Previous Article | Next Article ![]()
Applied and Environmental Microbiology, February 2003, p. 1067-1074, Vol. 69, No. 2
0099-2240/03/$08.00+0 DOI: 10.1128/AEM.69.2.1067-1074.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Research Institute of Agricultural Microbiology, St. Petersburg, Pushkin 196608,1 Department of Forage Sciences, N. I. Vavilov Russia Research Institute of Plant Genetic Resources, St. Petersburg 190000, Russia,3 Department of Applied Chemistry and Microbiology, Biocenter 1, Fin-0014 University of Helsinki, Finland,2 Department of Biology, University of York, York YO10 5YW, United Kingdom4
Received 10 June 2002/ Accepted 13 November 2002
|
|
|---|
|
|
|---|
The interesting symbiotic properties, together with the extensive taxonomic knowledge about R. galegae, prompted us to use the species to answer important questions related to the diversity and evolution of rhizobia and their symbiosis with legumes. A recent study indicated that R.galegae bv. orientalis strains have lower genetic diversity than R. galegae bv. officinalis strains (37). However, very few biovar orientalis strains were included in that study. In order to obtain a larger sample, R. galegae strains belonging to both biovars were collected during an expedition to the Caucasus, one of the centers of diversity of G. orientalis but not of G. officinalis. The collection was used for two different purposes.
First, we wanted to investigate the influence of the host plant on the genetic diversity of rhizobial populations. The host plant is an important factor shaping the genetic structure of a natural population of rhizobia (1, 2, 4, 9, 42, 45). The genomic diversity of the collected rhizobia was assessed by amplified fragment length polymorphism (AFLP) fingerprinting and by restriction fragment typing of the rRNA gene internal transcribed spacer (ITS). It is generally believed that insertion sequence (IS) elements in rhizobia promote genetic diversification through genomic rearrangements and recombination (8, 11, 17). The use of IS elements as probes for Southern hybridization can provide high-resolution fingerprints of rhizobial strains (26, 27, 31) and some information for understanding the function of these elements in rhizobial populations (2). An IS-related sequence, closely related to an Agrobacterium transposase from the widely studied plasmid pTiA6NC, has been found in the symbiotic region of R. galegae bv. orientalis strain HAMBI540 (24, 33), so we also used this IS sequence for fingerprinting the R. galegae strains.
Second, the collection was used to get a picture of the roles that the different parts of the symbiotic genome play in processes related to the evolutionary adaptation of rhizobia to different host plants. The genetic diversity in the symbiotic gene regions of R. galegae strains belonging to the two biovars was studied by PCR-restriction fragment length polymorphism (RFLP) analysis of five different symbiotic loci containing previously identified nod box sequences (32).
|
|
|---|
1 ha each) along a 300-km route through a mountainous part of the Krasnodar region in the northwest of the Caucasus (Fig. 1). An interesting observation was that the Galega species were found growing separately, with the exception of only two sites (Fig. 1, sites 104 and 126). Galega species in the area surveyed were readily distinguished by morphological traits of inflorescences and of the root system. G. orientalis was always part of well-established plant associations, whereas G. officinalis preferred disturbed or destroyed ecotopes and was part of the pioneer flora. The G. officinalis population was homogeneous morphologically, while the G. orientalis population was differentiated into distinct ecotypes, as is typical for the species. At least four well-differentiated ecotypes of G. orientalis were present in the region: mountain ecotypes A, from the Lagonaki tract (site 161), and B, from the mountain Otrub (site 97), were characterized by great height, good foliation, and a large inflorescence (up to 40 to 50 cm). The mountain ecotype C from Mezmay (site 160) was distinguished by dark foliage, and ecotype D from Khadyzhensk (site 176) and from Kamenoostrovsky (sites 104, 124, 126, and 127) was not as tall as the mountain ecotypes; had paler, shorter inflorescences; and sometime grew together with G. officinalis plants (sites 104 and 126).
![]() View larger version (37K): [in a new window] |
FIG. 1. Sampling locations. Solid circles, sites where G. orientalis grew; open circles, sites where G. officinalis grew; split circles, sites where both host plants grew.
|
|
View this table: [in a new window] |
TABLE 1. R. galegae strains used in the study
|
DNA isolation.
Total genomic DNAs of all R. galegae strains were obtained after lysozyme-sodium dodecly sulfate lysis followed by phenol-chloroform extraction and ethanol precipitation (12).
AFLP analyses.
The AFLP procedure was performed as described by Vos et al. (39), with some modifications. One hundred to 500 ng of DNA was digested with EcoRI and TruI (MseI) (MBI Fermentas) restriction enzymes in a PTC-200 thermal cycler (MJ Research, Watertown, Mass.) and, in the same step, ligated with double-stranded adapters specific for each restriction half-site. (7). Two sets of primers were used in separate PCRs: one set with two selective nucleotides (boldface) (MseI-gc [GAT GAG TCC TGA GTA AGC] and EcoRI-gc [GAC GTC GTA CCA ATT CGA GC]) and another set with a mixture of two and three selective nucleotides (MseI-gc and EcoRI-gag [GAC GTC GTA CCA ATT CGA GAG]). The PCR conditions were as described previously (37). Silver-stained polyacrylamide gels, prepared as described by Dresler-Nurmi et al. (7), were scanned with an AGFA scanner and analyzed using Bionumerics software, version 2.0 (Applied Maths, Kortrijk, Belgium). An unweighted pair group method with averaging (UPGMA) dendrogram was generated from all strains by using Pearson correlation coefficients. For further analyses, including PCR-RFLP and IS fingerprinting, a representative strain (Table 1) was chosen from each AFLP group. Based on AFLP data, a UPGMA dendrogram was constructed for each biovar. We then chose representatives from each AFLP cluster and loaded them on the same gel; there were 20 R. galegae bv. orientalis and 21 R. galegae bv. officinalis strains with distinct fingerprints. Two different sets of primers were used and two gels, containing all 41 chosen strains, were obtained. The final dendrogram was constructed from the combined gel. The HAMBI1141 and HAMBI540 strains were excluded from the final dendrogram construction to avoid distorting the comparison of the branching depths of the two biovars studied.
PCR amplification.
The PCRs were carried out with 5 to 10 ng of total DNA. Dynazyme DNA polymerase (Finnzymes) was used with the supplied buffer (1.5 mm MgCl2) and deoxynucleotide triphosphates. The primers fD1 and rD1 were used to amplify 16S rRNA genes (41). For ITS amplification, the primers FGPS1490-72 and FGPL132' (22) were used. The locations of primers for amplifying regions containing nod boxes are shown in Fig. 2. Primers for nod box regions were constructed using Jellyfish software (Biowire): nb1f, TCC ATC TCC TAG ATG CCT CA; nb1r, GGC TAC GCA GTG ATT TGA TG; nodDf, CAG ATG TTG CCT GGG TCG; nodDr, ATC AAC CTC AGT CAA CCG GC; nb2f, TCG ATT GTT TGA ATG CCA TA; nb2r, GGC AAC ACG TCC AAC TTT CT; nb3f, GCG TGA ATT GGT GTC ATC CT; nb3r, AAG CTC GTC GCC AAA ATA GA; nb4f, GCC GCG ATA TAC TGA CCA TT; nb4r, TAT TCG AGG CGT GGA TGT CT; nb5f, GGT CAT GCT CTC GAT GGT TT; and nb5r, CTT CAA ATC CAC GCA ATC AA. All PCRs were performed using a PTC-200 Peltier Thermal Cycler (MJ Research) with the following standard temperature profile: an initial denaturation at 95°C for 3 min; 30 cycles of denaturation (1 min at 94°C), annealing (1 min at 55°C), and extension (1.5 min at 72°C); and final extension at 72°C for 3 min. When amplification was weak, the PCR was repeated with a lower annealing temperature (52 or 53°C). The amplified DNA was examined by electrophoresis in 1% agarose gels.
|
View larger version (8K): [in a new window] |
FIG. 2. Maps of symbiotic gene regions of Rhizobium galegae. The locations of primers used in the study and of amplified fragments are shown.
|
xi 2)/(n - 1), where xi is the frequency of the ith RFLP type and n is the number of RFLP types (5). Similarities between the HAMBI540 and HAMBI1141 strains in symbiotic loci were estimated from the proportion of shared and unshared restriction fragments by using a similarity coefficient (20).
Southern hybridization.
An aliquot (0.5 to 1 µg) of the total DNA sample was digested to completion with HindIII (MBI Fermentas). Restriction fragments were separated by electrophoresis in 1% agarose. DNA of phage
digested with HindIII was included on the gel as a size marker. Southern blotting of DNA, hybridization, and detection were done according to instructions provided with the nonradioactive labeling and detection kit (Roche). The digoxigenin (DIG) labeling of the IS probe was performed using PCR amplification with the ISf (CGA TTG CTA CAA TGG CTT CA) and ISr (TCT TTC TTC CAC ATG CAA CC) primers, flanking an inner fragment of the unnamed IS element, constructed from sequence data of R.galegae strain HAMBI540 (34).
|
|
|---|
One hundred and one R. galegae strains were isolated from nodules of wild-growing G. officinalis and G. orientalis in this work. All the strains had 16S rRNA MspI restriction profile identical to that of the reference R. galegae strains HAMBI540 and HAMBI1141 (data not shown). AFLP fingerprinting (Fig. 3), ITS PCR-RFLP, and nod box PCR-RFLP (Table 1) clearly distinguished the strains isolated from each host plant. However, five strains originally isolated from G. orientalis nodules belong by their AFLP and ITS and nod box PCR-RFLP patterns to the G. officinalis cluster (Table 1). The symbiotic performances of these five isolates were studied. They produced relatively large pink nodules on G. officinalis plants and small white or light-green nodules on G. orientalis plants. Thus, we can conclude that these strains are R. galegae bv. officinalis even though they were isolated from nodules of the nontypical host plant G. orientalis.
![]() View larger version (30K): [in a new window] |
FIG. 3. UPGMA dendrogram constructed from AFLP data. The relative positions of R. galegae bv. orientalis HAMBI540 and R. galegae bv. officinalis HAMBI1141 are indicated by dashed lines.
|
The data obtained in the study by molecular characterization of rhizobial isolates correlates well with host plant diversity. AFLP analysis, which is based on whole-genome variability, shows that the Caucasian populations of R. galegae bv. orientalis are more diverse than the R.galegae bv. officinalis populations. The branching depth in the UPGMA dendrogram is deeper in the case of R. galegae bv. orientalis (Fig. 3). Interestingly, the R. galegae bv. orientalis strains group into two different AFLP clusters, A and B, the latter being closer to the R. galegae bv. officinalis cluster than to the main R. galegae bv. orientalis AFLP cluster (Fig. 3). Another interesting finding is that HAMBI1141 lies outside the main AFLP cluster of R. galegae bv. officinalis. This strain also has a unique ITS RFLP type, X (Table 1), and we believe that it is a rare genotype not typical of this region.
The results of the PCR-RFLP analysis, which allows the study of variability at particular loci (Table 1), are in agreement with those of the AFLP analyses. The corresponding values of heterogeneity in the ITS region for the two biovars were 0.87 for R. galegae bv. orientalis and 0.56 for R. galegae bv. officinalis. The average heterogeneity in the symbiotic regions analyzed was 0.42 for R. galegae bv. orientalis and 0.13 for R. galegae bv. officinalis. The only symbiotic region in which R. galegae bv. officinalis strains are more diverse is nod box 5. Overall, the strains of R. galegae bv. orientalis are more diverse than the R. galegae bv. officinalis strains in symbiotic regions as well as in the ITS region. This finding, together with the AFLP results, provides strong evidence that within the Caucasus region studied R. galegae bv. orientalis strains are more diverse than R. galegae bv. officinalis strains.
Interestingly, the results appear to contradict a previous study (37) in which G. officinalis strains were shown to be more diverse than G. orientalis strains. However, we believe that there is no contradiction: the population of G. officinalis studied by Terefework et al. (37) was geographically heterogeneous, representing five different locations (New Zealand, United Kingdom, Bulgaria, Italy, and Argentina). At least two of them (New Zealand and Bulgaria) are places where native G. officinalis is widely represented. At the same time, the population of G. orientalis analyzed by Terefework et al. represented two locations, Russia and Finland (where G. orientalis for the most part is an introduced plant) and included few strains isolated in the Caucasus. Therefore, the observations of Terefework et al. (37) are in concordance with the results reported here, emphasizing the importance of host plant diversity for the diversity of rhizobia.
Analysis of the symbiotic region is of prime interest when subtle genetic adaptation of rhizobia to particular host plants is under investigation. Five different fragments containing nod box sequences were studied in this work. The presence of nod boxes itself points to the great importance of these regions in symbiotic interaction. We constructed five pairs of primers (see Materials and Methods) to amplify the region adjacent to each of five nod boxes, which were identified by Suominen et al. (32) in the symbiotic regions of the R. galegae bv. orientalis reference strain HAMBI540 (Fig. 2). Regions nb3-nodD2, nb4-nodU, and nb5-nodE were amplified successfully for almost all R. galegae bv. officinalis and R. galegae bv. orientalis strains (Table 1). However, the nb2-dctA region was successfully amplified only for R. galegae bv. orientalis and not for R. galegae bv. officinalis strains (Table 1). The same was true for the nb1-nodB region (not shown), but when the nb1f primer was replaced by the nodDf primer, located inside the nodD1 gene (Fig. 2), the corresponding fragment comprising the nodD1-nodB region (which includes nb1) was amplified successfully for most of the strains of both biovars (Table 1). The last finding means that there is some difference between the nucleotide sequences of nod box 1 in strains belonged to different biovars.
The most conserved symbiotic genes seem to be in the nod box 3 and nod box 4 regions. The nod box 1 region is highly conserved among R. galegae bv. officinalis strains, and the nod box 5 region is highly conserved among R. galegae bv. orientalis strains. We expected to find high levels of divergence in the nod box regions because of the high proportion of noncoding sequence. However, they proved to be highly conserved despite the high noncoding content, e.g., 56% of noncoding sequence in the case of the most conserved nod box 3 region (unpublished sequence data for strain HAMBI540). We assume that there may be strict functional constraints on the noncoding regions adjacent to nod boxes (6).
The Nei similarity coefficient was calculated by comparing the restriction patterns of strains HAMBI540 and HAMBI1141 to give a measure of differences between different biovars in particular symbiotic loci. The coefficients were 0.77 for the nod box 3 region, 0.76 for nod box 4, 0.52 for nod box 5, and 0.47 for nod box 1. Thus, the greatest differences were found in the nod box 1 and nod box 5 regions, and the smallest differences were found in the nod box 3 and nod box 4 regions. The nod box 2 region has not been taken into account, as there was no amplified fragment for the R. galegae bv. officinalis strains. Thus, we can suggest that the nod box 1 and 5 regions with adjacent genes might play a more important role in host plant adaptation and symbiotic tuning. This suggestion is quite plausible because we know that the nodA and nodEF genes (located adjacent to nod boxes 1 and 5) are involved in synthesis of unsaturated fatty acids and their attachment to the Nod factor backbone (23). It is established that the unsaturated fatty acid strategy is the main factor for Rhizobium-Galega recognition (reviewed in reference 36). The presence and importance of different allelic forms of these host range genes (46) in different biovars of R. galegae cannot be excluded. However it seems there is no clear correlation between the genetic background and the Nod factors produced. It was shown that R. galegae strains belonging to different biovars produce nearly identical Nod factors, which are a mixture including a series of C18 and C20 fatty acids with carbonyl-conjugated double bonds (44).
To obtain the whole picture of genome adaptation to particular host plants, this comparative study should be expanded to include more symbiotic genes. Thus, further investigation is needed to complete the picture of symbiotic adaptation in the Rhizobium-Galega symbiosis. This work demonstrates an approach to trace evolutionary pressure along the symbiotic genomes of two closely related rhizobium biovars. All the analyses in this study, AFLP and RFLP of ITS and symbiotic regions, delineated the two biovars. The absence of shared restriction patterns or combined profiles from the ITS and symbiotic regions strongly suggests that there has been no recombination between the biovars. Thus, the major factor contributing to the divergence of this part of the symbiotic machinery could be the host plant.
We used a pair of primers which allowed us to amplify the inner fragment of an IS element closely related to an Agrobacterium transposase, which was revealed in the symbiotic region of R. galegae bv. orientalis strain HAMBI540 (33). Only 16 strains of R. galegae bv. orientalis were found to contain copies of the IS sequence (Table 1). Furthermore, we noticed that all IS-positive strains, including the reference strain HAMBI540, had closely similar rRNA ITS types: patterns O and G (Table 1) differ only in their HaeIII restriction patterns. Strains of these genotypes formed a clearly separated cluster (B) in the AFLP dendrogram (Fig. 3). It is particularly remarkable that 13 of these strains were isolated from the same site (Table 1). By using a DIG-labeled IS probe for Southern hybridization, it was shown that the strains contain from three to nine copies of the IS element; HAMBI540 contains two copies (Fig. 4). Each hybridizing band corresponds to at least one copy of the IS element, because the IS sequence does not contain HindIII restriction sites. Strains with identical IS hybridization patterns had identical AFLP fingerprints.
![]() View larger version (49K): [in a new window] |
FIG. 4. IS fingerprinting by Southern hybridization of R. galegae strains containing the IS element.
|
Our work shows that the study of rhizobial diversity in gene centers of the host is a way to answer some interesting questions in plant-microbe interaction and can shed some light on patterns and constraints of evolution in rhizobia.
|
|
|---|
,ß-unsaturated N-acyl substitutions. Mol. Microbiol. 34:227-237.[CrossRef][Medline]
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»