Using a weight of evidence approach to identify sources of microbiological contamination in a shellfish-growing area with "Restricted" classification.

Shellfish-growing areas in rural catchments are occasionally affected by elevated faecal contamination from diffuse sources and may be subject to frequent harvest closures/classification downgrades. We combined traditional risk management methods based on sanitary surveys and monitoring of Escherichia coli in seawater and shellfish with faecal source tracking, bacterial source apportionment, and hydrometeorological modelling to determine the causes of elevated E. coli concentrations contributing to harvest closures in Papanui Inlet (Aotearoa New Zealand). These multiple lines of evidence were used to inform a weight of evidence assessment of bacterial contamination in the inlet. Ruminant livestock was estimated to contribute 80% of the faecal coliform loading. Concentrations of E. coli in seawater were low (≤ 87 MPN 100 ml-1) whilst concentrations in tuaki/cockles/little neck clams (Austrovenus stutchburyi) occasionally exceeded the "Approved" classification limit (230 MPN 100 g-1). The most frequent positive genetic markers in seawater were the seagull (Catellicoccus marimammalium) (54% of seawater samples), and in shellfish, the bovine and seagull markers (both 12.5% of shellfish samples). Solar radiation was negatively correlated with E. coli in tuaki. We found that the growing area is affected by faecal inputs from animal and, to a lesser extent, human (septic tank discharges) sources which elevate contamination to levels detectable in tuaki but not in seawater, particularly in the summer months. The innovative approach can enhance the management of shellfish-growing areas affected by intermittent contamination and enables more targeted action to reduce pollution to improve shellfish water quality.

Growth conditions impact particulate absorption and pigment concentrations in two common bloom forming cyanobacterial species.

Remote sensing using satellite imagery has been promoted as a method to broaden the scale and frequency of cyanobacterial monitoring. This relies on the ability to establish relationships between the reflectance spectra of water bodies and the abundance of cyanobacteria. A challenge to achieving this comes from a limited understanding of the extent to which the optical properties of cyanobacteria vary according to their physiological state and growth environment. The aim of the present study was to determine how growth stage, nutrient status and irradiance affect pigment concentrations and absorption spectra in two common bloom forming cyanobacterial taxa: Dolichospermum lemmermannii and Microcystis aeruginosa. Each species was grown in laboratory batch culture under a full factorial design of low or high light intensity and low, medium, or high nitrate concentrations. Absorption spectra, pigment concentrations and cell density were measured throughout the growth phases. The absorption spectra were all highly distinguishable from each other, with greater interspecific than intraspecific differences, indicating that both D. lemmermannii and M. aeruginosa can be readily differentiated using hyperspectral absorption spectra. Despite this, each species exhibited different responses in the per-cell pigment concentrations with varying light intensity and nitrate exposure. Variability among treatments was considerably higher in D. lemmermannii than in M. aeruginosa, which exhibited smaller changes in pigment concentrations among the treatments. These results highlight the need to understand the physiology of the cyanobacteria and to take caution when estimating biovolumes from reflectance spectra when species composition and growth stage are unknown.

Growth at the limits: comparing trace metal limitation of a freshwater cyanobacterium (Dolichospermum lemmermannii) and a freshwater diatom (Fragilaria crotonensis).

Freshwater phytoplankton blooms are increasing in prevalence and there are conflicting views on whether trace metals limit growth of key species and thus bloom formation. The Taupo Volcanic Zone (TVZ), New Zealand, was formed by multiple eruptions of a super-volcano which emitted rhyolitic tephra leaving lakes depleted in trace metals. This provides an opportunity to test the potential of trace metal limitation on freshwater phytoplankton growth under nanomolar concentrations. Growth responses of two algal species isolated from Lake Taupo, Dolichospermum lemmermannii (cyanobacteria) and Fragilaria crotonensis (diatom), to six biologically important trace metals (manganese, iron, zinc, cobalt, copper and molybdenum) were examined in culture experiments. These were conducted at three trace metal concentrations: (1) ambient, (2) two-times ambient, and (3) ten-times ambient concentrations in Lake Taupo. Elevated concentrations of iron significantly increased growth rates and maximum cell densities in D. lemmermannii, whereas no significant concentration dependence was observed for other trace metals. Fragilaria crotonensis showed no significant growth response to elevated concentrations of trace metals. These results highlight the importance of iron as a growth limiting nutrient for cyanobacteria and indicate that even small (twofold) increases in Fe concentrations could enhance cyanobacteria growth rates in Lake Taupo, potentially causing cyanobacterial blooms.

Development of droplet digital Polymerase Chain Reaction assays for the detection of long-finned (Anguilla dieffenbachii) and short-finned (Anguilla australis) eels in environmental samples.

Freshwater eels are ecologically, and culturally important worldwide. The New Zealand long-finned eel (Anguilla dieffenbachii) and short-finned eel (Anguilla australis) are apex predators, playing an important role in ecosystem functioning of rivers and lakes. Recently, there has been a national decline in their populations due to habitat destruction and commercial harvest. The emergence of targeted environmental DNA detection methodologies provides an opportunity to enhance information about their past and present distributions. In this study we successfully developed species-specific droplet digital Polymerase Chain Reaction (ddPCR) assays to detect A. dieffenbachii and A. australis DNA in water and sediment samples. Assays utilized primers and probes designed for regions of the mitochondrial cytochrome b and 16S ribosomal RNA genes in A. dieffenbachii and A. australis, respectively. River water samples (n = 27) were analyzed using metabarcoding of fish taxa and were compared with the ddPCR assays. The presence of A. dieffenbachii and A. australis DNA was detected in a greater number of water samples using ddPCR in comparison to metabarcoding. There was a strong and positive correlation between gene copies (ddPCR analyses) and relative eel sequence reads (metabarcoding analyses) when compared to eel biomass. These ddPCR assays provide a new method for assessing spatial distributions of A. dieffenbachii and A. australis in a range of environments and sample types.

Phospholipid fatty acid (PLFA) analysis as a tool to estimate absolute abundances from compositional 16S rRNA bacterial metabarcoding data.

Microbial biodiversity monitoring through the analysis of DNA extracted from environmental samples is increasingly popular because it is perceived as being rapid, cost-effective, and flexible concerning the sample types studied. DNA can be extracted from diverse media before high-throughput sequencing of the prokaryotic 16S rRNA gene is used to characterize the taxonomic diversity and composition of the sample (known as metabarcoding). While sources of bias in metabarcoding methodologies are widely acknowledged, previous studies have focused mainly on the effects of these biases within a single substrate type, and relatively little is known of how these vary across substrates. We investigated the effect of substrate type (water, microbial mats, lake sediments, stream sediments, soil and a mock microbial community) on the relative performance of DNA metabarcoding in parallel with phospholipid fatty acid (PLFA) analysis. Quantitative estimates of the biomass of different taxonomic groups in samples were made through the analysis of PLFAs, and these were compared to the relative abundances of microbial taxa estimated from metabarcoding. Furthermore, we used the PLFA-based quantitative estimates of the biomass to adjust relative abundances of microbial groups determined by metabarcoding to provide insight into how the biomass of microbial taxa from PLFA analysis can improve understanding of microbial communities from environmental DNA samples. We used two sets of PLFA biomarkers that differed in their number of PLFAs to evaluate how PLFA biomarker selection influences biomass estimates. Metabarcoding and PLFA analysis provided significantly different views of bacterial composition, and these differences varied among substrates. We observed the most notable differences for the Gram-negative bacteria, which were overrepresented by metabarcoding in comparison to PLFA analysis. In contrast, the relative biomass and relative sequence abundances aligned reasonably well for Cyanobacteria across the tested freshwater substrates. Adjusting relative abundances of microbial taxa estimated by metabarcoding with PLFA-based quantification estimates of the microbial biomass led to significant changes in the microbial community compositions in all substrates. We recommend including independent estimates of the biomass of microbial groups to increase comparability among metabarcoding libraries from environmental samples, especially when comparing communities associated with different substrates.

Differential strain response in alkaline phosphatase activity to available phosphorus in Microcoleus autumnalis.

Toxic, benthic cyanobacterial proliferations have increased in frequency and severity globally and can have negative impacts on aquatic ecosystems, recreation and human health. Microcoleus autumnalis has been associated with numerous animal fatalities and is causing increasing concern. It tends to grow in systems with moderate dissolved inorganic nitrogen and very low dissolved reactive phosphorus. Acquisition of nutrients, particularly phosphorus, from organic sources may explain how M. autumnalis can reach the high biomass in these relatively nutrient deplete environments. In the present study the effect of phosphorus concentration and source on alkaline phosphatase activity was investigated in toxic and non-toxic M. autumnalis strains. Toxic strains exhibited significantly higher alkaline phosphatase activity than non-toxic strains (p < 0.05), and alkaline phosphatase activity increased in all strains under phosphorus-depleted conditions (p < 0.05). Alkaline phosphatase activity was also present in environmental M. autumnalis mats, though at lower levels than in laboratory experiments. The presence of alkaline phosphatase activity indicates that the acquisition of phosphorus from organic phosphorus sources may contribute to the ability of M. autumnalis to grow in systems with low dissolved reactive phosphorus.

Multiple cyanotoxin congeners produced by sub-dominant cyanobacterial taxa in riverine cyanobacterial and algal mats.

Benthic cyanobacterial proliferations in rivers are have been reported with increasing frequency worldwide. In the Eel and Russian rivers of California, more than a dozen dog deaths have been attributed to cyanotoxin toxicosis since 2000. Periphyton proliferations in these rivers comprise multiple cyanobacterial taxa capable of cyanotoxin production, hence there is uncertainty regarding which taxa are producing toxins. In this study, periphyton samples dominated by the cyanobacterial genera Anabaena spp. and Microcoleus spp. and the green alga Cladophora glomerata were collected from four sites in the Eel River catchment and one site in the Russian River. Samples were analysed for potential cyanotoxin producers using polymerase chain reaction (PCR) in concert with Sanger sequencing. Cyanotoxin concentrations were measured using liquid chromatography tandem-mass spectrometry, and anatoxin quota (the amount of cyanobacterial anatoxins per toxigenic cell) determined using droplet digital PCR. Sequencing indicated Microcoleus sp. and Nodularia sp. were the putative producers of cyanobacterial anatoxins and nodularins, respectively, regardless of the dominant taxa in the mat. Anatoxin concentrations in the mat samples varied from 0.1 to 18.6 µg g-1 and were significantly different among sites (p < 0.01, Wilcoxon test); however, anatoxin quotas were less variable (< 5-fold). Dihydroanatoxin-a was generally the most abundant variant in samples comprising 38% to 71% of the total anatoxins measured. Mats dominated by the green alga C. glomerata contained both anatoxins and nodularin-R at concentrations similar to those of cyanobacteria-dominated mats. This highlights that even when cyanobacteria are not the dominant taxa in periphyton, these mats may still pose a serious health risk and indicates that more widespread monitoring of all mats in a river are necessary.

Development and Application of a Quantitative PCR Assay to Assess Genotype Dynamics and Anatoxin Content in Microcoleus autumnalis-Dominated Mats.

Microcoleus is a filamentous cyanobacteria genus with a global distribution. Some species form thick, cohesive mats over large areas of the benthos in rivers and lakes. In New Zealand Microcoleus autumnalis is an anatoxin producer and benthic proliferations are occurring in an increasing number of rivers nationwide. Anatoxin content in M. autumnalis-dominated mats varies spatially and temporally, making understanding and managing proliferations difficult. In this study a M. autumnalis-specific TaqMan probe quantitative PCR (qPCR) assay targeting the anaC gene was developed. The assay was assessed against 26 non-M. autumnalis species. The assay had a detection range over seven orders of magnitude, with a limit of detection of 5.14 × 10-8 ng µL-1. The anaC assay and a cyanobacterial specific 16S rRNA qPCR were then used to determine toxic genotype proportions in 122 environmental samples collected from 19 sites on 10 rivers in New Zealand. Anatoxin contents of the samples were determined using LC-MS/MS and anatoxin quota per toxic cell calculated. The percentage of toxic cells ranged from 0 to 30.3%, with significant (p < 0.05) differences among rivers. The anatoxin content in mats had a significant relationship with the percentage of toxic cells (R² = 0.38, p < 0.001), indicating that changes in anatoxin content in M. autumnalis-dominated mats are primarily related to the dominance of toxic strains. When applied to more extensive samples sets the assay will enable new insights into how biotic and abiotic parameters influence genotype composition, and if applied to RNA assist in understanding anatoxin production.