Bioinformatics 2014 Jan 18;30(2):165-71. Epub 2013 Nov 18.
Institute of Bioengineering, School of Life Sciences, École Polytechnique Fédérale de Lausanne, Swiss Institute of Bioinformatics, 1015, Lausanne, Switzerland, Department of Genetic Medicine and Development, University of Geneva Medical School, Institute of Genetics and Genomics in Geneva, University of Geneva, 1211, Geneva, Switzerland and Center for Integrative Genomics, Faculty of Biology and Medicine, University of Lausanne, 1011, Lausanne, Switzerland.
Motivation: High-throughput sequencing technologies enable the genome-wide analysis of the impact of genetic variation on molecular phenotypes at unprecedented resolution. However, although powerful, these technologies can also introduce unexpected artifacts.
Results: We investigated the impact of library amplification bias on the identification of allele-specific (AS) molecular events from high-throughput sequencing data derived from chromatin immunoprecipitation assays (ChIP-seq). Putative AS DNA binding activity for RNA polymerase II was determined using ChIP-seq data derived from lymphoblastoid cell lines of two parent-daughter trios. We found that, at high-sequencing depth, many significant AS binding sites suffered from an amplification bias, as evidenced by a larger number of clonal reads representing one of the two alleles. To alleviate this bias, we devised an amplification bias detection strategy, which filters out sites with low read complexity and sites featuring a significant excess of clonal reads. This method will be useful for AS analyses involving ChIP-seq and other functional sequencing assays.
Availability: The R package abs filter for library clonality simulations and detection of amplification-biased sites is available from http://updepla1srv1.epfl.ch/waszaks/absfilter