Optimized design and data analysis of tag-based cytosine methylation assays
© Suzuki et al.; licensee BioMed Central Ltd. 2010
Received: 8 January 2010
Accepted: 1 April 2010
Published: 1 April 2010
Using the type III restriction-modification enzyme EcoP15I, we isolated sequences flanking sites digested by the methylation-sensitive HpaII enzyme or its methylation-insensitive MspI isoschizomer for massively parallel sequencing. A novel data transformation allows us to normalise HpaII by MspI counts, resulting in more accurate quantification of methylation at >1.8 million loci in the human genome. This HELP-tagging assay is not sensitive to sequence polymorphism or base composition and allows exploration of both CG-rich and depleted genomic contexts.
Epigenetic mechanisms of transcriptional regulation are increasingly being studied for their potential influences in human disease pathogenesis. Much of this interest is based on the paradigm of neoplastic transformation, in which epigenetic changes appear to be universal, widespread throughout the genome, causative of critical transcriptional changes and predictive of disease prognosis (reviewed in ). Furthermore, these epigenetic changes represent potential pharmacological targets for reversal and amelioration of the disease process .
Of the large number of regulatory processes referred to as epigenetic, there exist numerous assays to study chromatin component distribution, cytosine methylation and microRNA expression genome-wide. The chromatin components include a large number of post-translational modifications of histones, variant histones, DNA-binding proteins and associated complexes, all tested by chromatin immunoprecipitation (ChIP) approaches coupled with microarray hybridization or massively parallel sequencing (MPS). MicroRNAs can be identified and quantified by using microarrays and MPS, while cytosine methylation can be definitively studied by converting the DNA of the genome using sodium bisulfite, shotgun sequencing the product using MPS and mapping this back to the genome to count how frequently cytosines remain unconverted, indicating their methylation in the starting material, due to the resistance of methylcytosine to bisulfite conversion compared with unmethylated cytosines. This allows nucleotide resolution, strand-specific, quantitative assessment of cytosine methylation, with such studies performed in Arabidopsis [3–5] and human cells to date .
While this approach represents the ideal means of testing cytosine methylation, the amount of sequencing necessary (for the human genome, over 1 billion sequences of ~75 bp each ) to generate quantitative information genome-wide remains prohibitive in terms of cost, limiting these studies to the few referred to above. When studying human disease, the emphasis remains on cytosine methylation assays, as it is generally easier to collect clinical samples for DNA purification than for ChIP or even RNA assays. However, the cell populations harvested are rarely of high purity, and we generally do not know the degree of change in cytosine methylation in the disease of interest and thus the quantitative discrimination required for an assay, with some studies to date indicating that the changes may be quite subtle . These concerns emphasize the need for cytosine methylation assays that can detect methylation levels intermediate in value and changes in disease that are relatively modest in magnitude. Certain microarray-based assays to study cytosine methylation have addressed this issue, with the methylated DNA immunoprecipitation (meDIP) assay amenable to such quantification when used for CpG islands  and possibly also for less CG dinucleotide-rich regions . Restriction enzyme-based assays used with microarrays have also proven to be reasonably quantitative, whether based on methylation-sensitive (for example, the HELP assay ) or methylation-dependent (for example, MethylMapper ) enzymes. A promising new MPS-based assay is reduced representation bisulfite sequencing (RRBS), which is designed to study the CG-dense regions defined by short MspI fragments, and provides nucleotide resolution, quantitative data .
The use of MPS for what were previously microarray-based assays has been associated with improved performance , as we found when we modified our HELP (HpaII tiny fragment Enrichment by Ligation-mediated PCR) assay  for MPS, creating an assay similar to Methyl-Seq . The strength of the HELP assay involves the comparison of the HpaII with the methylation-insensitive MspI representation, allowing a normalization step that makes the assay semi-quantitative . The HELP representation approach was improved upon by Ball et al. , who developed the Methyl-Sensitive Cut Counting (MSCC) assay, which involves digesting DNA with HpaII, ligating an adapter to the cohesive end formed, using a restriction enzyme site within the adapter to digest at a flanking sequence and thus capturing the sequence immediate adjacent to the HpaII site. By adding a second MPS-compatible adapter, a library can be generated for MPS, allowing the counting of reads at these sites to represent the degree of methylation at the site. The authors demonstrated the assay to be reasonably quantitative, testing over 1.3 million sites in the human genome, representing not only HpaII sites clustered in CG-dense regions of the genome (approximately 12% of all HpaII sites are located in annotated CpG islands in the human genome ) but also the remaining majority of the genome in which CG dinucleotides are depleted, a genomic compartment not tested by RRBS as currently designed. A focus on the CG-dense minority of the genome will fail to observe changes such as those at CG-depleted promoters (such as OCT4 ) and CpG island shores , and within gene bodies where cytosine methylation has been found to be positively correlated with gene transcription . It is likely, therefore, that an assay system that can study both CG-dense and CG-depleted regions will acquire substantially more information about epigenomic states than those directed at the CG-dense compartment alone.
In the current study, we tested whether the use of an MspI control would improve MSCC assay performance, as we had found for microarray-based HELP, and whether we could develop an analytical pipeline for routine use of this assay in epigenome-wide association studies. We also explored the use of longer tags than those employed in the MSCC, and added T7 RNA polymerase and reverse transcription steps to allow the generation of libraries without contaminating products, thus obviating the need for gel extraction. The influence of base composition and fragment length parameters as potential sources of bias were also tested, using the H1 (WA01) human embryonic stem (ES) cell line. The outcome is a modified assay that combines the strengths of MSCC and HELP-seq/Methyl-seq, and the supporting analytical workflow that maximizes the quantitative capabilities of the data generated.
Results and discussion
Library preparation and sequencing
Data quality and reproducibility
Based on our experimental design, successfully generated products would be expected to possess a 5'-CGCTGCTG sequence at the 3' end of the read, the first two nucleotides (CG) representing the cohesive end for ligation of HpaII/MspI digestion products, the remaining six nucleotides the EcoP15I restriction enzyme recognition site. In order to evaluate the yield of desired products, we counted the number of reads containing this sequence and plotted the starting positions of this sequence within the reads obtained. We observed that approximately two-thirds of the reads contained the expected sequence, and found that the majority was located at base positions 25 and 26, consistent with the known digestion properties of the restriction enzyme . Removal of the approximately 30% of reads lacking the CG-EcoP15I sequence was performed to eliminate spurious sequences. In order to investigate sequence quality further, we also determined the number and relative position of Ns (ambiguous base calls) within the reads obtained. Overall, few reads were found to contain Ns, and where they were present, they were found to be evenly distributed by position within the sequence. To test data reproducibility, we compared the results of three experimental replicates against each other using the Pearson correlation coefficient metric. The results of this study showed that all replicates were highly correlated (all the r values exceed 0.9), which confirmed that the technical reproducibility of this assay was excellent (Table S2 in Additional file 1).
Distribution of MspI/HpaII sequence tags
We merged three lanes of MspI data and observed that approximately 80% of the 2,292,198 annotated HpaII sites in the human genome (hg18) were represented by at least one read, for a total of over 1.8 million loci throughout the genome. The mean numbers of reads per locus for MspI and HpaII were 3.94 and 1.82, respectively, and MspI counts were distributed evenly across all genomic compartments examined (Table S3 in Additional file 1). We hypothesize that a combination of incomplete genomic coverage and polymorphisms within some CCGG sites (as we have previously observed ) accounts for the 20% of HpaII sites that were not represented by any reads.
Normalization of HpaII by MspI counts and data transformation
Potential sources of bias: base composition and fragment length
Identification of polymorphic CCGG sequences
DNA methylation studies of human embryonic stem cells
The outcome of these studies was an improvement in the previously described MSCC  and HELP-seq  assays, not only by means of technical modifications such as the use of EcoP15I but also because of the concurrent use of MspI for normalization. The effect of these modifications was not only to increase the accuracy of the assay but also to enhance the ability to align sequences to the genome and thus identify polymorphic HpaII/MspI sites. The means of normalization of HpaII by MspI using an angular metric is an innovation that improved the data accuracy substantially and may have applications in other MPS assay normalization strategies. We were also able to discard reads that did not contain the expected adapter sequences, and created a straightforward data analytical pipeline that will facilitate processing of these HELP-tagging data by others.
The potential sources of systematic artifacts due to base composition or digestion product size were evaluated. Apart from a modest decrease in representation in regions above approximately 80% (G+C) content, base composition did not cause biases in representations, possibly in part due to our avoidance of a gel purification step in library preparation . Fragment length does influence the outcome, most likely due to effects on EcoP15I digestion , although the effects should be similar for both HpaII and MspI and should, therefore, largely cancel each other out in the normalization step. It is possible that endogeneous EcoP15I sites could influence the representations, but to have an effect they would have to be located within the 27 bp adjacent to HpaII/MspI sites and would cause digestion of the ligated adapter, causing those loci to be under-represented in both HpaII and MspI datasets. The most likely effect of these endogeneous sites is that they contribute to the proportion of loci at which we could not obtain sequence reads.
Our exploration of the distribution of cytosine methylation in the same human ES cell line studied by Lister et al.  showed consistent results, with hypomethylation of transcription start sites and methylation of transposable elements, as expected from long-standing observations in the field. We furthermore discovered a limited subset of transposable elements that is hypomethylated when in close proximity to transcription start sites. When this subset was studied to determine whether certain types of transposable elements were disproportionately over-represented, we found two broad classes, one of transposable element fossils with no innate capacity to replicate themselves (the ancient DNA, long interspersed repetitive elements (LINEs) and short interspersed repetitive elements (SINEs) shown in Figure 6c) and younger ERV1 long terminal repeat retroelements. Loss of methylation of functionally inactive transposable elements is likely to be of no negative consequence to the host genome, consistent with the host defense hypothesis , while the young ERV long terminal repeats represent a group of transposons whose function has been harnessed as promoters of endogeneous genes [26, 27]. This observation demonstrates the value of a high-resolution, genome-wide assay like HELP-tagging to define potential functional elements in an unbiased manner.
We propose that MPS-based assays such as RRBS , MSCC  and HELP-tagging will prove to be the assays of choice for epigenome-wide association studies in human disease, with the latter two preferable as we begin to explore the CG-depleted majority of the genome. It should not be necessary to run MspI assays every time a HELP-tagging assay is performed, suggesting that a common MspI dataset can serve as a universal reference for a species, allowing a single lane of Illumina sequencing of the HpaII library to provide the methylation data for that sample. The development of analytical pipelines to support analysis of these datasets will be critical to the success of these projects, while the careful ongoing assessment of potential sources of bias will also be essential for improving assay performance.
Materials and methods
Cell preparation and DNA purification
H1 human ES cells (NIH code WA01 from Wicell Research Institute, Madison, WI, USA) were cultured on matrigel (BD Biosciences, San Diego, CA, USA), at 37°C, 5% O2 and 5% CO2. Amplified human ES cell pluripotency was assessed by flow cytometry with SSEA4, CD24 and Oct4 markers. To extract DNA, the cells were suspended in 10 ml of a solution of 10 mM Tris-HCl (pH 8.0), 0.1 M EDTA and 1 ml of 10% SDS to which 10 μl of RNase A (20 mg/ml) was added. After incubation for 1 hour at 37°C, 50 μl of proteinase K (20 mg/ml) was added and the solution was gently mixed and incubated in a 50°C water bath overnight. To purify the lysate, it was extracted three times using saturated phenol, then twice with chloroform, and dialyzed for 16 hours at 4°C against three changes of 0.2× SSC. Following dialysis, the DNA was concentrated by coating the dialysis bags in polyethylene glycol (molecular weight 20,000). The purity and final concentration of the purified DNA was checked by spectrometry (Nanodrop, Wilmington, DE, USA).
Illumina library preparation
The sample preparation steps are illustrated in Figure 1. Two custom adapters were created for HELP-tagging, referred to as AE and AS. As well as an Illumina adapter sequence, adapter AE contains an EcoP15I recognition site and a T7 promoter sequence. Adapter AS contains an Illumina sequencing primer sequence. The adapter and primer sequences for library preparation are listed in Table S9 in Additional file 1. Genomic DNA (5 μg) was digested with HpaII and MspI in separate 200 μl reactions and purified by phenol/chloroform extraction followed by ethanol precipitation. The digested genomic DNA was ligated to adapter AE using a New England Biolabs Quick Ligation Kit (25 μl of 2× Quick ligase buffer, 3 μg of HpaII-digested DNA or 1 μg of MspI-digested DNA, 0.1 μl of Adapter AE (1 μM), 3 μl of Quick Ligase in a final volume of 50 μl). After AE ligation, the products were purified using Agencourt AMpure beads (Beckman Coulter, Brea CA, USA), then digested with EcoP15I (New England Biolabs). The restriction fragments were end-repaired to inhibit to dimerization of adapters, and tailed with a single dA, at the 3' end. After the dA tailing reaction, adapter AS was ligated to the dA-tailed fragments using a New England Biolabs Quick Ligation Kit (25 μl of 2× Quick ligase buffer, 2.5 μl of adapter AS (10 μM), 2.5 μl of Quick Ligase in a final volume 50 μl). After ligation, products were purified using the MinElute PCR purification kit (Qiagen, Hilden, Germany) and in vitro-transcribed using the Ambion MEGAshortscriptkit (Life Technologies, Carlsbad, CA, USA). Following in vitro transcription, products were purified with an RNeasy clean-up kit (Qiagen) before reverse transcription was performed using the Invitrogen SuperScript III kit (Life Technologies). The first strand cDNA produced was used as a template for PCR using the following conditions: 96°C for 2 minutes, then 18 cycles of 96°C for 15 seconds and 72°C for 15 seconds followed by 5 minutes at 72°C for the final extension. After PCR, the library was purified using a QIAQuick PCR clean-up kit (Qiagen).
Single-locus quantitative validation assays
Bisulfite conversion and MassArray (Sequenom, San Diego, CA, USA) were performed using an aliquot of the same sample of DNA as was used for the high-throughput assays described above. Bisulfite conversion was performed with an EZ DNA Methylation kit (Zymo Research, Orange, CA, USA). Bisulfite primers were designed using MethPrimer , specifying the desired product length (250 to 450 bp), primer length (23 to 29 bp) and primer Tm (56 to 62°C). PCR was performed using FastStart High Fidelity Taq polymerase (Roche, Basel, Switzerland) with the following conditions: 95°C for 10 minutes, then 42 cycles of 95°C for 30 seconds, primer-specific Tm for 30 seconds and 72°C for 1 minute, followed by 72°C for 10 minutes for the final extension. Primer-specific Tm and sequence information are provided in Table S6 in Additional file 1. Bisulfite MassArray assays were performed by the institutional Genomics Core Facility. The data were analyzed using the analytical pipeline we have previously described .
Four lanes of sequencing were performed using an Illumina GA IIx Sequencer at the institutional Epigenomics Shared Facility. Three lanes were used for technical replicates of MspI, for the methylation-insensitive reference dataset. Images generated by the Illumina sequencer were analyzed by Illumina pipeline software (versions 1.3 to 1.4). Initial data processing was performed using the default read length of 36 bp, after which we isolated the sequences in which we found adapter sequences on the 3'-end, replaced the adapter sequence with a poly(N) sequence of the same length, and re-ran the Illumina ELAND pipeline again on these sequences with the sequence length set at 27 bp (the 2 to 28 bp subsequence). The data within the ELAND_extended.txt files were used for counting the number of aligned sequences adjacent to each CCGG (HpaII/MspI) site annotated in the hg18 freeze of the human genome at the UCSC genome browser. We permitted up to two mismatches in each sequence, and allowed a sequence to align to up to a maximum of 10 locations within the genome. For non-unique alignments, a sequence was assigned a partial count for each alignment location amounting to 1/n, where n represents the total number of aligned positions. To normalize the data between experiments, the number of sequences associated with each HpaII site was divided by the total number of sequences (including partial counts) aligning to all HpaII sites in the same sample. We refer to this figure as the fixed count below.
To examine an influence of (G+C) mononucleotide content on counts of sequences obtained, we extracted the (G+C) annotation from the hg18 freeze of the human genome at UCSC and examined the distribution of sequence counts according to (G+C) content. Annotated percentages of (G+C) content were available for adjacent 5-bp windows. For each annotated HpaII site, we calculated the mean percentage (G+C) for a 50-bp region centered at the restriction site. Counts of sequences associated with HpaII sites were obtained for 50 sequential non-overlapping windows of 2% (G+C) (the minimum possible in a sample of 50-bp regions). These data were then normalized as a proportion of the total number of fragments. Comparisons were made to the expected frequencies, which for each 2% (G+C) bin was represented by the counts of HpaII sites falling within a range relative to the total number of HpaII sites in the genome. This analysis was performed on both HpaII and MspI-digested DNA for comparison.
The potential effect of distance between HpaII sites on sequences counts obtained at each HpaII site was measured by summing the counts of sequences aligning within each restriction fragment, and normalizing the result with respect to total sequence count. As with the (G+C) analysis above, this was performed for both MspI and HpaII digested restriction fragments. The data were compared with the expected distribution determined by performing virtual restriction digestion using genomic HpaII site coordinates, and normalizing the number of virtual fragments of each size with respect to the total number of these virtual fragments.
guanine and cytosine mononucleotides
HpaII tiny fragment Enrichment by Ligation-mediated PCR
methyl-sensitive cut counting
reduced representation bisulfite sequencing
University of California Santa Cruz.
This work is supported by a grant from the National Institute of Health (NIH, R01 HG004401) to JMG. The authors thank Shahina Maqbool PhD, Raul Olea and Gael Westby of the Einstein Epigenomics Shared Facility for their contributions, Drs Eric Bouhassira and Emmanuel Olivier (Einstein) for the WA01/H1 human ES cell line, and Einstein's Center for Epigenomics.
- Esteller M: Epigenetics in cancer. N Engl J Med. 2008, 358: 1148-1159. 10.1056/NEJMra072067.PubMedView ArticleGoogle Scholar
- Kopelovich L, Crowell JA, Fay JR: The epigenome as a target for cancer chemoprevention. J Natl Cancer Inst. 2003, 95: 1747-1757.PubMedView ArticleGoogle Scholar
- Zilberman D, Gehring M, Tran RK, Ballinger T, Henikoff S: Genome-wide analysis of Arabidopsis thaliana DNA methylation uncovers an interdependence between methylation and transcription. Nat Genet. 2007, 39: 61-69. 10.1038/ng1929.PubMedView ArticleGoogle Scholar
- Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, Pradhan S, Nelson SF, Pellegrini M, Jacobsen SE: Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature. 2008, 452: 215-219. 10.1038/nature06745.PubMedPubMed CentralView ArticleGoogle Scholar
- Lister R, O'Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, Ecker JR: Highly integrated single-base resolution maps of the epigenome in Arabidopsis. Cell. 2008, 133: 523-536. 10.1016/j.cell.2008.03.029.PubMedPubMed CentralView ArticleGoogle Scholar
- Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, Edsall L, Antosiewicz-Bourget J, Stewart R, Ruotti V, Millar AH, Thomson JA, Ren B, Ecker JR: Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009, 462: 315-322. 10.1038/nature08514.PubMedPubMed CentralView ArticleGoogle Scholar
- Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, Susser ES, Slagboom PE, Lumey LH: Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci USA. 2008, 105: 17046-17049. 10.1073/pnas.0806560105.PubMedPubMed CentralView ArticleGoogle Scholar
- Yamashita S, Hosoya K, Gyobu K, Takeshima H, Ushijima T: Development of a novel output value for quantitative assessment in methylated DNA immunoprecipitation-CpG island microarray analysis. DNA Res. 2009, 16: 275-286. 10.1093/dnares/dsp017.PubMedPubMed CentralView ArticleGoogle Scholar
- Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Graf S, Johnson N, Herrero J, Tomazou EM, Thorne NP, Backdahl L, Herberth M, Howe KL, Jackson DK, Miretti MM, Marioni JC, Birney E, Hubbard TJ, Durbin R, Tavare S, Beck S: A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nat Biotechnol. 2008, 26: 779-785. 10.1038/nbt1414.PubMedPubMed CentralView ArticleGoogle Scholar
- Oda M, Glass JL, Thompson RF, Mo Y, Olivier EN, Figueroa ME, Selzer RR, Richmond TA, Zhang X, Dannenberg L, Green RD, Melnick A, Hatchwell E, Bouhassira EE, Verma A, Suzuki M, Greally JM: High-resolution genome-wide cytosine methylation profiling with simultaneous copy number analysis and optimization for limited cell numbers. Nucleic Acids Res. 2009, 37: 3829-3839. 10.1093/nar/gkp260.PubMedPubMed CentralView ArticleGoogle Scholar
- Ordway JM, Bedell JA, Citek RW, Nunberg AN, Jeddeloh JA: MethylMapper: a method for high-throughput, multilocus bisulfite sequence analysis and reporting. Biotechniques. 2005, 39: 464, 466, 468 passim.Google Scholar
- Meissner A, Mikkelsen TS, Gu H, Wernig M, Hanna J, Sivachenko A, Zhang X, Bernstein BE, Nusbaum C, Jaffe DB, Gnirke A, Jaenisch R, Lander ES: Genome-scale DNA methylation maps of pluripotent and differentiated cells. Nature. 2008, 454: 766-770.PubMedPubMed CentralGoogle Scholar
- t Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: e141-10.1093/nar/gkn705.View ArticleGoogle Scholar
- Brunner AL, Johnson DS, Kim SW, Valouev A, Reddy TE, Neff NF, Anton E, Medina C, Nguyen L, Chiao E, Oyolu CB, Schroth GP, Absher DM, Baker JC, Myers RM: Distinct DNA methylation patterns characterize differentiated human embryonic stem cells and developing human fetal liver. Genome Res. 2009, 19: 1044-1056. 10.1101/gr.088773.108.PubMedPubMed CentralView ArticleGoogle Scholar
- Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009, 27: 361-368. 10.1038/nbt.1533.PubMedPubMed CentralView ArticleGoogle Scholar
- Fazzari MJ, Greally JM: Epigenomics: beyond CpG islands. Nat Rev Genet. 2004, 5: 446-455. 10.1038/nrg1349.PubMedView ArticleGoogle Scholar
- Hattori N, Nishino K, Ko YG, Hattori N, Ohgane J, Tanaka S, Shiota K: Epigenetic control of mouse Oct-4 gene expression in embryonic stem cells and trophoblast stem cells. J Biol Chem. 2004, 279: 17063-17069. 10.1074/jbc.M309002200.PubMedView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.PubMedPubMed CentralView ArticleGoogle Scholar
- Moncke-Buchner E, Rothenberg M, Reich S, Wagenfuhr K, Matsumura H, Terauchi R, Kruger DH, Reuter M: Functional characterization and modulation of the DNA cleavage efficiency of type III restriction endonuclease EcoP15I in its interaction with two sites in the DNA target. J Mol Biol. 2009, 387: 1309-1319. 10.1016/j.jmb.2009.02.047.PubMedView ArticleGoogle Scholar
- Greallylab. [http://greallylab.aecom.yu.edu/HELP-tagging/]
- Dohm JC, Lottaz C, Borodina T, Himmelbauer H: Substantial biases in ultra-short read data sets from high-throughput DNA sequencing. Nucleic Acids Res. 2008, 36: e105-10.1093/nar/gkn425.PubMedPubMed CentralView ArticleGoogle Scholar
- Quail MA, Kozarewa I, Smith F, Scally A, Stephens PJ, Durbin R, Swerdlow H, Turner DJ: A large genome center's improvements to the Illumina sequencing system. Nat Methods. 2008, 5: 1005-1010. 10.1038/nmeth.1270.PubMedPubMed CentralView ArticleGoogle Scholar
- Duncan BK, Miller JH: Mutagenic deamination of cytosine residues in DNA. Nature. 1980, 287: 560-561. 10.1038/287560a0.PubMedView ArticleGoogle Scholar
- Thompson RF, Suzuki M, Lau KW, Greally JM: A pipeline for the quantitative analysis of CG dinucleotide methylation using mass spectrometry. Bioinformatics. 2009, 25: 2164-2170. 10.1093/bioinformatics/btp382.PubMedPubMed CentralView ArticleGoogle Scholar
- Yoder JA, Walsh CP, Bestor TH: Cytosine methylation and the ecology of intragenomic parasites. Trends Genet. 1997, 13: 335-340. 10.1016/S0168-9525(97)01181-5.PubMedView ArticleGoogle Scholar
- Conley AB, Piriyapongsa J, Jordan IK: Retroviral promoters in the human genome. Bioinformatics. 2008, 24: 1563-1567. 10.1093/bioinformatics/btn243.PubMedView ArticleGoogle Scholar
- Cohen CJ, Lock WM, Mager DL: Endogenous retroviral LTRs as promoters for human genes: a critical assessment. Gene. 2009, 448: 105-114. 10.1016/j.gene.2009.06.020.PubMedView ArticleGoogle Scholar
- MethPrimer. [http://www.urogene.org/methprimer/]
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited