Trophic strategies explain the ocean niches of small eukaryotic phytoplankton.

A large fraction of marine primary production is performed by diverse small protists, and many of these phytoplankton are phagotrophic mixotrophs that vary widely in their capacity to consume bacterial prey. Prior analyses suggest that mixotrophic protists as a group vary in importance across ocean environments, but the mechanisms leading to broad functional diversity among mixotrophs, and the biogeochemical consequences of this, are less clear. Here we use isolates from seven major taxa to demonstrate a tradeoff between phototrophic performance (growth in the absence of prey) and phagotrophic performance (clearance rate when consuming Prochlorococcus). We then show that trophic strategy along the autotrophy-mixotrophy spectrum correlates strongly with global niche differences, across depths and across gradients of stratification and chlorophyll a. A model of competition shows that community shifts can be explained by greater fitness of faster-grazing mixotrophs when nutrients are scarce and light is plentiful. Our results illustrate how basic physiological constraints and principles of resource competition can organize complexity in the surface ocean ecosystem.

CoCoNet: an efficient deep learning tool for viral metagenome binning.

MOTIVATION: Metagenomic approaches hold the potential to characterize microbial communities and unravel the intricate link between the microbiome and biological processes. Assembly is one of the most critical steps in metagenomics experiments. It consists of transforming overlapping DNA sequencing reads into sufficiently accurate representations of the community's genomes. This process is computationally difficult and commonly results in genomes fragmented across many contigs. Computational binning methods are used to mitigate fragmentation by partitioning contigs based on their sequence composition, abundance or chromosome organization into bins representing the community's genomes. Existing binning methods have been principally tuned for bacterial genomes and do not perform favorably on viral metagenomes. RESULTS: We propose Composition and Coverage Network (CoCoNet), a new binning method for viral metagenomes that leverages the flexibility and the effectiveness of deep learning to model the co-occurrence of contigs belonging to the same viral genome and provide a rigorous framework for binning viral contigs. Our results show that CoCoNet substantially outperforms existing binning methods on viral datasets. AVAILABILITY AND IMPLEMENTATION: CoCoNet was implemented in Python and is available for download on PyPi (https://pypi.org/). The source code is hosted on GitHub at https://github.com/Puumanamana/CoCoNet and the documentation is available at https://coconet.readthedocs.io/en/latest/index.html. CoCoNet does not require extensive resources to run. For example, binning 100k contigs took about 4 h on 10 Intel CPU Cores (2.4 GHz), with a memory peak at 27 GB (see Supplementary Fig. S9). To process a large dataset, CoCoNet may need to be run on a high RAM capacity server. Such servers are typically available in high-performance or cloud computing settings. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Variable Freshwater Influences on the Abundance of Vibrio vulnificus in a Tropical Urban Estuary.

To better understand the controls on the opportunistic human pathogen Vibrio vulnificus in warm tropical waters, we conducted a year-long investigation in the Ala Wai Canal, a channelized estuary in Honolulu, HI. The abundance of V. vulnificus, as determined by quantitative PCR (qPCR) of the hemolysin gene (vvhA), varied spatially and temporally by nearly 4 orders of magnitude (≤3 to 14,000 mL-1). Unlike in temperate and subtropical systems, temperatures were persistently warm (19 to 31°C) and explained little of the variability in V. vulnificus abundance. Salinity (1 to 36 ppt) had a significant, but nonlinear, relationship with V. vulnificus abundance with the highest vvhA concentrations (>2,500 mL-1) observed only at salinities from 7 to 22 ppt. V. vulnificus abundances were lower on average during the summer dry season, when waters were warmer but more saline. The highest canal-wide average abundances were observed during a time of modest rainfall, when moderate salinities and elevated concentrations of reduced nitrogen species and silica suggested a groundwater influence. Parallel quantification of the vcgC gene suggested that C-type strains, which are responsible for most human infections, comprised 25% of the total V. vulnificus on average, but their relative contribution was greater at higher salinities, suggesting a broader salinity tolerance. Generalized regression models suggested that up to 67% of sample-to-sample variation (n = 202) in log-transformed V. vulnificus abundance was explained using the measured environmental variables, and up to 97% of the monthly variation in canal-wide average concentrations (n = 13) was explained with the best subset of four variables. IMPORTANCE Our data illustrate that, in the absence of strong seasonal variation in water temperature in the tropics, variation in salinity driven by rainfall becomes a primary controlling variable on V. vulnificus abundance. There is thus a tendency for a rainfall-driven seasonal cycle in V. vulnificus abundance which is inverted from the temperature-driven seasonal cycle at higher latitudes. However, stochasticity in rainfall and its nonlinear, indirect effects on V. vulnificus concentration means that high abundances can occur at any location in the canal at any time of year, making it challenging to predict concentrations of this pathogen at a high temporal or spatial resolution. Much of the variability in canal-wide average concentrations, on the other hand, was explained by a few variables that reflect the magnitude of freshwater input to the system, suggesting that relative risk of exposure to this pathogen could be predicted as an average for the system.

Making sense of virus size and the tradeoffs shaping viral fitness.

Viruses span an impressive size range, with genome length varying a thousandfold and virion volume nearly a millionfold. For cellular organisms the scaling of traits with size is a pervasive influence on ecological processes, but whether size plays a central role in viral ecology is unknown. Here, we focus on viruses of aquatic unicellular organisms, which exhibit the greatest known range of virus size. We outline hypotheses within a quantitative framework, and analyse data where available, to consider how size affects the primary components of viral fitness. We argue that larger viruses have fewer offspring per infection and slower contact rates with host cells, but a larger genome tends to increase infection efficiency, broaden host range, and potentially increase attachment success and decrease decay rate. These countervailing selective pressures may explain why a breadth of sizes exist and even coexist when infecting the same host populations. Oligotrophic ecosystems may be enriched in "giant" viruses, because environments with resource-limited phagotrophs at low concentrations may select for broader host range, better control of host metabolism, lower decay rate and a physical size that mimics bacterial prey. Finally, we describe where further research is needed to understand the ecology and evolution of viral size diversity.

Host Traits Drive Viral Life Histories across Phytoplankton Viruses.

Viruses are integral to ecological and evolutionary processes, but we have a poor understanding of what drives variation in key traits across diverse viruses. For lytic viruses, burst size, latent period, and genome size are primary characteristics controlling host-virus dynamics. Here we synthesize data on these traits for 75 strains of phytoplankton viruses, which play an important role in global biogeochemistry. We find that primary traits of the host (genome size, growth rate) explain 40%-50% of variation in burst size and latent period. Specifically, burst size and latent period both exhibit saturating relationships versus the host∶virus genome size ratio, with both traits increasing at low genome size ratios while showing no relationship at high size ratios. In addition, latent period declines as host growth rate increases. We analyze a model of latent period evolution to explore mechanisms that could cause these patterns. The model predicts that burst size may often be set by the host genomic resources available for viral construction, while latent period evolves to permit this maximal burst size, modulated by host metabolic rate. These results suggest that general mechanisms may underlie the evolution of diverse viruses. Future extensions of this work could help explain viral regulation of host populations, viral influence on community structure and diversity, and viral roles in biogeochemical cycles.

RNA viruses as major contributors to Antarctic virioplankton.

Early work on marine algal viruses focused exclusively on those having DNA genomes, but recent studies suggest that RNA viruses, especially those with positive-sense, single-stranded RNA (+ssRNA) genomes, are abundant in tropical and temperate coastal seawater. To test whether this was also true of polar waters, we estimated the relative abundances of RNA and DNA viruses using a mass ratio approach and conducted shotgun metagenomics on purified viral samples collected from a coastal site near Palmer Station, Antarctica on six occasions throughout a summer phytoplankton bloom (November-March). Our data suggest that RNA viruses contributed up to 65% of the total virioplankton (8-65%), and that, as observed previously in warmer waters, the majority of RNA viruses in these Antarctic RNA virus metagenomes had +ssRNA genomes most closely related to viruses in the order Picornavirales. Assembly of the metagenomic reads resulted in five novel, nearly complete genomes, three of which had features similar to diatom-infecting viruses. Our data are consistent with the hypothesis that RNA viruses influence diatom bloom dynamics in Antarctic waters.

Variables influencing extraction of nucleic acids from microbial plankton (viruses, bacteria, and protists) collected on nanoporous aluminum oxide filters.

Anodic aluminum oxide (AAO) filters have high porosity and can be manufactured with a pore size that is small enough to quantitatively capture viruses. These properties make the filters potentially useful for harvesting total microbial communities from water samples for molecular analyses, but their performance for nucleic acid extraction has not been systematically or quantitatively evaluated. In this study, we characterized the flux of water through commercially produced nanoporous (0.02 µm) AAO filters (Anotop; Whatman) and used isolates (a virus, a bacterium, and a protist) and natural seawater samples to test variables that we expected would influence the efficiency with which nucleic acids are recovered from the filters. Extraction chemistry had a significant effect on DNA yield, and back flushing the filters during extraction was found to improve yields of high-molecular-weight DNA. Using the back-flush protocol, the mass of DNA recovered from microorganisms collected on AAO filters was ≥ 100% of that extracted from pellets of cells and viruses and 94% ± 9% of that obtained by direct extraction of a liquid bacterial culture. The latter is a minimum estimate of the relative recovery of microbial DNA, since liquid cultures include dissolved nucleic acids that are retained inefficiently by the filter. In conclusion, we demonstrate that nucleic acids can be extracted from microorganisms on AAO filters with an efficiency similar to that achievable by direct extraction of microbes in suspension or in pellets. These filters are therefore a convenient means by which to harvest total microbial communities from multiple aqueous samples in parallel for subsequent molecular analyses.

The genome sequences of the marine diatom Epithemia pelagica strain UHM3201 (Schvarcz, Stancheva & Steward, 2022) and its nitrogen-fixing, endosymbiotic cyanobacterium.

We present the genome assembly of the pennate diatom Epithemia pelagica strain UHM3201 (Ochrophyta; Bacillariophyceae; Rhopalodiales; Rhopalodiaceae) and that of its cyanobacterial endosymbiont (Chroococcales: Aphanothecaceae). The genome sequence of the diatom is 60.3 megabases in span, and the cyanobacterial genome has a length of 2.48 megabases. Most of the diatom nuclear genome assembly is scaffolded into 15 chromosomal pseudomolecules. The organelle genomes have also been assembled, with the mitochondrial genome 40.08 kilobases and the plastid genome 130.75 kilobases in length. A number of other prokaryote MAGs were also assembled.

Broad phylogenetic and functional diversity among mixotrophic consumers of Prochlorococcus.

Small eukaryotic phytoplankton are major contributors to global primary production and marine biogeochemical cycles. Many taxa are thought to be mixotrophic, but quantitative studies of phagotrophy exist for very few. In addition, little is known about consumers of Prochlorococcus, the abundant cyanobacterium at the base of oligotrophic ocean food webs. Here we describe thirty-nine new phytoplankton isolates from the North Pacific Subtropical Gyre (Station ALOHA), all flagellates ~2-5 µm diameter, and we quantify their ability to graze Prochlorococcus. The mixotrophs are from diverse classes (dictyochophytes, haptophytes, chrysophytes, bolidophytes, a dinoflagellate, and a chlorarachniophyte), many from previously uncultured clades. Grazing ability varied substantially, with specific clearance rate (volume cleared per body volume) varying over ten-fold across isolates and six-fold across genera. Slower grazers tended to create more biovolume per prey biovolume consumed. Using qPCR we found that the haptophyte Chrysochromulina was most abundant among the isolated mixotrophs at Station ALOHA, with 76-250 cells mL-1 across depths in the upper euphotic zone (5-100 m). Our results show that within a single ecosystem the phototrophs that ingest bacteria come from many branches of the eukaryotic tree, and are functionally diverse, indicating a broad range of strategies along the spectrum from phototrophy to phagotrophy.

Refining real-time predictions of Vibrio vulnificus concentrations in a tropical urban estuary by incorporating dissolved organic matter dynamics.

The south shore of O'ahu, Hawai'i is one of the most visited coastal tourism areas in the United States with some of the highest instances of recreational waterborne disease. A population of the pathogenic bacterium Vibrio vulnificus lives in the estuarine Ala Wai Canal in Honolulu which surrounds the heavily populated tourism center of Waikiki. We developed a statistical model to predict V. vulnificus dynamics in this system using environmental measurements from moored oceanographic and atmospheric sensors in real time. During a year-long investigation, we analyzed water from 9 sampling events at 3 depths and 8 sites along the canal (n = 213) for 36 biogeochemical variables and V. vulnificus concentration using quantitative polymerase chain reaction (qPCR) of the hemolysin A gene (vvhA). The best multiple linear regression model of V. vulnificus concentration, explaining 80% of variation, included only six predictors: 5-day average rainfall preceding water sampling, daily maximum air temperature, water temperature, nitrate plus nitrite, and two metrics of humic dissolved organic matter (DOM). We show how real-time predictions of V. vulnificus concentration can be made using these models applied to the time series of water quality measurements from the Pacific Islands Ocean Observing System (PacIOOS) as well as the PacIOOS plume model based on the Waikiki Regional Ocean Modeling System (ROMS) products. These applications highlight the importance of including DOM variables in predictive modeling of V. vulnificus and the influence of rain events in elevating nearshore concentrations of V. vulnificus. Long-term climate model projections of locally downscaled monthly rainfall and air temperature were used to predict an overall increase in V. vulnificus concentration of approximately 2- to 3-fold by 2100. Improving these predictive models of microbial populations is critical for management of waterborne pathogen risk exposure, particularly in the wake of a changing global climate.

Overlooked and widespread pennate diatom-diazotroph symbioses in the sea.

Persistent nitrogen depletion in sunlit open ocean waters provides a favorable ecological niche for nitrogen-fixing (diazotrophic) cyanobacteria, some of which associate symbiotically with eukaryotic algae. All known marine examples of these symbioses have involved either centric diatom or haptophyte hosts. We report here the discovery and characterization of two distinct marine pennate diatom-diazotroph symbioses, which until now had only been observed in freshwater environments. Rhopalodiaceae diatoms Epithemia pelagica sp. nov. and Epithemia catenata sp. nov. were isolated repeatedly from the subtropical North Pacific Ocean, and analysis of sequence libraries reveals a global distribution. These symbioses likely escaped attention because the endosymbionts lack fluorescent photopigments, have nifH gene sequences similar to those of free-living unicellular cyanobacteria, and are lost in nitrogen-replete medium. Marine Rhopalodiaceae-diazotroph symbioses are a previously overlooked but widespread source of bioavailable nitrogen in marine habitats and provide new, easily cultured model organisms for the study of organelle evolution.

Towards an integrative view of virus phenotypes.

Understanding how phenotypes emerge from genotypes is a foundational goal in biology. As challenging as this task is when considering cellular life, it is further complicated in the case of viruses. During replication, a virus as a discrete entity (the virion) disappears and manifests itself as a metabolic amalgam between the virus and the host (the virocell). Identifying traits that unambiguously constitute a virus's phenotype is straightforward for the virion, less so for the virocell. Here, we present a framework for categorizing virus phenotypes that encompasses both virion and virocell stages and considers functional and performance traits of viruses in the context of fitness. Such an integrated view of virus phenotype is necessary for comprehensive interpretation of viral genome sequences and will advance our understanding of viral evolution and ecology.

Temporal and spatial variability in culturable pathogenic Vibrio spp. in Lake Pontchartrain, Louisiana, following hurricanes Katrina and Rita.

We investigated the abundance, distribution, and virulence gene content of Vibrio cholerae, V. parahaemolyticus, and V. vulnificus in the waters of southern Lake Pontchartrain in Louisiana on four occasions from October 2005 to September 2006, using selective cultivation and molecular assays. The three targeted pathogenic vibrios were generally below the detection level in January 2006, when the water was cold (13°C), and most abundant in September 2006, when the lake water was warmest (30°C). The maximum values for these species were higher than reported previously for the lake by severalfold to orders of magnitude. The only variable consistently correlated with total vibrio abundance within a single sampling was distance from shore (P = 0.000). Multiple linear regression of the entire data set revealed that distance from shore, temperature, and turbidity together explained 82.1% of the variability in total vibrio CFU. The log-transformed mean abundance of V. vulnificus CFU in the lake was significantly correlated with temperature (P = 0.014), but not salinity (P = 0.625). Virulence-associated genes of V. cholerae (ctx) and V. parahaemolyticus (trh and tdh) were not detected in any isolates of these species (n = 128 and n = 20, respectively). In contrast, 16S rRNA typing of V. vulnificus (n = 298) revealed the presence of both environmental (type A) and clinical (type B) strains. The percentage of the B-type V. vulnificus was significantly higher in the lake in October 2005 (35.8% of the total) than at other sampling times (P ≤ 0.004), consistent with the view that these strains represent distinct ecotypes.

Analysis of a viral metagenomic library from 200 m depth in Monterey Bay, California constructed by direct shotgun cloning.

BACKGROUND: Viruses have a profound influence on both the ecology and evolution of marine plankton, but the genetic diversity of viral assemblages, particularly those in deeper ocean waters, remains poorly described. Here we report on the construction and analysis of a viral metagenome prepared from below the euphotic zone in a temperate, eutrophic bay of coastal California. METHODS: We purified viruses from approximately one cubic meter of seawater collected from 200 m depth in Monterey Bay, CA. DNA was extracted from the virus fraction, sheared, and cloned with no prior amplification into a plasmid vector and propagated in E. coli to produce the MBv200m library. Random clones were sequenced by the Sanger method. Sequences were assembled then compared to sequences in GenBank and to other viral metagenomic libraries using BLAST analyses. RESULTS: Only 26% of the 881 sequences remaining after assembly had significant (E≤0.001) BLAST hits to sequences in the GenBank nr database, with most being matches to bacteria (15%) and viruses (8%). When BLAST analysis included environmental sequences, 74% of sequences in the MBv200m library had a significant match. Most of these hits (70%) were to microbial metagenome sequences and only 0.7% were to sequences from viral metagenomes. Of the 121 sequences with a significant hit to a known virus, 94% matched bacteriophages (Families Podo-, Sipho-, and Myoviridae) and 6% matched viruses of eukaryotes in the Family Phycodnaviridae (5 sequences) or the Mimivirus (2 sequences). The largest percentages of hits to viral genes of known function were to those involved in DNA modification (25%) or structural genes (17%). Based on reciprocal BLAST analyses, the MBv200m library appeared to be most similar to viral metagenomes from two other bays and least similar to a viral metagenome from the Arctic Ocean. CONCLUSIONS: Direct cloning of DNA from diverse marine viruses was feasible and resulted in a distribution of virus types and functional genes at depth that differed in detail, but were broadly similar to those found in surface marine waters. Targeted viral analyses are useful for identifying those components of the greater marine metagenome that circulate in the subcellular size fraction.

Morphological characterization of viruses in the stratified water column of alkaline, hypersaline Mono Lake.

Concentrations of viruses and prokaryotes in the alkaline, moderately hypersaline, seasonally stratified Mono Lake are among the highest reported for a natural aquatic environment. We used electron microscopy to test whether viral morphological characteristics differed among the epilimnion, metalimnion, and the anoxic hypolimnion of the lake and to determine how the properties of viruses in Mono Lake compare to other aquatic environments. Viral capsid size distributions were more similar in the metalimnion and hypolimnion of Mono Lake, while viral tail lengths were more similar in the epilimnion and metalimnion. The percentage of tailed viruses decreased with depth and the relative percentages of tailed phage families changed with depth. The presence of large (>125 nm capsid), untailed viruses in the metalimnion and hypolimnion suggests that eukaryotic viruses are produced in these suboxic and anoxic, hypersaline environments. Capsid diameters of viruses were larger on average in Mono Lake compared to other aquatic environments, and no lemon-shaped or filamentous viruses were found, in contrast to other high-salinity or high-altitude lakes and seas. Our data suggest that the physically and chemically distinct layers of Mono Lake harbor different viral assemblages, and that these assemblages are distinct from other aquatic environments that have been studied. Furthermore, we found that filtration of a sample through a 0.22-µm pore-size filter significantly altered the distribution of viral capsid diameters and tail lengths, resulting in a relative depletion of viruses having larger capsids and longer tails. This observation highlights the potential for bias in molecular surveys of viral diversity, which typically rely on filtration through 0.2- or 0.22-µm pore-size membrane filters to remove bacteria during sample preparation.

A Functional K+ Channel from Tetraselmis Virus 1, a Member of the Mimiviridae.

Potassium ion (K+) channels have been observed in diverse viruses that infect eukaryotic marine and freshwater algae. However, experimental evidence for functional K+ channels among these alga-infecting viruses has thus far been restricted to members of the family Phycodnaviridae, which are large, double-stranded DNA viruses within the phylum Nucleocytoviricota. Recent sequencing projects revealed that alga-infecting members of Mimiviridae, another family within this phylum, may also contain genes encoding K+ channels. Here we examine the structural features and the functional properties of putative K+ channels from four cultivated members of Mimiviridae. While all four proteins contain variations of the conserved selectivity filter sequence of K+ channels, structural prediction algorithms suggest that only two of them have the required number and position of two transmembrane domains that are present in all K+ channels. After in vitro translation and reconstitution of the four proteins in planar lipid bilayers, we confirmed that one of them, a 79 amino acid protein from the virus Tetraselmis virus 1 (TetV-1), forms a functional ion channel with a distinct selectivity for K+ over Na+ and a sensitivity to Ba2+. Thus, virus-encoded K+ channels are not limited to Phycodnaviridae but also occur in the members of Mimiviridae. The large sequence diversity among the viral K+ channels implies multiple events of lateral gene transfer.

Minimum Information about an Uncultivated Virus Genome (MIUViG).

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

A giant virus infecting green algae encodes key fermentation genes.

The family Mimiviridae contains uncommonly large viruses, many of which were isolated using a free-living amoeba as a host. Although the genomes of these and other mimivirids that infect marine heterokont and haptophyte protists have now been sequenced, there has yet to be a genomic investigation of a mimivirid that infects a member of the Viridiplantae lineage (green algae and land plants). Here we characterize the 668-kilobase complete genome of TetV-1, a mimivirid that infects the cosmopolitan green alga Tetraselmis (Chlorodendrophyceae). The analysis revealed genes not previously seen in viruses, such as the mannitol metabolism enzyme mannitol 1-phosphate dehydrogenase, the saccharide degradation enzyme alpha-galactosidase, and the key fermentation genes pyruvate formate-lyase and pyruvate formate-lyase activating enzyme. The TetV genome is the largest sequenced to date for a virus that infects a photosynthetic organism, and its genes reveal unprecedented mechanisms by which viruses manipulate their host's metabolism.

Variables influencing the efficiency and interpretation of reverse transcription quantitative PCR (RT-qPCR): An empirical study using Bacteriophage MS2.

Reverse transcription, quantitative PCR (RT-qPCR) is a sensitive method for quantification of specific RNA targets, but the first step of the assay, reverse transcription, is notoriously variable and sensitive to reaction conditions. In this study, we used purified Bacteriophage MS2 genomic RNA as a model virus target to test two different RT enzymes (SuperScript II and SuperScript III), two RT-priming strategies (gene-specific primers and random hexamers), and varying background RNA concentrations (0-50ngµl-1) to determine how these variables influence the efficiency of reverse transcription over a range of target concentrations (101-107 copies µl-1). The efficiency of the RT reaction was greatly improved by increasing both background RNA and primer concentrations, but the benefit provided by background RNA was source dependent. At a given target concentration, similar RT efficiencies were achieved with gene-specific primers and random hexamers, but the latter required much higher concentrations. With random hexamers, we observed a systematic variation in RT reaction efficiency as a function of target concentration. Using an RNA standard curve that was also subject to RT effectively normalized for this systematic variability, but the assay accuracy depended critically on the length of the standard RNA extending to the 3' end of the qPCR target site. Our results shed some light on previous contradictory conclusions in the literature, and provide insights that may aid in the design of RT-qPCR assays and the design of synthetic RNA standards when full-length material is not available.

Viruses in the Oceanic Basement.

Microbial life has been detected well into the igneous crust of the seafloor (i.e., the oceanic basement), but there have been no reports confirming the presence of viruses in this habitat. To detect and characterize an ocean basement virome, geothermally heated fluid samples (ca. 60 to 65°C) were collected from 117 to 292 m deep into the ocean basement using seafloor observatories installed in two boreholes (Integrated Ocean Drilling Program [IODP] U1362A and U1362B) drilled in the eastern sediment-covered flank of the Juan de Fuca Ridge. Concentrations of virus-like particles in the fluid samples were on the order of 0.2 × 105 to 2 × 105 ml-1 (n = 8), higher than prokaryote-like cells in the same samples by a factor of 9 on average (range, 1.5 to 27). Electron microscopy revealed diverse viral morphotypes similar to those of viruses known to infect bacteria and thermophilic archaea. An analysis of virus-like sequences in basement microbial metagenomes suggests that those from archaeon-infecting viruses were the most common (63 to 80%). Complete genomes of a putative archaeon-infecting virus and a prophage within an archaeal scaffold were identified among the assembled sequences, and sequence analysis suggests that they represent lineages divergent from known thermophilic viruses. Of the clustered regularly interspaced short palindromic repeat (CRISPR)-containing scaffolds in the metagenomes for which a taxonomy could be inferred (163 out of 737), 51 to 55% appeared to be archaeal and 45 to 49% appeared to be bacterial. These results imply that the warmed, highly altered fluids in deeply buried ocean basement harbor a distinct assemblage of novel viruses, including many that infect archaea, and that these viruses are active participants in the ecology of the basement microbiome.IMPORTANCE The hydrothermally active ocean basement is voluminous and likely provided conditions critical to the origins of life, but the microbiology of this vast habitat is not well understood. Viruses in particular, although integral to the origins, evolution, and ecology of all life on earth, have never been documented in basement fluids. This report provides the first estimate of free virus particles (virions) within fluids circulating through the extrusive basalt of the seafloor and describes the morphological and genetic signatures of basement viruses. These data push the known geographical limits of the virosphere deep into the ocean basement and point to a wealth of novel viral diversity, exploration of which could shed light on the early evolution of viruses.