Locality and diel cycling of viral production revealed by a 24 h time course cross-omics analysis in a coastal region of Japan

Viruses under the Antarctic Ice Shelf are active and potentially involved in global nutrient cycles

Viruses play an important role in the marine ecosystem. However, our comprehension of viruses inhabiting the dark ocean, and in particular, under the Antarctic Ice Shelves, remains limited. Here, we mine single-cell genomic, transcriptomic, and metagenomic data to uncover the viral diversity, biogeography, activity, and their role as metabolic facilitators of microbes beneath the Ross Ice Shelf. This is the largest Antarctic ice shelf with a major impact on global carbon cycle. The viral community found in the cavity under the ice shelf mainly comprises endemic viruses adapted to polar and mesopelagic environments. The low abundance of genes related to lysogenic lifestyle (<3%) does not support a predominance of the Piggyback-the-Winner hypothesis, consistent with a low-productivity habitat. Our results indicate a viral community actively infecting key ammonium and sulfur-oxidizing chemolithoautotrophs (e.g. Nitrosopumilus spp, Thioglobus spp.), supporting a "kill-the-winner" dynamic. Based on genome analysis, these viruses carry specific auxiliary metabolic genes potentially involved in nitrogen, sulfur, and phosphorus acquisition. Altogether, the viruses under Antarctic ice shelves are putatively involved in programming the metabolism of ecologically relevant microbes that maintain primary production in these chemosynthetically-driven ecosystems, which have a major role in global nutrient cycles.

Virus-Host Interactions Drive Contrasting Bacterial Diel Dynamics in the Ocean

Marine organisms perform a sea of diel rhythmicity. Planktonic diel dynamics have been shown to be driven by light, energy resources, circadian rhythms, and the coordinated coupling of photoautotrophs and heterotrophic bacterioplankton. Here, we explore the diel fluctuation of viral production and decay and their impact on the total and active bacterial community in the coastal and open seawaters of the South China Sea. The results showed that the night-production diel pattern of lytic viral production was concurrent with the lower viral decay at night, contributing to the accumulation of the viral population size during the night for surface waters. The diel variations in bacterial activity, community composition, and diversity were found highly affected by viral dynamics. This was revealed by the finding that bacterial community diversity was positively correlated to lytic viral production in the euphotic zone of the open ocean but was negatively related to lysogenic viral production in the coastal ocean. Such distinct but contrasting correlations suggest that viral life strategies can not only contribute to diversifying bacterial community but also potentially piggyback their host to dominate bacterial community, suggesting the tightly synchronized depth-dependent and habitat-specific diel patterns of virus-host interactions. It further implies that viruses serve as an ecologically important driver of bacterial diel dynamics across the ocean, highlighting the viral roles in bacterial ecological and biogeochemical processes in the ocean.

Benzo[a]pyrene stress impacts adaptive strategies and ecological functions of earthworm intestinal viromes

The earthworm gut virome influences the structure and function of the gut microbiome, which in turn influences worm health and ecological functions. However, despite its ecological and soil quality implications, it remains elusive how earthworm intestinal phages respond to different environmental stress, such as soil pollution. Here we used metagenomics and metatranscriptomics to investigate interactions between the worm intestinal phages and their bacteria under different benzo[a]pyrene (BaP) concentrations. Low-level BaP (0.1 mg kg^-1) stress stimulated microbial metabolism (1.74-fold to control), and enhanced the antiphage defense system (n = 75) against infection (8 phage-host pairs). Low-level BaP exposure resulted in the highest proportion of lysogenic phages (88%), and prophages expressed auxiliary metabolic genes (AMGs) associated with nutrient transformation (e.g., amino acid metabolism). In contrast, high-level BaP exposure (200 mg kg^-1) disrupted microbial metabolism and suppressed the antiphage systems (n = 29), leading to the increase in phage-bacterium association (37 phage-host pairs) and conversion of lysogenic to lytic phages (lysogenic ratio declined to 43%). Despite fluctuating phage-bacterium interactions, phage-encoded AMGs related to microbial antioxidant and pollutant degradation were enriched, apparently to alleviate pollution stress. Overall, these findings expand our knowledge of complex phage-bacterium interactions in pollution-stressed worm guts, and deepen our understanding of the ecological and evolutionary roles of phages.

Trade-Offs between Competitive Ability and Resistance to Top-Down Control in Marine Microbes

Trade-offs between competitive ability and resistance to top-down control manifest the "kill-the-winner" hypothesis that explains how mortality caused by protists and viruses can promote bacterial diversity. However, the existence of such trade-offs has rarely been investigated in natural marine bacterial communities. To address this question, we conducted on-board dilution experiments to manipulate top-down control pressure (protists only or protists plus viruses [protists+viruses] combined) and then applied 16S rRNA gene high-throughput sequencing techniques to assess the responses of each bacterial taxon. Dilution experiments enabled us to measure the top-down-control-free growth rate as the competitive ability and top-down-control-caused mortality as the reverse of resistance to top-down control. Overall, bacterial taxa with higher top-down-control-free growth rates were accompanied by lower top-down-control-caused resistance. Furthermore, competition-resistance trade-offs were stronger and more consistent when top-down control was caused by protists+viruses combined than by protists only. When protists+viruses were diluted, the bacterial rank abundance distribution became steepened and evenness and richness were decreased. However, when protists were diluted, only richness decreased. Our results indicate the existence of competition-resistance trade-offs in marine microbes and demonstrate the positive impacts of such trade-offs on bacterial diversity. Regardless, the strength of the competition-resistance trade-offs and the impacts on bacterial diversity were contingent on whether top-down control was caused by protists+viruses combined or protists only. IMPORTANCE We addressed the "kill-the-winner" hypothesis from the perspective of its principle (the competition-resistance trade-off) in marine bacterial communities incubated in situ. Our results supported the existence of competition-resistance trade-offs and the positive effect on bacterial community diversity. The study linked theoretical expectations and complex natural systems and provided new knowledge regarding how top-down controls and competition trade-offs shaped natural bacterial communities.

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

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

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

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

Viral Complexity

Although traditionally viewed as streamlined and simple, discoveries over the last century have revealed that viruses can exhibit surprisingly complex physical structures, genomic organization, ecological interactions, and evolutionary histories. Viruses can have physical dimensions and genome lengths that exceed many cellular lineages, and their infection strategies can involve a remarkable level of physiological remodeling of their host cells. Virus-virus communication and widespread forms of hyperparasitism have been shown to be common in the virosphere, demonstrating that dynamic ecological interactions often shape their success. And the evolutionary histories of viruses are often fraught with complexities, with chimeric genomes including genes derived from numerous distinct sources or evolved de novo. Here we will discuss many aspects of this viral complexity, with particular emphasis on large DNA viruses, and provide an outlook for future research.

Year-round dynamics of amplicon sequence variant communities differ among eukaryotes, Imitervirales, and prokaryotes in a coastal ecosystem

Coastal microbial communities are affected by seasonal environmental change, biotic interactions, and fluctuating nutrient availability. We investigated the seasonal dynamics of communities of eukaryotes, a major group of double-stranded DNA viruses that infect eukaryotes (order Imitervirales; phylum Nucleocytoviricota), and prokaryotes in the Uranouchi Inlet, Kochi, Japan. We performed metabarcoding using ribosomal RNA genes and viral polB genes as markers in 43 seawater samples collected over 20 months. Eukaryotes, prokaryotes, and Imitervirales communities characterized by the compositions of amplicon sequence variants (ASVs) showed synchronic seasonal cycles. However, the community dynamics showed intriguing differences in several aspects, such as the recovery rate after a year. We also showed that the differences in community dynamics were at least partially explained by differences in recurrence/persistence levels of individual ASVs among eukaryotes, prokaryotes, and Imitervirales. Prokaryotic ASVs were the most persistent, followed by eukaryotic ASVs and Imitervirales ASVs, which were the least persistent. We argue that the differences in the specificity of interactions (virus-eukaryote vs. prokaryote-eukaryote) as well as the niche breadth of community members were at the origin of the distinct community dynamics among eukaryotes, their viruses, and prokaryotes.

Spatio-temporal dynamics of bacterial communities in the shoreline of Laurentian great Lake Erie and Lake St. Clair's large freshwater ecosystems

Background

Long-term trends in freshwater bacterial community composition (BCC) and dynamics are not yet well characterized, particularly in large lake ecosystems. We addressed this gap by temporally (15 months) and spatially (6 sampling locations) characterizing BCC variation in lakes Erie and St. Clair; two connected ecosystems in the Laurentian Great Lakes.

Results

We found a spatial variation of the BCC between the two lakes and among the sampling locations (significant changes in the relative abundance of 16% of the identified OTUs at the sampling location level). We observed five distinct temporal clusters (UPGMA broad-scale temporal variation) corresponding to seasonal variation over the 15 months of sampling. Temporal variation among months was high, with significant variation in the relative abundance of 69% of the OTUs. We identified significant differences in taxonomic composition between summer months of 2016 and 2017, with a corresponding significant reduction in the diversity of BCC in summer 2017.

Conclusions

As bacteria play a key role in biogeochemical cycling, and hence in healthy ecosystem function our study defines the scope for temporal and spatial variation in large lake ecosystems. Our data also show that freshwater BCC could serve as an effective proxy and monitoring tool to access large lake health.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02306-y.

High Transcriptional Activity and Diverse Functional Repertoires of Hundreds of Giant Viruses in a Coastal Marine System

Viruses belonging to the Nucleocytoviricota phylum are globally distributed and include members with notably large genomes and complex functional repertoires. Recent studies have shown that these viruses are particularly diverse and abundant in marine systems, but the magnitude of actively replicating Nucleocytoviricota present in ocean habitats remains unclear. In this study, we compiled a curated database of 2,431 Nucleocytoviricota genomes and used it to examine the gene expression of these viruses in a 2.5-day metatranscriptomic time-series from surface waters of the California Current. We identified 145 viral genomes with high levels of gene expression, including 90 Imitervirales and 49 Algavirales viruses. In addition to recovering high expression of core genes involved in information processing that are commonly expressed during viral infection, we also identified transcripts of diverse viral metabolic genes from pathways such as glycolysis, the TCA cycle, and the pentose phosphate pathway, suggesting that virus-mediated reprogramming of central carbon metabolism is common in oceanic surface waters. Surprisingly, we also identified viral transcripts with homology to actin, myosin, and kinesin domains, suggesting that viruses may use these gene products to manipulate host cytoskeletal dynamics during infection. We performed phylogenetic analysis on the virus-encoded myosin and kinesin proteins, which demonstrated that most belong to deep-branching viral clades, but that others appear to have been acquired from eukaryotes more recently. Our results highlight a remarkable diversity of active Nucleocytoviricota in a coastal marine system and underscore the complex functional repertoires expressed by these viruses during infection. IMPORTANCE The discovery of giant viruses has transformed our understanding of viral complexity. Although viruses have traditionally been viewed as filterable infectious agents that lack metabolism, giant viruses can reach sizes rivalling cellular lineages and possess genomes encoding central metabolic processes. Recent studies have shown that giant viruses are widespread in aquatic systems, but the activity of these viruses and the extent to which they reprogram host physiology in situ remains unclear. Here, we show that numerous giant viruses consistently express central metabolic enzymes in a coastal marine system, including components of glycolysis, the TCA cycle, and other pathways involved in nutrient homeostasis. Moreover, we found expression of several viral-encoded actin, myosin, and kinesin genes, indicating viral manipulation of the host cytoskeleton during infection. Our study reveals a high activity of giant viruses in a coastal marine system and indicates they are a diverse and underappreciated component of microbial diversity in the ocean.

Seasonal and diel patterns of abundance and activity of viruses in the Red Sea

Virus-microbe interactions have been studied in great molecular details for many years in cultured model systems, yielding a plethora of knowledge on how viruses use and manipulate host machinery. Since the advent of molecular techniques and high-throughput sequencing, methods such as cooccurrence, nucleotide composition, and other statistical frameworks have been widely used to infer virus-microbe interactions, overcoming the limitations of culturing methods. However, their accuracy and relevance is still debatable as cooccurrence does not necessarily mean interaction. Here we introduce an ecological perspective of marine viral communities and potential interaction with their hosts, using analyses that make no prior assumptions on specific virus-host pairs. By size fractionating water samples into free viruses and microbes (i.e., also viruses inside or attached to their hosts) and looking at how viral group abundance changes over time along both fractions, we show that the viral community is undergoing a change in rank abundance across seasons, suggesting a seasonal succession of viruses in the Red Sea. We use abundance patterns in the different size fractions to classify viral clusters, indicating potential diverse interactions with their hosts and potential differences in life history traits between major viral groups. Finally, we show hourly resolved variations of intracellular abundance of similar viral groups, which might indicate differences in their infection cycles or metabolic capacities.

Predetermined clockwork microbial worlds: Current understanding of aquatic microbial diel response from model systems to complex environments

In the photic zone of aquatic ecosystems, microorganisms with different metabolisms and their viruses form complex interactions and food webs. Within these interactions, phototrophic microorganisms such as eukaryotic microalgae and cyanobacteria interact directly with sunlight, and thereby generate circadian rhythms. Diel cycling originally generated in microbial phototrophs is directly transmitted toward heterotrophic microorganisms utilizing the photosynthetic products as they are excreted or exuded. Such diel cycling seems to be indirectly propagated toward heterotrophs as a result of complex biotic interactions. For example, cell death of phototrophic microorganisms induced by viral lysis and protistan grazing provides additional resources of dissolved organic matter to the microbial community, and so generates diel cycling in other heterotrophs with different nutrient dependencies. Likewise, differences in the diel transmitting pathway via complex interactions among heterotrophs, and between heterotrophs and their viruses, may also generate higher variation and time lag diel rhythms in different heterotrophic taxa. Thus, sunlight and photosynthesis not only contribute energy and carbon supply, but also directly or indirectly control diel cycling of the microbial community through complex interactions in the photic zone of aquatic ecosystems.

A single-cell polony method reveals low levels of infected Prochlorococcus in oligotrophic waters despite high cyanophage abundances

Long-term stability of picocyanobacteria in the open oceans is maintained by a balance between synchronous division and death on daily timescales. Viruses are considered a major source of microbial mortality, however, current methods to measure infection have significant methodological limitations. Here we describe a method that pairs flow-cytometric sorting with a PCR-based polony technique to simultaneously screen thousands of taxonomically resolved individual cells for intracellular virus DNA, enabling sensitive, high-throughput, and direct quantification of infection by different virus lineages. Under controlled conditions with picocyanobacteria-cyanophage models, the method detected infection throughout the lytic cycle and discriminated between varying infection levels. In North Pacific subtropical surface waters, the method revealed that only a small percentage of Prochlorococcus (0.35–1.6%) were infected, predominantly by T4-like cyanophages, and that infection oscillated 2-fold in phase with the diel cycle. This corresponds to 0.35–4.8% of Prochlorococcus mortality daily. Cyanophages were 2–4-fold more abundant than Prochlorococcus, indicating that most encounters did not result in infection and suggesting infection is mitigated via host resistance, reduced phage infectivity and inefficient adsorption. This method will enable quantification of infection for key microbial taxa across oceanic regimes and will help determine the extent that viruses shape microbial communities and ecosystem level processes.

Differential Responses of a Coastal Prokaryotic Community to Phytoplanktonic Organic Matter Derived from Cellular Components and Exudates

The phytoplanktonic production and prokaryotic consumption of organic matter significantly contribute to marine carbon cycling. Organic matter released from phytoplankton via three processes (exudation of living cells, cell disruption through grazing, and viral lysis) shows distinct chemical properties. We herein investigated the effects of phytoplanktonic whole-cell fractions (WF) (representing cell disruption by grazing) and extracellular fractions (EF) (representing exudates) prepared from Heterosigma akashiwo, a bloom-forming Raphidophyceae, on prokaryotic communities using culture-based experiments. We analyzed prokaryotic community changes for two weeks. The shift in cell abundance by both treatments showed similar dynamics, reaching the first peak (~4.1×10⁶‍ ‍cells‍ ‍mL^–1) on day 3 and second peak (~1.1×10⁶‍ ‍cells‍ ‍mL^–1) on day 13. We classified the sequences obtained into operational taxonomic units (OTUs). A Bray-Curtis dissimilarity analysis revealed that the OTU-level community structure changed distinctively with the two treatments. Ten and 13 OTUs were specifically abundant in the WF and EF treatments, respectively. These OTUs were assigned as heterotrophic bacteria mainly belonging to the Alteromonadales (Gammaproteobacteria) and Bacteroidetes clades and showed successive dynamics following the addition of organic matter. We also analyzed the dynamics of these OTUs in the ocean using publicly available metagenomic data from a natural coastal bloom in Monterey Bay, USA. At least two WF treatment OTUs showed co-occurrence with H. akashiwo, indicating that the blooms of H. akashiwo also affect these OTUs in the ocean. The present results strongly suggest that the thriving and dead cells of uninfected phytoplankton differentially influence the marine prokaryotic community.

Temporal transcriptional patterns of cyanophage genes suggest synchronized infection of cyanobacteria in the oceans

BACKGROUND

Based on the peak expression times during infection, early, middle, and late genes have been characterized in viruses (cyanophages) that infect the unicellular cyanobacterium Prochlorococcus. Laboratory experiments show that some cyanophages can only replicate in the light and thus exhibit diurnal infection rhythms under light-dark cycles. Field evidence also suggests synchronized infection of Prochlorococcus by cyanophages in the oceans, which should result in progressive expression of cyanophage early, middle, and late genes. However, distinct temporal expression patterns have not been observed in cyanophage field populations.

RESULTS

In this study, we reanalyzed a previous metatranscriptomic dataset collected in the North Pacific Subtropical Gyre. In this dataset, it was previously shown that aggregate transcripts from cyanophage scaffolds display diurnal transcriptional rhythms with transcript abundances decreasing at night. By mapping metatranscriptomic reads to individual viral genes, we identified periodically expressed genes from putative viruses infecting the cyanobacteria Prochlorococcus and Synechococcus, heterotrophic bacteria, and algae. Of the 41 cyanophage genes, 35 were from cyanomyoviruses. We grouped the periodically expressed cyanomyovirus genes into early, middle, and late genes based on the conserved temporal expression patterns of their orthologs in cyanomyovirus laboratory cultures. We found that the peak expression times of late genes in cyanophage field populations were significantly later than those of early and middle genes, which were similar to the temporal expression patterns of synchronized cyanophage laboratory cultures.

CONCLUSIONS

The significantly later peak expression times of late genes in cyanomyovirus field populations suggest that cyanophage infection of Prochlorococcus is synchronized in the North Pacific Subtropical Gyre. The night-time peak expression of late genes also suggests synchronized lysis of Prochlorococcus at night, which might result in synchronized release of dissolved organic matter to the marine food web. Video abstract.

Tara Oceans: towards global ocean ecosystems biology

A planetary-scale understanding of the ocean ecosystem, particularly in light of climate change, is crucial. Here, we review the work of Tara Oceans, an international, multidisciplinary project to assess the complexity of ocean life across comprehensive taxonomic and spatial scales. Using a modified sailing boat, the team sampled plankton at 210 globally distributed sites at depths down to 1,000 m. We describe publicly available resources of molecular, morphological and environmental data, and discuss how an ecosystems biology approach has expanded our understanding of plankton diversity and ecology in the ocean as a planetary, interconnected ecosystem. These efforts illustrate how global-scale concepts and data can help to integrate biological complexity into models and serve as a baseline for assessing ecosystem changes and the future habitability of our planet in the Anthropocene epoch.

In silico Prediction of Virus-Host Interactions for Marine Bacteroidetes With the Use of Metagenome-Assembled Genomes

Bacteroidetes is one of the most abundant heterotrophic bacterial taxa in the ocean and play crucial roles in recycling phytoplankton-derived organic matter. Viruses of Bacteroidetes are also expected to have an important role in the regulation of host communities. However, knowledge on marine Bacteroidetes viruses is biased toward cultured viruses from a few species, mainly fish pathogens or Bacteroidetes not abundant in marine environments. In this study, we investigated the recently reported 1,811 marine viral genomes to identify putative Bacteroidetes viruses using various in silico host prediction techniques. Notably, we used microbial metagenome-assembled genomes (MAGs) to augment the marine Bacteroidetes reference genomic data. The examined viral genomes and MAGs were derived from simultaneously collected samples. Using nucleotide sequence similarity-based host prediction methods, we detected 31 putative Bacteroidetes viral genomes. The MAG-based method substantially enhanced the predictions (26 viruses) when compared with the method that is solely based on the reference genomes from NCBI RefSeq (7 viruses). Previously unrecognized genus-level groups of Bacteroidetes viruses were detected only by the MAG-based method. We also developed a host prediction method based on the proportion of Bacteroidetes homologs in viral genomes, which detected 321 putative Bacteroidetes virus genomes including 81 that were newly recognized as Bacteroidetes virus genomes. The majority of putative Bacteroidetes viruses were detected based on the proportion of Bacteroidetes homologs in both RefSeq and MAGs; however, some were detected in only one of the two datasets. Putative Bacteroidetes virus lineages included not only relatives of known viruses but also those phylogenetically distant from the cultured viruses, such as marine Far-T4 like viruses known to be widespread in aquatic environments. Our MAG and protein homology-based host prediction approaches enhanced the existing knowledge on the diversity of Bacteroidetes viruses and their potential interaction with their hosts in marine environments.

Diel cycling of the cosmopolitan abundant Pelagibacter virus 37-F6: one of the most abundant viruses on earth

The spatiotemporal dynamics for marine viral populations has only recently been explored. However, nothing is known about temporal activities of the uncultured Pelagibacter virus vSAG 37-F6, which was discovered by single-virus genomics as potentially the most abundant marine virus. Here, we investigate the diel cycling of 37-F6 virus and the putative SAR11 host using coastal and oceanic transcriptomic and viromic time-series data from Osaka Bay and North Pacific Subtropical Gyre. Virus 37-F6 and relatives displayed diel cycling of transcriptional activities synchronized with its putative host. In both virus and host, the lowest transcription rates were observed at 14:00-15:00, coinciding roughly with maximum solar irradiance, while higher transcriptional rates were detected during the night/early morning and afternoon. Diel abundance of free viruses of 37-F6 in seawater roughly mirrored the transcriptional activities of both virus and host. In Osaka Bay, among viral relatives (genus level), virus 37-F6 specifically showed the highest ratio of transcriptional activity to virome abundance, a proxy for viral transcriptional activity relative to free viral particle abundance. This high ratio suggests high infection rate efficiencies in vSAG 37-F6 virus compared to viral relatives. Thus, time-series data revealed temporal transcript activities in one of the most abundant viruses in Earth.

Absolute quantification of infecting viral particles by chip-based digital polymerase chain reaction

In silico and empirical quantification of viruses is paramount for obtaining information on viral populations that have a major impact on biogeochemical cycles. The uncultured Pelagibacter virus vSAG 37-F6 discovered via single-virus genomics is one of the most abundant and cosmopolitan marine viruses; however, little is understood about its temporal variation. Here, we estimated the absolute number of infecting 37-F6 viruses in coastal bacterioplankton from the Mediterranean Sea by using a novel, feasible SYBR Green I chip-based digital PCR (SYBR dPCR) technique, not implemented before for enumerating (uncultured) microbes. Quantitative SYBR dPCR estimated 450-3480 genome copies of virus 37-F6 in cells/mL (i.e. infecting viruses) and a total of ≈10-400 putative infected cells/mL with a potential C release of 0.12-4.9 pg/ml in the analysed samples. Considering that virus 37-F6 is ubiquitous and abundant in all Tara samples, an enormous amount of C could be transformed by this virus through the 'viral shunt'. Thus, this SYBR dPCR technique has enabled the absolute quantification of an ecologically relevant uncultured virus in nature and the estimation of its potential contribution on biogeochemical cycles. Overall, our study also shows that this approach has a broad applicability for quantifying any other target loci in Microbiology and Virology.

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22-5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.

Metabolic and biogeochemical consequences of viral infection in aquatic ecosystems

Ecosystems are controlled by 'bottom-up' (resources) and 'top-down' (predation) forces. Viral infection is now recognized as a ubiquitous top-down control of microbial growth across ecosystems but, at the same time, cell death by viral predation influences, and is influenced by, resource availability. In this Review, we discuss recent advances in understanding the biogeochemical impact of viruses, focusing on how metabolic reprogramming of host cells during lytic viral infection alters the flow of energy and nutrients in aquatic ecosystems. Our synthesis revealed several emerging themes. First, viral infection transforms host metabolism, in part through virus-encoded metabolic genes; the functions performed by these genes appear to alleviate energetic and biosynthetic bottlenecks to viral production. Second, viral infection depends on the physiological state of the host cell and on environmental conditions, which are challenging to replicate in the laboratory. Last, metabolic reprogramming of infected cells and viral lysis alter nutrient cycling and carbon export in the oceans, although the net impacts remain uncertain. This Review highlights the need for understanding viral infection dynamics in realistic physiological and environmental contexts to better predict their biogeochemical consequences.

Cooccurrence of Broad- and Narrow-Host-Range Viruses Infecting the Bloom-Forming Toxic Cyanobacterium Microcystis aeruginosa

Viruses play important roles in regulating the abundance and composition of bacterial populations in aquatic ecosystems. The bloom-forming toxic cyanobacterium Microcystis aeruginosa is predicted to interact with diverse cyanoviruses, resulting in Microcystis population diversification. However, current knowledge of the genomes from these viruses and their infection programs is limited to those of Microcystis virus Ma-LMM01. Here, we performed a time series sampling at a small pond in Japan during a Microcystis bloom and then investigated the genomic information and transcriptional dynamics of Microcystis-interacting viruses using metagenomic and metatranscriptomic approaches. We identified 15 viral genomic fragments classified into three groups, groups I (including Ma-LMM01), II (high abundance and transcriptional activity), and III (new lineages). According to the phylogenetic distribution of Microcystis strains possessing spacers against each viral group, the group II-original viruses interacted with all three phylogenetically distinct Microcystis population types (phylotypes), whereas the groups I and III-original viruses interacted with only one or two phylotypes, indicating the cooccurrence of broad- (group II) and narrow (groups I and III)-host-range viruses in the bloom. These viral fragments showed the highest transcriptional levels during daytime regardless of their genomic differences. Interestingly, M. aeruginosa expressed antiviral defense genes in the environment, unlike what was seen with an Ma-LMM01 infection in a previous culture experiment. Given that broad-host-range viruses often induce antiviral responses within alternative hosts, our findings suggest that such antiviral responses might inhibit viral multiplication, mainly that of broad-host-range viruses like those in group II.IMPORTANCE The bloom-forming toxic cyanobacterium Microcystis aeruginosa is thought to have diversified its population through the interactions between host and viruses in antiviral defense systems. However, current knowledge of viral genomes and infection programs is limited to those of Microcystis virus Ma-LMM01, which was a narrow host range in which it can escape from the highly abundant host defense systems. Our metagenomic approaches unveiled the cooccurrence of narrow- and broad-host-range Microcystis viruses, which included fifteen viral genomic fragments from Microcystis blooms that were classified into three groups. Interestingly, Microcystis antiviral defense genes were expressed against viral infection in the environment, unlike what was seen in a culture experiment with Ma-LMM01. Given that viruses with a broad host range often induce antiviral responses within alternative hosts, our findings suggest that antiviral responses inhibit viral reproduction, especially that of broad-range viruses like those in group II. This paper augments our understanding of the interactions between M. aeruginosa and its viruses and fills an important knowledge gap.

Tide driven microbial dynamics through virus-host interactions in the estuarine ecosystem

Microbes drive ecosystems and their viruses manipulate these processes, yet the importance of tidal functioning on the estuarine viruses and microbes remains poorly elucidated. Here, an integrative investigation on tidal patterns in viral and microbial communities and their inherent interactions over an entire spring-neap tidal cycle was conducted along a macrotidal subtropical estuary. The viral and microbial abundances oscillated significantly over the tidal cycle with relatively higher abundances observed at spring tide compared to neap tide. The distinct tidal dynamic patterns in bacterial production and community composition were tightly associated with the variations in viral infection, production and decay, revealing the tide-driven interactions between viruses and microbes. Concurrent with the higher viral decay but lower bacterial abundance and inhibited bacterial metabolism during the neap tide, lower gross viral production was coupled with a synchronous switching from viral lytic to lysogenic infection induced by the loss of viral infection efficiency and the transition from marine to freshwater bacterial populations triggered by tidal mixing. Our results highlighted the major tidal impact on the microbial dynamics through virus-host interactions, with cascading effects, neglected so far, on estuarine biogeochemical cycles.

Cyanobacterial viruses exhibit diurnal rhythms during infection

To adapt to the daily light–dark cycle, diurnal rhythms are used by the photosynthetic cyanobacteria Prochlorococcus and Synechococcus, which are the most abundant photosynthetic organisms on earth. Field studies revealed that cyanobacterial virus (cyanophage) populations in the oceans showed transcriptional rhythms. To explore the underlying mechanism, we used cyanophage laboratory cultures to find that some showed adsorption rhythms and all showed transcriptional rhythms. We discovered that the cyanophage transcriptional rhythm is partially caused by the photosynthetic activity of host cells, explaining transcriptional rhythms of field cyanophage populations. Our study shows that cultured viruses have diurnal infection rhythms which are critical for understanding how light–dark cycles shape the interaction of cyanophages and their hosts in the oceans.

As an adaptation to the daily light–dark (diel) cycle, cyanobacteria exhibit diurnal rhythms of gene expression and cell cycle. The light–dark cycle also affects the life cycle of viruses (cyanophages) that infect the unicellular picocyanobacteria Prochlorococcus and Synechococcus, which are the major primary producers in the oceans. For example, the adsorption of some cyanophages to the host cells depends on light, and the burst sizes of cyanophages are positively correlated to the length of light exposure during infection. Recent metatranscriptomic studies revealed transcriptional rhythms of field cyanophage populations. However, the underlying mechanism remains to be determined, as cyanophage laboratory cultures have not been shown to exhibit diurnal transcriptional rhythms. Here, we studied variation in infection patterns and gene expression of Prochlorococcus phages in laboratory culture conditions as a function of light. We found three distinct diel-dependent life history traits in dark conditions (diel traits): no adsorption (cyanophage P-HM2), adsorption but no replication (cyanophage P-SSM2), and replication (cyanophage P-SSP7). Under light–dark cycles, each cyanophage exhibited rhythmic transcript abundance, and cyanophages P-HM2 and P-SSM2 also exhibited rhythmic adsorption patterns. Finally, we show evidence to link the diurnal transcriptional rhythm of cyanophages to the photosynthetic activity of the host, thus providing a mechanistic explanation for the field observations of cyanophage transcriptional rhythms. Our study identifies that cultured viruses can exhibit diurnal rhythms during infection, which might impact cyanophage population-level dynamics in the oceans.

Unraveling the Contribution of High Temperature Stage to Jiang-Flavor Daqu, a Liquor Starter for Production of Chinese Jiang-Flavor Baijiu, With Special Reference to Metatranscriptomics

Jiang-flavor (JF) daqu is a liquor starter used for production of JF baijiu, a well-known distilled liquor in China. Although a high temperature stage (70°C) is necessary for qualifying JF daqu, little is known regarding its active microbial community and functional enzymes, along with its role in generating flavor precursors for JF baijiu aroma. In this investigation, based on metatranscriptomics, fungi, such as Aspergillus and Penicillium, were identified as the most active microbial members and 230 carbohydrate-active enzymes were identified as potential saccharifying enzymes at 70°C of JF daqu. Notably, most of enzymes in identified carbohydrate and energy pathways showed lower expression levels at 70°C of JF daqu than those at the high temperature stage (62°C) of Nong-flavor (NF) daqu, indicating lowering capacities of saccharification and fermentation by high temperature stage. Moreover, many enzymes, especially those related to the degradation of aromatic compounds, were only detected with low expression levels at 70°C of JF daqu albeit not at 62°C of NF daqu, indicating enhancing capacities of generating special trace aroma compounds in JF daqu by high temperature stage. Additionally, most of enzymes related to those capacities were highly expressed at 70°C by fungal genus of Aspergillus, Coccidioides, Paracoccidioides, Penicillium, and Rasamsonia. Therefore, this study not only sheds light on the crucial functions of high temperature stage but also paves the way to improve the quality of JF baijiu and provide active community and functional enzymes for other fermentation industries.

Adaptation to sub-optimal hosts is a driver of viral diversification in the ocean

Cyanophages of the Myoviridae family include generalist viruses capable of infecting a wide range of hosts including those from different cyanobacterial genera. While the influence of phages on host evolution has been studied previously, it is not known how the infection of distinct hosts influences the evolution of cyanophage populations. Here, using an experimental evolution approach, we investigated the adaptation of multiple cyanophage populations to distinct cyanobacterial hosts. We show that when infecting an "optimal" host, whose infection is the most efficient, phage populations accumulated only a few mutations. However, when infecting "sub-optimal" hosts, different mutations spread in the phage populations, leading to rapid diversification into distinct subpopulations. Based on our results, we propose a model demonstrating how shifts in microbial abundance, which lead to infection of "sub-optimal" hosts, act as a driver for rapid diversification of viral populations.

Degenerate PCR Primers to Reveal the Diversity of Giant Viruses in Coastal Waters

“Megaviridae” is a proposed family of giant viruses infecting unicellular eukaryotes. These viruses are ubiquitous in the sea and have impact on marine microbial community structure and dynamics through their lytic infection cycle. However, their diversity and biogeography have been poorly characterized due to the scarce detection of Megaviridae sequences in metagenomes, as well as the limitation of reference sequences used to design specific primers for this viral group. Here, we propose a set of 82 degenerated primers (referred to as MEGAPRIMER), targeting DNA polymerase genes (polBs) of Megaviridae. MEGAPRIMER was designed based on 921 Megaviridae polBs from sequenced genomes and metagenomes. By applying this primer set to environmental DNA meta-barcoding of a coastal seawater sample, we report 5595 non-singleton operational taxonomic units (OTUs) of Megaviridae at 97% nucleotide sequence identity. The majority of the OTUs were found to form diverse clades, which were phylogenetically distantly related to known viruses such as Mimivirus. The Megaviridae OTUs detected in this study outnumber the giant virus OTUs identified in previous individual studies by more than an order of magnitude. Hence, MEGAPRIMER represents a useful tool to study the diversity of Megaviridae at the population level in natural environments.

Contrasting biogeography and diversity patterns between diatoms and haptophytes in the central Pacific Ocean

Diatoms and haptophytes are two major phytoplankton groups, playing pivotal roles in global biogeochemical cycles and marine ecosystems. In general, diatoms have higher growth rates than haptophytes, whereas haptophytes tend to have higher nutrient uptake affinity. However, precise linkages between their ecological traits and geographical distributions remain poorly understood. Herein, we examined the basin-scale variability of the abundance and taxonomic composition of these two phytoplankton groups across 35 sites in the Pacific Ocean using DNA metabarcoding. The diatom community was generally dominated by a few genera at each sample site, whereas the haptophyte community consisted of a large number of genera in most of the sites. The coexistence of various haptophyte genera might be achieved by diversification of their ecophysiological traits such as mixotrophy. On the other hand, the diatom community might experience greater inter-genus competition due to the rapid uptake of nutrients. Our data further supports the notion that their distinct ecological strategies underlie the emergence of contrasting diversity patterns of these phytoplankton groups in the central Pacific at a basin scale.