Waitea circinata : a novel biocontrol agent against Meloidogyne enterolobii on tomato plants
Kategoria artykułu: Research Paper
Data publikacji: 14 mar 2025
Otrzymano: 17 cze 2024
DOI: https://doi.org/10.2478/jofnem-2025-0002
Słowa kluczowe
© 2025 Gianlucca de Urzêda Alves et al., published by Sciendo
This work is licensed under the Creative Commons Attribution 4.0 International License.
Tomato (
Several species of
The available strategies for managing
Given the environmental and practical limitations of chemical, genetic, and cultural control strategies, biological control has emerged as a promising alternative for sustainable nematode management. However, the biological control agents currently available in Brazil, including species from the genera
Recent research has explored novel biological agents targeting
This study explores the potential of
We hypothesize that
The study was conducted in two greenhouses at the Universidade Federal de Goiás, Brazil (Escola de Agronomia and Instituto de Ciências Biológicas IV), located at 16°35′46.4″S 49°16′48.4″W and 16°36′04.4″S 49°15′47.6″W. The climate is classified as Aw (tropical savanna) according to Köppen's classification (18). Greenhouse conditions were maintained with average temperatures of 27 ± 4°C and relative humidity at 70%. The study was conducted between December 2022 and July 2023.
An experiment was conducted to assess the ability of
A follow-up experiment assessed the enzymatic activity of the most suppressive treatment identified in the earlier experiments. It used four treatments (Water control, Wc: soil drenching with 15 g.L−1
The
The
The extraction procedure was an adaptation of an established method (19). Infected roots were fragmented into 2 cm pieces and blended with 0.5 % sodium hypochlorite for 30 seconds. The suspension was filtered through 100- and 500-mesh sieves. After washing under running water, the blended root material with water was placed in a 50 mL plastic tube, and 5g of kaolin was added. The filtrate was centrifuged at 1800 rpm for 5 minutes. The pellet was resuspended in 25 mL of a solution of 454 g.L−1 of sucrose, then centrifuged at 1800 rpm for 1 minute and filtered in 500-mesh sieves. The collected filtrate was derived from multiple plant samples. A 50 mL aliquot was obtained from the total filtrate volume to count the number of eggs and second-stage juveniles (J2) of
The immersion method consisted of taking seedlings with the substrate and immersing the root system completely in 70 mL of the mycelial suspension for 24 hours. After immersion, the seedlings were transplanted into 1L pots filled with a 1:1 sterilized soil-sand mixture pre-autoclaved at 120°C for 40 minutes. The drenching method consisted of taking seedlings, directly transplanting them into 1L pots, and applying 50 mL of the mycelial suspension via soil drenching. All seedlings received 5 g of NPK fertilizer (4-30-16) after transplantation. All plants were irrigated twice daily at 6:00 AM and 6:00 PM using an automated sprinkler system. Seven days post-treatment, the plants were inoculated with 2000 eggs + J2 nematodes per plant.
In all experiments, 35 days after nematode inoculation, the following traits were measured: Shoot Length (SL): Measured from the plant base to the last leaf node using a measuring tape; Root Fresh Weight (RFW): Roots were washed, and dried for 1 hour, and weighed using a digital scale (0.001 g precision); Shoot Fresh Weight (SFW): Weighed using the same digital scale; Final Nematode Population (FP): Total nematode count in roots used the method in the extraction for inoculum preparation; Population Density (DENS): Calculated as FP / RFW (eggs + J2 per gram of root) and Reproduction Factor (RF): Calculated as FP / initial inoculum (20).
Leaves were collected on days 0, 3, 5, 7, 8, 14, and 21 after transplantation and stored in a freezer under −20°C until extraction. For protein extraction, leaf tissue was ground in liquid nitrogen, and 0.1 g of grounded leaves were mixed with an extraction buffer (10 mM Tris-HCl, 150 mM NaCl, 2 mM EDTA, pH 7.0). Samples were centrifuged at 13,000 rpm for 30 minutes at 4°C, and the supernatant was used for enzyme assays.
Peroxidase (POX - EC 1.11.1.7) activity was measured by ABTS (2,2′-azino-bis(3-ethylbenzothiazoline-6- sulfonic acid) oxidation using its colorimetric property at 405 nm, and enzyme activity was defined as the activity that catalyzes the formation of ABTS that increases the absorbance of 1 unit per hour (21). Catalase (CAT - EC 1.11.1.6) activity was measured by the reduction of H2O2 at 240 nm, and enzyme activity was computed by calculating the amount of H2O2 decomposed (22). Chitinase (CHI - EC 3.2.1.14) activity was quantified using 3,5-dinitrosalicylic acid (DNS) at 540 nm using colloidal chitin and one unit of enzyme activity was defined as the amount of enzyme needed to release 1 mmol reducing sugar min−1 using a modified method (23). All essay readings were conducted on a spectrophotometer.
A Two-way ANOVA and Tukey post-hoc test (p < 0.05) were used to analyze the factorial interaction between application methods and dosages and separate the group means. Each variable residual was analyzed for normality and homoscedasticity using Shapiro-Wilk and Levene's tests. Box-Cox transformations (24) were applied in variables to meet Two-way ANOVA assumptions. A regression analysis was applied for nematological traits (FP, DENS, RF) to find the best descriptor of the dosage effect on nematode reproduction.
Student's t-test was used to compare nematological traits between treatments Me and Me+Wc in the enzymatic experiment to address the consistency of the nematode suppression. Enzyme activities were analyzed using ANOVA with Tukey's test (p < 0.05) to distinguish the group means. Principal component analysis assessed the relationship between enzymatic activity and nematode reproduction factor (RF). All analyses were performed using RStudio (version 4.3.1).
After performing the Two-way ANOVA, no interaction was detected between the application method and dosage for any traits; therefore, the two factors were analyzed independently (Table 1). An analysis of the application method revealed that the drenching method resulted in significantly higher SL, RFW, and SFW than the immersion method (p < 0.05, Table 2). The immersion method exhibited higher population density (DENS), indicating a less developed root system for a similar number of nematodes compared to drenching (p < 0.05, Table 2). The immersion method might affect plant development, reducing vegetative traits compared to the drenching method.
Analysis of variance of vegetative and nematological variables from tomato plants inoculated with
Application Method | SL(cm) | 1.303b | 0.257b | 100.229b | <0.001*b |
RFW(g) | 80.338 | <0.001* | 13.596 | <0.001* | |
SFW(g) | 0.929 | 0.338 | 53.982b | <0.001*b | |
FP | 1.641b | 0.208b | 0.939 | 0.337 | |
DENS | 2.468b | 0.120b | 14.793 | <0.001* | |
RF | 1.575b | 0.213b | 0.925 | 0.341 | |
Dosage | SL(cm) | 0.636b | 0.673b | 1.110b | 0.368b |
RFW(g) | 0.151 | 0.979 | 0.166 | 0.973 | |
SFW(g) | 0.906 | 0.482 | 1.318b | 0.272b | |
FP | 2.957b | 0.017*b | 2.263 | 0.063 | |
DENS | 2.355b | 0.048*b | 2.278 | 0.061 | |
RF | 2.949b | 0.017*b | 2.256 | 0.063 | |
Method * Dosage | SL(cm) | 0.524b | 0.757b | 0.712b | 0.617b |
RFW(g) | 1.317 | 0.266 | 1.930 | 0.106 | |
SFW(g) | 0.960 | 0.448 | 1.962b | 0.101b | |
FP | 0.507b | 0.769b | 0.138 | 0.982 | |
DENS | 0.506b | 0.770b | 1.078 | 0.384 | |
RF | 0.510b | 0.767b | 0.081 | 0.994 |
SL = Shoot length; RFW = Fresh root weight; SFW = Fresh weight of aerial part; FP = final population; DENS = Density of nematodes per gram of root; RF = Reproduction Factor.
The characters have statistically significant differences (p<0.05).
Analysis carried out with data transformed by the Box-Cox transformations (Box and Cox, 1964).
Vegetative and nematological traits of the first and second trials, differentiating the application methods of
1 | SL (cm) | 110.69 ± 20.01 | 107.58 ± 15.81 | 1.303b | 0.257b |
RFW(g) | 40.46 ± 10.18 | 62.26 ± 11.69 | 80.338 | <0.001* | |
SFW(g) | 119.18 ± 29.73 | 113.88 ± 19.29 | 0.929 | 0.338 | |
FP | 1150.00 ± 677.24 | 1555.95 ± 1187.65 | 1.641b | 0.208b | |
DENS | 29.90 ± 18.09 | 26.06 ± 20.63 | 2.468b | 0.120b | |
RF | 0.57 ± 0.33 | 0.78 ± 0.59 | 1.575b | 0.213b | |
2 | SL (cm) | 90.33 ± 10.11 | 114.10 ± 8.03 | 100.229b | <0.001*b |
RFW(g) | 29.81 ± 5.77 | 35.90 ± 7.01 | 13.596 | <0.001* | |
SFW(g) | 91.96 ± 15.30 | 112.38 ± 7.67 | 53.982b | <0.001*b | |
FP | 2443.33 ± 797.27 | 2248.33 ± 784.38 | 0.939 | 0.337 | |
DENS | 92.00 ± 26.37 | 68.72 ± 23.05 | 14.793 | <0.001* | |
RF | 1.22 ± 0.39 | 1.12 ± 0.39 | 0.925 | 0.341 |
SL = Shoot length; RFW = Fresh root weight; SFW = Fresh weight of aerial part; FP = final population; DENS = Density of nematodes per gram of root; RF = Reproduction Factor.
The characters have statistically significant differences between them (p<0.05).
Analysis carried out with data transformed by the Box-Cox transformations (Box and Cox, 1964).
The analysis of the dosage factor revealed that the FP, DENS, and FR variables had at least one dosage group mean different from the others (p < 0.05, Table 3). The Tukey test revealed that for FP, DENS, and FR variables, the 15 g.L−1 dosage was significantly different from the control treatment (p < 0.05). The regression analysis of the effect of doses of
Vegetative and nematological traits of the first and second trials, differentiating the dosages of
1 | SL (cm) | 101.83 ± 18.79 | 111.92 ± 13.39 | 108.64 ± 23.33 | 107.28 ± 23.35 | 112.92 ± 13.74 | 112.21 ± 12.68 | 0.636b | 0.673b |
RFW(g) | 51.31 ± 11.16 | 52.84 ± 19.36 | 52.20 ± 17.35 | 50.33 ± 15.50 | 51.74 ± 18.11 | 49.75 ± 12.22 | 0.151 | 0.979 | |
SFW(g) | 110.77 ± 21.76 | 120.73 ± 20.72 | 117.25 ± 28.44 | 110.7 ± 27.76 | 112.88 ± 28.34 | 126.85 ± 22.37 | 0.906 | 0.482 | |
FP | 2053.57 ± 1095.70 | 1285.71 ± 644.33 | 1171.42 ± 920.43 | 967.85 ± 825.46* | 1203.57 ± 974.37 | 1435.71 ± 1143.44 | 2.957b | 0.017*b | |
DENS | 40.59 ± 18.41 | 27.95 ± 17.37 | 26.13 ± 24.51 | 20.23 ± 16.33* | 24.38 ± 16.26 | 28.60 ± 19.31 | 2.355b | 0.048*b | |
RF | 1.02 ± 0.54 | 0.64 ± 0.32 | 0.58 ± 0.45 | 0.48 ± 0.41* | 0.60 ± 0.48 | 0.72 ± 0.57 | 2.949b | 0.017*b | |
2 | SL (cm) | 107.8 ± 16.51 | 103.7 ± 16.49 | 101.1 ± 12.63 | 100.9 ± 17.16 | 100.3 ± 14.63 | 99.5 ± 14.59 | 1.110b | 0.368b |
RFW(g) | 34.27 ± 8.99 | 33.00 ± 8.04 | 32.45 ± 5.94 | 32.10 ± 9.38 | 32.11 ± 5.20 | 33.18 ± 5.27 | 0.166 | 0.973 | |
SFW(g) | 95.77 ± 25.13 | 107.27 ± 11.45 | 106.48 ± 13.21 | 103.02 ± 17.83 | 101.01 ± 12.38 | 99.48 ± 11.20 | 1.318b | 0.272b | |
FP | 3015 ± 768.49 | 2470 ± 825.69 | 2120 ± 753.58 | 2035 ± 502.24 | 2070 ± 727.70 | 2365 ± 841.64 | 2.263 | 0.063 | |
DENS | 95.046 ± 33.35 | 79.43 ± 21.00 | 72.47 ± 23.12 | 67.94 ± 27.72 | 74.39 ± 27.22 | 92.90 ± 23.75 | 2.278 | 0.061 | |
RF | 1.50 ± 0.38 | 1.23 ± 0.41 | 1.06 ± 0.37 | 1.02 ± 0.25 | 1.03 ± 0.36 | 1.18 ± 0.42 | 2.256 | 0.063 |
SL = Shoot length; RFW = Fresh root weight; SFW = Fresh weight of aerial part; FP = final population; DENS = Density of nematodes per gram of root; RF = Reproduction Factor.
The characters have statistically significant differences (p<0.05).
Analysis carried out with data transformed by the Box-Cox transformations (Box and Cox, 1964).
ANOVA summary for the polynomial regression model of nematological traits based on mycorrhiza mycelial suspension.
Experiment 1 | PFa | Regression | 2 | 666569.667 | 333284.834 | 25.486 | *0.013 |
Residual | 3 | 39231.778 | 13077.259 | ||||
Total | 5 | 705801.446 | |||||
DENSa | Regression | 2 | 223.023 | 111.512 | 26.146 | *0.013 | |
Residual | 3 | 12.795 | 4.265 | ||||
Total | 5 | 235.818 | |||||
FRa | Regression | 2 | 0.166 | 0.083 | 25.698 | *0.013 | |
Residual | 3 | 0.010 | 0.003 | ||||
Total | 5 | 0.176 | |||||
Experiment 2 | PFa | Regression | 2 | 686036.548 | 343018.274 | 833.725 | *0.000 |
Residual | 3 | 1234.286 | 411.429 | ||||
Total | 5 | 687270.833 | |||||
DENSa | Regression | 2 | 611.832 | 305.916 | 65.484 | *0.003 | |
Residual | 3 | 14.015 | 4.672 | ||||
Total | 5 | 625.847 | |||||
FRa | Regression | 2 | 0.171 | 0.085 | 888.455 | *0.000 | |
Residual | 3 | 0.000 | 0.000 | ||||
Total | 5 | 0.171 |
FP = final population; DENS = Density of nematodes per gram of root; RF = Reproduction Factor.
The characters have statistically significant differences between them (p<0.05).
The regression results indicated that for FP in the first experiment, the 15 g.L−1 dosage reduced the control group's average of 2053.57 nematodes to 960.71 nematodes (53.21% reduction). The second experiment lowered from 3015 to 2035 nematodes (32.50% reduction). Based on the equations of the fitted curves, the lowest possible FP obtained by the curves were the dosages of 14.36 and 15.65 g.L−1, for experiments 1 and 2, respectively (Fig. 1A).

Regression analysis for the effect of
For DENS, in the first experiment, the 15 g.L−1 dosage reduced nematode densities from the control's average of 40.59 nematodes per gram of root to 20.53 (49.42% reduction). In the second experiment, the control average densities were 95.04, lowering to 67.94 nematodes per gram of root (28.51% reduction). The lowest possible DENS values were estimated at 14.71 and 13.03 g.L−1 (Fig. 1B).
For RF, when compared to the controls (1.03), the 15 g.L−1 dosage obtained the highest suppression, reaching 0.48 in the first experiment (53,39% reduction). In the second experiment, it was 1.51 for the controls, reducing to 1.02 (32.40% reduction). The lowest possible values for RF were estimated to be the dosages of 14.21 and 16 g.L−1 on experiments 1 and 2, respectively (Fig. 1C). For the first two experiments, the results converged on the use of dosage of 15g.L−1 and the method of soil drenching for enzymatic analysis due to better performance of vegetative traits averages (Tables 2, 3) and
The second part of the study focused on analyzing the enzyme activity in plants due to
Effect of
Me | 2944 ± 668,3 | 114,43 ± 28,67 | 1,47 ± 0,33 |
Wc+Me | 850 ± 463,1* | 26,38 ± 15,41* | 0,425 ± 0,23* |
FP = final population; DENS = Density of nematodes per gram of root; RF = Reproduction Factor.
The characters have statistically significant differences (p<0.05).
During the first seven days, significant differences in enzyme activity were observed primarily on the third and fifth days after mycorrhizal application. Three days post soil drenching; there was a marked increase in enzyme activity in the treatment with

Peroxidase, catalase, and chitinase enzymatic activity analyzed in tomato leaves under different treatments. Negative control (Water), inoculation of 2000 J2 +
However, CHI activity decreased by 20.02% in the treatment with
All enzyme activities showed differences among treatments on the fourteenth day after mycorrhizal soil drenching and seven days after nematode inoculation. The highest POX activity was recorded for the
On the 21st day after mycorrhizal soil drenching and fourteen days after nematode inoculation POX and CAT activities showed differences among treatments. The Wc+Me treatment (2.76 U.mg−1) presented a higher POX activity than the control (1.21 U.mg−1). The Me treatment (1.13 U.mg−1) showed no differences in POX activity compared to the control group. The Wc treatment (0.42 U.mg−1) presented the lowest POX activity of all treatments (Fig. 2A). For CAT activity, the treatment Wc+Me (76.91 U.mg−1) showed the highest activity. The Me treatment (49.99 U.mg−1) and Wc (61.52 U.mg−1) treatments were not different from the control treatment (46.14 U.mg−1) (Fig. 2B). No differences were observed for the CHI activity on the 21st day (Fig. 2C).
The first two principal components, Dim 1 and Dim 2 explain 38.92% and 23.34% of the total variance, respectively. Confidence ellipses at a 95% confidence level were added to visualize the spread and overlap of the treatment categories. The PCA plot (Fig. 3) displays the samples categorized by treatments, with confidence ellipses around each category. The black (Me) and red (Water) ellipses overlap slightly, indicating some similarity in the samples. Both treatments are clustered near the origin, suggesting limited variance in these treatments along the principal components. Meanwhile, both the green (Wc) and blue (Wc+Me) ellipses are distinctly separated from the other treatments, indicating a unique pattern in those categories (Fig. 3).

PCA plot of individuals with confidence ellipses around treatment categories and their contribution to Principal Components.
The PCA was performed on standardized data to ensure that each variable contributes equally to the analysis. The covariance matrix was computed, and eigenvalues and eigenvectors were extracted to determine the principal components. The vectors’ direction and length indicate each variable's contribution to the principal components. Longer vectors represent variables with higher contributions. The color coding represents the squared cosine (cos2) values, indicating the quality of the representation of the variables on the factor map. Higher cos2 values (red) indicate better representation. CAT.0D, CAT.7D, and POX.21D had longer vectors, indicating a higher contribution to Dim 1 and Dim 2. These variables were key contributors to the variance captured by the first two principal components. Variables related to similar time points or conditions tend to cluster together, suggesting underlying patterns. POX.3D and CAT.3D formed a cluster. Meanwhile, CHI.3D showed an opposite behavior, indicating a strong negative correlation. CAT.5D, CHI.5D, and POX.5D formed a cluster, indicating similar behavior, indicating a strong positive correlation. Both days 7 and 8 showed shorter vectors, except CAT.7D, indicating a lower quality of representation. No clusters were observed on day 14. On day 21, CHI, CAT, and POX showed approximately similar behavior, indicating similar variance characteristics and positive correlation (Fig. 4).

Principal Component Analysis (PCA) plot of enzymatic activities and collection day and their contribution to Principal Components.
In the present study, we tested two different methods of
Varying practices in tomato production drove the two-method approach. While soil drenching may suit most field applications, immersion could benefit seedling treatments in pre-planting conditions for fields with known nematode infestations (25). This exploratory investigation into immersion highlights the need for further research to understand its protective effects under different exposure durations. Both methods offer practical value, though cost, labor, and specific crop conditions must be evaluated to optimize applications across different production settings.
The 15 g·L−1 mycelial suspension dose was identified as the optimal concentration for suppressing nematode populations, consistently reducing the FP, RF, and DENS across three experiments. An average 50% reduction in FP, RF, and DENS was achieved, suggesting a consistent suppression effect at this dose. Despite conducting trials across two greenhouses at different times of the year, the nematode suppression results remained consistent, supported by regression analysis. This consistency reinforces the efficacy of the 15 g·L−1 dosage, although observed variations in suppression intensity may stem from environmental differences inherent in the distinct environmental conditions (26).
Previous studies used
The dose-response relationship between biocontrol agents and pathogens can vary depending on the genetic makeup of the plant, pathogen, and biocontrol agent itself, suggesting that genetic factors may influence the efficacy of mycorrhizal applications. This variability is crucial when developing commercial biocontrol products, as different plant-pathogen systems may require customized dosages for optimal effectiveness (27). The findings from this study provide a foundation for future research aimed at refining dosage guidelines for
The study also explored the enzymatic response associated with
In the combined treatment (Wc+Me), POX activity spiked by day 21, aligning with 14 days post-nematode inoculation and marking the midpoint of the ideal nematode life cycle (35,36). This response underscores the potential for
Additionally, CHI activity appeared suppressed during initial colonization. This suppression aligns with reports of mycorrhizal fungi, such as
The results parallel findings on other fungal biocontrol agents, such as
The study highlights W. circinata's potential as a biocontrol agent against root-knot nematodes, paving the way for further research on their application and resistance induction mechanisms. By elucidating these dynamics, we aim to enhance sustainable agricultural practices that reduce reliance on chemical treatments and improve plant resilience through symbiotic fungal associations.
To fully understand
Our study demonstrated that soil drenching with
The suppression rates achieved in this study are comparable to other biocontrol studies involving arbuscular mycorrhizal fungi (AMFs), highlighting