MBD-isolated Genome Sequencing provides a high-throughput and comprehensive survey of DNA methylation in the human genome

DNA methylation is an epigenetic modification involved in both normal developmental processes and disease states through the modulation of gene expression and the maintenance of genomic organization. Conventional methods of DNA methylation analysis, such as bisulfite sequencing, methylation sensitive restriction enzyme digestion and array-based detection techniques, have major limitations that impede high-throughput genome-wide analysis. We describe a novel technique, MBD-isolated Genome Sequencing (MiGS), which combines precipitation of methylated DNA by recombinant methyl-CpG binding domain of MBD2 protein and sequencing of the isolated DNA by a massively parallel sequencer. We utilized MiGS to study three isogenic cancer cell lines with varying degrees of DNA methylation. We successfully detected previously known methylated regions in these cells and identified hundreds of novel methylated regions. This technique is highly specific and sensitive and can be applied to any biological settings to identify differentially methylated regions at the genomic scale.     

MethylScreen: DNA methylation density monitoring using quantitative PCR

Aberrant gene silencing of genes through cytosine methylation has been demonstrated during the development of many types of cancers including prostate cancer Several genes including GSTP1 have been shown to be methylated in prostate cancer leading to the suggestion and demonstration that methylation status of such genes could be used as cancer diagnosis markers alone or in support of histology. We developed a bisulfite-free alternative, MethylScreen technology, an assay for DNA methylation detection utilizing combined restriction from both methylation-sensitive restriction enzymes (MSRE) and methylation-dependent restriction enzymes (MDRE). MethylScreen was used to analyze the 5' region of GSTP1 in cell lines, in vitro methylated DNA populations, and flash-frozen tissue samples in an effort to characterize the output and analytical performance characteristics of the assay. The output from the quantitative PCR assay suggested that it could not only detect fully methylated molecules in a mixed population below the 1% level, but it could also quantify the abundance of intermediately methylated molecules. Interestingly, the interpreted output from the four quantitative PCRs closely resembled the molecular population as described by clone-based bisulfite genomic sequencing.     

Quantifying DNA-protein interactions by double-stranded DNA arrays.

We have created double-stranded oligonucleotide arrays to perform highly parallel investigations of DNA-protein interactions. Arrays of single-stranded DNA oligonucleotides, synthesized by a combination of photolithography and solid-state chemistry, have been used for a variety of applications, including large-scale mRNA expression monitoring, genotyping, and sequence-variation analysis. We converted a single-stranded to a double-stranded array by synthesizing a constant sequence at every position on an array and then annealing and enzymatically extending a complementary primer. The efficiency of second-strand synthesis was demonstrated by incorporation of fluorescently labeled dNTPs (2'-deoxyribonucleoside 5'-triphosphates) and by terminal transferase addition of a fluorescently labeled ddNTP. The accuracy of second-strand synthesis was demonstrated by digestion of the arrayed double-stranded DNA (dsDNA) on the array with sequence-specific restriction enzymes. We showed dam methylation of dsDNA arrays by digestion with DpnI, which cleaves when its recognition site is methylated. This digestion demonstrated that the dsDNA arrays can be further biochemically modified and that the DNA is accessible for interaction with DNA-binding proteins. This dsDNA array approach could be extended to explore the spectrum of sequence-specific protein binding sites in genomes.     

DHPLC-based method for DNA methylation analysis of differential methylated regions from imprinted genes

The bisulfite genomic sequencing method is one of the most widely used techniques for methylation analysis in heterogeneous unbiased PCR, amplifying for both methylated and unmethylated alleles simultaneously. However, it requires labor-intensive and time-consuming cloning and sequencing steps. In the current study, we used a denaturing high-performance liquid chromatography (DHPLC) procedure in a complementary way with the bisulfite genomic sequencing to analyze the methylation of differentially methylated regions (DMRs) of imprinted genes. We showed reliable and reproducible results in distinguishing overall methylation profiles of DMRs regions of human SNRPN, H19, MEST/PEG1, LIT1, IGF2, TSSC5, WT1 antisense, and mouse H19, Mest/Peg1, Igf2R imprinted genes. These DHPLC profiles were in accordance with bisulfite genomic sequencing data and may serve as a type of "fingerprint," revealing the overall methylation status of DMRs associated with sample heterogeneity. We conclude that DHPLC analysis could be used to increase the throughput efficiency of methylation pattern analysis of imprinted genes after the bisulfite conversion of genomic DNA and unbiased PCR amplification.     

A simple algorithm for quantifying DNA methylation levels on multiple independent CpG sites in bisulfite genomic sequencing electropherograms

DNA methylation at cytosines is a widely studied epigenetic modification. Methylation is commonly detected using bisulfite modification of DNA followed by PCR and additional techniques such as restriction digestion or sequencing. These additional techniques are either laborious, require specialized equipment, or are not quantitative. Here we describe a simple algorithm that yields quantitative results from analysis of conventional four-dye-trace sequencing. We call this method Mquant and we compare it with the established laboratory method of combined bisulfite restriction assay (COBRA). This analysis of sequencing electropherograms provides a simple, easily applied method to quantify DNA methylation at specific CpG sites.     

Rolling circle DNA replication by extracts of herpes simplex virus type 1-infected human cells.

Whole-cell extracts of herpes simplex virus type 1-infected human cells (293 cells) can promote the rolling circle replication of circular duplex DNA molecules. The products of the reaction are longer than monomer unit length and are the result of semiconservative DNA replication by the following criteria: (i) resistance to DpnI and susceptibility to MboI restriction enzymes, (ii) shift in density on a CsCl gradient of the products synthesized in the presence of bromo-dUTP to a position on the gradient consistent with those of molecules composed mainly of one parental DNA strand and one newly synthesized DNA strand, and (iii) the appearance in the electron microscope of molecules consisting of duplex circles with multiunit linear appendages, a characteristic of a rolling circle mode of DNA replication. The reaction requires ATP and is dependent on herpes simplex virus type 1-encoded DNA polymerase.     

Microsoft Word macro for analysis of cytosine methylation by the bisulfite deamination reaction

Cytosine methylation at CpG dinucleotides is an important control mechanism in development, differentiation, and neoplasia. Bisulfite genomic sequencing and its modifications have been developed to examine methylation at these CpG dinucleotides. To use these methods, one has to (i) manually convert the sequence to that produced by bisulfite conversion and PCR amplification, taking into account that cytosine residues at CpG dinucleotides may or may not be converted depending on their methylation status, (ii) identify relevant restriction sites that may be used for methylation analysis, and (iii) conduct similar steps with the other DNA strand since the two strands of DNA are no longer complementary after bisulfite conversion. To automate these steps, we have developed a macro that can be used with Microsoft Word. This macro (i) converts genomic sequence to modified sequence that would result after bisulfite treatment facilitating primer design for bisulfite genomic sequencing and methylation-sensitive PCR assay and (ii) identifies restriction sites that are preserved in bisulfite-converted and PCR-amplified product only if cytosine residues at relevant CpG dinucleotides are methylated (and thereby not converted to uracil) in the genomic DNA.     

DNA forms of the geminivirus African cassava mosaic virus consistent with a rolling circle mechanism of replication.

We have analysed DNA from African cassava mosaic virus (ACMV)-infected Nicotiana benthamiana by two-dimensional agarose gel electrophoresis and detected ACMV-specific DNAs by blot-hybridisation. ACMV DNA forms including the previously characterised single-stranded, open-circular, linear and supercoiled DNAs along with five previously uncharacterised heterogeneous DNAs (H1-H5) were resolved. The heterogeneous DNAs were characterised by their chromatographic properties on BND-cellulose and their ability to hybridise to strand-specific and double-stranded probes. The data suggest a rolling circle mechanism of DNA replication, based on the sizes and strand specificity of the heterogeneous single-stranded DNA forms and their electrophoretic properties in relation to genome length single-stranded DNAs. Second-strand synthesis on a single-stranded virus-sense template is evident from the position of heterogeneous subgenomic complementary-sense DNA (H3) associated with genome-length virus-sense template (VT) DNA. The position of heterogeneous virus-sense DNA (H5), ranging in size from one to two genome lengths, is consistent with its association with genome-length complementary-sense template (CT) DNA, reflecting virus-sense strand displacement during replication from a double-stranded intermediate. The absence of subgenomic complementary-sense DNA associated with the displaced virus-sense strand suggests that replication proceeds via an obligate single-stranded intermediate. The other species of heterogeneous DNAs comprised concatemeric single-stranded virus-sense DNA (H4), and double-stranded or partially single-stranded DNA (H1 and H2).     

Rapid analysis of DNA methylation using new restriction enzyme sites created by bisulfite modification.

Bisulfite converts non-methylated cytosine in DNA to uracil leaving 5-methylcytosine unaltered. Here, predicted changes in restriction enzyme sites following reaction of genomic DNA with bisulfite and amplification of the product by the polymerase chain reaction (PCR) were used to assess the methylation of CpG sites. This procedure differs from conventional DNA methylation analysis by methylation-sensitive restriction enzymes because it does not rely on an absence of cleavage to detect methylated sites, the two strands of DNA produce different restriction enzyme sites and may be differentially analyzed, and closely related sequences may be separately analyzed by using specific PCR primers.     

Genome-wide analysis of DNA methylation patterns

Cytosine methylation is the most common covalent modification of DNA in eukaryotes. DNA methylation has an important role in many aspects of biology, including development and disease. Methylation can be detected using bisulfite conversion, methylation-sensitive restriction enzymes, methyl-binding proteins and anti-methylcytosine antibodies. Combining these techniques with DNA microarrays and high-throughput sequencing has made the mapping of DNA methylation feasible on a genome-wide scale. Here we discuss recent developments and future directions for identifying and mapping methylation, in an effort to help colleagues to identify the approaches that best serve their research interests.     

利用甲基化特异性引物高通量检测DNA甲基化

建立一种基于甲基化特异性引物和SAGE技术的高通量DNA甲基化定量检测新方法（MSP-SAGE）,首先利用亚硫酸氢钠对基因组DNA进行处理,使未甲基化的C转变为U,而甲基化的CpG不变.将处理和未处理的DNA双链变性后用随机引物PNNNNCG对存在含有CG的单链进行延伸,而无甲基化CG的单链处则不能延伸;将差异延伸的单链序列和频次信息经过系列分子操作后,引入PCR扩增模板;对中间带有未知序列的PCR扩增产物进行串连克隆测序.将来自于未处理组和处理组的某一CpG位点的序列出现的次数定义为［Tags］A和［Tags］B,将标准系列的实际甲基化水平和［Tags］B/［Tags］A之间建立线性回归方程.根据每一CpG位点的［Tags］B/［Tags］A比值可反推该位点的甲基化水平.MSP-SAGE具有良好的线性,基于标准系列的［Tags］B/［Tags］A与其实际甲基化水平的标准曲线方程为y＝1.455x（R²＝0.984,P<0.01）.MSP-SAGE的回收率在95%到110%之间,精确度位于4.2％和10.5%,检测限在3％左右,单次检测通量可达24个CpG位点.MSP-SAGE是一种很有应用前途的高通量DNA甲基化定量检测方法.     

Formamide as a denaturant for bisulfite conversion of genomic DNA: Bisulfite sequencing of the GSTPi and RARβ2 genes of 43 formalin-fixed paraffin-embedded prostate cancer specimens

Analysis of methylated DNA, which refers to 5-methycytosine (5mC) versus cytosine (C) at specific loci in genomic DNA (gDNA), has received increased attention in epigenomics, particularly in the area of cancer biomarkers. Many different methods for analysis of methylated DNA rely on initial reaction of gDNA with concentrated acidic sodium bisulfite to quantitatively convert C to uracil (U) via sulfonation of denatured, single-stranded gDNA under conditions where 5mC is resistant to analogous sulfonation leading to thymine (T). These methods typically employ polymerase chain reaction (PCR) amplification after bisulfite conversion, thereby leading to readily detectable amounts of amplicons where T and C are measured as surrogates for C and 5mC in the original unconverted gDNA. However, incomplete bisulfite conversion of C in gDNA has been reported to be a common source of error in analysis of methylated DNA. Incomplete conversion can be revealed during the course of bisulfite sequencing, which is the generally accepted “gold standard” for analysis of methylated DNA. Previous bisulfite sequencing investigations of conventional predenaturation of gDNA with NaOH followed by the use of bisulfite containing added urea to maintain denaturation and thus mitigate incomplete conversion of C have been reported to give conflicting results. The current study describes a new approach where conventional predenaturation of gDNA with NaOH is instead achieved with formamide and maintains denaturation during subsequent sample handling and sulfonation. This formamide-based method was applied to 46 formalin-fixed/paraffin-embedded (FFPE) biopsy tissue specimens from well-characterized patients with primary prostate cancer. These specimens were representative of difficult-to-analyze samples due to the chemically compromised nature of the gDNA, which was recovered by modifying the protocol for a commercially available total RNA/DNA extraction kit (RecoverALL). An additional novel aspect of this study was analysis of CpG-rich promoter regions of two prostate cancer-related genes: glutathione S-transferase pi (GSTPi) and retinoic acid receptor beta2 (RARβ2). High-quality bisulfite sequencing results were obtained for both genes in 43 of 46 (93%) specimens. Detection of methylated GSTPi and RARβ2 genes was significantly associated with primary prostate cancer as compared with the benign prostate (Fisher’s exact test, P < 0.001). The sensitivity and specificity of detection of methylated GSTPi and RARβ2 genes were 86% and 100% and 91% and 100%, respectively. Moreover, the presence of either methylated gene was detected in primary prostate cancer with sensitivity and specificity of 100% and 100%, respectively. The results demonstrated a high degree of reliability of formamide-based denaturation and bisulfite conversion that should extend, generally, to FFPE and other types of samples intended for any analytical method predicated on bisulfite conversion. This pilot study also demonstrated the efficacy of determining methylation of these two genes with high sensitivity and specificity in FFPE biopsy tissue specimens. Moreover, the results showed a highly significant association of methylated GSTPi and RARβ2 genes with primary prostate cancer. Finally, this improved procedure for determining these two methylated genes may allow the detection of prostate cancer cells in core biopsy specimens with insufficient numbers of cells and poor morphology.     

Maize chromomethylase Zea methyltransferase2 is required for CpNpG methylation

A cytosine DNA methyltransferase containing a chromodomain, Zea methyltransferase2 (Zmet2), was cloned from maize. The sequence of ZMET2 is similar to that of the Arabidopsis chromomethylases CMT1 and CMT3, with C-terminal motifs characteristic of eukaryotic and prokaryotic DNA methyltransferases. We used a reverse genetics approach to determine the function of the Zmet2 gene. Plants homozygous for a Mutator transposable element insertion into motif IX had a 13% reduction in methylated cytosines. DNA gel blot analysis of these plants with methylation-sensitive restriction enzymes and bisulfite sequencing of a 180-bp knob sequence showed reduced methylation only at CpNpG sites. No reductions in methylation were observed at CpG or asymmetric sites in heterozygous or homozygous mutant plants. Our research shows that chromomethylase Zmet2 is required for in vivo methylation of CpNpG sequences.     

CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing is widely used to measure cytosine methylation at single-base resolution in eukaryotes. It permits systems-level analysis of genomic methylation patterns associated with gene expression and chromatin structure. However, methods for large-scale identification of methylation patterns from bisulfite sequencing are lacking. We developed a comprehensive tool, CpG_MPs, for identification and analysis of the methylation patterns of genomic regions from bisulfite sequencing data. CpG_MPs first normalizes bisulfite sequencing reads into methylation level of CpGs. Then it identifies unmethylated and methylated regions using the methylation status of neighboring CpGs by hotspot extension algorithm without knowledge of pre-defined regions. Furthermore, the conservatively and differentially methylated regions across paired or multiple samples (cells or tissues) are identified by combining a combinatorial algorithm with Shannon entropy. CpG_MPs identified large amounts of genomic regions with different methylation patterns across five human bisulfite sequencing data during cellular differentiation. Different sequence features and significantly cell-specific methylation patterns were observed. These potentially functional regions form candidate regions for functional analysis of DNA methylation during cellular differentiation. CpG_MPs is the first user-friendly tool for identifying methylation patterns of genomic regions from bisulfite sequencing data, permitting further investigation of the biological functions of genome-scale methylation patterns.     

Evidence for cytosine methylation of non-symmetrical sequences in transgenic Petunia hybrida.

A considerable proportion of cytosine residues in plants are methylated at carbon 5. According to a well-accepted rule, cytosine methylation is confined to symmetrical sequences such as CpG and CpNpG, which provide the signal for faithful transmission of symmetrical methylation patterns by maintenance methylase. Using a genomic sequencing technique, we have analysed cytosine methylation patterns within a hypermethylated and a hypomethylated state of a transgene in Petunia hybrida. Examination of a part of the transgene promoter revealed that in both states m5C residues located within non-symmetrical sequences could be detected. Non-symmetrical C residues in the two states were methylated at frequencies of 5.9 and 31.9%, respectively. Methylation appeared to be distributed heterogeneously, but some DNA regions were more intensively methylated than others. Our results show that at least in a transgene, a heterogeneous methylation pattern, which does not depend on symmetry of target sequences, can be established and conserved.     

Solution structure of the endonuclease domain from the master replication initiator protein of the nanovirus faba bean necrotic yellows virus and comparison with the corresponding geminivirus and circovirus structures

Nanoviruses are a family of plant viruses that possess a genome of multiple circular single-stranded DNA (ssDNA) components and are strikingly similar in their replication mode to the plant geminiviruses and to the circoviruses that infect birds or mammals. These viruses multiply by rolling circle replication using virus-encoded multifunctional replication initiator proteins (Rep proteins) that catalyze the initiation of replication on a double-stranded DNA (dsDNA) intermediate and the resolution of the ssDNA into circles. Here we report the solution NMR three-dimensional structure of the endonuclease domain from the master Rep (M-Rep) protein of faba bean necrotic yellows virus (FBNYV), a representative of the nanoviruses. The domain comprises amino acids 2-95 (M-Rep2-95), and its global fold is similar to those previously described for the gemini- and circovirus Rep endonuclease domains, consisting of a central 5-stranded antiparallel beta-sheet covered on one side by an alpha-helix and irregular loops and on the other, more open side of the domain, by an alpha-helix containing the catalytic tyrosine residue (the catalytic helix). Longer domain constructs extending to amino acids 117 and 124 were also characterized. They contain an additional alpha-helix, are monomeric, and exhibit catalytic activity indistinguishable from that of M-Rep2-95. The binding site for the catalytic metal was identified by paramagnetic broadening and maps to residues on the exposed face of the central beta-sheet. A comparison with the previously determined Rep endonuclease domain structures of tomato yellow leaf curl Sardinia virus (TYLCSV), a geminivirus, and that of porcine circovirus type 2 (PCV2) Rep allows the identification of a positively charged surface that is most likely involved in dsDNA binding, and reveals common features shared by all endonuclease domains of nanovirus, geminivirus, and circovirus Rep proteins.     

Methylation assay by nucleotide incorporation: a quantitative assay for regional CpG methylation density

Although the aberrant methylation in CpG islands is of great interest as a causative role in human malignancies, it has been very difficult to accurately determine methylation density. Here we report a novel microplate-based quantitative methylation assay, designated MANIC, for a region containing a number of CpG sites based on incorporation of hapten-labeled dCTP at cytosine sites where the methylated cytosines have not been converted to uracil by the bisulfite treatment. Validation using control DNAs revealed that the method was sensitive enough to detect < 1.25% methylated DNA and that calibration curve was linear. With this approach, we determined relative methylation density of O6-methylguanine-DNA methyltransferase gene promoter containing 12 CpG sites among the 12 colorectal cancers and corresponding normal mucosal tissues. Consequently, MANIC showed a high concordance with results by a quantitative method, bisulfite PCR single-stranded conformational polymorphism (BiPS). MANIC is a technique that avoids cumbersome procedures such as electrophoresis or the use of radiolabeling and is applicable to any sequence regardless of the total number of CpG sites or heterogeneity in methylation status.