AcaFinder: genome mining for anti-CRISPR associated genes

Abstract Anti-CRISPR (Acr) proteins are encoded by (pro)viruses to inhibit their host’s CRISPR-Cas systems. Genes encoding Acr and Aca (Acr associated) proteins often co-localize to form acr-aca operons. Here, we present AcaFinder as the first Aca genome mining tool. AcaFinder can: (i) predict Acas and their associated acr-aca operons using guilt-by-association (GBA); (ii) identify homologs of known Acas using an HMM (Hidden Markov model) database; (iii) take input genomes for potential prophages, CRISPR-Cas systems, and self-targeting spacers (STSs); and (iv) provide a standalone program ( https://github.com/boweny920/AcaFinder ) and a web server ( http://aca.unl.edu/Aca ). AcaFinder was applied to mining over 16,000 prokaryotic and 142,000 gut phage genomes. After a multi-step filtering, 36 high-confident new Aca families were identified, which is three times of the 12 known Aca families. Seven new Aca families were from major human gut bacteria (Bacteroidota, Actinobacteria, Fusobacteria) and their phages, while most known Aca families were from Proteobacteria and Firmicutes. A complex association network between Acrs and Acas was revealed by analyzing their operonic co-localizations. It appears very common in evolution that the same aca genes can recombine with different acr genes and vice versa to form diverse acr-aca operon combinations. Importance At least four bioinformatics programs have been published for genome mining of Acrs since 2020. In contrast, no bioinformatics tools are available for automated Aca discovery. As the self-transcriptional repressor of acr-aca operons, Aca can be viewed as anti-anti-CRISPRs, with a great potential in the improvement of CRISPR-Cas technology. Although all the 12 known Aca proteins contain a conserved Helix-Turn-Helix (HTH) domain, not all HTH-containing proteins are Acas. However, HTH-containing proteins with an adjacent Acr homologs encoded in the same genetic operon are likely Aca proteins. AcaFinder implements this guilt-by-association (GBA) idea and the idea of using HMMs of known Acas for homologs into one software package. Applying AcaFinder in screening prokaryotic and gut phage genomes reveals a complex acr-aca operonic co-localization network between different families of Acrs and Acas.