Potential for evolutionary responses to climate change - evidence from tree populations

A multidimensional selective landscape drives adaptive divergence between and within closely related Phlox species

Selection causes local adaptation across populations within species and simultaneously divergence between species. However, it is unclear if either the force of or the response to selection is similar across these scales. We show that natural selection drives divergence between closely related species in a pattern that is distinct from local adaptation within species. We use reciprocal transplant experiments across three species of Phlox wildflowers to characterize widespread adaptive divergence. Using provenance trials, we also find strong local adaptation between populations within a species. Comparing divergence and selection between these two scales of diversity we discover that one suite of traits predicts fitness differences between species and that an independent suite of traits predicts fitness variation within species. Selection drives divergence between species, contributing to speciation, while simultaneously favoring extensive diversity that is maintained across populations within a species. Our work demonstrates how the selection landscape is complex and multidimensional.

Pollen thermotolerance of a widespread plant, Lotus corniculatus, in response to climate warming: possible local adaptation of populations from different elevations

One of the most vulnerable phases in the plant life cycle is sexual reproduction, which depends on effective pollen transfer, but also on the thermotolerance of pollen grains. Pollen thermotolerance is temperature-dependent and may be reduced by increasing temperature associated with global warming. A growing body of research has focused on the effect of increased temperature on pollen thermotolerance in crops to understand the possible impact of temperature extremes on yield. Yet, little is known about the effects of temperature on pollen thermotolerance of wild plant species. To fill this gap, we selected Lotus corniculatus s.l. (Fabaceae), a species common to many European habitats and conducted laboratory experiments to test its pollen thermotolerance in response to artificial increase in temperature. To test for possible local adaptation of pollen thermal tolerance, we compared data from six lowland (389-451 m a.s.l.) and six highland (841-1,030 m a.s.l.) populations. We observed pollen germination in vitro at 15 °C, 25 °C, 30 °C, and 40 °C. While lowland plants maintained a stable germination percentage across a broad temperature range (15-30 °C) and exhibited reduced germination only at extremely high temperatures (40 °C), highland plants experienced reduced germination even at 30 °C-temperatures commonly exceeded in lowlands during warm summers. This suggests that lowland populations of L. corniculatus may be locally adapted to higher temperature for pollen germination. On the other hand, pollen tube length decreased with increasing temperature in a similar way in lowland and highland plants. The overall average pollen germination percentage significantly differed between lowland and highland populations, with highland populations displaying higher germination percentage. On the other hand, the average pollen tube length was slightly smaller in highland populations. In conclusion, we found that pollen thermotolerance of L. corniculatus is reduced at high temperature and that the germination of pollen from plant populations growing at higher elevations is more sensitive to increased temperature, which suggests possible local adaptation of pollen thermotolerance.

Under Stress: Searching for Genes Involved in the Response of Abies pinsapo Boiss to Climate Change

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

Genotype-by-environment interactions and local adaptation shape selection in the US National Chip Processing Trial

KEY MESSAGE

We find evidence of selection for local adaptation and extensive genotype-by-environment interaction in the potato National Chip Processing Trial (NCPT). We present a novel method for dissecting the interplay between selection, local adaptation and environmental response in plant breeding schemes. Balancing local adaptation and the desire for widely adapted cultivars is challenging for plant breeders and makes genotype-by-environment interactions (GxE) an important target of selection. Selecting for GxE requires plant breeders to evaluate plants across multiple environments. One way breeders have accomplished this is to test advanced materials across many locations. Public potato breeders test advanced breeding material in the National Chip Processing Trial (NCPT), a public-private partnership where breeders from ten institutions submit advanced chip lines to be evaluated in up to ten locations across the country. These clones are genotyped and phenotyped for important agronomic traits. We used these data to interrogate the NCPT for GxE. Further, because breeders submitting clones to the NCPT select in a relatively small geographic range for the first 3 years of selection, we examined these data for evidence of incidental selection for local adaptation, and the alleles underlying it, using an environmental genome-wide association study (envGWAS). We found genomic regions associated with continuous environmental variables and discrete breeding programs, as well as regions of the genome potentially underlying GxE for yield.

Response of Extremely Small Populations to Climate Change-A Case of Trachycarpus nanus in Yunnan, China

Climate change affects the geographical distribution of plant species. Rare Trachycarpus nanus with a narrow distribution range, high medicinal value and extremely small population is facing increasing extinction risks under global climate change. In this study, 96 recorded occurrences and 23 environmental factors are used to predict the potential suitable area of T. nanus based on the optimized MaxEnt (3.4.4) model and ArcGIS (10.7) software. The results show that when the parameters are FC = LQ and RM = 1, the MaxEnt model is optimal and AUC = 0.946. The distribution patterns were predicted in the past, present, and four future phases, i.e., 2021-2040 (2030), 2041-2060 (2050), 2061-2080 (2070), and 2081-2100 (2090). The main factors are the annual precipitation (bio12), mean temperature of the coldest quarter (bio11), temperature seasonality (bio4), precipitation of the wettest quarter (bio16), and isothermality (bio3). The potential distribution of T. nanus is primarily concentrated in central Chuxiong, encompassing a total potential suitable area of 5.65 × 104 km2. In historical periods, the total habitat area is smaller than that in the present. In the future, the potential suitable area is generally increased. The centroid analysis shows that T. nanus will move to a high-altitude area and to the southeast. But its dispersal capacity may not keep up with the climate change rate. Therefore, additional protection sites for this species should be appropriately established and the habitat connectivity should be enhanced.

How useful are genomic data for predicting maladaptation to future climate?

Methods using genomic information to forecast potential population maladaptation to climate change or new environments are becoming increasingly common, yet the lack of model validation poses serious hurdles toward their incorporation into management and policy. Here, we compare the validation of maladaptation estimates derived from two methods-Gradient Forests (GFoffset) and the risk of non-adaptedness (RONA)-using exome capture pool-seq data from 35 to 39 populations across three conifer taxa: two Douglas-fir varieties and jack pine. We evaluate sensitivity of these algorithms to the source of input loci (markers selected from genotype-environment associations [GEA] or those selected at random). We validate these methods against 2- and 52-year growth and mortality measured in independent transplant experiments. Overall, we find that both methods often better predict transplant performance than climatic or geographic distances. We also find that GFoffset and RONA models are surprisingly not improved using GEA candidates. Even with promising validation results, variation in model projections to future climates makes it difficult to identify the most maladapted populations using either method. Our work advances understanding of the sensitivity and applicability of these approaches, and we discuss recommendations for their future use.

Weak evidence of provenance effects in spring phenology across Europe and North America

Forecasting the biological impacts of climate change requires understanding how species respond to warmer temperatures through interannual flexible variation vs through adaptation to local conditions. Yet, we often lack this information entirely or find conflicting evidence across studies, which is the case for spring phenology. We synthesized common garden studies across Europe and North America that reported spring event dates for a mix of angiosperm and gymnosperm tree species in the northern hemisphere, capturing data from 384 North American and 101 European provenances (i.e. populations) with observations from 1962 to 2019, alongside autumn event data when provided. Across continents, we found no evidence of provenance effects in spring phenology, but strong clines with latitude and mean annual temperature in autumn. These effects, however, appeared to diverge by continent and species type (gymnosperm vs angiosperm), with particularly pronounced clines in North America in autumn events. Our results suggest flexible, likely plastic responses, in spring phenology with warming, and potential limits - at least in the short term - due to provenance effects for autumn phenology. They also highlight that, after over 250 yr of common garden studies on tree phenology, we still lack a holistic predictive model of clines across species and phenological events.

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.

Spring wood phenology responds more strongly to chilling temperatures than bud phenology in European conifers

A comparative assessment of bud and wood phenology could aid a better understanding of tree growth dynamics. However, the reason for asynchronism or synchronism in leaf and cambial phenology remains unclear. To test the assumption that the temporal relationship between the budburst date and the onset date of wood formation is due to their common or different responses to environmental factors, we constructed a wood phenology dataset from previous literature, and compared it with an existing bud phenology dataset in Europe. We selected three common conifers (Larix decidua Mill., Picea abies (L.) H. Karst. and Pinus sylvestris L.) in both datasets and analyzed 909 records of the onset of wood formation at 47 sites and 238,720 records of budburst date at 3051 sites. We quantified chilling accumulation (CA) and forcing requirement (FR) of budburst and onset of wood formation based on common measures of CA and FR. We then constructed negative exponential CA-FR curves for bud and wood phenology separately. The results showed that the median, variance and probability distribution of CA-FR curves varied significantly between bud and wood phenology for three conifers. The different FR under the same chilling condition caused asynchronous bud and wood phenology. Furthermore, the CA-FR curves manifested that wood phenology was more sensitive to chilling than bud phenology. Thus, the FR of the onset of wood formation increases more than that of budburst under the same warming scenarios, explaining the stronger earlier trends in the budburst date than the onset date of woody formation simulated by the process-based model. Our work not only provides a possible explanation for asynchronous bud and wood phenology from the perspective of organ-specific responses to chilling and forcing, but also develops a phenological model for predicting both bud and wood phenology with acceptable uncertainties.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Alpine plant species converge towards adopting elevation-specific resource-acquisition strategy in response to experimental early snow-melting

Snow-melt is one of the important factors limiting growth and survival of alpine plants. Changes in snow-melt timing have profound effects on eco-physiological characteristics of alpine plant species through alterations in growing season length. Here, we conducted a field experiment and studied species response to experimentally induced early snow-melting (ES) (natural vs. early) at an alpine site (Rohtang) in the western Himalaya region. Eco-physiological response of eight snow-bed restricted alpine plant species from different elevations (lower: 3850 m and upper: 4150 m amsl) and belonging to contrasting resource acquisition strategies (conservative and acquisitive) were studied after 2-years (2019 & 2020) of initiating ES field experiment. We estimated the functional traits related to leaf economic spectrum and physiological performance and assessed their pattern of phenotypic plasticity. Analysis by linear mixed effect model showed that both the 'conservative' and 'acquisitive' species had responded to ES with significant effects on species specific leaf area, leaf dry matter content, leaf thickness, leaf water content and sugar content. Our results also revealed that ES treatment induced significant increase in leaf C/N ratio (10.57 % to 13.65 %) and protein content (15.85 % to 20.76 %) at both the elevations, irrespective of species groups. The phenotypic plasticity was found to be low and was essentially species-specific. However, for leaf protein content, the upper elevation species exhibited a higher phenotypic plasticity (0.43 ± 0.18) than the lower elevation species (0.31 ± 0.21). Interestingly, we found that irrespective of species unique functional strategy, species adapt to perform more conservative at lower elevation and more acquisitive at upper elevation, in response to ES. We conclude that plants occurring at contrasting elevations respond differentially to ES. However, species showed capacity for short-term acclimation to future environmental conditions, but may be vulnerable, if their niche is occupied by new species with greater phenotypic plasticity and a superior competitive ability.

Adaptive potential of maritime pine under contrasting environments

BACKGROUND

Predicting the adaptability of forest tree populations under future climates requires a better knowledge of both the adaptive significance and evolvability of measurable key traits. Phenotypic plasticity, standing genetic variation and degree of phenotypic integration shape the actual and future population genetic structure, but empirical estimations in forest tree species are still extremely scarce. We analysed 11 maritime pine populations covering the distribution range of the species (119 families and 8 trees/family, ca. 1300 trees) in a common garden experiment planted at two sites with contrasting productivity. We used plant height as a surrogate of fitness and measured five traits (mean and plasticity of carbon isotope discrimination, specific leaf area, needle biomass, Phenology growth index) related to four different strategies (acquisitive economics, photosynthetic organ size, growth allocation and avoidance of water stress).

RESULTS

Estimated values of additive genetic variation would allow adaptation of the populations to future environmental conditions. Overall phenotypic integration and selection gradients were higher at the high productivity site, while phenotypic integration within populations was higher at the low productivity site. Response to selection was related mainly to photosynthetic organ size and drought-avoidance mechanisms rather than to water use efficiency. Phenotypic plasticity of water use efficiency could be maladaptive, resulting from selection for height growth.

CONCLUSIONS

Contrary to the expectations in a drought tolerant species, our study suggests that variation in traits related to photosynthetic organ size and acquisitive investment of resources drive phenotypic selection across and within maritime pine populations. Both genetic variation and evolvability of key adaptive traits were considerably high, including plasticity of water use efficiency. These characteristics would enable a relatively fast micro-evolution of populations in response to the ongoing climate changes. Moreover, differentiation among populations in the studied traits would increase under the expected more productive future Atlantic conditions.

A species' response to spatial climatic variation does not predict its response to climate change

The dominant paradigm for assessing ecological responses to climate change assumes that future states of individuals and populations can be predicted by current, species-wide performance variation across spatial climatic gradients. However, if the fates of ecological systems are better predicted by past responses to in situ climatic variation through time, this current analytical paradigm may be severely misleading. Empirically testing whether spatial or temporal climate responses better predict how species respond to climate change has been elusive, largely due to restrictive data requirements. Here, we leverage a newly collected network of ponderosa pine tree-ring time series to test whether statistically inferred responses to spatial versus temporal climatic variation better predict how trees have responded to recent climate change. When compared to observed tree growth responses to climate change since 1980, predictions derived from spatial climatic variation were wrong in both magnitude and direction. This was not the case for predictions derived from climatic variation through time, which were able to replicate observed responses well. Future climate scenarios through the end of the 21st century exacerbated these disparities. These results suggest that the currently dominant paradigm of forecasting the ecological impacts of climate change based on spatial climatic variation may be severely misleading over decadal to centennial timescales.

Effects of freezing temperatures on early life stages of native trees of different elevational origin: implications for tree recruitment in seasonally dry mountain forests

In mountain forests, tree regeneration is limited by increasingly frequent frosts with increasing elevation. We investigated the effects of exposure to freezing temperature on early life stages of two native trees of different elevational origin in a seasonally dry mountain forest. We hypothesized that the negative effects of freezing exposure on performance of early life stages increases as freezing temperature decreases, and that frost resistance increases in plants of high elevational origin. We collected seeds of two tree species (Kageneckia lanceolata and Lithraea molleoides) from populations located at different elevations and grew seedlings and saplings in a greenhouse. Dry seeds, imbibed seeds and 1-month-old seedlings were exposed to seven temperature treatments ranging from 4 °C to -20 °C, while 12-month-old saplings were exposed to four temperature treatments from -8 °C to -20 °C. After freezing exposure in a climate chamber, we monitored seed germination and seedling and sapling survival. Germination of K. lanceolata decreased with decreasing temperature only for imbibed seeds from mid- and high elevations, whereas germination of L. molleoides slightly increased with decreasing temperature only for imbibed seeds from high elevations. For both species, seedling survival decreased with decreasing temperature. For K. lanceolata, the negative effects of freezing temperatures were weaker as elevational origin of seeds increased, whereas L. molleoides showed the opposite pattern. For both species, saplings only survived at the mildest applied freezing temperature (-8 °C). We conclude that effects of climatic variation associated with elevation depend on the study species and life stage. The observed patterns could be caused by maternal effects, which are absent at the sapling stage. Moreover, temperatures below -8 °C can limit recruitment since partial mortality of seedlings and saplings occurred at such values.

Projection of current and future distribution of adaptive genetic units in an alpine ungulate

Climate projections predict major changes in alpine environments by the end of the 21st century. To avoid climate-induced maladaptation and extinction, many animal populations will either need to move to more suitable habitats or adapt in situ to novel conditions. Since populations of a species exhibit genetic variation related to local adaptation, it is important to incorporate this variation into predictive models to help assess the ability of the species to survive climate change. Here, we evaluate how the adaptive genetic variation of a mountain ungulate-the Northern chamois (Rupicapra rupicapra)-could be impacted by future global warming. Based on genotype-environment association analyses of 429 chamois using a ddRAD sequencing approach, we identified genetic variation associated with climatic gradients across the European Alps. We then delineated adaptive genetic units and projected the optimal distribution of these adaptive groups in the future. Our results suggest the presence of local adaptation to climate in Northern chamois with similar genetic adaptive responses in geographically distant but climatically similar populations. Furthermore, our results predict that future climatic changes will modify the Northern chamois adaptive landscape considerably, with various degrees of maladaptation risk.

Low but significant evolutionary potential for growth, phenology and reproduction traits in European beech

Local survival of forest tree populations under climate change depends on existing genetic variation and their adaptability to changing environments. Responses to selection were studied in European beech (Fagus sylvatica) under field conditions. A total of 1087 adult trees, seeds, 1-year-old seedlings and established multiyear saplings were genotyped with 16 nuSSRs. Adult trees were assessed for phenotypic traits related to growth, phenology and reproduction. Parentage and paternity analyses were used to estimate effective female and male fecundity as a proxy of fitness and showed that few parents contributed to successful regeneration. Selection gradients were estimated from the relationship between traits and fecundity, while heritability and evolvability were estimated using mixed models and the breeder's equation. Larger trees bearing more fruit and early male flowering had higher total fecundity, while trees with longer growth season had lower total fecundity (directional selection). Stabilizing selection on spring phenology was found for female fecundity, highlighting the role of late frosts as a selection driver. Selection gradients for other traits varied between measurement years and the offspring cohort used to estimate parental fecundity. Compared to other studies in natural populations, we found low to moderate heritability and evolvability for most traits. Response to selection was higher for growth than for budburst, leaf senescence or reproduction traits, reflecting more consistent selection gradients across years and sex functions, and higher phenotypic variability in the population. Our study provides empirical evidence suggesting that populations of long-lived organisms such as forest trees can adapt locally, even at short-time scales.

Adaptation to climate change through dispersal and inherited timing in an avian migrant

Many organisms fail to adjust their phenology sufficiently to climate change. Studies have concentrated on adaptive responses within localities, but little is known about how latitudinal dispersal enhances evolutionary potential. Rapid adaptation is expected if dispersers from lower latitudes have improved synchrony to northern conditions, thereby gain fitness and introduce genotypes on which selection acts. Here we provide experimental evidence that dispersal in an avian migrant enables rapid evolutionary adaptation. We translocated Dutch female pied flycatchers (Ficedula hypoleuca) and eggs to Sweden, where breeding phenology is ~15 days later. Translocated females bred earlier, and their fitness was 2.5 times higher than local Swedish flycatchers. We show that between-population variation in timing traits is highly heritable, and hence immigration of southern genotypes promotes the necessary evolutionary response. We conclude that studies on adaptation to large-scale environmental change should not just focus on plasticity and evolution based on standing genetic variation but should also include phenotype-habitat matching through dispersal as a viable route to adjust.

A demo-genetic model shows how silviculture reduces natural density-dependent selection in tree populations

Biological production systems and conservation programs benefit from and should care for evolutionary processes. Developing evolution-oriented strategies requires knowledge of the evolutionary consequences of management across timescales. Here, we used an individual-based demo-genetic modelling approach to study the interactions and feedback between tree thinning, genetic evolution, and forest stand dynamics. The model combines processes that jointly drive survival and mating success-tree growth, competition and regeneration-with genetic variation of quantitative traits related to these processes. In various management and disturbance scenarios, the evolutionary rates predicted by the coupled demo-genetic model for a growth-related trait, vigor, fit within the range of empirical estimates found in the literature for wild plant and animal populations. We used this model to simulate non-selective silviculture and disturbance scenarios over four generations of trees. We characterized and quantified the effect of thinning frequencies and intensities and length of the management cycle on viability selection driven by competition and fecundity selection. The thinning regimes had a drastic long-term effect on the evolutionary rate of vigor over generations, potentially reaching 84% reduction, depending on management intensity, cycle length and disturbance regime. The reduction of genetic variance by viability selection within each generation was driven by changes in genotypic frequencies rather than by gene diversity, resulting in low-long-term erosion of the variance across generations, despite short-term fluctuations within generations. The comparison among silviculture and disturbance scenarios was qualitatively robust to assumptions on the genetic architecture of the trait. Thus, the evolutionary consequences of management result from the interference between human interventions and natural evolutionary processes. Non-selective thinning, as considered here, reduces the intensity of natural selection, while selective thinning (on tree size or other criteria) might reduce or reinforce it depending on the forester's tree choice and thinning intensity.

Comparing DNA isolation methods for forest trees: quality, plastic footprint, and time-efficiency

BACKGROUND

Genetic and genomic studies are seeing an increase in sample sizes together with a wider range of species investigated in response to environmental change concerns. In turn, these changes may come with challenges including the time and difficulty to isolate nucleic acids (DNA or RNA), the sequencing cost and environmental impacts of the growing amount of plastic waste generated in the process. Pseudotsuga menziesii var. menziesii (Mirbel) Franco (PM), Tsuga heterophylla (Raf.) Sarg. (TH) and Thuja plicata Donn ex D.Don (TP) are conifer species found in diverse woodlands both as natives and naturalized exotics. Our study was carried out whilst investigating their genetics to understand their population structure and potential for adaptation.

RESULTS

In the present study, we compared two different DNA isolation methods, i.e., spin-column DNeasy plant mini kit (QIAGEN), and temperature-driven enzymatic cocktail Plant DNA Extraction (MicroGEM). The quantity of recovered DNA and the quality of DNA were assessed along with the plastic footprint and time needed for three tree species. Both methods were optimised and proven to provide enough DNA for each studied species. The yield of DNA for each method depended on the species: QIAGEN showed higher yield in P. menziesii and T. heterophylla, while T. plicata recovered similar amount of DNA for both methods. The DNA quality was investigated using DNA barcoding techniques by confirming species identity and species discrimination. No difference was detected in the PCR amplification of the two barcoding loci, (rbcL and trnH-psbA), and the recovered sequences between DNA isolation methods. Measurement of the plastic use and the processing time per sample indicated that MicroGEM had a 52.64% lower plastic footprint and was 51.8% faster than QIAGEN.

CONCLUSIONS

QIAGEN gave higher yields in two of the species although both methods showed similar quality results across all species. However, MicroGEM was clearly advantageous to decrease the plastic footprint and improve the time efficiency. Overall, MicroGEM recovers sufficient and reliable DNA to perform common downstream analyses such as PCR and sequencing. Our findings illustrate the benefits of research and efforts towards developing more sustainable methods and techniques to reduce the environmental footprint of molecular analyses.

Identifying genetic variation associated with environmental gradients and drought-tolerance phenotypes in ponderosa pine

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.

Picturing local adaptation: Spectral and structural traits from drone remote sensing reveal clinal responses to climate transfer in common-garden trials of interior spruce (Picea engelmannii × glauca)

Common-garden trials of forest trees provide phenotype data used to assess growth and local adaptation; this information is foundational to tree breeding programs, genecology, and gene conservation. As jurisdictions consider assisted migration strategies to match populations to suitable climates, in situ progeny and provenance trials provide experimental evidence of adaptive responses to climate change. We used drone technology, multispectral imaging, and digital aerial photogrammetry to quantify spectral traits related to stress, photosynthesis, and carotenoids, and structural traits describing crown height, size, and complexity at six climatically disparate common-garden trials of interior spruce (Picea engelmannii × glauca) in western Canada. Through principal component analysis, we identified key components of climate related to temperature, moisture, and elevational gradients. Phenotypic clines in remotely sensed traits were analyzed as trait correlations with provenance climate transfer distances along principal components (PCs). We used traits showing clinal variation to model best linear unbiased predictions for tree height (R2  = .98-.99, root mean square error [RMSE] = 0.06-0.10 m) and diameter at breast height (DBH, R2  = .71-.97, RMSE = 2.57-3.80 mm) and generated multivariate climate transfer functions with the model predictions. Significant (p < .05) clines were present for spectral traits at all sites along all PCs. Spectral traits showed stronger clinal variation than structural traits along temperature and elevational gradients and along moisture gradients at wet, coastal sites, but not at dry, interior sites. Spectral traits may capture patterns of local adaptation to temperature and montane growing seasons which are distinct from moisture-limited patterns in stem growth. This work demonstrates that multispectral indices improve the assessment of local adaptation and that spectral and structural traits from drone remote sensing produce reliable proxies for ground-measured height and DBH. This phenotyping framework contributes to the analysis of common-garden trials towards a mechanistic understanding of local adaptation to climate.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Climate and Soil Microsite Conditions Determine Local Adaptation in Declining Silver Fir Forests

Ongoing climatic change is threatening the survival of drought-sensitive tree species, such as silver fir (Abies alba). Drought-induced dieback had been previously explored in this conifer, although the role played by tree-level genetic diversity and its relationship with growth patterns and soil microsite conditions remained elusive. We used double digest restriction-site-associated DNA sequencing (ddRADseq) to describe different genetic characteristics of five silver fir forests in the Spanish Pyrenees, including declining and non-declining trees. Single nucleotide polymorphisms (SNPs) were used to investigate the relationships between genetics, dieback, intraspecific trait variation (functional dendrophenotypic traits and leaf traits), local bioclimatic conditions, and rhizosphere soil properties. While there were no noticeable genetic differences between declining and non-declining trees, genome-environment associations with selection signatures were abundant, suggesting a strong influence of climate, soil physicochemical properties, and soil microbial diversity on local adaptation. These results provide novel insights into how genetics and diverse environmental factors are interrelated and highlight the need to incorporate genetic data into silver fir forest dieback studies to gain a better understanding of local adaptation.

Germination temperature sensitivity differs between co-occurring tree species and climate origins resulting in contrasting vulnerability to global warming

Climate change is shifting temperatures from historical patterns, globally impacting forest composition and resilience. Seed germination is temperature-sensitive, making the persistence of populations and colonization of available habitats vulnerable to warming. This study assessed germination response to temperature in foundation trees in south-western Australia's Mediterranean-type climate forests (Eucalyptus marginata (jarrah) and Corymbia calophylla (marri)) to estimate the thermal niche and vulnerability among populations. Seeds from the species' entire distribution were collected from 12 co-occurring populations. Germination thermal niche was investigated using a thermal gradient plate (5-40°C). Five constant temperatures between 9 and 33°C were used to test how the germination niche (1) differs between species, (2) varies among populations, and (3) relates to the climate of origin. Germination response differed among species; jarrah had a lower optimal temperature and thermal limit than marri (T _o 15.3°C, 21.2°C; ED₅₀ 23.4°C, 31°C, respectively). The thermal limit for germination differed among populations within both species, yet only marri showed evidence for adaptation to thermal origins. While marri has the capacity for germination at higher thermal temperatures, jarrah is more vulnerable to global warming exceeding safety margins. This discrepancy is predicted to alter species distributions and forest composition in the future.

Reduced within-population quantitative genetic variation is associated with climate harshness in maritime pine

How evolutionary forces interact to maintain genetic variation within populations has been a matter of extensive theoretical debates. While mutation and exogenous gene flow increase genetic variation, stabilizing selection and genetic drift are expected to deplete it. To date, levels of genetic variation observed in natural populations are hard to predict without accounting for other processes, such as balancing selection in heterogeneous environments. We aimed to empirically test three hypotheses: (i) admixed populations have higher quantitative genetic variation due to introgression from other gene pools, (ii) quantitative genetic variation is lower in populations from harsher environments (i.e., experiencing stronger selection), and (iii) quantitative genetic variation is higher in populations from heterogeneous environments. Using growth, phenological and functional trait data from three clonal common gardens and 33 populations (522 clones) of maritime pine (Pinus pinaster Aiton), we estimated the association between the population-specific total genetic variances (i.e., among-clone variances) for these traits and ten population-specific indices related to admixture levels (estimated based on 5165 SNPs), environmental temporal and spatial heterogeneity and climate harshness. Populations experiencing colder winters showed consistently lower genetic variation for early height growth (a fitness-related trait in forest trees) in the three common gardens. Within-population quantitative genetic variation was not associated with environmental heterogeneity or population admixture for any trait. Our results provide empirical support for the potential role of natural selection in reducing genetic variation for early height growth within populations, which indirectly gives insight into the adaptive potential of populations to changing environments.

Climate change decouples dominant tree species in African savannas

To understand how two dominant African savanna trees will continue to respond to climate changes, we examined their regeneration niche and adult tree distributions. Specifically, we wanted to (1) determine if distributional patterns were shifting, (2) predict future distributions under different climate change scenarios and (3) evaluate the realism of predicted future distributions. We randomly placed 40 grids into 6 strata across a climate gradient in the kingdom of Eswatini. Within these grids, we sampled adult and seedling marula (Scelerocarya birrea) and knobthorn (Senegalia nigrecens) trees and used the data to model their abundance. Next, we quantified shifts in distributional patterns (e.g., expansion or contraction) by measuring the current and projected areas of overlap between seedling and adult trees. Finally, we predicted future distributions of abundance based on predicted climate conditions. We found knobthorn seedlings within a small portion of the adult distribution, suggesting it was unlikely to track climate changes. Alternatively, finding marula seedlings on and beyond one edge of the adult distribution, suggested its range would shift toward cooler climates. Predicted future distributions suggest suitable climate for both species would transition out of savannas and into grasslands. Future projections (2041-2070) appeared consistent with observed distributions of marula, but knobthorn predictions were unrealistic given the lack of evidence for regeneration outside of its current range. The idiosyncratic responses of these species to climate change are likely to decouple these keystone structures in the coming decades and are likely to have considerable cascading effects including the potential rearrangement of faunal communities.

Microgeographic variation in early fitness traits of Pinus sylvestris from contrasting soils

PREMISE

The possibility of fine-scale intraspecific adaptive divergence under gene flow is established by theoretical models and has been confirmed empirically in tree populations distributed along steep altitudinal clines or across extreme edaphic discontinuities. However, the possibility of microgeographic adaptive divergence due to less severe but more frequent kinds of soil variation is unclear.

METHODS

In this study, we looked for evidence of local adaptation to calcareous versus siliceous soil types in two nearby Mediterranean Pinus sylvestris populations connected via pollen flow. Using a greenhouse experiment, we tested for variation in early (up to three years of age) seedling performance among open-pollinated maternal families originating from each edaphic provenance when experimentally grown on both types of natural local substrate.

RESULTS

Although seedlings were clearly affected by the edaphic environment, exhibiting lower and slower emergence as well as higher mortality on the calcareous than in the siliceous substrate, neither the performance on each substrate nor the plasticity among substrates varied significantly with seedling edaphic provenance.

CONCLUSIONS

We found no evidence of local adaptation to a non-extreme edaphic discontinuity over a small spatial scale, at least during early stages of seedling establishment. Future studies on microgeographic soil-driven adaptation should consider long-term experiments to minimize maternal effects and allow a potentially delayed expression of edaphic adaptive divergence.

Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.

Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

From common gardens to candidate genes: exploring local adaptation to climate in red spruce

Local adaptation to climate is common in plant species and has been studied in a range of contexts, from improving crop yields to predicting population maladaptation to future conditions. The genomic era has brought new tools to study this process, which was historically explored through common garden experiments. In this study, we combine genomic methods and common gardens to investigate local adaptation in red spruce and identify environmental gradients and loci involved in climate adaptation. We first use climate transfer functions to estimate the impact of climate change on seedling performance in three common gardens. We then explore the use of multivariate gene-environment association methods to identify genes underlying climate adaptation, with particular attention to the implications of conducting genome scans with and without correction for neutral population structure. This integrative approach uncovered phenotypic evidence of local adaptation to climate and identified a set of putatively adaptive genes, some of which are involved in three main adaptive pathways found in other temperate and boreal coniferous species: drought tolerance, cold hardiness, and phenology. These putatively adaptive genes segregated into two 'modules' associated with different environmental gradients. This study nicely exemplifies the multivariate dimension of adaptation to climate in trees.

Molecular insights into the dynamics of species invasion by hybridisation in Tasmanian eucalypts

In plants where seed dispersal is limited compared with pollen dispersal, hybridisation may enhance gene exchange and species dispersal. We provide genetic evidence of hybridisation contributing to the expansion of the rare Eucalyptus risdonii into the range of the widespread Eucalyptus amygdalina. These closely related tree species are morphologically distinct, and observations suggest that natural hybrids occur along their distribution boundaries and as isolated trees or in small patches within the range of E. amygdalina. Hybrid phenotypes occur outside the range of normal dispersal for E. risdonii seed, yet in some hybrid patches small individuals resembling E. risdonii occur and are hypothesised to be a result of backcrossing. Using 3362 genome-wide SNPs assessed from 97 individuals of E. risdonii and E. amygdalina and 171 hybrid trees, we show that (i) isolated hybrids match the genotypes expected of F₁ /F₂ hybrids, (ii) there is a continuum in the genetic composition among the isolated hybrid patches from patches dominated by F₁ /F₂ -like genotypes to those dominated by E. risdonii-backcross genotypes, and (iii) the E. risdonii-like phenotypes in the isolated hybrid patches are most-closely related to proximal larger hybrids. These results suggest that the E. risdonii phenotype has been resurrected in isolated hybrid patches established from pollen dispersal, providing the first steps in its invasion of suitable habitat by long-distance pollen dispersal and complete introgressive displacement of E. amygdalina. Such expansion accords with the population demographics, common garden performance data, and climate modelling which favours E. risdonii and highlights a role of interspecific hybridisation in climate change adaptation and species expansion.

Neutral and adaptive genetic diversity in plants: An overview

Genetic diversity is a prerequisite for evolutionary change in all kinds of organisms. It is generally acknowledged that populations lacking genetic variation are unable to evolve in response to new environmental conditions (e.g., climate change) and thus may face an increased risk of extinction. Although the importance of incorporating genetic diversity into the design of conservation measures is now well understood, less attention has been paid to the distinction between neutral (NGV) and adaptive (AGV) genetic variation. In this review, we first focus on the utility of NGV by examining the ways to quantify it, reviewing applications of NGV to infer ecological and evolutionary processes, and by exploring its utility in designing conservation measures for plant populations and species. Against this background, we then summarize the ways to identify and estimate AGV and discuss its potential use in plant conservation. After comparing NGV and AGV and considering their pros and cons in a conservation context, we conclude that there is an urgent need for a better understanding of AGV and its role in climate change adaptation. To date, however, there are only a few AGV studies on non-model plant species aimed at deciphering the genetic and genomic basis of complex trait variation. Therefore, conservation researchers and practitioners should keep utilizing NGV to develop relevant strategies for rare and endangered plant species until more estimates of AGV are available.

Tree-Level Growth Patterns and Genetic Associations Depict Drought Legacies in the Relict Forests of Abies marocana

The frequency and intensity of drought events are increasing worldwide, challenging the adaptive capacity of several tree species. Here, we evaluate tree growth patterns and climate sensitivity to precipitation, temperature, and drought in the relict Moroccan fir Abies marocana. We selected two study sites, formerly stated as harboring contrasting A. marocana taxa (A. marocana and A. tazaotana, respectively). For each tree, dendrochronological methods were applied to quantify growth patterns and climate-growth sensitivity. Further, ddRAD-seq was performed on the same trees and close saplings to obtain single nucleotide polymorphisms (SNPs) and related genotype-phenotype associations. Genetic differentiation between the two studied remnant populations of A. marocana was weak. Growth patterns and climate-growth relationships were almost similar at the two sites studied, supporting a negative effect of warming. Growth trends and tree size showed associations with SNPs, although there were no relationships with phenotypes related to climatic sensitivity. We found significant differences in the SNPs subjected to selection in the saplings compared to the old trees, suggesting that relict tree populations might be subjected to genetic differentiation and local adaptation to climate dryness. Our results illustrate the potential of tree rings and genome-wide analysis to improve our understanding of the adaptive capacity of drought-sensitive forests to cope with ongoing climate change.

Responses to Drought Stress in Poplar: What Do We Know and What Can We Learn?

Poplar (Populus spp.) is a high-value crop for wood and biomass production and a model organism for tree physiology and genomics. The early release, in 2006, of the complete genome sequence of P. trichocarpa was followed by a wealth of studies that significantly enriched our knowledge of complex pathways inherent to woody plants, such as lignin biosynthesis and secondary cell wall deposition. Recently, in the attempt to cope with the challenges posed by ongoing climate change, fundamental studies and breeding programs with poplar have gradually shifted their focus to address the responses to abiotic stresses, particularly drought. Taking advantage from a set of modern genomic and phenotyping tools, these studies are now shedding light on important processes, including embolism formation (the entry and expansion of air bubbles in the xylem) and repair, the impact of drought stress on biomass yield and quality, and the long-term effects of drought events. In this review, we summarize the status of the research on the molecular bases of the responses to drought in poplar. We highlight how this knowledge can be exploited to select more tolerant genotypes and how it can be translated to other tree species to improve our understanding of forest dynamics under rapidly changing environmental conditions.

Contrasting transcriptomic patterns reveal a genomic basis for drought resilience in the relict fir Abies pinsapo Boiss

Climate change challenges the adaptive capacity of several forest tree species in the face of increasing drought and rising temperatures. Therefore, understanding the mechanistic connections between genetic diversity and drought resilience is highly valuable for conserving drought-sensitive forests. Nonetheless, the post-drought recovery in trees from a transcriptomic perspective has not yet been studied by comparing contrasting phenotypes. Here, experimental drought treatments, gas-exchange dynamics and transcriptomic analysis (RNA-seq) were performed in the relict and drought-sensitive fir Abies pinsapo Boiss. to identify gene expression differences over immediate (24 h) and extended drought (20 days). Post-drought responses were investigated to define resilient and sensitive phenotypes. Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of A. pinsapo drought resilience. Weighted gene co-expression network analysis showed an activation of stomatal closing and an inhibition of plant growth-related genes during the immediate drought, consistent with an isohydric dynamic. During the extended drought, transcription factors, as well as cellular damage and homeostasis protection-related genes prevailed. Resilient individuals activate photosynthesis-related genes and inhibit aerial growth-related genes, suggesting a shifting shoot/root biomass allocation to improve water uptake and whole-plant carbon balance. About, 152 fixed SNPs were found between resilient and sensitive seedlings, which were mostly located in RNA-activity-related genes, including epigenetic regulation. Contrasting gene expression and SNPs were found between different post-drought resilience phenotypes for the first time in a forest tree, suggesting a transcriptomic and genomic basis for drought resilience. The obtained drought-related transcriptomic profile and drought-resilience candidate genes may guide conservation programs for this threatened tree species.

Weak genetic differentiation but strong climate-induced selective pressure toward the rear edge of mountain pine in north-eastern Spain

Local differentiation at distribution limits may influence species' adaptive capacity to environmental changes. However, drivers, such gene flow and local selection, are still poorly understood. We focus on the role played by range limits in mountain forests to test the hypothesis that relict tree populations are subjected to genetic differentiation and local adaptation. Two alpine treelines of mountain pine (Pinus uncinata Ram. ex DC) were investigated in the Spanish Pyrenees. Further, an isolated relict population forming the species' southernmost distribution limit in north-eastern Spain was also investigated. Using genotyping by sequencing, a genetic matrix conformed by single nucleotide polymorphisms (SNPs) was obtained. This matrix was used to perform genotype-environment and genotype-phenotype associations, as well as to model risk of non-adaptedness. Increasing climate seasonality appears as an essential element in the interpretation of SNPs subjected to selective pressures. Genetic differentiations were overall weak. The differences in leaf mass area and radial growth rate, as well as the identification of several SNPs subjected to selective pressures, exceeded neutral predictions of differentiation among populations. Despite genetic drift might prevail in the isolated population, the Fst values (0.060 and 0.066) showed a moderate genetic drift and Nm values (3.939 and 3.555) indicate the presence of gene flow between the relict population and both treelines. Nonetheless, the SNPs subjected to selection pressures provide evidences of possible selection in treeline ecotones. Persistence in range boundaries seems to involve several selective pressures in species' traits, which were significantly related to enhanced drought seasonality at the limit of P. uncinata distribution range. We conclude that gene flow is unlikely to constrain adaptation in the P. uncinata rear edge, although this species shows vulnerability to future climate change scenarios involving warmer and drier conditions.

Genomic signals of local adaptation and hybridization in Asian white birch

Disentangling the numerous processes that affect patterns of genome-wide diversity in widespread tree species has important implications for taxonomy, conservation, and forestry. Here, we investigate the population genomic structure of Asian white birch (Betula platyphylla) in China and seek to explain it in terms of hybridization, demography and adaptation. We generate whole genome sequence data from 83 individuals across the species range in China. Combining this with an existing data set for 79 European and Russian white birches, we show a clear distinction between B. pendula and B. platyphylla, which have sometimes been lumped taxonomically. Genomic diversity of B. platyphylla in north-western China and Central Russia is affected greatly by hybridization with B. pendula. Excluding these hybridized populations, B. platyphylla in China has a linear distribution from north-eastern to south-western China, along the edge of the inland mountainous region. Within this distribution, three genetic clusters are found, which we model as long diverged with subsequent episodes of gene flow. Patterns of covariation between allele frequencies and environmental variables in B. platyphylla suggest the role of natural selection in the distribution of diversity at 7609 SNPs of which 3767 were significantly differentiated among the genetic clusters. The putative adaptive SNPs are distributed throughout the genome and span 1633 genic regions. Of these genic regions, 87 were previously identified as candidates for selective sweeps in Eurasian B. pendula. We use the 7609 environmentally associated SNPs to estimate the risk of nonadaptedness for each sequenced B. platyphylla individual under a scenario of future climate change, highlighting areas where populations may be under future threat from rising temperatures.

De novo transcriptome sequencing and gene co-expression reveal a genomic basis for drought sensitivity and evidence of a rapid local adaptation on Atlas cedar (Cedrus atlantica)

Introduction

Understanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.

Methods

Here we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.

Results

De novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.

Discussion

This work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

Divergent selection predating the Last Glacial Maximum mainly acted on macro-phenotypes in Norway spruce

The current distribution and population structure of many species were, to a large extent, shaped by cycles of isolation in glacial refugia and subsequent population expansions. Isolation in and postglacial expansion through heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a constant interplay between evolutionary and demographic processes. We investigated population differentiation and adaptation of Norway spruce for juvenile growth, diameter of the stem, wood density, and tracheid traits at breast height. Data from 4461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of Q _ST vs. F _ST. We show that the macroscopic resultant trait (stem diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection that predated the Last Glacial Maximum. Altogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly preadapted, mainly through growth-related traits.

Local seed sourcing for sustainable forestry

Seed sourcing strategies are the basis for identifying genetic material meeting the requirements of future climatic conditions and social demands. Specifically, local seed sourcing has been extensively promoted, based on the expected adaptation of the populations to local conditions, but there are some limitations for the application. We analyzed Strict-sense local and Wide-sense local (based on climatic similarity) seed sourcing strategies. We determined species and genetic pools based on these strategies for 40 species and deployment zones in Spain. We also obtained the total number of seed sources and stands for these species in the EU countries. We analyzed the richness of the pools, the relationship with variables related to the use of the species in afforestation, and the availability of seed production areas approved for the production of reproductive material destined to be marketed. This study confirms the existence of extensive species and genetic local pools. Also, that the importance of these pools differs for different species, limitations being derived from the use of forest reproductive material and the existence of approved basic materials. Strategies derived from local seed sourcing approaches are the basis for the use of forest reproductive material because a large number of the species in the area considered in the study are under regulation. However, despite the extensive work done to approve basic materials, limitations based on the availability of seed production areas to provide local material for sustainable forestry are found in those species. Considering a Wide-sense local seed sourcing strategy we provide alternative pools in order to meet social demands under the actual regulations on marketing of reproductive materials.

Teasing apart the joint effect of demography and natural selection in the birth of a contact zone

Vast population movements induced by recurrent climatic cycles have shaped the genetic structure of plant species. During glacial periods species were confined to low-latitude refugia from which they recolonized higher latitudes as the climate improved. This multipronged recolonization led to many lineages that later met and formed large contact zones. We utilize genomic data from 5000 Picea abies trees to test for the presence of natural selection during recolonization and establishment of a contact zone in Scandinavia. Scandinavian P. abies is today made up of a southern genetic cluster originating from the Baltics, and a northern one originating from Northern Russia. The contact zone delineating them closely matches the limit between two major climatic regions. We show that natural selection contributed to its establishment and maintenance. First, an isolation-with-migration model with genome-wide linked selection fits the data better than a purely neutral one. Second, many loci show signatures of selection or are associated with environmental variables. These loci, regrouped in clusters on chromosomes, are often related to phenology. Altogether, our results illustrate how climatic cycles, recolonization and selection can establish strong local adaptation along contact zones and affect the genetic architecture of adaptive traits.

Modeling the effect of adaptation to future climate change on spring phenological trend of European beech (Fagus sylvatica L.)

Temperate trees could cope with climate change through phenotypic plasticity of phenological key events or adaptation in situ via selection on genetic variation. However, the relative contribution of local adaptation and phenotypic plasticity to phenological change is unclear for many ecologically important tree species. Here, we analyzed the leaf-out data of European beech (Fagus sylvatica L.) from 50 provenances planted in 7 trial sites. We first constructed a function between chilling accumulation (CA) and photoperiod-associated heat requirement (PHR) of leaf-out date for each provenance and quantified the relationship between parameters of the CA-PHR function and climatic variables at provenance origins by using the random forest model. Furthermore, we used the provenance-specific CA-PHR function to simulate future leaf-out dates under two climate change scenarios (RCP 4.5 and 8.5) and two assumptions (no adaptation and adaptation). The results showed that both CA, provenance, and their interactions affected the PHR of leaf-out. The provenances from southeastern Europe exhibited a stronger response of PHR to CA and thus flushed earlier than northwestern provenances. The parameters of the CA-PHR function were connected with climatic variables (e.g., mean diurnal temperature range, temperature seasonality) at the originating sites of each provenance. If only considering the phenotypic plasticity, the leaf-out date of European beech in 2070-2099 will advance by 6.8 and 9.0 days on average relative to 1951-2020 under RCP 4.5 and RCP 8.5, respectively. However, if F. sylvatica adapts to future climate change by adopting the current strategy, the advance of the leaf-out date will weaken by 1.4 and 3.4 days under RCP 4.5 and RCP 8.5, respectively. Our results suggest that the European beech could slow down its spring phenological advances and reduce its spring frost risk if it adopts the current strategy to adapt to future climate change.

Linking leaf economic and hydraulic traits with early-age growth performance and survival of Eucalyptus pauciflora

Selection on plant functional traits may occur through their direct effects on fitness (or a fitness component), or may be mediated by attributes of plant performance which have a direct impact on fitness. Understanding this link is particularly challenging for long-lived organisms, such as forest trees, where lifetime fitness assessments are rarely achievable, and performance features and fitness components are usually quantified from early-life history stages. Accordingly, we studied a cohort of trees from multiple populations of Eucalyptus pauciflora grown in a common-garden field trial established at the hot and dry end of the species distribution on the island of Tasmania, Australia. We related the within-population variation in leaf economic (leaf thickness, leaf area and leaf density) and hydraulic (stomatal density, stomatal length and vein density) traits, measured from two-year-old plants, to two-year growth performance (height and stem diameter) and to a fitness component (seven-year survival). When performance-trait relationships were modelled for all traits simultaneously, statistical support for direct effects on growth performance was only observed for leaf thickness and leaf density. Performance-based estimators of directional selection indicated that individuals with reduced leaf thickness and increased leaf density were favoured. Survival-performance relationships were consistent with size-dependent mortality, with fitness-based selection gradients estimated for performance measures providing evidence for directional selection favouring individuals with faster growth. There was no statistical support for an effect associated with the fitness-based quadratic selection gradient estimated for growth performance. Conditional on a performance measure, fitness-based directional selection gradients estimated for the leaf traits did not provide statistical support for direct effects of the focal traits on tree survival. This suggested that, under the environmental conditions of the trial site and time period covered in the current study, early-stage selection on the studied leaf traits may be mediated by their effects on growth performance, which in turn has a positive direct influence on later-age survival. We discuss the potential mechanistic basis of the direct effects of the focal leaf traits on tree growth, and the relevance of a putative causal pathway of trait effects on fitness through mediation by growth performance in the studied hot and dry environment.

Unveiling a hotspot of genetic diversity in southern Italy for the endangered Hermann’s tortoise Testudo hermanni

Background

Hotspots of intraspecific genetic diversity represent invaluable resources for species to cope with environmental changes, and their identification is increasingly recognized as a major goal of conservation ecology research. However, even for iconic and endangered species, conservation strategies are often planned without thorough information on the geographic patterns of genetic variation. Here, we investigated the spatial patterns of genetic variation of the endangered Hermann’s tortoise Testudo hermanni in the Italian Peninsula by genotyping 174 individuals at 7 microsatellite loci, with the aim to contribute to planning effective conservation strategies.

Results

Ordination-based and Bayesian clustering analyses consistently identified three main genetic clusters, one spread in the central and northern part of the peninsula, and two restricted to southern Italy and Sicily, respectively. The highest levels of genetic diversity were found in populations of the southern cluster and, in particular, at the northern edges of its distribution (He > 0.6, Ar > 2.8 ), that correspond to areas of putative secondary contact and admixture between distinct lineages. Our results clearly identify a hotspot of genetic diversity for the Hermann’s tortoise in southern Italy.

Conclusion

We inferred the evolutionary history and the spatial patterns of genetic variation of the Hermann’s tortoise in the Italian Peninsula. We identified three main genetic clusters along the peninsula and a hotspot of intraspecific diversity in southern Italy. Our results underline the urgent need for conservation actions to warrant the long-term persistence of viable tortoise populations in this area. Furthrmore, these data add further evidence to the role of southern Italy as a biodiversity hotspot for temperate fauna, claiming for higher consideration of this area in large scale conservation programs.

Dynamics of internal isoprenoid metabolites in young Picea abies (Norway spruce) shoots during drought stress conditions in springtime

Currently, large areas of Picea abies (Norway spruce) stands in Europe are increasingly affected by drought and heat waves. Moreover, early spring drought has occurred with much higher frequency. Our work focuses on physiological changes induced by drought in four-year-old spruce seedlings during shoot elongation. We investigated drought effect on photosynthetic rate, concentration of abscisic acid and its metabolites, amount and composition of monoterpenes in needles of seedlings from five different provenances (altitude range 550-1280 m above sea level) in Western Carpathians. Spruce seedlings subjected to one-month drought stress of moderate intensity (about 50% of soil water content at the end of experiment) showed significant reduction of CO₂ uptake and increased concentration of ABA related to untreated controls. Induced drought affected needle monoterpene content and composition. Observed changes in drought-induced physiological parameters were influenced by seedling provenance. The provenance from 920 m above sea level showed the greatest sensitivity to drought with significantly highest ABA content and, at the same time, a clear decline of CO₂ uptake and amounts of total monoterpenes. Our results indicating intra-specific provenance-related variability in physiological response of spruce seedlings to drought may provide a basis for improved reforestation strategies in drought risk areas.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint – the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.