Suxing Fu

Mountain river networks are composed of a mainstem and multiple tributaries, where complex hydrological processes and habitat heterogeneity provide diverse niches for fishes. Previous studies have documented the fish fauna and diversity in the Dadu River basin, providing baseline information for regional fish conservation. However, differences in fish assemblage composition among river units and their relationships with water environmental conditions require further examination through continuous multi-year surveys. From 2023 to 2025, a three-year continuous survey of fish assemblages was conducted across 28 sampling sites in the mainstem and nine major tributaries of the Shimian reach in the middle Dadu River (upper Yangtze River basin), along with concurrent monitoring of water physicochemical parameters. A total of 54 fish species belonging to three orders and seven families were recorded, with Cyprinidae and Cobitidae being the dominant groups. Fish assemblages were mainly composed of rheophilic, benthic, and demersal-spawning species, reflecting the typical ecological characteristics of fishes in high-gradient river habitats. Fish assemblage composition remained generally stable over the three survey years but showed marked spatial differences among river units. NMDS ordination and PERMANOVA indicated that functional habitat types (mainstem, upstream and confluences) explained a significant proportion of variation in assemblage structure, whereas the effect of survey year was not significant, indicating that spatial differences across habitat types were more pronounced than interannual variation. Distinct fish assemblages were formed among the mainstem, confluences, and individual tributaries, with several natural tributaries supporting higher proportions of endemic and cold-water fishes. Water environmental analyses revealed evident hydrochemical gradients among river units, and redundancy analysis (RDA) showed that water temperature, dissolved oxygen, and nutrient concentrations were closely associated with the spatial distribution of fish assemblages. Building on previous surveys of fish fauna in the Dadu River basin, this study adds three years of fish assemblage data at the mainstem–tributary river network scale and further describes spatial differences in assemblage composition among river units and their relationships with water environmental conditions. These results provide baseline information for understanding the spatial distribution of fish assemblages in alpine canyon rivers and for fish biodiversity conservation and tributary habitat management in the Dadu River basin.

Golden coloration can arise from changes in pigment production, pigment cell composition, or the function of pigment-containing organelles. Here, we characterized a golden variant (GV) of mandarin fish (Siniperca chuatsi) using morphological examination, skin pigment quantification, RNA sequencing, and untargeted metabolomics. Compared with wild-type (WT) fish, GV fish exhibited reduced visible melanization, lower melanin and pteridine contents, and higher β-carotene content. Transcriptomic analysis identified 1,163 differentially expressed genes, including 677 upregulated and 486 downregulated in GV relative to WT. Despite reduced melanin deposition, several genes associated with melanogenesis and melanosome biology were upregulated. Pathway enrichment analysis of upregulated genes and gene set enrichment analysis further supported changes in pigmentation-related transcriptional programs. Untargeted metabolomics yielded 2,995 putative metabolite annotations and identified 509 candidate differential features. However, no feature was significant at q < 0.05 after Benjamini-Hochberg correction across all tested annotations. A feature putatively annotated as L-tyrosine showed higher signal intensity in GV and mapped to pathways also represented in the transcriptomic dataset. Together, these findings show that reduced melanization and altered pigment composition in GV coexist with increased expression of several pigmentation-related genes. The candidate metabolic features require further validation, and the genetic basis and cellular processes underlying reduced melanin deposition remain unresolved.

Phytoplankton are sensitive indicators of environmental change in high-altitude aquatic ecosystems, but their community responses across river–lake confluence zones remain insufficiently understood. We investigated phytoplankton community characteristics and environmental associations in the Siling Co river–lake continuum. This endorheic system lies above 4500 m, where river inflows, lake-water mixing, and evaporative concentration generate marked hydrochemical gradients. Genus-level cell density and estimated biomass were assessed with 22 environmental variables from 14 sampling events at 10 lake and river sites. We recorded 54 phytoplankton genera, with Bacillariophyta dominant. Lake and River samples overlapped in composition, whereas sampling period accounted for more community variation than water-body type. Major-ion and salinity-related gradients were associated with cell-density community dissimilarity after accounting for water-body type. Cell-density and biomass composition also showed partly distinct associations with pH and SiO₂, and with Cl⁻, respectively. These observations indicate that phytoplankton communities in the Siling Co river–lake confluence zone are more closely associated with sampling-period variation and hydrochemical gradients than with Lake/River habitat classification. These findings provide an initial basis for integrated biological and hydrochemical monitoring.

River–lake confluence zones are shaped by inflow hydrology, evaporative concentration, and water mixing, providing natural systems for examining phytoplankton responses to environmental gradients. Here, we investigated Siling Co and adjacent river–lake confluence zones at elevations above 4,500 m. We integrated genus-level phytoplankton cell-density and biomass data with 22 environmental variables from 14 sampling events at 10 sites, including seven Lake and seven River samples. Community diversity, composition, and community–environment associations were evaluated using ordination, permutation tests, correlation analyses, distance-based redundancy analysis, Mantel tests, and multi-method evidence integration. In total, 54 genera were recorded, with Bacillariophyta representing the dominant phylum. Lake and River samples overlapped in community composition, whereas sampling period explained more genus-level variation than water-body type. Major-ion and mineralization gradients were consistently associated with community turnover. pH and SiO₂ were associated with cell-density composition, whereas Cl⁻ was associated with biomass composition. Cyanobacterial relative biomass was positively correlated with total nitrogen and negatively correlated with SO₄²⁻. These findings indicate that phytoplankton community structure in the Siling Co river–lake confluence zone cannot be attributed simply to two discrete habitat types, namely lakes and rivers, and provide a new perspective on community assembly in extremely high-altitude river–lake confluence zones.

The river ecosystems of the Qinghai–Tibet Plateau, recognized as a vital component of the “Asian Water Tower,” possess unique hydrological conditions and extreme environments that have shaped key indicator groups, most notably zooplankton. The community dynamics and structural characteristics of these zooplankton exhibit regular spatio-temporal distribution patterns across elevational gradients and seasonal successions. However, the intrinsic mechanisms underlying community succession and their correlations with environmental factors remain poorly understood, and the primary environmental drivers influencing community structure require further elucidation. Based on systematic zooplankton surveys and environmental data collection conducted across the Lhasa River basin from 2019 to 2021, this study established a comprehensive species inventory comprising 113 taxa across four major groups, alongside a multi-dimensional environmental dataset. We analyzed the spatio-temporal heterogeneities of zooplankton community structures—including abundance, biomass, and diversity indices—across different seasons and river reaches. The results revealed the composition and seasonal turnover of dominant taxa, with rotifers accounting for 39.82% of the total taxonomic richness. Mean zooplankton abundance and biomass across the basin were 1.18 ind./L and 343.60 × 10⁻⁵ mg/L, respectively, with peak values observed during autumn and within the Chabalang Wetland. The zooplankton community structure in the upstream, midstream, and downstream reaches, as well as associated wetlands, was significantly correlated with specific environmental factors (p < 0.05), including ammoniacal nitrogen (NH₄⁺-N), magnesium (Mg²⁺), total hardness (TH), potassium (K⁺), iron (Fe²⁺), sodium (Na⁺), sulfite (SO₃²⁻), nitrate ion (NO₃⁻), chloride ion (Cl⁻), total phosphorus (TP), and sulfide (S²⁻). Cl⁻, TH, Mg²⁺, SO₃²⁻, and elevation (Ele) were the key environmental drivers significantly influencing zooplankton abundance across seasons (p < 0.05). Furthermore, zooplankton abundance decreased significantly with increasing elevation during the winter. This research deepens our understanding of community assembly mechanisms in plateau river ecosystems and provides a scientific foundation for aquatic biodiversity conservation and ecological management in the Lhasa River basin.

Accurate sex identification is crucial for the artificial breeding of fish species, and identification based on phenotypic traits is the preferred method owing to its low cost and rapidity. This study aimed to investigate morphological differences between male and female Gymnocypris eckloni Herzenstein. A total of 166 G. eckloni individuals (83 males and 83 females) were included, and 20 morphological indicators, including head length, snout length, and body weight, were measured. Various analytical methods, such as correlation analysis, principal component analysis, discriminant analysis, and systematic cluster analysis, were used to compare the morphological characteristics of male and female G. eckloni. The result showed significant differences in morphological characteristics between male and female individuals (p ≤ 0.05). Principal component analysis revealed five main components that collectively contributed to 83.64% of the variance, primarily including body length, fork length, total length, and body weight, highlighting the significance of these components in sex differentiation. Subsequently, a mathematical model was constructed via discriminant analysis. This model demonstrated an accuracy of 79.5% and 73.5% in identifying male and female G. eckloni, respectively, with an overall accuracy of 76.50%. Furthermore, on analyzing the relationship between body length and body weight in G. eckloni, we found differences in growth patterns between sexes. In particular, female fish usually weighed more than male fish, which was consistent with the biological characteristics and energy distribution patterns of G. eckloni. These results suggest that sexual dimorphism in G. eckloni is primarily concentrated in body shape and head morphology, which provides a reliable basis for sex identification based on morphological differences.

The dynamic changes in plateau river ecosystems and the driving mechanisms of environmental factors have garnered significant attention. Phytoplankton, a core component of aquatic ecosystems, can directly reflect changes in the aquatic environment. This study focuses on the phytoplankton in the Lhasa River Basin, including the riverbed from the source to the river mouth, five largest tributaries, and two adjacent wetlands. We analyzed the spatial and temporal variation characteristics of phytoplankton and explored the environmental driving mechanisms based on four field surveys conducted between September 2019 and March 2021. Results showed that a total of 127 species of phytoplankton from six algal phyla were identified, including Cyanobacteria. Among these, Bacillariophyta was the dominant group, accounting for 41.7% of the identified species. Spatially, phytoplankton diversity showed a decreasing trend from upstream to downstream while temporally peaking in spring and autumn. Redundancy analysis revealed that upstream phytoplankton were driven by total hardness and altitude, midstream by pH and potassium ions, and downstream by nitrate and ammonium nitrogen. Classification and regression tree analysis showed total hardness, magnesium ions, and nitrite were key factors influencing phytoplankton abundance, diversity, and evenness. This study highlights the ecological dynamics and driving mechanisms of phytoplankton communities in the Lhasa River, demonstrating their high sensitivity to environmental factors. These findings could help to establish phytoplankton as critical indicators of aquatic ecosystem health and provide scientific guidance for the conservation and management of the plateau river ecosystems.

Schizothoracine fishes in saltwater lakes of the Tibetan Plateau are important models for studying the evolution and uplift of the Tibetan Plateau. Examining their adaptation to the high-salt environment is interesting. In this study, we first assembled the RNA-Seq data of each tissue of G. przewalskii, G. selincuoensis, and G. namensis from Qinghai Lake, Selincuo Lake, and Namtso Lake, respectively, obtained by the group previously. After obtaining reliable results, the adaptation of the gills, kidneys, and livers of the three species to the high-salinity environment was assessed by weighted gene co-expression network analysis (WGCNA). Using module eigengenes (ME), 21, 22, and 22 gene modules were identified for G. przewalskii, G. selincuoensis, and G. nemesis, respectively. Functional clustering analysis of genes in the significant association module identified several genes associated with osmolarity-regulated potential KEGG pathways in the gills of three species of Schizothoracine fish. Th17 cell differentiation pathway was up-regulated in the gills of all three species; histocompatibility class 2 II antigen and E alpha (h2-ea) were up-regulated genes in this pathway. Functional clustering analysis of genes in apparently related modules in the kidney unveiled several differential KEGG pathways. The pentose phosphate pathway was up-regulated in the three Schizothoracine fishes, and glucose-6-phosphate dehydrogenase (g6pd) was an up-regulated gene in this pathway. In the livers of the three Schizothorax species, the propanoate metabolism pathway was up-regulated, and succinate-CoA ligase GDP-forming subunit beta (suclg2) was an up-regulated gene in this pathway. The above analyses provide reference data for the adaptation of Schizothorax to high-salt environments and lay the foundation for future studies on the adaptive mechanism of Schizothorax in the plateau. These results partly fill the void in the knowledge gap in the survival adaptations of Schizothoracine fishes to highland saline lakes.

Two previously unknown species of Rhinogobius have been discovered in the streams of the Upper Youshui River, within the Yuan River Basin, Xiushan County, Chongqing, China. These new species are named as Rhinogobius sudoccidentalis and Rhinogobius lithopolychroma. Phylogenetic analysis based on mitochondrial genomes revealed that R. sudoccidentalis is genetically closest to R. reticulatus, while R. lithopolychroma shares the greatest genetic similarity with R. leavelli. Morphological distinctions allow for the clear differentiation of these species. Rhinogobius sudoccidentalissp. nov. is characterized by having VI–VII rays in the first dorsal fin and I, 8–9 rays in the second dorsal fin. The longitudinal scale series typically consists of 22–24 scales, while the transverse scale series comprises 7–8 scales. Notably, the predorsal scale series is absent and the total vertebrae count is 12+17=29. Rhinogobius lithopolychromasp. nov. can be distinguished from other species by the presence of 13–15 rays on the pectoral fin. Its longitudinal scale series ranges from 30 to 33 scales, with no scales in the predorsal area. The total vertebral count is 30, with 12 precaudal and 18 caudal vertebrae. The head and body of this species are light gray with irregular orange markings on the cheeks and opercle. Through morphological and molecular analyses, it has been confirmed that R. lithopolychroma and R. sudoccidentalis represent novel species within the Rhinogobius genus.

The Lhasa River is one of the five tributaries of the Yarlung Zangbo River. It is the main industrial, agricultural and animal husbandry area in Tibet. It plays an important ecological security barrier role in regulating regional climate and maintaining biodiversity. In order to understand the resource utilization ability and relationship of fish in Lhasa River, and provide basic data for the protection and sustainable utilization of fish diversity in Lhasa River, the fishery resources in September 2019 (autumn), June 2020 (summer), November 2020 (winter) and March 2021 (spring) were investigated for the first time. The main fish in the catch were determined by the Index of Relative Importance (IRI). Shannon index and Pianka’s index were used to analyze the niche breadth and overlap values of the main fish caught in Lhasa River Basin from three dimensions of time, space and space–time. The results show that a total of 1399 fish of 3 orders and 4 families were collected. Among them, there is one dominant species, Schizopygopsis younghusbandi, and five common species, Schizothorax waltoni, Schizothorax oconnori, Ptychobarbus dipogon, Triplophysa stenura and Carassius auratus. The spatial–temporal niche breadth ranges from 1.395 to 3.628. They are Schizopygopsis younghusbandi, Carassius auratus, Ptychobarbus dipogon, Schizothorax oconnori, Triplophysa stenura and Schizothorax waltoni in a declining order. The temporal niche breadth is generally small, indicating that the survival activities of the main fish are greatly affected by seasonal changes. The fluctuation range of spatial–temporal overlap value is 0–0.535, and only 20% of the total logarithm is meaningful. The species pairs with temporal niche overlap values exceeding 0.300 reached 80%, while the spatial niche overlap value is generally low. The results show that the spatial and temporal distributions of the main fish in Lhasa River were quite different, and the spatial and temporal niche overlap value of the relationship between species is obviously affected by the change in spatial distribution. This demonstrates that the competition relationship and intensity of the main fish are quite different in different spatial positions.