Measuring the Sensitivity of Single-locus “Neutrality Tests” Using a Direct Perturbation Approach

A large number of statistical tests have been proposed to detect natural selection based on a sample of variation at a single genetic locus. These tests measure the deviation of the allelic frequency distribution observed within populations from the distribution expected under a set of assumptions that includes both neutral evolution and equilibrium population demography. The present study considers a new way to assess the statistical properties of these tests of selection, by their behavior in response to direct perturbations of the steady-state allelic frequency distribution, unconstrained by any particular nonequilibrium demographic scenario. Results from Monte Carlo computer simulations indicate that most tests of selection are more sensitive to perturbations of the allele frequency distribution that increase the variance in allele frequencies than to perturbations that decrease the variance. Simulations also demonstrate that it requires, on average, 4N generations (N is the diploid effective population size) for tests of selection to relax to their theoretical, steady-state distributions following different perturbations of the allele frequency distribution to its extremes. This relatively long relaxation time highlights the fact that these tests are not robust to violations of the other assumptions of the null model besides neutrality. Lastly, genetic variation arising under an example of a regularly cycling demographic scenario is simulated. Tests of selection performed on this last set of simulated data confirm the confounding nature of these tests for the inference of natural selection, under a demographic scenario that likely holds for many species. The utility of using empirical, genomic distributions of test statistics, instead of the theoretical steady-state distribution, is discussed as an alternative for improving the statistical inference of natural selection

Worldwide genetic relationships of pigs as inferred from X chromosome SNPs

The phylogeography of the porcine X chromosome has not been studied despite the unique characteristics of this chromosome. Here, we genotyped 59 single nucleotide polymorphisms (SNPs) in 312 pigs from around the world, representing 39 domestic breeds and wild boars in 30 countries. Overall, widespread commercial breeds showed the highest heterozygosity values, followed by African and American populations. Structuring, as inferred from FST and analysis of molecular variance, was consistently larger in the non-pseudoautosomal (NPAR) than in the pseudoautosomal regions (PAR). Our results show that genetic relationships between populations can vary widely between the NPAR and the PAR, underscoring the fact that their genetic trajectories can be quite different. NPAR showed an increased commercial-like genetic component relative to the PAR, probably because human selection processes to obtain individuals with high productive parameters were mediated by introgressing boars rather than sows.WBP is funded by COLCIENCIAS (Departamento Administrativo de Ciencia, Tecnología e Innovación, Francisco José de Caldas fellowship 497/2009, Colombia). CAS was funded by a PhD grant from CAPES (Coordenaçao de Aperfeiçoamento de Pessoal de Nível Superior, Brazil) and Universidade Catolica de Brasilia, Brazil. This work is funded by grants AGL2010-14822 (MICINN, Spain) to MPE, CGL2009-09346 (MICINN, Spain) to SERO, and Consolider project (MICINN, Spain) to Centre of Research in Agricultural Genomics.Peer reviewe

Meta-analysis of northeast Atlantic marine taxa shows contrasting phylogeographic patterns following post-LGM expansions

Background. Comparative phylogeography enables the study of historical and evolutionary processes that have contributed to shaping patterns of contemporary genetic diversity across co-distributed species. In this study, we explored genetic structure and historical demography in a range of coastal marine species across the northeast Atlantic to assess whether there are commonalities in phylogeographic patterns across taxa and to evaluate whether the timings of population expansions were linked to the Last Glacial Maximum (LGM). Methods. A literature search was conducted using Web of Science. Search terms were chosen to maximise the inclusion of articles reporting on population structure and phylogeography from the northeast Atlantic; titles and abstracts were screened to identify suitable articles within the scope of this study. Given the proven utility of mtDNA in comparative phylogeography and the availability of these data in the public domain, a meta-analysis was conducted using published mtDNA gene sequences. A standardised methodology was implemented to ensure that the genealogy and demographic history of all mtDNA datasets were reanalysed in a consistent and directly comparable manner. Results. Mitochondrial DNA datasets were built for 21 species. The meta-analysis revealed significant population differentiation in 16 species and four main types of haplotype network were found, with haplotypes in some species unique to specific geographical locations. A signal of rapid expansion was detected in 16 species, whereas five species showed evidence of a stable population size. Corrected mutation rates indicated that the majority of expansions were estimated to have occurred after the earliest estimate for the LGM (similar to 26.5 Kyr), while few expansions were estimated to have pre-dated the LGM. Conclusion. This study suggests that post-LGM expansion appeared to be common in a range of marine taxa, supporting the concept of rapid expansions after the LGM as the ice sheets started to retreat. However, despite the commonality of expansion patterns in many of these taxa, phylogeographic patterns appear to differ in the species included in this study. This suggests that species-specific evolutionary processes, as well as historical events, have likely influenced the distribution of genetic diversity of marine taxa in the northeast Atlantic

Natural history and molecular evolution of demersal Mediterranean sharks and skates inferred by comparative phylogeographic and demographic analyses

Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.5560#supplemental-informationBackground. The unique and complex paleoclimatic and paleogeographic events which affected the Mediterranean Sea since late Miocene deeply influenced the distribution and evolution of marine organisms and shaped their genetic structure. Following the Messinian salinity crisis and the sea-level fluctuations during the Pleistocene, several Mediterranean marine species developed deep genetic differentiation, and some underwent rapid radiation. Here, we consider two of the most prioritized groups for conservation in the light of their evolutionary history: sharks and rays (elasmobranchs). This paper deals with a comparative multispecies analysis of phylogeographic structure and historical demography in two pairs of sympatric, phylogenetically- and ecologically-related elasmobranchs, two scyliorhinid catsharks (Galeus melastomus, Scyliorhinus canicula) and two rajid skates (Raja clavata, Raja miraletus). Sampling and experimental analyses were designed to primarily test if the Sicilian Channel can be considered as effective eco-physiological barrier for Mediterranean demersal sympatric elasmobranchs. Methods. The phylogeography and the historical demography of target species were inferred by analysing the nucleotide variation of three mitochondrial DNA markers (i.e., partial sequence of COI, NADH2 and CR) obtained from a total of 248 individuals sampled in the Western and Eastern Mediterranean Sea as well as in the adjacent northeastern Atlantic Ocean. Phylogeographic analysis was performed by haplotype networking and testing spatial genetic differentiation of samples (i.e., analysis of molecular variance and of principal components). Demographic history of Mediterranean populations was reconstructed using mismatch distribution and Bayesian Skyline Plot analyses. Results. No spatial genetic differentiation was identified in either catshark species, while phylogeographic structure of lineages was identified in both skates, with R. miraletus more structured than R. clavata. However, such structuring of skate lineages waSupplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.5560#supplemental-informatio

Diversification across an altitudinal gradient in the Tiny Greenbul (Phyllastrephus debilis) from the Eastern Arc Mountains of Africa

<p>Abstract</p> <p>Background</p> <p>The Eastern Arc Mountains of Africa have become one of the focal systems with which to explore the patterns and mechanisms of diversification among montane species and populations. One unresolved question is the extent to which populations inhabiting montane forest interact with those of adjacent lowland forest abutting the coast of eastern Africa. The Tiny Greenbul (<it>Phyllastephus debilis</it>) represents the only described bird species within the Eastern Arc/coastal forest mosaic, which is polytypic across an altitudinal gradient: the subspecies <it>albigula </it>(green head) is distributed in the montane Usambara and Nguru Mountains whereas the subspecies <it>rabai </it>(grey head) is found in Tanzanian lowland and foothill forest. Using a combination of morphological and genetic data, we aim to establish if the pattern of morphological differentiation in the Tiny Greenbul (<it>Phyllastrephus debilis</it>) is the result of disruptive selection along an altitudinal gradient or a consequence of secondary contact following population expansion of two differentiated lineages.</p> <p>Results</p> <p>We found significant biometric differences between the lowland (<it>rabai</it>) and montane (<it>albigula</it>) populations in Tanzania. The differences in shape are coupled with discrete differences in the coloration of the underparts. Using multi-locus data gathered from 124 individuals, we show that lowland and montane birds form two distinct genetic lineages. The divergence between the two forms occurred between 2.4 and 3.1 Myrs ago.</p> <p>Our coalescent analyses suggest that limited gene flow, mostly from the subspecies <it>rabai </it>to <it>albigula</it>, is taking place at three mid-altitude localities, where lowland and montane rainforest directly abut. The extent of this introgression appears to be limited and is likely a consequence of the recent expansion of <it>rabai </it>further inland.</p> <p>Conclusion</p> <p>The clear altitudinal segregation in morphology found within the Tiny Greenbul is the result of secondary contact of two highly differentiated lineages rather than disruptive selection in plumage pattern across an altitudinal gradient. Based on our results, we recommend <it>albigula </it>be elevated to species rank.</p

Fully Bayesian tests of neutrality using genealogical summary statistics

<p>Abstract</p> <p>Background</p> <p>Many data summary statistics have been developed to detect departures from neutral expectations of evolutionary models. However questions about the neutrality of the evolution of genetic loci within natural populations remain difficult to assess. One critical cause of this difficulty is that most methods for testing neutrality make simplifying assumptions simultaneously about the mutational model and the population size model. Consequentially, rejecting the null hypothesis of neutrality under these methods could result from violations of either or both assumptions, making interpretation troublesome.</p> <p>Results</p> <p>Here we harness posterior predictive simulation to exploit summary statistics of both the data and model parameters to test the goodness-of-fit of standard models of evolution. We apply the method to test the selective neutrality of molecular evolution in non-recombining gene genealogies and we demonstrate the utility of our method on four real data sets, identifying significant departures of neutrality in human influenza A virus, even after controlling for variation in population size.</p> <p>Conclusion</p> <p>Importantly, by employing a full model-based Bayesian analysis, our method separates the effects of demography from the effects of selection. The method also allows multiple summary statistics to be used in concert, thus potentially increasing sensitivity. Furthermore, our method remains useful in situations where analytical expectations and variances of summary statistics are not available. This aspect has great potential for the analysis of temporally spaced data, an expanding area previously ignored for limited availability of theory and methods.</p

A Role in Immunity for Arabidopsis Cysteine Protease RD21, the Ortholog of the Tomato Immune Protease C14

Secreted papain-like Cys proteases are important players in plant immunity. We previously reported that the C14 protease of tomato is targeted by cystatin-like EPIC proteins that are secreted by the oomycete pathogen Phytophthora infestans (Pinf) during infection. C14 has been under diversifying selection in wild potato species coevolving with Pinf and reduced C14 levels result in enhanced susceptibility for Pinf. Here, we investigated the role C14-EPIC-like interactions in the natural pathosystem of Arabidopsis with the oomycete pathogen Hyaloperonospora arabidopsidis (Hpa). In contrast to the Pinf-solanaceae pathosystem, the C14 orthologous protease of Arabidopsis, RD21, does not evolve under diversifying selection in Arabidopsis, and rd21 null mutants do not show phenotypes upon compatible and incompatible Hpa interactions, despite the evident lack of a major leaf protease. Hpa isolates express highly conserved EPIC-like proteins during infections, but it is unknown if these HpaEPICs can inhibit RD21 and one of these HpaEPICs even lacks the canonical cystatin motifs. The rd21 mutants are unaffected in compatible and incompatible interactions with Pseudomonas syringae pv. tomato, but are significantly more susceptible for the necrotrophic fungal pathogen Botrytis cinerea, demonstrating that RD21 provides immunity to a necrotrophic pathogen

Islands beneath islands: phylogeography of a groundwater amphipod crustacean in the Balearic archipelago

<p>Abstract</p> <p>Background</p> <p>Metacrangonyctidae (Amphipoda, Crustacea) is an enigmatic continental subterranean water family of marine origin (thalassoid). One of the species in the genus, <it>Metacrangonyx longipes</it>, is endemic to the Balearic islands of Mallorca and Menorca (W Mediterranean). It has been suggested that the origin and distribution of thalassoid crustaceans could be explained by one of two alternative hypotheses: (1) active colonization of inland freshwater aquifers by a marine ancestor, followed by an adaptative shift; or (2) passive colonization by stranding of ancestral marine populations in coastal aquifers during marine regressions. A comparison of phylogenies, phylogeographic patterns and age estimations of clades should discriminate in favour of one of these two proposals.</p> <p>Results</p> <p>Phylogenetic relationships within <it>M. longipes </it>based on three mitochondrial DNA (mtDNA) and one nuclear marker revealed five genetically divergent and geographically structured clades. Analyses of cytochrome oxidase subunit 1 (<it>cox1</it>) mtDNA data showed the occurrence of a high geographic population subdivision in both islands, with current gene flow occurring exclusively between sites located in close proximity. Molecular-clock estimations dated the origin of <it>M. longipes </it>previous to about 6 Ma, whereas major cladogenetic events within the species took place between 4.2 and 2.0 Ma.</p> <p>Conclusions</p> <p><it>M. longipes </it>displayed a surprisingly old and highly fragmented population structure, with major episodes of cladogenesis within the species roughly correlating with some of the major marine transgression-regression episodes that affected the region during the last 6 Ma. Eustatic changes (vicariant events) -not active range expansion of marine littoral ancestors colonizing desalinated habitats-explain the phylogeographic pattern observed in <it>M. longipes</it>.</p

Present and past climatic effects on the current distribution and genetic diversity of the Iberian spadefoot toad (Pelobates cultripes): an integrative approach

Aim: Predicting species responses to global change is one of the most pressing issues in Conservation Biogeography. A key part of the problem is understanding how organisms have reacted to climatic changes in the past. Here we use species distribution modelling to infer the effects of climate changes since the Last Interglacial (LIG, about 130,000 ybp) on patterns of genetic structure and diversity in the Western Spadefoot toad (Pelobates cultripes) in combination with spatially-explicit phylogeographic analyses. Location: Iberian Peninsula and mainland France. Methods: 524 individuals from 54 populations across the species range were sampled to document patterns of genetic diversity and infer their evolutionary history based on data from mtDNA and fourteen polymorphic microsatellites. Generalized linear models based on distribution data were used to infer climatic favourability for the species in the present and in projected scenarios for the LIG, the Mid Holocene and the last glacial maximum (LGM). Results: Estimates of genetic diversity show a decreasing trend from south to north, suggesting persistence of higher historical population sizes in the southern Iberian Peninsula. Species distribution models show differences in climatic favourability through time, with significant correlations between historically stable favourable areas and current patterns of genetic diversity. These results are corroborated by Bayesian Skyline Plots and continuous diffusion phylogeographic analyses. Main conclusions: The results indicate the presence of southern refugia, with moderate recent expansions at the northern end of the species’ range. Populations at the northern range margin exhibit the lowest genetic diversity and occupy historically unstable areas, classified as marginal in terms of favourability, rendering them most vulnerable to climate-mediated changes in the medium to long term

Low diversity Cryptococcus neoformans variety grubii multilocus sequence types from Thailand are consistent with an ancestral African origin.

Published versio