2
$\begingroup$

I was going over different sequencing library protocols and noticed that the adapter sequence can vary between protocols. In some adapter sequences the index sits between the primer binding site and the insert. In this case during sequencing the index will be sequenced and then the insert that comes right after it. In other cases, the adapter sequence is built so that the index is outside the primer binding site, i.e. index - primer binding site - insert (see two leftmost parts of the attached figure). Here the insert is sequenced using the Read1 primer and afterwards the index is sequenced in a separate reaction using the 17 index primer. figure 1

In this case are the index reads just added on to the ends of the insert reads in the FastQ files? Is there a way to tell from the FastQ files where the index was positioned within the adapter?

$\endgroup$
  • $\begingroup$ Software FastQC reports overrepresented kmers per position in read, perhaps that could help you figure out the order? Also, note that some sequencing providers are doing some of the trimming before they hand users the data. $\endgroup$ – Kamil S Jaron Nov 2 at 10:46
  • 1
    $\begingroup$ I do think this is a bioinformatics question. Ideally you could ignore this as an implementation detail that doesn't concern bioinformatics, but in reality I've had to stare at diagrams exactly like this when trying to troubleshoot strange illumina runs and interpret changes in base quality, for example. I'd like to see a clear answer to this myself. $\endgroup$ – Jesse Nov 2 at 16:35
0
$\begingroup$

This is definitely a bioinformatics question! Understanding how sequencing libraries are created in the wet-lab and how they are sequenced is fundamental to bioinformatics.

To answer the question, the index sequences are not included in the fastq file output from Illumina, or other machines. Illumina sequencers just use the index sequence to demultiplex the complex mixture of molecules in the sequencing pool. If a small "dot" on the glass sequencing slide has some sequence, but the indices don't match what the technician specified for a library, then that read is not included in the fastq output of that library since it could be contamination of some sort.

The sequences in the fasta file theoretically all should be the insert from the library that you prepared, but sometimes the molecule being sequenced is short, and there is read-through of the adapter sequence. This is the only circumstance under which the i7 or i5 index sequences might be in the normal fastq output sequence line. This is not ideal, and this is why people perform adapter trimming.

| improve this answer | |
$\endgroup$
  • 3
    $\begingroup$ The index sequences can be returned as separate fastqs, but most of the time, you don't need them. And cluster whose index doesn't match what's on the sample sheet does get outputted, in the "Undetermined" fastq. $\endgroup$ – swbarnes2 Nov 3 at 20:16
0
$\begingroup$

The order of sequencing in a dual index flow cell is:

  1. Read 1
  2. Index 1 (I7)
  3. Index 2 (I5)
  4. Read 2

Between each step there is a switch of the primer used (Read1-primer -> Index1-primer (I7) -> Index2-primer (I5) -> Read2-primer)
There are differences in how the I5 is dealt with between dual- and single-read flow cells as explained in the Illumina Indexed Sequencing Overview. In this document there is a lot more detailed explanation of the sequencing workflow that may be helpful to you.

But as explained by the other answers, none of this will matter in most cases. Even if you decide to output the Index sequences they will be exported as separate FASTQ files.

| improve this answer | |
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.