Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • mapping strand-specific RNA-seq reads using TopHat

    I am using ScriptSeq Complete Gold Kit to generate pair-end RNAseq data, and wished to map the reads in a strand-specific manner. I referred TopHat FAQ, which suggested to try 'library-type' - 'fr-firststrand' and 'fr-secondstrand' separately, and the one resulting in many more splicing junctions in junctions.bed should be the right library type.

    I checked the forum, and people said that the kit should go with 'fr-secondstrand', but I still wished to double-check it. So I ran 'unstranded','firststrand' and 'secondstrand' separately on my read pairs, but found the results are quite unexpectedly similar. Specifically the junctions.bed and acceptedhits.bam file are pretty much of the same size for all of them. As I am using pair end, presumably if I force strand-specific mapping, at least one option above should give really poor alignment. But the alignment was quite similar, as you can see from the attached figure (from 1st to 3rd panels are first-strand, second-strand and non-stranded, respectively).

    Can you help see what's wrong here? Thanks!

    Code:
    tophat --library-type fr-secondstrand -o $3 -r $4 -p 4 -g 1 -G $GTF $Ref $read_A $read_B
    Attached Files

  • #2
    Hi,
    Unstranded library of RNASeq can have strand bias, If yes then how can it be ruled out?
    do we really look into it in case of fusion identification using tophat fusion?
    Regards
    Bharati

    Comment


    • #3
      @caswater:
      Tophat by default aligns to the genome and retains discordant mappings as well. So that would explain why your three alignments are fairly similar. A better way to check the strand specificity and library-type is using RSeqC.

      see:
      Bridged amplification & clustering followed by sequencing by synthesis. (Genome Analyzer / HiSeq / MiSeq)

      Techniques and protocol discussions on sample preparation, library generation, methods and ideas

      Comment

      Latest Articles

      Collapse

      • seqadmin
        Quality Control Essentials for Next-Generation Sequencing Workflows
        by seqadmin




        Like all molecular biology applications, next-generation sequencing (NGS) workflows require diligent quality control (QC) measures to ensure accurate and reproducible results. Proper QC begins at nucleic acid extraction and continues all the way through to data analysis. This article outlines the key QC steps in an NGS workflow, along with the commonly used tools and techniques.

        Nucleic Acid Quality Control
        Preparing for NGS starts with isolating the...
        Yesterday, 01:58 PM
      • seqadmin
        An Introduction to the Technologies Transforming Precision Medicine
        by seqadmin


        In recent years, precision medicine has become a major focus for researchers and healthcare professionals. This approach offers personalized treatment and wellness plans by utilizing insights from each person's unique biology and lifestyle to deliver more effective care. Its advancement relies on innovative technologies that enable a deeper understanding of individual variability. In a joint documentary with our colleagues at Biocompare, we examined the foundational principles of precision...
        01-27-2025, 07:46 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by seqadmin, 02-07-2025, 09:30 AM
      0 responses
      19 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-05-2025, 10:34 AM
      0 responses
      37 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 02-03-2025, 09:07 AM
      0 responses
      37 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 01-31-2025, 08:31 AM
      0 responses
      41 views
      0 likes
      Last Post seqadmin  
      Working...
      X