Genetic diversity and population structure of three Schizothorax curvilabiatus populations in the Yarlung Zangbo River based on mitochondrial genes

Background: Schizothorax curvilabiatus is a commercially important freshwater fish endemic to the lower reaches of the Yarlung Zangbo River. However, its fishery resources have been severely depleted due to human activities, leading to a scarcity of data regarding its genetic characteristics. This study aims to assess the genetic diversity and population structure of this species to provide a scientific basis for its conservation.

Methods: Three mitochondrial DNA markers—Cytochrome b (Cyt b), Control Region (D-loop), and Cytochrome Oxidase Subunit I (COI)—were sequenced and analyzed. A total of 100 samples were collected from three natural populations in the lower Yarlung Zangbo River: Motuo (MT), Palong (PL), and Chayu (CY).

Results: High overall haplotype diversity (H_d: 0.889–0.951) but relatively low nucleotide diversity (P_i: 0.00450–0.0111) were observed. Notably, the PL population exhibited the lowest genetic diversity across all markers. Significant genetic differentiation was detected between the PL population and the other two populations (MT and CY), supported by F_ST indices, gene flow (N_m) estimates, and Analysis of Molecular Variance (AMOVA), which revealed that most variation resided within populations (>83%). Both Neighbor-Joining (NJ) phylogenetic trees and haplotype network analyses indicated a distinct geographical lineage pattern, with the PL and CY populations showing a distant genetic relationship.

Discussion: Although S. curvilabiatus maintains substantial overall genetic diversity, the low diversity within specific populations (particularly PL) and the significant inter-population differentiation highlight the vulnerability of its germplasm resources. These findings provide crucial insights into the genetic background of the species and offer valuable guidance for the conservation and sustainable management of its genetic resources.

1 Introduction

Schizothorax curvilabiatus is a freshwater fish species belonging to the order Cypriniformes, family Cyprinidae, subfamily Schizothorinae, and genus Schizothorax. This species is widely distributed in the lower reaches of the Yarlung Zangbo River, including the mainstream of Motuo and its tributaries, Palong Zangbo and Chagyu River, at elevations ranging from 500 to 2300 meters. Historically, it has served as an important contributor to the local economy and remains a dominant fish species in the region (Wu and Wu, 1992; Wu, 1985). Due to the unique conditions of its high-altitude habitat, indigenous fish in Tibet exhibit slow growth rates, prolonged resource replenishment cycles, and late sexual maturity (Liu et al., 2021), with females reaching maturity at approximately 12 years and males at around 6 years (Wang, 2018). In recent years, the wild population of S. curvilabiatus has experienced a significant decline due to various factors, including overfishing for economic purposes, habitat destruction caused by the construction of water and hydropower facilities, and the introduction of alien species (Wang et al., 2022; Zhang et al., 2023). This decline has been accompanied by a noticeable trend of resource depletion, as well as an increase in individual miniaturization and a shift toward younger age groups. There is an urgent need for in-depth basic research to establish a scientific basis for the conservation of its resources. Existing studies have primarily focused on the biological characteristics (Zhang et al., 2021; Liu and Wei, 2022) and muscle nutrient composition (Jin et al., 2022) of S. curvilabiatus, with limited research on the differentiation between geographical populations. Therefore, it is crucial to investigate the germplasm resources of S. curvilabiatus in the lower reaches of the Yarlung Zangbo River. This research should include an analysis of its genetic diversity and structure, as well as an evaluation of its historical dynamics. Such studies are essential for the conservation and sustainable utilization of S. curvilabiatus germplasm resources.

Genetic diversity is a core aspect of species conservation research and a prerequisite for survival, adaptation, and evolution (David et al., 2018; Cao et al., 2024). Mitochondrial DNA (mtDNA), characterized by maternal inheritance, a simple structure, and a rapid evolution rate, has been widely applied in phylogenetic analysis and population genetic studies of fish (He et al., 2022; Xu et al., 2022). Among mtDNA segments, Cytochrome b (Cyt b) evolves at a moderate rate due to selection pressure, making it a frequently used marker for analyzing phylogenetic origins, geographic lineage, and population historical dynamics (Sun et al., 2023; Purohit et al., 2023). The control region (D-loop), a non-coding sequence with the fastest evolution rate in mtDNA, contains abundant genetic variation, thereby providing more comprehensive insights in phylogenetic and population studies (Fang et al., 2022; Roy et al., 2024). The cytochrome c oxidase subunit I (COI) gene, known for its slow evolutionary rate and strong conservation in its 5’ end (500–700 bp), is frequently used in vertebrate species identification and population genetics studies due to its abundant base variation (Liu et al., 2020). Although the mtDNA sequence and structure of S. curvilabiatus have been described (Wang et al., 2017), its population genetic structure and diversity have not been systematically studied. Ma (2019) analyzed genetic differences among populations in the Motuo Section of Yarlung Zangbo River and Yigong Lake (a first-level tributary of the Palong Zangbo River) using the SLAF-seq method and developed SNP markers. However, due to insufficient sample size, a lack of analysis of inter-population gene flow, as well as the omission of the CY population, the genetic information available for S. curvilabiatus remains limited.

In this study, we analyzed the genetic diversity and population genetic structure of 100 samples collected from three sections of the lower reaches of the Yarlung Zangbo River (Motuo(MT), Palong(PL), and Chayu(CY)) using mitochondrial Cyt b, D-loop, and COI sequences. Our aim was to provide a scientific basis for the effective conservation strategies and germplasm resource managing of S. curvilabiatus.

2 Materials and methods

2.1 Sample collection and DNA extraction

In this study, samples of S. curvilabiatus were collected from the MT (35), PL (30) and CY (35) Rivers in the lower reaches of the Yarlung Zangbo River (Table 1; Figure 1). Fin clips were taken from the samples, preserved in anhydrous ethanol (95%), and stored at -20°C. DNA extraction was performed using a commercial DNA extraction kit (Sangon Biotech (Shanghai) Co., Ltd). The purity of the extracted DNA was verified via 1% agarose gel electrophoresis, and its concentration and A_260/280 ratio were measured using a Nanodrop 2000 spectrophotometer (Thermo Scientific, USA).

Table 1

Table 1. Sample collection information of S.curvilabiatus.

Figure 1

Figure 1. Sample sites of S. curvilabiatus. The red points represent MT population; the green points represent PL population; the black points represent CY population.

2.2 Cyt b, D-loop and COI gene sequence primer design and PCR amplification

Primers were designed using Primer 6 software based on the Cyt b, D-loop, and COI sequences (Genebank Accession no. NC_035994.1) from NCBI. Their specificity was verified prior to synthesis by Biomek, as detailed in Table 2. The PCR reaction consisted of a total volume of 20 μL, including 10 μL of 2× Taq PCR Master Mix (with dye), 6.5 μL of ddH₂O, 1 μL each of upstream and downstream primers (at a concentration of 10 μmol/L), and 1.5 μL of template DNA. The PCR cycling conditions were as follows: an initial denaturation at 95°C for 4 minutes, followed by 35 cycles of denaturation at 95°C for 30 seconds, annealing at 55°C for 30 seconds, and extension at 72°C for 40 seconds. A final extension was performed at 72°C for 10 minutes, followed by a hold at 15°C for 5 minutes. The PCR products were subsequently sent to Biomek for sequencing.

Table 2

Table 2. Mitochondrial gene primer sequence information.

2.3 Data analysis

The sequencing results were edited using BioEdit 7.0.5 (Hall, 1999) and manually proofread. Homology comparisons with the Cyt b, D-loop, and COI sequences with mtDNA (Gene bank Accession no. NC_035994.1) from NCBI were conducted using MEGA 7.0 (Kumar et al., 2016) to confirm sequence accuracy. Base composition, sequence variation, and inter-population genetic distances were analyzed. A phylogenetic tree was constructed using the Kimura 2-parameter (K2-P) model via the Neighbor-Joining (NJ) method, with branch confidence evaluated through Bootstrap analysis (1000 replicates). Genetic metrics, including the number of polymorphic sites (S), nucleotide diversity (P_i), number of haplotypes (H), and haplotype diversity index (H_d), were calculated using DnaSP 6.0 (Rozas et al., 2017). The haplotype network was constructed using PopART 1.7 (Leigh and Bryant, 2016). Additionally, Arlequin 3.5 (Excoffier et al., 2007) was used to perform an Analysis of Molecular Variance (AMOVA) to assess genetic variation within and between populations. To investigate potential population expansion, neutrality tests including Tajima’s D (Zhang et al., 2022) and Fu’s Fs (Zhao et al., 2021) were conducted, alongside with nucleotide mismatch distribution analysis.

3 Results

3.1 Sequence characteristics of mitochondrial genes

Following calibration and comparison of the sequencing results, Cyt b (1007 bp), D-loop (886 bp) and COI (1371 bp) sequences were selected for subsequent analysis. Analysis using MEGA revealed a significant AT bias in the mitochondrial gene sequences of the three populations of S. curvilabiatus, with the A+T content considerably exceeding the G+C content (Table 3). This observation is consistent with the typical characteristics of fish mtDNA. Additionally, the Cyt b sequence showed 107 polymorphic loci, consisting of 97 parsimony-informative sites and 10 singleton variable sites. The D-loop sequence contained 59 polymorphic sites, including 42 parsimony informative sites and 17 singleton variable sites. The COI sequence displayed 97 polymorphic sites, with 82 parsimony informative sites and 15 singleton variable sites (Table 4). The transition/transversion ratios for the three sequence regions were 11.4, 7.0 and 7.4, respectively.

Table 3

Table 3. Nucleotide composition of mtDNA Cyt b, D-loop and COI gene sequences of S. curvilabiatus.

Table 4

Table 4. Genetic diversity parameters in mtDNA Cyt b gene, D-loop region and COI gene sequence of three populations of S. curvilabiatus.

3.2 Analysis of population genetic diversity

The genetic diversity parameters of three populations of S. curvilabiatus were analyzed using DnaSP 6.0, and the results are presented in Table 4. Among the 100 samples, a total of 41 Cyt b haplotypes were identified. The H_d for the populations ranged from 0.813 to 0.929, with an overall value of 0.945. The P_i ranged from 0.00551 to 0.01684, with an overall value of 0.0111. The K ranged from 5.55347 to 16.95429, with an overall average of 11.16662. Notably, the CY population exhibited the highest haplotype diversity, while the MT population showed the highest nucleotide diversity and average number of nucleotide differences. For the D-loop sequence, 48 haplotypes were identified. The H_d, P_i and K were 0.951, 0.01191 and 9.353, respectively. The CY population exhibited the highest haplotype diversity, while the MT population had the highest nucleotide diversity and average number of nucleotide differences. Additionally, 31 haplotypes were detected in the COI sequence, with an overall H_d of 0.889 and P_i index of 0.00450. Notably, the CY population also showed the highest haplotype diversity, nucleotide diversity, and average number of nucleotide differences among the three populations.

Comparatively, the MT population had the highest nucleotide diversity and average number of nucleotide differences, while the CY population consistently demonstrated the highest haplotype diversity across all analyses. Furthermore, the genetic parameters derived from D-loop analysis were superior to those obtained from the COI sequence, which reflected the lowest values.

3.3 Analysis of population genetic structure

Statistical analyses of genetic distance and genetic differentiation indices revealed the following: based on the Cyt b gene, genetic distances ranged from 0.00864 to 0.0134, while genetic differentiation indices varied from 0.06293 to 0.20213. For the D-loop gene, genetic distances ranged from 0.0116 to 0.0134, with differentiation indices between 0.08634 and 0.24561. Regarding the COI gene, genetic distances ranged from 0.00238 to 0.00578, and differentiation indices ranged from 0.05183 to 0.16328 (Table 5). The N_m values for gene flow between the MT and PL populations ranged from 2.65 to 4.57, between the MT and CY populations from 1.48 to 2.74 and between the PL and CY populations from 0.77 to 1.28 (Table 6).

Table 5

Table 5. Population genetic differentiation and genetic distances.

Table 6

Table 6. The gene flow among S. curvilabiatus populations.

Results from the AMOVA analysis (Table 7) revealed that the contribution rates within population genetic variation were 88.59% (Cyt b), 83.72% (D-loop) and 91.58% (COI). Conversely, the contribution rates of among population genetic variation were 11.41%, 16.28% and 8.42%, respectively. These findings indicate that the primary source of genetic variation is within the populations. Furthermore, analysis of the genetic differentiation indices showed a significant degree of genetic differentiation among populations, with F_ST values of 0.11414 (Cyt b), 0.16276 (D-loop) and 0.08423 (COI) (P < 0.05).

Table 7

Table 7. Analysis of molecular variances (AMOVA).

3.4 Haplotype distribution and phylogenetic analysis

Based on multilocus sequence data, this study constructed a neighbor-joining (NJ) phylogenetic tree with S. oconnori and S. waltoni as out groups to investigate the population evolutionary relationships of S. curvilabiatus. The results revealed significant phylogeographic structure among populations (Figure 2), demonstrating that the MT population formed a monophyletic clade with the PL and CY populations, while the PL and CY populations each constituted independent subclades. To validate the phylogenetic topology, a haplotype network analysis was performed (Figure 3), showing that the distribution patterns of haplotypes in Cyt b (47 haplotypes), D-loop (48 haplotypes), and COI (31 haplotypes) highly corresponded with the NJ tree topology. Haplotype sharing analysis indicated that the MT-PL population shared 3 (Cyt b), 3 (D-loop) and 2 (COI) haplotypes, while the MT-CY population shared 2 (Cyt b), 1 (D-loop) and 1 (COI) haplotype. Notably, no shared haplotypes were detected between the PL and CY populations.

Figure 2

Figure 2. NJ phylogenetic trees based on Cyt b, D-loop and COI sequences using Kimura 2-parameter (K2-P).

Figure 3

Figure 3. Median-Joining network of Cyt b (A), D-loop (B) and COI (C) haplotypes from three populations of S. curvilabiatus.

3.5 Population historical dynamic information

The results of the neutral tests, Tajima’s D and Fu’s Fs, are presented in Table 8. For the Cyt b sequence, negative Tajima’s D values were observed across all populations, while a positive Fu’s Fs value was found only in the MT population. However, Tajima’s D values for all populations were not statistically significant (P > 0.05). In contrast, the Fu’s Fs values for the PL and CY populations were significantly different (P < 0.05). The overall values for Tajima’s D and Fu’s Fs were -1.48953 (P < 0.05) and -21.13796 (P < 0.01), respectively.

Table 8

Table 8. Neutrality test for three S. curvilabiatus populations.

Analysis using the D-loop sequence indicated that only the CY population had negative values for both Tajima’s D and Fu’s Fs, with the difference in Fu’s Fs being significant (P < 0.05). The D and Fs values for the PL population were positive but not significant, while the MT population exhibited positive D values and negative Fs values, with both results also being non-significant. The overall D and Fs values across the D-loop sequence were -0.35407 (P > 0.05) and -19.13379 (P < 0.05), respectively.

In the analysis of the COI sequence, all three populations had negative Tajima’s D values. The Fu’s Fs values were negative for the PL and CY populations, while the overall Fs values for all three populations were positive. Among these, only the CY population showed a significant difference (P < 0.05). The total values of D and Fs for the COI gene were -2.04701 (P < 0.01) and -13.29356 (P < 0.01), respectively.

Additionally, the nucleotide mismatch distribution curves based on Cyt b and D-loop sequences exhibited multi-peaks, indicating that the populations were in a stable state. The historical dynamic analysis based on these two genes showed that the MT population was in a relatively stable state without evidence of population expansion events (Figures 4, 5). Although there were significant differences in the Fu’s Fs values of the PL population, Tajima’s D values were not significantly different, and the nucleotide mismatch map showed multiple peaks, indicating that no population expansion event had occurred historically (Figures 4, 5).

Figure 4

Figure 4. Mismatch distribution of based on mtDNA Cyt b gene sequences. (A) MT; (B) PL; (C) CY; (D) Total.

Figure 5

Figure 5. Mismatch distribution of based on mtDNA D-loop region sequences. (A) MT; (B) PL; (C) CY; (D) Total.

The mismatch analysis for COI gene indicated population expansion (Table 7), with a single-peak overall nucleotide mismatch map (Figure 6). Among the three populations, only the CY population exhibited significantly negative Tajima’s D and Fu’s Fs values, alongside a double-peak nucleotide mismatch map (Figure 6), suggesting historical expansion followed by stabilization.

Figure 6

Figure 6. Mismatch distribution of based on mtDNA COI gene sequences. (A) MT; (B) PL; (C) CY; (D) Total.

Overall, the historical dynamics of S. curvilabiatus are complex. Based on neutral tests and nucleotide mismatch analysis, it is inferred that the recent expansion events are rare, and the species is currently in a relatively stable state.

4 Discussion

4.1 Genetic diversity of different geographical populations of S. curvilabiatus

Genetic diversity is a crucial prerequisite for assessing the status of population resources and informing conservation efforts. It serves as a key analytical method for studying the origin and evolution of species while providing insights into their evolutionary potential (Ward, 2000; Chen et al., 2020). Environmental selection pressures significantly influence population dynamics. Under low selection pressure, populations are more likely to survive, indicating higher adaptability to their environment. This adaptability is often associated with greater genetic diversity, which helps populations remain resilient to environmental changes. Conversely, low genetic diversity renders populations vulnerable to environmental pressures, potentially leading to bottlenecks or extinction (Manel et al., 2020; Williams and Johnson, 2022; Yao et al., 2022).

In this study, mitochondrial DNA sequences, including the Cyt b gene, D-loop region, and COI gene, were used as genetic markers to analyze the genetic diversity of three populations of S. curvilabiatus. Sequence analysis of the homologous fragments of the Cyt b, D-loop, and COI genes revealed a higher content of adenine (A) and thymine (T) compared to guanine (G) and cytosine (C), exhibiting a significant AT bias and a bias against G/C. This observation aligns with the known nucleotide distribution patterns in the mitochondrial genomes of fish (Tsuji et al., 2020; Choi et al., 2021; Tesfaye et al., 2021; Modeel et al., 2023).

The genetic diversity of fish populations is primarily assessed through H_d and P_i. Based on established classification criteria (H_d = 0.5, P_i = 0.005), genetic diversity can be categorized into four levels (Grant and Bowen, 1998; Bowen et al., 2001), the combination of low H_d - low P_i is indicative of the founder effect, whereas a high H_d- low P_i combination reflects post-bottleneck expansion. A low H_d- high P_i combination suggests geographic isolation and subsequent differentiation, while a high H_d- high P_i combination indicates a stable large population. In this study, the overall H_d and P_i values indicate a pattern of high haplotype diversity coupled with low nucleotide diversity. This pattern likely arises because haplotype diversity accumulates more rapidly, whereas nucleotide diversity increases over a longer time frame. It is speculated that a bottleneck effect or post-colonization expansion occurred within the S. curvilabiatus populations, resulting in limited nucleotide variation accumulation in a short period. Among the three genetic markers, the D-loop exhibited the highest H_d and P_i values, suggesting that it is more sensitive to changes than the Cyt b and COI genes. This increased sensitivity may be attributed to the location of the D-loop in the non-coding region, which tends to evolve faster than coding regions due to the absence of natural selection pressures. Therefore, the D-loop region is particularly well-suited for detecting fine-scale population structure. Nevertheless, owing to the inherent limitations in the amount of genetic information provided by any single locus, the protein-coding genes COI and Cyt b continue to hold significant value for phylogenetic reconstruction and species-level taxonomic assessments. Consequently, a comprehensive understanding of population genetic structure is best achieved through the integrated analysis of genetic variation across multiple, independently evolving loci. This finding is consistent with the mitochondrial gene patterns observed in Cyprinus carpio var. baisenensis and Labeo gonius (Roy et al., 2024; Sumana et al., 2024). When analyzing the genetic diversity among the three populations, both the MT and CY populations exhibited high genetic diversity, while the PL population displayed significantly lower genetic diversity. This pattern may align with the theory that environmental gradients, including altitude, climate, and habitat complexity, influence the dispersal and population genetic structure of plateau species. At the altitudinal limits of a species’ range, diversity and adaptation may be accelerated, potentially leading to a reduction in genetic diversity at higher altitudes (Schiffers et al., 2013; Halbritter et al., 2015).

4.2 Genetic structure and genetic differentiation of different geographical populations of S. curvilabiatus

The genetic structure and differentiation among fish populations were assessed using genetic distance, gene flow (N_m), and genetic differentiation index (F_ST). In this study, the three genes exhibited a consistent distribution trend, arranged from closest to farthest: MT and PL populations, MT and CY populations, and CY and PL populations. This pattern corresponds with the geographical distribution of S. curvilabiatus. Based on genetic distance classifications, distances greater than 0.9 were categorized as genus level, distances over 0.3 as species level, and those exceeding 0.05 as population level (Shaklee et al., 1982). None of the three populations reached the population classification standard, suggesting that they may have originated from a common ancestor. Geographical isolation and differing habitats within the same water restricted gene exchange among the populations, contributing significantly to their differentiation and leading to a genetic structure that mirrors the species distribution (Vilcot et al., 2022; Jia et al., 2024).

The genetic differentiation index (F_ST) is a critical measure for evaluating population differentiation. Classification criteria for F_ST are as follows: 0 < F_ST < 0.05: low differentiation; 0.05 < F_ST < 0.15: moderate differentiation; 0.15 < F_ST < 0.25: high differentiation; F_ST > 0.25: very high differentiation (Wright, 1931, 1978). Analysis of the three genes indicated that the F_ST values of the populations ranged from 0.05 to 0.25. Specifically, the F_ST values between the MT and PL populations were 0.05183 to 0.08634, indicating low differentiation. In contrast, the F_ST values for MT and CY populations ranged from 0.08427 to 0.14482, reflecting moderate differentiation. The values for the PL and CY populations were between 0.16328 to 0.24561, indicating a high level of genetic differentiation.

Gene flow (N_m) is another key metric for evaluating population genetic structure and differentiation. In this study, the N_m value of gene flow between the MT and PL populations ranged from 2.65 to 4.57, significantly exceeding 1, which suggests close geographical proximity and potential channels for gene exchange. This indicates that random mating is likely among individuals of the MT and PL populations (Petit and Excoffier, 2009). The N_m value between the MT and CY populations was 1.48 to 2.74, signifying relatively close geographical distance and the possibility of gene exchange. Conversely, the N_m value for the PL and CY populations was 0.77 to 1.28, indicating substantial genetic differentiation between these groups.

In summary, analyses of genetic distance, genetic differentiation index (F_ST), and gene flow (N_m) among three populations revealed a similar differentiation structure. The MT and PL populations were the closest, exhibiting the lowest F_ST values and the highest N_m values, suggesting low differentiation and a high probability of gene exchange. The MT and CY populations demonstrated moderate differentiation, accompanied by a lower likelihood of gene flow (N_m). In contrast, the PL and CY populations were highly differentiated, with minimal gene exchange. The results of the AMOVA analysis among the three populations indicated that over 83% of the genetic variation was within populations, suggesting that the within-population variation is the primary source of overall genetic diversity.

The phylogenetic tree constructed from three genes and the haplotype network indicated that the three populations formed two distinct geographic clusters: one composed of the MT and PL populations, with mixed haplotypes but no obvious phylogenetic relationship, and the other composed of the MT and CY populations, with no obvious phylogenetic relationship between their haplotypes. This suggests that the PL and CY populations formed distinct geographical populations and lineage structures, consistent with the above genetic differentiation pattern.

Based on the above analysis results, the genetic differentiation pattern of the three populations was closely related to their geographic location. The MT population was distributed in the main stream of the Yarlung Zangbo River, while the PL population was distributed in the Palong River basin, a first-level tributary of the Yarlung Zangbo River. The Palong River flows into the Yarlung Zangbo River in the Yarlung Zangbo Grand Canyon, providing opportunities for gene exchange. The CY population is distributed in the Chayu River, which flows into the Yarlung Zangbo River overseas, geographically distant. There is a theoretical possibility of gene exchange between the MT population in the main stream of the Yarlung Zangbo River and the PL population in the Palong River. However, the analysis results show that the genetic diversity of the CY population is high, and there is essentially no possibility of gene exchange. Thus, it became distinctly differentiated from the other populations.

4.3 Population historical dynamics of different geographical populations of S. curvilabiatus

This study employed neutrality tests for three genes and nucleotide mismatch analysis to infer whether the populations of S. curvilabiatus had experienced expansion events. According to the neutrality analysis, negative and statistically significant values of Tajima’s D and Fu’s Fs indicate that a population deviates significantly from the neutral model, suggesting a historical expansion event (Tajima, 1989). In conjunction with nucleotide mismatch analysis, a unimodal distribution, represented by a single peak in the Poisson distribution, further supports the occurrence of population expansion (Márcia et al., 2013). The overall Tajima’s D and Fu’s Fs values for the three genes varied, suggesting that the entire S. curvilabiatus population may have experienced expansion events. However, notable differences in these values indicated that only the CY population might have undergone distinct population expansion events.

The nucleotide mismatch analysis showed that only the overall distribution of the COI gene exhibited a unimodal pattern. In contrast, the remaining populations displayed multimodal distributions. Specifically, the MT and PL populations indicated a stable state without evidence of recent expansion, whereas the CY population may have experienced a recent expansion event or a bottleneck effect, followed by a period of relative stability. These findings align with results observed in Scophthalmus maximus and Pangasius spp (Firidin et al., 2020; Ha et al., 2020).

While this study has revealed the genetic differences among various geographical populations of S. curvilabiatus from a mitochondrial perspective, the limited genetic information available does not fully explain the reasons behind the current genetic structure of this species. Future research could expand the sample size or delve deeper into genetic information to analyze the timing and influencing factors of genetic differentiation among geographical populations at the genomic level.

In recent years, the fishery resources in the lower reaches of the Yarlung Zangbo River have drastically declined due to frequent geological disasters and human activities. As a representative fish species in this region, S. curvilabiatus is inevitably experiencing population decline and bottleneck effects, leading to a significant reduction in genetic diversity.

5 Conclusion

This study found that the genetic diversity of S. curvilabiatus in the MT and CY populations was relatively high overall, but significant differences existed in the genetic diversity within each population. The three populations exhibited two main phylogenetic relationships, revealing the genetic background of the germplasm resources of S. curvilabiatus in the lower reaches of the Yarlung Zangbo River. Moreover, the results provide valuable insights into the conservation, sustainable development, and utilization of the germplasm resources of S. curvilabiatus.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: NCBI, accession PX806090–PX806188.

Ethics statement

The animal study was approved by Animal Ethics Committee of the Academy of Agriculture and Animal Husbandry of Xizang Autonomous Region. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

QS: Conceptualization, Validation, Investigation, Data curation, Funding acquisition, Methodology, Writing – original draft, Formal analysis, Visualization. KH: Writing – original draft, Investigation, Supervision, Software, Resources, Methodology. PY: Writing – original draft, Investigation, Methodology. CZ: Formal analysis, Resources, Visualization, Validation, Writing – review & editing, Investigation, Supervision, Funding acquisition, Conceptualization.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work supported by the National Key R&D Program of China (Grant no. 2024YFD12007000); the Key Science and Technology Plan project of Tibet Autonomous Region (Grant no. XZ202301ZY0012N) and the financial special project of the Ministry of Agriculture and Rural Affairs (Grant no. CJW2023034).

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2025.1742741/full#supplementary-material

References

Bowen W., Bass A. L., Rocha L. A., Grant W. S., and Robertson D. R. (2001). Phylogeography of the Trumpet fishes (Aulostomus): ring species complex on a global scale. Evol. 55, 1029–1039. doi: 10.1111/j.0014-3820.2001.tb00619.x

PubMed Abstract | Crossref Full Text | Google Scholar

Cao L., Shao W. H., Yi W. J., and Zhang E. (2024). A review of conservation status of freshwater fish diversity in China. J. Fish. Biol. 104, 345–364. doi: 10.1111/jfb.15606

PubMed Abstract | Crossref Full Text | Google Scholar

Chen K., Wang W. F., Chen X. L., Zhu W. L., and Wang H. L. (2020). Analysis of genetic diversity among five populations of Lutraria maxima based on mtDNA Cyt b and D-Loop Sequences. Mar. Fisheries. 42, 149–160. doi: 10.13233/j.cnki.mar.fish.2020.02.003

Crossref Full Text | Google Scholar

Choi H. K., Jang J. E., Byeon S. Y., Kim Y. R., Maschette D., and Chung S. (2021). Genetic diversity and population structure of the Antarctic tooth fish, Dissostichus mawsoni, using mitochondrial and microsatellite DNA Markers. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.666417

Crossref Full Text | Google Scholar

David J. C., Margaret B., and Craig M. (2018). Genetic diversity and conservation units: dealing with the species-population continuum in the age of genomics. Front. Ecol. Evol. 6. doi: 10.3389/fevo.2018.00165

Crossref Full Text | Google Scholar

Excoffier L., Laval G., and Schneider S. (2007). Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform. 23, 47–50. doi: 10.1177/117693430500100003

PubMed Abstract | Crossref Full Text | Google Scholar

Fang D. A., Luo H., He M., Mao C. C., Kuang Z., and Qi H. F. (2022). Genetic diversity and population differentiation of Naked Carp (Gymnocypris przewalskii) revealed by cytochrome oxidase subunit I and D-Loop. Front. Ecol. Evol. 10. doi: 10.3389/fevo.2022.827654

Crossref Full Text | Google Scholar

Firidin S., Ozturk R. C., Alemdag M., Eroglu O., Terzi Y., and Kutlu I. (2020). Population genetic structure of Turbot (Scophthalmus maximus L. 1758) in the black sea. J. Fish Bio. 94, 1154–1164. doi: 10.1111/jfb.14487

PubMed Abstract | Crossref Full Text | Google Scholar

Grant W. A. S. and Bowen B. W. (1998). Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation. J. Hered. 89, 415–426. doi: 10.1093/jhered/89.5.415

Crossref Full Text | Google Scholar

Ha T. T. T., Nga T. T., Hang T. N. A., and Alam M. S. (2020). Genetic diversity in Pangasius spp. collected in Bangladesh based on mitochondrial Cytochrome b gene sequence analysis. Aquacul. Rep. 17, 100351. doi: 10.1016/j.aqrep.2020.100351

Crossref Full Text | Google Scholar

Halbritter A. H., Billeter R., Edwards P. J., and Alexander J. M. (2015). Local adaptation at range edges: comparing elevation and latitudinal gradients. J. Evol. Biol. 28, 1849–1860. doi: 10.1111/jeb.12701

PubMed Abstract | Crossref Full Text | Google Scholar

Hall T. A. (1999). BioEdit: A user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT. Nucleic. Acids Symposium. Ser. 41, 95–98.

Google Scholar

He W. P., Dong Z. L., Ma T. T., Yan H. G., Chen Z. X., and Yao W. Z. (2022). Three mitochondrial markers reveal genetic diversity and structure of Rock Carp (Procypris rabaudi) endemic to the upper Yangtze: implications for pre-release genetic assessment. Front. Mar. Sci. 9. doi: 10.3389/fmars.2022.939745

Crossref Full Text | Google Scholar

Jia M. X., Cheng F., Li J., Schmidt B. V., Li Y. C., and Xie S. G. (2024). Distribution patterns of the two genetic groups of Corbicula fluminea in A lotic-lentic System. Ecol. Evol. 14, e11339. doi: 10.1002/ece3.11339

PubMed Abstract | Crossref Full Text | Google Scholar

Jin H. Y., Li L., Qin D. L., Wu S., Wang N. M., and Jin X. (2022). Nutritional composition in muscle of Schizothorax curvilabiatus. Biotic. Resour. 44, 198–204. doi: 10.14188/j.ajsh.2022.02.009

Crossref Full Text | Google Scholar

Kumar S., Stecher G., and Tamura K. (2016). MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | Crossref Full Text | Google Scholar

Leigh J. W. and Bryant D. (2016). PopART: full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110–1116. doi: 10.1111/2041-210X.12410

Crossref Full Text | Google Scholar

Liu F., Li M. Z., Wang J., Gong Z., Liu M., and Liu H. Z. (2021). Species composition and longitudinal patterns of fish assemblages in the middle and lower Yarlung Zangbo River, Tibetan Plateau, China. Ecol. Indic. 125, 107542. doi: 10.1016/j.ecolind.2021.107542

Crossref Full Text | Google Scholar

Liu K. Y., Zhao S., Yu Z. C., Zhou Y. J., Yang J. Y., and Zhao R. (2020). Application of DNA barcoding in fish identification of supermarkets in Henan province, China: more and longer COI gene sequences were obtained by designing new primers. Food. Res. Int. 136, 109516. doi: 10.1016/j.foodres.2020.109516

PubMed Abstract | Crossref Full Text | Google Scholar

Liu Y. C. and Wei C. (2022). Age and growth characteristics of Schizothorax curilabiatus in the Chayu river, Tibet. Hubei. Agr. Sci. 61, 112–115. doi: 10.14088/j.cnki.issn0439-8114.2022.12.022

Crossref Full Text | Google Scholar

Ma H. X. (2019). Development of SNP markers and population genetics analysis of Schizothorax curvilabiatus based on SLAF-seq technology. Huazhong Agricultural University, Wuhan (IL.

Google Scholar

Manel S., Guerin P. E., Mouillot D., Blanchet S., Velez L., and Albouy C. (2020). Global determinants of freshwater and marine Fish genetic diversity. Nat. Commun. 11, 692–791. doi: 10.1038/s41467-020-14409-7

PubMed Abstract | Crossref Full Text | Google Scholar

Márcia B. A., Real R., Muñoz A. R., and Brown J. A. (2013). New measures for assessing model equilibrium and prediction mismatch in species distribution models. Divers. Distrib. 19, 1333–1338. doi: 10.1111/ddi.12100

Crossref Full Text | Google Scholar

Modeel S., Joshi B. D., Yadav S., Bharti M., and Negi R. K. (2023). Mitochondrial DNA reveals shallow population genetic structure in economically important Cyprinid fish Labeo rohita (Hamilton 1822) from south and southeast Asia. Mol. Biol. Rep. 50, 4759–4767. doi: 10.1007/s11033-023-08386-5

PubMed Abstract | Crossref Full Text | Google Scholar

Petit R. J. and Excoffier L. (2009). Gene flow and species delimitation. Trends. Ecol. Evol. 24, 386–393. doi: 10.1016/j.tree.2009.02.011

PubMed Abstract | Crossref Full Text | Google Scholar

Purohit S., Sharma P., Kothiyal S., Singh U., Nautiyal P., Singh M., et al. (2023). Resolving the phylogenetic relationship of himalayan snow trout Schizothorax plagiostomus with other species of Schizothoracine using mitochondrial COI and Cyt b genes. Mol. Biol. Rep. 50, 3927–3933. doi: 10.1007/s11033-023-08274-y

PubMed Abstract | Crossref Full Text | Google Scholar

Roy S., Behera B. K. V. L. R., Rout A. K., Kumar V., Parida P. K., and Jana A. K. (2024). Genetic characterization of Minor Carp (Labeo gonius) from Indian rivers revealed through mitochondrial ATPase 6/8 and D-loop region analysis: implications for conservation and management. Front. Mar. Sci. 11. doi: 10.3389/fmars.2024.1345649

Crossref Full Text | Google Scholar

Rozas J., Ferrer-Mata A., Sánchez-DelBarrio J. C., Guirao-Rico S., Librado P., and Ramos-Onsins S. E. (2017). DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299–3302. doi: 10.1093/molbev/msx248

PubMed Abstract | Crossref Full Text | Google Scholar

Schiffers K., Bourne E. C., Lavergne S., Thuiller W., and Travis J. M. (2013). Limited evolutionary rescue of locally adapted populations facing climate change. Philos. Trans. R. Soc Lond. 368, 20120083. doi: 10.1098/rstb.2012.0083

PubMed Abstract | Crossref Full Text | Google Scholar

Shaklee J. B., Tamaru C. S., and Waples R. S. (1982). Speciation and evolution of marine fishes studied by the electrophoretic analysis of proteins. Pac. Sci. 36, 141–157.

Google Scholar

Sumana S. L., Wang P. P., Zhang C. F., Jing X. J., Zhu J., Tang Y. K., et al. (2024). Genetic diversity of the common carp black strain population based on mtDNA (D-loop and cyt b). Heliyon 10, e30307. doi: 10.1016/j.heliyon.2024.e30307

PubMed Abstract | Crossref Full Text | Google Scholar

Sun C. H., Huang Q., Gu Y. L., Huang X. Y., Du H. W., and Lao Y. (2023). Genetic diversity and population structure of Saurida undosquamis in the south China sea based on Cyt b and the control region. Glob. Ecol. Conserv. 48, e02740. doi: 10.1016/j.gecco.2023.e02740

Crossref Full Text | Google Scholar

Tajima F. (1989). Statistical method for testing the neutral mutation hypothesis by DNA Polymorphism. Genetics. 123, 585–595. doi: 10.1093/genetics/123.3.585

PubMed Abstract | Crossref Full Text | Google Scholar

Tesfaye G., Curto M., Meulenbroek P., Englmaier G. K., Tibihika P. D., and Alemayehu E. (2021). Genetic diversity of Nile tilapia (Oreochromis niloticus) populations in Ethiopia: insights from nuclear DNA microsatellites and implications for conservation. BMC. Ecol. Evo. 21, 113. doi: 10.1186/s12862-021-01829-2

PubMed Abstract | Crossref Full Text | Google Scholar

Tsuji S., Maruyama A., Miya M., Ushio M., Sato H., and Minamoto T. (2020). Environmental DNA analysis shows high potential as a tool for estimating intra specific genetic diversity in a wild fish population. Mol. Ecol. Resour. 20, 1248–1258. doi: 10.1111/1755-0998.13165

PubMed Abstract | Crossref Full Text | Google Scholar

Vilcot M., Albouy C., Donati G. F. A., Claverie T., Julius P., and Manel S. (2022). Spatial genetic differentiation correlates with species assemblage turnover across tropical reef fish lineages. Global. Ecol. Biogeogr. 32, 535–547. doi: 10.1111/geb.13637

Crossref Full Text | Google Scholar

Wang J. (2018). Fish community structure in lower reach Yarlung Zangbo river and the biological trials of Schizothorax curvilabiatus. University of Chinese Academy of Sciences, Beijing (IL.

Google Scholar

Wang J., Zhang F. B., Hu H. M., Gong Z., Cao W. X., and Lin P. C. (2022). Characteristics of age and growth of Schizothorax curvilabiatus in the lower reaches of the Yarlung Zangbo river. Acta Hydrobiol. Sin. 46, 1770–1779. doi: 10.7541/2022.2021.060

Crossref Full Text | Google Scholar

Wang W. L., Zhang L. H., Pan Y. Z., Lhamo T. S., Zhang C., and Li B. H. (2017). The complete mitochondrial genome of the Schizothorax curilabiatus (Cypriniformes: Cyprinidae). Mitochondrial. DNA. B. 2, 683–684. doi: 10.1080/23802359.2017.1383202

PubMed Abstract | Crossref Full Text | Google Scholar

Ward R. D. (2000). Genetics in fisheries management//marine genetics (Dordrecht: Springer Netherlands Press).

Google Scholar

Williams T. J. and Johnson J. B. (2022). History predicts contemporary community diversity within a Biogeo graphic province of freshwater fish. J. Biogeogr. 49, 809–821. doi: 10.1111/jbi.14316

Crossref Full Text | Google Scholar

Wright S. (1931). Evolution in mendelian populations. genetics. 16, 97–159. doi: 10.1093/genetics/16.2.97

PubMed Abstract | Crossref Full Text | Google Scholar

Wright S. (1978). Evolution and the genetics of populations: variability within and among natural population (Chicago: University of Chicago Press).

Google Scholar

Wu Y. F. (1985). A preliminary analysis of fish fauna of Namjag barwa region in Xizang, China. Acta Bio. Plateau. Sin. 4, 61–70.

Google Scholar

Wu Y. F. and Wu C. Z. (1992). The fishes of Qinghai-Xizang plateau (Chengdu: Sichuan Science& Technology Press).

Google Scholar

Xu Z., Dou S. Z., Ding S. X., and Liu J. X. (2022). Temporal genetic stability despite decades of over exploitation for large yellow croaker in the East China Sea. Front. Mar. Sci. 9. doi: 10.3389/fmars.2022.861840

Crossref Full Text | Google Scholar

Yao M., Zhang S., Lu Q., Chen X. Y., Zhang S. Y., and Kong Y. Q. (2022). Fishing for fish environmental DNA: ecological applications, methodological considerations, surveying designs, and ways forward. Mol. Ecol. 31, 5132–5164. doi: 10.1111/mec.16659

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang J. L., Huang J. Q., Fang C., Xu T. Q., and Wang K. F. (2021). Study on the biological characteristics of Schizothorax curvilabiatus. Fish. Hebei. 1, 11–16. doi: 10.3969/j.issn.1004-6755.2021.01.003

Crossref Full Text | Google Scholar

Zhang W. Y., Jiang S. F., Salumy K. R., Xuan Z. Y., Xiong Y. W., and Jin S. B. (2022). Comparison of genetic diversity and population structure of eight Macrobrachium nipponense populations in China based on D-loop sequences. Aquacult. Rep. 23, 101086. doi: 10.1016/j.aqrep.2022.101086

Crossref Full Text | Google Scholar

Zhang C., Wang W. L., and Zhou J. S. (2023). Primary color map and mitochondrial genome of fish in Yarlung Zangbo river basin (Beijing: China Agriculture Press).

Google Scholar

Zhao Y., Zhu X., Jiang Y., Li Z., Li X., and Xu W. B. (2021). Genetic diversity and variation of seven Chinese grass shrimp (Palaemonetes sinensis) populations based on the mitochondrial COI gene. BMC. Ecol. Evo. 21, 167–176. doi: 10.1186/s12862-021-01893-8

PubMed Abstract | Crossref Full Text | Google Scholar