Molecular control of wood formation in trees

PagJAZ5 regulates cambium activity through coordinately modulating cytokinin concentration and signaling in poplar

Wood is resulted from the radial growth paced by the division and differentiation of vascular cambium cells in woody plants, and phytohormones play important roles in cambium activity. Here, we identified that PagJAZ5, a key negative regulator of jasmonate (JA) signaling, plays important roles in enhancing cambium cell division and differentiation by mediating cytokinin signaling in poplar 84K (Populus alba × Populus glandulosa). PagJAZ5 is preferentially expressed in developing phloem and cambium, weakly in developing xylem cells. Overexpression (OE) of PagJAZ5m (insensitive to JA) increased cambium activity and xylem differentiation, while jaz mutants showed opposite results. Transcriptome analyses revealed that cytokinin oxidase/dehydrogenase (CKXs) and type-A response regulators (RRs) were downregulated in PagJAZ5m OE plants. The bioactive cytokinins were significantly increased in PagJAZ5m overexpressing plants and decreased in jaz5 mutants, compared with that in 84K plants. The PagJAZ5 directly interact with PagMYC2a/b and PagWOX4b. Further, we found that the PagRR5 is regulated by PagMYC2a and PagWOX4b and involved in the regulation of xylem development. Our results showed that PagJAZ5 can increase cambium activity and promote xylem differentiation through modulating cytokinin level and type-A RR during wood formation in poplar.

PdRabG3f interfered with gibberellin-mediated internode elongation and xylem developing in poplar

As a member of the small GTPases family, Rab GTPases play a key role in specifying transport pathways in the intracellular membrane trafficking system and are involved in plant growth and development. By quantitative trait locus (QTL) mapping, PdRabG3f was identified as a candidate gene associated with shoot height in a hybrid offspring of Populus deltoides 'Danhong' × Populus simonii 'Tongliao1'. PdRabG3f localized to the nucleus, endoplasmic reticulum and tonoplast and was primarily expressed in the xylem and cambium. Overexpression of PdRabG3f in Populus alba × Populus glandulosa (84 K poplar) had inhibitory effects on vertical and radical growth. In the transgenic lines, there were evident changes in the levels of 15 gibberellin (GA) derivatives, and the application of exogenous GA3 partially restored the phenotypes mediated by GAs deficiency. The interaction between PdRabG3f and RIC4, which was the GA-responsive factor, provided additional explanation for PdRabG3f's inhibitory effect on poplar growth. RNA-seq analysis revealed differentially expressed genes (DEGs) associated with cell wall, xylem, and gibberellin response. PdRabG3f interfering endogenous GAs levels in poplar might involve the participation of MYBs and ultimately affected internode elongation and xylem development. This study provides a potential mechanism for gibberellin-mediated regulation of plant growth through Rab GTPases.

Blue light photoreceptor cryptochrome 1 promotes wood formation and anthocyanin biosynthesis in Populus

Blue light photoreceptor cryptochrome 1 (CRY1) in herbaceous plants plays crucial roles in various developmental processes, including cotyledon expansion, hypocotyl elongation and anthocyanin biosynthesis. However, the function of CRY1 in perennial trees is unclear. In this study, we identified two ortholog genes of CRY1 (PagCRY1a and PagCRY1b) from Populus, which displayed high sequence similarity to Arabidopsis CRY1. Overexpression of PagCRY1 substantially inhibited plant growth and promoted secondary xylem development in Populus, while CRISPR/Cas9-mediated knockout of PagCRY1 enhanced plant growth and delayed secondary xylem development. Moreover, overexpression of PagCRY1 dramatically increased anthocyanin accumulation. The further analysis supported that PagCRY1 functions specifically in response to blue light. Taken together, our results demonstrated that modulating the expression of blue light photoreceptor CRY1 ortholog gene in Populus could significantly influence plant biomass production and the process of wood formation, laying a foundation for further investigating the light-regulated tree growth.

WUSCHEL-RELATED HOMEOBOX genes are crucial for normal vascular organization and wood formation in poplar

Vascular cambium in tree species is a cylindrical domain of meristematic cells that are responsible for producing secondary xylem (also called wood) inward and secondary phloem outward. The poplar (Populus trichocarpa) WUSCHEL (WUS)-RELATED HOMEOBOX (WOX) family members, PtrWUSa and PtrWOX13b, were previously shown to be expressed in vascular cambium and differentiating xylem cells in poplar stems, but their functions remain unknown. Here, we investigated roles of PtrWUSa, PtrWOX13b and their close homologs in vascular organization and wood formation. Expression analysis showed that like PtrWUSa and PtrWOX13b, their close homologs, PtrWUSb, PtrWUS4a/b and PtrWOX13a/c, were also expressed in vascular cambium and differentiating xylem cells in poplar stems. PtrWUSa also exhibited a high level of expression in developing phloem fibers. Expression of PtrWUSa fused with the dominant EAR repression domain (PtrWUSa-DR) in transgenic poplar caused retarded growth of plants with twisted stems and curled leaves and a severe disruption of vascular organization. In PtrWUSa-DR stems, a drastic proliferation of cells occurred in the phloem region between vascular cambium and phloem fibers and they formed islands of ectopic vascular tissues or phloem fiber-like sclerenchyma cells. A similar proliferation of cells was also observed in PtrWUSa-DR leaf petioles and midveins. On the other hand, overexpression of PtrWOX4a-DR caused ectopic formation of vascular bundles in the cortical region, and overexpression of PtrWOX13a-DR and PtrWOX13b-DR led to a reduction in wood formation without affecting vascular organization in transgenic poplar plants. Together, these findings indicate crucial roles of PtrWUSa and PtrWOX13a/b in regulating vascular organization and wood formation, which furthers our understanding of the functions of WOX genes in regulating vascular cambium activity in tree species.

Analyses of high spatial resolution datasets identify genes associated with multi-layered secondary cell wall thickening in Pinus bungeana

BACKGROUND AND AIMS

Secondary cell wall (SCW) thickening is a major cellular developmental stage determining wood structure and properties. Although the molecular regulation of cell wall deposition during tracheary element differentiation has been well established in primary growth systems, less is known about the gene regulatory processes involved in the multi-layered SCW thickening of mature trees.

METHODS

Using third-generation [long-read single-molecule real-time (SMRT)] and second-generation [short-read sequencing by synthesis (SBS)] sequencing methods, we established a Pinus bungeana transcriptome resource with comprehensive functional and structural annotation for the first time. Using these approaches, we generated high spatial resolution datasets for the vascular cambium, xylem expansion regions, early SCW thickening, late SCW thickening and mature xylem tissues of 71-year-old Pinus bungeana trees.

KEY RESULTS

A total of 79 390 non-redundant transcripts, 31 808 long non-coding RNAs and 5147 transcription factors were annotated and quantified in different xylem tissues at all growth and differentiation stages. Furthermore, using this high spatial resolution dataset, we established a comprehensive transcriptomic profile and found that members of the NAC, WRKY, SUS, CESA and LAC gene families are major players in early SCW formation in tracheids, whereas members of the MYB and LBD transcription factor families are highly expressed during late SCW thickening.

CONCLUSIONS

Our results provide new molecular insights into the regulation of multi-layered SCW thickening in conifers. The high spatial resolution datasets provided can serve as important gene resources for improving softwoods.

Transcriptional dynamics reveals the asymmetrical events underlying graft union formation in pecan (Carya illinoinensis)

Grafting is a widely used technique for pecan propagation; however, the background molecular events underlying grafting are still poorly understood. In our study, the graft partners during pecan [Carya illinoinensis (Wangenh.) K. Koch] graft union formation were separately sampled for RNA-seq, and the transcriptional dynamics were described via weighted gene co-expression network analysis. To reveal the main events underlying grafting, the correlations between modules and grafting traits were analyzed. Functional annotation showed that during the entire graft process, signal transduction was activated in the scion, while messenger RNA splicing was induced in the rootstock. At 2 days after grafting, the main processes occurring in the scion were associated with protein synthesis and processing, while the primary processes occurring in the rootstock were energy release-related. During the period of 7-14 days after grafting, defense response was a critical process taking place in the scion; however, the main process functioning in the rootstock was photosynthesis. From 22 to 32 days after grafting, the principal processes taking place in the scion were jasmonic acid biosynthesis and defense response, whereas the highly activated processes associated with the rootstock were auxin biosynthesis and plant-type secondary cell wall biogenesis. To further prove that the graft partners responded asymmetrically to stress, hydrogen peroxide contents as well as peroxidase and β-1,3-glucanase activities were detected, and the results showed that their levels were increased in the scion not the rootstock at certain time points after grafting. Our study reveals that the scion and rootstock might respond asymmetrically to grafting in pecan, and the scion was likely associated with stress response, while the rootstock was probably involved in energy supply and xylem bridge differentiation during graft union formation.

Comparative physiological, biochemical, metabolomic, and transcriptomic analyses reveal the formation mechanism of heartwood for Acacia melanoxylon

Acacia melanoxylon is well known as a valuable commercial tree species owing to its high-quality heartwood (HW) products. However, the metabolism and regulatory mechanism of heartwood during wood development remain largely unclear. In this study, both microscopic observation and content determination proved that total amount of starches decreased and phenolics and flavonoids increased gradually from sapwood (SW) to HW. We also obtained the metabolite profiles of 10 metabolites related to phenolics and flavonoids during HW formation by metabolomics. Additionally, we collected a comprehensive overview of genes associated with the biosynthesis of sugars, terpenoids, phenolics, and flavonoids using RNA-seq. A total of ninety-one genes related to HW formation were identified. The transcripts related to plant hormones, programmed cell death (PCD), and dehydration were increased in transition zone (TZ) than in SW. The results of RT-PCR showed that the relative expression level of genes and transcription factors was also high in the TZ, regardless of the horizontal or vertical direction of the trunk. Therefore, the HW formation took place in the TZ for A. melanoxylon from molecular level, and potentially connected to plant hormones, PCD, and cell dehydration. Besides, the increased expression of sugar and terpenoid biosynthesis-related genes in TZ further confirmed the close connection between terpenoid biosynthesis and carbohydrate metabolites of A. melanoxylon. Furthermore, the integrated analysis of metabolism data and RNA-seq data showed the key transcription factors (TFs) regulating flavonoids and phenolics accumulation in HW, including negative correlation TFs (WRKY, MYB) and positive correlation TFs (AP2, bZIP, CBF, PB1, and TCP). And, the genes and metabolites from phenylpropanoid and flavonoid metabolism and biosynthesis were up-regulated and largely accumulated in TZ and HW, respectively. The findings of this research provide a basis for comprehending the buildup of metabolites and the molecular regulatory processes of HW formation in A. melanoxylon.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Methylation of microRNA genes and its effect on secondary xylem development of stem in poplar

MicroRNAs (miRNAs) and DNA methylation are both vital regulators of gene expression. DNA methylation can affect the transcription of miRNAs, just like coding genes, through methylating the CpG islands in the gene regions of miRNAs. Although previous studies have shown that DNA methylation and miRNAs can each be involved in the process of wood formation, the relationship between the two has been relatively little studied in plant wood formation. Studies have shown that the second internode (IN2) (from top to bottom) of 3-month-old poplar trees can represent the primary stage of poplar stem development and IN8 can represent the secondary stage. There were also significant differences in DNA methylation patterns and miRNA expression patterns obtained from PS and SS. In this study, we first interactively analyzed methylation and miRNA sequencing data to identify 43 differentially expressed miRNAs regulated by differential methylation from the primary stage and secondary stage, which were found to be involved in multiple biological processes related to wood formation by enrichment analysis. In addition, six miRNA/target gene modules were finally identified as potentially involved in secondary xylem development of poplar stems through degradome sequencing and functional analysis. In conclusion, this study provides important reference information on the mechanism of interaction between different regulatory pathways of wood formation.

Wood of trees: Cellular structure, molecular formation, and genetic engineering

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Growth-regulating factor 15-mediated vascular cambium differentiation positively regulates wood formation in hybrid poplar (Populus alba × P. glandulosa)

Introduction

Hybrid poplars are industrial trees in China. An understanding of the molecular mechanism underlying wood formation in hybrid poplars is necessary for molecular breeding. Although the division and differentiation of vascular cambial cells is important for secondary growth and wood formation, the regulation of this process is largely unclear.

Methods

In this study, mPagGRF15 OE and PagGRF15-SRDX transgenic poplars were generated to investigate the function of PagGRF15. RNA-seq and qRT-PCR were conducted to analyze genome-wide gene expression, while ChIP‒seq and ChIP-PCR were used to identified the downstream genes regulated by PagGRF15.

Results and discussion

We report that PagGRF15 from hybrid poplar (Populus alba × P. glandulosa), a growth-regulating factor, plays a critical role in the regulation of vascular cambium activity. PagGRF15 was expressed predominantly in the cambial zone of vascular tissue. Overexpression of mPagGRF15 (the mutated version of GRF15 in the miR396 target sequence) in Populus led to decreased plant height and internode number. Further stem cross sections showed that the mPagGRF15 OE plants exhibited significant changes in vascular pattern with an increase in xylem and a reduction in phloem. In addition, cambium cell files were decreased in the mPagGRF15 OE plants. However, dominant suppression of the downstream genes of PagGRF15 using PagGRF15-SRDX showed an opposite phenotype. Based on the RNA-seq and ChIP-seq results, combining qRT-PCR and ChIP-PCR analysis, candidate genes, such as WOX4b, PXY and GID1.3, were obtained and found to be mainly involved in cambial activity and xylem differentiation. Accordingly, we speculated that PagGRF15 functions as a positive regulator mediating xylem differentiation by repressing the expression of the WOX4a and PXY genes to set the pace of cambial activity. In contrast, PagGRF15 mediated the GA signaling pathway by upregulating GID1.3 expression to stimulate xylem differentiation. This study provides valuable information for further studies on vascular cambium differentiation mechanisms and genetic improvement of the specific gravity of wood in hybrid poplars.

Ectopic Expression of PtrLBD39 Retarded Primary and Secondary Growth in Populus trichocarpa

Primary and secondary growth of trees are needed for increments in plant height and stem diameter, respectively, affecting the production of woody biomass for applications in timber, pulp/paper, and related biomaterials. These two types of growth are believed to be both regulated by distinct transcription factor (TF)-mediated regulatory pathways. Notably, we identified PtrLBD39, a highly stem phloem-specific TF in Populus trichocarpa and found that the ectopic expression of PtrLBD39 in P. trichocarpa markedly retarded both primary and secondary growth. In these overexpressing plants, the RNA-seq, ChIP-seq, and weighted gene co-expression network analysis (WGCNA) revealed that PtrLBD39 directly or indirectly regulates TFs governing vascular tissue development, wood formation, hormonal signaling pathways, and enzymes responsible for wood components. This regulation led to growth inhibition, decreased fibrocyte secondary cell wall thickness, and reduced wood production. Therefore, our study indicates that, following ectopic expression in P. trichocarpa, PtrLBD39 functions as a repressor influencing both primary and secondary growth.

A cytokinin response factor PtCRF1 is involved in the regulation of wood formation in poplar

Wood formation is a complex developmental process under the control of multiple levels of regulatory transcriptional network and hormone signals in trees. It is well known that cytokinin (CK) signaling plays an important role in maintaining the activity of the vascular cambium. The CK response factors (CRFs) encoding a subgroup of AP2 transcription factors have been identified to mediate the CK-dependent regulation in different plant developmental processes. However, the functions of CRFs in wood development remain unclear. Here, we characterized the function of PtCRF1, a CRF transcription factor isolated from poplar, in the process of wood formation. The PtCRF1 is preferentially expressed in secondary vasculature, especially in vascular cambium and secondary phloem, and encodes a transcriptional activator. Overexpression of PtCRF1 in transgenic poplar plants led to a significant reduction in the cell layer number of vascular cambium. The development of wood tissue was largely promoted in the PtCRF1-overexpressing lines, while it was significantly compromised in the CRISPR/Cas9-generated double mutant plants of PtCRF1 and its closest homolog PtCRF2. The RNA sequencing (RNA-seq) and quantitative reverse transcription PCR (RT-qPCR) analyses showed that PtCRF1 repressed the expression of the typical CK-responsive genes. Furthermore, bimolecular fluorescence complementation assays revealed that PtCRF1 competitively inhibits the direct interactions between histidine phosphotransfer proteins and type-B response regulator by binding to PtHP protein. Collectively, these results indicate that PtCRF1 negatively regulates CK signaling and is required for woody cell differentiation in poplar.

The microRNA7833-AUX6 module plays a critical role in wood development by modulating cellular auxin influx in Populus tomentosa

The critical role of auxin on secondary vascular development in woody plants has been demonstrated. The concentration gradient of endogenous indole-3-acetic acid and the cellular and molecular pathways contributing to the auxin-directed vascular organization and wood growth have been uncovered in recent decades. However, our understanding of the roles and regulations of auxin influx in wood formation in trees remains limited. Here, we reported that a microRNA, miR7833, participates in the negative regulation of stem cambial cell division and secondary xylem development in Populus tomentosa. The miR7833 is mainly expressed in the vascular cambium during stem radical growth and specifically targets and represses two AUX/LAX family auxin influx carriers, AUX5 and AUX6, in poplar. We further revealed that poplar AUX6, the most abundant miR7833 target in the stem, is preferentially enriched in the developing xylem and is a positive regulator for cell division and differentiation events during wood formation. Moreover, inhibition of auxin influx carriers by 1-naphthoxyacetic acids abolished the regulatory effects of miR7833 and AUX6 on secondary xylem formation in poplar. Our results revealed the essential roles of the miR7833-AUX6 module in regulating cellular events in secondary xylem development and demonstrated an auxin influx-dependent mechanism for wood formation in poplar.

Trees need closure too: Wound-induced secondary vascular tissue regeneration

Trees play a pivotal role in terrestrial ecosystems as well as being an important natural resource. These attributes are primarily associated with the capacity of trees to continuously produce woody tissue from the vascular cambium, a ring of stem cells located just beneath the bark. Long-lived trees are exposed to a myriad of biological and environmental stresses that may result in wounding, leading to a loss of bark and the underlying vascular cambium. This affects both wood formation and the quality of timber arising from the tree. In addition, the exposed wound site is a potential entry point for pathogens that cause disease. In response to wounding, trees have the capacity to regenerate lost or damaged tissues at this site. Investigating gene expression changes associated with different stages of wound healing reveals complex and dynamic changes in the activity of transcription factors, signalling pathways and hormone responses. In this review we summarise these data and discuss how they relate to our current understanding of vascular cambium formation and xylem differentiation during secondary growth. Based on this analysis, a model for wound healing that provides the conceptual foundations for future studies aimed at understanding this intriguing process is proposed.

Effects of ethephon on heartwood formation and related physiological indices of Dalbergia odorifera T. Chen

Introduction

Dalbergia odorifera T. Chen, known as fragrant rosewood, is a rare and endangered tree species. Studies have shown that plant growth regulators can effectively promote heartwood formation. This study aimed to investigate the effects of ethephon (ETH) on heartwood formation and the influence of ethephon and hydrogen peroxide (H2O2) on the physiological characteristics in D. odorifera.

Methods

D. odorifera branches underwent treatment with 2.5% plant growth regulators, including ETH, jasmonic acid (JA), salicylic acid (SA), abscisic acid (ABA), H2O2, and inhibitors such as ascorbic acid (AsA) to inhibit H2O2 synthesis, and (S) -trans 2-amino-4 - (2-aminoethoxy) -3-butene (AVG) to inhibit ethylene synthesis. After a 14-day period, we conducted an analysis to evaluate the impact of these plant growth regulators on elongation distance, vessel occlusion percentage, and trans-nerol content. Additionally, the effects of ETH and H2O2 on endogenous plant hormones, H2O2 content, soluble protein content, and enzyme activity were investigated within 0-48 h of treatment.

Results

After treatment with ETH for 14 days, the extension distance of the heartwood material was 15 cm, while the trans-nerolol content was 15 times that of the ABA group. ETH and H2O2 promoted endogenous ethylene synthesis; Ethylene content peaked at 6 and 18 h. The peak ethylene content in the ETH group was 68.07%, 12.89%, and 20.87% higher than the initial value of the H2O2 group and ddH2O group, respectively, and 29.64% higher than that in the AVG group. The soluble protein content and activity of related enzymes were significantly increased following ETH treatment.

Discussion

ETH exhibited the most impact on heartwood formation while not hindering tree growth. This treatment effectively triggered the production of endogenous ethylene in plants and enhanced the activity of essential enzymes involved in heartwood formation. These findings serve as a valuable reference for future investigations into heartwood formation.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar

Introduction

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood.

Methods

In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031.

Results

Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex.

Discussion

Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

Expression Quantitative Trait Locus of Wood Formation-Related Genes in Salix suchowensis

Shrub willows are widely planted for landscaping, soil remediation, and biomass production, due to their rapid growth rates. Identification of regulatory genes in wood formation would provide clues for genetic engineering of willows for improved growth traits on marginal lands. Here, we conducted an expression quantitative trait locus (eQTL) analysis, using a full sibling F1 population of Salix suchowensis, to explore the genetic mechanisms underlying wood formation. Based on variants identified from simplified genome sequencing and gene expression data from RNA sequencing, 16,487 eQTL blocks controlling 5505 genes were identified, including 2148 cis-eQTLs and 16,480 trans-eQTLs. eQTL hotspots were identified, based on eQTL frequency in genomic windows, revealing one hotspot controlling genes involved in wood formation regulation. Regulatory networks were further constructed, resulting in the identification of key regulatory genes, including three transcription factors (JAZ1, HAT22, MYB36) and CLV1, BAM1, CYCB2;4, CDKB2;1, associated with the proliferation and differentiation activity of cambium cells. The enrichment of genes in plant hormone pathways indicates their critical roles in the regulation of wood formation. Our analyses provide a significant groundwork for a comprehensive understanding of the regulatory network of wood formation in S. suchowensis.

Genome-wide identification and expression analysis of GDP-D-mannose pyrophosphorylase and KATANIN in Corymbia citriodora

The GDP-D-mannose pyrophosphorylase (GMP) and microtubule severing enzyme KATANIN (KTN) are crucial for wood formation. Although functional identification has been performed in Arabidopsis, few comprehensive studies have been conducted in forest trees. In this study, we discovered 8 CcGMP and 4 CcKTN genes by analyzing the whole genome sequence of Corymbia citriodora. The chromosomal location, genome synteny, phylogenetic relationship, protein domain, motif identification, gene structure, cis-acting regulatory elements, and protein-interaction of CcGMP and CcKTN were all investigated. KTN has just one pair of segmentally duplicated genes, while GMP has no duplication events. According to gene structure, two 5' UTRs were identified in CcGMP4. Furthermore, there is no protein-interaction between KTN and GMP. Based on real-time PCR, the expression of most genes showed a positive connection with DBH diameters. In addition, the expression of CcGMP4 and CcKTN4 genes were greater in different size tree, indicating that these genes are important in secondary xylem production. Overall, this findings will enhance our comprehension of the intricacy of CcGMP&CcKTN across diverse DBHs and furnish valuable insights for future functional characterization of specific genes in C. citriodora.

Combined metabolome and transcriptome analysis reveals a critical role of lignin biosynthesis and lignification in stem-like pneumatophore development of the mangrove Avicennia marina

MAIN CONCLUSION

Transcriptional and metabolic regulation of lignin biosynthesis and lignification plays crucial roles in Avicennia marina pneumatophore development, facilitating its adaptation to coastal habitats. Avicennia marina is a pioneer mangrove species in coastal wetland. To cope with the periodic intertidal flooding and hypoxia environment, this species has developed a complex and extensive root system, with its most unique feature being a pneumatophore with a distinct above- and below-ground morphology and vascular structure. However, the characteristics of pneumatophore lignification remain unknown. Studies comparing the anatomy among above-ground pneumatophore, below-ground pneumatophore, and feeding root have suggested that vascular structure development in the pneumatophore is more like the development of a stem than of a root. Metabolome and transcriptome analysis illustrated that the accumulation of syringyl (S) and guaiacyl (G) units in the pneumatophore plays a critical role in lignification of the stem-like structure. Fourteen differentially accumulated metabolites (DAMs) and 10 differentially expressed genes involved in the lignin biosynthesis pathway were targeted. To identify genes significantly associated with lignification, we analyzed the correlation between 14 genes and 8 metabolites and further built a co-expression network between 10 transcription factors (TFs), including 5 for each of MYB and NAC, and 23 enzyme-coding genes involved in lignin biosynthesis. 4-Coumarate-CoA ligase, shikimate/quinate hydroxycinnamoyl transferase, cinnamyl alcohol dehydrogenase, caffeic acid 3-O-methyltransferase, phenylalanine ammonia-lyase, and peroxidase were identified to be strongly correlated with these TFs. Finally, we examined 9 key candidate genes through quantitative real-time PCR to validate the reliability of transcriptome data. Together, our metabolome and transcriptome findings reveal that lignin biosynthesis and lignification regulate pneumatophore development in the mangrove species A. marina and facilitate its adaptation to coastal habitats.

The MPK6-LTF1L1 module regulates lignin biosynthesis in rice through a distinct mechanism from Populus LTF1

Lignin is a complex polymer that provides structural support and defense to plants. It is synthesized in the secondary cell walls of specialized cells. Through regulates its stability, LTF1 acts as a switch to control lignin biosynthesis in Populus, a dicot plant. However, how lignin biosynthesis is regulated in rice, a monocot plant, remains unclear. By employing genetic, cellular, and chemical approaches, we discovered that LTF1L1, a rice homolog of LTF1, regulates lignin biosynthesis through a distinct mechanism from Populus LTF1. Knockout of LTF1L1 increased lignin synthesis in the sclerenchyma cells of rice stems, while overexpression of LTF1L1 decreased it. LTF1L1 is phosphorylated by OsMPK6 at Ser169, which did not affect its stability but impaired its ability to repress the expression of lignin biosynthesis genes. This was supported by the non-phosphorylated mutant of LTF1L1 (LTF1L1S169A), which displayed a stronger repressive effect on lignin biosynthesis in both rice and Populus. Our findings reveal that LTF1L1 acts as a negative regulator of lignin biosynthesis via a distinct mechanism from that of LTF1 in Populus and highlight the evolutionary diversity in the regulation of lignin biosynthesis in plants.

Transcriptome and miRNAs Profiles Reveal Regulatory Network and Key Regulators of Secondary Xylem Formation in "84K" Poplar

Secondary xylem produced by stem secondary growth is the main source of tree biomass and possesses great economic and ecological value in papermaking, construction, biofuels, and the global carbon cycle. The secondary xylem formation is a complex developmental process, and the underlying regulatory networks and potential mechanisms are still under exploration. In this study, using hybrid poplar (Populus alba × Populus glandulosa clone 84K) as a model system, we first ascertained three representative stages of stem secondary growth and then investigated the regulatory network of secondary xylem formation by joint analysis of transcriptome and miRNAs. Notably, 7507 differentially expressed genes (DEGs) and 55 differentially expressed miRNAs (DEMs) were identified from stage 1 without initiating secondary growth to stage 2 with just initiating secondary growth, which was much more than those identified from stage 2 to stage 3 with obvious secondary growth. DEGs encoding transcription factors and lignin biosynthetic enzymes and those associated with plant hormones were found to participate in the secondary xylem formation. MiRNA-target analysis revealed that a total of 85 DEMs were predicted to have 2948 putative targets. Among them, PagmiR396d-PagGRFs, PagmiR395c-PagGA2ox1/PagLHW/PagSULTR2/PagPolyubiquitin 1, PagmiR482d-PagLAC4, PagmiR167e-PagbHLH62, and PagmiR167f/g/h-PagbHLH110 modules were involved in the regulating cambial activity and its differentiation into secondary xylem, cell expansion, secondary cell wall deposition, and programmed cell death. Our results give new insights into the regulatory network and mechanism of secondary xylem formation.

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.

Chromosome-scale genomes of five Hongmu species in Leguminosae

The Legume family (Leguminosae or Fabaceae), is one of the largest and economically important flowering plants. Heartwood, the core of a tree trunk or branch, is a valuable and renewable resource employed for centuries in constructing sturdy and sustainable structures. Hongmu refers to a category of precious timber trees in China, encompassing 29 woody species, primarily from the legume genus. Due to the lack of genome data, detailed studies on their economic and ecological importance are limited. Therefore, this study generates chromosome-scale assemblies of five Hongmu species in Leguminosae: Pterocarpus santalinus, Pterocarpus macrocarpus, Dalbergia cochinchinensis, Dalbergia cultrata, and Senna siamea, using a combination of short-reads, long-read nanopore, and Hi-C data. We obtained 623.86 Mb, 634.58 Mb, 700.60 Mb, 645.98 Mb, and 437.29 Mb of pseudochromosome level assemblies with the scaffold N50 lengths of 63.1 Mb, 63.7 Mb, 70.4 Mb, 61.1 Mb and 32.2 Mb for P. santalinus, P. macrocarpus, D. cochinchinensis, D. cultrata and S. siamea, respectively. These genome data will serve as a valuable resource for studying crucial traits, like wood quality, disease resistance, and environmental adaptation in Hongmu.

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.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Brassinosteroids affect wood development and properties of Fraxinus mandshurica

Introduction

Xylem development plays a crucial role in wood formation in woody plants. In recent years, there has been growing attention towards the impact of brassinosteroids (BRs) on this xylem development. In the present study, we evaluated the dynamic variation of xylem development in Fraxinus mandshurica (female parent, M8) and a novel interspecific hybrid F. mandshurica × Fraxinus sogdiana (1601) from May to August 2020.

Methods

We obtained RNA-Seq transcriptomes of three tissue types (xylem, phloem, and leaf) to identify the differences in xylem-differentially expressed genes (X-DEGs) and xylem-specifically expressed genes (X-SEGs) in M8 and 1601 variants. We then further evaluated these genes via weighted gene co-expression network analysis (WGCNA) alongside overexpressing FmCPD, a BR biosynthesis enzyme gene, in transient transgenic F. mandshurica.

Results

Our results indicated that the xylem development cycle of 1601 was extended by 2 weeks compared to that of M8. In addition, during the later wood development stages (secondary wall thickening) of 1601, an increased cellulose content (14%) and a reduced lignin content (11%) was observed. Furthermore, vessel length and width increased by 67% and 37%, respectively, in 1601 compared with those of M8. A total of 4589 X-DEGs were identified, including enzymes related to phenylpropane metabolism, galactose metabolism, BR synthesis, and signal transduction pathways. WGCNA identified hub X-SEGs involved in cellulose synthesis and BR signaling in the 1601 wood formation-related module (CESA8, COR1, C3H14, and C3H15); in contrast, genes involved in phenylpropane metabolism were significantly enriched in the M8 wood formation-related module (CCoAOMT and CCR). Moreover, overexpression of FmCPD in transient transgenic F. mandshurica affected the expression of genes associated with lignin and cellulose biosynthesis signal transduction. Finally, BR content was determined to be approximately 20% lower in the M8 xylem than in the 1601 xylem, and the exogenous application of BRs (24-epi brassinolide) significantly increased the number of xylem cell layers and altered the composition of the secondary cell walls in F. mandshurica.

Discussion

Our findings suggest that BR biosynthesis and signaling play a critical role in the differing wood development and properties observed between M8 and 1601 F. mandshurica.

A method for analyzing programmed cell death in xylem development by flow cytometry

Programmed cell death (PCD) is a genetically regulated developmental process leading to the death of specific types of plant cells, which plays important roles in plant development and growth such as wood formation. However, an efficient method needs to be established to study PCD in woody plants. Flow cytometry is widely utilized to evaluate apoptosis in mammalian cells, while it is rarely used to detect PCD in plants, especially in woody plants. Here, we reported that the xylem cell protoplasts from poplar stem were stained with a combination of fluorescein annexin V-FITC and propidium iodide (PI) and then sorted by flow cytometry. As expected, living cells (annexin V-FITC negative/PI negative), early PCD cells (annexin V-FITC positive/PI negative), and late PCD cells (annexin V-FITC positive/PI positive) could be finely distinguished through this method and then subjected for quantitative analysis. The expression of cell-type- and developmental stages-specific marker genes was consistent with the cell morphological observation. Therefore, the newly developed fluorescence-activated cell sorting (FACS) method can be used to study PCD in woody plants, which will be beneficial for studying the molecular mechanisms of wood formation.

PagUNE12 encodes a basic helix-loop-helix transcription factor that regulates the development of secondary vascular tissue in poplar

Secondary growth in woody plants generates new cells and tissues via the activity of the vascular cambium and drives the radial expansion of stems and roots. It is regulated by a series of endogenous factors, especially transcription factors. Here, we cloned the basic helix-loop-helix (bHLH) transcription factor gene UNFERTILIZED EMBRYO SAC12 (UNE12) from poplar (Populus alba × Populus glandulosa Uyeki) and used biochemical, molecular, and cytological assays to investigate the biological functions and regulatory mechanism of PagUNE12. PagUNE12 mainly localized in the nucleus and possessed transcriptional activation activity. It was widely expressed in vascular tissues, including primary phloem and xylem and secondary phloem and xylem. Poplar plants overexpressing PagUNE12 showed significantly reduced plant height, shorter internodes, and curled leaves compared with wild-type plants. Optical microscopy and transmission electron microscopy revealed that overexpressing PagUNE12 promoted secondary xylem development, with thicker secondary cell walls than wild-type poplar. Fourier transform infrared spectroscopy, confocal Raman microscopy, and 2D Heteronuclear Single Quantum Correlation analysis indicated that these plants also had increased lignin contents, with a lower relative abundance of syringyl lignin units and a higher relative abundance of guaiacyl lignin units. Therefore, overexpressing PagUNE12 promoted secondary xylem development and increased the lignin contents of secondary xylem in poplar, suggesting that this gene could be used to improve wood quality in the future.

PagERF81 regulates lignin biosynthesis and xylem cell differentiation in poplar

Lignin is a major component of plant cell walls and is essential for plant growth and development. Lignin biosynthesis is controlled by a hierarchical regulatory network involving multiple transcription factors. In this study, we showed that the gene encoding an APETALA 2/ethylene-responsive element binding factor (AP2/ERF) transcription factor, PagERF81, from poplar 84 K (Populus alba × P. glandulosa) is highly expressed in expanding secondary xylem cells. Two independent homozygous Pagerf81 mutant lines created by gene editing, produced significantly more but smaller vessel cells and longer fiber cells with more lignin in cell walls, while PagERF81 overexpression lines had less lignin, compared to non-transgenic controls. Transcriptome and reverse transcription quantitative PCR data revealed that multiple lignin biosynthesis genes including Cinnamoyl CoA reductase 1 (PagCCR1), Cinnamyl alcohol dehydrogenase 6 (PagCAD6), and 4-Coumarate-CoA ligase-like 9 (Pag4CLL9) were up-regulated in Pagerf81 mutants, but down-regulated in PagERF81 overexpression lines. In addition, a transient transactivation assay revealed that PagERF81 repressed the transcription of these three genes. Furthermore, yeast one hybrid and electrophoretic mobility shift assays showed that PagERF81 directly bound to a GCC sequence in the PagCCR1 promoter. No known vessel or fiber cell differentiation related genes were differentially expressed, so the smaller vessel cells and longer fiber cells observed in the Pagerf81 lines might be caused by abnormal lignin deposition in the secondary cell walls. This study provides insight into the regulation of lignin biosynthesis, and a molecular tool to engineer wood with high lignin content, which would contribute to the lignin-related chemical industry and carbon sequestration.

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.

Global Transcriptomic Analyses Provide New Insight into the Molecular Mechanisms of Endocarp Formation and Development in Iron Walnut (Juglans sigillata Dode)

Iron walnut (Juglans sigillata Dode) is a native species in southwestern China that exhibits variation in both fruit morphology and shell thickness. However, the underlying molecular processes controlling hardened endocarp development in walnut has not yet been reported. Here, we generated transcriptional profiles of iron walnut endocarp at three developmental stages using “Dapao”, the most common commercial variety. Using pairwise comparisons between these three stages, a total of 8555 non-redundant differentially expressed genes (DEGs) were identified, and more than one-half of the total DEGs exhibited significant differential expression in stage I as compared with stage II or stage III, suggesting that the first stage may ultimately determine the final characteristics of the mature walnut shell. Furthermore, in the clustering analysis of the above DEGs, 3682, 2349, and 2388 genes exhibited the highest expression in stages I, II, and III, respectively. GO enrichment analysis demonstrated that the major transcriptional variation among the three developmental stages was caused by differences in cell growth, plant hormones, metabolic process, and phenylpropanoid metabolism. Namely, using the tissue-specific expression analysis and a gene co-expression network, we identified MADS-box transcription factor JsiFBP2 and bHLH transcription factor JsibHLH94 as candidate regulators of endocarp formation in the early stage, and JsiNAC56 and JsiMYB78 might play key roles in regulating the lignification process of endocarp in the late stage. This study provides useful information for further research to dissect the molecular mechanisms governing the shell formation and development of iron walnut.

SCL14 Inhibits the Functions of the NAC043-MYB61 Signaling Cascade to Reduce the Lignin Content in Autotetraploid Populus hopeiensis

Whole-genome duplication often results in a reduction in the lignin content in autopolyploid plants compared with their diploid counterparts. However, the regulatory mechanism underlying variation in the lignin content in autopolyploid plants remains unclear. Here, we characterize the molecular regulatory mechanism underlying variation in the lignin content after the doubling of homologous chromosomes in Populus hopeiensis. The results showed that the lignin content of autotetraploid stems was significantly lower than that of its isogenic diploid progenitor throughout development. Thirty-six differentially expressed genes involved in lignin biosynthesis were identified and characterized by RNA sequencing analysis. The expression of lignin monomer synthase genes, such as PAL, COMT, HCT, and POD, was significantly down-regulated in tetraploids compared with diploids. Moreover, 32 transcription factors, including MYB61, NAC043, and SCL14, were found to be involved in the regulatory network of lignin biosynthesis through weighted gene co-expression network analysis. We inferred that SCL14, a key repressor encoding the DELLA protein GAI in the gibberellin (GA) signaling pathway, might inhibit the NAC043-MYB61 signaling functions cascade in lignin biosynthesis, which results in a reduction in the lignin content. Our findings reveal a conserved mechanism in which GA regulates lignin synthesis after whole-genome duplication; these results have implications for manipulating lignin production.

Identification and Expression Profile of CLE41/44-PXY-WOX Genes in Adult Trees Pinus sylvestris L. Trunk Tissues during Cambial Activity

WUSCHEL (WUS)-related homeobox (WOX) protein family members play important roles in the maintenance and proliferation of the stem cells in the cambium, the lateral meristem that forms all the wood structural elements. Most studies have examined the function of these genes in angiosperms, and very little was known about coniferous trees. Pine is one of the most critical forest-forming conifers globally, and in this research, we studied the distribution of WOX4, WOX13, and WOXG genes expression in Pinus sylvestris L. trunk tissues. Further, we considered the role of TDIF(CLE41/44)/TDR(PXY) signaling in regulating Scots pine cambial activity. The distribution of CLE41/44-PXY-WOXs gene expression in Scots pine trunk tissues was studied: (1) depending on the stage of ontogenesis (the first group of objects); and (2) depending on the stage of cambial growth (the second group of objects). The first group of objects is lingonberry pine forests of different ages (30-, 80-, and 180-year-old stands) in the middle taiga subzone. At the time of selection, all the trees of the studied groups were at the same seasonal stage of development: the formation of late phloem and early xylem was occurring in the trunk. The second group of objects is 40-year-old pine trees that were selected growing in the forest seed orchard. We took the trunk tissue samples on 27 May 2022, 21 June 2022, and 21 July 2022. We have indicated the spatial separation expressed of PsCLE41/44 and PsPXY in pine trunk tissues. PsCLE41/44 was differentially expressed in Fraction 1, including phloem cells and cambial zone. Maximum expression of the PsPXY gene occurred in Fraction 2, including differentiating xylem cells. The maximum expression of the PsCLE41/44 gene occurred on 27 May, when the number of cells in the cambial zone was the highest, and then it decreased to almost zero. The PsPXY gene transcript level increased from May to the end of July. We found that the highest transcript level of the PsWOX4 gene was during the period of active cell proliferation in the cambial zone, and also in the trees with the cambial age 63 years, which were characterized by the largest number of cell layers in the cambial zone. In this study, we have examined the expression profiles of genes belonging to the ancient clade (PsWOXG and PsWOX13) in stem tissues in Scots pine for the first time. We found that, in contrast to PsWOX4 (high expression that was observed during the period of active formation of early tracheids), the expression of genes of the ancient clade of the WOX genes was observed during the period of decreased cambial activity in the second half of the growing season. We found that PsWOX13 expression was shifted to Fraction 1 in most cases and increased from the phloem side, while PsWOXG expression was not clearly bound to a certain fraction. Based on the data, the role of the CLE41/44-PXY-WOX signaling module in regulating P. sylvestris cambial growth is discussed.

Transcriptomic Evidence Reveals Low Gelatinous Layer Biosynthesis in Neolamarckia cadamba after Gravistimulation

Trees can control their shape and resist gravity by producing tension wood (TW), which is a special wood that results from trees being put under stress. TW is characterized by the presence of a gelatinous layer (G layer) and the differential distribution of cell wall polymers. In this study, we investigated whether or not gravistimulation in N. cadamba resulted in TW with an obvious G layer. The results revealed an absence of an obvious G layer in samples of the upper side of a leaning stem (UW), as well as an accumulation of cellulose and a decrease in lignin content. A negligible change in the content of these polymers was recorded and compared to untreated plant (NW) samples, revealing the presence of a G layer either in much lower concentrations or in a lignified form. A transcriptomic investigation demonstrated a higher expression of cell wall esterase- and hydrolase-related genes in the UW, suggesting an accumulation of noncellulosic sugars in the UW, similar to the spectroscopy results. Furthermore, several G-layer-specific genes were also downregulated, including fasciclin-like arabinogalactan proteins (FLA), beta-galactosidase (BGAL) and chitinase-like proteins (CTL). The gene coexpression network revealed a strong correlation between cell-wall-synthesis-related genes and G-layer-synthesis-specific genes, suggesting their probable antagonistic role during G layer formation. In brief, the G layer in N. cadamba was either synthesized in a very low amount or was lignified during an early stage of growth; further experimental validation is required to understand the exact mechanism and stage of G layer formation in N. cadamba during gravistimulation.

Overexpression of PtoMYB115 improves lignocellulose recalcitrance to enhance biomass digestibility and bioethanol yield by specifically regulating lignin biosynthesis in transgenic poplar

Background

Woody plants provide the most abundant biomass resource that is convertible for biofuels. Since lignin is a crucial recalcitrant factor against lignocellulose hydrolysis, genetic engineering of lignin biosynthesis is considered as a promising solution. Many MYB transcription factors have been identified to involve in the regulation of cell wall formation or phenylpropanoid pathway. In a previous study, we identified that PtoMYB115 contributes to the regulation of proanthocyanidin pathway, however, little is known about its role in lignocellulose biosynthesis and biomass saccharification in poplar.

Results

Here, we detected the changes of cell wall features and examined biomass enzymatic saccharification for bioethanol production under various chemical pretreatments in PtoMYB115 transgenic plants. We reported that PtoMYB115 might specifically regulate lignin biosynthesis to affect xylem development. Overexpression of PtoMYB115 altered lignin biosynthetic gene expression, resulting in reduced lignin deposition, raised S/G and beta-O-4 linkage, resulting in a significant reduction in cellulase adsorption with lignin and an increment in cellulose accessibility. These alterations consequently improved lignocellulose recalcitrance for significantly enhanced biomass saccharification and bioethanol yield in the PtoMYB115-OE transgenic lines. In contrast, the knockout of PtoMYB115 by CRISPR/Cas9 showed reduced woody utilization under various chemical pretreatments.

Conclusions

This study shows that PtoMYB115 plays an important role in specifically regulating lignin biosynthesis and improving lignocellulose features. The enhanced biomass saccharification and bioethanol yield in the PtoMYB115-OE lines suggests that PtoMYB115 is a candidate gene for genetic modification to facilitate the utilization of biomass.

Transcriptomic analysis and physiological characteristics of exogenous naphthylacetic acid application to regulate the healing process of oriental melon grafted onto squash

The plant graft healing process is an intricate development influenced by numerous endogenous and environmental factors. This process involves the histological changes, physiological and biochemical reactions, signal transduction, and hormone exchanges in the grafting junction. Studies have shown that applying exogenous plant growth regulators can effectively promote the graft healing process and improve the quality of grafted plantlets. However, the physiological and molecular mechanism of graft healing formation remains unclear. In our present study, transcriptome changes in the melon and cucurbita genomes were analyzed between control and NAA treatment, and we provided the first view of complex networks to regulate graft healing under exogenous NAA application. The results showed that the exogenous NAA application could accelerate the graft healing process of oriental melon scion grafted onto squash rootstock through histological observation, increase the SOD, POD, PAL, and PPO activities during graft union development and enhance the contents of IAA, GA₃, and ZR except for the IL stage. The DEGs were identified in the plant hormone signal-transduction, phenylpropanoid biosynthesis, and phenylalanine metabolism through transcriptome analysis of CK vs. NAA at the IL, CA, and VB stage by KEGG pathway enrichment analysis. Moreover, the exogenous NAA application significantly promoted the expression of genes involved in the hormone signal-transduction pathway, ROS scavenging system, and vascular bundle formation.

The Regulation of Xylem Development by Transcription Factors and Their Upstream MicroRNAs

Xylem, as a unique organizational structure of vascular plants, bears water transport and supports functions necessary for plant survival. Notably, secondary xylem in the stem (i.e., wood) also has important economic and ecological value. In view of this, the regulation of xylem development has been widely concerned. In recent years, studies on model plants Arabidopsis and poplar have shown that transcription factors play important regulatory roles in various processes of xylem development, including the directional differentiation of procambium and cambium into xylem, xylem arrangement patterns, secondary cell wall formation and programmed cell death. This review focuses on the regulatory roles of widely and thoroughly studied HD-ZIP, MYB and NAC transcription factor gene families in xylem development, and it also pays attention to the regulation of their upstream microRNAs. In addition, the existing questions in the research and future research directions are prospected.

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.

Changes in the Activity of the CLE41/PXY/WOX Signaling Pathway in the Birch Cambial Zone under Different Xylogenesis Patterns

The balance between cell proliferation and differentiation into other cell types is crucial for meristem indeterminacy, and both growth aspects are under genetic control. The peptide-receptor signaling module regulates the activity of the cambial stem cells and the differentiation of their derivatives, along with cytokinins and auxin. We identified the genes encoding the signaling module CLE41-PXY and the regulator of vascular cambium division WOX4 and studied their expression during the period of cambial growth in the radial row: the conducting phloem/cambial zone and the differentiating xylem in two forms of Betula pendula, silver birch and Karelian birch. We have shown that the expression maximum of the BpCLE41/44a gene precedes the expression maximum of the BpPXY gene. Non-figured Karelian birch plants with straight-grained wood are characterized by a more intensive growth and the high expression of CLE41/44-PXY-WOX4. Figured Karelian birch plants, where the disturbed ratio and spatial orientation of structural elements characterizes the wood, have high levels of BpWOX4 expression and a decrease in xylem growth as well as the formation of xylem with a lower vessel density. The mutual influences of CLE41-PXY signaling and auxin signaling on WOX4 gene activity and the proliferation of cambium stem cells are discussed.

Effect of the Longitudinal Tensile Creep on the Stiffness of Radiata Pine (Pinus radiata D. Don)

The influence of load on the cellulose microfibrils of single cells or thin wood foils is known. It can decrease the cellulose microfibril angles and, in turn, increase the stiffness. However, this modification of a piece of wood, which is made up of multiple cells, is unknown. The aim of this research was to study the effect of tensile creep on the longitudinal stiffness of radiata pine wood. The modulus of elasticity of each specimen was determined before and after being subjected to tensile creep. The samples were loaded at 1170 N and 1530 N for 20 min at 70 °C. The load was determined as a function of a percentage of the force at the proportional limit. The moduli of elasticity before and post-tensile creep showed no effect on the stiffness of wood at the macroscopic level, but neither were there damage to the cell structure. It can be assumed that there are changes at the microscopic level, but they are not enough to be reflected at the macro scale. It is also challenging to achieve the modifications that occur at the level of a single cell or in thin wood foils; however, the implications of this would be favorable for the development of stronger wood-based products.

Overexpression of EgrIAA20 from Eucalyptus grandis, a Non-Canonical Aux/IAA Gene, Specifically Decouples Lignification of the Different Cell-Types in Arabidopsis Secondary Xylem

Wood (secondary xylem) formation is regulated by auxin, which plays a pivotal role as an integrator of developmental and environmental cues. However, our current knowledge of auxin-signaling during wood formation is incomplete. Our previous genome-wide analysis of Aux/IAAs in Eucalyptus grandis showed the presence of the non-canonical paralog member EgrIAA20 that is preferentially expressed in cambium. We analyzed its cellular localization using a GFP fusion protein and its transcriptional activity using transactivation assays, and demonstrated its nuclear localization and strong auxin response repressor activity. In addition, we functionally tested the role of EgrIAA20 by constitutive overexpression in Arabidopsis to investigate for phenotypic changes in secondary xylem formation. Transgenic Arabidopsis plants overexpressing EgrIAA20 were smaller and displayed impaired development of secondary fibers, but not of other wood cell types. The inhibition in fiber development specifically affected their cell wall lignification. We performed yeast-two-hybrid assays to identify EgrIAA20 protein partners during wood formation in Eucalyptus, and identified EgrIAA9A, whose ortholog PtoIAA9 in poplar is also known to be involved in wood formation. Altogether, we showed that EgrIAA20 is an important auxin signaling component specifically involved in controlling the lignification of wood fibers.

Eucalyptus grandis AUX/INDOLE-3-ACETIC ACID 13 (EgrIAA13) is a novel transcriptional regulator of xylogenesis

Our Induced Somatic Sector Analysis and protein-protein interaction experiments demonstrate that Eucalyptus grandis IAA13 regulates xylem fibre and vessel development, potentially via EgrIAA13 modules involving ARF2, ARF5, ARF6 and ARF19. Auxin is a crucial phytohormone regulating multiple aspects of plant growth and differentiation, including regulation of vascular cambium activity, xylogenesis and its responsiveness towards gravitropic stress. Although the regulation of these biological processes greatly depends on auxin and regulators of the auxin signalling pathway, many of their specific functions remain unclear. Therefore, the present study aims to functionally characterise Eucalyptus grandis AUX/INDOLE-3-ACETIC ACID 13 (EgrIAA13), a member of the auxin signalling pathway. In Eucalyptus and Populus, EgrIAA13 and its orthologs are preferentially expressed in the xylogenic tissues and downregulated in tension wood. Therefore, to further investigate EgrIAA13 and its function during xylogenesis, we conducted subcellular localisation and Induced Somatic Sector Analysis experiments using overexpression and RNAi knockdown constructs of EgrIAA13 to create transgenic tissue sectors on growing stems of Eucalyptus and Populus. Since Aux/IAAs interact with Auxin Responsive Factors (ARFs), in silico predictions of IAA13-ARF interactions were explored and experimentally validated via yeast-2-hybrid experiments. Our results demonstrate that EgrIAA13 localises to the nucleus and that downregulation of EgrIAA13 impedes Eucalyptus xylem fibre and vessel development. We also observed that EgrIAA13 interacts with Eucalyptus ARF2, ARF5, ARF6 and ARF19A. Based on these results, we conclude that EgrIAA13 is a regulator of Eucalyptus xylogenesis and postulate that the observed phenotypes are likely to result from alterations in the auxin-responsive transcriptome via IAA13-ARF modules such as EgrIAA13-EgrARF5. Our results provide the first insights into the regulatory role of EgrIAA13 during xylogenesis.

Integrated mRNA and Small RNA Sequencing Reveals microRNAs Associated With Xylem Development in Dalbergia odorifera

Dalbergia odorifera is a rare and precious rosewood specie, whose wood is a very high-quality material for valuable furniture and carving crafts. However, limited information is available about the process of wood formation in D. odorifera. To determine genes that might be closely associated with the xylem differentiation process, we analyzed the differentially expressed genes (DEGs) and microRNAs (miRNAs) from specific xylem tissues of D. odorifera by RNA sequencing (RNA-seq) and small RNA sequencing (small RNA-seq). In total, we obtained 134,221,955 clean reads from RNA-seq and 90,940,761 clean reads from small RNA-seq. By comparing the transition zone (Dotz) and sapwood (Dosw) samples, a total of 395 DEGs were identified. Further analysis revealed that DEGs encoded for WRKY transcription factors (eight genes), lignin synthesis (PER47, COMT, CCR2), cell wall composition (UXS2), gibberellin synthesis (KAO2, GA20OX1), jasmonic acid synthesis (OPR2, CYP74A), and synthesis of flavonoids (PAL2) and terpenoids (CYP71A1). Subsequently, a preliminary analysis by small RNA-seq showed that the expressions of 14 miRNAs (such as miR168a-5p, miR167f-5p, miR167h-5p, miR167e, miR390a, miR156g, novel_52, and novel_9) were significantly different between Dotz and Dosw. Further analysis revealed that the target genes of these differentially expressed miRNAs were enriched in the GO terms "amino acid binding," "cellulase activity," and "DNA beta-glucosyltransferase activity". Further, KEGG pathway annotation showed significant enrichment in "fatty acid elongation" and "biosynthesis of unsaturated fatty acids". These processes might be participating in the xylem differentiation of D. odorifera. Next, expression correlation analysis showed that nine differentially expressed miRNAs were significantly negatively associated with 21 target genes, which encoded for proteins such as pyrH, SPL6, SPL12, GCS1, and ARF8. Overall, this is the first study on miRNAs and their potential functions in the xylem development of D. odorifera, which provides a stepping stone for a detailed functional investigation of D. odorifera miRNAs.

Mobile forms of carbon in trees: metabolism and transport

Plants constitute 80% of the biomass on earth, and almost two-thirds of this biomass is found in wood. Wood formation is a carbon (C)-demanding process and relies on C transport from photosynthetic tissues. Thus, understanding the transport process is of major interest for understanding terrestrial biomass formation. Here, we review the molecules and mechanisms used to transport and allocate C in trees. Sucrose is the major form in which C is transported in plants, and it is found in the phloem sap of all tree species investigated so far. However, in several tree species, sucrose is accompanied by other molecules, notably polyols and the raffinose family of oligosaccharides. We describe the molecules that constitute each of these transport groups, and their distribution across different tree species. Furthermore, we detail the metabolic reactions for their synthesis, the mechanisms by which trees load and unload these compounds in and out of the vascular system, and how they are radially transported in the trunk and finally catabolized during wood formation. We also address a particular C recirculation process between phloem and xylem that occurs in trees during the annual cycle of growth and dormancy. A search of possible evolutionary drivers behind the diversity of C-carrying molecules in trees reveals no consistent differences in C transport mechanisms between angiosperm and gymnosperm trees. Furthermore, the distribution of C forms across species suggests that climate-related environmental factors will not explain the diversity of C transport forms. However, the consideration of C-transport mechanisms in relation to tree-rhizosphere coevolution deserves further attention. To conclude the review, we identify possible future lines of research in this field.

PtrHAT22, as a higher hierarchy regulator, coordinately regulates secondary cell wall component biosynthesis in Populus trichocarpa

Homeodomain-leucine zipper (HD-Zip) II transcription factors (TFs) have been reported to play vital roles in diverse biological processes of plants. However, it remains unclear whether HD-Zip II TFs regulate secondary cell wall (SCW) in woody plants. In this study, we performed the functional characterization of a Populus trichocarpa HD-Zip II TF, PtrHAT22, which encodes a nuclear localized transcription repressor predominantly expressing in secondary developing tissues. Overexpression of PtrHAT22 showed arrested growths, including reduced heights and diameters above the ground, small leaves, and decreased biomass. Meanwhile, the contents of lignin, cellulose, and thickness of SCW significantly decreased, whilst the content of hemicellulose obviously increased in PtrHAT22 transgenic poplar. The expressions of some wood-associated TFs and structural genes significantly changed accordingly with the alternations of SCW characteristics in PtrHAT22 transgenic poplar. Furthermore, PtrHAT22 directly repressed the promoter activities of PtrMYB20, PtrMYB28, and PtrCOMT2, and bind two cis-acting elements that were specifically enriched in their promoter regions. Taken together, our results suggested that PtrHAT22, as a higher hierarchy TF like PtrWNDs, exerted coordination regulation of poplar SCW component biosynthesis through directly and indirectly regulating structural genes and different hierarchy TFs of SCW formation network.

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.

Genome-Wide Identification of Long Non-Coding RNAs and Their Potential Functions in Poplar Growth and Phenylalanine Biosynthesis

Poplar is an important bioenergy tree species. lncRNAs play important roles in various biological regulatory processes, and their expression pattern is more tissue-specific than mRNAs. In this study, P. deltoides “Danhong” (Pd) and P. simonii “Tongliao1” (Ps) with different growth rates and wood quality were used as experimental materials, and the transcriptomes of their shoot apical meristem, xylem, and phloem were sequenced. Furthermore, high-throughput RNA sequencing analysis revealed that the expression patterns of genes and lncRNAs are different between the two genotypes. 6,355 lncRNAs were identified. Based on target prediction, lncRNAs and target genes were involved in ADP binding, oxidoreductase activity, phenylpropanoid biosynthesis, and cyanoamino acid metabolism. The DElncRNAs in two poplars were co-expressed with transcription factors and structural genes of lignin and flavonoid pathways. In addition, we found the potential target lncRNAs of miRNA. This result provides basic evidence for a better understanding of the regulatory role of lncRNAs in regulating phenylalanine molecular pathways and wood formation.

PagGRF12a interacts with PagGIF1b to regulate secondary xylem development through modulating PagXND1a expression in Populus alba × P. glandulosa

Growth-regulating factors (GRFs) are important regulators of plant development and growth, but their possible roles in xylem development in woody plants remain unclear. Here, we report that Populus alba × Papulus glandulosa PagGRF12a negatively regulates xylem development in poplar. PagGRF12a is expressed in vascular tissues. Compared to non-transgenic control plants, transgenic poplar plants overexpressing PagGRF12a exhibited reduced xylem width and plants with repressed expression of PagGRF12a exhibited increased xylem width. Xylem NAC domain 1 (XND1) encodes a NAC domain transcription factor that regulates xylem development and transcriptional analyses revealed that PagXND1a is highly upregulated in PagGRF12a-overexpressing plants and downregulated in PagGRF12a-suppressed plants, indicating that PagGRF12a may regulate xylem development through PagXND1a. Transient transcriptional assays and chromatin immunoprecipitation-polymerase chain reaction assays confirmed that PagGRF12a directly upregulates PagXND1a. In addition, PagGRF12a interacts with the GRF-Interacting Factor (GIF) PagGIF1b, and this interaction enhances the effects of PagGRF12a on PagXND1a. Our results indicate that PagGRF12a inhibits xylem development by upregulating the expression of PagXND1a.

Molecular Genetic Characteristics of Different Scenarios of Xylogenesis on the Example of Two Forms of Silver Birch Differing in the Ratio of Structural Elements in the Xylem

Silver birch (Betula pendula Roth) is an economically important species in Northern Europe. The current research focused on the molecular background of different xylogenesis scenarios in the birch trunks. The study objects were two forms of silver birch, silver birch trees, and Karelian birch trees; the latter form is characterized by the formation of two types of wood, non-figured (straight-grained) and figured, respectively, while it is currently not clear which factors cause this difference. We identified VND/NST/SND genes that regulate secondary cell wall biosynthesis in the birch genome and revealed differences in their expression in association with the formation of xylem with different ratios of structural elements. High expression levels of BpVND7 accompanied differentiation of the type of xylem which is characteristic of the species. At the same time, the appearance of figured wood was accompanied by the low expression levels of the VND genes and increased levels of expression of NST and SND genes. We identified BpARF5 as a crucial regulator of auxin-dependent vascular patterning and its direct target—BpHB8. A decrease in the BpARF5 level expression in differentiating xylem was a specific characteristic of both Karelian birch with figured and non-figured wood. Decreased BpARF5 level expression in non-figured trees accompanied by decreased BpHB8 and VND/NST/SND expression levels compared to figured Karelian birch trees. According to the results obtained, we suggested silver birch forms differing in wood anatomy as valuable objects in studying the regulation of xylogenesis.