Responses of foraminifera communities to aquaculture‐derived organic enrichment as revealed by environmental DNA metabarcoding

Sediment DNA metabarcoding and morphology provide complementary insight into macrofauna and meiobenthos response to environmental gradients in an Arctic glacial fjord

Arctic fjords ecosystems are highly dynamic, with organisms exposed to various natural stressors along with productivity clines driven by advection of water masses from shelves. The benthic response to these environmental clines has been extensively studied using traditional, morphology-based approaches mostly focusing on macroinvertebrates. In this study we analyse the effects of glacially mediated disturbance on the biodiversity of benthic macrofauna and meiobenthos (meiofauna and Foraminifera) in a Svalbard fjord by comparing morphology and eDNA metabarcoding. Three genetic markers targeting metazoans (COI), meiofauna (18S V1V2) and Foraminifera (18S 37f) were analyzed. Univariate measures of alpha diversity and multivariate compositional dissimilarities were calculated and tested for similarities in response to environmental gradients using correlation analysis. Our study showed different taxonomic composition of morphological and molecular datasets for both macrofauna and meiobenthos. Some taxonomic groups while abundant in metabarcoding data were almost absent in morphology-based inventory and vice versa. In general, species richness and diversity measures in macrofauna morphological data were higher than in metabarcoding, and similar for the meiofauna. Both methodological approaches showed different patterns of response to the glacially mediated disturbance for the macrofauna and the meiobenthos. Macrofauna showed an evident distinction in taxonomic composition and a dramatic cline in alpha diversity indices between the outer and inner parts of fjord, while the meiobenthos showed a gradual change and more subtle responses to environmental changes along the fjord axis. The two methods can be seen as complementing rather than replacing each other. Morphological approach provides more accurate inventory of larger size species and more reliable quantitative data, while metabarcoding allows identification of inconspicuous taxa that are overlooked in morphology-based studies. As different taxa may show different sensitivities to environmental changes, both methods shall be used to monitor marine biodiversity in Arctic ecosystems and its response to dramatically changing environmental conditions.

Assessment of anthropogenic pollution in Guanabara Bay (SE Brazil) through biogeochemical data and stable isotope mixing models

This work intends to identify pollution sources along the margins of Guanabara Bay (GB; SE Brazil) through a multiproxy approach and Bayesian stable isotopic mixture model (BSIMM). For this purpose, 33 surface sediment samples were collected and analyzed for granulometry, geochemistry (heavy metals, total organic carbon-TOC, stable isotopes of carbon and nitrogen-δ13C and δ15N, Rock-Eval pyrolysis parameters-REPP), and physicochemical parameters. Metal concentrations (E) dissolved in water (EW), adsorbed by organic matter (EOM) and by Mn hydroxides (EMn), and total extracted concentrations (ET) were analyzed. Sampling was conducted in 2018 after an oil spill from Reduc Oil Refinery. Potential Ecological risk index (PERI), based on metals, classified 85% of the analyzed stations as having moderate to considerable ecological risk. The metals with the potential to cause the highest ecological risk were CdW, CdOM, PbOM, and HgOM. The combination of BSIMM and REPP data was an effective proxy for oil spill detection by indicating the presence of polycyclic aromatic hydrocarbons (PAHs). Relatively high TOC contents suggested that the analyzed stations are eutrophicated environments. BSIMM discriminated three groups of stations with different sources of organic matter (OM), endorsing the result previously shown by the cluster analysis: (A) Niterói region, Botafogo marina, Glória marina, Fiscal and Fundão islands with diffuse sources of OM, including marine phytoplankton and material of continental origin from highly polluted rivers and domestic sewage; (B) region near Fundão and Governador islands and Mangue Channel outlet with OM (≃70%) supplied by highly polluted streams and a small contribution of PAHs; (C) Duque de Caxias and Botafogo-Urca inlet with significant contributions of PAHs, materials from C-3 plants and rivers polluted by urban sewage. Results of linear regressions in conjunction with BSIMM indicate that HgMn and PbOM mainly affect Group A's stations. Although the eastern margin of GB (Niterói; Group A) showed greater oceanic interaction than the other groups, it presented substantial concentrations of metals, potentially harmful (i.e., Hg and Pb) to marine biota and human health.

Benthic foraminiferal community structure and its response to environmental factors revealed using high-throughput sequencing in the Zhoushan Fishing Ground, East China Sea

Benthic foraminifera are excellent tools for monitoring marine environments and reconstructing paleoenvironments. This study investigated the structure and diversity of benthic foraminiferal communities in 20 superficial sediment samples obtained from the Zhoushan Fishing Ground (ZFG) using high-throughput sequencing based on small subunit ribosomal DNA and RNA amplification. The results revealed Rotaliida as the most dominant group, with spatial heterogeneity in foraminiferal distribution. Total benthic foraminiferal communities exhibited higher species richness and diversity compared to active communities. While heavy metal pollution in the ZFG was moderate, areas with elevated concentrations of heavy metals exhibited low diversity and richness in foraminiferal communities. Total foraminiferal community structure was primarily influenced by factors such as water depth and Hg, Pb, Cd, and Zn levels. Notably, Hg levels emerged as a critical factor impacting the structure and diversity of the active foraminiferal community. The dominant species, Operculina, exhibited tolerance toward heavy metal pollution.

Benthic foraminifera as bioindicators of marine pollution: A bibliometric approach to unravel trends, patterns and perspectives

Benthic foraminifera, single-celled marine organisms, are known for their wide distribution, high abundance and species diversity, test (i.e., shell) preservation in the sedimentary (e.g., historical) record, and sensitivity to environmental changes. Because of these characteristics, they have been widely used as bioindicators in environmental monitoring and, more recently, as Biological Quality Elements (BQEs) in the Ecological Quality Status (EcoQS) evaluation. The global scientific literature on benthic foraminifera as bioindicators was gathered from the Scopus database (overall 966 papers from 1973 to 2022) and explored with scientometric software. The outcomes highlight that the investigation of benthic foraminiferal response to pollutants started over 50 years ago. Indeed, not only the number of published documents has recently peaked (i.e., 2021 and 2022) but there has been also a growth in the percentages of papers falling within the Decision Sciences category that deals with the application of foraminiferal indices for the EcoQS assessment.

Foraminifera and plastic pollution: Knowledge gaps and research opportunities

Plastic has become one of the most ubiquitous and environmentally threatening sources of pollution in the Anthropocene. Beyond the conspicuous visual impact and physical damages, plastics both carry and release a cocktail of harmful chemicals, such as monomers, additives and persistent organic pollutants. Here we show through a review of the scientific literature dealing with both plastic pollution and benthic foraminifera (Rhizaria), that despite their critical roles in the structure and function of benthic ecosystems, only 0.4% of studies have investigated the effects of micro- and nano-plastics on this group. Consequently, we urge to consider benthic foraminifera in plastic pollution studies via a tentative roadmap that includes (i) the use of their biological, physiological and behavioral responses that may unveil the effects of microplastics and nanoplastics and (ii) the evaluation of the indicative value of foraminiferal species to serve as proxies for the degree of pollution. This appears particularly timely in the context of the development of management strategies to restore coastal ecosystems.

Total and active benthic foraminiferal community and their response to heavy metals revealed by high throughput DNA and RNA sequencing in the Zhejiang coastal waters, East China Sea

Benthic foraminifera, large protists abundant in marine environments, have been widely used as bioindicators of environmental conditions. In this study, high-throughput sequencing based on small subunit rDNA and rRNA amplifications was used to investigate total and active benthic foraminifera community composition and diversity from nineteen and twelve superficial marine sediment samples in the Zhejiang coastal waters, respectively. The results showed that the dominant taxa of total foraminifera changed from Buliminellidae (hyaline) to Saccamminidae (agglutinated) from north to south along the coastal waters of Zhejiang Province. According to our survey, heavy metal contamination was moderate in Zhejiang coastal waters, and the potential ecological risks posed by Cd and Hg were higher. The contamination level of heavy metals at Yueqing Bay was the highest, followed by those at Sanmen Bay and Hangzhou Bay. Cd, Cu and grain size may be key factors affecting the distribution and composition of active foraminiferal communities.

Mitochondrial cytochrome c oxidase subunit I (COI) metabarcoding of Foraminifera communities using taxon-specific primers

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in most marine environments. Molecular methods such as metabarcoding have revealed a high, yet undescribed diversity of Foraminifera. However, so far only one molecular marker, the 18S ribosomal RNA, was available for metabarcoding studies on Foraminifera. Primers that allow amplification of foraminiferal mitochondrial cytochrome oxidase I (COI) and identification of Foraminifera species were recently published. Here we test the performance of these primers for the amplification of whole foraminiferal communities, and compare their performance to that of the highly degenerate LerayXT primers, which amplify the same COI region in a wide range of eukaryotes. We applied metabarcoding to 48 samples taken along three transects spanning a North Sea beach in the Netherlands from dunes to the low tide level, and analysed both sediment samples and meiofauna samples, which contained taxa between 42 µm and 1 mm in body size obtained by decantation from sand samples. We used single-cell metabarcoding (Girard et al., 2022) to generate a COI reference library containing 32 species of Foraminifera, and used this to taxonomically annotate our community metabarcoding data. Our analyses show that the highly degenerate LerayXT primers do not amplify Foraminifera, while the Foraminifera primers are highly Foraminifera- specific, with about 90% of reads assigned to Foraminifera and amplifying taxa from all major groups, i.e., monothalamids, Globothalamea, and Tubothalamea. We identified 176 Foraminifera ASVs and found a change in Foraminifera community composition along the beach transects from high tide to low tide level, and a dominance of single-chambered monothalamid Foraminifera. Our results highlight that COI metabarcoding can be a powerful tool for assessing Foraminiferal communities.

Combining eDNA and morphological approaches to reveal the impacts of long-term discharges of shale gas wastewaters on receiving waters

The potential threats of shale gas wastewater discharges to receiving waters is of great concern. In this study, chemical analyses and biomonitoring were performed three times in a small river that received treated wastewater over a two-year period. The results of chemical analyses showed that the concentrations of chloride, conductivity, barium, and strontium increased at the discharge site, but their concentrations decreased considerably farther downstream (≥500 m). The concentrations of toxic organic compounds (16 US EPA priority polycyclic aromatic hydrocarbons and 6 priority phthalates), trace metals (strontium, arsenic, zinc, copper, chromium, lead, cadmium, nickel, and neodymium), and natural radionuclides (⁴⁰K, ²³⁸U, ²²⁶Ra, and ²³²Th) were comparable to the corresponding background values or did not exhibit obvious accumulation in sediments with continued discharge. Morphological and environmental DNA approaches were used to reveal the potential effects of wastewater discharges on aquatic ecosystems. The results showed that the community structure of benthic invertebrates was not altered by the long-term discharges of shale gas wastewaters. However, the biodiversity indices (richness and Shannon) from the two approaches showed inconsistencies, which were caused by multiple reasons, and that substrates had a strong influence on the morphological biodiversity indices. A multimetric index was proposed to further analyze morphological and environmental DNA data, and the results showed no significant difference between the upstream and downstream sites. Generally, the chemical and biological results both demonstrated that the discharges of shale gas wastewaters had limited impacts on river ecosystems within two years.

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.

Assessing the ecological quality status of the highly polluted Bagnoli area (Tyrrhenian Sea, Italy) using foraminiferal eDNA metabarcoding

Morphology-based benthic foraminifera indices are increasingly used worldwide for biomonitoring the ecological quality of marine sediments. The recent development of foraminiferal eDNA metabarcoding offers a reliable, time-, and cost-effective alternative to morphology-based foraminiferal biomonitoring. However, the practical applications of these new tools are still highly limited. In the present study, we evaluate the response of benthic foraminifera and define the ecological quality status (EcoQS) in the Bagnoli area (Tyrrhenian Sea, Italy) based on a traditional morphology-based approach and eDNA metabarcoding. The geochemical data show that several sites in front of the former industrial plant contain higher concentrations of potentially toxic elements than the effect range median and are characterized by the highest total organic carbon (TOC) content, whereas the distantly located sites can be considered relatively low- to unpolluted. Significant differences (i.e., diversity and assemblage composition) in both morphological and molecular datasets were found between the relatively low- to unpolluted and the most polluted areas. Similarly, the selected ecological indices of both morphological and molecular datasets strikingly and congruently resulted in a clear separation following the environmental stress gradient. The molecular indices (i.e., g-exp(H'_bc), g-Foram AMBI, and g-Foram AMBI-MOTUs) reliably identified poor-to-bad EcoQS in the polluted area in front of the former industrial plant. On the other hand, the Foram-AMBI based on morphology well identified an overall trend but seemed to overestimate the EcoQS if the traditional class boundaries were considered. The congruent and complementary trends between morphological and metabarcoding data observed in the case of the Bagnoli site further support the application of foraminiferal metabarcoding in routine biomonitoring to assess the environmental impacts of heavily polluted marine areas.

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.

Biotic signals associated with benthic impacts of salmon farms from eDNA metabarcoding of sediments

Environmental DNA (eDNA) metabarcoding can rapidly characterize the composition and diversity of benthic communities, thus it has high potential utility for routine assessments of benthic impacts of marine finfish farming. In this study, 126 sediment grab samples from 42 stations were collected at six salmon farms in British Columbia, Canada. Benthic community changes were assessed by both eDNA metabarcoding of metazoans and macrofaunal polychaete surveys. The latter was done by analysing 11,466 individuals using a combination of morphology-based taxonomy and DNA barcoding. Study objectives were to: (i) compare biotic signals associated with benthic impacts of salmon farming in the two data sources, and (ii) identify potential eDNA indicators to facilitate monitoring in Canada. Alpha diversity parameters were consistently reduced near fish cage edge and negatively correlated with pore-water sulphide concentration, with coefficients ranging from -0.62 to -0.48. Although Polychaeta are a common indicator group, the negative correlation with pore-water sulphide concentration was much stronger for Nematoda OTU richness (correlation coefficient: -0.86) than for Polychaeta (correlation coefficient: -0.38). Presence/absence of Capitella generally agreed well between the two methods despite that they differed in the volume of sediments sampled and the molecular marker used. Multiple approaches were used to identify OTUs related to organic enrichment statuses. We demonstrate that eDNA metabarcoding generates biotic signals that could be leveraged for environmental assessment of benthic impacts of fish farms in multiple ways: both alpha diversity and Nematoda OTU richness could be used to assess the spatial extent of impact, and OTUs related to organic enrichment could be used to develop local biotic indices.

Benthic foraminiferal metabarcoding and morphology-based assessment around three offshore gas platforms: Congruence and complementarity

Since the 1960 s, there has been a rapid expansion of drilling activities in the central and northern Adriatic Sea to meet the increasing global energy demand. The discharges of organic and inorganic pollutants, as well as the alteration of the sediment substrate, are among the main impacts associated with these activities. In the present study, we evaluate the response of benthic foraminifera to the activities of three gas platforms in the northwestern Adriatic Sea, with a special focus on the Armida A platform for which extensive geochemical data (organic matter, trace elements, polycyclic aromatic hydrocarbons, other hydrocarbons, and volatile organic compounds) are available. The response to disturbance is assessed by analyzing the foraminiferal diversity using the traditional morphology-based approach and by 18S rDNA-based metabarcoding. The two methods give congruent results, showing relatively lower foraminiferal diversity and higher dominance values at stations closer to the platforms (<50 m). The taxonomic compositions of the morphological and metabarcoding datasets are very different, the latter being dominated by monothalamous, mainly soft-walled species. However, compositional changes consistently occur at 50 m from the platform and can be related to variations in sediment grain-size variation and higher concentrations of Ni, Zn, Ba, hydrocarbons and total organic carbon. Additionally, several morphospecies and Molecular Operational Taxonomic Units (MOTUs) show strong correlations with distance from the platform and with environmental parameters extracted from BIOENV analysis. Some of these MOTUs have the potential to become new bioindicators, complementing the assemblage of hard-shelled foraminiferal species detected through microscopic analyses. The congruence and complementarity between metabarcoding and morphological approaches support the application of foraminiferal metabarcoding in routine biomonitoring surveys as a reliable, time- and cost-effective methodology to assess the environmental impacts of marine industries.

Modeling the ecological status response of rivers to multiple stressors using machine learning: A comparison of environmental DNA metabarcoding and morphological data

Understanding the ecological status response of rivers to multiple stressors is a precondition for river restoration and management. However, this requires the collection of appropriate data, including environmental variables and the status of aquatic organisms, and analysis via a suitable model that captures the nonlinear relationships between ecological status and various stressors. The morphological approach has been the standard data collection method employed for establishing the status of aquatic organisms. However, this approach is very laborious and restricted to a specific set of organisms. Recently, an environmental DNA (eDNA) metabarcoding data approach has been developed that is far more efficient than the morphological approach and potentially applicable to an unlimited set of organisms. However, it remains unclear how well eDNA metabarcoding data reflects the impacts of environmental stressors on aquatic ecosystems compared with morphological data, which is essential for clarifying the potential applications of eDNA metabarcoding data in the ecological monitoring and management of rivers. The present work addresses this issue by modeling organism diversity based on three indices with respect to multiple environmental variables in both the catchment and reach scales. This is done by corresponding support vector machine (SVM) models constructed from eDNA metabarcoding and morphological data on 24 sampling locations in the Taizi River basin, China. According to the mean absolute percent error (MAPE) between the measured diversity index values and the index values predicted by the SVM models, the SVM models constructed from eDNA metabarcoding data (MAPE = 3.87) provide more accurate predictions than the SVM models constructed from morphological data (MAPE = 28.36), revealing that the eDNA metabarcoding data better reflects environmental conditions. In addition, the sensitivity of SVM model predictions of the ecological indices for both catchment-scale and reach-scale stressors is evaluated, and the stressors having the greatest impact on the ecological status of rivers are identified. The results demonstrate that the ecological status of rivers is more sensitive to environmental stressors at the reach scale than to stressors at the catchment scale. Therefore, our study is helpful in exploring the potential applications of eDNA metabarcoding data and SVM modeling in the ecological monitoring and management of rivers.