How can I perform DE analysis for GRO-seq data between conditions?
Entering edit mode
2.6 years ago
JJ ▴ 560

Hi all,

So GRO-seq is a clever method to measure nascent RNA and I have found a nice protocol to look for active enhancers.

However, it also measures nascent mRNA transcription. So now I am wondering how to do a differential expression (DE) analysis for genes between conditions. Can I simply use standard DE analysis tools? This has been done for GRO-seq apparently (from the suppl):

The hg19 RefSeq gene list was used for all transcription level analysis. Htseq was used to get gene level read counts number from Bowtie mapped bam files. The resultant gene read count table was subjected to DESeq2 for differential gene analysis

But what confuses me is that with GRO-seq one gets mainly incomplete transcripts, right? So is the standard DE analysis still an correct approach? Would it be better to use the transcripts found by e.g., groHMM or HOMER and feed these in a DE analysis?

Input is very much appreciated! Thanks

GRO-seq RNA-Seq • 1.2k views
Entering edit mode
2.6 years ago

If it's been published (hopefully more than once), you likely won't get any complaints about it. Plus DESeq2 is easy/quick enough to run that it's worth a shot just to see if the results make sense.


Login before adding your answer.

Traffic: 2789 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6