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

IRD, UMR 113, Laboratoire des Symbioses Tropicales et Méditerranéennes, Campus International de Baillarguet TA A-82/J, F-34398 Montpellier Cedex 5, France,1 Montpellier SupAgro, UMR 113, Laboratoire des Symbioses Tropicales et Méditerranéennes, Campus International de Baillarguet TA A-82/J, F-34398 Montpellier Cedex 5, France,2 INRA, USC1242, Laboratoire des Symbioses Tropicales et Méditerranéennes, Campus International de Baillarguet TA A-82/J, F-34398 Montpellier Cedex 5, France,3 CIRAD, UMR 113, Laboratoire des Symbioses Tropicales et Méditerranéennes, Campus International de Baillarguet TA A-82/J, F-34398 Montpellier Cedex 5, France,4 Université Montpellier II, UMR 113, Laboratoire des Symbioses Tropicales et Méditerranéennes, Campus International de Baillarguet TA A-82/J, F-34398 Montpellier Cedex 5, France5
Received 16 May 2008/ Accepted 14 July 2008
|
|
|---|
|
|
|---|
The genetic diversity of these two bacterial species, especially S. meliloti, has usually been described as quite important (2, 7, 11, 40). This bacterial diversity, assessed by phenotypic or genotypic analyses, has been reported to be influenced by several factors, including geographical location (32), soil factors (14), and plant species (19, 37). Furthermore, at an interspecies level, several studies on S. meliloti and S. medicae natural populations, trapped from various soils (using several Medicago species), reported either the balanced coexistence of the two bacterial species (3, 23) or the dominance of one of them (3, 31, 33, 43).
At a plant genotype level, host nodulation specificity could also influence the bacterial diversity in soil. Significant variations were detected among rhizobial populations associated with different plant lines of Vicia faba (38), Pisum sativum (10), or Medicago sativa (26). Symbiotic specificity studies among M. truncatula lines were based on single-strain inoculations. Some of them showed qualitative symbiotic polymorphism for nodulation (Nod+/Nod–) and nitrogen fixation (Fix+/Fix–) among a wide collection of M. truncatula fixed lines inoculated with diverse S. meliloti strains (35, 39). Other studies demonstrated that different plant-rhizobium combinations displayed different quantitative levels of effectiveness, described as the input that a bacterial nitrogen-fixing symbiosis makes to plant nitrogen metabolism, with S. meliloti (27) and S. medicae species (17). However, as a host plant usually interacts with several symbionts, these studies do not reflect possible quantitative variations of symbiotic effectiveness within a bacterial population encompassing diverse genotypes. Studying plant symbiotic specificity when put in contact with a single soil would reflect what occurs in a natural ecosystem compared to artificial single inoculations.
Nodulation competition among S. meliloti and S. medicae genotypes, at both intra- and interspecies levels, is another biotic factor which could also influence the genetic structure of bacterial populations recovered from soil. Strain competitiveness (defined as the ability to dominate nodulation) experiments were mostly performed to compare inoculated strains versus native bacteria, in order to enhance the agricultural yield by selecting the best strains for nodulation and nitrogen fixation efficiency (21, 42).
In the M. truncatula-Sinorhizobium symbiotic system, it is still unclear (3) if the plant cultivar and competition among bacterial genotypes acting at both inter- and intraspecies levels are true selective forces driving the partner choice in nitrogen fixation association. Furthermore, a study of the nitrogen fixation polymorphism of these bacteria within a single geographical site has never been performed. With the aim of improving our understanding of the bacterial population structure and diversity maintenance in this symbiotic association, we carried out a study whose objectives were (i) to explore the genetic diversity of symbiotic bacteria recovered from a single soil and trapped with four different M. truncatula lines, (ii) to test for significant symbiotic divergence among Sinorhizobium populations sampled with four different M. truncatula lines, (iii) to perform competition assays within and between the two bacterial species S. meliloti and S. medicae to evaluate strain competitiveness, and (iv) to compare the fixation efficiencies of single-strain inoculants and multistrain mixtures.
|
|
|---|
Aliquots of soil were put into contact with four M. truncatula fixed lines, originally sampled from different countries: (i) DZA315-16 (Algeria), (ii) F83-005 (France), (iii) PRT180-21 (Portugal), and (iv) the reference line, JemalongA17. These lines were chosen based on the data of Ronfort et al. (30), which showed large genetic distances among them using microsatellite data. After seed sterilization, five individuals per line were cultivated for 2 months on a soil aliquot, and nodules were collected from plant roots and stored at –80°C in 20% glycerol till they were used. From 9 to 14 nodules per plant were surface sterilized, and isolation of a single bacterial clone per nodule was performed.
DNA extraction and marker choice for genetic-diversity analyses.
Bacteria were grown in 20 ml yeast extract-mannitol liquid medium for at least 3 days at 27°C under agitation (41), and total genomic DNA was isolated (3).
The genome of S. meliloti Sm1021 comprises one chromosome (3.65 Mb) and two megaplasmids (pSymA and pSymB), whereas that of S. medicae strain WSM419 comprises in addition one small plasmid. We amplified five genomic regions evenly distributed along the genome and implicated in different bacterial functions. Markers are located either within genes or adjacent to them within an intergenic spacer (IGS). We used two markers involved in symbiosis: IGSNOD (on pSymA), framed by nodulation genes involved in symbiosis (nodE and nodG), and typA (on the chromosome), which is essential for symbiosis on M. truncatula JemalongA17 (20). We also used two paralogous genes, virB11 and avhB11, implicated in effector translocation and conjugation, respectively, in Agrobacterium tumefaciens. In this study, we named loci virB11 and avhB11 as described in A. tumefaciens. In the S. meliloti reference strain (Sm2011), a single gene, named virB11, has been found (on pSymA) with high similarity to avhB11 of A. tumefaciens. On the other hand, the S. medicae reference strain (WSM419) also contains a single locus, also named virB11, but with high similarity to virB11 of A. tumefaciens. Furthermore, we studied the 16S rRNA gene, which allows species discrimination between S. meliloti and S. medicae.
PCR-CE-SSCP analyses.
Allelic discrimination was performed by capillary electrophoresis-single-stranded conformation polymorphism (CE-SSCP) (36). For the PCRs used for the SSCP, the primers of each set were 5' labeled with two different fluorescent dyes, FAM (5-carboxyfluorescein) and HEX (5-carboxy-4',5',2',4',5',7'-hexachlorofluorescein) (Table 1), and we used PCR mixtures and cycling conditions without a touchdown program, as described previously (3). Precipitation of the PCR products was performed by adding 1/10 sodium acetate volume (sodium acetate-EDTA buffer [1.5 M sodium acetate, pH >8, and 250 mM EDTA]) and 4 volumes of 95% ethanol and centrifuging the mixture at 4,000 rpm and 10°C for 35 min. The pellet was washed in 100 µl 70% ethanol and then centrifuged again for 15 min. The DNA was resuspended in 50 µl sterile Tris-HCl (10 mM).
|
View this table: [in a new window] |
TABLE 1. Primers used in the study
|
Nucleotide sequencing and PCR-RFLP analyses of the 16S rRNA gene.
We sequenced each different CE-SSCP allele obtained. For similar alleles, we randomly selected and sequenced from 1 to 26 individuals each. PCR amplifications were performed as described above using unlabeled primers, except for virB11 and avhB11 (Table 1). Bands were extracted from 1% agarose gels and purified using the Qiaquick PCR purification kit (Qiagen). Sequencing reactions were performed by the Macrogen Company. All isolates were characterized by PCR-restriction fragment length polymorphism (RFLP) analysis of the 16S rRNA gene (Table 1) with the restriction enzyme RsaI, which differentiates the two Sinorhizobium species.
Single- and mixed-inoculant preparation.
We tested each of the 13 SSCP genotypes isolated from soil for effectiveness with the four plant lines and for competition for nodulation (using mixed-strain inoculants). We estimated the correlation between the optical density at 620 nm and cell numbers (determined by plate counts) for the 13 reference strains separately (data not shown). In the first step, we carried out single-strain inoculations. In the second step, two bacterial mixtures containing, respectively, the nine genotypes of S. meliloti and the four genotypes of S. medicae were prepared. In the third step, an interspecies bacterial competition was performed by mixing the dominant genotype for nodulation of each bacterial species within each plant line, according to results obtained from the intraspecies competition. Mixed inoculants were prepared by combining equal quantities (i.e., the same volume at the same optical density at 620 nm) of each single-strain inoculant. In all experiments, the inoculum was adjusted to approximately 5 x 106 cells ml–1.
Plant growth and data collection in the competition and single-strain inoculant experiments.
Plants were individually grown in tubes (215 by 20 mm) containing 20 ml of inclined sterile nutritive medium (15 g liter–1 of agar) and 50 ml of liquid sterile nutritive solution diluted four times. After sterilization and germination, each seedling was put into a tube and inoculated with 1.5 ml inoculum 2 days later. Negative control plants received only sterile water, and positive controls received 5 mM nitrogen [Ca(NO3)2]. The plants were grown for 7 weeks with a 16-h photoperiod (alternating 24°C and 16°C). Ten replicates were used per plant line and treatment. The total number of nodules and the aerial-shoot biomass on all plants of each experiment were obtained. We estimated plant fitness by the dry aerial biomass, which is strongly correlated with the fixation efficiency (18, 24, 39). Nodules were harvested in the competition plant assays, and one bacterial isolate was purified from each nodule. For the intraspecific-competition assay, bacterial isolates were genotyped as described above by PCR sequencing on IGSNOD for S. medicae and on the three loci IGSNOD, typA, and virB11 for S. meliloti, as these markers were sufficient to distinguish each genotype. In the interspecific-competition experiment, PCR-RFLP on the 16S rRNA gene was used for genotyping bacterial species.
Statistical analysis of genetic diversity and population differentiation.
Diversity indexes of S. meliloti and S. medicae bacterial populations obtained from trapping were calculated at each locus using two summary statistics, allelic richness and Nei's diversity index, using Genclone 2.0 (1). The number of polymorphic sites (S), the number of analyzed sites, the estimator of nucleotide diversity, and Watterson's estimate of the scaled mutation rate were computed at each locus. When applicable (i.e., S > 0), Tajima's D statistic was calculated and departure from the neutral evolutionary model was tested for each species at each locus using Arlequin 3.1 (12).
The genetic differentiation among bacterial populations isolated from each plant line was tested by exact tests of global genotypic differentiation using Genepop 3.4 (28). Significant departures from equal proportions of the two species in both trapping and interspecific-competition experiments were tested with a chi-square test using Statistica 6.
Plant fitness and nodulation ability analyses.
Based on the single-inoculation experiment results, we tested the effects of plant lines, bacterial species, bacterial genotypes nested within species, and the two plant line-bacterial species and plant line-bacterial genotype interactions using an analysis of variance performed with JMP 7.0.1 (SAS Institute Inc.). Based on the two intraspecific competitions (the strain mixture for each species), we used an analysis of variance to test the effects of plant lines, bacterial species, and their interaction. Least-significant-difference Fisher post hoc tests were performed to compare the plant biomasses induced by the mixtures to those of the single-strain inocula.
Nucleotide sequence accession numbers.
All sequences have been deposited in the GenBank database under the following accession numbers: IGSNOD, AM991101 to AM991108; typA, AM991109 to AM991113; virB11, AM991116, AM992458, AM992459, and AM993155; and avhB11, AM991114 and AM991115.
|
|
|---|
2 = 169.27; df = 7; P < 0.001).
![]() View larger version (31K): [in a new window] |
FIG. 1. Compositions of the bacterial populations in trapping and in the interspecific competition experiment.
|
|
View this table: [in a new window] |
TABLE 2. Genotypes of S. medicae and S. meliloti isolates from soil trappinga
|
Using all loci, 13 different multilocus genotypes were detected, four and nine, respectively, for S. medicae and S. meliloti (Table 2). At least for the bacterial populations isolated from plant lines F83-005 and PRT180-21, S. meliloti diversity appeared larger than that of S. medicae in terms of allelic richness and genotypic distribution. Nucleotide divergence among strains was low for each locus, especially for S. medicae (Table 3). Tajima's D value for the different loci of S. meliloti did not depart from a neutral hypothesis of evolution. In S. medicae, the most negative D value was found for IGSNOD, deviating significantly from neutrality.
|
View this table: [in a new window] |
TABLE 3. Genetic diversity of S. medicae and S. meliloti populations from soil trapping
|
|
View this table: [in a new window] |
TABLE 4. Percentages of genotypes of S. medicae and S. meliloti isolates from trapping and intraspecific competitions in four plant lines
|
Intraspecific competition with all species genotypes.
Each of the 13 different bacterial genotypes was tested independently for nodulation. Only the S. meliloti ml9 strain could not renodulate any of the four plant lines (Fig. 2). This strain might be an opportunist rhizobium (coinfecting a nodule), which might not harbor all the symbiotic genes.
![]() View larger version (49K): [in a new window] |
FIG. 2. Mean dry aerial biomass (± standard error) (a) and mean nodule number (± standard error) (b) obtained for each plant line-bacterial genotype combination. The four S. medicae (md1 to md4) and the eight S. meliloti (ml1 to ml8) strains were inoculated separately on each M. truncatula line. "S. medicae" and "S. meliloti" indicate the results of the two mixtures including all the strains for each species, and "interspecific" is the result of a mixture including one strain of each species (see the text). "control –" is the negative control (no inoculation), and "control +" is the positive control (5 mM nitrogen).
|
2 = 363.48; df = 15; P < 0.001). It is noteworthy that, conversely, the most frequent genotype, md4, had been recovered at a very low frequency in the trapping experiment.
Seven of the nine different S. meliloti genotypes isolated from the soil-trapping experiment were recovered from the 130 nodules genotyped after plant inoculation with an equal mixture of the nine strains (Table 4). Two genotypes (ml2 and ml9) were not reisolated in this experiment. The chi-square tests could not rule out the hypothesis of equal frequencies among strain genotypes in DZA315-16 (
2 = 8.69; df = 5; P = 0.12), whereas one genotype was found significantly more frequently in other plant lines (data not shown). Finally, the four S. meliloti populations showed highly significant genotypic differentiation (G-like test, P < 0.001).
Interspecific competition with two species genotypes.
Genotype md4 for S. medicae was systematically used and mixed with either S. meliloti genotype ml8 (for the F83-005 and JemalongA17 lines) or ml5 (for the PRT180-21 and DZA315-16 lines). One hundred and twenty-one nodules were genotyped. S. meliloti strains were recovered significantly more frequently than S. medicae strains in each of the four populations (
2 = 23.99; df = 7; P < 0.001) (Fig. 1).
Plant fitness statistical analyses.
Each plant line was inoculated with each single bacterial genotype isolated, with each intraspecific mixture, or with the mixture of the two dominant strains of each bacterial species. The means and variances of the aerial biomasses and nodule numbers are shown in Fig. 2. Measurements of DZA315-16 plants inoculated with genotype md4 are absent, due to technical problems.
(i) Single inoculants.
Direct estimation of plant line-bacterial genotype (species) interaction was not possible due to the nonequilibrated data we obtained because of the absence of the combination DZA315-16/md4. Therefore, we included in our data virtual data corresponding to this missing combination. We tested several virtual conservative measures (of the nodule number and aerial biomass), i.e., the mean of the md4 strain, the mean of the DZA315-16 line, and the mean of all S. medicae strains. The final results were independent of the choice of the mean (data not shown), suggesting the reliability of such an approach to compensate for missing data.
For both the dry aerial biomass and nodule number measures of all single genotypes, the plant line, bacterial species, bacterial genotype (species), and plant line genotype (species) interaction effects were highly significant. Plant line-bacterial species interaction was found to be barely significant for the aerial biomass and not significant for the nodule number (Table 5, all single genotypes).
|
View this table: [in a new window] |
TABLE 5. Analyses of variance for dry aerial biomasses and nodule numbers of plants
|
S. meliloti species induced a higher nodulation rate than S. medicae species (P = 0.01).
(ii) S. medicae mixture versus S. meliloti mixture effects.
When plant biomass means between S. medicae and S. meliloti mixtures were compared, a slight bacterial species effect was detected (P = 0.023) (Table 5, two intraspecific mixtures), with a higher biomass for plants inoculated with the mixed inoculant of S. medicae species.
(iii) Intraspecific mixtures versus single inoculants.
Within each species, dry aerial biomasses obtained with S. meliloti and S. medicae mixtures were either equivalent to or significantly higher than the mean of each single-strain inoculation (Fig. 2 and Table 5, S. medicae and S. meliloti). The S. medicae mixture and the md2 genotype induced similar plant biomasses, in contrast to the other three bacterial genotypes (Table 6). Furthermore, post hoc tests showed that the S. meliloti mixture induced greater aerial biomass than five genotypes (ml1, ml2, ml4, ml5, and ml7), whereas no significant difference was detected between the mixture and ml3 or ml6 for DZA315-16 and F83-005 and ml3 or ml8 for JemalongA17 (Table 6).
|
View this table: [in a new window] |
TABLE 6. Comparisons of plant aerial biomass inoculated with each single bacterial genotype and S. medicae or S. meliloti intraspecific mixtures
|
|
|
|---|
Nucleotidic polymorphism was observed for the two homologous genes implicated in translocation and conjugation, virB11 and avhB11, between the two species, but also within S. meliloti. Our results support the hypothesis that the common ancestor of the two Sinorhizobium species harbored the two copies virB11 and avhB11. As several authors (2, 5) have suggested that S. medicae might have emerged from a population of S. meliloti, the emerging species S. medicae would have kept one ortholog, whereas S. meliloti harbored either the two copies or only one. The role of each gene in S. meliloti and between the two species merits more study.
In the trapping, we observed a high predominance of S. medicae strains compared to S. meliloti. In previous literature (2), no marked symbiotic preference for S. medicae was described for any M. truncatula line. Thus, the proportions observed in our study certainly reflect a major disequilibrium between the two bacterial species in the soil.
Medicago species display a wide variability of specificity (4). Such variations imply that Medicago species might play a role in bacterial selection in the soil, as demonstrated by Jebara et al. (19). In our study, the proportions of the two bacterial species (S. meliloti and S. medicae) were slightly different among the four plant lines at both interspecific and intraspecific levels, suggesting that the plant lines might select their symbiotic partners differently. M. truncatula plant lines thus display differences in their specificities against bacterial populations, although the presence or absence of rare strains might influence this differentiation. Consequently, M. truncatula genetic diversity in a natural ecosystem might, at least in part, have an influence on rhizobial genetic diversity in the soil, because viable bacteria could be released into the plant rhizosphere (25).
Both intraspecies competition results and soil-trapping results showed one genotype highly predominant for S. medicae (but not the same genotype) and a large genetic diversity for S. meliloti. Thus, a strong competition might exist among the different S. medicae genotypes for nodulation, leading to one predominant genotype being recovered, whereas this interstrain competition is apparently much weaker in S. meliloti, for which several genotypes were recovered in similar frequencies. We have no clue to explain such proportion differences in the competition outputs between the two species, even if it is not an intrinsic specificity of each species, as Bailly et al. (3) did not observe this unequal distribution in similar trapping experiments. Moreover, in S. medicae, the most competitive strain (md4) in the competition experiment was not the genotype found to be dominant in soil trapping (md1). We suspect that the md1 genotype is quantitatively dominant in the soil used but is not the most competitive for plant nodulation under the tube conditions tested, because competition among strains seems to be dependent on the experimental culture conditions, as previously demonstrated in the Bradyrhizobium japonicum-soybean model (15).
At an interspecies level, the competition experiment resulted in a higher proportion of S. meliloti, whereas trapping from soil resulted in more than 90% S. medicae isolated. The initial unbalanced proportions obtained in soil trapping most certainly resulted from a "population effect" rather than a methodological bias or a much higher competitiveness of S. medicae in soil. At least in this local population, we may consider S. meliloti the best M. truncatula symbiont for nodulation. To generalize this hypothesis, other competition analyses should be done between the two bacterial species trapped from other natural populations.
In Sinorhizobium spp., several authors have shown variation in the ability to fix nitrogen among strains originating from different geographical areas for S. meliloti (22, 27) or for S. medicae (17). In this study, we demonstrated high polymorphism in terms of nitrogen fixation effectiveness and nodulation ability within a single soil, at both the interspecies and the intraspecies levels, even though weak genotypic and nucleotidic divergence were observed among strains.
Within each bacterial species, we detected a significant bacterial genotype-plant line interaction on the aerial biomass and nodule number data. Previous studies showed that the fixation effectiveness of three S. meliloti strains (27) and two S. medicae strains (17) was dependent on the M. truncatula lines. Our study suggests that such interaction is observed not only among strains sampled from various geographical locations, but also within a single site, and that this can lead to dramatically divergent phenotypes (Fig. 2a, strain STM 5472). Nitrogen fixation variation among plant lines had already been described, and ecotype-strain specificity is, at least partially, controlled by a single recessive allele (Mtsym6) for the plant partner (39), whereas succinoglycan oligosaccharide structure is involved, at least partially, in the bacterial partner specificity (34).
The bacterial mixture containing the four S. medicae genotypes gave the same plant aerial biomass as the genotype md2, which is the most efficient of the four genotypes. The effectiveness reached by the mixture could thus be due to this single strain (md2), even if it was barely found in nodules. However, genotype occurrence could be underestimated, as coinfected nodules, which may represent around 15% (13), were not detected. In S. meliloti, depending on the plant line, either several genotypes displayed the same effectiveness as the mixture or the effectiveness of the mixture was found to be intermediate among all the strains. By comparing the effectiveness of each genotype separately, S. meliloti appears to be a better nitrogen-fixing species than S. medicae. Conversely, the bacterial mixture encompassing all S. medicae genotypes induces a better aerial plant biomass than the mixture including all S. meliloti strains. These opposite results might be explained by synergistic effects between strains in mixtures, as suggested by Bromfield (6), leading to a fixation advantage for S. medicae populations. This beneficial effect is not always the rule, as Heath and Tiffin (17) observed an antagonistic effect with a mixture of strains. Discrimination between S. meliloti and S. medicae as the best nitrogen-fixing symbiont species is not trivial, considering that in soil, populations of bacteria, rather than single strains, interact with the plant.
The influence of the plant cover on the structure evolution of soil bacterial populations appears extremely complex. Fine symbiotic differences exist within M. truncatula lines and can be expressed even when they are put in contact with a single soil. We showed that the partner choice in the M. truncatula-Sinorhizobium association is controlled by the plant lines and by the bacterial competition for nodulation and not strictly by the fixation efficiencies of the bacterial genotypes. However, it is relevant to keep in mind that bacterial adaptation to the environment is essential for symbiosis.
We thank Lucette Mauré and Mélodie Bousquet, who provided technical help for strain isolation. We are very grateful to Christophe Mougel and Pierre Offre, who collected nodules from the soil-trapping experiment.
Published ahead of print on 25 July 2008. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»