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.

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

LongGF: Computational Algorithm and Software Tool for Fast and Accurate Detection of Gene Fusions by Long-Read Transcriptome Sequencing

A Novel Analytical Strategy To Identify Fusion Transcripts Between Repetitive Elements And Protein Coding-Exons Using Rna-Seq

CTAT-LR-fusion: accurate fusion transcript identification from long and short read isoform sequencing at bulk or single cell resolution

Comprehensive Assessment of Isoform Detection Methods for Third-Generation Sequencing Data

LncFinder: an Integrated Platform for Long Non-Coding RNA Identification Utilizing Sequence Intrinsic Composition, Structural Information and Physicochemical Property.

Comprehensive assessment of mRNA isoform detection methods for long-read sequencing data

Comprehensive Evaluation of Fusion Transcript Detection Algorithms and a Meta-Caller to Combine Top Performing Methods in Paired-End RNA-seq Data.

Enhancing novel isoform discovery: leveraging nanopore long-read sequencing and machine learning approaches

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

IsoDetect: Detection of Splice Isoforms from Third Generation Long Reads Based on Short Feature Sequences

LIQA: long-read isoform quantification and analysis

DEEP-LONG: A Fast and Accurate Aligner for Long RNA-Seq

FusionQ: a Novel Approach for Gene Fusion Detection and Quantification from Paired-End RNA-Seq.

An optimized workflow of full-length transcriptome sequencing for accurate fusion transcript identification

Flnc: Machine Learning Improves the Identification of Novel Long Noncoding RNAs from Stand-Alone RNA-Seq Data.

Abstract 7418: LongFuse: Detecting gene fusion transcripts from high throughput long-read single cell RNA sequencing data

Flnc: Machine Learning Improves the Identification of Novel Full-length Long Noncoding RNAs from RNA Sequencing Data Without Transcriptional Initiation Profiles

Lncrna-Mfdl: Identification of Human Long Non-Coding Rnas by Fusing Multiple Features and Using Deep Learning

LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions.

Comprehensive characterization of single-cell full-length isoforms in human and mouse with long-read sequencing