Announcement

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

  • Paired end miRNA-data, trimming question

    Hi!
    I'm new on this forum and in this field of research. I've tried to find an answer to a probably very simple question, but have failed. So I try here and see where I end up?

    I have paired end data for small RNA (sequenced on NextSeq500, library based on NEXTflex Small RNA-kit (Illumina compatible)). I realise most people have single end data, not paired end data for small RNAs. So I have problems with the adapter trimming. It's probably obvious how to handle that, but I'm new and I struggle...

    I've used cutadapt, as recommended by Bioo Scientific, which works fine for R1, but not R2. I obviously don't understand which sequence I need to state for the adapter for the second read. Coz I had expected both R1 and R2 to be processed at more or less the same level since its so short reads? Any help?

    Code:
    cutadapt -a TGGAATTCTCGGGTGCCAAGG -A <adapter???> -o /proj/Trimming/B1_trimmed1_fwd.fq -p /proj/Trimming/B1_trimmed1_rev.fq --minimum-length 23 /proj/FastQ/B1_fwd.fastq.gz /proj/FastQ/B1_rev.fastq.gz
    This is an example from when I ran the script, in this case I had stated the reverse complement for the adapter in the code above. I have also tried what results I get if I trim away the 5' adapter (but with the -g option instead) and don't get any wiser from that. *feel stupid*

    Total read pairs processed: 8,827,513
    Read 1 with adapter: 8,362,506 (94.7%)
    Read 2 with adapter: 45,734 (0.5%)
    Pairs that were too short: 91,927 (1.0%)
    Pairs written (passing filters): 8,735,586 (99.0%)

    Total basepairs processed: 1,341,781,976 bp
    Read 1: 670,890,988 bp
    Read 2: 670,890,988 bp
    Total written (filtered): 950,264,383 bp (70.8%)
    Read 1: 286,493,104 bp
    Read 2: 663,771,279 bp

    I appreciate any help and/or suggestions?
    //├ůsa

  • #2
    You can use the BBMap package to find your adapter sequences like this:

    bbmerge.sh in1=r1.fq.gz in2=r2.fq.gz outa=adapters.fa mininsert=17

    Once you have the adapter sequences, you can trim like this:

    bbduk.sh in1=r1.fq.gz in2=r2.fq.gz out1=trimmed1.fq.gz out2=trimmed2.fq.gz ref=adapters.fa tbo tpe ktrim=r k=21 mink=9 hdist=1

    You should check the contents of adapters.fa to make sure it contains sequence. If you have a mix of adapters it may not work correctly and just output 'N' for the adapter sequences, so in that case you could simply use the included adapters in /bbmap/resources/adapters.fa. But using the specific adapter sequences for your reads is more precise.

    Comment


    • #3
      Originally posted by Brian Bushnell View Post
      You can use the BBMap package to find your adapter sequences like this:

      bbmerge.sh in1=r1.fq.gz in2=r2.fq.gz outa=adapters.fa mininsert=17

      Once you have the adapter sequences, you can trim like this:

      bbduk.sh in1=r1.fq.gz in2=r2.fq.gz out1=trimmed1.fq.gz out2=trimmed2.fq.gz ref=adapters.fa tbo tpe ktrim=r k=21 mink=9 hdist=1

      You should check the contents of adapters.fa to make sure it contains sequence. If you have a mix of adapters it may not work correctly and just output 'N' for the adapter sequences, so in that case you could simply use the included adapters in /bbmap/resources/adapters.fa. But using the specific adapter sequences for your reads is more precise.
      I got a recommendation to use the following script, so I'm running that at the moment (first is the 3' adapter sequence, second is the reverse complement of the 5' adapter (I acctually had not tested that combination before )). Seems to do what I want so far:

      Code:
      cutadapt -a TGGAATTCTCGGGTGCCAAGG -A GATCGTCGGACTGTAGAACTCTGAAC -o R1.trim1.fq -p R2.trim1.fq R1.fastq.gz R2.fastq.gz
      But thank you for your suggestion, I'll try that if this one does not work as I think it does.

      Comment


      • #4
        If using NEXTFlex adapters you also nee to trim the first and last 4 bases from the adapter-clipped reads. Bioo adapters have four random bases at each end, so unless you remove them they will interfere with your mapping.

        Comment


        • #5
          Late addition - but was running into the same issue with paired-end reads using NEXTFlex prep for small RNAs. These posts were really helpful.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Advanced Methods for the Detection of Infectious Disease
            by seqadmin




            The recent pandemic caused worldwide health, economic, and social disruptions with its reverberations still felt today. A key takeaway from this event is the need for accurate and accessible tools for detecting and tracking infectious diseases. Timely identification is essential for early intervention, managing outbreaks, and preventing their spread. This article reviews several valuable tools employed in the detection and surveillance of infectious diseases.
            ...
            11-27-2023, 01:15 PM
          • seqadmin
            Strategies for Investigating the Microbiome
            by seqadmin




            Microbiome research has led to the discovery of important connections to human and environmental health. Sequencing has become a core investigational tool in microbiome research, a subject that we covered during a recent webinar. Our expert speakers shared a number of advancements including improved experimental workflows, research involving transmission dynamics, and invaluable analysis resources. This article recaps their informative presentations, offering insights...
            11-09-2023, 07:02 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, Yesterday, 09:55 AM
          0 responses
          11 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 11-30-2023, 10:48 AM
          0 responses
          17 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 11-29-2023, 08:26 AM
          0 responses
          14 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 11-29-2023, 08:12 AM
          0 responses
          14 views
          0 likes
          Last Post seqadmin  
          Working...
          X