CRISPRCasTyper: Automated Identification, Annotation, and Classification of CRISPR-Cas Loci

Complete genomes of Asgard archaea reveal diverse integrated and mobile genetic elements

Asgard archaea are of great interest as the progenitors of Eukaryotes, but little is known about the mobile genetic elements (MGEs) that may shape their ongoing evolution. Here, we describe MGEs that replicate in Atabeyarchaeia, a wetland Asgard archaea lineage represented by two complete genomes. We used soil depth-resolved population metagenomic data sets to track 18 MGEs for which genome structures were defined and precise chromosome integration sites could be identified for confident host linkage. Additionally, we identified a complete 20.67 kbp circular plasmid and two family-level groups of viruses linked to Atabeyarchaeia, via CRISPR spacer targeting. Closely related 40 kbp viruses possess a hypervariable genomic region encoding combinations of specific genes for small cysteine-rich proteins structurally similar to restriction-homing endonucleases. One 10.9 kbp integrative conjugative element (ICE) integrates genomically into the Atabeyarchaeum deiterrae-1 chromosome and has a 2.5 kbp circularizable element integrated within it. The 10.9 kbp ICE encodes an expressed Type IIG restriction-modification system with a sequence specificity matching an active methylation motif identified by Pacific Biosciences (PacBio) high-accuracy long-read (HiFi) metagenomic sequencing. Restriction-modification of Atabeyarchaeia differs from that of another coexisting Asgard archaea, Freyarchaeia, which has few identified MGEs but possesses diverse defense mechanisms, including DISARM and Hachiman, not found in Atabeyarchaeia. Overall, defense systems and methylation mechanisms of Asgard archaea likely modulate their interactions with MGEs, and integration/excision and copy number variation of MGEs in turn enable host genetic versatility.

Crosstalk between three CRISPR-Cas types enables primed type VI-A adaptation in Listeria seeligeri

CRISPR-Cas systems confer adaptive immunity to their prokaryotic hosts through the process of adaptation, where sequences are captured from foreign nucleic acids and integrated as spacers in the CRISPR array, and thereby enable crRNA-guided interference against new threats. While the Cas1-2 integrase is critical for adaptation, it is absent from many CRISPR-Cas loci, rendering the mechanism of spacer acquisition unclear for these systems. Here we show that the RNA-targeting type VI-A CRISPR system of Listeria seeligeri acquires spacers from DNA substrates through the action of a promiscuous Cas1-2 integrase encoded by a co-occurring type II-C system, in a transcription-independent manner. We further demonstrate that the type II-C integration complex is strongly stimulated by preexisting spacers in a third CRISPR system (type I-B) which imperfectly match phage targets and prime type VI-A adaptation. Altogether, our results reveal an unprecedented degree of communication among CRISPR-Cas loci encoded by a single organism.

The CRISPR-associated adenosine deaminase Cad1 converts ATP to ITP to provide antiviral immunity

Type III CRISPR systems provide immunity against genetic invaders through the production of cyclic oligo-adenylate (cAn) molecules that activate effector proteins that contain CRISPR-associated Rossman fold (CARF) domains. Here, we characterized the function and structure of an effector in which the CARF domain is fused to an adenosine deaminase domain, CRISPR-associated adenosine deaminase 1 (Cad1). We show that upon binding of cA4 or cA6 to its CARF domain, Cad1 converts ATP to ITP, both in vivo and in vitro. Cryoelectron microscopy (cryo-EM) structural studies on full-length Cad1 reveal an hexameric assembly composed of a trimer of dimers, with bound ATP at inter-domain sites required for activity and ATP/ITP within deaminase active sites. Upon synthesis of cAn during phage infection, Cad1 activation leads to a growth arrest of the host that prevents viral propagation. Our findings reveal that CRISPR-Cas systems employ a wide range of molecular mechanisms beyond nucleic acid degradation to provide adaptive immunity in prokaryotes.

Metagenomic characterization of viruses and mobile genetic elements associated with the DPANN archaeal superphylum

The archaeal superphylum DPANN (an acronym formed from the initials of the first five phyla discovered: Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanohaloarchaeota and Nanoarchaeota) is a group of ultrasmall symbionts able to survive in extreme ecosystems. The diversity and dynamics between DPANN archaea and their virome remain largely unknown. Here we use a metagenomic clustered regularly interspaced short palindromic repeats (CRISPR) screening approach to identify 97 globally distributed, non-redundant viruses and unclassified mobile genetic elements predicted to infect hosts across 8 DPANN phyla, including 7 viral groups not previously characterized. Genomic analysis suggests a diversity of viral morphologies including head-tailed, tailless icosahedral and spindle-shaped viruses with the potential to establish lytic, chronic or lysogenic infections. We also find evidence of a virally encoded Cas12f1 protein (probably originating from uncultured DPANN archaea) and a mini-CRISPR array, which could play a role in modulating host metabolism. Many metagenomes have virus-to-host ratios >10, indicating that DPANN viruses play an important role in controlling host populations. Overall, our study illuminates the underexplored diversity, functional repertoires and host interactions of the DPANN virome.

PlasmidScope: a comprehensive plasmid database with rich annotations and online analytical tools

Plasmids are extrachromosomal genetic molecules that replicate independent of chromosomes in bacteria, archaea, and eukaryotic organisms. They contain diverse functional elements and are capable of horizontal gene transfer among hosts. While existing plasmid databases have archived plasmid sequences isolated from individual microorganisms or natural environments, there is a need for a comprehensive, standardized, and annotated plasmid database to address the vast accumulation of plasmid sequences. Here, we propose PlasmidScope (https://plasmid.deepomics.org/), a plasmid database offering comprehensive annotations, automated online analysis, and interactive visualization. PlasmidScope harbors a substantial collection of 852 600 plasmids curated from 10 repositories. Along with consolidated background information, PlasmidScope utilizes 12 state-of-the-art tools and provides comprehensive annotations for the curated plasmids, covering genome completeness, topological structure, mobility, host source, tRNA, tmRNA, signal peptides, transmembrane proteins and CRISPR/Cas systems. PlasmidScope offers diverse functional annotations for its 25 231 059 predicted genes from 9 databases as well as corresponding protein structures predicted by ESMFold. In addition, PlasmidScope integrates online analytical modules and interactive visualization, empowering researchers to delve into the complexities of plasmids.

Metagenomic insights into microalgae-bacterium-virus interactions and viral functions in phycosphere facing environmental fluctuations

Despite the ecological and biotechnological significance of microalgae-bacterium symbionts, the response of host-virus interactions to external environmental fluctuations and the role of viruses in phycosphere remain largely unexplored. Herein, we employed algal-bacterial granular sludge (ABGS) with varying light intensity and organic carbon loading to investigate the mechanisms of microalgae-bacterium-virus symbionts in response to environmental fluctuations. Metagenomics revealed that enhanced light intensity decreased the diversity of microalgae, so did the diversity of symbiotic bacteria and viruses. As carbon sources decreased, bacteria prompted horizontal gene transfer in phycosphere by 12.76 %-157.40 %, increased the proportion of oligotrophs as keystone species (0.00 % vs 14.29 %) as well as viruses using oligotrophs as hosts (18.52 % vs 25.00 %). Furthermore, virus-carried auxiliary metabolic genes (AMGs) and biosynthetic gene clusters (BGCs) encoding vitamin B12 synthesis (e.g., cobS), antioxidation (e.g., queC), and microbial aggregation (e.g., cysE). Additionally, phylogenetic and similarity analysis further revealed the evolutionary origin and potential horizontal transfer of the AMGs and BGCs, which could potentially enhance the adaptability of bacteria and eukaryotic microalgae. Overall, our research demonstrates that environmental fluctuations have cascading effects on the microalgae-bacteria-virus interactions, and emphasizes the important role of viruses in maintaining the stability of the phycosphere symbiotic community.

Expanded genome and proteome reallocation in a novel, robust Bacillus coagulans strain capable of utilizing pentose and hexose sugars

Bacillus coagulans, a Gram-positive thermophilic bacterium, is recognized for its probiotic properties and recent development as a microbial cell factory. Despite its importance for biotechnological applications, the current understanding of B. coagulans' robustness is limited, especially for undomesticated strains. To fill this knowledge gap, we characterized the metabolic capability and performed functional genomics and systems analysis of a novel, robust strain, B. coagulans B-768. Genome sequencing revealed that B-768 has the largest B. coagulans genome known to date (3.94 Mbp), about 0.63 Mbp larger than the average genome of sequenced B. coagulans strains, with expanded carbohydrate metabolism and mobilome. Functional genomics identified a well-equipped genetic portfolio for utilizing a wide range of C5 (xylose, arabinose), C6 (glucose, mannose, galactose), and C12 (cellobiose) sugars present in biomass hydrolysates, which was validated experimentally. For growth on individual xylose and glucose, the dominant sugars in biomass hydrolysates, B-768 exhibited distinct phenotypes and proteome profiles. Faster growth and glucose uptake rates resulted in lactate overflow metabolism, which makes B. coagulans a lactate overproducer; however, slower growth and xylose uptake diminished overflow metabolism due to the high energy demand for sugar assimilation. Carbohydrate Transport and Metabolism (COG-G), Translation (COG-J), and Energy Conversion and Production (COG-C) made up 60%-65% of the measured proteomes but were allocated differently when growing on xylose and glucose. The trade-off in proteome reallocation, with high investment in COG-C over COG-G, explains the xylose growth phenotype with significant upregulation of xylose metabolism, pyruvate metabolism, and tricarboxylic acid (TCA) cycle. Strain B-768 tolerates and effectively utilizes inhibitory biomass hydrolysates containing mixed sugars and exhibits hierarchical sugar utilization with glucose as the preferential substrate.IMPORTANCEThe robustness of B. coagulans makes it a valuable microorganism for biotechnology applications; yet, this phenotype is not well understood at the cellular level. Through phenotypic characterization and systems analysis, this study elucidates the functional genomics and robustness of a novel, undomesticated strain, B. coagulans B-768, capable of utilizing inhibitory switchgrass biomass hydrolysates. The genome of B-768, enriched with carbohydrate metabolism genes, demonstrates high regulatory capacity. The coordination of proteome reallocation in Carbohydrate Transport and Metabolism (COG-G), Translation (COG-J), and Energy Conversion and Production (COG-C) is critical for effective cell growth, sugar utilization, and lactate production via overflow metabolism. Overall, B-768 is a novel, robust, and promising B. coagulans strain that can be harnessed as a microbial biomanufacturing platform to produce chemicals and fuels from biomass hydrolysates.

Differentiation of the Anammox core microbiome: Unraveling the evolutionary impetus of scalable gene flow

Anaerobic ammonium oxidation bacteria (AAOB), distinguished by their unique autotrophic nitrogen metabolism, hold pivotal positions in the global nitrogen cycle and environmental biotechnologies. However, the ecophysiology and evolution of AAOB remain poorly understood, attributed to the absence of monocultures. Hence, a comprehensive elucidation of the AAOB-dominated core microbiome, anammox core, is imperative to further completing the theory of engineered nitrogen removal and ecological roles of anammox. Performing taxonomic and phylogenetic analyses on collected genome repertoires, we show here that Candidatus Brocadia and Candidatus Kuenenia possesses a more compact core than Candidatus Jettenia, which partly explains why the latter has a less common ecological presence. Evidence of gene flow is particularly striking in functions related to biosynthesis and oxygen detoxification, underscoring the evolutionary forces driving lineage and core differentiation. Furthermore, CRISPR spacer traceback of the AAOB metagenome-assembled genomes (MAGs) reveals a series of genetic traces for the concealed phages. By reconceptualizing the functional divergence of AAOB with the historical role of phages, we ultimately propose a coevolutionary framework to understand the evolutionary trajectory of anammox microecology. The discoveries provided in this study offer new insights into understanding the evolution of AAOB and the ecology of anammox.

Targeted phage hunting to specific Klebsiella pneumoniae clinical isolates is an efficient antibiotic resistance and infection control strategy

Klebsiella pneumoniae is one of the most threatening multi-drug-resistant pathogens today, with phage therapy being a promising alternative for personalized treatments. However, the intrinsic capsule diversity in Klebsiella spp. poses a substantial barrier to the phage host range, complicating the development of broad-spectrum phage-based treatments. Here, we have isolated and genomically characterized phages capable of infecting each of the acquired 77 reference serotypes of Klebsiella spp., including capsular types widespread among high-risk K. pneumoniae clones causing nosocomial infections. We demonstrated the possibility of isolating phages for all capsular types in the collection, revealing high capsular specificity among taxonomically related phages, in contrast to a few phages that exhibited broad-spectrum infection capabilities. To decipher the determinants of the specificity of these phages, we focused on their receptor-binding proteins, with particular attention to depolymerases. We also explored the possibility of designing a broad-spectrum phage cocktail based on phages isolated in reference capsular-type strains and determining the ability to lyse relevant clinical isolates. A combination of 12 phages capable of infecting 55% of the reference Klebsiella spp. serotypes was tested on a panel of carbapenem-resistant K. pneumoniae clinical isolates. Thirty-one percent of isolates were susceptible to the phage cocktail. However, our results suggest that in a highly variable encapsulated bacterial host, phage hunting must be directed to the specific Klebsiella isolates. This work is a step forward in the understanding of the complexity of phage-host interactions and highlights the importance of implementing precise and phage-specific strategies to treat K. pneumoniae infections worldwide.IMPORTANCEThe emergence of resistant bacteria is a serious global health problem. In the absence of effective treatments, phages are a personalized and effective therapeutic alternative. However, little is still known about phage-host interactions, which are key to implementing effective strategies. Here, we focus on the study of Klebsiella pneumoniae, a highly pathogenic encapsulated bacterium. The complexity and variability of the capsule, where in most cases phage receptors are found, make it difficult for phage-based treatments. Here, we isolated a large collection of Klebsiella phages against all the reference strains and in a cohort of clinical isolates. Our results suggest that clinical isolates represent a challenge, especially high-risk clones. Thus, we propose targeted phage hunting as an effective strategy to implement phage-derived therapies. Our results are a step forward for new phage-based strategies to control K. pneumoniae infections, highlighting the importance of understanding phage-host interactions to design personalized treatments against Klebsiella spp.

Diverse viral cas genes antagonize CRISPR immunity

Prokaryotic CRISPR-Cas immunity is subverted by anti-CRISPRs (Acrs), which inhibit Cas protein activities when expressed during the phage lytic cycle or from resident prophages or plasmids1. Acrs often bind to specific cognate Cas proteins, and hence inhibition is typically limited to a single CRISPR-Cas subtype2. Furthermore, although acr genes are frequently organized together in phage-associated gene clusters3, how such inhibitors initially evolve has remained unclear. Here we investigated the Acr content and inhibition specificity of diverse Listeria isolates, which naturally harbour four CRISPR-Cas systems (types I-B, II-A, II-C and VI-A). We observed widespread antagonism of CRISPR, which we traced to 11 previously unknown and 4 known acr gene families encoded by endogenous mobile elements. Among these were two Acrs that possess sequence homology to type I-B Cas proteins, one of which assembles into a defective interference complex. Surprisingly, an additional type I-B Cas homologue did not affect type I immunity, but instead inhibited the RNA-targeting type VI CRISPR system by means of CRISPR RNA (crRNA) degradation. By probing viral sequence databases, we detected abundant orphan cas genes located within putative anti-defence gene clusters. Among them, we verified the activity of a particularly broad-spectrum cas3 homologue that inhibits type I-B, II-A and VI-A CRISPR immunity. Our observations provide direct evidence of Acr evolution by cas gene co-option, and new genes with potential for broad-spectrum control of genome editing technologies.

Methanotrophic Methanoperedens archaea host diverse and interacting extrachromosomal elements

Methane emissions are mitigated by anaerobic methane-oxidizing archaea, including Methanoperedens. Some Methanoperedens host huge extrachromosomal genetic elements (ECEs) called Borgs that may modulate their activity, yet the broader diversity of Methanoperedens ECEs is understudied. Here we report small enigmatic linear ECEs, circular viruses and unclassified ECEs that are predicted to replicate within Methanoperedens. Linear ECEs have inverted terminal repeats, tandem repeats and coding patterns that are strongly reminiscent of Borgs, but they are only 52-145 kb in length. As they share proteins with Borgs and Methanoperedens, we refer to them as mini-Borgs. Mini-Borgs are genetically diverse and can be assigned to at least five family-level groups. We identify eight families of Methanoperedens viruses, some of which encode multi-haem cytochromes, and circular ECEs encoding transposon-associated TnpB genes with proximal population-heterogeneous CRISPR arrays. These ECEs exchange genetic information with each other and with Methanoperedens, probably impacting their archaeal host activity and evolution.

Global marine microbial diversity and its potential in bioprospecting

The past two decades has witnessed a remarkable increase in the number of microbial genomes retrieved from marine systems1,2. However, it has remained challenging to translate this marine genomic diversity into biotechnological and biomedical applications3,4. Here we recovered 43,191 bacterial and archaeal genomes from publicly available marine metagenomes, encompassing a wide range of diversity with 138 distinct phyla, redefining the upper limit of marine bacterial genome size and revealing complex trade-offs between the occurrence of CRISPR-Cas systems and antibiotic resistance genes. In silico bioprospecting of these marine genomes led to the discovery of a novel CRISPR-Cas9 system, ten antimicrobial peptides, and three enzymes that degrade polyethylene terephthalate. In vitro experiments confirmed their effectiveness and efficacy. This work provides evidence that global-scale sequencing initiatives advance our understanding of how microbial diversity has evolved in the oceans and is maintained, and demonstrates how such initiatives can be sustainably exploited to advance biotechnology and biomedicine.

Metagenome-assembled microbial genomes from Parkinson's disease fecal samples

The human gut microbiome composition has been linked to Parkinson's disease (PD). However, knowledge of the gut microbiota on the genome level is still limited. Here we performed deep metagenomic sequencing and binning to build metagenome-assembled genomes (MAGs) from 136 human fecal microbiomes (68 PD samples and 68 control samples). We constructed 952 non-redundant high-quality MAGs and compared them between PD and control groups. Among these MAGs, there were 22 different genomes of Collinsella and Prevotella, indicating high variability of those genera in the human gut environment. Microdiversity analysis indicated that Ruminococcus bromii was statistically significantly (p < 0.002) more diverse on the strain level in the control samples compared to the PD samples. In addition, by clustering all genes and performing presence-absence analysis between groups, we identified several control-specific (p < 0.05) related genes, such as speF and Fe-S oxidoreductase. We also report detailed annotation of MAGs, including Clusters of Orthologous Genes (COG), Cas operon type, antiviral gene, prophage, and secondary metabolites biosynthetic gene clusters, which can be useful for providing a reference for future studies.

Unveiling the microbial realm with VEBA 2.0: a modular bioinformatics suite for end-to-end genome-resolved prokaryotic, (micro)eukaryotic and viral multi-omics from either short- or long-read sequencing

The microbiome is a complex community of microorganisms, encompassing prokaryotic (bacterial and archaeal), eukaryotic, and viral entities. This microbial ensemble plays a pivotal role in influencing the health and productivity of diverse ecosystems while shaping the web of life. However, many software suites developed to study microbiomes analyze only the prokaryotic community and provide limited to no support for viruses and microeukaryotes. Previously, we introduced the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite to address this critical gap in microbiome research by extending genome-resolved analysis beyond prokaryotes to encompass the understudied realms of eukaryotes and viruses. Here we present VEBA 2.0 with key updates including a comprehensive clustered microeukaryotic protein database, rapid genome/protein-level clustering, bioprospecting, non-coding/organelle gene modeling, genome-resolved taxonomic/pathway profiling, long-read support, and containerization. We demonstrate VEBA's versatile application through the analysis of diverse case studies including marine water, Siberian permafrost, and white-tailed deer lung tissues with the latter showcasing how to identify integrated viruses. VEBA represents a crucial advancement in microbiome research, offering a powerful and accessible software suite that bridges the gap between genomics and biotechnological solutions.

Pathogen dynamics and discovery of novel viruses and enzymes by deep nucleic acid sequencing of wastewater

Wastewater contains an extensive reservoir of genetic information, yet largely unexplored. Here, we analyzed by high-throughput sequencing total nucleic acids extracted from wastewater samples collected during a 17 month-period in Berlin, Germany. By integrating global wastewater datasets and applying a novel computational approach to accurately identify viral strains within sewage RNA-sequencing data, we demonstrated the emergence and global dissemination of a specific astrovirus strain. Astrovirus abundance and sequence variation mirrored temporal and spatial patterns of infection, potentially serving as footprints of specific timeframes and geographical locations. Additionally, we revealed more than 100,000 sequence contigs likely originating from novel viral species, exhibiting distinct profiles in total RNA and DNA datasets and including undescribed bunyaviruses and parvoviruses. Finally, we identified thousands of new CRISPR-associated protein sequences, including Transposase B (TnpB), a class of compact, RNA-guided DNA editing enzymes. Collectively, our findings underscore the potential of high-throughput sequencing of total nucleic acids derived from wastewater for a broad range of applications.

Characterization of Ligilactobacillus salivarius CRISPR-Cas systems

Ligilactobacillus is a diverse genus among lactobacilli with phenotypes that reflect adaptation to various hosts. CRISPR-Cas systems are highly prevalent within lactobacilli, and Ligilactobacillus salivarius, the most abundant species of Ligilactobacillus, possesses both DNA- and RNA-targeting CRISPR-Cas systems. In this study, we explore the presence and functional properties of I-B, I-C, I-E, II-A, and III-A CRISPR-Cas systems in over 500 Ligilactobacillus genomes, emphasizing systems found in L. salivarius. We examined the I-E, II-A, and III-A CRISPR-Cas systems of two L. salivarius strains and observed occurrences of split cas genes and differences in CRISPR RNA maturation in native hosts. This prompted testing of the single Cas9 and multiprotein Cascade and Csm CRISPR-Cas effector complexes in a cell-free context to demonstrate the functionality of these systems. We also predicted self-targeting spacers within L. salivarius CRISPR-Cas systems and found that nearly a third of L. salivarius genomes possess unique self-targeting spacers that generally target elements other than prophages. With these two L. salivarius strains, we performed prophage induction coupled with RNA sequencing and discovered that the prophages residing within these strains are inducible and likely active elements, despite targeting by CRISPR-Cas systems. These findings deepen our comprehension of CRISPR-Cas systems in L. salivarius, further elucidating their relationship with associated prophages and providing a functional basis for the repurposing of these Cas effectors for bacterial manipulation.

IMPORTANCE

Ligilactobacillus salivarius is a diverse bacterial species widely used in the food and dietary supplement industries. In this study, we investigate the occurrence and diversity of their adaptive immune systems, CRISPR-Cas, in over 500 genomes. We establish their function and provide insights into their role in the interplay between the bacterial host and the predatory phages that infect them. Such findings expand our knowledge about these important CRISPR-Cas immune systems widespread across the bacterial tree of life and also provide a technical basis for the repurposing of these molecular machines for the development of molecular biology tools and the manipulation and engineering of bacteria and other life forms.

Viral and Bacterial Community Dynamics in Food Waste and Digestate from Full-Scale Biogas Plants

Anaerobic digestion (AD) is commonly used in food waste treatment. Prokaryotic microbial communities in AD of food waste have been comprehensively studied. The role of viruses, known to affect microbial dynamics and metabolism, remains largely unexplored. This study employed metagenomic analysis and recovered 967 high-quality viral bins within food waste and digestate derived from 8 full-scale biogas plants. The diversity of viral communities was higher in digestate. In silico predictions linked 20.8% of viruses to microbial host populations, highlighting possible virus predators of key functional microbes. Lineage-specific virus-host ratio varied, indicating that viral infection dynamics might differentially affect microbial responses to the varying process parameters. Evidence for virus-mediated gene transfer was identified, emphasizing the potential role of viruses in controlling the microbiome. AD altered the specific process parameters, potentially promoting a shift in viral lifestyle from lysogenic to lytic. Viruses encoding auxiliary metabolic genes (AMGs) were involved in microbial carbon and nutrient cycling, and most AMGs were transcriptionally expressed in digestate, meaning that viruses with active functional states were likely actively involved in AD. These findings provided a comprehensive profile of viral and bacterial communities and expanded knowledge of the interactions between viruses and hosts in food waste and digestate.

Bioinformatic analysis of type III CRISPR systems reveals key properties and new effector families

Recognition of RNA from invading mobile genetic elements (MGE) prompts type III CRISPR systems to activate an HD nuclease domain and/or a nucleotide cyclase domain in the Cas10 subunit, eliciting an immune response. The cyclase domain can generate a range of nucleotide second messengers, which in turn activate a diverse family of ancillary effector proteins. These provide immunity by non-specific degradation of host and MGE nucleic acids or proteins, perturbation of membrane potentials, transcriptional responses, or the arrest of translation. The wide range of nucleotide activators and downstream effectors generates a complex picture that is gradually being resolved. Here, we carry out a global bioinformatic analysis of type III CRISPR loci in prokaryotic genomes, defining the relationships of Cas10 proteins and their ancillary effectors. Our study reveals that cyclic tetra-adenylate is by far the most common signalling molecule used and that many loci have multiple effectors. These typically share the same activator and may work synergistically to combat MGE. We propose four new candidate effector protein families and confirm experimentally that the Csm6-2 protein, a highly diverged, fused Csm6 effector, is a ribonuclease activated by cyclic hexa-adenylate.

An ecological-evolutionary perspective on the genomic diversity and habitat preferences of the Acidobacteriota

Members of the phylum Acidobacteriota inhabit a wide range of ecosystems including soils. We analyzed the global patterns of distribution and habitat preferences of various Acidobacteriota lineages across major ecosystems (soil, engineered, host-associated, marine, non-marine saline and alkaline, and terrestrial non-soil ecosystem) in 248,559 publicly available metagenomic datasets. Classes Terriglobia, Vicinamibacteria, Blastocatellia, and Thermoanaerobaculia were highly ubiquitous and showed clear preference to soil over non-soil habitats, class Polarisedimenticolia showed comparable ubiquity and preference between soil and non-soil habitats, while classes Aminicenantia and Holophagae showed preferences to non-soil habitats. However, while specific preferences were observed, most Acidobacteriota lineages were habitat generalists rather than specialists, with genomic and/or metagenomic fragments recovered from soil and non-soil habitats at various levels of taxonomic resolution. Comparative analysis of 1930 genomes strongly indicates that phylogenetic affiliation plays a more important role than the habitat from which the genome was recovered in shaping the genomic characteristics and metabolic capacities of the Acidobacteriota. The observed lack of strong habitat specialization and habitat transition driven lineage evolution in the Acidobacteriota suggest ready cross colonization between soil and non-soil habitats. We posit that such capacity is key to the successful establishment of Acidobacteriota as a major component in soil microbiomes post ecosystem disturbance events or during pedogenesis.

RNA-guided RNA silencing by an Asgard archaeal Argonaute

Argonaute proteins are the central effectors of RNA-guided RNA silencing pathways in eukaryotes, playing crucial roles in gene repression and defense against viruses and transposons. Eukaryotic Argonautes are subdivided into two clades: AGOs generally facilitate miRNA- or siRNA-mediated silencing, while PIWIs generally facilitate piRNA-mediated silencing. It is currently unclear when and how Argonaute-based RNA silencing mechanisms arose and diverged during the emergence and early evolution of eukaryotes. Here, we show that in Asgard archaea, the closest prokaryotic relatives of eukaryotes, an evolutionary expansion of Argonaute proteins took place. In particular, a deep-branching PIWI protein (HrAgo1) encoded by the genome of the Lokiarchaeon 'Candidatus Harpocratesius repetitus' shares a common origin with eukaryotic PIWI proteins. Contrasting known prokaryotic Argonautes that use single-stranded DNA as guides and/or targets, HrAgo1 mediates RNA-guided RNA cleavage, and facilitates gene silencing when expressed in human cells and supplied with miRNA precursors. A cryo-EM structure of HrAgo1, combined with quantitative single-molecule experiments, reveals that the protein displays structural features and target-binding modes that are a mix of those of eukaryotic AGO and PIWI proteins. Thus, this deep-branching archaeal PIWI may have retained an ancestral molecular architecture that preceded the functional and mechanistic divergence of eukaryotic AGOs and PIWIs.

Characteristics and immune functions of the endogenous CRISPR-Cas systems in myxobacteria

The clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) system widely occurs in prokaryotic organisms to recognize and destruct genetic invaders. Systematic collation and characterization of endogenous CRISPR-Cas systems are conducive to our understanding and potential utilization of this natural genetic machinery. In this study, we screened 39 complete and 692 incomplete genomes of myxobacteria using a combined strategy to dispose of the abridged genome information and revealed at least 19 CRISPR-Cas subtypes, which were distributed with a taxonomic difference and often lost stochastically in intraspecies strains. The cas genes in each subtype were evolutionarily clustered but deeply separated, while most of the CRISPRs were divided into four types based on the motif characteristics of repeat sequences. The spacers recorded in myxobacterial CRISPRs were in high G+C content, matching lots of phages, tiny amounts of plasmids, and, surprisingly, massive organismic genomes. We experimentally demonstrated the immune and self-target immune activities of three endogenous systems in Myxococcus xanthus DK1622 against artificial genetic invaders and revealed the microhomology-mediated end-joining mechanism for the immunity-induced DNA repair but not homology-directed repair. The panoramic view and immune activities imply potential omnipotent immune functions and applications of the endogenous CRISPR-Cas machinery.

IMPORTANCE

Serving as an adaptive immune system, clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) empower prokaryotes to fend off the intrusion of external genetic materials. Myxobacteria are a collective of swarming Gram-stain-negative predatory bacteria distinguished by intricate multicellular social behavior. An in-depth analysis of their intrinsic CRISPR-Cas systems is beneficial for our understanding of the survival strategies employed by host cells within their environmental niches. Moreover, the experimental findings presented in this study not only suggest the robust immune functions of CRISPR-Cas in myxobacteria but also their potential applications.

Type IV-A3 CRISPR-Cas systems drive inter-plasmid conflicts by acquiring spacers in trans

Plasmid-encoded type IV-A CRISPR-Cas systems lack an acquisition module, feature a DinG helicase instead of a nuclease, and form ribonucleoprotein complexes of unknown biological functions. Type IV-A3 systems are carried by conjugative plasmids that often harbor antibiotic-resistance genes and their CRISPR array contents suggest a role in mediating inter-plasmid conflicts, but this function remains unexplored. Here, we demonstrate that a plasmid-encoded type IV-A3 system co-opts the type I-E adaptation machinery from its host, Klebsiella pneumoniae (K. pneumoniae), to update its CRISPR array. Furthermore, we reveal that robust interference of conjugative plasmids and phages is elicited through CRISPR RNA-dependent transcriptional repression. By silencing plasmid core functions, type IV-A3 impacts the horizontal transfer and stability of targeted plasmids, supporting its role in plasmid competition. Our findings shed light on the mechanisms and ecological function of type IV-A3 systems and demonstrate their practical efficacy for countering antibiotic resistance in clinically relevant strains.

A panoramic view of the virosphere in three wastewater treatment plants by integrating viral-like particle-concentrated and traditional non-concentrated metagenomic approaches

Wastewater biotreatment systems harbor a rich diversity of microorganisms, and the effectiveness of biotreatment systems largely depends on the activity of these microorganisms. Specifically, viruses play a crucial role in altering microbial behavior and metabolic processes throughout their infection phases, an aspect that has recently attracted considerable interest. Two metagenomic approaches, viral-like particle-concentrated (VPC, representing free viral-like particles) and non-concentrated (NC, representing the cellular fraction), were employed to assess their efficacy in revealing virome characteristics, including taxonomy, diversity, host interactions, lifestyle, dynamics, and functional genes across processing units of three wastewater treatment plants (WWTPs). Our findings indicate that each approach offers unique insights into the viral community and functional composition. Their combined use proved effective in elucidating WWTP viromes. We identified nearly 50,000 viral contigs, with Cressdnaviricota and Uroviricota being the predominant phyla in the VPC and NC fractions, respectively. Notably, two pathogenic viral families, Asfarviridae and Adenoviridae, were commonly found in these WWTPs. We also observed significant differences in the viromes of WWTPs processing different types of wastewater. Additionally, various phage-derived auxiliary metabolic genes (AMGs) were active at the RNA level, contributing to the metabolism of the microbial community, particularly in carbon, sulfur, and phosphorus cycling. Moreover, we identified 29 virus-carried antibiotic resistance genes (ARGs) with potential for host transfer, highlighting the role of viruses in spreading ARGs in the environment. Overall, this study provides a detailed and integrated view of the virosphere in three WWTPs through the application of VPC and NC metagenomic approaches. Our findings enhance the understanding of viral communities, offering valuable insights for optimizing the operation and regulation of wastewater treatment systems.

Discarded diversity: Novel megaphages, auxiliary metabolic genes, and virally encoded CRISPR-Cas systems in landfills

Background

Viruses are the most abundant microbial entity on the planet, impacting microbial community structure and ecosystem services. Despite outnumbering Bacteria and Archaea by an order of magnitude, viruses have been comparatively underrepresented in reference databases. Metagenomic examinations have illustrated that viruses of Bacteria and Archaea have been specifically understudied in engineered environments. Here we employed metagenomic and computational biology methods to examine the diversity, host interactions, and genetic systems of viruses predicted from 27 samples taken from three municipal landfills across North America.

Results

We identified numerous viruses that are not represented in reference databases, including the third largest bacteriophage genome identified to date (~678 kbp), and note a cosmopolitan diversity of viruses in landfills that are distinct from viromes in other systems. Host-virus interactions were examined via host CRISPR spacer to viral protospacer mapping which captured hyper-targeted viral populations and six viral populations predicted to infect across multiple phyla. Virally-encoded auxiliary metabolic genes (AMGs) were identified with the potential to augment hosts' methane, sulfur, and contaminant degradation metabolisms, including AMGs not previously reported in literature. CRISPR arrays and CRISPR-Cas systems were identified from predicted viral genomes, including the two largest bacteriophage genomes to contain these genetic features. Some virally encoded Cas effector proteins appear distinct relative to previously reported Cas systems and are interesting targets for potential genome editing tools.

Conclusions

Our observations indicate landfills, as heterogeneous contaminated sites with unique selective pressures, are key locations for diverse viruses and atypical virus-host dynamics.

Exploring the Frozen Armory: Antiphage Defense Systems in Cold-Adapted Bacteria with a Focus on CRISPR-Cas Systems

Our understanding of the antiphage defense system arsenal in bacteria is rapidly expanding, but little is known about its occurrence in cold-adapted bacteria. In this study, we aim to shed light on the prevalence and distribution of antiphage defense systems in cold-adapted bacteria, with a focus on CRISPR-Cas systems. Using bioinformatics tools, Prokaryotic Antiviral Defense LOCator (PADLOC) and CRISPRCasTyper, we mapped the presence and diversity of antiphage defense systems in 938 available genomes of cold-adapted bacteria from diverse habitats. We confirmed that CRISPR-Cas systems are less frequent in cold-adapted bacteria, compared to mesophilic and thermophilic species. In contrast, several antiphage defense systems, such as dXTPases and DRTs, appear to be more frequently compared to temperate bacteria. Additionally, our study provides Cas endonuclease candidates with a potential for further development into cold-active CRISPR-Cas genome editing tools. These candidates could have broad applications in research on cold-adapted organisms. Our study provides a first-time map of antiphage defense systems in cold-adapted bacteria and a detailed overview of CRISPR-Cas diversity.

Specificities and commonalities of the Planctomycetes plasmidome

Plasmids, despite their critical role in antibiotic resistance and modern biotechnology, are understood in only a few bacterial groups in terms of their natural ecological dynamics. The bacterial phylum Planctomycetes, known for its unique molecular and cellular biology, has a largely unexplored plasmidome. This study offers a thorough exploration of the diversity of natural plasmids within Planctomycetes, which could serve as a foundation for developing various genetic research tools for this phylum. Planctomycetes plasmids encode a broad range of biological functions and appear to have coevolved significantly with their host chromosomes, sharing many homologues. Recent transfer events of insertion sequences between cohabiting chromosomes and plasmids were also observed. Interestingly, 64% of plasmid genes are distantly related to either chromosomally encoded genes or have homologues in plasmids from other bacterial groups. The planctomycetal plasmidome is composed of 36% exclusive proteins. Most planctomycetal plasmids encode a replication initiation protein from the Replication Protein A family near a putative iteron-containing replication origin, as well as active type I partition systems. The identification of one conjugative and three mobilizable plasmids suggests the occurrence of horizontal gene transfer via conjugation within this phylum. This comprehensive description enhances our understanding of the plasmidome of Planctomycetes and its potential implications in antibiotic resistance and biotechnology.

CoCas9 is a compact nuclease from the human microbiome for efficient and precise genome editing

The expansion of the CRISPR-Cas toolbox is highly needed to accelerate the development of therapies for genetic diseases. Here, through the interrogation of a massively expanded repository of metagenome-assembled genomes, mostly from human microbiomes, we uncover a large variety (n = 17,173) of type II CRISPR-Cas loci. Among these we identify CoCas9, a strongly active and high-fidelity nuclease with reduced molecular size (1004 amino acids) isolated from an uncultivated Collinsella species. CoCas9 is efficiently co-delivered with its sgRNA through adeno associated viral (AAV) vectors, obtaining efficient in vivo editing in the mouse retina. With this study we uncover a collection of previously uncharacterized Cas9 nucleases, including CoCas9, which enriches the genome editing toolbox.

Diverse and abundant phages exploit conjugative plasmids

Phages exert profound evolutionary pressure on bacteria by interacting with receptors on the cell surface to initiate infection. While the majority of phages use chromosomally encoded cell surface structures as receptors, plasmid-dependent phages exploit plasmid-encoded conjugation proteins, making their host range dependent on horizontal transfer of the plasmid. Despite their unique biology and biotechnological significance, only a small number of plasmid-dependent phages have been characterized. Here we systematically search for new plasmid-dependent phages targeting IncP and IncF plasmids using a targeted discovery platform, and find that they are common and abundant in wastewater, and largely unexplored in terms of their genetic diversity. Plasmid-dependent phages are enriched in non-canonical types of phages, and all but one of the 65 phages we isolated were non-tailed, and members of the lipid-containing tectiviruses, ssDNA filamentous phages or ssRNA phages. We show that plasmid-dependent tectiviruses exhibit profound differences in their host range which is associated with variation in the phage holin protein. Despite their relatively high abundance in wastewater, plasmid-dependent tectiviruses are missed by metaviromic analyses, underscoring the continued importance of culture-based phage discovery. Finally, we identify a tailed phage dependent on the IncF plasmid, and find related structural genes in phages that use the orthogonal type 4 pilus as a receptor, highlighting the evolutionarily promiscuous use of these distinct contractile structures by multiple groups of phages. Taken together, these results indicate plasmid-dependent phages play an under-appreciated evolutionary role in constraining horizontal gene transfer via conjugative plasmids.

Diverse Class 2 CRISPR Effectors as Active Nucleases with Expanded Targeting Capabilities

CRISPR-Cas systems have proven effective in a variety of applications due to their ease of use and relatively high editing efficiency. Yet, any individual CRISPR-Cas system has inherent limitations, necessitating a diversity of RNA-guided nucleases to suit applications with distinct needs. We searched through metagenomic sequences to identify RNA-guided nucleases and found enzymes from diverse CRISPR-Cas types and subtypes, the most promising of which we developed into gene-editing platforms. Based on prior annotations of the metagenomic sequences, we establish the likely taxa and sampling locations where Class 2 CRISPR-Cas systems active in eukaryotes may be found. The newly discovered systems show robust capabilities as gene editors and base editors.

Discovery and engineering of AiEvo2, a novel Cas12a nuclease for human gene editing applications

The precision of gene editing technology is critical to creating safe and effective therapies for treating human disease. While the programmability of CRISPR-Cas systems has allowed for rapid innovation of new gene editing techniques, the off-target activity of these enzymes has hampered clinical development for novel therapeutics. Here, we report the identification and characterization of a novel CRISPR-Cas12a enzyme from Acinetobacter indicus (AiCas12a). We engineer the nuclease (termed AiEvo2) for increased specificity, protospacer adjacent motif recognition, and efficacy on a variety of human clinical targets. AiEvo2 is highly precise and able to efficiently discriminate between normal and disease-causing alleles in Huntington's patient-derived cells by taking advantage of a single nucleotide polymorphism on the disease-associated allele. AiEvo2 efficiently edits several liver-associated target genes including PCSK9 and TTR when delivered to primary hepatocytes as mRNA encapsulated in a lipid nanoparticle. The enzyme also engineers an effective CD19 chimeric antigen receptor-T-cell therapy from primary human T cells using multiplexed simultaneous editing and chimeric antigen receptor insertion. To further ensure precise editing, we engineered an anti-CRISPR protein to selectively inhibit off-target gene editing while retaining therapeutic on-target editing. The engineered AiEvo2 nuclease coupled with a novel engineered anti-CRISPR protein represents a new way to control the fidelity of editing and improve the safety and efficacy of gene editing therapies.

Pangenome analysis of Shewanella xiamenensis revealed important genetic traits concerning genetic diversity, pathogenicity and antibiotic resistance

BACKGROUND

Shewanella xiamenensis, widely distributed in natural environments, has long been considered as opportunistic pathogen. Recently, significant changes in the resistance spectrum have been observed in S. xiamenensis, due to acquired antibiotic resistance genes. Therefore, a pan-genome analysis was conducted to illuminate the genomic changes in S. xiamenensis.

RESULTS

Phylogenetic analysis revealed three major clusters and three singletons, among which close relationship between several strains was discovered, regardless of their host and niches. The "open" genomes with diversity of accessory and strain-specific genomes took advantage towards diversity environments. The purifying selection pressure was the main force on genome evolution, especially in conservative genes. Only 53 gene families were under positive selection pressure. Phenotypic resistance analysis revealed 21 strains were classified as multi-drug resistance (MDR). Ten types of antibiotic resistance genes and two heavy metal resistance operons were discovered in S. xiamenensis. Mobile genetic elements and horizontal gene transfer increased genome diversity and were closely related to MDR strains. S. xiamenensis carried a variety of virulence genes and macromolecular secretion systems, indicating their important roles in pathogenicity and adaptability. Type IV secretion system was discovered in 15 genomes with various sequence structures, indicating it was originated from different donors through horizontal gene transfer.

CONCLUSIONS

This study provided with a detailed insight into the changes in the pan-genome of S. xiamenensis, highlighting its capability to acquire new mobile genetic elements and resistance genes for its adaptation to environment and pathogenicity to human and animals.

Increased mutations in lipopolysaccharide biosynthetic genes cause time-dependent development of phage resistance in Salmonella

Understanding how bacteria evolve resistance to phages has implications for phage-based therapies and microbial evolution. In this study, the susceptibility of 335 Salmonella isolates to the wide host range Salmonella phage BPSELC-1 was tested. Potentially significant gene sets that could confer resistance were identified using bioinformatics approaches based on phage susceptibility phenotypes; more than 90 potential antiphage defense gene sets, including those involved in lipopolysaccharide (LPS) biosynthesis, DNA replication, secretion systems, and respiratory chain, were found. The evolutionary dynamics of Salmonella resistance to phage were assessed through laboratory evolution experiments, which showed that phage-resistant mutants rapidly developed and exhibited genetic heterogeneity. Most representative Salmonella hosts (58.1% of 62) rapidly developed phage resistance within 24 h. All phage-resistant mutant clones exhibited genetic heterogeneity and observed mutations in LPS-related genes (rfaJ and rfaK) as well as other genes such as cellular respiration, transport, and cell replication-related genes. The study also identified potential trade-offs, indicating that bacteria tend to escape fitness trade-offs through multi-site mutations, all tested mutants increased sensitivity to polymyxin B, but this does not always affect their relative fitness or biofilm-forming capacity. Furthermore, complementing the rfaJ mutant gene could partially restore the phage sensitivity of phage-resistant mutants. These results provide insight into the phage resistance mechanisms of Salmonella and the complexity of bacterial evolution resulting from phage predation, which can inform future strategies for phage-based therapies and microbial evolution.

Comparative evaluation of bioinformatic tools for virus-host prediction and their application to a highly diverse community in the Cuatro Ciénegas Basin, Mexico

Due to the enormous diversity of non-culturable viruses, new viruses must be characterized using culture-independent techniques. The associated host is an important phenotypic feature that can be inferred from metagenomic viral contigs thanks to the development of several bioinformatic tools. Here, we compare the performance of recently developed virus-host prediction tools on a dataset of 1,046 virus-host pairs and then apply the best-performing tools to a metagenomic dataset derived from a highly diverse transiently hypersaline site known as the Archaean Domes (AD) within the Cuatro Ciénegas Basin, Coahuila, Mexico. Among host-dependent methods, alignment-based approaches had a precision of 66.07% and a sensitivity of 24.76%, while alignment-free methods had an average precision of 75.7% and a sensitivity of 57.5%. RaFAH, a virus-dependent alignment-based tool, had the best overall performance (F1_score = 95.7%). However, when predicting the host of AD viruses, methods based on public reference databases (such as RaFAH) showed lower inter-method agreement than host-dependent methods run against custom databases constructed from prokaryotes inhabiting AD. Methods based on custom databases also showed the greatest agreement between the source environment and the predicted host taxonomy, habitat, lifestyle, or metabolism. This highlights the value of including custom data when predicting hosts on a highly diverse metagenomic dataset, and suggests that using a combination of methods and qualitative validations related to the source environment and predicted host biology can increase the number of correct predictions. Finally, these predictions suggest that AD viruses infect halophilic archaea as well as a variety of bacteria that may be halophilic, halotolerant, alkaliphilic, thermophilic, oligotrophic, sulfate-reducing, or marine, which is consistent with the specific environment and the known geological and biological evolution of the Cuatro Ciénegas Basin and its microorganisms.

CRISPR genotyping methods: Tracing the evolution from spoligotyping to machine learning

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) systems provide prokaryotes with adaptive immunity defenses against foreign genetic invaders. The identification of CRISPR-Cas function is among the most impactful discoveries of recent decades that have shaped the development of genome editing in various organisms paving the way for a plethora of promising applications in biotechnology and health. Even before the discovery of CRISPR-Cas biological role, the particular structure of CRISPR loci has been explored for epidemiological genotyping of bacterial pathogens. CRISPR-Cas loci are arranged in CRISPR arrays of mostly identical direct repeats intercalated with invader-derived spacers and an operon of cas genes encoding the Cas protein components. Each small CRISPR RNA (crRNA) encoded within the CRISPR array constitutes a key functional unit of this RNA-based CRISPR-Cas defense system guiding the Cas effector proteins toward the foreign nucleic acids for their destruction. The information acquired from prior invader encounters and stored within CRISPR arrays turns out to be extremely valuable in tracing the microevolution and epidemiology of major bacterial pathogens. We review here the history of CRISPR-based typing strategies highlighting the first PCR-based methods that have set the stage for recent developments of high-throughput sequencing and machine learning-based approaches. A great amount of whole genome sequencing and metagenomic data accumulated in recent years opens up new avenues for combining experimental and computational approaches of high-resolution CRISPR-based typing.

Viral communities in millipede guts: Insights into the diversity and potential role in modulating the microbiome

Millipedes are important detritivores harbouring a diverse microbiome. Previous research focused on bacterial and archaeal diversity, while the virome remained neglected. We elucidated the DNA and RNA viral diversity in the hindguts of two model millipede species with distinct microbiomes: the tropical Epibolus pulchripes (methanogenic, dominated by Bacillota) and the temperate Glomeris connexa (non-methanogenic, dominated by Pseudomonadota). Based on metagenomic and metatranscriptomic assembled viral genomes, the viral communities differed markedly and preferentially infected the most abundant prokaryotic taxa. The majority of DNA viruses were Caudoviricetes (dsDNA), Cirlivirales (ssDNA) and Microviridae (ssDNA), while RNA viruses consisted of Leviviricetes (ssRNA), Potyviridae (ssRNA) and Eukaryotic viruses. A high abundance of subtypes I-C, I-B and II-C CRISPR-Cas systems was found, primarily from Pseudomonadota, Bacteroidota and Bacillota. In addition, auxiliary metabolic genes that modulate chitin degradation, vitamins and amino acid biosynthesis and sulphur metabolism were also detected. Lastly, we found low virus-to-microbe-ratios and a prevalence of lysogenic viruses, supporting a Piggyback-the-Winner dynamic in both hosts.

Dynamics of CRISPR-mediated virus-host interactions in the human gut microbiome

Arms races between mobile genetic elements and prokaryotic hosts are major drivers of ecological and evolutionary change in microbial communities. Prokaryotic defense systems such as CRISPR-Cas have the potential to regulate microbiome composition by modifying the interactions among bacteria, plasmids, and phages. Here, we used longitudinal metagenomic data from 130 healthy and diseased individuals to study how the interplay of genetic parasites and CRISPR-Cas immunity reflects on the dynamics and composition of the human gut microbiome. Based on the coordinated study of 80 000 CRISPR-Cas loci and their targets, we show that CRISPR-Cas immunity effectively modulates bacteriophage abundances in the gut. Acquisition of CRISPR-Cas immunity typically leads to a decrease in the abundance of lytic phages but does not necessarily cause their complete disappearance. Much smaller effects are observed for lysogenic phages and plasmids. Conversely, phage-CRISPR interactions shape bacterial microdiversity by producing weak selective sweeps that benefit immune host lineages. We also show that distal (and chronologically older) regions of CRISPR arrays are enriched in spacers that are potentially functional and target crass-like phages and local prophages. This suggests that exposure to reactivated prophages and other endemic viruses is a major selective pressure in the gut microbiome that drives the maintenance of long-lasting immune memory.

Deepdefense: annotation of immune systems in prokaryotes using deep learning

BACKGROUND

Due to a constant evolutionary arms race, archaea and bacteria have evolved an abundance and diversity of immune responses to protect themselves against phages. Since the discovery and application of CRISPR-Cas adaptive immune systems, numerous novel candidates for immune systems have been identified. Previous approaches to identifying these new immune systems rely on hidden Markov model (HMM)-based homolog searches or use labor-intensive and costly wet-lab experiments. To aid in finding and classifying immune systems genomes, we use machine learning to classify already known immune system proteins and discover potential candidates in the genome. Neural networks have shown promising results in classifying and predicting protein functionality in recent years. However, these methods often operate under the closed-world assumption, where it is presumed that all potential outcomes or classes are already known and included in the training dataset. This assumption does not always hold true in real-world scenarios, such as in genomics, where new samples can emerge that were not previously accounted for in the training phase.

RESULTS

In this work, we explore neural networks for immune protein classification, deal with different methods for rejecting unrelated proteins in a genome-wide search, and establish a benchmark. Then, we optimize our approach for accuracy. Based on this, we develop an algorithm called Deepdefense to predict immune cassette classes based on a genome. This design facilitates the differentiation between immune system-related and unrelated proteins by analyzing variations in model-predicted confidence values, aiding in the identification of both known and potentially novel immune system proteins. Finally, we test our approach for detecting immune systems in the genome against an HMM-based method.

CONCLUSIONS

Deepdefense can automatically detect genes and define cassette annotations and classifications using 2 model classifications. This is achieved by creating an optimized deep learning model to annotate immune systems, in combination with calibration methods, and a second model to enable the scanning of an entire genome.

The CRISPR-Cas system in clinical strains of Acinetobacter baumannii: an in-silico analysis

Acinetobacter baumannii is a relevant bacterium due to its high-resistance profile. It is well known that antimicrobial resistance is primarily linked to mutations and the acquisition of external genomic material, such as plasmids or phages, to which the Clustered Regularly Interspaced Short Palindromic Repeats associated with Cas proteins, or CRISPR-Cas, system is related. It is known that the system can influence the acquisition of foreign genetic material and play a role in various physiological pathways. In this study, we conducted an in-silico analysis using 91 fully assembled genomes of clinical strains obtained from the NCBI database. Among the analyzed genomes, the I-F1 subtype of the CRISPR-Cas system was detected showcasing variations in architecture and phylogeny. Using bioinformatic tools, we determined the presence, distribution, and specific characteristics of the CRISPR-Cas system. We found a possible association of the system with resistance genes but not with virulence determinants. Analysis of the system's components, including spacer sequences, suggests its potential role in protecting against phage infections, highlighting its protective function.

Diversity of Bathyarchaeia viruses in metagenomes and virus-encoded CRISPR system components

Bathyarchaeia represent a class of archaea common and abundant in sedimentary ecosystems. Here we report 56 metagenome-assembled genomes of Bathyarchaeia viruses identified in metagenomes from different environments. Gene sharing network and phylogenomic analyses led to the proposal of four virus families, including viruses of the realms Duplodnaviria and Adnaviria, and archaea-specific spindle-shaped viruses. Genomic analyses uncovered diverse CRISPR elements in these viruses. Viruses of the proposed family "Fuxiviridae" harbor an atypical Type IV-B CRISPR-Cas system and a Cas4 protein that might interfere with host immunity. Viruses of the family "Chiyouviridae" encode a Cas2-like endonuclease and two mini-CRISPR arrays, one with a repeat identical to that in the host CRISPR array, potentially allowing the virus to recruit the host CRISPR adaptation machinery to acquire spacers that could contribute to competition with other mobile genetic elements or to inhibit host defenses. These findings present an outline of the Bathyarchaeia virome and offer a glimpse into their counter-defense mechanisms.

Evaluation of Safety and Probiotic Traits from a Comprehensive Genome-Based In Silico Analysis of Ligilactobacillus salivarius P1CEA3, Isolated from Pigs and Producer of Nisin S

Ligilactobacillus salivarius is an important member of the porcine gastrointestinal tract (GIT). Some L. salivarius strains are considered to have a beneficial effect on the host by exerting different probiotic properties, including the production of antimicrobial peptides which help maintain a healthy gut microbiota. L. salivarius P1CEA3, a porcine isolated strain, was first selected and identified by its antimicrobial activity against a broad range of pathogenic bacteria due to the production of the novel bacteriocin nisin S. The assembled L. salivarius P1CEA3 genome includes a circular chromosome, a megaplasmid (pMP1CEA3) encoding the nisin S gene cluster, and two small plasmids. A comprehensive genome-based in silico analysis of the L. salivarius P1CEA3 genome reveals the presence of genes related to probiotic features such as bacteriocin synthesis, regulation and production, adhesion and aggregation, the production of lactic acid, amino acids metabolism, vitamin biosynthesis, and tolerance to temperature, acid, bile salts and osmotic and oxidative stress. Furthermore, the strain is absent of risk-related genes for acquired antibiotic resistance traits, virulence factors, toxic metabolites and detrimental metabolic or enzymatic activities. Resistance to common antibiotics and gelatinase and hemolytic activities have been discarded by in vitro experiments. This study identifies several probiotic and safety traits of L. salivarius P1CEA3 and suggests its potential as a promising probiotic in swine production.

A New RNA-Dependent Cas12g Nuclease

The development of RNA-targeting CRISPR-Cas systems represents a major step forward in the field of gene editing and regulation. RNA editing presents a viable alternative to genome editing in certain scenarios as it offers a reversible and manageable approach, reducing the likelihood of runaway mutant variants. One of the most promising applications is in the treatment of genetic disorders caused by mutations in RNA molecules. In this study, we investigate a previously undescribed Cas12g nuclease which was found in metagenomes from promising thermophilic microbial communities during the expedition to the Republic of North Ossetia-Alania in 2020. The method outlined in this study can be applied to other Cas orthologs and variants, leading to a better understanding of the CRISPR-Cas system and its enzymatic activities. The cis-cleavage activity of the new type V-G Cas effector was indicated by in vitro RNA cleavage experiments. While CRISPR-Cas systems are known for their high specificity, there is still a risk of unintended cleavage of nontargeted RNA molecules. Ultimately, the search for new genome editing tools and the study of their properties will remove barriers to research in this area. With continued research and development, we may be able to unlock their full potential.

In-silico whole-genome sequence analysis of a halotolerant filamentous mangrove cyanobacterium revealed CRISPR-Cas systems with unique properties

Novel CRISPR systems capable of cleaving both DNA and RNA are progressively emerging as attractive tools for genome manipulation of prokaryotic and eukaryotic organisms. We report specific characteristics of CRISPR systems present in Oxynema aestuarii AP17, a halotolerant, filamentous cyanobacterium and the second known member of the Oxynema genus. In-silico analyses of its whole-genome sequence revealed the presence of multiple Type I and Type III CRISPR loci with one Type I-G system previously unreported in cyanobacteria. We further identified the leader sequences at the 5' end of multiple CRISPR loci, many of which were distinct from previously reported cyanobacterial CRISPR leaders. Phylogenetic analyses of the O. aestuarii AP17 Cas1 proteins revealed two protein sequences that were unique and distantly related to other cyanobacterial Cas1 protein sequences. Our findings are significant because novel Class 1 CRISPR systems possess multi-subunit effectors and are highly flexible for repurposing by protein domain fusions made to the effector complex. Additionally, Type III CRISPRs are particularly useful for genome editing in certain extremophiles for which mesophilic Type II CRISPRs are ineffective.

Endogenous CRISPR-Cas mediated in situ genome editing: State-of-the-art and the road ahead for engineering prokaryotes

The CRISPR-Cas systems have shown tremendous promise as heterologous tools for genome editing in various prokaryotes. However, the perturbation of DNA homeostasis and the inherent toxicity of Cas9/12a proteins could easily lead to cell death, which led to the development of endogenous CRISPR-Cas systems. Programming the widespread endogenous CRISPR-Cas systems for in situ genome editing represents a promising tool in prokaryotes, especially in genetically intractable species. Here, this review briefly summarizes the advances of endogenous CRISPR-Cas-mediated genome editing, covering aspects of establishing and optimizing the genetic tools. In particular, this review presents the application of different types of endogenous CRISPR-Cas tools for strain engineering, including genome editing and genetic regulation. Notably, this review also provides a detailed discussion of the transposon-associated CRISPR-Cas systems, and the programmable RNA-guided transposition using endogenous CRISPR-Cas systems to enable editing of microbial communities for understanding and control. Therefore, they will be a powerful tool for targeted genetic manipulation. Overall, this review will not only facilitate the development of standard genetic manipulation tools for non-model prokaryotes but will also enable more non-model prokaryotes to be genetically tractable.

Characterization and diversity of CRISPR/Cas systems in Klebsiella oxytoca

CRISPR (clustered regularly interspaced short palindromic repeats)/Cas (CRISPR-associated protein) system is a crucial adaptive immune system for bacteria to resist foreign DNA infection. In this study, we investigated the prevalence and diversity of CRISPR/Cas systems in 175 Klebsiella oxytoca (K. oxytoca) strains. Specifically, 58.86% (103/175) of these strains possessed at least one confirmed CRISPR locus. Two CRISPR/Cas system types, I-F and IV-A3, were identified in 69 strains. Type I-F system was the most prevalent in this species, which correlated well with MLST. Differently, type IV-A3 system was randomly distributed. Moreover, the type IV-A3 system was separated into two subgroups, with subgroup-specific cas genes and repeat sequences. In addition, spacer origin analysis revealed that approximately one-fifth of type I-F spacers and one-third of type IV-A3 spacers had a significant match to MGEs. The phage tail tape measure protein and conjunctive transfer system protein were important targets of type I-F and IV-A3 systems in K. oxytoca, respectively. PAM sequences were inferred to be 5'-NCC-3' for type I-F, 5'-AAG-3' for subgroup IV-A3-a, and 5'-AAN-3' for subgroup IV-A3-b. Collectively, our findings will shed light on the prevalence, diversity, and functional effects of the CRISPR/Cas system in K. oxytoca.

Bacteriophages suppress CRISPR-Cas immunity using RNA-based anti-CRISPRs

Many bacteria use CRISPR-Cas systems to combat mobile genetic elements, such as bacteriophages and plasmids1. In turn, these invasive elements have evolved anti-CRISPR proteins to block host immunity2,3. Here we unveil a distinct type of CRISPR-Cas Inhibition strategy that is based on small non-coding RNA anti-CRISPRs (Racrs). Racrs mimic the repeats found in CRISPR arrays and are encoded in viral genomes as solitary repeat units4. We show that a prophage-encoded Racr strongly inhibits the type I-F CRISPR-Cas system by interacting specifically with Cas6f and Cas7f, resulting in the formation of an aberrant Cas subcomplex. We identified Racr candidates for almost all CRISPR-Cas types encoded by a diverse range of viruses and plasmids, often in the genetic context of other anti-CRISPR genes5. Functional testing of nine candidates spanning the two CRISPR-Cas classes confirmed their strong immune inhibitory function. Our results demonstrate that molecular mimicry of CRISPR repeats is a widespread anti-CRISPR strategy, which opens the door to potential biotechnological applications6.

Comparative Analysis and Phylogenetic Insights of Cas14-Homology Proteins in Bacteria and Archaea

Type-V-F Cas12f proteins, also known as Cas14, have drawn significant interest within the diverse CRISPR-Cas nucleases due to their compact size. This study involves analyzing and comparing Cas14-homology proteins in prokaryotic genomes through mining, sequence comparisons, a phylogenetic analysis, and an array/repeat analysis. In our analysis, we identified and mined a total of 93 Cas14-homology proteins that ranged in size from 344 aa to 843 aa. The majority of the Cas14-homology proteins discovered in this analysis were found within the Firmicutes group, which contained 37 species, representing 42% of all the Cas14-homology proteins identified. In archaea, the DPANN group had the highest number of species containing Cas14-homology proteins, a total of three species. The phylogenetic analysis results demonstrate the division of Cas14-homology proteins into three clades: Cas14-A, Cas14-B, and Cas14-U. Extensive similarity was observed at the C-terminal end (CTD) through a domain comparison of the three clades, suggesting a potentially shared mechanism of action due to the presence of cutting domains in that region. Additionally, a sequence similarity analysis of all the identified Cas14 sequences indicated a low level of similarity (18%) between the protein variants. The analysis of repeats/arrays in the extended nucleotide sequences of the identified Cas14-homology proteins highlighted that 44 out of the total mined proteins possessed CRISPR-associated repeats, with 20 of them being specific to Cas14. Our study contributes to the increased understanding of Cas14 proteins across prokaryotic genomes. These homologous proteins have the potential for future applications in the mining and engineering of Cas14 proteins.

Antiviral type III CRISPR signalling via conjugation of ATP and SAM

CRISPR systems are widespread in the prokaryotic world, providing adaptive immunity against mobile genetic elements1,2. Type III CRISPR systems, with the signature gene cas10, use CRISPR RNA to detect non-self RNA, activating the enzymatic Cas10 subunit to defend the cell against mobile genetic elements either directly, via the integral histidine-aspartate (HD) nuclease domain3-5 or indirectly, via synthesis of cyclic oligoadenylate second messengers to activate diverse ancillary effectors6-9. A subset of type III CRISPR systems encode an uncharacterized CorA-family membrane protein and an associated NrN family phosphodiesterase that are predicted to function in antiviral defence. Here we demonstrate that the CorA-associated type III-B (Cmr) CRISPR system from Bacteroides fragilis provides immunity against mobile genetic elements when expressed in Escherichia coli. However, B. fragilis Cmr does not synthesize cyclic oligoadenylate species on activation, instead generating S-adenosyl methionine (SAM)-AMP (SAM is also known as AdoMet) by conjugating ATP to SAM via a phosphodiester bond. Once synthesized, SAM-AMP binds to the CorA effector, presumably leading to cell dormancy or death by disruption of the membrane integrity. SAM-AMP is degraded by CRISPR-associated phosphodiesterases or a SAM-AMP lyase, potentially providing an 'off switch' analogous to cyclic oligoadenylate-specific ring nucleases10. SAM-AMP thus represents a new class of second messenger for antiviral signalling, which may function in different roles in diverse cellular contexts.

Metabolic versatility of Caldarchaeales from geothermal features of Hawai'i and Chile as revealed by five metagenome-assembled genomes

Members of the archaeal order Caldarchaeales (previously the phylum Aigarchaeota) are poorly sampled and are represented in public databases by relatively few genomes. Additional representative genomes will help resolve their placement among all known members of Archaea and provide insights into their roles in the environment. In this study, we analyzed 16S rRNA gene amplicons belonging to the Caldarchaeales that are available in public databases, which demonstrated that archaea of the order Caldarchaeales are diverse, widespread, and most abundant in geothermal habitats. We also constructed five metagenome-assembled genomes (MAGs) of Caldarchaeales from two geothermal features to investigate their metabolic potential and phylogenomic position in the domain Archaea. Two of the MAGs were assembled from microbial community DNA extracted from fumarolic lava rocks from Mauna Ulu, Hawai'i, and three were assembled from DNA obtained from hot spring sinters from the El Tatio geothermal field in Chile. MAGs from Hawai'i are high quality bins with completeness >95% and contamination <1%, and one likely belongs to a novel species in a new genus recently discovered at a submarine volcano off New Zealand. MAGs from Chile have lower completeness levels ranging from 27 to 70%. Gene content of the MAGs revealed that these members of Caldarchaeales are likely metabolically versatile and exhibit the potential for both chemoorganotrophic and chemolithotrophic lifestyles. The wide array of metabolic capabilities exhibited by these members of Caldarchaeales might help them thrive under diverse harsh environmental conditions. All the MAGs except one from Chile harbor putative prophage regions encoding several auxiliary metabolic genes (AMGs) that may confer a fitness advantage on their Caldarchaeales hosts by increasing their metabolic potential and make them better adapted to new environmental conditions. Phylogenomic analysis of the five MAGs and over 3,000 representative archaeal genomes showed the order Caldarchaeales forms a monophyletic group that is sister to the clade comprising the orders Geothermarchaeales (previously Candidatus Geothermarchaeota), Conexivisphaerales and Nitrososphaerales (formerly known as Thaumarchaeota), supporting the status of Caldarchaeales members as a clade distinct from the Thaumarchaeota.

Staphylococcus brunensis sp. nov. isolated from human clinical specimens with a staphylococcal cassette chromosome-related genomic island outside of the rlmH gene bearing the ccrDE recombinase gene complex

Novel species of coagulase-negative staphylococci, which could serve as reservoirs of virulence and antimicrobial resistance factors for opportunistic pathogens from the genus Staphylococcus, are recognized in human and animal specimens due to advances in diagnostic techniques. Here, we used whole-genome sequencing, extensive biotyping, MALDI-TOF mass spectrometry, and chemotaxonomy to characterize five coagulase-negative strains from the Staphylococcus haemolyticus phylogenetic clade obtained from human ear swabs, wounds, and bile. Based on the results of polyphasic taxonomy, we propose the species Staphylococcus brunensis sp. nov. (type strain NRL/St 16/872T = CCM 9024T = LMG 31872T = DSM 111349T). The genomic analysis revealed numerous variable genomic elements, including staphylococcal cassette chromosome (SCC), prophages, plasmids, and a unique 18.8 kb-long genomic island SbCIccrDE integrated into the ribosomal protein L7 serine acetyltransferase gene rimL. SbCIccrDE has a cassette chromosome recombinase (ccr) gene complex with a typical structure found in SCCs. Based on nucleotide and amino acid identity to other known ccr genes and the distinct integration site that differs from the canonical methyltransferase gene rlmH exploited by SCCs, we classified the ccr genes as novel variants, ccrDE. The comparative genomic analysis of SbCIccrDE with related islands shows that they can accumulate virulence and antimicrobial resistance factors creating novel resistance elements, which reflects the evolution of SCC. The spread of these resistance islands into established pathogens such as Staphylococcus aureus would pose a great threat to the healthcare system. IMPORTANCE The coagulase-negative staphylococci are important opportunistic human pathogens, which cause bloodstream and foreign body infections, mainly in immunocompromised patients. The mobile elements, primarily the staphylococcal cassette chromosome mec, which confers resistance to methicillin, are the key to the successful dissemination of staphylococci into healthcare and community settings. Here, we present a novel species of the Staphylococcus genus isolated from human clinical material. The detailed analysis of its genome revealed a previously undescribed genomic island, which is closely related to the staphylococcal cassette chromosome and has the potential to accumulate and spread virulence and resistance determinants. The island harbors a set of conserved genes required for its mobilization, which we recognized as novel cassette chromosome recombinase genes ccrDE. Similar islands were revealed not only in the genomes of coagulase-negative staphylococci but also in S. aureus. The comparative genomic study contributes substantially to the understanding of the evolution and pathogenesis of staphylococci.

Evidence of host specificity in Lactobacillus johnsonii genomes and its influence on probiotic potential in poultry

To date, the selection of candidate strains for probiotic development in production animals has been largely based upon screens for desired phenotypic traits. However, increasing evidence indicates that the use of host-specific strains may be important, because coevolution with the animal host better prepares a bacterial strain to colonize and succeed in its respective host animal species. This concept was applied to Lactobacillus johnsonii in commercial poultry production because of its previous correlation with enhanced bird performance. Using 204 naturally isolated chicken- and turkey-source L. johnsonii, we demonstrate that there is a strong phylogenetic signal for coevolution with the animal host. These isolates differ phenotypically, even within host source, and these differences can be correlated with certain L. johnsonii phylogenetic clades. In commercial turkey poults, turkey-specific strains with strong in vitro phenotypes performed better early in life than strains lacking those phenotypes. A follow-up performance trial in broiler chickens demonstrated that chicken-specific strains result in better overall bird performance than nonchicken-specific strains. Collectively, this work provides evidence for the impact of host adaptation on a probiotic strain's potential. Furthermore, this top-down approach is useful for screening larger numbers of isolates for probiotic candidates.