Single-cell genomics-based analysis of virus–host interactions in marine surface bacterioplankton

Viral infections dynamically alter the composition and metabolic potential of marine microbial communities and the evolutionary trajectories of host populations with resulting feedback on biogeochemical cycles. It is quite possible that all microbial populations in the ocean are impacted by viral infections. Our knowledge of virus–host relationships, however, has been limited to a minute fraction of cultivated host groups. Here, we utilized single-cell sequencing to obtain genomic blueprints of viruses inside or attached to individual bacterial and archaeal cells captured in their native environment, circumventing the need for host and virus cultivation. A combination of comparative genomics, metagenomic fragment recruitment, sequence anomalies and irregularities in sequence coverage depth and genome recovery were utilized to detect viruses and to decipher modes of virus–host interactions. Members of all three tailed phage families were identified in 20 out of 58 phylogenetically and geographically diverse single amplified genomes (SAGs) of marine bacteria and archaea. At least four phage–host interactions had the characteristics of late lytic infections, all of which were found in metabolically active cells. One virus had genetic potential for lysogeny. Our findings include first known viruses of Thaumarchaeota, Marinimicrobia, Verrucomicrobia and Gammaproteobacteria clusters SAR86 and SAR92. Viruses were also found in SAGs of Alphaproteobacteria and Bacteroidetes. A high fragment recruitment of viral metagenomic reads confirmed that most of the SAG-associated viruses are abundant in the ocean. Our study demonstrates that single-cell genomics, in conjunction with sequence-based computational tools, enable in situ, cultivation-independent insights into host–virus interactions in complex microbial communities.     

Efficient dilution-to-extinction isolation of novel virus–host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus–host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus–host systems in culture collections and genomic databases will improve both our understanding of virus–host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.Subject terms: Bacteriophages, Microbial ecology     

Top-down controls on bacterial community structure: microbial network analysis of bacteria,T4-like viruses and protists

Characterizing ecological relationships between viruses, bacteria and protists in the ocean are critical to understanding ecosystem function, yet these relationships are infrequently investigated together. We evaluated these relationships through microbial association network analysis of samples collected approximately monthly from March 2008 to January 2011 in the surface ocean (0–5 m) at the San Pedro Ocean Time series station. Bacterial, T4-like myoviral and protistan communities were described by Automated Ribosomal Intergenic Spacer Analysis and terminal restriction fragment length polymorphism of the gene encoding the major capsid protein (g23) and 18S ribosomal DNA, respectively. Concurrent shifts in community structure suggested similar timing of responses to environmental and biological parameters. We linked T4-like myoviral, bacterial and protistan operational taxonomic units by local similarity correlations, which were then visualized as association networks. Network links (correlations) potentially represent synergistic and antagonistic relationships such as viral lysis, grazing, competition or other interactions. We found that virus–bacteria relationships were more cross-linked than protist–bacteria relationships, suggestive of increased taxonomic specificity in virus–bacteria relationships. We also found that 80% of bacterial–protist and 74% of bacterial–viral correlations were positive, with the latter suggesting that at monthly and seasonal timescales, viruses may be following their hosts more often than controlling host abundance.     

Deep sea sediments associated with cold seeps are a subsurface reservoir of viral diversity

In marine ecosystems, viruses exert control on the composition and metabolism of microbial communities, influencing overall biogeochemical cycling. Deep sea sediments associated with cold seeps are known to host taxonomically diverse microbial communities, but little is known about viruses infecting these microorganisms. Here, we probed metagenomes from seven geographically diverse cold seeps across global oceans to assess viral diversity, virus–host interaction, and virus-encoded auxiliary metabolic genes (AMGs). Gene-sharing network comparisons with viruses inhabiting other ecosystems reveal that cold seep sediments harbour considerable unexplored viral diversity. Most cold seep viruses display high degrees of endemism with seep fluid flux being one of the main drivers of viral community composition. In silico predictions linked 14.2% of the viruses to microbial host populations with many belonging to poorly understood candidate bacterial and archaeal phyla. Lysis was predicted to be a predominant viral lifestyle based on lineage-specific virus/host abundance ratios. Metabolic predictions of prokaryotic host genomes and viral AMGs suggest that viruses influence microbial hydrocarbon biodegradation at cold seeps, as well as other carbon, sulfur and nitrogen cycling via virus-induced mortality and/or metabolic augmentation. Overall, these findings reveal the global diversity and biogeography of cold seep viruses and indicate how viruses may manipulate seep microbial ecology and biogeochemistry.Subject terms: Environmental microbiology, Microbial ecology     

Metagenomic Exploration of Viruses throughout the Indian Ocean

The characterization of global marine microbial taxonomic and functional diversity is a primary goal of the Global Ocean Sampling Expedition. As part of this study, 19 water samples were collected aboard the Sorcerer II sailing vessel from the southern Indian Ocean in an effort to more thoroughly understand the lifestyle strategies of the microbial inhabitants of this ultra-oligotrophic region. No investigations of whole virioplankton assemblages have been conducted on waters collected from the Indian Ocean or across multiple size fractions thus far. Therefore, the goals of this study were to examine the effect of size fractionation on viral consortia structure and function and understand the diversity and functional potential of the Indian Ocean virome. Five samples were selected for comprehensive metagenomic exploration; and sequencing was performed on the microbes captured on 3.0-, 0.8- and 0.1 µm membrane filters as well as the viral fraction (<0.1 µm). Phylogenetic approaches were also used to identify predicted proteins of viral origin in the larger fractions of data from all Indian Ocean samples, which were included in subsequent metagenomic analyses. Taxonomic profiling of viral sequences suggested that size fractionation of marine microbial communities enriches for specific groups of viruses within the different size classes and functional characterization further substantiated this observation. Functional analyses also revealed a relative enrichment for metabolic proteins of viral origin that potentially reflect the physiological condition of host cells in the Indian Ocean including those involved in nitrogen metabolism and oxidative phosphorylation. A novel classification method, MGTAXA, was used to assess virus-host relationships in the Indian Ocean by predicting the taxonomy of putative host genera, with Prochlorococcus, Acanthochlois and members of the SAR86 cluster comprising the most abundant predictions. This is the first study to holistically explore virioplankton dynamics across multiple size classes and provides unprecedented insight into virus diversity, metabolic potential and virus-host interactions.     

Links between viruses and prokaryotes throughout the water column along a North Atlantic latitudinal transect

Viruses are an abundant, diverse and dynamic component of marine ecosystems and have a key role in the biogeochemical processes of the ocean by controlling prokaryotic and phytoplankton abundance and diversity. However, most of the studies on virus–prokaryote interactions in marine environments have been performed in nearshore waters. To assess potential variations in the relation between viruses and prokaryotes in different oceanographic provinces, we determined viral and prokaryotic abundance and production throughout the water column along a latitudinal transect in the North Atlantic. Depth-related trends in prokaryotic and viral abundance (both decreasing by one order of magnitude from epi- to abyssopelagic waters), and prokaryotic production (decreasing by three orders of magnitude) were observed along the latitudinal transect. The virus-to-prokaryote ratio (VPR) increased from ∼19 in epipelagic to ∼53 in the bathy- and abyssopelagic waters. Although the lytic viral production decreased significantly with depth, the lysogenic viral production did not vary with depth. In bathypelagic waters, pronounced differences in prokaryotic and viral abundance were found among different oceanic provinces with lower leucine incorporation rates and higher VPRs in the North Atlantic Gyre province than in the provinces further north and south. The percentage of lysogeny increased from subpolar regions toward the more oligotrophic lower latitudes. Based on the observed trends over this latitudinal transect, we conclude that the viral–host interactions significantly change among different oceanic provinces in response to changes in the biotic and abiotic variables.     

Bacterial Growth Rate and Marine Virus–Host Dynamics

Abstract The dynamics of a marine virus–host system were investigated at different steady state growth rates in chemostat cultures and the data were analyzed using a simple model. The virus–host interactions showed strong dependence on host cell growth rate. The duration of the infection cycle and the virus burst size were found to depend on bacterial growth rate, and the rate of cell lysis and virus production were positively correlated with steady state growth rate in the cultures (r ² > 0.96, p < 0.05). At bacterial growth rates of 0.02 to 0.10 h⁻¹ in the chemostats the virus burst size increased from 12 ± 4 to 56 ± 4, and the latent period decreased from 2.0 to 1.7 h. Resistant clones of the host strain were present in the cultures from the beginning of the experiment and replaced the sensitive host cells following viral lysis in the cultures. Regrowth of resistant cells correlated significantly (r ²= 1.000, p < 0.02) with the lysis rate of sensitive cells, indicating that release of viral lysates stimulated growth of the non-infected, resistant cells. The constructed model was suitable for simulating the observed dynamics of the sensitive host cells, viruses and resistant clones in the cultures. The model was therefore used in an attempt to predict the dynamics of this virus–host interaction in a natural marine environment during a certain set of growth conditions. The simulation indicated that a steady state relationship between the specific viruses and sensitive and resistant bacterial clones may occur at densities that are reasonable to assume for natural environments. The study demonstrates that basic characterization and modeling of specific virus–host interactions may improve our understanding of the behavior of bacteria and viruses in natural systems. Received: 12 November 1999; Accepted: 2 May 2000; Online Publication: 11 August 2000     

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.

Comparative genomics and host range analysis reveals the remarkable diversity and evolution of tailed archaeal viruses of the order Caudoviricetes, which together with their bacterial relatives arguably represent the most abundant and widespread virus group on our planet.     

Indirect interactions in the microbial world: specificities and similarities to plant–insect systems

Trophic interactions between bacteria, viruses, and protozoan predators play crucial roles in structuring aquatic microbial communities and regulating microbe-mediated ecosystem functions (biogeochemical processes). In this microbial food web, protozoan predators and viruses share bacteria as a common resource, and protozoan predators can kill viruses [intraguild predation (IGP)] and vice versa, even though these latter processes are probably of less importance. However, protozoan predators (IG predator) and viruses (IG prey) generally occur together in various environments, and this cannot be fully explained by the classic IGP models. In addition, controlled experiments have often demonstrated that protozoan predators have apparently positive effects on viral activity. These surprising patterns can be explained by indirect interactions between them via induced trait changes in bacterial assemblages, which can be compared with trait-mediated indirect interactions (TMIIs) in terrestrial plant–insect systems. Here, we review some trait changes in bacterial assemblages that may positively affect the activities and abundance of viruses. It has been suggested that in bacterial assemblages, protozoan predation may enhance growth conditions for individual bacteria and induce both phenotypic trait changes at the individual (e.g., filament-forming bacteria) and group level as a result of changes in bacterial community composition (e.g., species dominance). We discuss the specificities of aquatic microbial systems and attempt find functional similarities between aquatic microbial systems and terrestrial plant–insect systems with regard to TMII function.     

High potential for temperate viruses to drive carbon cycling in chemoautotrophy‐dominated shallow‐water hydrothermal vents

Viruses are the most abundant life forms in the world's oceans and they are key drivers of biogeochemical cycles, but their impact on the microbial assemblages inhabiting hydrothermal vent ecosystems is still largely unknown. Here, we analysed the viral life strategies and virus‐host interactions in the sediments of a newly discovered shallow‐water hydrothermal field of the Mediterranean Sea. Our study reveals that temperate viruses, once experimentally induced to replicate, can cause large mortality of vent microbes, significantly reducing the chemoautotrophic carbon production, while enhancing the metabolism of microbial heterotrophs and the re‐cycling of the organic matter. These results provide new insights on the factors controlling primary and secondary production processes in hydrothermal vents, suggesting that the inducible provirus‐host interactions occurring in these systems can profoundly influence the functioning of the microbial food web and the efficiency in the energy transfer to the higher trophic levels.     

Viruses,cancer and non-self recognition

Virus–host interactions form an essential part of every aspect of life, and this review is aimed at looking at the balance between the host and persistent viruses with a focus on the immune system. The virus–host interaction is like a cat-and-mouse game and viruses have developed ingenious mechanisms to manipulate cellular pathways, most notably the major histocompatibility (MHC) class I pathway, to reside within infected cell while evading detection and destruction by the immune system. However, some of the signals sensing and responding to viral infection are derived from viruses and the fact that certain viruses can prevent the infection of others, highlights a more complex coexistence between the host and the viral microbiota. Viral immune evasion strategies also illustrate that processes whereby cells detect and present non-self genetic material to the immune system are interlinked with other cellular pathways. Immune evasion is a target also for cancer cells and a more detailed look at the interfaces between viral factors and components of the MHC class I peptide-loading complex indicates that these interfaces are also targets for cancer mutations. In terms of the immune checkpoint, however, viral and cancer strategies appear different.     

Transitional states in marine fisheries: adapting to predicted global change

Global climate change has the potential to substantially alter the production and community structure of marine fisheries and modify the ongoing impacts of fishing. Fish community composition is already changing in some tropical, temperate and polar ecosystems, where local combinations of warming trends and higher environmental variation anticipate the changes likely to occur more widely over coming decades. Using case studies from the Western Indian Ocean, the North Sea and the Bering Sea, we contextualize the direct and indirect effects of climate change on production and biodiversity and, in turn, on the social and economic aspects of marine fisheries. Climate warming is expected to lead to (i) yield and species losses in tropical reef fisheries, driven primarily by habitat loss; (ii) community turnover in temperate fisheries, owing to the arrival and increasing dominance of warm-water species as well as the reduced dominance and departure of cold-water species; and (iii) increased diversity and yield in Arctic fisheries, arising from invasions of southern species and increased primary production resulting from ice-free summer conditions. How societies deal with such changes will depend largely on their capacity to adapt--to plan and implement effective responses to change--a process heavily influenced by social, economic, political and cultural conditions.     

Characterization of Lysogens in Bacterioplankton Assemblages of the Southern California Borderland

Viruses cause significant mortality of marine microorganisms; however, their role in shaping the composition of microbial assemblages has not been fully elucidated. Because viruses may form lysogenic relationships with their hosts, temperate viruses may influence bacterial assemblage structures through direct lysis of hosts when induced by environmental stimuli or by homoimmunity (i.e., immunity to closely related viruses). We investigated the components of bacterioplankton assemblages that bore prophage using the lysogenic induction agent mitomycin C. Seawater was collected at two locations (the San Pedro Ocean Time Series Station and in the Santa Barbara Channel) in the Southern California Borderland and amended with mitomycin C. After 24-h incubation, the community structure of bacterioplankton was compared with unamended controls using automated rRNA intergenic spacer analysis. The addition of mitomycin C to seawater had effects on the community structure of bacterioplankton, stimulating detectable overall diversity and richness of fingerprints and causing the assemblages within incubations to become different to control assemblages. Most negatively impacted operational taxonomic units (OTU) in mitomycin C-amended incubations individually comprised a large fraction of total amplified DNA in initial seawater (5.3-23.3% of amplified DNA fluorescence) fingerprints, and data suggest that these include organisms putatively classified as members of the gamma-Proteobacteria, SAR11 cluster, and Synechococcus groups. The stimulation of assemblage richness by induction of lysogens, and the reduction in the contribution to total DNA of common OTU (and concomitant increase in rare OTU), suggests that temperate phage have the potential to strongly influence the diversity of bacterioplankton assemblages. Because lysogenic OTU may also be resistant to closely related lytic (i.e., free-living) viruses, the impact of lytic virioplankton on assemblages may only be pronounced transiently or when conditions causing lysogenic induction arise.     

Psychrophilic versus psychrotolerant bacteria--occurrence and significance in polar and temperate marine habitats.

The numerical dominance and ecological role of psychrophilic bacteria in bottom sediments, sea ice, surface water and melt pools of the polar oceans were investigated using isolates, colony forming units (CFU) and metabolic activities. All sediment samples of the Southern Ocean studied showed a clear numerical dominance of cold-loving bacteria. In Arctic sediments underlying the influence of cold polar water bodies psychrophiles prevailed also but they were less dominant in sediments influenced by the warm Atlantic Water. A predominance of psychrophiles was further found in consolidated Antarctic sea ice as well as in multiyear Arctic sea ice and in melt pools on top of Arctic ice floes. A less uniform adaptation response was, however, met in polar surface waters. In the very northern part of the Fram Strait (Arctic Ocean) we found bacterial counts and activities at 1 degree C exceeding those at 22 degrees C. In surface water of the Weddell Sea (Southern Ocean) psychrophiles also dominated numerically in early autumn but the dominance declined obviously with the onset of winter-water and a decrease of chlorphyll a. Otherwise in surface water of the Southern Ocean CFUs were higher at 22 degrees C than at 1 degree C while activities were vice versa indicating at least a functional dominance of psychrophiles. Even in the temperate sediments of the German Bight true psychrophiles were present and a clear shift towards cold adapted communities in winter observed. Among the polar bacteria a more pronounced cold adaptation of Antarctic in comparison with Arctic isolates was obtained. The results and literature data indicate that stenothermic cold adapted bacteria play a significant role in the global marine environment. On the basis of the temperature response of our isolates from different habitats it is suggested to expand the definition of Morita in order to meet the cold adaptation strategies of the bacteria in the various cold habitats.     

Abundance and function of bacteria in the Southern Ocean.

The very low water temperatures existing in polar oceans that experience seasonal advance and retreat of pack ice do not inhibit the presence of large bacterial populations. Bacteria may contribute significantly to the energy transfers within the Southern Ocean. In the last decades, notable progress has been made in the knowledge of the role of marine bacteria in the Southern Ocean. A short overview of the abundance and function ofAntarctic marine bacteria is given, with respect to metabolic activity. The importance of spatial and temporal variability is described. The ecological function of Antarctic marine bacterioplankton is discussed. Depending on food web structure, bacteria may be either a link in food webs supporting metazoan production, or a sink where bacterial production is metabolised by microorganisms. In the more oligotrophic areas and during certain periods of the year bacterial biomass dominates phytoplankton. The microbial food web is therefore the dominant pathway for carbon and energy flow in Antarctic seawater.     

Comparison of contemporary and fossil diatom assemblages from the western Antarctic Peninsula shelf

During the past ten years, the Antarctic Peninsula has been identified as the most rapidly warming region of the Southern Hemisphere and it is important to place this warming in the context of the natural climate and oceanographic variability of the recent geological past. Many biological proxies, such as marine diatom assemblages, have been used to determine Southern Ocean palaeoceanographic conditions during the Late Quaternary, however, few investigations have attempted to link observations of modern floras with the fossil record. In this study we examine a suite of modern austral spring (December 2003) and summer (February 2002) surface water samples from along the western Antarctic Peninsula (WAP) continental shelf and compare these to core-top, surface sediment samples. Using detrended correspondence analysis (DCA) and principal component analysis (PCA) of diatom abundance data we investigate the relationship of contemporary diatom floras with the fossil record. This multivariate analysis reveals that our modern assemblages can be divided into three groups: summer southern WAP sites, summer northern WAP sites, and spring WAP sites. Sea surface temperature (SST) is an important environmental variable for explaining seasonal differences in diatom assemblages between spring and summer, but sea surface salinity (SSS) is more important for understanding temporally-equivalent regional variations in assemblage. Our summer diatom samples are more reminiscent of early season assemblages, reflecting the unusually late sea ice retreat from the region that year. When the modern assemblages are compared to the fossil record, it is clear that most of the important diatoms from the summer assemblage are not preserved into the sediments, and that the fossil record more closely reflects spring assemblages. This observation is important for any future attempts to quantitatively reconstruct palaeoceanographic conditions along the WAP during the Late Quaternary and highlights the need for many more such studies in order to address longer timescales, such as interannual variability, in the context of the fossil record.     

Comprehensive analysis of endogenous bornavirus-like elements in eukaryote genomes

Bornaviruses are the only animal RNA viruses that establish a persistent infection in their host cell nucleus. Studies of bornaviruses have provided unique information about viral replication strategies and virus–host interactions. Although bornaviruses do not integrate into the host genome during their replication cycle, we and others have recently reported that there are DNA sequences derived from the mRNAs of ancient bornaviruses in the genomes of vertebrates, including humans, and these have been designated endogenous borna-like (EBL) elements. Therefore, bornaviruses have been interacting with their hosts as driving forces in the evolution of host genomes in a previously unexpected way. Studies of EBL elements have provided new models for virology, evolutionary biology and general cell biology. In this review, we summarize the data on EBL elements including what we have newly identified in eukaryotes genomes, and discuss the biological significance of EBL elements, with a focus on EBL nucleoprotein elements in mammalian genomes. Surprisingly, EBL elements were detected in the genomes of invertebrates, suggesting that the host range of bornaviruses may be much wider than previously thought. We also review our new data on non-retroviral integration of Borna disease virus.     

Metabolic reprogramming by viruses in the sunlit and dark ocean

Background

Marine ecosystem function is largely determined by matter and energy transformations mediated by microbial community interaction networks. Viral infection modulates network properties through mortality, gene transfer and metabolic reprogramming.

Results

Here we explore the nature and extent of viral metabolic reprogramming throughout the Pacific Ocean depth continuum. We describe 35 marine viral gene families with potential to reprogram metabolic flux through central metabolic pathways recovered from Pacific Ocean waters. Four of these families have been previously reported but 31 are novel. These known and new carbon pathway auxiliary metabolic genes were recovered from a total of 22 viral metagenomes in which viral auxiliary metabolic genes were differentiated from low-level cellular DNA inputs based on small subunit ribosomal RNA gene content, taxonomy, fragment recruitment and genomic context information. Auxiliary metabolic gene distribution patterns reveal that marine viruses target overlapping, but relatively distinct pathways in sunlit and dark ocean waters to redirect host carbon flux towards energy production and viral genome replication under low nutrient, niche-differentiated conditions throughout the depth continuum.

Conclusions

Given half of ocean microbes are infected by viruses at any given time, these findings of broad viral metabolic reprogramming suggest the need for renewed consideration of viruses in global ocean carbon models.