Previous Article | Next Article ![]()
Journal of Bacteriology, August 2006, p. 5570-5577, Vol. 188, No. 15
0021-9193/06/$08.00+0 doi:10.1128/JB.00335-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Soybean Genomics and Improvement Laboratory, Agricultural Research Service, U.S. Department of Agriculture, Beltsville, Maryland 20705,1 Penn State Berks College, Berks Campus, P.O. Box 7009, Reading, Pennsylvania 196102
Received 7 March 2006/ Accepted 25 May 2006
|
|
|---|
|
|
|---|
Over many decades there has been a concerted effort to collect and evaluate medics from various locations worldwide. Subsequent to their acquisition, plant introductions were usually evaluated for their potential application in agriculture. These programs were and are relatively easy since a plant's morphology, physiology, and performance are readily measured. By comparison, the collection of rhizobial microsymbionts, capable of nodulating these medics, has been more arbitrary since there is no efficient method to discriminate between genetically related clusters of these bacteria. Traditionally, a limited number of rhizobial cultures were isolated and then were tested on several varieties of a crop and those with the best overall performance were chosen for manufacture of inoculants. In the case of Medicago, USDA 1002 (strain 3Doa2 or ATCC 9930) was isolated in 1919 by USDA scientists in Virginia from a soil on the Arlington Farm (now the location of the Pentagon). This strain was subsequently chosen as the type for the species (then Rhizobium meliloti) based on its superior symbiotic performance.
With the development of techniques in molecular biology, it has become possible to examine newly isolated cultures for a much broader array of characters. These techniques have been used to investigate genetic diversity among strains, and the data have also been used to provide support for decisions in rhizobial taxonomic classification. One of the more popular techniques is DNA fingerprinting by PCR. This method depends on the use of primers for repeat sequences that reside throughout bacterial genomes (30). One such analysis was developed with rhizobia that form symbioses with Medicago (7). Other PCR-related techniques that have been applied to obtain DNA fingerprints of rhizobial genomes are random amplified polymorphic DNA (32) and amplified fragment length polymorphism (31). These tools have been used to examine genetic diversity among Medicago-nodulating rhizobia isolated in Croatia (3), Italy (5, 19, 20), Latin America (8), Tunisia (1, 33), and the Caucasian and Central Asian regions (2, 22). Even though with these methods the existence of genetic diversity among Medicago-nodulating rhizobia has been revealed, without sequence data it is not possible to obtain conclusive estimates of levels of their divergence. The data used to estimate divergence are inconclusive because PCR products of equal molecular size across lanes need not necessarily be identical. Besides the qualitative nature of fingerprint data, another limitation of the approach is that results are not portable and consequently numerous standard reference strains need to be included in every analysis.
DNA fingerprinting is one of a variety of methods, collectively referred to as a consensus or polyphasic approach, that are used in the characterization of rhizobia. However, there is lack of agreement among bacterial taxonomists and population geneticists on which of these methods are the most appropriate to define bacterial species limits (9, 25). As a consequence, rhizobial taxonomic decisions mostly have been subjective and arbitrary.
In two recent reviews, a theoretical basis for the use of a multilocus genotypic approach to define bacterial species limits was discussed (9, 25). The advantage of such an approach is that evidence can be provided to infer both founding genotypes and likely patterns of descent within clusters of closely related strains. The first attempt to apply a multilocus approach for genotyping Medicago-nodulating rhizobia was reported by Eardly et al. (10), who used multilocus enzyme electrophoresis (MLEE) to characterize 231 strains that originated predominantly from southwest Asia. Because results generated by MLEE also are not portable, comparisons of results with large numbers of strains in separate data sets are complex and require that multiple reference standards be included in each analysis. A significant improvement over MLEE is the multilocus sequence typing (MLST) method, where alleles of at least seven loci are identified by direct nucleotide sequence comparisons (16) rather than electrophoretic migration patterns of enzyme electromorphs. The combinations of alleles for the loci examined in each strain are then used to derive an allelic profile or sequence type (ST). The ST data are maintained in a database that can be manipulated to permit subsequent comparisons of large numbers of bacterial genotypes without the need to repeatedly analyze reference cultures. MLST results for numerous pathogenic bacterial species have already been reported (11), but to our knowledge this strategy has never been applied to rhizobia.
MLST also can be applied in the efficient management of genetic resource collections. For example clones of the same strain can be easily identified and eliminated.
Our primary goal was to use MLST to identify distinct groups of Sinorhizobium spp. and to estimate the relationships between these groups with the analysis of 10 chromosomal loci in 231 diverse strains of Medicago-nodulating rhizobia that had previously been characterized by MLEE (10). These results were then used in a follow-up study with several representative strains to examine any possibility of a relationship between chromosomal genotype and symbiotic phenotype.
|
|
|---|
Plant tests. Seeds of Medicago polymorpha were surface sterilized with concentrated H2SO4 for 3 min and were washed five times with sterile distilled water. The treated seeds were germinated on sterile water agar, seedlings were sown in sterile 50:50 (wt/wt) sand-vermiculite-filled Leonard jars (15), and 2 ml of modified arabinose-gluconate-grown late-log-phase broth cultures was used to inoculate each jar. The cultures tested for symbiosis were the type strains for the species S. meliloti (USDA 1002) and S. medicae (A321) and the MLEE group B strains CC169, M3, M16, M58, M75, M161, M173, and M254. Each treatment was prepared in three replications, and three jars without inoculated bacteria served as controls. The plants were grown in a greenhouse without supplemental lighting in two duplicate experiments for 30 and 42 days. The plants were uprooted, and the tops were cut off to determine nitrogenase activities as described by van Berkum et al. (29). Determinations for the concentration of ethylene in each chamber were as described by van Berkum and Sloger (27). The plant tops were dried at 60°C for 2 days to determine dry matter (6).
PCR primer design and PCR amplification of chromosomal loci. Loci for MLST analysis were chosen by referring to the complete genomic sequence of strain 1021 (4) to select 10 genes distributed across the chromosome (Table 1). The entire open reading frame of each locus was used to select primers that would amplify a portion of each gene between 200 and 500 bp in size by using the primer design software package Oligo Primer Analysis Software version 6.65 (Molecular Biology Insights, Inc., Cascade, CO). The oligonucleotides selected (Table 1) with Oligo were synthesized by Sigma-Genosys (The Woodlands, TX) and were received as dried preparations. Upon receipt the primers were dissolved with 10 mM Tris buffer, pH 8.0, to a final concentration of 1,000 pmol and were stored at 20°C. The PCRs for each locus were then optimized by using the FailSafe PCR PreMix selection kit (Epicentre, Madison, WI) and the thermal cycle protocol described by van Berkum and Fuhrmann (26) with an MJ Research PTC-225 Peltier thermal cycler (MJ Research Inc., Waltham, MA) using genomic DNAs of both USDA 1002 and A321 as templates. These 24 reactions were analyzed by horizontal agarose gel electrophoresis to select the FailSafe PCR system (Epicentre, Madison, WI) determined to be the most suitable for PCR amplification with the DNA preparations of all 231 strains used in this investigation. The presence of a single PCR product of the expected molecular size for each primer pair using each template was verified by horizontal gel electrophoresis. Each PCR product was then purified, especially to remove the PCR primers, by using the Ampure PCR purification system (Agincourt Bioscience Corporation, Beverly, MA).
|
View this table: [in a new window] |
TABLE 1. Primer sequences for the 10 loci used in MLST analysis of the chromosomes of 231 Medicago-nodulating rhizobia
|
Data analysis. A Microsoft Access database was created to compile the data collected with the 231 strains. The sequence length entered for each locus was identical, and the same alleles were identified using the software Sequence Comparator version 2.0.1, written by Keith Jolley. As additional alleles were identified they were assigned different allele numbers in the database. In each case, allelic variation was verified by confirming the sequence disparity using Genedoc (18) and then checking the electropherograms produced by the sequencing analysis to substantiate differences. In the case of ambiguities the sequencing analysis was repeated.
When the database was completed, the relationship function in Microsoft Access was used to create a "query" of the allelic allocation for each of the 10 loci across the 231 strains. The 10 entries across the strains were then listed according to an ascending value. The resulting tabulation of the data was then exported as a Microsoft Excel file to prepare data input files for the program Sequence Type Analysis and Recombinational Tests (START), written by Keith Jolley, University of Oxford (14); the statistics for population genetics computer program written by T. S. Whittam (23); and eBURST (12). START was used to determine both the allele and profile frequencies and to create an UPGMA (unweighted-pair group method average linkages) tree to portray the genetic relationships among the STs. Sequence types were classified as single-locus variants (SLVs), double-locus variants (DLVs), or singletons (STs differing at three or more loci). The computer program by T. S. Whittam was used to calculate the genetic distances among the alleles of each locus and to derive the mean value. The number of groups, the clonal complex of each group, and a population snapshot of the chromosomal variation were generated using eBURST. The null hypothesis of linkage equilibrium for the multilocus sequence data as defined by Maynard et al. (16) was evaluated using the program LIAN 3.0 (13).
Nucleotide sequence accession numbers. Sequences of the alleles for each locus were deposited with GenBank under accession numbers DQ423249 through DQ423361. The 16S rRNA gene sequences of strains 15B4, 74B15, 56A14, CC2013, S33, M7, and M58 were deposited with GenBank under accession numbers DQ423242 through DQ423248, respectively.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Numbers of alleles and genetic diversity across chromosomes of 231 Medicago-nodulating rhizobia
|
|
View this table: [in a new window] |
TABLE 2. STs with multiple strains
|
![]() View larger version (26K): [in a new window] |
FIG. 1. Linkage distance among 231 Medicago-nodulating rhizobia derived from allelic variation among 10 chromosomal loci. A matrix of the strain identification and the ST followed by the allele labels for each was used in the START (Sequence Type Analysis and Recombinational Tests, version 1.05) program to generate the UPGMA dendrogram. The program was written by Keith Jolley, University of Oxford (14). MLST group affiliation is indicated. Singletons are identified by an "s"; ST numbers followed by a blank belong to group 1.
|
![]() View larger version (32K): [in a new window] |
FIG. 2. Population snapshot of 231 Medicago-nodulating rhizobia derived from the allelic variation of 10 chromosomal loci. A matrix of the ST followed by the allele labels for each was used in eBURST (11) to generate a diagram of the evolutionary patterns among the strains. The snapshot was produced by setting the group definition to 0/10 genes. The sizes of the circles are related to the numbers of strains within each ST. The founder and cofounder genotypes are colored blue and red, respectively. Distances between STs indicated by connecting lines are arbitrary.
|
|
View this table: [in a new window] |
TABLE 4. Test of the null hypothesis of linkage equilibrium among 91 multilocus STs in a global collection of 231 Medicago-nodulating rhizobial strainsa
|
|
View this table: [in a new window] |
TABLE 5. Host plant and geographic origin of Medicago-nodulating rhizobia of the group 1 clonal complex
|
In addition to sequencing of the 10 loci for MLST, the 16S rRNA gene sequences were determined for a select set of strains. The 16S rRNA gene sequences of strains 56A14 (ST-5) and 74B15 (ST-51), representing groups 2 and 4, respectively, and of USDA 1002 (S. meliloti) were identical. The 16S rRNA gene sequences of group 6 strain S33 (ST-27) and of the singleton CC2013 (ST-71) were identical to each other, and both had a single nucleotide difference compared with the USDA 1002 allele. Strains M58 (ST-15) and M7 (ST-13) represent groups 3 and 5, respectively, and possessed a 16S rRNA gene sequence identical with that reported for strain A321 (21). The 16S rRNA gene sequence of these three strains had four base pair differences from the corresponding sequence of USDA 1002.
The symbiotic reactions of several strains that were placed in the cluster defined as S. medicae (21) were tested with M. polymorpha as the plant host. This was done because it has been reported that only S. medicae forms an effective symbiosis with this host legume species. If substantiated, this would represent a phenotypic trait that could be used to distinguish strains that were defined as S. medicae. The type strains for S. meliloti (USDA 1002) and S. medicae (A321) could be distinguished by the plant test since they produced an ineffective and an effective symbiotic response on M. polymorpha, respectively (Table 6 ). However, a mixed response was obtained with eight of the other strains placed in the S. medicae cluster: four were ineffective, one was intermediate in effectiveness, and the three remaining strains were fully effective, similar to strain A321 (Table 6). In a repeat test, M16, M161, M173, and M58 were again observed to nodulate and to establish an ineffective response with M. polymorpha while strains M254 and M3 were effective.
|
View this table: [in a new window] |
TABLE 6. Symbiotic reaction of several MLEE group B strains, identified as genospecies 2 or Sinorhizobium medicae (21), on Medicago polymorphaa
|
|
|
|---|
The presence of 91 STs among the 231 strains would indicate that there is extensive diversity in natural populations of Medicago-nodulating rhizobia. It seems likely that in this study only a small fraction of this diversity was revealed. It is anticipated that a much broader range of genotypes will be identified with an analysis of strains from other geographic regions.
In addition to its value for characterizing genetic diversity in bacterial populations, this method has potential practical applications as well. For example MLST data could be used for choosing unique strains for preservation in culture collections or for the selection of strains that might be useful in low-input sustainable agriculture.
One of the goals in this study was to determine if symbiotic host range and effectiveness might be associated with specific chromosomal genotypes. This question was based on the report of Rome et al. (21) and their conclusion that, while strains of S. medicae form effective symbioses with M. polymorpha, those of S. meliloti are ineffective with this legume species. The MLST analysis confirmed the separation of a cluster of strains that had previously been identified (by MLEE and other criteria) as belonging to the species S. medicae. However, from the subsequent plant tests it was apparent that there is no strict correlation between the chromosomal genotype of a strain and its symbiotic effectiveness on M. polymorpha. An explanation for this conclusion may be that the symbiotic determinants in S. meliloti and S. medicae reside on large extrachromosomal elements (e.g., pSymA) that may be subject to lateral transfer and recombination. The possibility that this may be the explanation for the results with the plant tests could be investigated with an MLST analysis of the extrachromosomal elements in the same strains.
|
|
|---|
, M., S. Sikora, S. Red
epovi
, and Z.
tafa. 2003. Genetic identification and symbiotic efficiency of an indigenous Sinorhizobium meliloti field population. Food Technol. Biotechnol. 41:69-75.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»