Efficient and accurate determination of genome-wide DNA methylation patterns in Arabidopsis thaliana with enzymatic methyl sequencing

Whole Genome Methylation Sequencing via Enzymatic Conversion (EM-seq): Protocol, Data Processing, and Analysis

Whole genome bisulfite sequencing (WGBS) has been the gold standard technique for base resolution analysis of DNA methylation for the last 15 years. It has been, however, associated with technical biases, which lead to overall overestimation of global and regional methylation values, and significant artifacts in extreme cytosine-rich DNA sequence contexts. Enzymatic conversion of cytosine is the newest approach, set to replace entirely the use of the damaging bisulfite conversion of DNA. The EM-seq technique utilizes TET2, T4-BGT, and APOBEC in a two-step conversion process, where the modified cytosines are first protected by oxidation and glucosylation, followed by deamination of all unmodified cytosines to uracil. As a result, EM-seq is degradation-free and bias-free, requires low DNA input, and produces high library yields with longer reads, little batch variation, less duplication, uniform genomic coverage, accurate methylation over a larger number of captured CpGs, and no sequence-specific artifacts.

Epigenetic gene regulation in plants and its potential applications in crop improvement

DNA methylation, also known as 5-methylcytosine, is an epigenetic modification that has crucial functions in plant growth, development and adaptation. The cellular DNA methylation level is tightly regulated by the combined action of DNA methyltransferases and demethylases. Protein complexes involved in the targeting and interpretation of DNA methylation have been identified, revealing intriguing roles of methyl-DNA binding proteins and molecular chaperones. Structural studies and in vitro reconstituted enzymatic systems have provided mechanistic insights into RNA-directed DNA methylation, the main pathway catalysing de novo methylation in plants. A better understanding of the regulatory mechanisms will enable locus-specific manipulation of the DNA methylation status. CRISPR-dCas9-based epigenome editing tools are being developed for this goal. Given that DNA methylation patterns can be stably transmitted through meiosis, and that large phenotypic variations can be contributed by epimutations, epigenome editing holds great promise in crop breeding by creating additional phenotypic variability on the same genetic material.

PTGS is dispensable for the initiation of epigenetic silencing of an active transposon in Arabidopsis

Transposable elements (TEs) are repressed in plants through transcriptional gene silencing (TGS), maintained epigenetic silencing marks such as DNA methylation. However, the mechanisms by which silencing is first installed remain poorly understood in plants. Small interfering (si)RNAs and post-transcriptional gene silencing (PTGS) are believed to mediate the initiation of TGS by guiding the first deposition of DNA methylation. To determine how this silencing installation works, we took advantage of ÉVADÉ (EVD), an endogenous retroelement in Arabidopsis, able to recapitulate true de novo silencing with a sequence of PTGS followed by a TGS. To test whether PTGS is required for TGS, we introduce active EVD into RNA-DEPENDENT-RNA-POLYMERASE-6 (RDR6) mutants, an essential PTGS component. EVD activity and silencing are monitored across several generations. In the absence of PTGS, silencing of EVD is still achieved through installation of RNA-directed DNA methylation (RdDM). Our study shows that PTGS is dispensable for de novo EVD silencing. Although we cannot rule out that PTGS might facilitate TGS, or control TE activity, initiation of epigenetic silencing can take place in its absence.

Perfluorooctanesulfonic acid (PFOS) induced cancer related DNA methylation alterations in human breast cells: A whole genome methylome study

DNA methylation plays a pivotal role in cancer. The ubiquitous contaminant perfluorooctanesulfonic acid (PFOS) has been epidemiologically associated with breast cancer, and can induce proliferation and malignant transformation of normal human breast epithelial cells (MCF-10A), but the information about its effect on DNA methylation is sparse. The aim of this study was to characterize the whole-genome methylome effects of PFOS in our breast cell model and compare the findings with previously demonstrated DNA methylation alterations in breast tumor tissues. The DNA methylation profile was assessed at single CpG resolution in MCF-10A cells treated with 1 μM PFOS for 72 h by using Enzymatic Methyl sequencing (EM-seq). We found 12,591 differentially methylated CpG-sites and 13,360 differentially methylated 100 bp tiles in the PFOS exposed breast cells. These differentially methylated regions (DMRs) overlapped with 2406 genes of which 494 were long non-coding RNA and 1841 protein coding genes. We identified 339 affected genes that have been shown to display altered DNA methylation in breast cancer tissue and several other genes related to cancer development. This includes hypermethylation of GACAT3, DELEC1, CASC2, LCIIAR, MUC16, SYNE1 and hypomethylation of TTN and KMT2C. DMRs were also found in estrogen receptor genes (ESR1, ESR2, ESRRG, ESRRB, GREB1) and estrogen responsive genes (GPER1, EEIG1, RERG). The gene ontology analysis revealed pathways related to cancer phenotypes such as cell adhesion and growth. These findings improve the understanding of PFOS's potential role in breast cancer and illustrate the value of whole-genome methylome analysis in uncovering mechanisms of chemical effects, identifying biomarker candidates, and strengthening epidemiological associations, potentially impacting risk assessment.

Improved detection of methylation in ancient DNA

Reconstructing premortem DNA methylation levels in ancient DNA has led to breakthrough studies such as the prediction of anatomical features of the Denisovan. These studies rely on computationally inferring methylation levels from damage signals in naturally deaminated cytosines, which requires expensive high-coverage genomes. Here, we test two methods for direct methylation measurement developed for modern DNA based on either bisulfite or enzymatic methylation treatments. Bisulfite treatment shows the least reduction in DNA yields as well as the least biases during methylation conversion, demonstrating that this method can be successfully applied to ancient DNA.

Cost-effective solutions for high-throughput enzymatic DNA methylation sequencing

Characterizing DNA methylation patterns is important for addressing key questions in evolutionary biology, geroscience, and medical genomics. While costs are decreasing, whole-genome DNA methylation profiling remains prohibitively expensive for most population-scale studies, creating a need for cost-effective, reduced representation approaches (i.e., assays that rely on microarrays, enzyme digests, or sequence capture to target a subset of the genome). Most common whole genome and reduced representation techniques rely on bisulfite conversion, which can damage DNA resulting in DNA loss and sequencing biases. Enzymatic methyl sequencing (EM-seq) was recently proposed to overcome these issues, but thorough benchmarking of EM-seq combined with cost-effective, reduced representation strategies has not yet been performed. To do so, we optimized Targeted Methylation Sequencing protocol (TMS)-which profiles ∼4 million CpG sites-for miniaturization, flexibility, and multispecies use at a cost of ∼$80. First, we tested modifications to increase throughput and reduce cost, including increasing multiplexing, decreasing DNA input, and using enzymatic rather than mechanical fragmentation to prepare DNA. Second, we compared our optimized TMS protocol to commonly used techniques, specifically the Infinium MethylationEPIC BeadChip (n=55 paired samples) and whole genome bisulfite sequencing (n=6 paired samples). In both cases, we found strong agreement between technologies (R² = 0.97 and 0.99, respectively). Third, we tested the optimized TMS protocol in three non-human primate species (rhesus macaques, geladas, and capuchins). We captured a high percentage (mean=77.1%) of targeted CpG sites and produced methylation level estimates that agreed with those generated from reduced representation bisulfite sequencing (R² = 0.98). Finally, we applied our protocol to profile age-associated DNA methylation variation in two subsistence-level populations-the Tsimane of lowland Bolivia and the Orang Asli of Peninsular Malaysia-and found age-methylation patterns that were strikingly similar to those reported in high income cohorts, despite known differences in age-health relationships between lifestyle contexts. Altogether, our optimized TMS protocol will enable cost-effective, population-scale studies of genome-wide DNA methylation levels across human and non-human primate species.

Comparing methylation levels assayed in GC-rich regions with current and emerging methods

DNA methylation is an epigenetic mechanism that regulates gene expression, and for mammals typically occurs on cytosines within CpG dinucleotides. A significant challenge for methylation detection methods is accurately measuring methylation levels within GC-rich regions such as gene promoters, as inaccuracies compromise downstream biological interpretation of the data. To address this challenge, we compared methylation levels assayed using four different Methods Enzymatic Methyl-seq (EM-seq), whole genome bisulphite sequencing (WGBS), Infinium arrays (Illumina MethylationEPIC, "EPIC"), and Oxford Nanopore Technologies nanopore sequencing (ONT) applied to human DNA. Overall, all methods produced comparable and consistent methylation readouts across the human genome. The flexibility offered by current gold standard WGBS in interrogating genome-wide cytosines is surpassed technically by both EM-seq and ONT, as their coverages and methylation readouts are less prone to GC bias. These advantages are tempered by increased laboratory time (EM-seq) and higher complexity (ONT). We further assess the strengths and weaknesses of each method, and provide recommendations in choosing the most appropriate methylation method for specific scientific questions or translational needs.

DNA methylation of exercise-responsive genes differs between trained and untrained men

BACKGROUND

Physical activity is well known for its multiple health benefits and although the knowledge of the underlying molecular mechanisms is increasing, our understanding of the role of epigenetics in long-term training adaptation remains incomplete. In this intervention study, we included individuals with a history of > 15 years of regular endurance or resistance training compared to age-matched untrained controls performing endurance or resistance exercise. We examined skeletal muscle DNA methylation of genes involved in key adaptation processes, including myogenesis, gene regulation, angiogenesis and metabolism.

RESULTS

A greater number of differentially methylated regions and differentially expressed genes were identified when comparing the endurance group with the control group than in the comparison between the strength group and the control group at baseline. Although the cellular composition of skeletal muscle samples was generally consistent across groups, variations were observed in the distribution of muscle fiber types. Slow-twitch fiber type genes MYH7 and MYL3 exhibited lower promoter methylation and elevated expression in endurance-trained athletes, while the same group showed higher methylation in transcription factors such as FOXO3, CREB5, and PGC-1α. The baseline DNA methylation state of those genes was associated with the transcriptional response to an acute bout of exercise. Acute exercise altered very few of the investigated CpG sites.

CONCLUSIONS

Endurance- compared to resistance-trained athletes and untrained individuals demonstrated a different DNA methylation signature of selected skeletal muscle genes, which may influence transcriptional dynamics following a bout of acute exercise. Skeletal muscle fiber type distribution is associated with methylation of fiber type specific genes. Our results suggest that the baseline DNA methylation landscape in skeletal muscle influences the transcription of regulatory genes in response to an acute exercise bout.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

Increased scalability and sequencing quality of an epigenetic age prediction assay

Epigenetic ageing in a human context, has been used to better understand the relationship between age and factors such as lifestyle and genetics. In an ecological setting, it has been used to predict the age of individual animals for wildlife management. Despite the importance of epigenetic ageing in a range of research fields, the assays to measure epigenetic ageing are either expensive on a large scale or complex. In this study, we aimed to improve the efficiency and sequencing quality of an existing epigenetic ageing assay for the Australian Lungfish (Neoceratodus forsteri). We used an enzyme-based alternative to bisulfite conversion to reduce DNA fragmentation and evaluated its performance relative to bisulfite conversion. We found the sequencing quality to be 12% higher with the enzymatic alternative compared to bisulfite treatment (p-value < 0.01). This new enzymatic based approach, although currently double the cost of bisulfite treatment can increases the throughput and sequencing quality. We envisage this assay setup being adopted increasingly as the scope and scale of epigenetic ageing research continues to grow.

Single-molecule long-read methylation profiling reveals regional DNA methylation regulated by Elongator Complex Subunit 2 in Arabidopsis roots experiencing spaceflight

BACKGROUND

The Advanced Plant Experiment-04 - Epigenetic Expression (APEX-04-EpEx) experiment onboard the International Space Station examined the spaceflight-altered cytosine methylation in two genetic lines of Arabidopsis thaliana, wild-type Col-0 and the mutant elp2-5, which is deficient in an epigenetic regulator Elongator Complex Subunit 2 (ELP2). Whole-genome bisulfite sequencing (WGBS) revealed distinct spaceflight associated methylation differences, presenting the need to explore specific space-altered methylation at single-molecule resolution to associate specific changes over large regions of spaceflight related genes. To date, tools of multiplexed targeted DNA methylation sequencing remain limited for plant genomes.

RESULTS

To provide methylation data at single-molecule resolution, Flap-enabled next-generation capture (FENGC), a novel targeted multiplexed DNA capture and enrichment technique allowing cleavage at any specified sites, was applied to survey spaceflight-altered DNA methylation in genic regions of interest. The FENGC capture panel contained 108 targets ranging from 509 to 704 nt within the promoter or gene body regions of gene targets derived from spaceflight whole-genome data sets. In addition to genes with significant changes in expression and average methylation levels between spaceflight and ground control, targets with space-altered distributions of the proportion of methylated cytosines per molecule were identified. Moreover, trends of co-methylation of different cytosine contexts were exhibited in the same DNA molecules. We further identified significant DNA methylation changes in three previously biological process-unknown genes, and loss-of-function mutants of two of these genes (named as EMO1 and EMO2 for ELP2-regulated Methylation in Orbit 1 and 2) showed enhanced root growth rate.

CONCLUSIONS

FENGC simplifies and reduces the cost of multiplexed, targeted, single-molecule profiling of methylation in plants, providing additional resolution along each DNA molecule that is not seen in population-based short-read data such as WGBS. This case study has revealed spaceflight-altered regional modification of cytosine methylation occurring within single DNA molecules of cell subpopulations, which were not identified by WGBS. The single-molecule survey by FENGC can lead to identification of novel functional genes. The newly identified EMO1 and EMO2 are root growth regulators which may be epigenetically involved in plant adaptation to spaceflight.

DNA methylation drives a new path in gastric cancer early detection: Current impact and prospects

Gastric cancer (GC) is one of the most common and deadly cancers worldwide. Early detection offers the best chance for curative treatment and reducing its mortality. However, the optimal population-based early screening for GC remains unmet. Aberrant DNA methylation occurs in the early stage of GC, exhibiting cancer-specific genetic and epigenetic changes, and can be detected in the media such as blood, gastric juice, and feces, constituting a valuable biomarker for cancer early detection. Furthermore, DNA methylation is a stable epigenetic alteration, and many innovative methods have been developed to quantify it rapidly and accurately. Nonetheless, large-scale clinical validation of DNA methylation serving as tumor biomarkers is still lacking, precluding their implementation in clinical practice. In conclusion, after a critical analysis of the recent existing literature, we summarized the evolving roles of DNA methylation during GC occurrence, expounded the newly discovered noninvasive DNA methylation biomarkers for early detection of GC, and discussed its challenges and prospects in clinical applications.

A comprehensive review on Gossypium hirsutum resistance against cotton leaf curl virus

Cotton (Gossypium hirsutum L.) is a significant fiber crop. Being a major contributor to the textile industry requires continuous care and attention. Cotton is subjected to various biotic and abiotic constraints. Among these, biotic factors including cotton leaf curl virus (CLCuV) are dominant. CLCuV is a notorious disease of cotton and is acquired, carried, and transmitted by the whitefly (Bemisia tabaci). A cotton plant affected with CLCuV may show a wide range of symptoms such as yellowing of leaves, thickening of veins, upward or downward curling, formation of enations, and stunted growth. Though there are many efforts to protect the crop from CLCuV, long-term results are not yet obtained as CLCuV strains are capable of mutating and overcoming plant resistance. However, systemic-induced resistance using a gene-based approach remained effective until new virulent strains of CLCuV (like Cotton Leaf Curl Burewala Virus and others) came into existence. Disease control by biological means and the development of CLCuV-resistant cotton varieties are in progress. In this review, we first discussed in detail the evolution of cotton and CLCuV strains, the transmission mechanism of CLCuV, the genetic architecture of CLCuV vectors, and the use of pathogen and nonpathogen-based approaches to control CLCuD. Next, we delineate the uses of cutting-edge technologies like genome editing (with a special focus on CRISPR-Cas), next-generation technologies, and their application in cotton genomics and speed breeding to develop CLCuD resistant cotton germplasm in a short time. Finally, we delve into the current obstacles related to cotton genome editing and explore forthcoming pathways for enhancing precision in genome editing through the utilization of advanced genome editing technologies. These endeavors aim to enhance cotton's resilience against CLCuD.

The Effects of DNA Methylation on Cytoplasmic Male Sterility in Sugar Beet

DNA methylation is widely found in higher plants and can control gene expression by regulation without changing the DNA sequence. In this study, the whole-genome methylation map of sugar beet was constructed by WGBS (whole-genome bisulfite sequencing) technology, and the results of WGBS were verified by bisulfite transformation, indicating that the results of WGBS technology were reliable. In addition, 12 differential methylation genes (DMGs) were identified, which were related to carbohydrate and energy metabolism, pollen wall development, and endogenous hormone regulation. Quantitative real-time PCR (qRT-PCR) showed that 75% of DMG expression levels showed negative feedback with methylation level, indicating that DNA methylation can affect gene expression to a certain extent. In addition, we found hypermethylation inhibited gene expression, which laid a foundation for further study on the molecular mechanism of DNA methylation at the epigenetic level in sugar beet male sterility.

Hypermethylation and small RNA expression are associated with increased age in almond (Prunus dulcis [Mill.] D.A. Webb) accessions

The focus of this study is to profile changes in DNA methylation and small RNA expression occurring with increased age in almond breeding germplasm to identify possible biomarkers of age that can be used to assess the potential of individuals to develop aging-related disorders. To profile DNA methylation in almond germplasm, 70 methylomes were generated from almond individuals representing three age cohorts (11, 7, and 2 years old) using an enzymatic methyl-seq approach followed by analysis to call differentially methylated regions (DMRs) within these cohorts. Small RNA (sRNA) expression was profiled in three breeding selections, each from two age cohorts (1 and 6 years old), using sRNA-Seq followed by differential expression analysis. Weighted chromosome-level methylation analysis reveals hypermethylation in 11-year-old almond breeding selections when compared to 2-year-old selections in the CG and CHH contexts. Seventeen consensus DMRs were identified in all age contrasts. sRNA expression differed significantly between the two age cohorts tested, with significantly decreased expression in sRNAs in the 6-year-old selections compared to the 1-year-old. Almond shows a pattern of hypermethylation and decreased sRNA expression with increased age. Identified DMRs and differentially expressed sRNAs could function as putative biomarkers of age following validation in additional age groups.

Nitric oxide controls shoot meristem activity via regulation of DNA methylation

Despite the importance of Nitric Oxide (NO) as signaling molecule in both plant and animal development, the regulatory mechanisms downstream of NO remain largely unclear. Here, we show that NO is involved in Arabidopsis shoot stem cell control via modifying expression and activity of ARGONAUTE 4 (AGO4), a core component of the RNA-directed DNA Methylation (RdDM) pathway. Mutations in components of the RdDM pathway cause meristematic defects, and reduce responses of the stem cell system to NO signaling. Importantly, we find that the stem cell inducing WUSCHEL transcription factor directly interacts with AGO4 in a NO dependent manner, explaining how these two signaling systems may converge to modify DNA methylation patterns. Taken together, our results reveal that NO signaling plays an important role in controlling plant stem cell homeostasis via the regulation of de novo DNA methylation.

DNA methylation profiles in bovine sperm are associated with daughter fertility

The current decline in dairy cattle fertility has resulted in significant financial losses for dairy farmers. In the past, most efforts to improve dairy cattle fertility have been focused on either management or genetics, while epigenetics have received less attention. In this study, 12 bulls were selected from a provided 100 bull list and studied (High daughter fertility = 6, Low daughter fertility = 6) for Enzymatic methylation sequencing in the Illumina HiSeq platform according to the Canadian daughter fertility index (DFI), sires with high and low daughter fertility have average DFI of 92 and 112.6, respectively. And the bull list provided shows a mean DFI of 103.4. 252 CpGs with methylation differences greater than 20% (q < 0.01) were identified, as well as the top 10 promising DMRs with a 15% methylation difference (q < 1.1e-26). Interestingly, the DMCs and DMRs were found to be distributed more on the X chromosome than on the autosome, and they were covered by gene clusters linked to germ cell formation and development. In conclusion, these findings could enhance our ability to make informed decisions when deciding on superior bulls and advance our understanding of paternal epigenetic inheritance.

ACD15, ACD21, and SLN regulate the accumulation and mobility of MBD6 to silence genes and transposable elements

DNA methylation mediates silencing of transposable elements and genes in part via recruitment of the Arabidopsis MBD5/6 complex, which contains the methyl-CpG binding domain (MBD) proteins MBD5 and MBD6, and the J-domain containing protein SILENZIO (SLN). Here, we characterize two additional complex members: α-crystalline domain (ACD) containing proteins ACD15 and ACD21. We show that they are necessary for gene silencing, bridge SLN to the complex, and promote higher-order multimerization of MBD5/6 complexes within heterochromatin. These complexes are also highly dynamic, with the mobility of MBD5/6 complexes regulated by the activity of SLN. Using a dCas9 system, we demonstrate that tethering the ACDs to an ectopic site outside of heterochromatin can drive a massive accumulation of MBD5/6 complexes into large nuclear bodies. These results demonstrate that ACD15 and ACD21 are critical components of the gene-silencing MBD5/6 complex and act to drive the formation of higher-order, dynamic assemblies at CG methylation (meCG) sites.

ACD15, ACD21 and SLN regulate accumulation and mobility of MBD6 to silence genes and transposable elements

DNA methylation mediates silencing of transposable elements and genes in part via recruitment of the Arabidopsis MBD5/6 complex, which contains the methyl-CpG-binding domain (MBD) proteins MBD5 and MBD6, and the J-domain containing protein SILENZIO (SLN). Here we characterize two additional complex members: α-crystalline domain containing proteins ACD15 and ACD21. We show that they are necessary for gene silencing, bridge SLN to the complex, and promote higher order multimerization of MBD5/6 complexes within heterochromatin. These complexes are also highly dynamic, with the mobility of complex components regulated by the activity of SLN. Using a dCas9 system, we demonstrate that tethering the ACDs to an ectopic site outside of heterochromatin can drive massive accumulation of MBD5/6 complexes into large nuclear bodies. These results demonstrate that ACD15 and ACD21 are critical components of gene silencing complexes that act to drive the formation of higher order, dynamic assemblies.

Whole genome methylation sequencing reveals epigenetic landscape and abnormal expression of FABP5 in extramammary Paget's disease

BACKGROUND

Extramammary Paget's disease (EMPD) is a rare cutaneous malignant tumor with a high recurrence rate after surgery. However, the genetic and epigenetic alterations underlying its pathogenesis remain unknown. DNA methylation is an important epigenetic modification involved in many biological processes.

METHODS

In this study, enzymatic methyl-sequencing (EM-seq) technique was used to investigate the landscape of genome-wide DNA methylation from three pairs of tumor tissues and adjacent tissues of patients with EMPD. Additionally, we conducted histopathological examinations to assess the expression of fatty acid-binding protein 5 (FABP5) in another three paired samples from EMPD patients.

RESULTS

The cluster analysis showed the good quality of the samples. A differential methylation region (DMR) heat map was used to quantitatively characterize genome-wide methylation differences between tumors and controls. Global DNA methylation level is lower in EMPD tissue compared to matched controls, indicating that DNA methylation discriminates between tumor and normal skin. And the top hypomethylation gene on the promoter region in tumor tissues was FABP5 on chromosome 8 with 38.44% decreased median methylation. We next identified the expression of FABP5 in paired tumors and adjacent tissues in three additional patients with EMPD. Immunofluorescence results showed FABP5 highly expressed in tumor tissues and co-located with CK7, CK20 and EMA. GO and KEGG enrichment analysis showed DMR genes on promoter are mainly enriched in the calcium ion transport, GTPase mediated signal transduction, Rap1 signaling pathway and GnRH signaling pathway.

CONCLUSION

Taken together, our findings provide the first description of the whole genome methylation map of EMPD and identify FABP5 as a pathogenic target of EMPD.

Navigating the pitfalls of mapping DNA and RNA modifications

Chemical modifications to nucleic acids occur across the kingdoms of life and carry important regulatory information. Reliable high-resolution mapping of these modifications is the foundation of functional and mechanistic studies, and recent methodological advances based on next-generation sequencing and long-read sequencing platforms are critical to achieving this aim. However, mapping technologies may have limitations that sometimes lead to inconsistent results. Some of these limitations are technical in nature and specific to certain types of technology. Here, however, we focus on common (yet not always widely recognized) pitfalls that are shared among frequently used mapping technologies and discuss strategies to help technology developers and users mitigate their effects. Although the emphasis is primarily on DNA modifications, RNA modifications are also discussed.

The plant response to high CO2 levels is heritable and orchestrated by DNA methylation

Plant responses to abiotic environmental challenges are known to have lasting effects on the plant beyond the initial stress exposure. Some of these lasting effects are transgenerational, affecting the next generation. The plant response to elevated carbon dioxide (CO₂ ) levels has been well studied. However, these investigations are typically limited to plants grown for a single generation in a high CO₂ environment while transgenerational studies are rare. We aimed to determine transgenerational growth responses in plants after exposure to high CO₂ by investigating the direct progeny when returned to baseline CO₂ levels. We found that both the flowering plant Arabidopsis thaliana and seedless nonvascular plant Physcomitrium patens continue to display accelerated growth rates in the progeny of plants exposed to high CO₂ . We used the model species Arabidopsis to dissect the molecular mechanism and found that DNA methylation pathways are necessary for heritability of this growth response. More specifically, the pathway of RNA-directed DNA methylation is required to initiate methylation and the proteins CMT2 and CMT3 are needed for the transgenerational propagation of this DNA methylation to the progeny plants. Together, these two DNA methylation pathways establish and then maintain a cellular memory to high CO₂ exposure.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

pycoMeth: a toolbox for differential methylation testing from Nanopore methylation calls

We present pycoMeth, a toolbox to store, manage and analyze DNA methylation calls from long-read sequencing data obtained using the Oxford Nanopore Technologies sequencing platform. Building on a novel, rapid-access, read-level and reference-anchored methylation storage format MetH5, we propose efficient algorithms for haplotype aware, multi-sample consensus segmentation and differential methylation testing. We show that MetH5 is more efficient than existing solutions for storing Oxford Nanopore Technologies methylation calls, and carry out benchmarking for pycoMeth segmentation and differential methylation testing, demonstrating increased performance and sensitivity compared to existing solutions designed for short-read methylation data.

Single-cell DNA methylation sequencing by combinatorial indexing and enzymatic DNA methylation conversion

BACKGROUND

DNA methylation is a critical molecular mark involved in cellular differentiation and cell-specific processes. Single-cell whole genome DNA methylation profiling methods hold great potential to resolve the DNA methylation profiles of individual cell-types. Here we present a method that couples single-cell combinatorial indexing (sci) with enzymatic conversion (sciEM) of unmethylated cytosines.

RESULTS

The sciEM method facilitates DNA methylation profiling of single-cells that is highly correlated with single-cell bisulfite-based workflows (r² > 0.99) whilst improving sequencing alignment rates, reducing adapter contamination and over-estimation of DNA methylation levels (CpG and non-CpG). As proof-of-concept we perform sciEM analysis of the temporal lobe, motor cortex, hippocampus and cerebellum of the human brain to resolve single-cell DNA methylation of all major cell-types.

CONCLUSION

To our knowledge sciEM represents the first non-bisulfite single-cell DNA methylation sequencing approach with single-base resolution.

Target-enriched enzymatic methyl sequencing: Flexible, scalable and inexpensive hybridization capture for quantifying DNA methylation

The increasing interest in studying DNA methylation to understand how traits or diseases develop requires new and flexible approaches for quantifying DNA methylation in a diversity of organisms. In particular, we need efficient yet cost-effective ways to measure CpG methylation states over large and complete regions of the genome. Here, we develop TEEM-Seq (target-enriched enzymatic methyl sequencing), a method that combines enzymatic methyl sequencing with a custom-designed hybridization capture bait set that can be scaled to reactions including large numbers of samples in any species for which a reference genome is available. Using DNA from a passerine bird, the superb starling (Lamprotornis superbus), we show that TEEM-Seq is able to quantify DNA methylation states similarly well to the more traditional approaches of whole-genome and reduced-representation sequencing. Moreover, we demonstrate its reliability and repeatability, as duplicate libraries from the same samples were highly correlated. Importantly, the downstream bioinformatic analysis for TEEM-Seq is the same as for any sequence-based approach to studying DNA methylation, making it simple to incorporate into a variety of workflows. We believe that TEEM-Seq could replace traditional approaches for studying DNA methylation in candidate genes and pathways, and be effectively paired with other whole-genome or reduced-representation sequencing approaches to increase project sample sizes. In addition, TEEM-Seq can be combined with mRNA sequencing to examine how DNA methylation in promoters or other regulatory regions is related to the expression of individual genes or gene networks. By maximizing the number of samples in the hybridization reaction, TEEM-Seq is an inexpensive and flexible sequence-based approach for quantifying DNA methylation in species where other capture-based methods are unavailable or too expensive, particularly for non-model organisms.

Performance comparisons of methylation and structural variants from low-input whole-genome methylation sequencing

Aim: Whole-genome methylation sequencing carries both DNA methylation and structural variant information (single nucleotide variant [SNV]; copy number variant [CNV]); however, limited data is available on the reliability of obtaining this information simultaneously from low-input DNA using various library preparation and sequencing protocols. Methods: A HapMap NA12878 sample was sequenced with three protocols (EM-sequencing, QIA-sequencing and Swift-sequencing) and their performance was compared on CpG methylation measurement and SNV and CNV detection. Results: At low DNA input (10-25 ng), EM-sequencing was superior in almost all metrics except CNV detection where all protocols were similar. EM-sequencing captured the highest number of CpGs and true SNVs. Conclusion: EM-sequencing is suitable to detect methylation, SNVs and CNVs from single sequencing with low-input DNA.

Identification of unique DNA methylation sites in Kabuki syndrome using whole genome bisulfite sequencing and targeted hybridization capture followed by enzymatic methylation sequencing

BACKGROUND

Kabuki syndrome (KS) is a congenital malformation syndrome caused by mutations in the KMT2D and KDM6A genes that encode histone modification enzymes. Although KS is considered a single gene disorder, its symptoms vary widely. Recently, disease-specific DNA methylation patterns, or episignatures, have been recognized and used as a diagnostic tool for KS. Because of various crosstalk mechanisms between histone modifications and DNA methylation, DNA methylation analysis may have high potential for investigations into the pathogenesis of KS.

RESULTS

In this study, we investigated altered CpG-methylation sites that were specific to KS to find important genes associated with the various phenotypes or pathogenesis of KS. Whole genome bisulfite sequencing (WGBS) was performed to select target CpG islands, and enzymatic conversion technology was applied after hybridization capture to confirm KS-specific episignatures of 130 selected differently methylated target regions (DMTRs) in DNA samples from the 65 participants, 31 patients with KS and 34 unaffected individuals, in this study. We identified 26 candidate genes in 22 DMTRs that may be associated with KS. Our results indicate that disease-specific methylation sites can be identified from a small number of WGBS samples, and hybridization capture followed by enzymatic methylation sequencing can simultaneously test the sites.

CONCLUSIONS

Although DNA methylation can be tissue-specific, our results suggest that methylation profiling of DNA extracted from peripheral blood may be a powerful approach to study the pathogenesis of diseases.

Conditional GWAS of non-CG transposon methylation in Arabidopsis thaliana reveals major polymorphisms in five genes

Genome-wide association studies (GWAS) have revealed that the striking natural variation for DNA CHH-methylation (mCHH; H is A, T, or C) of transposons has oligogenic architecture involving major alleles at a handful of known methylation regulators. Here we use a conditional GWAS approach to show that CHG-methylation (mCHG) has a similar genetic architecture-once mCHH is statistically controlled for. We identify five key trans-regulators that appear to modulate mCHG levels, and show that they interact with a previously identified modifier of mCHH in regulating natural transposon mobilization.

Epigenetic regulation of fetal brain development in pig

How fetal brain development is regulated at the molecular level is not well understood. Due to ethical challenges associated with research on the human fetus, large animals particularly pigs are increasingly used to study development and disorders of fetal brain. The pig fetal brain grows rapidly during the last ∼ 50 days before birth which is around day 60 (d60) of pig gestation. But what regulates the onset of accelerated growth of the brain is unknown. The current study tests the hypothesis that epigenetic alteration around d60 is involved in the onset of rapid growth of fetal brain of pig. To test this hypothesis, DNA methylation changes of fetal brain was assessed in a genome-wide manner by Enzymatic Methyl-seq (EM-seq) during two gestational periods (GP): d45 vs. d60 (GP1) and d60 vs. d90 (GP2). The cytosine-guanine (CpG) methylation data was analyzed in an integrative manner with the RNA-seq data generated from the same brain samples from our earlier study. A neural network based modeling approach was implemented to learn changes in methylation patterns of the differentially expressed genes, and then predict methylations of the brain in a genome-wide manner during rapid growth. This approach identified specific methylations that changed in a mutually informative manner during rapid growth of the fetal brain. These methylations were significantly overrepresented in specific genic as well as intergenic features including CpG islands, introns, and untranslated regions. In addition, sex-bias methylations of known single nucleotide polymorphic sites were also identified in the fetal brain ide during rapid growth.

Mutations in DNA polymerase δ subunit 1 co-segregate with CMD2-type resistance to Cassava Mosaic Geminiviruses

Cassava mosaic disease (CMD) suppresses cassava yields across the tropics. The dominant CMD2 locus confers resistance to cassava mosaic geminiviruses. It has been reported that CMD2-type landraces lose resistance after regeneration through de novo morphogenesis. As full genome bisulfite sequencing failed to uncover an epigenetic mechanism for this loss of resistance, whole genome sequencing and genetic variant analysis was performed and the CMD2 locus was fine-mapped to a 190 kilobase interval. Collectively, these data indicate that CMD2-type resistance is caused by a nonsynonymous, single nucleotide polymorphism in DNA polymerase δ subunit 1 (MePOLD1) located within this region. Virus-induced gene silencing of MePOLD1 in a CMD-susceptible cassava variety produced a recovery phenotype typical of CMD2-type resistance. Analysis of other CMD2-type cassava varieties identified additional candidate resistance alleles within MePOLD1. Genetic variation of MePOLD1, therefore, could represent an important genetic resource for resistance breeding and/or genome editing, and elucidating mechanisms of resistance to geminiviruses.

Whole-genome sequence and methylome profiling of the almond [Prunus dulcis (Mill.) D.A. Webb] cultivar 'Nonpareil'

Almond [Prunus dulcis (Mill.) D.A. Webb] is an economically important, specialty nut crop grown almost exclusively in the United States. Breeding and improvement efforts worldwide have led to the development of key, productive cultivars, including 'Nonpareil,' which is the most widely grown almond cultivar. Thus far, genomic resources for this species have been limited, and a whole-genome assembly for 'Nonpareil' is not currently available despite its economic importance and use in almond breeding worldwide. We generated a 571X coverage genome sequence using Illumina, PacBio, and optical mapping technologies. Gene prediction revealed 49,321 putative genes using MinION Oxford nanopore and Illumina RNA sequencing, and genome annotation found that 68% of predicted models are associated with at least one biological function. Furthermore, epigenetic signatures of almond, namely DNA cytosine methylation, have been implicated in a variety of phenotypes including self-compatibility, bud dormancy, and development of noninfectious bud failure. In addition to the genome sequence and annotation, this report also provides the complete methylome of several almond tissues, including leaf, flower, endocarp, mesocarp, exocarp, and seed coat. Comparisons between methylation profiles in these tissues revealed differences in genome-wide weighted % methylation and chromosome-level methylation enrichment.

Antioxidant Regulation and DNA Methylation Dynamics During Mikania micrantha Seed Germination Under Cold Stress

As a primary goal, adaptation to cold climate could expand an invasion range of exotic plants. Here, we aimed to explore the regulation strategy of M. micrantha seed development under cold stress through molecular physiology and multi-omics analysis. Significant increase of hydrogen peroxide, malondialdehyde, and electrolyte leakage observed under cold stress revealed that oxidative damage within M. micrantha seed cells was induced in the initial germination phase. Proteomic data underscored an activation of antioxidant activity to maintain redox homeostasis, with a cluster of antioxidant proteins identified. Genomic-wide transcriptome, in combination with time-series whole-genome bisulfite sequencing mining, elucidated that seven candidate genes, which were the target of DNA demethylation-dependent ROS scavenging, were possibly associated with an M. micrantha germ break. Progressive gain of CHH context DNA methylation identified in an early germination phrase suggested a role of a DNA methylation pathway, while an active DNA demethylation pathway was also initiated during late seed development, which was in line with the expression trend of methylation and demethylation-related genes verified through qRT-PCR. These data pointed out that cold-dependent DNA demethylation and an antioxidant regulatory were involved together in restoring seed germination. The expression level of total 441 genes presented an opposite trend to the methylation divergence, while the expression of total 395 genes was proved to be negatively associated with their methylation levels. These data provided new insights into molecular reprograming events during M. micrantha seed development.

Somatic DNA demethylation generates tissue-specific methylation states and impacts flowering time

Cytosine methylation is a reversible epigenetic modification of DNA. In plants, removal of cytosine methylation is accomplished by the four members of the DEMETER (DME) family of 5-methylcytosine DNA glycosylases, named DME, DEMETER-LIKE2 (DML2), DML3, and REPRESSOR OF SILENCING1 (ROS1) in Arabidopsis thaliana. Demethylation by DME is critical for seed development, preventing experiments to determine the function of the entire gene family in somatic tissues by mutant analysis. Here, we bypassed the reproductive defects of dme mutants to create somatic quadruple homozygous mutants of the entire DME family. dme; ros1; dml2; and dml3 (drdd) leaves exhibit hypermethylated regions compared with wild-type leaves and rdd triple mutants, indicating functional redundancy among all four demethylases. Targets of demethylation include regions co-targeted by RNA-directed DNA methylation and, surprisingly, CG gene body methylation, indicating dynamic methylation at these less-understood sites. Additionally, many tissue-specific methylation differences are absent in drdd, suggesting a role for active demethylation in generating divergent epigenetic states across wild-type tissues. Furthermore, drdd plants display an early flowering phenotype, which involves 5'-hypermethylation and transcriptional down-regulation of FLOWERING LOCUS C. Active DNA demethylation is therefore required for proper methylation across somatic tissues and defines the epigenetic landscape of intergenic and coding regions.

Analysis and Performance Assessment of the Whole Genome Bisulfite Sequencing Data Workflow: Currently Available Tools and a Practical Guide to Advance DNA Methylation Studies

DNA methylation is associated with transcriptional repression, genomic imprinting, stem cell differentiation, embryonic development, and inflammation. Aberrant DNA methylation can indicate disease states, including cancer and neurological disorders. Therefore, the prevalence and location of 5-methylcytosine in the human genome is a topic of interest. Whole-genome bisulfite sequencing (WGBS) is a high-throughput method for analyzing DNA methylation. This technique involves library preparation, alignment, and quality control. Advancements in epigenetic technology have led to an increase in DNA methylation studies. This review compares the detailed experimental methodology of WGBS using accessible and up-to-date analysis tools. Practical codes for WGBS data processing are included as a general guide to assist progress in DNA methylation studies through a comprehensive case study.

Whole Genome DNA Methylation Profiling of D2 Medium Spiny Neurons in Mouse Nucleus Accumbens Using Two Independent Library Preparation Methods

DNA methylation plays essential roles in various cellular processes. Next-generation sequencing has enabled us to study the functional implication of DNA methylation across the whole genome. However, this approach usually requires a substantial amount of genomic DNA, which limits its application to defined cell types within a discrete brain region. Here, we applied two separate protocols, Accel-NGS Methyl-Seq (AM-seq) and Enzymatic Methyl-seq (EM-seq), to profile the methylome of D2 dopamine receptor-expressing medium spiny neurons (D2-MSNs) in mouse nucleus accumbens (NAc). Using 40 ng DNA extracted from FACS-isolated D2-MSNs, we found that both methods yielded comparably high-quality methylome data. Additionally, we identified numerous unmethylated regions (UMRs) as cell type-specific regulatory regions. By comparing the NAc D2-MSN methylome with the published methylomes of mouse prefrontal cortex excitatory neurons and neural progenitor cells (NPCs), we identified numerous differentially methylated CpG and non-CpG regions. Our study not only presents a comparison of these two low-input DNA whole genome methylation profiling protocols, but also provides a resource of DNA methylome of mouse accumbal D2-MSNs, a neuron type that has critical roles in addiction and other neuropsychiatric disorders.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Epigenome guided crop improvement: current progress and future opportunities

Epigenomics encompasses a broad field of study, including the investigation of chromatin states, chromatin modifications and their impact on gene regulation; as well as the phenomena of epigenetic inheritance. The epigenome is a multi-modal layer of information superimposed on DNA sequences, instructing their usage in gene expression. As such, it is an emerging focus of efforts to improve crop performance. Broadly, this might be divided into avenues that leverage chromatin information to better annotate and decode plant genomes, and into complementary strategies that aim to identify and select for heritable epialleles that control crop traits independent of underlying genotype. In this review, we focus on the first approach, which we term ‘epigenome guided’ improvement. This encompasses the use of chromatin profiles to enhance our understanding of the composition and structure of complex crop genomes. We discuss the current progress and future prospects towards integrating this epigenomic information into crop improvement strategies; in particular for CRISPR/Cas9 gene editing and precision genome engineering. We also highlight some specific opportunities and challenges for grain and horticultural crops.

Quality control and evaluation of plant epigenomics data

Epigenomics is the study of molecular signatures associated with discrete regions within genomes, many of which are important for a wide range of nuclear processes. The ability to profile the epigenomic landscape associated with genes, repetitive regions, transposons, transcription, differential expression, cis-regulatory elements, and 3D chromatin interactions has vastly improved our understanding of plant genomes. However, many epigenomic and single-cell genomic assays are challenging to perform in plants, leading to a wide range of data quality issues; thus, the data require rigorous evaluation prior to downstream analyses and interpretation. In this commentary, we provide considerations for the evaluation of plant epigenomics and single-cell genomics data quality with the aim of improving the quality and utility of studies using those data across diverse plant species.

Comparison of EM-seq and PBAT methylome library methods for low-input DNA

DNA methylation is the most studied epigenetic mark involved in regulation of gene expression. For low input samples, a limited number of methods for quantifying DNA methylation genome-wide has been evaluated. Here, we compared a series of input DNA amounts (1–10ng) from two methylome library preparation protocols, enzymatic methyl-seq (EM-seq) and post-bisulfite adaptor tagging (PBAT) adapted from single-cell PBAT. EM-seq takes advantage of enzymatic activity while PBAT relies on conventional bisulfite conversion for detection of DNA methylation. We found that both methods accurately quantified DNA methylation genome-wide. They produced expected distribution patterns around genomic features, high C-T transition efficiency at non-CpG sites and high correlation between input amounts. However, EM-seq performed better in regard to library and sequencing quality, i.e. EM-seq produced larger insert sizes, higher alignment rates and higher library complexity with lower duplication rate compared to PBAT. Moreover, EM-seq demonstrated higher CpG coverage, better CpG site overlap and higher consistency between input series. In summary, our data suggests that EM-seq overall performed better than PBAT in whole-genome methylation quantification of low input samples.

Beyond the gene: epigenetic and cis-regulatory targets offer new breeding potential for the future

For millennia, natural and artificial selection has combined favourable alleles for desirable traits in crop species. While modern plant breeding has achieved steady increases in crop yields over the last century, on the current trajectory we will simply not meet demand by 2045. Novel breeding strategies and sources of genetic variation will be required to sustainably fill predicted yield gaps and meet new consumer preferences. Here, we highlight that stepping up to meet this grand challenge will increasingly require thinking 'beyond the gene'. Significant progress has been made in understanding the contributions of both epigenetic variation and cis-regulatory variation to plant traits. This non-genic variation has great potential in future breeding, synthetic biology and biotechnology applications.

Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA

Bisulfite sequencing detects 5mC and 5hmC at single-base resolution. However, bisulfite treatment damages DNA, which results in fragmentation, DNA loss, and biased sequencing data. To overcome these problems, enzymatic methyl-seq (EM-seq) was developed. This method detects 5mC and 5hmC using two sets of enzymatic reactions. In the first reaction, TET2 and T4-BGT convert 5mC and 5hmC into products that cannot be deaminated by APOBEC3A. In the second reaction, APOBEC3A deaminates unmodified cytosines by converting them to uracils. Therefore, these three enzymes enable the identification of 5mC and 5hmC. EM-seq libraries were compared with bisulfite-converted DNA, and each library type was ligated to Illumina adaptors before conversion. Libraries were made using NA12878 genomic DNA, cell-free DNA, and FFPE DNA over a range of DNA inputs. The 5mC and 5hmC detected in EM-seq libraries were similar to those of bisulfite libraries. However, libraries made using EM-seq outperformed bisulfite-converted libraries in all specific measures examined (coverage, duplication, sensitivity, etc.). EM-seq libraries displayed even GC distribution, better correlations across DNA inputs, increased numbers of CpGs within genomic features, and accuracy of cytosine methylation calls. EM-seq was effective using as little as 100 pg of DNA, and these libraries maintained the described advantages over bisulfite sequencing. EM-seq library construction, using challenging samples and lower DNA inputs, opens new avenues for research and clinical applications.

BiasCorrector: Fast and accurate correction of all types of experimental biases in quantitative DNA methylation data derived by different technologies

Quantification of DNA methylation in neoplastic cells is crucial both from mechanistic and diagnostic perspectives. However, such measurements are prone to different experimental biases. Polymerase chain reaction (PCR) bias results in an unequal recovery of methylated and unmethylated alleles at the sample preparation step. Post-PCR biases get introduced additionally by the readout processes. Correcting the biases is more practicable than optimising experimental conditions, as demonstrated previously. However, utilisation of our earlier developed algorithm strongly necessitates automation. Here, we present two R packages: rBiasCorrection, the core algorithms to correct biases; and BiasCorrector, its web-based graphical user interface frontend. The software detects and analyses experimental biases in calibration DNA samples at a single base resolution by using cubic polynomial and hyperbolic regression. The correction coefficients from the best regression type are employed to compensate for the bias. Three common technologies-bisulphite pyrosequencing, next-generation sequencing and oligonucleotide microarrays-were used to comprehensively test BiasCorrector. We demonstrate the accuracy of BiasCorrector's performance and reveal technology-specific PCR- and post-PCR biases. BiasCorrector effectively eliminates biases regardless of their nature, locus, the number of interrogated methylation sites and the detection method, thus representing a user-friendly tool for producing accurate epigenetic results.

Developmental remodelling of non-CG methylation at satellite DNA repeats

In vertebrates, DNA methylation predominantly occurs at CG dinucleotides however, widespread non-CG methylation (mCH) has been reported in mammalian embryonic stem cells and in the brain. In mammals, mCH is found at CAC trinucleotides in the nervous system, where it is associated with transcriptional repression, and at CAG trinucleotides in embryonic stem cells, where it positively correlates with transcription. Moreover, CAC methylation appears to be a conserved feature of adult vertebrate brains. Unlike any of those methylation signatures, here we describe a novel form of mCH that occurs in the TGCT context within zebrafish mosaic satellite repeats. TGCT methylation is inherited from both male and female gametes, remodelled during mid-blastula transition, and re-established during gastrulation in all embryonic layers. Moreover, we identify DNA methyltransferase 3ba (Dnmt3ba) as the primary enzyme responsible for the deposition of this mCH mark. Finally, we observe that TGCT-methylated repeats are specifically associated with H3K9me3-marked heterochromatin suggestive of a functional interplay between these two gene-regulatory marks. Altogether, this work provides insight into a novel form of vertebrate mCH and highlights the substrate diversity of vertebrate DNA methyltransferases.