SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Query regarding the cuffdiff command to get differential gene expression kumardeep RNA Sequencing 1 05-24-2012 09:07 AM
Does Cuffdiff normalize FPKMS before it does differential expression? greener Bioinformatics 0 02-25-2011 02:37 PM
Cufflinks/Cuffdiff significant differential expression memo Bioinformatics 5 01-25-2011 09:49 AM
Cuffdiff differential expression two groups jamminbeh Bioinformatics 0 12-09-2010 02:13 PM

Reply
 
Thread Tools
Old 02-05-2014, 10:12 AM   #1
gwilymh
Member
 
Location: Milwaukee

Join Date: Dec 2011
Posts: 72
Default Differential isoform expression in cuffdiff: how does it work?

I am running an analysis in CuffDiff using a GFF file of annotations from FlyBase (www.flybase.org). The GFF file contains annotation data for genes, mRNA, CDS and 5'- and 3'-UTR's.

Differential expression of isoforms is reported in some 1266 pairwise comparisons. However, most of the genes in the analysis contain only 1 transcript; and there is very little overlap between the genes showing differential expression (file gene_exp.diff) and the isoforms showing differential expression (file isoform_exp.diff). Can any one explain to me how this is possible?

My code is below:

cuffdiff \
-o "/gpfs/group1/f/flyinv/Outputs_CuffDiff/Exons_B_allFemales_Ch3/" \
-L AR,CH,PP,PX,ST,TL \
--total-hits-norm \
--dispersion-method per-condition \
--frag-bias-correct /gpfs/group1/f/flyinv/working_index/Dpse3_0_1.fa \
--multi-read-correct \
--library-norm-method geometric \
/gpfs/group1/f/flyinv/working_index/Ch3_FB_genes_exons_mRNA_CDS_UTR_PalmieriCh3Genes.gff3 \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/AR_DM1005_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/AR_KB945_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/AR_MSH126_Female/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/CH_JR272_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/CH_KB888_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/CH_MSH202_Female/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PP_DM1049_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PP_DM1054_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PP_DM1065_Female/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PX_SPE123_1-2_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PX_SPE123_4-1_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/PX_SPE123_5-2_Female/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/ST_JR138_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/ST_MSH177_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/ST_MSH49_Female/accepted_hits.bam \
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/TL_SCI_12-2_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/TL_SPE123_7-1_Female/accepted_hits.bam,\
/gpfs/group1/f/flyinv/Outputs_TopHat/transcriptomeSequence_exons/TL_SPE123_8-1_Female/accepted_hits.bam
gwilymh is offline   Reply With Quote
Reply

Tags
annotation, cuffdiff, gene expression, isoform

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


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