Both Hemibarbus labeo and H. medius (Cypriniformes: Cyprinidae: Gobioninae) are primary freshwater fishes and are widely distributed. As such, they provide an ideal model for phylogeographical studies. However, the similarity in morphological characters between these two species made the description of their distributions and the validation of species quite challenging. Here we employed variations in the DNA sequences of mitochondrial COI and ND5 genes (2151 bp) to solve this challenge and to study the population genetics structure of these two species. Among the 130 specimens belonging to 8 populations of H. labeo and 9 populations of H. medius from 17 drainage systems in southern China,196 variable sites (9.1% in the full sequences) falling into 50 haplotypes were identified. The haplotype diversity (h) and the nucleotide diversity (π) were 0.964 and 0.019, respectively, indicating a high level of genetic diversity and an evolutionary potential in both species. The result of neighbor-joining tree based on composite nucleotide sequences of the mtDNA COI and ND5 genes showed that the H. labeo and H. medius fell into two major clades (clade1and clade2): clade1was composed of some specimens of Oujiang River, all the specimens of Hanjiang River and Jiulongjiang River, whereas all remaining populations fell in clade2. The genetic distance between clade I and clade II was 0.036, while that between H. labeo and H. medius was 0.027. The haplotype network analyses indicated that the populations of Hanjiang River and Jiulongjiang River had relatively high genetic variation with the rest rivers. The po-pulations of Hainan Island migrated northward to Moyangjaing River. Haplotypes of the rivers of Hainan Island and Moyangjang River had relatively higher genetic variation with the Yangtze River than Pearl River. The populations of Xiangjiang River had no genetic variation with the populations of Guijiang River and Liujiang River. Analysis of molecular variance (AMOVA) indicated that the genetic variance mainly presented in individuals between geographical regions. The genetic variation of populations among regions was 71.2%, the genetic variation among populations within regions was 16.6%, and that within populations within the regions was 12.2%, indicating that most of the genetic variations resided in the populations among regions. The results of mismatch distribution and tests of neutrality suggested that in all populations, H. labeo, H. medius, clade1and clade2 were relatively stable.