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

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.

Phylogenomic analyses highlight innovation and introgression in the continental radiations of Fagaceae across the Northern Hemisphere

Northern Hemisphere forests changed drastically in the early Eocene with the diversification of the oak family (Fagaceae). Cooling climates over the next 20 million years fostered the spread of temperate biomes that became increasingly dominated by oaks and their chestnut relatives. Here we use phylogenomic analyses of nuclear and plastid genomes to investigate the timing and pattern of major macroevolutionary events and ancient genome-wide signatures of hybridization across Fagaceae. Innovation related to seed dispersal is implicated in triggering waves of continental radiations beginning with the rapid diversification of major lineages and resulting in unparalleled transformation of forest dynamics within 15 million years following the K-Pg extinction. We detect introgression at multiple time scales, including ancient events predating the origination of genus-level diversity. As oak lineages moved into newly available temperate habitats in the early Miocene, secondary contact between previously isolated species occurred. This resulted in adaptive introgression, which may have further amplified the diversification of white oaks across Eurasia.

Fagaceae are diverse family including trees of ecological and economic importance. This phylogenomic analysis of nuclear and plastid genomes reconstructs evolutionary history and finds evidence of multiple adaptive introgression events in this important plant family.

Small-scale genetic structure and mating patterns in an extensive sessile oak forest (Quercus petraea (Matt.) Liebl.)

Oaks (Quercus) are major components of temperate forest ecosystems in the Northern Hemisphere where they form intermediate or climax communities. Sessile oak (Quercus petraea) forests represent the climax vegetation in eastern Germany and western Poland. Here, sessile oak forms pure stands or occurs intermixed with Scots Pine (Pinus sylvestris). A large body of research is available on gene flow, reproduction dynamics, and genetic structure in fragmented landscapes and mixed populations. At the same time, our knowledge regarding large, contiguous, and monospecific populations is considerably less well developed. Our study is an attempt to further develop our understanding of the reproduction ecology of sessile oak as an ecologically and economically important forest tree by analyzing mating patterns and genetic structure within adult trees and seedlings originating from one or two reproduction events in an extensive, naturally regenerating sessile oak forest. We detected positive spatial genetic structure up to 30 meters between adult trees and up to 40 meters between seedlings. Seed dispersal distances averaged 8.4 meters. Pollen dispersal distances averaged 22.6 meters. In both cases, the largest proportion of the dispersal occurred over short distances. Dispersal over longer distances was more common for pollen but also appeared regularly for seeds. The reproductive success of individual trees was highly skewed. Only 41 percent of all adult trees produced any offspring while the majority did not participate in reproduction. Among those trees that contributed to the analyzed seedling sample, 80 percent contributed 1-3 gametes. Only 20 percent of all parent trees contributed four or more gametes. However, these relatively few most fertile trees contributed 51 percent of all gametes within the seedling sample. Vitality and growth differed significantly between reproducing and nonreproducing adult trees with reproducing trees being more vital and vigorous than nonreproducing individuals. Our study demonstrates that extensive, apparently homogenous oak forests are far from uniform on the genetic level. On the contrary, they form highly complex mosaics of remarkably small local neighborhoods. This counterbalances the levelling effect of long-distance dispersal and may increase the species' adaptive potential. Incorporating these dynamics in the management, conservation, and restoration of oak forests can support the conservation of forest genetic diversity and assist those forests in coping with environmental change.

Molecular Research on Stress Responses in Quercus spp.: From Classical Biochemistry to Systems Biology through Omics Analysis

The genus Quercus (oak), family Fagaceae, comprises around 500 species, being one of the most important and dominant woody angiosperms in the Northern Hemisphere. Nowadays, it is threatened by environmental cues, which are either of biotic or abiotic origin. This causes tree decline, dieback, and deforestation, which can worsen in a climate change scenario. In the 21st century, biotechnology should take a pivotal role in facing this problem and proposing sustainable management and conservation strategies for forests. As a non-domesticated, long-lived species, the only plausible approach for tree breeding is exploiting the natural diversity present in this species and the selection of elite, more resilient genotypes, based on molecular markers. In this direction, it is important to investigate the molecular mechanisms of the tolerance or resistance to stresses, and the identification of genes, gene products, and metabolites related to this phenotype. This research is being performed by using classical biochemistry or the most recent omics (genomics, epigenomics, transcriptomics, proteomics, and metabolomics) approaches, which should be integrated with other physiological and morphological techniques in the Systems Biology direction. This review is focused on the current state-of-the-art of such approaches for describing and integrating the latest knowledge on biotic and abiotic stress responses in Quercus spp., with special reference to Quercus ilex, the system on which the authors have been working for the last 15 years. While biotic stress factors mainly include fungi and insects such as Phytophthora cinnamomi, Cerambyx welensii, and Operophtera brumata, abiotic stress factors include salinity, drought, waterlogging, soil pollutants, cold, heat, carbon dioxide, ozone, and ultraviolet radiation. The review is structured following the Central Dogma of Molecular Biology and the omic cascade, from DNA (genomics, epigenomics, and DNA-based markers) to metabolites (metabolomics), through mRNA (transcriptomics) and proteins (proteomics). An integrated view of the different approaches, challenges, and future directions is critically discussed.

How does contemporary selection shape oak phenotypes?

Most existing forests are subjected to natural and human-mediated selection pressures, which have increased due to climate change and the increasing needs of human societies for wood, fibre and fuel resources. It remains largely unknown how these pressures trigger evolutionary changes. We address this issue here for temperate European oaks (Quercus petraea and Q. robur), which grow in mixed stands, under even-aged management regimes. We screened numerous functional traits for univariate selection gradients and for expected and observed genetic changes over two successive generations. In both species, growth, leaf morphology and physiology, and defence-related traits displayed significant selection gradients and predicted shifts, whereas phenology, water metabolism, structure and resilience-related traits did not. However, the direction of the selection response and the potential for adaptive evolution differed between the two species. Quercus petraea had a much larger phenotypic and genetic variance of fitness than Q. robur. This difference raises concerns about the adaptive response of Q. robur to contemporary selection pressures. Our investigations suggest that Q. robur will probably decline steadily, particularly in mixed stands with Q. petraea, consistent with the contrasting demographic dynamics of the two species.

Climate Change Synchronizes Growth and iWUE Across Species in a Temperate-Submediterranean Mixed Oak Forest

Tree species have good tolerance to a range of environmental conditions, though their ability to respond and persist to environmental changes is dramatically reduced at the rear-edge distribution limits. At those edges, gene flow conferring adaptation is impaired due to lack of populations at lower latitudes. Thus, trees mainly rely on phenotypic changes to buffer against long-term environmental changes. Interspecific hybridization may offer an alternative mechanism in the generation of novel genetic recombinants that could be particularly valuable to ensure persistence in geographically isolated forests. In this paper, we take advantage of the longevity of a temperate-submediterranean mixed-oak forest to explore the long-term impact of environmental changes on two different oak species and their hybrid. Individual trees were genetically characterized and classified into three groups: pure Quercus petraea (Matt.), Liebl, pure Q. pyrenaica Willd, and hybrids. We calculated basal area increment and intrinsic water-use efficiency (iWUE) from tree-ring width and δ13C per genetic group, respectively. Tree-growth drivers were assessed using correlation analyses and generalized linear mixed models for two contrasting climatic periods: (1880-1915, colder with [CO2] < 303 ppm; and 1980-2015, warmer with [CO2] > 338 ppm). The three genetic groups have increased radial growth and iWUE during the last decades, being the least drought-tolerant QuPe the most sensitive species to water stress. However, no significant differences were found among genetic groups neither in mean growth rate nor in mean iWUE. Furthermore, little differences were found in the response to climate among groups. Genetic groups only differed in the relationship between δ13C and temperature and precipitation during the earlier period, but such a difference disappeared during the recent decades. Climate change may have promoted species-level convergence as a response to environment-induced growth limitations, which translated in synchronized growth and response to climate as well as a tighter stomatal control and increased iWUE across coexisting oak species.

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.

A species-discriminatory single-nucleotide polymorphism set reveals maintenance of species integrity in hybridizing European white oaks (Quercus spp.) despite high levels of admixture

BACKGROUND AND AIMS

Hybridization and introgression play an important role in the evolution and diversification of plants. To assess the degree of past and current hybridization, the level of genetic admixture in populations needs to be investigated. Ongoing hybridization and blurred species separation have made it challenging to assign European white oak taxa based on leaf morphology and/or genetic markers and to assess the level of admixture. Therefore, there is a need for powerful markers that differentiate between taxa. Here, we established a condensed set of single-nucleotide polymorphism (SNP) markers to reliably differentiate between the three most common oak species in temperate European forests (Quercus robur, Q. petraea, Q. pubescens) and to assess the degree of admixture in a large set of selected Swiss populations.

METHODS

A training set of 194 presumably pure reference samples from Switzerland and Europe was used to assign 633 test individuals with two different approaches (population genetic-based/Bayesian vs. assumption-free/discriminative classifier) using 58 selected SNPs from coding regions. Admixture was calculated at the individual and population level with the Shannon diversity index based on individual assignment probabilities.

KEY RESULTS

Depending on the approach, 97.5-100 % of training individuals were assigned correctly, and additional analyses showed that the established SNP set could be further reduced while maintaining its discriminatory power. The two assignment approaches showed high overlap (99 %) in assigning training individuals and slightly less overlap in test individuals (84 %). Levels of admixture varied widely among populations. Mixed stands of Q. petraea and Q. pubescens revealed much higher degrees of admixture than mixed stands of the other two taxon pairs, accentuating high levels of gene flow between these two taxa in Switzerland.

CONCLUSIONS

Our set of SNPs warrants reliable taxon discrimination with great potential for further applications. We show that the three European white oak taxa have largely retained their species integrity in Switzerland despite high levels of admixture.

In situ estimation of genetic variation of functional and ecological traits in Quercus petraea and Q.robur

BACKGROUND

Predicting the evolutionary potential of natural tree populations requires the estimation of heritability and genetic correlations among traits on which selection acts, as differences in evolutionary success between species may rely on differences for these genetic parameters. In situ estimates are expected to be more accurate than measures done under controlled conditions which do not reflect the natural environmental variance.

AIMS

The aim of the current study was to estimate three genetic parameters (i.e. heritability, evolvability and genetic correlations) in a natural mixed oak stand composed of Quercus petraea and Quercus robur about 100 years old, for 58 traits of ecological and functional relevance (growth, reproduction, phenology, physiology, resilience, structure, morphology and defence).

METHODS

First we estimated genetic parameters directly in situ using realized genomic relatedness of adult trees and parentage relationships over two generations to estimate the traits additive variance. Secondly, we benefited from existing ex situ experiments (progeny tests and conservation collection) installed with the same populations, thus allowing comparisons of in situ heritability estimates with more traditional methods.

RESULTS

Heritability and evolvability estimates obtained with different methods varied substantially and showed large confidence intervals, however we found that in situ were less precise than ex situ estimates, and assessments over two generations (with deeper relatedness) improved estimates of heritability while large sampling sizes are needed for accurate estimations. At the biological level, heritability values varied moderately across different ecological and functional categories of traits, and genetic correlations among traits were conserved over the two species.

CONCLUSION

We identified limits for using realized genomic relatedness in natural stands to estimate the genetic variance, given the overall low variance of genetic relatedness and the rather low sampling sizes of currently used long term genetic plots in forestry. These limits can be overcome if larger sample sizes are considered, or if the approach is extended over the next generation.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco‐evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large‐scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco‐evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

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.

Common garden comparisons confirm inherited differences in sensitivity to climate change between forest tree species

The natural distribution, habitat, growth and evolutionary history of tree species are strongly dependent on ecological and genetic processes in ecosystems subject to fluctuating climatic conditions, but there have been few experimental comparisons of sensitivity between species. We compared the responses of two broadleaved tree species (Fagus sylvatica and Quercus petraea) and two conifer tree species (Pinus sylvestris and Picea abies) to climatic transfers by fitting models containing the same climatic variables. We used published data from European provenance test networks to model the responses of individual populations nested within species. A mixed model approach was applied to develop a response function for tree height over climatic transfer distance, taking into account the climatic conditions at both the seed source and the test location. The two broadleaved species had flat climatic response curves, indicating high levels of plasticity in populations, facilitating adaptation to a broader range of environments, and conferring a high potential for resilience in the face of climatic change. By contrast, the two conifer species had response curves with more pronounced slopes, indicating a lower resilience to climate change. This finding may reflect stronger genetic clines in P. sylvestris and P. abies, which constrain their climate responses to narrower climatic ranges. The response functions had maxima that deviated from the expected maximum productivity in the climate of provenance towards cooler/moister climate conditions, which we interpreted as an adaptation lag. Unilateral, linear regression analyses following transfer to warmer and drier sites confirmed a decline in productivity, predictive of the likely impact of ongoing climate change on forest populations. The responses to mimicked climate change evaluated here are of considerable interest for forestry and ecology, supporting projections of expected performance based on “real-time” field data.

Development of Target Sequence Capture and Estimation of Genomic Relatedness in a Mixed Oak Stand

Anticipating the evolutionary responses of long-lived organisms, such as trees, to environmental changes, requires the assessment of genetic variation of adaptive traits in natural populations. To this end, high-density markers are needed to calculate genomic relatedness between individuals allowing to estimate the genetic variance of traits in wild populations. We designed a targeted capture-based, next-generation sequencing assay based on the highly heterozygous pedunculate oak (Quercus robur) reference genome, for the sequencing of 3 Mb of genic and intergenic regions. Using a mixed stand of 293 Q. robur and Q. petraea genotypes we successfully captured over 97% of the target sequences, corresponding to 0.39% of the oak genome, with sufficient depth (97×) for the detection of about 190,000 SNPs evenly spread over the targeted regions. We validated the technique by evaluating its reproducibility, and comparing the genomic relatedness of trees with their known pedigree relationship. We explored the use of the technique on other related species and highlighted the advantages and limitations of this approach. We found that 92.07% of target sequences in Q. suber and 70.36% of sequences in Fagus sylvatica were captured. We used this SNP resource to estimate genetic relatedness in the mixed oak stand. Mean pairwise genetic relatedness was low within each species with a few values exceeding 0.25 (half-sibs) or 0.5 (full-sibs). Finally, we applied the technique to a long-standing issue in population genetics of trees regarding the relationship between inbreeding and components of fitness. We found very weak signals for inbreeding depression for reproductive success and no signal for growth within both species.