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.

Bellamya limnophila is a mollusk of significant medical and economic value in China. Understanding the complete mitochondrial genome of this species will better establish a foundation for systematic classification research on Viviparidae. Therefore, we sequenced the complete mitochondrial genome of B. limnophila, conducted a comprehensive analysis of its structural characteristics, and constructed a phylogenetic tree using maximum likelihood and Bayesian methods. The results showed that the genome sequence is 16,991 bp in length, including 13 protein-coding genes (PCGs), 20 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and 1 non-coding region (D-loop). In summary, the Ka/Ks ratios of all PCGs were <1, indicating that purifying selection dominated the evolutionary process of these snails. The entire genome structure exhibited conservative features, such as the majority of start codons being the standard ATG codon and the majority of tRNA genes having the standard cloverleaf secondary structure. B. limnophila and B. quadrata showed collinearity in terms of sequence homology. Phylogenetic analysis indicates that the clade formed by the genera Margarya, Cipangopaludina, and Bellamya is the sister group of the genus Viviparus; Bellamya limnophila is more closely related to B. quadrata than to other species. This study contributes to the mitochondrial genome database of the family Viviparidae and provides valuable insights into the phylogenetic relationships of related snails.