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.

A taxonomy-free diatom eDNA-based technique for assessing lake trophic level using lake sediments.

Anthropogenic eutrophication is one of the most pressing issues facing lakes globally. Our ability to manage lake eutrophication is hampered by the limited spatial and temporal extents of monitoring records, stemming from the time-consuming and expensive nature of physiochemical and biological monitoring. Diatom-based biomonitoring presents an alternative to traditional eutrophication monitoring, yet it is restricted by the high degree of taxonomic expertise required. Environmental DNA metabarcoding, while providing a promising substitute for diatom community enumeration, is plagued by inadequate taxonomic coverage of reference databases and methodological bias, limiting its use for biomonitoring. Here we show that taxonomy-free diatom-biomonitoring, in which environmental DNA metabarcoding data is utilised but not assigned to specific taxonomic classes, presents an accurate, fast, and relatively automated alternative to taxonomically assigned eutrophication biomonitoring. Our taxonomy-free index accounted for 85% of trophic level variability across 89 lakes and had the lowest average prediction error of the three approaches tested. By not relying on taxonomic identification or metabarcoding reference databases, taxonomy-free biomonitoring maintains diatom diversity that is lost in taxonomic assignment using molecular approaches. Furthermore, by utilising lake sediments, the approach outlined here presents a time-integrated estimation of lake trophic level and thus does not require time-consuming seasonal sampling. Taxonomy-free biomonitoring addresses the limitations of traditional physicochemical eutrophication monitoring and taxonomic biomonitoring alternatives and can be used to extend the spatial and temporal extents of eutrophication monitoring.

Resolving 500 Years of Anthropogenic Impacts in a Mesotrophic Lake: Nutrients Outweigh Other Drivers of Lake Change.

Interactions among multiple stressors, legacies of past perturbations, and the lack of historical information make it difficult to determine the influence of individual anthropogenic impacts on lakes and separate them from natural ecosystem variability. In the present study, we coupled paleolimnological approaches, historical data, and ecological experiments to disentangle the impacts of multiple long-term stressors on lake ecosystem structure and function. We found that the lake structure and function remained resistant to the impacts of catchment deforestation and erosion, and the introduction of several exotic fish species. Changes in ecosystem structure and function were consistent, with nutrient enrichment being the primary driver of change. Significant and sustained changes in the lake diatom community structure (and their nutrient requirements), bacterial community function, and paleolimnological proxies of ecosystem function coincided with nitrogen and phosphorus fertilizers in the catchment. The results highlight that the effects of increased nutrient inputs are much stronger than the influence of other, potentially significant, drivers of ecosystem change, and that the degree of nutrient impact can be underestimated by environmental monitoring due to its diffuse and accumulative nature. Delineating the effects of multiple anthropogenic drivers requires long-term records of both impacts and lake ecosystem change across multiple trophic levels.

Lake microbial communities are not resistant or resilient to repeated large-scale natural pulse disturbances.

Opportunities to study community-level responses to extreme natural pulse disturbances in unaltered ecosystems are rare. Lake sediment records that span thousands of years can contain well-resolved sediment pulses, triggered by earthquakes. These palaeorecords provide a means to study repeated pulse disturbances and processes of resistance (insensitivity to disturbance) and ecological resilience (capacity to regain structure, function and process). In this study, sedimentary DNA was extracted from a sediment core from Lake Paringa (New Zealand) that is situated in a near natural catchment. Metabarcoding and inferred functions were used to assess the lake microbial community over the past 1100 years - a period that included four major earthquakes. Microbial community composition and function differed significantly between highly perturbed (postseismic, ~50 years) phases directly after the earthquakes and more stable (interseismic, ~250 years) phases, indicating a lack of community resistance. Although community structure differed significantly in successive postseismic phases, function did not, suggesting potential functional redundancy. Significant differences in composition and function in successive interseismic phases demonstrate that communities are not resilient to large-scale natural pulse disturbances. The clear difference in structure and function, and high number of indicator taxa (responsible for driving differences in communities between phases) in the fourth interseismic phase probably represents a regime shift, possibly due to the two-fold increase in sediment and terrestrial biospheric organic carbon fluxes recorded following the fourth earthquake. Large pulse disturbances that enhance sediment inputs into lake systems may produce an underappreciated mechanism that destabilises lake ecosystem processes.

Blue Waters, Green Bottoms: Benthic Filamentous Algal Blooms Are an Emerging Threat to Clear Lakes Worldwide.

Nearshore (littoral) habitats of clear lakes with high water quality are increasingly experiencing unexplained proliferations of filamentous algae that grow on submerged surfaces. These filamentous algal blooms (FABs) are sometimes associated with nutrient pollution in groundwater, but complex changes in climate, nutrient transport, lake hydrodynamics, and food web structure may also facilitate this emerging threat to clear lakes. A coordinated effort among members of the public, managers, and scientists is needed to document the occurrence of FABs, to standardize methods for measuring their severity, to adapt existing data collection networks to include nearshore habitats, and to mitigate and reverse this profound structural change in lake ecosystems. Current models of lake eutrophication do not explain this littoral greening. However, a cohesive response to it is essential for protecting some of the world's most valued lakes and the flora, fauna, and ecosystem services they sustain.

A Microencapsulation Method for Delivering Tetrodotoxin to Bivalves to Investigate Uptake and Accumulation.

Most marine biotoxins are produced by microalgae. The neurotoxin tetrodotoxin (TTX) has been reported in many seafood species worldwide but its source is unknown, making accumulation and depuration studies in shellfish difficult. Tetrodotoxin is a water-soluble toxin and cannot be directly ingested by shellfish. In the present study, a method was developed which involved binding TTX to solid particles of humic acid and encapsulating them in agar-gelatin capsules. A controlled quantity of TTX-containing microcapsules (size range 20-280 µm) was fed to Paphies australis, a bivalve known to accumulate TTX in the wild. The TTX-containing microcapsules were fed to P. australis every second day for 13 days. Ten P. australis (including five controls fed non-toxic microalgae) were harvested after 7 days and ten after 13 days. Paphies australis accumulated TTX, reaching concentrations of up to 103 µg kg-1 by day 13, exceeding the European Food Safety Authority recommended concentration of 44 µg kg-1 in shellfish. This novel method will allow future studies to explore the effects, accumulation and depuration rates of TTX in different animals and document how it is transferred through food webs.

Elucidating Biodiversity Shifts in Ballast Water Tanks during a Cross-Latitudinal Transfer: Complementary Insights from Molecular Analyses.

In this study, the evolution of ballast water (BW) assemblages across different trophic levels was characterized over a 21 day cross-latitudinal vessel transit using a combination of molecular methods. Triplicate BW samples were collected every second day and size-fractionated (<2.7, 10, >50 µm). Measurements of adenosine triphosphate (ATP) and metabarcoding of environmental nucleic acid (DNA and RNA) analyses, complemented by microscopy and flow cytometry, were performed on each sample. Measured ATP concentrations exhibited high variance between replicates and a strong negative trend in the large (≥50 µm) fraction over the voyage. In concert with microscopy, the metabarcoding data indicated a die-off of larger metazoans during the first week of study and gradual reductions in dinoflagellates and ochrophytes. The ATP and metabarcoding data signaled persistent or increased cellular activity of heterotrophic bacteria and protists in the BW, which was supported by flow cytometry. The metabarcoding showed the presence of active bacteria in all size fractions, suggesting that the sequential filtration approach does not ensure taxonomical differentiation, which has implications for BW quality assessment. Although our data show that ATP and metabarcoding have potential for indicative BW screening for BW compliance monitoring, further research and technological development is needed to improve representativeness of sampling and deliver the unequivocal response criteria required by the international Ballast Water Management Convention.

Considerations for incorporating real-time PCR assays into routine marine biosecurity surveillance programmes: a case study targeting the Mediterranean fanworm (Sabella spallanzanii) and club tunicate (Styela clava) 1.

Molecular techniques may provide effective tools to enhance marine biosecurity surveillance. Prior to routine implementation, evidence-based consideration of their benefits and limitations is needed. In this study, we assessed the efficiency and practicality of visual diver surveys and real-time PCR assays (targeting DNA and RNA) for detecting two marine invasive species whose infestation levels varied between species and location: Sabella spallanzanii and Styela clava. Filtered water samples (n = 171) were collected in parallel with dive surveys at two locations as part of the New Zealand Marine High Risk Site Surveillance programme: Nelson Harbour (27 sites) and Waitemata Harbour (30 sites). Diver surveys resulted in a greater number of detections compared to real-time PCR: S. clava - 21 versus 5 sites in Nelson, 6 versus 1 in Auckland; S. spallanzanii - 18 versus 10 in Auckland, no detections in Nelson. Occupancy modelling derived detection probabilities for the real-time PCR for S. clava were low (14%), compared to S. spallanzanii (66%). This could be related to abundances, or species-specific differences in DNA shedding. Only one RNA sample was positive, suggesting that most detections were from extracellular DNA or non-viable fragments. While molecular methods cannot yet replace visual observations, this study shows they provide useful complementary information.

Metabarcoding and metabolomics offer complementarity in deciphering marine eukaryotic biofouling community shifts.

Metabarcoding and metabolomics were used to explore the taxonomic composition and functional diversity of eukaryotic biofouling communities on plates with antifouling paints at two French coastal sites: Lorient (North Eastern Atlantic Ocean; temperate and eutrophic) and Toulon (North-Western Mediterranean Sea; mesotrophic but highly contaminated). Four distinct coatings were tested at each site and season for one month. Metabarcoding showed biocidal coatings had less impact on eukaryotic assemblages compared to spatial and temporal effects. Ciliophora, Chlorophyceae or Cnidaria (mainly hydrozoans) were abundant at Lorient, whereas Arthropoda (especially crustaceans), Nematoda, and Ochrophyta dominated less diversified assemblages at Toulon. Seasonal shifts were observed at Lorient, but not Toulon. Metabolomics also showed clear site discrimination, but these were associated with a coating and not season dependent clustering. The meta-omics analysis enabled identifications of some associative patterns between metabolomic profiles and specific taxa, in particular those colonizing the plates with biocidal coatings at Lorient.

The Abundance of Toxic Genotypes Is a Key Contributor to Anatoxin Variability in Phormidium-Dominated Benthic Mats.

The prevalence of benthic proliferations of the anatoxin-producing cyanobacterium Phormidium are increasing in cobble-bed rivers worldwide. Studies to date have shown high spatial and temporal variability in anatoxin concentrations among mats. In this study we determined anatoxin quotas (toxins per cell) in field samples and compared these results to the conventionally-used concentrations (assessed per dry weight of mat). Three mats were selected at sites in two rivers and were sampled every 2-3 h for 24-26 h. The samples were lyophilized and ground to a fine homogenous powder. Two aliquots of known weights were analyzed for anatoxin congeners using liquid chromatography-mass spectrometry, or digital droplet PCR with Phormidium-specific anaC primers to measure absolute quantities of gene copies. Anatoxin concentrations in the mats varied 59- and 303-fold in the two rivers over the study periods. A similar pattern was observed among gene copies (53- and 2828-fold). When converted to anatoxin quotas there was markedly less variability (42- and 16-fold), but significantly higher anatoxin quotas were observed in mats from the second river (p < 0.001, Student's t-test). There were no obvious temporal patterns with high and low anatoxin concentrations or quotas measured at each sampling time and across the study period. These results demonstrate that variability in anatoxin concentrations among mats is primarily due to the abundance of toxic genotypes. No consistent modulation in anatoxin production was observed during the study, although significant differences in anatoxin quotas among rivers suggest that site-specific physiochemical or biological factors may influence anatoxin production.

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.

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.

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.

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.

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.

High levels of structural diversity observed in microcystins from Microcystis CAWBG11 and characterization of six new microcystin congeners.

Microcystins (MCs) are cyclic peptides produced by cyanobacteria, which can be harmful to humans and animals when ingested. Differences in the coding of the non­ribosomal peptide synthetase/polyketide synthase enzyme complex responsible for microcystin production have resulted in more than 100 microcystin variants being reported to date. The microcystin diversity of Microcystis CAWBG11 was investigated using matrix-assisted laser desorption/ionization-time of flight mass spectrometry and liquid chromatography-mass spectrometry. This revealed that CAWBG11 simultaneously produced 21 known microcystins and six new congeners: [Asp3] MC-RA, [Asp3] MC-RAba, [Asp3] MC-FA, [Asp3] MC-WA, MC-FAba and MC-FL. The new congeners were putatively characterized by tandem mass spectrometry and chemical derivatization. A survey of the microcystin congeners produced by 49 cyanobacterial strains documented in scientific literature showed that cyanobacteria generally produce four microcystin congeners, but strains which produce up to 47 microcystin congeners have been reported. Microcystis CAWBG11 (which produces at least 27 congeners) was positioned in the top ten percentile of the strains surveyed, and showed fluidity of the amino acids incorporated into both position two and position four.

First detection of tetrodotoxin in the bivalve Paphies australis by liquid chromatography coupled to triple quadrupole mass spectrometry with and without precolumn reaction.

Two methods for the determination of tetrodotoxin (TTX) in marine biota have been developed and validated using ultra-performance LC coupled to triple quadrupole MS. The direct analysis of TTX is completed in one method, while the other method detects the dehydration product of TTX after reaction with base. The methods were validated in a single-laboratory trial and used to test Paphies australis (pipi) samples collected from Whangapoua, New Zealand during April 2011. Pa. australis is a commonly eaten species of bivalve that was found to contain TTX at levels up to 0.80 mg/kg in this study. The methods exhibited recoveries ranging from 94 to 120%, and the within laboratory reproducibility ranged from 6 to 27% for Pleurobranchaea maculata (grey-side gilled sea slug) and bivalve matrixes. Use of the method using a dehydration step showed no evidence of TTX analogs in any of the samples.

An integrative approach to assess the impact of disturbance on native fish in lakes.

Freshwater fish are in a perilous state with more than 30% of species considered critically endangered. Yet significant ecological and methodological complexities constrain our ability to determine how disturbances are impacting native fish communities. We review current methods used to assess the responses of fish communities, especially native fish, to disturbances, with a focus on lakes. These methods include contemporary population surveys, manipulative experimental approaches, paleolimnological approaches and Indigenous Knowledge and social histories. We identify knowledge gaps, such as a lack of baseline data for native fish, an inability to assess the impact of historical disturbances, stressor response dynamics in contemporary multi-stressor environments, and natural disturbance regimes. Our assessment of the current methods highlights challenges to filling these knowledge gaps using the reviewed methods. We advocate strongly for the implementation of an integrative approach that combines emerging technologies (i.e. molecular-based techniques in contemporary surveys and paleolimnology) and underutilised knowledge streams (i.e. Indigenous Knowledge and social histories) which should be used in concert with conventional methods. This integrative approach will allow researchers to determine the key drivers of decline and the degree of change, which will enable more informed and successful management actions.

Sedimentary ancient DNA reveals the impact of anthropogenic land use disturbance and ecological shifts on fish community structure in small lowland lake.

Freshwater fish biodiversity and abundance are decreasing globally. The drivers of decline are primarily anthropogenic; however, the causative links between disturbances and fish community change are complex and challenging to investigate. We used a suite of sedimentary DNA methods (droplet digital PCR and metabarcoding) and traditional paleolimnological approaches, including pollen and trace metal analysis, ITRAX X-ray fluorescence and hyperspectral core scanning to explore changes in fish abundance and drivers over 1390 years in a small lake. This period captured a disturbance trajectory from pre-human settlement through subsistence living to intensive agriculture. Generalized additive mixed models explored the relationships between catchment inputs, internal drivers, and fish community structure. Fish community composition distinctly shifted around 1350 CE, with the decline of a sensitive Galaxias species concomitant with early land use changes. Total fish abundance significantly declined around 1950 CE related to increases in ruminant bacterial DNA (a proxy for ruminant abundance) and cadmium flux (a proxy for phosphate fertilizers), implicating land use intensification as a key driver. Concurrent shifts in phytoplankton and zooplankton suggested that fish communities were likely impacted by food web dynamics. This study highlights the potential of sedDNA to elucidate the long-term disturbance impacts on biological communities in lakes.

A paleoecological investigation of recent cyanobacterial blooms and their drivers in two contrasting lakes.

Cyanobacterial blooms are one of the most significant threats to global water security and freshwater biodiversity. Interactions among multiple stressors, including habitat degradation, species invasions, increased nutrient runoff, and climate change, are key drivers. However, assessing the role of anthropogenic activity on the onset of cyanobacterial blooms and exploring response variation amongst lakes of varying size and depth is usually limited by lack of historical records. In the present study we applied molecular, paleolimnological (trace metal, Itrax-µ-XRF and hyperspectral scanning, chronology), paleobotanical (pollen) and historical data to reconstruct cyanobacterial abundance and community composition and anthropogenic impacts in two dune lakes over a period of up to 1200 years. Metabarcoding and droplet digital PCR results showed very low levels of picocyanobacteria present in the lakes prior to about CE 1854 (1839-1870 CE) in the smaller shallow Lake Alice and CE 1970 (1963-1875 CE) in the larger deeper Lake Wiritoa. Hereafter bloom-forming cyanobacteria were detected and increased notably in abundance post CE 1984 (1982-1985 CE) in Lake Alice and CE 1997 (1990-2007 CE) in Lake Wiritoa. Currently, the magnitude of blooms is more pronounced in Lake Wiritoa, potentially attributable to hypoxia-induced release of phosphorus from sediment, introducing an additional source of nutrients. Generalized linear modelling was used to investigate the contribution of nutrients (proxy = bacterial functions), temperature, redox conditions (Mn:Fe), and erosion (Ti:Inc) in driving the abundance of cyanobacteria (ddPCR). In Lake Alice nutrients and erosion had a statistically significant effect, while in Lake Wiritoa nutrients and redox conditions were significant.