bcSeq: an R package for fast sequence mapping in high-throughput shRNA and CRISPR screens

Area of Science:

• Genomics
• Bioinformatics
• Molecular Biology

Background:

• CRISPR-Cas9 and shRNA screens are powerful tools in biological research.
• Accurate analysis of high-throughput sequencing data is crucial for reliable results.
• Existing methods may struggle with sequencing errors and ambiguous barcode mappings.

Purpose of the Study:

• To develop a computationally efficient and statistically robust method for analyzing CRISPR-Cas9 and shRNA high-throughput sequencing screens.
• To introduce bcSeq, an open-source R package designed for this purpose.

Main Methods:

• Implementation of a fast and parallelized algorithm for mapping sequencing reads to barcode libraries.
• Utilization of a Trie data structure for enhanced computational speed.
• Application of a statistical sequencing error model based on Phred scores to resolve ambiguous mappings.

Main Results:

• bcSeq provides a statistically principled approach to handle sequencing errors and ambiguous mappings.
• The package offers a computationally efficient solution for analyzing large-scale sequencing screen data.
• The algorithm demonstrates speed and accuracy in mapping reads to barcode libraries.

Conclusions:

• bcSeq is a valuable open-source tool for researchers utilizing CRISPR-Cas9 and shRNA high-throughput sequencing screens.
• The package improves the accuracy and efficiency of data analysis in basic and translational research.
• bcSeq facilitates more reliable interpretation of screening experiments.