CRISPR-Cas systems are widespread accessory elements across bacterial and archaeal plasmids.

Many prokaryotes encode CRISPR-Cas systems as immune protection against mobile genetic elements (MGEs), yet a number of MGEs also harbor CRISPR-Cas components. With a few exceptions, CRISPR-Cas loci encoded on MGEs are uncharted and a comprehensive analysis of their distribution, prevalence, diversity, and function is lacking. Here, we systematically investigated CRISPR-Cas loci across the largest curated collection of natural bacterial and archaeal plasmids. CRISPR-Cas loci are widely but heterogeneously distributed across plasmids and, in comparison to host chromosomes, their mean prevalence per Mbp is higher and their distribution is distinct. Furthermore, the spacer content of plasmid CRISPRs exhibits a strong targeting bias towards other plasmids, while chromosomal arrays are enriched with virus-targeting spacers. These contrasting targeting preferences highlight the genetic independence of plasmids and suggest a major role for mediating plasmid-plasmid conflicts. Altogether, CRISPR-Cas are frequent accessory components of many plasmids, which is an overlooked phenomenon that possibly facilitates their dissemination across microbiomes.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions.

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Streamlining CRISPR spacer-based bacterial host predictions to decipher the viral dark matter.

Thousands of new phages have recently been discovered thanks to viral metagenomics. These phages are extremely diverse and their genome sequences often do not resemble any known phages. To appreciate their ecological impact, it is important to determine their bacterial hosts. CRISPR spacers can be used to predict hosts of unknown phages, as spacers represent biological records of past phage-bacteria interactions. However, no guidelines have been established to standardize host prediction based on CRISPR spacers. Additionally, there are no tools that use spacers to perform host predictions on large viral datasets. Here, we developed a set of tools that includes all the necessary steps for predicting the hosts of uncharacterized phages. We created a database of >11 million spacers and a program to execute host predictions on large viral datasets. Our host prediction approach uses biological criteria inspired by how CRISPR-Cas naturally work as adaptive immune systems, which make the results easy to interpret. We evaluated the performance using 9484 phages with known hosts and obtained a recall of 49% and a precision of 69%. We also found that this host prediction method yielded higher performance for phages that infect gut-associated bacteria, suggesting it is well suited for gut-virome characterization.

CRISPRloci: comprehensive and accurate annotation of CRISPR-Cas systems.

CRISPR-Cas systems are adaptive immune systems in prokaryotes, providing resistance against invading viruses and plasmids. The identification of CRISPR loci is currently a non-standardized, ambiguous process, requiring the manual combination of multiple tools, where existing tools detect only parts of the CRISPR-systems, and lack quality control, annotation and assessment capabilities of the detected CRISPR loci. Our CRISPRloci server provides the first resource for the prediction and assessment of all possible CRISPR loci. The server integrates a series of advanced Machine Learning tools within a seamless web interface featuring: (i) prediction of all CRISPR arrays in the correct orientation; (ii) definition of CRISPR leaders for each locus; and (iii) annotation of cas genes and their unambiguous classification. As a result, CRISPRloci is able to accurately determine the CRISPR array and associated information, such as: the Cas subtypes; cassette boundaries; accuracy of the repeat structure, orientation and leader sequence; virus-host interactions; self-targeting; as well as the annotation of cas genes, all of which have been missing from existing tools. This annotation is presented in an interactive interface, making it easy for scientists to gain an overview of the CRISPR system in their organism of interest. Predictions are also rendered in GFF format, enabling in-depth genome browser inspection. In summary, CRISPRloci constitutes a full suite for CRISPR-Cas system characterization that offers annotation quality previously available only after manual inspection.

Casboundary: automated definition of integral Cas cassettes.

MOTIVATION: CRISPR-Cas are important systems found in most archaeal and many bacterial genomes, providing adaptive immunity against mobile genetic elements in prokaryotes. The CRISPR-Cas systems are encoded by a set of consecutive cas genes, here termed cassette. The identification of cassette boundaries is key for finding cassettes in CRISPR research field. This is often carried out by using Hidden Markov Models and manual annotation. In this article, we propose the first method able to automatically define the cassette boundaries. In addition, we present a Cas-type predictive model used by the method to assign each gene located in the region defined by a cassette's boundaries a Cas label from a set of pre-defined Cas types. Furthermore, the proposed method can detect potentially new cas genes and decompose a cassette into its modules. RESULTS: We evaluate the predictive performance of our proposed method on data collected from the two most recent CRISPR classification studies. In our experiments, we obtain an average similarity of 0.86 between the predicted and expected cassettes. Besides, we achieve F-scores above 0.9 for the classification of cas genes of known types and 0.73 for the unknown ones. Finally, we conduct two additional study cases, where we investigate the occurrence of potentially new cas genes and the occurrence of module exchange between different genomes. AVAILABILITY AND IMPLEMENTATION: https://github.com/BackofenLab/Casboundary. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Type IV CRISPR-Cas systems are highly diverse and involved in competition between plasmids.

CRISPR-Cas systems provide prokaryotes with adaptive immune functions against viruses and other genetic parasites. In contrast to all other types of CRISPR-Cas systems, type IV has remained largely overlooked. Here, we describe a previously uncharted diversity of type IV gene cassettes, primarily encoded by plasmid-like elements from diverse prokaryotic taxa. Remarkably, via a comprehensive analysis of their CRISPR spacer content, these systems were found to exhibit a strong bias towards the targeting of other plasmids. Our data indicate that the functions of type IV systems have diverged from those of other host-related CRISPR-Cas immune systems to adopt a role in mediating conflicts between plasmids. Furthermore, we find evidence for cross-talk between certain type IV and type I CRISPR-Cas systems that co-exist intracellularly, thus providing a simple answer to the enigmatic absence of type IV adaptation modules. Collectively, our results lead to the expansion and reclassification of type IV systems and provide novel insights into the biological function and evolution of these elusive systems.

Comprehensive search for accessory proteins encoded with archaeal and bacterial type III CRISPR-cas gene cassettes reveals 39 new cas gene families.

A study was undertaken to identify conserved proteins that are encoded adjacent to cas gene cassettes of Type III CRISPR-Cas (Clustered Regularly Interspaced Short Palindromic Repeats - CRISPR associated) interference modules. Type III modules have been shown to target and degrade dsDNA, ssDNA and ssRNA and are frequently intertwined with cofunctional accessory genes, including genes encoding CRISPR-associated Rossman Fold (CARF) domains. Using a comparative genomics approach, and defining a Type III association score accounting for coevolution and specificity of flanking genes, we identified and classified 39 new Type III associated gene families. Most archaeal and bacterial Type III modules were seen to be flanked by several accessory genes, around half of which did not encode CARF domains and remain of unknown function. Northern blotting and interference assays in Synechocystis confirmed that one particular non-CARF accessory protein family was involved in crRNA maturation. Non-CARF accessory genes were generally diverse, encoding nuclease, helicase, protease, ATPase, transporter and transmembrane domains with some encoding no known domains. We infer that additional families of non-CARF accessory proteins remain to be found. The method employed is scalable for potential application to metagenomic data once automated pipelines for annotation of CRISPR-Cas systems have been developed. All accessory genes found in this study are presented online in a readily accessible and searchable format for researchers to audit their model organism of choice: http://accessory.crispr.dk .

Toxin inhibition in C. crescentus VapBC1 is mediated by a flexible pseudo-palindromic protein motif and modulated by DNA binding.

Expression of bacterial type II toxin-antitoxin (TA) systems is regulated at the transcriptional level through direct binding of the antitoxin to pseudo-palindromic sequences on operator DNA. In this context, the toxin functions as a co-repressor by stimulating DNA binding through direct interaction with the antitoxin. Here, we determine crystal structures of the complete 90 kDa heterooctameric VapBC1 complex from Caulobacter crescentus CB15 both in isolation and bound to its cognate DNA operator sequence at 1.6 and 2.7 Å resolution, respectively. DNA binding is associated with a dramatic architectural rearrangement of conserved TA interactions in which C-terminal extended structures of the antitoxin VapB1 swap positions to interlock the complex in the DNA-bound state. We further show that a pseudo-palindromic protein sequence in the antitoxin is responsible for this interaction and required for binding and inactivation of the VapC1 toxin dimer. Sequence analysis of 4127 orthologous VapB sequences reveals that such palindromic protein sequences are widespread and unique to bacterial and archaeal VapB antitoxins suggesting a general principle governing regulation of VapBC TA systems. Finally, a structure of C-terminally truncated VapB1 bound to VapC1 reveals discrete states of the TA interaction that suggest a structural basis for toxin activation in vivo.

Characterizing leader sequences of CRISPR loci.

MOTIVATION: The CRISPR-Cas system is an adaptive immune system in many archaea and bacteria, which provides resistance against invading genetic elements. The first phase of CRISPR-Cas immunity is called adaptation, in which small DNA fragments are excised from genetic elements and are inserted into a CRISPR array generally adjacent to its so called leader sequence at one end of the array. It has been shown that transcription initiation and adaptation signals of the CRISPR array are located within the leader. However, apart from promoters, there is very little knowledge of sequence or structural motifs or their possible functions. Leader properties have mainly been characterized through transcriptional initiation data from single organisms but large-scale characterization of leaders has remained challenging due to their low level of sequence conservation. RESULTS: We developed a method to successfully detect leader sequences by focusing on the consensus repeat of the adjacent CRISPR array and weak upstream conservation signals. We applied our tool to the analysis of a comprehensive genomic database and identified several characteristic properties of leader sequences specific to archaea and bacteria, ranging from distinctive sizes to preferential indel localization. CRISPRleader provides a full annotation of the CRISPR array, its strand orientation as well as conserved core leader boundaries that can be uploaded to any genome browser. In addition, it outputs reader-friendly HTML pages for conserved leader clusters from our database. AVAILABILITY AND IMPLEMENTATION: CRISPRleader and multiple sequence alignments for all 195 leader clusters are available at http://www.bioinf.uni-freiburg.de/Software/CRISPRleader/ CONTACT: costa@informatik.uni-freiburg.de or backofen@informatik.uni-freiburg.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CRISPRstrand: predicting repeat orientations to determine the crRNA-encoding strand at CRISPR loci.

MOTIVATION: The discovery of CRISPR-Cas systems almost 20 years ago rapidly changed our perception of the bacterial and archaeal immune systems. CRISPR loci consist of several repetitive DNA sequences called repeats, inter-spaced by stretches of variable length sequences called spacers. This CRISPR array is transcribed and processed into multiple mature RNA species (crRNAs). A single crRNA is integrated into an interference complex, together with CRISPR-associated (Cas) proteins, to bind and degrade invading nucleic acids. Although existing bioinformatics tools can recognize CRISPR loci by their characteristic repeat-spacer architecture, they generally output CRISPR arrays of ambiguous orientation and thus do not determine the strand from which crRNAs are processed. Knowledge of the correct orientation is crucial for many tasks, including the classification of CRISPR conservation, the detection of leader regions, the identification of target sites (protospacers) on invading genetic elements and the characterization of protospacer-adjacent motifs. RESULTS: We present a fast and accurate tool to determine the crRNA-encoding strand at CRISPR loci by predicting the correct orientation of repeats based on an advanced machine learning approach. Both the repeat sequence and mutation information were encoded and processed by an efficient graph kernel to learn higher-order correlations. The model was trained and tested on curated data comprising >4500 CRISPRs and yielded a remarkable performance of 0.95 AUC ROC (area under the curve of the receiver operator characteristic). In addition, we show that accurate orientation information greatly improved detection of conserved repeat sequence families and structure motifs. We integrated CRISPRstrand predictions into our CRISPRmap web server of CRISPR conservation and updated the latter to version 2.0. AVAILABILITY: CRISPRmap and CRISPRstrand are available at http://rna.informatik.uni-freiburg.de/CRISPRmap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A novel interference mechanism by a type IIIB CRISPR-Cmr module in Sulfolobus.

Recent studies on CRISPR-based adaptive immune systems have revealed extensive structural and functional diversity of the interference complexes which often coexist intracellularly. The archaeon Sulfolobus islandicus REY15A encodes three interference modules, one of type IA and two of type IIIB. Earlier we showed that type IA activity eliminated plasmid vectors carrying matching protospacers with specific CCN PAM sequences. Here we demonstrate that interference-mediated by one type IIIB module Cmr-α, and a Csx1 protein, efficiently eliminated plasmid vectors carrying matching protospacers but lacking PAM motifs. Moreover, Cmr-α-mediated interference was dependent on directional transcription of the protospacer, in contrast to the transcription-independent activities of the type IA and type IIIA DNA interference. We infer that the interference mechanism involves transcription-dependent DNA targeting. A rationale is provided for the intracellular coexistence of the different interference systems in S. islandicus REY15A which cooperate functionally by sharing a single Cas6 protein for crRNA processing and utilize crRNA products from identical CRISPR spacers.

A novel single-tailed fusiform Sulfolobus virus STSV2 infecting model Sulfolobus species.

A newly isolated single-tailed fusiform virus, Sulfolobus tengchongensis spindle-shaped virus STSV2, from Hamazui, China, is characterised. It contains a double-stranded modified DNA genome of 76,107 bp and is enveloped by a lipid membrane structure. Virions exhibit a single coat protein that forms oligomers when isolated. STSV2 is related to the single-tailed fusiform virus STSV1 and, more distantly, to the two-tailed bicaudavirus ATV. The virus can be stably cultured over long periods in laboratory strains of Sulfolobus and no evidence was found for cell lysis under different stress conditions. Therefore, it constitutes an excellent model virus for archaeal virus-host studies.

CRISPR adaptive immune systems of Archaea.

CRISPR adaptive immune systems were analyzed for all available completed genomes of archaea, which included representatives of each of the main archaeal phyla. Initially, all proteins encoded within, and proximal to, CRISPR-cas loci were clustered and analyzed using a profile-profile approach. Then cas genes were assigned to gene cassettes and to functional modules for adaptation and interference. CRISPR systems were then classified primarily on the basis of their concatenated Cas protein sequences and gene synteny of the interference modules. With few exceptions, they could be assigned to the universal Type I or Type III systems. For Type I, subtypes I-A, I-B, and I-D dominate but the data support the division of subtype I-B into two subtypes, designated I-B and I-G. About 70% of the Type III systems fall into the universal subtypes III-A and III-B but the remainder, some of which are phyla-specific, diverge significantly in Cas protein sequences, and/or gene synteny, and they are classified separately. Furthermore, a few CRISPR systems that could not be assigned to Type I or Type III are categorized as variant systems. Criteria are presented for assigning newly sequenced archaeal CRISPR systems to the different subtypes. Several accessory proteins were identified that show a specific gene linkage, especially to Type III interference modules, and these may be cofunctional with the CRISPR systems. Evidence is presented for extensive exchange having occurred between adaptation and interference modules of different archaeal CRISPR systems, indicating the wide compatibility of the functionally diverse interference complexes with the relatively conserved adaptation modules.

Escherichia coli CRISPR arrays from early life fecal samples preferentially target prophages.

CRISPR-Cas systems are defense mechanisms against phages and other nucleic acids that invade bacteria and archaea. In Escherichia coli, it is generally accepted that CRISPR-Cas systems are inactive in laboratory conditions due to a transcriptional repressor. In natural isolates, it has been shown that CRISPR arrays remain stable over the years and that most spacer targets (protospacers) remain unknown. Here, we re-examine CRISPR arrays in natural E. coli isolates and investigate viral and bacterial genomes for spacer targets using a bioinformatics approach coupled to a unique biological dataset. We first sequenced the CRISPR1 array of 1769 E. coli isolates from the fecal samples of 639 children obtained during their first year of life. We built a network with edges between isolates that reflect the number of shared spacers. The isolates grouped into 34 modules. A search for matching spacers in bacterial genomes showed that E. coli spacers almost exclusively target prophages. While we found instances of self-targeting spacers, those involving a prophage and a spacer within the same bacterial genome were rare. The extensive search for matching spacers also expanded the library of known E. coli protospacers to 60%. Altogether, these results favor the concept that E. coli's CRISPR-Cas is an antiprophage system and highlight the importance of reconsidering the criteria use to deem CRISPR-Cas systems active.

Co-localization of antibiotic resistance genes is widespread in the infant gut microbiome and associates with an immature gut microbial composition.

BACKGROUND: In environmental bacteria, the selective advantage of antibiotic resistance genes (ARGs) can be increased through co-localization with genes such as other ARGs, biocide resistance genes, metal resistance genes, and virulence genes (VGs). The gut microbiome of infants has been shown to contain numerous ARGs, however, co-localization related to ARGs is unknown during early life despite frequent exposures to biocides and metals from an early age. RESULTS: We conducted a comprehensive analysis of genetic co-localization of resistance genes in a cohort of 662 Danish children and examined the association between such co-localization and environmental factors as well as gut microbial maturation. Our study showed that co-localization of ARGs with other resistance and virulence genes is common in the early gut microbiome and is associated with gut bacteria that are indicative of low maturity. Statistical models showed that co-localization occurred mainly in the phylum Proteobacteria independent of high ARG content and contig length. We evaluated the stochasticity of co-localization occurrence using enrichment scores. The most common forms of co-localization involved tetracycline and fluoroquinolone resistance genes, and, on plasmids, co-localization predominantly occurred in the form of class 1 integrons. Antibiotic use caused a short-term increase in mobile ARGs, while non-mobile ARGs showed no significant change. Finally, we found that a high abundance of VGs was associated with low gut microbial maturity and that VGs showed even higher potential for mobility than ARGs. CONCLUSIONS: We found that the phenomenon of co-localization between ARGs and other resistance and VGs was prevalent in the gut at the beginning of life. It reveals the diversity that sustains antibiotic resistance and therefore indirectly emphasizes the need to apply caution in the use of antimicrobial agents in clinical practice, animal husbandry, and daily life to mitigate the escalation of resistance. Video Abstract.

The infant gut virome is associated with preschool asthma risk independently of bacteria.

Bacteriophage (also known as phage) communities that inhabit the gut have a major effect on the structure and functioning of bacterial populations, but their roles and association with health and disease in early life remain unknown. Here, we analyze the gut virome of 647 children aged 1 year from the Copenhagen Prospective Studies on Asthma in Childhood2010 (COPSAC2010) mother-child cohort, all deeply phenotyped from birth and with longitudinally assessed asthma diagnoses. Specific temperate gut phage taxa were found to be associated with later development of asthma. In particular, the joint abundances of 19 caudoviral families were found to significantly contribute to this association. Combining the asthma-associated virome and bacteriome signatures had additive effects on asthma risk, implying an independent virome-asthma association. Moreover, the virome-associated asthma risk was modulated by the host TLR9 rs187084 gene variant, suggesting a direct interaction between phages and the host immune system. Further studies will elucidate whether phages, alongside bacteria and host genetics, can be used as preclinical biomarkers for asthma.

SMV1 virus-induced CRISPR spacer acquisition from the conjugative plasmid pMGB1 in Sulfolobus solfataricus P2.

Organisms of the crenarchaeal order Sulfolobales carry complex CRISPR (clustered regularly interspaced short palindromic repeats) adaptive immune systems. These systems are modular and show extensive structural and functional diversity, especially in their interference complexes. The primary targets are an exceptional range of diverse viruses, many of which propagate stably within cells and follow lytic life cycles without producing cell lysis. These properties are consistent with the difficulty of activating CRISPR spacer uptake in the laboratory, but appear to conflict with the high complexity and diversity of the CRISPR immune systems that are found among the Sulfolobales. In the present article, we re-examine the first successful induction of archaeal spacer acquisition in our laboratory that occurred exclusively for the conjugative plasmid pMGB1 in Sulfolobus solfataricus P2 that was co-infected with the virus SMV1 (Sulfolobus monocaudavirus 1). Although we reaffirm that protospacer selection is essentially a random process with respect to the pMGB1 genome, we identified single spacer sequences specific for each of CRISPR loci C, D and E that, exceptionally, occurred in many sequenced clones. Moreover, the same sequence was reproducibly acquired for a given locus in independent experiments, consistent with it being the first protospacer to be selected. There was also a small protospacer bias (1.6:1) to the antisense strand of protein genes. In addition, new experiments demonstrated that spacer acquisition in the previously inactive CRISPR locus A could be induced on freeze-thawing of the infected cells, suggesting that environmental stress can facilitate activation. Coincidentally with spacer acquisition, a mobile OrfB element was deleted from pMGB1, suggesting that interplay can occur between spacer acquisition and transposition.

Protospacer recognition motifs: mixed identities and functional diversity.

Protospacer adjacent motifs (PAMs) were originally characterized for CRISPR-Cas systems that were classified on the basis of their CRISPR repeat sequences. A few short 2-5 bp sequences were identified adjacent to one end of the protospacers. Experimental and bioinformatical results linked the motif to the excision of protospacers and their insertion into CRISPR loci. Subsequently, evidence accumulated from different virus- and plasmid-targeting assays, suggesting that these motifs were also recognized during DNA interference, at least for the recently classified type I and type II CRISPR-based systems. The two processes, spacer acquisition and protospacer interference, employ different molecular mechanisms, and there is increasing evidence to suggest that the sequence motifs that are recognized, while overlapping, are unlikely to be identical. In this article, we consider the properties of PAM sequences and summarize the evidence for their dual functional roles. It is proposed to use the terms protospacer associated motif (PAM) for the conserved DNA sequence and to employ spacer acqusition motif (SAM) and target interference motif (TIM), respectively, for acquisition and interference recognition sites.

The airway microbiota of neonates colonized with asthma-associated pathogenic bacteria.

Culture techniques have associated colonization with pathogenic bacteria in the airways of neonates with later risk of childhood asthma, whereas more recent studies utilizing sequencing techniques have shown the same phenomenon with specific anaerobic taxa. Here, we analyze nasopharyngeal swabs from 1 month neonates in the COPSAC2000 prospective birth cohort by 16S rRNA gene sequencing of the V3-V4 region in relation to asthma risk throughout childhood. Results are compared with previous culture results from hypopharyngeal aspirates from the same cohort and with hypopharyngeal sequencing data from the later COPSAC2010 cohort. Nasopharyngeal relative abundance values of Streptococcus pneumoniae, Haemophilus influenzae, and Moraxella catarrhalis are associated with the same species in the hypopharyngeal cultures. A combined pathogen score of these bacteria's abundance values is associated with persistent wheeze/asthma by age 7. No other taxa are associated. Compared to the hypopharyngeal aspirates from the COPSAC2010 cohort, the anaerobes Veillonella and Prevotella, which have previously been implicated in asthma development, are less commonly detected in the COPSAC2000 nasopharyngeal samples, but correlate with the pathogen score, hinting at latent community structures that bridge current and previous results. These findings have implications for future asthma prevention efforts.

Expanding known viral diversity in the healthy infant gut.

The gut microbiome is shaped through infancy and impacts the maturation of the immune system, thus protecting against chronic disease later in life. Phages, or viruses that infect bacteria, modulate bacterial growth by lysis and lysogeny, with the latter being especially prominent in the infant gut. Viral metagenomes (viromes) are difficult to analyse because they span uncharted viral diversity, lacking marker genes and standardized detection methods. Here we systematically resolved the viral diversity in faecal viromes from 647 1-year-olds belonging to Copenhagen Prospective Studies on Asthma in Childhood 2010, an unselected Danish cohort of healthy mother-child pairs. By assembly and curation we uncovered 10,000 viral species from 248 virus family-level clades (VFCs). Most (232 VFCs) were previously unknown, belonging to the Caudoviricetes viral class. Hosts were determined for 79% of phage using clustered regularly interspaced short palindromic repeat spacers within bacterial metagenomes from the same children. Typical Bacteroides-infecting crAssphages were outnumbered by undescribed phage families infecting Clostridiales and Bifidobacterium. Phage lifestyles were conserved at the viral family level, with 33 virulent and 118 temperate phage families. Virulent phages were more abundant, while temperate ones were more prevalent and diverse. Together, the viral families found in this study expand existing phage taxonomy and provide a resource aiding future infant gut virome research.