The present study examined genetic variability among Clonorchis sinensis isolates from four different geographical localities (Guangzhou, Nanning, Jiamusi and Daqing) and host species (cats, dogs, human and rabbits) in Mainland China by sequence analyses of two mitochondrial DNA (mtDNA) genes, namely NADH dehydrogenase subunits 2, 5 (nad2 and nad5) and ribosomal internal transcribed spacer 1 (ITS1). A portion of the ITS1, nad2 (pnad2) and nad5 (pnad5) was amplified by polymerase chain reaction separately from adult C. sinensis individuals and the amplicons were subjected to sequencing from both directions. The length of the sequences of ITS1, pnad2 and pnad5 was 643, 666 and 771 bp, respectively. The intraspecific sequence variations within C. sinensis were 0-1.7% for ITS1, 0-1.4% for pnad2 and 0-0.9% for pnad5. The interspecific sequence variations within other zoonotic trematodes, which were published previously, were 4.5-84.9% for ITS1, 21.9-43.6% for pnad2 and 19.2-48.9% for pnad5. The A+T contents of the sequences were 45.26-45.88% (ITS1), 62.91-63.51% (pnad2) and 58.24-58.63% (pnad5). Phylogenetic analyses using ribosomal and mitochondrial sequence data set, with three different computational algorithms (Bayesian inference, maximum parsimony and maximum likelihood), all revealed distinct groups with high statistical support. These findings demonstrated the existence of low-level intraspecific variations in ribosomal DNA (rDNA) and mtDNA sequences among C. sinensis isolates from four different regions and hosts in China and elucidated that mtDNA sequences and rDNA sequences provided reliable genetic markers for phylogenetic studies of zoonotic trematodes.