Comparative proteomic analysis of Populus trichocarpa early stem from primary to secondary growth

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.

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

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

Overexpression of miR390b promotes stem elongation and height growth in Populus

MicroRNA390 (miR390) is involved in plant growth and development by down-regulating the expression of the downstream genes trans-acting short interfering RNA3 (TAS3) and AUXIN RESPONSE FACTORs (ARFs). There is a scarcity of research on the involvement of the miR390-TAS3-ARFs pathway in the stem development of Populus. Here, differentially expressed miRNAs during poplar stem development were screened by small RNA sequencing analysis, and a novel function of miR390b in stem development was revealed. Overexpression of miR390b (OE-miR390b) resulted in a large increase in the number of xylem fiber cells and a slight decrease in the cell length at the longitudinal axis. Overall increases in stem elongation and plant height were observed in the OE-miR390b plants. According to transcriptome sequencing results and transient co-expression analysis, TAS3.1 and TAS3.2 were identified as the target genes of miR390 in poplar and were negatively regulated by miR390 in the apex. The transcription levels of ARF3.2 and ARF4 were significantly repressed in OE-miR390b plants and strongly negatively correlated with the number of xylem fiber cells along the longitudinal axis. These findings indicate that the conserved miR390-TAS3-ARFs pathway in poplar is involved in stem elongation and plant height growth.

Cas9/gRNA-Mediated Mutations in PtrFLA40 and PtrFLA45 Reveal Redundant Roles in Modulating Wood Cell Size and SCW Synthesis in Poplar

Fasciclin-like arabinogalactan proteins (FLAs) play an important role in plant development and adaptation to the environment. However, the roles of FLAs in wood formation remain poorly understood. Here, we identified a total of 50 PtrFLA genes in poplar. They were classified into four groups: A to D, among which group A was the largest group with 28 members clustered into four branches. Most PtrFLAs of group A were dominantly expressed in developing xylem based on microarray and RT-qPCR data. The roles of PtrFLA40 and PtrFLA45 in group A were investigated via the Cas9/gRNA-induced mutation lines. Loss of PtrFLA40 and PtrFLA45 increased stem length and diameter in ptrfla40ptrfla45 double mutants, but not in ptrfla40 or ptrfla45 single mutants. Further, our findings indicated that the ptrfla40ptrfla45 mutants enlarged the cell size of xylem fibers and vessels, suggesting a negative modulation in stem xylem cell size. In addition, wood lignin content in the ptrfla40fla45 mutants was increased by nearly 9%, and the lignin biosynthesis-related genes were significantly up-regulated in the ptrfla40fla45 mutants, in agreement with the increase in wood lignin content. Overall, Cas9/gRNA-mediated mutations in PtrFLA40 and PtrFLA45 reveal redundant roles in modulating wood cell size and secondary cell wall (SCW) synthesis in poplar.

Network-Based Analysis to Identify Hub Genes Involved in Spatial Root Response to Mechanical Constrains

Previous studies report that the asymmetric response, observed along the main poplar woody bent root axis, was strongly related to both the type of mechanical forces (compression or tension) and the intensity of force displacement. Despite a large number of targets that have been proposed to trigger this asymmetry, an understanding of the comprehensive and synergistic effect of the antistress spatially related pathways is still lacking. Recent progress in the bioinformatics area has the potential to fill these gaps through the use of in silico studies, able to investigate biological functions and pathway overlaps, and to identify promising targets in plant responses. Presently, for the first time, a comprehensive network-based analysis of proteomic signatures was used to identify functions and pivotal genes involved in the coordinated signalling pathways and molecular activities that asymmetrically modulate the response of different bent poplar root sectors and sides. To accomplish this aim, 66 candidate proteins, differentially represented across the poplar bent root sides and sectors, were grouped according to their abundance profile patterns and mapped, together with their first neighbours, on a high-confidence set of interactions from STRING to compose specific cluster-related subnetworks (I–VI). Successively, all subnetworks were explored by a functional gene set enrichment analysis to identify enriched gene ontology terms. Subnetworks were then analysed to identify the genes that are strongly interconnected with other genes (hub gene) and, thus, those that have a pivotal role in the bent root asymmetric response. The analysis revealed novel information regarding the response coordination, communication, and potential signalling pathways asymmetrically activated along the main root axis, delegated mainly to Ca2+ (for new lateral root formation) and ROS (for gravitropic response and lignin accumulation) signatures. Furthermore, some of the data indicate that the concave side of the bent sector, where the mechanical forces are most intense, communicates to the other (neighbour and distant) sectors, inducing spatially related strategies to ensure water uptake and accompanying cell modification. This information could be critical for understanding how plants maintain and improve their structural integrity—whenever and wherever it is necessary—in natural mechanical stress conditions.

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

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

Transcriptome-wide identification and characterization of microRNAs in diverse phases of wood formation in Populus trichocarpa

We applied miRNA expression profiling method to Populus trichocarpa stems of the three developmental stages, primary stem (PS), transitional stem (TS), and secondary stem (SS), to investigate miRNA species and their regulation on lignocellulosic synthesis and related processes. We obtained 892, 872, and 882 known miRNAs and 1727, 1723, and 1597 novel miRNAs, from PS, TS, and SS, respectively. Comparisons of these miRNA species among different developmental stages led to the identification of 114, 306, and 152 differentially expressed miRNAs (DE-miRNAs), which had 921, 2639, and 2042 candidate target genes (CTGs) in the three respective stages of the same order. Correlation analysis revealed 47, 439, and 71 DE-miRNA-CTG pairs of high negative correlation in PS, TS, and SS, respectively. Through biological process analysis, we finally identified 34, 6, and 76 miRNA-CTG pairs from PS, TS, and SS, respectively, and the miRNA target genes in these pairs regulate or participate lignocellulosic biosynthesis-related biological processes: cell division and differentiation, cell wall modification, secondary cell wall biosynthesis, lignification, and programmed cell death processes. This is the first report on an integrated analysis of genome-wide mRNA and miRNA profilings during multiple phases of poplar stem development. Our analysis results imply that individual miRNAs modulate secondary growth and lignocellulosic biosynthesis through regulating transcription factors and lignocellulosic biosynthetic pathway genes, resulting in more dynamic promotion, suppression, or regulatory circuits. This study advanced our understanding of many individual miRNAs and their essential, diversified roles in the dynamic regulation of secondary growth in woody tree species.

Functional understanding of secondary cell wall cellulose synthases in Populus trichocarpa via the Cas9/gRNA-induced gene knockouts

Plant cellulose is synthesized by a large plasma membrane-localized cellulose synthase (CesA) complex. However, an overall functional determination of secondary cell wall (SCW) CesAs is still lacking in trees, especially one based on gene knockouts. Here, the Cas9/gRNA-induced knockouts of PtrCesA4, 7A, 7B, 8A and 8B genes were produced in Populus trichocarpa. Based on anatomical, immunohistochemical and wood composition evidence, we gained a comprehensive understanding of five SCW PtrCesAs at the genetic level. Complete loss of PtrCesA4, 7A/B or 8A/B led to similar morphological abnormalities, indicating similar and nonredundant genetic functions. The absence of the gelatinous (G) layer, one-layer-walled fibres and a 90% decrease in cellulose in these mutant woods revealed that the three classes of SCW PtrCesAs are essential for multilayered SCW structure and wood G-fibre. In addition, the mutant primary and secondary phloem fibres lost the n(G + L)- and G-layers and retained the thicker S-layers (L, lignified; S, secondary). Together with polysaccharide immunolocalization data, these findings suggest differences in the role of SCW PtrCesAs-synthesized cellulose in wood and phloem fibre wall structures. Overall, this functional understanding of the SCW PtrCesAs provides further insights into the impact of lacking cellulose biosynthesis on growth, SCW, wood G-fibre and phloem fibre wall structures in the tree.

Molecular Changes Concomitant with Vascular System Development in Mature Galls Induced by Root-Knot Nematodes in the Model Tree Host Populus tremula × P. alba

One of the most striking features occurring in the root-knot nematode Meloidogyne incognita induced galls is the reorganization of the vascular tissues. During the interaction of the model tree species Populus and M. incognita, a pronounced xylem proliferation was previously described in mature galls. To better characterise changes in expression of genes possibly involved in the induction and the formation of the de novo developed vascular tissues occurring in poplar galls, a comparative transcript profiling of 21-day-old galls versus uninfected root of poplar was performed. Genes coding for transcription factors associated with procambium maintenance and vascular differentiation were shown to be differentially regulated, together with genes partaking in phytohormones biosynthesis and signalling. Specific signatures of transcripts associated to primary cell wall biosynthesis and remodelling, as well as secondary cell wall formation (cellulose, xylan and lignin) were revealed in the galls. Ultimately, we show that molecules derived from the monolignol and salicylic acid pathways and related to secondary cell wall deposition accumulate in mature galls.

A Molecular Blueprint of Lignin Repression

Although lignin is essential to ensure the correct growth and development of land plants, it may be an obstacle to the production of lignocellulosics-based biofuels, and reduces the nutritional quality of crops used for human consumption or livestock feed. The need to tailor the lignocellulosic biomass for more efficient biofuel production or for improved plant digestibility has fostered considerable advances in our understanding of the lignin biosynthetic pathway and its regulation. Most of the described regulators are transcriptional activators of lignin biosynthesis, but considerably less attention has been devoted to the repressors of this pathway. We provide a comprehensive overview of the molecular factors that negatively impact on the lignification process at both the transcriptional and post-transcriptional levels.

Data Integration in Poplar: 'Omics Layers and Integration Strategies

Populus trichocarpa is an important biofuel feedstock that has been the target of extensive research and is emerging as a model organism for plants, especially woody perennials. This research has generated several large ‘omics datasets. However, only few studies in Populus have attempted to integrate various data types. This review will summarize various ‘omics data layers, focusing on their application in Populus species. Subsequently, network and signal processing techniques for the integration and analysis of these data types will be discussed, with particular reference to examples in Populus.

Plant proteases during developmental programmed cell death

Proteases are among the key regulators of most forms of programmed cell death (PCD) in animals. Many PCD processes have also been associated with protease expression or activation in plants, However, functional evidence for the roles and actual modes of action of plant proteases in PCD remains surprisingly limited. In this review, we provide an update on protease involvement in the context of developmentally regulated plant PCD. To illustrate the diversity of protease functions, we focus on several prominent developmental PCD processes, including xylem and tapetum maturation, suspensor elimination, endosperm degradation, and seed coat formation, as well as plant senescence processes. Despite the substantial advances in the field, protease functions are often only correlatively linked to developmental PCD, and the specific molecular roles of proteases in many developmental PCD processes remain to be elucidated.

Genome-wide analysis of the Catalpa bungei caffeic acid O-methyltransferase (COMT) gene family: identification and expression profiles in normal, tension, and opposite wood

Caffeic acid O-methyltransferase (COMT) is an important protein that participates in lignin synthesis and is associated with the ratio of G-/S-type lignin in plants. COMTs are associated with the wood properties of forest trees; however, little known about the COMT family in Catalpa bungei, a valuable timber tree species in China . We performed a comprehensive analysis of COMT genes in the C. bungei genome by describing the gene structure and phylogenetic relationships of each family member using bioinformatics-based methods. A total of 23 putative COMT genes were identified using the conserved domain sequences and amino acid sequences of COMTs from Arabidopsis thaliana and Populus trichocarpa as probes. Phylogenetic analysis showed that 23 CbuCOMTs can be divided into three groups based on their structural characteristics; five conserved domains were found in the COMT family. Promoter analysis indicated that the CbuCOMT promoters included various cis-acting elements related to growth and development. Real-time quantitative polymerase chain reaction (PCR) analysis showed differential expression among CbuCOMTs. CbuCOMT2, 7, 8, 9, 10, 12, 13, 14, 21, and 23 were mainly expressed in xylem. Only CbuCOMT23 was significantly downregulated in tension wood and upregulated in opposite wood compared to normal wood. Our study provides new information about the CbuCOMT gene family and will facilitate functional characterisation in further research.

Exploring Cell Wall Composition and Modifications During the Development of the Gynoecium Medial Domain in Arabidopsis

In Arabidopsis, the gynoecium, the inner whorl of the flower, is the female reproductive part. Many tissues important for fertilization such as the stigma, style, transmitting tract, placenta, ovules, and septum, comprising the medial domain, arise from the carpel margin meristem. During gynoecium development, septum fusion occurs and tissues form continuously to prepare for a successful pollination and fertilization. During gynoecium development, cell wall modifications take place and one of the most important is the formation of the transmitting tract, having a great impact on reproductive competence because it facilitates pollen tube growth and movement through the ovary. In this study, using a combination of classical staining methods, fluorescent dyes, and indirect immunolocalization, we analyzed cell wall composition and modifications accompanying medial domain formation during gynoecium development. We detected coordinated changes in polysaccharide distribution through time, cell wall modifications preceding the formation of the transmitting tract, mucosubstances increase during transmitting tract formation, and a decrease of mannan distribution. Furthermore, we also detected changes in lipid distribution during septum fusion. Proper cell wall composition and modifications are important for postgenital fusion of the carpel (septum fusion) and transmitting tract formation, because these tissues affect plant reproductive competence.

Proteomics analysis reveals the molecular mechanism underlying the transition from primary to secondary growth of poplar

Wood is the most important natural source of energy and also provides fuel and fiber. Considering the significant role of wood, it is critical to understand how wood is formed. Integration of knowledge about wood development at the cellular and molecular levels will allow more comprehensive understanding of this complex process. In the present study, we used a comparative proteomic approach to investigate the differences in protein profiles between primary and secondary growth in young poplar stems using tandem mass tag (TMT)-labeling. More than 10,816 proteins were identified, and, among these, 3106 proteins were differentially expressed during primary to secondary growth. Proteomic data were validated using a combination of histochemical staining, enzyme activity assays, and quantitative real-time PCR. Bioinformatics analysis revealed that these differentially expressed proteins are related to various metabolic pathways, mainly including signaling, phytohormones, cell cycle, cell wall, secondary metabolism, carbohydrate and energy metabolism, and protein metabolism as well as redox and stress pathways. This large proteomics dataset will be valuable for uncovering the molecular changes occurring during the transition from primary to secondary growth. Further, it provides new and accurate information for tree breeding to modify wood properties.