Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2003, p. 884-893, Vol. 69, No. 2
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.2.884-893.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Laboratorio de Evolución Molecular y Experimental, Instituto de Ecología, Universidad Nacional Autónoma de México, México D. F. 04510,1 Programa de Ecología Molecular y Microbiana, Centro de Investigación sobre Fijación de Nitrógeno, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, Mexico2
Received 22 July 2002/ Accepted 14 November 2002
|
|
|---|
|
|
|---|
In Rhizobium a large proportion of the genome is composed of plasmids (16), which contribute significantly to the ecological fitness and symbiotic performance of rhizobial strains (6, 8, 16). The symbiotic genes for nodulation (nod) and nitrogen fixation (nif, fix) are located on the symbiotic plasmid (pSym). The pSym plasmids of biovar phaseoli strains that nodulate P. vulgaris have multiple copies of the nitrogenase reductase gene (nifH) and confer a restricted host range (28). In contrast, R. tropici and R. gallicum bv. gallicum pSym plasmids carry a single nifH copy and confer a broader host range that includes Leucaena spp. (3, 9). Segovia and collaborators (41) suggested that R. leguminosarum bv. phaseoli is the result of pSym transfer from R. etli bv. phaseoli in historic times. The R. gallicum and R. giardinii bv. phaseoli strains that nodulate P. vulgaris harbor pSym plasmids similar to those found in R. leguminosarum bv. phaseoli and R. etli bv. phaseoli, which led Amarger and colleagues (3) to propose that the acquisition of R. etli bv. phaseoli pSym took place via R. leguminosarum strains. The presence of viable R. etli bv. phaseoli strains on the testa of P. vulgaris seeds has been demonstrated (36). This finding could explain the geographical spread of R. etli bv. phaseoli and provides a scenario for the lateral transfer of pSym among indigenous Rhizobium species in historic and recent times along with the introduction of bean crops worldwide.
Besides the strong evidence for interspecies pSym transfer, population level analyses of Rhizobium species have also revealed the existence of lateral transfer of pSym within species in agricultural fields and pastures (25, 27, 40, 60). However, Wernegreen and colleagues (57) found a limited pattern of pSym transfer (48) in R. leguminosarum bv. trifolii populations associated with native Trifolium species growing in mountain meadows of California.
The effects of agricultural practices and plant host domestication on the genetic structure of Rhizobium populations have scarcely been addressed (30, 47). The domestication of wild P. vulgaris plants began around 4,000 years ago in Mesoamerica and the Andean region of South America (22). Beans were probably codomesticated with maize in Mesoamerica, since these two crops are grown in association in a traditional agrosystem called milpa (29). The milpa system is a prehispanic cultivation method, in which beans are intercropped with maize and squash, together with diverse other plant species that are locally used for medicinal and nutritional purposes (47). This cultural practice promotes bean nitrogen fixation, and its advantages have been recognized (5, 29). In a previous study we compared the genetic structures of R. etli bv. phaseoli strains associated with beans under different degrees of domestication (47). We found that the bacteria associated with milpa beans had a genetic structure intermediate between the genetic structures of the bacteria isolated from wild beans and monocultured beans (47). In a more recent study, we analyzed the genetic structure of a Rhizobium population obtained from nodules of P. vulgaris and Phaseolus coccineus plants from milpa plots in San Miguel, Puebla, Mexico, with the aim of understanding its spatial variation (45). To do this, six plots were sampled in a single year, and a hierarchical analysis of multilocus enzyme electrophoresis (MLEE) data revealed (i) the coexistence of two distantly related genetic groups, designated genetic divisions I and III, and (ii) the existence of numerically dominant genotypes within each division; in addition, linkage disequilibrium analyses indicated that recombination is frequent within each genetic division but not between the divisions, which led us to propose a reticulated and epidemic genetic structure (32, 45). In the present study, P. vulgaris plants from one of the previously studied milpa plots were sampled over three consecutive years to determine (i) the temporal stability of the population genetic structure, (ii) the structure and dynamics of the plasmidic compartment, and (iii) the taxonomic affiliations of the genetic groups identified.
|
|
|---|
MLEE.
Cell lysates of the isolates were obtained as previously described (45), and electrophoresis was performed on cellulose acetate membranes (20, 45). The following six enzymes were assessed: isocitrate dehydrogenase (EC 1.1.1.42), peptidase (EC 3.4.13), phosphoglucomutase (EC 5.4.2.2), glucose-6-phosphate dehydrogenase (EC 1.1.1.49), xanthine dehydrogenase (EC 1.1.1.204), and malate dehydrogenase (EC 1.1.1.37). Isocitrate dehydrogenase, peptidase, and phosphoglucomutase each exhibited one band of activity, glucose-6-phosphate dehydrogenase and xanthine dehydrogenase each exhibited two bands, and malate dehydrogenase exhibited three bands, which yielded a total of 10 loci for analysis.
Genetic diversity and cluster analysis.
Distinctive mobility variants of each enzyme, numbered in order of decreasing anodal mobility, were considered alleles at the corresponding locus (43). In the case of enzymes that produced more than one band, each band was considered a locus. The absence of enzyme activity was scored as a null allele and was treated as an ordinary allele. The combined allele profiles were defined as multilocus genotypes (electrophoretic types [ETs]). Based on allele frequencies for ETs, the genetic diversity for an enzyme locus (h) was calculated as follows: h = (1 -
xi2)[n/(n - 1)], where xi is the frequency of the ith allele and n is the number of ETs (43). The total mean genetic diversity (H) is the arithmetic mean of h values for all loci and represents the proportion of loci at which two randomly chosen genotypes can be expected to differ. To compute the H values, we used the program ETDIV, version 2.2 (58).
The genetic distance between each pair of different ETs was estimated by determining the mean character differences, and a similarity matrix was constructed by using the PAUP* program (53); the data were clustered by the unweighted pair group method with arithmetic averages (UPGMA) (46).
Genetic differentiation.
To estimate the relative genetic differentiation (Gst), we used Nei's equation, Gst = (Ht - Hs)/Ht (49), where Ht is the expected diversity in an equivalent randomly mating total population and Hs is the average diversity of the subpopulations. To compute the Gst values, we used the program ETDIV, version 2.2 (58). The indices range from 0, if there is no genetic differentiation at a given level, to 1, if there is maximal genetic differentiation (35). To test if the Gst values were significantly different from 0, we performed a chi-square test of independence as follows:
2 = nGst(a - 1), where n is the number of individuals and a is the total number of alleles. The degrees of freedom are (k - 1)(a - 1), where k is the number of subdivisions. The degrees of freedom and
2 values were summed across loci, and significance was examined at a P value of <0.05 (19, 59).
Linkage disequilibrium analyses.
To determine the extent to which populations exhibited nonrandom associations of alleles between loci, we used a multilocus index based on the distribution of allelic mismatches between pairs of isolates for all loci. The ratio of the variance in mismatches observed in a population (Vo) to the expected variance of the corresponding population at linkage equilibrium (random association of alleles) (Ve) provides a measure of linkage disequilibrium. If there is no linkage disequilibrium, Vo/Ve is 1. The significance of the difference between Vo and Ve was calculated by using a Monte Carlo procedure with 1,000 iterations, which was carried out with the LDV program (50).
Visualization and cluster analysis of plasmid profiles.
The plasmid contents of the isolates were visualized by using the Eckhardt procedure (12). Plasmid mobilities were determined in 0.7% agarose gels by using plasmids of R. etli bv. phaseoli CFN42 and Sinorhizobium meliloti 1021 as molecular size references. A plasmid profile similarity matrix was constructed by using the PAUP* program (53), and data were clustered by using the UPGMA algorithm (46).
PFGE of plasmids and pSym determination.
Pulsed-field gel electrophoresis (PFGE) was used to obtain accurate size estimates for plasmids from selected isolates and to identify the pSym plasmids. Intact genomic DNA was prepared in 0.8% agarose plugs and subjected to PFGE in a contour-clamped homogeneous electric field apparatus (CHEF-DRII; Bio-Rad) by following the manufacturer's instructions. Electrophoresis was carried out at 13°C with a constant voltage of 4.5 V cm-1 by using a two-block program as follows: block 1, 30-s initial switch, 60-s final switch, and 8-h run time; and block 2, 50-s initial switch, 100-s final switch, and 28-h run time. Saccharomyces cerevisiae chromosomes (Bio-Rad) were used as molecular size markers. Gels were stained with ethidium bromide, photographed, and transferred to nylon filters. Membrane hybridization with a nifH probe was performed as described below.
Analysis of nifH gene organization.
Genomic DNA from 34 isolates was digested with endonuclease BamHI, subjected to 1% agarose gel electrophoresis, stained with ethidium bromide, photographed, and transferred to nylon filters. An internal nifH fragment (
450 bp) from strain CFN42 was amplified with primers o1 and o3 (37), as described below. This fragment was labeled with digoxigenin-dUTP by using random primers, and detection was performed with anti-digoxigenin-alkaline phosphatase Fab fragments by using the chemiluminescence system and following the instructions of the manufacturer (Roche). Membrane hybridization and washing were performed under high-stringency conditions (65°C, 0.5x SSC [1x SSC is 0.15 M NaCl plus 0.015 sodium citrate]). Some membranes were also hybridized with a lambda probe to obtain accurate estimates of the hybridization signal sizes on gels normalized with HindIII-digested lambda DNA.
Plant nodulation tests and acetylene reduction assay.
Seeds of P. vulgaris cv. Negro Jamapa, Macroptilium atropurpureum, and Leucaena leucocephala cv. Peruvian were surface sterilized with sodium hypochlorite (54). Pregerminated seeds were placed in flasks filled with vermiculite, watered with an N-free plant nutrient solution (13), and inoculated with each of the 34 selected isolates. Plants were maintained in a growth chamber at 28°C with a photoperiod of 15 h. After 4 weeks, the numbers of nodules were counted, and nitrogen fixation was measured by the acetylene reduction assay (54).
PCR-restriction fragment length polymorphism (RFLP) analysis of the 16S rRNA gene.
The 16S rRNA genes of 34 isolates were PCR amplified by using primers fD1 and rD1 (55). The amplification products were digested with endonucleases Sau3AI, MspI, and PstI and visualized in agarose gels as described elsewhere (23). Type and reference strains of the five recognized species that nodulate common bean (R. etli bv. phaseoli CFN42, R. tropici CIAT899 and CFN299, R. leguminosarum bv. phaseoli USDA2671, R. gallicum bv. gallicum R602sp and FL27, and R. giardinii H152) were included for comparison.
PCR amplification and nucleotide sequencing.
Two isolates were selected for partial DNA sequencing of two chromosomal genes, rrs coding for 16S rRNA and glnII coding for glutamine synthetase, and two pSym genes, nifH and nodB, coding for the dinitrogenase reductase and N-acetylglucosamine deacetylase, respectively. For PCR amplification, a reaction mixture (50 µl) containing 1x PCR buffer (Gibco BRL), 1.5 mM MgCl2, each deoxynucleotide triphosphate at a concentration of 200 µM, each primer at a concentration of 0.2 µM, and 2 U of Taq polymerase was used. The following temperature profile was used for all amplifications: 3 min at 94°C, followed by 30 cycles of 1 min at 94°C, 1 min at the appropriate annealing temperature, and 1 min of extension at 72°C and a final extension cycle consisting of 72°C for 5 min. The annealing temperatures for the different primer pairs are indicated below. rrs gene fragments were amplified with primers PF2 (TACTGTCGATCTGGAGTATG) and PR1 (ATTGTAGCACGTGTGTAGCC) (annealing temperature, 60°C); glnII gene fragments were amplified with primers GSF1B and GSR2seq (56) (annealing temperature, 70°C); nifH sequences were amplified with primers o1 and o3 (37) (annealing temperature, 58°C); and almost complete nodB genes were amplified with primers nodB3F (56) and nodCRR (GAGACGGCGRCRRTGCTGGTTG) (annealing temperature, 65°C). The same primers were used for sequencing reactions with ABI Prism Big Dye chemistry, and the products were analyzed with an ABI377 sequencer (ABI, Foster City, Calif.).
Nucleotide sequence accession numbers.
The sequences obtained in this study have been deposited in the GenBank sequence database under accession numbers AF529011 to AF529022.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Distribution of ETs and plasmid profiles among 18 R. gallicum and 108 R. etli isolates from nodules of P. vulgaris plants over three consecutive yearsa
|
2 = 7.14; df = 4; P = 0.129), suggesting that there was temporally stable genotype dominance in the population. These dominant R. etli genotypes grouped in a tight cluster, as shown in Fig. 1 (mean genetic distance, <0.2), indicating that they form a clonal complex (14).
![]() View larger version (38K): [in a new window] |
FIG. 1. Dendrograms showing the genetic relatedness among chromosomal and plasmidic genotypes of R. etli and R. gallicum and the chromosome-plasmid profile combinations. (A) Genetic relatedness of the multilocus genotypes based on 10 isoenzymatic loci. The ET designations are indicated. (B) Genetic relatedness of the plasmid profiles based on the presence or absence of the different plasmid size classes. The plasmid profile designations are indicated. The number of isolates for each multiple ET or plasmid profile is given in parentheses. The genetic distance between each pair of ETs or plasmid profiles was estimated by determining mean character differences, and data were clustered by the UPGMA. The four most abundant ETs and their plasmid profile combinations are indicated by boldface type and solid lines; the other combinations are indicated by dashed lines.
|
|
View this table: [in a new window] |
TABLE 2. Genetic diversity, genetic differentiation, and linkage disequilibrium estimates for R. gallicum and R. etli populations associated with P. vulgaris plants in San Miguel Acuexcomac, Mexico
|
Linkage disequilibrium analyses.
A hierarchical linkage disequilibrium analysis was performed with the MLEE data set to estimate the extents of genetic exchange within and between the R. etli and R. gallicum populations. As shown in Table 2, both populations showed linkage equilibrium when only ETs were included in the analysis. When all isolates of each species were considered, the populations appeared to be in linkage disequilibrium, due to the presence of numerically dominant ETs (epidemic clones [32]). These results revealed an epidemic genetic structure for both species, in which the frequency of a few dominant genotypes increases to produce epidemic clones, but frequent genetic exchange occurs among the members of the population. Linkage disequilibrium was detected when the analysis was performed with either all 126 isolates or 43 ETs (Table 2), indicating that the extent of genetic exchange between R. etli and R. gallicum is negligible, if there is any exchange at all.
Diversity and cluster analysis of plasmid profiles.
Eleven plasmid size classes were identified among the 18 R. gallicum isolates; the sizes ranged from
50 to >1,500 kb and represented seven distinct profiles (Table 1 and Fig. 1). R. gallicum isolates harbored two to four plasmids (average, 3.1 plasmids), and all isolates contained a >1,500-kb megaplasmid, which was absent in the R. etli population (Table 3). Profile P4 was the dominant profile; it was observed in 44% of the isolates (Table 1 and Fig. 1) and was indistinguishable from the profile of the type strain of R. gallicum bv. gallicum R602sp (Table 3). Fifteen plasmid size classes were detected among the population of 108 R. etli isolates, and the plasmid sizes ranged from
125 to
700 kb. Twenty-six plasmid profiles were found, consisting of two to five plasmids (average, 3.5 plasmids). All of the R. etli profiles included the
700-kb plasmid, which was not found in the R. gallicum population (Table 3). The R. etli population displayed a high degree of plasmid profile dominance, with many isolates sharing a few plasmid profiles. Profiles P1, P2, and P3 were found in 57% of the R. etli isolates (Table 1 and Fig. 1); these profiles were the predominant profiles and were recovered in all 3 years. A chi-square test indicated that the frequencies of these three profiles varied significantly over the 3 years (
2= 19.32; df = 4; P < 0.001). Profile P1 was the most abundant profile in 1996, P2 was the most abundant profile in 1994, and P3 was the most abundant profile in 1995. These apparent temporal changes in plasmid profile dominance in the R. etli population could have been due simply to the small sample size or could reflect actual changes in the plasmid combinations selected in the crop seasons (for example, due to fluctuations in environmental conditions).
|
View this table: [in a new window] |
TABLE 3. Multilocus genotypes (ETs), plasmid profiles, 16S ribosomal DNA RFLP patterns, nifH hybridization patterns, and estimated sizes of plasmids and nifH hybridization bands for 34 selected P. vulgaris isolates and reference strains
|
Based on the MLEE and plasmid profile data, 34 isolates were selected for further molecular analyses; this selection included 10 isolates from the R. gallicum population and 24 isolates from the R. etli population. This selection included isolates from dominant ETs and plasmid profiles, as well as isolates displaying unique ETs or plasmid profiles (Table 3).
PFGE and pSym identification.
PFGE was used to obtain accurate plasmid size estimates for the 34 isolates selected for molecular analyses, and the values are shown in Table 3. These values were used to estimate the plasmid sizes for the whole sample of 126 isolates, based on Eckhardt plasmid profiles.
The pSym plasmids of 20 isolates (4 R. gallicum and 16 R. etli isolates) were identified by Southern analysis of plasmid profiles resolved by PFGE by using a nifH probe and are indicated in Table 3. Most R. gallicum isolates harbored 550-kb pSym plasmids; the single exception was isolate IE2751, which harbored a 525-kb pSym. The sizes of the pSym plasmids of R. etli isolates ranged from 350 to 550 kb. One-half of the R. etli isolates contained a 390-kb pSym, like the R. etli type strain CFN42. In general, pSym plasmids of the same size were associated with identical or similar plasmid profiles; for example, isolates with profiles P2 and P3 harbored 390-kb pSym plasmids. The only exception was isolate IE4874, which contained a 450-kb pSym, whereas the other two isolates with profile P1 contained 550-kb pSym plasmids (Table 3).
Determination of nifH copy number by Southern hybridization.
For the 34 isolates analyzed, six different hybridization patterns were recorded (Table 3). All R. gallicum isolates exhibited a single hybridization signal, and most of them shared an
8.3-kb hybridizing band (pattern C1) with R. gallicum bv. gallicum strains R602sp and FL27 (Table 3), supporting classification of our isolates as members of R. gallicum bv. gallicum. All R. etli isolates contained multiple (two or three) nifH copies, and two groups of related patterns were observed. Patterns A1 and A2 had two bands with sizes of
12.4 and
16.6 kb, while patterns B1 and B2 had two bands with sizes of
5.6 and
9.8 kb; the latter pattern was identical to that displayed by type strain R. etli CFN42 (Table 3).
Within the R. gallicum and R. etli populations, nifH hybridization patterns could be found scattered throughout the different branches of the MLEE and plasmid profile dendrograms, further supporting the conclusion derived from the comparative cluster analyses of MLEE and plasmid profile data that lateral transfer of plasmids occurs within each species but not between the species. However, in some instances an association between plasmid profile and nifH pattern within species was observed. For example, all the R. etli isolates that had a P1 plasmid profile displayed a type A1 nifH hybridization pattern, as summarized in Table 3. Notable associations between pSym size and nifH hybridization pattern were also found. For example, R. etli isolates with the A1 nifH pattern harbored 450- to 550-kb pSym plasmids, the A2 nifH patterns were associated with 390- to 425-kb pSym plasmids, and the B1 and B2 nifH patterns were strictly associated with 390-kb pSym plasmids (Table 3). R. gallicum isolates with 550-kb pSym plasmids yielded only the C1 nifH pattern, while a 525-kb pSym was associated with the distinct C3 nifH pattern. A 550-kb pSym was detected in the R. etli and R. gallicum populations, but the nifH hybridization patterns revealed that these plasmids were distinct.
Plant nodulation tests and symbiotic effectiveness.
All 34 selected isolates (Table 3) nodulated and fixed nitrogen with P. vulgaris plants (Table 4). However, R. gallicum isolates induced only one-third of the nodules induced by R. etli isolates. All R. gallicum isolates formed effective nodules on M. atropurpureum, whereas R. etli isolates induced less than one-third of the nodules induced by R. gallicum isolates on this host and were ineffective (Table 4). All R. gallicum isolates nodulated and fixed nitrogen on L. leucocephala plants, further supporting classification of these isolates as R. gallicum bv. gallicum, while all R. etli isolates were unable to nodulate this host (Table 4).
|
View this table: [in a new window] |
TABLE 4. Average numbers of nodules formed by 10 R. gallicum and 24 R. etli isolates on three leguminous hostsa
|
1.5 kb) for 34 isolates were PCR amplified and digested with endonucleases Sau3AI, MspI, and PstI. Four different composite restriction patterns were found (Table 3). R. gallicum isolates displayed three different patterns (AAA, AAB, and ABA); pattern AAB was identical to the pattern of R. gallicum FL27, while pattern AAA was identical the pattern of R. gallicum R602sp. All R. etli isolates displayed the same restriction pattern (BAB), which was identical to the pattern of R. etli CFN42 (Table 3). The restriction sites for each enzyme were mapped on the rrs sequence of R. gallicum R602sp (accession number U86343) and are indicated in Table 3. The patterns displayed by other Rhizobium reference species were different from those of our isolates (data not shown).
Nucleotide sequence analyses.
Partial sequences of two chromosomally encoded genes (rrs and glnII) and two pSym-encoded genes (nifH and nodB) were determined for R. gallicum isolate IE988 (ET12 in Table 1 and Fig. 1) and R. etli isolate IE2730 (ET1 in Table 1 and Fig. 1) and, when the sequences were not available in public sequence databases, also for R. gallicum strains R602sp and FL27.
The 558-bp rrs sequence segment of isolate IE988 was very similar to that of strain R602sp (Table 5). However, the level of similarity with the R602sp and FL27 sequences reported by Sessitsch et al. (44) was lower, due to the presence of several insertions which are absent in the R602sp sequence reported by Amarger et. al. (3) and in our IE988 sequence. The rrs sequence of isolate IE2730 differed only at two nucleotides from the rrs sequence of R. etli CFN42 (Table 5).
|
View this table: [in a new window] |
TABLE 5. Numbers of pairwise nucleotide differences and percentages of similarity for aligned rrs, glnII, nifH, and nodB gene segments of the R. gallicum IE988 and R. etli IE2730 isolates compared with reference strains
|
The 389-bp nifH sequence segment obtained for isolate IE988 was most similar to the nifH sequence segments of strains R602sp and FL27, while the nifH sequence of isolate IE2730 was identical to the sequence of the three paralogous nifH copies of R. etli CFN42 (Table 5).
An almost complete nodB sequence segment (571 bp) was obtained for isolates IE988 and IE2730 and for strains R602sp and FL27. The nodB sequence segment of IE988 was very similar to the nodB sequence segments of strains R602sp and FL27, while the nodB sequence of IE2730 was identical to the nodB sequence of R. etli CFN42 (Table 5).
Taken together, these results strongly support the taxonomic placement of isolate IE988 in R. gallicum bv. gallicum and the taxonomic placement of isolate IE2730 in R. etli bv. phaseoli.
|
|
|---|
The three dominant R. etli ETs identified in the present work correspond to those found in the previous study (45), indicating that the genetic composition of the populations at this site is stable both in terms of space (six plots sampled in a single year [45]) and in terms of time (one plot sampled over 3 years [this study]). Furthermore, the two distinct lineages that were consistently recovered over the 3 years sampled correspond to the previously described genetic divisions I and III (45). Linkage disequilibrium analyses confirmed that recombination takes place within but not between these lineages. Importantly, the set of molecular markers analyzed in this study, along with the host range experiments, unambiguously placed our division I population in R. gallicum bv. gallicum and the division III population in R. etli bv. phaseoli. Therefore, the two lineages correspond to two distinct species and are not compartments within an R. etli population, as suggested previously (45). As a consequence, we concluded that the genetic structure of the R. etli and R. gallicum populations at this site is epidemic but not reticulated, as defined by Maynard Smith and colleagues (32).
Strong evidence that there is genetic exchange within but not between the R. etli and R. gallicum populations was also provided by the analysis of the plasmid compartment. Both the plasmid profiles and the nifH hybridization patterns provided evidence that there is plasmid transfer among distinct multilocus genotypes within each species (Fig. 1 and Table 3). The panmictic plasmid transfer pattern (48) found for both species could have been the result of the selective pressure imposed by agricultural practices, as suggested by Wernegreen and colleagues (57), or may have been an intrinsic feature of these Rhizobium species. The evidence that recombination took place also in the chromosomal compartment, as demonstrated by the linkage disequilibrium analyses, favors the latter possibility. Taken together, these results indicate that genetic exchange is an important source of variation and cohesion in the ecology and evolution of the two species, although it is not great enough to prevent the emergence of epidemic clones that are recovered from nodules at a high frequency. This finding may be interpreted as evidence that there is strong selective pressure imposed by the host, which favors the maintenance of particular chromosome-plasmid associations. It will be interesting to test whether the epidemic clones are particularly competitive for nodulation or if they are simply numerically dominant clones in the soil. It is important that both species nodulate the same individual plants, indicating that mechanisms other than host, spatial, or temporal isolation may account for the genetic barrier between them.
The 34 strains selected for molecular and host range analyses revealed that the subsample of 10 R. gallicum isolates constitutes a diverse lineage that displays three of the four 16S PCR-RFLP patterns detected, three different nifH hybridization patterns, and two pSym sizes (Table 3). In addition to P. vulgaris, these isolates effectively nodulated the other two hosts tested (Table 4). On the other hand, the 24 R. etli isolates effectively nodulated only P. vulgaris. All of the isolates displayed the same 16S PCR-RFLPs, but they harbored six pSym size classes and displayed four nifH hybridization patterns. The comparison of plasmid profiles, nifH hybridization patterns, and pSym sizes revealed the complex and dynamic structure of the plasmidic compartment within species. Generally, pSym plasmids that were the same size and had identical nifH hybridization patterns were associated with identical plasmid profiles. However, there were some exceptions that might indicate the existence of genetic rearrangements which affect the plasmidic compartment. These rearrangements can operate within particular strains, as extensively documented for CFN42 and NGR234 (7, 15, 38, 39), and may be coupled with lateral transfer within populations, as proposed in the present study.
Only 14% of all of the isolates examined in this study are R. gallicum isolates. This is in good agreement with the results of a previous study (45), in which it was reported that 10% of the isolates obtained from P. vulgaris corresponded to genetic division I (R. gallicum). In that study, however, R. gallicum isolates were recovered from P. coccineus nodules at nearly the same frequency (54%) as R. etli isolates. These results indicate that R. gallicum and R. etli nodulate both hosts, although the latter species is clearly more competitive for P. vulgaris nodulation. Thus, although R. gallicum was originally isolated from P. vulgaris nodules in France (3), this plant may not be its primary host. This possibility is further supported by the report that R. gallicum strains were isolated in Canada from Onobrychis viciifolia and Oxytropis riparia (tribe Galegae), suggesting that not all potential hosts for this species have been identified yet (26).
Besides its broad host range, R. gallicum has a wide geographic distribution, which raises a question about its biogeography. R. gallicum populations have been isolated from bean nodules in several European countries and Tunisia (3, 21, 33), and only one American R. gallicum bean isolate (strain FL27) has been reported previously (44). The present study provides the first report of an R. gallicum population nodulating beans in America, and it is the first population genetic analysis for the species that was performed. Sessitsch and colleagues (44) reported that the European R. gallicum strains can be distinguished from the Mexican FL27 strain by an RFLP resulting from a single nucleotide difference in the 16S rRNA gene, which created a PstI restriction site that was absent in all the French and Austrian isolates which they analyzed but was present in the Mexican strain FL27. However, this restriction site was absent in 60% of our R. gallicum isolates (Table 3) and hence cannot be used to distinguish European R. gallicum strains from American R. gallicum strains. Furthermore, the plasmid profiles and nucleotide sequences of chromosomal and plasmidic genes of the Mexican population were very similar to those of R. gallicum type strain R602sp (in some instances even more similar than they were to those of Mexican strain FL27). These observations suggest that if R. gallicum was imported from America to Europe, as previously suggested (44), the lineages that gave rise to the European isolates remain in America and that the divergence between intercontinental populations is less pronounced than the divergence suggested by Sessitsch et al. (44). The presence of R. etli and R. gallicum in Tunisian soils recently cropped with beans and the absence of these organisms in fields which had not been cultivated with beans (33) further support the hypothesis that P. vulgaris microsymbionts were introduced by being carried along with bean seeds (36). Interestingly, R. gallicum bv. phaseoli was not present in the isolates from San Miguel Acuexcomac, even though R. etli bv. phaseoli and R. gallicum bv. gallicum may have coexisted at this site for centuries. One possible explanation for this is that R. gallicum bv. phaseoli may be outcompeted by R. etli bv. phaseoli and R. gallicum bv. gallicum strains. An alternative explanation is that since R. leguminosarum has not been found in Mexico, the hypothesized bridge for the conjugal transfer of pSym from R. etli bv. phaseoli to R. gallicum is missing (3). Perhaps the introduction of common beans along with their seed-borne symbionts to other continents exerted selective pressure and provided the ecological opportunity for lateral transfer of symbiotic information to resident Rhizobium populations, as in the case of the symbiotic island acquired by saprophytic Mesorhizobium spp. native soil populations in New Zealand (52).
In conclusion, our data are consistent with the biogeographical hypothesis of an American origin of the R. gallicum lineage (44), a view that is also consistent with the report of Canadian R. gallicum isolates obtained from O. viciifolia and O. riparia (26). Alternatively, however, R. gallicum could have a wide geographical distribution and a long evolutionary history of adaptation to different environments and leguminous hosts. Further research is clearly needed to examine these alternative hypotheses. We are currently analyzing the sequences of several chromosomal and plasmidic genes of R. gallicum isolates from different continents and hosts in order to elucidate the phylogeographic origin and dispersal pathways of this species.
We thank the community of San Miguel, especially the Aguilar family, for their kindness and help with fieldwork. B. Gaut is gratefully acknowledged for providing sequencing facilities at the University of California Irvine, and René Hernández is acknowledged for sequencing at IBT, UNAM. We thank Marco A. Rogel, Aldo Valera, Laura Espinosa, and Rafael Díaz for providing technical assistance and P. Burgos for designing primers PF2 and PR1. We acknowledge David Romero and Susana Brom for providing critical comments on the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»