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.
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. 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.
FDM is highly correlated with JSD* (r = 0.82) when average RNA-seq coverage of the transcripts is sufficiently deep.
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 qRT-PCR of several differential transcription sites identified by FDM, confirmed significant differences at these sites.
Singh D, Orellana CF, Hu Y, Jones CD, Liu Y, Chiang DY, Liu J, Prins JF. (2011) FDM: A Graph-based Statistical Method to Detect Differential Transcription using RNA-seq data. Bioinformatics [Epub ahead of print]. [abstract] [article]