Question: How can I perform DE analysis for GRO-seq data between conditions?
gravatar for JJ
2.1 years ago by
JJ520 wrote:

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

rna-seq gro-seq • 1.0k views
ADD COMMENTlink modified 2.1 years ago by jared.andrews077.2k • written 2.1 years ago by JJ520
gravatar for jared.andrews07
2.1 years ago by
Memphis, TN
jared.andrews077.2k wrote:

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.

ADD COMMENTlink written 2.1 years ago by jared.andrews077.2k
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: 972 users visited in the last hour