Characterization of the cork oak transcriptome dynamics during acorn development

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.

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.

Quercus suber Roots Activate Antioxidant and Membrane Protective Processes in Response to High Salinity

Cork oak (Quercus suber) is a species native to Mediterranean areas and its adaptation to the increasingly prevalent abiotic stresses, such as soil salinization, remain unknown. In sequence with recent studies on salt stress response in the leaf, it is fundamental to uncover the plasticity of roots directly exposed to high salinity to better understand how Q. suber copes with salt stress. In the present study we aimed to unveil the antioxidants and key-genes involved in the stress-responses (early vs. later responses) of Q. suber roots exposed to high salinity. Two-month-old Q. suber plants were watered with 300 mM NaCl solution and enzymatic and non-enzymatic antioxidants, lipid peroxidation and the relative expression of genes related to stress response were analysed 8 h and 6 days after salt treatment. After an 8 h of exposure, roots activated the expression of QsLTI30 and QsFAD7 genes involved in stress membrane protection, and QsRAV1 and QsCZF1 genes involved in tolerance and adaptation. As a result of the continued salinity stress (6 days), lipid peroxidation increased, which was associated with an upregulation of QsLTI30 gene. Moreover, other protective mechanisms were activated, such as the upregulation of genes related to antioxidant status, QsCSD1 and QsAPX2, and the increase of the antioxidant enzyme activities of superoxide dismutase, catalase, and ascorbate peroxidase, concomitantly with total antioxidant activity and phenols. These data suggest a response dependent on the time of salinity exposure, leading Q. suber roots to adopt protective complementary strategies to deal with salt stress.

The Relevance of a Physiological-Stage Approach Study of the Molecular and Environmental Factors Regulating Seed Germination in Wild Plants

Germination represents the culmination of the seed developmental program and is affected by the conditions prevailing during seed maturation in the mother plant. During maturation, the dormancy condition and tolerance to dehydration are established. These characteristics are modulated by the environment to which they are subjected, having an important impact on wild species. In this work, a review was made of the molecular bases of the maturation, the processes of dormancy imposition and loss, as well as the germination process in different wild species with different life histories, and from diverse habitats. It is also specified which of these species present a certain type of management. The impact that the domestication process has had on certain characteristics of the seed is discussed, as well as the importance of determining physiological stages based on morphological characteristics, to face the complexities of the study of these species and preserve their genetic diversity and physiological responses.

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

Roles of metabolic regulation in developing Quercus variabilis acorns at contrasting geologically-derived phosphorus sites in subtropical China

BACKGROUND

Phosphorus (P) -rich soils develop in phosphorite residing areas while P-deficient soils are ubiquitous in subtropical regions. Little has been reported that how metabolites participate in the seed development and the processes involved in their coping with contrasting-nutrient environments.

RESULTS

Here we quantified the metabolites of Quercus variabilis acorns in the early (July), middle (August), late (September) development stages, and determined element (C, H, O, N, P, K, Ca, Mg, S, Fe, Al, Mn, Na, Zn, and Cu) concentrations of acorns in the late stage, at geologically-derived contrasting-P sites in subtropical China. The primary metabolic pathways included sugar metabolism, the TCA cycle, and amino acid metabolism. Most metabolites (especially C- and N-containing metabolites) increased and then decreased from July to September. Acorns between the two sites were significantly discriminated at the three stages, respectively, by metabolites (predominantly sugars and organic acids). Concentrations of P, orthophosphoric acid and most sugars were higher; erythrose was lower in late-stage acorns at P-rich sites than those at P-deficient sites. No significant differences existed in the size and dry mass of individual acorns between oak populations at the two sites.

CONCLUSIONS

Oak acorns at the two sites formed distinct metabolic phenotypes related to their distinct geologically-derived soil conditions, and the late-stage acorns tended to increase P-use-efficiency in the material synthesis process at P-deficient sites, relative to those at P-rich sites.

Dissecting the Seed Maturation and Germination Processes in the Non-Orthodox Quercus ilex Species Based on Protein Signatures as Revealed by 2-DE Coupled to MALDI-TOF/TOF Proteomics Strategy

Unlike orthodox species, seed recalcitrance is poorly understood, especially at the molecular level. In this regard, seed maturation and germination were studied in the non-orthodox Quercus ilex by using a proteomics strategy based on two-dimensional gel electrophoresis coupled to matrix-assisted laser desorption ionization/time of flight (2-DE-MALDI-TOF).Cotyledons and embryo/radicle were sampled at different developmental stages, including early (M1–M3), middle (M4–M7), and late (M8–M9) seed maturation, and early (G1–G3) and late (G4–G5) germination. Samples corresponding to non-germinating, inviable, seeds were also included. Protein extracts were subjected to 2-dimensional gel electrophoresis (2-DE) and changes in the protein profiles were analyzed. Identified variable proteins were grouped according to their function, being the energy, carbohydrate, lipid, and amino acid metabolisms, together with protein fate, redox homeostasis, and response to stress are the most represented groups. Beyond the visual aspect, morphometry, weight, and water content, each stage had a specific protein signature. Clear tendencies for the different protein groups throughout the maturation and germination stages were observed for, respectively, cotyledon and the embryo axis. Proteins related to metabolism, translation, legumins, proteases, proteasome, and those stress related were less abundant in non-germinating seeds, it related to the loss of viability. Cotyledons were enriched with reserve proteins and protein-degrading enzymes, while the embryo axis was enriched with proteins of cell defense and rescue, including heat-shock proteins (HSPs) and antioxidants. The peaks of enzyme proteins occurred at the middle stages (M6–M7) in cotyledons and at late ones (M8–M9) in the embryo axis. Unlike orthodox seeds, proteins associated with glycolysis, tricarboxylic acid cycle, carbohydrate, amino acid and lipid metabolism are present at high levels in the mature seed and were maintained throughout the germination stages. The lack of desiccation tolerance in Q. ilex seeds may be associated with the repression of some genes, late embryogenesis abundant proteins being one of the candidates.

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.

In silico and expression analyses of fasciclin-like arabinogalactan proteins reveal functional conservation during embryo and seed development

KEY MESSAGE

The fasciclin-like arabinogalactan proteins organization into four groups is conserved and may be related to specific roles in developmental processes across angiosperms. Fasciclin-like arabinogalactan proteins (FLAs) are a subclass of arabinogalactan proteins (AGPs), which contain fasciclin-like domains in addition to typical AGP domains. FLAs are present across all embryophytes, and despite their low overall sequence similarity, conserved regions that define the fasciclin functional domain (FAS) have been identified, suggesting that the cell adhesion property is also conserved. FLAs in Arabidopsis have been organized into four subgroups according to the number and distribution of functional domains. Recent studies associated FLAs with cell wall-related processes where domain organization seemed to be related to functional roles. In Arabidopsis, FLAs containing a single FAS domain were found to be important for the integrity and elasticity of the plant cell wall matrix, and FLAs with two FAS domains and two AGP domains were found to be involved in maintaining proper cell expansion under salt stress conditions. The main purpose of the present work was to elucidate the expression pattern of selected FLA genes during embryo and seed development using RT-qPCR. AtFLA8 and AtFLA10, two Arabidopsis genes that stood out in previous microarray studies of embryo development, were further examined using promoter-driven gene reporter analyses. We also studied the expression of cork oak FLA genes and found that their expression partially parallels the expression patterns of the putative AtFLA orthologs. We propose that the functional organization of FLAs is conserved and may be related to fundamental aspects of embryogenesis and seed development across angiosperms. Phylogenetic studies were performed, and we show that the same basic four-subgroup organization described for Arabidopsis FLA gene classification is valid for most Arabidopsis FLA orthologs of several plant species, namely poplar, corn and cork oak.

Correlated metabolic and elemental variations between the leaves and seeds of oak trees at contrasting geologically derived phosphorus sites

The leaves and seeds of plants frequently function as the source and sink organs for distinct metabolites, which can interactively vary in response to adverse site conditions. Subtropical soils are typically characterized as having deficient phosphorus (P), calcium (Ca), and magnesium (Mg), with enriched aluminum (Al) and iron (Fe), while Al and manganese (Mn) are toxic at low pH. It remains largely unknown how leaf- and seed-sourced metabolites are synergistically linked to adapt to P-variable soils for trees in subtropical areas. Here we quantified the metabolic and elemental profiling in the mature leaves and immature seeds of Quercus variabilis at contrasting geologically-derived phosphorus sites in subtropical China. The results revealed that carbon (C) and nitrogen (N) based metabolites (primarily sugars and organic acids), as well as enzyme- and protein/nucleic acid-related elements (N, P, Mg, and Mn) played important roles toward characterizing the profiling of metabolites and ionomes in leaves and seeds at two site types, respectively. These metabolites (sugars, amino acids, and fatty acids) and elements (N, P, Mg, and Mn) of seeds were closely related to the sugars, organic acids, and elements (N, P, Mg, and Mn) of leaves at the two site types. For the most part, the content of N and P in the soil affected the accumulation of materials (such as, starchs and proteins) in seeds, as well as N and P assimilation in leaves, by influencing C- and N-containing metabolites in leaves. These results suggested that correlated disparities of C- and N-containing metabolites, along with enzyme- and protein/nucleic acid-related elements in both leaves and seeds played important roles in plants to facilitate their adaptation to nutrient-variable sites in subtropical zones.

Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos

BACKGROUND

There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis.

RESULTS

Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons.

CONCLUSIONS

Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.

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.

Characterization of the cork formation and production transcriptome in Quercus cerris × suber hybrids

Cork oak is the main cork-producing species worldwide, and plays a significant economic, ecological and social role in the Mediterranean countries, in particular in Portugal and Spain. The ability to produce cork is limited to a few species, hence it must involve specific regulation mechanisms that are unique to these species. However, to date, these mechanisms remain largely understudied, especially with approaches involving the use of high-throughput sequencing technology. In this study, the transcriptome of cork-producing and non-cork-producing Quercus cerris × suber hybrids was analyzed in order to elucidate the differences between the two groups of trees displaying contrasting phenotypes for cork production. The results revealed the presence of a significant number of genes exclusively associated with cork production, in the trees that developed cork. Moreover, several gene ontology subcategories, such as cell wall biogenesis, lipid metabolic processes, metal ion binding and apoplast/cell wall, were only detected in the trees with cork production. These results indicate the existence, at the transcriptome level, of mechanisms that seem to be unique and necessary for cork production, which is an advancement in our knowledge regarding the genetic regulation behind cork formation and production.

Transcriptome characterization and screening of molecular markers in ecologically important Himalayan species (Rhododendron arboreum)

Rhododendron arboreum is an ecologically prominent species, which also lends commercial and medicinal benefits in the form of palatable juices and useful herbal drugs. Local abundance and survival of the species under a highly fluctuating climate make it an ideal model for genetic structure and functional analysis. However, a lack of genomic data has hampered additional research. In the present study, cDNA libraries from floral and foliar tissues of the species were sequenced to provide a foundation for understanding the functional aspects of the genome and to construct an enriched repository that will promote genomics studies in the genera. Illumina’s platform facilitated the generation of ∼100 million high-quality paired-end reads. De novo assembly, clustering, and filtering out of shorter transcripts predicted 113 167 non-redundant transcripts with an average length of 1164.6 bases. Of these, 71 961 transcripts were categorized based on functional annotations in the Gene Ontology database, whereby 5710 were grouped into 141 pathways and 23 746 encoded for different transcription factors. Transcriptome screening further identified 35 419 microsatellite regions, of which, 43 polymorphic loci were characterized on 30 genotypes. Seven hundred and nineteen transcripts had 811 high-quality single-nucleotide polymorphic variants with a minimum coverage of 10, a total score of 20, and SNP% of 50.

The draft genome sequence of cork oak

Cork oak (Quercus suber) is native to southwest Europe and northwest Africa where it plays a crucial environmental and economical role. To tackle the cork oak production and industrial challenges, advanced research is imperative but dependent on the availability of a sequenced genome. To address this, we produced the first draft version of the cork oak genome. We followed a de novo assembly strategy based on high-throughput sequence data, which generated a draft genome comprising 23,347 scaffolds and 953.3 Mb in size. A total of 79,752 genes and 83,814 transcripts were predicted, including 33,658 high-confidence genes. An InterPro signature assignment was detected for 69,218 transcripts, which represented 82.6% of the total. Validation studies demonstrated the genome assembly and annotation completeness and highlighted the usefulness of the draft genome for read mapping of high-throughput sequence data generated using different protocols. All data generated is available through the public databases where it was deposited, being therefore ready to use by the academic and industry communities working on cork oak and/or related species.

De novo assembly, functional annotation, and analysis of the giant reed (Arundo donax L.) leaf transcriptome provide tools for the development of a biofuel feedstock

BACKGROUND

Arundo donax has attracted renewed interest as a potential candidate energy crop for use in biomass-to-liquid fuel conversion processes and biorefineries. This is due to its high productivity, adaptability to marginal land conditions, and suitability for biofuel and biomaterial production. Despite its importance, the genomic resources currently available for supporting the improvement of this species are still limited.

RESULTS

We used RNA sequencing (RNA-Seq) to de novo assemble and characterize the A. donax leaf transcriptome. The sequencing generated 1249 million clean reads that were assembled using single-k-mer and multi-k-mer approaches into 62,596 unique sequences (unitranscripts) with an N50 of 1134 bp. TransDecoder and Trinotate software suites were used to obtain putative coding sequences and annotate them by mapping to UniProtKB/Swiss-Prot and UniRef90 databases, searching for known transcripts, proteins, protein domains, and signal peptides. Furthermore, the unitranscripts were annotated by mapping them to the NCBI non-redundant, GO and KEGG pathway databases using Blast2GO. The transcriptome was also characterized by BLAST searches to investigate homologous transcripts of key genes involved in important metabolic pathways, such as lignin, cellulose, purine, and thiamine biosynthesis and carbon fixation. Moreover, a set of homologous transcripts of key genes involved in stomatal development and of genes coding for stress-associated proteins (SAPs) were identified. Additionally, 8364 simple sequence repeat (SSR) markers were identified and surveyed. SSRs appeared more abundant in non-coding regions (63.18%) than in coding regions (36.82%). This SSR dataset represents the first marker catalogue of A. donax. 53 SSRs (PolySSRs) were then predicted to be polymorphic between ecotype-specific assemblies, suggesting genetic variability in the studied ecotypes.

CONCLUSIONS

This study provides the first publicly available leaf transcriptome for the A. donax bioenergy crop. The functional annotation and characterization of the transcriptome will be highly useful for providing insight into the molecular mechanisms underlying its extreme adaptability. The identification of homologous transcripts involved in key metabolic pathways offers a platform for directing future efforts in genetic improvement of this species. Finally, the identified SSRs will facilitate the harnessing of untapped genetic diversity. This transcriptome should be of value to ongoing functional genomics and genetic studies in this crop of paramount economic importance.