Question: Differential expression analysis suggestions needed
0
gravatar for sierraallinone
2.6 years ago by
sierraallinone20 wrote:

We have a trinity transcriptome assembly for 6 muscles types, 3 of which have triplicate data. We generated a supertranscript but now I’m reading conflicting things on how to approach the differential analysis portion of this experiment.

We have no reference genome so if we go with HISAT2 —> cufflinks —> Cuffdiff —> cummerbund wed of course use the supertranscript as the reference for which the reads align back to and go from there.

We also toyed with the idea of using DESEq/EdgeR

Then it was suggested to do hisat2–> stringtie —> ballgown but then someone says that those programs are good for genome assisted analyses while salmon/kallisto —> sleuth is transcript level analysis. I’ve looked into this type of pseudo alignment but at this point I just really don’t know which pipeline to use. Furthermore, some people think transcript level analysis is better than gene level analysis but it seems that gene level analysis is more reliable.

rna-seq • 972 views
ADD COMMENTlink modified 2.6 years ago • written 2.6 years ago by sierraallinone20

https://cgatoxford.wordpress.com/2016/08/17/why-you-should-stop-using-featurecounts-htseq-or-cufflinks2-and-start-using-kallisto-salmon-or-sailfish/

ADD REPLYlink written 2.6 years ago by Arup Ghosh2.8k
1

Hmm interesting.. a hybrid approach. I noticed they mentioned the alignment independent method works really well for unique transcripts. Just how unique? Because my transcripts are most likely going to be pretty similar with a few different alternative splicing zones of course it’s hard to tell since there’s no reference genome and we really only have one muscle type sequences from one clone as our initial rna reference. So does that mean I should maybe stick to the cufflinks pipeline?

ADD REPLYlink written 2.6 years ago by sierraallinone20

You should avoid using cufflinks & tophat, any other analysis flow will be fine.

ADD REPLYlink written 2.6 years ago by Arup Ghosh2.8k
  1. Generate the reference transcriptome from the data using HISAT2 and StringTie (or Trinity)
  2. Use this reference transcriptome with your original FASTQs for the purposes of determining read count abundances over your identified transcripts using Kallisto / Salmon
  3. conduct differential expression analysis sing EdgeR / DESeq2

...or, from #1, use prepDE.py (bundled with StringTie, I believe) to immediately generate count data suitable for #3

ADD REPLYlink modified 2.6 years ago • written 2.6 years ago by Kevin Blighe70k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2895 users visited in the last hour
_