DNA Methylation: Bisulphite Modification and Analysis

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene regulation is complex, but essentially involves modifications to the DNA itself or the proteins with which DNA associates. The predominant epigenetic modification of DNA in mammalian genomes is methylation of cytosine nucleotides (5-MeC). DNA methylation provides instruction to gene expression machinery as to where and when the gene should be expressed. The primary target sequence for DNA methylation in mammals is 5''-CpG-3'' dinucleotides (Figure 1). CpG dinucleotides are not uniformly distributed throughout the genome, but are concentrated in regions of repetitive genomic sequences and CpG "islands" commonly associated with gene promoters (Figure 1). DNA methylation patterns are established early in development, modulated during tissue specific differentiation and disrupted in many disease states including cancer. To understand the biological role of DNA methylation and its role in human disease, precise, efficient and reproducible methods are required to detect and quantify individual 5-MeCs.This protocol for bisulphite conversion is the "gold standard" for DNA methylation analysis and facilitates identification and quantification of DNA methylation at single nucleotide resolution. The chemistry of cytosine deamination by sodium bisulphite involves three steps (Figure 2). (1) Sulphonation: The addition of bisulphite to the 5-6 double bond of cytosine (2) Hydrolic Deamination: hydrolytic deamination of the resulting cytosine-bisulphite derivative to give a uracil-bisulphite derivative (3) Alkali Desulphonation: Removal of the sulphonate group by an alkali treatment, to give uracil. Bisulphite preferentially deaminates cytosine to uracil in single stranded DNA, whereas 5-MeC, is refractory to bisulphite-mediated deamination. Upon PCR amplification, uracil is amplified as thymine while 5-MeC residues remain as cytosines, allowing methylated CpGs to be distinguished from unmethylated CpGs by presence of a cytosine "C" versus thymine "T" residue during sequencing.DNA modification by bisulphite conversion is a well-established protocol that can be exploited for many methods of DNA methylation analysis. Since the detection of 5-MeC by bisulphite conversion was first demonstrated by Frommer et al.¹ and Clark et al.², methods based around bisulphite conversion of genomic DNA account for the majority of new data on DNA methylation. Different methods of post PCR analysis may be utilized, depending on the degree of specificity and resolution of methylation required. Cloning and sequencing is still the most readily available method that can give single nucleotide resolution for methylation across the DNA molecule.     

The CHH motif in sugar beet satellite DNA: a modulator for cytosine methylation

Methylation of DNA is important for the epigenetic silencing of repetitive DNA in plant genomes. Knowledge about the cytosine methylation status of satellite DNAs, a major class of repetitive DNA, is scarce. One reason for this is that arrays of tandemly arranged sequences are usually collapsed in next‐generation sequencing assemblies. We applied strategies to overcome this limitation and quantified the level of cytosine methylation and its pattern in three satellite families of sugar beet (Beta vulgaris) which differ in their abundance, chromosomal localization and monomer size. We visualized methylation levels along pachytene chromosomes with respect to small satellite loci at maximum resolution using chromosome‐wide fluorescent in situ hybridization complemented with immunostaining and super‐resolution microscopy. Only reduced methylation of many satellite arrays was obtained. To investigate methylation at the nucleotide level we performed bisulfite sequencing of 1569 satellite sequences. We found that the level of methylation of cytosine strongly depends on the sequence context: cytosines in the CHH motif show lower methylation (44–52%), while CG and CHG motifs are more strongly methylated. This affects the overall methylation of satellite sequences because CHH occurs frequently while CG and CHG are rare or even absent in the satellite arrays investigated. Evidently, CHH is the major target for modulation of the cytosine methylation level of adjacent monomers within individual arrays and contributes to their epigenetic function. This strongly indicates that asymmetric cytosine methylation plays a role in the epigenetic modification of satellite repeats in plant genomes.     

Hi-Meth：特定位点DNA甲基化高通量检测平台

胞嘧啶甲基化是DNA表观遗传修饰的主要类型之一,在维持正常细胞功能和调控基因表达中具有重要作用。重亚硫酸盐测序法(bisulfite sequencing PCR,BSP)是特异性位点DNA甲基化检测的通用方法,能明确目的片段中每一个CpG位点的甲基化状态,但此方法需要大量的单克隆测序,操作过程较繁琐、成本昂贵。因此,开发准确、高效、便捷的DNA甲基化检测技术对提升表观遗传研究效率具有重要意义。基于本课题组开发的高通量突变类型检测平台Hi-TOM (high-throughput tracking of mutations),我们进一步建立了特定位点DNA甲基化高通量检测平台Hi-Meth (high-throughput detection of DNA methylation)。DNA样品通过重亚硫酸盐处理之后,仅需一轮PCR扩增即可通过Hi-Meth平台获得特定位点DNA甲基化分析结果。利用Hi-Meth平台,对水稻不同基因启动子区域进行了DNA甲基化检测分析,并与基于BSP方法获得的结果进行了比较。结果表明,Hi-Meth策略与BSP策略检测结果基本一致。而且通过Hi-Meth平台可以更准确、便捷地获得特异性位点DNA甲基化分析结果。综上所述,Hi-Meth为特定DNA区域提供了重要的甲基化检测平台,对表观遗传研究具有重要意义。     

Cooperative activity of DNA methyltransferases for maintenance of symmetrical and non-symmetrical cytosine methylation in Arabidopsis thaliana

Maintenance of cytosine methylation in plants is controlled by three DNA methyltransferases. MET1 maintains CG methylation, and DRM1/2 and CMT3 act redundantly to enforce non-CG methylation. RPS, a repetitive hypermethylated DNA fragment from Petunia hybrida, attracts DNA methylation when transferred into Petunia or other species. In Arabidopsis thaliana, which does not contain any RPS homologues, RPS transgenes are efficiently methylated in all sequence contexts. To test which DNA methylation pathways regulate RPS methylation, we examined maintenance of RPS methylation in various mutant backgrounds. Surprisingly, CG methylation was lost in a drm1/2/cmt3 mutant, and non-CG methylation was almost completely eliminated in a met1 mutant. An unusual cooperative activity of all three DNA methyltransferases is therefore required for maintenance of both CG and non-CG methylation in RPS. Other unusual features of RPS methylation are the independence of its non-CG methylation from the RNA-directed DNA methylation (RdDM) pathway and the exceptional maintenance of methylation at a CC(m)TGG site in some epigenetic mutants. This is indicative of activity of a methylation system in plants that may have evolved from the DCM methylation system that controls CC(m)WGG methylation in bacteria. Our data suggest that strict separation of CG and non-CG methylation pathways does not apply to all target regions, and that caution is required in generalizing methylation data obtained for individual genomic regions.     

Detection and measurement of PCR bias in quantitative methylation analysis of bisulphite-treated DNA.

Methylation analysis of individual cytosines in genomic DNA can be determined quantitatively by bisulphite treatment and PCR amplification of the target DNA sequence, followed by restriction enzyme digestion or sequencing. Methylated and unmethylated molecules, however, have different sequences after bisulphite conversion. For some sequences this can result in bias during the PCR amplification leading to an inaccurate estimate of methylation. PCR bias is sequence dependent and often strand-specific. This study presents a simple method for detection and measurement of PCR bias for any set of primers, and investigates parameters for overcoming PCR bias.     

Genome-wide high throughput analysis of DNA methylation in eukaryotes

Cytosine methylation is the quintessential epigenetic mark. Two well-established methods, bisulfite sequencing and methyl-DNA immunoprecipitation (MeDIP) lend themselves to the genome-wide analysis of DNA methylation by high throughput sequencing. Here we provide an overview and brief review of these methods. We summarize our experience with MeDIP followed by high throughput Illumina/Solexa sequencing, exemplified by the analysis of the methylated fraction of the Neurospora crassa genome ("methylome"). We provide detailed methods for DNA isolation, processing and the generation of in vitro libraries for Illumina/Solexa sequencing. We discuss potential problems in the generation of sequencing libraries. Finally, we provide an overview of software that is appropriate for the analysis of high throughput sequencing data generated by Illumina/Solexa-type sequencing by synthesis, with a special emphasis on approaches and applications that can generate more accurate depictions of sequence reads that fall in repeated regions of a chosen reference genome.     

DNA Methylation Analysis of Chromosome 21 Gene Promoters at Single Base Pair and Single Allele Resolution

Differential DNA methylation is an essential epigenetic signal for gene regulation, development, and disease processes. We mapped DNA methylation patterns of 190 gene promoter regions on chromosome 21 using bisulfite conversion and subclone sequencing in five human cell types. A total of 28,626 subclones were sequenced at high accuracy using (long-read) Sanger sequencing resulting in the measurement of the DNA methylation state of 580427 CpG sites. Our results show that average DNA methylation levels are distributed bimodally with enrichment of highly methylated and unmethylated sequences, both for amplicons and individual subclones, which represent single alleles from individual cells. Within CpG-rich sequences, DNA methylation was found to be anti-correlated with CpG dinucleotide density and GC content, and methylated CpGs are more likely to be flanked by AT-rich sequences. We observed over-representation of CpG sites in distances of 9, 18, and 27 bps in highly methylated amplicons. However, DNA sequence alone is not sufficient to predict an amplicon's DNA methylation status, since 43% of all amplicons are differentially methylated between the cell types studied here. DNA methylation in promoter regions is strongly correlated with the absence of gene expression and low levels of activating epigenetic marks like H3K4 methylation and H3K9 and K14 acetylation. Utilizing the single base pair and single allele resolution of our data, we found that i) amplicons from different parts of a CpG island frequently differ in their DNA methylation level, ii) methylation levels of individual cells in one tissue are very similar, and iii) methylation patterns follow a relaxed site-specific distribution. Furthermore, iv) we identified three cases of allele-specific DNA methylation on chromosome 21. Our data shed new light on the nature of methylation patterns in human cells, the sequence dependence of DNA methylation, and its function as epigenetic signal in gene regulation. Further, we illustrate genotype–epigenotype interactions by showing novel examples of allele-specific methylation.     

铜胁迫对不同抗性种群海州香薷酸性转化酶基因启动子甲基化的影响

采用亚硫酸氢盐测序技术,检测了Cu胁迫对海州香薷（Elsholtzia haichowensis Sun） Cu抗性和非抗性种群酸性转化酶基因启动子甲基化的影响。结果表明,海州香薷液泡转化酶基因（vINV）和细胞壁转化酶基因（cwINV）的启动子在CG位点分别表现出超甲基化和低甲基化的现象。酸性转化酶基因启动子CHG和CHH位点的甲基化状态受Cu胁迫影响较大。Cu胁迫下,vINV和cwINV启动子分别有1个CHG和6个CHH位点的甲基化状态在Cu抗性种群和非抗性种群之间表现出较大差异。抗Cu种群中这些甲基化差异位点对Cu胁迫不敏感,但在非抗性种群中这些位点的甲基化水平在Cu胁迫后出现大幅上升或下降。有些甲基化差异位点位于或者临近预测的启动子顺式作用元件区域,可能参与Cu胁迫下酸性转化酶基因的表达调控。     

DNA methylation in mouse A-repeats in DNA methyltransferase-knockout ES cells and in normal cells determined by bisulfite genomic sequencing

Mouse ES cells with a null mutation of the known DNA methyltransferase retain some residual DNA methylation and can methylate foreign sequences de novo. We have used bisulfite genomic sequencing to examine the sequence specificity and distributions of methylation of a hypermethylated CG island sequence, mouse A-repeats. There were 13 CG dinucleotides in the region examined, 12 of which were methylated to variable extents in all DNAs. We found that: (1) there is considerable residual DNA methylation in ES cells lacking the known DNA methyltransferase (29% of normal methylation in the complete knockout ES DNA); (2) this other activity methylates at exactly the same CG sites as the major methyltransferase; and (3) differences in the distribution of methylated sites between A-repeats in these DNAs are consistent with this other activity methylating in a random de novo fashion. Also, the lack of any methylation in non-CG sites argues that, in other studies where non-CG methylation sites have been found by bisulfite sequencing, detection of such sites of non-CG methylation is not an inherent artifact in this methodology.     

DNA methylation of retrotransposons,DNA transposons and genes in sugar beet (Beta vulgaris L.)

The methylation of cytosines shapes the epigenetic landscape of plant genomes, coordinates transgenerational epigenetic inheritance, represses the activity of transposable elements (TEs), affects gene expression and, hence, can influence the phenotype. Sugar beet (Beta vulgaris ssp. vulgaris), an important crop that accounts for 30% of worldwide sugar needs, has a relatively small genome size (758 Mbp) consisting of approximately 485 Mbp repetitive DNA (64%), in particular satellite DNA, retrotransposons and DNA transposons. Genome‐wide cytosine methylation in the sugar beet genome was studied in leaves and leaf‐derived callus with a focus on repetitive sequences, including retrotransposons and DNA transposons, the major groups of repetitive DNA sequences, and compared with gene methylation. Genes showed a specific methylation pattern for CG, CHG (H = A, C, and T) and CHH sites, whereas the TE pattern differed, depending on the TE class (class 1, retrotransposons and class 2, DNA transposons). Along genes and TEs, CG and CHG methylation was higher than that of adjacent genomic regions. In contrast to the relatively low CHH methylation in retrotransposons and genes, the level of CHH methylation in DNA transposons was strongly increased, pointing to a functional role of asymmetric methylation in DNA transposon silencing. Comparison of genome‐wide DNA methylation between sugar beet leaves and callus revealed a differential methylation upon tissue culture. Potential epialleles were hypomethylated (lower methylation) at CG and CHG sites in retrotransposons and genes and hypermethylated (higher methylation) at CHH sites in DNA transposons of callus when compared with leaves.     

A comparison of existing global DNA methylation assays to low-coverage whole-genome bisulfite sequencing for epidemiological studies

DNA methylation is an epigenetic mark at the interface of genetic and environmental factors relevant to human disease. Quantitative assessments of global DNA methylation levels have therefore become important tools in epidemiology research, particularly for understanding effects of environmental exposures in complex diseases. Among the available methods of quantitative DNA methylation measurements, bisulfite sequencing is considered the gold standard, but whole-genome bisulfite sequencing (WGBS) has previously been considered too costly for epidemiology studies with high sample numbers. Pyrosequencing of repetitive sequences within bisulfite-treated DNA has been routinely used as a surrogate for global DNA methylation, but a comparison of pyrosequencing to WGBS for accuracy and reproducibility of methylation levels has not been performed. This study compared the global methylation levels measured from uniquely mappable (non-repetitive) WGBS sequences to pyrosequencing assays of several repeat sequences and repeat assay-matched WGBS data and determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels. We determined sources of variation in repetitive pyrosequencing assays to be PCR amplification bias, PCR primer selection bias in methylation levels of targeted sequences, and inherent variability in methylation levels of repeat sequences. Low-coverage, uniquely mappable WGBS showed the strongest correlation between replicates of all assays. By using multiplexing by indexed bar codes, the cost of WGBS can be lowered significantly to improve the accuracy of global DNA methylation assessments for human studies.     

An Integrated Workflow for DNA Methylation Analysis

The analysis of cytosine methylation provides a new way to assess and describe epigenetic regulation at a whole-genome level in many eukaryotes. DNA methylation has a demonstrated role in the genome stability and protection, regulation of gene expression and many other aspects of genome function and maintenance. BS-seq is a relatively unbiased method for profiling the DNA methylation, with a resolution capable of measuring methylation at individual cytosines. Here we describe, as an example, a workflow to handle DNA methylation analysis, from BS-seq library preparation to the data visualization. We describe some applications for the analysis and interpretation of these data. Our laboratory provides public access to plant DNA methylation data via visualization tools available at our “Next-Gen Sequence” websites (http://mpss.udel.edu), along with small RNA, RNA-seq and other data types.