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

CrAss-like phages are suitable indicators of antibiotic resistance genes found in abundance in fecally polluted samples

Antibiotic resistance genes (ARGs) have been extensively observed in bacterial DNA, and more recently, in phage particles from various water sources and food items. The pivotal role played by ARG transmission in the proliferation of antibiotic resistance and emergence of new resistant strains calls for a thorough understanding of the underlying mechanisms. The aim of this study was to assess the suitability of the prototypical p-crAssphage, a proposed indicator of human fecal contamination, and the recently isolated crAssBcn phages, both belonging to the Crassvirales group, as potential indicators of ARGs. These crAss-like phages were evaluated alongside specific ARGs (blaTEM, blaCTX-M-1, blaCTX-M-9, blaVIM, blaOXA-48, qnrA, qnrS, tetW and sul1) within the total DNA and phage DNA fractions in water and food samples containing different levels of fecal pollution. In samples with high fecal load (>103 CFU/g or ml of E. coli or somatic coliphages), such as wastewater and sludge, positive correlations were found between both types of crAss-like phages and ARGs in both DNA fractions. The strongest correlation was observed between sul1 and crAssBcn phages (rho = 0.90) in sludge samples, followed by blaCTX-M-9 and p-crAssphage (rho = 0.86) in sewage samples, both in the phage DNA fraction. The use of crAssphage and crAssBcn as indicators of ARGs, considered to be emerging environmental contaminants of anthropogenic origin, is supported by their close association with the human gut. Monitoring ARGs can help to mitigate their dissemination and prevent the emergence of new resistant bacterial strains, thus safeguarding public health.

DNA polymerase swapping in Caudoviricetes bacteriophages

BACKGROUND

Viruses with double-stranded (ds) DNA genomes in the realm Duplodnaviria share a conserved structural gene module but show a broad range of variation in their repertoires of DNA replication proteins. Some of the duplodnaviruses encode (nearly) complete replication systems whereas others lack (almost) all genes required for replication, relying on the host replication machinery. DNA polymerases (DNAPs) comprise the centerpiece of the DNA replication apparatus. The replicative DNAPs are classified into 4 unrelated or distantly related families (A-D), with the protein structures and sequences within each family being, generally, highly conserved. More than half of the duplodnaviruses encode a DNAP of family A, B or C. We showed previously that multiple pairs of closely related viruses in the order Crassvirales encode DNAPs of different families.

METHODS

Groups of phages in which DNAP swapping likely occurred were identified as subtrees of a defined depth in a comprehensive evolutionary tree of tailed bacteriophages that included phages with DNAPs of different families. The DNAP swaps were validated by constrained tree analysis that was performed on phylogenetic tree of large terminase subunits, and the phage genomes encoding swapped DNAPs were aligned using Mauve. The structures of the discovered unusual DNAPs were predicted using AlphaFold2.

RESULTS

We identified four additional groups of tailed phages in the class Caudoviricetes in which the DNAPs apparently were swapped on multiple occasions, with replacements occurring both between families A and B, or A and C, or between distinct subfamilies within the same family. The DNAP swapping always occurs "in situ", without changes in the organization of the surrounding genes. In several cases, the DNAP gene is the only region of substantial divergence between closely related phage genomes, whereas in others, the swap apparently involved neighboring genes encoding other proteins involved in phage genome replication. In addition, we identified two previously undetected, highly divergent groups of family A DNAPs that are encoded in some phage genomes along with the main DNAP implicated in genome replication.

CONCLUSIONS

Replacement of the DNAP gene by one encoding a DNAP of a different family occurred on many independent occasions during the evolution of different families of tailed phages, in some cases, resulting in very closely related phages encoding unrelated DNAPs. DNAP swapping was likely driven by selection for avoidance of host antiphage mechanisms targeting the phage DNAP that remain to be identified, and/or by selection against replicon incompatibility.

Identification of mobile genetic elements with geNomad

Identifying and characterizing mobile genetic elements in sequencing data is essential for understanding their diversity, ecology, biotechnological applications and impact on public health. Here we introduce geNomad, a classification and annotation framework that combines information from gene content and a deep neural network to identify sequences of plasmids and viruses. geNomad uses a dataset of more than 200,000 marker protein profiles to provide functional gene annotation and taxonomic assignment of viral genomes. Using a conditional random field model, geNomad also detects proviruses integrated into host genomes with high precision. In benchmarks, geNomad achieved high classification performance for diverse plasmids and viruses (Matthews correlation coefficient of 77.8% and 95.3%, respectively), substantially outperforming other tools. Leveraging geNomad's speed and scalability, we processed over 2.7 trillion base pairs of sequencing data, leading to the discovery of millions of viruses and plasmids that are available through the IMG/VR and IMG/PR databases. geNomad is available at https://portal.nersc.gov/genomad .

Gut virome-wide association analysis identifies cross-population viral signatures for inflammatory bowel disease

BACKGROUND

The gut virome has been implicated in inflammatory bowel disease (IBD), yet a full understanding of the gut virome in IBD patients, especially across diverse geographic populations, is lacking.

RESULTS

In this study, we conducted a comprehensive gut virome-wide association study in a Chinese cohort of 71 IBD patients (15 with Crohn's disease and 56 with ulcerative colitis) and 77 healthy controls via viral-like particle (VLP) and bulk virome sequencing of their feces. By utilizing an integrated gut virus catalog tailored to the IBD virome, we revealed fundamental alterations in the gut virome in IBD patients. These characterized 139 differentially abundant viral signatures, including elevated phages predicted to infect Escherichia, Klebsiella, Enterococcus_B, Streptococcus, and Veillonella species, as well as IBD-depleted phages targeting Prevotella, Ruminococcus_E, Bifidobacterium, and Blautia species. Remarkably, these viral signatures demonstrated high consistency across diverse populations such as those in Europe and the USA, emphasizing their significance and broad relevance in the disease context. Furthermore, fecal virome transplantation experiments verified that the colonization of these IBD-characterized viruses can modulate experimental colitis in mouse models.

CONCLUSIONS

Building upon these insights into the IBD gut virome, we identified potential biomarkers for prognosis and therapy in IBD patients, laying the foundation for further exploration of viromes in related conditions. Video Abstract.

Mriyaviruses: small relatives of giant viruses

The phylum Nucleocytoviricota consists of large and giant viruses that range in genome size from about 100 kilobases (kb) to more than 2.5 megabases. Here, using metagenome mining followed by extensive phylogenomic analysis and protein structure comparison, we delineate a distinct group of viruses with double-stranded (ds) DNA genomes in the range of 35-45 kb that appear to be related to the Nucleocytoviricota. In phylogenetic trees of the conserved double jelly-roll major capsid proteins (MCPs) and DNA packaging ATPases, these viruses do not show affinity to any particular branch of the Nucleocytoviricota and accordingly would comprise a class which we propose to name "Mriyaviricetes" (after Ukrainian "mriya," dream). Structural comparison of the MCP suggests that, among the extant virus lineages, mriyaviruses are the closest one to the ancestor of the Nucleocytoviricota. In the phylogenetic trees, mriyaviruses split into two well-separated branches, the family Yaraviridae and proposed new family "Gamadviridae." The previously characterized members of these families, yaravirus and Pleurochrysis sp. endemic viruses, infect amoeba and haptophytes, respectively. The genomes of the rest of the mriyaviruses were assembled from metagenomes from diverse environments, suggesting that mriyaviruses infect various unicellular eukaryotes. Mriyaviruses lack DNA polymerase, which is encoded by all other members of the Nucleocytoviricota, and RNA polymerase subunits encoded by all cytoplasmic viruses among the Nucleocytoviricota, suggesting that they replicate in the host cell nuclei. All mriyaviruses encode a HUH superfamily endonuclease that is likely to be essential for the initiation of virus DNA replication via the rolling circle mechanism.

IMPORTANCE

The origin of giant viruses of eukaryotes that belong to the phylum Nucleocytoviricota is not thoroughly understood and remains a matter of major interest and debate. Here, we combine metagenome database searches with extensive protein sequence and structure analysis to describe a distinct group of viruses with comparatively small genomes of 35-45 kilobases that appear to comprise a distinct class within the phylum Nucleocytoviricota that we provisionally named "Mriyaviricetes." Mriyaviruses appear to be the closest identified relatives of the ancestors of the Nucleocytoviricota. Analysis of proteins encoded in mriyavirus genomes suggests that they replicate their genome via the rolling circle mechanism that is unusual among viruses with double-stranded DNA genomes and so far not described for members of Nucleocytoviricota.

VirRep: a hybrid language representation learning framework for identifying viruses from human gut metagenomes

Identifying viruses from metagenomes is a common step to explore the virus composition in the human gut. Here, we introduce VirRep, a hybrid language representation learning framework, for identifying viruses from human gut metagenomes. VirRep combines a context-aware encoder and an evolution-aware encoder to improve sequence representation by incorporating k-mer patterns and sequence homologies. Benchmarking on both simulated and real datasets with varying viral proportions demonstrates that VirRep outperforms state-of-the-art methods. When applied to fecal metagenomes from a colorectal cancer cohort, VirRep identifies 39 high-quality viral species associated with the disease, many of which cannot be detected by existing methods.

CrAss-Like Phages: From Discovery in Human Fecal Metagenome to Application as a Microbial Source Tracking Marker

CrAss-like phages are a diverse group of bacteriophages genetically similar to the prototypical crAssphage (p-crAssphage), which was discovered in the human gut microbiome through a metagenomics approach. It was identified as a ubiquitous and highly abundant bacteriophage group in the gut microbiome. Initial co-occurrence analysis postulated Bacteroides spp. as the prospective bacterial host. Subsequent studies have confirmed multiple host species under Phylum Bacteroidetes and some Firmicutes. Detection of crAss-like phages in sewage-contaminated environmental water and robust correlation with enteric viruses and bacteria has culminated in their adoption as a microbial source tracking (MST) marker. Polymerase chain reaction (PCR) and real-time PCR assays have been developed utilizing the conserved genes in the p-crAssphage genome to detect human fecal contamination of different water sources, with high specificity. Numerous investigations have examined the implications of crAss-like phages in diverse disease conditions, including ulcerative colitis, obesity and metabolic syndrome, autism spectrum disorders, rheumatoid arthritis, atopic eczema, and other autoimmune disorders. These studies have unveiled associations between certain diseases and diminished abundance and diversity of crAss-like phages. This review offers insights into the diverse aspects of research on crAss-like phages, including their discovery, genomic characteristics, structure, taxonomy, isolation, molecular detection, application as an MST marker, and role as a gut microbiome modulator with consequential health implications.

Assessing phage-host population dynamics by reintroducing virulent viruses to synthetic microbiomes

Microbiomes feature complex interactions between diverse bacteria and bacteriophages. Synthetic microbiomes offer a powerful way to study these interactions; however, a major challenge is obtaining a representative bacteriophage population during the bacterial isolation process. We demonstrate that colony isolation reliably excludes virulent viruses from sample sources with low virion-to-bacteria ratios such as feces, creating "virulent virus-free" controls. When the virulent dsDNA virome is reintroduced to a 73-strain synthetic gut microbiome in a bioreactor model of the human colon, virulent viruses target susceptible strains without significantly altering community structure or metabolism. In addition, we detected signals of prophage induction that associate with virulent predation. Overall, our findings indicate that dilution-based isolation methods generate synthetic gut microbiomes that are heavily depleted, if not devoid, of virulent viruses and that such viruses, if reintroduced, have a targeted effect on community assembly, metabolism, and prophage replication.

Decoding huge phage diversity: a taxonomic classification of Lak megaphages

High-throughput sequencing for uncultivated viruses has accelerated the understanding of global viral diversity and uncovered viral genomes substantially larger than any that have so far been cultured. Notably, the Lak phages are an enigmatic group of viruses that present some of the largest known phage genomes identified in human and animal microbiomes, and are dissimilar to any cultivated viruses. Despite the wealth of viral diversity that exists within sequencing datasets, uncultivated viruses have rarely been used for taxonomic classification. We investigated the evolutionary relationships of 23 Lak phages and propose a taxonomy for their classification. Predicted protein analysis revealed the Lak phages formed a deeply branching monophyletic clade within the class Caudoviricetes which contained no other phage genomes. One of the interesting features of this clade is that all current members are characterised by an alternative genetic code. We propose the Lak phages belong to a new order, the 'Grandevirales'. Protein and nucleotide-based analyses support the creation of two families, three sub-families, and four genera within the order 'Grandevirales'. We anticipate that the proposed taxonomy of Lak megaphages will simplify the future classification of related viral genomes as they are uncovered. Continued efforts to classify divergent viruses are crucial to aid common analyses of viral genomes and metagenomes.

Jumping DNA polymerases in bacteriophages

Viruses with double-stranded (ds) DNA genomes in the realm Duplodnaviria share a conserved structural gene module but show a broad range of variation in their repertoires of DNA replication proteins. Some of the duplodnaviruses encode (nearly) complete replication systems whereas others lack (almost) all genes required for replication, relying on the host replication machinery. DNA polymerases (DNAPs) comprise the centerpiece of the DNA replication apparatus. The replicative DNAPs are classified into 4 unrelated or distantly related families (A-D), with the protein structures and sequences within each family being, generally, highly conserved. More than half of the duplodnaviruses encode a DNAP of family A, B or C. We showed previously that multiple pairs of closely related viruses in the order Crassvirales encode DNAPs of different families. Here we identify four additional groups of tailed phages in the class Caudoviricetes in which the DNAPs apparently were swapped on multiple occasions, with replacements occurring both between families A and B, or A and C, or between distinct subfamilies within the same family. The DNAP swapping always occurs "in situ", without changes in the organization of the surrounding genes. In several cases, the DNAP gene is the only region of substantial divergence between closely related phage genomes, whereas in others, the swap apparently involved neighboring genes encoding other proteins involved in phage replication. We hypothesize that DNAP swapping is driven by selection for avoidance of host antiphage mechanisms targeting the phage DNAP that remain to be identified, and/or by selection against replicon incompatibility. In addition, we identified two previously undetected, highly divergent groups of family A DNAPs that are encoded in some phage genomes along with the main DNAP implicated in genome replication.

Why tRNA acquisition could be relevant to bacteriophages?

In this opinion, we discuss the role of tRNAs in phage biology and their importance in DNA replication and phage-host interactions. Phages are a diverse group of obligate bacterial viruses that possess genomes with a wide range of sizes. Among them, we find phages with few genes that depend entirely on their host's translational machinery for replication. However, some phages carry genes for all replication steps and even contain genes for their own translational synthesis. In these cases, the integration of tRNA genes in their genomes is not completely understood, generating different theories about their presence and function during the replication cycle. Although different studies have attempted to elucidate their role, additional studies are needed to clarify the presence and significance of tRNA genes in phages. Moreover, we highlight the importance of tRNA genes in phages from both ecological and therapeutic perspectives.

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.

Analysis and culturing of the prototypic crAssphage reveals a phage-plasmid lifestyle

The prototypic crAssphage (Carjivirus communis) is one of the most abundant, prevalent, and persistent gut bacteriophages, yet it remains uncultured and its lifestyle uncharacterized. For the last decade, crAssphage has escaped plaque-dependent culturing efforts, leading us to investigate alternative lifestyles that might explain its widespread success. Through genomic analyses and culturing, we find that crAssphage uses a phage-plasmid lifestyle to persist extrachromosomally. Plasmid-related genes are more highly expressed than those implicated in phage maintenance. Leveraging this finding, we use a plaque-free culturing approach to measure crAssphage replication in culture with Phocaeicola vulgatus, Phocaeicola dorei, and Bacteroides stercoris, revealing a broad host range. We demonstrate that crAssphage persists with its hosts in culture without causing major cell lysis events or integrating into host chromosomes. The ability to switch between phage and plasmid lifestyles within a wide range of hosts contributes to the prolific nature of crAssphage in the human gut microbiome.

Transmission and dynamics of mother-infant gut viruses during pregnancy and early life

Early development of the gut ecosystem is crucial for lifelong health. While infant gut bacterial communities have been studied extensively, the infant gut virome remains under-explored. To study the development of the infant gut virome over time and the factors that shape it, we longitudinally assess the composition of gut viruses and their bacterial hosts in 30 women during and after pregnancy and in their 32 infants during their first year of life. Using shotgun metagenomic sequencing applied to dsDNA extracted from Virus-Like Particles (VLPs) and bacteria, we generate 205 VLP metaviromes and 322 total metagenomes. With this data, we show that while the maternal gut virome composition remains stable during late pregnancy and after birth, the infant gut virome is dynamic in the first year of life. Notably, infant gut viromes contain a higher abundance of active temperate phages compared to maternal gut viromes, which decreases over the first year of life. Moreover, we show that the feeding mode and place of delivery influence the gut virome composition of infants. Lastly, we provide evidence of co-transmission of viral and bacterial strains from mothers to infants, demonstrating that infants acquire some of their virome from their mother's gut.

Mriyaviruses: Small Relatives of Giant Viruses

The phylum Nucleocytoviricota consists of large and giant viruses that range in genome size from about 100 kilobases (kb) to more than 2.5 megabases. Here, using metagenome mining followed by extensive phylogenomic analysis and protein structure comparison, we delineate a distinct group of viruses with double-stranded (ds) DNA genomes in the range of 35-45 kb that appear to be related to the Nucleocytoviricota. In phylogenetic trees of the conserved double jelly-roll major capsid proteins (MCP) and DNA packaging ATPases, these viruses do not show affinity to any particular branch of the Nucleocytoviricota and accordingly would comprise a class which we propose to name "Mriyaviricetes" (after Ukrainian Mriya, dream). Structural comparison of the MCP suggests that, among the extant virus lineages, mriyaviruses are the closest one to the ancestor of the Nucleocytoviricota. In the phylogenetic trees, mriyaviruses split into two well-separated branches, the family Yaraviridae and proposed new family "Gamadviridae". The previously characterized members of these families, Yaravirus and Pleurochrysis sp. endemic viruses, infect amoeba and haptophytes, respectively. The genomes of the rest of the mriyaviruses were assembled from metagenomes from diverse environments, suggesting that mriyaviruses infect various unicellular eukaryotes. Mriyaviruses lack DNA polymerase, which is encoded by all other members of the Nucleocytoviricota, and RNA polymerase subunits encoded by all cytoplasmic viruses among the Nucleocytoviricota, suggesting that they replicate in the host cell nuclei. All mriyaviruses encode a HUH superfamily endonuclease that is likely to be essential for the initiation of virus DNA replication via the rolling circle mechanism.

Assessment of crAssphage as a human fecal source tracking marker in the lower Great Lakes

CrAssphage or crAss-like phage ranks as the most abundant phage in the human gut and is present in human feces-contaminated environments. Due to its high human specificity and sensitivity, crAssphage is a potentially robust source tracking indicator that can distinguish human fecal contamination from agricultural or wildlife sources. Its suitability in the Great Lakes area, one of the world's most important water systems, has not been well tested. In this study, we tested a qPCR-based quantification method using two crAssphage marker genes (ORF18-mod and CPQ_064) at Toronto recreational beaches along with their adjacent river mouths. Our results showed a 71.4 % (CPQ_064) and 100 % (ORF18-mod) human sensitivity for CPQ_064 and ORF18-mod, and a 100 % human specificity for both marker genes. CrAssphage was present in 57.7 % or 71.2 % of environmental water samples, with concentrations ranging from 1.45 to 5.14 log10 gene copies per 100 mL water. Though concentrations of the two marker genes were strongly correlated, ORF18-mod features a higher human sensitivity and higher positive detection rates in environmental samples. Quantifiable crAssphage was mostly present in samples collected in June and July 2021 associated with higher rainfall. In addition, rivers had more frequent crAssphage presence and higher concentrations than their associated beaches, indicating more frequent and greater human fecal contamination in the rivers. However, crAssphage was more correlated with E. coli and Enterococcus at the beaches than in the rivers, suggesting human fecal sources may be more predominant in driving the increases in E. coli and Enterococcus at the beaches when impacted by river plumes.

Dynamic nature of viral and bacterial communities in human faeces

Bacteriophages are a major component of the gut microbiome and are believed to play a role in establishment and stabilization of microbial communities by influencing taxonomic and functional diversity. We show that the activity of lytic and temperate phages can also significantly affect bacterial community structure in a model of extended colonic retention. Intact fresh human feces were incubated anaerobically at 37°C without homogenization and subjected to metagenomic sequencing. We observed subject-specific blooms and collapses of selected bacteriophage and bacterial populations within some individuals. Most notable were striking collapses of Prevotella populations accompanied by increases in specific bacteriophages. In a number of cases, we even observed a shift from one bacterial "enterotype" to another within 48 h. These results confirm that intact feces represents a highly dynamic ecological system and suggests that colonic retention time could have a profound effect on microbiome composition, including a significant impact by bacteriophages.

Eco-evolutionary dynamics of gut phageome in wild gibbons (Hoolock tianxing) with seasonal diet variations

It has been extensively studied that the gut microbiome provides animals flexibility to adapt to food variability. Yet, how gut phageome responds to diet variation of wild animals remains unexplored. Here, we analyze the eco-evolutionary dynamics of gut phageome in six wild gibbons (Hoolock tianxing) by collecting individually-resolved fresh fecal samples and parallel feeding behavior data for 15 consecutive months. Application of complementary viral and microbial metagenomics recovers 39,198 virulent and temperate phage genomes from the feces. Hierarchical cluster analyses show remarkable seasonal diet variations in gibbons. From high-fruit to high-leaf feeding period, the abundances of phage populations are seasonally fluctuated, especially driven by the increased abundance of virulent phages that kill the Lachnospiraceae hosts, and a decreased abundance of temperate phages that piggyback the Bacteroidaceae hosts. Functional profiling reveals an enrichment through horizontal gene transfers of toxin-antitoxin genes on temperate phage genomes in high-leaf season, potentially conferring benefits to their prokaryotic hosts. The phage-host ecological dynamics are driven by the coevolutionary processes which select for tail fiber and DNA primase genes on virulent and temperate phage genomes, respectively. Our results highlight complex phageome-microbiome interactions as a key feature of the gibbon gut microbial ecosystem responding to the seasonal diet.

Infant gut DNA bacteriophage strain persistence during the first 3 years of life

Bacteriophages are key components of gut microbiomes, yet the phage colonization process in the infant gut remains uncertain. Here, we establish a large phage sequence database and use strain-resolved analyses to investigate DNA phage succession in infants throughout the first 3 years of life. Analysis of 819 fecal metagenomes collected from 28 full-term and 24 preterm infants and their mothers revealed that early-life phageome richness increases over time and reaches adult-like complexity by age 3. Approximately 9% of early phage colonizers, which are mostly maternally transmitted and infect Bacteroides, persist for 3 years and are more prevalent in full-term than in preterm infants. Although rare, phages with stop codon reassignment are more likely to persist than non-recoded phages and generally display an increase in in-frame reassigned stop codons over 3 years. Overall, maternal seeding, stop codon reassignment, host CRISPR-Cas locus prevalence, and diverse phage populations contribute to stable viral colonization.

Megataxonomy and global ecology of the virosphere

Nearly all organisms are hosts to multiple viruses that collectively appear to be the most abundant biological entities in the biosphere. With recent advances in metagenomics and metatranscriptomics, the known diversity of viruses substantially expanded. Comparative analysis of these viruses using advanced computational methods culminated in the reconstruction of the evolution of major groups of viruses and enabled the construction of a virus megataxonomy, which has been formally adopted by the International Committee on Taxonomy of Viruses. This comprehensive taxonomy consists of six virus realms, which are aspired to be monophyletic and assembled based on the conservation of hallmark proteins involved in capsid structure formation or genome replication. The viruses in different major taxa substantially differ in host range and accordingly in ecological niches. In this review article, we outline the latest developments in virus megataxonomy and the recent discoveries that will likely lead to reassessment of some major taxa, in particular, split of three of the current six realms into two or more independent realms. We then discuss the correspondence between virus taxonomy and the distribution of viruses among hosts and ecological niches, as well as the abundance of viruses versus cells in different habitats. The distribution of viruses across environments appears to be primarily determined by the host ranges, i.e. the virome is shaped by the composition of the biome in a given habitat, which itself is affected by abiotic factors.

Driving through stop signs: predicting stop codon reassignment improves functional annotation of bacteriophages

The majority of bacteriophage diversity remains uncharacterized, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and its subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome Catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome median) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased median coding capacity from 66.8% to 90.0% and from 69.0% to 89.8% for UHGV and INPHARED sequences predicted to use translation table 15. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

Structural and Genomic Diversity of Bacteriophages

Bacteriophage diversity is a relatively unknown frontier that is rapidly being explored, leading to a wealth of new information. New bacteriophages are being discovered at an astounding rate via both phage isolation studies and metagenomic analyses. In addition, a nucleotide sequence-based viral taxonomic system has been developed to better handle this wealth of new information. As a result of these developments, phage scientists are transitioning from knowing that there must be huge numbers of diverse kinds of phage particles in natural environments to identifying the actual abundance and phage diversity that is present in specific environments. This review documents the beginning of this transition, offering a glimpse into the magnitude of change unfolding in the field. It stands as a testament to the expanding frontiers of phage research, illuminating the remarkable progress made in unraveling the intricate world of bacteriophage diversity and advancing our understanding of these enigmatic viral entities.

Host translation machinery is not a barrier to phages that interact with both CPR and non-CPR bacteria

IMPORTANCE

Here, we profiled putative phages of Saccharibacteria, which are of particular importance as Saccharibacteria influence some human oral diseases. We additionally profiled putative phages of Gracilibacteria and Absconditabacteria, two Candidate Phyla Radiation (CPR) lineages of interest given their use of an alternative genetic code. Among the phages identified in this study, some are targeted by spacers from both CPR and non-CPR bacteria and others by both bacteria that use the standard genetic code as well as bacteria that use an alternative genetic code. These findings represent new insights into possible phage replication strategies and have relevance for phage therapies that seek to manipulate microbiomes containing CPR bacteria.

Predicting stop codon reassignment improves functional annotation of bacteriophages

The majority of bacteriophage diversity remains uncharacterised, and new intriguing mechanisms of their biology are being continually described. Members of some phage lineages, such as the Crassvirales, repurpose stop codons to encode an amino acid by using alternate genetic codes. Here, we investigated the prevalence of stop codon reassignment in phage genomes and subsequent impacts on functional annotation. We predicted 76 genomes within INPHARED and 712 vOTUs from the Unified Human Gut Virome catalogue (UHGV) that repurpose a stop codon to encode an amino acid. We re-annotated these sequences with modified versions of Pharokka and Prokka, called Pharokka-gv and Prokka-gv, to automatically predict stop codon reassignment prior to annotation. Both tools significantly improved the quality of annotations, with Pharokka-gv performing best. For sequences predicted to repurpose TAG to glutamine (translation table 15), Pharokka-gv increased the median gene length (median of per genome medians) from 287 to 481 bp for UHGV sequences (67.8% increase) and from 318 to 550 bp for INPHARED sequences (72.9% increase). The re-annotation increased mean coding density from 66.8% to 90.0%, and from 69.0% to 89.8% for UHGV and INPHARED sequences. Furthermore, the proportion of genes that could be assigned functional annotation increased, including an increase in the number of major capsid proteins that could be identified. We propose that automatic prediction of stop codon reassignment before annotation is beneficial to downstream viral genomic and metagenomic analyses.

Tentaclins-A Novel Family of Phage Receptor-Binding Proteins That Can Be Hypermutated by DGR Systems

Diversity-generating retroelements (DGRs) are prokaryotic systems providing rapid modification and adaptation of target proteins. In phages, the main targets of DGRs are receptor-binding proteins that are usually parts of tail structures and the variability of such host-recognizing structures enables phage adaptation to changes on the bacterial host surface. Sometimes, more than one target gene containing a hypermutated variable repeat (VR) can be found in phage DGRs. The role of mutagenesis of two functionally different genes is unclear. In this study, several phage genomes that contain DGRs with two target genes were found in the gut virome of healthy volunteers. Bioinformatics analysis of these genes indicated that they encode proteins with different topology; however, both proteins contain the C-type lectin (C-lec) domain with a hypermutated beta-hairpin on its surface. One of the target proteins belongs to a new family of proteins with a specific topology: N-terminal C-lec domain followed by one or more immunoglobulin domains. Proteins from the new family were named tentaclins after TENTACLe + proteIN. The genes encoding such proteins were found in the genomes of prophages and phages from the gut metagenomes. We hypothesized that tentaclins are involved in binding either to bacterial receptors or intestinal/immune cells.

Uncovering the morphological differences between SARS-CoV-2 and SARS-CoV based on transmission electron microscopy images

Comprehending the morphological disparities between SARS-CoV-2 and SARS-CoV viruses can shed light on the underlying mechanisms of infection and facilitate the development of effective diagnostic tools and treatments. Hence, this study aimed to conduct a comprehensive analysis and comparative assessment of the morphology of SARS-CoV-2 and SARS-CoV using transmission electron microscopy (TEM) images. The dataset encompassed 519 isolated SARS-CoV-2 images obtained from patients in Italy (INMI) and 248 isolated SARS-CoV images from patients in Germany (Frankfurt). In this paper, we employed TEM images to scrutinize morphological features, and the outcomes were contrasted with those of SARS-CoV viruses. The findings reveal disparities in the characteristics of SARS-CoV-2 and SARS-CoV, such as envelope protein (E) 98.6 and 102.2 nm, length of spike protein (S) 10.11 and 9.50 nm, roundness 0.86 and 0.88, circularity 0.78 and 0.76, and area sizes 25145.54 and 38591.35 pixels, respectively. In conclusion, these results will augment the identification of virus subtypes, aid in the study of antiviral medications, and enhance our understanding of disease progression and the virus life cycle. Moreover, these findings have the potential to assist in the development of more accurate epidemiological prediction models for COVID-19, leading to better outbreak management and saving lives.

Community Types of the Human Gut Virome are Associated with Endoscopic Outcome in Ulcerative Colitis

BACKGROUND

Inflammatory bowel disease [IBD] is a major debilitating disease. Recently, the gut microbiota has gained attention as an important factor involved in the pathophysiology of IBD. As a complement to the established bacterial 'enterotypes' associated with IBD, we focused here on viruses. We investigated the intestinal virome of IBD patients undergoing biological therapy for the presence of virome configurations associated with IBD, and to uncover how those configurations are associated with therapeutic success.

METHODS

Viral-like particle enrichment followed by deep sequencing was performed on 432 faecal samples from 181 IBD patients starting biological therapy. Redundancy analysis and Dirichlet Multinomial Mixtures were applied to determine covariates of the virome composition and to condense the gut virota into 'viral community types', respectively.

RESULTS

Patients were stratified based on unsupervised clustering into two viral community types. Community type CA showed a low α-diversity and a high relative abundance of Caudoviricetes [non-CrAss] phages and was associated with the dysbiotic Bact2-enterotype. Community type CrM showed a high α-diversity and a high relative abundance of Crassvirales and Malgrandaviricetes phages. During post-interventional analysis, endoscopic outcome was associated with gut virome composition. Remitting UC patients had a high percentage of community type CrM, a high Shannon diversity and a low lysogenic potential. Pre-interventional analyses also identified five novel phages associated with treatment success.

CONCLUSIONS

This study proposed two gut virome configurations that may be involved in the pathophysiology of IBD. Interestingly, those viral configurations are further associated with therapeutic success, suggesting a potential clinical relevance.

Genomes of a Novel Group of Phages That Use Alternative Genetic Code Found in Human Gut Viromes

Metagenomics provides detection of phage genome sequences in various microbial communities. However, the use of alternative genetic codes by some phages precludes the correct analysis of their genomes. In this study, the unusual phage genome (phAss-1, 135,976 bp) was found after the de novo assembly of the human gut virome. Genome analysis revealed the presence of the TAG stop codons in 41 ORFs, including characteristic phage ORFs, and three genes of suppressor tRNA. Comparative analysis indicated that no phages with similar genomes were described. However, two phage genomes (BK046881_ctckW2 and BK025033_ct6IQ4) with substantial similarity to phAss-1 were extracted from the human gut metagenome data. These two complete genomes demonstrated 82.7% and 86.4% of nucleotide identity, respectively, similar genome synteny to phAss-1, the presence of suppressor tRNA genes and suppressor TAG stop codons in many characteristic phage ORFs. These data indicated that phAss-1, BK046881_ctckW2, and BK025033_ct6IQ4 are distinct species within the proposed Phassvirus genus. Moreover, a monophyletic group of divergent phage genomes containing the proposed Phassvirus genus was found among metagenome data. Several phage genomes from the group also contain ORFs with suppressor TAG stop codons, indicating the need to use various translation tables when depositing phage genomes in GenBank.

Viruses in deep-sea cold seep sediments harbor diverse survival mechanisms and remain genetically conserved within species

Deep sea cold seep sediments have been discovered to harbor novel, abundant, and diverse bacterial and archaeal viruses. However, little is known about viral genetic features and evolutionary patterns in these environments. Here, we examined the evolutionary ecology of viruses across active and extinct seep stages in the area of Haima cold seeps in the South China Sea. A total of 338 viral operational taxonomic units are identified and linked to 36 bacterial and archaeal phyla. The dynamics of host-virus interactions are informed by diverse antiviral defense systems across 43 families found in 487 microbial genomes. Cold seep viruses are predicted to harbor diverse adaptive strategies to persist in this environment, including counter-defense systems, auxiliary metabolic genes, reverse transcriptases, and alternative genetic code assignments. Extremely low nucleotide diversity is observed in cold seep viral populations, being influenced by factors including microbial host, sediment depth, and cold seep stage. Most cold seep viral genes are under strong purifying selection with trajectories that differ depending on whether cold seeps are active or extinct. This work sheds light on the understanding of environmental adaptation mechanisms and evolutionary patterns of viruses in the sub-seafloor biosphere.

Factors Affecting Variation of the Human Gut Phageome

The gut microbiome is a dense and metabolically active consortium of microorganisms and viruses located in the lower gastrointestinal tract of the human body. Bacteria and their viruses (phages) are the most abundant members of the gut microbiome. Investigating their biology and the interplay between the two is important if we are to understand their roles in human health and disease. In this review, we summarize recent advances in resolving the taxonomic structure and ecological functions of the complex community of phages in the human gut-the gut phageome. We discuss how age, diet, and geography can all have a significant impact on phageome composition. We note that alterations to the gut phageome have been observed in several diseases such as inflammatory bowel disease, irritable bowel syndrome, and colorectal cancer, and we evaluate whether these phageome changes can directly or indirectly contribute to disease etiology and pathogenesis. We also highlight how lack of standardization in studying the gut phageome has contributed to variation in reported results.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Interrogating the viral dark matter of the rumen ecosystem with a global virome database

The diverse rumen virome can modulate the rumen microbiome, but it remains largely unexplored. Here, we mine 975 published rumen metagenomes for viral sequences, create a global rumen virome database (RVD), and analyze the rumen virome for diversity, virus-host linkages, and potential roles in affecting rumen functions. Containing 397,180 species-level viral operational taxonomic units (vOTUs), RVD substantially increases the detection rate of rumen viruses from metagenomes compared with IMG/VR V3. Most of the classified vOTUs belong to Caudovirales, differing from those found in the human gut. The rumen virome is predicted to infect the core rumen microbiome, including fiber degraders and methanogens, carries diverse auxiliary metabolic genes, and thus likely impacts the rumen ecosystem in both a top-down and a bottom-up manner. RVD and the findings provide useful resources and a baseline framework for future research to investigate how viruses may impact the rumen ecosystem and digestive physiology.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Impact statement

Bacteriophages play a crucial role in shaping microbial communities within the human gut. Among the most dominant bacteriophages in the human gut microbiome are Crassvirales phages, which infect Bacteroides. Despite being widely distributed, only a few Crassvirales genomes have been isolated, leading to a limited understanding of their biology, ecology, and evolution. This study isolated and characterized three novel Crassvirales genomes belonging to two different families, and three genera, but infecting one bacterial host, Bacteroides cellulosilyticus WH2. Notably, the observation confirmed the phages are not co-evolving with their bacterial hosts, rather have a shared ability to exploit similar features in their bacterial host. Additionally, the identification of a critical viral protein undergoing purifying selection and interacting with the bacterial receptors opens doors to targeted therapies against bacterial infections. Given Bacteroides role in polysaccharide degradation in the human gut, our findings advance our understanding of the phage-host interactions and could have important implications for the development of phage-based therapies. These discoveries may hold implications for improving gut health and metabolism to support overall well-being.

Data summary

The genomes used in this research are available on Sequence Read Archive (SRA) within the project, PRJNA737576. Bacteroides cellulosilyticus WH2, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. ' frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11 are all available on GenBank with accessions NZ_CP072251.1 ( B. cellulosilyticus WH2), QQ198717 (Bc01), QQ198718 (Bc03), and QQ198719 (Bc11), and we are working on making the strains available through ATCC. The 3D protein structures for the three Crassvirales genomes are available to download at doi.org/10.25451/flinders.21946034.

Characterization of crAss-like phage isolates highlights Crassvirales genetic heterogeneity and worldwide distribution

Crassvirales (crAss-like phages) are an abundant group of human gut-specific bacteriophages discovered in silico. The use of crAss-like phages as human fecal indicators is proposed but the isolation of only seven cultured strains of crAss-like phages to date has greatly hindered their study. Here, we report the isolation and genetic characterization of 25 new crAss-like phages (termed crAssBcn) infecting Bacteroides intestinalis, belonging to the order Crassvirales, genus Kehishuvirus and, based on their genomic variability, classified into six species. CrAssBcn phage genomes are similar to ΦCrAss001 but show genomic and aminoacidic differences when compared to other crAss-like phages of the same family. CrAssBcn phages are detected in fecal metagenomes around the world at a higher frequency than ΦCrAss001. This study increases the known crAss-like phage isolates and their abundance and heterogeneity open the question of what member of the Crassvirales group should be selected as human fecal marker.

MgCod: Gene Prediction in Phage Genomes with Multiple Genetic Codes

Massive sequencing of microbiomes has led to the discovery of a large number of phage genomes with intermittent stop codon recoding. We have developed a computational tool, MgCod, that identifies genomic regions (blocks) with distinct stop codon recoding simultaneously with the prediction of protein-coding regions. When MgCod was used to scan a large volume of human metagenomic contigs hundreds of viral contigs with intermittent stop codon recoding were revealed. Many of these contigs originated from genomes of known crAssphages. Further analyses had shown that intermittent recoding was associated with subtle patterns in the organization of protein-coding genes, such as 'single-coding' and 'dual-coding'. The dual-coding genes, clustered into blocks, could be translated by two alternative codes producing nearly identical proteins. It was observed that the dual-coded blocks were enriched with the early-stage phage genes, while the late-stage genes were residing in the single-coded blocks. MgCod can identify types of stop codon recoding in novel genomic sequences in parallel with gene prediction. It is available for download from https://github.com/gatech-genemark/MgCod.

The human gut phageome: composition, development, and alterations in disease

The human gastrointestinal tract is colonized by a large number of microorganisms, including bacteria, archaea, viruses, and eukaryotes. The bacterial community has been widely confirmed to have a significant impact on human health, while viruses, particularly phages, have received less attention. Phages are viruses that specifically infect bacteria. They are abundant in the biosphere and exist in a symbiotic relationship with their host bacteria. Although the application of high-throughput sequencing and bioinformatics technology has greatly improved our understanding of the genomic diversity, taxonomic composition, and spatio-temporal dynamics of the human gut phageome, there is still a large portion of sequencing data that is uncharacterized. Preliminary studies have predicted that the phages play a crucial role in driving microbial ecology and evolution. Prior to exploring the function of phages, it is necessary to address the obstacles that hinder establishing a comprehensive sequencing database with sufficient biological properties and understanding the impact of phage-bacteria interactions on human health. In this study, we provide an overview of the human gut phageome, including its composition, structure, and development. We also explore the various factors that may influence the phageome based on current research, including age, diet, ethnicity, and geographical location. Additionally, we summarize the relationship between the phageome and human diseases, such as IBD, IBS, obesity, diabetes, and metabolic syndrome.

The human gut virome: composition, colonization, interactions, and impacts on human health

The gut virome is an incredibly complex part of the gut ecosystem. Gut viruses play a role in many disease states, but it is unknown to what extent the gut virome impacts everyday human health. New experimental and bioinformatic approaches are required to address this knowledge gap. Gut virome colonization begins at birth and is considered unique and stable in adulthood. The stable virome is highly specific to each individual and is modulated by varying factors such as age, diet, disease state, and use of antibiotics. The gut virome primarily comprises bacteriophages, predominantly order Crassvirales, also referred to as crAss-like phages, in industrialized populations and other Caudoviricetes (formerly Caudovirales). The stability of the virome's regular constituents is disrupted by disease. Transferring the fecal microbiome, including its viruses, from a healthy individual can restore the functionality of the gut. It can alleviate symptoms of chronic illnesses such as colitis caused by Clostridiodes difficile. Investigation of the virome is a relatively novel field, with new genetic sequences being published at an increasing rate. A large percentage of unknown sequences, termed 'viral dark matter', is one of the significant challenges facing virologists and bioinformaticians. To address this challenge, strategies include mining publicly available viral datasets, untargeted metagenomic approaches, and utilizing cutting-edge bioinformatic tools to quantify and classify viral species. Here, we review the literature surrounding the gut virome, its establishment, its impact on human health, the methods used to investigate it, and the viral dark matter veiling our understanding of the gut virome.

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.

Gauge your phage: benchmarking of bacteriophage identification tools in metagenomic sequencing data

BACKGROUND

The prediction of bacteriophage sequences in metagenomic datasets has become a topic of considerable interest, leading to the development of many novel bioinformatic tools. A comparative analysis of ten state-of-the-art phage identification tools was performed to inform their usage in microbiome research.

METHODS

Artificial contigs generated from complete RefSeq genomes representing phages, plasmids, and chromosomes, and a previously sequenced mock community containing four phage species, were used to evaluate the precision, recall, and F1 scores of the tools. We also generated a dataset of randomly shuffled sequences to quantify false-positive calls. In addition, a set of previously simulated viromes was used to assess diversity bias in each tool's output.

RESULTS

VIBRANT and VirSorter2 achieved the highest F1 scores (0.93) in the RefSeq artificial contigs dataset, with several other tools also performing well. Kraken2 had the highest F1 score (0.86) in the mock community benchmark by a large margin (0.3 higher than DeepVirFinder in second place), mainly due to its high precision (0.96). Generally, k-mer-based tools performed better than reference similarity tools and gene-based methods. Several tools, most notably PPR-Meta, called a high number of false positives in the randomly shuffled sequences. When analysing the diversity of the genomes that each tool predicted from a virome set, most tools produced a viral genome set that had similar alpha- and beta-diversity patterns to the original population, with Seeker being a notable exception.

CONCLUSIONS

This study provides key metrics used to assess performance of phage detection tools, offers a framework for further comparison of additional viral discovery tools, and discusses optimal strategies for using these tools. We highlight that the choice of tool for identification of phages in metagenomic datasets, as well as their parameters, can bias the results and provide pointers for different use case scenarios. We have also made our benchmarking dataset available for download in order to facilitate future comparisons of phage identification tools. Video Abstract.

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.

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.

Strain-Level Dynamics Reveal Regulatory Roles in Atopic Eczema by Gut Bacterial Phages

The vast population of bacterial phages or viruses (virome) plays pivotal roles in the ecology of human microbial flora and health conditions. Obstacles, including poor viral sequence inference, strain-sensitive virus-host relationship, and the high diversity among individuals, hinder the in-depth understanding of the human virome. We conducted longitudinal studies of the virome based on constructing a high-quality personal reference metagenome (PRM). By applying long-read sequencing for representative samples, we could build a PRM of high continuity that allows accurate annotation and abundance estimation of viruses and bacterial species in all samples of the same individual by aligning short sequencing reads to the PRM. We applied this approach to a series of fecal samples collected for 6 months from a 2-year-old boy who had experienced a 2-month flare-up of atopic eczema (dermatitis) in this period. We identified 31 viral strains in the patient's gut microbiota and deciphered their strain-level relationship to their bacterial hosts. Among them, a lytic crAssphage developed into a dozen substrains and coordinated downregulation in the catabolism of aromatic amino acids (AAAs) in their host bacteria which govern the production of immune-active AAA derivates. The metabolic alterations confirmed based on metabolomic assays cooccurred with symptom remission. Our PRM-based analysis provides an easy approach for deciphering the dynamics of the strain-level human gut virome in the context of entire microbiota. Close temporal correlations among virome alteration, microbial metabolism, and disease remission suggest a potential mechanism for how bacterial phages in microbiota are intimately related to human health. IMPORTANCE The vast populations of viruses or bacteriophages in human gut flora remain mysterious. However, poor annotation and abundance estimation remain obstacles to strain-level analysis and clarification of their roles in microbiome ecology and metabolism associated with human health and diseases. We demonstrate that a personal reference metagenome (PRM)-based approach provides strain-level resolution for analyzing the gut microbiota-associated virome. When applying such an approach to longitudinal samples collected from a 2-year-old boy who has experienced a 2-month flare-up of atopic eczema, we observed thriving substrains of a lytic crAssphage, showing temporal correlation with downregulated catabolism of aromatic amino acids, lower production of immune-active metabolites, and remission of the disease. The PRM-based approach is practical and powerful for strain-centric analysis of the human gut virome, and the underlying mechanism of how strain-level virome dynamics affect disease deserves further investigation.

CrAssphage May Be Viable Markers of Contamination in Pristine and Contaminated River Water

Viruses are the most biologically abundant entities and may be ideal indicators of fecal pollutants in water. Anthropogenic activities have triggered drastic ecosystem changes in rivers, leading to substantial shifts in chemical and biological attributes. Here, we evaluate the viability of using the presence of crAssphage as indicators of fecal contamination in South African rivers. Shotgun analysis revealed diverse crAssphage viruses in these rivers, which are impacted by chemical and biological pollution. Overall, the diversity and relative abundances of these viruses was higher in contaminated sites compared to pristine locations. In contrast to fecal coliform counts, crAssphage sequences were detected in pristine rivers, supporting the assertion that the afore mentioned marker may be a more accurate indicator of fecal contamination. Our data demonstrate the presence of diverse putative hosts which includes members of the phyla Bacteroidota, Pseudomonadota, Verrucomicrobiota, and Bacillota. Phylogenetic analysis revealed novel subfamilies, suggesting that rivers potentially harbor distinct and uncharacterized clades of crAssphage. These data provide the first insights regarding the diversity, distribution, and functional roles of crAssphage in rivers. Taken together, the results support the potential application of crAssphage as viable markers for water quality monitoring. IMPORTANCE Rivers support substantial populations and provide important ecosystem services. Despite the application of fecal coliform tests and other markers, we lack rapid and reproducible approaches for determining fecal contamination in rivers. Waterborne viral outbreaks have been reported even after fecal indicator bacteria (FIB) were suggested to be absent or below regulated levels of coliforms. This indicates a need to develop and apply improved indicators of pollutants in aquatic ecosystems. Here, we evaluate the viability of crAssphage as indicators of fecal contamination in two South African rivers. We assess the abundance, distribution, and diversity of these viruses in sites that had been predicted pristine or contaminated by FIB analysis. We show that crAssphage are ideal and sensitive markers for fecal contamination and describe novel clades of crAss-like phages. Known crAss-like subfamilies were unrepresented in our data, suggesting that the diversity of these viruses may reflect geographic locality and dependence.

Alterations of the gut virome in patients with systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that has been linked to the dysbiosis of the gut microbiome and virome. However, the potential characterization of the gut virome in SLE patients needs to be explored more extensively.

Methods

Herein, we analyzed the gut viral community of 16 SLE patients and 31 healthy controls using both bulk and virus-like particle (VLP)-based metagenomic sequencing of their fecal samples. A total of 15,999 non-redundant viral operational taxonomic units (vOTUs) were identified from the metagenomic assembled contigs and used for gut virome profiling.

Results

SLE patients exhibited a significant decrease in gut viral diversity in the bulk metagenome dataset, but this change was not significant in the VLP metagenome dataset. Also, considerable alterations of the overall gut virome composition and remarkable changes in the viral family compositions were observed in SLE patients compared with healthy controls, as observed in both two technologies. We identified 408 vOTUs (177 SLE-enriched and 231 control-enriched) with significantly different relative abundances between patients and controls in the bulk virome, and 18 vOTUs (17 SLE-enriched in 1 control-enriched) in the VLP virome. The SLE-enriched vOTUs included numerous Siphoviridae, Microviridae, and crAss-like viruses and were frequently predicted to infect Bacteroides, Parabacteroides, and Ruminococcus_E, while the control-enriched contained numerous members of Siphoviridae and Myoviridae and were predicted to infect Prevotella and Lachnospirales_CAG-274. We explored the correlations between gut viruses and bacteria and found that some Lachnospirales_CAG-274 and Hungatella_A phages may play key roles in the virus-bacterium network. Furthermore, we explored the gut viral signatures for disease discrimination and achieved an area under the receiver operator characteristic curve (AUC) of above 0.95, suggesting the potential of the gut virome in the prediction of SLE.

Conclusion

Our findings demonstrated the alterations in viral diversity and taxonomic composition of the gut virome of SLE patients. Further research into the etiology of SLE and the gut viral community will open up new avenues for treating and preventing SLE and other autoimmune diseases.

Gut bacteriome, mycobiome and virome alterations in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic destructive autoimmune disease of the joints which causes significant pain, functional disability, and mortality. Although aberrant immune cell activation induced by the imbalance between T helper Th1/Th17 and Treg cells is implicated in the RA development, its etiopathogenesis remains unclear. The presence of mucosal inflammation and systemic IgA-isotype-autoantibodies (anti-citrullinated peptide antibodies and rheumatoid factor) in pre-clinical RA supports the mucosal origin hypothesis involving altered microbiota in disease development. The gut microbiota comprises diverse bacteria, fungal and viral components, which are critical in developing host immunity. Alterations in microbial abundance are known to exacerbate or attenuate immune responses in the gut microenvironment subsequently affecting the joints. Further, these changes can provide biomarkers for disease activity and outcome in RA. Most of the research till date has been focused on describing gut bacterial components in RA. Studies on gut mycobiome and virome components in RA are relatively new and burgeoning field. Given the paucity of mycobiome or virome specific studies in RA, this review, discusses the recent findings on alterations in gut bacterial, fungal, and viral components as well as their role in regulating the spectrum of immune-pathogenic events occurring in RA which might be explored in future as a potential therapeutic target. Further, we provide an overview on inter-kingdom interactions between bacteria, fungi, and viruses in RA. The current understanding on gut microbiota modulation for managing RA is also summarised.

Phage-encoded carbohydrate-interacting proteins in the human gut

In the human gastrointestinal tract, the gut mucosa and the bacterial component of the microbiota interact and modulate each other to accomplish a variety of critical functions. These include digestion aid, maintenance of the mucosal barrier, immune regulation, and production of vitamins, hormones, and other metabolites that are important for our health. The mucus lining of the gut is primarily composed of mucins, large glycosylated proteins with glycosylation patterns that vary depending on factors including location in the digestive tract and the local microbial population. Many gut bacteria have evolved to reside within the mucus layer and thus encode mucus-adhering and -degrading proteins. By doing so, they can influence the integrity of the mucus barrier and therefore promote either health maintenance or the onset and progression of some diseases. The viral members of the gut - mostly composed of bacteriophages - have also been shown to have mucus-interacting capabilities, but their mechanisms and effects remain largely unexplored. In this review, we discuss the role of bacteriophages in influencing mucosal integrity, indirectly via interactions with other members of the gut microbiota, or directly with the gut mucus via phage-encoded carbohydrate-interacting proteins. We additionally discuss how these phage-mucus interactions may influence health and disease states.

Pharokka: a fast scalable bacteriophage annotation tool

SUMMARY

In recent years, there has been an increasing interest in bacteriophages, which has led to growing numbers of bacteriophage genomic sequences becoming available. Consequently, there is a need for a rapid and consistent genomic annotation tool dedicated for bacteriophages. Existing tools either are not designed specifically for bacteriophages or are web- and email-based and require significant manual curation, which makes their integration into bioinformatic pipelines challenging. Pharokka was created to provide a tool that annotates bacteriophage genomes easily, rapidly and consistently with standards compliant outputs. Moreover, Pharokka requires only two lines of code to install and use and takes under 5 min to run for an average 50-kb bacteriophage genome.

AVAILABILITY AND IMPLEMENTATION

Pharokka is implemented in Python and is available as a bioconda package using 'conda install -c bioconda pharokka'. The source code is available on GitHub (https://github.com/gbouras13/pharokka). Pharokka has been tested on Linux-64 and MacOSX machines and on Windows using a Linux Virtual Machine.

Multiomic spatial analysis reveals a distinct mucosa-associated virome

The human gut virome has been increasingly explored in recent years. However, nearly all virome-sequencing efforts rely solely on fecal samples and few studies leverage multiomic approaches to investigate phage-host relationships. Here, we combine metagenomics, metaviromics, and metatranscriptomics to study virome-bacteriome interactions at the colonic mucosal-luminal interface in a cohort of three individuals with inflammatory bowel disease; non-IBD controls were not included in this study. We show that the mucosal viral population is distinct from the stool virome and houses abundant crAss-like phages that are undetectable by fecal sampling. Through viral protein prediction and metatranscriptomic analysis, we explore viral gene transcription, prophage activation, and the relationship between the presence of integrase and temperate phages in IBD subjects. We also show the impact of deep sequencing on virus recovery and offer guidelines for selecting optimal sequencing depths in future metaviromic studies. Systems biology approaches such as those presented in this report will enhance our understanding of the human virome and its interactions with our microbiome and our health.

The global virome: How much diversity and how many independent origins?

Viruses are considered to be the most abundant biological entities on earth. They also display striking genetic diversity as emphatically demonstrated by the recent advances of metagenomics and metatranscriptomics. But what are the limits of this diversity, that is, how many virus species in the earth virome? By combining the available estimates of the number of prokaryote species with those of the virome size, we obtain back-of-the-envelope estimates of the total number of distinct virus species, which come out astronomically large, from about 10⁷ to about 10⁹ . The route of virus origins apparently involved non-viral replicators capturing and exapting various cellular proteins to become virus capsid subunits. How many times in the history of life has this happened? In other words, how many realms of viruses, the highest rank taxa that are supposed to be monophyletic, comprise the global virome? We argue that viruses emerged on a number (even if far from astronomical) independent occasions, so the number of realms will considerably increase from the current 6, by splitting some of the current realms, giving the realm status to some of the currently unclassified groups of viruses and discovery of new distinct groups.

Prokaryotic and viral genomes recovered from 787 Japanese gut metagenomes revealed microbial features linked to diets, populations, and diseases

We reconstructed 19,084 prokaryotic and 31,395 viral genomes from 787 Japanese gut metagenomes as Japanese metagenome-assembled genomes (JMAG) and Japanese Virus Database (JVD), which are large microbial genome datasets for a single population. Population-specific enrichment of the Bacillus subtilis and β-porphyranase among the JMAG could derive from the Japanese traditional food natto (fermented soybeans) and nori (laver), respectively. Dairy-related Enterococcus_B lactis and Streptococcus thermophilus were nominally associated with the East Asian-specific missense variant rs671:G>A in ALDH2, which was associated with dairy consumption. Of the species-level viral genome clusters in the JVD, 62.9% were novel. The β crAss-like phage composition was low among the Japanese but relatively high among African and Oceanian peoples. Evaluations of the association between crAss-like phages and diseases showed significant disease-specific associations. Our large catalog of virus-host pairs identified the positive correlation between the abundance of the viruses and their hosts.