Multiplex padlock capture and sequencing reveal human hypermutable CpG variations

Abstract

Utilizing the full power of next generation sequencing often requires the ability to perform large scale multiplex enrichment of many specific genomic loci in multiple samples. Several technologies have been recently developed but await substantial improvements. We report the 10,000-fold improvement of a previously developed padlock-based approach, and extend the assay to identifying genetic variations in hypermutable CpG regions across human chromosome 21. From ~3 million reads derived from a single Illumina Genome Analyzer lane, ~94% (~50,500) target sites can be observed with at least one read. The uniformity of coverage was also greatly improved; up to 93% and 57% of all targets fell within a 100- and 10-fold coverage range, respectively. Alleles at >400,000 target base positions were determined across six subjects and examined for single nucleotide polymorphisms, and the concordance with independently obtained genotypes was 98.4-100%. We detected >500 SNPs not currently in dbSNP; 362 of which were in targeted CpG locations. Transitions in CpG sites were at least 13.7 times more abundant than non-CpG transitions. Fractions of polymorphic CpG sites are lower in CpG-rich regions and show higher correlation with human-chimpanzee divergence within CpG vs. non-CpG sites. This is consistent with the hypothesis that methylation rate heterogeneity along chromosomes contributes to mutation rate variation in humans. Our success suggests that targeted CpG resequencing is an efficient way for identifying common and rare genetic variations. In addition, the significantly improved padlock capture technology can be readily applied to other projects that require multiplex sample preparation.