Question: Differential quantification of splice junction
gravatar for Anil Kesarwani
3.9 years ago by
United States
Anil Kesarwani80 wrote:


From TopHat 2 alignment, I have obtained outputs containg splice junctions and the corresponding read counts. I am interested in calculating differential splice junction usage between wild type and mutant samples.

In a recent publication ( ), the DEXSeq function called "testForDEURT" was used to quantify differential splice junction usage between two samples. However, I am not able to figure out the important steps to use this function.

I have wild type and mutant samples in replicates, and therefire some statistics, for example, p-value and FDR, will be needed to define signficantly candidates.

Could somebody please guide me doing this analysis. Any response is highly appreciated.


rna-seq dexseq • 1.3k views
ADD COMMENTlink modified 3.9 years ago by SES8.1k • written 3.9 years ago by Anil Kesarwani80

The only actual important steps that are different from the normal DEXSeq workflow are (1) filtering the splice site counts (they mention the filter used in the methods) and (2) annotating splice sites with gene information to allow parsing by DEXseq. You'll have to write something to do (2). (1) can likely be done with awk.

ADD REPLYlink written 3.9 years ago by Devon Ryan88k

To me, the implementaiton of DEXSeq to calculate differential splice junction usage is not straightforward. Could you please explain me in more detail. I managed to prepare my read count (for splice junctions) file in the format of DEXSeq, but I am not able to figure out how the corresponding GTF file should be prepared.


ADD REPLYlink written 3.9 years ago by Anil Kesarwani80
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2502 users visited in the last hour