In Silico Processing of the Complete CRISPR-Cas Spacer Space for Identification of PAM Sequences

Abstract

Despite extensive exploration of the diversity of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) systems, biological applications have been mostly confined to Class 2 systems, specifically the Cas9 and Cas12 (formerly Cpf1) single effector proteins. A key limitation of exploring and utilizing other CRISPR-Cas systems with unique functionalities, particularly Class I types and their multi-protein effector complex, is the knowledge of the systems protospacer adjacent motif (PAM) sequence identity. In this work, we developed a systematic pipeline, named CASPERpam, that enables us to comprehensively assess the PAM sequences of all the available CRISPR-Cas systems in the NCBI database of bacterial genomes. The CASPERpam analysis revealed that within the 30,389 assemblies previously screen for CRISPR arrays, there exists 26,364 spacers that match somewhere in the viral, bacterial, and plasmid databases of NCBI, using the constraints of 95% sequence identity and 95% sequence coverage for blast hits. When grouping these results by species, we were able to identify putative PAM sequences for 1,049 among 1,493 unique species. The remaining species either have insufficient data or an undetermined result from the analysis. Finally, we were able to infer certain design principles that are relevant for understanding PAM diversity and a baseline for further experimental studies including PAM assays. We envision CASPERpam is a useful bioinformatic tool for understanding and harnessing the diversity of CRISPR systems.