The first part of this research was conducted in the chemistry and biology laboratory of the Agroengineering Institute. Subsequently, for agronomic management, a “vulture wing” type greenhouse was used. Both are located at the facilities of the Universidad del Papaloapan, Loma Bonita Campus, in the municipality of Loma Bonita, Oaxaca, Mexico, with coordinates 18°05′56.5″ N and 95°53′48.0″ W, at an altitude of 38 m above sea level.

A randomized block experimental design was conducted with four treatments, each with 20 repettions, equivalent to 20 plants per treatment. The treatments included were: Q (UNPA Chitosan + Nutrient Solution), BS (Bacillus subtilis), Q + BS (UNPA Chitosan + Bacillus subtilis), and T (Nutrient Solution), as detailed in Table 1.

The treatments were applied at the time of sowing and then at 15, 30, and 45 days after sowing, as well as at 15, 30, 45, 60, 75, 90, 105, and 120 days after transplanting (DAT). It is important to mention that the transplanting was done when the plants were 45 days from sowing. The bioproducts were applied at a dose of 15 mL per plant using the drench technique. As a solvent for chitosan, slightly acidified distilled water (pH 5.6) with 1% v/v glacial acetic acid was used to facilitate complete dissolution and ensure stability. In all treatments, the standard Hoagland and Arnon (1950) nutrient solution was used, composed of macronutrients (NO₃⁻, NH₄⁺, PO₄³⁻, K⁺, Ca²⁺, Mg²⁺, SO₄²⁻) and micronutrients (Fe, Mn, Zn, Cu, and B), adjusted to a pH of 6.0 ± 0.2.

The Bacillus subtilis strain was provided by the Bioprocess Laboratory of the University of Papaloapan (UNPA) at the Loma Bonita campus. The lactic acid strain was reactivated in Man Rogosa Sharpe (MRS) medium, incubating at 38°C for 48 h in an incubator (Felisa^®, FE132AD, Ciudad de México, Mexico). The resulting bacterial suspension was adjusted to a concentration of 1 × 10⁹ CFU·mL⁻¹ [13].

Two grams of chitosan were dissolved in one liter of distilled water acidified with 1% v/v glacial acetic acid and kept under constant agitation for 15 min to ensure complete dissolution. For Bacillus subtilis, 15 mL of the bacterial concentrate (1 × 10⁹ CFU·mL⁻¹) were added to one liter of distilled water and maintained under constant agitation for 15 min. Finally, for the combined product, two grams of chitosan and 15 mL of the Bacillus subtilis concentrate were added to one liter of acidified water, keeping constant agitation for 15 min.

Data were processed using InfoStat 2017 for Windows. For variables measured repeatedly over time (plant height, stem diameter, number of leaves, chlorophyll a and b, sap electrical conductivity, pH, and sap nutrient concentrations), a repeated-measures ANOVA was performed to evaluate the effects of treatment, time, and their interaction. For fruit-quality variables (morphological and physicochemical traits) evaluated across three fruit sizes (small, medium, and large), a two-way ANOVA (treatment × fruit size) was applied. When significant effects were detected (p ≤ 0.05), means were separated using Tukey’s post-hoc test. Graphs were generated with SigmaPlot 12.00 for Windows. Values are reported as means ± SE (n = 20 plants for morphological and yield traits; n = 3 replicates for physiological and chemical analyses).

The results of the analysis of shrimp exoskeleton samples, UNPA chitosan, and Sigma chitosan are presented in Table 2, showing the contents of minerals (% M), Demineralization (% DM), Deproteinization (% DP), proteins (% P), moisture (% H), and pH (hydrogen potential) in each sample.

From the analyses performed, differences were found between the samples in terms of their chemical composition and physical properties. The shrimp exoskeleton shows a high content of minerals (24.6%) and proteins (44.62%), making it a rich source of these components. This composition is expected, considering the structural and protective function of the exoskeleton in shrimp. Additionally, the low moisture level (6.5%) and slightly alkaline pH (8.8) are consistent with the typical characteristics of marine biological materials. On the other hand, both UNPA chitosan and SIGMA chitosan show much lower contents of minerals and proteins, indicating a process of demineralization and deproteinization during their production. These results suggest efficient removal of impurities and undesirable components to obtain high purity chitosan. However, a notable increase in moisture content is observed in the chitosans (11.5% and 11.33%, respectively), which can be attributed to the extraction and purification process. Furthermore, both chitosans exhibit an acidic pH (3.4), probably due to the acidic nature of the deproteinization processes. The differences in chemical composition and physical properties between the shrimp exoskeleton and the derived chitosans have important implications for their industrial applicability. While the shrimp exoskeleton can be valued as a natural source of minerals and proteins in the food and supplement industries, UNPA and SIGMA chitosans are more suitable for agricultural applications due to their purity and specific physicochemical properties. Finally, these results highlight the importance of understanding the chemical composition and properties of biological materials and their derivatives for effective use in various industrial and nutritional applications. The results of the colorimetric analysis for the shrimp exoskeleton, UNPA chitosan, and SIGMA chitosan are presented in Table 3. The color coordinates L* (luminance), A* (red/green axis), and B* (yellow/blue axis) are shown, along with the color description and corresponding color chart.

These results reveal significant differences in color profiles among the analyzed samples. The shrimp exoskeleton exhibits a “pale grayish medium tone” with L* (luminance), A* (red/green axis), and B* (yellow/blue axis) values of 62.23, 1.17, and 16.3, respectively. This color may be characteristic of the chemical composition and structure of the exoskeleton, with the presence of natural pigments that give it this specific hue. On the other hand, both UNPA chitosan and SIGMA chitosan show lighter and less saturated colors compared to the shrimp exoskeleton. UNPA chitosan presents a “grayish medium white tone”, with L*, A*, and B* values of 75.63, 0.56, and 3.46, respectively. Meanwhile, SIGMA chitosan exhibits a “medium brown tone”, with L*, A*, and B* values of 65.56, 4.8, and 12.1. These colors may be influenced by the extraction and purification process of chitosan, as well as by the chemical composition of the starting materials.

The difference in color profiles between shrimp exoskeleton and derived chitosans may be important in terms of their application in various industries. For example, the lighter and neutral color of chitosans may make them more suitable for use in applications where the color of the material is an important factor, such as in the food or pharmaceutical industries. Meanwhile, the characteristic color of shrimp exoskeleton may be useful in applications where the aim is to harness its natural and authentic properties, such as in the manufacturing of cosmetic products or dietary supplements. These results highlight the importance of considering the color profile along with other physical and chemical properties when selecting and using biological materials and their derivatives in various industrial and research applications.

An FT-IR analysis was performed to confirm the structure of UNPA chitosan, comparing it with Sigma chitosan by identifying the characteristic bands of this polymer. Fig. 1 shows the obtained FT-IR spectra, where the O–H and N–H groups are observed at 3445 and 3280 cm⁻¹, respectively. The band corresponding to the C–N group is found at 1554 cm⁻¹, while the amide III band is located at 1377 cm⁻¹. The bands for the C=O linkages of amide I are identified at 1630 cm⁻¹, and the stretching bands N–H for amide I are found at 3260 cm⁻¹. For Sigma chitosan, the most representative bands of the biopolymer are found at 3438 cm⁻¹, corresponding to the N–H stretching for amide I and the broad vibration of the hydroxyl groups. Another characteristic band of the carbonyl group is observed at 1633 cm⁻¹. It is important to note that these bands are presented with lower intensity compared to UNPA chitosan.

In Fig. 2, the results of the morphological parameters in seedlings are presented, evaluated at 15, 30, and 45 days after sowing (DAS), aiming to measure plant height (A), stem diameter (B), and number of leaves (C). The results revealed that at 15, 30, and 45 DAS, the combined treatment of Q + BS showed greater height, with means of 6.88 ± 0.10, 10.82 ± 0.26, and 17.54 ± 0.46 cm, respectively. However, at 15 DAS, no significant differences were observed among the treatments T, Q, and BS. Similarly, at 30 and 45 DAS, there were no significant differences between treatments Q and T, but significant differences were detected in the treatment with BS, which presented the lowest height, with means of 8.67 ± 0.21 and 14.22 ± 0.40 cm, respectively. Regarding stem diameter, the most effective treatment at 15, 30, and 45 DAS was the combination of Q + BS, with means of 1.49 ± 0.10, 1.55 ± 0.10, and 2.20 ± 0.10 mm, respectively. However, at 15 DAS, no significant differences were found between treatments Q and BS, with treatment BS presenting a lower diameter, showing no significant differences with treatment T. At 30 DAS, no significant differences were found between treatments Q and T either. Similarly, at 45 DAS, no significant differences were found among all treatments.

In relation to the number of leaves, no significant differences were observed among the four treatments at 15, 30, and 45 DAS. This phenomenon suggests that while the combination of Q + BS showed advantages in height and stem diameter at certain time points, it did not significantly influence the number of leaves. These results highlight the importance of considering the combination of different treatments to maximize growth and development in seedlings under specific conditions.

In Fig. 3, the morphological growth parameters of Soledad chili plants are presented, evaluated at 15, 30, 45, 60, 75, 90, and 105 days after transplanting (DAT). These parameters include plant height (A), stem diameter (B), leaf area (C), and the number of shoots (D). It was observed that for the height variable, no significant differences were found between treatments at 15, 30, 45, 90, and 105 DAT. However, at 60 DAT, the treatment with BS showed a greater height (120.19 ± 2.68 cm), distinguishing itself from the control and the combination of Q + BS, although it did not show significant differences with the application of Q, which reached a height of (114.8 ± 2.33 cm). Similarly, at 75 DAT, the combination of Q + BS presented the greatest height (137.50 ± 3.65 cm), with no significant differences with the BS treatment (133.60 ± 3.19 cm), but significant differences compared to the other treatments. Regarding stem diameter, the control showed significant differences compared to the other treatments (Q, BS, and Q + BS), with the combination of Q + BS exhibiting the greatest diameter (0.307 ± 0.006 cm). Furthermore, no significant differences were observed between treatments and the control at 30, 45, 60, 75, 90, and 105 DAT.

In relation to leaf area, the best treatments at 15 DAT were Q and Q + BS (33.65 ± 1.77 and 32.16 ± 2.01 cm², respectively), with no significant differences between them, but significant differences when compared to the control and BS. At 30 DAT, the greatest leaf area was obtained with the application of Q (59.57 ± 3.49 cm²), followed by the control, with no significant differences between these two treatments. However, the control did not show significant differences concerning the other treatments (BS and Q + BS). In the following time intervals, significant differences between treatments varied, highlighting the superiority of the combined treatment Q + BS at 75 DAT. Regarding the number of shoots, at 15 DAT, the application of Q and Q + BS showed better results (7.05 ± 0.29 and 6.25 ± 0.29 shoots, respectively), with no significant differences between them, but significant differences compared to the other treatments. In summary, variability was observed in the performance of the treatments throughout the study period, emphasizing the effectiveness of the combination of Q + BS in various aspects of the growth of Soledad chili plants. These results may be useful for optimizing cultivation practices and improving the yield of this particular crop.

In Fig. 4, the results obtained for photosynthetic pigments at 30, 60, 90, and 105 days after transplanting (DAT) are shown for chlorophyll A (A), B (B), and chlorophyll A y B (C). It was observed that for chlorophyll A at 30 DAT, no significant differences were found between the treatments compared to the control treatment (T). However, at 60, 90, and 105 DAT, the treatment with Q showed higher concentrations of chlorophyll A (5.97 ± 0.35, 2.41 ± 0.03, and 2.14 ± 0.07 μg mL⁻¹, respectively).

At 60 DAT, the control (T) did not show a significant difference compared to the application of Q, nor did the application of BS and the combination of Q + BS compared to the control, with the BS application presenting the lowest levels of chlorophyll A at a concentration of 2.15 ± 0.85 μg mL⁻¹. At 90 DAT, the application with BS showed no significant differences compared to T, with the combination of Q + BS presenting the lowest concentration of chlorophyll A (1.10 ± 0.17 μg mL⁻¹). At 105 DAT, no significant differences were found between the treatment with Q and the control; however, Q showed higher chlorophyll values but did present differences compared to the treatments with BS and the combination Q + BS.

Regarding chlorophyll B, the best results at 30 and 60 DAT were obtained with the application of Q, showing concentrations of (4.13 ± 0.42 and 6.74 ± 0.55 μg mL⁻¹, respectively). For 90 and 105 DAT, the highest concentration was observed in the combination of Q + BS (0.73 ± 0.08 and 13.87 ± 0.82 μg mL⁻¹, respectively), followed by the application of BS with 3.34 ± 0.20 μg mL⁻¹, which showed no significant differences. The application of BS showed no differences compared to T at 60 and 90 DAT, with a concentration of 2.69 ± 0.03 μg mL⁻¹, nor with the combination of Q + BS, which presented the lowest concentration of chlorophyll B (1.96 ± 0.33 μg mL⁻¹). It was observed that the combination of Q + BS and the individual application of Q showed promising results in the concentrations of chlorophyll A and B at different evaluation times. These findings suggest the potential of these treatments to enhance photosynthesis and, consequently, the growth and development of Chile Soledad plants under specific conditions.

The analysis of sap in the petiole allowed establishing the relationship between the application of chitosan, Bacillus subtilis, and the combination of UNPA chitosan with Bacillus subtilis in the variables of electrical conductivity (μs cm⁻¹) (A), pH (B), and Brix degrees (C) at 30, 60, 90, and 105 days after transplanting (DAT), as shown in Fig. 5. It was determined that at 30, 60, and 105 DAT, there were no significant differences in electrical conductivity among the treatments compared to the control (T); however, at 90 DAT, the treatment that presented the highest electrical conductivity was T (600 μs cm⁻¹), which showed a significant difference compared to the other treatments, followed by the application of Q (453 ± 17.55 μs cm⁻¹). Regarding the treatment with BS, it presented a conductivity of (364.66 ± 2.90 μs cm⁻¹), while the combined treatment Q + BS had the lowest electrical conductivity (246.33 ± 11.34 μs cm⁻¹). Concerning the pH values recorded (B), from 30 DAT, the pH of the crop did not show significant differences among the treatments, remaining above 5 ± 0.40; at 60 DAT, the treatments with the highest pH were BS (6 ± 0.40) and T (6.1 ± 0.05), with no significant differences between them. At 90 DAT, the BS and T treatments presented a pH of 6.1 ± 0.03, with no significant differences between them, while Q had a pH of (4.7 ± 0.05) and Q + BS (3.06 ± 0.03). Finally, at 105 DAT, the treatments with the highest pH were Q + BS and T (5.5 ± 0.05), with no significant differences between them, nor between Q and BS (4.76 ± 0.12) and (5.03 ± 0.03), respectively. Lastly, regarding the concentration of Brix degrees (Fig. 6C), at 30, 60, 90, and 105 DAT, the treatment with the highest concentration was observed with the application of Q, with values exceeding 3 ± 0.03. However, at 30 DAT, there were no significant differences compared to the combined treatment Q + BS, which showed a value of 2.53 ± 0.06. At 60 DAT, significant differences were observed among all treatments, where BS, T, and Q + BS had values of 3.76 ± 0.06, 3 ± 0.05, and 2.26 ± 0.03, respectively. At 90 DAT, the combined application Q + BS presented values of 3.1 ± 0.01, followed by the control with values of 2.7 ± 0.0, while BS had the lowest values of 2.1 ± 0.0. Similarly, at 105 DAT, it was determined that the highest value, followed by Q, was T with values of 3.76 ± 0.0; however, for the combined application Q + BS and BS, no significant differences were observed (2.16 ± 0.12 and 1.8 ± 0.05).

The results of the concentration of nitrates, phosphates, and potassium in sap (mg L⁻¹) evaluated at 30, 60, 90, and 105 DAT for each treatment are shown in Fig. 6. At the end of the evaluation, it was determined that the highest concentration of nitrates at 30, 60, 90, and 105 DAT was found in the treatment with Q, showing values exceeding (1000 ± 33.33 mg L⁻¹); however, at 30, 60, and 90 DAT, no significant differences were observed with the treatment with BS, which showed values exceeding (900 ± 33.33 mg L⁻¹). However, at 60 and 90 DAT, no significant differences were shown between treatments compared to the control, but differences were noted at 105 DAT. Regarding the concentration of phosphates present in the sap, it was determined that the highest concentration at 30, 60, 90, and 105 DAT was found with the application of BS, showing values of (820 ± 20.81, 430 ± 51.96, 992.66 ± 3.17, and 490 ± 40.41 mg L⁻¹), respectively; however, at 60 DAT, no significant difference was shown compared to the application with Q and the combination of Q + BS. Regarding the concentration of potassium, it was determined that at 30, 60, and 90 DAT, no significant differences were observed concerning the treatments, but at 90 DAT, the combination of Q + BS showed a higher concentration (1566.66 ± 33.33 mg L⁻¹); however, it did not present differences compared to treatments Q and BS.

Table 4 shows the repeated-measures ANOVA for the growth, physiological, and nutritional pa-rameters of Capsicum annuum L. evaluated over time. Significant main effects (p < 0.05) of both treatment and sampling time were detected for all analyzed variables, whereas the treatment × time interaction was not significant for most parameters. The combined treatment Q + BS (UNPA chitosan + Bacillus subtilis) resulted in the highest mean values of plant height, stem diameter, chloro-phyll a and b contents, and sap electrical conductivity. These results indicate a synergistic effect between the bioelicitor and the rhizospheric bacterium, which together enhanced photosynthetic efficiency and nutrient uptake throughout plant development.

The following presents the results obtained from the number of fruits per treatment for each cut made in the experiment. The first cut was made 69 days after transplanting (DAT), the second at 85 DAT, and the third and final cut at 106 DAT. The results are detailed in Table 5. In the first cut, the combination of Q + BS showed the highest number of fruits (22.0 ± 3.08 fruits), while the Q treatment presented the lowest number (16.0 ± 3.47). However, no significant differences were found between the treatments. In the second harvest, the highest number of fruits was recorded with the application of BS (49.80 ± 5.78 fruits), but no significant differences were observed between the Q and Q + BS treatments, although there were differences compared to the control (T), which had the lowest number of fruits. For the last harvest, the highest number of fruits was observed in the Q application (41.55 ± 4.11 fruits). Similarly, no significant differences were found between the Q + BS and BS treatments, but differences were noted compared to the control, which showed the lowest number of fruits per plant. The application of the biofertilizers did not show significant differences in the total weight of the fruit in cuts 1 and 3. In contrast, in cut 2, significant differences were observed between treatments compared to the control, with the Q application showing the highest fruit weight at 110.74 ± 11.07. These results are shown in Table 6.

The analysis of the effect of biofertilizers on the accumulated yield per plant from the three cuts performed is presented in Fig. 7. In the representation of the number of fruits by treatment (A), it is highlighted that the highest number of fruits was obtained with the application of Q (104.4 ± 6.02 fruits). However, no significant differences were found compared to the treatments with BS and Q + BS, which presented averages of (101.75 ± 7.56 and 91.89 ± 5.19 fruits), respectively.

Significant differences were observed regarding the treatment T, which showed the lowest number of fruits (57.75 ± 4.65). As for the accumulated weight per plant in grams (B), the highest amount was evidenced with the application of Q (246.86 ± 10.66 g). However, no significant differences were found compared to the application with Q + BS (228.04 ± 41.76 g). In contrast, the application with BS yielded (201.37 ± 15.39 g) per plant, with no significant differences compared to the control, which showed the lowest fruit weight. In summary, the results show that the application of Q led to a higher number and weight of fruits per plant compared to the treatments with BS and Q + BS. Additionally, a significant difference in fruit yield was observed between the treatments and the control, highlighting the positive impact of the application of biofertilizers, particularly Q, on the yield of Soledad chile.

Next, Fig. 8 presents the data obtained on the morphology for fruit quality from the first harvest. Variables such as fruit length (cm) (A), fruit diameter (cm) (B), fruit weight (g) (C), pericarp thickness (cm) (D), and number of locules (E) in three different sizes (small, medium, and large) were analyzed. It was determined that the fruit length in all three sizes was greater when the combined treatment Q + BS was applied, with averages of (5.32 ± 0.1, 6.29 ± 0.2, and 7.72 ± 0.2 cm) for the small, medium, and large sizes, respectively.

However, in the small size, there were no differences compared to the treatment with BS, and in the medium size, no significant differences were observed among the other treatments. Regarding fruit diameter, the largest diameter was shown in the application with Q + BS, with values of (8.17 ± 0.3, 8.36 ± 0.3, and 9.56 ± 0.3), respectively for small, medium, and large sizes. In terms of fruit weight, in the small size classification, the heaviest fruits were found in the application with Q, followed by BS (2.15 ± 0.1 and 1.87 ± 0.1 g), respectively, with no significant differences between both treatments. In the large size classification, the combination of Q + BS showed the greatest fruit weight, followed by the application with Q. Regarding pericarp thickness, in the small, medium, and large classifications, a greater thickness was observed with the combined treatment Q + BS; however, in the small size, no significant differences were found compared to the treatment with Q, but there were differences compared to the BS treatment and the control. In summary, the results indicate that the application of the combined treatment Q + BS showed significant improvements in various aspects of fruit morphology and quality, such as length, diameter, weight, and pericarp thickness, compared to the individual treatments and the control.

Fig. 9 presents the data obtained on the number of locules (E1) and the physicochemical quality variables in the fruits from the first harvest: degrees Brix (F1), percentage of citric acid (G1), and pH (H1). Regarding the number of locules, no significant differences were observed in the small and medium classifications, but significant differences were noted in the large category. The treatment with Q showed the highest number of locules, exceeding 5; however, it did not differ significantly from the treatments with BS and Q + BS, but did differ from the control treatment, which had an average of 4 locules.

In the physicochemical determinations of the fruit, it was observed that for the small size classification, the highest concentration of degrees Brix was found in the control (3.15), while Q and the combined treatment Q + BS showed lower concentrations (1.2 and 1.3, respectively), with no significant differences between both treatments. In the medium size, the highest concentration of degrees Brix was observed in the treatment with Q (3.8), showing significant differences compared to the other treatments, while BS and Q + BS had the lowest concentrations (1.9 and 2, respectively), with no significant differences between them. In the large classification, the highest concentration of degrees Brix was observed in Q and the control (2 and 2.3 ± 0.1, respectively), with no significant differences between them, but significant differences were observed compared to BS and the combined Q + BS.

Regarding the percentage of citric acid, in the small classification, treatment T showed the highest concentration (0.05%), with significant differences compared to the other treatments, while Q and Q + BS had the lowest concentration (0.01%), with no differences between them. For the medium size, the highest concentration was found in the combination Q + BS (0.04%), showing significant differences compared to the other treatments, but no differences between Q and BS (0.01 and 0.02%, respectively). In the large classification, the highest concentration of citric acid was observed with Q (0.049%), showing significant differences compared to the other treatments. When determining the pH, it was observed that in the small size classification, the treatment with Q presented a higher pH (5.9), with significant differences compared to the other treatments.

However, no differences were found among BS, Q + BS, and T. For the medium size, the highest pH was found in the treatments BS and Q (5.9), with no significant differences between them, but significant differences were noted compared to the other treatments, where Q + BS and T showed lower pH (5.7). In the large classification, the highest pH was presented in the combined treatment Q + BS (5.9); however, it did not differ significantly from the treatments with Q and BS (5.8), respectively, but it did differ from the T, which showed a pH of 5.6. The results indicate that the combined treatment Q + BS shows significant improvements in several physicochemical quality aspects compared to the individual treatments and the control, particularly highlighting the number of locules, degrees Brix concentration, percentage of citric acid, and pH, underscoring its effectiveness in improving fruit quality.

Fig. 10 presents morphological data on fruit quality from the second harvest. Variables analyzed include fruit length (cm) (A2), fruit diameter (cm) (B2), fruit weight (g) (C2), pericarp thickness (cm) (D2), and the number of locules (E) across three size classifications (small, medium, and large). For the fruit length variable, no differences were observed among treatments in the small and medium classifications. In the large classification, the combined treatment of Q + BS showed the greatest length (7.3 ± 0.2 cm), but no significant difference compared to the Q and BS treatments, which measured 6.9 and 6.7 ± 0.2 cm, respectively. However, both treatments differed from the control treatment (T), which had the shortest length (6.47 ± 0.2 cm).

Regarding fruit diameter, the highest diameter in the small classification was found in the Q treatment (0.802 ± 0.02 cm). There were no differences compared to the combined treatment Q + BS (0.692 ± 0.04 cm), but differences were noted between these and the BS and T treatments, although there were no significant differences between T and BS. In the medium size classification, no significant differences were found among the four treatments. Finally, in the large classification, the largest fruit diameter was in the combined treatment Q + BS (0.916 ± 0.03 cm), although it did not differ significantly from Q and BS treatments, but did differ from the control.

For fruit weight in the small classification, the BS treatment recorded the highest weight (2.02 ± 0.2 g), followed by the Q treatment (1.91 ± 0.1 g), with no significant differences between these treatments. In the medium classification, there were no significant differences among all treatments. In the large classification, the greatest fruit weight was noted in the Q + BS, Q, and BS treatments (4.11 ± 0.2, 4.02 ± 0.2, and 3.78 ± 0.1 g, respectively), showing no differences among these treatments but differing significantly from the control. Regarding pericarp thickness, no differences were found among treatments in the small and medium classifications. However, in the large classification, the BS treatment showed the greatest pericarp thickness, although it did not differ from the combined Q + BS and T treatments but did differ from the Q treatment.

Fig. 11 presents data on the number of locules (E2) and physicochemical quality variables from the second harvest: Brix degrees (F2), citric acid percentage (G2), and pH (H2). In the small and medium classifications, no significant differences in the number of locules were found among treatments. However, in the large classification, the combined Q + BS treatment exhibited the highest number of locules, followed by BS and Q, with no significant differences among these treatments but differing from T, which had the fewest locules. The concentration of Brix degrees in the three classifications (small, medium, and large) was highest in the control treatment, showing concentrations of (3.43, 3.66, and 2.33, respectively), with significant differences in the small and medium sizes compared to other treatments. In the medium size classification, the lowest concentration was observed in the BS treatment (1.3), while the large classification showed no significant differences among treatments. Regarding citric acid concentration in the small size classification, the highest concentration was found in the Q treatment, causing differences compared to other treat-ments. However, BS and T showed no differences, while the combined Q + BS treatment had the lowest concentration (0.03%). In the medium classification, the highest citric acid concentration was found in the Q treatment (0.115%), differing from the other treatments, while the BS and T treatments did not show differences, both having concentrations of (0.44%). In the large classification, the control treatment had the highest citric acid concentration (0.05%), showing differences from the other treatments, but Q and the combination with Q + BS showed no significant differences (both at 0.040%), while the lowest concentration was in the BS treatment (0.032%).

Finally, the pH concentration in the small classification was highest in the Q and BS treatments (5.8 and 5.9, respectively), with no significant differences between these treatments. Similarly, Q + BS and T treatments did not show significant differences, with pH values of (5.6 and 5.7, respectively). In the medium classification, the highest pH was found in the combined treatment (5.9), with no significant differences from T and Q treatments (both at 5.8). However, T and Q differed from the BS treatment, which had a pH of (5.7). In the large classification, no significant differences were found among treatments, with pH values ranging from (5.7 to 5.9).

The effects of the biofertilizer treatments on capsaicin content across fruit sizes and harvests are shown in Fig. 12. In Table 7, the two-way ANOVA revealed highly significant effects (p < 0.001) of treatment and fruit size on most morphological and physicochemical quality variables. The combined treatment Q + BS exhibited the highest values of fruit length, diameter, weight, pericarp thickness, and capsaicin content, indicating a positive interaction between both biological factors. The treatment × size interaction was significant only for capsaicin content (p = 0.017), suggesting that the response to the biofertilizers depends on the stage of fruit development.

A comparison of yields was made with various authors based on the results obtained in our procedure. As a result, a yield of 25.19% in UNPA chitosan was achieved, which was compared with the results obtained by [28,29,30], who achieved lower yields of 74%, 61%, and 53.75%, respectively. This was attributed to the increase in temperature during the deproteinization and demineralization process, as well as the higher concentration of the solvents used. On the other hand, Luo et al. [31] presented a yield of 14.5% because the extraction was performed using different solvents and concentrations. These differences can be explained by various factors, including sample processing, treatment type, solvent concentration, sample type, and origin. According to De Castro et al. [32], possible causes for obtaining a low yield in chitosan are that the treated sample contains a higher amount of proteins and fats, as the composition of chitosan from different sources (crustaceans, insects, and fungi) varies not only with the species in question but also with the characteristics used for chitosan extraction, which complicates the comparison among various authors.

Similarly, a comparison was made between commercial SIGMA chitosan and UNPA chitosan based on the percentage of minerals. The result was 2.71% for UNPA chitosan, which, when compared to the results obtained by Luo and Hernández-Cocoletzi et al. [31,33], who obtained chitosan from shrimp exoskeleton with 0.95% and 1.40% minerals, is superior. Conversely, when comparing our results with those presented by Martínez et al. [29], who evaluated chitosan from exoskeletons of estuarine and offshore shrimp, using the same methodology, which resulted in 14% and 24% final minerals; Torres et al. [34] evaluated chitosan from shrimp of different origins, finding 29% and 18% minerals. Additionally, Castro and Vidal found mineral values of 29.0% in chitosan from Emerita analoga exoskeletons [35]. The results obtained for each author differ because Martínez et al. [29] state that the physicochemical properties of chitosan depend on both the sources and the production conditions of the shrimp.

Regarding SIGMA chitosan, its percentage of minerals was 4.5%. Based on these results, it has been confirmed that UNPA chitosan presents superior quality compared to commercial SIGMA chitosan, according to the parameters established by Alvarenga [36]. These same authors mention that for high quality chitosan, the mineral percentage should be less than 1%, and a low ash content can be an essential factor for better solubility. Therefore, the percentage of ash is influenced by the presence of mineral impurities, such as calcium, CaCO₃ salt content, or even by metal contaminants. This result is largely subject to the origin, properties, and extraction conditions [37].

In relation to the percentage of demineralization, values of 97.28% and 95.46% were obtained for UNPA and SIGMA chitosan, respectively. These values contrast with the 98% obtained by [38]. The percentage of protein obtained from the exoskeleton is identical to that reported by Suárez-Villa et al. [38], at 44.7%, and lower than that obtained by Cisneros-Pérez et al. [39], at 42.33%, and higher than that reported by Luna [40] at 36.2%. A final protein percentage of 7.43% was obtained for chitosan, which is higher than the 2.83% achieved by Álvarez et al. [41]. On the other hand, Castro and Vidal [35] achieved a protein percentage of 0.83%, although it exceeds that achieved by Martínez et al. [29] in highseas and estuarine shrimp chitosan, which were 34.67% and 37.16%, respectively. The 7.43% of protein in UNPA chitosan is very low compared to the protein content in the shrimp exoskeleton we obtained, which is 44.62%. This confirms that the deproteinization process was effective, as a significant number of proteins were eliminated.

According to the results obtained for the percentage of deproteinization, percentages of 83.34% and 80.01% were achieved for UNPA and SIGMA chitosan, respectively. Compared to the 87.5% deproteinization obtained by Pérez et al. [37], the values obtained were lower. On the other hand, Castro and Vidal [35] obtained a percentage exceeding 89%. Regarding the moisture percentage, values of 11.5% and 11.3% were obtained for UNPA and SIGMA chitosan, showing a contrast with the 14% obtained by Martínez et al. [29], presenting lower results. In contrast to the results obtained by Torres et al. [34] at 0.96% and Luo et al. [31] at 2.7%, higher results were evidenced, while Velasco and Ibitoye et al. [42,43] analyzed two exoskeleton samples, obtaining values below 6.6% and 6.5%. The value obtained from UNPA chitosan is estimated to be equivalent to that of commercial chitosan at 11%, according to the findings of Barra et al. [44].

With these results, it is confirmed what Martínez et al. [29] reported, indicating higher moisture in shrimp chitosan, evidencing a decrease in the molecular mass of the polymer, given that a greater amount of free hydrophilic groups, such as hydroxyls, amines, and aldehydes, at the end of the polymer chain increases the hygroscopicity of the polymers after deacetylation. These results explain the increase in water absorption capacity as humidity following dehydration. Additionally, it is noted that moisture content is one of the fundamental factors for evaluating chitosan. Molecular water is adsorbed onto the polymer chains, exerting a plasticizing effect on thermal stability and structure. Likewise, the moisture percentage is related to the removal of acetyl groups as a consequence of the thermo alkaline deacetylation of chitin, which generates free amino groups in the polymer chain. This site is sensitive to the formation of hydrogen bonds with the oxygen of free radical –OH [45].

Finally, the pH obtained in the shrimp exoskeleton was 8.8, while that of UNPA chitosan and SIGMA chitosan is 3.4, values higher than 4 and 4.5 in highseas and estuarine shrimp [34]. The pH of chitosan is a quality parameter, as it constitutes a component of its chemical characteristics. It should be noted that there have been no studies conducted to evaluate the level of depigmentation in chitosan using the Hunter Lab method, which prevents comparisons from being made. However, it is maintained that this depigmentation is associated with the oxidative effect of sodium hypochlorite, and that the color obtained from chitosan depends on its commercial value, as established by Romero-Serrano and Pereira [45]. The chitosan that exhibits the highest commercial value is beige, obtained from fungi, followed by white to grayish chitosan, obtained from squid and shrimp. This is related to our results for UNPA chitosan, which has a whitishgray color with color number #BEBAB4. The price of chitosan ranges between $700,000 and $800,000 pesos per ton, depending on its purity [46].

The levels of deacetylation, deproteinization, demineralization, moisture percentage, pH, and color are considered, along with the factors influencing its commercial value [47]. It is required to point out that particle size has an effect on the processes of deproteinization and demineralization [48,49,50]. Chitosan is widely recognized for its ability to induce defense responses, modulate plant growth, and enhance nutrient absorption due to its polyelectrolytic properties and its interaction with receptors on plant cell membranes.

For UNPA chitosan, the O-H and N-H groups were reported at bands of 3440 and 3280 cm⁻¹, respectively, which fall within the range identified by Jeon et al. [51], who recognized these groups at bands of 3300 to 3500 cm⁻¹. Similarly, these findings align with those obtained by González and Silva [28], at 3430 cm⁻¹, while Martínez-Robinson et al. [52] reported a band at 3363 cm⁻¹. The C-N group was found in the band located at 1380 cm⁻¹ according to González and Silva [28], who also detected it at bands of 1581 cm⁻¹. Furthermore, the amide III band was located in the range obtained by González and Silva [28] at 1376 cm⁻¹.

The C=O bonds of amide I were found at bands of 1660 cm⁻¹, which are higher than those reported by González and Silva [28] at 1646 cm⁻¹. They mentioned that the most representative bands of chitosan correspond to the N-H stretching for amide I, which were found at bands of 3289 cm⁻¹, while UNPA chitosan was found at bands of 3260 cm⁻¹. In fact, the broad O-H stretching vibration band was located at 3440 cm⁻¹, lower than the 3471 cm⁻¹ reported by González and Silva [28]. Martínez-Robinson et al. [52] also showed strong typical bands at 1644 cm⁻¹ (amide I, C=O group, stretching vibration). It was confirmed that the N-H stretching bands belonging to amide I were found at 3260 cm⁻¹, while González and Silva [28] found them at 3289 cm⁻¹, reaffirming that these are the most representative bands of chitosan.

Seedling height was greater in the Q + BS treatment, differing significantly from other treatments. Similarly, stem diameter was positively influenced by Q and Q + BS. These results are explained by the synergistic interaction between chitosan and Bacillus subtilis. Chitosan acts as a carbon and nitrogen source, enhancing bacterial growth and metabolic activity in the rhizosphere. In turn, the bacteria improve chitosan degradation and nutrient mobilization, creating a positive feedback that promotes plant development. Bolaños Dircio et al. [53] corroborates that these interactions improve crop growth and yield. Luna [40] showed that Bacillus subtilis increases micronutrient availability, such as iron, which is chelated by siderophores and assimilated by the plant.

The biostimulant effect shown in this study could be attributed to its chemical composition, consisting of linear polymers with structural units of 2-amino-2-deoxy-D-glucopyranose linked by 1,4 glycosidic bonds [54,55,56,57]. Once the plant comes into contact with this biopolymer, a cascade of signaling is triggered, stimulating defense mechanisms mediated by SAR, which, in turn, produces a significant increase in the roots, allowing for greater nutrient absorption, thereby strengthening and enhancing their vigor [58].

When comparing height with the results obtained by Izquierdo-Oviedo et al. [59], who reported a height of 3.99 cm after applying chitosan to the roots and foliage of Capsicum annuum L. (bell pepper) at a dose of 1 mg L−1 15 days posttransplant, our results were superior at 6.88 cm. These findings align with previous reports indicating that chitosan can enhance plant growth and physiological responses due to its antioxidant properties and biostimulant activity [60]. Similarly, a comparison was made with the results obtained by Martínez et al. [29], who evaluated Lombrihumus + Kuntan + Lidavital^® (Lidavital^®, Ciudad de México, México) at 35 days after planting, obtaining a height of 8.75 cm, while our results reached 10.82 cm. Furthermore, when comparing our results with those reported by Álvarez-Pinedo et al. [61], at 45 DAT, who evaluated the application of chitosan combined with ecomic in Capsicum annuum L. (bell pepper), it was found that their results showed significant differences compared to the control at 13.75 ± 0.24 cm; however, our results remained superior at 17.54 ± 0.46 cm.

Regarding stem diameter, it can be inferred that the most suitable treatment at 15 and 45 DDS was Q + BS. According to Bolaños Dircio et al. [53], who evaluated Bacillus licheniformis M2-7 in poblano pepper (Capsicum annuum L.) at 15 DDS, a diameter of 0.102 cm was presented, while our results were superior with Q + BS at 0.149 ± 0.004 cm. Additionally, a comparison was made with the results obtained by Martínez et al. [29], who applied Lombrihumus + Lidavital^® foliarly, evaluating at 35 DDS, which showed a diameter of 0.285 cm, whereas our result was lower at 0.155 cm. Similarly, the results obtained resemble those shown by Lara-Capistrán et al. [13], who evaluated the application of Bacillus subtilis in sweet pepper (Capsicum annuum L.) and determined that it differed positively from the control.

Lastly, no significant differences in the number of leaves were observed among the treatments, consistent with the findings reported by other authors [62], who evaluated QuiToMax^® (Cuban Institute of Chitin and Chitosan Derivatives, Cuba) at a concentration of 1 g L⁻¹ on tomato seedlings of two varieties, ESEN and L-43, in which no significant differences were observed for this variable. On the other hand, Ref. [63] evaluated the effect of HMA-chitosan and HMA-bionutrient in foliar application at 35 days after emergence and showed a result exceeding 6 leaves compared to our results.

Overall, results have been reported indicating the potential of Bacillus, particularly CBRF12 and CBCC57, for the growth of habanero pepper seedlings [64], Other Bacillus spp. strains have reported increases in the growth of Capsicum annuum L. seedlings [65].

The results show an effect of the application of biofertilizers on plant growth involving height, leaf area, and number of shoots. Among the four treatments evaluated, Bacillus subtilis significantly promoted plant height at 60 DAT (120.19 cm) (p ≤ 0.05). Regarding the stem diameter variable, no differences were found among the treatments. The Capsicum annuum L. plants treated with Q, BS, and their combination experienced a significant increase in leaf area and number of shoots at different growth stages, confirming that the use of biofertilizers such as Bacillus subtilis is a viable alternative for plant growth, as it has the ability to produce low molecular weight organic acids that facilitate the release of soluble phosphorus from the soil. This macronutrient is utilized by plants for photosynthesis, energy, carbon source synthesis and degradation, enabling them to achieve greater growth.

Additionally, they can produce indole-3-acetic acid (IAA), a plant hormone that contributes to cell elongation and division, tissue differentiation, and responses to light and gravity [66]. Furthermore, they provide protection to plants against phytopathogens, positively influencing plant growth [67]. Ultimately, bacteria of the Bacillus genus have a high ability to adhere to and colonize plant roots, which can significantly impact their growth [67].

On the other hand, regarding the application of chitosan, various authors have attributed that chitosan stimulates chlorophyll production and photosynthesis in plants, leading to increased development [68]. Chitosan is also used as a biostimulant, not only to promote plant growth but also to induce tolerance to abiotic stress or resistance to pathogens. Among other functions, it stimulates the activity of enzymes essential for nitrogen metabolism and optimizes its transfer in the leaves, enhancing foliar function in growth and development [69]. Chitosan enhances photosynthesis, and its molecular activity varies among cells depending on their specific physiological and chemical characteristics, resulting in increased root development, flowering, and potential crop yield [70,71].

In relation to our research, it has been reported that the application of chitosan in peppers has improved vegetative growth in both height and stem diameter. Similarly, in tomato cultivation, it has been shown to increase growth and the number of flowers [72]. The findings of this study align with those obtained in our research, demonstrating that chitosan has advantages in the morphological characteristics of Capsicum annuum L. plants. Additionally, another study revealed that foliar application of chitosan increased vegetative growth and yield in cucumbers [73].

Furthermore, Gamboa-Angulo et al. [74] evaluated the application of Bacillus subtilis in Capsicum annuum L. at a concentration of 1 × 10¹² CFU mL⁻¹, and their results were superior to those presented in our study. When comparing our results with those reported by Reyes-Pérez et al. [75], who examined the application of chitosan in tomato cultivation, it was found that applying doses of 150, 200, 250, and 300 mg L⁻¹ had a significant impact on height, significantly increasing values (p ≤ 0.05) compared to the control.

Several authors have demonstrated growth stimulation and yield increases through the use of chitosan in various crops, including tobacco, beans, and maize [76,77]. It has been noted that microbial biofertilizers are a feasible option for integrating into production systems to reduce the use of chemical fertilizers, produce food without harming the environment, humans, or animals, and positively affect growth [74]. In the combination of chitosan and a bacterial consortium, which includes Bacillus subtilis, it has been reported that they are an alternative for grain sorghum production [78]. There is also reference to the use of Bacillus subtilis to evaluate the agronomic response of sweet peppers under greenhouse conditions, indicating it as an alternative for production [79].

According to the results obtained, the highest concentration of chlorophyll A was observed in plants treated with chitosan, while for chlorophyll B and A + B, the highest concentration was recorded with the combined Q + BS treatment. These results can be supported by studies showing that chitosan increases the diameter of chloroplasts in the orchid Dendrobium at concentrations of 10–50 ppm in leaves, as well as influencing chloroplast gene expression, potentially leading to higher chlorophyll concentrations [80]. Additionally, it has been confirmed that chitosan enhances photosynthesis (gas exchange and photochemical reactions) in maize crops [81,82]. This effect may be due to the action of chitosan molecules at physiological and molecular levels in leaves, such as enhancing defense mechanisms against reactive oxygen species (ROS) or acting on nuclear and chloroplast genes involved in the photosynthetic process [83].

Other studies indicate that establishing a response model to the effects of chitosan on photosynthetic pigments is complex, as the impact depends on the species and concentration [84]. Additionally, it has been demonstrated that Bacillus subtilis increases photosynthetic capacity, enhancing photosynthetic efficiency and chlorophyll content in Arabidopsis, attributed to an increase in the number of chloroplasts [85]. Similarly, Bacillus subtilis has been reported to delay leaf senescence, keeping leaves green for longer, which may explain the higher chlorophyll concentration in our treatments [86].

A comparison was made with studies that evaluated the application of biofertilizers, including Bacillus subtilis, in Capsicum annuum L. (bell pepper) cultivation, which demonstrated that the highest chlorophyll concentration was observed with Bacillus subtilis inoculation, similar to our results [86]. Likewise, other studies reported that inoculation with the LPM1 strain of Bacillus subtilis in VLK and SVN varieties of bell pepper resulted in higher chlorophyll content compared to the control [87]. Conversely, some studies found no significant differences in the application of Bacillus subtilis in Capsicum annuum L. plants [88].

Studies in white clover and maize crops have also shown an increase in chlorophyll content with the inoculation of Bacillus subtilis [89,90]. Among studies demonstrating increased chlorophyll content in various crops, the application of chitosan nanoparticles in triticale seedlings at a dose of 0.2 mg mL⁻¹ resulted in the highest chlorophyll concentration [91]. Similarly, the application of chitosan oligomers in coffee plants at various concentrations showed that increasing the dose also increased chlorophyll A and B concentrations, with the most effective dose being 60 mg L⁻¹ [92]. However, other studies reported that the application of chitosan oligosaccharides did not alter chlorophyll A and B concentrations in bean plants [93]. Likewise, the use of Quitomax^® (Cuban Institute of Chitin and Chitosan Derivatives, Cuba) in bean plants at a dose of 500 mg L⁻¹ showed no significant differences compared to the control [77].

For the variable electrical conductivity (EC) of the petiole sap, a difference was observed only in treatment T at 105 DAT, which exhibited the highest electrical conductivity. However, it is difficult to make various comparisons with other studies since there are no similar reports; nonetheless, the results may be attributed to the nutrients present in the substrate (vermicompost) being in a saturated or assimilable form, leading to greater absorption by the plants, as well as the content of humic acids in the vermicompost providing greater availability of magnesium (Mg), which interferes with protein transport, thereby improving electrical conductivity that indicates the concentration of dissolved salts [94].

In studies determining the quality and yield of bell pepper fruits (Capsicum annuum L.) in response to the conductivity of the nutrient solution, results showed that an electrical conductivity of 2.0 dS m⁻¹ or higher increased bioactive compounds but decreased fruit yield and firmness [95]. In this research, it was found that at 90 DAT, the treatment with the highest electrical conductivity was T (600 μs cm⁻¹), presenting significant differences compared to the other treatments.

It is also mentioned that electrical conductivity is a good indicator of the amount of fertilizer available to the plant. It has been reported that a greater exploration area for the roots provides adequate electrical conductivity and facilitates nutrient absorption [96]. It is important to note that this will depend on the edaphic conditions of the soil where the crop is established [97].

pH is an important factor that influences the growth of microorganisms. Some bacteria generally grow at low pH (3.0); the optimal pH range for bacteria is from 6.0 to 8.5, and only a few prefer pH 8.5 or higher [98]. This is consistent with the data obtained in the treatments of this research. It was generally observed that the concentration of N–NO₃⁻ in the cellular extract of the petiole (ECP) varied, with the maximum concentration presented in the application with chitosan at 1203 mg L⁻¹. These results are higher than those reported in studies evaluating nutrient solutions [99], and superior compared to our control treatment, which averaged 500.

Additionally, it has been reported that the concentration of nitrates in tomato plants varies according to the amounts of nitrates added to the nutrient solution [100]. These variations are detected during the vegetative stages and at the beginning of flowering. It should be kept in mind that when fertilizers with different elements are applied, compatibility between them must be considered. Moreover, the ionic interactions between different ions in the soil must be taken into account. An increase in potassium in soils reduces the absorption of calcium and magnesium, excess calcium insolubilizes soil phosphates, and excess chloride reduces the absorption of NO₃⁻. Conversely, an increase in nitrogen in soils improves potassium absorption [101].

When analyzing the effect of chitosan on biochemical indicators of fruit quality in tomato cultivation, the concentration of Brix degrees depending on the chitosan dose was found to range from 5.25 to 3.35 [102]. In the present study, the highest value was 105 DAT in the treatment with Q, presenting values higher than (3 ± 0.03), so there are no studies to compare regarding the fruits of Chile Soledad with the application of chitosan and Bacillus subtilis, making comparisons with various authors difficult. In a study on peppers grown in greenhouses, biostimulated plants based on algae showed a higher number of commercial fruits, along with a higher leaf area index compared to the non biostimulated control [103]. This coincides with the results obtained from the application of Q, Q + BS, and T. A significant increase in growth and yield of the crop was reported, with no statistically significant differences between the treatments compared to the control (p ≤ 0.000).

Bacillus subtilis applied to the root and salicylic acid applied to the foliage of A. solani in tomato plants showed similar outcomes [104]. On the other hand, results from another investigation indicated that in the yield of sorghum and protein quality of the grain, the most efficient treatment (with p < 0.05) was with QUI 100% + microbial consortium using 50% regional fertilization [105], consistent with the results obtained.

In a study on strawberry cultivation, the application of chitosan resulted in a significant increase in yield regarding fruit weight, as well as early and total yield. This was accompanied by vegetative growth (plant length, number of leaves/plant, leaf area, root growth, fresh and dry weight) [105].

Furthermore, an analysis on the application of chitosan and a microbial consortium in sorghum showed that in treatments with chitosan, plant length increased as the concentration increased by 50% and 100%, reaching a maximum length of 22.4 cm. Regarding grain yield (kg ha⁻¹), this occurred with 100% chitosan + BAC microbial consortium + 50% fertilization [106]. This indicates that the application of chitosan contributes to length and yield similar to what was obtained in our research.

Regarding fruit weight, it was determined that the application of Q resulted in a higher weight compared to the application of BS, which was somewhat similar to what was obtained in a study where chitosan was the second most important biostimulant for Capsicum annuum L., after humic acid, showing better results in fruit weight, as well as in the number of fruits per plant, fruit length, and fruit diameter [107].

Additionally, it has been found that chitosan exerts a stimulating effect on the following variables: plant height, clusters with fruits, fruits per plant, equatorial and polar diameter of the fruit, pericarp diameter, and fruit mass [108].

Although the results of this study demonstrated positive effects of chitosan and Bacillus subtilis on the growth, yield, and fruit quality of Soledad chili (Capsicum annuum L.) under non-controlled greenhouse conditions, several limitations should be considered. First, the results may vary in interaction with the natural soil microbiota, which could influence treatment efficacy and nutrient availability. Second, the scalability of chitosan and B. subtilis application at a commercial level requires optimization in terms of application volume, frequency, and associated costs. Additionally, practical considerations for farmer adoption, such as compatibility with existing agricultural practices and input availability, have not yet been fully evaluated. As next steps, experiments at a larger commercial greenhouse scale are recommended to validate the effectiveness of the treatments under real production conditions. It is also necessary to investigate interactions with soil microbiota and their impact on growth, yield, and fruit quality. Finally, optimizing the dose, frequency, and application method of chitosan and B. subtilis will help establish practical and economically viable protocols for adoption by growers.