Rxivist logo

QuASAR-MPRA: Accurate allele-specific analysis for massively parallel reporter assays

By Cynthia A Kalita, Gregory A. Moyerbrailean, Christopher Brown, Xiaoquan Wen, Francesca Luca, Roger Pique-Regi

Posted 03 Feb 2017
bioRxiv DOI: 10.1101/105627 (published DOI: 10.1093/bioinformatics/btx598)

The majority of the human genome is composed of non-coding regions containing regulatory elements such as enhancers, which are crucial for controlling gene expression. Many variants associated with complex traits are in these regions, and may disrupt gene regulatory sequences. Consequently, it is important to not only identify true enhancers but also to test if a variant within an enhancer affects gene regulation. Recently, allele-specific analysis in high-throughput reporter assays, such as massively parallel reporter assays (MPRA), have been used to functionally validate non-coding variants. However, we are still missing high-quality and robust data analysis tools for these datasets. We have further developed our method for allele-specific analysis QuASAR (quantitative allele-specific analysis of reads) to analyze allele-specific signals in barcoded read counts data from MPRAs. Using this approach, we can take into account the uncertainty on the original plasmid proportions, over-dispersion, and sequencing errors. The provided allelic skew estimate and its standard error also simplifies meta-analysis of replicate experiments. Additionally, we show that a beta-binomial distribution better models the variability present in the allelic imbalance of these synthetic reporters and results in a test that is statistically well calibrated under the null. Applying this approach to the MPRA data by Tewhey et al (2016), we find 602 SNPs with significant (FDR 10%) allele-specific regulatory function in LCLs. We also show that we can combine MPRA with QuASAR estimates to validate existing experimental and computational annotations of regulatory variants. Our study shows that by having the appropriate data analysis tools, we can greatly improve the power to detect allelic effects in high throughput reporter assays.

Download data

  • Downloaded 1,038 times
  • Download rankings, all-time:
    • Site-wide: 19,901
    • In bioinformatics: 2,366
  • Year to date:
    • Site-wide: 115,954
  • Since beginning of last month:
    • Site-wide: 98,461

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide


Sign up for the Rxivist weekly newsletter! (Click here for more details.)