Genomic insight into the taxonomy of Rhizobium genospecies that nodulate Phaseolus vulgaris

Due to the wide cultivation of bean (Phaseolus vulgaris L.), rhizobia associated with this plant have been isolated from many different geographical regions. In order to investigate the species diversity of bean rhizobia, comparative genome sequence analysis was performed in the present study for 69 Rhizobium strains mainly isolated from root nodules of bean and clover (Trifolium spp.). Based on genome average nucleotide identity, digital DNA:DNA hybridization, and phylogenetic analysis of 1,458 single-copy core genes, these strains were classified into 28 clusters, consistent with their species definition based on multilocus sequence analysis (MLSA) of atpD, glnII, and recA. The bean rhizobia were found in 16 defined species and nine putative novel species; in addition, 35 strains previously described as Rhizobium etli, Rhizobium phaseoli, Rhizobium vallis, Rhizobium gallicum, Rhizobium leguminosarum and Rhizobium spp. should be renamed. The phylogenetic patterns of symbiotic genes nodC and nifH were highly host-specific and inconsistent with the genomic phylogeny. Multiple symbiovars (sv.) within the Rhizobium species were found as a common feature: sv. phaseoli, sv. trifolii and sv. viciae in Rhizobium anhuiense; sv. phaseoli and sv. mimosae in Rhizobium sophoriradicis/R. etli/Rhizobium sp. III; sv. phaseoli and sv. trifolii in Rhizobium hidalgonense/Rhizobium acidisoli; sv. phaseoli and sv. viciae in R. leguminosarum/Rhizobium sp. IX; sv. trifolii and sv. viciae in Rhizobium laguerreae. Thus, genomic comparison revealed great species diversity in bean rhizobia, corrected the species definition of some previously misnamed strains, and demonstrated the MLSA a valuable and simple method for defining Rhizobium species.