Unconfigured Ad

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts
  • gwilymh
    Member
    • Dec 2011
    • 72

    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

Latest Articles

Collapse

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 06-05-2026, 10:09 AM
0 responses
11 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-04-2026, 08:59 AM
0 responses
23 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 12:03 PM
0 responses
28 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 06-02-2026, 11:40 AM
0 responses
22 views
0 reactions
Last Post SEQadmin2  
Working...