Extreme genome diversity and cryptic speciation in a harmful algal-bloom-forming eukaryote

Snow alga Sanguina aurantia as revealed through de novo genome assembly and annotation

To thrive on melting alpine and polar snow, some Chlorophytes produce an abundance of astaxanthin, causing red blooms, often dominated by genus Sanguina. The red cells have not been cultured, but we recently grew a green biciliate conspecific with Sanguina aurantia from a sample of watermelon snow. This culture provided source material for Oxford Nanopore Technology and Illumina sequencing. Our assembly pipeline exemplifies the value of a hybrid long- and short-read approach for the complexities of working with a culture grown from a field sample. Using bioinformatic tools, we separated assembled contigs into 2 genomic pools based on a difference in GC content (57.5 and 55.1%). We present the data as 2 assemblies of S. aurantia variants but explore other possibilities. High-throughput chromatin conformation capture analysis (Hi-C sequencing) was used to scaffold the assemblies into a 96-Mb genome designated as "A" and a 102-Mb genome designated as "B." Both assemblies are highly contiguous: genome A consists of 38 scaffolds with an N50 of 5.4 Mb, while genome B has 50 scaffolds with an N50 of 6.4 Mb. RNA sequencing was used to improve gene annotation.

Hormetic and transcriptomic responses of the toxic alga Prymnesium parvum to glyphosate

Growth of the toxic alga Prymnesium parvum is hormetically stimulated with environmentally relevant concentrations of glyphosate. The mechanisms of glyphosate hormesis in this species, however, are unknown. We evaluated the transcriptomic response of P. parvum to glyphosate at concentrations that stimulate maximum growth and where growth is not different from control values, the zero-equivalent point (ZEP). Maximum growth occurred at 0.1 mg l-1 and the ZEP was 2 mg l-1. At 0.1 mg l-1, upregulated transcripts outnumbered downregulated transcripts by one order of magnitude. Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that the upregulated transcriptome is primarily associated with metabolism and biosynthesis. Transcripts encoding heat shock proteins and co-chaperones were among the most strongly upregulated, and several others were associated with translation, Redox homeostasis, cell replication, and photosynthesis. Although most of the same transcripts were also upregulated at concentrations ≥ZEP, the proportion of downregulated transcripts greatly increased as glyphosate concentrations increased. At the ZEP, downregulated transcripts were associated with photosynthesis, cell replication, and anion transport, indicating that specific interference with these processes is responsible for the nullification of hormetic growth. Transcripts encoding the herbicidal target of glyphosate, 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), were upregulated at concentrations ≥ZEP but not at 0.1 mg l-1, indicating that disruption of EPSPS activity occurred at high concentrations and that nullification of hormetic growth involves the direct interaction of glyphosate with this enzyme. Results of this study may contribute to a better understanding of glyphosate hormesis and of anthropogenic factors that influence P. parvum biogeography and bloom formation.

The haplotype-resolved Prymnesium parvum (type B) microalga genome reveals the genetic basis of its fish-killing toxins

The catastrophic loss of aquatic life in the Central European Oder River in 2022, caused by a toxic bloom of the haptophyte microalga Prymnesium parvum (in a wide sense, s.l.), underscores the need to improve our understanding of the genomic basis of the toxin. Previous morphological, phylogenetic, and genomic studies have revealed cryptic diversity within P. parvum s.l. and uncovered three clade-specific (types A, B, and C) prymnesin toxins. Here, we used state-of-the-art long-read sequencing and assembled the first haplotype-resolved diploid genome of a P. parvum type B from the strain responsible for the Oder disaster. Comparative analyses with type A genomes uncovered a genome-size expansion driven by repetitive elements in type B. We also found conserved synteny but divergent evolution in several polyketide synthase (PKS) genes, which are known to underlie toxin production in combination with environmental cues. We identified an approximately 20-kbp deletion in the largest PKS gene of type B that we link to differences in the chemical structure of types A and B prymnesins. Flow cytometry and electron microscopy analyses confirmed diploidy in the Oder River strain and revealed differences to closely related strains in both ploidy and morphology. Our results provide unprecedented resolution of strain diversity in P. parvum s.l. and a better understanding of the genomic basis of toxin variability in haptophytes. The reference-quality genome will enable us to better understand changes in microbial diversity in the face of increasing environmental pressures and provides a basis for strain-level monitoring of invasive Prymnesium in the future.

Giant polyketide synthase enzymes in the biosynthesis of giant marine polyether toxins

Prymnesium parvum are harmful haptophyte algae that cause massive environmental fish kills. Their polyketide polyether toxins, the prymnesins, are among the largest nonpolymeric compounds in nature and have biosynthetic origins that have remained enigmatic for more than 40 years. In this work, we report the "PKZILLAs," massive P. parvum polyketide synthase (PKS) genes that have evaded previous detection. PKZILLA-1 and -2 encode giant protein products of 4.7 and 3.2 megadaltons that have 140 and 99 enzyme domains. Their predicted polyene product matches the proposed pre-prymnesin precursor of the 90-carbon-backbone A-type prymnesins. We further characterize the variant PKZILLA-B1, which is responsible for the shorter B-type analog prymnesin-B1, from P. parvum RCC3426 and thus establish a general model of haptophyte polyether biosynthetic logic. This work expands expectations of genetic and enzymatic size limits in biology.

From genes to toxins: Profiling Prymnesium parvum during a riverine harmful algal bloom

Blooms of Prymnesium parvum, a unicellular alga globally distributed in marine and brackish environments, frequently result in massive fish kills due to the production of toxins called prymnesins by this haptophyte. In August 2022, a harmful algal bloom (HAB) of this species occurred in the lower Oder River (Poland and Germany), which caused mass mortalities of fish and other organisms. This HAB was linked to low discharge of the Oder and mining activities that caused a significant increase in salinity. In this context, we report on the molecular detection and screening of this haptophyte and its toxins in environmental samples and clonal cultures derived thereof. Both conventional PCR and droplet digital PCR assays reliably detected P. parvum in environmental samples. eDNA metabarcoding using the V4 region of the 18S rRNA gene revealed a single Prymnesium sequence variant, but failed to identify it to species level. Four clonal cultures established from environmental samples were unambiguously identified as P. parvum by molecular phylogenetics (near full-length 18S rRNA gene) and light microscopy. Phylogenetic analysis (ITS1-5.8S-ITS2 marker region) placed the cultured phylotype within a clade containing other P. parvum strains known to produce B-type prymnesins. Toxin-screening of the cultures using liquid chromatography-electrospray ionization - time of flight mass spectrometry identified B-type prymnesins, which were also detected in extracts of filter residues from water samples of the Oder collected during the HAB. Overall, our investigation provides a detailed characterization of P. parvum, including their prymnesins, during this HAB in the Oder River, contributing valuable insights into this ecological disaster. In addition, the droplet digital PCR assay established here will be useful for future monitoring of low levels of P. parvum on the Oder River or any other salt-impacted and brackish water bodies.

Giant polyketide synthase enzymes biosynthesize a giant marine polyether biotoxin

Prymnesium parvum are harmful haptophyte algae that cause massive environmental fish-kills. Their polyketide polyether toxins, the prymnesins, are amongst the largest nonpolymeric compounds in nature, alongside structurally-related health-impacting "red-tide" polyether toxins whose biosynthetic origins have been an enigma for over 40 years. Here we report the 'PKZILLAs', massive P. parvum polyketide synthase (PKS) genes, whose existence and challenging genomic structure evaded prior detection. PKZILLA-1 and -2 encode giant protein products of 4.7 and 3.2 MDa with 140 and 99 enzyme domains, exceeding the largest known protein titin and all other known PKS systems. Their predicted polyene product matches the proposed pre-prymnesin precursor of the 90-carbon-backbone A-type prymnesins. This discovery establishes a model system for microalgal polyether biosynthesis and expands expectations of genetic and enzymatic size limits in biology.

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.

Review of Harmful Algal Blooms (HABs) Causing Marine Fish Kills: Toxicity and Mitigation

Extensive growth of microscopic algae and cyanobacteria results in harmful algal blooms (HABs) in marine, brackish, and freshwater environments. HABs can harm humans and animals through their toxicity or by producing ecological conditions such as oxygen depletion, which can kill fish and other economically or ecologically important organisms. This review summarizes the reports on various HABs that are able to bring about marine fish kills. The predominant HABs, their toxins, and their effects on fishes spread across various parts of the globe are discussed. The mechanism of HAB-driven fish kills is discussed based on the available reports, and existing mitigation methods are presented. Lapses in the large-scale implementation of mitigation methods demonstrated under laboratory conditions are projected. Clay-related technologies and nano-sorption-based nanotechnologies, although proven to make significant contributions, have not been put to use in real-world conditions. The gaps in the technology transfer of the accomplished mitigation prototypes are highlighted. Further uses of remote sensing and machine learning state-of-the-art techniques for the detection and identification of HABs are recommended.

Towards effective management of the marine-origin Prymnesium parvum (Haptophyta): A growing concern in freshwater reservoirs?

Freshwater ecosystems are highly susceptible to harmful algal blooms (HABs), which are often caused by monospecific dense blooms. Effective preventive management strategies are urgently needed to avoid wide-ranging and severe impacts often resulting in costly damage to resources and unsustainable management options. In this study, we utilized SDM techniques focused on Prymnesium parvum, one of the most notorious HABs species worldwide. We first compare the climatic space occupied by P. parvum in North America, Europe and Australia. Additionally, we use MaxEnt algorithm to infer, for the first time, the potentially suitable freshwater environments in the aforementioned ranges. We also discuss the risks of invasion in reservoirs - prone habitats to persistent blooms of pests and invasive phytoplanktonic species. Our results show populations with distinctive niches suggesting ecophysiological tolerances, perhaps reflecting different strains. Our model projections revealed that the potential extent for P. parvum invasions is much broader than its current geographic distribution. The spatial configuration of reservoirs, if not sustaining dense blooms due to non-optimal conditions, favors colonization of multiple basins and ecoregions not yet occupied by P. parvum. Our models can provide valuable insights to decision-makers and monitoring programs while reducing the resources required to control the spread of P. parvum in disturbed habitats. Lastly, as impact magnitude is influenced by toxicity which in turn varies between different strains, we suggest future studies to incorporate intraspecific genetic information and fine-scale environmental variables to estimate potential distribution of P. parvum.