Genetic diversity and phylogeography of Daphnia similoides sinensis located in the middle and lower reaches of the Yangtze River

Geographical patterns, climate, and environmental change have important influences on the distribution and spread of aquatic organisms. However, the relationships between the geographical pattern and phylogenetics of Daphnia as well as environmental change are not well known. The genetic diversity and phylogeography of seven D. similoides sinensis populations located in the middle and lower reaches of the Yangtze River were investigated based on the combination of mitochondrial (COI gene) and nuclear (14 microsatellite primers) markers. Based on the mitochondrial gene markers, D. similoides sinensis from the middle and lower reaches of the Yangtze River had one ancestral haplotype and two evolutionary clades. In addition, D. similoides sinensis population deviated from neutral evolution, showing signs of a bottleneck effect followed by population expansion. Based on the microsatellite markers, the seven D. similoides sinensis populations formed three main groups. The dendrogram (NJ/ME) showed that D. similoides sinensis based on the mitochondrial genes marker were obviously clustered two main clades, whereas there were three clades based on the microsatellite markers. Our results suggested that the habitat fragmentation due to the barrier of the dams and sluices promoted the genetic differentiation and phylogeography of D. similoides sinensis populations in the middle and lower reaches of the Yangtze River.     

Population genetics and sympatric divergence of the freshwater gudgeon,Gobiobotia filifer,in the Yangtze River inferred from mitochondrial DNA

The ecosystem and Pleistocene glaciations play important roles in population demography. The freshwater gudgeon, Gobiobotia filifer, is an endemic benthic fish in the Yangtze River and is a good model for ecological and evolutionary studies. This study aimed to decode the population structure of G. filifer in the Yangtze River and reveal whether divergence occurred before or after population radiation. A total of 292 specimens from eight locations in the upper and middle reaches of the Yangtze River were collected from 2014 to 2016 and analyzed via mitochondrial DNA Cyt b gene sequencing. A moderately high level of genetic diversity was found without structures among the population. However, phylogenetic and network topology showed two distinct haplotype groups, and each group contained a similar proportion of individuals from all sampled sites. This suggested the existence of two genetically divergent source populations in G. filifer. We deduced that a secondary contact of distinct glacial refugia was the main factor creating sympatric populations of G. filifer, and climate improvement promoted population expansion and colonization.     

Genetic variation and bidirectional gene flow in the riparian plant Miscanthus lutarioriparius,across its endemic range: implications for adaptive potential

Miscanthus lutarioriparius is an endemic species that grows along the middle and lower reaches of the Yangtze River and is a valuable source of germplasm for the development of second‐generation energy crops. The plant that propagates via seeds, stem nodes, and rhizomes shows high phenotypic variation and strong local adaptation. Here, we examined the magnitude and spatial distribution of genetic variation in M. lutarioriparius across its entire distributional range and tested underlying factors that shaped its genetic variation. Population genetic analyses were conducted on 644 individuals from 25 populations using 16 microsatellite markers. M. lutarioriparius exhibited a high level of genetic variation (H_E = 0.682–0.786; A_r = 4.74–8.06) and a low differentiation (F_ST = 0.063; D_est = 0.153). Of the total genetic variation, 10% was attributed to the differences among populations (df = 24, P < 0.0001), whereas 90% was attributed to the differences among individuals (df = 619, P ≤ 0.0001). Genetic diversity did not differ significantly across longitudes and did not increase in the populations growing downstream of the Yangtze River. However, significant associations were found between genetic differentiation and spatial distance. Six genetic discontinuities were identified, which mostly distributed among downstream populations. We conclude that anthropogenic factors and landscape features both contributed to shaping the pattern of gene flow in M. lutarioriparius, including long‐distance bidirectional dispersal. Our results explain the genetic basis of the high degree of adaptability in M. lutarioriparius and identify potential sources of new germplasm for the domestication of this potential second‐generation energy crop.     

Population genetic structure and colonization history of short ninespine sticklebacks (Pungitius kaibarae)

The contemporary distribution and genetic structure of a freshwater fish provide insight into its historical geodispersal and geographical isolation following Quaternary climate changes. The short ninespine stickleback, Pungitius kaibarae, is a small gasterosteid fish occurring in freshwater systems on the Korean Peninsula and in southeast Russia. On the Korean Peninsula, P. kaibarae populations are distributed in three geographically separated regions: the NE (northeast coast), SE (southeast coast), and a limited area in the ND (Nakdong River). In this study, we used mitochondrial loci and microsatellites to investigate the evolutionary history of P. kaibarae populations by assessing their pattern of genetic structure. Our analyses revealed a marked level of divergence among three regional populations, suggesting a long history of isolation following colonization, although ND individuals showed relatively higher genetic affinity to populations from SE than those from NE. The populations from NE showed a great degree of interpopulation differentiation, whereas populations from SE exhibited only weak genetic structuring. Upon robust phylogenetic analysis, P. kaibarae formed a monophyletic group with Russian P. sinensis and P. tymensis with strong node confidence values, indicating that P. kaibarae populations on the Korean Peninsula originated from the southward migration of its ancestral lineage around the middle Pleistocene.     

Effect of Yangtze River on population genetic structure of the relict plant Parrotia subaequalis in eastern China

Parrotia subaequalis (Hamamelidaceae) is a Tertiary relic species endemic in eastern China. We used inter‐simple sequence repeat (ISSR) markers to access genetic diversity and population genetic structure in natural five populations of P. subaequalis. The levels of genetic diversity were higher at species level (H = 0.2031) but lower at population level (H = 0.1096). The higher genetic diversity at species levels might be attributed to the accumulation of distinctive genotypes which adapted to the different habitats after Quaternary glaciations. Meanwhile, founder effects on the early stage, and subsequent bottleneck of population regeneration due to its biological characteristics, environmental features, and human activities, seemed to explain the low population levels of genetic diversity. The hierarchical AMOVA revealed high levels (42.60%) of among‐population genetic differentiation, which was in congruence with the high levels of Nei's genetic differentiation index (G_ST = 0.4629) and limited gene flow (N_m = 0.5801) among the studied populations. Mantel test showed a significant isolation‐by‐distance, indicating that geographic isolation has a significant effect on genetic structure in this species. Unweighted pair‐group method with arithmetic average clustering, PCoA, and Bayesian analyses uniformly recovered groups that matched the geographical distribution of this species. In particular, our results suggest that Yangtze River has served as a natural barrier to gene flow between populations occurred on both riversides. Concerning the management of P. subaequalis, the high genetic differentiation among populations indicates that preserving all five natural populations in situ and collecting enough individuals from these populations for ex situ conservation are necessary.     

Salinity appears to be the main factor shaping spatial COI diversity of Corbicula lineages within the Chinese Yangtze River Basin

Aim

To date, few studies have examined the phylogenetics of Corbicula clams in their native range and the environmental parameters influencing their distribution, although this could provide great insights into the biological adaptation and invasion dynamics of Corbicula clams. We sought to identify the genetic lineages of native Corbicula clams and elucidate the environmental factors shaping the distributions of identified lineages.

Location

China, mainly the Yangtze River Basin.

Methods

The alignment comprised 558 COI sequences including samples from China and 222 COI sequences from published studies. This dataset was used to generate phylogenetic trees and compare population diversity. We used dbRDA method to assess the relationship between these COI data and environmental factors measured to identify the important factors affecting Corbicula's distribution.

Results

The COI phylogenetic tree delineated the monophyly of 3 major COI clades and 77 distinct COI haplotypes in the Yangtze River Basin. The invasive lineage C/S (FW17) was not identified in our sampling in China, while invasive lineages A/R (FW5 = Hap6), B (FW1 = Hap17) and Rlc (FW4 = Hap43) were identified, abundant and widely distributed in the middle and lower reaches of Yangtze River. Focusing on populations from Yangtze River Basin, both the COI haplotype and nucleotide diversity in the lakes along the Yangtze River increased with longitude, except for the river mouth population. The AMOVA tests showed significant differentiation between the middle and lower reaches of the Yangtze River and among populations. The dbRDA results suggested that the parameter chloride explained most of the spatial COI haplotype distribution variation in the Yangtze River Basin, with the three invasive lineages tolerating broad fluctuations of salinity (chloride levels ranging from 4 to 60 mg/L).

Conclusions

The Corbicula COI haplotypes found within the Yangtze River Basin had distinct distribution preferences, with the invasive androgenetic lineages being the most abundant and widely distributed. Genetic diversity was higher in this native region than in invaded areas in Europe and America, while it decreased with increasing distance from the river mouth. Salinity appeared to be the main environmental factor shaping the COI haplotype distribution of Corbicula lineages within their native range.     

Genetic diversity and population structure of a cyprinid fish (Ancherythroculter nigrocauda) in a highly fragmented river

In this study, samples of Ancherythroculter nigrocauda, an endemic fish in the upper Yangtze River, were collected above and below dams in the Longxi River, a tributary of the upper Yangtze River, China, to investigate the genetic impacts of dams. The mitochondrial cytochrome b gene (cyt b) and 13 microsatellite (SSR) loci were used to analyze whether dams have resulted in loss of genetic diversity of the two fragmented populations or caused genetic differentiation between them. The results showed that the haplotype diversity (0.488; 0.486), nucleotide diversity (0.084%; 0.082%) and average expected heterozygosity (0.652; 0.676) of the two populations were all at a low level, and recent bottlenecks were detected. However, there was no genetic differentiation detected by the low genetic differentiation index (F_st, cyt b: ?0.1677, p = 0.99707; SSR: 0.00259, p = 0.81427). Besides, 11 pairs of half‐sibling relationship were found between the two populations indicating that there were individual movements and gene flow between them. This could be the larvae moving from upstream to downstream when water spilled over dams in flooding season. Therefore, our analysis showed that the dams have caused a loss of genetic diversity of the populations of A. nigrocauda in the Longxi River, blocked the active upstream movement but allowing passive downstream drift of larvae.     

Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Diverse processes shape deep phylogeographical divergence in Cobitis sinensis (Teleostei: Cobitidae) in East Asia

Little has been known about the impacts of past vicariance events on the phylogeography and population structure of freshwater fishes in East Asia. The aims of this study are to assess the genetic variability with extensive sampling throughout the range of Chinese spiny loach, Cobitis sinensis, and to infer the genetic structure and evolutionary history of populations. Cobitis sinensis in China may have initiated from two ancestral populations, namely Yangtze and Pearl Rivers, which diverged about 7.24 MYA likely due to drainage systems alteration. In the phylogroup I, a southward dispersal event occurred from East China (Yangtze River) to south ZheMin and Hainan subregions, followed by eastward dispersal from ZheMin to south Taiwan. In the phylogroup II, eastward colonization took place from Pearl River to north Taiwan in the late Pliocene, coupled with loss of genetic diversity in the island populations. This study showed that Cenozoic tectonic movements and climatic and sea‐level fluctuations may have shaped the genetic structure of C. sinensis in concert. Highly diverged mtDNA sequences suggest existence of cryptic species in morphospecies C. sinensis.     

Genetic diversity of common carp from two largest Chinese lakes and the Yangtze River revealed by microsatellite markers

Six polymorphic microsatellites (eight loci) were used to study the genetic diversity and population structure of common carp from Dongting Lake (DTC), Poyang Lake (PYC), and the Yangtze River (YZC) in China. The gene diversity was high among populations with values close to 1. The number of alleles per locus ranged from 2 to 11, and the average number of alleles among 3 populations ranged from 6.5 to 7.9. The mean observed (H _O) and expected (H _E) heterozygosity ranged from 0.4888 to 0.5162 and from 0.7679 to 0.7708, respectively. Significant deviations from Hardy–Weinberg Equilibrium expectation were found at majority of the loci and in all three populations in which heterozygote deficits were apparent. The analysis of molecular variance (AMOVA) indicated that the percent of variance among populations and within populations were 3.03 and 96.97, respectively. The Fst values between populations indicated that there were significant genetic differentiations for the common carp populations from the Yangtze River and two largest Chinese freshwater lakes. The factors that may result in genetic divergence and significant reduction of the observed heterozygosity were discussed.     

Population genetic structure and adaptive differentiation of iron walnut Juglans regia subsp. sigillata in southwestern China

Southwestern (SW) China is an area of active tectonism and erosion, yielding a dynamic, deeply eroded landscape that influences the genetic structure of the resident populations of plants and animals. Iron walnut (Juglans regia subsp. sigillata) is a deciduous tree species endemic to this region of China and cultivated there for its edible nuts. We sampled 36 iron walnut populations from locations throughout the species' range in SW China and genotyped a total of 765 individuals at five chloroplast DNA regions and 22 nuclear microsatellite loci. Species distribution models were produced to predict the evolution and historical biogeography of iron walnut and to estimate the impacts of climate oscillations and orographic environments on the species' demography. Our results indicated that J. regia subsp. sigillata had relatively low genetic diversity, high interpopulation genetic differentiation, and asymmetric interpopulation gene flow. Based on DIYABC analysis, we identified two lineages of J. sigillata in southwestern China. The lineages (subpopulations) diverge during the last glacial period (~1.34 Ma). Southwestern China was a glacial refuge during the last glacial period, but increasingly colder and arid climates might have fostered the fragmentation of J. regia subsp. sigillata within this refugium. Finally, we found that recent habitat fragmentation has led to a reduction in population connectivity and increased genetic differentiation by genetic drift in isolated populations. Our results support a conclusion that geological and climatic factors since the Miocene triggered the differentiation, evolutionary origin, and range shifts of J. sigillata in the studied region.     

Microsatellite diversity and population genetic structure of redfin culter ( Culter erythropterus ) in fragmented lakes of the Yangtze River

Redfin culter (Culter erythropterus) is a small lethic freshwater fish and widely distributed in the adjacent lakes of the Yangtze River of China. Five microsatellite loci were applied to investigate the genetic variation and population structure of redfin culter from seven lakes in the middle-and-lower reaches of the Yangtze River. The gene diversity was high among the populations (H > 0.9), the average number of alleles among seven populations was low with a range from 2.00 to 3.87. The mean observed (H _O) and expected (H _E) heterozygosity ranged from 0.111 to 0.419 and from 0.162 to 0.750, respectively. Significant deviations from Hardy-Weinberg Equilibrium expectation were found in 50% of the total locus-population combination tests in which heterozygote deficits were apparent. The analysis of molecular variance (AMOVA) indicated that the percentage of variance among and within these populations were 6.18 and 93.82, respectively. The Fst values (0.062, P < 0.001) among studied populations indicated that there were significant genetic differentiations among redfin culture populations from the scattered lakes with different connections to the Yangtze River. These results are useful for the evaluation and conservation of small freshwater fishes. The factors that may be involved in low intra-population polymorphism and the pattern of the population genetic structure of redfin culter from the Yangtze River were discussed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Handling editor: C. Sturmbauer     

中国栽培稻遗传多样性中心和起源研究

本研究对中国栽培稻6个地理分布群的700份古老地方栽培品种进行9个多态性等位酶基因位点的遗传多样性分析.结果表明籼稻和粳稻的平均基因多样性均以云南最大,淮河上游次之,黄河以北最小.中国栽培稻有3个遗传多样性中心:云南,长江中游-淮河上游,华南.长江中游-淮河上游可能是中国栽培稻的起源中心.结合考古学资料及前人的工作,认为云南不是中国栽培稻的起源中心,而可能是中国起源中心衍生的一个次生中心并受到南亚中心的强烈影响.     

Sequence Polymorphism and Geographical Variation at a Positively Selected MHC-DRB Gene in the Finless Porpoise (Neophocaena phocaenoides): Implication for Recent Differentiation of the Yangtze Finless Porpoise?

Sequence polymorphism at the MHC class II DRB locus was investigated in three finless porpoise (Neophocaena phocaenoides) populations in Chinese waters. Intragenic recombination and strong positive selection were the main forces in generating sequence diversity in the DRB gene. MHC sequence diversity changed significantly along the study period. Significant decrease in heterozygosity and lost alleles have been detected in the Yangtze River population and South China Sea population since 1990. Furthermore, there is a trend of increasing population differentiation over time. Especially, the genetic differentiation between the Yangtze River population and the Yellow Sea population was very low prior to 1990 (F _ST = 0.036, P = 0.009), but became very significant after 1990 (F _ST = 0.134, P < 0.001), suggesting a recent augmentation of genetic differentiation between both populations probably in a relatively short-term period. Porpoises from the Yangtze River displayed divergent frequencies of shared and private alleles from those displayed by two marine populations, which suggest that the former riverine population has been under a different selection regime (characteristic of a fresh water environment) than that of its marine counterparts.     

Genetic architecture of a body colour cline in Drosophila americana

Phenotypic variation within a species is often structured geographically in clines. In Drosophila americana, a longitudinal cline for body colour exists within North America that appears to be due to local adaptation. The tan and ebony genes have been hypothesized to contribute to this cline, with alleles of both genes that lighten body colour found in D. americana. These alleles are similar in sequence and function to the allele fixed in D. americana's more lightly pigmented sister species, Drosophila novamexicana. Here, we examine the frequency and geographic distribution of these D. novamexicana‐like alleles in D. americana. Among alleles from over 100 strains of D. americana isolated from 21 geographic locations, we failed to identify additional alleles of tan or ebony with as much sequence similarity to D. novamexicana as the D. novamexicana‐like alleles previously described. However, using genetic analysis of 51 D. americana strains derived from 20 geographic locations, we identified one new allele of ebony and one new allele of tan segregating in D. americana that are functionally equivalent to the D. novamexicana allele. An additional 5 alleles of tan also showed marginal evidence of functional similarity. Given the rarity of these alleles, however, we conclude that they are unlikely to be driving the pigmentation cline. Indeed, phenotypic distributions of the 51 backcross populations analysed indicate a more complex genetic architecture, with diversity in the number and effects of loci altering pigmentation observed both within and among populations of D. americana. This genetic heterogeneity poses a challenge to association studies and genomic scans for clinal variation, but might be common in natural populations.     

Microsatellite diversity and population genetic structure of yellowcheek, Elopichthys bambusa (Cyprinidae) in the Yangtze River

Yellowcheek carp (Elopichthys bambusa) is the only species of genus Elopichthys. It is widely distributed in Chinese freshwaters but currently its populations have declined to threatening level. We examined the genetic diversity and population structure of E. bambusa in the Yangtze River basin. A total of nine polymorphic microsatellite markers were employed to study five populations occurring in middle and lower reaches of the river. The results revealed low-to-moderate genetic diversity. The number of alleles per locus varied between 3 and 8 with an average of 4.8. Observed heterozygosity ranged from 0.15 to a maximum of 1.00. Significant deviations (P < 0.01) from Hardy–Weinberg equilibrium were observed for all the tested locus-population combinations with clear heterozygosity deficits. AMOVA indicated that majority of the variance lies within populations (93.81%) than among the populations (7.05%). Pairwise F_ST and unbiased genetic distance pointed out significant differentiation among the samples from populations with different connections to the Yangtze River. In the UPGMA dendrogram, clustering pattern of populations indicated that most of the populations are reproductively isolated due to anthropogenic interventions. Clustering of PYL and DTL populations shows ongoing gene flow through the mainstream. The recent hydrological alterations and overfishing are major factors shaping the current genetic structure. These results can be helpful for effective management and sustainable conservation of E. bambusa populations.     

Microsatellite‐based analysis of genetic structure and gene flow of Mythimna separata (Walker) (Lepidoptera: Noctuidae) in China

The oriental armyworm, Mythimna separata, is a serious agricultural pest in China. Seasonal and roundtrip migration has recently led to sudden, localized outbreaks and crop losses. To evaluate genetic differentiation between populations in eastern and western China and elucidate gene flow, the genetic structure of 20 natural populations from nine provinces was examined using seven microsatellite markers. The results indicated high genetic diversity. However, little to moderate (0 < F_ST < 0.15) genetic differentiation was detected, and there was no correlation between genetic distance and geographical distance. Bayesian clustering analysis identified three groups whereas discriminant analysis of principal components identified ten clusters that were considered as two clear‐cut clusters and one admixed group. Gene flow occurred frequently in most population pairs, and an asymmetrical migration rate was detected in several pairwise population comparisons. The bottleneck test showed that few populations had experienced recent bottlenecks. Correspondingly, large‐scale and long‐distance migration of M. separata has caused low genetic differentiation and frequent gene exchange. Our findings are important for studying genetic evolution and help to improve predictions of M. separata outbreaks in China.     

Mitochondrial introgressive hybridization following a demographic expansion in the tomato frogs of Madagascar,genus Dyscophus

Madagascar is a biodiversity hotspot with a unique fauna and flora largely endemic at the species level and highly threatened by habitat destruction. The processes underlying population‐level differentiation in Madagascar's biota are poorly understood and have been proposed to be related to Pleistocene climatic cycles, yet the levels of genetic divergence observed are often suggestive of ancient events. We combined molecular markers of different variability to assess the phylogeography of Madagascar's emblematic tomato frogs (Dyscophus guineti and D. antongilii) and interpret the observed pattern as resulting from ancient and recent processes. Our results suggest that the initial divergence between these taxa is probably old as reflected by protein‐coding nuclear genes and by a strong mitochondrial differentiation of the southernmost population. Dramatic changes in their demography appear to have been triggered by the end of the last glacial period and possibly by the short return of glacial conditions known as the 8K event. This dramatic change resulted in an approximately 50‐fold reduction of the effective population size in various populations of both species. We hypothesize these species' current mitochondrial DNA diversity distribution reflects a swamping of the mitochondrial genetic diversity of D. guineti by that of D. antongilii previous to the populations' bottlenecks during the Holocene, and probably as a consequence of D. antongilii demographic expansion approximately 1 million years ago. Our data support the continued recognition of D. antongilii and D. guineti as separate species and flag D. guineti as the more vulnerable species to past and probably also future environmental changes.     

Low population differentiation revealed in the highly threatened elongate loach (Leptobotia elongata,Bleeker), a species endemic to the fragmented upper reaches of the Yangtze River

The elongate loach (Leptobotia elongata) is a typical migratory fish that was formerly found in abundance throughout the upper and middle reaches of the Yangtze River but is now classified as a highly threatened fish on the China Species Red List. Here, a comparative study of microsatellite loci was performed to evaluate the genetic diversity and analyze the population structure in seven locations in the upper reaches of the Yangtze River, China. The results showed that genetic diversity of L. elongata is generally still at a high level. Populations from Xinshi, Panzhihua, Hechuan and Chishui, were found to have lower levels of genetic diversity than those that lived in Leshan, Yibin and Jiangjin. The IAM model indicated that populations from Panzhihua, Xinshi and Hechuan have undergone recent bottlenecks. AMOVA analysis revealed a small amount of differentiation among populations (3.6%), but most of the total variation occurred within populations (96.4%). Pairwise comparisons of populations confirmed a low population differentiation. STRUCTURE analysis showed no obvious genetic structure among populations, indicating that L. elongata living in the main stem of the river and its tributaries should be considered a single management unit.     

Genetic structure and historical demography of Schizothorax nukiangensis (Cyprinidae) in continuous habitat

Geographic distance, different living habitats or Pleistocene climatic oscillations have frequently been found to shape population genetic structure in many species. The genetic structure of Schizothorax nukiangensis, a high altitude, valuable fish species, which is distributed throughout the Nujiang River, was investigated by mitochondrial DNA sequence analysis. The cytochrome c oxidase subunit I (COI), cytochrome b (cytb), and the mitochondrial control region (MCR) of S. nukiangensis were concatenated for examination of population structure and demographic history. The concatenated data set (2405 bp) implied a pronounced genetic population structure (overall F_ST = 0.149) and defined two population units. Strong differentiation was detected between the Sanjiangkou (SJK) population and other populations due to environmental heterogeneity, dispersal ability, and/or glacial cycles. Additional DNA sequencing of the nuclear RAG2 gene also examined significant differentiation between two units and between SJK and the upstream populations (U‐unit). Recent expansion events suggest that S. nukiangensis may have undergone a rapid increase during warm interglacial periods. Surprisingly, S. nukiangensis appears to have undergone an obvious expansion during the last glaciations (LG) for cold hardiness and a sharp contraction from 1.5 ka to the present. However, two population units exhibited different reflections during the LG, which might be closely related to their living habitats and cold hardiness. A clear pattern of isolation by distance was detected in S. nukiangensis due to feeding habits, limited dispersal ability, and/or philopatry. It is vitally important that more attention be given to S. nukiangensis due to low genetic diversity, lack of gene flow, and recent population contraction.