Fagus japonica is generally known as one of two beech species native to Japan. It forms intermediate temperate forests that are vertically distributed from the lower part of the cool temperate zone into the upper part of the warm temperate zone, mostly along the Pacific side of Japan. We have studied the genetic diversity and population genetic structure in 16 F. japonica populations distributed throughout its range, using 13 nuclear microsatellite markers. Genetic diversity within populations was found to be high (average H E ¼ 0:659), whereas population differentiation was low (F ST ¼ 0:023; R ST ¼ 0:025). However, we detected significant isolation by distance. A neighbor-joining tree based on genetic distances among the populations reflected the geographic location of the populations and identified three distinct subdivisions, in the northeastern, central, and southwestern regions (the bootstrap values for internal branches between these regions are 87% and 91%). This finding was supported by boundaries detected by Monmonier's maximum-difference algorithm. In addition, Bayesian clustering analysis, using multilocus genotypes, revealed two distinct clusters that correlated with the geographical distributions of the populations, consistent with the results of the other analyses. The distribution of the two clusters suggests that, as a result of historical migrations, genetic mixing has occurred between two ancestral populations with different genetic backgrounds.