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

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Cyanobacteria-cyanophage interactions between freshwater and marine ecosystems based on large-scale cyanophage genomic analysis

The disparities in harmful algal blooms dynamics are largely attributed to variations in cyanobacteria populations within aquatic ecosystems. However, cyanobacteria-cyanophage interactions and their role in shaping cyanobacterial populations has been previously underappreciated. To address this knowledge gap, we isolated and sequenced 42 cyanophages from diverse water sources in China, with the majority (n = 35) originating from freshwater sources. We designated these sequences as the "Novel Cyanophage Genome sequence Collection" (NCGC). NCGC displayed notable genetic variations, with 95 % (40/42) of the sequences representing previously unidentified taxonomic ranks. By integrating NCGC with public data of cyanophages and cyanobacteria, we found evidence for more frequent historical cyanobacteria-cyanophage interactions in freshwater ecosystems. This was evidenced by a higher prevalence of prophage integrase-related genes in freshwater cyanophages (37.97 %) than marine cyanophages (7.42 %). In addition, freshwater cyanophages could infect a broader range of cyanobacteria orders (n = 4) than marine ones (n = 0). Correspondingly, freshwater cyanobacteria harbored more defense systems per million base pairs in their genomes, indicating more frequent phage infections. Evolutionary and cyanophage epidemiological studies suggest that interactions between cyanobacteria and cyanophages in freshwater and marine ecosystems are interconnected, and that brackish water can act as a transitional zone for freshwater and marine cyanophages. In conclusion, our research significantly expands the genetic information database of cyanophage, offering a wider selection of cyanophages to control harmful cyanobacterial blooms. Additionally, we represent a pioneering large-scale and comprehensive analysis of cyanobacteria and cyanophage sequencing data, and it provides theoretical guidance for the application of cyanophages in different environments.

Rheumatoid arthritis patients harbour aberrant enteric bacteriophages with autoimmunity-provoking potential: a paired sibling study

OBJECTIVES

Viruses have been considered as important participants in the development of rheumatoid arthritis (RA). However, the profile of enteric virome and its role in RA remains elusive. This study aimed to investigate the atlas and involvement of virome in RA pathogenesis.

METHODS

Faecal samples from 30 pairs of RA and healthy siblings that minimise genetic interferences were collected for metagenomic sequencing. The α and β diversity of the virome and the virome-bacteriome interaction were analysed. The differential bacteriophages were identified, and their correlations with clinical and immunological features of RA were analysed. The potential involvement of these differential bacteriophages in RA pathogenesis was further investigated by auxiliary metabolic gene annotation and molecular mimicry study. The responses of CD4+ T cells and B cells to the mimotopes derived from the differential bacteriophages were systemically studied.

RESULTS

The composition of the enteric bacteriophageome was distorted in RA. The differentially presented bacteriophages correlated with the immunological features of RA, including anti-CCP autoantibody and HLA-DR shared epitope. Intriguingly, the glycerolipid and purine metabolic genes were highly active in the bacteriophages from RA. Moreover, peptides of RA-enriched phages, in particular Prevotella phage and Oscillibacter phage could provoke the autoimmune responses in CD4+ T cells and plasma cells via molecular mimicry of the disease-associated autoantigen epitopes, especially those of Bip.

CONCLUSIONS

This study provides new insights into enteric bacteriophageome in RA development. In particular, the aberrant bacteriophages demonstrated autoimmunity-provoking potential that would promote the occurrence of the disease.

Tools and methodology to in silico phage discovery in freshwater environments

Freshwater availability is essential, and its maintenance has become an enormous challenge. Due to population growth and climate changes, freshwater sources are becoming scarce, imposing the need for strategies for its reuse. Currently, the constant discharge of waste into water bodies from human activities leads to the dissemination of pathogenic bacteria, negatively impacting water quality from the source to the infrastructure required for treatment, such as the accumulation of biofilms. Current water treatment methods cannot keep pace with bacterial evolution, which increasingly exhibits a profile of multidrug resistance to antibiotics. Furthermore, using more powerful disinfectants may affect the balance of aquatic ecosystems. Therefore, there is a need to explore sustainable ways to control the spreading of pathogenic bacteria. Bacteriophages can infect bacteria and archaea, hijacking their host machinery to favor their replication. They are widely abundant globally and provide a biological alternative to bacterial treatment with antibiotics. In contrast to common disinfectants and antibiotics, bacteriophages are highly specific, minimizing adverse effects on aquatic microbial communities and offering a lower cost-benefit ratio in production compared to antibiotics. However, due to the difficulty involving cultivating and identifying environmental bacteriophages, alternative approaches using NGS metagenomics in combination with some bioinformatic tools can help identify new bacteriophages that can be useful as an alternative treatment against resistant bacteria. In this review, we discuss advances in exploring the virome of freshwater, as well as current applications of bacteriophages in freshwater treatment, along with current challenges and future perspectives.

Characterization of CRISPR-Cas Systems in Shewanella algae and Shewanella haliotis: Insights into the Adaptation and Survival of Marine Pathogens

CRISPR-Cas systems are adaptive immune mechanisms present in most prokaryotes that play an important role in the adaptation of bacteria and archaea to new environments. Shewanella algae is a marine zoonotic pathogen with worldwide distribution, which accounts for the majority of clinical cases of Shewanella infections. However, the characterization of Shewanella algae CRISPR-Cas systems has not been well investigated yet. Through whole genome sequence analysis, we characterized the CRISPR-Cas systems in S. algae. Our results indicate that CRISPR-Cas systems are prevalent in S. algae, with the majority of strains containing the Type I-F system. This study provides new insights into the diversity and function of CRISPR-Cas systems in S. algae and highlights their potential role in the adaptation and survival of these marine pathogens.

Profiles of phage in global hospital wastewater: Association with microbial hosts, antibiotic resistance genes, metal resistance genes, and mobile genetic elements

Hospital wastewater (HWW) is known to host taxonomically diverse microbial communities, yet limited information is available on the phages infecting these microorganisms. To fill this knowledge gap, we conducted an in-depth analysis using 377 publicly available HWW metagenomic datasets from 16 countries across 4 continents in the NCBI SRA database to elucidate phage-host dynamics and phage contributions to resistance gene transmission. We first assembled a metagenomic HWW phage catalog comprising 13,812 phage operational taxonomic units (pOTUs). The majority of these pOTUs belonged to the Caudoviricetes order, representing 75.29 % of this catalog. Based on the lifestyle of phages, we found that potentially virulent phages predominated in HWW. Specifically, 583 pOTUs have been predicted to have the capability to lyse 81 potentially pathogenic bacteria, suggesting the promising role of HWW phages as a viable alternative to antibiotics. Among all pOTUs, 1.56 % of pOTUs carry 108 subtypes of antibiotic resistance genes (ARGs), 0.96 % of pOTUs carry 76 subtypes of metal resistance genes (MRGs), and 0.96 % of pOTUs carry 22 subtypes of non-phage mobile genetic elements (MGEs). Predictions indicate that certain phages carrying ARGs, MRGs, and non-phage MGEs could infect bacteria hosts, even potential pathogens. This suggests that phages in HWW may contribute to the dissemination of resistance-associated genes in the environment. This meta-analysis provides the first global catalog of HWW phages, revealing their correlations with microbial hosts and pahge-associated ARGs, MRG, and non-phage MGEs. The insights gained from this research hold promise for advancing the applications of phages in medical and industrial contexts.

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.

Efficient Recovery of Complete Gut Viral Genomes by Combined Short- and Long-Read Sequencing

Current metagenome assembled human gut phage catalogs contained mostly fragmented genomes. Here, comprehensive gut virome detection procedure is developed involving virus-like particle (VLP) enrichment from ≈500 g feces and combined sequencing of short- and long-read. Applied to 135 samples, a Chinese Gut Virome Catalog (CHGV) is assembled consisting of 21,499 non-redundant viral operational taxonomic units (vOTUs) that are significantly longer than those obtained by short-read sequencing and contained ≈35% (7675) complete genomes, which is ≈nine times more than those in the Gut Virome Database (GVD, ≈4%, 1,443). Interestingly, the majority (≈60%, 13,356) of the CHGV vOTUs are obtained by either long-read or hybrid assemblies, with little overlap with those assembled from only the short-read data. With this dataset, vast diversity of the gut virome is elucidated, including the identification of 32% (6,962) novel vOTUs compare to public gut virome databases, dozens of phages that are more prevalent than the crAssphages and/or Gubaphages, and several viral clades that are more diverse than the two. Finally, the functional capacities are also characterized of the CHGV encoded proteins and constructed a viral-host interaction network to facilitate future research and applications.

Advances in phage-host interaction prediction: in silico method enhances the development of phage therapies

Phages can specifically recognize and kill bacteria, which lead to important application value of bacteriophage in bacterial identification and typing, livestock aquaculture and treatment of human bacterial infection. Considering the variety of human-infected bacteria and the continuous discovery of numerous pathogenic bacteria, screening suitable therapeutic phages that are capable of infecting pathogens from massive phage databases has been a principal step in phage therapy design. Experimental methods to identify phage-host interaction (PHI) are time-consuming and expensive; high-throughput computational method to predict PHI is therefore a potential substitute. Here, we systemically review bioinformatic methods for predicting PHI, introduce reference databases and in silico models applied in these methods and highlight the strengths and challenges of current tools. Finally, we discuss the application scope and future research direction of computational prediction methods, which contribute to the performance improvement of prediction models and the development of personalized phage therapy.

Adaptive strategies and ecological roles of phages in habitats under physicochemical stress

Bacteriophages (phages) play a vital role in ecosystem functions by influencing the composition, genetic exchange, metabolism, and environmental adaptation of microbial communities. With recent advances in sequencing technologies and bioinformatics, our understanding of the ecology and evolution of phages in stressful environments has substantially expanded. Here, we review the impact of physicochemical environmental stress on the physiological state and community dynamics of phages, the adaptive strategies that phages employ to cope with environmental stress, and the ecological effects of phage-host interactions in stressful environments. Specifically, we highlight the contributions of phages to the adaptive evolution and functioning of microbiomes and suggest that phages and their hosts can maintain a mutualistic relationship in response to environmental stress. In addition, we discuss the ecological consequences caused by phages in stressful environments, encompassing biogeochemical cycling. Overall, this review advances an understanding of phage ecology in stressful environments, which could inform phage-based strategies to improve microbiome performance and ecosystem resilience and resistance in natural and engineering systems.

Forty-nine metagenomic-assembled genomes from an aquatic virome expand Caudoviricetes by 45 potential new families and the newly uncovered Gossevirus of Bamfordvirae

Twenty complete genomes (29-63 kb) and 29 genomes with an estimated completeness of over 90 % (30-90 kb) were identified for novel dsDNA viruses in the Yangshan Harbor metavirome. These newly discovered viruses contribute to the expansion of viral taxonomy by introducing 46 potential new families. Except for one virus, all others belong to the class Caudoviricetes. The exception is a novel member of the recently characterized viral group known as Gossevirus. Fifteen viruses were predicted to be temperate. The predicted hosts for the viruses appear to be involved in various aspects of the nitrogen cycle, including nitrogen fixation, oxidation and denitrification. Two viruses were identified to have a host of Flavobacterium and Tepidimonas fonticaldi, respectively, by matching CRISPR spacers with viral protospacers. Our findings provide an overview for characterizing and identifying specific viruses from Yangshan Harbor. The Gossevirus-like virus uncovered emphasizes the need for further comprehensive isolation and investigation of polinton-like viruses.

Metagenomic analysis of the honey bee queen microbiome reveals low bacterial diversity and Caudoviricetes phages

In eusocial insects, the health of the queens-the colony founders and sole reproductive females-is a primary determinant for colony success. Queen failure in the honey bee Apis mellifera, for example, is a major concern of beekeepers who annually suffer colony losses, necessitating a greater knowledge of queen health. Several studies on the microbiome of honey bees have characterized its diversity and shown its importance for the health of worker bees, the female non-reproductive caste. However, the microbiome of workers differs from that of queens, which, in comparison, is still poorly studied. Thus, direct investigations of the queen microbiome are required to understand colony-level microbiome assembly, functional roles, and evolution. Here, we used metagenomics to comprehensively characterize the honey bee queen microbiome. Comparing samples from different geographic locations and breeder sources, we show that the microbiome of queens is mostly shaped by the environment experienced since early life and is predicted to play roles in the breakdown of the diet and protection from pathogens and xenobiotics. We also reveal that the microbiome of queens comprises only four candidate core bacterial species, Apilactobacillus kunkeei, Lactobacillus apis, Bombella apis, and Commensalibacter sp. Interestingly, in addition to bacteria, we show that bacteriophages infect the queen microbiome, for which Lactobacillaceae are predicted to be the main reservoirs. Together, our results provide the basis to understand the honey bee colony microbiome assemblage, can guide improvements in queen-rearing processes, and highlight the importance of considering bacteriophages for queen microbiome health and microbiome homeostasis in eusocial insects.IMPORTANCEThe queen caste plays a central role in colony success in eusocial insects, as queens lay eggs and regulate colony behavior and development. Queen failure can cause colonies to collapse, which is one of the major concerns of beekeepers. Thus, understanding the biology behind the queen's health is a pressing issue. Previous studies have shown that the bee microbiome plays an important role in worker bee health, but little is known about the queen microbiome and its function in vivo. Here, we characterized the queen microbiome, identifying for the first time the present species and their putative functions. We show that the queen microbiome has predicted nutritional and protective roles in queen association and comprises only four consistently present bacterial species. Additionally, we bring to attention the spread of phages in the queen microbiome, which increased in abundance in failing queens and may impact the fate of the colony.

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

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

Biological and bioinformatic tools for the discovery of unknown phage-host combinations

The field of microbial ecology has been transformed by metagenomics in recent decades and has culminated in vast datasets that facilitate the bioinformatic dissection of complex microbial communities. Recently, attention has turned from defining the microbiota composition to the interactions and relationships that occur between members of the microbiota. Within complex microbiota, the identification of bacteriophage-host combinations has been a major challenge. Recent developments in artificial intelligence tools to predict protein structure and function as well as the relationships between bacteria and their infecting bacteriophages allow a strategic approach to identifying and validating phage-host relationships. However, biological validation of these predictions remains essential and will serve to improve the existing predictive tools. In this review, I provide an overview of the most recent developments in both bioinformatic and experimental approaches to predicting and experimentally validating unknown phage-host combinations.

Hiding in plain sight: The discovery of complete genomes of 11 hypothetical spindle-shaped viruses that putatively infect mesophilic ammonia-oxidizing archaea

The genome of a putative Nitrosopumilaceae virus with a hypothetical spindle-shaped particle morphology was identified in the Yangshan Harbour metavirome from the East China Sea through protein similarity comparison and structure analysis. This discovery was accompanied by a set of 10 geographically dispersed close relatives found in the environmental virus datasets from typical locations of ammonia-oxidizing archaeon distribution. Its host prediction was supported by iPHoP prediction and protein sequence similarity. The structure of the predicted major capsid protein, together with the overall N-glycosylation site, the transmembrane helices prediction, the hydrophilicity profile, and the docking simulation of the major capsid proteins, indicate that these viruses resemble spindle-shaped viruses. It suggests a similarly assembled structure and, consequently, a possibly spindle-shaped morphology of these newly discovered archaeal viruses.

Tailed bacteriophages (Caudoviricetes) dominate the microbiome of a diseased stingless bee

Bacteriophages, viruses that infect bacterial hosts, are known to rule the dynamics and diversity of bacterial populations in a number of ecosystems. Bacterial communities residing in the gut of animals, known as the gut microbiome, have revolutionized our understanding of many diseases. However, the gut phageome, while of apparent importance in this context, remains an underexplored area of research. Here we identify for the first time genomic sequences from tailed viruses (Caudoviricetes) that are associated with the microbiome of stingless bees (Melipona quadrifasciata). Both DNA and RNA were extracted from virus particles isolated from healthy and diseased forager bees, the latter showing symptoms from an annual syndrome that only affects M. quadrifasciata. Viral contigs from previously sequenced metagenomes of healthy and diseased forager bees were used for the analyses. Using conserved proteins deduced from their genomes, we found that Caudoviricetes were only present in the worker bee gut microbiome from diseased stingless bees. The most abundant phages are phylogenetically related to phages that infect Gram-positive bacteria from the order Lactobacillales and Gram-negative bacteria from the genus Gilliamella and Bartonella, that are common honey bee symbionts. The potential implication of these viruses in the M. quadrifasciata syndrome is discussed.

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.

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

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

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

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

Prediction and characterization of prophages of Stenotrophomonas maltophilia reveals a remarkable phylogenetic diversity of prophages

Prophages, which enables bacterial hosts to acquire novel traits, and increase genetic variation and evolutionary innovation, are considered to be one of the greatest drivers of bacterial diversity and evolution. Stenotrophomonas maltophilia is widely distributed and one of the most important multidrug resistant bacteria in hospitals. However, the distribution and genetic diversity of S. maltophilia prophages have not been elucidated. In this study, putative prophages were predicted in S. maltophilia genomes by using virus prediction tools, and the genetic diversity and phylogeny of S. maltophilia and the prophages they harbor were further analyzed. A total of 356 prophage regions were predicted from 88 S. maltophilia genomes. Among them, 144 were intact prophages, but 77.09% of the intact prophages did not match any known phage sequences in the public database. The number of prophage carried by S. maltophilia is related to its host habitat and is an important factor affecting the size of the host genome, but it is not related to the genetic diversity of the prophage. The prediction of auxiliary genes encoded by prophage showed that antibiotic resistance genes was not predicted for any of the prophages except for one questionable prophage, while 53 virulence genes and 169 carbohydrate active enzymes were predicted from 11.24 and 44.1% prophages, respectively. Most of the prophages (72.29%) mediated horizontal gene transfer of S. maltophilia genome, but only involved in 6.25% of the horizontal gene transfer events. In addition, CRISPR prediction indicated 97.75% S. maltophilia strains contained the CRISPR-Cas system containing 818 spacer sequences. However, these spacer sequences did not match any known S. maltophilia phages, and only a few S. maltophilia prophages. Comparative genomic analysis revealed a highly conserved and syntenic organization with genomic rearrangement between the prophages and the known related S. maltophilia phages. Our results indicate a high prevalence and genetic diversity of prophages in the genome of S. maltophilia, as well as the presence of a large number of uncharacterized phages. It provides an important complement to understanding the diversity and biological characteristics of phages, as well as the interactions and evolution between bacteria and phages.

Viral metagenomics reveals diverse virus-host interactions throughout the soil depth profile

IMPORTANCE

Soil viruses can moderate the roles that their host microbes play in global carbon cycling. However, given that most studies investigate the surface layer (i.e., top 20 cm) of soil, the extent to which this occurs in subsurface soil (i.e., below 20 cm) is unknown. Here, we leveraged public sequencing data to investigate the interactions between viruses and their hosts at soil depth intervals, down to 115 cm. While most viruses were detected throughout the soil depth profile, their adaptation to host microbes varied. Nonetheless, we uncovered evidence for the potential of soil viruses to encourage their hosts to recycle plant-derived carbon in both surface and subsurface soils. This work reasons that our understanding of soil viral functions requires us to continue to dig deeper and compare viruses existing throughout soil ecosystems.

Novel virulent and temperate cyanophages predicted to infect Microcoleus associated with anatoxin-producing benthic mats

Cyanophages are crucial for regulating cyanobacterial populations, but their influence on anatoxin-producing Microcoleus mat dynamics remains unexplored. Here, we use metagenomics to explore phage presence in benthic mats from the Wolastoq|Saint John River (New Brunswick, Canada) and the Eel River (California, USA). We recovered multiple viral-like sequences associated with different putative bacterial hosts, including two cyanophage genomes with apparently different replication strategies. A temperate cyanophage was found integrated in the genomes of Microcoleus sp. 3 recovered from the Eel River and is phylogenetically related to Phormidium phages. We also recovered novel virulent cyanophage genomes from Wolastoq and Eel River mats that were dominated by anatoxin-producing Microcoleus species predicted to be the host. Despite the geographical distance, these genomes have similar sizes (circa 239 kbp) and share numerous orthologous genes with high sequence identity. A considerable reduction of the anatoxin-producing Microcoleus species in Wolastoq mats following the emergence of the virulent phage suggests that phage infections have an important role in limiting the abundance of this toxigenic cyanobacterium and releasing anatoxins into the surrounding water. Our results constitute the first report of cyanophages predicted to infect mat-forming Microcoleus species associated with anatoxin production.

Genomic Analysis and Taxonomic Characterization of Seven Bacteriophage Genomes Metagenomic-Assembled from the Dishui Lake

Viruses in aquatic ecosystems exhibit remarkable abundance and diversity. However, scattered studies have been conducted to mine uncultured viruses and identify them taxonomically in lake water. Here, whole genomes (29-173 kbp) of seven uncultured dsDNA bacteriophages were discovered in Dishui Lake, the largest artificial lake in Shanghai. We analyzed their genomic signatures and found a series of viral auxiliary metabolic genes closely associated with protein synthesis and host metabolism. Dishui Lake phages shared more genes with uncultivated environmental viruses than with reference viruses based on the gene-sharing network classification. Phylogeny of proteomes and comparative genomics delineated three new genera within two known viral families of Kyanoviridae and Autographiviridae, and four new families in Caudoviricetes for these seven novel phages. Their potential hosts appeared to be from the dominant bacterial phyla in Dishui Lake. Altogether, our study provides initial insights into the composition and diversity of bacteriophage communities in Dishui Lake, contributing valuable knowledge to the ongoing research on the roles played by viruses in freshwater ecosystems.

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

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

Phages are unrecognized players in the ecology of the oral pathogen Porphyromonas gingivalis

BACKGROUND

Porphyromonas gingivalis (hereafter "Pg") is an oral pathogen that has been hypothesized to act as a keystone driver of inflammation and periodontal disease. Although Pg is most readily recovered from individuals with actively progressing periodontal disease, healthy individuals and those with stable non-progressing disease are also colonized by Pg. Insights into the factors shaping the striking strain-level variation in Pg, and its variable associations with disease, are needed to achieve a more mechanistic understanding of periodontal disease and its progression. One of the key forces often shaping strain-level diversity in microbial communities is infection of bacteria by their viral (phage) predators and symbionts. Surprisingly, although Pg has been the subject of study for over 40 years, essentially nothing is known of its phages, and the prevailing paradigm is that phages are not important in the ecology of Pg.

RESULTS

Here we systematically addressed the question of whether Pg are infected by phages-and we found that they are. We found that prophages are common in Pg, they are genomically diverse, and they encode genes that have the potential to alter Pg physiology and interactions. We found that phages represent unrecognized targets of the prevalent CRISPR-Cas defense systems in Pg, and that Pg strains encode numerous additional mechanistically diverse candidate anti-phage defense systems. We also found that phages and candidate anti-phage defense system elements together are major contributors to strain-level diversity and the species pangenome of this oral pathogen. Finally, we demonstrate that prophages harbored by a model Pg strain are active in culture, producing extracellular viral particles in broth cultures.

CONCLUSION

This work definitively establishes that phages are a major unrecognized force shaping the ecology and intra-species strain-level diversity of the well-studied oral pathogen Pg. The foundational phage sequence datasets and model systems that we establish here add to the rich context of all that is already known about Pg, and point to numerous avenues of future inquiry that promise to shed new light on fundamental features of phage impacts on human health and disease broadly. Video Abstract.

Characterization of prophages in bacterial genomes from the honey bee (Apis mellifera) gut microbiome

The gut of the European honey bee (Apis mellifera) possesses a relatively simple bacterial community, but little is known about its community of prophages (temperate bacteriophages integrated into the bacterial genome). Although prophages may eventually begin replicating and kill their bacterial hosts, they can also sometimes be beneficial for their hosts by conferring protection from other phage infections or encoding genes in metabolic pathways and for toxins. In this study, we explored prophages in 17 species of core bacteria in the honey bee gut and two honey bee pathogens. Out of the 181 genomes examined, 431 putative prophage regions were predicted. Among core gut bacteria, the number of prophages per genome ranged from zero to seven and prophage composition (the compositional percentage of each bacterial genome attributable to prophages) ranged from 0 to 7%. Snodgrassella alvi and Gilliamella apicola had the highest median prophages per genome (3.0 ± 1.46; 3.0 ± 1.59), as well as the highest prophage composition (2.58% ± 1.4; 3.0% ± 1.59). The pathogen Paenibacillus larvae had a higher median number of prophages (8.0 ± 5.33) and prophage composition (6.40% ± 3.08) than the pathogen Melissococcus plutonius or any of the core bacteria. Prophage populations were highly specific to their bacterial host species, suggesting most prophages were acquired recently relative to the divergence of these bacterial groups. Furthermore, functional annotation of the predicted genes encoded within the prophage regions indicates that some prophages in the honey bee gut encode additional benefits to their bacterial hosts, such as genes in carbohydrate metabolism. Collectively, this survey suggests that prophages within the honey bee gut may contribute to the maintenance and stability of the honey bee gut microbiome and potentially modulate specific members of the bacterial community, particularly S. alvi and G. apicola.

Large-scale phage cultivation for commensal human gut bacteria

Phages are highly abundant in the human gut, yet most of them remain uncultured. Here, we present a gut phage isolate collection (GPIC) containing 209 phages for 42 commensal human gut bacterial species. Genome analysis of the phages identified 34 undescribed genera. We discovered 22 phages from the Salasmaviridae family that have small genomes (∼10-20 kbp) and infect Gram-positive bacteria. Two phages from a candidate family, Paboviridae, with high prevalence in the human gut were also identified. Infection assays showed that Bacteroides and Parabacteroides phages are specific to a bacterial species, and strains of the same species also exhibit substantial variations in phage susceptibility. A cocktail of 8 phages with a broad host range for Bacteroides fragilis strains effectively reduced their abundance in complex host-derived communities in vitro. Our study expands the diversity of cultured human gut bacterial phages and provides a valuable resource for human microbiome engineering.

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.

High Level of Interaction between Phages and Bacteria in an Artisanal Raw Milk Cheese Microbial Community

Microbial starter cultures are used in the production of many cheeses around the world, such as Parmigiano-Reggiano, in Italy, Époisses, in France, and Canastra, in Brazil, providing many of the unique features of these cheeses. Bacteriophages (phages) are ubiquitous and well known to modulate the structure of bacterial communities, and recent data indicate that cheeses contain a high abundance of naturally occurring phages. Here, we analyze the viral and bacterial metagenomes of Canastra cheese: a traditional artisanal Brazilian cheese produced using an endogenous starter culture and raw milk. Over 1,200 viral operational taxonomic units were recovered using both isolated viral-like particles and complete metagenomic DNA. Common viral families identified included Siphoviridae and Myoviridae, with 40% of putative phage genomes unidentified at the family level of classification. We observed very high phage diversity, which varied greatly across different cheese producers, with 28% of phage genomes detected in only one producer. Several metagenome-assembled genomes were recovered for lactic acid-producing bacteria, as well as nonstarter bacterial species, and we identified several phage-bacterium interactions, at the strain level of resolution, varying across distinct cheese producers. We postulate that at least one bacterial strain detected could be endogenous and unique to the Canastra cheese-producing region in Brazil and that its growth seems to be modulated by autochthonous phages present in this artisanal production system. This phage-host relationship is likely to influence the fermentation dynamics and ultimately the sensorial profile of these cheeses, with implications for other similar cheese production systems around the world. IMPORTANCE Our work demonstrated a dynamic yet stable microbial ecosystem during cheese production using an endogenous starter culture. This was observed across several distinct producers and was marked by genomic evidence of continued phage-bacterium interactions, such as the presence of bacterial defense mechanisms. Furthermore, we provide evidence of unique microbial signatures for each individual cheese producer studied in the region, a fact that may have profound consequences on product traceability. This was the first effort to describe and understand the bacteriophage composition and ecological dynamics within the Brazilian Canastra cheese production system. The study of this prototypical backslopping production system provides a solid background for further mechanistic studies of the production of many cheeses around the world.

Functional ecology of bacteriophages in the environment

Bacteriophages are as ubiquitous as their bacterial hosts and often more abundant. Understanding how bacteriophages control their bacterial host populations requires a number of different approaches. Bacteriophages can control bacterial populations through lysis, drive evolution of bacterial immunity systems through infection, provide a conduit for horizontal gene transfer and alter host metabolism by carriage of auxiliary metabolic genes. Understanding and quantifying how bacteriophages drive these processes, requires both technological developments to take measurements in situ, and laboratory-based studies to understand mechanisms. Technological advances have allowed quantification of the number of infected cells in situ, revealing far-lower levels than expected. Understanding how observations in laboratory conditions relate to what occurs in the environment, and experimental confirmation of the predicted function of phage genes from observations in environmental omics data, remains challenging.

Advances in the field of phage-based therapy with special emphasis on computational resources

In the current era, one of the major challenges is to manage the treatment of drug/antibiotic-resistant strains of bacteria. Phage therapy, a century-old technique, may serve as an alternative to antibiotics in treating bacterial infections caused by drug-resistant strains of bacteria. In this review, a systematic attempt has been made to summarize phage-based therapy in depth. This review has been divided into the following two sections: general information and computer-aided phage therapy (CAPT). In the case of general information, we cover the history of phage therapy, the mechanism of action, the status of phage-based products (approved and clinical trials) and the challenges. This review emphasizes CAPT, where we have covered primary phage-associated resources, phage prediction methods and pipelines. This review covers a wide range of databases and resources, including viral genomes and proteins, phage receptors, host genomes of phages, phage-host interactions and lytic proteins. In the post-genomic era, identifying the most suitable phage for lysing a drug-resistant strain of bacterium is crucial for developing alternate treatments for drug-resistant bacteria and this remains a challenging problem. Thus, we compile all phage-associated prediction methods that include the prediction of phages for a bacterial strain, the host for a phage and the identification of interacting phage-host pairs. Most of these methods have been developed using machine learning and deep learning techniques. This review also discussed recent advances in the field of CAPT, where we briefly describe computational tools available for predicting phage virions, the life cycle of phages and prophage identification. Finally, we describe phage-based therapy's advantages, challenges and opportunities.

Prophage-Derived Regions in Curtobacterium Genomes: Good Things, Small Packages

Curtobacterium is a genus of Gram-positive bacteria within the order Actinomycetales. Some Curtobacterium species (C. flaccumfaciens, C. plantarum) are harmful pathogens of agricultural crops such as soybean, dry beans, peas, sugar beet and beetroot, which occur throughout the world. Bacteriophages (bacterial viruses) are considered to be potential curative agents to control the spread of harmful bacteria. Temperate bacteriophages integrate their genomes into bacterial chromosomes (prophages), sometimes substantially influencing bacterial lifestyle and pathogenicity. About 200 publicly available genomes of Curtobacterium species, including environmental metagenomic sequences, were inspected for the presence of sequences of possible prophage origin using bioinformatic methods. The comparison of the search results with several ubiquitous bacterial groups showed the relatively low level of the presence of prophage traces in Curtobacterium genomes. Genomic and phylogenetic analyses were undertaken for the evaluation of the evolutionary and taxonomic positioning of predicted prophages. The analyses indicated the relatedness of Curtobacterium prophage-derived sequences with temperate actinophages of siphoviral morphology. In most cases, the predicted prophages can represent novel phage taxa not described previously. One of the predicted temperate phages was induced from the Curtobacterium genome. Bioinformatic analysis of the modelled proteins encoded in prophage-derived regions led to the discovery of some 100 putative glycopolymer-degrading enzymes that contained enzymatic domains with predicted cell-wall- and cell-envelope-degrading activity; these included glycosidases and peptidases. These proteins can be considered for the experimental design of new antibacterials against Curtobacterium phytopathogens.

Bioinformatics approaches for unveiling virus-host interactions

The coronavirus disease-2019 (COVID-19) pandemic has elucidated major limitations in the capacity of medical and research institutions to appropriately manage emerging infectious diseases. We can improve our understanding of infectious diseases by unveiling virus-host interactions through host range prediction and protein-protein interaction prediction. Although many algorithms have been developed to predict virus-host interactions, numerous issues remain to be solved, and the entire network remains veiled. In this review, we comprehensively surveyed algorithms used to predict virus-host interactions. We also discuss the current challenges, such as dataset biases toward highly pathogenic viruses, and the potential solutions. The complete prediction of virus-host interactions remains difficult; however, bioinformatics can contribute to progress in research on infectious diseases and human health.

Metagenomic analysis reveals unexplored diversity of archaeal virome in the human gut

The human gut microbiome has been extensively explored, while the archaeal viruses remain largely unknown. Here, we present a comprehensive analysis of the archaeal viruses from the human gut metagenomes and the existing virus collections using the CRISPR spacer and viral signature-based approach. This results in 1279 viral species, of which, 95.2% infect Methanobrevibacteria_A, 56.5% shared high identity (>95%) with the archaeal proviruses, 37.2% have a host range across archaeal species, and 55.7% are highly prevalent in the human population (>1%). A methanogenic archaeal virus-specific gene for pseudomurein endoisopeptidase (PeiW) frequently occurs in the viral sequences (n = 150). Analysis of 33 Caudoviricetes viruses with a complete genome often discovers the genes (integrase, n = 29; mazE, n = 10) regulating the viral lysogenic-lytic cycle, implying the dominance of temperate viruses in the archaeal virome. Together, our work uncovers the unexplored diversity of archaeal viruses, revealing the novel facet of the human gut microbiome.

Phage Encounters Recorded in CRISPR Arrays in the Genus Oenococcus

The Oenococcus genus comprises four recognized species, and members have been found in different types of beverages, including wine, kefir, cider and kombucha. In this work, we implemented two complementary strategies to assess whether oenococcal hosts of different species and habitats were connected through their bacteriophages. First, we investigated the diversity of CRISPR-Cas systems using a genome-mining approach, and CRISPR-endowed strains were identified in three species. A census of the spacers from the four identified CRISPR-Cas loci showed that each spacer space was mostly dominated by species-specific sequences. Yet, we characterized a limited records of potentially recent and also ancient infections between O. kitaharae and O. sicerae and phages of O. oeni, suggesting that some related phages have interacted in diverse ways with their Oenococcus hosts over evolutionary time. Second, phage-host interaction analyses were performed experimentally with a diversified panel of phages and strains. None of the tested phages could infect strains across the species barrier. Yet, some infections occurred between phages and hosts from distinct beverages in the O. oeni species.

Revealing bacteria-phage interactions in human microbiome through the CRISPR-Cas immune systems

The human gut microbiome is composed of a diverse consortium of microorganisms. Relatively little is known about the diversity of the bacteriophage population and their interactions with microbial organisms in the human microbiome. Due to the persistent rivalry between microbial organisms (hosts) and phages (invaders), genetic traces of phages are found in the hosts' CRISPR-Cas adaptive immune system. Mobile genetic elements (MGEs) found in bacteria include genetic material from phage and plasmids, often resultant from invasion events. We developed a computational pipeline (BacMGEnet), which can be used for inference and exploratory analysis of putative interactions between microbial organisms and MGEs (phages and plasmids) and their interaction network. Given a collection of genomes as the input, BacMGEnet utilizes computational tools we have previously developed to characterize CRISPR-Cas systems in the genomes, which are then used to identify putative invaders from publicly available collections of phage/prophage sequences. In addition, BacMGEnet uses a greedy algorithm to summarize identified putative interactions to produce a bacteria-MGE network in a standard network format. Inferred networks can be utilized to assist further examination of the putative interactions and for discovery of interaction patterns. Here we apply the BacMGEnet pipeline to a few collections of genomic/metagenomic datasets to demonstrate its utilities. BacMGEnet revealed a complex interaction network of the Phocaeicola vulgatus pangenome with its phage invaders, and the modularity analysis of the resulted network suggested differential activities of the different P. vulgatus' CRISPR-Cas systems (Type I-C and Type II-C) against some phages. Analysis of the phage-bacteria interaction network of human gut microbiome revealed a mixture of phages with a broad host range (resulting in large modules with many bacteria and phages), and phages with narrow host range. We also showed that BacMGEnet can be used to infer phages that invade bacteria and their interactions in wound microbiome. We anticipate that BacMGEnet will become an important tool for studying the interactions between bacteria and their invaders for microbiome research.

Insights into the global freshwater virome

Viruses are by far the most abundant life forms on this planet. Yet, the full viral diversity remains mostly unknown, especially in environments like freshwater. Therefore, we aimed to study freshwater viruses in a global context. To this end, we downloaded 380 publicly available viral metagenomes (>1 TB). More than 60% of these metagenomes were discarded based on their levels of cellular contamination assessed by ribosomal DNA content. For the remaining metagenomes, assembled contigs were decontaminated using two consecutive steps, eventually yielding 273,365 viral contigs longer than 1,000 bp. Long enough contigs (≥ 10 kb) were clustered to identify novel genomes/genome fragments. We could recover 549 complete circular and high-quality draft genomes, out of which 10 were recognized as being novel. Functional annotation of these genomes showed that most of the annotated coding sequences are DNA metabolic genes or phage structural genes. On the other hand, taxonomic analysis of viral contigs showed that most of the assigned contigs belonged to the order Caudovirales, particularly the families of Siphoviridae, Myoviridae, and Podoviridae. The recovered viral contigs contained several auxiliary metabolic genes belonging to several metabolic pathways, especially carbohydrate and amino acid metabolism in addition to photosynthesis as well as hydrocarbon degradation and antibiotic resistance. Overall, we present here a set of prudently chosen viral contigs, which should not only help better understanding of freshwater viruses but also be a valuable resource for future virome studies.

A relationship between phages and organic carbon in wastewater treatment plant effluents

With stringent effluent requirements and the implementation of new processes for micropollutant removal, it is increasingly important for wastewater treatment plants (WWTPs) to understand the factors affecting effluent quality. Phages (viruses infecting prokaryotes) are abundant in the biological treatment processes. They can contribute to organic carbon in the treated effluent both because they are organic in nature and occur in the effluent and because they cause lysis of microorganisms. Today very little is known about the effects of phages on effluent quality. The goal of this study was, therefore, to determine the relationship between phages and organic carbon in WWTP effluents. We also examined the diversity, taxonomy, and host-association of DNA phages using metagenomics. Effluent samples were collected from four WWTPs treating municipal wastewater. Significant differences in both organic carbon and virus-like particle concentrations were observed between the plants and there was a linear relationship between the two parameters. The phage communities were diverse with many members being taxonomically unclassified. Putative hosts were dominated by bacteria known to be abundant in activated sludge systems such as Comamonadaceae. The composition of phages differed between the WWTPs, suggesting that local conditions shape the communities. Overall, our findings suggest that the abundance and composition of phages are related to effluent quality. Thus, there is a need for further research clarifying the association between phage dynamics and WWTP function.

Prophage Diversity Across Salmonella and Verotoxin-Producing Escherichia coli in Agricultural Niches of British Columbia, Canada

Prophages have long been regarded as an important contributor to the evolution of Salmonella and Verotoxin-producing E. coli (VTEC), members of the Enterobacteriaceae that cause millions of cases of foodborne illness in North America. In S. Typhimurium, prophages provide many of the genes required for invasion; similarly, in VTEC, the Verotoxin-encoding genes are located in cryptic prophages. The ability of prophages to quickly acquire and lose genes have driven their rapid evolution, leading to highly diversified populations of phages that can infect distantly-related bacterial hosts. To defend against foreign genetic materials (i.e., phages), bacteria have evolved Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) immunity, consisting of variable spacer regions that match short nucleic acid sequences of invaders previously encountered. The number of spacer regions varies widely amongst Enterobacteriaceae, and there is currently no clear consensus if the accumulation of spacers is linked to genomic prophage abundance. Given the immense prophage diversity and contribution to bacterial host phenotypes, we analyzed the prophage sequences within 118 strains of Salmonella and VTEC, 117 of which are of agricultural origin. Overall, 130 unique prophage sequences were identified and they were found to be remarkably diverse with <50% nucleotide similarity, particularly with the Gifsy-1 group which was identified in several Salmonella serovars and interestingly, a strain of VTEC. Additionally, we identified a novel plasmid-like phage that carried antibiotic resistance and bacteriocin resistance genes. The strains analyzed carried at least six distinct spacers which did not possess homology to prophages identified in the same genome. In fact, only a fraction of all identified spacers (14%) possessed significant homology to known prophages. Regression models did not discern a correlation between spacer and prophage abundance in our strains, although the relatively high number of spacers in our strains (an average of 27 in Salmonella and 19 in VTEC) suggest that high rates of infection may occur in agricultural niches and be a contributing driver in bacterial evolution. Cumulatively, these results shed insight into prophage diversity of Salmonella and VTEC, which will have further implications when informing development of phage therapies against these foodborne pathogens.

Advances and challenges in cataloging the human gut virome

The human gut virome, which is often referred to as the "dark matter" of the gut microbiome, remains understudied. A better understanding of the composition and variations of the gut virome across populations is critical for exploring its impact on diseases and health. A series of advances in the characterization of human gut virome have unveiled high genetic diversity and various functional potentials of gut viruses. Here, we summarize the recently available human gut virome databases and discuss their features, procedures, and challenges with the intention to provide a reference to researchers to use while choosing a profiling database. We also propose a "best practice" for cataloging the viral population.

Transplantation of bacteriophages from ulcerative colitis patients shifts the gut bacteriome and exacerbates the severity of DSS colitis

BACKGROUND

Inflammatory bowel diseases (IBDs) including Crohn's disease (CD) and ulcerative colitis (UC) are characterized by chronic and debilitating gut inflammation. Altered bacterial communities of the intestine are strongly associated with IBD initiation and progression. The gut virome, which is primarily composed of bacterial viruses (bacteriophages, phages), is thought to be an important factor regulating and shaping microbial communities in the gut. While alterations in the gut virome have been observed in IBD patients, the contribution of these viruses to alterations in the bacterial community and heightened inflammatory responses associated with IBD patients remains largely unknown.

RESULTS

Here, we performed in vivo microbial cross-infection experiments to follow the effects of fecal virus-like particles (VLPs) isolated from UC patients and healthy controls on bacterial diversity and severity of experimental colitis in human microbiota-associated (HMA) mice. Shotgun metagenomics confirmed that several phages were transferred to HMA mice, resulting in treatment-specific alterations in the gut virome. VLPs from healthy and UC patients also shifted gut bacterial diversity of these mice, an effect that was amplified during experimental colitis. VLPs isolated from UC patients specifically altered the relative abundance of several bacterial taxa previously implicated in IBD progression. Additionally, UC VLP administration heightened colitis severity in HMA mice, as indicated by shortened colon length and increased pro-inflammatory cytokine production. Importantly, this effect was dependent on intact VLPs.

CONCLUSIONS

Our findings build on recent literature indicating that phages are dynamic regulators of bacterial communities in the gut and implicate the intestinal virome in modulating intestinal inflammation and disease. Video Abstract.

Three families of Asgard archaeal viruses identified in metagenome-assembled genomes

Asgardarchaeota harbour many eukaryotic signature proteins and are widely considered to represent the closest archaeal relatives of eukaryotes. Whether similarities between Asgard archaea and eukaryotes extend to their viromes remains unknown. Here we present 20 metagenome-assembled genomes of Asgardarchaeota from deep-sea sediments of the basin off the Shimokita Peninsula, Japan. By combining a CRISPR spacer search of metagenomic sequences with phylogenomic analysis, we identify three family-level groups of viruses associated with Asgard archaea. The first group, verdandiviruses, includes tailed viruses of the class Caudoviricetes (realm Duplodnaviria); the second, skuldviruses, consists of viruses with predicted icosahedral capsids of the realm Varidnaviria; and the third group, wyrdviruses, is related to spindle-shaped viruses previously identified in other archaea. More than 90% of the proteins encoded by these viruses of Asgard archaea show no sequence similarity to proteins encoded by other known viruses. Nevertheless, all three proposed families consist of viruses typical of prokaryotes, providing no indication of specific evolutionary relationships between viruses infecting Asgard archaea and eukaryotes. Verdandiviruses and skuldviruses are likely to be lytic, whereas wyrdviruses potentially establish chronic infection and are released without host cell lysis. All three groups of viruses are predicted to play important roles in controlling Asgard archaea populations in deep-sea ecosystems.

Salmonella phage akira, infecting selected Salmonella enterica Enteritidis and Typhimurium strains, represents a new lineage of bacteriophages

Some serovars of Salmonella can cause life-threatening diarrhoeal diseases and bacteriemia. The emergence of multidrug-resistant strains has led to a need for alternative treatments such as phage therapy, which requires available, well-described, diverse, and suitable phages. Phage akira was found to lyse 19 out of 32 Salmonella enterica serovars and farm isolates tested, although plaque formation was observed with only two S. Enteritidis and one S. Typhimurium strain. Phage akira encodes anti-defence genes against type 1 R-M systems, is distinct (<65% nucleotide sequence identity) from related phages and has siphovirus morphology. We propose that akira represents a new genus in the class Caudoviricetes.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

Short- and long-read metagenomics expand individualized structural variations in gut microbiomes

In-depth profiling of genetic variations in the gut microbiome is highly desired for understanding its functionality and impacts on host health and disease. Here, by harnessing the long read advantage provided by Oxford Nanopore Technology (ONT), we characterize fine-scale genetic variations of structural variations (SVs) in hundreds of gut microbiomes from healthy humans. ONT long reads dramatically improve the quality of metagenomic assemblies, enable reliable detection of a large, expanded set of structural variation types (notably including large insertions and inversions). We find SVs are highly distinct between individuals and stable within an individual, representing gut microbiome fingerprints that shape strain-level differentiations in function within species, complicating the associations to metabolites and host phenotypes such as blood glucose. In summary, our study strongly emphasizes that incorporating ONT reads into metagenomic analyses expands the detection scope of genetic variations, enables profiling strain-level variations in gut microbiome, and their intricate correlations with metabolome.

Here, Wang and colleagues combine short and long sequencing reads to characterize structural variations, prophage and CRISPR spacer elements in human gut microbiomes, and reveal functional differences at a finer level of bacterial strains.

PHISDetector: A Tool to Detect Diverse In Silico Phage-host Interaction Signals for Virome Studies

Phage-microbe interactions are appealing systems to study coevolution, and have also been increasingly emphasized due to their roles in human health, disease, and the development of novel therapeutics. Phage-microbe interactions leave diverse signals in bacterial and phage genomic sequences, defined as phage-host interaction signals (PHISs), which include clustered regularly interspaced short palindromic repeats (CRISPR) targeting, prophage, and protein-protein interaction signals. In the present study, we developed a novel tool phage-host interaction signal detector (PHISDetector) to predict phage-host interactions by detecting and integrating diverse in silico PHISs, and scoring the probability of phage-host interactions using machine learning models based on PHIS features. We evaluated the performance of PHISDetector on multiple benchmark datasets and application cases. When tested on a dataset of 758 annotated phage-host pairs, PHISDetector yields the prediction accuracies of 0.51 and 0.73 at the species and genus levels, respectively, outperforming other phage-host prediction tools. When applied to on 125,842 metagenomic viral contigs (mVCs) derived from 3042 geographically diverse samples, a detection rate of 54.54% could be achieved. Furthermore, PHISDetector could predict infecting phages for 85.6% of 368 multidrug-resistant (MDR) bacteria and 30% of 454 human gut bacteria obtained from the National Institutes of Health (NIH) Human Microbiome Project (HMP). The PHISDetector can be run either as a web server (http://www.microbiome-bigdata.com/PHISDetector/) for general users to study individual inputs or as a stand-alone version (https://github.com/HIT-ImmunologyLab/PHISDetector) to process massive phage contigs from virome studies.

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

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

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

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

Organizing the Global Diversity of Microviruses

Microviruses encompass an astonishing array of small, single-stranded DNA phages that, due to the surge in metagenomic surveys, are now known to be prevalent in most environments. Current taxonomy concedes the considerable diversity within this lineage to a single family (the Microviridae), which has rendered it difficult to adequately and accurately assess the amount of variation that actually exists within this group. We amassed and curated the largest collection of microviral genomes to date and, through a combination of protein-sharing networks and phylogenetic analysis, discovered at least three meaningful taxonomic levels between the current ranks of family and genus. When considering more than 13,000 microviral genomes from recognized lineages and as-yet-unclassified microviruses in metagenomic samples, microviral diversity is better understood by elevating microviruses to the level of an order that consists of three suborders and at least 19 putative families, each with their respective subfamilies. These revisions enable fine-scale assessment of microviral dynamics: for example, in the human gut, there are considerable differences in the abundances of microviral families both between urban and rural populations and in individuals over time. In addition, our analysis of genome contents and gene exchange shows that microviral families carry no recognizable accessory metabolic genes and rarely, if ever, engage in horizontal gene transfer across microviral families or with their bacterial hosts. These insights bring microviral taxonomy in line with current developments in the taxonomy of other phages and increase the understanding of microvirus biology.

Global Composition of the Bacteriophage Community in Honey Bees

The microbial communities in animal digestive systems are critical for host development and health. They stimulate the immune system during development, synthesize important chemical compounds like hormones, aid in digestion, competitively exclude pathogens, etc. Compared to the bacterial and fungal components of the microbiome, we know little about the temporal and spatial dynamics of bacteriophage communities in animal digestive systems. Recently, the bacteriophages of the honey bee gut were characterized in two European bee populations. Most of the bacteriophages described in these two reports were novel, harbored many metabolic genes in their genomes, and had a community structure that suggests coevolution with their bacterial hosts. To describe the conservation of bacteriophages in bees and begin to understand their role in the bee microbiome, we sequenced the virome of Apis mellifera from Austin, TX, and compared bacteriophage compositions among three locations around the world. We found that most bacteriophages from Austin are novel, sharing no sequence similarity with anything in public repositories. However, many bacteriophages are shared among the three bee viromes, indicating specialization of bacteriophages in the bee gut. Our study, along with the two previous bee virome studies, shows that the bee gut bacteriophage community is simple compared to that of many animals, consisting of several hundred types of bacteriophages that primarily infect four of the dominant bacterial phylotypes in the bee gut. IMPORTANCE Viruses that infect bacteria (bacteriophages) are abundant in the microbial communities that live on and in plants and animals. However, our knowledge of the structure, dynamics, and function of these viral communities lags far behind our knowledge of their bacterial hosts. We sequenced the first bacteriophage community of honey bees from the United States and compared the U.S. honey bee bacteriophage community to those of samples from Europe. Our work is an important characterization of an economically critical insect species and shows how bacteriophage communities can contain highly conserved individuals and be highly variable in composition across a wide geographic range.