Shapiro–Wilk tests indicated that the distributions of all measured variables did not deviate significantly from normality in either season (p > 0.05). Analysis of variance revealed that the interaction between compost tea level (Ct) and foliar biostimulant treatment (F) significantly affected vegetative growth traits, umbel number, and oil-related traits (oil percentage, oil yield per plant, and oil yield per hectare) in both growing seasons.
Yield response relative to the control treatment
In the first season (Figure 2), oil yield per hectare (Oy ha⁻1) increased progressively with increasing compost tea level. All Ct3 (600 L ha⁻1) treatments produced higher oil yields than the control treatment (T1: Ct1 × F1, 360 L ha⁻1 without foliar spray). The highest oil yield (55.58 kg ha⁻1) was obtained from Ct3 combined with Spirulina at 2 cm3 L⁻1 (T13), which significantly exceeded the control and all Ct1 treatments, surpassing the LSD value (2.93 kg ha⁻1). Relative to the control, this treatment showed a marked increase in oil yield per hectare (Figure 3).

Figure 3: Biplot principal component analysis for the fifteen interacting treatments, based on the four vegetative traits, number of umbels per plant, and the three oil-related traits, using first and second components; a treatment is represented by each dot. Location inside the plot shows the degree of similarity or difference between the treatments according to the input factors (8 growth and oil characteristics); the closer together, the more similar, and the farther apart, the more distinct. PH: plant height (cm), NB: number of branches/plant, HFW: total herb fresh weight (g), HDW: total herb dry weight (g), O%: Essential oil content, Oy/p: oil yield per plant (cm3), Oy/ha: oil yield per hectare (L) and NU: Umbels number/plant T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4), and T15 (Ct3 × F5). Please click here to view a larger version of this figure.
Similar trends were observed for oil yield per plant and number of umbels, indicating that treatments producing higher oil yield per hectare were also associated with greater reproductive output and biomass accumulation. These relationships are presented as statistical associations among measured traits without inferring causal mechanisms.
In the second season, Ct3 × F3 (T13) again produced the highest oil yield per hectare (58.85 kg ha⁻1), significantly exceeding the control and most Ct2 treatments (LSD = 3.17 kg ha⁻1). The ranking of treatments relative to the control was largely consistent across seasons, indicating reproducible treatment effects under the tested environmental conditions.
Principal component analysis (PCA)
Principal component analysis was applied as an exploratory tool to summarize multivariate relationships among the eight measured traits across the 15 treatments. The first two principal components explained 98.4% of the total variance, with PC₁ accounting for 96.5% and PC₂ for 1.9% (Table 2; Figure 4). The dominance of PC₁ reflects strong covariance among growth- and yield-related traits in the dataset rather than distinct underlying biological processes.
| Parameter | PC1 | PC2 | PC3 | PC4 |
| NU | 0.3516 | 0.3064 | -0.3732 | 0.772 |
| O% | 0.349 | -0.6068 | -0.1873 | -0.0634 |
| Oy/p | 0.3591 | -0.11 | -0.0028 | -0.0974 |
| Oy/ha | 0.3583 | -0.1579 | -0.1046 | -0.2021 |
| PH | 0.3476 | 0.5757 | 0.3834 | -0.3296 |
| NB | 0.3495 | -0.26 | 0.7504 | 0.3693 |
| HFW | 0.3551 | 0.3138 | -0.1821 | -0.1535 |
| HDW | 0.358 | -0.0575 | -0.2675 | -0.2847 |
| Eigenvalue | 7.7201 | 0.1518 | 0.0721 | 0.0346 |
| % of Variance | 96.5 | 1.9 | 0.9 | 0.43 |
| Cumulative (%) | 96.5 | 98.4 | 99.3 | 99.73 |
| PH: plant height (cm), NB: number of branches/plant, HFW: total herb fresh weight (g), HDW: total herb dry weight (g), O%: Essential oil content, Oy/p: oil yield per plant (cm³), Oy/ha: oil yield per hectare (L) and NU: Umbels number/plant |
Table 2: Eigenvalues, eigenvectors, percentage variance, and cumulative variance of the first four principal components based on eight growth and oil-related traits across 15 treatment combinations. Trait abbreviations are as follows: PH (plant height, cm), NB (number of branches per plant), HFW (herb fresh weight, g), HDW (herb dry weight, g), O% (essential oil content), Oy/p (oil yield per plant, cm3), Oy/ha (oil yield per hectare, L), and NU (number of umbels per plant).

Figure 4: Loading plot of Eigen vector for 8 studied traits treated with the 15 interacted treatments: PC₁ bars use the left axis; PC₂, PC₃, PC₄ use distinct point shapes and colors on the right axis. PH: plant height (cm), NB: number of branches/plant, HFW: total herb fresh weight (g), HDW: total herb dry weight (g), O%: Essential oil content, Oy/p: oil yield per plant (cm3), Oy/ha: oil yield per hectare (L) and NU: Umbels number/plant. Error bars represent standard errors of the mean (SEM). Please click here to view a larger version of this figure.
All traits showed strong and similar positive loadings on PC₁, including number of umbels, oil percentage, oil yield per plant, oil yield per hectare, plant height, number of branches, herb fresh weight, and herb dry weight (Table 2; Figure 5). Accordingly, PC₁ can be interpreted as a general productivity gradient summarizing variation in biomass and oil yield traits across treatments. PC₂ explained only a small proportion of the variance and contributed limited additional discrimination among treatments.

Figure 5: PCA score plot illustrating the distribution of 15 treatment combinations based on the measured traits. Treatments are defined as follows: T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4), and T15 (Ct3 × F5). Error bars represent standard errors of the mean (SEM). Please click here to view a larger version of this figure.
PCA scores and comparison with the control
Treatments with high positive PC₁ scores (notably T13, T12, T15, and T14) also exhibited the greatest percentage increases relative to the control (T1) across most traits (Figure 5 and Figure 6). For example, T13 showed increases of 94.44% in the number of umbels and 98.18% in oil yield per hectare compared with the control.

Figure 6: Relative percentage increases of the measured traits compared to the control treatment (T1: Ct1 × F1). Treatments are defined as follows: T1 (Ct1 × F1), T2 (Ct1 × F2), T3 (Ct1 × F3), T4 (Ct1 × F4), T5 (Ct1 × F5), T6 (Ct2 × F1), T7 (Ct2 × F2), T8 (Ct2 × F3), T9 (Ct2 × F4), T10 (Ct2 × F5), T11 (Ct3 × F1), T12 (Ct3 × F2), T13 (Ct3 × F3), T14 (Ct3 × F4), and T15 (Ct3 × F5). Please click here to view a larger version of this figure.
Treatments with intermediate PC₁ scores (e.g., T8 and T10) showed moderate improvements over the control, whereas treatments with negative PC₁ scores, including the control and other Ct1 combinations, exhibited minimal relative changes. These patterns indicate that PC₁ scores were closely aligned with baseline-adjusted treatment responses, reinforcing the univariate ANOVA results. However, PCA was used strictly for data reduction and visualization of multivariate patterns and does not replace inferential statistical testing.
Trait associations in the PCA biplot
The PCA biplot (Figure 4) illustrates close associations among oil yield per plant, oil yield per hectare, and herb dry weight, reflecting strong correlations among these variables. Plant height, herb fresh weight, and number of umbels also displayed similar vector orientations, indicating coordinated variation among these traits within the dataset. The number of branches showed partial separation from other yield components, suggesting weaker associations under the present experimental conditions. These relationships describe observed correlations only and do not imply mechanistic linkages.
Treatments (T12–T15) receiving compost tea at 600 L ha⁻1 combined with foliar biostimulants showed superior multivariate performance across growth and oil traits. In contrast, treatments (T7–T11) exhibited comparatively lower values for most evaluated parameters. Treatments positioned in the opposite direction of most vectors, such as Ct1 (360 L ha⁻1) without foliar spray or combined with Nostoc, were associated with lower values for most measured traits.
Overall, PCA provided a concise multivariate summary that complements ANOVA by illustrating how treatments differed from the control across multiple traits simultaneously. The consistency between PCA patterns and baseline-adjusted univariate results supports the robustness of treatment effects without extending interpretation beyond the measured variables.
DATA AVAILABILITY:
The raw data for this manuscript are available within the supplementary files.
Supplementary Table 1: Raw experimental data for the results. Please click here to download this file.