Sudden peak in tetrodotoxin in French oysters during the summer of 2021: Source investigation using microscopy, metabarcoding and droplet digital PCR

Tetrodotoxin (TTX) is a potent neurotoxin causing human intoxications from contaminated seafood worldwide and is of emerging concern in Europe. Shellfish have been shown to contain varying TTX concentrations globally, with concentrations typically higher in Pacific oysters Crassostrea gigas in Europe. Despite many decades of research, the source of TTX remains unknown, with bacterial or algal origins having been suggested. The aim of this study was to identify potential source organisms causing TTX contamination in Pacific oysters in French coastal waters, using three different techniques. Oysters were deployed in cages from April to September 2021 in an estuary where TTX was previously detected. Microscopic analyses of water samples were used to investigate potential microalgal blooms present prior or during the peak in TTX. Differences in the bacterial communities from oyster digestive glands (DG) and remaining flesh were explored using metabarcoding, and lastly, droplet digital PCR assays were developed to investigate the presence of Cephalothrix sp., one European TTX-bearing species in the DG of toxic C. gigas. Oysters analysed by liquid chromatography-tandem mass spectrometry contained quantifiable levels of TTX over a three-week period (24 June-15 July 2021), with concentrations decreasing in the DG from 424 μg/kg for the first detection to 101 μg/kg (equivalent to 74 to 17 μg/kg of total flesh), and trace levels being detected until August 13, 2021. These concentrations are the first report of the European TTX guidance levels being exceeded in French shellfish. Microscopy revealed that some microalgae bloomed during the TTX peak, (e.g., Chaetoceros spp., reaching 40,000 cells/L). Prokaryotic metabarcoding showed increases in abundance of Rubritaleaceae (genus Persicirhabdus) and Neolyngbya, before and during the TTX peak. Both phyla have previously been described as possible TTX-producers and should be investigated further. Droplet digital PCR analyses were negative for the targeted TTX-bearing genus Cephalothrix.

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

Fatty acid production and associated gene pathways are altered by increased salinity and dimethyl sulfoxide treatments during cryopreservation of Symbiodinium pilosum (Symbiodiniaceae)

The Symbiodinium genus is ancestral among other Symbiodiniaceae lineages with species that are both symbiotic and free living. Changes in marine ecosystems threaten their existence and crucial ecological roles. Cryopreservation offers an avenue for their long-term storage for future habitat restoration after coral bleaching. In our previous study we demonstrated that high salinity treatments of Symbiodiniaceae isolates led to changes in their fatty acid (FA) profiles and higher cell viabilities after cryopreservation. In this study, we investigated the role of increased salinity on FA production and the genes involved in FA biosynthesis and degradation pathways during the cryopreservation of Symbiodinium pilosum. Overall, there was a twofold increase in mass of FAs produced by S. pilosum after being cultured in medium with increased salinity (54 parts per thousand; ppt). Dimethyl sulfoxide (Me2SO) led to a ninefold increase of FAs in standard salinity (SS) treatment, compared to a fivefold increase in increased salinity (IS) treatments. The mass of the FA classes returned to baseline during recovery. Transcriptomic analyses showed an acyl carrier protein gene was significantly upregulated after Me2SO treatment in the SS cultures. Cytochrome P450 reductase genes were significantly down regulated after Me2SO addition in SS treatment preventing FA degradation. These changes in the expression of FA biosynthesis and degradation genes contributed to more FAs in SS treated isolates. Understanding how increased salinity changes FA production and the roles of specific genes in regulating FA pathways will help improve current freezing protocols for Symbiodiniaceae and other marine microalgae.

Differential responses of selectively bred mussels (Perna canaliculus) to heat stress-survival, immunology, gene expression and microbiome diversity

New Zealand's green-lipped mussel (Perna canaliculus) is an ecologically and economically important species. Marine heatwaves are increasing in frequency around NZ's coastline, and these events are correlated with increased stress and mortality of some aquaculture species. This study aimed to identify general biomarkers of heat stress in P. canaliculus and to assess whether responses differed between genetically distinct selectively bred mussels. We exposed three families of selectively bred mussels (families A, B and C) to three seawater temperature regimes in the laboratory: 1) a "control" treatment (ambient 12°C), 2) a 26°C heat challenge with a subsequent recovery period, and 3) a sustained 26°C heat challenge with no recovery. We investigated whether the survival, immune response (hemocyte concentration and viability, oxidative stress and total antioxidant capacity), hemocyte gene expression and gill microbiome differed between the families during the temperature challenges. In the sustained heat-stress treatment, family A had the highest survival rate (42% compared with 25% and 5% for families C and B, respectively). Gene expression levels significantly shifted during thermal stress and differed between families, with family A more dissimilar than families B and C. Family C had substantially more genes impacted by temperature treatment and timepoint than the other families, while family B had very little genes/pathways that responded to thermal stress. Genes related to heat shock proteins and immune responses (e.g., AIF1, CTSC, TOLL8, CASP9, FNTA, AHCY, CRYAB, PPIF) were upregulated in all families during heat stress. Microbiome species-richness differed between families before and during heat-stress, with family A having a distinctly different microbiome flora than the other families. Microbial diversity changed similarly in all families exposed to prolonged heat-stress, with species of Vibrio and Campylobacter increasing in these mussels. Our study highlights the use of non-lethal sampling of hemocytes as a diagnostic tool to explore the immune response and gene expression of selectively bred mussels, to predict their response to ocean warming. This approach can identify potential thermotolerant candidates for further selective breeding, which may increase the resilience of the mussel aquaculture industry in a warming ocean.

Diversity and distribution of benthic dinoflagellates in Tonga include the potentially harmful genera Gambierdiscus and Fukuyoa

Benthic dinoflagellates that can cause illness, such as ciguatera poisoning (CP), are prevalent around the Pacific but are poorly described in many locations. This study represents the first ecological assessment of benthic harmful algae species in the Kingdom of Tonga, a country where CP occurs regularly. Surveys were conducted in June 2016 in the Tongatapu island group, and in June 2017 across three island groups: Ha'apai, Vava'u, and Tongatapu. Shallow subtidal coastal habitats were investigated by measuring water quality parameters and conducting quadrat surveys. Microalgae samples were collected using either macrophyte collection or the artificial substrate method. Benthic dinoflagellates (Gambierdiscus and/or Fukuyoa, Ostreopsis, and Prorocentrum) were counted using light microscopy, followed by molecular analyses (real-time PCR in 2016 and high throughput sequencing (metabarcoding) in 2017) to identify Gambierdiscus and Fukuyoa to species level. Six species were detected from the Tongatapu island group in 2016 (G. australes, G. carpenteri, G. honu, G. pacificus, F. paulensis, and F. ruetzleri) using real-time PCR. Using the metabarcoding approach in 2017, a total of eight species (G. australes, G. carpenteri, G. honu, G. pacificus, G. cheloniae, G. lewisii, G. polynesiensis, and F. yasumotoi) were detected. Species were detected in mixed assemblages of up to six species, with G. pacificus and G. carpenteri being the most frequently observed. Ha'apai had the highest diversity with eight species detected, which identifies this area as a Gambierdiscus diversity 'hotspot'. Vava'u and Tongatapu had three and six species found respectively. Gambierdiscus polynesiensis, a described ciguatoxin producer and proposed causative agent of CP was found only in Ha'apai and Vava'u in 2017, but not in Tongatapu in either year. Ostreopsis spp. and Prorocentrum spp. were also frequently observed, with Prorocentrum most abundant at the majority of sites. In 2016, the highest number of Gambierdiscus and/or Fukuyoa cells were observed on seagrass (Halodule uninervis) from Sopu, Tongatapu. In 2017, the highest numbers of Gambierdiscus and/or Fukuyoa from artificial substrate samples were recorded in the Halimeda dominant habitat at Neiafu Tahi, Vava'u, a low energy site. This raised the question of the effect of wave motion or currents on abundance measurements from artificial substrates. Differences in detection were noticed between macrophytes and artificial substrates, with higher numbers of species found on artificial substrates. This study provides a baseline of benthic dinoflagellate distributions and diversity for Tonga that may be used for future studies and the development of monitoring programmes.

Sub-tropical benthic/epiphytic dinoflagellates of Aotearoa New Zealand and Rangitāhua Kermadec Islands

Temperatures and temperature anomalies have been increasing in the sub-tropical regions of Aotearoa New Zealand and these changes may impact on harmful algal bloom (HAB) events. Benthic and epiphytic dinoflagellates, particularly the toxin producers, are the focus of this study as it is predicted that under future climate conditions they may produce more toxins or marine animals may become more susceptible to them. The results of past expeditions to Rangitāhua Kermadec Islands and sampling trips to Northland, Aotearoa New Zealand, are summarised and the results of the most recent trips to both regions are presented. The macroalgal habitats of the dinoflagellates are also characterised. Dinoflagellate species not previously identified in Rangitāhua include Coolia canariensis, C. palmyrensis, and C. tropicalis, all identified by DNA sequencing of the large subunit ribosomal RNA region. Gambierdiscus polynesiensis was again isolated and produced 44-methylgambierone and gambierone, and one isolate produced ciguatoxins, the cause of Ciguatera Poisoning. An Ostreopsis tairoto isolate, as analysed by the oxidative cleavage method, produced a palytoxin (PLTX)-like amine oxidation fragment, but when analysed for PLTX-like analogues using a new intact method none were detected indicating an 'unknown' PLTX-like compound is produced by this isolate. Isolates of O. cf. siamensis (Ostreopsis sp. 9), collected in Northland, were also analysed using the oxidative cleavage method, with the common PLTX-like amine fragment and the amide fragment corresponding to bishomoPLTX detected in all isolates. Again, the intact method indicated no detections in the isolates, again suggesting an unknown compound was being produced by these isolates. Prorocentrum hoffmannianum isolates produced okadaic acid (OA) and isoDTX-1 and P. lima isolates produced OA, DTX-1, and isoDTX-1. It is expected that new species of potentially harmful, benthic dinoflagellates will continue to be recorded in Aotearoa New Zealand and the results from Rangitāhua provide a guide to the HAB species to expect in sub-tropical Northland as the oceans continue to warm.

Exposure to Alexandrium spp. impairs the development of Green-lipped mussel (Perna canaliculus) embryos and larvae

The green-lipped mussel (GLM) Perna canaliculus is an economically, ecologically, and culturally important species in Aotearoa New Zealand. Since 2011, harmful algal blooms (HABs) of Alexandrium spp. have occurred annually in the Marlborough Sounds, the largest GLM aquaculture region in New Zealand. Across a similar timeframe, there has been a severe reduction in wild spat (juvenile mussel) catch. This research investigated the effects of Alexandrium pacificum (which produces paralytic shellfish toxins; PSTs) and A. minutum (a non-producer of PSTs) on the development of four GLM larval life stages (gametes, embryos, D-stage and settlement). Early life stages of GLM were exposed to environmentally relevant concentrations of Alexandrium spp. as whole cell, lysate and filtrate treatments. A 48-h exposure of embryos to whole A. pacificum cells at 500 cells mL-1 caused lysis of embryos, severe abnormalities, and reduced development through to veliger (D-stage) larvae by 85%. GLM growth was impaired at cell concentrations as low as 250 cells mL-1 during a 4-day exposure of D-stage larvae to both Alexandrium spp. Exposure of GLM to both whole and lysed treatments of Alexandrium spp. at 500 cells mL-1 resulted in halved larval growth rates (2.00 µm day-1 vs 4.48 µm day-1 in the control) and growth remained impeded during a 4-day recovery period. Both A. pacificum and A. minutum were found to negatively impact D-larvae. Both whole-cell and lysed-cell treatments of A. pacificum had similar negative effects, suggesting that Alexandrium spp. toxicity to D-larvae is independent of PSTs. Additionally, cell membrane-free treatments of A. pacificum had no negative effects on embryo development, indicating that cell surface-associated bioactive compounds may be responsible for the observed negative effects during this early life stage. Conversely, non-PST-producing A. minutum was toxic to D-stage larvae but not to embryos; larval growth was reduced following a brief 1 h exposure of sperm to cell membrane-free treatments of A. pacificum. No effects were recorded in GLM larvae exposed during settlement, highlighting the potential for differences in susceptibility of early life stages to Alexandrium spp. exposure and the influence of exposure durations. In the wild, blooms of Alexandrium spp. can persist for several months, reaching cell densities higher than those investigated in the present study, and as such may be detrimental to the vulnerable early life stages of GLM.

Comparison between droplet digital PCR and reverse transcription-quantitative PCR methods to measure ecotoxicology biomarkers

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is currently the gold-standard technique for detecting and quantifying messenger RNA. However, without proper validation, the method may produce artefactual and non-reproducible cycle threshold values generating poor-quality data. The newer droplet digital PCR (ddPCR) method allows for the absolute quantification of targeted nucleic acids providing more sensitive and accurate measurements without requiring external standards. This study compared these two PCR-based methods to measure the expression of well-documented genes used in ecotoxicology studies. We exposed Mediterranean mussels (Mytilus galloprovincialis) to copper and analyzed gene expression in gills and digestive glands using RT-qPCR and ddPCR assays. A step-by-step methodology to optimize and compare the two technologies is described. After ten-fold serial complementary DNA dilution, both RT-qPCR and ddPCR exhibited comparable linearity and efficiency and produced statistically similar results. We conclude that ddPCR is a suitable method to assess gene expression in an ecotoxicological context. However, RT-qPCR has a shorter processing time and remains more cost-effective.

Developing an integrated genomic selection approach beyond biomass for varietal protection and nutritive traits in perennial ryegrass (Lolium perenne L.)

Breeding target traits can be broadened to include nutritive value and plant breeder's rights traits in perennial ryegrass by using in-field regression-based spectroscopy phenotyping and genomic selection. Perennial ryegrass breeding has focused on biomass yield, but expansion into a broader set of traits is needed to benefit livestock industries whilst also providing support for intellectual property protection of cultivars. Numerous breeding objectives can be targeted simultaneously with the development of sensor-based phenomics and genomic selection (GS). Of particular interest are nutritive value (NV), which has been difficult and expensive to measure using traditional phenotyping methods, resulting in limited genetic improvement to date, and traits required to obtain varietal protection, known as plant breeder's rights (PBR) traits. In order to assess phenotyping requirements for NV improvement and potential for genetic improvement, in-field reflectance-based spectroscopy was assessed and GS evaluated in a single population for three key NV traits, captured across four timepoints. Using three prediction approaches, the possibility of targeting PBR traits using GS was evaluated for five traits recorded across three years of a breeding program. Prediction accuracy was generally low to moderate for NV traits and moderate to high for PBR traits, with heritability highly correlated with GS accuracy. NV did not show significant or consistent correlation between timepoints highlighting the need to incorporate seasonal NV into selection indexes and the value of being able to regularly monitor NV across seasons. This study has demonstrated the ability to implement GS for both NV and PBR traits in perennial ryegrass, facilitating the expansion of ryegrass breeding targets to agronomically relevant traits while ensuring necessary varietal protection is achieved.

Assessment of the recovery and photosynthetic efficiency of Breviolum psygmophilum and Effrenium voratum (Symbiodiniaceae) following cryopreservation

Many strains of Symbiodiniaceae have been isolated and their genetics, taxonomy, and metabolite production studied. Maintaining these cultures requires careful and regular sub-culturing that is costly with a high risk of species contamination or loss. Cryopreservation is a viable alternative for their long-term storage; however, there is uncertainty as to whether cryopreservation impacts the photosynthetic performance of Symbiodiniaceae. We investigated the growth rates and photosynthetic efficiency of two species, Breviolum psygmophilum and Effrenium voratum before and after cryopreservation. Rapid light curves (RLCs) produced using Pulse Amplitude Modulated (PAM) fluorometry were used to generate detailed information on the characteristics of photosystem II (PSII). The maximum electron transport rate (ETRmax) and the quantum yield (F_v/F_m) of the control (non-cryopreserved) and cryopreserved culture isolates were assessed across the growth cycle. The non-cryopreserved isolate of B. psygmophilum had a higher quantum yield than the cryopreserved isolate from day 12 to day 24, whereas there were no differences from day 28 to the late stationary phase. There were no significant differences in ETRmax. No significant differences were observed in quantum yield or ETRmax between the control and cryopreserved E. voratum isolates. The ability of cryopreserved strains to recover and regain their photosynthetic efficiency after freezing demonstrates the utility of this method for the long-term storage of these and other Symbiodiniaceae species.

Prokaryotic symbiont communities in three ascidian species introduced in both Ireland and New Zealand

Ascidians or sea squirts are among the marine taxa with the most introduced species worldwide. These animals have a suite of biological characteristics that contribute to their successful establishment, including long reproductive seasons, rapid growth rates, and resistance to pollution. Here, we sequenced a fragment of the 16S ribosomal RNA gene to characterize symbiont diversity and host-specificity in the solitary species Syela clava and Ascidiella aspersa, and the colonial species Didemnum vexillum. Samples were collected from introduced populations in several marinas and mussel facilities around Ireland, and a marina in New Zealand. Two additional colonial species Botrylloides violaceus and Didemnum sp. were collected in Ireland, and ambient seawater was sampled from both countries for comparison. Data revealed a strong effect of host species and location on prokaryote symbiont composition, consistent with recent ascidian microbiome literature. However, a location effect did not manifest in alpha diversity metrics (e.g., the same ascidian species at different locations exhibited similar diversity) but was evident in beta diversity metrics (greater intra-specific differences across locations than within locations). Location effects were stronger than species effects only for the solitary species (i.e., A. aspersa from New Zealand was more similar to S. clava from New Zealand than to A. aspersa from Ireland). D. vexillum and A. aspersa hosted a high abundance of prokaryotic symbionts that were previously found in other ascidian species, while S. clava symbiotic community was more closely related to bacteria common in the marine environment. Further studies should aim to unravel host-microbe coevolutionary patterns and the microbial role in facilitating host establishment in different habitats.

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.

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.

Development of a novel TaqMan qPCR assay for rapid detection and quantification of Gymnodinium catenatum for application to harmful algal bloom monitoring in coastal areas of Tunisia

Gymnodinium catenatum is a dinoflagellate known to cause paralytic shellfish poisoning (PSP), commonly associated with human muscular paralysis, neurological symptoms, and, in extreme cases, death. In the present work, we developed a real-time PCR-based assay for the rapid detection of the toxic microalgal species, G. catenatum, in environmental bivalve mollusc samples as well as seawater samples. G. catenatum-specific primers and probe were designed on the ITS1-5.8S-ITS2 rDNA region. Hydrolysis probe qPCR assay was optimized. ITS1-5.8S-ITS2 rDNA region copy numbers per G. catenatum cell genome were estimated to be 122.73 ± 5.54 copies/cell, allowing cell quantification. The application of the optimized qPCR assay for G. catenatum detection and quantification in field samples has been conducted, revealing high sensitivity (detection of around 1.310⁵ cells/L of seawater samples. Thus, the designed hydrolysis probe qPCR assay could be considered an efficient tool for phytoplankton monitoring whilst ensuring accuracy and sensitivity and providing cost and time savings.

Cryopreservation of six Symbiodiniaceae genera and assessment of fatty acid profiles in response to increased salinity treatments

Symbiodiniaceae are a diverse group of dinoflagellates, the majority of which are free-living and/or associated with a variety of protists and other invertebrate hosts. Maintenance of isolated cultures is labour-intensive and expensive, and cryopreservation provides an excellent avenue for their long-term storage. We aimed to cryopreserve 15 cultured isolates from six Symbiodiniaceae genera using dimethyl sulfoxide (DMSO) as the cryoprotectant agent (CPA). Under 15% DMSO, 10 isolates were successfully cryopreserved using either rapid freezing or controlled-rate freezing. Cultures that failed or had low survival, were subjected to (1) a reduction of CPA to 10%, or (2) increased salinity treatment before freezing. At 10% DMSO, three further isolates were successfully cryopreserved. At 15% DMSO there were high cell viabilities in Symbiodinium pilosum treated with 44 parts per thousand (ppt) and 54 ppt culture medium. An isolate of Fugacium sp. successfully cryopreserved after salinity treatments of 54 ppt and 64 ppt. Fatty acid (FA) analyses of S. pilosum after 54 ppt salinity treatment showed increased saturated FA levels, whereas Fugacium sp. had low poly-unsaturated FAs compared to normal salinity (34 ppt). Understanding the effects of salinity and roles of FAs in cryopreservation will help in developing protocols for these ecologically important taxa.

Unraveling the Gonyaulax baltica Species Complex: Cyst-theca Relationship of Impagidinium variaseptum, Spiniferites pseudodelicatus sp. nov. and S. ristingensis (Gonyaulacaceae, Dinophyceae), With Descriptions of Gonyaulax bohaiensis sp. nov, G. amoyensis sp. nov. and G. portimonensis sp. nov

The taxonomy of the extant dinoflagellate genus Gonyaulax is challenging since its thecate morphology is rather conservative. In contrast, cysts of Gonyaulax are varied in morphology and have been related with the fossil-based genera Spiniferites and Impagidinium. To better understand the systematics of Gonyaulax species, we performed germination experiments on cysts that can be identified as S. ristingensis, an unidentified Spiniferites with petaloid processes here described as Spiniferites pseudodelicatus sp. nov. and Impagidinium variaseptum from Chinese and Portuguese waters. Despite marked differences in cyst morphology, motile cells of S. pseudodelicatus and I. variaseptum are indistinguishable from Gonyaulax baltica. Motile cells hatched from S. ristingensis are morphologically similar to G. baltica as well but differ in the presence of one pronounced antapical spine. Three new species, Gonyaulax amoyensis (cyst equivalent S. pseudodelicatus), Gonyaulax bohaiensis (cyst equivalent I. variaseptum), and Gonyaulax portimonensis (cyst equivalent S. ristingensis), were erected. In addition, a new ribotype (B) of G. baltica was reported from South Korea and a bloom of G. baltica ribotype B is reported from New Zealand. Molecular phylogeny based on LSU and SSU rRNA gene sequences revealed that Gonyaulax species with minute or short antapical spines formed a well-resolved clade, whereas species with two pronounced antapical spines or lack of antapical spines formed the sister clade. Six strains of four above species were examined for yessotoxin production by liquid chromatography coupled with tandem mass spectrometry, and very low concentrations of yessotoxin were detected for one G. bohaiensis strain.

First transcriptome of the copepod Gladioferens pectinatus subjected to chronic contaminant exposures

Contaminants are often at low concentrations in ecosystems and their effects on exposed organisms can occur over long periods of time and across multiple generations. Alterations to subcellular mechanistic pathways in response to exposure to contaminants can provide insights into mechanisms of toxicity that methods measuring higher levels of biological may miss. Analysis of the whole transcriptome can identify novel mechanisms of action leading to impacts in exposed biota. The aim of this study was to characterise how exposures to copper, benzophenone and diclofenac across multiple generations altered molecular expression pathways in the marine copepod Gladioferens pectinatus. Results of the study demonstrated differential gene expression was observed in cultures exposure to diclofenac (569), copper (449) and benzophenone (59). Pathways linked to stress, growth, cellular and metabolic processes were altered by exposure to all three contaminants with genes associated with oxidative stress and xenobiotic regulation also impacted. Protein kinase functioning, cytochrome P450, transcription, skeletal muscle contraction/relaxation, mitochondrial phosphate translocator, protein synthesis and mitochondrial methylation were all differentially expressed with all three chemicals. The results of the study also suggested that using dimethyl sulfoxide as a dispersant influenced the transcriptome and future research may want to investigate it's use in molecular studies. Data generated in this study provides a first look at transcriptomic response of G. pectinatus exposed to contaminants across multiple generations, future research is needed to validate the identified biomarkers and link these results to apical responses such as population growth to demonstrate the predictive capacity of molecular tools.

q-PCR-based assay for the toxic dinoflagellate Karenia selliformis monitoring along the Tunisian coasts

Karenia selliformis is a marine dinoflagellate responsible for fish-kill events. Its presence has been reported along the Tunisian coasts (south-eastern Mediterranean Sea) since the 1990s. In the present study, a quantitative-PCR assay, based on the internal transcribed spacer (ITS) molecular marker, was developed to detect and quantify K. selliformis in environmental bivalve mollusk samples and in seawater samples. The assay was optimized, and its specificity was confirmed using cross-reactivity experiments against microalgal species commonly found on the Tunisian coasts and/or closely related to K. selliformis. Calibration curves were performed by tenfold dilutions of plasmid DNA harboring target sequence and genomic DNA, attaining a limit of detection of around 5 copies of target DNA per reaction, far below one K. selliformis cell per reaction. The field application of the developed assay showed a powerful detection capability. Thus, the designed assay could contribute to the deployment of in-field diagnostic tools for K. selliformis blooms monitoring.

Update of the Planktonic Diatom Genus Pseudo-nitzschia in Aotearoa New Zealand Coastal Waters: Genetic Diversity and Toxin Production

Domoic acid (DA) is produced by almost half of the species belonging to the diatom genus Pseudo-nitzschia and causes amnesic shellfish poisoning (ASP). It is, therefore, important to investigate the diversity and toxin production of Pseudo-nitzschia species for ASP risk assessments. Between 2018 and 2020, seawater samples were collected from various sites around Aotearoa New Zealand, and 130 clonal isolates of Pseudo-nitzschia were established. Molecular phylogenetic analysis of partial large subunit ribosomal DNA and/or internal transcribed spacer regions revealed that the isolates were divided into 14 species (Pseudo-nitzschia americana, Pseudo-nitzschia arenysensis, Pseudo-nitzschia australis, Pseudo-nitzschia calliantha, Pseudo-nitzschia cuspidata, Pseudo-nitzschia delicatissima, Pseudo-nitzschia fraudulenta, Pseudo-nitzschia galaxiae, Pseudo-nitzschia hasleana, Pseudo-nitzschia multiseries, Pseudo-nitzschia multistriata, Pseudo-nitzschia plurisecta, Pseudo-nitzschia pungens, and Pseudo-nitzschia cf. subpacifica). The P. delicatissima and P. hasleana strains were further divided into two clades/subclades (I and II). Liquid chromatography-tandem mass spectrometry was used to assess the production of DA and DA isomers by 73 representative strains. The analyses revealed that two (P. australis and P. multiseries) of the 14 species produced DA as a primary analogue, along with several DA isomers. This study is the first geographical distribution record of P. arenysensis, P.cuspidata, P. galaxiae, and P. hasleana in New Zealand coastal waters.

Development and application of a real-time PCR assay for the sensitive detection of diarrheic toxin producer Prorocentrum lima

Prorocentrum lima (P. lima) is a widely spread dinoflagellate in the Mediterranean Sea and it has become increasingly involved in harmful algal blooms. The purpose of this study is to develop a probe-based real-time polymerase chain reaction (PCR) targeting the ITS1-5.8S-ITS2 region for the detection and absolute quantification of P. lima based on linear and circular DNA standards. The results have shown that the quantitative PCR (q-PCR), using circular plasmid as a template, gave a threshold cycle number 1.79-5.6 greater than equimolar linear standards. When microalgae, commonly found in aquatic samples were tested, no cross-amplification was observed. The q-PCR brought about a good intra and inter-run reproducibility and a detection limit of 5 copies of linear plasmid per reaction. A quantitative relationship between the cell numbers and their corresponding plasmid copy numbers was attained. Afterwards, the effectiveness of the developed protocol was tested with 130 aquatic samples taken from 19 Tunisian sampling sites. The developed q-PCR had a detection sensitivity of up to 1 cell. All the positive samples were taken from three sampling sites of Medenine Governorate with cell abundances that ranged from 22 to 156,000 cells L-1 of seawater. The q-PCR assay revealed a high sensitivity in monitoring the aquatic samples in which the low concentrations of P. lima were not accurately detected by light microscopy. Indeed, this approach is at the same time rapid, specific and sensitive than the traditional microscopy techniques and it represents a great potential for the monitoring of P. lima blooms.

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.

Taxonomy and Molecular Phylogenetics of Ensiculiferaceae, fam. nov. (Peridiniales, Dinophyceae), with Consideration of their Life-history

In the current circumscription, the Thoracosphaeraceae comprise all dinophytes exhibiting calcified coccoid cells produced during their life-history. Species hitherto assigned to Ensiculifera and Pentapharsodinium are mostly based on the monadoid stage of life-history, while the link to the coccoid stage (occasionally treated taxonomically distinct) is not always resolved. We investigated the different life-history stages and DNA sequence data of Ensiculifera mexicana and other species occurring in samples collected from all over the world. Based on concatenated ribosomal RNA gene sequences Ensiculiferaceae represented a distinct peridinalean branch, which showed a distant relationship to other calcareous dinophytes. Both molecular and morphological data (particularly of the coccoid stage) revealed the presence of three distinct clades within Ensiculiferaceae, which may include other dinophytes exhibiting a parasitic life-history stage. At a higher taxonomic level, Ensiculiferaceae showed relationships to parasites and endosymbionts (i.e., Blastodinium and Zooxanthella) as well as to dinophytes harbouring diatoms instead of chloroplasts. These unexpected phylogenetic relationships are corroborated by the presence of five cingular plates in all such taxa, which differs from the six cingular plates of most other Thoracosphaeraceae. We herein describe Ensiculiferaceae, emend the descriptions of Ensiculifera and Pentapharsodinium, erect Matsuokaea and provide several new combinations at the species level.

Morphological and phylogenetic data do not support the split of Alexandrium into four genera

A recently published study analyzed the phylogenetic relationship between the genera Centrodinium and Alexandrium, confirming an earlier publication showing the genus Alexandrium as paraphyletic. This most recent manuscript retained the genus Alexandrium, introduced a new genus Episemicolon, resurrected two genera, Gessnerium and Protogonyaulax, and stated that: "The polyphyly [sic] of Alexandrium is solved with the split into four genera". However, these reintroduced taxa were not based on monophyletic groups. Therefore this work, if accepted, would result in replacing a single paraphyletic taxon with several non-monophyletic ones. The morphological data presented for genus characterization also do not convincingly support taxa delimitations. The combination of weak molecular phylogenetics and the lack of diagnostic traits (i.e., autapomorphies) render the applicability of the concept of limited use. The proposal to split the genus Alexandrium on the basis of our current knowledge is rejected herein. The aim here is not to present an alternative analysis and revision, but to maintain Alexandrium. A better constructed and more phylogenetically accurate revision can and should wait until more complete evidence becomes available and there is a strong reason to revise the genus Alexandrium. The reasons are explained in detail by a review of the available molecular and morphological data for species of the genera Alexandrium and Centrodinium. In addition, cyst morphology and chemotaxonomy are discussed, and the need for integrative taxonomy is highlighted.

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.

Assessment of Ciguatera and Other Phycotoxin-Related Risks in Anaho Bay (Nuku Hiva Island, French Polynesia): Molecular, Toxicological, and Chemical Analyses of Passive Samplers

Ciguatera poisoning is a foodborne illness caused by the consumption of seafood contaminated with ciguatoxins (CTXs) produced by dinoflagellates from the genera Gambierdiscus and Fukuyoa. The suitability of Solid Phase Adsorption Toxin Tracking (SPATT) technology for the monitoring of dissolved CTXs in the marine environment has recently been demonstrated. To refine the use of this passive monitoring tool in ciguateric areas, the effects of deployment time and sampler format on the adsorption of CTXs by HP20 resin were assessed in Anaho Bay (Nuku Hiva Island, French Polynesia), a well-known ciguatera hotspot. Toxicity data assessed by means of the mouse neuroblastoma cell-based assay (CBA-N2a) showed that a 24 h deployment of 2.5 g of resin allowed concentrating quantifiable amounts of CTXs on SPATT samplers. The CTX levels varied with increasing deployment time, resin load, and surface area. In addition to CTXs, okadaic acid (OA) and dinophysistoxin-1 (DTX1) were also detected in SPATT extracts using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS), consistent with the presence of Gambierdiscus and Prorocentrum species in the environment, as assessed by quantitative polymerase chain reaction (qPCR) and high-throughput sequencing (HTS) metabarcoding analyses conducted on passive window screen (WS) artificial substrate samples. Although these preliminary findings await further confirmation in follow-up studies, they highlight the usefulness of SPATT samplers in the routine surveillance of CP risk on a temporal scale, and the monitoring of other phycotoxin-related risks in ciguatera-prone areas.

Tetrodotoxin in marine bivalves and edible gastropods: A mini-review

Tetrodotoxin (TTX) is a potent neurotoxin responsible for countless human intoxications and deaths around the world. The distribution of TTX and its analogues is diverse and the toxin has been detected in organisms from both marine and terrestrial environments. Increasing detections seafood species, such as bivalves and gastropods, has drawn attention to the toxin, reinvigorating scientific interest and regulatory concerns. There have been reports of TTX in 21 species of bivalves and edible gastropods from ten countries since the 1980's. While TTX is structurally dissimilar to saxitoxin (STX), another neurotoxin detected in seafood, it has similar sodium channel blocking action and potency and both neurotoxins have been shown to have additive toxicities. The global regulatory level for the STX group toxins applied to shellfish is 800 μg/kg. The presence of TTX in shellfish is only regulated in one country; The Netherlands, with a regulatory level of 44 μg/kg. Due to the recent interest surrounding TTX in bivalves, the European Food Safety Authority established a panel to assess the risk and regulation of TTX in bivalves, and their final opinion was that a concentration below 44 μg of TTX per kg of shellfish would not result in adverse human effects. In this article, we review current knowledge on worldwide TTX levels in edible gastropods and bivalves over the last four decades, the different methods of detection used, and the current regulatory status. We suggest research needs that will assist with knowledge gaps and ultimately allow development of robust monitoring and management protocols.

Cryptic speciation in Protoceratium reticulatum (Dinophyceae): Evidence from morphological, molecular and ecophysiological data

The cosmopolitan, potentially toxic dinoflagellate Protoceratium reticulatum possesses a fossilizable cyst stage which is an important paleoenvironmental indicator. Slight differences in the internal transcribed spacer ribosomal DNA (ITS rDNA) sequences of P. reticulatum have been reported, and both the motile stage and cyst morphology of P. reticulatum display phenotypic plasticity, but how these morpho-molecular variations are related with ecophysiological preferences is unknown. Here, 55 single cysts or cells were isolated from localities in the Northern (Arctic to subtropics) and Southern Hemispheres (Chile and New Zealand), and in total 34 strains were established. Cysts and/or cells were examined with light microscopy and/or scanning electron microscopy. Large subunit ribosomal DNA (LSU rDNA) and/or ITS rDNA sequences were obtained for all strains/isolates. All strains/isolates of P. reticulatum shared identical LSU sequences except for one strain from the Mediterranean Sea that differs in one position, however ITS rDNA sequences displayed differences at eight positions. Molecular phylogeny was inferred using maximum likelihood and Bayesian inference based on ITS rDNA sequences. The results showed that P. reticulatum comprises at least three ribotypes (designated as A, B, and C). Ribotype A included strains from the Arctic and temperate areas, ribotype B included strains from temperate regions only, and ribotype C included strains from the subtropical and temperate areas. The average ratios of process length to cyst diameter of P. reticulatum ranged from 15% in ribotype A, 22% in ribotype B and 17% in ribotype C but cyst size could overlap. Theca morphology was indistinguishable among ribotypes. The ITS-2 secondary structures of ribotype A displayed one CBC (compensatory change on two sides of a helix pairing) compared to ribotypes B and C. Growth response of one strain from each ribotype to various temperatures was examined. The strains of ribotypes A, B and C exhibited optimum growth at 15 °C, 20 °C and 20-25 °C, respectively, thus corresponding to cold, moderate and warm ecotypes. The profiles of yessotoxins (YTXs) were examined for 25 strains using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The parent compound yessotoxin (YTX) was produced by strains of ribotypes A and B, but not by ribotype C strains, which only produced the structural variant homoyessotoxin (homoYTX). Our results support the notion that there is significant intra-specific variability in Protoceratium reticulatum and the biogeography of the different ribotypes is consistent with specific ecological preferences.

First description of the environmental niche of the epibenthic dinoflagellate species Coolia palmyrensis, C. malayensis, and C. tropicalis (Dinophyceae) from Eastern Australia

Environmental variables such as temperature, salinity, and irradiance are significant drivers of microalgal growth and distribution. Therefore, understanding how these variables influence fitness of potentially toxic microalgal species is particularly important. In this study, strains of the potentially harmful epibenthic dinoflagellate species Coolia palmyrensis, C. malayensis, and C. tropicalis were isolated from coastal shallow water habitats on the east coast of Australia and identified using the D1-D3 region of the large subunit (LSU) ribosomal DNA (rDNA). To determine the environmental niche of each taxon, growth was measured across a gradient of temperature (15-30°C), salinity (20-38), and irradiance (10-200 μmol photons · m^-2  · s^-1 ). Specific growth rates of Coolia tropicalis were highest under warm temperatures (27°C), low salinities (ca. 23), and intermediate irradiance levels (150 μmol photons · m^-2  · s^-1 ), while C. malayensis showed the highest growth at moderate temperatures (24°C) and irradiance levels (150 μmol photons · m^-2  · s^-1 ) and growth rates were consistent across the range of salinity levels tested (20-38). Coolia palmyrensis had the highest growth rate of all species tested and favored moderate temperatures (24°C), oceanic salinity (35), and high irradiance (>200 μmol photons · m^-2  · s^-1 ). This is the first study to characterize the environmental niche of species from the benthic harmful algal bloom genus Coolia and provides important information to help define species distributions and inform risk management.

High abundance of Amphidomataceae (Dinophyceae) during the 2015 spring bloom of the Argentinean Shelf and a new, non-toxigenic ribotype of Azadinium spinosum

Azaspiracids (AZA) are the most recently discovered group of lipophilic marine biotoxins of microalgal origin, and associated with human incidents of shellfish poisoning. They are produced by a few species of Amphidomataceae, but diversity and occurrence of the small-sized dinophytes remain poorly explored for many regions of the world. In order to analyze the presence and importance of Amphidomataceae in a highly productive area of Argentinean coastal waters (El Rincón area, SW Atlantic), a scientific cruise was performed in 2015 to sample the early spring bloom. In a multi-method approach, light microscopy was combined with real-time PCR molecular detection of Amphidomataceae, with chemical analysis of AZA, and with the establishment and characterization of amphidomatacean strains. Both light microscopy and PCR revealed that Amphidomataceae were widely present in spring plankton communities along the El Rincón area. They were particularly abundant offshore at the shelf front, reaching peak densities of 2.8 × 10⁵ cells L^-1, but no AZA were detected in field samples. In total, 31 new strains were determined as Az. dalianense and Az. spinosum, respectively. All Az. dalianense were non-toxigenic and shared the same rRNA sequences. The large majority of the new Az. spinosum strains revealed for the first time the presence of a non-toxigenic ribotype of this species, which is otherwise the most important AZA producer in European waters. One of the new Az. spinosum strains, with a particular slender shape and some other morphological peculiarities, clustered with toxigenic strains of Az. spinosum from Norway and, exceptionally for the species, produced only AZA-2 but not AZA-1. Results indicate a wide diversity within Az. spinosum, both in terms of sequence data and toxin profiles, which also will affect the qualitative and quantitative performance of the specific qPCR assay for this species. Overall, the new data provide a more differentiated perspective of diversity, toxin productivity and occurrence of Amphidomataceae in a poorly explored region of the global ocean.

Spatial variability and depuration of tetrodotoxin in the bivalve Paphies australis from New Zealand

Tetrodotoxin (TTX) is a potent neurotoxin responsible for many human intoxications globally. Despite its potency and widespread occurrence in taxonomically diverse species, the primary source of TTX remains uncertain. Paphies australis, an endemic clam found in New Zealand, has been found to contain TTX in several locations. However, it is unknown if this represents endogenous production or accumulation from an external source. To address this question, the concentrations of TTX in whole P. australis and dissected organs (siphons, foot, digestive gland and the ‘rest’) from thirteen sites around New Zealand were determined using liquid chromatography-tandem quadrupole mass spectrometry analysis (LC-MS/MS). Depuration rate of TTX was also investigated by harvesting and measuring concentrations in P. australis maintained in captivity on a toxin-free diet every three to 15 days for 150 days. The LC-MS/MS analyses of the spatial samples showed that TTX was present in P. australis from all regions tested, with significantly (p < 0.001) higher concentrations (15–50 μg kg⁻¹) observed at lower latitudes of the North Island compared with trace levels (0.5–3 μg kg⁻¹) in the South Island of New Zealand. Tetrodotoxin was detected in all the dissected organs but the siphons contained the highest concentrations of TTX at all sites analysed. A linear model of the depuration data identified a significant (p < 0.001) decline in total TTX concentrations in P. australis over the study period. The siphons maintained the highest amount of TTX across the entire depuration study. The digestive glands contained low concentrations at the start of the experiment, but this depurated rapidly and only traces remained after 21 days. These results provide evidence to suggest that P. australis does not produce TTX endogenously but obtains the neurotoxin from an exogenous source (e.g., diet) with the source more prevalent in warmer northern waters. The association of higher TTX concentrations in shellfish with warmer environments raises concerns that this toxin's distribution and abundance could become an increasing human health issue with global warming.

•

TTX-containing Paphies australis were maintained in captivity for 150 days and significantly depurated the toxin.

•

Thirteen populations of Paphies australis from around New Zealand were collected and tested for TTX.

•

All populations tested contained TTX but a significant latitudinal gradient was observed.

•

This study provides further evidence of an exogenous source of TTX in marine bivalves.

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

Paralytic shellfish toxin producing dinoflagellates have negatively impacted the shellfish aquaculture industry worldwide, including in Australia and New Zealand. Morphologically identical cryptic species of dinoflagellates that may differ in toxicity, in particular, species of the former Alexandrium tamarense species complex, co-occur in Australia, as they do in multiple regions in Asia and Europe. To understand the dynamics and the ecological drivers of the growth of each species in the field, accurate quantification at the species level is crucial. We have developed the first quantitative polymerase chain reaction (qPCR) primers for A. australiense, and new primers targeting A. ostenfeldii, A. catenella, and A. pacificum. We showed that our new primers for A. pacificum are more specific than previously published primer pairs. These assays can be used to quantify planktonic cells and cysts in the water column and in sediment samples with limits of detection of 2 cells/L for the A. catenella and A. australiense assays, 2 cells/L and 1 cyst/mg sediment for the A. pacificum assay, and 1 cells/L for the A. ostenfeldii assay, and efficiencies of >90%. We utilized these assays to discriminate and quantify co-occurring A. catenella, A. pacificum, and A. australiense in samples from the east coast of Tasmania, Australia.

Diversity, distribution, and azaspiracids of Amphidomataceae (Dinophyceae) along the Norwegian coast

Azaspiracids (AZA) are a group of lipophilic polyether compounds which have been implicated in shellfish poisoning incidents around Europe. They are produced by a few species of the dinophycean genera Azadinium and Amphidoma (Amphidomataceae). The presence of AZA toxins in Norway is well documented, but knowledge of the distribution and diversity of Azadinium and other Amphidomataceae along the Norwegian coast is rather limited and poorly documented. On a research survey along the Norwegian coast in 2015 from the Skagerrak in the South to Trondheimsfjorden in the North, plankton samples from 67 stations were analysed for the presence of Azadinium and Amphidoma and their respective AZA by on-board live microscopy, real-time PCR assays specific for Amphidomataceae, and liquid chromatography-tandem mass spectrometry (LC-MS/MS). Microscopy using live samples and positive real-time PCR assays using a general family probe and two species specific probes revealed the presence of Amphidomataceae distributed throughout the sampling area. Overall abundance was low, however, and was in agreement with a lack of detectable AZA in plankton samples. Single cell isolation and morphological and molecular characterisation of established strains revealed the presence of 7 amphidomatacean species (Azadiniun spinosum, Az. poporum, Az. obesum, Az. dalianense, Az. trinitatum, Az. polongum, Amphidoma languida) in the area. Azaspiracids were produced by the known AZA producing species Az. spinosum, Az. poporum and Am. languida only. LC-MS/MS analysis further revealed that Norwegian strains produce previously unreported AZA for Norway (AZA-11 by Az. spinosum, AZA-37 by Az. poporum, AZA-38 and AZA-39 by Am. languida), and also four novel compounds (AZA-50, -51 by Az. spinosum, AZA-52, -53 by Am. languida), whose structural properties are described and which now can be included in existing analytical protocols. A maximum likelihood analysis of concatenated rDNA regions (SSU, ITS1-ITS2, partial LSU) showed that the strains of Az. spinosum fell in two well supported clades, where most but not all new Norwegian strains formed the new Ribotype B. Ribotype differentiation was supported by a minor morphological difference with respect to the presence/absence of a rim around the pore plate, and was consistently reflected by different AZA profiles. Strains of Az. spinosum from ribotype A produce AZA-1, -2 and -33, whereas the new strains of ribotype B produce mainly AZA-11 and AZA-51. Significant sequence differences between both Az. spinosum ribotypes underline the need to redesign the currently used qPCR probes in order to detect all AZA producing Az. spinosum. The results generally underline the conclusion that for the Norwegian coast area it is important that amphidomatacean species are taken into account in future studies and monitoring programs.

Effects of chronic exposure to benzophenone and diclofenac on DNA methylation levels and reproductive success in a marine copepod

The UV-filter benzophenone and the anti-inflammatory diclofenac are commonly detected in the environment. The aim of this study was to assess the multigenerational effects of chronic exposure to low concentrations of these chemicals on toxicity and DNA methylation levels in the copepod Gladioferens pectinatus. Acute toxicity tests were conducted to determine the sensitivity of G. pectinatus to the chemicals. All chemicals impacted breeding, hatching and egg viability. Diclofenac (1 mg.L^-1) reduced the number of eggs per gravid female. Benzophenone (0.5 mg.L^-1) decreased egg hatching success. Exposure to the reference toxicant copper (0.02 mg.L^-1) led to unsuccessful hatching. Effects on DNA methylation was estimated by the percentage of 5- methylcytosine. The treatments resulted in strong differences in DNA methylation with increased methylation in the exposed animals. The two chemicals impacted both egg viability and the induction of differential DNA methylation, suggesting potential intra- and trans-generational evolutionary effects.

Distribution of Tetrodotoxin in the New Zealand Clam, Paphies australis, Established Using Immunohistochemistry and Liquid Chromatography-Tandem Quadrupole Mass Spectrometry

Tetrodotoxin (TTX) is one of the most potent neurotoxins known. It was originally thought to only occur in puffer fish but has now been identified in twelve different classes of freshwater and marine organisms, including bivalves. Despite being one of the world’s most studied biotoxins, its origin remains uncertain. There is contradictory evidence regarding the source of TTX and its pathway through food webs. To date, the distribution of TTX has not been examined in bivalves. In the present study, 48 Paphies australis, a TTX-containing clam species endemic to New Zealand, were collected. Thirty clams were dissected, and organs and tissues pooled into five categories (siphons, digestive gland, adductor muscles, and the ‘rest’) and analyzed for TTX using liquid chromatography-mass spectrometry (LC-MS). The micro-distribution of TTX was visualized in the remaining 18 individuals using an immunohistological technique incorporating a TTX-specific monoclonal antibody. The LC-MS analysis revealed that siphons contained the highest concentrations of TTX (mean 403.8 µg/kg). Immunohistochemistry analysis showed TTX in the outer cells of the siphons, but also in the digestive system, foot, and gill tissue. Observing TTX in organs involved in feeding provides initial evidence to support the hypothesis of an exogenous source in P. australis.

Effects of temperature and salinity stress on DNA methylation in a highly invasive marine invertebrate, the colonial ascidian Didemnum vexillum

Environmentally induced epigenetic changes may contribute to phenotypic plasticity, increase adaptive potential in changing environments, and play a key role in the establishment and spread of invasive species in new habitats. In this study, we used methylation-sensitive amplified polymorphism (MSAP) to assess environmentally induced DNA methylation changes in a globally invasive clonal ascidian, Didemnum vexillum. We tested the effect of increasing temperature (19, 25 and 27 °C) and decreasing salinity (34, 32, 30, 28 and 26 practical salinity units (PSU)) on global DNA methylation, growth and survival rates. Exposure to 27 °C resulted in significant changes in DNA methylation over time. Growth also decreased in colonies exposed to high temperatures, suggesting they were under thermal stress. In contrast, no differences in growth nor DNA methylation patterns were observed in colonies exposed to a decreasing salinity gradient, potentially due to prior adaptation. The results of this study show that environmental stress can induce significant global DNA methylation changes in an invasive marine invertebrate on very rapid timescales, and that this response varies depending on the type, magnitude, and duration of the stressor. Changes in genomic DNA methylation and the rate of growth may act to ‘buy survival time’ under stressful conditions, expanding the distribution limits of this globally invasive species.

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

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

The Epiphytic Genus Gambierdiscus (Dinophyceae) in the Kermadec Islands and Zealandia Regions of the Southwestern Pacific and the Associated Risk of Ciguatera Fish Poisoning

Species in the genus Gambierdiscus produce ciguatoxins (CTXs) and/or maitotoxins (MTXs), which may cause ciguatera fish poisoning (CFP) in humans if contaminated fish are consumed. Species of Gambierdiscus have previously been isolated from macroalgae at Rangitahua (Raoul Island and North Meyer Islands, northern Kermadec Islands), and the opportunity was taken to sample for Gambierdiscus at the more southerly Macauley Island during an expedition in 2016. Gambierdiscus cells were isolated, cultured, and DNA extracted and sequenced to determine the species present. Bulk cultures were tested for CTXs and MTXs by liquid chromatography-mass spectrometry (LC-MS/MS). The species isolated were G. australes, which produced MTX-1 (ranging from 3 to 36 pg/cell), and G. polynesiensis, which produced neither MTX-1 nor, unusually, any known CTXs. Isolates of both species produced putative MTX-3. The risk of fish, particularly herbivorous fish, causing CFP in the Zealandia and Kermadec Islands region is real, although in mainland New Zealand the risk is currently low. Both Gambierdiscus and Fukuyoa have been recorded in the sub-tropical northern region of New Zealand, and so the risk may increase with warming seas and shift in the distribution of Gambierdiscus species.

A new species of Gambierdiscus (Dinophyceae) from the south-west Pacific: Gambierdiscus honu sp. nov

Two isolates of a new tropical, epiphytic dinoflagellate species, Gambierdiscus honu sp. nov., were obtained from macroalgae sampled in Rarotonga, Cook Islands, and from North Meyer Island, Kermadec Islands. Gambierdiscus honu sp. nov. had the common Gambierdiscus Kofoidian plate formula: Po, 3', 6″, 6C?, 6 or 7S, 5‴, 1p and 2⁗. The characteristic morphological features of this species were its relatively small short dorsoventral length and width and the shape of individual plates, in particular the combination of the hatchet-shaped 2' and pentagonal 3' plates and the length to width ratio of the antapical 1p plate. The combination of these characteristics plus the smooth thecal surface and equal sized 1⁗ and 2⁗ plates differentiated this species from other Gambierdiscus species. The phylogenetic analyses supported the unique description. Both isolates of G. honu produced the putative maitotoxin (MTX)-3 analogue, but neither produced ciguatoxin (CTX) or MTX. Extracts of G. honu were shown to be highly toxic to mice by intraperitoneal injection (0.2mg/kg), although less toxic by gavage. It is possible that toxins other than putative MTX-3 are produced.

Role of Modular Polyketide Synthases in the Production of Polyether Ladder Compounds in Ciguatoxin-Producing Gambierdiscus polynesiensis and G. excentricus (Dinophyceae)

Gambierdiscus, a benthic dinoflagellate, produces ciguatoxins that cause the human illness Ciguatera. Ciguatoxins are polyether ladder compounds that have a polyketide origin, indicating that polyketide synthases (PKS) are involved in their production. We sequenced transcriptomes of Gambierdiscus excentricus and Gambierdiscus polynesiensis and found 264 contigs encoding single domain ketoacyl synthases (KS; G. excentricus: 106, G. polynesiensis: 143) and ketoreductases (KR; G. excentricus: 7, G. polynesiensis: 8) with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 24 contigs (G. excentricus: 3, G. polynesiensis: 21) encoding multiple PKS domains (forming typical type I PKSs modules) were found. The proposed structure produced by one of these megasynthases resembles a partial carbon backbone of a polyether ladder compound. Seventeen contigs encoding single domain KS, KR, s-malonyltransacylase, dehydratase and enoyl reductase with sequence similarity to type II fatty acid synthases (FAS) in plants were found. Type I PKS and type II FAS genes were distinguished based on the arrangement of domains on the contigs and their sequence similarity and phylogenetic clustering with known PKS/FAS genes in other organisms. This differentiation of PKS and FAS pathways in Gambierdiscus is important, as it will facilitate approaches to investigating toxin biosynthesis pathways in dinoflagellates.

A new Gambierdiscus species (Dinophyceae) from Rarotonga, Cook Islands: Gambierdiscus cheloniae sp. nov

Ciguatera fish poisoning (CFP) has been reported for many years in Rarotonga, Cook Islands, and has had the world's highest reported incidence of this illness for the last 20 years. Following intensive sampling to understand the distribution of the causative organisms of CFP, an undescribed Gambierdiscus species was isolated from the Rarotongan lagoon. Gambierdiscus cheloniae sp. nov. has the common Gambierdiscus Kofoidian plate formula (except for a variability in the number of precingular plates in aberrant cells): Po, 3', 6″ (7″), 6C?, 6 or 7S, 5'″, 1p and 2″″. The 2' plate is hatchet shaped and the dorsal end of 1p is pointed and the relatively narrow 1p plate. Morphologically G. cheloniae is similar to the genetically closely related species G. pacificus, G. toxicus and G. belizeanus, although smaller (depth and length) than G. toxicus. The apical pore plate varies from those of G. belizeanus and G. pacificus, which are shorter and narrower, and from G. toxicus, which is larger. G. cheloniae also differs from G. pacificus in the shape of the 2' plate. The description of this new species is supported by phylogenetic analyses using three different gene regions. G. cheloniae produced the putative maitotoxin-3 analogue, MTX-3, but neither maitotoxin or monitored ciguatoxin. Extracts of G. cheloniae were shown to be highly toxic to mice by intraperitoneal (i.p.) injection, although they were less toxic by gavage. It is possible that this species produces toxins other than putative MTX-3.

New scenario for speciation in the benthic dinoflagellate genus Coolia (Dinophyceae)

In this study, inter- and intraspecific genetic diversity within the marine harmful dinoflagellate genus Coolia Meunier was evaluated using isolates obtained from the tropics to subtropics in both Pacific and Atlantic Ocean basins. The aim was to assess the phylogeographic history of the genus and to clarify the validity of established species including Coolia malayensis. Phylogenetic analysis of the D1-D2 LSU rDNA sequences identified six major lineages (L1-L6) corresponding to the morphospecies Coolia malayensis (L1), C. monotis (L2), C. santacroce (L3), C. palmyrensis (L4), C. tropicalis (L5), and C. canariensis (L6). A median joining network (MJN) of C. malayensis ITS2 rDNA sequences revealed a total of 16 haplotypes; however, no spatial genetic differentiation among populations was observed. These MJN results in conjunction with CBC analysis, rDNA phylogenies and geographical distribution analyses confirm C. malayensis as a distinct species which is globally distributed in the tropical to warm-temperate regions. A molecular clock analysis using ITS2 rDNA revealed the evolutionary history of Coolia dated back to the Mesozoic, and supports the hypothesis that historical vicariant events in the early Cenozoic drove the allopatric differentiation of C. malayensis and C. monotis.

Rapid and accurate identification by real-time PCR of biotoxin-producing dinoflagellates from the family gymnodiniaceae

The identification of toxin-producing dinoflagellates for monitoring programmes and bio-compound discovery requires considerable taxonomic expertise. It can also be difficult to morphologically differentiate toxic and non-toxic species or strains. Various molecular methods have been used for dinoflagellate identification and detection, and this study describes the development of eight real-time polymerase chain reaction (PCR) assays targeting the large subunit ribosomal RNA (LSU rRNA) gene of species from the genera Gymnodinium, Karenia, Karlodinium, and Takayama. Assays proved to be highly specific and sensitive, and the assay for G. catenatum was further developed for quantification in response to a bloom in Manukau Harbour, New Zealand. The assay estimated cell densities from environmental samples as low as 0.07 cells per PCR reaction, which equated to three cells per litre. This assay not only enabled conclusive species identification but also detected the presence of cells below the limit of detection for light microscopy. This study demonstrates the usefulness of real-time PCR as a sensitive and rapid molecular technique for the detection and quantification of micro-algae from environmental samples.

Evaluating detection limits of next-generation sequencing for the surveillance and monitoring of international marine pests

Most surveillance programmes for marine invasive species (MIS) require considerable taxonomic expertise, are laborious, and are unable to identify species at larval or juvenile stages. Therefore, marine pests may go undetected at the initial stages of incursions when population densities are low. In this study, we evaluated the ability of the benchtop GS Junior™ 454 pyrosequencing system to detect the presence of MIS in complex sample matrices. An initial in-silico evaluation of the mitochondrial cytochrome c oxidase subunit I (COI) and the nuclear small subunit ribosomal DNA (SSU) genes, found that multiple primer sets (targeting a ca. 400 base pair region) would be required to obtain species level identification within the COI gene. In contrast a single universal primer set was designed to target the V1-V3 region of SSU, allowing simultaneous PCR amplification of a wide taxonomic range of MIS. To evaluate the limits of detection of this method, artificial contrived communities (10 species from 5 taxonomic groups) were created using varying concentrations of known DNA samples and PCR products. Environmental samples (water and sediment) spiked with one or five 160 hr old Asterias amurensis larvae were also examined. Pyrosequencing was able to recover DNA/PCR products of individual species present at greater than 0.64% abundance from all tested contrived communities. Additionally, single A. amurensis larvae were detected from both water and sediment samples despite the co-occurrence of a large array of environmental eukaryotes, indicating an equivalent sensitivity to quantitative PCR. NGS technology has tremendous potential for the early detection of marine invasive species worldwide.

Increased inter-colony fusion rates are associated with reduced COI haplotype diversity in an invasive colonial ascidian Didemnum vexillum

Considerable progress in our understanding of the population genetic changes associated with biological invasions has been made over the past decade. Using selectively neutral loci, it has been established that reductions in genetic diversity, reflecting founder effects, have occurred during the establishment of some invasive populations. However, some colonial organisms may actually gain an ecological advantage from reduced genetic diversity because of the associated reduction in inter-colony conflict. Here we report population genetic analyses, along with colony fusion experiments, for a highly invasive colonial ascidian, Didemnum vexillum. Analyses based on mitochondrial cytochrome oxidase I (COI) partial coding sequences revealed two distinct D. vexillum clades. One COI clade appears to be restricted to the probable native region (i.e., north-west Pacific Ocean), while the other clade is present in widely dispersed temperate coastal waters around the world. This clade structure was supported by 18S ribosomal DNA (rDNA) sequence data, which revealed a one base-pair difference between the two clades. Recently established populations of D. vexillum in New Zealand displayed greatly reduced COI genetic diversity when compared with D. vexillum in Japan. In association with this reduction in genetic diversity was a significantly higher inter-colony fusion rate between randomly paired New Zealand D. vexillum colonies (80%, standard deviation ±18%) when compared with colonies found in Japan (27%, standard deviation ±15%). The results of this study add to growing evidence that for colonial organisms reductions in population level genetic diversity may alter colony interaction dynamics and enhance the invasive potential of newly colonizing species.

Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars

Genetic map construction and identification of quantitative trait loci (QTLs) for blackleg resistance were performed for four mapping populations derived from five different canola source cultivars. Three of the populations were generated from crosses between single genotypes from the blackleg-resistant cultivars Caiman, Camberra and (AV)Sapphire and the blackleg-susceptible cultivar Westar(10). The fourth population was derived from a cross between genotypes from two blackleg resistant varieties (Rainbow and (AV)Sapphire). Different types of DNA-based markers were designed and characterised from a collection of 20,000 EST sequences generated from multiple Brassica species, including a new set of 445 EST-SSR markers of high value to the international community. Multiple molecular genetic marker systems were used to construct linkage maps with locus numbers varying between 219 and 468, and coverage ranging from 1173 to 1800 cM. The proportion of polymorphic markers assigned to map locations varied from 70 to 89% across the four populations. Publicly available simple sequence repeat markers were used to assign linkage groups to reference nomenclature, and a sub-set of mapped markers were also screened on the Tapidor x Ningyou (T x N) reference population to assist this process. QTL analysis was performed based on percentage survival at low and high disease pressure sites. Multiple QTLs were identified across the four mapping populations, accounting for 13-33% of phenotypic variance (V (p)). QTL-linked marker data are suitable for implementation in breeding for disease resistance in Australian canola cultivars. However, the likelihood of shifts in pathogen race structure across different geographical locations may have implications for the long-term durability of such associations.