Unravelling Metazoan and Fish Community Patterns in Yujiang River, China: Insights from Beta Diversity Partitioning and Co-Occurrence Network

Understanding the biodiversity of aquatic communities and the underlying mechanisms that shape biodiversity patterns and community dynamics is crucial for the effective conservation and management of freshwater ecosystems. However, traditional survey methods often fail to comprehensively capture species diversity, particularly for low-abundance taxa. Moreover, studies integrating both metazoan and fish communities at fine spatial scales remain limited. To address these gaps, we employed a multi-marker eDNA metabarcoding approach, targeting both the 12S and 18S rRNA gene regions, to comprehensively investigate the composition of metazoan and fish communities in the Yujiang River. A total of 12 metazoan orders were detected, encompassing 15 families, 21 genera, and 19 species. For the fish community, 32 species were identified, belonging to 25 genera, 10 families, and 7 orders. Among these, Adula falcatoides and Coptodon zillii were identified as the most prevalent and abundant metazoan and fish species, respectively. Notably, the most prevalent fish species, C. zillii and Oreochromis niloticus, are both recognized as invasive species. The Bray–Curtis distance of metazoa (average: 0.464) was significantly lower than that of fish communities (average: 0.797), suggesting higher community heterogeneity among fish assemblages. Beta-diversity decomposition indicated that variations in the metazoan and fish communities were predominantly driven by species replacement (turnover) (65.4% and 70.9% for metazoa and fish, respectively) rather than nestedness. Mantel tests further revealed that species turnover in metazoan communities was most strongly influenced by water temperature, while fish community turnover was primarily affected by water transparency, likely reflecting the physiological sensitivity of metazoans to thermal gradients and the dependence of fish on visual cues for foraging and habitat selection. In addition, a co-occurrence network of metazoan and fish species was constructed, highlighting potential predator-prey interactions between native species and Corbicula fluminea, which emerged as a potential keystone species. Overall, this study demonstrates the utility of multi-marker eDNA metabarcoding in characterizing aquatic community structures and provides new insights into the spatial dynamics and species interactions within river ecosystems.