Open access

Comparison of Soybean-Nodulating Bradyrhizobia Community Structures Along North Latitude Between Japan and USA

Written By

Yuichi Saeki and Sokichi Shiro

Submitted: 12 September 2013 Published: 29 January 2014

DOI: 10.5772/57165

From the Edited Volume

Advances in Biology and Ecology of Nitrogen Fixation

Edited by Takuji Ohyama

Chapter metrics overview

2,977 Chapter Downloads

View Full Metrics

1. Introduction

Soybean (Glycine max [L.] Merr.) establishes a symbiotic relationship by infection with soybean-nodulating bacteria and subsequent root nodule formation. Soybean acquires atmospheric nitrogen as ammonia through the symbiotic nitrogen fixation by soybean-nodulating rhizobia in the root nodules. The inoculation of soybean with bradyrhizobia that has high ability of nitrogen fixation is considered to be effective in increasing soybean production. However, the efficiency of the inoculum may be poor if the inoculum can not compete with indigenous soybean-nodulating rhizobia in the soil or can not establish an efficient symbiosis with the host plants because of the increased densities of indigenous rhizobia. To solve this problem, it is very important to understand the ecology of indigenous soybean-nodulating rhizobia in terms of their genetic diversity, community structure, geographical distribution, compatibility with the host soybean, and the ecosystems including environmental factors associated with the localization and dominance of the rhizobial strains in the soil. Therefore, analysis of the genetic diversity and field distribution of indigenous soybean-nodulating rhizobia is important to improve our understanding of rhizobial ecology as well as inoculation methodology under various environmental conditions. It is likely that the community of soybean-nodulating rhizobia vary from place to place because various wild soybeans are distributed and various soybean cultivars are cultivated in the northern to southern regions of the world. In Japan, Sawada et al. [1] isolated 85 indigenous soybean bradyrhizobia from soybean root nodules sampled from 46 soybean fields and the isolates were classified based on their serotype using rabbit antisera prepared against Bradyrhizobium USDA strains as antigens. Minamisawa et al. [2] also characterized 213 Japanese indigenous soybean bradyrhizobia isolated from six fields by analysis of their fingerprints with RSα, RSβ, nifDK and hupSL, and revealed diversity and endemism in their population structure. They suggested that bradyrhizobia might diversify in individual fields depending on the associated host plants and local soil conditions. Soybean-nodulating bradyrhizobia show physiological and genetic diversity and the bradyrhizobial community structures are constructed under the various environmental conditions. The major soybean-nodulating rhizobia that have been identified are Bradyrhizobium japonicum, Bradyrhizobium elkanii, and Sinorhizobium/Ensifer fredii [3-6]. Furthermore, additional species of soybean-nodulating rhizobia have been extensively discussed in the literature owing to the complexity of their taxonomical classification [7-11]. Soybean-nodulating bacteria are found over a wide region of the world, and their genetic diversity reflects geographical and climatic differences as well as host diversity.

In the host soybean, the genes related to nodulation, the Rj genes, are known as nodulation regulatory genes, and the Rj genotypes of rj1, Rj2, Rj3, Rj4 and non-Rj which lack these genetical phenotypes have been confirmed to exist in nature [12-17]. Specific rhizobial strains are incompatible with soybean cultivars harboring a particular Rj gene. In addition, indigenous soybean-nodulating rhizobia may show a preference for particular genotypes among the compatible genotypes, even among soybean plants cultivated in soil samples from the same field [18-23]. The ability of a soybean plant to host bradyrhizobia depends on the characteristics of Rj genes. Previous experimental results have also demonstrated that the community structure of soybean-nodulating bradyrhizobia depends on the host soybean Rj genotype and on the soybean cultivar, and it varies with cultivation temperature even in an identical soil sample [24, 25]. Since soybean cultivars harboring Rj genes are involved in the inhibition of effective nodulation by certain serogroups of rhizobia as well as in the preferential selection of appropriate rhizobia for nodulation, in the analysis of indigenous soybean-nodulating bacteria, it is important to use several kinds of Rj-genotypes of soybean cultivars for the isolation of rhizobia.

In our research group, Saeki et al. [26] investigated the genetic diversity and geographical distribution of indigenous soybean-nodulating rhizobia isolates from five sites in Japan (Hokkaido, Fukushima, Kyoto, Miyazaki, and Okinawa) by analyzing their restriction fragment length polymorphisms of polymerase chain reaction amplicons (PCR-RFLP) of the 16S-23S rRNA gene internal transcribed spacer (ITS) region, with 11 Bradyrhizobium strains that have USDA numbers as reference strains [27]. We reported that a geographical distribution of indigenous bradyrhizobia varied from northern to southern Japan. Furthermore, Saeki et al. [28] reported that the distribution of soybean-nodulating rhizobial niche in Japan was strongly correlated with latitude. The representative clusters of the isolated bradyrhizobia shifted from those of B. japonicum strains USDA 123, 110, and 6T to B. elkanii strain USDA 76T, moving from northern to southern Japan [29, 30].

Figure 1.

Schematic representation of the theme of this chapter. Are soybean-nodulating bradyrhizobial community structures similar between the same northern latitude in USA and Japan?

The United States of America (USA) is the world’s largest producer of soybeans. North latitudes between Japan and USA are similar and soybean cultivars are grown at latitudes similar to those of the soybean production areas of both countries. Understanding the geographical distribution of soybean-nodulating rhizobia in the USA therefore, would provide important knowledge about bradyrhizobial ecology and insights into appropriate inoculation techniques for soybean-nodulating rhizobia with high nitrogen fixation ability. We investigated the relationship between the genetic diversity of indigenous soybean-nodulating bradyrhizobia and their geographical distribution in the USA using nine communities of isolates from eight states [31]. We analyzed their genetic diversity and community structure by means of RFLP of PCR amplicons to target the 16S-23S rRNA gene ITS region, with 11 USDA Bradyrhizobium strains as reference strains. We also performed diversity analysis, multidimensional scaling analysis based on the mathematical dissimilarity index, and polar ordination analysis to describe the structure and geographical distribution of the soybean-nodulating bradyrhizobial community. The major clusters were B. japonicum belonging to the cluster Bj123 in the northern USA, and B. elkanii in the middle to southern regions. Dominance of bradyrhizobia in a community was generally larger for the cluster belonging to B. elkanii than for the cluster belonging to B. japonicum. The indigenous American soybean-nodulating bradyrhizobial community structure was also strongly correlated with latitude as well as that of Japan. Our results indicated that this community varies geographically. Adhikari et al. [32] revealed the genetic diversity of soybean-nodulating bradyrhizobia in relation to climate depending on altitude and soil properties, such as soil pH, in Nepal. Furthermore, our research group demonstrated that soybean-nodulating rhizobial communities differed greatly in nearby fields depending on whether the soil was acidic or basic, and it was found that S. fredii strains were dominant in the alkaline soils of Vietnam and Okinawa, Japan [23, 33]. These results suggest that a relationship exists between the geographic distribution of indigenous soybean-nodulating rhizobia, soil temperature (and its variations due to latitude and altitude), and soil pH. Therefore, knowledge of rhizobial ecology and biology in relation to numerous environmental factors and the environmental gradients is needed.

This chapter discusses the analysis of soybean-nodulating bradyrhizobial communities isolated from Rj-genotypes of soybean cultivars in relation to geographical differences including latitudinal gradients between Japan and USA (Fig.1). Analysis of RFLP of PCR amplicons of the 16S-23S rRNA gene ITS region and mathematical analysis of the PCR-RFLP results are demonstrated as possible approaches to the study of community diversity and ecosystem of soybean-nodulating rhizobia in relation to the rhizobial endemism in Japan and USA.

Advertisement

2. Classification of indigenous soybean-nodulating bradyrhizobia

Methods that are being developed and available for characterizing bradyrhizobial communities, include denaturing gradient gel electrophoresis (DGGE) analysis [34], terminal RFLP (T-RFLP) analysis [35], and automated ribosomal intergenic spacer analysis (ARISA) [36] using environmental DNA, and sequence polymorphisms targeting 16S rRNA gene (rDNA), the 16S–23S rDNA ITS region and other genomic and RNA sequences such as house-keeping genes and symbiotic functional genes [37-42]. In this section on experimental procedures, a relatively simple and reliable method for the study of indigenous soybean-nodulating bradyrhizobia isolated from nodules as described in detail previously [29] is demonstrated as one approach to the study of bradyrhizobial ecology.

2.1. Soil samples

Fresh soils for laboratory soybean cultivation were collected from some fields. We have analyzed soil samples from sixteen fields in Japan collected from 2004 to 2010 [30], and soil samples from nine experimental fields and farm fields in eight American states (US soils) in August 2010 [31] for isolation of soybean-nodulating bradyrhizobia (Table 1). These samples were weakly acidic-neutral soils collected from different regions along north latitude in these countries. At least, three soil samples were obtained from each field, to a depth of 10 cm, after removal of the surface litter, and the samples were homogenized to produce a single composite sample. Table 1 summarized the location, soil pH, and electric conductivity (EC) at these sites.

2.2. Isolation of indigenous soybean-nodulating bradyrhizobia

Since indigenous soybean-nodulating bacteria should be isolated from cultivars with different Rj genotypes, the soybean cultivars Akishirome, Bragg, or Orihime for non-Rj genotype, Bonminori, CNS, Hardee, or IAC-2 for Rj2Rj3 genotype, and Akisengoku, Fukuyutaka, or Hill for Rj4 genotype were cultivated in culture pots for 4 weeks in our laboratory. Soybean cultivars were planted in 1 L culture pots. The pots were filled with vermiculite and a 40% (v/v) N-free nutrient solution [21] and then autoclaved at 121°C for 20 min. The soybean seeds were sterilized by soaking for 30 s in 70% ethanol and 3 min in a dilute sodium hypochlorite solution (0.25% available chlorine), and then rinsed with sterile distilled water. A soil sample (2-3 g) was placed into the vermiculite at a depth of 2–3 cm, and the soybean seeds were sown into the soil. The plants were grown for 4 weeks in a growth chamber (day, 28°C for 16 h; night, 23°C for 8 h) with a weekly supply of sterile distilled water. After harvesting, the roots were washed thoroughly with tap water. The nodules were randomly collected and surface sterilized for 3 min in 70% ethanol and 30 min in a diluted sodium hypochlorite solution, and then rinsed with sterile distilled water. Each nodule was homogenized in sterile distilled water, and was streaked onto a yeast extract–mannitol agar (YMA) [43] plate and incubated for 5–7 days in the dark at 28˚C. To determine the genus of the isolates, a single colony was streaked onto YMA plates containing 0.002% (w/v) bromothymol blue (BTB) to determine whether the genus of the isolate was Bradyrhizobium or Sinorhizobium/Ensifer, based on change of BTB color [44], and incubated as described above. After incubation, each isolate was maintained on YMA slant medium at 4˚C for later analysis. As a negative control, soybean plants grown without soil were confirmed to form no nodules, eliminating the possibility of contamination with soybean-nodulating bacteria. Total number of soybean-nodulating bradyrhizobia isolated from each Rj-genotype soybean, for a sample soil, was considered as a soybean-nodulating rhizobial community in the soil sample.

Table 1.

Soil sample and the location of the sampling site, soil pH and EC in Japan and USA.

2.3. PCR-RFLP analysis of the 16S-23S rRNA gene ITS region

For DNA extraction, we cultured each isolate in 1.5 mL of HEPES-MES (HM) medium [45] supplemented with 0.1% L-arabinose [46] for 5 days at 28°C. Total DNA for the PCR template was extracted from the HM culture of the isolate as described by Hiraishi et al. [47]. Bacteria cells cultured in the HM medium were collected by centrifugation and washed with sterile distilled water. The cell pellet was suspended in 200 μL sterile distilled water. Then 40 μL of the suspension was mixed with 50 μL of BL buffer (40 mM Tris-HCl, 1% Tween 20, 0.5% Nonidet P-40, 1 mM EDTA, pH 8.0) and 10 μL of proteinase K (1 mg mL-1) and incubated at 60 °C for 30 min. Thereafter, the digested sample was incubated at 95 °C for 5 min. The sample was centrifuged at 15,000 × g for 10 min to remove undisrupted cells and large debris, and the supernatant was collected with a pipette. In the phylogenetic analysis, reference strains were used to classify the isolates, namely, eleven Bradyrhizobium USDA strains (B. japonicum USDA 4, 6T, 38, 110, 115, 123, 124, and 135, and B. elkanii USDA 46, 76T, and 94) were used in the RFLP analysis of the 16S-23S rRNA gene ITS region [27]. Total DNAs of the reference strains were extracted by means of the same procedure as that used for the isolates. In our study, PCR was carried out with Ex Taq DNA polymerase (TaKaRa Bio Inc., Otsu, Shiga, Japan). For ITS amplification, an ITS primer set (BraITS-F: 5’-GACTGGGGTGAAGTCGTAAC-3’, BraITS-R: 5’-ACGTCCTTCATCGCCTC-3’) designed for amplification of the 16S-23S rRNA gene ITS region of bradyrhizobia [26] was used for the PCR reaction. The PCR cycle consisted of a pre-run at 94 °C for 5 min, followed by 30 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 30 s, and extension at 72 °C for 1 min, with a final extension at 72 °C for 10 min. The RFLP analysis of the amplicon was performed with the restriction enzymes HaeIII, HhaI, MspI, and XspI (TaKaRa Bio Inc.). The restriction fragments were separated by a capillary electrophoresis apparatus (QIAxcel, Qiagen, UK) and submerged gel electrophoresis with ethidium bromide for visualization.

2.4. Cluster analysis

To construct a dendrogram based on the result of PCR-RFLP analysis for soybean-nodulating isolates, the fragment sizes on the electrophoresed gel were calculated by using appropriate fragment size markers and the fragment sizes deduced from the sequences of the reference strains. All reproducible fragments longer than 50 bp were used for the cluster analysis, but some irreproducible fragments were excluded. The genetic distance between pairs of isolates (D) was calculated using the following equation (1):

D= 1  [2NAB/(NA+NB)],E1

where NAB is the number of RFLP bands shared by strains A and B, and NA and NB are the numbers of restriction fragments of each of the two strains, respectively [48, 49]. The D values were calculated for all pairs of isolates, the cluster analysis was conducted by the unweighted pair group method using the arithmetic average (UPGMA), and a dendrogram was constructed with the PHYLIP software (J. Felsenstein, University of Washington, Seattle, WA, USA; http://evolution.genetics.washington.edu/phylip.html).

Community structures of soybean-nodulating bradyrhizobia from each soil samples were summarized into Table 2, as bradyrhizobial community structures from Japanese soils, and into Table 3, as those from US soils, respectively. Dendrograms constructed from these data are demonstrated in our reports [26, 30, 31, 33]. As reference strains, we used B. japonicum USDA strains 4, 6T, 38, 110, 115, 123, 124, and 135 and B. elkanii USDA strains 46, 76T, and 94. Although the topology can differ between dendrograms based on reference strain sequence data and those based on RFLP patterns, for analysis of soybean-nodulating bradyrhizobial communities, the dendrogram obtained as described is acceptable, because it is important to classify the numerous isolates into reliable clusters and/or groups based on reference strains. Furthermore, since the topology of a dendrogram can also differ depending on the electrophoretic resolution, fragment detection sensitivity, and analysis accuracy, it is important to classify isolates into appropriate RFLP patterns and clusters by using appropriate reference strains and restriction enzymes.

Table 2.

Cluster and the number of soybean-nodulating bradyrhizobia from Japan based on the report from Saeki et al. [30].

Upper number in each column indicates a whole number of isolates and lower indicates a number of isolates from non-Rj, Rj2Rj3, and Rj4, respectively.


Table 3.

Cluster and the number of soybean-nodulating bradyrhizobia from USA based on the report from Shiro et al. [31].

Upper number in each column indicates a whole number of isolates and lower indicates a number of isolates from non-Rj, Rj2Rj3, and Rj4, respectively.


2.5. Confirmation of nodule formation of isolate

Several representative isolates in each operational taxonomic unit (OTU) of the dendrogram were confirmed for their nodulation capability on host soybean by inoculation test. Each isolate was cultured in yeast extract-mannitol broth (YMB) culture [43] for 6 days at 28°C, and the cultures were then diluted with sterile distilled water to approximately 106 cells mL-1. The soybean seeds were sown into 500 mL prepared culture pots without soil, as described above, and inoculated with a 1 mL aliquot of each isolate per seed, with two or three replicates. We assessed nodule formation after 3 weeks in a growth chamber under the conditions described above.

Advertisement

3. Mathematical ecology analysis of soybean-nodulating bradyrhizobial communities

In previous section, we investigated the polymorphisms of 16S-23S rRNA gene ITS region of soybean-nodulating isolates and analyzed the community structures to elucidate bradyrhizobial ecology along latitude between Japan and USA. Several soybeans that are different in Rj-genotypes were cultivated for four weeks on the soil for isolation of soybean-nodulating bradyrhizobia. Polymorphisms of the isolates were detected by RFLP of 16S-23S rRNA gene ITS regions with HaeIII, HhaI, MspI, and XspI, with B. japonicum USDA 4, 6T, 38, 110, 115, 123, 124, 135, B. elkanii USDA 46, 76T, 94, as reference strains. As a result, 12 clusters from Japan and 9 clusters from USA were obtained (Tables 2 and 3). Isolates in this study were confirmed to belong to the genus Bradyrhizobium on the basis of their ability to form nodules on soybean roots, alkaline production on YMA medium containing BTB, and the PCR-RFLP analysis results. Among the clusters, Bj123, Bj110, Bj6, and Be76, from northern to southern regions, were isolated most frequently, followed by Be94, Bj115 from Japanese soils. In US soils, Bj123 was also dominant in the northern regions, whereas Be46, Be76, and Be94 were dominant in the central to southern regions, and Bj6 and Bj110 were moderately dominant in the central regions. In this section, mathematical approaches to the analysis of soybean-nodulating bradyrhizobial community diversity are demonstrated using the community structures of indigenous bradyrhizobia isolated from Japanese soils and US soils, presented in Tables 2 and 3.

3.1. Diversity indices analysis among field sites

To estimate differences among the diversities of the soybean-nodulating bradyrhizobial communities at different field sites, we used the Shannon-Wiener diversity index [50, 51]. Shannon-Wiener diversity index (H’) was calculated for each field site with equation (2):

H= PilnPi.E2

The Pi is the dominance of the isolate, expressed as (ni/N), where N and ni are the total number of isolates tested and the number of isolates belonging to a species, B. japonicum (Bj) or B. elkanii (Be), or a particular cluster at each field site, respectively. Thereafter, we calculated the alpha diversity (H’α), beta diversity (H’β), and gamma diversity (H’γ) to estimate the differences in the bradyrhizobial communities between pairs of soil samples [52, 53]. The H’α index represents a weighted average of the diversity indices of comparing two bradyrhizobial communities, the H’β index represents the differences between the bradyrhizobial communities from two soil samples (i.e., differences between sites), and the H’γ index represents the diversity of the total isolate communities from the two soil samples (Figure 2).

Figure 2.

Schematic representation of alpha, beta, and gamma diversities of soybean-nodulating rhizobial communities.

The relationship among these indices is expressed as the equation (3):

Hβ=Hγ Hα.E3

Thereafter, we estimated the differences among the compositions of the bradyrhizobial communities by comparing the ratio of beta to gamma diversity as the equation (4), taking into consideration the difference in gamma diversity in each pairwise comparison of bradyrhizobial communities.

Hβ/ Hγ.E4

3.2. Result of diversity indices analysis

At the northern sites, most isolates were classified into cluster Bj123, and the proportion of isolates in cluster Bj123 decreased southward. On the other hand, the frequency of isolation of B. elkanii increased southward, and the proportion of B. elkanii isolates was highest at the southern sites among the soil sample sites. As described above, Bj123, Bj110, Bj6, and Be76 clusters were isolated most frequently, from northern to southern regions in Japan. In US soils, Bj123 was also dominant in the northern regions, whereas Bj6 and Bj110 clusters were moderately dominant in the central regions, and Be46, Be76, and Be94 clusters were major in the central to southern regions in USA. In In comparing the northern and southern sites of both countries, beta diversities (H’β / H’γ) were higher than other comparisons (Tables 4 and 5, Figs. 3-6). These results mean that difference between bradyrhizobial community structures become larger in the case of comparison between northern and southern sites, and smaller in the case of similar latitudes. However, since varieties of beta diversities (H’β / H’γ) were larger in community structures of Japanese soils than those of US soils, tendency of beta diversity (H’β / H’γ), in comparison of bradyrhizobial community structures, is obscure as compared to those of US soil samples. This reason is due to the fact that some soil samples indicated high proportion of Bj110 cluster in Japanese bradyrhizobial community structures.

Table 4.

The ratio of beta to gamma diversity (H’β / H’γ) among soybean-nodulating bradyrhizobial communities in Japan.

Data in the upper triangle were based on the species (Bj and Be) data, and that in lower triangle were based on the cluster data.


Figure 3.

The ratio of beta to gamma diversity (H’β / H’γ) among pairs of soil sample sites based on the species (Bj and Be) data in Japan.

Figure 4.

The ratio of beta to gamma diversity (H’β / H’γ) among pairs of soil sample sites based on the cluster data in Japan.

Table 5.

The ratio of beta to gamma diversity (H’β / H’γ) among soybean-nodulating bradyrhizobial communities in USA.

Data in the upper triangle were based on the species (Bj and Be) data, and that in the lower triangle were based on the cluster data.


Figure 5.

The ratio of beta to gamma diversity (H’β / H’γ) among pairs of soil sample sites based on the species (Bj and Be) data in USA.

Figure 6.

The ratio of beta to gamma diversity (H’β / H’γ) among pairs of soil sample sites based on the cluster data in USA.

In the comparison between Japan and USA, bradyrhizobial community structures based on the species Bj and Be at similar geographical latitude indicated lower beta diversities and higher beta diversity at different geographical latitudes even in the comparison of different countries, Japan and USA (Table 6). On the other hand, in the comparison of beta diversities of community structures based on the cluster compositions, varieties of beta diversities (H’β / H’γ) were also larger due to high proportion of Bj110 cluster in some Japanese bradyrhizobial community structures, though tendency of results were similar to the result of H’β / H’γ based on the species data (Table 7).

Table 6.

The ratio of beta to gamma diversity (H’β / H’γ) among soybean-nodulating bradyrhizobial communities in the comparison between Japan and USA based on the species (Bj and Be) data.

3.3. Multidimensional scaling analysis

To describe the characteristics of the bradyrhizobial communities and the differences among field sample sites, we performed a multidimensional scaling (MDS) analysis based on the Bray-Curtis similarity measure. The Bray-Curtis similarity measure [54] has a robust monotonic relationship with ecological distance and a robust linear relationship with ecological distance up to large values of the distance. Thus, the Bray-Curtis similarity measure (BC) is one of the indices that best reflect the properties between communities [55]. The Bray-Curtis similarity measure was calculated using the following equation (5):

BCAB= [ |nAnB|] / [ (nA+nB)],E5

where BCAB is the dissimilarity between communities A and B, and nA and nB represent the total number of strains in B. japonicum (Bj) and B. elkanii (Be), or the number of strains in a particular cluster for communities A and B. The MDS analysis based on the Bray-Curtis similarity measure was conducted using the command “cmdscale” in the R software program version 2.15.1 (http://www.r-project.org/). Results of MDS analysis were indicated as two-dimensional (2-D) MDS analysis and three-dimensional (3-D) MDS analysis based on the species (Bj and Be) as independent variables or the clusters as independent variables.

Table 7.

The ratio of beta to gamma diversity (H’β / H’γ) among soybean-nodulating bradyrhizobial communities in the comparison between Japan and USA based on the cluster data.

3.4. Result of MDS analysis

As a result of MDS analysis of Bj and Be data set, the width of third dimension range was small in the 3-D MDS result, and it was considered that 2-D MDS result was able to explain their relationships among community structures (Figs. 7 and 8). As shown in Figure 7, communities mainly consisted of Bj were arranged at left part, and those mainly consisted of Be were arranged at right part of the MDS result. This result suggests that similar community compositions by Bj and Be exist in Japan and USA. Difference between Japan and USA is that frequency of isolation of Be of USA is higher than that of Japan. Thus, many bradyrhizobial communities of USA were positioned at right part of the 2-D MDS result.

In the case of characterization of community structures by MDS analysis based on Bj and Be as variables, 2-D MDS was reasonable to explain their relationships among communities because of only two variables, Bj and Be. However, in the case of characterization of community structure by MDS based on the cluster sets, 3-D MDS was more suitable to explain their relationships among communities (Figs. 9 and 10). This means that much number of clusters as variables are necessary for characterization of relationships among bradyrhizobial community structures. As described above, cluster Bj123, Bj110, Bj6, and Be76, from northern to southern regions in Japan, were isolated with high frequency, followed by Be94, Bj115 from Japanese soils. In US soils, Bj123 was also dominant in the northern regions, and Be46, Be76, and Be94 were dominant in the central to southern regions, and Bj6 and Bj110 were moderately dominant in the central regions. In the 3-D MDS result, similar coordinates were shown among bradyrhizobial communities isolated from soil samples, if the latitudes of the sample sites were near, without affecting the distance between Japan and USA. To inspect this result, a relationship between bradyrhizobial community structures and latitudes was investigated by a polar ordination of 3D-MDS coordinates of the community structures in the next analysis.

Figure 7.

Result of 2D-MDS analysis based on Bj and Be and the number of isolates as independent variables.

Figure 8.

Result of 3D-MDS analysis based on Bj and Be and the number of isolates as independent variables.

Figure 9.

Result of 2D-MDS analysis based on the clusters and the number of isolates as independent variables.

Figure 10.

Result of 3D-MDS analysis based on the clusters and the number of isolates as independent variables.

3.5. Polar ordination of community diversity and latitude

The 3-D MDS results were analyzed mathematically by comparing percentage differences between pairs of soybean-nodulating indigenous bradyrhizobial communities and by using a polar ordination analysis [56, 53] to examine the geographical distributions of soybean-nodulating bradyrhizobia between Japan and USA (Fig. 11). To determine the relative distances between the diversities based on the 3-D MDS plots of the communities in the 3-D Euclidean space as a function of latitude (°N), we calculated the Euclidean distances between the bradyrhizobial communities and poles. The MDS plot of Michigan was set as the northern pole, and that of Ishigaki was set as the southern pole, according to their latitudes of sample sites (Table 1). The distances between the MDS plots were calculated using the coordinates on the x-, y-, and z-axes as the Euclidean distance (Ed) using the equation (6):

EdAB= [(XAXB)2+(YAYB)2+(ZAZB)2]1/2,E6

where EdAB is the linear distance between communities A and B in the MDS plot and XA and XB, YA and YB, and ZA and ZB represent the x (axis 1 in 3-D MDS), y (axis 2 in 3-D MDS), and z (axis 3 in 3-D MDS) coordinates of communities A and B, respectively. The distances from each pole were converted into percent differences, D1 and D2, from the two polar communities (i.e., the Michigan and Ishigaki sites, which were considered to have a 100% difference). Simultaneous equations were constructed from the trigonometric figure using the Pythagorean theorem as described previously [28, 29].

Figure 11.

Schematic representation of polar ordination analysis.

Parameter x represents the polar difference (%) from the 0% pole (Ishigaki site) and is calculated as the following equation (7):

x= (L2+D12D22) / 2L,E7

where D1 and D2 are the percentage differences between a particular bradyrhizobial community and the communities at Ishigaki and at Michigan, respectively. Parameter L represents the 100% polar difference between the poles. Then, polar differences were calculated from the trigonometric diagram and plotted against the latitudinal difference between sites. The relationship of the polar ordination and the latitudes of the field sampling sites were estimated. This analysis was conducted for united data of sample soil site, based on the species (Bj and Be) and the clusters, for estimation of community structure distribution.

Figure 12.

Relationship between latitude and polar difference based on 3D-MDS from Bj-Be data set.

Figure 13.

Relationship between latitude and polar difference based on 3D-MDS from the cluster data set.

3.6. Result of polar ordination analysis

As shown in Figs. 12 and 13, the flora of indigenous soybean rhizobia changed gradually from north to south, and a distinctive flora was detected at each field site. The results of the polar ordination analysis showed that indigenous soybean-nodulating bradyrhizobial community structure was correlated with latitude. This result suggests that the community of indigenous bradyrhizobia at a particular geographical location might be affected by soil temperature associated with latitude or the diversity of the associated host plants acclimatized to that region’s climate. In our report, the higher dominance of localized B. elkanii strains in the soil was detected in USA [31]. In USA, the regression equations were shifted to lower than those of Japan. This is due to higher proportions of B. elkanii isolates in USA than those of Japan, though Bj123 was dominant in northern regions in both countries. The reason of dominance of B. elkanii in USA can be discussed based on other previous reports. Keyser et al. [57] examined the distribution of indigenous soybean-nodulating rhizobia in USA and found that serogroup of B. japonicum USDA123 was dominant in northern regions and B. elkanii in southern regions. Minamisawa et al. [58] investigated the preference of nodulation of soybean cultivars, a wild soybean progenitor (Glycine soja), and siratro (Macroptilium atropurpureum) by B. japonicum and B. elkanii strains. The result was that B. japonicum and B. elkanii preferentially nodulated G. max and M. atropurpureum, respectively, whereas both bradyrhizobial species formed nodules on G. soja with similar efficiency. Marr et al. [59] reported that Amphicarpaea bracteata performed nodule formation with B. japonicum and B. elkanii and performed nitrogen fixation with B. elkanii, though A. bracteata will not be the original host for B. elkanii. Furthermore, the microevolution and origins of Bradyrhizobium populations in eastern North America associated with soybean and native legumes (A. bracteata and Desmodium canadense) was investigated using genetic characterization by multilocus sequence typing of six core (housekeeping) gene sequences and two symbiotic gene sequences, and the results suggested that soybean-nodulating bacteria associated with native legumes represent a novel source of ecologically adapted bacteria for soybean inoculation [60]. Siratro is a major pasture legume that is cultivated in the tropics and subtropics, including parts of Australia, South and Central America, and some Pacific islands [61]. Additionally, A. bracteata is an annual legume that distributes widely in eastern North America [62, 63]. The presence of these legumes that have compatibility for nodulation with B. elkanii might be a reason to contribute strongly to the high dominance of B. elkanii in eastern North America.

In previous study, B. elkanii isolates were detected with high frequencies from southern regions of Japan [26, 33]. However, in those soil samples, isolation rates of B. elkanii were relatively low, and dominant isolates from fine-particle soils belonged to the Bj110 cluster (Table 2). In the United States, indigenous soybean-bradyrhizobia belonging to Bj123 cluster are dominant in northern regions, and Be clusters are dominant in central and southern regions [54, 31]. The Bj110 cluster is detected in central regions, but in association with few bradyrhizobia. These results suggest that soil chemical and/or physical properties determined by soil texture (e.g., silt versus clay) might affect indigenization and/or nodulation of soybean-nodulating bradyrhizobia. Fine-particle soils developed under paddy field conditions such as Gray Lowland soils might be suitable for indigenization of strains such as USDA110 strain under temperate climate regions. It has also been reported that soybean cultivation management practices affect the genetic variance of soybean-nodulating bacteria [64, 65]. Concerning the structure of soybean-nodulating bacterial communities, it is therefore necessary to consider soil types, cultivation conditions (e.g., sowing period), soil temperature, and soil management practices.

As a summary of these results, it is indicated that the composition of indigenous soybean-nodulating bradyrhizobial community was correlated with latitude in temperate regions. And it is suggested that the community structure of indigenous bradyrhizobia at a particular geographical location will be affected by soil temperature and/or the diversity of the associated host plants acclimatized to that region’s climate, with some exceptions in the case of fine-particle soils as discussed above, and alkaline-salinity soils, in which Bradyrhizobium serogroup 135 or Sinorhizobium/Ensifer fredii dominate [23, 33, 66-69].

Advertisement

4. Conclusion and future prospects

In this chapter, the RFLP analysis of the 16S–23S rRNA gene ITS region and mathematical analysis of the PCR-RFLP results were demonstrated as possible approaches to the study of community diversity and ecosystem of soybean-nodulating bradyrhizobia in relation to the rhizobial endemism in Japan and USA. As a result, generally, B. japonicum and B. elkanii generally indicated dominant existence from north to middle regions and from middle to south regions in both countries, respectively. Cluster Bj123 was dominant in northern regions, and cluster Be76 was dominant in southern regions in both countries. The bradyrhizobial community in USA was consisted from mainly Bj123, Bj110, Bj6, Be46, Be76, and Be94 clusters and diversity of Be cluster was higher than in Japan, and the bradyrhizobial community in Japan was consisted from mainly Bj123, Bj110, Bj6 and Be76 clusters and diversity of Bj clusters was higher than in USA. High coefficient of correlation was detected between community structures and north latitude. These results suggested that ecological niche of soybean-nodulating bradyrhizobial community will be detected along latitude, as a function of latitude and soil temperature.

The geographical distribution of bradyrhizobia along latitude reflects soil taxonomy such as zonal soils, the distribution of which on earth are affected by climate changes as a function of latitude. In contrast, the cluster of B. japonicum USDA 110 was dominant on fine-particle soils. The distribution of these strains reflects in part the distribution of soils such as intrazonal soils and is affected by water conditions and the oxidation-reduction potential in the soil. One of the reasons for the high occupancy of the Bj110 cluster in fine-particle soils might be strain capability for denitrification of bradyrhizobia. Though the end products of denitrification depend on the strain capability, B. japonicum strain USDA110 possesses a full set of functional denitrifying genes and reduces NO3 to N2 [70]. Furthermore, this strain evinces the denitrifying capability to reduce N2O surrounding the soybean root system [71]. Recently, Itakura et al. [72] demonstrated the mitigation of N2O emission from soils by inoculation of soils with B. japonicum USDA110 under field conditions. Therefore, utilization of useful bradyrhizobia that evince high N-fixing and full denitrifying capabilities is important not only for increasing yields but also for environmental conservation in agriculture and concerning the global warming.

Additionally, there are many reports on genetic diversity of soybean-nodulating rhizobia in subtropical-tropical regions. Appunu et al. [73] reported the genetic diversity of bradyrhizobia isolated from soybeans in India, and Jaiswal et al. [42] also reported the genetic diversity of soybean-nodulating rhizobia in India. They indicated the difference of soybean-nodulating rhizbial ecosystem from temperate regions and their broad host range. Yokoyama et al. [74] and Ando and Yokoyama [75] reported on Bradyrhizobium spp., which are different from B. japonicum and B. elkanii, based on genetic diversity of soybean-nodulating rhizobia in Thailand. Abaidoo et al. [76] reported heterogeneity of Bradyrhizobium spp. isolated from the new soybean cultivars in Africa as compared to bradyrhizobia from North American soybeans. These results suggest that diversity of soybean-nodulating rhizobia in subtropical-tropical regions and their ecosystems will be different from those in temperate regions. Further research on diversity and ecology of soybean-nodulating rhizobia in subtropical-tropical regions must be conducted for numerous environmental factors, containing soil types, climate conditions and soil managements to elucidate their ecology and to utilize their ecological traits for agriculture.

Because direct characterization of bradyrhizobial community structure in soil has so far been difficult, the characterization of rhizobial community structure has been limited with information coming only from analysis of soybean-nodulating rhizobial communities. It must be developed that the direct methods for the characterization of indigenous bradyrhizobial populations and community diversity in soils. Methods of characterizing indigenous rhizobial community structure from environmental DNA and use of media selective for bradyrhizobia from soil samples must therefore be developed to advance our understanding of indigenous rhizobial ecology and for construction of reliable models of soybean-nodulating rhizobial community structure.

References

  1. 1. Sawada Y, Miyashita K, Tanabe I, Kato K. Hup phenotype and serogroup identity of soybean nodulating bacteria isolated from Japanese soils. Soil Sci. Plant Nutr. 1989; 35: 281-288.
  2. 2. Minamisawa K, Nakatsuka Y, Isawa T. Diversity and field site variation of indigenous populations of soybean bradyrhizobia in Japan by fingerprints with repeated sequences RSα and RSβ. FEMS Microbiol. Ecol. 1999; 29: 171-178.
  3. 3. Jordan DC. Transfer of Rhizobium japonicum Buchanan 1980 to Bradyrhizobium gen. nov., a genus of slow-growing, root nodule bacteria from leguminous plants. Int. J. Syst. Bacteriol. 1982; 32: 136-139
  4. 4. Scholla MH, Elkan GH. Rhizobium fredii sp. nov., a fast-growing species that effectively nodulates soybeans. Int. J. Syst. Bacteriol. 1984; 34: 484-486.
  5. 5. Kuykendall LD, Saxena B, Cevine TE, Udell SE. Genetic diversity in Bradyrhizobium Jordan 1982 and a proposal for Bradyrhizobium elkanii sp. nov. Can. J. Microbiol. 1992; 38: 201-505.
  6. 6. Young JM. The genus name Ensifer Casida 1982 takes priority over Sinorhizobium Chen et al. 1988, and Sinorhizobium morelense Wang et al. 2002 is a later synonym of Ensifer adhaerens Casida 1982. Is the combination ‘Sinorhizobium adhaerens’ (Casida 1982) Willems et al. 2003 legitimate? Request for an opinion. Int. J. Syst. Evol. Microbiol. 2003; 53: 2107-2110.
  7. 7. Xu LM, Ge C, Cui Z, Li J, Fan H. Bradyrhizobium liaoningense sp. nov., isolated from the root nodules of soybeans. Int. J. Syst. Bacteriol. 1995; 45: 706-711.
  8. 8. Tan ZY, Xu XD, Wang ET, Gao JL, Martinez-Romero E, Chen WX. Phylogenetic and genetic relationships of Mesorhizobium tianshanense and related rhizobia. Int. J. Syst. Bacteriol. 1997; 47: 874-879.
  9. 9. Peng GX, Tan ZY, Wang ET, Reinhold-Hurek B, Chen WF, Chen WX. Identification of isolates from soybean nodules in Xinjiang region as Sinorhizobium xinjiangense and genetic differentiation of S. xinjiangense from Sinorhizobium fredii. Int. J. Syst. Evol. Microbiol. 2002; 52: 457-462.
  10. 10. Hungria M, Chueire LMO, Megías M, Lamrabet Y, Probanza A, Guttierrez-Mañero FJ, Campo RJ. Genetic diversity of indigenous tropical fast-growing rhizobia isolated from soybean nodules. Plant Soil 2006; 288: 343-356.
  11. 11. Vinuesa, P, Rojas-Jiménez K, Contreras-Moreira B, Mahna SK, Prasad BN, Moe H, Selvaraju SB, Thierfelder H, Werner D. Multilocus Sequence analysis for assessment of the biogeography and evolutionary genetics of four Bradyrhizobium species that nodulate soybeans on the Asiatic continent. Appl. Environ. Microbiol. 2008; 74: 6987-6996.
  12. 12. Williams LF, Lynch DL. Inheritance of a non-nodulating character in the soybean. Agron. J. 1954; 46: 28-29.
  13. 13. Caldwell BE. Inheritance of a strain-specific ineffective nodulation in soybeans, Crop Sci. 1966; 6: 427-428.
  14. 14. Caldwell BE, Hinson K, Johnson HW. A strain-specific ineffective nodulation reaction in the soybean Glycine max L. Merrill, Crop Sci. 1966; 6: 495-496.
  15. 15. Vest G. Rj3–A gene conditioning ineffective nodulation. Crop Sci., 10: 34–35. Vest, G. & Caldwell, B.E. (1972). Rj4–A gene conditioning ineffective nodulation in soybean. Crop Sci. 1970; 12: 692-693.
  16. 16. Vest G, Caldwell BE. Rj4–A gene conditioning ineffective nodulation in soybean. Crop Sci. 1972; 12: 692–693.
  17. 17. Devine TE, Kuykendall LD. Host genetic control of symbiosis in soybean (Glycine max L.), Plant Soil 1996; 186: 173-187.
  18. 18. Ishizuka J, Suemasu Y, Mizogami K. Preference of Rj-soybean cultivars for Bradyrhizobium japonicum for nodulation. Soil Sci. Plant Nutr. 1991; 37: 15-21.
  19. 19. Ishizuka J, Kim SD, Hussain AKM, Yamakawa T. Soybean preference for Bradyrhizobium japonicum for nodulation. Isolation of Rj2Rj4-lines from the cross of soybean cvs. IAC-2 (Rj2) and Hill (Rj4). Soil Sci. Plant Nutr. 1993; 39: 79-86.
  20. 20. Yamakawa T, Eriguchi M, Hussain AKMA, Ishizuka J. Soybean preference for Bradyrhizobium japonicum for nodulation. Nodulation by Rj2Rj3Rj4-genotypes isolated from the progenies of a cross between soybean cvs. IAC-2 (Rj2Rj3) and Hill (Rj4). Soil Sci. Plant Nutr. 1999; 45: 461-469.
  21. 21. Saeki Y, Akagi I, Takaki H, Nagatomo Y. Diversity of indigenous Bradyrhizobium strains isolated from three different Rj-soybean cultivars in terms of randomly amplified polymorphic DNA and intrinsic antibiotic resistance. Soil Sci. Plant Nutr. 2000; 46: 917-926.
  22. 22. Yamakawa T, Hussain AKMA, Ishizuka J. Soybean preference for Bradyrhizobium japonicum for nodulation. Occupation of serogroup USDA110 in nodules of soybean plants harboring various Rj-genes grown in a field. Soil Sci. Plant Nutr. 2003; 49: 835-841.
  23. 23. Saeki Y, Kaneko A, Hara T, Suzuki K, Yamakawa T, Nguyen MT, Nagatomo Y, Akao, S. Phylogenetic analysis of soybean-nodulating rhizobia isolated from alkaline soils in Vietnam. Soil Sci. Plant Nutr. 2005; 51: 1043-1052.
  24. 24. Minami M, Yamakawa T, Yamamoto A, Akao S, Saeki Y. Estimation of nodulation tendency among Rj-genotype soybeans using the bradyrhizobial community isolated from an Andosol. Soil Sci. Plant Nutr. 2009; 55: 65-72.
  25. 25. Shiro S, Yamamoto A, Umehara Y, Hayashi M, Yoshida N, Nishiwaki A, Yamakawa T, Saeki Y. Effect of Rj genotype and cultivation temperature on the community structure of soybean-nodulating bradyrhizobia. Appl. Environ. Microbiol. 2012; 78: 1243-1250.
  26. 26. Saeki Y, Aimi N, Tsukamoto S, Yamakawa T, Nagatomo Y, Akao S. Diversity and geographical distribution of indigenous soybean-nodulating bradyrhizobia in Japan. Soil Sci. Plant Nutr. 2006; 52: 418-426.
  27. 27. Saeki Y, Aimi N, Hashimoto M, Tsukamoto S, Kaneko A, Yoshida N, Nagatomo Y, Akao S. Grouping of Bradyrhizobium USDA strains by sequence analysis of 16S rDNA and 16S–23S rDNA internal transcribed spacer region. Soil Sci. Plant Nutr. 2004; 50: 517-525.
  28. 28. Saeki Y, Minami M, Yamamoto A, Akao S. Estimation of the bacterial community diversity of soybean-nodulating bradyrhizobia isolated from Rj genotype soybeans. Soil Sci. Plant Nutr. 2008; 54: 718-724.
  29. 29. Saeki Y. Characterization of soybean-nodulating rhizobial communities and diversity. In: Aleksandra S. (ed.) Soybean-molecular aspects of breeding. Rijeka: Intech; 2011. p163-184.
  30. 30. Saeki Y, Shiro S, Tajima T, Yamamoto A, Sameshima-Saito R, Sato T, Yamakawa T. Mathematical ecology analysis of geographical distribution of soybean-nodulating bradyrhizobia in Japan. Microbes Environ. 2013; in press.
  31. 31. Shiro S, Matsuura S, Saiki R, Sigua GC, Yamamoto A, Umehara Y, Hayashi M, Saeki Y. Genetic diversity and geographical distribution of indigenous soybean-nodulating bradyrhizobia in the United States. Appl. Environ. Microbiol. 2013; 79, 3610-3617.
  32. 32. Adhikari D, Kaneko M, Itoh K, Suyama K, Pokharel BB, Gaihre YK. Genetic diversity of soybean-nodulating rhizobia in Nepal in relation to climate and soil properties. Plant Soil 2012; 357: 131-145.
  33. 33. Suzuki K, Oguro H, Yamakawa T, Yamamoto A, Akao S, Saeki Y. Diversity and distribution of indigenous soybean- nodulating rhizobia in the Okinawa islands, Japan. Soil Sci. Plant Nutr. 2008; 54: 237-246.
  34. 34. Saeki Y, Ozumi S, Yamamoto A, Umehara Y, Hayashi M, Sigua GC. Changes in population occupancy of bradyrhizobia under different temperature regimes. Microbes Environ. 2010; 25: 309-312.
  35. 35. Wang Y, Wiatrowski HA, John R, Lin C-C, Young LY, Kerkhof LJ, Yee N, Barkay T. Impact of mercury on denitrification and denitrifying microbial communities in nitrate enrichments of subsurface sediments. Biodegradation 2013; 24: 33–46
  36. 36. Ikeda S, Rallos LEE, Okubo T, Eda S, Inaba S, Mitsui H, Minamisawa K. Microbial community analysis of field-grown soybeans with different nodulation phenotypes. Appl. Environ. Microbiol. 2008; 74: 5704-5709
  37. 37. Han TX, Wang ET, Han LL, Chen WF, Sui XH, Chen WX. Molecular diversity and phylogeny of rhizobia associated with wild legumes native to Xinjiang, China. Syst. Appl. Microbiol. 2008; 31: 287-301.
  38. 38. Han LL, Wang ET, Han TX, Liu J, Sui XH, Chen WF, Chen WX. Unique community structure and biogeography of soybean rhizobia in the saline-alkaline soils of Xinjiang, China. Plant Soil 2009; 324: 291-305.
  39. 39. Li QQ, Wang ET, Zhang YZ, Zhang YM, Tian CF, Sui XH, Chen WF, Chen WX. Diversity and biogeography of rhizobia isolated from root nodules of Glycine max grown in Hebei province, China. Microb. Ecol. 2011; 61: 917-931.
  40. 40. Wu LJ, Wang HQ, Wang ET, Chen WX, Tian CF. Genetic diversity of nodulating and non-nodulating thizobia associated with wild soybean (Glycine soja Siab. & Zucc.) in different ecoregions of China. FEMS Microbiol. Ecol. 2011; 76: 439-450.
  41. 41. Zhang YM, Li Y Jr, Chen WF, Wang ET, Tian CF, Li QQ, Zhang YZ, Sui XH, Chen WX. Biodiversity and biogeography of rhizobia associated with soybean plants grown in the north China plain. Appl. Environ. Microbiol. 2011; 77: 6331-6342.
  42. 42. Jaiswal SK, Anand A, Dhar B, Vaishampayan A. Genotypic characterization of phage-typed indigenous soybean bradyrhizobia and their host range symbiotic effectiveness. Microb. Ecol. 2012; 63: 116-126.
  43. 43. Vincent JM. A Manual for the Practical Study of the Root-Nodule Bacteria. International Biological Program, Blackwell Scientific, Oxford; 1970.
  44. 44. Keyser HH, Bohlool BB, Hu TS, Weber DF. Fast-growing rhizobia isolated from root nodules of soybean. Science 1982; 215: 1631-1632.
  45. 45. Cole MA, Elkan GH. Transmissible resistance to penicillin G, neomycin, and chloramphenicol in Rhizobium japonicum. Antimicrobial Agents & Chemofhempy 1973; 4: 248-253.
  46. 46. Sameshima R, Isawa T, Sadowsky MJ, Hamada T, Kasai H, Shutsrirung A, Mitsui H, Minamisawa K. Phylogeny and distribution of extra-slow-growing Bradyrhizobium japonicum harboring high copy numbers of RSα, RSβ and IS1631. FEMS Microbiol. Ecol. 2003; 44: 191-202.
  47. 47. Hiraishi A, Kamagata Y, Nakamura K. Polymerase chain reaction amplification and restriction fragment length polymorphism analysis of 16S rRNA genes from Methanogens. J. Ferment. Bioengineer. 1995; 79: 523-529.
  48. 48. Nei M, Li WH. Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc. Natl Acad. Sci. USA 1979; 79: 5269-5273.
  49. 49. Sakai M, Futamata H, Kim J-S, Matsuguchi T. Effect of soil salinity on population structure of fluorescent pseudomonads in spinach rhizosphere. Soil Sci. Plant Nutr. 1998; 44: 701-705.
  50. 50. MacArthur RH. Patterns of species diversity. Biol. Rev. 1965; 40: 510-533.
  51. 51. Pielou EC. Ecological diversity and its measurement. In: An Introduction to Mathematical Ecology. New York: Wiley/Interscience; 1969. p221-235.
  52. 52. Whittaker RH. Evolution and measurement of species diversity. Taxon 1972; 21: 213-251.
  53. 53. Bray JR, Curtis TJ. Anordinationoftheuplandforestcommunities of southern Wisconsin. Ecol. Monogr. 1957; 27: 325-349.
  54. 54. Faith PD, Minchin PR, Belbin L. Compositional dissimilarity as a robust measure of ecological distance. Vegetatio 1987; 69: 57-68.
  55. 55. Whittaker RH. Gradient analysis of vegetation. Biol. Rev. 1967; 49: 207-264.
  56. 56. Kobayashi S. Multivariate Analysis of Biological Communities. (in Japanese), Tokyo: Soju Shobo; 1995.
  57. 57. Keyser HH, Weber DF, Uratsu SL. Rhizobium japonicum serogroup and hydrogenase phenotype distribution in 12 states. Appl. Environ. Microbiol. 1984; 47: 613-615.
  58. 58. Minamisawa K, Onodera S, Tanimura Y, Kobayashi N, Yuhashi K, Kubota M. Preferential nodulation of Glycine max, Glycine soja and Macroptilium atropurpureum by two Bradyrhizobium species japonicum and elkanii. FEMS Microbiol. Ecol. 1997; 24: 49-56.
  59. 59. Marr DL, Devine TE, Parker MA. Nodulation restrictive genotypes of Glycine and Amphicarpaea: a comparative analysis. Plant Soil 1997; 189: 181-188.
  60. 60. Tang J, Bromfield ESP, Rodrigue N, Cloutier S, Tambong JT. Microevolution of symbiotic Bradyrhizobium populations associated with soybeans in east North America. Ecol. Evol. 2012; 2: 2943-2961.
  61. 61. Shaw HN, Whiteman PC. Siratro—a success story in breeding a tropical pasture legume. Trop. Grassl. 1977; 11: 7-14.
  62. 62. Parker MA. Local genetic differentiation for disease resistance in a selfing annual. Biol. J. Linn. Soc. 1991; 42: 337-349.
  63. 63. Turner BL, Fearing OS. A taxonomic study of the genus Amphicarpaea (Leguminosae). Southwest. Nat. 1964; 9: 207-218.
  64. 64. Bizarro MJ, Giongo A, Vargas LK, Roesch LFW, Gano KA, de Sá ELS, Passaglia LMP, Selbach PA. Genetic variability of soybean bradyrhizobia populations under different soil managements. Biol. Fertil. Soils 2011; 47: 357-362.
  65. 65. Grossman JM, Schipanski ME, Sooksanguan T, Seehaver S, Drinkwater LE. Diversity of rhizobia in soybean [Glycine max (Vinton)]nodules varies under organic and conventional management. Appl. Soil Ecol. 2011; 50: 14-20.
  66. 66. Damirgi SM, Frederick LR, Anderson IC. Serogroups of Rhizobium japonicum in soybean nodules as affected by soil types, Agron. J. 1967; 59: 10-12.
  67. 67. Ham GE, Frederick LR, Anderson IC. Serogroups of Rhizobium japonicum in soybean nodules sampled in Iowa. Agron. J. 1971; 63: 69-72.
  68. 68. Kowalski M, Ham GE, Frederick LR, Anderson IC. Relationship between strains of Rhizobium japonicum and their bacteriophages from soil and nodules of field- grown soybeans. Soil Sci. 1974; 118: 221-228.
  69. 69. Gross DC, Vidaver AK, Klucas RV. Plasmids, biological properties and efficacy of nitrogen fixation in Rhizobium japonicum strains indigenous to alkaline soils. J. Gen. Microbiol. 1979; 114: 257-266.
  70. 70. Sameshima-Saito R, Chiva K, Minamisawa K. Correlation of denitrifying capability with the existence of nap, nir, nor and nos genes in diverse strains of soybean bradyrhizobia. Microbes Environ. 2006; 21: 174-184.
  71. 71. Sameshima-Saito R, Chiba K, Hirayama J, Itakura M, Mitsui H, Eda S, Minamisawa K. Symbiotic Bradyrhizobium japonicum reduces N2O surrounding the soybean root system via nitrous oxide reductase. Appl. Environ. Microbiol. 2006; 72: 2526-2532.
  72. 72. Itakura M, Uchida Y, Akiyama H, Takada-Hoshino Y, Shimomura Y, Morimoto S, Tago K, Wang Y, Hayakawa C, Uetake Y, Sánchez C, Eda S, Hayatsu M, Minamisawa K. Mitigation of nitrous oxide emissions from soils by Bradyrhizobium japonicum inoculation. Nature Climate Change 2013; 3: 208-212.
  73. 73. Appunu C, N’Zoue A, Laguerre G. Genetic diversity of native bradyrhizobia isolated from soybeans (Glycine max L.) in different agricultural-ecological-climatic regions of India. Appl. Environ. Microbiol. 2008; 74: 5991-5996.
  74. 74. Yokoyama T, Ando S, Murakami T, Imai H. Genetic variability of the common nod gene in soybean bradyrhizobia isolated in Thailand and Japan. Can. J. Microbiol. 1996; 42: 1209-1218.
  75. 75. Ando S, Yokoyama T. Phylogenetic analyses of Bradyrhizobium strains nodulating soybean (Glycine max) in Thailand with reference to the USDA strains of Bradyrhizobium. Can. J. Microbiol. 1999; 45: 639-645.
  76. 76. Abaidoo RC, Keyser HH, Singleton PW, Borthakur D. Bradyrhizobium spp. (TGx) isolates nodulating the new soybean cultivars in Africa are diverse and distinct from bradyrhizobia that nodulate North American soybeans. Int. J. Syst. Evol. Microbiol. 2000; 50: 225-234.

Written By

Yuichi Saeki and Sokichi Shiro

Submitted: 12 September 2013 Published: 29 January 2014