qPCR Assays for the Detection and Quantification of Multiple Paralytic Shellfish Toxin-Producing Species of Alexandrium

Assessing methods for detecting Alexandrium catenella (Dinophyceae) and paralytic shellfish toxins in Southeast Alaska

Blooms of Alexandrium catenella threaten to disrupt subsistence, recreational, and commercial shellfish harvest in Alaska, as the paralytic shellfish toxins (PSTs) produced pose a serious public health risk and can lead to costly shutdowns for shellfish farmers. Current methods of PST detection in the region range from monitoring programs utilizing net tows to detect A. catenella to direct shellfish tissue testing via mouse bioassay (MBA) for commercial aquaculture harvest, as well as various optional testing methods for subsistence and recreational harvesters. The efficacy and feasibility of these methods vary, and they have not been directly compared in Southeast Alaska. In this study, we sought to assess and compare A. catenella and PST early detection methods to determine which can provide the most effective and accurate warning of A. catenella blooms or PST events. We found microscope counts to be variable and prone to missing lower numbers of A. catenella, which may be indicative of bloom formation. However, quantitative polymerase chain reaction (qPCR) significantly correlated with microscope counts and was able to effectively detect even low numbers of A. catenella on all sampling days. Paralytic shellfish toxin concentrations measured by enzyme-linked immunosorbent assay and MBA significantly correlated with each other, qPCR, and some microscope counts. These results show that qPCR is an effective tool for both monitoring A. catenella and serving as a proxy for PSTs. Further work is needed to refine qPCR protocols in this system to provide bloom warnings on an actionable timescale for the aquaculture industry and other shellfish harvesters. Integr Environ Assess Manag 2024;20:2189-2202. © 2024 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Biodiversity of microorganisms in the Baltic Sea: the power of novel methods in the identification of marine microbes

Until recently, the data on the diversity of the entire microbial community from the Baltic Sea were relatively rare and very scarce. However, modern molecular methods have provided new insights into this field with interesting results. They can be summarized as follows. (i) Although low salinity causes a reduction in the biodiversity of multicellular species relative to the populations of the North-East Atlantic, no such reduction occurs in bacterial diversity. (ii) Among cyanobacteria, the picocyanobacterial group dominates when considering gene abundance, while filamentous cyanobacteria dominate in means of biomass. (iii) The diversity of diatoms and dinoflagellates is significantly larger than described a few decades ago; however, molecular studies on these groups are still scarce. (iv) Knowledge gaps in other protistan communities are evident. (v) Salinity is the main limiting parameter of pelagic fungal community composition, while the benthic fungal diversity is shaped by water depth, salinity, and sediment C and N availability. (vi) Bacteriophages are the predominant group of viruses, while among viruses infecting eukaryotic hosts, Phycodnaviridae are the most abundant; the Baltic Sea virome is contaminated with viruses originating from urban and/or industrial habitats. These features make the Baltic Sea microbiome specific and unique among other marine environments.

An On-Farm Workflow for Predictive Management of Paralytic Shellfish Toxin-Producing Harmful Algal Blooms for the Aquaculture Industry

Paralytic shellfish toxins (PSTs) produced by marine dinoflagellates significantly impact shellfish industries worldwide. Early detection on-farm and with minimal training would allow additional time for management decisions to minimize economic losses. Here, we describe and test a standardized workflow based on the detection of sxtA4, an initial gene in the biosynthesis of PSTs. The workflow is simple and inexpensive and does not require a specialized laboratory. It consists of (1) water collection and filtration using a custom gravity sampler, (2) buffer selection for sample preservation and cell lysis for DNA, and (3) an assay based on a region of sxtA, DinoDtec lyophilized quantitative polymerase chain reaction (qPCR) assay. Water samples spiked with Alexandrium catenella showed a cell recovery of >90% when compared to light microscopy counts. The performance of the lysis method (90.3% efficient), Longmire's buffer, and the DinoDtec qPCR assay (tested across a range of Alexandrium species (90.7-106.9% efficiency; r2 > 0.99)) was found to be specific, sensitive, and efficient. We tested the application of this workflow weekly from May 2016 to 30th October 2017 to compare the relationship between sxtA4 copies L-1 in seawater and PSTs in mussel tissue (Mytilus galloprovincialis) on-farm and spatially (across multiple sites), effectively demonstrating an ∼2 week early warning of two A. catenella HABs (r = 0.95). Our tool provides an early, accurate, and efficient method for the identification of PST risk in shellfish aquaculture.

Development of saxitoxin biosynthesis gene sxtB-targeted qPCR assay for the quantification of toxic dinoflagellates Alexandrium catenella (group I) and A. pacificum (group IV) occurring in the Korean coast

Toxic dinoflagellate Alexandrium can produce saxitoxins (STXs) and cause paralytic shellfish poisoning (PSP), and thus they are monitored for environmental safety management. Microscopic discrimination of dinoflagellates is difficult to distinguish between toxic and non-toxic species due to their similar morphology. Meanwhile, an alternative quantitative PCR (qPCR) assay is sensitive, rapid, and cost-effective for harmful species monitoring. Herein, we developed a novel qPCR assay to detect the STXs biosynthesis gene sxtB of Alexandrium catenella and A. pacificum, the leading cause of PSP outbreaks in Asian coasts and worldwide. The newly designed sxtB TaqMan probes target the species without any positive signal in other relative dinoflagellates. Deming regression analysis revealed that the sxtB copy number of A. catenella and A. pacificum was 3.6 and 4.1 copies per cell, respectively. During the blooming periods (April 13th-14th, 2020), only A. catenella cells were detected through the qPCR assay, ranging from 5.0 × 10 to 2.5 × 104 eq cells L-1. In addition, sxtB qPCR quantified more accurately compared to large subunit (LSU) rRNA targeting qPCR assay that overestimate cell density. Besides, the sensitivity of sxtB was higher compared to the microscope when the species were rarely present (5.0 × 102 cells L-1). These suggest that the sxtB qPCR assay can be applied to toxic Alexandrium monitoring in the Korean coast, even in the early stage of bloomings.

Recovering sedimentary ancient DNA of harmful dinoflagellates accumulated over the last 9000 years off Eastern Tasmania, Australia

Harmful algal blooms (HABs) have had significant adverse impacts on the seafood industry along the Tasmanian east coast over the past 4 decades. To investigate the history of regional HABs, we performed analyses of sedimentary ancient DNA (sedaDNA) in coastal sediments up to ~9000 years old collected inshore and offshore of Maria Island, Tasmania. We used metagenomic shotgun sequencing and a hybridisation capture array ("HABbaits1") to target three harmful dinoflagellate genera, Alexandrium, Gymnodinium, and Noctiluca. Bioinformatic and DNA damage analyses verified the authenticity of the sedaDNA sequences. Our results show that dinoflagellates of Alexandrium genera have been present off eastern Tasmania during the last ~8300 years, and we sporadically detected and unambiguously verified sequences of Gymnodinium catenatum that were present offshore up to ~7600 years ago. We also recovered sedaDNA of the fragile, soft-bodied Noctiluca scintillans with increased relative abundance since 2010, consistent with plankton surveys. This study enabled us to identify challenges of sedaDNA sequence validation (in particular for G. catenatum, a microreticulate gymnodinoid species) and provided guidance for the development of tools to monitor past and present HAB species and improvement of future HAB event predictions.

Response of planktonic microbial assemblages to disturbance in an urban sub-tropical estuary

Microbes are sensitive indicators of estuarine processes because they respond rapidly to dynamic disturbance events. As most of the world's population lives in urban areas and climate change-related disturbance events are becoming more frequent, estuaries bounded by cities are experiencing increasing stressors, at the same time that their ecosystem services are required more than ever. Here, using a multidisciplinary approach, we determined the response of planktonic microbial assemblages in response to seasonality and a rainfall disturbance in an urban estuary bounded by Australia's largest city, Sydney. We used molecular barcoding (16S, 18S V4 rRNA) and microscopy-based identification to compare microbial assemblages at locations with differing characteristics and urbanisation histories. Across 142 samples, we identified 8,496 unique free-living bacterial zOTUs, 8,175 unique particle associated bacterial zOTUs, and 1,920 unique microbial eukaryotic zOTUs. Using microscopy, we identified only the top <10% abundant, larger eukaryotic taxa (>10 µm), however quantification was possible. The site with the greater history of anthropogenic impact showed a more even community of associated bacteria and eukaryotes, and a significant increase in dissolved inorganic nitrogen following rainfall, when compared to the more buffered site. This coincided with a reduced proportional abundance of Actinomarina and Synechococcus spp., a change in SAR 11 clades, and an increase in the eukaryotic microbial groups Dinophyceae, Mediophyceae and Bathyoccocaceae, including a temporary dominance of the harmful algal bloom dinoflagellate Prorocentrum cordatum (syn. P. minimum). Finally, a validated hydrodynamic model of the estuary supported these results, showing that the more highly urbanised and upstream location consistently experienced a higher magnitude of salinity reduction in response to rainfall events during the study period. The best abiotic variables to explain community dissimilarities between locations were TDP, PN, modelled temperature and salinity (r = 0.73) for the free living bacteria, TP for the associated bacteria (r = 0.43), and modelled temperature (r = 0.28) for the microbial eukaryotic communities. Overall, these results show that a minor disturbance such as a brief rainfall event can significantly shift the microbial assemblage of an anthropogenically impacted area within an urban estuary to a greater degree than a seasonal change, but may result in a lesser response to the same disturbance at a buffered, more oceanic influenced location. Fine scale research into the factors driving the response of microbial communities in urban estuaries to climate related disturbances will be necessary to understand and implement changes to maintain future estuarine ecosystem services.

Genomic copy number variability at the genus, species and population levels impacts in situ ecological analyses of dinoflagellates and harmful algal blooms

The application of meta-barcoding, qPCR, and metagenomics to aquatic eukaryotic microbial communities requires knowledge of genomic copy number variability (CNV). CNV may be particularly relevant to functional genes, impacting dosage and expression, yet little is known of the scale and role of CNV in microbial eukaryotes. Here, we quantify CNV of rRNA and a gene involved in Paralytic Shellfish Toxin (PST) synthesis (sxtA4), in 51 strains of 4 Alexandrium (Dinophyceae) species. Genomes varied up to threefold within species and ~7-fold amongst species, with the largest (A. pacificum, 130 ± 1.3 pg cell-1 /~127 Gbp) in the largest size category of any eukaryote. Genomic copy numbers (GCN) of rRNA varied by 6 orders of magnitude amongst Alexandrium (102- 108 copies cell-1) and were significantly related to genome size. Within the population CNV of rRNA was 2 orders of magnitude (105 - 107 cell-1) in 15 isolates from one population, demonstrating that quantitative data based on rRNA genes needs considerable caution in interpretation, even if validated against locally isolated strains. Despite up to 30 years in laboratory culture, rRNA CNV and genome size variability were not correlated with time in culture. Cell volume was only weakly associated with rRNA GCN (20-22% variance explained across dinoflagellates, 4% in Gonyaulacales). GCN of sxtA4 varied from 0-102 copies cell-1, was significantly related to PSTs (ng cell-1), displaying a gene dosage effect modulating PST production. Our data indicate that in dinoflagellates, a major marine eukaryotic group, low-copy functional genes are more reliable and informative targets for quantification of ecological processes than unstable rRNA genes.

A fast and accurate method for specific detection and quantification of the bloom-forming microalgae Karlodinium veneficum in the marine environment

Karlodinium veneficum is a toxic benthic globally distributed dinoflagellate which has direct impacts on human health and the environment. Early and accurate detection of this harmful algal bloom-forming species could be useful for potential risks monitoring and management. In the present work, a real-time PCR targeting the internal transcribed spacer ribosomal DNA region for the specific detection and absolute quantification of K. veneficum was designed. Then, the assay conditions were adjusted and validated. The developed qPCR was highly specific for the target species and displayed no cross-reactivity with closely related dinoflagellates and/or other microalgal species commonly distributed along the Tunisian coast. Its lowest detection limit was 5 rDNA copies per reaction, which is often considered satisfying. qPCR assay enumeration accuracy was evaluated using artificially inoculated environmental samples. The comparison of the cell abundance estimates obtained by qPCR assay with the theoretical estimates showed no statistically significant difference across a range of concentrations. We suggest that the qPCR approach developed in the present study may be a valuable tool to investigate the distribution and seasonal dynamics of K. veneficum in marine environments.

Apparent biogeographical trends in Alexandrium blooms for northern Europe: identifying links to climate change and effective adaptive actions

The marine dinoflagellate Alexandrium Halim represents perhaps the most significant and intensively studied genus with respect to species diversity, life history strategies, toxigenicity, biogeographical distribution, and global magnitude and consequences harmful algal blooms (HABs). The socioeconomic impacts, environmental and human health risks, and mitigation strategies for toxigenic Alexandrium blooms have also been explored in recent years. Human adaptive actions based on future scenarios of bloom dynamics and shifts in biogeographical distribution under climate-change parameters remain under development and not yet implemented on a regional scale. In the CoCliME (Co-development of climate services for adaptation to changing marine ecosystems) project these issues were addressed with respect to past, current and anticipated future status of key HAB genera and expected benefits of enhanced monitoring. Data on the distribution and frequency of Alexandrium blooms related to paralytic shellfish toxin (PST) events from key CoCliME Case Study areas, comprising the North Sea and adjacent Kattegat-Skagerrak, Norwegian Sea, and Baltic Sea, and eastern North Atlantic marginal seas, were evaluated in a contemporary and historical context over the past several decades. The first evidence of possible biogeographical expansion of Alexandrium taxa into eastern Arctic gateways was provided from DNA barcoding signatures. Various key climate change indicators, such as salinity, temperature, and water-column stratification, relevant to Alexandrium bloom initiation and development were identified. The possible influence of changing variables on bloom dynamics, magnitude, frequency and spatial and temporal distribution were interpreted in the context of regional ocean climate models. These climate change impact indicators may play key roles in selecting for the occurrence and diversity of Alexandrium species within the broader microeukaryote communities. For example, shifts to higher temperature and lower salinity regimes predicted for the southern North Sea indicate the potential for increased Alexandrium blooms, currently absent from this area. Ecological and socioeconomic impacts of Alexandrium blooms and effects on fisheries and aquaculture resources and coastal ecosystem function are evaluated, and, where feasible, effective adaptation strategies are proposed herein as emerging climate services.

In vitro selection and characterization of a DNA aptamer targeted to Prorocentrum minimum-A common harmful algae

Prorocentrum minimum is a common diarrhetic shellfish toxins-producing marine microalga that may seriously endanger marine resources and cause great economic losses. The development of a novel rapid detection technique is of great importance for the prevention and control of the damage caused by P. minimum. In this study, the aptamer against P. minimum was for the first time generated from an artificially synthesized single-stranded DNA library by systematic evolution of ligand by exponential enrichment (SELEX), using P. minimum and P. minimum-related species, including Prorocentrum donghaiense, Prorocentrum lima and Prorocentrum micans as target and counter-screening species, respectively. The aptamer library was successfully obtained at the end of 18 rounds of SELEX-screening by continuously monitoring the binding ratio of the resultant ssDNA from each round. Three sequences (Apt 1, Apt 2 and Apt 3) with the highest frequency in the aptamer library resulted from high-throughput sequencing were first selected as candidate aptamers. The secondary structure of these sequences was predicted and analyzed. In addition, the specificity and affinity of these candidate aptamers were determined by flow cytometry analysis. The results indicated that these aptamers had high specificity and affinity, with a K_D of (224.6 ± 8.8) nM (Apt 1), (286.6 ± 13.9) nM (Apt 2) and (388.5 ± 44.6) nM (Apt 3), respectively. Apt 1 was therefore chosen as the best aptamer against P. minimum. Finally, the fluorescence microscopic examination further confirmed that Apt 1 can well bind to P. minimum. In summary, Apt 1 may be promising for being used as a novel molecular recognition element for P. minimum.

Mapping the development of a Dinophysis bloom in a shellfish aquaculture area using a novel molecular qPCR assay

Diarrhetic shellfish toxins produced by certain species of the marine dinoflagellate Dinophysis can accumulate in shellfish in high concentrations, representing a significant food safety issue worldwide. This risk is routinely managed by monitoring programs in shellfish producing areas, however the methods used to detect these harmful marine microbes are not usually automated nor conducted onsite, and are often expensive and require specialized expertise. Here we designed a quantitative real-time polymerase chain reaction (qPCR) assay based on the ITS-5.8S ribosomal region of Dinophysis spp. and evaluated its specificity, efficiency, and sensitivity to detect species belonging to this genus. We designed and tested twenty sets of primers pairs using three species of Dinophysis - D. caudata, D. fortii and D. acuminata. We optimized a qPCR assay using the primer pair that sufficiently amplified each of the target species (Dacu_11F/Dacu_11R), and tested this assay for cross-reactivity with other dinoflagellates and diatoms in the laboratory (11 species) and in silico 8 species (15 strains) of Dinophysis, 3 species of Ornithocercus and 2 species of Phalacroma. The qPCR assay returned efficiencies of 92.4% for D. caudata, 91.3% for D fortii, and 91.5% for D. acuminata, while showing no cross-reactivity with other phytoplankton taxa. Finally, we applied this assay to a D. acuminata bloom which occurred in an oyster producing estuary in south eastern Australia, and compared cell numbers inferred by qPCR to those determined by microscopy counts (max abund. ∼6.3 × 10³ and 5.3 × 10³ cells L^-1 respectively). Novel molecular tools such as qPCR have the potential to be used on-farm, be automated, and provide an early warning for the management of harmful algal blooms.

Plankton biodiversity and species co-occurrence based on environmental DNA – a multiple marker study

Metabarcoding in combination with high-throughput sequencing (HTS) allows simultaneous detection of multiple taxa by targeting single or several taxonomically informative gene regions from environmental DNA samples. In this study, a multiple-marker HTS approach was applied to investigate the plankton diversity and seasonal succession in the Baltic Sea from winter to autumn. Four different markers targeting the 16S, 18S, and 28S ribosomal RNA genes were employed, including a marker for more efficient dinoflagellate detection. Typical seasonal changes were observed in phyto- and bacterioplankton communities. In phytoplankton, the appearance patterns of selected common, dominant, or harmful species followed the patterns also confirmed based on 20 years of phytoplankton monitoring data. In the case of zooplankton, both macro- and microzooplankton species were detected. However, no seasonal patterns were detected in their appearance. In total, 15 and 2 new zoo- and phytoplankton species were detected from the Baltic Sea. HTS approach was especially useful for detecting microzooplankton species as well as for investigating the co-occurrence and potential interactions of different taxa. The results of this study further exemplify the efficiency of metabarcoding for biodiversity monitoring and the advantage of employing multiple markers through the detection of species not identifiable based on a single marker survey and/or by traditional morphology-based methods.

Using qPCR and high-resolution sensor data to model a multi-species Pseudo-nitzschia (Bacillariophyceae) bloom in southeastern Australia

Harmful algal blooms, including those caused by the toxic diatom Pseudo-nitzschia, can have significant impacts on human health, ecosystem functioning and ultimately food security. In the current study we characterized a bloom of species of Pseudo-nitzschia that occurred in a south-eastern Australian oyster-growing estuary in 2019. Using light microscopy, combined with molecular (ITS/5.8S and LSU D1-D3 rDNA regions) and toxicological evidence, we observed the bloom to consist of multiple species of Pseudo-nitzschia including P. cf. cuspidata, P. hasleana, P. fraudulenta and P. multiseries, with P. cf. cuspidata being the only species that produced domoic acid (3.1 pg DA per cell). As several species of Pseudo-nitzschia co-occurred, only one of which produced DA, we developed a rapid, sensitive and efficient quantitative real-time polymerase chain reaction (qPCR) assay to detect only species belonging to the P. pseudodelicatissima complex Clade I, to which P. cf. cuspidata belongs, and this indicated that P. cuspidata or closely related strains may have dominated the Pseudo-nitzschia community at this time. Finally, using high resolution water temperature and salinity sensor data, we modeled the relationship between light microscopy determined abundance of P. delicatissima group and environmental variables (temperature, salinity, rainfall) at two sites within the estuary. A total of eight General Linear Models (GLMs) explaining between 9 and 54% of the deviance suggested that the temperature (increasing) and/or salinity (decreasing) data were generally more predictive of high cell concentrations than the rainfall data at both sites, and that overall, cell concentrations were more predictive at the more oceanic site than the more upstream site, using this method. We conclude that the combination of rapid molecular methods such as qPCR and real-time sensor data modeling, can provide a more rapid and effective early warning of harmful algal blooms of species of Pseudo-nitzschia, resulting in more beneficial regulatory and management outcomes.

Low Temperature and Cold Stress Significantly Increase Saxitoxins (STXs) and Expression of STX Biosynthesis Genes sxtA4 and sxtG in the Dinoflagellate Alexandrium catenella

Toxic dinoflagellate Alexandrium spp. produce saxitoxins (STXs), whose biosynthesis pathway is affected by temperature. However, the link between the regulation of the relevant genes and STXs' accumulation and temperature is insufficiently understood. In the present study, we evaluated the effects of temperature on cellular STXs and the expression of two core STX biosynthesis genes (sxtA4 and sxtG) in the toxic dinoflagellate Alexandrium catenella Alex03 isolated from Korean waters. We analyzed the growth rate, toxin profiles, and gene responses in cells exposed to different temperatures, including long-term adaptation (12, 16, and 20 °C) and cold and heat stresses. Temperature significantly affected the growth of A. catenella, with optimal growth (0.49 division/day) at 16 °C and the largest cell size (30.5 µm) at 12 °C. High concentration of STXs eq were detected in cells cultured at 16 °C (86.3 fmol/cell) and exposed to cold stress at 20→12 °C (96.6 fmol/cell) compared to those at 20 °C and exposed to heat stress. Quantitative real-time PCR (qRT-PCR) revealed significant gene expression changes of sxtA4 in cells cultured at 16 °C (1.8-fold) and cold shock at 20→16 °C (9.9-fold). In addition, sxtG was significantly induced in cells exposed to cold shocks (20→16 °C; 19.5-fold) and heat stress (12→20 °C; 25.6-fold). Principal component analysis (PCA) revealed that low temperature (12 and 16 °C) and cold stress were positively related with STXs' production and gene expression levels. These results suggest that temperature may affect the toxicity and regulation of STX biosynthesis genes in dinoflagellates.

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.

HPLC-UV, Metabarcoding and Genome Skims of Botanical Dietary Supplements: A Case Study in Echinacea

The use of DNA-based methods to authenticate botanical dietary supplements has been vigorously debated for a variety of reasons. More comparisons of DNA-based and chemical methods are needed, and concordant evaluation of orthogonal approaches on the same products will provide data to better understand the strengths and weaknesses of both approaches. The overall application of DNA-based methods is already firmly integrated into a wide array of continually modernizing stand alone and complementary authentication protocols. Recently, the use of full-length chloroplast genome sequences provided enhanced discriminatory capacity for closely related species of Echinacea compared to traditional DNA barcoding approaches (matK and rbcL). Here, two next-generation sequencing approaches were used: (1) genome skimming and (2) PCR amplicon (metabarcoding). The two genetic approaches were then combined with HPLC-UV to evaluate 20 commercially available dietary supplements of Echinacea representing "finished" products. The trade-offs involved in different DNA approaches were discussed, with a focus on how DNA methods support existing, accepted chemical methods. In most of the products (19/20), HPLC-UV suggested the presence of Echinacea spp. While metabarcoding was not useful with this genus and instead only resolved 7 products to the family level, genome skimming was able to resolve to species (9) or genus (1) with the 10/20 products where it was successful. Additional ingredients that HPLC-UV was unable to identify were also found in four products along with the relative sequence proportion of the constituents. Additionally, genome skimming was able to identify one product that was a different Echinacea species entirely.

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.

Assessing the Use of Molecular Barcoding and qPCR for Investigating the Ecology of Prorocentrum minimum (Dinophyceae), a Harmful Algal Species

Prorocentrum minimum is a species of marine dinoflagellate that occurs worldwide and can be responsible for harmful algal blooms (HABs). Some studies have reported it to produce tetrodotoxin; however, results have been inconsistent. qPCR and molecular barcoding (amplicon sequencing) using high-throughput sequencing have been increasingly applied to quantify HAB species for ecological analyses and monitoring. Here, we isolated a strain of P. minimum from eastern Australian waters, where it commonly occurs, and developed and validated a qPCR assay for this species based on a region of ITS rRNA in relation to abundance estimates from the cultured strain as determined using light microscopy. We used this tool to quantify and examine ecological drivers of P. minimum in Botany Bay, an estuary in southeast Australia, for over ~14 months in 2016–2017. We compared abundance estimates using qPCR with those obtained using molecular barcoding based on an 18S rRNA amplicon. There was a significant correlation between the abundance estimates from amplicon sequencing and qPCR, but the estimates from light microscopy were not significantly correlated, likely due to the counting method applied. Using amplicon sequencing, ~600 unique actual sequence variants (ASVs) were found, much larger than the known phytoplankton diversity from this region. P. minimum abundance in Botany Bay was found to be significantly associated with lower salinities and higher dissolved CO₂ levels.

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.

Harmful algal blooms and their effects in coastal seas of Northern Europe

Harmful algal blooms (HAB) are recurrent phenomena in northern Europe along the coasts of the Baltic Sea, Kattegat-Skagerrak, eastern North Sea, Norwegian Sea and the Barents Sea. These HABs have caused occasional massive losses for the aquaculture industry and have chronically affected socioeconomic interests in several ways. This status review gives an overview of historical HAB events and summarises reports to the Harmful Algae Event Database from 1986 to the end of year 2019 and observations made in long term monitoring programmes of potentially harmful phytoplankton and of phycotoxins in bivalve shellfish. Major HAB taxa causing fish mortalities in the region include blooms of the prymnesiophyte Chrysochromulina leadbeateri in northern Norway in 1991 and 2019, resulting in huge economic losses for fish farmers. A bloom of the prymesiophyte Prymnesium polylepis (syn. Chrysochromulina polylepis) in the Kattegat-Skagerrak in 1988 was ecosystem disruptive. Blooms of the prymnesiophyte Phaeocystis spp. have caused accumulations of foam on beaches in the southwestern North Sea and Wadden Sea coasts and shellfish mortality has been linked to their occurrence. Mortality of shellfish linked to HAB events has been observed in estuarine waters associated with influx of water from the southern North Sea. The first bloom of the dictyochophyte genus Pseudochattonella was observed in 1998, and since then such blooms have been observed in high cell densities in spring causing fish mortalities some years. Dinoflagellates, primarily Dinophysis spp., intermittently yield concentrations of Diarrhetic Shellfish Toxins (DST) in blue mussels, Mytilus edulis, above regulatory limits along the coasts of Norway, Denmark and the Swedish west coast. On average, DST levels in shellfish have decreased along the Swedish and Norwegian Skagerrak coasts since approximately 2006, coinciding with a decrease in the cell abundance of D. acuta. Among dinoflagellates, Alexandrium species are the major source of Paralytic Shellfish Toxins (PST) in the region. PST concentrations above regulatory levels were rare in the Skagerrak-Kattegat during the three decadal review period, but frequent and often abundant findings of Alexandrium resting cysts in surface sediments indicate a high potential risk for blooms. PST levels often above regulatory limits along the west coast of Norway are associated with A. catenella (ribotype Group 1) as the main toxin producer. Other Alexandrium species, such as A. ostenfeldii and A. minutum, are capable of producing PST among some populations but are usually not associated with PSP events in the region. The cell abundance of A. pseudogonyaulax, a producer of the ichthyotoxin goniodomin (GD), has increased in the Skagerrak-Kattegat since 2010, and may constitute an emerging threat. The dinoflagellate Azadinium spp. have been unequivocally linked to the presence of azaspiracid toxins (AZT) responsible for Azaspiracid Shellfish Poisoning (AZP) in northern Europe. These toxins were detected in bivalve shellfish at concentrations above regulatory limits for the first time in Norway in blue mussels in 2005 and in Sweden in blue mussels and oysters (Ostrea edulis and Crassostrea gigas) in 2018. Certain members of the diatom genus Pseudo-nitzschia produce the neurotoxin domoic acid and analogs known as Amnesic Shellfish Toxins (AST). Blooms of Pseudo-nitzschia were common in the North Sea and the Skagerrak-Kattegat, but levels of AST in bivalve shellfish were rarely above regulatory limits during the review period. Summer cyanobacteria blooms in the Baltic Sea are a concern mainly for tourism by causing massive fouling of bathing water and beaches. Some of the cyanobacteria produce toxins, e.g. Nodularia spumigena, producer of nodularin, which may be a human health problem and cause occasional dog mortalities. Coastal and shelf sea regions in northern Europe provide a key supply of seafood, socioeconomic well-being and ecosystem services. Increasing anthropogenic influence and climate change create environmental stressors causing shifts in the biogeography and intensity of HABs. Continued monitoring of HAB and phycotoxins and the operation of historical databases such as HAEDAT provide not only an ongoing status report but also provide a way to interpret causes and mechanisms of HABs.

Quantitative PCR assay for the simultaneous identification and enumeration of multiple Karenia species

Quantitative PCR (qPCR) is the method of choice for specific detection and quantification of harmful algal bloom (HAB) species. Development of qPCR assay for simultaneous enumeration of species that frequently co-exist in HABs is required. A high sensitivity TaqMan qPCR assay, using probe and primers, located at ITS1-5.8S-ITS2 rDNA region, detecting, specifically, Karenia selliformis, K. bidigitata, and K. mikimotoi, was designed. ITS1-5.8S-ITS2 rDNA region copy numbers per Karenia cell genome were estimated to 217.697 ± 67.904, allowing cell quantification. An application of the designed methodology in field samples has been conducted, and it showed high sensitivity (detection of around 10^-1 cell/100 mg of bivalve mollusk tissue, equivalent to about 20 copies of the target sequence). We suggest that the optimized method could contribute to early detection of three closely related Karenia species in seafood cultivating areas to promote control quality, guarantee a fast and effective intervention, and improve public health prevention.

Sierra Nevada sweep: metagenomic measurements of bioaerosols vertically distributed across the troposphere

To explore how airborne microbial patterns change with height above the Earth’s surface, we flew NASA’s C-20A aircraft on two consecutive days in June 2018 along identical flight paths over the US Sierra Nevada mountain range at four different altitudes ranging from 10,000 ft to 40,000 ft. Bioaerosols were analyzed by metagenomic DNA sequencing and traditional culturing methods to characterize the composition and diversity of atmospheric samples compared to experimental controls. The relative abundance of taxa changed significantly at each altitude sampled, and the diversity profile shifted across the two sampling days, revealing a regional atmospheric microbiome that is dynamically changing. The most proportionally abundant microbial genera were Mycobacterium and Achromobacter at 10,000 ft; Stenotrophomonas and Achromobacter at 20,000 ft; Delftia and Pseudoperonospora at 30,000 ft; and Alcaligenes and Penicillium at 40,000 ft. Culture-based detections also identified viable Bacillus zhangzhouensis, Bacillus pumilus, and Bacillus spp. in the upper troposphere. To estimate bioaerosol dispersal, we developed a human exposure likelihood model (7-day forecast) using general aerosol characteristics and measured meteorological conditions. By coupling metagenomics to a predictive atmospheric model, we aim to set the stage for field campaigns that monitor global bioaerosol emissions and impacts.

Distribution and abundance of azaspiracid-producing dinophyte species and their toxins in North Atlantic and North Sea waters in summer 2018

Representatives of the marine dinophyte family Amphidomataceae produce lipophilic phycotoxins called azaspiracids (AZA) which may cause azaspiracid shellfish poisoning (AZP) in humans after consumption of contaminated seafood. Three of the four known toxigenic species are observed frequently in the eastern North Atlantic. In 2018, a research survey was performed to strengthen knowledge on the distribution and abundance of toxigenic Amphidomataceae and their respective toxins in Irish coastal waters and in the North Sea. Species-specific quantification of the three toxigenic species (Azadinium spinosum, Azadinium poporum and Amphidoma languida) was based on recently developed qPCR assays, whose performance was successfully validated and tested with specificity tests and spike experiments. The multi-method approach of on-board live microscopy, qPCR assays and chemical AZA-analysis revealed the presence of Amphidomataceae in the North Atlantic including the three targeted toxigenic species and their respective AZA analogues (AZA-1, -2, -33, -38, -39). Azadinium spinosum was detected at the majority of Irish stations with a peak density of 8.3 x 10⁴ cells L^-1 and AZA (AZA-1, -2, -33) abundances up to 1,274 pg L^-1. Amphidoma languida was also present at most Irish stations but appeared in highest abundance in a bloom at a central North Sea station with a density of 1.2 x 10⁵ cells L^-1 and an AZA (AZA-38, -39) abundances of 618 pg L^-1. Azadinium poporum was detected sporadically at the Irish south coast and North Sea and was rather low in abundance during this study. The results confirmed the wide distribution and frequent occurrence of the target species in the North Atlantic area and revealed, for the first time, bloom abundances of toxigenic Amphidomataceae in this area. This emphasizes the importance of future studies and monitoring of amphidomatacean species and their respective AZA analogues in the North Atlantic.

First Detection of Paralytic Shellfish Toxins from Alexandrium pacificum above the Regulatory Limit in Blue Mussels (Mytilus galloprovincialis) in New South Wales, Australia

In 2016, 2017 and 2018, elevated levels of the species Alexandrium pacificum were detected within a blue mussel (Mytilus galloprovincialis) aquaculture area at Twofold Bay on the south coast of New South Wales, Australia. In 2016, the bloom persisted for at least eight weeks and maximum cell concentrations of 89,000 cells L-1 of A. pacificum were reported. The identity of A. pacificum was confirmed using molecular genetic tools (qPCR and amplicon sequencing) and complemented by light and scanning electron microscopy of cultured strains. Maximum reported concentrations of paralytic shellfish toxins (PSTs) in mussel tissue was 7.2 mg/kg PST STX equivalent. Elevated cell concentrations of A. pacificum were reported along the adjacent coastal shelf areas, and positive PST results were reported from nearby oyster producing estuaries during 2016. This is the first record of PSTs above the regulatory limit (0.8 mg/kg) in commercial aquaculture in New South Wales since the establishment of routine biotoxin monitoring in 2005. The intensity and duration of the 2016 A. pacificum bloom were unusual given the relatively low abundances of A. pacificum in estuarine and coastal waters of the region found in the prior 10 years.

New insights into the harmful algae inhibition by Spartina alterniflora: Cellular physiology and metabolism of extracellular secretion

This study investigated the response of three toxic harmful algal blooms (HABs) species, Heterosigma akashiwo, Chattonella marina, and Alexandrium tamarense to coastal invasive plant, Spartina alterniflora. In this study, the growth of three HABs species were suppressed significantly by S. alterniflora extracts, showing a dose-response relationship. The raphidophyte H. akashiwo and C. marina exhibited higher inhibitory response with EC_50,7d decreased by 14% and 75% as compared to the dinoflagellate A. tamarense. C. marina was the most sensitive among the three species. S. alterniflora extracts disrupted algal cellular integrity and photosynthesis. Furthermore, the extracellular organic matters were detected by fluorescence excitation-emission matrix. Algal metabolites, protein-like substances (tyrosine-like peak and tryptophan-like peak) decreased as time prolonged and the humic-like substances (UVA marine humic-like peak) increased when algal cells were exposed to S. alterniflora extracts. These results provide new insights to the inhibition mechanism of S. alterniflora extracts on HABs species.

Monitoring of the toxic dinoflagellate Alexandrium catenella in Osaka Bay, Japan using a massively parallel sequencing (MPS)-based technique

Since 2002, blooms of Alexandrium catenella sensu Fraga et al. (2015) and paralytic shellfish toxicity events have occurred almost yearly in Osaka Bay, Japan. To better understand the triggers for reoccurring A. catenella blooms in Osaka Bay, phytoplankton community was monitored during the spring seasons of 2012-2015. Monitoring was performed using massively parallel sequencing (MPS)-based technique on amplicon sequences of the 18S rRNA gene. Dense blooms of A. catenella occurred every year except in 2012, however, there was no significant correlation with the environmental parameters investigated. Plankton community diversity decreased before and middle of the A. catenella blooms, suggesting that the decline in diversity could be an indicator for the bloom occurrence. The yearly abundance pattern of A. catenella cells obtained by morphology-based counting coincided with the relative sequence abundances, which supports the effectiveness of MPS-based phytoplankton monitoring.

The Genetic Basis of Toxin Biosynthesis in Dinoflagellates

In marine ecosystems, dinoflagellates can become highly abundant and even dominant at times, despite their comparatively slow growth rates. One factor that may play a role in their ecological success is the production of complex secondary metabolite compounds that can have anti-predator, allelopathic, or other toxic effects on marine organisms, and also cause seafood poisoning in humans. Our knowledge about the genes involved in toxin biosynthesis in dinoflagellates is currently limited due to the complex genomic features of these organisms. Most recently, the sequencing of dinoflagellate transcriptomes has provided us with valuable insights into the biosynthesis of polyketide and alkaloid-based toxin molecules in dinoflagellate species. This review synthesizes the recent progress that has been made in understanding the evolution, biosynthetic pathways, and gene regulation in dinoflagellates with the aid of transcriptomic and other molecular genetic tools, and provides a pathway for future studies of dinoflagellates in this exciting omics era.

Environmental drivers of unprecedented Alexandrium catenella dinoflagellate blooms off eastern Tasmania, 2012-2018

Blooms of the highly toxic dinoflagellate Alexandrium catenella (previously referred to as tamarense group 1) were first detected off eastern Tasmania in 2012 and have since been responsible for incidences of human paralytic shellfish poisoning and extended closures (up to 25 weeks) of mussel, oyster, scallop, abalone and rock lobster industries (up to 150 mg/kg PST in mussels). Investigation of meteorological and oceanographic influences indicate that the annually recurrent winter-spring blooms (June-Oct) occur within a narrow water temperature window (10-15 °C) under two distinct sets of conditions: (1) following high rainfall and land run-off, under relatively light winds; and (2) following periods of anomalously low air temperatures and associated cooling of shallow coastal waters, again under relatively light winds. The common driver of blooms appears to be the development of stratification in coastal waters, via salinity and/or temperature gradients. We propose a framework for evaluating the risk of Alexandrium with the aim of developing a forecasting capability, and compare these environmental conditions with historic data to understand the recent advent of these blooms.