The East Himalaya-Hengduan Mountains region is the center of diversity of the genus Primula, and P. sikkimensis is one of the most common members of the genus in the region. In this study, the genetic diversity and structure of P. sikkimensis populations in China were assessed using inter-simple sequence repeat (ISSR) and chloroplast microsatellite markers. The 254 individuals analyzed represented 13 populations. High levels of genetic diversity were revealed by ISSR markers. At the species level, the expected heterozygosity and Shannon's index were 0.4032 and 0.5576, respectively. AMOVA analysis showed that 50.3% of the total genetic diversity was partitioned among populations. Three pairs of chloroplast microsatellite primers tested yielded a total of 12 size variants and 15 chloroplast haplotypes. Strong cpDNA genetic differentiation (G (ST) = 0.697) and evidence for phylogeographic structure were detected (N (ST )= 0.788, significantly higher than G (ST)). Estimated rates of pollen-mediated gene flow are approximately 27% greater than estimated rates of seed-mediated gene flow in P. sikkimensis. Both seed and pollen dispersal, however, are limited, and gene flow among populations appears to be hindered by the patchiness of the species' habitats and their geographic isolation. These features may have played important roles in shaping the genetic structure of P. sikkimensis. A minimum-spanning tree of chloroplast DNA haplotypes was constructed, and possible glacial refugia of P. sikkimensis were identified.