Genome-wide population genetic analysis identifies evolutionary forces establishing continuous population divergence

Elucidating the mechanism shaping the spatial variations of traits has long been a central concern of evolutionary biologists. Geographic clines of allele/morph frequencies along environmental gradients are suggested to be established and maintained by the balancing of two opposing evolutionary forces, namely selection that generates spatial differentiation in morph frequencies, and selection and/or stochastic factors that lead to the coexistence of multiple morphs within a population. Thus, testing for both selection and stochastic factors is necessary for a comprehensive understanding of the mechanism underlying clinal variation in morph/allele frequency in natural populations. Here, I identified the evolutionary forces responsible for clinal variation of color morph frequency in Ischnura senegalensis by comparing the population divergence of putatively neutral loci generated by high-throughput next-generation sequencing (F _STn) with that of the putative color locus (F _STc). No strong correlation was observed between F _STn and F _STc, suggesting that stochastic factors contribute less to color-locus population divergence. F _STc was less than F _STn between populations exposed to similar environmental conditions, but greater than F _STn between populations exposed to different environmental conditions, suggesting that both balancing selection and divergent selection act on the color locus. Therefore, two antagonistic selection factors rather than stochastic and historical factors contribute to establishing the clinal variation of morph frequency in I. senegalensis.     

The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura

Disentangling the relative importance and potential interactions of selection and genetic drift in driving phenotypic divergence of species is a classical research topic in population genetics and evolutionary biology. Here, we evaluate the role of stochastic and selective forces on population divergence of a colour polymorphism in seven damselfly species of the genus Ischnura, with a particular focus on I. elegans and I. graellsii. Colour-morph frequencies in Spanish I. elegans populations varied greatly, even at a local scale, whereas more similar frequencies were found among populations in eastern Europe. In contrast, I. graellsii and the other five Ischnura species showed little variation in colour-morph frequencies between populations. F(ST)-outlier analyses revealed that the colour locus deviated strongly from neutral expectations in Spanish populations of I. elegans, contrasting the pattern found in eastern European populations, and in I. graellsii, where no such discrepancy between morph divergence and neutral divergence could be detected. This suggests that divergent selection has been operating on the colour locus in Spanish populations of I. elegans, whereas processes such as genetic drift, possibly in combination with other forms of selection (such as negative frequency-dependent selection), appear to have been present in other regions, such as eastern Europe. Overall, the results indicate that both selective and stochastic processes operate on these colour polymorphisms, and suggest that the relative importance of factors varies between geographical regions.     

THE INTERPLAY BETWEEN LOCAL ECOLOGY,DIVERGENT SELECTION,AND GENETIC DRIFT IN POPULATION DIVERGENCE OF A SEXUALLY ANTAGONISTIC FEMALE TRAIT

Genetically polymorphic species offer the possibility to study maintenance of genetic variation and the potential role for genetic drift in population divergence. Indirect inference of the selection regimes operating on polymorphic traits can be achieved by comparing population divergence in neutral genetic markers with population divergence in trait frequencies. Such an approach could further be combined with ecological data to better understand agents of selection. Here, we infer the selective regimes acting on a polymorphic mating trait in an insect group; the dorsal structures (either rough or smooth) of female diving beetles. Our recent work suggests that the rough structures have a sexually antagonistic function in reducing male mating attempts. For two species (Dytiscus lapponicus and Graphoderus zonatus), we could not reject genetic drift as an explanation for population divergence in morph frequencies, whereas for the third (Hygrotus impressopunctatus) we found that divergent selection pulls morph frequencies apart across populations. Furthermore, population morph frequencies in H. impressopunctatus were significantly related to local bioclimatic factors, providing an additional line of evidence for local adaptation in this species. These data, therefore, suggest that local ecological factors and sexual conflict interact over larger spatial scales to shape population divergence in the polymorphism.     

Detecting population structure in a high gene-flow species,Atlantic herring (Clupea harengus): direct,simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.     

Molecular evidence for selection on female color polymorphism in the damselfly Ischnura graellsii

The significance of female color polymorphism in Odonata remains controversial despite many field studies. The importance of random factors (founder effects, genetic drift and migration) versus selective forces for the maintenance of this polymorphism is still discussed. In this study, we specifically test whether the female color polymorphism of Ischnura graellsii (Odonata, Coenagrionidae) is under selection in the wild. We compared the degree of genetic differentiation based on RAPD markers (assumed to be neutral) with the degree of differentiation based on color alleles. Weir and Cockerham's theta values showed a significant degree of population differentiation for both sets of loci (RAPD and color alleles) but the estimated degree of population differentiation (theta) was significantly greater for the set of RAPD loci. This result shows that some sort of selection contributes to the maintenance of similar color morph frequencies across the studied populations. Our results combined with those of previous field studies suggest that at least in some I. graellsii populations, density-dependent mechanisms might help to prevent the loss of this polymorphism but cannot explain the similarity in morph frequencies among populations.     

Spatial variation in colour morph, spotting and allozyme frequencies in the candy-stripe spider, Enoplognatha ovata (Theridiidae) on two Swedish archipelagos

The selective significance, if any, of many invertebrate visible polymorphisms is still not fully understood. Here we examine colour- and black spotting-morph frequencies in the spider Enoplognatha ovata in populations on two Swedish archipelagos with respect to different spatial scales and, in one archipelago, against the background of variation at four putative neutral allozyme marker loci. Every population studied was polymorphic for colour and 28 out of 30 contained all three colour morphs – lineata, redimita and ovata. We found no evidence for a breakdown in the traditional colour morph designation previously suggested for other northern European populations of this species. For colour there is no significant heterogeneity at spatial scales greater than between local populations within islands. Black spotting frequencies show a similar lack of pattern over larger spatial scales except that there are significant differences between the Stockholm and Göteborg archipelagos. Measures of population differentiation (θ) within the Stockholm islands for the two visible systems show them to be significantly more differentiated than the neutral markers, suggesting local selection acting on them in a population-specific manner. On the basis of previous observations and the distribution of spotting phenotypes on a European scale, it is argued that thermal selection might operate on black spotting during the juvenile stages favouring more spots in continental climates. It is not clear what selective forces act on colour.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

When can ecological speciation be detected with neutral loci?

It is not yet clear under what conditions empirical studies can reliably detect progress toward ecological speciation through the analysis of allelic variation at neutral loci. We use a simulation approach to investigate the range of parameter space under which such detection is, and is not, likely. We specifically test for the conditions under which divergent natural selection can cause a ‘generalized barrier to gene flow’ that is present across the genome. Our individual‐based numerical simulations focus on how population divergence at neutral loci varies in relation to recombination rate with a selected locus, divergent selection on that locus, migration rate and population size. We specifically test whether genetic differences at neutral markers are greater between populations in different environments than between populations in similar environments. We find that this expected signature of ecological speciation can be detected under part of the parameter space, most consistently when divergent selection is strong and migration is intermediate. By contrast, the expected signature of ecological speciation is not reliably detected when divergent selection is weak or migration is low or high. These findings provide insights into the strengths and weaknesses of using neutral markers to infer ecological speciation in natural systems.     

Island biology and morphological divergence of the Skyros wall lizard <Emphasis Type="Italic">Podarcis gaigeae</Emphasis>: a combined role for local selection and genetic drift on color morph frequency divergence?

Background  

Patterns of spatial variation in discrete phenotypic traits can be used to draw inferences about the adaptive significance of traits and evolutionary processes, especially when compared to patterns of neutral genetic variation. Population divergence in adaptive traits such as color morphs can be influenced by both local ecology and stochastic factors such as genetic drift or founder events. Here, we use quantitative color measurements of males and females of Skyros wall lizard, Podarcis gaigeae, to demonstrate that this species is polymorphic with respect to throat color, and the morphs form discrete phenotypic clusters with limited overlap between categories. We use divergence in throat color morph frequencies and compare that to neutral genetic variation to infer the evolutionary processes acting on islet- and mainland populations.     

THE EFFICACY OF DIVERGENCE HITCHHIKING IN GENERATING GENOMIC ISLANDS DURING ECOLOGICAL SPECIATION

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal “divergence hitchhiking” models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome‐wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.     

SELECTION ON THE COLOR POLYMORPHISM IN HAWAIIAN HAPPY-FACE SPIDERS: EVIDENCE FROM GENETIC STRUCTURE AND TEMPORAL FLUCTUATIONS

Throughout this century genetic polymorphisms for color have been widely used as a research tool to allow insights into key evolutionary processes. Although color variants can often be diverse within populations, frequencies of different morphs may be similar across populations, either as a result of balancing selection or gene flow. Under these circumstances selection can be extremely difficult to demonstrate. Here we test for balancing selection on the naturally occurring color forms of the Hawaiian happy-face spider, Theridion grallator with two approaches. First, allozyme loci are used to generate a null model against which to test selection. Frequencies of alleles involved in the color polymorphism of T. grallator are used to generate another estimate for comparison. The results suggest that statistically similar frequencies of color morphs among populations of T. grallator may be maintained by some form of balancing selection. Second, we make use of an unusual event in which the normally stable frequencies of unpatterned and patterned morphs within a population were found to have shifted toward an excess of unpatterned morphs. We scored offspring of all fertilized, unpatterned (bottom-recessive) females found during this period of skewed morph frequencies and also in a year when morph frequencies were normal to deduce paternal color phenotypes. Mating was found to be random in the normal year, but in the perturbed year females had mated with rare (patterned) males twice as frequently as expected on the basis of the frequency of this morph type in the population. Both of these results are consistent with selection operating on the color polymorphism, and we speculate that apostatic selection, perhaps mediated by bird predators, may provide the mechanism.     

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Allele frequency divergence reveals ubiquitous influence of positive selection in Drosophila

Resolving the role of natural selection is a basic objective of evolutionary biology. It is generally difficult to detect the influence of selection because ubiquitous non-selective stochastic change in allele frequencies (genetic drift) degrades evidence of selection. As a result, selection scans typically only identify genomic regions that have undergone episodes of intense selection. Yet it seems likely such episodes are the exception; the norm is more likely to involve subtle, concurrent selective changes at a large number of loci. We develop a new theoretical approach that uncovers a previously undocumented genome-wide signature of selection in the collective divergence of allele frequencies over time. Applying our approach to temporally resolved allele frequency measurements from laboratory and wild Drosophila populations, we quantify the selective contribution to allele frequency divergence and find that selection has substantial effects on much of the genome. We further quantify the magnitude of the total selection coefficient (a measure of the combined effects of direct and linked selection) at a typical polymorphic locus, and find this to be large (of order 1%) even though most mutations are not directly under selection. We find that selective allele frequency divergence is substantially elevated at intermediate allele frequencies, which we argue is most parsimoniously explained by positive—not negative—selection. Thus, in these populations most mutations are far from evolving neutrally in the short term (tens of generations), including mutations with neutral fitness effects, and the result cannot be explained simply as an ongoing purging of deleterious mutations.     

Isolation by elevation: genetic structure at neutral and putatively non-neutral loci in a dominant tree of subtropical forests, Castanopsis eyrei

Background

The distribution of genetic diversity among plant populations growing along elevational gradients can be affected by neutral as well as selective processes. Molecular markers used to study these patterns usually target neutral processes only, but may also be affected by selection. In this study, the effects of elevation and successional stage on genetic diversity of a dominant tree species were investigated controlling for neutrality of the microsatellite loci used.

Methodology/Principal Findings

Diversity and differentiation among 24 populations of Castanopsis eyrei from different elevations (251–920 m) and successional stages were analysed by eight microsatellite loci. We found that one of the loci (Ccu97H18) strongly deviated from a neutral model of differentiation among populations due to either divergent selection or hitchhiking with an unknown selected locus. The analysis showed that C. eyrei populations had a high level of genetic diversity within populations (A_R = 7.6, H_E = 0.82). Genetic variation increased with elevation for both the putatively selected locus Ccu97H18 and the neutral loci. At locus Ccu97H18 one allele was dominant at low elevations, which was replaced at higher elevations by an increasing number of other alleles. The level of genetic differentiation at neutral loci was similar to that of other Fagaceae species (F_ST = 0.032,  = 0.15). Population differentiation followed a model of isolation by distance but additionally, strongly significant isolation by elevation was found, both for neutral loci and the putatively selected locus.

Conclusions/Significance

The results indicate higher gene flow among similar elevational levels than across different elevational levels and suggest a selective influence of elevation on the distribution of genetic diversity in C. eyrei. The study underlines the importance to check the selective neutrality of marker loci in analyses of population structure.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

Patterns of selection and allele diversity of class I and class II major histocompatibility loci across the species range of sockeye salmon (Oncorhynchus nerka)

The major histocompatibility complex (MHC), an important component of the vertebrate immune system, provides an important suite of genes to examine the role of genetic diversity at non‐neutral loci for population persistence. We contrasted patterns of diversity at the two classical MHC loci in sockeye salmon (Oncorhynchus nerka), MHC class I (UBA) and MHC class II (DAB), and neutral microsatellite loci across 70 populations spanning the species range from Washington State to Japan. There was no correlation in allelic richness or heterozygosity between MHC loci or between MHC loci and microsatellites. The two unlinked MHC loci may be responding to different selective pressures; the distribution of F_ST values for the two loci was uncorrelated, and evidence for both balancing and directional selection on alleles and lineages of DAB and UBA was observed in populations throughout the species range but rarely on both loci within a population. These results suggest that fluctuating selection has resulted in the divergence of MHC loci in contemporary populations.     

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT₅₀) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.     

Local selection and population structure in a deep‐sea fish,the roundnose grenadier (Coryphaenoides rupestris)

Local populations within a species can become isolated by stochastic or adaptive processes, though it is most commonly the former that we quantify. Using presumably neutral markers we can assess the time‐dependent process of genetic drift, and thereby quantify patterns of differentiation in support of the effective management of diversity. However, adaptive differences can be overlooked in these studies, and these are the very characteristics that we hope to conserve by managing neutral diversity. In this study, we used 16 hypothetically neutral microsatellite markers to investigate the genetic structure of the roundnose grenadier in the North Atlantic. We found that one locus was a clear outlier under directional selection, with F_ST values much greater than at the remaining loci. Differentiation between populations at this locus was related to depth, suggesting directional selection, presumably acting on a linked locus. Considering only the loci identified as neutral, there remained significant population structure over the region of the North Atlantic studied. In addition to a weak pattern of isolation by distance, we identified a putative barrier to gene flow between sample sites either side of the Charlie‐Gibbs Fracture Zone, which marks the location where the sub‐polar front crosses the Mid‐Atlantic Ridge. This may reflect a boundary across which larvae are differentially distributed in separate current systems to some extent, promoting differentiation by drift. Structure due to both drift and apparent selection should be considered in management policy.