Plasmids impact on rhizobia-legumes symbiosis in diverse environments

Comparative genome analysis of Sesbania cannabina-nodulating Rhizobium spp. revealing the symbiotic and transferrable characteristics of symbiosis plasmids

Symbiotic nitrogen fixation between legumes and rhizobia makes a great contribution to the terrestrial ecosystem. The successful symbiosis between the partners mainly depends on the nod and nif genes in rhizobia, while the specific symbiosis is mainly determined by the structure of Nod factors and the corresponding secretion systems (type III secretion system; T3SS), etc. These symbiosis genes are usually located on symbiotic plasmids or a chromosomal symbiotic island, both could be transferred interspecies. In our previous studies, Sesbania cannabina-nodulating rhizobia across the world were classified into 16 species of four genera and all the strains, especially those of Rhizobium spp., harboured extraordinarily highly conserved symbiosis genes, suggesting that horizontal transfer of symbiosis genes might have happened among them. In order to learn the genomic basis of diversification of rhizobia under the selection of host specificity, we performed this study to compare the complete genome sequences of four Rhizobium strains associated with S. cannabina, YTUBH007, YTUZZ027, YTUHZ044 and YTUHZ045. Their complete genomes were sequenced and assembled at the replicon level. Each strain represents a different species according to the average nucleotide identity (ANI) values calculated using the whole-genome sequences; furthermore, except for YTUBH007, which was classified as Rhizobium binae, the remaining three strains were identified as new candidate species. A single symbiotic plasmid sized 345-402 kb containing complete nod, nif, fix, T3SS and conjugal transfer genes was detected in each strain. The high ANI and amino acid identity (AAI) values, as well as the close phylogenetic relationships among the entire symbiotic plasmid sequences, indicate that they have the same origin and the entire plasmid has been transferred among different Rhizobium species. These results indicate that S. cannabina stringently selects a certain symbiosis gene background of the rhizobia for nodulation, which might have forced the symbiosis genes to transfer from some introduced rhizobia to the related native or local-condition-adapted bacteria. The existence of almost complete conjugal transfer related elements, but not the gene virD, indicated that the self-transfer of the symbiotic plasmid in these rhizobial strains may be realized via a virD-independent pathway or through another unidentified gene. This study provides insight for the better understanding of high-frequency symbiotic plasmid transfer, host-specific nodulation and the host shift for rhizobia.

Nodules of wild legumes as unique natural hotspots of antibiotic resistance genes

Root nodules (RN) of legumes have distinct microenvironment from their symbiotic roots and surrounding soils. The rhizobia can withstand the host-produced phytoalexins and antimicrobial compounds. We thus hypothesize that the wild legume RN may develop unique natural resistome and be antibiotic resistance gene (ARG) hotspots. In this study, in comparison with rhizosphere soil (RS) and bulk soil (BS), we characterized the feature of antibiotic resistance in the RN of two wild legumes, Medicago polymorpha and Astragalus sinicus, by metagenomics. It was shown that the total relative abundance of ARGs followed the order of RN > RS > BS for both legumes. ARGs encoding antibiotic efflux pump predominated in all samples with increased proportion from BS to RN samples for both legumes. Totally 275 ARG subtypes were detected, and diversity of ARGs in RN was significantly lower than in BS samples for both legumes. 32 and 25 unique ARGs subtypes were detected in RN of both legumes. Bacterial community played a key role in shaping nodule-associated resistome because both ARG profiles and bacterial community differed greatly among BS, RS and RN. Rhizobia potentially hosted 10 and 15 ARGs subtypes for both legumes. The number and proportion of plasmid- and ARG-carrying contigs (ACCs) were higher in RN than in BS. Host tracking analysis of plasmid-ACCs suggests that proportion of rhizobial bacteria identified as their hosts decreased from BS to RN samples. No plasmid-ACCs with multiple ARGs were observed in BS samples, whereas they were detected in RN samples of both legumes. Our study showed that even wild legume nodules are unique natural ARG hotspots and enough attention should be paid to the dissemination risk of ARGs posed by globally produced legume crops.

Symbiosis Contribution of Non-nodulating <i>Bradyrhizobium cosmicum</i> S23321 after Transferal of the Symbiotic Plasmid pDOA9

The symbiotic properties of rhizobial bacteria are driven by the horizontal gene transfer of symbiotic genes, which are located in symbiosis islands or on plasmids. The symbiotic megaplasmid pDOA9 of Bradyrhizobium sp. DOA9, carrying the nod, nif, fix, and type three secretion system (T3SS) genes, has been conjugatively transferred to different Bradyrhizobium strains. In the present study, non-nodulating B. cosmicum S23321, which shows a close phylogenetic relationship with Bradyrhizobium sp. DOA9, but lacks symbiotic properties, was used to carry pDOA9 (annotated as chimeric S2:pDOA9). The results obtained showed that pDOA9 conferred symbiotic properties on S23321; however, nodulation phenotypes varied among the DOA9, chimeric ORS278:pDOA9, and S2:pDOA9 strains even though they all carried symbiotic pDOA9 plasmid. S23321 appeared to gain symbiotic nodulation from pDOA9 by processing nodulation genes and broadening the host range. The present results also showed the successful formation of active nodules in Arachis hypogaea (Dalbergoid) and Vigna radiata (Millitoid) by chimeric S2:pDOA9, while Crotalaria juncea (Genistoid) and Macroptilium atropurpureum (Millitoid) formed nodule-like structures. The formation of nodules and nodule-like structures occurred in a nod factor-dependent manner because the nod factor-lacking strain (S2:pDOA9ΩnodB) completely abolished nodulation in all legumes tested. Moreover, T3SS carried by S2:pDOA9 exerted negative effects on symbiosis with Crotalaria juncea, which was consistent with the results obtained on DOA9. T3SS exhibited symbiotic compatibility with V. radiata when nodulated by S23321. These outcomes implied that pDOA9 underwent changes during legume evolution that broadened host specificity and the compatibility of nodulation in a manner that was dependent on the chromosomal background of the recipient as well as legume host restrictions.

Phylogeographic distribution of rhizobia nodulating common bean (Phaseolus vulgaris L.) in Ethiopia

Rhizobia are soilborne bacteria that form symbiotic relations with legumes and fix atmospheric nitrogen. The nitrogen fixation potential depends on several factors such as the type of host and symbionts and on environmental factors that affect the distribution of rhizobia. We isolated bacteria nodulating common bean in Southern Ethiopia to evaluate their genetic diversity and phylogeography at nucleotide, locus (gene/haplotype) and species levels of genetic hierarchy. Phylogenetically, eight rhizobial genospecies (including previous collections) were determined that had less genetic diversity than found among reference strains. The limited genetic diversity of the Ethiopian collections was due to absence of many of the Rhizobium lineages known to nodulate beans. Rhizobium etli and Rhizobiumphaseoli were predominant strains of bean-nodulating rhizobia in Ethiopia. We found no evidence for a phylogeographic pattern in strain distribution. However, joint analysis of the current and previous collections revealed differences between the two collections at nucleotide level of genetic hierarchy. The differences were due to genospecies Rhizobium aethiopicum that was only isolated in the earlier collection.

Physiological disruption, structural deformation and low grain yield induced by neonicotinoid insecticides in chickpea: A long term phytotoxicity investigation

Arbitrary use of insecticides in agricultural practices cause severe environmental hazard that adversely affects the growth and productivity of edible crops. Considering theses, the aim of the present study was to evaluate the toxicological potential of two neonicotinoid insecticides, imidacloprid (IMID) and thiamethoxam (THIA) using chickpea as a test crop. Application of insecticides at three varying doses revealed a gradual decrease in biological performance of chickpea plants which however, varied noticeably among insecticides and their doses. Significant (P ≤ 0.05) decline in germination efficiency, length of plant organs under in vitro condition was observed in a dose related manner. Among insecticides, 300 μgIMIDkg^-1 (3X dose) maximally and significantly (P ≤ 0.05) inhibited germination efficiency, vigor index, length, dry matter accumulation, photosynthetic pigment formation, nodule volume and mass, nutrient uptake, grain yield and protein over untreated control. In contrast, 75 μgTHIAkg^-1 (3X dose) considerably declined the leghaemoglobin content, shoot phosphorus and root nitrogen. Enhanced expression of stress biomarkers including proline, malondialdehyde (MDA), and antioxidant defence enzymes was noticed in the presence of insecticides. For instance, at 3X IMID, shoot proline, MDA, ascorbate peroxidase (APX), guaiacol peroxidase (GPX) and peroxidase (POD) were increased significantly (P ≤ 0.05) by 66%, 81%, 36% and 35%, respectively. Additionally, electrolyte leakage was maximally (77%) increased at 3X dose of IMID, whereas, H₂O₂ in foliage was maximally accumulated (0.0156 μ moles min^-1 g^-1 fw) at 3X dose of THIA which was 58% greater than untreated foliage. A clear distortion/damage in tip and surface of roots and ultrastructural deformation in xylem and phloem of plant tissues as indication of insecticidal phytotoxicity was observed under scanning electron microscope (SEM). For oxidative stress and cytotoxicity assessment, root tips were stained with a combination of acridine orange and propidium iodide, and Evan blue dyes and examined. Confocal laser scanning microscopic (CLSM) images of roots revealed a 10-fold and 13.5-fold increase in red and blue fluorescence when 3X IMID treated roots were assessed quantitatively. Conclusively, the present investigation recommends that a careful and protected approach should be adopted before the application of insecticides in agricultural ecosystems.

Plasmids Related to the Symbiotic Nitrogen Fixation Are Not Only Cooperated Functionally but Also May Have Evolved over a Time Span in Family Rhizobiaceae

Rhizobia are soil bacteria capable of forming symbiotic nitrogen-fixing nodules associated with leguminous plants. In fast-growing legume-nodulating rhizobia, such as the species in the family Rhizobiaceae, the symbiotic plasmid is the main genetic basis for nitrogen-fixing symbiosis, and is susceptible to horizontal gene transfer. To further understand the symbioses evolution in Rhizobiaceae, we analyzed the pan-genome of this family based on 92 genomes of type/reference strains and reconstructed its phylogeny using a phylogenomics approach. Intriguingly, although the genetic expansion that occurred in chromosomal regions was the main reason for the high proportion of low-frequency flexible gene families in the pan-genome, gene gain events associated with accessory plasmids introduced more genes into the genomes of nitrogen-fixing species. For symbiotic plasmids, although horizontal gene transfer frequently occurred, transfer may be impeded by, such as, the host’s physical isolation and soil conditions, even among phylogenetically close species. During coevolution with leguminous hosts, the plasmid system, including accessory and symbiotic plasmids, may have evolved over a time span, and provided rhizobial species with the ability to adapt to various environmental conditions and helped them achieve nitrogen fixation. These findings provide new insights into the phylogeny of Rhizobiaceae and advance our understanding of the evolution of symbiotic nitrogen fixation.

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.

Phylogenomic Rhizobium Species Are Structured by a Continuum of Diversity and Genomic Clusters

The bacterial genus Rhizobium comprises diverse symbiotic nitrogen-fixing species associated with the roots of plants in the Leguminosae family. Multiple genomic clusters defined by whole genome comparisons occur within Rhizobium, but their equivalence to species is controversial. In this study we investigated such genomic clusters to ascertain their significance in a species phylogeny context. Phylogenomic inferences based on complete sets of ribosomal proteins and stringent core genome markers revealed the main lineages of Rhizobium. The clades corresponding to R. etli and R. leguminosarum species show several genomic clusters with average genomic nucleotide identities (ANI > 95%), and a continuum of divergent strains, respectively. They were found to be inversely correlated with the genetic distance estimated from concatenated ribosomal proteins. We uncovered evidence of a Rhizobium pangenome that was greatly expanded, both in its chromosomes and plasmids. Despite the variability of extra-chromosomal elements, our genomic comparisons revealed only a few chromid and plasmid families. The presence/absence profile of genes in the complete Rhizobium genomes agreed with the phylogenomic pattern of species divergence. Symbiotic genes were distributed according to the principal phylogenomic Rhizobium clades but did not resolve genome clusters within the clades. We distinguished some types of symbiotic plasmids within Rhizobium that displayed different rates of synonymous nucleotide substitutions in comparison to chromosomal genes. Symbiotic plasmids may have been repeatedly transferred horizontally between strains and species, in the process displacing and substituting pre-existing symbiotic plasmids. In summary, the results indicate that Rhizobium genomic clusters, as defined by whole genomic identities, might be part of a continuous process of evolutionary divergence that includes the core and the extrachromosomal elements leading to species formation.

Comparative analysis of rhizobial chromosomes and plasmids to estimate their evolutionary relationships

In the present study, complete genomic sequences retrieved from 57 rhizobial strains that covered four genera including 11 species were analyzed comprehensively. The four types of replicons: chromosomes, chromids, nonsymbiotic plasmids, and symbiotic plasmids were investigated and compared among these strains. Results showed that co-evolution occurred among these four replicons based on the similarities in average nucleotide identity. High correlation coefficient r values were observed between chromosomes and chromids, as well as between chromosomes and nonsymbiotic plasmids. Chromosomes and symbiotic plasmids showed different phylogenetic topology based on their core genes. Population structure analyses were performed to extrapolate the evolutionary histories of the test strains based on their chromosomal and symbiotic plasmid background. This resulted in seven ancestral types for chromosomal genes and three ancestral types for symbiotic plasmid genes. Rhizobial strains containing chromosome genes with ancestral type E tend to contain symbiotic plasmid genes with ancestral type II, while rhizobial strains containing chromosome genes with ancestral type G tend to contain symbiotic plasmid genes with ancestral type III. Seventeen strains associated with different host plant species which harbored the symbiotic genes with ancestral type I, exhibited high genetic diversity. In addition, Fu's test of the symbiotic plasmid genes with ancestral type III had undergone an expansion event, implying the influence of negative selection on these symbiotic plasmid genes.

Co-existence of Rhizobia and Diverse Non-rhizobial Bacteria in the Rhizosphere and Nodules of Dalbergia odorifera Seedlings Inoculated with Bradyrhizobium elkanii, Rhizobium multihospitium-Like and Burkholderia pyrrocinia-Like Strains

Rhizobia induce root nodules and fix atmospheric N₂ for most legume species in exchange for carbon. However, the diverse endophytic non-rhizobial bacteria in legume nodules that co-exist with rhizobia are often ignored because they are difficult to cultivate using routine cultivation approaches. To enhance our understanding of the incidence and diversity of legume–bacteria associations, a high-throughput sequencing analysis of bacterial 16S rRNA genes was used to examine the bacterial community in the rhizospheres and root nodules of Dalbergia odorifera seedlings that were uninoculated or inoculated with Bradyrhizobium elkanii H255, Rhizobium multihospitium–like HT221, or Burkholderia pyrrocinia–like H022238, in two growth media (nitrogen [N]-supplied soil or N-omitted potting mix). Seedlings inoculated with Bradyrhizobium had significantly more nodules than seedlings in the other inoculation conditions, regardless of growth media. Using the ¹⁵N natural abundance method, it was shown that the inoculated plants had significantly higher N₂ fixation efficiency (48–57%) and specific nodule activity [269–313 μg N mg⁻¹ of dry weight (dwt) nodule] compared to the uninoculated plants (203 μg N mg⁻¹ dwt nodule). The 16S rRNA gene analysis showed that there was generally a higher bacterial diversity in the rhizosphere than in the nodules in the corresponding condition. Both rhizobial inoculation and media status significantly altered the bacterial communities in the rhizospheres and nodules (P < 0.05), with the exception of the inoculated soil rhizospheres. Regarding non-rhizobial bacteria, three genera, i.e., Lactococcus, Bacillus, and Pseudomonas, were consistently enriched in the rhizosphere and Bradyrhizobium, Chloroplast norank (which belongs to Cyanobacteria), and Lactococcus were commonly found in the nodules. In contrast, common rhizobial genera (including Rhizobium, Mesorhizobium, and Burkholderia) were only present in the nodules at low relative abundances (0.01–3.41%). Regarding non-rhizobial bacteria, 32 genera were found in the nodules, with non-rhizobial bacteria being predominant in the N omitted potting mix (with a relative abundance of 56–87%). This study suggests that legume nodules are inhabited by a high diversity of non-rhizobial bacteria, which may play a vital role in nodulation and N₂ fixation in the host plants.

Retrieved 16S rRNA and nifH sequences reveal co-dominance of Bradyrhizobium and Ensifer (Sinorhizobium) strains in field-collected root nodules of the promiscuous host Vigna radiata (L.) R. Wilczek

In the present study, the relative distribution of endophytic rhizobia in field-collected root nodules of the promiscuous host mung bean was investigated by sequencing of 16S ribosomal RNA (rRNA) and nifH genes, amplified directly from the nodule DNA. Co-dominance of the genera Bradyrhizobium and Ensifer was indicated by 32.05 and 35.84% of the total retrieved 16S rRNA sequences, respectively, and the sequences of genera Mesorhizobium and Rhizobium comprised only 0.06 and 2.06% of the recovered sequences, respectively. Sequences amplified from rhizosphere soil DNA indicated that only a minor fraction originated from Bradyrhizobium and Ensifer strains, comprising about 0.46 and 0.67% of the total retrieved sequences, respectively. 16S rRNA gene sequencing has also identified the presence of several non-rhizobial endophytes from phyla Proteobacteria, Actinobacteria, Bacteroides, and Firmicutes. The nifH sequences obtained from nodules also confirmed the co-dominance of Bradyrhizobium (39.21%) and Ensifer (59.23%) strains. The nifH sequences of the genus Rhizobium were absent, and those of genus Mesorhizobium comprised only a minor fraction of the sequences recovered from the nodules and rhizosphere soil samples. Two bacterial isolates, identified by 16S rRNA gene sequence analysis as Bradyrhizobium strain Vr51 and Ensifer strain Vr38, successfully nodulated the original host (mung bean) plants. Co-dominance of Bradyrhizobium and Ensifer strains in the nodules of mung bean indicates the potential role of the host plant in selecting specific endophytic rhizobial populations. Furthermore, successful nodulation of mung bean by the isolates showed that strains of both the genera Bradyrhizobium and Ensifer can be used for production of inoculum.