Soil collection

The soil samples were obtained from a relatively undisturbed field in Bukidnon, Central Mindanao, Philippines (coordinates: 7.8324 N, 124.8829 E). It is located at an altitude of 790 m above sea level. The soil belongs to the La Castellana clay series in the order of Alfisols in the Hapludalf’s great group. The climate belongs to type III according to the Corona climate classification. Soils were collected from the top soil layer to a depth of 10 cm, and brought to the Molecular Ecology and Physiology (MEP) Laboratory of the Premier Research Institute of Science and Mathematics (PRISM), Mindanao State University - Iligan Institute of Technology (MSU-IIT), Iligan, Philippines. Soil properties such as pH, OM, N, Phosphorus (P) and Potassium (K) were measured at the beginning of soil collection and at the end of the experiment (45th day). At the time of soil collection soil pH was determined in a 1:1 soil soil-water suspension. The OM content was measured using the Walkley-Black titration method, which is based on the oxidation of OM by a potassium dichromate-sulphuric acid mixture. To extract the available P in the soil, the Olsen method was used, while K was measured using the cold sulphuric acid method (Labajo & Pabiona, 2022Labajo JRN, Pabiona MG, 2022. Physical and chemical properties of soil on selected sugarcane farms in Mt. Nebo, Valencia City, Bukidnon, Philippines. Asian J Agric, 6(2). 10.13057/asianjagric/g060204).

Experimental setup

The experiment was conducted in a screen house. A total of 20 plastic pots were prepared, consisting of four concentrations with five replicates for each treatment. Each plastic pot was filled with 5 kg of soil. Based on the chemical soil analysis of the Department of Agriculture of the CARAGA region, Philippines, the recommended amount of urea fertiliser in the field was 78.26 kg of urea per ha. Considering that 1 ha is equivalent to 2 x 10⁶ kg of soil, this recommended amount is approximately 196.65 mg of urea per 5 kg of soil or 39.33 mg/kg of soil. For the trial, urea with an NPK ratio of 46-0-0 was added to the soil pots at the following concentrations: 19.67 mg/kg soil (labelled as T1), which is 50% below the recommended value; 39.33 mg/kg soil (T2), which is the recommended value in this area based on the soil quality test; 78.66 mg/kg soil (T3), which is twice the recommended value and 0 mg/kg soil as a control. These amounts were applied according to the farmers’ usual schedule in two separate applications with the first 50% of the urea applied directly in solid form and mixed evenly into the soil, while the remaining urea was dissolved in water and sprayed evenly onto the soil 25 days after the first application. The solution was prepared by dissolving urea in tap water, resulting in mass concentrations of 4.66%, 8.9%, 17.8% and 0% for T1, T2, T3 and the control, respectively.

Nematode collection and processing

Immediately after soil sampling, 100 g of soil were taken to represent the original nematode communities and designated as day 0. While the rest of the homogenised soil was distributed among the pots, approximately 100 g of fresh soil were taken from each pot on day 15, 30 and 45. Following Martinez et al. (2019Martinez JG, Quiobe SP, Moens T, 2019. Effects of mercury (Hg) on soil nematodes: A microcosm approach. Arch Environ Con Tox, 77, pp. 421-431. 10.1007/s00244-019-00652-7), nematodes were extracted using a modified tray method and picked using a copper needle under a stereomicroscope (Olympus, BX 41). Nematodes were mounted on permanent slides (Siddiqi, 1997Siddiqi MR, 1997. Techniques and methodologies for nematode disease diagnosis and nematode identification. FAO Plant P, 144, pp. 21-44.) for identification to genus level based on morphological characteristics (Andrássy, 2005Andrássy I, 2005. Free-living Nematodes of Hungary, I (Nematoda errantia).: Vol. I. Hungarian Natural History Museum, Budapest., 2007Andrássy I, 2007. Free-living Nematodes of Hungary, II (Nematoda errantia).: Vol. II. Hungarian Natural History Museum, Budapest., 2009Andrássy I, 2009. Free-living Nematodes of Hungary, III (Nematoda errantia): Vol. III. Hungarian Natural History Museum, Budapest.) under a binocular Olympus CX 22® microscope (400–1000x magnification). The nematodes were assigned as the “oloniser-persister” (cp) scale (Ettema & Bongers, 1993Ettema CH, Bongers T, 1993. Characterization of nematode colonization and succession in disturbed soil using the Maturity Index. Biol Fert Soils, 16, pp. 79-85. 10.1007/BF00369407), and divided into trophic groups, namely bacterivores, fungivores, omnivores-predators, and plant parasites (Ferris & Bongers, 2009Ferris H, Bongers T, 2009. Indices developed specifically for analysis of nematode assemblages. Nematodes as environmental indicators, pp. 124-145. 10.1079/9781845933852.0124; Sieriebriennikov et al., 2014Sieriebriennikov B, Ferris H, de Goede RG, 2014. NINJA: An automated calculation system for nematode-based biological monitoring. Eur J Soil Biol, 61, pp. 90-93. 10.1016/j.ejsobi.2014.02.004).

Statistical analysis

For the urea-free experiment, an independent samples t-test was used to measure the difference in soil physico-chemical properties and community descriptors (i.e., abundance, no. of genera, Shannon-Weiner Index (H’), Simpson’s Evenness Index and Maturity Index (MI)) between day 0 and day 45. Non-metric multidimensional scaling (nMDS) based on the log (x+1) transformation was also used to visualise trends in nematode community composition between day 0 and day 45. Moreover, significant differences in nematode genus composition between day 0 and day 45 was determined using multivariate permutation analysis of variance (PERMANOVA) in PRIMER+ with 999 permutations at p<0.05. SIMPER (Similarity Percentages) analysis was used to calculate the contribution of each species (%) to the dissimilarity between the two groups. To determine the effects of urea on nematode communities, a two-way ANOVA was used to compare the nematode community descriptors, between treatments over time using the Statistica ® software package version 7.0. Prior to ANOVA, data were subjected to meet the normality assumptions. p-values < 0.05 were considered significant.

Suitability of the microcosm approach

In the urea-free experiment, the physicochemical properties of the soil were compared between the beginning (day 0) and the end (day 45) of the experiment. The results showed that the pH decreased significantly (p<0.05), while K increased significantly after 45 days (p<0.5). Other soil properties, such as OM, P and N, showed no significant changes. Nematode density (individuals/100g soil) and number of genera decreased significantly (both p<0.05) by 28% and 35%, respectively, at the end of the experiment. As shown in the nMDS, there was no overlap in nematode genera between day 0 and day 45, indicating distinct nematode assemblages at the start and end of the experiment (Fig. 1). In the same period of time, there was a significant decline in the H’ and the Simpson’s Eveness Index (both p<0.05), as well as the MI (Table 1). Among trophic groups, herbivores (32.4%) dominated the nematode assemblages on day 0, followed by bacterivores (25.2%), fungivores (11.9%), predators (16.3%) and omnivores (14.1%). A shift in the nematode communities was then observed on day 45, with the proportion of herbivores (7.3%), bacterivores (63.9%), fungivores (11.4%), predators (6.4%) and omnivores (11%) (Fig. 2).

  

Figure 1. Non-metric multidimensional scaling (nMDS) of the nematode genus composition of urea-free experiments (n=10). Black-filled circles and unfilled squares represent the nematode assemblages at the beginning (day 0) and at the end of the experiment (day 45), respectively. The two clusters mean that nematode genus composition between the two time periods (day 0 vs day 45) are significantly different from each other at p<0.05. 

  

Figure 2. Trophic grouping of nematodes across various urea concentrations per kg of soil (Control (0 mg/kg), T1 (19.67 mg/kg), T2 (39.33 mg/kg) and T3 (78.66 mg/kg) over different time intervals (day 0, 15, 30 and 45). 

The percentage dissimilarity between day 0 and day 45 was 49.93%, which can be attributed to the decrease of herbivores (Rotylenchus, Rotylenchulus and Belondira), predators (Mylonchulus, Iotonchus, Oxydirus and Eudorylaimus), bacterivores (Cephalobus) and omnivores (Crassolabium and Drepanodorylaimus). On the other hand, several genera also increased towards the end of the experiment, such as Rhabditis (bacterivore), Filenchus (fungivore) and Mesodorylaimus(omnivore) (Table 2).

Effect of urea on nematode communities

On day 15, the nematode abundance in T2 (39.33 mg/kg urea) was significantly higher than the other treatments (p<0.05) (Fig. 3A), while there was no significant difference (p>0,05) among the treatments on days 30 and 45. There was also no significant difference (p>0.05) in the number of genera (Fig. 3B), Shannon-Weiner Index (H’) (Fig. 3C), Simpson’s Evenness Index (Fig. 3D) and Maturity Index (MI) (Fig. 3E) between the treatments.

  

Figure 3. Descriptors of the nematode community per 100 g of soil across treatments over time. A) abundance B) richness expressed as number of genera, C) Shannon-Weiner Index (H’), D) Simpson’s Evenness Index and E) Maturity Index (MI). Data are mean values followed by the standard deviation (mean ± SD). Control (0 mg/kg), T1 (19.67 mg/kg), T2 (39.33 mg/kg) and T3 (78.66 mg/kg). 

Suitability of the microcosm approach

Field studies capture the complexity and variability of natural ecosystems, including interactions among multiple species and environmental factors. However, with such environmental variability, attributing observed effects caused by specific stressors can be challenging. On the other hand, microcosms offer a higher level of control over experimental conditions, allowing researchers to isolate and manipulate specific variables. By controlling specific factors, microcosm studies help in understanding the mechanisms underlying ecological processes. Thus, in this experiment, a microcosm approach was used to investigate the effects of urea on nematode communities.

The results showed a 28% reduction in abundance and 35% in the number of genera. These findings were consistent with previous studies involving nematodes in a terrestrial microcosm (Martikainen et al., 1998Martikainen E, Haimi J, Ahtiainen J, 1998. Effects of dimethoate and benomyl on soil organisms and soil processes-a microcosm study. Appl Soi Ecol, 9(1-3), 381-387. 10.1016/S0929-1393(98)00093-6; Martinez et al., 2019Martinez JG, Quiobe SP, Moens T, 2019. Effects of mercury (Hg) on soil nematodes: A microcosm approach. Arch Environ Con Tox, 77, pp. 421-431. 10.1007/s00244-019-00652-7). Although changes in physical and chemical properties have been shown to directly influence nematode abundance and diversity (Nisa et al., 2021Nisa RU, Tantray AY, Kouser N, Allie KA, Wani SM, Alamri SA, Alyemeni MN, Wijaya L, Shah AA, 2021. Influence of ecological and edaphic factors on biodiversity of soil nematodes. Saudi J Biol Sci, 28(5), pp. 3049-3059. 10.1016/j.sjbs.2021.02.046), the decline in abundance and genera, particularly herbivores, was influenced by the removal of plants (Martinez et al., 2019Martinez JG, Quiobe SP, Moens T, 2019. Effects of mercury (Hg) on soil nematodes: A microcosm approach. Arch Environ Con Tox, 77, pp. 421-431. 10.1007/s00244-019-00652-7; Ewald et al., 2020Ewald M, Glavatska O, Ruess L, 2020. Effects of resource manipulation on nematode community structure and metabolic footprints in an arable soil across time and depth. Nematology, 22(9), pp. 1025-1043. 10.1163/15685411-bja10009). Aside from providing nourishment to plant-feeding nematodes, plants can also alter soil conditions that can affect nematodes or can attract other organisms that serve as a food source of nematodes. Physical disturbance can also alter the nematode genera (Lazrova et al., 2021). A decline in abundance was observed among cp 4 and 5 genera, such as Mylonchulus, Iotonchus and Oxydirus (all predators) and Athernema, Crassolabium, Eudorylaimus, Drepanodorylaimus, Laimydorous and Paractinolaimus (all omnivores), suggesting that in addition to their high susceptibility to environmental disturbance (Sun et al., 2024Sun Z, Sun C, Feng X, Zhang T, Liu J, Wang X, Li S, Tang S, Jin K, 2024. Grazing alters the soil nematode communities in grasslands: A meta-analysis. J Environ Manage, 356, p. 120668. 10.1016/j.jenvman.2024.120668), predators and omnivores may also be sensitive to low food availability (Martinez et al., 2019Martinez JG, Quiobe SP, Moens T, 2019. Effects of mercury (Hg) on soil nematodes: A microcosm approach. Arch Environ Con Tox, 77, pp. 421-431. 10.1007/s00244-019-00652-7; Karuri, 2023Karuri H, 2023. Nematode community response to intensive tomato production in the tropics. Rhizosphere, 25, p. 100681. 10.1016/j.rhisph.2023.100681). The absence of a plant host may have caused the trophic groups to shift from herbivores to bacterivores after 45 days.

Nematodes are generally more abundant and diverse when the soil pH is close to neutral (approx. 6.5-7.5) (Li et al., 2023Li C, Wang X, Chen B, Wang L, Xie Z, Wang J, Yang Z, 2023. Fertilization restructures nematode assemblages by modifying soil pH in croplands of Northeast China. Front Environ Sci, 11, 1207379. 10.3389/fenvs.2023.1207379). Nisa et al. (2021Nisa RU, Tantray AY, Kouser N, Allie KA, Wani SM, Alamri SA, Alyemeni MN, Wijaya L, Shah AA, 2021. Influence of ecological and edaphic factors on biodiversity of soil nematodes. Saudi J Biol Sci, 28(5), pp. 3049-3059. 10.1016/j.sjbs.2021.02.046) reported that acidic soils have a lower abundance of nematodes compared to alkaline soils. Soil pH can also indirectly affect nematode diversity by influencing the abundance and diversity of other soil organisms that serve as a food source for nematodes. When pH decreases, soil nematodes increase the secretion of epidermal collagen to strengthen the cuticle which protect its body from osmotic pressure (Cong et al., 2020Cong Y, Yang H, Zhang P, Xie Y, Cao X, Zhang L, 2020. Transcriptome analysis of the nematode Caenorhabditis elegans in acidic stress environments. Front Psychol, 11, p. 1107. 10.3389/fphys.2020.01107). This suggests that nematodes can adapt to their environment and protect themselves from external stressors (e.g. osmotic pressure) and have the ability to strengthen their cuticle which can help them survive under adverse environmental conditions (Greiffer et al., 2022Greiffer L, Liebau E, Herrmann FC, Spiegler V, 2022. Condensed tannins act as anthelmintics by increasing the rigidity of the nematode cuticle. Sci Rep -UK, 12(1), p. 18850. 10.1038/s41598-022-23566-2). The decrease in pH by 0.5 after 45 days, along with physical disturbance, could also probably explain the decrease or disappearance of nematodes with high cp values. The omnivores and predators, e.g. Laimydorous and Paractinolaimus, were no longer found at the end of the experiment, however, there was an increase in Rhabditis and Filenchus. The response of Rhabditis in this study contradicts the results of Nisa et al. (2021Nisa RU, Tantray AY, Kouser N, Allie KA, Wani SM, Alamri SA, Alyemeni MN, Wijaya L, Shah AA, 2021. Influence of ecological and edaphic factors on biodiversity of soil nematodes. Saudi J Biol Sci, 28(5), pp. 3049-3059. 10.1016/j.sjbs.2021.02.046), who observed that bacterivorous (i.e. Rhabditis) nematodes are sensitive to decreasing pH. Other nematodes with similar cp values of 1 or 2 were found to be sensitive to changes in soil pH (Li et al., 2023Li C, Wang X, Chen B, Wang L, Xie Z, Wang J, Yang Z, 2023. Fertilization restructures nematode assemblages by modifying soil pH in croplands of Northeast China. Front Environ Sci, 11, 1207379. 10.3389/fenvs.2023.1207379), such as Cephalobus, whose abundance decreased on day 45. Cephalobus are known to feed on fungi and bacteria associated with plant roots, and their abundance can also be affected by the absence of plants. On the other hand, one genus (i.e. Ironus) with cp value of 4 was found at the end of the experiment but not recorded on day 0. This may have originated from an “egg bank” in the soil at the beginning of the experiment (Martinez et al., 2019Martinez JG, Quiobe SP, Moens T, 2019. Effects of mercury (Hg) on soil nematodes: A microcosm approach. Arch Environ Con Tox, 77, pp. 421-431. 10.1007/s00244-019-00652-7).

Effect of urea on nematode assemblages

The soil ecosystem can be disturbed by agricultural practices such as fertilisation (Ewald et al., 2022Ewald M, Rusch D, Rißmann C, Trost B, Theuerl S, Ruess L, 2022. Effects of irrigation and fertilization practice on soil nematode communities in arable land. Appl Soil Ecol, 177, p. 104546. 10.1016/j.apsoil.2022.104546). Adding fertilizers might make the soil more acidic or alkaline, and disrupt the soil ecosystem by altering nutrient availability.

Previous studies have shown a negative correlation between high levels of nitrogen fertiliser and soil nematode abundance (Liang et al., 2020Liang S, Kou X, Li Y, Lü X, Wang J, Li Q, 2020. Soil nematode community composition and stability under different nitrogen additions in a semiarid grassland. Global Ecol Conserv, 22, p. e00965. 10.1016/j.gecco.2020.e00965; Hu et al., 2022Hu J, Chen G, Hassan WM, Lan J, Si W, Wang W, Li G, Du G, 2022. The impact of fertilization intensity on soil nematode communities in a Tibetan Plateau grassland ecosystem. Appl Soil Ecol, 170, p. 104258. 10.1016/j.apsoil.2021.104258; Li et al., 2023Li C, Wang X, Chen B, Wang L, Xie Z, Wang J, Yang Z, 2023. Fertilization restructures nematode assemblages by modifying soil pH in croplands of Northeast China. Front Environ Sci, 11, 1207379. 10.3389/fenvs.2023.1207379; Wang et al., 2023aWang J, Wang H, Lin Q, Wu Y, He X, Chen X, Yan W, Zhao J, 2023a. Legume biological nitrogen fixation improves but chemical nitrogen fertilizer suppresses soil nematode communities in a Camellia oleifera plantation. Land Degrad Dev, 34(5), pp. 1403-1414. 10.1002/ldr.4542). As observed in all treatments (including the control) except T2 on day 15 (Fig. 3A), abundance decreased so this trend cannot be attributed to urea application alone. This may be caused by predation, resource competition or physical disturbance (Pothula et al., 2022Pothula SK, Phillips G, Bernard EC, 2022. Increasing levels of physical disturbance affect soil nematode community composition in a previously undisturbed ecosystem. J Nematol, 54(1). 10.2478/jofnem-2022-0022).

On the other hand, previous studies have reported an increase in nematode abundance and/or species richness with nitrogen application (Li et al., 2020Li J, Peng P, Zhao J, 2020. Assessment of soil nematode diversity based on different taxonomic levels and functional groups. Soil Ec Lett, 2(1), pp. 33-39. 10.1007/s42832-019-0019-5; Teshita et al., 2023Teshita A, Feng Y, Qian R, Wang X, Khan W, Gao Y, 2023. Alfalfa and maize intercropping enhances soil nematode structure and food web complexity in low-nitrogen soils. Appl Soil Ecol, 186, p. 104809. 10.1016/j.apsoil.2023.104809). This was observed in T2 (39.33 mg/kg), where a significantly higher abundance was recorded than in the control. Nitrogen fertiliser can stimulate microbial activity in the soil, leading to an increase in bacterial abundance which can serve as food for bacterivorous nematodes (Ni et al., 2024Ni X, Zhu X, Feng Q, Zhao D, Huang W, Pan F, 2024. Effect of Application Rates of N and P Fertilizers on Soil Nematode Community Structure in Mollisols. Agronomy, 14(3), p. 507. 10.3390/agronomy14030507). Hu et al. (2022Hu J, Chen G, Hassan WM, Lan J, Si W, Wang W, Li G, Du G, 2022. The impact of fertilization intensity on soil nematode communities in a Tibetan Plateau grassland ecosystem. Appl Soil Ecol, 170, p. 104258. 10.1016/j.apsoil.2021.104258) also reported that nitrogen input into the soil increases the abundance of nematodes. This was demonstrated in the present work in which an increase in the bacterivorous Rhabditis was observed.

The number of nematode genera and trophic groups tends to be affected by fertilisation (Hu et al., 2022Hu J, Chen G, Hassan WM, Lan J, Si W, Wang W, Li G, Du G, 2022. The impact of fertilization intensity on soil nematode communities in a Tibetan Plateau grassland ecosystem. Appl Soil Ecol, 170, p. 104258. 10.1016/j.apsoil.2021.104258). The effects of fertilisation on nematode genera and trophic groups can be complex and variable and depend on various factors, such as the type and amount of fertiliser applications, the soil type and the original composition of the nematode community. Bacterivorous nematodes tend to respond more to fertilisation than other trophic groups (Wang et al., 2023bWang J, Zhao X, Wei K, Tang J, Yuan C, Jin B, Sun X, Zhu B, 2023b. Correlations Between the Soil Bacterial-Feeding Nematodes, Bacteria, and Nitrogen in the Cropland of the Upper Yangtze River, China. J Soil Sci Plant Nut, 23(4), pp. 5840-5849. 10.1007/s42729-023-01443-9) (Fig. 2). On the other hand, fertilisation can have a negative impact on fungal populations due to the increased growth of microbial groups, leading to a decline in the abundance of fungivorous nematodes (Zhang et al., 2024Zhang H, Tian M, Jiang M, Yang J, Xu Q, Zhang Y, Ji M, Yao Y, Zhao C, Miao Y, 2024. Effects of nitrogen and phosphorus additions on soil nematode community of soybean farmland. Soil Ecol Lett. 6(2), p. 230200. 10.1007/s42832-023-0200-8). Moreover, the availability of nutrients, especially nitrogen, can also play a role in the shift from herbivorous to bacterivorous nematodes (Wilschut et al., 2021Wilschut RA, Geisen S, 2021. Nematodes as drivers of plant performance in natural systems. Trends Plant Sci, 26(3), pp. 237-247. 10.1016/j.tplants.2020.10.006). In contrast, bacterivorous nematodes can obtain nitrogen from the bacteria they consume, which may allow them to thrive in nutrient-poor environments (Mayrhofer et al., 2021Mayrhofer N, Velicer GJ, Schaal KA, Vasse M, 2021. Behavioral interactions between bacterivorous nematodes and predatory bacteria in a synthetic community. Microorganisms, 9(7), p. 1362.10.3390/microorganisms9071362; Shokoohi, 2024Shokoohi E, 2024. Interactions of Free-Living Nematodes and Associated Microorganisms with Plant-Parasitic Nematodes. Sustainable Management of Nematodes in Agriculture, Vol. 2: Role of Microbes-Assisted Strategies (pp. 127-147). 10.1007/978-3-031-52557-5_5). The results on day 0 showed that herbivores dominated the trophic groups, but as the days progressed, the dominant trophic groups shifted to bacterivores. In addition, the percentage of fungivores did not vary greatly from day 0 to day 45 in all treatments, including the control.

The average dissimilarity was 50%, due to the significant decline and disappearance of certain herbivorous and predatory species with cp 4-5 (Rotylenchus, Iotonchus, Belondira, Nygolaimus and Mylonchulus). These genera are known to be sensitive to disturbance (Krashevska et al., 2019Krashevska V, Kudrin AA, Widyastuti R, Scheu S, 2019. Changes in nematode communities and functional diversity with the conversion of rainforest into rubber and oil palm plantations. Front Ecol Evol7, p. 487. 10.3389/fevo.2019.00487) and the experimental set-up could be a disturbance to these nematodes. It was also observed that the bacterivorous community (Rhabditis) increased steadily across all treatments, including the control, as Rhabditis is known for its tolerance to pollutants (Allouche et al., 2020Allouche M, Nasri A, Harrath AH, Mansour L, Beyre H, Boufahj F, 2020. Migratory behavior of free-living marine nematodes surrounded by sediments experimentally contaminated by mixtures of polycyclic aromatic hydrocarbons. J King Saud Univ Sci32(2), pp. 1339-1345. 10.1016/j.jksus.2019.11.025). Ni et al. (2024Ni X, Zhu X, Feng Q, Zhao D, Huang W, Pan F, 2024. Effect of Application Rates of N and P Fertilizers on Soil Nematode Community Structure in Mollisols. Agronomy, 14(3), p. 507. 10.3390/agronomy14030507) found that fertilised plots were associated with an increase in the abundance of bacterivorous and suppressed the omnivorous nematodes. It suggests that the addition of nutrients (including NPK) can influence the nematode community by increasing the abundance and diversity of their microbial food sources in the soil. However, in this study the concentrations of urea used had no effect on the structure of the nematode community. This could be due to the soil’s natural ability to store, immobilise and degrade pollutants (Zhang et al., 2019Zhang B, Zhang L, Zhang X, 2019. Bioremediation of petroleum hydrocarbon-contaminated soil by petroleum-degrading bacteria immobilized on biochar. Rsc Adv, 9(60), pp. 35304-35311. 10.1039/C9RA06726D). In addition, since this is the first application of urea, the effects on the nematode community might not yet be noticeable. However, long-term and continuous application of urea over a longer period of time beyond the number of applications in this study, may alter the structure of the nematode community.