Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mBio
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mBio
    • Editor in Chief
    • Board of Editors
    • AAM Fellows
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article

Bacteriophage Distributions and Temporal Variability in the Ocean’s Interior

Elaine Luo, Frank O. Aylward, Daniel R. Mende, Edward F. DeLong
Mary Ann Moran, Editor
Elaine Luo
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Frank O. Aylward
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Daniel R. Mende
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Edward F. DeLong
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Edward F. DeLong
Mary Ann Moran
University of Georgia
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Eric Allen
Scripps Institution of Oceanography
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander Culley
Université Laval
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mBio.01903-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Bacteriophages are numerically the most abundant DNA-containing entities in the oligotrophic ocean, yet how specific phage populations vary over time and space remains to be fully explored. Here, we conducted a metagenomic time-series survey of double-stranded DNA phages throughout the water column in the North Pacific Subtropical Gyre, encompassing 1.5 years from depths of 25 to 1,000 m. Viral gene sequences were identified in assembled metagenomic samples, yielding an estimated 172,385 different viral gene families. Viral marker gene distributions suggested that lysogeny was more prevalent at mesopelagic depths than in surface waters, consistent with prior prophage induction studies using mitomycin C. A total of 129 ALOHA viral genomes and genome fragments from 20 to 108 kbp were selected for further study, which represented the most abundant phages in the water column. Phage genotypes displayed discrete population structures. Most phages persisted throughout the time-series and displayed a strong depth structure that mirrored the stratified depth distributions of co-occurring bacterial taxa in the water column. Mesopelagic phages were distinct from surface water phages with respect to diversity, gene content, putative life histories, and temporal persistence, reflecting depth-dependent differences in host genomic architectures and phage reproductive strategies. The spatiotemporal distributions of the most abundant open-ocean bacteriophages that we report here provide new insight into viral temporal persistence, life history, and virus-host-environment interactions throughout the open-ocean water column.

IMPORTANCE The North Pacific Subtropical Gyre represents one of the largest biomes on the planet, where microbial communities are central mediators of ecosystem dynamics and global biogeochemical cycles. Critical members of these communities are the viruses of marine bacteria, which can alter microbial metabolism and significantly influence their survival and productivity. To better understand these viral assemblages, we conducted genomic analyses of planktonic viruses over a seasonal cycle to ocean depths of 1,000 m. We identified 172,385 different viral gene families and 129 unique virus genotypes in this open-ocean setting. The spatiotemporal distributions of the most abundant open-ocean viruses that we report here provide new insights into viral temporal variability, life history, and virus-host-environment interactions throughout the water column.

INTRODUCTION

Viruses are abundant biotic entities that play critical roles in aquatic environments. Some of the most common among these viruses in the open ocean are double-stranded DNA (dsDNA) bacteriophages (phages) that infect many abundant and biogeochemically important groups of bacterioplankton, such as Prochlorococcus, Synechococcus, and numerous heterotrophic bacterial species in common genera such as Roseobacter, Alteromonas, Pelagibacter, and Puniceispirillum (1–6). Phages have been shown to kill hosts at rates of up to 20 to 40% of the total population per day, potentially strongly impacting bacterioplankton populations (7, 8). In addition, carbon flux through phage biomass is estimated at 145 gigatonnes per year, playing a substantial role in the global carbon cycle (9). Furthermore, phages can influence ocean biochemistry through microbial cell lysis leading to the production of dissolved organic matter (DOM) and via auxiliary metabolic genes (AMGs) that alter cellular carbon, sulfur, and nitrogen metabolism of their hosts during the phage replication cycle (10–14). Advancing fundamental knowledge of marine phages is therefore an important step toward developing a deeper understanding of marine ecosystem dynamics.

Phages have critical roles in the microbial ecology and biogeochemistry of the global ocean due to their tremendous abundance and diversity. While the genotypic diversity of marine phages has historically been difficult to ascertain, recent studies have provided new insights into viral genomic diversity in the oceans. Developments in high-throughput DNA sequencing have enabled the exploration of viral diversity in the environment at unprecedented scales (13–17). These frequent reports of large reservoirs of viral genetic diversity highlight the importance of further work of in situ characterization of novel marine phages using reference-independent metagenomic techniques that do not rely on a limited subset of cultivable hosts.

The majority of published marine viral metagenomic surveys to date have focused on cataloging the genomic diversity and geographic variability in surface water samples (14–16). The vertical and temporal distributions of environmental phage assemblages have received relatively less attention. To our knowledge, only two metagenomic studies have reported on phages recovered from deep-sea planktonic samples. The Pacific Ocean Virome data set included 12 samples from the deep Pacific revealing that aphotic zone viromes contained a unique set of AMGs that distinguish them from photic zone viromes (13), while the structure of 99 genomic fragments of bathypelagic phages has been reported from the Mediterranean Sea (17). These studies suggested that deep-ocean phages are largely novel and distinct from previously characterized surface phages, highlighting the need to explore the vast diversity of uncharacterized phages below the surface ocean. With respect to temporal variability, previous studies have focused on daily, weekly, or annual scales in surface waters (18–24). Additional metagenomic studies of depth profiles coupled with temporal variability should prove useful as well to provide further insight into the spatiotemporal dynamics of marine phage populations.

Elucidation of viral dynamics in space and time at well-defined study sites has potential to provide the environmental context for interpreting broader patterns and consequences of viral diversity. Here, we present a metagenomic depth profile time-series of phage genes, genomes, and genome fragments captured in cellular bacterioplankton fractions (>0.2 µm) from depths of 25 to 1,000 m over 1.5 years. We used two approaches to explore how phages vary through depth and time in the North Pacific Subtropical Gyre (NPSG), an oligotrophic habitat that is representative of the largest biome in the world (25). Using a genome-centric approach, we describe genomic fragments of abundant phage populations at Station ALOHA. Through implementation of a multistep reassembly workflow, we reconstructed viral population genomes to describe how the diversity, distribution, and genetic repertoire of phages vary through depth and time. In the second approach, we used a gene-centric methodology that leveraged a previous study of prokaryotic assemblages at Station ALOHA (26). Specifically, we used the previously reported nonredundant gene catalog constructed from Station ALOHA to analyze the vertical distribution of phage genes and examine how the diversity and functional repertoire of phages vary across depth profiles. We also used marker genes to explore how viral life-history strategies shift through the ocean’s water column. Our analyses characterize viral gene distributions, genotypes, and temporal dynamics across a range of depths and provide important insight into the genomic diversity and dynamics of viral assemblages in the open ocean.

RESULTS AND DISCUSSION

Viral genotypic diversity was inferred from metagenomic data collected at Station ALOHA from 25 to 1,000 m over 1.5 years, using both genome-centric and gene-centric approaches (see Fig. S2 in the supplemental material). The genome-centric methodologies used generated 129 viral contigs between 20 and 108 kbp in length, representing genomes or large genomic fragments of some of the most abundant phages at Station ALOHA. These were used to characterize the distributions of dominant viral populations and provide genomic insights into the ecology of specific phage groupings. In addition, the gene-centric approach of all Station ALOHA contigs assembled (26) captured a wide range of viral diversity from 177,713 nonredundant viral genes, which facilitated broader quantitative analyses of phage gene distributions in space and time, providing insight into the relationships between viral life history, environmental, and host variability.

FIG S1 

Our data set includes 83 samples from 7 depths and 12 time points from August 2010 to December 2011. Download FIG S1, PDF file, 0.5 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG S2 

Bioinformatic workflow from sequencing to gene- and genome-centric approaches. Download FIG S2, PDF file, 0.1 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

Novel ALOHA viral genomes, genome fragments, and AMGs.Our conservative viral genome assembly strategy yielded many novel viral populations different from previously characterized viruses (Fig. 1). Seventy-nine out of 129 ALOHA viral contigs shared relatively little sequence homology to any known phages or previously sequenced viral metagenomes. Of the 50 Station ALOHA viral contigs having some database homologues, 10 were related to known phage genomes in RefSeq75 (all cyanophages, Table S2), while 40 shared homology to phages in one of three previously sequenced viral metagenomes (14, 15, 27). These data suggested that at least some of the most abundant phage groups that we found were widespread across ocean basins.

TABLE S1 

Assembly statistics for the initial individual assemblies for each sample, number of contigs >3 kb in length used for VirSorter runs, and number of VirSorter identified contigs in all categories. Download TABLE S1, PDF file, 0.1 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

TABLE S2 

List of 10 out of 129 ALOHA viral contigs with hits to known phages in RefSeq75 database with more than half of its genes hitting an average amino acid identity of >60%. Download TABLE S2, PDF file, 0.1 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG 1 
  • Open in new tab
  • Download powerpoint
FIG 1 

Depth profiles of known and novel phage ALOHA viral contigs assigned using four reference databases: known phages in RefSeq75 protein database and three previously sequenced viral metagenomes. Phages are identified with a cutoff of >50% genes with >60% average amino acid identity to a reference database genome. (a) Proportion of contigs with hits to known phages in RefSeq75, other viromes, and novel phages. (b) Subset of hits to RefSeq75. (c) Subset of hits to other viromes: uvMED from the Mediterranean (15), Global Ocean Viromes (GOV) from Tara Oceans (14), and Earth viromes (EV) from human, terrestrial, and marine environments (27). The number of contigs in each category is shown in the legend. Proportion is normalized by total nucleotides mapped to each contig. Error bars show standard errors, which are summed among groups in stacked bars.

Mesopelagic phages appeared to be undersampled in current databases and were distinct from surface water phages. None of the ALOHA viral contigs from 770 m or deeper were similar to any previously reported phages in existing databases. Gene mapping revealed that some ALOHA viral contigs shared conserved genomic structure with known phages in reference databases, despite low amino acid similarity (Fig. S3a to e), supporting the existence of evolutionary constraints on gene order even among distantly related phages.

FIG S3 

Genome maps of abundant phages found at different depths in the water column at Station ALOHA. Download FIG S3, PDF file, 0.3 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

Among viral contigs of all lengths, we found 37 genes identified from 625 contigs colocated with phage structural genes that have functions not previously identified in marine phages (Table S3). These included 19 putative AMGs that provide insight into how phages can manipulate host metabolism. Some of these genes have putative functions in antibiotic synthesis (carbamoyltransferase C terminus, myoinositol-1-phosphate synthase), antimicrobial resistance (dolichyl-phosphate-mannose protein mannosyltransferase), antitoxin synthesis (antitoxin of toxin-antitoxin system), antigen synthesis (P83/100), transporters (sodium bile acid symporter), and superinfection immunity. These new phage-associated genes further expand our current knowledge of gene content in naturally occurring phages. Moreover, certain AMGs were found only in ALOHA viral contigs dominant at specific depths, such as myoinositol-1-phosphate synthase (25 and 75 m), dolichyl-phosphate-mannose protein mannosyltransferase (200 m), sodium bile acid symporter (200 m), P83/100 antigen proteins (200 m), and superinfection immunity protein (sporadic, 770 and 1,000 m). The functional roles of these depth-specific AMGs provide new insight into phage-host interactions in the open-ocean water column.

TABLE S3 

List of 37 novel genes in marine phages identified from 625 curated contigs colocated with viral structural genes (PFAM bit score of >10) that have not been reported in other previously sequenced viromes (Brum et al., 2015; Roux et al., 2016). Download TABLE S3, PDF file, 0.1 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

Phage genotype distributions in the Station ALOHA time-series depth profile.Most ALOHA viral contigs reached peak abundances at a single depth and could be broadly categorized into one of five groups based on abundance profiles: a surface group dominating 25 to 75 m, a deep-chlorophyll-maximum (DCM) group at 125 m, a 200-m group, a deeper mesopelagic group from 500 to 1,000 m, and a sporadic group of more temporally variable phages (Fig. 2). No evidence for eurybathic phages was evident in our data. Instead, the vertical distribution of phage contigs appeared to reflect the depth-stratified distributions of their potential bacterial and archaeal hosts (26). Overall, these results suggest that many dominant phage groups at Station ALOHA may have relatively narrow host ranges.

FIG 2 
  • Open in new tab
  • Download powerpoint
FIG 2 

Coverage profiles of 129 ALOHA viral contigs through depth and time. Each node on the top dendrogram and its associated column represent one contig, while each row represents one of 83 total samples. The height of the black column within a cell represents the mean coverage (second and third quartiles) of a contig in a given sample. The full bar height of each sample layer is log scaled to the maximum coverage in a given sample. Sample layers are ordered by depth from 25 to 1,000 m, shown on the right. Within each depth bin, samples are ordered by time from August 2010 to December 2011, shown on the left. Contigs are clustered by differential coverage, with manual binning into five groups on the bottom. Across the top, homology to known phages in the RefSeq75 protein database and other viral metagenomes was assigned based on >50% genes hitting at >60% average amino acid identity.

This depth structure in ALOHA viral contig distributions was also reflected in our gene-centric approach using the ALOHA gene catalog (26). Comparison of sample clustering patterns based on phage- versus bacterioplankton-specific genes indicated that viral gene distributions mirrored community-wide trends of their potential bacterioplankton hosts. Overall patterns of depth-specific gene content were highly similar between phage-specific and cellular genes (Fig. 3), reflecting depth clustering patterns in the genome-centric approach. Phage- and cell-associated gene distribution dendrograms were conserved across these depth clusters over time. This suggested that the overall pattern of gene distributions was conserved between phages and their cellular host community.

FIG 3 
  • Open in new tab
  • Download powerpoint
FIG 3 

Cluster analysis of phage-specific versus cellular genes identified in the Station ALOHA time-series nonredundant gene catalog. Each edge on the trees represents one sample, and corresponding samples are shown with connecting lines between the two dendrograms.

The pronounced spatial differences along the depth gradient were accompanied by depth-stratified differences in temporal variability of ALOHA viral contigs, with persistent phages dominating surface waters and more episodic phages dominating mesopelagic depths. Most phages displayed no clear trends of seasonality or shuffling of dominant phage groups, with the exception of a small group of phages in the sporadic group (Fig. 2). In the persistent groups (surface, DCM, 200 m, and mesopelagic), phage populations displayed remarkable temporal stability throughout the 1.5-year sampling period, similar to that previously reported for surface water viruses over shorter daily time scales (24). Additionally, four phages in our persistent group, captured in the 2010–2011 Hawaii Ocean Time-series (HOT) data set, shared high sequence similarity to phages found in a 2015 diel 15-m phage study in the NPSG (24), indicating that these populations were consistently present over multiple years (Table S4). In contrast to the patterns observed in surface waters, phages in the mesopelagic ocean exhibited a more-sporadic occurrence characteristic of boom-and-bust dynamics. Several ALOHA viral contigs were highly abundant at only a few time points and virtually absent at other times, as shown by the relative abundances of phages at 25 and 1,000 m (Fig. 4a). To confirm that the greater temporal variability at 1,000 m is not attributed to the lower number of viral contigs detected, mean-normalized variances show that mesopelagic phages were indeed more temporally variable than surface water phages (Fig. 4b).

TABLE S4 

List of 4 persistent ALOHA viral contigs with hits to surface phages in a 2015 data set near Station ALOHA (Aylward et al., https://doi.org/10.1073/pnas.1714821114 ) with more than half of its genes hitting an average amino acid identity of >60%. Download TABLE S4, PDF file, 0.1 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG 4 
  • Open in new tab
  • Download powerpoint
FIG 4 

Temporal variability of assembled viral contigs increased with depth. (a) Relative abundance of the most abundant ALOHA viral contigs at 25 and 1,000 m. Each row represents one of the 13 most abundant ALOHA viral contigs for each depth bin, with size indicated in kilobase pairs (kbp). Asterisks represent contigs containing one or more prophage markers (PFAM bit score of >30 to integrase, CI repressor, Cro, or excisionase). Relative abundance is scaled by area. (b) Mean-normalized variance of relative abundance through time calculated for each contig.

Genomic trends in lytic versus lysogenic viral life history.The phage depth distributions were also reflected in encoded viral life-history traits. Using prophage marker genes as indicators, we found that the genomic potential for lysogeny increased below the DCM, which hovers around 90 to 130 m throughout the year (25). We focus on three known prophage markers from studies on bacteriophage lambda as proxies for lysogeny, though it is worth noting that this is not the only means of lysogenic conversion. It should be noted that some putative prophage markers also exist in nonphage genetic elements, for example, the integrases associated with integrative conjugative elements. In our analyses, all three prophage markers (integrase, repressor protein cI [CI repressor], and excisionase) showed similar patterns of proportional increase with depth below 125 m (Fig. 5a to c). Collectively, the data suggested that the proportion of phages that are capable of lysogeny was ~5 times higher in the mesopelagic than in the surface ocean. Furthermore, the average number of prophage markers per cellular genome increased below 125 m, suggesting that the average number of integrated prophages per cell also increased with depth (Fig. S4a to c). The copy number of integrases appears to be anomalously high relative to other prophage and phage markers, suggesting that there may be representation from non-phage mobile genomic elements or that these gene families are more highly conserved and therefore more easily detected using homology-based methods. Despite the possibility of false positives in using integrase as a prophage marker, the overall increased abundance of other prophage markers provides evidence for an increased prevalence of lysogeny in the mesopelagic ocean. Lastly, the average copy number of 4 phage markers per cellular genome ranged from 1 to 5 at 125 m while values were 0.45, 0.31, 0.52, and 0.04 at 1,000 m for DNA polymerase, terminase, capsid, and tail, respectively (Fig. S4d to g). The strikingly similar decreases across all four phage markers suggest that the prevalence of actively replicating lytic phages within cells decreased greatly at mesopelagic depths. The potential inclusion of free phages adhered to particles in our samples is unlikely to significantly impact these results due to prefiltration of larger, >1.6-μm particles before bacterioplankton collection. Taken together, these marker gene data suggest a shift of increased lysogeny and overall decreased viral particle abundance per host, from the euphotic to the mesopelagic zone.

FIG S4 

Depth profile of phage marker proteins identified in the Station ALOHA nonredundant gene catalog (domain bit score of >50) using hidden Markov models generated with a manually curated set of proteins from NCBI. Download FIG S4, PDF file, 0.2 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

FIG 5 
  • Open in new tab
  • Download powerpoint
FIG 5 

Depth profile of prophage marker proteins identified in the Station ALOHA time-series nonredundant gene catalog (domain bit score of >50) using hidden Markov models generated with manually curated sets of viral marker proteins from NCBI. Each circle represents a sample mean, and each vertical bar represents a depth mean. Depth profile of prophage markers integrase (a), CI repressor (b), and excisionase (c). The fold increase of prophage markers with respect to surface mean is calculated using coverage of marker proteins normalized to coverage of capsid proteins.

Viral life-history strategies are important to consider for both phage and host ecology. For the host, a lytic cycle results in cell death and the release of cellular material into the environment, while a lysogenic cycle does not immediately kill the host but incurs a cost of carrying and reproducing foreign genetic material in the host genome. For the phage, a lytic cycle means a rapid increase in short-term fitness when a host is productive enough to support phage production, while lysogeny may represent an opportunity cost in reproduction but increased chances of survival. Characterizing which environments favor a certain viral life-history strategy is important when considering virus-host interactions and their resulting biogeochemical effects along the ocean’s interior.

Different theories exist to explain the ecologic factors that may influence the proportion of lytic versus lysogenic phages. For example, the “piggy-back the winner” hypothesis predicts that lysogeny should be more prevalent at higher host cell densities (28) and has been based on cell- and virus-like particle abundance patterns (9, 28, 29). Our metagenomic data do not directly support these results, since we found evidence for more lysogens in deeper waters that have lower cell densities (cell abundance data available at http://hahana.soest.hawaii.edu/hot/hot-dogs/ ), consistent with earlier results from mitomycin C phage induction studies (30).

An alternative hypothesis is that lysogeny will be favored in environments where there are fewer constraints on cellular genome size (31, 32). In surface oligotrophic oceans, bacteria tend to have smaller genomes, which has been hypothesized to be a consequence of low-nutrient conditions leading to metabolic constraints that favor smaller genomes, lower G+C contents, and decreased nitrogen content in peptides (25, 33). At Station ALOHA, bacterial genome size nearly doubles just below the DCM in the “genomic transition zone” (26), where we also observed a marked increase in prophage marker genes. This increase in genome size below the DCM is likely more permissive to foreign DNA accumulation, thereby enabling replicative elements such as prophages to persist in the more capacious bacterial genomes that comprise deeper-water microbial communities.

Lysogeny might also be advantageous in the deep ocean where low host productivity cannot support rapid phage production (34–36). The results reported here are consistent with a prior prophage induction experiment that showed greater potential for lysogeny in deep than in surface waters in the Mediterranean and Baltic Seas (30). Other field observations have found increased lysogeny in low-productivity environments such as oligotrophic rather than coastal oceans (37, 38) or in winter rather than spring (23, 39). Although this trend is not consistent across all studies (e.g., references 40 and 41), productivity appears to be a major correlate of phage-host interactions in some environments. At Station ALOHA, productivity declines sharply past 125 m near the DCM (25), coinciding with the sharp transition from lytic to lysogenic phage-host interaction in our marker gene profiles. Our study provides a genomic perspective and validation for induction experiments, as well as molecular evidence for a higher proportion of lysogens and decreased viral particles per host in the deep ocean. Our results also suggest that phage-induced mortality may be higher in more productive surface waters, shifting to more temperate phage-host interactions in the mesopelagic open ocean.

Ecology of surface and mesopelagic phages.Overall, the most abundant surface phage populations were well represented in the 1.6- to 0.2-μm size fraction throughout our sampling period (Fig. 2 and 4). This temporal persistence is intriguing, since dominant phages are expected to follow boom-and-bust dynamics in a negative density-dependent manner, according to some ecological models such as kill-the-winner or fluctuating selection dynamics model (reviewed in reference 42). On the other hand, phage resistance mechanisms allow coastal marine Synechococcus populations to coexist with abundant cyanophages over seasonal cycles (18), and persistent viral types have been observed in other viral time-series analyses (19–22). Moreover, the stability of phages in surface waters may reflect the overall stability of prokaryotic assemblages in the NSPG (43). Alternatively, phages could display increased temporal persistence due to a broader host range or longer viral residency time. A combination of these factors might contribute to the temporal stability of many different host-phage pairs, particularly those persistent in ocean surface waters.

The most temporally variable phage groups in our analyses occurred in the mesopelagic ocean, where specific phages dominate at certain time points only to disappear the next month (Fig. 2 and 4). This sporadic nature corresponds to high temporal variability in known phage and bacterial genes observed in the mesopelagic ocean, particularly at 1,000 m (Fig. S5a and b). Compared to the surface ocean, where sunlight drives consistent high productivity, productivity in the deep ocean is limited by sporadic rain of organic material exported from the surface (25). Temporally variable resources might select for copiotrophic or particle-attached bacteria (44) that can grow quickly when resources become available. These bursts of growth can lead to temporally sporadic phage distributions from prophage propagation with fast-growing hosts, prophage induction, or particularly successful lytic infections in terminating host blooms.

FIG S5 

Bubble plot of proportion of ALOHA gene catalog genes hitting four groups of known phages (a) and associated bacteria (b). Download FIG S5, PDF file, 0.3 MB.
Copyright © 2017 Luo et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license .

In this environment of temporally variable resources, we found one novel gene of interest that is specific to a mesopelagic sporadic phage: a superinfection immunity gene encoding a membrane-attached protein that confers host immunity to other phages (45) that has not been previously observed in marine phages. The associated contig (Fig. S3e) in our sampling period followed the distribution of Vibrio (Fig. S5b), copiotrophic bacteria that have been found on sinking particles and in deep samples at Station ALOHA (46). This putative lytic vibriophage is remarkably abundant in an environment where lysogeny is favored. By conferring host immunity to other competitor phages, this gene could contribute to the success of this abundant putative lytic phage and result in subsequent boom-and-bust dynamics of this phage and its particle-attached host in the mesopelagic ocean. It is worth noting here that such genomic characterization of phage populations in situ may be prone to generating false positives in functional capacity (47). Experimental verification of the function of carried phage genes will more reliably elucidate the genomic capacity of previously uncharacterized mesopelagic phages and the resulting phage-host interaction.

Conclusions.In summary, our time-series study in the NPSG leveraged both gene- and genome-centric approaches to provide insight into phage diversity, structure, and function across depth and time. We found that mesopelagic phages were distinct from surface phages and were largely novel and underrepresented in current viral reference and metagenomic databases. With respect to depth variability, discrete phage populations displayed strong depth structure, similar to that of putative bacterial hosts. There were virtually no eurybathic phages, suggesting that most dominant phage groups were adapted to relatively narrow host ranges. We also found unique AMGs in a mesopelagic phage suggestive of depth-specific adaptations to a more variable landscape of phage-host interactions in the aphotic zone. With respect to temporal variability, the most abundant phage groups were remarkably persistent, displaying little to no seasonality nor observable shuffling of dominant groups. With respect to coupled variability through space and time, mesopelagic phages were more sporadic in their distribution. Considering viral life history, we found five times more genes associated with lysogeny at greater depths and a sharp increase in lysogeny at and below the DCM, corresponding to the cell-associated genomic transition zone (26). Our observations, in addition to other recent studies (13–17), expand the realm of characterized phages in the surface to the mesopelagic ocean.

MATERIALS AND METHODS

Study site and sample collection.Bacterioplankton samples in the 0.2- to 1.6-µm size fraction were collected from 7 depths (25, 75, 125, 200, 500, 770, and 1,000 m) on 12 occasions in 2010–2011 at Station ALOHA (22°45′N, 158°W) in the North Pacific Subtropical Gyre (NPSG). Detailed sample collection has been previously described (26). As the study site of the Hawaii Ocean Time-series (HOT) program, Station ALOHA is one of the most sampled open-ocean systems in the world, with well-characterized biogeochemical gradients that provide context to our work (25) (examples in Fig. S1 in the supplemental material).

Station ALOHA metagenomic assembly and gene catalog.The methods for DNA extraction, library construction, and sequencing with the Illumina NextSeq and MiSeq platforms of these bacterioplankton samples have been previously described in detail (26). Briefly, metagenomic reads from each sample were assembled using MIRA v4.9.5_2, providing the basis for our genome-centric analyses. Forty million genes were predicted from all assemblies using Prodigal v2.6.0 (48) and clustered using CD-HIT v2.6 (49) to generate 8,966,703 nonredundant gene clusters. These nonredundant genes are here referred to as the ALOHA gene catalog, providing the basis for our gene-centric analyses. The relative abundance of each nonredundant gene in each sample was calculated by mapping reads using the BWA-MEM algorithm v0.7.15 (50) with default parameters and then dividing the resulting coverage by the total coverage of all genes (26).

Genome-centric analyses. (i) Assembling ALOHA viral contigs.Station ALOHA metagenomic assemblies were used to identify viral contigs (Fig. S2) using VirSorter v1.0.3 with a 3-kbp cutoff for improved recall (51). Sequence reads that were used to assemble these contigs were pooled across all samples and reassembled using default parameters in metaSPAdes (52). The 104,732 contigs were validated through a second VirSorter screen, in which 917 viral contigs 1.7 to 108 kbp in length from all VirSorter categories were retained. Most cellular contamination was removed at this stage, as the estimated cellular genomic completeness was only 1.2%, with only 35 out of 16,296 genes mapping to single-copy prokaryotic marker genes using the standard Anvi’o v.2.1.1 workflow (53–55). Of these 35 hits, 30 were restricted to recA and DNA helicase, which could be recombination and replication homologous phage genes. One contig found to contain a ribosomal bacterial marker gene was removed from our analyses.

Given that 20 kbp represents the lower end of genome size in DNA phages from marine systems (56), we retained 142 contigs of ≥20 kbp in size to focus on near-complete genomes or large genomic fragments. As a final quality control step, we removed contigs that did not contain any genes with distant homology to known phage structural proteins (PFAM bit score of >10 to terminase, portal, capsid, tail, base plate, spike, neck, and head genes). This step was conservative and eliminated 13 putative phage genomes from our analyses, 5 of which displayed significant homology to previously sequenced phages (identification described below). The resulting contigs represent a high-confidence subset of total phage diversity. We checked again for cellular contamination using cell genomic completion, which is now reduced to 0.6%. Only 12 out of a total of 5,877 putative phage genes mapped to bacterial or archaeal marker genes. Eleven of these hits are restricted to recA and DNA helicase, while none were ribosomal cell marker genes. Overall, these quality control steps successfully removed any detectable cellular genome contaminants, generating 129 ALOHA high-confidence viral contigs ranging between 20 and 108 kbp in length. Assembly and reassembly statistics are shown in Table S1.

(ii) Annotating ALOHA viral contigs.Contigs were annotated with a combined viral database from predicted proteins in assembled sequences from phages in RefSeq75 (57) and four other viral metagenomes: uvMED (15), uvDeep (17), GOV (14), and EV (27). Genes were predicted using Prodigal v2.6.3 (48). Contigs were identified using LAST (58) to this combined viral database with cutoffs of >50% of genes hit to a reference with >60% average amino acid identity. This relaxed cutoff is chosen to identify distant homology to limited reference phages based on a majority of genes on a given contig. To calculate the proportion of genomes assigned to each category in a given sample, contigs were normalized by their relative abundance of nucleotides mapped.

(iii) Depth and temporal distributions of ALOHA viral contigs.For reference-independent visualization of spatiotemporal distributions, we used BWA-MEM v0.7.15 (50) to map all reads from each sample to ALOHA viral contigs and generated coverage profiles using Anvi’o v.2.1.1 (55). To prevent the possibility of conserved phage genes inflating coverage, only the second and third coverage value quartiles across each contig were used to generate mean coverage profiles. ALOHA viral contigs were binned by differential coverage manually into five distribution groups. To examine temporal variability within each depth bin, we used the proportion of nucleotides mapped as a metric for relative abundance and plotted the 13 most abundant phages through time within each depth bin. To examine whether temporal variability in phage assemblages changes with depth, we calculated the variance of the relative abundance for each ALOHA viral contig within all depth bins through time. Given that more abundant contigs tend to have larger variances, we normalized the temporal variance of each contig by dividing by its mean abundance through time. This mean-normalized variance alleviates preferential weighting to high-abundance contigs and prevents changes in number of contigs among depths from influencing the variance metric.

(iv) Mapping to reference genomes.To examine the genomic structure of our ALOHA viral contigs, we generated genome maps of some of the most abundant or complete ALOHA viral contigs in each of five depth distribution groups. Phage genomes in the RefSeq75 database that were most similar to each ALOHA viral contig were included to visualize structural similarities using the GenomeDiagram module v0.2 in Python (59).

(v) Detection of novel phage genes and AMGs.To detect AMGs, we predicted and annotated 16,286 genes from 917 reassembled viral contigs (all lengths) and compared these genes to existing databases of marine phage genes including AMGs (13, 14). Genes were considered novel AMGs if they were both annotated with a PFAM bit score of >30 to a novel function not in these databases and colocated on the subset of 625 contigs with viral structural genes (PFAM bit score of >10 to terminase, portal, capsid, tail, base plate, spike, neck, and head).

Gene-centric analyses. (i) Phage- and cell-associated gene catalog assemblages.We clustered samples using the ALOHA gene catalog to examine how closely patterns of spatiotemporal diversity of phage assemblages mirrored those of broader bacterioplankton communities. We separated genes into two groups of interest: a phage group of 177,713 viral genes and a cell-associated group encompassing all remaining genes in the ALOHA gene catalog (8,788,990 in total). Viral genes were identified using the combined viral database described above with LAST at >90% amino acid identity. A total of 5,328 photosystem genes (PFAM bit score of >30) were removed from the phage group on the grounds that phage copies of this AMG (reviewed in reference 60) could not be accurately distinguished from bacterium-encoded copies (3). This approach identified 172,385 phage genes used in subsequent analyses. Independently for phage- and cell-associated groups, we clustered samples using gene coverages to generate Bray-Curtis distance matrices and subsequently average linkage hierarchical agglomerative clustering on the R vegan package (61). The resulting dendrograms were visualized using the R dendextend package (62).

To visualize the depth structure and distribution of known phages and bacteria, we identified 180,055 genes from the ALOHA gene catalog with >60% amino acid identity to known phages and 3,295,154 bacterial genes with >60% amino acid identity to known bacteria in RefSeq75. Using these gene coverages, we calculated relative abundances of three groups of well-characterized phages and their associated bacterial hosts for each sample.

(ii) Viral life history.To examine whether the prevalence of lysogeny changes with depth, we identified prophage markers in the ALOHA gene catalog and used their relative abundance as a proxy for the potential incidence of lysogeny. We also analyzed other sets of well-known phage structural proteins for normalization. We curated sets of annotated prophage and phage markers of interest from NCBI and identified 89 prophage integrases, 116 prophage CI repressors, 357 prophage excisionases, 1,780 phage DNA polymerases, 139 phage terminases, 971 phage capsids, and 104 phage tail fiber proteins. Curated proteins were selected from either full phage genomes or uncultured marine phages from previously sequenced viral metagenomes. We used these protein sets to generate hidden Markov models for each marker using MUSCLE v3.8.31 alignment and HMMER v3.1 (63, 64). We then used these models to identify marker proteins in the translated ALOHA gene catalog with a domain bit score of >50 across the whole sequence. The summed coverages of prophage markers (integrases, CI repressors, and excisionases) were normalized to the summed coverage of capsids as a proxy for the proportion of phages that had the potential to reproduce through lysogeny based on the presence of prophage marker genes. We further normalized this proportion by dividing by the surface ocean mean value to identify relative changes in lysogeny with depth.

To examine depth profiles of the number and nature of active infections, we calculated the average cell genomic copy number of phage markers (DNA polymerases, terminases, capsids, and tail fibers) and prophage markers (integrases, CI repressors, and excisionases). Given that prophages in cellular genomes are predominantly nonreplicative, the copy number of prophage markers was used to estimate the number of prophages per genome. The copy number of phage markers was used to estimate the number of total phages captured inside a bacterium. To generate copy number per cellular genome, coverages of prophage and phage markers were normalized to the average coverage of 10 universal single-copy bacterial marker genes, called mOTU profiling (26, 65).

Accession number(s).Sequence data are available at NCBI under BioProject no. PRJNA352737 and at https://www.imicrobe.us/#/projects/263 .

ACKNOWLEDGMENTS

We thank the captain and crew of the R/V Kilo Moana and the Hawaii Ocean Time-series marine operations team for sample collection and oceanographic data acquisition. We thank Tsultrim Palden and Anna Romano for DNA library preparation and sequencing. We thank John Eppley and Torben Nielson for generating initial assemblies and assistance with the ALOHA gene catalog. We thank Elisha Wood-Charlson, Jessica Bryant, Daniel Olson, and Murat Eren for expert advice in data analyses.

This research was supported by grants from the Simons Foundation to E.F.D. (SCOPE 329108) and the Gordon and Betty Moore Foundation to E.F.D. (GBMF 3777). Partial support for D.R.M. was provided by the European Molecular Biology Organization (ALTF 721-2015), and partial support for E.L. was provided by the Natural Sciences and Engineering Research Council of Canada (PGSD3-487490-2016). This work is a contribution of the Simons Collaboration on Ocean Processes and Ecology and the Center for Microbial Oceanography: Research and Education. The funders had no role in study design, sample collection, data analyses, or the decision to submit the work for publication.

FOOTNOTES

    • Received 17 October 2017
    • Accepted 25 October 2017
    • Published 28 November 2017
  • Copyright © 2017 Luo et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license .

REFERENCES

  1. 1.↵
    1. Wilson WH,
    2. Joint IR,
    3. Carr NG,
    4. Mann NH
    . 1993. Isolation and molecular characterization of five marine cyanophages propagated on Synechococcus sp. strain WH7803. Appl Environ Microbiol59:3736–3743.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Rohwer F,
    2. Segall A,
    3. Steward G,
    4. Seguritan V,
    5. Breitbart M,
    6. Wolven F,
    7. Azam F
    . 2000. The complete genomic sequence of the marine phage Roseophage SIO1 shares homology with nonmarine phages. Limnol Oceanogr45:408–418. doi:10.4319/lo.2000.45.2.0408.
    OpenUrlCrossRef
  3. 3.↵
    1. Lindell D,
    2. Sullivan MB,
    3. Johnson ZI,
    4. Tolonen AC,
    5. Rohwer F,
    6. Chisholm SW
    . 2004. Transfer of photosynthesis genes to and from Prochlorococcus viruses. Proc Natl Acad Sci U S A101:11013–11018. doi:10.1073/pnas.0401526101.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    1. Zhao Y,
    2. Temperton B,
    3. Thrash JC,
    4. Schwalbach MS,
    5. Vergin KL,
    6. Landry ZC,
    7. Ellisman M,
    8. Deerinck T,
    9. Sullivan MB,
    10. Giovannoni SJ
    . 2013. Abundant SAR11 viruses in the ocean. Nature494:357–360. doi:10.1038/nature11921.
    OpenUrlCrossRefPubMedWeb of Science
  5. 5.↵
    1. Kang I,
    2. Oh HM,
    3. Kang D,
    4. Cho JC
    . 2013. Genome of a SAR116 bacteriophage shows the prevalence of this phage type in the oceans. Proc Natl Acad Sci U S A110:12343–12348. doi:10.1073/pnas.1219930110.
    OpenUrlAbstract/FREE Full Text
  6. 6.↵
    1. Garcia-Heredia I,
    2. Rodriguez-Valera F,
    3. Martin-Cuadrado AB
    . 2013. Novel group of podovirus infecting the marine bacterium Alteromonas macleodii. Bacteriophage3:e24766. doi:10.4161/bact.24766.
    OpenUrlCrossRef
  7. 7.↵
    1. Weinbauer MG
    . 2004. Ecology of prokaryotic viruses. FEMS Microbiol Rev28:127–181. doi:10.1016/j.femsre.2003.08.001.
    OpenUrlCrossRefPubMedWeb of Science
  8. 8.↵
    1. Weinbauer MG,
    2. Rassoulzadegan F
    . 2004. Are viruses driving microbial diversification and diversity?Environ Microbiol6:1–11. doi:10.1046/j.1462-2920.2003.00539.x.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    1. Lara E,
    2. Vaqué D,
    3. Sà EL,
    4. Boras JA,
    5. Gomes A,
    6. Borrull E,
    7. Díez-Vives C,
    8. Teira E,
    9. Pernice MC,
    10. Garcia FC,
    11. Forn I,
    12. Castillo YM,
    13. Peiró A,
    14. Salazar F,
    15. Morán XAG,
    16. Massana R,
    17. Catalá TS,
    18. Luna GM,
    19. Agustí S,
    20. Estrada M,
    21. Gasol JM,
    22. Duarte CM
    . 2017. Unveiling the role and life strategies of viruses from the surface to the dark ocean. Sci Adv3:e1602565. doi:10.1126/sciadv.1602565.
    OpenUrlFREE Full Text
  10. 10.↵
    1. Lindell D,
    2. Jaffe JD,
    3. Johnson ZI,
    4. Church GM,
    5. Chisholm SW
    . 2005. Photosynthesis genes in marine viruses yield proteins during host infection. Nature438:86–89. doi:10.1038/nature04111.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Thompson LR,
    2. Zeng Q,
    3. Kelly L,
    4. Huang KH,
    5. Singer AU,
    6. Stubbe J,
    7. Chisholm SW
    . 2011. Phage auxiliary metabolic genes and the redirection of cyanobacterial host carbon metabolism. Proc Natl Acad Sci U S A108:E757–E764. doi:10.1073/pnas.1102164108.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Breitbart M
    . 2012. Marine viruses: truth or dare. Annu Rev Mar Sci4:425–448. doi:10.1146/annurev-marine-120709-142805.
    OpenUrlCrossRefPubMed
  13. 13.↵
    1. Hurwitz BL,
    2. Sullivan MB
    . 2013. The Pacific Ocean Virome (POV): a marine viral metagenomic dataset and associated protein clusters for quantitative viral ecology. PLoS One8:e57355. doi:10.1371/journal.pone.0057355.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Roux S,
    2. Brum JR,
    3. Dutilh BE,
    4. Sunagawa S,
    5. Duhaime MB,
    6. Loy A,
    7. Poulos BT,
    8. Solonenko N,
    9. Lara E,
    10. Poulain J,
    11. Pesant S,
    12. Kandels-Lewis S,
    13. Dimier C,
    14. Picheral M,
    15. Searson S,
    16. Cruaud C,
    17. Alberti A,
    18. Duarte CM,
    19. Gasol JM,
    20. Vaqué D,
    21. Tara Oceans Coordinators
    22. Bork P,
    23. Acinas SG,
    24. Wincker P,
    25. Sullivan MB
    . 2016. Ecogenomics and potential biogeochemical impacts of globally abundant ocean viruses. Nature537:689–693. doi:10.1038/nature19366.
    OpenUrlCrossRefPubMed
  15. 15.↵
    1. Mizuno CM,
    2. Rodriguez-Valera F,
    3. Kimes NE,
    4. Ghai R
    . 2013. Expanding the marine virosphere using metagenomics. PLoS Genet9:e1003987. doi:10.1371/journal.pgen.1003987.
    OpenUrlCrossRef
  16. 16.↵
    1. Brum JR,
    2. Ignacio-Espinoza JC,
    3. Roux S,
    4. Doulcier G,
    5. Acinas SG,
    6. Alberti A,
    7. Cruaud C,
    8. de Vargas C,
    9. Gasol JM,
    10. Gorsky G,
    11. Gregory AC,
    12. Guidi L,
    13. Hingamp P,
    14. Iudicone D,
    15. Not F,
    16. Ogata H,
    17. Pesant S,
    18. Poulos BT,
    19. Schwenck S,
    20. Speich S,
    21. Dimier C,
    22. Kandels-Lewis S,
    23. Picheral M,
    24. Searson S,
    25. Tara Oceans Coordinators
    26. Bork P,
    27. Bowler C,
    28. Sunagawa S,
    29. Wincker P,
    30. Karsenti E,
    31. Sullivan MB
    . 2015. Patterns and ecological drivers of ocean viral communities. Science348:1261498. doi:10.1126/science.1261498.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Mizuno CM,
    2. Ghai R,
    3. Saghaï A,
    4. López-García P,
    5. Rodriguez-Valera F
    . 2016. Genomes of abundant and widespread viruses from the deep ocean. mBio7:e00805-16. doi:10.1128/mBio.00805-16.
    OpenUrlAbstract/FREE Full Text
  18. 18.↵
    1. Waterbury JB,
    2. Valois FW
    . 1993. Resistance to co-occurring phages enables marine Synechococcus to co-exist with cyanophages abundant in seawater. Appl Environ Microbiol59:3393–3399.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Marston MF,
    2. Sallee JL
    . 2003. Genetic diversity and temporal variation in the cyanophage community infecting marine Synechococcus species in Rhode Island’s coastal waters. Appl Environ Microbiol69:4639–4647. doi:10.1128/AEM.69.8.4639-4647.2003.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Pagarete A,
    2. Chow CE,
    3. Johannessen T,
    4. Fuhrman JA,
    5. Thingstad TF,
    6. Sandaa RA
    . 2013. Strong seasonality and interannual recurrence in marine myovirus communities. Appl Environ Microbiol79:6253–6259. doi:10.1128/AEM.01075-13.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Chow CET,
    2. Kim DY,
    3. Sachdeva R,
    4. Caron DA,
    5. Fuhrman JA
    . 2014. Top-down controls on bacterial community structure: microbial network analysis of bacteria, T4-like viruses and protists. ISME J8:816–829. doi:10.1038/ismej.2013.199.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Goldsmith DB,
    2. Parsons RJ,
    3. Beyene D,
    4. Salamon P,
    5. Breitbart M
    . 2015. Deep sequencing of the viral phoH gene reveals temporal variation, depth-specific composition, and persistent dominance of the same viral phoH genes in the Sargasso Sea. PeerJ3:e997. doi:10.7717/peerj.997.
    OpenUrlCrossRef
  23. 23.↵
    1. Brum JR,
    2. Hurwitz BL,
    3. Schofield O,
    4. Ducklow HW,
    5. Sullivan MB
    . 2015. Seasonal time bombs: dominant temperate viruses affect Southern Ocean microbial dynamics. ISME J10:1–13. doi:10.1038/ismej.2015.125.
    OpenUrlCrossRef
  24. 24.↵
    1. Aylward FO,
    2. Boeuf D,
    3. Mende DR,
    4. Wood-Charlson EM,
    5. Vislova A,
    6. Eppley JM,
    7. Romano AE,
    8. DeLong EF
    . 2017. Diel cycling and long-term persistence of viruses in the ocean’s euphotic zone. Proc Natl Acad Sci U S A114:11446–11451. doi:10.1073/pnas.1714821114.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Karl DM,
    2. Church MJ
    . 2014. Microbial oceanography and the Hawaii Ocean Time-series programme. Nat Rev Microbiol12:699–713. doi:10.1038/nrmicro3333.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Mende DR,
    2. Bryant JA,
    3. Aylward FO,
    4. Eppley JM,
    5. Nielsen T,
    6. Karl DM,
    7. DeLong EF
    . 2017. Environmental drivers of a microbial genomic transition zone in the ocean’s interior. Nat Microbiol2:1367–1373. doi:10.1038/s41564-017-0008-3.
    OpenUrlCrossRef
  27. 27.↵
    1. Paez-Espino D,
    2. Eloe-Fadrosh EA,
    3. Pavlopoulos GA,
    4. Thomas AD,
    5. Huntemann M,
    6. Mikhailova N,
    7. Rubin E,
    8. Ivanova NN,
    9. Kyrpides NC
    . 2016. Uncovering Earth’s virome. Nature536:425–430. doi:10.1038/nature19094.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Knowles B,
    2. Silveira CB,
    3. Bailey BA,
    4. Barott K,
    5. Cantu VA,
    6. Cobián-Güemes AG,
    7. Coutinho FH,
    8. Dinsdale EA,
    9. Felts B,
    10. Furby KA,
    11. George EE,
    12. Green KT,
    13. Gregoracci GB,
    14. Haas AF,
    15. Haggerty JM,
    16. Hester ER,
    17. Hisakawa N,
    18. Kelly LW,
    19. Lim YW,
    20. Little M,
    21. Luque A,
    22. McDole-Somera T,
    23. McNair K,
    24. de Oliveira LS,
    25. Quistad SD,
    26. Robinett NL,
    27. Sala E,
    28. Salamon P,
    29. Sanchez SE,
    30. Sandin S,
    31. Silva GGZ,
    32. Smith J,
    33. Sullivan C,
    34. Thompson C,
    35. Vermeij MJA,
    36. Youle M,
    37. Young C,
    38. Zgliczynski B,
    39. Brainard R,
    40. Edwards RA,
    41. Nulton J,
    42. Thompson F,
    43. Rohwer F
    . 2016. Lytic to temperate switching of viral communities. Nature531:466–470. doi:10.1038/nature17193.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Wigington CH,
    2. Sonderegger DL,
    3. Brussaard CPD,
    4. Buchan A,
    5. Finke JF,
    6. Fuhrman J,
    7. Lennon JT,
    8. Middelboe M,
    9. Suttle CA,
    10. Stock C,
    11. Wilson WH,
    12. Wommack KE,
    13. Wilhelm SW,
    14. Weitz JS
    . 2016. Re-examining the relationship between virus and microbial cell abundances in the global oceans. Nat Microbiol1:1–21. doi:10.1038/nmicrobiol.2015.24.
    OpenUrlCrossRef
  30. 30.↵
    1. Weinbauer MG,
    2. Brettar I,
    3. Höfle MG
    . 2003. Lysogeny and virus-induced mortality of bacterioplankton in surface, deep, and anoxic marine waters. Limnol Oceanogr48:1457–1465. doi:10.4319/lo.2003.48.4.1457.
    OpenUrlCrossRef
  31. 31.↵
    1. Casjens S
    . 2003. Prophages and bacterial genomics: what have we learned so far?Mol Microbiol49:277–300. doi:10.1046/j.1365-2958.2003.03580.x.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    1. Touchon M,
    2. Bernheim A,
    3. Rocha EP
    . 2016. Genetic and life-history traits associated with the distribution of prophages in bacteria. ISME J10:2744–2754. doi:10.1038/ismej.2016.47.
    OpenUrlCrossRef
  33. 33.↵
    1. Batut B,
    2. Knibbe C,
    3. Marais G,
    4. Daubin V
    . 2014. Reductive genome evolution at both ends of the bacterial population size spectrum. Nat Rev Microbiol12:841–850. doi:10.1038/nrmicro3331.
    OpenUrlCrossRefPubMed
  34. 34.↵
    1. Stewart FM,
    2. Levin BR
    . 1984. The population biology of bacterial viruses: why be temperate. Theor Popul Biol26:93–117. doi:10.1016/0040-5809(84)90026-1.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Middelboe M
    . 2000. Bacterial growth rate and marine virus-host dynamics. Microb Ecol40:114–124.
    OpenUrlCrossRefPubMedWeb of Science
  36. 36.↵
    1. Moebus K
    . 1996. Marine bacteriophage reproduction under nutrient-limited growth of host bacteria. I. Investigations with six phage-host systems. Mar Ecol Prog Ser144:1–12. doi:10.3354/meps144001.
    OpenUrlCrossRef
  37. 37.↵
    1. Jiang SC,
    2. Paul JH
    . 1998. Significance of lysogeny in the marine environment: studies with isolates and a model of lysogenic phage production. Microb Ecol35:235–243. doi:10.1007/s002489900079.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Weinbauer MG,
    2. Suttle CA
    . 1999. Lysogeny and prophage induction in coastal and offshore bacterial communities. Aquat Microb Ecol18:217–225. doi:10.3354/ame018217.
    OpenUrlCrossRef
  39. 39.↵
    1. Williamson SJ,
    2. Houchin LA,
    3. McDaniel L,
    4. Paul JH
    . 2002. Seasonal variation in lysogeny as depicted by prophage induction in Tampa Bay, Florida. Appl Environ Microbiol68:4307–4314. doi:10.1128/AEM.68.9.4307-4314.2002.
    OpenUrlAbstract/FREE Full Text
  40. 40.↵
    1. Laybourn-Parry J,
    2. Marshall WA,
    3. Madan NJ
    . 2007. Viral dynamics and patterns of lysogeny in saline Antarctic lakes. Polar Biol30:351–358. doi:10.1007/s00300-006-0191-9.
    OpenUrlCrossRefWeb of Science
  41. 41.↵
    1. Payet JP,
    2. Suttle CA
    . 2013. To kill or not to kill: the balance between lytic and lysogenic viral infection is driven by trophic status. Limnol Oceanogr58:465–474. doi:10.4319/lo.2013.58.2.0465.
    OpenUrlCrossRefWeb of Science
  42. 42.↵
    1. Avrani S,
    2. Schwartz DA,
    3. Lindell D
    . 2012. Virus-host swinging party in the oceans: incorporating biological complexity into paradigms of antagonistic coexistence. Mob Genet Elements2:88–95. doi:10.4161/mge.20031.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Bryant JA,
    2. Aylward FO,
    3. Eppley JM,
    4. Karl DM,
    5. Church MJ,
    6. DeLong EF
    . 2016. Wind and sunlight shape microbial diversity in surface waters of the North Pacific Subtropical Gyre. ISME J10:1308–1322. doi:10.1038/ismej.2015.221.
    OpenUrlCrossRef
  44. 44.↵
    1. DeLong EF,
    2. Preston CM,
    3. Mincer T,
    4. Rich V,
    5. Hallam SJ,
    6. Frigaard NU,
    7. Martinez A,
    8. Sullivan MB,
    9. Edwards R,
    10. Brito BR,
    11. Chisholm SW,
    12. Karl DM
    . 2006. Community genomics among stratified microbial assemblages in the ocean’s interior. Science311:496–503. doi:10.1126/science.1120250.
    OpenUrlAbstract/FREE Full Text
  45. 45.↵
    1. Lu MJ,
    2. Henning U
    . 1989. The immunity (imm) gene of Escherichia coli bacteriophage T4. J Virol63:3472–3478.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Fontanez KM,
    2. Eppley JM,
    3. Samo TJ,
    4. Karl DM,
    5. DeLong EF
    . 2015. Microbial community structure and function on sinking particles in the North Pacific Subtropical Gyre. Front Microbiol6:469. doi:10.3389/fmicb.2015.00469.
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. Enault F,
    2. Briet A,
    3. Bouteille L,
    4. Roux S,
    5. Sullivan MB,
    6. Petit MA
    . 2017. Phages rarely encode antibiotic resistance genes: a cautionary tale for virome analyses. ISME J11:237–247. doi:10.1038/ismej.2016.90.
    OpenUrlCrossRef
  48. 48.↵
    1. Hyatt D,
    2. Chen GL,
    3. Locascio PF,
    4. Land ML,
    5. Larimer FW,
    6. Hauser LJ
    . 2010. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics11:119. doi:10.1186/1471-2105-11-119.
    OpenUrlCrossRefPubMed
  49. 49.↵
    1. Li W,
    2. Jaroszewski L,
    3. Godzik A
    . 2001. Clustering of highly homologous sequences to reduce the size of large protein databases. Bioinformatics17:282–283. doi:10.1093/bioinformatics/17.3.282.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Li H
    . 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXivarXiv:1303.3997v2 [q-bio.GN]. arXiv:1303.3997v2.
  51. 51.↵
    1. Roux S,
    2. Enault F,
    3. Hurwitz BL,
    4. Sullivan MB
    . 2015. VirSorter: mining viral signal from microbial genomic data. PeerJ3:e985. doi:10.7717/peerj.985.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Bankevich A,
    2. Nurk S,
    3. Antipov D,
    4. Gurevich AA,
    5. Dvorkin M,
    6. Kulikov AS,
    7. Lesin VM,
    8. Nikolenko SI,
    9. Pham S,
    10. Prjibelski AD,
    11. Pyshkin AV,
    12. Sirotkin AV,
    13. Vyahhi N,
    14. Tesler G,
    15. Alekseyev MA,
    16. Pevzner PA
    . 2012. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol19:455–477. doi:10.1089/cmb.2012.0021.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Rinke C,
    2. Schwientek P,
    3. Sczyrba A,
    4. Ivanova NN,
    5. Anderson IJ,
    6. Cheng JF,
    7. Darling A,
    8. Malfatti S,
    9. Swan BK,
    10. Gies EA,
    11. Dodsworth JA,
    12. Hedlund BP,
    13. Tsiamis G,
    14. Sievert SM,
    15. Liu WT,
    16. Eisen JA,
    17. Hallam SJ,
    18. Kyrpides NC,
    19. Stepanauskas R,
    20. Rubin EM,
    21. Hugenholtz P,
    22. Woyke T
    . 2013. Insights into the phylogeny and coding potential of microbial dark matter. Nature499:431–437. doi:10.1038/nature12352.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Campbell JH,
    2. O’Donoghue P,
    3. Campbell AG,
    4. Schwientek P,
    5. Sczyrba A,
    6. Woyke T,
    7. Söll D,
    8. Podar M
    . 2013. UGA is an additional glycine codon in uncultured SR1 bacteria from the human microbiota. Proc Natl Acad Sci U S A110:5540–5545. doi:10.1073/pnas.1303090110.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Eren AM,
    2. Esen ÖC,
    3. Quince C,
    4. Vineis JH,
    5. Morrison HG,
    6. Sogin ML,
    7. Delmont TO
    . 2015. Anvi’o: an advanced analysis and visualization platform for ’omics data. PeerJ3:e1319. doi:10.7717/peerj.1319.
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. Steward GF,
    2. Montiel JL,
    3. Azam F
    . 2000. Genome size distributions indicate variability and similarities among marine viral assemblages from diverse environments. Limnol Oceanogr45:1697–1706. doi:10.4319/lo.2000.45.8.1697.
    OpenUrlCrossRef
  57. 57.↵
    1. O’Leary NA,
    2. Wright MW,
    3. Brister JR,
    4. Ciufo S,
    5. Haddad D,
    6. McVeigh R,
    7. Rajput B,
    8. Robbertse B,
    9. Smith-White B,
    10. Ako-Adjei D,
    11. Astashyn A,
    12. Badretdin A,
    13. Bao Y,
    14. Blinkova O,
    15. Brover V,
    16. Chetvernin V,
    17. Choi J,
    18. Cox E,
    19. Ermolaeva O,
    20. Farrell CM,
    21. Goldfarb T,
    22. Gupta T,
    23. Haft D,
    24. Hatcher E,
    25. Hlavina W,
    26. Joardar VS,
    27. Kodali VK,
    28. Li W,
    29. Maglott D,
    30. Masterson P,
    31. McGarvey KM,
    32. Murphy MR,
    33. O’Neill K,
    34. Pujar S,
    35. Rangwala SH,
    36. Rausch D,
    37. Riddick LD,
    38. Schoch C,
    39. Shkeda A,
    40. Storz SS
    . 2016. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res44:D733–D745. doi:10.1093/nar/gkv1189.
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. Kiełbasa SM,
    2. Wan R,
    3. Sato K,
    4. Horton P,
    5. Frith MC
    . 2011. Adaptive seeds tame genomic sequence comparison. Genome Res21:487–493. doi:10.1101/gr.113985.110.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Pritchard L,
    2. White JA,
    3. Birch PRJ,
    4. Toth IK
    . 2006. GenomeDiagram: a python package for the visualization of large-scale genomic data. Bioinformatics22:616–617. doi:10.1093/bioinformatics/btk021.
    OpenUrlCrossRefPubMedWeb of Science
  60. 60.↵
    1. Hurwitz BL,
    2. U’Ren JM
    . 2016. Viral metabolic reprogramming in marine ecosystems. Curr Opin Microbiol31:161–168. doi:10.1016/j.mib.2016.04.002.
    OpenUrlCrossRef
  61. 61.↵
    1. Oksanen J,
    2. Blanchet FG,
    3. Friendly M,
    4. Kindt R,
    5. Legendre P,
    6. McGlinn D,
    7. Minchin PR,
    8. O’Hara RB,
    9. Simpson GL,
    10. Solymos P,
    11. Stevens MHH,
    12. Szoecs E,
    13. Wagner H
    . 2016. vegan: Community Ecology Package. R package version 2.4-0. https://CRAN.R-project.org/package=vegan .
  62. 62.↵
    1. Galili T
    . 2015. dendextend: an R package for visualizing, adjusting, and comparing trees of hierarchical clustering. Bioinformatics31:3718–3720. doi:10.1093/bioinformatics/btv428.
    OpenUrlCrossRefPubMed
  63. 63.↵
    1. Edgar RC
    . 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res32:1792–1797. doi:10.1093/nar/gkh340.
    OpenUrlCrossRefPubMedWeb of Science
  64. 64.↵
    1. Eddy SR
    . 2011. Accelerated profile HMM searches. PLoS Comput Biol7:e1002195. doi:10.1371/journal.pcbi.1002195.
    OpenUrlCrossRefPubMed
  65. 65.↵
    1. Sunagawa S,
    2. Mende DR,
    3. Zeller G,
    4. Izquierdo-Carrasco F,
    5. Berger SA,
    6. Kultima JR,
    7. Coelho LP,
    8. Arumugam M,
    9. Tap J,
    10. Nielsen HB,
    11. Rasmussen S,
    12. Brunak S,
    13. Pedersen O,
    14. Guarner F,
    15. de Vos WM,
    16. Wang J,
    17. Li J,
    18. Doré J,
    19. Ehrlich SD,
    20. Stamatakis A,
    21. Bork P
    . 2013. Metagenomic species profiling using universal phylogenetic marker genes. Nat Methods10:1196–1199. doi:10.1038/nmeth.2693.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Bacteriophage Distributions and Temporal Variability in the Ocean’s Interior
Elaine Luo, Frank O. Aylward, Daniel R. Mende, Edward F. DeLong
mBio Nov 2017, 8 (6) e01903-17; DOI: 10.1128/mBio.01903-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mBio article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Bacteriophage Distributions and Temporal Variability in the Ocean’s Interior
(Your Name) has forwarded a page to you from mBio
(Your Name) thought you would be interested in this article in mBio.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Bacteriophage Distributions and Temporal Variability in the Ocean’s Interior
Elaine Luo, Frank O. Aylward, Daniel R. Mende, Edward F. DeLong
mBio Nov 2017, 8 (6) e01903-17; DOI: 10.1128/mBio.01903-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS AND DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

bacteriophages
biodiversity
Oceans and Seas
Seawater
bacteriophages
bacterioplankton
marine microbiology
microbial ecology
microbial oceanography

Related Articles

Cited By...

About

  • About mBio
  • Editor in Chief
  • Board of Editors
  • AAM Fellows
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Article Types
  • Ethics
  • Contact Us

Follow #mBio

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2150-7511