Issue |
Aquat. Living Resour.
Volume 37, 2024
|
|
---|---|---|
Article Number | 5 | |
Number of page(s) | 12 | |
DOI | https://doi.org/10.1051/alr/2024003 | |
Published online | 11 April 2024 |
Research Article
Genetic structure of Macrobrachium nipponense, an important farmed freshwater shrimp in China, in the Three Gorges Reservoir
1
Hubei Engineering Research Center for Protection and Utilization of Special Biological Resources in the Hanjiang River Basin, School of Life Sciences, Jianghan University, Wuhan, China
2
Hubei Key Laboratory of Environmental and Health Effects of Persistent Toxic Substances, Jianghan University, Wuhan, China
3
College of Environmental Science and Engineering, China West Normal University, Nanchong, China
* Corresponding author: zhaidongdong@jhun.edu.cn
Received:
30
May
2023
Accepted:
7
March
2024
Macrobrachium nipponense is a major farmed species freshwater shrimp in China and its genetic diversity is of high value in aquaculture. The construction of the Three Gorges Project has slowed down the water velocity in the Three Gorges Reservoir (TGR) and made the water more clear, which is more favourable for the survival and reproduction of M. nipponense. However, there are lack of studies on the genetic diversity of M. nipponense populations in the TGR. In this study, mitochondrial COI gene sequences were used as molecular marker to analyze the genetic diversity, genetic differentiation and historical demography of eight M. nipponense populations in the TGR. The results showed that the M. nipponense populations in TGR have high genetic diversity. There was some genetic differentiation between the four populations in lower reaches of the TGR and four populations in the upper reaches of the TGR. We also found that there was a significant correlation between genetic distance and geographical distance among the M. nipponense populations in the TGR. The M. nipponense population in the TGR has experienced population reduction in the recent year, which might be related to the ice age movement in the Little Ice Age and human activities. This study provides a theoretical basis for the rational development and protection of M. nipponense resources in the TGR.
Key words: Three Gorges Reservoir / Macrobrachium nipponense / genetic diversity / population genetic structure
Handling Editor: Carlos Saavedra
© D. Wang et al., Published by EDP Sciences 2024
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Shrimp and crab are important aquaculture species in China, whose annual yield is second only to fish (Bureau of Fishery, Ministry of Agriculture and Affairs of the P.R.C, 2020). M. nipponense belongs to Macrobrachium in Palaemonidae, which is one of the most widely cultured freshwater shrimps in China. In 2019, the cultured yield of M. nipponense in China was up to 225,300 tons, second only to that of Litopenaeus Vannamei and Procambarus clarkii (Bureau of Fishery, Ministry of Agriculture and Affairs of the P.R.C, 2020). Improved variety play an important role in aquaculture (Sae-Lim et al., 2017). Continuous cultivation of new variety is an important guarantee for the healthy and sustainable development of aquaculture (Janssen et al., 2017). Wild germplasm resources carrying genes of potential value in improving yield, quality, and adaptation to the environment provide a wide range of raw materials for improved farmed varieties (Pathirana and Carimi, 2022). Therefore, protecting wild germplasm resources plays an important role in aquaculture (Hill et al., 2021; Janssen et al., 2018; Wang et al., 2021).
A comprehensive survey of the genetic diversity of M. nipponense can provide guidance for better utilization of germplasm resources for genetic improvement (Xiong et al., 2023). As an important farmed shrimp in China, there have been a large number of studies on the genetic diversity of M. nipponense populations in major producing areas or important water areas. For example, Feng et al. (2008) used mitochondrial COI gene sequence as a molecular marker study the genetic diversity of M. nipponense populations in the five freshwater lakes of China (i.e. Poyang Lake, Dongting Lake, Taihu Lake, Chaohu Lake, and Hongze Lake), they found that each populations has good genetic diversity, and there was a large genetic differentiation among each populations. And Feng et al. (2021) also used mitochondrial COI gene sequence as a molecular marker to investigate the genetic diversity of M. nipponense from four water systems in Henan Province, China. High haplotype diversity and low nucleotide diversity were observed, and there was no significant genetic differentiation for M. nipponense from the four water systems. The demographic history analysis showed that M. nipponense of Henan province did not experience recent population expansion. There are also many studies using microsatellite markers (SSR), for example, Ma et al. (2012) used nine polymorphic microsatellite markers to analyze the genetic diversity and structure of four M. nipponense populations in Qiandao Lake, China, and found the four populations displayed high genetic diversity. Mutation-drift equilibrium analysis showed no significant bottleneck effect. A moderate level of differentiation among these populations. However, most studies on the genetic diversity of M. nipponense were conducted in natural water areas, and few studies were conducted in areas that has been subjected to strong transformations by human activities, such as reservoirs.
The Yangtze River is the longest river in Asia and has abundant hydroenergy. The Three Gorges Dam was built in Sandouping, Yichang, in 1993 in order to solve the flood problem in the middle reaches of the Yangtze River and develop the hydroenergy resources of the Yangtze River (New and Xie, 2008). The Three Gorges Dam is 181 meters high. After the construction of the Three Gorges Dam, a 663 km long reservoir has been formed from Zigui to Jiangjin in the Yangtze River, with a water area of about 1,084 km2 (Xiang et al., 2021). After the impounding of the TGR, water flow velocity slowed down obviously, water surface area increased, providing more habitat for M. nipponense (Xiang et al., 2021b; Yang et al., 2018). At the same time, the slower flow rate of water makes sediment deposition, water clearing and algae productivity increased, which providing more food sources for M. nipponense and was conducive to the proliferation of M. nipponense (Long et al., 2011; Mirzajani et al., 2020; Xu et al., 2011). As we thought the change of hydrological environment had increased the population number of M. nipponense in the TGR (Yang et al., 2018). A comprehensive understanding of the genetic structure of M. nipponense populations in the TGR is conducive to the conservation and rational development of M. nipponense in the TGR (Alal et al., 2021).
At present, there are a few studies on M. nipponense in the TGR. Chen Sibao et al. (2020) studied the reproductive biology of M. nipponense population in the TGR and found that the spawning time of M. nipponense in the TGR was from late April to early October (Chen et al., 2020). In terms of genetic diversity, only Fu et al. (2010) used microsatellite (SSR) markers to study the genetic diversity of M. nipponense populations in the Yangtze River Basin in 2007, and designed two sample sites (Wanzhou and Chongqing) in the TGR. Fu et al. (2010) found that the genetic diversity of Wanzhou and Chongqing populations were at high level, and the degree of genetic differentiation between Chongqing and Wanzhou populations was low. But there are only two sample sites involved in Fu et al. (2010), which failed to fully reveal the genetic diversity pattern of M. nipponense populations in the TGR.
Due to the advantages of strict maternal inheritance and high mutation rate, mitochondrial DNA have been widely used in population genetics research. Therefore, in this study, partial sequences of mitochondrial cytochrome oxidase subunit I (COI) genes were used as molecular markers to analyze the genetic diversity, inter-population genetic differentiation, and historical demography of eight M. nipponense populations in the head, middle and tail of the TGR. In order to comprehensively understand the genetic diversity and genetic structure of M. nipponense populations in the TGR, and to provide theoretical basis for the protection and rational development of wild germplasm resources of M. nipponense populations in the TGR.
2 Materials and methods
2.1 Sample collection
M. nipponense samples were collected from August 2020 to June 2022 at eight sites in the TGR of the Yangtze River, including Gudongkou (GDK), Xiakou (XK), Shuitianba (STB), Wushan (WS), Wanzhou (WZ), Fengdu (FD), Mudong (MD) and Jiangjin (JJ). GDK, XK and STB are sample points in the tributary of the head area of the reservoir; WS is the sample point in the main stream of the head area of the reservoir; WZ and FD are the sample points in the main stream of the middle area of the reservoir; MD and JJ are sample points in the tributary of the tail area of the reservoir (Fig. 1 and Tab. 1). The experimental samples were collected using cage net, and species identification was carried out with reference to Fauna Sinica (Li et al., 2007). A total of 229 M. nipponense samples were collected from the eight sampling sites in the TGR. A small amount of abdominal muscle was cutted from M. nipponense and stored in 95% alcohol at −20°C in School of Life Sciences, Jianghan University.
![]() |
Fig. 1 Sampling sites of M. nipponense in the TGR. |
Sampling sites and number of individuals of M. nipponense in the TGR.
2.2 DNA extraction, amplification and sequencing
Total DNA was extracted using the Animal Tissues DNA Isolation Kit (Forge Biotechnology Co., Ltd, Chengdu China). The amplification primers for the mitochondrial gene COI fragment sequence were the LCO1490 5’-GGT CAA ATC ATA AAG ATA TTG G-3’ and HCO2198 5’-TAA ACT TCA GGG TGA CCA AAA AAT CA-3’ (Cui et al., 2018). The total reaction volume for the polymerase chain reaction (PCR) was 30 µL, containing PrimeSTAR Max Premix (2×) 15 µL (containing dNTPs, MgCl2, Taq DNA polymerase, and reaction buffer; Takara Bio., Ltd), 2.5 µL each of forward and reverse primers (10 mmol/L), 3 µL template DNA, and 7 µL ddH2O. The PCR amplification procedure was as follows: pre-denaturation at 97 °C for 3 min; followed by 35 cycles of denaturation at 98 °C for 15 s, annealing at 46 °C for 15 s, extension at 72 °C for 20 s; final extension at 72 °C for 5 min and storage at 4 °C. After the PCR products were subjected to agarose gel (1%) electrophoresis, the samples with accurate and clear target bands were sent to Wuhan Tianyi-Huayu Gen Sci Tech Co., Ltd. for purification, recovery, and sequence determination.
2.3 Sequence analysis
The raw sequences were checked and edited with reference to the .ab1 chromatogram files, and the sequences were aligned using MEGA 7 software (Kumar et al., 2016). The population genetic diversity indices calculated included the number of segregating sites (s), haplotype number (h), haplotype diversity (Hd), and nucleotide diversity (Pi) were calculated using DnaSP 5.10 software (Librado and Rozas, 2009).
Analysis of molecular variance (AMOVA) was performed on the M. nipponense population in the TGR to quantify the sources of genetic variation in Arlequin 3.0 software, 229 M. nipponense were divided into eight groups according to geographical location (Excoffier et al., 2005). Pairwise genetic differentiation index (FST values) among 8 M. nipponense populations with the significance level examined by 1000 permutations were calculated in Arlequin 3.0 (Excoffier et al., 2005). STRUCTURE analysis was applied to infer the number of genetically differentiated clusters (K) using STRUCTURE v2.3.4 (Pritchard et al., 2000). The parameters were set as follows: a burn-in period of 100,000 iterations followed by 500,000 recorded iterations for K = 1 to K = 8 clusters and 15 iterations per K values (Pritchard et al., 2000). The most probable number of clusters present in this dataset was determined using the Evanno’s ΔK approach using Structure Harvester online (Earl and Vonholdt, 2012; Evanno et al., 2005) (available at http://taylor0.biology.ucla.edu/structureHarvester/).
Additionally, the genetic distance based on Kimura 2-parameter (K2P) model was calculated using MEGA 7 software to construct a genetic distance matrix between M. nipponense populations (Kumar et al., 2016). The distance along the river between each sampling point was measured using the LocaSpace Viewer 4.33 to construct a geographic distance matrix between M. nipponense populations in the TGR. The Mantel test was performed using the “vegan” package in R software in order to identify the relationship between the genetic distance matrix and geographical distance matrix of M. nipponense populations in the TGR.
To detect whether there was distinct lineage differentiation among haplotypes or populations, a phylogenetic tree of haplotypes was reconstructed using the Bayesian (BI) approaches based on Hasegawa-Kishino-Yano (HKY) model plus gamma distribution rate (+G) plus evolutionarily invariable (+I) model (The best-fit model was calculated using the Bayesian Information Criterion (BIC) in jModeltest 2.1.7) in MrBayes 3.2.7 (Darriba et al., 2012; Ronquist et al., 2012), with Macrobrachium tratense (MW845636), Macrobrachium villosimanus (MW845638) and Macrobrachium yui (MW845641) as the outgroup. Four chains (three hot, one cold) were run for 108 generations, the last 75% generations were used do phylogenetic analyses. A haplotype network was constructed using Median-joining in Popart 1.7 software (Leigh and Bryant, 2015).
Mismatch distribution analyses using DnaSP 5.10 software and neutrality tests (Tajima’s D test and Fu’ FS test) using Arlequin 3.0 software were conducted to test for historical population expansion of the M. nipponense populations in the TGR (Excoffier et al., 2005; Librado and Rozas, 2009). The Bayesian skyline plot (BSP) was constructed using BEAST v2.7.0 software to estimate the historical demography of the M. nipponense population in the TGR (Suchard et al., 2018). The best-fit model for each population was calculated using the BIC in jModeltest 2.1.7 software (Darriba et al., 2012). Referring to the study by Feng et al. (2021), the base evolution rate of the COI gene use 1.4% per million years, which was estimated based on the divergence of sister species of the Palaemonidae caused by the Isthmus of Panama (Knowlton and Weigt 1998). And using a strict molecular clock model with a Markov chain Monte Carlo (MCMC) length of 1.5 × 109 generations, and finally, Tracer 1.5 software was used to detect the results of an effective sample size greater than 200 (ESS > 200), which were output as Bayesian skyline plots (Rambaut and Drummond, 2007).
3 Results
3.1 Sequence characterization and genetic diversity
A total of 229 mitochondrial COI gene sequences of M. nipponense from the TGR were obtained in this study, and the clearly sequenced 635-bp fragments were selected for analysis. There were no base insertions or deletions were observed. There were 51 variable sites in the 229 sequences, accounting for 8.03% of the total number of sites, including 10 singleton variable sites, and 41 parsimony informative sites. The mean total nucleotide composition of 229 sequences was A = 29.8%, T = 29.3%, G = 18.4%, C = 22.5%, that showed an AT bias, which was consistent with the base composition of other crustaceans mitochondrial genome (Wei, 2022).
A total of 46 haplotypes (GenBank accession numbers: OQ978506-OQ978551) were identified in this study, and each sampled population has 8–16 haplotypes. The overall haplotype diversity was 0.801, with the lowest haplotype diversity of 0.525 found in the GDK population and the highest haplotype diversity of 0.926 found in the FD population. The overall nucleotide diversity was 0.01540, with the lowest nucleotide diversity of 0.00932 found in the GDK population and the highest nucleotide diversity of 0.01829 found in the FD population (Tab. 2).
Genetic diversity of eight M. nipponense populations in TGR.
3.2 Population structure
Results of AMOVA demonstrated that genetic variations mainly originated from within populations (92.96%), while variation among populations was relatively low (7.04%) (Tab. 3). The fixation index Fst was 0.07 (0.01 < P < 0.05), indicating that there might be a slight degree of genetic differentiation among the eight M nipponensis populations in the TGR.
The pairwise FST values ranged from 0.0667 to 0.2529, among which the FST values between the GDK and FD populations were the largest, and the FST values between FD and JJ populations were the smallest (Sewall 1978) (Fig. 2). The FST values between populations showed that the four populations in the lower reaches of the TGR (GDK, XK, STB and WS) has widespread genetic differentiation among the four populations in the upper reaches of the TGR (WZ, FD, MD and JJ), and this suggests that the FST values may be related to the geographical distance between M. nipponense populations.
The results of STRUCTURE analysis of 229 M. nipponense in the TGR showed that ΔK reached a maximum value when K = 6. Figure 3 shows the distribution of genetic clusters of each M. nipponense population in the TGR when K = 2–8. No matter what the value of K is, we can clearly see that the proportion of genetic clusters between the four populations in lower reaches of the TGR and the four populations in the upper reaches of the TGR is different, but the number of genetic clusters in each population is the same. This is consistent with the results of pairwise FST values.
The genetic distance (GD) among the eight populations ranged from 0.96% to 1.92%, with the biggest genetic distances was observed between FD and XK (Fig. 2) and the result of the Mantel test was 0.6078 with P value of 0.014. This indicates that there is a significantly correlation between genetic distance and geographical distance of M. nipponense populations in the TGR.
AMOVA of eight M. nipponense populations based on mitochondrial COI gene.
![]() |
Fig. 2 Pairwise Fst values (above the diagonal) and genetic distance (GD, below the diagonal) based on mitochondrial COI gene among eight M. nipponense populations in the TGR. All results were tested for significance 1000 times, and the value marked by a red asterisk (*) represents a passing significance test. |
![]() |
Fig. 3 STRUCTURE analysis of eight M. nipponense populations in the TGR based on COI gene. |
3.3 Haplotype phylogenetic relationship
The haplotype BI phylogenetic tree of M. nipponense in the TGR did not show obvious clustering trend with high support rate, and the private haplotypes of each population did not cluster according to geographical distribution (Fig. 4).
The haplotype network of M. nipponense was roughly overlaps with the BI phylogenetic tree, and the distribution of haplotypes did not show obvious phylogeographic pattern. The haplotypes were closely related to each other, and the mutation steps between all haplotypes are very few (generally only 1–2 step). The network shows that Hap4, Hap10 and Hap17 are abundant and located at the base of the BI phylogenetic tree. It is implied that these haplotypes may be ancestral haplotypes of M. nipponense population in the TGR (Fig. 5).
![]() |
Fig. 4 BI phylogenetic tree of M. nipponense haplotypes constructed based on COI gene. |
![]() |
Fig. 5 The network diagram of M. nipponense haplotypes constructed based on COI gene. Each circle represents a haplotype, and the size of the circle is proportional to the number of individuals in the population. |
3.4 Historical demography
The P values of Tajima’s D-test and Fu’s FS-test for the eight M. nipponense populations in the TGR were not significantly (P < 0.05), indicating that the variation in all M. nipponense populations did not significantly deviate from neutral mutation (Tab. 4). The mismatch distribution analyse plots for each M. nipponense population and for the total population were all significantly multipeaked, again suggesting that no population expansion occurred historically in the M. nipponense population from the TGR (Fig. 6). The Bayesian skyline plot indicated that all populations have experienced population contraction in the recent years. Effective population size of all M. nipponense populations in GTR began to decrease by about 700 years ago, and it was reducing faster and faster, and before that, the effective population size of all populations was largely stable (Fig. 7).
Neutrality tests for M. nipponense population based on mitochondrial COI gene.
![]() |
Fig. 6 Mismatch distribution analyse for each M. nipponense population from different localities. |
![]() |
Fig. 7 Bayesian skyline plot of the effective population sizes through time for M. nipponense population. The X axis is the time scale in units of thousands of years, and the Y axis is the estimated effective population size. |
4 Discussion
4.1 Genetic diversity and historical demography
In broad terms, germplasm is the diversity of a farmed species and its wild relatives that can hybridise and produce fertile progeny (Pathirana and Carimi, 2022). And a comprehensive survey of genetic diversity can provide guidance for better utilization of germplasm resources for genetic improvement (Fu, 2015). Haplotype diversity and nucleotide diversity are important indicators that are often used to evaluate biological genetic diversity. The haplotype diversity of the M. nipponense population in the TGR was 0.801, and the nucleotide diversity was 0.01540, indicating that the genetic diversity is high (Grant and Bowen, 1998). The genetic diversity of the M. nipponense population in the TGR was not only higher than that of the M. nipponense population in the Henan Province (Hd = 0.78574, Pi = 0.01059) but also higher than that of the M. nipponense populations in freshwater lakes in China such as Dongting Lake (Hd = 0.716, Pi = 0.01151), Chaohu Lake (Hd = 0.795, Pi = 0.00494) and so on (Feng et al., 2008, 2021). The high genetic diversity of M. nipponense suggesting that there may be a large number of germplasm resources with excellent quality in the TGR, especially the WZ and FD populations, which are suitable for breeding work.
Grant and Bowen (1998) showed that if a population has high genetic diversity (Hd > 0.5, Pi > 0.005), then it may be due to a long evolutionary history in a large stable population. The haplotype network of M. nipponense in the TGR does not show a simple star-shaped structure. The neutrality tests and mismatch distribution analyses showed that the population size of M. nipponense in the TGR may remained stable in history, indicating that M. nipponense in the TGR may not have experienced severe bottleneck effect or population expansion (Wang et al., 2023). A large number of genetic mutations had been accumulated in M. nipponense population in the TGR, resulting in high genetic diversity in the long evolutionary history (Chen et al., 2017; Chiang and Schaal, 1999; Wang et al., 2000). Secondly, larger population size also favored the maintenance of high genetic diversity in M. nipponense population. Genetic drift can reduce population genetic diversity, and this effect is more significant in small populations. M. nipponense is the dominant species of Crustacea in the TGR with a large population, and little influence from genetic drift. The large population size of M. nipponense in the TGR is conducive to keep its high genetic diversity.
Bayesian skyline plot indicated that M. nipponense populations in the TGR began to shrink 700 years ago, which may be related to the Little Ice Age that began in the 13th century and ended in the early 20th century (Zhu, 1973). The Earth’s climate alternates between ice ages and warm periods, with aquatic organisms generally experiencing population contraction or bottlenecking during colder ice ages and population expansion during warmer interglacial periods (Hewitt, 2000, 1996; Li 2015). The Little Ice Age is the nearest ice age, and the average temperature during it was approximately 1–2 °C lower than present temperature, even in some places the minimum temperature may be 5–7 °C lower than the modern, which may have contributed to the contraction of M. nipponense populations in the TGR (Chang et al., 2020; Zhou et al., 2019; Zhu, 1973). However, in the 20th century, with the end of the Little Ice Age and the gradual warming of the global climate. The results of BSP show that effective population size of M. nipponense in the TGR was still shrinking. Fishing intensity increased may have been a major factor in the decline of effective populations of M. nipponense in the TGR after the 20th century (Jackson et al., 2001; Limburg and Waldman, 2009). Wild fishing was the main source of M. nipponense in the market before large-scale farming, and the fishing intensity of M. nipponense gradually increased with the improvement of fishing technology in the middle of the last century (Chen et al., 2011; Li, 2020). However, the genetic diversity of all M. nipponense populations in the TGR is still at high level, suggesting that population contraction at the Little Ice Age had a limited impact on the genetic diversity of M. nipponense in the TGR.
4.2 Population structure
Results of AMOVA suggest that genetic differentiation may exist among the eight M. nipponense populations in the TGR. The Mantel test indicated that there was a significantly positive correlation between genetic distance and geographical distance, and the pairwise FST values also indicate widespread genetic differentiation between the head area and tail area populations of the TGR. In addition, it is clearly show from the results of STRUCTURE analysis, that the proportion of genetic clusters between the four populations in lower reaches of the TGR and the four populations in the upper reaches of the TGR is a little bit different. This is similar to the study of Feng et al. (2008) on genetic diversity of M. nipponense in five major freshwater lakes in China and Cui et al. (2018) on genetic diversity of M. nipponense in the Huaihe River Basin in China. Genetic differentiation is common among M. nipponense populations whose distance are far, and the degree of genetic differentiation is related to geographical distance to a certain extent.
The weak mobility of M. nipponense may be the cause of genetic differentiation between the four populations in lower reaches of the TGR and four populations in the upper reaches of the TGR. Firstly, M. nipponense is small in size and poor in locomotion ability, which makes it difficult to carry out long-distance migration. Second, the female will stick the fertilized eggs to her abdomen after spawning (Ma et al., 2012). This behavior can improve the survival rate of fertilized eggs, but it prevents the M. nipponense from expanding gene exchange among different populations like species that produces drifting eggs such as Hypophthalmichthys molitrix and Aristichthys nobilis (Wang et al., 2003). In conclusion, because of the long geographical distance between head area populations and tail area populations in the TGR, there is little gene flow between them, so resulting in genetic differentiation (Joanna et al., 2015; Wang et al., 2003). Therefore, M nipponensis in the TGR cannot be simply regarded as a conservation unit, and the genetic structure of M nipponensis in the TGR needs to be continuously investigated.
The absence of genetic differentiation between the GDK population and other populations, despite the barrier imposed by the Gudongkou Hydropower Station, may be related to the relatively short period of time since the Gudongkou Hydropower Station was built. A sufficient number of isolated generations is an important factor in generating genetic differentiation, Ruzich et al. (2019) have found that the effects of dam barrier on the genetic diversity and differentiation of fish populations may not be reflected in a short period of time, but 40–60 generations are sufficient to show. The Gudongkou Hydropower Station cut off links between GDK groups and others in 1995, less than 30 years, and thus the GDK population may not have accumulated sufficient genetic differences from other populations in the head area of the TGR (Cui et al., 2018; Gong 2006).
5 Conclusions
In this study, COI gene sequences were employed to examine the genetic diversity, genetic differentiation and historical demography of eight M. nipponense populations in the TGR. The results showed that M. nipponense population in the TGR had high genetic diversity, so germplasm resources are abundant, and the screening work of excellent germplasm should be carried out to make sure the germplasm resources in the TGR can be used to breed, especially the FD and WZ populations because of their high genetic diversity. In addition, there was some genetic differentiation between the four populations in lower reaches of the TGR and the four populations in the upper reaches of the TGR. Therefore, we advise to segment two management units of M. nipponense in the TGR when we management or develop them. It should be noted that studies based only on mitochondrial genes as markers may have some limitations. Markers such as SSR and SNP should be used to analysis the genetic diversity and genetic structure of M. nipponense populations in the TGR in the future.
Supplementary material
Table 1. Pairwise Fst values (above the diagonal) and corrected P values (below the diagonal) based on mitochondrial COI gene among eight M. nipponense populations in the TGR.
Table 2. Genetic distance (GD) based on mitochondrial COI gene among eight M. nipponense populations in the TGR.
Table 3. Geographical distance (GD, above the diagonal) among eight M. nipponense populations in the TGR.
Access hereAcknowledgments
We are grateful to Yang He for his help in M. nipponense sampling. We are grateful to Wei Zhang for her help in data analysis. We are grateful to Yingxue Peng for her help in laboratory work. We thank AJE (www.aje.cn) for its linguistic assistance during the preparation of this manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (51979123), the Scientific Research Project of Jianghan University (Grant numbers 2021yb104), the Student Scientific Research Project of Jianghan University (Grant numbers 2022zd064), the Innovative Research Team Foundation of the Department of Education of Hubei Province, China (T2020034), and the Scientific Research Project of Jianghan University (2021kjzx006).
Conflicts of interest
All the authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Data availability statement
The original data of COI marker have been deposited in Figshare: Wang, Dong (2023). ALL COI sequence.fas. figshare. Dataset. https://doi.org/10.6084/m9.figshare.22817753.v1.
Author contribution statement
Dong Wang performed the laboratory work, and wrote the manuscript. Le Hu, Fubin Zhang and Hongyan Liu helped in data analysis. Fengqun Zheng and Mengyu Gong helped in laboratory work. Fei Xiong designed the study and revised the manuscript. Dongdong Zhai helped in revised the manuscript. All authors read and approved the final version of the manuscript. Dongdong Zhai also contributed to sampling. All authors read and approved the final manuscript.
References
- Alal GW, Barasa JE, Chemoiwa EJ, Kaunda‐Arara B, Akoll P, Masembe C. 2021. Genetic diversity and population structure of selected lacustrine and riverine populations of African catfish, Clarias gariepinus (Burchell, 1822), in Kenya. J Appl Ichthyol 37: 427–438. [CrossRef] [Google Scholar]
- Bureau of Fishery, Ministry of Agriculture and Affairs of the P.R.C. 2020. China fishery statistical yearbook. China Agricultural Press. [Google Scholar]
- Chang J, Zhang EL, Luo W, Liu EF, Yang XQ, James S. 2020. A quantitative temperature reconstruction of the ‘Little Ice Age’ in southern China. The Holocene 30: 709–720. [CrossRef] [Google Scholar]
- Chen LS, Li JL, Liu QG. 2011. Development history of Macrobrachium nipponense industry in China. Fish Sci Technol Inf 38: 35–39. [Google Scholar]
- Chen PC, Shih CH, Chu TJ, Lee YC, Tzeng TD. 2017. Phylogeography and genetic structure of the oriental river prawn Macrobrachium nipponense (Crustacea: Decapoda: Palaemonidae) in East Asia. PLOS ONE 12: e0173490. [CrossRef] [PubMed] [Google Scholar]
- Chen SB, Liao CS, Rypel AL., Lian YX, Ye SW, Zhang TL, Liu JS. 2020. Spatial and interspecific comparisons of the reproductive biology of two species of co-occurring freshwater shrimps (Decapoda: Caridea: Palaemonidae) in the Three Gorges Reservoir, China. J Crustac Biol 40: 316–324. [CrossRef] [Google Scholar]
- Chiang TY, Schaal BA. 1999. Phylogeography of North American populations of the moss species Hylocomium splendens based on the nucleotide sequence of internal transcribed spacer 2 of nuclear ribosomal DNA. Mol Ecol 8: 1037–1042. [CrossRef] [Google Scholar]
- Cui F, Yu Y, Bao F, Wang S, Xiao MS. 2018. Genetic diversity analysis of the oriental river prawn (Macrobrachium nipponense) in Huaihe River. Mitochondrial Dna Part A 29: 737–744. [CrossRef] [PubMed] [Google Scholar]
- Darriba D, Taboada GL, Doallo R, Posada D. 2012. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods 9: 772–772. [CrossRef] [PubMed] [Google Scholar]
- Earl D, Vonholdt BM. 2012. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour 4: 359–361. [CrossRef] [Google Scholar]
- Evanno G, Regnaut S, Goudet J. 2005. Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol. 14: 2611–2620. [CrossRef] [PubMed] [Google Scholar]
- Excoffier L, Laval G, Schneider S. 2005. Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinforma 1: 47–50. [Google Scholar]
- Feng JB, Sun YN, Cheng X, Li JL. 2008. Sequence analysis of mitochondrial COI gene of Macrobrachium nipponense from the five largest freshwater lakes in China. J Fish China 32: 517–525. [Google Scholar]
- Feng MX, Tang YT, Liu RY, Yan X, Hu B, Nie GX, Zhou CJ. 2021. Genetic diversity of Macrobrachium nipponense in Henan Province assessed using mitochondrial COI gene. J Fish Sci China 28: 1548–1558. [Google Scholar]
- Fu HT, Li FJ, Wu Y, Gong YS, Jiang SF, Xiong YW, Wang N. 2010. Genetic diversity of Macrobrachium nipponense on a regional scale in the Yangtze River. J Fish China 34: 204–212. [CrossRef] [Google Scholar]
- Fu YB. 2015. Understanding crop genetic diversity under modern plant breeding. Theor Appl Genet 128: 2131–2142. [CrossRef] [PubMed] [Google Scholar]
- Gong JH. 2006. Treatment for spillway foundation of step No.1 hydropower station of Gudongkou cascade. Hydropower New Energy 2006: 37–41. [Google Scholar]
- Grant WS, Bowen BW. 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. [CrossRef] [Google Scholar]
- Hewitt G. 1996. Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc 58: 247–276. [CrossRef] [Google Scholar]
- Hewitt G. 2000. The genetic legacy of the Quaternary ice ages. Nature 405: 907–913. [CrossRef] [PubMed] [Google Scholar]
- Hill CB, Angessa TT, Zhang X, Chen K, Zhou G, Tan C, Wang P, Westcott S, Li C. 2021. A global barley panel revealing genomic signatures of breeding in modern Australian cultivars. Plant J 106: 419–434. [CrossRef] [PubMed] [Google Scholar]
- Jackson JBC, Kirby MX, Berger WH, Bjorndal KA, Botsford LW, Bourque BJ, Bradbury RH, Cooke R, Erlandson J, Estes JA, Hughes TP, Kidwell S, Lange CB, Lenihan HS, Pandolfi JM, Peterson CH, Steneck RS, Tegner MJ, Warner RR. 2001. Historical Overfishing and the Recent Collapse of Coastal Ecosystems. Science 293: 629–637. [CrossRef] [PubMed] [Google Scholar]
- Janssen K, Chavanne H, Berentsen P, Komen H. 2017. Impact of selective breeding on European aquaculture. Aquaculture 472: 8–16. [CrossRef] [Google Scholar]
- Janssen K, Saatkamp H, Komen H. 2018.Cost-benefit analysis of aquaculture breeding programs. Genet. Sel Evol 50: 2. [CrossRef] [Google Scholar]
- Joanna RF, Heather K, Stephen DP. 2015. Molecular ecology. Higher Education Press. [Google Scholar]
- Knowlton N, Weigt LA. 1998. New dates and new rates for divergence across the Isthmus of Panama. Proc Royal Soc B: Biol Sci 265: 2257–2263. [CrossRef] [Google Scholar]
- Kumar S, Stecher G, Tamura K. 2016. MEG A7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol 33: 1870–1874. [CrossRef] [PubMed] [Google Scholar]
- Leigh JW, Bryant D. 2015. Popart: full‐feature software for haplotype network construction. Methods Ecol Evol 6: 1110–1116. [CrossRef] [Google Scholar]
- Li JJ. 2020. Comparative study on two cultivation models of macrobrachium nipponense (Thesis for the Degree of Master). Baoding, Hebei, China: Hebei University. [Google Scholar]
- Li XZ, Liu RY, Liang XQ, Chen GX. 2007. Fauna sinica, Crustacea Decapoda Palaemonoidea, Invertebrata. Beijing: Science Press. [Google Scholar]
- Li Y. 2015. Studies on Morphology and Genetics of Pampus Species (Doctor dissertation). Qingdao, Shandong, China: Ocean University of China. [Google Scholar]
- Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25: 1451–1452. [CrossRef] [PubMed] [Google Scholar]
- Limburg KE, Waldman JR. 2009. Dramatic declines in North Atlantic Diadromous Fishes. BioScience 59: 955–965. [CrossRef] [Google Scholar]
- Long T, Wu L, Meng G, Guo W. 2011. Numerical simulation for impacts of hydrodynamic conditions on algae growth in Chongqing Section of Jialing River, China. Ecol Model 222: 112–119. [CrossRef] [Google Scholar]
- Ma KY, Feng JB, Li JL. 2012. Genetic variation based on microsatellite analysis of the oriental river prawn, Macrobrachium nipponense from Qiandao Lake in China. Genet Mol Res 11: 4235–4244. [Google Scholar]
- Mirzajani A, Ghane A, Bagheri S, Abbasi K, Sayadrahim M, Salahi M, Lavajoo F. 2020. Diet survey and trophic position of Macrobrachium nipponense in the Food Web of Anzali Wetland. Wetlands 40: 1229–1239. [CrossRef] [Google Scholar]
- New T, Xie Z. 2008. Impacts of large dams on riparian vegetation: applying global experience to the case of China’s Three Gorges Dam. Biodivers Conserv 17: 3149–3163. [CrossRef] [Google Scholar]
- Pathirana R, Carimi F. 2022. Management and utilization of plant genetic resources for a sustainable agriculture. Plants 11: 2038. [CrossRef] [PubMed] [Google Scholar]
- Pritchard JK, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155: 945–959. [CrossRef] [PubMed] [Google Scholar]
- Rambaut A, Drummond AJ. Tracer v 1.5. 2007. http://beast.bio.ed.ac.uk/Tracer. [Google Scholar]
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. 2012. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61: 539–542. [Google Scholar]
- Ruzich J, Turnquist K, Nye N, Rowe D, Larson WA. 2019. Isolation by a hydroelectric dam induces minimal impacts on genetic diversity and population structure in six fish species. Conserv Genet 20: 1421–1436. [CrossRef] [Google Scholar]
- Sae-Lim P, Kause A, Mulder HA, Olesen I. 2017. BREEDING AND GENETICS SYMPOSIUM: Climate change and selective breeding in aquaculture. J Anim Sci 95: 1812. [Google Scholar]
- Sewall Wright. 1978. Evolution and the Genetics of Populations. 4: Variability within and among natural populations. The University of Chicago Press. [Google Scholar]
- Suchard MA, Lemey P, Baele G, Ayres DL, Drummond AJ, Rambaut A. 2018. Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10. Virus Evol 4: vey016. [CrossRef] [Google Scholar]
- Wang D, Luo JY, Xiong F, Liu HY, Zhai DD. 2023. Genetic structure of Pseudobrama simoni populations in the mainstream and tributaries of Three Gorges Reservoir. J Hydroecol 44: 124–131. [Google Scholar]
- Wang H, Jiang H, Duan G, Song G, Ling J, Pan T, Hu Y, Zhou H, Yang M. 2021. The genetic diversity of the rice-crayfish eco-farming Procambarus clarkii in Anhui Province, China. Turk J Fish Aquat Sci 22: TRJFAS19904. [Google Scholar]
- Wang JP, Hsu KC, Chiang TY. 2000. Mitochondrial DNA phylogeography of Acrossocheilus paradoxus (Cyprinidae) in Taiwan. Mol Ecol 9: 1483–1494. [CrossRef] [PubMed] [Google Scholar]
- Wang Z, Wu Q, Zhou J, Ye Y, Tong J. 2003. Silver Carp, Hypophthalmichthys molitrix, in the Poyang Lake belong to the Ganjiang River Population Rather than the Changjiang River Population. Environ. Biol Fishes 68: 261–267. [Google Scholar]
- Wei Z. 2022. The complete mitochondrial genomes of five Agrilinae (Coleoptera, Buprestidae) species and phylogenetic implications. ZooKeys 1092: 195–212. [CrossRef] [PubMed] [Google Scholar]
- Xiang R, Wang L, Li H, Tian Z, Zheng B. 2021. Water quality variation in tributaries of the Three Gorges Reservoir from 2000 to 2015. Water Res 195: 116993. [CrossRef] [PubMed] [Google Scholar]
- Xiang R, Wang L, Li H, Tian Z, Zheng B. 2021. Temporal and spatial variation in water quality in the Three Gorges Reservoir from 1998 to 2018. Sci Total Environ 768: 144866. [CrossRef] [PubMed] [Google Scholar]
- Xiong Y, Jiang S, Zhang L, Wang J, Zhang W, Jin S, Gong Y, Wu Y, Qiao H, Fu H. 2023. Genetic diversity and structure assessment of Macrobrachium nipponense populations: implications for the protection and management of genetic resources. Aquat Living Resour 36:5. [CrossRef] [EDP Sciences] [Google Scholar]
- Xu X, Tan Y, Yang G, Li H, Su W. 2011. Impacts of China’s Three Gorges Dam Project on net primary productivity in the reservoir area. Sci Total Environ 409: 4656–4662. [CrossRef] [PubMed] [Google Scholar]
- Yang HF, Yang SL, Xu KH, Milliman JD, Wang H, Yang Z, Chen Z, Zhang CY. 2018. Human impacts on sediment in the Yangtze River: A review and new perspectives. Glob Planet Change 162: 8–17. [CrossRef] [Google Scholar]
- Yang ZB, Liu YY, He RT, Yang K, Ayi QL, Xie SG, Kang YH, Zheng B. 2018. Responses of benthic macroinvertebrates in tributaries of different hydrological regimes in the Three Gorges Reservoir region to reservoir impoundment. Acta Ecol Sin 38: 7231–7241. [Google Scholar]
- Zhou X, Jiang D, Lang X. 2019. A multi-model analysis of ‘Little Ice Age’ climate over China. The Holocene 29: 592–605. [CrossRef] [Google Scholar]
- Zhu KZ. 1973. A preliminary study of climate change in China during the last five thousand years. Meteorol. Sci Technol S1: 1–23. [Google Scholar]
Cite this article as: Wang D, Hu L, Zhang F, Zheng F, Gong M, Xiong F, Liu H, Zhai D.. 2024. Genetic structure of Macrobrachium nipponense, an important farmed freshwater shrimp in China, in the Three Gorges Reservoir. Aquat. Living Resour. 37: 5
All Tables
All Figures
![]() |
Fig. 1 Sampling sites of M. nipponense in the TGR. |
In the text |
![]() |
Fig. 2 Pairwise Fst values (above the diagonal) and genetic distance (GD, below the diagonal) based on mitochondrial COI gene among eight M. nipponense populations in the TGR. All results were tested for significance 1000 times, and the value marked by a red asterisk (*) represents a passing significance test. |
In the text |
![]() |
Fig. 3 STRUCTURE analysis of eight M. nipponense populations in the TGR based on COI gene. |
In the text |
![]() |
Fig. 4 BI phylogenetic tree of M. nipponense haplotypes constructed based on COI gene. |
In the text |
![]() |
Fig. 5 The network diagram of M. nipponense haplotypes constructed based on COI gene. Each circle represents a haplotype, and the size of the circle is proportional to the number of individuals in the population. |
In the text |
![]() |
Fig. 6 Mismatch distribution analyse for each M. nipponense population from different localities. |
In the text |
![]() |
Fig. 7 Bayesian skyline plot of the effective population sizes through time for M. nipponense population. The X axis is the time scale in units of thousands of years, and the Y axis is the estimated effective population size. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.