Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Optimizing Rhizobium-legume symbioses by simultaneous measurement of rhizobial competitiveness and N2 fixation in nodules

View ORCID ProfileMarcela A. Mendoza-Suárez, View ORCID ProfileBarney A. Geddes, View ORCID ProfileCarmen Sánchez-Cañizares, View ORCID ProfileRicardo H. Ramírez-González, View ORCID ProfileCharlotte Kirchhelle, View ORCID ProfileBeatriz Jorrin, and View ORCID ProfilePhilip S. Poole
  1. aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  2. bJohn Innes Centre, Norwich Research Park, NR4 7UH Norwich, United Kingdom

See allHide authors and affiliations

PNAS May 5, 2020 117 (18) 9822-9831; first published April 21, 2020; https://doi.org/10.1073/pnas.1921225117
Marcela A. Mendoza-Suárez
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Marcela A. Mendoza-Suárez
Barney A. Geddes
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Barney A. Geddes
Carmen Sánchez-Cañizares
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Carmen Sánchez-Cañizares
Ricardo H. Ramírez-González
bJohn Innes Centre, Norwich Research Park, NR4 7UH Norwich, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ricardo H. Ramírez-González
Charlotte Kirchhelle
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Charlotte Kirchhelle
Beatriz Jorrin
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Beatriz Jorrin
Philip S. Poole
aDepartment of Plant Sciences, University of Oxford, OX1 3RB Oxford, United Kingdom;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Philip S. Poole
  • For correspondence: philip.poole@plants.ox.ac.uk
  1. Edited by Éva Kondorosi, Hungarian Academy of Sciences, Biological Research Centre, Szeged, Hungary, and approved March 12, 2020 (received for review December 3, 2019)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Increasing legume use in agriculture is a key requirement for the sustainable intensification of global farming systems. Maximizing the yield of legumes requires matching of the plant to elite rhizobia that are both competitive for nodulation and capable of high rates of nitrogen fixation. We have developed a high-throughput method that allows identification of elite strains at the single nodule level with the potential to revolutionize the search for elite indigenous rhizobia.

Abstract

Legumes tend to be nodulated by competitive rhizobia that do not maximize nitrogen (N2) fixation, resulting in suboptimal yields. Rhizobial nodulation competitiveness and effectiveness at N2 fixation are independent traits, making their measurement extremely time-consuming with low experimental throughput. To transform the experimental assessment of rhizobial competitiveness and effectiveness, we have used synthetic biology to develop reporter plasmids that allow simultaneous high-throughput measurement of N2 fixation in individual nodules using green fluorescent protein (GFP) and barcode strain identification (Plasmid ID) through next generation sequencing (NGS). In a proof-of-concept experiment using this technology in an agricultural soil, we simultaneously monitored 84 different Rhizobium leguminosarum strains, identifying a supercompetitive and highly effective rhizobial symbiont for peas. We also observed a remarkable frequency of nodule coinfection by rhizobia, with mixed occupancy identified in ∼20% of nodules, containing up to six different strains. Critically, this process can be adapted to multiple Rhizobium-legume symbioses, soil types, and environmental conditions to permit easy identification of optimal rhizobial inoculants for field testing to maximize agricultural yield.

  • nodulation
  • legume
  • rhizobium
  • competition
  • nitrogen-fixing

Legumes have evolved the remarkable ability to host N2 fixing bacteria, known as rhizobia, in specialized organs called root nodules. In a symbiotic partnership within the root nodule, legumes supply nutrients to rhizobia that fix N2 gas from the atmosphere into reduced forms that are supplied to the legume (1). These symbiotic associations have been estimated to fix ∼80% of the biologically fixed N2 in agricultural areas (2). Rhizobia have been studied for more than a 100 years because of their ability to increase yields of legumes crops. The first industrial production of rhizobial inoculants began by the end of the 19th century (3). Since then, farmers have applied rhizobial inoculants as an alternative to economically expensive chemical fertilizers and in organic crop systems. The bioavailable nitrogen that is generated in nodules of legumes benefits nonlegume crops grown in rotation or at the same time (3). In some cases, rhizobial inoculants that can significantly improve legume crop growth under controlled conditions fail to compete against native strains with inferior N2 fixing abilities, known as the “rhizobial competition problem” (4). In other cases, native rhizobia perform better at N2 fixation than commercial inoculants (5) or even N-fertilized controls (6), showing the importance of developing inoculants for local crops that take advantage of strains preadapted to thrive in local environments (tailor-made inoculants). A potential breakthrough in the use of inoculants is the recent development of a seed coating based on silk fibronin and trehalose that stabilizes rhizobia in saline and, possibly, arid environments (7). While this is an exciting development, to improve inoculants, the rhizobia coated onto seeds must be both competitive for nodule occupancy and effective in N2 fixation.

Since N2 fixation happens inside nodules, several indirect measurements of rhizobial effectiveness have been developed such as shoot dry weight (DW), acetylene reduction assays (ARA), total nitrogen content, 15N2 fixation, and many others (8). These techniques are not amenable for the analysis of individual nodules and, as a result, when screening for effective rhizobial strains, only one strain per plant is usually tested. In the case of nodule occupancy, different strategies have been used to assess the competitiveness of rhizobial strains. These include chromosomally encoded marker genes (9), naturally occurring or induced antibiotic-resistance genes (10), the use of PCR fingerprints which target specific genes (11), or metagenomic sequencing to monitor strain-level changes (12). Many of these methods are highly time-consuming, since they involve a complex procedure for chromosomal integration of genes. Currently available metagenomic sequencing approaches require very complex computational analyses and are dependent on the availability of reference genomes for strains of interest. Marker systems using plasmids are easier to introduce into rhizobia than chromosomal integrations and are not dependent on preexisting genome data. However, concerns about plasmid stability in the absence of antibiotic selection during plant assays have limited the use of plasmids to study competitiveness (13). Plasmids that show stability in the absence of antibiotic selection have now been developed, such as the stable promoter-probe and low-copy vector pJP2 (14), and recently adapted to Golden Gate cloning to enable high-throughput plasmid assembly (15).

Effectiveness at N2 fixation and competitiveness for nodulation are independent rhizobial traits (12, 16) and, due to the complexity of competition assays, there are few studies evaluating these two traits simultaneously. Therefore, strains that are both highly effective and competitive (elite strains) can take years to identify (17, 18). To improve the identification of elite rhizobial inoculants, it is crucial to increase the number of rhizobial strains that can be evaluated simultaneously in the presence of native rhizobia and concomitantly assess their N2 fixing ability (19, 20). To overcome these limitations, we present a plasmid-based library allowing simultaneous assessment of competitiveness and effectiveness for rhizobial strains. To achieve this goal, plasmids were developed with two key modules: 1) a module for evaluating rhizobial effectiveness that included a consensus nifH promoter (PsnifH) driving nodule-specific expression of green fluorescent protein (GFP) as reporter gene, and 2) a module for measuring bacterial competitiveness, consisting of a unique synthetic nucleotide sequence (Plasmid ID), allowing the identification of multiple rhizobial strains in a single experiment using next generation sequencing (NGS). Plasmids were assembled by high-throughput Golden Gate cloning (21) in the broad host-range vector pOGG026, which has excellent stability in rhizobia, even in the absence of antibiotics (15).

Results

A Synthetic nifH Promoter (PsnifH) Adapted for Use in Diverse Rhizobia.

When considering a molecular tool to measure effectiveness, we hypothesized that the expression of superfolder GFP (sfGFP) from a nifH promoter (which drives the expression of the genes encoding the nitrogenase complex, nifHDK) (22) would correlate with N2 fixation in individual nodules. The fluorescence signal from the nodules containing such a bioreporter should reflect the expression of nifH, as well as the size of the bacteroid population, which is related to nodule size that is associated with effectiveness of N2 fixation (23, 24). A promoter that is active only following nodule formation will also avoid fitness penalties prior to the invasion of the plant that might alter strain competitiveness. The well-studied region upstream of the genes encoding nitrogenase (nifHDK), which constitutes the nifH promoter (PnifH), has the characteristics of a type -24/-12 promoter recognized by the RNA polymerase sigma factor RpoN (σ54) (25). This promoter has been used for many years in studying Rhizobium competition and also as a symbiosis-induced promoter (15, 26, 27). To make a nifH promoter as broadly useful as possible, a 258-bp consensus promoter (PsnifH) was derived from alignment of 48 strains of the genus Rhizobium (SI Appendix, Fig. S1). It includes an upstream activator sequence (UAS) (TGT-N10-ACA), a binding sequence for RpoN σ54 (sigma factor for expression of all nif genes in nodules) (TGGCACG-N4-TTG), and a ribosome-binding site (RBS) (GAAGGAA) (SI Appendix, Fig. S2A) (25, 28, 29). The synthetic nifH promoter was synthesized as a “PU” module with BsaI restriction enzyme sites at either ends (SI Appendix, Fig. S2B), allowing Golden Gate assembly into the stable plasmid pOGG026 (15).

To evaluate the function of PsnifH in various rhizobial strains, the chromogenic markers gusA (β-glucuronidase) and celB (thermostable β-galactosidase) (26) were assembled downstream of the nifH promoter generating pOPS0262 and pOPS0263, respectively. They enable tracking of gene expression by staining GusA with Magenta-GlcA (pink) or CelB with X-Gal (blue). These plasmids were conjugated into 11 rhizobial strains and inoculated onto their respective host (Fig. 1A). Rhizobium leguminosarum bv. trifolii WSM1325 was used as a negative control since it has a different nifH promoter to the other R. leguminosarum strains (24). All mature nodules expected to show nifH promoter activity stained correctly after enzymatic reaction with Magenta-GlcA (for GusA) or X-Gal (for CelB) (Fig. 1A) in several Rhizobium and even members of the Sinorhizobium genus. Since GusA and CelB are useful markers for competition assays that are not dependent on expensive laboratory equipment, the efficacy of the plasmids in competition assays were compared to gusA or celB integrated in the chromosome (9). Pea plants were coinoculated at ratios: 10, 50, and 90% of RlvUPM791[pOPS0263] versus Rlv3841[pOPS0262]. At 21 days postinoculation (dpi), pea roots were sequentially double stained with Magenta-GlcA and X-Gal after thermal treatment (Fig. 1B and SI Appendix, Table S5). As found by ref. 9, RlvUPM791 was significantly more competitive than Rlv3841 (y intercept = 21.97; P < 0.0001 different from 0) (Fig. 1B), supporting the use of stable plasmids based on pOGG0026 to mark R. leguminosarum strains.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Design and validation of a PsnifH consensus promoter. (A) Common bean root system inoculated with CIAT899[pOPS0262] (celB under the control of PsnifH) stained with X-Gal after thermal treatment and results obtained from the full qualitative assessment of PsnifH performance assessing celB and gusA expression in nodules. (B) Competition assays for nodule formation, pea root system inoculated with UPM791[pOPS0263], and Rlv3841[pOPS0262] after sequential double-staining (n = 10). (C) Competition assays for nodule formation between Rlv3841[pOPS0381] (marked) and Rlv3841-WT (unmarked) strains on pea roots. GFP detection (scale, 3,000–30,000 cps). Analysis of competitiveness was performed by comparing percent of marked strain cells in the original mix inoculum with percent of the obtained marked nodules in each plant. Results (n = 5) showed that marked nodules were not significantly different from unmarked nodules. (D) Microscope images of pea nodules formed by Rlv3841sfGFP harboring a plasmid containing sfGFP under control of PsnifH: immature (i), mature (ii), and senescent nodule (iii). Bright-field (Left), fluorescence (GFP filter) (Center), and merged images (Right). (E) Maximum intensity projection of a confocal stack of GFP-expressing bacteroids formed by Rlv3841[pOPS0381] from crushed pea root nodules (Upper), and 3D-segmented bacteroids from the same confocal stack color-coded for volume (scale from 0 to 60 μm3) (Lower). Video of 3D structures in SI Appendix.

Next, the fluorescent proteins, sfGFP (30) and mCherry (31) were assembled with PsnifH, generating pOPS0381 and pOPS0382, respectively. The native Rlv3841 PnifH promoter driving sfGFP (pOPS0379) and mCherry (pOPS0380) were used as positive controls (15). These plasmids were conjugated into Rlv3841, inoculated onto pea plants and their roots analyzed using a NightOWL charge coupled device (CCD) camera to quantify green and red fluorescence expression. Nodule-specific expression of both sfGFP and mCherry occurred in the positive control and with PsnifH (SI Appendix, Fig. S3). Additionally, the Indigo software permitted automatic counting of nodules expressing either GFP or mCherry. However, as fluorescence of sfGFP was almost 10 times greater than mCherry, only sfGFP was used for subsequent experiments. Rlv3841 marked with sfGFP under control of PsnifH (Rlv3841[pOPS0381]), and unmarked Rlv3841 (WT) (Fig. 1C) were inoculated together on pea roots (Fig. 1C). The percentage of 21 dpi nodules containing marked strains was not significantly different from their percentage in the inoculum (Fig. 1C). The slope of the regression line was 0.9936 ± 0.096 (95% CI [0.7841, 1.203]), not significantly different from 1 (P = 0.9476), and the y intercept = −1.279 ± 5.475 (95% CI [13.21, 10.65]) is not significantly different from 0 (P = 0.8192), demonstrating pOPS0381 does not alter the competitiveness of Rlv3841.

PsnifH Promoting sfGFP Expression Is a Biosensor for Nitrogenase Activity.

Pea nodules formed by Rlv3841 marked with sfGFP under control of PsnifH (Fig. 1D) showed high GFP in the N2 fixing zone of mature nodules, but no expression in immature nodules. Furthermore, individual bacteroids of Rlv3841 expressing GFP from crushed nodules were visualized using confocal microscopy (Fig. 1E) and their three-dimensional (3D) shapes reconstructed. Individual bacteroids are coralloid in shape with volumes of approximately 10–60 μm3. These qualitative results confirmed the correlation of sfGFP expression from PsnifH with N2 fixation. To determine if nitrogenase activity correlates quantitatively with GFP expression from PsnifH, four R. leguminosarum strains with different levels of effectiveness at fixing nitrogen were used. WSM1475 and WSM1521 were previously identified with a low and high effectiveness through the analysis of biomass (32). Additionally, Rlv3841 and UPM791 were added as inoculants as they are well-studied strains and their behavior is broadly characterized. The plasmids pOPS0491, pOPS0503, pOPS0504, and pOPS0515 (differing only by Plasmid ID to pOPS0381) were conjugated into WSM1475, WSM1521, Rlv3841, and UPM791, and transconjugants inoculated on pea (Fig. 2A). Effectiveness for nitrogenase activity of the strains was determined by ARA at 28 dpi and 35 dpi (Fig. 2B) and shoot dry biomass (DW) at 35 dpi and 42 dpi (SI Appendix, Table S6). The high effectiveness of strain WSM1521 had been confirmed before (33) through glasshouse evaluation of dry matter yields. Therefore, WSM1521 was used as a positive control for nitrogenase activity, with the prediction that GFP fluorescence from WSM1521 nodules should be higher than WSM1475, Rlv3841, and UPM791. GFP detection, ARA, and shoot DW biomass were normalized by the number of nodules per plant (SI Appendix, Fig. S4). A two-tailed Pearson test shows a positive correlation for GFP fluorescence and ARA at 28 dpi of r = 0.8894 (95% CI [0.6444–0.9688]), P = 0.0001 and a R2 = 0.7910. At 35 dpi, the positive correlation was r = 0.8160 (95% CI [0.4234–0.9506]), P = 0.0022 and a R2 = 0.6659. Results from 28 dpi were grouped by strain in Fig. 2C. At both 28 dpi and 35 dpi, GFP fluorescence and nitrogenase activity correlated, with a peak at 28 dpi (Fig. 2C), in agreement with maximum ARA occurring at flowering in peas (34). Dry biomass measurements at 35 dpi and 42 dpi (SI Appendix, Table S6) showed that strain WSM1521 is highly effective and confirm that sfGFP fluorescence is an easy method for quantification of nitrogenase activity in individual nodules.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

PsnifH promoting sfGFP expression is a biosensor for nitrogenase activity. (A) Four strains with different effectiveness class were used to inoculate pea plants. (B) At 28 and 35 dpi, root systems were simultaneously analyzed for their GFP detection from nodules and for their nitrogenase activity by ARA. (C) Correlation of GFP (fluorescence, CPS) and ARA ([nmoles of C2H2 reduced][h−1]) at 28 dpi (n = 3).

Plasmid IDs Allow High-Throughput Identification of Rhizobial Strains.

Due to limitations in the number of strains able to be tested for competitiveness in traditional marker systems, we developed a system utilizing NGS technology to monitor the competitiveness of many strains simultaneously. Plasmids were therefore tagged with 96 unique 12-nt Golay barcodes (35), incorporated into a 57-bp Level 0 Golden Gate “T” module (SI Appendix, Table S3 and Fig. S5). The overall features are terminal BsaI enzyme recognition sites, Plasmid ID (Golay barcode) at nucleotides 12–23, and a universal primer binding site at nucleotides 24–46. Ninety-five individual Golden Gate Level 1 cloning one-tube one-step reactions were set up (BsaI digestion and ligation) (21) in a 96-well plate, where the only difference was the T module (Plasmid ID). (Fig. 3A and SI Appendix, Figs. S6 and S7). To test that Plasmid IDs allow multiple rhizobial strains to be marked and identified, strains UPM791gusA and Rlv3841celB, with marker genes integrated into their genomes, respectively, had Plasmid ID E10 (pOPS0548) and Plasmid ID G2 (pOPS0564) introduced by conjugation (Fig. 3B). UPM791gusA-E10 and Rlv3841celB-G2 were set to compete 1:1 for pea nodulation and the strain occupying each nodule identified by sequential double-staining with Magenta-GlcA and X-Gal after thermal treatment (Fig. 3C). Forty nodules were analyzed by NGS, 36 singly infected nodules and four mixed nodules as determined by Magenta-GlcA and X-Gal staining (Fig. 3D). The Plasmid IDs of all 40 samples matched the Magenta-GlcA and X-Gal staining. For the four mixed nodules, the percentages of IDs matched qualitatively with the relative colors of each nodule (Fig. 3E).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Plasmid IDs allow high-throughput identification of rhizobial strains. (A) Plasmid ID Library built by Golden Gate cloning. Promoter, marker gene, and backbone vector were identical, with the T module (Plasmid ID) as the unique identifier in each case. (B) Experimental design: Plasmid ID-E10 (pOPS0548) and Plasmid ID-G2 (pOPS0564) were conjugated into UPM791gusA and Rlv3841celB, respectively and used to coinoculate pea plants. (C) Pea roots that were sequentially double-stained with Magenta-GlcA and X-Gal after thermal treatment, resulting in pink nodules formed by UPM791gusA-E10[pOPS0548] (gusA constitutively expressed) and blue nodules formed by Rlv3841celB-G2[pOPS0564] (celB constitutively expressed). Four mixed nodules were identified (indicated by black arrows). (D) Forty nodules were selected and prepared for analysis by NGS with Ion Torrent, including the four mixed nodules on wells A1, B1, C1, and D1. (E) Sequencing results showing ID relative abundance per well sample.

Simultaneous High-Throughput Assessment of Nodulation Competitiveness and N2 Fixation Effectiveness.

Measuring effectiveness and competitiveness of Rhizobium strains in sterile conditions, deprived of N, limits the finding’s relevance to agricultural conditions. Therefore, a competition assay was performed by adding tagged-rhizobial strains to nonsterile nutrient-rich field soil containing indigenous rhizobia. Soil was from Yatesbury House Farm (Yatesbury, Calne, Wiltshire, United Kingdom). Soil was analyzed for chemical composition (SI Appendix, Tables S8 and S9) and by most probable number (MPN) (8) for rhizobial population size. Yatesbury soil contained ∼100 cells of indigenous rhizobia per gram, so that competition assays contained 6.20 × 104 rhizobial cells in 620.4 ± 1.96 g of Yatesbury soil. Ninety-five well-characterized R. leguminosarum bv. viciae (Rlv) strains were selected from different geographical locations (SI Appendix, Fig. S11) and tested for effective symbiosis with pea plants (SI Appendix, Fig. S12 A and B). The plasmid library (SI Appendix, Fig. S7) was grown in liquid culture in a 96-deepwell plate and then conjugated into the strains (SI Appendix, Fig. S11) (Fig. 4A), with 84 strains reaching the desired growth. The final Tagged-strains library was grown in liquid culture in a 96-deepwell plate, total cell density adjusted to ∼6.20 × 104 cells, and the combined Tagged-strains library inoculated onto pea seedlings in 0.5 L of Yatesbury soil (Fig. 4A). At 21 dpi (flowering time), pea plants were harvested (SI Appendix, Fig. S14A) and roots were analyzed (SI Appendix, Fig. S14 B and C). Plants (n = 30) formed an average of 54.3 nodules, with 57.4% (±2.6 SEM) expressing PsnifH GFP and 42.6% occupied by indigenous rhizobia (Fig. 4C). The IDs from GFP-expressing nodules were PCR-amplified using the conserved primer binding sites, sequenced by Ion Torrent and Plasmid IDs cross-referenced to strains (SI Appendix, Fig. S8). For a nodule to be scored as a mixed infection, a threshold for the second strain was set at 10%. After high-quality filtering (SI Appendix, Fig. S8), 769 sequenced nodules were identified, of which 611 nodules resulted from a single clonal infection and 158 contained more than one strain (i.e., mixed nodules). Even with such a conservative strain identification method (10% threshold), nodules with six strains were identified (SI Appendix, Table S10). Mixed nodules represent 20.6% of the total population (Fig. 4D), which is much higher than the ∼1% that we observed in competition assays when only two strains were inoculated in sterile vermiculite/sand (SI Appendix, Table S5). Furthermore, 20.6% is a lower estimate since indigenous rhizobia occupied 42.6% of nodules and mixtures of GFP-expressing tagged strains and indigenous rhizobia would not be identified. Remarkably, 15.2% of all sequenced nodules were occupied by strain G083, followed by strains V116 (7.02%), G067 (5.20%), G109 (5.07%) G0077 (4.42%), and Rlv3831 (4.29%).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

High-throughput competition assay among tagged-rhizobial strains in a nonsterile nutrient-rich field soil containing indigenous rhizobia. (A) Strategy for the competition assay among Tagged-strains in the presence of indigenous rhizobia. (B) Example of roots exposed to a blue-light transilluminator to check nodules expressing GFP belonged to Tagged-strains. (C) Average of nodules expressing GFP per plant (n = 30). (D) Percentage of single and mixed infection from total sequenced nodules. (E) Total occurrence of Tagged-strains identified from a total of 769 sequenced nodules.

Identification of a Supercompetitive and Highly Effective Rhizobial Strain.

We calculated a competitiveness index (CI) of each strain by comparing their own occurrence (the number of nodules in which the strains were identified from the 769 sequenced nodules) to the number of GFP-tagged nodules per plant:CI=[nAnT][100]CI=competitivenessnA=nodules occupied by strain AnT=total number of nodules analysedplant−1.Statistical analyses of the 10 strains with the highest CI (Fig. 5A) was performed and after ordinary one-way ANOVA, a Tukey’s multiple comparisons test showed that the CI of G083 (CI = 16.14) was higher (P < 0.0001) than any of the other strains (SI Appendix, Table S11). Additionally, G083 was the only strain that appeared in each of the 30 plants tested (Fig. 5A).

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Identification of a supercompetitive and highly effective rhizobial strain. (A) Competitiveness index (CI) of the top 10 strains. Dots represent CI of each strain per plant, and bars represent the SEM. (B) Comparison of the effectiveness of nodules occupied by different strains measured by GFP detection (cps) and their occurrence (number of nodules in which the strains were identified from the 240 most effective nodules); error bar represent SEM. (C) Effectiveness assessment of selected strains by ARA (n = 4). (D) Examples of 6-wk-old pea plants inoculated with selected strains. (E) Assessment of effectiveness by measurement of DW. WC; water control (n = 4).

All individual nodules were quantified for fluorescence while still on the plant (SI Appendix, Fig. S15A). The eight brightest nodules (i.e., the most effective) from each of the 30 plants analyzed (n = 240) were picked and gridded in microtiter plates, and their fluorescence was quantified (SI Appendix, Fig. S15B). DNA was then extracted, and Plasmid IDs determined by sequencing (Materials and Methods).

After ordinary one-way ANOVA and Tukey’s multiple comparisons, among the most effective nodules, G083 had the highest occurrence (n = 81) and showed significantly higher GFP fluorescence than nodules containing other strains (Fig. 5B), with the exception of L008 (n = 7), where the difference was marginally nonsignificant (P = 0.056) (SI Appendix, Table S12). These data suggested that in addition of being the most competitive strain, G083 is also a highly effective strain. While our high-throughput methodology represents a significant advance in simultaneous assessment of competitiveness and N2 fixing ability in Rhizobium-legume symbioses, can this be extrapolated to plant performance and growth yield? To answer this, the effectiveness of G083 was assessed by ARA at 21 dpi and by DW at 42 dpi. Along with G083, strains V116 and Rlv3841 were selected as controls—V116 as it had the second highest occurrence of the most effective nodules (n = 25) and Rlv3841 (n = 8) as the most well-studied strain of R. leguminosarum bv viciae. In both assays, ARA (Fig. 5C and SI Appendix, Table S13) and DW (Fig. 5 D and E and SI Appendix, Table S15) values for G083 were significantly higher than V116 and Rlv3841, while the latter strains did not show a significant difference from each other (ordinary one-way ANOVA and Dunnett’s multiple comparisons test) (SI Appendix, Tables S14 and S16). Additionally, G083 produced fewer but larger nodules and more pea pods per plant at 42 dpi (SI Appendix, Table S15). Both ARA and DW measurements show that the results obtained by high-throughput GFP detection in individual nodules correlates with plant performance and growth. Overall, the molecular toolkit developed here allows the identification of elite strains, such as G083 in pea, which is both supercompetitive and highly effective at N2 fixation.

Discussion

PsnifH as a Consensus Promoter for Rhizobium Strains.

We designed a synthetic consensus promoter (PsnifH) compatible with Golden Gate cloning to drive nodule-specific expression of CelB, GusA, mCherry, and sfGFP, which confirmed that PsnifH can be used in a broad range of R. leguminosarum strains and potentially other species of rhizobia (e.g., R. tropici, S. fredii) (Fig. 1A). Fluorescent proteins do not need cofactors or substrates (36), and additionally, sfGFP has a greater resistance to chemical denaturants than GFP (30). This could have been problematic if sfGFP fluorescence accumulated in the nodules over time. However, senescent areas of nodules did not show fluorescence, which indicates that the signal was restricted to the N2 fixing zone of nodules (Fig. 1 D, iii). Although we were able to show a strong positive correlation between fluorescence detection and ARA (r = 0.8894, Fig. 2C), this correlation is not strictly linear since the fluorescence detection from WSM1475 was lower compared to its results from ARA. These results do not invalidate the use of florescence detection as an easy and quantitative proxy for nitrogenase activity in individual nodules, but rather define the limit of the assay. There are more accurate measurements of total nitrogenase activity than GFP florescence, such as flow-through ARA systems and 15N2 assays, which will improve the linearity of correlation, but they are very laborious, time-consuming, and require expensive special equipment that is not accessible for low-income laboratories. Therefore, assessing a handful of strains is a major undertaking; moreover, assessing tens or hundreds of strains is almost impossible using these standard techniques.

For a rhizobial strain to be selected as a suitable inoculant, it must first be evaluated for its tolerance to environmental stresses, its N2 fixation effectiveness, its soil persistence, and most importantly, its ability to compete against indigenous strains for nodule formation in a directly relevant soil (37). The addition of the Plasmid IDs to the PsnifH sfGFP reporter plasmids provides the opportunity to assess the competitiveness of up to 95 strains simultaneously, although the number can be increased by building more Plasmid IDs since all parts are compatible with our recently developed bacterial expression vector archive (BEVA) bacterial cloning system (15).

Identification of a Supercompetitive and Highly Effective Rhizobial Strain.

Soil pH has been shown to affect the competitiveness of R. leguminosarum strains (38), which may explain the high competitiveness of G083 in our test soil. G083 was isolated from a soil with a pH 7.6 (39), which is very similar to that of Yatesbury soil (pH 7.75) (SI Appendix, Table S9). A key breakthrough of such high-throughput screening is that multiple soils, environmental conditions, and plant cultivars can be tested all at the same time (20). While G083 is undoubtedly an elite strain in our test soil, other strains may perform better under different conditions; e.g., of soil type, pH, moisture content, or nutrient status. High-throughput screening will enable matching of inoculant to field conditions.

One possible explanation for the competitive advantage of strain G083 is that it may migrate to root hairs faster and initiate infection much earlier than any other tested strains. This ability could be related to communication between the first strain that enters in the infection threads (ITs) and the plant itself, triggering an immune response that obstructs entry of further strains (40). For instance, necrosis of cortical cells of ITs in alfalfa is observed after S. meliloti infection (41). Alternatively, G083 may have active mechanisms of preventing other strains from infecting or of disabling them in a competitive race. Examples of such mechanisms include type VI secretion systems (T6SS), which secrete effector proteins (42) that may have antibacterial activity (43) as well a positive role in Rhizobium–legume symbioses (44). G083 was only found in nine mixed nodules (5.7% of mixed nodules), a very low occurrence compared to its presence in single-infection nodules (19.1% of single infection nodules). This suggests a supercompetitive strain may trigger autoregulation of nodulation (AON) much earlier. AON involves long-distance root-shoot signaling that inhibits further nodulation events (45⇓–47). We therefore presume that from the moment G083 began to form nodules, and fix N2 at high rates, it restricted less competitive strains from successful infection (45).

A remarkable observation from comparing nodulation of multiple strains in soil is that mixed nodules represent 20.6% of the total sequenced nodules. This percentage is much higher than the ∼1% observed in previous competition assays where only two strains were inoculated at a time in sterile vermiculite/sand (SI Appendix, Table S5). Furthermore, some of the nodules containing one tagged strain are likely to contain indigenous rhizobia that would not be detected. Tagged strain identification was also conservative with no mismatches allowed in either sample or Plasmid ID and a 10% threshold of total reads used for the presence of a strain. Although there are existing reports of high levels of mixed nodules, such as 36% reported in lentils (48), this tends to be the exception. Our results showing that inoculation of low numbers (i.e., same number as indigenous rhizobia in the soil) of multiple strains into an agricultural soil increases the mixed nodule occupancy from a few percent to over 20% has significant implications for the evolution and stability of Rhizobium-legume symbioses. Bacteroids of R. leguminosarum are terminally differentiated and do not regrow from nodules of peas due to the carefully orchestrated production of antimicrobial peptides by the plant (49⇓–51). Plants cannot select effective N2 fixing strains from the soil (40), instead they must sanction ineffective nodules, presumably by reducing resource allocation to a nodule and shutting it down. The result is that fewer undifferentiated bacteria are released from ineffective nodules compared to effective nodules (52). In legumes that produce antimicrobial peptides, this biased release of effective strains is thought to be a primary driver of the evolution of symbiosis. However, one way in which ineffective rhizobia can cheat on the symbiosis is to form mixed nodules with more effective strains. Therefore, the demonstration that the presence of a rich microbiota in soil promotes a large increase in mixed nodule occupancy suggests that less effective strains may benefit by “hitch hiking” with effective strains. Many strains also only appear in mixed nodules, suggesting there may be cooperative partnerships between strains. While the work here has only scratched the surface of such ecological investigations, it is clear that tagged rhizobial strain sets represent a powerful tool to address the mechanisms underlying competitive cooperation and exclusion.

Although we quantified GFP using the sensitive, but expensive, NightOWL CCD camera, it is possible to use low-cost equipment that is widely available. For example, by using a smart phone with a blue-light transilluminator and an orange filter, it is possible to obtain good quality photos (Fig. 4B) suitable for fluorescence image analysis with Fiji software (53). The only specialized equipment required is a PCR thermocycler to amplify the Plasmid IDs, and the samples can then be posted to one of the many low-cost sequencing centers around the world.

Until now, the search for elite strains has been a time-consuming process involving numerous studies in the laboratory before the strains could be assessed in greenhouse and field trials. The existing tagged PsnifH sfGFP plasmid library can be transferred to many other Rhizobium species for strain selection for improvement of plant performance or use in ecological studies of strain competition. By using our BEVA construction system, it is also easy to exchange the nifH promoter, IDs, or resistance markers, to make it more suitable for other rhizobia. While we have developed a plasmid system to facilitate its transfer, it could be adapted to Tn7 (54) for chromosomal integration in species of rhizobia where plasmid maintenance might be a problem.

In this work, we demonstrate a molecular approach to address the challenge of improving elite rhizobial inoculants for sustainable agriculture. It is further apparent that the tools developed here also have powerful applications to microbial ecology. We believe that the next crucial step in translating this research to a real-world benefit would be to integrate these findings into coordinated multidisciplinary studies with agronomists and breeders.

Materials and Methods

Bacterial Strains and Growth Conditions.

The bacterial strains and plasmids used in this study are listed in See SI Appendix, Table S1. Escherichia coli strains were grown in liquid or solid Luria–Bertani medium at 37 °C supplemented with appropriate antibiotics (µg·mL−1): ampicillin, 100; tetracycline, 10; and kanamycin, 20. For E. coli ST18 strains, 5‐aminolevulinic acid (ALA) was added at 50 µg·mL−1. All plasmids are available from Addgene (https://www.addgene.org) with the Addgene ID given in SI Appendix, Table S1. Rhizobial strains, listed in SI Appendix, Table S1, were grown at 28 °C in Tryptone-Yeast (TY) extract or Universal Minimal Salts (UMS) (55) with appropriate carbon and nitrogen sources at 10 mM unless otherwise stated. Antibiotics were used at the following concentrations (µg·mL−1) unless otherwise stated: gentamicin, 20; kanamycin, 20; neomycin, 40; nitrofurantoin, 10; rifampicin, 50; spectinomycin, 50; streptomycin 500; tetracycline (2 in UMS, 5 in TY).

Designing of a PsnifH Consensus Promoter and Plasmid ID.

To determine a consensus sequence for a nodule-specific promoter expressed in a wide variety of biovars and strains of R. leguminosarum and designed as a PU module for Golden Gate cloning, we followed the process described in SI Appendix, Fig. S1 and detailed in SI Appendix.

For the Plasmid ID, a unique 12-nt, error-correcting Golay barcode (35) was selected as suitable for strain identification (the 96 selected Golay barcodes are given in SI Appendix, Table S1) and the full design is in SI Appendix.

Plasmid ID Library.

A total of 95 individual and unique reporter plasmids were built by Golden Gate Level 1 cloning detailed in SI Appendix.

Primer Design for Multiplex Ion Torrent Sequencing Strategy.

We designed primers (SI Appendix, Table S4) for a two-step PCR. Full design is in SI Appendix.

Thermostable β-Galactosidase Activity from Cloned celB.

Roots were submerged upside down in a 50-mL Falcon tube containing 40 mL of phosphate buffer and kept at 70 °C for 1 h in order to destroy endogenous β-galactosidases. After cooling, 500 μL of X-Gal (5-bromo-4-chloro-3-indolyl-β-d-galactoside) solution (250 μg/mL) was added. Tubes were covered with aluminum foil and incubated overnight at 37 °C.

Sequential Double Staining to Detect both β-Glucuronidase and Thermostable β-Galactosidase Activity.

The gusA activity assay was run first with Magenta-GlcA (5-bromo-6-chloro-3-indolyl-β-d-glucuronide acid) solution (200 μg/mL) as substrate and, after overnight incubation at 28 °C, the thermostable β-galactosidase activity for celB was done overnight at 37 °C with X-Gal as indicated above (26).

Competition Assay.

Coinoculation with two strains.

Protocol is modified from ref. 9 and fully described in SI Appendix.

Coinoculation with multiple rhizobial strains.

Rhizobium cultures to be used as inoculum were grown to early stationary phase (OD600 < 0.6–0.8). The OD600 of a 200-μL sample from each strain was measured in 96-well plate using a FLUOStar Omega Microplate Reader (BMG). Cultures were adjusted to the desired OD and mixed in a flask. Inoculum was prepared by adding 1 mL of the Tagged-strain mix in 74 mL of 2.5× concentrated nitrogen-free rooting solution (55) and then distributed homogenously into 500-mL pots. Five-day-old similar-sized pea (Avolar) seedlings were transferred to preinoculated pots. Pea plants were grown without inoculum as negative controls (water control [WC]).

Analysis of Agricultural Soil for Competition Assays.

Agricultural soil was collected in mid-August 2017 from Yatesbury House Farm (Yatesbury, Calne, Wiltshire, UK). This is a 600-ha farm maintained under European Union organic regulations and certified by The Soil Association Certification Limited (SA Certification). The crops grown in this farm are wheat, barley, oats, field beans, and lentils are cultivated with legume-based mixtures of clover. Full details of the analysis are in SI Appendix. MPN of indigenous rhizobia from Yatesbury soil was carried out following the methodology described in ref. 8 and fully described in SI Appendix.

Image Acquisition.

For green and red fluorescence expression, fluorescence CCD images were acquired with a ring-light epi illumination accessory for NightOWL LB-983 camera and for full root system, a blue-light transilluminator (VWR) together with a Kodak Wratten Gelatin Filter no. 49 were used (specific parameters in SI Appendix). Dissecting microscope images were taken using a microscope Leica M165 FC and the digital camera Leica DFC310 FX (software LAS v4.5). Visible light filter, green fluorescent, and red fluorescent filters were used.

Confocal Microscopy for Bacteroids Analysis.

Confocal microscope images were acquired using an HCX PL APO CS 20x/0.7 IMM ultraviolet lens on a Leica SP 5. sfGFP was excited at 488 nm and detected at 500–600 nm. Three-dimensional segmentation was performed in MorphographX (56) (specific parameters in SI Appendix).

Bioinformatic Analysis and Availability of Materials.

Paired-end merging of amplicon reads was done with USEARCH 10.0.240 (57). Quality filtering of amplicon reads was carried out with Trimmomatic-0.36 (58) (steps fully described in SI Appendix), and the Plasmid ID was determined using a custom script available in: https://github.com/marcelamendoza/Plasmid-ID.

Data Manipulation and Statistical Analysis.

Data manipulation of PivotTables and conditional tables were done in Microsoft Excel, and data were processed in either the statistics program Prism8 (GraphPad Software) or The R project for Statistical computing (59).

Data Availability.

The data supporting the findings of the study are available in this article and its SI Appendix. Raw reads for this study have been deposited in the ENA at EMBL-EBI under accession no. PRJEB37539

Acknowledgments

This work was supported by the Biotechnology and Biological Sciences Research Council Grants BB/L011484/1, BB/N003608/1, and BB/N013387/1 (to P.S.P.). M.A.M.-S was recipient of the Mexican Government Grant Consejo Nacional de Ciencia y Tecnología 266954/399852. We thank John Howieson, Euan James, and Juan Imperial for rhizobial strains; Vanessa De Sousa Vieira for technical assistance funded by the University of Oxford Disability Advisory Service; Richard Ganlett for farm soil; and Alison East for critical revision of the manuscript.

Footnotes

  • ↵1To whom correspondence may be addressed. Email: philip.poole{at}plants.ox.ac.uk.
  • Author contributions: M.A.M.-S., B.A.G., C.S.-C., and P.S.P. designed research; M.A.M.-S., R.H.R.-G., and C.K. performed research; M.A.M.-S. contributed new reagents/analytic tools; M.A.M.-S., B.A.G., C.S.-C., R.H.R.-G., C.K., B.J., and P.S.P. analyzed data; and M.A.M.-S., B.A.G., C.S.-C., and P.S.P. wrote the paper.

  • The authors declare no competing interest.

  • This article is a PNAS Direct Submission.

  • Data deposition: Raw reads for this study have been deposited in the European Nucleotide Archive (ENA) at European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI) under accession no. PRJEB37539.

  • This article contains supporting information online at https://www.pnas.org/lookup/suppl/doi:10.1073/pnas.1921225117/-/DCSupplemental.

Published under the PNAS license.

References

  1. 1.↵
    1. J. I. Sprent,
    2. E. K. James
    , Legume evolution: Where do nodules and mycorrhizas fit in? Plant Physiol. 144, 575–581 (2007).
    OpenUrlFREE Full Text
  2. 2.↵
    1. G. W. O’Hara
    , The role of nitrogen fixation in crop production. J. Crop Prod. 1, 115–138 (1998).
    OpenUrl
  3. 3.↵
    1. J. J. Parnell et al
    ., From the lab to the farm: An industrial perspective of plant beneficial microorganisms. Front. Plant Sci. 7, 1110 (2016).
    OpenUrl
  4. 4.↵
    1. R. J. Yates,
    2. J. G. Howieson,
    3. W. G. Reeve,
    4. G. W. O’Hara
    , A re-appraisal of the biology and terminology describing rhizobial strain success in nodule occupancy of legumes in agriculture. Plant Soil 348, 255–267 (2011).
    OpenUrlCrossRef
  5. 5.↵
    1. A. M. Chibeba,
    2. S. Kyei-Boahen,
    3. M. F. Guimarães,
    4. M. A. Nogueira,
    5. M. Hungria
    , Isolation, characterization and selection of indigenous Bradyrhizobium strains with outstanding symbiotic performance to increase soybean yields in Mozambique. Agric. Ecosyst. Environ. 246, 291–305 (2017).
    OpenUrl
  6. 6.↵
    1. D. Mulas et al
    ., Distribution and efficiency of Rhizobium leguminosarum strains nodulating Phaseolus vulgaris in Northern Spanish soils: Selection of native strains that replace conventional N fertilization. Soil Biol. Biochem. 43, 2283–2293 (2011).
    OpenUrl
  7. 7.↵
    1. A. T. Zvinavashe,
    2. E. Lim,
    3. H. Sun,
    4. B. Marelli
    , A bioinspired approach to engineer seed microenvironment to boost germination and mitigate soil salinity. Proc. Natl. Acad. Sci. U.S.A. 116, 25555–25561 (2019).
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. J. G. Howieson,
    2. M. J. Dilworth
    1. J. G. Howieson,
    2. M. J. Dilworth
    , Working with Rhizobia, J. G. Howieson, M. J. Dilworth, Eds. (ACIAR, 2016), pp. 97–105.
  9. 9.↵
    1. C. Sánchez-Cañizares,
    2. J. Palacios
    , Construction of a marker system for the evaluation of competitiveness for legume nodulation in Rhizobium strains. J. Microbiol. Methods 92, 246–249 (2013).
    OpenUrl
  10. 10.↵
    1. B. A. Geddes,
    2. I. J. Oresnik
    , Inability to catabolize galactose leads to increased ability to compete for nodule occupancy in Sinorhizobium meliloti. J. Bacteriol. 194, 5044–5053 (2012).
    OpenUrlAbstract/FREE Full Text
  11. 11.↵
    1. S. Duodu,
    2. C. Brophy,
    3. J. Connolly,
    4. M. M. Svenning
    , Competitiveness of a native Rhizobium leguminosarum biovar trifolii strain for nodule occupancy is manifested during infection. Plant Soil 318, 117–126 (2009).
    OpenUrlCrossRef
  12. 12.↵
    1. L. T. Burghardt et al
    ., Select and resequence reveals relative fitness of bacteria in symbiotic and free-living environments. Proc. Natl. Acad. Sci. U.S.A. 115, 2425–2430 (2018).
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. V. Corich et al
    ., Aspects of marker/reporter stability and selectivity in soil microbiology. Microb. Ecol. 41, 333–340 (2001).
    OpenUrlPubMed
  14. 14.↵
    1. F. Pini et al
    ., Bacterial biosensors for in vivo spatiotemporal mapping of root secretion. Plant Physiol. 174, 1289–1306 (2017).
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. B. A. Geddes,
    2. M. A. Mendoza-Suárez,
    3. P. S. Poole
    , A bacterial expression vector archive (BEVA) for flexible modular assembly of golden gate-compatible vectors. Front. Microbiol. 9, 3345 (2019).
    OpenUrl
  16. 16.↵
    1. A. A. Cardoso,
    2. M. P. Andraus,
    3. T. C. Borba,
    4. C. C. Martin-Didonet,
    5. E. P. Ferreira
    , Characterization of rhizobia isolates obtained from nodules of wild genotypes of common bean. Braz. J. Microbiol. 48, 43–50 (2017).
    OpenUrl
  17. 17.↵
    1. V. Bourion et al
    ., Co-inoculation of a pea core-collection with diverse rhizobial strains shows competitiveness for nodulation and efficiency of nitrogen fixation are distinct traits in the interaction. Front. Plant Sci. 8, 2249 (2018).
    OpenUrlCrossRef
  18. 18.↵
    1. P. Irisarri et al
    ., Selection of competitive and efficient rhizobia strains for white clover. Front. Microbiol. 10, 768 (2019).
    OpenUrl
  19. 19.↵
    1. A. Checcucci,
    2. G. C. DiCenzo,
    3. M. Bazzicalupo,
    4. A. Mengoni
    , Trade, diplomacy, and warfare: The quest for elite rhizobia inoculant strains. Front. Microbiol. 8, 2207 (2017).
    OpenUrl
  20. 20.↵
    1. A. Greenlon et al
    ., Global-level population genomics reveals differential effects of geography and phylogeny on horizontal gene transfer in soil bacteria. Proc. Natl. Acad. Sci. U.S.A. 116, 15200–15209 (2019).
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. C. Engler,
    2. R. Kandzia,
    3. S. Marillonnet
    , A one pot, one step, precision cloning method with high throughput capability. PLoS One 3, e3647 (2008).
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. L. M. Rubio,
    2. P. W. Ludden
    , Biosynthesis of the iron-molybdenum cofactor of nitrogenase. Annu. Rev. Microbiol. 62, 93–111 (2008).
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. J. U. Regus et al
    ., Cell autonomous sanctions in legumes target ineffective rhizobia in nodules with mixed infections. Am. J. Bot. 104, 1299–1312 (2017).
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. V. J. Melino et al
    ., Identifying abnormalities in symbiotic development between Trifolium spp. and Rhizobium leguminosarum bv. trifolii leading to sub-optimal and ineffective nodule phenotypes. Ann. Bot. 110, 1559–1572 (2012).
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. B. Thöny,
    2. H. Hennecke
    , The -24/-12 promoter comes of age. FEMS Microbiol. Rev. 5, 341–357 (1989).
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. A. Sessitsch,
    2. K. J. Wilson,
    3. A. D. Akkermans,
    4. W. M. de Vos
    , Simultaneous detection of different Rhizobium strains marked with either the Escherichia coli gusA gene or the Pyrococcus furiosus celB gene. Appl. Environ. Microbiol. 62, 4191–4194 (1996).
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. M. D. Denton,
    2. W. G. Reeve,
    3. J. G. Howieson,
    4. D. R. Coventry
    , Competitive abilities of common field isolates and a commercial strain of Rhizobium leguminosarum bv. trifolii for clover nodule occupancy. Soil Biol. Biochem. 35, 1039–1048 (2003).
    OpenUrl
  28. 28.↵
    1. J. T. Sullivan,
    2. S. D. Brown,
    3. C. W. Ronson
    , The NifA-RpoN regulon of Mesorhizobium loti strain R7A and its symbiotic activation by a novel LacI/GalR-family regulator. PLoS One 8, e53762 (2013).
    OpenUrlCrossRef
  29. 29.↵
    1. B. Brito et al
    ., Hydrogenase genes from Rhizobium leguminosarum bv. viciae are controlled by the nitrogen fixation regulatory protein nifA. Proc. Natl. Acad. Sci. U.S.A. 94, 6019–6024 (1997).
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. J.-D. Pédelacq,
    2. S. Cabantous,
    3. T. Tran,
    4. T. C. Terwilliger,
    5. G. S. Waldo
    , Engineering and characterization of a superfolder green fluorescent protein. Nat. Biotechnol. 24, 79–88 (2006).
    OpenUrlCrossRefPubMed
  31. 31.↵
    1. N. C. Shaner et al
    ., Improved monomeric red, orange and yellow fluorescent proteins derived from Discosoma sp. red fluorescent protein. Nat. Biotechnol. 22, 1567–1572 (2004).
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. J. Howieson,
    2. A. Loi,
    3. S. Carr
    , Biserrula pelecinus L.-a legume pasture species with potential for acid, duplex soils which is nodulated by unique root-nodule bacteria. Aust. J. Agric. Res. 46, 997 (1995).
    OpenUrlCrossRef
  33. 33.↵
    1. J. Slattery,
    2. D. Pearce
    , Development of elite inoculant Rhizobium strains in Southeastern Australia. ACIAR Proc. 109e, 86–94 (2002).
    OpenUrl
  34. 34.↵
    1. G. Mulley et al
    ., Pyruvate is synthesized by two pathways in pea bacteroids with different efficiencies for nitrogen fixation. J. Bacteriol. 192, 4944–4953 (2010).
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. J. G. Caporaso et al
    ., Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 6, 1621–1624 (2012).
    OpenUrlCrossRefPubMed
  36. 36.↵
    1. M. Chalfie,
    2. Y. Tu,
    3. G. Euskirchen,
    4. W. W. Ward,
    5. D. C. Prasher
    , Green fluorescent protein as a marker for gene expression. Science 263, 802–805 (1994).
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. M. J. Dilworth,
    2. E. K. James,
    3. J. I. Sprent,
    4. W. E. Newton
    1. M. J. Dilworth,
    2. E. K. James,
    3. J. I. Sprent,
    4. W. E. Newton
    , Nitrogen-fixing Leguminous Symbioses, M. J. Dilworth, E. K. James, J. I. Sprent, W. E. Newton, Eds. (Springer Netherlands, 2008), p. 88.
  38. 38.↵
    1. S. D. Frey,
    2. L. K. Blum
    , Effect of pH on competition for nodule occupancy by type I and type II strains of Rhizobium leguminosarum bv. phaseoli. Plant Soil 163, 157–164 (1994).
    OpenUrl
  39. 39.↵
    1. B. Jorrin Rubio
    , “Genomics of Specificity in the Symbiotic Interaction between Rhizobium leguminosarum and Legumes,” PhD thesis, Universidad Politécnica de Madrid (2016). https://doi.org/10.20868/UPM.thesis.43380. Accessed 25 September 2018.
  40. 40.↵
    1. A. Westhoek et al
    ., Policing the legume-Rhizobium symbiosis: A critical test of partner choice. Sci. Rep. 7, 1419 (2017).
    OpenUrlCrossRef
  41. 41.↵
    1. J. Vasse,
    2. F. Billy,
    3. G. Truchet
    , Abortion of infection during the Rhizobium meliloti-alfalfa symbiotic interaction is accompanied by a hypersensitive reaction. Plant J. 4, 555–566 (1993).
    OpenUrlCrossRef
  42. 42.↵
    1. M. Gallique,
    2. M. Bouteiller,
    3. A. Merieau
    , The type VI secretion system: A dynamic system for bacterial communication? Front. Microbiol. 8, 1454 (2017).
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. A. B. Russell et al
    ., Diverse type VI secretion phospholipases are functionally plastic antibacterial effectors. Nature 496, 508–512 (2013).
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. A. Salinero-Lanzarote et al
    ., The Type VI secretion system of Rhizobium etli Mim1 has a positive effect in symbiosis. FEMS Microbiol. Ecol. 95, fiz054 (2019).
    OpenUrlCrossRef
  45. 45.↵
    1. V. Mortier,
    2. M. Holsters,
    3. S. Goormachtig
    , Never too many? How legumes control nodule numbers. Plant Cell Environ. 35, 245–258 (2012).
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. H. Ryu,
    2. H. Cho,
    3. D. Choi,
    4. I. Hwang
    , Plant hormonal regulation of nitrogen-fixing nodule organogenesis. Mol. Cells 34, 117–126 (2012).
    OpenUrlCrossRefPubMed
  47. 47.↵
    1. B. J. Ferguson et al
    ., Molecular analysis of legume nodule development and autoregulation. J. Integr. Plant Biol. 52, 61–76 (2010).
    OpenUrlCrossRefPubMed
  48. 48.↵
    1. S. N. May,
    2. B. B. Bohlool
    , Competition among Rhizobium leguminosarum strains for nodulation of Lentils (Lens esculenta). Appl. Environ. Microbiol. 45, 960–965 (1983).
    OpenUrlAbstract/FREE Full Text
  49. 49.↵
    1. P. Mergaert et al
    ., A novel family in Medicago truncatula consisting of more than 300 nodule-specific genes coding for small, secreted polypeptides with conserved cysteine motifs. Plant Physiol. 132, 161–173 (2003).
    OpenUrlAbstract/FREE Full Text
  50. 50.↵
    1. W. Van De Velde, et al
    ., Plant peptides govern terminal differentiation of bacteria in symbiosis. Science 327, 1122–1126 (2010).
    OpenUrlAbstract/FREE Full Text
  51. 51.↵
    1. J. Montiel et al
    ., Morphotype of bacteroids in different legumes correlates with the number and type of symbiotic NCR peptides. Proc. Natl. Acad. Sci. U.S.A. 114, 5041–5046 (2017).
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. B. Daubech et al
    ., Spatio-temporal control of mutualism in legumes helps spread symbiotic nitrogen fixation. eLife 6, e28683 (2017).
    OpenUrl
  53. 53.↵
    1. J. Schindelin et al
    ., Fiji: An open-source platform for biological-image analysis. Nat. Methods 9, 676–682 (2012).
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. L. Romero-Jiménez,
    2. D. Rodríguez-Carbonell,
    3. M. T. Gallegos,
    4. J. Sanjuán,
    5. D. Pérez-Mendoza
    , Mini-Tn7 vectors for stable expression of diguanylate cyclase PleD* in Gram-negative bacteria. BMC Microbiol. 15, 190 (2015).
    OpenUrl
  55. 55.↵
    1. P. S. Poole,
    2. N. A. Schofield,
    3. C. J. Reid,
    4. E. M. Drew,
    5. D. L. Walshaw
    , Identification of chromosomal genes located downstream of dctD that affect the requirement for calcium and the lipopolysaccharide layer of Rhizobium leguminosarum. Microbiology 140, 2797–2809 (1994).
    OpenUrlCrossRefPubMed
  56. 56.↵
    1. P. Barbier de Reuille et al
    ., MorphoGraphX: A platform for quantifying morphogenesis in 4D. eLife 4, 05864 (2015).
    OpenUrlCrossRefPubMed
  57. 57.↵
    1. R. C. Edgar
    , Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461 (2010).
    OpenUrlCrossRefPubMed
  58. 58.↵
    1. A. M. Bolger,
    2. M. Lohse,
    3. B. Usadel
    , Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. R Core Team
    , R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, Austria, 2014). http://www.R-project.org/. Accessed 14 April 2020.
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Optimizing Rhizobium-legume symbioses by simultaneous measurement of rhizobial competitiveness and N2 fixation in nodules
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Optimizing Rhizobium-legume symbioses by simultaneous measurement of rhizobial competitiveness and N2 fixation in nodules
Marcela A. Mendoza-Suárez, Barney A. Geddes, Carmen Sánchez-Cañizares, Ricardo H. Ramírez-González, Charlotte Kirchhelle, Beatriz Jorrin, Philip S. Poole
Proceedings of the National Academy of Sciences May 2020, 117 (18) 9822-9831; DOI: 10.1073/pnas.1921225117

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Optimizing Rhizobium-legume symbioses by simultaneous measurement of rhizobial competitiveness and N2 fixation in nodules
Marcela A. Mendoza-Suárez, Barney A. Geddes, Carmen Sánchez-Cañizares, Ricardo H. Ramírez-González, Charlotte Kirchhelle, Beatriz Jorrin, Philip S. Poole
Proceedings of the National Academy of Sciences May 2020, 117 (18) 9822-9831; DOI: 10.1073/pnas.1921225117
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Agricultural Sciences
Proceedings of the National Academy of Sciences: 117 (18)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Setting sun over a sun-baked dirt landscape
Core Concept: Popular integrated assessment climate policy models have key caveats
Better explicating the strengths and shortcomings of these models will help refine projections and improve transparency in the years ahead.
Image credit: Witsawat.S.
Model of the Amazon forest
News Feature: A sea in the Amazon
Did the Caribbean sweep into the western Amazon millions of years ago, shaping the region’s rich biodiversity?
Image credit: Tacio Cordeiro Bicudo (University of São Paulo, São Paulo, Brazil), Victor Sacek (University of São Paulo, São Paulo, Brazil), and Lucy Reading-Ikkanda (artist).
Syrian archaeological site
Journal Club: In Mesopotamia, early cities may have faltered before climate-driven collapse
Settlements 4,200 years ago may have suffered from overpopulation before drought and lower temperatures ultimately made them unsustainable.
Image credit: Andrea Ricci.
Steamboat Geyser eruption.
Eruption of Steamboat Geyser
Mara Reed and Michael Manga explore why Yellowstone's Steamboat Geyser resumed erupting in 2018.
Listen
Past PodcastsSubscribe
Birds nestling on tree branches
Parent–offspring conflict in songbird fledging
Some songbird parents might improve their own fitness by manipulating their offspring into leaving the nest early, at the cost of fledgling survival, a study finds.
Image credit: Gil Eckrich (photographer).

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490