Natalie Solonenko

Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses

Ocean microbes drive biogeochemical cycling on a global scale. However, this cycling is constrained by viruses that affect community composition, metabolic activity, and evolutionary trajectories. Owing to challenges with the sampling and cultivation of viruses, genome-level viral diversity remains poorly described and grossly understudied, with less than 1% of observed surface-ocean viruses known. Here we assemble complete genomes and large genomic fragments from both surface- and deep-ocean viruses sampled during the Tara Oceans and Malaspina research expeditions, and analyse the resulting ‘global ocean virome’ dataset to present a global map of abundant, double-stranded DNA viruses complete with genomic and ecological contexts. A total of 15,222 epipelagic and mesopelagic viral populations were identified, comprising 867 viral clusters (defined as approximately genus-level groups). This roughly triples the number of known ocean viral populations and doubles the number of candidate bacterial and archaeal virus genera, providing a near-complete sampling of epipelagic communities at both the population and viral-cluster level. We found that 38 of the 867 viral clusters were locally or globally abundant, together accounting for nearly half of the viral populations in any global ocean virome sample. While two-thirds of these clusters represent newly described viruses lacking any cultivated representative, most could be computationally linked to dominant, ecologically relevant microbial hosts. Moreover, we identified 243 viral-encoded auxiliary metabolic genes, of which only 95 were previously known. Deeper analyses of four of these auxiliary metabolic genes (dsrC, soxYZ, P-II (also known as glnB) and amoC) revealed that abundant viruses may directly manipulate sulfur and nitrogen cycling throughout the epipelagic ocean. This viral catalog and functional analyses provide a necessary foundation for the meaningful integration of viruses into ecosystem models where they act as key players in nutrient cycling and trophic networks.

Twelve previously unknown phage genera are ubiquitous in global oceans

Viruses are fundamental to ecosystems ranging from oceans to humans, yet our ability to study them is bottlenecked by the lack of ecologically relevant isolates, resulting in “unknowns” dominating culture-independent surveys. Here we present genomes from 31 phages infecting multiple strains of the aquatic bacterium Cellulophaga baltica (Bacteroidetes) to provide data for an underrepresented and environmentally abundant bacterial lineage. Comparative genomics delineated 12 phage groups that (i) each represent a new genus, and (ii) represent one novel and four well-known viral families. This diversity contrasts the few well-studied marine phage systems, but parallels the diversity of phages infecting human-associated bacteria. Although all 12 Cellulophaga phages represent new genera, the podoviruses and icosahedral, nontailed ssDNA phages were exceptional, with genomes up to twice as large as those previously observed for each phage type. Structural novelty was also substantial, requiring experimental phage proteomics to identify 83% of the structural proteins. The presence of uncommon nucleotide metabolism genes in four genera likely underscores the importance of scavenging nutrient-rich molecules as previously seen for phages in marine environments. Metagenomic recruitment analyses suggest that these particular Cellulophaga phages are rare and may represent a first glimpse into the phage side of the rare biosphere. However, these analyses also revealed that these phage genera are widespread, occurring in 94% of 137 investigated metagenomes. Together, this diverse and novel collection of phages identifies a small but ubiquitous fraction of unknown marine viral diversity and provides numerous environmentally relevant phage–host systems for experimental hypothesis testing.

Single-cell and population level viral infection dynamics revealed by phageFISH, a method to visualize intracellular and free viruses

Microbes drive the biogeochemical cycles that fuel planet Earth, and their viruses (phages) alter microbial population structure, genome repertoire, and metabolic capacity. However, our ability to understand and quantify phage–host interactions is technique-limited. Here, we introduce phageFISH – a markedly improved geneFISH protocol that increases gene detection efficiency from 40% to > 92% and is optimized for detection and visualization of intra- and extracellular phage DNA. The application of phageFISH to characterize infection dynamics in a marine podovirus–gammaproteobacterial host model system corroborated classical metrics (qPCR, plaque assay, FVIC, DAPI) and outperformed most of them to reveal new biology. PhageFISH detected both replicating and encapsidated (intracellular and extracellular) phage DNA, while simultaneously identifying and quantifying host cells during all stages of infection. Additionally, phageFISH allowed per-cell relative measurements of phage DNA, enabling single-cell documentation of infection status (e.g. early vs late stage infections). Further, it discriminated between two waves of infection, which no other measurement could due to population-averaged signals. Together, these findings richly characterize the infection dynamics of a novel model phage–host system, and debut phageFISH as a much-needed tool for studying phage–host interactions in the laboratory, with great promise for environmental surveys and lineage-specific population ecology of free phages.

Towards quantitative viromics for both double-stranded and single-stranded DNA viruses

Background Viruses strongly influence microbial population dynamics and ecosystem functions. However, our ability to quantitatively evaluate those viral impacts is limited to the few cultivated viruses and double-stranded DNA (dsDNA) viral genomes captured in quantitative viral metagenomes (viromes). This leaves the ecology of non-dsDNA viruses nearly unknown, including single-stranded DNA (ssDNA) viruses that have been frequently observed in viromes, but not quantified due to amplification biases in sequencing library preparations (Multiple Displacement Amplification, Linker Amplification or Tagmentation). Methods Here we designed mock viral communities including both ssDNA and dsDNA viruses to evaluate the capability of a sequencing library preparation approach including an Adaptase step prior to Linker Amplification for quantitative amplification of both dsDNA and ssDNA templates. We then surveyed aquatic samples to provide first estimates of the abundance of ssDNA viruses. Results Mock community experiments confirmed the biased nature of existing library preparation methods for ssDNA templates (either largely enriched or selected against) and showed that the protocol using Adaptase plus Linker Amplification yielded viromes that were ±1.8-fold quantitative for ssDNA and dsDNA viruses. Application of this protocol to community virus DNA from three freshwater and three marine samples revealed that ssDNA viruses as a whole represent only a minor fraction (<5%) of DNA virus communities, though individual ssDNA genomes, both eukaryote-infecting Circular Rep-Encoding Single-Stranded DNA (CRESS-DNA) viruses and bacteriophages from the Microviridae family, can be among the most abundant viral genomes in a sample. Discussion Together these findings provide empirical data for a new virome library preparation protocol, and a first estimate of ssDNA virus abundance in aquatic systems.

Regulation of infection efficiency in a globally abundant marine Bacteriodetes virus

Bacteria impact humans, industry and nature, but do so under viral constraints. Problematically, knowledge of viral infection efficiencies and outcomes derives from few model systems that over-represent efficient lytic infections and under-represent virus–host natural diversity. Here we sought to understand infection efficiency regulation in an emerging environmental Bacteroidetes–virus model system with markedly different outcomes on two genetically and physiologically nearly identical host strains. For this, we quantified bacterial virus (phage) and host DNA, transcripts and phage particles throughout both infections. While phage transcriptomes were similar, transcriptional differences between hosts suggested host-derived regulation of infection efficiency. Specifically, the alternative host overexpressed DNA degradation genes and underexpressed translation genes, which seemingly targeted phage DNA particle production, as experiments revealed they were both significantly delayed (by >30 min) and reduced (by >50%) in the inefficient infection. This suggests phage failure to repress early alternative host expression and stress response allowed the host to respond against infection by delaying phage DNA replication and protein translation. Given that this phage type is ubiquitous and abundant in the global oceans and that variable viral infection efficiencies are central to dynamic ecosystems, these data provide a critically needed foundation for understanding and modeling viral infections in nature.

Ecogenomics and biogeochemical impacts of uncultivated globally abundant ocean viruses

Ocean microbes drive global-scale biogeochemical cycling1, but do so under constraints imposed by viruses on host community composition, metabolism, and evolutionary trajectories2–5. Due to sampling and cultivation challenges, genome-level viral diversity remains poorly described and grossly understudied in nature such that <1% of observed surface ocean viruses, even those that are abundant and ubiquitous, are ‘known’5. Here we analyze a global map of abundant, double stranded DNA (dsDNA) viruses and viral-encoded auxiliary metabolic genes (AMGs) with genomic and ecological contexts through the Global Ocean Viromes (GOV) dataset, which includes complete genomes and large genomic fragments from both surface and deep ocean viruses sampled during the Tara Oceans and Malaspina research expeditions6,7. A total of 15,222 epi- and mesopelagic viral populations were identified that comprised 867 viral clusters (VCs, approximately genus-level groups8,9). This roughly triples known ocean viral populations10, doubles known candidate bacterial and archaeal virus genera9, and near-completely samples epipelagic communities at both the population and VC level. Thirty-eight of the 867 VCs were identified as the most impactful dsDNA viral groups in the oceans, as these were locally or globally abundant and accounted together for nearly half of the viral populations in any GOV sample. Most of these were predicted in silico to infect dominant, ecologically relevant microbes, while two thirds of them represent newly described viruses that lacked any cultivated representative. Beyond these taxon-specific ecological observations, we identified 243 viral-encoded AMGs in GOV, only 95 of which were known. Deeper analyses of 4 of these AMGs revealed that abundant viruses directly manipulate sulfur and nitrogen cycling, and do so throughout the epipelagic ocean. Together these data provide a critically-needed organismal catalog and functional context to begin meaningfully integrating viruses into ecosystem models as key players in nutrient cycling and trophic networks.

Life-Style and Genome Structure of Marine Pseudoalteromonas Siphovirus B8b Isolated from the Northwestern Mediterranean Sea

Marine viruses (phages) alter bacterial diversity and evolution with impacts on marine biogeochemical cycles, and yet few well-developed model systems limit opportunities for hypothesis testing. Here we isolate phage B8b from the Mediterranean Sea using Pseudoalteromonas sp. QC-44 as a host and characterize it using myriad techniques. Morphologically, phage B8b was classified as a member of the Siphoviridae family. One-step growth analyses showed that this siphovirus had a latent period of 70 min and released 172 new viral particles per cell. Host range analysis against 89 bacterial host strains revealed that phage B8b infected 3 Pseudoalteromonas strains (52 tested, >99.9% 16S rRNA gene nucleotide identity) and 1 non-Pseudoaltermonas strain belonging to Alteromonas sp. (37 strains from 6 genera tested), which helps bound the phylogenetic distance possible in a phage-mediated horizontal gene transfer event. The Pseudoalteromonas phage B8b genome size was 42.7 kb, with clear structural and replication modules where the former were delineated leveraging identification of 16 structural genes by virion structural proteomics, only 4 of which had any similarity to known structural proteins. In nature, this phage was common in coastal marine environments in both photic and aphotic layers (found in 26.5% of available viral metagenomes), but not abundant in any sample (average per sample abundance was 0.65% of the reads). Together these data improve our understanding of siphoviruses in nature, and provide foundational information for a new ‘rare virosphere’ phage–host model system.

Optimization of viral resuspension methods for carbon-rich soils along a permafrost thaw gradient.

Permafrost stores approximately 50% of global soil carbon (C) in a frozen form; it is thawing rapidly under climate change, and little is known about viral communities in these soils or their roles in C cycling. In permafrost soils, microorganisms contribute significantly to C cycling, and characterizing them has recently been shown to improve prediction of ecosystem function. In other ecosystems, viruses have broad ecosystem and community impacts ranging from host cell mortality and organic matter cycling to horizontal gene transfer and reprogramming of core microbial metabolisms. Here we developed an optimized protocol to extract viruses from three types of high organic-matter peatland soils across a permafrost thaw gradient (palsa, moss-dominated bog, and sedge-dominated fen). Three separate experiments were used to evaluate the impact of chemical buffers, physical dispersion, storage conditions, and concentration and purification methods on viral yields. The most successful protocol, amended potassium citrate buffer with bead-beating or vortexing and BSA, yielded on average as much as 2-fold more virus-like particles (VLPs) g−1 of soil than other methods tested. All method combinations yielded VLPs g−1 of soil on the 108 order of magnitude across all three soil types. The different storage and concentration methods did not yield significantly more VLPs g−1 of soil among the soil types. This research provides much-needed guidelines for resuspending viruses from soils, specifically carbon-rich soils, paving the way for incorporating viruses into soil ecology studies.

Large-scale maps of variable infection efficiencies in aquatic Bacteroidetes phage-host model systems.

Microbes drive ecosystem functioning and their viruses modulate these impacts through mortality, gene transfer and metabolic reprogramming. Despite the importance of virus-host interactions and likely variable infection efficiencies of individual phages across hosts, such variability is seldom quantified. Here, we quantify infection efficiencies of 38 phages against 19 host strains in aquatic Cellulophaga (Bacteroidetes) phage-host model systems. Binary data revealed that some phages infected only one strain while others infected 17, whereas quantitative data revealed that efficiency of infection could vary 10 orders of magnitude, even among phages within one population. This provides a baseline for understanding and modeling intrapopulation host range variation. Genera specific host ranges were also informative. For example, the Cellulophaga Microviridae, showed a markedly broader intra-species host range than previously observed in Escherichia coli systems. Further, one phage genus, Cba41, was examined to investigate nonheritable changes in plating efficiency and burst size that depended on which host strain it most recently infected. While consistent with host modification of phage DNA, no differences in nucleotide sequence or DNA modifications were detected, leaving the observation repeatable, but the mechanism unresolved. Overall, this study highlights the importance of quantitatively considering replication variations in studies of phage-host interactions.

Comparative Omics and Trait Analyses of Marine Pseudoalteromonas Phages Advance the Phage OTU Concept

Viruses influence the ecology and evolutionary trajectory of microbial communities. Yet our understanding of their roles in ecosystems is limited by the paucity of model systems available for hypothesis generation and testing. Further, virology is limited by the lack of a broadly accepted conceptual framework to classify viral diversity into evolutionary and ecologically cohesive units. Here, we introduce genomes, structural proteomes, and quantitative host range data for eight Pseudoalteromonas phages isolated from Helgoland (North Sea, Germany) and use these data to advance a genome-based viral operational taxonomic unit (OTU) definition. These viruses represent five new genera and inform 498 unaffiliated or unannotated protein clusters (PCs) from global virus metagenomes. In a comparison of previously sequenced Pseudoalteromonas phage isolates (n = 7) and predicted prophages (n = 31), the eight phages are unique. They share a genus with only one other isolate, Pseudoalteromonas podophage RIO-1 (East Sea, South Korea) and two Pseudoalteromonas prophages. Mass-spectrometry of purified viral particles identified 12–20 structural proteins per phage. When combined with 3-D structural predictions, these data led to the functional characterization of five previously unidentified major capsid proteins. Protein functional predictions revealed mechanisms for hijacking host metabolism and resources. Further, they uncovered a hybrid sipho-myovirus that encodes genes for Mu-like infection rarely described in ocean systems. Finally, we used these data to evaluate a recently introduced definition for virus populations that requires members of the same population to have >95% average nucleotide identity across at least 80% of their genes. Using physiological traits and genomics, we proposed a conceptual model for a viral OTU definition that captures evolutionarily cohesive and ecologically distinct units. In this trait-based framework, sensitive hosts are considered viral niches, while host ranges and infection efficiencies are tracked as viral traits. Quantitative host range assays revealed conserved traits within virus OTUs that break down between OTUs, suggesting the defined units capture niche and fitness differentiation. Together these analyses provide a foundation for model system-based hypothesis testing that will improve our understanding of marine copiotrophs, as well as phage–host interactions on the ocean particles and aggregates where Pseudoalteromonas thrive.