Bridging the gap between ecophysiological and genetic knowledge to assess the adaptive potential of European beech

Climate legacy in seed and seedling traits of European beech populations

Tree species' ability to persist within their current distribution ranges is determined by seed germination and seedling growth. Exploring variation in these traits in relation to climatic conditions helps to understand and predict tree population dynamics, and to support species management and conservation under future climate. We analyzed seeds and seedlings of 26 European beech populations from the northeastern boundary of the species range to test whether: 1) adaptation to climatic conditions is reflected in depth of dormancy and germination of seeds; 2) climatic characteristics of origin predictably affect seedling traits. The variation in seed dormancy and germination in a laboratory test, and seedling growth and morphology traits in a nursery common-garden test was examined. Populations originating from warmer and drier sites (mostly from the northern region), compared to those from the opposite end of climatic gradient, germinated later, with a lower success, and produced seedlings with shorter and tougher roots. They had deeper dormancy and poorer seed germination capacity, and are likely more vulnerable to environmental changes. The climatic conditions at the origin shape the intraspecific variation of seed germination and seedling traits, and may limit regeneration from seed and affect adaptation potential of beech to increasing temperatures and decreasing precipitation.

Demographic feedbacks during evolutionary rescue can slow or speed adaptive evolution

Populations declining toward extinction can persist via genetic adaptation in a process called evolutionary rescue. Predicting evolutionary rescue has applications ranging from conservation biology to medicine, but requires understanding and integrating the multiple effects of a stressful environmental change on population processes. Here we derive a simple expression for how generation time, a key determinant of the rate of evolution, varies with population size during evolutionary rescue. Change in generation time is quantitatively predicted by comparing how intraspecific competition and the source of maladaptation each affect the rates of births and deaths in the population. Depending on the difference between two parameters quantifying these effects, the model predicts that populations may experience substantial changes in their rate of adaptation in both positive and negative directions, or adapt consistently despite severe stress. These predictions were then tested by comparison to the results of individual-based simulations of evolutionary rescue, which validated that the tolerable rate of environmental change varied considerably as described by analytical results. We discuss how these results inform efforts to understand wildlife disease and adaptation to climate change, evolution in managed populations and treatment resistance in pathogens.

Forest Dynamics Models for Conservation, Restoration, and Management of Small Forests

Globally, there are myriad situations in which people aim to conserve, restore, or manage forest ecosystems at small spatial scales of 50 ha or less. To inform management, forest dynamics models provide an increasingly diverse and valuable portfolio of tools for projecting forest change under different management and environmental conditions. Yet, many models may not be appropriate or feasible to use in small forest management because of their design for larger-scale applications, the information needed to initialize models, or discrepancies between model outputs and information relevant for small forest management objectives. This review explores the suitability of 54 existing forest dynamics models to inform the management of small forests. We evaluated the characteristics of each model using five criteria with implications for small forest management: spatial resolution, number of species the model can simulate, inclusion of spatial structure, modeling approach, and mechanistic detail. While numerous models can be suitable under certain conditions, the review criteria led us to conclude that two models offered the broadest versatility and usability for small forest contexts, SORTIE and FORMIND. This review can help orient and guide small forest managers who wish to add modeling to their forest management efforts.

Genomic signatures of natural selection at phenology-related genes in a widely distributed tree species Fagus sylvatica L

BACKGROUND

Diversity among phenology-related genes is predicted to be a contributing factor in local adaptations seen in widely distributed plant species that grow in climatically variable geographic areas, such as forest trees. European beech (Fagus sylvatica L.) is widespread, and is one of the most important broadleaved tree species in Europe; however, its potential for adaptation to climate change is a matter of uncertainty, and little is known about the molecular basis of climate change-relevant traits like bud burst.

RESULTS

We explored single nucleotide polymorphisms (SNP) at candidate genes related to bud burst in beech individuals sampled across 47 populations from Europe. SNP diversity was monitored for 380 candidate genes using a sequence capture approach, providing 2909 unlinked SNP loci. We used two complementary analytical methods to find loci significantly associated with geographic variables, climatic variables (expressed as principal components), or phenotypic variables (spring and autumn phenology, height, survival). Redundancy analysis (RDA) was used to detect candidate markers across two spatial scales (entire study area and within subregions). We revealed 201 candidate SNPs at the broadest scale, 53.2% of which were associated with phenotypic variables. Additive polygenic scores, which provide a measure of the cumulative signal across significant candidate SNPs, were correlated with a climate variable (first principal component, PC1) related to temperature and precipitation availability, and spring phenology. However, different genotype-environment associations were identified within Southeastern Europe as compared to the entire geographic range of European beech.

CONCLUSIONS

Environmental conditions play important roles as drivers of genetic diversity of phenology-related genes that could influence local adaptation in European beech. Selection in beech favors genotypes with earlier bud burst under warmer and wetter habitats within its range; however, selection pressures may differ across spatial scales.

The within-population variability of leaf spring and autumn phenology is influenced by temperature in temperate deciduous trees

Leaf phenology is a major driver of ecosystem functioning in temperate forests and a robust indicator of climate change. Both the inter-annual and inter-population variability of leaf phenology have received much attention in the literature; in contrast, the within-population variability of leaf phenology has been far less studied. Beyond its impact on individual tree physiological processes, the within-population variability of leaf phenology can affect the estimation of the average budburst or leaf senescence dates at the population scale. Here, we monitored the progress of spring and autumn leaf phenology over 14 tree populations (9 tree species) in six European forests over the period of 2011 to 2018 (yielding 16 site-years of data for spring, 14 for autumn). We monitored 27 to 512 (with a median of 62) individuals per population. We quantified the within-population variability of leaf phenology as the standard deviation of the distribution of individual dates of budburst or leaf senescence (SD_BBi and SD_LSi, respectively). Given the natural variability of phenological dates occurring in our tree populations, we estimated from the data that a minimum sample size of 28 (resp. 23) individuals, are required to estimate SD_BBi (resp. SD_LSi) with a precision of 3 (resp. 7) days. The within-population of leaf senescence (average SD_LSi = 8.5 days) was on average two times larger than for budburst (average SD_BBi = 4.0 days). We evidenced that warmer temperature during the budburst period and a late average budburst date were associated with a lower SD_BBi, as a result of a quicker spread of budburst in tree populations, with a strong species effect. Regarding autumn phenology, we observed that later senescence and warm temperatures during the senescence period were linked with a high SD_LSi, with a strong species effect. The shares of variance explained by our models were modest suggesting that other factors likely influence the within-population variation in leaf phenology. For instance, a detailed analysis revealed that summer temperatures were negatively correlated with a lower SD_LSi.

"Green pointillism": detecting the within-population variability of budburst in temperate deciduous trees with phenological cameras

Phenological cameras have been used over a decade for identifying plant phenological markers (budburst, leaf senescence) and more generally the greenness dynamics of forest canopies. The analysis is usually carried out over the full camera field of view, with no particular analysis of the variability of phenological markers among trees. Here we show that images produced by phenological cameras can be used to quantify the within-population variability of budburst (WPVbb) in temperate deciduous forests. Using seven site-years of image analyses, we report a strong correlation (r² = 0.97) between the WPVbb determined with a phenological camera and its quantification through ground observation. We show that WPVbb varies strongly (by a factor of 4) from year to year in a given population and that those variations are linked with temperature conditions during the budburst period, with colder springs associated to a higher differentiation of budburst (higher WPVbb) among trees. Deploying our approach at the continental scale, i.e., throughout phenological cameras networks, would improve the understanding of the spatial (across populations) and temporal (across years) variations of WPVbb, which have strong implications on forest functioning, tree fitness and phenological modelling.

Towards a New Generation of Trait-Flexible Vegetation Models

Plant trait variability, emerging from eco-evolutionary dynamics that range from alleles to macroecological scales, is one of the most elusive, but possibly most consequential, aspects of biodiversity. Plasticity, epigenetics, and genetic diversity are major determinants of how plants will respond to climate change, yet these processes are rarely represented in current vegetation models. Here, we provide an overview of the challenges associated with understanding the causes and consequences of plant trait variability, and review current developments to include plasticity and evolutionary mechanisms in vegetation models. We also present a roadmap of research priorities to develop a next generation of vegetation models with flexible traits. Including trait variability in vegetation models is necessary to better represent biosphere responses to global change.

Forest carbon allocation modelling under climate change

Carbon allocation plays a key role in ecosystem dynamics and plant adaptation to changing environmental conditions. Hence, proper description of this process in vegetation models is crucial for the simulations of the impact of climate change on carbon cycling in forests. Here we review how carbon allocation modelling is currently implemented in 31 contrasting models to identify the main gaps compared with our theoretical and empirical understanding of carbon allocation. A hybrid approach based on combining several principles and/or types of carbon allocation modelling prevailed in the examined models, while physiologically more sophisticated approaches were used less often than empirical ones. The analysis revealed that, although the number of carbon allocation studies over the past 10 years has substantially increased, some background processes are still insufficiently understood and some issues in models are frequently poorly represented, oversimplified or even omitted. Hence, current challenges for carbon allocation modelling in forest ecosystems are (i) to overcome remaining limits in process understanding, particularly regarding the impact of disturbances on carbon allocation, accumulation and utilization of nonstructural carbohydrates, and carbon use by symbionts, and (ii) to implement existing knowledge of carbon allocation into defence, regeneration and improved resource uptake in order to better account for changing environmental conditions.

Local Adaptation and Response of Platycladus orientalis (L.) Franco Populations to Climate Change

Knowledge about the local adaptation and response of forest tree populations to the climate is important for assessing the impact of climate change and developing adaptive genetic resource management strategies. However, such information is not available for most plant species. Here, based on 69 provenances tested at 19 common garden experimental sites, we developed a universal response function (URF) for tree height at seven years of age for the important and wide-spread native Chinese tree species Platycladus orientalis (L.) Franco. URF was recently used to predict the potential growth response of a population originating from any climate and growing in any climate conditions. The developed model integrated both genetic and environmental effects, and explained 55% of the total variation in tree height observed among provenances and test sites in China. We found that local provenances performed better than non-local counterparts in habitats located in central, eastern, and southwestern China, showing the evidence of local adaptation as compared to other regions. In contrast, non-local provenances outperformed local ones in peripheral areas in northern and northwestern China, suggesting an adaptational lag in these areas. Future projections suggest that the suitable habitat areas of P. orientalis would expand by 15%–39% and shift northward by 0.8–3 degrees in latitude; however, the projected tree height of this species would decline by 4%–8% if local provenances were used. If optimal provenances were used, tree height growth could be improved by 13%–15%, along with 59%–71% suitable habitat expansion. Thus, assisted migration with properly selected seed sources would be effective in avoiding maladaptation in new plantations under a changing climate for P. orientalis.

Seed and pollen gene dispersal in Taxus baccata, a dioecious conifer in the face of strong population fragmentation

Background and Aims

Dispersal is crucial due to its direct impact on dynamics of a species' distribution as well as having a role in shaping adaptive potential through gene flow. In plants forming scarce and small populations, knowledge about the dispersal process is required to assess the potential for colonizing new habitats and connectivity of present and future populations. This study aimed to assess dispersal potential in Taxus baccata, a dioecious gymnosperm tree with a wide but highly fragmented distribution.

Methods

Seed and pollen dispersal kernels were estimated directly in the framework of the spatially explicit mating model, where genealogies of naturally established seedlings were reconstructed with the help of microsatellite markers. In this way, six differently shaped dispersal functions were compared.

Key Results

Seed dispersal followed a leptokurtic distribution, with the Exponential-Power, the Power-law and Weibull being almost equally best-fitting models. The pollen dispersal kernel appeared to be more fat-tailed than the seed dispersal kernel, and the Lognormal and the Exponential-Power function showed the best fit. The rate of seed immigration from the background sources was not significantly different from the rate of pollen immigration (13.1 % vs. 19.7 %) and immigration rates were in agreement with or below maximum predictions based on the estimated dispersal kernels. Based on the multimodel approach, 95 % of seeds travel <109 m, while 95 % of pollen travels <704 m from the source.

Conclusions

The results showed that, at a local spatial scale, yew seeds travel shorter distances than pollen, facilitating a rapid development of a kinship structure. At the landscape level, however, although yew exhibits some potential to colonize new habitats through seed dispersal, genetic connectivity between different yew remnants is strongly limited. Taking into account strong population fragmentation, the study suggests that gene dispersal may be a limiting factor of the adaptability of the species.

Risk of genetic maladaptation due to climate change in three major European tree species

Tree populations usually show adaptations to their local environments as a result of natural selection. As climates change, populations can become locally maladapted and decline in fitness. Evaluating the expected degree of genetic maladaptation due to climate change will allow forest managers to assess forest vulnerability, and develop strategies to preserve forest health and productivity. We studied potential genetic maladaptation to future climates in three major European tree species, Norway spruce (Picea abies), silver fir (Abies alba), and European beech (Fagus sylvatica). A common garden experiment was conducted to evaluate the quantitative genetic variation in growth and phenology of seedlings from 77 to 92 native populations of each species from across Switzerland. We used multivariate genecological models to associate population variation with past seed source climates, and to estimate relative risk of maladaptation to current and future climates based on key phenotypic traits and three regional climate projections within the A1B scenario. Current risks from climate change were similar to average risks from current seed transfer practices. For all three climate models, future risks increased in spruce and beech until the end of the century, but remained low in fir. Largest average risks associated with climate projections for the period 2061-2090 were found for spruce seedling height (0.64), and for beech bud break and leaf senescence (0.52 and 0.46). Future risks for spruce were high across Switzerland. However, areas of high risk were also found in drought-prone regions for beech and in the southern Alps for fir. Genetic maladaptation to future climates is likely to become a problem for spruce and beech by the end of this century, but probably not for fir. Consequently, forest management strategies should be adjusted in the study area for spruce and beech to maintain productive and healthy forests in the future.

Can phenological models predict tree phenology accurately in the future? The unrevealed hurdle of endodormancy break

The onset of the growing season of trees has been earlier by 2.3 days per decade during the last 40 years in temperate Europe because of global warming. The effect of temperature on plant phenology is, however, not linear because temperature has a dual effect on bud development. On one hand, low temperatures are necessary to break bud endodormancy, and, on the other hand, higher temperatures are necessary to promote bud cell growth afterward. Different process-based models have been developed in the last decades to predict the date of budbreak of woody species. They predict that global warming should delay or compromise endodormancy break at the species equatorward range limits leading to a delay or even impossibility to flower or set new leaves. These models are classically parameterized with flowering or budbreak dates only, with no information on the endodormancy break date because this information is very scarce. Here, we evaluated the efficiency of a set of phenological models to accurately predict the endodormancy break dates of three fruit trees. Our results show that models calibrated solely with budbreak dates usually do not accurately predict the endodormancy break date. Providing endodormancy break date for the model parameterization results in much more accurate prediction of this latter, with, however, a higher error than that on budbreak dates. Most importantly, we show that models not calibrated with endodormancy break dates can generate large discrepancies in forecasted budbreak dates when using climate scenarios as compared to models calibrated with endodormancy break dates. This discrepancy increases with mean annual temperature and is therefore the strongest after 2050 in the southernmost regions. Our results claim for the urgent need of massive measurements of endodormancy break dates in forest and fruit trees to yield more robust projections of phenological changes in a near future.

Bring in the genes: genetic-ecophysiological modeling of the adaptive response of trees to environmental change. With application to the annual cycle

The observation of strong latitudinal clines in the date of bud burst of tree species indicate that populations of these species are genetically adapted to local environmental conditions. Existing phenological models rarely address this clinal variation, so that adaptive responses of tree populations to changes in environmental conditions are not taken into account, e.g., in models on species distributions that use phenological sub-models. This omission of simulating adaptive response in tree models may over- or underestimate the effects of climate change on tree species distributions, as well as the impacts of climate change on tree growth and productivity. Here, we present an approach to model the adaptive response of traits to environmental change based on an integrated process-based eco-physiological and quantitative genetic model of adaptive traits. Thus, the parameter values of phenological traits are expressed in genetic terms (allele effects and-frequencies, number of loci) for individual trees. These individual trees thereby differ in their ability to acquire resources, grow and reproduce as described by the process-based model, leading to differential survival. Differential survival is thus the consequence of both differences in parameters values and their genetic composition. By simulating recombination and dispersal of pollen, the genetic composition of the offspring will differ from that of their parents. Over time, the distribution of both trait values and the frequency of the underlying alleles in the population change as a consequence of changes in environmental drivers leading to adaptation of trees to local environmental conditions. This approach is applied to an individual-tree growth model that includes a phenological model on the annual cycle of trees whose parameters are allowed to adapt. An example of the adaptive response of the onset of the growing season across Europe is presented.

Simulating local adaptation to climate of forest trees with a Physio-Demo-Genetics model

One challenge of evolutionary ecology is to predict the rate and mechanisms of population adaptation to environmental variations. The variations in most life history traits are shaped both by individual genotypic and by environmental variation. Forest trees exhibit high levels of genetic diversity, large population sizes, and gene flow, and they also show a high level of plasticity for life history traits. We developed a new Physio-Demo-Genetics model (denoted PDG) coupling (i) a physiological module simulating individual tree responses to the environment; (ii) a demographic module simulating tree survival, reproduction, and pollen and seed dispersal; and (iii) a quantitative genetics module controlling the heritability of key life history traits. We used this model to investigate the plastic and genetic components of the variations in the timing of budburst (TBB) along an elevational gradient of Fagus sylvatica (the European beech). We used a repeated 5 years climatic sequence to show that five generations of natural selection were sufficient to develop nonmonotonic genetic differentiation in the TBB along the local climatic gradient but also that plastic variation among different elevations and years was higher than genetic variation. PDG complements theoretical models and provides testable predictions to understand the adaptive potential of tree populations.

Within-population genetic structure in beech (Fagus sylvatica L.) stands characterized by different disturbance histories: does forest management simplify population substructure?

The fine-scale assessment of both spatially and non-spatially distributed genetic variation is crucial to preserve forest genetic resources through appropriate forest management. Cryptic within-population genetic structure may be more common than previously thought in forest tree populations, which has strong implications for the potential of forests to adapt to environmental change. The present study was aimed at comparing within-population genetic structure in European beech (Fagus sylvatica L.) plots experiencing different disturbance levels. Five plot pairs made up by disturbed and undisturbed plots having the same biogeographic history were sampled throughout Europe. Overall, 1298 individuals were analyzed using four highly polymorphic nuclear microsatellite markers (SSRs). Bayesian clustering within plots identified 3 to 11 genetic clusters (within-plot θ_ST ranged from 0.025 to 0.124). The proportion of within-population genetic variation due to genetic substructuring (F_CluPlot = 0.067) was higher than the differentiation among the 10 plots (F_PlotTot = 0.045). Focusing on the comparison between managed and unmanaged plots, disturbance mostly explains differences in the complexity of within-population genetic structure, determining a reduction of the number of genetic clusters present in a standardized area. Our results show that: i) genetic substructuring needs to be investigated when studying the within-population genetic structure in forest tree populations, and ii) indices describing subtle characteristics of the within-population genetic structure are good candidates for providing early signals of the consequences of forest management, and of disturbance events in general.

A road map for integrating eco-evolutionary processes into biodiversity models

The demand for projections of the future distribution of biodiversity has triggered an upsurge in modelling at the crossroads between ecology and evolution. Despite the enthusiasm around these so-called biodiversity models, most approaches are still criticised for not integrating key processes known to shape species ranges and community structure. Developing an integrative modelling framework for biodiversity distribution promises to improve the reliability of predictions and to give a better understanding of the eco-evolutionary dynamics of species and communities under changing environments. In this article, we briefly review some eco-evolutionary processes and interplays among them, which are essential to provide reliable projections of species distributions and community structure. We identify gaps in theory, quantitative knowledge and data availability hampering the development of an integrated modelling framework. We argue that model development relying on a strong theoretical foundation is essential to inspire new models, manage complexity and maintain tractability. We support our argument with an example of a novel integrated model for species distribution modelling, derived from metapopulation theory, which accounts for abiotic constraints, dispersal, biotic interactions and evolution under changing environmental conditions. We hope such a perspective will motivate exciting and novel research, and challenge others to improve on our proposed approach.

Long-distance gene flow and adaptation of forest trees to rapid climate change

Forest trees are the dominant species in many parts of the world and predicting how they might respond to climate change is a vital global concern. Trees are capable of long-distance gene flow, which can promote adaptive evolution in novel environments by increasing genetic variation for fitness. It is unclear, however, if this can compensate for maladaptive effects of gene flow and for the long-generation times of trees. We critically review data on the extent of long-distance gene flow and summarise theory that allows us to predict evolutionary responses of trees to climate change. Estimates of long-distance gene flow based both on direct observations and on genetic methods provide evidence that genes can move over spatial scales larger than habitat shifts predicted under climate change within one generation. Both theoretical and empirical data suggest that the positive effects of gene flow on adaptation may dominate in many instances. The balance of positive to negative consequences of gene flow may, however, differ for leading edge, core and rear sections of forest distributions. We propose future experimental and theoretical research that would better integrate dispersal biology with evolutionary quantitative genetics and improve predictions of tree responses to climate change.

Reconstruction of a beech population bottleneck using archival demographic information and Bayesian analysis of genetic data

Range expansion and contraction has occurred in the history of most species and can seriously impact patterns of genetic diversity. Historical data about range change are rare and generally appropriate for studies at large scales, whereas the individual pollen and seed dispersal events that form the basis of geneflow and colonization generally occur at a local scale. In this study, we investigated range change in Fagus sylvatica on Mont Ventoux, France, using historical data from 1838 to the present and approximate Bayesian computation (ABC) analyses of genetic data. From the historical data, we identified a population minimum in 1845 and located remnant populations at least 200 years old. The ABC analysis selected a demographic scenario with three populations, corresponding to two remnant populations and one area of recent expansion. It also identified expansion from a smaller ancestral population but did not find that this expansion followed a population bottleneck, as suggested by the historical data. Despite a strong support to the selected scenario for our data set, the ABC approach showed a low power to discriminate among scenarios on average and a low ability to accurately estimate effective population sizes and divergence dates, probably due to the temporal scale of the study. This study provides an unusual opportunity to test ABC analysis in a system with a well-documented demographic history and identify discrepancies between the results of historical, classical population genetic and ABC analyses. The results also provide valuable insights into genetic processes at work at a fine spatial and temporal scale in range change and colonization.