AspWood: High-Spatial-Resolution Transcriptome Profiles Reveal Uncharacterized Modularity of Wood Formation in Populus tremula

The transcription factor PagLBD4 represses cell differentiation and secondary cell wall biosynthesis in Populus

LBD (LATERAL ORGAN BOUNDARIES DOMAIN) transcription factors are key regulators of plant growth and development. In this study, we functionally characterized the PagLBD4 gene in Populus (Populus alba × Populus glandulosa). Overexpression of PagLBD4 (PagLBD4OE) significantly repressed secondary xylem differentiation and secondary cell wall (SCW) deposition, while CRISPR/Cas9-mediated PagLBD4 knockout (PagLBD4KO) significantly increased secondary xylem differentiation and SCW deposition. Consistent with the functional analysis, gene expression analysis revealed that SCW biosynthesis pathways were significantly down-regulated in PagLBD4OE plants but up-regulated in PagLBD4KO plants. We also performed DNA affinity purification followed by sequencing (DAP-seq) to identify genes bound by PagLBD4. Integration of RNA sequencing (RNA-seq) and DAP-seq data identified 263 putative direct target genes (DTGs) of PagLBD4, including important regulatory genes for SCW biosynthesis, such as PagMYB103 and PagIRX12. Together, our results demonstrated that PagLBD4 is a repressor of secondary xylem differentiation and SCW biosynthesis in Populus, which possibly lead to the dramatic growth repression in PagLBD4OE plants.

Systems genetic analysis of lignin biosynthesis in Populus tremula

The genetic control underlying natural variation in lignin content and composition in trees is not fully understood. We performed a systems genetic analysis to uncover the genetic regulation of lignin biosynthesis in a natural 'SwAsp' population of aspen (Populus tremula) trees. We analyzed gene expression by RNA sequencing (RNA-seq) in differentiating xylem tissues, and lignin content and composition using Pyrolysis-GC-MS in mature wood of 268 trees from 99 genotypes. Abundant variation was observed for lignin content and composition, and genome-wide association study identified proteins in the pentose phosphate pathway and arabinogalactan protein glycosylation among the top-ranked genes that are associated with these traits. Variation in gene expression and the associated genetic polymorphism was revealed through the identification of 312 705 local and 292 003 distant expression quantitative trait loci (eQTL). A co-expression network analysis suggested modularization of lignin biosynthesis and novel functions for the lignin-biosynthetic CINNAMYL ALCOHOL DEHYDROGENASE 2 and CAFFEOYL-CoA O-METHYLTRANSFERASE 3. PHENYLALANINE AMMONIA LYASE 3 was co-expressed with HOMEOBOX PROTEIN 5 (HB5), and the role of HB5 in stimulating lignification was demonstrated in transgenic trees. The systems genetic approach allowed linking natural variation in lignin biosynthesis to trees´ responses to external cues such as mechanical stimulus and nutrient availability.

C2H2 zinc finger protein PagIDD15A regulates secondary wall thickening and lignin biosynthesis in poplar

Wood production is largely determined by the activity of cambial cell proliferation, and the secondary cell wall (SCW) thickening of xylem cells determines the wood property. In this study, we identified an INDETERMINATE DOMAIN (IDD) type C2H2 zinc finger transcription factor PagIDD15A as a regulator of wood formation in Populus alba × Populus glandulosa. Downregulation of PagIDD15A expression by RNA interference (RNAi) inhibited xylem development and xylem cell secondary wall thickening. RNA-seq analysis showed that PagPAL1, PagCCR2 and PagCCoAOMT1 were downregulated in the differentiating xylem of the PagIDD15A-RNAi transgenic plants, showing that PagIDD15A may regulate SCW biosynthesis through inhibiting lignin biosynthesis. The downregulation of PagVND6-B2, PagMYB10 and PagMYC4 and upregulation of PagWRKY12 in the differentiating xylem of RNAi transgenic plants suggest that PagIDD15A may also regulate these transcription factor (TF) genes to affect SCW thickening. RT-qPCR analysis in the phloem-cambium of RNAi transgenic demonstrates that PagIDD15A may regulate the expression of the genes associated with cell proliferation, including, PagSHR (SHORTROOT), PagSCR (SCARECROW), PagCYCD3;1 (CYCLIN D3;1) and PagSMR4 (SIAMESE-RELATED4), to affect the cambial activity. This study provides the knowledge of the IDD-type C2H2 zinc finger protein in regulating wood formation.

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

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

Providing biological context for GWAS results using eQTL regulatory and co-expression networks in Populus

Our study utilized genome-wide association studies (GWAS) to link nucleotide variants to traits in Populus trichocarpa, a species with rapid linkage disequilibrium decay. The aim was to overcome the challenge of interpreting statistical associations at individual loci without sufficient biological context, which often leads to reliance solely on gene annotations from unrelated model organisms. We employed an integrative approach that included GWAS targeting multiple traits using three individual techniques for lignocellulose phenotyping, expression quantitative trait loci (eQTL) analysis to construct transcriptional regulatory networks around each candidate locus and co-expression analysis to provide biological context for these networks, using lignocellulose biosynthesis in Populus trichocarpa as a case study. The research identified three candidate genes potentially involved in lignocellulose formation, including one previously recognized gene (Potri.005G116800/VND1, a critical regulator of secondary cell wall formation) and two genes (Potri.012G130000/AtSAP9 and Potri.004G202900/BIC1) with newly identified putative roles in lignocellulose biosynthesis. Our integrative approach offers a framework for providing biological context to loci associated with trait variation, facilitating the discovery of new genes and regulatory networks.

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.

Vascular cambium stem cells: past, present and future

Secondary xylem and phloem originate from a lateral meristem called the vascular cambium that consists of one to several layers of meristematic cells. Recent lineage tracing studies have shown that only one of the cambial cells in each radial cell file functions as the stem cell, capable of producing both secondary xylem and phloem. Here, we first review how phytohormones and signalling peptides regulate vascular cambium formation and activity. We then propose how the stem cell concept, familiar from apical meristems, could be applied to cambium studies. Finally, we discuss how this concept could set the basis for future research.

PagARGOS promotes low-lignin wood formation in poplar

Wood formation, which occurs mainly through secondary xylem development, is important not only for supplying raw material for the 'ligno-chemical' industry but also for driving the storage of carbon. However, the complex mechanisms underlying the promotion of xylem formation remain to be elucidated. Here, we found that overexpression of Auxin-Regulated Gene involved in Organ Size (ARGOS) in hybrid poplar 84 K (Populus alba × Populus tremula var. glandulosa) enlarged organ size. In particular, PagARGOS promoted secondary growth of stems with increased xylem formation. To gain further insight into how PagARGOS regulates xylem development, we further carried out yeast two-hybrid screening and identified that the auxin transporter WALLS ARE THIN1 (WAT1) interacts with PagARGOS. Overexpression of PagARGOS up-regulated WAT1, activating a downstream auxin response promoting cambial cell division and xylem differentiation for wood formation. Moreover, overexpressing PagARGOS caused not only higher wood yield but also lower lignin content compared with wild-type controls. PagARGOS is therefore a potential candidate gene for engineering fast-growing and low-lignin trees with improved biomass production.

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.

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.

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.

PagEXPA1 combines with PagCDKB2;1 to regulate plant growth and the elongation of fibers in Populus alba × Populus glandulosa

Expansins are important plant cell wall proteins. They can loosen and soften the cell walls and lead to wall extension and cell expansion. To investigate their role in wood formation and fiber elongation, the PagEXPA1 that highly expressed in cell differentiation and expansion tissues was cloned from 84K poplar (Populus alba × P. glandulosa). The subcellular localization showed that PagEXPA1 located in the cell wall and it was highly expressed in primary stems and young leaves. Compared with non-transgenic 84K poplar, overexpression of PagEXPA1 can promote plant-growth, lignification, and fiber cell elongation, while PagEXPA1 Cas9-editing mutant lines exhibited the opposite phenotype. Transcriptome analysis revealed that DEGs were mainly enriched in some important processes, which are associated with cell wall formation and cellulose synthesis. The protein interaction prediction and expression analysis showed that PagCDKB2:1 and PagEXPA1 might have an interaction relationship. The luciferase complementary assay and bimolecular fluorescence complementary assay validated that PagEXPA1 can combined with PagCDKB2;1. So they promoted the expansion of xylem vascular tissues and the development of poplar though participating in the regulation of cell division and differentiation by programming the cell-cycle. It provides good foundation for molecular breeding of fast-growing and high-quality poplar varieties.

TARGET OF MONOPTEROS: key transcription factors orchestrating plant development and environmental response

Plants have an incredible ability to sustain root and vascular growth after initiation of the embryonic root and the specification of vascular tissue in early embryos. Microarray assays have revealed that a group of transcription factors, TARGET OF MONOPTEROS (TMO), are important for embryonic root initiation in Arabidopsis. Despite the discovery of their auxin responsiveness early on, their function and mode of action remained unknown for many years. The advent of genome editing has accelerated the study of TMO transcription factors, revealing novel functions for biological processes such as vascular development, root system architecture, and response to environmental cues. This review covers recent achievements in understanding the developmental function and the genetic mode of action of TMO transcription factors in Arabidopsis and other plant species. We highlight the transcriptional and post-transcriptional regulation of TMO transcription factors in relation to their function, mainly in Arabidopsis. Finally, we provide suggestions for further research and potential applications in plant genetic engineering.

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.

Seasonal dynamics of cell-to-cell transport in angiosperm wood

This study describes the seasonal changes in cell-to-cell transport in three selected angiosperm tree species, Acer pseudoplatanus (maple), Fraxinus excelsior (ash), and Populus tremula × tremuloides (poplar), with an emphasis on the living wood component, xylem parenchyma cells (XPCs). We performed anatomical studies, dye loading through the vascular system, measurements of non-structural carbohydrate content, immunocytochemistry, inhibitory assays and quantitative real-time PCR to analyse the transport mechanisms and seasonal variations in wood. The abundance of membrane dye in wood varied seasonally along with seasonally changing tree phenology, cambial activity, and non-structural carbohydrate content. Moreover, dyes internalized in vessel-associated cells and 'trapped' in the endomembrane system are transported farther between other XPCs via plasmodesmata. Finally, various transport mechanisms based on clathrin-mediated and clathrin-independent endocytosis, and membrane transporters, operate in wood, and their involvement is species and/or season dependent. Our study highlights the importance of XPCs in seasonally changing cell-to-cell transport in both ring-porous (ash) and diffuse-porous (maple, poplar) tree species, and demonstrates the involvement of both endocytosis and plasmodesmata in intercellular communication in angiosperm wood.

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.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

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.

The effect of nitrogen source and levels on hybrid aspen tree physiology and wood formation

Nitrogen can be taken up by trees in the form of nitrate, ammonium and amino acids, but the influence of the different forms on tree growth and development is poorly understood in angiosperm species like Populus. We studied the effects of both organic and inorganic forms of nitrogen on growth and wood formation of hybrid aspen trees in experimental conditions that allowed growth under four distinct steady-state nitrogen levels. Increased nitrogen availability had a positive influence on biomass accumulation and the radial dimensions of both xylem vessels and fibers, and a negative influence on wood density. An optimal level of nitrogen availability was identified where increases in biomass accumulation outweighed decreases in wood density. None of these responses depended on the source of nitrogen except for shoot biomass accumulation, which was stimulated more by treatments complemented with nitrate than by ammonium alone or the organic source arginine. The most striking difference between the nitrogen sources was the effect on lignin composition, whereby the abundance of H-type lignin increased only in the presence of nitrate. The differential effect of nitrate is possibly related to the well-known role of nitrate as a signaling compound. RNA-sequencing revealed that while the lignin-biosynthetic genes did not significantly (FDR <0.01) respond to added NO3 - , the expression of several laccases, catalysing lignin polymerization, was dependent on N-availability. These results reveal a unique role of nitrate in wood formation and contribute to the knowledge basis for decision-making in utilizing hybrid aspen as a bioresource.

Evidence for poplar PtaPLATZ18 in the regulation of plant growth and vascular tissues development

Introduction

Plant A/T-rich protein and zinc-binding protein (PLATZ) are plant-specific transcription factors playing a role in plant development and stress response. To assess the role of PLATZs in vascular system development and wood formation in poplar, a functional study for PtaPLATZ18, whose expression was associated with the xylem, was carried out.

Methods

Poplar dominant repressor lines for PtaPLATZ18 were produced by overexpressing a PtaPLATZ18-SRDX fusion. The phenotype of three independent transgenic lines was evaluated at morphological, biochemical, and molecular levels and compared to the wild type.

Results

The PtaPLATZ18-SRDX lines showed increased plant height resulting from higher internode length. Besides, a higher secondary xylem thickness was also evidenced in these dominant repression lines as compared to the wild type suggesting an activation of cambial activity. A higher amount of lignin was evidenced within wood tissue as compared to the wild type, indicating an alteration in cell wall composition within xylem cell types. This latter phenotype was linked to an increased expression of genes involved in lignin biosynthesis and polymerization.

Discussion

The phenotype observed in the PtaPLATZ18-SRDX lines argues that this transcription factor targets key regulators of plant growth and vascular tissues development.

Development of modified CTAB and Trizol protocols to isolate high molecular weight (HMW) RNA from polyphenol and polysaccharides rich pigeonpea (Cajanuscajan (L.) Millsp

Pigeonpea (Cajanuscajan L.) is a legume crop that contains high levels of polyphenolic compounds and polysaccharides that become a hindrance in extracting good-quality and enough amount of RNA from its tissues. With the existing methods of RNA isolation, the phenolic compounds may co-precipitate or bind to the RNA giving false results. Therefore, in the present study, we have modified conventional CTAB and Trizol-based methods which resulted in good quality with the absorbance A260/A280 ratios in the range of 1.83 to 1.98 and A260/230 ratios in the range of 2.0-2.23, revealed RNA to be of high purity and free of contaminants. Both of the proposed protocols yielded a good quantity of RNA ranging from 289 to 422μg per gram of tissue. Distinctly visible bands of 28S and 18S rRNA were observed without degradation or smear, which indicated the presence of intact RNA. RT-PCR analysis showed that isolated RNA was quantitatively sufficient and compliant for the subsequent gene expression analysis.

Fertilization Weakens the Ecological Succession of Dissolved Organic Matter in Paddy Rice Rhizosphere Soil at the Molecular Level

Dissolved organic matter (DOM) is involved in numerous biogeochemical processes, and understanding the ecological succession of DOM is crucial for predicting its response to farming (e.g., fertilization) practices. Although plentiful studies have examined how fertilization practice affects the content of soil DOM, it remains unknown how long-term fertilization drives the succession of soil DOM over temporal scales. Here, we investigated the succession of DOM in paddy rice rhizosphere soils subjected to different long-term fertilization treatments (CK: no fertilization; NPK: inorganic fertilization; OM: organic fertilization) along with plant growth. Our results demonstrated that long-term fertilization significantly promoted the molecular chemodiversity of DOM, but it weakened the correlation between DOM composition and plant development. Time-decay analysis indicated that the DOM composition had a shorter halving time under CK treatment (94.7 days), compared to NPK (337.4 days) and OM (223.8 days) treatments, reflecting a lower molecular turnover rate of DOM under fertilization. Moreover, plant development significantly affected the assembly process of DOM only under CK, not under NPK and OM treatments. Taken together, our results demonstrated that long-term fertilization, especially inorganic fertilization, greatly weakens the ecological succession of DOM in the plant rhizosphere, which has a profound implication for understanding the complex plant-DOM interactions.

Flexure wood formation via growth reprogramming in hybrid aspen involves jasmonates and polyamines and transcriptional changes resembling tension wood development

Stem bending in trees induces flexure wood but its properties and development are poorly understood. Here, we investigated the effects of low-intensity multidirectional stem flexing on growth and wood properties of hybrid aspen, and on its transcriptomic and hormonal responses. Glasshouse-grown trees were either kept stationary or subjected to several daily shakes for 5 wk, after which the transcriptomes and hormones were analyzed in the cambial region and developing wood tissues, and the wood properties were analyzed by physical, chemical and microscopy techniques. Shaking increased primary and secondary growth and altered wood differentiation by stimulating gelatinous-fiber formation, reducing secondary wall thickness, changing matrix polysaccharides and increasing cellulose, G- and H-lignin contents, cell wall porosity and saccharification yields. Wood-forming tissues exhibited elevated jasmonate, polyamine, ethylene and brassinosteroids and reduced abscisic acid and gibberellin signaling. Transcriptional responses resembled those during tension wood formation but not opposite wood formation and revealed several thigmomorphogenesis-related genes as well as novel gene networks including FLA and XTH genes encoding plasma membrane-bound proteins. Low-intensity stem flexing stimulates growth and induces wood having improved biorefinery properties through molecular and hormonal pathways similar to thigmomorphogenesis in herbaceous plants and largely overlapping with the tension wood program of hardwoods.

The AP2/ERF transcription factor PtoERF15 confers drought tolerance via JA-mediated signaling in Populus

Drought stress is one of the major limiting factors for the growth and development of perennial trees. Xylem vessels act as the center of water conduction in woody species, but the underlying mechanism of its development and morphogenesis under water-deficient conditions remains elucidation. Here, we identified and characterized an osmotic stress-induced ETHYLENE RESPONSE FACTOR 15 (PtoERF15) and its target, PtoMYC2b, which was involved in mediating vessel size, density, and cell wall thickness in response to drought in Populus tomentosa. PtoERF15 is preferentially expressed in differentiating xylem of poplar stems. Overexpression of PtoERF15 contributed to stem water potential maintaining, thus promoting drought tolerance. RNA-Seq and biochemical analysis further revealed that PtoERF15 directly regulated PtoMYC2b, encoding a switch of JA signaling pathway. Additionally, our findings verify that three sets of homologous genes from NAC (NAM, ATAF1/2, and CUC2) gene family: PtoSND1-A1/A2, PtoVND7-1/7-2, and PtoNAC118/120, as the targets of PtoMYC2b, are involved in the regulation of vessel morphology in poplar. Collectively, our study provides molecular evidence for the involvement of the PtoERF15-PtoMYC2b transcription cascade in maintaining stem water potential through the regulation of xylem vessel development, ultimately improving drought tolerance in poplar.

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.

Regulation of secondary growth by poplar BLADE-ON-PETIOLE genes in Arabidopsis

BLADE-ON-PETIOLE (BOP) genes are essential regulators of vegetative and reproductive development in land plants. First characterized in Arabidopsis thaliana (Arabidopsis), members of this clade function as transcriptional co-activators by recruiting TGACG-motif binding (TGA) basic leucine zipper (bZIP) transcription factors. Highly expressed at organ boundaries, these genes are also expressed in vascular tissue and contribute to lignin biosynthesis during secondary growth. How these genes function in trees, which undergo extensive secondary growth to produce wood, remains unclear. Here, we investigate the functional conservation of BOP orthologs in Populus trichocarpa (poplar), a widely-used model for tree development. Within the poplar genome, we identified two BOP-like genes, PtrBPL1 and PtrBPL2, with abundant transcripts in stems. To assess their functions, we used heterologous assays in Arabidopsis plants. The promoters of PtrBPL1 and PtrBPL2, fused with a β-glucuronidase (GUS) reporter gene showed activity at organ boundaries and in secondary xylem and phloem. When introduced into Arabidopsis plants, PtrBPL1 and PtrBPL2 complemented leaf and flower patterning defects in bop1 bop2 mutants. Notably, Arabidopsis plants overexpressing PtrBPL1 and PtrBPL2 showed defects in stem elongation and the lignification of secondary tissues in the hypocotyl and stem. Finally, PtrBPL1 and PtrBPL2 formed complexes with TGA bZIP proteins in yeast. Collectively, our findings suggest that PtrBPL1 and PtrBPL2 are orthologs of Arabidopsis BOP1 and BOP2, potentially contributing to secondary growth regulation in poplar trees. This work provides a foundation for functional studies in trees.

Natural variation in the prolyl 4-hydroxylase gene PtoP4H9 contributes to perennial stem growth in Populus

Perennial trees must maintain stem growth throughout their entire lifespan to progressively increase in size as they age. The overarching question of the molecular mechanisms that govern stem perennial growth in trees remains largely unanswered. Here we deciphered the genetic architecture that underlies perennial growth trajectories using genome-wide association studies (GWAS) for measures of growth traits across years in a natural population of Populus tomentosa. By analyzing the stem growth trajectory, we identified PtoP4H9, encoding prolyl 4-hydroxylase 9, which is responsible for the natural variation in the growth rate of diameter at breast height (DBH) across years. Quantifying the dynamic genetic contribution of PtoP4H9 loci to stem growth showed that PtoP4H9 played a pivotal role in stem growth regulation. Spatiotemporal expression analysis showed that PtoP4H9 was highly expressed in cambium tissues of poplars of various ages. Overexpression and knockdown of PtoP4H9 revealed that it altered cell expansion to regulate cell wall modification and mechanical characteristics, thereby promoting stem growth in Populus. We showed that natural variation in PtoP4H9 occurred in a BASIC PENTACYSTEINE transcription factor PtoBPC1-binding promoter element controlling PtoP4H9 expression. The geographic distribution of PtoP4H9 allelic variation was consistent with the modes of selection among populations. Altogether, our study provides important genetic insights into dynamic stem growth in Populus, and we confirmed PtoP4H9 as a potential useful marker for breeding or genetic engineering of poplars.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Recent genome resequencing paraded COBRA-Like gene family roles in abiotic stress and wood formation in Poplar

A cell wall determines the mechanical properties of a cell, serves as a barrier against plant stresses, and allows cell division and growth processes. The COBRA-Like (COBL) gene family encodes a putative glycosylphosphatidylinositol (GPI)-anchored protein that controls cellulose deposition and cell progression in plants by contributing to the microfibril orientation of a cell wall. Despite being studied in different plant species, there is a dearth of the comprehensive global analysis of COBL genes in poplar. Poplar is employed as a model woody plant to study abiotic stresses and biomass production in tree research. Improved genome resequencing has enabled the comprehensive exploration of the evolution and functional capacities of PtrCOBLs (Poplar COBRA-Like genes) in poplar. Phylogeny analysis has discerned and classified PtrCOBLs into two groups resembling the Arabidopsis COBL family, and group I genes possess longer proteins but have fewer exons than group II. Analysis of gene structure and motifs revealed PtrCOBLs maintained a rather stable motif and exon-intron pattern across members of the same group. Synteny and collinearity analyses exhibited that the evolution of the COBL gene family was heavily influenced by gene duplication events. PtrCOBL genes have undergone both segmental duplication and tandem duplication, followed by purifying selection. Promotor analysis flaunted various phytohormone-, growth- and stress-related cis-elements (e.g., MYB, ABA, MeJA, SA, AuxR, and ATBP1). Likewise, 29 Ptr-miRNAs of 20 families were found targeting 11 PtrCOBL genes. PtrCOBLs were found localized at the plasma membrane and extracellular matrix, while gene ontology analysis showed their involvement in plant development, plant growth, stress response, cellulose biosynthesis, and cell wall biogenesis. RNA-seq datasets depicted the bulk of PtrCOBL genes expression being found in plant stem tissues and leaves, rendering mechanical strength and rejoinders to environmental cues. PtrCOBL2, 3, 10, and 11 manifested the highest expression in vasculature and abiotic stress, and resemblant expression trends were upheld by qRT-PCR. Co-expression network analysis identified PtrCOBL2 and PtrCOBL3 as hub genes across all abiotic stresses and wood developing tissues. The current study reports regulating roles of PtrCOBLs in xylem differentiating tissues, tension wood formation, and abiotic stress latency that lay the groundwork for future functional studies of the PtrCOBL genes in poplar breeding.

Comprehensive analysis of lncRNA-mRNA regulatory network in Populus associated with xylem development

Long noncoding RNAs (lncRNAs) play essential roles in numerous biological processes in plants, such as regulating the gene expression. However, only a few studies have looked into their potential functions in xylem development. High-throughput sequencing of P. euramericana 'Zhonglin46' developing and mature xylem was performed in this study. Through sequencing analysis, 14,028 putative lncRNA transcripts were identified, including 4525 differentially expressed lncRNAs (DELs). Additional research revealed that in mature xylem, a total of 2320 DELs were upregulated and 2205 were downregulated compared to developing xylem. Meanwhile, there were a total of 8122 differentially expressed mRNAs (DEMs) that were upregulated and 16,424 that were downregulated in mature xylem compared with developing xylem. The cis- and trans-target genes of DELs were analyzed for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, which indicated that these DELs participate in controlling the phenylpropanoid and lignin biosynthesis pathway as well as the starch and sucrose metabolism pathway. Among the cis-regulated DELs, LNC_006291, LNC_006292, and LNC_006532 all participate in regulating multiple HCT gene family membranes. As targets, POPTR_001G045900v3 (CCR2) and POPTR_018G063500v3 (SUS) both have only one cis-regulatory lncRNA, referred to as LNC_000057 and LNC_006212, respectively. Moreover, LNC_004484 and two DELs named LNC_008014 and LNC_010781 were revealed to be important nodes in the co-expression network of trans-lncRNAs and mRNAs associated to the lignin biosynthesis pathway and cellulose and xylan biosynthetic pathways, respectively. Finally, quantitative real-time PCR (qRT-PCR) was used to confirme 34 pairs of lncRNA-mRNA. Taken together, these findings may help to clarify the regulatory role that lncRNAs play in xylem development and wood formation.

Impact of xylan on field productivity and wood saccharification properties in aspen

Xylan that comprises roughly 25% of hardwood biomass is undesirable in biorefinery applications involving saccharification and fermentation. Efforts to reduce xylan levels have therefore been made in many species, usually resulting in improved saccharification. However, such modified plants have not yet been tested under field conditions. Here we evaluate the field performance of transgenic hybrid aspen lines with reduced xylan levels and assess their usefulness as short-rotation feedstocks for biorefineries. Three types of transgenic lines were tested in four-year field tests with RNAi constructs targeting either Populus GT43 clades B and C (GT43BC) corresponding to Arabidopsis clades IRX9 and IRX14, respectively, involved in xylan backbone biosynthesis, GATL1.1 corresponding to AtGALT1 involved in xylan reducing end sequence biosynthesis, or ASPR1 encoding an atypical aspartate protease. Their productivity, wood quality traits, and saccharification efficiency were analyzed. The only lines differing significantly from the wild type with respect to growth and biotic stress resistance were the ASPR1 lines, whose stems were roughly 10% shorter and narrower and leaves showed increased arthropod damage. GT43BC lines exhibited no growth advantage in the field despite their superior growth in greenhouse experiments. Wood from the ASPR1 and GT43BC lines had slightly reduced density due to thinner cell walls and, in the case of ASPR1, larger cell diameters. The xylan was less extractable by alkali but more hydrolysable by acid, had increased glucuronosylation, and its content was reduced in all three types of transgenic lines. The hemicellulose size distribution in the GALT1.1 and ASPR1 lines was skewed towards higher molecular mass compared to the wild type. These results provide experimental evidence that GATL1.1 functions in xylan biosynthesis and suggest that ASPR1 may regulate this process. In saccharification without pretreatment, lines of all three constructs provided 8-11% higher average glucose yields than wild-type plants. In saccharification with acid pretreatment, the GT43BC construct provided a 10% yield increase on average. The best transgenic lines of each construct are thus predicted to modestly outperform the wild type in terms of glucose yields per hectare. The field evaluation of transgenic xylan-reduced aspen represents an important step towards more productive feedstocks for biorefineries.

Novel candidate genes for lignin structure identified through genome-wide association study of naturally varying Populus trichocarpa

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (-log₁₀(p-value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

Field testing of transgenic aspen from large greenhouse screening identifies unexpected winners

Trees constitute promising renewable feedstocks for biorefinery using biochemical conversion, but their recalcitrance restricts their attractiveness for the industry. To obtain trees with reduced recalcitrance, large-scale genetic engineering experiments were performed in hybrid aspen blindly targeting genes expressed during wood formation and 32 lines representing seven constructs were selected for characterization in the field. Here we report phenotypes of five-year old trees considering 49 traits related to growth and wood properties. The best performing construct considering growth and glucose yield in saccharification with acid pretreatment had suppressed expression of the gene encoding an uncharacterized 2-oxoglutarate-dependent dioxygenase (2OGD). It showed minor changes in wood chemistry but increased nanoporosity and glucose conversion. Suppressed levels of SUCROSE SYNTHASE, (SuSy), CINNAMATE 4-HYDROXYLASE (C4H) and increased levels of GTPase activating protein for ADP-ribosylation factor ZAC led to significant growth reductions and anatomical abnormalities. However, C4H and SuSy constructs greatly improved glucose yields in saccharification without and with pretreatment, respectively. Traits associated with high glucose yields were different for saccharification with and without pretreatment. While carbohydrates, phenolics and tension wood contents positively impacted the yields without pretreatment and growth, lignin content and S/G ratio were negative factors, the yields with pretreatment positively correlated with S lignin and negatively with carbohydrate contents. The genotypes with high glucose yields had increased nanoporosity and mGlcA/Xyl ratio, and some had shorter polymers extractable with subcritical water compared to wild-type. The pilot-scale industrial-like pretreatment of best-performing 2OGD construct confirmed its superior sugar yields, supporting our strategy.

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.

The PtoTCP20-miR396d-PtoGRF15 module regulates secondary vascular development in Populus

Secondary vascular development is a key biological characteristic of woody plants and the basis of wood formation. Our understanding of gene expression regulation and dynamic changes in microRNAs (miRNAs) during secondary vascular development is still limited. Here we present an integrated analysis of the miRNA and mRNA transcriptome of six phase-specific tissues-the shoot apex, procambium, primary vascular tissue, cambium, secondary phloem, and secondary xylem-in Populus tomentosa. Several novel regulatory modules, including the PtoTCP20-miR396d-PtoGRF15 module, were identified during secondary vascular development in Populus. A series of biochemical and molecular experiments confirmed that PtoTCP20 activated transcription of the miR396d precursor gene and that miR396d targeted PtoGRF15 to downregulate its expression. Plants overexpressing miR396d (35S:miR396d) showed enhanced secondary growth and increased xylem production. Conversely, during the transition from primary to secondary vascular development, plants with downregulated PtoTCP20expression (PtoTCP20-SRDX), downregulated miR396 expression (35S:STTM396), and PtoGRF15 overexpression (35S:PtoGRF15) showed delayed secondary growth. Novel regulatory modules were identified by integrated analysis of the miRNA and mRNA transcriptome, and the regulatory role of the PtoTCP20-miR396d-PtoGRF15 signaling cascade in secondary vascular development was validated in Populus, providing information to support improvements in forest cultivation and wood properties.

Different combinations of laccase paralogs nonredundantly control the amount and composition of lignin in specific cell types and cell wall layers in Arabidopsis

Vascular plants reinforce the cell walls of the different xylem cell types with lignin phenolic polymers. Distinct lignin chemistries differ between each cell wall layer and each cell type to support their specific functions. Yet the mechanisms controlling the tight spatial localization of specific lignin chemistries remain unclear. Current hypotheses focus on control by monomer biosynthesis and/or export, while cell wall polymerization is viewed as random and nonlimiting. Here, we show that combinations of multiple individual laccases (LACs) are nonredundantly and specifically required to set the lignin chemistry in different cell types and their distinct cell wall layers. We dissected the roles of Arabidopsis thaliana LAC4, 5, 10, 12, and 17 by generating quadruple and quintuple loss-of-function mutants. Loss of these LACs in different combinations led to specific changes in lignin chemistry affecting both residue ring structures and/or aliphatic tails in specific cell types and cell wall layers. Moreover, we showed that LAC-mediated lignification has distinct functions in specific cell types, waterproofing fibers, and strengthening vessels. Altogether, we propose that the spatial control of lignin chemistry depends on different combinations of LACs with nonredundant activities immobilized in specific cell types and cell wall layers.

Single-cell transcriptomics unveils xylem cell development and evolution

BACKGROUND

Xylem, the most abundant tissue on Earth, is responsible for lateral growth in plants. Typical xylem has a radial system composed of ray parenchyma cells and an axial system of fusiform cells. In most angiosperms, fusiform cells comprise vessel elements for water transportation and libriform fibers for mechanical support, while both functions are performed by tracheids in other vascular plants such as gymnosperms. Little is known about the developmental programs and evolutionary relationships of these xylem cell types.

RESULTS

Through both single-cell and laser capture microdissection transcriptomic profiling, we determine the developmental lineages of ray and fusiform cells in stem-differentiating xylem across four divergent woody angiosperms. Based on cross-species analyses of single-cell clusters and overlapping trajectories, we reveal highly conserved ray, yet variable fusiform, lineages across angiosperms. Core eudicots Populus trichocarpa and Eucalyptus grandis share nearly identical fusiform lineages, whereas the more basal angiosperm Liriodendron chinense has a fusiform lineage distinct from that in core eudicots. The tracheids in the basal eudicot Trochodendron aralioides, an evolutionarily reversed trait, exhibit strong transcriptomic similarity to vessel elements rather than libriform fibers.

CONCLUSIONS

This evo-devo framework provides a comprehensive understanding of the formation of xylem cell lineages across multiple plant species spanning over a hundred million years of evolutionary history.

Genome-Wide Identification and Characterization of Auxin Response Factor (ARF) Gene Family Involved in Wood Formation and Response to Exogenous Hormone Treatment in Populus trichocarpa

Auxin is a key regulator that virtually controls almost every aspect of plant growth and development throughout its life cycle. As the major components of auxin signaling, auxin response factors (ARFs) play crucial roles in various processes of plant growth and development. In this study, a total of 35 PtrARF genes were identified, and their phylogenetic relationships, chromosomal locations, synteny relationships, exon/intron structures, cis-elements, conserved motifs, and protein characteristics were systemically investigated. We also analyzed the expression patterns of these PtrARF genes and revealed that 16 of them, including PtrARF1, 3, 7, 11, 13-17, 21, 23, 26, 27, 29, 31, and 33, were preferentially expressed in primary stems, while 15 of them, including PtrARF2, 4, 6, 9, 10, 12, 18-20, 22, 24, 25, 28, 32, and 35, participated in different phases of wood formation. In addition, some PtrARF genes, with at least one cis-element related to indole-3-acetic acid (IAA) or abscisic acid (ABA) response, responded differently to exogenous IAA and ABA treatment, respectively. Three PtrARF proteins, namely PtrARF18, PtrARF23, and PtrARF29, selected from three classes, were characterized, and only PtrARF18 was a transcriptional self-activator localized in the nucleus. Moreover, Y2H and bimolecular fluorescence complementation (BiFC) assay demonstrated that PtrARF23 interacted with PtrIAA10 and PtrIAA28 in the nucleus, while PtrARF29 interacted with PtrIAA28 in the nucleus. Our results provided comprehensive information regarding the PtrARF gene family, which will lay some foundation for future research about PtrARF genes in tree development and growth, especially the wood formation, in response to cellular signaling and environmental cues.

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.

Unleashing the power within short-read RNA-seq for plant research: Beyond differential expression analysis and toward regulomics

RNA-seq has become a state-of-the-art technique for transcriptomic studies. Advances in both RNA-seq techniques and the corresponding analysis tools and pipelines have unprecedently shaped our understanding in almost every aspects of plant sciences. Notably, the integration of huge amount of RNA-seq with other omic data sets in the model plants and major crop species have facilitated plant regulomics, while the RNA-seq analysis has still been primarily used for differential expression analysis in many less-studied plant species. To unleash the analytical power of RNA-seq in plant species, especially less-studied species and biomass crops, we summarize recent achievements of RNA-seq analysis in the major plant species and representative tools in the four types of application: (1) transcriptome assembly, (2) construction of expression atlas, (3) network analysis, and (4) structural alteration. We emphasize the importance of expression atlas, coexpression networks and predictions of gene regulatory relationships in moving plant transcriptomes toward regulomics, an omic view of genome-wide transcription regulation. We highlight what can be achieved in plant research with RNA-seq by introducing a list of representative RNA-seq analysis tools and resources that are developed for certain minor species or suitable for the analysis without species limitation. In summary, we provide an updated digest on RNA-seq tools, resources and the diverse applications for plant research, and our perspective on the power and challenges of short-read RNA-seq analysis from a regulomic point view. A full utilization of these fruitful RNA-seq resources will promote plant omic research to a higher level, especially in those less studied species.

Aspartic proteases modulate programmed cell death and secondary cell wall synthesis during wood formation in poplar

Programmed cell death (PCD) is essential for wood development in trees. However, the determination of crucial factors involved in xylem PCD of wood development is still lacking. Here, two Populus trichocarpa typical aspartic protease (AP) genes, AP17 and AP45, modulate xylem maturation, especially fibre PCD, during wood formation. AP17 and AP45 were dominantly expressed in the fibres of secondary xylem, as suggested by GUS expression in APpro::GUS transgenic plants. Cas9/gRNA-induced AP17 or AP45 mutants delayed secondary xylem fibre PCD, and ap17ap45 double mutants showed more serious defects. Conversely, AP17 overexpression caused premature PCD in secondary xylem fibres, indicating a positive modulation in wood fibre PCD. Loss of AP17 and AP45 did not alter wood fibre wall thickness, whereas the ap17ap45 mutants showed a low lignin content in wood. However, AP17 overexpression led to a significant decrease in wood fibre wall thickness and lignin content, revealing the involvement in secondary cell wall synthesis during wood formation. In addition, the ap17ap45 mutant and AP17 overexpression plants resulted in a significant increase in saccharification yield in wood. Overall, AP17 and AP45 are crucial modulators in xylem maturation during wood development, providing potential candidate genes for engineering lignocellulosic wood for biofuel utilization.

PagDET2 promotes cambium cell division and xylem differentiation in poplar stem

Secondary growth of the woody tree stem is governed by meristematic cell division and differentiation in the vascular cambium. Multiple hormonal signals and endogenous developmental programs regulate vascular cambium activity. Brassinosteroids (BRs) significantly promote secondary stem growth and wood formation in poplar trees. However, the underlying regulatory mechanisms of BRs within the vascular tissue remain unclear. Genetic and anatomical approaches were used here to elucidate the role of PagDET2, the rate-limiting enzyme for BRs biosynthesis, in regulating secondary vascular cambium activity in Populus. This study showed that the elevated endogenous castasterone (CS) levels in tree stems through overexpressing PagDET2 could enhance cambium meristem cell activity and xylem (XY) differentiation to promote secondary stem growth. RNA-seq analysis revealed that genes involved in BRs response, vascular cambium cell division, XY differentiation, and secondary cell wall synthesis were up-regulated.

Reduced chemodiversity suppresses rhizosphere microbiome functioning in the mono-cropped agroecosystems

BACKGROUND

Rhizodeposits regulate rhizosphere interactions, processes, nutrient and energy flow, and plant-microbe communication and thus play a vital role in maintaining soil and plant health. However, it remains unclear whether and how alteration in belowground carbon allocation and chemodiversity of rhizodeposits influences microbiome functioning in the rhizosphere ecosystems. To address this research gap, we investigated the relationship of rhizosphere carbon allocation and chemodiversity with microbiome biodiversity and functioning during peanut (Arachis hypogaea) continuous mono-cropping. After continuously labeling plants with ¹³CO₂, we studied the chemodiversity and composition of rhizodeposits, along with the composition and diversity of active rhizosphere microbiome using metabolomic, amplicon, and shotgun metagenomic sequencing approaches based on DNA stable-isotope probing (DNA-SIP).

RESULTS

Our results indicated that enrichment and depletion of rhizodeposits and active microbial taxa varied across plant growth stages and cropping durations. Specifically, a gradual decrease in the rhizosphere carbon allocation, chemodiversity, biodiversity and abundance of plant-beneficial taxa (such as Gemmatimonas, Streptomyces, Ramlibacter, and Lysobacter), and functional gene pathways (such as quorum sensing and biosynthesis of antibiotics) was observed with years of mono-cropping. We detected significant and strong correlations between rhizodeposits and rhizosphere microbiome biodiversity and functioning, though these were regulated by different ecological processes. For instance, rhizodeposits and active bacterial communities were mainly governed by deterministic and stochastic processes, respectively. Overall, the reduction in carbon deposition and chemodiversity during peanut continuous mono-cropping tended to suppress microbial biodiversity and its functions in the rhizosphere ecosystem.

CONCLUSIONS

Our results, for the first time, provide the evidence underlying the mechanism of rhizosphere microbiome malfunctioning in mono-cropped systems. Our study opens new avenues to deeply disentangle the complex plant-microbe interactions from the perspective of rhizodeposits chemodiversity and composition and will serve to guide future microbiome research for improving the functioning and services of soil ecosystems. Video abstract.

AUXIN RESPONSE FACTOR7 integrates gibberellin and auxin signaling via interactions between DELLA and AUX/IAA proteins to regulate cambial activity in poplar

Cambial development in the stems of perennial woody species is rigorously regulated by phytohormones. Auxin and gibberellin (GA) play crucial roles in stimulating cambial activity in poplar (Populus spp.). In this study, we show that the DELLA protein REPRESSOR of ga1-3 Like 1 (RGL1), AUXIN RESPONSE FACTOR 7 (ARF7), and Aux/INDOLE-3-ACETIC ACID 9 (IAA9) form a ternary complex that mediates crosstalk between the auxin and GA signaling pathways in poplar stems during cambial development. Biochemical analysis revealed that ARF7 physically interacts with RGL1 and IAA9 through distinct domains. The arf7 loss-of-function mutant showed markedly attenuated responses to auxin and GA, whereas transgenic poplar plants overexpressing ARF7 displayed strongly improved cambial activity. ARF7 directly binds to the promoter region of the cambial stem cell regulator WOX4 to modulate its expression, thus integrating auxin and GA signaling to regulate cambial activity. Furthermore, the direct activation of PIN-FORMED 1 expression by ARF7 in the RGL1-ARF7-IAA9 module increased GA-dependent cambial activity via polar auxin transport. Collectively, these findings reveal that the crosstalk between auxin and GA signaling mediated by the RGL1-ARF7-IAA9 module is crucial for the precise regulation of cambial development in poplar.

Overexpression of PeNAC122 gene promotes wood formation and tolerance to osmotic stress in poplars

Finding the adequate balance between wood formation and abiotic stress resistance is still an important challenge for industrial woody crops. In this study, PeNAC122, a member of the NAC transcription factor (TF) family highly expressed in xylem, was cloned from Populus euphratica. Tissue expression and β-glucuronidase (GUS) staining showed that PeNAC122 was exclusively expressed in phloem fiber and secondary xylem of stems. Subcellular and yeast transactivation assays confirmed that PeNAC122 protein existed in the nucleus and did not have transcriptional activation and inhibitory activity. Overexpression of PeNAC122 poplar lines exhibited reduced plant height, thickened xylem, and accumulated lignin content in stems, and also upregulates the expression of secondary cell wall biosynthetic genes. Moreover, overexpression of PeNAC122 lines displayed more tolerance to PEG6000-induced osmotic stress, with stronger photosynthetic performance, higher antioxidant enzyme activity, and less accumulation of reactive oxygen species in leaves, and higher expression levels of stress response genes DREB2A, RD29, and NCED3. These results indicate that PeNAC122 plays a crucial role in wood formation and abiotic stress tolerance, which, in addition to potential use in improving wood quality, provides further insight into the role of NAC family TFs in balancing wood development and abiotic stress resistance.

Molecular basis of differential adventitious rooting competence in poplar genotypes

Recalcitrant adventitious root (AR) development is a major hurdle in propagating commercially important woody plants. Although significant progress has been made to identify genes involved in subsequent steps of AR development, the molecular basis of differences in apparent recalcitrance to form AR between easy-to-root and difficult-to-root genotypes remains unknown. To address this, we generated cambium tissue-specific transcriptomic data from stem cuttings of hybrid aspen, T89 (difficult-to-root) and hybrid poplar OP42 (easy-to-root), and used transgenic approaches to verify the role of several transcription factors in the control of adventitious rooting. Increased peroxidase activity was positively correlated with better rooting. We found differentially expressed genes encoding reactive oxygen species scavenging proteins to be enriched in OP42 compared with T89. A greater number of differentially expressed transcription factors in cambium cells of OP42 compared with T89 was revealed by a more intense transcriptional reprograming in the former. PtMYC2, a potential negative regulator, was less expressed in OP42 compared with T89. Using transgenic approaches, we demonstrated that PttARF17.1 and PttMYC2.1 negatively regulate adventitious rooting. Our results provide insights into the molecular basis of genotypic differences in AR and implicate differential expression of the master regulator MYC2 as a critical player in this process.

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.

PdeHCA2 affects biomass in Populus by regulating plant architecture, the transition from primary to secondary growth, and photosynthesis

PdeHCA2 regulates the transition from primary to secondary growth, plant architecture, and affects photosynthesis by targeting PdeBRC1 and controlling the anatomy of the mesophyll, and intercellular space, respectively. Branching, secondary growth, and photosynthesis are vital developmental processes of woody plants that determine plant architecture and timber yield. However, the mechanisms underlying these processes are unknown. Here, we report that the Populus transcription factor High Cambium Activity 2 (PdeHCA2) plays a role in the transition from primary to secondary growth, vascular development, and branching. In Populus, PdeHCA2 is expressed in undifferentiated provascular cells during primary growth, in phloem cells during secondary growth, and in leaf veins, which is different from the expression pattern of its homolog in Arabidopsis. Overexpression of PdeHCA2 has pleiotropic effects on shoot and leaf development; overexpression lines showed delayed growth of shoots and leaves, reduced photosynthesis, and abnormal shoot branching. In addition, auxin-, cytokinin-, and photosynthesis-related genes were differentially regulated in these lines. Electrophoretic mobility shift assays and transcriptome analysis indicated that PdeHCA2 directly up-regulates the expression of BRANCHED1 and the MADS-box gene PdeAGL9, which regulate plant architecture, by binding to cis-elements in the promoters of these genes. Taken together, our findings suggest that HCA2 regulates several processes in woody plants including vascular development, photosynthesis, and branching by affecting the proliferation and differentiation of parenchyma cells.