Assessment of RainDrop BS-seq as a method for large-scale,targeted bisulfite sequencing

We present a systematic assessment of RainDrop BS-seq, a novel method for large-scale, targeted bisulfite sequencing using microdroplet-based PCR amplification coupled with next-generation sequencing. We compared DNA methylation levels at 498 target loci (1001 PCR amplicons) in human whole blood, osteosarcoma cells and an archived tumor tissue sample. We assessed the ability of RainDrop BS-seq to accurately measure DNA methylation over a range of DNA quantities (from 10 to 1500 ng), both with and without whole-genome amplification (WGA) following bisulfite conversion. DNA methylation profiles generated using at least 100 ng correlated well (median R = 0.92) with those generated on Illumina Infinium HumanMethylation450 BeadChips, currently the platform of choice for epigenome-wide association studies (EWAS). WGA allowed for testing of samples with a starting DNA amount of 10 and 50 ng, although a reduced correlation was observed (median R = 0.79). We conclude that RainDrop BS-seq is suitable for measuring DNA methylation levels using nanogram quantities of DNA, and can be used to study candidate epigenetic biomarker loci in an accurate and high-throughput manner, paving the way for its application to routine clinical diagnostics.     

Differential methylation analysis for bisulfite sequencing using DSS

Bisulfite sequencing (BS-seq) technology measures DNA methylation at single nucleotide resolution. A key task in BS-seq data analysis is to identify differentially methylation (DM) under different conditions. Here we provide a tutorial for BS-seq DM analysis using Bioconductor package DSS. DSS uses a beta-binomial model to characterize the sequence counts from BS-seq, and implements rigorous statistical method for hypothesis testing. It provides flexible functionalities for a variety of DM analyses.     

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.     

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.     

Identification of differentially methylated regions using streptavidin bisulfite ligand methylation enrichment (SuBLiME), a new method to enrich for methylated DNA prior to deep bisulfite genomic sequencing

We have developed a method that enriches for methylated cytosines by capturing the fraction of bisulfite-treated DNA with unconverted cytosines. The method, called streptavidin bisulfite ligand methylation enrichment (SuBLiME), involves the specific labeling (using a biotin-labeled nucleotide ligand) of methylated cytosines in bisulfite-converted DNA. This step is then followed by affinity capture, using streptavidin-coupled magnetic beads. SuBLiME is highly adaptable and can be combined with deep sequencing library generation and/or genomic complexity-reduction. In this pilot study, we enriched methylated DNA from Csp6I-cut complexity-reduced genomes of colorectal cancer cell lines (HCT-116, HT-29 and SW-480) and normal blood leukocytes with the aim of discovering colorectal cancer biomarkers. Enriched libraries were sequenced with SOLiD-3 technology. In pairwise comparisons, we scored a total of 1,769 gene loci and 33 miRNA loci as differentially methylated between the cell lines and leukocytes. Of these, 516 loci were differently methylated in at least two promoter-proximal CpG sites over two discrete Csp6I fragments. Identified methylated gene loci were associated with anatomical development, differentiation and cell signaling. The data correlated with good agreement to a number of published colorectal cancer DNA methylation biomarkers and genomic data sets. SuBLiME is effective in the enrichment of methylated nucleic acid and in the detection of known and novel biomarkers.     

Design of efficient simplified genomic DNA and bisulfite sequencing in large plant populations

The next generation sequencing enables generation of high resolution and high throughput data for structure sequence of any genome at a fast declining cost. This opens opportunity for population based genetic and genomic analyses. In many applications, whole genome sequencing or re-sequencing is unnecessary or prohibited by budget limits. The Reduced Representation Genome Sequencing (RRGS), which sequences only a small proportion of the genome of interest, has been proposed to deal with the situations. Several forms of RRGS are proposed and implemented in the literature. When applied to plant or crop species, the current RRGS protocols shared a key drawback that a significantly high proportion (up to 60%) of sequence reads to be generated may be of non-genomic origin but attributed to chloroplast DNA or rRNA genes, leaving an exceptional low efficiency of the sequencing experiment. We recommended and discussed here the design of optimized simplified genomic DNA and bisulfite sequencing strategies, which may greatly improves efficiency of the sequencing experiments by bringing down the presentation of the undesirable sequencing reads to less than 10% in the whole sequence reads. The optimized RAD-seq and RRBS-seq methods are potentially useful for sequence variant screening and genotyping in large plant/crop populations.     

Comparison of bisulfite sequencing PCR with pyrosequencing for measuring differences in DNA methylation

DNA methylation strongly affects chromatin structure and the regulation of gene expression. For many years, bisulfite sequencing PCR (BSP) has served as the “gold standard” for measuring DNA methylation. However, with the evolution of pyrosequencing as a tool to evaluate DNA methylation, the need arises to compare the relative efficiencies of the two techniques in measuring DNA methylation. We provide for the first time a direct assessment of BSP and pyrosequencing to detect and quantify hypomethylation, hypermethylation, and mixed methylation of the ABCB1 promoter in various drug-sensitive and drug-resistant MCF-7 breast cancer cell lines through head-to-head experimentation. Our findings indicate that although both methods can reliably detect increased, decreased, and mixed methylation of DNA, BSP appears to be more sensitive than pyrosequencing at detecting strong hypermethylation of DNA. However, we also observed greater variability in the methylation of CpG sites by BSP, possibly due to the additional bacterial cloning step required by BSP over pyrosequencing. BSP and pyrosequencing equally detected hypomethylation and mixed methylation of DNA. The ability of pyrosequencing to reliably detect differences in DNA methylation across cell populations without requiring the cloning of bisulfite-treated DNA into bacterial expression vectors was seen as a major advantage of this technique.     

A systematic study of normalization methods for Infinium 450K methylation data using whole-genome bisulfite sequencing data

DNA methylation plays an important role in disease etiology. The Illumina Infinium HumanMethylation450 (450K) BeadChip is a widely used platform in large-scale epidemiologic studies. This platform can efficiently and simultaneously measure methylation levels at ∼480,000 CpG sites in the human genome in multiple study samples. Due to the intrinsic chip design of 2 types of chemistry probes, data normalization or preprocessing is a critical step to consider before data analysis. To date, numerous methods and pipelines have been developed for this purpose, and some studies have been conducted to evaluate different methods. However, validation studies have often been limited to a small number of CpG sites to reduce the variability in technical replicates. In this study, we measured methylation on a set of samples using both whole-genome bisulfite sequencing (WGBS) and 450K chips. We used WGBS data as a gold standard of true methylation states in cells to compare the performances of 8 normalization methods for 450K data on a genome-wide scale. Analyses on our dataset indicate that the most effective methods are peak-based correction (PBC) and quantile normalization plus β-mixture quantile normalization (QN.BMIQ). To our knowledge, this is the first study to systematically compare existing normalization methods for Illumina 450K data using novel WGBS data. Our results provide a benchmark reference for the analysis of DNA methylation chip data, particularly in white blood cells.     

Single-molecule polymerase chain reaction reduces bias: application to DNA methylation analysis by bisulfite sequencing

The treatment of DNA with bisulfite, which converts C to U but leaves 5-methyl-C unchanged, forms the basis of many analytical techniques for DNA methylation analysis. Many techniques exist for measuring the methylation state of a single CpG but, for analysis of an entire region, cloning and sequencing remains the gold standard. However, biases in polymerase chain reaction (PCR) amplification and in cloning can skew the results. We hypothesized that single-molecule PCR (smPCR) amplification would eliminate the PCR amplification bias because competition between templates that amplify at different efficiencies no longer exists. The amplified products can be sequenced directly, thus eliminating cloning bias. We demonstrated this accurate and unbiased approach by analyzing a sample that was expected to contain a 50:50 ratio of methylated to unmethylated molecules: a region of the X-linked FMR1 gene from a human female cell line. We compared traditional cloning and sequencing to smPCR and sequencing. Sequencing smPCR products gave an expected methylated to unmethylated ratio of 48:52, whereas conventional cloning and sequencing gave a biased ratio of 72:28. Our results show that smPCR sequencing can eliminate both PCR and cloning bias and represents an attractive approach to bisulfite sequencing.     

DNA甲基化分析中重亚硫酸盐处理DNA转化效率的评估方法

为建立一种评估重亚硫酸盐处理DNA样本后胞嘧啶转化效率的有效方法,以两组不同的TaqMan qPCR检测梯度稀释的重亚硫酸盐处理和未处理的DNA标准品,建立转化与未转化的DNA Ct值以及对应的DNA拷贝数的标准曲线。使用相同的探针定量检测重亚硫酸盐处理后的DNA样本评估转化效率。结果显示该方法应用两组探针,根据相应的标准曲线,精确评估样本经重亚硫酸盐处理的转化效率。使用已知转化和未转化拷贝数的混合DNA作为模板,证实了该方法的可靠性。同时也对不同重亚硫酸盐试剂盒处理DNA的转化效率进行了评估,结果显示,该方法能够有效地评估DNA样品重亚硫酸盐的转化效率,为DNA甲基化准确分析提供了可靠快捷的方法。     

Methylation status of Sillago japonica satellite DNA examined by bisulfite modification

A member of Sillago japonica satellite DNA contained internal subrepeats in its 174 bp unit. S. Japonica genomic DNA isolated from liver tissue was subjected to bisulfite modification, and the DNA sequences of about 40 bp flanked by both subrepeats were amplified by polymerase chain reaction (PCR). This protocol, combination of bisulfite reaction and PCR, converts cytosines in the genomic DNA to thymines in the amplified DNA, whereas 5-methylcytosines in the genomic DNA remain as cytosines. Sequence analysis of the amplified DNA fragments revealed that most of the cytosine residues at CpG were methylated in this region.     

Q-RRBS: a quantitative reduced representation bisulfite sequencing method for single-cell methylome analyses

Reduced representation bisulfite sequencing (RRBS) is a powerful method of DNA methylome profiling that can be applied to single cells. However, no previous report has described how PCR-based duplication-induced artifacts affect the accuracy of this method when measuring DNA methylation levels. For quantifying the effects of duplication-induced artifacts on methylome profiling when using ultra-trace amounts of starting material, we developed a novel method, namely quantitative RRBS (Q-RRBS), in which PCR-induced duplication is excluded through the use of unique molecular identifiers (UMIs). By performing Q-RRBS on varying amounts of starting material, we determined that duplication-induced artifacts were more severe when small quantities of the starting material were used. However, through using the UMIs, we successfully eliminated these artifacts. In addition, Q-RRBS could accurately detect allele-specific methylation in absence of allele-specific genetic variants. Our results demonstrate that Q-RRBS is an optimal strategy for DNA methylation profiling of single cells or samples containing ultra-trace amounts of cells.     

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.     

An Optimized Method for the Construction of a DNA Methylome from Small Quantities of Tissue or Purified DNA from Arabidopsis Embryo

DNA methylation is an important epigenetic mechanism affecting genome structure, gene regulation, and the silencing of transposable elements. Cell- and tissue-specific methylation patterns are critical for differentiation and development in eukaryotes. Dynamic spatiotemporal methylation data in these cells or tissues is, therefore, of great interest. However, the construction of bisulfite sequencing libraries can be challenging if the starting material is limited or the genome size is small, such as in Arabidopsis. Here, we describe detailed methods for the purification of Arabidopsis embryos at all stages, and the construction of comprehensive bisulfite libraries from small quantities of input. We constructed bisulfite libraries by releasing embryos from intact seeds, using a different approach for each developmental stage, and manually picking single-embryo with microcapillaries. From these libraries, reliable Arabidopsis methylome data were collected allowing, on average, 11-fold coverage of the genome using as few as five globular, heart, and torpedo embryos as raw input material without the need for DNA purification step. On the other hand, purified DNA from as few as eight bending torpedo embryos or a single mature embryo is sufficient for library construction when RNase A is treated before DNA extraction. This method can be broadly applied to cells from different tissues or cells from other model organisms. Methylome construction can be achieved using a minimal amount of input material using our method; thereby, it has the potential to increase our understanding of dynamic spatiotemporal methylation patterns in model organisms.     

Cross-platform compatibility of Hi-Plex,a streamlined approach for targeted massively parallel sequencing

Although per-base sequencing costs have decreased during recent years, library preparation for targeted massively parallel sequencing remains constrained by high reagent cost, limited design flexibility, and protocol complexity. To address these limitations, we previously developed Hi-Plex, a polymerase chain reaction (PCR) massively parallel sequencing strategy for screening panels of genomic target regions. Here, we demonstrate that Hi-Plex applied with hybrid adapters can generate a library suitable for sequencing with both the Ion Torrent and the TruSeq chemistries and that adjusting primer concentrations improves coverage uniformity. These results expand Hi-Plex capabilities as an accurate, affordable, flexible, and rapid approach for various genetic screening applications.