Ecological and evolutionary dynamics of CRISPR-Cas systems in Clostridium botulinum: Insights from genome mining and comparative analysis

Understanding the prevalence and distribution of CRISPR-Cas systems across different strains can illuminate the ecological and evolutionary dynamics of Clostridium botulinum populations. In this study, we conducted genome mining to characterize the CRISPR-Cas systems of C. botulinum strains. Our analysis involved retrieving complete genome sequences of these strains and assessing the diversity, prevalence, and evolution of their CRISPR-Cas systems. Subsequently, we performed an analysis of homology in spacer sequences from identified CRISPR arrays to investigate and characterize the range of targeted phages and plasmids. Additionally, we investigated the evolutionary trajectory of C. botulinum strains under selective pressures from foreign invasive DNA. Our findings revealed that 306 strains possessed complete CRISPR-Cas structures, comprising 58% of the studied C. botulinum strains. Secondary structure prediction of consensus repeats indicated that subtype II-C, with longer stems compared to subtypes ID and IB, tended to form more stable RNA secondary structures. Moreover, protospacer motif analysis demonstrated that strains with subtype IB CRISPR-Cas systems exhibited 5′-CGG-3′, 5′-CC-3′, and 5′-CAT-3′ motifs in the 3′ flanking regions of protospacers. The diversity observed in CRISPR-Cas systems indicated their classification into subtypes IB, ID, II-C, III-B, and III-D. Furthermore, our results showed that systems with subtype ID and III-D frequently harbored similar spacer patterns. Moreover, analysis of spacer sequences homology with phage and prophage genomes highlighted the specific activities exhibited by subtype IB and III-B against phages and plasmids, providing valuable insights into the functional specialization within these systems.