Chrysochromulina: Genomic assessment and taxonomic diagnosis of the type species for an oleaginous algal clade

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.

Phylogeny and genetic variations of the three genome compartments in haptophytes shed light on the rapid evolution of coccolithophores

Haptophyte algae, including coccolithophores, play key roles in global carbon cycling and ecosystem. They exhibit exceptional morphological and functional diversity. However, their phylogeny is mostly based on short markers and genome researches are always limited to few species, hindering a better understanding about their evolution and diversification. In this study, by assembling 69 new plastid genomes, 65 new mitochondrial genomes, and 55 nuclear drafts, we systematically analyzed their genome variations and built the most comprehensive phylogenies in haptophytes and Noelaerhabdaceae, with the latter is the family of the model coccolithophore Emiliania huxleyi. The haptophyte genomes vary significantly in size, gene content, and structure. We detected phylogenetic incongruence of Prymnesiales between genome compartments. In Noelaerhabdaceae, by including Reticulofenestra sessilis and a proper outgroup, we found R. sessilis was not the basal taxon of this family. Noelaerhabdaceae strains have very similar genomic features and conserved sequences, but different gene content and dynamic structure. We speculate that was caused by DNA double-strand break repairs. Our results provide valuable genetic resources and new insights into the evolution of haptophytes, especially coccolithophores.

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

Harmful algal blooms of the toxic haptophyte Prymnesium parvum are a recurrent problem in many inland and estuarine waters around the world. Strains of P. parvum vary in the toxins they produce and in other physiological traits associated with harmful algal blooms, but the genetic basis for this variation is unknown. To investigate genome diversity in this morphospecies, we generated genome assemblies for 15 phylogenetically and geographically diverse strains of P. parvum, including Hi-C guided, near-chromosome-level assemblies for two strains. Comparative analysis revealed considerable DNA content variation between strains, ranging from 115 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 the strains identified several functional categories associated with metabolic and genome size variation in P. parvum, including genes for the biosynthesis of toxic metabolites and proliferation of transposable elements. Together, our results indicate that P. parvum comprises 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 in P. parvum and demonstrate the need for similar resources for other harmful algal-bloom-forming morphospecies.

Isolation and infection cycle of a polinton-like virus virophage in an abundant marine alga

Virophages are small double stranded DNA (dsDNA) viruses that can only replicate in a host by co-infecting with another virus. Marine algae are commonly associated with virophage-like elements such as Polinton-like viruses (PLVs) that remain largely uncharacterized. Here we isolated a PLV that co-infects the alga Phaeocystis globosa with the Phaeocystis globosa virus-14T (PgV-14T), a close relative of the "Phaeocystis globosa virus-virophage" genomic sequence. We name this PLV 'Gezel-14T. Gezel is phylogenetically distinct from the Lavidaviridae family where all known virophages belong. Gezel-14T co-infection decreases the fitness of its viral host by reducing burst sizes of PgV-14T, yet insufficiently to spare the cellular host population. Genomic screens show Gezel-14T-like PLVs integrated into Phaeocystis genomes, suggesting that these widespread viruses are capable of integration into cellular host genomes. This system presents an opportunity to better understand the evolution of eukaryotic dsDNA viruses as well as the complex dynamics and implications of viral parasitism.

High-resolution metagenomic reconstruction of the freshwater spring bloom

BACKGROUND

The phytoplankton spring bloom in freshwater habitats is a complex, recurring, and dynamic ecological spectacle that unfolds at multiple biological scales. Although enormous taxonomic shifts in microbial assemblages during and after the bloom have been reported, genomic information on the microbial community of the spring bloom remains scarce.

RESULTS

We performed a high-resolution spatio-temporal sampling of the spring bloom in a freshwater reservoir and describe a multitude of previously unknown taxa using metagenome-assembled genomes of eukaryotes, prokaryotes, and viruses in combination with a broad array of methodologies. The recovered genomes reveal multiple distributional dynamics for several bacterial groups with progressively increasing stratification. Analyses of abundances of metagenome-assembled genomes in concert with CARD-FISH revealed remarkably similar in situ doubling time estimates for dominant genome-streamlined microbial lineages. Discordance between quantitations of cryptophytes arising from sequence data and microscopic identification suggested the presence of hidden, yet extremely abundant aplastidic cryptophytes that were confirmed by CARD-FISH analyses. Aplastidic cryptophytes are prevalent throughout the water column but have never been considered in prior models of plankton dynamics. We also recovered the first metagenomic-assembled genomes of freshwater protists (a diatom and a haptophyte) along with thousands of giant viral genomic contigs, some of which appeared similar to viruses infecting haptophytes but owing to lack of known representatives, most remained without any indication of their hosts. The contrasting distribution of giant viruses that are present in the entire water column to that of parasitic perkinsids residing largely in deeper waters allows us to propose giant viruses as the biological agents of top-down control and bloom collapse, likely in combination with bottom-up factors like a nutrient limitation.

CONCLUSION

We reconstructed thousands of genomes of microbes and viruses from a freshwater spring bloom and show that such large-scale genome recovery allows tracking of planktonic succession in great detail. However, integration of metagenomic information with other methodologies (e.g., microscopy, CARD-FISH) remains critical to reveal diverse phenomena (e.g., distributional patterns, in situ doubling times) and novel participants (e.g., aplastidic cryptophytes) and to further refine existing ecological models (e.g., factors affecting bloom collapse). This work provides a genomic foundation for future approaches towards a fine-scale characterization of the organisms in relation to the rapidly changing environment during the course of the freshwater spring bloom. Video Abstract.

The complete mitochondrial genome of Isochrysis galbana harbors a unique repeat structure and a specific trans-spliced cox1 gene

The haptophyte Isochrysis galbana is considered as a promising source for food supplements due to its rich fucoxanthin and polyunsaturated fatty acids content. Here, the I. galbana mitochondrial genome (mitogenome) was sequenced using a combination of Illumina and PacBio sequencing platforms. This 39,258 bp circular mitogenome has a total of 46 genes, including 20 protein-coding genes, 24 tRNA genes and two rRNA genes. A large block of repeats (~12.7 kb) was segregated in one region of the mitogenome, accounting for almost one third of the total size. A trans-spliced gene cox1 was first identified in I. galbana mitogenome and was verified by RNA-seq and DNA-seq data. The massive expansion of tandem repeat size and cis- to trans-splicing shift could be explained by the high mitogenome rearrangement rates in haptophytes. Strict SNP calling based on deep transcriptome sequencing data suggested the lack of RNA editing in both organelles in this species, consistent with previous studies in other algal lineages. To gain insight into haptophyte mitogenome evolution, a comparative analysis of mitogenomes within haptophytes and among eight main algal lineages was performed. A core gene set of 15 energy and metabolism genes is present in haptophyte mitogenomes, consisting of 1 cob, 3 cox, 7 nad, 2 atp and 2 ribosomal genes. Gene content and order was poorly conserved in this lineage. Haptophyte mitogenomes have lost many functional genes found in many other eukaryotes including rps/rpl, sdh, tat, secY genes, which make it contain the smallest gene set among all algal taxa. All these implied the rapid-evolving and more recently evolved mitogenomes of haptophytes compared to other algal lineages. The phylogenetic tree constructed by cox1 genes of 204 algal mitogenomes yielded well-resolved internal relationships, providing new evidence for red-lineages that contained plastids of red algal secondary endosymbiotic origin. This newly assembled mitogenome will add to our knowledge of general trends in algal mitogenome evolution within haptophytes and among different algal taxa.

Phylogenetic analyses and reclassification of the oleaginous marine species Nannochloris sp. "desiccata" (Trebouxiophyceae, Chlorophyta), formerly Chlorella desiccata, supported by a high-quality genome assembly

Microalgae are diverse, with many gaps remaining in phylogenetic and physiological understanding. Thus, studying new microalgae species increases our broader comprehension of biological diversity, and evaluation of new candidates as algal production platforms can lead to improved productivity under a variety of cultivation conditions. Chlorella is a genus of fast-growing species often isolated from freshwater habitats and cultivated as a source of nutritional supplements. However, the use of freshwater increases competition with other freshwater needs. We identified Chlorella desiccata to be worthy of further investigation as a potential algae production strain, due to its isolation from a marine environment and its promising growth and biochemical composition properties. Long-read genomic sequencing was conducted for C. desiccata UTEX 2526, resulting in a high-quality, near chromosome level, diploid genome with an assembly length of 21.55 Mbp in only 18 contigs. We also report complete circular mitochondrial and chloroplast genomes. Phylogenomic and phylogenetic analyses using nuclear, chloroplast, 18S rRNA, and actin sequences revealed that this species clades within strains currently identified as Nannochloris (Trebouxiophyceae, Chlorophyta), leading to its reclassification as Nannochloris sp. "desiccata" UTEX 2526. The mode of cell division for this species is autosporulation, differing from the type species N. bacillaris. As has occurred across multiple microalgae genera, there are repeated examples of Nannochloris species reclassification in the literature. This high-quality genome assembly and phylogenetic analysis of the potential algal production strain Nannochloris sp. "desiccata" UTEX 2526 provides an important reference and useful tool for further studying this region of the phylogenetic tree.

The Genome of the Haptophyte Diacronema lutheri (Pavlova lutheri, Pavlovales): A Model for Lipid Biosynthesis in Eukaryotic Algae

Haptophytes are biogeochemically and industrially important protists with underexplored genomic diversity. We present a nuclear genome assembly for the class Pavlovales, which was assembled with PacBio long-read data into highly contiguous sequences. We sequenced strain Diacronema lutheri NIVA-4/92, formerly known as Pavlova lutheri, because it has established roles in aquaculture and has been a key organism for studying microalgal lipid biosynthesis. Our data show that D. lutheri has the smallest and most streamlined haptophycean genome assembled to date, with an assembly size of 43.503 Mb and 14,446 protein-coding genes. Together with its high nuclear GC content, Diacronema is an important genus for investigating selective pressures on haptophyte genome evolution, contrasting with the much larger and more repetitive genome of the coccolithophore Emiliania huxleyi. The D. lutheri genome will be a valuable resource for resolving the genetic basis of algal lipid biosynthesis and metabolic remodeling that takes place during adaptation and stress response in natural and engineered environments.

Large Differences in the Haptophyte Phaeocystis globosa Mitochondrial Genomes Driven by Repeat Amplifications

The haptophyte Phaeocystis globosa is a well-known species for its pivotal role in global carbon and sulfur cycles and for its capability of forming harmful algal blooms (HABs) with serious ecological consequences. Its mitochondrial genome (mtDNA) sequence has been reported in 2014 but it remains incomplete due to its long repeat sequences. In this study, we constructed the first full-length mtDNA of P. globosa, which was a circular genome with a size of 43,585 bp by applying the PacBio single molecular sequencing method. The mtDNA of this P. globosa strain (CNS00066), which was isolated from the Beibu Gulf, China, encoded 19 protein-coding genes (PCGs), 25 tRNA genes, and two rRNA genes. It contained two large repeat regions of 6.7 kb and ∼14.0 kb in length, respectively. The combined length of these two repeat regions, which were missing from the previous mtDNA assembly, accounted for almost half of the entire mtDNA and represented the longest repeat region among all sequenced haptophyte mtDNAs. In this study, we tested the hypothesis that repeat unit amplification is a driving force for different mtDNA sizes. Comparative analysis of mtDNAs of five additional P. globosa strains (four strains obtained in this study, and one strain previously published) revealed that all six mtDNAs shared identical numbers of genes but with dramatically different repeat regions. A homologous repeat unit was identified but with hugely different numbers of copies in all P. globosa strains. Thus, repeat amplification may represent an important driving force of mtDNA evolution in P. globosa.

Pyrenoids: CO2-fixing phase separated liquid organelles

Pyrenoids are non-membrane bound organelles found in chloroplasts of algae and hornwort plants that can be seen by light-microscopy. Pyrenoids are formed by liquid-liquid phase separation (LLPS) of Rubisco, the primary CO₂ fixing enzyme, with an intrinsically disordered multivalent Rubisco-binding protein. Pyrenoids are the heart of algal and hornwort biophysical CO₂ concentrating mechanisms, which accelerate photosynthesis and mediate about 30% of global carbon fixation. Even though LLPS may underlie the apparent convergent evolution of pyrenoids, our current molecular understanding of pyrenoid formation comes from a single example, the model alga Chlamydomonas reinhardtii. In this review, we summarise current knowledge about pyrenoid assembly, regulation and structural organization in Chlamydomonas and highlight evidence that LLPS is the general principle underlying pyrenoid formation across algal lineages and hornworts. Detailed understanding of the principles behind pyrenoid assembly, regulation and structural organization within diverse lineages will provide a fundamental understanding of this biogeochemically important organelle and help guide ongoing efforts to engineer pyrenoids into crops to increase photosynthetic performance and yields.².