Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Beginner Questions

    I'm a beginner to the field. I'm trying to align my reads to a reference genome. However, I'm not quite sure I understand what tophat is trying to do. In fact, I want to make sure my understanding to RNA sequencing is correct.

    Let's say,

    "tophat -p 8 -G genes.gtf -o C1_R1_thout genome C1_R1_1.fq C1_R1_2.fq"

    From my understanding, the *.fq files are raw sequence reads generated by experiments. Those raw reads are aligned to sequences defined in a reference genome. The reference genome in this example is "genome". genes.gtf lists some known transcripts.

    The point of doing it is to analyse unknown transcripts. Say, if some of my reads can't be aligned to the reference genome, we know it might mean that there's a new mutation. Am I correct?

    Another question, why do we need to provide an annotation file (genes.gtf) to tophat? Tophat already has the information from the reference genome to align with?

    Thanks,

  • #2
    Originally posted by tedwong View Post
    The point of doing it is to analyse unknown transcripts. Say, if some of my reads can't be aligned to the reference genome, we know it might mean that there's a new mutation. Am I correct?
    Tophat helps serve a few different and sometimes overlapping roles. Firstly, it will align the reads to known (annotated) transcripts. Secondly, it will align reads directly to the genome. Thirdly, it will attempt to align reads in such a way that finding novel (unannotated) transcripts is possible. For actually finding these, tools like cufflinks are useful.

    There are also many reasons a given read won't align to the the genome/transcriptome. Firstly, it could be from a contaminant (bacteria, etc.). This is extremely common. Secondly, it could be from a portion of the genome that doesn't exist in the reference sequence (e.g., the 45S rRNA sequence isn't in the mouse genome and there are a number of sequence gaps in most (all?) mammalian genomes). Thirdly, there may be variants in the sample that was sequenced that aren't present in the reference sequence (there's almost always some of this). Finally, this can be caused by sequencing errors.

    There are undoubtedly other reasons in addition to those I listed above!

    Another question, why do we need to provide an annotation file (genes.gtf) to tophat? Tophat already has the information from the reference genome to align with?
    Providing tophat with a reference annotation allows it to more easily align reads to the known transcriptome. Most of you reads will align either to that or to the genome, so this decreases the number of reads that need to be used to find novel splice sites. Finding these is somewhat computationally demanding and not without its issues.

    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
    11 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
    68 views
    0 likes
    Last Post seqadmin  
    Working...
    X