# Unearthing Optimal Symbiotic Rhizobia Partners from the Main Production Area of Phaseolus vulgaris in Yunnan

Junjie Zhang, Jingqi Wang, Yufeng Feng, Brigitte Brunel, Xuxiao Zong

► To cite this version:

Junjie Zhang, Jingqi Wang, Yufeng Feng, Brigitte Brunel, Xuxiao Zong. Unearthing Optimal Symbiotic Rhizobia Partners from the Main Production Area of Phaseolus vulgaris in Yunnan. International Journal of Molecular Sciences, 2024, 25, pp.8511. 10.3390/ijms25158511. hal-04678104

HAL Id: hal-04678104

https://hal.inrae.fr/hal-04678104v1

Submitted on 26 Aug 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Unearthing Optimal Symbiotic Rhizobia Partners from the Main Production Area of *Phaseolus vulgaris* in Yunnan

Junjie Zhang<sup>1,2,*</sup>, Jingqi Wang<sup>1</sup>, Yufeng Feng<sup>1</sup>, Brigitte Brunel<sup>3</sup> and Xuxiao Zong<sup>4 ID</sup>

**Abstract:** *Phaseolus vulgaris* is a globally important legume cash crop, which can carry out symbiotic nitrogen fixation with rhizobia. The presence of suitable rhizobia in cultivating soils is crucial for legume cropping, especially in areas beyond the plant-host native range, where soils may lack efficient symbiotic partners. We analyzed the distribution patterns and traits of native rhizobia associated with *P. vulgaris* in soils of Yunnan, where the common bean experienced a recent expansion. A total of 608 rhizobial isolates were tracked from soils of fifteen sampling sites using two local varieties of *P. vulgaris*. The isolates were discriminated into 43 genotypes as defined by IGS PCR-RFLP. Multiple locus sequence analysis based on *recA*, *atpD* and *rpoB* of representative strains placed them into 11 rhizobial species of *Rhizobium* involving *Rhizobium sophorae*, *Rhizobium acidisoli*, *Rhizobium ecuadorense*, *Rhizobium hidalgonense*, *Rhizobium vallis*, *Rhizobium sophoriradicis*, *Rhizobium croatiense*, *Rhizobium anhuiense*, *Rhizobium phaseoli*, *Rhizobium chutanense* and *Rhizobium etli*, and five unknown *Rhizobium* species; *Rhizobium* genosp. I~V. R. phaseoli* and *R. anhuiense* were the dominant species (28.0% and 28.8%) most widely distributed, followed by *R. croatiense* (14.8%). The other rhizobial species were less numerous or site-specific. Phylogenies of *nodC* and *nifH* markers, were divided into two specific symbiovars, sv. phaseoli regardless of the species affiliation and sv. viciae associated with *R. vallis*. Through symbiotic effect assessment, all the tested strains nodulated both *P. vulgaris* varieties, often resulting with a significant greenness index (91–98%). However, about half of them exhibited better plant biomass performance, at least on one common bean variety, and two isolates (CYAH-6 and BLYH-15) showed a better symbiotic efficiency score. Representative strains revealed diverse abiotic stress tolerance to NaCl, acidity, alkalinity, temperature, drought and glyphosate. One strain efficient on both varieties and exhibiting stress abiotic tolerance (BLYH-15) belonged to *R. genosp*. IV sv. phaseoli, a species first found as a legume symbiont.

**Keywords:** common bean; *Rhizobium*; genetic diversity; symbiovar; biogeographic distribution; stress tolerance

## 1. Introduction

The common bean (*Phaseolus vulgaris* L.), known as green bean, is a leguminous plant of the genus *Phaseolus* [1,2]. *P. vulgaris* is one of the frequently consumed vegetables, native to Mexico and Argentina. *P. vulgaris* was first domesticated in the Americas, including Mexico, Colombia, Ecuador and northern Peru, as well as in the Andean center from southern Peru to northern Argentina [3]. *P. vulgaris* prefers warmth and is not tolerant to frost. It was only at the end of the 16th century that China began to introduce its cultivation. When *P. vulgaris* was introduced from the Americas, it was first introduced and domesticated in Guizhou, Yunnan, Sichuan and the surrounding provinces, before spreading to the northeast [4]. According to reports, beans are grown in more than 90 countries worldwide,

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Zhang, J.; Wang, J.; Feng, Y.; Brunel, B.; Zong, X. Unearthing Optimal Symbiotic Rhizobia Partners from the Main Production Area of *Phaseolus vulgaris* in Yunnan. *Int. J. Mol. Sci.* **2024**, *25*, 8511. https://doi.org/10.3390/ijms25158511

Academic Editor: Anna Gałązka

Received: 24 May 2024

Revised: 25 July 2024

Accepted: 31 July 2024

Published: 4 August 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

<footnote n="1">¹ College of Food and Bioengineering, Zhengzhou University of Light Industry, Zhengzhou 450000, China</footnote>

<footnote n="2">² Collaborative Innovation Center for Food Production and Safety of Henan Province,
Zhengzhou 450002, China</footnote>

<footnote n="3">³ Eco&Sols, University Montpellier, CIRAD, INRAE, Institut Agro, IRD, F-34398 Montpellier, France;
brigitte.brunel@supagro.fr</footnote>

<footnote n="4">⁴ Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China</footnote>

<footnote>* Correspondence: junjiezh@zzuli.edu.cn; Tel.: +86-371-8660-8670</footnote>
with a total planting area of 36.5 million km² and a total output of 31.4 million tons. This accounts for 50% of the total output of all edible beans, second only to the legume crop soybean. The planting area of *P. vulgaris* in China is 10 million km², and the total output is 1.3 million tons [5,6]. At present, *P. vulgaris* is widely cultivated in major producing areas such as northeast and north Xinjiang, and southwest China. Shandong Province and Hebei Province in north China, Heilongjiang Province and Jilin Province in northeast China and Yunnan, Guizhou and Sichuan provinces in southwest China have very rich varieties and types of beans. Yunnan is the province with the largest total planting area and output of edible beans in China, serving as the production base for high-quality fresh bean raw materials in the country. *P. vulgaris* has emerged as an important local cash crop because of its easy cultivation and low planting cost, resulting in a relatively extensive planting area. In some dam areas of Chuxiong, Dali, Baoshan and other places within Yunnan Province, the fresh bean industry is thriving. Cultivation of soft pod vine fresh beans has yielded substantial economic benefits, providing a new way for farmers in alpine mountainous areas to alleviate poverty and achieve prosperity. *P. vulgaris* is not only a vegetable, but also a grain or cash crop for foreign exchange export. Therefore, the vigorous development of bean production holds great significance for China's agriculture.

Nitrogen-fixing bacteria—known as rhizobia—are a group of Gram-negative bacteria characterized by their ability to fix nitrogen in symbiosis with legumes [7,8]. When rhizobia infect leguminous plants, they induce the formation of root nodules, and within root nodules, rhizobia ultimately differentiate in bacteroids that coexist with plant cells. Bacteroids draw nutrients from plants and obtain an environment suitable for nitrogen fixation, while plants obtain nitrogen nutrients through nitrogen fixation by bacteroids, which reduce atmospheric nitrogen to ammonia forms that are assimilable by the plants [9–13].

At present, there are several international studies on the diversity of rhizobia associated with *P. vulgaris*, and different rhizobial species forming symbiosis with *P. vulgaris* have been reported. These include *Rhizobium phaseoli* and *Rhizobium tropici* [14], *Rhizobium acidisoli* [15], *Rhizobium croatiense* and *Rhizobium redzepovicii* [16], *Rhizobium freirei* [17], *Rhizobium leguminosarum* [18], *Rhizobium etli* [19], *Rhizobium hidalgonense* [20], *Rhizobium ecuadorense* [21], *Rhizobium giardinii* and *Rhizobium gallicum* [22], *Rhizobium leucaenae* [23], *Rhizobium lusitanum*, *Pararhizobium giardinii* and *Rhizobium pisi* [24], *Rhizobium azibense* [25], *Rhizobium aethiopicum* [26] and *Rhizobium mesoamericanum* [27]. In addition, Chinese sinorhizobionodulating *P. vulgaris* have been found in species of *Sinorhizobium meliloti* [28], *Sinorhizobium americanum* [29] and *Sinorhizobium fredii* [30].

There have been several studies on the diversity of rhizobia associated with *P. vulgaris* in China, revealing various rhizobia populations and distribution. For instance, *Rhizobium vallis* strains were isolated from Yunnan leguminous plants in China [31] while *Rhizobium chutanense* strains were isolated from *P. vulgaris* in Jiangxi Province [32]. In addition, *P. vulgaris*-nodulating rhizobia resources from Heilongjiang, Liaoning and Jiangxi provinces in China were collected, and most of the isolated strains belonged to the genus *Rhizobium*. Among these, the species *R. leguminosarum*, *Rhizobium laguerreae*, *R. phaseoli* and *R. vallis* were identified, along with members of the genus *Bradyrhizobium* [33]. Furthermore, the genera *Sinorhizobium*, *Rhizobium*, *Bradyrhizobium*, *Ochrobactrum* and *Agrobacterium* were isolated from *P. vulgaris* in Shaanxi Province. These included species such as *R. phaseoli*, *R. vallis*, *R. giardinii*, *Rhizobium yanglingense*, *R. leguminosarum*, *Sinorhizobium adhaerens*, *S. fredii*, *Sinorhizobium kunmerowiae*, *Bradyrhizobium liaoningense*, *Ochrobactrum anthropic* and *Agrobacterium radiobacter*, totaling 11 species. These findings highlight the very rich diversity of rhizobia associated with *P. vulgaris* [34].

Across different other countries and continents, *P. vulgaris* can also coexist with a wide variety of rhizobia belonging to different genotypes, genera (*Rhizobium*, *Sinorhizobium* and *Bradyrhizobium*) and species with dominant strains or taxa [35]. In its center of origin, the dominant species among *P. vulgaris* nodules is *R. etli*. However, in some regions of Latin America, the prevalence of native strains (*R. leguminosarum*, *R. gallicum* and *R. giardinii*) can hinder the effectiveness of inoculation with *R. etli*, which is generally used to promote
nodulation and nitrogen fixation. In areas where pulses have been cultivated, legume nodules harbor many species other than *R. etli* [24]. *R. tropici*, characterized by high degree of genetic stability, dominates in acidic soils and regions with high-temperatures [36]. In Africa, *R. phaseoli*, *R. etli* and a new *Rhizobium* taxon had a great advantage in forming symbiotic relationships with *P. vulgaris*. Meanwhile, *R. leguminosarum* sv. phaseoli, isolated from Moroccan soil, was found more tolerant to acidic conditions in culture media or sterile soil [18].

With the sustainable development of agriculture in China, the frequent use of chemical fertilizers and pesticides has led to environmental contamination and declining soil fertility. Nitrogen is one of the most restricted nutrients for plant growth, and nitrogen fertilizers incur significant costs in crop production [37,38]. In agriculture, the symbiotic relationship between leguminous plants and rhizobia is the key to agricultural and environmental sustainability [39,40]. This relationship provides a natural and renewable nitrogen resource for crops, which is economical and environmentally friendly [41,42]. The rhizobia-leguminous symbiosis system, renowned for its strongest nitrogen fixation capacity, plays an important role in promoting the ecological restoration of contaminated land, including areas affected by saline conditions, heavy metals and pesticide pollution [43]. Yu and Liu [44] found that enhancing soybean salt tolerance could lead to improved soybean yield in saline soils, thus highlighting the potential for biological improvement. Han [45] collected germplasm resources of elite rhizobia from broad beans in the Qinghai area, and finally screened out the strains with high drought and salt-alkali tolerance, offering a good application prospect. Chi et al. [46] conducted NaCl and drought tolerance tests on 58 peanut-nodulating rhizobia isolated from different regions of Shandong Province, and found diversity in salt and drought tolerance among rhizobia not only throughout the province but also within a same region. Wu et al. [47] screened out excellent rhizobia from soybean nodules with strong resistance to high temperatures, salt, antibiotic, acidic and alkaline conditions, as well as strong resistance to dyes and chemical drugs. Cheng et al. [48] simulated cultivated land with a history of glyphosate use by spraying different concentrations of glyphosate solution on the soil before sowing alfalfa. The results showed that different concentrations of glyphosate had inhibitory effects on the growth and nitrogen fixation of alfalfa, with the inhibitory effect strengthening with increasing glyphosate concentration. Therefore, improving crop salt tolerance, acidity and alkalinity tolerance, high temperature tolerance, drought tolerance and other traits, as well as the comprehensive development of biological treatment of difficult and constrained soils, are major issues for future agricultural development. In this study, an optimal matching experiment between rhizobia and *P. vulgaris* was carried out to screen out rhizobia with high efficiency in nitrogen fixation and strong abiotic stress tolerance. The aim was to improve the yield of *P. vulgaris*, mitigate environmental contamination caused by chemical fertilizers and pesticides and improve the quality of *P. vulgaris* [49,50].

Considering all the aforementioned aspects, and the fact that common bean-nodulating rhizobia in Yunnan (China) have not been systematically studied, we conducted the present study. The aim of this work was to evaluate the diversity, relative abundance and geographic distribution of native rhizobia that nodulate *P. vulgaris* in Yunnan Province. Traditionally, common beans have been grown largely in Yunnan as an economic crop for fresh vegetables and grains, and the cultivation area has rapidly expanded to meet consumer demand. Thus, the taxonomic status of the isolated strains from root nodules was determined through ribosomal intergenic typing, phylogenetic analyses of housekeeping genes (*recA*, *atpD* and *rpoB*), the 16S rRNA gene and symbiotic marker genes (*nodC* and *nifH*). Additionally, the distribution of rhizobia in relation to soil properties and environmental factors, as well as the potential of representative strains to induce effective symbiosis and tolerate abiotic stress, were investigated.
## 2. Results

### 2.1. Physicochemical Characteristics of Soils and the Environment

All 15 sites differed in pH, as well as in their levels of organic matter (OM), alkaline hydrolyzable nitrogen (AN), available phosphorus (AP), available potassium and total salts (EC) (Supplementary Table S1). Soil samples ranged from acidic (five sites) to neutral (five sites) to slightly alkaline (five sites). Most of them contained an average OM level (2–4% across 11 sites), while two had a slightly low OM level (1.7–1.9%) and one a high level (>5%). The field soil at DLEY contained the highest contents of OM (62.6 g/kg soil) and AN (388.9 mg/kg soil). BSSD contained the highest content of AP (388.9 mg/kg soil) while soil from site CXMD contained the highest AK (358.3 mg/kg soil). At the opposite end, site BSCN exhibited the lowest proportions of AP and OM in soil. The field soil at CXDY had the highest salinity (EC = 819.3 µS/kg soil) and site DLMD had the lowest (69.1). The sites belong to the subtropical highland climate (Cwb, Kôppen classification) with mild temperatures and dry winters. Precipitations varied from moderate (nine sites with 870–985 mm/year) to high (six sites with 1020–1350 mm/year). The altitude of the sampling sites was relatively high (1033–2967 m).

### 2.2. Rhizobial Collection and IGS PCR-RFLP Analysis

In total, 608 rhizobial isolates were obtained from the 15 locations (varying from nine (CXWD) to fifty-two (CXDM and CXLF) isolates/site): 303 were obtained by host trapping with variety 1 (black seeds) and the remaining 305 isolates by variety 2 (white seeds) (Supplementary Table S2). The 608 isolates were distinguished into 25, 23 and 22 RFLP patterns obtained from the restriction enzymes *Hae*III, *Msp*I and *Hha*I, respectively. By combining all RFLP patterns, the IGS PCR-RFLP typing allowed the classification of all isolates into 43 IGS types (Table 1). Among the 43 IGS types, IGS type 1 represented the most abundant population (with 153 isolates), followed by IGS type 2 (63 isolates), type 3 (44 isolates), type 4 (38 isolates), type 5 (26 isolates) and type 6 (19 isolates), while the 37 remaining IGS types represented 265 isolates (Supplementary Table S2). Isolates of IGS type 1 were distributed over twelve sampling sites, with the CXLF and DLMD sites having the highest number (68 isolates, 11.2%). IGS type 2 isolates were detected at sites CXMD (nineteen isolates, 3.1%), CXLF (two isolates, 0.3%), CXDH (nine isolates, 1.5%), CXSB (thirteen isolates, 2.1%), CXNH (four isolates, 0.7%), DLXY (nine isolates, 1.5%), BSLY (four isolates, 0.7%) and BSSD (three isolates, 0.5%) (Supplementary Table S2). Sites CXLF, DLXY and CXNH were inhabited by the most diverse *P. vulgaris*-nodulating rhizobial community with thirteen to seventeen IGS types, while sites CXWD and BSCN both harbored only three to four IGS types, indicating that the richness and evenness of rhizobial IGS types varied across the sampling area (Supplementary Table S2).

### 2.3. Identification of Species by Phylogenetic Analysis of Core Genes

Nearly full-length 16S rRNA genes were successfully amplified and sequenced for 46 rhizobial isolates representing all 43 IGS types and sites of origin (Table 1). The representative isolates divided into three groups (Groups 1–3) in the phylogenetic tree (Supplementary Figure S1). First of all, fifteen representative isolates clustered together with several defined *Rhizobium* species in Group 1 that showed 99.7–100% similarity in their 16S rRNA gene sequences. This clade comprises the type strains of *Rhizobium acidisoli* FH13T, *R. anhuiense* CCBAU 23252T, *Rhizobium hidagonense* FH14T and *Rhizobium sophorae* CCBAU 03386T. Secondly, five representative isolates clustered in Group 2, which also included the type strains of *Rhizobium dioscoreae* S-93T, and *Rhizobium vallis* CCBAU 65647T. They shared 98.1–100% similarities with each other. Finally, Group 3 contained 26 representative isolates which shared 99.6–100% similarities with *Rhizobium phaseoli* ATCC 14482T, *Rhizobium ecuadorense* CNPSo 671T, *Rhizobium chutanense* C5T, *Rhizobium bangladeshense* BLR175T, *Rhizobium aethiopicum* HBR26T, *Rhizobium etli* CFN 42T, *Rhizobium sophoriradicis* CCBAU 03470T and *Rhizobium croatiense* 13T. Thus, all the representative isolates had been identified as belonging to *Rhizobium*.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Maximum likelihood phylogenetic tree based on concatenated *recA-atpD-rpoB* gene sequences (1185 base pairs) showing the relationships of rhizobia isolated from *Phaseolus vulgaris* L. in Yunnan Province of China. The tree was constructed under the best-fit model (GTR + G + I). Scale bar indicates 0.02 nt substitution per site. Bootstrap confidence values (%) calculated for 500 replications > 50% are indicated at the internodes.</figcaption></figure>
<table><caption>**Table 1.** Genetic groupings of *Rhizobium* isolates associated with *Phaseolus vulgaris* and their geographical distribution in the different sampling sites.</caption><thead><tr><th rowspan="2">IGS RFLP Type</th><th rowspan="2">Isolate Number</th><th rowspan="2">Representative Isolate WYCCWR No./Origin Site <sup>a</sup></th><th colspan="11">MLSA Similarity (%) with <sup>b</sup>:</th><th rowspan="2">Species Identification (Clade in Figure 1)</th></tr><tr><th>Rso</th><th>Rac</th><th>Rec</th><th>Rhi</th><th>Rva</th><th>Rsr</th><th>Rcr</th><th>Ran</th><th>Rph</th><th>Rch</th><th>Ret</th></tr></thead><tbody><tr><td>20</td><td>9</td><td>WYCCWR13672/CDYB-22</td><td>99.1</td><td>95.8</td><td>96.3</td><td>94.9</td><td>95.1</td><td>95.4</td><td>95.1</td><td>95.1</td><td>94.9</td><td>94.2</td><td>93.4</td><td>Rhizobium sophorae (C1)</td></tr><tr><td>18</td><td>9</td><td>WYCCWR13777/CNHH-10</td><td>95.2</td><td>98.6</td><td>96.3</td><td>95.5</td><td>95.3</td><td>95</td><td>94.9</td><td>95</td><td>95</td><td>95.7</td><td>94.3</td><td>Rhizobium acidisoli (C2)</td></tr><tr><td>13</td><td>11</td><td>WYCCWR14027/DEYH-16</td><td>95.9</td><td>96.8</td><td>99</td><td>95.7</td><td>95.1</td><td>94.6</td><td>95</td><td>94.9</td><td>94.7</td><td>95.1</td><td>93.7</td><td>Rhizobium ecuadorense (C3)</td></tr><tr><td>10</td><td>16</td><td>WYCCWR13936/DMDH-21</td><td>94.3</td><td>94.6</td><td>93.9</td><td>95.3</td><td>94</td><td>93.6</td><td>93.9</td><td>93.6</td><td>94.4</td><td>95.6</td><td>94.8</td><td rowspan="2">Rhizobium genosp. I (C4)</td></tr><tr><td>30</td><td>5</td><td>WYCCWR14133/BSDB-6</td><td>94.2</td><td>94.7</td><td>94</td><td>95.2</td><td>94</td><td>93.5</td><td>93.9</td><td>94</td><td>94.5</td><td>95.6</td><td>94.6</td></tr><tr><td>43</td><td>1</td><td>WYCCWR14114/BSDH-6</td><td>94.7</td><td>95.3</td><td>95.2</td><td>97.7</td><td>95</td><td>95.3</td><td>95.7</td><td>94.1</td><td>96.5</td><td>95.2</td><td>93.4</td><td rowspan="7">Rhizobium hidalgonense (C5)</td></tr><tr><td>25</td><td>6</td><td>WYCCWR13998/DWSB-10</td><td>94</td><td>95.2</td><td>95.1</td><td>98.5</td><td>94.5</td><td>94.1</td><td>93.9</td><td>95.7</td><td>94</td><td>94.7</td><td>93.3</td></tr><tr><td>42</td><td>1</td><td>WYCCWR13867/DXYH-4</td><td>94.5</td><td>95.6</td><td>95.7</td><td>99.8</td><td>95.1</td><td>94.4</td><td>94.3</td><td>94.3</td><td>94.2</td><td>94.9</td><td>93.8</td></tr><tr><td>37</td><td>3</td><td>WYCCWR13972/DWSH-6</td><td>94.6</td><td>95.7</td><td>95.8</td><td>99.8</td><td>95.2</td><td>94.6</td><td>94.5</td><td>94.5</td><td>94.4</td><td>95.1</td><td>93.6</td></tr><tr><td>41</td><td>1</td><td>WYCCWR13551/CMDH-27</td><td>94.6</td><td>95.7</td><td>95.8</td><td>100</td><td>95.2</td><td>94.6</td><td>94.5</td><td>94.5</td><td>94.4</td><td>95.1</td><td>93.6</td></tr><tr><td>8</td><td>16</td><td>WYCCWR13528/CMDH-4</td><td>94.6</td><td>95.7</td><td>95.8</td><td>100</td><td>95.2</td><td>94.6</td><td>94.5</td><td>94.5</td><td>94.4</td><td>95.1</td><td>93.6</td></tr><tr><td>32</td><td>4</td><td>WYCCWR13680/CLFH-5</td><td>94.6</td><td>95.7</td><td>95.8</td><td>100</td><td>95.2</td><td>94.6</td><td>94.5</td><td>94.5</td><td>94.4</td><td>95.1</td><td>93.6</td></tr><tr><td>39</td><td>2</td><td>WYCCWR13741/CDHH-14</td><td>95.1</td><td>95.1</td><td>94.6</td><td>95</td><td>98.9</td><td>95.1</td><td>96.3</td><td>94.4</td><td>95.1</td><td>95.5</td><td>93.4</td><td rowspan="2">Rhizobium vallis (C6)</td></tr><tr><td>12</td><td>12</td><td>WYCCWR13728/CDHH-1</td><td>95.1</td><td>95.2</td><td>94.9</td><td>95</td><td>99.3</td><td>95.1</td><td>96.3</td><td>94.3</td><td>95.1</td><td>95.5</td><td>93.4</td></tr><tr><td>28</td><td>5</td><td>WYCCWR13814/CWDB-3</td><td>95.6</td><td>95.4</td><td>95</td><td>95.1</td><td>95.3</td><td>98.5</td><td>96.3</td><td>94.9</td><td>95.2</td><td>95.4</td><td>93.5</td><td>Rhizobium sophoriradicis (C7)</td></tr><tr><td>23</td><td>8</td><td>WYCCWR13632/CDYH-7</td><td>95.1</td><td>95.4</td><td>94.6</td><td>93.9</td><td>96.2</td><td>96.8</td><td>98.6</td><td>94</td><td>95.8</td><td>94.8</td><td>93.9</td><td rowspan="5">Rhizobium croatiense (C8)</td></tr><tr><td>3</td><td>44</td><td>WYCCWR14106/BCNB-10</td><td>95.1</td><td>95.4</td><td>94.6</td><td>93.9</td><td>96.2</td><td>96.8</td><td>98.6</td><td>94</td><td>95.8</td><td>94.8</td><td>93.9</td></tr><tr><td>15</td><td>10</td><td>WYCCWR13582/CYAH-6</td><td>95.1</td><td>95.3</td><td>94.6</td><td>93.9</td><td>96.2</td><td>96.6</td><td>98.4</td><td>94</td><td>95.8</td><td>94.8</td><td>93.9</td></tr><tr><td>6</td><td>19</td><td>WYCCWR13581/CYAH-5</td><td>95.1</td><td>95.1</td><td>94.6</td><td>93.8</td><td>96.2</td><td>96.9</td><td>98.6</td><td>94</td><td>95.9</td><td>94.6</td><td>93.9</td></tr><tr><td>21</td><td>9</td><td>WYCCWR13601/CYAH-25</td><td>95.2</td><td>95.4</td><td>94.7</td><td>94</td><td>96.3</td><td>96.7</td><td>98.5</td><td>94</td><td>95.9</td><td>94.7</td><td>94</td></tr><tr><td></td><td></td><td>WYCCWR13663/CDYB-13</td><td>95.1</td><td>95.3</td><td>94.6</td><td>93.9</td><td>96.2</td><td>96.6</td><td>98.4</td><td>94</td><td>95.8</td><td>94.8</td><td>93.9</td></tr><tr><td>1</td><td>153</td><td>WYCCWR13630/CDYH-5</td><td>95.1</td><td>95.8</td><td>95.1</td><td>95.1</td><td>94.9</td><td>95.3</td><td>95.1</td><td>98.6</td><td>95.1</td><td>95</td><td>94.3</td><td rowspan="2">Rhizobium anhuiense (C9)</td></tr><tr><td></td><td></td><td>WYCCWR13646/CDYH-21</td><td>95.1</td><td>95.7</td><td>95.1</td><td>94.7</td><td>94.6</td><td>95.1</td><td>94.7</td><td>98.9</td><td>94.7</td><td>94.6</td><td>94.1</td></tr><tr><td>29</td><td>5</td><td>WYCCWR13999/DWSB-11</td><td>94.2</td><td>95.1</td><td>94.9</td><td>96.3</td><td>94.8</td><td>94.5</td><td>94</td><td>98.2</td><td>94.5</td><td>94.7</td><td>93.4</td><td rowspan="3">Rhizobium anhuiense (C9)</td></tr><tr><td>11</td><td>13</td><td>WYCCWR13604/CYAB-3</td><td>94.9</td><td>95.3</td><td>94.9</td><td>94.4</td><td>94.5</td><td>94.2</td><td>94.1</td><td>99.7</td><td>94.5</td><td>94.4</td><td>93.7</td></tr><tr><td>34</td><td>4</td><td>WYCCWR13828/CSBH-8</td><td>95</td><td>95.3</td><td>95</td><td>94.5</td><td>94.5</td><td>94.6</td><td>94.3</td><td>100</td><td>94.6</td><td>94.6</td><td>93.9</td></tr><tr><td>14</td><td>10</td><td>WYCCWR13549/CMDH-25</td><td>94.6</td><td>95.6</td><td>95.3</td><td>97</td><td>95.2</td><td>94.6</td><td>94.6</td><td>95.1</td><td>96.3</td><td>95.8</td><td>94.6</td><td>Rhizobium genosp. II (C10)</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th rowspan="2">IGS RFLP Type</th><th rowspan="2">Isolate Number</th><th rowspan="2">Representative Isolate WYCCWR No./Origin Site <sup>a</sup></th><th colspan="11">MLSA Similarity (%) with <sup>b</sup>:</th><th rowspan="2">Species Identification (Clade in Figure 1)</th></tr><tr><th>Rso</th><th>Rac</th><th>Rec</th><th>Rhi</th><th>Rva</th><th>Rsr</th><th>Rcr</th><th>Ran</th><th>Rph</th><th>Rch</th><th>Ret</th></tr></thead><tbody><tr><td>2</td><td>63</td><td>WYCCWR13564/CMDB-13</td><td>94.6</td><td>95.1</td><td>94.7</td><td>94.7</td><td>95</td><td>95.3</td><td>95.8</td><td>94.7</td><td>98.4</td><td>96</td><td>94.4</td><td rowspan="11">Rhizobium phaseoli (C11)</td></tr><tr><td>40</td><td>2</td><td>WYCCWR13810/CNHB-19</td><td>94.5</td><td>95.1</td><td>94.6</td><td>94.6</td><td>95</td><td>95.1</td><td>95.8</td><td>94.7</td><td>98.2</td><td>96</td><td>94.2</td></tr><tr><td>4</td><td>38</td><td>WYCCWR13550/CMDH-26</td><td>94.5</td><td>94.7</td><td>94.2</td><td>93.9</td><td>95</td><td>95.5</td><td>96.3</td><td>94.6</td><td>98.9</td><td>95.3</td><td>94.2</td></tr><tr><td>38</td><td>3</td><td>WYCCWR13893/DXYB-3</td><td>94.5</td><td>94.7</td><td>94.2</td><td>93.9</td><td>95</td><td>95.5</td><td>96.3</td><td>94.6</td><td>98.9</td><td>95.3</td><td>94.2</td></tr><tr><td>22</td><td>9</td><td>WYCCWR14123/BSDH-15</td><td>94.7</td><td>95.1</td><td>94.5</td><td>94.1</td><td>95.2</td><td>95.7</td><td>96.6</td><td>94.6</td><td>99.5</td><td>95.4</td><td>94.4</td></tr><tr><td>36</td><td>3</td><td>WYCCWR13915/DXYB-25</td><td>94.6</td><td>95.1</td><td>94.7</td><td>94.7</td><td>95</td><td>95.3</td><td>95.8</td><td>94.7</td><td>98.4</td><td>96</td><td>94.4</td></tr><tr><td>7</td><td>17</td><td>WYCCWR13909/DXYB-19</td><td>94.6</td><td>94.7</td><td>94.4</td><td>94</td><td>95</td><td>95.7</td><td>96.3</td><td>94.6</td><td>99.1</td><td>95.3</td><td>94.4</td></tr><tr><td>16</td><td>10</td><td>WYCCWR13839/CSBH-19</td><td>94.6</td><td>95.1</td><td>94.7</td><td>94.7</td><td>95</td><td>95.3</td><td>95.8</td><td>94.7</td><td>98.4</td><td>96</td><td>94.4</td></tr><tr><td>33</td><td>4</td><td>WYCCWR13823/CSBH-3</td><td>94.6</td><td>94.7</td><td>94.4</td><td>94</td><td>95</td><td>95.7</td><td>96.3</td><td>94.6</td><td>99.1</td><td>95.3</td><td>94.4</td></tr><tr><td>26</td><td>5</td><td>WYCCWR13529/CMDH-5</td><td>94.6</td><td>95.1</td><td>94.7</td><td>94.7</td><td>95</td><td>95.3</td><td>95.8</td><td>94.7</td><td>98.4</td><td>96</td><td>94.4</td></tr><tr><td>9</td><td>16</td><td>WYCCWR13697/CLFH-22</td><td>94.6</td><td>95.1</td><td>94.7</td><td>94.7</td><td>95</td><td>95.3</td><td>95.8</td><td>94.7</td><td>98.4</td><td>96</td><td>94.4</td></tr><tr><td>24</td><td>6</td><td>WYCCWR13654/CDYB-4</td><td>93.9</td><td>94</td><td>94</td><td>94.3</td><td>95.1</td><td>94.2</td><td>94.5</td><td>94.1</td><td>95</td><td>96.3</td><td>95.1</td><td>Rhizobium genosp. III (C12)</td></tr><tr><td>27</td><td>5</td><td>WYCCWR14006/DWSB-18</td><td>94</td><td>95.4</td><td>94.6</td><td>95.1</td><td>95.8</td><td>94.3</td><td>94.3</td><td>94</td><td>94.7</td><td>98.6</td><td>94.7</td><td rowspan="3">Rhizobium chutanense (C13)</td></tr><tr><td>35</td><td>3</td><td>WYCCWR13808/CNHB-17</td><td>94</td><td>95.5</td><td>94.6</td><td>95.1</td><td>95.3</td><td>94.2</td><td>94.2</td><td>93.9</td><td>95</td><td>98.6</td><td>94.5</td></tr><tr><td>5</td><td>26</td><td>WYCCWR13754/CDHB-7</td><td>94.5</td><td>95.9</td><td>95.1</td><td>95.1</td><td>95.2</td><td>94.6</td><td>94.6</td><td>94.6</td><td>95.1</td><td>98.6</td><td>94.8</td></tr><tr><td>31</td><td>4</td><td>WYCCWR13563/CMDB-12</td><td>93</td><td>93.9</td><td>93</td><td>92.8</td><td>93.4</td><td>93.1</td><td>93.8</td><td>92.8</td><td>94.3</td><td>94.7</td><td>97.1</td><td>Rhizobium etli (C14)</td></tr><tr><td>17</td><td>9</td><td>WYCCWR14090/BLYB-15</td><td>93.3</td><td>93.3</td><td>92.9</td><td>92.4</td><td>93.4</td><td>94.5</td><td>93.9</td><td>92.3</td><td>94.5</td><td>94.6</td><td>94.2</td><td>Rhizobium genosp. IV (C15)</td></tr><tr><td>19</td><td>9</td><td>WYCCWR14071/BLYH-17</td><td>91.8</td><td>91.5</td><td>91.7</td><td>91.1</td><td>91.7</td><td>92.3</td><td>91.9</td><td>91.1</td><td>92.9</td><td>92.9</td><td>92.6</td><td>Rhizobium genosp. V (C16)</td></tr><tr><td>IGS type number: 43</td><td>Total isolate number: 608</td><td>Total representatives: 46</td><td colspan="11">Total species: 16</td></tr></tbody></table>

a. The collection center number of (WYCCWR)/original isolated number; b. Similarities in MLSA associated with the following type strains-Rso: Rhizobium sophorae; Rac: Rhizobium acidisoli; Rec: Rhizobium ecuadorense; Rhi: Rhizobium hidalgonense; Rva: Rhizobium vallis; Rsr: Rhizobium sophoriradicis; Rcr: Rhizobium croatiense; Ran: Rhizobium anthuiense; Rph: Rhizobium phaseoli; Rch: Rhizobium chutanense; Ret: Rhizobium etli.
The representative isolates were divided into thirteen clades (C1–C16) in the phylogenetic tree based on their concatenated *recA-atpD-rpoB* sequences (Figure 1). C1 includes nine isolates from IGS type 20 with the type strain of *R. sophorae*, sharing 99.1% similarities (Table 1). Thus, cluster C1 was identified as *R. sophorae*. Cluster C2, identified as *R. acidisoli*, contained *R. acidisoli* FH13T and IGS type 18 representing nine isolates with 98.6% similarities. C3, identified as *R. ecuadorense*, contained 11 isolates of IGS type 13, showing 99% similarities with *R. ecuadorense* CNPSo 671T. C5 was identified as *R. hidalgonense* and comprised the type strain of *R. hidalgonense* with IGS types 8, 25, 32, 37, 41, 42 and 43 (representing together 32 isolates (5.3% in total) sharing 97.7–100% similarities with each other). Cluster C6, identified as *R. vallis*, contained 14 isolates from IGS types 12 and 39, showing 98.9–99.3% similarities. C7 comprised only the IGS type 28, covering five isolates, and was identified as belonging to the species *R. sophoriradicis* with 98.5% similarity to type strain CCBAU 03470 T. C8, identified as *R. croatiense*, contained *R. croatiense* 13T and IGS types 3, 6, 15, 21 and 23, totaling 90 isolates (14.8%) with 98.4–98.6% similarities. C9 included four representative isolates from IGS types 1, 11, 29 and 34 (total of 175 isolates, 28.8% in total) identified as *R. anhuiense* with 98.2–100% similarities. C11, identified as *R. phaseoli*, contained *R. phaseoli* ATCC 14482T and IGS types 2, 4, 7, 9, 16, 22, 26, 33, 36, 38 and 40, representing 170 isolates (28.0% in total) with 98.2–99.5% similarities. C13 included three representative isolates from IGS types 5, 27 and 35, identified as *R. chutanense* with 98.6% similarities. C14, identified as *R. etli*, contained four isolates from IGS type 31, showing 97.1% similarities with *R. etli* CFN 42 T. In addition, C4, 10, 12, 15 and 16 (Clade 4, 10, 12, 15, 16) contained six representative strains, and the similarity with the sequences of all known population strains was less than 97%, the model strain with the highest similarity in C4, C12 and C15 was *R. chutanense*, and the similarity was only 94.6–96.3%, and the model bacteria with the highest similarity between C10 and C16 was *R. phaseoli*, with a similarity of only 92.9–96.3%, is suspected to be a new population of rhizobia strains of such groups C4, 10, 12, 15 and 16. The phylogenetic analyses of the single gene of *recA*, *atpD* or *rpoB* are shown in Supplementary Figures S2–S4.

Group 1 isolates in the 16S rRNA gene phylogenetic tree (Figure S1) contained strains belonging to clades 1, 2, 4 7 and 13 (15 strains). Group 2 (five strains) included strains from clades 7 and 13 along with a strain (CDYB-4) from clade 10. Group 3 (26 strains) encompassed strains from clades 3, 6, 11 and 12, as well as clade 10 except for the strain CDYB-4. Apart from CDYB-4, the phylogenetic 16S rRNA gene tree was consistent with the MLSA classification in clades (Figure 1). Specifically, the CDYB-4 sequence exhibited low bootstrap support in clade 10 (<50%) compared to the other C10 sequences (98%), forming a single branch outside of the defined clades.

Populations of *R. anhuiense* were identified in 13 out of the 15 sampling sites (28.8% in total), with the highest number of isolates, representing the dominant species, observed in CXDY (59%), CXLF (71%), DLWS (53%), DLMD (66%) and DLEY (55%) (Figure 2). *R. phaseoli* was also found dominant (28.0%) and recovered in 12 sampling sites, with populations of CXMD (67%), CXSB (67%), DLXY (66%), CXNH (50%) and CXDH (35%) being the most abundant, establishing it as the dominant species in these five sites. Both *R. hidalgonense* (5.3%) and *R. croatiense* (14.8%) were recovered in 15 sampling sites, with dominance observed at CXYA (73%), CXWD (33%), BSCN (63%) and BSLY (39%) respectively. *R. chutanense* was detected in seven sampling sites, with CXDH (20%) and CXNH (16%) hosting the highest proportion of its members. *R. sophorae*, *R. acidisoli*, *R. ecuadorense*, *R. vallis*, *R. sophoriradicis*, *R. etli* and *Rhizobium* genosp. I~V were found in minority (<2%) in a limited number of sites (one to five). No differences were detected between rhizobia trapped with black beans of variety 1 compared to white beans of variety 2 within clades and sampling sites (Supplementary Figure S5, Table S2), indicating that the population distribution regardless of the common bean variety used for trapping was similar.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Distribution of core gene clades identifying species of rhizobia isolated from *Phaseolus vulgaris* L. per site and across all sites in Yunnan Province of China. The number of isolates is indicated under the site code. CXSB, Chuxiong Shuang bai, CXWD: Chuxiong Wu ding, CXYA: Chuxiong Yao an, CXMD: Chuxiong Mu ding, CXDH: Chuxiong Dong hua, CXDY: Chuxiong Dayao, CXLF: Chuxiong Lufeng, CXNH: Chuxiong Nanhua, DLEY: Dali Er yuan, DLXY: Dali Xiang yun, DLWS: Dali Wei shan, DLMD: Dali Mi du, BSLY: Baoshan Long yang, BSCN: Baoshan Chang ning, BSSD: Baoshan Shi dian.</figcaption></figure>

## 2.4. Identification of Symbiovars by Phylogenetic Analysis of Symbiotic Genes

Besides the core genes studied above to taxonomically identify rhizobia, it is complementary to type rhizobia according to their symbiovar related to the legume-host spectrum, which does not follow taxonomic phylogeny. To do so, common genes involved in nodulation (*nodC* coding for an N-acetyltransferase) and fixation (*nifH* coding for a subunit of nitrogenase) were investigated in the 46 representative strains.

The *nodC* genes were amplified from all the 46 representative isolates, confirming their genetic basis for rhizobial symbiosis. Phylogenetic analysis defined two strongly supported *nodC* groups (N1-N2) among the representatives (Figure 3). First, the symbiotype N1 (44 representatives) clustered together with *nodC* sequences from strains of several species, all from sv. phaseoli (such as *R. fabae* NC1, *R. croatiense* 13T<sup>T</sup>, *R. sophorae* CCBAU03386, *R. acidisoli* FH13<sup>T</sup>, *R. sophoriradicis* CCBAU 03470<sup>T</sup>, *R. chutanense* C16, *R. phaseoli* ATCC 14482<sup>T</sup>, *R. hidalgonense* FH14<sup>T</sup>, *R. etli* CFN L9<sup>T</sup> and *R. leguminosarum* bv. phaseoli LCS0313, showing a similarity of 97.5–100%. Second, symbiotype N2 (strains CDHH-1 and CDHH-14) encompassed *nodC* sequences 96% similar to the *nodC* of *R. ruizarguesonis* UMP1133<sup>T</sup>. Thus, all representatives belonged to the sv. phaseoli or sv. viciae (Figure 3).
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Maximum likelihood phylogenetic tree based on symbiotic gene *nodC* (376 base pairs) showing the relationships of the rhizobia isolated from nodules of *Phaseolus vulgaris* L. in Yunnan Province of China. The two *nodC* groups found among isolates are named N1 and N2. The tree was constructed using the maximum likelihood method under the best-fit model (T92 + G + I). Scale bar indicates 0.05 nt per site. Bootstrap confidence values (%) calculated for 500 replications > 70% are indicated at the internodes.</figcaption></figure>
Phylogenetic analysis of *nifH* symbiotic gene sequences divided the 46 studied representative strains of *Phaseolus* rhizobia into two distinct groups (Supplementary Figure S6). The larger group clustered 44 representatives with strains *R. chutanense* C5T, *R. aethiopicum* HBR26T, *R. ecuadorense* CNPSo 671T, *R. phaseoli* ATCC 14482T, *R. acidisoli* FH23T, *R. sophoriradicis* CCBAU 03470T and *R. croatiense* 13T<sup>T</sup>, sharing a sequence similarity of 95.3–100%. The smaller group consisted of two representative strains (CDHH-1 and CDHH-14) associated with *R. lentis* BLR27<sup>T</sup>. These two groups correspond to the *nodC* groups N1 (sv. phaseoli) and N2 (sv. viciae), respectively (Figure 3).

## 2.5. Correlation Analysis of Rhizobial Distribution with Soil and Environmental Properties

PCA was used to explore the relationships between soil and environment properties, and the rhizobial community composition based on IGS genotypes. The PCA results (Figure 4) showed that the soil chemical factors had different effects on the distribution of the rhizobia populations and IGS types. The IGS types 3, 6, 15, 17, 19, 20, 23, 30 and 43 (left upper part of Figure 4) included 114 isolates recovered mainly in CXDY, CXYA, BSSD, BSCN and BSLY (Supplementary Table S2). Their distribution was positively correlated with pH, EC and average rainfall (AvePrecp), and negatively related with soil OM contents. In particular, IGS20 and IGS15 showed a positive correlation with pH and EC, while IGS3 and IGS20 were associated with the average rainfall (Figure 4). The IGS types 1, 7, 11, 22, 25, 27, 29, 32, 34, 37 and 38 (lower left part of Figure 4) gathered the majority of isolates (222) distributed across four main sites (CXLF, DLEY, DLMD and DLWS) (Supplementary Table S2); they presented a positive association with AN, pH and OM, but were negatively correlated with AK values. In particular, the IGS27 and IGS29 were positively correlated with AN, while IGS37 and IGS7 correlated with higher OM values (Figure 4). The IGS types 5, 8, 10, 13, 24, 35, 36 and 42 (including 82 strains) (lower right part of Figure 4) tended to be associated with DLXY (Supplementary Table S2). Meanwhile, these strains were positively correlated with OM and negatively correlated with EC, pH value and average rainfall. Notably, the IGS10 was positively correlated with OM values (Figure 4). Finally, the IGS types 2, 4, 21, 12, 9, 31, 39, 18, 28, 33, 16, 14, 26, 40 and 41 (middle and upper right part of Figure 4), representing 190 isolates mainly from CXNH, CXDH, CXSB, CXMD and CXWD, were positively distributed with AK, but lower AN, EC and pH values. In particular, IGS28 was positively correlated with AK values. Thus, we found that the different rhizobial populations were influenced differently by edaphic factors, showing potential positive or negative significant correlations.

*R. sophoriradicis*, *R. acidisoli*, *R. vallis*, *R. etli* and *R. genosp*. II were mainly distributed in CXWD, CXSB, CXDH, CXNH and CXMD. Their distribution was positively correlated with the available potassium content in the soil and negatively correlated with electrical conductivity and pH value. *R. phaseoli*, *R. chutanense*, *R. hidalgonense*, *R. ecuadorense*, *R. genosp*. I and *R. genosp*. III were primarily distributed in the sampling points of DLXY, DLWS, DLMD and DLEY, with a positive correlation with organic matter content and negatively correlated with average precipitation. *R. sophorae* and *R. anhuiense* were mainly distributed in CXDY, CXLF and BSSD, with a positive correlation with alkaline nitrogen, electrical conductivity, and pH value, and a negative correlation with available potassium content. *R. croatiense*, *R. genosp*. IV and *R. genosp*. V were primarily distributed in CXYA, BSCN and BSLY. Their distribution was positively correlated with average precipitation and negatively correlated with the content of organic matter and alkaline nitrogen in the soil. Furthermore, *R. anhuiense* was widely distributed in all 13 sampling points, indicating that these strains were with strong adaptability to different conditions(Supplementary Figure S10).
<figure><img src="image_8.png" /><figcaption>**Figure 4.** PCA to relate the distribution of the 42 IGS types of isolates (n > 1) to physicochemical factors of soils and environment collected from the different sites. The blue arrows indicate IGS types of rhizobia, green indicate the sampling sites and red arrows represent soil properties and environmental factors. The longer the arrow was, the greater the influence of the soil property and environmental factor presents on the distribution of the IGS types. The smaller the angle between the arrow and the IGS type was, the stronger the effect of the soil property or environmental factor on distribution of the IGS type.</figcaption></figure>

## 2.6. Symbiotic Efficiency of Rhizobial Representative Strains

All representative isolates from the 42 different IGS types successfully nodulated *P. vulgaris*, with each plant forming more than 50 nodules, and the control group (CK) treated without bacterial treatment did not have nodulation. However, significant variations in symbiotic efficiency, were observed among the nodulated plants (Supplementary Figure S7). The chlorophyll content of *P. vulgaris* plants in the inoculated group was increased (Supplementary Figure S8), showing a significant difference compared to CK, except for strains CDHH-1 and DXYB-19 associated with variety v1 and 4 strains with v2 (CMDH-4, CSBH-8, BSDH-6 and BSDB-6). The average chlorophyll content per black *P. vulgaris* plant of the inoculated strains CDYB-4, CDYH-5, DWSB-18, CDYH-7 and CLFH-5 was among the highest with more than 1.33 times higher than that of CK. Similarly, the chlorophyll content of white *P. vulgaris* inoculated with strains CMDH-27, CYAH-6, BLYH-17, CDYB-13 and DMDH-21 were among the highest. In particular, the average chlorophyll content per plant of CDYB-13 and CNHB-17 was over 1.38 times higher than that of CK. The total mass of black *P. vulgaris* plants inoculated with rhizobia of *P. vulgaris*
was significantly improved (Supplementary Figure S9) compared to the control (CK) for 11 strains (CMDH-4, CYAH-6, CDYH-7, CDHB-7, CWDB-3, DXYH-4, DWSB-10, DWSB-11, DWSB-18, DEYH-16 and BLYB-15). Notably, DWSB-18, BLYB-15, CDYH-7, DWSB-11, CMDH-4, DEYH-16, CYAH-6 and DWSB-10 showed dry biomass more than 2.8 times that of CK. Similarly, the total dry weight mass of all white *P. vulgaris* plants inoculated with *P. vulgaris* rhizobia was significantly improved, particularly with strains BLYH-17, CYAH-6, CMDH-25, CDHB-7 and CYAH-25 exceeding 2.5 times that of CK.

The representative strains were evaluated based on the combined score of the three symbiotic indices mentioned above: chlorophyll (25%) + number of nodules (25%) + total plant dry weight (50%) (Figure 5). Among these, five strains were identified among the high-efficiency nitrogen-fixing group suitable for variety 1 (CMDH-4, CYAH-6, CDYH-7, DWSB-18 and BLYB-15), while another five strains could be designated among the high-efficiency nitrogen-fixing group suitable for variety 2 (CYAH-6, CMDH-27, CMDB-12, BLYH-17 and BLYB-15). Notably, two strains (CYAH-6 and BLYB-15) were identified as high nitrogen-fixing bacteria in both variety 1 and variety 2 (Figure 5).

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Symbiotic performance score of 46 representative strains based on three symbiotic indexes: chlorophyll index (25%) + number of nodules (25%) + dry plant weight (50%). CK is negative control (uninoculated plants). (Left): Scores with black *Phaseolus vulgaris*. (Right): Score with white *P. vulgaris*. Treatments (CK or strains) were conducted in triplicates. Bars indicate mean ± standard error. Bars with a same letter are not significantly different (ANOVA + LSD test).</figcaption></figure>

## 2.7. Abiotic Stress Tolerance of Representative Strains

The experimental results obtained from growth media tests showed that among the 46 strains tested, none could grow at pH 5. Twenty-one strains could grow normally at a pH of 6 (Supplementary Table S3). All strains thrived under pH 7 and pH 8 conditions. Forty strains grew normally at pH 9, while only twenty strains could tolerate pH 10. Remarkably, at pH 11, only nine strains demonstrated normal growth (CMDB-12, CYAH-5, CDYH-7, CDYB-13, CDHH-14, DXYH-4, DXYB-25, BLYH-17 and BLYB-15), indicating their strong alkali tolerance (Supplementary Table S3). All strains exhibited normal grow at a NaCl concentration of 0.01% (control), but as the NaCl concentration increased, strain growth was gradually inhibited. At 1.0% NaCl, only 17 strains could grow normally. When NaCl reached 2.0%, strains CYAH-25 and BLYH-17 showed normal grow and stronger salt tolerance (Supplementary Table S4). None of the strains grew at 4 °C and 10 °C, but all were able to grow at 28 °C. At 37 °C, 12 strains grew normally. At 45 °C, only three strains (CDYH-7, BLYH-17 and BLYB-15) did, indicating their wide temperature adaptation (Supplementary Table S5). All strains were able to grow in presence of 3% and 5% of polyethylene glycol. This number decreased to 31 strains at 7%, and to 11 at 10%. The
11 strains (CMDH-25, CYAH-6, CYAB-3, CDYB-4, CDYB-13, CDYB-22, CDHH-1, CDHH-14, DXYB-25, BLYH-17 and BLYB-15) showed the highest drought tolerance (Supplementary Table S6). The 46 representative rhizobial strains could grow at the lowest concentration of glyphosate but slowly, revealing some growth inhibition in the glyphosate-containing media. This number reached 38 tolerant strains growing in a glyphosate concentration of 1.2 mL/L and 27 strains in 1.8 mL/L, indicating the highest glyphosate tolerance (Supplementary Table S7).

The determination of abiotic stress tolerance across 46 symbiotic strains of *P. vulgaris* showed a rich diversity in terms of adaptability among the different strains. The comprehensive experimental results, encompassing all the abiotic stress tested, i.e., acidity, alkalinity, NaCl, temperature, PEG and glyphosate tolerance, revealed CWDB-3, DXYH-4, BLYH-17 and BLYB-15 as the most tolerant rhizobial strains (Supplementary Table S8). Combined with symbiotic assay results (Section 2.6), BLYB-15 (belonging to *R. dioscoreae*, IGS type 17 according to Section 2.3.) emerged as the strain with both high potential symbiotic efficiency scores on v1 and v2 local bean varieties, together with notable environmental stress tolerance potential (growth at pH 6, pH 11, 45 °C, with 1% NaCl, 10% PEG or 1.8 mL/L glyphosate).

### 3. Discussion

The study of *P. vulgaris* nodulating rhizobia in Yunnan, combined with the examination of soil physicochemical properties and associated environmental factors, revealed the bio-geographic distribution of rhizobia in the area. Unlike previous studies, we systematically investigated the diversity of *P. vulgaris* symbionts at 15 sampling sites in Chuxiong, Dali, and Baoshan areas of Yunnan Province. A total of 608 rhizobial isolates were thus obtained from root nodules and characterized genetically and symbiotically. Taking all the results from the IGS PCR-RFLP typing, phylogenies of 16S rRNA gene sequences, concatenated *recA-atpD-rpoB* sequences, as well as *nodC* and *nifH* sequences, the isolates obtained in this study were classified as a diverse community consisting of 43 IGS types within 16 species all within the genus *Rhizobium*. Most of the representatives selected by IGS type and site belonged to the sv. phaseoli and a minority to sv. viciae. The dominant species among common bean-nodulating rhizobia were *R. anhuiense* (28.8% in relative abundance) and *R. phaseoli* (28.0%), followed by *R. croatiense* (14.8%) and *R. hidalgonense* (5.3%). These species were found to be widely distributed across most of the fifteen sites (11–13 sites). The twelve other species (*R. sophorae*, *R. acidisoli*, *R. ecuadorense*, *R. vallis*, *R. sophoriradicis*, *R. etli*, *R. chutanense* and *Rhizobium* genosp. I~V) were less represented with a narrower distribution (<5.6% and restricted to 1–7 sites).

The data provide new insights on the diversity structure and geographic distribution of *P. vulgaris*-symbiotic rhizobia in China and areas outside of the plant-host native range. Among the recovered species in Yunnan areas, *R. sophorae* was initially isolated from effective nodules of the shrubby Sophora (Sophora flavescens) in Changzhi City (Shanxi Province, China), and was found to be able to effectively nodulate not only Sophora but also *P. vulgaris* [51]. This rhizobial species was later recovered from root nodules of *Vicia faba* L. grown in Panxi (southwest China) [52], southwest China [53], and Hebei Province (northeast China) [54,55]. These findings underscore the widespread existence of *R. sophorae*, emphasizing its adaptability to diverse soil and environmental conditions in China. *R. acidisoli*, initially isolated from *P. vulgaris* nodules in acidic soils of Mexico and forming symbiosis with *P. vulgaris* [15], was later discovered in Morocco [56]. *R. ecuadorense*, isolated from *P. vulgaris* in northern and central Ecuador, had a valuable effective nitrogen fixation with *P. vulgaris* [21]. *R. hidalgonense* was isolated from *P. vulgaris* nodules in the State of Mexico [20]. *R. anhuiense* was first described in rhizobia of pea and *V. faba* in Anhui and Jiangxi Provinces [57], and has been widely recorded in several other provinces in China (provinces of Shandong [58], Sichuan [59] and Hebei [55]). The present study, with a total of 28% relative abundance detected over 13 Yunnan sites, suggests a wide geographic distribution of *R. anhuiense* in China. *R. croatiense* was isolated from *P. vulgaris* landraces
from soils of northeast Croatia [16]. *R. vallis*, isolated from nodules of three legume plants (*P. vulgaris*, *Mimosa mimosa* and *Indigofera spicata*) grown in the Yunnan province of China, effectively nodulated *P. vulgaris* but did not nodulate *M. mimosa* and *I. spicata*, suggesting that *R. vallis* could be endophytic in some root nodules [31]. *R. sophoriradicis*, originally isolated from Sophora (S. flavescens) [51], has been extensively studied in *P. vulgaris* in Iran [60], South Africa [61] and Peru [62], indicating its widespread distribution among plant species. *R. phaseoli*, which predominates in native areas of common beans, was isolated in Mexico [63], but also from non-native countries such as Ethiopia [26], Brazil [64] and Eswatini [65]. *R. chutanense* was isolated from *P. vulgaris* in Jiangxi Province for the first time, and can effectively trigger nodulation with both *P. vulgaris* and soybean [32]. *R. etli*, rarely found in this study (<1%), was initially discovered in root nodules of the Mexican leguminous plant *Mimosa affinis* and was able to form nodules and fix nitrogen on *M. mimosa* [66]. This species was subsequently found in *P. vulgaris* nodules collected across different agro-ecological zones in Senegal, Gambia (West Africa) [67] and northwestern Argentina. It emerged as a predominant species in common bean nodules from *P. vulgaris* origin areas [24,68] and from different regions of Jordan [69] and Brazil [70]. Moreover, *R. etli* was widely found to be tolerant to high salinity and pH levels in northwestern Morocco [71], Ethiopia [26,72] and Egypt [73,74]. These results show the wide distribution of rhizobia in time and space, as well as its high adaptability to different soil and environmental factors. Adding to the diversity distribution. Notably, *R. tropici* was not observed in this study, similar to findings in the Shaanxi Province of China [34]. However, *R. tropici*, along with *R. etli* and *R. phaseoli*, is another predominant species among *P. vulgaris* symbionts widely recovered in continents and countries such as Columbia [14], Argentina [68], Brazil [64,75] and Iran [60], as well as in north, west or South Africa [61,67,75]. Bacteria of the genus *Rhizobium* have been isolated from a variety of sources and widely utilized for their nitrogen-fixing abilities, both in environment and agriculture. Our findings highlight the broad geographic distribution of these species worldwide, suggesting that *P. vulgaris* plants have selected chromosomal backgrounds or communities suited to local conditions over long-term legume cultivation. Ultimately, the detection of highly conserved *nodC* genes across the 13 rhizobial species identified in this study provides further evidence of *P. vulgaris* plants' stringent selection of symbiosis genes in their microsymbionts. This selective pressure may favor symbiosis genes toward the most adapted indigenous rhizobia, as reported in other cases [76]. To sum up, our results underscore the necessity of screening and selecting high-quality rhizobial strains with strong adaptability to local conditions for inoculant production aimed at enhancing *P. vulgaris* and other legume yields.

The values of OM, AN, pH and AK varied significantly among the field soils sampled in this study, and the rhizobial distribution at the IGS type level was found related to these soil characteristics. In addition, we observed that the same species isolated from different soil types had different abundance responses according to soil physicochemical characteristics, as evidenced by a distinct pattern observed at the IGS genotypes within species (Figure 4). Soil chemical and climatic factors exerted varying effects on the distribution of the rhizobial species and IGS types. For instance, we noted a positive correlation between IGS types 3, 6, 15 and 23 (*R. croatiense*), 17 (*R. genosp*. IV), 19 (*R. genosp*. V), 20 (*R. sophorae*), 30 (*R. genosp*. I) and 43 (*R. hidalgonense*) and a negative correlation of IGS types 5, 35 (*R. chutanense*), 24 (*R. genosp*. III), 8, 42, (*R. hidalgonense*), 10 (*R. genosp*. I), 13 (*R. ecuadorense*) and 36 (*R. phaseoli*) with average precipitation (Figure 4). Similarly, a positive correlation was observed between IGS types 2, 4, 9, 16, 26, 33 and 40 (*R. phaseoli*), 12 and 39 (*R. vallis*), 14 (*R. genosp*. II), 41 (*R. hidalgonense*), 18 (*R. acidisoli*), 21 (*R. croatiense*), 31 (*R. etli*) and 28 (*R. sophoriradicis*) and a negative correlation of IGS types 1, 11, 29 and 34 (*R. anhuiense*), 7, 22 and 38 (*R. phaseoli*), 25, 32 and 37 (*R. hidalgonense*) and 27 (*R. chutanense*), with AK, influencing their separation in two distinct distribution patterns. These relationships were consistent with previous reports indicating that soil traits were selective abiotic factors shaping the biogeography of rhizobial species or IGS types [77].
In conclusion, our study demonstrated the existence of indigenous rhizobia that form an effective symbiosis with *P. vulgaris* cultivated in southwest China, in which *R. anhuiense* was newly recorded as *P. vulgaris*-nodulating rhizobia. *R. phaseoli*, *R. anhuiense*, *R. croatiense* and *R. hidalgonense* were the most abundant and widely distributed in the studied soils (28.0%, 28.8%, 14.8% and 5.6%) and formed unique species assemblage of *P. vulgaris*-nodulating rhizobia along with twelve less frequent species (*R. sophorae*, *R. acidisoli*, *R. ecuadorense*, *R. vallis*, *R. sophoriradicis*, *R. etli*, *R. chutanense* and *Rhizobium* genosp. I~V). Moreover *R. anhuiense* occurred in thirteen tested soil types, indicating that strains of these species could be better competitors and adapted to different soil conditions. Finally, all strains tested belonged to the sv. phaseoli regardless of their species affiliation or sv. viciae associated with *R. vallis*. Through symbiotic experiments, the number of nodules per plant was good with chlorophyll index higher than controls for most strains. Dry weights of host plants were significantly improved under inoculation for about half of the *P. vulgaris*-nodulating strains. Therefore, inoculation with rhizobia associated with *P. vulgaris* can significantly promote the growth of *P. vulgaris* plants and achieve effective symbiosis. In parallel, representative strains to IGS types and sites were screened for abiotic stress tolerance, aiming to potentially assist the host legume to cope with soil and environmental stresses. A total of four strains with strong comprehensive tolerance were identified. Combined with the results of symbiotic assay on both *P. vulgaris* varieties, the strain BLYB-15 (*R. genosp*. IV sv. phaseoli, IGS type 17) was found to exhibit high symbiotic efficiency and strong stress tolerance. The resource collection of common bean-associated rhizobia in Yunnan Province, analyzed for genetic and phenotypic diversity, provides significant guidelines for increasing local and sustainable *P. vulgaris* production, with the potential for site-specific selection of efficient rhizobial genetic types.

# 4. Materials and Methods

## 4.1. Field Soil Sampling and Soil and Environmental Characteristics

Soils were sampled from fields cultivated with local common beans (*Phaseolus vulgaris* L.) located in the cities of Chuxiong (CX), Dali (DL) and Baoshan (BS), all in the province of Yunnan, southwest China. A total of 15 sites were studied: 8 at CX (CX-MD, YA, DY LF, DH, WD, SB, NH), 4 at DL (DL-XY, WS, MD, EY) and 3 at BS (BS-LY, CN, SD) (Supplementary Table S1). At each site, soil was collected from the main production area of *P. vulgaris* to a depth of 10–20 cm near the *P. vulgaris* roots during its flowering stage in June 2022. For each site, three randomly taken soil sub-samples of equal volume were crushed to a uniform state, and transported to the laboratory in an ice-filled cooler [77]. Part of each soil sample was chemically analyzed for pH, electrical conductivity (EC), organic matter (OM), available phosphorus (AP), available potassium (AK) and alkaline hydrolyzable nitrogen (AN), as described previously [78]. Using the location information of the sampled sites, annual climate data were collected for each site using DIVA-GIS software Version 7.5 (University of California, Los Angeles, CA, USA). This includes the altitude of the sampling site (Alt), the amount of rainfall throughout the year (AvePrecp), the maximum (AveTmax) and minimum temperatures (AveTmin) and the analysis of climate data [79].

## 4.2. Rhizobial Isolation and Conservation

Surface-sterilized seeds (2.5% *w*/*v* NaClO solution for 5 min) of the common bean (variety 1 of black *P. vulgaris*: garden bean and variety 2 of white *P. vulgaris*: edible pod bean from the local areas of Yunnan Province) were germinated and the seedings were sown in surface-sterilized plastic pots (15 cm high × 10 cm diameter) filled with each sampled soil mixed with sterilized vermiculite (1/5 *v*/*v*). The three randomly taken soil sub-samples were thoroughly mixed to constitute a representative soil sample per field site. For each representative soil sample, 10 repetitions of tracking experiments were performed. All plants were grown under greenhouse conditions of 25/20 °C (day/night) with a 16 h photoperiod. Sterilized water was added to the pots throughout the experiment as required. After 45 days, all plants in all soils were uprooted and bacterial strains were isolated from
nodules according to the standard protocol [53]. Five plants from each sampled soil were selected for excision of root nodules and isolation of rhizobia. Root nodules were surface-sterilized, then each individual sterilized root nodule was crushed in sterile water and the bacterial suspension was streaked onto a Yeast extract Mannitol Agar (YMA) plate. After incubation at 28 °C for 2 to 3 days, single colonies representing the dominant bacteria on each plate were picked up and purified by cross-streaking on new YMA plates until pure cultures were visually obtained. All purified isolates were conserved in Tryptone Yeast (TY) broth (tryptone 5 g; yeast extract 3 g; CaCl₂ 0.6 g; distilled water 1 L, pH 7.0) supplied with glycerol (20%, *v/v*) and stored at -80 °C for long-term storage. Additionally, they were maintained on YMA slants at 4 °C for temporary storage.

### 4.3. Genomic Characterization of Rhizobial Isolates

The genomic DNA of each isolate was purified according to Terefwork et al. [80,81]. DNA was used as template for PCR amplifications of the 16S-23S rRNA intergenic spacer (IGS) region with primers IGS1490 (forward, TGCGGCTGGATCACCTCCTT) and IGS132' (reverse, CCGGGTTTCCCCATTCGG) [82]. PCR amplification was carried out in a standard 50 µL reaction mixture, including 1 µL of DNA template and 5 U of Taq DNA polymerase (Sangon Biotech (Shanghai) Co., Ltd., Shanghai, China). Aliquots of amplified PCR products (900 base pairs) were visualized after electrophoresis in a 1.0% (*w/v*) agarose gel labeled with GoldView type I. Then, PCR products were digested separately with the endonucleases HaeIII, MspI and Hhal [83] at 37 °C for 10 h. The 16S-23S rRNA gene IGS type of each strain was designated after separation and visualization of restriction fragments by electrophoresis in 2.5% (*w/v*) agarose gel and UV-illumination.

### 4.4. Molecular and Phylogenetic Identification of the Isolates, Alpha-Diversity Estimation

Isolates sharing the same RFLP pattern of 16S-23S rRNA gene IGS in this study were designed as an IGS type. One representative strain for each IGS type in each sample site (total of 46 strains) was selected for amplification of the 16S rRNA gene using the forward primer P1 (CGGGATCCAGAGTTTGATCCTGGTCAGAACGCT) and reverse primer P6 (CGGGATCCTACGGCTACCTTGTTAC GACTTCACCCC3) [82]. The PCR products were verified as mentioned above, and were sent for commercial sequencing based on the Sanger method (Sangon Biotech (Shanghai, China) Co., Ltd.). The acquired sequences were compared against the NCBI database using the online BLASTN tool, and sequences for type strains of defined *Rhizobium* species sharing similarities greater than 97.0% with the new isolates were extracted. The phylogenetic analysis was conducted in the MEGA 7.0 software [84]. Sequences were aligned using Clustal W and the best model of sequence evolution was selected. Then, the phylogenetic tree was inferred using the maximum likelihood (ML) and the non-parametric bootstrap (500 pseudo-replications) methods.

DNA fragments of *recA* (coding for DNA recombination protein), *atpD* (encoding for ATP synthase beta chain) and *rpoB* (encoding the RNA polymerase β subunit) were amplified separately by PCR using the primer pairs recA41F/recA640R (TTCGGCAAGGGMTC-GRTSATG/ACATSACRCCGATCTTCATGC), atpD255F/atpD782R (GCTSGGCCGCATCM-TSAACGTC/GCCGACACTTCMGAACCNGCCTG) and rpoB83F/rpoB1061R (CCTSATC-GAGGTTCACAGAAGGC/AGCGTGTTGCGGATATAGGCG), respectively [85,86]. The symbiotic genes *nodC* and *nifH* were amplified using the primer pairs nodC540F/nodC1160R (TGATYGAYATG GARTAYTGGCT/CGYGACARCCARTCGCTRTTG) and nifHF/nifHR (TACGGNAARGGSGGNATCGGCAA/AGCATGTCYTCSAGY TCNTCCA) [87,88]. The PCR conditions were adopted from previous reports [88]. PCR products verification, sequencing and tree construction of each sequenced gene were performed as mentioned above. Furthermore, sequences of *atpD* (297 base pairs), *recA* (289 base pairs) and *rpoB* (599 base pairs) genes were concatenated and aligned using Clustal W [89]. Distance calculation, construction of ML trees based on each simple gene and the concatenated housekeeping genes, bootstrap analysis, were performed in MEGA 7.0 as described above.
The sequences have been deposited in the NCBI database (accession numbers are indicated on the trees (Figures 1 and 2 and Supplementary Figures S1–S5)).

Alpha-diversity was calculated to estimate the effective numbers of species (ENS) using three diversity indices (numbers of clusters, exponential Shannon index and the inverse of the Simpson index [90]).

### 4.5. Correlation Analysis of Soil Properties and Environmental Factors with Rhizobial Communities

The principal component analysis (PCA) was performed using CANOCO version 5.0 [91] to investigate the relationships between soil properties (AN, AP, AK, OM, EC and pH) and environmental factors (Alt, AveTmin, AveTmax and AvePrecp) and the rhizobial community composition based on IGS genotypes and rhizobial species. The distance matrix generated from the response variable (i.e., rhizobial composition data) were based on percentage dissimilarity (i.e., Bray-Curtis dissimilarity) obtained from the 43 IGS genotypes (used to identify rhizobia). This matrix was then correlated to environmental and soil factors.

### 4.6. Symbiotic Efficiency Measurements

Symbiotic efficiency of the representative strains was evaluated on both variety 1 and variety 2 of the local common bean. Briefly, surface-sterilized seedlings were aseptically transferred in pots (1 plant/pot) containing sterile vermiculite as substrate and inoculated with 1 mL of rhizobial suspension (OD<sub>600</sub> = 1.0). Plants were grown under greenhouse conditions and watered with sterile N-free nutrient solution as required [92]. Symbiotic performance was evaluated 45 days after inoculation for both bean varieties. Common bean growth was estimating by weighting their dry root and shoot biomass, counting the number of nodules and measuring leaf chlorophyll contents (SPAD chlorophyll meter). Uninoculated plants were included as a negative control, and all treatments were performed in three replicates. A symbiotic performance score (%) was calculated for each treatment, considering 25% of the nodulation index (number of nodules per plant), 25% of the leaf chlorophyll index (SPAD value per plant) and 50% of the plant biomass index (dry roots and leaves per plant). Each performance index was standardized by dividing it by the best average value obtained in the experiment (n = 3). Data were analyzed by one-way ANOVA followed by an LSD post hoc test (p = 0.001).

### 4.7. Measurements of Abiotic Stress on Rhizobial Strain Growth

Acid and alkali tolerance test: The 46 representative strains of common bean-nodulating rhizobia were inoculated on YMA solid medium at different pH (5, 6, 7, 8, 9, 10 and 11). Each strain was repeated on 3 plates and incubated at 28 °C. Growth results were recorded after 2 and 3 days in comparison with the control at pH 7.

Salt tolerance test: Representative rhizobial strains were inoculated on YMA solid medium containing different concentrations of NaCl (0.01%, 1%, 2%, 3% and 4% (w/v)), with 0.01% NaCl plates as control. Each strain was tested in triplicate, and growth was assessed after 2 and 3 days of incubation at 28 °C.

Temperature growth tolerance test: Each representative rhizobium was inoculated in liquid TY medium and cultured at different temperatures (4 °C, 10 °C, 28 °C, 37 °C and 45 °C). Each treatment was replicated 3 times. Optical density (O.D.) was measured at 600 nm after 2 and 3 days of incubation.

Drought tolerance test: Polyethylene glycol 6000 (PEG 6000) was used to artificially simulate drought conditions. Different amounts of PEG were added to YM broth medium, resulting in final concentration of 0%, 3%, 5%, 7%, 10% and 15% (w/v). Each rhizobium was inoculated into these media (three replicates), and cultures were incubated in a rotary incubator at 28 °C for 2 and 3 days.

Glyphosate tolerance test: Glyphosate was sourced from a commercial product containing 41% (w/v) glyphosate isopropylamine saline solution in H<sub>2</sub>O, at an effective concentration of 2.16 M. Each rhizobial strain was inoculated on YMA medium supplemented
with different concentrations of glyphosate (by adding 0, 0.6 mL/L, 1.2 mL/L and 1.8 mL/L of the 41% glyphosate isopropylamine saline solution). Plates were then incubated at 28 °C for 2 and 3 days.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25158511/s1.

**Author Contributions:** Conceptualization, J.Z.; Methodology, J.Z.; Validation, J.Z. and B.B.; Formal analysis, J.W.; Investigation, J.W. and Y.F.; Resources, X.Z.; Writing—original draft, J.W.; Writing—review & editing, J.Z., B.B. and X.Z.; Supervision, J.Z.; Project administration, J.Z.; Funding acquisition, J.Z. and X.Z. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was financed by the National Key R&D Program of China (2021YFD1600602-02/YFD1600602) to X.X. Zong, and the Science and Technology Innovation Talents in Universities of Henan Province (22HASTIT035) and Central Plains Youth Top Talent Project (Yuzutong[2023]No.11) to J.J. Zhang. The APC was funded by 22HASTIT035 and Yuzutong[2023]No.11.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Guo, Z. High yield cultivation and management technology of kidney beans in North China. *Agric. Technol. Equip.* **2017**, *11*, 43–45.

2. Luo, A.; Li, M.; Liu, F.; Fu, Q.; Dong, J. Research progress in genetic breeding of common bean. *Vegetables* **2023**, 36–42.

3. Grange, L.; Hungria, M.; Graham, P.; Martinez-Romero, E. New insights into the origins and evolution of rhizobia that nodulate common bean (*Phaseolus vulgaris*) in Brazil. *Soil Biol. Biochem.* **2007**, *39*, 867–876. [CrossRef]

4. Zhang, X.; Blair, M.W.; Wang, S. Genetic diversity of Chinese common bean (*Phaseolus vulgaris* L.) landraces assessed with simple sequence repeat markers. *Theor. Appl. Genet.* **2008**, *117*, 629–640. [CrossRef]

5. Ma, X.; Li, L. 16S rDNA sequence analysis of Rhizobium strains isolated from Soybeans and kidney beans in Rizhao. *Hubei Agric. Sci.* **2017**, *56*, 2938–2941. [CrossRef]

6. Wu, J. Research progress on common bean genomics study. *J. Sichuan Agric. Univ.* **2021**, *39*, 4–10. [CrossRef]

7. Yao, Y.; Jie, W.; Du, Y.; Zhao, D.; Yan, X. Taxonomy, Identification and application of Rhizobium. *Chin. Agric. Sci. Bull.* **2020**, *36*, 100–105.

8. Khalid, R.; Zhang, Y.; Ali, S.; Sui, X.; Zhang, X.; Amara, U.; Chen, W.; Hayat, R. *Rhizobium pakistanensis* sp. nov., isolated from groundnut (*Arachis hypogaea*) nodules grown in rainfed Pothwar, Pakistan. *Antonie Van Leeuwenhoek* **2014**, *107*, 281–290. [CrossRef] [PubMed]

9. Liu, X.; Li, Y.; Li, Z. The Effects of Intercropping and Dual Inoculation on Maize and Common Bean (*Phaseolus vulgaris*) Growth. *Chin. Agric. Sci. Bull.* **2014**, *30*, 182–188.

10. Sena, P.T.S.; do Nascimento, T.R.; Lino, J.d.O.S.; Oliveira, G.S.; Ferreira Neto, R.A.; de Freitas, A.D.S.; Fernandes-Júnior, P.I.; Martins, L.M.V. Molecular, physiological, and symbiotic characterization of cowpea rhizobia from soils under different agricultural systems in the semiarid region of Brazil. *J. Soil Sci. Plant Nutr.* **2020**, *20*, 1178–1192. [CrossRef]

11. Kumar, N.; Srivastava, P.; Vishwakarma, K.; Kumar, R.; Kuppala, H.; Maheshwari, S.K.; Vats, S. The *Rhizobium*-plant symbiosis: State of the art. In *Plant Microbe Symbiosis*; Varma, A., Tripathi, S., Prasad, R., Eds.; Springer: Cham, Switzerland, 2020.

12. Yan, Y.; Tian, C.; Yang, J.; Wang, Y.; Lin, M. Establishment of artificial efficiency biological nitrogen fixation system and its agricultural application. *Chin. Bull. Life Sci.* **2021**, *33*, 1532–1543. [CrossRef]

13. Xiang, Z.; Song, Y. Breeding of *Rh. phaseoli* and its effective application technology. *J. Microbiol.* **1988**, 40–44.

14. Martínez-Romero, E.; Segovia, L.; Mercante, F.M.; Franco, A.A.; Graham, P.; Pardo, M.A. *Rhizobium tropici*, a novel species nodulating *Phaseolus vulgaris* L. beans and *Leucaena* sp. trees. *Int. J. Syst. Bacteriol.* **1991**, *41*, 417–426. [CrossRef] [PubMed]

15. Roman-Ponce, B.; Zhang, Y.J.; Soledad Vasquez-Murrieta, M.; Sui, X.H.; Chen, W.F.; Alberto Padilla, J.C.; Guo, X.W.; Gao, J.L.; Yan, J.; Wei, G.H.; et al. *Rhizobium acidisoli* sp nov., isolated from root nodules of *Phaseolus vulgaris* in acid soils. *Int. J. Syst. Evol. Microbiol.* **2016**, *66*, 398–406. [CrossRef] [PubMed]

16. Rajnovic, I.; Ramírez-Bahena, M.-H.; Kajic, S.; Igual, J.M.; Peix, Á.; Velázquez, E.; Sikora, S. *Rhizobium croatiense* sp. nov. and *Rhizobium redzepovicii* sp. nov., two new species isolated from nodules of *Phaseolus vulgaris* in Croatia. *Syst. Appl. Microbiol.* **2022**, *45*, 126317. [CrossRef] [PubMed]
17. Dall'Agnol, R.F.; Ribeiro, R.A.; Ormeno-Orrillo, E.; Rogel, M.A.; Hungria, M. *Rhizobium freirei* sp. nov., a symbiont of *Phaseolus vulgaris* that is very effective at fixing nitrogen. *Int. J. Syst. Evol. Microbiol.* **2013**, *63* Pt 11, 4167–4173. [CrossRef] [PubMed]

18. Aurag, J.; Sasson, A. Tolerance of Rhizobium leguminosarum bv. phaseoli to acidity and drought. *World J. Microbiol. Biotechnol.* **1992**, *8*, 532–533. [CrossRef]

19. Segovia, L.; Young, J.P.W.; Martinez-Romero, E. Reclassification of American *Rhizobium leguminosarum* Biovar Phaseoli Type I Strains as *Rhizobium etli* sp. nov. *Int. J. Syst. Bacteriol.* **1993**, *43*, 374–377. [CrossRef] [PubMed]

20. Yan, J.; Yan, H.; Liu, L.; Chen, W.; Zhang, X.; Verástegui-Valdés, M.M.; Wang, E.; Han, X. *Rhizobium hidagonense* sp. nov., a nodule endophytic bacterium of *Phaseolus vulgaris* in acid soil. *Arch. Microbiol.* **2016**, *199*, 97–104. [CrossRef]

21. Ribeiro, R.A.; Martins, T.B.; Ormeno-Orrillo, E.; Marcon Delamuta, J.R.; Rogel, M.A.; Martinez-Romero, E.; Hungria, M. *Rhizobium ecuadorense* sp nov., an indigenous N₂-fixing symbiont of the Ecuadorian common bean (*Phaseolus vulgaris* L.) genetic pool. *Int. J. Syst. Evol. Microbiol.* **2015**, *65* Pt 9, 3162–3169. [CrossRef]

22. Amarger, N.; Macheret, V.; Laguerre, G. *Rhizobium gallicum* sp. nov. and *Rhizobium giardinii* sp. nov., from *Phaseolus vulgaris* Nodules. *Int. J. Syst. Bacteriol.* **1997**, *47*, 996–1006. [CrossRef] [PubMed]

23. Ribeiro, R.A.; Rogel, M.A.; Lopez-Lopez, A.; Ormeno-Orrillo, E.; Barcellos, F.G.; Martinez, J.; Thompson, F.L.; Martinez-Romero, E.; Hungria, M. Reclassification of *Rhizobium tropici* type A strains as *Rhizobium leucaenae* sp nov. *Int. J. Syst. Evol. Microbiol.* **2012**, *62*, 1179–1184. [CrossRef] [PubMed]

24. Martínez-Romero, E. Diversity of *Rhizobium-Phaseolus vulgaris* symbiosis: Overview and perspectives. *Plant Soil* **2003**, *252*, 11–23. [CrossRef]

25. Mnasri, B.; Liu, T.Y.; Saidi, S.; Chen, W.F.; Chen, W.X.; Zhang, X.X.; Mhamdi, R. *Rhizobium azibense* sp nov., a nitrogen fixing bacterium isolated from root-nodules of *Phaseolus vulgaris*. *Int. J. Syst. Evol. Microbiol.* **2014**, *64* Pt 5, 1501–1506. [CrossRef] [PubMed]

26. Aserse, A.A.; Rasanen, L.A.; Assefa, F.; Hailemariam, A.; Lindstrom, K. Phylogeny and genetic diversity of native rhizobia nodulating common bean (*Phaseolus vulgaris* L.) in Ethiopia. *Syst. Appl. Microbiol.* **2012**, *35*, 120–131. [CrossRef] [PubMed]

27. Lopez-Lopez, A.; Rogel-Hernandez, M.A.; Barois, I.; Ortiz Ceballos, A.I.; Martinez, J.; Ormeno-Orrillo, E.; Martinez-Romero, E. *Rhizobium grahamii* sp nov., from nodules of *Dalea leporina*, *Leucaena leucocephala* and *Clitoria ternatea*, and *Rhizobium mesoamericanum* sp nov., from nodules of *Phaseolus vulgaris*, *siratro*, *cowpea* and *Mimosa pudica*. *Int. J. Syst. Evol. Microbiol.* **2012**, *62*, 2264–2271. [CrossRef] [PubMed]

28. Luis Zurdo-Pineiro, J.; Garcia-Fraile, P.; Rivas, R.; Peix, A.; Leon-Barrios, M.; Willems, A.; Francisco Mateos, P.; Martinez-Molina, E.; Velazquez, E.; van Berkum, P. Rhizobia from Lanzarote, the Canary Islands, that nodulate *Phaseolus vulgaris* have characteristics in common with *Sinorhizobium meliloti* Isolates from Mainland Spain. *Appl. Environ. Microbiol.* **2009**, *75*, 2354–2359. [CrossRef] [PubMed]

29. Mnasri, B.; Saïdi, S.; Chihaoui, S.-A.; Mhamdi, R. *Sinorhizobium americanum* symbiovar mediterranense is a predominant symbiont that nodulates and fixes nitrogen with common bean (*Phaseolus vulgaris* L.) in a Northern Tunisian field. *Syst. Appl. Microbiol.* **2012**, *35*, 263–269. [CrossRef] [PubMed]

30. Sadowsky, M.J.; Cregan, P.B.; Keyser, H.H. Nodulation and nitrogen fixation efficacy of *Rhizobium fredii* with *Phaseolus vulgaris* genotypes. *Appl. Environ. Microbiol.* **1988**, *54*, 1907–1910. [CrossRef]

31. Wang, F.; Wang, E.T.; Wu, L.J.; Sui, X.H.; Li, Y.; Chen, W.X. *Rhizobium vallis* sp. nov., isolated from nodules of three leguminous species. *Int. J. Syst. Evol. Microbiol.* **2011**, *61*, 2582–2588. [CrossRef]

32. Huo, Y.; Tong, W.; Wang, J.; Wang, F.; Bai, W.; Wang, E.; Shi, P.; Chen, W.; Wei, G. *Rhizobium chutanense* sp. nov., isolated from root nodules of *Phaseolus vulgaris* in China. *Int. J. Syst. Evol. Microbiol.* **2019**, *69*, 2049–2056. [CrossRef] [PubMed]

33. Cao, Y. Diversity and Population Structures of Rhizobia Nodulated with Phaseolus vulgaris in Two Ecoregions of China. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2015.

34. Wang, L.; Cao, Y.; Wang, E.T.; Qiao, Y.J.; Jiao, S.; Liu, Z.S.; Zhao, L.; Wei, G.H. Biodiversity and biogeography of rhizobia associated with common bean (*Phaseolus vulgaris* L.) in Shaanxi Province. *Syst. Appl. Microbiol.* **2016**, *39*, 211–219. [CrossRef]

35. Tong, W. Comparative Genomics and Phylogeny of Rhizobia Isolated from *Phaseolus vulgaris*. Ph.D. Thesis, Northwest A&F University, Yangling, China, 2019.

36. Anyango, B.; Wilson, K.J.; Beynon, J.L.; Giller, K.E. Diversity of rhizobia nodulating *Phaseolus vulgaris* L. in two Kenyan Soils with contrasting pHs. *Appl. Environ. Microbiol.* **1995**, *61*, 4016–4021.

37. Efstathiadou, E.; Savvas, D.; Tampakaki, A.P. Genetic diversity and phylogeny of indigenous rhizobia nodulating faba bean (*Vicia faba* L.) in Greece. *Syst. Appl. Microbiol.* **2020**, *43*, 126149. [CrossRef] [PubMed]

38. Xia, C. Identification of Common *Phaseolus vulgaris* Germplasm Resources Based on Phenotype and SSR Molecular Markers. Master's Thesis, Heilongjiang Bayi Agricultural Reclamation University, Daqing, China, 2020.

39. Satyanandam, T.; Babu, K.; Yellamanda, B.; Kumar, K.V.; Rosaiah, G.; Vijayalakshmi, M. Diversity of indigenous symbiotic nitrogen fixing bacteria from blackgram *Vigna mungo* (L.) Hepper cultivated in rice fallows. *Legume Res.* **2022**, *45*, 994–999. [CrossRef]

40. Wang, X.; Yang, J.; Xu, J.; Az, G. Isolation and back-grafting test of rhizobia from kidney bean. *Xinjiang Agric. Sci.* **2001**, *76–77*. [CrossRef]

41. Savy, D.; Cozzolino, V. Novel fertilising products from lignin and its derivatives to enhance plant development and increase the sustainability of crop production. *J. Clean. Prod.* **2022**, *366*, 132832. [CrossRef]

42. Xiang, Z.; Song, Y. Breeding and application of rhizobia vulgaris. *Liaoning Agric. Sci.* **1985**, 35–39.
43. Wei, G.; Ma, Z. Application of rhizobia-legume symbiosis for remediation of heavy-metal contaminated soils. *Acta Microbiol. Sin.* **2010**, *50*, 1421–1430. [CrossRef]

44. Yu, B.; Liu, Y. Recent progress in study on salt tolerance in soybean. *Soybean Sci.* **2000**, *19*, 154–159.

45. Han, M. Drought and saline tolerance of *Vicia faba* Rhizobium. *J. Qinghai Univ.* **2019**, *37*, 35–41. [CrossRef]

46. Chi, Y.C.; Wang, J.H.; Fan, T.Q.; Yu, S.L.; Jiang, C.J.; Zhou, G.; Li, H.Z. Preliminary study on NaCl-tolerance and drought-tolerance of bradyrhizobial strains (*Arachis*) isolated from Shandong. *J. Peanut Sci.* **2008**, *21*–25.

47. Wu, P.; He, Q.y.; Li, Z.p.; Shi, J.; Zhu, C.w.; Sheng, W. Screening of *Rhizobium* from soybean in North Anhui for strong resistance to stresses. *J. Anhui Sci. Technol. Univ.* **2012**, *26*, 15–22.

48. Cheng, B.; Wang, J.; Shi, H.b.; Zhang, H.m.; Zhang, R.q.; Zhen, C.; Zhai, B. Effects of glyphosate and Rhizobium on the quality traits and nitrogen fixation of *Medicago sativa*. *Chin. J. Grassl.* **2021**, *43*, 47–53. [CrossRef]

49. Chen, Z.; Jin, Y.; Yao, X.; Wei, X.; Li, X.; Li, C.; White, J.F.; Nan, Z. Gene analysis reveals that leaf litter from Epichloe endophyte-infected perennial ryegrass alters diversity and abundance of soil microbes involved in nitrification and denitrification. *Soil Biol. Biochem.* **2021**, *154*, 108123. [CrossRef]

50. Xu, J. Screening and Molecular Marker of High Efficient Symbiotic Rhizobia on Kidney Bean. Master's Thesis, China Agricultural University, Beijing, China, 2004.

51. Jiao, Y.S.; Yan, H.; Ji, Z.J.; Liu, Y.H.; Sui, X.H.; Wang, E.T.; Guo, B.L.; Chen, W.X.; Chen, W.F. *Rhizobium sophorae* sp nov and *Rhizobium sophoriradicis* sp nov., nitrogen-fixing rhizobial symbionts of the medicinal legume *Sophora flavescens*. *Int. J. Syst. Evol. Microbiol.* **2015**, *65*, 497–503. [CrossRef]

52. Chen, Y.X.; Zou, L.; Penttinen, P.; Chen, Q.; Li, Q.Q.; Wang, C.Q.; Xu, K.W. Faba Bean (*Vicia faba* L.) nodulating rhizobia in Panxi, China, are diverse at species, plant growth promoting ability, and symbiosis related gene levels. *Front. Microbiol.* **2018**, *9*, 1338. [CrossRef]

53. Zhang, J.; Shang, Y.; Peng, S.; Chen, W.; Wang, E.; de Lajudie, P.; Li, B.; Guo, C.; Liu, C. *Rhizobium sophorae*, *Rhizobium laguerreae*, and two novel *Rhizobium genospecies* associated with *Vicia sativa* L. in Northwest China. *Plant Soil* **2019**, *442*, 113–126. [CrossRef]

54. Zhang, J.; Li, S.; Wang, N.; Yang, T.; Brunel, B.; Andrews, M.; Zong, X.; Wang, E. *Rhizobium sophorae* is the dominant rhizobial symbiont of *Vicia faba* L. In North China. *Syst. Appl. Microbiol.* **2022**, *45*, 126291. [CrossRef]

55. Zhang, J.; Wang, N.; Li, S.; Brunel, B.; Wang, J.; Feng, Y.; Yang, T.; Zong, X. Genotypic composition and performance of pea-nodulating rhizobia from soils outside the native plant-host range. *Front. Microbiol.* **2023**, *14*, 1201140. [CrossRef]

56. El Idrissi, M.M.; Lamin, H.; Bouhnik, O.; Lamrabet, M.; Alami, S.; Jabrone, Y.; Bennis, M.; Bedmar, E.J.; Abdelmoumen, H. Characterization of *Pisum sativum* and *Vicia faba* microsymbionts in Morocco and definition of symbiovar viciae in *Rhizobium* acidisoli. *Syst. Appl. Microbiol.* **2020**, *43*, 126084. [CrossRef] [PubMed]

57. Zhang, Y.J.; Zheng, W.T.; Everall, I.; Young, J.P.W.; Zhang, X.X.; Tian, C.F.; Sui, X.H.; Wang, E.T.; Chen, W.X. *Rhizobium anhuiense* sp nov., isolated from effective nodules of *Vicia faba* and *Pisum sativum*. *Int. J. Syst. Evol. Microbiol.* **2015**, *65*, 2960–2967. [CrossRef] [PubMed]

58. Li, Y.; Wang, E.T.; Liu, Y.; Li, X.; Yu, B.; Ren, C.; Liu, W.; Li, Y.; Xie, Z. *Rhizobium anhuiense* as the predominant microsymbionts of Lathyrus maritimus along the Shandong Peninsula seashore line. *Syst. Appl. Microbiol.* **2016**, *39*, 384–390. [CrossRef]

59. Zou, L.; Chen, Y.X.; Penttinen, P.; Lan, Q.; Wang, K.; Liu, M.; Peng, D.; Zhang, X.; Chen, Q.; Zhao, K.; et al. Genetic diversity and symbiotic efficiency of nodulating rhizobia isolated from root nodules of faba bean in one field. *PLoS ONE* **2016**, *11*, e0167804. [CrossRef]

60. Rouhazi, K.; Khodakaramian, G.; Velazquez, E. Phylogenetic diversity of rhizobial species and symbiovars nodulating *Phaseolus vulgaris* in Iran. *FEMS Microbiol. Lett.* **2016**, *363*, fnw024. [CrossRef]

61. Zinga, M.K.; Jaiswal, S.K.; Dakora, F.D. Presence of diverse rhizobial communities responsible for nodulation of common bean (*Phaseolus vulgaris*) in South African and Mozambican soils. *FEMS Microbiol. Ecol.* **2017**, *93*, fiw236. [CrossRef] [PubMed]

62. Ormeno-Orrillo, E.; Aguilar-Cuba, Y.; Zuniga-Davila, D. Draft Genome Sequence of *Rhizobium sophoriradicis* H4, a nitrogen-fixing bacterium associated with the leguminous plant *Phaseolus vulgaris* on the coast of Peru. *Genome Announc.* **2018**, *6*, e00241-18. [CrossRef]

63. Ribeiro, R.A.; Ormeno-Orrillo, E.; Dall'Agnol, R.F.; Graham, P.H.; Martinez-Romero, E.; Hungria, M. Novel *Rhizobium* lineages isolated from root nodules of the common bean (*Phaseolus vulgaris* L.) in Andean and Mesoamerican areas. *Res. Microbiol.* **2013**, *164*, 740–748. [CrossRef]

64. Costa, M.R.; Chibeba, A.M.; Mercante, F.M.; Hungria, M. Polyphasic characterization of rhizobia microsymbionts of common bean *Phaseolus vulgaris* (L.) isolated in Mato Grosso do Sul, a hotspot of Brazilian biodiversity. *Symbiosis* **2018**, *76*, 163–176. [CrossRef]

65. Gunununu, R.P.; Mohammed, M.; Jaiswal, S.K.; Dakora, F.D. Phylogeny and symbiotic effectiveness of indigenous rhizobial microsymbionts of common bean (*Phaseolus vulgaris* L.) in Malkerns, Eswatini. *Sci. Rep.* **2023**, *13*, 17029. [CrossRef]

66. Wang, E.T.; Rogel, M.A.; Garcia-de los Santos, A.; Martinez-Romero, J.; Cevallos, M.A.; Martinez-Romero, E. *Rhizobium etli* bv. mimosae, a novel biovar isolated from Mimosa affinis. *Int. J. Syst. Bacteriol.* **1999**, *49* Pt 4, 1479–1491. [CrossRef] [PubMed]

67. Diouf, A.; de Lajudie, P.; Neyra, M.; Kersters, K.; Gillis, M.; Martinez-Romero, E.; Gueye, M. Polyphasic characterization of rhizobia that nodulate *Phaseolus vulgaris* in West Africa (Senegal and Gambia). *Int. J. Syst. Evol. Microbiol.* **2000**, *50* Pt 1, 159–170. [CrossRef]

68. Aguilar, O.M.; Lopez, M.V.; Riccillo, P.M. The diversity of rhizobia nodulating beans in Northwest Argentina as a source of more efficient inoculant strains. *J. Biotechnol.* **2001**, *91*, 181–188. [CrossRef] [PubMed]
69. Tamimi, S.A.; Young, J.P.W. *Rhizobium etli* is the dominant common bean nodulating rhizobia in cultivated soils from different locations in Jordan. *Appl. Soil Ecol.* **2004**, *26*, 193–200. [CrossRef]

70. Moura, F.T.; Helene, L.C.F.; Ribeiro, R.A.; Nogueira, M.A.; Hungria, M. The outstanding diversity of rhizobia microsymbionts of common bean (*Phaseolus vulgaris* L.) in Mato Grosso do Sul, central-western Brazil, revealing new *Rhizobium* species. *Arch. Microbiol.* **2023**, *205*, 325. [CrossRef] [PubMed]

71. Mouhsine, B.; Prell, J.; Filali-Maltouf, A.; Priefer, U.B.; Aurag, J. Diversity, phylogeny and distribution of bean rhizobia in salt-affected soils of North-West Morocco. *Symbiosis* **2007**, *43*, 83–96.

72. Beyene, D.; Kassa, S.; Ampy, F.; Asseffa, A.; Gebremedhin, T.; van Berkum, P. Ethiopian soils harbor natural populations of rhizobia that form symbioses with common bean (*Phaseolus vulgaris* L.). *Arch. Microbiol.* **2004**, *181*, 129–136. [CrossRef] [PubMed]

73. Shamseldin, A.; Werner, D. Presence of *Rhizobium etli* bv. phaseoli and Rhizobium gallicum bv. gallicum in Egyptian soils. *World J. Microbiol. Biotechnol.* **2007**, *23*, 285–289. [CrossRef]

74. Youseif, S.H.; Abd El-Megeed, F.H.; Ageez, A.; Cocking, E.C.; Saleh, S.A. Phylogenetic multilocus sequence analysis of native rhizobia nodulating faba bean (*Vicia faba* L.) in Egypt. *Syst. Appl. Microbiol.* **2014**, *37*, 560–569. [CrossRef]

75. Grange, L.; Hungria, M. Genetic diversity of indigenous common bean (*Phaseolus vulgaris*) rhizobia in two Brazilian ecosystems. *Soil Biol. Biochem.* **2004**, *36*, 1389–1398. [CrossRef]

76. He, Y.; Guo, L.; Zhang, H.; Huang, G. Symbiotic effectiveness of pea-rhizobia associations and the implications for farming systems in the western Loess Plateau, China. *Afr. J. Biotechnol.* **2011**, *10*, 3540–3548.

77. Zhang, J.; Shang, Y.; Wang, E.; Chen, W.; de Lajudie, P.; Li, B.; Guo, C.; Yang, X.; Zheng, J.; Liu, C. *Mesorhizobium jarvisii* sv. astragali as predominant microsymbiont for *Astragalus sinicus* L. in acidic soils, Xinyang, China. *Plant Soil* **2018**, *433*, 201–212. [CrossRef]

78. Han, L.L.; Wang, E.T.; Han, T.X.; Liu, J.; Sui, X.H.; Chen, W.F.; Chen, W.X. Unique community structure and biogeography of soybean rhizobia in the saline-alkaline soils of Xinjiang, China. *Plant Soil* **2009**, *324*, 291–305. [CrossRef]

79. Appunu, C.; Sasirekha, N.; Prabavathy, V.R.; Nair, S. A significant proportion of indigenous rhizobia from India associated with soybean (*Glycine max* L.) distinctly belong to Bradyrhizobium and Ensifer genera. *Biol. Fertil. Soils* **2009**, *46*, 57–63. [CrossRef]

80. Shang, Y. The Diversity Analysis of Legume Green Manure Symbiotic Rhizobia and Mining of Dominant Species Resources in China. Master's Thesis, Zhengzhou University of Light Industry, Zhengzhou, China, 2020.

81. Terefework, Z.; Kaijalainen, S.; Lindström, K. AFLP fingerprinting as a tool to study the genetic diversity of Rhizobium galegae isolated from Galega orientalis and Galega officinalis. *J. Biotechnol.* **2001**, *91*, 169–180. [CrossRef] [PubMed]

82. Laguerre, G.; Mavingui, P.; Allard, M.R.; Charnay, M.P.; Louvrier, P.; Mazurier, S.I.; Rigottier-Gois, L.; Amarger, N. Typing of rhizobia by PCR DNA fingerprinting and PCR-restriction fragment length polymorphism analysis of chromosomal and symbiotic gene regions: Application to *Rhizobium leguminosarum* and its different biovars. *Appl. Environ. Microbiol.* **1996**, *62*, 2029–2036. [CrossRef] [PubMed]

83. Laguerre, G.; Allard, M.R.; Revoy, F.; Amarger, N. Rapid identification of rhizobia by restriction fragment length polymorphism analysis of PCR-amplified 16S rRNA Genes. *Appl. Environ. Microbiol.* **1994**, *60*, 56–63. [CrossRef] [PubMed]

84. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. *Mol. Biol. Evol.* **2016**, *33*, 1870–1874. [CrossRef] [PubMed]

85. Vinuesa, P.; Silva, C.; Werner, D.; Martínez-Romero, E. Population genetics and phylogenetic inference in bacterial molecular systematics: The roles of migration and recombination in Bradyrhizobium species cohesion and delineation. *Mol. Phylogenetics Evol.* **2005**, *34*, 29–54. [CrossRef]

86. Sijilmassi, B.; Filalimaltouf, A.; Boulahyaoui, H.; Kricha, A.; Boubekri, K.; Udupa, S.; Kumar, S.; Amri, A. Assessment of Genetic Diversity and Symbiotic Efficiency of Selected Rhizobia Strains Nodulating Lentil (*Lens culinaris* Medik.). *Plants* **2021**, *10*, 15. [CrossRef] [PubMed]

87. Laguerre, G.; Nour, S.M.; Macheret, V.; Sanjuan, J.; Drouin, P.; Amarger, N. Classification of rhizobia based on *nodC* and *nifH* gene analysis reveals a close phylogenetic relationship among *Phaseolus vulgaris* symbionts. *Microbiology* **2001**, *147*, 981–993. [CrossRef] [PubMed]

88. Sarita, S.; Sharma, P.K.; Priefer, U.B.; Prell, J. Direct amplification of rhizobial *nodC* sequences from soil total DNA and comparison to *nodC* diversity of root nodule isolates. *FEMS Microbiol. Ecol.* **2005**, *54*, 1–11. [CrossRef] [PubMed]

89. Higgins, D.G.; Jeanmougin, F.; Gibson, T.J.; Plewniak, F.; Thompson, J.D. *The CLUSTAL_X Windows Interface: Flexible Strategies for Multiple Sequence Alignment Aided by Quality Analysis Tools*; Oxford University Press: Oxford, UK, 1997.

90. Jost, L. Entropy and diversity. *Oikos* **2006**, *113*, 363–375. [CrossRef]

91. Braak, T.C.J.F.; Smilauer, P. *Canoco Reference Manual and User's Guide: Software for Ordination, Version 5.0*; Microcomputer Power: Ithaca, NY, USA, 2012; p. 2012.

92. Zhang, F.; Cui, Z.; Chen, X.; Ju, X.; Shen, J.; Chen, Q.; Liu, X.; Zhang, W.; Mi, G.; Fan, M.; et al. Chapter one—Integrated Nutrient Management for Food Security and Environmental Quality in China. In *Advances in Agronomy*; Sparks, D.L., Ed.; Elsevier: Amsterdam, The Netherlands, 2012; Volume 116, pp. 1–40.

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.