This study explores the genetic structure of Siberian indigenous populations on the basis of standard blood group and protein markers and DNA variable number of tandem repeats (VNTR) variation. Four analytical methods were utilized in this study: Harpending and Jenkin's R-matrix; Harpending and Ward's method of correlating genetic heterozygosity (H) to the distance from the centroid of the gene frequency array (r ii ); spatial autocorrelation, and Mantel tests. Because of the underlying assumptions of the various methods, the numbers of populations used in the analyses varied from 15 to 62. Since spatial autocorrelation is based upon separate correlations between alleles, a larger number of standard blood markers and populations were used. Fewest Siberian populations have been sampled for VNTRs, thus, only a limited comparison was possible. The four analytical procedures employed in this study yielded complementary results suggestive of the effects of unique historical events, evolutionary forces, and geography on the distribution of alleles in Siberian indigenous populations. The principal components analysis of the R-matrix demonstrated the presence of populational clusters that reflect their phylogenetic relationship. Mantel comparisons of matrices indicate that an intimate relationship exists between geography, languages, and genetics of Siberian populations. Spatial autocorrelation patterns reflect the isolation-by-distance model of Malecot and the possible effects of long-distance migration.