CHH Methylation Research Articles

ZmKTF1 promotes salt tolerance by mediating RNA-directed DNA methylation in maize.

The epigenetic process of RNA-directed DNA methylation (RdDM) regulates the expression of genes and transposons. However, little is known about the involvement of RdDM in the response of maize (Zea mays) to salt stress. Here, we isolated a salt-sensitive maize mutant and cloned the underlying gene, which encodes KOW DOMAIN-CONTAINING TRANSCRIPTION FACTOR1 (KTF1), an essential component of the RdDM pathway. Evolutionary analysis identified two homologs of KTF1 (ZmKTF1A and ZmKTF1B) with highly similar expression patterns. Whole-genome bisulfite sequencing revealed that mutations in ZmKTF1 substantially decrease genome-wide CHH (H = A, C, or T) methylation levels. Moreover, our findings suggest that ZmKTF1-mediated DNA methylation regulates the expression of multiple key genes involved in oxidoreductase activity upon exposure to salt, concomitant with increased levels of reactive oxygen species. In addition, insertion-deletion mutations (InDels) in the promoter of ZmKTF1 affect its expression, thereby altering Na+ concentrations in seedlings in a natural maize population. Therefore, ZmKTF1 might represent an untapped epigenetic resource for improving salt tolerance in maize. Overall, our work demonstrates the critical role of ZmKTF1 involved in the RdDM pathway in maize salt tolerance.

Single-Base Methylome Analysis of Sweet Cherry (Prunus avium L.) on Dwarfing Rootstocks Reveals Epigenomic Differences Associated with Scion Dwarfing Conferred by Grafting

Plant grafting using dwarfing rootstocks is one of the important cultivation measures in the sweet cherry (Prunus avium) industry. In this work, we aimed to explore the effects of the dwarfing rootstock “Pd1” (Prunus tomentosa) on sweet cherry ‘Shuguang2’ scions by performing morphological observations using the paraffin slice technique, detecting GA (gibberellin) and IAA (auxin) contents using UPLC-QTRAP-MS (ultra-performance liquid chromatography coupled with a hybrid triple quadrupole-linear ion trap mass spectrometer), and implementing integration analyses of the epigenome and transcriptome using whole-genome bisulfite sequencing and transcriptome sequencing. Anatomical analysis indicated that the cell division ability of the SAM (shoot apical meristem) in dwarfing plants was reduced. Pd1 rootstock significantly decreased the levels of GAs and IAA in sweet cherry scions. Methylome analysis showed that the sweet cherry genome presented 15.2~18.6%, 59.88~61.55%, 28.09~33.78%, and 2.99~5.28% methylation at total C, CG, CHG, and CHH sites, respectively. Shoot tips from dwarfing plants exhibited a hypermethylated pattern mostly due to increased CHH methylation, while leaves exhibited a hypomethylated pattern. According to GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, DMGs (differentially methylated genes) and DEGs (differentially expressed genes) were enriched in hormone-related GO terms and KEGG pathways. Global correlation analysis between methylation and transcription revealed that mCpG in the gene body region enhanced gene expression and mCHH in the region near the TSS (transcription start site) was positively correlated with gene expression. Next, we found some hormone-related genes and TFs with significant changes in methylation and transcription, including SAURs, ARF, GA2ox, ABS1, bZIP, MYB, and NAC. This study presents a methylome map of the sweet cherry genome, revealed widespread DNA methylation alterations in scions caused by dwarfing rootstock, and obtained abundant genes with methylation and transcription alterations that are potentially involved in rootstock-induced growth changes in sweet cherry scions. Our findings can lay a good basis for further epigenetic studies on sweet cherry dwarfing and provide valuable new insight into understanding rootstock–scion interactions.

Divergent Flowering Time Responses to Increasing Temperatures Are Associated With Transcriptome Plasticity and Epigenetic Modification Differences at FLC Promoter Region of Arabidopsis thaliana.

Understanding the genetic, and transcriptomic changes that drive the phenotypic plasticity of fitness traits is a central question in evolutionary biology. In this study, we utilised 152 natural Swedish Arabidopsis thaliana accessions with re-sequenced genomes, transcriptomes and methylomes and measured flowering times (FTs) under two temperature conditions (10°C and 16°C) to address this question. We revealed that the northern accessions exhibited advanced flowering in response to decreased temperature, whereas the southern accessions delayed their flowering, indicating a divergent flowering response. This contrast in flowering responses was associated with the isothermality of their native ranges, which potentially enables the northern accessions to complete their life cycle more rapidly in years with shorter growth seasons. At the transcriptome level, we observed extensive rewiring of gene co-expression networks, with the expression of 25 core genes being associated with the mean FT and its plastic variation. Notably, variations in FLC expression sensitivity between northern and southern accessions were found to be associated with the divergence FT response. Further analysis suggests that FLC expression sensitivity is associated with differences in CG, CHG and CHH methylation at the promoter region. Overall, our study revealed the association between transcriptome plasticity and flowering time plasticity among different accessions, providing evidence for its relevance in ecological adaptation. These findings offer deeper insights into the genetics of rapid responses to environmental changes and ecological adaptation.

Functional Characterization of Accessible Chromatin in Common Wheat.

Eukaryotic gene transcription is fine-tuned by precise spatiotemporal interactions between cis-regulatory elements (CREs) and trans-acting factors. However, how CREs individually or coordinated with epigenetic marks function in regulating homoeolog bias expression is still largely unknown in wheat. In this study, through comprehensively characterizing open chromatin coupled with DNA methylation in the seedling and spikelet of common wheat, we observed that differential chromatin openness occurred between the seedling and spikelet, which plays important roles in tissue development through regulating the expression of related genes or through the transcription factor (TF)-centered regulatory network. Moreover, we found that CHH methylation may act as a key determinant affecting the differential binding of TFs, thereby resulting in differential expression of target genes. In addition, we found that sequence variations in MNase hypersensitive sites (MHSs) result in the differential expression of key genes responsible for important agronomic traits. Thus, our study provides new insights into the roles of CREs in regulating tissue or homoeolog bias expression, and controlling important agronomic traits in common wheat. It also provides potential CREs for genetic and epigenetic manipulation toward improving desirable traits for wheat molecule breeding.

Cashmere cyclic growth affected by different photoperiods alters DNA methylation patterns

Whole-genome bisulfite sequencing was used to study epigenetic differences in hair follicles from Inner Mongolia Albas cashmere goats during long and short photoperiods. Results showed earlier anagen entry in short photoperiod. Genome-wide cytosine methylation was 3.23%, with CG, CHG, and CHH methylation rates at 48.00%, 0.19%, and 0.17%, respectively. Methylation was higher for CG than for CHG and CHH, and higher for introns than exons and downstream regions. Both groups had similar levels of CG, CHG, and CHH methylation, but CHH methylation levels differed slightly between photoperiods. A total of 9,700 differentially methylated regions (DMRs) and 4,187 genes (DMGs) were identified. Gene Ontology (GO) and KEGG analyses revealed significant enrichment in 67 GO terms and 41 signaling pathways related to circadian rhythms, bringing new insights in the role of DNA methylation in cashmere growth and development.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana.

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Substrate specificity and protein stability drive the divergence of plant-specific DNA methyltransferases.

DNA methylation is an important epigenetic mechanism essential for transposon silencing and genome integrity. Across evolution, the substrates of DNA methylation have diversified between kingdoms to account for genome complexity. In plants, Chromomethylase3 (CMT3) and CMT2 are the major methyltransferases mediating CHG and CHH methylation, respectively. However, how these two enzymes diverge on substrate specificities during evolution remains unknown. Here, we reveal that CMT2 originates from a duplication of the evolutionarily more ancient CMT3 in flowering plants. Lacking a key arginine residue recognizing CHG in CMT2 impairs its CHG methylation activity in most flowering plants. An engineered V1200R mutation empowers CMT2 to restore both CHG and CHH methylation in Arabidopsis cmt2cmt3 mutant, testifying a loss-of-function effect for CMT2 after ∼200 million years of evolution. Interestingly, CMT2 has evolved a long and unstructured N-terminus critical for balancing protein stability, especially under heat stress. Furthermore, CMT2 N-terminus is plastic and can be tolerant to various natural mutations. Together, this study reveals the mechanism of chromomethylase divergence for context-specific DNA methylation in plants and sheds important lights on DNA methylation evolution and function.

Transcriptomic and epigenetic responses shed light on soybean resistance to Phytophthora sansomeana.

Phytophthora root rot, caused by oomycete pathogens in the Phytophthora genus, poses a significant threat to soybean productivity. While resistance mechanisms against Phytophthora sojae have been extensively studied in soybean, the molecular basis underlying immune responses to Phytophthora sansomeana remains unclear. In this study, we investigated transcriptomic and epigenetic responses of two resistant (Colfax and NE2701) and two susceptible (Williams 82 and Senaki) soybean lines at four time points (2, 4, 8, and 16h post inoculation [hpi]) after P. sansomeana inoculation. Comparative transcriptomic analyses revealed a greater number of differentially expressed genes (DEGs) upon pathogen inoculation in resistant lines, particularly at 8 and 16 hpi. These DEGs were predominantly associated with defense response, ethylene, and reactive oxygen species-mediated defense pathways. Moreover, DE transposons were predominantly upregulated after inoculation, and more of them were enriched near genes in Colfax than other soybean lines. Notably, we identified a long non-coding RNA (lncRNA) within the mapped region of the resistance gene that exhibited exclusive upregulation in the resistant lines after inoculation, potentially regulating two flanking LURP-one-related genes. Furthermore, DNA methylation analysis revealed increased CHH (where H=A, T, or C) methylation levels in lncRNAs after inoculation, with delayed responses in Colfax compared to Williams 82. Overall, our results provide comprehensive insights into soybean responses to P. sansomeana, highlighting potential roles of lncRNAs and epigenetic regulation in plant defense.

Genome-wide chromatin accessibility reveals transcriptional regulation of heterosis in inter-subspecific hybrid rice.

The utilization of rice heterosis is essential for ensuring global food security; however, its molecular mechanism remains unclear. In this study, comprehensive analyses of accessible chromatin regions (ACRs), DNA methylation, and gene expression in inter-subspecific hybrid and its parents were performed to determine the potential role of chromatin accessibility in rice heterosis. The hybrid exhibited abundant ACRs, in which the gene ACRs and proximal ACRs were directly related to transcriptional activation rather than the distal ACRs. Regarding the dynamic accessibility contribution of the parents, paternal ZHF1015 transmitted a greater number of ACRs to the hybrid. Accessible genotype-specific target genes were enriched with overrepresented transcription factors, indicating a unique regulatory network of genes in the hybrid. Compared with its parents, the differentially accessible chromatin regions with upregulated chromatin accessibility were much greater than those with downregulated chromatin accessibility, reflecting a stronger regulation in the hybrid. Furthermore, DNA methylation levels were negatively correlated with ACR intensity, and genes were strongly affected by CHH methylation in the hybrid. Chromatin accessibility positively regulated the overall expression level of each genotype. ACR-related genes with maternal Z04A-bias allele-specific expression tended to be enriched during carotenoid biosynthesis, whereas paternal ZHF1015-bias genes were more active in carbohydrate metabolism. Our findings provide a new perspective on the mechanism of heterosis based on chromatin accessibility in inter-subspecific hybrid rice.

Difference in DNA methylation pattern and expression of active ingredients between garden ginseng and ginseng under forest

Panax ginseng is widely used for its medicinal value. With the decline in wild ginseng resources, garden ginseng (GG) and ginseng under forest (FG) have become the main sources. Their different growth conditions and ages may account for the differences in epigenetics. The genomic DNA of Panax ginseng was digested using the methylation-insensitive enzyme MseⅠ combined with the methylation-sensitive enzymes AciⅠ, PstⅠ, and EcoT22Ⅰ. Then, a sequencing library of CG, CHG, and CHH methylation contexts was constructed and sequenced through amplicon sequencing. Differences in DNA methylation were observed using the methylation context-sensitive enzyme ddRAD (MCSeEd), revealing eight differentially methylated regions (DMRs). The bisulfite sequencing technology revealed that the methylation level of GG was higher than that of FG. The Plant ISOform Sequencing Database was used to predict the regulatory function of CpG island–coding region so as to estimate the relationship between methylation and effective components. The quantitative real-time polymerase chain reaction (qPCR) showed that the expression of functional genes in FG was higher than that in GG. These studies showed that the increase in GG methylation level was associated with the decrease in gene expression, providing a new perspective on epigenetic differences and gene expression in ginseng.

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.

Reduced CHH methylation levels reveal a critical role of aging pathway genes in Moso bamboo flowering

Moso bamboo holds significant economic importance in China, serving various purposes, such as food, material, ornamentation, and greenery. Despite its versatility, the occurrence of flowering in Moso bamboo poses a threat to bamboo forests, resulting in substantial losses. The underlying cause of bamboo flowering remains elusive. Dynamic fluctuations in DNA methylation govern the transcriptional levels of crucial genes pivotal for plant growth and development. In this study, we conducted comprehensive DNA methylation (by whole-genome bisulfite sequencing) and transcriptome (by RNA-seq) analyses on non-flowering leaves, flowering leaves, and spikelets of Moso bamboo. Our findings revealed a notable reduction in the overall DNA methylation level, particularly CHH methylation, from leaves to spikelets, influencing the expression of differentially regulated genes. Notably, we identified DNA methylation as a regulatory mechanism for numerous flowering-related genes, including SPLs, FT, and SOC1. Specifically, the SPL3f gene, a key regulator of the aging pathway, exhibited hypomethylation and a high expression level in spikelets. Conversely, SOC1c displayed transcriptional silencing attributed to hypermethylation in the CHH context in the leaves of non-flowering plants. DNA methylation may affect the flowering mechanism of Moso bamboo by regulating the expression of key genes. In summary, our results shed light on the dynamic changes in DNA methylation between leaves and spikelets, unraveling an important epigenetic modification mechanism for flowering in Moso bamboo.

Abstract 7030: Reference materials for analysis of DNA methylation in cell-free circulating tumor DNA

Abstract Aberrant DNA methylation is associated with many cancers including breast, liver, bone, and colon and has been identified as a biomarker for early cancer screening. Liquid biopsies are beginning to screen for DNA methylation in cancer-derived DNA, but accurate assessment of methylation is not trivial, and methods like bisulfite sequencing can result in biased data. Reference materials for establishing the analytical validity of measurements of CpG methylation are not widely available. In order to prepare such reference materials, genomic DNA (gDNA) was extracted from the well characterized human reference cell line GM24385 and fragmented to a circulating tumor DNA (ctDNA)-like size distribution. This DNA was processed further to prepare CpG-methylated and unmethylated versions. The size distribution was determined by using the Bioanalyzer High Sensitivity DNA Kit (Agilent, Santa Clara, CA). Methylation status in the reference materials was confirmed by both digital PCR (dPCR) and NGS based methods. A digital PCR method was developed to analyze heterozygous single nucleotide polymorphisms (SNPs) using the Bio-Rad QX200™ Droplet Digital PCR (ddPCR) System. To determine methylation levels by NGS, the NEBNext™ Enzymatic Methyl-seq Kit was used (EM-seq). Paired-end sequencing of the libraries was performed on a NextSeq™ 2000 (Illumina, San Diego, CA). Data alignment to hg19 was done with the bwa-meth (Pedersen, BS., et al.) aligner, and assessment of methylation in the aligned data was done with MethylDackel (Ryan, D.). Each reference material had an average peak size of ~160 bp in length, similar to patient cfDNA. By dPCR, the unmethylated reference material had an average of 0.73% DNA methylation, while the methylated reference material had an average of 100.10% DNA methylation. Values can be above 100% by dPCR due to unequal representation of both alleles. By NGS, the unmethylated reference material had an average of 0.93 % CpG methylation, which was similar to the background of 1.13 % and 1.31 % for CHG and CHH methylation, respectively. The methylated reference material had an average of 93.95 % CpG methylation, with a background of 0.69 % and 0.74 % for CHG and CHH methylation, respectively. Because EM-seq relies on enzymatic protection and conversion steps, it may be that these reactions did not allow for full protection and full conversion, which led to values above 0 % for unmethylated and below 100 % for methylated materials. We have developed reference materials for the assessment of DNA methylation in ctDNA. One unmethylated and one methylated reference material was created. These reference materials can be mixed together to the level of methylation required to help aid in assay design, optimization, and validation. Citation Format: Alysha Higgs, Angelique Ryan, Jayanthi Ramprakash, Dana Ruminski Lowe, Yves Konigshofer, Catherine Huang, Andrew Anfora, Russell Garlick, Bharathi Anekella. Reference materials for analysis of DNA methylation in cell-free circulating tumor DNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7030.

Comparative and functional analysis unveils the contribution of photoperiod to DNA methylation, sRNA accumulation, and gene expression variations in short-day and long-day grasses.

Photoperiod employs complicated networks to regulate various developmental processes in plants, including flowering transition. However, the specific mechanisms by which photoperiod affects epigenetic modifications and gene expression variations in plants remain elusive. In this study, we conducted a comprehensive analysis of DNA methylation, small RNA (sRNA) accumulation, and gene expressions under different daylengths in facultative long-day (LD) grass Brachypodium distachyon and short-day (SD) grass rice. Our results showed that while overall DNA methylation levels were minimally affected by different photoperiods, CHH methylation levels were repressed under their favorable light conditions, particularly in rice. We identified numerous differentially methylated regions (DMRs) that were influenced by photoperiod in both plant species. Apart from differential sRNA clusters, we observed alterations in the expression of key components of the RNA-directed DNA methylation pathway, DNA methyltransferases, and demethylases, which may contribute to the identified photoperiod-influenced CHH DMRs. Furthermore, we identified many differentially expressed genes in response to different daylengths, some of which were associated with the DMRs. Notably, we discovered a photoperiod-responsive gene MYB11 in the transcriptome of B.distachyon, and further demonstrated its role as a flowering inhibitor by repressing FT1 transcription. Together, our comparative and functional analysis sheds light on the effects of daylength on DNA methylation, sRNA accumulation, and gene expression variations in LD and SD plants, thereby facilitating better designing breeding programs aimed at developing high-yield crops that can adapt to local growing seasons.

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer.

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Dynamic DNA Methylation Regulates Season-Dependent Secondary Metabolism in the New Shoots of Tea Plants.

Plant secondary metabolites are critical quality-conferring compositions of plant-derived beverages, medicines, and industrial materials. The accumulations of secondary metabolites are highly variable among seasons; however, the underlying regulatory mechanism remains unclear, especially in epigenetic regulation. Here, we used tea plants to explore an important epigenetic mark DNA methylation (5mC)-mediated regulation of plant secondary metabolism in different seasons. Multiple omics analyses were performed on spring and summer new shoots. The results showed that flavonoids and theanine metabolism dominated in the metabolic response to seasons in the new shoots. In summer new shoots, the genes encoding DNA methyltransferases and demethylases were up-regulated, and the global CG and CHG methylation reduced and CHH methylation increased. 5mC methylation in promoter and gene body regions influenced the seasonal response of gene expression; the amplitude of 5mC methylation was highly correlated with that of gene transcriptions. These differentially methylated genes included those encoding enzymes and transcription factors which play important roles in flavonoid and theanine metabolic pathways. The regulatory role of 5mC methylation was further verified by applying a DNA methylation inhibitor. These findings highlight that dynamic DNA methylation plays an important role in seasonal-dependent secondary metabolism and provide new insights for improving tea quality.

The lemon genome and DNA methylome unveil epigenetic regulation of citric acid biosynthesis during fruit development.

Citric acid gives lemons their unique flavor, which impacts their sensory traits and market value. However, the intricate process of citric acid accumulation during lemon fruit growth remains incompletely understood. Here, we achieved a chromosomal-level genome assembly for the 'Xiangshui' lemon variety, spanning 364.85Mb across nine chromosomes. This assembly revealed 27 945 genes and 51.37% repetitive sequences, tracing the divergence from citron 2.85 million years ago. DNA methylome analysis of lemon fruits across different developmental stages revealed significant variations in DNA methylation. We observed decreased CG and CHG methylation but increased CHH methylation. Notably, the expression of RdDM pathway-related genes increased with fruit development, suggesting a connection with elevated CHH methylation, which is potentially influenced by the canonical RdDM pathway. Furthermore, we observed that elevated CHH DNA methylation within promoters significantly influenced the expression of key genes, critically contributing to vital biological processes, such as citric acid accumulation. In particular, the pivotal gene phosphoenolpyruvate carboxykinase (ClPEPCK), which regulates the tricarboxylic acid cycle, was strikingly upregulated during fruit development, concomitant with increased CHH methylation in its promoter region. Other essential genes associated with citric acid accumulation, such as the MYB transcription factor (ClPH1/4/5) and ANTHOCYANIN 1 (ClAN1), were strongly correlated with DNA methylation levels. These results strongly indicate that DNA methylation crucially orchestrates the metabolic synthesis of citric acid. In conclusion, our study revealed dynamic changes in DNA methylation during lemon fruit development, underscoring the significant role of DNA methylation in controlling the citric acid metabolic pathway.

Heritable changes of epialleles near genes in maize can be triggered in the absence of CHH methylation.

Trans-chromosomal interactions resulting in changes in DNA methylation during hybridization have been observed in several plant species. However, little is known about the causes or consequences of these interactions. Here, we compared DNA methylomes of F1 hybrids that are mutant for a small RNA biogenesis gene, Mop1 (Mediator of paramutation1), with that of their parents, wild-type siblings, and backcrossed progeny in maize (Zea mays). Our data show that hybridization triggers global changes in both trans-chromosomal methylation (TCM) and trans-chromosomal demethylation (TCdM), most of which involved changes in CHH methylation. In more than 60% of these TCM differentially methylated regions (DMRs) in which small RNAs are available, no significant changes in the quantity of small RNAs were observed. Methylation at the CHH TCM DMRs was largely lost in the mop1 mutant, although the effects of this mutant varied depending on the location of these DMRs. Interestingly, an increase in CHH at TCM DMRs was associated with enhanced expression of a subset of highly expressed genes and suppressed expression of a small number of lowly expressed genes. Examination of the methylation levels in backcrossed plants demonstrates that both TCM and TCdM can be maintained in the subsequent generation, but that TCdM is more stable than TCM. Surprisingly, although increased CHH methylation in most TCM DMRs in F1 plants required Mop1, initiation of a new epigenetic state of these DMRs did not require a functional copy of this gene, suggesting that initiation of these changes is independent of RNA-directed DNA methylation.

RNA m6A modification facilitates DNA methylation during maize kernel development.

N6-methyladenosine (m6A) in mRNA and 5-methylcytosine (5mC) in DNA have critical functions for regulating gene expression and modulating plant growth and development. However, the interplay between m6A and 5mC is an elusive territory and remains unclear mechanistically in plants. We reported an occurrence of crosstalk between m6A and 5mC in maize (Zea mays) via the interaction between mRNA adenosine methylase (ZmMTA), the core component of the m6A methyltransferase complex, and decrease in DNA methylation 1 (ZmDDM1), a key chromatin-remodeling factor that regulates DNA methylation. Genes with m6A modification were coordinated with a much higher level of DNA methylation than genes without m6A modification. Dysfunction of ZmMTA caused severe arrest during maize embryogenesis and endosperm development, leading to a significant decrease in CHH methylation in the 5' region of m6A-modified genes. Instead, loss of function of ZmDDM1 had no noteworthy effects on ZmMTA-related activity. This study establishes a direct link between m6A and 5mC during maize kernel development and provides insights into the interplay between RNA modification and DNA methylation.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context, but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.