Fine-mapping cellular QTLs with RASQUAL and ATAC-seq
Open Access
- 14 December 2015
- journal article
- research article
- Published by Springer Science and Business Media LLC in Nature Genetics
- Vol. 48 (2), 206-213
- https://doi.org/10.1038/ng.3467
Abstract
Natsuhiko Kumasaka, Andrew Knights and Daniel Gaffney develop a new statistical approach for association mapping that models genetic effects and accounts for biases in sequencing data in a single probabilistic framework. They apply this method to generate a map of chromatin accessibility QTLs and show how it can be used to fine-map regulatory variants and link distal regulatory elements with genes. When cellular traits are measured using high-throughput DNA sequencing, quantitative trait loci (QTLs) manifest as fragment count differences between individuals and allelic differences within individuals. We present RASQUAL (Robust Allele-Specific Quantitation and Quality Control), a new statistical approach for association mapping that models genetic effects and accounts for biases in sequencing data using a single, probabilistic framework. RASQUAL substantially improves fine-mapping accuracy and sensitivity relative to existing methods in RNA-seq, DNase-seq and ChIP-seq data. We illustrate how RASQUAL can be used to maximize association detection by generating the first map of chromatin accessibility QTLs (caQTLs) in a European population using ATAC-seq. Despite a modest sample size, we identified 2,707 independent caQTLs (at a false discovery rate of 10%) and demonstrated how RASQUAL and ATAC-seq can provide powerful information for fine-mapping gene-regulatory variants and for linking distal regulatory elements with gene promoters. Our results highlight how combining between-individual and allele-specific genetic signals improves the functional interpretation of noncoding variation.Keywords
This publication has 43 references indexed in Scilit:
- TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusionsGenome Biology, 2013
- Matrix eQTL: ultra fast eQTL analysis via large matrix operationsBioinformatics, 2012
- Fast gapped-read alignment with Bowtie 2Nature Methods, 2012
- DNase I sensitivity QTLs are a major determinant of human expression variationNature, 2012
- False positive peaks in ChIP-seq and other sequencing-based functional assays caused by unannotated high copy number regionsBioinformatics, 2011
- Transcriptome genetics using second generation sequencing in a Caucasian populationNature, 2010
- Understanding mechanisms underlying human gene expression variation with RNA sequencingNature, 2010
- Genome-wide analysis of allelic expression imbalance in human primary cells by high-throughput transcriptome resequencingHuman Molecular Genetics, 2009
- Effect of read-mapping biases on detecting allele-specific expression from RNA-sequencing dataBioinformatics, 2009
- Fast and accurate short read alignment with Burrows–Wheeler transformBioinformatics, 2009