Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Bowtie2: rRNA data from Silva fail to align with Illumina PE reads

    Hi!

    Could you help me with my long list of problems? I have Illumina HiSeq paired-end RNA data from an unmodel organism. The reads have a lot of rRNA and other sources of contamination and they have quite low quality scores.

    I'm trying to filter out the rRNA out of the samples using the whole length shotgun sequences from a closest model species. I can't use the species whose rRNA genes we have found from the data because they have only partial sequences available. So I downloaded the rRNA SSU and LSU r117 sequences from Silva as FASTA without gaps and decompressed the .tar.gz files into FASTA before usage. Should I use some other format than FASTA without gaps?

    I created the index out of the Silva-files and aligned the reads using Bowtie2 using simplest possible commands:

    bowtie2-build C_brenneriSSU.fasta C_brenneriSSU_index
    bowtie2-align -x C_brenneriSSU_index -1 ove2_R1.fastq -2 ove2_R2.fastq -S result.fastq
    But nothing is aligning. Why is that? The reads are untrimmed because I've read that it's best for Bowtie2 to use the original reads.

    We have also tried to align some example rRNA-sequencefiles and the Silva-files with consensuses of contigs constructed from combined samples gained from CLC assembly but they won't align with the rRNA datas with Bowtie2 either. Any ideas why not?

    And most importantly, what should I do to find the rRNA sequences? I'm familiar with extracting the unmapped reads by using samtools so it's basically just the alignment that seems to be the problem.
    Last edited by jyuems; 05-07-2014, 04:12 AM.

  • #2
    This is the result from Bowtie2:

    36906390 pairs aligned 0 times concordantly or discordantly; of these:
    73812780 mates make up the pairs; of these:
    73812780 (100.00%) aligned 0 times
    0 (0.00%) aligned exactly 1 time
    0 (0.00%) aligned >1 times
    What could be going wrong here?

    Comment


    • #3
      Hi Jyumes...

      I am also new but all i can suggest that use MGRAST for rRNA Annotation. There is a paper in plos for this study.
      Also can download rRNA database for prokaryotes and eukaryote and perform Blast (all vs all) for tabular out put.
      this will solve the rRNA reads.

      Also i would like to ask have you used any rRNA cancellation method during RNA isolation otherwise there is uptp70% rRNA reads contamination in your sequences.

      Comment


      • #4
        Thank you guptavipin142! I will give the MGRAST a go.

        Yes, the samples were attempted to be purified from rRNA twice during isolation but according to BLAST results there are still a many rRNA sequences. Also according to FastQC report there are a lot of overrepresented sequences present in the samples.

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