Entering edit mode
3.8 years ago
yaghoub.amraei
▴
10
Hi all. I have 24 samples RNA-seq, paired-end data including 2 genotypes, 3 Zone from root, 2 conditions and each condition including two replicate. How can I do cuffdiff for this expriments. thanks alot
Thanks, Can't use Caddiff with this number of samples?
Yes, you can, but again, that is deprecated
Were did you hear/read that Cufflinks/Cuffdiff is deprecated?
The authors of those tools say so themselves (see their lab websites; their twitters; their previous posts on biostars; etc.).
I have not been able to find it. Could you point to something more recent than this (which states cufflinks it is not deprecated).
Sure.
Most recent: This mentions it as "retired software" and "no longer under active development and are not being supported". It does not list it as "now maintained/developed elsewhere".
Older posts: See this tweet (referring to it as obsolete). And see this.
I would also not use this pipeline. Stringtie is overly complicated unless you want to assemble the transcriptome. For a standard DE analysis (i a proper reference transcriptome exists) there is (imho) no need to bother with it. ALso leightweight pseudo- or selective aligners such as kallisto and salmon typically outperform traditional alignment for RNA-seq quantification (check any of the recent benchmark papers). This ballgown tool does not really provide features which other, better documented and actively maintained tools such as DESeq2 or edgeR do not offer.
I would go with quantification by salmon, importing data into R with tximport, then DESeq2, as described e.g. here: https://www.bioconductor.org/packages/devel/workflows/vignettes/rnaseqGene/inst/doc/rnaseqGene.html
my goal is LncRNA detection...
Recommendations for lncRNA: https://academic.oup.com/gigascience/article/8/12/giz145/5663671
Still recommend pseudoalignment, e.g. kallisto, for this.