# Unearthing the Plant Growth-Promoting Traits of Bacillus megaterium RmBm31, an Endophytic Bacterium Isolated From Root Nodules of Retama monosperma

Malika Affaf Dahmani, Antoine Desrut, Bouziane Moumen, Julien Verdon, Lamia Mermouri, Mourad Kacem, Pierre Coutos-Thévenot, Meriem Kaid-Harche, Thierry Bergès, Cécile Vriet

► To cite this version:

Malika Affaf Dahmani, Antoine Desrut, Bouziane Moumen, Julien Verdon, Lamia Mermouri, et al. Unearthing the Plant Growth-Promoting Traits of Bacillus megaterium RmBm31, an Endophytic Bacterium Isolated From Root Nodules of Retama monosperma. Frontiers in Plant Science, 2020, 11 (27), pp.124. 10.3389/fpls.2020.00124. hal-05517079

HAL Id: hal-05517079

https://hal.science/hal-05517079v1

Submitted on 18 Feb 2026

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 CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Unearthing the Plant Growth-Promoting Traits of *Bacillus megaterium* RmBm31, an Endophytic Bacterium Isolated From Root Nodules of *Retama monosperma*

OPEN ACCESS

**Edited by:**
Essaid Ait Barka,
Université de Reims Champagne-
Ardenne, France

**Reviewed by:**
Claire Prigent-Combaret,
UMR5557 Ecologie Microbienne
(LEM), France
José López-Bucio,
Universidad Michoacana de San
Nicolás de Hidalgo, Mexico

*Correspondence:
Cécile Vriet
cecile.vriet@univ-poitiers.fr

†These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Plant Microbe Interactions,
a section of the journal
Frontiers in Plant Science

Received: 08 November 2019

Accepted: 28 January 2020

Published: 27 February 2020

**Citation:**
Dahmani MA, Desrut A, Moumen B,
Verdon J, Mermouri L, Kacem M,
Coutos-Thévenot P, Kaid-Harche M,
Bergès T and Vriet C (2020)
Unearthing the Plant Growth-
Promoting Traits of Bacillus
megaterium RmBm31, an Endophytic
Bacterium Isolated From Root
Nodules of Retama monosperma.
Front. Plant Sci. 11:124.
doi: 10.3389/fpls.2020.00124

Malika Affaf Dahmani<fnref n="1" />,²†, Antoine Desrut<fnref n="3" />†, Bouziane Moumen<fnref n="3" />, Julien Verdon<fnref n="3" />, Lamia Mermouri<fnref n="1" />, Mourad Kacem<fnref n="4" />, Pierre Coutos-Thévenot<fnref n="3" />, Meriem Kaid-Harche<fnref n="1" />, Thierry Bergès² and Cécile Vriet<fnref n="3" />*

Plants live in association with complex populations of microorganisms, including Plant Growth-Promoting Rhizobacteria (PGPR) that confer to plants an improved growth and enhanced stress tolerance. This large and diverse group includes endophytic bacteria that are able to colonize the internal tissues of plants. In the present study, we have isolated a nonrhizobial species from surface sterilized root nodules of *Retama monosperma*, a perennial leguminous species growing in poor and high salinity soils. Sequencing of its genome reveals this endophytic bacterium is a *Bacillus megaterium* strain (RmBm31) that possesses a wide range of genomic features linked to plant growth promotion. Furthermore, we show that RmBm31 is able to increase the biomass and positively modify the root architecture of seedlings of the model plant species *Arabidopsis thaliana* both in physical contact with its roots and via the production of volatile organic compounds. Lastly, we investigated the molecular mechanisms implicated in RmBm31 plant beneficial effects by carrying out a transcriptional analysis on a comprehensive set of phytohormone-responsive marker genes. Altogether, our results demonstrate that RmBm31 displays plant growth-promoting traits of potential interest for agricultural applications.

**Keywords:** plant growth-promoting rhizobacteria, endophyte, root nodules, *Bacillus megaterium*, *Retama monosperma*, *Arabidopsis thaliana*, volatile compounds

**Abbreviations:** IAA, Indole-3-Acetic Acid; PGPR, plant growth-promoting rhizobacteria; PGP, plant growth promotion; VOC, volatile organic compounds.

<footnote n="1">¹Laboratoire des Productions, Valorisation végétales et microbiennes (LP2VM), Département de biotechnologies, Faculté SNV, Université des Sciences et de la Technologie d'Oran-Mohammed BOUDIAF (USTO M.B), Oran, Algéria, ²Laboratoire Signalisation et Transports Ioniques Membranaires (STIM), CNRS EA7349, Université de Poitiers, Poitiers, France,</footnote>

<footnote n="3">³Laboratoire Ecologie et Biologie des Interactions (EBI), UMR CNRS 7267, Université de Poitiers, Poitiers, France,</footnote>

<footnote n="4">⁴Département de Biotechnologie, Faculté SNV, Université d'Oran Ahmed Ben Bella, Oran, Algéria</footnote>
INTRODUCTION

*Retama monosperma* (L.) Boiss. is a perennial leguminous species that plays an important role in stabilizing sand dunes and in soil fertility due to its developed root system and ability to form symbiosis with root-nodulating nitrogen-fixing rhizobia (Selami et al., 2014). Relatively recently, several nonrhizobial bacteria have also been isolated from root nodules of leguminous (*Fabaceae*) plant species (Muresu et al., 2008; Peix et al. 2014; Martínez-Hidalgo and Hirsch, 2017). The function and biological importance of these endophytic bacteria in nodules remain largely unknown (Martínez-Hidalgo and Hirsch, 2017).

Endophytic bacteria belong to a larger and diverse set of soil bacteria referred to as plant growth-promoting rhizobacteria (PGPR) that are able to confer to the host plant an improved growth and/or tolerance to diverse biotic and abiotic stresses (Lugtenberg and Kamilova, 2009; Vacheron et al., 2013; Hardoim et al., 2015; Santoyo et al., 2016). The exact mechanisms by which PGPR promote plant growth are not fully understood, but evidence in the literature suggests the existence of preferential associations between genes contributing to phytobeneficial traits (Bruto et al., 2014). Their plant beneficial functions involve the bacterial ability to improve plant nutrient acquisition, in particular *via* their capacity to carry out mineralization/solubilization of organic/insoluble phosphates, fix nitrogen, and to produce siderophores in iron-limiting conditions that chelate ferric ions and make them available to the plant (Hardoim et al., 2015). Some PGPR also stimulate the plant's iron (Fe) acquisition machinery and thereby improve its uptake of Fe (Zhang et al., 2009; Zhou et al., 2016). In addition, some PGPR are able to produce plant growth regulators such as auxins, gibberellins, and cytokinins and/or to change their concentration *in planta* (Vacheron et al., 2013). Indeed, synthesis of the phytohormone auxin (Indole-3 Acetic Acid; IAA) by bacteria has been well characterized (Spaepen and Vanderleyden, 2011). The precursor that modulates the level of IAA biosynthesis is the amino acid tryptophan found in plant root exudates (Lugtenberg and Kamilova, 2009; Vacheron et al., 2013). IAA controls many plant growth and development processes including root architecture traits such as the formation of lateral roots and root hairs and the primary root length. Furthermore, some PGPR strains release Volatile Organic Compounds (VOC) that are able to promote plant growth (Gutiérrez-Luna et al., 2010; Hernández-Calderón et al., 2018). Among the substances identified figure the volatile compounds 2,3-butanediol, acetoine, 1-hexanol, pentadecane, and indole (Ryu et al., 2003; Blom et al., 2011). Indole is known to have an effect on the root system architecture, beyond its role as auxin precursor, by interfering with the auxin signaling pathway (Bailly et al., 2014).

PGPR also possess the ability to promote plant growth indirectly, notably by improving the plant tolerance against diverse abiotic stresses or pathogens. This biocontrol activity is achieved through direct antagonism with soil-borne pathogens (i.e. via the production of antibiotics, toxins and/or antimicrobial enzyme activities) and/or by competing with them for the acquisition of nutrients (eg. via siderophore-mediated competition for iron) (Lugtenberg and Kamilova, 2009; Vacheron et al., 2013). Besides, some PGPR strains that are able to enter inside the roots and establish themselves as endophytes are also able to reach the vascular

system and migrate to the stem and leaves and thereby exert a direct biocontrol activity against leaf pathogens (Miotto-Vilanova et al., 2016). Moreover, some PGPR confer to plants an enhanced tolerance against a wide range of aboveground pathogens by stimulating the plant defenses. This process known as Induced Systemic Resistance (ISR) involves a priming phenomenon that allows PGPR-inoculated plants to react faster and/or more efficiently to subsequent pathogen attacks (Pieterse et al., 2014). Lastly, some PGPR possess the ability to promote plant growth by affecting the plant interaction with other beneficial microorganisms including N₂-fixing rhizobia and arbuscular mycorrhizal fungi (Barea et al., 2005). Future researches on the interaction between plants and PGPR are needed before these beneficial bacteria can be widely and efficiently exploited in agriculture as biocontrol and/or biofertilizer agents.

In the present study we have isolated an endophytic bacterium from surface-sterilized root nodules of *Retama monosperma*. This bacterium was identified as a *Bacillus megaterium* strain using whole genome sequencing and named RmBm31. RmBm31 harbors numerous genomic features related to plant growth promotion, hence we have explored its capacity to enhance plant growth and development using the model plant species *Arabidopsis thaliana*. Results from our *in vitro* assays reveal that RmBm31 is able to increase the biomass and improve the root system architecture of the seedlings both in physical contact with their roots and solely *via* the emission of volatile compounds. Lastly, we have explored the effects of RmBm31 on the expression of a comprehensive set of marker genes of the plant phytohormone pathways in order to shed light on the molecular mechanisms involved in its plant growth-promoting ability. Remarkably, we did not identify major RmBm31-triggered transcriptional changes in the pathways of the plant growth-related hormones auxin, cytokinin and gibberellins. Instead, we found that this PGPR induces the expression of genes involved in the signaling pathways of ethylene, jasmonic acid, and salicylic acid, three phytohormones commonly implicated in plant response to biotic stress. In agreement with previous studies (Bordiec et al., 2011; Drogue et al., 2014; Pieterse et al., 2014; Hernández-Calderón et al., 2018; Kandaswamy et al., 2019), these findings highlight the importance of plant defense responses in beneficial plant-PGPR interactions and suggest yet unknown molecular mechanisms may be at play in RmBm31-triggered growth-promoting effects.

MATERIALS AND METHODS

Isolation of *Bacillus megaterium* RmBm31
From Surface-Sterilized Root Nodules of
*Retama monosperma* (L.) Boiss

Root nodules were collected from *Retama monosperma* plants growing in natural conditions in the area of Oran, Northwest Algeria (N 35°42 W 00°34; 162 m altitude). The nodules were teared off from the roots, washed once with tap water, and then thoroughly washed with sterile, distilled water. Their surfaces were sterilized by immersion in 70% ethyl alcohol for 20 to 30 s
and then in 2.6% sodium hypochlorite solution for 5 min. The nodules were then washed 7 to 8 times with sterile distilled water. Next, the surface-sterilized nodules were crushed with sterile forceps and plated on YMA solid medium (mannitol 10 g.L⁻¹; yeast extract 0.12 g.L⁻¹; NaCl 1 g.L⁻¹; MgSO₄ 0.2 g.L⁻¹; K₂HPO₄ 0.5 g.L⁻¹; agar 20 g.L⁻¹; pH 6.8) and incubated at 28°C. Single colonies were picked up, purified through repetitive streaking on YMA medium and stored in 30-40% glycerol at -80°C for subsequent molecular identification.

# Production of Indole-3-Acetic Acid (IAA) by RmBm31

## Detection of IAA by Thin-Layer Chromatography

For IAA identification assay, 10 µl of extract and 5 µl of IAA standard solution (50 µg.ml⁻¹ in water/methanol) were spotted on a glass capillary plate coated with a thin layer of silica gel 60 (TLC F 254 plates; Merck, Germany). The solvent system used was N-butanol:ammonia:water (10:1:10) (Sethia et al., 2015). The plates were developed by spraying Salkowski reagent (0.01 M of FeCl₃ in 35% perchloric acid) and heating at 40°C for 10 min. Spots were visualized under light and UV light (254 nm) and identified by comparing Retention factor (Rf) value of the IAA produced by RmBm31 with the standard.

## IAA Extraction, Concentration and Identification by Liquid Chromatography-Mass Spectrometry

In order to confirm the production of IAA by RmBm31, 100 ml of YMA medium supplemented with L-tryptophan (0.2 mg.ml⁻¹) was inoculated with RmBm31 and incubated at 28°C with orbital shaking at 150 rpm. After 3 days of culture, bacterial cells were eliminated by centrifugation at 19,000 g for 15 min at 4°C. The pH of the supernatant was adjusted to pH 9 with 1 M KOH to keep IAA ionized and extracted with an equal volume of ethyl acetate. The lower aqueous phase was transferred to an Erlenmeyer and the pH of the solution was lowered below pH 3 to preserve IAA in protonated form. The acidic sample was extracted with an equal volume of ethyl acetate and the upper organic phase was evaporated under nitrogen. The dried extract was dissolved in 1 ml of methanol and kept at 4°C prior to be analyzed by LC-MS.

Settings were first established with IAA standard (1 mg.mL⁻¹) by Reverse Phase-HPLC analysis. Ten µl was injected into a reverse phase C18 column (5.25 µm x 0.46 cm) equipped with a binary pump and a UV detector. The sample was eluted with methanol and 0.3% acetic acid at 0.8 mL.min⁻¹ flow rate. For LC-MS, the extracts were first filtered through a 0.2 µm membrane and 10 µl was injected. The molecular mass of the samples was determined by Electrospray Ionization-Mass Spectrometry (ESI-MS) with a XevoQ-TOF (Waters, Milford, MA, USA) mass spectrometer coupled to a fused C18 HPLC column (5.25 m x 0.46 cm) and analyzed in a positive mode. The spray voltage was set to 3.0 kV, the source temperature to 120 °C, and the desolvation temperature to 250°C. The MS spectra of IAA were identified through retention time and mass spectra comparison with IAA standard.

# Sequencing of RmBm31 Genome and Analysis of Encoded Functions

## RmBm31 DNA Extraction

Glycerol stock of *Bacillus megaterium* RmBm31 strain was streaked on LB medium and grown 24 h at 28°C. One colony forming unit (CFU) was used to inoculate LB liquid medium (Tryptone 10 g.L⁻¹; Yeast extract 5g.L⁻¹; and NaCl 5g.L⁻¹) without amendments of antibiotic and incubated at 28°C under orbital agitation until the OD₆₀₀nm reached 0,7-1.

Total DNA was isolated and purified using NucleoSpin Microbial® DNA kit (MACHEREY-NAGEL). RNAse treatment was performed according to manufacturer instructions. DNA integrity was checked by 1.5% agarose gel electrophoresis and quantification was performed on a Qubit fluorometer using a Qubit DNA HS Assay Kit (Thermo Fischer Scientific). The absence of RNA contamination was checked using a Qubit RNA Assay Kit (Thermo Fischer Scientific).

## Identification of RmBm31 by DNA-Sequencing and Phylogenomic Analysis

The genome of *Bacillus megaterium* RmBm31 was sequenced by GATC Biotech using the Illumina HiSeq technology in paired-end mode. A total number of 24.393.100 reads were obtained representing a theoretical coverage of 500×.

Genome size was estimated by counting k-mer frequencies of quality-filtered paired-end reads with Jellyfish tool (version 2.2.6) (Marçais and Kingsford, 2011). De novo assembly of the genome was performed using the Shovill pipeline (version 0.2; https://github.com/tseemann/shovill). This pipeline includes all the necessary tools to assemble a bacterial genome: getting statistics for raw reads (seqtk), estimating genome size (kmc), trimming reads (Trimmomatic), correcting reads (lighter), overlapping reads (flash), assembling reads (SPAdes, version 3.9), and correcting errors from the assembly (BWA and samtools and pilon). Annotation of the assembled genome was done using the pipeline prokka (version 1.11) (Seemann, 2014). Lastly, classification of CDS into functional COG (Clusters of Orthologous Groups) categories was performed using WebMGA (http://weizhongli-lab.org/webMGA/) (Wu et al., 2011). The circular representation of the genome was drawn with cgview (Stothard and Wishart, 2005). Like the closely related *Bacillus* strains, RmBm31 might also possess plasmids (not investigated in detail in the present study).

RmBm31 was identified as a *Bacillus megaterium* strain by performing a phylogenomic analysis using a representative set of *Bacillus* genomes available in NCBI Genbank. Orthofinder2 (Emms and Kelly, 2015) was used to find their core genome constituted of 38 single copy genes. A phylogenomic tree was generated from the alignment of these genes using IQ-TREE (Nguyen et al., 2015) (maximum-likelihood method) with bootstrap test (1,000 replicates) (Hoang et al., 2018). As in a previous study (Bhattacharyya et al., 2017), the tree was rooted with *Sulfolobus acidocaldarius* DSM 639.
## Genomic Analysis of Plant Growth-Promoting Features

A comprehensive set of 105 genes coding for known PGP-related genes were selected from searches in the literature. Using the NCBI Protein and KEGG databases, their protein sequence IDs were retrieved preferentially from the *Bacillus megaterium* genome QM B1551, or in the other *Bacillus* species genomes *Bacillus cereus* ATCC 14579 and *Bacillus* subtilis QB928 (phylogenetic tree of Zou et al., 2010). Failing that, protein sequence IDs were retrieved from the gram-negative proteobacteria genomes *Burkholderia* sp. KJ006, *Pseudomonas fluorescens* F113, or *Pseudomonas aeruginosa* PAO1. Sequences of these proteins were then retrieved from NCBI Batch Entrez. Their homologs in RmBm31 genome were searched using Blastp. Only best hits with an E-value cutoff of 1e-05 (≤0.00001) were retained. Moreover, proteins were considered as homologs only if they shared the same Pfam functional protein domains. Protein domain assignment was performed using InterProScan.

## Inoculation of *Arabidopsis*
*Thaliana* Seedlings With *B. megaterium*
RmBm31 and In Vitro Cocultivation Assays
Plant Material and Growth Conditions

*Arabidopsis thaliana* (Arabidopsis) ecotype Columbia (Col-0) was used as model plant species is this study. Arabidopsis seeds were surface sterilized for 7 min in 5% (v/v) hypochlorite and 0.1% Tween®20 detergent solution and washed five times in sterile water. Seeds were sown in sterile Petri-dishes (140 × 40 × 20,6 mm) on half strength (0,5×) Murashige and Skoog (MS) medium (M0222, Duchefa Biochemie, Haarlem, The Netherlands), without sucrose, supplemented with 0.5% of MES (Morpholino-Ethane-Sulfonic acid monohydrate; MW = 213,2 gmol⁻¹) (ACROS Organics™, 172591000) and 0.8% (w/v) of plant Agar Type A (Sigma-Aldrich®, A4550), pH 5.8. The plates containing 10 seeds each were sealed with two layers of parafilm. After two days of seed stratification at 4°C in the dark, the plates were positioned vertically in a growth chamber under a long-day photoperiod (16 h of light), with a light intensity of 120–140 µmol.m⁻².s⁻¹, and at a temperature of 22°C. Seedlings were grown for five days in these growth conditions prior their cocultivation with the rhizobacteria.

## Preparation of RmBm31 Inoculum and Inoculation Treatment

For preparation of the inoculum, an aliquot of a glycerol stock of BmRm31 was streaked on solid LB medium (10 g.L⁻¹ Bacto-tryptone, 5 g.L⁻¹ Yeast Extract, 5 g.L⁻¹ NaCl, and 15 g.L⁻¹ Agar, pH7). After 24 h at 28°C, bacterial cells were collected in 10 mM MgSO₄, washed twice with 50 ml of 10 mM MgSO₄ by centrifugation for 5 min at 5,000 g, and resuspended in 50 ml of 10 mM MgSO₄. The bacterial titer was adjusted to an OD₆₀₀nm of 0.01 to obtain an inoculum with a bacterial density of 2.10⁶ colony-forming units.ml⁻¹ (CFU.ml⁻¹). For all experiments, this bacterial density was confirmed by counting the number of CFU on LB medium.

For experiments with bacteria in physical contact with plant roots, five days after seed sowing, Arabidopsis seedlings (10 per plates) were inoculated with 10 µl of the inoculum at 2.10⁶ CFU.ml⁻¹ in physical contact with their roots, at 1 cm under the shoot-root junction. For experiments involving the effects of RmBm31 volatile compounds only, 100 µl of bacterial suspension was spotted on LB medium into a small Petri dish (35 × 10mm) itself placed inside a larger Petri dish (140 × 140 × 20,6 mm) containing the 5-day-old Arabidopsis seedlings (10 per plate). In this experimental setting, the seedlings were physically separated from the inoculum, but gas exchange was allowed between the two compartments. For the mock treatment, seedlings were treated the same way with a 10 mM MgSO₄ solution. The plates were sealed with two layers of parafilm and transferred back to the plant growth chamber.

## Phenotypic Analyses

For shoot and root fresh weight measurements, seedlings were sectioned at the root/shoot junction, and the weight (pool of three plants for the roots) immediately measured on an analytical balance.

For primary and lateral root analyses, images of Arabidopsis seedlings were taken using a digital camera (D5200 Nikkon equipped with a AF 90 mm 1:2.8 Macro 1:1 TAMRON lens). For root hair analyses, images of the primary root tips were acquired using a macroscope (Macro View MVX10 Olympus) with a 20-fold magnification. The ImageJ software (Schneider et al., 2012) equipped with the SmartRoot plugin (Lobet et al., 2011) was used to measure the primary and lateral root lengths and root hair length and density on a 1 mm² segment located at 1 mm above the root tip.

Data are means ± the standard error of the mean (SEM) of nine biological replicates from three independent experiments.

## Gene Expression Profiling
Total RNA Extraction

Plant samples for gene expression analysis were harvested at midday (8 h of light, 16 h photoperiod), 7 days post inoculation. Roots and shoots of Arabidopsis seedlings were harvested separately by sectioning the root-shoot junction, immediately frozen in liquid nitrogen, and stored at -80°C.

Total RNA was extracted from 25–100 mg of shoot and root tissues using a phenol/chloroform extraction procedure adapted from Box et al. (2011). Briefly, frozen tissues were grinded using a TissueLyser II (Qiagen) bead mill (30 Hz; 30 s). Total RNA was extracted with 0.5 ml of RNA extraction buffer solution (25 mM Tris-HCl (pH 8); 25 mM EDTA; 75 mM NaCl; 2% (w/v) SDS; 7.8% (v/v) β-Mercaptoethanol) and 0.5 ml of Phenol/Chloroform/Isoamyl Alcohol solution (25:24:1, v/v/v). Samples were vortexed for 30 s and centrifuged 15 min at 13,200 g at 4°C. The aqueous upper phase was recovered, transferred to new tubes, and an equal volume (0.4 ml) of Chloroform/Isoamyl Alcohol (24:1; v/v) was added. Samples were then centrifuged 5 min at 13,200 g at 4°C and the aqueous layer was transferred to a new tube containing 0.4 ml of 4 M LiCl₂. Samples were precipitated overnight at 4°C and centrifuged 20 min at 20,600 g at 4°C.
The supernatant was discarded, and the RNA pellet was resuspended in 50 µl of nuclease-free water (Promega). After resuspension, 5 µl of sodium-acetate (3 M) and 125 µl of ice-cold 96% (v/v) ethanol were added. Samples were then left 1 h at -20°C for RNA precipitation. After centrifugation 20 min at 20,600 g at 4°C, the supernatant was discarded and the pellet washed in ice-cold 70% (v/v) ethanol by centrifugation 5 min at 20,600 g at 4°C and dried under a fume hood. The RNA pellet was finally resuspended with 6-30 µl of nuclease-free water and stored at -80°C. RNA was quantified by spectrometry at 260 nm using the Multiskan GO Thermo Scientific spectrophotometer equipped with a µDrop™ Thermo Scientific Plate. RNA integrity was checked by gel electrophoresis on an agarose gel and by a measure of the 260/280 nm ratio.

## Relative Gene Expression Analysis by Real Time, Quantitative RT-PCR

Primers for qRT-PCR were designed using the NCBI Primer-Blast software (Ye et al., 2012), with the following criteria: a primer size comprised between 18 and 25 bp, a GC % of 45-60%, a melting temperature (T<sub>m</sub>) between 58 and 63°C, and a PCR product size of 50 to 200 bp. Moreover, preferences were given to primer pairs that were exon-exon shuffling or intron spanning.

Gene expression analyses were performed by qRT-PCR using the GoTaq qPCR MasterMix (Promega) according to the manufacturer instructions (1× GoTaq® qPCR Master Mix, 0.33 µM of forward and reverse primer, and 5 µl of 10-fold diluted cDNA per well). The amplification reaction was carried out using a 96-well plate thermal cycler (Mastercycler Realplex2, Eppendorf). The amplification program consisted of an initial denaturation step (2 min at 95°C) followed by 40 cycles of amplification with two stages (15 s at 95°C, 1 min at 60°C). The specificity of the amplification was checked with a melting curve analysis consisting of an initial denaturation step of 15 s at 95°C followed by an increasing temperature with a step of 0.1°C, from 60 to 95°C, for 20 min. Target gene expression was normalized using the reference gene At4g26410 (Czechowski et al., 2005; Lemonnier et al., 2014) whose expression remained stable in all condition evaluated (in the different tissues, time points, and following inoculation with the PGPR strain) according to the results obtained with a second reference gene: AtUPL7 (At3g53090). Results were expressed as relative gene expression values using the 2<sup>-ΔCt</sup> method (Schmittgen and Livak, 2008).

Data are mean ± SEM of 4 biological replicates, each from an independent experiment.

## Statistical Analysis

Statistical analyses of differences for morphological traits and relative gene expressions were carried out using a nonparametric Mann-Whitney-Wilcoxon test (n < 30) unless otherwise indicated. Tests were performed using the software GraphPad Prism® version 7.0.

# RESULTS

## Isolation and Identification of the Endophytic and IAA Producing Bacterium *Bacillus megaterium* RmBm31

*Bacillus megaterium* strain RmBm31 was isolated from surface sterilized root-nodules of *Retama monosperma* growing in Algeria following an approach already used in a previous study to identify nodule-associated microbial symbionts (Muresu et al., 2008). Among the endophytic bacterial strains identified, isolate RmBm31 was selected as the best producer of IAA according to a Thin Layer Chromatography assay (TLC) (Figure S1).

By sequencing its genome (Figure 1) and performing a phylogenomic analysis, we determined that RmBm31 is a *Bacillus megaterium* strain (Figure 2 and Table S1). The genome of RmBm31 has a total of 5,817,483 bp in size with an average G+C content of 37.64% and contains 5,817 predicted coding sequences (CDS) with an average length of 817 bp. Coding regions covered 81.7% of the genome. Biological roles (annotation) were assigned to 4,318 (74.2%) genes of the predicted CDS based on similarity searches with the NCBI nonredundant protein database (nr). The remaining CDS (1,499; 25.8%) code for proteins with unknown function (17) and hypothetical proteins (1,482). Lastly, a total of 20 rRNAs (comprising eighteen 5S rRNAs, one 16S rRNA, and one 23S rRNA) together with 119 tRNAs genes were identified in RmBm31 genome (Table S2). The low number of rRNA genes can be explained by the difficulty to assemble these genes from short read technology, since the assembler collapses these genes in one contig.

The classification of CDS into functional categories according to the COG (Clusters of Orthologous Groups) (Tatusov et al., 2000; Galperin et al., 2015) database is summarized in Table S3. Poorly characterized genes, including those with only a general function prediction or an unknown function represented the majority of protein-coding genes in the annotated genome of RmBm31 (35.23%), whereas genes not assigned in COG represented 11.24% of this total number of proteins. Out of the 3,114 proteins assigned to a COG functional category, “Amino acid transport and metabolism” [E], “Energy production and conversion” [C], and “Translation, ribosomal structure and biogenesis” [J] were the three most represented across the 25 categories with 4.64, 4.44, and 4.21% of protein encoding genes, respectively. These categories were closely followed by the categories “Replication, recombination, and repair” [L] (4.09%), “Transcription” [K] (3.97%), and “Carbohydrate transport and metabolism” [G] (3.95%) (Table S3).

## The Genome of RmBm31 Contains a Large Set of Genes Contributing to Plant Beneficial Functions

Endophytic nonrhizobial species often display beneficial effects on the host plant growth and development and/or stress resilience.
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Graphical representation of the genome structure of *Bacillus megaterium* RmBm31. Circles from the outside to the inside show the positions of protein-coding genes (blue), tRNA genes (red) and rRNA genes (pink) and others (grey) on the positive (circle 1), and negative (circle 2) strands. Circles 3 and 4 show plots of AT content and AT skew plotted as the deviation from the average content of the entire sequence. Circles 5 and 6 show plots of GC content and GC skew plotted as the deviation from the average content of the entire sequence.</figcaption></figure>

Genes commonly associated with the beneficial effects of rhizospheric and endophytic PGPR on plant productivity encode for proteins involved in nutrient acquisition, including phosphate solubilization, nitrogen fixation and metabolism, and siderophore production, as well as the production of phytohormones, the emission of volatile organic and inorganic compounds, and proteins involved in root colonization traits (Santoyo et al., 2016; Liu et al., 2017). In order to determine whether RmBm31 may display one or several of these genetic traits known to be associated with Plant Growth-Promotion (PGP), we took advantage of its sequenced genome. Briefly, we established a comprehensive list of

PGP-related genes and their associated protein sequences from literature and databases searches. Homologous proteins encoded by the genome of RmBm31 were then retrieved by performing an NCBI Protein BLAST (Blastp) search (**Table 1** and **Table S4**).

Results from this analysis revealed that the genome of RmBm31 encodes proteins linked to phosphate solubilization, including an exopolyphosphatase (Bmeg_02393), and an alkaline phosphatase (Bmeg_02494, Bmeg_04574). Genes involved in transport and assimilation of inorganic phosphate (Pho regulon) were also identified (Bmeg_04258 and Bmeg_04259). In addition, RmBm31 genome contains genes for nitrate transport (Bmeg_02938), nitrate/
<figure><img src="image_5.png" /><figcaption>**FIGURE 2** | Continued</figcaption></figure>
**FIGURE 2** | Phylogenomic tree based on 38 single copy genes obtained by the maximum-likelihood method showing the phylogenetic relationship of *Bacillus megaterium* RmBm31 with the related species. Orthofinder2 and IQ-tree were used to determine the position of *Bacillus megaterium* RmBm31 strain in relation to other species within the genus *Bacillus*. The tree was rooted with *Sulfolobus acidocaldarius* DSM639. The numbers above the branches are support values obtained from 1,000 bootstrap replicates.

nitrite reduction (Bmeg_02939) and nitrite/ammonium reduction (Bmeg_02563, Bmeg_02940), but seems to lack *nif* genes that are necessary for nitrogen fixation, and therefore, most likely does not have the ability to perform this biological process unlike nitrogen-fixing rhizobia species and other rhizobacteria (Bruto et al., 2014) (Table 1 and Table S4).

Besides improving the acquisition of nutrients by plants, rhizobacteria also display phytobeneficial effects through their ability to produce phytohormones, including auxin. A large proportion of bacteria isolated from the rhizosphere are able to synthetize Indole Acetic Acid (IAA) from tryptophan via at least one of the five alternative pathways identified to date (Spaepen and Vanderleyden, 2011). In agreement with the demonstrated ability to produce IAA (Figure S1), RmBm31 genome encodes proteins known to be required for the biosynthesis of tryptophan from chorismate (Bmeg_03801, Bmeg_03802, Bmeg_03803, Bmeg_03804, Bmeg_03805) as well as the four enzymes involved in the indole-3-pyruvate (IPA) pathway of auxin synthesis, including the indole-3-pyruvate decarboxylase (encoded by the gene *ipdC*) (Bmeg_02539), and one enzyme (encoded by the gene *iaaH*) involved in the Indole-3-Acetamide (IAM) pathway (Bmeg_04800). In addition, we also identified the presence of genes involved in the synthesis of ACC deaminase, a bacterial enzyme (encoded by the gene *acdS*) able to promote plant growth by lowering the plant endogenous ethylene levels (Bmeg_03384) (Bruto et al., 2014).

Moreover, our study identified genes implicated in the production of indole, acetoin, and 2,3-butanediol, three VOC known to enhance plant growth and/or positively modify the

root system architecture (Ryu et al., 2003; Bailly et al., 2014). RmBm31 genome also possesses a gene involved in the production of nitric oxide (nos; Bmeg_04865). Lastly, the capacity of rhizobacteria to migrate toward the roots and colonize them in response to root exudates also requires a set of genes, including some involved in flagellar biosynthesis/ regulation (Cole et al., 2017) identified in the genome of RmBm31 (Table 1 and Table S4). Altogether, these results suggest that RmBm31 is an endophytic bacterium that possesses a large set of genetic features generally associated with plant growth-promoting traits of potential interest for its host *R. monosperma* and more widely for agricultural applications.

## RmBm31 Promotes the Growth and Development of Arabidopsis Seedlings

In order to determine whether the endophytic rhizobacteria strain RmBm31 displays plant growth-promoting activities, we used an *in vitro* experimental system in which *Arabidopsis* seedlings were grown axenically prior to their inoculation with RmBm31 in physical contact with their roots, or prior to their exposure to RmBm31 volatile compounds. Seven days post-inoculation, we performed shoot and root fresh weight measurements and analyzed a set of root architecture parameters (Figure 3 and Table S5). In physical contact with roots, in which condition both the actions of diffusible and volatile substances may be involved, RmBm31 significantly enhanced the root biomass of the seedlings (79% increase), the lateral root length and number (320 and 210% increase,

<table><caption>**TABLE 1** | Bacterial genes associated with plant growth-promoting traits and with rhizosphere competence identified in *Bacillus megaterium* RmBm31.</caption><thead><tr><th>Function</th><th>Gene name</th><th>Best hit (B_meg_)</th></tr></thead><tbody><tr><td>Phosphate solubilization and mineralization</td><td>pqqF, phoR, phoP, phoA, phoD, ppx</td><td>03619, 04259, 04258, 02494, 04574, 02393</td></tr><tr><td>Nitrogen assimilation and reduction</td><td>narK, nasB, nasC, nasD</td><td>02938, 02940, 02939, 02563</td></tr><tr><td>Siderophore synthesis/Fe-uptake</td><td>pvdL, pvdJ, pchA, pchC, pchD, pchE, yfmC, yfmD, yfmE, yfmF</td><td>05807, 04387, 05808, 04929, 02837, 02836, 02835, 04489</td></tr><tr><td>L-tryptophane, indole synthesis</td><td>trpA, trpB, trpC, trpD, trpF</td><td>03801, 03802, 03804, 03805, 03803</td></tr><tr><td>Auxin (Indole-3-Acetic acid) synthesis</td><td>iaaH, ipdC</td><td>04800, 02539</td></tr><tr><td>ACC (1-Aminocyclopropane-1-Carboxylate)-deamination</td><td>acdS</td><td>3384</td></tr><tr><td>Spermidine synthesis</td><td>speA, speB, speH, speE</td><td>02358, 04672, 00344, 02847, 04252, 04673, 02937</td></tr><tr><td>Acetoin, Butanediol synthesis</td><td>alsD, alsS, ilvH, ilvB, bdhA</td><td>02929, 02928, 04187, 04188, 01788</td></tr><tr><td>Nitric oxide synthesis</td><td>nos</td><td>4865</td></tr><tr><td>Hydrogen cyanide synthesis</td><td>hcnC</td><td>4991</td></tr><tr><td>2,4-Diacetylphloroglucinol synthesis</td><td>phlA, phlD</td><td>05629, 02253</td></tr><tr><td>Flagellar assembly</td><td>flgB, flgC, flgD, flgE, flgK, flgL, flgM, flgN, flhA, flhC, flhD, flhE, flhF, flhG, flhH, flhI, flhJ, flhK, flhM, flhN, flhP, flhQ, flhR, flhS, motA, motB</td><td>03672, 03671, 03662, 03661, 04591, 04590, 04593, 04592, 03650, 02625, 04581, 03670, 03669, 03668, 03667, 03666, 03665, 03663, 03658, 03657, 03654, 03653, 03652, 04599, 04580, 00123, 00124</td></tr><tr><td>Bacterial chemotaxis</td><td>cheA, cheB, cheD, cheR, cheY, cheW</td><td>03646, 03647, 03644, 03809, 03656, 03645</td></tr></tbody></table>
respectively) and the root hair length (36% increase) in comparison to the control condition (mock treatment) (Figure 3 and Table S5). These data suggest RmBm31 may be defined as a PGPR strain.

To investigate the putative contribution of volatile compounds produced by RmBm31 on the plant growth-promoting effects described above, we designed and used an *in vitro* experimental system in which the strain RmBm31 was physically separated from the *Arabidopsis* seedlings, so that only the volatile compounds potentially produced by RmBm31 could be responsible for the phenotypic effects observed. In these conditions, beneficial effects of RmBm31 were observed on the seedling shoot and root fresh weights (196 and 118% increase, respectively). In addition, the volatile compounds produced by this strain strongly modified the root system architecture (Figure 4 and Table S5). More particularly, RmBm31 volatile compounds had a large influence on the lateral root length and number (217 and 256% increase, respectively) and significantly increased (84%) the root hair length of RmBm31-

treated seedlings in comparison to the mock treatment (Figure 4 and Table S5). Overall, the most prominent changes triggered by RmBm31 on root morphology were the stimulation of the lateral root formation and development and an increase in root hair length in both experimental conditions tested (with RmBm31 in physical contact with the seedlings roots, and via the emission of volatile compounds only). These results also clearly demonstrate a strong implication of the volatile compounds produced by RmBm31 on its beneficial effects on *Arabidopsis* seedlings grown in controlled *in vitro* conditions.

## B. megaterium RmBm31 Induces Transcriptional Changes in Marker Genes Responsive to Defense-Associated Phytohormones

To study the contribution of auxin and other phytohormones on the positive changes in seedling biomass and root-system architecture

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

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

**FIGURE 3** | Phenotypic effects of *Bacillus megaterium* RmBm31 on *Arabidopsis thaliana* Col-0 seedlings 7 days post inoculation (dpi). Five-day-old seedlings were either mock-treated ("Mock") or inoculated with 10 µl of RmBm31 at 2.10⁶ CFU.ml⁻¹ in MgSO₄ 10 mM ("RmBm31"). (A) Pictures of whole seedlings (scale bar, 1 cm) and macroscopy pictures of root tips (scale bar, 1 mm) at 7 dpi. (B) Root and shoot biomass of the seedlings and quantitative phenotypic analysis of their root system architecture at 7 dpi. Data are means ± SEM of nine biological replicates (n) from three independent experiments. Stars indicate statistically significant differences between the mock- and RmBm31-treated conditions according to a Mann-Whitney-Wilcoxon test (ns, nonsignificant; \*, P < 0.05; \*,\*, P < 0.01; \*,\*,\*, P < 0.001; \*,\*,\*,\*, P < 0.0001).
<figure><img src="image_8.png" /><figcaption>B</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**FIGURE 4** | Phenotypic effects of *Bacillus megaterium* RmBm31 volatile compounds on *Arabidopsis thaliana* Col-0 seedlings 7 days post inoculation (dpi). Five-day-old seedlings were either mock-treated ("Mock") or exposed to RmBm31 volatile compounds by spotting 100 µl of RmBm31 at 2.10⁶ CFU.ml⁻¹ in MgSO₄ 10 mM on LB medium, per plate of 10 seedlings, and without physical contact with them ("RmBm31"). (A) Pictures of whole seedlings (scale bar, 1cm) and macroscopy pictures of root tips (scale bar, 1 mm) at 7 dpi. (B) Root and shoot biomass of the seedlings and quantitative phenotypic analysis of their root system architecture at 7 dpi. Data are means ± SEM of nine biological replicates (n) from three independent experiments. Stars indicate statistically significant differences between the mock- and RmBm31-treated conditions according to a Mann-Whitney-Wilcoxon test (ns, nonsignificant ; *, P < 0,05; **, P < 0,01; ***, P < 0,001; ****, P < 0,0001).</figcaption></figure>

induced by RmBm31, we examined by RT-qPCR the expression of phytohormone-responsive genes in the shoots and roots of *Arabidopsis* seedlings either inoculated with RmBm31 or exposed solely to its volatile compounds (Figure 5 and Table S7). This gene expression analysis was performed using a comprehensive set of 33 marker genes for the biosynthesis, transport or signaling pathways of the nine plant hormones (Table S6). Preference were given to marker genes that showed specific responsiveness to single exogenous hormone applications based on previous studies (Nemhauser et al., 2006; Paponov et al., 2008; De Smet et al., 2015; Pandey et al., 2016; Papadopoulou et al., 2018) and according to data in gene expression databases (BAR eFP Browser and Genevestigator).

Unexpectedly, RmBm31 did not induce any major transcriptional changes in the biosynthesis or signaling pathways of auxin, a growth-related hormone known to be implicated in PGPR-induced beneficial effects on plant growth and the root

system architecture (Zamioudis et al., 2013). Instead, several marker genes responsive to the defense-related hormone were found relatively strongly (fold change > 2) and statistically significantly induced or repressed, either in root or in shoot, in at least one of the two experimental conditions tested. More specifically, when the seedlings were inoculated with the PGPR in physical contact with their roots, we observed the induction of the expression of *ERF1B* (Ethylene), *WRKY38* (Jasmonic acid), *PR1* and *WRKY70* (Salicylic acid) in the root and/or the shoot. In the second experimental condition, in which the seedlings were solely exposed to the volatile compounds emitted by RmBm31, only *PDF1.2* (Jasmonic acid/ Ethylene), *WRKY38*, and *PR1* were significantly induced in the shoots. Hence, only two genes, *WRKY38* and *PR1*, exhibited a differential expression in RmBm31- and Mock-treated seedlings shoots, in both experimental conditions tested (Figure 5 and Table S7). Altogether, these results suggest that the plant-
<table><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="2">ROOT</th><th colspan="2">SHOOT</th></tr><tr><th>PC</th><th>VC</th><th>PC</th><th>VC</th></tr></thead><tbody><tr><td rowspan="6">Auxin (IAA)</td><td>IAA19</td><td>1,2</td><td>1,2</td><td>1,2</td><td>1,8</td></tr><tr><td>IAA29</td><td>1,5</td><td>1,4</td><td>1,2</td><td>1,8</td></tr><tr><td>LBD29</td><td>1,7</td><td>2,1</td><td>2,2</td><td>3,8</td></tr><tr><td>AXR4</td><td>1,2</td><td>1,2</td><td>1,2</td><td>1,1</td></tr><tr><td>AUX1</td><td>1,6</td><td>1,5</td><td>1,5</td><td>1,5</td></tr><tr><td>PAT1</td><td>1,1</td><td>1,3</td><td>1,5</td><td>1,1</td></tr><tr><td>NIT1</td><td>2,1</td><td>1,8</td><td>1,3</td><td>1,1</td></tr><tr><td rowspan="3">Cytokinin (CK)</td><td>ARR9</td><td>1,6</td><td>1,3</td><td>1,4</td><td>1,7</td></tr><tr><td>CRE1</td><td>2,5</td><td>2,0</td><td>1,3</td><td>1,4</td></tr><tr><td>IPT5</td><td>1,3</td><td>3,1</td><td>1,3</td><td>5,4*</td></tr><tr><td rowspan="3">Gibberellin (GA)</td><td>RGA</td><td>1,3</td><td>1,3</td><td>1,0</td><td>1,3</td></tr><tr><td>GID1A</td><td>1,2</td><td>1,5</td><td>1,1</td><td>1,8</td></tr><tr><td>GA20OX1</td><td>1,1</td><td>1,0</td><td>1,1</td><td>1,1</td></tr><tr><td rowspan="6">Brassinosteroid (BR)</td><td>MSBP1</td><td>1,1</td><td>1,1</td><td>1,2</td><td>1,1</td></tr><tr><td>BAK1</td><td>1,1</td><td>1,4</td><td>1,3</td><td>1,2</td></tr><tr><td>BRI1</td><td>1,1</td><td>1,6</td><td>1,1</td><td>1,2</td></tr><tr><td>SAUR-Ac</td><td>1,3</td><td>1,7</td><td>1,4</td><td>1,1</td></tr><tr><td>BZR1</td><td>1,1</td><td>1,5</td><td>1,3</td><td>1,1</td></tr><tr><td>DWF4</td><td>1,1</td><td>1,3</td><td>1,0</td><td>1,4</td></tr><tr><td rowspan="2">Strigolactone (SL)</td><td>D14</td><td>1,1</td><td>1,2</td><td>1,0</td><td>1,8</td></tr><tr><td>MAX1</td><td>1,3</td><td>1,2</td><td>1,6</td><td>1,1</td></tr><tr><td rowspan="2">Abscissic acid (ABA)</td><td>PYL8</td><td>1,8</td><td>1,1</td><td>1,7</td><td>1,4</td></tr><tr><td>ABA1</td><td>1,8</td><td>2,6</td><td>1,1</td><td>1,5</td></tr><tr><td rowspan="3">Ethylene (ET)</td><td>ERS2</td><td>1,1</td><td>2,0</td><td>1,2</td><td>2,1</td></tr><tr><td>ERF1B</td><td>3,6*</td><td>2,3</td><td>2,8*</td><td>1,6</td></tr><tr><td>ACS5</td><td>1,4</td><td>1,8</td><td>1,1</td><td>1,3</td></tr><tr><td rowspan="3">Jasmonic acid (JA)</td><td>PDF1.2</td><td>4,3</td><td>38,1</td><td>1,5</td><td>37,0*</td></tr><tr><td>WRKY38</td><td>3,9*</td><td>4,2</td><td>5,9*</td><td>7,1*</td></tr><tr><td>OPR3</td><td>1,5</td><td>1,1</td><td>1,3</td><td>1,2</td></tr><tr><td rowspan="3">Salicylic acid (SA)</td><td>PR1</td><td>1,7</td><td>3,9</td><td>3,0*</td><td>3,5*</td></tr><tr><td>WRKY70</td><td>3,1*</td><td>2,0</td><td>2,9*</td><td>1,9</td></tr><tr><td>MES1</td><td>1,3</td><td>1,1</td><td>1,7*</td><td>1,9</td></tr><tr><td>Camalexin (CAM)</td><td>PAD3</td><td>1,1</td><td>1,8</td><td>2,6*</td><td>2,3</td></tr></tbody></table>

**FIGURE 5** | Effects of *Bacillus megaterium* RmBm31 on the relative expression of phytohormone marker genes in *Arabidopsis thaliana* Col-0 seedlings roots and shoots. Heatmap of gene expression changes in roots and in shoots induced by RmBm31 in physical contact with the seedling roots 7 days post inoculation ("PC"), or by RmBm31 volatile compounds at 7 dpi ("VC"). Five-day-old seedlings were either mock-treated or treated with RmBm31 inoculum. Seven days post inoculation, root and shoot tissues were harvested (at midday, 8 h light) and transcript levels were quantified by qRT-PCR. The expression level of each gene was normalized to the reference gene At4g26410. Data are mean ± SEM of four biological replicates, each from an independent experiment. Fold changes (Fc) of gene expression are in red for induction and in blue for repression in comparison to the mock. Stars indicate statistically significant differences according to a Mann-Whitney-Wilcoxon test (ns, nonsignificant; *, P < 0,05).
growth promotion and root architecture changes triggered by RmBm31 might involve molecular mechanisms independent of auxin and the phytohormonal pathways well known to date to be involved in such biological processes.

DISCUSSION

**B. megaterium RmBm31 is an Endophytic Bacterium Isolated From Root Nodules of R. monosperma That Possesses Numerous Genetic Features Linked to Plant Growth-Promotion**

*Retama* species are leguminous plants of agronomic interest for biofertilization, dune stabilization, and revegetation of semiarid and arid ecosystems, notably due to their capacity to form symbiosis with N-fixing rhizobial species (Selami et al, 2014). Despite this importance, the symbiotic interaction between *Retama* and rhizobia species and with nonrhizobial endophytes remains largely understudied (Ahnia et al., 2018; Lamin et al., 2019). Here, we isolated an endophytic bacterium from *Retama monosperma* root nodules and identified it as a *Bacillus megaterium* strain by whole genome sequencing. To the best of our knowledge, this is the first report of an endophytic nonrhizobial rhizobacteria being isolated from root nodules of *Retama monosperma*. Nevertheless, other *Bacillus megaterium* strains have already been identified in nodules of several leguminous species like the common bean (*Phaseolus vulgaris* L.) (Korir et al., 2017), pigeon pea (*Cajanus cajan*) (Rajendran et al., 2008), *Medicago polymorpha* (Chinnaswamy et al., 2018), and *Medicago sativa* L. (Khalifa and Almalki, 2015). More generally, a large and diverse set of nonrhizobial species live in legume root nodules in association with N₂-fixing rhizobial species (Peix et al., 2014; Martínez-Hidalgo and Hirsch, 2017).

Similarly to rhizospheric PGPR, endophytic bacteria are able to promote plant growth and development through direct or indirect mechanisms that include an improved plant nutrient acquisition (typically nitrogen, phosphorus, and iron), the synthesis of phytohormones, or the emission of VOC (Vacheron et al., 2013). We determined whether RmBm31 possesses genes associated with these traits. Results from our analysis showed that RmBm31 possesses many of the genetic features known to be associated with the PGPR-induced plant growth-promoting traits, including the ability to improve the plant nutrient acquisition and to synthetize phytohormones and phytobeneficial volatile compounds. However, although screening for the bacterial PGP genes and PGP traits have been successfully employed as selection criteria in many studies, other studies have shown that the abundance of these traits does not always correlate with the PGPR plant growth-promoting efficiency (Compant et al., 2019). These findings highlight the usefulness of performing plant-PGPR interaction assays when inferring the agronomical potential of newly isolated PGPR.

**B. megaterium RmBm31 Promotes Plant Growth and Root Development in *Arabidopsis* Through the Production of Volatile Compounds**

In order to evaluate the plant growth-promoting capacities of RmBm31, we set up an *in vitro* experimental system in which the model plant species *Arabidopsis thaliana* was co-cultivated with this PGPR. Moreover, two different conditions were used to investigate solely the effects of the bacterial volatile compounds in comparison to a second condition in which bacterial diffusible substances and the physical contact by itself could also be involved. Our results revealed that RmBm31 display marked plant growth-promoting activities (increase of the shoot and root biomass, lateral root length and number, and root hair length), both in physical contact with seedling roots (putatively including the production of diffusible and/or volatile compounds) and *via* the production of volatile compounds only. Based on these results, RmBm31 can be considered as a PGPR. Remarkably, this endophytic strain was originally isolated in root nodules of *Retama monosperma* but displays strong plant growth-promoting properties on the nonhost, non nodulating species *Arabidopsis thaliana*. This finding suggests that RmBm31 may have the potential to efficiently improve the growth of a wide range of plant species, and therefore be of interest for agronomical applications. However, its PGP activities in soil condition and on crop species remain to be investigated.

Results from our phenotypic analyses also demonstrated that colonization of the plant roots by RmBm31 is not a prerequisite for its phytobeneficial activities since these latter are also observed in condition where the effects are induced solely by volatile compounds. These results are fully in agreement with previous studies which demonstrated that microbial volatile organic compounds are important inducers of plant growth promotion (Ryu et al., 2003; Blom et al., 2011). Noteworthy, RmBm31 grew better on LB medium than in physical contact with the seedlings roots on MS medium. Therefore volatile compounds produced may be different or in different concentrations in the two experimental conditions (Blom et al., 2011). This could explain the difference observed in the intensity of the PGP effects between the two cocultivation conditions. Nevertheless, similar PGPR-induced plant growth-promoting activities and root developmental changes were observed in both systems. Altogether, our results suggest that the volatile compounds emitted by RmBm31 contribute to a large extent to its plant growth-promoting activities.

Major RmBm31-induced changes on root morphology are the stimulation of lateral root formation and root hair development. These root system architecture modifications are reminiscent of those observed with other PGPR strains (Ortíz-Castro et al., 2011; Zamioudis et al., 2013; Verbon and Liberman, 2016; Palacio-Rodríguez et al., 2017). Noteworthy, auxin signaling and transport have been shown to play an important role in these PGPR-induced effects (Ortíz-Castro et al., 2011; Zamioudis et al., 2013; Zamioudis et al., 2015). The major contribution of volatile
compounds in the plant growth-promoting effects observed allow to rule out the possibility bacterial auxin or other nonvolatile phytohormones produced by RmBm31 could be the main sources of its phytobeneficial properties. However, a possible regulation by the PGPR volatile compounds of the plant auxin biosynthesis and/or signaling pathways leading to the PGP effects observed could still be at play. Indole might represent a good candidate since treatment of *Arabidopsis* seedlings with this volatile molecule lead to similar beneficial effects on the shoot and root biomass and root system architecture than those triggered by RmBm31 (Bailly et al., 2014). Remarkably, indole is known to participate in auxin-controlled lateral root initiation and development by interfering with auxin signaling, beyond its role as auxin synthesis precursor (Bailly et al., 2014).

## B. megaterium RmBm31 Plant Growth-Promoting Effects are Accompanied With Transcriptional Changes in Defense-Related Pathways

Results from our studies demonstrate that the production of IAA by RmBm31 is not the major contributor of its PGP activities since these latter were observed in the *in vitro* experimental system in which only volatile compounds effects were involved. However, these findings do not exclude a possible regulation of the plant IAA biosynthesis and/or signaling pathways leading to the PGP effects observed by the volatile compounds emitted by RmBm31. Evidence in the literature indeed reveals a key role played by these pathways in PGPR-triggered root system architecture changes similar to those observed in response to RmBm31 inoculation (Zamioudis et al., 2013). However, the beneficial effects of other strains of *Bacillus megaterium* have been shown to involve auxin-independent molecular mechanisms (Kishore et al., 2005; López-Bucio et al., 2007). These results are also in agreement with a study performed on a strain of *Bacillus aryabhattai*, a subspecies closely related to *Bacillus megaterium* (Park et al., 2017).

Nevertheless, the growth-promotion effects of RmBm31 are likely to involve changes in the host plant hormonal pathways. To test this hypothesis, one could use mutants defecives in particular steps of the pathways and investigate whether mutations affect the growth promotion of the plant by the PGPR. However, interpretation of the results may be ambiguous, notably given the strong and/or pleiotropic phenotypic effects sometimes conferred by the mutation(s) (eg. auxin and cytokinin perception triple mutants *tir1afb2afb3* and *cre1-12/ahk2-2tk/ahk3-3*) (Ortíz-Castro et al., 2008; Zamioudis et al., 2013). In recent years, genome-wide gene expression profiling and targeted gene expression analyses with pathway-specific marker genes have emerged as useful complementary or alternative strategies (Papadopoulou et al., 2018) to screen for molecular mechanisms that may be involved in a particular biological process. We chose to apply such strategy by investigating the capacity of RmBm31 to induce transcriptional changes in auxin response/signaling, as well as other phytohormonal pathways in the seedling roots and shoots. Interestingly, the gene expression analysis we performed failed to identify major RmBm31-triggered changes in the expression of genes involved in the pathways of the plant

growth-related hormones auxin, cytokinin and gibberellins. Instead, the only two marker genes that were found to be transcriptionally regulated by RmBm31, in the two experimental conditions tested (and therefore could be associated with the PGP properties of this PGPR strain) were *WRKY38* and *PR1*, two genes commonly associated with the plant defense pathways (Van Loon et al., 2006; Kim et al., 2008; Xie et al., 2010).

## CONCLUSION AND PERSPECTIVES

We report here on the identification and characterization of a *Bacillus megaterium* endophytic strain, RmBm31, isolated from root nodules of the legume species *Retama monosperma*, a perennial leguminous species growing in poor and high salinity soils. Our study reveals that RmBm31 is an IAA-producing endophytic bacterium that possesses a large set of genes associated with plant growth-promoting traits. Using the model plant species *Arabidopsis thaliana*, we demonstrate this strain display beneficial effects on plant growth and root development *via* the production of volatile compounds. Remarkably, our study suggests these effects might involve auxin-independent signaling mechanisms. Future work should identify the volatile compounds and elucidate the molecular circuitry involved in the efficiency of this plant-PGPR interaction. Whether RmBm31 is endophytic in *Arabidopsis* also merits to be investigated. In addition, the potential beneficial effects on RmBm31 in soil assays and on the plant abiotic stress tolerance remains to be investigated. In longer term this new knowledge may help develop PGPR-based strategies for crop protection and productivity improvement in sustainable agriculture.

## DATA AVAILABILITY STATEMENT

The datasets generated for this study can be found in the The Whole Genome Shotgun project of *Bacillus megaterium* RmBm31 has been deposited at DDBJ/ENA/GenBank under the accession VYTX00000000 (Bioproject PRJNA566080, Biosample SAMN12779970). The version described in this paper is version VYTX01000000.

## AUTHOR CONTRIBUTIONS

CV, MD, AD, and BM designed the experiments. MD, AD, BM, CV, LM, and MK-H performed the experiments and data analyses. CV, AD, and MD wrote the manuscript. CV, AD, MD, BM, JV, PC-T, MH, and TB discussed the data and revised the manuscript. All authors approved the final manuscript.

## FUNDING

This work was funded by the French Ministry of Higher Education, Research and Innovation (AD, PhD grant) and the by Erasmus mundus-Alidrissi A scholarship scheme for exchange and cooperation between Europe and North Africa (MAD). This study was also supported by the 2015-2020 State-Region Planning
Contracts (CPER), the European Regional Development Fund (FEDER), the Centre National de la Recherche Scientifique (CNRS), the University of Poitiers and USTO-MB University.

## ACKNOWLEDGMENTS

The authors gratefully acknowledge the ImageUP platform and Bruno Merceron for technical assistance in macroscopy analyses. We also want to express our appreciation to students who

contributed to this work, especially Laurie Taylor, Baptiste Biet and David Landry.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00124/full#supplementary-material

## REFERENCES

Ahnia, H., Bourebaba, Y., Duran, D., Boulila, F., Palacios, J. M., Rey, L., et al. (2018). *Bradyrhizobium algeriense* sp. nov., a novel species isolated from effective nodules of *Retama sphaerocarpa* from Northeastern Algeria. *Syst. Appl. Microbiol.* 41, 333–339. doi: 10.1016/j.syapm.2018.03.004

Bailly, A., Groenhagen, U., Schulz, S., Geisler, M., Eberl, L., and Weisskopf, L. (2014). The inter-kingdom volatile signal indole promotes root development by interfering with auxin signalling. *Plant J.* 80, 758–771. doi: 10.1111/tpj.12666

Barea, J. M., Pozo, M. J., Azcón, R., and Azcón-Aguilar, C. (2005). Microbial co-operation in the rhizosphere. *J. Exp. Bot.* 56, 1761–1778. doi: 10.1093/jxb/eri197

Bhattacharyya, C., Bakshi, U., Mallick, I., Mukherji, S., Bera, B., and Ghosh, A. (2017). Genome-guided insights into the plant growth promotion capabilities of the physiologically versatile bacillus aryabhattai strain AB211. *Front. Microbiol.* 8, 411. doi: 10.3389/fmicb.2017.00411

Blom, D., Fabbri, C., Connor, E. C., Schiestl, F. P., Klauser, D. R., Boller, T., et al. (2011). Production of plant growth modulating volatiles is widespread among rhizosphere bacteria and strongly depends on culture conditions. *Environ. Microbiol.* 13, 3047–3058. doi: 10.1111/j.1462-2920.2011.02582.x

Bordiec, S., Paquis, S., Lacroix, H., Dhondt, S., Ait Barka, E., Kauffmann, S., et al. (2011). Comparative analysis of defence responses induced by the endophytic plant growth-promoting rhizobacterium *Burkholderia phytofirmans* strain PsJN and the non-host bacterium Pseudomonas syringae pv. pisi in grapevine cell suspensions. *J. Exp. Bot.* 62, 595–603. doi: 10.1093/jxb/erq291

Box, M. S., Coustham, V., Dean, C., and Mylne, J. S. (2011). Protocol: A simple phenol-based method for 96-well extraction of high quality RNA from *Arabidopsis*. *Plant Methods* 7, 7. doi: 10.1186/1746-4811-7-7

Bruto, M., Prigent-Combaret, C., Muller, D., and Moenne-Locoz, Y. (2014). Analysis of genes contributing to plant-beneficial functions in plant growth-promoting rhizobacteria and related proteobacteria. *Sci. Rep.* 4, 6261. doi: 10.1038/srep06261

Chinnaswamy, A., Coba de la Peña, T., Stoll, A., de la Peña Rojo, D., Bravo, J., Rincón, A., et al. (2018). A nodule endophytic Bacillus megaterium strain isolated from Medicago polymorpha enhances growth, promotes nodulation by *Ensifer medicae* and alleviates salt stress in alfalfa plants. *Ann. Appl. Biol.* 172, 295–308. doi: 10.1111/aab.12420

Cole, B. J., Feltcher, M. E., Waters, R. J., Wetmore, K. M., Mucyn, T. S., Ryan, E. M., et al. (2017). Genome-wide identification of bacterial plant colonization genes. *PloS Biol.* 15, e2002860. doi: 10.1371/journal.pbio.2002860

Compant, S., Samad, A., Faist, H., and Sessitsch, A. (2019). A review on the plant microbiome: Ecology, functions, and emerging trends in microbial application. *J. Adv. Res.* 19, 29–37. doi: 10.1016/j.jare.2019.03.004

Czechowski, T., Stitt, M., Altmann, T., Udvardi, M. K., and Scheible, W. R. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization in *Arabidopsis*. *Plant Physiol.* 139, 5–17.

De Smet, S., Cuypers, A., Vangronsveld, J., and Remans, T. (2015). Gene Networks Involved in Hormonal Control of Root Development in *Arabidopsis thaliana*: A Framework for Studying Its Disturbance by Metal Stress. *Int. J. Mol. Sci.* 16, 19195–19224. doi: 10.3390/ijms160819195

Drogue, B., Sanguin, H., Chamam, A., Mozar, M., Llauro, C., Panaud, O., et al. (2014). Plant root transcriptome profiling reveals a strain-dependent response during *Azospirillum*-rice cooperation. *Front. Plant Sci.* 5, 607. doi: 10.3389/fpls.2014.00607

Emms, D. M., and Kelly, S. (2015). OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. *Genome Biol.* 16, 157. doi: 10.1186/s13059-015-0721-2

Galperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2015). Expanded microbial genome coverage and improved protein family annotation in the COG database. *Nucleic Acids Res.* 43, D261–D269. doi: 10.1093/nar/gku1223

Gutiérrez-Luna, F. M., López-Bucio, J., Altamirano-Hernández, J., Valencia-Cantero, E., Reyes de la Cruz, H., and Macías-Rodríguez, L. (2010). Plant growth-promoting rhizobacteria modulate root-system architecture in *Arabidopsis thaliana* through volatile organic compound emission. *Symbiosis* 51, 75–83. doi: 10.1007/s13199-010-0066-2

Hardoim, P. R., Van Overbeek, L. S., Berg, G., Pirtila, A. M., Compant, S., Campisano, A., et al. (2015). The hidden world within plants: ecological and evolutionary considerations for defining functioning of microbial endophytes. *Microbiol. Mol. Biol. Rev.* 79, 293–320. doi: 10.1128/MMBR.00050-14

Hernández-Calderón, E., Aviles-Garcia, M. E., Castulo-Rubio, D. Y., Macías-Rodríguez, L., Ramírez, V. M., Santoyo, G., et al. (2018). Volatile compounds from beneficial or pathogenic bacteria differentially regulate root exudation, transcription of iron transporters, and defense signaling pathways in Sorghum bicolor. *Plant Mol. Biol.* 96, 291–304. doi: 10.1007/s11103-017-0694-5

Hoang, D. T., Chernomor, O., Von Haeseler, A., Minh, B. Q., and Vinh, L. S. (2018). UFBoot2: Improving the ultrafast bootstrap approximation. *Mol. Biol. Evol.* 35, 518–522. doi: 10.1093/molbev/msx281

Kandaswamy, R., Ramasamy, M. K., Palanivel, R., and Balasundaram, U. (2019). Impact of *pseudomonas putida* RRF3 on the root transcriptome of rice plants: Insights into defense response, secondary metabolism and root exudation. *J. Biosci.* 44, 98. doi: 10.1007/s12038-019-9922-2

Khalifa, A. Y. Z., and Almalki, M. A. (2015). Isolation and characterization of an endophytic bacterium, *Bacillus megaterium* BMN1, associated with root-nodules of *Medicago sativa* L. growing in Al-Ahsaa region, Saudi Arabia. *Ann. Microbiol.* 65, 1017–1026. doi: 10.1007/s13213-014-0946-4

Kim, K. C., Lai, Z., Fan, B., and Chen, Z. (2008). Arabidopsis WRKY38 and WRKY62 transcription factors interact with histone deacetylase 19 in basal defense. *Plant Cell* 20, 2357–2371. doi: 10.1105/tpc.107.055566

Kishore, G. K., Pande, S., and Podile, A. R. (2005). Phylloplane bacteria increase seedling emergence, growth and yield of field-grown groundnut (*Arachis hypogaea* L.). *Lett. Appl. Microbiol.* 40, 260–268. doi: 10.1111/j.1472-765X.2005.01664.x

Korir, H., Mungai, N. W., Thuita, M., Hamba, Y., and Masso, C. (2017). Co-inoculation effect of rhizobia and plant growth promoting rhizobacteria on common bean growth in a low phosphorus soil. *Front. Plant Sci.* 8, 141. doi: 10.3389/fpls.2017.00141

López-Bucio, J., Campos-Cuevas, J. C., Hernández-Calderón, E., Velásquez-Becerra, C., Farías-Rodríguez, R., Macías-Rodríguez, L. I., et al. (2007). *Bacillus megaterium* rhizobacteria promote growth and alter root-system architecture through an auxin- and ethylene-independent signaling mechanism in *Arabidopsis thaliana*. *Mol. Plant Microbe Interact.* 20, 207–217. doi: 10.1094/MPMI-20-2-0207

Lamin, H., Alami, S., Bouhnik, O., Elfaik, S., Abdelmoumen, H., Bedmar, E. J., et al. (2019). Nodulation of *retama monosperma* by *Ensifer aridii* in an abandoned lead mine soils in eastern morocco. *Front. Microbiol.* 10, 1456. doi: 10.3389/fmicb.2019.01456

Lemonnier, P., Gaillard, C., Veillet, F., Verbeke, J., Lemoine, R., Coutos-Thevenot, P., et al (2014). Expression of Arabidopsis sugar transport protein STP13 differentially affects glucose transport activity and basal resistance to Botrytis cinerea. *Plant Mol. Biol.* 85, 473–484.

Liu, H., Carvalhais, L. C., Crawford, M., Singh, E., Dennis, P. G., Pieterse, C. M. J., et al. (2017). Inner Plant Values: Diversity, Colonization and Benefits from Endophytic Bacteria. *Front. Microbiol.* 8, 2552. doi: 10.3389/fmicb.2017.02552
Lobet, G., Pagès, L., and Draye, X. (2011). A novel image-analysis toolbox enabling quantitative analysis of root system architecture. *Plant Physiol.* 157, 29–39. doi: 10.1104/pp.111.179895

Lugtenberg, B., and Kamilova, F. (2009). Plant-growth-promoting rhizobacteria. *Annu. Rev. Microbiol.* 63, 541–556. doi: 10.1146/annurev.micro.62.081307.162918

Marçais, G., and Kingsford, C. (2011). A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. *Bioinformatics* 27, 764–770. doi: 10.1093/bioinformatics/btr011

Martínez-Hidalgo, P., and Hirsch, A. M. (2017). The nodule microbiome: N2-fixing rhizobia do not live alone. *Phytobiomes* 1, 70–82. doi: 10.1094/PBIOMES-12-16-0019-RVW

Miotto-Vilanova, L., Jacquard, C., Courteaux, B., Wortham, L., Michel, J., Clement, C., et al. (2016). *Burkholderia phytofirmans* PsJN confers grapevine resistance against botrytis cinerea via a direct antimicrobial effect combined with a better resource mobilization. *Front. Plant Sci.* 7, 1236. doi: 10.3389/fpls.2016.01236

Muresu, R., Polone, E., Sulas, L., Baldan, B., Tondello, A., Delogu, G., et al. (2008). Coexistence of predominantly nonculturable rhizobia with diverse, endophytic bacterial taxa within nodules of wild legumes. *FEMS Microbiol. Ecol.* 63, 383–400. doi: 10.1111/j.1574-6941.2007.00424.x

Nemhauser, J. L., Hong, F., and Chory, J. (2006). Different plant hormones regulate similar processes through largely nonoverlapping transcriptional responses. *Cell* 126, 467–475. doi: 10.1016/j.cell.2006.05.050

Nguyen, L. T., Schmidt, H. A., Von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. *Mol. Biol. Evol.* 32, 268–274. doi: 10.1093/molbev/msu300

Ortíz-Castro, R., Valencia-Cantero, E., and López-Bucio, J. (2008). Plant growth promotion by *Bacillus megaterium* involves cytokinin signaling. *Plant Signal Behav.* 3, 263–265. doi: 10.4161/psb.3.4.5204

Ortíz-Castro, R., Diaz-Perez, C., Martinez-Trujillo, M., Del Río, R. E., Campos-García, J., and López-Bucio, J. (2011). Transkingdom signaling based on bacterial cyclodipeptides with auxin activity in plants. *Proc. Natl. Acad. Sci. U. S. A.* 108, 7253–7258. doi: 10.1073/pnas.1006740108

Palacio-Rodríguez, R., Coria-Arellano, J. L., Lizbeth, J., López-Bucio, J., Sánchez-Salas, J., Muro-Pérez, G., et al. (2017). Halophilic rhizobacteria from *Distichlis spicata* promote growth and improve salt tolerance in heterologous plant hosts. *Symbiosis* 73, 179–189. doi: 10.1007/s13199-017-0481-8

Pandey, A., Sharma, M., and Pandey, G. K. (2016). Emerging roles of strigolactones in plant responses to stress and development. *Front. Plant Sci.* 7, 434. doi: 10.3389/fpls.2016.00434

Papadopoulou, G. V., Maedicke, A., Grosser, K., Van Dam, N. M., and Martinez-Medina, A. (2018). Defence signalling marker gene responses to hormonal elicitation differ between roots and shoots. *AoB Plants* 10, ply031. doi: 10.1093/aobpla/ply031

Paponov, I. A., Paponov, M., Teale, W., Menges, M., Chakrabortee, S., Murray, J. A., et al. (2008). Comprehensive transcriptome analysis of auxin responses in Arabidopsis. *Mol. Plant* 1, 321–337. doi: 10.1093/mp/ssm021

Park, Y. G., Mun, B. G., Kang, S. M., Hussain, A., Shahzad, R., Seo, C. W., et al. (2017). *Bacillus aryabhattai* SRB02 tolerates oxidative and nitrosative stress and promotes the growth of soybean by modulating the production of phytohormones. *PloS One* 12, e0173203. doi: 10.1371/journal.pone.0173203

Peix, A., Ramírez-Bahena, M. H., Velázquez, E., and Bedmar, E. J. (2014). Bacterial associations with legumes. *CRC Crit. Rev. Plant Sci.* 34, 17–42. doi: 10.1080/07352689.2014.897899

Pieterse, C. M., Zamioudis, C., Berendsen, R. L., Weller, D. M., Van Wees, S. C., and Bakker, P. A. (2014). Induced systemic resistance by beneficial microbes. *Annu. Rev. Phytopathol.* 52, 347–375. doi: 10.1146/annurev-phyto-082712-102340

Rajendran, G., Sing, F., Desai, A. J., and Archana, G. (2008). Enhanced growth and nodulation of pigeon pea by co-inoculation of *Bacillus* strains with Rhizobium spp. *Bioresour. Technol.* 99, 4544–4550. doi: 10.1016/j.biortech.2007.06.057

Ryu, C. M., Farag, M. A., Hu, C. H., Reddy, M. S., Wei, H. X., Pare, P. W., et al. (2003). Bacterial volatiles promote growth in Arabidopsis. *Proc. Natl. Acad. Sci. U. S. A.* 100, 4927–4932. doi: 10.1073/pnas.0730845100

Santoyo, G., Moreno-Hagelsieg, G., Orozco-Mosqueda Mdel, C., and Glick, B. R. (2016). Plant growth-promoting bacterial endophytes. *Microbiol. Res.* 183, 92–99. doi: 10.1016/j.micres.2015.11.008

Schmittgen, T. D., and Livak, K. J. (2008). Analyzing real-time PCR data by the comparative C(T) method. *Nat. Protoc.* 3, 1101–1108. doi: 10.1038/nprot.2008.73

Schneider, C. A., Rasband, W. S., and Eliceiri, K. W. (2012). NIH Image to imageJ: 25 years of image analysis. *Nat. Methods* 9, 671–675. doi: 10.1038/nmeth.2089

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. *Bioinformatics* 30, 2068–2069. doi: 10.1093/bioinformatics/btu153

Selami, N., Auriac, M.-C., Catrice, O., Capela, D., Kaid-Harche, M., and Timmers, T. (2014). Morphology and anatomy of root nodules of *Retama monosperma* (L.)Boiss. *Plant Soil* 379, 109–119. doi: 10.1007/s11104-014-2045-5

Sethia, B., Mustafa, M., Manohar, S., Patil, S. V., Jayamohan, N. S., and Kumudini, B. S. (2015). Indole acetic acid production by fluorescent Pseudomonas spp. from the rhizosphere of Plectranthus amboinicus (Lour.) Spreng. and their variation in extragenic repetitive DNA sequences. *Indian J. Exp. Biol.* 53, 342–349.

Spaepen, S., and Vanderleyden, J. (2011). Auxin and plant-microbe interactions. *Cold Spring Harb. Perspect. Biol.* 3, a001438. doi: 10.1101/cshperspect.a001438

Stothard, P., and Wishart, D. S. (2005). Circular genome visualization and exploration using CGView. *Bioinformatics* 21, 537–539. doi: 10.1093/bioinformatics/bti054

Tatusov, R. L., Galperin, M. Y., Natale, D. A., and Koonin, E. V. (2000). The COG database: a tool for genome-scale analysis of protein functions and evolution. *Nucleic Acids Res.* 28, 33–36. doi: 10.1093/nar/28.1.33

Vacheron, J., Desbrosses, G., Bouffaud, M. L., Touraine, B., Moenne-Loccoz, Y., Muller, D., et al. (2013). Plant growth-promoting rhizobacteria and root system functioning. *Front. Plant Sci.* 4, 356. doi: 10.3389/fpls.2013.00356

Van Loon, L. C., Rep, M., and Pieterse, C. M. (2006). Significance of inducible defense-related proteins in infected plants. *Annu. Rev. Phytopathol.* 44, 135–162. doi: 10.1146/annurev.phyto.44.070505.143425

Verbon, E. H., and Liberman, L. M. (2016). Beneficial Microbes Affect Endogenous Mechanisms Controlling Root Development. *Trends Plant Sci.* 21, 218–229. doi: 10.1016/j.tplants.2016.01.013

Wu, S., Zhu, Z., Fu, L., Niu, B., and Li, W. (2011). WebMGA: a customizable web server for fast metagenomic sequence analysis. *BMC Genomics* 12, 444. doi: 10.1186/1471-2164-12-444

Xie, C., Zhou, X., Deng, X., and Guo, Y. (2010). PKS5, a SNF1-related kinase, interacts with and phosphorylates NPR1, and modulates expression of WRKY38 and WRKY62. *J. Genet. Genomics* 37, 359–369. doi: 10.1016/S1673-8527(09)60054-0

Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., and Madden, T. L. (2012). Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. *BMC Bioinf.* 13, 134. doi: 10.1186/1471-2105-13-134

Zamioudis, C., Mastranesti, P., Dhonukshe, P., Blilou, I., and Pieterse, C. M. (2013). Unraveling root developmental programs initiated by beneficial *Pseudomonas* spp. bacteria. *Plant Physiol.*, 162, 304–318. doi: 10.1104/pp.112.212597

Zamioudis, C., Korteland, J., Van Pelt, J. A., Van Hamersveld, M., Dombrowski, N., Bai, Y., et al. (2015). Rhizobacterial volatiles and photosynthesis-related signals coordinate MYB72 expression in Arabidopsis roots during onset of induced systemic resistance and iron-deficiency responses. *Plant J.* 84, 309–322. doi: 10.1111/tpj.12995

Zhang, H., Sun, Y., Xie, X., Kim, M. S., Dowd, S. E., and Pare, P. W. (2009). A soil bacterium regulates plant acquisition of iron via deficiency-inducible mechanisms. *Plant J.* 58, 568–577. doi: 10.1111/j.1365-313X.2009.03803.x

Zhou, C., Guo, J., Zhu, L., Xiao, X., Xie, Y., Zhu, J., et al. (2016). *Paenibacillus polymyxa* BFKC01 enhances plant iron absorption via improved root systems and activated iron acquisition mechanisms. *Plant Physiol. Biochem.* 105, 162–173. doi: 10.1016/j.plaphy.2016.04.025

Zou, C., Li, Z., and Yu, D. (2010). *Bacillus megaterium* strain XTBG34 promotes plant growth by producing 2-pentylfuran. *J. Microbiol.* 48, 460–466. doi: 10.1007/s12275-010-0068-z

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Dahmani, Desrut, Moumen, Verdon, Mermouri, Kacem, Coutos-Thévenot, Kaid-Harche, Bergès and Vriet. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.