A Comprehensive Curation Shows the Dynamic Evolutionary Patterns of Prokaryotic CRISPRs

[Why so rare if so essentiel: the determinants of the sparse distribution of CRISPR-Cas systems in bacterial genomes]

CRISPR-Cas (Cluster of Regularly Interspaced Short Palindromic Repeats) systems confer bacteria and archaea an adaptative immunity against phages and other invading genetic elements playing an important role in bacterial evolution. However, despite the protection they generate and high rate of horizontal transfer, less than 50% of bacterial genomes harbor a CRISPR-Cas system. As a comparison, 90% of archaea encode a CRISPR-Cas system and a bacterial genome codes for two restriction modification systems on average. This review describes CRISPR-Cas systems distribution in bacterial genomes and then details the different hypotheses put forward to explain the relative scarcity of CRISPR-Cas systems. More specifically, phage escape mechanisms, ecological factors such as phage diversity and abundance and intrinsic costs, such as maintenance or autoimmunity, are discussed. Overall, a better understanding of the downsides of encoding CRISPR-Cas systems is essential to explain their evolutionary dynamics and their relative success in different environments and clades.

CRISPRdigger: detecting CRISPRs with better direct repeat annotations

Clustered regularly interspaced short palindromic repeats (CRISPRs) are important genetic elements in many bacterial and archaeal genomes, and play a key role in prokaryote immune systems’ fight against invasive foreign elements. The CRISPR system has also been engineered to facilitate target gene editing in eukaryotic genomes. Using the common features of mis-annotated CRISPRs in prokaryotic genomes, this study proposed an accurate de novo CRISPR annotation program CRISPRdigger, which can take a partially assembled genome as its input. A comprehensive comparison with the three existing programs demonstrated that CRISPRdigger can recover more Direct Repeats (DRs) for CRISPRs and achieve a higher accuracy for a query genome. The program was implemented by Perl and all the parameters had default values, so that a user could annotate CRISPRs in a query genome by supplying only a genome sequence in the FASTA format. All the supplementary data are available at http://www.healthinformaticslab.org/supp/.