Symbiotic plasmid rearrangement in Rhizobium leguminosarum bv. viciae VF39SM

Closed genomes of commercial inoculant rhizobia provide a blueprint for management of legume inoculation

Rhizobia are soil bacteria capable of establishing symbiosis within legume root nodules, where they reduce atmospheric N2 into ammonia and supply it to the plant for growth. Australian soils often lack rhizobia compatible with introduced agricultural legumes, so inoculation with exotic strains has become a common practice for over 50 years. While extensive research has assessed the N2-fixing capabilities of these inoculants, their genomics, taxonomy, and core and accessory gene phylogeny are poorly characterized. Furthermore, in some cases, inoculant strains have been developed from isolations made in Australia. It is unknown whether these strains represent naturalized exotic organisms, native rhizobia with a capacity to nodulate introduced legumes, or recombinant strains arising from horizontal transfer between introduced and native bacteria. Here, we describe the complete, closed genome sequences of 42 Australian commercial rhizobia. These strains span the genera, Bradyrhizobium, Mesorhizobium, Methylobacterium, Rhizobium, and Sinorhizobium, and only 23 strains were identified to species level. Within inoculant strain genomes, replicon structure and location of symbiosis genes were consistent with those of model strains for each genus, except for Rhizobium sp. SRDI969, where the symbiosis genes are chromosomally encoded. Genomic analysis of the strains isolated from Australia showed they were related to exotic strains, suggesting that they may have colonized Australian soils following undocumented introductions. These genome sequences provide the basis for accurate strain identification to manage inoculation and identify the prevalence and impact of horizontal gene transfer (HGT) on legume productivity.

IMPORTANCE

Inoculation of cultivated legumes with exotic rhizobia is integral to Australian agriculture in soils lacking compatible rhizobia. The Australian inoculant program supplies phenotypically characterized high-performing strains for farmers but in most cases, little is known about the genomes of these rhizobia. Horizontal gene transfer (HGT) of symbiosis genes from inoculant strains to native non-symbiotic rhizobia frequently occurs in Australian soils and can impact the long-term stability and efficacy of legume inoculation. Here, we present the analysis of reference-quality genomes for 42 Australian commercial rhizobial inoculants. We verify and classify the genetics, genome architecture, and taxonomy of these organisms. Importantly, these genome sequences will facilitate the accurate strain identification and monitoring of inoculants in soils and plant nodules, as well as enable detection of horizontal gene transfer to native rhizobia, thus ensuring the efficacy and integrity of Australia's legume inoculation program.

Diversification of plasmids in a genus of pathogenic and nitrogen-fixing bacteria

Members of the agrobacteria-rhizobia complex (ARC) have multiple and diverse plasmids. The extent to which these plasmids are shared and the consequences of their interactions are not well understood. We extracted over 4000 plasmid sequences from 1251 genome sequences and constructed a network to reveal interactions that have shaped the evolutionary histories of oncogenic virulence plasmids. One newly discovered type of oncogenic plasmid is a mosaic with three incomplete, but complementary and partially redundant virulence loci. Some types of oncogenic plasmids recombined with accessory plasmids or acquired large regions not known to be associated with pathogenicity. We also identified two classes of partial virulence plasmids. One class is potentially capable of transforming plants, but not inciting disease symptoms. Another class is inferred to be incomplete and non-functional but can be found as coresidents of the same strain and together are predicted to confer pathogenicity. The modularity and capacity for some plasmids to be transmitted broadly allow them to diversify, convergently evolve adaptive plasmids and shape the evolution of genomes across much of the ARC. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Why are rhizobial symbiosis genes mobile?

Rhizobia are one of the most important and best studied groups of bacterial symbionts. They are defined by their ability to establish nitrogen-fixing intracellular infections within plant hosts. One surprising feature of this symbiosis is that the bacterial genes required for this complex trait are not fixed within the chromosome, but are encoded on mobile genetic elements (MGEs), namely plasmids or integrative and conjugative elements. Evidence suggests that many of these elements are actively mobilizing within rhizobial populations, suggesting that regular symbiosis gene transfer is part of the ecology of rhizobial symbionts. At first glance, this is counterintuitive. The symbiosis trait is highly complex, multipartite and tightly coevolved with the legume hosts, while transfer of genes can be costly and disrupt coadaptation between the chromosome and the symbiosis genes. However, horizontal gene transfer is a process driven not only by the interests of the host bacterium, but also, and perhaps predominantly, by the interests of the MGEs that facilitate it. Thus understanding the role of horizontal gene transfer in the rhizobium-legume symbiosis requires a 'mobile genetic element's-eye view' on the ecology and evolution of this important symbiosis. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Genetic Variation in Host-Specific Competitiveness of the Symbiont Rhizobium leguminosarum Symbiovar viciae

Legumes of the Fabeae tribe form nitrogen-fixing root nodules resulting from symbiotic interaction with the soil bacteria Rhizobium leguminosarum symbiovar viciae (Rlv). These bacteria are all potential symbionts of the Fabeae hosts but display variable partner choice when co-inoculated in mixture. Because partner choice and symbiotic nitrogen fixation mostly behave as genetically independent traits, the efficiency of symbiosis is often suboptimal when Fabeae legumes are exposed to natural Rlv populations present in soil. A core collection of 32 Rlv bacteria was constituted based on the genomic comparison of a collection of 121 genome sequences, representative of known worldwide diversity of Rlv. A variable part of the nodD gene sequence was used as a DNA barcode to discriminate and quantify each of the 32 bacteria in mixture. This core collection was co-inoculated on a panel of nine genetically diverse Pisum sativum, Vicia faba, and Lens culinaris genotypes. We estimated the relative Early Partner Choice (EPC) of the bacteria with the Fabeae hosts by DNA metabarcoding on the nodulated root systems. Comparative genomic analyses within the bacterial core collection identified molecular markers associated with host-dependent symbiotic partner choice. The results revealed emergent properties of rhizobial populations. They pave the way to identify genes related to important symbiotic traits operating at this level.

Defining the requirements for the conjugative transfer of Rhizobium leguminosarum plasmid pRleVF39b

Rhizobium leguminosarum strain VF39 contains a plasmid, pRleVF39b, which encodes a distinctive type of conjugation system (rhizobial type IVa) that is relatively widespread among rhizobial genomes. The cluster of genes encoding the transfer functions lacks orthologs to genes such as traCD, traF and traB, but contains 15 conserved genes of unknown function. We determined the importance of these genes in conjugation by constructing marked and unmarked mutations in each gene, and established that six genes, now designated trcA-F, played a significant role in plasmid transfer. Like the relaxase gene, traA, and the genes encoding the MPF system (trb genes), five of these genes, located in two divergently transcribed operons, are regulated by the Xre family repressor TrbR. The other gene, trcF encodes a protein with similarity to histidinol phosphatases, and its role in conjugation is unclear, but mutations in trcF are severely impaired for conjugation. TrcF does not play a role in regulation of other conjugation genes.

Exploring the evolutionary dynamics of Rhizobium plasmids through bipartite network analysis

The genus Rhizobium usually has a multipartite genome architecture with a chromosome and several plasmids, making these bacteria a perfect candidate for plasmid biology studies. As there are no universally shared genes among typical plasmids, network analyses can complement traditional phylogenetics in a broad-scale study of plasmid evolution. Here, we present an exhaustive analysis of 216 plasmids from 49 complete genomes of Rhizobium by constructing a bipartite network that consists of two classes of nodes, the plasmids and homologous protein families that connect them. Dissection of the network using a hierarchical clustering strategy reveals extensive variety, with 34 homologous plasmid clusters. Four large clusters including one cluster of symbiotic plasmids and two clusters of chromids carrying some truly essential genes are widely distributed among Rhizobium. In contrast, the other clusters are quite small and rare. Symbiotic clusters and rare accessory clusters are exogenetic and do not appear to have co-evolved with the common accessory clusters; the latter ones have a large coding potential and functional complementarity for different lifestyles in Rhizobium. The bipartite network also provides preliminary evidence of Rhizobium plasmid variation and formation including genetic exchange, plasmid fusion and fission, exogenetic plasmid transfer, host plant selection, and environmental adaptation.

Rhizobial extrachromosomal replicon variability, stability and expression in natural niches

In bacteria, niche adaptation may be determined by mobile extrachromosomal elements. A remarkable characteristic of Rhizobium and Ensifer (Sinorhizobium) but also of Agrobacterium species is that almost half of the genome is contained in several large extrachromosomal replicons (ERs). They encode a plethora of functions, some of them required for bacterial survival, niche adaptation, plasmid transfer or stability. In spite of this, plasmid loss is common in rhizobia upon subculturing. Rhizobial gene-expression studies in plant rhizospheres with novel results from transcriptomic analysis of Rhizobium phaseoli in maize and Phaseolus vulgaris roots highlight the role of ERs in natural niches and allowed the identification of common extrachromosomal genes expressed in association with plant rootlets and the replicons involved.

Rhizobial communities in symbiosis with legumes: genetic diversity, competition and interactions with host plants

The term ‘Rhizobium-legume symbiosis’ refers to numerous plant-bacterial interrelationships. Typically, from an evolutionary perspective, these symbioses can be considered as species-to-species interactions, however, such plant-bacterial symbiosis may also be viewed as a low-scale environmental interplay between individual plants and the local microbial population. Rhizobium-legume interactions are therefore highly important in terms of microbial diversity and environmental adaptation thereby shaping the evolution of plant-bacterial symbiotic systems. Herein, the mechanisms underlying and modulating the diversity of rhizobial populations are presented. The roles of several factors impacting successful persistence of strains in rhizobial populations are discussed, shedding light on the complexity of rhizobial-legume interactions.

Recovery of nonpathogenic mutant bacteria from tumors caused by several Agrobacterium tumefaciens strains: a frequent event?

We have evaluated the interaction that bacterial genotypes and plant hosts have with the loss of pathogenicity in tumors, using seven Agrobacterium tumefaciens strains inoculated on 12 herbaceous and woody hosts. We performed a screening of the agrobacteria present inside the tumors, looking for nonpathogenic strains, and found a high variability of those strains in this niche. To verify the origin of the putative nonpathogenic mutant bacteria, we applied an efficient, reproducible, and specific randomly amplified polymorphic DNA analysis method. In contrast with previous studies, we recovered a very small percentage (0.01%) of nonpathogenic strains that can be considered true mutants. Of 5,419 agrobacterial isolates examined, 662 were nonpathogenic in tomato, although only 7 (from pepper and tomato tumors induced by two A. tumefaciens strains) could be considered to derive from the inoculated strain. Six mutants were affected in the transferred DNA (T-DNA) region; one of them contained IS426 inserted into the iaaM gene, whereas the whole T-DNA region was apparently deleted in three other mutants, and the virulence of the remaining two mutants was fully restored with the T-DNA genes as well. The plasmid profile was altered in six of the mutants, with changes in the size of the Ti plasmid or other plasmids and/or the acquisition of new plasmids. Our results also suggest that the frequent occurrence of nonpathogenic clones in the tumors is probably due to the preferential growth of nonpathogenic agrobacteria, of either endophytic or environmental origin, but different from the bacterial strain inducing the tumor.

Occurrence of sep insecticidal toxin complex genes in Serratia spp. and Yersinia frederiksenii

Some strains of Serratia entomophila and S. proteamaculans cause amber disease of the grass grub Costelytra zealandica (Coleoptera: Scarabaeidae). Three genes required for virulence, sepABC, are located on a large plasmid, pADAP. Sequence analysis suggests that the sepABC gene cluster may be part of a horizontally mobile region. This study presents evidence for the putative mobility of the sep genes of pADAP. Southern blot analysis showed that orthologues of the sep genes reside on plasmids within S. entomophila, S. liquefaciens, S. proteamaculans, and a plasmid from Yersinia frederiksenii. Three plasmids hybridized to the pADAP sep virulence-associated region but not the pADAP replication and conjugation regions. Subsequent DNA sequence analysis of the Y. frederiksenii sep-like genes, designated tcYF1 and tcYF2, showed that they had 88% and 87% DNA identity to sepA and sepB, respectively. These results indicate that the sep genes are part of a discrete horizontally mobile region.

Phase variation and genomic architecture changes in Azospirillum

The plant growth-promoting rhizobacterium Azospirillum lipoferum 4B generates in vitro at high frequency a stable nonswimming phase variant designated 4V(I), which is distinguishable from the wild type by the differential absorption of dyes. The frequency of variants generated by a recA mutant of A. lipoferum 4B was increased up to 10-fold. The pleiotropic modifications characteristic of the phase variant are well documented, but the molecular processes involved are unknown. Here, the objective was to assess whether genomic rearrangements take place during phase variation of strain 4B. The random amplified polymorphic DNA (RAPD) profiles of strains 4B and 4V(I) differed. RAPD fragments observed only with the wild type were cloned, and three cosmids carrying the corresponding fragments were isolated. The three cosmids hybridized with a 750-kb plasmid and pulse-field gel electrophoresis analysis revealed that this replicon was missing in the 4V(I) genome. The same rearrangements took place during phase variation of 4BrecA. Large-scale genomic rearrangements during phase variation were demonstrated for two additional strains. In Azospirillum brasilense WN1, generation of stable variants was correlated with the disappearance of a replicon of 260 kb. For Azospirillum irakense KBC1, the variant was not stable and coincided with the formation of a new replicon, whereas the revertant recovered the parental genomic architecture. This study shows large-scale genomic rearrangements in Azospirillum strains and correlates them with phase variation.

Screening possible mechanisms mediating cadmium resistance in Rhizobium leguminosarum bv. viciae isolated from contaminated Portuguese soils

Environment heavy-metal contamination is now widespread. Soils may become contaminated from a variety of anthropogenic sources, such as smelters, mining, industry, and application of metal-containing pesticides and fertilizers. Soil microorganisms are very sensitive to moderate heavy-metal concentrations. Therefore, the present work was designed to screen possible mechanisms involved in Rhizobium's Cd resistance; with this purpose, we determined the tolerance levels of several isolates originated from sites with different heavy-metal contamination. Whole-cell-soluble proteins and plasmid profiles were analyzed. We also determined Cd cell concentrations and lipopolysaccharide (LPS) amounts. Results showed different tolerances among Rhizobium isolates; according to their maximum resistance level, isolates were divided in four groups: sensitive (0-125 microM CdCl(2)), moderately tolerant (125-210 microM CdCl(2)), tolerant (250-500 microM CdCl(2)), and extremely tolerant (> or =750 microM CdCl(2)). Intracellular Cd concentrations were lower when compared to wall-bound Cd. Unexpectedly, extremely tolerant isolates accumulated higher levels of metal, suggesting the presence of intracellular agents that prevent metal interfering with important metabolic pathways. The electrophoretic patterns of whole-cell-soluble proteins evidenced cadmium as an inducer of protein metabolism alterations, which were more evident in some polypeptides. Plasmid profiles also showed differences; most tolerant isolates presented two plasmids with molecular weights of 485 and 415 kb, indicating that extrachromosomal DNA may be involved in cadmium resistance. LPS showed to be a common mechanism of resistance. However, the degree of tolerance conferred by LPS is not enough to support tolerance to the higher levels of stress imposed. Presence of other resistance mechanisms is currently being investigated.

Apparent incompatibility of plasmid pSfrYC4b of Sinorhizobium fredii with two different plasmids in another strain

Sinorhizobium fredii YC4B is a spontaneous mutant derivative of strain YC4 that is unable to nodulate soybeans. The second-largest plasmid of strain YC4B, termed pSfrYC4b (810 kb), was transferred to S. fredii HN01SR, a strain which contains three large indigenous plasmids (pSfrHN01a, pSfrHN01b and pSfrHN01c). Surprisingly, two stable indigenous plasmids (pSfrHN01a and pSfrHN01b) of strain HN01SR were cured simultaneously by the introduction of pSfrYC4b. Furthermore, a novel, unstable plasmid (pHY4) became visible in agarose gels. The electrophoretic mobility of plasmid pHY4 was slower than that shown by the cured plasmids, indicating that the molecular weight of the former is higher than that of plasmids pSfrYC4b and pSfrHN01b. Replication gene repC-like sequences were detected by polymerase chain reaction (PCR) on pSfrHN01a and pSfrYC4b, but not on pSfrHN01b. Sau3AI and PstI restriction patterns of the PCR-amplified repC-like sequences from HN01SR and YC4B were very similar.

Transfer of the symbiotic plasmid of Rhizobium etli CFN42 requires cointegration with p42a, which may be mediated by site-specific recombination

Plasmid p42a from Rhizobium etli CFN42 is self-transmissible and indispensable for conjugative transfer of the symbiotic plasmid (pSym). Most pSym transconjugants also inherit p42a. pSym transconjugants that lack p42a always contain recombinant pSyms, which we designated RpSyms*. RpSyms* do not contain some pSym segments and instead have p42a sequences, including the replication and transfer regions. These novel recombinant plasmids are compatible with wild-type pSym, incompatible with p42a, and self-transmissible. The symbiotic features of derivatives simultaneously containing a wild-type pSym and an RpSym* were analyzed. Structural analysis of 10 RpSyms* showed that 7 shared one of the two pSym-p42a junctions. Sequencing of this common junction revealed a 53-bp region that was 90% identical in pSym and p42a, including a 5-bp central region flanked by 9- to 11-bp inverted repeats reminiscent of bacterial and phage attachment sites. A gene encoding an integrase-like protein (intA) was localized downstream of the attachment site on p42a. Mutation or the absence of intA abolished pSym transfer from a recA mutant donor. Complementation with the wild-type intA gene restored transfer of pSym. We propose that pSym-p42a cointegration is required for pSym transfer; cointegration may be achieved either through homologous recombination among large reiterated sequences or through IntA-mediated site-specific recombination between the attachment sites. Cointegrates formed through the site-specific system but resolved through RecA-dependent recombination or vice versa generate RpSyms*. A site-specific recombination system for plasmid cointegration is a novel feature of these large plasmids and implies that there is unique regulation which affects the distribution of pSym in nature due to the role of the cointegrate in conjugative transfer.

Differentiation of plasmids in marine diazotroph assemblages determined by randomly amplified polymorphic DNA analysis

Nitrogen fixation by diazotrophic bacteria is a significant source of new nitrogen in salt marsh ecosystems. Recent studies have characterized the physiological and phylogenetic diversity of oxygen-utilizing diazotrophs isolated from the rhizoplanes of spatially separated intertidal macrophyte habitats. However, there is a paucity of information regarding the traits encoded by and the diversity of plasmids occurring in this key ecological functional group. Five-hundred and twenty-one isolates cultivated from the rhizoplanes of Juncus roemarianus, Spartina patens and different growth forms (short-form and tall-form) of Spartina alterniflora were screened for the presence of plasmids. One-hundred and thirty-four diazotrophs carrying plasmids that ranged in size from 2 to >100 kbp were identified. The majority of the marine bacteria contained one plasmid. Diazotrophs from the short-form S. alterniflora rhizoplane contained significantly fewer plasmids relative to isolates from tall-form S. alterniflora, J. roemarianus and S. patens. Although some plasmids exhibited homology to a nifH gene probe, the majority of the plasmids were classified as cryptic. Two oligonucleotide primers were developed to facilitate genotypic typing of the endogenously isolated marine plasmids by the randomly amplified polymorphic DNA (RAPD)-PCR technique. These primers proved to be more effective than 21 commercially available primers tested to generate RAPD-PCR patterns. Analysis of the RAPD-PCR patterns indicated as many as 71 different plasmid genotypes occurring in diazotroph bacterial assemblages within and between the four different salt marsh grass rhizoplane habitats investigated in this study.

Genetic diversity of indigenous Rhizobium leguminosarum bv. viciae isolates nodulating two different host plants during soil restoration with alfalfa

A total of 360 Rhizobium leguminosarum bv. viciae strains was isolated from three brown-coal mining restoration fields of different age and plant cover (without and in the first and second year of alfalfa, Medicago sativa, cultivation) using two host species (Vicia hirsuta and Pisum sativum) as capture plants. The strains were genetically typed by restriction fragment length polymorphism analysis of polymerase chain reaction (PCR)-generated 16S-23S ribosomal DNA intergenic spacer regions (IGS-RFLP) and characterized by plasmid profiles and RFLP analysis of amplified nodABC genes. The R. leguminosarum bv. viciae population was dominated by the same group of strains (irrespective of the trap plant used). According to type richness, the genetic diversity of indigenous R. leguminosarum in the second year of restoration was lower than in the first year and it resembled that of the fallow field, except for plasmid types, in which it was higher than that of the fallow field. Some of the less frequent nodABC genotypes were associated with distinct chromosomal IGS genotypes and symbiotic plasmids (pSyms) of different sizes, indicating that horizontal transfer and rearrangements of pSym can occur in natural environments. However, the dominant pSym and chromosomal genotypes were strictly correlated suggesting a genetically stable persistence of the prevailing R. leguminosarum bv. viciae genotypes in the absence of its host plant.