Spatial Distribution of Nematodes to Improve their Sampling and Management Decision

Mahfouz M.M. Abd-Elgawad

Plant Pathology Department, National Research Centre, Elbehoth St, Dokki 12622, Giza, Egypt.

Received | August 16, 2023; Accepted | October 22, 2023; Published | November 15, 2023

*Correspondence | Mahfouz M.M. Abd-Elgawad, Plant Pathology Department, National Research Centre, Elbehoth St, Dokki 12622, Giza, Egypt; Email: [email protected]

Citation | Abd-Elgawad, M.M.M., 2023. Spatial distribution of nematodes to improve their sampling and management decision. Pakistan Journal of Nematology, 41(2): 144-152.

DOI | https://dx.doi.org/10.17582/journal.pjn/2023/41.2.144.152

Keywords | Spatial distribution, Index of dispersion, Sampling and management, Ecological modeling, Sample size, Count transformation

Copyright: 2023 by the authors. Licensee ResearchersLinks Ltd, England, UK.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Introduction

Challenges for advanced agricultural systems comprise producing sufficient food while decreasing negative effects on the environment and utilizing our resources as best we can in durable approaches. Such challenges necessitate sound solving of problems linked to the full spectrum of losses and limitations in crop yields. Plant-parasitic nematodes (PPNs) can inflict substantial agricultural losses globally. A recent revision of such losses caused by PPNs for the 40 life-sustaining crops with significant values for food and export reported an average of 13.5% losses; equalled USD 358.24 billion annually (Abd-Elgawad and Askary, 2015; Abd-Elgawad, 2020a; 2022). Although various aspects of PPN biology and ecology are basic to optimize nematode management, major challenges are materialized in interpreting the intricacy of agroecosystem structures and dynamics into quantifiable metrics. Therefore, in the last few decades, novel technologies and spatial/ temporal statistics were established to provide definite tools for analyzing point distributional patterns. They can upgrade detection and identification of spatial heterogeneity such as clustering or gradient (Perry et al., 2002). Campos-Herrera et al. (2013) reviewed a few of these patterns that have general applications; for nematodes and other organisms. The reviewed patterns are also implied in issues related to catching where multitrophic reactions disclose, and how to characterize and explicate them.

The present study focuses on nematode distribution patterns in terms of their reflections on the economic significance of these nematodes especially from the standpoint of pest management. The study briefly addresses the basic models of the frequency distribution of nematode counts from a sample as a rationale for practical applications of additional related models. Other developed models to tackle ecological and biological issues of nematode distribution along with their merits and demerits are also presented. Related modifications of applying nematode distribution models to derive more accurate and appropriate conclusions for nematode management are discussed herein.

Using basic mathematical models in nematology

Although sampling precision and accuracy are keys to sound nematode research and effective management decision, much experimentation pass over these terms without giving satisfactory attention. For sound nematological research, Salama and Abd-Elgawad (2010) reviewed basic models of the three possible relations between the variance (Ơ2) of nematode counts and their arithmetic mean (µ) of a population. These are: (1) Positive binomial as a comparative model when the variance is significantly less than the mean (Ơ2 < µ), (2) Poisson series as the basic framework when the variance is almost equal to the mean (Ơ2 = µ), and (3) Negative binomial as the most flexible of several possible equations when the variance is significantly larger than the mean (Ơ2 > µ). These models are rigidly probability distributions and are converted to frequency distributions by multiplying each probability by the nematode-sample size. Advantageously, when the frequency distribution of nematode numbers in a sample fits one of these models, then: (1) the spatial dispersion of the population (ecological meaning) can be depicted in mathematical terms, (2) errors of population parameters can be assessed, (3) spatial and temporal shifts in nematode population density can be compared, and (4) the effect of ecological factors can be estimated (Herve et al., 2005; Duncan and Phillips, 2009; Salama and Abd-Elgawad, 2010; Rumiani et al., 2021). In this context, Wheeler et al. (1994) found that Meloidogyne hapla populations were best depicted by the negative binomial distribution, showing clumped distribution of its populations. Also, Poisson distribution was applicable to two Pratylenchus penetrans populations, indicating their random spatial distribution (Campbell and Noe, 1985; Wheeler et al., 1994). As nematodes are usually not evenly distributed, the related positive binomial is not of relevant usage for characterizing natural nematode distribution patterns. Nonetheless, positive binomial may be used for setting probability levels in detecting a nematode-infested unit for a range of lot/sample sizes over several infestation levels (Salama and Abd-Elgawad, 2003).

Models applied to measure spatial structure of nematodes: their merits and demerits

Sampling plans reviewed by Duncan and Phillips (2009) remains chiefly current as general sampling and extraction techniques have shifted little in the subsequent years (Abd-Elgawad, 2017). Therefore, employing models applied to other organisms to characterize also the nematode distribution patterns have been reported and developed (Been and Schomaker, 2013; Campos-Herrera et al., 2013; Berg et al., 2014; Gorny et al., 2020; Abd-Elgawad, 2021). Such applicable models in nematology can help in unraveling the complexities of ecosystem dynamics into quantifiable metrics via their indices of nematode dispersion. Also, they may open new advantageous avenues such as developing nematode-optimum sample size and use nematode-data transformation to meet assumptions necessary for parametric statistical analyses (McSorley et al., 1985; Duncan et al., 1989; Abd-Elgawad, 1992; Abd-Elgawad and Hasabo, 1995; Ghaderi et al., 2012; Rumiani et al., 2021).

Statistical methods for quantifying the spatial/temporal ordering of ongoing nematode population data rely mainly on fitting the above-mentioned frequency distributions, but applying other fine-tuning models have been a continuous process in nematology. These may comprise using definite models to establish the relationship between mean and variance of the sampled nematodes such as Taylor’s power law, TPL (Taylor, 1961), mapping and quantifying spatial autocorrelation (Campbell and Noe, 1985; Cressie, 1989, 1993), developing two-dimensional maps of spatial patterns for rather immobile organisms (Campos-Herrera et al., 2013; Diggle, 2013), and applying X, Y coordinates in a two dimensional space (Perry, 1995, 1998) to quantify nematodes-distribution patterns (Spiridonov et al., 2007; Kabir et al., 2018; Abd-Elgawad, 2020). Gorny et al. (2020) stressed that assumptions are often made on spatial and spatio-temporal distributions, lags, and related nematode-population parameters. Such assumptions are intended to circumvent the potential for high spatial variability and the considerable expenditure of resources as well as time devoted to gather, transfer, and analyze nematode samples to identify and quantify their populations. Because these assumptions can influence both accuracy and precision of assessing nematode populations and, eventually, management decisions (Robertson and Freckman, 1995; Steinberg and Kareiva, 1997; Duncan and Phillips, 2009), several options were adopted to solve such issues. From these options, TPL model was utilized to assess an index of nematode distribution. Consequently, this model could compute the optimal sampling intensities of a polyspecific nematode community (6-8 samples per 2 × 2 m row plot) for estimating populations within 50% of the coefficient of variation of the mean (Wheeler et al., 1994). These intensities comprised nematode species with both clumped (M. hapla) and random (P. penetrans) distributions of their populations. In another study, the equation of TPL was regrouped to solve for sampling accuracy via defining several levels of accuracy and precision linked to a fixed, cost-determined, sample size (Abd-Elgawad, 2017). Such levels may furnish the decision-maker with more informative inputs on sampling techniques. On the other hand, information on the spatiotemporal attributes of soil-borne disease epidemics is functional for deducing pathogen introduction circumstances and epidemic development. Thus, this information may comprise shifts in nematode population levels and patterns as a reflection of field management as well (Campbell and Noe, 1985; Abd-Elgawad et al., 2016; Contina et al., 2018). Furthermore, analysis of spatial patterns may test assumptions about links among pathogen taxa such as affinity or co-occurrence (Pethybridge and Turechek, 2003), ecological biology, or impacts of abiotic factors such as edaphic elements on spatial patterns, crop losses, and multitrophic interactions among the existing taxa (Campos-Herrera et al., 2013).

Notably, TPL was published (Taylor, 1961, 1970) as rule of thumb. It is used as a model primed for nematology by Ferris (1984) to offer certain services based on calculating nematode dispersion indices (McSorley et al., 1985; Duncan et al., 1989; Abd-Elgawad, 1992; Abd-Elgawad and Hasabo, 1995; Ghaderi et al., 2012; Rumiani et al., 2021). It states that: the variance (S2) of a population is proportional to a fractional power (b) of the arithmetic mean (X): S2 = aXb or log S2 = log a + b log X (i), where a and b are population parameters, a depends chiefly on the sample size and b is an index of nematode dispersion (McSorley et al., 1985). When parameters a and b of TPL are known, the sample size N (or N*) can be calculated from:

N = (1/E)2 a (x̄)b-2 or N* = (tα[n-1]/D)2 a (x̄)b-2

(ii) where N or N* is the number of samples, tα[n-1] is the appropriate Student’s t value for confidence limits of 1- α and n - 1 degrees of freedom, and sampling reliability is defined in terms of the standard error to mean ratio (E) or the ratio of the half-width of the confidence interval to the mean (D) of the samples. Thus, TPL is useful in determining transformations (Taylor, 1970) and developing nematode sampling plans (Ferris, 1984; McSorley et al., 1985; Duncan et al., 1989). Insights on applying TPL could offer further conclusions for the outputs of nematode data. Duncan et al. (1989) distinguished between two groups of citrus orchards based on the slope value of TPL. This value for the orchards with smaller patches of Tylenchulus semipenetrans-infected trees was different (P < 0.05) from that from large-patch orchards. For optimum sample size, TPL calculations revealed that 12 citrus trees (samples) are adequate in the geographic survey and 11 trees in the temporal for T. semipenetrans levels of 1,000 juveniles and males/100 cm3 soil (at E = 0.20). On the contrary, such sample size in sporadically infested orchards was assessed to be 69 trees (Duncan et al., 1989). Also, Abd-Elgawad (2016) introduced improvements via iteration process for accurate calculation of nematode-optimum sample size (Table 1). Consequently, minimum number of Meloidogyne incognita samples, for example, needed to achieve a 25% level of reliability as defined in terms of confidence interval half-width to mean ratio (D) was 402 and 389 samples before (McSorley et al., 1985) and after (Abd-Elgawad, 2016) iteration, respectively (Table 1). Moreover, rearranging the TPL formulae to solve for the ratio of the half-width of the confidence interval to the mean of the nematode numbers rather than sample size was suggested to select an acceptable level of sampling reliability (Table 2).

The mechanism through which TPL comes out is still deliberately talked over or argued in the literature (Campos-Herrera et al., 2013; Xiao et al., 2015; Xu and Cohen, 2021). On the other hand, TPL usage is confirmed to be compatible with orders of many organisms and applicable to their distribution patterns in nature (Taylor, 2019). For instance, application of TPL revealed that its exponents varied systematically with potential drivers concerning distribution patterns in time and space of fishes in the North Sea and therefore gains useful ecological acquaintance of the targeted system (Cobain et al., 2019). Hence, these authors stressed that its usage backed the notion of detecting community spatio-temporal structures and dynamics as impacted by driving biotic/abiotic factors. In addition to being a descriptive model, there have been other practical utilizations of TPL in biology and ecology (Taylor, 2019).

 

Table 1: Minimum number of nematode samples needed to achieve a 25% level of reliability as defined in terms of standard error to mean ratio (E) or confidence interval half-width to mean ratio (D) with iteration*.

Mean count per sample	Number of samples via E	Student’s t-value	Number of samples via D	Reference
Criconemella spp.: The power law parameters a = 3.076, b = 1.218
10	8	2 (assumed)	33	Abd-Elgawad and Hasabo, 1995
		2.037 (n=33)	34
		2.035 (n=34)	34
100	1	2 (assumed)	5
		2.776 (n=5)	10
		2.262 (n=10)	7
		2.447 (n=7)	8
		2.365 (n=8)	8
Meliodogyne incognita: The power law parameters a = 4.77, b = 2.12
10	101	2 (assumed)	402	McSorley et al., 1985
		1.966 (n=402)	389
		1.966 (n=389)	389
5	93	2 (assumed)	370
		1.966 (n=370)	358
		1.967 (n=358)	358

 

*The t-value is either assumed as 2 for 95% confidence interval or iterated using its tabulated value from: http://www.danielsoper.com/statcalc3/calc.aspx?id=10 (Abd-Elgawad, 2016).

 

Table 2: Percentage level of accuracy as defined in terms of the standard error to mean ratio (E) and the ratio of the half-width of the confidence intervalx to the mean (D) for stratified random sampling of Meloidogyne spp. in fields of berseem clover in Egypt.

Cost of samples (US $)	Finance based number of samples	Mean nematode count per sampley	Level of accuracy/reliabilityz
			D	E
Taylor’s power law parameters: a = 3.483, b = 1.729
150	15	40	63%	29%
200	20	30	55%	26%
250	25	20	51%	25%
300	30	20	48%	23%
400	40	10	44%	22%
200	20	1	87%	42%
250	25	1	77%	37%
750	150	1	30%	15%

 

xThe t-value at 95% confidence interval was obtained from: http://www.danielsoper.com/statcalc3/calc.aspx?id=10. yBased on a sample size of 100 gm soil (Abd-Elgawad and Hasabo, 1995; Abd-Elgawad, 2016, 2017). zThe fractional values rounded up to nearest two decimals (i.e. percentage).

Such a reliability/ accuracy of nematode sampling is especially significant for decision-maker in case of limited fund. Therefore, to define exact levels of accuracy linked to a fixed, cost-determined, sample size could be calculated (Abd-Elgawad, 2016). These levels may also furnish the decision-maker with more informative data on sampling approaches and their accuracies (Abd-Elgawad, 2017). However, such statistical manipulations should lead to identical confidence profile that is below, at, or above the economic threshold level, in most cases. As identification of new nematode species/strains/races is tremendously continuing via accurate and rapid morphological, biochemical and molecular tools (Shao et al., 2023) where modern trend is further facilitated by setting deep learning patterns for their automatic identification (Shabrina et al., 2023), quite high confidence interval should be assured. Hence, Abd-Elgawad (2017) reported a 99% confidence band above action threshold of Tylenchulus semipenetrans population with mean= 2500 at a precision of 80% or more as attainable probability levels due to emerging trends in gathering/processing samples and identifying their nematodes. Ultimately, such high levels of sampling precision and accuracy may be tailored with the current advances and evolutions in processing soil sampling rate to extra sizes of nematode samples.

 

Table 3: Comparison of index of aggregation (Ia)* values of four different studies on entomopathogenic nematode (EPN) distributions using different approaches in various regions.

EPN studied population	Form of measured EPNs	Ia value	Comments (location)	Reference
Heterorhabditis bacteriophora-infective stage juveniles (IJs) applied uniformly, in one, or nine patches on Kentucky bluegrass	EPN-infected Galleria mellonella larvae over time	All mean values were less than one but differed (P ≤ 0.05) until 20 weeks, no more, after EPN application	The values suggest a more even distribution than a random one (New Jersey/ USA)	Wilson et al., 2003
Natural populations of Steinernema feltiae and S. affine in grassland plots	IJs assigned to one of 4 groups of increasing physiological age	The values ranged 1.27-1.45 with group II as the most aggregated one.	All values indicate aggregated distribution (Merelbeke/ Belgium)	Spiridonov et al., 2007
H. bacteriophora or S. carpocapsae-infected G. mellonella larvae applied within 24 h of initial IJ emergence to cultivated fields and adjoining grassy border plots	H. bacteriophora and S. carpocapsae-IJs recovered from G. mellonella larvae baits applied several times after the cadavers	Range < 1 to > 2. Mean values differed between EPN species in bait traps and between soil management regimes at 48 h and 16 days after placing the cadavers, respectively	Spatial distributions after dispersing from a grassy border into the adjacent planted field plots were more aggregated for H. bacteriophora than for S. carpocapsae (Ohio/ USA)	Bal et al., 2017
Natural populations of H. indica in citrus and mango grove	EPN-infected G. mellonella larvae	0.913	Ia refers to even, close to random, distribution (Giza/ Egypt)	Abd-Elgawad, 2020

 

*Ia = the observed value of distance to regularity/ the mean randomized value (Perry, 1995).

 

Campos-Herrera et al. (2013) reviewed geostatistical models as a different set that can simultaneously utilize sample values and locations to illustrate spatial patterns and assess values at unsampled locations (Clark, 2001). Geostatistical models can effectively use global positioning systems (GPS) to define the location/site of each sample unit (Goovaerts, 1999). Unfortunately, this set of models does not supply tests to estimate the statistical significance of the established patterns. On the contrary, Analysis by Distance IndicEs (SADIE) manifested in a software package (Perry, 1995, 1998; Winder et al., 2019) can estimate the statistical probability level of spatial links between organisms or the same organism at different times (Perry, 1995). Although SADIE is more realistic measure than the empirical TYL, it is still of more limited use than TYL in nematology. Therefore, examples of its useful applications are given (Table 3).

Admittedly, these distribution models with their related indices can better be applied to complement each other and demonstrate more aspects of their established patterns. In this context, Gorny et al. (2020) simultaneously used two indices to fix adequate sampling protocols and to define certain sites for nematicidal application. For nematode interaction with other soil microorganisms, Wu et al. (2019) utilized SADIE to record regulation of entomopathogenic nematodes (EPNs) by the soil fungus Fusarium solani. Thus, using such indices of dispersion to boost biocontrol potential or to save unnecessary costs in nematicidal usage is ideally consistent with IPM. Moreover, spatial patterns are more depictive of nematode distributions in samples collected far apart which will be more influenced by different microhabitats than samples taken near to each other. This notion could be supported by utilizing SADIE analyses and semi-variogram together to better understand spatial patterns and spatiotemporal dynamics of PPNs (Gorny et al., 2020) and EPNs (Stuart et al., 2015). To ease its applications, SADIE software program including its main indices and graphical demonstrations were made available and reviewed (Winder et al., 2019). SADIE was also combined with other techniques to examine soil food webs in citrus groves for developing novel biocontrol approach within an IPM plans (Campos-Herrera et al., 2012, 2013).

Complementary approaches (Lavigne et al., 2010; Li et al., 2012) are also important to optimize detection of spatial heterogeneity especially when SADIE cannot tackle it effectively. Such approaches are useful when clusters of the targeted organism are located on elongated or square domains and close to the edges of the studied areas. In such cases, these techniques can distinguish clusters with small radius and in smaller sample size than that of SADIE. They can also be adapted to suit the absolute position or the value of the counts.

Overall advantages of determining nematode distribution patterns

Initially, it goes without saying that the above-mentioned pros and cons for different models of nematode distribution should be taken into account. Additionally, distribution patterns of plant nematodes are always derived from sound sampling to generally achieve one or more of the following advantages: (i) Select plant species/genotype that best fit to certain sites/locations, (ii) Manipulate variable rate techniques for applying the exact volume of nematicides, (iii) Unravel the complexities of nematode-ecosystem dynamics into quantifiable metrics, (iv) Distinguish relationships between various organisms in the sampled area, (v) Exploit such relationships in time and space for alert integrated pest management (IPM) schemes (Abd-Elgawad, 2021). Stuart and Gaugler (1994) speculated that EPN aggregated distribution can have substantial ramifications at the natural polyspecific community levels by shifting the dynamics of predation, parasitism, and competition. Moreover, the related temporal and spatial (horizontal/vertical) distributions may be compared with other relevant biotic and physical force(s) for better design of IPM schemes, and (vi) Certain models can function for nematode-count transformation to meet not only assumptions for parametric statistical analysis but also to achieve sound and accurate treatment comparisons.

Conclusions and Recommendations

Plant-parasitic nematodes deserve sound management approaches. Research on spatial and temporal distribution patterns of the nematodes based on their authentic dispersion indices will optimize such approaches. It offers realistic marks to unravel the complexities of ecosystem dynamics into quantifiable metrics. These patterns can enable to leverage variable rate techniques for the used nematicides and accurately depict relationships between organisms in space and time for careful IPM while basically serve in the precise nematode-count transformation to perfect treatment comparisons.

Acknowledgements

The author acknowledges the financial support and the facilities offered by The National Research Centre via in house project no.13050112 entitled Pesticide alternatives against soil-borne pathogens and pests attacking economically significant export crops.

Novelty Statement

Advanced agricultural systems can manipulate spatial distributions of nematodes to improve their sampling and management decision. The application of concepts and models related to nematode distribution patterns presented herein may be harnessed to tackle ecological and biological issues of nematodes and lower unwanted environmental effects.

Author Contribution

Mahfouz M. M. Abd-Elgawad: Suggested the idea, collected the literature necessary for the manuscript, wrote and reviewed the manuscript. Also, the author read and approved the manuscript.

Abbreviations

IPM, Integrated pest management; RKNs, Root-knot nematodes; E, the standard error to mean ratio; D, the ratio of the half-width of the confidence interval to the mean; Ia, the observed value of distance to regularity divided by the mean randomized value; TPL, Taylor’s power law; S2, variance of a population; X, arithmetic mean; a and b, population parameters of Taylor’s power law; PPNs, Plant-parasitic nematodes; EPNs, Entomopathogenic nematodes.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Funding

Financial support was partially made by The National Research Centre, Egypt. The facilities offered by The National Research Centre via in-house project no.13050112 are appreciated.

Conflict of interests

The author has declared no conflict of interest.

References

Abd-Elgawad, M.M., 1992. Spatial distribution of the phytonematode community in Egyptian citrus groves. Rev. Nématologie, 14: 367-373.

Abd-Elgawad, M.M., and Hasabo, S.A., 1995. Spatial distribution of the phytonematode community in Egyptian berseem clover fields. Fund. Appl. Nematol., 18(4): 329-334.

Abd-Elgawad, M.M.M., 2016. Use of Taylor’s power law parameters in nematode sampling. Int. J. Pharm. Tech. Res., 9(12): 999-1004.

Abd-Elgawad, M.M., 2017. Accuracy and precision of phytonematode sampling plans. Agric. Eng. Int. CIGR J., 19: 6–15.

Abd-Elgawad, M.M.M., 2020. Can rational sampling maximise isolation and fix distribution measure of entomopathogenic nematodes? Nematology, 22(8): 907-916. https://doi.org/10.1163/15685411-00003350

Abd-Elgawad, M.M.M., 2021. Optimizing sampling and extraction methods for plant-parasitic and entomopathogenic nematodes. Plants, 10(4): 629. https://doi.org/10.3390/plants10040629

Abd-Elgawad, M.M.M., 2022. Understanding molecular plant-nematode interactions to develop alternative approaches for nematode control. Plants, 11: 2141. https://doi.org/10.3390/plants11162141

Abd-Elgawad, M.M.M., 2020a. Biological control agents in the integrated nematode management of pepper in Egypt. Egypt. J. Biol. Pest Cont., 30: 70. https://doi.org/10.1186/s41938-020-00273-6.

Abd-Elgawad, M.M.M., and Askary, T.H., 2015. Impact of phytonematodes on agriculture economy. In: Askary, T.H., and Martinelli, P.R.P. (Eds). Biocontrol Agents of Phytonematodes. CAB International, Wallingford, UK, pp. 3-49. https://doi.org/10.1079/9781780643755.0003

Abd-Elgawad, M.M.M., Koura, F.F.H., Montasser, S.A., and Hammam, M.M.A., 2016. Distribution and losses of Tylenchulus semipenetrans in citrus orchards on reclaimed land in Egypt. Nematology, 18: 1141–1150. https://doi.org/10.1163/15685411-00003020

Bal, H.K., Acosta, N., Cheng, Z., Grewal, P.S., and Hoy, C.W., 2017. Effect of habitat and soil management on dispersal and distribution patterns of entomopathogenic nematodes. Appl. Soil Ecol., 121: 48-59. https://doi.org/10.1016/j.apsoil.2017.08.018

Been, T.H., and Schomaker, C.H., 2013. Distribution patterns and sampling. In: Perry, R.N., and Moens, M. (Eds). Plant Nematology, 2nd ed. CAB International, St. Albans, UK, pp. 331-358. https://doi.org/10.1079/9781780641515.0331

Berg, W.V.D., Hartsema, O., and Nijs, L.J.M.F.D., 2014. Statistical analysis of nematode counts from inter-laboratory proficiency tests. Nematology, 16(2): 229–243. https://doi.org/10.1163/15685411-00002761

Campbell, C.L., and Noe, J.P., 1985. The spatial analysis of soilborne pathogens and root diseases. Ann. Rev. Phytopathol., 23: 129-148. https://doi.org/10.1146/annurev.py.23.090185.001021

Campos-Herrera, R., Ali, J.G., Diaz, B.M., and Duncan, L.W., 2013. Analyzing spatial patterns linked to the ecology of herbivores and their natural enemies in the soil. Front. Plant Sci., 4: 378. https://doi.org/10.3389/fpls.2013.00378

Campos-Herrera, R., El-Borai, F.E., and Duncan, L.W., 2012. Wide interguild relationships among entomopathogenic and free-living nematodes in soil as measured by real time qPCR. J. Invertebr. Pathol., 111: 126–135. https://doi.org/10.1016/j.jip.2012.07.006

Clark, I., 2001. Practical Geostatistics. Available on line at: www.uni.tuebingen.de/egpinfo/Pedometrie/ PracticalGeostatistics1979.pdf

Cobain, M.R.D., Brede, M., and Trueman, C.N., 2019. Taylor’s power law captures the effects of environmental variability on community structure: An example from fishes in the North Sea. J. Anim. Ecol., 88: 290–301. https://doi.org/10.1111/1365-2656.12923

Contina, J.B., Dandurand, L.M., and Knudsen, G.R., 2018. A spatial analysis of the potato cyst nematode Globodera pallida in Idaho. Phytopathology, 108: 988-1001. https://doi.org/10.1094/PHYTO-11-17-0388-R

Cressie, N., 1989. Geostatistics. Am. Statist., 43: 197-202. https://doi.org/10.1080/00031305.1989.10475658

Cressie, N., 1993. Statistics for spatial data. John Wiley and Sons, New York, USA, pp. 424. https://doi.org/10.1002/9781119115151

Diggle, P.J., 2013. Statistical analysis of spatial and spatio-temporal point patterns. Academic press, London, UK, pp. 300. https://doi.org/10.1201/b15326

Duncan, L.W., and Phillips, M.S., 2009. Sampling root-knot nematodes. In: Perry, R.N., Moens, M., and Starr, J.L. (Eds). Root-knot Nematodes. CAB International, St. Albans, UK, pp. 275–300. https://doi.org/10.1079/9781845934927.0275

Duncan, L.W., Ferguson, J.J., Dunn, R.A., and Noling, J.W., 1989. Application of Taylor’s power law to sample statistics of Tylenchulus semipenetrans in Florida citrus. J. Nematol. (Suppl.), 21: 707-711.

Ferris, H., 1984. Probability range in damage predictions as related to sampling decisions. J. Nematol., 16: 246-251.

Ghaderi, R., Hamzehzarghani, H., and Karegar, A., 2012. Sampling optimization for root lesion nematodes in the irrigated wheat fields of Marvdasht region, Fars, Iran. Nematol. Medit., 40: 3-10.

Goovaerts, P., 1999. Geostatistics in soil science: State-of-the-art and perspectives. Geoderma, 89: 1-45. https://doi.org/10.1016/S0016-7061(98)00078-0

Gorny, A.M., Hay, F.S., Esker, P., and Pethybridge, S.J., 2020. Spatial and spatiotemporal analysis of Meloidogyne hapla and Pratylenchus penetrans populations in commercial potato fields in New York, USA. Nematology, 23: 139-151. https://doi.org/10.1163/15685411-bja10034

Herve, G., Bertrand, B., Villain, L., Licardie, D., and Cilasand, C., 2005. Distribution analyses of Meloidogyne spp. and Pratylenchus coffeae sensu lato in coffee plots in Costa Rica and Guatemala. Plant Pathol., 54: 471-475. https://doi.org/10.1111/j.1365-3059.2005.01206.x

Kabir, M.F., Mwamula, A.O., Lee, J.K., Jeong, M., Lee, D.W., and Park, J.J., 2018. Spatial distribution of Heterodera trifolii in Chinese cabbage fields. J. Asia-Pacific Entomol., 21(2): 688-694. https://doi.org/10.1016/j.aspen.2018.04.011

Lavigne, C., Ricci, B., Franck, P., and Senoussi, R., 2010. Spatial analyses of ecological count data: A density map comparison approach. Basic Appl. Ecol., 11: 734-742. https://doi.org/10.1016/j.baae.2010.08.011

Li, B., Madden, L.V., and Xu, X., 2012. Spatial analysis by distance indices: An alternative local clustering index for studying spatial patterns. Methods Ecol. Evol., 3: 368-377. https://doi.org/10.1111/j.2041-210X.2011.00165.x

McSorley, R., Dankers, W.H., Parrado, J.L., and Reynolds, J.S., 1985. Spatial distribution of the nematode commun­ity on Perrine marl soil. Nematropica, 15: 77-92.

Perry, J.N., 1995. Spatial analysis by distance indices. J Anim. Ecol., 64: 303-314. https://doi.org/10.2307/5892

Perry, J.N., 1998. Measures of spatial patterns for counts. Ecology, 79: 1008-1017. https://doi.org/10.1890/0012-9658(1998)079[1008:MOSPFC]2.0.CO;2

Perry, J.N., Liebhold, A.M., Rosenberg, M.S., Dungan, J., Miriti, M., Jakomulska, A., Citron-Pousty, S., 2002. Illustrations and guidelines for selecting statistical methods for quantifying spatial patterns in ecological data. Ecography, 25: 578–600. https://doi.org/10.1034/j.1600-0587.2002.250507.x

Pethybridge, S.J., and Turechek, W.W., 2003. Analysis of the association among three viruses infecting hop in Australia. Plant Pathol., 52: 158-167. https://doi.org/10.1046/j.1365-3059.2003.00812.x

Robertson, G.P., and Freckman, D.W., 1995. The spatial distribution of nematode trophic groups across a cultivated ecosystem. Ecology, 76: 1425-1432. https://doi.org/10.2307/1938145

Rumiani, M., Habiballah Hamhehzarghani, H., Karegar, A., and Ghaderi, R., 2021. Optimization of citrus tree sampling pattern for estimating population of citrus nematode in the soil of infested orchards in Fars province, southern Iran. Crop Prot., 142: 105523. https://doi.org/10.1016/j.cropro.2020.105523

Salama, H.S., and Abd-Elgawad, M.M.M., 2003. Quarantine problems: An analytical approach with special reference to palm weevils and phytonematodes. Arch. Phytopathol. Plant Prot., 36: 41-46. https://doi.org/10.1080/0323540031000080137

Salama, H.S., and Abd-Elgawad, M.M.M., 2010. Spatial patterns of the red palm weevil and applied entomopathogenic nematode Heterorhabditis bacteriophora. Arch. Phytopathol. Plant Prot., 43(7): 689-699. https://doi.org/10.1080/03235400802144397

Shabrina, NH., Lika, RA., and Indarti, S., 2023. Deep learning models for automatic identification of plant-parasitic nematode. Artif. Intell. Agric. 7: 1-12. https://doi.org/10.1016/j.aiia.2022.12.002

Shao, H., Zhang, P., Peng, D., Huang, W., Kong, L-a., Li, C., Liu, E., and Peng, H., 2023. Current advances in the identification of plant nematode diseases: From lab assays to in-field diagnostics. Front. Plant Sci., 14: 1106784. https://doi.org/10.3389/fpls.2023.1106784

Spiridonov, S.E., Moens, M., and Wilson, M.J., 2007. Fine scale spatial distributions of two entomopathogenic nematodes in agrassland soil. Appl. Soil Ecol., 37: 192-201. https://doi.org/10.1016/j.apsoil.2007.06.009

Steinberg, K., and Kareiva, P., 1997. Challenges and opportunities for empirical evaluation of spatial theory. In: Kareiva, P., and Tilman, D. (Eds). Spatial ecology: The role of space in population dynamics and interspecific interactions. Princeton University Press, Princeton, Massachusetts, USA, pp. 318-331. https://doi.org/10.2307/j.ctv36zpzm.22

Stuart, R., Gaugler, R., 1994. Patchiness in Populations of Entomopathogenic Nematodes. J. Invertebr. Pathol., 64: 39–45. https://doi.org/10.1006/jipa.1994.1066

Stuart, R.J., Barbercheck, M.E., Grewal, P.S., 2015. Entomopathogenic nematodes in the soil environment: Distributions, interactions and the influence of biotic and abiotic factors. In: Campos-Herrera, R. (Ed.). Nematode pathogenesis of insects and other pests. Springer, Cham, Switzerland, pp. 97–138. https://doi.org/10.1007/978-3-319-18266-7_4

Taylor, L.R., 1961. Aggregation, variance and the mean. Nature, 189: 732-735. https://doi.org/10.1038/189732a0

Taylor, L.R., 1970. Aggregation and the transformation of counts of Aphis fabae Scop. On beans. Ann. Appl. Biol., 65: 181-189. https://doi.org/10.1111/j.1744-7348.1970.tb04577.x

Taylor, R.A.J., 2019. Taylor’s power law: Order and pattern in nature. Elsevier Academic Press, Cambridge, UK, pp. 1286.

Wheeler, T.A., Madden, L.V., Riedel, R.M., and Rowe, R.C., 1994. Distribution and yield-loss relations of Verticillium dahliae, Pratylenchus penetrans, P. crenatus, and Meloidogyne hapla in commercial potato fields. Phytopathology, 84: 843-852. https://doi.org/10.1094/Phyto-84-843

Wilson, M.J., Lewis, E.E., Yoder, F., and Gaugler, R., 2003. Application pattern and persistence of the entomopathogenic nematode Heterorhabditis bacteriophora. Biol. Cont., 26: 180-188. https://doi.org/10.1016/S1049-9644(02)00125-1

Winder, L., Alexander, C., Griffiths, G., Holland, J., Wooley, C., and Perry, J., 2019. Twenty years and counting with SADIE: Spatial analysis by distance indices software and review of its adoption and use. Rethink. Ecol., 4: 1-16. https://doi.org/10.3897/rethinkingecology.4.30890

Wu, S.-Y., El-Borai, F.E., Graham, J.H., and Duncan, L.W., 2019. Geospatial relationships between native entomopathogenic nematodes and Fusarium solani in a Florida citrus orchard. Appl. Soil Ecol., 140: 108–114. https://doi.org/10.1016/j.apsoil.2019.04.009

Xiao, X., Locey, K., and White, E.P., 2015. A process-independent explanation for the general form of Taylor’s law. Am. Natur., 186(2): 51-60. https://doi.org/10.1086/682050

Xu, M., and Cohen, J.E., 2021. Spatial and temporal autocorrelations affect Taylor’s law for US county populations: Descriptive and predictive models. PLoS One, 16(1): e0245062. https://doi.org/10.1371/journal.pone.0245062