SEQanswers

Go Back   SEQanswers > Applications Forums > RNA Sequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
Cufflinks 0 FPKM values herstein Bioinformatics 2 07-24-2013 10:21 PM
Cufflinks - Odd MIR FPKM values. ramma RNA Sequencing 2 07-16-2013 08:23 AM
Different FPKM values of cufflinks and cuffdiff mrfox Bioinformatics 5 10-17-2012 01:10 PM
Cufflinks and cuffdiff FPKM values combiochem Bioinformatics 12 10-13-2012 11:37 PM
cufflinks-1.0.3 produces very high FPKM values when compared to cufflinks-0.9.3. Why? pinki999 Bioinformatics 5 06-09-2012 06:48 AM

Reply
 
Thread Tools
Old 12-31-2013, 01:07 PM   #1
gwilymh
Member
 
Location: Milwaukee

Join Date: Dec 2011
Posts: 72
Default Cufflinks returning elevated FPKM values for 'new' transcripts

I am running RNA-Seq in a Tuxedo pipeline. In the final Cufflinks output, I am getting massively elevated FPKM values for some transcripts. These inflated values are only returned for the transcripts newly discovered by Cufflinks/TopHat (i.e. transcripts that were not previously annotated). The transcripts from annotated genes seemed normal.

In one analysis, for example, the known genes had an average FPKM value of 214.5, and a maximum FPKM value of 147,473. The newly discovered transcripts, however, returned a mean FPKM of 139,234.6 and a maximum of 74,769,200; 2.5% of the new transcripts had vales greater than the maximum FPKM detected for the annotated genes. The results for the new transcripts clearly contain artifacts.

My pipeline is below
tophat \
--min-anchor-length 10 \
--splice-mismatches 1 \
--min-intron-length 5 \
--microexon-search \
--fusion-search \
--transcriptome-index=/gpfs/group1/f/flyinv/working_index/Dpse3_0_1_exons_ \
-o "/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/AR_MSH126_Male" \
/gpfs/group1/f/flyinv/working_index/Dpse3_0_1 \
/gpfs/group1/f/flyinv/RNASeq/AR_MSH126_Male_1_sequence.txt \
/gpfs/group1/f/flyinv/RNASeq/AR_MSH126_Male_2_sequence.txt \

cufflinks \
--output-dir "/gpfs/group1/f/flyinv/Outputs_CuffLinks/transcriptomeSequence_exons/AR_MSH126_Male" \
--GTF-guide /gpfs/group1/f/flyinv/working_index/Dpse3_0_Exons.gff3 \
--upper-quartile-norm \
--min-intron-length 5 \
--quiet \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/AR_MSH126_Male/accepted_hits.bam \


Does anyone know how or why this occurs? How can it be prevented? How can such artifacts be screened out of downstream analyses?
gwilymh is offline   Reply With Quote
Reply

Tags
cufflinks, fpkm, rna-seq, tophat

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:13 PM.


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