PagPXYs improve drought tolerance by regulating reactive oxygen species homeostasis in the cambium of Populus alba × P. glandulosa.

PXY (Phloem intercalated with xylem) is a receptor kinase required for directional cell division during the development of plant vascular tissue. Drought stress usually affects plant stem cell division and differentiation thereby limiting plant growth. However, the role of PXY in cambial activities of woody plants under drought stress is unclear. In this study, we analyzed the biological functions of two PXY genes (PagPXYa and PagPXYb) in poplar growth and development and in response to drought stress in a hybrid poplar (Populus alba × P. glandulosa, '84K'). Expression analysis indicated that PagPXYs, similar to their orthologs PtrPXYs in Populus trichocarpa, are mainly expressed in the stem vascular system, and related to drought. Interestingly, overexpression of PagPXYa and PagPXYb in poplar did not have a significant impact on the growth status of transgenic plants under normal condition. However, when treated with 8 % PEG6000 or 100 mM H2O2, PagPXYa and PagPXYb overexpressing lines consistently exhibited more cambium cell layers, fewer xylem cell layers, and enhanced drought tolerance compared to the non-transgenic control '84K'. In addition, PagPXYs can alleviate the damage caused by H2O2 to the cambium under drought stress, thereby maintaining the cambial division activity of poplar under drought stress, indicating that PagPXYs play an important role in plant resistance to drought stress. This study provides a new insight for further research on the balance of growth and drought tolerance in forest trees.

Studying plant vascular development using single-cell approaches.

Vascular cells form a highly complex and heterogeneous tissue. Its composition, function, shape, and arrangement vary with the developmental stage and between organs and species. Understanding the transcriptional regulation underpinning this complexity thus requires a high-resolution technique that is capable of capturing rapid events during vascular cell formation. Single-cell and single-nucleus RNA sequencing (sc/snRNA-seq) approaches provide powerful tools to extract transcriptional information from these lowly abundant and dynamically changing cell types, which allows the reconstruction of developmental trajectories. Here, we summarize and reflect on recent studies using single-cell transcriptomics to study vascular cell types and discuss current and future implementations of sc/snRNA-seq approaches in the field of vascular development.

Transcription factor PagMYB31 positively regulates cambium activity and negatively regulates xylem development in poplar.

Wood formation involves consecutive developmental steps, including cell division of vascular cambium, xylem cell expansion, secondary cell wall (SCW) deposition, and programmed cell death. In this study, we identified PagMYB31 as a coordinator regulating these processes in Populus alba × Populus glandulosa and built a PagMYB31-mediated transcriptional regulatory network. PagMYB31 mutation caused fewer layers of cambial cells, larger fusiform initials, ray initials, vessels, fiber and ray cells, and enhanced xylem cell SCW thickening, showing that PagMYB31 positively regulates cambial cell proliferation and negatively regulates xylem cell expansion and SCW biosynthesis. PagMYB31 repressed xylem cell expansion and SCW thickening through directly inhibiting wall-modifying enzyme genes and the transcription factor genes that activate the whole SCW biosynthetic program, respectively. In cambium, PagMYB31 could promote cambial activity through TRACHEARY ELEMENT DIFFERENTIATION INHIBITORY FACTOR (TDIF)/PHLOEM INTERCALATED WITH XYLEM (PXY) signaling by directly regulating CLAVATA3/ESR-RELATED (CLE) genes, and it could also directly activate WUSCHEL HOMEOBOX RELATED4 (PagWOX4), forming a feedforward regulation. We also observed that PagMYB31 could either promote cell proliferation through the MYB31-MYB72-WOX4 module or inhibit cambial activity through the MYB31-MYB72-VASCULAR CAMBIUM-RELATED MADS2 (VCM2)/PIN-FORMED5 (PIN5) modules, suggesting its role in maintaining the homeostasis of vascular cambium. PagMYB31 could be a potential target to manipulate different developmental stages of wood formation.

Network Analysis of Metabolome and Transcriptome Revealed Regulation of Different Nitrogen Concentrations on Hybrid Poplar Cambium Development.

Secondary development is a key biological characteristic of woody plants and the basis of wood formation. Exogenous nitrogen can affect the secondary growth of poplar, and some regulatory mechanisms have been found in the secondary xylem. However, the effect of nitrogen on cambium has not been reported. Herein, we investigated the effects of different nitrogen concentrations on cambium development using combined transcriptome and metabolome analysis. The results show that, compared with 1 mM NH4NO3 (M), the layers of hybrid poplar cambium cells decreased under the 0.15 mM NH4NO3 (L) and 0.3 mM NH4NO3 (LM) treatments. However, there was no difference in the layers of hybrid poplar cambium cells under the 3 mM NH4NO3 (HM) and 5 mM NH4NO3 (H) treatments. Totals of 2365, 824, 649 and 398 DEGs were identified in the M versus (vs.) L, M vs. LM, M vs. HM and M vs. H groups, respectively. Expression profile analysis of the DEGs showed that exogenous nitrogen affected the gene expression involved in plant hormone signal transduction, phenylpropanoid biosynthesis, the starch and sucrose metabolism pathway and the ubiquitin-mediated proteolysis pathway. In M vs. L, M vs. LM, M vs. HM and M vs. H, differential metabolites were enriched in flavonoids, lignans, coumarins and saccharides. The combined analysis of the transcriptome and metabolome showed that some genes and metabolites in plant hormone signal transduction, phenylpropanoid biosynthesis and starch and sucrose metabolism pathways may be involved in nitrogen regulation in cambium development, whose functions need to be verified. In this study, from the point of view that nitrogen influences cambium development to regulate wood formation, the network analysis of the transcriptome and metabolomics of cambium under different nitrogen supply levels was studied for the first time, revealing the potential regulatory and metabolic mechanisms involved in this process and providing new insights into the effects of nitrogen on wood development.

Gibberellin promotes cambium reestablishment during secondary vascular tissue regeneration after girdling in an auxin-dependent manner in Populus.

Secondary vascular tissue (SVT) development and regeneration are regulated by phytohormones. In this study, we used an in vitro SVT regeneration system to demonstrate that gibberellin (GA) treatment significantly promotes auxin-induced cambium reestablishment. Altering GA content by overexpressing or knocking down ent-kaurene synthase (KS) affected secondary growth and SVT regeneration in poplar. The poplar DELLA gene GIBBERELLIC ACID INSENSITIVE (PtoGAI) is expressed in a specific pattern during secondary growth and cambium regeneration after girdling. Overexpression of PtoGAI disrupted poplar growth and inhibited cambium regeneration, and the inhibition of cambium regeneration could be partially restored by GA application. Further analysis of the PtaDR5:GUS transgenic plants, the localization of PIN-FORMED 1 (PIN1) and the expression of auxin-related genes found that an additional GA treatment could enhance the auxin response as well as the expression of PIN1, which mediates auxin transport during SVT regeneration. Taken together, these findings suggest that GA promotes cambium regeneration by stimulating auxin signal transduction.

Genes expression profiles in vascular cambium of Eucalyptus urophylla × Eucalyptus grandis at different ages.

BACKGROUND: Wood is a secondary xylem generated by vascular cambium. Vascular cambium activities mainly include cambium proliferation and vascular tissue formation through secondary growth, thereby producing new secondary phloem inward and secondary xylem outward and leading to continuous tree thickening and wood formation. Wood formation is a complex biological process, which is strictly regulated by multiple genes. Therefore, molecular level research on the vascular cambium of different tree ages can lead to the identification of both key and related genes involved in wood formation and further explain the molecular regulation mechanism of wood formation. RESULTS: In the present study, RNA-Seq and Pac-Bio Iso-Seq were used for profiling gene expression changes in Eucalyptus urophylla × Eucalyptus grandis (E. urograndis) vascular cambium at four different ages. A total of 59,770 non-redundant transcripts and 1892 differentially expressed genes (DEGs) were identified. The expression trends of the DEGs related to cell division and differentiation, cell wall biosynthesis, phytohormone, and transcription factors were analyzed. The DEGs encoding expansin, kinesin, cycline, PAL, GRP9, KNOX, C2C2-dof, REV, etc., were highly expressed in E. urograndis at three years old, leading to positive effects on growth and development. Moreover, some gene family members, such as NAC, MYB, HD-ZIP III, RPK, and RAP, play different regulatory roles in wood formation because of their sophisticated transcriptional network and function redundantly. CONCLUSIONS: These candidate genes are a potential resource to further study wood formation, especially in fast-growing and adaptable eucalyptus. The results may also serve as a basis for further research to unravel the molecular mechanism underlying wood formation.

Combining single-cell RNA sequencing with spatial transcriptome analysis reveals dynamic molecular maps of cambium differentiation in the primary and secondary growth of trees.

Primary and secondary growth of the tree stem are responsible for corresponding increases in trunk height and diameter. However, our molecular understanding of the biological processes that underlie these two types of growth is incomplete. In this study, we used single-cell RNA sequencing and spatial transcriptome sequencing to reveal the transcriptional landscapes of primary and secondary growth tissues in the Populus stem. Comparison between the cell atlas and differentiation trajectory of primary and secondary growth revealed different regulatory networks involved in cell differentiation from cambium to xylem precursors and phloem precursors. These regulatory networks may be controlled by auxin accumulation and distribution. Analysis of cell differentiation trajectories suggested that vessel and fiber development followed a sequential pattern of progressive transcriptional regulation. This research provides new insights into the processes of cell identity and differentiation that occur throughout primary and secondary growth of tree stems, increasing our understanding of the cellular differentiation dynamics that occur during stem growth in trees.

The LBD11-ROS feedback regulatory loop modulates vascular cambium proliferation and secondary growth in Arabidopsis.

Vascular cambium produces the phloem and xylem, vascular tissues that transport resources and provide mechanical support, making it an ideal target for crop improvement. However, much remains unknown about how vascular cambium proliferates. In this study, through pharmaceutical and genetic manipulation of reactive oxygen species (ROS) maxima, we demonstrate a direct link between levels of ROS and activity of LATERAL ORGAN BOUNDARIES DOMAIN 11 (LBD11) in maintaining vascular cambium activity. LBD11 activates the transcription of several key ROS metabolic genes, including PEROXIDASE 71 and RESPIRATORY BURST OXIDASE HOMOLOGS D and F, to generate local ROS maxima in cambium, which in turn enhance the proliferation of cambial cells. In a negative feedback mechanism, higher ROS levels then repress LBD11 expression and maintain the balance of cambial cell proliferation. Our findings thus reveal the role of a novel LBD11/ROS-dependent feedback regulatory system in maintaining vascular cambium-specific redox homeostasis and radial growth in plants.

Sucrose Signaling Contributes to the Maintenance of Vascular Cambium by Inhibiting Cell Differentiation.

Plants produce sugars by photosynthesis and use them for growth and development. Sugars are transported from source-to-sink organs via the phloem in the vasculature. It is well known that vascular development is precisely controlled by plant hormones and peptide hormones. However, the role of sugars in the regulation of vascular development is poorly understood. In this study, we examined the effects of sugars on vascular cell differentiation using a vascular cell induction system named 'Vascular Cell Induction Culture System Using Arabidopsis Leaves' (VISUAL). We found that sucrose has the strongest inhibitory effect on xylem differentiation, among several types of sugars. Transcriptome analysis revealed that sucrose suppresses xylem and phloem differentiation in cambial cells. Physiological and genetic analyses suggested that sucrose might function through the BRI1-EMS-SUPPRESSOR1 transcription factor, which is the central regulator of vascular cell differentiation. Conditional overexpression of cytosolic invertase led to a decrease in the number of cambium layers due to an imbalance between cell division and differentiation. Taken together, our results suggest that sucrose potentially acts as a signal that integrates environmental conditions with the developmental program.

High-resolution anatomical and spatial transcriptome analyses reveal two types of meristematic cell pools within the secondary vascular tissue of poplar stem.

The secondary vascular tissue emanating from meristems is central to understanding how vascular plants such as forest trees evolve, grow, and regulate secondary radial growth. However, the overall molecular characterization of meristem origins and developmental trajectories from primary to secondary vascular tissues in woody tree stems is technically challenging. In this study, we combined high-resolution anatomic analysis with a spatial transcriptome (ST) technique to define features of meristematic cells in a developmental gradient from primary to secondary vascular tissues in poplar stems. The tissue-specific gene expression of meristems and derived vascular tissue types were accordingly mapped to specific anatomical domains. Pseudotime analyses were used to track the origins and changes of meristems throughout the development from primary to secondary vascular tissues. Surprisingly, two types of meristematic-like cell pools within secondary vascular tissues were inferred based on high-resolution microscopy combined with ST, and the results were confirmed by in situ hybridization of, transgenic trees, and single-cell sequencing. The rectangle shape procambium-like (PCL) cells develop from procambium meristematic cells and are located within the phloem domain to produce phloem cells, whereas fusiform shape cambium zone (CZ) meristematic cells develop from fusiform metacambium meristematic cells and are located inside the CZ to produce xylem cells. The gene expression atlas and transcriptional networks spanning the primary transition to secondary vascular tissues generated in this work provide new resources for studying the regulation of meristem activities and the evolution of vascular plants. A web server (https://pgx.zju.edu.cn/stRNAPal/) was also established to facilitate the use of ST RNA-seq data.

Cell-type-specific PtrWOX4a and PtrVCS2 form a regulatory nexus with a histone modification system for stem cambium development in Populus trichocarpa.

Stem vascular cambium cells in forest trees produce wood for materials and energy. WOX4 affects the proliferation of such cells in Populus. Here we show that PtrWOX4a is the most highly expressed stem vascular-cambium-specific (VCS) gene in P. trichocarpa, and its expression is controlled by the product of the second most highly expressed VCS gene, PtrVCS2, encoding a zinc finger protein. PtrVCS2 binds to the PtrWOX4a promoter as part of a PtrWOX13a-PtrVCS2-PtrGCN5-1-PtrADA2b-3 protein tetramer. PtrVCS2 prevented the interaction between PtrGCN5-1 and PtrADA2b-3, resulting in H3K9, H3K14 and H3K27 hypoacetylation at the PtrWOX4a promoter, which led to fewer cambium cell layers. These effects on cambium cell proliferation were consistent across more than 20 sets of transgenic lines overexpressing individual genes, gene-edited mutants and RNA interference lines in P. trichocarpa. We propose that the tetramer-PtrWOX4a system may coordinate genetic and epigenetic regulation to maintain normal vascular cambium development for wood formation.

A Constitutively Active Cytokinin Receptor Variant Increases Cambial Activity and Stem Growth in Poplar.

The cambial meristem is responsible for bark and wood formation in woody plants. The activity of the cambial meristem is controlled by various factors; one of them is the plant hormone cytokinin. Here, we have explored different approaches to genetically engineering cambial activity in poplar plants by the ectopic expression of a cytokinin biosynthesis gene with enhanced activity (named ROCK4) or of a gene encoding a constitutively active cytokinin receptor variant (ROCK3). Both genes are derived from Arabidopsis thaliana and were expressed in poplar trees under the control of their own promoter or the cambium-specific pHB8 promoter. pIPT3:ROCK4- and pHB8:ROCK4-expressing plants were smaller than wild-type plants and formed more lateral branches; pHB8:ROCK4 transgenic plants additionally showed an increased stem diameter. In contrast, pAHK3:ROCK3- and pHB8:ROCK3-expressing plants grew taller than wild type without an altered branching pattern and formed more cambial cells, leading to increased radial stem growth. The effectivity of ROCK3 when expressed in either secondary phloem cells or in cambial cells is consistent with a dual, tissue-autonomous and non-autonomous activity of cytokinin in regulating cambial activity. We propose ROCK3 as a novel gene to enhance biomass formation in woody plants.

Control of cambium initiation and activity in Arabidopsis by the transcriptional regulator AHL15.

Plant secondary growth, which is the basis of wood formation, includes the production of secondary xylem, which is derived from meristematic cambium cells embedded in vascular tissue. Here, we identified an important role for the Arabidopsis thaliana (Arabidopsis) AT-HOOK MOTIF CONTAINING NUCLEAR LOCALIZED 15 (AHL15) transcriptional regulator in controlling vascular cambium activity. The limited secondary xylem development in inflorescence stems of herbaceous Arabidopsis plants was significantly reduced in ahl15 loss-of-function mutants, whereas constitutive or vascular meristem-specific AHL15 overexpression produced woody inflorescence stems. AHL15 was required for enhanced secondary xylem formation in the woody suppressor of overexpression of constans 1 (soc1) fruitfull (ful) double loss-of-function mutant. Moreover, we found that AHL15 induces vascular cambium activity downstream of the repressing SOC1 and FUL transcription factors, most likely similar to how it enhances lateral branching by promoting biosynthesis of the hormone cytokinin. Our results uncover a novel pathway driving cambium development, in which AHL15 expression levels act in parallel to and are dependent on the well-established TDIF-PXY-WOX pathway to differentiate between herbaceous and woody stem growth.

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.

Non-Coding RNA Analyses of Seasonal Cambium Activity in Populus tomentosa.

Non-coding RNA, known as long non-coding RNA (lncRNA), circular RNA (circRNA) and microRNA (miRNA), are taking part in the multiple developmental processes in plants. However, the roles of which played during the cambium activity periodicity of woody plants remain poorly understood. Here, lncRNA/circRNA-miRNA-mRNA regulatory networks of the cambium activity periodicity in Populus tomentosa was constructed, combined with morphologic observation and transcriptome profiling. Light microscopy and Periodic Acid Schiff (PAS) staining revealed that cell walls were much thicker and number of cell layers was increased during the active-dormant stage, accompanied by abundant change of polysaccharides. The novel lncRNAs and circRNAs were investigated, and we found that 2037 lncRNAs and 299 circRNAs were differentially expression during the vascular cambium period, respectively. Moreover, 1046 genes were identified as a target gene of 2037 novel lncRNAs, and 89 of which were the miRNA precursors or targets. By aligning miRNA precursors to the 7655 lncRNAs, 21 lncRNAs were identified as precursors tof 19 known miRNAs. Furthermore, the target mRNA of lncRNA/circRNA-miRNA network mainly participated in phytohormone, cell wall alteration and chlorophyll metabolism were analyzed by GO enrichment and KEGG pathway. Especially, circRNA33 and circRNA190 taking part in the phytohormone signal pathway were down-regulated during the active-dormant transition. Xyloglucan endotransglucosylase/hydrolase protein 24-like and UDP-glycosyltransferase 85A1 involved in the cell wall modification were the targets of lncRNA MSTRG.11198.1 and MSTRG.1050.1. Notably, circRNA103 and MSTRG.10851.1 regulate the cambium periodicity may interact with the miR482. These results give a new light into activity-dormancy regulation, associated with transcriptional dynamics and non-coding RNA networks of potential targets identification.

The role of senescence-associated gene101 (PagSAG101a) in the regulation of secondary xylem formation in poplar.

Wood is produced by the accumulation of secondary xylem via proliferation and differentiation of the cambium cells in woody plants. Identifying the regulators involved in this process remains a challenging task. In this study, we isolated PagSAG101a, the homolog of Arabidopsis thaliana SAG101, from a hybrid poplar (Populus alba × Populus glandulosa) clone 84K and investigated its role in secondary xylem development. PagSAG101a was expressed predominantly in lignified stems and localized in the nucleus. Compared with non-transgenic 84K plants, transgenic plants overexpressing PagSAG101a displayed increased plant height, internode number, stem diameter, xylem width, and secondary cell wall thickness, while opposite phenotypes were observed for PagSAG101a knock-out plants. Transcriptome analyses revealed that differentially expressed genes were enriched for those controlling cambium cell division activity and subsequent secondary cell wall deposition during xylem formation. In addition, the tandem CCCH zinc finger protein PagC3H17, which positively regulates secondary xylem width and secondary wall thickening in poplar, could bind to the promoter of PagSAG101a and mediate the regulation of xylem differentiation. Our results support that PagSAG101a, downstream of PagC3H17, functions in wood development.

Strigolactone signaling regulates cambial activity through repression of WOX4 by transcription factor BES1.

During secondary growth, meristematic cells in the cambium can either proliferate to maintain the stem cell population or differentiate into xylem or phloem. The balance between these two developmental trajectories is tightly regulated by many environmental and endogenous cues. Strigolactones (SLs), a class of plant hormones, were previously reported to regulate secondary growth by promoting cambium activity. However, the underlying molecular mechanisms of SL action in plant secondary growth are not well understood. We performed histological, genetic, and biochemical analyses using genetic materials in Arabidopsis (Arabidopsis thaliana) with altered activity of the transcription factors BRI1-EMS-SUPPRESSOR1 (BES1) or WUSCHEL-related HOMEOBOX4 (WOX4) or lacking MORE AXILLARY SHOOT2 (MAX2), a key positive component in the SL signaling pathway. We found that BES1, a downstream regulator in the SL signaling pathway that promotes shoot branching and xylem differentiation, also inhibits WOX4 expression, a key regulator of cambium cell division in the intercellular TRACHEARY ELEMENT DIFFERENTIATION INHIBITORY FACTOR (TDIF)-TDIF RECEPTOR (TDR) signaling pathway. The antagonistic roles of BES1 and WOX4 in the regulation of cambium activity may integrate intercellular TDIF signals to efficiently and bidirectionally modulate cambium cell proliferation and differentiation. As both BES1 and WOX4 are widely involved in various endogenous signals and responses to environmental stimuli, these findings may provide insight into the dynamic regulation of cambium development.

Transcriptional landscape of highly lignified poplar stems at single-cell resolution.

BACKGROUND: Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood. RESULTS: Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees. CONCLUSIONS: These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees.

The transcription factor PagLBD3 contributes to the regulation of secondary growth in Populus.

LATERAL ORGAN BOUNDARIES DOMAIN (LBD) genes encode plant-specific transcription factors that participate in regulating various developmental processes. In this study, we genetically characterized PagLBD3 encoding an important regulator of secondary growth in poplar (Populus alba × Populus glandulosa). Overexpression of PagLBD3 increased stem secondary growth in Populus with a significantly higher rate of cambial cell differentiation into phloem, while dominant repression of PagLBD3 significantly decreased the rate of cambial cell differentiation into phloem. Furthermore, we identified 1756 PagLBD3 genome-wide putative direct target genes (DTGs) through RNA sequencing (RNA-seq)-coupled DNA affinity purification followed by sequencing (DAP-seq) assays. Gene Ontology analysis revealed that genes regulated by PagLBD3 were enriched in biological pathways regulating meristem development, xylem development, and auxin transport. Several central regulator genes for vascular development, including PHLOEM INTERCALATED WITH XYLEM (PXY), WUSCHEL RELATED HOMEOBOX4 (WOX4), Secondary Wall-Associated NAC Domain 1s (SND1-B2), and Vascular-Related NAC-Domain 6s (VND6-B1), were identified as PagLBD3 DTGs. Together, our results indicate that PagLBD3 and its DTGs form a complex transcriptional network to modulate cambium activity and phloem/xylem differentiation.

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.