Community structure changes in sediment bacterial communities along an organic enrichment gradient associated with fish farming

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.

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.

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.

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.

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.

Seasonal dynamics of the coastal bacterioplankton at intensive fish-farming areas of the Yellow Sea, China revealed by high-throughput sequencing

Marine aquaculture areas are facing stressed environmental challenges, especially the degradation of coastal ecosystems. Here a coordinated time-series study was used to investigate the coastal bacterioplankton biodiversity dynamics of the Yellow Sea, China. Bacterial 16S rRNA gene sequencing revealed a temporal pattern of decreasing of diversity in summer. Functional prediction indicated that metabolic pathways related to the adenosine triphosphate (ATP)-binding cassette transporters and other membrane transporters were significantly enriched in May, while the genetic information processing category was most abundant in March. The May microbiomes showed most significant positive correlation with phosphate concentration, while the August and November microbiomes correlated with temperature and chemical oxygen demand (COD) most, and the March microbiomes showed significant correlation with Cu²⁺ level, pH and salinity. The correlations between representative bacteria and environmental parameters revealed in this study may provide insights into the potential influences of human aquaculture activities, on the biodiversity of coastal bacterioplankton.

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.

Detecting the presence of fish farm-derived organic matter at the seafloor using stable isotope analysis of phospholipid fatty acids

The expansion of global aquaculture activities is important for the wellbeing of future generations in terms of employment and food security. Rearing animals in open-exchange cages permits the release of organic wastes, some of which ultimately reaches the underlying sediments. The development of rapid, quantitative and objective monitoring techniques is therefore central to the environmentally sustainable growth of the aquaculture industry. Here, we demonstrate that fish farm-derived organic wastes can be readily detected at the seafloor by quantifying sediment phospholipid fatty acids (PLFAs) and their carbon stable isotope signatures. Observations across five farms reveal that farm size and/or distance away from it influence the spatial distribution of the generated organic wastes and their effect on benthic bacterial biomass. Comparison to the isotopic signatures of fish feed-derived PLFAs indicates that 16:0 and 18:1(n-9) are potential biomarkers for fish farm-derived organic wastes. Our results suggest that stable isotope analysis of sediment PLFAs has potential for monitoring the environmental performance of aquaculture activities, particularly given the increasing prevalence of terrigenous organic matter in aquaculture feed stocks because it is isotopically district to marine organic matter.

Nutrient Limitation in Surface Waters of the Oligotrophic Eastern Mediterranean Sea: an Enrichment Microcosm Experiment

The growth rates of planktonic microbes in the pelagic zone of the Eastern Mediterranean Sea are nutrient limited, but the type of limitation is still uncertain. During this study, we investigated the occurrence of N and P limitation among different groups of the prokaryotic and eukaryotic (pico-, nano-, and micro-) plankton using a microcosm experiment during stratified water column conditions in the Cretan Sea (Eastern Mediterranean). Microcosms were enriched with N and P (either solely or simultaneously), and the PO4 turnover time, prokaryotic heterotrophic activity, primary production, and the abundance of the different microbial components were measured. Flow cytometric and molecular fingerprint analyses showed that different heterotrophic prokaryotic groups were limited by different nutrients; total heterotrophic prokaryotic growth was limited by P, but only when both N and P were added, changes in community structure and cell size were detected. Phytoplankton were N and P co-limited, with autotrophic pico-eukaryotes being the exception as they increased even when only P was added after a 2-day time lag. The populations of Synechococcus and Prochlorococcus were highly competitive with each other; Prochlorococcus abundance increased during the first 2 days of P addition but kept increasing only when both N and P were added, whereas Synechococcus exhibited higher pigment content and increased in abundance 3 days after simultaneous N and P additions. Dinoflagellates also showed opportunistic behavior at simultaneous N and P additions, in contrast to diatoms and coccolithophores, which diminished in all incubations. High DNA content viruses, selective grazing, and the exhaustion of N sources probably controlled the populations of diatoms and coccolithophores.

Biosurfactant production from marine hydrocarbon-degrading consortia and pure bacterial strains using crude oil as carbon source

Biosurfactants (BSs) are "green" amphiphilic molecules produced by microorganisms during biodegradation, increasing the bioavailability of organic pollutants. In this work, the BS production yield of marine hydrocarbon degraders isolated from Elefsina bay in Eastern Mediterranean Sea has been investigated. The drop collapse test was used as a preliminary screening test to confirm BS producing strains or mixed consortia. The community structure of the best consortia based on the drop collapse test was determined by 16S-rDNA pyrotag screening. Subsequently, the effect of incubation time, temperature, substrate and supplementation with inorganic nutrients, on BS production, was examined. Two types of BS - lipid mixtures were extracted from the culture broth; the low molecular weight BS Rhamnolipids and Sophorolipids. Crude extracts were purified by silica gel column chromatography and then identified by thin layer chromatography and Fourier transform infrared spectroscopy. Results indicate that BS production yield remains constant and low while it is independent of the total culture biomass, carbon source, and temperature. A constant BS concentration in a culture broth with continuous degradation of crude oil (CO) implies that the BS producing microbes generate no more than the required amount of BSs that enables biodegradation of the CO. Isolated pure strains were found to have higher specific production yields than the complex microbial marine community-consortia. The heavy oil fraction of CO has emerged as a promising substrate for BS production (by marine BS producers) with fewer impurities in the final product. Furthermore, a particular strain isolated from sediments, Paracoccus marcusii, may be an optimal choice for bioremediation purposes as its biomass remains trapped in the hydrocarbon phase, not suffering from potential dilution effects by sea currents.