Question 14 of 20 Fundamentals Practical
When and why should you trim adapters from sequencing reads?
Short answer
Trim when fragments can be shorter than the read length, because the sequencer then reads into the adapter at the 3' end. Leftover adapter bases do not match the genome, so they lower mapping rates and can create false mismatches near read ends.
Why
Every library fragment is flanked by adapters. If the insert is longer than the read, the read never reaches the adapter. If the insert is shorter, the read runs past it into adapter sequence.
Illumina's guidance says adapter sequence can appear at the 3' end of reads when read length exceeds insert size, and lists the sequences to trim for each kit, such as AGATCGGAAGAGC at the start of TruSeq adapters and CTGTCTCTTATACACATCT for many tagmentation kits. Short inserts are common in small RNA libraries, degraded samples and cell-free DNA, so trimming matters most there. Tools such as cutadapt find and remove the adapter, can trim low-quality 3' ends, and can drop reads that become too short.
For two-color instruments such as NextSeq, cutadapt notes that a read falling off the end of its fragment shows a run of high-quality but wrong G calls, which ordinary quality trimming misses, and offers a dedicated option. Trimming is not always needed. Many modern aligners soft-clip unmatched ends, and Illumina's FASTQ generation can trim adapters itself, so check what was already done before trimming twice.
How to work it through
- Check the library kit and expected insert size against the read length.
- Run a QC report and look at adapter content along the read.
- Look up the adapter sequence for the kit from the vendor's list.
- Trim with a tool such as cutadapt, using paired mode for paired reads.
- Discard reads shorter than a minimum length that can still map uniquely.
- Rerun QC to confirm adapter content has gone and record the settings.
An example
A cell-free DNA library with fragments around 160 bp is sequenced with 2 x 150 bp reads. The QC report shows adapter content rising from 20 percent of reads beyond position 120. After trimming the TruSeq adapter in paired mode and discarding reads shorter than 30 bases, the mapping rate rises from 88 to 97 percent and a cluster of false mismatches near read ends disappears.
Where the answer stops holding
Aggressive quality trimming can shorten reads enough to make them map less uniquely. Trimming with the wrong adapter sequence removes nothing useful. Some pipelines, such as those relying on fixed read lengths, can be affected by variable read lengths after trimming.
A strong answer covers
- Adapter read-through happens when inserts are shorter than reads
- Use the vendor's adapter sequence for the kit
- Trim paired reads together and set a minimum length
- Two-color chemistry needs special handling of poly-G ends
- Check whether trimming was already done upstream
Common mistakes
- Trimming single files of a pair separately so the pairs fall out of order
- Assuming a high mapping rate means adapters are absent
- Applying a fixed number of trimmed bases instead of finding the adapter
The interviewer might ask next
How would you estimate the insert size distribution from aligned paired reads?