Unravelling the effects of contemporary and historical range expansion on the distribution of genetic diversity in the damselfly Coenagrion scitulum

Although genetic diversity provides the basic substrate for evolution, there are a limited number of studies that assess the impact of recent climate change on intraspecific genetic variation. This study aims to unravel the degree to which historical and contemporary factors shape genetic diversity and structure across a large part of the range of the range‐expanding damselfly Coenagrion scitulum (Rambur, 1842). A total of 525 individuals from 31 populations were genotyped at nine microsatellites, and a subset was sequenced at two mitochondrial genes. We inferred the importance of geography, environmental factors, and recent range expansion on genetic diversity and structure. Genetic diversity decreased going westwards, suggesting a signature of historical post‐glacial expansion from east to west and the presence of eastern refugia. Although genetic differentiation decreased going northwards, it increased in the northern edge populations, suggesting a role of contemporary range expansion on the genetic make‐up of populations. The phylogeographical context was proven to be essential in understanding and identifying the genetic signatures of local contemporary processes. Within this framework, our results highlight that recent range expansion of a good disperser can decrease genetic diversity and increase genetic differentiation which should be considered when devising suitable conservation strategies.     

Rapid range expansion increases genetic differentiation while causing limited reduction in genetic diversity in a damselfly

Many ectothermic species are currently expanding their geographic range due to global warming. This can modify the population genetic diversity and structure of these species because of genetic drift during the colonization of new areas. Although the genetic signatures of historical range expansions have been investigated in an array of species, the genetic consequences of natural, contemporary range expansions have received little attention, with the only studies available focusing on range expansions along a narrow front. We investigate the genetic consequences of a natural range expansion in the Mediterranean damselfly Coenagrion scitulum, which is currently rapidly expanding along a broad front in different directions. We assessed genetic diversity and genetic structure using 12 microsatellite markers in five centrally located populations and five recently established populations at the edge of the geographic distribution. Our results suggest that, although a marginal significant decrease in the allelic richness was found in the edge populations, genetic diversity has been preserved during the range expansion of this species. Nevertheless, edge populations were genetically more differentiated compared with core populations, suggesting genetic drift during the range expansion. The smaller effective population sizes of the edge populations compared with central populations also suggest a contribution of genetic drift after colonization. We argue and document that range expansion along multiple axes of a broad expansion front generates little reduction in genetic diversity, yet stronger differentiation of the edge populations.     

Local adaptation at range edges: comparing elevation and latitudinal gradients

Local adaptation at range edges influences species’ distributions and how they respond to environmental change. However, the factors that affect adaptation, including gene flow and local selection pressures, are likely to vary across different types of range edge. We performed a reciprocal transplant experiment to investigate local adaptation in populations of Plantago lanceolata and P. major from central locations in their European range and from their latitudinal and elevation range edges (in northern Scandinavia and Swiss Alps, respectively). We also characterized patterns of genetic diversity and differentiation in populations using molecular markers. Range‐centre plants of P. major were adapted to conditions at the range centre, but performed similarly to range‐edge plants when grown at the range edges. There was no evidence for local adaptation when comparing central and edge populations of P. lanceolata. However, plants of both species from high elevation were locally adapted when compared with plants from high latitude, although the reverse was not true. This asymmetry was associated with greater genetic diversity and less genetic differentiation over the elevation gradient than over the latitudinal gradient. Our results suggest that adaptation in some range‐edge populations could increase their performance following climate change. However, responses are likely to differ along elevation and latitudinal gradients, with adaptation more likely at high‐elevation. Furthermore, based upon these results, we suggest that gene flow is unlikely to constrain adaptation in range‐edge populations of these species.     

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.     

Crossing the threshold: gene flow,dominance and the critical level of standing genetic variation required for adaptation to novel environments

Genetic architecture plays an important role in the process of adaptation to novel environments. One example is the role of allelic dominance, where advantageous recessive mutations have a lower probability of fixation than advantageous dominant mutations. This classic observation, termed ‘Haldane's sieve’, has been well explored theoretically for single isolated populations adapting to new selective regimes. However, the role of dominance is less well understood for peripheral populations adapting to novel environments in the face of recurrent and maladaptive gene flow. Here, we use a combination of analytical approximations and individual‐based simulations to explore how dominance influences the likelihood of adaptation to novel peripheral environments. We demonstrate that in the face of recurrent maladaptive gene flow, recessive alleles can fuel adaptation only when their frequency exceeds a critical threshold within the ancestral range.     

The effect of geographic range and dichogamy on genetic variability and population genetic structure in Tricyrtis section Flavae

Populations of each of the four species of Tricyrtis sect. Flavae were sampled using enzyme electrophoresis to examine the effect of geographic range and dichogamy on the genetic diversity of the species. The most widespread species, T. nana, had the lowest level of genetic diversity at both the population and the species level. The depauperate genetic diversity at the population level of T. nana appears to result from the high self-fertilization of the species. The low genetic diversity at the species level of T. nana probably resulted from the bottleneck effect during the speciation process in which this species diverged from the progenitor species, T. flava. Genetic differentiation among populations was high in both adichogamous T. nana and protandrous T. flava. High self-fertilization in T. nana and the colonizing nature of T. flava are likely the main factors causing the differentiated population genetic structure. In contrast to a previous study on chloroplast DNA (cpDNA) variation in Tricyrtis sect. Flavae, T. nana was most closely related to T. flava, which corresponds to the morphological resemblance of both species.     

Strong genetic differentiation but not local adaptation toward the range limit of a coastal dune plant

All species have limited geographic distributions; but the ecological and evolutionary mechanisms causing range limits are largely unknown. That many species’ geographic range limits are coincident with niche limits suggests limited evolutionary potential of marginal populations to adapt to conditions experienced beyond the range. We provide a test of range limit theory by combining population genetic analysis of microsatellite polymorphisms with a transplant experiment within, at the edge of, and 60 km beyond the northern range of a coastal dune plant. Contrary to expectations, lifetime fitness increased toward the range limit with highest fitness achieved by most populations at and beyond the range edge. Genetic differentiation among populations was strong, with very low, nondirectional gene flow suggesting range limitation via constraints to dispersal. In contrast, however, local adaptation was negligible, and a distance‐dependent decline in fitness only occurred for those populations furthest from home when planted beyond the range limit. These results challenge a commonly held assumption that stable range limits match niche limits, but also raise questions about the unique value of peripheral populations in expanding species’ geographical ranges.     

Population genetic analysis of a recent range expansion: mechanisms regulating the poleward range limit in the volcano barnacle Tetraclita rubescens

As range shifts coincident with climate change have become increasingly well documented, efforts to describe the causes of range boundaries have increased. Three mechanisms—genetic impoverishment, migration load, or a physical barrier to dispersal—are well described theoretically, but the data needed to distinguish among them have rarely been collected. We describe the distribution, abundance, genetic variation, and environment of Tetraclita rubescens, an intertidal barnacle that expanded its northern range limit by several hundreds of kilometres from San Francisco, CA, USA, since the 1970s. We compare geographic variation in abundance with abiotic and biotic patterns, including sea surface temperatures and the distributions of 387 co‐occurring species, and describe genetic variation in cytochrome c oxidase subunit I, mitochondrial noncoding region, and nine microsatellite loci from 27 locations between Bahia Magdalena (California Baja Sur, Mexico) and Cape Mendocino (CA, USA). We find very high gene flow, high genetic diversity, and a gradient in physical environmental variation coincident with the range limit. We infer that the primary cause of the northern range boundary in T. rubescens is migration load arising from flow of maladapted alleles into peripheral locations and that environmental change, which could have reduced selection against genotypes immigrating into the newly colonized portion of the range, is the most likely cause of the observed range expansion. Because environmental change could similarly affect all taxa in a region whose distributional limits are established by migration load, these mechanisms may be common causes of range boundaries and largely synchronous multi‐species range expansions.     

Range margin populations show high climate adaptation lags in European trees

How populations of long‐living species respond to climate change depends on phenotypic plasticity and local adaptation processes. Marginal populations are expected to have lags in adaptation (i.e. differences between the climatic optimum that maximizes population fitness and the local climate) because they receive pre‐adapted alleles from core populations preventing them from reaching a local optimum in their climatically marginal habitat. Yet, whether adaptation lags in marginal populations are a common feature across phylogenetically and ecologically different species and how lags can change with climate change remain unexplored. To test for range‐wide patterns of phenotypic variation and adaptation lags of populations to climate, we (a) built model ensembles of tree height accounting for the climate of population origin and the climate of the site for 706 populations monitored in 97 common garden experiments covering the range of six European forest tree species; (b) estimated populations' adaptation lags as the differences between the climatic optimum that maximizes tree height and the climate of the origin of each population; (c) identified adaptation lag patterns for populations coming from the warm/dry and cold/wet margins and from the distribution core of each species range. We found that (a) phenotypic variation is driven by either temperature or precipitation; (b) adaptation lags are consistently higher in climatic margin populations (cold/warm, dry/wet) than in core populations; (c) predictions for future warmer climates suggest adaptation lags would decrease in cold margin populations, slightly increasing tree height, while adaptation lags would increase in core and warm margin populations, sharply decreasing tree height. Our results suggest that warm margin populations are the most vulnerable to climate change, but understanding how these populations can cope with future climates depend on whether other fitness‐related traits could show similar adaptation lag patterns.     

High genetic diversity declines towards the geographic range periphery of Adonis vernalis,a Eurasian dry grassland plant

Genetic diversity is important for species' fitness and evolutionary processes but our knowledge on how it varies across a species' distribution range is limited. The abundant centre hypothesis (ACH) predicts that populations become smaller and more isolated towards the geographic range periphery – a pattern that in turn should be associated with decreasing genetic diversity and increasing genetic differentiation. We tested this hypothesis in Adonis vernalis, a dry grassland plant with an extensive Eurasian distribution. Its life‐history traits and distribution characteristics suggest a low genetic diversity that decreases and a high genetic differentiation that increases towards the range edge. We analysed AFLP fingerprints in 28 populations along a 4698‐km transect from the geographic range core in Russia to the western range periphery in Central and Western Europe. Contrary to our expectation, our analysis revealed high genetic diversity (range of proportion of polymorphic bands = 56–81%, He = 0.168–0.238) and low genetic differentiation across populations (Φ_ST = 0.18). However, in congruence with the genetic predictions of the ACH, genetic diversity decreased and genetic differentiation increased towards the range periphery. Spanish populations were genetically distinct, suggesting a divergent post‐glacial history in this region. The high genetic diversity and low genetic differentiation in the remaining A. vernalis populations is surprising given the species' life‐history traits and points to the possibility that the species has been widely distributed in the studied region or that it has migrated from a diverse source in an East–West direction, in the past.     

Local adaptation at the range peripheries of Sitka spruce

High‐dispersal rates in heterogeneous environments and historical rapid range expansion can hamper local adaptation; however, we often see clinal variation in high‐dispersal tree species. To understand the mechanisms of the species’ distribution, we investigated local adaptation and adaptive plasticity in a range‐wide context in Sitka spruce, a wind‐pollinated tree species that has recently expanded its range after glaciations. Phenotypic traits were observed using growth chamber experiments that mimicked temperature and photoperiodic regimes from the limits of the species realized niche. Bud phenology exhibited parallel reaction norms among populations; however, putatively adaptive plasticity and strong divergent selection were seen in bud burst and bud set timing respectively. Natural selection appears to have favoured genotypes that maximize growth rate during available frost‐free periods in each environment. We conclude that Sitka spruce has developed local adaptation and adaptive plasticity throughout its range in response to current climatic conditions despite generally high pollen flow and recent range expansion.     

Evolutionary shifts to self‐fertilisation restricted to geographic range margins in North American Arabidopsis lyrata

Cross‐fertilisation predominates in eukaryotes, but shifts to self‐fertilisation are common and ecologically and evolutionarily important. Reproductive assurance under outcross gamete limitation is one eco‐evolutionary process held responsible for the shift to selfing. Although small effective population size is a situation where selfing plants could theoretically benefit from reproductive assurance, empirical tests of the role of population size are rare. Here, we show that selfing evolved repeatedly at range margins, where historical demographic processes produced low effective population sizes. Outcrossing populations of North American Arabidopsis lyrata have low genetic diversity at geographic margins, with a signature of post‐glacial range expansion in the north and rear‐edge isolation in the south. Selfing populations occur at the margins of two genetic groups and never in their interior. These results corroborate small effective population size as the promoter of self‐fertilisation and have important implications for our understanding of species turnover, range limits and range dynamics.     

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.     

Experimental studies of adaptation in Clarkia xantiana. II. Fitness variation across a subspecies border

Because the range boundary is the locale beyond which a taxon fails to persist, it provides a unique opportunity for studying the limits on adaptive evolution. Adaptive constraints on range expansion are perplexing in view of widespread ecotypic differentiation by habitat and region within a species' range (regional adaptation) and rapid evolutionary response to novel environments. In this study of two parapatric subspecies, Clarkia xantiana ssp. xantiana and C. x. ssp. parviflora, we compared the fitness of population transplants within their native region, in a non-native region within the native range, and in the non-native range to assess whether range expansion might be limited by a greater intensity of selection on colonists of a new range versus a new region within the range. The combined range of the two subspecies spans a west-to-east gradient of declining precipitation in the Sierra Nevada of California, with ssp. xantiana in the west being replaced by ssp. parviflora in the east. Both subspecies had significantly higher fitness in the native range (range adaptation), whereas regional adaptation was weak and was found only in the predominantly outcrossing ssp. xantiana but was absent in the inbreeding ssp. parvifilora. Because selection intensity on transplants was much stronger in the non-native range relative to non-native regions, there is a larger adaptive barrier to range versus regional expansion. Three of five sequential fitness components accounted for regional and range adaptation, but only one of them, survivorship from germination to flowering, contributed to both. Flower number contributed to regional adaptation in ssp. xantiana and fruit set (number of fruits per flower) to range adaptation. Differential survivorship of the two taxa or regional populations of ssp. xantiana in non-native environments was attributable, in part, to biotic interactions, including competition, herbivory, and pollination. For example, low fruit set in ssp. xantiana in the east was likely due to the absence of its principal specialist bee pollinators in ssp. parviflora's range. Thus, convergence on self-fertilization may be necessary for ssp. xantiana to invade ssp. parviflora's range, but the evolution of outcrossing would not be required for ssp. parviflora to invade ssp. xantiana's range.     

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θ_sil=∼0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima''s D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (π_tot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4N_e generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role.     

Conservation implications of genetic variation between spatially and temporally distinct colonies of the endangered damselfly Coenagrion mercuriale

Abstract.  1. Good conservation management is underpinned by a thorough understanding of species' historical and contemporary dispersal capabilities along with the possible adaptive or neutral processes behind any spatio-temporal genetic structuring. These issues are investigated with respect to the rare damselfly Coenagrion mercuriale (Charpentier) – the only odonate species currently listed in the U.K.'s Biodiversity Action Plan – in east Devon where its distribution has become fragmented.
2. The two east Devon C. mercuriale populations, only 3.5 km apart, have accumulated strong differences in frequencies of alleles at 14 microsatellite loci as a consequence of poor adult dispersal and drift. There is no contemporary migration between sites.
3. A genetic signature of population decline at both sites corresponds with known demographic reductions. Coenagrion mercuriale in east Devon are now significantly genetically less diverse than those from a population stronghold in the Itchen Valley.
4. Colonies would benefit from improved connectivity between areas and possibly by a transfer of individuals from other ecologically similar areas.
5. Because C. mercuriale has a semivoltine life cycle throughout the U.K., the possibility that alternate-year cohorts are reproductively isolated is explored. Genetic differentiation among cohorts is an order of magnitude less than between sites, suggesting that some larvae delay their development into adults for a year and recruit to a different cohort.
6. To our knowledge, this is the first study to document migration and gene flow between alternate-year cohorts in a species of odonate. From a conservation standpoint, the cohorts do not require separate management.     

The performance of range maps and species distribution models representing the geographic variation of species richness at different resolutions

Aim The method used to generate hypotheses about species distributions, in addition to spatial scale, may affect the biodiversity patterns that are then observed. We compared the performance of range maps and MaxEnt species distribution models at different spatial resolutions by examining the degree of similarity between predicted species richness and composition against observed values from well‐surveyed cells (WSCs). Location Mexico. Methods We estimated amphibian richness distributions at five spatial resolutions (from 0.083° to 2°) by overlaying 370 individual range maps or MaxEnt predictions, comparing the similarity of the spatial patterns and correlating predicted values with the observed values for WSCs. Additionally, we looked at species composition and assessed commission and omission errors associated with each method. Results MaxEnt predictions reveal greater geographic differences in richness between species rich and species poor regions than the range maps did at the five resolutions assessed. Correlations between species richness values estimated by either of the two procedures and the observed values from the WSCs increased with decreasing resolution. The slopes of the regressions between the predicted and observed values indicate that MaxEnt overpredicts observed species richness at all of the resolutions used, while range maps underpredict them, except at the finest resolution. Prediction errors did not vary significantly between methods at any resolution and tended to decrease with decreasing resolution. The accuracy of both procedures was clearly different when commission and omission errors were examined separately. Main conclusions Despite the congruent increase in the geographic richness patterns obtained from both procedures as resolution decreases, the maps created with these methods cannot be used interchangeably because of notable differences in the species compositions they report.     

Lack of concordance between genetic diversity estimates at the molecular and quantitative-trait levels

In many applications of population genetics, particularly in the field of conservation biology, estimates of molecular diversity are used as surrogate indicators of less easily acquired measures of genetic variation for quantitative traits. The general validity of this approach to inferring levels of quantitative genetic variation within populations is called into question by the demonstration that estimates of molecular and quantitative-genetic variation are essentially uncorrelated in natural populations of Daphnia, one of the few organisms for which multiple estimates of both quantities are available. On the other hand, molecular measures of population subdivision seem to give conservatively low estimates of the degree of genetic subdivision at the level of quantitative traits. This suggests that although molecular markers provide little information on the level of genetic variation for quantitative traits within populations, they may be valid indicators of population subdivision for such characters.