Research Article
Functional and Numerical Response of the Predatory Mite, Typhlodromus Setubali on Immature Stages of Panonychus Ulmi (Phytoseiidae: Tetranychidae)
Said Ouassat* and Latifa Allam
Corresponding Author: Said Ouassat, Laboratory of Health and Environment, Crop protection team, Department of Biology, Faculty of Sciences, Moulay Ismail University, Zitoune, 11201, Meknes, Morocco
Received: April 24, 2021; Revised: May 15, 2021; Accepted: May 17, 2021 Available Online: May 19, 2021
Citation: Ouassat S & Allam L. (2021) Functional and Numerical Response of the Predatory Mite, Typhlodromus Setubali on Immature Stages of Panonychus Ulmi (Phytoseiidae: Tetranychidae). Int J Biopro Biotechnol Advance, 7(3): 326-331.
Copyrights: ©2021 Ouassat S & Allam L. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Share :
  • 419

    Views & Citations

The European red spider mite, Panonychus ulmi (Koch) (Acari: Tetranychidae), continues to be the most serious mite pest of apple, Malus pumila L. Borkh (Rosacea) and can be very damaging to other crops. For a successful biological control using the predatory mite, T. setubali Dosse (Acari: Phytoseiidae), the predation (functional response), as well as the oviposition of gravid females of predator (numerical response), were investigated at different densities of P. ulmi immatures. The logistic regression revealed a type II functional response of T. setubali. The attack rate (α) and handling time (Th) of predator, were 0.029 prey/day and 0.573 day/prey, respectively. The oviposition of adult females increased above the prey density of 6.54 mites/cm2, with the highest daily fecundity of 2.13 eggs/ day. The prey density needed to achieve half the maximum response was estimated to be 15.18 prey, corresponding to 2.15 mites/ cm2. At lower densities of prey, the proportion Ne/N0 was higher. Therefore, T. setubali can be accurately a successful biological control in an augmentative release program. Further mite-based investigations are needed to refine the findings.

Keywords: Attack rate, Handling time, Typhlodromus setubali, Functional response


Panonychus ulmi Koch (Acari: Tetranychidae), commonly named European red spider mite, is the major phytophagous mite of various crops in the world [1]. Like other mite pests, P. ulmi feed on the cellular contents of leaves, compromising the canopy’s photosynthetic function and causing premature leaves abscission [2]. Based-chemical management programs assigned to insect and mite pests are a double edged, most spider mites of the family Tetranychidae, including P. ulmi, can quickly develop resistance against frequently used acaricides [3]. Facing limitation of the use of chemical control, the use of natural enemies, especially predatory mites of the family Phytoseiidae, has made an increased great interest. The use of phytoseiids is not recent [4], and this taxon groups more than 2 500 valid species, widely distributed around the world [5].

Biological control of the red spider mite, has been established on many cultures, with different phytoseiid species. The predatory mite, T. setubali Dosse (Acari: Phytoseiidae) is widely distributed in Morocco. It was reported on different crops, in many surrounding countries of the Mediterranean basin [6]. Based-field experiments during three seasons of biological control in Morocco, T. setubali is a valuable biological control agent of P. ulmi on different apple varieties. The main objective of this study is to apprehend the functional response of the predaceous mite, T. setubali to various densities of P. ulmi, whether it is of type I, II, or type III. Although lacking bibliographic information related to the development and feeding style of T. setubali, functional response experiments conducted on increasing densities of P. ulmi, can be helpful to improve P. ulmi management programs through understanding the specific behaviour of this species.

 The functional response is one of the most commonly used mathematical models to describe predator-prey interactions [7-9]. It is defined as the relationship between the prey density available and the number of preys attacked by a single predator during a given time interval. In biological control, the decision about which natural enemy to use is strongly based on the results of such an analysis, which leads to draw biological conclusions in combination with independent sources of information [10].

Holling [8] has defined three types of functional responses. In Type I, the number of consumed preys increases linearly to a maximum, then remains constant as increasing prey density. In Type II, the number approaches the asymptote hyperbolically as prey densities increases. Whereas in Type III, the number of preys consumed approaches it as a sigmoid curve. Therefore, a clear distinction between type II and III based on predation is rather difficult [11]. Many previously published results reported that type II is the commonly frequent type among predators, especially, among predaceous mites belonging to the family Phytoseiidae [12,13].

Although several regression-based models describing the predation rate of a biological control agent, have been published, the choice of one model to another constitutes a theoretical debate in terms of a complete analysis of data set [14]. Whenever possible, the functional response investigation remains a crucial method in population ecology, it helps to improve the knowledge of the change in prey- predator system, for a successful biological control. Overall, Modeling of the relationship between prey and its predator, has progressed considerably to more mathematically fitted models, which correct errors related to certain classical aspects of calculation [15].


Rearing of mites

Panonychus ulmi Koch was collected from apple trees in Oulmes region, Morocco, and reared on green bean plants Phaseolus vulgaris L. (Leguminosae). To obtain the seven densities that have been using in functional response experiments, the oviposition of gravid adult females of mite was taken into account. Prey numbers of 2, 4, 8, 16, 32, 64, and 128 immatures, at the same age, were obtained by transferring 1, 3, 5, 10, 16, 34, and 54 females, respectively, onto the leaf discs. Females were allowed to lay eggs for 24 hours, then were removed. The eggs were followed daily once until hatching and emergence of immature stages, which will be used during experiments as food for predators.

The initial population of T. setubali Dosse was obtained from Riyad-Fruit orchard, located at Tiddas, Morocco. The colony of predator was kept in rearing unit and P. ulmi infected bean leaves were daily added. The rearing unit consists of black plastic tiles, placed on a floating sponge into a plastic dish full water (40×27×17 cm). Borders of the black plastic support were covered with a wet wide band of Kleenex, to prevent mite escape.

Rearing of prey and predator colonies were carried out in the culture chamber at 26 ± 1°C, 65 ± 5 % RH, and 16:8 h (L: D) as photoperiod, for four generations before the beginning of the experiments.

To obtain same age-females, a total of fifty gravid females were individually transferred from the stock colony to bean leaf discs (diameter of 3 cm), which containing P. ulmi immatures as food diet. The females were allowed to lay eggs for 12h and then removed. After hatching, newly emerged larvae of T. setubali were followed until adult stage. Subsequently, males and females were paired onto new bean leaf discs until copulation success, and immediately separated. Homogenous cohorts of gravid females were then obtained for stating functional response tests.

Functional response study

To assess the functional response of T. setubali on P. ulmi immatures, 24 h starved adult females were assigned individually to one of seven prey densities (2, 4, 8, 16, 32, 64 and 128 prey), corresponding to 0.28, 0.56, 1.13, 2.26, 4.52, 9.05 and 18.11 mites/ cm2, respectively. After 24h, the predator females were removed and the number of preys killed was counted.

Because of predatory mites spend more time searching and handling prey, even if confined into a smaller space such as a leaf disc, the model does not detect the initial increase in attacks at the lower prey densities [14]. To overcome this model limitation, our experimental design involved ten replicates for mite densities of 0.28, 0.56, 1.13 mites/ cm2, seven replicates for 2.26, and 4.52 mites/cm2, and five for mite densities of 9.05 and 18.11 mites/cm2. Overall, seven prey density levels were tested for one predator female. Controls (without predators) were set up for each prey density.

Numerical response study

The numerical response is the number of offspring produced by female, fed with a given prey density during a time interval [16]. To quantify the numerical response of T. setubali, the oviposition of females was daily followed. The number of eggs laid by females at seven constant densities of P. ulmi, was recorded every 24 h, for 8 consecutive days.

Statistical analysis

Data analysis was performed in R program ver. 3.5.3. [17]. Functional response was determined according to Juliano’s procedure [18]. The logistic regression adjusting a polynomial equation (1) of the proportion of prey attacked (Ne), as a function of the initial prey density (N0), was used to estimate the linear, quadratic and cubic coefficients and therefore, determined the shape of the functional response curve of T. setubali. The type of response was determined by the signs of the linear and quadratic coefficients (P1 and P2). If the linear coefficient is negative (P1 < 0), it describes a type II functional response. If P1 > 0 and P2 < 0, it presents a type III functional response [19].
Where N0 is the initial number of prey, Ne is the number of prey eaten, Ne/ N0 is the probability of being attacked and P0, P1, P2, and P3 are the intercept, linear, quadratic, and cubic coefficients, respectively.

The attack rate (α) and handling time (h) can be determined by using the Holing disc equation or Rogers's random predator equations of type II (2) and type III (3) [20], namely RRPE-II and RRPE-III, respectively.

Rogers’s random predator equations include an attack exponent (q), to describe the per capita prey consumption in low prey densities and overcome the prey depletion at the end of experiments. Even if the number of attacked prey (Ne) appears on both sides of equation (2), the fit of data was performed by using iterative Newton’s method, as an alternative to LambertW function, which allows an explicit solution of the implicit RRPE-II [20]. For the functional response type III, a fitted version (3) of the original model has presented [22]. Both equations make it possible to predict how prey will be depleted over time during functional response experiments.

To describe the numerical response, the relationship between the fecundity of T. setubali females and the prey density available, was fitted by using a hyperbolic model based on this equation:

With (y) is the daily oviposition of females of T. setubali at the prey densities (x), (a) is the maximum daily oviposition, whereas (b) is the number of preys needed to reach half the maximum response.


Functional response

The logistic regression delivered a negative linear coefficient (P1 < 0) for P. ulmi densities, suggesting that T. setubali performed a type II functional response, assuming that the predation rate of adult females of T. setubali increased as prey density increased and levelling off to a maximum of 18 prey (Figure 1a). At the low prey densities ranged from 2 to 16, more than 90 % of P. ulmi immatures was consumed, and as its density increased, the proportion of prey attacked decreased; the rate of consumption tended to stabilize beyond a prey density of 64 mites (Figure 1b).

Estimates of the logistic regression coefficients are given in Table 1. The results revealed that P0, P1, P2, and P3 are significantly different from 0. Further, as the linear coefficient, P1 is negative (P = 0.007) and a positive quadratic coefficient (P = 0.016), T. setubali exhibited a type II functional response. The logistic regression adjusting a polynomial equation of the proportion of prey attacked as function of prey density can be written as:

Functional response data of T. setubali were successfully fitted with the Rogers´s model [20]. Both parameters, the attack rate (α) and handling time (Th), were found to be significantly different from 0. The model explained 95 % of the total variation in data (Table 2). The highest predation rate was found by 34.48 prey/ day, and half of the maximum predation was reached at density of 46.24 prey, corresponding to 6.54 prey/ cm2.

The attack rate was 0.573 prey/ day and the handling time was found be Th = 0.029 days/ prey. However, the Th estimated for this species seems to be short, meaning that predators spent less time in capturing, killing, and digesting the prey. These findings were consistent with those obtained for Euseius concordis Chant, reared on Oligonychus ilicis McGregor at 25 °C (0.033 days/ prey) [13]. Similarly, Neoseiulus cucumeris Oudemans had attack rate and handling time on Thrips flavidulus Bagnall, whereas some predatory mites, as Neoseiulus barkeri Hughes and Euseius nicholsi Ehara & Lee, had higher parameters at 26 °C (0.042 prey/ day and 0.057 days/ prey, respectively) [23].

Feeding with P. ulmi, the predatory mites, Euseius finlandicus Oudemans and Amblyseius andersoni Chant, displayed a functional response of type II, when both predators have been exposed to 25°C and increasing P. ulmi densities [24]. Chileseius camposi González y Schuster fed with P. ulmi at 20 ± 2ºC, showed likewise a type II [25]. Moreover, the functional response of a predatory mite, can change from type II to III, and vice versa, when different factors described above intervene over time of experiments.

The functional response of a predatory mite alone, cannot explain the success or failure of a biological control [26]. Several determining factors can affect the behavior of a natural enemy, such as interactions with potential predators [27], properties of host plants [28], as well as abiotic factors, such as temperature [29]. In this study, although the functional response has been studied at restricted conditions of food, functional response parameters of T. setubali could change in presence of pollen during spring, then affects the capacity of predator. For example, the cattail pollen of Typha latifolia L., which reduced significantly the predation rate and search efficiency of Typhlodromus pyri Scheuten, over time when it was tested on P. ulmi [30].

Numerical response
To assess the numerical response of gravid adult females of T. setubali, the hyperbolic model (4) explained 89 % of the total variation in daily oviposition (Table 2). The maximum daily oviposition (m) of predator female was estimated to be 1.86 eggs/ day, while the prey density needed to achieve half the maximum response (n) was estimated to be 2.15 P. ulmi immatures (Figure 2).

The results showed that the oviposition of T. setubali females, increased curvilinearly with prey density until reaching a plateau at prey densities from 32 to 128 P. ulmi immatures, with approximately 2.15 eggs/ female (Figure 2). The prey density-oviposition relationship was positively correlated with the predation rate, because the transformation rate in phytoseiid mites depends strongly on the number of resources invested per egg [31].

Both functional and numerical responses of T. setubali to P. ulmi can be useful for an elementary knowledge, for use as biological control agent. In the case of two or more concurrent predators, the determination of the functional response allows comparing specific effectiveness of each one [32]. Therefore, these results are exempt to field conditions and other mortality factors, like pesticide treatment and potential presence of other natural enemies. Further experiments are needed to evaluate the functional and numerical response of T. setubali, for wide biological control of pest mites on various crops around the Mediterranean.


This study suggests a type II functional response of T. setubali on P. ulmi immatures, which can be an appropriate food of this species. In addition, gravid females of predator showed a significant oviposition rate at 26 °C.
Although the experimental conditions may not correspond to the situation in field, the present results can be helpful for a first evaluation of the functional response of T. setubali and its incorporating as a biological control agent on various crops.
  1. Thistlewood HMA, Bostanian NJ, Hardman JM (2013) Panonychus ulmi (Koch) European Red Mite (Trombidiformes: Tetranychidae), in Biological Control Programmes in Canada 2001-2012, Mason P. G. and Gillespie D. R., CABI. pp: 238-243.
  2. Candolfi MP, Wermelinger B, Boller EF (1993) Photosynthesis and transpiration of “Riesling × Sylvaner” grapevine leaves as affected by European red mite (Panonychus ulmi Koch) (Acari, Tetranychidae) feeding. J Appl Entomol 115(1-5): 233-239.
  3. Kramer T, Nauen R (2011) Monitoring of spirodiclofen susceptibility in field populations of European red mites, Panonychus ulmi (Koch) (Acari: Tetranychidae), and the cross-resistance pattern of a laboratory-selected strain. Pest Manag Sci 67(10): 1285-1293.
  4. McMurtry JA, Sourassou NF, Demite PR (2015) The Phytoseiidae (Acari: Mesostigmata) as Biological Control Agents, in Prospects for Biological Control of Plant Feeding Mites and Other Harmful Organisms, D. Carrillo, G. J. de Moraes, and J. E. Peña, Eds. Cham: Springer International Publishing. pp: 133-149.
  5. Demite PR, de Moraes GJ, McMurtry JA, Denmark HA, Castilho RC (2019) Phytoseiidae Database. 2019, Accessed: March 7, 2021. Available online at:
  6. Tixier M-S, Allam L, Douin M, Kreiter S (2016) Phytoseiidae (Acari: Mesostigmata) of Morocco: New records, descriptions of five new species, re-descriptions of two species, and key for identification. Zootaxa 4067(5): 501-551.
  7. Solomon ME (1949) The Natural Control of Animal Populations. J Anim Ecol 8(1): 1-35.
  8. Holling CS ()1959) The Components of Predation as Revealed by a Study of Small-Mammal Predation of the European Pine Sawfly. Can Entomol 91(5): 293-320.
  9. Chiou J-M, Müller H-G, Wang J-L (2004) Functional Response Models. Stat Sin 14(3): 675-693.
  10. Casas J, Hulliger B (1994) Statistical analysis of functional response experiments. Biocontrol Sci Technol 4(2): 133-145.
  11. Rosenbaum B, Rall BC (2018) Fitting functional responses: Direct parameter estimation by simulating differential equations. Methods Ecol Evol 9(10): 2076-2090.
  12. Bazgir F, Shakarami J, Jafari S (2020) Functional response of the predatory mite Amblyseius swirskii (Acari: Phytoseiidae) to Eotetranychus frosti (Tetranychidae) and Cenopalpus irani (Tenuipalpidae). Acarologia 60(1): 30-39.
  13. da Silveira EC, Reis PR, Siqueira MF, Toledo MA, Liska GR, et al. (2020) Functional response of Euseius concordis feeding on Oligonychus ilicis (Acari: Phytoseiidae, Tetranychidae). Exp Appl Acarol 14(14): 371-382.
  14. Okuyama T (2013) On selection of functional response models: Holling’s models and more. BioControl 58(3): 293-298.
  15. Zhang JF, Papanikolaou NE, Kypraios T, Drovandi CC (2018) Optimal experimental design for predator–prey functional response experiments. J R Soc Interface 15(144): 20180186.
  16. Reeve JD, Ayres MP, Lorio PL (1995) Chapter 17 - Host Suitability, Predation, and Bark Beetle Population Dynamics, in Population Dynamics, N. Cappuccino and P. W. Price, Eds. San Diego: Academic Press. pp: 339-357.
  17. R Core Team (2019) A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
  18. Juliano SA (2001) Nonlinear curve fitting: predation and functional response curves, In: Scheiner, S.M. and J. Gurevitch (eds.)-Design and Analysis of Ecological Experiments, New York: Oxford University Press. pp: 178-196.
  19. Rutherford A (2000) Introducing Anova and Ancova: A GLM Approach. SAGE.
  20. Rogers D (1972) Random Search and Insect Population Models J Anim Ecol 41(2): 369-383.
  21. Pritchard D (2017) Frair: Tools for Functional Response Analysis.
  22. Hassell MP, Lawton JH, Beddington JR (1977) Sigmoid Functional Responses by Invertebrate Predators and Parasitoids. J Anim Ecol 46(1): 249-262.
  23. Yao H, Zheng W, Tariq K, Zhang H (2014) Functional and Numerical Responses of Three Species of Predatory Phytoseiid Mites (Acari: Phytoseiidae) to Thrips flavidulus (Thysanoptera: Thripidae). Neotrop Entomol 43(5): 437-445.
  24. Koveos DS, Broufas GD (2000) Functional response of Euseius finlandicus and Amblyseius andersoni to Panonychus ulmi on apple and peach leaves in the laboratory. Exp Appl Acarol 24(4): 247-256.
  25. Sepúlveda F, Carrillo R (2008) Functional response of the predatory mite Chileseius camposi (Acarina: Phytoseiidae) on densities of its prey, Panonychus ulmi (Acarina: Tetranychidae). Rev Biol Trop 56(3): 1255-1260.
  26. Gangisaffat F, Perring TM (2015) Prey stage preference and functional response of the predatory mite Galendromus flumenis to Oligonychus pratensis. Biol Control 82: 40-45.
  27. Sentis A, Hemptinne JL, Brodeur J (2013) How functional response and productivity modulate intraguild predation. Ecosphere 4(4): 46.
  28. Barbosa MFC, Poletti M, Poletti EC (2019) Functional response of Amblyseius tamatavensis Blommers (Mesostigmata: Phytoseiidae) to eggs of Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) on five host plants. Biol Control 138: 104030.
  29. Jalali MA, Tirry L, De Clercq P (2010) Effect of temperature on the functional response of Adalia bipunctata to Myzus persicae. BioControl 55(2): 261-269.
  30. Wei Q, Walde SJ (1997) The functional response of Typhlodromus pyri to its prey, Panonychus ulmi: the effect of pollen. Exp Appl Acarol 21(10): 677-684.
  31. de Moraes GJ, McMurtry JA (1985) Comparison of Tetranychus evansi and T. urticae [Acari: Tetranychidae] as prey for eight species of phytoseiid mites. Entomophaga 30(4): 393-397.