The present study examined sequence variation in four mitochondrial (mt) genes, namely cytochrome c oxidase subunits 1 (cox1) and 2 (cox2), and NADH dehydrogenase subunits 1 and 2 (nad1 and nad2) among Clonorchis sinensis isolates from different endemic regions in China, and their phylogenetic relationships with other zoonotic trematodes were reconstructed. A portion of the cox1 and cox2 genes (pcox1 and pcox2), and nad1 and nad2 genes (pnad1 and pnad2) were amplified separately from individual liver flukes by polymerase chain reaction (PCR) and the amplicons were subjected to sequencing from both directions. The intra-specific sequence variations within C. sinensis were 0–1.6% for pcox1, 0–1.4% for pcox2, 0–0.9% for pnad1 and 0–1.0% for pnad2. Phylogenetic analyses based on the combined sequences of pcox1, pcox2, pnad1 and pnad2 revealed that all the C. sinensis isolates grouped together and were closely related to Opisthorchis felineus. These findings revealed the existence of intra-specific variation in mitochondrial DNA (mtDNA) sequences among C. sinensis isolates from different geographic regions, and demonstrated that mtDNA sequences provide reliable genetic markers for phylogenetic studies of zoonotic trematodes.