The conservation biology field underscores the importance of understanding genetic diversity and gene flow within plant populations and the factors that influence them. This study employs Simple Sequence Repeat (SSR) molecular markers to investigate the genetic diversity of the endangered plant species Saussurea involucrata, offering a theoretical foundation for its conservation efforts. Utilizing sequencing results to screen SSR loci, we designed and scrutinized 18 polymorphic microsatellite primers across 112 samples from 11 populations in the Bayinbuluke region. Our findings reveal high genetic diversity (I = 0.837, He = 0.470) and substantial gene flow (Nm = 1.390) among S. involucrata populations (China, Xinjiang), potentially attributed to efficient pollen and seed dispersal mechanisms. Principal Coordinate Analysis (PCoA) indicates a lack of distinct genetic structuring within the Bayinbuluke populations. The cluster analysis using STRUCTURE reflected the genetic structure of S. involucrata to a certain extent compared with PCoA. The results showed that all samples were divided into four groups. To safeguard this species, we advocate for the in situ conservation of all S. involucrata populations in the area. The SSR markers developed in this study provide a valuable resource for future genetic research on S. involucrata.

Rare and vulnerable endemic plants represent different evolutionary units that occur at different times, and protecting these species is a key issue in biological protection. Understanding the impact of the history of endangered plant populations on their genetic diversity helps to reveal evolutionary history and is crucial for guiding conservation efforts. Saussurea involucrata, a perennial alpine species mainly distributed in the Tianshan Mountains, is famous for its medicinal value but has become endangered due to over-exploitation. In the present study, we employed both nuclear and chloroplast DNA sequences to investigate the genetic distribution pattern and evolutionary history of S. involucrata. A total of 270 individuals covering nine S. involucrata populations were sampled for the amplification and sequencing of nrDNA Internal Transcribed Spacer (ITS) and chloroplast trnL-trnF, matK and ndhF-rpl32 sequences. Via calculation, we identified 7 nuclear and 12 plastid haplotypes. Among the nine populations, GL and BA were characterized by high haplotype diversity, whereas BG revealed the lowest haplotype diversity. Molecular dating estimations suggest that divergence among S. involucrata populations occurred around 0.75 Ma, coinciding with the uplift of Tianshan Mountains. Our results reveal that both isolation-by-distance (IBD) and isolation-by-resistance (IBR) have promoted genetic differentiation among populations of S. involucrata. The results from the ecological niche modeling analyses show a more suitable habitat for S. involucrata in the past than at present, indicating a historical distribution contraction of the species. This study provides new insight into understanding the genetic differentiation of S. involucrata, as well as the theoretical basis for conserving this species.