Commonalities and differences among symbiosis islands of three Mesorhizobium loti strains

Identification and functional analysis of recent IS transposition events in rhizobia

Rhizobia are alpha- and beta- Proteobacteria that, through the establishment of symbiotic interactions with leguminous plants, are able to fix atmospheric nitrogen as ammonium. The successful establishment of a symbiotic interaction is highly dependent on the availability of nitrogen sources in the soil, and on the specific rhizobia strain. Insertion sequences (ISs) are simple transposable genetic elements that can move to different locations within the host genome and are known to play an important evolutionary role, contributing to genome plasticity by acting as recombination hot-spots, and disrupting coding and regulatory sequences. Disruption of coding sequences may have occurred either in a common ancestor of the species or more recently. By means of ISComapare, we identified Differentially Located ISs (DLISs) in nearly related rhizobial strains of the genera Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium. Our results revealed that recent IS transposition could have a role in adaptation by enabling the activation and inactivation of genes that could dynamically affect the competition and survival of rhizobia in the rhizosphere.

Rhizobial diversity is associated with inoculation history at a two-continent scale

A total of 120 Mesorhizobium strains collected from the central dry zone of Myanmar were analyzed in a pot experiment to evaluate nodulation and symbiotic effectiveness (SE%) in chickpea plants. Phylogenetic analyses revealed all strains belonged to the genus Mesorhizobium according to 16–23S rDNA IGS and the majority of chickpea nodulating rhizobia in Myanmar soils were most closely related to M. gobiense, M. muleiense, M. silamurunense, M. tamadayense and M. temperatum. Around two-thirds of the Myanmar strains (68%) were most closely related to Indian strain IC-2058 (CA-181), which is also most closely related to M. gobiense. There were no strains that were closely related to the cognate rhizobial species to nodulate chickpea: M. ciceri and M. mediterraneum. Strains with diverse 16S–23S rDNA IGS shared similar nodC and nifH gene sequences with chickpea symbionts. Detailed sequence analysis of nodC and nifH found that the strains in Myanmar were somewhat divergent from the group including M. ciceri and were more closely related to M. muleiense and IC-2058. A cross-continent analysis between strains isolated in Australia compared with Myanmar found that there was little overlap in species, where Australian soils were dominated with M. ciceri, M. temperatum and M. huakuii. The only co-occurring species found in both Myanmar and Australia were M. tamadayense and M. silumurunense. Continued inoculation with CC1192 may have reduced diversity of chickpea strains in Australian soils. Isolated strains in Australian and Myanmar had similar adaptive traits, which in some cases were also phylogenetically related. The genetic discrepancy between chickpea nodulating strains in Australia and Myanmar is not only due to inoculation history but to adaptation to soil conditions and crop management over a long period, and there has been virtually no loss of symbiotic efficiency over this time in strains isolated from soils in Myanmar.

Evolution of rhizobial symbiosis islands through insertion sequence-mediated deletion and duplication

Symbiosis between organisms influences their evolution via adaptive changes in genome architectures. Immunity of soybean carrying the Rj2 allele is triggered by NopP (type III secretion system [T3SS]-dependent effector), encoded by symbiosis island A (SymA) in B. diazoefficiens USDA122. This immunity was overcome by many mutants with large SymA deletions that encompassed T3SS (rhc) and N2 fixation (nif) genes and were bounded by insertion sequence (IS) copies in direct orientation, indicating homologous recombination between ISs. Similar deletion events were observed in B. diazoefficiens USDA110 and B. japonicum J5. When we cultured a USDA122 strain with a marker gene sacB inserted into the rhc gene cluster, most sucrose-resistant mutants had deletions in nif/rhc gene clusters, similar to the mutants above. Some deletion mutants were unique to the sacB system and showed lower competitive nodulation capability, indicating that IS-mediated deletions occurred during free-living growth and the host plants selected the mutants. Among 63 natural bradyrhizobial isolates, 2 possessed long duplications (261-357 kb) harboring nif/rhc gene clusters between IS copies in direct orientation via homologous recombination. Therefore, the structures of symbiosis islands are in a state of flux via IS-mediated duplications and deletions during rhizobial saprophytic growth, and host plants select mutualistic variants from the resultant pools of rhizobial populations. Our results demonstrate that homologous recombination between direct IS copies provides a natural mechanism generating deletions and duplications on symbiosis islands.

Symbiosis islands of Loteae-nodulating Mesorhizobium comprise three radiating lineages with concordant nod gene complements and nodulation host-range groupings

Mesorhizobium is a genus of soil bacteria, some isolates of which form an endosymbiotic relationship with diverse legumes of the Loteae tribe. The symbiotic genes of these mesorhizobia are generally carried on integrative and conjugative elements termed symbiosis islands (ICESyms). Mesorhizobium strains that nodulate Lotus spp. have been divided into host-range groupings. Group I (GI) strains nodulate L. corniculatus and L. japonicus ecotype Gifu, while group II (GII) strains have a broader host range, which includes L. pedunculatus. To identify the basis of this extended host range, and better understand Mesorhizobium and ICESym genomics, the genomes of eight Mesorhizobium strains were completed using hybrid long- and short-read assembly. Bioinformatic comparison with previously sequenced mesorhizobia genomes indicated host range was not predicted by Mesorhizobium genospecies but rather by the evolutionary relationship between ICESym symbiotic regions. Three radiating lineages of Loteae ICESyms were identified on this basis, which correlate with Lotus spp. host-range grouping and have lineage-specific nod gene complements. Pangenomic analysis of the completed GI and GII ICESyms identified 155 core genes (on average 30.1 % of a given ICESym). Individual GI or GII ICESyms carried diverse accessory genes with an average of 34.6 % of genes unique to a given ICESym. Identification and comparative analysis of NodD symbiotic regulatory motifs – nod boxes – identified 21 branches across the NodD regulons. Four of these branches were associated with seven genes unique to the five GII ICESyms. The nod boxes preceding the host-range gene nodZ in GI and GII ICESyms were disparate, suggesting regulation of nodZ may differ between GI and GII ICESyms. The broad host-range determinant(s) of GII ICESyms that confer nodulation of L. pedunculatus are likely present amongst the 53 GII-unique genes identified.

Assessment of Polygala paniculata (Polygalaceae) characteristics for evolutionary studies of legume-rhizobia symbiosis

Root nodule (RN) symbiosis is a mutualistic interaction observed between nitrogen-fixing soil bacteria and nodulating plants, which are scattered in only four orders of angiosperms called nitrogen-fixing clade. Most of legumes engage in RN symbiosis with rhizobia. Molecular genetic analyses with legumes and non-leguminous nodulating plants revealed that RN symbiosis utilizes early signalling components that are required for symbiosis with arbuscular mycorrhizal (AM) fungi. However detailed evolutionary processes are still largely unknown. Comparative analyses with non-nodulating species phylogenetically related to legumes could be better strategies to study the evolution of RN symbiosis in legumes. Polygala paniculata is a non-leguminous species that belongs to a family different from legumes but that is classified into the same order, Fabales. It has appropriate characteristics for cultivation in laboratories: small body size, high fertility and short lifecycles. Therefore, we further assessed whether this species is suitable as a model species for comparative studies with legumes. We first validated that the plant we obtained in Palau was truly P. paniculata by molecular phylogenetic analysis using rbcL sequences. The estimated genome size of this species was less than those of two model legumes, Lotus japonicus and Medicago truncatula. We determined conditions for cultivation in vitro and for hairy root formation from P. paniculata seedlings. It would facilitate to investigate gene functions in this species. The ability of P. paniculata to interact with AM fungi was confirmed by inoculation with Rhizophagus irregularis, suggesting the presence of early signalling factors that might be involved in RN symbiosis. Unexpectedly, branching of root hairs was observed when inoculated with Mesorhizobium loti broad host range strain NZP2037, indicating that P. paniculata has the biological potential to respond to rhizobia. We propose that P. paniculata is used as a model plant for the evolutionary study of RN symbiosis.

traG Gene Is Conserved across Mesorhizobium spp. Able to Nodulate the Same Host Plant and Expressed in Response to Root Exudates

Evidences for an involvement of the bacterial type IV secretion system (T4SS) in the symbiotic relationship between rhizobia and legumes have been pointed out by several recent studies. However, information regarding this secretion system in Mesorhizobium is still very scarce. The aim of the present study was to investigate the phylogeny and expression of the traG gene, which encodes a substrate receptor of the T4SS. In addition, the occurrence and genomic context of this and other T4SS genes, namely, genes from tra/trb and virB/virD4 complexes, were also analyzed in order to unveil the structural and functional organization of T4SS in mesorhizobia. The location of the T4SS genes in the symbiotic region of the analyzed rhizobial genomes, along with the traG phylogeny, suggests that T4SS genes could be horizontally transferred together with the symbiosis genes. Regarding the T4SS structural organization in Mesorhizobium, the virB/virD4 genes were absent in all chickpea (Cicer arietinum L.) microsymbionts and in the Lotus symbiont Mesorhizobium japonicum MAFF303099^T. Interestingly, the presence of genes belonging to another secretion system (T3SS) was restricted to these strains lacking the virB/virD4 genes. The traG gene expression was detected in M. mediterraneum Ca36^T and M. ciceri LMS-1 strains when exposed to chickpea root exudates and also in the early nodules formed by M. mediterraneum Ca36^T, but not in older nodules. This study contributes to a better understanding of the importance of T4SS in mutualistic symbiotic bacteria.

Potentially Mobile Denitrification Genes Identified in Azospirillum sp. Strain TSH58

Denitrification ability is sporadically distributed among diverse bacteria, archaea, and fungi. In addition, disagreement has been found between denitrification gene phylogenies and the 16S rRNA gene phylogeny. These facts have suggested potential occurrences of horizontal gene transfer (HGT) for the denitrification genes. However, evidence of HGT has not been clearly presented thus far. In this study, we identified the sequences and the localization of the nitrite reductase genes in the genomes of 41 denitrifying Azospirillum sp. strains and searched for mobile genetic elements that contain denitrification genes. All Azospirillum sp. strains examined in this study possessed multiple replicons (4 to 11 replicons), with their sizes ranging from 7 to 1,031 kbp. Among those, the nitrite reductase gene nirK was located on large replicons (549 to 941 kbp). Genome sequencing showed that Azospirillum strains that had similar nirK sequences also shared similar nir-nor gene arrangements, especially between the TSH58, Sp7T, and Sp245 strains. In addition to the high similarity between nir-nor gene clusters among the three Azospirillum strains, a composite transposon structure was identified in the genome of strain TSH58, which contains the nir-nor gene cluster and the novel IS6 family insertion sequences (ISAz581 and ISAz582). The nirK gene within the composite transposon system was actively transcribed under denitrification-inducing conditions. Although not experimentally verified in this study, the composite transposon system containing the nir-nor gene cluster could be transferred to other cells if it is moved to a prophage region and the phage becomes activated and released outside the cells. Taken together, strain TSH58 most likely acquired its denitrification ability by HGT from closely related Azospirillum sp. denitrifiers.IMPORTANCE The evolutionary history of denitrification is complex. While the occurrence of horizontal gene transfer has been suggested for denitrification genes, most studies report circumstantial evidences, such as disagreement between denitrification gene phylogenies and the 16S rRNA gene phylogeny. Based on the comparative genome analyses of Azospirillum sp. denitrifiers, we identified denitrification genes, including nirK and norCBQD, located on a mobile genetic element in the genome of Azospirillum sp. strain TSH58. The nirK was actively transcribed under denitrification-inducing conditions. Since this gene was the sole nitrite reductase gene in strain TSH58, this strain most likely benefitted by acquiring denitrification genes via horizontal gene transfer. This finding will significantly advance our scientific knowledge regarding the ecology and evolution of denitrification.

The Rhizobia-Lotus Symbioses: Deeply Specific and Widely Diverse

The symbiosis between Lotus and rhizobia has been long considered very specific and only two bacterial species were recognized as the microsymbionts of Lotus: Mesorhizobium loti was considered the typical rhizobia for the L. corniculatus complex, whereas Bradyrhizobium sp. (Lotus) was the symbiont for L. uliginosus and related species. As discussed in this review, this situation has dramatically changed during the last 15 years, with the characterization of nodule bacteria from worldwide geographical locations and from previously unexplored Lotus spp. Current data support that the Lotus rhizobia are dispersed amongst nearly 20 species in five genera (Mesorhizobium, Bradyrhizobium, Rhizobium, Ensifer, and Aminobacter). As a consequence, M. loti could be regarded an infrequent symbiont of Lotus, and several plant-bacteria compatibility groups can be envisaged. Despite the great progress achieved with the model L. japonicus in understanding the establishment and functionality of the symbiosis, the genetic and biochemical bases governing the stringent host-bacteria compatibility pairships within the genus Lotus await to be uncovered. Several Lotus spp. are grown for forage, and inoculation with rhizobia is a common practice in various countries. However, the great diversity of the Lotus rhizobia is likely squandered, as only few bacterial strains are used as inoculants for Lotus pastures in very different geographical locations, with a great variety of edaphic and climatic conditions. The agroecological potential of the genus Lotus can not be fully harnessed without acknowledging the great diversity of rhizobia-Lotus interactions, along with a better understanding of the specific plant and bacterial requirements for optimal symbiotic nitrogen fixation under increasingly constrained environmental conditions.

Sequential induction of three recombination directionality factors directs assembly of tripartite integrative and conjugative elements

Tripartite integrative and conjugative elements (ICE³) are a novel form of ICE that exist as three separate DNA regions integrated within the genomes of Mesorhizobium spp. Prior to conjugative transfer the three ICE³ regions of M. ciceri WSM1271 ICEMcSym¹²⁷¹ combine and excise to form a single circular element. This assembly requires three coordinated recombination events involving three site-specific recombinases IntS, IntG and IntM. Here, we demonstrate that three excisionases–or recombination directionality factors—RdfS, RdfG and RdfM are required for ICE³ excision. Transcriptome sequencing revealed that expression of ICE³ transfer and conjugation genes was induced by quorum sensing. Quorum sensing activated expression of rdfS, and in turn RdfS stimulated transcription of both rdfG and rdfM. Therefore, RdfS acts as a “master controller” of ICE³ assembly and excision. The dependence of all three excisive reactions on RdfS ensures that ICE³ excision occurs via a stepwise sequence of recombination events that avoids splitting the chromosome into a non-viable configuration. These discoveries expose a surprisingly simple control system guiding molecular assembly of these novel and complex mobile genetic elements and highlight the diverse and critical functions of excisionase proteins in control of horizontal gene transfer.

Bacteria evolve and adapt quickly through the horizontal transfer of DNA. A major mechanism facilitating this transfer is conjugation. Conjugative DNA elements that integrate into the chromosome are termed ‘Integrative and Conjugative Elements’ (ICE). We recently discovered a unique form of ICE that undergoes a complex series of recombination events with the host chromosome to split itself into three separate parts. This tripartite ICE must also precisely order its recombination when leaving the current host to avoid splitting the host chromosome and the ICE into non-viable parts. In this work, we show that the tripartite ICEs use chemical cell-cell communication to stimulate recombination and that recombination events are specifically ordered through cascaded transcriptional activation of small DNA-binding proteins called recombination directionality factors. Despite the inherent complexity of tripartite ICEs this work exposes a surprisingly simple system to stimulate their precise and ordered molecular assembly prior to horizontal transfer.

Symbiosis island shuffling with abundant insertion sequences in the genomes of extra-slow-growing strains of soybean bradyrhizobia

Extra-slow-growing bradyrhizobia from root nodules of field-grown soybeans harbor abundant insertion sequences (ISs) and are termed highly reiterated sequence-possessing (HRS) strains. We analyzed the genome organization of HRS strains with the focus on IS distribution and symbiosis island structure. Using pulsed-field gel electrophoresis, we consistently detected several plasmids (0.07 to 0.4 Mb) in the HRS strains (NK5, NK6, USDA135, 2281, USDA123, and T2), whereas no plasmids were detected in the non-HRS strain USDA110. The chromosomes of the six HRS strains (9.7 to 10.7 Mb) were larger than that of USDA110 (9.1 Mb). Using MiSeq sequences of 6 HRS and 17 non-HRS strains mapped to the USDA110 genome, we found that the copy numbers of ISRj1, ISRj2, ISFK1, IS1632, ISB27, ISBj8, and IS1631 were markedly higher in HRS strains. Whole-genome sequencing showed that the HRS strain NK6 had four small plasmids (136 to 212 kb) and a large chromosome (9,780 kb). Strong colinearity was found between 7.4-Mb core regions of the NK6 and USDA110 chromosomes. USDA110 symbiosis islands corresponded mainly to five small regions (S1 to S5) within two variable regions, V1 (0.8 Mb) and V2 (1.6 Mb), of the NK6 chromosome. The USDA110 nif gene cluster (nifDKENXSBZHQW-fixBCX) was split into two regions, S2 and S3, where ISRj1-mediated rearrangement occurred between nifS and nifB. ISs were also scattered in NK6 core regions, and ISRj1 insertion often disrupted some genes important for survival and environmental responses. These results suggest that HRS strains of soybean bradyrhizobia were subjected to IS-mediated symbiosis island shuffling and core genome degradation.

Genome sequence of the Lotus spp. microsymbiont Mesorhizobium loti strain NZP2037

Mesorhizobium loti strain NZP2037 was isolated in 1961 in Palmerston North, New Zealand from a Lotus divaricatus root nodule. Compared to most other M. loti strains, it has a broad host range and is one of very few M. loti strains able to form effective nodules on the agriculturally important legume Lotus pedunculatus. NZP2037 is an aerobic, Gram negative, non-spore-forming rod. This report reveals that the genome of M. loti strain NZP2037 does not harbor any plasmids and contains a single scaffold of size 7,462,792 bp which encodes 7,318 protein-coding genes and 70 RNA-only encoding genes. This rhizobial genome is one of 100 sequenced as part of the DOE Joint Genome Institute 2010 Genomic Encyclopedia for Bacteria and Archaea-Root Nodule Bacteria (GEBA-RNB) project.