Comprehensive transcriptome analysis of developing xylem responding to artificial bending and gravitational stimuli in Betula platyphylla

Genome-wide identification, expression and salt stress tolerance analysis of the GRAS transcription factor family in Betula platyphylla

The GRAS gene family is a plant-specific family of transcription factors and play a vital role in many plant growth processes and abiotic stress responses. Nevertheless, the functions of the GRAS gene family in woody plants, especially in Betula platyphylla (birch), are hardly known. In this study, we performed a genome-wide analysis of 40 BpGRAS genes (BpGRASs) and identified typical GRAS domains of most BpGRASs. The BpGRASs were unevenly distributed on 14 chromosomes of birch and the phylogenetic analysis of six species facilitated the clustering of 265 GRAS proteins into 17 subfamilies. We observed that closely related GRAS homologs had similar conserved motifs according to motif analysis. Besides, an analysis of the expression patterns of 26 BpGRASs showed that most BpGRASs were highly expressed in the leaves and responded to salt stress. Six BpGRASs were selected for cis-acting element analysis because of their significant upregulation under salt treatment, indicating that many elements were involved in the response to abiotic stress. This result further confirmed that these BpGRASs might participate in response to abiotic stress. Transiently transfected birch plants with transiently overexpressed 6 BpGRASs and RNAi-silenced 6 BpGRASs were generated for gain- and loss-of-function analysis, respectively. In addition, overexpression of BpGRAS34 showed phenotype resistant to salt stress, decreased the cell death and enhanced the reactive oxygen species (ROS) scavenging capabilities and proline content under salt treatment, consistent with the results in transiently transformed birch plants. This study is a systematic analysis of the GRAS gene family in birch plants, and the results provide insight into the molecular mechanism of the GRAS gene family responding to abiotic stress in birch plants.

Bouncing back stronger: Diversity, structure, and molecular regulation of gelatinous fiber development

Gelatinous fibers (G-fibers) are specialized contractile cells found in a diversity of vascular plant tissues, where they provide mechanical support and/or facilitate plant mobility. G-fibers are distinct from typical fibers by the presence of an innermost thickened G-layer, comprised mainly of axially oriented cellulose microfibrils. Despite the disparate developmental origins-tension wood fibers from the vascular cambium or primary phloem fibers from the procambium-G-fiber development, composition, and molecular signatures are remarkably similar; however, important distinctions do exist. Here, we synthesize current knowledge of the phylogenetic diversity, compositional makeup, and the molecular profiles that characterize G-fiber development and highlight open questions for future investigation.

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.

Comprehensive Transcriptome Analysis of Stem-Differentiating Xylem Upon Compression Stress in Cunninghamia Lanceolata

Compression wood (CW) in gymnosperm brings great difficulties to wood industry using wood as raw materials since CW presents special wood structure and have different physical and chemical properties from those of normal wood (NW). Chinese fir (Cunninghamia lanceolata) is widely distributed in China. However, global transcriptome profiling of coding and long non-coding RNA in response to compression stress has not been reported in the gymnosperm species. In this study, we revealed that CW in Chinese fir exhibited distinct morphology and cytology properties compared with those of NW, including high lignin content, thick and round tracheid cells. Furthermore, we combined both PacBio long-read SMRT sequencing (Iso-Seq) and Illumina short-read RNA-Seq to reveal the transcriptome in stem-differentiating xylem (SDX) under different time points (2, 26, and 74 h) upon compression stress in NW, CW, and OW (opposite wood), respectively. Iso-Seq was successfully assembled into 41,253 de-novo full-length transcriptome reference (average length 2,245 bp). Moreover, there were striking differences in expression upon compression stress, which were involved 13 and 7 key enzyme genes in the lignin and cellulose synthesis, respectively. Especially, we revealed 11 secondary growth-related transcription factors show differential expression under compression stress, which was further validated by qRT-PCR. Finally, the correlation between 6,533 differentially expressed coding genes and 372 differentially expressed long non-coding RNAs (lncRNAs) indicates that these lncRNAs may affect cell wall biogenesis and xyloglucan metabolism. In conclusion, our results provided comprehensive cytology properties and full-length transcriptome profiling of wood species upon compression stress. Especially we explored candidate genes, including both coding and long non-coding genes, and provided a theoretical basis for further research on the formation mechanism of CW in gymnosperm Chinese fir.

Transcriptional reprogramming of xylem cell wall biosynthesis in tension wood

Tension wood (TW) is a specialized xylem tissue developed under mechanical/tension stress in angiosperm trees. TW development involves transregulation of secondary cell wall genes, which leads to altered wood properties for stress adaptation. We induced TW in the stems of black cottonwood (Populus trichocarpa, Nisqually-1) and identified two significantly repressed transcription factor (TF) genes: class B3 heat-shock TF (HSFB3-1) and MYB092. Transcriptomic analysis and chromatin immunoprecipitation (ChIP) were used to identify direct TF-DNA interactions in P. trichocarpa xylem protoplasts overexpressing the TFs. This analysis established a transcriptional regulatory network in which PtrHSFB3-1 and PtrMYB092 directly activate 8 and 11 monolignol genes, respectively. The TF-DNA interactions were verified for their specificity and transactivator roles in 35 independent CRISPR-based biallelic mutants and overexpression transgenic lines of PtrHSFB3-1 and PtrMYB092 in P. trichocarpa. The gene-edited trees (mimicking the repressed PtrHSFB3-1 and PtrMYB092 under tension stress) have stem wood composition resembling that of TW during normal growth and under tension stress (i.e., low lignin and high cellulose), whereas the overexpressors showed an opposite effect (high lignin and low cellulose). Individual overexpression of the TFs impeded lignin reduction under tension stress and restored high levels of lignin biosynthesis in the TW. This study offers biological insights to further uncover how metabolism, growth, and stress adaptation are coordinately regulated in trees.

Transcriptomics and Proteomics Reveal the Cellulose and Pectin Metabolic Processes in the Tension Wood (Non-G-Layer) of Catalpa bungei

: Catalpa bungei is an economically important tree with high-quality wood and highly valuable to the study of wood formation. In this work, the xylem microstructure of C. bungei tension wood (TW) was observed, and we performed transcriptomics, proteomics and Raman spectroscopy of TW, opposite wood (OW) and normal wood (NW). The results showed that there was no obvious gelatinous layer (G-layer) in the TW of C. bungei and that the secondary wall deposition in the TW was reduced compared with that in the OW and NW. We found that most of the differentially expressed mRNAs and proteins were involved in carbohydrate polysaccharide synthesis. Raman spectroscopy results indicated that the cellulose and pectin content and pectin methylation in the TW were lower than those in the OW and NW, and many genes and proteins involved in the metabolic pathways of cellulose and pectin, such as galacturonosyltransferase (GAUT), polygalacturonase (PG), endoglucanase (CLE) and β-glucosidase (BGLU) genes, were significantly upregulated in TW. In addition, we found that the MYB2 transcription factor may regulate the pectin degradation genes PG1 and PG3, and ARF, ERF, SBP and MYB1 may be the key transcription factors regulating the synthesis and decomposition of cellulose. In contrast to previous studies on TW with a G-layer, our results revealed a change in metabolism in TW without a G-layer, and we inferred that the change in the pectin type, esterification and cellulose characteristics in the TW of C. bungei may contribute to high tensile stress. These results will enrich the understanding of the mechanism of TW formation.

The BpMYB4 Transcription Factor From Betula platyphylla Contributes Toward Abiotic Stress Resistance and Secondary Cell Wall Biosynthesis

The MYB (v-myb avian myeloblastosis viral oncogene homolog) family is one of the largest transcription factor families in plants, and is widely involved in the regulation of plant metabolism. In this study, we show that a MYB4 transcription factor, BpMYB4, identified from birch (Betula platyphylla Suk.) and homologous to EgMYB1 from Eucalyptus robusta Smith and ZmMYB31 from Zea mays L. is involved in secondary cell wall synthesis. The expression level of BpMYB4 was higher in flowers relative to other tissues, and was induced by artificial bending and gravitational stimuli in developing xylem tissues. The expression of this gene was not enriched in the developing xylem during the active season, and showed higher transcript levels in xylem tissues around sprouting and near the dormant period. BpMYB4 also was induced express by abiotic stress. Functional analysis indicated that expression of BpMYB4 in transgenic Arabidopsis (Arabidopsis thaliana) plants could promote the growth of stems, and result in increased number of inflorescence stems and shoots. Anatomical observation of stem sections showed lower lignin deposition, and a chemical contents test also demonstrated increased cellulose and decreased lignin content in the transgenic plants. In addition, treatment with 100 mM NaCl and 200 mM mannitol resulted in the germination rate of the over-expressed lines being higher than that of the wild-type seeds. The proline content in transgenic plants was higher than that in WT, but MDA content was lower than that in WT. Further investigation in birch using transient transformation techniques indicated that overexpression of BpMYB4 could scavenge hydrogen peroxide and O₂ ^.- and reduce cell damage, compared with the wild-type plants. Therefore, we believe that BpMYB4 promotes stem development and cellulose biosynthesis as an inhibitor of lignin biosynthesis, and has a function in abiotic stress resistance.

Effects of exogenous 24-epibrassinolide and brassinazole on negative gravitropism and tension wood formation in hybrid poplar (Populus deltoids × Populus nigra)

Exogenous 24-epibrassinolide (BL) and brassinazole (BRZ) have regulatory roles in G-fiber cell wall development and secondary xylem cell wall carbohydrate biosynthesis during tension wood formation in hybrid poplar. Brassinosteroids (BRs) play important roles in regulating gravitropism and vasculature development. Here, we report the effect of brassinosteroids on negative gravitropism and G-fiber cell wall development of the stem in woody angiosperms. We applied exogenous 24-epibrassinolide (BL) or its biosynthesis inhibitor brassinazole (BRZ) to slanted hybrid poplar trees (Populus deltoids × Populus nigra) and measured the morphology of gravitropic stems, anatomy and chemistry of secondary cell wall. We furthermore analyzed the expression levels of auxin transport and cellulose biosynthetic genes after 24-epibrassinolide (BL) or brassinazole (BRZ) application. The BL-treated seedlings showed no negative gravitropism bending, whereas application of BRZ dramatically enhanced negative gravitropic bending. BL treatment stimulated secondary xylem fiber elongation and G-fiber formation on the upper side of stems but delayed G-fiber maturation. BRZ inhibited xylem fiber elongation but induced the production of more mature G-fibers on the upper side of stems. Wood chemistry analyses and immunolocalization demonstrated that BL and BRZ treatments increased the cellulose content and modified the deposition of cell wall carbohydrates including arabinose, galactose and rhamnose in the secondary xylem. The expression of cellulose biosynthetic genes, especially those related to cellulose microfibril deposition (PtFLA12 and PtCOBL4) was significantly upregulated in BL- and BRZ-treated TW stems compared with control stems. The significant differences of G-fibers development and negative gravitropism bending between 24-epibrassinolide (BL) and brassinazole (BRZ) application suggest that brassinosteroids are important for secondary xylem development during tension wood formation. Our findings provide potential insights into the mechanism by which BRs regulate G-fiber cell wall development to accomplish negative gravitropism in TW formation.

BpNAC012 Positively Regulates Abiotic Stress Responses and Secondary Wall Biosynthesis

NAC (NAM, ATAF1/2, and CUC2) transcription factors play important roles in plant biological processes and stress responses. Here, we characterized the functional roles of BpNAC012 in white birch (Betula platyphylla). We found that BpNAC012 serves as a transcriptional activator. Gain- and loss-of-function analyses revealed that the transcript level of BpNAC012 was positively associated with salt and osmotic stress tolerance. BpNAC012 activated the core sequence CGT[G/A] to induce the expression of abiotic stress-responsive downstream genes, including Δ-1-pyrroline-5-carboxylate synthetase, superoxide dismutase, and peroxidase, resulting in enhanced salt and osmotic stress tolerance in BpNAC012 overexpression transgenic birch lines. We also showed that BpNAC012 is expressed predominantly in mature stems and that RNA interference-induced suppression of BpNAC012 caused a drastic reduction in the secondary wall thickening of stem fibers. Overexpression of BpNAC012 activated the expression of secondary wall-associated downstream genes by directly binding to the secondary wall NAC-binding element sites, resulting in ectopic secondary wall deposition in the stem epidermis. Moreover, salt and osmotic stresses elicited higher expression levels of lignin biosynthetic genes and elevated lignin accumulation in BpNAC012 overexpression lines. These findings provide insight into the functions of NAC transcription factors.

RNA-Seq analysis of differential gene expression in Betula luminifera xylem during the early stages of tension wood formation

Background

Betula luminifera H. Winkler, which is widely distributed in southern China, is an economically important broadleaf tree species. However, little genomic information of B. luminifera is available, and little is known about the molecular mechanisms of wood formation in this species. Meanwhile, few efforts have focused on investigating the early transcriptional changes during tension wood formation in woody plants.

Results

A reference transcriptome dataset was first generated containing 45,700 Unigenes, and 35,135 (76.9%) Unigenes were annotated by a BLAST similarity search against four public databases. Then, based on an anatomical investigation, the global gene expression changes during the early stages of tension wood formation were analyzed. Gene expression profiling showed that a total of 13,273 Unigenes were differentially regulated during the early stages of tension wood formation. Most genes involved in cellulose and lignin biosynthesis were highlighted to reveal their biological importance in tension wood formation. In addition, the transcription levels of many genes involved in the auxin response pathway were significantly changed during the early stages of tension wood formation. Furthermore, 18 TFs co-expressed with key enzymes of cellulose synthesis were identified.

Conclusions

Our results revealed the transcriptional changes associated with TW formation and identified potential key genes in the regulation of this process. These results will help to dissect the molecular mechanism of wood formation and provide key candidate genes for marker-assisted selection in B. luminifera.

Protein expression in tension wood formation monitored at high tissue resolution in Populus

Tension wood (TW) is a specialized tissue with contractile properties that is formed by the vascular cambium in response to gravitational stimuli. We quantitatively analysed the proteomes of Populus tremula cambium and its xylem cell derivatives in stems forming normal wood (NW) and TW to reveal the mechanisms underlying TW formation. Phloem-, cambium-, and wood-forming tissues were sampled by tangential cryosectioning and pooled into nine independent samples. The proteomes of TW and NW samples were similar in the phloem and cambium samples, but diverged early during xylogenesis, demonstrating that reprogramming is an integral part of TW formation. For example, 14-3-3, reactive oxygen species, ribosomal and ATPase complex proteins were found to be up-regulated at early stages of xylem differentiation during TW formation. At later stages of xylem differentiation, proteins involved in the biosynthesis of cellulose and enzymes involved in the biosynthesis of rhamnogalacturonan-I, rhamnogalacturonan-II, arabinogalactan-II and fasciclin-like arabinogalactan proteins were up-regulated in TW. Surprisingly, two isoforms of exostosin family proteins with putative xylan xylosyl transferase function and several lignin biosynthesis proteins were also up-regulated, even though xylan and lignin are known to be less abundant in TW than in NW. These data provided new insight into the processes behind TW formation.

Transcriptional and computational study of expansins differentially expressed in response to inclination in radiata pine

Plants have the ability to reorient their vertical growth when exposed to inclination. This response can be as quick as 2 h in inclined young pine (Pinus radiata D. Don) seedlings, with over accumulation of lignin observed after 9 days s. Several studies have identified expansins involved in cell expansion among other developmental processes in plants. Six putative expansin genes were identified in cDNA libraries isolated from inclined pine stems. A differential transcript abundance was observed by qPCR analysis over a time course of inclination. Five genes changed their transcript accumulation in both stem sides in a spatial and temporal manner compared with non-inclined stem. To compare these expansin genes, and to suggest a possible mechanism of action at molecular level, the structures of the predicted proteins were built by comparative modeling methodology. An open groove on the surface of the proteins composed of conserved zresidues was observed. Using a cellulose polymer as ligand the protein-ligand interaction was evaluated, with the results showing differences in the protein-ligand interaction mode. Differences in the binding energy interaction can be explained by changes in some residues that generate differences in electrostatic surface in the open groove region, supporting the participation of six members of multifamily proteins in this specific process. The data suggests participation of different expansin proteins in the dissembling and remodeling of the complex cell wall matrix during the reorientation response to inclination.

Expression of the MYB transcription factor gene BplMYB46 affects abiotic stress tolerance and secondary cell wall deposition in Betula platyphylla

Plant MYB transcription factors control diverse biological processes, such as differentiation, development and abiotic stress responses. In this study, we characterized BplMYB46, an MYB gene from Betula platyphylla (birch) that is involved in both abiotic stress tolerance and secondary wall biosynthesis. BplMYB46 can act as a transcriptional activator in yeast and tobacco. We generated transgenic birch plants with overexpressing or silencing of BplMYB46 and subjected them to gain- or loss-of-function analysis. The results suggest that BplMYB46 improves salt and osmotic tolerance by affecting the expression of genes including SOD, POD and P5CS to increase both reactive oxygen species scavenging and proline levels. In addition, BplMYB46 appears to be involved in controlling stomatal aperture to reduce water loss. Overexpression of BplMYB46 increases lignin deposition, secondary cell wall thickness and the expression of genes in secondary cell wall formation. Further analysis indicated that BplMYB46 binds to MYBCORE and AC-box motifs and may directly activate the expression of genes involved in abiotic stress responses and secondary cell wall biosynthesis whose promoters contain these motifs. The transgenic BplMYB46-overexpressing birch plants, which have improved salt and osmotic stress tolerance, higher lignin and cellulose content and lower hemicellulose content than the control, have potential applications in the forestry industry.

An Ethylene-responsive Factor BpERF11 Negatively Modulates Salt and Osmotic Tolerance in Betula platyphylla

Ethylene responsive factors (ERFs) play important roles in the abiotic stress; however, only a few ERF genes from woody plants have been functionally characterized. In the present study, an ERF gene from Betula platyphylla (birch), BpERF11, was functionally characterized in response to abiotic stress. BpERF11 is a nuclear protein, which could specifically bind to GCC boxes and DRE motifs. BpERF11-overexpressing and BpERF11 RNA interference (RNAi) knockdown plants were generated for gain- and loss-of-function analysis. BpERF11 negatively regulates resistance to salt and severe osmotic stress, and the transgenic birch plants overexpressing BpERF11 shows increased electrolyte leakage and malondialdehyde (MDA) contents. BpERF11 inhibits the expression of an AtMYB61 homologous gene, resulting in increased stomatal aperture, which elevated the transpiration rate. Furthermore, BpERF11 downregulates the expression of P5CS, SOD and POD genes, but upregulates the expression of PRODH and P5CDH, which results in reduced proline levels and increased reactive oxygen species (ROS) accumulation. BpERF11 also significantly inhibits the expression of LEA and dehydrin genes that involve in abiotic stress tolerance. Therefore, BpERF11 serves as a transcription factor that negatively regulates salt and severe osmotic tolerance by modulating various physiological processes.

Exogenous GA₃ Application Enhances Xylem Development and Induces the Expression of Secondary Wall Biosynthesis Related Genes in Betula platyphylla

Gibberellin (GA) is a key signal molecule inducing differentiation of tracheary elements, fibers, and xylogenesis. However the molecular mechanisms underlying the effect of GA on xylem elongation and secondary wall development in tree species remain to be determined. In this study, Betula platyphylla (birch) seeds were treated with 300 ppm GA₃ and/or 300 ppm paclobutrazol (PAC), seed germination was recorded, and transverse sections of hypocotyls were stained with toluidine blue; the two-month-old seedlings were treated with 50 μM GA₃ and/or 50 μM PAC, transverse sections of seedling stems were stained using phloroglucinol–HCl, and secondary wall biosynthesis related genes expression was analyzed by real-time quantitative PCR. Results indicated that germination percentage, energy and time of seeds, hypocotyl height and seedling fresh weight were enhanced by GA₃, and reduced by PAC; the xylem development was wider in GA₃-treated plants than in the control; the expression of NAC and MYB transcription factors, CESA, PAL, and GA oxidase was up-regulated during GA₃ treatment, suggesting their role in GA₃-induced xylem development in the birch. Our results suggest that GA₃ induces the expression of secondary wall biosynthesis related genes to trigger xylogenesis in the birch plants.

Quantitative Proteomic and Phosphoproteomic Approaches for Deciphering the Signaling Pathway for Tension Wood Formation in Poplar

Trees adjust their growth following forced changes in orientation to re-establish a vertical position. In angiosperms, this adjustment involves the differential regulation of vascular cambial activity between the lower (opposite wood) and upper (tension wood) sides of the leaning stem. We investigated the molecular mechanisms leading to the formation of differential wood types through a quantitative proteomic and phosphoproteomic analysis on poplar subjected to a gravitropic stimulus. We identified and quantified 675 phosphopeptides, corresponding to 468 phosphoproteins, and 3 763 nonphosphorylated peptides, corresponding to 1 155 proteins, in the differentiating xylem of straight-growing trees (control) and trees subjected to a gravitational stimulus during 8 weeks. About 1% of the peptides were specific to a wood type (straight, opposite, or tension wood). Proteins quantified in more than one type of wood were more numerous: a mixed linear model showed 389 phosphopeptides and 556 proteins to differ in abundance between tension wood and opposite wood. Twenty-one percent of the phosphoproteins identified here were described in their phosphorylated form for the first time. Our analyses revealed remarkable developmental molecular plasticity, with wood type-specific phosphorylation events, and highlighted the involvement of different proteins in the biosynthesis of cell wall components during the formation of the three types of wood.