Analyzing whole genome bisulfite sequencing data from highly divergent genotypes

Phillip Wulfridge, Ben Langmead, Andrew P. Feinberg, Kasper D. Hansen

Research output: Contribution to journalArticlepeer-review

3 Scopus citations


In the study of DNA methylation, genetic variation between species, strains or individuals can result in CpG sites that are exclusive to a subset of samples, and insertions and deletions can rearrange the spatial distribution of CpGs. How to account for this variation in an analysis of the interplay between sequence variation and DNA methylation is not well understood, especially when the number of CpG differences between samples is large. Here, we use wholegenome bisulfite sequencing data on two highly divergent mouse strains to study this problem. We show that alignment to personal genomes is necessary for valid methylation quantification. We introduce a method for including strain-specific CpGs in differential analysis, and show that this increases power. We apply our method to a human normalcancer dataset, and show this improves accuracy and power, illustrating the broad applicability of our approach. Our method uses smoothing to impute methylation levels at strain-specific sites, thereby allowing strain-specific CpGs to contribute to the analysis, while accounting for differences in the spatial occurrences of CpGs. Our results have implications for joint analysis of genetic variation and DNA methylation using bisulfite-converted DNA, and unlocks the use of personal genomes for addressing this question.

Original languageEnglish (US)
Pages (from-to)E117
JournalNucleic acids research
Issue number19
StatePublished - Nov 4 2019

ASJC Scopus subject areas

  • Genetics


Dive into the research topics of 'Analyzing whole genome bisulfite sequencing data from highly divergent genotypes'. Together they form a unique fingerprint.

Cite this