Genome-wide investigation of DNA methylation dynamics reveals a critical role of DNA demethylation during the early somatic embryogenesis of Dimocarpus longan Lour

Transcriptomic, epigenomic and physiological comparisons reveal key factors for different manganese tolerances in three Chenopodium ambrosioides L. populations

Chenopodium ambrosioides is a manganese (Mn) hyperaccumulator that can be used for Mn-polluted soil phytoremediation. However, the mechanism of Mn tolerance of C. ambrosioides remains largely unknown. In this study, the key factors for Mn tolerance of C. ambrosioides was investigated from the aspects of DNA methylation pattern, gene expression regulation and physiological function. We found that the two genotypes of C. ambrosioides populations have differentiated tolerance to Mn stress (Mn-tolerant: CS and XC, Mn-sensitive: WH). Although there was no difference in Mn accumulation between two types under excess Mn, the biomass and photosynthetic systems were more severely inhibited in Mn-sensitive type, as well as suffering more serious oxidative damage. More differentially expressed genes (DEGs) were downregulated in the Mn-tolerant type, indicating that the Mn-tolerant type tends to inhibit gene expression to cope with Mn stress. DEGs related to metal transport, antioxidant system, phytohormone and transcription factors contribute to the tolerance of C. ambrosioides to Mn, and account for difference in Mn stress sensitivities between the Mn-sensitive and tolerant types. We also found that DNA methylation variation may help to cope with Mn stress. The global DNA methylation level in C. ambrosioides increased under Mn stress, especially in the Mn-sensitive type. Dozens of methylated loci were significantly associated with the Mn accumulation trait of C. ambrosioides, and some critical DEGs were regulated by DNA methylation. Our study comprehensively demonstrated the Mn tolerance mechanism of C. ambrosioides for the first time, and highlighted the roles of epigenetic modification in C. ambrosioides response to Mn stress. Our findings may contribute to elucidating the adaptation mechanism of hyperaccumulator to the heavy metal toxicity.

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.

Genome-wide circular RNA profiling and competing endogenous RNA regulatory network analysis provide new insights into the molecular mechanisms underlying early somatic embryogenesis in Dimocarpus longan Lour

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the 'non-homologous end-joining' (NHEJ) and 'butanoate metabolism' pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways 'mitogen-activated protein kinase signaling' and 'nitrogen metabolism'. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Triticale doubled haploid plant regeneration factors linked by structural equation modeling

Triticale regeneration via anther culture faces many difficulties, e.g., a low percentage of regenerated plants and the presence of albinos. Plant regeneration may be affected by abiotic stresses and by ingredients added to the induction medium. The latter influences biochemical pathways and plant regeneration efficiency. Among such ingredients, copper and silver ions acting as cofactors for enzymatic reactions are of interest. However, their role in plant tissue cultures and relationships with biochemical pathways has not been studied yet.

The study evaluated relationships between DNA methylation, changes in DNA sequence variation, and green plant regeneration efficiency influenced by copper and silver ions during triticale plant regeneration. For this purpose, a biological model based on donor plants and their regenerants, a methylation-sensitive amplified fragment length polymorphism, and structural equation modeling were employed.

The green plant regeneration efficiency varied from 0.71 to 6.06 green plants per 100 plated anthers. The values for the components of tissue culture-induced variation related to cytosine methylation in a CHH sequence context (where H is A, C, or T) were 8.65% for sequence variation, 0.76% for DNA demethylation, and 0.58% for de novo methylation. The proposed model states that copper ions affect the regeneration efficiency through cytosine methylation and may induce mutations through, e.g., oxidative processes, which may interfere with the green plant regeneration efficiency. The linear regression confirms that the plant regeneration efficiency rises with increasing copper ion concentration in the absence of Ag ions in the induction medium. The least absolute shrinkage and selection operator regression shows that de novo methylation, demethylation, and copper ions may be involved in the green plant regeneration efficiency. According to structural equation modeling, copper ions play a central role in the model determining the regeneration efficiency.

Weighted Gene Correlation Network Analysis (WGCNA) of Arabidopsis Somatic Embryogenesis (SE) and Identification of Key Gene Modules to Uncover SE-Associated Hub Genes

Somatic embryogenesis (SE), which occurs naturally in many plant species, serves as a model to elucidate cellular and molecular mechanisms of embryo patterning in plants. Decoding the regulatory landscape of SE is essential for its further application. Hence, the present study was aimed at employing Weighted Gene Correlation Network Analysis (WGCNA) to construct a gene coexpression network (GCN) for Arabidopsis SE and then identifying highly correlated gene modules to uncover the hub genes associated with SE that may serve as potential molecular targets. A total of 17,059 genes were filtered from a microarray dataset comprising four stages of SE, i.e., stage I (zygotic embryos), stage II (proliferating tissues at 7 days of induction), stage III (proliferating tissues at 14 days of induction), and stage IV (mature somatic embryos). This included 1,711 transcription factors and 445 EMBRYO DEFECTIVE genes. GCN analysis identified a total of 26 gene modules with the module size ranging from 35 to 3,418 genes using a dynamic cut tree algorithm. The module-trait analysis revealed that four, four, seven, and four modules were associated with stages I, II, III, and IV, respectively. Further, we identified a total of 260 hub genes based on the degree of intramodular connectivity. Validation of the hub genes using publicly available expression datasets demonstrated that at least 78 hub genes are potentially associated with embryogenesis; of these, many genes remain functionally uncharacterized thus far. In silico promoter analysis of these genes revealed the presence of cis-acting regulatory elements, "soybean embryo factor 4 (SEF4) binding site," and "E-box" of the napA storage-protein gene of Brassica napus; this suggests that these genes may play important roles in plant embryo development. The present study successfully applied WGCNA to construct a GCN for SE in Arabidopsis and identified hub genes involved in the development of somatic embryos. These hub genes could be used as molecular targets to further elucidate the molecular mechanisms underlying SE in plants.

Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.

Genome-wide identification and characterization of DEAD-box helicase family associated with early somatic embryogenesis in Dimocarpus longan Lour

DEAD-box (DDX) proteins belong to the largest subfamily of RNA helicase SF2, which contributes to all biological processes of RNA metabolism in the plant kingdom. Till now, no significant data are available regarding studies on DDX in Somatic Embryogenesis (SE) of woody plants. It is important to investigate the biological function of the DlDDX family in longan SE. Thus, a comprehensive analysis of 58 longan DEAD-box (DlDDX) genes characterization was performed by genome-wide identification and transcript abundance validation analysis. Homologous evolution has revealed that some DlDDXs in longan had high sequence similarity with Mus musculus, Citrus and Saccharomyces cerevisiae, indicating that DlDDXs were highly conservative in the animal, plant, and microorganism. Remarkably, gene duplication, purifying selection, and alternative splicing events, and new auxiliary domains have likely contributed to the functional evolution of DlDDX, indicating that DlDDX appeared neofunctionalization in longan. Besides, DlDDX3, 15, 28, 36 might interact with protein complex (MAC3A, MAC3B, CDC5, CBP20) of miRNA biosynthesis. Notably, DlDDX28 contained a novel auxiliary domain (CAF-1 p150), which might contribute to DNA demethylation in longan early SE. 4 DlDDX genes significantly expressed not only in early SE and zygotic embryogenesis (ZE) but also up-regulated at high levels in 'Honghezi' and 'Quanlongbaihe' with abortive seeds, which are of great significance. Moreover, some DlDDXs presented abiotic stress-response dynamic expression patterns by ABA, SA, JA, and NaCl treatments during early SE. Hence, DEAD-box is essential to SE development and seed abortive in longan.