| Home | E-Submission | Sitemap | Editorial Office |  
Korean J. Pl. Taxon > Volume 53(1); 2023 > Article
/home/virtual/kjpt/journal//../xmls/kjpt-53-1-14.xml NYAMGEREL, BAASANMUNKH, OYUNTSETSEG, OYUNTSETSEG, NORRIS, CHOI, and BAYARMAA: Genetic diversity of the threatened Saussurea dorogostaiskii (Asteraceae) in the Khuvsgul region of Mongolia


Saussurea dorogostaiskii Palib. (Asteraceae) is a critically endangered medicinal plant in Mongolia and Russia. We studied the genetic variation of S. dorogostaiskii from three mountains of northern Mongolia. The genetic profile was assessed in 70 individuals from eight populations using five inter-simple sequence repeat markers, producing 53 loci with 96.4% polymorphism across all bands. Shannon’s index (I) and Nei’s gene diversity (H) value at the species level of S. dorogostaiskii are 0.25 and 0.17, respectively. An AMOVA showed high genetic variation among the populations (22% of populations and 32% of mountains), consistent with the high genetic differentiation (GST = 0.49) and low gene flow (Nm = 0.51) in S. dorogostaiskii populations. Eight populations were clustered into two groups, corresponding to their geographic locations. The low within-population genetic diversity and high genetic differentiation among S. dorogostaiskii populations factor into their endangered designation. This genetic analysis reveals that all populations are equally threatened, and community-based conservation is appropriate for these species.


Saussurea dorogostaiskii Palib. (Asteraceae) is one of the rare and threatened plants among the genus Saussurea DC. This species grows on the alpine talus of the southern Siberian to northern Mongolian mountains at altitudes from 1,500 to 2,750m a.s.l. (Von Raab-Straube, 2017; Norris, 2020; Baasanmunkh et al., 2022). In general, S. dorogostaiskii resembles S. orgaadayi Khanm. & Krasnob. and S. involucrata (Kar. & Kir.) Sch.Bip, but can be distinguished by its elongated racemiform synflorescence, oblong-elliptic to obovate leaves, and purple corolla (Von Raab-Straube, 2017). All species are referred to as ‘snow lotus’ in traditional medicine. Saussurea dorogostaiskii has long been used to treat lung diseases, cough, and cold symptoms (Krasnoborov, 1984; Norris, 2020). Recently, this species has been classified as critically endangered globally, with the area of occupancy is estimated at 68 km2 (Nyamgerel et al., unpubl. data) but the species population size over time is undocumented. This situation requires immediate action to sustainably manage its natural populations by developing and adapting appropriate conservation strategies.
Biodiversity conservation approaches focus on the geographic areas, ecosystems, ecological communities, and species, with less attention on genetic diversity and the evolutionary change among populations. Indeed, most conservation studies used species as a metric for biodiversity, focusing on intraspecific units within intensively managed threatened species (Coates et al., 2018). Intraspecific variation can determine the adaptive potential of species and influence its persistence over time (Barrett and Schluter, 2008). Most recent studies of plant species variation have identified high levels of diversity using DNA fingerprinting techniques (Godhe and Rynearson, 2017).
Inter simple sequence repeats (ISSR) markers are widely used DNA-based marker for population genetic studies because it generates more reliable and reproductive bands, as a compared than random amplification of polymorphic DNA (RAPD) (Nagaoka and Ogihara, 1997; Zhang and Dai, 2010). ISSR are technically simpler compared to simple sequence repeats and amplified fragment length polymorphism markers since they do not require prior sequence information to generate DNA application products (Joshi et al., 2000; Bornet and Branchard, 2001; Wang et al., 2008). Recently, ISSR markers are used to determine genetic variability and population relationship of plant species such as Prangos fedtschenkoi (Regel & Schmalh.) Korovin (Mustafina et al., 2017), Etlingera elatior (Jack) R.M.Sm. (Ismail et al., 2019), Mallotus oblongifolius (Miq.) Mull.Arg. (Yan et al., 2019) and Eleusine coracana (L.) Gaertn. (Venkatesan et al., 2021).
The genetic variability in Saussurea genus has been successfully investigated using genetic markers (Yuan et al., 2009; Gailīte et al., 2011; Jeong et al., 2012; Kreuzer et al., 2014; Dashzeveg et al., 2017; Sharma et al., 2019; Semwal et al., 2020; Yun and Kim, 2021). However, studies on genetic variation of S. dorogostaiskii are not studied to date. Therefore, the present study aimed to assess the genetic diversity and differentiation of S. dorogostaiskii in northern Mongolia using the ISSR markers. The knowledge of genetic variation will be helpful in the development of S. dorogostaiskii conservation strategy.



A total of 70 samples of S. dorogostaiskii from eight populations were collected from three isolated mountains in Khuvsgul province, namely, the Khoridol Saridag, Ulaan Taiga and Tengis-Shishget (Fig. 1A). Samples from protected areas were collected with permissions and supports from the relevant authorities. All individuals sampled within each population were separated by at least 10 m to avoid collecting individuals connected by a rhizome. Seven to ten individuals were randomly collected from each population due to its limited population size. The voucher specimens were deposited at the Herbarium of the National University of Mongolia (UBU). Detailed information on each population is given in Table 1.

DNA extraction and PCR amplification

Genomic DNA from the silica-gel-dried leaves was isolated by the CTAB methods (Doyle and Doyle, 1987). The quality and quantity of the genomic DNA were measured with spectrophotometry using a BioPhotometer (Eppendorf AG, Hamburg, Germany), and the integrity of the DNA was checked visually by ethidium bromide staining of a 0.8% agarose gel. DNA samples were diluted to a concentration of 50 ng/μL using a TE buffer. A total of 100 ISSR primers were initially screened with one individual from each population for reproducible amplification. Genetic profile study of S. dorogostaiskii was conducted by five ISSR primers, the remaining primers showed no amplification. Selected primers produced distinct and well-resolved fragments, which enable effective scoring for genetic diversity study within and among the populations.
The DNA amplifications were carried out in 50 μL volumes consisting of 2.0 μL DNA; 25 μL AmpliTaq Gold 360 Master Mix (Thermo Fisher Scientific Ltd., Waltham, MA, USA); 4 μL primer (800 nM); and the final volume was completed with ultrapure water. The PCR amplification was performed in a thermal cycler under the following conditions: initial denaturation at 94oC for 5 min; 35 cycles of 15 s at 94ºC for denaturation, 30 s at 47ºC for primer annealing and 2 min at 72ºC for extension, with a final extension for 5 min at 72ºC. The PCR products were separated by horizontal electrophoresis in agarose gel (2.0% w/v) immersed in a 1× TAE (Tris-Acetate-EDTA) buffer at constant voltage (150 V) for approximately 1 h and then stained with 0.5 μg/mL ethidium bromide. The DNA bands were then visualized with UV light and documented manually.

Data analysis

Distinct, reproducible fragments were scored as present (1) or absent (0) for each ISSR-PCR reaction and presented as part of a binary matrix. The number of total bands (TB), the number of polymorphic bands (PB), and the percentage of polymorphic bands were obtained. The following parameters were used to measure the efficiency of ISSR markers for genetic profile assessment. The polymorphic information content (PIC) is the potential of a primer to detect the polymorphism among individuals in the population, and it is calculated for each primer as described by Roldán-Ruiz et al. (2000): PICi = 2fi (1-fi), where PICi is the PIC of the locus i, fi is the frequency of the amplified allele, and 1-fi is the frequency of the non-amplified allele. The frequency was calculated as the proportion between the number of amplified alleles at each locus and the total number of accessions. The PIC of each primer was calculated using the average PIC value from all loci of each primer. The general usefulness of the markers used was determined by the marker index (MI) for each ISSR primer. The MI for each primer was calculated as the product of the PIC and effective multiplex ratio (EMR) according to Varshney et al. (2007); MI = PIC × EMR. The EMR is influenced by the overall utility of each marker, it was estimated as EMR = n × β, where n is the average number of the alleles amplified by accession to a specific system marker (multiplex ratio) and β was estimated after considering the number of polymorphic loci (PB) and the non-polymorphic loci (MB); β = PB/(PB + MB). The resolving power (RP) is the ability of the primer to differentiate between genotypes, which was calculated as RP = ΣIb, where Ib represents the informative bands (Prevost and Wilkinson, 1999). The Ib can be represented into a 0/1 scale by using the following formula: Ib = 1 – (2 × |0.5 – pi|), where pi is the ratio of accessions comprising the ith band.
Genetic parameters such as the observed number of alleles (Na), the effective number of alleles (Ne), the Shannon’s information index (I) and the genetic diversity or expected heterozygosity (H), were determined with GenALEx version 6.5 (Peakall and Smouse, 2012). Genetic divergence between populations was investigated by total genetic diversity (HT), genetic diversity within populations (HS), Nei’s (1978) genetic differentiation index among populations (GST), and gene flow estimates between populations (Nm) with POPGENE version 1.31 (Yeh et al., 1991).
An analysis of molecular variance (AMOVA) was performed to estimate the partitioning of genetic variation within and among populations using GenALEx 6.5 (Peakall and Smouse, 2012). A study on the correlation (Mantel test) between the genetic and geographic distance among the populations using XLSTAT software (Addinsoft, 2021). The unbiased genetic distance was utilized for the construction of a dendrogram using unweighted pair group method with arithmetic mean (UPGMA) in POPGENE program version 1.31. In addition, Nei’s genetic distance matrix was used to construct a Neighbor-Joining dendrogram, using XLSTAT, for the 70 individuals belonging to the eight populations of S. dorogostaiskii. Also, the principal coordinate analysis (PCoA) was also performed using XLSTAT to determine the relative genetic distance between individuals and to check the consistency of population genetic differentiation as defined by cluster analysis.


ISSR marker analysis

Five primers that displayed clear and reproductive bands were obtained through the screening of a total of 100 primers. A total of 53 amplified ISSR fragments (loci) were found from 70 S. dorogostaiskii samples with 96.4% polymorphism and an average of 10.6 loci by primer. The band’s range was from 150 bp to 1,000 bp. The information on the genetic profile of each sample obtained using the five ISSR primers was used to estimate marker performance using the following parameters (Table 2). A high PIC value was detected for primer UBC891 at 0.41 and a low PIC value of 0.17 for UBC810. The average PIC value per primer is 0.24. The maximum EMR (7.52) was detected with the primer UBC855 and the lowest was shown by the primer UBC840 (1.18), with a mean EMR of 5.14 per primer. The average MI was 3.04 per ISSR primers. The highest MI value was detected with the primer UBC810 (4.44) and the lowest with UBC840 (1.5). The average RP was 5.75 per ISSR primer. The highest RP value was detected with the primer UBC836 (9.28), and the lowest with the primer UBC855 (1.94).

Genetic diversity

Estimates of genetic diversity parameters varied among S. dorogostaiskii populations (Table 3). The percentage of polymorphic loci for each population ranged from 0.30% (HG; Mt. Khug, Ulaan Taiga, Ulaan-Uul) to 0.60% (UR; Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe) with an average of 0.46%. The UR population included the highest number of effective alleles per locus (1.40) and the lowest alleles were observed in the HG population (0.85). The UR population from the Khoridol Saridag Mountains proved to be the most genetically diverse (I = 0.33, H = 0.23), while the HG population from the Ulaan Taiga Mountains had the lowest (I = 0.16, H = 0.11) of all populations. The average genetic diversity of S. dorogostaiskii was low (I = 0.25, H = 0.17).

Genetic differentiation

AMOVA revealed a statistically significant difference (p < 0.01) in variance among eight populations of S. dorogostaiskii. Of the total genetic diversity, 22% occurred among regions, 32% among populations, and 46% within populations (Table 4). The additional genetic differentiation parameters support AMOVA result (Table 5). The genetic diversity within populations (HS = 0.16) was equal to half of the total genetic diversity (HT = 0.32). The observed genetic differentiation among populations (GST) was 0.49 demonstrating that 49% of genetic variation is among the populations. This result was further substantiated by the low level of gene flow among the populations, as showed by the low Nm value (0.51). The mantel test showed no significant correlation between population genetic differentiation and geographical distance (r = 0.098, p = 0.632) (Fig. 2). This indicates that genetic variation among S. dorogostaiskii is not limited by the remoteness of populations.

Relationship among populations

The pairwise genetic distance and genetic identity among eight populations of S. dorogostaiskii ranged from 0.008 to 0.40 and 0.67 to 0.92, respectively (Table 6). The maximum genetic distance was observed between BJ (Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh)–UR pair. The least genetic distance was found in TL (Mt. Tag, Khoridol Saridag, Ulaan-Uul sum)–HT (Mt. Khorin, Khoridol Saridag, Renchinlkhumbe) and MO (Mt. Mongorog, Ulaan Taiga, Ulaan-Uul)–HG pairs. The dendrogram of the eight populations of S. dorogostaiskii was based on the genetic distances of Nei’s using UPGMA cluster analysis (Fig. 3). UPGMA analysis separated out two main clusters. Cluster I included TL, TS (Mt. Toshin, Khoridol Saridag, Renchinlkhumbe), HT, and UR collected from the Khoridol Saridag Mountains, and SO (Mt. Soyo, Tenger Shishget, Tsaaannuur) from the Tengis-Shishget Mountains. Cluster II was from the HG, MO, and BJ from Ulaan Taiga Mountains.
PCoA was obtained for the individuals and the differentiation was divided into three main clusters (Fig. 4): Cluster 1 included populations TL, TS, HT, and SO, consistent with cluster I of the UPGMA analysis; Ulaan Taiga Mountains populations were divided into Cluster 2 (HG and MO) and Cluster 3 (BJ). The Neighbor-Joining tree for individuals based on Nei’s genetic distance showed that most individuals from the same population origin were formed together (Fig. 5). However, the individuals from the Khoridol Saridag Mountains were formed in partially mixed clusters. Additionally, the SO population of the Tengis-Shishget region is more similar to the populations of the Khoridol Saridag Mountains than to the populations of the Ulaan Taiga Mountains.


The primary goal in conserving endangered plant species is maintaining genetic diversity (Margules and Pressey, 2000). The S. dorogostaiskii, a globally threatened species, has low genetic diversity, which was similar to the previous findings in the Saussurea genus that also used the ISSR markers (Jeong et al., 2012; Semwal et al., 2020). In addition, the average genetic diversity of S. dorogostaiskii was lower (H = 0.17, I = 0.25) than the short-lived perennial plants (H = 0.32) and endemic species (H = 0.20) estimated using RAPD markers (Nybom, 2004). The proportion of genetic diversity within the populations was low in S. dorogostaiskii (46%) compared to most Saussurea species (59–92%) (Gailīte et al., 2011; Jeong et al., 2012; Kreuzer et al., 2014; Dashzeveg et al., 2017; Semwal et al., 2020; Yun and Kim, 2021).
Genetic diversity within the population is influenced by many factors such as population size, breeding system, natural selection, and life history traits (Li et al., 2018). Previous studies reported that S. dorogostaiskii is a monocarpic perennial that reproduces exclusively by seed; it has been found with more rocks and less vegetation cover in the middle elevation of mountains up to 2,735 m, with a relatively low number of mature individuals (Von Raab-Straube, 2017; Norris, 2020). In addition, Boikov (2009) discovered that S. dorogostaiskii, as a rare plant species adapted to ecological stress conditions in Siberia, is not very competitive, grows slowly, and prefers an open canopy area. For a species already low in number, these traits of reproducing slowly and in highstress niche environments, indicate that gene flow would be limited. In our findings, the gene flow (Nm) between S. dorogostaiskii populations was 0.51. Generally, Nm < 1 is insufficient to counteract the effects of genetic drift (Hamrick et al., 1995). Genetic drift affects the distribution of genetic variation in two ways: (1) the decrease of variation within populations (loss of heterozygosity and eventual fixation of alleles), and (2) the increase of differentiation among populations (Ellstrand and Elam, 1993). Additionally, we calculated the GST values to define genetic differences among populations. This value is often very low, but statistically significant, it is classified as low (< 0.05), medium (0.05–0.15), or high (> 0.15) (Li et al., 2018). S. dorogostaiskii populations have a high level (0.49) of differentiation, which occurs when gene flow between populations is restricted.
Looking at the clustering of population genetics can help reveal what important factors restrict gene flow. In this study, all Ulaan-Taiga populations formed one cluster, and other populations formed a second main cluster. The Tengis-Shishget Mountains is linked by connecting terrain and proximity to the northern edge of the Khoridol Saridag Mountains. This connection explains why the SO population is more genetically alike to the Khoridol Saridag populations than the Ulaan-Taiga populations. The Ulaan Taiga and Khoridol Saridag Mountains are separated by the Darkhad Depression (Fig. 1A) (Clark, 2014), is a significant area and entirely different ecosystem which could limit gene exchange between S. dorogostaiskii populations. This trend of genetic isolation among snow lotus species similarly occurred at high elevations in China (Chen and Wang, 2018). These species are distributed in the Tianshan Mountains and the Altai Mountains, which are separated by the Junggar Basin. Furthermore, the BJ population is genetically differentiated from the other two populations in the same area. This difference may be because the BJ population is located in very different ecological conditions, e.g., in a canyon on the top of a mountain range, and that this population had only three adults and less than 50 individuals. In small populations (e.g., < 100 individuals), allele frequencies may undergo large and unpredictable fluctuations due to random genetic drift (Barrett and Kohn, 1991).
In the natural populations of S. dorogostaiskii, significant loss of genetic diversity within the populations due to genetic drift caused by small population sizes, ecological stress conditions and geographical isolation in northern Mongolia. The population is the fundamental unit of species evolution, low genetic diversity within-population may reduce adaptability to changing environments (Li et al., 2018). Therefore, the protection of all existing populations is necessary, required to maintain genetic variability and protect it particularly in the face of alpine climate change. Furthermore, continuing studies will require assessing the genetic diversity in southern Siberian S. dorogostaiskii populations and comparing them to northern Mongolian populations to see if these patterns of genetic diversity and isolation are similar.
Most populations of S. dorogostaiskii are found in federally protected areas and harvesting is prohibited by law in Mongolia. However, because it grows at high altitudes in vast areas (approximately 1,528,000 hectares), rangers have rarely been able to enforce this law and catch the harvesters (Norris et al., 2020). Harvesting during the flowering season is reported to be a significant cause of S. dorogostaiskii endangerment. According to this previously mentioned study, the most effective ways to protect S. dorogostaiskii are to educate locals about its rarity and threats, increase safeguarding time during the flowering season, and consider community cultivation programs for the species harvest. Furthermore, harvesting of the Himalayan snow lotus correlated with decreased plant height and status over a 100 years period. This situation suggests that humans can unconsciously drive evolution by selecting the largest individuals for harvest; this behavioral pattern, should be considered when managing threatened species (Law and Salick, 2005).


This is the first study to consider the genetic profile of S. dorogostaiskii in northern Mongolia. The ISSR markers revealed relatively lower genetic diversity for S. dorogostaiskii than for other Saussurea species and rare medicinal plants. This species has quite low diversity within the population and high genetic differentiation among the populations. This is likely because, the populations of S. dorogostaiskii are affected by genetic drift due to high elevation isolation and the expanse of the Darkhad Depression. Therefore, all populations of S. dorogostaiskii equally warrant an active conservation strategy to maintain their natural diversity. However, forthcoming studies with genome wide genotyping-by- sequencing (Elshire et al., 2011) or MIG-Seq (multiplexed ISSR genotyping by sequencing) (Suyama and Matsuki, 2015) will be valuable in better understanding of the genetic diversity and distribution patterns of S. dorogostaiskii with related Saussurea species known as Snow Lotus group.


Our study was supported by a research project (Population genetic study of Saussurea dorogostaiskii in Khuvsgul region; Grant P2019-3758) of National University of Mongolia, Mongolia, Korea National Arboretum, Republic of Korea (A study on the distribution of vascular plants in Mongolia; Grant KNA1-2-38, 20-5), and the University of Vermont, USA (Ian Worley Award and the Simon Family Public Research Fellowship). We are grateful to Tsegmed Z. for his assistance with the geographical map. We also thank Tumursukh J., Davaasuren B., Byambaa, Batzaya, Nyam-Ochir D., Tsogoo, Darmaabazar U. and UTPAA employees and rangers for their helpful communication assistance during our field study.


The authors declare that there are no conflicts of interest.

Fig. 1.
Saussurea dorogostaiskii in Khuvsgul region, Mongolia. A. Distribution of eight populations. B. Mature plant in Tag Mountain, Khoridol Saridag. TL, Mt. Tag, Khoridol Saridag, Ulaan-Uul sum; TS, Mt. Toshin, Khoridol Saridag, Renchinlkhumbe; HT, Mt. Khorin, Khoridol Saridag, Renchinlkhumbe; SO, Mt. Soyo, Tenger Shishget, Tsaaannuur; UR, Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe; BJ, Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh; HG, Mt. Khug, Ulaan Taiga, Ulaan-Uul; MO, Mt. Mongorog, Ulaan Taiga, Ulaan-Uul.
Fig. 2.
A graph comparing geographic distance and genetic differentiation among eight Saussurea dorogostaiskii populations.
Fig. 3.
UPGMA dendrogram (with length of branch) based on Nei’s unbiased measurements of genetic distance between eight populations of Saussurea dorogostaiskii. TL, Mt. Tag, Khoridol Saridag, Ulaan-Uul sum; HT, Mt. Khorin, Khoridol Saridag, Renchinlkhumbe; TS, Mt. Toshin, Khoridol Saridag, Renchinlkhumbe; UR, Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe; SO, Mt. Soyo, Tenger Shishget, Tsaaannuur; BJ, Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh; HG, Mt. Khug, Ulaan Taiga, Ulaan-Uul; MO, Mt. Mongorog, Ulaan Taiga, Ulaan-Uul.
Fig. 4.
Two-dimensional plot of principal coordinate analysis showing clustering of individuals belonging to eight populations of Saussurea dorogostaiskii. TL, Mt. Tag, Khoridol Saridag, Ulaan-Uul sum; TS, Mt. Toshin, Khoridol Saridag, Renchinlkhumbe; HT, Mt. Khorin, Khoridol Saridag, Renchinlkhumbe; UR, Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe; SO, Mt. Soyo, Tenger Shishget, Tsaaannuur; BJ, Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh; HG, Mt. Khug, Ulaan Taiga, Ulaan-Uul; MO, Mt. Mongorog, Ulaan Taiga, Ulaan-Uul.
Fig. 5.
Neighbor-Joining clustering pattern for 70 individuals belonging to eight populations of Saussurea dorogostaiskii.
Table 1.
The detailed sample information of eight populations of S. dorogostaiskii in Khuvsgul province, Mongolia, with abbreviation of location, collection site, and sample size (Ns).
No. Abbreviation Collection site Altitude (m) Ns
1 TL Mt. Tag, Khoridol Saridag, Ulaan-Uul sum 2,538 9
2 TS Mt. Toshin, Khoridol Saridag, Renchinlkhumbe 2,396 7
3 HT Mt. Khorin, Khoridol Saridag, Renchinlkhumbe 2,083 9
4 SO Mt. Soyo, Tenger Shishget, Tsaaannuur 1,708 7
5 UR Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe 2,134 9
6 BJ Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh 2,322 9
7 HG Mt. Khug, Ulaan Taiga, Ulaan-Uul 1,609 10
8 MO Mt. Mongorog, Ulaan Taiga, Ulaan-Uul 2,731 10
Table 2.
Sequences of primers, number of total bands (TB), polymorphic bands (PB), approximate size range (in base pairs), and genetic characterization of 5 ISSR loci for 70 accessions of Saussurea dorogostaiskii.
Primer Sequence 5′ → 3′ Amplified bands PIC EMR MI RP

TB PB (%) Size range (bp)
1 UBC 855 (AC)8YT 11 9 (81.8) 150–450 0.23 7.52 1.71 1.94
2 UBC 891 HVH(TG)8 9 9 (100) 200–1,000 0.41 5.64 3.32 6.74
3 UBC 836 (AG)8YA 12 12 (100) 150–1,000 0.21 5.28 4.20 9.28
4 UBC 840 (GA)8YT 6 6 (100) 150–500 0.17 1.18 1.50 2.12
5 UBC 810 (GA)8T 15 15 (100) 300–1,000 0.16 5.44 4.44 8.66

Total 53
Mean 10.6 10.2 (96.4) 0.24 5.14 3.04 5.75

Y, single letter abbreviation for mixed base positions for C and T; H, single letter abbreviation for mixed base positions for A, C and T; V, single letter abbreviation for mixed base positions for G, C and A; PIC, polymorphic information content; EMR, effective multiplex ratio; MI, marker index; RP, resolving power.

Table 3.
Genetic diversity within populations parameters of eight populations of Saussurea dorogostaiskii.
Populations % P I H Ne Na
TL 0.55 0.31 0.21 1.34 1.38
TS 0.51 0.30 0.21 1.30 1.36
HT 0.51 0.28 0.19 1.30 1.34
UR 0.60 0.33 0.23 1.40 1.41
SO 0.50 0.27 0.19 1.30 1.34
BJ 0.36 0.19 0.13 1.09 1.23
HG 0.30 0.16 0.11 0.85 1.19
MO 0.36 0.17 0.11 0.98 1.17
Mean 0.46 0.25 0.17 1.20 1.30

% P, percentage of polymorphic loci; I, Shannon’s Information Index; H, expected heterozygosity; Ne, effective number of alleles; Na, observed number of alleles; TL, Mt. Tag, Khoridol Saridag, Ulaan-Uul sum; TS, Mt. Toshin, Khoridol Saridag, Renchinlkhumbe; HT, Mt. Khorin, Khoridol Saridag, Renchinlkhumbe; UR, Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe; SO, Mt. Soyo, Tenger Shishget, Tsaaannuur; BJ, Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh; HG, Mt. Khug, Ulaan Taiga, Ulaan-Uul; MO, Mt. Mongorog, Ulaan Taiga, Ulaan-Uul.

Table 4.
Analysis of molecular variance (AMOVA) for five ISSR markers among Saussurea dorogostaiskii populations.
Source of variance df VC Variation (%)
Among mountains 2 2.230 22**
Among populations 5 3.345 32**
Within populations 62 4.724 46**

Total 69 10.300 100

ISSR, inter simple sequence repeats; df, degree of freedom; VC, variance components.

** p < 0.01.

Table. 5
Genetic differentiation parameters among the eight Saussurea dorogostaiskii populations.
S. dorogostaiskii HT HS GST Nm*
Total 0.32 0.16 0.49 0.51

HT, total genetic diversity; HS, genetic diversity within populations; GST, relative magnitude of genetic differentiation among populations; Nm, estimate of gene flow from GST.

Table 6.
Nei’s unbiased measures of genetic identity (above right) and genetic distance (below left) of eight Saussurea dorogostaiskii populations.
Populations TL TS HT UR SO BJ HG MO
TL 0.89 0.92 0.91 0.82 0.73 0.75 0.76
TS 0.11 0.84 0.87 0.84 0.75 0.79 0.76
HT 0.08 0.17 0.88 0.79 0.68 0.74 0.73
UR 0.09 0.14 0.12 0.78 0.67 0.74 0.71
SO 0.20 0.18 0.24 0.25 0.74 0.76 0.72
BJ 0.32 0.29 0.39 0.40 0.30 0.80 0.80
HG 0.29 0.23 0.30 0.30 0.27 0.22 0.92
MO 0.27 0.27 0.32 0.34 0.33 0.22 0.08

TL, Mt. Tag, Khoridol Saridag, Ulaan-Uul sum; TS, Mt. Toshin, Khoridol Saridag, Renchinlkhumbe; HT, Mt. Khorin, Khoridol Saridag, Renchinlkhumbe; UR, Mt. Urtrag, Khoridol Saridag, Renchinlkhumbe; SO, Mt. Soyo, Tenger Shishget, Tsaaannuur; BJ, Mt. Baruun Jarga, Ulaan Taiga, Bayanzurkh; HG, Mt. Khug, Ulaan Taiga, Ulaan-Uul; MO, Mt. Mongorog, Ulaan Taiga, Ulaan-Uul.


Addinsoft. 2021. XLSTAT statistical and data analysis solution, version 7.0. Long Island, New York, USA.

Akshay, K. R. Sudharani, N. Anjali, K. B and Deepak, T. M. 2014. Biodiversity and strategies for conservation of rare, endangered and threatened medicinal plants. Research & Reviews: Journal of Pharmacognosy and Phytochemistry 2: 12-20.

Baasanmunkh, S. Urgamal, M. Oyuntsetseg, B. Sukhorukov, A. P. Tsegmed, Z. Son, D. C. Erst, A. Oyundelger, K. Kechaykin, A. A. Norris, J. Kosachev, P. Ma, J.-S. Chang, K. S and Choi, H. J. 2022. Flora of Mongolia: annotated checklist of native vascular plants. PhytoKeys 192: 63-169.
pmid pmc
Barrett, S. C. H and Kohn, J. R. 1991. Genetic and evolutionary consequences of small population size in plants: implications for conservation. Genetics and Conservation of Rare Plants. Falk, D. A. Holsinger, K. E (eds.), Oxford University Press, New York. 3-30.

Barrett, R. D. H and Schluter, D. 2008. Adaptation from standing genetic variation. Trends in Ecology and Evolution 23: 38-44.
crossref pmid
Boikov, T. G. 2009. Functioning of rare plant populations of Central Siberia. Contemporary Problems of Ecology 2: 444-450.
crossref pdf
Bornet, B and Branchard, M. 2001. Nonanchored inter simple sequence repeat (ISSR) markers: reproducible and specific tools for genome fingerprinting. Plant Molecular Biology Reporter 19: 209-215.
crossref pdf
Chen, J and Wang, Y.-J. 2018. New Saussurea (Asteraceae) species from Bogeda Mountain, eastern Tianshan, China, and inference of its evolutionary history and medical usage. PLoS ONE 13: e0199416.
crossref pmid pmc
Clark, J. K. 2014. Modeling late prehistoric and early historic pastoral adaptations in northern Mongolia’s Darkhad Depression. PhD dissertation,. University of Pittsburgh, Pittsburgh, PA, USA. 275 pp.

Coates, D.J. Byrne, M and Moritz, C. 2018. Genetic diversity and conservation units: Dealing with the species-population continuum in the age of genomics. Frontiers in Ecology and Evolution 6: 165.
Dashzeveg, N. Buerkert, A and Wiehle, M. 2017. Morphological and genetic diversity and seed germination behavior of a snow lotus (Saussurea involucrata, Asteraceae) from the Mongolian Altay Mountains, Western Mongolia. Genetic Resources and Crop Evolution 64: 927-934.
crossref pdf
Doyle, J. J and Doyle, J. L. 1987. A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 19: 11-15.

Ellstrand, N. C and Elam, D. R. 1993. Population genetic consequences of small population size: Implications for plant conservation. Annual review of Ecology and Systematics 24: 217-242.
Elshire, R. J. Glaubitz, J. C. Sun, Q. Poland, J. A. Kawamoto, K. Buckler, E. S and Mitchell, S. E. 2011. A robust, simple geno-typing- by-sequencing (GBS) approach for high diversity species. PLoS ONE 6: e19379.
Gailīte, A. Ievinsh, G and Ruņģis, D. 2011. Genetic diversity analysis of Latvian and Estonian Saussurea esthonica populations. Environmental and Experimental Biology 9: 115-119.

Godhe, A and Rynearson, T. 2017. The role of intraspecific variation in the ecological and evolutionary success of diatoms in changing environments. Philosophical Transactions of the Royal Society B: Biological Sciences 372: 20160399.
crossref pmid pmc pdf
Hamrick, J. L. Godt, M. J. W and Sherman-Broyles, S. L. 1995. Gene flow among plant population: Evidence from genetic markers. Experimental and Molecular Approaches to Plant Biosystematics. Hoch, P. C. Stephenson, A. G (eds.), Missouri Botanical Garden Press, St. Louis, MO. 215-232.

Ismail, N. A. Rafii, M. Y. Mahmud, T. M. M. Hanafi, M. M and Miah, G. 2019. Genetic diversity of torch ginger (Etlingera elatior) germplasm revealed by ISSR and SSR markers. BioMed Research International 2019. 5904804:
crossref pmid pdf
Jeong, J. H. Lee, B. C. Yoo, K. O. Jang, S. K and Kim, Z. S. 2012. Influence of small-scale habitat patchiness on the genetic diversity of the Korean endemic species Saussurea chabyoungsanica (Asteraceae). Biochemical Systematics and Ecology 43: 14-24.
Joshi, S. P. Gupta, V. S. Aggarwal, R. K. Ranjekar, P. K and Brar, D. S. 2000. Genetic diversity and phylogenetic relationship as revealed by inter simple sequence repeat (ISSR) polymorphism in the genus Oryza . Theoretical and Applied Genetics 100: 1311-1320.
crossref pdf
Krasnoborov, I. M. 1984. Key to Plants of the Tuva ASSR. Nauka Publisher, Novosibirsk. 335 (in Russian).

Kreuzer, M. Tribsch, A and Nyffeler, R. 2014. Ecological and genetic differentiation of two subspecies of Saussurea alpina in the Western Alps. Alpine Botany 124: 49-58.
crossref pdf
Law, W and Salick, J. 2005. Human-induced dwarfing of Himalayan snow lotus, Saussurea laniceps (Asteraceae). Proceedings of the National Academy of Sciences of the United States of America 102: 10218-10220.
crossref pmid pmc
Li, S. Gan, X. Han, H. Zhang, X and Tian, Z. 2018. Low withinpopulation genetic diversity and high genetic differentiation among populations of the endangered plant Tetracentron sinense Oliver revealed by inter-simple sequence repeat analysis. Annals of Forest Science 75: 74.
crossref pdf
Margules, C. R and Pressey, R. L. 2000. Systematic conservation planning. Nature 405: 243-253.
crossref pmid pdf
Mustafina, F. U. Kim, E. H. Son, S.-W. Turginov, O. T. Chang, K. S and Choi, K. 2017. Assessment of genetic diversity of Prangos fedtschenkoi (Apiaceae) and its conservation status based on ISSR markers. Korean Journal of Plant Taxonomy 47: 11-22.

Nagaoka, T and Ogihara, Y. 1997. Applicability of inter-simple sequence repeat polymorphisms in wheat for use as DNA markers in comparison to RFLP and RAPD markers. Theoretical and Applied Genetics 94: 597-602.
crossref pdf
Nei, M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89: 583-590.
crossref pmid pmc pdf
Norris, J. Nudkhuu, N. Jal, T. Davaajii, L. Buyanjargal, D. Dvvjii, N.-O. Ornobish, B. Erdenebat, B and Ornobish, T. 2020. Synthesis of the local knowledge and ecology of vansemberuu (Saussurea dorogostaiskii Palib.) for conservation, Darkhad Valley, Khuvsgul, Mongolia. UVM Honors College Senior Theses. (358): UVM ScholarWorks; 75.

Nybom, H. 2004. Comparison of different nuclear DNA markers for estimating intraspecific genetic diversity in plants. Molecula Ecology 13: 1143-1155.
crossref pmid
Peakall, R and Smouse, P. E. 2012. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research: An update. Bioinformatics 28: 2537-2539.
crossref pmid pmc pdf
Prevost, A and Wilkinson, M. J. 1999. A new system of comparing PCR primers applied to ISSR fingerprinting of potato cultivars. Theoretical and Applied Genetics 98: 107-112.
crossref pdf
Roldán-Ruiz, I. Dendauw, J. Van Bockstaele, E. Depicker, A and De Loose, M. 2000. AFLP markers reveal high polymorphic rates in ryegrasses (Lolium spp.). Molecular Breeding 6: 125-134.

Semwal, P. Painuli, S. Jugran, A. Palni, L. M. S. Verma, S. Thapliyal, A and Tewari, D. 2020. Genetic diversity of scanty available Himalayan Saussurea obvallata (DC.) Edgew. Iranin Journal of Science and Technology, Transactions A: Science 44: 587-594.
crossref pdf
Sharma, S. Sharma, R. Sharma, P. Thakur, K and Dutt, B. 2019. Direct shoot organogenesis from seedling derived shoot tip explants of endangered medicinal plant Saussurea costus (Falc.) Lipsch. Proceedings of the National Academy of Sciences India Section B: Biological Sciences 89: 755-764.
crossref pdf
Suyama, Y and Matsuki, Y. 2015. MIG-seq: An effective PCRbased method for genome-wide single-nucleotide polymorphism genotyping using the next-generation sequencing platform. Scientific Reports 5: 16963.
crossref pdf
Varshney, R. K. Chabane, K. Hendre, P. S. Aggarwal, R. K and Graner, A. 2007. Comparative assessment of EST-SSR, EST-SNP and AFLP markers for evaluation of genetic diversity and conservation of genetic resources using wild, cultivated and elite barleys. Plant Science 173: 638-649.
Venkatesan, J. Ramu, V. Sethuraman, T. Sivagnanam, C and Doss, G. 2021. Molecular marker for characterization of traditional and hybrid derivatives of Eleusine coracana (L.) using ISSR marker. Journal of Genetic Engineering and Biotechnology 19: 178.
crossref pmid pmc pdf
Von Raab-Straube, E. 2017. Taxonomic revision of Saussurea subgenus Amphilaena (Compositae, Cardueae). Englera 34: 1-274.

Wang, M. X. Zhang, H. L. Zhang, D. L. Qi, Y. W. Fan, Z. L. Li, D. Y. Pan, D. J. Cao, Y. S. Qiu, Z. E. Yu, P. Yang, Q. W. Wang, X. K and Li, Z. C. 2008. Genetic structure of Oryza rufipogon Griff in China. Heredity 101: 527-535.
crossref pmid pdf
Yan, W. Li, J. Zheng, D. Friedman, C and Wang, H. 2019. Analysis of genetic population structure and diversity in Mallotus oblongifolius using ISSR and SRAP markers. PeerJ 7: e7173.
crossref pmid pmc pdf
Yeh, C. F. Yang, R. C and Boyle, T. 1991. POPGENE: Microsoft Windows-Based Freeware for Population Genetics Analysis, Version 1.31. University of Alberta and Centre for International Forestry Research, Alberta, Canada.

Yun, S.A and Kim, S.-C. 2021. Genetic diversity and structure of Saussurea polylepis (Asteraceae) on continental islands of Korea: Implications for conservation strategies and management. PLoS ONE 16: e0249752.
crossref pmid pmc
Yuan, X. F. Dai, Z. H. Wang, X. D and Zhao, B. 2009. Assessment of genetic stability in tissue-cultured products and seedlings of Saussurea involucrata by RAPD and ISSR markers. Biotechnology Letters 31: 1279-1287.
crossref pmid pdf
Zhang, L.-J and Dai, S.-L. 2010. Genetic variation within and among populations of Orychophragmus violaceus (Cruciferae) in China as detected by ISSR analysis. Genetic Resources and Crop Evolution 57: 55-64.
crossref pdf
PDF Links  PDF Links
PubReader  PubReader
ePub Link  ePub Link
Full text via DOI  Full text via DOI
Download Citation  Download Citation
Web of Science
Related article
Editorial Office
Korean Journal of Plant Taxonomy
Department of Biology, Daejeon University, Daejeon 34520, Korea
TEL: +82-42-280-2434   E-mail: kjpt1968@gmail.com
About |  Browse Articles |  Current Issue |  For Authors and Reviewers
Copyright © Korean Society of Plant Taxonomists.                 Developed in M2PI