This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Cufflinks On Targeted Transcripts

From RNASeq data, I'm interesting in investigating the differential expression of transcripts of only a small number of genes (a hypothesis-based approach rather than hypothesis-free looking across the whole transcriptome).
Would it be a valid approach to run cufflinks using a gtf file containing only the 10-20 transcripts I'm interested in and then cuffdiff on the output?

Thanks, Dave

cufflinks rnaseq rna-seq

1 answer

Cufflinks can assemble novel transcripts, so make sure to also supply the -G option. One thing to think about is the size normalization and dispersion estimation steps in cuffdiff. With only 10-20 transcripts, I would worry that at least the library size normalization step wouldn't work very well. If you have enough replicates that you don't need to share information across genes, then the dispersion calculations should be fine. The simpler method, btw, would be to just use the full dataset for all of the estimates and then just subset it to the transcripts of interest afterward (you can then readjust the p.values).

Log in to answer this question.