Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Is it possible to pull out exon sequences from cuffmerge out? TPH Bioinformatics 0 03-01-2016 12:42 PM
Cuffmerge Exon vs Transcript output from Cufflinks Problem. tonup69 Bioinformatics 0 02-18-2016 07:58 AM
Cuffmerge outputs only exon annotation data gwilymh Bioinformatics 0 02-03-2015 10:06 AM
cuffmerge: same tss_id assigned to transcripts with different starting sites mra Bioinformatics 0 04-02-2014 11:42 PM
cuffmerge tss_id gik Bioinformatics 0 03-25-2013 08:22 AM

Thread Tools
Old 02-04-2017, 08:56 PM   #1
Senior Member
Location: 63130

Join Date: Apr 2012
Posts: 125
Default cuffmerge:the same tss_id, but different leftmost exon

Hello, there,

I did de novo RNA-seq cufflinks assembly. I then merged transcripts.gtf files from different replicates to obtain merged.gtf in the absence of reference annotation.

This merged.gtf showed that the transcripts with the same tss_id have different leftmost exons:

Here is the example:

chr2 Cufflinks exon 25289899 25290661 . . . gene_id "XLOC_005458"; transcript_id "TCONS_00010739"; exon_number "1"; oId "CUFF.5451.1"; tss_id "TSS7438";
chr2 Cufflinks exon 25290738 25290883 . . . gene_id "XLOC_005458"; transcript_id "TCONS_00010739"; exon_number "2"; oId "CUFF.5451.1"; tss_id "TSS7438";
chr2 Cufflinks exon 25290976 25291190 . . . gene_id "XLOC_005458"; transcript_id "TCONS_00010739"; exon_number "3"; oId "CUFF.5451.1"; tss_id "TSS7438";
chr2 Cufflinks exon 25289938 25290082 . . . gene_id "XLOC_005458"; transcript_id "TCONS_00010740"; exon_number "1"; oId "CUFF.5451.2"; tss_id "TSS7438";
chr2 Cufflinks exon 25290388 25291177 . . . gene_id "XLOC_005458"; transcript_id "TCONS_00010740"; exon_number "2"; oId "CUFF.5451.2"; tss_id "TSS7438";

this gene: XLOC_005458, has two transcript isoforms: TCONS_00010739 and TCONS_00010740, and those two isoforms starts with the same transcription start site since the tss_ids are the same: TSS7438

However, transcript TCONS_00010739 apparently starts at 25289899 (exon1 left coordinate) on chr2, and transcript TCONS_00010740 starts at 25289938 (exon1 left coordinate) on chr2

How should I explain this? I mean if the tss_ids are the same, their leftmost exons should start at the same coordinates.

Thank you very much.

capricy 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 11:51 PM.

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