Nosema ceranae is a widespread honeybee parasite, considered to be one of the pathogens involved in the colony losses phenomenon. To date, little is known about its intraspecific genetic variability. The few studies on N. ceranae variation have focused on the subunits of ribosomal DNA, which are not ideal for this purpose and have limited resolution. Here we characterized three single copy loci (Actin, Hsp70 and RPB1) in three N. ceranae isolates from Hungary and Hawaii. Our results provide evidence of unexpectedly high levels of intraspecific polymorphism, the coexistence of a wide variety of haplotypes within each bee colony, and the occurrence of genetic recombination in RPB1. Most haplotypes are not shared across isolates and derive from a few frequent haplotypes by a reduced number of singletons (mutations that appear usually just once in the sample), which suggest that they have a fairly recent origin. Overall, our data indicate that this pathogen has experienced a recent population expansion. The presence of multiple haplotypes within individual isolates could be explained by the existence of different strains of N. ceranae infecting honeybee colonies in the field which complicates, and must not be overlooked, further analysis of host–parasite interactions.