Adaptation of lodgepole pine and interior spruce to climate: implications for reforestation in a warming world

Weak local adaptation to drought in seedlings of a widespread conifer

Tree seedlings from populations native to drier regions are often assumed to be more drought tolerant than those from wetter provenances. However, intraspecific variation in drought tolerance has not been well-characterized despite being critical for developing climate change mitigation and adaptation strategies, and for predicting the effects of drought on forests. We used a large-scale common garden drought-to-death experiment to assess range-wide variation in drought tolerance, measured by decline of photosynthetic efficiency, growth, and plastic responses to extreme summer drought in seedlings of 73 natural populations of the two main varieties of Douglas-fir (Pseudotsuga menziesii var. menziesii and var. glauca). Local adaptation to drought was weak in var. glauca and nearly absent in menziesii. Var. glauca showed higher tolerance to drought but slower growth than var. menziesii. Clinal variation in drought tolerance and growth species-wide was mainly associated with temperature rather than precipitation. A higher degree of plasticity for growth was observed in var. menziesii in response to extreme drought. Genetic variation for drought tolerance in seedlings within varieties is maintained primarily within populations. Selective breeding within populations may facilitate adaptation to drought more than assisted gene flow.

Local adaptation: Causal agents of selection and adaptive trait divergence

Divergent selection across the landscape can favor the evolution of local adaptation in populations experiencing contrasting conditions. Local adaptation is widely observed in a diversity of taxa, yet we have a surprisingly limited understanding of the mechanisms that give rise to it. For instance, few have experimentally confirmed the biotic and abiotic variables that promote local adaptation, and fewer yet have identified the phenotypic targets of selection that mediate local adaptation. Here, we highlight critical gaps in our understanding of the process of local adaptation and discuss insights emerging from in-depth investigations of the agents of selection that drive local adaptation, the phenotypes they target, and the genetic basis of these phenotypes. We review historical and contemporary methods for assessing local adaptation, explore whether local adaptation manifests differently across life history, and evaluate constraints on local adaptation.

Using landscape genomics to delineate seed and breeding zones for lodgepole pine

Seed and breeding zones traditionally are delineated based on local adaptation of phenotypic traits associated with climate variables, an approach requiring long-term field experiments. In this study, we applied a landscape genomics approach to delineate seed and breeding zones for lodgepole pine. We used a gradient forest (GF) model to select environment-associated single nucleotide polymorphisms (SNPs) using three SNP datasets (full, neutral and candidate) and 20 climate variables for 1906 lodgepole pine (Pinus contorta) individuals in British Columbia and Alberta, Canada. The two GF models built with the full (28 954) and candidate (982) SNPs were compared. The GF models identified winter-related climate as major climatic factors driving genomic patterns of lodgepole pine's local adaptation. Based on the genomic gradients predicted by the full and candidate GF models, lodgepole pine distribution range in British Columbia and Alberta was delineated into six seed and breeding zones. Our approach is a novel and effective alternative to traditional common garden approaches for delineating seed and breeding zone, and could be applied to tree species lacking data from provenance trials or common garden experiments.

Leaf Economic and Hydraulic Traits Signal Disparate Climate Adaptation Patterns in Two Co-Occurring Woodland Eucalypts

With climate change impacting trees worldwide, enhancing adaptation capacity has become an important goal of provenance translocation strategies for forestry, ecological renovation, and biodiversity conservation. Given that not every species can be studied in detail, it is important to understand the extent to which climate adaptation patterns can be generalised across species, in terms of the selective agents and traits involved. We here compare patterns of genetic-based population (co)variation in leaf economic and hydraulic traits, climate–trait associations, and genomic differentiation of two widespread tree species (Eucalyptus pauciflora and E. ovata). We studied 2-year-old trees growing in a common-garden trial established with progeny from populations of both species, pair-sampled from 22 localities across their overlapping native distribution in Tasmania, Australia. Despite originating from the same climatic gradients, the species differed in their levels of population variance and trait covariance, patterns of population variation within each species were uncorrelated, and the species had different climate–trait associations. Further, the pattern of genomic differentiation among populations was uncorrelated between species, and population differentiation in leaf traits was mostly uncorrelated with genomic differentiation. We discuss hypotheses to explain this decoupling of patterns and propose that the choice of seed provenances for climate-based plantings needs to account for multiple dimensions of climate change unless species-specific information is available.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of climate change.

Non-Native Forest Tree Species in Europe: The Question of Seed Origin in Afforestation

Non-native forest tree species have been introduced in Europe since the 16th century, but only in the second half of the 20th century the significance of the seed source origin for their economic use was recognized, resulting in the establishment of numerous provenance trials at a national, regional, European and International level, as those led by IUFRO. Breeding programs have also been launched in the continent for the most economically important species. Aim of this work is the formulation of provenance recommendations for planting of five non-native tree species in Europe (Douglas fir, grand fir, Sitka spruce, lodgepole pine and black locust), based on the information obtained from twenty countries, in the frame of the EU FP-1403 NNEXT Cost Action. The survey revealed that official and non-official national recommendations, based on provenance research results, have been elaborated and followed at a different level and extend for the above five species, but only for Douglas fir recommendations exist in almost all the participating to the survey countries. The compilation of provenance recommendations across Europe for each species is presented in the current work. Besides the recommended introduced seed sources, European seed sources are also preferred for planting, due to ease of access and high availability of forest reproductive material. European breeding programs yielding genetic material of high productivity and quality constitute currently the seed source of choice for several species and countries. Consolidation of trial data obtained across countries will allow the joint analysis that is urgently needed to draw solid conclusions, and will facilitate the development of ‘Universal-Response-Functions’ for the species of interest, rendering possible the identification of the genetic material suitable for global change. New provenance trial series that will test seed sources from the entire climatic range of the species, established in sites falling within and outside the environmental envelopes of their natural ranges, are urgently needed to pinpoint and understand the species-specific climate constraints, as well as to correlate functional traits to the seed origin and the environmental conditions of the test sites, so that the selection of suitable forest reproductive material of non-native tree species in the face of climate change can be feasible.

Optimizing sampling design and sequencing strategy for the genomic analysis of quantitative traits in natural populations

Mapping the genes underlying ecologically relevant traits in natural populations is fundamental to develop a molecular understanding of species adaptation. Current sequencing technologies enable the characterization of a species' genetic diversity across the landscape or even over its whole range. The relevant capture of the genetic diversity across the landscape is critical for a successful genetic mapping of traits and there are no clear guidelines on how to achieve an optimal sampling and which sequencing strategy to implement. Here we determine, through simulation, the sampling scheme that maximizes the power to map the genetic basis of a complex trait in an outbreeding species across an idealized landscape and draw genomic predictions for the trait, comparing individual and pool sequencing strategies. Our results show that quantitative trait locus detection power and prediction accuracy are higher when more populations over the landscape are sampled and this is more cost-effectively done with pool sequencing than with individual sequencing. Additionally, we recommend sampling populations from areas of high genetic diversity. As progress in sequencing enables the integration of trait-based functional ecology into landscape genomics studies, these findings will guide study designs allowing direct measures of genetic effects in natural populations across the environment.

Silica Particles Trigger the Exopolysaccharide Production of Harsh Environment Isolates of Growth-Promoting Rhizobacteria and Increase Their Ability to Enhance Wheat Biomass in Drought-Stressed Soils

In coming decades, drought is expected to expand globally owing to increased evaporation and reduced rainfall. Understanding, predicting, and controlling crop plants’ rhizosphere has the potential to manipulate its responses to environmental stress. Our plant growth-promoting rhizobacteria (PGPR) are isolated from a natural laboratory, ‘The Evolution Canyon’, Israel, (EC), from the wild progenitors of cereals, where they have been co-habituating with their hosts for long periods of time. The study revealed that commercial TM50 silica particles (SN) triggered the PGPR production of exopolysaccharides (EPS) containing D-glucuronate (D-GA). The increased EPS content increased the PGPR water-holding capacity (WHC) and osmotic pressure of the biofilm matrix, which led to enhanced plant biomass in drought-stressed growth environments. Light- and cryo-electron- microscopic studies showed that, in the presence of silica (SN) particles, bacterial morphology is changed, indicating that SNs are associated with significant reprogramming in bacteria. The findings encourage the development of large-scale methods for isolate formulation with natural silicas that ensure higher WHC and hyperosmolarity under field conditions. Osmotic pressure involvement of holobiont cohabitation is also discussed.

Annual aboveground carbon uptake enhancements from assisted gene flow in boreal black spruce forests are not long-lasting

Assisted gene flow between populations has been proposed as an adaptive forest management strategy that could contribute to the sequestration of carbon. Here we provide an assessment of the mitigation potential of assisted gene flow in 46 populations of the widespread boreal conifer Picea mariana, grown in two 42-year-old common garden experiments and established in contrasting Canadian boreal regions. We use a dendroecological approach taking into account phylogeographic structure to retrospectively analyse population phenotypic variability in annual aboveground net primary productivity (NPP). We compare population NPP phenotypes to detect signals of adaptive variation and/or the presence of phenotypic clines across tree lifespans, and assess genotype-by-environment interactions by evaluating climate and NPP relationships. Our results show a positive effect of assisted gene flow for a period of approximately 15 years following planting, after which there was little to no effect. Although not long lasting, well-informed assisted gene flow could accelerate the transition from carbon source to carbon sink after disturbance.

Adaptation to drought is coupled with slow growth, but independent from phenology in marginal silver fir (Abies alba Mill.) populations

Drought is one of the most important selection pressures for forest trees in the context of climate change. Yet, the different evolutionary mechanisms, and their environmental drivers, by which certain populations become more drought tolerant than others is still little understood. We studied adaptation to drought in 16 silver fir (Abies alba Mill.) populations from the French Mediterranean Alps by combining observations on seedlings from a greenhouse experiment (N = 8,199) and on adult tress in situ (N = 315). In the greenhouse, we followed half-sib families for four growing seasons for growth and phenology traits, and tested their water stress response in a "drought until death" experiment. Adult trees in the field were assessed for δ 13C, a proxy for water use efficiency, and genotyped at 357 SNP loci. SNP data was used to generate a null expectation for seedling trait divergence between populations in order to detect the signature of selection, and 31 environmental variables were used to identify the selective environment. We found that seedlings originating from populations with low soil water capacity grew more slowly, attained a smaller stature, and resisted water stress for a longer period of time in the greenhouse. Additionally, adult trees of these populations exhibited a higher water use efficiency as evidenced by their δ 13C. These results suggest a correlated evolution of the growth-drought tolerance trait complex. Population divergence in bud break phenology was adaptive only in the second growing season, and evolved independently from the growth-drought tolerance trait complex. Adaptive divergence in bud break phenology was principally driven by the inter- and intra-annual variation in temperature at the geographic origin of the population. Our results illustrate the different evolutionary strategies used by populations to cope with drought stress at the range limits across a highly heterogeneous landscape, and can be used to inform assisted migration programs.

Variation in the phenology of photosynthesis among eastern white pine provenances in response to warming

In higher-latitude trees, temperature and photoperiod control the beginning and end of the photosynthetically active season. Elevated temperature (ET) has advanced spring warming and delayed autumn cooling while photoperiod remains unchanged. We assessed the effects of warming on the length of the photosynthetically active season of three provenances of Pinus strobus L. seedlings from different latitudes, and evaluated the accuracy of the photochemical reflectance index (PRI) and the chlorophyll/carotenoid index (CCI) for tracking the predicted variation in spring and autumn phenology of photosynthesis among provenances. Seedlings from northern, local and southern P. strobus provenances were planted in a temperature-free-air-controlled enhancement (T-FACE) experiment and exposed to ET (+1.5/3°C; day/night). Over 18 months, we assessed photosynthetic phenology by measuring chlorophyll fluorescence, gas exchange, leaf spectral reflectance and pigment content. During autumn, all seedlings regardless of provenance followed the same sequence of phenological events with the initial downregulation of photosynthesis, followed by the modulation of non-photochemical quenching and associated adjustments of zeaxanthin pool sizes. However, the timing of autumn downregulation differed between provenances, with delayed onset in the southern provenance (SP) and earlier onset in the northern relative to the local provenance, indicating that photoperiod at the provenance origin is a dominant factor controlling autumn phenology. Experimental warming further delayed the downregulation of photosynthesis during autumn in the SP. A provenance effect during spring was also observed but was generally not significant. The vegetation indices PRI and CCI were both effective at tracking the seasonal variations of energy partitioning in needles and the differences of carotenoid pigments indicative of the stress status of needles. These results demonstrate that PRI and CCI can be useful tools for monitoring conifer phenology and for the remote monitoring of the length of the photosynthetically active season of conifers in a changing climate.

Local Adaptation to Biotic Interactions: A Meta-analysis across Latitudes

Adaptation to local conditions can increase species' geographic distributions and rates of diversification, but which components of the environment commonly drive local adaptation-particularly the importance of biotic interactions-is unclear. Biotic interactions should drive local adaptation when they impose consistent divergent selection; if this is common, we expect transplant experiments to detect more frequent and stronger local adaptation when biotic interactions are left intact. We tested this hypothesis using a meta-analysis of transplant experiments from >125 studies (mostly of plants). Overall, local adaptation was common, and biotic interactions affected fitness. Nevertheless, local adaptation was neither more common nor stronger when biotic interactions were left intact, either between experimental treatments within studies (control vs. biotic interactions experimentally manipulated) or between studies that used natural versus biotically altered transplant environments. However, the effect of ameliorating negative interactions varied with latitude, suggesting that interactions may promote local adaptation more often in tropical than in temperate ecosystems, although few tropical studies were available to test this. Our results suggest that biotic interactions often fail to drive local adaptation even though they strongly affect fitness, perhaps because temperate biotic environments are unpredictable at the spatiotemporal scales required for local adaptation.

Evaluating genomic data for management of local adaptation in a changing climate: A lodgepole pine case study

We evaluate genomic data, relative to phenotypic and climatic data, as a basis for assisted gene flow and genetic conservation. Using a seedling common garden trial of 281 lodgepole pine (Pinus contorta) populations from across western Canada, we compare genomic data to phenotypic and climatic data to assess their effectiveness in characterizing the climatic drivers and spatial scale of local adaptation in this species. We find that phenotype-associated loci are equivalent or slightly superior to climate data for describing local adaptation in seedling traits, but that climate data are superior to genomic data that have not been selected for phenotypic associations. We also find agreement between the climate variables associated with genomic variation and with 20-year heights from a long-term provenance trial, suggesting that genomic data may be a viable option for identifying climatic drivers of local adaptation where phenotypic data are unavailable. Genetic clines associated with the experimental traits occur at broad spatial scales, suggesting that standing variation of adaptive alleles for this and similar species does not require management at scales finer than those indicated by phenotypic data. This study demonstrates that genomic data are most useful when paired with phenotypic data, but can also fill some of the traditional roles of phenotypic data in management of species for which phenotypic trials are not feasible.

Northern forest tree populations are physiologically maladapted to drought

Northern forests at the leading edge of their distributions may not show increased primary productivity under climate warming, being limited by climatic extremes such as drought. Looking beyond tree growth to underlying physiological mechanisms is fundamental for accurate predictions of forest responses to climate warming and drought stress. Within a 32-year genetic field trial, we analyze relative contributions of xylem plasticity and inferred stomatal response to drought tolerance in regional populations of a widespread conifer. Genetic adaptation leads to varying responses under drought. Trailing-edge tree populations produce fewer tracheids with thicker cell walls, characteristic of drought-tolerance. Stomatal response explains the moderate drought tolerance of tree populations in central areas of the species range. Growth loss of the northern population is linked to low stomatal responsiveness combined with the production of tracheids with thinner cell walls. Forests of the western boreal may therefore lack physiological adaptations necessary to tolerate drier conditions.

Modularity of genes involved in local adaptation to climate despite physical linkage

BACKGROUND

Linkage among genes experiencing different selection pressures can make natural selection less efficient. Theory predicts that when local adaptation is driven by complex and non-covarying stresses, increased linkage is favored for alleles with similar pleiotropic effects, with increased recombination favored among alleles with contrasting pleiotropic effects. Here, we introduce a framework to test these predictions with a co-association network analysis, which clusters loci based on differing associations. We use this framework to study the genetic architecture of local adaptation to climate in lodgepole pine, Pinus contorta, based on associations with environments.

RESULTS

We identify many clusters of candidate genes and SNPs associated with distinct environments, including aspects of aridity and freezing, and discover low recombination rates among some candidate genes in different clusters. Only a few genes contain SNPs with effects on more than one distinct aspect of climate. There is limited correspondence between co-association networks and gene regulatory networks. We further show how associations with environmental principal components can lead to misinterpretation. Finally, simulations illustrate both benefits and caveats of co-association networks.

CONCLUSIONS

Our results support the prediction that different selection pressures favor the evolution of distinct groups of genes, each associating with a different aspect of climate. But our results went against the prediction that loci experiencing different sources of selection would have high recombination among them. These results give new insight into evolutionary debates about the extent of modularity, pleiotropy, and linkage in the evolution of genetic architectures.

Genetic differentiation in life history traits and thermal stress performance across a heterogeneous dune landscape in Arabidopsis lyrata

Background and Aims

Over very short spatial scales, the habitat of a species can differ in multiple abiotic and biotic factors. These factors may impose natural selection on several traits and can cause genetic differentiation within a population. We studied multivariate genetic differentiation in a plant species of a sand dune landscape by linking environmental variation with differences in genotypic trait values and gene expression levels to find traits and candidate genes of microgeographical adaptation.

Methods

Maternal seed families of Arabidopsis lyrata were collected in Saugatuck Dunes State Park, Michigan, USA, and environmental parameters were recorded at each collection site. Offspring plants were raised in climate chambers and exposed to one of three temperature treatments: regular occurrence of frost, heat, or constant control conditions. Several traits were assessed: plant growth, time to flowering, and frost and heat resistance.

Key Results

The strongest trait-environment association was between a fast switch to sexual reproduction and weaker growth under frost, and growing in the open, away from trees. The second strongest association was between the trait combination of small plant size and early flowering under control conditions combined with large size under frost, and the combination of environmental conditions of growing close to trees, at low vegetation cover, on dune bottoms. Gene expression analysis by RNA-seq revealed candidate genes involved in multivariate trait differentiation.

Conclusions

The results support the hypothesis that in natural populations, many environmental factors impose selection, and that they affect multiple traits, with the relative direction of trait change being complex. The results highlight that heterogeneity in the selection environment over small spatial scales is a main driver of the maintenance of adaptive genetic variation within populations.

Adaptations of white spruce to climate: strong intraspecific differences in cold hardiness linked to survival

Understanding local adaptation of tree populations to climate allows the development of assisted migration guidelines as a tool for forest managers to address climate change. Here, we study the relationship among climate, a wide range of physiological traits, and field performance of selected white spruce provenances originating from throughout the species range. Tree height, survival, cold hardiness, hydraulic, and wood anatomical traits were measured in a 32-year-old common garden trial, located in the center of the species range. Provenance performance included all combinations of high versus low survival and growth, with the most prevalent population differentiation for adaptive traits observed in cold hardiness. Cold hardiness showed a strong association with survival and was associated with cold winter temperatures at the site of seed origin. Tree height was mostly explained by the length of the growing season at the origin of the seed source. Although population differentiation was generally weak in wood anatomical and hydraulic traits, within-population variation was substantial in some traits, and a boundary analysis revealed that efficient water transport was associated with vulnerable xylem and low wood density, indicating that an optimal combination of high water transport efficiency and high cavitation resistance is not possible. Our results suggest that assisted migration prescriptions may be advantageous under warming climate, but pronounced trade-offs between survival and cold hardiness require a careful consideration of the distances of these transfers.

Bioinformatically predicted deleterious mutations reveal complementation in the interior spruce hybrid complex

BACKGROUND

Mutation load is expected to be reduced in hybrids via complementation of deleterious alleles. While local adaptation of hybrids confounds phenotypic tests for reduced mutation load, it may be possible to assess variation in load by analyzing the distribution of putatively deleterious alleles. Here, we use this approach in the interior spruce (Picea glauca x P. engelmannii) hybrid complex, a group likely to suffer from high mutation load and in which hybrids exhibit local adaptation to intermediate conditions. We used PROVEAN to bioinformatically predict whether non-synonymous alleles are deleterious, based on conservation of the position and abnormality of the amino acid change.

RESULTS

As expected, we found that predicted deleterious alleles were at lower average allele frequencies than alleles not predicted to be deleterious. We were unable to detect a phenotypic effect on juvenile growth rate of the many rare alleles predicted to be deleterious. Both the proportion of alleles predicted to be deleterious and the proportion of loci homozygous for predicted deleterious alleles were higher in P. engelmannii (Engelmann spruce) than in P. glauca (white spruce), due to higher diversity and frequencies of rare alleles in Engelmann. Relative to parental species, the proportion of alleles predicted to be deleterious was intermediate in hybrids, and the proportion of loci homozygous for predicted deleterious alleles was lowest.

CONCLUSION

Given that most deleterious alleles are recessive, this suggests that mutation load is reduced in hybrids due to complementation of deleterious alleles. This effect may enhance the fitness of hybrids.

Growth gains from selective breeding in a spruce hybrid zone do not compromise local adaptation to climate

Hybrid zones contain extensive standing genetic variation that facilitates rapid responses to selection. The Picea glauca × Picea engelmannii hybrid zone in western Canada is the focus of tree breeding programs that annually produce ~90 million reforestation seedlings. Understanding the direct and indirect effects of selective breeding on adaptive variation is necessary to implement assisted gene flow (AGF) polices in Alberta and British Columbia that match these seedlings with future climates. We decomposed relationships among hybrid ancestry, adaptive traits, and climate to understand the implications of selective breeding for climate adaptations and AGF strategies. The effects of selection on associations among hybrid index estimated from ~6,500 SNPs, adaptive traits, and provenance climates were assessed for ~2,400 common garden seedlings. Hybrid index differences between natural and selected seedlings within breeding zones were small in Alberta (average +2%), but larger and more variable in BC (average −7%, range −24% to +1%), slightly favoring P. glauca ancestry. The average height growth gain of selected seedlings over natural seedlings within breeding zones was 36% (range 12%–86%). Clines in growth with temperature‐related variables were strong, but differed little between selected and natural populations. Seedling hybrid index and growth trait associations with evapotranspiration‐related climate variables were stronger in selected than in natural seedlings, indicating possible preadaptation to drier future climates. Associations among cold hardiness, hybrid ancestry, and cold‐related climate variables dominated signals of local adaptation and were preserved in breeding populations. Strong hybrid ancestry–phenotype–climate associations suggest that AGF will be necessary to match interior spruce breeding populations with shifting future climates. The absence of antagonistic selection responses among traits and maintenance of cold adaptation in selected seedlings suggests breeding populations can be safely redeployed using AGF prescriptions similar to those of natural populations.