Application of Multilocus Sequence Typing To Study the Genetic Structure of Megaplasmids in Medicago-Nodulating Rhizobia

A multilocus sequence typing (MLST) analysis was used to examine the genetic structure and diversity within the two large extrachromosomal replicons in Medicago-nodulating rhizobia (Sinorhizobium meliloti and Sinorhizobium medicae). The allelic diversity within these replicons was high compared to the reported diversity within the corresponding chromosomes of the same strains (P. van Berkum et al., J. Bacteriol. 188:5570-5577, 2006). Also, there was strong localized linkage disequilibrium (LD) between certain pSymA loci: e.g., nodC and nifD. Although both of these observations could be explained by positive (or diversifying) selection by plant hosts, results of tests for positive selection did not provide consistent support for this hypothesis. The strong LD observed between the nodC and nifD genes could also be explained by their close proximity on the pSymA replicon. Evidence was obtained that some nodC alleles had a history of intragenic recombination, while other alleles of this locus had a history of intergenic recombination. Both types of recombination were associated with a decline in symbiotic competence with Medicago sativa as the host plant. The combined observations of LD between the nodC and nifD genes and intragenic recombination within one of these loci indicate that the symbiotic gene region on the pSymA plasmid has evolved as a clonal segment, which has been laterally transferred within the natural populations.Plants of the genus Medicago are legumes that often benefit from a mutualistic symbiosis with rhizobia. The most agriculturally significant species of rhizobia that nodulate these plants are Sinorhizobium meliloti (9) and Sinorhizobium medicae (22). Previously reported population genetic analyses of these bacteria have focused on the study of how allelic variants at multiple loci are distributed within and among natural populations (2, 3, 10, 26, 31, 32). This was also the focus of the present study, but it was extended by examining more loci in many more strains of both species of Sinorhizobium coupled with an analysis having a range of symbiotic genotypes. One goal was to determine if there were any obvious correlations between the megaplasmid genotypes observed and their symbiotic competence. A second goal was to determine if selection by their host plants may have influenced the evolution of their symbiotic relationships.The genes for symbiosis reside on the extrachromosomal replicons pSymA (1,354,226 nucleotides [nt]) and pSMED02 (1,245,408 nt) in the genomes of S. meliloti Rm1021 and S. medicae WSM419, respectively (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"AE006469","term_id":"25168258","term_text":"AE006469"}}AE006469 and {"type":"entrez-nucleotide","attrs":{"text":"CP000740","term_id":"150031715","term_text":"CP000740"}}CP000740, respectively). Besides these two plasmids, these two strains each harbor one other large extrachromosomal replicon, pSymB (1,683,333 nt) and pSMED01 (1,570,951 nt), respectively (GenBank accession no. {"type":"entrez-nucleotide","attrs":{"text":"AL591985","term_id":"30407151","term_text":"AL591985"}}AL591985 and {"type":"entrez-nucleotide","attrs":{"text":"CP000739","term_id":"150030273","term_text":"CP000739"}}CP000739, respectively).Multilocus sequence typing (MLST) (16) is a form of genomic indexing that is commonly used to study the population genetic structure and phylogenetic relatedness within diverse groups of bacteria. In this method, nucleotide sequences of a fixed set of common loci are obtained from a collection of strains, and polymorphic sites among these sequences are used to derive an allelic profile or sequence type (ST) for each genome. Comparisons of the resulting data can be used to infer phylogenetic relationships among the organisms in the sample population, and they also can be used to infer how evolutionary processes, such as recombination and selection, have shaped the genetic structure of the population. For example, levels of intergenic recombination among chromosomal genes in natural populations of Neisseria meningitidis reportedly are relatively high, while corresponding levels within populations of Staphylococcus aureus were low (28). Depending on the specific pairs of loci examined, the levels of linkage disequilibrium (LD) (a lack of intergenic recombination) among several chromosomally carried core genes of S. meliloti were reported to be generally moderate to high (26).The MLST approach has been used to confirm that the chromosomes of S. meliloti and S. medicae are sexually isolated (2, 3, 31) and to provide evidence that horizontal gene transfer (HGT) does occur between the symbiotic megaplasmids of these species (3, 32). It has also been used to demonstrate that levels of intergenic recombination, as indicated by linkage disequilibrium, differ between the three replicons of S. meliloti (26). Levels of intergenic recombination within the pSymB replicons of these strains are generally high, unlike the chromosomes and pSymA replicons within the same strains (26). Bailly et al. (3) hypothesized that the region of the pSymA plasmid that contains the nodulation (nod) genes is frequently transferred in natural populations. They also suggested that selective pressures from the host plant may have influenced both nod gene diversity and patterns of polymorphism across the entire nod gene region.In the present study, multilocus allelic variation of the two megaplasmids was examined among 231 Medicago-nodulating rhizobia that originated primarily from southwest Asia (10). Previously, 91 different chromosomal sequence types (STs) were identified among the same strains from sequence variation in 10 loci (31). This collection of strains had earlier been divided into two closely related groups based on results of multilocus enzyme electrophoresis (10), and this result was subsequently cited in support of separating the Medicago-nodulating rhizobia into the two species S. meliloti and S. medicae (22).The objectives of this study were (i) to use MLST to examine the genetic relationships within and among the large extrachromosomal replicons in S. meliloti and S. medicae, (ii) to estimate levels of intergenic and intragenic recombination in these replicons, (iii) to evaluate the nitrogen-fixing competence of representative symbiotic genotypes with Medicago sativa, and (iv) to determine whether positive (or diversifying) selection may have influenced the genetic structure of the megaplasmids.