Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Low mapping TopHat

    Hi guys,

    another problem to solve. I am running TopHat on some human cells samples I have submitted to our facility for RNAseq. The bioanalyzer values were great (all above RIN 8.5) and thus a library based exclusively on ribo depletion has been created.
    I have multiplexed 4 samples in a lane (75bp, single-end) and obtained roughly 50 mil reads for each sample.
    This is my TopHat summary:

    Reads kept by tophat2 for mapping: 48777778
    Reads discarded by tophat2: 321249

    Number of unique reads that mapped: 5922988
    Percentage of unique reads that mapped: 12.06%

    I have used transcriptome for Tophat2...but I cannot understand why such a low mapping percentage (not the same across the samples, ranges across 12-40%).
    Do you have any idea???
    Other groups in our lab (doing polyA selection, on same cells, using same RNAseq pipeline) map almost 90% of the reads.

    Help me please!!!
    Thanks
    Manu

  • #2
    Perhaps a lot of rRNA left in the samples, leading to many multimapping reads due to repetitiveness, or "duplicates" which are really rRNA genes with very high coverage. It could also be adapter contamination but I have seen such low mapping rates often in e g RiboMinus samples due to poor rRNA removal.

    Comment


    • #3
      It might be useful to grab some random read sequences that didn't align (tophat makes a file called unmapped.bam) and try aligning them with BLAST and see what it comes up with. not very efficient but if you repeatedly see poor / no alignments from BLAST you can be pretty sure that Tophat isn't messing up.

      You can pull some reads from it like this:

      Code:
      samtools view unmapped.bam | \
        perl -slane 'if(rand() < 0.000001) { print ">$F[0]\n$F[9]"; }' > random_reads.fa
      You can tweak that value in the if test to control how much gets extracted. What I did here should pull on average 1 read per million.
      /* Shawn Driscoll, Gene Expression Laboratory, Pfaff
      Salk Institute for Biological Studies, La Jolla, CA, USA */

      Comment


      • #4
        Thanks for both suggestions!
        I have run the script and I BLASTed the unmapped reads of one of my samples...
        apparently most of my reads match with a pre-rRNA (or rRNA in general) and also incredibly to some exotic fish, which I assume is a mere coincidence and an artifact.
        So my question is: shouldn't rRNA map anyhow to the transcriptome?
        More: is there any strategy to limit such a huge contamination? (I really need to do ribo depletion rather than polyA selection)

        Thanks a lot!!!
        Manu

        Comment


        • #5
          Hi,

          I have a quick follow up question......I am using that perl script to see why i am getting sooo low % of mapped reads (60%). When i run it i only pull ~ 11 samples of ~ 120 bp (not really a million reads). Am i missing something?

          Comment


          • #6
            Originally posted by Gonza View Post
            Hi,

            I have a quick follow up question......I am using that perl script to see why i am getting sooo low % of mapped reads (60%). When i run it i only pull ~ 11 samples of ~ 120 bp (not really a million reads). Am i missing something?
            @sdriscoll's code pulls one read from a million sampled so you must have about 11 million reads that are in unmapped.bam. You can change the "0.000001" value to sample more reads.

            While it is possible that your samples are contaminated with something (either because of something you did or the sequencing facility did) you will want to take a more careful look before concluding anything.

            Brian asked a set of questions in your other post so let us see what you have to say about those.

            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, 03-27-2024, 06:37 PM
            0 responses
            12 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-27-2024, 06:07 PM
            0 responses
            11 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-22-2024, 10:03 AM
            0 responses
            53 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-21-2024, 07:32 AM
            0 responses
            69 views
            0 likes
            Last Post seqadmin  
            Working...
            X