Kelussia odoratissima is well known for its medicinal importance. It has been announced as an endangered species. Thus, examining the genetic variation and conservation of this plant is necessary. In the present study, inter simple sequence repeat (ISSR) and sequence-related amplified polymorphism (SRAP) molecular markers were employed for the first time to access the genetic diversity and relationships of 77 wild individual plants of K. odoratissima collected from seven populations in Central Zagros region of Iran. A total of 146 bands were amplified by 12 ISSR primers, of which 129 (87.80 %) were polymorphic, while 69 polymorphic bands (83.30 %) were observed among 86 bands amplified by 11 SRAP primers. Polymorphic information content (PIC = 0.32), resolving power (Rp = 7.80), and marker informativeness (MI = 3.48) generated by ISSR primers were higher than that of SRAP analysis (PIC = 0.30, Rp = 5.61, and MI = 1.88). The study indicated that ISSR were more effective than SRAP markers for assessing the degree of genetic variation of K. odoratissima. In both UPGMA dendrograms of ISSR and SRAP, in most cases, individuals from each population were clustered in various groups without clear separation, which demonstrates the high variability of this germplasm in Iran. UPGMA cluster analysis revealed inconsistencies in the clustering patterns, as the Mantel’s test between the dendrograms for ISSR and SRAP data indicated a poor fit for the ISSR and SRAP data types (r = 0.10). Besides, principal coordinate analysis results showed that the first three principal coordinates account for 65.57 % of the total variation and studied seven populations were separated from each other and placed into five groups. These results have an important implication for K. odoratissima germplasm characterization, improvement, and conservation.