Rxivist logo

A junction coverage compatibility score to quantify the reliability of transcript abundance estimates and annotation catalogs

By Charlotte Soneson, Michael I Love, Rob Patro, Shobbir Hussain, Dheeraj Malhotra, Mark D Robinson

Posted 28 Jul 2018
bioRxiv DOI: 10.1101/378539 (published DOI: 10.26508/lsa.201800175)

Most methods for statistical analysis of RNA-seq data take a matrix of abundance estimates for some type of genomic features as their input, and consequently the quality of any obtained results are directly dependent on the quality of these abundances. Here, we present the junction coverage compatibility (JCC) score, which provides a way to evaluate the reliability of transcript-level abundance estimates as well as the accuracy of transcript annotation catalogs. It works by comparing the observed number of reads spanning each annotated splice junction in a genomic region to the predicted number of junction-spanning reads, inferred from the estimated transcript abundances and the genomic coordinates of the corresponding annotated transcripts. We show that while most genes show good agreement between the observed and predicted junction coverages, there is a small set of genes that do not. Genes with poor agreement are found regardless of the method used to estimate transcript abundances, and the corresponding transcript abundances should be treated with care in any downstream analyses.

Download data

  • Downloaded 635 times
  • Download rankings, all-time:
    • Site-wide: 50,093
    • In bioinformatics: 5,106
  • Year to date:
    • Site-wide: 120,012
  • Since beginning of last month:
    • Site-wide: 116,187

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide