Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • valei
    Member
    • Sep 2010
    • 14

    #1

    TopHat2 with genes lacking introns

    Dear all,
    I'm working with RNA-seq paired-end data and I'm using TopHat-BowTie in order to align my reads against the Human reference genome.
    My problem is the following: I' have some reads coming from genes that are duplicated in the genome and they are very similar in terms of sequence, the only difference between them is that the second gene doesn't have the introns (for example this is the case of GLUD and GLUD2, or NACA and NACA2, etc..). When TopHat works with these reads it prefer align with mismatches in the contiguous sequence (GLUD2 or NACA2) instead of splitting the read in two parts and align it in two separated exon without any mismatch (GLUD or NACA).
    Obviously the correct alignment is the second one, I'm sure of this because I got the same mismatches in all the samples that I analyzed.
    How can I overcome this problem?
    Can I use some option or particular parameter? I would like to constrain TopHat to prefer no-mismatch alignments even if they are splitted in two of more segments.
    I hope my question is clear....
    Many many thanks for your help!
    Last edited by valei; 06-08-2012, 02:16 AM.
  • Dario1984
    Senior Member
    • Jun 2011
    • 166

    #2
    You have to upgrade to TopHat version 2. You supply your transcriptome annotation with the --GTF option. It maps the reads first to the spliced transcript models, which will take away your problem.

    The name for the genes that don't have introns and are the same sequence as ones that do is pseudogenes.

    Comment

    Latest Articles

    Collapse

    • SEQadmin2
      Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
      by SEQadmin2



      Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

      Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
      ...
      10-14-2026, 09:55 AM
    • SEQadmin2
      New Genomics Technologies Take Aim at Long-Standing Limits
      by SEQadmin2


      Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

      We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
      ...
      09-28-2026, 10:25 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by SEQadmin2, 10-22-2026, 11:13 AM
    0 responses
    30 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-22-2026, 10:48 AM
    0 responses
    22 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-15-2026, 10:52 AM
    0 responses
    39 views
    0 reactions
    Last Post SEQadmin2  
    Started by SEQadmin2, 10-14-2026, 05:39 AM
    0 responses
    45 views
    0 reactions
    Last Post SEQadmin2  
    Working...