Impact of Plant Growth Promoting Rhizobacteria in the Orchestration of Lycopersicon esculentum Mill. Resistance to Plant Parasitic Nematodes: A Metabolomic Approach to Evaluate Defense Responses Under Field Conditions

The present study deals with biological control of Meloidogyne incognita in 45-days old Lycopersicon esculentum, inoculated with Pseudomonas aeruginosa(M1) and Burkholderia gladioli (M2). The improved plant growth and biomass of nematode infested Plant growth promoting rhizobacteria (PGPR) inoculated plants was observed. Remarkable reduction in the numbers of second stage juvenile (J2s), root galls was recorded after treatment of microbes relative to experimental controls. Moreover, the lowered activities of oxidative stress markers (H2O2 (hydrogen peroxide), O2− (superoxide anion), malondialdehyde (MDA)) was estimated in plants after rhizobacterial supplementation. Higher activities of enzymatic (SOD (Superoxide dismutase), POD (Guaiacol peroxidase), CAT (Catalase), GPOX (Glutathione peroxidase), APOX (Ascorbate peroxidase), GST (Glutathione-S-transferase), GR (Glutathione reductase), DHAR (Dehydroascorbate reductase), PPO (Polyphenol oxidase)) and non-enzymatic (glutathione, ascorbic acid, tocopherol) antioxidants were further determined in nematode infected plants following the addition of bacterial strains. The upregulation of photosynthetic activities were depicted by evaluating plant pigments and gas exchange attributes. An increase in the levels of phenolic compounds (total phenols, flavonoids, anthocyanins), osmoprotectants (total osmolytes, carbohydrates, reducing sugars, trehalose, proline, glycine betaine, free amino acids) and organic acids (fumaric, succinic, citric, malic acid) were reflected in infected plants, showing further enhancement after application of biocontrol agents. The study revealed the understanding of plant metabolism, along with the initiative to commercially exploit the biocontrol agents as an alternative to chemical nematicides in infected fields for sustainable agriculture.

The M. incognita culture was raised and maintained in green house by collection of samples from infected sites. M. incognita was recognized, screened and isolated from the sites to use it further for experimental purpose. They were identified by assessing perineal patterns of females that are formed by expansion or alteration of juvenile bodies, followed by retained tail tips, lateral patterns or presence of phasmids [33]. The roots were washed and cleaned properly (0.9% saline solution) to remove the soil particles. The egg masses were carefully taken off from infected roots and suspended into 45% lactic acid solution. They were allowed to hatch in distilled water into second stage juveniles which were further used to inoculate the tomato and brinjal seedlings raised for maintaining the pure cultures. Second stage juveniles obtained from cultures were then used for experimentation.

Raising of Plants
Authoritative seeds of L. esculentum Mill. (var. Pusa Ruby) were purchased from market, sterilized by 0.01% mercuric chloride and rinsed for 4, 5 times with double distilled water. Soil was prepared by mixing garden soil, sand and organic manure in 3:1:1 ratio followed by autoclaving. Pots were washed and cleaned thoroughly and filled with 300 g soil mixture to which the seeds were sown. After the appearance of true leaves, inoculation of M. incognita was given near roots (300 J2s). Later, the soil was given rhizobacterial treatment after the onset of germination and inoculation of M. incognita. The pots were marked and kept under field conditions and irrigated properly. The harvesting was done carefully after 45-days of inoculation/germination by washing the roots with distilled water to get rid of soil particles. The plants were stored to carry out morphological, physiological and biochemical investigations.

Growth Attributes
The growth attributes of 45-days old L. esculentum plants were measured through assessing root length and shoot length. Further, plant biomass was also determined in terms of fresh weight and dry weight of samples.

Total Chlorophyll and Carotenoid Content
The estimation of total chlorophyll and carotenoid content was done by Arnon [34] and Maclachlan and Zalik [35] method. For this, 1 g sample was grounded in 80% acetone, followed by centrifugation in cooling centrifuge with 4 • C at 10,000 rpm for 15 min. After that, absorbance of supernatant was read at 645 and 663 nm respectively for total chl and 480 and 510 nm for carotenoid content.

Total Xanthophyll Content
Estimation of total xanthophyll content was done by Lawrence [36] protocol. For this purpose, 50 mg plant sample completely dried and crushed was taken in the flask. The extractant (30 mL) was prepared by combining hexane (10 mL): absolute alcohol (6 mL): toluene (7 mL): acetone (7 mL) and transferred and meticulously mixed for 15 min. To this, 2 mL 40% Met. KOH was mixed and heated to 58 • C in water bath. Sample incubation was done under dark conditions for an hour after which 30 mL of hexane in addition to 10% Na 2 SO 4 was combined in flask to make it to 100 mL volume. Again, incubations for an hour were done under dark conditions after vigorously mixing the contents. Later, the upper layer was removed to other flask and hexane was added to make final contents of 50 mL after which the optical density at 474 nm was read.

Gaseous Exchange Parameters
The measurement of gas exchange attributes were done in 45-days old plants. Estimation of Net photosynthetic rate (Pn), Transpiration rate (E), Stomatal conductance (Gs) and intercellular CO 2 rate (Ci) was done by Portable Photosynthesis Measuring System Unit (Li COR-6400, LiCOR Instruments, Lincoln, NE, USA) during sunny day. Maintenance of relative humidity (80-90%); temperature (25 • C); photon flux density (1000 µmol m −2 g −1 ) and CO 2 rate (400 µmol mol −1 ) were the essential conditions for the standard working of the instrument.

Superoxide Anion Content
The measurement of superoxide anion content was done after following Wu et al. [37] method. Crushing of 1 g plant material in 2% polyvinylpyrrolidone (PVP) and 0.5% Triton X-100 containing, 4 mL 50 mM phosphate buffer (pH-7.8) was done. Then, they were spun at 10,000 rpm for 15 min.
After removing the 0.5 mL supernatant in separate tube, 0.5 mL of phosphate buffer and 0.1 mL hydroxylamine hydrochloride was mixed. Incubations at 25 • C for 30 min were done. To this, 1 mL each of 1-naphthylamine and 3-aminobenzenesulphonic acid was supplemented and placed for 15 min at 26 • C. Optical density at 530 nm was further read.

Hydrogen Peroxide Content
The estimation of hydrogen peroxide content by Velikova et al. [38] proposed method was done. Homogenisation of 100 mg of plant tissue in freshly prepared 0.1% Trichloroacetic acid (TCA) was done in cold motor pestle. They were then subjected to centrifuge at temperature of 4 • C, at 10,000 rpm for 20-25 min. The supernatant was removed in the separate tube and mixed with 0.5 mL of potassium phosphate buffer and 0.5 mL KI. Optical density at 390 nm was then recorded.

Malondialdehyde Content
Estimation of malondialdehyde (MDA) accumulation was conducted by Heath and Packer [39] method. To start with the measurement, 1g of fresh plant material was macerated using freshly prepared 0.1% 5 mL TCA. Samples were centrifuged at 8000 rpm at 4 • C for 15 min. 1 mL supernatant was removed in a separate tube to which 4 mL of 0.5% thio-butyric acid (TBA) (made with 20% TCA) was combined. The incubations for 30-35 min at 90 • C were given by placing the tubes in water bath. After cooling the tubes, the optical density at 532 nm was read.

Non-Enzymatic Antioxidants
To estimate the levels of non-enzymatic antioxidants, fresh plant tissue was grounded in 50 mM of 4 mL Tris Buffer with pH 10 using cold, sterile and autoclaved motor pestle. After that, centrifugation in cold centrifuge (4 • C) at 10,000 rpm was done for 15-20 min. Collection of supernatant in other tubes was done to quantify antioxidants (glutathione, ascorbic acid and tocopherol).

Glutathione Content
Estimation of glutathione content was done through standard protocol offered by Noctor and Foyer [49]. For this, combination of following contents, l00 µL supernatant, 50 µL DTNB, 4 mL absolute alcohol and l000 µL Tris buffer were taken in a tube and placed under room temperature for 10 min. The tubes were vortexed for 10 min until a yellow coloured product is formed. Optical density at 412 nm was read.

Ascorbic Acid Content
Estimation of ascorbic acid was conducted through procedure opted by Roe and Kuether [50]. For this, combination of 1 mL supernatant, 1 mL freshly prepared 50% TCA and 0.5 mL DPNH was done. The tubes were shaken vigorously after which they were incubated at room temperature for 4 h. The osazones crystals produced were dissolved in cold sulphuric acid. Incubation in ice for 20 min was given and absorbance at 540 nm was read.

Tocopherol Content
Estimation of tocopherol content was done by Martinek [51] established protocol. For this, 500 µL homogenate was combined to 500 µL distilled water, 500 µL absolute ethanol prepared by adding 0.12% FeCl 3 6H 2 O. The tubes were shaken vigorously until the protein crystals were developed. Crystals were dissolved in 500 µL xylene and vortexed for 2 min. The layer formed on upper phase was combined to TPTZ in separate tube and optical density at 600 nm was read.
2.9. Phenolic Compounds 2.9.1. Total Phenols Estimation of total phenol content was conducted through method examined by Singleton and Rossi [52]. For this, 500 mg oven dried sample was grounded in 60% ethanol, after which it was warmed at 65 • C using water bath for 20 min. Samples were filtered by Whatmann (GE Healthcare Life Sciences, Maidstone, UK) filter paper till the residue was left behind. 100 mL of total sample volume was maintained by 60% ethanol. There was 2 mL of this sample mixed in another tube containing 8 mL FC reagent and 8 mL sodium carbonate. Incubation for 4 h was given followed by measuring optical density at 765 nm. (Standard curve prepared by gallic acid).

Total Flavonoids
Estimation of total flavonoids was conducted by standard procedure given by [53]. For this, 100 mg of dried sample was macerated in 4 mL absolute alcohol. The preparation of extractant was done by combining 3 mL each of 5% NaNO 2 , 5% AlCl 3 and distilled water. 1 mL of above prepared extractant was then combined to 2 mL sodium hydroxide and 2.5 mL of distilled water in sample tube. Incubation was done for 10 min followed by measuring optical density at 510 nm (Standard used: Rutin).

Total Anthocyanins
Estimation of total anthocyanins was done through procedure followed by Mancinelli [54]. For this, 1000 mg of plant material was crushed in reagent prepared by combining absolute methanol, distilled water and hydrochloric acid in the ratio of 78:20:2. The sample was incubated under cold conditions (4 • C) for 12 h and spun at 14,000 rpm for 20 min. Absorbance at 530 and 657 nm was read to calculate the anthocyanin content.

Total Osmolytes
Estimation of total osmolytes was done with Vapour Pressure Osmometer (Vapro 5600, ELITechGroup, Paris, France). The fresh plant samples were stored in liquid nitrogen and thawed for sap extraction using syringe. The sap was collected for further analysis. Only 10 µL of sap was used for analysis and loaded directly on filter paper discs kept upon the instrument. Calibration of instrument was done using Sodium chloride standard kits with 1000, 290 and 100 mOsm osmolarities that come along with the instrument. Values displayed on the instrument were noted after maintaining room temperature conditions of 22, 23 • C.

Total Carbohydrates
Total carbohydrates were calculated through the Hedge et al. [55] procedure. For this, 100 mg sample was warmed in 4 mL of 2.6 N HCl for 5 h. Incubations at room temperature were given followed by the addition of sodium carbonate to offset the reaction. A total volume of 30 mL was maintained by addition of 25 mL distilled water. There was 4 mL of anthrone reagent taken in the separate tube to which 1 mL of above sample was added and subjected to boil for 15 min. Samples were cooled and optical density at 630 nm was read (d-glucose: standard).

Total Reducing Sugars
The analysis of total reducing sugars was done by method provided by Miller [56]. For this, 100 mg dried and finely powdered sample was crushed in 80% absolute alcohol. Preparation of DNSA was done through combination of 0.2 g phenol crystals and 0.05 g Na 2 SO 3 under cold conditions. In total, 3 mL of this freshly prepared reagent was added to 1 mL of above supernatant and 40% KNaC 4 H 4 O 6 ·4H 2 O. Incubation for 2 min was given at room temperature to record the optical density at 510 nm.

Trehalose Content
Estimation of trehalose content was done through Trevelyan and Harrison [57] assay. The analysis was begun by crushing 500 mg oven dried sample in 80% ethanol and further centrifuged under cold conditions of 4 • C at 10,000 rpm for 20 min. The components of reaction were combined by adding 100 µL of above extractant, 4 mL anthrone reagent and 3 mL freshly prepared TCA. Optical density at 620 nm was read after the formation of deep yellow coloured product (d-glucose: standard).

Glycine Betaine Content
The measurement of glycine betaine (GB) content was predicted by Grieve and Grattan [58] protocol. To begin with the analysis, 500 mg dried sample was crushed in 0.05% toluene and 4 mL double distilled water which was then incubated for 24 h. Samples were filtered through filters with 0.2 µm pore size. The reaction components were mixed by combining 500 µL supernatant, 1000 µL 3 N HCl and 100 µL potassium iodide. Incubations of samples were done in ice for 4 h. To the reaction tubes, 3 mL chilled water and 8 mL 1,2-dichloroethane was mixed. The tubes were vortexed till the formation of two different layers. The above layer was thrown, and bottom pink layer was kept for analysis. Optical density at 365 nm was then read (Standard: betaine hydrochloride).

Proline Content
Estimation of proline content was done through Bates et al. [59] analysis method. For this, 500 mg sample was crushed in 10 mL 3% sulphosalicylic acid. Centrifugation of sample tubes at 10,000 rpm for 15-20 min was done. Sample tubes containing 2 mL supernatant was separated and combined to 3 mL ninhydrin reagent and 3 mL glacial acetic acid. The tubes were heated to 100 • C and placed in ice after 5 min to counterbalance the reaction. Later, 5 mL toluene was further mixed to it and tubes were vortexed for 3 min. The red toluene layer was formed of which the optical density at 520 nm was read (Standard: l-proline).

Free Amino Acid Content
Estimation of free amino acid content was done via Lee and Takahashi [60] method. To begin with, 100 mg dried sample was crushed using 4 mL absolute ethanol and warmed at 50 • C for 10 min using water bath. The reaction tubes were vortexed for 5 min. Later, freshly formed ninhydrin reagent was combined with 200 µL sample and heated for 5 min. Optical density at 570 nm was taken after the formation of blue colour.

Organic Acid Profiling
Organic acid profiling was done through following Chen et al. [61] method with slight modifications. To begin with the analysis, 50 mg dried crushed sample was taken to which 500 µL absolute methanol and 500 µL 2.5 N HCl was pooled. Sample tubes were then placed in shaker for 5 h, and centrifuged under cold conditions (4 • C) at 10,000 rpm for 10 min. The residue was left behind and above liquid was transferred into another tube to which 0.3 mL methanol and 0.1 mL 50% conc. sulphuric acid was combined, vortexed and incubated for 10 h at 65 • C in water bath. Sample tubes were cooled at 28 • C and 800 µL chloroform and 400 µL double distilled water mixed into it. Vortexing was done till the formation of two different layers. The basal layer was collected for organic acid analysis by GC-MS. The system (Shimadzu GC-MS-Q2010Plus, Japan) required 2 µL sample for analysis. Standard conditions of the instrument were: (i) for GC: helium carrier gas, column with flow rate of 1.7 mL/min, analytical column with length of 30 m, ID-0.025 mm and DB, 5 ms, temperature of 125 • C and 25 • C/min; Biomolecules 2019, 9, 676 9 of 30 (ii) for MS: ion source with temp 200 • C and detection mode is relative. The organic acid quantification was carried out by comparing mass spectra or peaks with National Institute of Standard and Technology and Wiley 7 Library.

Statistical Analysis
The morphological, physiological and biochemical data was statistically tested through the self-built program in Microsoft Excel (Microsoft Office Excel 2007, Albuquerque, New Mexico, USA). Results in the form of mean ± standard deviation (SD) (level of significance checked at p ≤ 0.05 and 0.01) were presented. Null hypothesis (H 0 ) was verified by two-way ANOVA (analysis of variance). Data was also confirmed with Tukey's multiple comparison test and HSD (honestly significant difference).

Effect of Biocontrol Agents on Growth Parameters of PPN Infected L. esculentum Plants
The ramifications of microbial strains (P. aeruginosa (M1) and B. gladioli (M2) on morphological parameters of 45-days old L. esculentum plants raised under PPN infestation was investigated in terms of root length, shoot length, biomass (fresh weight, dry weight) and number of galls. The results indicated a remarkable reduction in the root length and shoot length of plants by 43.2 and 21.5% relative to control plants. The decline in the plant biomass of infected plants (fresh weight (35.1%) and dry weight (47.6%)) was also reported. However, after the inoculation of P. aeruginosa (M1) the enhancement in the root length, shoot length, fresh weight and dry weight was noted by 73.3, 54.3, 33.2, and 105.4%, respectively. The supplementation of B. gladioli (M2) in infected plants also maximised the root length (51.6%), shoot length (69.7%), fresh weight (49.6%) and dry weight (110.5%) when compared to PPN experimental plants. Furthermore, the number of galls were also measured in PPN inoculated plants. The gall numbers (27) were counted in plants inhabited with nematodes, and their significant reduction was noted in the plants supplemented with rhizobacterial strains. Reduction of 28.06% was observed after P. aeruginosa (M1) addition while reduction to 19.5% was seen after B. gladioli (B2) treatment in contrast to experimental controls (Table 1).

Effect of Biocontrol Agents on Photosynthetic Pigments of PPNPPN Infected L. esculentumPlants
The plant pigments were determined in terms of total chlorophyll, carotenoid and xanthophyll. It was evaluated that level of total chlorophyll was reduced to 19.7% in plants raised with nematode infestation as compared to experimental controls. But, supplementation of P. aeruginosa (M1) to nematode infected plants, led to uplifted levels of total chlorophyll from 19.7 to 68.5% respectively. Also, the inoculation of B. gladioli (M2) in infected plants stimulated the total chlorophyll content in plants by 91.1% respectively. The levels of carotenoids estimated, showed lowered levels by 60.9% in plants infected with nematodes relative to controls. After the application of P. aeruginosa (M1), a sharp enhancement in the level of carotenoids by 183% was noticed in nematode infested plants and by 220% after the inoculation of B. gladioli (M2). In addition, the levels of xanthophyll were also minimised by 33% in plants raised under nematode treated soils in contrast to control. A remarkable stimulation in xanthophyll content was measured by 56.4 and 86.7%, after the inoculation of P. aeruginosa (M1) and B. gladioli (M2) respectively, in nematode infested plants when compared with experimental controls (Figure 1).

Effect of Biocontrol Agents on Photosynthetic Pigments of PPNPPN Infected L. esculentumPlants
The plant pigments were determined in terms of total chlorophyll, carotenoid and xanthophyll. It was evaluated that level of total chlorophyll was reduced to 19.7% in plants raised with nematode infestation as compared to experimental controls. But, supplementation of P. aeruginosa (M1) to nematode infected plants, led to uplifted levels of total chlorophyll from 19.7 to 68.5% respectively. Also, the inoculation of B. gladioli (M2) in infected plants stimulated the total chlorophyll content in plants by 91.1% respectively. The levels of carotenoids estimated, showed lowered levels by 60.9% in plants infected with nematodes relative to controls. After the application of P. aeruginosa (M1), a sharp enhancement in the level of carotenoids by 183% was noticed in nematode infested plants and by 220% after the inoculation of B. gladioli (M2). In addition, the levels of xanthophyll were also minimised by 33% in plants raised under nematode treated soils in contrast to control. A remarkable stimulation in xanthophyll content was measured by 56.4 and 86.7%, after the inoculation of P. aeruginosa (M1) and B. gladioli (M2) respectively, in nematode infested plants when compared with experimental controls (Figure 1

Effect of Biocontrol Agents on Gas Exchange Parameters of PPN Infected L. esculentum Plants
The results depicted an alteration in the gas exchange parameters of plants raised under the influence of nematode infection and bacterial strains. The results displayed a noticeable reduction in net photosynthetic rate, stomatal conductance, intercellular CO2 and transpiration rate by 33.9, 31.6, 11.3 and 40.3% respectively, in 45-days old PPN infected plants. Supplementation of P. aeruginosa (M1) in infected plants led to an upregulation in all these gas exchange parameters. The inclination in the net photosynthetic rate, stomatal conductance, intercellular CO2 and transpiration rate was revealed by 32.7, 80.2, 8.64 and 30% respectively, after inoculations in the infected plants. The inoculation with B. gladioli (M2) also enhanced the activities of net photosynthetic rate, stomatal conductance, intercellular CO2 and transpiration rate by 87, 88.5, 24.4 and 46% respectively, in infected plants in comparison to experimental controls ( Table 2).

Effect of Biocontrol Agents on Gas Exchange Parameters of PPN Infected L. esculentum Plants
The results depicted an alteration in the gas exchange parameters of plants raised under the influence of nematode infection and bacterial strains. The results displayed a noticeable reduction in net photosynthetic rate, stomatal conductance, intercellular CO 2 and transpiration rate by 33.9, 31.6, 11.3 and 40.3% respectively, in 45-days old PPN infected plants. Supplementation of P. aeruginosa (M1) in infected plants led to an upregulation in all these gas exchange parameters. The inclination in the net photosynthetic rate, stomatal conductance, intercellular CO 2 and transpiration rate was revealed by 32.7, 80.2, 8.64 and 30% respectively, after inoculations in the infected plants. The inoculation with B. gladioli (M2) also enhanced the activities of net photosynthetic rate, stomatal conductance, intercellular CO 2 and transpiration rate by 87, 88.5, 24.4 and 46% respectively, in infected plants in comparison to experimental controls (Table 2).  that mean values of treatments are significantly different at p < 0.5 according to Tukey's multiple comparison.

Effect of Biocontrol Agents on Oxidative Damage in PPN Infected L. esculentum Plants
The oxidative damage was assessed by estimating the levels of oxidative stress markers namely, superoxide anion, H2O2 and MDA. It was revealed in the present investigation that PPN infestation in 45-

Effect of Biocontrol Agents on Antioxidative Enzymes Activities of PPN Infected L. esculentum Plants
The  (Table 3).  (Table 3).

Effect of Biocontrol Agents on Non-Enzymatic Antioxidants of PPNInfected L. esculentum Plants
The results pertaining to non-enzymatic antioxidants, glutathione, ascorbic acid, and tocopherol were found to be modulated in plants subjected to nematode infection and microbial treatment. A remarkable stimulation in the non-enzymatic antioxidants was estimated in the nematode infected plants. The glutathione content was raised by 51.2% in nematode infested plants relative to experimental controls. This was further accelerated after the addition of P. aeruginosa (M1) in infected plants and by 7.11 and 16.9% after the treatment of B. gladioli (M2) respectively. A similar observation was recorded for ascorbic acid levels where an escalation in the ascorbic acid levels was noted by 36.9% in infested plants while that of controls. This was further augmented by P. aeruginosa (M1) and B. gladioli (M2) independently, after which an enhancement in the ascorbic acid levels by 10.9% and 17.2% were reported. Another antioxidant, tocopherol was also estimated and found to be stimulated by 37.3% in PPN infested plants. Moreover, an inoculation of P. aeruginosa (M1) raised the tocopherol content in PPN affected plants by 9    observation was recorded for ascorbic acid levels where an escalation in the ascorbic acid levels was noted by 36.9% in infested plants while that of controls. This was further augmented by P. aeruginosa(M1) and B. gladioli (M2) independently, after which an enhancement in the ascorbic acid levels by 10.9% and 17.2% were reported. Another antioxidant, tocopherol was also estimated and found to be stimulated by 37.3% in PPN infested plants. Moreover, an inoculation of P. aeruginosa(M1) raised the tocopherol content in PPN affected plants by 9

L. esculentum Plants
The present study depicted that measurement of total phenols, flavonoid and anthocyanin levels were enhanced in the plants stressed with nematodes. The upliftment in the levels of total phenols occurred by 25.56%, flavonoids by 20.93% and anthocyanins by 58.46%, respectively. However, the further enhancement in the levels of these phenolic compounds were observed after the treatment of P. aeruginosa (M1) by 40.37% (total phenols), 40 23.09% respectively, in PPN stressed plants in comparison to experimental control plants. Increase in the total osmolytes was recorded by 25.97 and 33.26% in stressed plants after treating with P. aeruginosa (M1) and B. gladioli (M2) respectively. Besides, the total carbohydrates were also raised by 23.26 and 35.21% in infected plants in the presence of P. aeruginosa (M1) and B. gladioli (M2) respectively. The reducing sugar content was also improved by 27.7% after the amendment with P. aeruginosa (M1) and 40.87% with B. gladioli (M2) treatment ( Figure 5).

Effect of biocontrol Agents on Organic Acid Profiling in PPN Infected L. esculentum Plants
Organic acid profiling was done to assess the levels of fumaric acid, succinic acid, citric acid and malic acid in plants exposed to nematode infestation. An upregulation in the levels of organic acids, fumaric acid (9.7%), succinic acid (29.67%), citric acid (42.26%) and malic acid (57.52%) was noticed in PPN exposed plants in contrast to experimental controls. The drastic elevation in fumaric, succinic, citric and malic acid was recorded by 27.07, 12.02, 43.9 and 118.8% respectively, after the PPN infected plants were amended with P. aeruginosa (M1). Similarly, organic acid levels of infected plants inoculated with B. gladioli (M2) were also enhanced by 38.05% (fumaric acid), by 7.16% (succinic acid), by 62.66% (citric acid) and by 140.56% (malic acid) respectively (Table 4).

Discussion
Microbial community present in the rhizosphere have been explored to be useful as natural biocontrol agents in managing PPN infection in L. esculentum and improving their growth conditions under normal as well as infectious stages [32]. Our studies indicated that PPN infection declined the growth of 45-days old infected plants in terms of root length, shoot length, plant biomass (fresh and dry weight) and induced gall formation in plants. Earlier also, reduction in the growth attributes have been reported by Gupta et al. [62], Bacopa monnieri plants raised under the influence of nematodes. Moreover, reduction in growth of Salvado rapersica L. and S. oleoides Decne was also reported in terms of decreased plant height, root length, root and shoot biomass, due to PPN infection in the plant species [28]. The reduction in the growth parameters of plants occurs due to infestation of PPN into the roots followed by their further proliferation into interior region leading to development of galls over the root surfaces [63]. Furthermore, the development of galls in root zone of infected plants have also been recorded in the present study. Similar to our observation, the formation of galls have also been noticed in tomatoes and carrots due to the infection of M. incognita [19]. They illustrated that nematodes are quite effective in penetrating the root system both fibrous and tap root through secretion of array of compounds to dissolute cell wall for invasion [19]. The present study also depicted that inoculation of P. aeruginosa(M1) and B. gladioli (M2) improved the morphological characters of the plants inhabited with PPNs. Our studies are in the agreement to previous studies conducted by Abd-El-Khair et al. [64], who reported that biocontrol agents Bacillus subtilis, B. pumilus and P. fluorescens had a protective role against PPN M. incognita in improving root length, shoot length, fresh and dry weights in cowpea. According to their assumptions, the biocontrol agents have the stimulating effects on plants by inducing growth, plant nutrition, producing plant hormones and suppressing various pathogens in the surrounding zones [64]. Moreover, they can also survive under severe conditions to act aggressively against PPNs by producing growth promoting materials in order to cease the negative effects of pathogens [65]. The reduction in the gall number, egg masses and gall index have also been reported by Viljoen et al. [19] in tomato and carrot plants infected with PPNs on supplementation of B. firmus, B. aryabhattai, B. cereus, Paenibacillus baricinonensis and P. alvei respectively. They concluded that the reduction in the gall formation in the plant species is mainly because of antagonistic effects of these microbial strains on PPNs each possessing its own mode of action to inhibit their proliferation around the root zones [19]. It was also reported that Bacillus and Paenibacillus sp. are involved in the production of different nematicidal substances that directly inhibit the growth activities of PPNs [66]. Similarly, plant growth promotion has also been studied in nematode infected cucumber plants after the supplementation of Photorhabdus luminescens in terms of leaf number, root and shoot length, plant biomass, and fruit yield [67]. The improvement in the growth parameters and reduction in the gall formation of plants in current study might be due to production of extracellular enzymes (serine alkaline protease, chitinases, glucanases etc.), antibiotics and toxins to induce mortality of nematodes [19,32,67]. The main mechanism associated with escalated resistance is the colonisation of the plant species to mediate systemic resistance against parasitic nematodes. This tends to barricade the biting of nematodes to plant roots thereby, inhibiting their penetration inside roots to the maximum extent.
The results of the present study also demonstrated the reduction in the plant pigments in plants infected with nematodes. Our results concur with many similar findings, as it is extensively reported that Meloidogyne sp. reduced chlorophyll pigment in tomato plants [68]. They found that disruption in the plant pigments occur mainly due to hindrance in the reactions associated with photosynthesis due to nematode infection, that caused the loss of plant pigments derived during photosynthetic processes [68]. Moreover, nematode damage the root and leaf tissues, chloroplasts and disturb the translocation of photosynthates in plants [69]. The studies conducted in Capsicum annuum L. showed a significant reduction in the total carotenoid levels after the inoculation of M. incognita [70]. Our results also demonstrated an elevation in the levels of plant pigments in nematode infected L. esculentum plants after the rhizobacterial treatment. Our findings coincide with the studies reported by Abd-El-Khair et al. [64], who found the raised levels of chla, chlb, total chlorophyll and carotenoid content in nematode infested cowpea plants. Furthermore, a study reported in nematode infected Pinus pinaster revealed a decreased incidence of nematode proliferation after microbial treatment along with the prevention in the degradation of plant photosynthetic pigments. They found a significant stimulation in the pigment content in infected plants due to protection of key enzymes involved in photosynthesis from nematode attack [71]. The improved plant pigment levels in the present study are most likely due to protection of photosynthetic apparatus by microbes along with their role in improving translocation of nutrients in plants that ultimately regulates the effective photosynthesis. Moreover, the mechanism entails microbe mediated limitation of PPN migration inside the vascular tissues that ramifies in effective transport of essential minerals required for photosynthesis.
A down regulation in the activities of gas exchange parameters have also been noted in nematode infested plants in the present work. Our studies corroborate with the findings of Goulart, et al. [70], who suggested that nematode infected coffee plants showed remarkable reduction in the transpiration rate, stomatal conductance and CO 2 concentration. Moreover, they found that overall rate of photosynthesis was lowered in infected plants with reduced starch synthesis [72]. Similarly, the reduction in gaseous exchange parameters, mainly transpiration and assimilation rate in nematode infected Solanum tuberosum plants was also recorded [73]. The decrease in the gaseous exchange mechanisms in plants under nematode infection are mainly due to the impairment in the water transport pathways because of root damage by nematodes [74]. Moreover, the reduction in the net photosynthesis accompanied by stomatal regulation and photochemical efficiency was observed in C. annuum under the infection of root-feeding parasitic nematodes [70]. The increase in the gaseous exchange parameters was noticed in the present study when an inoculation of P. aeruginosa and B. gladioli was done. The results of present study coincide with the earlier studies conducted by Mioranza et al. [75] in tomato plants infested by M. incognita that showed an elevation in photosynthetic rate, stomatal conductance, internal CO 2 concentration, leaf temperature and transpiration upon supplementation of Thuya occidentalis to infected plants. Furthermore, it was revealed that M. incognita infected tomato plants led to enhanced effects upon stomatal conductance and transpiration rate in the presence of B. lichenoformis [76]. The acceleration of gaseous exchange parameters in the current study is most likely due to higher synthesis of enzymes involved in photosynthetic processes.
The current study also reported that nematode infection uplifted the levels of oxidative stress markers in plants. Our studies find support from our earlier studies, Khanna et al. [32], where levels of superoxide anions, H 2 O 2 and MDA were found raised in 10-days old seedlings of tomato under the influence of PPNs. This is directly linked to oxidative damage in plants due to pathogenicity of parasitic nematodes in the form of ROS (reactive oxygen species) generation. Moreover, studies have been reported where levels of total reactive oxygen species, superoxide anions, H 2 O 2 and PCD (programmed cell death) were elevated in B. monnieri L. infected with PPNs [77]. The accumulation of oxidative stress markers and free radicals occurs as a result of toxic by-products of metabolic pathways in plants, when their levels exceed the optimum limits on exposure to parasitic nematodes [78]. Our studies also showed a significant reduction in the levels of all the stress markers in the presence of biocontrol agents. Similar to our studies, the pathogen induced overproduction of ROS has been reported to damage lipid synthesis and biomolecules, that has been considerably alleviated in the presence of microbial agents, implying their role as true antagonists against PPNs [79]. However, the reduction in the levels of MDA have also been observed in Cicerarietinum plants, prominently after their treatment with microbial consortium [80]. Another study revealed that nematode infested B. monnieri L. plants when inoculated with biocontrol agents (Chitinphilus sp. and Streptomyces sp.) lowered lipid peroxidation, superoxide anions and hydrogen peroxide levels [77]. They speculated that these inoculating agents had a stimulatory effect in defense mechanisms through reducing the free radicals generated and strengthening their immune system. In addition, they also strengthen the cell membrane of plant tissues that prevent their negative effects on biomolecules and also in preventing their movement towards other parts of the plant [81]. The mechanism associated with the reduced levels of oxidative stress markers in the present study by the aid of microbes is most likely due to elevation in the chain of resistance processes and antioxidative defense genes in plants along with their abilities to directly scavenge the ROS initiated during stress. Along with this, they also restricts their invasion inside the plants so as to balance the redox homeostasis with the plant inner zone.
The up regulation in the activities of enzymatic as well as non-enzymatic antioxidants have also been observed in the present study in plants associated with nematode infection. A study reported in barley infected with PPNs showed enhanced activities of SOD, POD, CAT and PPO [82]. They reported that increase in antioxidative enzymes are relative to resistance mechanisms developed by plants. Moreover, enhanced activities of peroxidases (APOX, GPOX and POD) are directly linked to their involvement in providing structural rigidity to plant tissues through enhanced lignin synthesis in cell walls in order to prevent the nematode penetration in plants [83]. Increases in the activities of peroxidases have also been reported in Ipomoea batatas subjected to PPN infection as assessed through transcriptomic analysis [84]. Moreover, an elevation in GPOX, CAT and SOD activities was observed in tomato plants infected by M. incognita [85]. The increase in the activities of GR, GST and other antioxidants in plants exposed to nematode infection are due to higher glutathione pools that further activates the enzymes DHAR and GPOX, key enzymes of GR cycle. Also, the level of non-enzymatic antioxidants glutathione, ascorbic acid and tocopherol were enhanced in nematode affected plants in the present study. AsA-GSH pathway is the most important pathway that gets activated during nematode infection in plants. It has been reported to activate the antioxidants as a defense strategy in plants followed by maintaining the redox state during plant-nematode interactions [86]. These antioxidants have been known to directly participate in free radical scavenging, against ROS accumulated during nematode attack [86]. Additionally, tocopherol have been seen to scavenge lipid peroxy radicals accumulated within cell membrane via free radical reactions [87]. Accumulation of different antioxidants have been reported in grape vine by Kesba and El-Beltagi [88] during nematode infection. Our results also showed further upliftment of enzymatic and non-enzymatic antioxidants following treatment with rhizobacterial strains. Our findings corroborate with the studies conducted by Molinari and Leonetti [89], who reported modulation in the activities of defense genes encoding antioxidative system CAT, APOX in nematode infected tomato plants. They also reported that SAR (Systemic Acquired Resistance) is the proposed mechanism that activates molecular signaling at different stages of nematode infection within plants, thereby preventing their penetration and movement inside roots [89]. Furthermore, another study reported in tomato plants infected with PPNs showed higher activities of SOD and POD after the inoculation of Drechslerella dactyloides and Dactylaria brochopaga [90]. According to their research, the microbes participate in the defense signaling during the presence of pathogens, that tend to up regulate the activities of specific biomolecules and antioxidants [90]. Moreover, SOD being a first line of defense in dismutating H 2 O 2 into non-toxic forms is foremost activated by microbes against ROS generated in cytosol during respiration and photosynthesis [91]. Elevation in the SOD activity and ascorbic acid content have also been reported in nematode infested B. monnieri L. plants under the influence of Chitophilus sp. and Streptomyces sp. [62,77]. Also, the synthesis of non-enzymatic antioxidants were enhanced in these plants. The mechanism associated with stimulatory effects of microbes on antioxidative defense system is most probably due to up-regulation of genes and protein transcripts of antioxidant molecules. And this further attributes to the adaptive responses of plants with protection from PPN infection. The microbes have the capacity to incline the immune responses in plants against toxic environment created within them by advent of PPN proliferation.
Present results also demonstrated the accumulation of phenolic compounds in nematode infested plants. The increased levels of total phenolics have also been reported by Gálvez et al. [70] in C. annuum after PPN inoculation in order to induce the resistance in plants. Furthermore, a study conducted in Oryza sativa associated with PPN showed higher levels of total phenols (TP) in plants. They depicted that TP induced lignification in epidermal regions of plants, as a constitutive and inducible post penetration mechanism for providing resistance against nematode infection [92]. Flavonoids also act as signaling molecules during plant defense towards PPNs in rhizosphere in modulating auxin transport during gall formation. A study depicting role in flavonoids in nematode parasitized Medicago truncatula was conducted, that revealed an induction in the flavonoid pathways to accumulate different flavonoids such as chalcones and flavonones (kaempferol, flavone, isoflavonoids, medicarpin, pterocarpans etc.) to protect them against PPNs and inhibiting gall development in plants [93]. An increase in the levels of anthocyanins was also observed by Labudda et al. [94], in Arabidopsis thaliana infected with nematodes as an important antioxidant in protecting the plants from infection. Our study showed a further enhancement in the levels of phenolic compounds when rhizobacterial inoculations were done in PPN infected plants. It was revealed by Nunes da Silva et al. [71] in nematode stressed P. pinaster that biosynthesis of phenolics was dramatically increased with the reduction of nematode population at that specific site, being a crucial metabolite in controlling the pathogen prevalence. This is mainly attributed to activation of shikimate biosynthetic pathway by microbes to overproduce metabolites involved in the resistance mechanisms [71]. It was further suggested that inoculation of D. dactyloides and D. brochopaga in nematode infected tomato plants led to accumulation of phenolic compounds and enzymes involved in its synthesis [90]. They reported that PGPR activated PAL (Phenylalanine ammonia lyase) pathway in roots to improve plant growth, trap nematode population and initiate ISR in plants. Increase in the total phenolic content was also recorded in nematode infected cowpea plants under the influence of different rhizobacterial strains such as B. subtilis, P. fluorescens and B. pumilus respectively [64]. This phenomenon is quintessential to eradicate the impact of infection from the plants. It can be relatable to direct involvement of phenolic compounds in preventing the impregnation of nematodes inside plant tissues. The phenol deposition at cell walls makes intricate for PPN to biting and invading further inside plants.
Secondary metabolites also known as defense-related metabolites are activated in plants against pathogen attack. They play an important role in defense by triggering direct toxic effects on pathogen or by acting as signaling molecule to activate defense mechanisms in plants [95]. Present study showed a stimulation in the secondary metabolites (osmoprotectants and organic acids) in nematode infested plants. A remarkable increase in the total carbohydrates was observed in A. thaliana after the inoculation of PPN [94]. They speculated that an enhancement in the carbohydrates is due to alteration in the photosynthetic efficacy, in order to maintain the permanent inflow of sugars for nematodes prevailing in roots. The amino acid levels were also reported to be lessened in C. annuum plants infected with PPNs [70]. They formulated that Leu, Gly, Phe, Cys, Gly, Thr, Val and Ala were found to be declined in infected plants due to the alterations in the expression levels of different genes involved in amino acid metabolism [70]. Moreover, a study conducted in A. thaliana roots also showed a higher accumulation of different sugars such as raffinoses, trehaloses, ketoses and kestoses and amino acids (methionine, glutamine, aspartic acid, glutamic acid) in nematode infected plants [96]. According to their analysis, the higher levels of these metabolites produced under nematode infection are mainly due to the fact that they interfere with nematode-plant signaling by reducing the incidence of nematode associations with plants [96]. An elevation in organic acid levels has also been observed in nematode infested C. annuum plants [97]. The higher levels of organic acid synthesis in plants have been linked to disease resistance against different pathogens [98]. Our results also demonstrated a further rise in the secondary metabolic profiles in nematode infested plants after amendment with PGPR. Studies have been reported that showed a modulation in the genes involved in amino acid metabolism in nematode infested A. thaliana plants inoculated with Serratia plymuthica M24T3 [99]. Furthermore, stimulation in the activities of secondary metabolic profiles have also been reported by Oh et al. [100] by Enterobacter arsburiae HK169 in nematode infected tomato plants in order to encounter the oxidative stress generated by pathogen attack. However, an upliftment in the secondary metabolic profiles have also been found in tomatoes and carrots exposed to nematode infection after the inoculation of PGPR, B. firmus T11, B. aryabhattai AO8, B. cereus, P. alvei T30 and P. barcinonensis A10 respectively [19]. The accumulation of organic acids has also been reported by inoculation of Aspergillus niger and A. candidus in plants exposed to nematode infection to form an antagonistic associations with them [101]. The metabolomic study showed exaggerated levels of secondary metabolite profiles in plants with respect to infection alone as well as with microbes. The main mechanism behind these triggered levels are attributed to activation of genes encoding enzymes involved in synthesis of different secondary metabolites.

Conclusions
Present work suggested the synergistic role of PGPR in providing rhizosphere community resilience under PPN stress. The association of microbes with plants in promoting growth, metabolic activities and upregulation of defense system to improve plant health and management of PPNs infestation has been explored. The use of microbes as biocontrol agents is the most recent and environment friendly strategy, as they improved the morphological, biochemical, physiological, and secondary metabolic profiles within plants. Moreover, they have also been known to induce ISR in plants along with the production of different toxins, antibiotics, enzymes etc. to control proliferation and penetration of nematodes within plants. They also alter the root exudations that further hinder the mobility of PPN into the soil so that they could not enter into the plant system. Consequently, it has been the most innovative approach to develop strategies for fabrication of most recent technologies in agricultural productivity in ecological manner. Nevertheless, further analysis is essential for better understanding of the mechanisms associated with plant-microbes-PPNs interactions affecting rhizospheric community, where the functions of rhizodeposition is primarily important.