I have used Deseq for estimating the differentially expressed introns for the RNA seq data using the UCSC intron bed file. can Deseq be used for differentially expressed intron analysis?
-
$\begingroup$ Welcome to the site. You write Deseq in the question but you tagged with deseq2 (of the DEseq2 package). Which one do you really use? Also if you did it can be used, but it seems like you have problems with the results or something, if this is the case, you will need to provide more of information like the relevant part of the code used $\endgroup$– llrsCommented May 2, 2018 at 13:49
1 Answer
Typical tools to use for this would be the DEXseq packages (as done here) or rMATS, which will calculate intron retention rates.
I would generally suggest not summing intronic counts by gene and using that in DESeq2 (I assume that's what you meant by DESeq, since there's no reason to use the original DESeq for anything). What you'd be getting if you did this is less expressed introns and more sample-specific bias in enriching for mature versus pre-mRNA.
-
$\begingroup$ Thank you Ryan. Yeah, I will use either DEXseq or rMATS in that case. $\endgroup$ Commented May 3, 2018 at 8:16