Drought Sensitivity of Norway Spruce at the Species' Warmest Fringe: Quantitative and Molecular Analysis Reveals High Genetic Variation Among and Within Provenances

Circadian rhythm and redox homeostasis candidate genes showed association with shallow elevation in Norway spruce

The analysis of genetic variation underlying local adaptation in natural populations, together with the response to different external stimuli, is currently a hot topic in forest sciences, with the aim of identifying genetic markers controlling key phenotypic traits of interest for their inclusion in restoration and breeding programs. In Europe, one of the main tree species is Norway spruce (Picea abies (L.) H.Karst.). Using the MassARRAY® platform, 568 trees from North Rhine-Westphalia (Germany) were genotyped with 94 single nucleotide polymorphisms (SNPs) related to circadian and growth rhythms, and to stress response. The association analysis of the selected markers with health status and elevation was performed using three different methods, and those identified by at least two of these were considered as high confidence associated SNPs. While just five markers showed a weak association with health condition, 32 SNPs were correlated with elevation, six of which were considered as high confidence associated SNPs, as indicated by at least two different association methods. Among these genes, thioredoxin and pseudo response regulator 1 (PRR1) are involved in redox homeostasis and ROS detoxification, APETALA2-like 3 (AP2L3), a transcription factor, is involved in seasonal apical growth, and a RPS2-like is a disease resistance gene. The function of some of these genes in controlling light-dependent reactions and metabolic processes suggests signatures of adaptation to local photoperiod and the synchronization of the circadian rhythm. This work provides new insights into the genetic basis of local adaptation over a shallow elevation gradient in Norway spruce.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Growth, drought response, and climate-associated genomic structure in whitebark pine in the Sierra Nevada of California

Whitebark pine (Pinus albicaulis Engelm.) has experienced rapid population declines and is listed as threatened under the Endangered Species Act in the United States. Whitebark pine in the Sierra Nevada of California represents the southernmost end of the species' distribution and, like other portions of its range, faces threats from an introduced pathogen, native bark beetles, and a rapidly warming climate. Beyond these chronic stressors, there is also concern about how this species will respond to acute stressors, such as drought. We present patterns of stem growth from 766 large (average diameter at breast height >25 cm), disease-free whitebark pine across the Sierra Nevada before and during a recent period of drought. We contextualize growth patterns using population genomic diversity and structure from a subset of 327 trees. Sampled whitebark pine generally had positive to neutral stem growth trends from 1970 to 2011, which was positively correlated with minimum temperature and precipitation. Indices of stem growth during drought years (2012 to 2015) relative to a predrought interval were mostly positive to neutral at our sampled sites. Individual tree growth response phenotypes appeared to be linked to genotypic variation in climate-associated loci, suggesting that some genotypes can take better advantage of local climatic conditions than others. We speculate that reduced snowpack during the 2012 to 2015 drought years may have lengthened the growing season while retaining sufficient moisture to maintain growth at most study sites. Growth responses may differ under future warming, however, particularly if drought severity increases and modifies interactions with pests and pathogens.

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.

Genetic basis of growth reaction to drought stress differs in contrasting high-latitude treeline ecotones of a widespread conifer

Climate change is increasing the frequency and intensity of drought events in many boreal forests. Trees are sessile organisms with a long generation time, which makes them vulnerable to fast climate change and hinders fast adaptations. Therefore, it is important to know how forests cope with drought stress and to explore the genetic basis of these reactions. We investigated three natural populations of white spruce (Picea glauca) in Alaska, located at one drought-limited and two cold-limited treelines with a paired plot design of one forest and one treeline plot. We obtained individual increment cores from 458 trees and climate data to assess dendrophenotypes, in particular the growth reaction to drought stress. To explore the genetic basis of these dendrophenotypes, we genotyped the individual trees at 3000 single nucleotide polymorphisms in candidate genes and performed genotype-phenotype association analysis using linear mixed models and Bayesian sparse linear mixed models. Growth reaction to drought stress differed in contrasting treeline populations. Therefore, the populations are likely to be unevenly affected by climate change. We identified 40 genes associated with dendrophenotypic traits that differed among the treeline populations. Most genes were identified in the drought-limited site, indicating comparatively strong selection pressure of drought-tolerant phenotypes. Contrasting patterns of drought-associated genes among sampled sites and in comparison to Canadian populations in a previous study suggest that drought adaptation acts on a local scale. Our results highlight genes that are associated with wood traits which in turn are critical for the establishment and persistence of future forests under climate change.

Revealing the Complex Relationship Among Hyperspectral Reflectance, Photosynthetic Pigments, and Growth in Norway Spruce Ecotypes

Norway spruce has a wide natural distribution range, harboring substantial physiological and genetic variation. There are three altitudinal ecotypes described in this species. Each ecotype has been shaped by natural selection and retains morphological and physiological characteristics. Foliar spectral reflectance is readily used in evaluating the physiological status of crops and forest ecosystems. However, underlying genetics of foliar spectral reflectance and pigment content in forest trees has rarely been investigated. We assessed the reflectance in a clonal bank comprising three ecotypes in two dates covering different vegetation season conditions. Significant seasonal differences in spectral reflectance among Norway spruce ecotypes were manifested in a wide-ranging reflectance spectrum. We estimated significant heritable variation and uncovered phenotypic and genetic correlations among growth and physiological traits through bivariate linear models utilizing spatial corrections. We confirmed the relative importance of the red edge within the context of the study site's ecotypic variation. When interpreting these findings, growth traits such as height, diameter, crown length, and crown height allowed us to estimate variable correlations across the reflectance spectrum, peaking in most cases in wavelengths connected to water content in plant tissues. Finally, significant differences among ecotypes in reflectance and other correlated traits were detected.

Integrating genomic information and productivity and climate-adaptability traits into a regional white spruce breeding program

Tree improvement programs often focus on improving productivity-related traits; however, under present climate change scenarios, climate change-related (adaptive) traits should also be incorporated into such programs. Therefore, quantifying the genetic variation and correlations among productivity and adaptability traits, and the importance of genotype by environment interactions, including defense compounds involved in biotic and abiotic resistance, is essential for selecting parents for the production of resilient and sustainable forests. Here, we estimated quantitative genetic parameters for 15 growth, wood quality, drought resilience, and monoterpene traits for Picea glauca (Moench) Voss (white spruce). We sampled 1,540 trees from three open-pollinated progeny trials, genotyped with 467,224 SNP markers using genotyping-by-sequencing (GBS). We used the pedigree and SNP information to calculate, respectively, the average numerator and genomic relationship matrices, and univariate and multivariate individual-tree models to obtain estimates of (co)variance components. With few site-specific exceptions, all traits examined were under genetic control. Overall, higher heritability estimates were derived from the genomic- than their counterpart pedigree-based relationship matrix. Selection for height, generally, improved diameter and water use efficiency, but decreased wood density, microfibril angle, and drought resistance. Genome-based correlations between traits reaffirmed the pedigree-based correlations for most trait pairs. High and positive genetic correlations between sites were observed (average 0.68), except for those pairs involving the highest elevation, warmer, and moister site, specifically for growth and microfibril angle. These results illustrate the advantage of using genomic information jointly with productivity and adaptability traits, and defense compounds to enhance tree breeding selection for changing climate.

Genome-wide association identifies candidate genes for drought tolerance in coast redwood and giant sequoia

Drought is a major limitation for survival and growth in plants. With more frequent and severe drought episodes occurring due to climate change, it is imperative to understand the genomic and physiological basis of drought tolerance to be able to predict how species will respond in the future. In this study, univariate and multitrait multivariate genome-wide association study methods were used to identify candidate genes in two iconic and ecosystem-dominating species of the western USA, coast redwood and giant sequoia, using 10 drought-related physiological and anatomical traits and genome-wide sequence-capture single nucleotide polymorphisms. Population-level phenotypic variation was found in carbon isotope discrimination, osmotic pressure at full turgor, xylem hydraulic diameter, and total area of transporting fibers in both species. Our study identified new 78 new marker × trait associations in coast redwood and six in giant sequoia, with genes involved in a range of metabolic, stress, and signaling pathways, among other functions. This study contributes to a better understanding of the genomic basis of drought tolerance in long-generation conifers and helps guide current and future conservation efforts in the species.

Genetic diversity of Norway spruce ecotypes assessed by GBS-derived SNPs

We investigated the genetic structure of three phenotypically distinct ecotypic groups of Norway spruce (Picea abies) belonging to three elevational classes; namely, low- (acuminata), medium- (europaea), and high-elevation (obovata) form, each represented by 150 trees. After rigorous filtering, we used 1916 Genotyping-by-Sequencing generated SNPs for analysis. Outputs from three multivariate analysis methods (Bayesian clustering algorithm implemented in STRUCTURE, Principal Component Analysis, and the Discriminant Analysis of Principal Components) indicated the presence of a distinct genetic cluster representing the high-elevation ecotypic group. Our findings bring a vital message to forestry practice affirming that artificial transfer of forest reproductive material, especially for stands under harsh climate conditions, should be considered with caution.

Low resistance but high resilience to drought of flushing Norway spruce seedlings

Spring drought episodes are becoming more frequent and intensive in European temperate forests. To study tree resilience to spring drought, Norway spruce seedlings were exposed to three levels of drought stress (well-watered (W), moderately stressed (M) and severely stressed (S)) for 42 days and then fully irrigated for 14 days. Drought strongly reduced gas exchange parameters for both M and S seedlings. After 42 days, stomatal conductance was lower by 83 and 97% in M and S, respectively, than in W seedlings. Respiration prevailed over photosynthesis in S seedlings at the end of the drought period. Drought mostly reduced longitudinal growth, especially in shoots and needles. Xylem growth reduction was caused mainly by a lower number of newly produced tracheids, not by changes in their size. Norway spruce seedlings showed good resilience to spring drought, as the observed physiological parameters started to recover after rewatering and seedlings started to sprout and form new tracheids. In M seedlings, all physiological traits recovered to the level of W seedlings during the 14-day irrigation period but the recovery took longer in S seedlings. Shoots and needles did not regrow in length but leaf mass per area increased during the recovery phase. To conclude, Norway spruce seedlings showed good resilience to spring single-drought event, but time necessary to full recovery from stress could make seedlings more vulnerable to recurrent drought events.

Post-glacial re-colonization and natural selection have shaped growth responses of silver fir across Europe

Warmer climate and more frequent extreme droughts will pose major threats to forest ecosystems. Past demography processes due to post-glacial recolonization and adaptation to local environmental conditions are among the main contributors to genetic differentiation processes among provenances. Assessing the intra-specific variability of tree growth responses to such changes is crucial to explore a species' potential to cope with climate warming. We combined growth-related traits derived from tree-ring width series with neutral genetic information of 18 European provenances of silver fir (Abies alba Mill.) growing in two common garden experiments in Switzerland. Analyses based on neutral single nucleotide polymorphisms revealed that the studied provenances grouped into three longitudinal clusters. These three genetic clusters showed differences in growth traits (height and DBH), with the provenances from the eastern cluster exhibiting the highest growth. The Pyrenees cluster showed significantly lower recovery and resilience to the extreme drought of 2003 as well as lower values of growth autocorrelation. Q_ST-F_ST and correlation analyses with climate of provenance origin suggest that the differences among provenances found in some traits result from natural selection. Our study suggests that the last post-glacial re-colonization and natural selection are the major drivers explaining the intra-specific variability in growth of silver fir across Europe. These findings highlight the importance of combining dendroecology and genetic analyses on fitness-related traits to assess the potential of a species to cope with global environmental change and provide insights to support assisted gene flow to ensure the persistence of the species in European forests.

Candidate regulators and target genes of drought stress in needles and roots of Norway spruce

Drought stress impacts seedling establishment, survival and whole-plant productivity. Molecular responses to drought stress have been most extensively studied in herbaceous species, mostly considering only aboveground tissues. Coniferous tree species dominate boreal forests, which are predicted to be exposed to more frequent and acute drought as a result of ongoing climate change. The associated impact at all stages of the forest tree life cycle is expected to have large-scale ecological and economic impacts. However, the molecular response to drought has not been comprehensively profiled for coniferous species. We assayed the physiological and transcriptional response of Picea abies (L.) H. Karst seedling needles and roots after exposure to mild and severe drought. Shoots and needles showed an extensive reversible plasticity for physiological measures indicative of drought-response mechanisms, including changes in stomatal conductance (gs), shoot water potential and abscisic acid (ABA). In both tissues, the most commonly observed expression profiles in response to drought were highly correlated with the ABA levels. Still, root and needle transcriptional responses contrasted, with extensive root-specific down-regulation of growth. Comparison between previously characterized Arabidopsis thaliana L. drought-response genes and P. abies revealed both conservation and divergence of transcriptional response to drought. In P. abies, transcription factors belonging to the ABA responsive element(ABRE) binding/ABRE binding factors ABA-dependent pathway had a more limited role. These results highlight the importance of profiling both above- and belowground tissues, and provide a comprehensive framework to advance the understanding of the drought response of P. abies. The results demonstrate that a short-term, severe drought induces severe physiological responses coupled to extensive transcriptome modulation and highlight the susceptibility of Norway spruce seedlings to such drought events.

Looking for the needle in a downsized haystack: Whole-exome sequencing unravels genomic signals of climatic adaptation in Douglas-fir (Pseudotsuga menziesii)

Conifers often occur along steep gradients of diverse climates throughout their natural ranges, which is expected to result in spatially varying selection to local climate conditions. However, signals of climatic adaptation can often be confounded, because unraveled clines covary with signals caused by neutral evolutionary processes such as gene flow and genetic drift. Consequently, our understanding of how selection and gene flow have shaped phenotypic and genotypic differentiation in trees is still limited.A 40-year-old common garden experiment comprising 16 Douglas-fir (Pseudotsuga menziesii) provenances from a north-to-south gradient of approx. 1,000 km was analyzed, and genomic information was obtained from exome capture, which resulted in an initial genomic dataset of >90,000 single nucleotide polymorphisms. We used a restrictive and conservative filtering approach, which permitted us to include only SNPs and individuals in environmental association analysis (EAA) that were free of potentially confounding effects (LD, relatedness among trees, heterozygosity deficiency, and deviations from Hardy-Weinberg proportions). We used four conceptually different genome scan methods based on FST outlier detection and gene-environment association in order to disentangle truly adaptive SNPs from neutral SNPs.We found that a relatively small proportion of the exome showed a truly adaptive signal (0.01%-0.17%) when population substructuring and multiple testing was accounted for. Nevertheless, the unraveled SNP candidates showed significant relationships with climate at provenance origins, which strongly suggests that they have featured adaptation in Douglas-fir along a climatic gradient. Two SNPs were independently found by three of the employed algorithms, and one of them is in close proximity to an annotated gene involved in circadian clock control and photoperiodism as was similarly found in Populus balsamifera. Synthesis. We conclude that despite neutral evolutionary processes, phenotypic and genomic signals of adaptation to climate are responsible for differentiation, which in particular explain disparity between the well-known coastal and interior varieties of Douglas-fir.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Piceaglauca [Moench] Voss), by combining genotype–environment associations, genotype–phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought‐related traits (i.e., anatomical, growth response and climate‐sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse‐controlled conditions. The interlinked phenotype–genotype–environment network revealed eight high‐confidence genes involved in white spruce adaptation to drought, of which four were drought‐responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

Evolvability of Drought Response in Four Native and Non-native Conifers: Opportunities for Forest and Genetic Resource Management in Europe

Intraspecific genetic variation in drought response is expected to play an important role in determining the persistence of tree populations in global change as it (1) allows for spontaneous selection and local adaptation of tree populations, (2) supports assisted seed transfer of less-drought-sensitive provenance, and (3) enables the integration of drought-sensitivity traits into tree breeding. Estimating the potential of such adaptation options requires quantitative genetic knowledge of drought sensitivity across significant parts of species distributions and a comparative assessment of genetic variation within economically and ecologically important tree species. We quantified genetic variation within and among populations of four conifers growing within common garden experiments in the drought-prone eastern Austria. This region experienced three strong drought periods between 1980 and 2010 that resulted in significant reductions in radial growth. Among the four tested species, Douglas-fir revealed the highest resistance during drought and silver fir the best recovery after drought, while European larch and Norway spruce showed the lowest resistance. High genetic variation among populations and phenotypic stability across all three drought events was found for Norway spruce and silver fir, but not for the other species. Heritability and evolvability of drought traits, both approximated via genetic repeatability, revealed strong differences among populations of all four species. Repeatability and evolvability for resistance were highest in Norway spruce and, for recovery, highest in European larch. Our comparison indicates that the mean drought sensitivity of a species is not related to the intraspecific genetic variation in drought response. Thus, also highly drought-sensitive species, such as Norway spruce and European larch, harbor significant genetic variation in drought response within and among populations to justify targeted tree breeding, assisted gene flow, and supportive forest management to foster local adaptations to future conditions.

Variability among Sites and Climate Models Contribute to Uncertain Spruce Growth Projections in Denmark

Projecting trees species growth into future climate is subject to large uncertainties and it is of importance to quantify the different sources (e.g., site, climate model) to prioritize research efforts. This study quantifies and compares sites and climate model-induced uncertainties in projected Norway spruce growth from Denmark. We analyzed tree-rings from 340 Norway spruce trees sampled in 14 planted stands (1. Plantation; period 1950–1987) and additionally 36 trees from six trials in a common garden experiment (2. Common garden; period 1972–2012). Growth-climate correlations were estimated and multiple linear and nonlinear regression models relating growth with climate were tested. Tree growth was projected up to 2100 applying multiple linear or quadratic regression models based on the 15 Atmosphere-Ocean General Circulation Models (AOGCMs) of the Coupled Model Inter-comparison Project Phase 5 (CMIP5). The climate-growth models showed that summer drought and warm previous-year late-summer and early-autumn constrain growth. In some stands, warm springs affected growth positively. The projections of growth under future climates on average showed from no to slightly negative changes in growth compared to present growth rates. However, projections showed a very large variation, ranging from highly positive to highly negative growth changes. The uncertainties due to variation in site responses and in climate models were substantial. A lesser degree of uncertainty was related to the emission scenarios. Even though our projections on average suggest that Norway spruce may experience a growth reduction in the future, the tremendous variation in growth predictions due to differences between stands and climate models calls for further research and caution when projections are interpreted. These results also suggest that forest managers in general should avoid the use of Norway spruce on exposed and drought prone sites and as an additional resilience measure primarily use it in mixtures with other more climate tolerant species.

The dehydrins gene expression differs across ecotypes in Norway spruce and relates to weather fluctuations

Norway spruce has a broad natural distribution range, which results in a substantial variety of its physiological and genetic variation. There are three distinct altitudinal ecotypes described in this tree species. The physiological optimum of each ecotype may be shifted due to ongoing climate change, especially in traits associated with water demand that might be crucial for adaptation. Dehydrins are proteins that help to mitigate the adverse effects of dehydration. Dehydrin gene expression patterns appeared to be a suitable marker for plant stress assessment. Genetically determined differences in response between individuals and populations were formerly studied, however, mainly in controlled conditions. We evaluated ecotypic variation in dehydrin gene expression in a clonal bank comprised of all three ecotypes. A genetic relationship among targeted trees was uncovered utilizing GBS (Genotyping by Sequencing) platform. We sampled 4-6 trees of each ecotype throughout 15 months period. Subsequently, we assessed the RNA expression of dehydrin genes by qRT-PCR. For this study, we deliberately selected dehydrins from different categories. Our findings detected significant differences among ecotypes in dehydrin expression. The association of recorded climatic variables and individual gene expression across the study period was evaluated and revealed, for certain genes, a correlation between dehydrin gene expression and precipitation, temperature, and day-length.

Shaping Durum Wheat for the Future: Gene Expression Analyses and Metabolites Profiling Support the Contribution of BCAT Genes to Drought Stress Response

Global climate change, its implications for agriculture, and the complex scenario presented by the scientific community are of worldwide concern. Drought is a major abiotic stress that can restrict plants growth and yields, thus the identification of genotypes with higher adaptability to drought stress represents one of the primary goals in breeding programs. During abiotic stress, metabolic adaptation is crucial for stress tolerance, and accumulation of specific amino acids and/or as secondary metabolites deriving from amino acid metabolism may correlate with the increased tolerance to adverse environmental conditions. This work, focused on the metabolism of branched chain-amino acids (BCAAs) in durum wheat and the role of branched-chain amino acid aminotransferases (BCATs) in stress response. The role of BCATs in plant response to drought was previously proposed for Arabidopsis, where the levels of BCAAs were altered at the transcriptional level under drought conditions, triggering the onset of defense response metabolism. However, in wheat the role of BCAAs as a trigger of the onset of the drought defense response has not been elucidated. A comparative genomic approach elucidated the composition of the BCAT gene family in durum wheat. Here we demonstrate a tissue and developmental stage specificity of BCATs regulation in the drought response. Moreover, a metabolites profiling was performed on two contrasting durum wheat cultivars Colosseo and Cappelli resulting in the detection of a specific pattern of metabolites accumulated among genotypes and, in particular, in an enhanced BCAAs accumulation in the tolerant cv Cappelli further supporting a role of BCAAs in the drought defense response. The results support the use of gene expression and target metabolomic in modern breeding to shape new cultivars more resilient to a changing climate.

Adaptive genetic variation to drought in a widely distributed conifer suggests a potential for increasing forest resilience in a drying climate

Drought intensity and frequency are increasing under global warming, with soil water availability now being a major factor limiting tree growth in circumboreal forests. Still, the adaptive capacity of trees in the face of future climatic regimes remains poorly documented. Using 1481 annually resolved tree-ring series from 29-yr-old trees, we evaluated the drought sensitivity of 43 white spruce (Picea glauca (Moench) Voss) populations established in a common garden experiment. We show that genetic variation among populations in response to drought plays a significant role in growth resilience. Local genetic adaptation allowed populations from drier geographical origins to grow better, as indicated by higher resilience to extreme drought events, compared with populations from more humid geographical origins. The substantial genetic variation found for growth resilience highlights the possibility of selecting for drought resilience in boreal conifers. As a major research outcome, we showed that adaptive genetic variation in response to changing local conditions can shape drought vulnerability at the intraspecific level. Our findings have wide implications for forest ecosystem management and conservation.

Individual reproductive success in Norway spruce natural populations depends on growth rate, age and sensitivity to temperature

Quantifying the individual reproductive success and understanding its determinants is a central issue in evolutionary research for the major consequences that the transmission of genetic variation from parents to offspring has on the adaptive potential of populations. Here, we propose to distil the myriad of information embedded in tree-ring time series into a set of tree-ring-based phenotypic traits to be investigated as potential drivers of reproductive success in forest trees. By using a cross-disciplinary approach that combines parentage analysis and a thorough dendrophenotypic characterisation of putative parents, we assessed sex-specific relationships between such dendrophenotypic traits (i.e., age, growth rate and parameters describing sensitivity to climate and to extreme climatic events) and reproductive success in Norway spruce. We applied a full probability method for reconstructing parent-offspring relationships between 604 seedlings and 518 adult trees sampled within five populations from southern and central Europe. We found that individual female and male reproductive success was positively associated with tree growth rate and age. Female reproductive success was also positively influenced by the correlation between growth and the mean temperature of the previous vegetative season. Overall, our results showed that Norway spruce individuals with the highest fitness are those who are able to keep high-growth rates despite potential growth limitations caused by reproductive costs and climatic limiting conditions. Identifying such functional links between the individual ecophysiological behaviour and its evolutionary gain would increase our understanding on how natural selection shapes the genetic composition of forest tree populations over time.

Climate sensitivity of Cryptomeria japonica in two contrasting environments: Perspectives from QTL mapping

Long-lived forest tree species experience a wide range of environmental conditions throughout their lifespan. Evaluation of the underlying growth and development mechanisms of these species is essential to predict tree growth under climate change. This study investigated climate sensitivity to temperature, precipitation, dry periods, and the associated genomic regions in Cryptomeria japonica, Japan’s most commercially important tree. We used tree rings and common garden experiments with three clonal replicates planted in two contrasting environments in Kyushu (Kumamoto site) and Honshu (Chiba site), Japan. Tree growth showed a significant negative correlation with the dry period (>4 days) in March of the year of tree-ring formation at the Chiba site. In contrast, temperature and precipitation had little influence on tree growth. Quantitative trait locus (QTL) analysis was performed to investigate climate sensitivity to dry periods at the Chiba site, revealing 13 significant QTLs. One QTL showed a substantially large contribution to the overall climate sensitivity, accounting for 12.4% of the total phenotypic variation. The phenotypic variance explained (PVE) by other QTLs ranged from 0.9% to 2.9%, and the total PVE by all QTLs was 35.6%. These findings indicate that the tree population at the Chiba site could be vulnerable to drought in early spring and that the QTL showing the greatest impact on climate sensitivity may be closely related to genes associated with tolerance or adaptation to drought stress. The QTLs identified in this study could be useful for molecular breeding, forest management, and predicting the growth of C. japonica under a changing climate.

The GenTree Dendroecological Collection, tree-ring and wood density data from seven tree species across Europe

The dataset presented here was collected by the GenTree project (EU-Horizon 2020), which aims to improve the use of forest genetic resources across Europe by better understanding how trees adapt to their local environment. This dataset of individual tree-core characteristics including ring-width series and whole-core wood density was collected for seven ecologically and economically important European tree species: silver birch (Betula pendula), European beech (Fagus sylvatica), Norway spruce (Picea abies), European black poplar (Populus nigra), maritime pine (Pinus pinaster), Scots pine (Pinus sylvestris), and sessile oak (Quercus petraea). Tree-ring width measurements were obtained from 3600 trees in 142 populations and whole-core wood density was measured for 3098 trees in 125 populations. This dataset covers most of the geographical and climatic range occupied by the selected species. The potential use of it will be highly valuable for assessing ecological and evolutionary responses to environmental conditions as well as for model development and parameterization, to predict adaptability under climate change scenarios.

Adaptive limitations of white spruce populations to drought imply vulnerability to climate change in its western range

A cost-effective climate change adaptation strategy for the forestry sector is to move seed sources to more northern and higher elevation planting sites as part of ongoing reforestation programs. This is meant to match locally adapted populations with anticipated environments, but adaptive traits do not always show population differences suitable to mitigate climate change impacts. For white spruce, drought tolerance is a critical adaptive trait to prevent mortality and productivity losses. Here, we use a 40-year-old provenance experiment that has been exposed to severe drought periods in 1999 and 2002 to retrospectively investigate drought response and the adaptive capacity of white spruce populations across their boreal range. Relying on dendrochronological analysis under experimentally controlled environments, we evaluate population differences in resistance, resilience, and recovery to these extreme events. Results showed evidence for population differentiation in resistance and recovery parameters, but provenances conformed to approximately the same growth rates under drought conditions and had similar resilience metrics. The lack of populations with better growth rates under drought conditions is contrary to expectations for a wide-ranging species with distinct regional climates. Populations from the wettest environments in the northeastern boreal were surprisingly drought-tolerant, suggesting that these populations would readily resist water deficits projected for the 2080s, and supporting the view that northeastern Canada will provide a refugium for boreal species under climate change. The findings also suggest that white spruce is sensitive to growth reductions under climate change in the western boreal. The study highlights that population differentiation in adaptive capacity is species- and trait-specific, and we provide a counterexample for drought tolerance traits, where assisted migration prescriptions may be ineffective to mitigate climate change impacts. For resource managers and policy makers, we provide maps where planning for widespread declines of boreal white spruce forests may be unavoidable.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Disentangling the role of climate and soil on tree growth and its interaction with seed origin

When considering options for adapting forests under climate change, climate is treated as the dominant driver of forest growth, while soil properties are often ignored mainly due to shortage of accurate data. The effects of climate and soil on forest growth may vary due to local adaptation to both climate and soil, and these local adaptations might need to be considered when transferring seed provenances under climate change. Data from 29 provenance trials of Norway spruce (Picea abies (L.) Karst.) across a wide gradient of planting conditions in Austria was used to develop Structural Equation Models (SEMs) to quantified the role of climatic and soil drivers and their interactions on juvenile growth performance and to test if provenance origin affects the relative importance of these drivers. Climate and soil of the planting site location were found to have similar direct effects on juvenile tree growth, however, climate was found to be more important because of additional indirect effects via interactions with soil parameters. Notably, the relative effects of climate and soil vary among different provenance groups. Climate constraints are dominant for seed sources originating from colder and/or high altitude locations, while test site climate and soil are equally important contributors of growth for provenances originating from warmer origin and lower elevation sites. Together with the better growth performance of the latter provenance group their plasticity allows them to utilize a wide range of soil conditions.

Long-Term Studies Reveal Differential Responses to Climate Change for Trees Under Soil- or Herbivore-Related Stress

Worldwide, trees are confronting increased temperature and aridity, exacerbating susceptibility to herbivory. Long-term studies comparing patterns of plant performance through drought can help identify variation among and within populations in vulnerability to climate change and herbivory. We use long-term monitoring data to examine our overarching hypothesis that the negative impacts of poor soil and herbivore susceptibility would be compounded by severe drought. We studied pinyon pine, Pinus edulis, a widespread southwestern tree species that has suffered extensive climate-change related mortality. We analyzed data on mortality, growth, male reproduction, and herbivory collected for 14-32 years in three areas with distinct soil-types. We used standardized precipitation-evapotranspiration index (SPEI) as a climate proxy that summarizes the impacts of drought due to precipitation and temperature variation on semi-arid forests. Several key findings emerged: (1) Plant performance measurements did not support our hypothesis that trees growing in stressful, coarse-textured soils would suffer more than trees growing in finer-textured soils. Stem growth at the area with coarse, young cinder soils (area one) responded only weakly to drought, while stem growth on more developed soils with sedimentary (area two) and volcanic (area three) substrates, was strongly negatively affected by drought. Male reproduction declined less with drought at area one and more at areas two and three. Overall mortality was 30% on coarse cinder soils (area one) and averaged 55% on finer soil types (areas two and three). (2) Although moth herbivore susceptible trees were hypothesized to suffer more with drought than moth resistant trees, the opposite occurred. Annual stem growth was negatively affected by drought for moth resistant trees, but much less strongly for moth susceptible trees. (3) In contrast to our hypothesis, moths declined with drought. Overall, chronically water-stressed and herbivore-susceptible trees had smaller declines in performance relative to less-stressed trees during drought years. These long-term findings support the idea that stressed trees might be more resistant to drought since they may have adapted or acclimated to resist drought-related mortality.

High Morphological Differentiation in Crown Architecture Contrasts with Low Population Genetic Structure of German Norway Spruce Stands

High elevation sites in the low mountain ranges in Germany are naturally covered by Norway spruce (Picea abies (Karst.) L.) stands. Historically, large scale anthropogenic range expansion starting in the mid to late 18th century had a huge impact on the forest composition throughout Germany. Utilisation and exploitation often led to artificial regeneration, mostly carried out using seeds from allochthonous provenances. Usually, autochthonous (natural) high elevation Norway spruce trees have narrow crown phenotypes, whereas lowland trees have broader crowns. Narrow crown phenotypes are likely the result of adaptation to heavy snow loads combined with high wind speeds. In the present study, neighbouring stand pairs of putative autochthonous and allochthonous origin with contrasting phenotypes in high elevation sites were investigated with 200 samples each. These stands are located in the Ore Mountains, the Thuringian Forest, and the Harz Mountains. Additionally, a relict population with the typical narrow high elevation phenotypes was sampled in Thuringia, known as “Schlossbergfichte”. The objective of the study was to quantify supposedly adaptive phenotypic differences in crown architecture and the genetic differentiation of 11 putatively neutral nuclear microsatellite markers (i.e., simple sequence repeats (nSSRs)). The high differentiation of morphological traits (PST = 0.952–0.989) between the neighbouring autochthonous and allochthonous stands of similar age contrasts with the very low neutral genetic differentiation (FST = 0.002–0.007; G″ST = 0.002–0.030), suggesting that directional selection at adaptive gene loci was involved in phenotypic differentiation. Comparing the regions, a small isolation by distance effect for the Harz Mountains was detected, suggesting landscape resistance restricting gene flow. Finally, the differentiation of the very old autochthonous (up to 250 years) stand “Schlossbergfichte” with typical high elevation phenotypes could cohere with the sampling of a relict genepool.

Influence of Spatiotemporal Dynamics on the Fine-Scale Spatial Genetic Structure of Differently Managed Picea abies Stands

The tandem analysis of dendrochronological and genetic data is piquing forest ecologists’ interest and represents a promising approach for studying the temporal development of genetic structure in forest tree populations. Such multidisciplinary approach can help elucidate to what extent different management practices have impacted the fine-scale spatial genetic structure of forest stands through time. In this study, we jointly analysed spatial, age and genetic data from three differently managed Norway spruce permanent plots to assess: (1) possible differences among plots in the spatial distribution of individuals and their genetic structure due to different management practices, and (2) whether modifications in the age structure influenced the fine-scale spatial genetic structure within each permanent plot. With these aims, we genetically characterized at five nuclear microsatellite markers a large subset (328) of all the trees for which spatial and age data were collected (1472). We found that different management practices determined a similar spatial structure in terms of trees’ ages (r < 25 m in all plots) and neutral genetic diversity (Sp ranging from 0.002 to 0.004). Hot spots and cold spots of trees’ age were not statistically different in terms of genetic diversity, and trees’ age was not statistically different among the genetic clusters detected. On the other hand, the spatial distribution of individuals was significantly clustered up to 22 m only in the wooded pasture plot. Our main findings show that forest land use and management can indeed determine markedly different spatial layouts of Norway spruce individuals but do not produce strong distortions in the spatial structure of age and genetic parameters.