Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood.

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

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.

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.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits.

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Oaks: an evolutionary success story.

The genus Quercus is among the most widespread and species-rich tree genera in the northern hemisphere. The extraordinary species diversity in America and Asia together with the continuous continental distribution of a limited number of European species raise questions about how macro- and microevolutionary processes made the genus Quercus an evolutionary success. Synthesizing conclusions reached during the past three decades by complementary approaches in phylogenetics, phylogeography, genomics, ecology, paleobotany, population biology and quantitative genetics, this review aims to illuminate evolutionary processes leading to the radiation and expansion of oaks. From opposing scales of time and geography, we converge on four overarching explanations of evolutionary success in oaks: accumulation of large reservoirs of diversity within populations and species; ability for rapid migration contributing to ecological priority effects on lineage diversification; high rates of evolutionary divergence within clades combined with convergent solutions to ecological problems across clades; and propensity for hybridization, contributing to adaptive introgression and facilitating migration. Finally, we explore potential future research avenues, emphasizing the integration of microevolutionary and macroevolutionary perspectives.

Oak symbolism in the light of genomics.

Throughout the Northern Hemisphere, human societies, political systems, and religions have appropriated oaks in symbolic representations. In this review, we explore the possible associations between recent genetic and genomic findings and the symbolic representations of oaks. We first consider the ways in which evolutionary history during the Holocene has tightened links between humans and oaks in Europe, and how this may have led to symbolic representations. We then show how recent findings concerning the structure and evolution of the oak genome have provided additional knowledge about symbolic representations, such as longevity, cohesiveness, and robustness.

Adaptive introgression as a driver of local adaptation to climate in European white oaks.

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

Massive postglacial gene flow between European white oaks uncovered genes underlying species barriers.

Oaks are dominant forest tree species widely distributed across the Northern Hemisphere, where they constitute natural resources of economic, ecological, social and historical value. Hybridisation and adaptive introgression have long been thought to be major drivers of their ecological success. Therefore, the maintenance of species barriers remains a key question, given the extent of interspecific gene flow. In this study, we made use of the tremendous genetic variation among four European white oak species (31 million single nucleotide polymorphisms (SNPs)) to infer the evolutionary history of these species, study patterns of genetic differentiation and identify reproductive barriers. We first analysed the ecological and historical relationships among these species and inferred a long-term strict isolation followed by a recent and extensive postglacial contact using approximate Bayesian computation. Assuming this demographic scenario, we then performed backward simulations to generate the expected distributions of differentiation under neutrality to scan their genomes for reproductive barriers. We finally identified important intrinsic and ecological functions driving the reproductive isolation. We discussed the importance of identifying the genetic basis for the ecological preferences between these oak species and its implications for the renewal of European forests under global warming.

Genomic landscape of the global oak phylogeny.

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

Range margin populations show high climate adaptation lags in European trees.

How populations of long-living species respond to climate change depends on phenotypic plasticity and local adaptation processes. Marginal populations are expected to have lags in adaptation (i.e. differences between the climatic optimum that maximizes population fitness and the local climate) because they receive pre-adapted alleles from core populations preventing them from reaching a local optimum in their climatically marginal habitat. Yet, whether adaptation lags in marginal populations are a common feature across phylogenetically and ecologically different species and how lags can change with climate change remain unexplored. To test for range-wide patterns of phenotypic variation and adaptation lags of populations to climate, we (a) built model ensembles of tree height accounting for the climate of population origin and the climate of the site for 706 populations monitored in 97 common garden experiments covering the range of six European forest tree species; (b) estimated populations' adaptation lags as the differences between the climatic optimum that maximizes tree height and the climate of the origin of each population; (c) identified adaptation lag patterns for populations coming from the warm/dry and cold/wet margins and from the distribution core of each species range. We found that (a) phenotypic variation is driven by either temperature or precipitation; (b) adaptation lags are consistently higher in climatic margin populations (cold/warm, dry/wet) than in core populations; (c) predictions for future warmer climates suggest adaptation lags would decrease in cold margin populations, slightly increasing tree height, while adaptation lags would increase in core and warm margin populations, sharply decreasing tree height. Our results suggest that warm margin populations are the most vulnerable to climate change, but understanding how these populations can cope with future climates depend on whether other fitness-related traits could show similar adaptation lag patterns.

High-Throughput DNA sequencing of ancient wood.

Reconstructing the colonization and demographic dynamics that gave rise to extant forests is essential to forecasts of forest responses to environmental changes. Classical approaches to map how population of trees changed through space and time largely rely on pollen distribution patterns, with only a limited number of studies exploiting DNA molecules preserved in wooden tree archaeological and subfossil remains. Here, we advance such analyses by applying high-throughput (HTS) DNA sequencing to wood archaeological and subfossil material for the first time, using a comprehensive sample of 167 European white oak waterlogged remains spanning a large temporal (from 550 to 9,800 years) and geographical range across Europe. The successful characterization of the endogenous DNA and exogenous microbial DNA of 140 (~83%) samples helped the identification of environmental conditions favouring long-term DNA preservation in wood remains, and started to unveil the first trends in the DNA decay process in wood material. Additionally, the maternally inherited chloroplast haplotypes of 21 samples from three periods of forest human-induced use (Neolithic, Bronze Age and Middle Ages) were found to be consistent with those of modern populations growing in the same geographic areas. Our work paves the way for further studies aiming at using ancient DNA preserved in wood to reconstruct the micro-evolutionary response of trees to climate change and human forest management.

Assessing inter- and intraspecific variability of xylem vulnerability to embolism in oaks.

The genus Quercus comprises important species in forestry not only for their productive value but also for their ability to withstand drought. Hence an evaluation of inter- and intraspecific variation in drought tolerance is important for selecting the best adapted species and provenances for future afforestation. The presence of long vessels makes it difficult to assess xylem vulnerability to embolism in oak. Thanks to the development of an in situ flow centrifuge equipped with a large rotor, we quantified (i) the between species variability of embolism resistance in four native and two exotic species of oaks in Europe and (ii) the within species variability in Quercus petraea. Embolism resistance varied significantly among species, with the pressure inducing 50% loss of hydraulic conductivity (P50 ) ranging between - 7.0 and -4.2 MPa. Species native to the Mediterranean region were more resistant than pan-European species. In contrast, intraspecific variability in embolism resistance in Q. petraea was low within provenances and null among provenances. A positive correlation between P50 and vessel diameter among the six oak species indicates that the more embolism resistant species had narrower xylem vessels. However, this tradeoff between hydraulic efficiency and safety was not observed between Q. petraea provenances.

Fine-scale species distribution changes in a mixed oak stand over two successive generations.

Large-scale tree distribution changes have received considerable attention but underlying demo-genetic mechanisms are less well documented. We used a diachronic approach to track species shifts in a mixed oak stand (Quercus petraea-Quercus robur) at a fine spatiotemporal scale. Species assignment was made using single nucleotide polymorphism (SNP) fingerprints employing clustering and parentage analysis. Mating patterns and reproductive success were assessed by parentage analysis. Plot-based inventories of soil parameters and sapling densities provided ecological and demographic information, respectively. Sapling density and reproductive success was higher in Q. petraea than in Q. robur, and were correlated with a spatial expansion of Q. petraea (50% to 67% of the area). Admixed trees resulting from hybridization and backcrossing between the two species were more frequent under the Q. robur canopy. We suspect that species' differential responses to ongoing environmental changes and interspecific competition are the predominant factors accounting for the recruitment success of Q. petraea, while human interference, differential reproduction and hybridization (and backcrossings) are probably of more limited importance. We anticipate in mixed Q. petraea-Q. robur stands, under current ongoing environmental change, that these processes will be enhanced, at least in the western part of the distribution of the two species.

Extensive recent secondary contacts between four European white oak species.

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.

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.

Back to America: tracking the origin of European introduced populations of Quercus rubra L.

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.

Forest tree genomics: growing resources and applications.

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

Signatures of local adaptation in candidate genes of oaks (Quercus spp.) with respect to present and future climatic conditions.

Testing how populations are locally adapted and predicting their response to their future environment is of key importance in view of climate change. Landscape genomics is a powerful approach to investigate genes and environmental factors involved in local adaptation. In a pooled amplicon sequencing approach of 94 genes in 71 populations, we tested whether >3500 single nucleotide polymorphisms (SNPs) in the three most common oak species in Switzerland (Quercus petraea, Q. pubescens, Q. robur) show an association with abiotic factors related to local topography, historical climate and soil characteristics. In the analysis including all species, the most frequently associated environmental factors were those best describing the habitats of the species. In the species-specific analyses, the most important environmental factors and associated SNPs greatly differed among species. However, we identified one SNP and seven genes that were associated with the same environmental factor across all species. We finally used regressions of allele frequencies of the most strongly associated SNPs along environmental gradients to predict the risk of nonadaptedness (RONA), which represents the average change in allele frequency at climate-associated loci theoretically required to match future climatic conditions. RONA is considerable for some populations and species (up to 48% in single populations) and strongly differs among species. Given the long generation time of oaks, some of the required allele frequency changes might not be realistic to achieve based on standing genetic variation. Hence, future adaptedness requires gene flow or planting of individuals carrying beneficial alleles from habitats currently matching future climatic conditions.

The oak gene expression atlas: insights into Fagaceae genome evolution and the discovery of genes regulated during bud dormancy release.

BACKGROUND: Many northern-hemisphere forests are dominated by oaks. These species extend over diverse environmental conditions and are thus interesting models for studies of plant adaptation and speciation. The genomic toolbox is an important asset for exploring the functional variation associated with natural selection. RESULTS: The assembly of previously available and newly developed long and short sequence reads for two sympatric oak species, Quercus robur and Quercus petraea, generated a comprehensive catalog of transcripts for oak. The functional annotation of 91 k contigs demonstrated the presence of a large proportion of plant genes in this unigene set. Comparisons with SwissProt accessions and five plant gene models revealed orthologous relationships, making it possible to decipher the evolution of the oak genome. In particular, it was possible to align 9.5 thousand oak coding sequences with the equivalent sequences on peach chromosomes. Finally, RNA-seq data shed new light on the gene networks underlying vegetative bud dormancy release, a key stage in development allowing plants to adapt their phenology to the environment. CONCLUSION: In addition to providing a vast array of expressed genes, this study generated essential information about oak genome evolution and the regulation of genes associated with vegetative bud phenology, an important adaptive traits in trees. This resource contributes to the annotation of the oak genome sequence and will provide support for forward genetics approaches aiming to link genotypes with adaptive phenotypes.

Chilling and heat requirements for leaf unfolding in European beech and sessile oak populations at the southern limit of their distribution range.

With global warming, an advance in spring leaf phenology has been reported worldwide. However, it is difficult to forecast phenology for a given species, due to a lack of knowledge about chilling requirements. We quantified chilling and heat requirements for leaf unfolding in two European tree species and investigated their relative contributions to phenological variations between and within populations. We used an extensive database containing information about the leaf phenology of 14 oak and 10 beech populations monitored over elevation gradients since 2005. In parallel, we studied the various bud dormancy phases, in controlled conditions, by regularly sampling low- and high-elevation populations during fall and winter. Oak was 2.3 times more sensitive to temperature for leaf unfolding over the elevation gradient and had a lower chilling requirement for dormancy release than beech. We found that chilling is currently insufficient for the full release of dormancy, for both species, at the lowest elevations in the area studied. Genetic variation in leaf unfolding timing between and within oak populations was probably due to differences in heat requirement rather than differences in chilling requirement. Our results demonstrate the importance of chilling for leaf unfolding in forest trees and indicate that the advance in leaf unfolding phenology with increasing temperature will probably be less pronounced than forecasted. This highlights the urgent need to determine experimentally the interactions between chilling and heat requirements in forest tree species, to improve our understanding and modeling of changes in phenological timing under global warming.