Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bug? duplicated genes in cufflinks output genes.expr silin284 Bioinformatics 3 05-18-2014 12:19 AM
Cufflinks, differentially expressed genes statsteam Bioinformatics 5 11-15-2013 12:28 PM
Cufflinks cufflinks v1.0.3 - segmentation fault bias correction chrNT annotations adrian Bioinformatics 0 06-08-2011 02:28 PM
Cufflinks merging nearby genes xinchen Bioinformatics 2 01-24-2011 07:28 PM
inconsistent gene names in genes.expr - Cufflinks Boel Bioinformatics 2 04-14-2010 06:16 AM

Thread Tools
Old 11-19-2010, 02:37 AM   #1
Location: France

Join Date: Sep 2010
Posts: 11
Default Cufflinks merges adjacent genes

Hi all,

I am working on the RNA-Seq of an organism with very short intergenic distance (e.g. 50bp). I am using the TopHat/Cufflinks package for my analysis, and I found out Cufflinks tends to merge adjacent genes as a single trasfrag, when the intergenic distance is short. I wonder whether there's an option to adjust the cutoff for the coverage gap size for the definition of a transfrag.

Thanks in advance.

proteomania is offline   Reply With Quote
Old 11-20-2010, 03:58 PM   #2
Location: Pasadena, CA

Join Date: May 2009
Posts: 45

There's a coverage option for filtering out transcripts with "retained introns" that are assembled due to complete coverage of the intron with reads from pre-spliced mRNAs. However, this isn't going to work in your case as the intergenic regions aren't going to be spanned by junctions and I don't think there is any option that does what you want.

Genes that are on different strands should be assembled relatively correctly, as there the orientation of the junctions will tell Cufflinks what to do. But for genes in the same orientation, you are out of luck, especially if they are at about the same expression levels.

What you could do is look at you transcript models and try to find those that have two tandem ORFs, those are likely candidates for merged genes. Also, you could try doing ChIP-Seq for histone marks associated with promoters or preinitiation complex components.
GKM is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 08:56 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO