Multisubstance Indicators Based on Caged Gammarus Bioaccumulation Reveal the Influence of Chemical Contamination on Stream Macroinvertebrate Abundances across France

Identification and organ-specific patterns of expression of two metallothioneins in the sentinel species Gammarus fossarum

Metal pollution is a major concern for aquatic environments. Widespread contamination by various trace metal ions has been described in freshwater streams as well as their subsequent bioaccumulation, potentially leading to toxicity and trophic transfer. Metallothioneins constitute an evolutionary conserved family of low molecular weight, cysteine-rich, metal-chelating proteins, whose known physiological functions are the maintenance of the homeostasis of essential metals, the detoxification of non-essential metals, and the protection against oxidative stress and free radicals. In this study, we identified two metallothionein-coding transcripts, mt1 and mt2, in the transcriptome of the amphipod Gammarus fossarum, a sentinel species widely used to assess the quality of watersheds. For the first time, we investigated the organ-specific patterns of expression of these two mt transcripts at the individual level in the gills and the caeca of this small crustacean. In silico analysis and experimental exposures to environmentally relevant concentrations of cadmium, zinc and silver showed that G. fossarum mt1 induction is stronger after Cd exposure compared to the other tested metals. G. fossarum mt1 was more significantly induced in the caeca than in the gills of exposed organisms for any metal exposure, while G. fossarum mt2 was, at least at the individual level, more inducible in the gills than in the caeca of G. fossarum exposed to Cd and Zn. Our results provide new genetic resources that will help to improve the understanding of metal homeostasis in this sentinel species.

Contrasted Antibiotics and Pesticides Occurrence in Fish Exposed In Situ to Urban Effluents: A 20-Day Caging Experiment

Urban freshwater ecosystems receive a wide array of organic pollutants through wastewater-treatment plant (WWTP) discharges and agricultural runoff. Evaluating the fate and effects of antibiotics and pesticides can be a challenging task, especially the effects on freshwater vertebrates because of their abilities to metabolize and excrete these chemicals and because of their high mobility and escape behavior when exposed to stressful environmental conditions. In the present study, 37 wild gudgeons (Gobio gobio) were caged for a period of up to 20 days, upstream and downstream of a WWTP effluent discharge in the Orge River (a tributary of the Seine River, France). Levels of pesticides and antibiotics in fish muscles were monitored weekly and compared with environmental contamination (water and sediments). Our results highlighted a slight bioaccumulation of pesticides in the gudgeon muscles at the downstream site after 20 days of exposure. Concerning antibiotics, ofloxacin was the most detected compound in fish muscles (85% of occurrence) and ranged from undetectable to 8 ng g-1  dry weight. Antibiotic levels in fish muscle were not higher at the downstream site and did not increase with exposure duration, despite high levels in the water (up to 29 times greater than upstream). Potential ecotoxicological effects were also evaluated: Body condition did not differ between the caging location and exposure time. Three oxidative status markers in the fish livers showed significant shifts after 14 days of caging. Our results suggest a high clearance rate of antibiotics and, to a lesser extent, of pesticides in wild gudgeons, which could be explained by changes in xenobiotic metabolism with pollutant exposure. Environ Toxicol Chem 2024;00:1-11. © 2023 SETAC.

Acclimation and transgenerational plasticity support increased cadmium tolerance in Gammarus populations exposed to natural metal contamination in headwater streams

Considering long-term population effects of chronic exposure to contaminants remains limited in ecological risk assessment. Field evidence that multigenerational exposure influences organisms' sensitivity is still scarce, and mechanisms have yet to be elucidated in the environmental context. This study focuses on the crustacean Gammarus fossarum, for which an increased tolerance to cadmium (Cd) has previously been reported in a naturally low-contaminated headwater stream. Our objectives were to investigate whether Cd tolerance is a common phenomenon in headwater populations, and to elucidate the nature of the tolerance and its intergenerational transmission. For this, we carried out an in-depth in situ characterization of Cd exposure (gammarids' caging) and levels of tolerance in nine populations on a regional scale, as well as laboratory maintenance and cross-breeding of contaminated and uncontaminated populations. Acute tolerance levels correlate positively with bioavailable Cd contamination levels among streams. The contaminated and non-contaminated populations differ about two-fold in sensitivity to Cd. Tolerance was found in all age classes of contaminated populations, it can be transiently lost during the year, and it was transmissible to offspring. In addition, tolerance levels dropped significantly when organisms were transferred to a Cd-free environment for two months. These organisms also ceased producing tolerant offspring, confirming a non-genetic transmission of Cd tolerance between generations. These findings support that Cd tolerance corresponds to non-genetic acclimation combined with transgenerational plasticity. Moreover, cross-breeding revealed that tolerance transmission to offspring is not limited to maternal effect. We suggest epigenetics as a plausible mechanism for the plasticity of Cd sensitivity observed in the field. Our results therefore highlight the neglected role of plasticity and non-genetic transmission of modified sensitivities during the long-term exposure of natural populations to environmental contamination.

Avoidance behaviour of aquatic macroinvertebrates for real-time detection of micropollutant surge in wastewater effluents

Micropollutants are regularly detected at the outlets of wastewater treatment plants (WWTPs). Across urban and industrial WWTPs, monitoring directives only require assessment for a handful of chemicals via sampling methods that fail to capture the temporal variability in micropollutant discharge. In this study, we develop a biotest for real-time on-line monitoring of micropollutant discharge dynamics in WWTPs effluents. The selected biomonitoring device ToxMate uses videotracking of invertebrate movement, which was used to deduce avoidance behaviour of the amphipod Gammarus fossarum. Organism conditioning was set up to induce a state of minimal locomotor activity in basal conditions to maximise avoidance signal sensitivity to micropollutant spikes. We showed that with a standardised protocol, it was possible to minimise both overall movement and sensitivity to physio-chemical variations typical to WWTP effluents, as well as capture the spikes of two micropollutants upon exposure (copper and methomyl). Spikes in avoidance behaviour were consistently seen for the two chemicals, as well as a strong correlation between avoidance intensity and spiked concentration. A two-year effluent monitoring case study also illustrates how this biomonitoring method is suitable for real-time on-site monitoring, and shows a promising non-targeted approach for characterising complex micropollutant discharge variability at WWTP effluents, which today remains poorly understood.

The retinoid metabolism of Gammarus fossarum is disrupted by exogenous all-trans retinoic acid, citral, and methoprene but not by the technical formulation of glyphosate

Over the last decade, fluctuations of retinoids (RETs), also known as vitamin A and derivatives, have proved to be useful biomarkers to assess the environmental chemical pressure on a wide variety of non-target vertebrates. This use of RET-based biomarkers is of particular interest in the non-target sentinel species Gammarus fossarum in which RETs were shown to influence crucial physiological functions. To study and probe this metabolism in this crustacean model, a UHPLC-MS/MS method was developed to 1) identify and 2) monitor several endogenous RETs in unexposed females throughout their reproductive cycle. Then, females were exposed in controlled conditions to exogenous all-trans retinoic acid (atRA) and citral (CIT), a RA synthesis inhibitor, to simulate an excess or deficiency in RA. Perturbation of vitamin A metabolism by pesticides was further studied in response to methoprene (MET), a juvenile hormone analog as well as glyphosate (GLY). The developed method allowed, for the first time in this model, the identification of RA metabolites (all-trans 4-oxo and 13-cis 4-oxo RA), RA isomers (all-trans and 13-cis RA) as well as retinaldehyde (RALD) isomers (all-trans, 11-cis, and 13-cis RALD) and showed two distinct phases in the reproductive cycle. Retinoic acid successfully increased the tissular concentration of both RA isomers and CIT proved to be efficient at perturbating the conversion from RALD to RA. Methoprene perturbed the ratios between RA isomers whereas GLY had no observed effects on the RET system of G. fossarum females. We were able to discriminate different dynamics of RET perturbations by morphogens (atRA or CIT) or MET which highlights the plausible mediation of RETs in MET-induced disorders. Ultimately, our study shows that RETs are influenced by exposure to MET and strengthen their potential to assess aquatic ecosystem chemical status.

Identifying the impact of toxicity on stream macroinvertebrate communities in a multi-stressor context based on national ecological and ecotoxicological monitoring databases

In situ bioassays are used to measure the harmful effects induced by mixtures of toxic chemicals in watercourses. In France, national-scale biomonitoring data are available including invertebrate surveys and in-field chemical toxicity measures with caged gammarids to assess environmental toxicity of mixtures of chemicals. The main objective of our study is to present a proof-of-concept approach identifying possible links between in-field chemical toxicity, stressors and the ecological status. We used two active biomonitoring databases comprising lethal toxicity (222 in situ measures of gammarid mortality) and sublethal toxicity (101 in situ measures of feeding inhibition). We measured the ecological status of each active biomonitoring site using the I₂M₂ metric (macroinvertebrate-based multimetric index), accounted for known stressors of nutrients and organic matter, hydromorphology and chemical toxicity. We observed a negative relationship between stressors (hydromorphology, nutrients and organic matter, and chemical toxicity) and the good ecological status. This relationship was aggravated in watercourses where toxicity indicators were degraded. We validated this hypothesis for instance with nutrients and organic matter like nitrates or hydromorphological conditions like percentage of vegetation on banks. Future international assesments concerning the role of in-field toxic pollution on the ecological status in a multi-stressor context are now possible via the current methodology.

Influence of the exposure concentration of dissolved cadmium on its organotropism, toxicokinetic and fate in Gammarus fossarum

Information on the relationship between the exposure concentrations of metals and their biodistribution among organs remained scarce in invertebrates. The objective of this study was to investigate the effects of Cd concentration on the organotropism, toxico-kinetic and fate of this metal in different organs of gammarids exposed to dissolved radioisotope ¹⁰⁹Cd. Gammarids male were exposed for 7 days to three environmental Cd concentrations (i.e. 4, 52 and 350 ng.L^-1) before being placed in depuration conditions (i.e. uncontaminated water). At several sampling times, Cd concentrations were determined by ¹⁰⁹Cd γ-counting in water, caeca, cephalon, gills, intestine and remaining tissues. Bioconcentration Factors (BCF) and Cd relative proportions in organs were calculated to assess the exposure concentration effect on the bioaccumulation capacities. The dependance of the organ-specific kinetic parameters to Cd water concentrations were estimated by fitting nested one-compartment toxico-kinetic (TK) models to both the accumulation and depuration data, by Bayesian inference. Then, for each Cd concentrations, the metal exchanges among organs over time were formalized by a multi-compartments TK model fitted to all organ data simultaneously. Our results highlighted that, at the end of the exposure phase, BCF and Cd relative proportions, in each organ, were not significantly modulated by water concentrations. Kinetically, Cd accumulation rates in all organs (except intestines) were depended on the exposure concentration, but not the elimination rates. The in vivo management of Cd (i.e. metal exchanges among organs) remained qualitatively unchanged according to exposure concentration. All these results also highlighted key role of that organs in the management of Cd: bioconcentration by caeca, storage by gills and main entry pathway by intestine. This study shows the interest of implementing TK approaches to test the effect of environmental factors on bioaccumulation, inter-organ exchanges and fate of contaminants in invertebrate body to enhance the understanding of the toxicity risk.

Multiple reaction monitoring mass spectrometry for the discovery of environmentally modulated proteins in an aquatic invertebrate sentinel species, Gammarus fossarum

Multiple reaction monitoring (MRM) mass spectrometry is emerging as a relevant tool for measuring customized molecular markers in freshwater sentinel species. While this technique is typically used for the validation of protein molecular markers preselected from shotgun experiments, recent gains of MRM multiplexing capacity offer new possibilities to conduct large-scale screening of animal proteomes. By combining the strength of active biomonitoring strategies and MRM technologies, this study aims to propose a new strategy for the discovery of candidate proteins that respond to environmental variability. For this purpose, 249 peptides derived from 147 proteins were monitored by MRM in 273 male gammarids caged in 56 environmental sites, representative of the diversity of French water bodies. A methodology is here proposed to identify a set of customized housekeeping peptides (HKPs) used to correct analytical batch effects and allow proper comparison of peptide levels in gammarids. A comparative analysis performed on HKPs-normalized data resulted in the identification of peptides highly modulated in the environment and derived from proteins likely involved in the environmental stress response. Overall, this study proposes a breakthrough approach to screen and identify potential proteins responding to relevant environmental conditions in sentinel species.

Metal bioavailable contamination engages richness decline, species turnover but unchanged functional diversity of stream macroinvertebrates at the scale of a French region

Freshwater ecosystems are the main source of water for sustaining life on earth, and the biodiversity they support is the main source of valuable goods and services for human populations. Despite growing recognition of the impairment of freshwater ecosystems by micropollutant contamination, different conceptual and methodological considerations can newly be addressed to improve our understanding of the ecological impact into these ecosystems. Here, we originally combined in situ ecotoxicology and community ecology concepts to unveil the mechanisms structuring macroinvertebrate communities along a regional contamination gradient. The novelty of our study lies in the use of an innovative biomonitoring approach (measurement of metal contents in caged crustaceans) allowing to quantify and compare on a regional scale the levels of bioavailable metal contamination to which stream communities are exposed. We were hence able to identify 23 streams presenting a significant gradient of bioavailable metal contamination within the same catchment area in the South West of France, from which we also obtained data on the composition of resident macroinvertebrate communities. Analyses of structural and functional integrity of communities revealed an unexpected decoupling between taxonomic and functional diversity of communities in response to bioavailable metal contamination. We show that despite the negative impact of bioavailable metal contamination exposure on taxonomic diversity (with an average species loss of 17% in contaminated streams), functional diversity is maintained through a process of non-random species replacement by functional redundant species at the regional scale. Such unanticipated findings call for a deeper characterization of metal-tolerant communities' ability to cope with environmental variability in multi-stressed ecosystems.

In-situ biological effects, bioaccumulation, and multi-media distribution of organic contaminants in a shallow lake

Anthropogenic activity-impacted aquatic environment contains a complex mixture of contaminants, and ecological risk assessment solely based on chemical analysis is insufficient and biological assessment is required. However, traditional assessment heavily relies on laboratory bioassays, which may cause uncertainty due to inevitable laboratory-related artifact. A self-designed in-situ bioassay system was successfully applied to simultaneously evaluate water and sediment toxicity by co-exposure of two native species, Chinese rare minnows (Gobiocypris rarus) and Asian clams (Corbicula fluminea) in Tai Lake Basin, China. In-situ exposure caused pronouncedly lethal and sublethal effects (i.e., metabolic and oxidative stress, neurotoxicity, reproductive toxicity) on both fish and clams. Meanwhile, multi-media distribution of organic contaminants in water-sediment-biota system was analyzed. Besides hydrophobicity, metabolism was recognized as an influential factor on phase distribution of contaminants in water-fish and sediment-clam systems. Traditional hazard quotient (HQ) method based on environmental concentrations of 98 contaminants showed bias in risk assessment. Instead, a weight of evidence method by integrating three lines of evidence, including in-situ survival, enhanced integrated biomarker response values and environmental concentrations, successfully differentiate high- and moderate-risk sites in the shallow lakes. The present study incorporated in-situ bioassays into risk assessment using a weight of evidence approach, which reduced uncertainty in decision-making.

Quantification of multi-scale links of anthropogenic pressures with PAH and PCB bioavailable contamination in French freshwaters

Aquatic ecosystems are exposed to multiple environmental pressures including chemical contamination. Polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) are persistent organic pollutants (POPs) known as preoccupying substances for the environment. Active biomonitoring (ABM) is a surveillance method for polluted aquatic ecosystems measuring bioavailable contamination. In this work, the aim was to quantify the total links between environmental pressures and bioavailable contamination (for PAHs and PCBs) at the French national scale. Based on 245 sites experimented by ABM from 2017 to 2019, environmental pressures (anthropogenic pressures and environmental parameters) were defined (point source landfill density, point source urban density, point source industry density, point source road density, nonpoint source industry density, nonpoint source road density, nonpoint source urban density, nutrients and organic matter, slope, dams, straightness, coarse sediment, summer precipitation, hydrographic network density and watershed size) and characterized by one or a combination of measures called stressor indicators. The links between environmental pressures and bioavailable POPs contamination (ABM measure) at a large spatial scale were defined and quantified via structural equation modeling. Point source urban density, nutrients and organic matter, summer precipitation, straightness and point source industry density are correlated positively with PAH bioavailable contamination. In contrast, nonpoint source urban density, nonpoint source industry density, nonpoint source road density and watershed size are positively correlated with PCB bioavailable contamination. The dominant pressures linked to PAHs and PCBs were different, respectively local and large-scale pressures were linked to PAH bioavailable contamination, and only large-scale pressures were linked to PCB bioavailable contamination.

Co-expression network analysis identifies novel molecular pathways associated with cadmium and pyriproxyfen testicular toxicity in Gammarus fossarum

Omics approaches are continuously providing new clues on the mechanisms of action of contaminants in species of environmental relevance, contributing to the emergence of molecular ecotoxicology. Co-expression network approaches represent a suitable methodological framework for studying the rich content of omics datasets. This study aimed to find evidence of key pathways and proteins related to the testicular toxicity in the sentinel crustacean species Gammarus fossarum exposed to endocrine disruptors using a weighted protein co-expression network analysis. From a shotgun proteomics dataset of male gonads of G. fossarum organisms exposed to cadmium (Cd), pyriproxyfen (Pyr) and methoxyfenozide (Met) in laboratory conditions, four distinct modules were identified as significantly correlated to contaminants' exposure. Protein set enrichment analysis identified modules involved in cytoskeleton organization and oxidative stress response associated with the Cd exposure. The module associated with Pyr exposure was associated with endoplasmic reticulum stress (ER) response, and the module correlated with Met exposure was characterized by a significant proportion of amphipod-restricted proteins whose functions are still not characterized. Our results show that co-expression networks are efficient and adapted tools to identify new potential mode of actions from environmental sentinel species, such as G. fossarum, using a proteogenomic approach, even without an annotated genome.

How to quantify the links between bioavailable contamination in watercourses and pressures of anthropogenic land cover, contamination sources and hydromorphology at multiple scales?

Active biomonitoring permits the quantification of biological exposure to chemicals through measurements of bioavailable concentrations in biota and biological markers of toxicity in organisms. It enables respective comparison of the levels of contamination between sites and sampling campaigns. Caged gammarids are recently proposed as relevant probes for measuring bioavailable contamination in freshwater systems. The purpose of the present study was to develop a multi-pressure and multiscale approach, considering metallic contamination levels (from data based on active biomonitoring) as a response to pressures (combination of individual stressors). These pressures were anthropogenic land cover, industry density, wastewater treatment plant density, pressures on stream hydromorphological functioning, riverside vegetation and bioavailability factors. A dataset combining active biomonitoring and potentially related pressures was established at the French national scale, with 196 samplings from 2009 to 2016. The links between pressures and metallic contamination were defined and modelled via structural equation modeling (more specifically partial least squares - path modeling). The model enabled the understanding of the respective influences of pressures on metallic bioconcentration in caged sentinel organisms. Beyond validating the local influence of industries and wastewater treatment plants on metallic contamination, this model showed a complementary effect of driving forces of anthropogenic land cover (leading to human activities). It also quantified a significant influence of pressures on stream hydromorphological functioning, presence of vegetation and physico-chemical parameters on metal bioconcentration. This hierarchical multi-pressure approach could serve as a concept on how pressures and contamination (assessed by active biomonitoring) can be connected. Its future application will enable better understanding of environmental pressures leading to contamination in freshwater ecosystems.

High-multiplexed monitoring of protein biomarkers in the sentinel Gammarus fossarum by targeted scout-MRM assay, a new vision for ecotoxicoproteomics

Ecotoxicoproteomics employs mass spectrometry-based approaches centered on proteins of sentinel organisms to assess for instance, chemical toxicity in fresh water. In this study, we combined proteogenomics experiments and a novel targeted proteomics approach free from retention time scheduling called Scout-MRM. This methodology will enable the measurement of simultaneously changes in the relative abundance of multiple proteins involved in key physiological processes and potentially impacted by contaminants in the freshwater sentinel Gammarus fossarum. The development and validation of the assay were performed to target 157 protein biomarkers of this non-model organism. We carefully chose and validated the transitions to monitor using conventional parameters (linearity, repeatability, LOD, LOQ). Finally, the potential of the methodology is illustrated by measuring 277-peptide-plex assay (831 transitions) in sentinel animals exposed in natura to different agricultural sites potentially exposed to pesticide contamination. Multivariate data analyses highlighted the modulation of several key proteins involved in feeding and molting. This multiplex-targeted proteomics assay paves the way for the discovery and the use of a large panel of novel protein biomarkers in emergent ecotoxicological models for environmental monitoring in the future. BIOLOGICAL SIGNIFICANCE: The study contributed to the development of Scout-MRM for the high-throughput quantitation of a large panel of proteins in the Gammarus fossarum freshwater sentinel. Increasing the number of markers in ecotoxicoproteomics is of most interest to assess the impact of pollutants in freshwater organisms. The development and validation of the assay enabled the monitoring of a large panel of reporter peptides of exposed gammarids. To illustrate the applicability of the methodology, animals from different agricultural sites were analysed. The application of the assay highlighted the modulation of some biomarker proteins involved in key physiological pathways, such as molting, feeding and general stress response. Increasing multiplexing capabilities and field test will provide the development of diagnostic protein biomarkers for emergent ecotoxicological models in future environmental biomonitoring programs.

A "Population Dynamics" Perspective on the Delayed Life-History Effects of Environmental Contaminations: An Illustration with a Preliminary Study of Cadmium Transgenerational Effects over Three Generations in the Crustacean Gammarus

We explore the delayed consequences of parental exposure to environmentally relevant cadmium concentrations on the life-history traits throughout generations of the freshwater crustacean Gammarus fossarum. We report the preliminary results obtained during a challenging one-year laboratory experiment in this environmental species and propose the use of population modeling to interpret the changes in offspring life-history traits regarding their potential demographic impacts. The main outcome of this first long-term transgenerational assay is that the exposure of spawners during a single gametogenesis cycle (3 weeks) could result in severe cascading effects on the life-history traits along three unexposed offspring generations (one year). Indeed, we observed a decrease in F1 reproductive success, an early onset of F2 offspring puberty with reduced investment in egg yolk reserves, and finally a decrease in the growth rate of F3 juveniles. However, the analysis of these major transgenerational effects by means of a Lefkovitch matrix population model revealed only weak demographic impacts. Population compensatory processes mitigating the demographic consequences of parental exposure seem to drive the modification of life-history traits in offspring generations. This exploratory study sheds light on the role of population mechanisms involved in the demographic regulation of the delayed effects of environmental toxicity in wild populations.

In Situ Reproductive Bioassay with Caged Gammarus fossarum (Crustacea): Part 2-Evaluating the Relevance of Using a Molt Cycle Temperature-Dependent Model as a Reference to Assess Toxicity in Freshwater Monitoring

Active biomonitoring approaches are now recognized as relevant for monitoring water contamination and toxicity. Nevertheless, due to the confounding influence of variable and uncontrolled environmental conditions such as temperature, biological markers measured on transplanted individuals to assess water quality are difficult to interpret. The purpose of the present study is to propose a methodology for adapting a laboratory test of chronic sublethal toxicity based on the molting cycle of Gammarus fossarum to in situ assays. To this end, we 1) adapted the molt cycle temperature-dependent model developed in Part 1 (Chaumot et al. 2020, this issue) to the fluctuating temperatures measured in the field; 2) assessed the predictive power of our approach as a "reference value" from gammarids caged in 9 nonimpacted sites at different seasons; and 3) tested the relevance of our tool to interpret in situ reproductive bioassays from 5 upstream/downstream studies and a large-scale deployment in 12 sites. Our approach based on modeling the progress of gammarid molting cycle as a function of temperature appeared to be a relevant and robust tool for interpreting in situ observations in different environmental contexts in time and space. By avoiding using a "reference" or upstream situation as a baseline from which water quality could be assessed, this approach provides a real added value to water quality diagnosis in biomonitoring programs. Environ Toxicol Chem 2020;39:678-691. © 2019 SETAC.