Structure and function of virion RNA polymerase of a crAss-like phage

Genome Analysis of Epsilon CrAss-like Phages

CrAss-like phages play an important role in maintaining ecological balance in the human intestinal microbiome. However, their genetic diversity and lifestyle are still insufficiently studied. In this study, a novel CrAssE-Sib phage genome belonging to the epsilon crAss-like phage genomes was found. Comparative analysis indicated that epsilon crAss-like phages are divided into two putative genera, which were proposed to be named Epsilonunovirus and Epsilonduovirus; CrAssE-Sib belongs to the former. The crAssE-Sib genome contains a diversity-generating retroelement (DGR) cassette with all essential elements, including the reverse transcriptase (RT) and receptor binding protein (RBP) genes. However, this RT contains the GxxxSP motif in its fourth domain instead of the usual GxxxSQ motif found in all known phage and bacterial DGRs. RBP encoded by CrAssE-Sib and other Epsilonunoviruses has an unusual structure, and no similar phage proteins were found. In addition, crAssE-Sib and other Epsilonunoviruses encode conserved prophage repressor and anti-repressors that could be involved in lysogenic-to-lytic cycle switches. Notably, DNA primase sequences of epsilon crAss-like phages are not included in the monophyletic group formed by the DNA primases of all other crAss-like phages. Therefore, epsilon crAss-like phage substantially differ from other crAss-like phages, indicating the need to classify these phages into a separate family.

Solving protein structures by combining structure prediction, molecular replacement and direct-methods-aided model completion

Highly accurate protein structure prediction can generate accurate models of protein and protein-protein complexes in X-ray crystallography. However, the question of how to make more effective use of predicted models for completing structure analysis, and which strategies should be employed for the more challenging cases such as multi-helical structures, multimeric structures and extremely large structures, both in the model preparation and in the completion steps, remains open for discussion. In this paper, a new strategy is proposed based on the framework of direct methods and dual-space iteration, which can greatly simplify the pre-processing steps of predicted models both in normal and in challenging cases. Following this strategy, full-length models or the conservative structural domains could be used directly as the starting model, and the phase error and the model bias between the starting model and the real structure would be modified in the direct-methods-based dual-space iteration. Many challenging cases (from CASP14) have been tested for the general applicability of this constructive strategy, and almost complete models have been generated with reasonable statistics. The hybrid strategy therefore provides a meaningful scheme for X-ray structure determination using a predicted model as the starting point.

Structure of the recombinant RNA polymerase from African Swine Fever Virus

African Swine Fever Virus is a Nucleo-Cytoplasmic Large DNA Virus that causes an incurable haemorrhagic fever in pigs with a high impact on global food security. ASFV replicates in the cytoplasm of the infected cell and encodes its own transcription machinery that is independent of cellular factors, however, not much is known about how this system works at a molecular level. Here, we present methods to produce recombinant ASFV RNA polymerase, functional assays to screen for inhibitors, and high-resolution cryo-electron microscopy structures of the ASFV RNAP in different conformational states. The ASFV RNAP bears a striking resemblance to RNAPII with bona fide homologues of nine of its twelve subunits. Key differences include the fusion of the ASFV assembly platform subunits RPB3 and RPB11, and an unusual C-terminal domain of the stalk subunit vRPB7 that is related to the eukaryotic mRNA cap 2´-O-methyltransferase 1. Despite the high degree of structural conservation with cellular RNA polymerases, the ASFV RNAP is resistant to the inhibitors rifampicin and alpha-amanitin. The cryo-EM structures and fully recombinant RNAP system together provide an important tool for the design, development, and screening of antiviral drugs in a low biosafety containment environment.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

Tail-tape-fused virion and non-virion RNA polymerases of a thermophilic virus with an extremely long tail

Thermus thermophilus bacteriophage P23-45 encodes a giant 5,002-residue tail tape measure protein (TMP) that defines the length of its extraordinarily long tail. Here, we show that the N-terminal portion of P23-45 TMP is an unusual RNA polymerase (RNAP) homologous to cellular RNAPs. The TMP-fused virion RNAP transcribes pre-early phage genes, including a gene that encodes another, non-virion RNAP, that transcribes early and some middle phage genes. We report the crystal structures of both P23-45 RNAPs. The non-virion RNAP has a crab-claw-like architecture. By contrast, the virion RNAP adopts a unique flat structure without a clamp. Structure and sequence comparisons of the P23-45 RNAPs with other RNAPs suggest that, despite the extensive functional differences, the two P23-45 RNAPs originate from an ancient gene duplication in an ancestral phage. Our findings demonstrate striking adaptability of RNAPs that can be attained within a single virus species.

Structural atlas of a human gut crassvirus

CrAssphage and related viruses of the order Crassvirales (hereafter referred to as crassviruses) were originally discovered by cross-assembly of metagenomic sequences. They are the most abundant viruses in the human gut, are found in the majority of individual gut viromes, and account for up to 95% of the viral sequences in some individuals1-4. Crassviruses are likely to have major roles in shaping the composition and functionality of the human microbiome, but the structures and roles of most of the virally encoded proteins are unknown, with only generic predictions resulting from bioinformatic analyses4,5. Here we present a cryo-electron microscopy reconstruction of Bacteroides intestinalis virus ΦcrAss0016, providing the structural basis for the functional assignment of most of its virion proteins. The muzzle protein forms an assembly about 1 MDa in size at the end of the tail and exhibits a previously unknown fold that we designate the 'crass fold', that is likely to serve as a gatekeeper that controls the ejection of cargos. In addition to packing the approximately 103 kb of virus DNA, the ΦcrAss001 virion has extensive storage space for virally encoded cargo proteins in the capsid and, unusually, within the tail. One of the cargo proteins is present in both the capsid and the tail, suggesting a general mechanism for protein ejection, which involves partial unfolding of proteins during their extrusion through the tail. These findings provide a structural basis for understanding the mechanisms of assembly and infection of these highly abundant crassviruses.

Bacteriophages of the Order Crassvirales: What Do We Currently Know about This Keystone Component of the Human Gut Virome?

The order Crassvirales comprises dsDNA bacteriophages infecting bacteria in the phylum Bacteroidetes that are found in a variety of environments but are especially prevalent in the mammalian gut. This review summarises available information on the genomics, diversity, taxonomy, and ecology of this largely uncultured viral taxon. With experimental data available from a handful of cultured representatives, the review highlights key properties of virion morphology, infection, gene expression and replication processes, and phage-host dynamics.

Function of the bacteriophage P2 baseplate central spike Apex domain in the infection process

The contractile tail of bacteriophage P2 functions to drive the tail tube across the outer membrane of its host bacterium, a prerequisite event for subsequent translocation of phage genomic DNA into the host cell. The tube is equipped with a spike-shaped protein (product of P2 gene V , gpV or Spike) that contains a membrane-attacking Apex domain carrying a centrally positioned Fe ion. The ion is enclosed in a histidine cage that is formed by three symmetry-related copies of a conserved HxH (histidine, any residue, histidine) sequence motif. Here, we used solution biophysics and X-ray crystallography to characterize the structure and properties of Spike mutants in which the Apex domain was either deleted or its histidine cage was either destroyed or replaced with a hydrophobic core. We found that the Apex domain is not required for the folding of full-length gpV or its middle intertwined β-helical domain. Furthermore, despite its high conservation, the Apex domain is dispensable for infection in laboratory conditions. Collectively, our results show that the diameter of the Spike but not the nature of its Apex domain determines the efficiency of infection, which further strengthens the earlier hypothesis of a drill bit-like function of the Spike in host envelope disruption.

Pectobacterium versatile Bacteriophage Possum: A Complex Polysaccharide-Deacetylating Tail Fiber as a Tool for Host Recognition in Pectobacterial Schitoviridae

Novel, closely related phages Possum and Horatius infect Pectobacterium versatile, a phytopathogen causing soft rot in potatoes and other essential plants. Their properties and genomic composition define them as N4-like bacteriophages of the genus Cbunavirus, a part of a recently formed family Schitoviridae. It is proposed that the adsorption apparatus of these phages consists of tail fibers connected to the virion through an adapter protein. Tail fibers possess an enzymatic domain. Phage Possum uses it to deacetylate O-polysaccharide on the surface of the host strain to provide viral attachment. Such an infection mechanism is supposed to be common for all Cbunavirus phages and this feature should be considered when designing cocktails for phage control of soft rot.

Structural insights into the dual activities of the two-barrel RNA polymerase QDE-1

Neurospora crassa protein QDE-1, a member of the two-barrel polymerase superfamily, possesses both DNA- and RNA-dependent RNA polymerase (DdRP and RdRP) activities. The dual activities are essential for the production of double-stranded RNAs (dsRNAs), the precursors of small interfering RNAs (siRNAs) in N. crassa. Here, we report five complex structures of N-terminal truncated QDE-1 (QDE-1ΔN), representing four different reaction states: DNA/RNA-templated elongation, the de novo initiation of RNA synthesis, the first step of nucleotide condensation during de novo initiation and initial NTP loading. The template strand is aligned by a bridge-helix and double-psi beta-barrels 2 (DPBB2), the RNA product is held by DPBB1 and the slab domain. The DNA template unpairs with the RNA product at position –7, but the RNA template remains paired. The NTP analog coordinates with cations and is precisely positioned at the addition site by a rigid trigger loop and a proline-containing loop in the active center. The unique C-terminal tail from the QDE-1 dimer partner inserts into the substrate-binding cleft and plays regulatory roles in RNA synthesis. Collectively, this work elucidates the conserved mechanisms for DNA/RNA-dependent dual activities by QDE-1 and other two-barrel polymerase superfamily members.

Genomics and Geographic Diversity of Bacteriophages Associated With Endosymbionts in the Guts of Workers and Alates of Coptotermes Species (Blattodea: Rhinotermitidae)

Subterranean termites depend nutritionally on their gut microbiota, which includes protozoa as well as taxonomically and functionally diverse bacteria. Our previous metavirome study revealed a high diversity and novel families of bacteriophages in the guts of Coptotermes formosanus workers from New Orleans, Louisiana, United States. Two assembled bacteriophage genomes (Phages TG-crAlp-04 and 06, family Podoviridae) existed in all colonies and showed similarity to a prophage (ProJPt-Bp1) previously sequenced from a bacterial endosymbiont (Candidatus Azobacteroides pseudotrichonymphae, CAP) of protozoa in the gut of a termite species of the genus Prorhinotermes from Taiwan. In this study the genomes of Phage TG-crAlp-04 and 06 were subjected to detailed functional annotation. Both phage genomes contained conserved genes for DNA packaging, head and tail morphogenesis, and phage replication. Approximately 30% of the amino acid sequences derived from genes in both genomes matched to those of ProJPt-Bp1 phage or other phages from the crAss-like phage group. No integrase was identified; the lack of a lysogeny module is a characteristic of crAss-like phages. Primers were designed to sequence conserved genes of the two phages and their putative host bacterium (CAP) to detect their presence in different termite species from native and introduced distribution ranges. Related strains of the host bacterium were found across different termite genera and geographic regions. Different termite species had separate CAP strains, but intraspecific geographical variation was low. These results together with the fact that CAP is an important intracellular symbiont of obligate cellulose-digesting protozoa, suggest that CAP is a core gut bacterium and co-evolved across several subterranean termite species. Variants of both crAss-like phages were detected in different Coptotermes species from the native and introduced range, but they did not differentiate by species or geographic region. Since similar phages were detected in different termite species, we propose the existence of a core virome associated with core bacterial endosymbionts of protozoa in the guts of subterranean termites. This work provides a strong basis for further study of the quadripartite relationship of termites, protozoa, bacteria, and bacteriophages.

Recruitment of Mobile Genetic Elements for Diverse Cellular Functions in Prokaryotes

Prokaryotic genomes are replete with mobile genetic elements (MGE) that span a continuum of replication autonomy. On numerous occasions during microbial evolution, diverse MGE lose their autonomy altogether but, rather than being quickly purged from the host genome, assume a new function that benefits the host, rendering the immobilized MGE subject to purifying selection, and resulting in its vertical inheritance. This mini-review highlights the diversity of the repurposed (exapted) MGE as well as the plethora of cellular functions that they perform. The principal contribution of the exaptation of MGE and their components is to the prokaryotic functional systems involved in biological conflicts, and in particular, defense against viruses and other MGE. This evolutionary entanglement between MGE and defense systems appears to stem both from mechanistic similarities and from similar evolutionary predicaments whereby both MGEs and defense systems tend to incur fitness costs to the hosts and thereby evolve mechanisms for survival including horizontal mobility, causing host addiction, and exaptation for functions beneficial to the host. The examples discussed demonstrate that the identity of an MGE, overall mobility and relationship with the host cell (mutualistic, symbiotic, commensal, or parasitic) are all factors that affect exaptation.

Discovery, diversity, and functional associations of crAss-like phages in human gut metagenomes from four Dutch cohorts

The crAss-like phages are a diverse group of related viruses that includes some of the most abundant viruses of the human gut. To explore their diversity and functional role in human population and clinical cohorts, we analyze gut metagenomic data collected from 1,950 individuals from the Netherlands. We identify 1,556 crAss-like phage genomes, including 125 species-level and 32 genus-level clusters absent from the reference databases used. Analysis of their genomic features shows that closely related crAss-like phages can possess strikingly divergent regions responsible for transcription, presumably acquired through recombination. Prediction of crAss-like phage hosts points primarily to bacteria of the phylum Bacteroidetes, consistent with previous reports. Finally, we explore the temporal stability of crAss-like phages over a 4-year period and identify associations between the abundance of crAss-like phages and several human phenotypes, including depletion of crAss-like phages in inflammatory bowel disease patients.

Implications of AlphaFold2 for crystallographic phasing by molecular replacement

The AlphaFold2 results in the 14th edition of Critical Assessment of Structure Prediction (CASP14) showed that accurate (low root-mean-square deviation) in silico models of protein structure domains are on the horizon, whether or not the protein is related to known structures through high-coverage sequence similarity. As highly accurate models become available, generated by harnessing the power of correlated mutations and deep learning, one of the aspects of structural biology to be impacted will be methods of phasing in crystallography. Here, the data from CASP14 are used to explore the prospects for changes in phasing methods, and in particular to explore the prospects for molecular-replacement phasing using in silico models.

Structure and RNA template requirements of Arabidopsis RNA-DEPENDENT RNA POLYMERASE 2

RNA-dependent RNA polymerases play essential roles in RNA-mediated gene silencing in eukaryotes. In Arabidopsis, RNA-DEPENDENT RNA POLYMERASE 2 (RDR2) physically interacts with DNA-dependent NUCLEAR RNA POLYMERASE IV (Pol IV) and their activities are tightly coupled, with Pol IV transcriptional arrest, induced by the nontemplate DNA strand, somehow enabling RDR2 to engage Pol IV transcripts and generate double-stranded RNAs. The double-stranded RNAs are then released from the Pol IV-RDR2 complex and diced into short-interfering RNAs that guide RNA-directed DNA methylation and silencing. Here we report the structure of full-length RDR2, at an overall resolution of 3.1 Å, determined by cryoelectron microscopy. The N-terminal region contains an RNA-recognition motif adjacent to a positively charged channel that leads to a catalytic center with striking structural homology to the catalytic centers of multisubunit DNA-dependent RNA polymerases. We show that RDR2 initiates 1 to 2 nt internal to the 3' ends of its templates and can transcribe the RNA of an RNA/DNA hybrid, provided that 9 or more nucleotides are unpaired at the RNA's 3' end. Using a nucleic acid configuration that mimics the arrangement of RNA and DNA strands upon Pol IV transcriptional arrest, we show that displacement of the RNA 3' end occurs as the DNA template and nontemplate strands reanneal, enabling RDR2 transcription. These results suggest a model in which Pol IV arrest and backtracking displaces the RNA 3' end as the DNA strands reanneal, allowing RDR2 to engage the RNA and synthesize the complementary strand.

Protein inter-residue contact and distance prediction by coupling complementary coevolution features with deep residual networks in CASP14

This article reports and analyzes the results of protein contact and distance prediction by our methods in the 14th Critical Assessment of techniques for protein Structure Prediction (CASP14). A new deep learning-based contact/distance predictor was employed based on the ensemble of two complementary coevolution features coupling with deep residual networks. We also improved our multiple sequence alignment (MSA) generation protocol with wholesale meta-genome sequence databases. On 22 CASP14 free modeling (FM) targets, the proposed model achieved a top-L/5 long-range precision of 63.8% and a mean distance bin error of 1.494. Based on the predicted distance potentials, 11 out of 22 FM targets and all of the 14 FM/template-based modeling (TBM) targets have correctly predicted folds (TM-score >0.5), suggesting that our approach can provide reliable distance potentials for ab initio protein folding.

Physics-based protein structure refinement in the era of artificial intelligence

Protein structure refinement is the last step in protein structure prediction pipelines. Physics-based refinement via molecular dynamics (MD) simulations has made significant progress during recent years. During CASP14, we tested a new refinement protocol based on an improved sampling strategy via MD simulations. MD simulations were carried out at an elevated temperature (360 K). An optimized use of biasing restraints and the use of multiple starting models led to enhanced sampling. The new protocol generally improved the model quality. In comparison with our previous protocols, the CASP14 protocol showed clear improvements. Our approach was successful with most initial models, many based on deep learning methods. However, we found that our approach was not able to refine machine-learning models from the AlphaFold2 group, often decreasing already high initial qualities. To better understand the role of refinement given new types of models based on machine-learning, a detailed analysis via MD simulations and Markov state modeling is presented here. We continue to find that MD-based refinement has the potential to improve AI predictions. We also identified several practical issues that make it difficult to realize that potential. Increasingly important is the consideration of inter-domain and oligomeric contacts in simulations; the presence of large kinetic barriers in refinement pathways also continues to present challenges. Finally, we provide a perspective on how physics-based refinement could continue to play a role in the future for improving initial predictions based on machine learning-based methods.

Mechanisms and clinical importance of bacteriophage resistance

We are in the midst of a golden age of uncovering defense systems against bacteriophages. Apart from the fundamental interest in these defense systems, and revolutionary applications that have been derived from them (e.g. CRISPR-Cas9 and restriction endonucleases), it is unknown how defense systems contribute to resistance formation against bacteriophages in clinical settings. Bacteriophages are now being reconsidered as therapeutic agents against bacterial infections due the emergence of multidrug resistance. However, bacteriophage resistance through defense systems and other means could hinder the development of successful phage-based therapies. Here, we review the current state of the field of bacteriophage defense, highlight the relevance of bacteriophage defense for potential clinical use of bacteriophages as therapeutic agents and suggest new directions of research.

Long-term persistence of crAss-like phage crAss001 is associated with phase variation in Bacteroides intestinalis

BACKGROUND

The crAss-like phages are ubiquitous and highly abundant members of the human gut virome that infect commensal bacteria of the order Bacteroidales. Although incapable of lysogeny, these viruses demonstrate long-term persistence in the human gut microbiome, dominating the virome in some individuals.

RESULTS

Here we show that rapid phase variation of alternate capsular polysaccharides in Bacteroides intestinalis cultures plays an important role in a dynamic equilibrium between phage sensitivity and resistance, allowing phage and bacteria to multiply in parallel. The data also suggests the role of a concomitant phage persistence mechanism associated with delayed lysis of infected cells, similar to carrier state infection. From an ecological and evolutionary standpoint, this type of phage-host interaction is consistent with the Piggyback-the-Winner model, which suggests a preference towards lysogenic or other "benign" forms of phage infection when the host is stably present at high abundance.

CONCLUSION

Long-term persistence of bacteriophage and host could result from mutually beneficial mechanisms driving bacterial strain-level diversity and phage survival in complex environments.

Target classification in the 14th round of the critical assessment of protein structure prediction (CASP14)

An evolutionary-based definition and classification of target evaluation units (EUs) is presented for the 14th round of the critical assessment of structure prediction (CASP14). CASP14 targets included 84 experimental models submitted by various structural groups (designated T1024-T1101). Targets were split into EUs based on the domain organization of available templates and performance of server groups. Several targets required splitting (19 out of 25 multidomain targets) due in part to observed conformation changes. All in all, 96 CASP14 EUs were defined and assigned to tertiary structure assessment categories (Topology-based FM or High Accuracy-based TBM-easy and TBM-hard) considering their evolutionary relationship to existing ECOD fold space: 24 family level, 50 distant homologs (H-group), 12 analogs (X-group), and 10 new folds. Principal component analysis and heatmap visualization of sequence and structure similarity to known templates as well as performance of servers highlighted trends in CASP14 target difficulty. The assigned evolutionary levels (i.e., H-groups) and assessment classes (i.e., FM) displayed overlapping clusters of EUs. Many viral targets diverged considerably from their template homologs and thus were more difficult for prediction than other homology-related targets. On the other hand, some targets did not have sequence-identifiable templates, but were predicted better than expected due to relatively simple arrangements of secondary structural elements. An apparent improvement in overall server performance in CASP14 further complicated traditional classification, which ultimately assigned EUs into high-accuracy modeling (27 TBM-easy and 31 TBM-hard), topology (23 FM), or both (15 FM/TBM).

Highly accurate protein structure prediction with AlphaFold

Proteins are essential to life, and understanding their structure can facilitate a mechanistic understanding of their function. Through an enormous experimental effort1-4, the structures of around 100,000 unique proteins have been determined5, but this represents a small fraction of the billions of known protein sequences6,7. Structural coverage is bottlenecked by the months to years of painstaking effort required to determine a single protein structure. Accurate computational approaches are needed to address this gap and to enable large-scale structural bioinformatics. Predicting the three-dimensional structure that a protein will adopt based solely on its amino acid sequence-the structure prediction component of the 'protein folding problem'8-has been an important open research problem for more than 50 years9. Despite recent progress10-14, existing methods fall far short of atomic accuracy, especially when no homologous structure is available. Here we provide the first computational method that can regularly predict protein structures with atomic accuracy even in cases in which no similar structure is known. We validated an entirely redesigned version of our neural network-based model, AlphaFold, in the challenging 14th Critical Assessment of protein Structure Prediction (CASP14)15, demonstrating accuracy competitive with experimental structures in a majority of cases and greatly outperforming other methods. Underpinning the latest version of AlphaFold is a novel machine learning approach that incorporates physical and biological knowledge about protein structure, leveraging multi-sequence alignments, into the design of the deep learning algorithm.

Atomistic basis of force generation, translocation, and coordination in a viral genome packaging motor

Double-stranded DNA viruses package their genomes into pre-assembled capsids using virally-encoded ASCE ATPase ring motors. We present the first atomic-resolution crystal structure of a multimeric ring form of a viral dsDNA packaging motor, the ATPase of the asccφ28 phage, and characterize its atomic-level dynamics via long timescale molecular dynamics simulations. Based on these results, and previous single-molecule data and cryo-EM reconstruction of the homologous φ29 motor, we propose an overall packaging model that is driven by helical-to-planar transitions of the ring motor. These transitions are coordinated by inter-subunit interactions that regulate catalytic and force-generating events. Stepwise ATP binding to individual subunits increase their affinity for the helical DNA phosphate backbone, resulting in distortion away from the planar ring towards a helical configuration, inducing mechanical strain. Subsequent sequential hydrolysis events alleviate the accumulated mechanical strain, allowing a stepwise return of the motor to the planar conformation, translocating DNA in the process. This type of helical-to-planar mechanism could serve as a general framework for ring ATPases.

Fishing for phages in metagenomes: what do we catch, what do we miss?

Metagenomics and metatranscriptomics have become the principal approaches for discovery of novel bacteriophages and preliminary characterization of their ecology and biology. Metagenomic sequencing dramatically expanded the known diversity of tailed and non-tailed phages with double-stranded DNA genomes and those with single-stranded DNA genomes, whereas metatranscriptomics led to the discovery of thousands of new single-stranded RNA phages. Apart from expanding phage diversity, metagenomics studies discover major novel groups of phages with unique features of genome organization, expression strategy and virus-host interaction, such as the putative order 'crAssvirales', which includes the most abundant human-associated viruses. The continued success of metagenomics hinges on the combination of the most powerful computational methods for phage genome assembly and analysis including harnessing CRISPR spacers for the discovery of novel phages and host assignment. Together, these approaches could make a comprehensive characterization of the earth phageome a realistic goal.

Analysis of metagenome-assembled viral genomes from the human gut reveals diverse putative CrAss-like phages with unique genomic features

CrAssphage is the most abundant human-associated virus and the founding member of a large group of bacteriophages, discovered in animal-associated and environmental metagenomes, that infect bacteria of the phylum Bacteroidetes. We analyze 4907 Circular Metagenome Assembled Genomes (cMAGs) of putative viruses from human gut microbiomes and identify nearly 600 genomes of crAss-like phages that account for nearly 87% of the DNA reads mapped to these cMAGs. Phylogenetic analysis of conserved genes demonstrates the monophyly of crAss-like phages, a putative virus order, and of 5 branches, potential families within that order, two of which have not been identified previously. The phage genomes in one of these families are almost twofold larger than the crAssphage genome (145-192 kilobases), with high density of self-splicing introns and inteins. Many crAss-like phages encode suppressor tRNAs that enable read-through of UGA or UAG stop-codons, mostly, in late phage genes. A distinct feature of the crAss-like phages is the recurrent switch of the phage DNA polymerase type between A and B families. Thus, comparative genomic analysis of the expanded assemblage of crAss-like phages reveals aspects of genome architecture and expression as well as phage biology that were not apparent from the previous work on phage genomics.