Wanted dead or alive? Using metabarcoding of environmental DNA and RNA to distinguish living assemblages for biosecurity applications

High-throughput sequencing metabarcoding studies in marine biosecurity have largely focused on targeting environmental DNA (eDNA). DNA can persist extracellularly in the environment, making discrimination of living organisms difficult. In this study, bilge water samples (i.e., water accumulating on-board a vessel during transit) were collected from 15 small recreational and commercial vessels. eDNA and eRNA molecules were co-extracted and the V4 region of the 18S ribosomal RNA gene targeted for metabarcoding. In total, 62.7% of the Operational Taxonomic Units (OTUs) were identified at least once in the corresponding eDNA and eRNA reads, with 19.5% unique to eDNA and 17.7% to eRNA. There were substantial differences in diversity between molecular compartments; 57% of sequences from eDNA-only OTUs belonged to fungi, likely originating from legacy DNA. In contrast, there was a higher percentage of metazoan (50.2%) and ciliate (31.7%) sequences in the eRNA-only OTUs. Our data suggest that the presence of eRNA-only OTUs could be due to increased cellular activities of some rare taxa that were not identified in the eDNA datasets, unusually high numbers of rRNA transcripts in ciliates, and/or artefacts produced during the reverse transcriptase, PCR and sequencing steps. The proportions of eDNA/eRNA shared and unshared OTUs were highly heterogeneous within individual bilge water samples. Multiple factors including boat type and the activities performed on-board, such as washing of scientific equipment, may play a major role in contributing to this variability. For some marine biosecurity applications analysis, eDNA-only data may be sufficient, however there are an increasing number of instances where distinguishing the living portion of a community is essential. For these circumstances, we suggest only including OTUs that are present in both eDNA and eRNA data. OTUs found only in the eRNA data need to be interpreted with caution until further research provides conclusive evidence for their origin.

Development of a real-time polymerase chain reaction assay for the detection of the invasive Mediterranean fanworm, Sabella spallanzanii, in environmental samples

The Mediterranean fanworm, Sabella spallanzanii Gmelin 1791, was first detected in the Southern Hemisphere in the 1990s and is now abundant in many parts of southern Australia and in several locations around northern New Zealand. Once established, it can proliferate rapidly, reaching high densities with potential ecological and economic impacts. Early detection of new S. spallanzanii incursions is important to prevent its spread, guide eradication or control efforts and to increase knowledge on the species' dispersal pathways. In this study, we developed a TaqMan probe real-time polymerase chain reaction assay targeting a region of the mitochondrial cytochrome oxidase I gene. The assay was validated in silico and in vitro using DNA from New Zealand and Australian Sabellidae with no cross-reactivity detected. The assay has a linear range of detection over seven orders of magnitude with a limit of detection reached at 12.4 × 10^-4 ng/μL of DNA. We analysed 145 environmental (water, sediment and biofouling) samples and obtained positive detections only from spiked samples and those collected at a port where S. spallanzanii is known to be established. This assay has the potential to enhance current morphological and molecular-based methods, through its ability to rapidly and accurately identify S. spallanzanii in environmental samples.

Metabarcoding monitoring analysis: the pros and cons of using co-extracted environmental DNA and RNA data to assess offshore oil production impacts on benthic communities

Sequencing environmental DNA (eDNA) is increasingly being used as an alternative to traditional morphological-based identification to characterize biological assemblages and monitor anthropogenic impacts in marine environments. Most studies only assess eDNA which, compared to eRNA, can persist longer in the environment after cell death. Therefore, eRNA may provide a more immediate census of the environment due to its relatively weaker stability, leading some researchers to advocate for the use of eRNA as an additional, or perhaps superior proxy for portraying ecological changes. A variety of pre-treatment techniques for screening eDNA and eRNA derived operational taxonomic units (OTUs) have been employed prior to statistical analyses, including removing singleton taxa (i.e., OTUs found only once) and discarding those not present in both eDNA and eRNA datasets. In this study, we used bacterial (16S ribosomal RNA gene) and eukaryotic (18S ribosomal RNA gene) eDNA- and eRNA-derived data from benthic communities collected at increasing distances along a transect from an oil production platform (Taranaki, New Zealand). Macro-infauna (visual classification of benthic invertebrates) and physico-chemical data were analyzed in parallel. We tested the effect of removing singleton taxa, and removing taxa not present in the eDNA and eRNA libraries from the same environmental sample (trimmed by shared OTUs), by comparing the impact of the oil production platform on alpha- and beta-diversity of the eDNA/eRNA-based biological assemblages, and by correlating these to the morphologically identified macro-faunal communities and the physico-chemical data. When trimmed by singletons, presence/absence information from eRNA data represented the best proxy to detect changes on species diversity for both bacteria and eukaryotes. However, assessment of quantitative beta-diversity from read abundance information of bacteria eRNA did not, contrary to eDNA, reveal any impact from the oil production activity. Overall, the data appeared more robust when trimmed by shared OTUs, showing a greater effect of the platform on alpha- and beta-diversity. Trimming by shared OTUs likely removes taxa derived from legacy DNA and technical artefacts introduced through reverse transcriptase, polymerase-chain-reaction and sequencing. Findings from our scoping study suggest that metabarcoding-based biomonitoring surveys should, if funds, time and expertise allow, be assessed using both eDNA and eRNA products.

Adsorption of Ten Microcystin Congeners to Common Laboratory-Ware Is Solvent and Surface Dependent

Cyanobacteria can produce heptapetides called microcystins (MC) which are harmful to humans due to their ability to inhibit cellular protein phosphatases. Quantitation of these toxins can be hampered by their adsorption to common laboratory-ware during sample processing and analysis. Because of their structural diversity (>100 congeners) and different physico-chemical properties, they vary in their adsorption to surfaces. In this study, the adsorption of ten different MC congeners (encompassing non-arginated to doubly-arginated congeners) to common laboratory-ware was assessed using different solvent combinations. Sample handling steps were mimicked with glass and polypropylene pipettes and vials with increasing methanol concentrations at two pH levels, before analysis by liquid chromatography-tandem mass spectrometry. We demonstrated that MC adsorb to polypropylene surfaces irrespective of pH. After eight successive pipet actions using polypropylene tips ca. 20% of the MC were lost to the surface material, which increased to 25%–40% when solutions were acidified. The observed loss was alleviated by changing the methanol (MeOH) concentration in the final solvent. The required MeOH concentration varied depending on which congener was present. Microcystins only adsorbed to glass pipettes (loss up to 30% after eight pipet actions) when in acidified aqueous solutions. The latter appeared largely dependent on the presence of ionizable groups, such as arginine residues.

Trophic state and geographic gradients influence planktonic cyanobacterial diversity and distribution in New Zealand lakes

Cyanobacteria are commonly associated with eutrophic lakes, where they often form blooms and produce toxins. However, they are a ubiquitous component of phytoplankton in lakes of widely varying trophic status. We hypothesised that cyanobacterial diversity would vary among lakes of differing trophic status, but that the relative importance of geographical and hydromorphological characteristics driving these patterns would differ across trophic groups. DNA from 143 New Zealand lakes that spanned a range of geographic, hydromorphological and trophic gradients was analysed using automated rRNA intergenic spacer analysis and screened for genes involved in cyanotoxin production. Statistical analysis revealed significant delineation among cyanobacterial communities from different trophic classes. Multivariate regression indicated that geographical features (latitude, longitude and altitude) were significant in driving cyanobacterial community structure; however, partitioning of their effects varied among trophic categories. High-throughput sequencing was undertaken on selected samples to investigate their taxonomic composition. The most abundant and diverse (71 operational taxonomic units) taxon across all lake types was the picocyanobacteria genus Synechococcus Cyanotoxins (microcystins n = 23, anatoxins n = 1) were only detected in eutrophic lowland lakes. Collectively, these data infer that increasing eutrophication of lakes will have broad-scale impacts on planktonic cyanobacteria diversity and the prevalence of cyanotoxins.

Molecular Characterisation and Co-cultivation of Bacterial Biofilm Communities Associated with the Mat-Forming Diatom Didymosphenia geminata

Didymosphenia geminata (Lyngbye) M. Schmidt is a stalked freshwater diatom that is expanding its range globally. In some rivers, D. geminata forms thick and expansive polysaccharide-dominated mats. Like other stalked diatoms, D. geminata cells attach to the substratum with a pad of adhesive extracellular polymeric substance. Research on D. geminata and other diatoms suggests that bacterial biofilm composition may contribute to successful attachment. The aim of this study was to investigate the composition and role of bacterial biofilm communities in D. geminata attachment and survival. Bacterial biofilms were collected at four sites in the main stem of a river (containing D. geminata) and in four tributaries (free of D. geminata). Samples were characterised using automated rRNA intergenic spacer analysis and high-throughput sequencing (HTS). Mat-associated bacteria were isolated and their effect on the early establishment of D. geminata cells assessed using co-culturing experiments. ARISA and HTS data showed differences in bacterial communities between samples with and without D. geminata at two of the four sites. Samples with D. geminata had a higher relative abundance of Sphingobacteria (p < 0.01) and variability in community composition was reduced. Analysis of the 76 bacteria isolated from the mat revealed 12 different strains representing 8 genera. Co-culturing of a Carnobacterium sp. with D. geminata reduced survival (p < 0.001) and attachment (p < 0.001) of D. geminata. Attachment was enhanced by Micrococcus sp. and Pseudomonas sp. (p < 0.001 and p < 0.01, respectively). These data provide evidence that bacteria play a role in the initial attachment and on-going survival of D. geminata, and may partly explain observed distribution patterns.

Application of a spectrofluorimetric tool (bbe BenthoTorch) for monitoring potentially toxic benthic cyanobacteria in rivers

Over the last decade reports of animal poisoning following accidental consumption of neurotoxin-producing benthic cyanobacteria (mainly Phormidium spp.) have increased. There is a need for rapid and cost-effective tools to survey benthic cyanobacteria. In this study we assessed the performance of the BenthoTorch, a fluorometric probe that provides in situ estimations of cyanobacteria, diatoms and green algae biomass in biofilms. Biofilms (n = 288) were analysed from two rivers in France and eight in New Zealand. Correlations between chlorophyll-a measured using the BenthoTorch and spectrophotometry were higher for thin (<2 mm) compared to thick (>2 mm) biofilms (r(2) = 0.58 and 0.27 respectively; p < 0.001). When cyanobacteria represented less than 50% of the total biomass (based on biovolumes), microscopic and BenthoTorch compositional estimations were significantly correlated (r(2) = 0.53, p < 0.001). Conversely, there was no correlation when cyanobacteria exceeded 50% of the total biomass. Under this scenario diatoms were overestimated. Our results suggest that the observed biases occur because the BenthoTorch only measures the upper biofilm layer and it underestimates the biomass of phycoerythrin-containing cyanobacteria. To improve the performance of this sensor and render it a useful tool for a rapid evaluation of benthic cyanobacterial biomass in rivers, we propose that: (i) the algorithms based on the LEDs responses currently available on this tool need revision, (ii) new excitation wavelengths should be included that allow the fingerprints of phycoerythrin-containing cyanobacteria to be discriminated, and (iii) a sensor that penetrates the biofilms is needed to obtain more accurate estimates of cyanobacterial biomass.

First evaluation of foraminiferal metabarcoding for monitoring environmental impact from an offshore oil drilling site

At present, environmental impacts from offshore oil and gas activities are partly determined by measuring changes in macrofauna diversity. Morphological identification of macrofauna is time-consuming, expensive and dependent on taxonomic expertise. In this study, we evaluated the applicability of using foraminiferal-specific metabarcoding for routine monitoring. Sediment samples were collected along distance gradients from two oil platforms off Taranaki (New Zealand) and their physico-chemical properties, foraminiferal environmental DNA/RNA, and macrofaunal composition analyzed. Macrofaunal and foraminiferal assemblages showed similar shifts along impact gradients, but responded differently to environmental perturbations. Macrofauna were affected by hypoxia, whereas sediment grain size appeared to drive shifts in foraminifera. We identified eight foraminiferal molecular operational taxonomic units that have potential to be used as bioindicator taxa. Our results show that metabarcoding represents an effective tool for assessing foraminiferal communities near offshore oil and gas platforms, and that it can be used to complement current monitoring techniques.

Geographically conserved microbiomes of four temperate water tunicates

Tunicates are useful models for exploring microbiomes because they have an innate immune system resembling that of chordates. Automated ribosomal RNA intergenic spacer analysis and High-Throughput Sequencing were used to compare the tunic microbiomes of Ciona robusta (formerly Ciona intestinalis type A), Ciona savignyi, Botrylloides leachi and Botryllus schlosseri sampled from three distinct locations with limited genetic connectivity. Bacterial phylotype profiles were conserved within each species, and there were no detectable differences between tunic and tunic + cuticle subsamples from an individual. Bacterial operational taxonomic unit (OTU) diversity was lowest for C. savignyi (320 ± 190 OTUs) and highest for B. schlosseri (1260 ± 190 OTUs). Each species had a distinct set of bacterial OTUs (pseudo-F = 3.0, p > 0.001), with the exception of B. leachi and B. schlosseri from one sampling location (t = 1.2, p = 0.09). Of note were OTUs assigned to Alphaproteobacteria from C. robusta plus Phyllobacteriaceae and Endozoicomonas from C. savignyi. These OTUs contributed 51, 22 and 10% of sequence reads, respectively, and are related to known bacterial symbionts. The within-species conservation of core OTUs across three distinct and co-occurring populations of tunicates provides compelling evidence that these tunicates foster defined microbiomes.

Deubiquitylating enzyme, USP9X, regulates proliferation of cells of head and neck cancer lines

OBJECTIVES

Truncating mutations in USP9X have been identified in oral squamous cell carcinoma patients. The aim of this study was to determine USP9X's functional role, if any, in head and neck cancer cells.

MATERIALS AND METHODS

USP9X was depleted/overexpressed in head and neck cancer cell line: SCC15 (tongue), CAL27 (tongue), FaDu (pharynx) and Detroit 562 (pharynx). Cell proliferation was monitored using the CyQUANT assay, and cell cycle distribution was determined by flow cytometry. Immunoblot assays were conducted to assess protein levels. RT-qPCR was performed to determine Notch and Wnt pathway target gene expression.

RESULTS

Our data showed a direct correlation between USP9X protein levels and proliferation, as well as Notch pathway activity in head and neck cancer cells. However, at least in FaDu, USP9X did not appear to regulate proliferation through the Notch pathway. Immunoblotting revealed a dramatic reduction in downstream targets of mTOR complex 1, namely total ribosomal protein (S6) and its phosphorylated form (pS6), when USP9X was depleted in FaDu cells. In contrast, in immortalized but non-tumorigenic HaCaT keratinocytes, USP9X depletion led to increase in cell proliferation, maintaining direct regulation of Notch activity.

CONCLUSIONS

The functional role of USP9X was found to be context dependent. USP9X possibly promotes head and neck cancer cell proliferation through the mTOR pathway.

Metabarcoding improves detection of eukaryotes from early biofouling communities: implications for pest monitoring and pathway management

In this experimental study the patterns in early marine biofouling communities and possible implications for surveillance and environmental management were explored using metabarcoding, viz. 18S ribosomal RNA gene barcoding in combination with high-throughput sequencing. The community structure of eukaryotic assemblages and the patterns of initial succession were assessed from settlement plates deployed in a busy port for one, five and 15 days. The metabarcoding results were verified with traditional morphological identification of taxa from selected experimental plates. Metabarcoding analysis identified > 400 taxa at a comparatively low taxonomic level and morphological analysis resulted in the detection of 25 taxa at varying levels of resolution. Despite the differences in resolution, data from both methods were consistent at high taxonomic levels and similar patterns in community shifts were observed. A high percentage of sequences belonging to genera known to contain non-indigenous species (NIS) were detected after exposure for only one day.

The rise of toxic benthic Phormidium proliferations: A review of their taxonomy, distribution, toxin content and factors regulating prevalence and increased severity

There has been a marked increase in the distribution, intensity and frequency of proliferations of some species of the benthic mat-forming, toxin-producing genus Phormidium in rivers globally over the last decade. This review summarises current knowledge on their taxonomy, distribution, toxin content, environmental drivers of proliferations, and monitoring and management strategies in New Zealand. Although toxic Phormidium proliferation occurs in rivers worldwide little is known about these factors in most countries. Proliferations, defined as >20% cover of a riverbed, have been identified in 103 rivers across New Zealand. Morphological and molecular data indicate the main species responsible is Phormidium autumnale. In New Zealand Phormidium produces anatoxins (anatoxin-a, homoanatoxin-a, dihydroanatoxin-a, and dihydrohomoanatoxin-a) and these were detected in 67% of 771 samples from 40 rivers. The highest concentration measured was 712mgkg^-1 dried weight (Oreti River, Southland), with considerable spatial and temporal variability in anatoxin concentrations between and within rivers. A synthesis of field based studies suggests that Phormidium proliferations are most likely when there is some enrichment of dissolved inorganic nitrogen but when water-column dissolved reactive phosphorus is less than 0.01mgL^-1. Once established Phormidium-dominated mats trap sediment and internal mat biogeochemistry can mobilise sediment-bound phosphorus, which is then available for growth. Removal of Phormidium-dominated mats is primarily due to shear stress and substrate disturbance, although there is also evidence for autogenic detachment. A combination of factors including; changes to riparian margins, increased nitrate and fine sediment loads, and alterations in flow regimes are likely to have contributed to the rise in Phormidium proliferations.

A review of the global ecology, genomics, and biogeography of the toxic cyanobacterium, Microcystis spp

This review summarizes the present state of knowledge regarding the toxic, bloom-forming cyanobacterium, Microcystis, with a specific focus on its geographic distribution, toxins, genomics, phylogeny, and ecology. A global analysis found documentation suggesting geographic expansion of Microcystis, with recorded blooms in at least 108 countries, 79 of which have also reported the hepatatoxin microcystin. The production of microcystins (originally "Fast-Death Factor") by Microcystis and factors that control synthesis of this toxin are reviewed, as well as the putative ecophysiological roles of this metabolite. Molecular biological analyses have provided significant insight into the ecology and physiology of Microcystis, as well as revealed the highly dynamic, and potentially unstable, nature of its genome. A genetic sequence analysis of 27 Microcystis species, including 15 complete/draft genomes are presented. Using the strictest biological definition of what constitutes a bacterial species, these analyses indicate that all Microcystis species warrant placement into the same species complex since the average nucleotide identity values were above 95%, 16S rRNA nucleotide identity scores exceeded 99%, and DNA-DNA hybridization was consistently greater than 70%. The review further provides evidence from around the globe for the key role that both nitrogen and phosphorus play in controlling Microcystis bloom dynamics, and the effect of elevated temperature on bloom intensification. Finally, highlighted is the ability of Microcystis assemblages to minimize their mortality losses by resisting grazing by zooplankton and bivalves, as well as viral lysis, and discuss factors facilitating assemblage resilience.

Accurate assessment of the impact of salmon farming on benthic sediment enrichment using foraminiferal metabarcoding

Assessing the environmental impact of salmon farms on benthic systems is traditionally undertaken using biotic indices derived from microscopic analyses of macrobenthic infaunal (MI) communities. In this study, we tested the applicability of using foraminiferal-specific high-throughput sequencing (HTS) metabarcoding for monitoring these habitats. Sediment samples and physico-chemical data were collected along an enrichment gradient radiating out from three Chinook salmon (Oncorhynchus tshawytscha) farms in New Zealand. HTS of environmental DNA and RNA (eDNA/eRNA) resulted in 1,875,300 sequences that clustered into 349 Operational Taxonomic Units. Strong correlations were observed among various biotic indices calculated from MI data and normalized fourth-root transformed HTS data. Correlations were stronger using eRNA compared to eDNA data. Quantile regression spline analyses identified 12 key foraminiferal taxa that have potential to be used as bioindicator species. This study demonstrates the huge potential for using this method for biomonitoring of fish-farming and other marine industrial activities.

The effect of cyanobacterial biomass enrichment by centrifugation and GF/C filtration on subsequent microcystin measurement

Microcystins are cyclic peptides produced by multiple cyanobacterial genera. After accumulation in the liver of animals they inhibit eukaryotic serine/threonine protein phosphatases, causing liver disease or death. Accurate detection/quantification of microcystins is essential to ensure safe water resources and to enable research on this toxin. Previous methodological comparisons have focused on detection and extraction techniques, but have not investigated the commonly used biomass enrichment steps. These enrichment steps could modulate toxin production as recent studies have demonstrated that high cyanobacterial cell densities cause increased microcystin levels. In this study, three microcystin-producing strains were processed using no cell enrichment steps (by direct freezing at three temperatures) and with biomass enrichment (by centrifugation or GF/C filtration). After extraction, microcystins were analyzed using liquid chromatography-tandem mass spectrometry. All processing methods tested, except GF/C filtration, resulted in comparable microcystin quotas for all strains. The low yields observed for the filtration samples were caused by adsorption of arginine-containing microcystins to the GF/C filters. Whilst biomass enrichment did not affect microcystin metabolism over the time-frame of normal sample processing, problems associated with GF/C filtration were identified. The most widely applicable processing method was direct freezing of samples as it could be utilized in both field and laboratory environments.

Further characterization of glycine-containing microcystins from the McMurdo dry Valleys of Antarctica

Microcystins are hepatotoxic cyclic peptides produced by several cyanobacterial genera worldwide. In 2008, our research group identified eight new glycine-containing microcystin congeners in two hydro-terrestrial mat samples from the McMurdo Dry Valleys of Eastern Antarctica. During the present study, high-resolution mass spectrometry, amino acid analysis and micro-scale thiol derivatization were used to further elucidate their structures. The Antarctic microcystin congeners contained the rare substitution of the position-1 ᴅ-alanine for glycine, as well as the acetyl desmethyl modification of the position-5 Adda moiety (3S-amino-9S-methoxy-2S,6,8S-trimethyl-10-phenyldeca-4E,6E-dienoic acid). Amino acid analysis was used to determine the stereochemistry of several of the amino acids and conclusively demonstrated the presence of glycine in the microcystins. A recently developed thiol derivatization technique showed that each microcystin contained dehydrobutyrine in position-7 instead of the commonly observed N-methyl dehydroalanine.

Successional change in microbial communities of benthic Phormidium-dominated biofilms

Benthic cyanobacterial blooms are increasing worldwide and can be harmful to human and animal health if they contain toxin-producing species. Microbial interactions are important in the formation of benthic biofilms and can lead to increased dominance and/or toxin production of one or few taxa. This study investigated how microbial interactions contribute to proliferation of benthic blooms dominated by the neurotoxin-producing Phormidium autumnale. Following a rainfall event that cleared the substrate, biofilm succession was characterised at a site on the Hutt River (New Zealand) by sampling every 2-3 days over 32 days. A combination of morphological and molecular community analyses (automated ribosomal intergenic spacer analysis and Illumina™ MiSeq sequencing) identified three distinct phases of succession in both the micro-algal and bacterial communities within P. autumnale-dominated biofilms. Bacterial composition shifted between the phases, and these changes occurred several days before those of the micro-algal community. Alphaproteobacteria and Betaproteobacteria dominate in the early phase; Alphaproteobacteria, Betaproteobacteria, Sphingobacteria and Flavobacteria in the mid-phase; and Sphingobacteria and Flavobacteria in the late phase. Collectively, the results suggest that succession is driven by bacteria in the early stages but becomes dependent on micro-algae in the mid- and late stages of biofilm formation.

Intracellular immunohistochemical detection of tetrodotoxin in Pleurobranchaea maculata (Gastropoda) and Stylochoplana sp. (Turbellaria)

Tetrodotoxin (TTX), is a potent neurotoxin targeting sodium channels that has been identified in multiple marine and terrestrial organisms. It was recently detected in the Opisthobranch Pleurobranchaea maculata and a Platyhelminthes Stylochoplana sp. from New Zealand. Knowledge on the distribution of TTX within these organisms is important to assist in elucidating the origin and ecological role of this toxin. Intracellular micro-distribution of TTX was investigated using a monoclonal antibody-based immunoenzymatic technique. Tetrodotoxin was strongly localized in neutral mucin cells and the basement membrane of the mantle, the oocytes and follicles of the gonad tissue, and in the digestive tissue of P. maculata. The ova and pharynx were the only two structures to contain TTX in Stylochoplana sp. Using liquid chromatography-mass spectrometry, TTX was identified in the larvae and eggs, but not the gelatinous egg cases of P. maculata. Tetrodotoxin was present in egg masses of Stylochoplana sp. These data suggest that TTX has a defensive function in adult P. maculata, who then invest this in their progeny for protection. Localization in the digestive tissue of P. maculata potentially indicates a dietary source of TTX. Stylochoplana sp. may use TTX in prey capture and for the protection of offspring.

Early detection of eukaryotic communities from marine biofilm using high-throughput sequencing: an assessment of different sampling devices

Marine biofilms are precursors for colonization by larger fouling organisms, including non-indigenous species (NIS). In this study, high-throughput sequencing (HTS) of 18S rRNA metabarcodes was used to investigate four sampling methods (modified syringe, sterilized sponge, underwater tape and sterilized swab) for characterizing eukaryotic communities in marine biofilms. Perspex™ plates were sampled in and out of water. DNA collected with tape did not amplify. Otherwise, there were no statistical differences in communities among the remaining three sampling devices or between the two environments. Sterilized sponges are recommended for ease of use underwater. In-depth HTS analysis identified diverse eukaryotic communities, dominated by Metazoa and Chromoalveolata. Among the latter, diatoms (Bacillariophyceae) were particularly abundant (33% of reads assigned to Chromalveolata). The NIS Ciona savignyi was detected in all samples. The application of HTS in marine biofilm surveillance could facilitate early detection of NIS, improving the probability of successful eradication.

First identification of tetrodotoxin (TTX) in the flatworm Stylochoplana sp.; a source of TTX for the sea slug Pleurobranchaea maculata

High concentrations of the neurotoxin tetrodotoxin (TTX) were detected by liquid chromatography-mass spectrometry (LC-MS) in the Platyhelminthes Stylochoplana sp. from Pilot Bay (Tauranga, New Zealand). This is the first detection of TTX in this genus. Concentrations were monitored from March to November (2013) and found to significantly decrease from a peak in July (avg. 551 mg kg(-1)) to November (avg. 140 mg kg(-1)). Stylochoplana sp. co-occurred with TTX-containing Pleurobranchaea maculata (Opisthobranchia). A Stylochoplana sp.-specific real-time PCR assay was developed targeting the mitochondrial cytochrome c oxidase subunit I gene to determine if P. maculata consumed Stylochoplana sp. Positive Stylochoplana sp. signals were obtained for 7 of 19 P. maculata tested. Mass calculations indicate Stylochoplana sp. could supply Pilot Bay P. maculata with the TTX required to account for the concentrations reported in previous studies (ca. 1.04 mg TTX per individual) based on an ingestion rate of one individual every 2-3 days throughout their lifetime. However, due to the lack of Stylochoplana sp. in areas with dense P. maculata populations, and high concentration (ca. 1400 mg kg(-1)) of TTX detected in some individuals, it is unlikely that Stylochoplana sp. represent the sole source of TTX in P. maculata.