Studies were performed to assess the genetic variation amongst isolates of the aphid-pathogenic fungus Pandora neoaphidis (syn. Erynia neoaphidis). 37 isolates were examined, from a range of pest and non-pest aphid species, as well as 21 from eight other entomophthoralean species. Universal primers were used to amplify the ITS rDNA regions and all of the species tested produced discrete ITS groups, with the exception of Conidiobolus spp. Neighbour-joining analysis of the ITS2 regions from P. neoaphidis, P. kondoiensis and Zoophthora radicans demonstrated that these three species formed distinct groups with sequence identities of 58–82% between the groups. An ITS size of ca 1100 bp was diagnostic for P. neoaphidis, while ca 1450 bp was characteristic of P. kondoiensis. ITS-RFLP analysis failed to yield intraspecific polymorphisms in any of the P. neoaphidis isolates screened, although it was useful in distinguishing between different entomophthoralean species. Some intraspecific variation in the ITS region was detected in a number of isolates of Z. radicans and Conidiobolus spp. We propose that two isolates previously identified as P. neoaphidis based on conidia morphology, are actually P. kondoiensis based on molecular studies. Sequencing analysis of the complete ITS region from P. neoaphidis and P. kondoiensis allowed species-specific primers to be developed for P. neoaphidis and P. kondoiensis. These were used to screen aphids infected in laboratory bioassays and from field-collected samples, without prior isolation of the fungus. The primers are useful tools for quantifying the epizootiology of P. neoaphidis in aphid populations, as well as assessing competitive interactions between these two species.