Deciphering the Factors for Nodulation and Symbiosis of Mesorhizobium Associated with Cicer arietinum in Northwest India

: The compatibility between rhizobia and legumes for nitrogen ‐ fixing nodules and the stages of root hair curling, formation of infection thread, and nodulation initiation have been vitally studied, but the factors for the sustainable root surface colonization and efficient symbiosis within chickpea and rhizobia have been poorly investigated. Hence, we aimed to analyze phenotypic properties and phylogenetic relationships of root ‐ nodule bacteria associated with chickpea ( Cicer arietinum ) in the north ‐ west Indo Gangetic Plains (NW ‐ IGP) region of Uttar Pradesh, India. In this study, 54 isolates were recovered from five agricultural locations. Strains exhibited high exopolysaccharide production and were capable of survival at 15–42 °C. Assays for phosphate solubilization, catalase, oxidase, Indole acetic acid (IAA) production, and 1 ‐ aminocyclopropane ‐ 1 ‐ carboxylic acid (ACC) deaminase activity revealed that all the tested isolates possessed plant growth ‐ promoting potential. Metabolic profiling using Biolog plates indicated that patterns of substrate utilization differed considerably among isolates. A biofilm formation assay showed that isolates displayed a nearly four ‐ fold range in their capacity for biofilm development. Inoculation experiments indicated that all isolates formed nodules on chickpea, but they exhibited more than a two ‐ fold range in symbiotic efficiency. No nodules were observed on four other legumes ( Phaseolus vulgaris , Pisum sativum , Lens culinaris , and Vigna mungo ). Concatenated sequences from six loci ( gap , edD , glnD , gnD , rpoB , and nodC ) supported the assignment of all isolates to the species Mesorhizobium ciceri , with strain M. ciceri Ca181 as their closest relative


Introduction
For years, plant-associated bacteria have been speculated to be the organisms capable of promoting growth and/or suppressing diseases when present in the rhizosphere and as endophytes within healthy plant tissues [1][2]. Leguminous plants have the capability to associate with rhizobia by the nodulation process and fix atmospheric nitrogen (N) in the root nodule. This can help plants in soil that is low in nitrogen (N) and condition the soil itself, if other factors are favorable for growth. Bacteria have typically been associated with or adhere around the legumes plant roots and act as a plant growth promoter by several secretory or substrate capabilities, such as phosphate, ACC deaminase, siderophore, IAA, catalase, oxidase, and NH3 production [3][4]. Moreover, they also act as endophytic or nodule-forming bacteria (found within the tissues of the plant) as well as saprophytic bacteria (found free-living in the soil).
Rhizobia may colonize soil that envelops the roots [5] or adhere to the root surface [6], which is directly influenced by their attachments on their desired surface. Biofilm formation by bacteria offers conductive sites for an appropriate, effective, and reproductive environment for the bacteria to adhere to the surface. Biofilm thus provides the space for growth, a protective degree of homeostasis, and allows the bacteria to overcome biotic and abiotic stresses-this is done by a complex extrapolymeric substance (EPS) matrix sheath [7][8][9]. During infection, thread formation, and root colonization, EPS plays an important role in the infection and the active nodules [8]. Bacteria adheres to the legume roots with the help of EPS and can fix the nitrogen for plant growth, act as biocontrol against pathogens, produce phytohormone, and mobilize the nutrients to enriched the soil and the environment [10][11]. Rhizobia have been used in agricultural practices mainly for nitrogen fixation and plant growth promotion (PGP) [12][13][14] due to their wide distribution. The first criterion for a Rhizobium strain to be used in legume inocula is that it must be highly effective when fixing nitrogen [15].
Chickpea is one of the major pulse crops throughout the world and ranks second (The Food and Agriculture Organization of United Nation) amongst food legumes in terms of world production and is cultivated on a large scale in arid and semi-arid environments [16][17]. Chickpea is the most important legume in the Indo Gangetic Plain (IGP), and its cultivation is completely dependent on rhizobia, which have an impact on agriculture as well as the environment. Currently, chickpea rhizobia are included in the genus Mesorhizobium (gram-negative bacterium), with two species that form functional nodules after interacting with chickpea roots [18][19]. These species are described as specific microsymbionts-M. ciceri [20] and M. mediterraneum [21]. When not in symbiosis, they are part of the habitat and behave as saprophytes. M. ciceri cells sense the root exudates in the soil and invade after adhering to the root hairs of C. arietinum, which are dependent on the exopolysaccharides of rhizobia. Which in turn, act as signal molecules for the establishment of symbiosis between rhizobia and legumes [22][23]. After invading the root, the formation of specific root organs, called nodules, takes place. During the maturation of the nodules, invading rhizobial cells differentiate into endo-symbiotic structures called bacteroides. Fixed nitrogen (NH4 + ) is provided to plants by bacteroides. It is hypothesized that, in return, the plants deliver organic acids and several carbon and energy sources to the bacteroides as root exudates [24][25][26]. The overall procedure is completely based on saprophytic survival and adherence during chickpea cropping.
The absence of a high population of potential nitrogen-fixing rhizobia in the soil and their sustainability in symbiosis has been one of the main limiting factors for legume production. Plant growth promotion, resistance to stress, and utilization of various carbon sources support the rhizobia in their saprophytic survival (while they are in a free-living state), effective nodulation with legumes, and symbiotic efficiency. Moreover, the ability to survive against stress and the formation of biofilm can shape the abundance and efficiency of rhizobia. Hence, we aimed to investigate the major factors required for efficient symbiosis and PGP for mesorhizobia, associated symbiotically with C. arietinum. Further, sequencing of six gene loci was done to characterize the diversity and relationships of Mesorhizobium strains associated with chickpea in the IGP region of India.

Sampling Sites, Isolation, and Culture Condition
This study was carried out on chickpea nodules from 32 different representatives from five semiarid or sub-humid alkaline field sites of the NW-IGP region Uttar Pradesh, India ( Figure S1). Isolation and purification of rhizobia from root nodules was done on yeast extract-mannitol agar (YMA) by using the standard procedure [27]. All procured rhizobia were incubated on YMA slants at 28 °C and maintained at 4 °C for routine use and in 20% (w/v) glycerol at −80 °C for long-term storage [28]. Standard culture of M. ciceri strains IC-2058 and IC-2018 were taken from the microbial genomics laboratory of Indian council of Agriculture research-National Bureau of Agriculturally Important Microorganisms (ICAR-NBAIM) Mau, India.

Validation, Nodulation, and Host-Specificity Test of Chickpea Nodule Rhizobia
Studied isolates were validated as rhizobia by testing the ability to induce root nodules on chickpea according to previous reported methods [29]. Briefly, 100 μL log phase culture of tested isolates was taken and inoculated on to pre-surface sterilized chickpea seeds (Avarodhi variety). Further, the sterilized clay pots (30 cm high 20 cm diameter) were filled with equal amount of autoclaved soil (approx. 4.0 kg). After that, inoculated seeds were sown in clay pots for growth. The ideal condition for plant growth was adjusted (28/19 °C (day/night) with a 12-h photoperiod). After 40 days, plants were uprooted gently from the pots for observation of various parameters such as nodule number (NN), shoot dry weight (SDW), chlorophyll content, and symbiotic efficiency (SE) according to Gibson A.H. [30]. Nitrogen supplemented (T2) plant was treated as a control while without N (T1) acted as a negative control. The whole experiment was performed in 5 replicates.
Host specificity was examined by cross nodulation test as described by Nandasena et al. [31]. Briefly, all the tested rhizobia were inoculated in Phaseolus vulgaris, Pisum sativum, Lens culinaris, and Vigna mungo separately in pots, as described for nodulation assay in this study. During the flowering stage, plants were uprooted gently and the root nodules were examined.

PGP Traits
Phosphate solubilization index (S.I.) of the tested strains was determined using Pikovaskya agar plates spot method [32], calculated as the ratio of halo zone diameter to colony diameter for positive isolates [33]. Qualitative and quantitative estimation of IAA production was estimated as described by Khanmna et al. [34]. Ammonia (NH3) production of each isolate was evaluated in peptone water. Briefly, 100 μL of log phase cultures (10 6 cells mL −1 ) were inoculated in 15 ml nutrient broth (NB from Hi-Media Co.) and incubated at 125 rpm and 28 ± 2 °C in rotator shaker for 4-5 days. Development of deep yellow to brown color by addition of Nessler's reagent (Hi-Media Co.) (0.5 mL) was noted as a positive test for NH3 production [34]. Catalase and oxidase production tests were carried out according to Smibert et al. [35]. For ACC deaminase, minimal medium (ACC as sole source of nitrogen) was inoculated with 100 μL log phase culture, and growth was observed at different time intervals through the absorbance measurement (at 600 nm), as well as viable count by serial dilutionspread plate method. Multiplication capabilities of inoculated strains were considered as ACC deaminase positive. Quantification of ACC deaminase was carried out according to Glick et al. [36].
All morpho-phenotypic and PGP tests were performed in five replicates at three different time intervals for their reproducibility and consistency.

Biofilm Assay
Qualitative and quantitative estimation of biofilm formation ability of all the tested rhizobia were evaluated as described by Mirani J. A. and Jamil N. [37] with partial modification. Briefly, log phase culture of isolates was pelleted out, washed 2 times with ddH2O, and resuspended into sterile ddH2O to adjust the optical density (OD) 0.2 at A600 with spectrophotometer. Further, 100 μL of dissolved culture was inoculated into the small size culture vial that contained 250 μL of M63 minimal medium supplemented with 0.2% glucose (Hi Media Co. In,), 1 mM magnesium sulfate (Hi Media Co. In.), and 0.5% casamino acids (Hi Media Co. In.). Further, the cover slip was placed vertically and then incubated for 48 h at 28 °C. After that, the cover slip was put carefully in 0.1% crystal violet (CV) for 30 min to stain the developed biofilm. Then, it was solubilized in 1 ml 30% glacial acetic acid and the absorbance was recorded (Shimadzu Co. Ltd. In.) at 563 nm wave length by using at least five replicates for each assay to avoid the deviation and error.

Biolog Carbon Source Utilization and Numerical Analysis
Based on the SE test, morpho-phenotypic characteristics, and biofilm assay, 16 most prominent isolates were selected for further characterization. Isolates were cultivated separately in YEM broth, harvested at log phase, washed twice with autoclave ddH2O, and maintained the OD600 1.0 in autoclave ddH2O for Biolog plate inoculation. Individual strains were tested for sole carbon source utilization on Biolog GN microtiter plates (Microlog2, Version4.2, Biolog Inc., Hayward, CA, USA). GN plate wells were inoculated with 150 μl of cell suspension of selected isolates and incubated at 28 °C for 24 h. BIOLOG GN microplates used in this study contained 95 wells each with a separate sole C source, and a control well without a C source. Complete information regarding the plates has been given by Garland and Mills [38]. The average well color development (AWCD) was calculated for each microplate. AWCD analysis allows a comparison of plates that have achieved the same degree of color development regardless of differences in cell inoculum density that may cause differences in rates of color development [39]. OD readings were first truncated to lie in the range {0, 2}, since values below 0 are clearly erroneous and values above 2 have been shown to be dominated by measurement error.
The AWCD for each microplate was calculated by subtracting the control well optical density from the substrate well OD (blanked substrate wells), setting any resultant blanked substrate wells with negative values to zero and taking the mean of the 95 blanked substrate wells. The mean of the AWCD for each set of triplicate plates was calculated. If OD represents the corrected OD for well, i of replicate j at time t, then the AWCD for replicate j at time t is given as: The standardized OD values are then given by: The readings obtained from the selected Biolog microplates were analyzed and compared by principal component analysis (PCA), using SPSS v16. Interpretation of the principal components was based on significant factor loading of the individual substrates on each of the principal components [39].

Molecular Identification, Sequence Typing, Genetic Differentiation, and Gene Flow of (gap, edD, gnD, glnD, rpoB, and nodC)
Genomic DNA of selected strains were isolated through Wizard ® Genomic DNA Purification Kit (Promega) and quality was measured in NanoDrop microvolume spectrophotometers. Gene amplicons of 16S rRNA, glyceraldehyde3 -phosphate dehydrogenase (gap), phosphogluconate dehydratase (edD), 6-phosphogluconate dehydrogenase (gnD) and protein-PII uridylyl transferase (glnD) genes, ribosomal polymerase B subunit (rpoB), and nodulation gene (nodC) were amplified by primer pairs in in an ABI. PCR system. Primer sequences, locus tag and locus location of selected primers were mentioned in supplementary Table S2. also [40] Amplified genes were examined by horizontal electrophoresis in 1.5% agarose with 1 μL aliquots of PCR product and further purified by gel elution method (Nucleopore gel elution kit). Amplified gene products were sequenced in automated sequencer (model ABI 3130xl), using ABI cycle sequencing kit version 3.1. Further, BLAST searches were done to identify similar sequences in the NCBI database. Phylogeny, robustness of the tree topology, and distance calculation of 16S rRNA, gap, edD, gnD, glnD, rpoB, and nodC gene sequences were carried out by using the minimum evolution bootstrap phylogeny method and bootstrapping algorithms confined (1000 replication) in the MEGA version 6.0 [41]. Further, sequences were concatenated and aligned using ClustalW with the manually concatenated sequences of the same housekeeping genes from type strains of the defined mesorhizobium species obtained from the NCBI database.
The genetic differentiation and gene flows among the tested strains from different MLSA clades and the 9 representatives of mesorhizobial type strains sequence (M. ciceri Ca181, M. ciceri bio. Bis. amorphae CCNNGS0123, M. australicum WSM2073, and M. opportunistum WSM2075) were predicted by using the DnaSP 5 software [42]. The parameter used for nucleotide sequence was set as DNA, genomic state was haploid, and chromosomal location was prokaryotic. Further, the concatenated sequence tree of gap, edD, gnD, glnD, rpoB, and nodC was constructed as described previously [43,44]. Briefly, the concatenation of sequences of the same gene of all strains were aligned using MEGA 6.0 [41] and saved as "fasta", then merged using online Fasta alignment joiner tool [45] to form a matrix of concatenated gene sequences.

Statistical Analysis
Pearson correlation test is performed by BioVinci 1.1.5 software, for PGP traits and symbiotic efficiency test. Mean values were compared by Duncan's multiple-range test [46].

Morph-Phenotypic Characters
Phenotypically, all the studied 54 isolates were mucilaginous, round with convex colony, and failed to absorb congo red in the medium. Most of the isolates were moderately fast growing with a generation time (GT) between 4 and 6 h, however, CPN3 and CPN16 were fast growing with GT lesser than 3 h. The bromothymol blue (BTB) test displayed the yellow color (acid production) in most of the isolates, after incubation at 28 °C on YEMA incorporated with BTB while CPN4, CPN13, CPN37, CPN38, and CPN44 produced blue color colonies. In this study, all of the tested symbiont strains showed growth between temperatures 15 to 42 °C. However, only CPN29, CPN32, and CPN34 displayed the average growth at 45 °C.
The P-solubilization index value (PSI) of tested isolates revealed it positive and were found between 2.87 mm to 4.33 mm with respect to reference strain (Table 1) which proved them as good P-solubilizers. In addition, the PSI value of strain CPN8, CPN9, CPN29, CPN32, CPN34, and CPN52 were signified as strong P-solubilizer. The catalase, oxidase, NH3, and IAA tests revealed the potentiality of isolates as PGP enzyme producers (Table 1). IAA production of all the isolates varied from 79.31 to 107.41, 100 μg mL −1 after 72 h incubation (in the presence of tryptophan). All the tested rhizobia were found catalase and NH3 positive. Meanwhile, CPN5, CPN11, CPN13, and CPN14 were found negative to oxidase ( Table 1). The Pearson correlation test of PGP traits (PSI, IAA, ACC-Deaminase, and Biofilm) and SE of tested strains revealed the positive correlation among them. Biofilm efficiency was corelated positively with SE while ACC deaminase is somewhat positive to SE. PSI and IAA have partial or negatively corelated with SE ( Figure 1)
The cross-nodulation results of all 54 rhizobial isolates were strictly considered as monospecific to Cicer arietinum due to failure of nodulation in Phaseolus vulgaris, Pisum sativum, Lens culinaris, and in Vigna mungo.

Carbon Utilization
Supplementary Table S1 shows the loading scores of 95 carbon sources in Biolog GN plates in the first two principal components. The principal component analysis (PCA) showed that first and the second principal component (PC1 and PC2) accounted for 20.2% and 9.6% of data variance respectively, the 2 principal component factors explained 29.9% of total variance. The most important carbon sources in differentiating among our isolates were defined as those that had high variance explained by PC1 or PC2 (>0.70). It was reported that samples located in different spaces were relevant to the ability of carbon substrates utilization. The higher the loading scores; the larger were the effects of carbon source on the principal components. Taking the microbial metabolic pathway of three major nutrients as the basic delineation standards, the carbon sources in Biolog-GN plates were divided into following types: carbohydrates and their derivatives, amino acid and their derivatives, carboxylic acids, amino acids, and miscellaneous. As mentioned in Table S1, 15 carbon sources had an impact on the first principal components (PC1) and 19 carbon sources had impacts on the second principal component. PCA analysis aggregated the isolates primarily into 2 areas that very roughly corresponded to the areas of sampling and alkalotolerance. The variability between the two groups was accounted for by the utilization of different carbon sources ( Figure 2). All the isolates from the Agathuwa region distributed along the PC1 axis, and those from Darshan Nagar and Maznai distributed along the PC2 axis. Most of the isolates, showing very high alkalotolerant properties, correlated with PC1. Group 1 was formed mainly by CPN29, CPN34, CPN35 (from Agathuwa); CPN3, CPN7 (from Mau), and CPN21 (from Raunahi), which shared similar overall substrate utilization patterns, but differed from the isolates CPN6, CPN16, and CPN32 along PC 1. These stations showed a higher utilization of xylitol, L-fucose, D-psicose, hydroxy-L-proline, phenyethylamine, 2-aminoethanol, and L-pyroglutamic than the others. Group 2 was formed by isolates from CPN39, CPN45 (Darshan Nagar), CPN52 (Maznai), CPN1 (Mau), and CPN22 (Raunahi), which had similar carbon utilization patterns, but different from CPN10 (Mau) and CPN16 (Raunahi), along the PC2. Most of them were grouped due to their D, L-lactic acid, D-saccharic acid, turanose, bromosuccinic acid, D-galactonic acid lactone, D-gluconic acid and D, L-lactic acid utilization. Mostly all the highly alkalotolerant isolates distributed along PC1, indicating their similar carbon utilization patterns.

Molecular Identification and Sequence Typing
16S rRNA gene amplification and sequencing of the best SE strains was done. Results designating them as genus Mesorhizobium sp. ciceri with the 98% to 100% identical report, after the annotation and Blast analysis at the NCBI database. Minimum evolution bootstrap phylogeny of the homogenous gene sequences displayed the common clustering among the tested as well as type strains (Figure 3). Moreover, the sequences for four genes coding for gap, edD, gnD, glnD, rpoB, and nodC were simply provided for the further understanding of studied strains relationship ( Figures S2, S3, S4, S4, S6, and S7). Repeated attempts to amplify and sequence the tested gene of screened isolates of chickpea-nodulating rhizobia were successful. The blast report of the above amplified gene fragments found 98% to 100% sequence similarity with genus Mesorhizobium. The gap gene sequences of all the selected strains of M. ciceri revealed great variation in their alignment. The constructed phylogenetic tree was clustered into three groups which differed from 16S rRNA phylogeny. Strains CPN21, CPN29, CPN34, and CPN39 were clustered with Mesorhizobium sp. B7, while the rest of the 12 isolates were clustered with designated Mesorhizobium species (Figure 3). Moreover, other representative rhizobium sp., such as R. undicola MG11875, R. gallicum R.4387, R. giardinni R-4385, R. tropici LGM9503, and R. leguminosarum WSM2304 was separately aligned in the tree. A similar approach was used in the edD gene phylogeny which revealed the comprehensive divergence with the various type strains of Mesorhizobium.
Moreover, the concatenated phylogenetic tree of six different genes of studied isolates and representative type strains reflected the precise diversification ( Figure 4)

Discussion
Bio-inoculation of C. arietinum with rhizobia is now indispensable for agriculture and environment due to its high nutritional value as well as nitrogen fixation and is entirely based on the symbiotic relationship with rhizobia which can be disturbed by several exogenous as well as indigenous factors. Hence, we investigated the actual influencing factors that are required for establishing the symbiotic relationship between C. arietinum and rhizobia. The results obtained from Gram staining, growth on YMA-congo red (CR), and YMA (BTB) as utilization of CR, and from moderately yellow to deep yellow due to their acid producing ability which is a common characteristic of fast growing Mesorhizobium species. Koskey et al. [47] reported similar observations in terms of reaction on the YMA (BTB). Indigenous soybean root-nodulating bacteria were also categorized by the use of YMA-BTB medium in Japan [48]. In light of this characteristic feature, the tested strains were further evaluated for physiological and biochemical factors. The optimum temperature range for ideal growth and activity of rhizobia is between 25 to 31 °C, which is in agreement with previous studies in chickpea rhizobia and other species [49][50]. However, some reports indicate that the maximum growth temperature is 40 °C for both M. ciceri and M. mediterraneum [51][52]. PGP attributes of test strains indicated their competence to increase the soil fertility nutritional value, in which, essentiality of phosphorus (P) as a key element for plant growth and development has been reported widely [52][53]. Phosphorous is one of the most limiting nutrients in tropical soil and only 0.1% of the total P present is available to the plants because of its chemical bonding and low solubility [52]. Interestingly, all the tested strains were able to solubilized the phosphate which proved that rhizobia had the potential to accelerate the solubilization of unavailable soil phosphate and thus could account for high efficiency of phosphorus usage. Similar research was reported previously [54][55]. It could be revealed that Mesorhizobium is useful as a reducing agent of the negative effects of soil calcification on soil P nutrition. The ability to produce IAA by tested strains can have direct effects on production of phytohormones, enzymatic activities, nitrogen accumulation in shoots and seeds, and plant biomass [56][57][58][59]. The isolated strains showed high catalase activity, similar to the observations made by Nakayama et al. [60], which indicated that these isolates could possibly be resistant to environmental, mechanical, and chemical stress. Downie A. J. [61] has observed that high activity oxidase points towards the enhanced ability of strains to effect nodulation. However, CPN11, CPM13, and CPN1 gave negative results for oxidase activity. Each isolate was found positive in NH3 production that might be responsible for the facilitation of limited bacterial growth within the nodules which is controlled by the legumes by diverse ways. Hence, NH3 worked as a major factor for controlled growth of bacteroides inside the nodule [60]. It is established well that plant associated-bacteria lowers the ethylene levels of plants through the ACC deaminase activity and helps to grow under biotic and abiotic stresses [61]. In this study, all the strains were found positive for ACC deaminase, but it was lower than previous reports of ACC-D production by Pseudomonas fluorescens and Streptomyces djakartensis TB4 [61][62][63].
Mono-specificity of tested strains were established by cross-nodulation assay against Cicer arietinum, Pisum sativum, Lens culinaris, and in Vigna mungo. The results revealed that the studied strains nodulated only C. arietinum and no cross-nodulation compatibility was observed between chickpea nodule rhizobia and other mesorhizobia leguminous plants. Thiswas in congruence with previous reports of cross nodulation tests [31,64]. Biofilm formation ability of all 54 tested strains can be considered as the key factor for symbiotic efficiency and it supports the possibility for their commercial use in stress condition [66]. Among all, 16 isolates (CPN1, CPN3, CPN6, CPN7, CPN8, CPN9, CPN10, CPN16, CPN21, CPN22, CPN29, CPN31, CPN32, CPN34, CPN35, CPN39, CPN45, and CPN52) were efficient in biofilm formation, while the rest of the 37 isolates were found moderate or weak in quantitative estimation of biofilm (Table 1). This ability to form biofilms can be positively correlated to symbiotic efficiency test of 16 most promising isolates selected in previous studies [65][66][67]. Wang et al. [68] also indicated that M. huakuii and M. tianshanense are efficient to biofilm formation and were reviewed by Rinaudi L.V. and Giordano W. [10]. Hence, it can be speculated that biofilm is an important factor of an efficient nodulation in chickpea.
Correlation map among the PGP traits (phosphate solubilization, IAA, ACC-deaminase production), biofilm formation, and SE revealed that biofilm is the positive factor efficient nodulation, while ACC-D is partially corelated to SE. Phosphate solubilization and IAA might be useful for the bacterial survival in saprophytic survival.
Based on the above and previous results, CPN1, CPN3, CPN6, CPN7, CPN8, CPN10, CPN16, CPN21, CPN22, CPN29, CPN32, CPN34, CPN35, CPN39, CPN45, and CPN52 were used to study the carbon utilization pattern and functional genes by Biolog and gene sequencing, respectively. Community-level carbon source utilization (Biolog) profiles have recently been introduced as means of classifying microbial communities on the basis of heterotrophic metabolism. Such a classification system might allow microbial ecologists to compare the metabolic roles of microbial communities from different environments without involving tedious isolation and identification of community members. The consumption of carbon substrates present in the Biolog system is a sensitive indicator of short-term changes in the microbial functional diversity [69]. The use of carbon substrates present in the Biolog plate was sensitive enough to detect short-term changes in the microbial functional diversity. Biolog analyses of the sixteen tested isolates revealed differences between soil microbial functional diversity from the region. Prior results showed the rhizosphere to have its own unique soil microbial community capable of producing distinctive metabolic diversity patterns of different composition, and suggested that differences in carbon source utilization can be linked to differences in carbon source availability [70]. There seems to have been no previous studies of metabolic diversity patterns in soil microbial communities in alkaline soils of the IGP region. However, a perusal of our data indicates that the driving environmental factor causing shifts in the metabolic diversity patterns seemed to be the varying amount of alkalinity in the soil.
The accurate identification is mandatory for the strains for further implication and utilization. In this study, the 16S rRNA gene was studied to identify the tested strains. Results revealed that tested strains were identical to Mesorhizobium sp. ciceri and hence, it was proposed that the NW-IGP region of India have positive abundance of C. arietinum nodulating rhizobia. Though, 16S rRNA genes are highly similar within the genus which might be a constrain for their accurate differentiation and identification [71][72]. The usefulness of housekeeping gene sequences for taxonomy and phylogeny analysis has been demonstrated widely [73]. The report displayed the currently published 30 validly species of genus Mesorhizobium in which about 50% were added in the last ten years and has been supported mainly by recA (homologous recombination protein A), atpD (ATP synthase Beta-subunit), and glnII (glutamine synthetase II) phylogenies [19]. Hence, gap, edD, gnD, glnD, rpoB, and nodC gene sequences have been amplified, sequenced, and characterized for the most stable classification. In our study, screened rhizobia were designated as Mesorhizobium on the basis of the 16S rRNA gene, while the species definition and molecular correlation among them was based on the polyphasic functional gene study. This apparent link between gap, edD, gnD, glnD, rpoB, and nodC gene sequences of tested strains and field site may be at least in part due to adaptation of the bacteria to local conditions outside the host plant, and this warrants further study.

Conclusions
In conclusion, the plant growth promotion and symbiotic performances of the isolates revealed the causes and highlighted the factors for efficient nodulation in chickpea. Further, carbon utilization pattern and their biofilm formation ability validated them as best bioprospecting strains in saprophytic as well as endophytic environment. The cross nodulation test has confirmed the strain specificity for chickpea as the restricted host in the NW-IGP region of India and its molecular typing support them as genus Mesorhizobium, which has been revalidated well by the concatenation approached phylogeny.
Supplementary Materials: Figure S1: Geographical location of sampling sites of North west IGP Region of India. Figure S2. Phylogeny showing the relationships between Mesorhizobia and related rhizobia-based aligned gap gene sequences. Scale bars represents 0.1 substitutions per site. Significant bootstrap probability values are indicated at the branching points (only values greater than 80% over 1000 replicates are shown). Figure S3. Phylogenetic tree based on the edD gene sequences, showing the relationships between chickpea mesorhizobia isolates from the NW-IGP region and type strains. Percentage bootstrap is indicated on internal branches (1000 replicates). Scale bars represent the average number of substitutions per site. Figure S4. Phylogenetic trees based on gnD gene sequences, showing the relationships between chickpea mesorhizobia isolates from the NW-IGP region and type strains. Percentage bootstrap is indicated on internal branches (1000 replicates). Scale bars represent the average number of substitutions per site. Figure S5. Phylogenetic tree based on the glnD sequences, showing the relationships between chickpea mesorhizobia isolates from the NW-IGP region and type strains. Percentage bootstrap is indicated on internal branches (1000 replicates). Scale bars represent the average number of substitutions per site. Figure S6. Phylogenetic trees based on rpoB gene sequences. Percentage bootstrap is indicated on internal branches (1000 replicates). Scale bars represent the average number of substitutions per site. Figure S7. Phylogenetic trees based on nodC gene sequences. Percentage bootstrap is indicated on internal branches (1000 replicates). Scale bars represent the average number of substitutions per site. Table S1: Carbon substrates utilized by selected isolates in Biolog plate, significantly correlated to PC1 and PC2 (R > 0.70). Table  S2

Conflicts of Interest:
The authors declare no conflict of interest.