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.

Long non-coding RNA-mediated competing endogenous RNA regulatory network during flower development and color formation in Melastoma candidum.

M. candidum, an evergreen shrubby flower known for its superior adaptation ability in South China, has gained increased attention in garden applications. However, scant attention has been paid to its flower development and color formation process at the non-coding RNA level. To fill this gap, we conducted a comprehensive analysis based on long non-coding RNA sequencing (lncRNA-seq), RNA-seq, small RNA sequencing (sRNA-seq), and widely targeted metabolome detection of three different flower developmental stages of M. candidum. After differentially expressed lncRNAs (DElncRNAs), differentially expressed mRNAs (DEmRNAs), differentially expressed microRNAs (DEmiRNAs), and differentially synthesized metabolites (DSmets) analyses between the different flower developmental stages, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were conducted to identify some key genes and metabolites in flavonoid, flavone, anthocyanin, carotenoid, and alkaloid-related GO terms and biosynthetic pathways. Three direct-acting models, including antisense-acting, cis-acting, and trans-acting between lncRNAs and mRNAs, were detected to illustrate the direct function of lncRNAs on target genes during flower development and color formation. Based on the competitive endogenous RNA (ceRNA) regulatory theory, we constructed a lncRNA-mediated regulatory network composed of DElncRNAs, DEmiRNAs, DEmRNAs, and DSmets to elucidate the indirect role of lncRNAs in the flower development and color formation of M. candidum. By utilizing correlation analyses between DERNAs and DSmets within the ceRNA regulatory network, alongside verification trials of the ceRNA regulatory mechanism, the study successfully illustrated the significance of lncRNAs in flower development and color formation process. This research provides a foundation for improving and regulating flower color at the lncRNA level in M. candidum, and sheds light on the potential applications of non-coding RNA in studies of flower development.

PeGSTU58, a Glutathione S-Transferase from Populus euphratica, Enhances Salt and Drought Stress Tolerance in Transgenic Arabidopsis.

Glutathione S-transferases (GSTs) play a crucial role in responding to abiotic stress and are an important target for research on plant stress tolerance mechanisms. Populus euphratica is a promising candidate species for investigating the abiotic tolerance mechanisms in woody plants. In our previous study, PeGSTU58 was identified as being associated with seed salinity tolerance. In the present study, PeGSTU58 was cloned from P. euphratica and functionally characterized. PeGSTU58 encodes a Tau class GST and is located in both the cytoplasm and nucleus. Transgenic Arabidopsis overexpressing PeGSTU58 displayed enhanced tolerance to salt and drought stress. Under salt and drought stress, the transgenic plants exhibited significantly higher activities of antioxidant enzymes, including SOD, POD, CAT, and GST, compared to the wild-type (WT) plants. Additionally, the expression levels of several stress-responsive genes, including DREB2A, COR47, RD22, CYP8D11, and SOD1 were upregulated in PeGSTU58 overexpression lines compared to those in WT Arabidopsis under salt and drought stress conditions. Furthermore, yeast one-hybrid assays and luciferase analysis showed that PebHLH35 can directly bind to the promoter region of PeGSTU58 and activate its expression. These results indicated that PeGSTU58 was involved in salt and drought stress tolerances by maintaining ROS homeostasis, and its expression was positively regulated by PebHLH35.

Pollen-specific gene SKU5-SIMILAR 13 enhances growth of pollen tubes in the transmitting tract in Arabidopsis.

Pollen tube growth and penetration in female tissues are essential for the transfer of sperm to the embryo sac during plant pollination. Despite its importance during pollination, little is known about the mechanisms that mediate pollen tube growth in female tissues. In this study, we identified an Arabidopsis thaliana pollen/pollen tube-specific gene, SKU5-SIMILAR 13 (SKS13), which was critical for the growth of pollen tubes in the transmitting tract. The SKS13 protein was distributed throughout the cytoplasm and pollen tube walls at the apical region. In comparison with wild-type pollen tubes, those of the sks13 mutants burst more frequently when grown in vitro. Additionally, the growth of sks13 pollen tubes was retarded in the transmitting tract, thereby resulting in decreased male fertility. The accumulation of pectin and cellulose in the cell wall of sks13 pollen tubes was altered, and the content of jasmonic acid (JA) in sks13 pollen was reduced. The pollen tubes treated with an inhibitor of JA biosynthesis grew much more slowly and had an altered distribution of pectin, which is similar to the pollen tube phenotypes of the SKS13 mutation. Our results suggest that SKS13 is essential for pollen tube growth in the transmitting tract by mediating the biosynthesis of JA that modifies the components of pollen tube cell walls.

Corrigendum: BEL1-like Homeodomain Protein BLH6a Is a Negative Regulator of CAld5H2 in Sinapyl Alcohol Monolignol Biosynthesis in Poplar.

[This corrects the article DOI: 10.3389/fpls.2021.695223.].

Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus.

High-throughput single-cell RNA sequencing (scRNA-seq) has advantages over traditional RNA-seq to explore spatiotemporal information on gene dynamic expressions in heterogenous tissues. We performed Drop-seq, a method for the dropwise sequestration of single cells for sequencing, on protoplasts from the differentiating xylem of Populus alba × Populus glandulosa. The scRNA-seq profiled 9,798 cells, which were grouped into 12 clusters. Through characterization of differentially expressed genes in each cluster and RNA in situ hybridizations, we identified vessel cells, fiber cells, ray parenchyma cells and xylem precursor cells. Diffusion pseudotime analyses revealed the differentiating trajectory of vessels, fiber cells and ray parenchyma cells and indicated a different differentiation process between vessels and fiber cells, and a similar differentiation process between fiber cells and ray parenchyma cells. We identified marker genes for each cell type (cluster) and key candidate regulators during developmental stages of xylem cell differentiation. Our study generates a high-resolution expression atlas of wood formation at the single cell level and provides valuable information on wood formation.

Arabidopsis AtPRP17 functions in embryo development by regulating embryonic patterning.

MAIN CONCLUSION: Arabidopsis AtPRP17, a homolog of yeast splicing factor gene PRP17, is expressed in siliques and embryos and functions in embryo development via regulating embryonic patterning. Yeast splicing factor PRP17/CDC40 is essential for cell growth through involvement in cell cycle regulation. Arabidopsis genome encodes a homolog of PRP17, AtPRP17; however, its function in Arabidopsis development is unknown. This study showed that AtPRP17 was highly expressed in siliques and embryos, and the protein was localized in the nucleus. The loss-of-function mutation of AtPRP17 led to shrunken seeds in Arabidopsis mature siliques. Further analysis revealed that the defective mature seeds of the mutant resulted from abnormal embryos with shriveled cotyledons, unequal cotyledons, swollen and shortened hypocotyls, or shortened radicles. During embryogenesis, mutant embryos showed delayed development and defective patterning of the apical and base domains, such as inhibited cotyledons and disorganized quiescent center cells and columella. Our results suggested that AtPRP17 functions in Arabidopsis embryo development via regulating embryonic patterning.

BEL1-like Homeodomain Protein BLH6a Is a Negative Regulator of CAl5H2 in Sinapyl Alcohol Monolignol Biosynthesis in Poplar.

Lignin is one of the major components of xylem cell walls in tree stems. The lignin in the wood of most flowering plants (dicotyledonous angiosperms) is typically polymerized from three monolignol precursors, coniferyl alcohol, sinapyl alcohol, and p-coumaroyl alcohol, resulting in guaiacyl (G), syringyl (S), and hydroxyphenyl (H) subunits, respectively. In this study, we focus on the transcriptional regulation of a coniferaldehyde 5-hydroxylase (CAld5H2) gene, which encodes a key enzyme for sinapyl alcohol biosynthesis. We carried out a yeast one-hybrid (Y1H) screen to identify candidate upstream transcription factors (TFs) regulating CAld5H2. We obtained 12 upstream TFs as potential regulators of CAld5H2. One of these TF genes, BLH6a, encodes a BEL1-like homeodomain (BLH) protein and negatively regulated the CAld5H2 promoter activity. The direct regulation of CAld5H2 promoter by BLH6a was supported by chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) and dominant repression of BLH6a in transgenic plants. Luciferase complementation imaging analyses showed extensive protein-protein interactions among these 12 TFs. We propose that BLH6a is a negative regulator of CAld5H2, which acts through combinatorial regulation of multiple TFs for sinapyl alcohol (S monolignol) biosynthesis in poplar.

Investigation Into Different Wood Formation Mechanisms Between Angiosperm and Gymnosperm Tree Species at the Transcriptional and Post-transcriptional Level.

Enormous distinctions of the stem structure and cell types between gymnosperms and angiosperms tree species are expected to cause quite different wood physical and mechanical attributes, however, the molecular mechanisms underlying the differing wood morphology are still unclear. In this study, we compared the transcriptomes obtained by RNA-Seq between Populus alba × P. glandulosa clone 84K, and Larix kaempferi (Lamb.) Carr trees. Available genome resource served as reference for P. alba × P. glandulosa and the Iso-Seq results of a three-tissues mixture (xylem, phloem, and leaf) were used as the reference for L. kaempferi to compare the xylem-specifically expressed genes and their alternative splicing model. Through screening, we obtained 13,907 xylem-specifically expressed genes (5,954 up-regulated, 7,953 down-regulated) in the xylem of P. alba × P. glandulosa, and 2,596 xylem-specifically expressed genes (1,648 up-regulated, 948 down-regulated) in the xylem of L. kaempferi. From the GO and KEGG analyses, some genes associated with two wood formation-related pathways, namely those for phenylpropanoid biosynthesis, and starch and sucrose metabolism, were successfully screened. Then the distributions and gene expression models between P. alba × P. glandulosa and L. kaempferi in those pathways were compared, which suggested differential wood formation processes between the angiosperm and gymnosperm trees. Furthermore, a Weight Gene Co-expression Network Analysis (WGCNA) for total xylem-specifically expressed genes in two species was conducted, from which wood formation-related modules were selected to build a co-expression network for the two tree species. The genes within this co-expression network showed different co-expression relationships between the angiosperm and gymnosperm woody species. Comparing the alternative splicing events for wood formation-related genes suggests a different post-transcriptional regulation process exists between the angiosperm and gymnosperm trees. Our research thus provides the foundation for the in-depth investigation of different wood formation mechanisms of angiosperm and gymnosperm species.

MicroRNA comparison between poplar and larch provides insight into the different mechanism of wood formation.

KEY MESSAGE: MiRNA transcriptome analysis of different tissues in poplar and larch suggests variant roles of miRNAs in regulating wood formation between two kinds of phyla. Poplar and larch belong to two different phyla. Both are ecological woody species and major resources for wood-related industrial applications. However, wood properties are different between these two species and the molecular basis is largely unknown. In this study, we performed high-throughput sequencing of microRNAs (miRNAs) in the three tissues, xylem, phloem and leaf of Populus alba × Populus glandulosa and Larix kaempferi. Differentially expressed miRNA (DEmiRNA) analysis identified 85 xylem-specific miRNAs in P. alba × P. glandulosa and 158 xylem-specific miRNAs in L. kaempferi. Among 36 common miRNAs, 12 were conserved between the two species. GO and KEGG analyses of the miRNA target genes showed similar metabolism in two species. Through KEGG and BLASTN, we predicted target genes of xylem differentially expressed (DEmiRNA) in the wood formation-related pathways and located DEmiRNAs in these pathways. A network was built for wood formation-related DEmiRNAs, their target genes and orthologous genes in Arabidopsis thaliana. Comparison of DEmiRNA and target gene annotation between P. alba × P. glandulosa and L. kaempferi suggested the different functions of DEmiRNAs and divergent mechanism in wood formation between two species, providing knowledge to understand wood formation mechanism in gymnosperm and angiosperm woody plants.

Involvement of CesA4, CesA7-A/B and CesA8-A/B in secondary wall formation in Populus trichocarpa wood.

Cellulose synthase A genes (CesAs) are responsible for cellulose biosynthesis in plant cell walls. In this study, functions of secondary wall cellulose synthases PtrCesA4, PtrCesA7-A/B and PtrCesA8-A/B were characterized during wood formation in Populus trichocarpa (Torr. & Gray). CesA RNAi knockdown transgenic plants exhibited stunted growth, narrow leaves, early necrosis, reduced stature, collapsed vessels, thinner fiber cell walls and extended fiber lumen diameters. In the RNAi knockdown transgenics, stems exhibited reduced mechanical strength, with reduced modulus of rupture (MOR) and modulus of elasticity (MOE). The reduced mechanical strength may be due to thinner fiber cell walls. Vessels in the xylem of the transgenics were collapsed, indicating that water transport in xylem may be affected and thus causing early necrosis in leaves. A dramatic decrease in cellulose content was observed in the RNAi knockdown transgenics. Compared with wildtype, the cellulose content was significantly decreased in the PtrCesA4, PtrCesA7 and PtrCesA8 RNAi knockdown transgenics. As a result, lignin and xylem contents were proportionally increased. The wood composition changes were confirmed by solid-state NMR, two-dimensional solution-state NMR and sum-frequency-generation vibration (SFG) analyses. Both solid-state nuclear magnetic resonance (NMR) and SFG analyses demonstrated that knockdown of PtrCesAs did not affect cellulose crystallinity index. Our results provided the evidence for the involvement of PtrCesA4, PtrCesA7-A/B and PtrCesA8-A/B in secondary cell wall formation in wood and demonstrated the pleiotropic effects of their perturbations on wood formation.

Characterization of Cellulose synthase-like D (CSLD) family revealed the involvement of PtrCslD5 in root hair formation in Populus trichocarpa.

Cellulose synthase-like D (CSLD) family was characterized for their expression and functions in Populus trichocarpa. Ten members, PtrCslD1-10, were identified in the P. trichocarpa genome, and they belong to 4 clades by phylogenetic tree analysis. qRT-PCR and promoter:GUS assays in Arabidopsis and P. trichocarpa displayed divergent expression patterns of these 10 PtrCSLD genes in root hairs, root tips, leaves, vascular tissues, xylem and flowers. Among PtrCslD2, PtrCslD4, PtrCslD5, PtrCslD6, and PtrCslD8 that all exhibited expression in root hairs, only PtrCslD5 could restore the root hairless phenotype of the atcsld3 mutant, demonstrating that PtrCslD5 is the functional ortholog of AtCslD3 for root hair formation. Our results suggest more possible functions for other PtrCslD genes in poplar.