Peganum harmala L. is a perennial herbaceous plant and can be a future drug due to its wide medicinal purposes. Despite its economic importance, the molecular genetics of P. harmal have not yet been studied in detail. Genetic diversity of 12 P. harmala genotypes were investigated by using Inter-Simple Sequence Repeats (ISSR), PCR-RFLP of rDNA-ITS, PCR-SSCP of rDNA-ITS and Simple Sequence Repeat (SSR) markers. The level of polymorphism revealed by ITS-SSCP is the lowest, followed by ITS-RFLP then ISSR and the highest polymorphism level was reported for SSR marker. The AMOVA analysis implied that most of the variation occurred within the Populations. A value of inbreeding coefficient Fis estimated by the three co-dominant markers was nearly equal and offer an indication of the partial out-crossing reproductive system of P. harmala. Principal Coordinate Analysis (PCOA) plot revealed a clear pattern of clustering based on the locations of collected plants which coincide with the isolation by distance. The study revealed that ITS-SSCP and ISSR markers respectively were more informative than the other used markers in the assessment of genetic diversity of P. harmala. The results reflect the great diversity of P. harmala and data obtained from this study can be used for future collecting missions.