Environmental Viral Genomes Shed New Light on Virus-Host Interactions in the Ocean

Potential of a newly isolated lytic bacteriophage to control Pseudomonas coronafaciens pv. garcae in coffee plants: Molecular characterization with in vitro and ex vivo experiments

Traditionally, control of coffee plant bacterial halo blight (BHB) caused by the phytopathogen Pseudomonas coronafaciens pv. garcae (Pcg) involves frequent spraying of coffee plantations with non-environmentally friendly and potentially bacterial resistance-promoting copper products or with kasugamycin hydrochloride. In this study we report a leap forward in the quest for a new ecofriendly approach, characterizing (both physicochemically and biologically) and testing both in vitro and ex vivo a new lytic phage for Pcg. An in-depth molecular (genomic and DNA structural features) characterization of the phage was also undertaken. Phage PcgS01F belongs to the class Caudoviricetes, Drexlerviridae family and genus Guelphvirus, and presents a siphovirus-like morphotype. Phage PcgS01F showed a latency period of 40 min and a burst size of 46 PFU/host cell, allowing to conclude that it replicates well in Pcg IBSBF-158. At Multiplicity Of Infection (MOI, or the ratio of phage to bacteria) 1000, the performance of phage PcgS01F was much better than at MOI 10, promoting increasing bacterial reductions until the end of the in vitro inactivation assays, stabilizing at a significant 82 % bacterial load reduction. Phage PcgS01F infected and killed Pcg cells ex vivo in coffee plant leaves artificially contaminated, with a maximum of Pcg inactivation of 7.66 log CFU/mL at MOI 1000 after 36 h of incubation. This study provides evidence that the isolated phage is a promising candidate against the causative agent of BHB in coffee plants.

Isolation and characterization of new lytic bacteriophage PSA-KC1 against Pseudomonas aeruginosa isolates from cystic fibrosis patients

A novel lytic bacteriophage, PSA-KC1, was isolated from wastewater. In this study, the whole genome of the bacteriophage PSA-KC1 was analyzed, and its lytic properties were assessed. PSA-KC1 has a linear double-stranded DNA genome with a total length of 43,237 base pairs and a GC content of 53.6%. In total, 65 genes were predicted, 46 of which were assigned functions as structural proteins involved in genome replication, packaging or phage lysis. PSA-KC1 belongs to the genus Septimatrevirus under the Caudoviricetes class. The aim of this study was to investigate the efficacy of the lytic bacteriophage PSA-KC1 and compare it with that of the Pyophage phage cocktail on 25 multi drug resistant (MDR) Pseudomonas aeruginosa strains isolated from sputum samples of cystic fibrosis patients. Seventeen of these strains were susceptible (68%) to the PSA-KC1 lytic phage we isolated, whereas eight clinical strains were resistant. However, 22 (88%) of the P. aeruginosa strains were susceptible to the Pyophage cocktail, and three (12%) were resistant to the Phage cocktail. At the end of our study, a new lytic phage active against multidrug-resistant P. aeruginosa strains from CF patients was isolated, and its genome was characterized. Since the PSA-KC1 phage does not contain virulence factors, toxins or integrase genes, it can be expected to be a therapeutic candidate with the potential to be used safely in phage therapy.

A new bacteriophage infecting Staphylococcus epidermidis with potential for removing biofilms by combination with chimeric lysin CHAPSH3b and vancomycin

Staphylococcus epidermidis is the cause of serious skin and prosthetic joint infections despite being a common inhabitant of human body surfaces. However, both the rise in antibiotic resistance in this species and its ability to form biofilms are increasingly limiting the available therapeutic options against these illnesses. In this landscape, phage therapy stands out as an interesting alternative to antibiotics. In the present study, we report the isolation and characterization of a novel virulent phage infecting S. epidermidis (Staphylococcus phage IPLA-AICAT), which belongs to the Herelleviridae family. The estimated genome size of this virus is 139,941 bp, and sequence analysis demonstrated the absence of antibiotic resistance genes and virulence factors. This phage infects a high proportion (79%) of clinically relevant S. epidermidis strains and exhibits antibiofilm activity. Moreover, a combination of this phage with other antimicrobials, i.e., vancomycin and the lytic protein CHAPSH3b, further improved the reduction in surface-attached bacteria. Notably, the combination of Staphylococcus phage IPLA-AICAT (109 PFU/mL) and CHAPSH3b (8 µM), originally designed to target Staphylococcus aureus, was able to reduce the number of viable cells by 3.06 log units in 5-h-old biofilms. In 24-h-old biofilms, the reduction was also significant after 6 h of treatment (2.06 log units) and 24 h of treatment (2.52 log units). These results confirm our previous data regarding the potential of phage-lysin mixtures against staphylococcal biofilms.IMPORTANCEStaphylococcus epidermidis is one of the main causes to device-associated infections mostly due to its ability to form stable biofilms attached to human tissues. Besides the inherent antimicrobial tolerance of biofilms, this microorganism is also increasingly becoming resistant to standard-of-care antibiotics. To fight against this problem, phage therapy is a viable option to complement the available antibiotics in the treatment of recalcitrant infections. This work describes a new phage infecting S. epidermidis clinical strains that is a member of the Herelleviridae family and the combination with other antimicrobials to further improve the reduction of biofilms. Together with the significant progress achieved in the development of diagnostic tools, phages and their derived proteins will bring us much closer to a therapeutic landscape in which we are not so heavily reliant on antibiotics to combat bacterial pathogens.

The influence of depth on the global deep-sea plasmidome

Plasmids play a crucial role in facilitating genetic exchange and enhancing the adaptability of microbial communities. Despite their importance, environmental plasmids remain understudied, particularly those in fragile and underexplored ecosystems such as the deep-sea. In this paper we implemented a bioinformatics pipeline to study the composition, diversity, and functional attributes of plasmid communities (plasmidome) in 81 deep-sea metagenomes from the Tara and Malaspina expeditions, sampled from the Pacific, Atlantic, and Indian Oceans at depths ranging from 270 to 4005 m. We observed an association between depth and plasmid traits, with the 270-1000 m range (mesopelagic samples) exhibiting the highest number of plasmids and the largest plasmid sizes. Plasmids of Alphaproteobacteria and Gammaproteobacteria were predominant across the oceans, particularly in this depth range, which also showed the highest species diversity and abundance of metabolic pathways, including aromatic compound degradation. Surprisingly, relatively few antibiotic resistance genes were found in the deep-sea ecosystem, with most being found in the mesopelagic layer. These included classes such as beta-lactamase, biocide resistance, and aminoglycosides. Our study also identified the MOBP and MOBQ relaxase families as prevalent across various taxonomic classes. This research underscores the importance of studying the plasmidome independently from the chromosomal context. Our limited understanding of the deep-sea's microbial ecology, especially its plasmidome, necessitates caution in human activities like mining. Such activities could have unforeseen impacts on this largely unexplored ecosystem.

Taxonomic difference in marine bloom-forming phytoplanktonic species affects the dynamics of both bloom-responding prokaryotes and prokaryotic viruses

The production of dissolved organic matter during phytoplankton blooms and consumption by heterotrophic prokaryotes promote marine carbon biogeochemical cycling. Although prokaryotic viruses presumably affect this process, their dynamics during blooms are not fully understood. Here, we investigated the effects of taxonomic difference in bloom-forming phytoplankton on prokaryotes and their viruses. We analyzed the dynamics of coastal prokaryotic communities and viruses under the addition of dissolved intracellular fractions from taxonomically distinct phytoplankton, the diatom Chaetoceros sp. (CIF) and the raphidophycean alga Heterosigma akashiwo (HIF), using microcosm experiments. Ribosomal RNA gene amplicon and viral metagenomic analyses revealed that particular prokaryotes and prokaryotic viruses specifically increased in either CIF or HIF, indicating that taxonomic difference in bloom-forming phytoplankton promotes distinct dynamics of not only the prokaryotic community but also prokaryotic viruses. Furthermore, combining our microcosm experiments with publicly available environmental data mining, we identified both known and novel possible host-virus pairs. In particular, the growth of prokaryotes associating with phytoplanktonic organic matter, such as Bacteroidetes (Polaribacter and NS9 marine group), Vibrio spp., and Rhodobacteriales (Nereida and Planktomarina), was accompanied by an increase in viruses predicted to infect Bacteroidetes, Vibrio, and Rhodobacteriales, respectively. Collectively, our findings suggest that changes in bloom-forming species can be followed by an increase in a specific group of prokaryotes and their viruses and that elucidating these tripartite relationships among specific phytoplankton, prokaryotes, and prokaryotic viruses improves our understanding of coastal biogeochemical cycling in blooms.IMPORTANCEThe primary production during marine phytoplankton bloom and the consumption of the produced organic matter by heterotrophic prokaryotes significantly contribute to coastal biogeochemical cycles. While the activities of those heterotrophic prokaryotes are presumably affected by viral infection, the dynamics of their viruses during blooms are not fully understood. In this study, we experimentally demonstrated that intracellular fractions of taxonomically distinct bloom-forming phytoplankton species, the diatom Chaetoceros sp. and the raphidophycean alga Heterosigma akashiwo, promoted the growth of taxonomically different prokaryotes and prokaryotic viruses. Based on their dynamics and predicted hosts of those viruses, we succeeded in detecting already-known and novel possible host-virus pairs associating with either phytoplankton species. Altogether, we propose that the succession of bloom-forming phytoplankton would change the composition of the abundant prokaryotes, resulting in an increase in their viruses. These changes in viral composition, depending on bloom-forming species, would alter the dynamics and metabolism of prokaryotes, affecting biogeochemical cycling in blooms.

Isolation of viruses, including mollivirus, with the potential to infect Acanthamoeba from a Japanese warm temperate zone

Acanthamoeba castellanii is infected with diverse nucleocytoplasmic large DNA viruses. Here, we report the co-isolation of 12 viral strains from marine sediments in Uranouchi Inlet, Kochi, Japan. Based on the morphological features revealed by electron microscopy, these isolates were classified into four viral groups including Megamimiviridae, Molliviridae, Pandoraviridae, and Pithoviridae. Genomic analyses indicated that these isolates showed high similarities to the known viral genomes with which they are taxonomically clustered, and their phylogenetic relationships were also supported by core gene similarities. It is noteworthy that Molliviridae was isolated from the marine sediments in the Japanese warm temperate zone because other strains have only been found in the subarctic region. Furthermore, this strain has 19 and 4 strain-specific genes found in Mollivirus sibericum and Mollivirus kamchatka, respectively. This study extends our knowledge about the habitat and genomic diversity of Molliviridae.

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

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

Dinucleotide biases in the genomes of prokaryotic and eukaryotic dsDNA viruses and their hosts

The genomes of cellular organisms display CpG and TpA dinucleotide composition biases. Such biases have been poorly investigated in dsDNA viruses. Here, we show that in dsDNA virus, bacterial, and eukaryotic genomes, the representation of TpA and CpG dinucleotides is strongly dependent on genomic G + C content. Thus, the classical observed/expected ratios do not fully capture dinucleotide biases across genomes. Because a larger portion of the variance in TpA frequency was explained by G + C content, we explored which additional factors drive the distribution of CpG dinucleotides. Using the residuals of the linear regressions as a measure of dinucleotide abundance and ancestral state reconstruction across eukaryotic and prokaryotic virus trees, we identified an important role for phylogeny in driving CpG representation. Nonetheless, phylogenetic ANOVA analyses showed that few host associations also account for significant variations. Among eukaryotic viruses, most significant differences were observed between arthropod-infecting viruses and viruses that infect vertebrates or unicellular organisms. However, an effect of viral DNA methylation status (either driven by the host or by viral-encoded methyltransferases) is also likely. Among prokaryotic viruses, cyanobacteria-infecting phages resulted to be significantly CpG-depleted, whereas phages that infect bacteria in the genera Burkolderia and Staphylococcus were CpG-rich. Comparison with bacterial genomes indicated that this effect is largely driven by the general tendency for phages to resemble the host's genomic CpG content. Notably, such tendency is stronger for temperate than for lytic phages. Our data shed light into the processes that shape virus genome composition and inform manipulation strategies for biotechnological applications.

Virus-prokaryote infection pairs associated with prokaryotic production in a freshwater lake

Viruses infect and kill prokaryotic populations in a density- or frequency-dependent manner and affect carbon cycling. However, the effects of the stratification transition, including the stratified and de-stratified periods, on the changes in prokaryotic and viral communities and their interactions remain unclear. We conducted a monthly survey of the surface and deep layers of a large and deep freshwater lake (Lake Biwa, Japan) for a year and analyzed the prokaryotic production and prokaryotic and viral community composition. Our analysis revealed that, in the surface layer, 19 prokaryotic species, accounting for approximately 40% of the total prokaryotic abundance, could potentially contribute to the majority of prokaryotic production, which is the highest during the summer and is suppressed by viruses. This suggests that a small fraction of prokaryotes and phages were the key infection pairs during the peak period of prokaryotic activity in the freshwater lake. We also found that approximately 50% of the dominant prokaryotic and viral species in the deep layer were present throughout the study period. This suggests that the "kill the winner" model could explain the viral impact on prokaryotes in the surface layer, but other dynamics may be at play in the deep layer. Furthermore, we found that annual vertical mixing could result in a similar rate of community change between the surface and deep layers. These findings may be valuable in understanding how communities and the interaction among them change when freshwater lake stratification is affected by global warming in the future.IMPORTANCEViral infection associated with prokaryotic production occurs in a density- or frequency-dependent manner and regulates the prokaryotic community. Stratification transition and annual vertical mixing in freshwater lakes are known to affect the prokaryotic community and the interaction between prokaryotes and viruses. By pairing measurements of virome analysis and prokaryotic production of a 1-year survey of the depths of surface and deep layers, we revealed (i) the prokaryotic infection pairs associated with prokaryotic production and (ii) the reset in prokaryotic and viral communities through annual vertical mixing in a freshwater lake. Our results provide a basis for future work into changes in stratification that may impact the biogeochemical cycling in freshwater lakes.

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

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

Molecular Characterization and Genome Mechanical Features of Two Newly Isolated Polyvalent Bacteriophages Infecting Pseudomonas syringae pv. garcae

Coffee plants have been targeted by a devastating bacterial disease, a condition known as bacterial blight, caused by the phytopathogen Pseudomonas syringae pv. garcae (Psg). Conventional treatments of coffee plantations affected by the disease involve frequent spraying with copper- and kasugamycin-derived compounds, but they are both highly toxic to the environment and stimulate the appearance of bacterial resistance. Herein, we report the molecular characterization and mechanical features of the genome of two newly isolated (putative polyvalent) lytic phages for Psg. The isolated phages belong to class Caudoviricetes and present a myovirus-like morphotype belonging to the genuses Tequatrovirus (PsgM02F) and Phapecoctavirus (PsgM04F) of the subfamilies Straboviridae (PsgM02F) and Stephanstirmvirinae (PsgM04F), according to recent bacterial viruses' taxonomy, based on their complete genome sequences. The 165,282 bp (PsgM02F) and 151,205 bp (PsgM04F) genomes do not feature any lysogenic-related (integrase) genes and, hence, can safely be assumed to follow a lytic lifestyle. While phage PsgM02F produced a morphogenesis yield of 124 virions per host cell, phage PsgM04F produced only 12 virions per host cell, indicating that they replicate well in Psg with a 50 min latency period. Genome mechanical analyses established a relationship between genome bendability and virion morphogenesis yield within infected host cells.

DiGAlign: Versatile and Interactive Visualization of Sequence Alignment for Comparative Genomics

With the explosion of available genomic information, comparative genomics has become a central approach to understanding microbial ecology and evolution. We developed DiGAlign (https://www.genome.jp/digalign/), a web server that provides versatile functionality for comparative genomics with an intuitive interface. It allows the user to perform the highly customizable visualization of a synteny map by simply uploading nucleotide sequences of interest, ranging from a specific region to the whole genome landscape of microorganisms and viruses. DiGAlign will serve a wide range of biological researchers, particularly experimental biologists, with multifaceted features that allow the rapid characterization of genomic sequences of interest and the generation of a publication-ready figure.

Exploring the Archaeal Virosphere by Metagenomics

During the past decade, environmental research has demonstrated that archaea are abundant and widespread in nature and play important ecological roles at a global scale. Currently, however, the majority of archaeal lineages cannot be cultivated under laboratory conditions and are known exclusively or nearly exclusively through metagenomics. A similar trend extends to the archaeal virosphere, where isolated representatives are available for a handful of model archaeal virus-host systems. Viral metagenomics provides an alternative way to circumvent the limitations of culture-based virus discovery and offers insight into the diversity, distribution, and environmental impact of uncultured archaeal viruses. Presently, metagenomics approaches have been successfully applied to explore the viromes associated with various lineages of extremophilic and mesophilic archaea, including Asgard archaea (Asgardarchaeota), ANME-1 archaea (Methanophagales), thaumarchaea (Nitrososphaeria), altiarchaea (Altiarchaeota), and marine group II archaea (Poseidoniales). Here, we provide an overview of methods widely used in archaeal virus metagenomics, covering metavirome preparation, genome annotation, phylogenetic and phylogenomic analyses, and archaeal host assignment. We hope that this summary will contribute to further exploration and characterization of the enigmatic archaeal virome lurking in diverse environments.

Genomic and transcriptomic insights into complex virus-prokaryote interactions in marine biofilms

Marine biofilms are complex communities of microorganisms that play a crucial ecological role in oceans. Although prokaryotes are the dominant members of these biofilms, little is known about their interactions with viruses. By analysing publicly available and newly sequenced metagenomic data, we identified 2446 virus-prokaryote connections in 84 marine biofilms. Most of these connections were between the bacteriophages in the Uroviricota phylum and the bacteria of Proteobacteria, Cyanobacteria and Bacteroidota. The network of virus-host pairs is complex; a single virus can infect multiple prokaryotic populations or a single prokaryote is susceptible to several viral populations. Analysis of genomes of paired prokaryotes and viruses revealed the presence of 425 putative auxiliary metabolic genes (AMGs), 239 viral genes related to restriction-modification (RM) systems and 38,538 prokaryotic anti-viral defence-related genes involved in 15 defence systems. Transcriptomic evidence from newly established biofilms revealed the expression of viral genes, including AMGs and RM, and prokaryotic defence systems, indicating the active interplay between viruses and prokaryotes. A comparison between biofilms and seawater showed that biofilm prokaryotes have more abundant defence genes than seawater prokaryotes, and the defence gene composition differs between biofilms and the surrounding seawater. Overall, our study unveiled active viruses in natural biofilms and their complex interplay with prokaryotes, which may result in the blooming of defence strategists in biofilms. The detachment of bloomed defence strategists may reduce the infectivity of viruses in seawater and result in the emergence of a novel role of marine biofilms.

ICTV Virus Taxonomy Profile: Aoguangviridae 2023

The family Aoguangviridae includes dsDNA viruses that have been associated with marine archaea. Currently, members of this virus family are known through metagenomics. Virions are predicted to consist of an icosahedral capsid and a helical tail, characteristic of members in the class Caudoviricetes. Aoguangviruses have some of the largest genomes among archaeal viruses and possess most of the components of the DNA replication machinery as well as auxiliary functions. The family Aoguangviridae includes the species Aobingvirus yangshanense. Many unclassified relatives of this virus group, referred to as 'magroviruses', have been discovered by metagenomics in globally distributed marine samples. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Aoguangviridae, which is available at ictv.global/report/aoguangviridae.

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.

Emergence of Raoultella ornithinolytica in human infections from different hospitals in Ecuador with OXA-48-producing resistance

Purpose

The purpose of this study was to highlight the clinical and molecular features of 13 Raoultella ornithinolytica strains isolated from clinical environments in Ecuador, and to perform comparative genomics with previously published genomes of Raoultella spp. As Raoultella is primarily found in environmental, clinical settings, we focused our work on identifying mechanisms of resistance that can provide this bacterium an advantage to establish and persist in hospital environments.

Methods

We analyzed 13 strains of Raoultella ornithinolytica isolated from patients with healthcare associated infections (HAI) in three hospitals in Quito and one in Santo Domingo de Los Tsáchilas, Ecuador, between November 2017 and April 2018. These isolates were subjected to phenotypic antimicrobial susceptibility testing, end-point polymerase chain reaction (PCR) to detect the presence of carbapenemases and whole-genome sequencing.

Results

Polymerase chain reaction revealed that seven isolates were positive isolates for blaOXA-48 and one for blaKPC-2 gene. Of the seven strains that presented the blaOXA-48 gene, six harbored it on an IncFII plasmid, one was inserted into the bacterial chromosome. The blaKPC gene was detected in an IncM2/IncR plasmid. From the bioinformatics analysis, nine genomes had the gene blaOXA-48, originating from Ecuador. Moreover, all R. ornithinolytica strains contained the ORN-1 gene, which confers resistance for β-lactams, such as penicillins and cephalosporins. Comparative genome analysis of the strains showed that the pangenome of R. ornithinolytica is considered an open pangenome, with 27.77% of core genes, which could be explained by the fact that the antibiotic resistance genes in the ancestral reconstruction are relatively new, suggesting that this genome is constantly incorporating new genes.

Conclusion

These results reveal the genome plasticity of R. ornithinolytica, particularly in acquiring antibiotic-resistance genes. The genomic surveillance and infectious control of these uncommon species are important since they may contribute to the burden of antimicrobial resistance and human health.

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.

Prevalence of Viral Frequency-Dependent Infection in Coastal Marine Prokaryotes Revealed Using Monthly Time Series Virome Analysis

Viruses infecting marine prokaryotes have a large impact on the diversity and dynamics of their hosts. Model systems suggest that viral infection is frequency dependent and constrained by the virus-host encounter rate. However, it is unclear whether frequency-dependent infection is pervasive among the abundant prokaryotic populations with different temporal dynamics. To address this question, we performed a comparison of prokaryotic and viral communities using 16S rRNA amplicon and virome sequencing based on samples collected monthly for 2 years at a Japanese coastal site, Osaka Bay. Concurrent seasonal shifts observed in prokaryotic and viral community dynamics indicated that the abundance of viruses correlated with that of their predicted host phyla (or classes). Cooccurrence network analysis between abundant prokaryotes and viruses revealed 6,423 cooccurring pairs, suggesting a tight coupling of host and viral abundances and their "one-to-many" correspondence. Although stable dominant species, such as SAR11, showed few cooccurring viruses, a fast succession of their viruses suggests that viruses infecting these populations changed continuously. Our results suggest that frequency-dependent viral infection prevails in coastal marine prokaryotes regardless of host taxa and temporal dynamics. IMPORTANCE There is little room for doubt that viral infection is prevalent among abundant marine prokaryotes regardless of their taxa or growth strategy. However, comprehensive evaluations of viral infections in natural prokaryotic communities are still technically difficult. In this study, we examined viral infection in abundant prokaryotes by monitoring the monthly dynamics of prokaryotic and viral communities at a eutrophic coastal site, Osaka Bay. We compared the community dynamics of viruses with those of their putative hosts based on genome-based in silico host prediction. We observed frequent cooccurrence among the predicted virus-host pairs, suggesting that viral infection is prevalent in abundant prokaryotes regardless of their taxa or temporal dynamics. This likely indicates that frequent lysis of the abundant prokaryotes via viral infection has a considerable contribution to the biogeochemical cycling and maintenance of prokaryotic community diversity.

Diverse viruses of marine archaea discovered using metagenomics

During the past decade, metagenomics became a method of choice for the discovery of novel viruses. However, host assignment for uncultured viruses remains challenging, especially for archaeal viruses, which are grossly undersampled compared to viruses of bacteria and eukaryotes. Here, we assessed the utility of CRISPR spacer targeting, tRNA gene matching and homology searches for viral signature proteins, such as major capsid proteins, for the assignment of archaeal hosts and validated these approaches on metaviromes from Yangshan Harbor (YSH). We report 35 new genomes of viruses which could be confidently assigned to hosts representing diverse lineages of marine archaea. We show that the archaeal YSH virome is highly diverse, with some viruses enriching the previously described virus groups, such as magroviruses of Marine Group II Archaea (Poseidoniales), and others representing novel groups of marine archaeal viruses. Metagenomic recruitment of Tara Oceans datasets on the YSH viral genomes demonstrated the presence of YSH Poseidoniales and Nitrososphaeria viruses in the global oceans, but also revealed the endemic YSH-specific viral lineages. Furthermore, our results highlight the relationship between the soil and marine thaumarchaeal viruses. We propose three new families within the class Caudoviricetes for the classification of the five complete viral genomes predicted to replicate in marine Poseidoniales and Nitrososphaeria, two ecologically important and widespread archaeal groups. This study illustrates the utility of viral metagenomics in exploring the archaeal virome and provides new insights into the diversity, distribution and evolution of marine archaeal viruses.

PhaTYP: predicting the lifestyle for bacteriophages using BERT

Bacteriophages (or phages), which infect bacteria, have two distinct lifestyles: virulent and temperate. Predicting the lifestyle of phages helps decipher their interactions with their bacterial hosts, aiding phages' applications in fields such as phage therapy. Because experimental methods for annotating the lifestyle of phages cannot keep pace with the fast accumulation of sequenced phages, computational method for predicting phages' lifestyles has become an attractive alternative. Despite some promising results, computational lifestyle prediction remains difficult because of the limited known annotations and the sheer amount of sequenced phage contigs assembled from metagenomic data. In particular, most of the existing tools cannot precisely predict phages' lifestyles for short contigs. In this work, we develop PhaTYP (Phage TYPe prediction tool) to improve the accuracy of lifestyle prediction on short contigs. We design two different training tasks, self-supervised and fine-tuning tasks, to overcome lifestyle prediction difficulties. We rigorously tested and compared PhaTYP with four state-of-the-art methods: DeePhage, PHACTS, PhagePred and BACPHLIP. The experimental results show that PhaTYP outperforms all these methods and achieves more stable performance on short contigs. In addition, we demonstrated the utility of PhaTYP for analyzing the phage lifestyle on human neonates' gut data. This application shows that PhaTYP is a useful means for studying phages in metagenomic data and helps extend our understanding of microbial communities.

Diversity, distribution, and functional potentials of magroviruses from marine and brackish waters

Marine group II (MGII) archaea (Ca. Poseidoniales) are among the most abundant microbes in global oceanic surface waters and play an important role in driving marine biogeochemical cycles. Magroviruses - the viruses of MGII archaea have been recently found to occur ubiquitously in surface ocean. However, their diversity, distribution, and potential ecological functions in coastal zones especially brackish waters are unknown. Here we obtained 234 non-redundant magroviral genomes from brackish surface waters by using homology searches for viral signature proteins highlighting the uncovered vast diversity of this novel viral group. Phylogenetic analysis based on these brackish magroviruses along with previously reported marine ones identified six taxonomic groups with close evolutionary connection to both haloviruses and the viruses of Marine Group I archaea. Magroviruses were present abundantly both in brackish and open ocean samples with some showing habitat specification and others having broad spectrums of distribution between different habitats. Genome annotation suggests they may be involved in regulating multiple metabolic pathways of MGII archaea. Our results uncover the previously overlooked diversity and ecological potentials of a major archaeal virial group in global ocean and brackish waters and shed light on the cryptic evolutionary history of archaeal viruses.

Significant Differences in Planktonic Virus Communities Between "Cellular Fraction" (0.22 ~ 3.0 µm) and "Viral Fraction" (< 0.22 μm) in the Ocean

Compared to free-living viruses (< 0.22 m) in the ocean, planktonic viruses in the "cellular fraction" (0.22 ~ 3.0 μm) are now far less well understood, and the differences between them remain largely unexplored. Here, we revealed that even in the same seawater samples, the "cellular fraction" comprised significantly distinct virus communities from the free virioplankton, with only 13.87% overlap in viral contigs at the species level. Compared to the viral genomes deposited in NCBI RefSeq database, 99% of the assembled viral genomes in the "cellular fraction" represented novel genera. Notably, the assembled (near-) complete viral genomes within the "cellular fraction" were significantly larger than that in the "viral fraction," and the "cellular fraction" contained three times more species of giant viruses or jumbo phages with genomes > 200 kb than the "viral fraction." The longest complete genomes of jumbo phage (~ 252 kb) and giant virus (~ 716 kb) were both detected only in the "cellular fraction." Moreover, a relatively higher proportion of proviruses were predicted within the "cellular fraction" than "viral fraction." Besides the substantial divergence in viral community structure, the different fractions also contained their unique viral auxiliary metabolic genes; e.g., those potentially participating in inorganic carbon fixation in deep sea were detected only in the "cellular-fraction" viromes. In addition, there was a considerable divergence in the community structure of both "cellular fraction" and "viral fraction" viromes between the surface and deep-sea habitats, suggesting that they might have similar environmental adaptation properties. The findings deepen our understanding of the complexity of viral community structure and function in the ocean.

Novel Virus Identification through Metagenomics: A Systematic Review

Metagenomic Next Generation Sequencing (mNGS) allows the evaluation of complex microbial communities, avoiding isolation and cultivation of each microbial species, and does not require prior knowledge of the microbial sequences present in the sample. Applications of mNGS include virome characterization, new virus discovery and full-length viral genome reconstruction, either from virus preparations enriched in culture or directly from clinical and environmental specimens. Here, we systematically reviewed studies that describe novel virus identification through mNGS from samples of different origin (plant, animal and environment). Without imposing time limits to the search, 379 publications were identified that met the search parameters. Sample types, geographical origin, enrichment and nucleic acid extraction methods, sequencing platforms, bioinformatic analytical steps and identified viral families were described. The review highlights mNGS as a feasible method for novel virus discovery from samples of different origins, describes which kind of heterogeneous experimental and analytical protocols are currently used and provides useful information such as the different commercial kits used for the purification of nucleic acids and bioinformatics analytical pipelines.

Widespread yet persistent low abundance of TIM5-like cyanophages in the oceans

Ocean ecosystems are inhabited by a diverse set of viruses that impact microbial mortality and evolution. However, the distribution and abundances of specific viral lineages, particularly those from the large bank of rare viruses, remains largely unknown. Here, we assessed the diversity and abundance of the TIM5-like cyanophages. The sequencing of three new TIM5-like cyanophage genomes and environmental amplicons of a signature gene from the Red Sea revealed highly conserved gene content and sequence similarity. We adapted the polony method, a solid-phase polymerase chain reaction assay, to quantify TIM5-like cyanophages during three 2000 km expeditions in the Pacific Ocean and four annual cycles in the Red Sea. TIM5-like cyanophages were widespread, detected at all latitudes and seasons surveyed throughout the photic zone. Yet they were generally rare, ranging between <100 and 4000 viruses·ml^-1 . Occasional peaks in abundance of 10- to 100-fold were observed, reaching 71,000 viruses·ml^-1 . These peaks were ephemeral and seasonally variable in the Red Sea. Infection levels, quantified during one such peak, were very low. These characteristics of low diversity and abundance, as well as variable outbreaks, distinguishes the TIM5-like lineage from other major cyanophage lineages and illuminates that rare virus lineages can be persistent and widespread in the oceans.

Increased phage resistance through lysogenic conversion accompanying emergence of monophasic Salmonella Typhimurium ST34 pandemic strain

Salmonella enterica serovar Typhimurium (S. Typhimurium) comprises a group of closely related human and animal pathogens that account for a large proportion of all Salmonella infections globally. The epidemiological record of S. Typhimurium in Europe is characterized by successive waves of dominant clones, each prevailing for approximately 10-15 years before replacement. Succession of epidemic clones may represent a moving target for interventions aimed at controlling the spread and impact of this pathogen on human and animal health. Here, we investigate the relationship of phage sensitivity and population structure of S. Typhimurium using data from the Anderson phage typing scheme. We observed greater resistance to phage predation of epidemic clones circulating in livestock over the past decades compared to variants with a restricted host range implicating increased resistance to phage in the emergence of epidemic clones of particular importance to human health. Emergence of monophasic S. Typhimurium ST34, the most recent dominant multidrug-resistant clone, was accompanied by increased resistance to phage predation during clonal expansion, in part by the acquisition of the mTmII prophage that may have contributed to the fitness of the strains that replaced ancestors lacking this prophage.

Isolation and Molecular Characterization of Two Novel Lytic Bacteriophages for the Biocontrol of Escherichia coli in Uterine Infections: In Vitro and Ex Vivo Preliminary Studies in Veterinary Medicine

E. coli is one of the etiological agents responsible for pyometra in female dogs, with conventional treatment involving ovariohysterectomy. Here, we report the isolation and full characterization of two novel lytic phages, viz. vB_EcoM_Uniso11 (ph0011) and vB_EcoM_Uniso21 (ph0021). Both phages belong to the order Caudovirales and present myovirus-like morphotypes, with phage ph0011 being classified as Myoviridae genus Asteriusvirus and phage ph0021 being classified as Myoviridae genus Tequatrovirus, based on their complete genome sequences. The 348,288 bp phage ph0011 and 165,222 bp phage ph0021 genomes do not encode toxins, integrases or antimicrobial resistance genes neither depolymerases related sequences. Both phages were shown to be effective against at least twelve E. coli clinical isolates in in vitro antibacterial activity assays. Based on their features, both phages have potential for controlling pyometra infections caused by E. coli. Phage ph0011 (reduction of 4.24 log CFU/mL) was more effective than phage ph0021 (reduction of 1.90 log CFU/mL) after 12 h of incubation at MOI 1000. As a cocktail, the two phages were highly effective in reducing the bacterial load (reduction of 5.57 log CFU/mL) at MOI 100, after 12 h of treatment. Both phages were structurally and functionally stabilized in vaginal egg formulations.

Viral spillover risk increases with climate change in High Arctic lake sediments

The host spectrum of viruses is quite diverse, as they can sustainedly infect a few species to several phyla. When confronted with a new host, a virus may even infect it and transmit sustainably in this new host, a process called ‘viral spillover’. However, the risk of such events is difficult to quantify. As climate change is rapidly transforming environments, it is becoming critical to quantify the potential for spillovers. To address this issue, we resorted to a metagenomics approach and focused on two environments, soil and lake sediments from Lake Hazen, the largest High Arctic freshwater lake in the world. We used DNA and RNA sequencing to reconstruct the lake’s virosphere in both its sediments and soils, as well as its range of eukaryotic hosts. We then estimated the spillover risk by measuring the congruence between the viral and the eukaryotic host phylogenetic trees, and show that spillover risk increases with runoff from glacier melt, a proxy for climate change. Should climate change also shift species range of potential viral vectors and reservoirs northwards, the High Arctic could become fertile ground for emerging pandemics.

Ecogenomics reveals viral communities across the Challenger Deep oceanic trench

Despite the environmental challenges and nutrient scarcity, the geographically isolated Challenger Deep in Mariana trench, is considered a dynamic hotspot of microbial activity. Hadal viruses are the least explored microorganisms in Challenger Deep, while their taxonomic and functional diversity and ecological impact on deep-sea biogeochemistry are poorly described. Here, we collect 13 sediment cores from slope and bottom-axis sites across the Challenger Deep (down to ~11 kilometers depth), and identify 1,628 previously undescribed viral operational taxonomic units at species level. Community-wide analyses reveals 1,299 viral genera and distinct viral diversity across the trench, which is significantly higher at the bottom-axis vs. slope sites of the trench. 77% of these viral genera have not been previously identified in soils, deep-sea sediments and other oceanic settings. Key prokaryotes involved in hadal carbon and nitrogen cycling are predicted to be potential hosts infected by these viruses. The detected putative auxiliary metabolic genes suggest that viruses at Challenger Deep could modulate the carbohydrate and sulfur metabolisms of their potential hosts, and stabilize host's cell membranes under extreme hydrostatic pressures. Our results shed light on hadal viral metabolic capabilities, contribute to understanding deep sea ecology and on functional adaptions of hadal viruses for future research.

Viruses in Subsurface Environments

Over the past 20 years, our knowledge of virus diversity and abundance in subsurface environments has expanded dramatically through application of quantitative metagenomic approaches. In most subsurface environments, viral diversity and abundance rival viral diversity and abundance observed in surface environments. Most of these viruses are uncharacterized in terms of their hosts and replication cycles. Analysis of accessory metabolic genes encoded by subsurface viruses indicates that they evolved to replicate within the unique features of their environments. The key question remains: What role do these viruses play in the ecology and evolution of the environments in which they replicate? Undoubtedly, as more virologists examine the role of viruses in subsurface environments, new insights will emerge.

Extensive gut virome variation and its associations with host and environmental factors in a population-level cohort

Indigenous bacteriophage communities (virome) in the human gut have a huge impact on the structure and function of gut bacterial communities (bacteriome), but virome variation at a population scale is not fully investigated yet. Here, we analyse the gut dsDNA virome in the Japanese 4D cohort of 4198 deeply phenotyped individuals. By assembling metagenomic reads, we discover thousands of high-quality phage genomes including previously uncharacterised phage clades with different bacterial hosts than known major ones. The distribution of host bacteria is a strong determinant for the distribution of phages in the gut, and virome diversity is highly correlated with anti-viral defence mechanisms of the bacteriome, such as CRISPR-Cas and restriction-modification systems. We identify 97 various intrinsic/extrinsic factors that significantly affect the virome structure, including age, sex, lifestyle, and diet, most of which showed consistent associations with both phages and their predicted bacterial hosts. Among the metadata categories, disease and medication have the strongest effects on the virome structure. Overall, these results present a basis to understand the symbiotic communities of bacteria and their viruses in the human gut, which will facilitate the medical and industrial applications of indigenous viruses.

The Synthesis and Assembly of a Truncated Cyanophage Genome and Its Expression in a Heterogenous Host

Cyanophages play an important role in regulating the dynamics of cyanobacteria communities in the hydrosphere, representing a promising biological control strategy for cyanobacterial blooms. Nevertheless, most cyanophages are host-specific, making it difficult to control blooming cyanobacteria via single or multiple cyanophages. In order to address the issue, we explore the interaction between cyanophages and their heterologous hosts, with the aim of revealing the principles of designing and constructing an artificial cyanophage genome towards multiple cyanobacterial hosts. In the present study, we use synthetic biological approaches to assess the impact of introducing a fragment of cyanophage genome into a heterologous cyanobacterium under a variety of environmental conditions. Based on a natural cyanophage A-4L genome (41,750 bp), a truncated cyanophage genome Syn-A-4-8 is synthesized and assembled in Saccharomyces cerevisiae. We found that a 351-15,930 bp area of the A-4L genome has a fragment that is lethal to Escherichia coli during the process of attempting to assemble the full-length A-4L genome. Syn-A-4-8 was successfully introduced into E. coli and then transferred into the model cyanobacterium Synechococcus elongatus PCC 7942 (Syn7942) via conjugation. Although no significant phenotypes of Syn7942 carrying Syn-A-4-8 (LS-02) could be observed under normal conditions, its growth exhibited a prolonged lag phase compared to that of the control strain under 290-millimolar NaCl stress. Finally, the mechanisms of altered salt tolerance in LS-02 were revealed through comparative transcriptomics, and ORF25 and ORF26 on Syn-A-4-8 turned out to be the key genes causing the phenotype. Our research represents an important attempt in designing artificial cyanophages towards multiple hosts, and offers new future insights into the control of cyanobacterial blooms.

Accurate identification of bacteriophages from metagenomic data using Transformer

Motivation

Bacteriophages are viruses infecting bacteria. Being key players in microbial communities, they can regulate the composition/function of microbiome by infecting their bacterial hosts and mediating gene transfer. Recently, metagenomic sequencing, which can sequence all genetic materials from various microbiome, has become a popular means for new phage discovery. However, accurate and comprehensive detection of phages from the metagenomic data remains difficult. High diversity/abundance, and limited reference genomes pose major challenges for recruiting phage fragments from metagenomic data. Existing alignment-based or learning-based models have either low recall or precision on metagenomic data.

Results

In this work, we adopt the state-of-the-art language model, Transformer, to conduct contextual embedding for phage contigs. By constructing a protein-cluster vocabulary, we can feed both the protein composition and the proteins’ positions from each contig into the Transformer. The Transformer can learn the protein organization and associations using the self-attention mechanism and predicts the label for test contigs. We rigorously tested our developed tool named PhaMer on multiple datasets with increasing difficulty, including quality RefSeq genomes, short contigs, simulated metagenomic data, mock metagenomic data and the public IMG/VR dataset. All the experimental results show that PhaMer outperforms the state-of-the-art tools. In the real metagenomic data experiment, PhaMer improves the F1-score of phage detection by 27%.

The OceanDNA MAG catalog contains over 50,000 prokaryotic genomes originated from various marine environments

Marine microorganisms are immensely diverse and play fundamental roles in global geochemical cycling. Recent metagenome-assembled genome studies, with particular attention to large-scale projects such as Tara Oceans, have expanded the genomic repertoire of marine microorganisms. However, published marine metagenome data is still underexplored. We collected 2,057 marine metagenomes covering various marine environments and developed a new genome reconstruction pipeline. We reconstructed 52,325 qualified genomes composed of 8,466 prokaryotic species-level clusters spanning 59 phyla, including genomes from the deep-sea characterized as deeper than 1,000 m (n = 3,337), low-oxygen zones of <90 μmol O2 per kg water (n = 7,884), and polar regions (n = 7,752). Novelty evaluation using a genome taxonomy database shows that 6,256 species (73.9%) are novel and include genomes of high taxonomic novelty, such as new class candidates. These genomes collectively expanded the known phylogenetic diversity of marine prokaryotes by 34.2%, and the species representatives cover 26.5-42.0% of prokaryote-enriched metagenomes. Thoroughly leveraging accumulated metagenomic data, this genome resource, named the OceanDNA MAG catalog, illuminates uncharacterized marine microbial 'dark matter' lineages.

Complete Genome Sequence of a Phage Infecting Sphingomonadaceae

We isolated a phage infecting a member of the Sphingomonadaceae family from a freshwater lake. The phage has a DNA genome of 41,771 bp, with a GC content of 61.7%. The genome harbors 50 predicted protein-coding genes and an auxiliary metabolic gene, which encodes a protein belonging to the radical S-adenosylmethionine superfamily.

Diversity and origins of bacterial and archaeal viruses on sinking particles reaching the abyssal ocean

Sinking particles and particle-associated microbes influence global biogeochemistry through particulate matter export from the surface to the deep ocean. Despite ongoing studies of particle-associated microbes, viruses in these habitats remain largely unexplored. Whether, where, and which viruses might contribute to particle production and export remain open to investigation. In this study, we analyzed 857 virus population genomes associated with sinking particles collected over three years in sediment traps moored at 4000 m in the North Pacific Subtropical Gyre. Particle-associated viruses here were linked to cellular hosts through matches to bacterial and archaeal metagenome-assembled genome (MAG)-encoded prophages or CRISPR spacers, identifying novel viruses infecting presumptive deep-sea bacteria such as Colwellia, Moritella, and Shewanella. We also identified lytic viruses whose abundances correlated with particulate carbon flux and/or were exported from the photic to abyssal ocean, including cyanophages. Our data are consistent with some of the predicted outcomes of the viral shuttle hypothesis, and further suggest that viral lysis of both autotrophic and heterotrophic prokaryotes may play a role in carbon export. Our analyses revealed the diversity and origins of prevalent viruses found on deep-sea sinking particles and identified prospective viral groups for future investigation into processes that govern particle export in the open ocean.

Newly identified HMO-2011-type phages reveal genomic diversity and biogeographic distributions of this marine viral group

Viruses play critical roles in influencing biogeochemical cycles and adjusting host mortality, population structure, physiology, and evolution in the ocean. Marine viral communities are composed of numerous genetically distinct subfamily/genus-level viral groups. Among currently identified viral groups, the HMO-2011-type group is known to be dominant and broadly distributed. However, only four HMO-2011-type cultivated representatives that infect marine SAR116 and Roseobacter strains have been reported to date, and the genetic diversity, potential hosts, and ecology of this group remain poorly elucidated. Here, we present the genomes of seven HMO-2011-type phages that were isolated using four Roseobacter strains and one SAR11 strain, as well as additional 207 HMO-2011-type metagenomic viral genomes (MVGs) identified from various marine viromes. Phylogenomic and shared-gene analyses revealed that the HMO-2011-type group is a subfamily-level group comprising at least 10 discernible genus-level subgroups. Moreover, >2000 HMO-2011-type DNA polymerase sequences were identified, and the DNA polymerase phylogeny also revealed that the HMO-2011-type group contains diverse subgroups and is globally distributed. Metagenomic read-mapping results further showed that most HMO-2011-type phages are prevalent in global oceans and display distinct geographic distributions, with the distribution of most HMO-2011-type phages being associated with temperature. Lastly, we found that members in subgroup IX, represented by pelagiphage HTVC033P, were among the most abundant HMO-2011-type phages, which implies that SAR11 bacteria are crucial hosts for this viral group. In summary, our findings substantially expand current knowledge regarding the phylogenetic diversity, evolution, and distribution of HMO-2011-type phages, highlighting HMO-2011-type phages as major ecological agents that can infect certain key bacterial groups.

Viruses affect picocyanobacterial abundance and biogeography in the North Pacific Ocean

The photosynthetic picocyanobacteria Prochlorococcus and Synechococcus are models for dissecting how ecological niches are defined by environmental conditions, but how interactions with bacteriophages affect picocyanobacterial biogeography in open ocean biomes has rarely been assessed. We applied single-virus and single-cell infection approaches to quantify cyanophage abundance and infected picocyanobacteria in 87 surface water samples from five transects that traversed approximately 2,200 km in the North Pacific Ocean on three cruises, with a duration of 2-4 weeks, between 2015 and 2017. We detected a 550-km-wide hotspot of cyanophages and virus-infected picocyanobacteria in the transition zone between the North Pacific Subtropical and Subpolar gyres that was present in each transect. Notably, the hotspot occurred at a consistent temperature and displayed distinct cyanophage-lineage composition on all transects. On two of these transects, the levels of infection in the hotspot were estimated to be sufficient to substantially limit the geographical range of Prochlorococcus. Coincident with the detection of high levels of virally infected picocyanobacteria, we measured an increase of 10-100-fold in the Synechococcus populations in samples that are usually dominated by Prochlorococcus. We developed a multiple regression model of cyanophages, temperature and chlorophyll concentrations that inferred that the hotspot extended across the North Pacific Ocean, creating a biological boundary between gyres, with the potential to release organic matter comparable to that of the sevenfold-larger North Pacific Subtropical Gyre. Our results highlight the probable impact of viruses on large-scale phytoplankton biogeography and biogeochemistry in distinct regions of the oceans.

Novel Viral Communities Potentially Assisting in Carbon, Nitrogen, and Sulfur Metabolism in the Upper Slope Sediments of Mariana Trench

Viruses are ubiquitous in the oceans. Even in the deep sediments of the Mariana Trench, viruses have high productivity. However, little is known about their species composition and survival strategies in that environment. Here, we uncovered novel viral communities (3,206 viral scaffolds) in the upper slope sediments of the Mariana Trench via metagenomic analysis of 15 sediment samples. Most (99%) of the viral scaffolds lack known viral homologs, and ca. 59% of the high-quality viral genomes (total of 111 with completeness of >90%) represent novel genera, including some Phycodnaviridae and jumbo phages. These viruses contain various auxiliary metabolic genes (AMGs) potentially involved in organic carbon degradation, inorganic carbon fixation, denitrification, and assimilatory sulfate reduction, etc. This study provides novel insight into the almost unknown benthic viral communities in the Mariana Trench. IMPORTANCE The Mariana Trench harbors a substantial number of infective viral particles. However, very little is known about the identity, survival strategy, and potential functions of viruses in the trench sediments. Here, through metagenomic analysis, unusual benthic viral communities with high diversity and novelty were discovered. Among them, 59% of the viruses with a genome completeness of >90% represent novel genera. Various auxiliary metabolic genes carried by these viruses reflect the potential adaptive characteristics of viruses in this extreme environment and the biogeochemical cycles that they may participate in. This study gives us a deeper understanding of the peculiarities of viral communities in deep-sea/hadal sediments.

Highly diverse flavobacterial phages isolated from North Sea spring blooms

It is generally recognized that phages are a mortality factor for their bacterial hosts. This could be particularly true in spring phytoplankton blooms, which are known to be closely followed by a highly specialized bacterial community. We hypothesized that phages modulate these dense heterotrophic bacteria successions following phytoplankton blooms. In this study, we focused on Flavobacteriia, because they are main responders during these blooms and have an important role in the degradation of polysaccharides. A cultivation-based approach was used, obtaining 44 lytic flavobacterial phages (flavophages), representing twelve new species from two viral realms. Taxonomic analysis allowed us to delineate ten new phage genera and ten new families, from which nine and four, respectively, had no previously cultivated representatives. Genomic analysis predicted various life styles and genomic replication strategies. A likely eukaryote-associated host habitat was reflected in the gene content of some of the flavophages. Detection in cellular metagenomes and by direct-plating showed that part of these phages were actively replicating in the environment during the 2018 spring bloom. Furthermore, CRISPR/Cas spacers and re-isolation during two consecutive years suggested that, at least part of the new flavophages are stable components of the microbial community in the North Sea. Together, our results indicate that these diverse flavophages have the potential to modulate their respective host populations.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

Arms race in a cell: genomic, transcriptomic, and proteomic insights into intracellular phage-bacteria interplay in deep-sea snail holobionts

BACKGROUND

Deep-sea animals in hydrothermal vents often form endosymbioses with chemosynthetic bacteria. Endosymbionts serve essential biochemical and ecological functions, but the prokaryotic viruses (phages) that determine their fate are unknown.

RESULTS

We conducted metagenomic analysis of a deep-sea vent snail. We assembled four genome bins for Caudovirales phages that had developed dual endosymbiosis with sulphur-oxidising bacteria (SOB) and methane-oxidising bacteria (MOB). Clustered regularly interspaced short palindromic repeat (CRISPR) spacer mapping, genome comparison, and transcriptomic profiling revealed that phages Bin1, Bin2, and Bin4 infected SOB and MOB. The observation of prophages in the snail endosymbionts and expression of the phage integrase gene suggested the presence of lysogenic infection, and the expression of phage structural protein and lysozyme genes indicated active lytic infection. Furthermore, SOB and MOB appear to employ adaptive CRISPR-Cas systems to target phage DNA. Additional expressed defence systems, such as innate restriction-modification systems and dormancy-inducing toxin-antitoxin systems, may co-function and form multiple lines for anti-viral defence. To counter host defence, phages Bin1, Bin2, and Bin3 appear to have evolved anti-restriction mechanisms and expressed methyltransferase genes that potentially counterbalance host restriction activity. In addition, the high-level expression of the auxiliary metabolic genes narGH, which encode nitrate reductase subunits, may promote ATP production, thereby benefiting phage DNA packaging for replication.

CONCLUSIONS

This study provides new insights into phage-bacteria interplay in intracellular environments of a deep-sea vent snail. Video Abstract.

Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome

Anelloviruses are a ubiquitous component of healthy human viromes and remain highly prevalent after being acquired early in life. The full extent of "anellome" diversity and its evolutionary dynamics remain unexplored. We employed in-depth sequencing of blood-transfusion donor(s)-recipient pairs coupled with public genomic resources for a large-scale assembly of anellovirus genomes and used the data to characterize global and personal anellovirus diversity through time. The breadth of the anellome is much greater than previously appreciated, and individuals harbor unique anellomes and transmit lineages that can persist for several months within a diverse milieu of endemic host lineages. Anellovirus sequence diversity is shaped by extensive recombination at all levels of divergence, hindering traditional phylogenetic analyses. Our findings illuminate the transmission dynamics and vast diversity of anelloviruses and set the foundation for future studies to characterize their biology.

Expanding standards in viromics: in silico evaluation of dsDNA viral genome identification, classification, and auxiliary metabolic gene curation

BACKGROUND

Viruses influence global patterns of microbial diversity and nutrient cycles. Though viral metagenomics (viromics), specifically targeting dsDNA viruses, has been critical for revealing viral roles across diverse ecosystems, its analyses differ in many ways from those used for microbes. To date, viromics benchmarking has covered read pre-processing, assembly, relative abundance, read mapping thresholds and diversity estimation, but other steps would benefit from benchmarking and standardization. Here we use in silico-generated datasets and an extensive literature survey to evaluate and highlight how dataset composition (i.e., viromes vs bulk metagenomes) and assembly fragmentation impact (i) viral contig identification tool, (ii) virus taxonomic classification, and (iii) identification and curation of auxiliary metabolic genes (AMGs).

RESULTS

The in silico benchmarking of five commonly used virus identification tools show that gene-content-based tools consistently performed well for long (≥3 kbp) contigs, while k-mer- and blast-based tools were uniquely able to detect viruses from short (≤3 kbp) contigs. Notably, however, the performance increase of k-mer- and blast-based tools for short contigs was obtained at the cost of increased false positives (sometimes up to ∼5% for virome and ∼75% bulk samples), particularly when eukaryotic or mobile genetic element sequences were included in the test datasets. For viral classification, variously sized genome fragments were assessed using gene-sharing network analytics to quantify drop-offs in taxonomic assignments, which revealed correct assignations ranging from ∼95% (whole genomes) down to ∼80% (3 kbp sized genome fragments). A similar trend was also observed for other viral classification tools such as VPF-class, ViPTree and VIRIDIC, suggesting that caution is warranted when classifying short genome fragments and not full genomes. Finally, we highlight how fragmented assemblies can lead to erroneous identification of AMGs and outline a best-practices workflow to curate candidate AMGs in viral genomes assembled from metagenomes.

CONCLUSION

Together, these benchmarking experiments and annotation guidelines should aid researchers seeking to best detect, classify, and characterize the myriad viruses 'hidden' in diverse sequence datasets.

BACPHLIP: predicting bacteriophage lifestyle from conserved protein domains

Bacteriophages are broadly classified into two distinct lifestyles: temperate and virulent. Temperate phages are capable of a latent phase of infection within a host cell (lysogenic cycle), whereas virulent phages directly replicate and lyse host cells upon infection (lytic cycle). Accurate lifestyle identification is critical for determining the role of individual phage species within ecosystems and their effect on host evolution. Here, we present BACPHLIP, a BACterioPHage LIfestyle Predictor. BACPHLIP detects the presence of a set of conserved protein domains within an input genome and uses this data to predict lifestyle via a Random Forest classifier that was trained on a dataset of 634 phage genomes. On an independent test set of 423 phages, BACPHLIP has an accuracy of 98% greatly exceeding that of the previously existing tools (79%). BACPHLIP is freely available on GitHub (https://github.com/adamhockenberry/bacphlip) and the code used to build and test the classifier is provided in a separate repository (https://github.com/adamhockenberry/bacphlip-model-dev) for users wishing to interrogate and re-train the underlying classification model.

Seasonal dynamics of natural Ostreococcus viral infection at the single cell level using VirusFISH

Ostreococcus is a cosmopolitan marine genus of phytoplankton found in mesotrophic and oligotrophic waters, and the smallest free-living eukaryotes known to date, with a cell diameter close to 1 μm. Ostreococcus has been extensively studied as a model system to investigate viral-host dynamics in culture, yet the impact of viruses in naturally occurring populations is largely unknown. Here, we used Virus Fluorescence in situ Hybridization (VirusFISH) to visualize and quantify viral-host dynamics in natural populations of Ostreococcus during a seasonal cycle in the central Cantabrian Sea (Southern Bay of Biscay). Ostreococcus were predominantly found during summer and autumn at surface and 50 m depth, in coastal, mid-shelf and shelf waters, representing up to 21% of the picoeukaryotic communities. Viral infection was only detected in surface waters, and its impact was variable but highest from May to July and November to December, when up to half of the population was infected. Metatranscriptomic data available from the mid-shelf station unveiled that the Ostreococcus population was dominated by the species O. lucimarinus. This work represents a proof of concept that the VirusFISH technique can be used to quantify the impact of viruses on targeted populations of key microbes from complex natural communities.

A persistent giant algal virus, with a unique morphology, encodes an unprecedented number of genes involved in energy metabolism

Viruses have long been viewed as entities possessing extremely limited metabolic capacities. Over the last decade, however, this view has been challenged, as metabolic genes have been identified in viruses possessing large genomes and virions-the synthesis of which is energetically demanding. Here, we unveil peculiar phenotypic and genomic features of Prymnesium kappa virus RF01 (PkV RF01), a giant virus of the Mimiviridae family. We found that this virus encodes an unprecedented number of proteins involved in energy metabolism, such as all four succinate dehydrogenase (SDH) subunits (A-D) as well as key enzymes in the β-oxidation pathway. The SDHA gene was transcribed upon infection, indicating that the viral SDH is actively used by the virus- potentially to modulate its host's energy metabolism. We detected orthologous SDHA and SDHB genes in numerous genome fragments from uncultivated marine Mimiviridae viruses, which suggests that the viral SDH is widespread in oceans. PkV RF01 was less virulent compared with other cultured prymnesioviruses, a phenomenon possibly linked to the metabolic capacity of this virus and suggestive of relatively long co-evolution with its hosts. It also has a unique morphology, compared to other characterized viruses in the Mimiviridae family. Finally, we found that PkV RF01 is the only alga-infecting Mimiviridae virus encoding two aminoacyl-tRNA synthetases and enzymes corresponding to an entire base-excision repair pathway, as seen in heterotroph-infecting Mimiviridae These Mimiviridae encoded-enzymes were found to be monophyletic and branching at the root of the eukaryotic tree of life. This placement suggests that the last common ancestor of Mimiviridae was endowed with a large, complex genome prior to the divergence of known extant eukaryotes.IMPORTANCE Viruses on Earth are tremendously diverse in terms of morphology, functionality, and genomic composition. Over the last decade, the conceptual gap separating viruses and cellular life has tightened because of the detection of metabolic genes in viral genomes that express complex virus phenotypes upon infection. Here, we describe Prymnesium kappa virus RF01, a large alga-infecting virus with a unique morphology, an atypical infection profile, and an unprecedented number of genes involved in energy metabolism (such as the tricarboxylic (TCA) cycle and the β-oxidation pathway). Moreover, we show that the gene corresponding to one of these enzymes (the succinate dehydrogenase subunit A) is transcribed during infection and is widespread among marine viruses. This discovery provides evidence that a virus has the potential to actively regulate energy metabolism with its own gene.

Culturing novel and abundant pelagiphages in the ocean

Viruses play a key role in biogeochemical cycling and host mortality, metabolism, physiology and evolution in the ocean. Viruses that infect the globally abundant SAR11 bacteria (pelagiphages) were reported to be an important component of the marine viral communities. Our current knowledge of pelagiphages is based on a few studies and therefore is limited. In this study, 10 new pelagiphages were isolated and genomically characterized. These pelagiphages represent the first cultivated representatives of four viral lineages only found in metagenomic sequencing datasets previously. Many abundant environmental viral sequences, i.e., single-virus vSAG 37-F6 and several Global Ocean Viromes (GOV) viral populations, are now further confirmed with these pelagiphages. Viromic read mapping reveals that these new pelagiphages are globally distributed in the ocean and can be detected throughout the water column. Remarkably, isolation of these pelagiphages contributed up to 12% of all viromic reads annotated in the analysed viromes. Altogether, this study has greatly broadened our understanding of pelagiphages regarding their morphology, genetic diversity, infection strategies, and distribution pattern. The availability of these newly isolated pelagiphages and their genome sequences will allow us to further explore their infectivities and ecological strategies.