Multilocus Sequence Analysis for Assessment of the Biogeography and Evolutionary Genetics of Four Bradyrhizobium Species That Nodulate Soybeans on the Asiatic Continent

A highly supported maximum-likelihood species phylogeny for the genus Bradyrhizobium was inferred from a supermatrix obtained from the concatenation of partial atpD, recA, glnII, and rpoB sequences corresponding to 33 reference strains and 76 bradyrhizobia isolated from the nodules of Glycine max (soybean) trap plants inoculated with soil samples from Myanmar, India, Nepal, and Vietnam. The power of the multigene approach using multiple strains per species was evaluated in terms of overall tree resolution and phylogenetic congruence, representing a practical and portable option for bacterial molecular systematics. Potential pitfalls of the approach are highlighted. Seventy-five of the isolates could be classified as B. japonicum type Ia (USDA110/USDA122-like), B. liaoningense, B. yuanmingense, or B. elkanii, whereas one represented a novel Bradyrhizobium lineage. Most Nepalese B. japonicum Ia isolates belong to a highly epidemic clone closely related to strain USDA110. Significant phylogenetic evidence against the monophyly of the of B. japonicum I and Ia lineages was found. Analysis of their DNA polymorphisms revealed high population distances, significant genetic differentiation, and contrasting population genetic structures, suggesting that the strains in the Ia lineage are misclassified as B. japonicum. The DNA polymorphism patterns of all species conformed to the expectations of the neutral mutation and population equilibrium models and, excluding the B. japonicum Ia lineage, were consistent with intermediate recombination levels. All species displayed epidemic clones and had broad geographic and environmental distribution ranges, as revealed by mapping climate types and geographic origins of the isolates on the species tree.