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$ – llrs May 2 '18 at 13:49

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.

| improve this answer | |
  • $\begingroup$ Thank you Ryan. Yeah, I will use either DEXseq or rMATS in that case. $\endgroup$ – jagat shrestha May 3 '18 at 8:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.