Arabidopsis VIM proteins regulate epigenetic silencing by modulating DNA methylation and histone modification in cooperation with MET1

GhVIM28, a negative regulator identified from VIM family genes, positively responds to salt stress in cotton

The VIM (belonged to E3 ubiquitin ligase) gene family is crucial for plant growth, development, and stress responses, yet their role in salt stress remains unclear. We analyzed phylogenetic relationships, chromosomal localization, conserved motifs, gene structure, cis-acting elements, and gene expression patterns of the VIM gene family in four cotton varieties. Our findings reveal 29, 29, 17, and 14 members in Gossypium hirsutum (G.hirsutum), Gossypium barbadense (G.barbadense), Gossypium arboreum (G.arboreum), and Gossypium raimondii (G. raimondii), respectively, indicating the maturity and evolution of this gene family. motifs among GhVIMs genes were observed, along with the presence of stress-responsive, hormone-responsive, and growth-related elements in their promoter regions. Gene expression analysis showed varying patterns and tissue specificity of GhVIMs genes under abiotic stress. Silencing GhVIM28 via virus-induced gene silencing revealed its role as a salt-tolerant negative regulator. This work reveals a mechanism by which the VIM gene family in response to salt stress in cotton, identifying a potential negative regulator, GhVIM28, which could be targeted for enhancing salt tolerance in cotton. The objective of this study was to explore the evolutionary relationship of the VIM gene family and its potential function in salt stress tolerance, and provide important genetic resources for salt tolerance breeding of cotton.

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

How transposable elements are recognized and epigenetically silenced in plants?

Plant genomes are littered with transposable elements (TEs). Because TEs are potentially highly mutagenic, host organisms have evolved a set of defense mechanisms to recognize and epigenetically silence them. Although the maintenance of TE silencing is well studied, our understanding of the initiation of TE silencing is limited, but it clearly involves small RNAs and DNA methylation. Once TEs are silent, the silent state can be maintained to subsequent generations. However, under some circumstances, such inheritance is unstable, leading to the escape of TEs to the silencing machinery, resulting in the transcriptional activation of TEs. Epigenetic control of TEs has been found to be closely linked to many other epigenetic phenomena, such as genomic imprinting, and is known to contribute to regulation of genes, especially those near TEs. Here we review and discuss the current models of TE silencing, its unstable inheritance after hybridization, and the effects of epigenetic regulation of TEs on genomic imprinting.

PHYTOCHROME-INTERACTING FACTORS are involved in starch degradation adjustment via inhibition of the carbon metabolic regulator QUA-QUINE STARCH in Arabidopsis

As sessile organisms, plants encounter dynamic and challenging environments daily, including abiotic/biotic stresses. The regulation of carbon and nitrogen allocations for the synthesis of plant proteins, carbohydrates, and lipids is fundamental for plant growth and adaption to its surroundings. Light, one of the essential environmental signals, exerts a substantial impact on plant metabolism and resource partitioning (i.e., starch). However, it is not fully understood how light signaling affects carbohydrate production and allocation in plant growth and development. An orphan gene unique to Arabidopsis thaliana, named QUA-QUINE STARCH (QQS) is involved in the metabolic processes for partitioning of carbon and nitrogen among proteins and carbohydrates, thus influencing leaf, seed composition, and plant defense in Arabidopsis. In this study, we show that PHYTOCHROME-INTERACTING bHLH TRANSCRIPTION FACTORS (PIFs), including PIF4, are required to suppress QQS during the period at dawn, thus preventing overconsumption of starch reserves. QQS expression is significantly de-repressed in pif4 and pifQ, while repressed by overexpression of PIF4, suggesting that PIF4 and its close homologs (PIF1, PIF3, and PIF5) act as negative regulators of QQS expression. In addition, we show that the evening complex, including ELF3 is required for active expression of QQS, thus playing a positive role in starch catabolism during night-time. Furthermore, QQS is epigenetically suppressed by DNA methylation machinery, whereas histone H3 K4 methyltransferases (e.g., ATX1, ATX2, and ATXR7) and H3 acetyltransferases (e.g., HAC1 and HAC5) are involved in the expression of QQS. This study demonstrates that PIF light signaling factors help plants utilize optimal amounts of starch during the night and prevent overconsumption of starch before its biosynthesis during the upcoming day.

Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Genetic and environmental drivers of large-scale epigenetic variation in Thlaspi arvense

Natural plant populations often harbour substantial heritable variation in DNA methylation. However, a thorough understanding of the genetic and environmental drivers of this epigenetic variation requires large-scale and high-resolution data, which currently exist only for a few model species. Here, we studied 207 lines of the annual weed Thlaspi arvense (field pennycress), collected across a large latitudinal gradient in Europe and propagated in a common environment. By screening for variation in DNA sequence and DNA methylation using whole-genome (bisulfite) sequencing, we found significant epigenetic population structure across Europe. Average levels of DNA methylation were strongly context-dependent, with highest DNA methylation in CG context, particularly in transposable elements and in intergenic regions. Residual DNA methylation variation within all contexts was associated with genetic variants, which often co-localized with annotated methylation machinery genes but also with new candidates. Variation in DNA methylation was also significantly associated with climate of origin, with methylation levels being lower in colder regions and in more variable climates. Finally, we used variance decomposition to assess genetic versus environmental associations with differentially methylated regions (DMRs). We found that while genetic variation was generally the strongest predictor of DMRs, the strength of environmental associations increased from CG to CHG and CHH, with climate-of-origin as the strongest predictor in about one third of the CHH DMRs. In summary, our data show that natural epigenetic variation in Thlaspi arvense is significantly associated with both DNA sequence and environment of origin, and that the relative importance of the two factors strongly depends on the sequence context of DNA methylation. T. arvense is an emerging biofuel and winter cover crop; our results may hence be relevant for breeding efforts and agricultural practices in the context of rapidly changing environmental conditions.

Cellular differentiation, hormonal gradient, and molecular alternation between the division zone and the elongation zone of bamboo internodes

Bamboo differentiates a cell division zone (DZ) and a cell elongation zone (EZ) to promote internode elongation during rapid growth. However, the biological mechanisms underlying this sectioned growth behavior are still unknown. Using histological, physiological, and genomic data, we found that the cell wall and other subcellular organelles such as chloroplasts are more developed in the EZ. Abundant hydrogen peroxide accumulated in the pith cells of the EZ, and stomata formed completely in the EZ. In contrast, most cells in the DZ were in an undifferentiated state with wrinkled cell walls and dense cytoplasm. Hormone detection revealed that the levels of gibberellin, auxin, cytokinin, and brassinosteroid were higher in the DZ than in the EZ. However, the levels of salicylic acid and jasmonic acid were higher in the EZ than in the DZ. Transcriptome analysis with qRT-PCR quantification revealed that the transcripts for cell division and primary metabolism had higher expression in the DZ, whereas the genes for photosynthesis, cell wall growth, and secondary metabolism were dramatically upregulated in the EZ. Overexpression of a MYB transcription factor, BmMYB83, promotes cell wall lignification in transgenic plants. BmMYB83 is specifically expressed in cells that may have lignin deposits, such as protoxylem vessels and fiber cells. Our results indicate that hormone gradient and transcriptome reprogramming, as well as specific expression of key genes such as BmMYB83, may lead to differentiation of cell growth in the bamboo internode.

The Arabidopsis APOLO and human UPAT sequence-unrelated long noncoding RNAs can modulate DNA and histone methylation machineries in plants

BACKGROUND

RNA-DNA hybrid (R-loop)-associated long noncoding RNAs (lncRNAs), including the Arabidopsis lncRNA AUXIN-REGULATED PROMOTER LOOP (APOLO), are emerging as important regulators of three-dimensional chromatin conformation and gene transcriptional activity.

RESULTS

Here, we show that in addition to the PRC1-component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1), APOLO interacts with the methylcytosine-binding protein VARIANT IN METHYLATION 1 (VIM1), a conserved homolog of the mammalian DNA methylation regulator UBIQUITIN-LIKE CONTAINING PHD AND RING FINGER DOMAINS 1 (UHRF1). The APOLO-VIM1-LHP1 complex directly regulates the transcription of the auxin biosynthesis gene YUCCA2 by dynamically determining DNA methylation and H3K27me3 deposition over its promoter during the plant thermomorphogenic response. Strikingly, we demonstrate that the lncRNA UHRF1 Protein Associated Transcript (UPAT), a direct interactor of UHRF1 in humans, can be recognized by VIM1 and LHP1 in plant cells, despite the lack of sequence homology between UPAT and APOLO. In addition, we show that increased levels of APOLO or UPAT hamper VIM1 and LHP1 binding to YUCCA2 promoter and globally alter the Arabidopsis transcriptome in a similar manner.

CONCLUSIONS

Collectively, our results uncover a new mechanism in which a plant lncRNA coordinates Polycomb action and DNA methylation through the interaction with VIM1, and indicates that evolutionary unrelated lncRNAs with potentially conserved structures may exert similar functions by interacting with homolog partners.

miR778 mediates gene expression, histone modification, and DNA methylation during cyst nematode parasitism

Despite the known critical regulatory functions of microRNAs, histone modifications, and DNA methylation in reprograming plant epigenomes in response to pathogen infection, the molecular mechanisms underlying the tight coordination of these components remain poorly understood. Here, we show how Arabidopsis (Arabidopsis thaliana) miR778 coordinately modulates the root transcriptome, histone methylation, and DNA methylation via post-transcriptional regulation of the H3K9 methyltransferases SU(var)3-9 homolog 5 (SUVH5) and SUVH6 upon infection by the beet cyst nematode Heterodera schachtii. miR778 post-transcriptionally silences SUVH5 and SUVH6 upon nematode infection. Manipulation of the expression of miR778 and its two target genes significantly altered plant susceptibility to H. schachtii. RNA-seq analysis revealed a key role of SUVH5 and SUVH6 in reprograming the transcriptome of Arabidopsis roots upon H. schachtii infection. In addition, chromatin immunoprecipitation (ChIP)-seq analysis established SUVH5 and SUVH6 as the main enzymes mediating H3K9me2 deposition in Arabidopsis roots in response to nematode infection. ChIP-seq analysis also showed that these methyltransferases possess distinct DNA binding preferences in that they are targeting transposable elements under noninfected conditions and protein-coding genes in infected plants. Further analyses indicated that H3K9me2 deposition directed by SUVH5 and SUVH6 contributes to gene expression changes both in roots and in nematode feeding sites and preferentially associates with CG DNA methylation. Together, our results uncovered multi-layered epigenetic regulatory mechanisms coordinated by miR778 during Arabidopsis-H. schachtii interactions.

A tale of two players: the role of phosphate in iron and zinc homeostatic interactions

This minireview details the impact of iron-phosphate and zinc-phosphate interactions in plants and provides perspectives for further areas of research regarding nutrient homeostasis. Iron (Fe) and zinc (Zn) are among the most important micronutrients for plant growth and have numerous implications for human health and agriculture. While plants have developed efficient uptake and transport mechanisms for Fe and Zn, emerging research has shown that the availability of other nutrients in the environment influences the homeostasis of Fe and Zn within plants. In this minireview, we present the current knowledge regarding homeostatic interactions of Fe and Zn with the macronutrient phosphorous (P) and the resulting physiological responses to combined deficiencies of these nutrients. Fe and P interactions have been shown to influence root development, photosynthesis, and biological processes aiding Fe uptake. Zn and P interactions also influence root growth, and coordination of Zn-dependent transcriptional regulation contributes to phosphate (Pi) transport in the plant. Understanding homeostatic interactions among these different nutrients is of critical importance to obtain a more complete understanding of plant nutrition in complex soil environments.

GWAS on multiple traits identifies mitochondrial ACONITASE3 as important for acclimation to submergence stress

Flooding causes severe crop losses in many parts of the world. Genetic variation in flooding tolerance exists in many species; however, there are few examples for the identification of tolerance genes and their underlying function. We conducted a genome-wide association study (GWAS) in 387 Arabidopsis (Arabidopsis thaliana) accessions. Plants were subjected to prolonged submergence followed by desubmergence, and seven traits (score, water content, Fv/Fm, and concentrations of nitrate, chlorophyll, protein, and starch) were quantified to characterize their acclimation responses. These traits showed substantial variation across the range of accessions. A total of 35 highly significant single-nucleotide polymorphisms (SNPs) were identified across the 20 GWA datasets, pointing to 22 candidate genes, with functions in TCA cycle, DNA modification, and cell division. Detailed functional characterization of one candidate gene, ACONITASE3 (ACO3), was performed. Chromatin immunoprecipitation followed by sequencing showed that a single nucleotide polymorphism in the ACO3 promoter co-located with the binding site of the master regulator of retrograde signaling ANAC017, while subcellular localization of an ACO3-YFP fusion protein confirmed a mitochondrial localization during submergence. Analysis of mutant and overexpression lines determined changes in trait parameters that correlated with altered submergence tolerance and were consistent with the GWAS results. Subsequent RNA-seq experiments suggested that impairing ACO3 function increases the sensitivity to submergence by altering ethylene signaling, whereas ACO3 overexpression leads to tolerance by metabolic priming. These results indicate that ACO3 impacts submergence tolerance through integration of carbon and nitrogen metabolism via the mitochondrial TCA cycle and impacts stress signaling during acclimation to stress.

Multi-omics data integration reveals link between epigenetic modifications and gene expression in sugar beet (Beta vulgaris subsp. vulgaris) in response to cold

Background

DNA methylation is thought to influence the expression of genes, especially in response to changing environmental conditions and developmental changes. Sugar beet (Beta vulgaris ssp. vulgaris), and other biennial or perennial plants are inevitably exposed to fluctuating temperatures throughout their lifecycle and might even require such stimulus to acquire floral competence. Therefore, plants such as beets, need to fine-tune their epigenetic makeup to ensure phenotypic plasticity towards changing environmental conditions while at the same time steering essential developmental processes. Different crop species may show opposing reactions towards the same abiotic stress, or, vice versa, identical species may respond differently depending on the specific kind of stress.

Results

In this study, we investigated common effects of cold treatment on genome-wide DNA methylation and gene expression of two Beta vulgaris accessions via multi-omics data analysis. Cold exposure resulted in a pronounced reduction of DNA methylation levels, which particularly affected methylation in CHH context (and to a lesser extent CHG) and was accompanied by transcriptional downregulation of the chromomethyltransferase CMT2 and strong upregulation of several genes mediating active DNA demethylation.

Conclusion

Integration of methylomic and transcriptomic data revealed that, rather than methylation having directly influenced expression, epigenetic modifications correlated with changes in expression of known players involved in DNA (de)methylation. In particular, cold triggered upregulation of genes putatively contributing to DNA demethylation via the ROS1 pathway. Our observations suggest that these transcriptional responses precede the cold-induced global DNA-hypomethylation in non-CpG, preparing beets for additional transcriptional alterations necessary for adapting to upcoming environmental changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08312-2.

The CLASSY family controls tissue-specific DNA methylation patterns in Arabidopsis

DNA methylation shapes the epigenetic landscape of the genome, plays critical roles in regulating gene expression, and ensures transposon silencing. As is evidenced by the numerous defects associated with aberrant DNA methylation landscapes, establishing proper tissue-specific methylation patterns is critical. Yet, how such differences arise remains a largely open question in both plants and animals. Here we demonstrate that CLASSY1-4 (CLSY1-4), four locus-specific regulators of DNA methylation, also control tissue-specific methylation patterns, with the most striking pattern observed in ovules where CLSY3 and CLSY4 control DNA methylation at loci with a highly conserved DNA motif. On a more global scale, we demonstrate that specific clsy mutants are sufficient to shift the epigenetic landscape between tissues. Together, these findings reveal substantial epigenetic diversity between tissues and assign these changes to specific CLSY proteins, elucidating how locus-specific targeting combined with tissue-specific expression enables the CLSYs to generate epigenetic diversity during plant development.

CLASSY (CLSY) proteins regulate DNA methylation at specific loci in the Arabidopsis genome. Here the authors show that the CLSYs also control tissue-specific DNA methylation, including at siren loci in ovules, and that the lack of an individual CLSYs can shift the epigenetic landscape between tissues.

The UHRF protein family in epigenetics, development, and carcinogenesis

The UHRF protein family consists of multidomain regulatory proteins that sense modification status of DNA and/or proteins and catalyze the ubiquitylation of target proteins. Through their functional domains, they interact with other molecules and serve as a hub for regulatory networks of several important biological processes, including maintenance of DNA methylation and DNA damage repair. The UHRF family is conserved in vertebrates and plants but is missing from fungi and many nonvertebrate animals. Mammals commonly have UHRF1 and UHRF2, but, despite their high structural similarity, the two paralogues appear to have distinct functions. Furthermore, UHRF1 and UHRF2 show different expression patterns and different outcomes in gene knockout experiments. In this review, we summarize the current knowledge on the molecular function of the UHRF family in various biological pathways and discuss their roles in epigenetics, development, gametogenesis, and carcinogenesis, with a focus on the mammalian UHRF proteins.

Expression Patterns of DNA Methylation and Demethylation Genes during Plant Development and in Response to Phytohormones

DNA methylation and demethylation precisely and effectively modulate gene expression during plant growth and development and in response to stress. However, expression profiles of genes involved in DNA methylation and demethylation during plant development and their responses to phytohormone treatments remain largely unknown. We characterized the spatiotemporal expression patterns of genes involved in de novo methylation, methyl maintenance, and active demethylation in roots, shoots, and reproductive organs using β-glucuronidase (GUS) reporter lines. Promoters of DNA demethylases were generally more highly active at the mature root tissues, whereas the promoters of genes involved in DNA methylation were more highly active at fast-growing root tissues. The promoter activity also implies that methylation status in shoot apex, leaf primordia, floral organs, and developing embryos is under tight equilibrium through the activity of genes involved in DNA methylation and demethylation. The promoter activity of DNA methylation and demethylation-related genes in response to various phytohormone treatments revealed that phytohormones can alter DNA methylation status in specific and redundant ways. Overall, our results illustrate that DNA methylation and demethylation pathways act synergistically and antagonistically in various tissues and in response to phytohormone treatments and point to the existence of hormone-linked methylome regulation mechanisms that may contribute to tissue differentiation and development.

Helical and β-Turn Conformations in the Peptide Recognition Regions of the VIM1 PHD Finger Abrogate H3K4 Peptide Recognition

The PHD finger-containing VARIANT IN METHYLATION/ORTHRUS (VIM/ORTH) family of proteins in Arabidopsis consists of functional homologues of mammalian UHRF1 and is required for the maintenance of DNA methylation. Comparison of the sequence with those of other PHD fingers implied that VIM1 and VIM3 PHD could recognize lysine 4 of histone H3 (H3K4) through interactions mediated by a conserved aspartic acid. However, our calorimetric and modified histone peptide array binding studies suggested that neither H3K4 nor other histone marks are recognized by VIM1 and VIM3 PHD fingers. Here, we report a 2.6 Å resolution crystal structure of the VIM1 PHD finger and demonstrate significant structural changes in the putative H3 recognition segments in contrast to canonical H3K4 binding PHD fingers. These changes include (i) the H3A1 binding region, (ii) strand β1 that forms an intermolecular β-sheet with the H3 peptide, and (iii) an aspartate-containing motif involved in salt bridge interaction with H3K4, which together appear to abrogate recognition of H3K4 by the VIM1 PHD finger. To understand the significance of the altered structural features in the VIM1 PHD that might prevent histone H3 recognition, we modeled a chimeric VIM1 PHD (chmVIM1 PHD) by grafting the peptide binding structural features of the BHC80 PHD onto the VIM1 PHD. Molecular dynamics simulation and metadynamics analyses revealed that the chmVIM1 PHD-H3 complex is stable and also showed a network of intermolecular interactions similar to those of the BHC80 PHD-H3 complex. Collectively, this study reveals that subtle structural changes in the peptide binding region of the VIM1 PHD abrogate histone H3 recognition.

A single-cell view of the transcriptome during lateral root initiation in Arabidopsis thaliana

Root architecture is a major determinant of plant fitness and is under constant modification in response to favorable and unfavorable environmental stimuli. Beyond impacts on the primary root, the environment can alter the position, spacing, density, and length of secondary or lateral roots. Lateral root development is among the best-studied examples of plant organogenesis, yet there are still many unanswered questions about its earliest steps. Among the challenges faced in capturing these first molecular events is the fact that this process occurs in a small number of cells with unpredictable timing. Single-cell sequencing methods afford the opportunity to isolate the specific transcriptional changes occurring in cells undergoing this fate transition. Using this approach, we successfully captured the transcriptomes of initiating lateral root primordia in Arabidopsis thaliana and discovered many upregulated genes associated with this process. We developed a method to selectively repress target gene transcription in the xylem pole pericycle cells where lateral roots originate and demonstrated that the expression of several of these targets is required for normal root development. We also discovered subpopulations of cells in the pericycle and endodermal cell files that respond to lateral root initiation, highlighting the coordination across cell files required for this fate transition.

The multi-functionality of UHRF1: epigenome maintenance and preservation of genome integrity

During S phase, the cooperation between the macromolecular complexes regulating DNA synthesis, epigenetic information maintenance and DNA repair is advantageous for cells, as they can rapidly detect DNA damage and initiate the DNA damage response (DDR). UHRF1 is a fundamental epigenetic regulator; its ability to coordinate DNA methylation and histone code is unique across proteomes of different species. Recently, UHRF1’s role in DNA damage repair has been explored and recognized to be as important as its role in maintaining the epigenome. UHRF1 is a sensor for interstrand crosslinks and a determinant for the switch towards homologous recombination in the repair of double-strand breaks; its loss results in enhanced sensitivity to DNA damage. These functions are finely regulated by specific post-translational modifications and are mediated by the SRA domain, which binds to damaged DNA, and the RING domain. Here, we review recent studies on the role of UHRF1 in DDR focusing on how it recognizes DNA damage and cooperates with other proteins in its repair. We then discuss how UHRF1’s epigenetic abilities in reading and writing histone modifications, or its interactions with ncRNAs, could interlace with its role in DDR.

F-box protein CFK1 interacts with and degrades de novo DNA methyltransferase in Arabidopsis

DNA methylation plays crucial roles in cellular development and stress responses through gene regulation and genome stability control. Precise regulation of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), the de novo Arabidopsis DNA methyltransferase, is crucial to maintain DNA methylation homeostasis to ensure genome integrity. Compared with the extensive studies on DRM2 targeting mechanisms, little information is known regarding the quality control of DRM2 itself. Here, we conducted yeast two-hybrid screen assay and identified an E3 ligase, COP9 INTERACTING F-BOX KELCH 1 (CFK1), as a novel DRM2-interacting partner and targets DRM2 for degradation via the ubiquitin-26S proteasome pathway in Arabidopsis thaliana. We also performed whole genome bisulfite sequencing (BS-seq) to determine the biological significance of CFK1-mediated DRM2 degradation. Loss-of-function CFK1 leads to increased DRM2 protein abundance and overexpression of CFK1 showed reduced DRM2 protein levels. Consistently, CFK1 overexpression induces genome-wide CHH hypomethylation and transcriptional de-repression at specific DRM2 target loci. This study uncovered a distinct mechanism regulating de novo DNA methyltransferase by CFK1 to control DNA methylation level.

Genome-Wide Association Study on Seminal and Nodal Roots of Wheat Under Different Growth Environments

The root of wheat consists of seminal and nodal roots. Comparatively speaking, fewer studies have been carried out on the nodal root system because of its disappearance at the early seedling stage under indoor environments. In this study, 196 accessions from the Huanghuai Wheat Region (HWR) were used to identify the characteristics of seminal and nodal root traits under different growth environments, including indoor hydroponic culture (IHC), outdoor hydroponic culture (OHC), and outdoor pot culture (OPC), for three growing seasons. The results indicated that the variation range of root traits in pot environment was larger than that in hydroponic environment, and canonical coefficients were the greatest between OHC and OPC (0.86) than those in other two groups, namely, IHC vs. OPC (0.48) and IHC vs. OHC (0.46). Most root traits were negatively correlated with spikes per area (SPA), grains per spike (GPS), and grain yield (GY), while all the seminal root traits were positively correlated with thousand-kernel weight (TKW). Genome-wide association study (GWAS) was carried out on root traits by using a wheat 660K SNP array. A total of 35 quantitative trait loci (QTLs)/chromosomal segments associated with root traits were identified under OPC and OHC. In detail, 11 and 24 QTLs were significantly associated with seminal root and nodal root traits, respectively. Moreover, 13 QTLs for number of nodal roots per plant (NRP) containing 14 stable SNPs, were distributed on chromosomes 1B, 2B, 3A, 4B, 5D, 6D, 7A, 7B, and Un. Based on LD and bioinformatics analysis, these QTLs may contain 17 genes closely related to NRP. Among them, TraesCS2B02G552500 and TraesCS7A02G428300 were highly expressed in root tissues. Moreover, the frequencies of favorable alleles of these 14 SNPs were confirmed to be less than 70% in the natural population, suggesting that the utilization of these superior genes in wheat root is still improving.

Small DNA Methylation, Big Player in Plant Abiotic Stress Responses and Memory

DNA methylation is a conserved epigenetic mark that plays important roles in maintaining genome stability and regulating gene expression. As sessile organisms, plants have evolved sophisticated regulatory systems to endure or respond to diverse adverse abiotic environmental challenges, i.e., abiotic stresses, such as extreme temperatures (cold and heat), drought and salinity. Plant stress responses are often accompanied by changes in chromatin modifications at diverse responsive loci, such as 5-methylcytosine (5mC) and N 6-methyladenine (6mA) DNA methylation. Some abiotic stress responses are memorized for several hours or days through mitotic cell divisions and quickly reset to baseline levels after normal conditions are restored, which is referred to as somatic memory. In some cases, stress-induced chromatin marks are meiotically heritable and can impart the memory of stress exposure from parent plants to at least the next stress-free offspring generation through the mechanisms of transgenerational epigenetic inheritance, which may offer the descendants the potential to be adaptive for better fitness. In this review, we briefly summarize recent achievements regarding the establishment, maintenance and reset of DNA methylation, and highlight the diverse roles of DNA methylation in plant responses to abiotic stresses. Further, we discuss the potential role of DNA methylation in abiotic stress-induced somatic memory and transgenerational inheritance. Future research directions are proposed to develop stress-tolerant engineered crops to reduce the negative effects of abiotic stresses.

DNA Geminivirus Infection Induces an Imprinted E3 Ligase Gene to Epigenetically Activate Viral Gene Transcription

Flowering plants and mammals contain imprinted genes that are primarily expressed in the endosperm and placenta in a parent-of-origin manner. In this study, we show that early activation of the geminivirus genes C2 and C3 in Arabidopsis (Arabidopsis thaliana) plants, encoding a viral suppressor of RNA interference and a replication enhancer protein, respectively, is correlated with the transient vegetative expression of VARIANT IN METHYLATION5 (VIM5), an endosperm imprinted gene that is conserved in diverse plant species. VIM5 is a ubiquitin E3 ligase that directly targets the DNA methyltransferases MET1 and CMT3 for degradation by the ubiquitin-26S proteasome proteolytic pathway. Infection with Beet severe curly top virus induced VIM5 expression in rosette leaf tissues, possibly via the expression of the viral replication initiator protein, leading to the early activation of C2 and C3 coupled with reduced symmetric methylation in the C2-3 promoter and the onset of disease symptoms. These findings demonstrate how this small DNA virus recruits a host imprinted gene for the epigenetic activation of viral gene transcription. Our findings reveal a distinct strategy used by plant pathogens to exploit the host machinery in order to inhibit methylation-mediated defense responses when establishing infection.

Conservation of centromeric histone 3 interaction partners in plants

The loading and maintenance of centromeric histone 3 (CENH3) at the centromere are critical processes ensuring appropriate kinetochore establishment and equivalent segregation of the homologous chromosomes during cell division. CENH3 loss of function is lethal, whereas mutations in the histone fold domain are tolerated and lead to chromosome instability and chromosome elimination in embryos derived from crosses with wild-type pollen. A wide range of proteins in yeast and animals have been reported to interact with CENH3. The histone fold domain-interacting proteins are potentially alternative targets for the engineering of haploid inducer lines, which may be important when CENH3 mutations are not well supported by a given crop. Here, we provide an overview of the corresponding plant orthologs or functional homologs of CENH3-interacting proteins. We also list putative CENH3 post-translational modifications that are also candidate targets for modulating chromosome stability and inheritance.

DREAM complex suppresses DNA methylation maintenance genes and precludes DNA hypermethylation

The DNA methyltransferases MET1 and CMT3 are known to be responsible for maintenance of DNA methylation at symmetric CG and CHG sites, respectively, in Arabidopsis thaliana. However, it is unknown how the expression of methyltransferase genes is regulated in different cell states and whether change in expression affects DNA methylation at the whole-genome level. Using a reverse genetic screen, we identified TCX5, a tesmin/TSO1-like CXC domain-containing protein, and demonstrated that it is a transcriptional repressor of genes required for maintenance of DNA methylation, which include MET1, CMT3, DDM1, KYP and VIMs. TCX5 functions redundantly with its paralogue TCX6 in repressing the expression of these genes. In the tcx5 tcx6 double mutant, expression of these genes is markedly increased, thereby leading to markedly increased DNA methylation at CHG sites and, to a lesser extent, at CG sites at the whole-genome level. Furthermore, our whole-genome DNA methylation analysis indicated that the CG and CHG methylation level is lower in differentiated quiescent cells than in dividing cells in the wild type but is comparable in the tcx5/6 mutant, suggesting that TCX5/6 are required for maintenance of the difference in DNA methylation between the two cell types. We identified TCX5/6-containing multi-subunit complexes, which are known as DREAM in other eukaryotes, and demonstrated that the Arabidopsis DREAM components function as a whole to preclude DNA hypermethylation. Given that the DREAM complexes are conserved from plants to animals, the preclusion of DNA hypermethylation by DREAM complexes may represent a conserved mechanism in eukaryotes.

Epigenetic memory marks determine epiallele stability at loci targeted by de novo DNA methylation

It is generally assumed that DNA methylation changes at genomic regions targeted by the de novo RNA-directed DNA methylation (RdDM) pathway are unstable. Here, we show that RdDM targets in Arabidopsis can be classified into two groups on the basis of whether there is remethylation following the restoration of NRPD1 function in nrpd1 mutant plants-remethylable loci and non-remethylable loci. In contrast to the remethylable loci, the non-remethylable loci contain higher levels of the euchromatic marks of trimethylation at Lys 4 of histone H3 (H3K4me3), which interferes with the recruitment of the RdDM molecular machinery, and acetylation at Lys 18 of histone H3 (H3K18ac), which helps to recruit the DNA demethylase ROS1 to antagonize RdDM. Here, using targeted methylation erasure by CRISPR-dCas9-TET1, we demonstrate that methylated CG (mCG) and mCHG (where H represents A, C or T) are memory marks that are required for targeting the RdDM machinery to remethylable loci. Our results show that histone and DNA methylation marks are critical in determining the ability of RdDM target loci to form stable epialleles, and contribute to understanding the formation and transmission of epialleles.

Examining Short-Term Responses to a Long-Term Problem: RNA-Seq Analyses of Iron Deficiency Chlorosis Tolerant Soybean

Iron deficiency chlorosis (IDC) is a global crop production problem, significantly impacting yield. However, most IDC studies have focused on model species, not agronomically important crops. Soybean is the second largest crop grown in the United States, yet the calcareous soils across most of the upper U.S. Midwest limit soybean growth and profitability. To understand early soybean iron stress responses, we conducted whole genome expression analyses (RNA-sequencing) of leaf and root tissue from the iron efficient soybean (Glycine max) cultivar Clark, at 30, 60 and 120 min after transfer to iron stress conditions. We identified over 10,000 differentially expressed genes (DEGs), with the number of DEGs increasing over time in leaves, but decreasing over time in roots. To investigate these responses, we clustered our expression data across time to identify suites of genes, their biological functions, and the transcription factors (TFs) that regulate their expression. These analyses reveal the hallmarks of the soybean iron stress response (iron uptake and homeostasis, defense, and DNA replication and methylation) can be detected within 30 min. Furthermore, they suggest root to shoot signaling initiates early iron stress responses representing a novel paradigm for crop stress adaptations.

DNA Methylation Readers in Plants

In plants, DNA methylation occurs in distinct sequence contexts, including CG, CHG, and CHH. Thus, plants have developed a surprisingly diverse set of DNA methylation readers to cope with an extended repertoire of methylated sites. The Arabidopsis genome contains twelve Methyl-Binding Domain proteins (MBD), and nine SET and RING finger-associated (SRA) domain containing proteins belonging to the SUVH clade, in addition to three homologs of UHRF1, namely VIM1-3, all containing SRA domains. In this review, we will highlight several research questions that remain unresolved with respect to the function of plant DNA methylation readers, which can have both de novo demethylase and maintenance activity. We argue that maintenance of CG methylation in plants likely involved actors not found in their mammalian counterparts, and that new evidence suggests significant reprogramming of DNA methylation during plant reproduction as an important new development in the field.

DNA Methylation of Farnesyl Pyrophosphate Synthase, Squalene Synthase, and Squalene Epoxidase Gene Promoters and Effect on the Saponin Content of Eleutherococcus Senticosus

Eleutherococcus senticosus (Ruper. et Maxim.) Maxim is a traditional Chinese medicine. The saponin components of E. senticosus have several biological effects, including reduction of blood lipids; protection against liver, heart, and vascular disease; and antitumor activity. The DNA methylation of E. senticosus farnesyl pyrophosphate synthase (FPS), squalene synthase (SS), and squalene epoxidase (SE) gene promoters and the mechanism of the influence of these enzymes on saponin synthesis and accumulation in E. senticosus were explored using bisulfite sequencing technology, real-time PCR, the vanillin-concentrated sulfuric acid chromogenic method, and LC-MS. There are 19 DNA methylation sites and 8 methylation types in the FPS gene. The SS gene has nine DNA methylation sites and two DNA methylation types. The SE gene has 16 DNA methylation sites and 7 methylation types. The total saponin content in the high and low DNA methylation groups were 1.07 ± 0.12 and 2.92 ± 0.32 mg/g, respectively. Statistical analysis indicated that the gene expression of the FPS, SS, and SE genes was significantly positively correlated with the saponin content (p < 0.05), and that the methylation ratio was significantly negatively correlated with the saponin content (p < 0.01), while the expression of the SS and SE genes was significantly positively correlated (p < 0.01). A total of 488 metabolites were detected from E. senticosus and 100 different metabolites were screened out by extensive targeted metabolomics. The amount of most metabolites related to the mevalonate pathway was higher in the low DNA methylation group than in the high DNA methylation group. It was demonstrated that there are DNA methylation sites in the promoter regions of the FPS, SS, and SE genes of E. senticosus, and DNA methylation in this region could significantly inhibit synthesis in the mevalonate pathway, thus reducing the content of the final product E. senticosus saponin.

Changes in grapevine DNA methylation and polyphenols content induced by solar ultraviolet-B radiation, water deficit and abscisic acid spray treatments

Environment and crop management shape plant's phenotype. Argentinean high-altitude vineyards are characterized by elevated solar ultraviolet-B radiation (UVB) and water deficit (D) that enhance enological quality for red winemaking. These signals promote phenolics accumulation in leaves and berries, being the responses mediated by abscisic acid (ABA). DNA methylation is an epigenetic mechanism that regulates gene expression and may affect grapevine growth, development and acclimation, since methylation patterns are mitotically heritable. Berry skins low molecular weight polyphenols (LMWP) were characterized in field grown Vitis vinifera L. cv. Malbec plants exposed to contrasting UV-B, D, and ABA treatments during one season. The next season early fruit shoots were epigenetically (methylation-sensitive amplification polymorphism; MSAP) and biochemically (LMWP) characterized. Unstable epigenetic patterns and/or stochastic stress-induced methylation changes were observed. UV-B and D were the treatments that induced greater number of DNA methylation changes respect to Control; and UV-B promoted global hypermethylation of MSAP epiloci. Sequenced MSAP fragments associated with UV-B and ABA showed similarities with transcriptional regulators and ubiquitin ligases proteins activated by light. UV-B was associated with flavonols accumulation in berries and with hydroxycinnamic acids in the next season fruit shoots, suggesting that DNA methylation could regulate the LMWP accumulation and participate in acclimation mechanisms.

EFFECTOR OF TRANSCRIPTION factors are novel plant-specific regulators associated with genomic DNA methylation in Arabidopsis

Plant-specific EFFECTORS OF TRANSCRIPTION (ET) are characterised by a variable number of highly conserved ET repeats, which are involved in zinc and DNA binding. In addition, ETs share a GIY-YIG domain, involved in DNA nicking activity. It was hypothesised that ETs might act as epigenetic regulators. Here, methylome, transcriptome and phenotypic analyses were performed to investigate the role of ET factors and their involvement in DNA methylation in Arabidopsis thaliana. Comparative DNA methylation and transcriptome analyses in flowers and seedlings of et mutants revealed ET-specific differentially expressed genes and mostly independently characteristic, ET-specific differentially methylated regions. Loss of ET function results in pleiotropic developmental defects. The accumulation of cyclobutane pyrimidine dimers after ultraviolet stress in et mutants suggests an ET function in DNA repair.

BAH domains and a histone-like motif in DNA methyltransferase 1 (DNMT1) regulate de novo and maintenance methylation in vivo

DNA methyltransferase 1 (DNMT1) is a multidomain protein believed to be involved only in the passive transmission of genomic methylation patterns via maintenance methylation. The mechanisms that regulate DNMT1 activity and targeting are complex and poorly understood. We used embryonic stem (ES) cells to investigate the function of the uncharacterized bromo-adjacent homology (BAH) domains and the glycine-lysine (GK) repeats that join the regulatory and catalytic domains of DNMT1. We removed the BAH domains by means of a CRISPR/Cas9-mediated deletion within the endogenous Dnmt1 locus. The internally deleted protein failed to associate with replication foci during S phase in vivo and lost the ability to mediate maintenance methylation. The data indicate that ablation of the BAH domains causes DNMT1 to be excluded from replication foci even in the presence of the replication focus-targeting sequence (RFTS). The GK repeats resemble the N-terminal tails of histones H2A and H4 and are normally acetylated. Substitution of lysines within the GK repeats with arginines to prevent acetylation did not alter the maintenance activity of DNMT1 but unexpectedly activated de novo methylation of paternal imprinting control regions (ICRs) in mouse ES cells; maternal ICRs remained unmethylated. We propose a model under which DNMT1 deposits paternal imprints in male germ cells in an acetylation-dependent manner. These data reveal that DNMT1 responds to multiple regulatory inputs that control its localization as well as its activity and is not purely a maintenance methyltransferase but can participate in the de novo methylation of a small but essential compartment of the genome.

Heterologous Expression of the Transcription Factor EsNAC1 in Arabidopsis Enhances Abiotic Stress Resistance and Retards Growth by Regulating the Expression of Different Target Genes

Heterologous expression of a transcription factor (TF) gene in a related species is a useful method for crop breeding and the identification of gene function. The differences in phenotype and target gene expression between HE lines (with the heterologous expression of an ortholog) and OX lines (with an overexpressed native gene) must be understood. EsNAC1, encoding a NAC protein and the ortholog of RD26 in Arabidopsis, was cloned from Eutrema salsugineum and introduced into Arabidopsis. The heterologous expression of EsNAC1 retarded the vegetative growth of Arabidopsis, and the transgenic plants (HE lines) showed much greater resistance to salt and oxidative stress than the wild type, Col-0. The HE lines accumulated 2.8-fold (8-h light) of starch, 1.42-fold of Chlorophyll a and 1.31-fold of Chlorophyll b than Col-0 during the light period, with obvious differences compared to the RD26OX line. A genome-wide ChIP (chromatin immunoprecipitation analysis)-on-chip assay revealed that EsNAC1 targeted promoters of different genes compared to RD26. In HE lines, EsNAC1 could specifically upregulate the expression level of TF genes NAC DOMAIN CONTAINING PROTEIN 62 (ANAC062), INTEGRASE-TYPE DNA-BINDING PROTEIN (TINY2), and MYB HYPOCOTYL ELONGATION-RELATED (MYBH) to show more effective abiotic stress resistance than RD26OX lines. Moreover, DELTA1-PYRROLINE-5-CARBOXYLATE SYNTHASE 1 (P5CS1), TRYPTOPHAN BIOSYNTHESIS 2 (TRP2) or GALACTINOL SYNTHASE 2 (GOLS2), was also specifically regulated by EsNAC1 to retard the vegetative growth of HE lines, but not the brassinosteroid singling pathway in RD26OX lines. These differences in phenotypes and metabolism between the HE lines and the RD26OX line implied that the differential features could be produced from the diversity of target genes in the transgenic plants when the ortholog was introduced.

ORESARA15, a PLATZ transcription factor, mediates leaf growth and senescence in Arabidopsis

Plant leaves undergo a series of developmental changes from leaf primordium initiation through growth and maturation to senescence throughout their life span. Although the mechanisms underlying leaf senescence have been intensively elucidated, our knowledge of the interrelationship between early leaf development and senescence is still fragmentary. We isolated the oresara15-1Dominant (ore15-1D) mutant, which had an extended leaf longevity and an enlarged leaf size, from activation-tagged lines of Arabidopsis. Plasmid rescue identified that ORE15 encodes a PLANT A/T-RICH SEQUENCE- AND ZINC-BINDING PROTEIN family transcription factor. Phenotypes of ore15-1D and ore15-2, a loss-of-function mutant, were evaluated through physiological and anatomical analyses. Microarray, quantitative reverse transcription polymerase chain reaction, and chromatin immunoprecipitation as well as genetic analysis were employed to reveal the molecular mechanism of ORE15 in the regulation of leaf growth and senescence. ORE15 enhanced leaf growth by promoting the rate and duration of cell proliferation in the earlier stage and suppressed leaf senescence in the later stage by modulating the GROWTH-REGULATING FACTOR (GRF)/GRF-INTERACTING FACTOR regulatory pathway. Our study highlighted a molecular conjunction through ORE15 between growth and senescence, which are two temporally separate developmental processes during leaf life span.

Arabidopsis Zinc-Finger-Like Protein ASYMMETRIC LEAVES2 (AS2) and Two Nucleolar Proteins Maintain Gene Body DNA Methylation in the Leaf Polarity Gene ETTIN (ARF3)

Arabidopsis ASYMMETRIC LEAVES2 (AS2) plays a critical role in leaf adaxial-abaxial partitioning by repressing expression of the abaxial-determining gene ETTIN/AUXIN RESPONSE FACTOR3 (ETT/ARF3). We previously reported that six CpG dinucleotides in its exon 6 are thoroughly methylated by METHYLTRASFERASE1, that CpG methylation levels are inversely correlated with ETT/ARF3 transcript levels and that methylation levels at three out of the six CpG dinucleotides are decreased in as2-1. All these imply that AS2 is involved in epigenetic repression of ETT/ARF3 by gene body DNA methylation. The mechanism of the epigenetic repression by AS2, however, is unknown. Here, we tested mutations of NUCLEOLIN1 (NUC1) and RNA HELICASE10 (RH10) encoding nucleolus-localized proteins for the methylation in exon 6 as these mutations enhance the level of ETT/ARF3 transcripts in as2-1. Methylation levels at three specific CpGs were decreased in rh10-1, and two of those three overlapped with those in as2-1. Methylation levels at two specific CpGs were decreased in nuc1-1, and one of those three overlapped with that in as2-1. No site was affected by both rh10-1 and nuc1-1. One specific CpG was unaffected by these mutations. These results imply that the way in which RH10, NUC1 and AS2 are involved in maintaining methylation at five CpGs in exon 6 might be through at least several independent pathways, which might interact with each other. Furthermore, we found that AS2 binds specifically the sequence containing CpGs in exon 1 of ETT/ARF3, and that the binding requires the zinc-finger-like motif in AS2 that is structurally similar to the zinc finger-CxxC domain in vertebrate DNA methyltransferase1.

Enhancing gene regulatory network inference through data integration with markov random fields

A gene regulatory network links transcription factors to their target genes and represents a map of transcriptional regulation. Much progress has been made in deciphering gene regulatory networks computationally. However, gene regulatory network inference for most eukaryotic organisms remain challenging. To improve the accuracy of gene regulatory network inference and facilitate candidate selection for experimentation, we developed an algorithm called GRACE (Gene Regulatory network inference ACcuracy Enhancement). GRACE exploits biological a priori and heterogeneous data integration to generate high- confidence network predictions for eukaryotic organisms using Markov Random Fields in a semi-supervised fashion. GRACE uses a novel optimization scheme to integrate regulatory evidence and biological relevance. It is particularly suited for model learning with sparse regulatory gold standard data. We show GRACE’s potential to produce high confidence regulatory networks compared to state of the art approaches using Drosophila melanogaster and Arabidopsis thaliana data. In an A. thaliana developmental gene regulatory network, GRACE recovers cell cycle related regulatory mechanisms and further hypothesizes several novel regulatory links, including a putative control mechanism of vascular structure formation due to modifications in cell proliferation.

The founder-cell transcriptome in the Arabidopsis apetala1 cauliflower inflorescence meristem

BACKGROUND

Although the pattern of lateral organ formation from apical meristems establishes species-specific plant architecture, the positional information that confers cell fate to cells as they transit to the meristem flanks where they differentiate, remains largely unknown. We have combined fluorescence-activated cell sorting and RNA-seq to characterise the cell-type-specific transcriptome at the earliest developmental time-point of lateral organ formation using DORNRÖSCHEN-LIKE::GFP to mark founder-cell populations at the periphery of the inflorescence meristem (IM) in apetala1 cauliflower double mutants, which overproliferate IMs.

RESULTS

Within the lateral organ founder-cell population at the inflorescence meristem, floral primordium identity genes are upregulated and stem-cell identity markers are downregulated. Additional differentially expressed transcripts are involved in polarity generation and boundary formation, and in epigenetic and post-translational changes. However, only subtle transcriptional reprogramming within the global auxin network was observed.

CONCLUSIONS

The transcriptional network of differentially expressed genes supports the hypothesis that lateral organ founder-cell specification involves the creation of polarity from the centre to the periphery of the IM and the establishment of a boundary from surrounding cells, consistent with bract initiation. However, contrary to the established paradigm that sites of auxin response maxima pre-pattern lateral organ initiation in the IM, auxin response might play a minor role in the earliest stages of lateral floral initiation.

Mobile small RNAs regulate genome-wide DNA methylation

RNA silencing at the transcriptional and posttranscriptional levels regulates endogenous gene expression, controls invading transposable elements (TEs), and protects the cell against viruses. Key components of the mechanism are small RNAs (sRNAs) of 21-24 nt that guide the silencing machinery to their nucleic acid targets in a nucleotide sequence-specific manner. Transcriptional gene silencing is associated with 24-nt sRNAs and RNA-directed DNA methylation (RdDM) at cytosine residues in three DNA sequence contexts (CG, CHG, and CHH). We previously demonstrated that 24-nt sRNAs are mobile from shoot to root in Arabidopsis thaliana and confirmed that they mediate DNA methylation at three sites in recipient cells. In this study, we extend this finding by demonstrating that RdDM of thousands of loci in root tissues is dependent upon mobile sRNAs from the shoot and that mobile sRNA-dependent DNA methylation occurs predominantly in non-CG contexts. Mobile sRNA-dependent non-CG methylation is largely dependent on the DOMAINS REARRANGED METHYLTRANSFERASES 1/2 (DRM1/DRM2) RdDM pathway but is independent of the CHROMOMETHYLASE (CMT)2/3 DNA methyltransferases. Specific superfamilies of TEs, including those typically found in gene-rich euchromatic regions, lose DNA methylation in a mutant lacking 22- to 24-nt sRNAs (dicer-like 2, 3, 4 triple mutant). Transcriptome analyses identified a small number of genes whose expression in roots is associated with mobile sRNAs and connected to DNA methylation directly or indirectly. Finally, we demonstrate that sRNAs from shoots of one accession move across a graft union and target DNA methylation de novo at normally unmethylated sites in the genomes of root cells from a different accession.