Effect of migration and environmental heterogeneity on the maintenance of quantitative genetic variation: a simulation study

The paradox of high genetic variation observed in traits under stabilizing selection is a long‐standing problem in evolutionary theory, as mutation rates appear too low to explain observed levels of standing genetic variation under classic models of mutation–selection balance. Spatially or temporally heterogeneous environments can maintain more standing genetic variation within populations than homogeneous environments, but it is unclear whether such conditions can resolve the above discrepancy between theory and observation. Here, we use individual‐based simulations to explore the effect of various types of environmental heterogeneity on the maintenance of genetic variation (V_A) for a quantitative trait under stabilizing selection. We find that V_A is maximized at intermediate migration rates in spatially heterogeneous environments and that the observed patterns are robust to changes in population size. Spatial environmental heterogeneity increased variation by as much as 10‐fold over mutation–selection balance alone, whereas pure temporal environmental heterogeneity increased variance by only 45% at max. Our results show that some combinations of spatial heterogeneity and migration can maintain considerably more variation than mutation–selection balance, potentially reconciling the discrepancy between theoretical predictions and empirical observations. However, given the narrow regions of parameter space required for this effect, this is unlikely to provide a general explanation for the maintenance of variation. Nonetheless, our results suggest that habitat fragmentation may affect the maintenance of V_A and thereby reduce the adaptive capacity of populations.     

Environmental Heterogeneity Explains the Genetic Structure of Continental and Mediterranean Populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.     

Quantifying the climatic niche of symbiont partners in a lichen symbiosis indicates mutualist‐mediated niche expansions

The large distributional areas and ecological niches of many lichenized fungi may in part be due to the plasticity in interactions between the fungus (mycobiont) and its algal or cyanobacterial partners (photobionts). On the one hand, broad‐scale phylogenetic analyses show that partner compatibility in lichens is rather constrained and shaped by reciprocal selection pressures and codiversification independent of ecological drivers. On the other hand, sub‐species‐level associations among lichen symbionts appear to be environmentally structured rather than phylogenetically constrained. In particular, switching between photobiont ecotypes with distinct environmental preferences has been hypothesized as an adaptive strategy for lichen‐forming fungi to broaden their ecological niche. The extent and direction of photobiont‐mediated range expansions in lichens, however, have not been examined comprehensively at a broad geographic scale. Here we investigate the population genetic structure of Lasallia pustulata symbionts at sub‐species‐level resolution across the mycobiont's Europe‐wide range, using fungal MCM7 and algal ITS rDNA sequence markers. We show that variance in occurrence probabilities in the geographic distribution of genetic diversity in mycobiont‐photobiont interactions is closely related to changes in climatic niches. Quantification of niche extent and overlap based on species distribution modeling and construction of Hutchinsonian climatic hypervolumes revealed that combinations of fungal–algal interactions change at the sub‐species level along latitudinal temperature gradients and in Mediterranean climate zones. Our study provides evidence for symbiont‐mediated niche expansion in lichens. We discuss our results in the light of symbiont polymorphism and partner switching as potential mechanisms of environmental adaptation and niche evolution in mutualisms.     

Environmental heterogeneity generates opposite gene‐by‐environment interactions for two fitness‐related traits within a population

Theory predicts that environmental heterogeneity offers a potential solution to the maintenance of genetic variation within populations, but empirical evidence remains sparse. The live‐bearing fish Xiphophorus variatus exhibits polymorphism at a single locus, with different alleles resulting in up to five distinct melanistic “tailspot” patterns within populations. We investigated the effects of heterogeneity in two ubiquitous environmental variables (temperature and food availability) on two fitness‐related traits (upper thermal limits and body condition) in two different tailspot types (wild‐type and upper cut crescent). We found gene‐by‐environment (G × E) interactions between tailspot type and food level affecting upper thermal limits (UTL), as well as between tailspot type and thermal environment affecting body condition. Exploring mechanistic bases underlying these G × E patterns, we found no differences between tailspot types in hsp70 gene expression despite significant overall increases in expression under both thermal and food stress. Similarly, there was no difference in routine metabolic rates between the tailspot types. The reversal of relative performance of the two tailspot types under different environmental conditions revealed a mechanism by which environmental heterogeneity can balance polymorphism within populations through selection on different fitness‐related traits.     

Evolutionary history of Scinax treefrogs on land‐bridge islands in south‐eastern Brazil

Aim We investigated how Pleistocene refugia and recent (c. 12,000 years ago) sea level incursions shaped genetic differentiation in mainland and island populations of the Scinax perpusillus treefrog group. Location Brazilian Atlantic Forest, São Paulo state, south‐eastern Brazil. Methods Using mitochondrial and microsatellite loci, we examined population structure and genetic diversity in three species from the S. perpusillus group, sampled from three land‐bridge islands and five mainland populations, in order to understand the roles of Pleistocene forest fragmentation and sea level incursions on genetic differentiation. We calculated metrics of relatedness and genetic diversity to assess whether island populations exhibit signatures of genetic drift and isolation. Two of the three island populations in this study have previously been described as new species based on a combination of distinct morphological and behavioural characters, thus we used the molecular datasets to determine whether phenotypic change is consistent with genetic differentiation. Results Our analyses recovered three distinct lineages or demes composed of northern mainland São Paulo populations, southern mainland São Paulo populations, and one divergent island population. The two remaining island populations clustered with samples from adjacent mainland populations. Estimates of allelic richness were significantly lower, and estimates of relatedness were significantly higher, in island populations relative to their mainland counterparts. Main conclusions Fine‐scale genetic structure across mainland populations indicates the possible existence of local refugia within São Paulo state, underscoring the small geographic scale at which populations diverge in this species‐rich region of the Atlantic Coastal Forest. Variation in genetic signatures across the three islands indicates that the populations experienced different demographic processes after marine incursions fragmented the distribution of the S. perpusillus group. Genetic signatures of inbreeding and drift in some island populations indicate that small population sizes, coupled with strong ecological selection, may be important evolutionary forces driving speciation on land‐bridge islands.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.     

Inter and intra‐population phenotypic and genotypic structuring in the European whitefish Coregonus lavaretus,a rare freshwater fish in Scotland

This study revealed between‐lake genetic structuring between Coregonus lavaretus collected from the only two native populations of this species in Scotland, U.K. (Lochs Eck and Lomond) evidenced by the existence of private alleles (12 in Lomond and four in Eck) and significant genetic differentiation (F_ST = 0·056) across 10 microsatellite markers. Juvenile C. lavaretus originating from eggs collected from the two lakes and reared in a common‐garden experiment showed clear phenotypic differences in trophic morphology (i.e. head and body shape) between these populations indicating that these characteristics were, at least partly, inherited. Microsatellite analysis of adults collected from different geographic regions within Loch Lomond revealed detectable and statistically significant but relatively weak genetic structuring (F_ST = 0·001–0·024) and evidence of private alleles related to the basin structure of the lake. Within‐lake genetic divergence patterns suggest three possibilities for this observed pattern: (1) differential selection pressures causing divergence into separate gene pools, (2) a collapse of two formerly divergent gene pools and (3) a stable state maintained by balancing selection forces resulting from spatial variation in selection and lake heterogeneity. Small estimates of effective population sizes for the populations in both lakes suggest that the capacity of both populations to adapt to future environmental change may be limited.     

Geography is more important than host plant use for the population genetic structure of a generalist insect herbivore

Population divergence can occur due to mechanisms associated with geographic isolation and/or due to selection associated with different ecological niches. Much of the evidence for selection‐driven speciation has come from studies of specialist insect herbivores that use different host plant species; however, the influence of host plant use on population divergence of generalist herbivores remains poorly understood. We tested how diet breadth, host plant species and geographic distance influence population divergence of the fall webworm (Hyphantria cunea; FW). FW is a broadly distributed, extreme generalist herbivore consisting of two morphotypes that have been argued to represent two different species: black‐headed and red‐headed. We characterized the differentiation of FW populations at two geographic scales. We first analysed the influence of host plant and geographic distance on genetic divergence across a broad continental scale for both colour types. We further analysed the influence of host plant, diet breadth and geographic distance on divergence at a finer geographic scale focusing on red‐headed FW in Colorado. We found clear genetic and morphological distinction between red‐ and black‐headed FW, and Colorado FW formed a genetic cluster distinct from other locations. Although both geographic distance and host plant use were correlated with genetic distance, geographic distance accounted for up to 3× more variation in genetic distance than did host plant use. As a rare study investigating the genetic structure of a widespread generalist herbivore over a broad geographic range (up to 3,000 km), our study supports a strong role for geographic isolation in divergence in this system.     

Genomic differentiation between temperate and tropical Australian populations of Drosophila melanogaster

Determining the genetic basis of environmental adaptation is a central problem of evolutionary biology. This issue has been fruitfully addressed by examining genetic differentiation between populations that are recently separated and/or experience high rates of gene flow. A good example of this approach is the decades-long investigation of selection acting along latitudinal clines in Drosophila melanogaster. Here we use next-generation genome sequencing to reexamine the well-studied Australian D. melanogaster cline. We find evidence for extensive differentiation between temperate and tropical populations, with regulatory regions and unannotated regions showing particularly high levels of differentiation. Although the physical genomic scale of geographic differentiation is small--on the order of gene sized--we observed several larger highly differentiated regions. The region spanned by the cosmopolitan inversion polymorphism In(3R)P shows higher levels of differentiation, consistent with the major difference in allele frequencies of Standard and In(3R)P karyotypes in temperate vs. tropical Australian populations. Our analysis reveals evidence for spatially varying selection on a number of key biological processes, suggesting fundamental biological differences between flies from these two geographic regions.     

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long‐standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter‐population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.     

Genetically informed ecological niche models improve climate change predictions

We examined the hypothesis that ecological niche models (ENMs) more accurately predict species distributions when they incorporate information on population genetic structure, and concomitantly, local adaptation. Local adaptation is common in species that span a range of environmental gradients (e.g., soils and climate). Moreover, common garden studies have demonstrated a covariance between neutral markers and functional traits associated with a species’ ability to adapt to environmental change. We therefore predicted that genetically distinct populations would respond differently to climate change, resulting in predicted distributions with little overlap. To test whether genetic information improves our ability to predict a species’ niche space, we created genetically informed ecological niche models (gENMs) using Populus fremontii (Salicaceae), a widespread tree species in which prior common garden experiments demonstrate strong evidence for local adaptation. Four major findings emerged: (i) gENMs predicted population occurrences with up to 12‐fold greater accuracy than models without genetic information; (ii) tests of niche similarity revealed that three ecotypes, identified on the basis of neutral genetic markers and locally adapted populations, are associated with differences in climate; (iii) our forecasts indicate that ongoing climate change will likely shift these ecotypes further apart in geographic space, resulting in greater niche divergence; (iv) ecotypes that currently exhibit the largest geographic distribution and niche breadth appear to be buffered the most from climate change. As diverse agents of selection shape genetic variability and structure within species, we argue that gENMs will lead to more accurate predictions of species distributions under climate change.     

The experimental evolution of specialists,generalists, and the maintenance of diversity

Environmental heterogeneity may be a general explanation for both the quantity of genetic variation in populations and the ecological niche width of individuals. To evaluate this hypothesis, I review the literature on selection experiments in heterogeneous environments. The niche width usually – but not invariably – evolves to match the amount of environmental variation, specialists evolving in homogeneous environments and generalists evolving in heterogeneous environments. The genetics of niche width are more complex than has previously been recognized, particularly with respect to the magnitude of costs of adaptation and the putative constraints on the evolution of generalists. Genetic variation in fitness is more readily maintained in heterogeneous environments than in homogeneous environments and this diversity is often stably maintained through negative frequency‐dependent selection. Moreover environmental heterogeneity appears to be a plausible mechanism for at least two well‐known patterns of species diversity at the landscape scale. I conclude that environmental heterogeneity is a plausible and possibly very general explanation for diversity across the range of scales from individuals to landscapes.     

Genomics,environment and balancing selection in behaviourally bimodal populations: The caribou case

Selection forces that favour different phenotypes in different environments can change frequencies of genes between populations along environmental clines. Clines are also compatible with balancing forces, such as negative frequency‐dependent selection (NFDS), which maintains phenotypic polymorphisms within populations. For example, NFDS is hypothesized to maintain partial migration, a dimorphic behavioural trait prominent in species where only a fraction of the population seasonally migrates. Overall, NFDS is believed to be a common phenomenon in nature, yet a scarcity of studies were published linking naturally occurring allelic variation with bimodal or multimodal phenotypes and balancing selection. We applied a Pool‐seq approach and detected selection on alleles associated with environmental variables along a North–South gradient in western North American caribou, a species displaying partially migratory behaviour. On 51 loci, we found a signature of balancing selection, which could be related to NFDS and ultimately the maintenance of the phenotypic polymorphisms known within these populations. Yet, remarkably, we detected directional selection on a locus when our sample was divided into two behaviourally distinctive groups regardless of geographic provenance (a subset of GPS‐collared migratory or sedentary individuals), indicating that, within populations, phenotypically homogeneous groups were genetically distinctive. Loci under selection were linked to functional genes involved in oxidative stress response, body development and taste perception. Overall, results indicated genetic differentiation along an environmental gradient of caribou populations, which we found characterized by genes potentially undergoing balancing selection. We suggest that the underlining balancing force, NFDS, plays a strong role within populations harbouring multiple haplotypes and phenotypes, as it is the norm in animals, plants and humans too.     

Woody climbers show greater population genetic differentiation than trees: Insights into the link between ecological traits and diversification

The climbing habit is a key innovation in plants: climbing taxa have higher species richness than nonclimbing sister groups. We evaluated the hypothesis that climbing plant species show greater among‐population genetic differentiation than nonclimber species. We compared the among‐population genetic distance in woody climbers (eight species, 30 populations) and trees (seven species, 29 populations) coexisting in nine communities in a temperate rainforest. We also compared within‐population genetic diversity in co‐occurring woody climbers and trees in two communities. Mean genetic distance between populations of climbers was twice that of trees. Isolation by distance (increase in genetic distance with geographic distance) was greater for climbers. Climbers and trees showed similar within‐population genetic diversity. Our longevity estimate suggested that climbers had shorter generation times, while other biological features often associated with diversification (dispersal and pollination syndromes, mating system, size, and metabolic rate) did not show significant differences between groups. We hypothesize that the greater population differentiation in climbers could result from greater evolutionary responses to local selection acting on initially higher within‐population genetic diversity, which could be driven by neutral processes associated with shorter generation times. Increased population genetic differentiation could be incorporated as another line of evidence when testing for key innovations.     

Genetic attributes of midwife toad (Alytes obstetricans) populations do not correlate with degree of species decline

Genetic diversity is crucial for long‐term population persistence. Population loss and subsequent reduction in migration rate among the most important processes that are expected to lead to a reduction in genetic diversity and an increase in genetic differentiation. While the theory behind this is well‐developed, empirical evidence from wild populations is inconsistent. Using microsatellite markers, we compared the genetic structure of populations of an amphibian species, the midwife toad (Alytes obstetricans), in four Swiss regions where the species has suffered variable levels of subpopulation extirpation. We also quantified the effects of several geographic factors on genetic structure and used a model selection approach to ascertain which of the variables were important for explaining genetic variation. Although subpopulation pairwise F_ST‐values were highly significant even over small geographic scales, neither any of the geographic variables nor loss of subpopulations were important factors for predicting spatial genetic structure. The absence of a signature of subpopulation loss on genetic differentiation may suggest that midwife toad subpopulations function as relatively independent units.     

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.     

THE EVOLUTION OF ENVIRONMENTAL TOLERANCE AND RANGE SIZE: A COMPARISON OF GEOGRAPHICALLY RESTRICTED AND WIDESPREAD MIMULUS

The geographic ranges of closely related species can vary dramatically, yet we do not fully grasp the mechanisms underlying such variation. The niche breadth hypothesis posits that species that have evolved broad environmental tolerances can achieve larger geographic ranges than species with narrow environmental tolerances. In turn, plasticity and genetic variation in ecologically important traits and adaptation to environmentally variable areas can facilitate the evolution of broad environmental tolerance. We used five pairs of western North American monkeyflowers to experimentally test these ideas by quantifying performance across eight temperature regimes. In four species pairs, species with broader thermal tolerances had larger geographic ranges, supporting the niche breadth hypothesis. As predicted, species with broader thermal tolerances also had more within‐population genetic variation in thermal reaction norms and experienced greater thermal variation across their geographic ranges than species with narrow thermal tolerances. Species with narrow thermal tolerance may be particularly vulnerable to changing climatic conditions due to lack of plasticity and insufficient genetic variation to respond to novel selection pressures. Conversely, species experiencing high variation in temperature across their ranges may be buffered against extinction due to climatic changes because they have evolved tolerance to a broad range of temperatures.