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

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

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.

OakRootRNADB-a consolidated RNA-seq database for coding and noncoding RNA in roots of pedunculate oak (Quercus robur)

The degree to which roots elongate is determined by the expression of genes that regulate root growth in each developmental zone of a root. Most studies have, however, focused on the molecular factors that regulate primary root growth in annual plants. In contrast, the relationship between gene expression and a specific pattern of taproot development and growth in trees is poorly understood. However, the presence of a deeply located taproot, with branching lateral roots, can especially mitigate the effect of insufficient water availability in long-lived trees, such as pedunculated oak. In the present article, we integrated the ribonucleic acid (RNA) sequencing data on roots of oak trees into a single comprehensive database, named OakRootRNADB that contains information on both coding and noncoding RNAs. The sequences in the database also enclose information pertaining to transcription factors, transcriptional regulators and chromatin regulators, as well as a prediction of the cellular localization of a transcript. OakRootRNADB has a user-friendly interface and functional tools that increase access to genomic information. Integrated knowledge of molecular patterns of expression, specifically occurring within and between root zones and within root types, can elucidate the molecular mechanisms regulating taproot growth and enhanced root soil exploration. Database URL https://oakrootrnadb.idpan.poznan.pl/.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

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.

Landscape genomics of Quercus lobata reveals genes involved in local climate adaptation at multiple spatial scales

Understanding how the environment shapes genetic variation provides critical insight about the evolution of local adaptation in natural populations. At multiple spatial scales and multiple geographic contexts within a single species, such information could address a number of fundamental questions about the scale of local adaptation and whether or not the same loci are involved at different spatial scales or geographic contexts. We used landscape genomic approaches from three local elevational transects and rangewide sampling to (a) identify genetic variation underlying local adaptation to environmental gradients in the California endemic oak, Quercus lobata; (b) examine whether putatively adaptive SNPs show signatures of selection at multiple spatial scales; and (c) map putatively adaptive variation to assess the scale and pattern of local adaptation. Of over 10 k single-nucleotide polymorphisms (SNPs) generated with genotyping-by-sequencing, we found signatures of natural selection by climate or local environment at over 600 SNPs (536 loci), some at multiple spatial scales across multiple analyses. Candidate SNPs identified with gene-environment tests (LFMM) at the rangewide scale also showed elevated associations with climate variables compared to the background at both rangewide and elevational transect scales with gradient forest analysis. Some loci overlap with those detected in other oak species, raising the question of whether the same loci might be involved in local climate adaptation in different congeneric species that inhabit different geographic contexts. Mapping landscape patterns of adaptive versus background genetic variation identified regions of marked local adaptation and suggests nonlinear association of candidate SNPs and environmental variables. Taken together, our results offer robust evidence for novel candidate genes for local climate adaptation at multiple spatial scales.

Genome-Wide Identification of Epigenetic Regulators in Quercus suber L

Modifications of DNA and histones, including methylation and acetylation, are critical for the epigenetic regulation of gene expression during plant development, particularly during environmental adaptation processes. However, information on the enzymes catalyzing all these modifications in trees, such as Quercus suber L., is still not available. In this study, eight DNA methyltransferases (DNA Mtases) and three DNA demethylases (DDMEs) were identified in Q. suber. Histone modifiers involved in methylation (35), demethylation (26), acetylation (8), and deacetylation (22) were also identified in Q. suber. In silico analysis showed that some Q. suber DNA Mtases, DDMEs and histone modifiers have the typical domains found in the plant model Arabidopsis, which might suggest a conserved functional role. Additional phylogenetic analyses of the DNA and histone modifier proteins were performed using several plant species homologs, enabling the classification of the Q. suber proteins. A link between the expression levels of each gene in different Q. suber tissues (buds, flowers, acorns, embryos, cork, and roots) with the functions already known for their closest homologs in other species was also established. Therefore, the data generated here will be important for future studies exploring the role of epigenetic regulators in this economically important species.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

From bud formation to flowering: transcriptomic state defines the cherry developmental phases of sweet cherry bud dormancy

Background

Bud dormancy is a crucial stage in perennial trees and allows survival over winter to ensure optimal flowering and fruit production. Recent work highlighted physiological and molecular events occurring during bud dormancy in trees. However, they usually examined bud development or bud dormancy in isolation. In this work, we aimed to further explore the global transcriptional changes happening throughout bud development and dormancy onset, progression and release.

Results

Using next-generation sequencing and modelling, we conducted an in-depth transcriptomic analysis for all stages of flower buds in several sweet cherry (Prunus avium L.) cultivars that are characterized for their contrasted dates of dormancy release. We find that buds in organogenesis, paradormancy, endodormancy and ecodormancy stages are defined by the expression of genes involved in specific pathways, and these are conserved between different sweet cherry cultivars. In particular, we found that DORMANCY ASSOCIATED MADS-box (DAM), floral identity and organogenesis genes are up-regulated during the pre-dormancy stages while endodormancy is characterized by a complex array of signalling pathways, including cold response genes, ABA and oxidation-reduction processes. After dormancy release, genes associated with global cell activity, division and differentiation are activated during ecodormancy and growth resumption. We then went a step beyond the global transcriptomic analysis and we developed a model based on the transcriptional profiles of just seven genes to accurately predict the main bud dormancy stages.

Conclusions

Overall, this study has allowed us to better understand the transcriptional changes occurring throughout the different phases of flower bud development, from bud formation in the summer to flowering in the following spring. Our work sets the stage for the development of fast and cost effective diagnostic tools to molecularly define the dormancy stages. Such integrative approaches will therefore be extremely useful for a better comprehension of complex phenological processes in many species.

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.

Proteomics, Holm Oak (Quercus ilex L.) and Other Recalcitrant and Orphan Forest Tree Species: How do They See Each Other?

Proteomics has had a big impact on plant biology, considered as a valuable tool for several forest species, such as Quercus, Pines, Poplars, and Eucalyptus. This review assesses the potential and limitations of the proteomics approaches and is focused on Quercus ilex as a model species and other forest tree species. Proteomics has been used with Q. ilex since 2003 with the main aim of examining natural variability, developmental processes, and responses to biotic and abiotic stresses as in other species of the genus Quercus or Pinus. As with the progress in techniques in proteomics in other plant species, the research in Q. ilex moved from 2-DE based strategy to the latest gel-free shotgun workflows. Experimental design, protein extraction, mass spectrometric analysis, confidence levels of qualitative and quantitative proteomics data, and their interpretation are a true challenge with relation to forest tree species due to their extreme orphan and recalcitrant (non-orthodox) nature. Implementing a systems biology approach, it is time to validate proteomics data using complementary techniques and integrate it with the -omics and classical approaches. The full potential of the protein field in plant research is quite far from being entirely exploited. However, despite the methodological limitations present in proteomics, there is no doubt that this discipline has contributed to deeper knowledge of plant biology and, currently, is increasingly employed for translational purposes.

Ion Torrent and lllumina, two complementary RNA-seq platforms for constructing the holm oak (Quercus ilex) transcriptome

Transcriptome analysis is widely used in plant biology research to explore gene expression across a large variety of biological contexts such as those related to environmental stress and plant-pathogen interaction. Currently, next generation sequencing platforms are used to obtain a high amount of raw data to build the transcriptome of any plant. Here, we compare Illumina and Ion Torrent sequencing platforms for the construction and analysis of the holm oak (Quercus ilex) transcriptome. Genomic analysis of this forest tree species is a major challenge considering its recalcitrant character and the absence of previous molecular studies. In this study, Quercus ilex raw sequencing reads were obtained from Illumina and Ion Torrent and assembled by three different algorithms, MIRA, RAY and TRINITY. A hybrid transcriptome combining both sequencing technologies was also obtained in this study. The RAY-hybrid assembly generated the most complete transcriptome (1,116 complete sequences of which 1,085 were single copy) with a E90N50 of 1,122 bp. The MIRA-Illumina and TRINITY-Ion Torrent assemblies annotated the highest number of total transcripts (62,628 and 74,058 respectively). MIRA-Ion Torrent showed the highest number of shared sequences (84.8%) with the oak transcriptome. All the assembled transcripts from the hybrid transcriptome were annotated with gene ontology grouping them in terms of biological processes, molecular functions and cellular components. In addition, an in silico proteomic analysis was carried out using the translated assemblies as databases. Those from Ion Torrent showed more proteins compared to the Illumina and hybrid assemblies. This new generated transcriptome represents a valuable tool to conduct differential gene expression studies in response to biotic and abiotic stresses and to assist and validate the ongoing Q. ilex whole genome sequencing.

Gene expression analysis of bud burst process in European hazelnut (Corylus avellana L.) using RNA-Seq

The control of bud burst process depending on temperature is crucial factor in woody perennial plants to survive in unfavorable ecological conditions. Although it has important economic and agronomic values, little information is available on the molecular mechanism of the bud burst process in Corylus avellana. Here for the first time, we conducted a de novo transcriptome-based experiment using eco-dormant leaf bud tissues. Four transcriptome libraries were constructed from the leaf bud tissues and sequenced via Illumina platform. Transcriptome analysis revealed 86,394 unigenes with a mean length of 1189 nt and an N50 of 1916 nt. Among these unigenes, 63,854 (73.78%) of them were annotated by at least one database. De novo assembled transcripts were enriched in phenylpropanoid metabolism, phytohormone biosynthesis and signal transduction pathways. Analyses of phytohormone-associated genes revealed important changes during bud burst, in response to gibberellic acid, auxin, and brassinosteroids. Approximately 2163 putative transcription factors were predicted, of which the largest number of unique transcripts belonged to the MYB transcription factor family. These results contribute to a better understanding of the regulation of bud burst genes in perennial plants.

Diversification, adaptation, and community assembly of the American oaks (Quercus), a model clade for integrating ecology and evolution

Contents Summary 669 I. Model clades for the study and integration of ecology and evolution 670 II. Oaks: an important model clade 671 III. Insights from the history of the American oaks for understanding community assembly and ecosystem dominance 673 IV. Bridging the gap between micro- and macroevolutionary processes relevant to ecology 679 V. How do we reconcile evidence for adaptive evolution with niche conservatism and long-term stasis? 682 VI. High plasticity and within-population genetic variation contribute to population persistence 683 VII. Emerging technologies for tracking functional change 685 VIII. Conclusions 685 Acknowledgements 686 References 686 SUMMARY: Ecologists and evolutionary biologists are concerned with explaining the diversity and composition of the natural world and are aware of the inextricable linkages between ecological and evolutionary processes that maintain the Earth's life support systems. Yet examination of these linkages remains challenging due to the contrasting nature of focal systems and research approaches. Model clades provide a critical means to integrate ecology and evolution, as illustrated by the oaks (genus Quercus), an important model clade, given their ecological dominance, remarkable diversity, and growing phylogenetic, genomic, and ecological data resources. Studies of the clade reveal that their history of sympatric parallel adaptive radiation continues to influence community assembly today, highlighting questions on the nature and extent of coexistence mechanisms. Flexible phenology and hydraulic traits, despite evolutionary stasis, may have enabled adaptation to a wide range of environments within and across species, contributing to their high abundance and diversity. The oaks offer fundamental insights at the intersection of ecology and evolution on the role of diversification in community assembly processes, on the importance of flexibility in key functional traits in adapting to new environments, on factors contributing to persistence of long-lived organisms, and on evolutionary legacies that influence ecosystem function.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

How Quercus ilex L. saplings face combined salt and ozone stress: a transcriptome analysis

BACKGROUND

Similar to other urban trees, holm oaks (Quercus ilex L.) provide a physiological, ecological and social service in the urban environment, since they remove atmospheric pollution. However, the urban environment has several abiotic factors that negatively influence plant life, which are further exacerbated due to climate change, especially in the Mediterranean area. Among these abiotic factors, increased uptake of Na + and Cl - usually occurs in trees in the urban ecosystem; moreover, an excess of the tropospheric ozone concentration in Mediterranean cities further affects plant growth and survival. Here, we produced and annotated a de novo leaf transcriptome of Q. ilex as well as transcripts over- or under-expressed after a single episode of O3 (80 nl l-1, 5 h), a salt treatment (150 mM for 15 days) or a combination of these treatments, mimicking a situation that plants commonly face, especially in urban environments.

RESULTS

Salinity dramatically changed the profile of expressed transcripts, while the short O3 pulse had less effect on the transcript profile. However, the short O3 pulse had a very strong effect in inducing over- or under-expression of some genes in plants coping with soil salinity. Many differentially regulated genes were related to stress sensing and signalling, cell wall remodelling, ROS sensing and scavenging, photosynthesis and to sugar and lipid metabolism. Most differentially expressed transcripts revealed here are in accordance with a previous report on Q. ilex at the physiological and biochemical levels, even though the expression profiles were overall more striking than those found at the biochemical and physiological levels.

CONCLUSIONS

We produced for the first time a reference transcriptome for Q. ilex, and performed gene expression analysis for this species when subjected to salt, ozone and a combination of the two. The comparison of gene expression between the combined salt + ozone treatment and salt or ozone alone showed that even though many differentially expressed genes overlap all treatments, combined stress triggered a unique response in terms of gene expression modification. The obtained results represent a useful tool for studies aiming to investigate the effects of environmental stresses in urban-adapted tree species.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

A molecular approach to drought-induced reduction in leaf CO2 exchange in drought-resistant Quercus ilex

Drought-induced reduction of leaf gas exchange entails a complex regulation of the plant leaf metabolism. We used a combined molecular and physiological approach to understand leaf photosynthetic and respiratory responses of 2-year-old Quercus ilex seedlings to drought. Mild drought stress resulted in glucose accumulation while net photosynthetic CO₂ uptake (P_n ) remained unchanged, suggesting a role of glucose in stress signaling and/or osmoregulation. Simple sugars and sugar alcohols increased throughout moderate-to-very severe drought stress conditions, in parallel to a progressive decline in P_n and the quantum efficiency of photosystem II; by contrast, minor changes occurred in respiration rates until drought stress was very severe. At very severe drought stress, 2-oxoglutarate dehydrogenase complex gene expression significantly decreased, and the abundance of most amino acids dramatically increased, especially that of proline and γ-aminobutyric acid (GABA) suggesting enhanced protection against oxidative damage and a reorganization of the tricarboxylic cycle acid cycle via the GABA shunt. Altogether, our results point to Q. ilex drought tolerance being linked to signaling and osmoregulation by hexoses during early stages of drought stress, and enhanced protection against oxidative damage by polyols and amino acids under severe drought stress.

Novel Insights into Tree Biology and Genome Evolution as Revealed Through Genomics

Reference genome sequences are the key to the discovery of genes and gene families that determine traits of interest. Recent progress in sequencing technologies has enabled a rapid increase in genome sequencing of tree species, allowing the dissection of complex characters of economic importance, such as fruit and wood quality and resistance to biotic and abiotic stresses. Although the number of reference genome sequences for trees lags behind those for other plant species, it is not too early to gain insight into the unique features that distinguish trees from nontree plants. Our review of the published data suggests that, although many gene families are conserved among herbaceous and tree species, some gene families, such as those involved in resistance to biotic and abiotic stresses and in the synthesis and transport of sugars, are often expanded in tree genomes. As the genomes of more tree species are sequenced, comparative genomics will further elucidate the complexity of tree genomes and how this relates to traits unique to trees.

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.

First Draft Assembly and Annotation of the Genome of a California Endemic Oak Quercus lobata Née (Fagaceae)

Oak represents a valuable natural resource across Northern Hemisphere ecosystems, attracting a large research community studying its genetics, ecology, conservation, and management. Here we introduce a draft genome assembly of valley oak (Quercus lobata) using Illumina sequencing of adult leaf tissue of a tree found in an accessible, well-studied, natural southern California population. Our assembly includes a nuclear genome and a complete chloroplast genome, along with annotation of encoded genes. The assembly contains 94,394 scaffolds, totaling 1.17 Gb with 18,512 scaffolds of length 2 kb or longer, with a total length of 1.15 Gb, and a N50 scaffold size of 278,077 kb. The k-mer histograms indicate an diploid genome size of ∼720–730 Mb, which is smaller than the total length due to high heterozygosity, estimated at 1.25%. A comparison with a recently published European oak (Q. robur) nuclear sequence indicates 93% similarity. The Q. lobata chloroplast genome has 99% identity with another North American oak, Q. rubra. Preliminary annotation yielded an estimate of 61,773 predicted protein-coding genes, of which 71% had similarity to known protein domains. We searched 956 Benchmarking Universal Single-Copy Orthologs, and found 863 complete orthologs, of which 450 were present in > 1 copy. We also examined an earlier version (v0.5) where duplicate haplotypes were removed to discover variants. These additional sources indicate that the predicted gene count in Version 1.0 is overestimated by 37–52%. Nonetheless, this first draft valley oak genome assembly represents a high-quality, well-annotated genome that provides a tool for forest restoration and management practices.

Implication of the suberin pathway in adaptation to waterlogging and hypertrophied lenticels formation in pedunculate oak (Quercus robur L.)

Waterlogging causes stressful conditions for perennial species. The temporary overabundance of water in waterlogged soil can induce hypoxia in the rhizosphere, leading to root death, tree decline and even dieback. Two closely related members of the European white oak complex, pedunculate (Quercus robur L.) and sessile (Quercus petraea Matt. Liebl.) oaks, have different ecological characteristics, especially regarding their adaptation to soil waterlogging. The tolerance of waterlogging observed in pedunculate oak is driven principally by its ability to produce adaptive structures, hypertrophied lenticels and adventitious roots, and to switch rapidly its metabolism to the fermentative pathway. This study had two objectives: (i) to identify genes important for adaptation to waterlogging and (ii) to gain insight into the molecular mechanisms involved in hypertrophied lenticel formation in pedunculate oak. We subjected seedlings of the two species to hypoxia by maintaining the water level 2 cm above the collar. The immersed part of the stem (i.e., containing hypertrophied lenticels in pedunculate oak) was sampled after 9 days of waterlogging stress and its gene expression was investigated by RNA-seq. Genes displaying differential expression between the two species were identified with the DESeq R package and a false discovery rate of 0.001. We found that 3705 contigs were differentially regulated between the two species. Twenty-two differentially expressed genes were validated by real-time quantitative polymerase chain reaction. The suberin biosynthesis pathway was found to be upregulated in pedunculate oak, consistent with molecular mechanisms analogous to those operating in the radial oxygen loss barrier in waterlogging-tolerant species.

A comprehensive view of the web-resources related to sericulture

Recent progress in the field of sequencing and analysis has led to a tremendous spike in data and the development of data science tools. One of the outcomes of this scientific progress is development of numerous databases which are gaining popularity in all disciplines of biology including sericulture. As economically important organism, silkworms are studied extensively for their numerous applications in the field of textiles, biomaterials, biomimetics, etc. Similarly, host plants, pests, pathogens, etc. are also being probed to understand the seri-resources more efficiently. These studies have led to the generation of numerous seri-related databases which are extremely helpful for the scientific community. In this article, we have reviewed all the available online resources on silkworm and its related organisms, including databases as well as informative websites. We have studied their basic features and impact on research through citation count analysis, finally discussing the role of emerging sequencing and analysis technologies in the field of seri-data science. As an outcome of this review, a web portal named SeriPort, has been created which will act as an index for the various sericulture-related databases and web resources available in cyberspace.Database URL: http://www.seriport.in/.

High-density linkage mapping and distribution of segregation distortion regions in the oak genome

We developed the densest single-nucleotide polymorphism (SNP)-based linkage genetic map to date for the genus Quercus An 8k gene-based SNP array was used to genotype more than 1,000 full-sibs from two intraspecific and two interspecific full-sib families of Quercus petraea and Quercus robur A high degree of collinearity was observed between the eight parental maps of the two species. A composite map was then established with 4,261 SNP markers spanning 742 cM over the 12 linkage groups (LGs) of the oak genome. Nine genomic regions from six LGs displayed highly significant distortions of segregation. Two main hypotheses concerning the mechanisms underlying segregation distortion are discussed: genetic load vs. reproductive barriers. Our findings suggest a predominance of pre-zygotic to post-zygotic barriers.

Endogenous rhythmic growth in oak trees is regulated by internal clocks rather than resource availability

Common oak trees display endogenous rhythmic growth with alternating shoot and root flushes. To explore the mechanisms involved, microcuttings of the Quercus robur L. clone DF159 were used for (13)C/(15)N labelling in combination with RNA sequencing (RNASeq) transcript profiling of shoots and roots. The effect of plant internal resource availability on the rhythmic growth of the cuttings was tested through inoculation with the ectomycorrhizal fungus Piloderma croceum. Shoot and root flushes were related to parallel shifts in above- and below-ground C and, to a lesser extent, N allocation. Increased plant internal resource availability by P. croceum inoculation with enhanced plant growth affected neither the rhythmic growth nor the associated resource allocation patterns. Two shifts in transcript abundance were identified during root and shoot growth cessation, and most concerned genes were down-regulated. Inoculation with P. croceum suppressed these transcript shifts in roots, but not in shoots. To identify core processes governing the rhythmic growth, functions [Gene Ontology (GO) terms] of the genes differentially expressed during the growth cessation in both leaves and roots of non-inoculated plants and leaves of P. croceum-inoculated plants were examined. Besides genes related to resource acquisition and cell development, which might reflect rather than trigger rhythmic growth, genes involved in signalling and/or regulated by the circadian clock were identified. The results indicate that rhythmic growth involves dramatic oscillations in plant metabolism and gene regulation between below- and above-ground parts. Ectomycorrhizal symbiosis may play a previously unsuspected role in smoothing these oscillations without modifying the rhythmic growth pattern.

Large scale transcriptome analysis reveals interplay between development of forest trees and a beneficial mycorrhiza helper bacterium

BACKGROUND

Pedunculate oak, Quercus robur is an abundant forest tree species that hosts a large and diverse community of beneficial ectomycorrhizal fungi (EMFs), whereby ectomycorrhiza (EM) formation is stimulated by mycorrhiza helper bacteria such as Streptomyces sp. AcH 505. Oaks typically grow rhythmically, with alternating root flushes (RFs) and shoot flushes (SFs). We explored the poorly understood mechanisms by which oaks integrate signals induced by their beneficial microbes and endogenous rhythmic growth at the level of gene expression. To this end, we compared transcript profiles of oak microcuttings at RF and SF during interactions with AcH 505 alone and in combination with the basidiomycetous EMF Piloderma croceum.

RESULTS

The local root and distal leaf responses to the microorganisms differed substantially. More genes involved in the recognition of bacteria and fungi, defence and cell wall remodelling related transcription factors (TFs) were differentially expressed in the roots than in the leaves of oaks. In addition, interaction with AcH 505 and P. croceum affected the expression of a higher number of genes during SF than during RF, including AcH 505 elicited defence response, which was attenuated by co-inoculation with P. croceum in the roots during SF. Genes encoding leucine-rich receptor-like kinases (LRR-RLKs) and proteins (LRR-RLPs), LRR containing defence response regulators, TFs from bZIP, ERF and WRKY families, xyloglucan cell wall transglycolases/hydrolases and exordium proteins were differentially expressed in both roots and leaves of plants treated with AcH 505. Only few genes, including specific RLKs and TFs, were induced in both AcH 505 and co-inoculation treatments.

CONCLUSION

Treatment with AcH 505 induces and maintains the expression levels of signalling genes encoding candidate receptor protein kinases and TFs and leads to differential expression of cell wall modification related genes in pedunculate oak microcuttings. Local gene expression response to AcH 505 alone and in combination with P. croceum are more pronounced when roots are in resting stages, possibly due to the fact that non growing roots re-direct their activity towards plant defence rather than growth.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole-genome shotgun (WGS) approach, without the use of costly and time-consuming methods, such as fosmid or BAC clone-based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS-FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired-end and mate-pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired-end reads and contaminants detected, resulting in a total of 17,910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.