Abstract:MOTIVATION: In eukaryotic cells, alternative splicing expands the diversity of RNA transcripts and plays an important role in tissue-specific differentiation, and can be misregulated in disease. To understand these processes, there is a great need for methods to detect differential transcription between samples. Our focus is on samples observed using short-read RNA sequencing (RNA-seq).METHODS: We characterize differential transcription between two samples as the difference in the relative abundance of the transcript isoforms present in the samples. The magnitude of differential transcription of a gene between two samples can be measured by the square root of the Jensen Shannon Divergence (JSD*) between the gene's transcript abundance vectors in each sample. We define a weighted splice-graph representation of RNA-seq data, summarizing in compact form the alignment of RNA-seq reads to a reference genome. The flow difference metric (FDM) identifies regions of differential RNA transcript expression between pairs of splice graphs, without need for an underlying gene model or catalog of transcripts. We present a novel non-parametric statistical test between splice graphs to assess the significance of differential transcription, and extend it to group-wise comparison incorporating sample replicates.RESULTS: Using simulated RNA-seq data consisting of four technical replicates of two samples with varying transcription between genes, we show that (i) the FDM is highly correlated with JSD* (r=0.82) when average RNA-seq coverage of the transcripts is sufficiently deep; and (ii) the FDM is able to identify 90% of genes with differential transcription when JSD* >0.28 and coverage >7. This represents higher sensitivity than Cufflinks (without annotations) and rDiff (MMD), which respectively identified 69 and 49% of the genes in this region as differential transcribed. Using annotations identifying the transcripts, Cufflinks was able to identify 86% of the genes in this region as differentially transcribed. Using experimental data consisting of four replicates each for two cancer cell lines (MCF7 and SUM102), FDM identified 1425 genes as significantly different in transcription. Subsequent study of the samples using quantitative real time polymerase chain reaction (qRT-PCR) of several differential transcription sites identified by FDM, confirmed significant differences at these sites.AVAILABILITY: http://csbio-linux001.cs.unc.edu/nextgen/software/FDM CONTACT: darshan@email.unc.eduSUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

FDM: a graph-based statistical method to detect differential transcription using RNA-seq data

PennDiff: Detecting Differential Alternative Splicing and Transcription by RNA Sequencing.

Degps is a Powerful Tool for Detecting Differential Expression in RNA-sequencing Studies

Rseqdiff: Detecting Differential Isoform Expression from RNA-Seq Data Using Hierarchical Likelihood Ratio Test.

RNA-Seq Analysis of Differential Splice Junction Usage and Intron Retentions by DEXSeq

Identifying Differentially Spliced Genes from Two Groups of Rna-Seq Samples

Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation

A comprehensive benchmarking of differential splicing tools for RNA-seq analysis at the event level

Comprehensive benchmark of differential transcript usage analysis for static and dynamic conditions

Differential analysis in Transcriptomic: The strength of randomly picking 'reference' genes

IFDlong: an isoform and fusion detector for accurate annotation and quantification of long-read RNA-seq data

Detecting differential usage of exons from RNA-seq data

Detection and Visualization of Differential Splicing in RNA-Seq Data with JunctionSeq

DifferentialRegulation: a Bayesian hierarchical approach to identify differentially regulated genes

DRME: Count-based Differential RNA Methylation Analysis at Small Sample Size Scenario

Differential analysis of RNA methylome with improved spatial resolution

Gene Ranking of RNA-Seq Data Via Discriminant Non-Negative Matrix Factorization

Selecting Differential Splicing Methods: Practical Considerations

: intuitive analysis and visualisation of differential alternative splicing using beta distributions

RADPRE: A Computational Program for Identification of Differential Mrna Processing Including Alternative Polyadenylation

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