Rewiring Host Lipid Metabolism by Large Viruses Determines the Fate of Emiliania huxleyi, a Bloom-Forming Alga in the Ocean

Structure and replication cycle of a virus infecting climate-modulating alga Emiliania huxleyi

The globally distributed marine alga Emiliania huxleyi has cooling effect on the Earth's climate. The population density of E. huxleyi is restricted by Nucleocytoviricota viruses, including E. huxleyi virus 201 (EhV-201). Despite the impact of E. huxleyi viruses on the climate, there is limited information about their structure and replication. Here, we show that the dsDNA genome inside the EhV-201 virion is protected by an inner membrane, capsid, and outer membrane. EhV-201 virions infect E. huxleyi by using fivefold vertices to bind to and fuse the virus' inner membrane with the cell plasma membrane. Progeny virions assemble in the cytoplasm at the surface of endoplasmic reticulum-derived membrane segments. Genome packaging initiates synchronously with the capsid assembly and completes through an aperture in the forming capsid. The genome-filled capsids acquire an outer membrane by budding into intracellular vesicles. EhV-201 infection induces a loss of surface protective layers from E. huxleyi cells, which enables the continuous release of virions by exocytosis.

The consequences of viral infection on protists

Protists encompass a vast widely distributed group of organisms, surpassing the diversity observed in metazoans. Their diverse ecological niches and life forms are intriguing characteristics that render them valuable subjects for in-depth cell biology studies. Throughout history, viruses have played a pivotal role in elucidating complex cellular processes, particularly in the context of cellular responses to viral infections. In this comprehensive review, we provide an overview of the cellular alterations that are triggered in specific hosts following different viral infections and explore intricate biological interactions observed in experimental conditions using different host-pathogen groups.

Algal blooms in the ocean: hot spots for chemically mediated microbial interactions

The cycling of major nutrients in the ocean is affected by large-scale phytoplankton blooms, which are hot spots of microbial life. Diverse microbial interactions regulate bloom dynamics. At the single-cell level, interactions between microorganisms are mediated by small molecules in the chemical crosstalk that determines the type of interaction, ranging from mutualism to pathogenicity. Algae interact with viruses, bacteria, parasites, grazers and other algae to modulate algal cell fate, and these interactions are dependent on the environmental context. Recent advances in mass spectrometry and single-cell technologies have led to the discovery of a growing number of infochemicals - metabolites that convey information - revealing the ability of algal cells to govern biotic interactions in the ocean. The diversity of infochemicals seems to account for the specificity in cellular response during microbial communication. Given the immense impact of algal blooms on biogeochemical cycles and climate regulation, a major challenge is to elucidate how microscale interactions control the fate of carbon and the recycling of major elements in the ocean. In this Review, we discuss microbial interactions and the role of infochemicals in algal blooms. We further explore factors that can impact microbial interactions and the available tools to decipher them in the natural environment.

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.

Genomic adaptation of giant viruses in polar oceans

Despite being perennially frigid, polar oceans form an ecosystem hosting high and unique biodiversity. Various organisms show different adaptive strategies in this habitat, but how viruses adapt to this environment is largely unknown. Viruses of phyla Nucleocytoviricota and Mirusviricota are groups of eukaryote-infecting large and giant DNA viruses with genomes encoding a variety of functions. Here, by leveraging the Global Ocean Eukaryotic Viral database, we investigate the biogeography and functional repertoire of these viruses at a global scale. We first confirm the existence of an ecological barrier that clearly separates polar and nonpolar viral communities, and then demonstrate that temperature drives dramatic changes in the virus-host network at the polar-nonpolar boundary. Ancestral niche reconstruction suggests that adaptation of these viruses to polar conditions has occurred repeatedly over the course of evolution, with polar-adapted viruses in the modern ocean being scattered across their phylogeny. Numerous viral genes are specifically associated with polar adaptation, although most of their homologues are not identified as polar-adaptive genes in eukaryotes. These results suggest that giant viruses adapt to cold environments by changing their functional repertoire, and this viral evolutionary strategy is distinct from the polar adaptation strategy of their hosts.

A bacterial quorum sensing signal is a potent inhibitor of de novo pyrimidine biosynthesis in the globally abundant Emiliania huxleyi

Interactions between marine phytoplankton, viruses, and bacteria drive biogeochemical cycling, shape marine trophic structures, and impact global climate. Microbially produced compounds have emerged as key players in influencing eukaryotic organismal physiology, and in turn, remodel microbial community structure. This work aimed to reveal the molecular mechanism by which the bacterial quorum sensing molecule 2-heptyl-4-quinolone (HHQ), produced by the marine gammaproteobacterium Pseudoalteromonas spp., arrests cell division and confers protection from virus-induced mortality in the bloom-forming coccolithophore Emiliania huxleyi. Previous work has established alkylquinolones as inhibitors of dihydroorotate dehydrogenase (DHODH), a fundamental enzyme catalyzing the fourth step in pyrimidine biosynthesis and a potential antiviral drug target. An N-terminally truncated version of E. huxleyi DHODH was heterologously expressed in E. coli, purified, and kinetically characterized. Here, we show HHQ is a potent inhibitor (Ki of 2.3 nM) of E. huxleyi DHODH. E. huxleyi cells exposed to brequinar, the canonical human DHODH inhibitor, experienced immediate, yet reversible cellular arrest, an effect which mirrors HHQ-induced cellular stasis previously observed. However, brequinar treatment lacked other notable effects observed in HHQ-exposed E. huxleyi including significant changes in cell size, chlorophyll fluorescence, and protection from virus-induced lysis, indicating HHQ has additional as yet undiscovered physiological targets. Together, these results suggest a novel and intricate role of bacterial quorum sensing molecules in tripartite interdomain interactions in marine ecosystems, opening new avenues for exploring the role of microbial chemical signaling in algal bloom regulation and host-pathogen dynamics.

Virologs, viral mimicry, and virocell metabolism: the expanding scale of cellular functions encoded in the complex genomes of giant viruses

The phylum Nucleocytoviricota includes the largest and most complex viruses known. These "giant viruses" have a long evolutionary history that dates back to the early diversification of eukaryotes, and over time they have evolved elaborate strategies for manipulating the physiology of their hosts during infection. One of the most captivating of these mechanisms involves the use of genes acquired from the host-referred to here as viral homologs or "virologs"-as a means of promoting viral propagation. The best-known examples of these are involved in mimicry, in which viral machinery "imitates" immunomodulatory elements in the vertebrate defense system. But recent findings have highlighted a vast and rapidly expanding array of other virologs that include many genes not typically found in viruses, such as those involved in translation, central carbon metabolism, cytoskeletal structure, nutrient transport, vesicular trafficking, and light harvesting. Unraveling the roles of virologs during infection as well as the evolutionary pathways through which complex functional repertoires are acquired by viruses are important frontiers at the forefront of giant virus research.

Identification and Verification of Candidate miRNA Biomarkers with Application to Infection with Emiliania huxleyi Virus

The interactions of Emiliania huxleyi and its specific lytic virus (EhV) have a profound influence on marine biogeochemical carbon-sulfur cycles and play a prominent role in global climate change. MicroRNAs (miRNAs) have emerged as promising candidates with extensive diagnostic potential due to their role in virus-host interactions. However, the application of miRNA signatures as diagnostic markers in marine viral infection has made limited progress. Based on our previous small-RNA sequencing data, one host miRNA biomarker that is upregulated in early infection and seven viral miRNA biomarkers that are upregulated in late infection were identified and verified using qRT-PCR and a receiver operating characteristic curve analysis in pure culture, mixed culture, and natural seawater culture. The host ehx-miR20-5p was able to significantly differentiate infection groups from the control in the middle (24 h post-infection, hpi) and late infection (48 hpi) phases, while seven virus-derived miRNA biomarkers could diagnose the early and late stages of EhV infection. Functional enrichment analysis showed that these miRNAs participated in numerous essential metabolic pathways, including gene transcription and translation, cell division-related pathways, protein-degradation-related processes, and lipid metabolism. Additionally, a dual-luciferase reporter assay confirmed the targeted relationship between a viral ehv-miR7-5p and the host dihydroceramide desaturase gene (hDCD). This finding suggests that the virus-derived miRNA has the ability to inhibit the host sphingolipid metabolism, which is a specific characteristic of EhV infection during the late stage. Our data revealed a cluster of potential miRNA biomarkers with significant regulatory functions that could be used to diagnose EhV infection, which has implications for assessing the infectious activity of EhV in a natural marine environment.

Aerobic bacteria produce nitric oxide via denitrification and promote algal population collapse

Microbial interactions govern marine biogeochemistry. These interactions are generally considered to rely on exchange of organic molecules. Here we report on a novel inorganic route of microbial communication, showing that algal-bacterial interactions between Phaeobacter inhibens bacteria and Gephyrocapsa huxleyi algae are mediated through inorganic nitrogen exchange. Under oxygen-rich conditions, aerobic bacteria reduce algal-secreted nitrite to nitric oxide (NO) through denitrification, a well-studied anaerobic respiratory mechanism. The bacterial NO is involved in triggering a cascade in algae akin to programmed cell death. During death, algae further generate NO, thereby propagating the signal in the algal population. Eventually, the algal population collapses, similar to the sudden demise of oceanic algal blooms. Our study suggests that the exchange of inorganic nitrogen species in oxygenated environments is a potentially significant route of microbial communication within and across kingdoms.

Virus infection of phytoplankton increases average molar mass and reduces hygroscopicity of aerosolized organic matter

Viral infection of phytoplankton is a pervasive mechanism of cell death and bloom termination, which leads to the production of dissolved and colloidal organic matter that can be aerosolized into the atmosphere. Earth-observing satellites can track the growth and death of phytoplankton blooms on weekly time scales but the impact of viral infection on the cloud forming potential of associated aerosols is largely unknown. Here, we determine the influence of viral-derived organic matter, purified viruses, and marine hydrogels on the cloud condensation nuclei activity of their aerosolized solutions, compared to organic exudates from healthy phytoplankton. Dissolved organic material derived from exponentially growing and infected cells of well-characterized eukaryotic phytoplankton host-virus systems, including viruses from diatoms, coccolithophores and chlorophytes, was concentrated, desalted, and nebulized to form aerosol particles composed of primarily of organic matter. Aerosols from infected phytoplankton cultures resulted in an increase in critical activation diameter and average molar mass in three out of five combinations evaluated, along with a decrease in organic kappa (hygroscopicity) compared to healthy cultures and seawater controls. The infected samples also displayed evidence of increased surface tension depression at realistic cloud water vapor supersaturations. Amending the samples with xanthan gum to simulate marine hydrogels increased variability in organic kappa and surface tension in aerosols with high organic to salt ratios. Our findings suggest that the pulses of increased dissolved organic matter associated with viral infection in surface waters may increase the molar mass of dissolved organic compounds relative to surface waters occupied by healthy phytoplankton or low phytoplankton biomass.

Functional Profiling and Evolutionary Analysis of a Marine Microalgal Virus Pangenome

Phycodnaviridae are large double-stranded DNA viruses, which facilitate studies of host-virus interactions and co-evolution due to their prominence in algal infection and their role in the life cycle of algal blooms. However, the genomic interpretation of these viruses is hampered by a lack of functional information, stemming from the surprising number of hypothetical genes of unknown function. It is also unclear how many of these genes are widely shared within the clade. Using one of the most extensively characterized genera, Coccolithovirus, as a case study, we combined pangenome analysis, multiple functional annotation tools, AlphaFold structural modeling, and literature analysis to compare the core and accessory pangenome and assess support for novel functional predictions. We determined that the Coccolithovirus pangenome shares 30% of its genes with all 14 strains, making up the core. Notably, 34% of its genes were found in at most three strains. Core genes were enriched in early expression based on a transcriptomic dataset of Coccolithovirus EhV-201 algal infection, were more likely to be similar to host proteins than the non-core set, and were more likely to be involved in vital functions such as replication, recombination, and repair. In addition, we generated and collated annotations for the EhV representative EhV-86 from 12 different annotation sources, building up information for 142 previously hypothetical and putative membrane proteins. AlphaFold was further able to predict structures for 204 EhV-86 proteins with a modelling accuracy of good-high. These functional clues, combined with generated AlphaFold structures, provide a foundational framework for the future characterization of this model genus (and other giant viruses) and a further look into the evolution of the Coccolithovirus proteome.

Sugarcane mosaic virus orchestrates the lactate fermentation pathway to support its successful infection

Viruses often establish their own infection by altering host metabolism. How viruses co-opt plant metabolism to support their successful infection remains an open question. Here, we used untargeted metabolomics to reveal that lactate accumulates immediately before and after robust sugarcane mosaic virus (SCMV) infection. Induction of lactate-involved anaerobic glycolysis is beneficial to SCMV infection. The enzyme activity and transcriptional levels of lactate dehydrogenase (LDH) were up-regulated by SCMV infection, and LDH is essential for robust SCMV infection. Moreover, LDH relocates in viral replicase complexes (VRCs) by interacting with SCMV-encoded 6K2 protein, a key protein responsible for inducing VRCs. Additionally, lactate could promote SCMV infection by suppressing plant defense responses. Taken together, we have revealed a viral strategy to manipulate host metabolism to support replication compartment but also depress the defense response during the process of infection.

Diversity in Infection Specificity between the Bloom-forming Microalga Heterosigma akashiwo and Its dsDNA Virus, Heterosigma akashiwo Virus

Heterosigma akashiwo virus (HaV) is a dsDNA virus that infects the bloom-forming raphidoflagellate Heterosigma akashiwo. Both the host and its virus are phenotypically diverse in terms of infection specificity. Their relationships have been examined based on the occurrence or absence of algal lysis following virus inoculation; however, variations in the strain-level host-virus relationship regarding infectivity and lysis rates remain unclear. Therefore, we performed a series of cross-infectivity tests using 60 H. akashiwo and 22 HaV strains isolated from the coastal waters of western Japan. The host strains were divided into 5 different groups and viruses into 4 groups. Using a representative strain from each group, algal lysis was observed in 14 of the (5×4=) 20 host-virus combinations; the concentration of infectious units in each HaV suspension was then assessed using the most probable number (MPN) assay on the five host strains. Virus titers ranged between 1.1×10¹ and 2.1×10⁷ infectious units mL^-1; the titer of each viral lysate was differently estimated using distinct H. akashiwo strains as hosts. These results suggest that (1) a clonal viral lysate comprises virions with different intraspecific infection specificities and/or (2) the efficiency and error rates of each intracellular replication process vary in each host-virus combination.

MicroRNA-mediated regulation of lipid metabolism in virus-infected Emiliania huxleyi

The interactions between Emiliania huxleyi and E. huxleyi virus (EhV) regulate marine carbon and sulfur biogeochemical cycles and play a prominent role in global climate change. As a large DNA virus, EhV has developed a novel "virocell metabolism" model to meet its high metabolic needs. Although it has been widely demonstrated that EhV infection can profoundly rewire lipid metabolism, the epigenetic regulatory mechanisms of lipid metabolism are still obscure. MicroRNAs (miRNAs) can regulate biological pathways by targeting hub genes in the metabolic processes. In this study, the transcriptome, lipidome, and miRNAome were applied to investigate the epigenetic regulation of lipid metabolism in E. huxleyi cells during a detailed time course of viral infection. Combined transcriptomic, lipidomic, and physiological experiments revealed reprogrammed lipid metabolism, along with mitochondrial dysfunction and calcium influx through the cell membrane. A total of 69 host miRNAs (including 1 known miRNA) and 7 viral miRNAs were identified, 27 of which were differentially expressed. Bioinformatic prediction revealed that miRNAs involved in the regulation of lipid metabolism and a dual-luciferase reporter assay suggested that phosphatidylinositol 3-kinase (PI3K) gene might be a target of ehx-miR5. Further qPCR and western blot analysis showed a significant negative correlation between the expression of ehx-miR5 and its target gene PI3K, along with the lower activity of its downstream components (p-Akt, p-TOR, SREBP), indicating that lipid metabolism might be regulated by ehx-miR5 through the PI3K-Akt-TOR signaling pathway. Our findings reveal several novel mechanisms of viral strategies to manipulate host lipid metabolism and provide evidence that ehx-miR5 negatively modulates the expression of PI3K and disturbs lipid metabolism in the interactions between E. huxleyi and EhV.

Giant virus biology and diversity in the era of genome-resolved metagenomics

The discovery of giant viruses, with capsids as large as some bacteria, megabase-range genomes and a variety of traits typically found only in cellular organisms, was one of the most remarkable breakthroughs in biology. Until recently, most of our knowledge of giant viruses came from ~100 species-level isolates for which genome sequences were available. However, these isolates were primarily derived from laboratory-based co-cultivation with few cultured protists and algae and, thus, did not reflect the true diversity of giant viruses. Although virus co-cultures enabled valuable insights into giant virus biology, many questions regarding their origin, evolution and ecological importance remain unanswered. With advances in sequencing technologies and bioinformatics, our understanding of giant viruses has drastically expanded. In this Review, we summarize our understanding of giant virus diversity and biology based on viral isolates as laboratory cultivation has enabled extensive insights into viral morphology and infection strategies. We then explore how cultivation-independent approaches have heightened our understanding of the coding potential and diversity of the Nucleocytoviricota. We discuss how metagenomics has revolutionized our perspective of giant viruses by revealing their distribution across our planet's biomes, where they impact the biology and ecology of a wide range of eukaryotic hosts and ultimately affect global nutrient cycles.

Phytosterol metabolism in plant positive-strand RNA virus replication

The genome of most plant viruses consists of a single positive-strand of RNA (+ ssRNA). Successful replication of these viruses is fully dependent on the endomembrane system of the infected cells, which experiences a massive proliferation and a profound reshaping that enables assembly of the macromolecular complexes where virus genome replication occurs. Assembly of these viral replicase complexes (VRCs) requires a highly orchestrated interplay of multiple virus and co-opted host cell factors to create an optimal microenvironment for efficient assembly and functioning of the virus genome replication machinery. It is now widely accepted that VRC formation involves the recruitment of high levels of sterols, but the specific role of these essential components of cell membranes and the precise molecular mechanisms underlying sterol enrichment at VRCs are still poorly known. In this review, we intend to summarize the most relevant knowledge on the role of sterols in ( +)ssRNA virus replication and discuss the potential of manipulating the plant sterol pathway to help plants fight these infectious agents.

Temporal and Spatial Signaling Mediating the Balance of the Plankton Microbiome

The annual patterns of plankton succession in the ocean determine ecological and biogeochemical cycles. The temporally fluctuating interplay between photosynthetic eukaryotes and the associated microbiota balances the composition of aquatic planktonic ecosystems. In addition to nutrients and abiotic factors, chemical signaling determines the outcome of interactions between phytoplankton and their associated microbiomes. Chemical mediators control essential processes, such as the development of key morphological, physiological, behavioral, and life-history traits during algal growth. These molecules thus impact species succession and community composition across time and space in processes that are highlighted in this review. We focus on spatial, seasonal, and physiological dynamics that occur during the early association of algae with bacteria, the exponential growth of a bloom, and its decline and recycling. We also discuss how patterns from field data and global surveys might be linked to the actions of metabolic markers in natural phytoplankton assemblages.

Resolving the microalgal gene landscape at the strain level: A novel hybrid transcriptome of Emiliania huxleyi CCMP3266

Microalgae are key ecological players with a complex evolutionary history. Genomic diversity, in addition to limited availability of high-quality genomes, challenge studies that aim to elucidate molecular mechanisms underlying microalgal ecophysiology. Here, we present a novel and comprehensive transcriptomic hybrid approach to generate a reference for genetic analyses, and resolve the microalgal gene landscape at the strain level. The approach is demonstrated for a strain of the coccolithophore microalga Emiliania huxleyi, which is a species complex with considerable genome variability. The investigated strain is commonly studied as a model for algal-bacterial interactions, and was therefore sequenced in the presence of bacteria to elicit the expression of interaction-relevant genes. We applied complementary PacBio Iso-Seq full-length cDNA, and poly(A)-independent Illumina total RNA sequencing, which resulted in a de novo assembled, near complete hybrid transcriptome. In particular, hybrid sequencing improved the reconstruction of long transcripts and increased the recovery of full-length transcript isoforms. To use the resulting hybrid transcriptome as a reference for genetic analyses, we demonstrate a method that collapses the transcriptome into a genome-like dataset, termed "synthetic genome" (sGenome). We used the sGenome as a reference to visually confirm the robustness of the CCMP3266 gene assembly, to conduct differential gene expression analysis, and to characterize novel E. huxleyi genes. The newly-identified genes contribute to our understanding of E. huxleyi genome diversification, and are predicted to play a role in microbial interactions. Our transcriptomic toolkit can be implemented in various microalgae to facilitate mechanistic studies on microalgal diversity and ecology. Importance Microalgae are key players in the ecology and biogeochemistry of our oceans. Efforts to implement genomic and transcriptomic tools in laboratory studies involving microalgae suffer from the lack of published genomes. In the case of coccolithophore microalgae, the problem has long been recognized; the model species Emiliania huxleyi is a species complex with genomes composed of a core, and a large variable portion. To study the role of the variable portion in niche adaptation, and specifically in microbial interactions, strain-specific genetic information is required. Here we present a novel transcriptomic hybrid approach, and generated strain-specific genome-like information. We demonstrate our approach on an E. huxleyi strain that is co-cultivated with bacteria. By constructing a "synthetic genome", we generated comprehensive gene annotations that enabled accurate analyses of gene expression patterns. Importantly, we unveiled novel genes in the variable portion of E. huxleyi that play putative roles in microbial interactions.

An Emiliania huxleyi pan-transcriptome reveals basal strain specificity in gene expression patterns

Emiliania huxleyi is a cosmopolitan coccolithophore widespread in temperate oceans. This unicellular photoautotroph forms massive recurring blooms that play an important role in large biogeochemical cycles of carbon and sulfur, which play a role in climate change. The mechanism of bloom formation and demise, controlled by giant viruses that routinely infect these blooms, is poorly understood. We generated a pan-transcriptome of E. huxleyi, derived from three strains with different susceptibility to viral infection. Expression profiling of E. huxleyi sensitive and resistant strains showed major basal differences, including many genes that are induced upon viral infection. This suggests that basal gene expression can affect the host metabolic state and the susceptibility of E. huxleyi to viruses. Due to its ecological importance, the pan-transcriptome and its protein translation, applicable to many E. huxleyi strains, is a powerful resource for investigation of eukaryotic microbial communities.

Adaptive evolution of viruses infecting marine microalgae (haptophytes), from acute infections to stable coexistence

Collectively known as phytoplankton, photosynthetic microbes form the base of the marine food web, and account for up to half of the primary production on Earth. Haptophytes are key components of this phytoplankton community, playing important roles both as primary producers and as mixotrophs that graze on bacteria and protists. Viruses influence the ecology and diversity of phytoplankton in the ocean, with the majority of microalgae-virus interactions described as 'boom and bust' dynamics, which are characteristic of acute virus-host systems. Most haptophytes are, however, part of highly diverse communities and occur at low densities, decreasing their chance of being infected by viruses with high host specificity. Viruses infecting these microalgae have been isolated in the laboratory, and there are several characteristics that distinguish them from acute viruses infecting bloom-forming haptophytes. Herein we synthesise what is known of viruses infecting haptophyte hosts in the ocean, discuss the adaptive evolution of haptophyte-infecting viruses -from those that cause acute infections to those that stably coexist with their host - and identify traits of importance for successful survival in the ocean.

Ecological significance of extracellular vesicles in modulating host-virus interactions during algal blooms

Extracellular vesicles are produced by organisms from all kingdoms and serve a myriad of functions, many of which involve cell-cell signaling, especially during stress conditions and host-pathogen interactions. In the marine environment, communication between microorganisms can shape trophic level interactions and population succession, yet we know very little about the involvement of vesicles in these processes. In a previous study, we showed that vesicles produced during viral infection by the ecologically important model alga Emiliania huxleyi, could act as a pro-viral signal, by expediting infection and enhancing the half-life of the virus in the extracellular milieu. Here, we expand our laboratory findings and show the effect of vesicles on natural populations of E. huxleyi in a mesocosm setting. We profile the small-RNA (sRNA) cargo of vesicles that were produced by E. huxleyi during bloom succession, and show that vesicles applied to natural assemblages expedite viral infection and prolong the half-life of this major mortality agent of E. huxleyi. We subsequently reveal that exposure of the natural assemblage to E. huxleyi-derived vesicles modulates not only host-virus dynamics, but also other components of the microbial food webs, thus emphasizing the importance of extracellular vesicles to microbial interactions in the marine environment.

Viral infection of algal blooms leaves a unique metabolic footprint on the dissolved organic matter in the ocean

Algal blooms are hotspots of primary production in the ocean, forming the basis of the marine food web and fueling the dissolved organic matter (DOM) pool. Viruses are key players in controlling algal demise, thereby diverting biomass from higher trophic levels to the DOM pool, a process termed the "viral shunt." To decode the metabolic footprint of the viral shunt in the environment, we induced a bloom of Emiliania huxleyi and followed its succession using untargeted exometabolomics. We show that bloom succession induces dynamic changes in the exometabolic landscape. We found a set of chlorine-iodine-containing metabolites that were induced by viral infection and released during bloom demise. These metabolites were further detected in virus-infected oceanic E. huxleyi blooms. Therefore, we propose that halogenation with both chlorine and iodine is a distinct hallmark of the virus-induced DOM of E. huxleyi, providing insights into the metabolic consequences of the viral shunt.

Metatranscriptome Library Preparation Influences Analyses of Viral Community Activity During a Brown Tide Bloom

There is growing interest in the use of metatranscriptomics to study virus community dynamics. We used RNA samples collected from harmful brown tides caused by the eukaryotic alga Aureococcus anophagefferens within New York (United States) estuaries and in the process observed how preprocessing of libraries by either selection for polyadenylation or reduction in ribosomal RNA (rRNA) influenced virus community analyses. As expected, more reads mapped to the A. anophagefferens genome in polyadenylation-selected libraries compared to the rRNA-reduced libraries, with reads mapped in each sample correlating to one another regardless of preprocessing of libraries. Yet, this trend was not seen for reads mapping to the Aureococcus anophagefferens Virus (AaV), where significantly more reads (approximately two orders of magnitude) were mapped to the AaV genome in the rRNA-reduced libraries. In the rRNA-reduced libraries, there was a strong and significant correlation between reads mappings to AaV and A. anophagefferens. Overall, polyadenylation-selected libraries produced fewer viral contigs, fewer reads mapped to viral contigs, and different proportions across viral realms and families, compared to their rRNA-reduced pairs. This study provides evidence that libraries generated by rRNA reduction and not selected for polyadenylation are more appropriate for quantitative characterization of viral communities in aquatic ecosystems by metatranscriptomics.

Phytosterol Profiles, Genomes and Enzymes - An Overview

The remarkable diversity of sterol biosynthetic capacities described in living organisms is enriched at a fast pace by a growing number of sequenced genomes. Whereas analytical chemistry has produced a wealth of sterol profiles of species in diverse taxonomic groups including seed and non-seed plants, algae, phytoplanktonic species and other unicellular eukaryotes, functional assays and validation of candidate genes unveils new enzymes and new pathways besides canonical biosynthetic schemes. An overview of the current landscape of sterol pathways in the tree of life is tentatively assembled in a series of sterolotypes that encompass major groups and provides also peculiar features of sterol profiles in bacteria, fungi, plants, and algae.

A persistent giant algal virus, with a unique morphology, encodes an unprecedented number of genes involved in energy metabolism

Viruses have long been viewed as entities possessing extremely limited metabolic capacities. Over the last decade, however, this view has been challenged, as metabolic genes have been identified in viruses possessing large genomes and virions-the synthesis of which is energetically demanding. Here, we unveil peculiar phenotypic and genomic features of Prymnesium kappa virus RF01 (PkV RF01), a giant virus of the Mimiviridae family. We found that this virus encodes an unprecedented number of proteins involved in energy metabolism, such as all four succinate dehydrogenase (SDH) subunits (A-D) as well as key enzymes in the β-oxidation pathway. The SDHA gene was transcribed upon infection, indicating that the viral SDH is actively used by the virus- potentially to modulate its host's energy metabolism. We detected orthologous SDHA and SDHB genes in numerous genome fragments from uncultivated marine Mimiviridae viruses, which suggests that the viral SDH is widespread in oceans. PkV RF01 was less virulent compared with other cultured prymnesioviruses, a phenomenon possibly linked to the metabolic capacity of this virus and suggestive of relatively long co-evolution with its hosts. It also has a unique morphology, compared to other characterized viruses in the Mimiviridae family. Finally, we found that PkV RF01 is the only alga-infecting Mimiviridae virus encoding two aminoacyl-tRNA synthetases and enzymes corresponding to an entire base-excision repair pathway, as seen in heterotroph-infecting Mimiviridae These Mimiviridae encoded-enzymes were found to be monophyletic and branching at the root of the eukaryotic tree of life. This placement suggests that the last common ancestor of Mimiviridae was endowed with a large, complex genome prior to the divergence of known extant eukaryotes.IMPORTANCE Viruses on Earth are tremendously diverse in terms of morphology, functionality, and genomic composition. Over the last decade, the conceptual gap separating viruses and cellular life has tightened because of the detection of metabolic genes in viral genomes that express complex virus phenotypes upon infection. Here, we describe Prymnesium kappa virus RF01, a large alga-infecting virus with a unique morphology, an atypical infection profile, and an unprecedented number of genes involved in energy metabolism (such as the tricarboxylic (TCA) cycle and the β-oxidation pathway). Moreover, we show that the gene corresponding to one of these enzymes (the succinate dehydrogenase subunit A) is transcribed during infection and is widespread among marine viruses. This discovery provides evidence that a virus has the potential to actively regulate energy metabolism with its own gene.

Lipidome disturbances in preadipocyte differentiation associated with bisphenol A and replacement bisphenol S exposure

Bisphenol S (BPS) is the major substitute for the production of bisphenol A (BPA)-free products and detected in both food and environment. Although the relationship between BPA exposure and increased risk of obesity and diabetes has been noted, the potential influence of BPS is not fully understood. Herein, a non-targeted lipidomic study was performed to explore BPA/BPS exposure actions using the 3T3-L1 preadipocyte differentiation model, and revealed the comprehensive lipidome disturbance induced by either BPA or BPS exposure at different doses of 0.01, 1 and 100 μM. BPA was more potent than BPS in disturbance of lipid metabolism. A considerable similarity of BPS exposure to BPA was discovered. The key lipid remodeling in response to exposure was found to involve the cardiolipins, phosphatidylglycerols and fatty acids metabolic pathways, providing novel clues of potential mechanism in which both BPA and BPS exposure could be associated with increased risk of insulin resistance. Our study supplies the perspective into the lipidome response to environmental stress induced by BPA/BPS, and shows that BPA-free products are not necessarily safer. Substitution of BPA by its structural analog BPS should be therefore performed with caution.

Validation of superior reference genes for qRT-PCR and Western blot analyses in marine Emiliania huxleyi-virus model system

AIMS

To search for a set of reference genes for reliable gene expression analysis in the globally important marine coccolithophore Emiliania huxleyi-virus model system.

METHODS AND RESULTS

Fifteen housekeeping genes (CDKA, CYP15, EFG3, POLAI, RPL30, RPL13, SAMS, COX1, GPB1-2, HSP90, TUA, TUB, UBA1, CAM3 and GAPDH) were evaluated for their stability as potential reference genes for qRT-PCR using ΔCt, geNorm, NormFinder, Bestkeeper and RefFinder software. CDKA, TUA and TUB genes were tested as loading controls for Western blot in the same sample panel. Additionally, target genes associated with cell apoptosis, that is metacaspase genes, were applied to validate the selection of reference genes. The analysis results demonstrated that putative housekeeping genes exhibited significant variations in both mRNA and protein content during virus infection. After a comprehensive analysis with all the algorithms, CDKA and GAPDH were recommended as the most stable reference genes for E huxleyi virus (EhV) infection treatments. For Western blot, significant variation was seen for TUA and TUB, whereas CDKA was stably expressed, consistent with the results of qRT-PCR.

CONCLUSIONS

CDKA and GAPDH are the best choice for gene and protein expression analysis than the other candidate reference genes under EhV infection conditions.

SIGNIFICANCE AND IMPACT OF THE STUDY

The stable internal control genes identified in this work will help to improve the accuracy and reliability of gene expression analysis and gain insight into complex E. huxleyi-EhV interaction regulatory networks.

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.

High CO2 concentration and iron availability determine the metabolic inventory in an Emiliania huxleyi-dominated phytoplankton community

Ocean acidification (OA), a consequence of anthropogenic carbon dioxide (CO₂ ) emissions, strongly impacts marine ecosystems. OA also influences iron (Fe) solubility, affecting biogeochemical and ecological processes. We investigated the interactive effects of CO₂ and Fe availability on the metabolome response of a natural phytoplankton community. Using mesocosms we exposed phytoplankton to ambient (390 μatm) or future CO₂ levels predicted for the year 2100 (900 μatm), combined with ambient (4.5 nM) or high (12 nM) dissolved iron (dFe). By integrating over the whole phytoplankton community, we assigned functional changes based on altered metabolite concentrations. Our study revealed the complexity of phytoplankton metabolism. Metabolic profiles showed three stages in response to treatments and phytoplankton dynamics. Metabolome changes were related to the plankton group contributing respective metabolites, explaining bloom decline and community succession. CO₂ and Fe affected metabolic profiles. Most saccharides, fatty acids, amino acids and many sterols significantly correlated with the high dFe treatment at ambient pCO₂ . High CO₂ lowered the abundance of many metabolites irrespective of Fe. However, sugar alcohols accumulated, indicating potential stress. We demonstrate that not only altered species composition but also changes in the metabolic landscape affecting the plankton community may change as a consequence of future high-CO₂ oceans.

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.

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.

Internal Nitrogen Pools Shape the Infection of Aureococcus anophagefferens CCMP 1984 by a Giant Virus

The pelagophyte Aureococcus anophagefferens blooms annually in shallow bays around the world, where it is hypothesized to outcompete other phytoplankton in part by using alternative nitrogen sources. The high proportion of natural populations that are infected during the late stages of the bloom suggest viruses cause bloom collapse. We hypothesized that the Aureococcus anophagefferens Virus (AaV) infection cycle would be negatively influenced in cultures acclimated to decreasing external nitrogen conditions, but that the real-time external nitrogen concentration would not influence the infection cycle. Cultures acclimated in NO 3 - concentrations (0.0147 mM; N:P = 0.1225) that showed reduced end point cell abundances, forward scatter (a proxy for size) and red fluorescence (a proxy for chlorophyll a), also produced fewer viruses per cell at a slower rate. Decreasing the external concentration of nitrogen post infection did not alter burst size or time to lysis. These data suggest that the nitrogen used for new viral progeny is present within host cells at the time of infection. Flow cytometric data of an infection cycle showed a reduction in red fluorescence around twelve hours post infection, consistent with degradation of nitrogen-rich chloroplasts during the infection cycle. Using cell and virus quota estimates, we determined that A. anophagefferens cells had sufficient nitrogen and carbon for the lower ranges of burst sizes determined but did not contain enough phosphorous. Consistent with this observation, expression of nitrate and sugar transporters did not increase in the publicly available transcriptome data of the infection cycle, while several phosphorus transporters were. Our data demonstrate that dynamics of viruses infecting Aureococcus over the course of a bloom is dictated by the host cell state upon infection, which is set a priori by external nutrient supplies.

Extrachromosomal Genetic Engineering of the Marine Diatom Phaeodactylum tricornutum Enables the Heterologous Production of Monoterpenoids

Geraniol is a commercially relevant plant-derived monoterpenoid that is a main component of rose essential oil and used as insect repellent. Geraniol is also a key intermediate compound in the biosynthesis of the monoterpenoid indole alkaloids (MIAs), a group of over 2000 compounds that include high-value pharmaceuticals. As plants naturally produce extremely small amounts of these molecules and their chemical synthesis is complex, industrially sourcing these compounds is costly and inefficient. Hence, microbial hosts suitable to produce MIA precursors through synthetic biology and metabolic engineering are currently being sought. Here, we evaluated the suitability of a eukaryotic microalga, the marine diatom Phaeodactylum tricornutum, for the heterologous production of monoterpenoids. Profiling of endogenous metabolism revealed that P. tricornutum, unlike other microbes employed for industrial production of terpenoids, accumulates free pools of the precursor geranyl diphosphate. To evaluate the potential for larger synthetic biology applications, we engineered P. tricornutum through extrachromosomal, episome-based expression, for the heterologous biosynthesis of the MIA intermediate geraniol. By profiling the production of geraniol resulting from various genetic and cultivation arrangements, P. tricornutum reached the maximum geraniol titer of 0.309 mg/L in phototrophic conditions. This work provides (i) a detailed analysis of P. tricornutum endogenous terpenoid metabolism, (ii) a successful demonstration of extrachromosomal expression for metabolic pathway engineering with potential gene-stacking applications, and (iii) a convincing proof-of-concept of the suitability of P. tricornutum as a novel production platform for heterologous monoterpenoids, with potential for complex pathway engineering aimed at the heterologous production of MIAs.

Modeling the avoidance behavior of zooplankton on phytoplankton infected by free viruses

In any ecosystem, chaotic situations may arise from equilibrium state for different reasons. To overcome these chaotic situations, sometimes the system itself exhibits some mechanisms of self-adaptability. In this paper, we explore an eco-epidemiological model consisting of three aquatic groups: phytoplankton, zooplankton, and marine free viruses. We assume that the phytoplankton population is infected by external free viruses and zooplankton get affected on consumption of infected phytoplankton; also, the infected phytoplankton do not compete for resources with the susceptible one. In addition, we model a mechanism by which zooplankton recognize and avoid infected phytoplankton, at least when susceptible phytoplankton are present. The zooplankton extinction chance increases on increasing the force of infection or decreasing the intensity of avoidance. Further, when the viral infection triggers chaotic dynamics, high zooplankton avoidance intensity can stabilize again the system. Interestingly, for high avoidance intensity, nutrient enrichment has a destabilizing effect on the system dynamics, which is in line with the paradox of enrichment. Global sensitivity analysis helps to identify the most significant parameters that reduce the infected phytoplankton in the system. Finally, we compare the dynamics of the system by allowing the infected phytoplankton also to share resources with the susceptible phytoplankton. A gradual increase of the virus replication factor turns the system dynamics from chaos to doubling state to limit cycle to stable state and the system finally settles down to the zooplankton-free equilibrium point. Moreover, on increasing the intensity of avoidance, the system shows a transcritical bifurcation from the zooplankton-free equilibrium to the coexistence steady state and remains stable thereafter.

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.

The oomycete Lagenisma coscinodisci hijacks host alkaloid synthesis during infection of a marine diatom

Flagellated oomycetes frequently infect unicellular algae, thus limiting their proliferation. Here we show that the marine oomycete Lagenisma coscinodisci rewires the metabolome of the bloom-forming diatom Coscinodiscus granii, thereby promoting infection success. The algal alkaloids β-carboline and 4-carboxy-2,3,4,9-tetrahydro-1H-β-carboline are induced during infection. Single-cell profiling with AP-MALDI-MS and confocal laser scanning microscopy reveals that algal carbolines accumulate in the reproductive form of the parasite. The compounds arrest the algal cell division, increase the infection rate and induce plasmolysis in the host. Our results indicate that the oomycete manipulates the host metabolome to support its own multiplication.

Flagellated oomycetes frequently infect unicellular algae, thus limiting their proliferation. Here, the authors show that an oomycete rewires the metabolome of a marine bloom-forming diatom, thereby promoting infection success.

Emerging lipidome patterns associated with marine Emiliania huxleyi-virus model system

Emiliania huxleyi (Coccolithophore) plays a prominent role in the global carbon cycle and in climate processes. The annual collapse of massive E. huxleyi blooms in the marine environment has been shown to be frequently linked to viral control. These host-virus interactions shape the evolution and dynamics of oceanic microscale ecosystems, yet we still understand little of the molecular mechanism of these virus-mediated processes. Here, we present a detailed characterization of the lipidome of E. huxleyi BOF92 strain, both of uninfected cells and those infected with its specific lytic virus EhV-99B1. Non-targeted lipidomics analysis was performed in order to evaluate the dynamic alterations underlying virus-induced metabolic remodeling. The host lipidome (both lipid content and composition) significantly changed in response to the viral infection. The most statistically significant differential lipids were screened as potential biomarkers for assessing E. huxleyi population sensitivity to EhV infection. Our results reveal that the remodeling of lipid metabolism that underlies the pathogenesis of this infection primarily involved sphingolipid, glycerolipid and fatty acid metabolic pathways. Our study provides insights into how viruses shape their hosts metabolism to support their unique life cycle and a lipid-based chemical arms race during host-virus dynamic interactions in a marine environment.

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.

Physiological responses of Oxyrrhis marina to a diet of virally infected Emiliania huxleyi

The coccolithophore Emiliania huxleyi forms some of the largest phytoplankton blooms in the ocean. The rapid demise of these blooms has been linked to viral infections. E. huxleyi abundance, distribution, and nutritional status make them an important food source for the heterotrophic protists which are classified as microzooplankton in marine food webs. In this study we investigated the fate of E. huxleyi (CCMP 374) infected with virus strain EhV-86 in a simple predator-prey interaction. The ingestion rates of Oxyrrhis marina were significantly lower (between 26.9 and 50.4%) when fed virus-infected E. huxleyi cells compared to non-infected cells. Despite the lower ingestion rates, O. marina showed significantly higher growth rates (between 30 and 91.3%) when fed infected E. huxleyi cells, suggesting higher nutritional value and/or greater assimilation of infected E. huxleyi cells. No significant differences were found in O. marina cell volumes or fatty acids profiles. These results show that virally infected E. huxleyi support higher growth rates of single celled heterotrophs and in addition to the “viral shunt” hypothesis, viral infections may also divert more carbon to mesozooplankton grazers.

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.

Phaeobacter inhibens induces apoptosis-like programmed cell death in calcifying Emiliania huxleyi

The model coccolithophore, Emiliania huxleyi, forms expansive blooms dominated by the calcifying cell type, which produce calcite scales called coccoliths. Blooms last several weeks, after which the calcified algal cells rapidly die, descending into the deep ocean. E. huxleyi bloom collapse is attributed to E. huxleyi viruses (EhVs) that infect and kill calcifying cells, while other E. huxleyi pathogens, such as bacteria belonging to the roseobacter clade, are overlooked. EhVs kill calcifying E. huxleyi by inducing production of bioactive viral-glycosphingolipids (vGSLs), which trigger algal programmed cell death (PCD). The roseobacter Phaeobacter inhibens was recently shown to interact with and kill the calcifying cell type of E. huxleyi, but the mechanism of algal death remains unelucidated. Here we demonstrate that P. inhibens kills calcifying E. huxleyi by inducing a highly specific type of PCD called apoptosis-like-PCD (AL-PCD). Host death can successfully be abolished in the presence of a pan-caspase inhibitor, which prevents the activation of caspase-like molecules. This finding differentiates P. inhibens and EhV pathogenesis of E. huxleyi, by demonstrating that bacterial-induced AL-PCD requires active caspase-like molecules, while the viral pathogen does not. This is the first demonstration of a bacterium inducing AL-PCD in an algal host as a killing mechanism.

In plaque-mass spectrometry imaging of a bloom-forming alga during viral infection reveals a metabolic shift towards odd-chain fatty acid lipids

Tapping into the metabolic cross-talk between a host and its virus can reveal unique strategies employed during infection. Viral infection is a dynamic process that generates an evolving metabolic landscape. Gaining a continuous view into the infection process is highly challenging and is limited by current metabolomics approaches, which typically measure the average of the entire population at various stages of infection. Here, we took an innovative approach to study the metabolic basis of host-virus interactions between the bloom-forming alga Emiliania huxleyi and its specific virus. We combined a classical method in virology, the plaque assay, with advanced mass spectrometry imaging (MSI), an approach we termed ‘in plaque-MSI’. Taking advantage of the spatial characteristics of the plaque, we mapped the metabolic landscape induced during infection in a high spatiotemporal resolution, unfolding the infection process in a continuous manner. Further unsupervised spatially-aware clustering, combined with known lipid biomarkers, revealed a systematic metabolic shift during infection towards lipids containing the odd-chain fatty acid pentadecanoic acid (C15:0). Applying ‘in plaque-MSI’ might facilitate the discovery of bioactive compounds that mediate the chemical arms race of host-virus interactions in diverse model systems.

Light regulation of coccolithophore host-virus interactions

Viruses that infect photoautotrophs have a fundamental relationship with light, given the need for host resources. We investigated the role of light on Coccolithovirus (EhV) infection of the globally distributed coccolithophore, Emiliania huxleyi. Light was required for EhV adsorption, and viral production was highest when host cultures were maintained in continuous light or at irradiance levels of 150-300 μmol m^-2  s^-1 . During the early stages of infection, photosynthetic electron transport remained high, while RuBisCO expression decreased concomitant with an induction of the pentose phosphate pathway, the primary source of de novo nucleotides. A mathematical model developed and fitted to the laboratory data supported the hypothesis that EhV replication was controlled by a trade-off between host nucleotide recycling and de novo synthesis, and that photoperiod and photon flux could toggle this switch. Laboratory results supported field observations that light was the most robust driver of EhV replication within E. huxleyi populations collected across a 2000 nautical mile transect in the North Atlantic. Collectively, these findings demonstrate that light can drive host-virus interactions through a mechanistic interplay between host metabolic processes, which serve to structure infection and phytoplankton mortality in the upper ocean.

Nitric oxide production and antioxidant function during viral infection of the coccolithophore Emiliania huxleyi

Emiliania huxleyi is a globally important marine phytoplankton that is routinely infected by viruses. Understanding the controls on the growth and demise of E. huxleyi blooms is essential for predicting the biogeochemical fate of their organic carbon and nutrients. In this study, we show that the production of nitric oxide (NO), a gaseous, membrane-permeable free radical, is a hallmark of early-stage lytic infection in E. huxleyi by Coccolithoviruses, both in culture and in natural populations in the North Atlantic. Enhanced NO production was detected both intra- and extra-cellularly in laboratory cultures, and treatment of cells with an NO scavenger significantly reduced viral production. Pre-treatment of exponentially growing E. huxleyi cultures with the NO donor S-nitroso-N-acetylpenicillamine (SNAP) prior to challenge with hydrogen peroxide (H₂O₂) led to greater cell survival, suggesting that NO may have a cellular antioxidant function. Indeed, cell lysates generated from cultures treated with SNAP and undergoing infection displayed enhanced ability to detoxify H₂O₂. Lastly, we show that fluorescent indicators of cellular ROS, NO, and death, in combination with classic DNA- and lipid-based biomarkers of infection, can function as real-time diagnostic tools to identify and contextualize viral infection in natural E. huxleyi blooms.

Algal Oxylipins Mediate the Resistance of Diatoms against Algicidal Bacteria

Algicidal bacteria can lyse microalgal blooms and trigger shifts within plankton communities. Resistant algal species can escape lysis, and have the opportunity to dominate the phytoplankton after a bacterial infection. Despite their important function in ecosystem regulation, little is known about mechanisms of resistance. Here, we show that the diatom Chaetoceros didymus releases eicosanoid oxylipins into the medium, and that the lytic algicidal bacterium, Kordia algicida, induces the production of several wound-activated oxylipins in this resistant diatom. Neither releases nor an induction occurs in the susceptible diatom Skeletonema costatum that is lysed by the bacterium within a few days. Among the upregulated oxylipins, hydroxylated eicosapentaenoic acids (HEPEs) dominate. However, also, resolvins, known lipid mediators in mammals, increase upon exposure of the algae to the algicidal bacteria. The prevailing hydroxylated fatty acid, 15-HEPE, significantly inhibits growth of K. algicida at a concentration of approximately 1 µM. The oxylipin production may represent an independent line of defense of the resistant alga, acting in addition to the previously reported upregulation of proteases.