Insights into the harmful algal flora in northwestern Mediterranean coastal lagoons revealed by pyrosequencing metabarcodes of the 28S rRNA gene

Diversity of eukaryotic and prokaryotic microbiota revealed by metabarcoding in Neotropical floodplain lakes

The diversity of eukaryotic and prokaryotic communities has been assessed by morphological and genetic approaches, which are used to characterize the microbiota in different environments. Here, planktonic prokaryotic and eukaryotic communities of the Araguaia River, located in the Central region of Brazil, were analyzed based on metabarcoding analysis of rRNA genes to evaluate the diversity of these groups in tropical floodplain lakes. Also, we tested their spatial concordance throughout the Araguaia river. Water samples were collected from 8 floodplain lakes in Araguaia River. The 16S and 18S rRNA genes were amplified and sequenced using Illumina MiSeq. For eukaryotes, 34,242 merged reads were obtained and 225 distinct OTUs were delineated, of which 106 OTUs were taxonomically classified. For prokaryotes, 26,426 sequences were obtained and 351 OTUs were detected. Of them, 231 were classified in at least one taxonomic category. The most representative eukaryotes belonged to Ciliophora, Chlorophyta and Charophyta. The prokaryotic phylum with the most OTUs classified were Proteobacteria, Actinobacteria and Bacteroidetes. The lakes did not show spatial concordance when comparing the similarity between their microbiota. The knowledge of freshwater biodiversity using DNA sequencing for important rivers, such as Araguaia River, can improve microbiota inventories of tropical biodiversity hotspots.

Monitoring of benthic eukaryotic communities in two tropical coastal lagoons through eDNA metabarcoding: a spatial and temporal approximation

Tropical coastal lagoons are important ecosystems that support high levels of biodiversity and provide several goods and services. Monitoring of benthic biodiversity and detection of harmful or invasive species is crucial, particularly in relation to seasonal and spatial variation of environmental conditions. In this study, eDNA metabarcoding was used in two tropical coastal lagoons, Chacahua (CH) and Corralero (C) (Southern Mexican Pacific), to describe the benthic biodiversity and its spatial–temporal dynamics. The distribution of benthic diversity within the lagoons showed a very particular pattern evidencing a transition from freshwater to seawater. Although the two lagoon systems are similar in terms of the species composition of metazoans and microeukaryotes, our findings indicate that they are different in taxa richness and structure, resulting in regional partitioning of the diversity with salinity as the driving factor of community composition in CH. Harmful, invasive, non-indigenous species, bioindicators and species of commercial importance were detected, demonstrating the reach of this technique for biodiversity monitoring along with the continued efforts of building species reference libraries.

Diversity and distribution of harmful microalgae in the Gulf of Thailand assessed by DNA metabarcoding

Information on the diversity and distribution of harmful microalgae in the Gulf of Thailand is very limited and mainly based on microscopic observations. Here, we collected 44 water samples from the Gulf of Thailand and its adjacent water (Perhentian Island, Malaysia) for comparison in 2018. DNA metabarcoding was performed targeting the partial large subunit ribosomal RNA gene (LSU rDNA D1-D3) and the internal transcribed spacers (ITS1 and ITS2). A total of 50 dinoflagellate genera (made up of 72 species) were identified based on the LSU rDNA dataset, while the results of ITS1 and ITS2 datasets revealed 33 and 32 dinoflagellate genera comprising 69 and 64 species, respectively. Five potentially toxic Pseudo-nitzschia (Bacillariophyceae) species were detected, with four as newly recorded species in the water (Pseudo-nitzschia americana/brasilliana, Pseudo-nitzschia simulans/delicatissima, P. galaxiae and P. multistriata). The highest relative abundances of P. galaxiae and P. multistriata were found in Trat Bay and Chumphon (accounting for 0.20% and 0.06% of total ASVs abundance, respectively). Three paralytic shellfish toxin producing dinoflagellate species were detected: Alexandrium tamiyavanichii, Alexandrium fragae, and Gymnodinium catenatum. The highest abundance of A. tamiyavanichii was found in the surface sample of Chumphon (CHO7 station), accounting for 1.95% of total ASVs abundance. Two azaspiracid producing dinoflagellate species, Azadinium poporum ribotype B, Azadinium spinosum ribotype A, and a pinnatoxin producing dinoflagellate species Vulcanodinium rugosum, with two ribotypes B and C, were revealed from the datasets although with very low abundances. Six fish killing dinoflagellate species, including Margalefidinium polykrikoides group IV, Margalefidinium fulvescens, Karenia mikimotoi, Karenia selliformis ribotype B, Karlodinium australe, and Karlodinium digitatum were detected and all representing new records in this area. The findings of numerous harmful microalgal species in the Gulf of Thailand highlight the potential risk of human intoxication and fish killing events.

Environmental DNA Metabarcoding for Simultaneous Monitoring and Ecological Assessment of Many Harmful Algae

Harmful algae can have profound economic, environmental, and social consequences. As the timing, frequency, and severity of harmful algal blooms (HABs) change alongside global climate, efficient tools to monitor and understand the current ecological context of these taxa are increasingly important. Here we employ environmental DNA metabarcoding to identify patterns in a wide variety of potentially harmful algae and associated ecological communities in the Hood Canal of Puget Sound in Washington State, USA. Tracking trends of occurrence in a series of water samples over a period of 19 months, we find algal sequences from genera with harmful members in a majority of samples, suggesting that these groups are routinely present in local waters. We report patterns in variants of the economically important genus Pseudo-nitzschia (of which some members produce domoic acid; family Bacillariaceae), as well as multiple potentially harmful algal taxa previously unknown or poorly documented in the region, including a cold-water variant from the genus Alexandrium (of which some members produce saxitoxin; family Gonyaulacaceae), two variants from the genus Karlodinium (of which some members produce karlotoxins; family Kareniaceae), and one variant from the parasitic genus Hematodinium (family Syndiniaceae). We then use data on environmental variables and the biological community surrounding each algal taxon to illustrate the ecological context in which they are commonly found. Environmental DNA metabarcoding thus simultaneously (1) alerts us to potential new or cryptic occurrences of algae from harmful genera, (2) expands our knowledge of the co-occurring conditions and species associated with the growth of these organisms in changing marine environments, and (3) suggests a pathway for multispecies monitoring and management moving forward.

Development of an absolute quantification method for ribosomal RNA gene copy numbers per eukaryotic single cell by digital PCR

Recent increase of Harmful Algal Blooms (HAB) causes world-wide ecological, economical, and health issues, and more attention is paid to frequent coastal monitoring for the early detection of HAB species to prevent or reduce such impacts. Use of molecular tools in addition to traditional microscopy-based observation has become one of the promising methodologies for coastal monitoring. However, as ribosomal RNA (rRNA) genes are commonly targeted in molecular studies, variability in the rRNA gene copy number within and between species must be considered to provide quantitative information in quantitative PCR (qPCR), digital PCR (dPCR), and metabarcoding analyses. Currently, this information is only available for a limited number of species. The present study utilized a dPCR technology to quantify copy numbers of rRNA genes per single cell in 16 phytoplankton species, the majority of which are toxin-producers, using a newly developed universal primer set accompanied by a labeled probe with a fluorophore and a double-quencher. In silico PCR using the newly developed primers allowed the detection of taxa from 8 supergroups, demonstrating universality and broad coverage of the primer set. Chelex buffer was found to be suitable for DNA extraction to obtain DNA fragments with suitable size to avoid underestimation of the copy numbers. The study successfully demonstrated the first comparison of absolute quantification of 18S rRNA copy numbers per cell from 16 phytoplankton species by the dPCR technology.

Toxic marine microalgae and noxious blooms in the Mediterranean Sea: A contribution to the Global HAB Status Report

We review the spatial distribution of toxic marine microalgal species and the impacts of all types of harmful algal events (Harmful Algal Blooms, HABs) in the Mediterranean Sea (MS), including the Black Sea, the Sea of Marmara, coastal lagoons and transitional waters, based on two databases compiled in the Ocean Biogeographic Information System (OBIS). Eighty-four potentially toxic species have been detected in the MS (2,350 records), of which 16 described from these waters between 1860 and 2014 and a few suspected to have been introduced. More than half of these species (46) produce toxins that may affect human health, the remainders ichthyotoxic substances (29) or other types of toxins (9). Nevertheless, toxicity-related events are not frequent in the MS (308 records in 31 years), and mainly consist of impacts on aquaculture, caused by the dinoflagellates Dinophysis and Alexandrium, along with a few actual shellfish poisoning cases. Pseudo-nitzschia blooms are widespread, but domoic acid in shellfish rarely exceeds regulatory levels. Fish kills are probably less sporadic than reported, representing a problem at a few places along the southern MS coasts and in the Ebro River Delta. Since the last decade of the 20^th century, blooms of the benthic dinoflagellates Ostreopsis cf. ovata have regularly occurred all along rocky shores of the MS, at times with human health problems caused by toxic aerosol. New records of Gambierdiscus and Fukuyoa, until now reported for the westernmost and easternmost MS coasts, raise concerns about the risk of ciguatera, a syndrome so far known only for subtropical and tropical areas. Recent discoveries are the dinoflagellates Vulcanodinium rugosum, responsible for the presence of pinnatoxins in French lagoons' shellfish, and the azaspiracid-producers Azadinium spp. Mucilages and discolorations have a major impact on tourism in summer. Reports of toxic species and HABs have apparently increased in the MS over the last half century, which is likely related to the increased awareness and monitoring operations rather than to an actual increase of these phenomena. Indeed, while the case of Ostreopsis appears as a sudden upsurge rather than a trend, no actual increase of toxic or noxious events has so far emerged in intensively studied areas, such as the French and Spanish coasts or the Adriatic Sea. Moreover, some cases of decrease are reported, e.g., for Alexandrium minutum blooms disappearing from the Harbour of Alexandria. Overall, main HAB risks derive from cases of massive development of microalgal biomass and consequent impacts of reduced coastal water quality on tourism, which represents the largest part of the marine economy along the MS coasts.

Seasonal and Spatial Variations in Bacterial Communities From Tetrodotoxin-Bearing and Non-tetrodotoxin-Bearing Clams

Tetrodotoxin (TTX) is one of the most potent naturally occurring compounds and is responsible for many human intoxications worldwide. Paphies australis are endemic clams to New Zealand which contain varying concentrations of TTX. Research suggests that P. australis accumulate the toxin exogenously, but the source remains uncertain. The aim of this study was to identify potential bacterial TTX-producers by exploring differences in bacterial communities in two organs of P. australis: the siphon and digestive gland. Samples from the digestive glands of a non-toxic bivalve Austrovenus stutchburyi that lives amongst toxic P. australis populations were also analyzed. Bacterial communities were characterized using 16S ribosomal RNA gene metabarcoding in P. australis sourced monthly from the Hokianga Harbor, a site known to have TTX-bearing clams, for 1 year, from ten sites with varying TTX concentrations around New Zealand, and in A. stutchburyi from the Hokianga Harbor. Tetrodotoxin was detected in P. australis from sites all around New Zealand and in all P. australis collected monthly from the Hokianga Harbor. The toxin averaged 150 μg kg^-1 over the year of sampling in the Hokianga Harbor but no TTX was detected in the A. stutchburyi samples from the same site. Bacterial species diversity differed amongst sites (p < 0.001, F = 5.9) and the diversity in siphon samples was significantly higher than in digestive glands (p < 0.001, F = 65.8). Spirochaetaceae (4-60%) and Mycoplasmataceae (16-78%) were the most abundant families in the siphons and the digestive glands, respectively. The bacterial communities were compared between sites with the lowest TTX concentrations and the Hokianga Harbor (site with the highest TTX concentrations), and the core bacterial communities from TTX-bearing individuals were analyzed. The results from both spatial and temporal studies corroborate with previous hypotheses that Vibrio and Bacillus could be responsible for the source of TTX in bivalves. The results from this study also indicate that marine cyanobacteria, in particular picocyanobacteria (e.g., Cyanobium, Synechococcus, Pleurocapsa, and Prochlorococcus), should be investigated further as potential TTX producers.

Seasonal variations of phytoplankton assemblages in relation to environmental factors in Mediterranean coastal waters of Morocco, a focus on HABs species

Studies on phytoplankton and in particular Harmful Algal Blooms (HABs) species in southern Mediterranean waters are scarce. We performed from April 2008 to June 2009 weekly investigations on microphytoplankton community structure and abundance in two contrasted marine ecosystems located in the western Moroccan Mediterranean coast, M'diq Bay and Oued Laou Estuary. Simultaneously, we measured the main physico-chemical parameters. Globally, the two studied areas showed comparable values of the assessed abiotic environmental factors. Temperature and salinity followed seasonal variation with values ranging from 13.5 °C to 21.4 °C and 31 to 36.8, respectively. Average nutrient values in surface water ranged from 0.7 to 45.76 μM for dissolved inorganic nitrogen, 0.02-2.10 μM for PO₄ and 0.23-17.46 μM for SiO₄ in the study areas. A total of 92 taxa belonging to 8 taxonomic classes were found. The highest number of microphytoplankton abundance reached 1.2 × 10⁶ cells L^-1 with diatoms being the most abundant taxa. Factorial Discriminant Analysis (FDA) and Spearman correlation test showed a significant seasonal discrimination of dominant microphytoplankton species. These micro-organisms were associated with different environmental variables, in particular temperature and salinity. Numerous HABs species were encountered regularly along the year. Although Dinophysis species and Prorocentrum lima were present in both sites, no Lipophilic Shellfish Poisoning was detected for the analyzed bivalve mollusks. Domoic acid (DA), produced by toxic species of Pseudo-nitzschia was found with concentrations up to 18 µg DA g^-1 in the smooth clam Callista chione. Data showed that the observed persistent and dramatic Paralytic Shellfish Poisoning (PSP) intoxication of mollusks resulted probably of Gymnodinium catenatum proliferations in both studied areas. Contrary to C. chione, the cockle Achanthocardia tuberculatum showed a permanent and extremely high toxicity level during the 15 months survey with up to 7545 µg Equivalent Saxitoxin kg^-1 flesh (ten times higher than the sanitary threshold of 800 µg eqSTX Kg^-1flesh). The present work highlights for the first time the dynamic of microphytoplankton including HABs species and their associated toxin accumulation in the commercially exploited shellfish in the southern western Mediterranean waters of Morocco. Furthermore, the acquired data will help us to improve the monitoring of HABs species and related toxins in these coastal marine systems.

Patterns in Alpha and Beta Phytoplankton Diversity along a Conductivity Gradient in Coastal Mediterranean Lagoons

Understanding the diversity patterns of phytoplankton assemblages in coastal lagoons is clearly important for water management. In this study, we explored alpha and beta diversity patterns in phytoplankton communities across five Mediterranean lagoons hydrologically connected to Vistonikos Gulf. We examined the phytoplankton community composition and biomass on a monthly basis from November 2018 to October 2019. For this, water samples were collected from seven inshore, brackish and coastal waters, sampling sites covering a wide range of conductivity. We found significant spatial and temporal differences in phytoplankton alpha diversity and in phytoplankton biomass metrics explained by the high variation of conductivity. Evenness remained low throughout the study period, reflecting significant dominance of several phytoplankton blooms. Harmful algal blooms of Prorocentrum minimum, Alexandrium sp., Rhizosolenia setigera and Cylindrotheca closterium occurred. The system’s species pool was characterized by relatively high phytoplankton beta diversity (average ~0.7) resulting from high temporal species turnover (90%). Overall, alpha and beta diversity components were indicative of rather heterogeneous phytoplankton communities which were associated with the high differences in conductivity among the sampling sites.

Comparing metabarcoding and morphological approaches to identify phytoplankton taxa associated with harmful algal blooms

Molecular techniques are expected to be highly useful in detecting taxa causing harmful algal blooms (HABs). This is the first report in Canada evaluating HABs-related species identification using a combination of morphological and molecular approaches. Microscopy, quantitative polymerase chain reaction (qPCR), and metabarcoding with multiple markers (i.e., 16S, 18S-dinoflagellate and 18S-diatom, large subunit (28S) rDNA) were applied on samples ( n = 54) containing suspected harmful algae (e.g., Alexandrium spp., Chattonella sp., Chrysochromulina spp., Dictyocha spp., Heterosigma akashiwo, Protoceratium reticulatum, Pseudochattonella verruculosa, Pseudo-nitzschia spp., Pseudopedinella sp.). Owing to methodology limitations, qPCR result interpretation was limited, although good detectability occurred using previously published assays for Alexandrium tamarense, H. akashiwo, and P. verruculosa. Overall, the multiple-marker metabarcoding results were superior to the morphology-based methods, with the exception of taxa from the silicoflagellate group. The combined results using both 18S markers and the 28S marker together closely corresponded with morphological identification of targeted species, providing the best overall taxonomic coverage and resolution. The most numerous unique taxa were identified using the 18S-dinoflagellate amplicon, and the best resolution to the species level occurred using the 28S amplicon. Molecular techniques are therefore promising for HABs taxa detection but currently depend on deploying multiple markers for metabarcoding.

Evaluation of FluoroProbe® performance for the phytoplankton-based assessment of the ecological status of Mediterranean coastal lagoons

The European Water Framework Directive and several other legislations worldwide have selected phytoplankton for monitoring the ecological status of surface waters. This assessment is a complicated task in coastal lagoons due to their intrinsic variability, prompting moves to use real-time measurements. Here, we tested the ability of the submersible spectrofluorometer FluoroProbe® to accurately estimate the phytoplankton biomass and to efficiently discriminate spectral groups in Mediterranean coastal lagoons, by using sub-surface water samples (n = 107) collected at Biguglia lagoon (Corsica) in different environmental situations (salinity and trophic state) from March 2012 to December 2014. We compared the estimates of biomass and phytoplankton group composition obtained with the FluoroProbe® (in situ and lab measurements) with the spectrofluorimetrically measured biomass and HPLC-derived quantifications of pigment concentrations. FluoroProbe® provided good estimates of the total phytoplankton biomass (particularly, the lab measurements). The FluoroProbe® data were significantly correlated with the HPLC results, except for the in situ measurements of very weak concentrations of blue-green and red algae. Our findings indicate that factory-calibrated FluoroProbe® is an efficient and easy-to-use real-time phytoplankton monitoring tool in coastal lagoons, especially as an early warning system for the detection of potentially harmful algal blooms. Practical instructions dedicated to non-specialist field operators are provided. A simple and efficient method for discarding in situ measurement outliers is also proposed.

Pseudo-nitzschia, Nitzschia, and domoic acid: New research since 2011

Some diatoms of the genera Pseudo-nitzschia and Nitzschia produce the neurotoxin domoic acid (DA), a compound that caused amnesic shellfish poisoning (ASP) in humans just over 30 years ago (December 1987) in eastern Canada. This review covers new information since two previous reviews in 2012. Nitzschia bizertensis was subsequently discovered to be toxigenic in Tunisian waters. The known distribution of N. navis-varingica has expanded from Vietnam to Malaysia, Indonesia, the Philippines and Australia. Furthermore, 15 new species (and one new variety) of Pseudo-nitzschia have been discovered, bringing the total to 52. Seven new species were found to produce DA, bringing the total of toxigenic species to 26. We list all Pseudo-nitzschia species, their ability to produce DA, and show their global distribution. A consequence of the extended distribution and increased number of toxigenic species worldwide is that DA is now found more pervasively in the food web, contaminating new marine organisms (especially marine mammals), affecting their physiology and disrupting ecosystems. Recent findings highlight how zooplankton grazers can induce DA production in Pseudo-nitzschia and how bacteria interact with Pseudo-nitzschia. Since 2012, new discoveries have been reported on physiological controls of Pseudo-nitzschia growth and DA production, its sexual reproduction, and infection by an oomycete parasitoid. Many advances are the result of applying molecular approaches to discovering new species, and to understanding the population genetic structure of Pseudo-nitzschia and mechanisms used to cope with iron limitation. The availability of genomes from three Pseudo-nitzschia species, coupled with a comparative transcriptomic approach, has allowed advances in our understanding of the sexual reproduction of Pseudo-nitzschia, its signaling pathways, its interactions with bacteria, and genes involved in iron and vitamin B₁₂ and B₇ metabolism. Although there have been no new confirmed cases of ASP since 1987 because of monitoring efforts, new blooms have occurred. A massive toxic Pseudo-nitzschia bloom affected the entire west coast of North America during 2015-2016, and was linked to a 'warm blob' of ocean water. Other smaller toxic blooms occurred in the Gulf of Mexico and east coast of North America. Knowledge gaps remain, including how and why DA and its isomers are produced, the world distribution of potentially toxigenic Nitzschia species, the prevalence of DA isomers, and molecular markers to discriminate between toxigenic and non-toxigenic species and to discover sexually reproducing populations in the field.

dinoref: A curated dinoflagellate (Dinophyceae) reference database for the 18S rRNA gene

Dinoflagellates are a heterogeneous group of protists present in all aquatic ecosystems where they occupy various ecological niches. They play a major role as primary producers, but many species are mixotrophic or heterotrophic. Environmental metabarcoding based on high-throughput sequencing is increasingly applied to assess diversity and abundance of planktonic organisms, and reference databases are definitely needed to taxonomically assign the huge number of sequences. We provide an updated 18S rRNA reference database of dinoflagellates: dinoref. Sequences were downloaded from genbank and filtered based on stringent quality criteria. All sequences were taxonomically curated, classified taking into account classical morphotaxonomic studies and molecular phylogenies, and linked to a series of metadata. dinoref includes 1,671 sequences representing 149 genera and 422 species. The taxonomic assignation of 468 sequences was revised. The largest number of sequences belongs to Gonyaulacales and Suessiales that include toxic and symbiotic species. dinoref provides an opportunity to test the level of taxonomic resolution of different 18S barcode markers based on a large number of sequences and species. As an example, when only the V4 region is considered, 374 of the 422 species included in dinoref can still be unambiguously identified. Clustering the V4 sequences at 98% similarity, a threshold that is commonly applied in metabarcoding studies, resulted in a considerable underestimation of species diversity.