Characterization of Bradyrhizobium spp. Nodulating Lupinus cosentinii and L. luteus Microsymbionts in Morocco

In this work, we analyzed the diversity of the nodule-forming bacteria associated with Lupinus luteus and Lupinus cosentinii grown in the Maamora Cork oak forest acidic soils in Morocco. The phenotypic analysis showed the high diversity of the strains nodulating the two lupine's species. The strains were not tolerant to acidity or high alkalinity. They do not tolerate salinity or high temperatures either. The strains isolated from L. luteus were more tolerant to antibiotics and salinity than those isolated from L. cosentinii. The plant growth promoting (PGP) activities of our strains are modest, as among the 28 tested isolates, only six produced auxins, six produced siderophores, whereas three solubilized phosphates. Only two strains possess the three activities. The rrs gene sequences from eight representative strains selected following ARDRA and REP-PCR results revealed that they were members of the genus Bradyrhizobium. Six strains were then retained for further molecular analysis. The glnII, recA, gyrB, dnaK, and rpoB housekeeping gene sequence phylogeny showed that some strains were close to B. lupini LMG28514T whereas others may constitute new genospecies in the genus Bradyrhizobium. The strains were unable to nodulate Glycine max and Phaseolus vulgaris and effectively nodulated L. luteus, L. cosentinii, L. angustifolius, Chamaecytisus albidus, and Retama monosperma. The nodC and nodA symbiotic gene phylogenies showed that the strains are members of the genistearum symbiovar.


INTRODUCTION
Legumes (Fabaceae) are of great importance in different aspects, whether in agriculture, medicine, or ecology. Some species of this plant family can establish a specific symbiosis with rhizobia. The latter are soil bacteria normally found in the rhizosphere and in large numbers on the root surface, because of the nutrients secreted by the plant as root exudates (Sessitsch et al., 2002). In the rhizobia-legume symbiosis, bacteria provide plants with nitrogen compounds and plants provide them with carbon sources. A successful symbiotic interaction is based on a molecular dialogue between the two partners and generally involves rhizobial lipo-chitooligosaccharide signals called Nod factors (NF) (Wang et al., 2018;Lindström and Mousavi, 2020). These nodulation factors secreted by the rhizobia and perceived by the root cell receptors will initiate molecular and physiological responses in the plant, leading to the formation of the nodule, which is the site of fixation of the nitrogen (Clúa et al., 2018). Nitrogen fixation and uptake contributes to the high seed protein content and also provides residual N for the following crops, particularly cereals (Kermah et al., 2018).
Rhizobia are also collectively called root or shoot nodule bacteria, comprise more than 238 species in 18 genera regrouped in two clades; however, the description of rhizobial species included only 23% of the 19,000 legume species that exist throughout the world (Shamseldin et al., 2017).
Lupinus is one of the most important genera in Fabaceae with a rich diversity of species divided into Mediterranean "Old World" species and American "New World" species. Wild lupines are distributed over a wide range of climates around the world, from the coldest to the warmest and from the driest to the most humid (Wolko et al., 2011). On the other hand, the centers of diversity of Lupinus are the Americas, the Mediterranean region, and the regions of North and East Africa (Susek et al., 2016).
Wild lupines are toxic to humans and unpalatable to animals because their seeds contain alkaloids. Lupine cultivation improved in the 1930's, when breeders discovered sweet varieties of white lupine. Thus, its culture was stimulated in different European countries (Kohajdova et al., 2011). Nowadays, lupines are mainly consumed by ruminants. They are used in human food for the manufacture of ingredients used in the bakery and pastry industries in Europe and Australia (Hall et al., 2005;Smith et al., 2006). The current cultivated varieties of lupines, known as sweet varieties, are devoid of bitterness such as white and yellow lupines, which are used for human consumption. The white lupine is exploited generally as flour, recommended in low glycemic index protein, and gluten-free diets (Kohajdova et al., 2011;Yorgancilar and Bilgiçli, 2014).
The extensive and ongoing development of molecular biology tools over the last 20 years has facilitated the identification of new nodulating bacteria and has resulted in significant changes in the classification and proposition of new and different species. Lupines were first thought to be mainly nodulated by slowgrowing rhizobia of the genus Bradyrhizobium, although fastgrowing strains associated with lupines have been identified since 1988 (Miller and Pepper, 1988). Lupinus is currently reported as a promiscuous host that can be nodulated by different symbiotic bacteria belonging to the genera Bradyrhizobium, Ochrobactrum, Microvirga, Phyllobacterium, Neorhizobium, and Rhizobium (Jarabo-Lorenzo et al., 2003;Trujillo et al., 2005;Velázquez et al., 2010;Ardley et al., 2012;Bourebaba et al., 2016;Msaddak et al., 2018;Tounsi-Hammami et al., 2019;Missbah El Idrissi et al., 2020).
There has probably been a horizontal transfer of symbiotic genes between different Bradyrhizobium spp. associated with Lupinus spp. and other legumes of the Genisteae tribe, and it is very likely that the Microvirga and Ochrobactrum species isolated from different lupins have obtained their symbiotic genes from other more common rhizobial genera (Andrews et al., 2018).
Lateral transfer of specific symbiosis genes into rhizobial genera is an important mechanism for legumes to form symbiosis with selected rhizobia more adapted to particular soils and a legume-specific rhizobium strain symbiosis can develop in specific habitats.
In this work, we isolated 36 bacteria from the root nodules of two Lupine species, Lupinus cosentinii and Lupinus luteus, grown in soils of the Maamora forest, one of the most important Cork Oak Forest in the world (Aafi et al., 2005). The bacteria were first screened for their nodulation gene possession as a marker for their belonging to symbiotic microsymbionts and then characterized by molecular and phenotypic analysis.

Bacterial Strains
All the 36 strains were isolated from different soils in Maamora forest in the region of Rabat known by its mild climate (Aafi et al., 2005) (Supplementary Table 1). Twenty isolates were isolated from L. cosentinii and 16 were isolated from L. luteus plants root nodules, according to the method of Howieson and Dilworth (2016). The root segments bearing nodules were collected, washed under running tap water and surface sterilized by immersion in 5% sodium hypochlorite for 3 min, and finally washed seven times with sterile distilled water. Each nodule was then crushed, and the extract was streaked onto plates of YMA. Plates were incubated for 15 days at 28 • C, and then the single colonies obtained were checked for purity by repeated streaking on YMA medium, supplemented with Congo red dye. The incubation period of 2 weeks was sufficient to check for the growth of any slow-growing symbiotic bacteria on YEM medium. Pure isolates were maintained at −20 • C in 50% (m/v) glycerol and at 4 • C in YMA that was routinely used for rhizobial culture.

Nodulation and Host Range Experiments
The seeds of the two lupines were surface disinfected with 70% ethanol for 1 min and later sterilized with 5% sodium hypochlorite for 15 min. The seeds were scarified with concentrated sulfuric acid for 10 min and germinated in water agar (0.6% w/v) containing Petri plates. The seedlings were then transferred to Gibson tubes (Howieson and Dilworth, 2016) and inoculated with 1 ml per plant of a suspension containing ∼10 8 isolates CFU ml −1 . The plants were cultivated at 26 • C for 60 days under a 16.0/8.0 h light/dark photoperiod. Eight weeks after inoculation, the plants were checked for nodules appearance and plant aspect. Indirect effectiveness of the nodules for N 2 fixation was estimated by visual assay of red leghemoglobin presence in cross-sections and by the dark green intensity of the leaves compared to uninoculated control plants.

DNA Extraction and REP-PCR Fingerprinting
Bacteria were grown on solid Tryptone-Yeast Extract (TY) slants (Beck et al., 1993) for 96 h at 28 • C, and colonies were suspended in 2 ml of distilled water. The suspensions were then centrifuged twice at 3,000 r/min for 15 min. The resulting pellet was treated with proteinase K (20 mg/ml), and total DNA was extracted as previously described by Guerrouj et al. (2013). Rep-PCR (Repetitive Extragenic Palindromic Polymerase Chain Reaction) using ERIC1R and ERIC2 primers (De Bruijn, 1992) was used for amplification of the DNA. The DNA template was denatured for 5 min at 95 • C, and PCR was carried out for 35 cycles (94 • C for 30 s, 52 • C for 1 min, and 72 • C for 1 min), with an elongation step at 72 • C for 7 min. PCR products were analyzed by horizontal electrophoresis in 2% agarose gels in TAE buffer at 55 V for 4 h. Gels were stained with ethidium bromide, visualized under UV radiation, and photographed with a digital camera. Cluster analysis of the fingerprints obtained was performed with Gelcompar II program (version 2.15) (Supplementary Figure 1).
All the amplifications consisted in an initial 3-min denaturation at 94 • C and 35 cycles of 30-s denaturation at 94 • C, 30 s of annealing depending on the primers used, and an extension at 72 • C for 45 s. The final extension is performed at 72 • C for 3 min. PCR amplification of 16S rRNA gene fragments was done using the two opposing primers fD1 and rD1 and the primer pairs nodCFn-nodCI and nodA1f-nodAb1r were used for amplification of the nodC and nodA gene (Chaintreuil et al., 2001;Laguerre et al., 2001).
The PCR products and their concentration were assessed by electrophoresis of 6 µl of product on a 1% agarose gel and staining with ethidium bromide. A molecular mass marker (Promega 1-kb ladder) was included to estimate the length of the amplicons.
The amplification products were purified using a Qiaquick PCR purification kit (Qiagen) and sequenced using the same primers as for PCR amplification using the dideoxynucleotide chain-termination method with fluorescent ddNTPs (Applied Biosystems), with ABI Prism Dye Chemistry, and analyzed with a 3130xl automatic sequencer at the sequencing facilities of UATRS (Unités d'Appui Technique à la Recherche Scientifique, CNRST, Rabat, Morocco) and the CSIC sequencing facilities in Zaidin experimental station, CSIC, Granada, Spain. The sequences obtained were compared with those from GenBank using the BLASTN program (Altschul et al., 1990).

Sequences Analysis and Phylogeny
Sequences were aligned using MEGA 7 software (Kumar et al., 2016), and the distances calculated according to Kimura's two-parameter model (Kimura, 1980) were used to infer phylogenetic trees with the neighbor-joining analysis (Saitou and Nei, 1987) with MEGA 7 software (Kumar et al., 2016). Phylogenetic trees were subjected to 1000 bootstrap replications, and preferred topologies were plotted.

Phenotypic Characterization
The analysis was pursued by determining some phenotypic properties of the isolated strains. Tests were performed in broth tubes or agar plates inoculated with an exponentially growing liquid culture. The tolerance of the rhizobial isolates to high temperatures was tested on TY broth medium (Beringer, 1974) at 30, 35, 40, and 45 • C. The ability of isolates to grow in acidic or basic media was determined on YMA Petri dishes for which the pH has been adjusted and buffered to 5.0, 5.5, 6.0, 7.0, 8.0, 9.0, 9.5, or 10.0, as described by Zerhari et al. (2000). The salt tolerance of the isolates was tested at 0, 86, 171, 342, 510, and 685 mM.
Utilization of 12 amino acids as sole nitrogen source was investigated on a modified YMB medium in which the yeast extract was replaced by the amino acid to test. The amino acids were sterilized separately by filtration and added to the medium at a concentration of at 1% (w/v). Carbohydrate assimilation was carried out on solid YEM containing 0.005% YEM and 0.4% (w/v) of glucose, fructose, sucrose, maltose, lactose or starch. The carbohydrates were sterilized separately by filtration.
To evaluate oxidase activity, a loopful of bacterial mass from the colonies of the isolates was spread out over cellulose paper soaked in a solution containing 1% N,N-dimethyl-pphenylenediamine oxalate 98%. The appearance of blue color is indicative of oxidase activity.

PGPR Characteristics
We used the PVK solid medium containing 2.5 g L −1 of (Ca 3 PO 4 ) 2 to determine the ability of the isolates to solubilize phosphates. To this purpose, 5 µl of the bacterial suspension was plated onto PVK medium in Petri dishes and incubated at 28 • C for 10 days. The size of the halos that appeared around the colonies and the size of the colonies were measured to estimate the ability of the isolates to solubilize phosphate. Siderophoreproducing isolates were identified by the development of orange halos around the colonies according to Lakshmanan et al. (2015). To assess the ability of the isolates to produce indole acetic acid (IAA) and IAA-related compounds, a volume of 5 µl of the bacterial suspension was placed on the solid YEM medium supplemented with 0.5 g L −1 tryptophan, and the plates were then incubated at 28 • C for 10 days. The production of IAA was detected using the Salkowski reagent followed by the appearance of a pink halo around the colonies after incubation for 30 to 60 min at 28 • C.

Strain Isolation, REP-PCR, and ARDRA Analysis
Thirty-six bacterial cells were isolated from root nodules of L. cosentinii (20) and L. luteus (16) grown in Maamora cork oak forest soils. Among the strains, 28 were able to re-nodulate their original host under axenic conditions using Gibson's tubes. REP-PCR using primers REP1R-I and REP2-I grouped the strains into six clusters and some single strains (Supplementary Figure 1). The REP PCR technique is employed for analysis of the genetic diversity of rhizobia (Laguerre et al., 1997;Menna et al., 2009).
Eight strains representing the different groups were randomly selected and used for in vitro ARDRA, based on 1,450 nucleotides of rrs gene, with the restriction enzymes cfoI, HaeIII, and MspI, which distributed the strains in three ribotypes (Supplementary Table 3). Furthermore, the in silico ARDRA of the eight representative strains, together with other 45 bradyrhizobial species and the outgroup R. gallicum R602, using 30 restriction enzymes, including MspI, HaeIII, and RsaI, showed that the strains from the two lupines clustered in three main groups, with B. lupini-B. canariense and B. cytisi (Supplementary Figure 2).
The strains Lcos102 and Llut5 have similarity values of 99.79 and 99.86% with B. canariense BTA-1 T . They have similarities of 99.08-99.85% with the North African strains tested (Supplementary Table 4). The strains Lcos8.1 and Llut8 share similarities of 99.64 and 99.79% with B. cytisi CTAW11 T and share also similarities of 98.93-99.85% with the strains isolated from the Algerian and Tunisian lupines.
The phylogenetic tree based on the 16S rRNA sequences showed that the strains Lut6, Lcos6, and Lcos7.2 clustered together with different strains isolated from North African lupines, in Algeria and Tunisia in group I, along with B. lupini USDA 3051 T . The strains Lcos102 and Llut5 clustered in group II with B. canariense LMG 22265 T , whereas strains Lcos8.1 and Llut8 clustered in group III with B. cytisi CTAW T and some strains isolated from Algerian and Tunisian lupines. The strain Llut4 grouped further into the cluster composed of groups 1 and 2 (Figure 1). The rrs sequences of strains Lcos6, Lcos7.2, and Llut6 are 100% similar, and thus, we selected one strain (Llcos6) to represent the three. The strains retained as representatives for further molecular analyses are Lcos6, Llut4, Llut5, Lcos102, Lcos8.1, and Llut8.

Multilocus Sequence Analysis
Five housekeeping genes, dnaK, glnII, gyrB, recA, and rpoB of the six representative strains were analyzed in this study. The recA and glnII sequences were analyzed and compared with other strains isolated from North African lupines root nodules. All the sequence accession numbers are reported on the different figures. We could not compare our sequences with the dnaK, gyrB, and rpoB sequences of Algerian and Tunisian strains since they are not published.
The sequence analyses showed that the recA genes of strains Llut5 (MK346962) and Lcos6 (MK346959) are 100% similar and have 99.80% of similarity with B. lupini 3051 T . The two strains share 99.58-100% similarity with strain LaO10, LaO17, and Lat46 isolated from L. angustifolius in Algeria (Mellal et al., 2019), and 99.59% with strain LmiT3 isolated from the root nodules of L. micranthus in Tunisia (Msaddak et al., 2017). The strains Lcos8.1 (MW527003), Llut8 (MK346964), and Lcos102 (MK346960) share 99.23 to 99.60% and have similarity percentages of 94.43 to 95.23% with Llut4 (MK346961). They have percentages of similarities ranging from 93.51 to 94.30 with the different Algerian and Tunisian strains. The phylogenetic tree based on the recA sequences (Supplementary Figure 3) shows that strains Lcos6 and Llut5 clustered with strains nodulating L. micranthus and L. angustifolius in Algeria all together with the type strain B. lupini USDA3051. The strain Llut4 was closer to B. cytisi LMG25866 T with which it constituted a group containing B. japonicum LMG6138 T and B. diazoefficiens USDA110 T . The strains Llut8, Lcos8.1, and Lcos102 were regrouped in a different cluster with B. ganzhouense RITF806 T .
The glnII neighbor-joining phylogenetic tree (Supplementary Figure 3) shows that strains Lcos6 and Llut5 are close to B. lupini 3051 T with which they were regrouped along with strains LaO10, LaO17, and LaT46 isolated from L. angustifolius in Algeria and strain LmiT3 isolated from L. micranthus in Tunisia.
Similar results were obtained with the gyrB, dnaK, and rpoB sequences (Supplementary Figures 4-8). However, we could not compare our sequences with those of other North African lupines' microsymbionts.
The concatenation of the glnII, gyrB, and recA genes showed that the strains Llut8, Llcos81, and Llcos102 have identity percentages of 95.8, 95.73, and 95.59% with B. ganzhouense RITF806 T and share similarities of 99.43-99.57%. The strains Lcos6 and Llut5 share 99.70 similarities and have identity percentages of 99.43 and 99.50% with B. lupini LMG28514 T , respectively. The strain Llut4 has low similarities with all the strains tested and shares a maximum identity of 94% with B. cytisi LMG25866 T (Supplementary Table 4). The phylogenetic tree based on the concatenation of the three genes (Figure 2) revealed also that B. cytisi LMG25866 T is the closest parent to strain Llut4, whereas the strains Lut5 and Lcos6 clustered with B. lupini LMG28514 T , and strains Llut8, Llcos81, and Llcos102 were regrouped with B. ganzhouense RITF806 T .

nodA and nodC Symbiotic Genes Sequences Analysis
The nodC sequences of the strains Lut5, Llut8, and Lcos6 share 99.9-100% similarity. The nodC-based phylogenetic tree (Figure 3) showed that they were regrouped in a same cluster with strain LAM15, isolated from the nodules of L. angustifolius in Algeria (Mellal et al., 2019), with which they have a similarity percentage of 92.1%, within the genistearum symbiovar. The strain Llut4 has a 95.9% percentage similarity with the three other representatives and was regrouped with B. cytisi LMG25866T with which it shares 97% similarity.
The analysis of the nodA gene sequences showed also that our strains were very close and share 99.6 to 100% similarity. The phylogenetic tree based on the nodA sequences (Supplementary Figure 9) showed also that they formed a single sub-cluster in the symbiovar genistearum. The phylogenetic analysis of the two symbiotic genes' concatenated sequences also produced similar results (Supplementary Figure 10).
The Moroccan as well as the Algerian and the Tunisian strains were clustered together in the same clade with B. lupini, B. canariense, B. cytisi, and B. rifense, members of the genistearum symbiovar.

Phenotypic Tests
All the isolated strains were slow growing, as their colonies appeared on YEM solid medium after 7 days incubation time at 28 • C. The phenotypic results reported in Supplementary Table 6 show that the strains grow in media with different pH values ranging from 6 to 8. All strains nodulating L. cosentinii use trehalose as sole carbon source; 93% assimilate fructose and myoinositol, and 64% grew in the presence of arabinose, while only 7% can assimilate maltose. The latter cannot be assimilated by strains nodulating L. luteus. The results also showed that all strains of L. luteus assimilate arabinose and fructose, while 79% of these strains use myo-inositol and 36% develop in the presence of trehalose. None of the two lupine species microsymbionts is capable to use lactose, starch, or carboxymethyl-cellulose (CMC).
The strains nodulating lupines in the Maâmora forest are able to use a wide range of amino acids as sole nitrogen source. All the strains nodulating L. luteus assimilate 10 amino acids, arginine, tyrosine, proline, asparagine, histidine, valine, serine, tryptophan, leucine, and phenylalanine, while 93% of these strains use methionine, and 64% use alanine, whereas only 14% are able to use glycine, lysine, and aspartic acid. On the other hand, all strains that nodulate L. cosentinii are able to use tyrosine, serine, and phenylalanine as the sole source of nitrogen, while 93% assimilate arginine, proline, valine, tryptophan, and leucine, whereas 86% can assimilate asparagine and histidine while 57% are able to assimilate methionine. Glycine, lysine, and aspartic acid are used by only 21% of these strains.
We also tested some enzymatic activities in the strains studied. Thus, 18% of strains nodulating L. cosentinii and 9% among strains of L. luteus produce gelatinase, and all the strains isolated from the two lupine species produce catalase.
Tolerance to NaCl varies considerably between strains nodulating the two species. The results show that all strains tested can grow in the presence of 85 mM NaCl, while 21% of L. cosentinii strains and 14% of L. luteus strains develop in salt concentrations between 170 mM and 680 mM, while only 7% of the strains nodulating L. cosentinii and 7% of those of L. luteus continue to develop at 854 mM NaCl.
Strains react differently when grown in YEM media with different pH values. All the strains nodulating the two species of lupine have the ability to grow in a slightly acidic to neutral pH (6 to 7). All the strains nodulating L. luteus grow in a slightly alkaline pH (pH 7.5) and 93% grow at pH 8. They are sensitive to more alkaline pH, with only 50% of the strains able to grow at pH 8.5 and 14% at pH 9. The sensitivity of these strains was also noted with respect to acidic pH; thus, no strain tolerated pH 5; however, 14% of these strains were able to grow at pH 5.5. The results also show that the majority of strains nodulating L. cosentinii (86%) can develop at pH 7.5 and that with the increase in pH, the number of strains that manage to grow there decreases, 36% grow at pH 8, 29% at pH 8.5, and only 21% at pH 9. In acidic environments, 21% of the strains manage to develop at pH 5.5 and 7% at pH 5.
Among all the heavy metals tested, barium, manganese chloride, manganese sulfate, aluminum, and cadmium (at 100 µg/ml) were the least harmful for the growth of the strains nodulating L. luteus with 79-93% that resist. On the other hand, the presence of cobalt and copper inhibited the growth of half of the strains. The heavy metals most toxic to strains of L. luteus are mercury, magnesium, zinc, and lead, and only 21% could grow in their presence. The results showed also that 50% of strains nodulating L. cosentinii develop in the presence of manganese chloride and 43% are resistant to aluminum. On the other hand, in presence of barium, cobalt, and manganese sulfate, 36% of the strains were able to resist, and 28% of the strains were able to develop in the presence of copper. Mercury, magnesium, zinc at 500 µg/ml, and lead at 1,000 µg/ml were found to be more toxic with only 21% of strains that were able to tolerate their presence.
The tolerance to antibiotic is one of the defining criteria in the analysis of rhizobial biodiversity. In this work, we tested 11 antibiotics on the different strains isolated from the two lupines and the results showed that there is variability in tolerance. Thus, all the strains nodulating L. luteus were resistant to several families of antibiotics: beta-lactams (ampicillin and penicillin), quinolones (nalidixic acid), glucopeptides (vancomycin), aminoglycosides (spectinomycin), gramicidins (bacitracin), and trimethoprims. They grew also in presence of tetracyclines and aminoglycosides (gentamycins and kanamycin) with high percentages between 70 and 86%. The chloramphenicol affects the growth of these bacteria since only 36% were FIGURE 2 | Maximum likelihood phylogenetic tree based on concatenated gyrB + recA + glnII sequences (1,430 nt) of strains from nodules of L. cosentinii and L. luteus and phylogenetically related species within the genus Bradyrhizobium. The analysis was based on 1,405 nucleotides. Isolates are denoted in bold. Bootstrap values are indicated as percentages derived from 1,000 replications. Values lower than 50 are not shown. Bar, 2 nucleotides substitution per 100 nucleotides. The tree is rooted with Rhizobium pisi DSM 30132 T . able to grow in its presence, at the concentrations used. The results also showed that 43% of the strains nodulating L. cosentinii tolerated the majority of antibiotics: nalidixic acid, vancomycin, spectinomycin, bacitracin, trimethoprim, and beta-lactams, and 36% developed in the presence of chloramphenicol. These strains are sensitive to aminoglycosides, since only 21% are resistant to kanamycin and 14% are resistant to gentamycin.
Temperature also has a differential effect on strains' growth. Thus, all the strains grew at 30 • C, whereas at 35 • C, the strains showed a fluctuation in their responses, with 79% of the strains nodulating L. cosentinii, and 71% for the strains nodulating L. luteus, which were able to grow at this temperature. No strain can grow above 40 • C.

PGPR Growth Promoting Activities Test Results
Out of the 28 strains tested, six strains produced the IAA, four strains from L. cosentinii and two from L. luteus. Only three strains were able to solubilize the inorganic phosphates, two strains from L. cosentinii and one strain among L. luteus microsymbionts. Two strains were hence very effective with a solubilization index of 275% for lut5 and 375% for cos10.2. Six strains produced siderophores, four strains from L. cosentinii and two from L. luteus. The two strains Lcos10.2 and Llut5 have the capacity to solubilize phosphates and produce siderophores and IAA.

DISCUSSION
In this work, we aimed to characterize and identify some bacterial strains isolated from the root nodules of two lupines grown in different soils of the Maamora forest, in the vicinity of Rabat.
We started the molecular characterization by the analysis of the eight representative strains' rrs gene sequences, and their phylogeny revealed that they are more related to B. cytisi and B. lupini. The rrs gene was considered as a preliminary genetic marker to situate a strain taxonomic position at the genus level, although it does not provide clear phylogenetic diversity for species of the Bradyrhizobium genus (Martinez-Romero and Ormeño-Orrillo, 2019). To this purpose, an MLSA of different housekeeping genes as additional phylogenetic markers is frequently used for accurate identification of Bradyrhizobium species and strains (Tampakaki et al., 2017;Wójcik and Kalita, 2019).
Six strains were then selected from the eight representatives, as the rrs sequences of strains Llut6 and Lcos7.2 were identical (100%) to strain Lcos6 rrs sequence. The phylogenies of the three individual genes glnII, gyrB, and recA showed that strains Lcos6 and Llut5 were closely related to B. lupini LMG 28514 T . In fact, we sequenced also the dnaK and rpoB genes of most strains, but we failed to amplify them in some. Furthermore, data concerning these two genes in B. lupini are not available in the databases, and the dnaK (JACJNR010000003.1) and rpoB (JACJNR010000001.1) sequences of the strain B. lupini DSM 30140 extracted from the whole genome sequence available in the NCBI database are more related to B. japonicum LMG6138 T (Supplementary Figures 5,  6). Consequently, we preferred to continue with only three genes.
The concatenated housekeeping genes' sequences phylogeny confirmed that the closest parent of the strains Lcos6 and Llut5 is B. lupini LMG 28514 T , whereas strain Llut4 was more related to B. cytisi LMG 25866 T . The three strains Lcos8.1, Lcos102 and Llut8 were regrouped in a single cluster more related to B. ganzhouense RITF806 T , isolated from the rot nodules of Acacia melanoxylon in China (Lu et al., 2014). The two strains Lcos8.1 and Lcos102 were isolated from L. cosentinii, whereas strain Llut8 was isolated from L. luteus. None of the three strains grew in YEMA medium in the presence of 1% NaCl. The strain Lcos6 grows between pH 6 and 7.5, while cos8.1 grows between pH 6 and 8, and the Lut8 strain grows in pH 6-8.5, but do not grow at 40 • C. They use glucose, fructose, mannose, and sucrose as sole carbon sources. Inversely, the type strain of B. ganzhouense does not assimilate sucrose; it grows at pH 5-12 and develops in the presence of 3% NaCl (Lu et al., 2014). Their low percentages of similarity with their closest parent suggest that they may constitute a putative new genospecies in the genus Bradyrhizobium.
Plant species of the genus Lupinus were considered as nodulated exclusively by Bradyrhizobium spp. for a long time until Ardley et al. (2012) reported the isolation and identification of three symbiotic Microvirga species from L. texensis and Listia angolensis root nodules. More recently, lupines were found to be nodulated by Microvirga strains in Tunisian, Moroccan, and American soils (Beligala et al., 2017;Msaddak et al., 2017Msaddak et al., , 2019Rejili et al., 2019;Missbah El Idrissi et al., 2020). Msaddak et al. (2017) suggested that lupines would prefer Bradyrhizobium spp. strains as microsymbionts in acidic soils and Microvirga spp. in the alkaline soils. Lamrabet et al. (2020) also reported that in alkaline soils, some Genistae such as Retama spp. would prefer members of Microvirga as microsymbiont instead of Bradyrhizobium, although Tounsi-Hammami et al. (2019) found that strains isolated from L. albus in alkaline soils were affiliated to Rhizobium and Neorhizobium. Furthermore, Durán et al. (2014) reported the nodulation of the endemic Lupinus mariaejosephae by Bradyrhizobium valentinum in alkaline-limed soils in Eastern Spain.
The phenotypic analysis showed that all the strains assimilate different carbon sources including glucose, mannose, galactose, and sucrose, while no strain used starch or lactose; this result agrees with previously reported results (Rome et al., 1996;Missbah El Idrissi et al., 2020). Bradyrhizobium canariense BTA-1T does not use sucrose or lactose either (Vinuesa et al., 2005).
There was diversity in the utilization of amino acids, but all the strains were able to assimilate tyrosine and phenylalanine as sole nitrogen source, which contradicts results of Miller and Pepper (1988) who reported that the asparagine was the most preferred amino acid for the growth of rhizobia isolated from lupines in the desert of Mexico. Vinuesa et al. (2005) reported also that B. canariense and B. yuanmingense do not use glycine. Some strains possess the gelatinase activity, an enzyme that hydrolyzes collagen into amino acids or peptides, which allows them to persist in in the rhizosphere while waiting for the plant partner.
We analyzed the effect of environmental factors including salinity, drought, acidity, and alkalinity as well as the behavior of strains toward heavy metals and antibiotics, factors that compromise survival, growth, and nitrogen fixing ability of rhizobia strains (Zahran et al., 1994). We noticed a low tolerance to salinity in the majority of the strains as reported by Vinuesa et al. (2005), but some strains tolerated salt concentrations as high as 854 mM. Many studies have reported very interesting resistance profiles in strains isolated from lupines in Egypt that can tolerate up to 1,700 mM NaCl (Zahran et al., 1994). However, the limits of salinity tolerance between rhizobia can vary considerably from one species to another (Elsheikh and Wood, 1989), and even between strains of the same species (Boukhatem et al., 2012). Peix et al. (2015) found that among 11 type strains of Bradyrhizobium species tested, only B. betae, B. japonicum, and B. ganzhouense grew at 171 mM of NaCl. In fact, in saline environments, the symbiosis between rhizobia and lupines depends not only on the bacteria but on the host plant also. Some strains were able to grow at pH 9, while few strains tolerated pH 5. Our strains tolerate more alkaline than acidic pHs. The type strains of B. cytisi and B. rifense are reported to grow at pH 4.5 (Peix et al., 2015). All the isolates grow up to 35 • C, while no strain grew at 40 • C, which is very common in bradyrhizobia (Peix et al., 2015). Bradyrhizobium sp. strains isolated from L. luteus in Tunisia were not able to grow at 37 • C, pH 12, and 1% NaCl (Msaddak et al., 2017). Other studies have shown that certain rhizobial strains can tolerate higher temperatures, but it has been suggested that rhizobia may be protected against different soil constraints by living in particular niches in the rhizosphere and inside nodules (Boukhatem et al., 2012).
L. luteus nodulating strains are more resistant to antibiotics (ATB) and heavy metals than strains isolated from L. cosentinii, among which only 50% are able to resist to some antibiotics. The resistance level varies between strains, depending on the type of ATB and its mechanism of action. Similar results were obtained by (Msaddak et al., 2017).
The majority of our strains showed modest activities of phosphate solubilization and production of siderophores or auxins. Sarkar and Laha (2013) reported that 80% of rhizobia have the capacity to produce IAA and IAA-related compounds. However; siderophore production is usually low and phosphate solubilization is variable in rhizobia (Lamrabet et al., 2020). Few rhizobia may possess all the plant growth promoting (PGP) activities and the plants need hence to interact with different beneficial strains in their rhizosphere. To sustain their development, the most healthy plants attract the best PGPRs (Backer et al., 2018;Gouda et al., 2018).
It is evident that environmental factors such as salinity, drought, acidity, alkalinity, heavy metals, and the presence of antibiotics compromise the survival, growth, and the ability to fix nitrogen of strains of rhizobia (De la Peña and Pueyo, 2012), which shows the interest of knowing their phenotypic and biochemical characteristics in advance before being introduced and inoculated to the plants in the field, to ensure their survival, their adaptability, and their competitiveness (Boukhatem et al., 2012).
The nodC symbiotic gene (encoding the N-acetylglucosaminyl transferase protein) is essential for nodulation of compatible host legumes (Laguerre et al., 2001). However, although its phylogeny is not congruent with the core genes' phylogeny (Andrews and Andrews, 2017), it is used to characterize the new isolates at the symbiovar level (Rogel et al., 2011).
Analysis of the individual and concatenated nodA and nodC gene sequences and their phylogeny showed that the strains Llut5, Lcos6, Llut4, and Llut8 are members of the genistearum symbiovar.
All the strains are able to nodulate R. monosperma, L. luteus, L. albus, L. cosentinii, C. albidus, and C. monspessulanus but not G. max or Phaseolus vulgaris. It is known that the Genisteae microsymbionts do not nodulate soybean (G. max or G. soya: Phaseoleae tribe). However, they are able to nodulate different genera and species of the tribe Genisteae (Vinuesa et al., 2005). The fact that the strains are able to nodulate and fix nitrogen with different Genisteae is an advantage to members of this tribe. It is important to note that we observed no significant visible differences between the two lupine plants inoculated with the Lcos strains (isolated from L. cosentinii) or the Llut strains (isolated from L. luteus). All the plants were in good health as shown by the leaves' green color after 3 months. This represents an advantage for these legume plants, because they can establish a nitrogen fixing symbiosis with different strains with different PGP potentialities, which contributes to their resilience in the case of any edaphoclimatic changes.

CONCLUSION
In this work, we report on the characterization of some strains isolated from the root nodules of two lupines grown in the acidic soils of the Maamora cork forest, in Morocco.
We found that L. cosentinii is nodulated by members of the genus Bradyrhizobum, whereas the species has been reported as nodulated by Microvirga sp. in alkaline soils of the same forest. As already reported, the selection of lupine microsymbionts in the south Mediterranean would depend on the type and pH of the soil. Some Bradyrhizobium spp. strains newly isolated from L. luteus and L. cosentinii may constitute a new genospecies, and this needs confirmation by more sophisticated research such as whole genome sequencing or DNA-DNA hybridization.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI GenBank [accession: MW526996-MW527060, MW494669-MW494676, MW547016, and MW565871-MW565872].

AUTHOR CONTRIBUTIONS
MM participated in the conception and design of the experimentations and discussion of the results. OB participated in the analysis of genomic data and phylogenies. SE, HL, and SA participated in the isolation and characterization of isolates. EB participated in the sequencing experiments. HA participated in the conception of the project. All authors contributed to the article and approved the submitted version.