Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences

Abstract

High-throughput sequencing of cDNA (RNA-seq) is used extensively to characterize the transcriptome of cells. Many transcriptomic studies aim at comparing either abundance levels or the transcriptome composition between given conditions, and as a first step, the sequencing reads must be used as the basis for abundance quantification of transcriptomic features of interest, such as genes or transcripts. Various quantification approaches have been proposed, ranging from simple counting of reads that overlap given genomic regions to more complex estimation of underlying transcript abundances. In this paper, we show that gene-level abundance estimates and statistical inference offer advantages over transcript-level analyses, in terms of performance and interpretability. We also illustrate that the presence of differential isoform usage can lead to inflated false discovery rates in differential gene expression analyses on simple count matrices but that this can be addressed by incorporating offsets derived from transcript-level abundance estimates. We also show that the problem is relatively minor in several real data sets. Finally, we provide an R package ( tximport ) to help users integrate transcript-level abundance estimates from common quantification pipelines into count-based statistical inference engines.

Cited by papers we hold 62

• Alternatively spliced isoforms of IRF7 differentially regulate interferon expression to tune response to viral infection · 2025
  Not checked
• edgeR v4: powerful differential analysis of sequencing data with expanded functionality and improved support for small counts and larger datasets · 2025
  No notice
• Aging clocks based on accumulating stochastic variation · 2024
  No notice
• Tumour-selective activity of RAS-GTP inhibition in pancreatic cancer · 2024
  Corrected
• Safety, efficacy and determinants of response of allogeneic CD19-specific CAR-NK cells in CD19+ B cell tumors: a phase 1/2 trial · 2024
  No notice
• Vaccine-boosted CAR T crosstalk with host immunity to reject tumors with antigen heterogeneity · 2023
  No notice
• Platelet-derived exerkine CXCL4/platelet factor 4 rejuvenates hippocampal neurogenesis and restores cognitive function in aged mice · 2023
  No notice
• Systematic transcriptional analysis of human cell lines for gene expression landscape and tumor representation · 2023
  No notice
• A CRISPRi/a platform in human iPSC-derived microglia uncovers regulators of disease states · 2022
  No notice
• Single-cell RNA sequencing reveals evolution of immune landscape during glioblastoma progression · 2022
  No notice
• The complete sequence of a human genome · 2022
  No notice
• Strain-specific impacts of probiotics are a significant driver of gut microbiome development in very preterm infants · 2022
  No notice
• Acquired mutations and transcriptional remodeling in long-term estrogen-deprived locoregional breast cancer recurrences · 2021
  Not checked
• Serine residues 726 and 780 have nonredundant roles regulating STAT5a activity in luminal breast cancer · 2021
  Not checked
• Differential gene expression in Drosophila melanogaster and D. nigrosparsa infected with the same Wolbachia strain · 2021
  Not checked
• Autosomal recessive loci contribute significantly to quantitative variation of male fertility in a dairy cattle population · 2021
  Not checked
• RNA-Binding Protein Expression Alters Upon Differentiation of Human B Cells and T Cells · 2021
  Not checked
• Coordination of microbe–host homeostasis by crosstalk with plant innate immunity · 2021
  No notice
• TPM, FPKM, or Normalized Counts? A Comparative Study of Quantification Measures for the Analysis of RNA-seq Data from the NCI Patient-Derived Models Repository · 2021
  No notice
• Altered lipid metabolism marks glioblastoma stem and non-stem cells in separate tumor niches · 2021
  No notice

and 42 more

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.