Systematic assessment of reduced representation bisulfite sequencing to human blood samples: A promising method for large-sample-scale epigenomic studies

Complementary to the time- and cost-intensive direct bisulfite sequencing, we applied reduced representation bisulfite sequencing (RRBS) to the human peripheral blood mononuclear cells (PBMC) from YH, the Asian individual whose genome and epigenome has been deciphered in the YH project and systematically assessed the genomic coverage, coverage depth and reproducibility of this technology as well as the concordance of DNA methylation levels measured by RRBS and direct bisulfite sequencing for the detected CpG sites. Our result suggests that RRBS can cover more than half of CpG islands and promoter regions with a good coverage depth and the proportion of the CpG sites covered by the biological replicates reaches 80-90%, indicating good reproducibility. Given a smaller data quantity, RRBS enjoys much better coverage depth than direct bisulfite sequencing and the concordance of DNA methylation levels between the two methods is high. It can be concluded that RRBS is a time and cost-effective sequencing method for unbiased DNA methylation profiling of CpG islands and promoter regions in a genome-wide scale and it is the method of choice to assay certain genomic regions for multiple samples in a rapid way.     

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.     

MethPrimer: designing primers for methylation PCRs

MOTIVATION: DNA methylation is an epigenetic mechanism of gene regulation. Bisulfite- conversion-based PCR methods, such as bisulfite sequencing PCR (BSP) and methylation specific PCR (MSP), remain the most commonly used techniques for methylation mapping. Existing primer design programs developed for standard PCR cannot handle primer design for bisulfite-conversion-based PCRs due to changes in DNA sequence context caused by bisulfite treatment and many special constraints both on the primers and the region to be amplified for such experiments. Therefore, the present study was designed to develop a program for such applications. RESULTS: MethPrimer, based on Primer 3, is a program for designing PCR primers for methylation mapping. It first takes a DNA sequence as its input and searches the sequence for potential CpG islands. Primers are then picked around the predicted CpG islands or around regions specified by users. MethPrimer can design primers for BSP and MSP. Results of primer selection are delivered through a web browser in text and in graphic view.     

Regionally Specific and Genome-Wide Analyses Conclusively Demonstrate the Absence of CpG Methylation in Human Mitochondrial DNA

Although CpG methylation clearly distributes genome-wide in vertebrate nuclear DNA, the state of methylation in the vertebrate mitochondrial genome has been unclear. Several recent reports using immunoprecipitation, mass spectrometry, and enzyme-linked immunosorbent assay methods concluded that human mitochondrial DNA (mtDNA) has much more than the 2 to 5% CpG methylation previously estimated. However, these methods do not provide information as to the sites or frequency of methylation at each CpG site. Here, we have used the more definitive bisulfite genomic sequencing method to examine CpG methylation in HCT116 human cells and primary human cells to independently answer these two questions. We found no evidence of CpG methylation at a biologically significant level in these regions of the human mitochondrial genome. Furthermore, unbiased next-generation sequencing of sodium bisulfite treated total DNA from HCT116 cells and analysis of genome-wide sodium bisulfite sequencing data sets from several other DNA sources confirmed this absence of CpG methylation in mtDNA. Based on our findings using regionally specific and genome-wide approaches with multiple human cell sources, we can definitively conclude that CpG methylation is absent in mtDNA. It is highly unlikely that CpG methylation plays any role in direct control of mitochondrial function.     

Improved reduced representation bisulfite sequencing for epigenomic profiling of clinical samples

Background

DNA methylation plays crucial roles in epigenetic gene regulation in normal development and disease pathogenesis. Efficient and accurate quantification of DNA methylation at single base resolution can greatly advance the knowledge of disease mechanisms and be used to identify potential biomarkers. We developed an improved pipeline based on reduced representation bisulfite sequencing (RRBS) for cost-effective genome-wide quantification of DNA methylation at single base resolution. A selection of two restriction enzymes (Taq^αI and MspI) enables a more unbiased coverage of genomic regions of different CpG densities. We further developed a highly automated software package to analyze bisulfite sequencing results from the Solexa GAIIx system.

Results

With two sequencing lanes, we were able to quantify ~1.8 million individual CpG sites at a minimum sequencing depth of 10. Overall, about 76.7% of CpG islands, 54.9% of CpG island shores and 52.2% of core promoters in the human genome were covered with at least 3 CpG sites per region.

Conclusions

With this new pipeline, it is now possible to perform whole-genome DNA methylation analysis at single base resolution for a large number of samples for understanding how DNA methylation and its changes are involved in development, differentiation, and disease pathogenesis.     

Reduced representation bisulfite sequencing design for assessing the methylation of human CpG islands in large samples

The reduced representation bisulfite sequencing (RRBS) method has been developed for the high-throughput analysis of DNA methylation based on the sequencing of genomic libraries treated with sodium bisulfite by next-generation approaches. In contrast to whole-genome sequencing, the RRBS approach elaborates specific endonucleases to prepare libraries in order to produce pools of CpG-rich DNA fragments. The original RRBS technology based on the use of the MspI libraries allows one to increase the relative number of CpG islands in the pools of genomic fragments compared to whole-genome bisulfite sequencing. Nevertheless, this technology is rarely used due to the high cost compared with bisulfite methylation analysis with hybridization microarrays and significant residual amount of data represented by the sequences of genomic repeats that complicates the alignment and is not of particular interest for developing DNA methylation markers, which is often the main goal of biomedical research. We have developed an algorithm for estimating the likelihood that recognition sites of restriction endonucleases will be represented in CpG islands and present a method of reducing the effective size of the RRBS library without a significant loss of the CpG islands based on the use of the XmaI endonuclease for library preparation. In silico analysis demonstrates that the optimum range of the XmaI-RRBS fragment lengths is 110–200 base pairs. The sequencing of this library allows one to assess the methylation status of over 125000 CpG dinucleotides, of which over 90000 belong to CpG islands.

     

COHCAP: an integrative genomic pipeline for single-nucleotide resolution DNA methylation analysis

COHCAP (City of Hope CpG Island Analysis Pipeline) is an algorithm to analyze single-nucleotide resolution DNA methylation data produced by either an Illumina methylation array or targeted bisulfite sequencing. The goal of the COHCAP algorithm is to identify CpG islands that show a consistent pattern of methylation among CpG sites. COHCAP is currently the only DNA methylation package that provides integration with gene expression data to identify a subset of CpG islands that are most likely to regulate downstream gene expression, and it can generate lists of differentially methylated CpG islands with ∼50% concordance with gene expression from both cell line data and heterogeneous patient data. For example, this article describes known breast cancer biomarkers (such as estrogen receptor) with a negative correlation between DNA methylation and gene expression. COHCAP also provides visualization for quality control metrics, regions of differential methylation and correlation between methylation and gene expression. This software is freely available at https://sourceforge.net/projects/cohcap/.     

High density DNA methylation array with single CpG site resolution

We have developed a new generation of genome-wide DNA methylation BeadChip which allows high-throughput methylation profiling of the human genome. The new high density BeadChip can assay over 480K CpG sites and analyze twelve samples in parallel. The innovative content includes coverage of 99% of RefSeq genes with multiple probes per gene, 96% of CpG islands from the UCSC database, CpG island shores and additional content selected from whole-genome bisulfite sequencing data and input from DNA methylation experts. The well-characterized Infinium® Assay is used for analysis of CpG methylation using bisulfite-converted genomic DNA. We applied this technology to analyze DNA methylation in normal and tumor DNA samples and compared results with whole-genome bisulfite sequencing (WGBS) data obtained for the same samples. Highly comparable DNA methylation profiles were generated by the array and sequencing methods (average R² of 0.95). The ability to determine genome-wide methylation patterns will rapidly advance methylation research.     

Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

CpG PatternFinder: a Windows-based utility program for easy and rapid identification of the CpG methylation status of DNA

The bisulfite genomic sequencing technique is one of the most widely used techniques to study sequence-specific DNA methylation because of its unambiguous ability to reveal DNA methylation status to the order of a single nucleotide. One characteristic feature of the bisulfite genomic sequencing technique is that a number of sample sequence files will be produced from a single DNA sample. The PCR products of bisulfite-treated DNA samples cannot be sequenced directly because they are heterogeneous in nature; therefore they should be cloned into suitable plasmids and then sequenced. This procedure generates an enormous number of sample DNA sequence files as well as adding extra bases belonging to the plasmids to the sequence, which will cause problems in the final sequence comparison. Finding the methylation status for each CpG in each sample sequence is not an easy job. As a result CpG PatternFinder was developed for this purpose. The main functions of the CpG PatternFinder are: (i) to analyze the reference sequence to obtain CpG and non-CpG-C residue position information. (ii) To tailor sample sequence files (delete insertions and mark deletions from the sample sequence files) based on a configuration of ClustalW multiple alignment. (iii) To align sample sequence files with a reference file to obtain bisulfite conversion efficiency and CpG methylation status. And, (iv) to produce graphics, highlighted aligned sequence text and a summary report which can be easily exported to Microsoft Office suite. CpG PatternFinder is designed to operate cooperatively with BioEdit, a freeware on the internet. It can handle up to 100 files of sample DNA sequences simultaneously, and the total CpG pattern analysis process can be finished in minutes. CpG PatternFinder is an ideal software tool for DNA methylation studies to determine the differential methylation pattern in a large number of individuals in a population. Previously we developed the CpG Analyzer program; CpG PatternFinder is our further effort to create software tools for DNA methylation studies.     

Plasma total homocysteine is associated with DNA methylation in patients with schizophrenia

Schizophrenia (SCZ) is a devastating psychiatric disorder with a median lifetime prevalence rate of 0.7–0.8%. Elevated plasma total homocysteine has been suggested as a risk factor for SCZ, and various biological effects of hyperhomocysteinemia have been proposed to be relevant to the pathophysiology of SCZ. As increased attention is paid to aberrant DNA methylation in SCZ, homocysteine is attracting additional interest as a potential key substance. Homocysteine is formed in the methionine cycle, which is involved in one-carbon methyl group-transfer metabolism, and it acts as a methyl donor when it is converted to S-adenosyl-methionine. To date, no studies have examined the relationship between homocysteine and genome-wide DNA methylation in SCZ. We examined the relationship between plasma total homocysteine and DNA methylation patterns in the peripheral leukocytes of patients with SCZ (n = 42) using a quantitative high-resolution DNA methylation array (485,764 CpG sites). Significant homocysteine-related changes in DNA methylation were observed at 1,338 CpG sites that were located across whole gene regions, including promoters, gene bodies and 3′-untranslated regions. Of the 1,338 sites, 758 sites (56.6%) were located in the CpG islands (CGIs) and in the regions flanking CGIs (CGI: 15.8%; CGI shore: 28.2%; CGI shelf: 12.6%), and positive correlations between plasma total homocysteine and DNA methylation were observed predominantly at CpG sites in the CGIs. Our results suggest that homocysteine might play a role in the pathogenesis of SCZ via a molecular mechanism that involves alterations to DNA methylation.