Comparative metagenomic analysis of microbial community compositions and functions in cage aquaculture and its nearby non-aquaculture environments

In the context of burgeoning global aquaculture, its environmental repercussions, particularly in marine ecosystems, have gained significant attentions. Cage aquaculture, a prominent method, has been observed to significantly influence marine environments by discharging substantial amounts of organic materials and pollutants. It is also one of the important reasons for water eutrophication. This study investigated the impacts of cage aquaculture on microbial diversity and functional potential using metagenomics. Specifically, a comparison was made of the physicochemical indicators and microbial diversity between three grouper aquaculture cage nets in Lingshui Xincun Port and three nearby non-aquaculture area surface waters. We found that compared to non-aquaculture areas, the eutrophication indicators in aquaculture environments significantly increased, and the abundances of Vibrio and Pseudoalteromonas in aquaculture environments significantly rose. Additionally, microbial functional genes related to carbon, nitrogen, and sulfur metabolisms were also found to be significantly affected by aquaculture activities. The correlation analysis between microbial populations and environmental factors revealed that the abundances of most microbial taxa showed positive correlations with dissolved inorganic nitrogen, soluble reactive phosphorus, NH4+, and negative correlations with dissolved oxygen. Overall, this study elucidated the significant impacts of aquaculture-induced eutrophication on the diversity and functions of planktonic bacterial communities.

Environmental DNA methods for biosecurity and invasion biology in terrestrial ecosystems: Progress, pitfalls, and prospects

Analysis of environmental DNA (eDNA) enables indirect detection of species without the need to directly observe and sample them. For biosecurity and invasion biology, eDNA-based methods are useful to address biological invasions at all phases, from detecting arrivals to confirming eradication of past invasions. We conducted a systematic review of the literature and found that in biosecurity and invasion biology, eDNA has primarily been used to detect new incursions and monitor spread in marine and freshwater ecosystems, with much slower uptake in terrestrial ecosystems, reflecting a broader trend common to the usage of eDNA tools. In terrestrial ecosystems, eDNA research has mostly focussed on the use of eDNA metabarcoding to characterise biodiversity, rather than targeting biosecurity threats or non-native populations. We discuss how eDNA-based methods are being applied to terrestrial ecosystems for biosecurity and managing non-native populations at each phase of the invasion continuum: transport, introduction, establishment, and spread; across different management options: containment, control, and eradication; and for detecting the impact of non-native organisms. Finally, we address some of the current technical issues and caveats of eDNA-based methods, particularly for terrestrial ecosystems, and how these might be solved. As eDNA-based methods improve, they will play an increasingly important role in the early detection and adaptive management of biological invasions, and the implementation of effective biosecurity controls.

Sources and contamination routes of seafood with human pathogenic Vibrio spp.: A Farm-to-Fork approach

Vibrio spp., known human foodborne pathogens, thrive in freshwater, estuaries, and marine settings, causing vibriosis upon ingestion. The rising global vibriosis cases due to climate change necessitate a deeper understanding of Vibrio epidemiology and human transmission. This review delves into Vibrio contamination in seafood, scrutinizing its sources and pathways. We comprehensively assess the contamination of human-pathogenic Vibrio in the seafood chain, covering raw materials to processed products. A "Farm-to-Fork" approach, aligned with the One Health concept, is essential for grasping the complex nature of Vibrio contamination. Vibrio's widespread presence in natural and farmed aquatic environments establishes them as potential entry points into the seafood chain. Environmental factors, including climate, human activities, and wildlife, influence contamination sources and routes, underscoring the need to understand the origin and transmission of pathogens in raw seafood. Once within the seafood chain, the formation of protective biofilms on various surfaces in production and processing poses significant food safety risks, necessitating proper cleaning and disinfection to prevent microbial residue. In addition, inadequate seafood handling, from inappropriate processing procedures to cross-contamination via pests or seafood handlers, significantly contributes to Vibrio food contamination, thus warranting attention to reduce risks. Information presented here support the imperative for proactive measures, robust research, and interdisciplinary collaboration in order to effectively mitigate the risks posed by human pathogenic Vibrio contamination, safeguarding public health and global food security. This review serves as a crucial resource for researchers, industrials, and policymakers, equipping them with the knowledge to develop biosecurity measures associated with Vibrio-contaminated seafood.

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.

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.

A First Insight into the Microbial and Viral Communities of Comau Fjord—A Unique Human-Impacted Ecosystem in Patagonia (42∘ S)

While progress has been made in surveying the oceans to understand microbial and viral communities, the coastal ocean and, specifically, estuarine waters, where the effects of anthropogenic activity are greatest, remain partially understudied. The coastal waters of Northern Patagonia are of interest since this region experiences high-density salmon farming as well as other disturbances such as maritime transport of humans and cargo. Here, we hypothesized that viral and microbial communities from the Comau Fjord would be distinct from those collected in global surveys yet would have the distinctive features of microbes from coastal and temperate regions. We further hypothesized that microbial communities will be functionally enriched in antibiotic resistance genes (ARGs) in general and in those related to salmon farming in particular. Here, the analysis of metagenomes and viromes obtained for three surface water sites showed that the structure of the microbial communities was distinct in comparison to global surveys such as the Tara Ocean, though their composition converges with that of cosmopolitan marine microbes belonging to Proteobacteria, Bacteroidetes, and Actinobacteria. Similarly, viral communities were also divergent in structure and composition but matched known viral members from North America and the southern oceans. Microbial communities were functionally enriched in ARGs dominated by beta-lactams and tetracyclines, bacitracin, and the group macrolide–lincosamide–streptogramin (MLS) but were not different from other communities from the South Atlantic, South Pacific, and Southern Oceans. Similarly, viral communities were characterized by exhibiting protein clusters similar to those described globally (Tara Oceans Virome); however, Comau Fjord viromes displayed up to 50% uniqueness in their protein content. Altogether, our results indicate that microbial and viral communities from the Comau Fjord are a reservoir of untapped diversity and that, given the increasing anthropogenic impacts in the region, they warrant further study, specifically regarding resilience and resistance against antimicrobials and hydrocarbons.

Environmental nucleic acids: A field-based comparison for monitoring freshwater habitats using eDNA and eRNA

Nucleic acids released by organisms and isolated from environmental substrates are increasingly being used for molecular biomonitoring. While environmental DNA (eDNA) has received much attention, the potential of environmental RNA as a biomonitoring tool remains under-explored. Several recent studies using paired DNA and RNA metabarcoding of bulk samples suggest that RNA might better reflect "metabolically active" parts of the community. However, such studies mainly capture organismal eDNA and eRNA. For larger eukaryotes, isolation of extra-organismal RNA will be important, but viability needs to be examined in a field-based setting. In this study we evaluate (a) whether extra-organismal eRNA release from macroeukaryotes can be detected given its supposedly rapid degradation, and (b) if the same field collection methods for eDNA can be applied to eRNA. We collected eDNA and eRNA from water in lakes where fish community composition is well documented, enabling a comparison between the two nucleic acids in two different seasons with monitoring using conventional methods. We found that eRNA is released from macroeukaryotes and can be filtered from water and metabarcoded in a similar manner as eDNA to reliably provide species composition information. eRNA had a small but significantly greater true positive rate than eDNA, indicating that it correctly detects more species known to exist in the lakes. Given relatively small differences between the two molecules in describing fish community composition, we conclude that if eRNA provides significant advantages in terms of lability, it is a strong candidate to add to the suite of molecular monitoring tools.

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.

Environmental DNA metabarcoding for benthic monitoring: A review of sediment sampling and DNA extraction methods

Environmental DNA (eDNA) metabarcoding (parallel sequencing of DNA/RNA for identification of whole communities within a targeted group) is revolutionizing the field of aquatic biomonitoring. To date, most metabarcoding studies aiming to assess the ecological status of aquatic ecosystems have focused on water eDNA and macroinvertebrate bulk samples. However, the eDNA metabarcoding has also been applied to soft sediment samples, mainly for assessing microbial or meiofaunal biota. Compared to classical methodologies based on manual sorting and morphological identification of benthic taxa, eDNA metabarcoding offers potentially important advantages for assessing the environmental quality of sediments. The methods and protocols utilized for sediment eDNA metabarcoding can vary considerably among studies, and standardization efforts are needed to improve their robustness, comparability and use within regulatory frameworks. Here, we review the available information on eDNA metabarcoding applied to sediment samples, with a focus on sampling, preservation, and DNA extraction steps. We discuss challenges specific to sediment eDNA analysis, including the variety of different sources and states of eDNA and its persistence in the sediment. This paper aims to identify good-practice strategies and facilitate method harmonization for routine use of sediment eDNA in future benthic monitoring.

Aquaculture drives distinct patterns of planktonic and sedimentary bacterial communities: insights into co-occurrence pattern and assembly processes

Aquaculture would change the environmental condition in the lake ecosystem, affecting the structure and function of the aquatic ecosystem. However, little is known about the underlying mechanisms controlling the distribution patterns of bacterial community respond to aquaculture in water column and sediment. Here, we investigated the composition, co-occurrence patterns, and assembly processes of planktonic and sedimentary bacterial communities (PBC vs. SBC) from an aquaculture-influenced zone of the Eastern Lake Taihu, China. We found that aquaculture activity greatly influenced the diversity and composition of SBC by inducing excess nitrogen into the sediments. Meanwhile, network analysis revealed that aquaculture activity strengthened species interactions within the SBC network but weakened the species interactions within the PBC network. Aquaculture activity also increased the importance of deterministic processes governing the assembly of SBC by heightening the importance of environmental filtering, whereas it decreased the relative importance of deterministic processes within the assembly of PBC. In addition, ecological restoration with macrophytes increased the diversity of PBC and formed a more stable PBC network by increasing the number of network keystones. Overall, our results indicated that aquaculture drove distinct co-occurrence patterns and assembly mechanisms of PBC and SBC. This study has fundamental implications in the lake ecosystem for evaluating the microbially mediated ecological consequences of aquaculture.

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.

Deciphering the molecular signal from past and alive bacterial communities in aquatic sedimentary archives

Lake sediments accumulate information on biological communities thus acting as natural archives. Traditionally paleolimnology has focussed on fossilized remains of organisms, however, many organisms do not leave fossil evidence, meaning major ecosystem components are missing from environmental reconstructions. Many paleolimnology studies now incorporate molecular methods, including investigating microbial communities using environmental DNA (eDNA), but there is uncertainty about the contribution of living organisms to molecular inventories. In the present study, we obtained DNA and RNA inventories from sediment spanning 700 years to investigate the contribution of past and active communities to the molecular signal from sedimentary archives. Additionally, a droplet digital PCR (ddPCR) targeting the 16S ribosomal RNA (16S rRNA) gene of the photosynthetic cyanobacterial genera Microcystis was used to explore if RNA signals were from legacy RNA. We posit that the RNA signal is a mixture of legacy RNA, dormant cells, living bacteria and modern-day trace level contaminants that were introduced during sampling and preferentially amplified. The presence of legacy RNA was confirmed by the detection of Microcystis in sediments aged to ~200 years ago. Recent comparisons between 16S rRNA gene metabarcoding and traditional paleo proxies showed that past changes in bacterial communities can be reconstructed from sedimentary archives. The recovery of RNA in the present study has provided new insights into the origin of these signals. However, caution is required during analysis and interpretation of 16S rRNA gene metabarcoding data especially in recent sediments were there are potentially active bacteria.

Interaction of environmental eukaryotic microorganisms and fungi in the pond-cultured carps: new insights into the potential pathogenic fungi in the freshwater aquaculture

The quality and safety of the aquatic products have gradually become the focus of global attention. In this study, the environmental eukaryotic and fungi communities in pond-cultured grass carp (Ctenopharyngodon idellus) and the koi carp (Cyprinus carpio L.) were investigated. For comparative analysis, the alpha diversity shows that the environmental microbial abundance in the koi carp groups were higher than that in the grass carp groups, while beta diversity reveals that the differences of the microbial community composition and structures in the grass carp groups were significantly higher than those in the koi carp groups. Meanwhile, the environmental microbial diversity of grass carp groups was higher than that of koi carp groups at phylum level, but showed no significant difference at genus level. Additionally, the dominant total phyla were Opisthokonta, Stramenopiles plusAlveolates plusRhizaria, Archaeplastida, Cryptophyceae, and Centrohelida for the 18S rRNA gene and Ciliophora, Chlorophyta, and Ascomycota for the ITS2 rRNA gene in both of the two carp groups. Additionally, annotation analysis showed that the biomarkers in the grass carp groups are significantly higher than those of the koi carp groups. Furthermore, the functional prediction of Funguild showed significant difference in outputs, while similarity in trophic modes and guild types between the two carp groups. Meanwhile, the total relative abundances of animal pathogen, fungal parasite, and plant pathogen were extremely similar between the two carp groups. Surprisingly, one pathogenic fungus of genus Fusarium was identified in both the environments of two carp groups based on filtered operational taxonomic unit tables. Overall, this is the first robust report to understand the characteristics of environmental eukaryotic microorganisms and fungi in the edible and ornamental carps. Our results also provide the basic data for the prevention of fungal diseases and the healthy culture of the carps.

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.

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.

A step towards the validation of bacteria biotic indices using DNA metabarcoding for benthic monitoring

Environmental genomics is a promising field for monitoring biodiversity in a timely fashion. Efforts have increasingly been dedicated to the use of bacteria DNA derived data to develop biotic indices for benthic monitoring. However, a substantial debate exists about whether bacteria-derived data using DNA metabarcoding should follow, for example, a taxonomy-based or a taxonomy-free approach to marine bioassessments. Here, we showcase the value of DNA-based monitoring using the impact of fish farming as an example of anthropogenic disturbances in coastal areas and compare the performance of taxonomy-based and taxonomy-free approaches in detecting environmental alterations. We analysed samples collected near to the farm cages and along distance gradients from two aquaculture installations, and at control sites, to evaluate the effect of this activity on bacterial assemblages. Using the putative response of bacterial taxa to stress we calculated the taxonomy-based biotic index microgAMBI. The distribution of individual amplicon sequence variants (ASVs), as a function of a gradient in sediment acid volatile sulphides, was then used to derive a taxonomy-free bacterial biotic index specific for this data set using a de novo approach based on quantile regression splines. Our results show that microgAMBI revealed a organically enriched environment along the gradient. However, the de novo biotic index outperformed microgAMBI by providing a higher discriminatory power in detecting changes in abiotic factors directly related to fish production, whilst allowing the identification of new ASVs bioindicators. The de novo strategy applied here represents a robust method to define new bioindicators in regions or habitats where no previous information about the response of bacteria to environmental stressors exists.

Towards the Optimization of eDNA/eRNA Sampling Technologies for Marine Biosecurity Surveillance

The field of eDNA is growing exponentially in response to the need for detecting rare and invasive species for management and conservation decisions. Developing technologies and standard protocols within the biosecurity sector must address myriad challenges associated with marine environments, including salinity, temperature, advective and deposition processes, hydrochemistry and pH, and contaminating agents. These approaches must also provide a robust framework that meets the need for biosecurity management decisions regarding threats to human health, environmental resources, and economic interests, especially in areas with limited clean-laboratory resources and experienced personnel. This contribution aims to facilitate dialogue and innovation within this sector by reviewing current approaches for sample collection, post-sampling capture and concentration of eDNA, preservation, and extraction, all through a biosecurity monitoring lens.

Metabarcoding on both environmental DNA and RNA highlights differences between fungal communities sampled in different habitats

In recent years, metabarcoding has become a key tool to describe microbial communities from natural and artificial environments. Thanks to its high throughput nature, metabarcoding efficiently explores microbial biodiversity under different conditions. It can be performed on environmental (e)DNA to describe so-called total microbial community, or from environmental (e)RNA to describe active microbial community. As opposed to total microbial communities, active ones exclude dead or dormant organisms. For what concerns Fungi, which are mostly filamentous microorganisms, the relationship between DNA-based (total) and RNA-based (active) communities is unclear. In the present study, we evaluated the consequences of performing metabarcoding on both soil and wood-extracted eDNA and eRNA to delineate molecular operational taxonomic units (MOTUs) and differentiate fungal communities according to the environment they originate from. DNA and RNA-based communities differed not only in their taxonomic composition, but also in the relative abundances of several functional guilds. From a taxonomic perspective, we showed that several higher taxa are globally more represented in either “active” or “total” microbial communities. We also observed that delineation of MOTUs based on their co-occurrence among DNA and RNA sequences highlighted differences between the studied habitats that were overlooked when all MOTUs were considered, including those identified exclusively by eDNA sequences. We conclude that metabarcoding on eRNA provides original functional information on the specific roles of several taxonomic or functional groups that would not have been revealed using eDNA alone.

Environmental DNA metabarcoding reveals estuarine benthic community response to nutrient enrichment - Evidence from an in-situ experiment

Nutrient loading is a major threat to estuaries and coastal environments worldwide, therefore, it is critical that we have good monitoring tools to detect early signs of degradation in these ecologically important and vulnerable ecosystems. Traditionally, bottom-dwelling macroinvertebrates have been used for ecological health assessment but recent advances in environmental genomics mean we can now characterize less visible forms of biodiversity, offering a more holistic view of the ecosystem and potentially providing early warning signals of disturbance. We carried out a manipulative nutrient enrichment experiment (0, 150 and 600 g N fertilizer m^-2) in two estuaries in New Zealand to assess the effects of nutrient loading on benthic communities. After seven months of enrichment, environmental DNA (eDNA) metabarcoding was used to examine the response of eukaryotic (18S rRNA), diatom only (rbcL) and bacterial (16S rRNA) communities. Multivariate analyses demonstrated changes in eukaryotic, diatom and bacterial communities in response to nutrient enrichment at both sites, despite differing environmental conditions. These patterns aligned with changes in macrofaunal communities identified using traditional morphological techniques, confirming concordance between disturbance indicators detected by eDNA and current monitoring approaches. Clear shifts in eukaryotic and bacterial indicator taxa were seen in response to nutrient loading while changes in diatom only communities were more subtle. Community changes were discernible between 0 and 150 g N m^-2 treatments, suggesting that estuary health assessment tools could be developed to detect early signs of degradation. Increasing variation in community structure associated with nutrient loading could also be used as an indicator of stress or approaching tipping points. This work represents a first step towards the development of molecular-based estuary monitoring tools, which could provide a more holistic and standardized approach to ecosystem health assessment with faster turn-around times and lower costs.

Spatiotemporal Dynamics of Vibrio Communities and Abundance in Dongshan Bay, South of China

The Vibrio genus inhabit estuarine and marine ecosystem throughout the world and can cause severe infections in humans and animals. Previous studies have demonstrated the dynamics of Vibrio at both community and population levels and assessed the close relationship between environmental factors and Vibrio diversity and abundance, such as temperature, salinity, and dissolved oxygen. It is also generally believed that aquaculture is the fastest-growing food sector, which is also applying great environmental impacts on microbial communities in aquatic ecosystems. However, our understanding of the spatiotemporal quantification of Vibrio throughout the four seasons in the aquaculture zone and response to environmental factors remains poor. To explore the spatiotemporal distribution and abundance of the Vibrio community with their related environmental factors and detect the relationships among them, we collected 10 seawater sites spanning four seasons across the whole year in Dongshan Bay for investigating the Vibrio community dynamics. Marked differences in diversity and abundance of the Vibrio community were observed between seasons, which were mainly driven by temperature, dissolved oxygen, nitrate, and nitrite. qPCR analysis showed that Vibrio abundance was most abundant in the summer (5.37 × 10⁶ copies/L), compared with the autumn (4.58 × 10⁶ copies/L), spring (1.18 × 10⁶ copies/L), and winter (1.55 × 10⁴ copies/L). A total of 22 Vibrio operational taxonomic units (OTUs) and 28 species were identified by universal bacteria 16S rRNA gene and cultivation methods, with Vibrio fortis the dominant in these aquaculture areas. To summarize, our present study is one of the few studies to research the occurrence of Vibrio in marine aquaculture of South China, and the results indicate that Vibrio are widely distributed in aquaculture environment and that a further risk assessment is needed to be conducted.

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.

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.

A microbial mandala for environmental monitoring: Predicting multiple impacts on estuarine prokaryote communities of the Bay of Biscay

Routine monitoring of benthic biodiversity is critical for managing and understanding the anthropogenic impacts on marine, transitional and freshwater ecosystems. However, traditional reliance on morphological identification generally makes it cost-prohibitive to increase the scale of monitoring programmes. Metabarcoding of environmental DNA has clear potential to overcome many of the problems associated with traditional monitoring, with prokaryotes and other microorganisms showing particular promise as bioindicators. However, due to the limited knowledge regarding the ecological roles and responses of environmental microorganisms to different types of pressure, the use of de novo approaches is necessary. Here, we use two such approaches for the prediction of multiple impacts present in estuaries and coastal areas of the Bay of Biscay based on microbial communities. The first (Random Forests) is a machine learning method while the second (Threshold Indicator Taxa Analysis and quantile regression splines) is based on de novo identification of bioindicators. Our results show that both methods overlap considerably in the indicator taxa identified, but less for sequence variants. Both methods also perform well in spite of the complexity of the studied ecosystem, providing predictive models with strong correlation to reference values and fair to good agreement with ecological status groups. The ability to predict several specific types of pressure is especially appealing. The cross-validated models and biotic indices developed can be directly applied to predict the environmental status of estuaries in the same geographical region, although more work is needed to evaluate and improve them for use in new regions or habitats.

Supervised machine learning is superior to indicator value inference in monitoring the environmental impacts of salmon aquaculture using eDNA metabarcodes

Increasing anthropogenic impact and global change effects on natural ecosystems has prompted the development of less expensive and more efficient bioassessments methodologies. One promising approach is the integration of DNA metabarcoding in environmental monitoring. A critical step in this process is the inference of ecological quality (EQ) status from identified molecular bioindicator signatures that mirror environmental classification based on standard macroinvertebrate surveys. The most promising approaches to infer EQ from biotic indices (BI) are supervised machine learning (SML) and the calculation of indicator values (IndVal). In this study we compared the performance of both approaches using DNA metabarcodes of bacteria and ciliates as bioindicators obtained from 152 samples collected from seven Norwegian salmon farms. Results from standard macroinvertebrate-monitoring of the same samples were used as reference to compare the accuracy of both approaches. First, SML outperformed the IndVal approach to infer EQ from eDNA metabarcodes. The Random Forest (RF) algorithm appeared to be less sensitive to noisy data (a typical feature of massive environmental sequence data sets) and uneven data coverage across EQ classes (a typical feature of environmental compliance monitoring scheme) compared to a widely used method to infer IndVals for the calculation of a BI. Second, bacteria allowed for a more accurate EQ assessment than ciliate eDNA metabarcodes. For the implementation of DNA metabarcoding into routine monitoring programmes to assess EQ around salmon aquaculture cages, we therefore recommend bacterial DNA metabarcodes in combination with SML to classify EQ categories based on molecular signatures.

Marine mammals are natural hosts of Oceanivirga salmonicida, a bacterial pathogen of Atlantic salmon

During 1992 and 1993, a bacterial disease occurred in a seawater Atlantic salmon Salmo salar farm, causing serious mortalities. The causative agent was subsequently named as Oceanivirga salmonicida, a member of the Leptotrichiaceae. Searches of 16S rRNA gene sequence databases have shown sequence similarities between O. salmonicida and uncultured bacterial clones from the digestive tracts of marine mammals. In the current study, oral samples were taken from stranded dolphins (common dolphin Delphinus delphis, striped dolphin Stenella coeruleoalba) and healthy harbour seals Phoca vitulina. A bacterium with growth characteristics consistent with O. salmonicida was isolated from a common dolphin. The isolate was confirmed as O. salmonicida, by comparisons to the type strain, using 16S rRNA gene, gyrB, groEL, and recA sequence analyses, average nucleotide identity analysis, and MALDI-TOF mass spectrometry. Metagenomic analysis indicated that the genus Oceanivirga represented a significant component of the oral bacterial microbiomes of the dolphins and seals. However, sequences consistent with O. salmonicida were only found in the dolphin samples. Analyses of marine mammal microbiome studies in the NCBI databases showed sequences consistent with O. salmonicida from the common dolphin, striped dolphin, bottlenose dolphin Tursiops truncatus, humpback whale Megaptera novaeangliae, and harbour seal. Sequences from marine environmental studies in the NCBI databases showed no sequences consistent with O. salmonicida. The findings suggest that several species of marine mammals are natural hosts of O. salmonicida.

Release and degradation of environmental DNA and RNA in a marine system

Over the last decade, there has been growing interest in the analysis of environmental DNA (eDNA) to infer the presence of organisms in aquatic environments. The efficacy of eDNA/eRNA based tools are highly depend on the turnover rate of the molecule (their release and degradation). Environmental DNA has been shown to persist for days, weeks or years in environmental samples. Environmental RNA (eRNA) is thought to degrade faster than eDNA, however to our knowledge, no experimental studies have explored this. Here we present an aquarium study to investigate eDNA and eRNA shedding rates and degradation for two sessile marine invertebrates. The copy numbers for eDNA and eRNA were assessed using droplet digital PCR targeting the mitochondrial Cytochrome c Oxidase subunit 1 (COI) gene. Environmental RNA persisted after organism removal for much longer than expected with detections for up to 13 h. In contrast, eDNA was detected is samples collected up to 94 h after organism removal. There was no evidence that the decay rates constants for eDNA and eRNA were different (p = 0.6, Kruskal-Wallis tests). Both eDNA and eRNA was detected in biofilms collected at the end of the experiment (day 21). This suggests binding with organic or inorganic compounds or stabilization of these molecules in the biofilm matrix. The finding of the prolonged persistence of eRNA may provide new opportunities for improved biodiversity surveys through reducing false positives caused by legacy DNA and could also facilitate new research on environmental transcriptomics.

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 m² 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.

Cumulative impact of long-term intensive mariculture on total and active bacterial communities in the core sediments of the Ailian Bay, North China

The exponential growth of off-shore mariculture worldwide over the last 20 years has had significant impact on coastal sediment biogeochemistry. However, there are no long-term records of the cumulative impacts of mariculture on the benthic bacterial community. Here, total (DNA) and active (RNA) bacterial community compositions were characterized using MiSeq sequencing of 16S rRNA gene in four core sediments of the Ailian Bay, one of the typical intensive mariculture areas in China with more than fifty-year history of kelp and scallop cultivation. The γ-Proteobacteria, δ-Proteobacteria, Acidobacteria and Acitinobacteria were more abundant in the total bacterial communities, while β-Proteobacteria, Anaerolineae, Clostridia, Spirochaetes and Cyanobacteria were enriched in the active bacterial communities. Significant differences were observed between total and active benthic bacterial communities. The influences of different mariculture modes on the total bacterial communities were more significant than those on the active bacterial communities. Only limited groups of the total bacterial communities were significant influenced by the cumulative effects of the long-term mariculture. The bacterial genera with the function in the sulfide cycling and organic consumption were enriched in the total bacterial population of the integrated multi-trophic aquaculture (IMTA) areas. The variations of both total and active bacterial communities were significantly influenced by grain sizes, total organic carbon and nutrients. Both total and active bacterial communities exhibited a slightly stronger response to environmental factors than to spatial (distance) factors. The effects of mutualism might dominate the total and active bacterial networks in the Ailian Bay. The present study demonstrated that the cumulative influences of the long-term and intensive IMTA mariculture on total benthic bacterial communities in the sub-surface sediments of the Ailian Bay were stronger than those on the active benthic bacterial communities, which provided some insights into the potential ecological roles of specific taxa in the sediments of the IMTA ecosystems.

Seasonality overwhelms aquacultural activity in determining the composition and assembly of the bacterial community in Lake Taihu, China

Aquaculture accounts for an extremely valuable and rapidly expanding sector of global food production, yet its environmental impacts on aquatic ecosystems have attracted much concern. In this study, we collected water samples from eastern Lake Taihu, China. We targeted sites varying in their intensity of aquacultural activities and sampled them over multiple seasons. For each sample, we measured physicochemical variables, and we sequenced the 16S rRNA gene of the respective bacterial communities using an Illumina second-generation sequencing platform. Marked differences in diversity and bacterial community composition were observed between seasons, whereas we observed relatively weak differences between sites. Remarkable differences in the abundance of the bacterial community were observed at the phylum and genus levels across the different seasons. Stochastic processes dominated the assembly of bacterial communities in the aquaculture-influenced systems, and the assembly processes of bacterial community differed between seasons. Our observations highlight the effect of seasonality on bacterial communities and provide a more complete knowledge base for the proper assessment of the effects of aquacultural activities on freshwater ecosystems.

Influence of suspended mariculture on vertical distribution profiles of bacteria in sediment from Daya Bay, Southern China

Mariculture is known to contribute to oxygen depletion, pH decline and accumulation of nutrients and organic matter in sediments. However, studies on the bacterial vertical distribution of mariculture area are very limited. The bacterial abundance in the non-culture site (3.8 ± 0.8 × 10⁹ copies g^-1) was significantly higher than that in the three mariculture sites (1.2 ± 0.2 × 10⁹ copies g^-1), and bacterial diversity in the non-culture site was significantly higher than that in fish cage-TF (p < 0.05). The vertical distribution profiles of bacteria in non-culture and oyster culture sites were similar but very different from that of fish cage-TF. In addition, significant downward trends in bacterial abundance and diversity were observed as sediment depth increased (p < 0.05), and the most relevant environmental factors were moisture content, total nitrogen, total organic carbon and carbon/nitrogen. The dominant bacterial phyla in sediment were Proteobacteria, Chloroflexi and Bacteroidetes.

Multi-marker eDNA metabarcoding survey to assess the environmental impact of three offshore gas platforms in the North Adriatic Sea (Italy)

The environmental DNA (eDNA) metabarcoding represents a new promising tool for biomonitoring and environmental impact assessment. One of the main advantages of eDNA metabarcoding, compared to the traditional morphotaxonomy-based methods, is to provide a more holistic biodiversity information that includes inconspicuous morphologically non-identifiable taxa. Here, we use eDNA metabarcoding to survey marine biodiversity in the vicinity of the three offshore gas platforms in North Adriatic Sea (Italy). We isolated eDNA from 576 water and sediment samples collected at 32 sampling sites situated along four axes at increasing distances from the gas platforms. We obtained about 46 million eDNA sequences for 5 markers from nuclear 18S V1V2, 18S V4, 18S 37F and mitochondrial 16S and COI genes that cover a wide diversity of benthic and planktonic eukaryotes. Our results showed some impact of platform activities on benthic and pelagic communities at very close distance (<50 m), while communities for intermediate (125 m, 250 m, 500 m) and reference (1000 m, 2000 m) sites did not show any particular biodiversity changes that could be related to platforms activities. The most significant community change along the distance gradient was obtained with the 18S V1V2 marker targeting benthic eukaryotes, even though other markers showed similar trends, but to a lesser extent. These results were congruent with the AMBI index inferred from the eDNA sequences assigned to benthic macrofauna. We finally explored the relation between various physicochemical parameters, including hydrocarbons, on benthic community in the case of one of the platforms. Our results showed that these communities were not significantly impacted by most of hydrocarbons, but rather by macro-elements and sediment texture.

Temporal Bacterial Surveillance of Salmon Aquaculture Sites Indicates a Long Lasting Benthic Impact With Minimal Recovery

Coastal aquaculture has experienced substantial growth in the last few decades and associated impacts on natural environments are of increasing importance. Understanding both the effects of aquaculture on marine ecosystems and the processes of recovery during fallowing periods is crucial for the development of a more environmentally sustainable industry. Because bacteria are sensitive to environmental change, surveying fluctuations in bacterial communities is a promising tool for monitoring the status of benthic environments. Here, we used 16S rRNA gene high-throughput sequencing to characterize bacterial communities in flocculent matter samples collected over a period of 3 years and at various distances from cages (0–200 meters) at production and fallow (3–35 months) salmon aquaculture sites in southern Newfoundland to evaluate the environmental impact of aquaculture on predominantly hard-bottom substrates. Bacterial composition analysis revealed four clusters, three of which (defined as “recently disturbed,” “intermediate impact,” and “high impact”) differed markedly from a fourth “low impact” cluster that contained far-field samples collected >500 m from cages. Samples within the high impact group were most often collected directly under cages, whereas those in the intermediate impact group were mainly sampled from 20 to 40 m from cages. Large scale phylum shifts (increases of Bacteroidetes, Firmicutes, Spirochaetes, and decreases in Proteobacteria and Epsilonbacteraeota) and a decline in bacterial diversity were observed in the high impact cluster, indicating significant ecological change. Samples from sites of different fallow duration were found in the high impact cluster, indicating a lack of recovery, even after 35 months of fallowing. Finally, we identified 28 genera as bacterial biomarkers, specific to one or more clusters, including genera associated with organically enriched environments and previously reported in the context of aquaculture impacts. Tracking the relative abundance of biomarkers in relation to different lengths of fallowing in the three more impacted clusters showed that these markers remained significantly above low impact cluster levels at all times, further pointing toward incomplete recovery. Our results suggest that coastal aquaculture on hard-bottom substrates is prone to long lasting impacts on bacterial communities, especially below cages, and that effects can be accurately tracked using bacterial community profiles or specific biomarkers.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Dynamics of the bacterial community in a channel catfish nursery pond with a cage-pond integration system

The changes in the bacterial community composition in a channel catfish nursery pond with a cage-pond integration system were investigated by sequencing of the 16S rRNA gene through Illumina MiSeq sequencing platforms. A total of 1 362 877 sequences and 1440 operational taxonomic units were obtained. Further analysis showed that the dominant phyla in the cage and pond groups were similar, including Actinobacteria, Cyanobacteria, Proteobacteria, and Bacteroidetes, although a significant difference was detected between them by ANOSIM (P < 0.05). Temporal changes and site variation were significantly related to the variation of the bacterial community. A comprehensive analysis of the diversity and evenness of the bacterial 16S rRNA gene, redundancy analysis (RDA), and partial Mantel test showed that the bacterial community composition in a cage-pond integration system was shaped more by temporal variation than by site variation. RDA also indicated that water temperature, total dissolved solids, and Secchi depth had the largest impact on bacterial populations.

Environmental DNA metabarcoding of benthic bacterial communities indicates the benthic footprint of salmon aquaculture

We evaluated benthic bacterial communities as bioindicators in environmental impact assessments of salmon aquaculture, a rapidly growing sector of seafood industry. Sediment samples (n=72) were collected from below salmon cages towards distant reference sites. Bacterial community profiles inferred from DNA metabarcodes were compared to reference data from standard macrofauna biomonitoring surveys of the same samples. Deltaproteobacteria were predominant in immediate vicinity of the salmon cages. Along the transect, significant shifts in bacterial community structures were observed with Gammaproteobacteria dominating the less-impacted sites. Alpha- and beta-diversity measures of bacterial communities correlated significantly with macrofauna diversity metrics and with five ecological status indices. Benthic bacterial communities mirror the reaction of macrofauna bioindicators to environmental disturbances caused by salmon farming. The implementation of bacterial eDNA metabarcoding in future Strategic Framework Directives is an alternative cost-effective high-throughput biomonitoring solution, providing a basis for management strategies in a matter of days rather than months.

A cross-taxa study using environmental DNA/RNA metabarcoding to measure biological impacts of offshore oil and gas drilling and production operations

Standardized ecosystem-based monitoring surveys are critical for providing information on marine ecosystem health. Environmental DNA/RNA (eDNA/eRNA) metabarcoding may facilitate such surveys by quickly and effectively characterizing multi-trophic levels. In this study, we assessed the suitability of eDNA/eRNA metabarcoding to evaluate changes in benthic assemblages of bacteria, Foraminifera and other eukaryotes along transects at three offshore oil and gas (O&G) drilling and production sites, and compared these to morphologically characterized macro-faunal assemblages. Bacterial communities were the most responsive to O&G activities, followed by Foraminifera, and macro-fauna (the latter assessed by morphology). The molecular approach enabled detection of hydrocarbon degrading taxa such as the bacteria Alcanivorax and Microbulbifer at petroleum impacted stations. Most identified indicator taxa, notably among macro-fauna, were highly specific to site conditions. Based on our results we suggest that eDNA/eRNA metabarcoding can be used as a stand-alone method for biodiversity assessment or as a complement to morphology-based monitoring approaches.

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.

Bacterial communities 16S rDNA fingerprinting as a potential tracing tool for cultured seabass Dicentrarchus labrax

Traceability of seafood has become crucial with market globalization and consumer’s awareness. The present study used PCR-DGGE and 454 pyrosequencing to assess if bacterial communities fingerprint associated to seabass (Dicentrarchus labrax) skin mucus can be used to discriminate the geographic origin of fishes cultured in three semi-intensive fish farms. PCR-DGGE and pyrosequencing results were congruent and suggested that this molecular approach has the potential to trace fish farms with a spatial resolution <500 m. Pyrosequencing results provided a detailed insight into the bacterial community composition of seabass skin mucus and revealed the existence of a core of bacterial communities within family Pseudomonadaceae and Rhodobacteraceae. This approach also allowed to recognized key OTUs that are potentially relevant to discriminate the geographic origin of the fish being surveyed. Overall, the present study increased our knowledge on farmed seabass microbiome and demonstrated that specific and unique bacterial taxa can act as natural signatures that allow us to trace fish to its respective geographic origin. Our study provides valuable clues that should be more investigated in future studies as a way to fulfill current traceability needs in the global trade of seafood.

Dispersal and assimilation of an aquaculture waste subsidy in a low productivity coastal environment

To understand dispersal and assimilation of aquaculture waste subsidies in a naturally low-productivity environment, we applied a novel, rapid transmethylation technique to analyse sediment and biota fatty acid composition. This technique was initially validated at Atlantic salmon farms in Macquarie Harbour, Australia, where sediments were collected at farm and control locations. Subsequently, sediment, benthic polychaete and zooplankton were sampled at sites 0, 50, 250, 500 and 1000m distant from multiple cages. Results demonstrated an acute deposition zone up to 50m from cages and a diffuse zone extending 500m from cages. Changes in sediment concentration of linoleic acid, oleic acid and total fatty acids were effective tracers of farm deposition. Bacterial biomarkers indicated that aquaculture waste stimulates bacterial productivity in sediments, with elevated biomarker concentrations also detected in benthic polychaetes. Overall, fatty acid analysis was a sensitive technique to characterize the benthic footprint of aquaculture influence.

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.

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.

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.