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

Single-cell RNA-seq of the rare virosphere reveals the native hosts of giant viruses in the marine environment

Giant viruses (phylum Nucleocytoviricota) are globally distributed in aquatic ecosystems. They play fundamental roles as evolutionary drivers of eukaryotic plankton and regulators of global biogeochemical cycles. However, we lack knowledge about their native hosts, hindering our understanding of their life cycle and ecological importance. In the present study, we applied a single-cell RNA sequencing (scRNA-seq) approach to samples collected during an induced algal bloom, which enabled pairing active giant viruses with their native protist hosts. We detected hundreds of single cells from multiple host lineages infected by diverse giant viruses. These host cells included members of the algal groups Chrysophycae and Prymnesiophycae, as well as heterotrophic flagellates in the class Katablepharidaceae. Katablepharids were infected with a rare Imitervirales-07 giant virus lineage expressing a large repertoire of cell-fate regulation genes. Analysis of the temporal dynamics of these host-virus interactions revealed an important role for the Imitervirales-07 in controlling the population size of the host Katablepharid population. Our results demonstrate that scRNA-seq can be used to identify previously undescribed host-virus interactions and study their ecological importance and impact.

Homing in on the rare virosphere reveals the native host of giant viruses

Giant viruses (phylum Nucleocytoviricota) are globally distributed in aquatic ecosystems1,2. They play major roles as evolutionary drivers of eukaryotic plankton3 and regulators of global biogeochemical cycles4. Recent metagenomic studies have significantly expanded the known diversity of marine giant viruses1,5-7, but we still lack fundamental knowledge about their native hosts, thereby hindering our understanding of their lifecycle and ecological importance. Here, we aim to discover the native hosts of giant viruses using a novel, sensitive single-cell metatranscriptomic approach. By applying this approach to natural plankton communities, we unraveled an active viral infection of several giant viruses, from multiple lineages, and identified their native hosts. We identify a rare lineage of giant virus (Imitervirales-07) infecting a minute population of protists (class Katablepharidaceae) and revealed the prevalence of highly expressed viral-encoded cell-fate regulation genes in infected cells. Further examination of this host-virus dynamics in a temporal resolution suggested this giant virus controls its host population demise. Our results demonstrate how single-cell metatranscriptomics is a sensitive approach for pairing viruses with their authentic hosts and studying their ecological significance in a culture-independent manner in the marine environment.

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

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

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.

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.

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.

Closely related viruses of the marine picoeukaryotic alga Ostreococcus lucimarinus exhibit different ecological strategies

In marine ecosystems, viruses are major disrupters of the direct flow of carbon and nutrients to higher trophic levels. Although the genetic diversity of several eukaryotic phytoplankton virus groups has been characterized, their infection dynamics are less understood, such that the physiological and ecological implications of their diversity remain unclear. We compared genomes and infection phenotypes of the two most closely related cultured phycodnaviruses infecting the widespread picoprasinophyte Ostreococcus lucimarinus under standard- (1.3 divisions per day) and limited-light (0.41 divisions per day) nutrient replete conditions. OlV7 infection caused early arrest of the host cell cycle, coinciding with a significantly higher proportion of infected cells than OlV1-amended treatments, regardless of host growth rate. OlV7 treatments showed a near-50-fold increase of progeny virions at the higher host growth rate, contrasting with OlV1's 16-fold increase. However, production of OlV7 virions was more sensitive than OlV1 production to reduced host growth rate, suggesting fitness trade-offs between infection efficiency and resilience to host physiology. Moreover, although organic matter released from OlV1- and OlV7-infected hosts had broadly similar chemical composition, some distinct molecular signatures were observed. Collectively, these results suggest that current views on viral relatedness through marker and core gene analyses underplay operational divergence and consequences for host ecology.

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.

Ancestrality and Mosaicism of Giant Viruses Supporting the Definition of the Fourth TRUC of Microbes

Giant viruses of amoebae were discovered in 2003. Since then, their diversity has greatly expanded. They were suggested to form a fourth branch of life, collectively named ‘TRUC’ (for “Things Resisting Uncompleted Classifications”) alongside Bacteria, Archaea, and Eukarya. Their origin and ancestrality remain controversial. Here, we specify the evolution and definition of giant viruses. Phylogenetic and phenetic analyses of informational gene repertoires of giant viruses and selected bacteria, archaea and eukaryota were performed, including structural phylogenomics based on protein structural domains grouped into 289 universal fold superfamilies (FSFs). Hierarchical clustering analysis was performed based on a binary presence/absence matrix constructed using 727 informational COGs from cellular organisms. The presence/absence of ‘universal’ FSF domains was used to generate an unrooted maximum parsimony phylogenomic tree. Comparison of the gene content of a giant virus with those of a bacterium, an archaeon, and a eukaryote with small genomes was also performed. Overall, both cladistic analyses based on gene sequences of very central and ancient proteins and on highly conserved protein fold structures as well as phenetic analyses were congruent regarding the delineation of a fourth branch of microbes comprised by giant viruses. Giant viruses appeared as a basal group in the tree of all proteomes. A pangenome and core genome determined for Rickettsia bellii (bacteria), Methanomassiliicoccus luminyensis (archaeon), Encephalitozoon intestinalis (eukaryote), and Tupanvirus (giant virus) showed a substantial proportion of Tupanvirus genes that overlap with those of the cellular microbes. In addition, a substantial genome mosaicism was observed, with 51, 11, 8, and 0.2% of Tupanvirus genes best matching with viruses, eukaryota, bacteria, and archaea, respectively. Finally, we found that genes themselves may be subject to lateral sequence transfers. In summary, our data highlight the quantum leap between classical and giant viruses. Phylogenetic and phyletic analyses and the study of protein fold superfamilies confirm previous evidence of the existence of a fourth TRUC of life that includes giant viruses, and highlight its ancestrality and mosaicism. They also point out that best evolutionary representations for giant viruses and cellular microorganisms are rhizomes, and that sequence transfers rather than gene transfers have to be considered.

Viruses of Eukaryotic Algae: Diversity, Methods for Detection, and Future Directions

The scope for ecological studies of eukaryotic algal viruses has greatly improved with the development of molecular and bioinformatic approaches that do not require algal cultures. Here, we review the history and perceived future opportunities for research on eukaryotic algal viruses. We begin with a summary of the 65 eukaryotic algal viruses that are presently in culture collections, with emphasis on shared evolutionary traits (e.g., conserved core genes) of each known viral type. We then describe how core genes have been used to enable molecular detection of viruses in the environment, ranging from PCR-based amplification to community scale "-omics" approaches. Special attention is given to recent studies that have employed network-analyses of -omics data to predict virus-host relationships, from which a general bioinformatics pipeline is described for this type of approach. Finally, we conclude with acknowledgement of how the field of aquatic virology is adapting to these advances, and highlight the need to properly characterize new virus-host systems that may be isolated using preliminary molecular surveys. Researchers can approach this work using lessons learned from the Chlorella virus system, which is not only the best characterized algal-virus system, but is also responsible for much of the foundation in the field of aquatic virology.

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.

Prasinovirus Attack of Ostreococcus Is Furtive by Day but Savage by Night

Prasinoviruses are large DNA viruses that infect diverse genera of green microalgae worldwide in aquatic ecosystems, but molecular knowledge of their life cycles is lacking. Several complete genomes of both these viruses and their marine algal hosts are now available and have been used to show the pervasive presence of these species in microbial metagenomes. We have analyzed the life cycle of Ostreococcus tauri virus 5 (OtV5), a lytic virus, using transcriptome sequencing (RNA-Seq) from 12 time points of healthy or infected Ostreococcus tauri cells over a day/night cycle in culture. In the day, viral gene transcription remained low while host nitrogen metabolism gene transcription was initially strongly repressed for two successive time points before being induced for 8 h, but during the night, viral transcription increased steeply while host nitrogen metabolism genes were repressed and many host functions that are normally reduced in the dark appeared to be compensated either by genes expressed from the virus or by increased expression of a subset of 4.4% of the host's genes. Some host cells underwent lysis progressively during the night, but a larger proportion were lysed the following morning. Our data suggest that the life cycles of algal viruses mirror the diurnal rhythms of their hosts.IMPORTANCE Prasinoviruses are common in marine environments, and although several complete genomes of these viruses and their hosts have been characterized, little is known about their life cycles. Here we analyze in detail the transcriptional changes occurring over a 27-h-long experiment in a natural diurnal rhythm, in which the growth of host cells is to some extent synchronized, so that host DNA replication occurs late in the day or early in the night and cell division occurs during the night. Surprisingly, viral transcription remains quiescent over the daytime, when the most energy (from light) is available, but during the night viral transcription activates, accompanied by expression of a few host genes that are probably required by the virus. Although our experiment was accomplished in the lab, cyclical changes have been documented in host transcription in the ocean. Our observations may thus be relevant for eukaryotic phytoplankton in natural environments.

Cyanophage-encoded lipid desaturases: oceanic distribution, diversity and function

Cyanobacteria are among the most abundant photosynthetic organisms in the oceans; viruses infecting cyanobacteria (cyanophages) can alter cyanobacterial populations, and therefore affect the local food web and global biochemical cycles. These phages carry auxiliary metabolic genes (AMGs), which rewire various metabolic pathways in the infected host cell, resulting in increased phage fitness. Coping with stress resulting from photodamage appears to be a central necessity of cyanophages, yet the overall mechanism is poorly understood. Here we report a novel, widespread cyanophage AMG, encoding a fatty acid desaturase (FAD), found in two genotypes with distinct geographical distribution. FADs are capable of modulating the fluidity of the host’s membrane, a fundamental stress response in living cells. We show that both viral FAD (vFAD) families are Δ9 lipid desaturases, catalyzing the desaturation at carbon 9 in C16 fatty acid chains. In addition, we present a comprehensive fatty acid profiling for marine cyanobacteria, which suggests a unique desaturation pathway of medium- to long-chain fatty acids no longer than C16, in accordance with the vFAD activity. Our findings suggest that cyanophages are capable of fiddling with the infected host’s membranes, possibly leading to increased photoprotection and potentially enhancing viral-encoded photosynthetic proteins, resulting in a new viral metabolic network.