Evolutionary processes and cellular functions underlying divergence in Alexandrium minutum

Rapid diversification underlying the global dominance of a cosmopolitan phytoplankton

Marine phytoplankton play important roles in the global ecosystem, with a limited number of cosmopolitan keystone species driving their biomass. Recent studies have revealed that many of these phytoplankton are complexes composed of sibling species, but little is known about the evolutionary processes underlying their formation. Gephyrocapsa huxleyi, a widely distributed and abundant unicellular marine planktonic algae, produces calcified scales (coccoliths), thereby significantly affects global biogeochemical cycles via sequestration of inorganic carbon. This species is composed of morphotypes defined by differing degrees of coccolith calcification, the evolutionary ecology of which remains unclear. Here, we report an integrated morphological, ecological and genomic survey across globally distributed G. huxleyi strains to reconstruct evolutionary relationships between morphotypes in relation to their habitats. While G. huxleyi has been considered a single cosmopolitan species, our analyses demonstrate that it has evolved to comprise at least three distinct species, which led us to formally revise the taxonomy of the G. huxleyi complex. Moreover, the first speciation event occurred before the onset of the last interglacial period (~140 ka), while the second followed during this interglacial. Then, further rapid diversifications occurred during the most recent ice-sheet expansion of the last glacial period and established morphotypes as dominant populations across environmental clines. These results suggest that glacial-cycle dynamics contributed to the isolation of ocean basins and the segregations of oceans fronts as extrinsic drivers of micro-evolutionary radiations in extant marine phytoplankton.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Genetic association of toxin production in the dinoflagellate Alexandrium minutum

Dinoflagellates of the genus Alexandrium are responsible for harmful algal blooms and produce paralytic shellfish toxins (PSTs). Their very large and complex genomes make it challenging to identify the genes responsible for toxin synthesis. A family-based genomic association study was developed to determine the inheritance of toxin production in Alexandrium minutum and identify genomic regions linked to this production. We show that the ability to produce toxins is inheritable in a Mendelian way, while the heritability of the toxin profile is more complex. We developed the first dinoflagellate genetic linkage map. Using this map, several major results were obtained: 1. A genomic region related to the ability to produce toxins was identified. 2. This region does not contain any polymorphic sxt genes, known to be involved in toxin production in cyanobacteria. 3. The sxt genes, known to be present in a single cluster in cyanobacteria, are scattered on different linkage groups in A. minutum. 4. The expression of two sxt genes not assigned to any linkage group, sxtI and sxtG, may be regulated by the genomic region related to the ability to produce toxins. Our results provide new insights into the organization of toxicity-related genes in A. minutum, suggesting a dissociated genetic mechanism for the production of the different analogues and the ability to produce toxins. However, most of the newly identified genes remain unannotated. This study therefore proposes new candidate genes to be further explored to understand how dinoflagellates synthesize their toxins.

Insights into Alexandrium minutum Nutrient Acquisition, Metabolism and Saxitoxin Biosynthesis through Comprehensive Transcriptome Survey

Simple Summary

Alexandrium minutum is one of the causing organisms for the occurrence of harmful algae bloom (HABs) in marine ecosystems. This species produces saxitoxin, one of the deadliest neurotoxins which can cause human mortality. However, molecular information such as genes and proteins catalog on this species is still lacking. Therefore, this study has successfully characterized several new molecular mechanisms regarding A. minutum environmental adaptation and saxitoxin biosynthesis. Ultimately, this study provides a valuable resource for facilitating future dinoflagellates’ molecular response to environmental changes.

Abstract

The toxin-producing dinoflagellate Alexandrium minutum is responsible for the outbreaks of harmful algae bloom (HABs). It is a widely distributed species and is responsible for producing paralytic shellfish poisoning toxins. However, the information associated with the environmental adaptation pathway and toxin biosynthesis in this species is still lacking. Therefore, this study focuses on the functional characterization of A. minutum unigenes obtained from transcriptome sequencing using the Illumina Hiseq 4000 sequencing platform. A total of 58,802 (47.05%) unigenes were successfully annotated using public databases such as NCBI-Nr, UniprotKB, EggNOG, KEGG, InterPRO and Gene Ontology (GO). This study has successfully identified key features that enable A. minutum to adapt to the marine environment, including several carbon metabolic pathways, assimilation of various sources of nitrogen and phosphorus. A. minutum was found to encode homologues for several proteins involved in saxitoxin biosynthesis, including the first three proteins in the pathway of saxitoxin biosynthesis, namely sxtA, sxtG and sxtB. The comprehensive transcriptome analysis presented in this study represents a valuable resource for understanding the dinoflagellates molecular metabolic model regarding nutrient acquisition and biosynthesis of saxitoxin.

Characterization of toxin-producing strains of Dinophysis spp. (Dinophyceae) isolated from French coastal waters, with a particular focus on the D. acuminata-complex

Dinoflagellates of the genus Dinophysis are the most prominent producers of Diarrhetic Shellfish Poisoning (DSP) toxins which have an impact on public health and on marine aquaculture worldwide. In particular, Dinophysis acuminata has been reported as the major DSP agent in Western Europe. Still, its contribution to DSP events in the regions of the English Channel and the Atlantic coast of France, and the role of the others species of the Dinophysis community in these areas are not as clear. In addition, species identification within the D. acuminata complex has proven difficult due to their highly similar morphological features. In the present study, 30 clonal strains of the dominant Dinophysis species have been isolated from French coasts including the English Channel (3 sites), the Atlantic Ocean (11 sites) and the Mediterranean Sea (6 sites). Morphologically, strains were identified as three species: D. acuta, D. caudata, D. tripos, as well as the D. acuminata-complex. Sequences of the ITS and LSU rDNA regions confirmed these identifications and revealed no genetic difference within the D. acuminata-complex. Using the mitochondrial gene cox1, two groups of strains differing by only one substitution were found in the D. acuminata-complex, but SEM analysis of various strains showed a large range of morphological variations. Based on geographical origin and morphology, strains of the subclade A were ascribed to 'D. acuminata' while those of the subclade B were ascribed to 'D. sacculus'. Nevertheless, the distinction into two separate species remains questionable and was not supported by our genetic data. The considerable variations observed in cultured strains suggest that physiological factors might influence cell contour and bias identification. Analyses of Dinophysis cultures from French coastal waters using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) revealed species-conserved toxin profiles for D. acuta (dinophysistoxin 2 (DTX2), okadaic acid (OA), pectenotoxin 2 (PTX2)), D. caudata (PTX2) and D. tripos (PTX2), irrespective of geographical origin (Atlantic Ocean or Mediterranean Sea). Within the D. acuminata-complex, two different toxin profiles were observed: the strains of 'D. acuminata' (subclade A) from the English Channel and the Atlantic Ocean contained only OA while strains of 'D. sacculus' (subclade B) from Mediterranean Sea/Atlantic Ocean contained PTX2 as the dominant toxin, with OA and C9-esters also being present, albeit in lower proportions. The same difference in toxin profiles between 'D. sacculus' and 'D. acuminata' was reported in several studies from Galicia (NW- Spain). This difference in toxin profiles has consequences in terms of public health, and consequently for monitoring programs. While toxin profile could appear as a reliable feature separating 'D. acuminata' from 'D. sacculus' on both French and Spanish coasts, this does not seem consistent with observations on a broader geographical scale for the D. acuminata complex, possibly due to the frequent lack of genetic characterization.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

Strong constitutive expression divergence among strains but no evidence of differential expression associated with sexual reproduction in Alexandrium minutum

Sexual reproduction remains poorly characterized in dinoflagellates. This is especially the case at the molecular level. Here crossing experiments were performed among strains of the toxic dinoflagellate Alexandrium minutum belonging to two genetically divergent groups. Gene expression was compared between sexually compatible and incompatible crosses at the time of gamete fusion and resting cyst (~zygote) formation. Not a single transcript was identified as differentially expressed between compatible and incompatible crosses at these two crucial time points of the dinoflagellate life cycle. However, several thousands of transcripts displayed constitutive expression differences between strains. This was especially the case between the strains belonging to the genetically divergent groups. A few hundreds of transcripts were also identified as differentially expressed between strains belonging to opposite mating types. Some of these transcripts displayed homology with the SxtA protein, known to be involved in saxitoxin production in cyanobacteria, as well as with proteins potentially involved in mating in fungi.

Alexandrium pacificum and Alexandrium minutum: Harmful or environmentally friendly?

Alexandrium minutum and Alexandrium pacificum are representatives of the dinoflagellate genus that regularly proliferate on the French coasts and other global coastlines. These harmful species may threaten shellfish harvest and human health due to their ability to synthesize neurotoxic alkaloids of the saxitoxin group. However, some dinoflagellates such as A. minutum, and as reported here A. pacificum as well, may also have a beneficial impact on the environment by producing dimethylsulfoniopropionate-DMSP, the precursor of dimethylsulfur-DMS and sulfate aerosols involved in climate balance. However, environmental conditions might influence Alexandrium physiology towards the production of harmful or environmentally friendly compounds. After assessing the influence of two salinity regimes (33 and 38) relative to each species origin (Atlantic French coast and Mediterranean Lagoon respectively), it appears that DMSP and toxin content was variable between the three experimented strains and that higher salinity disadvantages toxin production and tends to favor the production of the osmolytes DMSP and glycine betaine. Hence, this key metabolite production is strain and species-dependent and is influenced by environmental conditions of salinity which in turn, can diversely affect the environment. Widespread coastal blooms of A. minutum and A. pacificum, although being a risk for seafood contamination with toxins, are also a DMSP and DMS source that potentially contribute to the ecosystem structuration and climate. Regarding recent advances in DMSP biosynthesis pathway, 3 dsyB homologs were found in A. minutum but no homolog of the diatom sequence TpMMT.

Production of Paralytic Shellfish Toxins (PSTs) in Toxic Alexandrium catenella is Intertwined with Photosynthesis and Energy Production

To investigate the mechanism for the production of paralytic shellfish toxins (PST) in toxic dinoflagellates, with a 2D-gel based approach, we had made two sets of proteomic comparisons: (a) between a toxic Alexandrium catenella (AC-T) and a phylogenetically closely related non-toxic strain (AC-N), (b) between toxic AC-T grown in a medium with 10% normal amount of phosphate (AC-T-10%P) known to induce higher toxicity and AC-T grown in normal medium. We found that photosynthesis and energy production related proteins were up-regulated in AC-T when compared to AC-N. However, the same group of proteins was down-regulated in AC-T-10%P when compared to normal AC-T. Examining the relationship of photosynthesis and toxin content of AC-T upon continuous photoperiod experiment revealed that while growth and associated toxin content increased after 8 days of continuous light, toxin content maintained constant when cells were shifted from continuous light to continuous dark for 3 days. This emphasized the cruciality of light availability on toxin biosynthesis in AC-T, while another light-independent mechanism may be responsible for higher toxicity in AC-T-10%P compared to normal AC-T. Taken all together, it is believed that the interplay between “illumination”, “photosynthesis”, “phosphate availability”, and “toxin production” is much more complicated than what we had previously anticipated.

Error, noise and bias in de novo transcriptome assemblies

De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.

Species specific gene expression dynamics during harmful algal blooms

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Biosynthesis of Saxitoxin in Marine Dinoflagellates: An Omics Perspective

Saxitoxin is an alkaloid neurotoxin originally isolated from the clam Saxidomus giganteus in 1957. This group of neurotoxins is produced by several species of freshwater cyanobacteria and marine dinoflagellates. The saxitoxin biosynthesis pathway was described for the first time in the 1980s and, since then, it was studied in more than seven cyanobacterial genera, comprising 26 genes that form a cluster ranging from 25.7 kb to 35 kb in sequence length. Due to the complexity of the genomic landscape, saxitoxin biosynthesis in dinoflagellates remains unknown. In order to reveal and understand the dynamics of the activity in such impressive unicellular organisms with a complex genome, a strategy that can carefully engage them in a systems view is necessary. Advances in omics technology (the collective tools of biological sciences) facilitated high-throughput studies of the genome, transcriptome, proteome, and metabolome of dinoflagellates. The omics approach was utilized to address saxitoxin-producing dinoflagellates in response to environmental stresses to improve understanding of dinoflagellates gene–environment interactions. Therefore, in this review, the progress in understanding dinoflagellate saxitoxin biosynthesis using an omics approach is emphasized. Further potential applications of metabolomics and genomics to unravel novel insights into saxitoxin biosynthesis in dinoflagellates are also reviewed.

Uptake of Inorganic and Organic Nitrogen Sources by Dinophysis acuminata and D. acuta

Dinoflagellate species of Dinophysis are obligate mixotrophs that require light, nutrients, and prey for sustained growth. Information about their nitrogenous nutrient preferences and their uptake kinetics are scarce. This study aimed to determine the preferred nitrogen sources in cultures of D. acuminata and D. acuta strains from the Galician Rías Baixas (NW Spain) and to compare their uptake kinetics. Well-fed versus starved cultures of D. acuminata and D. acuta were supplied with N¹⁵ labeled inorganic (nitrate, ammonium) and organic (urea) nutrients. Both species showed a preference for ammonium and urea whereas uptake of nitrate was negligible. Uptake rates by well-fed cells of D. acuminata and D. acuta were 200% and 50% higher, respectively, than by starved cells. Uptake of urea by D. acuminata was significantly higher than that of ammonium in both nutritional conditions. In contrast, similar uptake rates of both compounds were observed in D. acuta. The apparent inability of Dinophysis to take up nitrate suggests the existence of incomplete nitrate-reducing and assimilatory pathways, in line with the paucity of nitrate transporter homologs in the D. acuminata reference transcriptome. Results derived from this study will contribute to understand Harmful Algal Blooms succession and differences in the spatio-temporal distribution of the two Dinophysis species when they co-occur in stratified scenarios.

A review of the global distribution of Alexandrium minutum (Dinophyceae) and comments on ecology and associated paralytic shellfish toxin profiles, with a focus on Northern Europe

Alexandrium minutum is a globally distributed harmful algal bloom species with many strains that are known to produce paralytic shellfish toxins (PSTs) and consequently represent a concern to human and ecosystem health. This review highlights that A. minutum typically occurs in sheltered locations, with cell growth occurring during periods of stable water conditions. Sediment characteristics are important in the persistence of this species within a location, with fine sediments providing cyst deposits for ongoing inoculation to the water column. Toxic strains of A. minutum do not produce a consistent toxin profile, different populations produce a range of PSTs in differing quantities. Novel cluster analysis of published A. minutum toxin profiles indicates five PST profile clusters globally. Some clusters are grouped geographically (Northern Europe) while others are widely spread. Isolates from Taiwan have a range of toxin profile clusters and this area appears to have the most diverse set of PST producing A. minutum populations. These toxin profiles indicate that within the United Kingdom there are two populations of A. minutum grouping with strains from Northern France and Southern Ireland. There is a degree of interconnectivity in this region due to oceanic circulation and a high level of shipping and recreational boating. Further research into the interrelationships between the A. minutum populations in this global region would be of value.

Robust Dinoflagellata phylogeny inferred from public transcriptome databases

Dinoflagellates are dominant members of the plankton and play key roles in ocean ecosystems as primary producers, predators, parasites, coral photobionts, and causative agents of algal blooms that produce toxins harmful to humans and commercial fisheries. These unicellular protists exhibit remarkable trophic and morphological diversity and include species with some of the largest reported nuclear genomes. Despite their high ecological and economic importance, comprehensive genome (or transcriptome) based dinoflagellate trees of life are few in number. To address this issue, we used recently generated public sequencing data, including from the Moore Microbial Eukaryote Transcriptome Sequencing Project, to identify dinoflagellate-specific ortholog groups. These orthologs were combined to create a broadly sampled and highly resolved phylogeny of dinoflagellates. Our results emphasize the scope and utility of public sequencing databases in creating broad and robust phylogenies for large and complex taxonomic lineages, while also providing unique insights into the evolution of thecate dinoflagellates.

Can genomic data alone tell us whether speciation happened with gene flow?

The allopatric model, which requires a period of geographical isolation for speciation to complete, has been the standard model in the modern era. Recently, "speciation with gene flow" has been widely discussed in relation to the model of "strict allopatry" and the level of DNA divergence across genomic regions. We wish to caution that genomic data by themselves may only permit the rejection of the simplest form of allopatry. Even a slightly more complex and realistic model that starts with subdivided populations would be impossible to reject by the genomic data alone. To resolve this central issue of speciation, other forms of observations such as the sequencing of reproductive isolation genes or the identification of geographical barrier(s) will be necessary.