SEQanswers (
-   Bioinformatics (
-   -   cufflinks : analysis comparison with and without a gtf reference file (

sohnic 06-09-2011 06:54 AM

cufflinks : analysis comparison with and without a gtf reference file
I have many questions about cufflinks output. Here one of them :
First I've used tophat to map my RNAseq (100pb) to obtain a accepted_hits.bam file.
Then I've used cufflinks in two ways :
- simply : cufflinks accepted_hits.bam
- with a gtf file, that is the actually annotation of my genome (eucalyptus) :
cufflinks -g annotGenome.gtf accepted_hits.bam
Note that I've used the g and not the G option.

One example result :
- without reference gtf :
one gene / one isoform : 12110-17714
The first part of this isoform 12-16530 has the same structure intron/exon than isoforms formed with the reference.
Then I have a last exon 16530-17714.

-with reference gtf
one gene/two isoforms
+ transcript 1: 12024-17350 = exact transcript from the reference
full_read_support "no";
The corresponding no reference last exon is now :
16530 - 16561
16597 - 17350
That's my reference, but in my run this intron is mapped. There is no read that split in two parts. A few reads begin at position 16595. I've checked no read ending at 16561.
I thing this RNA doesn't exist in my transcriptome.
+ transcript 2 : 12024-17714
full_read_support "yes";
last exon : 16530-17714, the same exon than the no reference version
Why this transcript2 contains the 12024-12109 portion that is not mapped with RNAseq (instead the reference=transcript1 begin with this sequence) ?

for the two isoforms, I have FPKM values (4 for transcript1 that doesn't seem to exist in my transcriptome and 13 for the transcript2). How cufflinks attributes those values ?

With the version without gtf reference, I have a FPKM=36, that is the double comparing with the version with reference (13+4=17) while the mapping file is the same.

At least, note that those transcripts are located on the forward strand of the genome and that there is nothing in gtf and cufflink results on the opposite strand at this location.

Many thanks for your suggestion,


Stefano Manzini 12-17-2013 01:45 AM

Hello sohnic,

I have no answer to your question, but I would like to ask you one because you're posting about the very same issue I am wondering about.

Can you just give me some hints about the difference of using (or not) a local reference to run cufflinks?
I would like to understand what does cufflinks use to assembly the information contained in target.bam file when you don't provide it a local reference. Further, I'd also like to know whether using it will make cufflinks name the transcripts with the gene names instead of CUFF.65279 and the like.

Thank you for your help.

ahalfpen 11-03-2016 06:28 PM

My experience with cufflinks
Cufflinks will use the mapped reads and the reference genome to create a parsimonious assembly (a minimum spanning tree-esque structure) that explains the reads. The logic of the algorithm is based on the Burset and Guigo paper (Burset, Guigo, Evaluation of Gene Structure Prediction, 1996) and their observations regarding high false positive rates of gene finding programs that are guided by annotation. This is my understanding of the matter based on the past year I have spent working on similar issues.

brojee 07-07-2019 06:40 AM

Sleeve fasteners will utilize the mapped peruses and the reference genome to make a closefisted get together (a base spreading over tree-esque structure) that clarifies the peruses. The rationale of the calculation depends on the Burset and Guigo paper (Burset, Guigo, Evaluation of Gene Structure Prediction, 1996) and their perceptions in regards to high false positive rates of quality discovering programs that are guided by explanation. This is my comprehension of the issue dependent on the previous year I have spent taking a shot at comparable issues.

All times are GMT -8. The time now is 12:10 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.