Reconstructing the suberin pathway in poplar by chemical and transcriptomic analysis of bark tissues

Expression of dehydroshikimate dehydratase in poplar induces transcriptional and metabolic changes in the phenylpropanoid pathway

Modification of lignin in feedstocks via genetic engineering aims to reduce biomass recalcitrance to facilitate efficient conversion processes. These improvements can be achieved by expressing exogenous enzymes that interfere with native biosynthetic pathways responsible for the production of the lignin precursors. In planta expression of a bacterial 3-dehydroshikimate dehydratase in poplar trees reduced lignin content and altered the monomer composition, which enabled higher yields of sugars after cell wall polysaccharide hydrolysis. Understanding how plants respond to such genetic modifications at the transcriptional and metabolic levels is needed to facilitate further improvement and field deployment. In this work, we acquired fundamental knowledge on lignin-modified poplar expressing 3-dehydroshikimate dehydratase using RNA-seq and metabolomics. The data clearly demonstrate that changes in gene expression and metabolite abundance can occur in a strict spatiotemporal fashion, revealing tissue-specific responses in the xylem, phloem, or periderm. In the poplar line that exhibited the strongest reduction in lignin, we found that 3% of the transcripts had altered expression levels and ~19% of the detected metabolites had differential abundance in the xylem from older stems. The changes affected predominantly the shikimate and phenylpropanoid pathways as well as secondary cell wall metabolism, and resulted in significant accumulation of hydroxybenzoates derived from protocatechuate and salicylate.

MYB24, MYB144, and MYB168 positively regulate suberin biosynthesis at potato tuber wounds during healing

The essence of wound healing is the accumulation of suberin at wounds, which is formed by suberin polyphenolic (SPP) and suberin polyaliphatic (SPA). The biosynthesis of SPP and SPA monomers is catalyzed by several enzyme classes related to phenylpropanoid metabolism and fatty acid metabolism, respectively. However, how suberin biosynthesis is regulated at the transcriptional level during potato (Solanum tuberosum) tuber wound healing remains largely unknown. Here, 6 target genes and 15 transcription factors related to suberin biosynthesis in tuber wound healing were identified by RNA-seq technology and qRT-PCR. Dual luciferase and yeast one-hybrid assays showed that StMYB168 activated the target genes StPAL, StOMT, and St4CL in phenylpropanoid metabolism. Meanwhile, StMYB24 and StMYB144 activated the target genes StLTP, StLACS, and StCYP in fatty acid metabolism, and StFHT involved in the assembly of SPP and SPA domains in both native and wound periderms. More importantly, virus-induced gene silencing in S. tuberosum and transient overexpression in Nicotiana benthamiana assays confirmed that StMYB168 regulates the biosynthesis of free phenolic acids, such as ferulic acid. Furthermore, StMYB24/144 regulated the accumulation of suberin monomers, such as ferulates, α, ω-diacids, and ω-hydroxy acids. In conclusion, StMYB24, StMYB144, and StMYB168 have an elaborate division of labor in regulating the synthesis of suberin during tuber wound healing.

The Analysis, Description, and Examination of the Maize LAC Gene Family's Reaction to Abiotic and Biotic Stress

Laccase (LAC) is a diverse group of genes found throughout the plant genome essential for plant growth and the response to stress by converting monolignin into intricate lignin formations. However, a comprehensive investigation of maize laccase has not yet been documented. A bioinformatics approach was utilized in this research to conduct a thorough examination of maize (Zea mays L.), resulting in the identification and categorization of 22 laccase genes (ZmLAC) into six subfamilies. The gene structure and motifs of each subgroup were largely consistent. The distribution of the 22 LAC genes was uneven among the maize chromosomes, with the exception of chromosome 9. The differentiation of the genes was based on fragment replication, and the differentiation time was about 33.37 million years ago. ZmLAC proteins are primarily acidic proteins. There are 18 cis-acting elements in the promoter sequences of the maize LAC gene family associated with growth and development, stress, hormones, light response, and stress response. The analysis of tissue-specific expression revealed a high expression of the maize LAC gene family prior to the V9 stage, with minimal expression at post-V9. Upon reviewing the RNA-seq information from the publicly available transcriptome, it was discovered that ZmLAC5, ZmLAC10, and ZmLAC17 exhibited significant expression levels when exposed to various biotic and abiotic stress factors, suggesting their crucial involvement in stress responses and potential value for further research. This study offers an understanding of the functions of the LAC genes in maize's response to biotic and abiotic stress, along with a theoretical basis for comprehending the molecular processes at play.

Suberin deficiency and its effect on the transport physiology of young poplar roots

The precise functions of suberized apoplastic barriers in root water and nutrient transport physiology have not fully been elucidated. While lots of research has been performed with mutants of Arabidopsis, little to no data are available for mutants of agricultural crop or tree species. By employing a combined set of physiological, histochemical, analytical, and transport physiological methods as well as RNA-sequencing, this study investigated the implications of remarkable CRISPR/Cas9-induced suberization defects in young roots of the economically important gray poplar. While barely affecting overall plant development, contrary to literature-based expectations significant root suberin reductions of up to 80-95% in four independent mutants were shown to not evidently affect the root hydraulic conductivity during non-stress conditions. In addition, subliminal iron deficiency symptoms and increased translocation of a photosynthesis inhibitor as well as NaCl highlight the involvement of suberin in nutrient transport physiology. The multifaceted nature of the root hydraulic conductivity does not allow drawing simplified conclusions such as that the suberin amount must always be correlated with the water transport properties of roots. However, the decreased masking of plasma membrane surface area could facilitate the uptake but also leakage of beneficial and harmful solutes.

CsMYB67 participates in the flavonoid biosynthesis of summer tea leaves

Flavonoids are important compounds in tea leaves imparting bitter and astringent taste, which also play key roles in tea plants responding to environmental stress. Our previous study showed that the expression level of CsMYB67 was positively correlated with the accumulation of flavonoids in tea leaves as exposed to sunlight. Here, we newly reported the function of CsMYB67 in regulating flavonoid biosynthesis in tea leaves. CsMYB67 was localized in the nucleus and responded to temperature. The results of transient expression assays showed the co-transformation of CsMYB67 and CsTTG1 promoted the transcription of CsANS promoter in the tobacco system. CsTTG1 was bound to the promoter of CsANS based on the results of yeast one-hybrid (Y1H) and transient expression assays, while CsMYB67 enhanced the transcription of CsANS through protein interaction with CsTTG1 according to the results of yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC). Thus, CsMYB67-CsTTG1 module enhanced the anthocyanin biosynthesis through up-regulating the transcription of CsANS. Besides, CsMYB67 also enhanced the transcription of CsFLS and CsUFGT through forming transcription factor complexes. The function of CsMYB67 on flavonoid biosynthesis in tea leaves was validated by gene suppression assay. As CsMYB67 was suppressed, the transcriptional level of CsFLS was greatly reduced, leading to a significant increase in the contents of total catechins and total anthocyanidins. Hence, CsMYB67 plays an important role in regulating the downstream pathway of flavonoid biosynthesis in summer tea leaves.

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.

Comparative Analysis of Cd Uptake and Tolerance in Two Mangrove Species (Avicennia marina and Rhizophora stylosa) with Distinct Apoplast Barriers

Mangrove plants demonstrate an impressive ability to tolerate environmental pollutants, but excessive levels of cadmium (Cd) can impede their growth. Few studies have focused on the effects of apoplast barriers on heavy metal tolerance in mangrove plants. To investigate the uptake and tolerance of Cd in mangrove plants, two distinct mangrove species, Avicennia marina and Rhizophora stylosa, are characterized by unique apoplast barriers. The results showed that both mangrove plants exhibited the highest concentration of Cd2+ in roots, followed by stems and leaves. The Cd2+ concentrations in all organs of R. stylosa consistently exhibited lower levels than those of A. marina. In addition, R. stylosa displayed a reduced concentration of apparent PTS and a smaller percentage of bypass flow when compared to A. marina. The root anatomical characteristics indicated that Cd treatment significantly enhanced endodermal suberization in both A. marina and R. stylosa roots, and R. stylosa exhibited a higher degree of suberization. The transcriptomic analysis of R. stylosa and A. marina roots under Cd stress revealed 23 candidate genes involved in suberin biosynthesis and 8 candidate genes associated with suberin regulation. This study has confirmed that suberized apoplastic barriers play a crucial role in preventing Cd from entering mangrove roots.

Effects of drought and rehydration on root gene expression in seedlings of Pinus massoniana Lamb

The mechanisms underlying plant response to drought involve the expression of numerous functional and regulatory genes. Transcriptome sequencing based on the second- and/or third-generation high-throughput sequencing platforms has proven to be powerful for investigating the transcriptional landscape under drought stress. However, the full-length transcriptomes related to drought responses in the important conifer genus Pinus L. remained to be delineated using the third-generation sequencing technology. With the objectives of identifying the candidate genes responsible for drought and/or rehydration and clarifying the expression profile of key genes involved in drought regulation, we combined the third- and second-generation sequencing techniques to perform transcriptome analysis on seedling roots under drought stress and rewatering in the drought-tolerant conifer Pinus massoniana Lamb. A sum of 294,114 unique full-length transcripts were produced with a mean length of 3217 bp and N50 estimate of 5075 bp, including 279,560 and 124,438 unique full-length transcripts being functionally annotated and Gene Ontology enriched, respectively. A total of 4076, 6295 and 18,093 differentially expressed genes (DEGs) were identified in three pair-wise comparisons of drought-treatment versus control transcriptomes, including 2703, 3576 and 8273 upregulated and 1373, 2719 and 9820 downregulated DEGs, respectively. Moreover, 157, 196 and 691 DEGs were identified as transcription factors in the three transcriptome comparisons and grouped into 26, 34 and 44 transcription factor families, respectively. Gene Ontology enrichment analysis revealed that a remarkable number of DEGs were enriched in soluble sugar-related and cell wall-related processes. A subset of 75, 68 and 97 DEGs were annotated to be associated with starch, sucrose and raffinose metabolism, respectively, while 32 and 70 DEGs were associated with suberin and lignin biosynthesis, respectively. Weighted gene co-expression network analysis revealed modules and hub genes closely related to drought and rehydration. This study provides novel insights into root transcriptomic changes in response to drought dynamics in Masson pine and serves as a fundamental work for further molecular investigation on drought tolerance in conifers.

Endosperm Persistence in Arabidopsis Results in Seed Coat Fractures and Loss of Seed Longevity

Seeds are specialized plant organs that carry, nurture, and protect plant offspring. Developmental coordination between the three genetically distinct seed tissues (the embryo, endosperm, and seed coat) is crucial for seed viability. In this study, we explore the relationship between the TFs AtHB25 and ICE1. Previous results identified ICE1 as a target gene of AtHB25. In seeds, a lack of ICE1 (ice1-2) suppresses the enhanced seed longevity and impermeability of the overexpressing mutant athb25-1D, but surprisingly, seed coat lipid polyester deposition is not affected, as shown by the double-mutant athb25-1D ice1-2 seeds. zou-4, another mutant lacking the transcriptional program for proper endosperm maturation and for which the endosperm persists, also presents a high sensitivity to seed aging. Analysis of gso1, gso2, and tws1-4 mutants revealed that a loss of embryo cuticle integrity does not underlie the seed-aging sensitivity of ice1-2 and zou-4. However, scanning electron microscopy revealed the presence of multiple fractures in the seed coats of the ice1 and zou mutants. Thus, this study highlights the importance of both seed coat composition and integrity in ensuring longevity and demonstrates that these parameters depend on multiple factors.

Cutinized and suberized barriers in leaves and roots: Similarities and differences

Anatomical, histochemical, chemical, and biosynthetic similarities and differences of cutinized and suberized plant cell walls are presented and reviewed in brief. Based on this, the functional properties of cutinized and suberized plant cell walls acting as transport barriers are compared and discussed in more detail. This is of general importance because fundamental misconceptions about relationships in plant-environment water relations are commonly encountered in the scientific literature. It will be shown here, that cuticles represent highly efficient apoplastic transport barriers significantly reducing the diffusion of water and dissolved compounds. The transport barrier of cuticles is mainly established by the deposition of cuticular waxes. Upon wax extraction, with the cutin polymer remaining, cuticular permeability for water and dissolved non-ionized and lipophilic solutes are increasing by 2-3 orders of magnitude, whereas polar and charged substances (e.g., nutrient ions) are only weakly affected (2- to 3-fold increases in permeability). Suberized apoplastic barriers without the deposition of wax are at least as permeable as the cutin polymer matrix without waxes and hardly offer any resistance to the free movement of water. Only upon the deposition of significant amounts of wax, as it is the case with suberized periderms exposed to the atmosphere, an efficient transport barrier for water can be established by suberized cell walls. Comparing the driving forces (gradients between water potentials inside leaves and roots and the surrounding environment) for water loss acting on leaves and roots, it is shown that leaves must have a genetically pre-defined highly efficient transpiration barrier fairly independent from rapidly changing environmental influences. Roots, in most conditions facing a soil environment with relative humidities very close to 100%, are orders of magnitude more permeable to water than leaf cuticles. Upon desiccation, the permanent wilting point of plants is defined as -1.5 MPa, which still corresponds to nearly 99% relative humidity in soil. Thus, the main reason for plant water stress leading to dehydration is the inability of root tissues to decrease their internal water potential to values more negative than -1.5 MPa and not the lack of a transport barrier for water in roots and leaves. Taken together, the commonly mentioned concepts that a drought-induced increase of cuticular wax or root suberin considerably strengthens the apoplastic leaf or root transport barriers and thus aids in water conservation appears highly questionable.

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.

MdMYB52 regulates lignin biosynthesis upon the suberization process in apple

Our previous studies, comparing russeted vs. waxy apple skin, highlighted a MYeloBlastosys (Myb) transcription factor (MdMYB52), which displayed a correlation with genes associated to the suberization process. The present article aims to assess its role and function in the suberization process. Phylogenetic analyses and research against Arabidopsis thaliana MYBs database were first performed and the tissue specific expression of MdMYB52 was investigated using RT-qPCR. The function of MdMYB52 was further investigated using Agrobacterium-mediated transient overexpression in Nicotiana benthamiana leaves. An RNA-Seq analysis was performed to highlight differentially regulated genes in response MdMYB52. Transcriptomic data were supported by analytical chemistry and microscopy. A massive decreased expression of photosynthetic and primary metabolism pathways was observed with a concomitant increased expression of genes associated with phenylpropanoid and lignin biosynthesis, cell wall modification and senescence. Interestingly key genes involved in the synthesis of suberin phenolic components were observed. The analytical chemistry displayed a strong increase in the lignin content in the cell walls during MdMYB52 expression. More specifically, an enrichment in G-Unit lignin residues was observed, supporting transcriptomic data as well as previous work describing the suberin phenolic domain as a G-unit enriched lignin-like polymer. The time-course qPCR analysis revealed that the observed stress response, might be explain by this lignin biosynthesis and by a possible programmed senescence triggered by MdMYB52. The present work supports a crucial regulatory role for MdMYB52 in the biosynthesis of the suberin phenolic domain and possibly in the fate of suberized cells in russeted apple skins.

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.

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.

Chemical and Molecular Characterization of Wound-Induced Suberization in Poplar (Populus alba × P. tremula) Stem Bark

Upon mechanical damage, plants produce wound responses to protect internal tissues from infections and desiccation. Suberin, a heteropolymer found on the inner face of primary cell walls, is deposited in specific tissues under normal development, enhanced under abiotic stress conditions and synthesized by any tissue upon mechanical damage. Wound-healing suberization of tree bark has been investigated at the anatomical level but very little is known about the molecular mechanisms underlying this important stress response. Here, we investigated a time course of wound-induced suberization in poplar bark. Microscopic changes showed that polyphenolics accumulate 3 days post wounding, with aliphatic suberin deposition observed 5 days post wounding. A wound periderm was formed 9 days post wounding. Chemical analyses of the suberin polyester accumulated during the wound-healing response indicated that suberin monomers increased from 0.25 to 7.98 mg/g DW for days 0 to 28, respectively. Monomer proportions varied across the wound-healing process, with an overall ratio of 2:1 (monomers:glycerol) found across the first 14 days post wounding, with this ratio increasing to 7:2 by day 28. The expression of selected candidate genes of poplar suberin metabolism was investigated using qRT-PCR. Genes queried belonging to lipid polyester and phenylpropanoid metabolism appeared to have redundant functions in native and wound-induced suberization. Our data show that, anatomically, the wounding response in poplar bark is similar to that described in periderms of other species. It also provides novel insight into this process at the chemical and molecular levels, which have not been previously studied in trees.

An ARF1-binding factor triggering programmed cell death and periderm development in pear russet fruit skin

Plants have a cuticular membrane (CM) and periderm membrane (PM), which act as barriers to terrestrial stresses. The CM covers primary organs with a continuous hydrophobic layer of waxes embedded in cutin, while the PM stacks with suberized cells outermost to the secondary tissues. The formation of native periderm is regulated by a postembryonic meristem phellogen that produces suberized phellem (cork) outwardly. However, the mechanism controlling phellogen differentiation to phellem remains to be clarified. Here, map-based cloning in a pear F1 population with segregation for periderm development in fruit skin facilitated the identification of an aspartic acid repeat deletion in Pyrus Periderm Programmed Cell Death 1.1 (PyPPCD1.1) that triggers phellogen activity for cork formation in pear russet fruit skin. PyPPCD1.1 showed preferential expression in pear fruit skin, and the encoded protein shares a structural similarity to that of the viral capsid proteins. Asp deletion in PyPPCD1.1 weakened its nuclear localization but increased its accumulation in the chloroplast. Both PyPPCD1.1 and its recessive allele directly interact with ADP-ribosylation factor 1 (ARF1). PyPPCD1.1 triggered PCD in an ARF1-dependent manner. Thus, this study identified the switch gene for PCD and periderm development and provided a new molecular regulatory mechanism underlying the development of this trait.

The metabolic and proteomic repertoires of periderm tissue in skin of the reticulated Sikkim cucumber fruit

Suberized and/or lignified (i.e. lignosuberized) periderm tissue appears often on surface of fleshy fruit skin by mechanical damage caused following environmental cues or developmental programs. The mechanisms underlying lignosuberization remain largely unknown to date. Here, we combined an assortment of microscopical techniques with an integrative multi-omics approach comprising proteomics, metabolomics and lipidomics to identify novel molecular components involved in fruit skin lignosuberization. We chose to investigate the corky Sikkim cucumber (Cucumis sativus var. sikkimensis) fruit. During development, the skin of this unique species undergoes massive cracking and is coated with a thick corky layer, making it an excellent model system for revealing fundamental cellular machineries involved in fruit skin lignosuberization. The large-scale data generated provides a significant source for the field of skin periderm tissue formation in fleshy fruit and suberin metabolism.

Hydroponic cultivation conditions allowing the reproducible investigation of poplar root suberization and water transport

BACKGROUND

With increasing joint research cooperation on national and international levels, there is a high need for harmonized and reproducible cultivation conditions and experimental protocols in order to ensure the best comparability and reliability of acquired data. As a result, not only comparisons of findings of different laboratories working with the same species but also of entirely different species would be facilitated. As Populus is becoming an increasingly important genus in modern science and agroforestry, the integration of findings with previously gained knowledge of other crop species is of high significance.

RESULTS

To ease and ensure the comparability of investigations of root suberization and water transport, on a high degree of methodological reproducibility, we set up a hydroponics-based experimental pipeline. This includes plant cultivation, root histochemistry, analytical investigation, and root water transport measurement. A 5-week-long hydroponic cultivation period including an optional final week of stress application resulted in a highly consistent poplar root development. The poplar roots were of conical geometry and exhibited a typical Casparian band development with subsequent continuously increasing suberization of the endodermis. Poplar root suberin was composed of the most frequently described suberin substance classes, but also high amounts of benzoic acid derivatives could be identified. Root transport physiology experiments revealed that poplar roots in this developmental stage have a two- to tenfold higher hydrostatic than osmotic hydraulic conductivity. Lastly, the hydroponic cultivation allowed the application of gradually defined osmotic stress conditions illustrating the precise adjustability of hydroponic experiments as well as the previously reported sensitivity of poplar plants to water deficits.

CONCLUSIONS

By maintaining a high degree of harmonization, we were able to compare our results to previously published data on root suberization and water transport of barley and other crop species. Regarding hydroponic poplar cultivation, we enabled high reliability, reproducibility, and comparability for future experiments. In contrast to abiotic stress conditions applied during axenic tissue culture cultivation, this experimental pipeline offers great advantages including the growth of roots in the dark, easy access to root systems before, during, and after stress conditions, and the more accurate definition of the developmental stages of the roots.

Effects of Salinity and Abscisic Acid on Lipid Transfer Protein Accumulation, Suberin Deposition and Hydraulic Conductance in Pea Roots

Lipid transfer proteins (LTPs) participate in many important physiological processes in plants, including adaptation to stressors, e.g., salinity. Here we address the mechanism of this protective action of LTPs by studying the interaction between LTPs and abscisic acid (ABA, a "stress" hormone) and their mutual participation in suberin deposition in root endodermis of salt-stressed pea plants. Using immunohistochemistry we show for the first time NaCl induced accumulation of LTPs and ABA in the cell walls of phloem paralleled by suberin deposition in the endoderm region of pea roots. Unlike LTPs which were found localized around phloem cells, ABA was also present within phloem cells. In addition, ABA treatment resulted in both LTP and ABA accumulation in phloem cells and promoted root suberization. These results suggested the importance of NaCl-induced accumulation of ABA in increasing the abundance of LTPs and of suberin. Using molecular modeling and fluorescence spectroscopy we confirmed the ability of different plant LTPs, including pea Ps-LTP1, to bind ABA. We therefore hypothesize an involvement of plant LTPs in ABA transport (unloading from phloem) as part of the salinity adaptation mechanism.

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.

Genetic Architecture and Genome-Wide Adaptive Signatures Underlying Stem Lenticel Traits in Populus tomentosa

The stem lenticel is a highly specialized tissue of woody plants that has evolved to balance stem water retention and gas exchange as an adaptation to local environments. In this study, we applied genome-wide association studies and selective sweeping analysis to characterize the genetic architecture and genome-wide adaptive signatures underlying stem lenticel traits among 303 unrelated accessions of P. tomentosa, which has significant phenotypic and genetic variations according to climate region across its natural distribution. In total, we detected 108 significant single-nucleotide polymorphisms, annotated to 88 candidate genes for lenticel, of which 9 causative genes showed significantly different selection signatures among climate regions. Furthermore, PtoNAC083 and PtoMYB46 showed significant association signals and abiotic stress response, so we overexpressed these two genes in Arabidopsis thaliana and found that the number of stem cells in all three overexpression lines was significantly reduced by PtoNAC083 overexpression but slightly increased by PtoMYB46 overexpression, suggesting that both genes are involved in cell division and expansion during lenticel formation. The findings of this study demonstrate the successful application of an integrated strategy for dissecting the genetic basis and landscape genetics of complex adaptive traits, which will facilitate the molecular design of tree ideotypes that may adapt to future climate and environmental changes.

An integrated metabolic and transcriptomic analysis reveals the mechanism through which fruit bagging alleviates exocarp semi-russeting in pear fruit

Fruit semi-russeting is an undesirable quality trait that occurs in fruit production. It is reported that preharvest fruit bagging could effectively alleviate fruit exocarp semi-russeting, but the physiological and molecular mechanisms remain unclear. In the present study, we performed an in-depth investigation into pear fruit semi-russeting from morphologic, metabolic and transcriptomic perspectives by comparing control (semi-russeted) and bagged (non-russeted) 'Cuiguan' pear fruits. The results showed that significant changes in cutin and suberin resulted in pear fruit semi-russeting. Compared with the skin of bagged fruits, the skin of the control fruits presented reduced cutin contents accompanied by an accumulation of suberin, which resulted in fruit semi-russeting; α, ω-dicarboxylic acids accounted for the largest proportion of typical suberin monomers. Moreover, combined transcriptomic and metabolic analysis revealed a series of genes involved in cutin and suberin biosynthesis, transport and polymerization differentially expressed between the two groups. Furthermore, the expression levels of genes involved in the stress response and in hormone biosynthesis and signaling were significantly altered in fruits with contrasting phenotypes. Finally, a number of transcription factors, including those of the MYB, NAC, bHLH and bZIP families, were differentially expressed. Taken together, the results suggest that the multilayered mechanism through which bagging alleviates pear fruit semi-russeting is complex, and the large number of candidate genes identified provides a good foundation for future functional studies.

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.

Differential regulation of triterpene biosynthesis induced by an early failure in cuticle formation in apple

Waxy apple cuticles predominantly accumulate ursane-type triterpenes, but the profile shifts with the induction of skin russeting towards lupane-type triterpenes. We previously characterised several key enzymes in the ursane-type and lupane-type triterpene pathways, but this switch in triterpene metabolism associated with loss of cuticle integrity is not fully understood. To analyse the relationship between triterpene biosynthesis and russeting, we used microscopy, RNA-sequencing and metabolite profiling during apple fruit development. We compared the skin of three genetically-close clones of 'Golden Delicious' (with waxy, partially russeted and fully russeted skin). We identified a unique molecular profile for the russet clone, including low transcript abundance of multiple cuticle-specific metabolic pathways in the early stages of fruit development. Using correlation analyses between gene transcription and metabolite concentration we found MYB transcription factors strongly associated with lupane-type triterpene biosynthesis. We showed how their transcription changed with the onset of cuticle cracking followed by russeting and that one factor, MYB66, was able to bind the promoter of the oxidosqualene cyclase OSC5, to drive the production of lupeol derivatives. These results provide insights into the breakdown of cuticle integrity leading to russet and how this drives MYB-regulated changes to triterpene biosynthesis.

Importance of suberin biopolymer in plant function, contributions to soil organic carbon and in the production of bio-derived energy and materials

Suberin is a hydrophobic biopolymer of significance in the production of biomass-derived materials and in biogeochemical cycling in terrestrial ecosystems. Here, we describe suberin structure and biosynthesis, and its importance in biological (i.e., plant bark and roots), ecological (soil organic carbon) and economic (biomass conversion to bioproducts) contexts. Furthermore, we highlight the genomics and analytical approaches currently available and explore opportunities for future technologies to study suberin in quantitative and/or high-throughput platforms in bioenergy crops. A greater understanding of suberin structure and production in lignocellulosic biomass can be leveraged to improve representation in life cycle analysis and techno-economic analysis models and enable performance improvements in plant biosystems as well as informed crop system management to achieve economic and environmental co-benefits.

Drought activates MYB41 orthologs and induces suberization of grapevine fine roots

The permeability of roots to water and nutrients is controlled through a variety of mechanisms and one of the most conspicuous is the presence of the Casparian strips and suberin lamellae. Roots actively regulate the creation of these structures developmentally, along the length of the root, and in response to the environment, including drought. In the current study, we characterized the suberin composition along the length of grapevine fine roots during development and in response to water deficit, and in the same root systems we quantified changes in expression of suberin biosynthesis- and deposition-related gene families (via RNAseq) allowing the identification of drought-responsive suberin-related genes. Grapevine suberin composition did not differ between primary and lateral roots, and was similar to that of other species. Under water deficit there was a global upregulation of suberin biosynthesis which resulted in an increase of suberin specific monomers, but without changes in their relative abundances, and this upregulation took place across all the developmental stages of fine roots. These changes corresponded to the upregulation of numerous suberin biosynthesis- and export-related genes which included orthologs of the previously characterized AtMYB41 transcriptional factor. Functional validation of two grapevine MYB41 orthologs, VriMYB41 and VriMYB41-like, confirmed their ability to globally upregulate suberin biosynthesis, export, and deposition. This study provides a detailed characterization of the developmental and water deficit induced suberization of grapevine fine roots and identifies important orthologs responsible for suberin biosynthesis, export, and its regulation in grape.

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.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

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.

PRX2 and PRX25, peroxidases regulated by COG1, are involved in seed longevity in Arabidopsis

Permeability is a crucial trait that affects seed longevity and is regulated by different polymers including proanthocyanidins, suberin, cutin and lignin located in the seed coat. By testing mutants in suberin transport and biosynthesis, we demonstrate the importance of this biopolymer to cope with seed deterioration. Transcriptomic analysis of cog1-2D, a gain-of-function mutant with increased seed longevity, revealed the upregulation of several peroxidase genes. Reverse genetics analysing seed longevity uncovered redundancy within the seed coat peroxidase gene family; however, after controlled deterioration treatment, seeds from the prx2 prx25 double and prx2 prx25 prx71 triple mutant plants presented lower germination than wild-type plants. Transmission electron microscopy analysis of the seed coat of these mutants showed a thinner palisade layer, but no changes were observed in proanthocyanidin accumulation or in the cuticle layer. Spectrophotometric quantification of acetyl bromide-soluble lignin components indicated changes in the amount of total polyphenolics derived from suberin and/or lignin in the mutant seeds. Finally, the increased seed coat permeability to tetrazolium salts observed in the prx2 prx25 and prx2 prx25 prx71 mutant lines suggested that the lower permeability of the seed coats caused by altered polyphenolics is likely to be the main reason explaining their reduced seed longevity.

Investigation of Salt Tolerance Mechanisms Across a Root Developmental Gradient in Almond Rootstocks

The intensive use of groundwater in agriculture under the current climate conditions leads to acceleration of soil salinization. Given that almond is a salt-sensitive crop, selection of salt-tolerant rootstocks can help maintain productivity under salinity stress. Selection for tolerant rootstocks at an early growth stage can reduce the investment of time and resources. However, salinity-sensitive markers and salinity tolerance mechanisms of almond species to assist this selection process are largely unknown. We established a microscopy-based approach to investigate mechanisms of stress tolerance in and identified cellular, root anatomical, and molecular traits associated with rootstocks exhibiting salt tolerance. We characterized three almond rootstocks: Empyrean-1 (E1), Controller-5 (C5), and Krymsk-86 (K86). Based on cellular and molecular evidence, our results show that E1 has a higher capacity for salt exclusion by a combination of upregulating ion transporter expression and enhanced deposition of suberin and lignin in the root apoplastic barriers, exodermis, and endodermis, in response to salt stress. Expression analyses revealed differential regulation of cation transporters, stress signaling, and biopolymer synthesis genes in the different rootstocks. This foundational study reveals the mechanisms of salinity tolerance in almond rootstocks from cellular and structural perspectives across a root developmental gradient and provides insights for future screens targeting stress response.

Analysis of Orthologous SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) Promotor Activity in Herbaceous and Woody Angiosperms

SECONDARY WALL-ASSOCIATED NAC DOMAIN1 (SND1) is a master regulator of fibre secondary wall deposition in Arabidopsis thaliana (Arabidopsis), with homologs in other angiosperms and gymnosperms. However, it is poorly understood to what extent the fibre-specific regulation of the SND1 promoter, and that of its orthologs, is conserved between diverged herbaceous and woody lineages. We performed a reciprocal reporter gene analysis of orthologous SND1 promoters from Arabidopsis (AthSND1), Eucalyptus grandis (EgrNAC61) and Populus alba × P. grandidentata (PagWND1A) relative to secondary cell wall-specific Cellulose Synthase4 (CesA4) and CesA7 promoters, in both a non-woody (Arabidopsis) and a woody (poplar) system. β-glucuronidase (GUS) reporter analysis in Arabidopsis showed that the SND1 promoter was active in vascular tissues as previously reported and showed interfascicular and xylary fibre-specific expression in inflorescence stems, while reporter constructs of the woody plant-derived promoters were partial to the (pro)cambium-phloem and protoxylem. In transgenic P. tremula × P. alba plants, all three orthologous SND1 promoters expressed the GUS reporter similarly and preferentially in developing secondary xylem, ray parenchyma and cork cambium. Ours is the first study to reciprocally test orthologous SND1 promoter specificity in herbaceous and woody species, revealing diverged regulatory functions in the herbaceous system.

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.

ChIP-Seq reveals that QsMYB1 directly targets genes involved in lignin and suberin biosynthesis pathways in cork oak (Quercus suber)

BACKGROUND

Gene activity is largely controlled by transcriptional regulation through the action of transcription factors and other regulators. QsMYB1 is a member of the R2R3-MYB transcription factor family related to secondary growth, and in particular, with the cork development process. In order to identify the putative gene targets of QsMYB1 across the cork oak genome we developed a ChIP-Seq strategy.

RESULTS

Results provide direct evidence that QsMY1B targets genes encoding for enzymes involved in the lignin and suberin pathways as well as gene encoding for ABCG transporters and LTPs implicated in the transport of monomeric suberin units across the cellular membrane. These results highlight the role of QsMYB1 as a regulator of lignin and suberin biosynthesis, transport and assembly.

CONCLUSION

To our knowledge, this work constitutes the first ChIP-Seq experiment performed in cork oak, a non-model plant species with a long-life cycle, and these results will contribute to deepen the knowledge about the molecular mechanisms of cork formation and differentiation.

Transcriptional profiling of cork oak phellogenic cells isolated by laser microdissection

The phenylpropanoid pathway impacts the cork quality development. In cork of bad quality, the flavonoid route is favored, whereas in good quality, cork lignin and suberin production prevails. Cork oaks develop a thick cork tissue as a protective shield that results of the continuous activity of a secondary meristem, the cork cambium, or phellogen. Most studies applied to developmental processes do not consider the cell types from which the samples were extracted. Here, laser microdissection (LM) coupled with transcript profiling using RNA sequencing (454 pyrosequencing) was applied to phellogen cells of trees producing low- and good quality cork. Functional annotation and functional enrichment analyses showed that stress-related genes are enriched in samples extracted from trees producing good quality cork (GQC). This process is under tight transcriptional (transcription factors, kinases) regulation and also hormonal control involving ABA, ethylene, and auxins. The phellogen cells collected from trees producing bad quality cork (BQC) show a consistent up-regulation of genes belonging to the flavonoid pathway as a response to stress. They also display a different modulation of cell wall genes resulting into a thinner cork layer, i.e., less meristematic activity. Based on the analysis of the phenylpropanoid pathway regulating genes, in GQC, the synthesis of lignin and suberin is promoted, whereas in BQC, the same pathway favors the biosynthesis of free phenolic compounds. This study provided new insights of how cell-specific gene expression can determine tissue and organ morphology and physiology and identified robust candidate genes that can be used in breeding programs aiming at improving cork quality.

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.