Genetic differentiation across a steep and narrow environmental gradient: Quantitative genetic and genomic insights into Lake Superior populations of Quercus rubra

Adaptive differentiation of traits and underlying loci can occur at a small geographical scale if natural selection is stronger than countervailing gene flow and drift. We investigated this hypothesis using coupled quantitative genetic and genomic approaches for a wind-pollinated tree species, Quercus rubra, along the steep, narrow gradient of the Lake Superior coast that encompasses four USDA Hardiness Zones within 100 km. For the quantitative genetic component of this study, we examined phenotypic differentiation among eight populations in a common garden, measuring seed mass, germination, height, stem diameter, leaf number, specific leaf area and survival. For the genomic component, we quantified genetic differentiation for 26 populations from the same region using RAD-seq. Because hybridisation with Quercus ellipsoidalis occurs in other parts of the species' range, we included two populations of this congener for comparison. In the common garden study, we found a strong signal of population differentiation that was significantly associated with at least one climate factor for nine of 10 measured traits. In contrast, we found no evidence of genomic differentiation among populations based on FST or any other measures. However, both distance-based and genotype-environment association analyses identified loci showing the signature of selection, with one locus in common across five analyses. This locus was associated with the minimum temperature of the coldest month, a factor that defines the climate zones and was also significant in the common garden analyses. In addition, we documented introgression from Q. ellipsoidalis into Q. rubra, with rates of introgression correlated with the climate gradient. In sum, this study reveals signatures of selection at the quantitative trait and genomic level consistent with climate adaptation, a pattern that is more often documented at a much broader geographical scale, especially in long-lived wind-pollinated species.

Clinal variation in quantitative traits but not in evolutionary potential along elevational and latitudinal gradients in the widespread Anthyllis vulneraria

PREMISE

Strong elevational and latitudinal gradients allow the study of genetic differentiation in response to similar environmental changes. However, it is uncertain whether the environmental changes along the two types of gradients result in similar genetically based changes in quantitative traits. Peripheral arctic and alpine populations are thought to have less evolutionary potential than more central populations do.

METHODS

We studied quantitative traits of the widespread Anthyllis vulneraria in a common garden. Plants originated from 20 populations along a 2000-m elevational gradient from the lowlands to the elevational limit of the species in the Alps, and from 20 populations along a 2400-km latitudinal gradient from the center of the distribution of the species in Central Europe to its northern distributional margin.

RESULTS

Most traits showed similar clinal variations with elevation and latitude of origin, and the magnitude of all measured traits in relation to mean annual temperature was similar. Higher QST values than FST values in several traits indicated diversifying selection, but for others QST was smaller than FST. Genetic diversity of quantitative traits and neutral molecular markers was not correlated. Plasticity in response to favorable conditions declined with elevation and less strongly with latitude of origin, but the evolvability of traits did not.

CONCLUSIONS

The clinal variation suggests adaptive differentiation of quantitative traits along the two gradients. The evolutionary potential of peripheral populations is not necessarily reduced, but lower plasticity may threaten their survival under rapidly changing climatic conditions.

Impacts of anthropogenic climate change on tropical montane forests: an appraisal of the evidence

In spite of their small global area and restricted distributions, tropical montane forests (TMFs) are biodiversity hotspots and important ecosystem services providers, but are also highly vulnerable to climate change. To protect and preserve these ecosystems better, it is crucial to inform the design and implementation of conservation policies with the best available scientific evidence, and to identify knowledge gaps and future research needs. We conducted a systematic review and an appraisal of evidence quality to assess the impacts of climate change on TMFs. We identified several skews and shortcomings. Experimental study designs with controls and long-term (≥10 years) data sets provide the most reliable evidence, but were rare and gave an incomplete understanding of climate change impacts on TMFs. Most studies were based on predictive modelling approaches, short-term (<10 years) and cross-sectional study designs. Although these methods provide moderate to circumstantial evidence, they can advance our understanding on climate change effects. Current evidence suggests that increasing temperatures and rising cloud levels have caused distributional shifts (mainly upslope) of montane biota, leading to alterations in biodiversity and ecological functions. Neotropical TMFs were the best studied, thus the knowledge derived there can serve as a proxy for climate change responses in under-studied regions elsewhere. Most studies focused on vascular plants, birds, amphibians and insects, with other taxonomic groups poorly represented. Most ecological studies were conducted at species or community levels, with a marked paucity of genetic studies, limiting understanding of the adaptive capacity of TMF biota. We thus highlight the long-term need to widen the methodological, thematic and geographical scope of studies on TMFs under climate change to address these uncertainties. In the short term, however, in-depth research in well-studied regions and advances in computer modelling approaches offer the most reliable sources of information for expeditious conservation action for these threatened forests.

Antagonistic Effects of Assortative Mating on the Evolution of Phenotypic Plasticity along Environmental Gradients

AbstractPrevious theory has shown that assortative mating for plastic traits can maintain genetic divergence across environmental gradients despite high gene flow. Yet these models did not examine how assortative mating affects the evolution of plasticity. We here describe patterns of genetic variation across elevation for plasticity in a trait under assortative mating, using multiple-year observations of budburst date in a common garden of sessile oaks. Despite high gene flow, we found significant spatial genetic divergence for the intercept, but not for the slope, of reaction norms to temperature. We then used individual-based simulations, where both the slope and the intercept of the reaction norm evolve, to examine how assortative mating affects the evolution of plasticity, varying the intensity and distance of gene flow. Our model predicts the evolution of either suboptimal plasticity (reaction norms with a slope shallower than optimal) or hyperplasticity (slopes steeper than optimal) in the presence of assortative mating when optimal plasticity would evolve under random mating. Furthermore, a cogradient pattern of genetic divergence for the intercept of the reaction norm (where plastic and genetic effects are in the same direction) always evolves in simulations with assortative mating, consistent with our observations in the studied oak populations.

Plant traits poorly predict winner and loser shrub species in a warming tundra biome

Climate change is leading to species redistributions. In the tundra biome, shrubs are generally expanding, but not all tundra shrub species will benefit from warming. Winner and loser species, and the characteristics that may determine success or failure, have not yet been fully identified. Here, we investigate whether past abundance changes, current range sizes and projected range shifts derived from species distribution models are related to plant trait values and intraspecific trait variation. We combined 17,921 trait records with observed past and modelled future distributions from 62 tundra shrub species across three continents. We found that species with greater variation in seed mass and specific leaf area had larger projected range shifts, and projected winner species had greater seed mass values. However, trait values and variation were not consistently related to current and projected ranges, nor to past abundance change. Overall, our findings indicate that abundance change and range shifts will not lead to directional modifications in shrub trait composition, since winner and loser species share relatively similar trait spaces.

Key triggers of adaptive genetic variability of sessile oak [Q. petraea (Matt.) Liebl.] from the Balkan refugia: outlier detection and association of SNP loci from ddRAD-seq data

Knowledge on the genetic composition of Quercus petraea in south-eastern Europe is limited despite the species' significant role in the re-colonisation of Europe during the Holocene, and the diverse climate and physical geography of the region. Therefore, it is imperative to conduct research on adaptation in sessile oak to better understand its ecological significance in the region. While large sets of SNPs have been developed for the species, there is a continued need for smaller sets of SNPs that are highly informative about the possible adaptation to this varied landscape. By using double digest restriction site associated DNA sequencing data from our previous study, we mapped RAD-seq loci to the Quercus robur reference genome and identified a set of SNPs putatively related to drought stress-response. A total of 179 individuals from eighteen natural populations at sites covering heterogeneous climatic conditions in the southeastern natural distribution range of Q. petraea were genotyped. The detected highly polymorphic variant sites revealed three genetic clusters with a generally low level of genetic differentiation and balanced diversity among them but showed a north-southeast gradient. Selection tests showed nine outlier SNPs positioned in different functional regions. Genotype-environment association analysis of these markers yielded a total of 53 significant associations, explaining 2.4-16.6% of the total genetic variation. Our work exemplifies that adaptation to drought may be under natural selection in the examined Q. petraea populations.

Oak stands along an elevation gradient have different molecular strategies for regulating bud phenology

BACKGROUND

Global warming raises serious concerns about the persistence of species and populations locally adapted to their environment, simply because of the shift it produces in their adaptive landscape. For instance, the phenological cycle of tree species may be strongly affected by higher winter temperatures and late frost in spring. Given the variety of ecosystem services they provide, the question of forest tree adaptation has received increasing attention in the scientific community and catalyzed research efforts in ecology, evolutionary biology and functional genomics to study their adaptive capacity to respond to such perturbations.

RESULTS

In the present study, we used an elevation gradient in the Pyrenees Mountains to explore the gene expression network underlying dormancy regulation in natural populations of sessile oak stands sampled along an elevation cline and potentially adapted to different climatic conditions mainly driven by temperature. By performing analyses of gene expression in terminal buds we identified genes displaying significant dormancy, elevation or dormancy-by-elevation interaction effects. Our Results highlighted that low- and high-altitude populations have evolved different molecular strategies for minimizing late frost damage and maximizing the growth period, thereby increasing potentially their respective fitness in these contrasting environmental conditions. More particularly, population from high elevation overexpressed genes involved in the inhibition of cell elongation and delaying flowering time while genes involved in cell division and flowering, enabling buds to flush earlier were identified in population from low elevation.

CONCLUSION

Our study made it possible to identify key dormancy-by-elevation responsive genes revealing that the stands analyzed in this study have evolved distinct molecular strategies to adapt their bud phenology in response to temperature.

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

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

Chasing genetic correlation breakers to stimulate population resilience to climate change

Global climate change introduces new combinations of environmental conditions, which is expected to increase stress on plants. This could affect many traits in multiple ways that are as yet unknown but will likely require the modification of existing genetic relationships among functional traits potentially involved in local adaptation. Theoretical evolutionary studies have determined that it is an advantage to have an excess of recombination events under heterogeneous environmental conditions. Our study, conducted on a population of radiata pine (Pinus radiata D. Don), was able to identify individuals that show high genetic recombination at genomic regions, which potentially include pleiotropic or collocating QTLs responsible for the studied traits, reaching a prediction accuracy of 0.80 in random cross-validation and 0.72 when whole family was removed from the training population and predicted. To identify these highly recombined individuals, a training population was constructed from correlation breakers, created through tandem selection of parents in the previous generation and their consequent mating. Although the correlation breakers showed lower observed heterogeneity possibly due to direct selection in both studied traits, the genomic regions with statistically significant differences in the linkage disequilibrium pattern showed higher level of heretozygosity, which has the effect of decomposing unfavourable genetic correlation. We propose undertaking selection of correlation breakers under current environmental conditions and using genomic predictions to increase the frequency of these ’recombined’ individuals in future plantations, ensuring the resilience of planted forests to changing climates. The increased frequency of such individuals will decrease the strength of the population-level genetic correlations among traits, increasing the opportunity for new trait combinations to be developed in the future.

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

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

Directional Selection on Tree Seedling Traits Driven by Experimental Drought Differs Between Mesic and Dry Populations

We evaluated population differences and drought-induced phenotypic selection on four seedling traits of the Australian forest tree Eucalyptus pauciflora using a glasshouse dry-down experiment. We compared dry and mesic populations and tested for directional selection on lamina length (reflecting leaf size), leaf shape, the node of ontogenetic transition to the petiolate leaf (reflecting the loss of vegetative juvenility), and lignotuber size (reflecting a recovery trait). On average, the dry population had smaller and broader leaves, greater retention of the juvenile leaf state and larger lignotubers than the mesic population, but the populations did not differ in seedling survival. While there was statistical support for directional selection acting on the focal traits in one or other population, and for differences between populations in selection gradient estimates for two traits, only one trait—lamina length—exhibited a pattern of directional selection consistent with the observed population differences being a result of past adaptation to reduce seedling susceptibility to acute drought. The observed directional selection for lamina length in the mesic population suggests that future increases in drought risk in the wild will shift the mean of the mesic population toward that of the dry population. Further, we provide evidence suggesting an early age trade-off between drought damage and recovery traits, with phenotypes which develop larger lignotubers early being more susceptible to drought death. Such trade-offs could have contributed to the absence of population mean differences in survival, despite marked differentiation in seedling traits.

Drone-based physiological index reveals long-term acclimation and drought stress responses in trees

Monitoring early tree physiological responses to drought is key to understanding progressive impacts of drought on forests and identifying resilient species. We combined drone-based multispectral remote sensing with measurements of tree physiology and environmental parameters over two growing seasons in a 100-y-old Pinus sylvestris forest subject to 17-y of precipitation manipulation. Our goal was to determine if drone-based photochemical reflectance index (PRI) captures tree drought stress responses and whether responses are affected by long-term acclimation. PRI detects changes in xanthophyll cycle pigment dynamics, which reflect increases in photoprotective non-photochemical quenching activity resulting from drought-induced photosynthesis downregulation. Here, PRI of never-irrigated trees was up to 10 times lower (higher stress) than PRI of irrigated trees. Long-term acclimation to experimental treatment, however, influenced the seasonal relationship between PRI and soil water availability. PRI also captured diurnal decreases in photochemical efficiency, driven by vapour pressure deficit. Interestingly, 5 years after irrigation was stopped for a subset of the irrigated trees, a positive legacy effect persisted, with lower stress responses (higher PRI) compared with never-irrigated trees. This study demonstrates the ability of remotely sensed PRI to scale tree physiological responses to an entire forest and the importance of long-term acclimation in determining current drought stress responses.

Adaptive responses to temperature and precipitation variation at the early-life stages of Pinus sylvestris

Early-stage fitness variation has been seldom evaluated at broad scales in forest tree species, despite the long tradition of studying climate-driven intraspecific genetic variation. In this study, we evaluated the role of climate in driving patterns of population differentiation at early-life stages in Pinus sylvestris and explored the fitness and growth consequences of seed transfer within the species range. We monitored seedling emergence, survival and growth over a 2-yr period in a multi-site common garden experiment which included 18 European populations and spanned 25° in latitude and 1700 m in elevation. Climate-fitness functions showed that populations exhibited higher seedling survival and growth at temperatures similar to their home environment, which is consistent with local adaptation. Northern populations experienced lower survival and growth at warmer sites, contrary to previous studies on later life stages. Seed mass was higher in populations from warmer areas and was positively associated with survival and growth at more southern sites. Finally, we did not detect a survival-growth trade-off; on the contrary, bigger seedlings exhibited higher survival probabilities under most climatic conditions. In conclusion, our results reveal that contrasting temperature regimes have played an important role in driving the divergent evolution of P. sylvestris populations at early-life stages.

Counter-gradient variation of reproductive effort in a widely distributed temperate oak (Quercus petraea)

The genetic and phenotypic variability of life history traits determines the demographic attributes of tree populations and, thus, their responses to anthropogenic climate change. Growth- and survival-related traits have been widely studied in forest ecology, but little is known about the determinism of reproductive traits.Using an elevation gradient experiment in the Pyrenees we assessed the degree to which variations in reproductive effort along climatic gradients are environmentally or genetically driven, by comparing oak populations (Quercus petraea) growing under field and common garden conditions.In situ monitoring revealed a decline in reproductive effort with increasing elevation and decreasing temperature. In common garden conditions, significant genetic differentiation was observed between provenances for reproduction and growth: trees from cold environments (high elevations) grew more slowly, and produced larger acorns in larger numbers. Our observations show that genetic and phenotypic clines for reproductive traits have opposite signs (counter-gradient) along the environmental gradient as opposed to growth, for which genetic variation parallels phenotypic variation (co-gradient).The counter-gradient found here for reproductive effort reveals that genetic variation partly counteracts the phenotypic effect of temperature, moderating the change in reproductive effort according to temperature. We consider the possible contribution to this counter-gradient in reproductive effort as an evolutionary trade-off between reproduction and growth.

Ecological Interactions, Environmental Gradients, and Gene Flow in Local Adaptation

Despite long-standing interest in local adaptation of plants to their biotic and abiotic environment, existing theory, and many case studies, little work to date has addressed within-species evolution of concerted strategies and how these might contrast with patterns across species. Here we consider the interactions between pollinators, herbivores, and resource availability in shaping plant local adaptation, how these interactions impact plant phenotypes and gene flow, and the conditions where multiple traits align along major environmental gradients such as latitude and elevation. Continued work in emerging model systems will benefit from the melding of classic experimental approaches with novel population genetic analyses to reveal patterns and processes in plant local adaptation.

Comparing Genetic Diversity in Three Threatened Oaks

Genetic diversity is a critical resource for species’ survival during times of environmental change. Conserving and sustainably managing genetic diversity requires understanding the distribution and amount of genetic diversity (in situ and ex situ) across multiple species. This paper focuses on three emblematic and IUCN Red List threatened oaks (Quercus, Fagaceae), a highly speciose tree genus that contains numerous rare species and poses challenges for ex situ conservation. We compare the genetic diversity of three rare oak species—Quercus georgiana, Q. oglethorpensis, and Q. boyntonii—to common oaks; investigate the correlation of range size, population size, and the abiotic environment with genetic diversity within and among populations in situ; and test how well genetic diversity preserved in botanic gardens correlates with geographic range size. Our main findings are: (1) these three rare species generally have lower genetic diversity than more abundant oaks; (2) in some cases, small population size and geographic range correlate with genetic diversity and differentiation; and (3) genetic diversity currently protected in botanic gardens is inadequately predicted by geographic range size and number of samples preserved, suggesting non-random sampling of populations for conservation collections. Our results highlight that most populations of these three rare oaks have managed to avoid severe genetic erosion, but their small size will likely necessitate genetic management going forward.

Revisiting the Functional Zoning Concept under Climate Change to Expand the Portfolio of Adaptation Options

Climate change is threatening our ability to manage forest ecosystems sustainably. Despite strong consensus on the need for a broad portfolio of options to face this challenge, diversified management options have yet to be widely implemented. Inspired by functional zoning, a concept aimed at optimizing biodiversity conservation and wood production in multiple-use forest landscapes, we present a portfolio of management options that intersects management objectives with forest vulnerability to better address the wide range of goals inherent to forest management under climate change. Using this approach, we illustrate how different adaptation options could be implemented when faced with impacts related to climate change and its uncertainty. These options range from establishing ecological reserves in climatic refuges, where self-organizing ecological processes can result in resilient forests, to intensive plantation silviculture that could ensure a stable wood supply in an uncertain future. While adaptation measures in forests that are less vulnerable correspond to the traditional functional zoning management objectives, forests with higher vulnerability might be candidates for transformative measures as they may be more susceptible to abrupt changes in structure and composition. To illustrate how this portfolio of management options could be applied, we present a theoretical case study for the eastern boreal forest of Canada. Even if these options are supported by solid evidence, their implementation across the landscape may present some challenges and will require good communication among stakeholders and with the public.

Low Population Differentiation but High Phenotypic Plasticity of European Beech in Germany

Drought is increasingly impairing the vitality of European beech (Fagus sylvatica L.) in several regions of its distribution range. In times of climate change, adaptive traits such as plant phenology and frost tolerance are also becoming more important. Adaptive patterns of European beech seem to be complex, as contrasting results regarding the relative effect of phenotypic plasticity and genetic variation in trait variation have been reported. Here, we used a large translocation experiment comprising more than 15,500 seedlings in three regions of Germany to investigate local adaptation and phenotypic plasticity in beech. We found low population differentiation regarding plant survival, and plant height increment, but high phenotypic plasticity for these traits. Survival showed a positive correlation with temperature variables and a less pronounced and negative correlation with precipitation-related variables. This suggests a predominant effect of temperature and growing degree days on the survival of beech seedlings under moderate drought stress. The high phenotypic plasticity may help beech to cope with changing environmental conditions, albeit increasing drought stress may make adaptive changes necessary in the long term.

Juniperus communis populations exhibit low variability in hydraulic safety and efficiency

The performance and distribution of woody species strongly depend on their adjustment to environmental conditions based on genotypic and phenotypic properties. Since more intense and frequent drought events are expected due to climate change, xylem hydraulic traits will play a key role under future conditions, and thus, knowledge of hydraulic variability is of key importance. In this study, we aimed to investigate the variability in hydraulic safety and efficiency of the conifer shrub Juniperus communis based on analyses along an elevational transect and a common garden approach. We studied (i) juniper plants growing between 700 and 2000 m a.s.l. Innsbruck, Austria, and (ii) plants grown in the Innsbruck botanical garden (Austria) from seeds collected at different sites across Europe (France, Austria, Ireland, Germany and Sweden). Due to contrasting environmental conditions at different elevation and provenance sites and the wide geographical study area, pronounced variation in xylem hydraulics was expected. Vulnerability to drought-induced embolisms (hydraulic safety) was assessed via the Cavitron and ultrasonic acoustic emission techniques, and the specific hydraulic conductivity (hydraulic efficiency) via flow measurements. Contrary to our hypothesis, relevant variability in hydraulic safety and efficiency was neither observed across elevations, indicating a low phenotypic variation, nor between provenances, despite expected genotypic differences. Interestingly, the provenance from the most humid and warmest site (Ireland) and the northernmost provenance (Sweden) showed the highest and the lowest embolism resistance, respectively. The hydraulic conductivity was correlated with plant height, which indicates that observed variation in hydraulic traits was mainly related to morphological differences between plants. We encourage future studies to underlie anatomical traits and the role of hydraulics for the broad ecological amplitude of J. communis.

What do you mean "functional" in ecology? Patterns versus processes

Use of the term "functional" trait has increased exponentially in ecology. Although accounting for numerous ecological questions, this concept raises several issues. We propose that the term "functional" could be misleading because (1) no rigorous criteria exist to identify "functional" traits and (2) it suggests that only some traits ("functional" ones) can inform our understanding of species functioning, whatever the scale or discipline. Hence, the concept of "functional" trait in ecology is starting to be challenged and it remains unclear why some traits should be considered functional, whereas other traits should not. We argue that the most used "functional" traits are meaningful because they reflect important differences between populations or species, based on synchronic comparisons, that is, irrespective of time (hereafter "pattern" traits). Hence, they are useful for identifying trade-offs and strategies across large numbers of observations, usually at rather coarse scales, and are most often used in analyses of "big data." However, given that many ecological processes occur across short time scales and narrow gradients of climate and resource availability, the efficacy of these traits to inform us about these ecological processes appears questionable. We show that trait measurements that take time explicitly into account (hereafter "process" traits) differ from pattern traits because they quantify the flows of material and energy within a given environment across a defined period of time. Although pattern traits and process traits are both functional, it is important to understand the differences between the approaches. Moreover, better accounting of ontogeny, life form, plasticity, and genetic variability is required to enhance the convergence between pattern and process approaches. This revised framework allows more explicit connections between trait ecology and other biological sciences. It should enhance the study of processes at all scales in order to investigate efficiently the adaptive responses of biological organisms to climate change.

Monuments Unveiled: Genetic Characterization of Large Old Chestnut (Castanea sativa Mill.) Trees Using Comparative Nuclear and Chloroplast DNA Analysis

Large old trees are extraordinary organisms. They not only represent a historical, landscape and environmental heritage of inestimable value, but they also witness a long history of environmental changes and human interventions, and constitute an as yet poorly known reserve of genetic variability which can be considered a great resource for management programs of forest species. This is the first genetic study on Italian, large, old chestnut trees (Castanea sativa Mill.). Ninety-nine trees were surveyed and analysed. For each tree, more than one sample from canopy and root suckers was collected to test for the genetic integrity of the individuals. All samples were genotyped using nine nuclear microsatellite markers (nSSRs) and 106 unique genetic profiles were identified. A Bayesian analysis performed with the software STRUCTURE revealed the occurrence of two main gene pools and unveiled the genetic relationships existing among the genotyped individuals, and with the natural chestnut populations living in proximity. A phylogeographic structure of the plastid diversity was also obtained by the use of DNA sequence variation at two marker regions, revealing different origins and probable connections of the old trees with different glacial refugia. Our results contribute to an improved evaluation of the European chestnut genetic resources and provide useful insights into the species’ history and domestication in Italy. The importance of carefully targeted conservation strategies for these invaluable organisms is reaffirmed.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Varying Levels of Genetic Control and Phenotypic Plasticity in Timing of Bud Burst, Flower Opening, Leaf Senescence and Leaf Fall in Two Common Gardens of Prunus padus L.

Several phenological phases mark the seasonal growth pattern in temperate woody perennials. To gain further insight into the way these phases react on an altering growth environment, we tested whether vegetative and reproductive phenophases in a shrub species respond differentially among different genetic entities and between two different planting sites. We scored leaf bud burst, flower opening, leaf senescence and leaf fall on 267 ramets of Prunus padus L. belonging to 53 genotypes that were sampled in 9 local populations, and that were planted in 2 common gardens in the northern part of Belgium. The data were processed with cumulative logistic regression. The contribution of genetic and non-genetic components to the total variability varied between the four studied seasonal phenophases. The timing of flower opening displayed the smallest relative amount of intragenotypic variance (between ramets), suggesting a stronger genetic control and a lesser need at the individual plant level for plastic fine tuning to the micro-environment. In addition, whereas leaf bud burst showed the highest relative variance at the interpopulation level among all phenophases, probably at least partly attributable to local adaptation, flower opening displayed the highest intergenotypic variance, which may have been promoted more by assortative mating. Spring phenophases were strongly correlated (r = 0.89) as well as the autumnal phenophases (r = 0.72). Flower opening was not correlated with the autumnal phenophases. Timing of leaf bud burst and leaf senescence were negatively correlated, demonstrating that the length of the growing season enlarged or diminished among the studied genotypes. Although the two planting sites were only 24 km apart, all phenophases were advanced at the less exposed site, indicating a phenotypic plastic response. Together, our results suggest that in P. padus, flowering is less sensitive to environmental variation than leaf bud break and may show a lesser impact of a changing environment on this reproductive phenophase.

Genetic variation in early fitness traits across European populations of silver birch (Betula pendula)

Given that the ecological niche of tree species is typically narrower for earlier life stages, intraspecific genetic variation at early fitness traits may greatly influence the adaptive response of tree populations to changing environmental conditions. In this study, we evaluated genetic variation in early fitness traits among 12 populations of Betula pendula from a wide latitudinal range in Europe (41-55°N). We first conducted a chamber experiment to test for population differences in germination and the effect of pre-chilling treatment on seed dormancy release. We then established three common gardens spread across the species latitudinal range in order to evaluate levels of quantitative genetic variation and genotype-by-environment interaction at different early life traits. Our results showed significant variation in chamber germination rates among populations (0-60 %), with southern populations exhibiting lower germination. Pre-chilling treatments did not generally improve germination success. Population seedling emergence rates in the field were correlated with chamber germination rates, though being an order of magnitude lower, with an average ranging from 0 to 1.3 % across gardens. Highly significant variation was found in field emergence rates among populations, and between seed-crop years within populations, but not among families within populations. Populations differed in seedling height, diameter, slenderness and budburst date, with significant among-family variation. Population latitude was positively associated with chamber germination rate and with seedling emergence rate in one of the central field sites. Overall, genetic, environmental and demographic factors seem to influence the observed high levels of variation in early fitness traits among B. pendula populations. Our results suggest limited regeneration capacity for the study species under drier conditions, but further field trials with sufficient replication over environments and seed crops will improve our understanding of its vulnerability to climate change.

Genetic Analysis by nuSSR Markers of Silver Birch (Betula pendula Roth) Populations in Their Southern European Distribution Range

In the main distribution area the genetic pattern of silver birch is dominated by two haplotypes: haplotype A located in the western and north-western Europe, and haplotype C in eastern and southeastern Europe, characterized by high levels of neutral genetic variability within populations, and low differentiation among populations. Information about the amount and structure of genetic variation in the southern marginal areas, representing rear populations left during the expansion of this species from southern glacial refugia, are lacking. The general aim of the study was to investigate the existence of the climatic characteristics typical of the environmental niche of the species, jointly to genetic organization, variation and gene flow, in marginal populations on the Italian Apennines and Greek Southern Rhodope and compare them with populations of the southern part of the main distribution range on the Alps and Balkans. Genetic analysis was performed using nuclear microsatellites loci on 311 trees sampled from 14 populations. Environmental analysis was performed on the multivariate analysis of derived climatic variables. The allelic pattern was analyzed to assess genetic diversity, population diversity and differentiation, population structure and gene flow. The geographic and environmental peripherality did not always match, with some Apennine sites at higher elevation enveloped in the environmental niche. In the peripheral populations on the Apennines, we observed a lower genetic diversity and higher differentiation, with evident genetic barriers detected around these sites. These characteristics were not shown in the marginal Greek populations. Unexpectedly, the southern Italian marginal populations showed genetic links with the Greek and central area of the distribution range. The Greek populations also showed evident gene flow with the Alpine and Balkan areas. The disparity of results in these two marginal areas show that it is not the geographic peripherality or even the ecological marginality that may shape the genetic diversity and structure of marginal populations, but primarily their position as part of the continuous range or as disjunct populations. This outcome suggests different considerations on how to manage their gene pools and the role that these rear populations can play in maintaining the biodiversity of this species.

Adaptive fitness of Sapindus emarginatus Vahl populations towards future climatic regimes and the limiting factors of its distribution

Sapindus emarginatus Vahl (Sapindaceae) also known as ‘Indian Soap nut’ is significantly important for saponin content in its fruits. However, its current population in India is heavily fragmented due to a lack of sustainable harvesting practices. Moreover, changing climatic regimes may further limit its distribution and possibly compromise the survival of the species in nature. The aim of the present study was to: predict the future distribution range of S. emarginatus; identify the bioclimatic variables limiting this distribution and to evaluate its adaptive fitness and genomic resilience towards these variables. To determine future species distribution range and identify limiting bioclimatic variables, we applied two different ecological niche models (ENMs; BioClim and MaxEnt) on real occurrence data (n = 88 locations). The adaptive fitness of the species was evaluated by quantifying the genetic variability with AFLP markers and marker-environmental associations, using AFLP-associated Bayesian statistics. We found 77% overlap between the baseline (2030) and predicted (2100) species distribution ranges, which were primarily determined by maximum temperature (T_MAX) and mean annual precipitation (MAP). The T_MAX and MAP contributed 43.1% and 27.1%, respectively to ENM model prediction. Furthermore, AFLP loci significantly associated with bioclimatic variables, and T_MAX and MAP represent the lowest proportion (6.15%), confirming to the severe response of the species genome towards these variables. Nevertheless, the very low Linkage disequilibrium (LD) in these loci (4.54%) suggests that the current sensitivity to T_MAX and MAP is subject to change during recombination. Moreover, a combination of high heterozygosity (0.40–0.43) and high within-population variability (91.63 ± 0.31%) confirmed high adaptive fitness to maintain reproductive success. Therefore, the current populations of S. emarginatus have substantial genomic resilience towards future climate change, albeit significant conservation efforts (including mass multiplication) are warranted to avoid future deleterious impacts of inbreeding depression on the fragmented populations.

Genotype-by-environment interaction in coast redwood outside natural distribution - search for environmental cues

Background

Effective matching of genotypes and environments is required for the species to reach optimal productivity and act effectively for carbon sequestration. A common garden experiment across five different environments was undertaken to assess genotype x environment interaction (GxE) of coast redwood in order to understand the performance of genotypes across environments.

Results

The quantitative genetic analysis discovered no GxE between investigated environments for diameter at breast height (DBH). However, no genetic component was detected at one environment possibly due to stressful conditions. The implementation of universal response function allowed for the identification of important environmental factors affecting species productivity. Additionally, this approach enabled us to predict the performance of species across the New Zealand environmental conditions.

Conclusions

In combination with quantitative genetic analysis which identified genetically superior material, the URF model can directly identify the optimal geographical regions to maximize productivity. However, the finding of ideally uncorrelated climatic variables for species with narrow ecological amplitude is rather challenging, which complicates construction of informative URF model. This, along with a small number of tested environments, tended to overfit a prediction model which resulted in extreme predictions in untested environments.

Temperature and Rainfall Are Separate Agents of Selection Shaping Population Differentiation in a Forest Tree

Research highlights: We present evidence indicating that covariation of functional traits among populations of a forest tree is not due to genetic constraints, but rather selective covariance arising from local adaptation to different facets of the climate, namely rainfall and temperature. Background and Aims: Traits frequently covary among natural populations. Such covariation can be caused by pleiotropy and/or linkage disequilibrium, but also may arise when the traits are genetically independent as a direct consequence of natural selection, drift, mutation and/or gene flow. Of particular interest are cases of selective covariance, where natural selection directly generates among-population covariance in a set of genetically independent traits. We here studied the causes of population-level covariation in two key traits in the Australian tree Eucalyptus pauciflora. Materials and Methods: We studied covariation in seedling lignotuber size and vegetative juvenility using 37 populations sampled from throughout the geographic and ecological ranges of E. pauciflora on the island of Tasmania. We integrated evidence from multiple sources: (i) comparison of patterns of trait covariation within and among populations; (ii) climate-trait modelling using machine-learning algorithms; and (iii) selection analysis linking trait variation to field growth in an arid environment. Results: We showed strong covariation among populations compared with the weak genetic correlation within populations for the focal traits. Population differentiation in these genetically independent traits was correlated with different home-site climate variables (lignotuber size with temperature; vegetative juvenility with rainfall), which spatially covaried. The role of selection in shaping the population differentiation in lignotuber size was supported by its relationship with fitness measured in the field. Conclusions: Our study highlights the multi-trait nature of adaptation likely to occur as tree species respond to spatial and temporal changes in climate.

Interannual variation in season length is linked to strong co-gradient plasticity of phenology in a montane annual plant

Species are commonly distributed along latitudinal and elevational gradients of growing season length to which they might respond via phenotypic plasticity and/or adaptive genetic differentiation. However, the relative contribution of these processes and whether plasticity, if it occurs, facilitates expansion along season-length gradients remain unclear, but are important for predicting species fates during anthropogenic change. We quantified phenological trait variation in the montane annual Rhinanthus minor for three generations at 12 sites across 900 m of elevation in the Canadian Rocky Mountains and conducted a reciprocal transplant experiment for two generations among nine sites. We compared clines and interannual variation of phenological traits between natural and transplanted individuals. Season length declined by c. 37% along our elevational gradient and, as expected, plants emerged, reached first flower and made their first seed in c. 41% fewer growing degree days under shorter growing seasons. Although reciprocal transplants revealed modest genetic differentiation across elevation, trait clines primarily were due to striking co-gradient plasticity that paralleled genetic differentiation. Co-gradient plasticity likely evolved in response to considerable interannual variation in season length across our elevational transect, and should prepare R. minor to make adaptive changes to phenology in response to ongoing climate change predicted for montane environments.

Modelling of population structure through contemporary groups in genetic evaluation

Background

Forest trees can occupy extensive geography and environmentally highly variable areas which result in high genetic variability in the direction of pressure from natural selection. At the same time, the majority of conifer species are wind-pollinated from both short and long distances, resulting in wide-spread gene flow, which can lead to maladaptation to local conditions. Quantitative analyses of provenance/progeny tests correct for genetic differences between populations to ensure unbiased genetic parameters are obtained. Commonly, the provenance effect is fitted as a fixed term or can be implemented as a contemporary group in the pedigree.

Results

The use of a provenance effect, either as a fixed term or as the same contemporary groups in both maternal and paternal sides of the pedigree, resulted in fairly similar precision of genetic parameters in our case. However, when we developed a phantom contemporary group for the paternal side of the pedigree that considered a different genetic quality of pollen compared with the maternal contribution from trees in the local environment, the model fit and accuracy of breeding values increased.

Conclusion

Consideration of the mating dynamics and the vector of gene flow are important factors in modelling contemporary genetic groups, particularly when implementing pedigrees within a mixed model framework to obtain unbiased estimates of genetic parameters. This approach is especially important in traits involved in local adaptation.

Abiotic genetic adaptation in the Fagaceae

Fagaceae can be found in tropical and temperate regions and contain species of major ecological and economic importance. In times of global climate change, tree populations need to adapt to rapidly changing environmental conditions. The predicted warmer and drier conditions will potentially result in locally maladapted populations. There is evidence that major genera of the Fagaceae are already negatively affected by climate change-related factors such as drought and associated biotic stressors. Therefore, knowledge of the mechanisms underlying adaptation is of great interest. In this review, we summarise current literature related to genetic adaptation to abiotic environmental conditions. We begin with an overview of genetic diversity in Fagaceae species and then summarise current knowledge related to drought stress tolerance, bud burst timing and frost tolerance in the Fagaceae. Finally, we discuss the role of hybridisation, epigenetics and phenotypic plasticity in adaptation.

Mediterranean and Northern Iberian gene pools of wild Castanea sativa Mill. are two differentiated ecotypes originated under natural divergent selection

Nine wild Iberian provenances of Castanea sativa Mill. grouped in two gene pools, North Iberian Peninsula and Mediterranean, were evaluated for several adaptive traits in two provenance-progeny trials with the aim of evaluating the role of natural selection in shaping adaptive variation and increasing our understanding of the genetic structure of this species, as well as reporting complete information on the genetic variation among and within the studied populations. An annual growth rhythm experiment was evaluated during the first 3 years after establishment for phenology, growth, stem form and survival, and a periodic drought-stress experiment was evaluated for dry weight, growth, survival and other related drought traits in both well-watered and drought-stress treatments. The high genetic variability reported in both trials is largely due to the genetic variation among populations. The significant differences reported between quantitative genetic and neutral marker differentiation indicated the local adaptation of these populations through directional selection, mainly for phenology, growth and biomass allocation. A clinal variation among populations was determined through correlations of phenology with latitude and xerothermic index of the provenances, showing that central and southern Mediterranean populations had earlier phenology than northern populations and that drought played a relevant role in this differentiation. The significant correlation between phenological traits and the ancestry values in the Mediterranean gene pool supported the different pattern of behavior between both gene pools and also indicated the existence of two ecotypes: xeric and mesophytic ecotypes, corresponding to Mediterranean and North Iberian gene pools, respectively. The results obtained in the drought-stress experiment confirmed that, in general terms, xeric populations showed a greater adaptability to drought, with more developed root systems and higher survival than northern populations. Moreover, the genetic variability observed within populations indicated the potential response capacity of Iberian C. sativa populations to undergo fast adaptive evolution.

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.

New insights into adaptation and population structure of cork oak using genotyping by sequencing

Species respond to global climatic changes in a local context. Understanding this process, including its speed and intensity, is paramount due to the pace at which such changes are currently occurring. Tree species are particularly interesting to study in this regard due to their long generation times, sedentarism, and ecological and economic importance. Quercus suber L. is an evergreen forest tree species of the Fagaceae family with an essentially Western Mediterranean distribution. Despite frequent assessments of the species' evolutionary history, large-scale genetic studies have mostly relied on plastidial markers, whereas nuclear markers have been used on studies with locally focused sampling strategies. In this work, "Genotyping by sequencing" is used to derive 1,996 single nucleotide polymorphism markers to assess the species' evolutionary history from a nuclear DNA perspective, gain insights into how local adaptation is shaping the species' genetic background, and to forecast how Q. suber may respond to global climatic changes from a genetic perspective. Results reveal (a) an essentially unstructured species, where (b) a balance between gene flow and local adaptation keeps the species' gene pool somewhat homogeneous across its distribution, but still allowing (c) variation clines for the individuals to cope with local conditions. "Risk of Non-Adaptedness" (RONA) analyses suggest that for the considered variables and most sampled locations, (d) the cork oak should not require large shifts in allele frequencies to survive the predicted climatic changes. Future directions include integrating these results with ecological niche modeling perspectives, improving the RONA methodology, and expanding its use to other species. With the implementation presented in this work, the RONA can now also be easily assessed for other organisms.

The extended consequences of genetic conductivity: Mating distance affects community phenotypes in Norway spruce

Anthropogenic landscape-level alterations such as habitat fragmentation and long distance translocation of genetic material are currently altering the genetic connectivity and structure of forest tree populations globally. As the susceptibility of individual trees to dependent organisms is often genetically determined, it is possible that these genetic changes may extend beyond individuals to affect associated communities. To test this, we examined how variation in crossing distance among the progeny of 18 controlled crosses of Norway spruce (Picea abies) populations occurring across central Sweden affected chemical defense, and subsequently, a small community of galling Adelges aphids infecting planted trees at two common garden trails. Although crossing distance did not influence growth, vitality or reproduction in the studied population, it did influence the expression of one candidate defensive chemical compound, apigenin, which was found in higher concentrations within outcrossed trees. We also show that this variation in apigenin induced by crossing distance correlated with susceptibility to one member of the galling community but not the other. Furthermore, the effect of crossing distance on galling communities and the general susceptibility of Norway spruce to infection also varied with environment. Specifically, in the more benign environment, inbred trees suffered greater gall infection than outcrossed trees, which is contrary to general predictions that the effects of inbreeding should be more pronounced in harsher environments. These findings suggest that the effects of variation in crossing distance in forest trees can extend beyond the individual to influence whole communities.

The future of subalpine forests in the Southern Rocky Mountains: Trajectories for Pinus aristata genetic lineages

Like many other high elevation alpine tree species, Rocky Mountain bristlecone pine (Pinus aristata Engelm.) may be particularly vulnerable to climate change. To evaluate its potential vulnerability to shifts in climate, we defined the suitable climate space for each of four genetic lineages of bristlecone pine and for other subalpine tree species in close proximity to bristlecone pine forests. Measuring changes in the suitable climate space for lineage groups is an important step beyond models that assume species are genetically homogenous. The suitable climate space for bristlecone pine in the year 2090 is projected to decline by 74% and the proportional distribution of suitable climate space for genetic lineages shifts toward those associated with warmer and wetter conditions. The 2090 climate space for bristlecone pine exhibits a bimodal distribution along an elevation gradient, presumably due to the persistence of the climate space in the Southern Rocky Mountains and exclusion at mid-elevations by conditions that favor the climate space of other species. These shifts have implications for changes in fire regimes, vulnerability to pest and pathogens, and altered carbon dynamics across the southern Rockies, which may reduce the likelihood of bristlecone pine trees achieving exceptional longevity in the future. The persistence and expansion of climate space for southern bristlecone pine genetic lineage groups in 2090 suggests that these sources may be the least vulnerable in the future. While these lineages may be more likely to persist and therefore present opportunities for proactive management (e.g., assisted migration) to maintain subalpine forest ecosystem services in a warmer world, our findings also imply heighted conservation concern for vulnerable northern lineages facing range contractions.

Phenotypic plasticity vs. local adaptation in quantitative traits differences of Stipa grandis in semi-arid steppe, China

Whether plants are able to adapt to environmental changes depends on their genetic characteristics and phenotypic plastic responses. We investigated the phenotypic responses of 7 populations of an important dominant species in semi-arid steppe of China - Stipa grandis, and then distinguished which adaptive mechanism(s), phenotypic plasticity or local adaptation, was/were involved in this species to adapt to environmental changes. (1) All traits were significantly influenced by the interaction of population and growth condition and by population in each condition, and inter-population variability (CVinter) was larger in the field than in the common garden for 8/9 traits, indicating that both phenotypic plasticity and genetic differentiation controlled the phenotypic differences of S. grandis. (2) From a functional standpoint, the significant relationships between the values of traits in the common garden and the environmental variables in their original habitats couldn't support local habitat adaptation of these traits. (3) Low CVintra, low quantitative differentiation among populations (Q ST ), and low plasticity shown in the western populations indicated the very low adaptive potential of S. grandis to environmental changes. (4) From the original habitats to the common garden which is far away from S. grandis distribution region, positive phenotypic responses were found in several populations, indicating that some original habitats have become unfavorable for S. grandis.

Evolutionary dynamics of the leaf phenological cycle in an oak metapopulation along an elevation gradient

It has been predicted that environmental changes will radically alter the selective pressures on phenological traits. Long-lived species, such as trees, will be particularly affected, as they may need to undergo major adaptive change over only one or a few generations. The traits describing the annual life cycle of trees are generally highly evolvable, but nothing is known about the strength of their genetic correlations. Tight correlations can impose strong evolutionary constraints, potentially hampering the adaptation of multivariate phenological phenotypes. In this study, we investigated the evolutionary, genetic and environmental components of the timing of leaf unfolding and senescence within an oak metapopulation along an elevation gradient. Population divergence, estimated from in situ and common-garden data, was compared to expectations under neutral evolution, based on microsatellite markers. This approach made it possible (1) to evaluate the influence of genetic correlation on multivariate local adaptation to elevation and (2) to identify traits probably exposed to past selective pressures due to the colder climate at high elevation. The genetic correlation was positive but very weak, indicating that genetic constraints did not shape the local adaptation pattern for leaf phenology. Both spring and fall (leaf unfolding and senescence, respectively) phenology timings were involved in local adaptation, but leaf unfolding was probably the trait most exposed to climate change-induced selection. Our data indicated that genetic variation makes a much smaller contribution to adaptation than the considerable plastic variation displayed by a tree during its lifetime. The evolutionary potential of leaf phenology is, therefore, probably not the most critical aspect for short-term population survival in a changing climate.

Estimating the ability of plants to plastically track temperature-mediated shifts in the spring phenological optimum

One consequence of rising spring temperatures is that the optimum timing of key life-history events may advance. Where this is the case, a population's fate may depend on the degree to which it is able to track a change in the optimum timing either via plasticity or via adaptation. Estimating the effect that temperature change will have on optimum timing using standard approaches is logistically challenging, with the result that very few estimates of this important parameter exist. Here we adopt an alternative statistical method that substitutes space for time to estimate the temperature sensitivity of the optimum timing of 22 plant species based on >200 000 spatiotemporal phenological observations from across the United Kingdom. We find that first leafing and flowering dates are sensitive to forcing (spring) temperatures, with optimum timing advancing by an average of 3 days °C^-1 and plastic responses to forcing between -3 and -8 days °C^-1 . Chilling (autumn/winter) temperatures and photoperiod tend to be important cues for species with early and late phenology, respectively. For most species, we find that plasticity is adaptive, and for seven species, plasticity is sufficient to track geographic variation in the optimum phenology. For four species, we find that plasticity is significantly steeper than the optimum slope that we estimate between forcing temperature and phenology, and we examine possible explanations for this countergradient pattern, including local adaptation.

Trait variation along elevation gradients in a dominant woody shrub is population-specific and driven by plasticity

Elevation gradients are frequently used as space-for-time substitutions to infer species' trait responses to climate change. However, studies rarely investigate whether trait responses to elevation are widespread or population-specific within a species, and the relative genetic and plastic contributions to such trait responses may not be well understood. Here, we examine plant trait variation in the dominant woody shrub, Rhododendron maximum, along elevation gradients in three populations in the South Central Appalachian Mountains, USA, in both field and common garden environments. We ask the following: (i) do plant traits vary along elevation? (ii) do trait responses to elevation differ across populations, and if so, why? and (iii) does genetic differentiation or phenotypic plasticity drive trait variation within and among populations? We found that internode length, shoot length, leaf dry mass, and leaf area varied along elevation, but that these responses were generally unique to one population, suggesting that trait responses to environmental gradients are population-specific. A common garden experiment identified no genetic basis to variation along elevation or among populations in any trait, suggesting that plasticity drives local and regional trait variation and may play a key role in the persistence of plant species such as R. maximum with contemporary climate change. Overall, our findings highlight the importance of examining multiple locations in future elevation studies and indicate that, for a given plant species, the magnitude of trait responses to global climate change may vary by location.

Adaptive and plastic responses of Quercus petraea populations to climate across Europe

How temperate forests will respond to climate change is uncertain; projections range from severe decline to increased growth. We conducted field tests of sessile oak (Quercus petraea), a widespread keystone European forest tree species, including more than 150 000 trees sourced from 116 geographically diverse populations. The tests were planted on 23 field sites in six European countries, in order to expose them to a wide range of climates, including sites reflecting future warmer and drier climates. By assessing tree height and survival, our objectives were twofold: (i) to identify the source of differential population responses to climate (genetic differentiation due to past divergent climatic selection vs. plastic responses to ongoing climate change) and (ii) to explore which climatic variables (temperature or precipitation) trigger the population responses. Tree growth and survival were modeled for contemporary climate and then projected using data from four regional climate models for years 2071-2100, using two greenhouse gas concentration trajectory scenarios each. Overall, results indicated a moderate response of tree height and survival to climate variation, with changes in dryness (either annual or during the growing season) explaining the major part of the response. While, on average, populations exhibited local adaptation, there was significant clinal population differentiation for height growth with winter temperature at the site of origin. The most moderate climate model (HIRHAM5-EC; rcp4.5) predicted minor decreases in height and survival, while the most extreme model (CCLM4-GEM2-ES; rcp8.5) predicted large decreases in survival and growth for southern and southeastern edge populations (Hungary and Turkey). Other nonmarginal populations with continental climates were predicted to be severely and negatively affected (Bercé, France), while populations at the contemporary northern limit (colder and humid maritime regions; Denmark and Norway) will probably not show large changes in growth and survival in response to climate change.

Genetic structure and diversity of the selfing model grass Brachypodium stacei (Poaceae) in Western Mediterranean: out of the Iberian Peninsula and into the islands

Annual Mediterranean species of the genus Brachypodium are promising model plants for energy crops since their selfing nature and short-life cycles are an advantage in breeding programs. The false brome, B. distachyon, has already been sequenced and new genomic initiatives have triggered the de-novo genome sequencing of its close relatives such as B. stacei, a species that was until recently mistaken for B. distachyon. However, the success of these initiatives hinges on detailed knowledge about the distribution of genetic variation within and among populations for the effective use of germplasm in a breeding program. Understanding population genetic diversity and genetic structure is also an important prerequisite for designing effective experimental populations for genomic wide studies. However, population genetic data are still limited in B. stacei. We therefore selected and amplified 10 nuclear microsatellite markers to depict patterns of population structure and genetic variation among 181 individuals from 19 populations of B. stacei occurring in its predominant range, the western Mediterranean area: mainland Iberian Peninsula, continental Balearic Islands and oceanic Canary Islands. Our genetic results support the occurrence of a predominant selfing system with extremely high levels of homozygosity across the analyzed populations. Despite the low level of genetic variation found, two different genetic clusters were retrieved, one clustering all SE Iberian mainland populations and the island of Minorca and another one grouping all S Iberian mainland populations, the Canary Islands and all Majorcan populations except one that clustered with the former group. These results, together with a high sharing of alleles (89%) suggest different colonization routes from the mainland Iberian Peninsula into the islands. A recent colonization scenario could explain the relatively low levels of genetic diversity and low number of alleles found in the Canary Islands populations while older colonization events are hypothesized to explain the high genetic diversity values found in the Majorcan populations. Our study provides widely applicable information about geographical patterns of genetic variation in B. stacei. Among others, the genetic pattern and the existence of local alleles will need to be adequately reflected in the germplasm collection of B. stacei for efficient genome wide association studies.

Local Adaptation in European Firs Assessed through Extensive Sampling across Altitudinal Gradients in Southern Europe

BACKGROUND

Local adaptation is a key driver of phenotypic and genetic divergence at loci responsible for adaptive traits variations in forest tree populations. Its experimental assessment requires rigorous sampling strategies such as those involving population pairs replicated across broad spatial scales.

METHODS

A hierarchical Bayesian model of selection (HBM) that explicitly considers both the replication of the environmental contrast and the hierarchical genetic structure among replicated study sites is introduced. Its power was assessed through simulations and compared to classical 'within-site' approaches (FDIST, BAYESCAN) and a simplified, within-site, version of the model introduced here (SBM).

RESULTS

HBM demonstrates that hierarchical approaches are very powerful to detect replicated patterns of adaptive divergence with low false-discovery (FDR) and false-non-discovery (FNR) rates compared to the analysis of different sites separately through within-site approaches. The hypothesis of local adaptation to altitude was further addressed by analyzing replicated Abies alba population pairs (low and high elevations) across the species' southern distribution range, where the effects of climatic selection are expected to be the strongest. For comparison, a single population pair from the closely related species A. cephalonica was also analyzed. The hierarchical model did not detect any pattern of adaptive divergence to altitude replicated in the different study sites. Instead, idiosyncratic patterns of local adaptation among sites were detected by within-site approaches.

CONCLUSION

Hierarchical approaches may miss idiosyncratic patterns of adaptation among sites, and we strongly recommend the use of both hierarchical (multi-site) and classical (within-site) approaches when addressing the question of adaptation across broad spatial scales.

Intraspecific Variation in Wood Anatomical, Hydraulic, and Foliar Traits in Ten European Beech Provenances Differing in Growth Yield

In angiosperms, many studies have described the inter-specific variability of hydraulic-related traits and little is known at the intra-specific level. This information is however mandatory to assess the adaptive capacities of tree populations in the context of increasing drought frequency and severity. Ten 20-year old European beech (Fagus sylvatica L.) provenances representing the entire distribution range throughout Europe and differing significantly in aboveground biomass increment (ABI) by a factor of up to four were investigated for branch wood anatomical, hydraulic, and foliar traits in a provenance trial located in Northern Europe. We quantified to which extend xylem hydraulic and leaf traits are under genetic control and tested whether the xylem hydraulic properties (hydraulic efficiency and safety) trades off with yield and wood anatomical and leaf traits. Our results showed that only three out of 22 investigated ecophysiological traits showed significant genetic differentiations between provenances, namely vessel density (VD), the xylem pressure causing 88% loss of hydraulic conductance and mean leaf size. Depending of the ecophysiological traits measured, genetic differentiation between populations explained 0-14% of total phenotypic variation, while intra-population variability was higher than inter-population variability. Most wood anatomical traits and some foliar traits were additionally related to the climate of provenance origin. The lumen to sapwood area ratio, vessel diameter, theoretical specific conductivity and theoretical leaf-specific conductivity as well as the C:N-ratio increased with climatic aridity at the place of origin while the carbon isotope signature (δ(13)C) decreased. Contrary to our assumption, none of the wood anatomical traits were related to embolism resistance but were strong determinants of hydraulic efficiency. Although ABI was associated with both VD and δ(13)C, both hydraulic efficiency and embolism resistance were unrelated, disproving the assumed trade-off between hydraulic efficiency and safety. European beech seems to compensate increasing water stress with growing size mainly by adjusting vessel number and not vessel diameter. In conclusion, European beech has a high potential capacity to cope with climate change due to the high degree of intra-population genetic variability.