The draft genome sequence of cork oak

Rhytidome- and cork-type barks of holm oak, cork oak and their hybrids highlight processes leading to cork formation

BACKGROUND

The periderm is basic for land plants due to its protective role during radial growth, which is achieved by the polymers deposited in the cell walls. In most trees, like holm oak, the first periderm is frequently replaced by subsequent internal periderms yielding a heterogeneous outer bark made of a mixture of periderms and phloem tissues, known as rhytidome. Exceptionally, cork oak forms a persistent or long-lived periderm which results in a homogeneous outer bark of thick phellem cell layers known as cork. Cork oak and holm oak distribution ranges overlap to a great extent, and they often share stands, where they can hybridize and produce offspring showing a rhytidome-type bark.

RESULTS

Here we use the outer bark of cork oak, holm oak, and their natural hybrids to analyse the chemical composition, the anatomy and the transcriptome, and further understand the mechanisms underlying periderm development. We also include a unique natural hybrid individual corresponding to a backcross with cork oak that, interestingly, shows a cork-type bark. The inclusion of hybrid samples showing rhytidome-type and cork-type barks is valuable to approach cork and rhytidome development, allowing an accurate identification of candidate genes and processes. The present study underscores that abiotic stress and cell death are enhanced in rhytidome-type barks whereas lipid metabolism and cell cycle are enriched in cork-type barks. Development-related DEGs showing the highest expression, highlight cell division, cell expansion, and cell differentiation as key processes leading to cork or rhytidome-type barks.

CONCLUSION

Transcriptome results, in agreement with anatomical and chemical analyses, show that rhytidome and cork-type barks are active in periderm development, and suberin and lignin deposition. Development and cell wall-related DEGs suggest that cell division and expansion are upregulated in cork-type barks whereas cell differentiation is enhanced in rhytidome-type barks.

Genomic basis of the distinct biosynthesis of β-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that β-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for β-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different β-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding β-glucogallin and HT biosynthesis in closely related oak species.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Wounding-Related Signaling Is Integrated within the Auxin-Response Framework to Induce Adventitious Rooting in Chestnut

Wounding and exogenous auxin are needed to induce adventitious roots in chestnut microshoots. However, the specific inductive role of wounding has not been characterized in this species. In the present work, two main goals were established: First, we prompted to optimize exogenous auxin treatments to improve the overall health status of the shoots at the end of the rooting cycle. Second, we developed a time-series transcriptomic analysis to compare gene expression in response to wounding alone and wounding plus auxin, focusing on the early events within the first days after treatments. Results suggest that the expression of many genes involved in the rooting process is under direct or indirect control of both stimuli. However, specific levels of expression of relevant genes are only attained when both treatments are applied simultaneously, leading to the successful development of roots. In this sense, we have identified four transcription factors upregulated by auxin (CsLBD16, CsERF113, Cs22D and CsIAA6), with some of them also being induced by wounding. The highest expression levels of these genes occurred when wounding and auxin treatments were applied simultaneously, correlating with the rooting response of the shoots. The results of this work clarify the genetic nature of the wounding response in chestnut, its relation to adventitious rooting, and might be helpful in the development of more specific protocols for the vegetative propagation of this species.

The genome assembly of Island Oak (Quercus tomentella), a relictual island tree species

Island oak (Quercus tomentella) is a rare relictual island tree species that exists only on six islands off the coast of California and Mexico, but was once widespread throughout mainland California. Currently, this species is endangered by threats such as non-native plants, grazing animals, and human removal. Efforts for conservation and restoration of island oak currently underway could benefit from information about its range-wide genetic structure and evolutionary history. Here we present a high-quality genome assembly for Q. tomentella, assembled using PacBio HiFi and Omni-C sequencing, developed as part of the California Conservation Genomics Project (CCGP). The resulting assembly has a length of 781 Mb, with a contig N50 of 22.0 Mb and a scaffold N50 of 63.4 Mb. This genome assembly will provide a resource for genomics-informed conservation of this rare oak species. Additionally, this reference genome will be the first one available for a species in Quercus section Protobalanus, a unique oak clade present only in western North America.

ContScout: sensitive detection and removal of contamination from annotated genomes

Contamination of genomes is an increasingly recognized problem affecting several downstream applications, from comparative evolutionary genomics to metagenomics. Here we introduce ContScout, a precise tool for eliminating foreign sequences from annotated genomes. It achieves high specificity and sensitivity on synthetic benchmark data even when the contaminant is a closely related species, outperforms competing tools, and can distinguish horizontal gene transfer from contamination. A screen of 844 eukaryotic genomes for contamination identified bacteria as the most common source, followed by fungi and plants. Furthermore, we show that contaminants in ancestral genome reconstructions lead to erroneous early origins of genes and inflate gene loss rates, leading to a false notion of complex ancestral genomes. Taken together, we offer here a tool for sensitive removal of foreign proteins, identify and remove contaminants from diverse eukaryotic genomes and evaluate their impact on phylogenomic analyses.

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness

Quercus variabilis is a deciduous woody species with high ecological and economic value, and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

Identification of Functional Brassinosteroid Receptor Genes in Oaks and Functional Analysis of QmBRI1

Brassinosteroids (BRs) play important regulatory roles in plant growth and development, with functional BR receptors being crucial for BR recognition or signaling. Although functional BR receptors have been extensively studied in herbaceous plants, they remain largely under-studied in forest tree species. In this study, nine BR receptors were identified in three representative oak species, of which BRI1s and BRL1s were functional BR receptors. Dispersed duplications were a driving force for oak BR receptor expansion, among which the Brassinosteroid-Insensitive-1 (BRI1)-type genes diverged evolutionarily from most rosids. In oak BRI1s, we identified that methionine in the conserved Asn-Gly-Ser-Met (NGSM) motif was replaced by isoleucine and that the amino acid mutation occurred after the divergence of Quercus and Fagus. Compared with QmBRL1, QmBRI1 was relatively highly expressed during BR-induced xylem differentiation and in young leaves, shoots, and the phloem and xylem of young stems of Quercus mongolica. Based on Arabidopsis complementation experiments, we proved the important role of QmBRI1 in oak growth and development, especially in vascular patterning and xylem differentiation. These findings serve as an important supplement to the findings of the structural, functional and evolutionary studies on functional BR receptors in woody plants and provide the first example of natural mutation occurring in the conserved BR-binding region (NGSM motif) of angiosperm BRI1s.

Telomere-to-telomere and haplotype-resolved genome assembly of the Chinese cork oak (Quercus variabilis)

The Quercus variabilis, a deciduous broadleaved tree species, holds significant ecological and economical value. While a chromosome-level genome for this species has been made available, it remains riddled with unanchored sequences and gaps. In this study, we present a nearly complete comprehensive telomere-to-telomere (T2T) and haplotype-resolved reference genome for Q. variabilis. This was achieved through the integration of ONT ultra-long reads, PacBio HiFi long reads, and Hi-C data. The resultant two haplotype genomes measure 789 Mb and 768 Mb in length, with a contig N50 of 65 Mb and 56 Mb, and were anchored to 12 allelic chromosomes. Within this T2T haplotype-resolved assembly, we predicted 36,830 and 36,370 protein-coding genes, with 95.9% and 96.0% functional annotation for each haplotype genome. The availability of the T2T and haplotype-resolved reference genome lays a solid foundation, not only for illustrating genome structure and functional genomics studies but also to inform and facilitate genetic breeding and improvement of cultivated Quercus species.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa"

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

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

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

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

A chromosome-scale genome assembly of Castanopsis hystrix provides new insights into the evolution and adaptation of Fagaceae species

Fagaceae species dominate forests and shrublands throughout the Northern Hemisphere, and have been used as models to investigate the processes and mechanisms of adaptation and speciation. Compared with the well-studied genus Quercus, genomic data is limited for the tropical-subtropical genus Castanopsis. Castanopsis hystrix is an ecologically and economically valuable species with a wide distribution in the evergreen broad-leaved forests of tropical-subtropical Asia. Here, we present a high-quality chromosome-scale reference genome of C. hystrix, obtained using a combination of Illumina and PacBio HiFi reads with Hi-C technology. The assembled genome size is 882.6 Mb with a contig N50 of 40.9 Mb and a BUSCO estimate of 99.5%, which are higher than those of recently published Fagaceae species. Genome annotation identified 37,750 protein-coding genes, of which 97.91% were functionally annotated. Repeat sequences constituted 50.95% of the genome and LTRs were the most abundant repetitive elements. Comparative genomic analysis revealed high genome synteny between C. hystrix and other Fagaceae species, despite the long divergence time between them. Considerable gene family expansion and contraction were detected in Castanopsis species. These expanded genes were involved in multiple important biological processes and molecular functions, which may have contributed to the adaptation of the genus to a tropical-subtropical climate. In summary, the genome assembly of C. hystrix provides important genomic resources for Fagaceae genomic research communities, and improves understanding of the adaptation and evolution of forest trees.

Laser Microdissection of Woody and Suberized Plant Tissues for RNA-Seq Analysis

An accurate profile of gene expression at a cellular level can contribute to a better understanding of biological processes and complexities involved in regulatory mechanism of woody plants. Laser microdissection is one technique that allows isolation of specific, target cells or tissue from a heterogeneous cell population. This technique entails microscopic visualization of the selected tissue and use a laser beam to separate the desired cells from surrounding tissue. Initial identification of these cells is made based on morphology and/or histological staining. Some works have been made in several tissues and plant models. However, there are few studies of laser microdissection application in woody species, particularly, lignified and suberized cells. Moreover, the presence of high level of suberin in cell walls can be a big challenge for the application of this approach. In our study it was developed a technique for tissue isolation, using laser microdissection of four different plant cell types (phellogen, lenticels, cortex and xylem) from woody tissues of cork oak (Quercus suber), followed by RNA extraction and RNA-Seq. We tested several methodologies regarding laser microdissection, cryostat equipments, fixation treatments, duration of single-cells collection and number of isolated cells by laser microdissection and RNA extraction procedures. A simple and efficient protocol for tissue isolation by laser microdissection and RNA purification was obtained, with a final method validation of RNA-Seq analysis. The optimized methodology combining RNA-Seq for expression analysis will contribute to elucidate the molecular pathways associated with different development processes of the xylem and phellem in oaks, including the lenticular channels formation.

Germplasm Resources of Oaks (Quercus L.) in China: Utilization and Prospects

Oaks exhibit unique biological characteristics and high adaptability to complex climatic and soil conditions. They are widely distributed across various regions, spanning 40 degrees latitude and 75 degrees longitude. The total area of oak forest in China is 16.72 million hm². There are 60 lineages of Quercus in China, including 49 species, seven varieties, and four subgenera. Archaeological data indicate that oaks were already widely distributed in ancient times, and they are dominant trees in vast regions of China's forests. In addition, the acorn was an important food for ancestral humans, and it has accompanied human civilization since the early Paleolithic. Diverse oak species are widely distributed and have great functional value, such as for greening, carbon sequestration, industrial and medicinal uses, and insect rearing. Long-term deforestation, fire, diseases, and pests have led to a continuous decline in oak resources. This study discusses the Quercus species and their distribution in China, ecological adaptation, and the threats facing the propagation and growth of oaks in a changing world. This will give us a better understanding of Quercus resources, and provide guidance on how to protect and better utilize germplasm resources in China. The breeding of new varieties, pest control, and chemical and molecular research also need to be strengthened in future studies.

GlPS1 overexpression accumulates coumarin secondary metabolites in transgenic Arabidopsis

The dried root of Glehnia littoralis is a traditional Chinese herbal medicine mainly used to treat lung diseases and plays an important role in fighting coronavirus disease 2019 pneumonia in China. This study focused on the key enzyme gene GlPS1 for furanocoumarin synthesis in G. littoralis. In the 35S:GlPS1 transgenic Arabidopsis study, the Arabidopsis thaliana-overexpressing GlPS1 gene was more salt-tolerant than Arabidopsis in the blank group. Metabolomics analysis showed 30 differential metabolites in Arabidopsis, which overexpressed the GlPS1 gene. Twelve coumarin compounds were significantly upregulated, and six of these coumarin compounds were not detected in the blank group. Among these differential coumarin metabolites, isopimpinellin and aesculetin have been annotated by the Kyoto Encyclopedia of Genes and Genomes and isopimpinellin was not detected in the blank group. Through structural comparison, imperatorin was formed by dehydration and condensation of zanthotoxol and a molecule of isoprenol, and the difference between them was only one isoprene. Results showed that the GlPS1 gene positively regulated the synthesis of coumarin metabolites in A. thaliana and at the same time improved the salt tolerance of A. thaliana.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11240-022-02427-w.

Nitrogen Acquisition and Transport in the Ectomycorrhizal Symbiosis-Insights from the Interaction between an Oak Tree and Pisolithus tinctorius

In temperate forests, the roots of various tree species are colonized by ectomycorrhizal fungi, which have a key role in the nitrogen nutrition of their hosts. However, not much is known about the molecular mechanisms related to nitrogen metabolism in ectomycorrhizal plants. This study aimed to evaluate the nitrogen metabolic response of oak plants when inoculated with the ectomycorrhizal fungus Pisolithus tinctorius. The expression of candidate genes encoding proteins involved in nitrogen uptake and assimilation was investigated in ectomycorrhizal roots. We found that three oak ammonium transporters were over-expressed in root tissues after inoculation, while the expression of amino acid transporters was not modified, suggesting that inorganic nitrogen is the main form of nitrogen transferred by the symbiotic fungus into the roots of the host plant. Analysis by heterologous complementation of a yeast mutant defective in ammonium uptake and GFP subcellular protein localization clearly confirmed that two of these genes encode functional ammonium transporters. Structural similarities between the proteins encoded by these ectomycorrhizal upregulated ammonium transporters, and a well-characterized ammonium transporter from E. coli, suggest a similar transport mechanism, involving deprotonation of NH4+, followed by diffusion of uncharged NH3 into the cytosol. This view is supported by the lack of induction of NH4+ detoxifying mechanisms, such as the GS/GOGAT pathway, in the oak mycorrhizal roots.

Transcriptomics Analysis Reveals a Putative Role for Hormone Signaling and MADS-Box Genes in Mature Chestnut Shoots Rooting Recalcitrance

Maturation imposes several changes in plants, which are particularly drastic in the case of trees. In recalcitrant woody species, such as chestnut (Castanea sativa Mill.), one of the major maturation-related shifts is the loss of the ability to form adventitious roots in response to auxin treatment as the plant ages. To analyze the molecular mechanisms underlying this phenomenon, an in vitro model system of two different lines of microshoots derived from the same field-grown tree was established. While juvenile-like shoots root readily when treated with exogenous auxin, microshoots established from the crown of the tree rarely form roots. In the present study, a transcriptomic analysis was developed to compare the gene expression patterns in both types of shoots 24 h after hormone and wounding treatment, matching the induction phase of the process. Our results support the hypothesis that the inability of adult chestnut tissues to respond to the inductive treatment relies in a deep change of gene expression imposed by maturation that results in a significant transcriptome modification. Differences in phytohormone signaling seem to be the main cause for the recalcitrant behavior of mature shoots, with abscisic acid and ethylene negatively influencing the rooting ability of the chestnut plants. We have identified a set of related MADS-box genes whose expression is modified but not suppressed by the inductive treatment in mature shoots, suggesting a putative link of their activity with the rooting-recalcitrant behavior of this material. Overall, distinct maturation-derived auxin sensibility and homeostasis, and the related modifications in the balance with other phytohormones, seem to govern the outcome of the process in each type of shoots.

Strategies to mitigate shifts in red oak (Quercus sect. Lobatae) distribution under a changing climate

Red oaks (Quercus sect. Lobatae) are a taxonomic group of hardwood trees, which occur in swamp forests, subtropical chaparral and savannahs from Columbia to Canada. They cover a wide range of ecological niches, and many species are thought to be able to cope with current trends in climate change. Genus Quercus encompasses ca. 500 species, of which ca. 80 make up sect. Lobatae. Species diversity is greatest within the southeastern USA and within the northern and eastern regions of Mexico. This review discusses the weak reproductive barriers between species of red oaks and the effects this has on speciation and niche range. Distribution and diversity have been shaped by drought adaptations common to the species of sect. Lobatae, which enable them to fill various xeric niches across the continent. Drought adaptive traits of this taxonomic group include deciduousness, deep tap roots, ring-porous xylem, regenerative stump sprouting, greater leaf thickness and smaller stomata. The complex interplay between these anatomical and morphological traits has given red oaks features of drought tolerance and avoidance. Here, we discuss physiological and genetic components of these adaptations to address how many species of sect. Lobatae reside within xeric sites and/or sustain normal metabolic function during drought. Although extensive drought adaptation appears to give sect. Lobatae a resilience to climate change, aging tree stands, oak life history traits and the current genetic structures place many red oak species at risk. Furthermore, oak decline, a complex interaction between abiotic and biotic agents, has severe effects on red oaks and is likely to accelerate species decline and fragmentation. We suggest that assisted migration can be used to avoid species fragmentation and increase climate change resilience of sect. Lobatae.

A high-quality genome assembly and annotation of Quercus acutissima Carruth

INTRODUCTION

Quercus acutissima is an economic and ecological tree species often used for afforestation of arid and semi-arid lands and is considered as an excellent tree for soil and water conservation.

METHODS

Here, we combined PacBio long reads, Hi-C, and Illumina short reads to assemble Q. acutissima genome.

RESULTS

We generated a 957.1 Mb genome with a contig N50 of 1.2 Mb and scaffold N50 of 77.0 Mb. The repetitive sequences constituted 55.63% of the genome, among which long terminal repeats were the majority and accounted for 23.07% of the genome. Ab initio, homology-based and RNA sequence-based gene prediction identified 29,889 protein-coding genes, of which 82.6% could be functionally annotated. Phylogenetic analysis showed that Q. acutissima and Q. variabilis were differentiated around 3.6 million years ago, and showed no evidence of species-specific whole genome duplication.

CONCLUSION

The assembled and annotated high-quality Q. acutissima genome not only promises to accelerate the species molecular biology studies and breeding, but also promotes genome level evolutionary studies.

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/.

Transcriptome Analysis of Persian Oak (Quercus brantii L.) Decline Using RNA-seq Technology

Since the late 1980s, the oak decline has affected the Zagros oak forests in western Iran. Persian oak (Quercus brantii L.) the most important tree species of these forests has been damaged more than any other plant species. In the present study, the RNA sequencing technique was used for the first time to identify key genes and molecular mechanisms involved in Persian oak decline. The RNA was extracted from the leaves of healthy and declined oak trees, and sequenced using the Illumina HiSeq 2500 platform (2 × 150 bp paired-end reads). De novo transcriptome assembly of Persian oak revealed 56,743 unigenes and 6049 differentially expressed genes (DEGs) between declined and control samples. The results of gene ontology analysis showed that most of the DEGs involved in oak decline belong to the group of stress-responsive genes. In general, oak decline samples showed significant reductions in gene expression associated with "photosynthesis and storage of sugar" and "protein synthesis and related processes." Additionally, DEGs related to the starch degradation pathway were up-regulated, whereas DEGs associated with acetate-mevalonate (MVA), biosynthesis of lignin, and lignases pathways were down-regulated. The present study's findings can be an effective step in identifying the genes involved in oak decline and deciphering the relationship between this phenomenon and biotic and abiotic stresses.

A chromosome-scale genome assembly of Quercus gilva: Insights into the evolution of Quercus section Cyclobalanopsis (Fagaceae)

Quercus gilva is an ecologically and economically important species of Quercus section Cyclobalanopsis and is a dominant species in evergreen broad-leaved forests in subtropical regions of East Asia. In the present study, we reported a high-quality chromosome-scale genome assembly of Q. gilva, the first reference genome for section Cyclobalanopsis, using the combination of Illumina and PacBio sequencing with Hi-C technologies. The assembled genome size of Q. gilva was 889.71 Mb, with a contig number of 773 and a contig N50 of 28.32 Mb. Hi-C scaffolding anchored 859.07 Mb contigs (96.54% of the assembled genome) onto 12 pseudochromosomes, with a scaffold N50 of 70.35 Mb. A combination of de novo, homology-based, and transcript-based predictions predicted a final set of 36,442 protein-coding genes distributed on 12 pseudochromosomes, and 97.73% of them were functionally annotated. A total of 535.64 Mb (60.20%) of repetitive sequences were identified. Genome evolution analysis revealed that Q. gilva was most closely related to Q. suber and they diverged at 40.35 Ma, and Q. gilva did not experience species-specific whole-genome duplication in addition to the ancient gamma (γ) whole-genome triplication event shared by core eudicot plants. Q. gilva underwent considerable gene family expansion and contraction, with 598 expanded and 6,509 contracted gene families detected. The first chromosome-scale genome of Q. gilva will promote its germplasm conservation and genetic improvement and provide essential resources for better studying the evolution of Quercus section Cyclobalanopsis.

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.

Symbiotic Fungi of an Ambrosia Beetle Alter the Volatile Bouquet of Cork Oak Seedlings

In Portugal, fungal symbionts of the ambrosia beetle Platypus cylindrus affect tree vigor of cork oak (Quercus suber) and are linked with the cork oak decline process. Fungal symbionts play crucial roles in the life history of bark and ambrosia beetles and recent work indicates complex interactions on the fungal and plant metabolic level. Colonized trees may respond with an array of currently unknown volatile metabolites being indicative of such interactions, acting as infochemicals with their environment. In this study, we examined volatile organic compounds (VOCs) of cork oak seedlings wound inoculated with strains of three fungal associates of P. cylindrus (Raffaelea montetyi, R. quercina, and Ceratocystiopsis sp. nov.) over a 45-day period by means of thermodesorption gas chromatography-mass spectrometry techniques. Fungal strains induced largely quantitative but species-specific changes among the 58 VOCs characterized. Overall, monoterpenes-the major volatiles of cork oak foliage-were significantly reduced, possibly a result of fungal biotransformation. Acetophenone, sulcatone, and nonanal-volatiles known for mediating ambrosia beetle behavior-increased in response to fungal inoculation. Qualitative VOC profiles of excised tissue of wood lesions (21 VOCs) and pure fungal cultures (60 VOCs) showed little overlap with seedling VOCs, indicating their plant-derived but fungal-induced origin. This chemoecological study expands on the limited knowledge of VOCs as infochemicals emitted from oak trees threatened by oak decline in relation to beetle-vectored ophiostomatoid fungi. It opens new avenues of research to clarify mutualistic or pathogenic aspects of these complex symbiotic interactions and develop new control strategies for P. cylindrus, including its mycobiota.

Genome-wide identification, phylogeny, and gene duplication of the epigenetic regulators in Fagaceae

Epigenetic regulators are proteins involved in controlling gene expression. Information about the epigenetic regulators within the Fagaceae, a relevant family of trees and shrubs of the northern hemisphere ecosystems, is scarce. With the intent to characterize these proteins in Fagaceae, we searched for orthologs of DNA methyltransferases (DNMTs) and demethylases (DDMEs) and Histone modifiers involved in acetylation (HATs), deacetylation (HDACs), methylation (HMTs), and demethylation (HDMTs) in Fagus, Quercus, and Castanea genera. Blast searches were performed in the available genomes, and freely available RNA-seq data were used to de novo assemble transcriptomes. We identified homologs of seven DNMTs, three DDMEs, six HATs, 11 HDACs, 32 HMTs, and 21 HDMTs proteins. Protein analysis showed that most of them have the putative characteristic domains found in these protein families, which suggests their conserved function. Additionally, to elucidate the evolutionary history of these genes within Fagaceae, paralogs were identified, and phylogenetic analyses were performed with DNA and histone modifiers. We detected duplication events in all species analyzed with higher frequency in Quercus and Castanea and discuss the evidence of transposable elements adjacent to paralogs and their involvement in gene duplication. The knowledge gathered from this work is a steppingstone to upcoming studies concerning epigenetic regulation in this economically important family of Fagaceae.

Spatiotemporal development of suberized barriers in cork oak taproots

The longevity and high activity of the cork cambium (or phellogen) from Quercus suber L. (cork oak) are the cornerstones for the sustainable exploitation of a unique raw material. Cork oak is a symbolic model to study cork development and cell wall suberization, yet most genetic and molecular studies on these topics have targeted other model plants. In this study, we explored the potential of taproots as a model system to study phellem development and suberization in cork oak, thereby avoiding the time constraints imposed when studying whole plants. In roots, suberin deposition is found in mature endodermis cells during primary development and in phellem cells during secondary development. By investigating the spatiotemporal characteristics of both endodermis and phellem suberization in young seedling taproots, we demonstrated that secondary growth and phellogen activity are initiated very early in cork oak taproots (approx. 8 days after sowing). We further compared the transcriptomic profile of root segments undergoing primary (PD) and secondary development (SD) and identified multiple candidate genes with predicted roles in cell wall modifications, mainly lignification and suberization, in addition to several regulatory genes, particularly transcription factor- and hormone-related genes. Our results indicate that the molecular regulation of suberization and secondary development in cork oak roots is relatively conserved with other species. The provided morphological characterization creates new opportunities to allow a faster assessment of phellogen activity (as compared with studies using stem tissues) and to tackle fundamental questions regarding its regulation.

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.

High-quality genome and methylomes illustrate features underlying evolutionary success of oaks

The genus Quercus, which emerged ∼55 million years ago during globally warm temperatures, diversified into ∼450 extant species. We present a high-quality de novo genome assembly of a California endemic oak, Quercus lobata, revealing features consistent with oak evolutionary success. Effective population size remained large throughout history despite declining since early Miocene. Analysis of 39,373 mapped protein-coding genes outlined copious duplications consistent with genetic and phenotypic diversity, both by retention of genes created during the ancient γ whole genome hexaploid duplication event and by tandem duplication within families, including numerous resistance genes and a very large block of duplicated DUF247 genes, which have been found to be associated with self-incompatibility in grasses. An additional surprising finding is that subcontext-specific patterns of DNA methylation associated with transposable elements reveal broadly-distributed heterochromatin in intergenic regions, similar to grasses. Collectively, these features promote genetic and phenotypic variation that would facilitate adaptability to changing environments.

A chromosome-scale genome assembly of the Mongolian oak (Quercus mongolica)

Mongolian oak (Quercus mongolica Fisch.) is an ecologically and economically important white oak species native to and widespread in the temperate zone of East Asia. Here, we present a chromosome-scale reference genome assembly of Q. mongolica, a representative white oak species, by combining Illumina and PacBio data with Hi-C mapping technologies that is the first reference genome created for an Asian oak. Our results showed that the PacBio draft genome size was 809.84 Mb, with a BUSCO complete gene percentage of 92.71%. Hi-C scaffolding anchored 774.59 Mb contigs (95.65% of draft assembly) onto 12 pseudochromosomes. The contig N50 and scaffold N50 were 2.64 Mb and 66.74 Mb, respectively. Of the 36,553 protein-coding genes predicted in the study, approximately 95% had functional annotations in public databases. A total of 435.34 Mb (53.75% of the genome) of repetitive sequences were predicted in the assembled genome. Genome evolution analysis showed that Q. mongolica is closely related to Q. robur from Europe, and they shared a common ancestor ~11.8 million years ago. Gene family evolution analysis of Q. mongolica revealed that the nucleotide-binding site (NBS)-encoding gene family related to disease resistance was significantly contracted, whereas the ECERIFERUM 1 (CER1) homologous genes related to cuticular wax biosynthesis was significantly expanded. This pioneering Asian oak genome resource represents an important supplement to the oak genomics community and will improve our understanding of Asian white oak biology and evolution.

Kinship Analysis and Pedigree Reconstruction of a Natural Regenerated Cork Oak (Quercus suber) Population

Cork oak (Quercus suber L.) is a valuable forest species in the western Mediterranean Basin due to its ecological value and the production of cork (a renewable natural material). Cork quality depends on the genetic background and cork oak environment, which has long been recognized. As no cork oak genetic trials with pedigree information were available, the inference of the genetic relatedness between individuals from molecular markers can potentially be applied to natural populations. This work aimed to investigate the potential of performing kinship prediction and pedigree reconstruction by SNP genotyping a natural cork oak population. A total of 494 trees located in Portugal were genotyped with 8K SNPs. The raw SNP set was filtered differently, producing four SNP sets that were further filtered by missing data, genotype frequency, and minor allele frequency. For each set, an identity by descent (IBD) matrix was generated to perform the relationship prediction, revealing from 22,114 to 23,859 relationships. Familial categories from the first to the third degree were able to be assigned. The feasibility of SNP genotyping for future studies on the kinship analysis and pedigree reconstruction of cork oak populations was demonstrated. The information produced may be used in further breeding and conservation programs for cork oak.

Identification of Proteases and Protease Inhibitors in Seeds of the Recalcitrant Forest Tree Species Quercus ilex

Proteases and protease inhibitors have been identified in the recalcitrant species Quercus ilex using in silico and wet methods, with focus on those present in seeds during germination. In silico analyses showed that the Q. ilex transcriptome database contained 2,240 and 97 transcripts annotated as proteases and protease inhibitors, respectively. They belonged to the different families according to MEROPS, being the serine and metallo ones the most represented. The data were compared with those previously reported for other Quercus species, including Q. suber, Q. lobata, and Q. robur. Changes in proteases and protease inhibitors alongside seed germination in cotyledon and embryo axis tissues were assessed using proteomics and in vitro and in gel activity assays. Shotgun (LC-MSMS) analysis of embryo axes and cotyledons in nonviable (NV), mature (T1) and germinated (T3) seeds allowed the identification of 177 proteases and 12 protease inhibitors, mostly represented by serine and metallo types. Total protease activity, as determined by in vitro assays using azocasein as substrate, was higher in cotyledons than in embryo axes. There were not differences in activity among cotyledon samples, while embryo axis peaked at germinated T4 stage. Gel assays revealed the presence of protease activities in at least 10 resolved bands, in the Mr range of 60-260 kDa, being some of them common to cotyledons and embryo axes in either nonviable, mature, and germinated seeds. Bands showing quantitative or qualitative changes upon germination were observed in embryo axes but not in cotyledons at Mr values of 60-140 kDa. Proteomics shotgun analysis of the 10 bands with protease activity supported the results obtained in the overall proteome analysis, with 227 proteases and 3 protease inhibitors identified mostly represented by the serine, cysteine, and metallo families. The combined use of shotgun proteomics and protease activity measurements allowed the identification of tissue-specific (e.g., cysteine protease inhibitors in embryo axes of mature acorns) and stage-specific proteins (e.g., those associated with mobilization of storage proteins accumulated in T3 stage). Those proteins showing differences between nonviable and viable seeds could be related to viability, and those variables between mature and germinated could be associated with the germination process. These differences are observed mostly in embryo axes but not in cotyledons. Among them, those implicated in mobilization of reserve proteins, such as the cathepsin H cysteine protease and Clp proteases, and also the large number of subunits of the CNS and 26S proteasome complex differentially identified in embryos of the several stages suggests that protein degradation via CNS/26S plays a major role early in germination. Conversely, aspartic proteases such as nepenthesins were exclusively identified in NV seeds, so their presence could be used as indicator of nonviability.

A chromosome-level genome assembly of the Chinese cork oak (Quercus variabilis)

Quercus variabilis (Fagaceae) is an ecologically and economically important deciduous broadleaved tree species native to and widespread in East Asia. It is a valuable woody species and an indicator of local forest health, and occupies a dominant position in forest ecosystems in East Asia. However, genomic resources from Q. variabilis are still lacking. Here, we present a high-quality Q. variabilis genome generated by PacBio HiFi and Hi-C sequencing. The assembled genome size is 787 Mb, with a contig N50 of 26.04 Mb and scaffold N50 of 64.86 Mb, comprising 12 pseudo-chromosomes. The repetitive sequences constitute 67.6% of the genome, of which the majority are long terminal repeats, accounting for 46.62% of the genome. We used ab initio, RNA sequence-based and homology-based predictions to identify protein-coding genes. A total of 32,466 protein-coding genes were identified, of which 95.11% could be functionally annotated. Evolutionary analysis showed that Q. variabilis was more closely related to Q. suber than to Q. lobata or Q. robur. We found no evidence for species-specific whole genome duplications in Quercus after the species had diverged. This study provides the first genome assembly and the first gene annotation data for Q. variabilis. These resources will inform the design of further breeding strategies, and will be valuable in the study of genome editing and comparative genomics in oak species.

PANNZER-A practical tool for protein function prediction

The facility of next-generation sequencing has led to an explosion of gene catalogs for novel genomes, transcriptomes and metagenomes, which are functionally uncharacterized. Computational inference has emerged as a necessary substitute for first-hand experimental evidence. PANNZER (Protein ANNotation with Z-scoRE) is a high-throughput functional annotation web server that stands out among similar publically accessible web servers in supporting submission of up to 100,000 protein sequences at once and providing both Gene Ontology (GO) annotations and free text description predictions. Here, we demonstrate the use of PANNZER and discuss future plans and challenges. We present two case studies to illustrate problems related to data quality and method evaluation. Some commonly used evaluation metrics and evaluation datasets promote methods that favor unspecific and broad functional classes over more informative and specific classes. We argue that this can bias the development of automated function prediction methods. The PANNZER web server and source code are available at http://ekhidna2.biocenter.helsinki.fi/sanspanz/.

Studying the genetic basis of masting

The mechanisms underlying mast seeding have traditionally been studied by collecting long-term observational data on seed crops and correlating seedfall with environmental variables. Significant progress in ecological genomics will improve our understanding of the evolution of masting by clarifying the genetic basis of masting traits and the role of natural selection in shaping those traits. Here, we summarize three important aspects in studying the evolution of masting at the genetic level: which traits govern masting, whether those traits are genetically regulated, and which taxa show wide variation in these traits. We then introduce recent studies on the molecular mechanisms of masting. Those studies measure seasonal changes in gene expression in natural conditions to quantify how multiple environmental factors combine to regulate floral initiation, which in many masting plant species is the single largest contributor to among-year variation in seed crops. We show that Fagaceae offers exceptional opportunities for evolutionary investigations because of its diversity at both the phenotypic and genetic levels and existing documented genome sequences. This article is part of the theme issue 'The ecology and evolution of synchronized seed production in plants'.

Chromosome-scale genome assembly of Castanopsis tibetana provides a powerful comparative framework to study the evolution and adaptation of Fagaceae trees

Fagaceae species are increasingly used as models to elucidate the process and mechanism of adaptation and speciation by integrating ecology, evolution and genomics. The genus Castanopsis belongs to the family Fagaceae and is mainly distributed across subtropical and tropical Asia. In the present study, we reported the first chromosome-scale genome assembly of Castanopsis tibetana, a common species of evergreen broadleaved forests in subtropical China. The combination of Nanopore sequencing and Hi-C technologies enabled a high-quality genome assembly. The final assembled genome size of C. tibetana was 878.6 Mb (97.6% of the estimated genome size), consisting of 477 contigs with an N50 length of 3.3 Mb. The benchmarking universal single-copy orthologue (BUSCO) assessment indicated a completeness of 93.0%. Hi-C scaffolding generated 12 pseudochromosomes, representing 98.7% of the assembled genome. Subsequently, 40,937 protein-coding genes were predicted and 90.04% of them were functionally annotated. More than 476.9 Mb of repetitive sequences (54.3% of the genome) were identified, and the percentage of the genome covered by TE elements was 39.98%. Comparative genomics analysis revealed that C. tibetana was most closely related to Castanea mollissima and diverged at 18.48 Ma, and that C. tibetana has undergone considerable gene family expansion and contraction. Evidence of positive selection was detected in 53 genes, which showed different arrangement pattern compared to Quercus robur. The chromosome-scale genome assembly of C. tibetana will expand Fagaceae genome resources across the family and provide a powerful comparative framework to study the adaptation and evolution of Fagaceae trees.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

Comparative Transcriptome Analysis Revealed Candidate Genes Potentially Related to Desiccation Sensitivity of Recalcitrant Quercus variabilis Seeds

Chinese cork oak (Quercus variabilis) is a widely distributed and highly valuable deciduous broadleaf tree from both ecological and economic perspectives. Seeds of this species are recalcitrant, i.e., sensitive to desiccation, which affects their storage and long-term preservation of germplasm. However, little is known about the underlying molecular mechanism of desiccation sensitivity of Q. variabilis seeds. In this study, the seeds were desiccated with silica gel for certain days as different treatments from 0 (Control) to 15 days (T15) with a gradient of 1 day. According to the seed germination percentage, four key stages (Control, T2, T4, and T11) were found. Then the transcriptomic profiles of these four stages were compared. A total of 4,405, 4,441, and 5,907 differentially expressed genes (DEGs) were identified in T2 vs. Control, T4 vs. Control, and T11 vs. Control, respectively. Among them, 2,219 DEGs were overlapped in the three comparison groups. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these DEGs were enriched into 124 pathways, such as "Plant hormone signal transduction" and "Glycerophospholipid metabolism". DEGs related to hormone biosynthesis and signal transduction (ZEP, YUC, PYR, ABI5, ERF1B, etc.), stress response proteins (LEA D-29, HSP70, etc.), and phospholipase D (PLD1) were detected during desiccation. These genes and their interactions may determine the desiccation sensitivity of seeds. In addition, group specific DEGs were also identified in T2 vs. Control (PP2C62, UNE12, etc.), T4 vs. Control (WRKY1-like, WAK10, etc.), and T11 vs. Control (IBH1, bZIP44, etc.), respectively. Finally, a possible work model was proposed to show the molecular regulation mechanism of desiccation sensitivity in Q. variabilis seeds. This is the first report on the molecular regulation mechanism of desiccation sensitivity of Q. variabilis seeds using RNA-Seq. The findings could make a great contribution to seed storage and long-term conservation of recalcitrant seeds in the future.

MicroRNAs in Woody Plants

MicroRNAs (miRNAs) are small (∼21-nucleotides) non-coding RNAs found in plant and animals. MiRNAs function as critical post-transcriptional regulators of gene expression by binding to complementary sequences in their target mRNAs, leading to mRNA destabilization and translational inhibition. Plant miRNAs have some distinct characteristics compared to their animal counterparts, including greater evolutionary conservation and unique miRNA processing methods. The lifecycle of a plant begins with embryogenesis and progresses through seed germination, vegetative growth, reproductive growth, flowering and fruiting, and finally senescence and death. MiRNAs participate in the transformation of plant growth and development and directly monitor progression of these processes and the expression of certain morphological characteristics by regulating transcription factor genes involved in cell growth and differentiation. In woody plants, a large and rapidly increasing number of miRNAs have been identified, but their biological functions are largely unknown. In this review, we summarize the progress of miRNA research in woody plants to date. In particular, we discuss the potential roles of these miRNAs in growth, development, and biotic and abiotic stresses responses in woody plants.

Transcriptomic analysis of cork during seasonal growth highlights regulatory and developmental processes from phellogen to phellem formation

The phellogen or cork cambium stem cells that divide periclinally and outwardly specify phellem or cork. Despite the vital importance of phellem in protecting the radially-growing plant organs and wounded tissues, practically only the suberin biosynthetic process has been studied molecularly so far. Since cork oak (Quercus suber) phellogen is seasonally activated and its proliferation and specification to phellem cells is a continuous developmental process, the differentially expressed genes during the cork seasonal growth served us to identify molecular processes embracing from phellogen to mature differentiated phellem cell. At the beginning of cork growth (April), cell cycle regulation, meristem proliferation and maintenance and processes triggering cell differentiation were upregulated, showing an enrichment of phellogenic cells from which phellem cells are specified. Instead, at maximum (June) and advanced (July) cork growth, metabolic processes paralleling the phellem cell chemical composition, such as the biosynthesis of suberin, lignin, triterpenes and soluble aromatic compounds, were upregulated. Particularly in July, polysaccharides- and lignin-related secondary cell wall processes presented a maximal expression, indicating a cell wall reinforcement in the later stages of cork formation, presumably related with the initiation of latecork development. The putative function of relevant genes identified are discussed in the context of phellem ontogeny.

Proteomics Data Analysis for the Identification of Proteins and Derived Proteotypic Peptides of Potential Use as Putative Drought Tolerance Markers for Quercus ilex

Drought is one of the main causes of mortality in holm oak (Quercus ilex) seedlings used in reforestation programs. Although this species shows high adaptability to the extreme climate conditions prevailing in Southern Spain, its intrinsic genetic variability may play a role in the differential response of some populations and individuals. The aim of this work was to identify proteins and derived proteotypic peptides potentially useful as putative markers for drought tolerance in holm oak by using a targeted post-acquisition proteomics approach. For this purpose, we used a set of proteins identified by shotgun (LC-MSMS) analysis in a drought experiment on Q. ilex seedlings from four different provenances (viz. the Andalusian provinces Granada, Huelva, Cadiz and Seville). A double strategy involving the quantification of proteins and target peptides by shotgun analysis and post-acquisition data analysis based on proteotypic peptides was used. To this end, an initial list of proteotypic peptides from proteins highly represented under drought conditions was compiled that was used in combination with the raw files from the shotgun experiment to quantify the relative abundance of the fragment’s ion peaks with the software Skyline. The most abundant peptides under drought conditions in at least two populations were selected as putative markers of drought tolerance. A total of 30 proteins and 46 derived peptides belonging to the redox, stress-related, synthesis,-folding and degradation, and primary and secondary metabolism functional groups were thus identified. Two proteins (viz., subtilisin and chaperone GrpE protein) were found at increased levels in three populations, which make them especially interesting for validation drought tolerance markers in subsequent experiments.

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.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach

The pathological interaction between oak trees and Phytophthora cinnamomi has implications in the cork oak decline observed over the last decades in the Iberian Peninsula. During host colonization, the phytopathogen secretes effector molecules like elicitins to increase disease effectiveness. The objective of this study was to unravel the proteome changes associated with the cork oak immune response triggered by P. cinnamomi inoculation in a long-term assay, through SWATH-MS quantitative proteomics performed in the oak leaves. Using the Arabidopis proteome database as a reference, 424 proteins were confidently quantified in cork oak leaves, of which 80 proteins showed a p-value below 0.05 or a fold-change greater than 2 or less than 0.5 in their levels between inoculated and control samples being considered as altered. The inoculation of cork oak roots with P. cinnamomi increased the levels of proteins associated with protein-DNA complex assembly, lipid oxidation, response to endoplasmic reticulum stress, and pyridine-containing compound metabolic process in the leaves. In opposition, several proteins associated with cellular metabolic compound salvage and monosaccharide catabolic process had significantly decreased abundances. The most significant abundance variations were observed for the Ribulose 1,5-Bisphosphate Carboxylase small subunit (RBCS1A), Heat Shock protein 90–1 (Hsp90-1), Lipoxygenase 2 (LOX2) and Histone superfamily protein H3.3 (A8MRLO/At4G40030) revealing a pertinent role for these proteins in the host-pathogen interaction mechanism. This work represents the first SWATH-MS analysis performed in cork oak plants inoculated with P. cinnamomi and highlights host proteins that have a relevant action in the homeostatic states that emerge from the interaction between the oomycete and the host in the long term and in a distal organ.

CorkOakDB-The Cork Oak Genome Database Portal

Quercus suber (cork oak) is an evergreen tree native to the Mediterranean basin, which plays a key role in the ecology and economy of this area. Over the last decades, this species has gone through an observable decline, mostly due to environmental factors. Deciphering the mechanisms of cork oak's response to the environment and getting a deep insight into its biology are crucial to counteract biotic and abiotic stresses compromising the stability of a unique ecosystem. In the light of these setbacks, the publication of the genome in 2018 was a major step towards understanding the genetic make-up of this species. In an effort to integrate this information in a comprehensive, accessible and intuitive format, we have developed The Cork Oak Genome Database Portal (CorkOakDB). The CorkOakDB is supported by the BioData.pt e-infrastructure, the Portuguese ELIXIR node for biological data. The portal gives public access to search and explore the curated genomic and transcriptomic data on this species. Moreover, CorkOakDB provides a user-friendly interface and functional tools to help the research community take advantage of the increased accessibility to genomic information. A study case is provided to highlight the functionalities of the portal. CorkOakDB guarantees the update, curation and data collection, aiming to collect data besides the genetic/genomic information, in order to become the main repository in cork oak research. Database URL: http://corkoakdb.org/.

Adaptive Evolution of Chalcone Isomerase Superfamily in Fagaceae

Chalcone Isomerase (CHI) catalyzes the biosynthesis of flavonoids and secondary metabolism in plants. Currently, there is no systematic analysis of CHIs gene family in Fagaceae which is available. In this study, twenty-two CHI proteins were identified in five species of the Fagaceae family. The CHI superfamily in Fagaceae can be classified into three subfamilies and five groups using phylogenetic analysis, analysis of physicochemical properties, and structural prediction. Results indicated that serine (Ser) and isoleucine (Ile) residues determine the substrate preferred by active Type I Fagaceae CHI, and the chalcone isomerase-like (CHIL) of Fagaceae had active site residues. Adaptive analysis of CHIs showed that CHIs are subject to selection pressure. The active CHI gene of Fagaceae was located in the cytoplasm, and it had the typical gene structure of CHI and contains four exons. All the twenty-two identified CHIs had the conserved domain motif 3, and the different groups had their own structural characteristics. In the process of fatty acid binding protein (FAP) evolution to CHIL and CHI, the physical and chemical properties of proteins also had significant differences in addition to changes in protein functions.

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.

Oxidosqualene cyclases involved in the biosynthesis of triterpenoids in Quercus suber cork

Cork is a water-impermeable, suberin-based material harboring lignin, (hemi)cellulose, and extractable small molecules (primarily triterpenoids). Extractables strongly influence the properties of suberin-based materials. Though these previous findings suggest a key role for triterpenoids in cork material quality, directly testing this idea is hindered in part because it is not known which genes control cork triterpenoid biosynthesis. Here, we used gas chromatography and mass spectrometry to determine that the majority (>85%) of non-polar extractables from cork were pentacyclic triterpenoids, primarily betulinic acid, friedelin, and hydroxy-friedelin. In other plants, triterpenoids are generated by oxidosqualene cyclases (OSCs). Accordingly, we mined Quercus suber EST libraries for OSC fragments to use in a RACE PCR-based approach and cloned three full-length OSC transcripts from cork (QsOSC1-3). Heterologous expression in Saccharomyces cerevisiae revealed that QsOSC1-3 respectively encoded enzymes with lupeol synthase, mixed α- and β-amyrin synthase, and mixed β-amyrin and friedelin synthase activities. These activities together account for the backbone structures of the major cork triterpenoids. Finally, we analyzed the sequences of QsOSC1-3 and other plant OSCs to identify residues associated with specific OSC activities, then combined this with analyses of Q. suber transcriptomic and genomic data to evaluate potential redundancies in cork triterpenoid biosynthesis.