Enhanced Methylation Analysis by Recovery of Unsequenceable Fragments

Bisulfite sequencing is a valuable tool for mapping the position of 5-methylcytosine in the genome at single base resolution. However, the associated chemical treatment causes strand scission, which depletes the number of sequenceable DNA fragments in a library and thus necessitates PCR amplification. The AT-rich nature of the library generated from bisulfite treatment adversely affects this amplification, resulting in the introduction of major biases that can confound methylation analysis. Here, we report a method that enables more accurate methylation analysis, by rebuilding bisulfite-damaged components of a DNA library. This recovery after bisulfite treatment (ReBuilT) approach enables PCR-free bisulfite sequencing from low nanogram quantities of genomic DNA. We apply the ReBuilT method for the first whole methylome analysis of the highly AT-rich genome of Plasmodium berghei. Side-by-side comparison to a commercial protocol involving amplification demonstrates a substantial improvement in uniformity of coverage and reduction of sequence context bias. Our method will be widely applicable for quantitative methylation analysis, even for technically challenging genomes, and where limited sample DNA is available.

[1]  S. Henikoff,et al.  DNA Methylation Profiling Identifies CG Methylation Clusters in Arabidopsis Genes , 2005, Current Biology.

[2]  F. Miura,et al.  Amplification-free whole-genome bisulfite sequencing by post-bisulfite adaptor tagging , 2012, Nucleic acids research.

[3]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[4]  Zuhong Lu,et al.  A novel method to quantify local CpG methylation density by regional methylation elongation assay on microarray , 2008, BMC Genomics.

[5]  Lee E. Edsall,et al.  Human DNA methylomes at base resolution show widespread epigenomic differences , 2009, Nature.

[6]  Robert J. Schmitz,et al.  Methylated DNA is over-represented in whole-genome bisulfite sequencing data , 2014, Front. Genet..

[7]  M. T. McCabe,et al.  Cancer DNAMethylation : Molecular Mechanisms and Clinical Implications , 2009 .

[8]  M. T. McCabe,et al.  Cancer DNA Methylation: Molecular Mechanisms and Clinical Implications , 2009, Clinical Cancer Research.

[9]  Rudolf Jaenisch,et al.  Role for DNA methylation in genomic imprinting , 1993, Nature.

[10]  N. Heintz,et al.  The Nuclear DNA Base 5-Hydroxymethylcytosine Is Present in Purkinje Neurons and the Brain , 2009, Science.

[11]  M. Duraisingh,et al.  Epigenetics in Plasmodium: What Do We Really Know? , 2010, Eukaryotic Cell.

[12]  R. Sager,et al.  Tissue specificity and clustering of methylated cystosines in bovine satellite I DNA. , 1982, Proceedings of the National Academy of Sciences of the United States of America.

[13]  G. Ast,et al.  DNA-methylation effect on cotranscriptional splicing is dependent on GC architecture of the exon–intron structure , 2013, Genome research.

[14]  David R. Liu,et al.  Conversion of 5-Methylcytosine to 5- Hydroxymethylcytosine in Mammalian DNA by the MLL Partner TET1 , 2009 .

[15]  Roland Eils,et al.  Tagmentation-based whole-genome bisulfite sequencing , 2013, Nature Protocols.

[16]  David R. Liu,et al.  The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing , 2010, PloS one.

[17]  Kees-Jan Françoijs,et al.  Linear amplification for deep sequencing , 2011, Nature Protocols.

[18]  Fidel Ramírez,et al.  deepTools: a flexible platform for exploring deep-sequencing data , 2014, Nucleic Acids Res..

[19]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[20]  L. E. McDonald,et al.  A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[21]  M. Ehrlich,et al.  Comparison of bisulfite modification of 5-methyldeoxycytidine and deoxycytidine residues. , 1980, Nucleic acids research.

[22]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[23]  Zbynek Bozdech,et al.  Dynamic Epigenetic Regulation of Gene Expression during the Life Cycle of Malaria Parasite Plasmodium falciparum , 2013, PLoS pathogens.

[24]  Akimitsu Okamoto,et al.  Degradation of DNA by bisulfite treatment. , 2007, Bioorganic & medicinal chemistry letters.

[25]  Bing Ren,et al.  Tet-assisted bisulfite sequencing of 5-hydroxymethylcytosine , 2012, Nature Protocols.

[26]  S. Lonardi,et al.  Genome-wide mapping of DNA methylation in the human malaria parasite Plasmodium falciparum. , 2013, Cell host & microbe.

[27]  Lee T. Sam,et al.  Deep sequencing reveals distinct patterns of DNA methylation in prostate cancer. , 2011, Genome research.

[28]  Marcel Martin.  Cutadapt removes adapter sequences from high-throughput sequencing reads , 2011 .

[29]  N. A. Malmquist,et al.  Small-molecule histone methyltransferase inhibitors display rapid antimalarial activity against all blood stage forms in Plasmodium falciparum , 2012, Proceedings of the National Academy of Sciences.

[30]  T. Fennell,et al.  Analyzing and minimizing PCR amplification bias in Illumina sequencing libraries , 2011, Genome Biology.

[31]  O. Stegle,et al.  Single-Cell Genome-Wide Bisulfite Sequencing for Assessing Epigenetic Heterogeneity , 2014, Nature Methods.

[32]  S. Balasubramanian,et al.  Quantitative Sequencing of 5-Methylcytosine and 5-Hydroxymethylcytosine at Single-Base Resolution , 2012, Science.

[33]  Brent S. Pedersen,et al.  Fast and accurate alignment of long bisulfite-seq reads , 2014, 1401.1129.

[34]  Aaron R. Quinlan,et al.  Bioinformatics Applications Note Genome Analysis Bedtools: a Flexible Suite of Utilities for Comparing Genomic Features , 2022 .