Usage of prepDE.py with stringtie and deseq2
1
0
Entering edit mode
Nancy ▴ 10
@nancy-14435
Last seen 5.4 years ago

 I have been following the usage of prepDE.py script to convert stringtie quantifications to counts to input to DESeq2, exactly as described here: [http://ccb.jhu.edu/software/stringtie/index.shtml?t=manual#deseq]

We typically run Stringtie in the simple reference-based mode (using the –eB arguments, bypassing the merging and assembling steps, simply quantifying genes found in the reference). We then supply to prepDE.py the abundances output from this to get counts.

Over a series of analyses, performed on different data sets, different experiments, etc. I have consistently found the results to show very few, if any, differential expressed genes, and those that I do identify are of very small magnitude fold-change (no greater than 2 - 2.5 fold)

Typical experimental set ups comprise of a minimum or 3 replicates, bulk-tissue RNA-seq and these are in cell-lines or animal models. (human/mouse/rat). Inter-replicate variation (looking at individual gene fpkms or counts) does exist in some cases, whereas replicates are more tight in other cases.

While this is possible in the biological sense - most of these are discovery experiments so we don't have a "positive control" - I am concerned whether the analytical pipeline is responsible for what we see. (I am not entirely clear how prepDE works, to compute counts from coverage - could there be any over-normalizing, or ?? ).

I am curious if there are others who use this approach routinely and whether there are any other things to pay attention to, that I may not have considered.

Thanks for any inputs.

 

prepDE STRINGTIE deseq2 • 4.5k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 14 hours ago
United States

 I don’t know of any reason why you would have decreased sensitivity with this pipeline. You can also try asking the StringTie authors here:

http://ccb.jhu.edu/software/stringtie/index.shtml#contact

 

ADD COMMENT

Login before adding your answer.

Traffic: 595 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6