A comparative transcriptomic approach to understanding the formation of cork

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.

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.

Beyond width and density: stable carbon and oxygen isotopes in cork-rings provide insights of physiological responses to water stress in Quercus suber L

As climate change increasingly affects forest ecosystems, detailed understanding of major effects is important to anticipate their consequences under future climate scenarios. The Mediterranean region is a prominent climate change hotspot, and evergreen cork oak (Quercus suber L.) woodlands are particularly climatically sensitive due to cork (bark) harvesting. Cork oak's drought avoidance strategy is well-known and includes structural and physiological adaptations that maximise soil water uptake and transport and limit water use, potentially leading to reduced stem and cork growth. Trees' responses to cope with water-limited conditions have been extensively described based on cork-rings width and, more recently, on cork-rings density, in dendroecological studies. However, so far, tree functional attributes and physiological strategies, namely photosynthetic metabolism adjustments affecting cork formation, have never been addressed and/or integrated on these previous cork-rings-based studies. In this study, we address the relation between carbon and oxygen stable isotopes of cork rings and precipitation and temperature, in two distinct locations of southwestern Portugal-the (wetter) Tagus basin peneplain and the (drier) Grândola mountains. We aimed at assessing whether the two climatic factors affect cork-ring isotopic composition under contrasting conditions of water availability, and, therefore, if carbon and oxygen signatures in cork can reflect tree functional (physiological and structural) responses to stressful conditions, which might be aggravated by climate change. Our results indicate differences between the study areas. At the drier site, the stronger statistically significant negative cork δ ¹³C correlations were found with mean temperature, whereas strong positive cork δ ¹⁸O correlations were fewer and found only with precipitation. Moreover, at the wetter site, cork rings are enriched in ¹⁸O and depleted in ¹³C, indicating, respectively, shallow groundwater as the water source for physiological processes related with biosynthesis of non-photosynthetic secondary tissues, such as suberin, and a weak stomatal regulation under high water availability, consistent with non-existent water availability constrains. In contrast, at the drier site, trees use water from deeper ground layers, depleted in ¹⁸O, and strongly regulate stomatal conductance under water stress, thus reducing photosynthetic carbon uptake and probably relying on stored carbon reserves for cork ring formation. These results suggest that although stable isotopes signatures in cork rings are not proxies for net growth, they may be (fairly) robust indicators of trees' physiological and structural adjustments to climate and environmental changes in Mediterranean environments.

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.

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.

Peridermal fruit skin formation in Actinidia sp. (kiwifruit) is associated with genetic loci controlling russeting and cuticle formation

BACKGROUND

The skin (exocarp) of fleshy fruit is hugely diverse across species. Most fruit types have a live epidermal skin covered by a layer of cuticle made up of cutin while a few create an outermost layer of dead cells (peridermal layer).

RESULTS

In this study we undertook crosses between epidermal and peridermal skinned kiwifruit, and showed that epidermal skin is a semi-dominant trait. Furthermore, backcrossing these epidermal skinned hybrids to a peridermal skinned fruit created a diverse range of phenotypes ranging from epidermal skinned fruit, through fruit with varying degrees of patches of periderm (russeting), to fruit with a complete periderm. Quantitative trait locus (QTL) analysis of this population suggested that periderm formation was associated with four loci. These QTLs were aligned either to ones associated with russet formation on chromosome 19 and 24, or cuticle integrity and coverage located on chromosomes 3, 11 and 24.

CONCLUSION

From the segregation of skin type and QTL analysis, it appears that skin development in kiwifruit is controlled by two competing factors, cuticle strength and propensity to russet. A strong cuticle will inhibit russeting while a strong propensity to russet can create a continuous dead skinned 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.

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

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.

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.

Tissue-specific study across the stem reveals the chemistry and transcriptome dynamics of birch bark

Tree bark is a highly specialized array of tissues that plays important roles in plant protection and development. Bark tissues develop from two lateral meristems; the phellogen (cork cambium) produces the outermost stem-environment barrier called the periderm, while the vascular cambium contributes with phloem tissues. Although bark is diverse in terms of tissues, functions and species, it remains understudied at higher resolution. We dissected the stem of silver birch (Betula pendula) into eight major tissue types, and characterized these by a combined transcriptomics and metabolomics approach. We further analyzed the varying bark types within the Betulaceae family. The two meristems had a distinct contribution to the stem transcriptomic landscape. Furthermore, inter- and intraspecies analyses illustrated the unique molecular profile of the phellem. We identified multiple tissue-specific metabolic pathways, such as the mevalonate/betulin biosynthesis pathway, that displayed differential evolution within the Betulaceae. A detailed analysis of suberin and betulin biosynthesis pathways identified a set of underlying regulators and highlighted the important role of local, small-scale gene duplication events in the evolution of metabolic pathways. This work reveals the transcriptome and metabolic diversity among bark tissues and provides insights to its development and evolution, as well as its biotechnological applications.

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.