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

Abstract

ABSTRACTHarmful algal blooms (HABs) of the toxic haptophytePrymnesium parvumare a recurrent problem in many inland and estuarine waters around the world. Strains ofP. parvumvary in the toxins they produce and in other physiological traits associated with HABs, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for fifteen phylogenetically and geographically diverse strains ofP. parvumincluding Hi-C guided, near-chromosome level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 Mbp to 845 Mbp. Strains included haploids, diploids, and polyploids, but not all differences in DNA content were due to variation in genome copy number. Haploid genome size between strains of different chemotypes differed by as much as 243 Mbp. Syntenic and phylogenetic analyses indicate that UTEX 2797, a common laboratory strain from Texas, is a hybrid that retains two phylogenetically distinct haplotypes. Investigation of gene families variably present across strains identified several functional categories associated with metabolism, including candidates for the biosynthesis of toxic metabolites, as well as genome size variation, including recent proliferations of transposable elements. Together, our results indicate thatP. parvumis comprised of multiple cryptic species. These genomes provide a robust phylogenetic and genomic framework for investigations into the eco-physiological consequences of the intra- and inter-specific genetic variation present inP. parvumand demonstrate the need for similar resources for other HAB-forming morphospecies.SIGNIFICANCE STATEMENTHarmful algal blooms (HABs) are a global concern. Efforts to understand the genetic basis of traits associated with the success of HAB-forming species are limited by a dearth of genomic resources. In this paper we present genomes for fifteen strains ofPrymnesium parvum, a toxic alga that causes ecosystem and societally disruptive HABs around the world. We uncover an unprecedented amount of sequence-level, gene family, and genome architecture evolution inP. parvumand provide evidence for both cryptic speciation and hybridization. These results illustrate how both inter- and intraspecific genetic variation can be dramatically underestimated in a protist morphospecies. More work is needed to understand the eco-physiological consequences of hidden genetic diversity inP. parvumand HAB-forming species more generally.