Experimental Evolution of Legume Symbionts: What Have We Learnt?

Resistance to carbapenems in the urban soil isolate Cupriavidus taiwanensis S2-1-W is associated with OXA-1206, a newly discovered carbapenemase

AIMS

Cupriavidus isolates are found in environmental and clinical samples and are often resistant to carbapenems, which are last-resort antibiotics. However, their carbapenem-resistance molecular mechanisms remain unknown. This study aimed to (i) characterize and sequence the carbapenem-resistant soil isolate Cupriavidus taiwanensis S2-1-W to uncover its antibiotic resistance determinants; and (ii) clone and characterize a putative novel carbapenemase gene identified in this isolate.

METHODS AND RESULTS

Antibiotic susceptibility testing of C. taiwanensis S2-1-W revealed that it was resistant to most carbapenems, other β-lactams, and aminoglycosides tested. Genome sequencing of this isolate revealed a complex chromosomal resistome that included multidrug efflux pump genes, one aminoglycoside transferase gene, and three β-lactamase genes. Among them, we identified a novel putative class D β-lactamase gene (blaOXA-1206) that is highly conserved among other sequenced C. taiwanensis isolates. Cloning and characterization of blaOXA-1206 confirmed that it encodes for a newly discovered carbapenemase (OXA-1206) that confers resistance to carbapenems and other β-lactams.

CONCLUSION

Carbapenem-resistance in C. taiwanensis S2-1-W is associated with a newly discovered carbapenemase, OXA-1206.

The role of microbial interactions on rhizobial fitness

Rhizobia are soil bacteria that can establish a nitrogen-fixing symbiosis with legume plants. As horizontally transmitted symbionts, the life cycle of rhizobia includes a free-living phase in the soil and a plant-associated symbiotic phase. Throughout this life cycle, rhizobia are exposed to a myriad of other microorganisms that interact with them, modulating their fitness and symbiotic performance. In this review, we describe the diversity of interactions between rhizobia and other microorganisms that can occur in the rhizosphere, during the initiation of nodulation, and within nodules. Some of these rhizobia-microbe interactions are indirect, and occur when the presence of some microbes modifies plant physiology in a way that feeds back on rhizobial fitness. We further describe how these interactions can impose significant selective pressures on rhizobia and modify their evolutionary trajectories. More extensive investigations on the eco-evolutionary dynamics of rhizobia in complex biotic environments will likely reveal fascinating new aspects of this well-studied symbiotic interaction and provide critical knowledge for future agronomical applications.

Comparative phylotranscriptomics reveals ancestral and derived root nodule symbiosis programmes

Symbiotic interactions such as the nitrogen-fixing root nodule symbiosis (RNS) have structured ecosystems during the evolution of life. Here we aimed at reconstructing ancestral and intermediate steps that shaped RNS observed in extant flowering plants. We compared the symbiotic transcriptomic responses of nine host plants, including the mimosoid legume Mimosa pudica for which we assembled a chromosome-level genome. We reconstructed the ancestral RNS transcriptome composed of most known symbiotic genes together with hundreds of novel candidates. Cross-referencing with transcriptomic data in response to experimentally evolved bacterial strains with gradual symbiotic proficiencies, we found the response to bacterial signals, nodule infection, nodule organogenesis and nitrogen fixation to be ancestral. By contrast, the release of symbiosomes was associated with recently evolved genes encoding small proteins in each lineage. We demonstrate that the symbiotic response was mostly in place in the most recent common ancestor of the RNS-forming species more than 90 million years ago.

Evaluation of qPCR to Detect Shifts in Population Composition of the Rhizobial Symbiont Mesorhizobium japonicum during Serial in Planta Transfers

Microbial symbionts range from mutualistic to commensal to antagonistic. While these roles are distinct in their outcome, they are also fluid in a changing environment. Here, we used the Lotus japonicus-Mesorhizobium japonicum symbiosis to investigate short-term and long-term shifts in population abundance using an effective, fast, and low-cost tracking methodology for M. japonicum. We use quantitative polymerase chain reaction (qPCR) to track previously generated signature-tagged M. japonicum mutants targeting the Tn5 transposon insertion and the flanking gene. We used a highly beneficial wild type and moderately beneficial and non-beneficial mutants of M. japonicum sp. nov. to demonstrate the specificity of these primers to estimate the relative abundance of each genotype within individual nodules and after serial transfers to new hosts. For the moderate and non-beneficial genotypes, qPCR allowed us to differentiate genotypes that are phenotypically indistinguishable and investigate host control with suboptimal symbionts. We consistently found the wild type increasing in the proportion of the population, but our data suggest a potential reproductive trade-off between the moderate and non-beneficial genotypes. The multi-generation framework we used, coupled with qPCR, can easily be scaled up to track dozens of M. japonicum mutants simultaneously. Moreover, these mutants can be used to explore M. japonicum genotype abundance in the presence of a complex soil community.

Adaptive Evolution of Rhizobial Symbiosis beyond Horizontal Gene Transfer: From Genome Innovation to Regulation Reconstruction

There are ubiquitous variations in symbiotic performance of different rhizobial strains associated with the same legume host in agricultural practices. This is due to polymorphisms of symbiosis genes and/or largely unexplored variations in integration efficiency of symbiotic function. Here, we reviewed cumulative evidence on integration mechanisms of symbiosis genes. Experimental evolution, in concert with reverse genetic studies based on pangenomics, suggests that gain of the same circuit of key symbiosis genes through horizontal gene transfer is necessary but sometimes insufficient for bacteria to establish an effective symbiosis with legumes. An intact genomic background of the recipient may not support the proper expression or functioning of newly acquired key symbiosis genes. Further adaptive evolution, through genome innovation and reconstruction of regulation networks, may confer the recipient of nascent nodulation and nitrogen fixation ability. Other accessory genes, either co-transferred with key symbiosis genes or stochastically transferred, may provide the recipient with additional adaptability in ever-fluctuating host and soil niches. Successful integrations of these accessory genes with the rewired core network, regarding both symbiotic and edaphic fitness, can optimize symbiotic efficiency in various natural and agricultural ecosystems. This progress also sheds light on the development of elite rhizobial inoculants using synthetic biology procedures.

A Novel Module Promotes Horizontal Gene Transfer in Azorhizobium caulinodans ORS571

Azorhizobium caulinodans ORS571 contains an 87.6 kb integrative and conjugative element (ICE^Ac) that conjugatively transfers symbiosis genes to other rhizobia. Many hypothetical redundant gene fragments (rgfs) are abundant in ICE^Ac, but their potential function in horizontal gene transfer (HGT) is unknown. Molecular biological methods were employed to delete hypothetical rgfs, expecting to acquire a minimal ICE^Ac and consider non-functional rgfs as editable regions for inserting genes related to new symbiotic functions. We determined the significance of rgf4 in HGT and identified the physiological function of genes designated rihF1a (AZC_3879), rihF1b (AZC_RS26200), and rihR (AZC_3881). In-frame deletion and complementation assays revealed that rihF1a and rihF1b work as a unit (rihF1) that positively affects HGT frequency. The EMSA assay and lacZ-based reporter system showed that the XRE-family protein RihR is not a regulator of rihF1 but promotes the expression of the integrase (intC) that has been reported to be upregulated by the LysR-family protein, AhaR, through sensing host’s flavonoid. Overall, a conservative module containing rihF1 and rihR was characterized, eliminating the size of ICE^Ac by 18.5%. We propose the feasibility of constructing a minimal ICE^Ac element to facilitate the exchange of new genetic components essential for symbiosis or other metabolic functions between soil bacteria.

Host-Associated Rhizobial Fitness: Dependence on Nitrogen, Density, Community Complexity, and Legume Genotype

The environmental context of the nitrogen-fixing mutualism between leguminous plants and rhizobial bacteria varies over space and time. Variation in resource availability, population density, and composition likely affect the ecology and evolution of rhizobia and their symbiotic interactions with hosts. We examined how host genotype, nitrogen addition, rhizobial density, and community complexity affected selection on 68 rhizobial strains in the Sinorhizobium meliloti-Medicago truncatula mutualism. As expected, host genotype had a substantial effect on the size, number, and strain composition of root nodules (the symbiotic organ). The understudied environmental variable of rhizobial density had a stronger effect on nodule strain frequency than the addition of low nitrogen levels. Higher inoculum density resulted in a nodule community that was less diverse and more beneficial but only in the context of the more selective host genotype. Higher density resulted in more diverse and less beneficial nodule communities with the less selective host. Density effects on strain composition deserve additional scrutiny as they can create feedback between ecological and evolutionary processes. Finally, we found that relative strain rankings were stable across increasing community complexity (2, 3, 8, or 68 strains). This unexpected result suggests that higher-order interactions between strains are rare in the context of nodule formation and development. Our work highlights the importance of examining mechanisms of density-dependent strain fitness and developing theoretical predictions that incorporate density dependence. Furthermore, our results have translational relevance for overcoming establishment barriers in bioinoculants and motivating breeding programs that maintain beneficial plant-microbe interactions across diverse agroecological contexts. IMPORTANCE Legume crops establish beneficial associations with rhizobial bacteria that perform biological nitrogen fixation, providing nitrogen to plants without the economic and greenhouse gas emission costs of chemical nitrogen inputs. Here, we examine the influence of three environmental factors that vary in agricultural fields on strain relative fitness in nodules. In addition to manipulating nitrogen, we also use two biotic variables that have rarely been examined: the rhizobial community's density and complexity. Taken together, our results suggest that (i) breeding legume varieties that select beneficial strains despite environmental variation is possible, (ii) changes in rhizobial population densities that occur routinely in agricultural fields could drive evolutionary changes in rhizobial populations, and (iii) the lack of higher-order interactions between strains will allow the high-throughput assessments of rhizobia winners and losers during plant interactions.

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'.

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.

Response of Sugarcane Rhizosphere Bacterial Community to Drought Stress

Sugarcane is an important sugar and energy crop, and its yield is greatly affected by drought. Although a large number of studies have shown that rhizosphere microorganisms can help improve the adaptability of plants to biotic or abiotic stresses, there is a lack of studies on the adaptability of sugarcane rhizosphere microbial communities to host plants. Therefore, we conducted drought stress treatment and normal irrigation treatment on three sugarcane varieties GT21, GT31, and GT42 widely cultivated in Guangxi. Using 16S rDNA sequencing technology to analyze the changes in abundance of the sugarcane rhizosphere bacterial community under different treatments, combined with the determination of soil enzyme activity, soil nutrient content, and sugarcane physiological characteristics, we explored the sugarcane rhizosphere bacterial community response to drought stress. In addition, we used the structural equation model to verify the response path of sugarcane rhizosphere bacteria. The results show that the bacterial community structure in the rhizosphere of sugarcane is stable under normal water conditions. The change in the bacterial community structure under drought stress has a 25.2% correlation with the drought adaptability of sugarcane, but the correlation with drought stress is as high as 42.17%. The changes in abundance of rhizosphere bacteria under drought stress are mainly concentrated in the phylum Rhizobiales and Streptomycetales. This change is directly related to the physiological state of the host plant under drought stress, soil available phosphorus, soil urease and soil acid protease. We investigated the response species of rhizosphere microorganisms and their response pathways under drought stress, providing a scientific basis for rhizosphere microorganisms to assist host plants to improve drought adaptability.

Minimal gene set from Sinorhizobium (Ensifer) meliloti pSymA required for efficient symbiosis with Medicago

Reduction of N₂ gas to ammonia in legume root nodules is a key component of sustainable agricultural systems. Root nodules are the result of a symbiosis between leguminous plants and bacteria called rhizobia. Both symbiotic partners play active roles in establishing successful symbiosis and nitrogen fixation: while root nodule development is mostly controlled by the plant, the rhizobia induce nodule formation, invade, and perform N₂ fixation once inside the plant cells. Many bacterial genes involved in the rhizobia-legume symbiosis are known, and there is much interest in engineering the symbiosis to include major nonlegume crops such as corn, wheat, and rice. We sought to identify and combine a minimal bacterial gene complement necessary and sufficient for symbiosis. We analyzed a model rhizobium, Sinorhizobium (Ensifer) meliloti, using a background strain in which the 1.35-Mb symbiotic megaplasmid pSymA was removed. Three regions representing 162 kb of pSymA were sufficient to recover a complete N₂-fixing symbiosis with alfalfa, and a targeted assembly of this gene complement achieved high levels of symbiotic N₂ fixation. The resulting gene set contained just 58 of 1,290 pSymA protein-coding genes. To generate a platform for future synthetic manipulation, the minimal symbiotic genes were reorganized into three discrete nod, nif, and fix modules. These constructs will facilitate directed studies toward expanding the symbiosis to other plant partners. They also enable forward-type approaches to identifying genetic components that may not be essential for symbiosis, but which modulate the rhizobium's competitiveness for nodulation and the effectiveness of particular rhizobia-plant symbioses.

Experimental Evolution in Plant-Microbe Systems: A Tool for Deciphering the Functioning and Evolution of Plant-Associated Microbial Communities

In natural environments, microbial communities must constantly adapt to stressful environmental conditions. The genetic and phenotypic mechanisms underlying the adaptive response of microbial communities to new (and often complex) environments can be tackled with a combination of experimental evolution and next generation sequencing. This combination allows to analyse the real-time evolution of microbial populations in response to imposed environmental factors or during the interaction with a host, by screening for phenotypic and genotypic changes over a multitude of identical experimental cycles. Experimental evolution (EE) coupled with comparative genomics has indeed facilitated the monitoring of bacterial genetic evolution and the understanding of adaptive evolution processes. Basically, EE studies had long been done on single strains, allowing to reveal the dynamics and genetic targets of natural selection and to uncover the correlation between genetic and phenotypic adaptive changes. However, species are always evolving in relation with other species and have to adapt not only to the environment itself but also to the biotic environment dynamically shaped by the other species. Nowadays, there is a growing interest to apply EE on microbial communities evolving under natural environments. In this paper, we provide a non-exhaustive review of microbial EE studies done with systems of increasing complexity (from single species, to synthetic communities and natural communities) and with a particular focus on studies between plants and plant-associated microorganisms. We highlight some of the mechanisms controlling the functioning of microbial species and their adaptive responses to environment changes and emphasize the importance of considering bacterial communities and complex environments in EE studies.

Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping

Rhizobium–legume symbioses serve as paradigmatic examples for the study of mutualism evolution. The genus Ensifer (syn. Sinorhizobium) contains diverse plant-associated bacteria, a subset of which can fix nitrogen in symbiosis with legumes. To gain insights into the evolution of symbiotic nitrogen fixation (SNF), and interkingdom mutualisms more generally, we performed extensive phenotypic, genomic, and phylogenetic analyses of the genus Ensifer. The data suggest that SNF likely emerged several times within the genus Ensifer through independent horizontal gene transfer events. Yet, the majority (105 of 106) of the Ensifer strains with the nodABC and nifHDK nodulation and nitrogen fixation genes were found within a single, monophyletic clade. Comparative genomics highlighted several differences between the “symbiotic” and “nonsymbiotic” clades, including divergences in their pangenome content. Additionally, strains of the symbiotic clade carried 325 fewer genes, on average, and appeared to have fewer rRNA operons than strains of the nonsymbiotic clade. Initial characterization of a subset of ten Ensifer strains identified several putative phenotypic differences between the clades. Tested strains of the nonsymbiotic clade could catabolize 25% more carbon sources, on average, than strains of the symbiotic clade, and they were better able to grow in LB medium and tolerate alkaline conditions. On the other hand, the tested strains of the symbiotic clade were better able to tolerate heat stress and acidic conditions. We suggest that these data support the division of the genus Ensifer into two main subgroups, as well as the hypothesis that pre-existing genetic features are required to facilitate the evolution of SNF in bacteria.

Transfer of the Symbiotic Plasmid of Rhizobium etli CFN42 to Endophytic Bacteria Inside Nodules

Conjugative transfer is one of the mechanisms allowing diversification and evolution of bacteria. Rhizobium etli CFN42 is a bacterial strain whose habitat is the rhizosphere and is able to form nodules as a result of the nitrogen-fixing symbiotic relationship it may establish with the roots of Phaseolus vulgaris. R. etli CFN42 contains one chromosome and six large plasmids (pRet42a - pRet42f). Most of the genetic information involved in the establishment of the symbiosis is localized on plasmid pRet42d, named as the symbiotic plasmid (pSym). This plasmid is able to perform conjugation, using pSym encoded transfer genes controlled by the RctA/RctB system. Another plasmid of CFN42, pRet42a, has been shown to perform conjugative transfer not only in vitro, but also on the surface of roots and inside nodules, using other rhizobia as recipients. In addition to the rhizobia involved in the formation of nodules, these structures have been shown to contain endophytic bacteria from different genera and species. In this work, we have explored the conjugative transfer of the pSym (pRet42d) from R. etli CFN42 to endophytic bacteria as putative recipients, using as donor a CFN42 derivative labeled with GFP in the pRet42d and RFP in the chromosome. We were able to isolate some transconjugants, which inherit the GFP, but not the RFP marker. Some of them were identified, analyzed and evaluated for their ability to nodulate. We found transconjugants from genera such as Stenotrophomonas, Achromobacter, and Bacillus, among others. Although all the transconjugants carried the GFP marker, and nod, fix, and nif genes from pRet42d, not all were able to nodulate. Ultrastructure microscopy analysis showed some differences in the structure of the nodules of one of the transconjugants. A replicon of the size of pRet42d (371 Kb) could not be visualized in the transconjugants, suggesting that the pSym or a segment of the plasmid is integrated in the chromosome of the recipients. These findings strengthen the proposal that nodules constitute a propitious environment for exchange of genetic information among bacteria, in addition to their function as structures where nitrogen fixation and assimilation takes place.

Genome-Wide Identification of the CrRLK1L Subfamily and Comparative Analysis of Its Role in the Legume-Rhizobia Symbiosis

The plant receptor-like-kinase subfamily CrRLK1L has been widely studied, and CrRLK1Ls have been described as crucial regulators in many processes in Arabidopsis thaliana (L.), Heynh. Little is known, however, about the functions of these proteins in other plant species, including potential roles in symbiotic nodulation. We performed a phylogenetic analysis of CrRLK1L subfamily receptors of 57 different plant species and identified 1050 CrRLK1L proteins, clustered into 11 clades. This analysis revealed that the CrRLK1L subfamily probably arose in plants during the transition from chlorophytes to embryophytes and has undergone several duplication events during its evolution. Among the CrRLK1Ls of legumes and A. thaliana, protein structure, gene structure, and expression patterns were highly conserved. Some legume CrRLK1L genes were active in nodules. A detailed analysis of eight nodule-expressed genes in Phaseolus vulgaris L. showed that these genes were differentially expressed in roots at different stages of the symbiotic process. These data suggest that CrRLK1Ls are both conserved and underwent diversification in a wide group of plants, and shed light on the roles of these genes in legume–rhizobia symbiosis.