Viruses are the most abundant biological entities on earth: in most environments, from ocean water to the content of animal guts, the number of detected virus particles exceeds that of cells by one to two orders of magnitude ^[@R2]^. Among these viruses, more than 90% are tailed bacteriophages ^[@R1]^. More than 99% of the prokaryotic diversity in the biosphere is represented by bacteria and archaea that fail to grow in laboratory cultures, and accordingly, the great majority of the viruses are thought to infect these uncultivated microbes ^[@R1]^. Moreover, analysis of the human gut virome shows that most of the sequences, in contrast to the bacterial and archaeal sequences, have no matches in the current sequence databases, suggesting a vast virome consisting primarily of 'dark matter' ^[@R5]--[@R7]^.

The crAssphage is the utmost manifestation of this trend. The complete crAssphage (after Cross Assembly) genome was assembled by joining contigs obtained from several human fecal viral metagenomes as a circular double-stranded (ds) DNA molecule of approximately 97 kilobase (kb)^[@R4]^; the circular genome map apparently results from the terminal redundancy and/or circular permutation. The crAssphage is extremely abundant, accounting for up to 90% of the reads in the virus-like particle-enriched fraction of the gut metagenome and about 22% of the reads in the total metagenome. Numerous reads matching the crAssphage genome have been identified in numerous gut metagenomes collected in diverse geographic locations, indicating that crAssphage is not only the most abundant virus in the human gut microbiome but also a (nearly) ubiquitous one ^[@R4],[@R8],[@R9]^. Read co-occurrence analysis points to bacteria of the phylum Bacteroidetes as the host(s) of the crAssphage ^[@R4],[@R10]^. This assignment is compatible with the presence, in the crAssphage genome, of a protein containing carbohydrate-binding domains (BACON domains) that is highly similar to a homologous protein from *Bacteroides* and with partial matches between two crAssphage sequences and CRISPR spacers from two species of *Bacteroides* ^[@R4]^. Members of the Bacteroidetes dominate the gut microbiome but most of these bacteria so far have not been grown in culture ^[@R11],[@R12]^. Thus, it is hardly surprising that the most abundant -- but never isolated -- phage from this environment appears to be a parasite of Bacteroidetes. Analysis of the protein sequences encoded in the crAssphage genome failed to identify specific relationships with other bacteriophages ^[@R4]^. Several proteins implicated in phage genome replication have been identified including a family B DNA polymerase (DNAP), a primase and a flavin-dependent thymidylate synthase, but neither the major capsid protein nor other structural and morphogenetic proteins were detected. In an attempt to clarify the provenance of this most abundant but enigmatic human-associated virus, we re-analyzed the crAssphage genome using the most sensitive available methods for protein sequence analysis and taking advantage of the database growth since the time of crAssphage discovery. The result is the identification of a previously unknown, expansive bacteriophage family that appears to be associated with diverse members of Bacteroidetes and for which we now recognize the structural, replication and expression gene modules.

The sequences of the crAssphage proteins were compared, using PSI-BLAST, to the non-redundant protein sequence database (nr) and the Whole Genome Shotgun (WGS) databases (NCBI, NIH, Bethesda) containing microbial genomic and metagenomic sequences. Sequences with significant similarity to crAssphage proteins were detected in four genomes of previously identified bacteriophages and numerous contigs assigned to bacterial genomes (possibly, prophages) and metagenomic contigs. These sequences were highly diverse, and most were not closely related (despite the statistical significance of the detected similarity) to the crAssphage proteins ([Supplementary Table 1](#SD1){ref-type="supplementary-material"}). Thus, crAssphage relatives identified here might comprise a previously unidentified, large, diverse family of bacteriophages (henceforth crAss-like family), or potentially, even two or more families. Altogether, we identified several hundred putative representatives of the crAss-like phage family ([Supplementary Figure 1](#SD1){ref-type="supplementary-material"}); 37 diverse representatives, for which (nearly) complete genomes were available, were selected for in-depth analysis ([Supplementary Table 1](#SD1){ref-type="supplementary-material"}). We then constructed multiple alignments of the crAssphage proteins and their homologs, and used these alignments as queries for profile-profile searches against a comprehensive collection of protein families using the HHPred software, one of the most sensitive current methods for protein sequence analysis (see Methods for details). We identified a block of 5 genes that appear to comprise the structural module of the emerging family of bacteriophages. These genes encode a predicted major capsid protein (MCP) of the HK97 fold, portal protein, large terminase subunit, and two uncharacterized proteins that are conserved throughout the crAss-like family and, given the consistent adjacency to the MCP, could be components of the virion or the morphogenetic machinery ([Table 1](#T1){ref-type="table"}; [Figure 1](#F1){ref-type="fig"}; [Supplementary Table 2](#SD1){ref-type="supplementary-material"} and [Supplementary Note 1](#SD1){ref-type="supplementary-material"}).

Despite the low sequence conservation, even within the family, and remote similarity to proteins from other phages, the gene order in the capsid module of the crAss-like family is nearly invariant ([Figure 1](#F1){ref-type="fig"}), suggesting congruent evolution of these genes. A concatenated alignment of all 5 genes was used to construct a phylogenetic tree of the crAss-like family which includes a strongly supported clade of crAssphage relatives and several other distinct groups of (predicted) bacteriophages ([Figure 1](#F1){ref-type="fig"}).. Three of these groups included previously identified phages, namely *Azobacteroides* phage ProJPt-Bp ^[@R13]^, *Flavobacterium psychrophilum* phage Fpv3 and *Cellulophaga* phage phi14:2, widespread, although apparently not highly abundant phages in the oceans ^[@R14]^. The same group that included the *Cellulophaga* phage also contained the genome of the IAS (immunodeficiency-associated stool) virus that is highly abundant in gut viromes of HIV-infected individuals^[@R15]^. Most of the other members of the crAss-like family are unassigned metagenomic sequences but several come from bacterial genome assemblies and might represent prophages. All experimentally characterized crAss-like phages are associated with Bacteroidetes. Among the other sequences included in the family, several are assigned to other bacteria, in particular, *Chlamydia trachomatis*, as well as several members of the recently identified candidate phyla radiation (CPR) ^[@R16]^. However, in the phylogenetic trees of the predicted MCP, these sequences are embedded within groups consisting of sequences associated with Bacteroidetes ([Figure 1](#F1){ref-type="fig"} and [Supplementary Figure 1](#SD1){ref-type="supplementary-material"}), and none of these contigs contained genes that could be linked to the host bacteria. Thus, the available data appear compatible with exclusive association between the crAss-like phages and Bacteroidetes.

In the linear genome maps shown in [Figure 2](#F2){ref-type="fig"} (see [Supplementary Table 1](#SD1){ref-type="supplementary-material"} and [Supplementary Note 1](#SD1){ref-type="supplementary-material"} for the crAssphage and IAS phage gene annotations), the capsid structural module occupies about 10 kb near one end of the crAss-like phage genomes. Downstream of this module, are the genes encoding predicted tail proteins and two proteins homologous to bacterial Integration Host Factor (IHF) that is essential for chromatin packaging in bacteria and some phages ^[@R17]^ ([Figure 2](#F2){ref-type="fig"}; [Table 1](#T1){ref-type="table"}; [Supplementary Table 2](#SD1){ref-type="supplementary-material"}). The two most conserved tail proteins encoded by the crAss-like group are homologous to tail components gp4 (tubular tail protein) and gp10 (tail stabilization protein) of bacteriophage P22 ^[@R18]^. In the same putative operon, crAssphage also encodes a homolog of the tail needle protein gp26 of phage P22. These three proteins are sufficient for the formation of a short tail similar to that of bacteriophage P22 ^[@R19]^.The gp10 homologs of crAss-like phages are large (\>1400 aa), apparently multidomain proteins in which the gp10-homologous region accounts only for about 150 aa. The additional domains of this protein could be involved in host recognition similarly to the tail spike protein of P22-like phages ^[@R20]^. Some crAss-like phages encode additional, auxiliary proteins in the tail module, e.g. an IAS protein homologous to the tail-associated lysozyme gp13 of short-tailed phage phi29 ^[@R21]^. Thus, the crAss-like phages can be predicted to possess short, stubby tails, a hallmark of the family *Podoviridae*. One of the isolated crAss-like phages, *Cellulophaga* phage phi14:2, is indeed a typical podophage ^[@R14]^. Unlike phages with long tails (families *Myoviridae* and *Syphoviridae*), P22-like phages do not encode maturation proteases ^[@R22]^, consistent with the apparent absence of such a protease among gene products of the crAss-like phages. However, in the midst of the genes for predicted tail components, some of the crAss-like phages, including the crAssphage group, encode a predicted Zn-dependent protease ([Figure 2](#F2){ref-type="fig"}) that might be involved in processing of the tail and/or capsid proteins. In crAssphage group genomes, the protease gene is embedded within a block of genes encoding putative additional tail components which are highly similar to homologs from uncharacterized prophages integrated in Bacteroides genomes that are otherwise unrelated to crAss-like phages ([Figure 2](#F2){ref-type="fig"}; also see [Supplementary Table 2](#SD1){ref-type="supplementary-material"}). Thus, evolution of the crAssphage group apparently involved relatively recent recombination with an unrelated (pro)phage from the same host(s).

The lytic replication module genes occupy about 30 kb on the opposite end of the genome from the capsid module and are transcribed towards the middle of genome ([Figure 2](#F2){ref-type="fig"}). This module encodes a versatile suite of proteins implicated in DNA replication and repair, and shows a patchy gene distribution, without a single universally conserved gene, and a much greater variability within the crAss-like family than the structural modules ([Table 1](#T1){ref-type="table"}; [Figure 3](#F3){ref-type="fig"}). The most conserved replicative gene is a predicted DnaG family primase; many crAss-like family members also encode a putative superfamily 2 (SNF2 family) helicase implicated in DNA replication (this helicase is inactivated in the crAssphage subfamily as indicated by mutiple amino acid replacements in the catalytic sites), an ATP-dependent DNA ligase, a uracyl-DNA glycosylase (UDG), a flavin-dependent thymidylate synthase (ThyX), and one or two diverged single-stranded DNA-binding proteins (SSB) ([Table 1](#T1){ref-type="table"}; [Figure 3](#F3){ref-type="fig"}). The crAss-like family viruses encode one of the two distinct DNA polymerases (DNAPs): the crAssphage clade has a family B DNAP, whereas the other family members have either a family B or a family A DNAP, or no DNAP at all ([Table 1](#T1){ref-type="table"}; [Figure 3](#F3){ref-type="fig"}). Evolutionary reconstruction suggests that the family A DNAP is ancestral in crAss-like phages and was lost or replaced with the family B DNAP on several occasions ([Supplementary Figure 2](#SD1){ref-type="supplementary-material"}). The only near universal replicative protein, DNA primase, appears to be monophyletic in the crAss-like family and forms a strongly supported clade with the primases of Bacteroidetes ([Supplementary Figure 2](#SD1){ref-type="supplementary-material"}). This is the only conserved crAss-like family gene that shows a deep connection to Bacteroidetes, indicating that a founder crAss-like phage acquired this gene from a Bacteroidetes host and implying that the virus-host link is evolutionarily ancient. An unusual evolutionary connection was detected for the phage ligase: in the phylogenetic tree, the crAss-like family ligases clustered with those of eukaryotic giant dsDNA viruses suggesting that these viruses acquired the ligase from crAss-like phages ([Supplementary Figure 2](#SD1){ref-type="supplementary-material"}). Portions of the replicative gene block of the IAS virus and its closest relatives show high similarity to homologs from putative prophages of Bacteroidetes, suggesting another recombination event ([Figure 2](#F2){ref-type="fig"} and see [Supplementary Table 2](#SD1){ref-type="supplementary-material"}).

The structural and replicative gene blocks of the crAss-like phages are separated by an array of uncharacterized genes that are transcribed in the same direction as the structural genes and are universally conserved across the crAss-like family ([Figure 2](#F2){ref-type="fig"} and [Table 1](#T1){ref-type="table"}). Several of these genes encode giant proteins, up to 6000 amino acids in size. An HHpred search initiated with the multiple alignment of the homologs of one of these large proteins (crAssphage gene 46 product) identified a small region of similarity with the β′-subunit of the bacterial RNA polymerase (RNAP), which contained the signature catalytic loop with the metal-binding DxDxD motif ^[@R23]^ ([Supplementary Figure 3](#SD1){ref-type="supplementary-material"}). Detailed sequence analysis resulted in the identification of two additional conserved motifs typical of the RNAP β subunit ^[@R23]^, suggesting that the two subunits are fused in this protein ([Figure 4](#F4){ref-type="fig"} and [Supplementary Figure 3](#SD1){ref-type="supplementary-material"}). Although the similarity between these crAss-like family protein sequences and the large RNAP subunits was limited, the strict conservation of several predicted key motifs that comprise the RNAP catalytic site across the crAss-like family, the fusion of the putative homologs of two RNAP subunits in a single large protein compatible with the combined lengths of the β and β′ RNAP subunits, and the compatibility of the predicted secondary structure elements with the RNAP core structure ([Supplementary Figure 3](#SD1){ref-type="supplementary-material"}) strongly suggest that the crAss-like family phages encode an active RNAP.

The putative β--β′ RNAP fusion protein contains another large region of similarity with other phages, which in some of them resides in a separate protein (e.g. gene_53 product of Azobacteroides phage ProJPt-Bp1; [Figure 4](#F4){ref-type="fig"}). Another protein conserved in most of the crAss-like phages (crAssphage gene 47) is typically encoded next to or is fused to the β--β′ RNAP ([Figure 4](#F4){ref-type="fig"}). Most likely, all three proteins are functionally linked and form a multisubunit RNAP; although no homologs of the gene 47 product were detected, the size and association with the RNAP subunits suggest that this could be a highly diverged α subunit. In most crAss-like phages that encode fused subunits of the predicted RNAP, a putative zincin-like protease domain is encoded in the vicinity, e.g. within the gene 45--46 fusion product of the crAssphage, whereas most of the phages that encode the RNAP subunits in separate genes lack the predicted protease ([Figure 4](#F4){ref-type="fig"} and [Supplementary Figure 3](#SD1){ref-type="supplementary-material"}). Thus, the fused RNAP subunits might be cleaved by the Zn-dependent protease to produce the mature proteins. Fusions of zincin family proteases are typical of different multidomain phage proteins ^[@R24]^, which is compatible with the RNAP cleavage hypothesis. Nonetheless, it cannot be ruled out that the fusion protein is the active form of the phage RNAP. Only a few groups of phages encode their own multisubunit RNAPs, including *Lactococcus* phage 1706, *Rhodococcus* phage ReqiPepy6, *Bacillus* phage SPbeta, *Thermus* phages P74 and P23, and giant phages of *Pseudomonas* ^[@R25]--[@R28]^. In most of these phages, β and β′ subunits are fused ^[@R25],[@R29]^, whereas in *Pseudomonas* group phages, each subunit is split into two proteins ^[@R28]^. The phage RNAPs belong to diverged families that can be considered signatures of each respective phage group ^[@R26],[@R28]^. The crAss-like family RNAPs are even more extremely divergent than those of other phages. To our knowledge, processing of RNAP polyproteins by a dedicated protease so far has not been identified in viruses or cellular organisms. Phage RNAPs transcribe either early (replicative) ^[@R28]^ or late (structural) ^[@R27]^ genes; in the former case, the RNAP presumably is packaged into the virion. We attempted to identify promoters of early and late genes of crAss-like phages by searching the sequences upstream of the genes for potential conserved nucleotide motif, but no such motifs were detected.

The discovery of crAssphage, the most abundant virus in the human gut virome, appeared particularly striking because the genome was *terra incognita*, with few homologs detected in other viruses or bacteria, and the virion proteins not identified ^[@R4]^. The present analysis changes this by showing that crAssphage belongs to an expansive phage family that is only distantly related to other known phages and has unusual predicted features, in particular, a previously unknown putative mechanism of RNAP maturation via polyprotein processing. The MCP of these phages, a distinct form of the HK97 class of icosahedral capsid proteins ^[@R30]^, is now confidently predicted and amenable for experiments aimed at direct identification and characterization of the phage. Altogether, homologs with characterized functions have been detected for 53% (48 out of 91) of MetaGeneMark-predicted crAssphage proteins (compared to 26% in the original analysis and 14% in the current RefSeq annotation; [Supplementary Table 2](#SD1){ref-type="supplementary-material"}).

The crAss-like family includes at least one additional phage, IAS, that can reach high abundance in the human gut ^[@R15]^. Generally, the crAss-like family appears to be abundant and widespread in diverse habitats, both animal-associated and environmental. Various bacteria of the phylum *Bacteroidetes* appear to be the primary hosts of crAss-like phages as indicated by the presence of several genes apparently derived from these bacteria including the DNA primase that is ancestral in the family and the BACON domain protein implicated in phage adhesion to mucus that could increase the frequency of the encounters with the host bacteria ^[@R31]^. This virus-host association is supported by CRISPR spacer analysis; in addition to the previously reported imperfect matches to *Bacteroides* genomes ^[@R4]^, we detected perfect matches of crAssphage sequences to two spacers from *Porhyromonas* sp. ([Supplementary Note 2](#SD1){ref-type="supplementary-material"}). It seems likely that crAssphage has a broad host range among the Bacteroidetes, which could contribute to the (near) ubiquity of this phage in humans. For the IAS virus, spacers with partial matches were detected in CRISPR arrays of *Prevotella* ([Supplementary Note 2](#SD1){ref-type="supplementary-material"}), again indicating a *Bacteroidetes* host. Although some of the crAss-like family sequences identified here are assigned to genomes of other bacteria, these assignments could be erroneous (see above), suggesting that the crAss-like virus family is Bacteroidetes-specific. Our results indicate that some of the crAss-like phages are temperate and lysogenize their hosts by integrating into their genomes with the aid of phage-encoded tyrosine integrases. The temperate life style increases the opportunities for recombination with other phages and could account for the \\presence of regions with high similarity to otherwise unrelated Bacteroidetes prophages.

Our analysis of the predicted tail proteins indicates that the crAss-like phages possess short, podovirus-like tails. Thus, under the classical morphology-guided classification scheme, these phages would be classified into the family *Podoviridae*. However, given that phage taxonomy is moving towards sequenced-based approaches ^[@R32]^, crAss-like phages are likely to become a family within the order *Caudovirales*. The general lesson from this study is that, with the current proliferation of genomic and especially metagenomic sequence databases and advances in database search approaches, any discovered abundant virus or microbe is likely to become a prototype of a previously undetected, often highly diverse group of organisms. Such advanced analyses can guide experimental study of viruses and microbes that are currently known only through genomic sequences but could be major players in the microbiota.

Methods {#S1}
=======

The search for crAssphage structural proteins was performed as follows. The sequences of the crAssphage proteins were first used as queries in a PSI-BLAST search ^[@R33]^ of the NCBI non-redundant (nr) database. Proteins that produced no statistically significant hits to nr proteins with predicted functions were considered as candidates for crAssphage structural proteins. For each of these proteins, homologs from both nr and and environmental (env_nr) protein sequence databases were collected (using PSI-BLAST with default parameters until convergence or for 5--6 iterations). The homologs detected in this search were aligned with the query crAssphage protein, and the alignments were used as queries for HHPred ^[@R34]^ searches. Three of these queries produced statistically significant hits to phage structural proteins (terminase large subunit, portal, and HK97 family MCP).

The sequences of the predicted crAssphage MCP and its homologs ("initial MCP set") were used as queries for translating blast (TBLASTN) searches against the wgs and nr nucleotide databases. Nucleotide sequences of the hits were retrieved and translated in 6 frames, and the sequences of the putative MCP homologs were validated by comparison to the initial MCP set. The extracted protein sequences of the putative MCP homologs (either complete proteins or fragments longer than 150 aa) were clustered with blastclust at 90% identity, aligned with MUSCLE ^[@R35]^, and used for the phylogenetic reconstruction shown in [Supplementary Figure 1](#SD1){ref-type="supplementary-material"}.

In addition, other conserved crAssphage proteins were searched against GenBank databases by protein blast (BLASTP) and TBLASTN ^[@R33]^. These searches led to the identification of several complete and partial crAss-like viral genomes in the nr database as well as several hundred crAss-like contigs in the wgs databases ([Supplementary Figure 1](#SD1){ref-type="supplementary-material"}; see <ftp://ftp.ncbi.nih.gov/pub/yutinn/crassphage_2017/> for additional data).

For in-depth analysis, 37 representative genomes (contigs) were selected among the crAss-like family members identified in the nr and wgs databases. In some cases (namely, CVNZ01000019ext, CVNZ0100007ext, contig0001, contig0002, and contig0005), representative contigs were obtained by additional assembly using the Geneious software. The original contigs used for these assemblies are listed in [Supplementary Table 1](#SD1){ref-type="supplementary-material"}.

All representative genomes were translated using MetaGeneMark ^[@R36]^. The set of open reading frames produced by this translation for crAssphage was virtually identical to the original set that was produced using the Glimmer software ^[@R37]^; all genes detected by MetaGeneMark, for which a function was predicted, were represented also in the original Glimmer translation ^[@R4]^ ([Supplementary Table 2](#SD1){ref-type="supplementary-material"}). Homologous protein sequences were aligned using MUSCLE ^[@R35]^. In order to identify putative homologs outside the crAss-like family, the alignments of all conserved proteins were used as queries for PSI-BLAST ^[@R33]^ and HHPred ^[@R34]^ searches.

The analysis reported here differed from the original analysis of the crAssphage genome in the following respects: i) expanded databases of genomic and metagenomics sequences were searched, ii) the updated versions of the PSI-BLAST and HHPred software were used, iii) multiple alignments of crAssphage proteins and their homologs, rather than individual crAssphage protein sequences, were used to initiate the searches. Furthermore, no pre-set cut-offs were used in database searches, and all search results and alignments were examined individually for conservation of diagnostic sequence motifs. Together, these amendments to the sequence analysis protocol yielded substantially enhanced search results.

For phylogenetic reconstruction of crAss-like family MCP, PolA, PolB, primase, and ligase ([Supplementary Figures 1 and 2](#SD1){ref-type="supplementary-material"}), gapped columns (more than 30% of gaps) and columns with low information content were removed from the alignments ^[@R38]^; filtered alignments were used for tree reconstructions using the FastTree program ^[@R39]^. For the phylogenetic tree shown on [Figure 1](#F1){ref-type="fig"}, the alignments of five conserved proteins of the capsid module were concatenated and used for phylogenetic analysis with PhyML program (<http://www.atgc-montpellier.fr/phyml-sms/>) ^[@R40]^. The best model identified by PhyML was LG +G+I+F (LG substitution model, gamma distributed site rates with gamma shape parameter estimated from the alignment; fraction of invariable sites estimated from the alignment; and empirical equilibrium frequencies).

Search for nucleotide sequence motifs was performed using the MEME ^[@R41]^ and Gibbs Centroid Sampler ^[@R42]^ programs.

Data availability {#S2}
-----------------

All the data used for the analysis reported in this work are publicly available through GenBank. Genbank Accession Numbers for the representative set of contigs that have been analyzed in detail are given in the [Supplementary Dataset 1](#SD1){ref-type="supplementary-material"}, and the complete set of Accession Numbers for all crAss-like contigs is available at <ftp://ftp.ncbi.nih.gov/pub/yutinn/crassphage_2017/MCP_containing_contig_sources.xlsx>. The annotation of the crAssphage and IAS virus genes is given in the [Supplementary Dataset 1](#SD1){ref-type="supplementary-material"}. Further [supporting information](#SD1){ref-type="supplementary-material"} is available in [Supplementary Note 1](#SD1){ref-type="supplementary-material"} and at <ftp://ftp.ncbi.nih.gov/pub/yutinn/crassphage_2017/>.

Supplementary Material {#S3}
======================

###### 

**Supplementary Table 1**. The selected representative set of crAss-like family members analyzed in detail in this work; conserved genes in the crAss-like phage family (an extended version of [Table 1](#T1){ref-type="table"}); BLAST scores of conserved crAss-like family proteins

**Supplementary Table 2**. Annotation of the crAssphage and IAS phage genes and comparison of the MetaGeneMark annotation with the current RefSeq (Glimmer) crAssphage annotations

**Supplementary Figure 1**. Phylogenetic tree of the MCP for all identified members of the crAss-like family

**Supplementary Figure 2**. Phylogenetic trees for crAss-like family PolA, PolB, primase, and ligase

**Supplementary Figure 3**. Multiple alignment of conserved regions of β and β′ RNA polymerase subunits of the crAss-like family

**Supplementary Note 1**. Supporting evidence for remote sequence similarities used for crAssphage and IAS phage gene annotation

**Supplementary Note 2.** CRISPR spacers matching crAss-like family genomes

The authors thank Yuri I. Wolf and Sergey Shmakov for invaluable technical help, and Koonin group members for helpful discussion. NY, KSM, ABG, and EVK are supported by intramural funds of the US Department of Health and Human Services (to the National Library of Medicine).

**Authors' contributions**

EVK conceived of the study; NY, KSM and MK performed research; NY, KSM, ABG, MK, AS, RE and EVK analyzed the data; EVK wrote the manuscript that was read, edited and approved by all authors.

**Competing interests**

The authors declare no competing financial interests.

![Architecture and evolution of the capsid gene module of the crAss-like phage family\
The phylogenetic tree was constructed from concatenated multiple alignments of the 5 proteins of the capsid module. The genomic maps of the capsid gene block are shown for each branch. The 5 genes of the capsid module are color-coded, and uncharacterized adjacent genes are shown by empty block arrows. The colors of labels and branches indicate host or metagenome source: Red, human gut or fecal metagenome; Green, termite gut metagenome; Purple, terrestrial/groundwater; Brown, Soda lake (hypersaline brine); Turquoise, Marine sediment; Orange, populus root microbiome; Black, *Flavobacterium psychrophilum* 950106-1/1 (fish pathogen). Tree branch colors indicate the DNA polymerase family represented in the respective genome (contig): Purple, family B DNAP; Red, family A DNAP; Green, no DNAP; Black, unknown (incomplete genomes). Support values were obtained using 100 bootstrap replications; values greater than 50% are shown. The scale bar represents the number of amino acid substitutions per site. No outgroup was included due to the low (or absent) similarity between the crAss-like family protein to homologs from other phages.](nihms910991f1){#F1}

![Whole genome maps of crAssphage and IAS virus, the two members of the crAss-like family that are abundant in the human gut virome\
Conserved crAss-like family genes are color-coded. Dashed boxes highlight capsid, tail, replication, and transcription gene blocks. Genome regions encoding proteins with strong similarity to Bacteroidetes are shaded in pale green. Gene numbers are according to the crAssphage and IAS virus MetaGeneMark translations. Abbreviations: ssb, single-stranded DNA-binding protein; SF1, SF1 helicase; UDG, uracyl-DNA glycosylase; PolB, DNA polymerase family B; SF2, SNF2-family helicase; RecT, phage RecT recombinase; primase, DnaG family primase; ligase, ATP-dependent DNA ligase; dNK, deoxynucleoside monophosphate kinase; ThyX, flavin-dependent thymidylate synthase; Gp157, Siphovirus Gp157; dUTP, dUTPase; N4_gp49, phage protein of N4_gp49/Sf6_gp66 family; RepL, plasmid replication initiation protein RepL; IHF, integration host factor IHF subunit; PD-(D/E)XK, PD-(D/E)XK family nuclease; Rep_Org, putative replisome organizer protein; DnaB, DnaB replicative DNA helicase; AAA, AAA domain ATPase; rIIA, rIIA-like protector protein; rIIB, rIIB-like protector protein; NRDD, anaerobic ribonucleoside-triphosphate reductase; RNR, anaerobic ribonucleoside-triphosphate reductase activating protein; PolA, DNA polymerase family A; DprA, DNA processing protein DprA. For further details on the annotation, see [Supplementary Table 2](#SD1){ref-type="supplementary-material"}.](nihms910991f2){#F2}

###### The replicative gene module of the crAss-like phage family

A. The crAssphage group

B. The rest of the crAss-like family

Homologous genes are marked by the same colors and labels. Genes with no predicted function are numbered according to crAssphage translation, OBV-13 virus translation (suffix 'b'), Cellulophaga phage phi14:2 translation (suffix 'c'), or IAS virus translation (suffix 'i'). Abbreviations:: RNRm, class II ribonucleotide reductase; RNRa, ribonucleoside reductase alpha chain; RNRb, ribonucleoside reductase beta chain; GGCT, Gamma-glutamyl cyclotransferase; Gn_AT, glucosamine-fructose-6-phosphate aminotransferase; NTP-PPase, nucleoside triphosphate pyrophosphohydrolase. The rest of the abbreviations are as in [Figure 2](#F2){ref-type="fig"}.

![](nihms910991f3a)

![](nihms910991f3b)

![The genome expression gene module of the crass-like phage family\
The predicted RNAP subunits as well as the RNAP and protease motifs are color-coded as shown at the bottom of the figure. The PD-DxK nucleases are most likely encoded in Group I introns.](nihms910991f4){#F4}

###### 

Conserved genes in the crAss-like phage family

  gene                                                                   crAssphage gene number   IAS virus gene number   crAssphage   CDZK01015469   CEAR01029167   LSPZ01000006   LAZR01000126   CESO01030555   BCSF01000013   LSPY01000004   Azobact_ph_AP017903   Chlamydia_CVNZ01000019ex   contig0001   FUFK010039141   Woesebacteria_MGFQ01000035   Chitinophaga_FOJF01000001   Flavobact_ph_NC_031904   MDTC01014143   CEUT01009082/   Cellulophaga_ph_NC_021806   IAS_virus_KJ003983
  ---------------------------------------------------------------------- ------------------------ ----------------------- ------------ -------------- -------------- -------------- -------------- -------------- -------------- -------------- --------------------- -------------------------- ------------ --------------- ---------------------------- --------------------------- ------------------------ -------------- --------------- --------------------------- --------------------
  ***Core genes of crAss-like family***[a](#TFN1){ref-type="table-fn"}                                                                                                                                                                                                                                                                                                                                                                                                     
  *Terminase*                                                            *79*                     *15*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *Portal*                                                               *78*                     *16*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *putative structural protein*                                          *77*                     *32*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *MCP*                                                                  *76*                     *33*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *putative structural protein*                                          *75*                     *34*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *putative structural protein*                                          *74*                     *35,36*                 **Y**        **Y**          **Y**          **Y**          Y              Y              Y                             Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *putative structural protein*                                          *73*                     *-*                     **Y**        **Y**          **Y**          **Y**                         Y              Y              Y                                                               Y            Y               Y                            Y                           Y                        Y              Y               Y                           
  *IHF subunit*                                                          *54*                     *40*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                                                Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *Tail tubular protein (P22 gp4-like)*                                  *52*                     *34*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                                                Y            Y                                            Y                           Y                        Y              Y               Y                           Y
  *Tail stabilization protein (P22 gp10-like)*                           *51*                     *44*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                                                Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *Putative RNAP subunit*                                                *47*                     *47*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *β--β′ RNAP*                                                           *46N*                    *49*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *Putative RNAP subunit*                                                *46C*                    *48*                    **Y**        **Y**          **Y**          **Y**          ?              ?              Y              ?              Y                     Y                          Y            Y               Y                            Y                           Y                        Y              Y               Y                           Y
  *DnaG family primase*                                                  *22*                     *80_79*                 **Y**        **Y**          **Y**          **Y**          Y              Y              Y              Y                                    Y                                                                                    Y                           Y                        Y              Y               Y                           Y
  ***Conserved genes in the crAssphage group***                                                                                                                                                                                                                                                                                                                                                                                                                            
  *Uncharacterized protein*                                              *86*                     *13*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y                                                                              Y            Y                                                                        Y                        Y              Y               Y                           Y
  *Uncharacterized protein*                                              *49*                     *45*                    **Y**        **Y**          **Y**          **Y**                         Y              Y                                                                                                                                        Y                           Y                                                       Y                           Y
  *Uncharacterized protein*                                              *34*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                       Y            Y                                                                                                                Y                                           
  *Siphovirus Gp157 homolog*                                             *32*                     *-*                     **Y**        **Y**          **Y**          **Y**          Y                                                                                                                                                                                                                                                                                      
  *Uncharacterized protein*                                              *23*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                                                                                                                                                                                                 
  *ssb*                                                                  *21*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                                                                                                                                                                                                 
  *RecT*                                                                 *19*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                                                                                                                                                                                                 
  *SNF2 helicase*[b](#TFN2){ref-type="table-fn"}                         *18*                     *87*                    **Y**        **Y**          **Y**          **Y**          Y              Y              Y                                                                                                                                                                    Y                        Y              Y               Y                           Y
  *PolB*                                                                 *17*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                            Y                                                                                    Y                           Y                                                       Y                           
  *Uncharacterized protein*                                              *20*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                                                                                                                                                                                                 
  *UDG*                                                                  *16*                     *75*                    **Y**        **Y**          **Y**          **Y**          Y                             Y                                                                                                                                        Y                           Y                                       Y               Y                           Y
  *SF1 helicase*                                                         *15*                     *-*                     **Y**        **Y**          **Y**          **Y**                                                                                                                                                                                                                                                                                                 

The core genes include those that are represented in all 5 major branches of the family according to the MCP tree ([Figure 1](#F1){ref-type="fig"} and [Supplementary figure S1](#SD1){ref-type="supplementary-material"}). The crAssphage group is highlighted by bold type.

The predicted SNF2 family helicase is inactivated in the crAssphage group but active in the other phages.

Other crAss-like family members encode distant members of the PD-(D/E)XK nuclease family.
