Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • drunk_coder
    Junior Member
    • Dec 2011
    • 9

    #1

    cuffdif bug?

    Recently, our lab followed tophat->cuffdif protocol to analyze RNA-seq gene expression, when we look into the cuffdif gene_exp.diff output, we find a problem that
    one gene with low reads mapped has extremely high FPKM while other genes' expression values seem fine.

    Here is the Detail:
    1 RNA-Seq library: a) A1 (100bp paired, no replicate) b) C1-2 (100bp paired,no replicate)
    2 Genome data: Homo_sapiens.GRCh37.69.dna.toplevel.fa
    3 GTF : Homo_sapiens.GRCh37.69.gtf
    4 Gene with above problem in cuffdif "gene_exp.diff" file
    test_id gene_id gene locus sample_1 sample_2 status value_1 value_2 log2(fold_change) test_stat p_value q_value significant
    ENSG00000200170 ENSG00000200170 Y_RNA 10:73856277-73995472 A1 C1-2 OK 4.75944e+06 0 -inf 0 0.1507 0.999526 no


    5 Gene with above problem in cuffdif "gene.count_tracking" file
    tracking_id A1_count A1_count_variance A1_count_uncertainty_var A1_count_dispersion_var A1_status C1-2_count C1-2_count_variance C1-2_count_uncertainty_var C1-2_count_dispersion_var C1-2_status
    ENSG00000200170 1.19092 31 0 30.2743 OK 0 1 0 0 OK


    6 tophat version: 2.0.4
    7 cufflinks version: 2.0.0 or 2.10 (similar result)
    8 Commands: (follow the nature protocol):
    tophat: tophat-2.0.4.Linux_x86_64/tophat -r 500 -o /tophatmapping/A1_Remapping -p 8 -G /tophatmapping/gtf/Homo_sapiens.GRCh37.69.gtf /tophatmapping/index/Homo_sapiens.GRCh37.69.dna.toplevel /dou_reads/A1_1.fq /dou_reads/A1_2.fq
    tophat-2.0.4.Linux_x86_64/tophat -r 500 -o /tophatmapping/C1_Remapping -p 8 -G /tophatmapping/gtf/Homo_sapiens.GRCh37.69.gtf /tophatmapping/index/Homo_sapiens.GRCh37.69.dna.toplevel /dou_reads/C1_1.fq /dou_reads/C1_2.fq

    cuffdif: ~/bin/cufflinks-2.1.0.Linux_x86_64/cuffdiff -o /tophatmapping/A1_C1-2_Recuffdiff/diff_out -b /tophatmapping/index/Homo_sapiens.GRCh37.69.dna.toplevel.fa -p 8 -L A1,C1-2 -u /tophatmapping/gtf/Homo_sapiens.GRCh37.69.gtf /tophatmapping/output_A1/accepted_hits.bam ~/fudan/maobin2/tophatmapping/C1_Remapping/accepted_hits.bam

    Question:
    1 So how can a gene with only 1 fragment mapped has 4.75944e+06 FPKM?
    2 Locus of this gene (ENSG00000200170) in gene_exp.diff file is different from locus in Homo_sapiens.GRCh37.69.gtf ,how could this happen ?

Latest Articles

Collapse

  • kinpict.com
    Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by kinpict.com
    The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
    Today, 06:37 PM
  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2





    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
    10-30-2026, 06:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 11-12-2026, 07:59 AM
0 responses
27 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-05-2026, 09:25 AM
0 responses
35 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-03-2026, 11:54 AM
2 responses
56 views
0 reactions
Last Post kinpict.com  
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
175 views
0 reactions
Last Post SEQadmin2  
Working...