Evidence map›Paper›PMID 29444201›Full record

ArticleBioinformatics (Oxford, England)2018

Hierarchical analysis of RNA-seq reads improves the accuracy of allele-specific expression.

Narayanan Raghupathy, Kwangbom Choi, Matthew J Vincent, Glen L Beane, Keith S Sheppard, Steven C Munger, Ron Korstanje, Fernando Pardo-Manual de Villena, Gary A Churchill

Abstract read
In one paragraph

Article in Bioinformatics (Oxford, England), 2018. The graph could read no effect estimate from its abstract, so it casts no vote on the map. Cited by 67 papers.

0numbers the graph read from it
0cells of the map it votes in
67citing papers in PubMed
–field-weighted citation impact
1 · What the graph read from it

What it found

Each row is one number read from the abstract, on the scale the paper reported it, with its interval. Left of the dashed line favours the treatment, right favours the comparator. Under each row is the sentence it came from. New to these charts? A ten-minute tutorial.

The abstract states no effect estimate the extractor could read, or names no intervention and outcome on the map, so this paper lights no cell and moves no belief. It is still indexed, cited and linked below.

2 · The registry

The trial behind it

Trials whose registry record cites this paper, or whose number appears in the abstract. A trial that started after this paper was published is citing it as background, not reporting it.

Neither the registry nor the abstract names a trial number. If this is a trial report, that itself is worth knowing.

3 · Its place in the literature

Who cites it

67 citing papers in PubMed.

  1. Article
  2. Article
  3. Article
  4. Article
  5. Article
  6. Article
  7. Article
  8. Article
  9. Low-coverage whole-genome sequencing facilitates accurate and cost-effective haplotype reconstruction in complex mouse crosses.Mammalian genome : official journal of the International Mammalian Genome Society · 2025
    Article
  10. Article
  11. Article
  12. Article
  13. Article
  14. Article
  15. Article
  16. Article
  17. The regulatory architecture of gene expression variation inbioRxiv : the preprint server for biology · 2025
    Article
  18. Article
  19. Article
  20. Article

7 more citing papers are in PubMed but not listed here.

4 · The record

Corrections and comments

PubMed lists nothing against this paper. Absence here is not a guarantee, only a check that was made.

5 · Who and what money

Authors and funding

9 authors.

Narayanan RaghupathyThe Jackson Laboratory, Bar Harbor, USA.
Kwangbom ChoiThe Jackson Laboratory, Bar Harbor, USA.
Matthew J VincentThe Jackson Laboratory, Bar Harbor, USA.
Glen L BeaneThe Jackson Laboratory, Bar Harbor, USA.
Keith S SheppardThe Jackson Laboratory, Bar Harbor, USA.
Steven C MungerThe Jackson Laboratory, Bar Harbor, USA.
Ron KorstanjeThe Jackson Laboratory, Bar Harbor, USA.
Fernando Pardo-Manual de VillenaDepartment of Genetics, The University of North Carolina, Chapel Hill, USA.
Gary A ChurchillThe Jackson Laboratory, Bar Harbor, USA.

Funding

Shared Resource ManagementP30CA034196 · NCI · JACKSON LABORATORY · PI William Carl Skarnes · 1985 to 2026
$61.9M
UNC and PERLEGEN, IncP50GM076468 · NIGMS · JACKSON LABORATORY · PI PETKOV, PETKO M · 2006 to 2015
$31.0M
NCI NIH HHS P30 CA034196NIGMS NIH HHS P50 GM076468
6 · The paper itself

Abstract

Motivation: Allele-specific expression (ASE) refers to the differential abundance of the allelic copies of a transcript. RNA sequencing (RNA-seq) can provide quantitative estimates of ASE for genes with transcribed polymorphisms. When short-read sequences are aligned to a diploid transcriptome, read-mapping ambiguities confound our ability to directly count reads. Multi-mapping reads aligning equally well to multiple genomic locations, isoforms or alleles can comprise the majority (>85%) of reads. Discarding them can result in biases and substantial loss of information. Methods have been developed that use weighted allocation of read counts but these methods treat the different types of multi-reads equivalently. We propose a hierarchical approach to allocation of read counts that first resolves ambiguities among genes, then among isoforms, and lastly between alleles. We have implemented our model in EMASE software (Expectation-Maximization for Allele Specific Expression) to estimate total gene expression, isoform usage and ASE based on this hierarchical allocation. Results: Methods that align RNA-seq reads to a diploid transcriptome incorporating known genetic variants improve estimates of ASE and total gene expression compared to methods that use reference genome alignments. Weighted allocation methods outperform methods that discard multi-reads. Hierarchical allocation of reads improves estimation of ASE even when data are simulated from a non-hierarchical model. Analysis of RNA-seq data from F1 hybrid mice using EMASE reveals widespread ASE associated with cis-acting polymorphisms and a small number of parent-of-origin effects. Availability and implementation: EMASE software is available at https://github.com/churchill-lab/emase. Supplementary information: Supplementary data are available at Bioinformatics online.

Indexed as

AllelesAlternative SplicingSoftwareTranscriptomeAnimalsGenomicsMaleMiceSequence Analysis, RNA

Identifiers

PMID29444201
PMCPMC6022640

What OpenQuestion holds

Textmetadata
Read underepoch 390

Registered trials

None linked

Read under generation 80e0d062 · epoch 390. Bibliography from PubMed, PubMed Central and OpenAlex; grants from NIH RePORTER; trial links from ClinicalTrials.gov; estimates, votes and beliefs from the OpenQuestion graph.