Linking dendroecology and association genetics in natural populations: Stress responses archived in tree rings associate with SNP genotypes in silver fir (Abies albaMill.)

Diverging growth trends and climate sensitivities of individual pine trees after the 1976 extreme drought

Summer droughts are affecting the productivity and functioning of central European forests, with potentially lasting consequences for species composition and carbon sequestration. Long-term recovery rates and individual growth responses that may diverge from species-specific and population-wide behaviour are, however, poorly understood. Here, we present 2052 pine (Pinus sylvestris) ring width series from 19 forest sites in south-west Germany to investigate growth responses of individual trees to the exceptionally hot and dry summer of 1976. This outstanding drought event presents a distinctive test case to examine long-term post-drought recovery dynamics. We have proposed a new classification approach to identify a distinct sub-population of trees, referred to as "temporarily affected trees", with a prevalence ranging from 9 to 33 % across the forest stands. These trees exhibited an exceptionally prolonged growth suppression, lasting over a decade, indicating significantly lower resilience to the 1976 drought and a 50 % reduced capacity to recover to pre-drought states. Furthermore, shifts in resilience and recovery dynamics are accompanied by changing climate sensitivities, notably an increased response to maximum temperatures and summer droughts in post-1976 affected pines. Our findings underscore the likely interplay between individual factors and micro-site conditions that contribute to divergent tree responses to droughts. Assessing these factors at the individual tree level is recommended to advancing our understanding of forest responses to extreme drought events. By analyzing sub-population growth patterns, our study provides valuable insights into the impacts of summer droughts on central European forests in context of increasing drought events.

Genome-wide family prediction unveils molecular mechanisms underlying the regulation of agronomic traits in Urochloa ruziziensis

Tropical forage grasses, particularly those belonging to the Urochloa genus, play a crucial role in cattle production and serve as the main food source for animals in tropical and subtropical regions. The majority of these species are apomictic and tetraploid, highlighting the significance of U. ruziziensis, a sexual diploid species that can be tetraploidized for use in interspecific crosses with apomictic species. As a means to support breeding programs, our study investigates the feasibility of genome-wide family prediction in U. ruziziensis families to predict agronomic traits. Fifty half-sibling families were assessed for green matter yield, dry matter yield, regrowth capacity, leaf dry matter, and stem dry matter across different clippings established in contrasting seasons with varying available water capacity. Genotyping was performed using a genotyping-by-sequencing approach based on DNA samples from family pools. In addition to conventional genomic prediction methods, machine learning and feature selection algorithms were employed to reduce the necessary number of markers for prediction and enhance predictive accuracy across phenotypes. To explore the regulation of agronomic traits, our study evaluated the significance of selected markers for prediction using a tree-based approach, potentially linking these regions to quantitative trait loci (QTLs). In a multiomic approach, genes from the species transcriptome were mapped and correlated to those markers. A gene coexpression network was modeled with gene expression estimates from a diverse set of U. ruziziensis genotypes, enabling a comprehensive investigation of molecular mechanisms associated with these regions. The heritabilities of the evaluated traits ranged from 0.44 to 0.92. A total of 28,106 filtered SNPs were used to predict phenotypic measurements, achieving a mean predictive ability of 0.762. By employing feature selection techniques, we could reduce the dimensionality of SNP datasets, revealing potential genotype-phenotype associations. The functional annotation of genes near these markers revealed associations with auxin transport and biosynthesis of lignin, flavonol, and folic acid. Further exploration with the gene coexpression network uncovered associations with DNA metabolism, stress response, and circadian rhythm. These genes and regions represent important targets for expanding our understanding of the metabolic regulation of agronomic traits and offer valuable insights applicable to species breeding. Our work represents an innovative contribution to molecular breeding techniques for tropical forages, presenting a viable marker-assisted breeding approach and identifying target regions for future molecular studies on these agronomic traits.

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.

Effects of drought stress memory on the accumulation of stress-protective compounds in naturally grown pine and spruce

Forest trees are subjected to multiple stressors during their long lifetime and therefore require effective and finely regulated stress-protective systems. Stressors can induce protective systems either directly or with the involvement of stress memory mechanisms. Stress memory has only begun to be uncovered in model plants and is unexplored in coniferous species. Therefore, we studied the possible role of stress memory in the regulation of the accumulation of stress-protective compounds (heat shock proteins, dehydrins, proline) in the needles of naturally grown Scots pine and Norway spruce trees subjected to the subsequent action of long-term (multiyear) and short-term (seasonal) water shortages. Although the water deficit was relatively mild, it significantly influenced the pattern of expression of stress memory-related heat shock factor (HSF) and SWI/SNF genes, indicating the formation of stress memory in both species. In spruce, dehydrin accumulation was increased by water shortage in a manner compatible with Type II stress memory. The accumulation of HSP40 in spruce needles was positively influenced by long-term water shortage, but this increase was unlikely to be of biological importance due to the concomitant decrease in HSP70, HSP90 and HSP101 accumulation. Finally, proline accumulation was negatively influenced by short-term water deficit in spruce. In pine, no one protective compound accumulated in response to water stress. Taken together, the results indicate that the accumulation of stress-protective compounds was generally independent of stress memory effects both in pine and in spruce.

Fire frequency, as well as stress response and developmental gene control serotiny level variation in a widespread pioneer Mediterranean conifer, Pinus halepensis

Many plants undergo adaptation to fire. Yet, as global change is increasing fire frequency worldwide, our understanding of the genetics of adaptation to fire is still limited. We studied the genetic basis of serotiny (the ability to disseminate seeds exclusively after fire) in the widespread, pioneer Mediterranean conifer Pinus halepensis Mill., by linking individual variation in serotiny presence and level to fire frequency and to genetic polymorphism in natural populations. After filtering steps, 885 single nucleotide polymorphisms (SNPs) out of 8000 SNPs used for genotyping were implemented to perform an in situ association study between genotypes and serotiny presence and level. To identify serotiny‐associated loci, we performed random forest analyses of the effect of SNPs on serotiny levels, while controlling for tree size, frequency of wildfires, and background environmental parameters. Serotiny showed a bimodal distribution, with serotinous trees more frequent in populations exposed to fire in their recent history. Twenty‐two SNPs found in genes involved in stress tolerance were associated with the presence‐absence of serotiny while 37 found in genes controlling for flowering were associated with continuous serotiny variation. This study shows the high potential of P. halepensis to adapt to changing fire regimes, benefiting from a large and flexible genetic basis of trait variation.

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.

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.

The role of nutritional impairment in carbon-water balance of silver fir drought-induced dieback

Rear-edge populations at the xeric distribution limit of tree species are particularly vulnerable to forest dieback triggered by drought. This is the case of silver fir (Abies alba) forests located in Southwestern Europe. While silver fir drought-induced dieback patterns have been previously explored, information on the role played by nutritional impairment is lacking despite its potential interactions with tree carbon-water balances. We performed a comparative analysis of radial growth, intrinsic water-use efficiency (iWUE), oxygen isotopes (δ¹⁸ O) and nutrient concentrations in leaves of declining (DD) and non-declining (ND) trees in silver fir in four forests in the Spanish Pyrenees. We also evaluated the relationships among dieback predisposition, intraspecific trait variation (wood density and leaf traits) and rhizosphere soil physical-chemical properties. The onset of growth decline in DD trees occurred more than two decades ago, and they subsequently showed low growth resilience against droughts. The DD trees presented consistently lower foliar concentrations of nutrients such as P, K, Cu and Ni than ND trees. The strong effects of foliar nutrient status on growth resilience indices support the key role played by mineral nutrition in tree functioning and growth before, during and after drought. In contrast, variability in wood density and leaf morphological traits, as well as soil properties, showed weak relationships with tree nutritional status and drought performance. At the low elevation, warmer sites, DD trees showed stronger climate-growth relationships and lower δ¹⁸ O than ND trees. The uncoupling between iWUE and δ¹⁸ O, together with the positive correlations between P and K leaf concentrations and δ¹⁸ O, point to deeper soil/bedrock water sources and vertical decoupling between nutrient and water uptake in DD trees. This study provides novel insights into the mechanisms driving silver fir dieback and highlights the need to incorporate tree nutrition into forest dieback studies.

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.

Fine-scale spatial genetic structure across the species range reflects recent colonization of high elevation habitats in silver fir (Abies alba Mill.)

Variation in genetic diversity across species ranges has long been recognized as highly informative for assessing populations’ resilience and adaptive potential. The spatial distribution of genetic diversity within populations, referred to as fine‐scale spatial genetic structure (FSGS), also carries information about recent demographic changes, yet it has rarely been connected to range scale processes. We studied eight silver fir (Abies alba Mill.) population pairs (sites), growing at high and low elevations, representative of the main genetic lineages of the species. A total of 1,368 adult trees and 540 seedlings were genotyped using 137 and 116 single nucleotide polymorphisms (SNPs), respectively. Sites revealed a clear east‐west isolation‐by‐distance pattern consistent with the post‐glacial colonization history of the species. Genetic differentiation among sites (F_CT = 0.148) was an order of magnitude greater than between elevations within sites (F_SC = 0.031), nevertheless high elevation populations consistently exhibited a stronger FSGS. Structural equation modelling revealed that elevation and, to a lesser extent, post‐glacial colonization history, but not climatic and habitat variables, were the best predictors of FSGS across populations. These results suggest that high elevation habitats have been colonized more recently across the species range. Additionally, paternity analysis revealed a high reproductive skew among adults and a stronger FSGS in seedlings than in adults, suggesting that FSGS may conserve the signature of demographic changes for several generations. Our results emphasize that spatial patterns of genetic diversity within populations provide information about demographic history complementary to non‐spatial statistics, and could be used for genetic diversity monitoring, especially in forest trees.

Genome-wide approaches for the identification of markers and genes associated with sugarcane yellow leaf virus resistance

Sugarcane yellow leaf (SCYL), caused by the sugarcane yellow leaf virus (SCYLV) is a major disease affecting sugarcane, a leading sugar and energy crop. Despite damages caused by SCYLV, the genetic base of resistance to this virus remains largely unknown. Several methodologies have arisen to identify molecular markers associated with SCYLV resistance, which are crucial for marker-assisted selection and understanding response mechanisms to this virus. We investigated the genetic base of SCYLV resistance using dominant and codominant markers and genotypes of interest for sugarcane breeding. A sugarcane panel inoculated with SCYLV was analyzed for SCYL symptoms, and viral titer was estimated by RT-qPCR. This panel was genotyped with 662 dominant markers and 70,888 SNPs and indels with allele proportion information. We used polyploid-adapted genome-wide association analyses and machine-learning algorithms coupled with feature selection methods to establish marker-trait associations. While each approach identified unique marker sets associated with phenotypes, convergences were observed between them and demonstrated their complementarity. Lastly, we annotated these markers, identifying genes encoding emblematic participants in virus resistance mechanisms and previously unreported candidates involved in viral responses. Our approach could accelerate sugarcane breeding targeting SCYLV resistance and facilitate studies on biological processes leading to this trait.

Tackling the challenges of evolutionary forest research with multidata approaches

Many forest tree species have characteristics that make the study of their evolutionary ecology complex. For example, they are long-lived and thus have long generation times, and their often large, complex genomes have hampered establishing genomic resources. One way to tackle this challenge is to access multiple complementary data sources and analytical approaches when attempting to infer patterns of adaptive evolution. In the cover article of this issue of Molecular Ecology, Depardieu et al. (2021) combine large amounts of environmental, genomic, dendrochronological, and gene expression data in a common garden to explore the polygenic basis of drought resistance in white spruce (Picea glauca), a long-lived conifer. They identify candidate genes involved in growth and resistance to extreme drought events and show how multiple data sets may deliver complementary evidence to circumvent the manifold challenges of the research field.

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.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Intraspecific Growth Response to Drought of Abies alba in the Southeastern Carpathians

The intensity and frequency of drought have increased considerably during the last decades in southeastern Europe, and projected scenarios suggest that southern and central Europe will be affected by more drought events by the end of the 21st century. In this context, assessing the intraspecific genetic variation of forest tree species and identifying populations expected to be best adapted to future climate conditions is essential for increasing forest productivity and adaptability. Using a tree-ring database from 60 populations of 38-year-old silver fir (Abies alba) in five trial sites established across Romania, we studied the variation of growth and wood characteristics, provenance-specific response to drought, and climate-growth relationships during the period 1997–2018. The drought response of provenances was determined by four drought parameters: resistance, recovery, resilience, and relative resilience. Based on the standardized precipitation index, ten years with extreme and severe drought were identified for all trial sites. Considerable differences in radial growth, wood characteristics, and drought response parameters among silver fir provenances have been found. The provenances’ ranking by resistance, recovery, and resilience revealed that a number of provenances from Bulgaria, Italy, Romania, and Czech Republic placed in the top ranks in almost all sites. Additionally, there are provenances that combine high productivity and drought tolerance. The correlations between drought parameters and wood characters are positive, the most significant correlations being obtained between radial growth and resilience. Correlations between drought parameters and wood density were non-significant, indicating that wood density cannot be used as indicator of drought sensitivity. The negative correlations between radial growth and temperature during the growing season and the positive correlations with precipitation suggest that warming and water deficit could have a negative impact on silver fir growth in climatic marginal sites. Silvicultural practices and adaptive management should rely on selection and planting of forest reproductive material with high drought resilience in current and future reforestation programs.

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

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

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

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

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

The Pattern of Genetic Variation, Survival and Growth in the Abies alba Mill. Population within the Introgression Zone of Two Refugial Lineages in the Carpathians

Background and Objectives: The contact zones of different refugial lineages, where mixing of genetic backgrounds leads to new gene combinations or pre-adaptations, represent hotspots of genetic diversity. The aim of the study was to compare patterns in the genetic structure of the Abies alba Mill. population in the Eastern and Western Carpathians (Eastern Europe) within the introgression zone of two refugial lineages and the growth response of provenances located in a gradient of pollen-mediated gene fluxes. Materials and Methods: The mitochondrial nad5-4 marker and five polymorphic microsatellite nuclear markers (nSSR) were analyzed in 56 subpopulations from Romania, Ukraine, Slovakia and Poland. The survival rate and height growth up to an age of 15 years were compared for 33 subpopulations, forming a distance gradient between 170 and 470 km from the meeting zone of the refugial lineages. Results: The results of the analysis of molecular variance indicated that 8.2% of the total genetic variation is attributable to the between-subpopulation level and 1.7% to the between-lineage level. The pollen-mediated influence of the eastern lineage was detectable at a distance of at least 300 km in the western direction. Eastern provenances with origin sites closer to the meeting zone of the refugial lineages were characterized by lower survival rate and lower heights (about 8% lower than the average tree height) compared to subpopulations from the central and western part of the studied region. Conclusions: Pollen-mediated gene flow between lineages appears to have been sufficient to cause a significant change in phenotypic traits related to tree growth. Subpopulations from the central and western parts of the studied region are better adapted to current climatic conditions. Nonetheless, given the increasing aridity of the regional climate, a safe guideline is to increase genetic mixing.

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.

Reviewing the Use of Resilience Concepts in Forest Sciences

PURPOSE OF REVIEW

Resilience is a key concept to deal with an uncertain future in forestry. In recent years, it has received increasing attention from both research and practice. However, a common understanding of what resilience means in a forestry context and how to operationalise it is lacking. Here, we conducted a systematic review of the recent forest science literature on resilience in the forestry context, synthesizing how resilience is defined and assessed.

RECENT FINDINGS

Based on a detailed review of 255 studies, we analysed how the concepts of engineering resilience, ecological resilience and social-ecological resilience are used in forest sciences. A clear majority of the studies applied the concept of engineering resilience, quantifying resilience as the recovery time after a disturbance. The two most used indicators for engineering resilience were basal area increment and vegetation cover, whereas ecological resilience studies frequently focus on vegetation cover and tree density. In contrast, important social-ecological resilience indicators used in the literature are socioeconomic diversity and stock of natural resources. In the context of global change, we expected an increase in studies adopting the more holistic social-ecological resilience concept, but this was not the observed trend.

SUMMARY

Our analysis points to the nestedness of these three resilience concepts, suggesting that they are complementary rather than contradictory. It also means that the variety of resilience approaches does not need to be an obstacle for operationalisation of the concept. We provide guidance for choosing the most suitable resilience concept and indicators based on the management, disturbance and application context.

Functional Relationships of Wood Anatomical Traits in Norway Spruce

The quantitative assessment of wood anatomical traits offers important insights into those factors that shape tree growth. While it is known that conduit diameter, cell wall thickness, and wood density vary substantially between and within species, the interconnection between wood anatomical traits, tree-ring width, tree height and age, as well as environment effects on wood anatomy remain unclear. Here, we measure and derived 65 wood anatomical traits in cross-sections of the five outermost tree rings (2008-2012) of 30 Norway spruce [Picea abies (L.) H. Karst.] trees growing along an altitudinal gradient (1,400-1,750 m a.s.l.) in the northern Apennines (Italy). We assess the relationship among each anatomical trait and between anatomical trait groups according to their function for (i) tree-ring growth, (ii) cell growth, (iii) hydraulic traits, and (iv) mechanical traits. The results show that tree height significantly affects wood hydraulic traits, as well as number and tangential diameter of tracheids, and ultimately the total ring width. Moreover, the amount of earlywood and latewood percentage influence wood hydraulic safety and efficiency, as well as mechanical traits. Mechanically relevant wood anatomical traits are mainly influenced by tree age, not necessarily correlated with tree height. An additional level of complexity is also indicated by some anatomical traits, such as latewood lumen diameter and the cell wall reinforcement index, showing large inter-annual variation as a proxy of phenotypic plasticity. This study unravels the complex interconnection of tree-ring tracheid structure and identifies anatomical traits showing a large inter-individual variation and a strong interannual coherency. Knowing and quantifying anatomical variation in cells of plant stem is crucial in ecological and biological studies for an appropriate interpretation of abiotic drivers of wood formation often related to tree height and/or tree age.

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.

Influence of starch deficiency on photosynthetic and post-photosynthetic carbon isotope fractionations

Carbon isotope (13C) fractionations occurring during and after photosynthetic CO2 fixation shape the carbon isotope composition (δ13C) of plant material and respired CO2. However, responses of 13C fractionations to diel variation in starch metabolism in the leaf are not fully understood. Here we measured δ13C of organic matter (δ13COM), concentrations and δ13C of potential respiratory substrates, δ13C of dark-respired CO2 (δ13CR), and gas exchange in leaves of starch-deficient plastidial phosphoglucomutase (pgm) mutants and wild-type plants of four species (Arabidopsis thaliana, Mesembryanthemum crystallinum, Nicotiana sylvestris, and Pisum sativum). The strongest δ13C response to the pgm-induced starch deficiency was observed in N. sylvestris, with more negative δ13COM, δ13CR, and δ13C values for assimilates (i.e. sugars and starch) and organic acids (i.e. malate and citrate) in pgm mutants than in wild-type plants during a diel cycle. The genotype differences in δ13C values could be largely explained by differences in leaf gas exchange. In contrast, the PGM-knockout effect on post-photosynthetic 13C fractionations via the plastidic fructose-1,6-bisphosphate aldolase reaction or during respiration was small. Taken together, our results show that the δ13C variations in starch-deficient mutants are primarily explained by photosynthetic 13C fractionations and that the combination of knockout mutants and isotope analyses allows additional insights into plant metabolism.

Disentangling the effects of spatial proximity and genetic similarity on individual growth performances in Norway spruce natural populations

Cambial growth is a phenotypic trait influenced by various physiological processes, numerous biotic and abiotic drivers, as well as by the genetic background. By archiving the outcome of such complex interplay, tree-rings are an exceptional resource for addressing individual long-term growth responses to changing environments and climate. Disentangling the effects of the different drivers of tree growth, however, remains challenging because of the lack of multidisciplinary data. Here, we combine individual dendrochronological, genetic and spatial data to assess the relative importance of genetic similarity and spatial proximity on Norway spruce (Picea abies (L.) Karst.) growth performances. We intensively sampled five plots from two populations in southern and central Europe, characterizing a total of 482 trees. A two-step analytical framework was developed. First, the effects of climate and tree age on tree-ring width (TRW) were estimated for each individual using a random slope linear mixed-effects model. Individual parameters were then tested against genetic and spatial variables by Mantel tests, partial redundancy analyses and variance partitioning. Our modelling approach successfully captured a large fraction of variance in TRW (conditional R² values up to 0.94) which was largely embedded in inter-individual differences. All statistical approaches consistently showed that genetic similarity was not related to variation in the individual parameters describing growth responses. In contrast, up to 29% of the variance of individual parameters was accounted by spatial variables, revealing that microenvironmental features are more relevant than genetic similarity in determining similar growth patterns. Our study highlights both the advantages of modelling dendrochronological data at the individual level and the relevance of microenvironmental variation on individual growth patterns. These two aspects should be carefully considered in future multidisciplinary studies on growth dynamics in natural populations.

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.