Draft genome sequence of the coccolithovirus EhV-84

Long-read genome sequencing provides novel insights into the harmful algal bloom species Prymnesium parvum

Prymnesium parvum is a toxin-producing haptophyte that causes harmful algal blooms worldwide, which are often associated with massive fish-kills and subsequent economic losses. In here, we present nuclear and plastid genome assemblies using PacBio HiFi long reads and DNBseq short reads for the two P. parvum strains UTEX 2797 and CCMP 3037, representing producers of type A prymnesins. Our results show that the P. parvum strains have a moderate haptophyte genome size of 97.56 and 107.32 Mb. The genome assemblies present one of highest contiguous assembled contig sequences to date consisting of 463 and 362 contigs with a contig N50 of 596.99 kb and 968.39 kb for strain UTEX 2797 and CCMP 3037, respectively. The assembled contigs of UTEX 2797 and CCMP 3037 were anchored to 34 scaffolds, with a scaffold N50 of 5.35 Mb and 3.61 Mb, respectively, accounting for 93.2 % and 97.9 % of the total length. Each plastid genome comprises a circular contig. A total of 20,578 and 19,426 protein-coding genes were annotated for UTEX 2797 and CCMP 3037. The expanded gene family analysis showed that starch and sucrose metabolism, sulfur metabolism, energy metabolism and ABC transporters are involved in the evolution of P. parvum. Polyketide synthase (PKS) genes responsible for the production of secondary metabolites such as prymnesins displayed different expression patterns under nutrient limitation. Overlap with repeats and horizontal gene transfer may be two contributing factors to the high number of PKS genes found in this species. The two high quality P. parvum genomes will serve as valuable resources for ecological, genetic, and toxicological studies of haptophytes that can be used to monitor and potentially manage harmful blooms of ichthyotoxic P. parvum in the future.

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.

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.

Emerging Interaction Patterns in the Emiliania huxleyi-EhV System

Viruses are thought to be fundamental in driving microbial diversity in the oceanic planktonic realm. That role and associated emerging infection patterns remain particularly elusive for eukaryotic phytoplankton and their viruses. Here we used a vast number of strains from the model system Emiliania huxleyi/Emiliania huxleyi Virus to quantify parameters such as growth rate (µ), resistance (R), and viral production (Vp) capacities. Algal and viral abundances were monitored by flow cytometry during 72-h incubation experiments. The results pointed out higher viral production capacity in generalist EhV strains, and the virus-host infection network showed a strong co-evolution pattern between E. huxleyi and EhV populations. The existence of a trade-off between resistance and growth capacities was not confirmed.

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

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

Metagenomic Investigation of Viral Communities in Ballast Water

Ballast water is one of the most important vectors for the transport of non-native species to new aquatic environments. Due to the development of new ballast water quality standards for viruses, this study aimed to determine the taxonomic diversity and composition of viral communities (viromes) in ballast and harbor waters using metagenomics approaches. Ballast waters from different sources within the North America Great Lakes and paired harbor waters were collected around the Port of Duluth-Superior. Bioinformatics analysis of over 550 million sequences showed that a majority of the viral sequences could not be assigned to any taxa associated with reference sequences, indicating the lack of knowledge on viruses in ballast and harbor waters. However, the assigned viruses were dominated by double-stranded DNA phages, and sequences associated with potentially emerging viral pathogens of fish and shrimp were detected with low amino acid similarity in both ballast and harbor waters. Annotation-independent comparisons showed that viromes were distinct among the Great Lakes, and the Great Lakes viromes were closely related to viromes of other cold natural freshwater systems but distant from viromes of marine and human designed/managed freshwater systems. These results represent the most detailed characterization to date of viruses in ballast water, demonstrating their diversity and the potential significance of the ship-mediated spread of viruses.

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

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

Complete genome sequence of Roseophage vB_DshP-R1, which infects Dinoroseobacter shibae DFL12

The Roseophages, a group of marine viruses that uniquely infect the Roseobacter clade of bacteria, play a significant role in marine ecosystems. Here we present a complete genomic sequence of an N4 phage ‘vB_DshP-R1’, which infects Dinoroseobacter shibae DFL12, together with its structural and genomic features. vB_DshP-R1 has an ~ 75 nm diameter icosahedral structure and a complete genome of 75,028 bp. This is the first genome sequence of a lytic phage of the genus Dinoroseobacter.

Complete genome sequence of Roseophage vB_DshP-R1, which infects Dinoroseobacter shibae DFL12

The Roseophages, a group of marine viruses that uniquely infect the Roseobacter clade of bacteria, play a significant role in marine ecosystems. Here we present a complete genomic sequence of an N4 phage ‘vB_DshP-R1’, which infects Dinoroseobacter shibae DFL12, together with its structural and genomic features. vB_DshP-R1 has an ~ 75 nm diameter icosahedral structure and a complete genome of 75,028 bp. This is the first genome sequence of a lytic phage of the genus Dinoroseobacter.

Dip in the gene pool: metagenomic survey of natural coccolithovirus communities

Despite the global oceanic distribution and recognised biogeochemical impact of coccolithoviruses (EhV), their diversity remains poorly understood. Here we employed a metagenomic approach to study the occurrence and progression of natural EhV community genomic variability. Analysis of EhV metagenomes from the early and late stages of an induced bloom led to three main discoveries. First, we observed resilient and specific genomic signatures in the EhV community associated with the Norwegian coast, which reinforce the existence of limitations to the capacity of dispersal and genomic exchange among EhV populations. Second, we identified a hyper-variable region (approximately 21kbp long) in the coccolithovirus genome. Third, we observed a clear trend for EhV relative amino-acid diversity to reduce from early to late stages of the bloom. This study validated two new methodological combinations, and proved very useful in the discovery of new genomic features associated with coccolithovirus natural communities.

Permanent draft genomes of four new coccolithoviruses: EhV-18, EhV-145, EhV-156 and EhV-164

Coccolithoviruses infect the marine coccolithophorid microalga Emiliania huxleyi. Here, we describe the genomes of four new coccolithoviruses isolated from UK coastal locations. Of particular interest, EhV-18 and EhV-145 encode serine palmitoyltransferase function via two distinct genes, whereas all other coccolithoviruses have SPT as a gene fusion of LCB1/LCB2 domains.

Functional inferences of environmental coccolithovirus biodiversity

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

Draft genome sequence of four coccolithoviruses: Emiliania huxleyi virus EhV-88, EhV-201, EhV-207, and EhV-208

The Coccolithoviridae are a group of viruses which infect the marine coccolithophorid microalga Emiliania huxleyi. The Emiliania huxleyi viruses (known as EhVs) described herein have 160- to 180-nm diameter icosahedral structures, have genomes of approximately 400 kbp, and consist of more than 450 predicted coding sequences (CDSs). Here, we describe the genomic features of four newly sequenced coccolithoviruses (EhV-88, EhV-201, EhV-207, and EhV-208) together with their draft genome sequences and their annotations, highlighting the homology and heterogeneity of these genomes to the EhV-86 model reference genome.

Draft genome sequence of the coccolithovirus Emiliania huxleyi virus 202

Emiliania huxleyi virus 202 (EhV-202) is a member of the Coccolithoviridae, a group of viruses that infect the marine coccolithophorid Emiliania huxleyi. EhV-202 has a 160- to 180-nm-diameter icosahedral structure and a genome of approximately 407 kbp, consisting of 485 coding sequences (CDSs). Here we describe the genomic features of EhV-202, together with a draft genome sequence and its annotation, highlighting the homology and heterogeneity of this genome in comparison with the EhV-86 reference genome.