Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • help with differential gene expression with cufflinks and tophat

    Hi All,
    We are looking for differential gene expression of two different samples of the same species.

    After running tophat, for two different samples of the same species, we got accepted_hits.bam files for both. Then we ran cufflinks with the following command for both of the samples:

    $./cufflinks -G refGene.gtf accepted_hits1.bam
    $./cufflinks -G refGene.gtf accepted_hits2.bam
    where refGene.gtf is the main gtf file of the species downloaded from UCSC.

    This provided us with three files, genes.expr, transcripts.expr, transcripts.gtf (for two samples, all in all, we got 6 files, genes1.expr, transcripts1.expr, transcripts1.gtf, genes2.expr, transcripts2.expr, transcripts2.gtf)

    Then, we carried out cuffcompare with the following command but we are not quite sure as to this is the right way to do it or not.

    $./cuffcompare -r refGene.gtf /sample1/transcripts1.gtf /sample2/transcripts2.gtf (Is this how it should be?)

    Output: stdout.combined.gtf, stdout.loci, stdout.tracking

    Then we carried out cuffdiff with the following command but we are not quite sure as to this is the right way to do it or not.

    $./cuffdiff stdout.combined.gtf accepted_hits1.bam accepted_hits2.bam (Is this how it should be?)


    Now we want to do the differential gene expression for which we guess we need to use cuffcompare and cuffdiff, right? It will be highly appreciated if anyone can guide us on the same. Also, please tell us as to what files would be input for cuffcompare and cuffdiff to do the differential gene expression.

    Thank you very much for the helping hand.

  • #2
    Your commands appear to be correct.

    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
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    51 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    67 views
    0 likes
    Last Post seqadmin  
    Working...
    X