The distribution of Triplophysa yarkandensis is restricted to Xinjiang’s Tarim River basin. We collected 119 T. yarkandensis samples from nine geographic populations in the Tarim River basin and utilized the RAD-seq method for SNP genotyping. In this study, a total of 164.81 Gb bases were generated with the Illumina platform, and 129,873 candidate SNPs were obtained with the Stacks pipeline for population genetic analyses. High levels of genetic diversity were detected among nine populations. The AMOVA results showed that the majority of genetic variations originated from among populations (FST = 0.67), and the pairwise FST values ranged from 0.4579 to 0.8736, indicating high levels of genetic differentiation among these populations. The discriminate analysis of principal components (DAPCs) and neighbor joining (NJ) tree revealed that the nine populations could be separated into two clusters (i.e., south and north populations), and modest genetic differentiation between south and north populations was observed, while the individuals from several populations were not clustered together by geographical location. The evidence of two genetic boundaries between south and north populations (except TTM) was supported by barrier analysis. The Bayesian skyline plotting indicated that T. yarkandensis populations in the Tarim River basin had not experienced genetic bottlenecks, and the effective population size remained stable. This study first clarified the genetic diversity and differentiation of T. yarkandensis populations in the Tarim River basin, and it provided valuable molecular data for conservation and management of natural populations.