Question: Mask gtf file for cufflinks
0
gravatar for DVA
2.5 years ago by
DVA550
United States
DVA550 wrote:

Is the "-M" option necessary when using Cufflinks? I understand it would prevent downstream analysis focusing on regions like rRNA, but could this be done at the very end after differential expression? Also, I could not find a mask gtf to download from USCS genome browser... Any advice please? Thank you.

rna-seq • 900 views
ADD COMMENTlink written 2.5 years ago by DVA550
2

You should know that the old 'Tuxedo' pipeline of Tophat(2) and Cufflinks is no longer the "advisable" tool for RNA-seq analysis. The software is deprecated/ in low maintenance and should be replaced by HISAT2, StringTie and ballgown. See this paper: Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown.

ADD REPLYlink written 2.5 years ago by WouterDeCoster44k

Thanks a lot. I was trying to repeat a previous analysis done by a colleague.

ADD REPLYlink written 2.5 years ago by DVA550

Alright, makes sense, but probably also interesting to see if you can replicate the findings with a more recent method.

ADD REPLYlink written 2.5 years ago by WouterDeCoster44k
1

Yes absolutely. That is what I will do next. Thank you!

ADD REPLYlink written 2.5 years ago by DVA550

Do you know if my library is 3' rna seq library, what parameters should I alter with StringTie? Ideally, for 3' rna seq lib, transcript assembly is not needed for our DE analysis. Some default setting such as 200bps in -m, minimum length allowed for predicted transcripts,will be not relevant.

What I have in mind is the alternative workflow on : http://ccb.jhu.edu/software/stringtie/index.shtml?t=manual (second workflow under "Differential Expression Analysis).

ADD REPLYlink modified 2.5 years ago • written 2.5 years ago by DVA550

So you don't have full RNA-seq? Then I doubt whether stringtie is going to be much use, there is not a lot to assemble. I would suggest you use featureCounts or salmon for quantification and go directly to DESeq2 or edgeR.

ADD REPLYlink written 2.5 years ago by WouterDeCoster44k
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: 937 users visited in the last hour