Evolutionary insights from studies of geographic variation: Contemporary variation and looking to the future

Future carbon sequestration potential in a widespread transcontinental boreal tree species: Standing genetic variation matters!

Climate change (CC) necessitates reforestation/afforestation programs to mitigate its impacts and maximize carbon sequestration. But comprehending how tree growth, a proxy for fitness and resilience, responds to CC is critical to maximize these programs' effectiveness. Variability in tree response to CC across populations can notably be influenced by the standing genetic variation encompassing both neutral and adaptive genetic diversity. Here, a framework is proposed to assess tree growth potential at the population scale while accounting for standing genetic variation. We applied this framework to black spruce (BS, Picea mariana [Mill] B.S.P.), with the objectives to (1) determine the key climate variables having impacted BS growth response from 1974 to 2019, (2) examine the relative roles of local adaptation and the phylogeographic structure in this response, and (3) project BS growth under two Shared Socioeconomic Pathways while taking standing genetic variation into account. We modeled growth using a machine learning algorithm trained with dendroecological and genetic data obtained from over 2600 trees (62 populations divided in three genetic clusters) in four 48-year-old common gardens, and simulated growth until year 2100 at the common garden locations. Our study revealed that high summer and autumn temperatures negatively impacted BS growth. As a consequence of warming, this species is projected to experience a decline in growth by the end of the century, suggesting maladaptation to anticipated CC and a potential threat to its carbon sequestration capacity. This being said, we observed a clear difference in response to CC within and among genetic clusters, with the western cluster being more impacted than the central and eastern clusters. Our results show that intraspecific genetic variation, notably associated with the phylogeographic structure, must be considered when estimating the response of widespread species to CC.

Forest tree species adaptation to climate across biomes: Building on the legacy of ecological genetics to anticipate responses to climate change

Intraspecific variation plays a critical role in extant and future forest responses to climate change. Forest tree species with wide climatic niches rely on the intraspecific variation resulting from genetic adaptation and phenotypic plasticity to accommodate spatial and temporal climate variability. A centuries-old legacy of forest ecological genetics and provenance trials has provided a strong foundation upon which to continue building on this knowledge, which is critical to maintain climate-adapted forests. Our overall objective is to understand forest trees intraspecific responses to climate across species and biomes, while our specific objectives are to describe ecological genetics models used to build our foundational knowledge, summarize modeling approaches that have expanded the traditional toolset, and extensively review the literature from 1994 to 2021 to highlight the main contributions of this legacy and the new analyzes of provenance trials. We reviewed 103 studies comprising at least three common gardens, which covered 58 forest tree species, 28 of them with range-wide studies. Although studies using provenance trial data cover mostly commercially important forest tree species from temperate and boreal biomes, this synthesis provides a global overview of forest tree species adaptation to climate. We found that evidence for genetic adaptation to local climate is commonly present in the species studied (79%), being more common in conifers (87.5%) than in broadleaf species (67%). In 57% of the species, clines in fitness-related traits were associated with temperature variables, in 14% of the species with precipitation, and in 25% of the species with both. Evidence of adaptation lags was found in 50% of the species with range-wide studies. We conclude that ecological genetics models and analysis of provenance trial data provide excellent insights on intraspecific genetic variation, whereas the role and limits of phenotypic plasticity, which will likely determine the fate of extant forests, is vastly understudied.

Geographic variation in offspring size: Long- and short-term climate affect mean seed mass of Streptanthus populations

Offspring size is a key functional trait that can affect subsequent life history stages; in many species, it exhibits both local adaptation and phenotypic plasticity. Variation among populations in offspring size may be explained by various factors, including local climatic conditions. However, geographic variation in climate may be partitioned into long‐term and interannual sources of variation, which may differ in their effects on population mean offspring size. To assess environmental correlates of offspring size, we evaluated geographic variation in seed mass among 88 populations representing 6 species of Streptanthus (Brassicaceae) distributed across a broad climatic gradient in California. We examined the effects of temperature‐mediated growing season length and precipitation on population mean seed mass to determine whether it is best explained by (1) long‐term mean climatic conditions; (2) interannual climate anomalies (i.e., deviations in climate from long‐term means) during the year of seed development, or (3) interactions between climate variables. Both long‐term mean climate and climate anomalies in the year of collection were associated with population mean seed mass, but their effects differed in direction and magnitude. Relatively large seeds were produced at chronically wet sites but also during drier‐than‐average years. This contrast indicates that these associations may be generated by different mechanisms (i.e., adaptive evolution vs. phenotypic plasticity) and may be evidence of countergradient plasticity in seed mass. In addition, populations occurring in locations characterized by relatively long growing seasons produced comparatively large seeds, particularly among chronically dry sites. This study highlights the need to consider that the responses of seed mass to long‐term versus recent climatic conditions may differ and that climate variables may interact to predict seed mass. Such considerations are especially important when using these patterns to forecast the long‐ and short‐term responses of seed mass to climate change. The results presented here also contribute to our broader understanding of how climate drives long‐term (e.g., local adaptation) and short‐term (e.g., phenotypic plasticity) variation in functional traits, such as offspring size across landscapes.

Ecological genetics of Juglans nigra: Differences in early growth patterns of natural populations

Many boreal and temperate forest tree species distributed across large geographic ranges are composed of populations adapted to the climate they inhabit. Forestry provenance studies and common gardens provide evidence of local adaptation to climate when associations between fitness traits and the populations' home climates are observed. Most studies that evaluate tree height as a fitness trait do so at a specific point in time. In this study, we elucidate differences in early growth patterns in black walnut (Juglans nigra L.) populations by modeling height growth from seed up to age 11. The data comprise tree height measurements between ages 2 and 11 for 52 natural populations of black walnut collected through its geographic range and planted in one or more of 3 common gardens. We use the Chapman-Richards growth model in a mixed effects framework and test whether populations differ in growth patterns by incorporating populations' home climate into the model. In addition, we evaluate differences in populations' absolute growth and relative growth based on the fitted model. Models indicated that populations from warmer climates had the highest cumulative growth through time, with differences in average tree height between populations from home climates with a mean annual temperature (MAT) of 13°C and of 7°C estimated to be as high as 80% at age 3. Populations from warmer climates were also estimated to have higher and earlier maximum absolute growth rate than populations from colder climates. In addition, populations from warm climates were predicted to have higher relative growth rates at any given tree size. Results indicate that natural selection may shape early growth patterns of populations within a tree species, suggesting that fast early growth rates are likely selected for in relatively mild environments where competition rather than tolerance to environmental stressors becomes the dominant selection pressure.

Population climatic history predicts phenotypic responses in novel environments for Arabidopsis thaliana in North America

PREMISE

Determining how species perform in novel climatic environments is essential for understanding (1) responses to climate change and (2) evolutionary consequences of biological invasions. For the vast majority of species, the number of population characteristics that will predict performance and patterns of natural selection in novel locations in the wild remains limited.

METHODS

We evaluated phenological, vegetative, architectural, and fitness-related traits in experimental gardens in contrasting climates (Ontario, Canada, and South Carolina, USA) in the North American non-native distribution of Arabidopsis thaliana. We assessed the effects of climatic distance, geographic distance, and genetic features of history on performance and patterns of natural selection in the novel garden settings.

RESULTS

We found that plants had greater survivorship, flowered earlier, were larger, and produced more fruit in the south, and that genotype-by-environment interactions were significant between gardens. However, our analyses revealed similar patterns of natural selection between gardens in distinct climate zones. After accounting for genetic ancestry, we also detected that population climatic distance best predicted performance within gardens.

CONCLUSIONS

These data suggest that colonization success in novel, non-native environments is determined by a combination of climate and genetic history. When performance at novel sites was assessed with seed sources from geographically and genetically disparate, established non-native populations, proximity to the garden alone was insufficient to predict performance. Our study highlights the need to evaluate seed sources from diverse origins to describe comprehensively phenotypic responses to novel environments, particularly for taxa in which many source populations may contribute to colonization.

Genetic data reveals a complex history of multiple admixture events in presently allopatric wild gingers (Asarum spp.) showing intertaxonomic clinal variation in calyx lobe length

Clinal variation is a major pattern of observed phenotypic diversity and identifying underlying demographic processes is a necessary step to understand the establishment of clinal variation. The wild ginger series Sakawanum (genus Asarum) comprises four taxa, which exhibit intertaxonomic clinal variation in calyx lobe length across two continental islands isolated by a sea strait. To test alternative hypotheses of the evolutionary history and to determine the implications for the formation of clinal variation, we conducted approximate Bayesian computation (ABC) analysis and ecological niche modeling (ENM). ABC analysis indicated that the scenario assuming multiple admixture events was strongly supported. This scenario assumed two admixture events occurred between morphologically distinct taxa, likely leading to the generation of intermediate taxa. One of the admixture events was estimated to have occurred during the last glacial maximum (LGM), during which the taxa were estimated to have formed a common refugia in southern areas by ENM analysis. Although four taxa are currently distributed allopatrically on different islands and trans-oceanic dispersal appears unlikely, the formation of a land bridge and the geographic range shift to refugia would have allowed secondary contact between previously isolated taxa. This study suggests that clinal variation can be shaped by demographic history including multiple admixtures due to climatic oscillations.

Capturing genetic variation in crop wild relatives: An evolutionary approach

Crop wild relatives (CWRs) offer novel genetic resources for crop improvement. To assist in the urgent need to collect and conserve CWR germplasm, we advance here the concept of an "evolutionary" approach. Central to this approach is the predictive use of spatial proxies of evolutionary processes (natural selection, gene flow and genetic drift) to locate and capture genetic variation. As a means to help validate this concept, we screened wild-collected genotypes of woodland strawberry (Fragaria vesca) in a common garden. A quantitative genetic approach was then used to test the ability of two such proxies-mesoclimatic variation (a proxy of natural selection) and landscape isolation and geographic distance between populations (proxies of gene flow potential)-to predict spatial genetic variation in three quantitative traits (plant size, early season flower number and flower frost tolerance). Our results indicated a significant but variable effect of mesoclimatic conditions in structuring genetic variation in the wild, in addition to other undetermined regional scale processes. As a proxy of gene flow potential, landscape isolation was also a likely determinant of observed patterns-as opposed to, and regardless of, geographic distance between populations. We conclude that harnessing proxies of adaptive and nonadaptive evolutionary processes could provide a robust and valuable means to identify genetic variation in CWRs. We thus advocate wider use and development of this approach amongst researchers, breeders and practitioners, to expedite the capture and in situ conservation of genetic resources provided by crop wild relatives.

Spatio-temporal variation in fitness responses to contrasting environments in Arabidopsis thaliana

The evolutionary response of organisms to global climate change is expected to be strongly conditioned by preexisting standing genetic variation. In addition, natural selection imposed by global climate change on fitness-related traits can be heterogeneous over time. We estimated selection of life-history traits of an entire genetic lineage of the plant Arabidopsis thaliana occurring in north-western Iberian Peninsula that were transplanted over multiple years into two environmentally contrasting field sites in southern Spain, as southern environments are expected to move progressively northwards with climate change in the Iberian Peninsula. The results indicated that natural selection on flowering time prevailed over that on recruitment. Selection favored early flowering in six of eight experiments and late flowering in the other two. Such heterogeneity of selection for flowering time might be a powerful mechanism for maintaining genetic diversity in the long run. We also found that north-western A. thaliana accessions from warmer environments exhibited higher fitness and higher phenotypic plasticity for flowering time in southern experimental facilities. Overall, our transplant experiments suggested that north-western Iberian A. thaliana has the means to cope with increasingly warmer environments in the region as predicted by trends in global climate change models.

Polymorphic nuclear markers for coastal plant species with dynamic geographic distributions, the rock samphire (Crithmum maritimum) and the vulnerable dune pansy (Viola tricolor subsp. curtisii)

Identifying spatial patterns of genetic differentiation across a species range is critical to set up conservation and restoration decision-making. This is especially timely, since global change triggers shifts in species' geographic distribution and in the geographical variation of mating system and patterns of genetic differentiation, with varying consequences at the trailing and leading edges of a species' distribution. Using 454 pyrosequencing, we developed nuclear microsatellite loci for two plant species showing a strictly coastal geographical distribution and contrasting range dynamics: the expanding rock samphire (Crithmum maritimum, 21 loci) and the highly endangered and receding dune pansy (Viola tricolor subsp. curtisii, 12 loci). Population genetic structure was then assessed by genotyping more than 100 individuals from four populations of each of the two target species. Rock samphire displayed high levels of genetic differentiation (F_ST = 0.38), and a genetic structure typical of a mostly selfing species (F_IS ranging from 0.16 to 0.58). Populations of dune pansy showed a less pronounced level of population structuring (F_ST = 0.25) and a genotypic structure more suggestive of a mixed-mating system when excluding two loci with heterozygote excess. These results demonstrate that the genetic markers developed here are useful to assess the mating system of populations of these two species. They will be tools of choice to investigate phylogeographical patterns and variation in mating system over the geographical distribution ranges for two coastal plant species that are subject to dynamic evolution due to rapid contemporary global change.

Relative contributions of neutral and non-neutral processes to clinal variation in calyx lobe length in the series Sakawanum (Asarum: Aristolochiaceae)

Background and Aims

Clines, the gradual variation in measurable traits along a geographical axis, play a major role in evolution and can contribute to our understanding of the relative roles of selective and neutral process in trait variation. Using genetic and morphological analyses, the relative contributions of neutral and non-neutral processes were explored to infer the evolutionary history of species of the series Sakawanum (genus Asarum), which shows significant clinal variation in calyx lobe length.

Methods

A total of 27 populations covering the natural geographical distribution of the series Sakawanum were sampled. Six nuclear microsatellite markers were used to investigate genetic structure and genetic diversity. The lengths of calyx lobes of multiple populations were measured to quantify their geographical and taxonomic differentiation. To detect the potential impact of selective pressure, morphological differentiation was compared with genetic differentiation (QCT-FST comparison).

Key Results

Average calyx lobe length of A. minamitanianum was 124.11 mm, while that of A. costatum was 13.80 mm. Though gradually changing along the geographical axis within series, calyx lobe lengths were significantly differentiated among the taxa. Genetic differentiation between taxa was low (FST = 0.099), but a significant geographical structure along the morphological cline was detected. Except for one taxon pair, pairwise QCT values were significantly higher than the neutral genetic measures of FST and G'ST.

Conclusions

Divergent selection may have driven the calyx lobe length variation in series Sakawanum taxa, although the underlying mechanism is still not clear. The low genetic differentiation indicates recent divergence and/or gene flows between geographically close taxa. These neutral processes would also affect the clinal variation in calyx lobe lengths. Overall, this study implies the roles of population history and divergent selection in shaping the current cline of a flower trait in the series Sakawanum.

Biotic and abiotic drivers of intraspecific trait variation within plant populations of three herbaceous plant species along a latitudinal gradient

Background

The importance of intraspecific trait variation (ITV) is increasingly acknowledged among plant ecologists. However, our understanding of what drives ITV between individual plants (ITV_BI) at the population level is still limited. Contrasting theoretical hypotheses state that ITV_BI can be either suppressed (stress-reduced plasticity hypothesis) or enhanced (stress-induced variability hypothesis) under high abiotic stress. Similarly, other hypotheses predict either suppressed (niche packing hypothesis) or enhanced ITV_BI (individual variation hypothesis) under high niche packing in species rich communities. In this study we assess the relative effects of both abiotic and biotic niche effects on ITV_BI of four functional traits (leaf area, specific leaf area, plant height and seed mass), for three herbaceous plant species across a 2300 km long gradient in Europe. The study species were the slow colonizing Anemone nemorosa, a species with intermediate colonization rates, Milium effusum, and the fast colonizing, non-native Impatiens glandulifera.

Results

Climatic stress consistently increased ITV_BI across species and traits. Soil nutrient stress, on the other hand, reduced ITV_BI for A. nemorosa and I. glandulifera, but had a reversed effect for M. effusum. We furthermore observed a reversed effect of high niche packing on ITV_BI for the fast colonizing non-native I. glandulifera (increased ITV_BI), as compared to the slow colonizing native A. nemorosa and M. effusum (reduced ITV_BI). Additionally, ITV_BI in the fast colonizing species tended to be highest for the vegetative traits plant height and leaf area, but lowest for the measured generative trait seed mass.

Conclusions

This study shows that stress can both reduce and increase ITV_BI, seemingly supporting both the stress-reduced plasticity and stress-induced variability hypotheses. Similarly, niche packing effects on ITV_BI supported both the niche packing hypothesis and the individual variation hypothesis. These results clearly illustrates the importance of simultaneously evaluating both abiotic and biotic factors on ITV_BI. This study adds to the growing realization that within-population trait variation should not be ignored and can provide valuable ecological insights.

Electronic supplementary material

The online version of this article (10.1186/s12898-017-0151-y) contains supplementary material, which is available to authorized users.

Landscape genomic insights into the historic migration of mountain hemlock in response to Holocene climate change

PREMISE OF THE STUDY

Untangling alternative historic dispersal pathways in long-lived tree species is critical to better understand how temperate tree species may respond to climatic change. However, disentangling these alternative pathways is often difficult. Emerging genomic technologies and landscape genetics techniques improve our ability to assess these pathways in natural systems. We address the question to what degree have microrefugial patches and long-distance dispersal been responsible for the colonization of mountain hemlock (Tsuga mertensiana) on the Alaskan Kenai Peninsula.

METHODS

We used double-digest restriction-associated DNA sequencing (ddRADseq) to identify genetic variants across eight mountain hemlock sample sites on the Kenai Peninsula, Alaska. We assessed genetic diversity and linkage disequilibrium using landscape and population genetics approaches. Alternative historic dispersal pathways were assessed using discriminant analysis of principle components and electrical circuit theory.

KEY RESULTS

A combination of decreasing diversity, high gene flow, and landscape connectivity indicates that mountain hemlock colonization on the Kenai Peninsula is the result of long-distance dispersal. We found that contemporary climate best explained gene flow patterns and that isolation by resistance was a better model explaining genetic variation than isolation by distance.

CONCLUSIONS

Our findings support the conclusion that mountain hemlock colonization is the result of several long-distance dispersal events following Pleistocene glaciation. The high dispersal capability suggests that mountain hemlock may be able to respond to future climate change and expand its range as new habitat opens along its northern distribution.