Evolutionary Dynamics of Insertion Sequences in Relation to the Evolutionary Histories of the Chromosome and Symbiotic Plasmid Genes of Rhizobium etli Populations

Insertion sequences (IS) are mobile genetic elements that are distributed in many prokaryotes. In particular, in the genomes of the symbiotic nitrogen-fixing bacteria collectively known as rhizobia, IS are fairly abundant in plasmids or chromosomal islands that carry the genes needed for symbiosis. Here, we report an analysis of the distribution and genetic conservation of the IS found in the genome of Rhizobium etli CFN42 in a collection of 87 Rhizobium strains belonging to populations with different geographical origins. We used PCR to generate presence/absence profiles of the 39 IS found in R. etli CFN42 and evaluated whether the IS were located in consistent genomic contexts. We found that the IS from the symbiotic plasmid were frequently present in the analyzed strains, whereas the chromosomal IS were observed less frequently. We then examined the evolutionary dynamics of these strains based on a population genetic analysis of two chromosomal housekeeping genes (glyA and dnaB) and three symbiotic sequences (nodC and the two IS elements). Our results indicate that the IS contained within the symbiotic plasmid have a higher degree of genomic context conservation, lower nucleotide diversity and genetic differentiation, and fewer recombination events than the chromosomal housekeeping genes. These results suggest that the R. etli populations diverged recently in Mexico, that the symbiotic plasmid also had a recent origin, and that the IS elements have undergone a process of cyclic infection and expansion.Insertion sequences (IS) are the smallest transposable elements found in prokaryotes (usually less than 3 kb in size). They encode a transposase and may also encode small hypothetical proteins (4, 9). IS are distinguished by their ability to move within prokaryotic replicons, including both the chromosome and plasmids, and to copy themselves into various genomic sites. In this manner, IS elements can inactivate or alter the expression of adjacent genes (4). When IS occur in two or more identical copies within a genome, they can participate in various types of genetic rearrangements (e.g., duplications, inversions, and deletions), suggesting that IS may play an important role in the evolution of their hosts by promoting genomic plasticity (34). Due to these evolutionary dynamics, the diversity and distribution of IS elements differ greatly between taxa and even within strains of the same species (27).Various theories have been proposed to explain the evolution of IS elements in laboratory model strains and environmental bacterial populations (8, 18, 25, 29). Two main hypotheses seek to explain how these elements are maintained over the long term in their host genomes. The first proposes that they occasionally generate beneficial mutations and therefore may represent a selective advantage to their hosts (34). The second suggests that IS elements are genomic parasites that are maintained by their high rate of transposition and might be disseminated among different bacterial lineages by horizontal gene transfer (HGT). Data supporting the second hypothesis have shown that some IS elements may transpose at high rates upon entering a new host (42). After the initial infection, however, purifying selection may continuously remove these elements from the genome. Thus, IS may undergo an infection-expansion-extinction cycle that allows them to remain in different bacterial populations within the gene pool (42). These two hypotheses are not contradictory, and the evolutionary dynamics and distribution of IS may differ greatly depending on several factors, including (most notably) the rate of transposition and HGT, as well as selective pressures, population size, and the host''s habitat (6, 18, 21, 25, 27, 29).In the nitrogen-fixing symbiotic bacteria of the genera Rhizobium, Sinorhizobium, Mesorhizobium, Bradyrhizobium (of the alphaproteobacteria), Cupriavidus, and Burkholderia (of the betaproteobacteria), IS are particularly abundant in symbiotic plasmids (pSym) and symbiotic chromosomal islands (SI) (2, 12, 14, 19, 20, 43). SI and pSym include most of the genes needed to establish symbiosis in the roots of leguminous plants through nodule formation and nitrogen fixation (11). It is generally believed that these elements entered the rhizobial genomes through HGT (39, 40). Comparative genomic analyses have shown that both pSym and SI are highly variable, with the exception of a common set of genes encoding factors critical to nitrogen fixation (nif) and nodulation (nod) (5, 14). SI and pSym have been found to have lower GC contents and different codon usages than the corresponding chromosomal and nonsymbiotic plasmid sequences, suggesting that they were recently acquired by HGT.Some of these symbiotic elements, such as in pSym of Rhizobium etli CFN42 and the SI of Mesorhizobium loti, are conjugative and mobile (30, 32). Genomic analysis of R. etli CFN42 revealed the presence of 39 IS belonging to different families (14); they were found in the chromosome (11 IS); the 371-kb symbiotic plasmid (13 IS); the smaller 192-kb conjugative plasmid, p42a (13 IS); and two other plasmids, p42b and p42c (2 IS). Interestingly, this particular strain shows no evidence of IS disrupting open reading frames (ORFs) or having transpositional activity. However, another 42 incomplete IS may be found in the chromosome, pSym, and the conjugative plasmid; these incomplete sequences are truncated or contain stop codons in their coding sequences.Here, we focused on the dynamics and distribution of IS in different populations of the nitrogen-fixing symbiont R. etli. Since the maintenance of IS in bacterial species might depend on their transpositional activities and horizontal transfer rates, the identification of IS in the same genomic contexts across different strains of the same species could provide new insights into their persistence and divergence over short evolutionary periods. To examine the evolutionary dynamics of IS in natural populations of R. etli, we characterized the distributions, genomic contexts, and sequence variations of IS in isolates of R. etli from three populations with different origins, as well as in some other Rhizobium species. More specifically, we used PCR to generate presence/absence profiles of the 39 IS found in R. etli CFN42 in a collection of 87 strains representing different geographical sites and a gradient of domestication of the bacterial host, the common bean (Phaseolus vulgaris). We also evaluated whether the IS were conserved in the same genomic context relative to their position in R. etli CFN42 and determined the nucleotide sequences of two IS found in most of the isolates. Several population genetic tests applied to these IS, another pSym gene (nodC), and two chromosomal housekeeping genes (dnaB and glyA) suggest that these two IS elements have been inherited vertically and represent recent components of the R. etli gene pool. Finally, the present study strongly suggests that symbiotic plasmids have a recent origin within the R. etli populations.     

Conserved Plasmid Hydrogen-Uptake (hup)-Specific Sequences within Hup+Rhizobium leguminosarum Strains

Thirteen Rhizobium leguminosarum strains previously reported as H₂-uptake hydrogenase positive (Hup⁺) or negative (Hup⁻) were analyzed for the presence and conservation of DNA sequences homologous to cloned Bradyrhizobium japonicum hup-specific DNA from cosmid pHU1 (M. A. Cantrell, R. A. Haugland, and H. J. Evans, Proc. Natl. Acad. Sci. USA 80:181-185, 1983). The Hup phenotype of these strains was reexamined by determining hydrogenase activity induced in bacteroids from pea nodules. Five strains, including H₂ oxidation-ATP synthesis-coupled and -uncoupled strains, induced significant rates of H₂-uptake hydrogenase activity and contained DNA sequences homologous to three probe DNA fragments (5.9-kilobase [kb] HindIII, 2.9-kb EcoRI, and 5.0-kb EcoRI) from pHU1. The pattern of genomic DNA HindIII and EcoRI fragments with significant homology to each of the three probes was identical in all five strains regardless of the H₂-dependent ATP generation trait. The restriction fragments containing the homology totalled about 22 kb of DNA common to the five strains. In all instances the putative hup sequences were located on a plasmid that also contained nif genes. The molecular sizes of the identified hup-sym plasmids ranged between 184 and 212 megadaltons. No common DNA sequences homologous to B. japonicum hup DNA were found in genomic DNA from any of the eight remaining strains showing no significant hydrogenase activity in pea bacteroids. These results suggest that the identified DNA region contains genes essential for hydrogenase activity in R. leguminosarum and that its organization is highly conserved within Hup⁺ strains in this symbiotic species.     

Characterization of IntA,a Bidirectional Site-Specific Recombinase Required for Conjugative Transfer of the Symbiotic Plasmid of Rhizobium etli CFN42

Site-specific recombination occurs at short specific sequences, mediated by the cognate recombinases. IntA is a recombinase from Rhizobium etli CFN42 and belongs to the tyrosine recombinase family. It allows cointegration of plasmid p42a and the symbiotic plasmid via site-specific recombination between attachment regions (attA and attD) located in each replicon. Cointegration is needed for conjugative transfer of the symbiotic plasmid. To characterize this system, two plasmids harboring the corresponding attachment sites and intA were constructed. Introduction of these plasmids into R. etli revealed IntA-dependent recombination events occurring at high frequency. Interestingly, IntA promotes not only integration, but also excision events, albeit at a lower frequency. Thus, R. etli IntA appears to be a bidirectional recombinase. IntA was purified and used to set up electrophoretic mobility shift assays with linear fragments containing attA and attD. IntA-dependent retarded complexes were observed only with fragments containing either attA or attD. Specific retarded complexes, as well as normal in vivo recombination abilities, were seen even in derivatives harboring only a minimal attachment region (comprising the 5-bp central region flanked by 9- to 11-bp inverted repeats). DNase I-footprinting assays with IntA revealed specific protection of these zones. Mutations that disrupt the integrity of the 9- to 11-bp inverted repeats abolish both specific binding and recombination ability, while mutations in the 5-bp central region severely reduce both binding and recombination. These results show that IntA is a bidirectional recombinase that binds to att regions without requiring neighboring sequences as enhancers of recombination.     

Effect of Sym Plasmid Curing on Symbiotic Effectiveness in Rhizobium fredii

A mutant, USDA 206C, of Rhizobium fredii USDA 206 was obtained by passage on acridine plates. This mutant was cured of its 197-megadalton Sym plasmid but retained its symbiotic effectiveness. Multiple plasmid and chromosomally borne nif gene copies have previously been shown in R. fredii USDA 206. HindIII and EcoRI restriction enzyme digests of plasmid and total DNA showed that at least two nif gene copies are probably missing in USDA 206C. To compare the symbiotic effectiveness of USDA 206 and USDA 206C, plant tests were carried out. Statistically significant differences were obtained for nodule number, nodule mass, nitrogenase activity per plant, nitrogenase specific activity, and total plant dry weight. There was an apparent correlation between loss of Sym plasmidborne nif gene copies and reduction of overall symbiotic effectiveness. Delayed nodulation by strain USDA 206C relative to strain USDA 206 also indicated an association with the loss of plasmidborne nodulation functions and the reduced symbiotic effectiveness of strain USDA 206C.     

Diversification of DNA sequences in the symbiotic genome of Rhizobium etli

Bacteria of the genus Rhizobium and related genera establish nitrogen-fixing symbioses with the roots of leguminous plants. The genetic elements that participate in the symbiotic process are usually compartmentalized in the genome, either as independent replicons (symbiotic plasmids) or as symbiotic regions or islands in the chromosome. The complete nucleotide sequence of the symbiotic plasmid of Rhizobium etli model strain CFN42, symbiont of the common bean plant, has been reported. To better understand the basis of DNA sequence diversification of this symbiotic compartment, we analyzed the distribution of single-nucleotide polymorphisms in homologous regions from different Rhizobium etli strains. The distribution of polymorphisms is highly asymmetric in each of the different strains, alternating regions containing very few changes with regions harboring an elevated number of substitutions. The regions showing high polymorphism do not correspond with discrete genetic elements and are not the same in the different strains, indicating that they are not hypervariable regions of functional genes. Most interesting, some highly polymorphic regions share exactly the same nucleotide substitutions in more than one strain. Furthermore, in different regions of the symbiotic compartment, different sets of strains share the same substitutions. The data indicate that the majority of nucleotide substitutions are spread in the population by recombination and that the contribution of new mutations to polymorphism is relatively low. We propose that the horizontal transfer of homologous DNA segments among closely related organisms is a major source of genomic diversification.     

Effect of Plasmid pIJ1008 from Rhizobium leguminosarum on Symbiotic Function of Rhizobium meliloti

Plasmid pIJ1008, which carries determinants for uptake hydrogenase (Hup) activity, was transferred from Rhizobium leguminosarum to Rhizobium meliloti without impairing the capacity of the latter species to form root nodules on alfalfa. The plasmid was still present in rhizobia reisolated from the root nodules of 12 different alfalfa cultivars, but only low levels of Hup activity were detected in alfalfa.     

Genetic Determinants of Swarming in Rhizobium etli

Swarming motility is considered to be a social phenomenon that enables groups of bacteria to move coordinately atop solid surfaces. The differentiated swarmer cell population is embedded in an extracellular slime layer, and the phenomenon has previously been linked with biofilm formation and virulence. The gram-negative nitrogen-fixing soil bacterium Rhizobium etli CNPAF512 was previously shown to display swarming behavior on soft agar plates. In a search for novel genetic determinants of swarming, a detailed analysis of the swarming behavior of 700 miniTn5 mutants of R. etli was performed. Twenty-four mutants defective in swarming or displaying abnormal swarming patterns were identified and could be divided into three groups based on their swarming pattern. Fourteen mutants were completely swarming deficient, five mutants showed an atypical swarming pattern with no completely smooth edge and local extrusions, and five mutants displayed an intermediate swarming phenotype. Sequence analysis of the targeted genes indicated that the mutants were likely affected in quorum-sensing, polysaccharide composition or export, motility, and amino acid and polyamines metabolism. Several of the identified mutants displayed a reduced symbiotic nitrogen fixation activity.     

Three Replicons of Rhizobium sp. Strain NGR234 Harbor Symbiotic Gene Sequences

Rhizobium sp. strain NGR234 contains three replicons: the symbiotic plasmid or pNGR234a, a megaplasmid (pNGR234b), and the chromosome. Symbiotic gene sequences not present in pNGR234a were analyzed by hybridization. DNA sequences homologous to the genes fixLJKNOPQGHIS were found on the chromosome, while sequences homologous to nodPQ and exoBDFLK were found on pNGR234b.     

Conserved Plasmid Hydrogen-Uptake (hup)-Specific Sequences within HupRhizobium leguminosarum Strains

Thirteen Rhizobium leguminosarum strains previously reported as H(2)-uptake hydrogenase positive (Hup) or negative (Hup) were analyzed for the presence and conservation of DNA sequences homologous to cloned Bradyrhizobium japonicum hup-specific DNA from cosmid pHU1 (M. A. Cantrell, R. A. Haugland, and H. J. Evans, Proc. Natl. Acad. Sci. USA 80:181-185, 1983). The Hup phenotype of these strains was reexamined by determining hydrogenase activity induced in bacteroids from pea nodules. Five strains, including H(2) oxidation-ATP synthesis-coupled and -uncoupled strains, induced significant rates of H(2)-uptake hydrogenase activity and contained DNA sequences homologous to three probe DNA fragments (5.9-kilobase [kb] HindIII, 2.9-kb EcoRI, and 5.0-kb EcoRI) from pHU1. The pattern of genomic DNA HindIII and EcoRI fragments with significant homology to each of the three probes was identical in all five strains regardless of the H(2)-dependent ATP generation trait. The restriction fragments containing the homology totalled about 22 kb of DNA common to the five strains. In all instances the putative hup sequences were located on a plasmid that also contained nif genes. The molecular sizes of the identified hup-sym plasmids ranged between 184 and 212 megadaltons. No common DNA sequences homologous to B. japonicum hup DNA were found in genomic DNA from any of the eight remaining strains showing no significant hydrogenase activity in pea bacteroids. These results suggest that the identified DNA region contains genes essential for hydrogenase activity in R. leguminosarum and that its organization is highly conserved within Hup strains in this symbiotic species.     

Genomic basis of symbiovar mimosae in Rhizobium etli

Background

Symbiosis genes (nod and nif) involved in nodulation and nitrogen fixation in legumes are plasmid-borne in Rhizobium. Rhizobial symbiotic variants (symbiovars) with distinct host specificity would depend on the type of symbiosis plasmid. In Rhizobium etli or in Rhizobium phaseoli, symbiovar phaseoli strains have the capacity to form nodules in Phaseolus vulgaris while symbiovar mimosae confers a broad host range including different mimosa trees.

Results

We report on the genome of R. etli symbiovar mimosae strain Mim1 and its comparison to that from R. etli symbiovar phaseoli strain CFN42. Differences were found in plasmids especially in the symbiosis plasmid, not only in nod gene sequences but in nod gene content. Differences in Nod factors deduced from the presence of nod genes, in secretion systems or ACC-deaminase could help explain the distinct host specificity. Genes involved in P. vulgaris exudate uptake were not found in symbiovar mimosae but hup genes (involved in hydrogen uptake) were found. Plasmid pRetCFN42a was partially contained in Mim1 and a plasmid (pRetMim1c) was found only in Mim1. Chromids were well conserved.

Conclusions

The genomic differences between the two symbiovars, mimosae and phaseoli may explain different host specificity. With the genomic analysis presented, the term symbiovar is validated. Furthermore, our data support that the generalist symbiovar mimosae may be older than the specialist symbiovar phaseoli.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-575) contains supplementary material, which is available to authorized users.     

Fermentative and aerobic metabolism in Rhizobium etli.

Strains of Rhizobium etli, Rhizobium meliloti, and Rhizobium tropici decreased their capacity to grow after successive subcultures in minimal medium, with a pattern characteristic for each species. During the growth of R. etli CE 3 in minimal medium (MM), a fermentation-like response was apparent: the O2 content was reduced and, simultaneously, organic acids and amino acids were excreted and poly-beta-hydroxybutyrate (PHB) was accumulated. Some of the organic acids excreted into the medium were tricarboxylic acid (TCA) cycle intermediates, and, concomitantly, the activities of several TCA cycle and auxiliary enzymes decreased substantially or became undetectable. Optimal and sustained growth and a low PHB content were found in R. etli CE 3 when it was grown in MM inoculated at a low cell density with O2 maintained at 20% or with the addition of supplements that have an effect on the supply of substrates for the TCA cycle. In the presence of supplements such as biotin or thiamine, no amino acids were excreted and the organic acids already excreted into the medium were later reutilized. Levels of enzyme activities in cells from supplemented cultures indicated that carbon flux through the TCA cycle was maintained, which did not happen in MM. It is proposed that the fermentative state in Rhizobium species is triggered by a cell density signal that results in the regulation of some of the enzymes responsible for the flux of carbon through the TCA cycle and that this in turn determines how much carbon is available for the synthesis and accumulation of PHB. The fermentative state of free-living Rhizobium species may be closely related to the metabolism that these bacteria express during symbiosis.     

Multiresistance genes of Rhizobium etli CFN42

Multidrug efflux pumps of bacteria are involved in the resistance to various antibiotics and toxic compounds. In Rhizobium etli, a mutualistic symbiont of Phaseolus vulgaris (bean), genes resembling multidrug efflux pump genes were identified and designated rmrA and rmrB. rmrA was obtained after the screening of transposon-generated fusions that are inducible by bean-root released flavonoids. The predicted gene products of rmrAB shared significant homology to membrane fusion and major facilitator proteins, respectively. Mutants of rmrA formed on average 40% less nodules in bean, while mutants of rmrA and rmrB had enhanced sensitivity to phytoalexins, flavonoids, and salicylic acid, compared with the wild-type strain. Multidrug resistance genes emrAB from Escherichia coli complemented an rmrA mutant from R. etli for resistance to high concentrations of naringenin.     

Attenuation of Symbiotic Effectiveness by Rhizobium meliloti SAF22 Related to the Presence of a Cryptic Plasmid

Several wild-type strains of Rhizobium meliloti isolated from alfalfa nodules exhibited different plasmid profiles, yet did not differ in growth rate in yeast-mannitol medium, utilization of 43 different carbon sources, intrinsic resistance to 14 antibiotics, or detection of 16 enzyme activities. In contrast, three measures of effectiveness in symbiotic nitrogen fixation with alfalfa (shoot length, dry weight, and nitrogen content) indicated that R. meliloti SAF22, whose plasmid profile differs from those of the other strains tested, is significantly less effective than other wild-type strains in symbiotic nitrogen fixation. Light microscopy of nodules infected with strain SAF22 showed an abnormal center of nitrogen fixation zone III, with bacteria occupying a smaller portion of the infected host cells and vacuoles occupying a significantly larger portion of adjacent uninfected host cells. In contrast, the effective nodules infected with other wild types or plasmid pRmSAF22c-cured segregants of SAF22 did not display this cytological abnormality.     

豌豆根瘤菌共生基因pJB5JI与华癸中生根瘤菌7653R共生质粒的相互关系

华癸中生根瘤菌(Mesorhizobium huakuii)7653R是分离自我国南方水稻田的一株根瘤菌,含有2个内源质粒:p7653Ra和p7653Rb,其中7653Rb是共生质粒.通过Tn5-sacB的插入方法来消除质粒,获得7653Rb消除的突变株7653RD.将豌豆根瘤菌T83K3的共生质粒pJB5JI导入7653R和7653RD中,盆栽结果表明含有pJB5JI的转移接合子7653R-197的竞争结瘤能力和共生固氮能力均高于7653R.pJB5JI不能恢复7653RD在紫云英上的结瘤能力.含有pJB5JI的7653RD可以在豌豆上结无效瘤,表明pJB5JI可以在7653R的染色体背景下表达其功能.对转移接合子中的质粒稳定性进行检测,结果表明pJB5JI在人工传代的情况下可以稳定存在,但经过共生之后发生了遗传分离,对转移接合子和出发菌株及分离菌株进行kan基因的PCR扩增,除了受体菌外其他菌株都可得到PCR产物,由此推测,pJB5JI可能部分或全部整合到了受体菌的染色体基因组中.