Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Strand specificity issues

    Hi,

    Lately we have been working with bacterial paired-end Illumina HISeq2000 RNAseq data. All the alignments were done in the core facility where we had our RNAseq performed. All I know is that the samples were sequenced in strand-specific manner, rRNA was removed, the length of the reads was 90 bp, and they were aligned using TopHat.

    While analyzing the results I came across some problems. Depending on the way I try to visualize the alignments I keep on getting different results concerning the strand specificity. Therefore I have two main questions:

    Question 1:
    Attached are two files that show the same region viewed with IGV and Artemis. The difference is that Artemis was used to view the original BAM file (Artemis_XS_tagging.png), while IGV was used for viewing processed BAM files, which were first split according to the strand specificity (IGV after spliting.png). BAM files were split using the split.sh script. I’d like to know if the script is a reliable and which results to trust more, the ones seen with Artemis or the ones seen with IGV?

    Question 2:
    Viewing the BAM file with Artemis shows the antisense transcript originating always from the intergenic regions. It appears basically between each end every gene (Artemis_intergenic_regions.png). I assume that this is just an artifact due to technical problems. But where does this problem come from. How can I distinguish the true asRNA from these artifacts?

    * IGV after spliting.png - Snapshot taken from IGV. The BAM file was first split using the script (split.sh) into two separate BAM files based on the strand specificity. Subsequently, both the files were converted into a TDF files. Two replicas of samples from bacteria grown at 22°C (blue bars) and two replicas of samples from bacteria grown at 37°C (red bars) were overlaid.
    * Artemis_XS_tagging.png - Snapshot taken from Artemis. The original BAM file (bacteria grown at 37°C) was viewed with Artemis v16. Reads were colored by the RNAseq strand specific tag (XS). Alignment was viewed as a strand stack.
    * Artemis_intergenic_regions.png - Snapshot taken from Artemis taken in the same way as Artemis_XS_tagging.png.
    Attached Files

  • #2
    Checking your original and processed BAM files using the approach here might help:
    Bridged amplification & clustering followed by sequencing by synthesis. (Genome Analyzer / HiSeq / MiSeq)

    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
    13 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-27-2024, 06:07 PM
    0 responses
    12 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