Taxonomic assignment of uncultivated prokaryotic virus genomes is enabled by gene-sharing networks

Viromes vs. mixed community metagenomes: choice of method dictates interpretation of viral community ecology

Background

Viruses, the majority of which are uncultivated, are among the most abundant biological entities on Earth. From altering microbial physiology to driving community dynamics, viruses are fundamental members of microbiomes. While the number of studies leveraging viral metagenomics (viromics) for studying uncultivated viruses is growing, standards for viromics research are lacking. Viromics can utilize computational discovery of viruses from total metagenomes of all community members (hereafter metagenomes) or use physical separation of virus-specific fractions (hereafter viromes). However, differences in the recovery and interpretation of viruses from metagenomes and viromes obtained from the same samples remain understudied.

Results

Here, we compare viral communities from paired viromes and metagenomes obtained from 60 diverse samples across human gut, soil, freshwater, and marine ecosystems. Overall, viral communities obtained from viromes were more abundant and species rich than those obtained from metagenomes, although there were some exceptions. Despite this, metagenomes still contained many viral genomes not detected in viromes. We also found notable differences in the predicted lytic state of viruses detected in viromes vs metagenomes at the time of sequencing. Other forms of variation observed include genome presence/absence, genome quality, and encoded protein content between viromes and metagenomes, but the magnitude of these differences varied by environment.

Conclusions

Overall, our results show that the choice of method can lead to differing interpretations of viral community ecology. We suggest that the choice of whether to target a metagenome or virome to study viral communities should be dependent on the environmental context and ecological questions being asked. However, our overall recommendation to researchers investigating viral ecology and evolution is to pair both approaches to maximize their respective benefits.

Investigating the Human Intestinal DNA Virome and Predicting Disease-Associated Virus-Host Interactions in Severe Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS)

Understanding how the human virome, and which of its constituents, contributes to health or disease states is reliant on obtaining comprehensive virome profiles. By combining DNA viromes from isolated virus-like particles (VLPs) and whole metagenomes from the same faecal sample of a small cohort of healthy individuals and patients with severe myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS), we have obtained a more inclusive profile of the human intestinal DNA virome. Key features are the identification of a core virome comprising tailed phages of the class Caudoviricetes, and a greater diversity of DNA viruses including extracellular phages and integrated prophages. Using an in silico approach, we predicted interactions between members of the Anaerotruncus genus and unique viruses present in ME/CFS microbiomes. This study therefore provides a framework and rationale for studies of larger cohorts of patients to further investigate disease-associated interactions between the intestinal virome and the bacteriome.

Gut virome profiling identifies an association between temperate phages and colorectal cancer promoted by Helicobacter pylori infection

Colorectal cancer (CRC) is one of the most commonly diagnosed cancers worldwide. While a close correlation between chronic Helicobacter pylori infection and CRC has been reported, the role of the virome has been overlooked. Here, we infected Apc-mutant mouse models and C57BL/6 mice with H. pylori and conducted a comprehensive metagenomics analysis of H. pylori-induced changes in lower gastrointestinal tract bacterial and viral communities. We observed an expansion of temperate phages in H. pylori infected Apc+/1638N mice at the early stage of carcinogenesis. Some of the temperate phages were predicted to infect bacteria associated with CRC, including Enterococcus faecalis. We also observed a high prevalence of virulent genes, such as flgJ, cwlJ, and sleB, encoded by temperate phages. In addition, we identified phages associated with pre-onset and onset of H. pylori-promoted carcinogenesis. Through co-occurrence network analysis, we found strong associations between the viral and bacterial communities in infected mice before the onset of carcinogenesis. These findings suggest that the expansion of temperate phages, possibly caused by prophage induction triggered by H. pylori infection, may have contributed to the development of CRC in mice by interacting with the bacterial community.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

Characterizations of the multi-kingdom gut microbiota in Chinese patients with gouty arthritis

OBJECTIVE

The gut microbial composition has been linked to metabolic and autoimmune diseases, including arthritis. However, there is a dearth of knowledge on the gut bacteriome, mycobiome, and virome in patients with gouty arthritis (GA).

METHODS

We conducted a comprehensive analysis of the multi-kingdom gut microbiome of 26 GA patients and 28 healthy controls, using whole-metagenome shotgun sequencing of their stool samples.

RESULTS

Profound alterations were observed in the gut bacteriome, mycobiome, and virome of GA patients. We identified 1,117 differentially abundant bacterial species, 23 fungal species, and 4,115 viral operational taxonomic units (vOTUs). GA-enriched bacteria included Escherichia coli_D GENOME144544, Bifidobacterium infantis GENOME095938, Blautia_A wexlerae GENOME096067, and Klebsiella pneumoniae GENOME147598, while control-enriched bacteria comprised Faecalibacterium prausnitzii_G GENOME147678, Agathobacter rectalis GENOME143712, and Bacteroides_A plebeius_A GENOME239725. GA-enriched fungi included opportunistic pathogens like Cryptococcus neoformans GCA_011057565, Candida parapsilosis GCA_000182765, and Malassezia spp., while control-enriched fungi featured several Hortaea werneckii subclades and Aspergillus fumigatus GCA_000002655. GA-enriched vOTUs mainly attributed to Siphoviridae, Myoviridae, Podoviridae, and Microviridae, whereas control-enriched vOTUs spanned 13 families, including Siphoviridae, Myoviridae, Podoviridae, Quimbyviridae, Phycodnaviridae, and crAss-like. A co-abundance network revealed intricate interactions among these multi-kingdom signatures, signifying their collective influence on the disease. Furthermore, these microbial signatures demonstrated the potential to effectively discriminate between patients and controls, highlighting their diagnostic utility.

CONCLUSIONS

This study yields crucial insights into the characteristics of the GA microbiota that may inform future mechanistic and therapeutic investigations.

Characterization of a Vibriophage Infecting Pathogenic Vibrio harveyi

Bacterial diseases caused by Vibrio spp. are prevalent in aquaculture and can lead to high mortality rates among aquatic species and significant economic losses. With the increasing emergence of multidrug-resistant Vibrio strains, phage therapy is being explored as a potential alternative to antibiotics for biocontrol of infectious diseases. Here, a new lytic phage named vB_VhaS_R21Y (R21Y) was isolated against Vibrio harveyi BVH1 obtained from seawater from a scallop-farming area in Rongcheng, China. Its morphology, infection cycle, lytic profile, phage stability, and genetic features were characterized. Transmission electronic microscopy indicated that R21Y is siphovirus-like, comprising an icosahedral head (diameter 73.31 ± 2.09 nm) and long noncontractile tail (205.55 ± 0.75 nm). In a one-step growth experiment, R21Y had a 40-min latent period and a burst size of 35 phage particles per infected cell. R21Y was highly species-specific in the host range test and was relatively stable at pH 4-10 and 4-55 °C. Genomic analysis showed that R21Y is a double-stranded DNA virus with a genome size of 82,795 bp and GC content of 47.48%. Its high tolerance and lytic activity indicated that R21Y may be a candidate for phage therapy in controlling vibriosis in aquacultural systems.

Characterization of phage vB_EcoS-EE09 infecting E. coli DSM613 Isolated from Wastewater Treatment Plant Effluent and Comparative Proteomics of the Infected and Non-Infected Host

Phages influence microbial communities, can be applied in phage therapy, or may serve as bioindicators, e.g., in (waste)water management. We here characterized the Escherichia phage vB_EcoS-EE09 isolated from an urban wastewater treatment plant effluent. Phage vB_EcoS-EE09 belongs to the genus Dhillonvirus, class Caudoviricetes. It has an icosahedral capsid with a long non-contractile tail and a dsDNA genome with an approximate size of 44 kb and a 54.6% GC content. Phage vB_EcoS-EE09 infected 12 out of the 17 E. coli strains tested. We identified 16 structural phage proteins, including the major capsid protein, in cell-free lysates by protein mass spectrometry. Comparative proteomics of protein extracts of infected E. coli cells revealed that proteins involved in amino acid and protein metabolism were more abundant in infected compared to non-infected cells. Among the proteins involved in the stress response, 74% were less abundant in the infected cultures compared to the non-infected controls, with six proteins showing significant less abundance. Repressing the expression of these proteins may be a phage strategy to evade host defense mechanisms. Our results contribute to diversifying phage collections, identifying structural proteins to enable better reliability in annotating taxonomically related phage genomes, and understanding phage-host interactions at the protein level.

Biochemical characterization of mRNA capping enzyme from Faustovirus

The mammalian mRNA 5' cap structures play important roles in cellular processes such as nuclear export, efficient translation, and evading cellular innate immune surveillance and regulating 5'-mediated mRNA turnover. Hence, installation of the proper 5' cap is crucial in therapeutic applications of synthetic mRNA. The core 5' cap structure, Cap-0, is generated by three sequential enzymatic activities: RNA 5' triphosphatase, RNA guanylyltransferase, and cap N7-guanine methyltransferase. Vaccinia virus RNA capping enzyme (VCE) is a heterodimeric enzyme that has been widely used in synthetic mRNA research and manufacturing. The large subunit of VCE D1R exhibits a modular structure where each of the three structural domains possesses one of the three enzyme activities, whereas the small subunit D12L is required to activate the N7-guanine methyltransferase activity. Here, we report the characterization of a single-subunit RNA capping enzyme from an amoeba giant virus. Faustovirus RNA capping enzyme (FCE) exhibits a modular array of catalytic domains in common with VCE and is highly efficient in generating the Cap-0 structure without an activation subunit. Phylogenetic analysis suggests that FCE and VCE are descended from a common ancestral capping enzyme. We found that compared to VCE, FCE exhibits higher specific activity, higher activity toward RNA containing secondary structures and a free 5' end, and a broader temperature range, properties favorable for synthetic mRNA manufacturing workflows.

Top-heavy trophic structure within benthic viral dark matter

A better understanding of system-specific viral ecology in diverse environments is needed to predict patterns of virus-host trophic structure in the Anthropocene. This study characterised viral-host trophic structure within coral reef benthic cyanobacterial mats-a globally proliferating cause and consequence of coral reef degradation. We employed deep longitudinal multi-omic sequencing to characterise the viral assemblage (ssDNA, dsDNA, and dsRNA viruses) and profile lineage-specific host-virus interactions within benthic cyanobacterial mats sampled from Bonaire, Caribbean Netherlands. We recovered 11,012 unique viral populations spanning at least 10 viral families across the orders Caudovirales, Petitvirales, and Mindivirales. Gene-sharing network analyses provided evidence for extensive genomic novelty of mat viruses from reference and environmental viral sequences. Analysis of coverage ratios of viral sequences and computationally predicted hosts spanning 15 phyla and 21 classes revealed virus-host abundance (from DNA) and activity (from RNA) ratios consistently exceeding 1:1, suggesting a top-heavy intra-mat trophic structure with respect to virus-host interactions. Overall, our article contributes a curated database of viral sequences found in Caribbean coral reef benthic cyanobacterial mats (vMAT database) and provides multiple lines of field-based evidence demonstrating that viruses are active members of mat communities, with broader implications for mat functional ecology and demography.

A compendium of viruses from methanogenic archaea reveals their diversity and adaptations to the gut environment

Methanogenic archaea are major producers of methane, a potent greenhouse gas and biofuel, and are widespread in diverse environments, including the animal gut. The ecophysiology of methanogens is likely impacted by viruses, which remain, however, largely uncharacterized. Here we carried out a global investigation of viruses associated with all current diversity of methanogens by assembling an extensive CRISPR database consisting of 156,000 spacers. We report 282 high-quality (pro)viral and 205 virus-like/plasmid sequences assigned to hosts belonging to ten main orders of methanogenic archaea. Viruses of methanogens can be classified into 87 families, underscoring a still largely undiscovered genetic diversity. Viruses infecting gut-associated archaea provide evidence of convergence in adaptation with viruses infecting gut-associated bacteria. These viruses contain a large repertoire of lysin proteins that cleave archaeal pseudomurein and are enriched in glycan-binding domains (Ig-like/Flg_new) and diversity-generating retroelements. The characterization of this vast repertoire of viruses paves the way towards a better understanding of their role in regulating methanogen communities globally, as well as the development of much-needed genetic tools.

Viral but not bacterial community successional patterns reflect extreme turnover shortly after rewetting dry soils

As central members of soil trophic networks, viruses have the potential to drive substantial microbial mortality and nutrient turnover. Pinpointing viral contributions to terrestrial ecosystem processes remains a challenge, as temporal dynamics are difficult to unravel in the spatially and physicochemically heterogeneous soil environment. In Mediterranean grasslands, the first rainfall after seasonal drought provides an ecosystem reset, triggering microbial activity during a tractable window for capturing short-term dynamics. Here, we simulated precipitation in microcosms from four distinct dry grassland soils and generated 144 viromes, 84 metagenomes and 84 16S ribosomal RNA gene amplicon datasets to characterize viral, prokaryotic and relic DNA dynamics over 10 days. Vastly different viral communities in each soil followed remarkably similar successional trajectories. Wet-up triggered a significant increase in viral richness, followed by extensive compositional turnover. Temporal succession in prokaryotic communities was much less pronounced, perhaps suggesting differences in the scales of activity captured by viromes (representing recently produced, ephemeral viral particles) and total DNA. Still, differences in the relative abundances of Actinobacteria (enriched in dry soils) and Proteobacteria (enriched in wetted soils) matched those of their predicted phages, indicating viral predation of dominant bacterial taxa. Rewetting also rapidly depleted relic DNA, which subsequently reaccumulated, indicating substantial new microbial mortality in the days after wet-up, particularly of the taxa putatively under phage predation. Production of abundant, diverse viral particles via microbial host cell lysis appears to be a conserved feature of the early response to soil rewetting, and results suggest the potential for 'Cull-the-Winner' dynamics, whereby viruses infect and cull but do not decimate dominant host populations.

Diversity and function of mountain and polar supraglacial DNA viruses

Mountain and polar glaciers cover 10% of the Earth's surface and are typically extreme environments that challenge life of all forms. Viruses are abundant and active in supraglacial ecosystems and play a crucial role in controlling the supraglacial microbial communities. However, our understanding of virus ecology on glacier surfaces and their potential impacts on downstream ecosystems remains limited. Here, we present the supraglacial virus genome (SgVG) catalog, a 15-fold expanded genomic inventory of 10,840 DNA-virus species from 38 mountain and polar glaciers, spanning habitats such as snow, ice, meltwater, and cryoconite. Supraglacial DNA-viruses were highly specific compared to viruses in other ecosystems yet exhibited low public health risks. Supraglacial viral communities were primarily constrained by habitat, with cryoconite displaying the highest viral activity levels. We observed a prevalence of lytic viruses in all habitats, especially in cryoconite, but a high level of lysogenic viruses in snow and ice. Additionally, we found that supraglacial viruses could be linked to ∼83% of obtained prokaryotic phyla/classes and possessed the genetic potential to promote metabolism and increase cold adaptation, cell mobility, and phenolic carbon use of hosts in hostile environmental conditions using diverse auxiliary metabolic genes. Our results provide the first systematic characterization of the diversity, function, and public health risks evaluation of mountain and polar supraglacial DNA viruses. This understanding of glacial viruses is crucial for function assessments and ecological modeling of glacier ecosystems, especially for the Tibetan Plateau's Mountain glaciers, which support ∼20% of the human populations on Earth.

Virus diversity and activity is driven by snowmelt and host dynamics in a high-altitude watershed soil ecosystem

BACKGROUND

Viruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics' seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes.

RESULTS

We identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages (Leviviricetes), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a "seed-bank" viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: "Kill-The-Winner" which proposes that lytic phages are actively infecting abundant bacteria, and "Piggyback-The-Persistent" which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth.

CONCLUSION

Taken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry. Video Abstract.

Identification of hidden N4-like viruses and their interactions with hosts

IMPORTANCE

The findings of this study are significant, as N4-like viruses represent a unique viral lineage with a distinct replication mechanism and a conserved core genome. This work has resulted in a comprehensive global map of the entire N4-like viral lineage, including information on their distribution in different biomes, evolutionary divergence, genomic diversity, and the potential for viral-mediated host metabolic reprogramming. As such, this work significantly contributes to our understanding of the ecological function and viral-host interactions of bacteriophages.

Global diversity and biogeography of DNA viral communities in activated sludge systems

BACKGROUND

Activated sludge (AS) systems in wastewater treatment plants (WWTPs) harbor enormous viruses that regulate microbial metabolism and nutrient cycling, significantly influencing the stability of AS systems. However, our knowledge about the diversity of viral taxonomic groups and functional traits in global AS systems is still limited. To address this gap, we investigated the global diversity and biogeography of DNA viral communities in AS systems using 85,114 viral operational taxonomic units (vOTUs) recovered from 144 AS samples collected across 54 WWTPs from 13 different countries.

RESULTS

AS viral communities and their functional traits exhibited distance-decay relationship (DDR) at the global scale and latitudinal diversity gradient (LDG) from equator to mid-latitude. Furthermore, it was observed that AS viral community and functional gene structures were largely driven by the geographic factors and wastewater types, of which the geographic factors were more important. Carrying and disseminating auxiliary metabolic genes (AMGs) associated with the degradation of polysaccharides, sulfate reduction, denitrification, and organic phosphoester hydrolysis, as well as the lysis of crucial functional microbes that govern biogeochemical cycles were two major ways by which viruses could regulate AS functions. It was worth noting that our study revealed a high abundance of antibiotic resistance genes (ARGs) in viral genomes, suggesting that viruses were key reservoirs of ARGs in AS systems.

CONCLUSIONS

Our results demonstrated the highly diverse taxonomic groups and functional traits of viruses in AS systems. Viral lysis of host microbes and virus-mediated HGT can regulate the biogeochemical and nutrient cycles, thus affecting the performance of AS systems. These findings provide important insights into the viral diversity, function, and ecology in AS systems on a global scale. Video Abstract.

Isolation and molecular characterization of the Salmonella Typhimurium orphan phage Arash

The current threat of multidrug resistant strains necessitates development of alternatives to antibiotics such as bacteriophages. This study describes the isolation and characterization of a novel Salmonella Typhimurium phage 'Arash' from hospital wastewater in Leuven, Belgium. Arash has a myovirus morphology with a 95 nm capsid and a 140 nm tail. The host range of Arash is restricted to its isolation host. Approximately 86% of the phage particles are adsorbed to a host cell within 10 min. Arash has latent period of 65 min and burst size of 425 PFU/cell. Arash has a dsDNA genome of 180,819 bp with GC content of 53.02% with no similarities to any characterized phages, suggesting Arash as a novel species in the novel 'Arashvirus' genus. Arash carries no apparent lysogeny-, antibiotic resistance- nor virulence-related genes. Proteome analysis revealed 116 proteins as part of the mature phage particles of which 27 could be assigned a function. Therefore, the present findings shed light on the morphological, microbiological and genomic characteristics of Arash and suggest its potential application as therapeutic and/or biocontrol agent.

Metagenomic analysis unveils the underexplored roles of prokaryotic viruses in a full-scale landfill leachate treatment plant

Enormous viral populations have been identified in activated sludge systems, but their ecological and biochemical roles in landfill leachate treatment plants remain poorly understood. To address this knowledge gap, we conducted an in-depth analysis using 36 metagenomic datasets that we collected and sequenced during a half-year time-series sampling campaign at six sites in a full-scale landfill leachate treatment plant (LLTP), elucidating viral distribution, virus‒host dynamics, virus-encoded auxiliary metabolic genes (AMGs), and viral contributions to the spread of virulence and antibiotic resistance genes. Our findings demonstrated that viral and prokaryotic communities differed widely among different treatment units, with stability over time. LLTP viruses were linked to various prokaryotic hosts, spanning 35 bacterial phyla and one archaeal phylum, which included the core microbes involved in biological treatments, as well as some of the less well-characterized microbial dark matter phyla. By encoding 2364 auxiliary metabolic genes (AMGs), viruses harbored the potential to regulate microbial nucleotide metabolism, facilitate the biodegradation of complex organic matter, and enhance flocculation and settling in biological treatment plants. The abundance distribution of AMGs varied considerably across treatment units and showed a lifestyle-dependent pattern with temperate virus-associated AMGs exhibiting a higher average abundance in downstream biological treatment units and effluent water. Meanwhile, temperate viruses tended to carry a higher load of virulence factor genes (VFGs), antibiotic resistance genes (ARGs), and biotic and metal resistance genes (BMRGs), and engaged in more frequent gene exchanges with prokaryotes than lytic viruses, thus acting as a pivotal contributor to the dissemination of pathogenicity and resistance genes in downstream LLTP units. This study provided a comprehensive profile of viral and prokaryotic communities in the LLTP and unveiled the varying roles of different-lifestyle viruses in biochemical processes and water quality safety.

Characterization and genomic analysis of Stutzerimonas stutzeri phage vB_PstS_ZQG1, representing a novel viral genus

Stutzerimonas stutzeri is an opportunistic pathogenic bacterium belonging to the Gammaproteobacteria, exhibiting wide distribution in the environment and playing significant ecological roles such as nitrogen fixation or pollutant degradation. Despite its ecological importance, only two S. stutzeri phages have been isolated to date. Here, a novel S. stutzeri phage, vB_PstS_ZQG1, was isolated from the surface seawater of Qingdao, China. Transmission electron microscopy analysis indicates that vB_PstS_ZQG1 has a morphology characterized by a long non-contractile tail. The genomic sequence of vB_PstS_ZQG1 contains a linear, double-strand 61,790-bp with the G+C content of 53.24% and encodes 90 putative open reading frames. Two auxiliary metabolic genes encoding TolA protein and nucleotide pyrophosphohydrolase were identified, which are likely involved in host adaptation and phage reproduction. Phylogenetic and comparative genomic analyses demonstrated that vB_PstS_ZQG1 exhibits low similarity with previously isolated phages or uncultured viruses (average nucleotide identity values range from 21.7 to 29.4), suggesting that it represents a novel viral genus by itself, here named as Fuevirus. Biogeographic analysis showed that vB_PstS_ZQG1 was only detected in epipelagic and mesopelagic zone with low abundance. In summary, our findings of the phage vB_PstS_ZQG1 will provide helpful insights for further research on the interactions between S. stutzeri phages and their hosts, and contribute to discovering unknown viral sequences in the metagenomic database.

Distinct adaptive strategies and microbial interactions of soil viruses under different metal(loid) contaminations

Viruses, as the most abundant organisms, significantly influence ecological function and microbial survival in soils, yet little was known about how viruses and virus-microbe interactions respond to environmental stresses induced by metal(loid) contaminations. Here, we conducted the metagenomic analysis to investigate the adaptative mechanisms of soil viruses under different metal(loid) contamination levels. By capturing a catalogue of 23,066 viruses, we found that viral communities exhibited the increased richness, diversity, and the temperate to lytic ratio in facing the highest metal(loid) contaminations. Meanwhile, viruses displayed obvious lineage-specific infection modes to distinct dominant hosts under different pollution levels. Viral functions linking to the inhibition of transcription and the enhancement of DNA repairment as well as multiple resistance not only contributed to coping with elevated multiple metal(loid) stresses, but also facilitated the adaptation and functioning of viral hosts. Moreover, the harmonious coexistence of viruses and resistant/pathogenic bacteria under the heaviest contaminations potentially exacerbated disseminating resistance and pathogenicity, while viruses under the lightest contaminations might be natural predators of resistant/pathogenic bacteria through lysing host cells. Overall, this study highlights the ecological importance of viral adaptation and the interactions between viruses and resistant/pathogenic bacteria in contaminated environments, contributing to developing virus-based approaches to soil restoration.

Marine viruses disperse bidirectionally along the natural water cycle

Marine viruses in seawater have frequently been studied, yet their dispersal from neuston ecosystems at the air-sea interface towards the atmosphere remains a knowledge gap. Here, we show that 6.2% of the studied virus population were shared between air-sea interface ecosystems and rainwater. Virus enrichment in the 1-mm thin surface microlayer and sea foams happened selectively, and variant analysis proved virus transfer to aerosols collected at ~2 m height above sea level and rain. Viruses detected in rain and these aerosols showed a significantly higher percent G/C base content compared to marine viruses. CRISPR spacer matches of marine prokaryotes to foreign viruses from rainwater prove regular virus-host encounters at the air-sea interface. Our findings on aerosolization, adaptations, and dispersal support transmission of viruses along the natural water cycle.

Unveiling the hidden role of aquatic viruses in hydrocarbon pollution bioremediation

Hydrocarbon pollution poses substantial environmental risks to water and soil. Bioremediation, which utilizes microorganisms to manage pollutants, offers a cost-effective solution. However, the role of viruses, particularly bacteriophages (phages), in bioremediation remains unexplored. This study examines the diversity and activity of hydrocarbon-degradation genes encoded by environmental viruses, focusing on phages, within public databases. We identified 57 high-quality phage-encoded auxiliary metabolic genes (AMGs) related to hydrocarbon degradation, which we refer to as virus-encoded hydrocarbon degradation genes (vHYDEGs). These genes are encoded by taxonomically diverse aquatic phages and highlight the under-characterized global virosphere. Six protein families involved in the initial alkane hydroxylation steps were identified. Phylogenetic analyses revealed the diverse evolutionary trajectories of vHYDEGs across habitats, revealing previously unknown biodegraders linked evolutionarily with vHYDEGs. Our findings suggest phage AMGs may contribute to alkane and aromatic hydrocarbon degradation, participating in the initial, rate-limiting hydroxylation steps, thereby aiding hydrocarbon pollution bioremediation and promoting their propagation. To support future research, we developed vHyDeg, a database containing identified vHYDEGs with comprehensive annotations, facilitating the screening of hydrocarbon degradation AMGs and encouraging their bioremediation applications.

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.

Metagenomic stable isotope probing reveals bacteriophage participation in soil carbon cycling

Soil viruses are important components of the carbon (C) cycle, yet we still know little about viral ecology in soils. We added diverse 13 C-labelled carbon sources to soil and we used metagenomic-SIP to detect 13 C assimilation by viruses and their putative bacterial hosts. These data allowed us to link a 13 C-labelled bacteriophage to its 13 C-labelled Streptomyces putative host, and we used qPCR to track the dynamics of the putative host and phage in response to C inputs. Following C addition, putative host numbers increased rapidly for 3 days, and then more gradually, reaching maximal abundance on Day 6. Viral abundance and virus:host ratio increased dramatically over 6 days, and remained high thereafter (8.42 ± 2.94). From Days 6 to 30, virus:host ratio remained high, while putative host numbers declined more than 50%. Putative host populations were 13 C-labelled on Days 3-30, while 13 C-labelling of phage was detected on Days 14 and 30. This dynamic suggests rapid growth and 13 C-labelling of the host fueled by new C inputs, followed by extensive host mortality driven by phage lysis. These findings indicate that the viral shunt promotes microbial turnover in soil following new C inputs, thereby altering microbial community dynamics, and facilitating soil organic matter production.

Revisiting evolutionary trajectories and the organization of the Pleolipoviridae family

Archaeal pleomorphic viruses belonging to the Pleolipoviridae family represent an enigmatic group as they exhibit unique genomic features and are thought to have evolved through recombination with different archaeal plasmids. However, most of our understanding of the diversity and evolutionary trajectories of this clade comes from a handful of isolated representatives. Here we present 164 new genomes of pleolipoviruses obtained from metagenomic data of Australian hypersaline lakes and publicly available metagenomic data. We perform a comprehensive analysis on the diversity and evolutionary relationships of the newly discovered viruses and previously described pleolipoviruses. We propose to classify the viruses into five genera within the Pleolipoviridae family, with one new genus represented only by virus genomes retrieved in this study. Our data support the current hypothesis that pleolipoviruses reshaped their genomes through recombining with multiple different groups of plasmids, which is reflected in the diversity of their predicted replication strategies. We show that the proposed genus Epsilonpleolipovirus has evolutionary ties to pRN1-like plasmids from Sulfolobus, suggesting that this group could be infecting other archaeal phyla. Interestingly, we observed that the genome size of pleolipoviruses is correlated to the presence or absence of an integrase. Analyses of the host range revealed that all but one virus exhibit an extremely narrow range, and we show that the predicted tertiary structure of the spike protein is strongly associated with the host family, suggesting a specific adaptation to the host S-layer glycoprotein organization.

Phage-bacteria dynamics during the first years of life revealed by trans-kingdom marker gene analysis

Humans are colonized with commensal bacteria soon after birth, and, while this colonization is affected by lifestyle and other factors, bacterial colonization proceeds through well-studied phases. However, less is known about phage communities in early human development due to small study sizes, inability to leverage large databases, and lack of appropriate bioinformatics tools. In this study, whole genome shotgun sequencing data from the TEDDY study, composed of 12,262 longitudinal samples from 887 children in 4 countries, is reanalyzed to assess phage and bacterial dynamics simultaneously. Reads from these samples were mapped to marker genes from both bacteria and a new database of tens of thousands of phage taxa from human microbiomes. We uncover that each child is colonized by hundreds of different phages during the early years, and phages are more transitory than bacteria. Participants' samples continually harbor new phage species over time whereas the diversification of bacterial species begins to saturate. Phage data improves the ability for machine learning models to discriminate samples by country. Finally, while phage populations were individual-specific, striking patterns arose from the larger dataset, showing clear trends of ecological succession amongst phages, which correlated well with putative host bacteria. Improved understanding of phage-bacterial relationships may reveal new means by which to shape and modulate the microbiome and its constituents to improve health and reduce disease, particularly in vulnerable populations where antibiotic use and/or other more drastic measures may not be advised.

A systematic analysis of marine lysogens and proviruses

Viruses are ubiquitous in the oceans, exhibiting high abundance and diversity. Here, we systematically analyze existing genomic sequences of marine prokaryotes to compile a Marine Prokaryotic Genome Dataset (MPGD, consisting of over 12,000 bacterial and archaeal genomes) and a Marine Temperate Viral Genome Dataset (MTVGD). At least 40% of the MPGD genomes contain one or more proviral sequences, indicating that they are lysogens. The MTVGD includes over 12,900 viral contigs or putative proviruses, clustered into 10,897 viral genera. We show that lysogens and proviruses are abundant in marine ecosystems, particularly in the deep sea, and marine lysogens differ from non-lysogens in multiple genomic features and growth properties. We reveal several virus-host interaction networks of potential ecological relevance, and identify proviruses that appear to be able to infect (or to be transferred between) different bacterial classes and phyla. Auxiliary metabolic genes in the MTVGD are enriched in functions related to carbohydrate metabolism. Finally, we experimentally demonstrate the impact of a prophage on the transcriptome of a representative marine Shewanella bacterium. Our work contributes to a better understanding of the ecology of marine prokaryotes and their viruses.

PhaGenus: genus-level classification of bacteriophages using a Transformer model

MOTIVATION

Bacteriophages (phages for short), which prey on and replicate within bacterial cells, have a significant role in modulating microbial communities and hold potential applications in treating antibiotic resistance. The advancement of high-throughput sequencing technology contributes to the discovery of phages tremendously. However, the taxonomic classification of assembled phage contigs still faces several challenges, including high genetic diversity, lack of a stable taxonomy system and limited knowledge of phage annotations. Despite extensive efforts, existing tools have not yet achieved an optimal balance between prediction rate and accuracy.

RESULTS

In this work, we develop a learning-based model named PhaGenus, which conducts genus-level taxonomic classification for phage contigs. PhaGenus utilizes a powerful Transformer model to learn the association between protein clusters and support the classification of up to 508 genera. We tested PhaGenus on four datasets in different scenarios. The experimental results show that PhaGenus outperforms state-of-the-art methods in predicting low-similarity datasets, achieving an improvement of at least 13.7%. Additionally, PhaGenus is highly effective at identifying previously uncharacterized genera that are not represented in reference databases, with an improvement of 8.52%. The analysis of the infants' gut and GOV2.0 dataset demonstrates that PhaGenus can be used to classify more contigs with higher accuracy.

Characterization of 29 newly isolated bacteriophages as a potential therapeutic agent against IMP-6-producing Klebsiella pneumoniae from clinical specimens

Carbapenemase-producing Enterobacteriaceae (CPE) are one of the most detrimental species of antibiotic-resistant bacteria globally. Phage therapy has emerged as an effective strategy for the treatment of CPE infections. In western Japan, the rise of Klebsiella pneumoniae strains harboring the pKPI-6 plasmid encoding bla IMP-6 is of increasing concern. To address this challenge, we isolated 29 phages from Japanese sewage, specifically targeting 31 K. pneumoniae strains and one Escherichia coli strain harboring the pKPI-6 plasmid. Electron microscopy analysis revealed that among the 29 isolated phages, 21 (72.4%), 5 (17.2%), and 3 (10.3%) phages belonged to myovirus, siphovirus, and podovirus morphotypes, respectively. Host range analysis showed that 18 Slopekvirus strains within the isolated phages infected 25-26 K. pneumoniae strains, indicating that most of the isolated phages have a broad host range. Notably, K. pneumoniae strain Kp21 was exclusively susceptible to phage øKp_21, whereas Kp22 exhibited susceptibility to over 20 phages. Upon administering a phage cocktail composed of 10 phages, we observed delayed emergence of phage-resistant bacteria in Kp21 but not in Kp22. Intriguingly, phage-resistant Kp21 exhibited heightened sensitivity to other bacteriophages, indicating a "trade-off" for resistance to phage øKp_21. Our proposed phage set has an adequate number of phages to combat the K. pneumoniae strain prevalent in Japan, underscoring the potential of a well-designed phage cocktail in mitigating the occurrence of phage-resistant bacteria. IMPORTANCE The emergence of Klebsiella pneumoniae harboring the bla IMP-6 plasmid poses an escalating threat in Japan. In this study, we found 29 newly isolated bacteriophages that infect K. pneumoniae strains carrying the pKPI-6 plasmid from clinical settings in western Japan. Our phages exhibited a broad host range. We applied a phage cocktail treatment composed of 10 phages against two host strains, Kp21 and Kp22, which displayed varying phage susceptibility patterns. Although the phage cocktail delayed the emergence of phage-resistant Kp21, it was unable to hinder the emergence of phage-resistant Kp22. Moreover, the phage-resistant Kp21 became sensitive to other phages that were originally non-infective to the wild-type Kp21 strains. Our study highlights the potential of a well-tailored phage cocktail in reducing the occurrence of phage-resistant bacteria.

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.

Characterization and Comparative Genomic Analysis of a Deep-Sea Bacillus Phage Reveal a Novel Genus

As the most abundant biological entities, viruses are the major players in marine ecosystems. However, our knowledge on virus diversity and virus-host interactions in the deep sea remains very limited. In this study, vB_BteM-A9Y, a novel bacteriophage infecting Bacillus tequilensis, was isolated from deep-sea sediments in the South China Sea. vB_BteM-A9Y has a hexametric head and a long, complex contractile tail, which are typical features of myophages. vB_BteM-A9Y initiated host lysis at 60 min post infection with a burst size of 75 PFU/cell. The phage genome comprises 38,634 base pairs and encodes 54 predicted open reading frames (ORFs), of which 27 ORFs can be functionally annotated by homology analysis. Interestingly, abundant ORFs involved in DNA damage repair were identified in the phage genome, suggesting that vB_BteM-A9Y encodes multiple pathways for DNA damage repair, which may help to maintain the stability of the host/phage genome. A BLASTn search of the whole genome sequence of vB_BteM-A9Y against the GenBank revealed no existing homolog. Consistently, a phylogenomic tree and proteome-based phylogenetic tree analysis showed that vB_BteM-A9Y formed a unique branch. Further comparative analysis of genomic nucleotide similarity and ORF homology of vB_BteM-A9Y with its mostly related phages showed that the intergenomic similarity between vB_BteM-A9Y and these phages was 0-33.2%. Collectively, based on the comprehensive morphological, phylogenetic, and comparative genomic analysis, we propose that vB_BteM-A9Y belongs to a novel genus under Caudoviricetes. Therefore, our study will increase our knowledge on deep-sea virus diversity and virus-host interactions, as well as expanding our knowledge on phage taxonomy.

Phage-microbe dynamics after sterile faecal filtrate transplantation in individuals with metabolic syndrome: a double-blind, randomised, placebo-controlled clinical trial assessing efficacy and safety

Bacteriophages (phages) are bacterial viruses that have been shown to shape microbial communities. Previous studies have shown that faecal virome transplantation can decrease weight gain and normalize blood glucose tolerance in diet-induced obese mice. Therefore, we performed a double-blind, randomised, placebo-controlled pilot study in which 24 individuals with metabolic syndrome were randomised to a faecal filtrate transplantation (FFT) from a lean healthy donor (n = 12) or placebo (n = 12). The primary outcome, change in glucose metabolism, and secondary outcomes, safety and longitudinal changes within the intestinal bacteriome and phageome, were assessed from baseline up to 28 days. All 24 included subjects completed the study and are included in the analyses. While the overall changes in glucose metabolism are not significantly different between both groups, the FFT is well-tolerated and without any serious adverse events. The phage virion composition is significantly altered two days after FFT as compared to placebo, which coincides with more virulent phage-microbe interactions. In conclusion, we provide evidence that gut phages can be safely administered to transiently alter the gut microbiota of recipients.

An extended catalog of integrated prophages in the infant and adult fecal microbiome shows high prevalence of lysogeny

Background and aims

The acquisition and gradual maturation of gut microbial communities during early childhood is central to an individual's healthy development. Bacteriophages have the potential to shape the gut bacterial communities. However, the complex ecological interactions between phages and their bacterial host are still poorly characterized. In this study, we investigated the abundance and diversity of integrated prophages in infant and adult gut bacteria by detecting integrated prophages in metagenome assembled genomes (MAGs) of commensal bacteria.

Methods

Our study included 88 infants sampled at 3 weeks, 3 months, 6 months, and 12 months (n = 323 total samples), and their parents around delivery time (n = 138 total samples). Fecal DNA was extracted and characterized by using shotgun metagenomic sequencing, and a collection of prokaryotic MAGs was generated. The MAG collection was screened for the presence of integrated bacteriophage sequences, allowing their taxonomic and functional characterization.

Results

A large collection of 6,186 MAGs from infant and adult gut microbiota was obtained and screened for integrated prophages, allowing the identification of 7,165 prophage sequences longer than 10 kb. Strikingly, more than 70% of the near-complete MAGs were identified as lysogens. The prevalence of prophages in MAGs varied across bacterial families, with a lower prevalence observed among Coriobacteriaceae, Eggerthellaceae, Veillonellaceae and Burkholderiaceae, while a very high prevalence of lysogen MAGs were observed in Oscillospiraceae, Enterococcaceae, and Enterobacteriaceae. Interestingly for several bacterial families such as Bifidobacteriaceae and Bacteroidaceae, the prevalence of prophages in MAGs was higher in early infant time point (3 weeks and 3 months) than in later sampling points (6 and 12 months) and in adults. The prophage sequences were clustered into 5,616 species-like vOTUs, 77% of which were novel. Finally, we explored the functional repertoire of the potential auxiliary metabolic genes carried by these prophages, encoding functions involved in carbohydrate metabolism and degradation, amino acid metabolism and carbon metabolism.

Conclusion

Our study provides an enhanced understanding of the diversity and prevalence of lysogens in infant and adult gut microbiota and suggests a complex interplay between prophages and their bacterial hosts.

A predicted CRISPR-mediated symbiosis between uncultivated archaea

CRISPR-Cas systems defend prokaryotic cells from invasive DNA of viruses, plasmids and other mobile genetic elements. Here, we show using metagenomics, metatranscriptomics and single-cell genomics that CRISPR systems of widespread, uncultivated archaea can also target chromosomal DNA of archaeal episymbionts of the DPANN superphylum. Using meta-omics datasets from Crystal Geyser and Horonobe Underground Research Laboratory, we find that CRISPR spacers of the hosts Candidatus Altiarchaeum crystalense and Ca. A. horonobense, respectively, match putative essential genes in their episymbionts' genomes of the genus Ca. Huberiarchaeum and that some of these spacers are expressed in situ. Metabolic interaction modelling also reveals complementation between host-episymbiont systems, on the basis of which we propose that episymbionts are either parasitic or mutualistic depending on the genotype of the host. By expanding our analysis to 7,012 archaeal genomes, we suggest that CRISPR-Cas targeting of genomes associated with symbiotic archaea evolved independently in various archaeal lineages.

Genome Characterization of the Novel Lytic Enterobacter cloacae Phage vB_EclM_Q7622

Enterobacter cloacae is a widespread opportunistic pathogen that causes urinary tract infection. The abuse of antibiotics enabled multidrug-resistant strains to spread. Bacteriophage therapy is a naturally, safe, and efficient alternative treatment technology against multi-resistant bacteria. In this study, a virulent phage vB_EclM_Q7622 (Q7622) was isolated from the sewage of Jiangcun poultry market in Guangzhou city. Transmission electron microscopy indicated that Q7622 had an icosahedral head (97.8 ± 5.6 nm in diameter) and a short, contractile tail (113.7 ± 4.5 nm). Its double-stranded DNA genome is composed of 173,871 bp with a GC content of 40.02%. It possesses 297 open reading frames and 9 tRNAs. No known virulence and resistance genes were detected, indicated that phage Q7622 could be used for pathogens prevention and control safely. Comparative genomic and phylogenetic analysis showed that Q7622 was highly similar to the phages vB_EclM_CIP9 and vB_EhoM-IME523. The highest nucleotide similarity between Q7622 and the similar phages in NCBI calculated by pyANI and VIRIDIC were 94.9% and 89.1% with vB_EhoM-IME523 respectively, below 95%. Thus, according to the result of nucleotide similarity calculation results, Q7622 was a novel virulent Enterobacter cloacae phage strain of the genus Kanagawavirus.

Comparative genomic analysis of alloherpesviruses: Exploring an available genus/species demarcation proposal and method

The family Alloherpesviridae contains herpesviruses of fish and amphibians. Due to the significant economic losses to aquaculture that herpesviruses can cause, the primary areas of research interest are concerning their pathogenesis and prevention. Despite alloherpesvirus genomic sequences becoming more widely accessible, methods regarding their genus/species classification are still relatively unexplored. In the present study, the phylogenetic relationships between 40 completely sequenced alloherpesviruses were illustrated by the viral proteomic tree (ViPTree), which was divided into three monophyletic groups, namely Cyprinivirus, Ictalurivirus and Batrachovirus. Additionally, average nucleotide identity (ANI) and average amino acid identity (AAI) analyses were performed across all available sequences and clearly displayed species boundaries with the threshold value of ANI/AAI set at 90%. Subsequently, core-pan analysis uncovered 809 orthogroups and 11 core genes shared by all 40 alloherpesvirus genome sequences. For the former, a 15 percent identity depicts a clear genus boundary; for the latter, 8 of them may be qualified for phylogenetic analysis based on amino acid or nucleic acid sequences after being verified using maximum likelihood (ML) or neighbor-joining (NJ) phylogenetic trees. Finally, although the dot plot analysis was valid for the members within Ictalurivirus, it was unsuccessful for Cyprinivirus and Batrachovirus. Taken together, the comparison of individual methodologies provides a wide range of alternatives for alloherpesviruses classification under various circumstances.

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.

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.

MArVD2: a machine learning enhanced tool to discriminate between archaeal and bacterial viruses in viral datasets

Our knowledge of viral sequence space has exploded with advancing sequencing technologies and large-scale sampling and analytical efforts. Though archaea are important and abundant prokaryotes in many systems, our knowledge of archaeal viruses outside of extreme environments is limited. This largely stems from the lack of a robust, high-throughput, and systematic way to distinguish between bacterial and archaeal viruses in datasets of curated viruses. Here we upgrade our prior text-based tool (MArVD) via training and testing a random forest machine learning algorithm against a newly curated dataset of archaeal viruses. After optimization, MArVD2 presented a significant improvement over its predecessor in terms of scalability, usability, and flexibility, and will allow user-defined custom training datasets as archaeal virus discovery progresses. Benchmarking showed that a model trained with viral sequences from the hypersaline, marine, and hot spring environments correctly classified 85% of the archaeal viruses with a false detection rate below 2% using a random forest prediction threshold of 80% in a separate benchmarking dataset from the same habitats.

Comparative Genomics of Staphylococcus rostri, an Undescribed Bacterium Isolated from Dairy Mastitis

This study characterizes 81 S. rostri isolates from bovine mastitis (of which 80 were subclinical). The isolates were first identified as S. microti by MALDI-TOF MS, but later whole genome sequencing analysis allowed reclassification as S. rostri. The isolates were derived from 52 cows and nine dairy herds in Denmark. To describe the pathogenicity of S. rostri, we used whole genome sequencing to infer the distribution of genes associated with virulence, antibiotic resistance, and mobile genetic elements. Also, we performed a core-genome phylogeny analysis to study the genetic relatedness among the isolates. All 81 isolates expressed the same virulence profile comprising two putative virulence genes, clpP and clpC. Three isolates carried a resistance gene encoding streptomycin (str) or lincomycin (lnuA) resistance. The distribution of plasmids suggested the detected antibiotic resistance genes to be plasmid-mediated. Phages were abundant among the isolates, and the single isolate from clinical mastitis acquired a phage disparate from the rest, which potentially could be involved with virulence in S. rostri. The core genome phylogeny revealed a strong genetic intra-herd conservation, which indicates the source of introduction being herd-specific and might further imply the ability of S. rostri to adapt to the bovine niche and spread from cow-to-cow in a contagious manner. With this study, we aim to acquaint clinicians and professionals with the existence of S. rostri which might have been overlooked so far.

Characterisation of the Novel Filamentous Phage PMBT54 Infecting the Milk Spoilage Bacteria Pseudomonas carnis and Pseudomonas lactis

Filamentous bacteriophages are lysogenic and pseudo-lysogenic viruses that do not lyse their host but are often continuously secreted from the infected cell. They belong to the order Tubulavirales, which encompasses three families, with the Inoviridae being the largest. While the number of identified inoviral sequences has greatly increased in recent years due to metagenomic studies, morphological and physiological characterisation is still restricted to only a few members of the filamentous phages. Here, we describe the novel filamentous phage PMBT54, which infects the spoilage-relevant Pseudomonas species P. carnis and P. lactis. Its genome is 7320 bp in size, has a mol% GC content of 48.37, and codes for 13 open-reading frames, two of which are located on the (-) strand. The virion exhibits a typical filamentous morphology and is secreted from the host cell at various lengths. The phage was shown to promote biofilm formation in both host strains and, therefore, has potential implications for milk spoilage, as biofilms are a major concern in the dairy industry.

Genomic diversity and biogeographic distributions of a novel lineage of bacteriophages that infect marine OM43 bacteria

The marine methylotrophic OM43 clade is considered an important bacterial group in coastal microbial communities. OM43 bacteria, which are closely related to phytoplankton blooms, have small cell sizes and streamlined genomes. Bacteriophages profoundly shape the evolutionary trajectories, population dynamics, and physiology of microbes. The prevalence and diversity of several phages that infect OM43 bacteria have been reported. In this study, we isolated and sequenced two novel OM43 phages, MEP401 and MEP402. These phages share 90% of their open reading frames (ORFs) and are distinct from other known phage isolates. Furthermore, a total of 99 metagenomic viral genomes (MVGs) closely related to MEP401 and MEP402 were identified. Phylogenomic analyses suggest that MEP401, MEP402, and these identified MVGs belong to a novel subfamily in the family Zobellviridae and that they can be separated into two groups. Group I MVGs show conserved whole-genome synteny with MEP401, while group II MVGs possess the MEP401-type DNA replication module and a distinct type of morphogenesis and packaging module, suggesting that genomic recombination occurred between phages. Most members in these two groups were predicted to infect OM43 bacteria. Metagenomic read-mapping analysis revealed that the phages in these two groups are globally ubiquitous and display distinct biogeographic distributions, with some phages being predominant in cold regions, some exclusively detected in estuarine stations, and others displaying wider distributions. This study expands our knowledge of the diversity and ecology of a novel phage lineage that infects OM43 bacteria by describing their genomic diversity and global distribution patterns. IMPORTANCE OM43 phages that infect marine OM43 bacteria are important for host mortality, community structure, and physiological functions. In this study, two OM43 phages were isolated and characterized. Metagenomic viral genome (MVG) retrieval using these two OM43 phages as baits led to the identification of two phage groups of a new subfamily in the family Zobellviridae. We found that group I MVGs share similar genomic content and arrangement with MEP401 and MEP402, whereas group II MVGs only possess the MEP401-type DNA replication module. Metagenomic mapping analysis suggests that members in these two groups are globally ubiquitous with distinct distribution patterns. This study provides important insights into the genomic diversity and biogeography of the OM43 phages in the global ocean.

Psychrobacter Phage Encoding an Antibiotics Resistance Gene Represents a Novel Caudoviral Family

Psychrobacter is an important bacterial genus that is widespread in Antarctic and marine environments. However, to date, only two complete Psychrobacter phage sequences have been deposited in the NCBI database. Here, the novel Psychrobacter phage vB_PmaS_Y8A, infecting Psychrobacter HM08A, was isolated from sewage in the Qingdao area, China. The morphology of vB_PmaS_Y8A was characterized by transmission electron microscopy, revealing an icosahedral head and long tail. The genomic sequence of vB_PmaS_Y8A is linear, double-stranded DNA with a length of 40,226 bp and 44.1% G+C content, and encodes 69 putative open reading frames. Two auxiliary metabolic genes (AMGs) were identified, encoding phosphoadenosine phosphosulfate reductase and MarR protein. The first AMG uses thioredoxin as an electron donor for the reduction of phosphoadenosine phosphosulfate to phosphoadenosine phosphate. MarR regulates multiple antibiotic resistance mechanisms in Escherichia coli and is rarely found in viruses. No tRNA genes were identified and no lysogeny-related feature genes were detected. However, many similar open reading frames (ORFs) were found in the host genome, which may indicate that Y8A also has a lysogenic stage. Phylogenetic analysis based on the amino acid sequences of whole genomes and comparative genomic analysis indicate that vB_PmaS_Y8A contains a novel genomic architecture similar only to that of Psychrobacter phage pOW20-A, although at a low similarity. vB_PmaS_Y8A represents a new family-level virus cluster with 22 metagenomic assembled viral genomes, here named Minviridae. IMPORTANCE Although Psychrobacter is a well-known and important bacterial genus that is widespread in Antarctic and marine environments, genetic characterization of its phages is still rare. This study describes a novel Psychrobacter phage containing an uncharacterized antibiotic resistance gene and representing a new virus family, Minviridae. The characterization provided here will bolster current understanding of genomes, diversity, evolution, and phage-host interactions in Psychrobacter populations.

Long-term CRISPR locus dynamics and stable host-virus co-existence in subsurface fractured shales

Viruses are the most ubiquitous biological entities on Earth. Even so, elucidating the impact of viruses on microbial communities and associated ecosystem processes often requires identification of unambiguous host-virus linkages-an undeniable challenge in many ecosystems. Subsurface fractured shales present a unique opportunity to first make these strong linkages via spacers in CRISPR-Cas arrays and subsequently reveal complex long-term host-virus dynamics. Here, we sampled two replicated sets of fractured shale wells for nearly 800 days, resulting in 78 metagenomes from temporal sampling of six wells in the Denver-Julesburg Basin (Colorado, USA). At the community level, there was strong evidence for CRISPR-Cas defense systems being used through time and likely in response to viral interactions. Within our host genomes, represented by 202 unique MAGs, we also saw that CRISPR-Cas systems were widely encoded. Together, spacers from host CRISPR loci facilitated 2,110 CRISPR-based viral linkages across 90 host MAGs spanning 25 phyla. We observed less redundancy in host-viral linkages and fewer spacers associated with hosts from the older, more established wells, possibly reflecting enrichment of more beneficial spacers through time. Leveraging temporal patterns of host-virus linkages across differing well ages, we report how host-virus co-existence dynamics develop and converge through time, possibly reflecting selection for viruses that can evade host CRISPR-Cas systems. Together, our findings shed light on the complexities of host-virus interactions as well as long-term dynamics of CRISPR-Cas defense among diverse microbial populations.

Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines

BACKGROUND

Climate change threatens Earth's ice-based ecosystems which currently offer archives and eco-evolutionary experiments in the extreme. Arctic cryopeg brine (marine-derived, within permafrost) and sea ice brine, similar in subzero temperature and high salinity but different in temporal stability, are inhabited by microbes adapted to these extreme conditions. However, little is known about their viruses (community composition, diversity, interaction with hosts, or evolution) or how they might respond to geologically stable cryopeg versus fluctuating sea ice conditions.

RESULTS

We used long- and short-read viromics and metatranscriptomics to study viruses in Arctic cryopeg brine, sea ice brine, and underlying seawater, recovering 11,088 vOTUs (~species-level taxonomic unit), a 4.4-fold increase of known viruses in these brines. More specifically, the long-read-powered viromes doubled the number of longer (≥25 kb) vOTUs generated and recovered more hypervariable regions by >5-fold compared to short-read viromes. Distribution assessment, by comparing to known viruses in public databases, supported that cryopeg brine viruses were of marine origin yet distinct from either sea ice brine or seawater viruses, while 94% of sea ice brine viruses were also present in seawater. A virus-encoded, ecologically important exopolysaccharide biosynthesis gene was identified, and many viruses (~half of metatranscriptome-inferred "active" vOTUs) were predicted as actively infecting the dominant microbial genera Marinobacter and Polaribacter in cryopeg and sea ice brines, respectively. Evolutionarily, microdiversity (intra-species genetic variations) analyses suggested that viruses within the stable cryopeg brine were under significantly lower evolutionary pressures than those in the fluctuating sea ice environment, while many sea ice brine virus-tail genes were under positive selection, indicating virus-host co-evolutionary arms races.

CONCLUSIONS

Our results confirmed the benefits of long-read-powered viromics in understanding the environmental virosphere through significantly improved genomic recovery, expanding viral discovery and the potential for biological inference. Evidence of viruses actively infecting the dominant microbes in subzero brines and modulating host metabolism underscored the potential impact of viruses on these remote and underexplored extreme ecosystems. Microdiversity results shed light on different strategies viruses use to evolve and adapt when extreme conditions are stable versus fluctuating. Together, these findings verify the value of long-read-powered viromics and provide foundational data on viral evolution and virus-microbe interactions in Earth's destabilized and rapidly disappearing cryosphere. Video Abstract.

Presence and role of viruses in anaerobic digestion of food waste under environmental variability

BACKGROUND

The interaction among microorganisms in the anaerobic digestion of food waste (ADFW) reactors lead to the degradation of organics and the recycling of energy. Viruses are an important component of the microorganisms involved in ADFW, but are rarely investigated. Furthermore, little is known about how viruses affect methanogenesis.

RESULTS

Thousands of viral sequences were recovered from five full-scale ADFW reactors. Gene-sharing networks indicated that the ADFW samples contained substantial numbers of unexplored anaerobic-specific viruses. Moreover, the viral communities in five full-scale reactors exhibited both commonalities and heterogeneities. The lab-scale dynamic analysis of typical ADFW scenarios suggested that the viruses had similar kinetic characteristics to their prokaryotic hosts. By associating with putative hosts, a majority of the bacteria and archaea phyla were found to be infected by viruses. Viruses may influence prokaryotic ecological niches, and thus methanogenesis, by infecting key functional microorganisms, such as sulfate-reducing bacteria (SRB), syntrophic acetate-oxidizing bacteria (SAOB), and methanogens. Metabolic predictions for the viruses suggested that they may collaborate with hosts at key steps of sulfur and long-chain fatty acid (LCFA) metabolism and could be involved in typical methanogenesis pathways to participate in methane production.

CONCLUSIONS

Our results expanded the diversity of viruses in ADFW systems and suggested two ways that viral manipulated ADFW biochemical processes. Video Abstract.

Genetic characteristics and integration specificity of Salmonella enterica temperate phages

Introduction

Temperate phages can engage in the horizontal transfer of functional genes to their bacterial hosts. Thus, their genetic material becomes an intimate part of bacterial genomes and plays essential roles in bacterial mutation and evolution. Specifically, temperate phages can naturally transmit genes by integrating their genomes into the bacterial host genomes via integrases. Our previous study showed that Salmonella enterica contains the largest number of temperate phages among all publicly available bacterial species. S. enterica is an important pathogen that can cause serious systemic infections and even fatalities.

Methods

Initially, we extracted all S. enterica temperate phages from the extensively developed temperate phage database established in our previous study. Subsequently, we conducted an in-depth analysis of the genetic characteristics and integration specificity exhibited by these S. enterica temperate phages.

Results

Here we identified 8,777 S. enterica temperate phages, all of which have integrases in their genomes. We found 491 non-redundant S. enterica temperate phage integrases (integrase entries). S. enterica temperate phage integrases were classified into three types: intA, intS, and phiRv2. Correlation analysis showed that the sequence lengths of S. enterica integrase and core regions of attB and attP were strongly correlated. Further phylogenetic analysis and taxonomic classification indicated that both the S. enterica temperate phage genomes and the integrase gene sequences were of high diversities.

Discussion

Our work provides insight into the essential integration specificity and genetic diversity of S. enterica temperate phages. This study paves the way for a better understanding of the interactions between phages and S. enterica. By analyzing a large number of S. enterica temperate phages and their integrases, we provide valuable insights into the genetic diversity and prevalence of these elements. This knowledge has important implications for developing targeted therapeutic interventions, such as phage therapy, to combat S. enterica infections. By harnessing the lytic capabilities of temperate phages, they can be engineered or utilized in phage cocktails to specifically target and eradicate S. enterica strains, offering an alternative or complementary approach to traditional antibiotic treatments. Our study has implications for public health and holds potential significance in combating clinical infections caused by S. enterica.

ViWrap: A modular pipeline to identify, bin, classify, and predict viral-host relationships for viruses from metagenomes

Viruses are increasingly being recognized as important components of human and environmental microbiomes. However, viruses in microbiomes remain difficult to study because of the difficulty in culturing them and the lack of sufficient model systems. As a result, computational methods for identifying and analyzing uncultivated viral genomes from metagenomes have attracted significant attention. Such bioinformatics approaches facilitate the screening of viruses from enormous sequencing datasets originating from various environments. Though many tools and databases have been developed for advancing the study of viruses from metagenomes, there is a lack of integrated tools enabling a comprehensive workflow and analyses platform encompassing all the diverse segments of virus studies. Here, we developed ViWrap, a modular pipeline written in Python. ViWrap combines the power of multiple tools into a single platform to enable various steps of virus analysis, including identification, annotation, genome binning, species- and genus-level clustering, assignment of taxonomy, prediction of hosts, characterization of genome quality, comprehensive summaries, and intuitive visualization of results. Overall, ViWrap enables a standardized and reproducible pipeline for both extensive and stringent characterization of viruses from metagenomes, viromes, and microbial genomes. Our approach has flexibility in using various options for diverse applications and scenarios, and its modular structure can be easily amended with additional functions as necessary. ViWrap is designed to be easily and widely used to study viruses in human and environmental systems. ViWrap is publicly available via GitHub (https://github.com/AnantharamanLab/ViWrap). A detailed description of the software, its usage, and interpretation of results can be found on the website.

Astrovirology: how viruses enhance our understanding of life in the Universe

Viruses are the most numerically abundant biological entities on Earth. As ubiquitous replicators of molecular information and agents of community change, viruses have potent effects on the life on Earth, and may play a critical role in human spaceflight, for life-detection missions to other planetary bodies and planetary protection. However, major knowledge gaps constrain our understanding of the Earth's virosphere: (1) the role viruses play in biogeochemical cycles, (2) the origin(s) of viruses and (3) the involvement of viruses in the evolution, distribution and persistence of life. As viruses are the only replicators that span all known types of nucleic acids, an expanded experimental and theoretical toolbox built for Earth's viruses will be pivotal for detecting and understanding life on Earth and beyond. Only by filling in these knowledge and technical gaps we will obtain an inclusive assessment of how to distinguish and detect life on other planetary surfaces. Meanwhile, space exploration requires life-support systems for the needs of humans, plants and their microbial inhabitants. Viral effects on microbes and plants are essential for Earth's biosphere and human health, but virus-host interactions in spaceflight are poorly understood. Viral relationships with their hosts respond to environmental changes in complex ways which are difficult to predict by extrapolating from Earth-based proxies. These relationships should be studied in space to fully understand how spaceflight will modulate viral impacts on human health and life-support systems, including microbiomes. In this review, we address key questions that must be examined to incorporate viruses into Earth system models, life-support systems and life detection. Tackling these questions will benefit our efforts to develop planetary protection protocols and further our understanding of viruses in astrobiology.

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.

Diversity, Relationship, and Distribution of Virophages and Large Algal Viruses in Global Ocean Viromes

Virophages are a group of small double-stranded DNA viruses that replicate and proliferate with the help of the viral factory of large host viruses. They are widely distributed in aquatic environments but are more abundant in freshwater ecosystems. Here, we mined the Global Ocean Viromes 2.0 (GOV 2.0) dataset for the diversity, distribution, and association of virophages and their potential host large viruses in marine environments. We identified 94 virophage sequences (>5 kbp in length), of which eight were complete genomes. The MCP phylogenetic tree showed that the GOV virophages were widely distributed on the global virophage tree but relatively clustered on three major branches. The gene-sharing network divided GOV virophages into 21 outliers, 2 overlaps, and 14 viral clusters, of which 4 consisted of only the GOV virophages. We also identified 45 large virus sequences, 8 of which were >100 kbp in length and possibly involved in cell-virus-virophage (C-V-v) trisome relationships. The potential eukaryotic hosts of these eight large viruses and the eight virophages with their complete genomes identified are likely to be algae, based on comparative genomic analysis. Both homologous gene and codon usage analyses support a possible interaction between a virophage (GOVv18) and a large algal virus (GOVLV1). These results indicate that diverse and novel virophages and large viruses are widespread in global marine environments, suggesting their important roles and the presence of complicated unknown C-V-v relationships in marine ecosystems.