The genus Bradyrhizobium has been considered to be a taxonomically difficult group. In this study, phylogenetics and evolutionary genetics analyses were used to investigate divergence levels among Bradyrhizobium strains nodulating soybeans in China. Eleven genospecies were identified by sequence analysis of three phylogenetic and taxonomic markers (SMc00019, thrA, and truA). This was also supported by analyses of eight genes outside the symbiosis island ("off-island" genes; SMc00019, thrA, truA, fabB, glyA, phyR, exoN, and hsfA). However, seven genes inside the symbiosis island ("island" genes; nifA, nifH, nodC, nodV, fixA, trpD, and rhcC2) showed contrasting lower levels of nucleotide diversity and recombination rates than did off-island genes. Island genes had significantly incongruent gene phylogenies compared to the species tree. Four phylogenetic clusters were observed in island genes, and the epidemic cluster IV (harbored by Bradyrhizobium japonicum, Bradyrhizobium diazoefficiens, Bradyrhizobium huanghuaihaiense, Bradyrhizobium liaoningense, Bradyrhizobium daqingense, Bradyrhizobium sp. I, Bradyrhizobium sp. III, and Bradyrhizobium sp. IV) was not found in Bradyrhizobium yuanmingense, Bradyrhizobium sp. II, or Bradyrhizobium elkanii. The gene flow level of island genes among genospecies is discussed in the context of the divergence level of off-island genes.
Soybeans (Glycine max L.) were first domesticated in China and then introduced into different parts of the planet (1), now with an annual harvest area of 100 million hectares around the world (FAO, 2011). Ninety percent of their production comes from the United States, Brazil, Argentina, China, and India (FAO, 2007(FAO, to 2011. One of the key features of soybean is its ability to form symbiotic nitrogen-fixing nodules with diverse rhizobial species (2, 3), implying its important role in sustainable agriculture. It has been recurrently reported that Bradyrhizobium japonicum, Bradyrhizobium elkanii, Bradyrhizobium liaoningense, Bradyrhizobium yuanmingense, and Sinorhizobium fredii could nodulate soybeans (2-5). Recently, Bradyrhizobium huanghuaihaiense, Bradyrhizobium daqingense, Sinorhizobium sojae, and several unnamed species were also found to be effective microsymbionts of soybeans (2, 3, 6-8). Strain USDA110 represents a widely distributed type formerly known as B. japonicum Ia, but it has recently been proposed as a member of the new species Bradyrhizobium diazoefficiens (9).Recent studies not only suggested differences in the biogeographic distribution of rhizobial species nodulating soybeans but also demonstrated a biased selection of rhizobial species by different genotypes of soybeans (2,4,10,11). Consistent with these findings, comparative genomics of rhizobia revealed that the phyletic distribution of rhizobial functional genes involved in environmental adaptations and symbiotic interactions generally agrees with the phylogeny of rhizobial species (7). Therefore, it is important to distinguish different rhizobial species and even subdivis...