Effects of Large-Scale Releases on the Genetic Structure of Red Sea Bream (Pagrus major,Temminck et Schlegel) Populations in Japan

Large-scale hatchery releases are carried out for many marine fish species worldwide; nevertheless, the long-term effects of this practice on the genetic structure of natural populations remains unclear. The lack of knowledge is especially evident when independent stock enhancement programs are conducted simultaneously on the same species at different geographical locations, as occurs with red sea bream (Pagrus major, Temminck et Schlegel) in Japan. In this study, we examined the putative effects of intensive offspring releases on the genetic structure of red sea bream populations along the Japanese archipelago by genotyping 848 fish at fifteen microsatellite loci. Our results suggests weak but consistent patterns of genetic divergence (F_ST = 0.002, p < 0.001). Red sea bream in Japan appeared spatially structured with several patches of distinct allelic composition, which corresponded to areas receiving an important influx of fish of hatchery origin, either released intentionally or from unintentional escapees from aquaculture operations. In addition to impacts upon local populations inhabiting semi-enclosed embayments, large-scale releases (either intentionally or from unintentional escapes) appeared also to have perturbed genetic structure in open areas. Hence, results of the present study suggest that independent large-scale marine stock enhancement programs conducted simultaneously on one species at different geographical locations may compromise native genetic structure and lead to patchy patterns in population genetic structure.     

Spatial and temporal genetic structure in a hybrid cordgrass invasion

Invasive hybrids and their spread dynamics pose unique opportunities to study evolutionary processes. Invasive hybrids of native Spartina foliosa and introduced S. alterniflora have expanded throughout San Francisco Bay intertidal habitats within the past 35 years by deliberate plantation and seeds floating on the tide. Our goals were to assess spatial and temporal scales of genetic structure in Spartina hybrid populations within the context of colonization history. We genotyped adult and seedling Spartina using 17 microsatellite loci and mapped their locations in three populations. All sampled seedlings were hybrids. Bayesian ordination analysis distinguished hybrid populations from parent species, clearly separated the population that originated by plantation from populations that originated naturally by seed and aligned most seedlings within each population. Population genetic structure estimated by analysis of molecular variance was substantial (F_ST=0.21). Temporal genetic structure among age classes varied highly between populations. At one population, the divergence between adults and 2004 seedlings was low (F_ST=0.02) whereas at another population this divergence was high (F_ST=0.26). This latter result was consistent with local recruitment of self-fertilized seed produced by only a few parental plants. We found fine-scale spatial genetic structure at distances less than ∼200 m, further supporting local seed and/or pollen dispersal. We posit a few self-fertile plants dominating local recruitment created substantial spatial genetic structure despite initial long-distance, human dispersal of hybrid Spartina through San Francisco Bay. Fine-scale genetic structure may more strongly develop when local recruits are dominated by the offspring of a few self-fertile plants.     

Coalescence Times and FST Under a Skewed Offspring Distribution Among Individuals in a Population

Estimates of gene flow between subpopulations based on F_ST (or N_ST) are shown to be confounded by the reproduction parameters of a model of skewed offspring distribution. Genetic evidence of population subdivision can be observed even when gene flow is very high, if the offspring distribution is skewed. A skewed offspring distribution arises when individuals can have very many offspring with some probability. This leads to high probability of identity by descent within subpopulations and results in genetic heterogeneity between subpopulations even when Nm is very large. Thus, we consider a limiting model in which the rates of coalescence and migration can be much higher than for a Wright–Fisher population. We derive the densities of pairwise coalescence times and expressions for F_ST and other statistics under both the finite island model and a many-demes limit model. The results can explain the observed genetic heterogeneity among subpopulations of certain marine organisms despite substantial gene flow.     

Characterizing genetic integrity of rear-edge trout populations in the southern Appalachians

Vertebrate populations at the periphery of their range can show pronounced genetic drift and isolation, and therefore offer unique challenges for conservation and management. These populations are often candidates for management actions such as translocations that are designed to improve demographic and genetic integrity. This is particularly true of coldwater species like brook trout (Salvelinus fontinalis), whose numbers have declined greatly across its historic range. At the southern margin, remnant wild populations persist in isolated headwater streams, and many have a history of receiving translocated individuals through either stocking of hatchery reared fish, relocation of wild fish, or both during restoration attempts. To determine current genetic integrity and resolve the genetic effects of past management actions for brook trout populations in SC, USA, we genetically assessed all 18 documented remaining brook trout populations along with individuals acquired from six hatcheries with recorded stocking events in SC. Our results indicated that six of the 18 streams showed signs of hatchery admixture (range 57–97%) and restored patches retained genetic signatures from multiple source populations. Populations had among the lowest genetic diversity (min average H_E?=?0.147) and effective number of breeders (mean N_b?=?31.2) estimates observed throughout the native brook trout range. Populations were highly differentiated (mean pair-wise F_ST?=?0.396), and substantial genetic divergence was evident across major river drainages (max pair-wise F_ST?=?0.773). The lowest local genetic diversity and highest genetic differentiation ever reported for this species make its conservation a challenging task, particularly when combined with other threats such as climate change and non-native species. We offer recommendations on managing peripheral populations with depleted genetic characteristics and provide a reference for determining which existing populations will best serve as sources for future translocation efforts aimed at enhancing or restoring wild brook trout genetic integrity.     

Little to no genetic structure in the ectomycorrhizal basidiomycete Suillus spraguei (Syn. S. pictus) across parts of the northeastern USA

Population genetic studies of ectomycorrhizal fungi at scales larger than 100 km are still relatively rare with highly variable results. In this study, we determined the population genetic structure of the ectomycorrhizal basidiomycete Suillus spraguei over distances up to 600 km in northeastern USA forests. S. spraguei associates exclusively with five-needled pines and only with white pine (Pinus strobus) in the eastern USA. We used six microsatellite loci to assess the genetic structure between eight sites sampled in the Adirondack Park of New York and seven sites sampled in other forests of New York, Pennsylvania, and Massachusetts. Except for one site, little to no genetic differentiation was detected in pairwise comparisons of the sites (F _ST?=?0 to 0.05). Only one site was moderately differentiated from most other sites (F _ST?=?0.02 to 0.15). The Mantel test showed no significant correlation between genetic and geographic distances (isolation by distance; R ²?=?0.003, P?=?0.3). The STRUCTURE analysis also supported the presence of a single cluster (K?=?1).     

Is the genetic diversity of small scattered forest tree populations at the southern limits of their range more prone to stochastic events? A wild cherry case study by microsatellite-based markers

Wild cherry (Prunus avium L.) is a widespread, partially asexual, noble hardwood European species characterized by a scattered distribution, small population sizes, and human exploitation for its valuable wood. These characteristics, especially at the southern limits of the species natural distribution where additional varying stresses may occur, render P. avium populations prone to potential stochastic, genetic, and demographic events. In this study, we used dominant inter simple sequence repeat (ISSR) and codominant simple sequence repeat (SSR) markers to infer the genetic structure of P. avium. Five populations from northern Greece were evaluated based on 46 ISSR and 11 SSR loci. Populations presented a relatively high level of genetic variation, with a mean genetic diversity of H _e?=?0.166 and H _e?=?0.740 regarding ISSR and SSR analysis, respectively. We observed moderate population differentiation for ISSR (G _ST?=?0.113) and SSR (F _ST?=?0.097) markers. AMOVA also detected significant differentiation among populations for ISSRs (?? _ST?=?0.338) and SRRs (?? _ST?=?0.162). According to linkage disequilibrium analysis, estimates of effective population size were generally sufficient for maintaining extant genetic variability and evolutionary potential. A possible bottleneck was detected for only one population. In general, it appears that despite the particular characteristics of the P. avium populations studied, genetic stochasticity events were not apparent. The studied populations, located at the rear edge of the species European distribution, reveal a wealth of genetic variation that is very valuable for the genetic conservation of local adaptive gene complexes, especially under contemporary climatic change scenarios.     

Analysis of Genetic Variation in Natural Populations of Medicago truncatula of Southern Tunisian Ecological Areas, Using Morphological Traits and SSR Markers

We used 19 quantitative traits and 14 microsatellite markers (SSRs) to analyze the genetic variation in four natural populations of the model legume Medicago truncatula sampled in southern Tunisia. The greatest genetic variation of quantitative traits and molecular markers occurred within populations (>71%). In contrast to quantitative population differentiation (Q _ST ?=?0.09), a high level of molecular differentiation (F _ST ?=?0.23) was found among populations. The majority of quantitative traits exhibited Q _ST values significantly less than F _ST values, suggesting that selection may be acting to suppress differentiation for these traits. There was no significant correlation between genetic variation of quantitative traits and molecular markers within populations. On the other hand, significant correlations were found between measured quantitative characters and the site-of-origin environmental factors. The eco-geographical factors with the greatest influence on the variation of measured traits among populations were altitude, followed by soil texture, assimilated phosphorus (P₂O₅) and organic matter. Nevertheless, there were no consistent patterns of associations between gene diversity (He) and eco-geographical factors.     

Patterns of Genetic Diversity in Platanthera bifolia (Orchidaceae) with Respect to Life History Traits and Recent Range Expansion

For evolutionary and ecological analyses, genetic diversity at different scales needs to be studied in terms of biological properties, habitat, population size and population history. We surveyed Platanthera bifolia populations from six regions in northeastern Poland to determine the impact of the mating system and population history on genetic diversity. Based on variation at allozyme markers, genetic variation was relatively moderate (P?=?22.3%, A?=?1.48, H _O?=?0.083, F _IS?=??0.015) and similar to other Platanthera species. These parameters varied between populations (P?=?13.3%–26.6%, A?=?1.26–1.66, H _O?=?0.055–0.111, F _IS?=??0.262–0.147). The genetic diversity patterns were shaped by different proportions of facilitated selfing and/or outcrossing, resulting in positive and negative F _IS values, respectively. No relationship between inbreeding coefficient and population size, however, and no impact of apomixis on the level of genetic diversity of P. bifolia were found. The relatively low level of genetic differentiation among the investigated regions (F _CT?=?0.002, P?>?0.05) and among populations (F _ST?=?0.048, P?<?0.001), and the lack of a significant relationship between genetic and geographical distance, are discussed in the context of possible scenaria of postglacial expansion.     

Geographic Influence on Genetic Structure in the Widespread Neotropical Tree Simarouba amara (Simaroubaceae)

This study aimed to assess the population genetic structure of a widespread Neotropical tree species, Simarouba amara, at local, regional and continental spatial scales. We used five microsatellite loci to examine genetic variation in 14 natural populations (N?=?478 individuals) of this vertebrate-dispersed rain forest tree species in Panama, Ecuador, and French Guiana. Estimates of genetic differentiation (F _st and R _st) were significant among all but one population pair and global differentiation was moderate (F _st?=?0.25, R _st?=?0.33) with 94% of genetic variation ascribed to differences among three main geographic regions (Central America, Western Ecuador, Amazon basin). There was no evidence of isolation by distance within regions. Allele-size mutations contributed significantly (R _ST > F _ST) to the divergences between cis- and trans-Andean populations, highlighting the role of the northern Andean cordilleras as an important geographic barrier for this species.     

Spatial Scales of Pollen and Seed-Mediated Gene Flow in Tropical Rain Forest Trees

Gene flow via seed and pollen is a primary determinant of genetic and species diversity in plant communities at different spatial scales. This paper reviews studies of gene flow and population genetic structure in tropical rain forest trees and places them in ecological and biogeographic context. Although much pollination is among nearest neighbors, an increasing number of genetic studies report pollination ranging from 0.5–14 km for canopy tree species, resulting in extensive breeding areas in disturbed and undisturbed rain forest. Direct genetic measures of seed dispersal are still rare; however, studies of fine scale spatial genetic structure (SGS) indicate that the bulk of effective seed dispersal occurs at local scales, and we found no difference in SGS (Sp statistic) between temperate (N?=?24 species) and tropical forest trees (N?=?15). Our analysis did find significantly higher genetic differentiation in tropical trees (F _ST?=?0.177; N?=?42) than in temperate forest trees (F _ST?=?0.116; N?=?82). This may be due to the fact that tropical trees experience low but significant rates of self-fertilization and bi-parental inbreeding, whereas half of the temperate tree species in our survey are wind pollinated and are more strictly allogamous. Genetic drift may also be more pronounced in tropical trees due to the low population densities of most species.     

Genetic diversity and population differentiation in the yellow drum Nibea albiflora along the coast of the China Sea

Nibea albiflora (yellow drum) is an important seafood fish species in East Asia. We explored the population genetic variation of N. albiflora along the coastal waters of the China Sea using microsatellite markers to facilitate a selective breeding programme that is undertaken in China. A total of 256 alleles were detected at 12 loci in four wild populations. A high level of genetic diversity was observed with the mean number of alleles and the observed and expected heterozygosity in each population ranging from 7.917 to 14.083, 0.701 to 0.764 and 0.765 to 0.841, respectively. Pairwise fixation index (F_ST) analysis indicated significant but weak genetic differentiation among populations from four localities (F_ST?=?0.030, P?<?0.01), which was also confirmed by analysis of molecular variance (AMOVA). Significant genetic differentiation was detected between Ningde and the other populations (F_ST?=?0.047–0.056, P?<?0.01). Structure analysis suggested that N. albiflora within the examined range might be composed of two stocks. The data of the present study revealed high genetic diversity and low genetic differentiation among the N. albiflora populations along the coast of the China Sea. This baseline information could be valuable for future selective breeding programmes of N. albiflora.     

Genetic diversity and relationships among wild and cultivated Stenocereus queretaroensis populations in western Mexico

Stenocereus queretaroensis is an endemic, chiropterophilous, columnar cactus of economic importance in Mexico. To investigate the effect of artificial selection on genetic diversity, inter-simple sequence repeat (ISSR) markers were used to estimate the genetic variation of wild, orchard, and backyard populations of S. queretaroensis from two regions in western Mexico. Six primers were used to generate 62 bands, of which 39 were polymorphic (62.9%). The total genetic diversity was similar in the wild (H_T = 0.296) and orchard (H_T = 0.291) groups, and slightly lower in the backyard group (H_T = 0.281). Wild populations (F_ST = 0.13) were less differentiated than backyard (F_ST = 0.17) and orchard (F_ST = 0.21) populations. The wild and backyard groups were genetically closer (0.015) than the wild and orchard (0.018) groups. A Mantel test revealed a positive correlation between genetic and geographic distances (r = 0.433, p = 0.002). In conclusion, gene flow and the prevailing management system have efficiently maintained genetic diversity and facilitated inter-population differentiation in S. queretaroensis.     

The genetic population structure and temporal genetic stability of gilthead sea bream Sparus aurata populations in the Aegean and Ionian Seas,using microsatellite DNA markers

We examined 662 gilthead sea bream Sparus aurata from wild samples of the species in the Aegean and Ionian Seas, using 20 EST-linked microsatellite markers, in three multiplex panels, as well as seven anonymous loci. Most of the markers were revealed to be highly polymorphic. We found low genetic differentiation between the sampling stations/areas with total F_ST 0.002 (P < 0.05). Based on comparison of five temporal samples, our results indicate genetic data consistency over time for all tested samples, pointing to stable populations, despite reported repeated escape events. Our results confirm the genetic population structure previously observed in these specific areas, using by far more markers than in previous studies in both coding and non-coding DNA loci. The limited genetic structure and the temporal genetic stability indicate neither major genetic differentiation of local populations by geographic isolation nor influence from anthropogenic factors. These results provide a baseline for future reference in any management programme of both wild and farmed population of S. aurata as well as of other aquaculture species with a potential introgression among farmed and wild populations.     

Morphological Plasticity and Phylogeny in a Monogenean Parasite Transferring between Wild and Reared Fish Populations

It is widely accepted that disease interactions between cultured and wild fish occur repeatedly, although reported cases have mainly relied just on the observation of similar symptoms in affected populations. Whether there is an explicit pathogen transfer between fish stocks, or each develops its own pathogen population, has been insufficiently studied and rarely supported by molecular tools. In this study, we used population dynamics and genetic structure of the monogenean Furnestinia echeneis in reared and neighbouring wild sea bream to indicate pathogen transfer, characterized by the phenotypic plasticity of the parasite attachment apparatus and the lack of phylogenetic differentiation. The observed pattern of genetic variation inferred by nuclear DNA Internal Transcribed Spacer 1 (ITS1) and mtDNA cytochrome C oxidase 1 (COI), between parasite populations is most likely caused by a recent shared demographic history like a reduced species area in the last glacial period. In spite of such recent expansion that populations underwent, F. echeneis shows differentiation in haptor morphometry as an adaptive trait in closely related populations at the aquaculture site. This suggests that differentiation in morphology may occur relatively rapidly in this species and that adaptive forces, not the speciation process, drives this monogenean parasitation. On the other hand, the observed phylogenetic inertia suggests a low to moderate gene flow (based on F_ST) between parasites in cultured and wild fish, evidencing for the first time the transfer of pathogens at the aquaculture site inferred by a molecular tool.     

QST–FST comparisons with unbalanced half‐sib designs

Q_ST, a measure of quantitative genetic differentiation among populations, is an index that can suggest local adaptation if Q_ST for a trait is sufficiently larger than the mean F_ST of neutral genetic markers. A previous method by Whitlock and Guillaume derived a simulation resampling approach to statistically test for a difference between Q_ST and F_ST, but that method is limited to balanced data sets with offspring related as half‐sibs through shared fathers. We extend this approach (i) to allow for a model more suitable for some plant populations or breeding designs in which offspring are related through mothers (assuming independent fathers for each offspring; half‐sibs by dam); and (ii) by explicitly allowing for unbalanced data sets. The resulting approach is made available through the R package QstFstComp.     

Inheritance of a darkened caudal peduncle and yellow body coloration in the gilthead sea bream (Sparus aurata L.)

This study presents data on inheritance of a darkened caudal peduncle (ebony) and yellow body coloration (yellow) in the gilthead sea bream (Sparus aurata). Fifteen progeny groups, obtained by crossing fish with three color phenotypes and of known origin, were analyzed. Analyzes of segregation in F₁ progeny involving groups from parental crosses of wild‐type colored × wild‐type colored; ebony × ebony; yellow × yellow, showed that the parents produced the offspring only with the same phenotypes (true breeding). Crosses involving F₁ wild‐type colored parents (that resulted from crosses of wild‐type parents with either ebony or yellow fish) showed in their F₂ progeny groups of which their phenotypic segregations did not differ significantly from a 3 : 1 Mendelian ratio. The progeny of back‐cross of ebony × (F₁ wild‐type colored × ebony) showed phenotypic segregations that did not differ significantly from the 1 : 1 Mendelian ratio. Overall, the results of the crossing experiments demonstrated that, similar to albinism described in a number of aquacultured species, ebony and yellow body coloration in S. aurata are both due to a single recessive allele. However, the yellow mutation of a gene controlling yellow pigment synthesis affects the yellow color of the whole fish body, whereas the ebony mutation causes production of melanin only in a specific area of the fish body, resulting in the development of a black coloration of the caudal peduncle. Experiments to assess culture performance showed that the color genes controlling ebony and yellow coloration had significant detrimental pleiotropic effects on growth, survival and body shape. Color mutations in the gilthead sea bream may be used as models for the study of: (i) genetic and physiological mechanisms of sterility, (ii) stress and disease resistance, (iii) effects of heterosis, (iv) genetic polymorphism in populations, and (v) methods of genetic protection in selected sea bream strains as well as in experiments on chromosome set manipulation.     

Identification of the causative gene of a transparent phenotype of juvenile red sea bream Pagrus major

Deformities in cultured fish species may be genetic, and identifying causative genes is essential to expand production and maintain farmed animal welfare. We previously reported a genetic deformity in juvenile red sea bream, designated a transparent phenotype. To identify its causative gene, we conducted genome-wide linkage analysis and identified two single nucleotide polymorphisms (SNP) located on LG23 directly linked to the transparent phenotype. The scaffold on which the two SNPs were located contained two candidate genes, duox and duoxa, which are related to thyroid hormone synthesis. Four missense mutations were found in duox and one in duoxa, with that in duoxa showing perfect association with the transparent phenotype. The mutation of duoxa was suggested to affect the transmembrane structure and thyroid-related traits, including an enlarged thyroid gland and immature erythrocytes, and lower thyroxine (T₄) concentrations were observed in the transparent phenotype. The transparent phenotype was rescued by T₄ immersion. Loss-of-function of duoxa by CRISPR–Cas9 induced the transparent phenotype in zebrafish. Evidence suggests that the transparent phenotype of juvenile red sea bream is caused by the missense mutation of duoxa and that this mutation disrupts thyroid hormone synthesis. The newly identified missense mutation will contribute to effective selective breeding of red sea bream to purge the causative gene of the undesirable phenotype and improve seed production of red sea bream as well as provide basic information of the mechanisms of thyroid hormones and its related diseases in fish and humans.Subject terms: Agricultural genetics, Animal breeding     

Population genetic structure of Penaeus monodon, in relation to monsoon current patterns in Southwest, East and Andaman coastal waters of India

The black tiger shrimp (Penaeus monodon), a commercially important penaeid species, is widely distributed across the Indo-Pacific region. Genetic diversity in P. monodon collected from eight geographical regions in Southwest, East and Andaman coastal waters of India (N = 418) was investigated using 10 polymorphic microsatellite loci. Average observed heterozygosity at sampled loci were high, ranging from 0.643 (Coromandel Coast) to 0.753 (South Andaman). Pairwise F_ST (ranged from 0.005 to 0.078) and R_ST (ranged from 0.005 to 0.171) estimates revealed surprisingly strong and statistically significant genetic structure among tiger shrimp populations. A synthetic map generated by multidimensional scaling shows an apparent cline in allele frequencies paralleling the roughly circular flow of surface currents in the Bay of Bengal. Significant heterozygote deficiencies were noted in most population samples at most loci. Andaman Island sites showed the highest diversity. Recognition of high genetic diversity and distinct population structuring of P. monodon in Indian seas has important implications for future domestication of this species in India, for two reasons: identification of the best wild founding stocks for aquaculture and, subsequently, the potential impacts of release of domesticates to the wild, either accidentally or deliberately (i.e. for stock enhancement).     

Evidence for the Formation of the Male Gynogenetic Fish

The females and unexpected males of gynogenetic red crucian carps (GRCC) with the 1:1 sex ratio were found in the progeny of the distant crossing of red crucian carp (RCC; ♀, 2n?=?100) × blunt snout bream (BSB; ♂, 2n?=?48). The females and males of GRCC were fertile, and they mated each other to generate the red crucian carps (GRCC₁) and another variational gray crucian carps (GGCC). The GRCC and their offspring were proved to be diploids (2n?=?100) with one to three microchromosomes by examining the chromosomal metaphases. The evidences for the male’s genetic effect in GRCC were provided by means of fluorescence in situ hybridization, Sox-HMG DNA markers, and microsatellite DNA markers. The genotypic variances of GRCC resulted in their phenotypic variances which were quite different from their maternal parent. It was concluded that the formation of the male gynogenetic fish in GRCC resulted from the genetic leakage of the paternal fish in the form of the microchromosomes including the paternal male-determining gene. After being activated by the sperm of BSB, which was inactivated and finally degraded but left the microchromosomes, the egg of RCC, in which the 50 chromosomes were spontaneously doubled to 100 chromosomes, developed into the diploid male gynogenetic fish. The formation of the bisexual GRCC and their progeny indicated that the distant hybridization could generate the bisexual diploid gynogenetic fish with genetic variation derived from the paternal fish, which is of great significance in both fish genetic breeding and evolutionary biology.     

High genetic diversity in an endangered medicinal plant, <Emphasis Type="Italic">Saussurea involucrata</Emphasis> (<Emphasis Type="Italic">Saussurea</Emphasis>, Asteraceae), in western Tianshan Mountains,China

Saussurea involucrata (Asteraceae) is a medicinal and second-degree national priority endangered plant that is mainly distributed in the high latitude region of the western Tianshan Mountains. The population is fragmented and isolated, and extensive human impact merits a suitable and specific conservation strategy, which can be compiled based on the genetic diversity, population structure, and demographic history. Phylogeographic studies were conducted on a total of five natural populations and 150 individuals were sampled. Data from three cpDNA intergenic spacer regions (trnL-F, matK, and ndhF-rpl32) and nrDNA ITS sequences showed that twelve haplotypes in cpDNA and five haplotypes in nrDNA indicated high genetic diversity among populations sampled (H _T?=?0.820 and 0.756) and within populations sampled (H _S?=?0.792 and 0.721). Additionally, the high genetic diversity did not mirror genetic structure in either cpDNA (F _ST?=?0.03153, G _ST?>?N _ST, p?<?0.05) or nrDNA (F _ST?=?0.03666, meaningless G _ST and N _ST). Two groups (north and south) were determined for a SAMOVA analysis. Based on this analysis, the demographic history was conducted with a Bayesian Skyline Plot and Isolation with Migration analysis, which showed sustainable and stable extension without a marked bottleneck. Divergence time was indicated at c. 6.25 Mya (90%HPD: 15.30–0.22 Mya) in the Miocene, which is consistent with the formation of the Kelasu section of Tianshan. The southern populations in the Bayanbulak and Gonglu regions require additional attention and transplanting would be an effective way to restore rare cpDNA haplotypes, increase effective population size, and migration rate. Our results suggested that in situ conservation of S. involucrata in western Tianshan should be the main strategy for protection and recovery of the species.