This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Does Cuffmerge discard single-exon transcripts?

I'm using cuffmerge to combine the transcripts built using cufflinks on different samples. I've just noticed that the annotation I get from cuffmerge contains no transcripts with only a single exon, they all have at least 2. Looking at the cufflinks output for the individual libraries, they do contain transcripts with a single exon. Does cuffmerge just automatically discard any single-exon transcripts, or there some other metric it uses and it just happens that none of the single-exon transcripts pass the metric? I don't see any options for cuffmerge to adjust this behavior.

cufflinks rna-seq
• 2,599 views • link
updated by Satyajeet Khare Satyajeet Khare • written by colin.kern colin.kern

1 answer

cuffmerge does not exclude single exon transcripts. You can check canonical histone genes which are mostly mono-exonic. In merged.gtf file, they should have only one exon marked as "exon 1".

• 1 views • link
written by Satyajeet Khare Satyajeet Khare

You are correct, those genes are there. I forgot to mention that I have filtered out all transcripts with class_code "=" since I'm interested in finding novel long non-coding RNA. This leaves me with no single-exon transcripts. Could this be because cuffmerge "automatically filters a number of transfrags that are probably artfifacts"?

• 1 views • link
written by colin.kern colin.kern

Yes, cuffmerge may have filtered them out. Btw, many long non-coding RNAs are polyexonic. And therefore you may find novel lncRNAs in your sample.

• 1 views • link
written by Satyajeet Khare Satyajeet Khare

Log in to answer this question.

More posts like this