Exploiting immune cell receptor information to quantify index switching in single cell transcriptome sequencing experiment

By offering high sequencing speed and ultra-high-throughput at a low price, Illumina next-generation sequencing platforms has been widely adopted in recent years. However, an experiment with multiplexed library could be at risk of molecular recombination, known as “index switching”, which causes certain amount of reads to be assigned to an incorrect library. It is reported that a new advance, exclusion amplification (ExAmp) in conjunction with the patterned flow cell technology introduced on HiSeq 3000/HiSeq 4000/HiSeq X sequencing systems, potentially suffers from a higher rate of index switching than conventional bridge amplification. We took advantage of the diverse but highly cell-unique expression of immune cell receptors to quantify index switching on single cell RNA-seq data that were sequenced on HiSeq 3500 and HiSeq 4000. By utilizing the unique T-cell receptor (TCR) expression, we could quantify the spread-of-signal from many different wells (n=51 from total three batches) due to index switching. We used TraCer to reconstruct full-length TCR from all samples, and then used Kallisto to quantify TCR gene expression. We found index switching in all three batches of experiments investigated. The median percentage of incorrectly detected markers was estimated to be 4.2% (interquartile range (IQR): 2.0%-8.7%). We did not detect any consistent pattern of certain indices to be more prone for switching than others, suggesting that index switching is a stochastic process. We confirm that index switching is a problem that affects all samples run in multiplexed libraries on Illumina HiSeq 3500 and HiSeq 4000 platforms.

• Type: Other
• Archiver: European Genome-Phenome Archive (EGA)