The use of multilocus sequence analysis (MLSA) for the taxonomy of Bradyrhizobium was assessed. We compared partial sequences for atpD, recA, gyrB, rpoB and dnaK for a set of reference strains representing named species and genospecies, and a number of new isolates from Lupinus albus, Arachis hypogaea and Ornithopus compressus from Spain. The phylogenies of the individual genes were compared with previous DNA-DNA hybridization results. High hybridization values were well reflected, but intermediary hybridization values were less clearly apparent. However, the phylogeny of a concatenated dataset of the five genes did reflect all values and thus is more informative of overall genome similarity. Our results indicate that only for the genes gyrB, rpoB and dnaK there is a small gap between the interspecies sequence similarities and the intraspecies similarity, and therefore cut-off levels for species delineation cannot be set, although high sequence similarity (>99%) does permit identification. In a few instances, a reference strain did not group as expected for one of the five genes tested. This may be a result of horizontal gene transfer and recombination events occasionally involving housekeeping genes. This observation indicates it is best to consider more than one gene for taxonomic inferences. The majority of the new isolates from the three host species was identified as Bradyrhizobium canariense. Four strains from L. albus from León, Spain, formed a separate group close to Bradyrhizobium japonicum.