Dynamics and genotypic composition of Emiliania huxleyi and their co-occurring viruses during a coccolithophore bloom in the North Sea

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.

Macroalgal virosphere assists with host-microbiome equilibrium regulation and affects prokaryotes in surrounding marine environments

The microbiomes in macroalgal holobionts play vital roles in regulating macroalgal growth and ocean carbon cycling. However, the virospheres in macroalgal holobionts remain largely underexplored, representing a critical knowledge gap. Here we unveil that the holobiont of kelp (Saccharina japonica) harbors highly specific and unique epiphytic/endophytic viral species, with novelty (99.7% unknown) surpassing even extreme marine habitats (e.g. deep-sea and hadal zones), indicating that macroalgal virospheres, despite being closest to us, are among the least understood. These viruses potentially maintain microbiome equilibrium critical for kelp health via lytic-lysogenic infections and the expression of folate biosynthesis genes. In-situ kelp mesocosm cultivation and metagenomic mining revealed that kelp holobiont profoundly reshaped surrounding seawater and sediment virus-prokaryote pairings through changing surrounding environmental conditions and virus-host migrations. Some kelp epiphytic viruses could even infect sediment autochthonous bacteria after deposition. Moreover, the presence of ample viral auxiliary metabolic genes for kelp polysaccharide (e.g. laminarin) degradation underscores the underappreciated viral metabolic influence on macroalgal carbon cycling. This study provides key insights into understanding the previously overlooked ecological significance of viruses within macroalgal holobionts and the macroalgae-prokaryotes-virus tripartite relationship.

Altered growth and death in dilution-based viral predation assays

Viral lysis of phytoplankton is one of the most common forms of death on Earth. Building on an assay used extensively to assess rates of phytoplankton loss to predation by grazers, lysis rates are increasingly quantified through dilution-based techniques. In this approach, dilution of viruses and hosts are expected to reduce infection rates and thus increase host net growth rates (i.e., accumulation rates). The difference between diluted and undiluted host growth rates is interpreted as a measurable proxy for the rate of viral lytic death. These assays are usually conducted in volumes ≥ 1 L. To increase throughput, we implemented a miniaturized, high-throughput, high-replication, flow cytometric microplate dilution assay to measure viral lysis in environmental samples sourced from a suburban pond and the North Atlantic Ocean. The most notable outcome we observed was a decline in phytoplankton densities that was exacerbated by dilution, instead of the increased growth rates expected from lowered virus-phytoplankton encounters. We sought to explain this counterintuitive outcome using theoretical, environmental, and experimental analyses. Our study shows that, while die-offs could be partly explained by a 'plate effect' due to small incubation volumes and cells adhering to walls, the declines in phytoplankton densities are not volume-dependent. Rather, they are driven by many density- and physiology-dependent effects of dilution on predation pressure, nutrient limitation, and growth, all of which violate the original assumptions of dilution assays. As these effects are volume-independent, these processes likely occur in all dilution assays that our analyses show to be remarkably sensitive to dilution-altered phytoplankton growth and insensitive to actual predation pressure. Incorporating altered growth as well as predation, we present a logical framework that categorizes locations by the relative dominance of these mechanisms, with general applicability to dilution-based assays.

Spatiotemporal Dynamics of Coastal Viral Community Structure and Potential Biogeochemical Roles Affected by an Ulva prolifera Green Tide

The world's largest macroalgal green tide, caused by Ulva prolifera, has resulted in serious consequences for coastal waters of the Yellow Sea, China. Although viruses are considered to be one of the key factors in controlling microalgal bloom demise, understanding of the relationship between viral communities and the macroalgal green tide is still poor. Here, a Qingdao coastal virome (QDCV) time-series data set was constructed based on the metagenomic analysis of 17 DNA viromes along three coastal stations of the Yellow Sea, covering different stages of the green tide from Julian days 165 to 271. A total of 40,076 viral contigs were detected and clustered into 28,058 viral operational taxonomic units (vOTUs). About 84% of the vOTUs could not be classified, and 62% separated from vOTUs in other ecosystems. Green tides significantly influenced the spatiotemporal dynamics of the viral community structure, diversity, and potential functions. For the classified vOTUs, the relative abundance of Pelagibacter phages declined with the arrival of the bloom and rebounded after the bloom, while Synechococcus and Roseobacter phages increased, although with a time lag from the peak of their hosts. More than 80% of the vOTUs reached peaks in abundance at different specific stages, and the viral peaks were correlated with specific hosts at different stages of the green tide. Most of the viral auxiliary metabolic genes (AMGs) were associated with carbon and sulfur metabolism and showed spatiotemporal dynamics relating to the degradation of the large amount of organic matter released by the green tide. IMPORTANCE To the best of our knowledge, this study is the first to investigate the responses of viruses to the world's largest macroalgal green tide. It revealed the spatiotemporal dynamics of the unique viral assemblages and auxiliary metabolic genes (AMGs) following the variation and degradation of Ulva prolifera. These findings demonstrate a tight coupling between viral assemblages, and prokaryotic and eukaryotic abundances were influenced by the green tide.

Aureococcus anophagefferens (Pelagophyceae) genomes improve evaluation of nutrient acquisition strategies involved in brown tide dynamics

The pelagophyte Aureococcus anophagefferens causes harmful brown tide blooms in marine embayments on three continents. Aureococcus anophagefferens was the first harmful algal bloom species to have its genome sequenced, an advance that evidenced genes important for adaptation to environmental conditions that prevail during brown tides. To expand the genomic tools available for this species, genomes for four strains were assembled, including three newly sequenced strains and one assembled from publicly available data. These genomes ranged from 57.11 to 73.62 Mb, encoding 13,191-17,404 potential proteins. All strains shared ~90% of their encoded proteins as determined by homology searches and shared most functional orthologs as determined by KEGG, although each strain also possessed coding sequences with unique functions. Like the original reference genome, the genomes assembled in this study possessed genes hypothesized to be important in bloom proliferation, including genes involved in organic compound metabolism and growth at low light. Cross-strain informatics and culture experiments suggest that the utilization of purines is a potentially important source of organic nitrogen for brown tides. Analyses of metatranscriptomes from a brown tide event demonstrated that use of a single genome yielded a lower read mapping percentage (~30% of library reads) as compared to a database generated from all available genomes (~43%), suggesting novel information about bloom ecology can be gained from expanding genomic space. This work demonstrates the continued need to sequence ecologically relevant algae to understand the genomic potential and their ecology in the environment.

Viruses in Extreme Environments, Current Overview, and Biotechnological Potential

Virus research has advanced significantly since the discovery of the tobacco mosaic virus (TMV), the characterization of its infection mechanisms and the factors that determine their pathogenicity. However, most viral research has focused on pathogenic viruses to humans, animals and plants, which represent only a small fraction in the virosphere. As a result, the role of most viral genes, and the mechanisms of coevolution between mutualistic viruses, their host and their environment, beyond pathogenicity, remain poorly understood. This review focuses on general aspects of viruses that interact with extremophile organisms, characteristics and examples of mechanisms of adaptation. Finally, this review provides an overview on how knowledge of extremophile viruses sheds light on the application of new tools of relevant use in modern molecular biology, discussing their value in a biotechnological context.

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.

Temporal Variability of Virioplankton during a Gymnodinium catenatum Algal Bloom

Viruses are key biogeochemical engines in the regulation of the dynamics of phytoplankton. However, there has been little research on viral communities in relation to algal blooms. Using the virMine tool, we analyzed viral information from metagenomic data of field dinoflagellate (Gymnodinium catenatum) blooms at different stages. Species identification indicated that phages were the main species. Unifrac analysis showed clear temporal patterns in virioplankton dynamics. The viral community was dominated by Siphoviridae, Podoviridae, and Myoviridae throughout the whole bloom cycle. However, some changes were observed at different phases of the bloom; the relatively abundant Siphoviridae and Myoviridae dominated at pre-bloom and peak bloom stages, while at the post-bloom stage, the members of Phycodnaviridae and Microviridae were more abundant. Temperature and nutrients were the main contributors to the dynamic structure of the viral community. Some obvious correlations were found between dominant viral species and host biomass. Functional analysis indicated some functional genes had dramatic response in algal-associated viral assemblages, especially the CAZyme encoding genes. This work expands the existing knowledge of algal-associated viruses by characterizing viral composition and function across a complete algal bloom cycle. Our data provide supporting evidence that viruses participate in dinoflagellate bloom dynamics under natural conditions.

Biochemical diversity of glycosphingolipid biosynthesis as a driver of Coccolithovirus competitive ecology

Coccolithoviruses (EhVs) are large, double-stranded DNA-containing viruses that infect the single-celled, marine coccolithophore Emiliania huxleyi. Given the cosmopolitan nature and global importance of E. huxleyi as a bloom-forming, calcifying, photoautotroph, E. huxleyi-EhV interactions play a key role in oceanic carbon biogeochemistry. Virally-encoded glycosphingolipids (vGSLs) are virulence factors that are produced by the activity of virus-encoded serine palmitoyltransferase (SPT). Here, we characterize the dynamics, diversity and catalytic production of vGSLs in an array of EhV strains in relation to their SPT sequence composition and explore the hypothesis that they are a determinant of infectivity and host demise. vGSL production and diversity was positively correlated with increased virulence, virus replication rate and lytic infection dynamics in laboratory experiments, but they do not explain the success of less-virulent EhVs in natural EhV communities. The majority of EhV-derived SPT amplicon sequences associated with infected cells in the North Atlantic derived from slower infecting, less virulent EhVs. Our lab-, field- and mathematical model-based data and simulations support ecological scenarios whereby slow-infecting, less-virulent EhVs successfully compete in North Atlantic populations of E. huxleyi, through either the preferential removal of fast-infecting, virulent EhVs during active infection or by having access to a broader host range.

Dynamics of transparent exopolymer particle production and aggregation during viral infection of the coccolithophore, Emiliania huxleyi

Emiliania huxleyi produces calcium carbonate (CaCO₃ ) coccoliths and transparent exopolymer particles (TEP), sticky, acidic carbohydrates that facilitate aggregation. E. huxleyi's extensive oceanic blooms are often terminated by coccolithoviruses (EhVs) with the transport of cellular debris and associated particulate organic carbon (POC) to depth being facilitated by TEP-bound 'marine snow' aggregates. The dynamics of TEP production and particle aggregation in response to EhV infection are poorly understood. Using flow cytometry, spectrophotometry and FlowCam visualization of alcian blue (AB)-stained aggregates, we assessed TEP production and the size spectrum of aggregates for E. huxleyi possessing different degrees of calcification and cellular CaCO₃ :POC mass ratios, when challenged with two EhVs (EhV207 and EhV99B1). FlowCam imaging also qualitatively assessed the relative amount of AB-stainable TEP (i.e., blue:red ratio of each particle). We show significant increases in TEP during early phase EhV207-infection (∼ 24 h) of calcifying strains and a shift towards large aggregates following EhV99B1-infection. We also observed the formation of large aggregates with low blue:red ratios, suggesting that other exopolymer substances contribute towards aggregation. Our findings show the potential for virus infection and the associated response of their hosts to impact carbon flux dynamics and provide incentive to explore these dynamics in natural populations.

Viruses in Marine Ecosystems: From Open Waters to Coral Reefs

Viruses infect all kingdoms of marine life from bacteria to whales. Viruses in the world's oceans play important roles in the mortality of phytoplankton, and as drivers of evolution and biogeochemical cycling. They shape host population abundance and distribution and can lead to the termination of algal blooms. As discoveries about this huge reservoir of genetic and biological diversity grow, our understanding of the major influences viruses exert in the global marine environment continues to expand. This chapter discusses the key discoveries that have been made to date about marine viruses and the current direction of this field of research.

Expression profiling of host and virus during a coccolithophore bloom provides insights into the role of viral infection in promoting carbon export

The cosmopolitan coccolithophore Emiliania huxleyi is a unicellular eukaryotic alga that forms vast blooms in the oceans impacting large biogeochemical cycles. These blooms are often terminated due to infection by the large dsDNA virus, E. huxleyi virus (EhV). It was recently established that EhV-induced modulation of E. huxleyi metabolism is a key factor for optimal viral infection cycle. Despite the huge ecological importance of this host-virus interaction, the ability to assess its spatial and temporal dynamics and its possible impact on nutrient fluxes is limited by current approaches that focus on quantification of viral abundance and biodiversity. Here, we applied a host and virus gene expression analysis as a sensitive tool to quantify the dynamics of this interaction during a natural E. huxleyi bloom in the North Atlantic. We used viral gene expression profiling as an index for the level of active infection and showed that the latter correlated with water column depth. Intriguingly, this suggests a possible sinking mechanism for removing infected cells as aggregates from the E. huxleyi population in the surface layer into deeper waters. Viral infection was also highly correlated with induction of host metabolic genes involved in host life cycle, sphingolipid, and antioxidant metabolism, providing evidence for modulation of host metabolism under natural conditions. The ability to track and quantify defined phases of infection by monitoring co-expression of viral and host genes, coupled with advance omics approaches, will enable a deeper understanding of the impact that viruses have on the environment.

Genomic exploration of individual giant ocean viruses

Viruses are major pathogens in all biological systems. Virus propagation and downstream analysis remains a challenge, particularly in the ocean where the majority of their microbial hosts remain recalcitrant to current culturing techniques. We used a cultivation-independent approach to isolate and sequence individual viruses. The protocol uses high-speed fluorescence-activated virus sorting flow cytometry, multiple displacement amplification (MDA), and downstream genomic sequencing. We focused on 'giant viruses' that are readily distinguishable by flow cytometry. From a single-milliliter sample of seawater collected from off the dock at Boothbay Harbor, ME, USA, we sorted almost 700 single virus particles, and subsequently focused on a detailed genome analysis of 12. A wide diversity of viruses was identified that included Iridoviridae, extended Mimiviridae and even a taxonomically novel (unresolved) giant virus. We discovered a viral metacaspase homolog in one of our sorted virus particles and discussed its implications in rewiring host metabolism to enhance infection. In addition, we demonstrated that viral metacaspases are widespread in the ocean. We also discovered a virus that contains both a reverse transcriptase and a transposase; although highly speculative, we suggest such a genetic complement would potentially allow this virus to exploit a latency propagation mechanism. Application of single virus genomics provides a powerful opportunity to circumvent cultivation of viruses, moving directly to genomic investigation of naturally occurring viruses, with the assurance that the sequence data is virus-specific, non-chimeric and contains no cellular contamination.

Coccolithoviruses: A Review of Cross-Kingdom Genomic Thievery and Metabolic Thuggery

Coccolithoviruses (Phycodnaviridae) infect and lyse the most ubiquitous and successful coccolithophorid in modern oceans, Emiliania huxleyi. So far, the genomes of 13 of these giant lytic viruses (i.e., Emiliania huxleyi viruses—EhVs) have been sequenced, assembled, and annotated. Here, we performed an in-depth comparison of their genomes to try and contextualize the ecological and evolutionary traits of these viruses. The genomes of these EhVs have from 444 to 548 coding sequences (CDSs). Presence/absence analysis of CDSs identified putative genes with particular ecological significance, namely sialidase, phosphate permease, and sphingolipid biosynthesis. The viruses clustered into distinct clades, based on their DNA polymerase gene as well as full genome comparisons. We discuss the use of such clustering and suggest that a gene-by-gene investigation approach may be more useful when the goal is to reveal differences related to functionally important genes. A multi domain “Best BLAST hit” analysis revealed that 84% of the EhV genes have closer similarities to the domain Eukarya. However, 16% of the EhV CDSs were very similar to bacterial genes, contributing to the idea that a significant portion of the gene flow in the planktonic world inter-crosses the domains of life.

Intragenus competition between coccolithoviruses: an insight on how a select few can come to dominate many

Viruses are a major cause of coccolithophore bloom demise in both temperate and sub-temperate oceanic regions. Most infection studies on coccolithoviruses have been conducted with a single virus strain, and the effect of intragenus competition by closely related coccolithoviruses has been ignored. Here we conducted combined infection experiments, infecting Emiliania huxleyi CCMP 2090 with two coccolithoviruses: EhV-86 and EhV-207 both simultaneously and independently. EhV-207 displayed a shorter lytic cycle and increased production potential than EhV-86 and was remarkably superior under competitive conditions. Although the viruses displayed identical adsorption kinetics in the first 2 h post infection, EhV-207 gained a numerical advantage as early as 8 h post infection. Quantitative polymerase chain reaction (PCR) revealed that when infecting in combination, EhV-207 was not affected by the presence of EhV-86, whereas EhV-86 was quickly out-competed, and a significant reduction in free and cell-associated EhV-86 was seen as early as 2 days after the initial infection. The observation of such clear phenotypic differences between genetically distinct, yet similar, coccolithovirus strains, by flow cytometry and quantitative real-time PCR allowed tentative links to the burgeoning genomic, transcriptomic and metabolic data to be made and the factors driving their selection, in particular to the de novo coccolithovirus-encoded sphingolipid biosynthesis pathway. This work illustrates that, even within a family, not all viruses are created equally, and the potential exists for relatively small genetic changes to infer disproportionately large competitive advantages for one coccolithovirus over another, ultimately leading to a few viruses dominating the many.

Seasonal diversity and dynamics of haptophytes in the Skagerrak, Norway, explored by high-throughput sequencing

Microalgae in the division Haptophyta play key roles in the marine ecosystem and in global biogeochemical processes. Despite their ecological importance, knowledge on seasonal dynamics, community composition and abundance at the species level is limited due to their small cell size and few morphological features visible under the light microscope. Here, we present unique data on haptophyte seasonal diversity and dynamics from two annual cycles, with the taxonomic resolution and sampling depth obtained with high-throughput sequencing. From outer Oslofjorden, S Norway, nano- and picoplanktonic samples were collected monthly for 2 years, and the haptophytes targeted by amplification of RNA/cDNA with Haptophyta-specific 18S rDNA V4 primers. We obtained 156 operational taxonomic units (OTUs), from c. 400.000 454 pyrosequencing reads, after rigorous bioinformatic filtering and clustering at 99.5%. Most OTUs represented uncultured and/or not yet 18S rDNA-sequenced species. Haptophyte OTU richness and community composition exhibited high temporal variation and significant yearly periodicity. Richness was highest in September–October (autumn) and lowest in April–May (spring). Some taxa were detected all year, such as Chrysochromulina simplex, Emiliania huxleyi and Phaeocystis cordata, whereas most calcifying coccolithophores only appeared from summer to early winter. We also revealed the seasonal dynamics of OTUs representing putative novel classes (clades HAP-3–5) or orders (clades D, E, F). Season, light and temperature accounted for 29% of the variation in OTU composition. Residual variation may be related to biotic factors, such as competition and viral infection. This study provides new, in-depth knowledge on seasonal diversity and dynamics of haptophytes in North Atlantic coastal waters.

Development of DNA mismatch repair gene, MutS, as a diagnostic marker for detection and phylogenetic analysis of algal Megaviruses

Megaviruses are generically defined as giant viruses with genomes up to 1.26Mb that infect eukaryotic unicellular protists; they are clearly delineated in DNA polymerase B phylogenetic trees; in addition, common features often include an associated virophage observed during infection; the presence of an amino acyl tRNA synthetase gene; and a nucleic acid mismatch repair protein, MutS gene. The archetypal representative of this evolving putative family is Mimivirus, an opportunistic pathogen of Acanthamoeba spp. originally thought to be a bacterium until its genome sequence was published in 2004. Subsequent analysis of marine metagenomic data revealed Megaviruses are likely ubiquitous on the surface ocean. Analysis of genome sequences of giant viruses isolated from naturally occurring marine protists such as microalgae and a microflagellate grazer, started the expansion of the Megaviridae. Here, we explored the possibility of developing Megavirus specific markers for mutS that could be used in virus molecular ecology studies. MutS is split into 15 different clades representing a wide range of cellular life, and two that contain Megaviruses, clade MutS7 and clade MutS8. We developed specific PCR primers that recognized Megavirus clade MutS8, a clade that we propose discriminates most of the algal Megaviruses. Analysis of seawater off the coast of Maine, US, revealed novel groups of algal Megaviruses that were present in all samples tested. The Megavirus clade MutS8 marker should be considered as a tool to reveal new diversity and distribution of this enigmatic group of viruses.

Functional inferences of environmental coccolithovirus biodiversity

The cosmopolitan calcifying alga Emiliania huxleyi is one of the most abundant bloom forming coccolithophore species in the oceans and plays an important role in global biogeochemical cycling. Coccolithoviruses are a major cause of coccolithophore bloom termination and have been studied in laboratory, mesocosm and open ocean studies. However, little is known about the dynamic interactions between the host and its viruses, and less is known about the natural diversity and role of functionally important genes within natural coccolithovirus communities. Here, we investigate the temporal and spatial distribution of coccolithoviruses by the use of molecular fingerprinting techniques PCR, DGGE and genomic sequencing. The natural biodiversity of the virus genes encoding the major capsid protein (MCP) and serine palmitoyltransferase (SPT) were analysed in samples obtained from the Atlantic Meridional Transect (AMT), the North Sea and the L4 site in the Western Channel Observatory. We discovered nine new coccolithovirus genotypes across the AMT and L4 site, with the majority of MCP sequences observed at the deep chlorophyll maximum layer of the sampled sites on the transect. We also found four new SPT gene variations in the North Sea and at L4. Their translated fragments and the full protein sequence of SPT from laboratory strains EhV-86 and EhV-99B1 were modelled and revealed that the theoretical fold differs among strains. Variation identified in the structural distance between the two domains of the SPT protein may have an impact on the catalytic capabilities of its active site. In summary, the combined use of 'standard' markers (i.e. MCP), in combination with metabolically relevant markers (i.e. SPT) are useful in the study of the phylogeny and functional biodiversity of coccolithoviruses, and can provide an interesting intracellular insight into the evolution of these viruses and their ability to infect and replicate within their algal hosts.

Host-virus dynamics and subcellular controls of cell fate in a natural coccolithophore population

Marine viruses are major evolutionary and biogeochemical drivers in marine microbial foodwebs. However, an in-depth understanding of the cellular mechanisms and the signal transduction pathways mediating host-virus interactions during natural bloom dynamics has remained elusive. We used field-based mesocosms to examine the "arms race" between natural populations of the coccolithophore Emiliania huxleyi and its double-stranded DNA-containing coccolithoviruses (EhVs). Specifically, we examined the dynamics of EhV infection and its regulation of cell fate over the course of bloom development and demise using a diverse suite of molecular tools and in situ fluorescent staining to target different levels of subcellular resolution. We demonstrate the concomitant induction of reactive oxygen species, caspase-specific activity, metacaspase expression, and programmed cell death in response to the accumulation of virus-derived glycosphingolipids upon infection of natural E. huxleyi populations. These subcellular responses to viral infection simultaneously resulted in the enhanced production of transparent exopolymer particles, which can facilitate aggregation and stimulate carbon flux. Our results not only corroborate the critical role for glycosphingolipids and programmed cell death in regulating E. huxleyi-EhV interactions, but also elucidate promising molecular biomarkers and lipid-based proxies for phytoplankton host-virus interactions in natural systems.