Phylogenetic relationships among six species of Epistylis (i. e. E. plicatilis, E. urceolata, E. chrysemydis, E. hentscheli, E. wenrichi, and E. galea) were investigated using sequences of the first internal transcribed spacer region (ITS-1) of ribosomal DNA (rDNA). Amplified rDNA fragment sequences consisted of 215 or 217 bases of the flanking 18S and 5.8S regions, and the entire ITS-1 region (from 145 to 155 bases). There were more than 33 variable bases between E. galea and the other five species in both the 18S region and the ITS-1 region. The affiliation of them was assessed using Neighbor-joining (NJ), maximum parsimony (MP) and maximum likelihood (ML) analyses. In all the NJ, MP and ML analyses E. galea, whose macronucleic position and shape are distinctly different from those of the other five species, was probably diverged from the ancestor of Epistylis earlier than the other five species. The topology in which E. plicatilis and E. hentscheli formed a strongly supported sister clade to E. urceolata, E. chrysemydis, and E. wenrichi was consistent with variations in the thickness of the peristomial lip. We concluded that the macronucleus and peristomial lip might be the important phylogenetic characteristics within the genus Epistylis.