Thuja sutchuenensis Franch. is a critically endangered plant endemic to the North-East Chongqing, China. Genetic variation was studied to assess the distribution of genetic diversity within and among seven populations from the single remnant locations, using inter-simple sequence repeat (ISSR) markers. A total of 15 primers generated 310 well defined bands, with an average of 20.7 bands per primer. The seven populations revealed a relatively high level of genetic diversity in the species. The percentage of polymorphic bands, Nei's gene diversity and Shannon's information index at the population and species level were 76.1%, 0.155, 0.252 and 100%, 0.165, 0.295, respectively. A low level of genetic differentiation among populations (G(ST) = 0.102), in line with the results of Analyses of Molecular Variance (AMOVA), and a high level of gene flow (N(m) = 4.407) were observed. Both the Unweighted Pair Group Method with Arithmatic Mean (UPGMA) cluster analysis and Principal Coordinates Analysis (PCoA) supported the grouping of all seven populations into two groups. In addition, Mantel test revealed no significant correlation between genetic and geographical distances (r = 0.329, p = 0.100). The low genetic differentiation among populations implies that the conservation efforts should aim to preserve all the extant populations of this endangered species.