Global Trends of Benthic Bacterial Diversity and Community Composition Along Organic Enrichment Gradients of Salmon Farms

An approach for quantifying the influence of fish farm waste on hard-bottom habitats

Sea-based fish farms release significant quantities of fish feaces, waste feed and associated contaminants into the surrounding environment, necessitating ecological effects-based monitoring. Traditionally, this has been achieved by measuring geochemical and biological properties of sediments obtained by a benthic grab. However, many modern high-capacity farms are situated over hard substrates, rendering conventional sampling methods ineffective and hindering environmental management. Fortunately, new possibilities have emerged with the advent of high-throughput sequencing, and in-particular, the use of microbial eDNA as an indicator of benthic enrichment, and new methods have been developed for obtaining the source material irrespective of substrate. Here, we first demonstrate how applying the default approach used by the world's largest salmon farming industry to hard substrates is essentially qualitative and can lead to inappropriate conclusions. We then show the results from a substrate independent benthic sampling (SIBS) method coupled with a structured quantitative visual assessment of seabed images. The approach utilizes three novel indicators that quantify the degree of visual organic loading (VOL) and ecological effects (VEE) and elucidate the magnitude and spatial extent of the waste influence field via an index (bMBI) derived from 16S eDNA. VOL and VEE were useful for quantifying and classifying effects beneath farms, but were of limited use further away, whereas the bMBI was both highly quantitative and sensitive in the immediate and more distant receiving environments. This dual approach provides a basis for the development of a quantitative hard-bottom monitoring system and should therefore permit effective environment management in the future.

A comparison of visual and molecular methods for inferring biological communities in aquaculture enriched sediments - Impact assessment and cost-benefit analysis

Nutrients introduced to the environment by finfish aquaculture pose environmental risks, which can be mitigated by robust environmental monitoring. Biological communities in soft sediments are good indicators of aquaculture derived environmental changes. Traditionally, monitoring programs have visually surveyed macrofauna communities. However, DNA metabarcoding is a potentially more efficient alternative. We compared alpha diversity, multivariate dispersion and taxonomic composition of macrofauna communities with metabarcoding derived bacterial and eukaryote communities along an organic enrichment gradient at a salmon farm in Tasmania, Australia. Additionally, we conducted a cost-benefit analysis comparing the approaches. All methods identified indicator taxa that changed in abundance over the enrichment gradient. Macrofauna analysis was the most sensitive method for detecting changes in alpha diversity, while metabarcoding was most sensitive for multivariate dispersion. Taxonomic composition of animal communities derived from the two methods differed drastically. Metabarcoding was cheaper than macrofauna for ≥93 samples and quicker for ≥14 samples.

The Core of the Matter-Importance of Identification Method and Biological Replication for Benthic Marine Monitoring

Benthic macrofauna are important and widely used biological indicators of marine ecosystems as they have limited mobility and therefore integrate the effects of local environmental stressors over time. Recently, environmental DNA (eDNA) analysis has provided a potentially more resource-efficient approach for benthic biomonitoring than traditional morphology-based methods. Several studies have compared eDNA with morphology-based monitoring, but few have compared the two approaches using the exact same sediment cores. In addition, the meiofauna and pelagic organisms obtained as 'bycatch' using eDNA have largely been disregarded from comparisons. Here, we address these shortcomings through comparative invertebrate analyses of six sediment sample replicates from each of four stations in Denmark, using eDNA metabarcoding and morphological identification. Our results revealed large variation between the six replicates for both methods and little overlap in taxon compositions between methods. While the morphological dataset was dominated by molluscs and annelids, the eDNA dataset was dominated by arthropods and annelids. Using community composition data, we found that sampling stations could be distinguished both with eDNA and morphology. Finally, we inferred expected total richness from extrapolated accumulation curves of detected taxa from each method. This indicated that eDNA metabarcoding requires less replication than morphology for maximum coverage of diversity to be reached. However, both methods required high levels of replication, and our results on taxonomic composition add to the evidence that morphological and eDNA-based methods should preferably be used as complimentary tools for marine bioassessment.

Descriptive analyses of bacterial communities in marine sediment microcosms spiked with fish wastes, emamectin benzoate, and oxytetracycline

In marine sediments surrounding salmon aquaculture sites, organic matter (OM) enrichment has been shown to influence resident bacterial community composition; however, additional effects on these communities due to combined use of the sea-lice therapeutant emamectin benzoate (EMB) and the widely used antibiotic oxytetracycline (OTC) are unknown. Here, we use sediment microcosms to assess the influence of OM, EMB, and OTC on benthic bacterial communities. Microcosms consisted of mud or sand sediments enriched with OM (fish and feed wastes) and spiked with EMB and OTC at environmentally-relevant concentrations. Samples were collected from initial matrices at the initiation of the trial and after 110 days for 16 S rRNA gene sequencing of the V3-V4 region and microbiome profiling. The addition of OM in both mud and sand sediments reduced alpha diversities; for example, an average of 1106 amplicon sequence variants (ASVs) were detected in mud with no OM addition, while only 729 and 596 ASVs were detected in mud with low OM and high OM, respectively. Sediments enriched with OM had higher relative abundances of Spirochaetota, Firmicutes, and Bacteroidota. For instance, Spirochaetota were detected in sediments with no OM with a relative abundance range of 0.01-1.2%, while in sediments enriched with OM relative abundance varied from 0.16% to 26.1%. In contrast, the addition of EMB (60 ng/g) or OTC (150 ng/g) did not result in distinct taxonomic shifts in the bacterial communities compared to un-spiked sediments during the timeline of this experiment. EMB and OTC concentrations may have been below effective inhibitor concentrations for taxa in these communities; further work should explore gene content and the presence of antibiotic resistance genes (ARGs) in sediment-dwelling bacteria.

Simultaneous analysis of seven 16S rRNA hypervariable gene regions increases efficiency in marine bacterial diversity detection

Environmental DNA sequencing is the gold standard to reveal microbial community structures. In most applications, a one-fragment PCR approach is applied to amplify a taxonomic marker gene, usually a hypervariable region of the 16S rRNA gene. We used a new reverse complement (RC)-PCR-based assay that amplifies seven out of the nine hypervariable regions of the 16S rRNA gene, to interrogate bacterial communities in sediment samples collected from different coastal marine sites with an impact gradient. In parallel, we employed a traditional one-fragment analysis of the hypervariable V3-V4 region to investigate whether the RC-PCR reveals more of the 'unseen' diversity obtained by the one-fragment approach. As a benchmark for the full deck of diversity, we subjected the samples to PCR-free metagenomic sequencing. None of the two PCR-based approaches recorded the full taxonomic repertoire obtained from the metagenomics datasets. However, the RC-PCR approach detected 2.8 times more bacterial genera compared to the near-saturation sequenced V3-V4 samples. RC-PCR is an ideal compromise between the standard one-fragment approach and metagenomics sequencing and may guide future environmental sequencing studies, in which bacterial diversity is a central subject.

Host Species and Environment Shape the Gut Microbiota of Cohabiting Marine Bivalves

Pacific oysters (Crassostrea gigas) and Mediterranean mussels (Mytilus galloprovincialis) are commercially important marine bivalves that frequently coexist and have overlapping feeding ecologies. Like other invertebrates, their gut microbiota is thought to play an important role in supporting their health and nutrition. Yet, little is known regarding the role of the host and environment in driving these communities. Here, bacterial assemblages were surveyed from seawater and gut aspirates of farmed C. gigas and co-occurring wild M. galloprovincialis in summer and winter using Illumina 16S rRNA gene sequencing. Unlike seawater, which was dominated by Pseudomonadata, bivalve samples largely consisted of Mycoplasmatota (Mollicutes) and accounted for >50% of the total OTU abundance. Despite large numbers of common (core) bacterial taxa, bivalve-specific species (OTUs) were also evident and predominantly associated with Mycoplasmataceae (notably Mycoplasma). An increase in diversity (though with varied taxonomic evenness) was observed in winter for both bivalves and was associated with changes in the abundance of core and bivalve-specific taxa, including several representing host-associated and environmental (free-living or particle-diet associated) organisms. Our findings highlight the contribution of the environment and the host in defining the composition of the gut microbiota in cohabiting, intergeneric bivalve populations.

Comparing quantile regression spline analyses and supervised machine learning for environmental quality assessment at coastal marine aquaculture installations

Organic enrichment associated with marine finfish aquaculture is a local stressor of marine coastal ecosystems. To maintain ecosystem services, the implementation of biomonitoring programs focusing on benthic diversity is required. Traditionally, impact-indices are determined by extracting and identifying benthic macroinvertebrates from samples. However, this is a time-consuming and expensive method with low upscaling potential. A more rapid, inexpensive, and robust method to infer the environmental quality of marine environments is eDNA metabarcoding of bacterial communities. To infer the environmental quality of coastal habitats from metabarcoding data, two taxonomy-free approaches have been successfully applied for different geographical regions and monitoring goals, namely quantile regression splines (QRS) and supervised machine learning (SML). However, their comparative performance remains untested for monitoring the impact of organic enrichment introduced by aquaculture on marine coastal environments. We compared the performance of QRS and SML using bacterial metabarcoding data to infer the environmental quality of 230 aquaculture samples collected from seven farms in Norway and seven farms in Scotland along an organic enrichment gradient. As a measure of environmental quality, we used the Infaunal Quality Index (IQI) calculated from benthic macrofauna data (reference index). The QRS analysis plotted the abundance of amplicon sequence variants (ASVs) as a function to the IQI from which the ASVs with a defined abundance peak were assigned to eco-groups and a molecular IQI was subsequently calculated. In contrast, the SML approach built a random forest model to directly predict the macrofauna-based IQI. Our results show that both QRS and SML perform well in inferring the environmental quality with 89% and 90% accuracy, respectively. For both geographic regions, there was high correspondence between the reference IQI and both the inferred molecular IQIs (p < 0.001), with the SML model showing a higher coefficient of determination compared to QRS. Among the 20 most important ASVs identified by the SML approach, 15 were congruent with the good quality spline ASV indicators identified via QRS for both Norwegian and Scottish salmon farms. More research on the response of the ASVs to organic enrichment and the co-influence of other environmental parameters is necessary to eventually select the most powerful stressor-specific indicators. Even though both approaches are promising to infer environmental quality based on metabarcoding data, SML showed to be more powerful in handling the natural variability. For the improvement of the SML model, addition of new samples is still required, as background noise introduced by high spatio-temporal variability can be reduced. Overall, we recommend the development of a powerful SML approach that will be onwards applied for monitoring the impact of aquaculture on marine ecosystems based on eDNA metabarcoding data.

Maximising signal-to-noise ratios in environmental DNA-based monitoring

Man's impacts on global ecosystems are increasing and there is a growing demand that these activities be appropriately monitored. Monitoring requires measurement of a response metric ('signal') that changes maximally and consistently in response to the monitored activity irrespective of other factors ('noise'), thus maximising the signal-to-noise ratio. Indices derived from time-consuming morphology-based taxonomic identification of organisms are a core part of many monitoring programmes. Metabarcoding is an alternative to morphology-based identification and involves the sequencing of short fragments of DNA ('markers') from multiple taxa simultaneously. DNA suitable for metabarcoding includes that extracted from environmental samples (eDNA). Metabarcoding outputs DNA sequences that can be identified (annotated) by matching them against archived annotated sequences. However, sequences from most organisms are not archived - preventing annotation and potentially limiting metabarcoding in monitoring applications. Consequently, there is growing interest in using unannotated sequences as response metrics in monitoring programmes. We compared the sequences from three commonly used markers (16S (V3/V4 regions), 18S (V1/V2 regions) and COI) and, sampling along steep impact gradients, showed that the 16S and COI sequences were associated with the largest and smallest signal-to-noise ratio respectively. We trialled four separate, intuitive, noise-reduction approaches and demonstrated that removing less frequent sequences improved the signal-to-noise ratio, partitioning an additional 25 % from noise to explanatory factors in non-parametric ANOVA (NPA) and reducing dispersion in the data. For the 16S marker, retaining only the most frequently observed sequence, per sample, resulting in nine sequences across 150 samples, generated a near-maximal signal-to-noise ratio (95 % of the variance explained in NPA). We recommend that NPA, combined with rigorous elimination of less frequent sequences, be used to pre-filter sequences/taxa being used in monitoring applications. Our approach will simplify downstream analysis, for example the identification of key taxa and functional associations.

Environmental DNA-based profiling of benthic bacterial and eukaryote communities along a crude oil spill gradient in a coral reef in the Persian Gulf

Coral reef ecosystems in the Persian Gulf are frequently exposed to crude oil spills. We investigated benthic bacterial and eukaryote community structures at such coral reef sites subjected to different degrees of polycyclic aromatic hydrocarbon (PAH) pollution using environmental DNA (eDNA) metabarcoding. Both bacterial and eukaryote communities responded with pronounced shifts to crude oil pollution and distinguished control sites, moderately and heavily impacted sites with significant confidentiality. The observed community patterns were predominantly driven by Alphaproteobacteria and metazoans. Among these, we identified individual genera that were previously linked to oil spill stress, but also taxa, for which a link to hydrocarbon still remains to be established. Considering the lack of an early-warning system for the environmental status of coral reef ecosystems exposed to frequent crude-oil spills, our results encourage further research towards the development of an eDNA-based biomonitoring tool that exploits benthic bacterial and eukaryote communities as bioindicators.

Composition and functionality of bacterioplankton communities in marine coastal zones adjacent to finfish aquaculture

Finfish aquaculture is a fast-growing primary industry and is increasingly common in coastal ecosystems. Bacterioplankton is ubiquitous in marine environment and respond rapidly to environmental changes. Changes in bacterioplankton community are not well understood in semi-enclosed stratified embayments. This study aims to examine aquaculture effects in the composition and functional profiles of the bacterioplankton community using amplicon sequencing along a distance gradient from two finfish leases in a marine embayment. Results revealed natural stratification in bacterioplankton associated to NOx, conductivity, salinity, temperature and PO₄. Among the differentially abundant bacteria in leases, we found members associated with nutrient enrichment and aquaculture activities. Abundant predicted functions near leases were assigned to organic matter degradation, fermentation, and antibiotic resistance. This study provides a first effort to describe changes in the bacterioplankton community composition and function due to finfish aquaculture in a semi-enclosed and highly stratified embayment with a significant freshwater input.

Comparing sediment preservation methods for genomic biomonitoring of coastal marine ecosystems

To avoid loss of genetic information in environmental DNA (eDNA) field samples, the preservation of nucleic acids during field sampling is a critical step. In the development of standard operating procedures (SOPs) for eDNA-based compliance monitoring, the effect of different routinely used sediment preservations on biological community structures serving as bioindicators has gone untested. We compared eDNA metabarcoding results of marine bacterial communities from sample aliquots that were treated with a nucleic acid preservation solution (treated samples) and aliquots that were frozen without further treatment (non-treated samples). Sediment samples were obtained from coastal locations subjected to different stressors (aquaculture, urbanization, industry). DNA extraction efficiency, bacterial community profiles, and measures of alpha- and beta-diversity were highly congruent between treated and non-treated samples. As both preservation methods provide the same relevant information to environmental managers and regulators, we recommend the inclusion of both methods into SOPs for biomonitoring in marine coastal environments.