Light regulation of coccolithophore host-virus interactions

Expansion of the genomic and functional diversity of global ocean giant viruses

Giant viruses (GVs) play crucial roles in the global ocean microbial food web and biogeochemistry. Recent metagenomic advances have uncovered >1800 new GV genomes from the world's oceans. While this rapid increase in genomic information is impressive, it is nowhere close to the extensive genomic information available for other marine entities-e.g., prokaryotes and their "virome". We present 230 new high-quality GV genomes (genomes with 4 or more marker genes) and 398 partial genomes from nine global ocean datasets. Notably, we identified numerous GV genomes from the Baltic Sea, offering insights into their phylogenomics, metabolic potential, and environmental drivers in one of the largest brackish water ecosystems. We discovered new GV functions and identified a significant functional divide between the Imitervirales and Algavirales orders. Additionally, we evaluated factors affecting GV abundance through a case study on the Baltic Sea dataset. Our study significantly expands the marine GV genomic and functional diversity, broadening our understanding of their roles in the food web and biogeochemistry.

Phylogenetic proximity drives temporal succession of marine giant viruses in a five-year metagenomic time-series

Nucleocytoplasmic Large DNA Viruses (NCLDVs, also called giant viruses) are widespread in marine systems and infect a broad range of microbial eukaryotes (protists). Recent biogeographic work has provided global snapshots of NCLDV diversity and community composition across the world's oceans, yet little information exists about the guiding 'rules' underpinning their community dynamics over time. We leveraged a five-year monthly metagenomic time-series to quantify the community composition of NCLDVs off the coast of Southern California and characterize these populations' temporal dynamics. NCLDVs were dominated by Algavirales (Phycodnaviruses, 59%) and Imitervirales (Mimiviruses, 36%). We identified clusters of NCLDVs with distinct classes of seasonal and non-seasonal temporal dynamics. Overall, NCLDV population abundances were often highly dynamic with a strong seasonal signal. The Imitervirales group had highest relative abundance in the more oligotrophic late summer and fall, while Algavirales did so in winter. Generally, closely related strains had similar temporal dynamics, suggesting that evolutionary history is a key driver of the temporal niche of marine NCLDVs. However, a few closely-related strains had drastically different seasonal dynamics, suggesting that while phylogenetic proximity often indicates ecological similarity, occasionally phenology can shift rapidly, possibly due to host-switching. Finally, we identified distinct functional content and possible host interactions of two major NCLDV orders-including connections of Imitervirales with primary producers like the diatom Chaetoceros and widespread marine grazers like Paraphysomonas and Spirotrichea ciliates. Together, our results reveal key insights on season-specific effect of phylogenetically distinct giant virus communities on marine protist metabolism, biogeochemical fluxes and carbon cycling.

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. We find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

Cell-to-cell heterogeneity drives host-virus coexistence in a bloom-forming alga

Algal blooms drive global biogeochemical cycles of key nutrients and serve as hotspots for biological interactions in the ocean. The massive blooms of the cosmopolitan coccolithophore Emiliania huxleyi are often infected by the lytic E. huxleyi virus, which is a major mortality agent triggering bloom demise. This multi-annual "boom and bust" pattern of E. huxleyi blooms suggests that coexistence is essential for these host-virus dynamics. To investigate host-virus coexistence, we developed a new model system from an E. huxleyi culture that recovered from viral infection. The recovered population coexists with the virus, as host cells continue to divide in parallel to viral production. By applying single-molecule fluorescence in situ hybridization (smFISH) to quantify the fraction of infected cells, and assessing infection-specific lipid biomarkers, we identified a small subpopulation of cells that were infected and produced new virions, whereas most of the host population could resist infection. To further assess population heterogeneity, we generated clonal strain collections using single-cell sorting and subsequently phenotyped their susceptibility to E. huxleyi virus infection. This unraveled substantial cell-to-cell heterogeneity across a continuum of susceptibility to resistance, highlighting that infection outcome may vary depending on the individual cell. These results add a new dimension to our understanding of the complexity of host-virus interactions that are commonly assessed in bulk and described by binary definitions of resistance or susceptibility. We propose that phenotypic heterogeneity drives the host-virus coexistence and demonstrate how the coexistence with a lytic virus provides an ecological advantage for the host by killing competing strains.

Daily turnover of active giant virus infection during algal blooms revealed by single-cell transcriptomics

Giant viruses infect many unicellular eukaryotes, including algae that form massive oceanic blooms. Despite the major impact of viruses on the marine ecosystem, the ability to quantify and assess active viral infection in nature remains a major challenge. We applied single-cell RNA sequencing, to profile virus and host transcriptomes of 12,000 single algal cells from a coccolithophore bloom. Viral infection was detected already at early exponential bloom phase, negatively correlating with the bloom intensity. A consistent percent of infected coccolithophores displayed the early phase of viral replication for several consecutive days, indicating a daily turnover and continuous virocell-associated metabolite production, potentially affecting the surrounding microbiome. Linking single-cell infection state to host physiology revealed that infected cells remained calcified even in the late infection stage. These findings stress the importance of studying host-virus dynamics in natural populations, at single-cell resolution, to better understand virus life cycle and its impact on microbial food webs.

Marine phytoplankton downregulate core photosynthesis and carbon storage genes upon rapid mixed layer shallowing

Marine phytoplankton are a diverse group of photoautotrophic organisms and key mediators in the global carbon cycle. Phytoplankton physiology and biomass accumulation are closely tied to mixed layer depth, but the intracellular metabolic pathways activated in response to changes in mixed layer depth remain less explored. Here, metatranscriptomics was used to characterize the phytoplankton community response to a mixed layer shallowing (from 233 to 5 m) over the course of two days during the late spring in the Northwest Atlantic. Most phytoplankton genera downregulated core photosynthesis, carbon storage, and carbon fixation genes as the system transitioned from a deep to a shallow mixed layer and shifted towards catabolism of stored carbon supportive of rapid cell growth. In contrast, phytoplankton genera exhibited divergent transcriptional patterns for photosystem light harvesting complex genes during this transition. Active virus infection, taken as the ratio of virus to host transcripts, increased in the Bacillariophyta (diatom) phylum and decreased in the Chlorophyta (green algae) phylum upon mixed layer shallowing. A conceptual model is proposed to provide ecophysiological context for our findings, in which integrated light limitation and lower division rates during transient deep mixing are hypothesized to disrupt resource-driven, oscillating transcript levels related to photosynthesis, carbon fixation, and carbon storage. Our findings highlight shared and unique transcriptional response strategies within phytoplankton communities acclimating to the dynamic light environment associated with transient deep mixing and shallowing events during the annual North Atlantic bloom.

Marine viruses and climate change: Virioplankton, the carbon cycle, and our future ocean

Interactions between marine viruses and microbes are a critical part of the oceanic carbon cycle. The impacts of virus-host interactions range from short-term disruptions in the mobility of microbial biomass carbon to higher trophic levels through cell lysis (i.e., the viral shunt) to long-term reallocation of microbial biomass carbon to the deep sea through accelerating the biological pump (i.e., the viral shuttle). The biogeochemical backdrop of the ocean-the physical, chemical, and biological landscape-influences the likelihood of both virus-host interactions and particle formation, and the fate and flow of carbon. As climate change reshapes the oceanic landscape through large-scale shifts in temperature, circulation, stratification, and acidification, virus-mediated carbon flux is likely to shift in response. Dynamics in the directionality and magnitude of changes in how, where, and when viruses mediate the recycling or storage of microbial biomass carbon is largely unknown. Integrating viral infection dynamics data obtained from experimental models and field systems, with particle motion microphysics and global observations of oceanic biogeochemistry, into improved ecosystem models will enable viral oceanographers to better predict the role of viruses in marine carbon cycling in the future ocean.

Proton-transporting heliorhodopsins from marine giant viruses

Rhodopsins convert light into signals and energy in animals and microbes. Heliorhodopsins (HeRs), a recently discovered new rhodopsin family, are widely present in archaea, bacteria, unicellular eukaryotes, and giant viruses, but their function remains unknown. Here, we report that a viral HeR from Emiliania huxleyi virus 202 (V2HeR3) is a light-activated proton transporter. V2HeR3 absorbs blue-green light, and the active intermediate contains the deprotonated retinal Schiff base. Site-directed mutagenesis study revealed that E191 in TM6 constitutes the gate together with the retinal Schiff base. E205 and E215 form a PAG of the Schiff base, and mutations at these positions converted the protein into an outward proton pump. Three environmental viral HeRs from the same group as well as a more distantly related HeR exhibited similar proton-transport activity, indicating that HeR functions might be diverse similarly to type-1 microbial rhodopsins. Some strains of E. huxleyi contain one HeR that is related to the viral HeRs, while its viruses EhV-201 and EhV-202 contain two and three HeRs, respectively. Except for V2HeR3 from EhV-202, none of these proteins exhibit ion transport activity. Thus, when expressed in the E. huxleyi cell membranes, only V2HeR3 has the potential to depolarize the host cells by light, possibly to overcome the host defense mechanisms or to prevent superinfection. The neuronal activity generated by V2HeR3 suggests that it can potentially be used as an optogenetic tool, similarly to type-1 microbial rhodopsins.

Visualizing active viral infection reveals diverse cell fates in synchronized algal bloom demise

Marine viruses are the most abundant biological entity in the ocean and are considered as major evolutionary drivers of microbial life [C. A. Suttle, Nat. Rev. Microbiol. 5, 801-812 (2007)]. Yet, we lack quantitative approaches to assess their impact on the marine ecosystem. Here, we provide quantification of active viral infection in the bloom forming single-celled phytoplankton Emiliania huxleyi infected by the large virus EhV, using high-throughput single-molecule messenger RNA in situ hybridization (smFISH) of both virus and host transcripts. In natural samples, viral infection reached only 25% of the population despite synchronized bloom demise exposing the coexistence of infected and noninfected subpopulations. We prove that photosynthetically active cells chronically release viral particles through nonlytic infection and that viral-induced cell lysis can occur without viral release, thus challenging major assumptions regarding the life cycle of giant viruses. We could also assess active infection in cell aggregates linking viral infection and carbon export to the deep ocean [C. P. Laber et al., Nat. Microbiol. 3, 537-547 (2018)] and suggest a potential host defense strategy by enrichment of infected cells in sinking aggregates. Our approach can be applied to diverse marine microbial systems, opening a mechanistic dimension to the study of biotic interactions in the ocean.

Spring Accumulation Rates in North Atlantic Phytoplankton Communities Linked to Alterations in the Balance Between Division and Loss

For nearly a century, phytoplankton spring blooms have largely been explained in the context of abiotic factors regulating cellular division rates (e.g., mixed-layer light levels). However, the accumulation of new phytoplankton biomass represents a mismatch between phytoplankton division and mortality rates. The balance between division and loss, therefore, has important implications for marine food webs and biogeochemical cycles. A large fraction of phytoplankton mortality is due to the combination of microzooplankton grazing and viral lysis, however, broad scale simultaneous measurements of these mortality processes are scarce. We applied the modified dilution assay along a West-to-East diagonal transect in the North Atlantic during spring. Our results demonstrate positive accumulation rates with losses dominated by microzooplankton grazing. Considering the dynamic light environment phytoplankton experience in the mixed surface layer, particularly in the spring, we tested the potential for incubation light conditions to affect observed rates. Incubations acted as short-term 'light' perturbations experiments, in which deeply mixed communities are exposed to elevated light levels. These "light perturbations" increased phytoplankton division rates and resulted in proportional changes in phytoplankton biomass while having no significant effect on mortality rates. These results provide experimental evidence for the Disturbance-Recovery Hypothesis, supporting the tenet that biomass accumulation rates co-vary with the specific rate of change in division.

Nutrient Loading and Viral Memory Drive Accumulation of Restriction Modification Systems in Bloom-Forming Cyanobacteria

The mechanisms driving cyanobacterial harmful algal blooms (HABs) like those caused by Microcystis aeruginosa remain elusive, but improved defense against viral predation has been implicated for success in eutrophic environments. Our genus-level analyses of 139,023 genomes revealed that HAB-forming cyanobacteria carry vastly more restriction modification systems per genome (RMPG) than nearly all other prokaryotic genera, suggesting that viral defense is a cornerstone of their ecological success. In contrast, picocyanobacteria that numerically dominate nutrient-poor systems have the fewest RMPG within the phylum Cyanobacteria. We used classic resource competition models to explore the hypothesis that nutrient enrichments drive ecological selection for high RMPG due to increased host-phage contact rate. These classic models, agnostic to the mechanism of defense, explain how nutrient loading can select for increased RMPG but, importantly, fail to explain the extreme accumulation of these defense systems. However, extreme accumulation of RMPG can be achieved in a novel “memory” model that accounts for a unique activity of restriction modification systems: the accidental methylation of viral DNA by the methyltransferase. The methylated virus “remembers” the RM defenses of its former host and can evade these defenses if they are present in the next host. This viral memory leads to continual RM system devaluation; RMs accumulate extensively because the benefit of each addition is diminished. Our modeling leads to the hypothesis that nutrient loading and virion methylation drive the extreme accumulation of RMPG in HAB-forming cyanobacteria. Finally, our models suggest that hosts with different RMPG values can coexist when hosts have unique sets of RM systems.

Host Range and Coding Potential of Eukaryotic Giant Viruses

Giant viruses are a group of eukaryotic double-stranded DNA viruses with large virion and genome size that challenged the traditional view of virus. Newly isolated strains and sequenced genomes in the last two decades have substantially advanced our knowledge of their host diversity, gene functions, and evolutionary history. Giant viruses are now known to infect hosts from all major supergroups in the eukaryotic tree of life, which predominantly comprises microbial organisms. The seven well-recognized viral clades (taxonomic families) have drastically different host range. Mimiviridae and Phycodnaviridae, both with notable intrafamilial genome variation and high abundance in environmental samples, have members that infect the most diverse eukaryotic lineages. Laboratory experiments and comparative genomics have shed light on the unprecedented functional potential of giant viruses, encoding proteins for genetic information flow, energy metabolism, synthesis of biomolecules, membrane transport, and sensing that allow for sophisticated control of intracellular conditions and cell-environment interactions. Evolutionary genomics can illuminate how current and past hosts shape viral gene repertoires, although it becomes more obscure with divergent sequences and deep phylogenies. Continued works to characterize giant viruses from marine and other environments will further contribute to our understanding of their host range, coding potential, and virus-host coevolution.

Making sense of virus size and the tradeoffs shaping viral fitness

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

Temperate infection in a virus-host system previously known for virulent dynamics

The blooming cosmopolitan coccolithophore Emiliania huxleyi and its viruses (EhVs) are a model for density-dependent virulent dynamics. EhVs commonly exhibit rapid viral reproduction and drive host death in high-density laboratory cultures and mesocosms that simulate blooms. Here we show that this system exhibits physiology-dependent temperate dynamics at environmentally relevant E. huxleyi host densities rather than virulent dynamics, with viruses switching from a long-term non-lethal temperate phase in healthy hosts to a lethal lytic stage as host cells become physiologically stressed. Using this system as a model for temperate infection dynamics, we present a template to diagnose temperate infection in other virus-host systems by integrating experimental, theoretical, and environmental approaches. Finding temperate dynamics in such an established virulent host-virus model system indicates that temperateness may be more pervasive than previously considered, and that the role of viruses in bloom formation and decline may be governed by host physiology rather than by host-virus densities.

A single-cell view on alga-virus interactions reveals sequential transcriptional programs and infection states

The discovery of giant viruses infecting eukaryotes from diverse ecosystems has revolutionized our understanding of the evolution of viruses and their impact on protist biology, yet knowledge on their replication strategies and transcriptome regulation remains limited. Here, we profile single-cell transcriptomes of the globally distributed microalga Emiliania huxleyi and its specific giant virus during infection. We detected profound heterogeneity in viral transcript levels among individual cells. Clustering single cells based on viral expression profiles enabled reconstruction of the viral transcriptional trajectory. Reordering cells along this path unfolded highly resolved viral genetic programs composed of genes with distinct promoter elements that orchestrate sequential expression. Exploring host transcriptome dynamics across the viral infection states revealed rapid and selective shutdown of protein-encoding nuclear transcripts, while the plastid and mitochondrial transcriptomes persisted into later stages. Single-cell RNA-seq opens a new avenue to unravel the life cycle of giant viruses and their unique hijacking strategies.

A Primer for Microbiome Time-Series Analysis

Time-series can provide critical insights into the structure and function of microbial communities. The analysis of temporal data warrants statistical considerations, distinct from comparative microbiome studies, to address ecological questions. This primer identifies unique challenges and approaches for analyzing microbiome time-series. In doing so, we focus on (1) identifying compositionally similar samples, (2) inferring putative interactions among populations, and (3) detecting periodic signals. We connect theory, code and data via a series of hands-on modules with a motivating biological question centered on marine microbial ecology. The topics of the modules include characterizing shifts in community structure and activity, identifying expression levels with a diel periodic signal, and identifying putative interactions within a complex community. Modules are presented as self-contained, open-access, interactive tutorials in R and Matlab. Throughout, we highlight statistical considerations for dealing with autocorrelated and compositional data, with an eye to improving the robustness of inferences from microbiome time-series. In doing so, we hope that this primer helps to broaden the use of time-series analytic methods within the microbial ecology research community.

Linking Light-Dependent Life History Traits with Population Dynamics for Prochlorococcus and Cyanophage

Prochlorococcus cyanobacteria grow in diurnal rhythms driven by diel cycles. Their ecology depends on light, nutrients, and top-down mortality processes, including lysis by viruses. Cyanophage, viruses that infect cyanobacteria, are also impacted by light. For example, the extracellular viability and intracellular infection kinetics of some cyanophage vary between light and dark conditions. Nonetheless, it remains unclear whether light-dependent viral life history traits scale up to influence population-level dynamics. Here, we examined the impact of diel forcing on both cellular- and population-scale dynamics in multiple Prochlorococcus-phage systems. To do so, we developed a light-driven population model, including both cellular growth and viral infection dynamics. We then tested the model against measurements of experimental infection dynamics with diel forcing to examine the extent to which population level changes in both viral and host abundances could be explained by light-dependent life history traits. Model-data integration reveals that light-dependent adsorption can improve fits to population dynamics for some virus-host pairs. However, light-dependent variation alone does not fully explain realized host and virus population dynamics. Instead, we show evidence consistent with lysis saturation at relatively high virus-to-cell ratios. Altogether, our study represents a quantitative approach to integrate mechanistic models to reconcile Prochlorococcus-virus dynamics spanning cellular-to-population scales.IMPORTANCE The cyanobacterium Prochlorococcus is an essential member of global ocean ecosystems. Light rhythms drive Prochlorococcus photosynthesis, ecology, and interactions with potentially lethal viruses. At present, the impact of light on Prochlorococcus-virus interactions is not well understood. Here, we analyzed Prochlorococcus and virus population dynamics with a light-driven population model and compared our results with experimental data. Our approach revealed that light profoundly drives both cellular- and population-level dynamics for some host-virus systems. However, we also found that additional mechanisms, including lysis saturation, are required to explain observed host-virus dynamics at the population scale. This study provides the basis for future work to understand the intertwined fates of Prochlorococcus and associated viruses in the surface ocean.

Influence of light on the infection of Aureococcus anophagefferens CCMP 1984 by a "giant virus"

The pelagophyte Aureococcus anophagefferens has caused recurrent brown tide blooms along the northeast coast of the United States since the mid-1980's, and more recently spread to other regions of the globe. These blooms, due to the high cell densities, are associated with severe light attenuation that destroys the sea grass beds which provide the basis for many fisheries. Data collected by transmission electron microscopy, PCR, and metatranscriptomic studies of the blooms, support the hypothesis that large dsDNA viruses play a role in bloom dynamics. While a large (~140 nm) icosahedral virus, with a 371 kbp genome, was first isolated more than a decade ago, the constraints imposed by environmental parameters on bloom infection dynamics by Aureococcus anophagefferens Virus, (AaV) remain unknown. To investigate the role light plays in infection by this virus, we acclimated A. anophagefferens to light intensities of 30 (low), 60 (medium) or 90 μmol photons m-2 s-1 (high) and infected cultures at these irradiance levels. Moreover, we completed light shift experiments where acclimated cultures were exposed to even lower light intensities (0, 5, and 15 μmol photons m-2 s-1) consistent with irradiance found during the peak of the bloom when cell concentrations are highest. The abundance of viruses produced per lytic event (burst size) was lower in the low irradiance acclimated cultures compared to the medium and high acclimated cultures. Transferring infected cultures to more-limiting light availabilities further decreased burst size and increased the length of time it took for cultures to lyse, regardless of acclimation irradiance level. A hypothetical mechanism for the reduced efficiency of the infection cycle in low light due to ribosome biogenesis was predicted from pre-existing transcriptomes. Overall, these studies provide a framework for understanding light effects on infection dynamics over the course of the summer months when A. anophagefferens blooms occur.

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.

Unmasking cellular response of a bloom-forming alga to viral infection by resolving expression profiles at a single-cell level

Infection by large dsDNA viruses can lead to a profound alteration of host transcriptome and metabolome in order to provide essential building blocks to support the high metabolic demand for viral assembly and egress. Host response to viral infection can typically lead to diverse phenotypic outcome that include shift in host life cycle and activation of anti-viral defense response. Nevertheless, there is a major bottleneck to discern between viral hijacking strategies and host defense responses when averaging bulk population response. Here we study the interaction between Emiliania huxleyi, a bloom-forming alga, and its specific virus (EhV), an ecologically important host-virus model system in the ocean. We quantified host and virus gene expression on a single-cell resolution during the course of infection, using automatic microfluidic setup that captures individual algal cells and multiplex quantitate PCR. We revealed high heterogeneity in viral gene expression among individual cells. Simultaneous measurements of expression profiles of host and virus genes at a single-cell level allowed mapping of infected cells into newly defined infection states and allowed detection specific host response in a subpopulation of infected cell which otherwise masked by the majority of the infected population. Intriguingly, resistant cells emerged during viral infection, showed unique expression profiles of metabolic genes which can provide the basis for discerning between viral resistant and susceptible cells within heterogeneous populations in the marine environment. We propose that resolving host-virus arms race at a single-cell level will provide important mechanistic insights into viral life cycles and will uncover host defense strategies.

Almost all of our current understanding of the molecular mechanisms that govern host-pathogen interactions in the ocean is derived from experiments carried out at the population level, neglecting any heterogeneity. Here we used a single cell approach to unmask the phenotypic heterogeneity produced within infected populations of the cosmopolitan bloom-forming alga Emiliania huxleyi by its specific lytic virus. We found high variability in expression of viral genes among individual cells. This heterogeneity was used to map cells into their infection state and allowed to uncover a yet unrecognized host response. We also provide evidence that variability in host metabolic states provided a sensitive tool to decipher between susceptible and resistant cells.