The genomic content and context of auxiliary metabolic genes in roseophages

Novel phage infecting the Roseobacter CHUG lineage reveals a diverse and globally distributed phage family

Bacteriophages play an essential role in shaping the diversity and metabolism of bacterial communities. Marine Roseobacter group is an abundant heterotrophic bacterial group that is involved in many major element cycles, especially carbon and sulfur. Members of the Roseobacter CHUG (Clade Hidden and Underappreciated Globally) lineage are globally distributed and are activated in pelagic marine environments. In this study, we isolated and characterized a phage, CRP-810, that infects the CHUG strain FZCC0198. The genome of CRP-810 was dissimilar to those of other known phages. Additionally, 251 uncultured viral genomes (UViGs) closely related to CRP-810 were obtained from the uncultivated marine viral contig databases. Comparative genomic and phylogenetic analyses revealed that CRP-810 and these related UViGs exhibited conserved genome synteny, representing a new phage family with at least eight subgroups. Most of the CRP-810-type phages contain an integrase gene, and CRP-810 can be integrated into the host genome. Further analysis revealed that three CRP-810-type members were prophages found in the genomes of marine SAR11, Poseidonocella, and Sphingomonadaceae. Finally, viromic read-mapping analysis showed that CRP-810-type phages were globally distributed and displayed distinct biogeographic patterns related to temperature and latitude. Many members with a lower G + C content were mainly distributed in the trade station, whereas members with a higher G + C content were mainly distributed in polar and westerlies station, indicating that the niche differentiation of phages was subject to host adaptation. Collectively, these findings identify a novel phage family and expand our understanding of phylogenetic diversity, evolution, and biogeography of marine phages.

IMPORTANCE

The Roseobacter CHUG lineage, affiliated with the Pelagic Roseobacter Cluster (PRC), is widely distributed in the global oceans and is active in oligotrophic seawater. However, knowledge of the bacteriophages that infect CHUG members is limited. In this study, a CHUG phage, CRP-810, that infects the CHUG strain FZCC0198, was isolated and shown to have a novel genomic architecture. In addition, 251 uncultured viral genomes closely related to CRP-810 were recovered and included in the analyses. Phylogenomic analyses revealed that the CRP-810-type phages represent a new phage family containing at least eight genus-level subgroups. Members of this family were predicted to infect various marine bacteria. We also demonstrated that the CRP-810-type phages are widely distributed in global oceans and display distinct biogeographic patterns related to latitude. Collectively, this study provides important insights into the genomic organization, diversity, and ecology of a novel phage family that infect ecologically important bacteria in the global ocean.

Transfer of modified gut viromes improves symptoms associated with metabolic syndrome in obese male mice

Metabolic syndrome encompasses amongst other conditions like obesity and type-2 diabetes and is associated with gut microbiome (GM) dysbiosis. Fecal microbiota transplantation (FMT) has been explored to treat metabolic syndrome by restoring the GM; however, concerns on accidentally transferring pathogenic microbes remain. As a safer alternative, fecal virome transplantation (FVT, sterile-filtrated feces) has the advantage over FMT in that mainly bacteriophages are transferred. FVT from lean male donors have shown promise in alleviating the metabolic effects of high-fat diet in a preclinical mouse study. However, FVT still carries the risk of eukaryotic viral infections. To address this, recently developed methods are applied for removing or inactivating eukaryotic viruses in the viral component of FVT. Modified FVTs are compared with unmodified FVT and saline in a diet-induced obesity model on male C57BL/6 N mice. Contrasted with obese control, mice administered a modified FVT (nearly depleted for eukaryotic viruses) exhibits enhanced blood glucose clearance but not weight loss. The unmodified FVT improves liver pathology and reduces the proportions of immune cells in the adipose tissue with a non-uniform response. GM analysis suggests that bacteriophage-mediated GM modulation influences outcomes. Optimizing these approaches could lead to the development of safe bacteriophage-based therapies targeting metabolic syndrome through GM restoration.

High niche specificity and host genetic diversity of groundwater viruses

Viruses are key members of microbial communities that exert control over host abundance and metabolism, thereby influencing ecosystem processes and biogeochemical cycles. Aquifers are known to host taxonomically diverse microbial life, yet little is known about viruses infecting groundwater microbial communities. Here, we analysed 16 metagenomes from a broad range of groundwater physicochemistries. We recovered 1571 viral genomes that clustered into 468 high-quality viral operational taxonomic units. At least 15% were observed to be transcriptionally active, although lysis was likely constrained by the resource-limited groundwater environment. Most were unclassified (95%), and the remaining 5% were Caudoviricetes. Comparisons with viruses inhabiting other aquifers revealed no shared species, indicating substantial unexplored viral diversity. In silico predictions linked 22.4% of the viruses to microbial host populations, including to ultra-small prokaryotes, such as Patescibacteria and Nanoarchaeota. Many predicted hosts were associated with the biogeochemical cycling of carbon, nitrogen, and sulfur. Metabolic predictions revealed the presence of 205 putative auxiliary metabolic genes, involved in diverse processes associated with the utilization of the host's intracellular resources for biosynthesis and transformation reactions, including those involved in nucleotide sugar, glycan, cofactor, and vitamin metabolism. Viruses, prokaryotes overall, and predicted prokaryotic hosts exhibited narrow spatial distributions, and relative abundance correlations with the same groundwater parameters (e.g. dissolved oxygen, nitrate, and iron), consistent with host control over viral distributions. Results provide insights into underexplored groundwater viruses, and indicate the large extent to which viruses may manipulate microbial communities and biogeochemistry in the terrestrial subsurface.

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.

Structure and proposed DNA delivery mechanism of a marine roseophage

Tailed bacteriophages (order, Caudovirales) account for the majority of all phages. However, the long flexible tail of siphophages hinders comprehensive investigation of the mechanism of viral gene delivery. Here, we report the atomic capsid and in-situ structures of the tail machine of the marine siphophage, vB_DshS-R4C (R4C), which infects Roseobacter. The R4C virion, comprising 12 distinct structural protein components, has a unique five-fold vertex of the icosahedral capsid that allows genome delivery. The specific position and interaction pattern of the tail tube proteins determine the atypical long rigid tail of R4C, and further provide negative charge distribution within the tail tube. A ratchet mechanism assists in DNA transmission, which is initiated by an absorption device that structurally resembles the phage-like particle, RcGTA. Overall, these results provide in-depth knowledge into the intact structure and underlining DNA delivery mechanism for the ecologically important siphophages.

The persistence and stabilization of auxiliary genes in the human skin virome

BACKGROUND

The human skin contains a diverse microbiome that provides protective functions against environmental pathogens. Studies have demonstrated that bacteriophages modulate bacterial community composition and facilitate the transfer of host-specific genes, potentially influencing host cellular functions. However, little is known about the human skin virome and its role in human health. Especially, how viral-host relationships influence skin microbiome structure and function is poorly understood.

RESULTS

Population dynamics and genetic diversity of bacteriophage communities in viral metagenomic data collected from three anatomical skin locations from 60 subjects at five different time points revealed that cutaneous bacteriophage populations are mainly composed of tailed Caudovirales phages that carry auxiliary genes to help improve metabolic remodeling to increase bacterial host fitness through antimicrobial resistance. Sequence variation in the MRSA associated antimicrobial resistance gene, erm(C) was evaluated using targeted sequencing to further confirm the presence of antimicrobial resistance genes in the human virome and to demonstrate how functionality of such genes may influence persistence and in turn stabilization of bacterial host and their functions.

CONCLUSIONS

This large temporal study of human skin associated viruses indicates that the human skin virome is associated with auxiliary metabolic genes and antimicrobial resistance genes to help increase bacterial host fitness.

Genetic diversity of virus auxiliary metabolism genes associated with phosphorus metabolism in Napahai plateau wetland

Viruses play important roles in ecosystems by interfering with the central metabolic pathways of the host during infection via the expression of auxiliary metabolic genes (AMGs), altering the productivity of ecosystems and thus affecting geochemical cycling. In this study, the genetic diversity of phosphorus metabolism AMGs phoH, phoU and pstS was investigated by phylogenetic analysis, PCoA analysis, and alpha diversity analysis based on metagenomic data. It was found that the majority of the sequences were unique to Napahai plateau wetland. It was shown that the genetic diversity of phoH, phoU and pstS genes was independent of both habitats and host origins. In addition, the metabolic pathway of AMGs associated with the phosphorus cycling was identified based on metagenomic data. When phosphorus is deficient, virus utilizes AMGs to affect the metabolic pathway, contributing to higher phosphorus levels in the host and facilitating virus survival, replication, and propagation in the host cell.

Biogeochemical sulfur cycling of virus auxiliary metabolic genes involved in Napahai plateau wetland

Virus plays important roles in regulating microbial community structure, horizontal gene transfer, and promoting biological evolution, also augmenting host metabolism during infection via the expression of auxiliary metabolic genes (AMGs), and thus affect biogeochemical cycling in the oceans. As the "kidney of the earth," wetlands have rich biodiversity and abundant resources. Based on metagenomic data, 10 AMGs associated with sulfur cycling, i.e., tusA, moaD, dsrE, soxA, soxB, soxC, soxD, soxX, soxY, and soxZ, were analyzed in Napahai plateau wetland. The phylogenetic trees of AMGs involved in sulfur metabolism from different habitats and host origins were constructed. Combined with principal coordinate analysis, it revealed that most AMGs associated with sulfur metabolism clustered separately, indicating the abundance and uniqueness in this region. The sulfur metabolism pathways involved by AMGs were mainly SOX systems, among which sulfur oxidation was associated with moaD and dsrE genes, while sulfur transport was related to tusA genes. It provides an insight into the biogeochemical sulfur cycling in plateau wetlands and lays the foundation for further study on the co-evolution of virus and host.

Structural analysis of the housecleaning nucleoside triphosphate pyrophosphohydrolase MazG from Mycobacterium tuberculosis

The housecleaning enzyme of Mycobacterium tuberculosis (Mtb), MazG, is a nucleoside triphosphate pyrophosphohydrolase (NTP-PPase) and can hydrolyze all canonical or non-canonical NTPs into NMPs and pyrophosphate. The Mycobacterium tuberculosis MazG (Mtb-MazG) contributes to antibiotic resistance in response to oxidative or nitrosative stress under dormancy, making it a promising target for treating TB in latent infection patients. However, the structural basis of Mtb-MazG is not clear. Here we describe the crystal structure of Mtb-MazG (1-185) at 2.7 Å resolution, composed of two similar folded spherical domains in tandem. Unlike other all-α NTP pyrophosphatases, Mtb-MazG has an N-terminal extra region composed of three α-helices and five β-strands. The second domain is global, with five α-helices located in the N-terminal domain. Gel-filtration assay and SAXS analysis show that Mtb-MazG forms an enzyme-active dimer in solution. In addition, the metal ion Mg²⁺ is bound with four negative-charged residues Glu119, Glu122, Glu138, and Asp141. Different truncations and site-directed mutagenesis revealed that the full-length dimeric form and the metal ion Mg²⁺ are indispensable for the catalytic activity of Mtb-MazG. Thus, our work provides new insights into understanding the molecular basis of Mtb-MazG.

Distribution of rare N4-like viruses in temperate estuaries unveiled by viromics

The relative abundance of N4-like viruses in two temperate estuaries was assessed using four different methods, read mapping to known N4-like virus isolates, read mapping to native viral contigs, reciprocal blast search based on core genes, and read taxonomy classification using Kaiju. Overall, N4-like viruses were found to be of low abundance in the estuarine viromes. When mapping reads to only known N4-like virus genomes, high occurrences of N4-like viruses infecting Roseobacter were found, with their diversity consisting mostly of locally isolated Roseobacter N4-like virus species. Both contig-based methods and Kaiju classification showed similar seasonal patterns for N4-like viruses, and redundancy analysis revealed a negative correlation between N4-like viruses and temperature, suggesting that N4-like viruses may be more abundant in colder water. The discrepancy of relative abundance estimates using different methods indicates that N4-like viruses are best represented by native viral sequences. Our study indicates that N4-like viruses are rare in the marine environment and also provide insight into the importance of including local viral sequences in reference databases.

Ecological and functional roles of bacteriophages in contrasting environments: marine, terrestrial and human gut

While they are the most abundant biological entities on the planet, the role of bacteriophages (phages) in the microbiome remains enigmatic and understudied. With a rise in the number of metagenomics studies and the publication of highly efficient phage mining programmes, we now have extensive data on the genomic and taxonomic diversity of (mainly) DNA bacteriophages in a wide range of environments. In addition, the higher throughput and quality of sequencing is allowing for strain-level reconstructions of phage genomes from metagenomes. These factors will ultimately help us to understand the role these phages play as part of specific microbial communities, enabling the tracking of individual virus genomes through space and time. Using lessons learned from the latest metagenomic studies, we focus on two explicit aspects of the role bacteriophages play within the microbiome, their ecological role in structuring bacterial populations, and their contribution to microbiome functioning by encoding auxiliary metabolism genes.

Isolation and characterization of Brochothrix phage ADU4

B. thermosphacta is a psychrotrophic bacterium that often forms the predominant part of the spoilage microflora of aerobically and anaerobically stored meats. Bacteriophages are natural enemies of bacteria and their potential for use in environmentally friendly biocontrol of specific pathogens in food is being intensively studied. In this study, we reported the isolation and characterization of the newly isolated lytic Brochothrix phage ADU4, which is capable of infecting the B. thermosphacta bacterium. For the characterization of Brochothrix phage ADU4; host range, multiplicity of infection values (MOI), phage growth parameters (latent period and burst size), stability at various temperatures and pH, reduction growth of bacteria, effect on biofilm, and molecular characterization were investigated. The spot-test analysis showed positivity with B. thermosphacta strains, while no infection was observed in any other species and genera of bacteria tested. The optimal MOI value of the phage was determined as 0.1. The phage latent period and burst sizes were 40-50 min and 311 PFU/ml per infected host cell, respectively by one-step growth curve analysis. Brochothrix phage ADU4 reduced bacteria immediately after infection, which is shown by optical density (OD) measurement and colony counting (<10 CFU/ml) for 3 days. The degradation of B. thermosphacta in biofilm by Brochothrix phage ADU4 was analyzed and it was found that high titer phage breakdown the existing biofilm and also persistently inhibited biofilm formation. Brochothrix phage ADU4 genome was found to be 127,819 bp, and GC content 41.65%. The genome contains 217 putative open reading frames (ORFs), 4 tRNAs, and additionally, no known virulence and antibiotic resistance genes (AMR) were identified. Brochothrix phage ADU4 showed a high identity (96.09%) to the A9 phage that belongs to the Herelleviridae family. Nevertheless, the assembly module and its around appeared less conserved, and some DNA fragments in Brochothrix phage ADU4 genome were not found in A9 genome and vice versa. A9 contains TnpB, a transposase accessory protein involved in lysogenicity which is absent in Brochothrix phage ADU4. In contrary Brochothrix phage ADU4 had auxiliary metabolic genes (AMG) mostly carried by lytic phages. All these results showed that the Brochothrix phage ADU4 has excellent properties such as strong antibacterial activity, short latent period, high burst size, stability in different conditions, inhibition of biofilms, and absence of virulence and AMR genes. Based on all these features, this newly isolated phage is promising to control B. thermosphacta contamination in meat and meat products, and has the potential to be used alone or in combination with phage cocktails.

Metagenomics of Virus Diversities in Solid-State Brewing Process of Traditional Chinese Vinegar

Traditional Chinese vinegar offers an exceptional flavor and rich nutrients due to its unique solid-state fermentation process, which is a multiple microbial fermentation system including various bacteria, fungi and viruses. However, few studies on the virus diversities in traditional Chinese vinegar have been reported. In this paper, using Zhenjiang aromatic vinegar as a model system, we systemically explored the viral communities in the solid-state brewing process of traditional Chinese vinegar using bacterial and viral metagenomes. Results showed that the viral diversity in vinegar Pei was extensive and the virus communities varied along with the fermentation process. In addition, there existed some interactions between viral and bacterial communities. Moreover, abundant antibiotic resistance genes were found in viromes, indicating that viruses might protect fermentation bacteria strains from the stress of antibiotics in the fermentation environment. Remarkably, we identified abundant auxiliary carbohydrate metabolic genes (including alcohol oxidases, the key enzymes for acetic acid synthesis) from viromes, implying that viruses might participate in the acetic acid synthesis progress of the host through auxiliary metabolic genes. Taken together, our results indicated the potential roles of viruses in the vinegar brewing process and provided a new perspective for studying the fermentation mechanisms of traditional Chinese vinegar.

A snapshot of the global drinking water virome: Diversity and metabolic potential vary with residual disinfectant use

Viruses are important drivers of microbial community ecology and evolution, influencing microbial mortality, metabolism, and horizontal gene transfer. However, the effects of viruses remain largely unknown in many environments, including in drinking water systems. Drinking water metagenomic studies have offered a whole community perspective of bacterial impacts on water quality, but have not yet considered the influences of viruses. In this study, we address this gap by mining viral DNA sequences from publicly available drinking water metagenomes from distribution systems in six countries around the world. These datasets provide a snapshot of the taxonomic diversity and metabolic potential of the global drinking water virome; and provide an opportunity to investigate the effects of geography, climate, and drinking water treatment practices on viral diversity. Both environmental conditions and differences in sample processing were found to influence the viral composition. Using free chlorine as the residual disinfectant was associated with clear differences in viral taxonomic diversity and metabolic potential, with significantly fewer viral populations and less even viral community structures than observed in distribution systems without residual disinfectant. Additionally, drinking water viruses carry antibiotic resistance genes (ARGs), as well as genes to survive oxidative stress and nitrogen limitation. Through this study, we have demonstrated that viral communities are diverse across drinking water systems and vary with the use of residual disinfectant. Our findings offer directions for future research to develop a more robust understanding of how virus-bacteria interactions in drinking water distribution systems affect water quality.

Isolation and Characterization of a Novel Cyanophage Encoding Multiple Auxiliary Metabolic Genes

As significant drivers of cyanobacteria mortality, cyanophages have been known to regulate the population dynamics, metabolic activities, and community structure of this most important marine autotrophic picoplankton and, therefore, influence the global primary production and biogeochemical cycle in aquatic ecosystems. In the present study, a lytic Synechococcus phage, namely S-SZBM1, was isolated and identified. Cyanophage S-SZBM1 has a double-stranded DNA genome of 177,834 bp with a G+C content of 43.31% and contains a total of 218 predicted ORFs and six tRNA genes. Phylogenetic analysis and nucleotide-based intergenomic similarity suggested that cyanophage S-SZBM1 belongs to a new genus under the family Kyanoviridae. A variety of auxiliary metabolic genes (AMGs) that have been proved or speculated to relate to photosynthesis, carbon metabolism, nucleotide synthesis and metabolism, cell protection, and other cell metabolism were identified in cyanophage S-SZBM1 genome and may affect host processes during infection. In addition, 24 of 32 predicted structural proteins were identified by a high-throughput proteome analysis which were potentially involved in the assembly processes of virion. The genomic and proteomic analysis features of cyanophage S-SZBM1 offer a valuable insight into the interactions between cyanophages and their hosts during infection.

Genomic Characterization of Two Novel RCA Phages Reveals New Insights into the Diversity and Evolution of Marine Viruses

Viruses are the most abundant living entities in marine ecosystems, playing critical roles in altering the structure and function of microbial communities and driving ocean biogeochemistry. Phages that infect Roseobacter clade-affiliated (RCA) cluster strains are an important component of marine viral communities. Here, we characterize the genome sequences of two new RCA phages, CRP-9 and CRP-13, which infect RCA strain FZCC0023. Genomic analysis reveals that CRP-9 and CRP-13 represent a novel evolutionary lineage of marine phages. They both have a DNA replication module most similar to those in Cobavirus group phages. In contrast, their morphogenesis and packaging modules are distinct from those in cobaviruses but homologous to those in HMO-2011-type phages. The genomic architecture of CRP-9 and CRP-13 suggests a genomic recombination event between distinct phage groups. Metagenomic data sets were examined for metagenome-assembled viral genomes (MAVGs) with similar recombinant genome architectures. Fifteen CRP-9-type MAVGs were identified from marine viromes. Additionally, 158 MAVGs were identified containing HMO-2011-type morphogenesis and packaging modules with other types of DNA replication genes, providing more evidence that recombination between different phage groups is a major driver of phage evolution. Altogether, this study significantly expands the understanding of diversity and evolution of marine roseophages. Meanwhile, the analysis of these novel RCA phages and MAVGs highlights the critical role of recombination in shaping phage diversity. These phage sequences are valuable resources for inferring the evolutionary connection of distinct phage groups. IMPORTANCE Diversity and evolution of phages that infect the relatively slow-growing but dominant Roseobacter lineages are largely unknown. In this study, RCA phages CRP-9 and CRP-13 have been isolated on a Roseobacter RCA strain and shown to have a unique genomic architecture, which appears to be the result of a recombination event. CRP-9 and CRP-13 have a DNA replication module most similar to those in Cobavirus group phages and morphogenesis and packaging modules most similar to those in HMO-2011-type phages. HMO-2011-type morphogenesis and packaging modules are found in combination with distinct types of DNA replication genes, suggesting compatibility with various DNA replication modules. Altogether, this study contributes toward a better understanding of marine viral diversity and evolution.

Prophage Genomics and Ecology in the Family Rhodobacteraceae

Roseobacters are globally abundant bacteria with critical roles in carbon and sulfur biogeochemical cycling. Here, we identified 173 new putative prophages in 79 genomes of Rhodobacteraceae. These prophages represented 1.3 ± 0.15% of the bacterial genomes and had no to low homology with reference and metagenome-assembled viral genomes from aquatic and terrestrial ecosystems. Among the newly identified putative prophages, 35% encoded auxiliary metabolic genes (AMGs), mostly involved in secondary metabolism, amino acid metabolism, and cofactor and vitamin production. The analysis of integration sites and gene homology showed that 22 of the putative prophages were actually gene transfer agents (GTAs) similar to a GTA of Rhodobacter capsulatus. Twenty-three percent of the predicted prophages were observed in the TARA Oceans viromes generated from free viral particles, suggesting that they represent active prophages capable of induction. The distribution of these prophages was significantly associated with latitude and temperature. The prophages most abundant at high latitudes encoded acpP, an auxiliary metabolic gene involved in lipid synthesis and membrane fluidity at low temperatures. Our results show that prophages and gene transfer agents are significant sources of genomic diversity in roseobacter, with potential roles in the ecology of this globally distributed bacterial group.