Impact of demography on linked selection in two outcrossing Brassicaceae species

Genetic diversity is shaped by mutation, genetic drift, gene flow, recombination, and selection. The dynamics and interactions of these forces shape genetic diversity across different parts of the genome, between populations and species. Here, we have studied the effects of linked selection on nucleotide diversity in outcrossing populations of two Brassicaceae species, Arabidopsis lyrata and Capsella grandiflora, with contrasting demographic history. In agreement with previous estimates, we found evidence for a modest population size expansion thousands of generations ago, as well as efficient purifying selection in C. grandiflora. In contrast, the A. lyrata population exhibited evidence for very recent strong population size decline and weaker efficacy of purifying selection. Using multiple regression analyses with recombination rate and other genomic covariates as explanatory variables, we can explain 47% of the variance in neutral diversity in the C. grandiflora population, while in the A. lyrata population, only 11% of the variance was explained by the model. Recombination rate had a significant positive effect on neutral diversity in both species, suggesting that selection at linked sites has an effect on patterns of neutral variation. In line with this finding, we also found reduced neutral diversity in the vicinity of genes in the C. grandiflora population. However, in A. lyrata no such reduction in diversity was evident, a finding that is consistent with expectations of the impact of a recent bottleneck on patterns of neutral diversity near genes. This study thus empirically demonstrates how differences in demographic history modulate the impact of selection at linked sites in natural populations.     

High levels of population subdivision in a morphologically conserved Mediterranean toad (Alytes cisternasii) result from recent,multiple refugia: evidence from mtDNA,microsatellites and nuclear genealogies

Pleistocene glaciations often resulted in differentiation of taxa in southern European peninsulas, producing the high levels of endemism characteristic of these regions (e.g. the Iberian Peninsula). Despite their small ranges, endemic species often exhibit high levels of intraspecific differentiation as a result of a complex evolutionary history dominated by successive cycles of fragmentation, expansion and subsequent admixture of populations. Most evidence so far has come from the study of species with an Atlantic distribution in northwestern Iberia, and taxa restricted to Mediterranean‐type habitats remain poorly studied. The Iberian Midwife toad (Alytes cisternasii) is a morphologically conserved species endemic to southwestern and central Iberia and a typical inhabitant of Mediterranean habitats. Applying highly variable genetic markers from both mitochondrial and nuclear genomes to samples collected across the species’ range, we found evidence of high population subdivision within A. cisternasii. Mitochondrial haplotypes and microsatellites show geographically concordant patterns of genetic diversity, suggesting population fragmentation into several refugia during Pleistocene glaciations followed by subsequent events of geographical and demographic expansions with secondary contact. In addition, the absence of variation at the nuclear β‐fibint7 and Ppp3caint4 gene fragments suggests that populations of A. cisternasii have been recurrently affected by episodes of extinction and recolonization, and that documented patterns of population subdivision are the outcome of recent and multiple refugia. We discuss the evolutionary history of the species with particular interest in the increasing relevance of Mediterranean refugia for the survival of genetically differentiated populations during the Pleistocene glaciations as revealed by studies in co‐distributed taxa.     

Multilocus nuclear DNA markers reveal population structure and demography of Anopheles minimus

Utilization of multiple putatively neutral DNA markers for inferring evolutionary history of species population is considered to be the most robust approach. Molecular population genetic studies have been conducted in many species of Anopheles genus, but studies based on single nucleotide polymorphism (SNP) data are still very scarce. Anopheles minimus is one of the principal malaria vectors of Southeast (SE) Asia including the Northeastern (NE) India. Although population genetic studies with mitochondrial genetic variation data have been utilized to infer phylogeography of the SE Asian populations of this species, limited information on the population structure and demography of Indian An. minimus is available. We herewith have developed multilocus nuclear genetic approach with SNP markers located in X chromosome of An. minimus in eight Indian and two SE Asian population samples (121 individual mosquitoes in total) to infer population history and test several hypotheses on the phylogeography of this species. While the Thai population sample of An. minimus presented the highest nucleotide diversity, majority of the Indian samples were also fairly diverse. In general, An. minimus populations were moderately substructured in the distribution range covering SE Asia and NE India, largely falling under three distinct genetic clusters. Moreover, demographic expansion events could be detected in the majority of the presently studied populations of An. minimus. Additional DNA sequencing of the mitochondrial COII region in a subset of the samples (40 individual mosquitoes) corroborated the existing hypothesis of Indian An. minimus falling under the earlier reported mitochondrial lineage B.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Population genomic footprints of selection and associations with climate in natural populations of Arabidopsis halleri from the Alps

Natural genetic variation is essential for the adaptation of organisms to their local environment and to changing environmental conditions. Here, we examine genomewide patterns of nucleotide variation in natural populations of the outcrossing herb Arabidopsis halleri and associations with climatic variation among populations in the Alps. Using a pooled population sequencing (Pool‐Seq) approach, we discovered more than two million SNPs in five natural populations and identified highly differentiated genomic regions and SNPs using F_ST‐based analyses. We tested only the most strongly differentiated SNPs for associations with a nonredundant set of environmental factors using partial Mantel tests to identify topo‐climatic factors that may underlie the observed footprints of selection. Possible functions of genes showing signatures of selection were identified by Gene Ontology analysis. We found 175 genes to be highly associated with one or more of the five tested topo‐climatic factors. Of these, 23.4% had unknown functions. Genetic variation in four candidate genes was strongly associated with site water balance and solar radiation, and functional annotations were congruent with these environmental factors. Our results provide a genomewide perspective on the distribution of adaptive genetic variation in natural plant populations from a highly diverse and heterogeneous alpine environment.     

TESTING MODELS OF MIGRATION AND ISOLATION AMONG POPULATIONS OF CHINOOK SALMON (ONCORHYNCHUS TSCHAWYTSCHA)

The chinook salmon (Oncorhynchus tschawytscha) is a behaviorally, morphologically, and ecologically variable species distributed over a large geographic range. Although previous genetic surveys have revealed considerable genetic differences among populations with different life history types and from different major river drainages, it is not clear to what degree these genetically distinct populations are connected by low levels of gene flow. The work described in this paper addresses this question by surveying DNA restriction site variation at six nuclear genes from nine populations encompassing most of the species's North American range, and then attempting to fit the patterns of variation observed at these genes to five different evolutionary models using computer simulations of the coalescent process. Two commonly used constant population size models, one hypothesizing no migration among populations and one hypothesizing equal rates of migration among populations, were found to be statistically inconsistent with the observed patterns of variation. The other three models, which involved either recent divergence among populations coupled with large changes in populations size, unequal migration rates among populations, or selection, were all found to be consistent with the observed patterns of variation. Assuming selective neutrality, these results suggest that either the populations have all descended from a common ancestral population within the last ~50,000 years and have all suffered large declines in effective population size since that time, or that they have a more ancient divergence time but are connected by low levels of gene flow. These conclusions rest on the assumption of selective neutrality. With the methods employed, it was not possible to simultaneously test hypotheses of both selective neutrality and population structure.     

Evolutionary dynamics of quantitative variation in an adaptive trait at the regional scale: The case of zinc hyperaccumulation in Arabidopsis halleri

Metal hyperaccumulation in plants is an ecological trait whose biological significance remains debated, in particular because the selective pressures that govern its evolutionary dynamics are complex. One of the possible causes of quantitative variation in hyperaccumulation may be local adaptation to metalliferous soils. Here, we explored the population genetic structure of Arabidopsis halleri at fourteen metalliferous and nonmetalliferous sampling sites in southern Poland. The results were integrated with a quantitative assessment of variation in zinc hyperaccumulation to trace local adaptation. We identified a clear hierarchical structure with two distinct genetic groups at the upper level of clustering. Interestingly, these groups corresponded to different geographic subregions, rather than to ecological types (i.e., metallicolous vs. nonmetallicolous). Also, approximate Bayesian computation analyses suggested that the current distribution of A. halleri in southern Poland could be relictual as a result of habitat fragmentation caused by climatic shifts during the Holocene, rather than due to recent colonization of industrially polluted sites. In addition, we find evidence that some nonmetallicolous lowland populations may have actually derived from metallicolous populations. Meanwhile, the distribution of quantitative variation in zinc hyperaccumulation did separate metallicolous and nonmetallicolous accessions, indicating more recent adaptive evolution and diversifying selection between metalliferous and nonmetalliferous habitats. This suggests that zinc hyperaccumulation evolves both ways—towards higher levels at nonmetalliferous sites and lower levels at metalliferous sites. Our results open a new perspective on possible evolutionary relationships between A. halleri edaphic types that may inspire future genetic studies of quantitative variation in metal hyperaccumulation.     

Patterns of nucleotide sequence variation in <Emphasis Type="Italic">ICAM1</Emphasis> and <Emphasis Type="Italic">TNF</Emphasis> genes in twelve ethnic groups of India: roles of demographic history and natural selection

We have studied DNA sequence variation in and around the genes ICAM1 and TNF, which play functional and correlated roles in inflammatory processes and immune cell responses, in 12 diverse ethnic groups of India, with a view to investigating the relative roles of demographic history and natural selection in shaping the observed patterns of variation. The total numbers of single nucleotide polymorphisms (SNPs) detected at the ICAM1 and TNF loci were 29 and 12, respectively. Haplotype and allele frequencies differed significantly across populations. The site frequency spectra at these loci were significantly different from those expected under neutrality, and showed an excess of intermediate-frequency variants consistent with balancing selection. However, as expected under balancing selection, there was no significant reduction of F _ST values compared to neutral autosomal loci. Mismatch distributions were consistent with population expansion for both loci. On the other hand, the phylogenetic network among haplotypes for the TNF locus was similar to expectations under population expansion, while that for the ICAM1 was as expected under balancing selection. Nucleotide diversity at the ICAM1 locus was an order of magnitude lower in the promoter region, compared to the introns or exons, but no such difference was noted for the TNF gene. Thus, we conclude that the pattern of nucleotide variation in these genes has been modulated by both demographic history and selection. This is not surprising in view of the known allelic associations of several polymorphisms in these genes with various diseases, both infectious and noninfectious.     

Genetic diversity and differentiation among Lagopus lagopus populations in Scandinavia and Scotland: evolutionary significant units confirmed by SNP markers

Single Nucleotide Polymorphism in four Scandinavian populations of willow grouse (Lagopus lagopus) and two Scottish populations of red grouse (Lagopus lagopus scoticus) were assessed at 13 protein‐coding loci. We found high levels of diversity, with one substitution every 55 bp as an average and a total of 76 unlinked parsimony informative SNPs. Different estimators of genetic diversity such as: number of synonymous and non‐synonymous sites, average number of alleles, number and percentage of polymorphic loci, mean nucleotide diversity (π_s, π_a) and gene diversity at synonymous and non‐synonymous sites showed higher diversity in the northern populations compared to southern ones. Strong levels of purifying selection found in all the populations together with neutrality tests conforming to neutral expectations agree with large effective population sizes. Assignment tests reported a clear distinction between Scandinavian and Scottish grouse suggesting the existence of two different evolutionary significant units. The divergence time between willow and red grouse ranging between 12 500 and 125 000 years, in conjunction with the presence of ‘specific’ markers for each subspecies prompt a reassessment of the taxonomical status of the Scottish red grouse.     

Low genetic diversity and strong but shallow population differentiation suggests genetic homogenization by metapopulation dynamics in a social spider

Mating systems and population dynamics influence genetic diversity and structure. Species that experience inbreeding and limited gene flow are expected to evolve isolated, divergent genetic lineages. Metapopulation dynamics with frequent extinctions and colonizations may, on the other hand, deplete and homogenize genetic variation, if extinction rate is sufficiently high compared to the effect of drift in local demes. We investigated these theoretical predictions empirically in social spiders that are highly inbred. Social spiders show intranest mating, female‐biased sex ratio, and frequent extinction and colonization events, factors that deplete genetic diversity within nests and populations and limit gene flow. We characterized population genetic structure in Stegodyphus sarasinorum, a social spider distributed across the Indian subcontinent. Species‐wide genetic diversity was estimated over approximately 2800 km from Sri Lanka to Himalayas, by sequencing 16 protein‐coding nuclear loci. We found 13 SNPs in 6592 bp (π = 0.00045) indicating low species‐wide nucleotide diversity. Three genetic lineages were strongly differentiated; however, only one fixed difference among them suggests recent divergence. This is consistent with a scenario of metapopulation dynamics that homogenizes genetic diversity across the species' range. Ultimately, low standing genetic variation may hamper a species' ability to track environmental change and render social inbreeding spiders ‘evolutionary dead‐ends’.     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Unexpected high genetic diversity in small populations suggests maintenance by associative overdominance

The effective population size (N_e) is a central factor in determining maintenance of genetic variation. The neutral theory predicts that loss of variation depends on N_e, with less genetic drift in larger populations. We monitored genetic drift in 42 Drosophila melanogaster populations of different adult census population sizes (10, 50 or 500) using pooled RAD sequencing. In small populations, variation was lost at a substantially lower rate than expected. This observation was consistent across two ecological relevant thermal regimes, one stable and one with a stressful increase in temperature across generations. Estimated ratios between N_e and adult census size were consistently higher in small than in larger populations. The finding provides evidence for a slower than expected loss of genetic diversity and consequently a higher than expected long‐term evolutionary potential in small fragmented populations. More genetic diversity was retained in areas of low recombination, suggesting that associative overdominance, driven by disfavoured homozygosity of recessive deleterious alleles, is responsible for the maintenance of genetic diversity in smaller populations. Consistent with this hypothesis, the X‐chromosome, which is largely free of recessive deleterious alleles due to hemizygosity in males, fits neutral expectations even in small populations. Our experiments provide experimental answers to a range of unexpected patterns in natural populations, ranging from variable diversity on X‐chromosomes and autosomes to surprisingly high levels of nucleotide diversity in small populations.     

Genetic diversity and population structure in the mtDNA control region of Liza haematocheilus (Temminck & Schlegel, 1845)

The redlip mullet (Liza haematocheilus) is a commercial marine species widespread in the Northwestern Pacific. The genetic diversity and population structure of the redlip mullet collected from 2007 to 2010 along Chinese and Japanese coastal regions, were investigated by mitochondrial DNA. A total of 82 L. haematocheilus individuals were collected from five locations and 41 haplotypes were obtained. A pattern of distribution of genetic variability with high level of haplotype diversity (0.9642 ± 0.0104) and moderate nucleotide diversity (0.0162 ± 0.0085) was detected. The minimum spanning tree constructed with 41 haplotypes showed three divergent clades, corresponding to those defined in the NJ tree. Analyses of molecular variance and the population statistic F_ST also revealed significant genetic structures among populations of L. haematocheilus. The demographic history of L. haematocheilus was examined using neutrality tests and mismatch distribution analysis, which indicated a Pleistocene epoch population expansion. Comparisons of the Hakodate and Qingdao populations in different years demonstrated that marginal but significant (P < 0.05) differences exist in the two Qingdao populations, but with no significant difference (P > 0.05) in the Hakodate populations. Knowledge of their genetic diversity and genetic structure is crucial to establishing appropriate fishery management for the species.     

Phenotypic differentiation of the Solidago virgaurea complex along an elevational gradient: Insights from a common garden experiment and population genetics

Plant species distributed along wide elevational or latitudinal gradients show phenotypic variation due to their heterogeneous habitats. This study investigated whether phenotypic variation in populations of the Solidago virgaurea complex along an elevational gradient is caused by genetic differentiation. A common garden experiment was based on seeds collected from nine populations of the S. virgaurea complex growing at elevations from 1,597 m to 2,779 m a.s.l. on Mt. Norikura in central Japan. Population genetic analyses with microsatellite markers were used to infer the genetic structure and levels of gene flow between populations. Leaf mass per area was lower, while leaf nitrogen and chlorophyll concentrations were greater for higher elevations at which seeds were originally collected. For reproductive traits, plants derived from higher elevations had larger flower heads on shorter stems and flowering started earlier. These elevational changes in morphology were consistent with the clines in the field, indicating that phenotypic variation along the elevational gradient would have been caused by genetic differentiation. However, population genetic analysis using 16 microsatellite loci suggested an extremely low level of genetic differentiation of neutral genes among the nine populations. Analysis of molecular variance also indicated that most genetic variation was partitioned into individuals within a population, and the genetic differentiation among the populations was not significant. This study suggests that genome regions responsible for adaptive traits may differ among the populations despite the existence of gene flow and that phenotypic variation of the S. virgaurea complex along the elevational gradient is maintained by strong selection pressure.     

Population genetics of tandem trypsin inhibitor genes in Arabidopsis species with contrasting ecology and life history

Duplicated genes are important in the evolution and ecology of plant-defences because herbivore and pathogen attack can be countered via functional diversification at two levels: among duplicated loci and within loci. We explore molecular sequence variation for two members of a defence-related gene family, Arabidopsis thaliana trypsin inhibitors (ATTI), in A. thaliana and a closely related species, A. lyrata subspp. petraea. A worldwide sample of the inbreeding annual A. thaliana had less genetic variation at two ATTI loci (piTOTAL 相似文献   

Population structure and migration of the Tobacco Blue Mold Pathogen,Peronospora tabacina,into North America and Europe

Tobacco blue mold, caused by Peronospora tabacina, is an oomycete plant pathogen that causes yearly epidemics in tobacco (Nicotiana tabacum) in the United States and Europe. The genetic structure of P. tabacina was examined to understand genetic diversity, population structure and patterns of migration. Two nuclear loci, Igs2 and Ypt1, and one mitochondrial locus, cox2, were amplified, cloned and sequenced from fifty‐four isolates of P. tabacina from the United States, Central America–Caribbean–Mexico (CCAM), Europe and the Middle East (EULE). Cloned sequences from the three genes showed high genetic variability across all populations. Nucleotide diversity and the population mean mutation parameter per site (Watterson's theta) were higher in EULE and CCAM and lower in U.S. populations. Neutrality tests were significant and the equilibrium model of neutral evolution was rejected, indicating an excess of recent mutations or rare alleles. Hudson's S_nn tests were performed to examine population subdivision and gene flow among populations. An isolation‐with‐migration analysis (IM) supported the hypothesis of long‐distance migration of P. tabacina from the Caribbean region, Florida and Texas into other states in the United States. Within the European populations, the model documented migration from North Central Europe into western Europe and Lebanon, and migration from western Europe into Lebanon. The migration patterns observed support historical observations about the first disease introductions and movement in Europe. The models developed are applicable to other aerial dispersed emerging pathogens and document that high‐evolutionary‐risk plant pathogens can move over long distances to cause disease due to their large effective population size, population expansion and dispersal.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus,Osmeridae), an anadromous Pacific smelt

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction‐site‐associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single‐nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F_ST outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N_e) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome‐wide neutral and adaptive genetic variation for the applications of conservation and management.