Swimming Downstream: Statistical Analysis of Differential Transcript Usage Following Salmon Quantification

Overview

Journal F1000Res

Specialties Biomedical Engineering
Science

Date 2018 Oct 30

PMID 30356428

Citations 74

Authors

Michael I Love

Charlotte Soneson

Rob Patro

Affiliations

Soon will be listed here.

Abstract

Detection of differential transcript usage (DTU) from RNA-seq data is an important bioinformatic analysis that complements differential gene expression analysis. Here we present a simple workflow using a set of existing R/Bioconductor packages for analysis of DTU. We show how these packages can be used downstream of RNA-seq quantification using the Salmon software package. The entire pipeline is fast, benefiting from inference steps by Salmon to quantify expression at the transcript level. The workflow includes live, runnable code chunks for analysis using DRIMSeq and DEXSeq, as well as for performing two-stage testing of DTU using the stageR package, a statistical framework to screen at the gene level and then confirm which transcripts within the significant genes show evidence of DTU. We evaluate these packages and other related packages on a simulated dataset with parameters estimated from real data.

Citing Articles

A systematic benchmark of Nanopore long-read RNA sequencing for transcript-level analysis in human cell lines.

Chen Y, Davidson N, Kei Wan Y, Yao F, Su Y, Gamaarachchi H Nat Methods. 2025; .

PMID: 40082608 DOI: 10.1038/s41592-025-02623-4.

Systematic review and meta-analysis of bulk RNAseq studies in human Alzheimer's disease brain tissue.

Heberle B, Fox K, Libermann L, Xavier S, Dallarosa G, Santos R Alzheimers Dement. 2025; 21(3):e70025.

PMID: 40042520 PMC: 11881636. DOI: 10.1002/alz.70025.

Temporal phosphoproteomics reveals circuitry of phased propagation in insulin signaling.

Turewicz M, Skagen C, Hartwig S, Majda S, Thedinga K, Herwig R Nat Commun. 2025; 16(1):1570.

PMID: 39939313 PMC: 11821911. DOI: 10.1038/s41467-025-56335-6.

A fast hybridization chain reaction method in embryos and ovaries.

Mikami K, Kozono Y, Masukawa M, Kobayashi S Fly (Austin). 2024; 19(1):2428499.

PMID: 39639000 PMC: 11633216. DOI: 10.1080/19336934.2024.2428499.

Systematic review and meta-analysis of bulk RNAseq studies in human Alzheimer's disease brain tissue.

Heberle B, Fox K, Libermann L, Xavier S, Dallarosa G, Santos R bioRxiv. 2024; .

PMID: 39574617 PMC: 11580990. DOI: 10.1101/2024.11.07.622520.

References

Heller R, Manduchi E, Grant G, Ewens W . A flexible two-stage procedure for identifying gene sets that are differentially expressed. Bioinformatics. 2009; 25(8):1019-25. DOI: 10.1093/bioinformatics/btp076. View

McCarthy D, Chen Y, Smyth G . Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012; 40(10):4288-97. PMC: 3378882. DOI: 10.1093/nar/gks042. View

Reyes A, Huber W . Alternative start and termination sites of transcription drive most transcript isoform differences across human tissues. Nucleic Acids Res. 2017; 46(2):582-592. PMC: 5778607. DOI: 10.1093/nar/gkx1165. View

Van den Berge K, Soneson C, Robinson M, Clement L . stageR: a general stage-wise method for controlling the gene-level false discovery rate in differential expression and differential transcript usage. Genome Biol. 2017; 18(1):151. PMC: 5547545. DOI: 10.1186/s13059-017-1277-0. View

Battle A, Brown C, Engelhardt B, Montgomery S . Genetic effects on gene expression across human tissues. Nature. 2017; 550(7675):204-213. PMC: 5776756. DOI: 10.1038/nature24277. View

Soneson C, Love M, Robinson M . Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. 2016; 4:1521. PMC: 4712774. DOI: 10.12688/f1000research.7563.2. View

Nowicka M, Robinson M . DRIMSeq: a Dirichlet-multinomial framework for multivariate count outcomes in genomics. F1000Res. 2017; 5:1356. PMC: 5200948. DOI: 10.12688/f1000research.8900.2. View

Collado-Torres L, Nellore A, Kammers K, Ellis S, Taub M, Hansen K . Reproducible RNA-seq analysis using recount2. Nat Biotechnol. 2017; 35(4):319-321. PMC: 6742427. DOI: 10.1038/nbt.3838. View

Scotti M, Swanson M . RNA mis-splicing in disease. Nat Rev Genet. 2015; 17(1):19-32. PMC: 5993438. DOI: 10.1038/nrg.2015.3. View

10.

Trapnell C, Hendrickson D, Sauvageau M, Goff L, Rinn J, Pachter L . Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2012; 31(1):46-53. PMC: 3869392. DOI: 10.1038/nbt.2450. View

11.

Pimentel H, Bray N, Puente S, Melsted P, Pachter L . Differential analysis of RNA-seq incorporating quantification uncertainty. Nat Methods. 2017; 14(7):687-690. DOI: 10.1038/nmeth.4324. View

12.

Love M, Soneson C, Patro R . Swimming downstream: statistical analysis of differential transcript usage following Salmon quantification. F1000Res. 2018; 7:952. PMC: 6178912. DOI: 10.12688/f1000research.15398.3. View

13.

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

14.

Bray N, Pimentel H, Melsted P, Pachter L . Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016; 34(5):525-7. DOI: 10.1038/nbt.3519. View

15.

Blume J, DAgostino McGowan L, Dupont W, Greevy Jr R . Second-generation p-values: Improved rigor, reproducibility, & transparency in statistical analyses. PLoS One. 2018; 13(3):e0188299. PMC: 5863943. DOI: 10.1371/journal.pone.0188299. View

16.

Soneson C, Matthes K, Nowicka M, Law C, Robinson M . Isoform prefiltering improves performance of count-based methods for analysis of differential transcript usage. Genome Biol. 2016; 17:12. PMC: 4729156. DOI: 10.1186/s13059-015-0862-3. View

17.

Soneson C, Robinson M . iCOBRA: open, reproducible, standardized and live method benchmarking. Nat Methods. 2016; 13(4):283. DOI: 10.1038/nmeth.3805. View

18.

Soneson C, Robinson M . Towards unified quality verification of synthetic count data with countsimQC. Bioinformatics. 2017; 34(4):691-692. PMC: 5860609. DOI: 10.1093/bioinformatics/btx631. View

19.

Vitting-Seerup K, Sandelin A . IsoformSwitchAnalyzeR: analysis of changes in genome-wide patterns of alternative splicing and its functional consequences. Bioinformatics. 2019; 35(21):4469-4471. DOI: 10.1093/bioinformatics/btz247. View

20.

Koster J, Rahmann S . Snakemake--a scalable bioinformatics workflow engine. Bioinformatics. 2012; 28(19):2520-2. DOI: 10.1093/bioinformatics/bts480. View