Lack of genetic variation prevents adaptation at the geographic range margin in a damselfly

What limits a species' distribution in the absence of physical barriers? Genetic load due to asymmetric gene flow and the absence of genetic variation due to lack of gene flow are hypothesized to constrain adaptation to novel environments in marginal populations, preventing range expansion. Here, we examined the genetic structure and geographic variation in morphological traits in two damselflies (Ischnura asiatica and I. senegalensis) along a latitudinal gradient in Japan, which is the distribution centre of I. asiatica and the northern limit of I. senegalensis. Genomewide genetic analyses found a loss of genetic diversity at the edge of distribution in I. senegalensis but consistently high diversity in I. asiatica. Gene flow was asymmetric in a south–north direction in both species. Although body size and wing loading showed decreasing latitudinal clines (smaller in north) in I. asiatica in Japan, increasing latitudinal clines (larger in north) in these phenotypic markers were observed in I. senegalensis, particularly near the northern boundary, which coincided well with the location where genetic diversity began a sharp decline. In ectothermic animals, increasing latitudinal cline in these traits was suggested to be established when they failed to adapt to thermal gradient. Therefore, our findings support the possibility that a lack of genetic variation rather than geneflow swamping is responsible for the constraint of adaptation at the margin of geographic distribution.     

Geographical patterns of adaptation within a species' range: interactions between drift and gene flow

We use individual-based stochastic simulations and analytical deterministic predictions to investigate the interaction between drift, natural selection and gene flow on the patterns of local adaptation across a fragmented species' range under clinally varying selection. Migration between populations follows a stepping-stone pattern and density decreases from the centre to the periphery of the range. Increased migration worsens gene swamping in small marginal populations but mitigates the effect of drift by replenishing genetic variance and helping purge deleterious mutations. Contrary to the deterministic prediction that increased connectivity within the range always inhibits local adaptation, simulations show that low intermediate migration rates improve fitness in marginal populations and attenuate fitness heterogeneity across the range. Such migration rates are optimal in that they maximize the total mean fitness at the scale of the range. Optimal migration rates increase with shallower environmental gradients, smaller marginal populations and higher mutation rates affecting fitness.     

Evidence for adaptation from standing genetic variation on an antimicrobial peptide gene in the mussel Mytilus edulis

Genome scans of population differentiation identify candidate loci for adaptation but provide little information on how selection has influenced the genetic structure of these loci. Following a genome scan, we investigated the nature of the selection responsible for the outlying differentiation observed between populations of the marine mussel Mytilus edulis at a leucine/arginine polymorphism (L31R) in the antimicrobial peptide MGD2. We analysed DNA sequence polymorphisms, allele frequencies and population differentiation of polymorphisms closely linked to L31R, and pairwise and third‐order linkage disequilibria. An outlying level of population differentiation was observed at L31R only, while no departure from panmixia was observed at linked loci surrounding L31R, as in most of the genome. Selection therefore seems to affect L31R directly. Three hypotheses can explain the lack of differentiation in the chromosomal region close to L31R: (i) hitchhiking has occurred but migration and recombination subsequently erased the signal, (ii) selection was weak enough and recombination strong enough to limit the hitchhiking effect to a very small chromosomal region or (iii) selection acted on a pre‐existing polymorphism (i.e. standing variation) at linkage equilibrium with its background. Linkage equilibrium was observed between L31R and linked polymorphisms in every population analysed, as expected under the three hypotheses. However, linkage disequilibrium was observed in some populations between pairs of loci located upstream and downstream to L31R, generating a complex pattern of third‐order linkage disequilibria which is best explained by the hypothesis of selection on a pre‐existing polymorphism. We hypothesise that selection could be either balanced, maintaining alleles at different frequencies depending on the pathogen community encountered locally by mussels, or intermittent, resulting in sporadic fluctuations in allele frequency.     

Quantifying the effects of migration and mutation on adaptation and demography in spatially heterogeneous environments

How do mutation and gene flow influence population persistence, niche expansion and local adaptation in spatially heterogeneous environments? In this article, we analyse a demographic and evolutionary model of adaptation to an environment containing two habitats in equal frequencies, and we bridge the gap between different theoretical frameworks. Qualitatively, our model yields four qualitative types of outcomes: (i) global extinction of the population, (ii) adaptation to one habitat only, but also adaptation to both habitats with, (iii) specialized phenotypes or (iv) with generalized phenotypes, and we determine the conditions under which each equilibrium is reached. We derive new analytical approximations for the local densities and the distributions of traits in each habitat under a migration–selection–mutation balance, compute the equilibrium values of the means, variances and asymmetries of the local distributions of phenotypes, and contrast the effects of migration and mutation on the evolutionary outcome. We then check our analytical results by solving our model numerically, and also assess their robustness in the presence of demographic stochasticity. Although increased migration results in a decrease in local adaptation, mutation in our model does not influence the values of the local mean traits. Yet, both migration and mutation can have dramatic effects on population size and even lead to metapopulation extinction when selection is strong. Niche expansion, the ability for the population to adapt to both habitats, can also be prevented by small migration rates and a reduced evolutionary potential characterized by rare mutation events of small effects; however, niche expansion is otherwise the most likely outcome. Although our results are derived under the assumption of clonal reproduction, we finally show and discuss the links between our model and previous quantitative genetics models.     

The genetics of adaptation to discrete heterogeneous environments: frequent mutation or large‐effect alleles can allow range expansion

Range expansions are complex evolutionary and ecological processes. From an evolutionary standpoint, a populations' adaptive capacity can determine the success or failure of expansion. Using individual‐based simulations, we model range expansion over a two‐dimensional, approximately continuous landscape. We investigate the ability of populations to adapt across patchy environmental gradients and examine how the effect sizes of mutations influence the ability to adapt to novel environments during range expansion. We find that genetic architecture and landscape patchiness both have the ability to change the outcome of adaptation and expansion over the landscape. Adaptation to new environments succeeds via many mutations of small effect or few of large effect, but not via the intermediate between these cases. Higher genetic variance contributes to increased ability to adapt, but an alternative route of successful adaptation can proceed from low genetic variance scenarios with alleles of sufficiently large effect. Steeper environmental gradients can prevent adaptation and range expansion on both linear and patchy landscapes. When the landscape is partitioned into local patches with sharp changes in phenotypic optimum, the local magnitude of change between subsequent patches in the environment determines the success of adaptation to new patches during expansion.     

Effects of founding genetic variation on adaptation to a novel resource

Population genetic theory predicts that adaptation in novel environments is enhanced by genetic variation for fitness. However, theory also predicts that under strong selection, demographic stochasticity can drive populations to extinction before they can adapt. We exposed wheat-adapted populations of the flour beetle (Tribolium castaneum) to a novel suboptimal corn resource, to test the effects of founding genetic variation on population decline and subsequent extinction or adaptation. As previously reported, genetically diverse populations were less likely to go extinct. Here, we show that among surviving populations, genetically diverse groups recovered faster after the initial population decline. Within two years, surviving populations significantly increased their fitness on corn via increased fecundity, increased egg survival, faster larval development, and higher rate of egg cannibalism. However, founding genetic variation only enhanced the increase in fecundity, despite existing genetic variation-and apparent lack of trade-offs-for egg survival and larval development time. Thus, during adaptation to novel habitats the positive impact of genetic variation may be restricted to only a few traits, although change in many life-history traits may be necessary to avoid extinction. Despite severe initial maladaptation and low population size, genetic diversity can thus overcome the predicted high extinction risk in new habitats.     

Evidence for evolutionary change associated with the recent range expansion of the British butterfly, Aricia agestis, in response to climate change

Poleward range expansions are widespread responses to recent climate change and are crucial for the future persistence of many species. However, evolutionary change in traits such as colonization history and habitat preference may also be necessary to track environmental change across a fragmented landscape. Understanding the likelihood and speed of such adaptive change is important in determining the rate of species extinction with ongoing climate change. We conducted an amplified fragment length polymorphism (AFLP)‐based genome scan across the recently expanded UK range of the Brown Argus butterfly, Aricia agestis, and used outlier‐based (DFDIST and BayeScan) and association‐based (Isolation‐By‐Adaptation) statistical approaches to identify signatures of evolutionary change associated with range expansion and habitat use. We present evidence for (i) limited effects of range expansion on population genetic structure and (ii) strong signatures of selection at approximately 5% AFLP loci associated with both the poleward range expansion of A. agestis and differences in habitat use across long‐established and recently colonized sites. Patterns of allele frequency variation at these candidate loci suggest that adaptation to new habitats at the range margin has involved selection on genetic variation in habitat use found across the long‐established part of the range. Our results suggest that evolutionary change is likely to affect species’ responses to climate change and that genetic variation in ecological traits across species’ distributions should be maximized to facilitate range shifts across a fragmented landscape, particularly in species that show strong associations with particular habitats.     

The consequences of dominance and gene flow for local adaptation and differentiation at two linked loci

For a subdivided population the consequences of dominance and gene flow for the maintenance of multilocus polymorphism, local adaptation, and differentiation are investigated. The dispersing population inhabits two demes in which selection acts in opposite direction. Fitness is determined additively by two linked diallelic loci with arbitrary intermediate dominance (no over- or underdominance). For weak as well as strong migration, the equilibrium structure is derived. As a special case, a continuous-time continent–island model (CI model) is analyzed, with one-way migration from the continent to the island. For this CI model, the equilibrium and stability configuration is obtained explicitly for weak migration, for strong migration, for independent loci, and for complete linkage. For independent loci, the possible bifurcation patterns are derived as functions of the migration rate. These patterns depend strongly on the degree of dominance. The effects of dominance, linkage, and migration on the amount of linkage disequilibrium (LD) and the degree of local adaptation are explored. Explicit formulas are obtained for $D$   $(=x_1{x}_4-x_2{x}_3)$ and $r^2$ (the squared correlation in allelic state). They demonstrate that dominant island alleles increase $D$ and decrease $r^2$. Local adaptation is elevated by dominance of the locally adaptive alleles. The effective migration rate at a linked neutral locus is calculated. If advantageous alleles are dominant, it is decreased only slightly below the actual migration rate. For a quantitative trait that is determined by two additive loci, the influence of dominance on measures of differentiation is studied. Explicit expressions for $Q_{ST}$ and two types of $F_{ST}$ at equilibrium are deduced and their relation is discussed.     

Triangulating the genetic basis of adaptation to multifarious selection

Understanding how natural populations adapt to their local environments is a major research theme for ecological genomics. This endeavour begins by sleuthing for shared genetic similarities among unrelated natural populations sharing adaptive traits to documented selective pressures. When the selective pressures have low dimensionality, and the genetic response is localized to a few genes of major effect, this detective work is relatively straightforward. However, in the real world, populations face a complex mixture of selective pressures and many adaptive responses are the result of changes in quantitative traits that have a polygenic genetic basis. This complex relationship between environment and adaptation presents a significant challenge. How can we begin to identify drivers of adaptation in natural settings? In this issue of Molecular Ecology, Orsini et al. (2012) take advantage of the biological attributes of the freshwater microcrustacean Daphnia ( Fig. 1 ) to disentangle multidimensional selection’s signature on the genome of populations that have repeatedly evolved adaptive responses to isolated selective pressures including predation, parasitism and anthropogenic changes in land use. Orsini et al. (2012) leverage a powerful combination of spatially structured populations in a geographic mosaic of environmental stressors, the historical archive of past genotypes preserved in lake‐bottom sediments and selection experiments to identify sets of candidate genomic regions associated with adaptation in response to these three environmental stressors. This study provides a template for future investigation in ecological genomics, combining multiple experimental approaches with the genomic investigation of a well‐studied ecological model species.

Figure 1 Open in figure viewer PowerPoint Adult Daphnia magna carrying a resting egg in the brood pouch. The water flea Daphnia is a renowned ecological model system and rapidly developing as an ecological and environmental genomics model species. Photo credit Joachim Mergeay.     

Mate choice evolution, dominance effects, and the maintenance of genetic variation

Female mate choice influences the maintenance of genetic variation by altering the mating success of males with different genotypes. The evolution of preferences themselves, on the other hand, depends on genetic variation present in the population. Few models have tracked this feedback between a choice gene and its effects on genetic variation, in particular when genes that determine offspring viability and attractiveness have dominance effects. Here we build a population genetic model that allows comparing the evolution of various choice rules in a single framework. We first consider preferences for good genes and show that focused preferences for homozygotes evolve more easily than broad preferences, which allow heterozygous males high mating success too. This occurs despite better maintenance of genetic diversity in the latter scenario, and we discuss why empirical findings of superior mating success of heterozygous males consequently do not immediately lead to a better understanding of the lek paradox. Our results thus suggest that the mechanisms that help maintain genetic diversity also have a flipside of making female choice an inaccurate means of producing the desired kind of offspring. We then consider preferences for heterozygosity per se, and show that these evolve only under very special conditions. Choice for compatible genotypes can evolve but its selective advantage diminishes quickly due to frequency-dependent selection. Finally, we show that our model reproduces earlier results on selfing, when the female choice strategy produces assortative mating. Overall, our model indicates that various forms of heterozygote-favouring (or variable) female choice pose a problem for the theory of sexual ornamentation based on indirect benefits, rather than a solution.     

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.     

Recency, range expansion, and unsorted lineages: implications for interpreting neutral genetic variation in the sharp-tailed grouse (Tympanuchus phasianellus)

Both current and historical patterns of variation are relevant to understanding and managing ecological diversity. Recently derived species present a challenge to the reconstruction of historical patterns because neutral molecular data for these taxa are more likely to exhibit effects of recent and ongoing demographic processes. We studied geographical patterns of neutral molecular variation in a species thought to be of relatively recent origin, Tympanuchus phasianellus (sharp-tailed grouse), using mitochondrial control region sequences (CR-I), amplified fragment length polymorphisms (AFLP), and microsatellites. For historical context, we also analysed CR-I in all species of Tympanuchus. Within T. phasianellus, we found evidence for restricted gene flow between eastern and western portions of the species range, generally corresponding with the range boundary of T. p. columbianus and T. p. jamesi. The mismatch distribution and molecular clock estimates from the CR-I data suggested that all Tympanuchus underwent a range expansion prior to sorting of mitotypes among the species, and that sorting may have been delayed as a result of mutation-drift disequilibrium. This study illustrates the challenge of using genetic data to detect historical divergence in groups that are of relatively recent origin, or that have a history dominated by nonequilibrium conditions. We suggest that in such cases, morphological, ecological, and behavioural data may be particularly important adjuncts to molecular data for the recognition of historically or adaptively divergent groups.     

Landscape genetics and hierarchical genetic structure in Atlantic salmon: the interaction of gene flow and local adaptation

Disentangling evolutionary forces that may interact to determine the patterns of genetic differentiation within and among wild populations is a major challenge in evolutionary biology. The objective of this study was to assess the genetic structure and the potential influence of several ecological variables on the extent of genetic differentiation at multiple spatial scales in a widely distributed species, the Atlantic salmon, Salmo salar . A total of 2775 anadromous fish were sampled from 51 rivers along the North American Atlantic coast and were genotyped using 13 microsatellites. A Bayesian analysis clustered these populations into seven genetically and geographically distinct groups, characterized by different environmental and ecological factors, mainly temperature. These groups were also characterized by different extent of genetic differentiation among populations. Dispersal was relatively high and of the same magnitude within compared to among regional groups, which contrasted with the maintenance of a regional genetic structure. However, genetic differentiation was lower among populations exchanging similar rates of local as opposed to inter-regional migrants, over the same geographical scale. This raised the hypothesis that gene flow could be constrained by local adaptation at the regional scale. Both coastal distance and temperature regime were found to influence the observed genetic structure according to landscape genetic analyses. The influence of other factors such as latitude, river length and altitude, migration tactic, and stocking was not significant at any spatial scale. Overall, these results suggested that the interaction between gene flow and thermal regime adaptation mainly explained the hierarchical genetic structure observed among Atlantic salmon populations.     

The effects of migration load,selfing, inbreeding depression,and the genetics of adaptation on autotetraploid versus diploid establishment in peripheral habitats

The distribution and abundance of polyploids has intrigued biologists since their discovery in the early 20th century. A pattern in nature that may give insight to processes that shape the distribution and abundance of polyploids is that polyploid populations are sometimes associated with peripheral habitats within the range of a species of mixed ploidy. Here, adaptation and competition of a diploid versus an autotetraploid population in a peripheral habitat are examined theoretically. It is shown that a nascent autotetraploid population adapts to and outcompetes a diploid population in the periphery when the rate of gamete dispersal is high, and when the mode of gene action is recessive for moderate to high rates of selfing. With additive or dominant modes of gene action, the conditions for an autotetraploid to outcompete a diploid in the periphery appear determined more by the rate of selfing and less by gamete dispersal. All of these results are based on empirical work that suggests inbreeding depression is higher in diploids versus autotetraploids. Generally, the results indicate that, although autotetraploids incur minority cytotype exclusion, diploids face burdens themselves. In the case of adaptation to a peripheral habitat, this burden is migration load from gamete and propagule dispersal.     

Indirect genetic effects and the genetic bases of social dominance: evidence from cattle

Genetic studies of social behaviour have currently received new impetus from models including indirect genetic effects (IGEs) of social partners. This study aimed at investigating the contribution of conspecifics in social dominance, considered as response of dyadic interaction that is, winning (dominant individual) or losing (subordinate). A genetic correlation of −1 is expected between the attitude to win and the attitude to loose, and because a population always accounts for half winners and half losers, the heritability of the dominant status should be close to zero. Specifically, social dominance was studied in Aosta Chestnut and Aosta Black Pied (Bos taurus) breeds, alpine rustic cattle famous for traditional tournaments where pairs of cows assess dominant status in bloodless fights. The outcomes of 25 590 dyadic interactions performed by 8159 individuals in 11 years were analysed by applying a classical quantitative model and models including indirect effects. Data were analysed via Bayesian approach on a threshold trait. The assessment of variances revealed a genetic correlation of −0.976 between direct and indirect genetic components. The heritability measured on a liability scale was 0.122 for direct phenotype, but decreased to 0.014 when the total heritable variance (TBV) was considered. The trend of estimated breeding values showed that the total TBV was constant over the years, even though its direct component increased and the indirect part decreased. This result confirms the relevance of IGEs on social behaviour and the assumption that the mean individual social dominance cannot evolve within a population, due to the evolutionary constraints imposed by the ‘social environment''.     

Cryptic genetic variation and adaptation to waterlogging in Caledonian Scots pine,Pinus sylvestris L.

Local adaptation occurs as the result of differential selection among populations. Observations made under common environmental conditions may reveal phenotypic differences between populations with an underlying genetic basis; however, exposure to a contrasting novel environment can trigger release of otherwise unobservable (cryptic) genetic variation. We conducted a waterlogging experiment on a common garden trial of Scots pine, Pinus sylvestris (L.), saplings originating from across a steep rainfall gradient in Scotland. A flood treatment was maintained for approximately 1 year; physiological responses were gauged periodically in terms of photochemical capacity as measured via chlorophyll fluorescence. During the treatment, flooded individuals experienced a reduction in photochemical capacity, F_v/F_m, this reduction being greater for material originating from drier, eastern sites. Phenotypic variance was increased under flooding, and this increase was notably smaller in saplings originating from western sites where precipitation is substantially greater and waterlogging is more common. We conclude that local adaptation has occurred with respect to waterlogging tolerance and that, under the flooding treatment, the greater increase in variability observed in populations originating from drier sites is likely to reflect a relative absence of past selection. In view of a changing climate, we note that comparatively maladapted populations may possess considerable adaptive potential, due to cryptic genetic variation, that should not be overlooked.     

Inferring local adaptation from QST-FST comparisons: neutral genetic and quantitative trait variation in European populations of great snipe

We applied a phenotypic QST (PST) vs. FST approach to study spatial variation in selection among great snipe (Gallinago media) populations in two regions of northern Europe. Morphological divergence between regions was high despite low differentiation in selectively neutral genetic markers, whereas populations within regions showed very little neutral divergence and trait differentiation. QST > FST was robust against altering assumptions about the additive genetic proportions of variance components. The homogenizing effect of gene flow (or a short time available for neutral divergence) has apparently been effectively counterbalanced by differential natural selection, although one trait showed some evidence of being under uniform stabilizing selection. Neutral markers can hence be misleading for identifying evolutionary significant units, and adopting the PST-FST approach might therefore be valuable when common garden experiments is not an option. We discuss the statistical difficulties of documenting uniform selection as opposed to divergent selection, and the need for estimating measurement error. Instead of only comparing overall QST and FST values, we advocate the use of partial matrix permutation tests to analyse pairwise QST differences among populations, while statistically controlling for neutral differentiation.