Cuffdiff Results ZERO / 0 FPKM for all genes

gopal_botany

Junior Member

Join Date: Jan 2014

Posts: 3
- Share
- Tweet
#1

Cuffdiff Results ZERO / 0 FPKM for all genes

01-11-2014, 02:22 AM

Dear All,

I am new to this community and please accept my greetings....

I have performed denovo assembly of plant transcriptome by using SOAP-denovo and I have used Tophat ---> Cufflink -----> Cuffdiff pipeline for downstream analysis.

Respective commands are given below -

tophat --coverage-search --microexon-search -r 200 -o /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/Flower /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/SAF_TOPHAT_REF /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Flower_adapter_trimmed_R1.fastq /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Flower_adapter_trimmed_R2.fastq

tophat --coverage-search --microexon-search -r 200 -o /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/Root /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/SAF_TOPHAT_REF /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Root_adapter_trimmed_R1.fastq /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Root_adapter_trimmed_R2.fastq

tophat --coverage-search --microexon-search -r 200 -o /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/Shoot /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/SAF_TOPHAT_REF /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/shoot_adapter_trimmed_R1.fastq /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/shoot_adapter_trimmed_R2.fastq

tophat --coverage-search --microexon-search -r 200 -o /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/Leaf /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/SAF_TOPHAT_REF /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Leaf_adapter_trimmed_R1.fastq /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Leaf_adapter_trimmed_R2.fastq

tophat --coverage-search --microexon-search -r 200 -o /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/Leaf /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Tophat/SAF_TOPHAT_REF /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Leaf_fastx_trimmer_R1_adapter_trimmed.fq /home/gopal/SAF/Transcriptome/New_Analysis_24_Dec_2013/Leaf_fastx_trimmer_R2_adapter_trimmed.fq

cufflinks accepted_hits.bam

cuffmerge path.txt

cuffdiff -N -p 40 -o ../cuffdiff/new/ merged.gtf accepted_hits_FL.bam accepted_hits_LF.bam accepted_hits_RT.bam accepted_hits_ST.bam

My query is why cuffdiff results zero FPKM in all genes and each sample, while FPKM values were there in clufflink results.

Please help me.....

Thanking you.
Tags: None
TiborNagy

Senior Member

Join Date: Mar 2010

Posts: 329
- Share
- Tweet
#2

01-13-2014, 05:02 AM

The parameters looks OK. Maybe you mixed accepted_hits_*.bam in the command line. It is hard to tell, because we do not know how do you created accepted_hits_FL.bam, accepted_hits_LF.bam, accepted_hits_RT.bam and accepted_hits_ST.bam
Comment
gopal_botany

Junior Member

Join Date: Jan 2014

Posts: 3
- Share
- Tweet
#3

01-13-2014, 10:47 PM

Thanks TiborNagy,

As I written above the step wise step execution of pipeline... it will help you to understand how these files generated.

Thanks
Comment

Previous template Next

From Collection to Sequencing: Why Sample Preparation and Preservation Define Sequencing Data

by SEQadmin2

Data variability is still an issue in sequencing technologies despite the advances in reproducibility and accuracy of these platforms. But the problem does not originate in the sequencing itself, but in the previous steps, before the sample reaches the sequencer.

The first step is collection, followed by preservation and sample preparation for analysis. Most scientists overlook those steps, but not being careful might just be skewing the experiment’s results.
...
- Channel: Articles
06-02-2026, 10:05 AM
Single-Cell Sequencing at an Inflection Point: Early Impacts of New Platforms and Emerging Trends

by SEQadmin2

With the launch of new single-cell sequencing platforms in 2026, the field stands at an exciting inflection point. This article surveys the most impactful advances in the field and discusses how they’re reshaping research in cancer, immunology, and beyond.

Introduction

Single-cell sequencing technologies have undergone remarkable advances over the past decade, transitioning from low-throughput experimental approaches to highly scalable platforms capable of...
- Channel: Articles
05-22-2026, 06:42 AM

Topics	Statistics	Last Post
Sequencing the Two-Toed Sloth Genome Reveals Jumping Genes Tied to Its Extreme Metabolism by SEQadmin2 Started by SEQadmin2, 06-09-2026, 11:58 AM	0 responses 19 views 0 reactions	Last Post by SEQadmin2 06-09-2026, 11:58 AM
A New Method Makes Hantavirus Genome Analysis Faster and More Accessible by SEQadmin2 Started by SEQadmin2, 06-05-2026, 10:09 AM	0 responses 27 views 0 reactions	Last Post by SEQadmin2 06-05-2026, 10:09 AM
A New Single-Cell Method Maps DNA-Protein Interactions by SEQadmin2 Started by SEQadmin2, 06-04-2026, 08:59 AM	0 responses 38 views 0 reactions	Last Post by SEQadmin2 06-04-2026, 08:59 AM
Long-Read RNA Sequencing Uncovers a Hidden Layer of Immune Cell Regulation by SEQadmin2 Started by SEQadmin2, 06-02-2026, 12:03 PM	0 responses 61 views 0 reactions	Last Post by SEQadmin2 06-02-2026, 12:03 PM

Unconfigured Ad

Cuffdiff Results ZERO / 0 FPKM for all genes

Comment

Comment

Latest Articles

ad_right_rmr

News