Hypermethylation of Auxin-Responsive Motifs in the Promoters of the Transcription Factor Genes Accompanies the Somatic Embryogenesis Induction in Arabidopsis

Protein profile changes during priming explants to embryogenic response in Coffea canephora: identification of the RPN12 proteasome subunit involved in the protein degradation

Plant somatic embryogenesis encompasses somatic cells switch into embryogenic cells that can later produce somatic embryos with the ability to produce plantlets. Previously, we defined in vitro culture settings for the somatic embryogenesis process of Coffea canephora that comprise adequate plantlets with auxin plus cytokinin followed by cut-leaf explant cultivation with cytokinin, producing embryos with the ability to regenerate plantlets. Here, we confirmed that cultivating cut-leaf explants with cytokinin is sufficient to promote somatic embryos proliferation and the high yield of somatic embryos in the protocol requires adequate plantlets with auxin plus cytokinin. Two-dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis gels reveal auxin-plus cytokinin-dependent regulated proteins in plantlets with up and down abundance. Chitinase A class III, proteins involved in the metabolism and folding of proteins, photosynthesis, antioxidant activity, and chromatin organization were identified. The RPN12 protein, which is a subunit of the proteasome 26S, has an abundance that is not associated with transcript changes, suggesting post-translational regulation.

Transcriptomic profiling reveals histone acetylation-regulated genes involved in somatic embryogenesis in Arabidopsis thaliana

BACKGROUND

Somatic embryogenesis (SE) exemplifies the unique developmental plasticity of plant cells. The regulatory processes, including epigenetic modifications controlling embryogenic reprogramming of cell transcriptome, have just started to be revealed.

RESULTS

To identify the genes of histone acetylation-regulated expression in SE, we analyzed global transcriptomes of Arabidopsis explants undergoing embryogenic induction in response to treatment with histone deacetylase inhibitor, trichostatin A (TSA). The TSA-induced and auxin (2,4-dichlorophenoxyacetic acid; 2,4-D)-induced transcriptomes were compared. RNA-seq results revealed the similarities of the TSA- and auxin-induced transcriptomic responses that involve extensive deregulation, mostly repression, of the majority of genes. Within the differentially expressed genes (DEGs), we identified the master regulators (transcription factors - TFs) of SE, genes involved in biosynthesis, signaling, and polar transport of auxin and NITRILASE-encoding genes of the function in indole-3-acetic acid (IAA) biosynthesis. TSA-upregulated TF genes of essential functions in auxin-induced SE, included LEC1/LEC2, FUS3, AGL15, MYB118, PHB, PHV, PLTs, and WUS/WOXs. The TSA-induced transcriptome revealed also extensive upregulation of stress-related genes, including those related to stress hormone biosynthesis. In line with transcriptomic data, TSA-induced explants accumulated salicylic acid (SA) and abscisic acid (ABA), suggesting the role of histone acetylation (Hac) in regulating stress hormone-related responses during SE induction. Since mostly the adaxial side of cotyledon explant contributes to SE induction, we also identified organ polarity-related genes responding to TSA treatment, including AIL7/PLT7, RGE1, LBD18, 40, HB32, CBF1, and ULT2. Analysis of the relevant mutants supported the role of polarity-related genes in SE induction.

CONCLUSION

The study results provide a step forward in deciphering the epigenetic network controlling embryogenic transition in somatic cells of plants.

Genome-wide methylation landscape during somatic embryogenesis in Medicago truncatula reveals correlation between Tnt1 retrotransposition and hyperactive methylation regions

Medicago truncatula is a model legume for fundamental research on legume biology and symbiotic nitrogen fixation. Tnt1, a retrotransposon from tobacco, was used to generate insertion mutants in M. truncatula R108. Approximately 21 000 insertion lines have been generated and publicly available. Tnt1 retro-transposition event occurs during somatic embryogenesis (SE), a pivotal process that triggers massive methylation changes. We studied the SE of M. truncatula R108 using leaf explants and explored the dynamic shifts in the methylation landscape from leaf explants to callus formation and finally embryogenesis. Higher cytosine methylation in all three contexts of CG, CHG, and CHH patterns was observed during SE compared to the controls. Higher methylation patterns were observed in assumed promoter regions (~2-kb upstream regions of transcription start site) of the genes, while lowest was recorded in the untranslated regions. Differentially methylated promoter region analysis showed a higher CHH methylation in embryogenesis tissue samples when compared to CG and CHG methylation. Strong correlation (89.71%) was identified between the differentially methylated regions (DMRs) and the site of Tnt1 insertions in M. truncatula R108 and stronger hypermethylation of genes correlated with higher number of Tnt1 insertions in all contexts of CG, CHG, and CHH methylation. Gene ontology enrichment and KEGG pathway enrichment analysis identified genes and pathways enriched in the signal peptide processing, ATP hydrolysis, RNA polymerase activity, transport, secondary metabolites, and nitrogen metabolism pathways. Combined gene expression analysis and methylation profiling showed an inverse relationship between methylation in the DMRs (regions spanning genes) and the expression of genes. Our results show that a dynamic shift in methylation happens during the SE process in the context of CG, CHH and CHG methylation, and the Tnt1 retrotransposition correlates with the hyperactive methylation regions.

The roles of epigenetic regulators in plant regeneration: Exploring patterns amidst complex conditions

Plants possess remarkable capability to regenerate upon tissue damage or optimal environmental stimuli. This ability not only serves as a crucial strategy for immobile plants to survive through harsh environments, but also made numerous modern plant improvements techniques possible. At the cellular level, this biological process involves dynamic changes in gene expression that redirect cell fate transitions. It is increasingly recognized that chromatin epigenetic modifications, both activating and repressive, intricately interact to regulate this process. Moreover, the outcomes of epigenetic regulation on regeneration are influenced by factors such as the differences in regenerative plant species and donor tissue types, as well as the concentration and timing of hormone treatments. In this review, we focus on several well-characterized epigenetic modifications and their regulatory roles in the expression of widely studied morphogenic regulators, aiming to enhance our understanding of the mechanisms by which epigenetic modifications govern plant regeneration.

Suspension cell cultures of Panax vietnamensis as a biotechnological source of ginsenosides: growth, cytology, and ginsenoside profile assessment

Introduction

Panax vietnamensis is a valuable medicinal plant and a source of a broad spectrum of biologically active ginsenosides of different structural groups. Overexploitation and low adaptability to planation cultivation have made this species vulnerable to human pressure and prompted the development of cell cultivation in vitro as a sustainable alternative to harvesting wild plants for their bioactive components. Despite high interest in biotechnological production, little is known about the main factors affecting cell growth and ginsenoside biosynthesis of this species under in vitro conditions. In this study, the potential of cell cultures of P. vietnamensis as a biotechnological source of ginsenosides was was assessed.

Methods

Six suspension cell lines that were developed from different sections of a single rhizome through a multi-step culture optimization process and maintained for over 3 years on media with different mineral salt base and varying contents of auxins and cytokinins. These cell lines were evaluated for productivity parameters and cytological characteristics. Ginsenoside profiles were assessed using a combination of the reversed-phase ultra-high-performance liquid chromatography-Orbitrap-tandem mass spectrometry (UHPLC-Orbitrap-MS/MS) and ultra-performance liquid chromatography-time of flight-mass spectrometry (UPLC-TOF-MS).

Results

All lines demonstrated good growth with a specific growth rate of 0.1-0.2 day-1, economic coefficient of 0.31-0.70, productivity on dry weight (DW) of 0.30-0.83 gDW (L·day)-1, and maximum biomass accumulation varying from 10 to 22 gDW L-1. Ginsenosides of the protopanaxadiol (Rb1, Rb2/Rb3, malonyl-Rb1, and malonyl-Rb2/Rb3), oleanolic acid (R0 and chikusetsusaponin IV), and ocotillol (vinaginsenoside R1) groups and their isomers were identified in cell biomass extracts. Chikusetsusaponin IV was identified in P. vietnamensis cell culture for the first time.

Discussion

These results suggest that suspension cell cultures of Vietnamese ginseng have a high potential for the biotechnological production of biomass containing ginsenosides, particularly of the oleanolic acid and ocotillol groups.

DNA methylome analysis provides insights into gene regulatory mechanism for better performance of rice under fluctuating environmental conditions: epigenomics of adaptive plasticity

MAIN CONCLUSION

Roots play an important role in adaptive plasticity of rice under dry/direct-sown conditions. However, hypomethylation of genes in leaves (resulting in up-regulated expression) complements the adaptive plasticity of Nagina-22 under DSR conditions. Rice is generally cultivated by transplanting which requires plenty of water for irrigation. Such a practice makes rice cultivation a challenging task under global climate change and reducing water availability. However, dry-seeded/direct-sown rice (DSR) has emerged as a resource-saving alternative to transplanted rice (TPR). Though some of the well-adapted local cultivars are used for DSR, only limited success has been achieved in developing DSR varieties mainly because of a limited knowledge of adaptability of rice under fluctuating environmental conditions. Based on better morpho-physiological and agronomic performance of Nagina-22 (N-22) under DSR conditions, N-22 and IR-64 were grown by transplanting and direct-sowing and used for whole genome methylome analysis to unravel the epigenetic basis of adaptive plasticity of rice. Comparative methylome and transcriptome analyses indicated a large number (4078) of genes regulated through DNA methylation/demethylation in N-22 under DSR conditions. Gene × environment interactions play important roles in adaptive plasticity of rice under direct-sown conditions. While genes for pectinesterase, LRK10, C2H2 zinc-finger protein, splicing factor, transposable elements, and some of the unannotated proteins were hypermethylated, the genes for regulation of transcription, protein phosphorylation, etc. were hypomethylated in CG context in the root of N-22, which played important roles in providing adaptive plasticity to N-22 under DSR conditions. Hypomethylation leading to up-regulation of gene expression in the leaf complements the adaptive plasticity of N-22 under DSR conditions. Moreover, differential post-translational modification of proteins and chromatin assembly/disassembly through DNA methylation in CHG context modulate adaptive plasticity of N-22. These findings would help developing DSR cultivars for increased water-productivity and ecological efficiency.

DNA methylation mediates overgrazing-induced clonal transgenerational plasticity

Overgrazing generally induces dwarfism in grassland plants, and these phenotypic traits could be transmitted to clonal offspring even when overgrazing is excluded. However, the dwarfism-transmitted mechanism remains largely unknown, despite generally thought to be enabled by epigenetic modification. To clarify the potential role of DNA methylation on clonal transgenerational effects, we conducted a greenhouse experiment with Leymus chinensis clonal offspring from different cattle/sheep overgrazing histories via the demethylating agent 5-azacytidine. The results showed that clonal offspring from overgrazed (by cattle or sheep) parents were dwarfed and the auxin content of leaves significantly decreased compared to offspring from no-grazed parents'. The 5-azaC application generally increased the auxin content and promoted the growth of overgrazed offspring while inhibited no-grazed offspring growth. Meanwhile, there were similar trends in the expression level of genes related to auxin-responsive target genes (ARF7, ARF19), and signal transduction gene (AZF2). These results suggest that DNA methylation leads to overgrazing-induced plant transgenerational dwarfism via inhibiting auxin signal pathway.

Interplay of phytohormones and epigenetic regulation: A recipe for plant development and plasticity

Both phytohormone signaling and epigenetic mechanisms have long been known to play crucial roles in plant development and plasticity in response to ambient stimuli. Indeed, diverse signaling pathways mediated by phytohormones and epigenetic processes integrate multiple upstream signals to regulate various plant traits. Emerging evidence indicates that phytohormones and epigenetic processes interact at multiple levels. In this review, we summarize the current knowledge of the interplay between phytohormones and epigenetic processes from the perspective of phytohormone biology. We also review chemical regulators used in epigenetic studies and propose strategies for developing novel regulators using multidisciplinary approaches.

Regulatory Small RNAs for a Sustained Eco-Agriculture

Small RNA (sRNA) has become an alternate biotechnology tool for sustaining eco-agriculture by enhancing plant solidity and managing environmental hazards over traditional methods. Plants synthesize a variety of sRNA to silence the crucial genes of pests or plant immune inhibitory proteins and counter adverse environmental conditions. These sRNAs can be cultivated using biotechnological methods to apply directly or through bacterial systems to counter the biotic stress. On the other hand, through synthesizing sRNAs, microbial networks indicate toxic elements in the environment, which can be used effectively in environmental monitoring and management. Moreover, microbes possess sRNAs that enhance the degradation of xenobiotics and maintain bio-geo-cycles locally. Selective bacterial and plant sRNA systems can work symbiotically to establish a sustained eco-agriculture system. An sRNA-mediated approach is becoming a greener tool to replace xenobiotic pesticides, fertilizers, and other chemical remediation elements. The review focused on the applications of sRNA in both sustained agriculture and bioremediation. It also discusses limitations and recommends various approaches toward future improvements for a sustained eco-agriculture system.

Towards a hierarchical gene regulatory network underlying somatic embryogenesis

Genome-editing technologies have advanced in recent years but designing future crops remains limited by current methods of improving somatic embryogenesis (SE) capacity. In this Opinion, we provide an update on the molecular event by which the phytohormone auxin promotes the acquisition of plant cell totipotency through evoking massive changes in transcriptome and chromatin accessibility. We propose that the chromatin states and individual totipotency-related transcription factors (TFs) from disparate gene families organize into a hierarchical gene regulatory network underlying SE. We conclude with a discussion of the practical paths to probe the cellular origin of the somatic embryo and the epigenetic landscape of the totipotent cell state in the era of single-cell genomics.

Integrated transcriptomic and metabolic analyses provide insights into the maintenance of embryogenic potential and the biosynthesis of phenolic acids and flavonoids involving transcription factors in Larix kaempferi (Lamb.) Carr

Somatic embryogenesis (SE) techniques have been established for micropropagation or basic research related to plant development in many conifer species. The frequent occurrence of non-embryogenic callus (NEC) during SE has impose constraints on the application of somatic embryogenesis SE in Larix kaempferi (Lamb.) Carr, but the potential regulatory mechanisms are poorly understood. In this study, integrated transcriptomic and metabolomic analyses were performed in embryogenic callus (EC) and NEC originating from a single immature zygotic embryo to better decipher the key molecular and metabolic mechanisms required for embryogenic potential maintenance. The results showed that a total of 13,842 differentially expressed genes (DEGs) were found in EC and NEC, among which many were enriched in plant hormone signal transduction, starch and sucrose metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis, and the biosynthesis of amino acids pathways. Metabolite profiling showed that 441 differentially accumulated metabolites (DAMs) were identified in EC and NEC. Both EC and NEC had vigorous primary metabolic activities, while most secondary metabolites were upregulated in NEC. Many totipotency-related transcription factor (TF) genes such as BBMs, WUSs, and LEC1 showed higher expression levels in EC compared with NEC, which may result in the higher accumulation of indole 3-acetic acid (IAA) in EC. NEC was characterized by upregulation of genes and metabolites associated with stress responses, such as DEGs involved in jasmonic acid (JA) and ethylene (ETH) biosynthesis and signal transduction pathways, and DEGs and DAMs related to phenylpropanoid and flavonoid biosynthesis. We predicted and analyzed TFs that could target several key co-expressed structural DEGs including two C4H genes, two CcoAOMT genes and three HCT genes involved in phenylpropanoid and flavonoid biosynthesis. Based on the targeted relationship and the co-expression network, two ERFs (Lk23436 and Lk458687), one MYB (Lk34626) and one C2C2-dof (Lk37167) may play an important role in regulating phenolic acid and flavonoid biosynthesis by transcriptionally regulating the expression of these structural genes. This study shows an approach involving integrated transcriptomic and metabolic analyses to obtain insights into molecular events underlying embryogenic potential maintenance and the biosynthesis mechanisms of key metabolites involving TF regulation, which provides valuable information for the improvement of SE efficiency in L. kaempferi.

A Comparative Transcriptome Analysis Reveals the Molecular Mechanisms That Underlie Somatic Embryogenesis in Peaonia ostii ‘Fengdan’

Low propagation rate is the primary problem that limits industry development of tree peony. In this study, a highly efficient regeneration system for tree peony using somatic embryogenesis (SE) was established. The transcriptomes of zygotic embryo explants (S0), non-embryonic callus (S1), embryonic callus (S2), somatic embryos (S3), and regenerated shoots (S4) were analyzed to determine the regulatory mechanisms that underlie SE in tree peony. The differentially expressed genes (DEGs) were identified in the pairwise comparisons of S1-vs-S2 and S1-vs-S3, respectively. The enriched DEGs were primarily involved in hormone signal transduction, stress response and the nucleus (epigenetic modifications). The results indicated that cell division, particularly asymmetric cell division, was enhanced in S3. Moreover, the genes implicated in cell fate determination played central roles in S3. Hormone signal pathways work in concert with epigenetic modifications and stress responses to regulate SE. SERK, WOX9, BBM, FUS3, CUC, and WUS were characterized as the molecular markers for tree peony SE. To our knowledge, this is the first study of the SE of tree peony using transcriptome sequencing. These results will improve our understanding of the molecular mechanisms that underly SE in tree peony and will benefit the propagation and genetic engineering of this plant.

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.

HCV Core protein represses DKK3 expression via epigenetic silencing and activates the Wnt/β-catenin signaling pathway during the progression of HCC

The Wnt/β-catenin signaling pathway is frequently activated in hepatocellular carcinoma (HCC). A number of studies have focused on the aberrant hypermethylation of the DKK family proteins and its role in regulating the activation of specific signaling pathways. However, the exact way by which DKK regulates the signaling pathway caused by Core protein of HCV has not been reported. In the present study, we evaluated the expression level of DKK and its aberrant promoter methylation to investigate the involvement of epigenetic regulation in hepatoma cell lines. The transcription and protein expression of DKK1 was significantly increased, whereas the transcription and protein expression levels of DKK2, DKK3, and DKK4 were significantly decreased following overexpression of Core protein. Pyrosequencing indicated that hypermethylation of DKK3 was increased. This was associated with increased expression of Dnmt1. The investigation of the molecular mechanism indicated that HCV Core protein interacted with Dnmt1, which combined with the promoter of DKK3, leading to methylation of DKK3. Functional studies indicated that Core protein promoted the growth, migration and invasion of cancer cells. However, upregulation of the expression of DKK3 and/or the knockdown of the expression of Dnmt1 inhibited the growth, migration and invasion of cancer cells. Taken together, the data indicated that epigenetic silencing of DKK3 caused by Dnmt1 activated the Wnt/β-catenin pathway in HCV Core-mediated HCC. Therefore, DKK3 may be a potential diagnostic and therapeutic target for HCC.

AGL15 Promotion of Somatic Embryogenesis: Role and Molecular Mechanism

Plants have amazing regenerative properties with single somatic cells, or groups of cells able to give rise to fully formed plants. One means of regeneration is somatic embryogenesis, by which an embryonic structure is formed that "converts" into a plantlet. Somatic embryogenesis has been used as a model for zygotic processes that are buried within layers of maternal tissues. Understanding mechanisms of somatic embryo induction and development are important as a more accessible model for seed development. We rely on seed development not only for most of our caloric intake, but also as a delivery system for engineered crops to meet agricultural challenges. Regeneration of transformed cells is needed for this applied work as well as basic research to understand gene function. Here we focus on a MADS-domain transcription factor, AGAMOUS-Like15 (AGL15) that shows a positive correlation between accumulation levels and capacity for somatic embryogenesis. We relate AGL15 function to other transcription factors, hormones, and epigenetic modifiers involved in somatic embryo development.

Insights into the Histone Acetylation-Mediated Regulation of the Transcription Factor Genes That Control the Embryogenic Transition in the Somatic Cells of Arabidopsis

Somatic embryogenesis (SE), which is a process that involves the in vitro-induced embryogenic reprogramming of plant somatic cells, requires dynamic changes in the cell transcriptome. These changes are fine-tuned by many genetic and epigenetic factors, including posttranslational histone modifications such as histone acetylation. Antagonistically acting enzymes, histone acetyltransferases (HATs) and deacetylases (HDACs), which control histone acetylation in many developmental processes, are believed to control SE. However, the function of specific HAT/HDACs and the genes that are subjected to histone acetylation-mediated regulation during SE have yet to be revealed. Here, we present the global and gene-specific changes in histone acetylation in Arabidopsis explants that are undergoing SE. In the TSA (trichostatin A)-induced SE, we demonstrate that H3 and H4 acetylation might control the expression of the critical transcription factor (TF) genes of a vital role in SE, including LEC1, LEC2 (LEAFY COTYLEDON 1; 2), FUS3 (FUSCA 3) and MYB118 (MYB DOMAIN PROTEIN 118). Within the HATs and HDACs, which mainly positively regulate SE, we identified HDA19 as negatively affecting SE by regulating LEC1, LEC2 and BBM. Finally, we provide some evidence on the role of HDA19 in the histone acetylation-mediated regulation of LEC2 during SE. Our results reveal an essential function of histone acetylation in the epigenetic mechanisms that control the TF genes that play critical roles in the embryogenic reprogramming of plant somatic cells. The results implicate the complexity of Hac-related gene regulation in embryogenic induction and point to differences in the regulatory mechanisms that are involved in auxin- and TSA-induced SE.

Pup1 QTL Regulates Gene Expression Through Epigenetic Modification of DNA Under Phosphate Starvation Stress in Rice

Cytosine methylation, epigenetic DNA modification, is well known to regulate gene expression. Among the epigenetic modifications, 5-methylcytosine (5-mC) has been one of the extensively studied epigenetic changes responsible for regulating gene expression in animals and plants. Though a dramatic change in 5-mC content is observed at the genome level, the variation in gene expression is generally less than that it is expected. Only less is understood about the significance of 5-mC in gene regulation under P-starvation stress in plants. Using whole-genome bisulfite sequencing of a pair of rice [Pusa-44 and its near-isogenic line (NIL)-23 harboring Pup1 QTL] genotypes, we could decipher the role of Pup1 on DNA (de)methylation-mediated regulation of gene expression under P-starvation stress. We observed 13-15% of total cytosines to be methylated in the rice genome, which increased significantly under the stress. The number of differentially methylated regions (DMRs) for hypomethylation (6,068) was higher than those (5,279) for hypermethylated DMRs under the stress, particularly in root of NIL-23. Hypomethylation in CHH context caused upregulated expression of 489 genes in shoot and 382 genes in root of NIL-23 under the stress, wherein 387 genes in shoot and 240 genes in root were upregulated exclusively in NIL-23. Many of the genes for DNA methylation, a few for DNA demethylation, and RNA-directed DNA methylation were upregulated in root of NIL-23 under the stress. Methylation or demethylation of DNA in genic regions differentially affected gene expression. Correlation analysis for the distribution of DMRs and gene expression indicated the regulation of gene mainly through (de)methylation of promoter. Many of the P-responsive genes were hypomethylated or upregulated in roots of NIL-23 under the stress. Hypermethylation of gene body in CG, CHG, and CHH contexts caused up- or downregulated expression of transcription factors (TFs), P transporters, phosphoesterases, retrotransposon proteins, and other proteins. Our integrated transcriptome and methylome analyses revealed an important role of the Pup1 QTL in epigenetic regulation of the genes for transporters, TFs, phosphatases, carbohydrate metabolism, hormone-signaling, and chromatin architecture or epigenetic modifications in P-starvation tolerance. This provides insights into the molecular function of Pup1 in modulating gene expression through DNA (de)methylation, which might be useful in improving P-use efficiency or productivity of rice in P-deficient soil.

Biogenesis, Trafficking, and Function of Small RNAs in Plants

Small RNAs (sRNAs) encoded by plant genomes have received widespread attention because they can affect multiple biological processes. Different sRNAs that are synthesized in plant cells can move throughout the plants, transport to plant pathogens via extracellular vesicles (EVs), and transfer to mammals via food. Small RNAs function at the target sites through DNA methylation, RNA interference, and translational repression. In this article, we reviewed the systematic processes of sRNA biogenesis, trafficking, and the underlying mechanisms of its functions.