The genus Harpagophytum has two species: H. procumbens which is an important medicinal plant in southern Africa, and H. zeyheri. Genetic diversity in 96 samples, obtained by germinating seeds collected from Botswana, was assessed using six inter-simple sequence repeat (ISSR) and 10 random amplified polymorphic DNA (RAPD) primers. These DNA markers yielded a total of 138 polymorphic bands. Polymorphism information content (PIC) ranged from 0.06 to 0.39 for ISSR primers, and from 0.09 to 0.43 for RAPD primers. Jaccard's similarity coefficients were highest when seedlings derived from the same fruit capsule were compared, while seedlings from different fruits on the same plant had intermediate values. The lowest values were recorded among seedlings from different plants. These results were consistent with an outcrossing breeding system in Harpagophytum. Analysis of molecular variance revealed significant differentiation (P < 0.01) between taxonomic units within Harpagophytum. About 39% of the variability occurred between the two species, H. procumbens and H. zeyheri. Plants with an intermediate morphology, i.e. putative hybrids (PH), showed 21% differentiation when compared with H. procumbens ssp. procumbens (PP), and 19% when compared with H. procumbens ssp. transvaalense (PT) or with H. zeyheri (ZZ). In addition, a deviating variant of PT was identified, here termed 'procumbens new variety' (PN). PN showed only 9% differentiation when compared with PT, 22% when compared with PP or with PH, and 41% when compared with ZZ. Considerable differentiation between the two Harpagophytum species was revealed also by a cluster analysis. Introgression was, however, suggested by the intermediate position of the putative hybrid plants in a principal component analysis while inter-specific gene flow was shown by a Bayesian genetic structure analysis.