# Integrated roles of BclA and DD-carboxypeptidase 1 in Bradyrhizobium differentiation within NCR-producing and NCR-lacking root nodules

Quentin Barrière, Ibtissem Guefrachi, Djamel Gully, Florian Lamouche, Olivier Pierre, Joël Fardoux, Clémence Chaintreuil, Benoît Alunni, Tatiana Timchenko, Eric Giraud, et al.

► To cite this version:

Quentin Barrière, Ibtissem Guefrachi, Djamel Gully, Florian Lamouche, Olivier Pierre, et al. Integrated roles of BclA and DD-carboxypeptidase 1 in Bradyrhizobium differentiation within NCR-producing and NCR-lacking root nodules. Scientific Reports, 2017, 7 (1), pp.9063. 10.1038/s41598-017-08830-0. hal-02187766

HAL Id: hal-02187766

https://hal.science/hal-02187766

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# Integrated roles of BclA and DD-carboxypeptidase 1 in *Bradyrhizobium* differentiation within NCR-producing and NCR-lacking root nodules

Quentin Barrière<fnref n="1" />, Ibtissem Guefrachi<fnref n="1" />,²,⁴, Djamel Gully³, Florian Lamouche<fnref n="1" />, Olivier Pierre<fnref n="1" />,⁵, Joël Fardoux³, Clémence Chaintreuil³, Benoît Alunni<fnref n="1" />, Tatiana Timchenko<fnref n="1" />, Eric Giraud³ & Peter Mergaert<fnref n="1" />

Legumes harbor in their symbiotic nodule organs nitrogen fixing rhizobium bacteria called bacteroids. Some legumes produce Nodule-specific Cysteine-Rich (NCR) peptides in the nodule cells to control the intracellular bacterial population. NCR peptides have antimicrobial activity and drive bacteroids toward terminal differentiation. Other legumes do not produce NCR peptides and their bacteroids are not differentiated. Bradyrhizobia, infecting NCR-producing *Aeschynomene* plants, require the peptide uptake transporter BclA to cope with the NCR peptides as well as a specific peptidoglycan-modifying DD-carboxypeptidase, DD-CPase1. We show that *Bradyrhizobium diazoefficiens* strain USDA110 forms undifferentiated bacteroids in NCR-lacking soybean nodules. Unexpectedly, in *Aeschynomene afraspera* nodules the nitrogen fixing USDA110 bacteroids are hardly differentiated despite the fact that this host produces NCR peptides, suggesting that USDA110 is insensitive to the host peptide effectors and that nitrogen fixation can be uncoupled from differentiation. In agreement with the absence of bacteroid differentiation, USDA110 does not require its *bclA* gene for nitrogen fixing symbiosis with these two host plants. Furthermore, we show that the BclA and DD-CPase1 act independently in the NCR-induced morphological differentiation of bacteroids. Our results suggest that BclA is required to protect the rhizobia against the NCR stress but not to induce the terminal differentiation pathway.

Nodules on the roots of legume plants - and in certain cases also on the stems - are symbiotic organs that house within their cells thousands of nitrogen-fixing rhizobium bacteria, called bacteroids<fnref n="1" />,². This intimate co-habitation of eukaryotic cells with bacteria requires adaptations of both the host and the micro-symbiont to maintain a cooperative state to avoid one partner succumbing to the aggression or immune responses of the other. The production of effector peptides, called Nodule-specific Cysteine-Rich (NCR) peptides by the Inverted Repeat Lacking Clade (IRLC) and Dalbergoid legumes is such a host adaptation that helps control the endosymbiotic bacterial population³⁻⁵. These peptides, related to innate immunity antimicrobial peptides, target the bacteria and induce them into a terminally differentiated state which is characterized by a pronounced cell enlargement resulting in either elongated or spherical bacteroids, depending on the host species. This bacterial cell enlargement is accompanied by a strong polyploidization of the bacterial genome.

Received: 18 May 2017

Accepted: 12 July 2017

Published online: 22 August 2017

<footnote n="1">¹Institute for Integrative Biology of the Cell, UMR9198, CNRS, Université Paris-Sud, CEA, Gif-sur-Yvette, France.
²Research Unit Biodiversity & Valorization of Arid Areas Bioressources (BVBA), Faculty of Sciences, Gabès University, Erriadh-Zrig, 6072, Gabès, Tunisia.
³Laboratoire des Symbioses Tropicales et Méditerranéennes, Institut pour la Recherche et le Développement, UMR IRD/SupAgro/INRA/UM2/CIRAD, Campus International de Baillarguet, TA A-82/J, 34398, Montpellier Cedex 5, France.
⁴Present address: Université de Pau et des Pays de l'Adour, Pau, France.
⁵Present address: Institut Sophia AgroBiotech, Sophia-Antipolis, France. Quentin Barrière and Ibtissem Guefrachi contributed equally to this work. Correspondence and requests for materials should be addressed to P.M. (email: peter.mergaert@i2bc.paris-saclay.fr)</footnote>
Some NCR peptides, mostly cationic ones, have the capacity to kill rhizobia, as well as other bacteria, in agree-
ment with their similarity to antimicrobial peptides such as defensins. Similarly to defensins, the NCR peptides
affect bacterial envelope permeability which may lead to bacterial death³. In addition to their effects on the bacte-
rial membranes, NCR peptides also have diverse intracellular targets⁶, including the cell cycle machinery, which
leads to the bacteroid genome amplification⁷.

NCR-induced bacterial elongation and polyploidization and bacterial survival in the face of NCR toxicity require a SbmA_BacA transmembrane domain-containing peptide importer protein which is able to transport the NCR peptides within the bacterial cytosol. It was proposed that NCR peptide uptake reduces their toxicity by preventing their deleterious effects on the bacterial membrane⁸,⁹. Depending on the rhizobium species, two types of NCR transporters have been identified. In *Sinorhizobium meliloti*, the symbiont of the IRLC legumes *Medicago sativa* and *Medicago truncatula*, the NCR transporter, called BacA, consists solely of the SbmA_BacA trans-membrane domain and is highly homologous to *E. coli* SbmA which drives peptide transport across the inner membrane by the membrane electrochemical gradient¹⁰. On the other hand, *Bradyrhizobium* spp. nodulating *Aeschynomene* spp. (Dalbergoids) use for NCR uptake the SbmA_BacA transmembrane domain-containing ABC transporter BclA which has, unlike BacA, a C-terminal ATPase domain and drives transport by ATP hydrolysis⁹. The SbmA_BacA transmembrane domain transporters are wide-spread in bacteria and are not symbiosis-specific proteins. Their bacterial household function remains poorly defined but in *E. coli* or *S. meliloti* they can import a diversity of peptides or peptide derivatives, including antimicrobial peptides that have intracellular targets¹¹⁻²².

Many legumes do not produce NCR peptides and the bacteroids in those legumes do not enlarge, do not become polyploid and are not terminally differentiated<sup>23, 24</sup>. This is the case for legumes of the Robinioids and Millettioid clades such as *Lotus japonicus*, *Phaseolus vulgaris*, *Lablab purpureus*, *Tephrosia vogelii*, *Vigna unguiculata* and soybean (*Glycine max*). Consequently, the *bacA* gene of their symbionts, *Mesorhizobium loti*, *Rhizobium leguminosarum* bv. *phaseoli*, *Rhizobium etli* and *Sinorhizobium fredii*, respectively, can be inactivated without consequences for the symbiosis<sup>1, 8, 25-28</sup>. Based on these case studies, it was proposed that the *bacA* gene in rhizobia has only a symbiotic role when the host plant produces NCR peptides<sup>8</sup>. However, the role of the bradyrhizobial *bclA* gene in interactions with host plants that do not produce NCR peptides has not been studied.

Besides the SbmA_BacA transporters BacA or BclA, a peptidoglycan-modifying enzyme, a specific DD-carboxypeptidase (DD-CPase) enzyme with a peptidoglycan-binding SPOR domain and encoded by the *DD-CPase1* gene, was recently identified in bradyrhizobia as a second rhizobial factor involved in terminal bacteroid differentiation in legumes producing NCR peptides<sup>29</sup>. DD-CPases are periplasmic enzymes that regulate the degree of reticulation in the peptidoglycan layer and thereby affect the physical strength of the bacterial wall. In the NCR-producing legumes *Aeschynomene indica* and *Aeschynomene afraspera*, the bacteroids of the *DD-CPase1* mutant are malformed and hypertrophied. However, in the NCR-lacking soybean, the mutation does not affect the bacteroids<sup>29</sup>.

In this study, we investigated whether *Bradyrhizobium diazoefficiens* (previously *Bradyrhizobium japonicum*) strain USDA110 requires its *bclA* gene in symbiosis with NCR-lacking and NCR-producing legume hosts. We also analyzed how the *bclA* and *DD-CPase1* determinants of terminal bacteroid differentiation interfere with each other, by creating and analyzing double mutants in the *bclA* and *DD-CPase1* genes in *B. diazoefficiens* strain USDA110 and *Bradyrhizobium* sp. strain ORS285.

Results

**The Bradyrhizobium diazoefficiens USDA110 blr7537-encoded protein is a BclA homologue.** Similar to other analyzed *Bradyrhizobium* strains⁹, no *bacA* homologous gene was identified in the genome of *B. diazoefficiens* USDA110. However, the *B. diazoefficiens* gene blr7537 was identified as a homolog of the *bclA* genes of *Bradyrhizobium* spp. ORS285 and ORS278⁹. The encoded protein has 70–72% identity and 80–81% similarity to the BclA proteins of strains ORS285 and ORS278. The protein has the same structure with an N-terminal SbmA-BacA transmembrane domain and a C-terminal ATPase domain. It thus encodes a potentially functional ABC transporter. Moreover, the genes in the 3 *Bradyrhizobium* species are located in syntenic regions that extend to over 200 kb. Therefore we designated blr7537 as *bclA*. Similarly as in strains ORS285 and ORS278, the *B. diazoefficiens* locus lacks genes that potentially encode additional components of the ABC transporter, such as a periplasmic binding protein for substrate binding and delivery to the transporter.

**Bradyrhizobium diazoefficiens USDA110 BclA is an NCR peptide transporter and functional in symbiosis.** To test the *in vitro* and *in vivo* activity of the *B. diazoefficiens* USDA110 BclA protein, a deletion mutant of the *bclA* gene was constructed. The gene was also cloned into the plasmids pMG103 and pRF771, downstream of the *trp* promoter, and introduced into the *bclA*, *bacA* and *sbmA* mutants of *Bradyrhizobium* sp. ORS285, *S. meliloti* strain Sm1021 and *E. coli* strain BW25113, respectively.

Bleomycin and Bac7 are antimicrobial compounds which have intracellular targets, DNA and the ribosomes respectively, and they require active transport to be taken up in the bacterial cells. In *E. coli*, *S. meliloti* and *Bradyrhizobium* spp. ORS285 and ORS278, the uptake is mediated by the SbmA/BacA/BclA transporters<sup>9, 17, 30</sup>. Thus strains expressing one of these transporters display a significantly increased sensitivity to bleomycin or Bac7. We find that wild type *B. diazoefficiens* strain USDA110 or the *E. coli* sbmA and *S. meliloti* bacA mutants expressing USDA110 *bclA* from the pRF771 plasmid are more sensitive to bleomycin or Bac7 than the corresponding strains lacking *bclA* (Fig. 1), in agreement with BclA of USDA110 being able to transport these peptides.

Contrary to bleomycin and Bac7, the sensitivity to antimicrobial NCR peptides is reduced in the presence of *sbmA*, *bacA* or *bclA*⁸,⁹. This opposite response is likely because the toxicity of the NCR peptides resides in their potential to provoke membrane permeability and loss of membrane potential, rather than in the inhibition of some intracellular process. Also the *bclA* gene of USDA110 is able to reduce the sensitivity to NCR peptides in the *S. meliloti* *bacA* mutant (Fig. 2a). In addition, similarly as shown before for the *S. meliloti* *bacA* or *Bradyrhizobium*
<figure><img src="image_3.png" /><figcaption>**Figure 1.** The USDA110 *bclA* gene confers sensitivity to the antibiotic bleomycin and the antimicrobial peptide Bac7. (a) Bleomycin sensitivity in *B. diazoefficiens* strain USDA110 and its *bclA* mutant derivative and in *E. coli* strain BW25113Δ*sbmA* and *S. meliloti* strain Sm1021Δ*bacA* carrying an empty vector or a vector expressing USDA110 *bclA*. Blue bars are for strains expressing *bclA* and red bars are for strains lacking the gene. Bleomycin concentrations were applied as indicated and growth was determined by optical density measurement at 600 nm with a plate reader after 72 h, 48 h or 24 h incubation for USDA110, Sm1021 and BW25113 respectively. (b) Survival of *E. coli* strain BW25113Δ*sbmA* (top) and *S. meliloti* strain Sm1021Δ*bacA* (bottom) derivatives carrying the empty pRF771 plasmid or the USDA110 *bclA* gene located on plasmid pRF771 after treatment with the peptide Bac7 at the indicated concentration. The surviving bacteria were counted and expressed as % from the water control treatment. Error bars in all panels are standard deviations.</figcaption></figure>

strain ORS285 *bclA* genes⁹, the expression of USDA110 *bclA* promotes the uptake of an FITC-modified NCR247 peptide into the *S. meliloti* strain Sm1021Δ*bacA* (Fig. 2b), while FITC alone is not taken up (Fig. 2c).

Finally, the *bclA* gene of USDA110 can complement the ORS285Δ*bclA* mutant for bacteroid differentiation in *A. indica* nodules (Fig. 3a–f). This mutant induces small nodules in which bacteria do not differentiate and die as revealed by live/dead staining of nodule sections⁹ (Fig. 3c,d). The USDA110 *bclA* gene, similarly as the ORS285 *bclA* gene, and introduced into this mutant on the pMG103 plasmid, restores the wild type phenotype with the formation of large nodules inhabited with well-formed spherical bacteroids that remain alive (Fig. 3a,b,e,f). Plants inoculated with both complemented strains showed vigorous growth, contrary to those inoculated with the un-complemented mutant indicating that nitrogen fixation defect of the mutant is restored by the USDA110 *bclA* gene. Moreover, the USDA110 *bclA* gene also complements in part the Sm1021Δ*bacA* strain for nodulation of *M. sativa* (Fig. 3g–l). The nodules induced by the complemented strain become elongated and pinkish, compared to the small white nodules induced by the *bacA* mutant (Fig. 3g,i,k). The *bacA* mutant bacteria die as
<figure><img src="image_4.png" /><figcaption>a</figcaption></figure>

Sm1021ΔbacA

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

**Figure 2.** The USDA110 *bclA* gene confers resistance to antimicrobial NCR peptides and mediates the uptake of FITC-NCR247. (a) *S. meliloti* strain Sm1021Δ*bacA* derivatives expressing no *bacA*-related gene (empty vector) or the USDA110 *bclA* gene were incubated with NCR247, NCR335 or NCR035 or with water (control) and the surviving bacteria were counted and expressed as % from the control treatment. Error bars in all panels are standard deviations. (b) FITC-NCR247 uptake by *S. meliloti* strain Sm1021Δ*bacA* or its derivative expressing the USDA110 *bclA* gene was measured by flow cytometry in the presence of trypan blue to quench extracellular fluorescence. FITC-positive cells are marked with a green box. (c) FITC is not taken up by *S. meliloti* derivatives. FITC-treated bacteria were analyzed by flow cytometry in the presence of trypan blue to quench extracellular fluorescence. No FITC-positive cells were detected.
<figure><img src="image_11.png" /><figcaption>**Figure 3.** The *bclA* gene of *B. diazoefficiens* USDA110 complements the *bclA* mutation in *Bradyrhizobium* sp. ORS285Δ*bclA* and the *bacA* mutation in *S. meliloti* Sm1021Δ*bacA*. (a,c,e) Phenotype of *A. indica* nodules at 14 dpi, infected with the indicated strains. Scale bars are 1 mm. (b,d,f) Bacteroid viability determined by live/dead and calcofluor white staining of nodule sections and confocal microscopy in *A. indica* nodules induced by the indicated strains. Scale bars are 10 µm. (g,i,k) Phenotype of *M. sativa* nodules at 28 dpi, infected with the indicated strains. Scale bars are 1 mm. (h,j,l) Bacteroid viability determined by live/dead staining of nodule sections and confocal microscopy in *M. sativa* nodules induced by the indicated strains. Scale bars are 10 µm.</figcaption></figure>

soon as they are released inside the nodule cells and exposed to the NCR peptides⁸ (Fig. 3j). The bacteroids of the complemented strain however are viable within the symbiotic cells indicating that the hypersensitivity of the *bacA* mutant to the NCR peptides is suppressed by the *bclA* gene of USDA110 (Fig. 3l). Nevertheless, the defect of the *S. meliloti* *bacA* mutation is only partially repaired by *bclA* of USDA110 because nodules do not fix nitrogen
and do not support plant growth (data not shown). This phenotype is similar to the one obtained with the S. *meliloti* *bacA* complementation by *bclA* of strain ORS285⁹ or even by more similar *bacA* genes of *Sinorhizobium* or *Rhizobium* species³¹. This suggests that although rhizobial *bacA* and *bclA* genes have overlapping specificity for peptide uptake, they may also display differences in the set and/or amount of NCR peptides they can handle, probably because they evolved in the context of specific interactions with host plants, each producing its specific arsenal of NCR peptides. Nevertheless, the USDA110 *bclA* gene seems to be capable to treat the *Aeschynomene* NCR peptides.

Together, the bleomycin, Bac7 and NCR peptide sensitivity assays as well as the *in vivo* complementation of the *Bradyrhizobium* sp. ORS285 *bclA* mutation or the S. *meliloti* *bacA* mutation indicate that the *bclA* gene of *B. diazoefficiens* strain USDA110 is functional and has a similar activity as the *bclA* gene of *Bradyrhizobium* strain ORS285 and the S. *meliloti* *bacA* gene.

**Bradyrhizobium diazoefficiens USDA110 BclA is not required for symbiosis with NCR-lacking soybean.** In agreement with the taxonomic position of soybean within the Millettioids, the *B. diazoefficiens* strain USDA110 bacteroids in soybean nodules are undifferentiated and display no cell enlargement and polyphloidy as revealed by microscopy of nodule sections (Fig. 4c) and flow cytometry analysis of purified bacteroids (Fig. 4e). We tested whether the *bclA* gene in strain USDA110 is required for symbiosis with soybean. The *bclA* mutant of USDA110 was undistinguishable from the wild type strain for all parameters analyzed, including nodule tissue structure and bacterial occupation, bacterial viability, morphology, size and DNA content as well as nitrogen fixation (Fig. 4a-e). Thus the *Bradyrhizobium bclA* gene, similarly to *bacA* in other rhizobium species, is not required for symbiosis when bacteroids are not constrained by the host plant-produced NCR peptides to differentiate into an elongated and polyploid morphotype.

**Bradyrhizobium diazoefficiens USDA110 BclA is also not required for symbiosis with NCR-producing Aeschynomene afraspera.** *B. diazoefficiens* strain USDA110, which is a natural soybean symbiont, can also form functional nodules on *A. afraspera*²⁹,³². The strain forms nodules composed of a central zone with fully infected symbiotic cells and a cortex layer surrounding the infected cells. This nodule organization is very similar to the histology of nodules induced by the natural *Aeschynomene* symbiont *Bradyrhizobium* strain ORS285 (Fig. S1). Nevertheless, USDA110 is a less efficient nitrogen fixer, supporting lower plant biomass production and nitrogen accumulation than ORS285³² (Fig. S2). The latter strain forms elongated and polyploid bacteroids on *A. afraspera* and requires the *bclA* gene for elongated bacteroid formation⁹. Therefore, we analyzed the bacteroid type formed by strain USDA110 within *A. afraspera* nodules and the role of the USDA110 *bclA* gene. Unexpectedly, observations by confocal microscopy showed that USDA110 bacteroids in *A. afraspera* nodules were not or only very slightly elongated (Fig. 5b), contrary to bacteroids of strain ORS285 which are strongly elongated⁴,⁹. To confirm this unpredicted observation we used flow cytometry analysis of the bacteroid cell size, determined by the forward scatter (a measure for cell size), and the DNA content, measured by DAPI fluorescence. A slight increase in size of the bacteroids compared to the bacteria in culture was measured but this was not accompanied with an increase in the DNA content of the bacteroids (Fig. 5c). Thus the USDA110 bacteroids are much less or hardly differentiated compared to ORS285 bacteroids which have, besides the strong cell enlargement, also a marked increase in DNA content⁴,⁹. The absence of a pronounced differentiation of USDA110 bacteroids is not likely resulting from a defect in NCR gene expression in USDA110-infected nodules since five tested NCR genes were expressed at similar or even higher levels in USDA110-infected nodules compared to ORS285-infected nodules (Fig. S3).

In agreement with the absence of differentiation of the wild type USDA110, the USDA110$\Delta bclA$ mutant was not affected in symbiosis with *A. afraspera*: nodules infected with wild type or mutant were similar, supported plant growth and fixed nitrogen to the same extent (Fig. 5a,d), both types of nodules contained symbiotic cells which were completely infected with bacteroids that seemed not or only slightly elongated (Fig. 5b; Fig. S1) which was confirmed by flow cytometry (Fig. 5c).

**BclA is not required for the formation of differentiated bacteroids in *Aeschynomene afraspera* nodules by the USDA110 DD-CPase1 mutant.** Contrary to the wild type USDA110 strain, the UDSA110 DD-CPase1 mutant forms strongly elongated and polyploid bacteroids in the *A. afraspera* nodules²⁹ (Fig. 5b,c), indicating that by affecting their cell wall strength, the bacteria become sensitive to the NCR differentiation signals produced by the nodule cells. Nevertheless, the nitrogenase activity of plants infected with the UDSA110 DD-CPase1 mutant is strongly reduced, in large part because the mutant induces much less nodules than the wild type²⁹ (Fig. 5d). We created a *bclA/DD-CPase1* double mutant to determine whether the cell wall-determined bacteroid differentiation is depending on the BclA peptide transporter. Unexpectedly, the double mutant exhibited a similar symbiotic phenotype than the DD-CPase1 single mutant (Fig. 5b-d). This result indicates that the differentiation of USDA110, made possible by the inactivation of the DD-CPase1 gene, is independent of BclA.

A possible explanation could be that the DD-CPase1 mutation increases the permeability of cells for peptides, including NCR peptides, rendering the BclA transporter superfluous. To test this possibility, we measured sensitivity of strains against the peptide bleomycin which needs to be internalized to target the bacterial DNA. We found that the DD-CPase1 mutant strain displays a sensitivity to bleomycin similar to the wild type strain (Fig. 5e), indicating that the peptidoglycan remodeling, induced by the mutation, does not interfere with peptide uptake. Similarly, the double mutant is just as much or even slightly more resistant to bleomycin than the *bclA* mutant (Fig. 5e).
<figure><img src="image_12.png" /><figcaption>**Figure 4.** The *B. diazoefficiens* bclA mutant is not affected in symbiosis with soybean. (a) Plant growth and nodule phenotype of *G. max* plants inoculated with *B. diazoefficiens* USDA110 wild type and the bclA mutant at 14 dpi. Scale bars are 3.5 cm (top panels) and 1 mm (bottom panels). (b) The nitrogen fixation activity of *G. max* plants inoculated with *B. diazoefficiens* USDA110 wild type and the bclA mutant, measured by the acetylene reduction assay per mg of fresh nodule weight (nmol ethylene produced per hour incubation and per mg nodule weight). Box-plots represent in the rectangle the first quartile to the third quartile, divided by the median value, whiskers above and below the box show the minimum and maximum measured values. (c,d) Toluidine blue stained thin sections (two top rows) and confocal microscopy of bacteroid viability determination by live/dead staining (bottom row) of *G. max* nodules induced by USDA110 or USDA110ΔbclA. Scale bars are 100 µm for the top panels and 10 µm for the other panels. (e) Flow cytometry analysis of DNA content by DAPI fluorescence (top) and of cell size by forward scatter (bottom) in free-living *B. diazoefficiens* USDA110 (pink) or bacteroids isolated from *G. max* nodules infected with USDA110 wild type (blue) or USDA110ΔbclA (orange). The forward scatter and DAPI fluorescence profiles are completely overlapping for the 3 samples.</figcaption></figure>

**BclA and DD-CPase1 in Bradyrhizobium strain ORS285 act independently in bacteroid differentiation.** To further explore the interdependence of BclA and DD-CPase1 in bacteroid differentiation, we created the double mutant also in *Bradyrhizobium* strain ORS285. This strain forms nodules on *A. afraespera* as well as on *A. indica* in which it differentiates into either elongated polyploid or spherical polyploid bacteroids respectively⁴. A *bclA* mutation in this strain blocks the differentiation process in both hosts⁹ (Fig. 6a) while a *DD-CPase* mutation induces hypertrophied bacteroids²⁹ (Fig. 6a). In a similar way as in USDA110, the ORS285 double mutant was still capable to induce strongly enlarged bacteroids in nodules of both *A. afraespera* and *A. indica* (Fig. 6a). In *A. indica*, the *bclA* mutation induced death of the bacteria as revealed by the red fluorescence in the live/dead staining procedure of nodule sections (Fig. 6a). Even if in the double mutant, many bacteroids were strongly enlarged (Fig. 6a), others remained undifferentiated. This may be related to bacterial death induced by the *bclA* mutation in such a way that many bacteria die before having the chance to enlarge. Thus, these results indicate that the *bclA* and *DD-CPase1* mutations have a cumulative effect and that the *bclA* function is not upstream of the bacterial differentiation provoked by the *DD-CPase1* mutation. The cumulative effect of the two mutations is also observed when measuring with the acetylene reduction assay the nitrogenase activity of *A. afraespera* nodules (Fig. 6b). The *bclA* mutation has a stronger impact on nitrogen fixation than the *DD-CPase1*
<figure><img src="image_13.png" /><figcaption>**Figure 5.** The *B. diazoefficiens* bclA, DD-CPase1, bclA/DD-CPase1 mutant phenotypes in symbiosis with *A. afraspera*. (a) Plant growth and nodule phenotype of *A. afraspera* plants inoculated with *B. diazoefficiens* USDA110 wild type and indicated mutants at 14 dpi. Scale bars are 2 cm (top panels) and 1 mm (bottom panels). (b) Confocal microscopy of bacteroid viability determination by live/dead and calcofluor white staining of *A. afraspera* nodules induced by USDA110 wild type and indicated mutants. Scale bars are 10 µm. (c) Flow cytometry analysis of DNA content by DAPI fluorescence and of cell size by forward scatter (FS) in free living *B. diazoefficiens* USDA110 (blue) or bacteroids isolated from *A. afraspera* nodules infected with USDA110 wild type and indicated mutants (pink). (d) Relative nitrogen fixation activity (% of maximum/h/plant) of *A. afraspera* plants inoculated with *B. diazoefficiens* USDA110 wild type and indicated mutants, measured by the acetylene reduction assay per plant. Box-plots represent in the rectangle the first quartile to the third quartile, divided by the median value, whiskers above and below the box show the minimum and maximum measured values. (e) Bleomycin sensitivity in *B. diazoefficiens* strain USDA110 and its bclA, DD-CPase1, bclA/DD-CPase1 mutant derivatives at 0,0375 µg bleomycin/ml. Growth was determined by optical density measurement at 600 nm with a plate reader after 72 h. Error bars in panels (d) and (e) are standard deviations.</figcaption></figure>

mutation and the double mutant has the same low nitrogenase activity as the bclA single mutant while on the other hand bacteroids of the double mutant resemble morphologically more the DD-CPase1 bacteroids.

Discussion

The BacA transporter is widespread but not universal in bacteria. It has been recruited by rhizobia that infect host plants producing NCR peptides in their infected nodule cells. In such host cells, the BacA protein is essential for the rhizobia to differentiate and establish a chronic infection. However, BacA is not critical for symbiosis in host cells of legumes that do not produce NCR peptides⁸. Rhizobia such as *Bradyrhizobium* spp., which lack a *bona fide* *bacA* gene but infect legume hosts that do produce NCR peptides in the infected nodule cells, use the different but related transporter BclA⁹. Here we show that the BclA function in *B. diazoefficiens*, in analogy with the BacA function in other rhizobia, is superfluous for chronic infection of soybean symbiotic cells that do not produce NCR peptides or when the bacterial strain is inherently non-responsive to the peptides produced by *A. afraspera* symbiotic cells. In agreement with our current findings, we previously demonstrated that the *bacA* gene of the broad host range bacterium *Sinorhizobium fredii* strain HH103 is not required for symbiosis with soybean²⁶. These results, together with previously published studies support the view that the major, if not sole symbiotic role of SbmA_BacA domain transporters is the import of NCR peptides.

*B. diazoefficiens* strain USDA110 is not a natural symbiont of *Aeschynomene* host plants. The strain produces the same major Nod factors as the genuine *Aeschynomene* symbiont *Bradyrhizobium* strain ORS285 and is able to induce with equal efficiency nodules on *A. afraspera* roots which have a normal histological organization but a reduced nitrogen fixation and plant growth promoting activity³² (Fig. S2). However a striking difference between the two strains is the absence of marked differentiation of the USDA110 bacteroids (this work) while the ORS285
<figure><img src="image_14.png" /><figcaption>**Figure 6.** The *Bradyrhizobium* strain ORS285 Δ*bclA*, DD-CPase1, Δ*bclA/DDS-CPase1* mutant phenotypes in symbiosis with *A. afraspera* and *A. indica*. (a) Confocal microscopy of bacteroid viability determination by live/dead staining of *A. afraspera* and *A. indica* nodules induced by ORS285 wild type and the indicated mutants. Scale bars are 10 µm. (b) Relative nitrogen fixation activity (% of maximum/h/plant) of *A. afraspera* and *A. indica* plants inoculated with ORS285 wild type and the indicated mutants, measured by the acetylene reduction assay per plant. Box-plots represent in the rectangle the first quartile to the third quartile, divided by the median value, whiskers above and below the box show the minimum and maximum measured values.</figcaption></figure>

bacteroids strongly elongate and become polyploid⁴. The absence of differentiation of *B. diazoefficiens* strain USDA110 in *A. afraspera* nodules is unlikely to be the result of its *bclA* gene being a defective allele because we found that the USDA110 *bclA* mutant has a similar response to bleomycin as the corresponding mutants in strains ORS285 and ORS278. The USDA110 *bclA* gene fully complements the symbiotic defect of the ORS285 *bclA* mutation for bacteroid differentiation in *Aeschynomene*. It also complements *S. meliloti* *bacA* and *E. coli* *sbmA* mutants, similar to the corresponding genes of ORS285 and ORS278⁹, for *Bac*7- and NCR-uptake and -sensitivity and for symbiosis with *M. sativa*. Together, these findings demonstrate that USDA110 *BclA* has the same capacity
as BclA from ORS285 or BacA from S. meliloti to transport a broad range of antimicrobial peptides, including Bac7, bleomycin and NCR peptides. They also tell us that USDA110 BclA does not lack the ability to uptake some specific peptide or set of peptides produced in *Aeschynomene* nodules which would explain the inability of the strain to differentiate within these nodules. Moreover, transcriptome analysis by RNA-seq (data not shown) revealed that *bclA* of USDA110 is weakly expressed, in culture as well as in soybean or *A. afraspera* nodules, but the orthologous gene of ORS285 is similarly weakly expressed. It is also unlikely that the lack of differentiation of USDA110 bacteroids results from a defect in NCR gene expression in USDA110-infected nodules since five tested NCR genes were expressed at similar or even higher levels in USDA110-infected nodules compared to ORS285-infected nodules.

A candidate factor for the strain USDA110 recalcitrance to bacteroid differentiation in *A. afraspera* nodules is the production of hopanoids³³. Hopanoids are a bacterial class of lipids, produced by certain bacteria, notably the bradyrhizobia but not the other rhizobal genera. These lipids, which are similar to eukaryotic steroid lipids, render membranes more rigid and resistant to membrane stresses, including to antimicrobial peptides³⁴, ³⁵. Strain USDA110 is remarkable among the bradyrhizobia because of its very high accumulation of hopanoids, up to 40% of the total cellular lipid fraction³³. The amount of hopanoids accumulating in the membranes of strain ORS285 is currently unknown but it could be sufficiently low for NCR responsiveness while the high hopanoid content in USDA110 could make this strain insensitive to the NCR signals of the host. In agreement with this hypothesis, an *hpnH* mutant of USDA110, affected in the synthesis of a specific type of hopanoids, the C35 hopanoids, is drastically impaired in symbiosis with *A. afraspera* but not with soybean; the mutant, in contrast to the wild type, is also sensitive to the antimicrobial activity of NCR peptides³⁶. Moreover, even if the *Bradyrhizobium* USDA110 and ORS285 strains have a high proportion of shared genes, they have nevertheless also a substantial number of unique genes and they belong to phylogenetically distinct clades in the *Bradyrhizobium* genus³⁷, ³⁸. A recent comparative genome analysis of ten *Bradyrhizobium* strains, including strains USDA110 and ORS285, revealed that USDA110 has over 3000 unique genes while ORS285 had almost 500 genes not shared with the other strains included in the study³⁸. Probably the different genetic repertoires of the USDA110 and ORS285 strains contribute to their differential response to the bacteroid differentiation factors produced by the symbiotic host cells. This difference could be exploited to identify new bacterial factors involved in the differentiation process and the response to NCR peptides. Whatever its molecular basis, the weak or non-responsiveness of strain USDA110 to the host-produced bacteroid differentiation factors most likely reflects the fact that this strain has co-evolved with NCR-lacking soybean and not with NCR-producing *Aeschynomene* or other NCR-producing hosts.

It is interesting to note that although the root nodules infected with the strain USDA110 fix nitrogen, they do so less efficiently when compared to nodules infected with strain ORS285; they supported less plant growth with the accumulation of lower amounts of nitrogen in the plant biomass, indicating that the undifferentiated bacteroids of USDA110 are less efficient than the differentiated bacteroids of ORS285. We recently demonstrated an analogous correlation between the level of bacteroid differentiation and symbiotic performance in the interactions of two *M. truncatula* accessions with four different *Sinorhizobium* strains³⁹. These concordant observations in different legumes suggest that the differentiation of bacteroids may contribute - by a mechanism that still has to be clarified - to the nitrogen fixation efficiency of rhizobium-legume symbioses⁴⁰⁻⁴². Variability in the response to NCR peptides can be at the basis of such variations in nitrogen fixation efficiency of *Rhizobium* strains. This has been nicely illustrated recently in the *S. meliloti*-M. truncatula interaction⁴³, ⁴⁴. Allelic variations in single NCR genes between *M. truncatula* accessions can determine the nitrogen fixation effectiveness of a particular *S. meliloti* strain while not affecting other strains. Thus the genome of rhizobial symbionts should be optimally aligned with the spectrum of NCR peptides produced by the nodule cells for efficient bacteroid differentiation and nitrogen fixation.

The BclA/BacA transporters provide protection against the antimicrobial activity of NCR peptides, both *in vitro* and *in vivo*, and in addition, the transporters are also required for bacterial elongation and polyploidization⁸, ⁹ (this work). This implies two possibilities concerning the molecular function of these transporters: either their sole function is providing protection against the antimicrobial stress in the nodule cells and the lack of bacteroid differentiation in the corresponding mutants is an indirect consequence of the reduced fitness of these mutants; alternatively, the transporters are, in addition to providing protection, also directly involved in the activation of the differentiation mechanism leading to cell enlargement and polyploidization, for example by transporting the NCR peptide signals into the bacterial cytoplasm and so allowing the peptides to interfere with the bacterial cell cycle machinery and other intracellular targets⁶, ⁷. Our phenotypic analyses of the *bclA/DD-CPase1* double mutants in the *Bradyrhizobium* strains USDA110 and ORS285 are in favor of the first option; they suggest that the sole role of BclA is the protection against the NCR stress. The abnormal and excessive differentiation provoked by the *DD-CPase1* mutation uncovered the differentiation potential of the *bclA* mutants. This interpretation of our data is in agreement with our previous observation that *in vitro* the *bacA* mutant of *S. meliloti* responds similarly as the wild type strain to NCR treatment by elongation and DNA amplification at peptide concentrations below their toxic level⁸. This implies that the NCR peptides are perceived at the outer or inner membrane, activating a signaling process that induces the differentiation process or alternatively, that NCR peptides can enter in a BclA/BacA independent manner the bacterial cells to reach their cytosolic targets for differentiation.

## Methods

**Bacterial growth conditions.** *Bradyrhizobium* spp. strains were grown in yeast mannitol⁴⁵ (YM) or buffered nodulation medium for bacteria (BNM-B)⁹ at 30 °C, *E. coli* and *S. meliloti* strains in Luria-Bertani (LB) medium at 37 °C or 30 °C respectively. Antibiotics for strain and plasmid selection were used at the following concentrations: streptomycin (500 µg/ml); spectinomycin (100 µg/ml); chloramphenicol (12.5 µg/ml); gentamycin (15 or 50 µg/ml); carbenicilin (50 µg/ml); tetracycline (10 µg/ml); kanamycin (50 or 200 µg/ml); neomycin (120 µg/ml); and nalidixic acid (25 µg/ml).
**Construction of bacterial mutants and complemented strains.** Standard molecular biology techniques were used for all cloning work. The construction of the *B. diazoefficiens* USDA110 *bclA* deletion mutant (USDA110Δ*bclA* ) was performed with the same strategy as the previously described mutant in *Bradyrhizobium* strain ORS285⁹. The following primers were used for the amplification of the upstream and downstream region of the *bclA* gene: GATAGAAAGCTTAAGCGTCCGGTGGTCACCGTCACCTG (forward, upstream region); GTCTTGCGCCGGATCCAGAGCTCTGTCTCCTGAGGGGATG (reverse, upstream region); ACAGAGCTCTGGATCCGGCGCAAGACGATCGCTATCGTAG (forward, downstream region); GATAGAACTAGTGTTCGGCACCTCGGACGCCTTCTAC (reverse, downstream region).

The mutant strains ORS285Δ*bclA*, ORS285Δ*DD-CPase1* and USDA110Δ*DD-CPase1* were obtained from previous studies⁹,²⁹. The double mutants ORS285Δ*bclADD-CPase1* and USDA110 Δ*bclADD-CPase1* were constructed from the corresponding single mutant in *bclA* and using the respective *DD-CPase1* pVO155-npt2-GFP mutagenesis plasmids which were introduced by triparental mating, as described²⁹.

The *bclA* gene of *B. diazoefficiens* USDA110 was amplified by PCR using the following primer couple: TTATCGTCTAGACCCTCAGGAGACAGAGCTCTGTGAAG (forward) and CATGATGGATCCGATCGTCTTGCGCCTCAGCGCGCCACGGTC (reverse). The PCR fragment was cloned into the *Xba*I and *Bam*HI sites of plasmid pRF771⁴⁶. The resulting plasmid was introduced into *S. meliloti* Sm1021Δ*bacA* by tri-parental conjugation and into *E. coli* BW25113Δ*sbmA* by electroporation. The same USDA110 *bclA*-encoding PCR fragment was also introduced in plasmid pMG103⁴⁷. First, the *trp* promoter was amplified with primers GTGCCGAATTCGGCAAATATACTGAAATAGGTGTTG (forward) and GAGTGCATGCGGTA CCGGATCCATGGAATCTAGATTTAAAGTACTTCGAA (reverse) and introduced into the *Eco*RI and *Sph*I sites of pMG103. The *bclA* PCR product was subsequently cloned in the *Xba*I and *Bam*HI sites of this pMG103 derivative. The resulting construct was introduced into *Bradyrhizobium* sp. ORS285Δ*bclA* by electroporation.

**Peptide sensitivity assays, single-cell NCR peptide uptake assay and flow cytometry.** All sensitivity assays were as described before⁹. The peptide uptake assay by flow cytometry was as described⁴⁸.

**Real-Time Quantitative PCR Expression Analyses.** The relative expression level of *NCR* genes was determined by RT-qPCR using methods for RNA extraction, cDNA synthesis, quantitative PCR, normalization and with primers previously described⁴.

**Plant growth and nodulation.** The soybean variety used in this study was *Glycine max* Williams 82. Seeds were cleaned with 100% ethanol for 30 sec and then surface-sterilized with bleach for 15 min. *A. indica* and *A. afraespera* seeds were surface sterilized using sulfuric acid (30 min) and bleach (30 min). Sterilized seeds were then transferred at 28 °C in the darkness on agar plates (0.8%) with tap water for germination. One-day old seedlings were transferred to test tubes containing liquid BNM⁴⁹. Seedlings were grown at 28 °C with a 16 h light regime and 70% humidity. Seven days after germination, each seedling was inoculated with 1 ml of a bacterial suspension adjusted to an OD₆₀₀ = 0.1.

**Symbiotic analyses.** Nitrogen-fixation assays were performed with ten plants per condition at 14 days post-inoculation, as previously described⁹. Plant biomass production was determined from the fresh weight of shoots that were removed from whole plants. The nitrogen and carbon content was determined on pools of three de-nodulated whole plants that were dried in an oven at 60 °C for 48 h. The dried samples were analyzed with an Isoprime element analyzer (Elementar).

For microscopy analysis, nodules were harvested, embedded with agarose (6%), and then freshly sectioned with a Leica VT1200S vibratome (Leica Microsystems GmbH) into 50 µm tissue slices. Slices were incubated for 20 min in Live/Dead BacLight (Molecular Probes) staining solution in Tris buffer (50 mM), pH 7.0, containing 0.01% CalcofluorWhiteM2R (Sigma). Sections were washed to remove excess of dye and observed using a Leica TCS SP8x confocal microscope.

Toluidine blue staining of thin sections was made as follows. Nodules were fixed in 1% glutaraldehyde/4% paraformaldehyde. After washing, nodules were dehydrated in ethanol series and embedded in Technovit 7100 resin (KulzerHistoTechnik) according to the manufacturer's instructions. Five µm sections were obtained with a Leica RM2155 microtome and stained with toluidine blue (0.005%). Bright field microscopy was performed with an Eclipse 80i microscope (Nikon).

**Bacteroid purification and flow cytometry analysis.** Bacteroids were purified as described²⁴. Bacteria were then fixed by heat treatment (70 °C, 10 min) and stained by DAPI (20 µg/ml) before flow cytometry analysis using a MoFlo ASTRIOS flow cytometer (Beckman Coulter). The bacterial DNA content was assessed by the DAPI fluorescence with a 355-nm laser line. Each single event was recorded and analyzed with the Summit 6.2 software (Beckman Coulter).

**Data Availability.** All data generated or analyzed during this study are included in this published article (and its Supplementary Information files).

**References**

1. Kondorosi, E., Mergaert, P. & Kereszt, A. A paradigm for endosymbiotic life: cell differentiation of *Rhizobium* bacteria provoked by host plant factors. *Annu. Rev. Microbiol.* **67**, 611–28 (2013).

2. Alunni, B. & Gourion, B. Terminal bacteroid differentiation in the legume–rhizobium symbiosis: nodule-specific cysteine-rich peptides and beyond. *New Phytol.* **211**, 411–417 (2016).

3. Van de Velde, W. et al. Plant peptides govern terminal differentiation of bacteria in symbiosis. *Science* **327**, 1122–1126 (2010).
4. Czernic, P. et al. Convergent evolution of endosymbiont differentiation in Dalbergioid and Inverted Repeat-Lacking Clade legumes mediated by nodule-specific cysteine-rich peptides. *Plant Physiol.* **169**, 1254-1265 (2015).

5. Montiel, J. et al. Morphotype of bacteroids in different legumes correlates with the number and type of symbiotic NCR peptides. *Proc. Natl. Acad. Sci. USA* **114**, 5041–5046 (2017).

6. Farkas, A. et al. *Medicago truncatula* symbiotic peptide NCR247 contributes to bacteroid differentiation through multiple mechanisms. *Proc. Natl. Acad. Sci. USA* **111**, 5183–5188 (2014).

7. Penterman, J. et al. Host plant peptides elicit a transcriptional response to control the *Sinorhizobium meliloti* cell cycle during symbiosis. *Proc. Natl. Acad. Sci. USA* **111**, 3561–3566 (2014).

8. Haag, A. F. et al. Protection of *Sinorhizobium* against host cysteine-rich antimicrobial peptides is critical for symbiosis. *PLoS Biol.* **9**, e1001169, doi:10.1371/journal.pbio.1001169 (2011).

9. Guefrachi, I. et al. *Bradyrhizobium* BclA is a peptide transporter required for bacterial differentiation in symbiosis with *Aeschynomene* legumes. *Mol. Plant. Microbe. Interact.* **28**, 1155–1166 (2015).

10. Runti, G. et al. Functional characterization of SbmA, a bacterial inner membrane transporter required for importing the antimicrobial peptide Bac7(1-35). *J. Bacteriol.* **195**, 5343-5351 (2013).

11. Laviña, M., Pugsley, A. P. & Moreno, F. Identification, mapping, cloning and characterization of a gene (sbmA) required for microcin B17 action on *Escherichia coli* K12. *J. Gen. Microbiol.* **132**, 1685-1693 (1986).

12. Yorgey, P. et al. Posttranslational modifications in microcin B17 define an additional class of DNA gyrase inhibitor. *Proc. Natl. Acad. Sci. USA* **91**, 4519–4523 (1994).

13. Salomón, R. A. & Farias, R. N. The peptide antibiotic microcin 25 is imported through the TonB pathway and the SbmA protein. *J. Bacteriol.* **177**, 3323–3325 (1995).

14. Mattiuzzo, M. et al. Role of the *Escherichia coli* SbmA in the antimicrobial activity of proline-rich peptides. *Mol. Microbiol.* **66**, 151–163 (2007).

15. Pränting, M., Negrea, A., Rhen, M. & Andersson, D. I. Mechanism and fitness costs of PR-39 resistance in *Salmonella enterica* serovar typhimurium LT2. *Antimicrob. Agents Chemother.* **52**, 2734–2741 (2008).

16. Marlow, V. L. et al. Essential role for the BacA protein in the uptake of a truncated eukaryotic peptide in *Sinorhizobium meliloti*. *J. Bacteriol.* **191**, 1519–1527 (2009).

17. Wehmeier, S. et al. Internalization of a thiazole-modified peptide in *Sinorhizobium meliloti* occurs by BacA-dependent and -independent mechanisms. *Microbiology* **156**, 2702–2713 (2010).

18. Ghosal, A., Vitali, A., Stach, J. E. M. & Nielsen, P. E. Role of SbmA in the uptake of Peptide Nucleic Acid (PNA) -peptide conjugates in *E. coli*. *ACS Chem. Biol.* **8**, 360–367 (2012).

19. Puckett, S. E. et al. Bacterial resistance to antisense peptide phosphorodiamidate morpholino oligomers. *Antimicrob. Agents Chemother.* **56**, 6147–6153 (2012).

20. Arnold, M. F. F. et al. Partial complementation of *Sinorhizobium meliloti bacA* mutant phenotypes by the *Mycobacterium tuberculosis* BacA protein. *J. Bacteriol.* **195**, 389–398 (2013).

21. Narayanan, S. et al. Mechanism of *Escherichia coli* resistance to pyrrhocoricin. *Antimicrob. Agents Chemother.* **58**, 2754–2762 (2014).

22. Guefrachi, I., Verly, C., Kondorosi, E., Alunni, B., & Mergaert, P. Role of the bacterial BacA ABC-transporter in chronic infection of nodule cells by *Rhizobium* bacteria. In *Biological Nitrogen Fixation* (ed. de Bruijn, F. J.) 315-324 (John Wiley & Sons, Inc., 2015).

23. Mergaert, P. 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).

24. Mergaert, P. et al. Eukaryotic control on bacterial cell cycle and differentiation in the *Rhizobium*-legume symbiosis. *Proc. Natl. Acad. Sci. USA* **103**, 5230–5235 (2006).

25. Maruya, J. & Saeki, K. The *bacA* Gene Homolog, *mlr7400*, in *Mesorhizobium loti* MAFF303099 is dispensable for symbiosis with *Lotus japonicus* but partially capable of supporting the symbiotic function of *bacA* in *Sinorhizobium meliloti*. *Plant Cell Physiol.* **51**, 1443–1452 (2010).

26. Crespo-Rivas, J. C. et al. *Sinorhizobium fredii* HH103 bacteroids are not terminally differentiated and show altered O-antigen in nodules of the Inverted Repeat-Lacking Clade legume *Glycyrrhiza uralensis*. *Environ. Microbiol.* **18**, 2392–2404 (2016).

27. Ardissoone, S. et al. Role of BacA in lipopolysaccharide synthesis, peptide transport, and nodulation by *Rhizobium* sp. strain NGR234. *J. Bacteriol.* **193**, 2218–2228 (2011).

28. Karunakaran, R. et al. BacA is essential for bacteroid development in nodules of galeoid, but not phaseoloid, legumes. *J. Bacteriol.* **192**, 2920–2928 (2010).

29. Gully, D. et al. A peptidoglycan-remodeling enzyme is critical for bacteroid differentiation in *Bradyrhizobium* during legume symbiosis. *Mol. Plant. Microbe. Interact.* **29**, 447–457 (2016).

30. Mardirossian, M. et al. The host antimicrobial peptide Bac7<sub>1-35</sub> binds to bacterial ribosomal proteins and inhibits protein synthesis. *Chem. Biol.* **21**, 1639–1647 (2014).

31. diCenzo, G. C., Zamani, M., Ludwig, H. N. & Finan, T. M. Heterologous complementation reveals a specialized activity for BacA in the *Medicago - Sinorhizobium meliloti* symbiosis. *Mol. Plant Microbe Interact.* **30**, 312–324 (2017).

32. Renier, A. et al. Photosynthetic Bradyrhizobium sp. strain ORS285 synthesizes 2-O-methylfucosylated lipichitoooligosaccharides for nod gene-dependent interaction with *Aeschynomene* plants. *Mol. Plant Microbe Interact.* **24**, 1440–1447 (2011).

33. Kannenberg, E. L., Perzl, M. & Härtner, T. The occurrence of hopanoid lipids in *Bradyrhizobium* bacteria. *FEMS Microbiol. Lett.* **127**, 255–261 (1995).

34. Schmerk, C. L., Bernards, M. A. & Valvano, M. A. Hopanoid production is required for low-pH tolerance, antimicrobial resistance, and motility in *Burkholderia cenocepacia*. *J. Bacteriol.* **193**, 6712–6723 (2011).

35. Malott, R. J., Steen-Kinnaird, B. R., Lee, T. D. & Speert, D. P. Identification of hopanoid biosynthesis genes involved in polymyxin resistance in *Burkholderia multivorans*. *Antimicrob. Agents Chemother.* **56**, 464–471 (2012).

36. Kulkarni, G. et al. Specific hopanoid classes differentially affect free-living and symbiotic states of *Bradyrhizobium diazoefficiens*. *MBio* **6**, e01251–15, doi:10.1128/mBio.01251-15 (2015).

37. Giraud, E. et al. Legumes symbioses: absence of nod genes in photosynthetic bradyrhizobia. *Science* **316**, 1307–1312 (2007).

38. Mornico, D. et al. Comparative genomics of *Aeschynomene* symbionts: insights into the ecological lifestyle of nod-independent photosynthetic bradyrhizobia. *Genes* **3**, 35–61, doi:10.3390/genes3010035 (2012).

39. Kazmierczak, T. et al. Specific host-responsive associations between *Medicago truncatula* accessions and *Sinorhizobium* strains. *Mol. Plant Microbe Interact.* **30**, 399–409 (2017).

40. Oono, R. & Denison, R. F. Comparing symbiotic efficiency between swollen versus nonswollen rhizobial bacteroids. *Plant Physiol.* **154**, 1541–1548 (2010).

41. Kereszt, A., Mergaert, P. & Kondorosi, E. Bacteroid development in legume nodules: evolution of mutual benefit or of sacrificial victims? *Mol. Plant-Microbe Interact.* **24**, 1300–1309 (2011).

42. Sen, D., Weaver, R. W. & Bal, A. K. Structure and organization of effective peanut and cowpea root nodules induced by rhizobial strain 32H1. *J. Exp. Bot.* **37**, 356–363 (1986).

43. Yang, S. et al. Microsymbiont discrimination mediated by a host-secreted peptide in *Medicago truncatula*. *Proc. Natl. Acad. Sci. USA* **114**, 6848–6853 (2017).

44. Wang, Q. et al. Host-secreted antimicrobial peptide enforces symbiotic selectivity in *Medicago truncatula*. *Proc. Natl. Acad. Sci. USA* **114**, 6854–6859 (2017).
45. Giraud, E., Hannibal, L., Fardoux, J., Verméglio, A. & Dreyfus, B. Effect of *Bradyrhizobium* photosynthesis on stem nodulation of *Aeschynomene sensitiva*. *Proc. Natl. Acad. Sci. USA* **97**, 14795–14800 (2000).

46. Wells, D. H. & Long, S. R. The *Sinorhizobium meliloti* stringent response affects multiple aspects of symbiosis. *Mol Microbiol.* **43**, 1115–1127 (2002).

47. Inui, M., Roh, J. H., Zahn, K. & Yukawa, H. Sequence analysis of the cryptic plasmid pMG101 from *Rhodopseudomonas palustris* and construction of stable cloning vectors. *Appl Environ Microbiol.* **66**, 54–63 (2000).

48. Benincasa, M. et al. Single cell flow cytometry assay for peptide uptake by bacteria. *Bio-Protocol* **6**, e2038, doi:10.21769/BioProtoc.2038 (2016).

49. Ehrhardt, D. W. et al. Depolarization of alfalfa root hair membrane potential by *Rhizobium meliloti* Nod factors. **256**, 998–1000 (1992).

**Acknowledgements**

We thank Christien Merrifield (Institute for Integrative Biology of the Cell) for comments on the manuscript. Q.B. and F.L. were supported by PhD fellowships from the Université Paris-Sud and I.G. by PhD fellowships from the Tunisian Ministry of Higher Education and Scientific Research, the International Relations Office of Université Paris-Sud and the Agence Universitaire de la Francophonie (AUF). The present work has benefited from the core facilities of Imagerie-Gif (http://www.i2bc.paris-saclay.fr), member of IBiSA (http://www.ibisa.net), supported by 'France-BioImaging' (ANR-10-INBS-04-01), and the Labex 'Saclay Plant Science' (ANR-11-IDEX-0003-02). This work was funded by the Agence Nationale de la Recherche, grant n° ANR-13-BSV7-0013.

**Author Contributions**

Q.B., I.G., E.G. and P.M. conceived the experiments; Q.B., I.G., D.G., F.L., O.P., J.F. and C.C. performed the experiments; Q.B., I.G., T.T., B.A., E.G. and P.M. analyzed the data. Q.B. and P.M. wrote the paper and all authors reviewed and corrected the manuscript.

**Additional Information**

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-08830-0

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_15.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017