A total of 41 Clostridium botulinum serotype E strains from different geographic regions, including Canada, Denmark, Finland, France, Greenland, Japan, and the United States, were compared by multilocus sequence typing (MLST), amplified fragment length polymorphism (AFLP) analysis, variable-number tandem-repeat (VNTR) analysis, and botulinum neurotoxin (bont) E gene sequencing. The strains, representing environmental, food-borne, and infant botulism samples collected from 1932 to 2007, were analyzed to compare serotype E strains from different geographic regions and types of botulism and to determine whether each of the strains contained the transposon-associated recombinase rarA, involved with bont/E insertion. MLST examination using 15 genes clustered the strains into several clades, with most members within a cluster sharing the same BoNT/E subtype (BoNT/E1, E2, E3, or E6). Sequencing of the bont/E gene identified two new variants (E7, E8) that showed regions of recombination with other E subtypes. The AFLP dendrogram clustered the 41 strains similarly to the MLST dendrogram. Strains that could not be differentiated by AFLP, MLST, or bont gene sequencing were further examined using three VNTR regions. Both intact and split rarA genes were amplified by PCR in each of the strains, and their identities were confirmed in 11 strains by amplicon sequencing. The findings suggest that (i) the C. botulinum serotype E strains result from the targeted insertion of the bont/E gene into genetically conserved bacteria and (ii) recombination events (not random mutations) within bont/E result in toxin variants or subtypes within strains.