Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Advice Needed on Tuxedo Suite

    Hi Everyone,

    I am working on a project with two groups (A & B) each containing 5 samples. Within each sample I have two RNA-seq FASTQ reads which I downloaded. So as of now, I am planning to run tophat for each of the samples (for a total of 10 accepted_hits.bam files). From here I plan on running cufflinks for each sample with its respective .bam file to produce a total of 10 transcripts.gtf files.

    My goal is to compare the differential gene expression of group A and group B. I am using the hg19 reference genome.

    I have two questions. First, is there a way I can run tophat/cufflinks on all the samples all at once or do I have to run one sample, then wait until its done, then do the next sample and so on?

    Second, once I have the transcripts.gtf files for each of the samples for groups A & B, is it necessary to merge them into one single .gtf file? I assume this would be done with cuffmerge right? What exactly does cuffmerge do and what are the advantages in using it. In addition, is there a need to use cuffcompare since I am comparing the differential gene expression of group A and group B or can I go straight to cuffdiff.

    In general, my second question is essentially asking the pipeline for using software after tophat & cufflinks based on what I am trying to do. It would be great if someone could give me a suggestion as to the order of how I should approach this.

    THANKS!

  • #2
    If you have a cluster or single computer with sufficient RAM and cores then you can run as many instances of tophat/cufflinks simultaneously as you want.

    If you don't use the -G option, then not using cuffmerge would cause you to compare apples and oranges (suppose one sample had a modified gene model and another didn't, the resulting FPKM comparison would be meaningless). If you use cuffmerge then you don't need to use cuffcompare. I think they had a methods paper with a step-by-step walk-through of an analysis, which I'd recommend you read through.

    Comment


    • #3
      Is part of your objective to identify transcripts which have not yet been annotated in the human genome? If not then there is really no reason to use Cufflinks at all.

      To perform DE analysis between groups A & B considering only already annotated genes in hg19 I would:

      1. Download the hg19 reference set from the Illumina iGenomes site.

      2. Prepare the known genes bowtie index as described in the TopHat manual.

      3. Align your reads to hg19 with Tophat 2.x.x using the "-T" option (align to transcriptome only).

      4. Count the reads aligned to genes from each of your 10 accepted_hits.bam (e.g. with feature counts or htseq-count).

      5. Use the raw read counts matrix as input to your preferred DE module (e.g. EdgeR, DESeq2).

      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