Swarm and UNOISE outperform DADA2 and Deblur for denoising high-diversity marine seafloor samples.

The performance of sequence variant resolution analytic tools for metabarcoding has not yet been adequately benchmarked for high-diversity environmental samples. We therefore evaluated the sequence variant tools DADA2, Deblur, Swarm, and UNOISE, using high-diversity seafloor samples, resulting in comparisons of 1800 sequence variant tables. The evaluation was based on 30 sediment grab samples, for which 3 replica samples were collected. Each replica sample was extracted using 5 common DNA extraction kits, resulting in 450 DNA extracts which were 16S rRNA gene sequenced (V3-V4), using Illumina. Assessments included variation across replica samples, extraction kits, and denoising methods, in addition to applying prior knowledge about alpha diversity correlations toward the cosmopolitan marine archaeon Nitrosopumilus with high diversity and the sulfide oxidizing Sulfurovum with low diversity. DADA2 displayed the highest variance between replicates (Manhattan distance 1.14), while Swarm showed the lowest variance (Manhattan distance 0.93). For the analysis based on prior biological knowledge, UNOISE displayed the highest alpha diversity (Simpson's D) correlation toward Nitrosopumilus (Spearman rho = 0.85), while DADA2 showed the lowest (Spearman rho = 0.10). Deblur completely eliminated Nitrosopumilus from the dataset. For Sulfurovum, on the other hand, all the methods showed comparable results. In conclusion, our evaluations show that Swarm and UNOISE performed better than DADA2 and Deblur for high-diversity seafloor samples.

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

The analysis of benthic bacterial community structure has emerged as a powerful alternative to traditional microscopy-based taxonomic approaches to monitor aquaculture disturbance in coastal environments. However, local bacterial diversity and community composition vary with season, biogeographic region, hydrology, sediment texture, and aquafarm-specific parameters. Therefore, without an understanding of the inherent variation contained within community complexes, bacterial diversity surveys conducted at individual farms, countries, or specific seasons may not be able to infer global universal pictures of bacterial community diversity and composition at different degrees of aquaculture disturbance. We have analyzed environmental DNA (eDNA) metabarcodes (V3-V4 region of the hypervariable SSU rRNA gene) of 138 samples of different farms located in different major salmon-producing countries. For these samples, we identified universal bacterial core taxa that indicate high, moderate, and low aquaculture impact, regardless of sampling season, sampled country, seafloor substrate type, or local farming and environmental conditions. We also discuss bacterial taxon groups that are specific for individual local conditions. We then link the metabolic properties of the identified bacterial taxon groups to benthic processes, which provides a better understanding of universal benthic ecosystem function(ing) of coastal aquaculture sites. Our results may further guide the continuing development of a practical and generic bacterial eDNA-based environmental monitoring approach.

Beyond taxonomy: Validating functional inference approaches in the context of fish-farm impact assessments.

Characterization of microbial assemblages via environmental DNA metabarcoding is increasingly being used in routine monitoring programs due to its sensitivity and cost-effectiveness. Several programs have recently been developed which infer functional profiles from 16S rRNA gene data using hidden-state prediction (HSP) algorithms. These might offer an economic and scalable alternative to shotgun metagenomics. To date, HSP-based methods have seen limited use for benthic marine surveys and their performance in these environments remains unevaluated. In this study, 16S rRNA metabarcoding was applied to sediment samples collected at 0 and ≥1,200 m from Norwegian salmon farms, and three metabolic inference approaches (Paprica, Picrust2 and Tax4Fun2) evaluated against metagenomics and environmental data. While metabarcoding and metagenomics recovered a comparable functional diversity, the taxonomic composition differed between approaches, with genera richness up to 20× higher for metabarcoding. Comparisons between the sensitivity (highest true positive rates) and specificity (lowest true negative rates) of HSP-based programs in detecting functions found in metagenomic data ranged from 0.52 and 0.60 to 0.76 and 0.79, respectively. However, little correlation was observed between the relative abundance of their specific functions. Functional beta-diversity of HSP-based data was strongly associated with that of metagenomics (r ≥ 0.86 for Paprica and Tax4Fun2) and responded similarly to the impact of fish farm activities. Our results demonstrate that although HSP-based metabarcoding approaches provide a slightly different functional profile than metagenomics, partly due to recovering a distinct community, they represent a cost-effective and valuable tool for characterizing and assessing the effects of fish farming on benthic ecosystems.

Ecosystems monitoring powered by environmental genomics: A review of current strategies with an implementation roadmap.

A decade after environmental scientists integrated high-throughput sequencing technologies in their toolbox, the genomics-based monitoring of anthropogenic impacts on the biodiversity and functioning of ecosystems is yet to be implemented by regulatory frameworks. Despite the broadly acknowledged potential of environmental genomics to this end, technical limitations and conceptual issues still stand in the way of its broad application by end-users. In addition, the multiplicity of potential implementation strategies may contribute to a perception that the routine application of this methodology is premature or "in development", hence restraining regulators from binding these tools into legal frameworks. Here, we review recent implementations of environmental genomics-based methods, applied to the biomonitoring of ecosystems. By taking a general overview, without narrowing our perspective to particular habitats or groups of organisms, this paper aims to compare, review and discuss the strengths and limitations of four general implementation strategies of environmental genomics for monitoring: (a) Taxonomy-based analyses focused on identification of known bioindicators or described taxa; (b) De novo bioindicator analyses; (c) Structural community metrics including inferred ecological networks; and (d) Functional community metrics (metagenomics or metatranscriptomics). We emphasise the utility of the three latter strategies to integrate meiofauna and microorganisms that are not traditionally utilised in biomonitoring because of difficult taxonomic identification. Finally, we propose a roadmap for the implementation of environmental genomics into routine monitoring programmes that leverage recent analytical advancements, while pointing out current limitations and future research needs.

Comparing sediment DNA extraction methods for assessing organic enrichment associated with marine aquaculture.

Marine sediments contain a high diversity of micro- and macro-organisms which are important in the functioning of biogeochemical cycles. Traditionally, anthropogenic perturbation has been investigated by identifying macro-organism responses along gradients. Environmental DNA (eDNA) analyses have recently been advocated as a rapid and cost-effective approach to measuring ecological impacts and efforts are underway to incorporate eDNA tools into monitoring. Before these methods can replace or complement existing methods, robustness and repeatability of each analytical step has to be demonstrated. One area that requires further investigation is the selection of sediment DNA extraction method. Environmental DNA sediment samples were obtained along a disturbance gradient adjacent to a Chinook (Oncorhynchus tshawytscha) salmon farm in Otanerau Bay, New Zealand. DNA was extracted using four extraction kits (Qiagen DNeasy PowerSoil, Qiagen DNeasy PowerSoil Pro, Qiagen RNeasy PowerSoil Total RNA/DNA extraction/elution and Favorgen FavorPrep Soil DNA Isolation Midi Kit) and three sediment volumes (0.25, 2, and 5 g). Prokaryotic and eukaryotic communities were amplified using primers targeting the 16S and 18S ribosomal RNA genes, respectively, and were sequenced on an Illumina MiSeq. Diversity and community composition estimates were obtained from each extraction kit, as well as their relative performance in established metabarcoding biotic indices. Differences were observed in the quality and quantity of the extracted DNA amongst kits with the two Qiagen DNeasy PowerSoil kits performing best. Significant differences were observed in both prokaryotes and eukaryotes (p < 0.001) richness among kits. A small proportion of amplicon sequence variants (ASVs) were shared amongst the kits (~3%) although these shared ASVs accounted for the majority of sequence reads (prokaryotes: 59.9%, eukaryotes: 67.2%). Differences were observed in the richness and relative abundance of taxonomic classes revealed with each kit. Multivariate analysis showed that there was a significant interaction between "distance" from the farm and "kit" in explaining the composition of the communities, with the distance from the farm being a stronger determinant of community composition. Comparison of the kits against the bacterial and eukaryotic metabarcoding biotic index suggested that all kits showed similar patterns along the environmental gradient. Overall, we advocate for the use of Qiagen DNeasy PowerSoil kits for use when characterizing prokaryotic and eukaryotic eDNA from marine farm sediments. We base this conclusion on the higher DNA quality values and richness achieved with these kits compared to the other kits/amounts investigated in this study. The additional advantage of the PowerSoil Kits is that DNA extractions can be performed using an extractor robot, offering additional standardization and reproducibility of results.

Translational Molecular Ecology in practice: Linking DNA-based methods to actionable marine environmental management.

Molecular-based approaches can provide timely biodiversity assessments, showing an immense potential to facilitate decision-making in marine environmental management. However, the uptake of molecular data into environmental policy remains minimal. Here, we showcase a selection of local to global scale studies applying molecular-based methodologies for environmental management at various stages of implementation. Drawing upon lessons learned from these case-studies, we provide a roadmap to facilitate applications of DNA-based methods to marine policies and to overcome the existing challenges. The main impediment identified is the need for standardized protocols to guarantee data comparison across spatial and temporal scales. Adoption of Translational Molecular Ecology - the sustained collaboration between molecular ecologists and stakeholders, will enhance consensus with regards to the objectives, methods, and outcomes of environmental management projects. Establishing a sustained dialogue among stakeholders is key to accelerating the adoption of molecular-based approaches for marine monitoring and assessment.

Mixed-habitat assimilation of organic waste in coastal environments - It's all about synergy!

Fish farms are increasingly situated in strong current sites above or near to mixed-bottom habitats that include organisms not normally considered in the context of organic enrichment. This study takes a holistic view of the benthic enrichment process by combining different survey techniques on complimentary spatial scales: conventional macrofaunal cores, larger-scale visual quantification of epibiota and environmental-DNA metabarcoding of microbial communities. A large tube forming polychaete (Arenicola marina), normally found intertidally and living too deep for conventional sampling, was observed occupying an opportunistic niche in areas of high deposition and in very close association with Capitellid worm complexes. The surface-dwelling brittlestar, Ophiocomina nigra, was abundant at distances of 250-1000 m from Farm-B, suggesting a positive response to enrichment, but was displaced where sedimentation exceed 5 g m2 d-1. A corresponding gradient was evident within the sediment microbial communities, supporting established theories about ecosystem engineering and multi-species synergies for organic waste assimilation. Many of the bacteria present in the near-farm sediments were linked to the farmed fish and fish health issues suggesting one or two-way inoculation pressures. These functionally different benthic organisms are intrinsically linked and the resulting synergy has the potential to assimilate significant quantities of anthropogenically produced organic waste contributing to environmental sustainability.

Assessing the effects of salmon farming seabed enrichment using bacterial community diversity and high-throughput sequencing.

Aquaculture is an extremely valuable and rapidly expanding sector of the seafood industry. The sediment below active aquaculture farms receives inputs of organic matter from uneaten food and faecal material and this has led to concerns related to environmental sustainability. The impacts of organic enrichment on macrobenthic infauna are well characterized; however, much less is known about effect on bacterial communities. In this study, sediment, macrobenthic infauna samples and environmental data were collected along an enrichment gradient radiating out from a Chinook salmon (Oncorhynchus tshawytscha) farm (Marlborough Sounds; New Zealand). DNA and RNA were extracted and 16S rRNA metabarcodes from bacterial communities characterized using high-throughput sequencing. Desulfobacterales dominated at the cage (DNA and RNA), and at sites 50 m (DNA and RNA) and 150 m (RNA) from the farm. In contrast, unclassified bacteria from the class Gammaproteobacteria were the most abundant taxa at control sites (625 and 4000 m). Pronounced differences among DNA and RNA samples occurred at the cage site where Desulfobacterales abundance was markedly higher in RNA samples. There were strong correlations between shifts in bacterial communities and total organic matter and redox. This suggests that bacterial composition is strongly influenced by organic enrichment, a trait that may make them useful for assessing impacts associated with aquaculture farms.

Spatial and temporal dynamics in macrobenthos during recovery from salmon farm induced organic enrichment: when is recovery complete?

This study documents eight years of benthic recovery at a highly impacted salmon farm. Substantial recovery occurred in the first 2 years, and was assessed to be complete after ∼5 years. However, minor differences were still evident, along with some on-going benthic instability, attributable to medium-scale spatial movements and successional patterns of macrobenthos. Quantifying the endpoint of 'recovery' proved challenging due to: lack of a widely accepted definition, inherent variability in recovering sediments, differing trajectories of impact and reference sites, and statistical challenges. More complex biotic indices and metrics incorporating multiple variables were the most robust indicators. Statistical tests for 'parallelism' in the trajectories of Cage and Reference sites proved useful, but results were contingent upon how the method was applied, and should therefore be used in conjunction with data-visualisation methods. The study highlights the importance of a predetermined recovery endpoint, and using multiple indicators and a weight-of-evidence assessment approach.

Novel observations of benthic enrichment in contrasting flow regimes with implications for marine farm monitoring and management.

We examine macrofaunal and physico-chemical responses to organic enrichment beneath salmon farms in contrasting flow environments, and reveal pronounced flow-related differences in the magnitude and spatial extent of effects. Total macrofaunal abundances at high flow sites were nearly an order of magnitude greater than at comparable low flow sites, representing a significant benthic biomass. These very high abundances occurred in conjunction with moderate-to-high species richness, and were evident in the absence of appreciable organic matter accumulation. Biological responses to increasing sulfide were variable; however a significant biological threshold was evident at 1500 µM. Macrofaunal responses at high flow sites differed substantially from the Pearson-Rosenberg model. The atypical ecological conditions were attributed to (i) limited accumulation of fine sediments, (ii) maintenance of aerobic conditions in near-surface sediments, and (iii) an abundant food supply. Thus, enhanced resilience to organic waste at well-flushed sites appears related to both biological and physical processes.