Transcriptional profiling of cork oak phellogenic cells isolated by laser microdissection

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.

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.

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.

Potato Periderm Development and Tuber Skin Quality

The periderm is a corky tissue that replaces the epidermis when the latter is damaged, and is critical for preventing pathogen invasion and water loss. The periderm is formed through the meristematic activity of phellogen cells (cork cambium). The potato skin (phellem cells) composes the outer layers of the tuber periderm and is a model for studying cork development. Early in tuber development and following tuber expansion, the phellogen becomes active and produces the skin. New skin layers are continuously added by division of the phellogen cells until tuber maturation. Some physiological disorders of the potato tuber are related to abnormal development of the skin, including skinning injuries and russeting of smooth-skinned potatoes. Thus, characterizing the potato periderm contributes to modeling cork development in plants and helps to resolve critical agricultural problems. Here, we summarize the data available on potato periderm formation, highlighting tissue characteristics rather than the suberization processes.

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.

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.

Laser Microdissection of Specific Stem-Base Tissue Types from Olive Microcuttings for Isolation of High-Quality RNA

Simple Summary

Only a small portion of the stem cells participate in the process of adventitious root formation and the cells/tissues types involved in this process is species-dependent. In olive, it is still unclear which type of cells acquire competence for rooting. Regardless, the entire stem nodal segment (containing a mixture of distinct cell types) continues to be used in studies related to the molecular mechanisms underlying this process. Laser microdissection (LM) technology has been applied to isolate specific tissue and cell types. However, it is difficult to find a standard LM protocol suitable for all plant species and cell types and, thus, LM procedures must be developed and optimized for each particular tissue. In this study, we aimed to evaluate the efficiency of a LM protocol in olive microcuttings stem-base samples. This work presents a simple, rapid and efficient LM procedure for harvesting specific tissue types used for further high-quality RNA isolation. This will encourage future cell type-specific transcriptomic studies, contributing at deciphering rooting-competent cells in olive stems and to better understand the molecular mechanisms underlying the process of adventitious root formation.

Abstract

Higher plants are composed of different tissue and cell types. Distinct cells host different biochemical and physiological processes which is reflected in differences in gene expression profiles, protein and metabolite levels. When omics are to be carried out, the information provided by a specific cell type can be diluted and/or masked when using a mixture of distinct cells. Thus, studies performed at the cell- and tissue-type level are gaining increasing interest. Laser microdissection (LM) technology has been used to isolate specific tissue and cell types. However, this technology faces some challenges depending on the plant species and tissue type under analysis. Here, we show for the first time a LM protocol that proved to be efficient for harvesting specific tissue types (phloem, cortex and epidermis) from olive stem nodal segments and obtaining RNA of high quality. This is important for future transcriptomic studies to identify rooting-competent cells. Here, nodal segments were flash-frozen in liquid nitrogen-cooled isopentane and cryosectioned. Albeit the lack of any fixatives used to preserve samples’ anatomy, cryosectioned sections showed tissues with high morphological integrity which was comparable with that obtained with the paraffin-embedding method. Cells from the phloem, cortex and epidermis could be easily distinguished and efficiently harvested by LM. Total RNA isolated from these tissues exhibited high quality with RNA Quality Numbers (determined by a Fragment Analyzer System) ranging between 8.1 and 9.9. This work presents a simple, rapid and efficient LM procedure for harvesting specific tissue types of olive stems and obtaining high-quality RNA.

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.

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

The sugarcane ShMYB78 transcription factor activates suberin biosynthesis in Nicotiana benthamiana

KEY MESSAGE

A sugarcane MYB present in the culm induces suberin biosynthesis and is involved both with fatty acid and phenolics metabolism. Few transcription factors have been described as regulators of cell wall polymers deposition in C4 grasses. Particularly, regulation of suberin biosynthesis in this group of plants remains poorly understood. Here, we showed that the sugarcane MYB transcription factor ShMYB78 is an activator of suberin biosynthesis and deposition. ShMYB78 was identified upon screening genes whose expression was upregulated in sugarcane internodes undergoing suberization during culm development or triggered by wounding. Agrobacterium-mediated transient expression of ShMYB78 in Nicotiana benthamiana leaves induced the ectopic deposition of suberin and its aliphatic and aromatic monomers. Further, the expression of suberin-related genes was induced by ShMYB78 heterologous expression in Nicotiana benthamiana leaves. ShMYB78 was shown to be a nuclear protein based on its presence in sugarcane internode nuclear protein extracts, and protoplast transactivation assays demonstrated that ShMYB78 activates the promoters of the sugarcane suberin biosynthetic genes β-ketoacyl-CoA synthase (ShKCS20) and caffeic acid-O-methyltransferase (ShCOMT). Our results suggest that ShMYB78 may be involved in the transcriptional regulation of suberin deposition, from fatty acid metabolism to phenylpropanoid biosynthesis, in sugarcane internodes.

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.

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

Cork cambium (or phellogen) is a secondary meristem responsible for the formation of phelloderm and phellem/cork, which together compose the periderm. In Quercus suber L., the phellogen is active throughout the entire life of the tree, producing a continuous and renewable outer bark of cork. To identify specific candidate genes associated with cork cambium activity and phellem differentiation, we performed a comparative transcriptomic study of Q. suber secondary growth tissues (xylem and phellogen/phellem) using RNA-seq. The present work provides a high-resolution map of all the transcripts identified in the phellogen/phellem tissues. A total of 6013 differentially expressed genes were identified, with 2875 of the transcripts being specifically enriched during the cork formation process versus secondary xylem formation. Furthermore, cork samples originating from the original phellogen (`virgin' cork) and from a traumatic phellogen (`amadia' cork) were also compared. Our results point to a shortlist of potentially relevant candidate genes regulating phellogen activity and phellem differentiation, including novel genes involved in the suberization process, as well as genes associated to ethylene and jasmonate signaling and to meristem function. The future functional characterization of some of the identified candidate genes will help to elucidate the molecular mechanisms underlying cork cambium activity and phellem differentiation.

The development of the periderm: the final frontier between a plant and its environment

The periderm acts as the first line of defence for a plant, protecting wood and phloem from abiotic and biotic stresses. During secondary growth, through the increase in girth of plant organs, the periderm replaces the epidermis as the outermost tissue. The phellogen, a bifacial post-embryonic meristem, forms the phelloderm inwards (toward the vasculature) and the suberized phellem outwards (toward the environment). These three tissues are collectively referred to as the periderm. Here, we summarize recent findings on the molecular mechanisms of periderm development by describing periderm formation in connection to the fate of the surrounding tissues, by discussing common regulatory hubs between the vascular cambium and the phellogen, and by highlighting transcription factors (TFs) controlling phellem differentiation.

Phellem Cell-Wall Components Are Discriminants of Cork Quality in Quercus suber

Cork is a renewable, non-wood high valued forest product, with relevant ecological and economic impact in the Mediterranean-type ecosystems. Currently, cork is ranked according to its commercial quality. The most valuable planks are chosen for cork stoppers production. Cork planks with adequate thickness and porosity are classified as stoppable quality cork (SQC). The chemical composition of cork is known, but the regulation of metabolic pathways responsible of cork production and composition, hence of cork quality, is largely unknown. Here, we tested the hypothesis that post-genomic events may be responsible for the development of SQC and N-SQC (non-stoppable quality cork). Here, we show that combined proteomics and targeted metabolomics (namely soluble and cell wall bound phenolics) analyzed on recently formed phellem allows discriminate cork planks of different quality. Phellem cells of SQC and N-SQC displayed different reducing capacity, with consequential impact on both enzymatic pathways (e.g., glycolysis) and other cellular functions, including cell wall assembly and suberization. Glycolysis and respiration related proteins were abundant in both cork quality groups, whereas the level of several proteins associated to mitochondrial metabolism was higher in N-SQC. The soluble and cell wall-bound phenolics in recently formed phellem clearly discriminated SQC from N-SCQ. In our study, SQC was characterized by a high incorporation of aromatic components of the phenylpropanoid pathway in the cell wall, together with a lower content of hydrolysable tannins. Here, we propose that the level of hydrolysable tannins may represent a valuable diagnostic tool for screening recently formed phellem, and used as a proxy for the quality grade of cork plank produced by each tree.

The transcriptome of potato tuber phellogen reveals cellular functions of cork cambium and genes involved in periderm formation and maturation

The periderm is a protective corky tissue that is formed through the cambial activity of phellogen cells, when the outer epidermis is damaged. Timely periderm formation is critical to prevent pathogen invasion and water loss. The outer layers of the potato periderm, the tuber skin, serves as a model to study cork development. Early in tuber development the phellogen becomes active and produces the skin. During tuber maturation it becomes inactive and the skin adheres to the tuber flesh. The characterization of potato phellogen may contribute to the management of costly agricultural problems related to incomplete skin-set and the resulting skinning injuries, and provide us with new knowledge regarding cork development in planta. A transcriptome of potato tuber phellogen isolated by laser capture microdissection indicated similarity to vascular cambium and the cork from trees. Highly expressed genes and transcription factors indicated that phellogen activation involves cytokinesis and gene reprograming for the establishment of a dedifferentiation state; whereas inactivation is characterized by activity of genes that direct organ identity in meristem and cell-wall modifications. The expression of selected genes was analyzed using qPCR in native and wound periderm at distinct developmental stages. This allowed the identification of genes involved in periderm formation and maturation.

A comparative transcriptomic approach to understanding the formation of cork

The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.

Cork Oak Young and Traumatic Periderms Show PCD Typical Chromatin Patterns but Different Chromatin-Modifying Genes Expression

Plants are subjected to adverse conditions being outer protective tissues fundamental to their survival. Tree stems are enveloped by a periderm made of cork cells, resulting from the activity of the meristem phellogen. DNA methylation and histone modifications have important roles in the regulation of plant cell differentiation. However, studies on its involvement in cork differentiation are scarce despite periderm importance. Cork oak periderm development was used as a model to study the formation and differentiation of secondary protective tissues, and their behavior after traumatic wounding (traumatic periderm). Nuclei structural changes, dynamics of DNA methylation, and posttranslational histone modifications were assessed in young and traumatic periderms, after cork harvesting. Lenticular phellogen producing atypical non-suberized cells that disaggregate and form pores was also studied, due to high impact for cork industrial uses. Immunolocalization of active and repressive marks, transcription analysis of the corresponding genes, and correlations between gene expression and cork porosity were investigated. During young periderm development, a reduction in nuclei area along with high levels of DNA methylation occurred throughout epidermis disruption. As cork cells became more differentiated, whole nuclei progressive chromatin condensation with accumulation in the nuclear periphery and increasing DNA methylation was observed. Lenticular cells nuclei were highly fragmented with faint 5-mC labeling. Phellogen nuclei were less methylated than in cork cells, and in lenticular phellogen were even lower. No significant differences were detected in H3K4me3 and H3K18ac signals between cork cells layers, although an increase in H3K4me3 signals was found from the phellogen to cork cells. Distinct gene expression patterns in young and traumatic periderms suggest that cork differentiation might be under specific silencing regulatory pathways. Significant correlations were found between QsMET1, QsMET2, and QsSUVH4 gene expression and cork porosity. This work evidences that DNA methylation and histone modifications play a role in cork differentiation and epidermis induced tension-stress. It also provides the first insights into chromatin dynamics during cork and lenticular cells differentiation pointing to a distinct type of remodeling associated with cell death.