Demographic processes shaping genetic variation

Genetic analyses reveal panmixia in Indian waters and population subdivision across Indian Ocean and Indo-Malay Archipelago for Decapterus russelli

The Indian Scad, Decapterus russelli is an important pelagic carangid widely distributed throughout the Indian Ocean and the Indo-West Pacific. Despite being widely distributed in the Indian Ocean, the information regarding genetic structuring and diversity has been lacking compared to its Indo Malay Archipelago counterparts. The present study was conducted to investigate the genetic stock structure of D. russelli based on mitochondrial (Cyt b) and nuclear (DrAldoB1) markers along Indian waters. The results indicated the presence of a single panmictic stock across the Indian Ocean region. High haplotype diversity associated with low nucleotide diversity suggested a population bottleneck followed by rapid population growth. Phylogenetic analysis revealed the absence of geographical clustering of lineages with the most common haplotype distributed globally. The pelagic life style, migratory capabilities, and larval dispersal may be the contributing factors to the observed spatial homogeneity of D. russelli. However, significant genetic differentiation was observed between the populations from Indian Ocean and Indo-Malay Archipelago. Hierarchical molecular variance analysis (AMOVA), pairwise FST comparisons and SAMOVA showed existence of two distinct genetic stocks of D. russelli in the Indian Ocean and IMA. The observed interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the Indo-Malay Archipelago and the Indian Ocean regions suggest the influence of historic isolation, ocean surface currents and biotic and abiotic features of the ocean. Also, there was a significant relationship between genetic distance and geographical distance between population pairs in a manner consistent with isolation-by-distance. These resulted in the evolution of a phylogeographic break for this species between these regions. The findings of these results suggest that D. russelli from the Indian Ocean shall be managed in its entire area of distribution as a single stock. Further, the Indian Ocean and Indo-Malayan stocks can be managed separately.

Reduced within-population quantitative genetic variation is associated with climate harshness in maritime pine

How evolutionary forces interact to maintain genetic variation within populations has been a matter of extensive theoretical debates. While mutation and exogenous gene flow increase genetic variation, stabilizing selection and genetic drift are expected to deplete it. To date, levels of genetic variation observed in natural populations are hard to predict without accounting for other processes, such as balancing selection in heterogeneous environments. We aimed to empirically test three hypotheses: (i) admixed populations have higher quantitative genetic variation due to introgression from other gene pools, (ii) quantitative genetic variation is lower in populations from harsher environments (i.e., experiencing stronger selection), and (iii) quantitative genetic variation is higher in populations from heterogeneous environments. Using growth, phenological and functional trait data from three clonal common gardens and 33 populations (522 clones) of maritime pine (Pinus pinaster Aiton), we estimated the association between the population-specific total genetic variances (i.e., among-clone variances) for these traits and ten population-specific indices related to admixture levels (estimated based on 5165 SNPs), environmental temporal and spatial heterogeneity and climate harshness. Populations experiencing colder winters showed consistently lower genetic variation for early height growth (a fitness-related trait in forest trees) in the three common gardens. Within-population quantitative genetic variation was not associated with environmental heterogeneity or population admixture for any trait. Our results provide empirical support for the potential role of natural selection in reducing genetic variation for early height growth within populations, which indirectly gives insight into the adaptive potential of populations to changing environments.

Genetic signatures of small effective population sizes and demographic declines in an endangered rattlesnake, Sistrurus catenatus

Endangered species that exist in small isolated populations are at elevated risk of losing adaptive variation due to genetic drift. Analyses that estimate short-term effective population sizes, characterize historical demographic processes, and project the trajectory of genetic variation into the future are useful for predicting how levels of genetic diversity may change. Here, we use data from two independent types of genetic markers (single nucleotide polymorphisms [SNPs] and microsatellites) to evaluate genetic diversity in 17 populations spanning the geographic range of the endangered eastern massasauga rattlesnake (Sistrurus catenatus). First, we use SNP data to confirm previous reports that these populations exhibit high levels of genetic structure (overall Fst = 0.25). Second, we show that most populations have contemporary Ne estimates <50. Heterozygosity-fitness correlations in these populations provided no evidence for a genetic cost to living in small populations, though these tests may lack power. Third, model-based demographic analyses of individual populations indicate that all have experienced declines, with the onset of many of these declines occurring over timescales consistent with anthropogenic impacts (<200 years). Finally, forward simulations of the expected loss of variation in relatively large (Ne = 50) and small (Ne = 10) populations indicate they will lose a substantial amount of their current standing neutral variation (63% and 99%, respectively) over the next 100 years. Our results argue that drift has a significant and increasing impact on levels of genetic variation in isolated populations of this snake, and efforts to assess and mitigate associated impacts on adaptive variation should be components of the management of this endangered reptile.

Fluctuating selection across years and phenotypic variation in food-deceptive orchids

Nectarless flowers that deceive pollinators offer an opportunity to study asymmetric plant-insect interactions. Orchids are a widely used model for studying these interactions because they encompass several thousand species adopting deceptive pollination systems. High levels of intra-specific phenotypic variation have been reported in deceptive orchids, suggesting a reduced consistency of pollinator-mediated selection on their floral traits. Nevertheless, several studies report on widespread directional selection mediated by pollinators even in these deceptive orchids. In this study we test the hypothesis that the observed selection can fluctuate across years in strength and direction thus likely contributing to the phenotypic variability of this orchid group. We performed a three-year study estimating selection differentials and selection gradients for nine phenotypic traits involved in insect attraction in two Mediterranean orchid species, namely Orchis mascula and O. pauciflora, both relying on a well-described food-deceptive pollination strategy. We found weak directional selection and marginally significant selection gradients in the two investigated species with significant intra-specific differences in selection differentials across years. Our data do not link this variation with a specific environmental cause, but our results suggest that pollinator-mediated selection in food-deceptive orchids can change in strength and in direction over time. In perennial plants, such as orchids, different selection differentials in the same populations in different flowering seasons can contribute to the maintenance of phenotypic variation often reported in deceptive orchids.

Genetic diversity and structure of Lolium perenne ssp. multiflorum in California vineyards and orchards indicate potential for spread of herbicide resistance via gene flow

Management of agroecosystems with herbicides imposes strong selection pressures on weedy plants leading to the evolution of resistance against those herbicides. Resistance to glyphosate in populations of Lolium perenne L. ssp. multiflorum is increasingly common in California, USA, causing economic losses and the loss of effective management tools. To gain insights into the recent evolution of glyphosate resistance in L. perenne in perennial cropping systems of northwest California and to inform management, we investigated the frequency of glyphosate resistance and the genetic diversity and structure of 14 populations. The sampled populations contained frequencies of resistant plants ranging from 10% to 89%. Analyses of neutral genetic variation using microsatellite markers indicated very high genetic diversity within all populations regardless of resistance frequency. Genetic variation was distributed predominantly among individuals within populations rather than among populations or sampled counties, as would be expected for a wide-ranging outcrossing weed species. Bayesian clustering analysis provided evidence of population structuring with extensive admixture between two genetic clusters or gene pools. High genetic diversity and admixture, and low differentiation between populations, strongly suggest the potential for spread of resistance through gene flow and the need for management that limits seed and pollen dispersal in L. perenne.

Escape to Ferality: The Endoferal Origin of Weedy Rice from Crop Rice through De-Domestication

Domestication is the hallmark of evolution and civilization and harnesses biodiversity through selection for specific traits. In regions where domesticated lines are grown near wild relatives, congeneric sources of aggressive weedy genotypes cause major economic losses. Thus, the origins of weedy genotypes where no congeneric species occur raise questions regarding management effectiveness and evolutionary mechanisms responsible for weedy population success. Since eradication in the 1970s, California growers avoided weedy rice through continuous flood culture and zero-tolerance guidelines, preventing the import, presence, and movement of weedy seeds. In 2003, after decades of no reported presence in California, a weedy rice population was confirmed in dry-seeded fields. Our objectives were to identify the origins and establishment of this population and pinpoint possible phenotypes involved. We show that California weedy rice is derived from a different genetic source among a broad range of AA genome Oryzas and is most recently diverged from O. sativa temperate japonica cultivated in California. In contrast, other weedy rice ecotypes in North America (Southern US) originate from weedy genotypes from China near wild Oryza, and are derived through existing crop-wild relative crosses. Analyses of morphological data show that California weedy rice subgroups have phenotypes like medium-grain or gourmet cultivars, but have colored pericarp, seed shattering, and awns like wild relatives, suggesting that reversion to non-domestic or wild-like traits can occur following domestication, despite apparent fixation of domestication alleles. Additionally, these results indicate that preventive methods focused on incoming weed sources through contamination may miss burgeoning weedy genotypes that rapidly adapt, establish, and proliferate. Investigating the common and unique evolutionary mechanisms underlying global weed origins and subsequent interactions with crop relatives sheds light on how weeds evolve and addresses broader questions regarding the stability of selection during domestication and crop improvement.

Human-modified habitats change patterns of population genetic structure and group relatedness in Peter's tent-roosting bats

Although coloniality is widespread among mammals, it is still not clear what factors influence composition of social groups. As animals need to adapt to multiple habitat and environmental conditions throughout their range, variation in group composition should be influenced by adaptive adjustment to different ecological factors. Relevant to anthropogenic disturbance, increased habitat modification by humans can alter species' presence, density, and population structure. Therefore, it is important to understand the consequences of changes to landscape composition, in particular how habitat modification affects social structure of group-forming organisms. Here, we combine information on roosting associations with genetic structure of Peter's tent-roosting bats, Uroderma bilobatum to address how different habitat characteristics at different scales affect structure of social groups. By dividing analyses by age and sex, we determined that genetic structure was greater for adult females than adult males or offspring. Habitat variables explained 80% of the variation in group relatedness (mainly influenced by female relatedness) with roost characteristics contributing the most explained variation. This suggests that females using roosts of specific characteristics exhibit higher relatedness and seem to be philopatric. These females mate with more males than do more labile female groups. Results describe ecological and microevolutionary processes, which affect relatedness and social structure; findings are highly relevant to species distributions in both natural and human-modified environments.

Reconstructing the invasion history of Heracleum persicum (Apiaceae) into Europe

Sparse, incomplete and inappropriate historical records of invasive species often hamper invasive species management interventions. Population genetic analyses of invaders might provide a suitable context for the identification of their source populations and possible introduction routes. Here, we describe the population genetics of Heracleum persicum Desf. ex Fisch and trace its route of introduction into Europe. Microsatellite markers revealed a significantly higher genetic diversity of H. persicum in its native range, and the loss of diversity in the introduced range may be attributed to a recent genetic bottleneck. Bayesian cluster analysis on regional levels identified three and two genetic clusters in the native and the introduced ranges, respectively. A global structure analysis revealed two worldwide distinct genetic groups: one primarily in Iran and Denmark, the other primarily in Norway. There were also varying degrees of admixture in England, Sweden, Finland and Latvia. Approximate Bayesian computation indicated two independent introductions of H. persicum from Iran to Europe: the first one in Denmark and the second one in England. Finland was subsequently colonized by English populations. In contrast to the contemporary hypothesis of English origin of Norwegian populations, we found Finland to be a more likely source for Norwegian populations, a scenario supported by higher estimated historical migration from Finland to Norway. Genetic diversity per se is not a primary determinant of invasiveness in H. persicum. Our results indicate that, due to either pre-adaptations or rapid local adaptations, introduced populations may have acquired invasiveness after subsequent introductions, once a suitable environment was encountered.

Genetic models reveal historical patterns of sea lamprey population fluctuations within Lake Champlain

The origin of sea lamprey (Petromyzon marinus) in Lake Champlain has been heavily debated over the past decade. Given the lack of historical documentation, two competing hypotheses have emerged in the literature. First, it has been argued that the relatively recent population size increase and concomitant rise in wounding rates on prey populations are indicative of an invasive population that entered the lake through the Champlain Canal. Second, recent genetic evidence suggests a post-glacial colonization at the end of the Pleistocene, approximately 11,000 years ago. One limitation to resolving the origin of sea lamprey in Lake Champlain is a lack of historical and current measures of population size. In this study, the issue of population size was explicitly addressed using nuclear (nDNA) and mitochondrial DNA (mtDNA) markers to estimate historical demography with genetic models. Haplotype network analysis, mismatch analysis, and summary statistics based on mtDNA noncoding sequences for NCI (479 bp) and NCII (173 bp) all indicate a recent population expansion. Coalescent models based on mtDNA and nDNA identified two potential demographic events: a population decline followed by a very recent population expansion. The decline in effective population size may correlate with land-use and fishing pressure changes post-European settlement, while the recent expansion may be associated with the implementation of the salmonid stocking program in the 1970s. These results are most consistent with the hypothesis that sea lamprey are native to Lake Champlain; however, the credibility intervals around parameter estimates demonstrate that there is uncertainty regarding the magnitude and timing of past demographic events.

Models of selection, isolation, and gene flow in speciation

Many marine ecologists aspire to use genetic data to understand how selection and demographic history shape the evolution of diverging populations as they become reproductively isolated species. I propose combining two types of genetic analysis focused on this key early stage of the speciation process to identify the selective agents directly responsible for population divergence. Isolation-with-migration (IM) models can be used to characterize reproductive isolation between populations (low gene flow), while codon models can be used to characterize selection for population differences at the molecular level (especially positive selection for high rates of amino acid substitution). Accessible transcriptome sequencing methods can generate the large quantities of data needed for both types of analysis. I highlight recent examples (including our work on fertilization genes in sea stars) in which this confluence of interest, models, and data has led to taxonomically broad advances in understanding marine speciation at the molecular level. I also highlight new models that incorporate both demography and selection: simulations based on these theoretical advances suggest that polymorphisms shared among individuals (a key source of information in IM models) may lead to false-positive evidence of selection (in codon models), especially during the early stages of population divergence and speciation that are most in need of study. The false-positive problem may be resolved through a combination of model improvements plus experiments that document the phenotypic and fitness effects of specific polymorphisms for which codon models and IM models indicate selection and reproductive isolation (such as genes that mediate sperm-egg compatibility at fertilization).

Maximum-likelihood inference of population size contractions from microsatellite data

Understanding the demographic history of populations and species is a central issue in evolutionary biology and molecular ecology. In this work, we develop a maximum-likelihood method for the inference of past changes in population size from microsatellite allelic data. Our method is based on importance sampling of gene genealogies, extended for new mutation models, notably the generalized stepwise mutation model (GSM). Using simulations, we test its performance to detect and characterize past reductions in population size. First, we test the estimation precision and confidence intervals coverage properties under ideal conditions, then we compare the accuracy of the estimation with another available method (MSVAR) and we finally test its robustness to misspecification of the mutational model and population structure. We show that our method is very competitive compared with alternative ones. Moreover, our implementation of a GSM allows more accurate analysis of microsatellite data, as we show that the violations of a single step mutation assumption induce very high bias toward false contraction detection rates. However, our simulation tests also showed some limits, which most importantly are large computation times for strong disequilibrium scenarios and a strong influence of some form of unaccounted population structure. This inference method is available in the latest implementation of the MIGRAINE software package.

Bioclimatic regions influence genetic structure of four Jordanian Stipa species

Strong environmental gradients can affect the genetic structure of plant populations, but little is known as to whether closely related species respond similarly or idiosyncratically to ecogeographic variation. We analysed the extent to which gradients in temperature and rainfall shape the genetic structure of four Stipa species in four bioclimatic regions in Jordan. Genetic diversity, differentiation and structure of Stipa species were investigated using amplified fragment length polymorphism (AFLP) molecular markers. For each of the four study species, we sampled 120 individuals from ten populations situated in distinct bioclimatic regions and assessed the degree of genetic diversity and genetic differentiation within and among populations. The widespread ruderals Stipa capensis and S. parviflora had higher genetic diversity than the geographically restricted semi-desert species S. arabica and S. lagascae. In three of the four species, genetic diversity strongly decreased with precipitation, while genetic diversity increased with temperature in S. capensis. Most genetic diversity resided among populations in the semi-desert species (Φ(ST) = 0.572/0.595 in S. arabica/lagascae) but within populations in the ruderal species (Φ(ST) = 0.355/0.387 S. capensis/parviflora). Principal coordinate analysis (PCoA) and STRUCTURE analysis showed that Stipa populations of all species clustered ecogeographically. A genome scan revealed that divergent selection at particular AFLP loci contributed to genetic differentiation. Irrespective of their different life histories, Stipa species responded similarly to the bioclimatic gradient in Jordan. We conclude that, in addition to predominant random processes, steep climatic gradients might shape the genetic structure of plant populations.

Signatures of demography and recombination at coding genes in naturally-distributed populations of Arabidopsis lyrata subsp. petraea

Demography impacts the observed standing level of genetic diversity present in populations. Distinguishing the relative impacts of demography from selection requires a baseline of expressed gene variation in naturally occurring populations. Six nuclear genes were sequenced to estimate the patterns and levels of genetic diversity in natural Arabidopsis lyrata subsp. petraea populations that differ in demographic histories since the Pleistocene. As expected, northern European populations have genetic signatures of a strong population bottleneck likely due to glaciation during the Pleistocene. Levels of diversity in the northern populations are about half of that in central European populations. Bayesian estimates of historical population size changes indicate that central European populations also have signatures of population size change since the last glacial maxima, suggesting that these populations are not as stable as previously thought. Time since divergence amongst northern European populations is higher than amongst central European populations, suggesting that the northern European populations were established before the Pleistocene and survived glaciation in small separated refugia. Estimates of demography based on expressed genes are complementary to estimates based on microsatellites and transposable elements, elucidating temporal shifts in population dynamics and confirming the importance of marker selection for tests of demography.

History of expansion and anthropogenic collapse in a top marine predator of the Black Sea estimated from genetic data

Two major ecological transitions marked the history of the Black Sea after the last Ice Age. The first was the postglacial transition from a brackish-water to a marine ecosystem dominated by porpoises and dolphins once this basin was reconnected back to the Mediterranean Sea (ca. 8,000 y B.P.). The second occurred during the past decades, when overfishing and hunting activities brought these predators close to extinction, having a deep impact on the structure and dynamics of the ecosystem. Estimating the extent of this decimation is essential for characterizing this ecosystem's dynamics and for formulating restoration plans. However, this extent is poorly documented in historical records. We addressed this issue for one of the main Black Sea predators, the harbor porpoise, using a population genetics approach. Analyzing its genetic diversity using an approximate Bayesian computation approach, we show that only a demographic expansion (at most 5,000 y ago) followed by a contemporaneous population collapse can explain the observed genetic data. We demonstrate that both the postglacial settlement of harbor porpoises in the Black Sea and the recent anthropogenic activities have left a clear footprint on their genetic diversity. Specifically, we infer a strong population reduction (~90%) that occurred within the past 5 decades, which can therefore clearly be related to the recent massive killing of small cetaceans and to the continuing incidental catches in commercial fisheries. Our study thus provides a quantitative assessment of these demographically catastrophic events, also showing that two separate historical events can be inferred from contemporary genetic data.

The emergence and maintenance of sickle cell hotspots in the Mediterranean

Genetic disorders of haemoglobin (haemoglobinopathies), including the thalassaemias and sickle cell anaemia, abound in historically malarious regions, due to the protection they provide against death from severe malaria. Despite the overall spatial correlation between malaria and these disorders, inter-population differences exist in the precise combinations of haemoglobinopathies observed. Greece and Italy present a particularly interesting case study: their high frequencies of beta thalassaemia speak to a history of intense malaria selection, yet they possess very little of the strongly malaria protective mutation responsible for sickle cell anaemia, despite historical migrational links with Africa where high frequencies of sickle cell occur. Twentieth century surveys of beta thalassaemia and sickle cell in Greece, Sicily and Sardinia have revealed striking sickle cell ‘hotspots’ – places where the frequency of sickle cell approaches that seen in Africa while neighbouring populations remain relatively sickle cell free. It remains unclear how these hotspots have been maintained over time without sickle cell spreading throughout the region. Here we use a metapopulation model to show that (i) epistasis between the alpha and beta forms of thalassaemia can restrict the spread of sickle cell through a network of linked subpopulations and (ii) the emergence of sickle cell hotspots requires relatively low levels of gene flow, but the aforementioned epistasis increases the chances of hotspots forming.

Inferring population decline and expansion from microsatellite data: a simulation-based evaluation of the Msvar method

Reconstructing the demographic history of populations is a central issue in evolutionary biology. Using likelihood-based methods coupled with Monte Carlo simulations, it is now possible to reconstruct past changes in population size from genetic data. Using simulated data sets under various demographic scenarios, we evaluate the statistical performance of Msvar, a full-likelihood Bayesian method that infers past demographic change from microsatellite data. Our simulation tests show that Msvar is very efficient at detecting population declines and expansions, provided the event is neither too weak nor too recent. We further show that Msvar outperforms two moment-based methods (the M-ratio test and Bottleneck) for detecting population size changes, whatever the time and the severity of the event. The same trend emerges from a compilation of empirical studies. The latest version of Msvar provides estimates of the current and the ancestral population size and the time since the population started changing in size. We show that, in the absence of prior knowledge, Msvar provides little information on the mutation rate, which results in biased estimates and/or wide credibility intervals for each of the demographic parameters. However, scaling the population size parameters with the mutation rate and scaling the time with current population size, as coalescent theory requires, significantly improves the quality of the estimates for contraction but not for expansion scenarios. Finally, our results suggest that Msvar is robust to moderate departures from a strict stepwise mutation model.

Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests

Deviations of the site frequency spectrum of mutations (SFS) from neutral expectations may be caused by natural selection or by demographic processes such as population subdivision or temporal changes in population size. As most widespread temperate and boreal tree species have expanded from glacial refugia in the past 13,000 years, colonization bottlenecks associated with this migration may have left variable demographic signatures among geographic populations corresponding to distance from the refugia. To determine whether the signature of postglacial re-colonization has skewed the SFS in the widely distributed conifer Sitka spruce (Picea sitchensis (Bong.) Carr.), we re-sequenced 153 nuclear genes in six populations from across the species range. We found that while the SFS for the pooled sample produced negative values for Tajima's D and Fay and Wu's H, these statistics exhibited strong clinal variation when populations were analysed separately (R(2) = 0.84, P = 0.007 for Tajima's D and R(2)  = 0.65, P = 0.033 for Fay and Wu's H). When historical bottlenecks of varying age were simulated using approximate Bayesian computation, distance of populations from the southern range limit explained most of the variation in bottleneck timing among populations (R(2) = 0.89, P = 0.003). These data suggest that sequential population bottlenecks during postglacial re-colonization have resulted in diverse among-population signatures within the contemporary SFS in Sitka spruce, with rare variants more common in the south, and medium-frequency variants more common in the north. Our results also emphasize the need to consider sampling strategy and to explore population-specific null demographic models in surveys of nucleotide variation in widely distributed species.

Demographic and genetic patterns of variation among populations of Arabidopsis thaliana from contrasting native environments

BACKGROUND

Understanding the relationship between environment and genetics requires the integration of knowledge on the demographic behavior of natural populations. However, the demographic performance and genetic composition of Arabidopsis thaliana populations in the species' native environments remain largely uncharacterized. This information, in combination with the advances on the study of gene function, will improve our understanding on the genetic mechanisms underlying adaptive evolution in A. thaliana.

METHODOLOGY/PRINCIPAL FINDINGS

We report the extent of environmental, demographic, and genetic variation among 10 A. thaliana populations from Mediterranean (coastal) and Pyrenean (montane) native environments in northeast Spain. Geographic, climatic, landscape, and soil data were compared. Demographic traits, including the dynamics of the soil seed bank and the attributes of aboveground individuals followed over a complete season, were also analyzed. Genetic data based on genome-wide SNP markers were used to describe genetic diversity, differentiation, and structure. Coastal and montane populations significantly differed in terms of environmental, demographic, and genetic characteristics. Montane populations, at higher altitude and farther from the sea, are exposed to colder winters and prolonged spring moisture compared to coastal populations. Montane populations showed stronger secondary seed dormancy, higher seedling/juvenile mortality in winter, and initiated flowering later than coastal populations. Montane and coastal regions were genetically differentiated, montane populations bearing lower genetic diversity than coastal ones. No significant isolation-by-distance pattern and no shared multilocus genotypes among populations were detected.

CONCLUSIONS/SIGNIFICANCE

Between-region variation in climatic patterns can account for differences in demographic traits, such as secondary seed dormancy, plant mortality, and recruitment, between coastal and montane A. thaliana populations. In addition, differences in plant mortality can partly account for differences in the genetic composition of coastal and montane populations. This study shows how the interplay between variation in environmental, demographic, and genetic parameters may operate in natural A. thaliana populations.