An assessment on the intrapopulational and intraindividual genetic diversity in LSU rDNA in the harmful algal blooms-forming dinoflagellate Margalefidinium (= Cochlodinium) fulvescens based on clonal cultures and bloom samples from Jiaozhou Bay, China

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

eDNA metabarcoding reveals biodiversity and depth stratification patterns of dinoflagellate assemblages within the epipelagic zone of the western Coral Sea

BACKGROUND

Dinoflagellates play critical roles in the functioning of marine ecosystems but also may pose a hazard to human and ecosystem health by causing harmful algal blooms (HABs). The Coral Sea is a biodiversity hotspot, but its dinoflagellate assemblages in pelagic waters have not been studied by modern sequencing methods. We used metabarcoding of the 18 S rRNA V4 amplicon to assess the diversity and structure of dinoflagellate assemblages throughout the water column to a depth of 150 m at three stations in the Western Coral Sea. Additionally, at one station we compared metabarcoding with morphological methods to optimise identification and detection of dinoflagellates.

RESULTS

Stratification of dinoflagellate assemblages was evident in depth-specific relative abundances of taxonomic groups; the greatest difference was between the 5-30 m assemblages and the 130-150 m assemblages. The relative abundance of Dinophyceae (photosynthetic and heterotrophic) decreased with increasing depth, whereas that of Syndiniales (parasitic) increased with increasing depth. The composition of major taxonomic groups was similar among stations. Taxonomic richness and diversity of amplicon sequence variants (ASVs) were similar among depths and stations; however, the abundance of dominant taxa was highest within 0-30 m, and the abundance of rare taxa was highest within 130-150 m, indicating adaptations to specific depth strata. The number of unclassified ASVs at the family and species levels was very high, particularly for Syndinian representatives.

CONCLUSIONS

Dinoflagellate assemblages in open water of the Coral Sea are highly diverse and taxonomically stratified by depth; patterns of relative abundance along the depth gradient reflect environmental factors and ecological processes. Metabarcoding detects more species richness than does traditional microscopical methods of sample analysis, yet the methods are complementary, with morphological analysis revealing additional richness. The large number of unclassified dinoflagellate-ASVs indicates a need for improved taxonomic reference databases and suggests presence of dinoflagellate-crypto and-morphospecies.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

Diversity and Distribution of Harmful Algal Bloom Species from Seamount to Coastal Waters in the South China Sea

Mount Xianbei is one of the largest shallow seamounts located in the middle of the South China Sea (SCS), which might play a role in shaping the biodiversity of surrounding continental coastal waters, particularly the diversity of phytoplankton species causing frequent harmful algal blooms (HABs) in northern SCS. However, the diversity, composition, and distribution of phytoplankton species in the seamount regions of Xianbei remain largely unexplored. In this study, samples around and outside the seamount regions were collected during a late summer cruise of 2021 to test whether seamounts play a role in HAB species propagation. In total, we identified 19 HAB species across all samples using the ASV-based DNA metabarcoding approach, 6 of which had not been reported previously in the SCS, suggesting a diverse HAB species in the SCS. Specifically, 16 HAB species were found in the seamount region of Xianbei, and 5 of them were also found in the coastal waters, indicating a close connection between seamount and coastal waters. This study was the first attempt to explore HAB species' spatial diversity and vertical distribution in the seamount region of Xianbei at single-nucleotide resolution, which provides a novel explanation for the coastal HAB occurrence in the northern SCS. IMPORTANCE There are a number of seamounts under the water of the South China Sea (SCS). The seamounts might play a role in shaping the biodiversity of surrounding continental coastal waters. However, there is no direct evidence revealing the relationship of the biodiversity of phytoplankton between seamounts and coastal waters in the SCS, especially those species having the potential to form harmful algal blooms (HABs). Some HAB species might proliferate in certain geographic locations, while others may be broadly distributed across oceanic provinces. In this study, we provided a detailed analysis of phytoplankton composition and molecular detection of HAB species from seamount to coastal waters in the SCS, which suggested a strong interaction in the HAB species between the two areas. This finding provides new insights into the diversity and distribution of HABs in seamounts and their role in shaping the composition and the occurrence of HABs in coastal water.

DNA barcoding, an effective tool for species identification: a review

DNA barcoding is a powerful taxonomic tool to identify and discover species. DNA barcoding utilizes one or more standardized short DNA regions for taxon identification. With the emergence of new sequencing techniques, such as Next-generation sequencing (NGS), ONT MinION nanopore sequencing, and Pac Bio sequencing, DNA barcoding has become more accurate, fast, and reliable. Rapid species identification by DNA barcodes has been used in a variety of fields, including forensic science, control of the food supply chain, and disease understanding. The Consortium for Barcode of Life (CBOL) presents various working groups to identify the universal barcode gene, such as COI in metazoans; rbcL, matK, and ITS in plants; ITS in fungi; 16S rRNA gene in bacteria and archaea, and creating a reference DNA barcode library. In this article, an attempt has been made to analyze the various proposed DNA barcode for different organisms, strengths & limitations, recent advancements in DNA barcoding, and methods to speed up the DNA barcode reference library construction. This study concludes that constructing a reference library with high species coverage would be a major step toward identifying species by DNA barcodes. This can be achieved in a short period of time by using advanced sequencing and data analysis methods.

Geographical distribution, molecular and toxin diversity of the dinoflagellate species Gambierdiscus honu in the Pacific region

An increase in cases of ciguatera poisoning (CP) and expansion of the causative species in the South Pacific region highlight the need for baseline data on toxic microalgal species to help identify new areas of risk and manage known hot spots. Gambierdiscus honu is a toxin producing and potential CP causing dinoflagellate species, first described in 2017. Currently no high-resolution geographical distribution, intraspecific genetic variation or toxin production diversity data is available for G. honu. This research aimed to further characterize G. honu by investigating its distribution using species-specific real-time polymerase chain reaction assays at 25 sites in an area spanning ∼8000 km of the Coral Sea/Pacific Ocean, and assessing intraspecific genetic variation, toxicity and toxin production of isolated strains. Assessment of genetic variation of the partial rRNA operon of isolates demonstrated no significant intraspecific population structure, in addition to a lack of adherence to isolation by distance (IBD) model of evolution. The detected distribution of G. honu in the Pacific region was within the expected tropical to temperate latitudinal ranges of 10° to -30° and extended from Australia to French Polynesia. In the lipophilic fractions, the neuroblastoma cell-based assay (CBA-N2a) showed no ciguatoxin (CTX)-like activity for nine of the 10 isolates, and an atypical pattern for CAWD233 isolate which showed cytotoxic activity in OV- and OV+ conditions. In the same way, liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis confirmed no Pacific-CTXs (CTX-3B, CTX-3C, CTX-4A, CTX-4B) were produced by the ten strains. The CBA-N2a assessment of the hydrophilic fractions showed moderate to high cytotoxicity in both OV- and OV+ condition for all the strains showing a cytotoxic profile similar to that of gambierone. Indeed, this study is the first to show the cytotoxic activity of gambierone on mouse neuroblastoma cells while no cytotoxicity was observed when 44-MG was analysed at the same concentrations using the CBA-N2a. Analysis of the hydrophilic via LC-MS/MS confirmed production of gambierone in all isolates, ranging from 2.1 to 38.1 pg/cell, with 44-methylgambierone (44-MG) also produced by eight of the isolates, ranging from 0.3 to 42.9 pg/cell. No maitotoxin-1 was detected in any of the isolates. Classification of the G. honu strains according to the quantities of gambierone produced aligned with the classification of their cytotoxicity using the CBA-N2a. Finally, no maitotoxin-1 (MTX) was detected in any of the isolates. This study shows G. honu is widely distributed within the Pacific region with no significant intraspecific population structure present. This aligns with the view of microalgal populations as global metapopulations, however more in-depth assessment with other genetic markers could detect further structure. Toxicity diversity across 10 isolates assessed did not display any geographical patterns.

Biodiversity and Spatial-Temporal Dynamics of Margalefidinium Species in Jiaozhou Bay, China

Many Margalefidinium species are cosmopolitan harmful algal bloom (HAB) species that have caused huge economic and ecological damage. Despite extensive research on Margalefidinium species, the biodiversity and spatial-temporal dynamics of these species remain obscure. Jiaozhou Bay is an ideal area for HAB research, being one of the earliest marine survey areas in China. In this study, we carried out the first metabarcoding study on the temporal and spatial dynamics of Margalefidinium species using the 18S rDNA V4 region as the molecular marker and samples collected monthly at 12 sampling sites in Jiaozhou Bay in 2019. Two harmful Margalefidinium species (M. polykrikoides and M. fulvescens) were identified with potentially high genetic diversity (although we cannot rule out the possibility of intra-genome sequence variations). Both M. polykrikoides and M. fulvescens demonstrated strong temporal preference with a sharp peak of abundance in early autumn (September), but without showing strong location preference in Jiaozhou Bay. Our results revealed that temperature might be the main driver for their temporal dynamics. Knowledge of biodiversity and spatial-temporal dynamics of the Margalefidinium species may shed light on the understanding of mechanisms underlying strongly biased occurrences of Margalefidinium blooms recorded globally.

Dependence of genome size and copy number of rRNA gene on cell volume in dinoflagellates

Dinoflagellates are an ecologically important group of protists in aquatic environment and have evolved many unusual and enigmatic genomic features such as immense genome sizes, high repeated genes, and a large portion of hydroxymethyluracil in DNA. Although previous studies have observed positive correlations between the large subunit (LSU) rRNA gene copy number and genome size of a variety of eukaryotic organisms (e.g. higher plants and animals), or between cell volume and LSU rRNA gene copy number, and/or between genome size and cell size, which suggests a possible co-evolution among these three features in different lineages of life, it remains an open question regarding the relationships among these three parameters in dinoflagellates. For the first time, we estimated the copy numbers of the LSU rRNA gene, the genome sizes, and cell volumes within a broad range of dinoflagellates (covering 15 species of 11 genera) using single-cell qPCR-based assay (determining LSU rRNA gene copy number), FlowCAM (cell volume measurement), and ultraviolet spectrophotometry (genome size estimation). The measured copy number of LSU rRNA gene ranged from 398 ± 184 (Prorocentrum minimum) to 152,078 ± 33,555 copies•cell^-1 (Alexandrium pacificum), while the genome size and the cell volume ranged from 5.6 ± 0.2 (Karlodinium veneficum) to 853 ± 19.9 pg•cell^-1 (Pseliodinium pirum), and from 1,070 ± 225 (Kar. veneficum) to 168,474 ± 124,180 μm³ (Ps. pirum), respectively. Together with the three parameters measured in literature, there are significant positive linear correlations between LSU rRNA gene copy numbers and genome sizes, cell volumes and LSU rRNA gene copy numbers, and between genome sizes and cell volumes via comparisons of multi-model regression analyses, suggesting a dependence of genome size and rRNA gene copy number on the cell volumes of dinoflagellates. Validation of the measurement methods was conducted via comparisons between reported data in the literature and that predicted using the linear equations we obtained, and between genome size measured by flow cytometry (FCM) and ultraviolet spectrophotometry (Nanodrop). These results provide insightful understandings of dinoflagellate evolution in terms of the relationships among genomes, gene copy number, and cell volume, and of rRNA gene-based studies in intra-populational and intra-individual genetic diversity, taxonomy, and diversity assessment in the environment of dinoflagellates. The results also provide a dataset useful for reads calibration in environmental metabarcoding studies of dinoflagellates and selection of candidate species for whole genome sequencing.

Molecular characterization of harmful algal blooms in the Bohai Sea using metabarcoding analysis

Although the occurrences of harmful algal blooms (HABs) have been intensifying, many HABs in coastal waters may have been neglected despite their damaging impact directly on ecology and indirectly on human and animal health. The current detection of HABs depends primarily on the water coloration, chlorophyll intensity, cell density, and mortality due to HAB toxicity. Such methods may not be adequately sensitive to detecting HABs that are relatively transient or small scale. The Bohai Sea is the largest inlet of the Yellow Sea located on the northeast coast of China and famous for shipping and marine aquacultures. HABs frequently occur in the Bohai Sea. In this study, we explored the composition, diversity, and distribution of HAB species using the metabarcoding approach. Through sequencing and the analyzing the 18S rDNA V4 region of 15 samples collected from spatially isolated sites in the Bohai Sea during an expedition in the summer of 2019, we identified 74 potential HAB species including 34 that had not been reported in the Bohai Sea in previous studies. This project provided a detailed analysis of phytoplankton composition, and molecular detection of HAB species in the Bohai Sea. In particular, these analyses revealed extremely high relative abundances of the ichthyotoxic phytoplankton species Vicicitus globosus (Dictyochophyceae) at multiple adjacent sampling sites in the Bohai Bay, which were close to the Yellow River Estuary during the expedition. The results revealed the occurrence of a potential HAB event that would be otherwise undetected using conventional methods, highlighting the sensitivity and power of metabarcoding analysis in detecting HABs and HAB species. This research suggested the value for routine and long-term monitoring of HAB species as an approach for monitoring HABs.