0
$\begingroup$

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?

$\endgroup$
1
  • $\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
    Commented May 2, 2018 at 13:49

1 Answer 1

2
$\begingroup$

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.

$\endgroup$
1
  • $\begingroup$ Thank you Ryan. Yeah, I will use either DEXseq or rMATS in that case. $\endgroup$ Commented May 3, 2018 at 8:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

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