Phellem versus xylem: genome-wide transcriptomic analysis reveals novel regulators of cork formation in cork oak

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.

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.

Comprehensive analysis of lncRNA-mRNA regulatory network in Populus associated with xylem development

Long noncoding RNAs (lncRNAs) play essential roles in numerous biological processes in plants, such as regulating the gene expression. However, only a few studies have looked into their potential functions in xylem development. High-throughput sequencing of P. euramericana 'Zhonglin46' developing and mature xylem was performed in this study. Through sequencing analysis, 14,028 putative lncRNA transcripts were identified, including 4525 differentially expressed lncRNAs (DELs). Additional research revealed that in mature xylem, a total of 2320 DELs were upregulated and 2205 were downregulated compared to developing xylem. Meanwhile, there were a total of 8122 differentially expressed mRNAs (DEMs) that were upregulated and 16,424 that were downregulated in mature xylem compared with developing xylem. The cis- and trans-target genes of DELs were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, which indicated that these DELs participate in controlling the phenylpropanoid and lignin biosynthesis pathway as well as the starch and sucrose metabolism pathway. Among the cis-regulated DELs, LNC_006291, LNC_006292, and LNC_006532 all participate in regulating multiple HCT gene family membranes. As targets, POPTR_001G045900v3 (CCR2) and POPTR_018G063500v3 (SUS) both have only one cis-regulatory lncRNA, referred to as LNC_000057 and LNC_006212, respectively. Moreover, LNC_004484 and two DELs named LNC_008014 and LNC_010781 were revealed to be important nodes in the co-expression network of trans-lncRNAs and mRNAs associated to the lignin biosynthesis pathway and cellulose and xylan biosynthetic pathways, respectively. Finally, quantitative real-time PCR (qRT-PCR) was used to confirme 34 pairs of lncRNA-mRNA. Taken together, these findings may help to clarify the regulatory role that lncRNAs play in xylem development and wood formation.

Priorities for Bark Anatomical Research: Study Venues and Open Questions

The bark fulfils several essential functions in vascular plants and yields a wealth of raw materials, but the understanding of bark structure and function strongly lags behind our knowledge with respect to other plant tissues. The recent technological advances in sampling and preparation of barks for anatomical studies, along with the establishment of an agreed bark terminology, paved the way for more bark anatomical research. Whilst datasets reveal bark's taxonomic and functional diversity in various ecosystems, a better understanding of the bark can advance the understanding of plants' physiological and environmental challenges and solutions. We propose a set of priorities for understanding and further developing bark anatomical studies, including periderm structure in woody plants, phloem phenology, methods in bark anatomy research, bark functional ecology, relationships between bark macroscopic appearance, and its microscopic structure and discuss how to achieve these ambitious goals.

Cork Development: What Lies Within

The cork layer present in all dicotyledonous plant species with radial growth is the result of the phellogen activity, a secondary meristem that produces phellem (cork) to the outside and phelloderm inwards. These three different tissues form the periderm, an efficient protective tissue working as a barrier against external factors such as environmental aggressions and pathogen attacks. The protective function offered by cork cells is mainly due to the abundance of suberin in their cell walls. Chemically, suberin is a complex aliphatic network of long chain fatty acids and alcohols with glycerol together with aromatic units. In most woody species growing in temperate climates, the first periderm is replaced by a new functional periderm upon a few years after being formed. One exception to this bark development can be found in cork oak (Quercus suber) which display a single periderm that grows continuously. Quercus suber stands by its thick cork layer development with continuous seasonal growth. Cork raw material has been exploited by man for centuries, especially in Portugal and Spain. Nowadays, its applications have widened vastly, from the most known product, stoppers, to purses or insulating materials used in so many industries, such as construction and car production. Research on how cork develops, and the effect environmental factors on cork oak trees is extremely important to maintain production of good-quality cork, and, by maintaining cork oak stands wealthy, we are preserving a very important ecosystem both by its biodiversity and its vital social and economic role in areas already showing a population declination.

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.

Populus × canescens root suberization in reaction to osmotic and salt stress is limited to the developing younger root tip region

Populus is a valuable and fast-growing tree species commonly cultivated for economic and scientific purposes. But most of the poplar species are sensitive to drought and salt stress. Thus, we compared the physiological effects of osmotic stress (PEG8000) and salt treatment (NaCl) on poplar roots to identify potential strategies for future breeding or genetic engineering approaches. We investigated root anatomy using epifluorescence microscopy, changes in root suberin composition and amount using gas chromatography, transcriptional reprogramming using RNA sequencing, and modifications of root transport physiology using a pressure chamber. Poplar roots reacted to the imposed stress conditions, especially in the developing younger root tip region, with remarkable differences between both types of stress. Overall, the increase in suberin content was surprisingly small, but the expression of key suberin biosynthesis genes was strongly induced. Significant reductions of the radial water transport in roots were only observed for the osmotic and not the hydrostatic hydraulic conductivity. Our data indicate that the genetic enhancement of root suberization processes in poplar might be a promising target to convey increased tolerance, especially against toxic sodium chloride.

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.

The Making of Plant Armor: The Periderm

The periderm acts as armor protecting the plant's inner tissues from biotic and abiotic stress. It forms during the radial thickening of plant organs such as stems and roots and replaces the function of primary protective tissues such as the epidermis and the endodermis. A wound periderm also forms to heal and protect injured tissues. The periderm comprises a meristematic tissue called the phellogen, or cork cambium, and its derivatives: the lignosuberized phellem and the phelloderm. Research on the periderm has mainly focused on the chemical composition of the phellem due to its relevance as a raw material for industrial processes. Today, there is increasing interest in the regulatory network underlying periderm development as a novel breeding trait to improve plant resilience and to sequester CO₂. Here, we discuss our current understanding of periderm formation, focusing on aspects of periderm evolution, mechanisms of periderm ontogenesis, regulatory networks underlying phellogen initiation and cork differentiation, and future challenges of periderm research.

Translational profile of developing phellem cells in Arabidopsis thaliana roots

The phellem is a specialized boundary tissue providing the first line of defense against abiotic and biotic stresses in organs undergoing secondary growth. Phellem cells undergo several differentiation steps, which include cell wall suberization, cell expansion, and programmed cell death. Yet, the molecular players acting particularly in phellem cell differentiation remain poorly described, particularly in the widely used model plant Arabidopsis thaliana. Using specific marker lines we followed the onset and progression of phellem differentiation in A. thaliana roots and further targeted the translatome of newly developed phellem cells using translating ribosome affinity purification followed by mRNA sequencing (TRAP-SEQ). We showed that phellem suberization is initiated early after phellogen (cork cambium) division. The specific translational landscape was organized in three main domains related to energy production, synthesis and transport of cell wall components, and response to stimulus. Novel players in phellem differentiation related to suberin monomer transport and assembly as well as novel transcription regulators were identified. This strategy provided an unprecedented resolution of the translatome of developing phellem cells, giving a detailed and specific view on the molecular mechanisms acting on cell differentiation in periderm tissues of the model plant Arabidopsis.

Cork cells in cork oak periderms undergo programmed cell death and proanthocyanidin deposition

Vascular plants with secondary growth develop a periderm mostly composed of dead suberized cork cells to face environmental hostile conditions. Cork oak has a highly active and long-living phellogen forming a remarkably thick periderm that is periodically debarked for industrial purposes. This wounding originates the quick formation of a new traumatic periderm, making cork oak an exceptional model to study the first periderm differentiation during normal development in young sprigs and traumatic (wound) periderm formation after debarking. Here, we studied the poorly known first periderm differentiation steps that involve cell wall suberization, polyphenolic accumulation and programmed cell death (PCD) by combining transmission electron microscopy, histochemical and molecular methods in periderms from young sprigs. These processes were further compared with traumatic periderms formed after wounding using molecular and histochemical techniques, such as the polyphenolic accumulation. In the first periderms from young sprigs, four distinct differentiation stages were defined according to the presence of PCD morphological features. First young and traumatic periderms showed an upregulation of genes related to suberin biosynthesis, proanthocyanidins biosynthesis and transport, autophagy, and PCD. Traumatic periderms revealed an overall upregulation of these genes, likely resulting from ontogeny differences and distinct phellogen origin associated with a faster metabolism, highlighting the impact of wounding on phellogen activity after debarking. First periderms from young sprigs showed gradual accumulation of proanthocyanidins in the vacuoles throughout PCD stages until total filled lumens, whereas in traumatic periderms, these compounds were found cell wall linked in already empty cells. This work enabled a comprehensive overview of the cork cells differentiation processes contributing to deepening the knowledge of the fundamental ontogenic program of this protective tissue, which is also a unique forest product, constituting the basis of a sustainable and profitable industry.

Silencing of StRIK in potato suggests a role in periderm related to RNA processing and stress

BACKGROUND

The periderm is a protective barrier crucial for land plant survival, but little is known about genetic factors involved in its development and regulation. Using a transcriptomic approach in the cork oak (Q. suber) periderm, we previously identified an RS2-INTERACTING KH PROTEIN (RIK) homologue of unknown function containing a K homology (KH)-domain RNA-binding protein, as a regulatory candidate gene in the periderm.

RESULTS

To gain insight into the function of RIK in the periderm, potato (S. tuberosum) tuber periderm was used as a model: the full-length coding sequence of RIK, hereafter referred to as StRIK, was isolated, the transcript profile analyzed and gene silencing in potato performed to analyze the silencing effects on periderm anatomy and transcriptome. The StRIK transcript accumulated in all vegetative tissues studied, including periderm and other suberized tissues such as root and also in wounded tissues. Downregulation of StRIK in potato by RNA interference (StRIK-RNAi) did not show any obvious effects on tuber periderm anatomy but, unlike Wild type, transgenic plants flowered. Global transcript profiling of the StRIK-RNAi periderm did show altered expression of genes associated with RNA metabolism, stress and signaling, mirroring the biological processes found enriched within the in silico co-expression network of the Arabidopsis orthologue.

CONCLUSIONS

The ubiquitous expression of StRIK transcript, the flower associated phenotype and the differential expression of StRIK-RNAi periderm point out to a general regulatory role of StRIK in diverse plant developmental processes. The transcriptome analysis suggests that StRIK might play roles in RNA maturation and stress response in the periderm.

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.

Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.

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