Previous Article | Next Article ![]()
Applied and Environmental Microbiology, October 2008, p. 5991-5996, Vol. 74, No. 19
0099-2240/08/$08.00+0 doi:10.1128/AEM.01320-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Angèle N'Zoue,4,5,6,7,8 and
Gisèle Laguerre2,3,4,5,6,7,8*
Microbial Genetics Laboratory, Department of Genetics and Plant Breeding, Institute of Agricultural Sciences, Banaras Hindu University, Varanasi 221 005, India,1 INRA, UMR 1229 Microbiologie et Géochimie des Sols, F-21065 Dijon, France,2 Université de Bourgogne, UMR 1229, F-21065 Dijon, France,3 IRD, UMR 113 Symbioses Tropicales et Méditerranéennes, F-34398 Montpellier, France,4 INRA, UMR 113, F-34398 Montpellier, France,5 CIRAD, UMR 113, F-34398 Montpellier, France,6 SupAgro, UMR 113, F-34398 Montpellier, France,7 Université de Montpellier II, UMR 113, F-34398 Montpellier, France8
Received 13 June 2008/ Accepted 28 July 2008
|
|
|---|
|
|
|---|
China was the first center of domestication of the soybean plant, about 4,000 years ago, and is the largest producer in Asia. However, although commercial cultivation of soybeans started in India in the 1960s, this crop was probably introduced to this country as soon as it was domesticated in China. Therefore, India is considered a secondary center of domestication of soybeans (2, 10). This country is the fifth producer of soybeans in the world today (8.3 Mt in 2006 [http://faostat.fao.org]). However, the average yield of about 1.1 t ha–1 is low compared to the world average of 1.8 t ha–1, and improving the crop performance is a major challenge for India (6). The effectiveness of symbiotic N2 fixation may be an important factor to take into consideration through successful management of symbiosis between soybean cultivars and native rhizobia.
Soybean-nodulating rhizobia are genetically diverse and are classified into different genera and species, with all the species described so far having been identified in China. The slow growers are distributed in three species of the Bradyrhizobium genus, namely, Bradyrhizobium japonicum (9), Bradyrhizobium liaoningense (30), and Bradyrhizobium elkanii (11). Fast growers belong to Sinorhizobium fredii and S. xinjiangense and also include other unclassified rhizobia (3, 5, 16). Soybean rhizobia with a variable generation time were classified into Mesorhizobium tianshanense (4). However, despite the importance of soybean cultivation in India, little is known about the genetic resources in native rhizobia. In particular, the genetic diversity of soybean rhizobia has not been described so far. Therefore, in this study we examine the genetic diversity of a core collection of 50 isolates from soybean nodules sampled in various fields in different agricultural-ecological-climatic regions of India, including the main cropping areas. Four loci were analyzed by restriction fragment length polymorphism analysis of PCR-amplified DNA fragments (PCR-RFLP), including the 16S rRNA gene, the intergenic spacer (IGS) between the 16S and 23S rRNA genes, and the symbiotic genes nifH and nodC. The IGS and housekeeping genes are currently used as markers for molecular systematics and for estimations of phylogenetic relationships among bradyrhizobia (20, 23, 26-28), with the 16S rRNA gene having few polymorphisms within the Bradyrhizobium genus (23, 27, 28). Genetic relationships of the soybean isolates to other rhizobia were therefore investigated by sequencing the IGS and the housekeeping genes dnaK, glnII, and recA. Phylogenetic analysis of nifH was also performed.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Origins of soybean rhizobial isolates and agricultural-ecological-climatic characteristics of sampled regions
|
|
View this table: [in a new window] |
TABLE 2. Distribution of rhizobia isolated from root nodules of soybean plants into genotypes
|
Nucleotide sequence accession numbers.
All sequences from soybean isolates have been deposited in the GenBank database under accession numbers EU357926 through EU357931 for the IGS, EU357921 through EU357925 for nifH, EU818928 through EU818931 for dnaK, EU818932 through EU818935 for glnII, and EU818936 through EU818939 for recA. The nifH sequences of the B. liaoningense, Bradyrhizobium canariense, and Bradyrhizobium yuanmingense type strains have also been deposited, under accession numbers EU818925, EU818926, and EU818927, respectively.
|
|
|---|
Geographical distribution of genotypes.
Three to six genotypes were identified in each region, and genotypes I and VII were detected in all geographic areas (Fig. 1). An exact test of sample differentiation based on haplotype frequencies according to geographic region was performed using Arlequin software (18). No differentiation was found (P = 0.2) according to the global test as well as in testing differentiation between all pairs of samples, except between the populations from regions A and C (P = 0.02). One genotype (VI) was predominant in three populations (40 to 60% of isolates). Because the sampling in each region was very small, diversity was presumably underestimated. However, it should be noticed that the three isolates sampled at the Chengalpattu site (region D) were classified into three predominant genotypes, I, III, and VII. This result suggests that the diversity of soybean bradyrhizobia might be especially high in this geographic area, which is not a traditional site for soybean cultivation, in contrast to the other sampled regions, which are considered soybean hubs of India.
![]() View larger version (20K): [in a new window] |
FIG. 1. Distribution of bradyrhizobial genotypes (I to VIII) defined by the combined data from PCR-RFLP analysis of the 16S rRNA gene, IGS, nifH, and nodC among agricultural-ecological-climatic regions of India.
|
70%. The three representatives of IGS types I, II, and III formed a cluster (97 to 100% similarity between the aligned parts of the sequences) with B. yuanmingense and other unclassified bradyrhizobia, including strain MSDJ G17, isolated from a soybean plant in Argentina (A. Hartmann, personal communication). The IGS type V and VI sequences were closely related (99 to 100% similarity) to the type strain of B. liaoningense and to strain USDA 135 from soybean, which was formerly assigned to B. japonicum prior to the description of B. liaoningense as a novel species. The IGS type IV sequence (strain SR83) and unclassified Bradyrhizobium strains isolated from Pachyrhizus erosus nodules (17) formed a tight cluster that could not be related to any described species (only 76 to 80% similarity between the IGS sequences) but which might have a closer relationship with Bradyrhizobium genospecies IV (93 to 94% similarity). The sequence of strain SR83 did not show high percent identities with those of the various B. japonicum serotypes and other bradyrhizobial accessions of the USDA ARS National Rhizobium Germplasm Collection (23, 24). Therefore, our results show that 74% of isolates (IGS types I to IV) could not be related to a species (or genospecies) that was described so far as including soybean symbionts.
![]() View larger version (25K): [in a new window] |
FIG. 2. Phylogenetic ML tree based on 938-bp alignment of nucleotide sequences of the IGS between the 16S and 23S rRNA genes. Only bootstrap probability values of 70% (over 100 replicates) are indicated at the branching points. The scale bar indicates the number of substitutions per site. Soybean isolates are shown in bold, with their IGS haplotypes indicated. The letter "T" indicates the type strain of the species. Accession numbers of sequences extracted from GenBank are as follows: PAC53, AY628088; type XIX, AY493859; LMG R16434T, AJ534605; MSDJ G17, AF338851; PAC40, AY628087; LMG 6138T, AJ279264; BC-P22, AY386706; LMG 21987T, AJ631967; USDA 110, BA000040; LMG 18230T, AJ279301; USDA 135, AF208511; LMG 11944, AJ279287; LMG 10726, AJ279281; EC550-1, AY628085; PAC43, AY628089; PAC18, AY628093; LMG 10677, AJ534597; LMG 10671, AJ279272; LMG 10705, AJ534592; LMG 6134T, AJ279308; LMG 11951, AJ534594; and LMG 8293, AJ279311.
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Phylogenetic ML tree based on 1,493-bp alignment of concatenated nucleotide sequences of dnaK (489 bp), glnII (519 bp), and recA (482 bp). Only bootstrap probability values of 70% (over 100 replicates) are indicated at the branching points. The scale bar indicates the number of substitutions per site. Soybean isolates are shown in bold, with their IGS haplotypes indicated. The letter "T" indicates the type strain of the species. Accession numbers of sequences extracted from GenBank are as follows: CCBAU 10071T, AY923039 (dnaK), AY386780 (glnII), and AY591566 (recA); LMG 18230T, AY923041 (dnaK), AY386775 (glnII), and AY591564 (recA); USDA 6T, AM168362 (dnaK), AF169582 (glnII), and AM168341 (recA); USDA 110, BA000040 (complete genome); LMG 21987T, AY923046 (dnaK), AB353733 (glnII), and AB353734 (recA); BTA-1T, AY923047 (dnaK), AY386765 (glnII), and AY591553 (recA); USDA 76T, AM168363 (dnaK), AY599117 (glnII), and AM168342 (recA); ORS 278, CU234118 (complete genome); and BTAI1, CP000495 (complete genome).
|
700 bp) of representatives of the three nifH haplotypes and of the type strains of B. liaoningense, B. canariense, and B. yuanmingense were determined. The three nifH haplotypes were placed in two distinct clusters (Fig. 4). nif types I and II were grouped with B. yuanmingense reference strains, including Chinese soybean isolates (e.g., GenBank accession number EU146010). As expected, nif type III sequences were very similar to those of the B. japonicum bv. glycinearum and B. liaoningense reference strains.
![]() View larger version (23K): [in a new window] |
FIG. 4. Phylogenetic ML tree based on 612-bp alignment of nucleotide sequences of the nifH gene. Only bootstrap probability values of 70% (over 100 replicates) are indicated at the branching points. The scale bar indicates the number of substitutions per site. Soybean isolates are shown in bold, with their nifH haplotypes indicated. The letter "T" indicates the type strain of the species. Accession numbers of sequences extracted from GenBank are as follows: CCBAU 83698, EU146010; CCBAU 23230, AY934871; CCBAU 61071, EU113228; USDA 110, BA000040; S 127, DQ485701; IFO 14791T, AB094963; CPAC 402, DQ485714; MAFF303099, BA000012; 3841, AM236084; and 1021, AE007235.
|
|
|
|---|
Only 26% of the Indian soybean isolates (genotypes VII and VIII) could be assigned to one valid species known to include soybean symbionts, i.e., B. liaoningense. Our current knowledge of the genetic diversity of rhizobia associated with soybeans refers to various studies in different continents. However, all of the species described so far could be found in China (32), which is considered the center of origin of the soybean plant. Therefore, assignment of a significant proportion of Indian soybean isolates to B. liaoningense was not surprising, since this species is widespread in Asia (19, 30, 32).
The finding that 36% of isolates (genotypes I to V) were closely related to B. yuanmingense was firstly unexpected, since spontaneous nodulation of soybean plants by this species has not been reported. However, in searching for sequence similarities in databases, we found highly similar sequences suggesting that soybean bradyrhizobia closely related to B. yuanmingense also exist in China and South America. The analysis of nifH and nodC corroborated the close relationship of these isolates to B. yuanmingense. This result argues for vertical acquisition of the symbiotic genes rather than relatively recent horizontal gene transfer of symbiotic genes from B. liaoningense or B. japonicum bv. glycinearum. The B. yuanmingense species was originally created based on the description of nodule isolates from Lespedeza spp. (33). Since then, several studies have reported nodulation of various legumes by B. yuanmingense-related rhizobia (8, 15, 17, 30). Studies of host specificity showed that these B. yuanmingense isolates were not capable of effective nodulation of soybean plants (8, 17, 31, 33). Two representative B. yuanmingense soybean isolates were evaluated on soybean plants and were effective in symbiotic N2 fixation (data not shown). Therefore, our results suggest the occurrence of different biovars within the B. yuanmingense species.
The third genomic group (38% of isolates; IGS type IV, genotype VI) characterized in this study could not be assigned to previously described Bradyrhizobium species but harbors the same symbiotic genetic information as B. japonicum bv. glycinearum and B. liaoningense, which might result from acquisition by horizontal transfer of a B. liaoningense or B. japonicum symbiotic island. Such an event was recently reported to have occurred in Brazilian soybean fields (1). IGS type IV showed significant similarity with unclassified bradyrhizobia from Pachyrhizus erosus plants native to Central and South America, but these were not able to nodulate soybean plants (17). We also evaluated one representative of IGS type IV on soybean plants and found that the strain was fully effective for N2 fixation (data not shown). These results again suggest the occurrence of different biovars of the same genomic species. Further phylogenic and taxonomic analyses are in progress to clarify the classification of soybean isolates grouped in genotype VI.
In conclusion, the present study provides the first analysis of the genetic diversity of soybean-nodulating Bradyrhizobium isolates from India, a major soybean producer. The results show that the diversity is wider than expected based on previous studies in various geographic areas and on the current taxonomy of soybean rhizobia. Notably, the diversity of the soybean symbionts appears to be conserved across the agricultural-ecological-climatic regions sampled. Further studies should investigate the link between genetic diversity and functional variability of the soybean-rhizobium symbiosis.
Chinnaswamy Appunu was supported by a research fellowship from the Indian Council of Scientific and Industrial Research—University Grant Commission and by a short-term training fellowship from the European Molecular Biology Organization. This work was supported by grant 45 from the French governmental organization Bureau des Ressources Génétiques.
Published ahead of print on 1 August 2008. ![]()
Present address: Division of Microbiology, M. S. Swaminathan Research Foundation, Taramani Institutional Area, Third Cross Road, Chennai 620113, India. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2010 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»