Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • rRNA and adapter removal bondaries in RNA-seq HiSeq 3000 Illumina

    Hello

    I am new in RNA-seq analysis and it is difficult for me to find out info about quality assessment in the matter. I did an analysis with fastqc and MultiQC in a Paired End, 125nts (2x125), HiSeq3000 experiment. I obtained these two plots (attached) and I am not sure if a rRNA or adapter removal is required prior to alignment.

    In the case of the adapter content plot it looks like all samples have adapters included, in fact none of the samples passed the QA. In the case of the plot with %GC content (I read it can be used to detect a possible rRNA enrichment) I see a proportion of reads having a GC content higher than 60%, but I am not sure if it is enough to confirm a high presence of rRNA that could affect to posterior analysis.

    Any help will be useful!

    Thanks

    Magda
    Attached Files

  • #2
    You should scan/trim your sequences with a suitable program (bbduk.sh from BBMap, trimmomatic, cutadapt etc.). While some rely on aligners to soft-clip adapter sequences (and they likely will), trimming them beforehand ensures that your initial data is clean.

    Presence of rRNA is a separate issue. If you feel that your data has rRNA reads in it then the best solution is to align a sample (or whole data) to rDNA repeat for your organism. While rRNA reads would not directly interfere with counts if you have uneven levels of them across samples then they can indirectly affect your results.

    Comment


    • #3
      Thanks for your quick answer, it was very useful! Just a last question, could you recommend me a tool or database I could use to make the alignment with the rDNA repeats of an organism?

      Thanks

      Comment


      • #4
        What organism are you working with?

        Comment


        • #5
          We are working with human

          Comment


          • #6
            You can get the human rDNA repeat from here. If you do see contamination you can use bbsplit.sh from BBMap suite to separate those reads from rest of your data, if you wish.

            Comment

            Latest Articles

            Collapse

            • 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
            • seqadmin
              The Impact of AI in Genomic Medicine
              by seqadmin



              Artificial intelligence (AI) has evolved from a futuristic vision to a mainstream technology, highlighted by the introduction of tools like OpenAI's ChatGPT and Google's Gemini. In recent years, AI has become increasingly integrated into the field of genomics. This integration has enabled new scientific discoveries while simultaneously raising important ethical questions1. Interviews with two researchers at the center of this intersection provide insightful perspectives into...
              02-26-2024, 02:07 PM

            ad_right_rmr

            Collapse

            News

            Collapse

            Topics Statistics Last Post
            Started by seqadmin, 03-14-2024, 06:13 AM
            0 responses
            32 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-08-2024, 08:03 AM
            0 responses
            71 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-07-2024, 08:13 AM
            0 responses
            80 views
            0 likes
            Last Post seqadmin  
            Started by seqadmin, 03-06-2024, 09:51 AM
            0 responses
            68 views
            0 likes
            Last Post seqadmin  
            Working...
            X