Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • 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.

  • #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

    • seqadmin
      Essential Discoveries and Tools in Epitranscriptomics
      by seqadmin


      The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
      Yesterday, 07:01 AM
    • seqadmin
      Current Approaches to Protein Sequencing
      by seqadmin


      Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
      04-04-2024, 04:25 PM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 04-11-2024, 12:08 PM
    0 responses
    44 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 10:19 PM
    0 responses
    43 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-10-2024, 09:21 AM
    0 responses
    38 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 04-04-2024, 09:00 AM
    0 responses
    55 views
    0 likes
    Last Post seqadmin  
    Working...
    X