Multiple in-stream stressors degrade biological assemblages in five U.S. regions

Using DNA metabarcoding to characterize national scale diatom-environment relationships and to develop indicators in streams and rivers of the United States

Recent advancements in DNA techniques, metabarcoding, and bioinformatics could help expand the use of benthic diatoms in monitoring and assessment programs by providing relatively quick and increasingly cost-effective ways to quantify diatom diversity in environmental samples. However, such applications of DNA-based approaches are relatively new, and in the United States, unknowns regarding their applications at large scales exist because only a few small-scale studies have been done. Here, we present results from the first nationwide survey to use DNA metabarcoding (rbcL) of benthic diatoms, which were collected from 1788 streams and rivers across nine ecoregions spanning the conterminous USA. At the national scale, we found that diatom assemblage structure (1) was strongly associated with total phosphorus and total nitrogen concentrations, conductivity, and pH and (2) had clear patterns that corresponded with differences in these variables among the nine ecoregions. These four variables were strong predictors of diatom assemblage structure in ecoregion-specific analyses, but our results also showed that diatom-environment relationships, the importance of environmental variables, and the ranges of these variables within which assemblage changes occurred differed among ecoregions. To further examine how assemblage data could be used for biomonitoring purposes, we used indicator species analysis to identify ecoregion-specific taxa that decreased or increased along each environmental gradient, and we used their relative abundances of gene reads in samples as metrics. These metrics were strongly correlated with their corresponding variable of interest (e.g., low phosphorus diatoms with total phosphorus concentrations), and generalized additive models showed how their relationships compared among ecoregions. These large-scale national patterns and nine sets of ecoregional results demonstrated that diatom DNA metabarcoding is a robust approach that could be useful to monitoring and assessment programs spanning the variety of conditions that exist throughout the conterminous United States.

Linking wastewater treatment plant effluents to water quality and hydrology: Effects of multiple stressors on fish communities

Wastewater treatment plants (WWTPs) are essential for maintaining a good water quality of surface waters. However, WWTPs are also associated with water quality deterioration and hydro-morphological alteration. Riverine communities respond to these stressors with changes in their community structure, abundance and diversity. In this study, we used a dataset of 94 monitoring sites across North Rhine-Westphalia, Germany to investigate the influence of WWTPs on the water quality and hydro-morphological quality in river sections downstream of WWTP effluents. More specifically, we analyzed the effects of the percentage of WWTP effluents (in relation to median base flow) on four stressor groups (physico-chemistry, micropollutants, hydrological and morphological alteration) using Linear Mixed Models (LMM). Furthermore, we assessed the impact of a selection of twelve ecologically relevant stressor variables reflecting water quality deterioration and hydro-morphological alteration on reference fish communities using Canonical Correspondence Analysis (CCA). The percentage of WWTP effluents was correlated with water quality, especially with toxic units of a wide range of pharmaceuticals including diclofenac, venlafaxine and sulfamethoxazole (R² up to 0.54) as well as specific pesticides (e.g., terbutryn: R² = 0.33). The correlation of percent WWTP effluents with hydro-morphological alteration was weaker and most pronounced for the frequency of high flow (R² = 0.24) and flow variability (R² = 0.19). About 40 % of the variance in the fish community structure were explained by 12 stressor variables in the CCA models. Water quality and hydrological, but not morphological stressors showed strong albeit highly variable effects on individual fish species. The results indicate that water quality degradation and hydrological alteration are important factors determining the ecological status of fish communities. In this context, WWTP effluents can impose relevant point sources of pollution that affect water quality but also cause alterations of the hydrological regime. Further management measures addressing both stressor groups are needed to improve the ecological status.

Predicting macroinvertebrate average score per taxon (ASPT) at water quality monitoring sites in Japanese rivers

Biomonitoring with bioindicators such as river macroinvertebrates is fundamental for assessing the status of freshwater ecosystems. In Japan, water quality and biomonitoring surveys are conducted separately, leading to a lack of nationwide information on their relationships and the biological status of water quality monitoring (WQM) sites. To understand the biological status of WQM sites across Japan, we developed a multiple linear regression model to estimate the average score per taxon (ASPT) using river macroinvertebrate data surveyed at a total of 237 "aligned" sites based on the co-occurrence of biomonitoring and WQM sites. The resulting regression model with eight predictors, such as biological oxygen demand, the proportion of urban areas in the catchment, could predict ASPT with reasonable accuracy (e.g., an error of ±1 for 96% of the aligned data). Using this model, we estimated ASPT values at 2925 WQM sites in rivers nationwide, categorizing them into four levels of river environment quality: "very good" (29% of WQM sites), "good" (50%), "fairly good" (14%), and "not good" (8%). Furthermore, we observed statistically significant correlations (p < 0.05; 0.4 ≤ r ≤ 0.7) between ASPT and all eight macroinvertebrate metrics examined, such as mayfly and stonefly richness, providing ecological implications of changes in ASPT.

Multiple lines of evidence point to pesticides as stressors affecting invertebrate communities in small streams in five United States regions

Multistressor studies were performed in five regions of the United States to assess the role of pesticides as stressors affecting invertebrate communities in wadable streams. Pesticides and other chemical and physical stressors were measured in 75 to 99 streams per region for 4 weeks, after which invertebrate communities were surveyed (435 total sites). Pesticides were sampled weekly in filtered water, and once in bed sediment. The role of pesticides as a stressor to invertebrate communities was assessed by evaluating multiple lines of evidence: toxicity predictions based on measured pesticide concentrations, multivariate models and other statistical analyses, and previously published mesocosm experiments. Toxicity predictions using benchmarks and species sensitivity distributions and statistical correlations suggested that pesticides were present at high enough concentrations to adversely affect invertebrate communities at the regional scale. Two undirected techniques-boosted regression tree models and distance-based linear models-identified which pesticides were predictors of (respectively) invertebrate metrics and community composition. To put insecticides in context with known, influential covariates of invertebrate response, generalized additive models were used to identify which individual pesticide(s) were important predictors of invertebrate community condition in each region, after accounting for natural covariates. Four insecticides were identified as stressors to invertebrate communities at the regional scale: bifenthrin, chlordane, fipronil and its degradates, and imidacloprid. Fipronil was particularly important in the Southeast region, and imidacloprid, bifenthrin, and chlordane were important in multiple regions. For imidacloprid, bifenthrin, and fipronil, toxicity predictions were supported by mesocosm experiments that demonstrated adverse effects on naïve aquatic communities when dosed under controlled conditions. These multiple lines of evidence do not prove causality-which is challenging in the field under multistressor conditions-but they make a strong case for the role of insecticides as stressors adversely affecting invertebrate communities in streams within the five sampled regions.

Improving ecosystem health in highly altered river basins: a generalized framework and its application to the Mississippi-Atchafalaya River Basin

Continued large-scale public investment in declining ecosystems depends on demonstrations of "success". While the public conception of "success" often focuses on restoration to a pre-disturbance condition, the scientific community is more likely to measure success in terms of improved ecosystem health. Using a combination of literature review, workshops and expert solicitation we propose a generalized framework to improve ecosystem health in highly altered river basins by reducing ecosystem stressors, enhancing ecosystem processes and increasing ecosystem resilience. We illustrate the use of this framework in the Mississippi-Atchafalaya River Basin (MARB) of the central United States (U.S.), by (i) identifying key stressors related to human activities, and (ii) creating a conceptual ecosystem model relating those stressors to effects on ecosystem structure and processes. As a result of our analysis, we identify a set of landscape-level indicators of ecosystem health, emphasizing leading indicators of stressor removal (e.g., reduced anthropogenic nutrient inputs), increased ecosystem function (e.g., increased water storage in the landscape) and increased resilience (e.g., changes in the percentage of perennial vegetative cover). We suggest that by including these indicators, along with lagging indicators such as direct measurements of water quality, stakeholders will be better able to assess the effectiveness of management actions. For example, if both leading and lagging indicators show improvement over time, then management actions are on track to attain desired ecosystem condition. If, however, leading indicators are not improving or even declining, then fundamental challenges to ecosystem health remain to be addressed and failure to address these will ultimately lead to declines in lagging indicators such as water quality. Although our model and indicators are specific to the MARB, we believe that the generalized framework and the process of model and indicator development will be valuable in an array of altered river basins.

A synthesis of anthropogenic stress effects on emergence-mediated aquatic-terrestrial linkages and riparian food webs

Anthropogenic stress alters the linkage between aquatic and terrestrial ecosystems in various ways. Here, we review the contemporary literature on how alterations in aquatic systems through environmental pollution, invasive species and hydromorphological changes carry-over to terrestrial ecosystems and the food webs therein. We consider both the aquatic insect emergence and flooding as pathways through which stressors can propagate from the aquatic to the terrestrial system. We specifically synthesize and contextualize results on the roles of pollutants in the emergence pathway and their top-down consequences. Our review revealed that the emergence and flooding pathway are only considered in isolation and that the overall effects of invasive species or pollutants on food webs at the water-land interface require further attention. While very few recent studies looked at invasive species, a larger number of studies focused on metal transfer compared to pesticides, pharmaceuticals or PCBs, and multiple stress studies up to now left aquatic-terrestrial linkages unconsidered. Recent research on pollutants and emergence used aquatic-terrestrial mesocosms to elucidate the effects of aquatic stressors such as the mosquito control agent Bti, metals or pesticides to understand the effects on riparian spiders. Quality parameters, such as the structural and functional composition of emergent insect communities, the fatty acid profiles, yet also the composition of pollutants transferred to land prove to be important for the effects on riparian spiders. Process-based models including quality of emergence are useful to predict the resulting top-down directed food web effects in the terrestrial recipient ecosystem. In conclusion, we present and recommend a combination of empirical and modelling approaches in order to understand the complexity of aquatic-terrestrial stressor propagation and its spatial and temporal variation.

Water quality deterioration remains a major stressor for macroinvertebrate, diatom and fish communities in German rivers

About 60 % of Europe's rivers fail to meet ecological quality standards derived from biological criteria. The causes are manifold, but recent reports suggest a dominant role of hydro-morphological and water quality-related stressors. Yet, in particular micropollutants and hydrological stressors often tend to be underrepresented in multiple-stressor studies. Using monitoring data from four Federal States in Germany, this study investigated the effects of 19 stressor variables from six stressor groups (nutrients, salt ions, dissolved oxygen/water temperature, mixture toxicity of 51 micropollutants, hydrological alteration and morphological habitat quality) on three biological assemblages (fishes, macroinvertebrates, benthic diatoms). Biological effects were analyzed for 35 community metrics and quantified using Random Forest (RF) analyses to put the stressor groups into a hierarchical context. To compare metric responses, metrics were grouped into categories reflecting important characteristics of biological communities, such as sensitivity, functional traits, diversity and community composition as well as composite indices that integrate several metrics into one single index (e.g., ecological quality class). Water quality-related stressors - but not micropollutants - turned out to dominate the responses of all assemblages. In contrast, the effects of hydro-morphological stressors were less pronounced and stronger for hydrological stressors than for morphological stressors. Explained variances of RF models ranged 23-64 % for macroinvertebrates, 16-40 % for benthic diatoms and 18-48 % for fishes. Despite a high variability of responses across assemblages and stressor groups, sensitivity metrics tended to reveal stronger responses to individual stressors and a higher explained variance in RF models than composite indices. The results of this study suggest that (physico-chemical) water quality deterioration continues to impact biological assemblages in many German rivers, despite the extensive progress in wastewater treatment during the past decades. To detect water quality deterioration, monitoring schemes need to target relevant physico-chemical stressors and micropollutants. Furthermore, monitoring needs to integrate measures of hydrological alteration (e.g., flow magnitude and dynamics). At present, hydro-morphological surveys rarely address the degree of hydrological alteration. In order to achieve a good ecological status, river restoration and management needs to address both water quality-related and hydro-morphological stressors. Restricting analyses to just one single organism group (e.g., macroinvertebrates) or only selected metrics (e.g., ecological quality class) may hamper stressor identification and its hierarchical classification and, thus may mislead river management.

Chemical Mixtures and Multiple Stressors: Same but Different?

Ecosystems are strongly influenced by multiple anthropogenic stressors, including a wide range of chemicals and their mixtures. Studies on the effects of multiple stressors have largely focussed on nonchemical stressors, whereas studies on chemical mixtures have largely ignored other stressors. However, both research areas face similar challenges and require similar tools and methods to predict the joint effects of chemicals or nonchemical stressors, and frameworks to integrate multiple chemical and nonchemical stressors are missing. We provide an overview of the research paradigms, tools, and methods commonly used in multiple stressor and chemical mixture research and discuss potential domains of cross-fertilization and joint challenges. First, we compare the general paradigms of ecotoxicology and (applied) ecology to explain the historical divide. Subsequently, we compare methods and approaches for the identification of interactions, stressor characterization, and designing experiments. We suggest that both multiple stressor and chemical mixture research are too focused on interactions and would benefit from integration regarding null model selection. Stressor characterization is typically more costly for chemical mixtures. While for chemical mixtures comprehensive classification systems at suborganismal level have been developed, recent classification systems for multiple stressors account for environmental context. Both research areas suffer from rather simplified experimental designs that focus on only a limited number of stressors, chemicals, and treatments. We discuss concepts that can guide more realistic designs capturing spatiotemporal stressor dynamics. We suggest that process-based and data-driven models are particularly promising to tackle the challenge of prediction of effects of chemical mixtures and nonchemical stressors on (meta-)communities and (meta-)food webs. We propose a framework to integrate the assessment of effects for multiple stressors and chemical mixtures. Environ Toxicol Chem 2023;42:1915-1936. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Effects of multiple stressors on benthic invertebrates using Water Framework Directive monitoring data

Multiple stressors affect freshwater systems and cause a deficient ecological status according to the European Water Framework Directive (WFD). To select effective mitigation measures and improve the ecological status, knowledge on the stressor hierarchy and individual and joined effects is necessary. However, compared to common stressors like nutrient enrichment and morphological degradation, the relative importance of micropollutants such as pesticides and pharmaceuticals is largely unaddressed. We used WFD monitoring data from Saxony (Germany) to investigate the importance of 85 environmental variables (including 34 micropollutants) for 18 benthic invertebrate metrics at 108 sites. The environmental variables were assigned to five groups (natural factors, nutrient enrichment, metals, micropollutants and morphological degradation) and were ranked according to their relative importance as group and individually within and across groups using Principal Component Analyses (PCAs) and Boosted Regression Trees (BRTs). Overall, natural factors contributed the most to the total explained deviance of the models. This variable group represented not only typological differences between sampling sites but also a gradient of human impact by strongly anthropogenically influenced variables such as electric conductivity and dissolved oxygen. These large-scale effects can mask the individual importance of the other variable groups, which may act more specifically at a subset of sites. Accordingly, micropollutants were not represented by a few dominant variables but rather a diverse palette of different chemicals with similar contribution. As a group, micropollutants contributed similarly as metals, nutrient enrichment and morphological degradation. However, the importance of micropollutants might be underestimated due to limitations of the current chemical monitoring practices.

Effects of environment and metacommunity delineation on multiple dimensions of stream fish beta diversity

Introduction

Beta diversity represents changes in community composition among locations across a landscape. While the effects of human activities on beta diversity are becoming clearer, few studies have considered human effects on the three dimensions of beta diversity: taxonomic, functional, and phylogenetic. Including anthropogenic factors and multiple dimensions of biodiversity may explain additional variation in stream fish beta diversity, providing new insight into how metacommunities are structured within different spatial delineations.

Methods

In this study, we used a 350 site stream fish abundance dataset from South Carolina, United States to quantify beta diversity explainable by spatial, natural environmental, and anthropogenic variables. We investigated three spatial delineations: (1) a single whole-state metacommunity delineated by political boundaries, (2) two metacommunities delineated by a natural geomorphic break separating uplands from lowlands, and (3) four metacommunities delineated by natural watershed boundaries. Within each metacommunity we calculated taxonomic, functional, and phylogenetic beta diversity and used variation partitioning to quantify spatial, natural environmental, and anthropogenic contributions to variations in beta diversity.

Results

We explained 25–81% of the variation in stream fish beta diversity. The importance of these three factors in structuring metacommunities differed among the diversity dimensions, providing complementary perspectives on the processes shaping beta diversity in fish communities. The effect of spatial, natural environmental, and anthropogenic factors varied among the spatial delineations, which indicate conclusions drawn from variation partitioning may depend on the spatial delineation chosen by researchers.

Discussion

Our study highlights the importance of considering human effects on metacommunity structure, quantifying multiple dimensions of beta diversity, and careful consideration of user-defined metacommunity boundaries in beta diversity analyses.

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.

Identifying Key Stressors Driving Biological Impairment in Freshwater Streams in the Chesapeake Bay Watershed, USA

Biological communities in freshwater streams are often impaired by multiple stressors (e.g., flow or water quality) originating from anthropogenic activities such as urbanization, agriculture, or energy extraction. Restoration efforts in the Chesapeake Bay watershed, USA seek to improve biological conditions in 10% of freshwater tributaries and to protect the biological integrity of existing healthy watersheds. To achieve these goals, resource managers need to better understand which stressors are most likely driving biological impairment. Our study addressed this knowledge gap through two approaches: 1) reviewing and synthesizing published multi-stressor studies, and 2) examining 303(d) listed impairments linked to biological impairment as identified by jurisdiction regulatory agencies (the states within the watershed and the District of Columbia). Results identified geomorphology (i.e., physical habitat), salinity, and toxic contaminants as important for explaining variability in benthic community metrics in the literature review. Geomorphology (i.e., physical habitat and sediment), salinity, and nutrients were the most reported stressors in the jurisdictional impairment analysis. Salinity is likely a major stressor in urban and mining settings, whereas geomorphology was commonly reported in agricultural settings. Toxic contaminants, such as pesticides, were rarely measured; more research is needed to quantify the extent of their effects in the region. Flow alteration was also highlighted as an important urban stressor in the literature review but was rarely measured in the literature or reported by jurisdictions as a cause of impairment. These results can be used to prioritize stressor monitoring by managers, and to improve stressor identification methods for identifying causes of biological impairment.

Characterizing temporal variability in streams supports nutrient indicator development using diatom and bacterial DNA metabarcoding

Interest in developing periphytic diatom and bacterial indicators of nutrient effects continues to grow in support of the assessment and management of stream ecosystems and their watersheds. However, temporal variability could confound relationships between indicators and nutrients, subsequently affecting assessment outcomes. To document how temporal variability affects measures of diatom and bacterial assemblages obtained from DNA metabarcoding, we conducted weekly periphyton and nutrient sampling from July to October 2016 in 25 streams in a 1293 km² mixed land use watershed. Measures of both diatom and bacterial assemblages were strongly associated with the percent agriculture in upstream watersheds and total phosphorus (TP) and total nitrogen (TN) concentrations. Temporal variability in TP and TN concentrations increased with greater amounts of agriculture in watersheds, but overall diatom and bacterial assemblage variability within sites-measured as mean distance among samples to corresponding site centroids in ordination space-remained consistent. This consistency was due in part to offsets between decreasing variability in relative abundances of taxa typical of low nutrient conditions and increasing variability in those typical of high nutrient conditions as mean concentrations of TP and TN increased within sites. Weekly low and high nutrient diatom and bacterial metrics were more strongly correlated with site mean nutrient concentrations over the sampling period than with same day measurements and more strongly correlated with TP than with TN. Correlations with TP concentrations were consistently strong throughout the study except briefly following two major precipitation events. Following these events, biotic relationships with TP reestablished within one to three weeks. Collectively, these results can strengthen interpretations of survey results and inform monitoring strategies and decision making. These findings have direct applications for improving the use of diatoms and bacteria, and the use of DNA metabarcoding, in monitoring programs and stream site assessments.