Dynamics, phylogeny and phyto-stimulating potential of chitinase synthesizing bacterial root endosymbiosiome of North Western Himalayan Brassica rapa L.

The less phytopathogen susceptibility in Himalayan Brassica rapa L. has made it an exceptional crop eluding synthetic pesticide inputs, thereby guarantying economically well-founded and ecologically sustainable agriculture. The relevance of niche microflora of this crop has not been deliberated in this context, as endosymbiosiome is more stable than their rhizosphere counterparts on account of their restricted acquaintance with altering environment; therefore, the present investigation was carried out to study the endophytic microfloral dynamics across the B. rapa germplasm in context to their ability to produce chitinase and to characterize the screened microflora for functional and biochemical comportments in relevance to plant growth stimulation. A total of 200 colonies of bacterial endophytes were isolated from the roots of B. rapa across the J&K UT, comprising 66 locations. After morphological, ARDRA, and sequence analysis, eighty-one isolates were selected for the study, among the isolated microflora Pseudomonas sp. Bacillus sp. dominated. Likewise, class γ-proteobacteria dominated, followed by Firmicutes. The diversity studies have exposed changing fallouts on all the critical diversity indices, and while screening the isolated microflora for chitinase production, twenty-two strains pertaining to different genera produced chitinase. After carbon source supplementation to the chitinase production media, the average chitinase activity was significantly highest in glycerol supplementation. These 22 strains were further studied, and upon screening them for their fungistatic behavior against six fungal species, wide diversity was observed in this context. The antibiotic sensitivity pattern of the isolated strains against chloramphenicol, rifampicin, amikacin, erythromycin, and polymyxin-B showed that the strains were primarily sensitive to chloramphenicol and erythromycin. Among all the strains, only eleven produced indole acetic acid, ten were able to solubilize tricalcium phosphate and eight produced siderophores. The hydrocyanic acid and ammonia production was observed in seven strains each. Thus, the present investigation revealed that these strains could be used as potential plant growth promoters in sustainable agriculture systems besides putative biocontrol agents.


Results
Physio-biochemical characterization. A total of 200 colonies were isolated from the germplasm (roots) of B. rapa across the North-Western Himalayan zone of India (J&K) based on morphology (texture, color, elevation, margin, etc.). After ARDRA analysis, 81 strains were selected for this study. However, no colonies were on sterilized uncut roots close to the surface; the selected bacterial strains belonged to various sampling sites (Table S1). Out of the total of 81 isolated strains (Fig. 1a, Table S2) Pseudomonas sp. (n = 40, 49.38%) dominated followed by Bacillus sp. (n = 23, 28. 39%) (Fig. 1b), likewise, class γ-proteobacteria dominated (n = 47, 58.02%) followed by Firmicutes (n = 23, 28.39%) (Fig. 1c). The isolated bacterial strains showed a wide diversity across the germplasm collected from various sites. The Dominance diversity index was highest in Kansipora germplasm (1) and lowest in germplasm collected from Zakura (0.25), Simpson diversity index was 0 to 0.6939 in Kansipora and Akura, respectively. Similarly, the Shannon index was observed to be 0 and 1.475 in the sampled germplasm of Kansipora and Panzmulla, respectively. The Margalef index was 2.164 and 0 in the sampled germplasm from Zakura and Kansipora, respectively, similar to all other diversity indices [Brillouin, Menhinick, Fisher_alpha, etc.] varied widely among the isolated microflora (Fig. 1d).
The phylogenetic relationship among the isolated strains is represented in Fig. 2. A total of 22 chitinase-producing strains belonging to eight genera were further studied in context to their diversity in the collected samples and plant growth-promoting potential; among these strains, Bacillus sp. dominated, followed by Pseudomonas sp. (Fig. 3a) likewise, class γ-proteobacteria dominated (Fig. 3b). The chitinase producing bacterial strains also showed a wide diversity; the Dominance diversity index was observed to vary from 0.5 to 1, the Simpson diversity index varied from 0 to 0.5, the Shannon index was observed to be varying from 0 to 0.693, Margalef index was as high as 1.443 and as low as 0 in the sampled germplasm, in the similar way all other diversity indices [Brillouin, Menhinick, Fisher_alpha, etc.] varied widely among the isolated microflora (Fig. 3c). The population of bacterial endophytes in the sampled tissues varied from Log 10 8.386 × 10 5 to Log 10 4.227 × 10 5 (CFU/g FW) at Hugam and Nowgam, respectively (Table S1).
In the quantitative estimation of the chitinase production experiment, the highest chitinase production was observed in isolate smppsap2 (33.34 units/mL), while SKUAST-K31 was observed to produce the least (9.04 units/mL) chitinase. The CS: CZ ratio varied from 4.08 in strain smppsap2 to 1.63 in skuast2 (Fig. 4a). The species  The evolutionary history was anecdotal via the Neighbor-Joining method. The ideal tree with the sum of branch length = 26.06144331 is shown. The bootstrap consensus tree anecdotal from 1000 replicates is kept to and divisions consistent to partitions replicated in less than 50 per cent bootstrap replicates are warped. The evolutionary detachments were computed by means of the p-distance method and are in the units of the number of base variances per site. The confidence likelihood that the interior branch length is larger than 0, as projected using the bootstrap test (revealed next to the kindling). This analysis involved 81 nucleotide sequences. Codon sites encompassed were 1st + 2nd + 3rd + Noncoding. All sites with fewer than 95% site coverage were rejected, i.e., ambiguous bases were allowed at any position (partial deletion option), missing data and fewer than 5% alignment gaps. There were a total of 1427 positions in the final dataset. Evolutionary studies were conducted in MEGA X (v 10.2.6).  (Table 2). Based on 16S rRNA gene sequences of the chitinase-producing microflora, the details of their closest neighbor are given in Table 3.
Chitinase production. Based on their ability to produce chitinase, the top ten isolates were further studied to make chitinase when augmented with various carbon sources (maltose, glycerol, glucose, and sucrose) supplementations under different incubation periods (Fig. 4b). The average chitinase activity was significantly highest in glycerol supplementation, followed by glucose. The minor chitinase activity was observed in maltose. The incubation period greatly influenced chitinase activity, with the average highest chitinase activity at 192 h followed by 240-h incubation (Fig. 4c). With the increase in the incubation period, there was a corresponding increase in chitinase activity, and maximum activity was observed at the 192-h incubation period, beyond which a decline was noticed. The significantly most minor average chitinase activity among all the incubations was honored at the incubation period of 48 h, and all the isolates followed this trend. The chitinase production by all the isolates increased significantly with the incubation period from 48 to 192 h; there was a decrease beyond 192 h incubation.
Indole-3-acetic acid and other plant growth associated metabolite production. Among all the strains, only eleven produced IAA, and the highest production was observed in smppsap6 (16.55 µg/mL) and the lowest (3.13 mg/mL) in skuast2 strain (Table 4). A total of ten strains were able to solubilize tricalcium phosphate varying from 3.07 to 9.04 mm; the solubilization index ranged from 1.24 in strain SKUAST-K11 to 2.53 in strain SKUAST-K43. The phosphorus released expressed as mg/L varied from 38.23 in strain SKUAST-K67 to 264.80 mg/L in strain SKUAST-K46. The siderophore production expressed as percent unit was produced by eight strains and varied from 9.93 in strain SKUAST-K35 to 18.83 in smppsap4 strain. The HCN production described as absorbance (λ 625 nm) was created by seven strains and ranged from 0.012 in strain SKUAST-K46 to   www.nature.com/scientificreports/ 0.217 in strain smppsap4. Similarly, ammonia production was observed in seven isolates and ranged from 9.14 in strain SKUAST-K54 to 52.17 in strain SKUAST-K50 (Table 4).

Antifungal activity.
To assess the fungistatic property of the bacterial strains, six fungal species, namely Fusarium oxysporum, Fusarium solani, Dematophora necatrix 1 , Dematophora necatrix 2 Rhizoctonia solani, and Pythium amphanidermatum, were used, and fungistatic activity was expressed by calculating percent growth inhibition. All the strains exhibited fungistatic activity (Fig. 5a). A total of 19 isolates were found to be fungistatic against F. oxysporium, and the percentage growth inhibition ranged from 11.28 to 78.72 percent. Similarly, 19 isolates were observed to be fungistatic against F. solani. The percentage growth inhibition ranged from 11.39 to 59.05 percent; 21 isolates were honored to be fungistatic against R. solani. The percentage growth inhibition ranged from 11.15 to 77.54 percent; 13 isolates were observed to be fungistatic against P. amphanidermatum. The percentage growth inhibition ranged from 8.66 to 63.23 percent. 19 isolates were honored to be fungistatic against D. nectrix strain 1, and the percentage growth inhibition ranged from 7.04 to 58.90 percent. A total of 20 isolates were found to be fungistatic against D. nectrix strain 2, and the percentage growth inhibition ranged from 12.55 to 67.06 percent (Fig. 5a). The correlation studies represented by both matrix and map in Fig. 5b,c depict that chitinase activity and antifungal behavior in most isolated strains seem to be positively correlated.   Fig. 6a-e indicate that B. rapa endophytic chitinase producing microflora was primarily vulnerable to rifampicin, erythromycin, polymyxin-B, and amikacin, followed by chloramphenicol and erythromycin. While as the resistance in the isolated strains was highest towards amikacin, followed by polymyxin-B and rifampicin, respectively. Among all the strains, none was found resistant to chloramphenicol and erythromycin.
Stress tolerance. All the bacterial strains grew at all temperatures. At the temperature of 25 °C, the lowest growth was found in (Log 10 6.32 × 10 6 ) in the SKUAST-K10 strain, whereas the highest growth of (Log 10 9.41 × 10 6 ) was recorded in the smppsap4 strain against the control. At 35 °C temperature, the lowest growth was found (Log 10 6.10 × 10 6 ) in SKUAST-K53 strain, and the highest increase was recorded (Log 10 9.10 × 10 6 ) in smpp-sap4 strain, in contrast, to control (Fig. 6f). All the bacterial strains grew at all concentrations of PEG6000. At 25% PEG concentration, the smppsap4 showed the lowest (Log 10 0.2 × 10 6 ), and SKUAST-K43 (Log 10 10.7 × 10 6 ) had the highest population respectively, against the control (Fig. 6g).
In vivo studies. The inoculation experiments under pot house conditions using B. rapa as a host revealed that all the bacterial strains caused a considerable increase in both growths, yield and disease containment of the host plant (Fig. 7). The germination percentage was highest observed in the smppsap2 inoculated host (95.31 ± 0.21) against the control (83.12 ± 0.05); all the isolates caused a significant increase in root and shoot dry weights with the highest source and shoot dry weight (mg/plant) in 267.18 ± 0.03 and 1068.01 ± 0.11 against the control with 137.21 ± 0.41 and 561.02 ± 0.27 respectively. Similarly, the yield per plant and oil content (%) were significantly elevated in all the inoculated plants. The oil content (%) and seed yield (g/plant) was highest in smmppsap2 inoculated plants (34.19 ± 0.18 and 10.18 respectively) against the control with an oil content of 21.04 ± 0.14 (%) and 5.76 ± 0.09 (g/plant) respectively. All the strains prevalently reduced the disease intensity in all the inoculated plants against the control. The disease intensity was as low as 10.44% in the SKUAST-K46 inoculated strain against the control (48.09%).

Discussion
Over the past few decades, extensive use of chemical fertilizers, herbicides, and pesticides has raised environmental concerns to prevent food scarcity globally and escalate crop yield. It is one of the burning issues in agriculture and oilseed production 23 . It has become imperative to seek alternative means to contain these ecologically flawed stratagems. The beneficial bacteria that inhabit the plants can be instrumental in improving the growth of plants under both usual and challenging circumstances 24 . They can improve the plant nutrient uptake, amend the soil health, and prime plant defense to develop quality bio inoculants for a stable yield and sustainability in agriculture 1 . Hence, understanding the endophytic microfloral community associated with B. rapa has become significant, besides their role in facilitating plant growth. Plant growth-promoting bacteria proliferate the endosphere after rhizosphere colonization to endeavor beneficial effects on plant growth. Many hypotheses, including motility, plant-polymer degradation, production of phytohormones, solubilization of phosphates, enzyme production, siderophore production, and evasion of plant defenses, have been described as involved in growth promotion by these biological entities 25 . B. rapa, together with other family members, forms the world's second most important oilseed crop nutritionally; its protein profile is comparable with the soya bean and contains more S-amino acids than any other oilseed meal 26 . The huge diversity of the bacterial microflora in the Himalayan germplasm of the plants is attributed to Spatio-temporal scales besides the wide climate variation gradient stretched across the Himaliyas 1 . Because of this huge diversity on account of the above predisposing factors in the area, we investigated the endophytic bacterial microflora of the crop (Brassica rapa), and the dynamics of this microflora in chitinase producing ability as no literature about the crop is available in this context. The present investigation was undertaken to document the diversity of bacterial endosymbiosiome associated with this crop to fill this knowledge gap. The screening for chitinase production and the chitinase-producing microflora was further investigated for its potential to increase plant growth, disease containment, and stress tolerance. The sampling area was divided into three parts (Baramulla in the North, Srinagar in the center, and Anantnag in the south) depending on the presence of maximum germplasm footprints of B. rapa L. These three different sites were explored for collecting the samples to study the bacterial endosymbiosiome followed by its chitinase producing ability. The PCR amplification followed by 16S rRNA gene sequencing and phylogenetic analysis determined these isolates' molecular identity and genetic diversity. A similarity search using NCBI GenBank identified all the isolates showing 91-100% identity with reference sequences. The 16S rRNA gene sequencing is a robust technique for identifying bacteria at the genus level, differentiating between closely related bacterial species 27      www.nature.com/scientificreports/ investigations 11 , Which reported the existence of thirty bacterial endophytes in pearl millet with varying degree of biochemical and morphological attributes. It could be attributed to studying bacteria's endosymbiotic niche colonization lifestyle 1 . Chitinases with similar activity ranges have been characterized in many endophytic microfloral strains of bacteria 39,40 ; these strains produced chitinase as low as 0.5 U/mL up to as high as 30 U/mL, thereby supporting our findings strongly. In the present investigation, 22 endophytic bacterial root endophytes associated with B. rapa produced chitinase enzyme. The CS: CZ ratio varied from 1.63 to 4.08. The significance of bacterial endophytes delivering chitinase in a colloidal chitin-containing environment with similar CS: CZ has been previously investigated 41,42 , wherein the isolated microflora from multiple agronomic crop studies had the CS: CZ ratio ranging from 1-6. After optimizing various carbon source supplementations under 5 incubation periods, the average chitinase activity was significantly highest in glycerol supplementation and least in maltose. With the increase in the incubation period, there was a corresponding increase in chitinase activity, and maximum activity was observed during the 192-h incubation period. A decline was noticed at least at 48 h. Our findings indicate that bacterial endophyte chitinase is inducible and repressed by additional sugar catabolites. The increased chitinase production with an incubation period up to 192 h and then declines could be attributed to autolysis of the enzyme or catabolite repression when the threshold limit is reached; these findings are in line with various investigations 43,44 , although contrary to our findings there was huge variation among the bacterial strains, some showing peak production with maltose as well, which could be attributed to the diversity in the metabolism of the endophytic bacterial microflora. The chitinase producing bacterial strains also showed a wide variety; the Dominance diversity index was observed to vary from 0.5 to 1, the Simpson diversity index changed from 0 to 0.5, the Shannon index was honored to be ranging from 0 to 0.693, Margalef index was as high as 1.443 and as low as 0 in the sampled germplasm, in the similar way all other diversity indices [Brillouin, Menhinick, Fisher_alpha, etc.] varied widely among the isolated microflora. The variation in the diversity indices of chitinase-producing bacterial endophytic populations has been observed in a similar range. A similar magnitude of diversity indices viz., Dominance, Simpson diversity index, Brillouin, Menhinick, Fisher_alpha, and other indices were also reported in the endophytic strains of Himalayan origin Padder et al. 1 . It is attributed to selection pressure in disease escape 4 . Of all the isolated strains, only eleven produced IAA, HCN, and ammonia production in seven isolates. Dey and Raghuwanshi 45 , while investigating the comprehensive assessment of growth parameters for screening endophytic bacterial strains in Solanum Lycopersicum, observed that IAA production varied from 8.818 to 67.636 μg/ml in the isolated strains. Similarly, HCN and ammonia production was detected among the strains. A total of ten endophytic strains were able to solubilize tricalcium phosphate; the solubilization index and phosphorous released in liquid assay varied considerably among the isolates. Dey and Raghuwanshi 46 also observed the solubilization index as high as 3.714, and solubilization of tricalcium phosphate was caused by all the isolated www.nature.com/scientificreports/ endophytic bacterial strains except the one strain RR11. Phosphate solubilization has been observed in various endophytic bacteria strains in a similar range 47,48 . Furthermore, it was observed that eight strains produced siderophores, Lipková et al. 49 while visualizing the growth promotion and disease suppressing ability of endophytic strains of Brassica napus L. observed that the majority (37 out of 38) strains produced siderophores, thus validating our supporting our findings that siderophore producing endophytic bacterial strains exist in the brassicacae family. Iron-chelating microbes in complicated environments have made plants eligible for growth 46 . These plants have been found to provide appropriate niches to support the development of microflora 50 . The present investigation revealed that the isolated bacterial microflora behaved differently to the various antibiotics; it was observed that the isolated strains were primarily vulnerable to rifampicin, erythromycin, polymyxin-B, and amikacin, followed by chloramphenicol and erythromycin. The endophytic microflora has a massive diversity in context to antibiotic sensitivity. Even the candidates from the sole genus behave erratically to a specific antibiotic 51,52 . T heir behavior towards antibiotics like rifampicin, erythromycin, polymyxin-B, and amikacin has been observed similarly earlier 53 . All isolated bacterial endophytes grew in induced drought and elevated temperature mediated stress. The endophytes have been reported to withstand elevated temperatures even as high as 55 °C as low as subzero by Padder et al. 1 . The stress tolerance in the isolated microflora has revealed the more significant role of these endosymbiosiome entities shortly in changing global climate scenarios [54][55][56][57][58][59] . The stress tolerance in the bacteria can be attributed to many factors, like temperature, water, and other indigenous factors of the bacterial niche from which they are isolated 60 . Therefore, the possible reason to stress tolerance in the isolated microflora of the present investigation could be that the isolation in the present study was done where the temperature goes as high as 35 °C in the active growing season of B. rapa and during winters, the frozen environment causes the water deficit conditions and as such the associated microflora has evolved to challenge this stress.
To assess the fungistatic property of the bacterial strains, six fungal species, namely F. oxysporum, F. solani, D. necatrix 1 , D. necatrix 2 , R. solani, and P. amphanidermatum, were used. All the strains exhibited fungistatic activity except strain smppsap3. 19,19,21,13,19, and 20 isolates were found fungistatic against F. oxysporium, F. solani, R. solani, P. amphanidermatum, D. nectrix 1 , and D. nectrix 2 , respectively. Endophytic microflora, predominantly the species of Bacillus and Pseudomonas, have been found to cause the percentage growth inhibition to the extent of 90% and more in dual culture experiments with many pathogens like Fusarium sp., Pythium sp. etc. and thus protect plants against various pathogens and increase seedling emergence 61 . ' Antibiotic associations among fungi and bacteria are significant for soil and ecological functions 62 . The antifungal activity of the isolated strains can be attributed to lytic enzymes, chitinases, antibiotic synthesis, ammonia, and siderophores, among other secretory metabolites.
The stress tolerance in the isolated bacterial microflora could be autonomous entities to facilitate crop growth under exceptional and marginal environments. The attributes of endophytic microflora in this context have been reported previously 63 . The inoculation experiments under pot house conditions revealed that all the bacterial strains caused a considerable increase in the host plant's growth and yield. The germination percentage, root, and shoot dry weights increased significantly in all the inoculations against the control. Similar observations were recorded regarding the product per plant and oil content (%) in all the bacterial inoculations. The disease intensity was as low as 10.44% in the SKUAST-K46 inoculated strain against the control (48.09%). Similar to our results, Dawwan et al. 64 reported significant improvements in vegetative growth parameters of inoculated plants as compared to uninoculated control. Shi et al. 65 also found a substantial increase in fresh and dry weight plant height when inoculated beetroots with selected endophytic bacterial isolates. Padder et al. 66 reported an increase in germination percentage up to 95%, disease incidence containment [17% against the control (76.09%)], shoot length, and root length as high as 25-50% in Dalbergia sissoo inoculated with bacterial isolates due to their ability to produce a plethora of antifungal compounds like chitinase. The reduced disease incidence and intensity and stress alleviation have been the prominent attribute of bacterial endophytic formulations due to their ability to enhance the different gene expressions coupled with elevated biochemical response in the host 1 . The present investigation's increased germination percentage, yield, and growth could inhibit some prevalent seed-associated phytopathogens by cascading bacteria-associated metabolites. Chitinase seems to play a central role in phytopathogen containment 66 . This could be why the highest chitinase-producing strain smmppsap2 resulted in the highest germination percentage, ultimately alleviating the host's biotic stress and ultimately good yield and increased plant growth. Similar findings were reported by padder et al. [67][68][69] with Pseudomonas fluorescence strain smmppsap5.

Conclusion
This study offers acumens on the endophytic microbiome of B. rapa, which has not been characterized previously. The conclusions of the present investigation could be summarized as (a) B. rapa plant sustains a diverse endophytic microbiome which varies mainly between the sampling sites across the Himalayan germplasm. (b) The chitinase-producing ability in the sampled germplasm reflected dynamism in the endophytic community. (c) The presence of plant growth-promoting traits and chitinase production in the isolated culturable microfloral strains suggests that there is the possibility that these endophytes contribute to the growth, development and disease standing of the crop. (e) The stress tolerance attributes in the isolated microfloral strains suggest that these candidates can be instrumental in helping the crop to sustain the changing climate-mediated stress, both temperature and drought. (f) The controlled experiment studies reflect that there is the possibility that these strains can be applied as bioinoculant formulations at the field level. These concepts suggest that plants and endophyes evolve together within a particular niche. The various factors shape the endophytic symbiosome; these factors need to be studied more extensively to decisively determine the fundamental elements that determine the dynamics of endophytic microflora across the germplasm of B. rapa. Isolation. The surface sterilization technique was employed with some modifications 15 . Before separating the endophytes, the collected samples were washed under tap water, then dried by soaking in 70% ethanol for one minute. These dried samples were then treated with three percent sodium hypochlorite for 30 s. Lastly, subjected to washing with sterile distilled water (thrice) and then dried using tissue paper (sterilized). The samples, after getting through the process of surface disinfection, were placed for grinding with a sterile mortar pestle, and then the ground samples were dissolved in phosphate buffer saline (1000 ml containing 1.44 g Na 2 HPO 4 ; 2.4 g KH 2 HPO 4 ; 0.2 g KCl; 8.0 g NaCl and a pH of 7.4). Ground plant material was then serially diluted, and a volume of 100 µl from apiece dilution of 10 -2 , 10 -3 , 10 -6 , 10 -9 were spread plated on solid growth media (trypticase soy agar) followed by incubation at 37 °C for 24-72 h. The sterilizers' (sodium hypochlorite) surface sterilization performance was tested by inoculating surface sterilizing the whole root on tryptic soy/agar (TSA; Merck Co., Germany) media plates; the no growth indicated surface sterilization.

Scientific Reports
The morphological appearances, such as shape, color, and growth features of all the colonies, were examined. Dissimilar bacterial colonies were selected and sub-cultured. Each bacterial endophytic culture was assigned a specific code. All the isolated cultures were preserved at − 80 °C in a solution of 20% glycerol for further studies.
Diversity analysis and molecular description with 16S rRNA and ARDRA analysis. The CTAB procedure (N-Cetyl-N, N-trimethyl-ammonium bromide) was used to isolate the total genetic Material [70][71][72] . For 16S rDNA amplification, both forward and reverse primers viz., fD1 (5′ AGA GTT TGA TCC TGG CTC AG 3′) and rD1 (5′ AAG GAG GTG ATC CAG CCG CA 3′) 73 were used in a thermocycler (BioRad, USA). Amplification materials were resolved on a 1.5 percent agarose gel and interpreted using a gel doc framework (Alfa Imager, Alfa Innotech Corporation, USA). Genie PureTM fast PCR purification kit (GeNeiTM, Bengaluru, India) was used to purify the amplicons, then quantified at 260 nm using a spectrophotometer using calf thymus DNA as a control. Restriction fragment length polymorphism of amplified 16S ribosomal DNA with RsaI and HhaI (Thermo scientific) employed the method 1 to avoid duplicity among the isolated bacterial strains. In an Applicable Biosystems 3130 Genomic Analysis tool, the purified partial 16S rDNA amplification were sequenced (Applied Biosystems, CA, and the USA).

Analysis of 16S rDNA sequences. The Nucleotide Basic Local Search Tool (BLAST N) Program (www.
ncbi. nlm. nih. gov/ BLAST)' was used to compare the partial nucleotide sequences with already existing arrangements in the NCBI database, the lines with > 99% similarity were retrieved. All the deals were trimmed for the Neighbor-Joining Tree construction, engaging trim primers for mapped reads (the reads hitting the primer parameters using primer track) in the CLC package (QUAGEN Bioinformatics). All the sequences were deposited in NCBI at https:// submit. ncbi. nlm. nih. gov/ about/ genba nk. The datasets generated in the current investigation are available in NCBI at https:// www. ncbi. nlm. nih. gov with accession numers (Table S2). Phylogenetic analysis and multiple cluster alignment were carried out using the MEGA package (v.10.2.6) founded on the neighbor-joining method using a one-kilo repetition bootstrap. The diversity analysis of the isolated taxonomic entities was estimated employing Past 3 software (v 4.03) and EstimateS (v.9.1.1).
Chitinase activity. To screen the isolated bacterial microflora for chitinase activity, the bacterial strains after serial dilution were streak plated on (15μL of 3 × 10 9 CFU/mL) fresh medium [2.0 g colloidal chitin (crustacean); 0.1 g KH 2 PO 4 ; 0.01 g MgSO 4 .7H 2 O; 3.0 g NaCl; 0.7 g (NH 4 ) 2 SO 4 ; 0.05 g yeast extract; 2.0 g agar and 50 mM of sodium phosphate buffer, pH 6.0 (all the components are given on per liter basis)] at 30 °C over three days, after incubation, the bacterial strains were examined for the development of clearance zone on the solid media, and its presence indicated the presence of chitinase activity. The clearance zone and colony diameter measurements were made to calculate the CS: CZ ratio. To determine chitinase activity (U/mL), Berger and Reynold's method calculated chitinase activity 74 . The same media (100 mL) except the agar was taken in 250 ml Erlenmeyer flasks and inoculated with each chitinase positive bacterial strain (15μL of 3 × 10 9 CFU/mL log-phase cells) and incubated at 30 °C for 72 (hrs.). After incubation, the enzyme mixture was collected by centrifuging the mixture for 20 min at 12,000 rpm at room temperature. The substrate and enzyme mixture [0.5 ml (1% w/v of colloidal chitin in minimal media) and 0.5 ml of enzyme solution] was incubated at 45 °C for 1 h. The reaction was then stopped by 3 ml of a 5-dinitrosalicylic acid reagent and then 'heated for five minutes at 100 °C to cease the reaction completely. ' The supernatant (collected by centrifuging the heated mixture at 10,000 rpm for 15 min at 4 °C) was collected and subjected to absorbance measurement with a UV (530 nm) spectrometer. For determination of enzyme unit, serial dilutions of Nacetylglosamine (from 0 to 50 mM) were prepared. One unit (U) of the chitinase activity was defined as the amount of enzyme required to release one mmol of N-acetyl D glucosamine (as a standard) from chitin/min. The experimental protocol was modified as per the desired set of conditions to carry out experiments on the effect of various carbohydrate sources and different incubation periods.
Biochemical characterization of bacterial isolates. All chitinase-producing endophytic bacterial isolates were characterized based on biochemical features as per Bergey's manual of determinative Bacteriology 75 .
Indole-3-acetic acid (IAA) production. The Salkowski method of color development was employed for calculating IAA 76 ."The bacterial endophytic strains were inoculated into 25 mL nutrient broth fortified with www.nature.com/scientificreports/ 0.2 mg/mL L-tryptophan and incubated in a BOD incubator at 28 °C. Following incubation of 72 h, 2 mL per cultivation broth was centrifuged for two minutes at 7,000 rpm, followed by the addition of the same amount of Salkowski's reagent (1 mL 0.05 M FeCl 3 in 50 mL 35 percent HClO 4 ) to cultivation supernatant for the colored complex production which was estimated calorimetrically at 500 nm wavelength using a spectrophotometer.
Phosphate solubilization. The TSA broth-developed bacterial strains were placed [3 × 10 9 CFU/mL log phase developing cells (15μL)] on Pikovskaya's medium plates for the detection of solubilization zone 77 . These plates were incubated for 4 days at 28 °C. The established solubility region and colony diameter measured were used to calculate the solubilization index (SI = diameter of the colony/ diameter of halozone/colony diameter) 78 . Bray and Kartz were employed using Pikovskaya's broth 79 . The cultures inoculated in Pikovskaya's broth [3 × 10 9 CFU/mL log phase developing cells (15μL)] were incubated for 4 days at 28 °C. The soluble phosphorus released in the broth was estimated calorimetrically.
Production of ammonia and hydrocyanic acid. The log phase developing cells (15μL of [3 × 10 9 CFU/ mL] were inoculated in peptone water (10 mL) and incubated for (72 h at 30 °C). Approximately (0.5 mL) Nessler reagent was added, and the ammonia output's yellow/brown pigment appearance was recorded 80 . The culture was inoculated after incubation was centrifuged for 15 min at 10,000 rpm for quantitative estimation. To the collected supernatant (1 mL), a volume of 1 mL of the Nessler reagent was added. To this mixture, purified water was added for serial dilution because of the sensitivity of the calculation method to higher ammonia concentrations. Brown/yellow color was developed, and absorbance was measured at 630 nm' . The ammonia volume was determined by the ammonium sulfate level (0.1-1 mol/mL) curve 81 . The King's B medium was mixed with 4.4 g/L of glycine to measure the hydrocyanic acid (HCN) 82,83 . 'The degree of color shift from yellow to dark brown on filter paper dipped into 2% sodium carbonate (prepared in picric acid) was determined. The bacteria (15μL [3 × 10 9 CFU/mL were grown in the glycine-modified King's B broth, and the regular filter paper strip (100 × 15 mm 2 ) in alkaline picrates was soaked and hanged for quantitative testing in the conical flask 84 . After the filter paper was inserted into a test tube containing 10 mL of filtered water, the absorption of this colored water was estimated at 625 nm.
Siderophore estimation. The liquid Chrome Azurol S (CAS) assay was used for siderophore estimation 85 .
A 0.5 mL of Chrome Azurol-S (CAS) was added to 0.5 mL of supernatant (1 mL) of the cell-free extract [obtained by inoculating 15μL of 3 × 10 9 CFU/mL log-phase cells into the nutrient broth and after incubation was centrifuged for 15 min at 10,000 rpm] besides 10 µL shuttle solution (0.2 M 5-Sulfosalicylic acid). The absorption range was estimated at 630 nm after 10 min of incubation at room temperature. The media (minimal) was taken as the reference (made with the same ingredients as the minimal medium, except for the addition of cell-free culture-supernatant extract). Siderophore units were obtained from the following equation: where [Ar is the reference absorption at (630 nm)].
Antifungal activity. All of the chitinase-producing bacterial strains were tested against 6 phytopathogenic strains :(Fusarium oxysporum, Fusarium solani, Pythium amphidermatum, Dematophora necatrix 1 , Rhizoctonia solani, and Dematophora necatrix 2 ). The bacterial endophytes (48 h. old culture) full-grown on nutrient agar Petri dishes were placed (as line streaks) on PDA at an equal distance of 2.2 cm from the center of the 100 × 15 mm Petri dish. A 5 mm disc taken from the 10-day old fungal culture was placed at the center of the Petri dish as per the protocol 86 with some modifications. The plates were incubated at 23 °C. The plates without the bacterial endophyte were kept under control. The experiment was repeated per the required replications set for the experimental design. 'The formula calculated the percentage growth inhibition: Antibiotic sensitivity. Antibiotic susceptibility testing was carried out using discs impregnated with antibiotics (6 mm diameter). The strains' antibiotic susceptibility was assessed using the Kirby Bauer disc-diffusion assay system against chloramphenicol, rifampicin, amikacin, erythromycin, and polymyxin-B (analytical grade) 87 . ' Antibiotics were used in amounts of 35.0 µg per disk. According to the DIFCO Manual 10th edition, strains were classified as resistant or responsive based on the inhibition zone reported 88 .
Stress tolerance. The bacterial isolates were grown at different PEG concentrations (10, 15, 20, and 25%) and at different temperature ranges (15, 25, and 35 °C) in nutrient broth as per the local climatic conditions to observe their ability to grow in possible extreme conditions of the zone; the CFU/mL was measured. Formulation Development.
In vivo studies. The effect of isolated cultures on the growth of B. rapa was ascertained under pot culture conditions by using it as a host plant (B. rapa). Seeds of B. rapa were surface sterilized by 0.2% sodium hypochlorite. They were sown in pots of 4 kg capacity containing 2 kg silty clay loam soil and inoculated with individual Percent Siderophore Unit = (Ar − As/As = 100) .
Control − Treatment/Control × 100 www.nature.com/scientificreports/ isolates. All the pots were inoculated with 3 ml inoculum of log-phase bacterial isolates every 30 days till the harvesting of the crop. A control was also maintained by inoculating with broth devoid of bacterial culture. The experiment was conducted in a completely randomized design, with each treatment replicated four times. The pots were irrigated as and when needed. At the crop's maturity, the plants were uprooted, and observations were recorded for various growth and yield attributes, viz., oil content, yield per plant, dry root weight, and shoot dry weight. The yield per plant was calculated by collecting the seeds per plant and weighing them. The germination percentage was calculated as per the method adopted by 89 . Oil content was determined by following the solvent extraction technique 90 , for which 3 g of B. rapa seeds were first crushed in 3 g Na 2 SO 4, and the resultant powder containing oil was taken into test tubes; then 20 ml hexane was poured into the test tubes as mobile phase. Elute-containing oil was stored in a vile, and hexane was evaporated in a hot water bath. The remaining oil was weighed, and its percentage was calculated using the following formulae: Measurement of disease intensity. The Sclerotinia sclerotiorum infected plant samples of B. rapa were collected from the plains of J&K, and the cultures were introduced in the plants grown under controlled conditions; the disease incidence was recorded every three weeks by counting the number of infected leaves, and visually estimating their per-centage diseased area. Observations were made for five randomized branches per plant, and ten older leaves were evaluated to assess the disease incidence and severity using the following formula: