Analysis of polycomb repressive complex 2 (PRC2) subunits in Picea abies with a focus on embryo development.

BACKGROUND: Conserved polycomb repressive complex 2 (PRC2) mediates H3K27me3 to direct transcriptional repression and has a key role in cell fate determination and cell differentiation in both animals and plants. PRC2 subunits have undergone independent multiplication and functional divergence in higher plants. However, relevant information is still absent in gymnosperms. RESULTS: To launch gymnosperm PRC2 research, we identified and cloned the PRC2 core component genes in the conifer model species Picea abies, including one Esc/FIE homolog PaFIE, two p55/MSI homologs PaMSI1a and PaMSI1b, two E(z) homologs PaKMT6A2 and PaKMT6A4, a Su(z)12 homolog PaEMF2 and a PaEMF2-like fragment. Phylogenetic and protein domain analyses were conducted. The Esc/FIE homologs were highly conserved in the land plant, except the monocots. The other gymnospermous PRC2 subunits underwent independent evolution with angiospermous species to different extents. The relative transcript levels of these genes were measured in endosperm and zygotic and somatic embryos at different developmental stages. The obtained results proposed the involvement of PaMSI1b and PaKMT6A4 in embryogenesis and PaKMT6A2 and PaEMF2 in the transition from embryos to seedlings. The PaEMF2-like fragment was predominantly expressed in the endosperm but not in the embryo. In addition, immunohistochemistry assay showed that H3K27me3 deposits were generally enriched at meristem regions during seed development in P. abies. CONCLUSIONS: This study reports the first characterization of the PRC2 core component genes in the coniferous species P. abies. Our work may enable a deeper understanding of the cell reprogramming process during seed and embryo development and may guide further research on embryonic potential and development in conifers.

Integrated Lipidomic and Transcriptomic Analysis Reveals Phospholipid Changes in Somatic Embryos of Picea asperata in Response to Partial Desiccation.

Partial desiccation treatment (PDT) is an effective technology for promoting the germination and conversion of conifer somatic embryos (SEs). PDT, as a drought stress, induces intensive physiological responses in phospholipid metabolism, which are not well understood in the conifer SEs. Here, we integrated lipidomics, transcriptomics and proteomics analyses to reveal the molecular basis of lipid remodeling under PDT in Picea asperata SEs. Among the 82 lipid molecular species determined by mass spectrometry, phosphatidic acid (PA) had a significant effect after PDT and was the most critical lipid in the response to PDT. The transcriptomics results showed that multiple transcripts in the glycerolipid and glycerophospholipid metabolism pathways were differentially expressed, and these included five PLDα1 transcripts that catalyze the conversion of phosphatidylcholine (PC) to PA. Furthermore, the enzyme activity of this phospholipase D (PLD) was significantly enhanced in response to PDT, and PDT also significantly increased the protein level of PLDα1 (MA_10436582g0020). In addition, PA is a key factor in gibberellin, abscisic acid and ethylene signal transduction. One GDI1, one DELLA, three ABI1s, two SnRK2s, one CTR and 12 ERFs showed significantly differential expression between SEs before and after PDT in this study. Our data suggest that the observed increases in the PA contents might result from the activation of PLDα by PDT. PA not only affects the physical and chemical properties of the cell membrane but also participates in plant hormone signal transduction. Our work provides novel insight into the molecular mechanism through which PDT promotes the germination of SEs of coniferous tree species and fills the gap in the understanding of the mechanism of somatic embryo lipid remodeling in response to PDT.

Genome-wide analysis of lncRNA and mRNA expression and endogenous hormone regulation during tension wood formation in Catalpa bungei.

BACKGROUND: Phytohormones are the key factors regulating vascular development in plants, and they are also involved in tension wood (TW) formation. Although the theory of hormone distribution in TW formation is widely supported, the effects of endogenous hormones on TW formation have not yet been assessed. In this study, TW formation was induced in Catalpa bungei by artificial bending. The phytohormone content of TW, opposite wood (OW) and normal wood (NW) was determined using liquid chromatography-mass spectrometry (LC-MS), and transcriptome sequencing was performed. The hormone content and related gene expression data were comprehensively analyzed. RESULTS: The results of analyses of the plant hormone contents indicated significantly higher levels of cis-zeatin (cZ), indoleacetic acid (IAA) and abscisic acid (ABA) in TW than in OW. Genes involved in the IAA and ABA synthesis pathways, such as ALDH (evm. MODEL: group5.1511) and UGT (evm. MODEL: scaffold36.20), were significantly upregulated in TW. and the expression levels of ARF (evm. MODEL: group5.1332), A-ARR (evm. MODEL: group0.1600), and TCH4 (evm. MODEL: group2.745), which participate in IAA, cZ and Brassinolide (BR) signal transduction, were significantly increased in TW. In particular, ARF expression may be regulated by long noncoding RNAs (lncRNAs) and the HD-ZIP transcription factor ATHB-15. CONCLUSIONS: We constructed a multiple hormone-mediated network of C. bungei TW formation based on hormone levels and transcriptional expression profiles were identified during TW formation.

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.

Complete plastome sequences of Picea asperata and P. crassifolia and comparative analyses with P. abies and P. morrisonicola.

Picea asperata and P. crassifolia have sympatric ranges and are closely related, but the differences between these species at the plastome level are unknown. To better understand the patterns of variation among Picea plastomes, the complete plastomes of P. asperata and P. crassifolia were sequenced. Then, the plastomes were compared with the complete plastomes of P. abies and P. morrisonicola, which are closely and distantly related to the focal species, respectively. We also used these sequences to construct phylogenetic trees to determine the relationships among and between the four species as well as additional taxa from Pinaceae and other gymnosperms. Analysis of our sequencing data allowed us to identify 438 single nucleotide polymorphism (SNPs) point mutation events, 95 indel events, four inversion events, and seven highly variable regions, including six gene spacer regions (psbJ-petA, trnT-psaM, trnS-trnD, trnL-rps4, psaC-ccsA, and rps7-trnL) and one gene (ycf1). The highly variable regions are appropriate targets for future use in the phylogenetic reconstructions of closely related, sympatric species of Picea as well as Pinaceae in general.

Transcriptomic, Proteomic, and Metabolic Profiles of Catalpa bungei Tension Wood Reveal New Insight Into Lignin Biosynthesis Involving Transcription Factor Regulation.

Lignin is a complex polymer in plant cell walls whose proportion is second only to that of cellulose and plays an important role in the mechanical properties of wood and stress resistance of plants. Here, we induced tension wood (TW) formation in Catalpa bungei by artificial bending and analyzed the lignin metabolism of the TW. LC-MS analysis showed that a significantly higher content of coniferyl aldehyde was observed in the TW cell wall than in the opposite wood (OW) and normal wood (NW) cell walls. TW had significantly lower contents of coniferyl alcohol than OW and NW. Raman spectroscopy results indicated that TW had lower total lignin than OW and NW. The transcription and translation levels of most of the differentially expressed genes (DEGs) involved in lignin monomer biosynthesis indicated upregulation in TW/OW and TW/NW. We found no significant difference in the transcription levels of three collision gases (CADs) between TW and OW or between NW, but their translation levels were significantly downregulated in TW, suggesting post-transcriptional control for CAD. We predicted and analyzed transcription factors that could target DEGs involved in lignin monomer biosynthesis in TW. Based on the analysis of the relationships of targeting and coexpression, we found that NAC (evm.model.group1.695) could potentially target 4CLs and CCoAOMT, that HD-Zip (evm.model.group7.1157) had potential targeting relationships with CCoAOMT, F5H, and CCR, and that their expression levels were significantly positive. It is speculated that the upregulation of NAC and HD-ZIP transcription factors activates the expression of downstream target genes, which leads to a significant increase in coniferyl aldehyde in TW. However, the decrease in total lignin in TW may be caused by the significant downregulation of CAD translation and the significant decrease in precursors (coniferyl alcohol). Whether the expression of CAD genes is regulated by post-transcriptional control and affects TW lignin metabolism needs further study.

PICEAdatabase: a web database for Picea omics and phenotypic information.

Picea belongs to the Pinaceae family and is a famous commercial tree species because of its straight trunk and excellent timber traits. Recently, omics have been widely used for fundamental and mechanism studies on Picea plants. To improve the accessibility to omics and phenotypic data and facilitate further studies, we compiled the sequences of 2 chloroplast genomes (Picea crassifolia and Picea asperata) and 32 complete omics data sets, including 20 transcriptomes, 4 proteomes, 2 degradomes and 6 microRNAs from P. crassifolia, P. asperata, Picea balfouriana and Picea abies tissues under different treatments, in PICEAdatabase. In addition, phenotypic data on plant growth and wood property traits were collected from two field trials of P. crassifolia. PICEAdatabase also includes useful analysis tools, such as BLAST, DESeq2 and JBrowse, to assist with analyses.

ConTEdb: a comprehensive database of transposable elements in conifers.

Conifers are the largest and most ubiquitous group of gymnosperms and have significant ecological significance and economic importance. However, the huge and complex genomes have hindered the sequencing and mining of conifer genomes. In this study, we identified 413 423 transposable elements (TEs) from Picea abies, Picea glauca and Pinus taeda using a combination of multiple approaches and classified them into 11 133 families. A comprehensive web-based database, ConTEdb, was constructed and served for researchers. ConTEdb enables users to browse, retrieve and download the TE sequences from the database. Several analysis tools are integrated into ConTEdb to help users mine the TE data easily and effectively. In summary, ConTEdb provides a platform to study TE biology and functional genomics in conifers.