Global patterns of pelagic dinoflagellate diversity across protist size classes unveiled by metabarcoding

Molecular Identification of Gambierdiscus and Fukuyoa (Dinophyceae) from Environmental Samples

Ciguatera Fish Poisoning (CFP) is increasing across the Pacific and the distribution of the causative dinoflagellates appears to be expanding. Subtle differences in thecal plate morphology are used to distinguish dinoflagellate species, which are difficult to determine using light microscopy. For these reasons we sought to develop a Quantitative PCR assay that would detect all species from both Gambierdiscus and Fukuyoa genera in order to rapidly screen environmental samples for potentially toxic species. Additionally, a specific assay for F. paulensis was developed as this species is of concern in New Zealand coastal waters. Using the assays we analyzed 31 samples from three locations around New Zealand and the Kingdom of Tonga. Fourteen samples in total were positive for Gambierdiscus/Fukuyoa and two samples were also positive using the F. paulensis assay. Samples from the Kermadec Islands were further characterized using high-throughput sequencing metabarcoding. The majority of reads corresponded to Gambierdiscus species with three species identified at all sites (G. australes, G. honu and G. polynesiensis). This is the first confirmed identification of G. polynesiensis, a known ciguatoxin producer, in New Zealand waters. Other known toxin-producing genera were also detected, included Alexandrium, Amphidinium, Azadinium, Dinophysis, Ostreopsis, and Prorocentrum.

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.

De novo assembly and characterization of the transcriptome of the newly described dinoflagellate Ansanella granifera: Spotlight on flagellum-associated genes

Many dinoflagellates are known to cause red tides and often outgrow non-motile diatoms and motile small flagellates through active vertical migration between well-lit surface and eutrophic deep waters and/or by locating and ingesting prey cells. Their flagella play important roles in these two critical behaviors. However, the structural and functional genes of dinoflagellate flagella are very little known. Thus, a de novo assembly and characterization of the transcriptome of the fast-swimming dinoflagellate Ansanella granifera were conducted and its flagellum genes were compared with those of other dinoflagellates, motile small flagellates, and non-motile protist species. Based on assembled data using Trinity/CLC combined strategy, 83,652 transcripts of A. granifera were identified. The assembled consensus sequences were annotated to the NCBI non-redundant (nr), InterProScan, Gene Ontology (GO), and KEGG pathway analyses. Moreover, 71 structural and 35 functional flagellum-associated genes expressed were identified. The number of expressed flagellar structural and functional genes of A. granifera was not markedly different from those of other dinoflagellates or motile small flagellates, but much greater than those of non-motile species. Furthermore, in both phylogenetic trees based on the outer dynein arm (ODA1, ODA9, and DLC1) and inner dynein arm (IDA4, IDA7, and BOP5) flagellum genes of dinoflagellates, the problem of the long-branch attraction artifacts of Oxyrrhis marina which has been reported in the phylogenetic trees based on ribosomal DNA was removed. Moreover, in both phylogenetic trees based on the ODA and IDA flagellum genes, the species in the order Peridiniales or Gymnodiniales were revealed to belong to a big clade of each order. Therefore, the phylogenetic tree based on the flagellum genes is likely to give a clue to resolve the problem of separation in a big clade of a dinoflagellate order which has also been reported in the phylogenetic trees based on ribosomal DNA.

Molecular and phylogenetic characterization of Ostreopsis (Dinophyceae) and the description of a new species, Ostreopsis rhodesae sp. nov., from a subtropical Australian lagoon

Cryptic and pseudo-cryptic species are common amongst marine phytoplankton, and may cause misleading inferences of ecological and physiological data of plankton community studies. Deciphering the diversity and distribution of species of the benthic dinoflagellate Ostreopsis is one example, as there are many morphologically indistinct clades that differ greatly genetically and toxicologically from one another. In this study, a new species, Ostreopsis rhodesae from the southern Great Barrier Reef was described. While it initially appeared to be highly similar to several other Ostreopsis species, we found O. rhodesae can be distinguished based on the relative size of the second apical plate (2'), which is twice as long as the APC plate, and separates the third apical (3') from the third precingular (3'') plate. Phylogenetic trees based on the SSU, ITS/5.8S and D1-D2 and D8-D10 regions of the LSU rRNA were well supported, and showed a clear difference to other Ostreopsis clades. Compensatory base changes (CBCs) were identified in helices of the ITS2 between O. rhodesae and O. cf. ovata and O. cf. siamensis, which were also present in the same habitat. Fish gill cell lines were toxic to O. rhodesae, cell extracts but no palytoxin-like analogues were found in them. The findings highlight a case of pseudo-cryptic speciation, found in sympatry with closely related and morphologically similar species, but biologically and functionally distinct.

Development of a cost-effective metabarcoding strategy for analysis of the marine phytoplankton community

We developed a cost-effective metabarcoding strategy to analyze phytoplankton community structure using the Illumina MiSeq system. The amplicons (404-411 bp) obtained by end-pairing of two reads were sufficiently long to distinguish algal species and provided barcode data equivalent to those generated with the Roche 454 system, but at less than 1/20th of the cost. The original universal primer sequences targeting the 23S rDNA region and the PCR strategy were both modified, and this resulted in higher numbers of eukaryotic algal sequences by excluding non-photosynthetic proteobacterial sequences supporting effectiveness of this strategy. The novel strategy was used to analyze the phytoplankton community structure of six water samples from the East/Japan Sea: surface and 50 m depths at coastal and open-sea sites, with collections in May and July 2014. In total, 345 operational taxonomic units (OTUs) were identified, which covered most of the prokaryotic and eukaryotic algal phyla, including Dinophyta, Rhodophyta, Ochrophyta, Chlorophyta, Streptophyta, Cryptophyta, Haptophyta, and Cyanophyta. This highlights the importance of plastid 23S primers, which perform better than the currently used 16S primers for phytoplankton community surveys. The findings also revealed that more efforts should be made to update 23S rDNA sequences as well as those of 16S in the databases. Analysis of algal proportions in the six samples showed that community structure differed depending on location, depth and season. Across the six samples evaluated, the numbers of OTUs in each phylum were similar but their relative proportions varied. This novel strategy would allow laboratories to analyze large numbers of samples at reasonable expense, whereas this has not been possible to date due to cost and time. In addition, we expect that this strategy will generate a large amount of novel data that could potentially change established methods and tools that are currently used in the realms of oceanography and marine ecology.