Spatio-temporal patterns of multi-trophic biodiversity and food-web characteristics uncovered across a river catchment using environmental DNA

Ecological impacts of treated effluent on multitrophic biodiversity and their interactions

The reuse of water, particularly treated effluent from wastewater treatment plants (WWTPs), is a crucial and sustainable strategy for mitigating water scarcity, especially in megacities with high water demand and limited resources. However, the ecological risks associated with effluent discharge into receiving waterbodies have gained significant global attention. Understanding the dynamic effects of WWTP effluent on multi-trophic groups and their interactions is essential for assessing ecological impacts in aquatic ecosystems and informing management strategies. In this study, we examined five taxonomic groups representing different trophic levels of the freshwater food web - bacteria (decomposers), algae (primary producers), zooplankton (primary consumers), and benthic macroinvertebrates and fish (predators) - across two rivers to elucidate ecological responses to WWTP effluent from a multi-trophic perspective. Our results revealed significant but variable biological responses among these groups, depending on river conditions and trophic level. In the nutrient-rich river, primary consumers (zooplankton) were most affected, whereas in the nutrient-poor river, primary producers (algae) exhibited the strongest responses primarily derived from environmental disturbances. Notably, interactions between environmental variables and taxa were highly diverse, with trophic dynamics influenced by both bottom-up and top-down processes in the nutrient-rich river, whereas bottom-up effects dominated in the nutrient-poor river. Furthermore, niche overlap in algae-zooplankton networks was higher in the nutrient-rich river than in the nutrient-poor river. This study underscores the importance of integrating multi-trophic biodiversity profiling and trophic interaction analyses to comprehensively assess the ecological effects of WWTP effluent in receiving aquatic ecosystems with contrasting environmental contexts. Our findings highlight the importance of conservation and sustainable management practices, especially in urban aquatic ecosystems located in (semi-)arid regions that experience prolonged periods of low precipitation.

Riverine network size determined major driving factors of the composition and diversity of aquatic invertebrate communities in a multi-tributary mountain river basin

Revealing the spatial variation of aquatic invertebrates and their response to biotic and abiotic factors, from headwaters to estuaries, is crucial for understanding their successional patterns and protecting watershed ecosystems. This study aimed to explore the biogeographic patterns and identify the primary drivers of invertebrate community structure across river networks of varying sizes using environmental DNA (eDNA) technology. To assess the contribution of biotic and abiotic factors to invertebrate communities, we collected six categories of abiotic factors: geography, climate, hydro-morphology, human footprint index, land use, and water quality. For biotic factors, four microbial groups including archaea, bacteria, fungi, and protists were identified using eDNA techniques. Water samples were collected from a total of 187 sample sites in the upper Hanjiang River basin (China) during two seasons (Spring and Autumn), covering the transition from the headwater tributaries to the lower reaches of the main channel. The results revealed that environmental factors explained approximately 6.5 times more variation in invertebrate eDNA communities than geographic factors. Water quality and biotic factors had strong explanatory power for invertebrate eDNA diversity. Ecological succession of invertebrate eDNA communities along the river continuum showed a shift from Arthropoda-dominated communities in the headwaters to a co-dominance of Arthropoda, Rotifera, and Cnidaria downstream. The cumulative dendritic distance upstream, representing the location of each sampling site within the river network, emerged as the most predictive spatial feature. Significant differences were observed in the dominant environmental factors influencing community diversity across different river network sizes. In small river networks, invertebrate eDNA diversity was primarily influenced by biotic factors, while in medium-sized networks, it was shaped by a combination of biotic factors and water quality. In large river networks, water quality emerged as the primary driver. These findings suggest that invertebrate communities throughout the Hanjiang River basin undergo ecological succession along the river continuum, primarily shaped by environmental factors related to river network size.

Ecological restoration enhanced the stability of epiphytic microbial food webs of submerged macrophytes: Insights from predation characteristics of epiphytic predators

The application of various submerged macrophytes for ecological restoration has gained increasing attention in urban lake ecosystems. The multitrophic microbial communities that colonized in various submerged macrophytes constitute microbial food webs through trophic cascade effects, which affect the biogeochemical cycles of the lake ecosystem and directly determine the effects of ecological restoration. Therefore, it is essential to reveal the diversity, composition, assembly processes, and stability of the microbial communities within epiphytic food webs of diverse submerged macrophytes under eutrophication and ecological restoration scenarios. In this study, we explored the epiphytic microbial food webs of Vallisneria natans and Hydrilla verticillata in both eutrophic and ecological restoration regions. The obtained results indicated that the two regions with different nutrient levels remarkably affected the diversity and composition of epiphytic multitrophic microbial communities of submerged macrophytes, among them, the community composition of epiphytic predators were more prone to change. Secondly, environmental filtering effects played a more important role in driving the community assembly of epiphytic predators than that of prey. Furthermore, the generality and intraguild predation of epiphytic predators were significantly improved within ecological restoration regions, which increased the stability of epiphytic microbial food webs. Additionally, compared with Hydrilla verticillata, the epiphytic microbial food webs of Vallisneria natans exhibited higher multitrophic diversity and higher network stability regardless of regions. Overall, this study focused on the role of the epiphytic microbial food webs of submerged macrophytes in ecological restoration and uncovered the potential of epiphytic predators to enhance the stability of microbial food webs, which may provide new insights into the development of ecological restoration strategies.

Evaluating eDNA and eRNA metabarcoding for aquatic biodiversity assessment: From bacteria to vertebrates

The monitoring and management of aquatic ecosystems depend on precise estimates of biodiversity. Metabarcoding analyses of environmental nucleic acids (eNAs), including environmental DNA (eDNA) and environmental RNA (eRNA), have garnered attention for their cost-effective and non-invasive biomonitoring capabilities. However, the accuracy of biodiversity estimates obtained through eNAs can vary among different organismal groups. Here we evaluate the performance of eDNA and eRNA metabarcoding across nine organismal groups, ranging from bacteria to terrestrial vertebrates, in three cross-sections of the Yangtze River, China. We observe robust complementarity between eDNA and eRNA data. The relative detectability of eNAs was notably influenced by major taxonomic groups and organismal sizes, with eDNA providing more robust signals for larger organisms. Both eDNA and eRNA exhibited similar cross-sectional and longitudinal patterns. However, the detectability of larger organisms declined in eRNA metabarcoding, possibly due to differential RNA release and decay among different organismal groups or sizes. While underscoring the potential of eDNA and eRNA in large river biomonitoring, we emphasize the need for differential interpretation of eDNA versus eRNA data. This highlights the importance of careful method selection and interpretation in biomonitoring studies.

Measuring the state of aquatic environments using eDNA-upscaling spatial resolution of biotic indices

Aquatic macroinvertebrates, including many aquatic insect orders, are a diverse and ecologically relevant organismal group yet they are strongly affected by anthropogenic activities. As many of these taxa are highly sensitive to environmental change, they offer a particularly good early warning system for human-induced change, thus leading to their intense monitoring. In aquatic ecosystems there is a plethora of biotic monitoring or biomonitoring approaches, with more than 300 assessment methods reported for freshwater taxa alone. Ultimately, monitoring of aquatic macroinvertebrates is used to calculate ecological indices describing the state of aquatic systems. Many of the methods and indices used are not only hard to compare, but especially difficult to scale in time and space. Novel DNA-based approaches to measure the state and change of aquatic environments now offer unprecedented opportunities, also for possible integration towards commonly applicable indices. Here, we first give a perspective on DNA-based approaches in the monitoring of aquatic organisms, with a focus on aquatic insects, and how to move beyond traditional point-based biotic indices. Second, we demonstrate a proof-of-concept for spatially upscaling ecological indices based on environmental DNA, demonstrating how integration of these novel molecular approaches with hydrological models allows an accurate evaluation at the catchment scale. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Distinct multitrophic biodiversity composition and community organization in a freshwater lake and a hypersaline lake on the Tibetan Plateau

Alpine lakes play pivotal roles in plateau hydrological processes but are highly sensitive to climate change, yet we lack comprehensive knowledge of their multitrophic biodiversity patterns. Here, we compared the biodiversity characteristics of diverse taxonomic groups across water depths and in surface sediments from a freshwater lake and a hypersaline lake on the northwestern Tibetan Plateau. Using multi-marker environmental DNA metabarcoding, we detected 134 cyanobacteria, 443 diatom, 1,519 invertebrate, and 28 vertebrate taxa. Each group had a substantially different community composition in the two lakes, and differences were also found between water and sediments within each lake. Cooccurrence network analysis revealed higher network complexity, lower modularity, and fewer negative cohesions in the hypersaline lake, suggesting that high salinity may destabilize ecological networks. Our results provide the first holistic view of Tibetan lake biodiversity under contrasting salinity levels and reveal structural differences in the ecological networks that may impact ecosystem resilience.

Environmental DNA: The next chapter

Molecular tools are an indispensable part of ecology and biodiversity sciences and implemented across all biomes. About a decade ago, the use and implementation of environmental DNA (eDNA) to detect biodiversity signals extracted from environmental samples opened new avenues of research. Initial eDNA research focused on understanding population dynamics of target species. Its scope thereafter broadened, uncovering previously unrecorded biodiversity via metabarcoding in both well-studied and understudied ecosystems across all taxonomic groups. The application of eDNA rapidly became an established part of biodiversity research, and a research field by its own. Here, we revisit key expectations made in a land-mark special issue on eDNA in Molecular Ecology in 2012 to frame the development in six key areas: (1) sample collection, (2) primer development, (3) biomonitoring, (4) quantification, (5) behaviour of DNA in the environment and (6) reference database development. We pinpoint the success of eDNA, yet also discuss shortfalls and expectations not met, highlighting areas of research priority and identify the unexpected developments. In parallel, our retrospective couples a screening of the peer-reviewed literature with a survey of eDNA users including academics, end-users and commercial providers, in which we address the priority areas to focus research efforts to advance the field of eDNA. With the rapid and ever-increasing pace of new technical advances, the future of eDNA looks bright, yet successful applications and best practices must become more interdisciplinary to reach its full potential. Our retrospect gives the tools and expectations towards concretely moving the field forward.

An integrated spatio-temporal view of riverine biodiversity using environmental DNA metabarcoding

Anthropogenically forced changes in global freshwater biodiversity demand more efficient monitoring approaches. Consequently, environmental DNA (eDNA) analysis is enabling ecosystem-scale biodiversity assessment, yet the appropriate spatio-temporal resolution of robust biodiversity assessment remains ambiguous. Here, using intensive, spatio-temporal eDNA sampling across space (five rivers in Europe and North America, with an upper range of 20-35 km between samples), time (19 timepoints between 2017 and 2018) and environmental conditions (river flow, pH, conductivity, temperature and rainfall), we characterise the resolution at which information on diversity across the animal kingdom can be gathered from rivers using eDNA. In space, beta diversity was mainly dictated by turnover, on a scale of tens of kilometres, highlighting that diversity measures are not confounded by eDNA from upstream. Fish communities showed nested assemblages along some rivers, coinciding with habitat use. Across time, seasonal life history events, including salmon and eel migration, were detected. Finally, effects of environmental conditions were taxon-specific, reflecting habitat filtering of communities rather than effects on DNA molecules. We conclude that riverine eDNA metabarcoding can measure biodiversity at spatio-temporal scales relevant to species and community ecology, demonstrating its utility in delivering insights into river community ecology during a time of environmental change.

Towards (better) fluvial meta-ecosystem ecology: a research perspective

Rivers are an important component of the global carbon cycle and contribute to atmospheric carbon exchange disproportionately to their total surface area. Largely, this is because rivers efficiently mobilize, transport and metabolize terrigenous organic matter (OM). Notably, our knowledge about the magnitude of globally relevant carbon fluxes strongly contrasts with our lack of understanding of the underlying processes that transform OM. Ultimately, OM processing en route to the oceans results from a diverse assemblage of consumers interacting with an equally diverse pool of resources in a spatially complex network of heterogeneous riverine habitats. To understand this interaction between consumers and OM, we must therefore account for spatial configuration, connectivity, and landscape context at scales ranging from local ecosystems to entire networks. Building such a spatially explicit framework of fluvial OM processing across scales may also help us to better predict poorly understood anthropogenic impacts on fluvial carbon cycling, for instance human-induced fragmentation and changes to flow regimes, including intermittence. Moreover, this framework must also account for the current unprecedented human-driven loss of biodiversity. This loss is at least partly due to mechanisms operating across spatial scales, such as interference with migration and habitat homogenization, and comes with largely unknown functional consequences. We advocate here for a comprehensive framework for fluvial networks connecting two spatially aware but disparate lines of research on (i) riverine metacommunities and biodiversity, and (ii) the biogeochemistry of rivers and their contribution to the global carbon cycle. We argue for a research agenda focusing on the regional scale-that is, of the entire river network-to enable a deeper mechanistic understanding of naturally arising biodiversity-ecosystem functioning coupling as a major driver of biogeochemically relevant riverine carbon fluxes.

Environmental DNA and remote sensing datasets reveal the spatial distribution of aquatic insects in a disturbed subtropical river system

Biodiversity datasets with high spatial resolution are critical prerequisites for river protection and management decision-making. However, traditional morphological biomonitoring is inefficient and only provides several site estimates, and there is an urgent need for new approaches to predict biodiversity on fine spatial scales throughout the entire river systems. Here, we combined the environmental DNA (eDNA) and remote sensing (RS) technologies to develop a novel approach for predicting the spatial distribution of aquatic insects with high spatial resolution in a disturbed subtropical Dongjiang River system of southeast China. First, we screened thirteen RS-based vegetation indices that significantly correlated with the eDNA-inferred richness of aquatic insects. In particular, the green normalized difference vegetation index (GNDVI) and normalized difference red-edge2 (NDRE2) were closely related to eDNA-inferred richness. Second, using the gradient boosting decision tree, our data showed that the spatial pattern of eDNA-inferred richness could achieve a high spatial resolution to 500 m reach and accurate prediction of more than 80%, and the prediction efficiency of the headwater streams (Strahler stream order = 1) was slightly higher than the downstream (Strahler stream order >1). Third, using the random forest algorithm, the spatial distribution of aquatic insects could reach a prediction rate of over 70% for the presence or absence of specific genera. Overall, this study provides a new approach to achieving high spatial resolution prediction of the distribution of aquatic insects, which supports decision-making on river diversity protection under climate changes and human impacts.

Hypoxia and salinity constrain the sediment microbiota-mediated N removal potential in an estuary: A multi-trophic interrelationship perspective

Reactive nitrogen (N) enrichment is a common environmental problem in estuarine ecosystems, while the microbial-mediated N removal process is complicated for other multi-environmental factors. Therefore, A systematic investigation is necessary to understand the multi-trophic microbiota-mediated N removal characteristics under various environmental factors in estuaries. Here, we studied how multiple factors affect the multi-trophic microbiota-mediated N removal potential (denitrification and anammox) and N2O emission along a river-estuary-bay continuum in southeastern China using the environmental DNA (eDNA) approach. Results suggested that hypoxia and salinity were the dominant environmental factors affecting multi-trophic microbiota-mediated N removal in the estuary. The synergistic effect of hypoxia and salinity contributed to the loss of taxonomic (MultiTaxa) and phylogenetic (MultiPhyl) diversity across multi-trophic microbiota and enhanced the interdependence among multi-trophic microbiota in the estuary. The N removal potential calculated as the activities of key N removal enzymes was also significantly constrained in the estuary (0.011), compared with the river (0.257) and bay (0.461). Structural equation modeling illustrated that metazoans were central to all sediment N removal potential regulatory pathways. The top-down forces (predation by metazoans) restrained the growth of heterotrophic bacteria, which may affect microbial N removal processes in the sediment. Furthermore, we found that the hypoxia and salinity exacerbated the N2O emission in the estuary. This study clarifies that hypoxia and salinity constrain estuarine multi-trophic microbiota-mediated N removal potential and highlights the important role of multi-trophic interactions in estuarine N removal, providing a new perspective on mitigating estuarine N accumulation.

eDNA-based diversity and multitrophic network reveal the effects of land use and pollutants on the subtropical Dongjiang River systems

Biodiversity and its constituted multitrophic network in rivers are accelerating change under human land use and pollutants. However, due to the lack of complete datasets across taxa limited by traditional morphological biomonitoring, the change patterns of biodiversity and multitrophic networks are still unclear. Here, we used the eDNA approach to capture multitrophic communities (including fish, aquatic insects, protozoa, diatom and bacteria) in the Dongjiang River, a typical subtropical river in southeast China, and analyzed the changing patterns of biodiversity and multitrophic networks in relation to land use and water pollution. First, our data showed that the eDNA approach provided a snapshot of the multitrophic communities in the Dongjiang River, and the monitored 5833 OTUs were annotated to 55 phyla, 144 classes, 329 orders, 521 families, 945 genera and 406 species. Second, the multitrophic diversity index had similar patterns on the longitudinal scale of rivers, with significant decreases from the upstream to the downstream, while individual taxonomic groups exhibited variable spatial patterns. While there were similar spatial patterns between network metrics and diversity index, the former had stronger relationships with the spatial distance. Third, the multitrophic diversity and networks were significantly negatively correlated with land use and water pollution (e.g., CODMn), and network structures often had stronger and non-linear responses. Overall, this study highlights that eDNA biomonitoring of multitrophic communities and networks can provide deeper insights into ecosystem changes and help develop more targeted management strategies.

metanetwork: A R package dedicated to handling and representing trophic metanetworks

Trophic networks describe interactions between species at a given location and time. Due to environmental changes, anthropogenic perturbations or sampling effects, trophic networks may vary in space and time. The collection of network time series or networks in different sites thus constitutes a metanetwork. We present here the R package metanetwork, which will ease the representation, the exploration and analysis of trophic metanetwork data sets that are increasingly available. Our main methodological advance consists in suitable layout algorithm for trophic networks, which is based on trophic levels and dimension reduction in a graph diffusion kernel. In particular, it highlights relevant features of trophic networks (trophic levels, energetic channels). In addition, we developed tools to handle, compare visually and quantitatively and aggregate those networks. Static and dynamic visualisation functions have been developed to represent large networks. We apply our package workflow to several trophic network data sets.

Modelling environmental DNA transport in rivers reveals highly resolved spatio-temporal biodiversity patterns

The ever-increasing threats to riverine ecosystems call for novel approaches for highly resolved biodiversity assessments across taxonomic groups and spatio-temporal scales. Recent advances in the joint use of environmental DNA (eDNA) data and eDNA transport models in rivers (e.g., eDITH) allow uncovering the full structure of riverine biodiversity, hence elucidating ecosystem processes and supporting conservation measures. We applied eDITH to a metabarcoding dataset covering three taxonomic groups (fish, invertebrates, bacteria) and three seasons for a catchment sampled for eDNA at 73 sites. We upscaled eDNA-based biodiversity predictions to approximately 1900 reaches, and assessed α- and β-diversity patterns across seasons and taxonomic groups over the whole network. Genus richness predicted by eDITH was generally higher than values from direct eDNA analysis. Both predicted α- and β-diversity varied depending on season and taxonomic group. Predicted fish α-diversity increased downstream in all seasons, while invertebrate and bacteria α-diversity either decreased downstream or were unrelated to network position. Spatial β-diversity mostly decreased downstream, especially for bacteria. The eDITH model yielded a more refined assessment of freshwater biodiversity as compared to raw eDNA data, both in terms of spatial coverage, diversity patterns and effect of covariates, thus providing a more complete picture of freshwater biodiversity.

Destabilizing Effects of Environmental Stressors on Aquatic Communities and Interaction Networks across a Major River Basin

Human-driven environmental stressors are increasingly threatening species survival and diversity of river systems worldwide. However, it remains unclear how the stressors affect the stability changes across aquatic multiple communities. Here, we used environmental DNA (eDNA) data sets from a human-dominated river in China over 3 years and analyzed the stability changes in multiple communities under persistent anthropogenic stressors, including land use and pollutants. First, we found that persistent stressors significantly reduced multifaceted species diversity (e.g., species richness, Shannon's diversity, and Simpson's diversity) and species stability but increased species synchrony across multiple communities. Second, the structures of interaction networks inferred from an empirical meta-food web were significantly changed under persistent stressors, for example, resulting in decreased network modularity and negative/positive cohesion. Third, piecewise structural equation modeling proved that the persistent stress-induced decline in the stability of multiple communities mainly depended upon diversity-mediated pathways rather than the direct effects of stress per se; specifically, the increase of species synchrony and the decline of interaction network modularity were the main biotic drivers of stability variation. Overall, our study highlights the destabilizing effects of persistent stressors on multiple communities as well as the mechanistic dependencies, mainly through reducing species diversity, increasing species synchrony, and changing interaction networks.

Multi-group biodiversity distributions and drivers of metacommunity organization along a glacial-fluvial-limnic pathway on the Tibetan plateau

Extensive global glacial retreats are threatening cryosphere ecosystem functioning and the associated biota in glacier-fed water systems. Understanding multi-group biodiversity distributions and compositional variation across diverse but hydrologically linked habitats under varying glacial influences will help explain the mechanisms underlying glacial community organization and ecosystem processes. However, such data are generally lacking due to the difficulty of obtaining biodiversity information across wide taxonomic ranges. Here, we used a multi-marker environmental DNA metabarcoding approach to simultaneously investigate the spatial patterns of community compositions and assembly mechanisms of four taxonomic groups (cyanobacteria, diatoms, invertebrates, and vertebrates) along the flowpaths of a tributary of Lake Nam Co on the Tibetan Plateau-from its glacier headwaters, through its downstream river and wetlands, to its estuary. We detected 869 operational taxonomic units: 119 cyanobacterial, 395 diatom, 269 invertebrate, and 86 vertebrate. Taxonomic richnesses consistently increased from upstream to downstream, and although all groups showed community similarity distance decay patterns, the trend for vertebrates was the weakest. Cyanobacteria, diatom, and invertebrate community compositions were significantly correlated with several environmental factors, while the vertebrate community was only correlated with waterway width. Variation partitioning analysis indicated that varying extents of environmental conditions and spatial factors affected community organizations for different groups. Furthermore, stochastic processes contributed prominently to the microorganisms' community assembly (Sloan's neutral model R² = 0.77 for cyanobacteria and 0.73 for diatoms) but were less important for macroorganisms (R² = 0.21 for invertebrates and 0.15 for vertebrates). That trend was further substantiated by modified stochasticity ratio analyses. This study provides the first holistic picture of the diverse biotic communities residing in a series of hydrologically connected glacier-influenced habitats. Our results both uncovered the distinct mechanisms that underlie the metacommunity organizations of different glacial organisms and helped comprehensively predict the ecological impacts of the world's melting glaciers.

Changes of In Situ Prokaryotic and Eukaryotic Communities in the Upper Sanya River to the Sea over a Nine-Hour Period

The transition areas of riverine, estuarine, and marine environments are particularly valuable for the research of microbial ecology, biogeochemical processes, and other physical-chemical studies. Although a large number of microbial-related studies have been conducted within such systems, the vast majority of sampling have been conducted over a large span of time and distance, which may lead to separate batches of samples receiving interference from different factors, thus increasing or decreasing the variability between samples to some extent. In this study, a new in situ filtration system was used to collect membrane samples from six different sampling sites along the Sanya River, from upstream freshwater to the sea, over a nine-hour period. We used high-throughput sequencing of 16S and 18S rRNA genes to analyze the diversity and composition of prokaryotic and eukaryotic communities. The results showed that the structures of these communities varied according to the different sampling sites. The α-diversity of the prokaryotic and eukaryotic communities both decreased gradually along the downstream course. The structural composition of prokaryotic and eukaryotic communities changed continuously with the direction of river flow; for example, the relative abundances of Rhodobacteraceae and Flavobacteriaceae increased with distance downstream, while Sporichthyaceae and Comamonadaceae decreased. Some prokaryotic taxa, such as Phycisphaeraceae and Chromobacteriaceae, were present nearly exclusively in pure freshwater environments, while some additional prokaryotic taxa, including the SAR86 clade, Clade I, AEGEAN-169 marine group, and Actinomarinaceae, were barely present in pure freshwater environments. The eukaryotic communities were mainly composed of the Chlorellales X, Chlamydomonadales X, Sphaeropleales X, Trebouxiophyceae XX, Annelida XX, and Heteroconchia. The prokaryotic and eukaryotic communities were split into abundant, common, and rare communities for NCM analysis, respectively, and the results showed that assembly of the rare community assembly was more impacted by stochastic processes and less restricted by species dispersal than that of abundant and common microbial communities for both prokaryotes and eukaryotes. Overall, this study provides a valuable reference and new perspectives on microbial ecology during the transition from freshwater rivers to estuaries and the sea.

Environmental DNA Biomonitoring Reveals the Interactive Effects of Dams and Nutrient Enrichment on Aquatic Multitrophic Communities

Dam construction and nutrient enrichment are two pervasive stressors in rivers worldwide, which trigger a sharp decline in biodiversity and ecosystem services. However, the interactive effects of both stressors on multitrophic taxonomic groups remain largely unclear. Here, we used the multitrophic datasets captured by the environmental DNA (eDNA) approach to reveal the interactions between dams and nutrient enrichment on aquatic communities from the aspects of taxonomic α diversity, β diversity, and food webs. First, our data showed that dams and nutrient enrichment jointly shaped a unique spatial pattern of aquatic communities across the four river systems, and the dissimilarity of community structure significantly declined (i.e., structural homogenization) under both stressors. Second, dams and nutrients together explained 40-50% of the variations in aquatic communities, and dams had a stronger impact on fish, aquatic insects, and bacteria, yet nutrients had a stronger power to drive protozoa, fungi, and eukaryotic algae. Finally, we found that additive, synergistic, and antagonistic interactions of dams and nutrient enrichment were common and coexisted in river systems and led to significantly simplified aquatic food webs, with decreases in modularity (synergistic) and robustness (additive) and an increase in coherence (synergistic). Overall, our study highlights that eDNA-based datasets can provide multitrophic perspectives for fostering the understanding of the interactive effects of multiple stressors on rivers.