Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • NOTEST everywhere with cuffdiff

    Hello, I am running the following pipeline:

    tophat -p 4 -G genes.gtf -o sample1TopHat genome sample1R1.fastq sample1R2.fastq
    tophat -p 4 -G genes.gtf -o sample2TopHat genome sample2R1.fastq sample2R2.fastq

    cufflinks -o sample1Cufflinks -p 4 -G genes.gtf sample1TopHat/accepted_hits.bam
    cufflinks -o sample1Cufflinks -p 4 -G genes.gtf sample1TopHat/accepted_hits.bam

    cuffmerge -o sampleCuffMerge -p 8 -g genes.gtf -s genome.fa assemblylist.txt

    cuffdiff -o sampleCuffdiff -p 8 -b genome.fa -L s1,s2 -u sampleCuffMerge/merged.gtf sample1TopHat/accepted_hits.bam sample2TopHat/accepted_hits.bam

    However, I am getting NOTEST on the splicing.diff, whereas the isoform_exp.diff does show some information. I am just not sure how to debug this. Any suggestions?

    I have messed with the -c flag in cuffdiff. I am using cufflinks 2.2.1, and I am using the UCSC Homo Sapiens HG19 genome and annotation (genes.gtf). Any suggestions appreciated.

    Thank you in advance,
    Ramiro
    Last edited by ramiro2k; 06-22-2016, 09:18 AM. Reason: Added more information on software used.

  • #2
    So, you're expecting significant results from an experiment with no replicates?

    Also see : https://www.biostars.org/p/163535/

    Comment


    • #3
      Ran it with replicates with the same result

      Hello,

      Thank you so much for your reply. I ran it with two replicates and got a similar result, nothing significant, nothing expressed in the splicing.diff. I had looked at the link you mentioned in the past but I am surprised that nothing comes out as expressed or significant. On the differentialy expressed results I do get many genes as significant, etc. I am guessing this means that no splicing events ocurred, yet if you or anyone has any suggestions about anything I might be doing wrong please let me know.

      Thank you,
      Ramiro

      tophat -p 4 -G genes.gtf -o sample1r1TopHat genome sample1r1R1.fastq sample1r1R2.fastq
      tophat -p 4 -G genes.gtf -o sample1r2TopHat genome sample1r2R1.fastq sample1r2R2.fastq
      tophat -p 4 -G genes.gtf -o sample2r1TopHat genome sample2r1R1.fastq sample2r1R2.fastq
      tophat -p 4 -G genes.gtf -o sample2r2TopHat genome sample2r2R1.fastq sample2r2R2.fastq

      cufflinks -o sample1r1Cufflinks -p 4 -G genes.gtf sample1r1TopHat/accepted_hits.bam
      cufflinks -o sample1r2Cufflinks -p 4 -G genes.gtf sample1r2TopHat/accepted_hits.bam
      cufflinks -o sample2r1Cufflinks -p 4 -G genes.gtf sample2r1TopHat/accepted_hits.bam
      cufflinks -o sample2r2Cufflinks -p 4 -G genes.gtf sample2r2TopHat/accepted_hits.bam

      cuffmerge -o sampleCuffMerge -p 8 -g genes.gtf -s genome.fa assemblylist.txt

      cuffdiff -o sampleCuffdiff -p 8 -b genome.fa -L s1,s2 -u sampleCuffMerge/merged.gtf sample1r1TopHat/accepted_hits.bam,sample1r2TopHat/accepted_hits.bam sample2r1TopHat/accepted_hits.bam,sample2r2TopHat/accepted_hits.bam

      Comment


      • #4
        I actually found the answer in another post from 2011:

        One of the updates going from version 1.3 to version 2 (quoting from the cufflinks website here) was:

        "After substantial performance testing, we have determined that the false discovery rate of Jensen-Shannon-based tests (differential splicing, CDS switching, and promoter switching) can be unnacceptably high when used with fewer than three replicates in the conditions being compared. Cuffdiff now refrains from performing significance tests when one of the conditions involved has fewer than three replicates. You can change this behavior with the new --min-reps-for-js-test option. Cuffdiff still produces splicing.diff, cds.diff, and promoters.diff regardless of how many replicates you have. These files will include the JS distance scores, but none of the genes will be marked significant if you have fewer than the required number of replicates."

        This explained the situation.

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Strategies for Sequencing Challenging Samples
          by seqadmin


          Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
          03-22-2024, 06:39 AM
        • seqadmin
          Techniques and Challenges in Conservation Genomics
          by seqadmin



          The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

          Avian Conservation
          Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
          03-08-2024, 10:41 AM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, Yesterday, 06:37 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, Yesterday, 06:07 PM
        0 responses
        8 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-22-2024, 10:03 AM
        0 responses
        49 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 03-21-2024, 07:32 AM
        0 responses
        66 views
        0 likes
        Last Post seqadmin  
        Working...
        X