Whole methylome analysis by ultra-deep sequencing using two-base encoding.

Methylation, the addition of methyl groups to cytosine (C), plays an important role in the regulation of gene expression in both normal and dysfunctional cells. During bisulfite conversion and subsequent PCR amplification, unmethylated Cs are converted into thymine (T), while methylated Cs will not be converted. Sequencing of this bisulfite-treated DNA permits the detection of methylation at specific sites. Through the introduction of next-generation sequencing technologies (NGS) simultaneous analysis of methylation motifs in multiple regions provides the opportunity for hypothesis-free study of the entire methylome. Here we present a whole methylome sequencing study that compares two different bisulfite conversion methods (in solution versus in gel), utilizing the high throughput of the SOLiD System. Advantages and disadvantages of the two different bisulfite conversion methods for constructing sequencing libraries are discussed. Furthermore, the application of the SOLiD bisulfite sequencing to larger and more complex genomes is shown with preliminary in silico created bisulfite converted reads.

Formamide as a denaturant for bisulfite conversion of genomic DNA: Bisulfite sequencing of the GSTPi and RARbeta2 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 (RARbeta2). High-quality bisulfite sequencing results were obtained for both genes in 43 of 46 (93%) specimens. Detection of methylated GSTPi and RARbeta2 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 RARbeta2 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 RARbeta2 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.

Preparation of genome-wide DNA fragment libraries using bisulfite in polyacrylamide gel electrophoresis slices with formamide denaturation and quality control for massively parallel sequencing by oligonucleotide ligation and detection.

Bisulfite sequencing is widely used for analysis of DNA methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich or other loci in genomic DNA (gDNA). Such methods typically involve reaction of gDNA with bisulfite followed by polymerase chain reaction (PCR) amplification of specific regions of interest that, overall, converts C-->T (thymine) and 5mC-->C and then capillary sequencing to measure C versus T composition at CpG sites. Massively parallel sequencing by oligonucleotide ligation and detection (SOLiD) has recently enabled relatively low-cost whole genome sequencing, and it would be highly desirable to apply such massively parallel sequencing to bisulfite-converted whole genomes to determine DNA methylation status of an entire genome, which has heretofore not been reported. As an initial step toward achieving this goal, we have extended our ongoing interest in improving bisulfite conversion sample preparation to include a human genome-wide fragment library for SOliD. The current article features novel use of formamide denaturant during bisulfite conversion of a suitably constructed library directly in a band slice from polyacryamide gel electrophoresis (PAGE). To validate this new protocol for 5mC-protected fragment library conversion, which we refer to as Bis-PAGE, capillary-based size analysis and Sanger sequencing were carried out for individual amplicons derived from single-molecule PCR (smPCR) of randomly selected library fragments. smPCR/Capillary Sanger sequencing of approximately 200 amplicons unambiguously demonstrated greater than 99% C-->T conversion. All of these approximately 200 Sanger sequences were analyzed with a previously published web-accessible bioinformatics tool (methBLAST) for mapping to human chromosomes, the results of which indicated random distribution of analyzed fragments across all chromosomes. Although these particular Bis-PAGE conversion and quality control methods were exemplified in the context of a fragment library for SOLiD, the concepts can be generalized to include other genome-wide library constructions intended for DNA methylation analysis by alternative high-throughput or massively parallelized methods that are currently available.

Capillary electrophoretic analysis of methylation status in CpG-rich regions by single-base extension of primers modified with N6-methoxy-2,6-diaminopurine.

Polymerase-mediated single-base extension (SBE) of primers using a fluorescently labeled 2',3'-dideoxynucleotide triphosphate terminator was originally commercialized as SNaPshot for analysis of single-nucleotide polymorphisms by capillary electrophoresis (CE). Application of this general method to bisulfite-converted/PCR-amplified genomic DNA (gDNA) to analytically infer polymorphic methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich regions of gDNA has been noted previously by others to be limited by CE mobility-shifted peaks for SBE products derived from guanine (G)/adenine (A)-mixed-base primers used to hybridize to possible polymorphic sites (i.e., C vs. thymine [T] resulting from 5mC vs. C, respectively). This limitation is precluded in the current study by using novel SNaPshot primers modified with N(6)-methoxy-2,6-diaminopurine (K), which was originally described in 1991 by Brown and Lin as a unique adenine-guanine analog capable of participating in three H-bonds with C or T in DNA. Oligonucleotides modified by K as a bispecific complement for C/T are commercially available or can be readily synthesized, and they may have utility in other assay formats used to analyze CpG methylation status.

Methylation-dependent fragment separation: novel analysis of 5-methyl cytosine by capillary electrophoresis of amplified DNA using PCR incorporation of chemically modified dCTP.

Methylation of the cytosine (C) ring to form 5-methyl cytosine (MeC) in normally unmethylated CpG-rich regions of promoters in genes is associated with transcriptional silencing. Quantification of MeC is of current interest in findining new biomarkers for cancer. To this end, and for basic research in epigenomics, we have investigated a new method for relatively simple measurement of MeC content by capillary electrophoresis (CE). PCR amplicons for CE analysis are generated from bisulfite-converted DNA [C --> uracil (U)] using fluorescently labeled primers that anneal independent of methylation status. Resultant incorporation of C vs. T at original MeC vs. C positions can lead to separate CE peaks for signal integration that is proportional to MeC content. Furthermore, these PCR products are suitable for additional methylation analyses by sequencing, single-base extension, or TaqMan. Interestingly, PCR using alpha -thio-dCTP led to greater CE separations.

Methylation-dependent fragment separation: direct detection of DNA methylation by capillary electrophoresis of PCR products from bisulfite-converted genomic DNA.

Fundamental to understanding the role of cytosine (C) methylation in genomic DNA (gDNA) is the need for robust analysis methods to determine the location and degree of this modification. We report a novel method for methylation detection by denaturing capillary electrophoresis (CE) using standard fragment analysis conditions. Bisulfite treatment of gDNA will selectively deaminate C but not 5-methylcytosine (5mC). Amplicons generated from bisulfite-converted gDNA are analyzed immediately after PCR using a 6-carboxy fluorescein (6-FAM) dye-labeled primer. The amplicons from methylated and unmethylated gDNA separate based solely on base composition due to the presence of multiple C versus thymine (T) differences. By direct detection of PCR amplicons following PCR using primers that anneal independent of methylation status, the overall workflow from gDNA sample input to data analysis is relatively simple. Furthermore, the same PCR product is suitable for additional analyses such as direct sequencing, cloning and sequencing, single-base extension, and post-PCR incorporation of a modified dCTP, the latter of which allows resolution of amplicons with as little as a single C/T difference. We show the utility of this novel CE detection assay by analyzing the hypermethylated region of the fragile-X FMR1 locus.

Assessing gene expression variation in normal human tissues using GeneTag, a novel, global, sensitive profiling method.

GeneTag is a novel expression profiling method that allows the visualization, quantification and identification of expressed genes-whether known or novel-in any species, tissue or cell type, independent of knowledge of the underlying sequence. Here we describe the application of this method to determine variation of gene expression in individual human liver samples and the identification of tissue-specific genes by comparing expression patterns across several human organs. Expression data are stored in a database for future reference and data analysis relies on proprietary software, which allows complex comparisons to be performed. Differentially expressed genes are quickly identified through a link to a sequence database. The results from our study underscore the importance of knowledge of individual variation of gene expression for the design and interpretation of transcript profiling experiments in the context of any biological question.