We determined partial ND4 gene sequences of mitochondrial DNA from 15 heterorhabditid nematode isolates, representing 5 species collected from different regions of the world, by using polymerase chain reaction (PCR) and direct-sequencing of PCR products. Aligned nucleotide as well as amino acid sequences were used to differentiate nematode species by comparing sequence divergence and to infer phylogeny of the nematodes by using maximum parsimony and likelihood methods. Robustness of our phylogenetic trees was checked by bootstrap tests. The 15 nematode isolates can be divided into 7 haplotypes based on DNA sequences. On a larger scale, the sequence divergence revealed 4 distinct groups corresponding to 4 described species. No sequence divergence was detected from 5 isolates of Heterorhabditis bacteriophora or between Heterorhabditis marelatus to Heterorhabditis hepialius. Our sequence data yielded phylogenetic trees with identical topologies when different tree-building methods were used. Most relationships were also confirmed by using amino acid sequences in maximum parsimony analysis. Our molecular phylogeny of Heterorhabditis species support an existing taxonomy that is based largely on morphology and the sequence divergence of the ND4 gene permits species identification.