https://orcid.org/0000-0002-6519-1069
Figures
Abstract
Amblyseius eharai is a generalist predatory mite that consumes spider mites, rust mites, thrips, and pollen, with a high adaptability to various plants. To better understand ecological and behavioral aspects of this species, we investigated its functional response to different stages of two-spotted spider mite, Tetranychus urticae. Furthermore, we compared its environmental adaptability with that of other referenced phytoseiids using a temperature-dependent model of the intrinsic rate of increase. We were able to calculate the functional response parameters of both sexes of A. eharai when preying on eggs or larvae of T. urticae and, for females only, when preying on the deutonymph of T. urticae. Among the various combinations tested herein, A. eharai females preying on T. urticae larvae had the highest attack rate and shortest handling time. For eggs of T. urticae, A. eharai showed a lower attack rate; however, its handling time for eggs was significantly shorter compared to other phytoseiids. Using T. urticae larva as a prey, the attack rate of female A. eharai was higher and the handling time of both sexes of this species was shorter than those of other phytoseiid mites. Amblyseius eharai populations can show maximum performance quickly due to this species’ lower optimal temperature for population growth (28.1°C) compared to other phytoseiid mites. Thus, we provided evidence that this predatory mite has the potential to be a new, effective biological control agent of greenhouse pests such as T. urticae due to its high predation capacity and low optimal temperature.
Citation: Park Y-G, Lee J-H, Lim UT (2021) Functional response of Amblyseius eharai (Acari: Phytoseiidae) on Tetranychus urticae (Acari: Tetranychidae). PLoS ONE 16(12): e0260861. https://doi.org/10.1371/journal.pone.0260861
Editor: Ramzi Mansour, University of Carthage, TUNISIA
Received: August 1, 2021; Accepted: November 17, 2021; Published: December 2, 2021
Copyright: © 2021 Park et al. 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.
Data Availability: All relevant data are within the paper and its Supporting information files.
Funding: This work was carried out with the support of Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ013362042021), Rural Development Administration, Republic of Korea. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Predatory mites are very effective biological control agents and are commonly used to control small-bodied pests such as spider mites, thrips, and whiteflies in agricultural crops [1, 2]. Phytoseiid mites are among the most extensively commercialized and frequently applied natural enemies for augmentative biological control in greenhouses [3–6]. Among them, Amblyseius swirskii (Athias-Henriot), Neoseiulus californicus (McGregor), and Phytoseiulus persimilis (Athias-Henriot) (Acari: Phytoseiidae) are widely commercialized species.
Amblyseius eharai (Amitai et Swirski) (Acari: Phytoseiidae) is a native species in East Asia and is a type III generalist predator [2, 3, 7]. It feeds on spider mites, eriophyid mites, thrips, and pollen [2, 3, 8–11] and has a lower optimal temperature for development and oviposition compared to other predatory mites such as A. swirskii, N. californicus, N. longispinosus (Evans), and N. womersleyi (Schicha) [11–15]. Park and Lee [11] reported that A. eharai did not oviposit at 33.2°C nor were its larvae able to develop at 35.9°C, while in comparison, A. swirskii, N. californicus, N. longispinosus, and N. womersleyi all showed better development and reproduction than A. eharai at 36.0°C [12–15]. Despite these ecological adaptations for lower temperature conditions, A. eharai is potentially a valuable biological control agent since it has a high adaptability to many types of plants such as deciduous trees, conifers, shrubs, herbs, and vine [16]. Especially, this predator would adapt easily to East Asian agricultural land to which it is native. However, there is limited information available on other ecological aspects of A. eharai, especially on its functional response, which can be important when developing a new species for use as a biological control agent. Knowledge of a species’ type of functional response to a particular pest can help evaluate its potential control efficacy in a particular environment [17]. Type II and III are the most common functional response types in predatory arthropods [18–23], but these patterns can be affected by the physiological status of either the predator or prey species [20, 21]. From the functional response type, it is possible to estimate a predator’s attack rate and handling time. Although promising parameter values for a predator’s functional response do not guarantee its success for use in a biological control program, they can be used to assess a predator’s potential relative to other species in the same predator guild.
The two-spotted spider mite, Tetranychus urticae (Koch) (Acari: Tetranychidae), is a worldwide pest that can cause severe damage to many horticultural crops in greenhouses and open fields [24–26]. Spider mites reduce crop photosynthesis by feeding on plant cells and are hard to control due to their high population growth rates on many crops [24, 27, 28] and their rapid development of resistance to pesticides [29, 30]. As an alternative to chemical insecticides, biological control agents including predatory mites, anthocorid bugs, and pathogenic fungi can be applied [31, 32]. For the successful biological control of T. urticae in farmland, agents used should be well adapted to the area’s environmental and ecological conditions [31, 33].
In this study, we evaluated functional response of A. eharai, a phytoseiid native to East Asia, on different development stages of T. urticae to explore the potential of this predator as a biological control agent. Additionally, the attack rate and handling time of A. eharai on eggs and larvae of T. urticae were calculated and compared with those of other phytoseiid mites used as biocontrol agents. Finally, to estimate where A. eharai might function well for biological control, we compared its environmental adaptability with other predatory mites using a temperature-dependent model of intrinsic rate of increase values obtained from the scientific literature.
Materials and methods
Laboratory rearing of predator and prey mites
Tetranychus urticae was reared in the Laboratory of Insect Ecology, Seoul National University using kidney bean plants as food. Two to three kidney beans were seeded in a pot (85 × 70 mm; diameter × height) and grown for two weeks before being used to rear T. urticae. Six to seven such potted plants were inoculated with mixed stages of T. urticae (over 1000 mites) and then placed on a small steel tray (365 × 320 × 50 mm; width × length × height). This small tray was then placed in a bigger steel tray (550 × 455 × 60 mm; width × length × height) that was flooded with water to prevent mite escape. About half of the plants in a small tray were replaced with new ones every 3–4 days. The stems and leaves of plants being removed were cut off and draped over the new plants to transfer their T. urticae individuals. Rearing was done at 25°C, 60–80% RH, and a photoperiod of 16:8 (L:D) h in an insectary.
Amblyseius eharai used in our experiments were from a colony started with mites collected from the overwintering grape buds of vineyards in Hwaseong, Korea (N 37° 10’ 23.5", E 126° 41’ 39.1"). A total 150 A. eharai females were collected from the buds and kept at 26–8°C, 60–80% RH, and a photoperiod of 16:8 (L:D) h in an incubator (50 × 50 × 50 cm; width × length × height; Hanbaek Scientific Technology, Bucheon, Korea). To provide food for the predatory mites, mixed stages of T. urticae on an infested kidney bean leaf from the rearing colony were placed in each Petri dish (100 × 42 mm; diameter × height) (SPL life science, Pocheon, Korea) used to rear predators. The Petri dishes with the predatory mites each contained an excised kidney bean leaf disc (70 mm diameter, upside down) (infested with T. urticae) on water-saturated cotton. Each Petri dish was inoculated with 10 to 15 A. eharai females. The kidney bean leaves infested with T. urticae in these dishes were replaced with new ones from the spider mite rearing colony once a week. When leaves with prey were replaced, all detectable A. eharai stages including eggs were transferred onto the new prey-infested leaves with a fine brush. At least 30 to 40 generations of A. eharai were produced in this way before the colony was used in this experiment.
Functional response experiment
To assess the functional response of A. eharai, the experiments were conducted at 26.6 ± 0.39°C, 72.2 ± 2.51% RH, and a photoperiod of 16:8 (L:D) h using adults of A. eharai of our rearing colony as described above. To obtain freshly emerged A. eharai adults, deutonymphs were collected from the rearing colony and isolated in a Petri dish under the same conditions as used for rearing. The immature development of these mites was checked daily, and newly molted adults of both sexes were moved to another Petri dish (held under the same conditions) and allowed to feed, mate, and lay eggs for three days. After three days, adults (now 3-day-old) of each sex were placed individually in Petri dishes with a single kidney bean leaf with no prey mites and starved for 24 h.
Smaller Petri dishes (50 × 15 mm; diameter × height; SPL Life Science, Pocheon, Korea) were used as test arenas for the experiments. In each test Petri dish, one kidney bean leaf disc (35 mm diameter) was placed upside down on water-saturated cotton. Preys (T. urticae) in various densities were released in these Petri dishes with predators (A. eharai). The densities of each prey stage and its replication numbers for the functional response experiments are presented in Table 1.
To obtain T. urticae eggs for experiments, five to ten adult female mites collected from the rearing colony were placed in each intended test dish, where they were allowed to lay eggs for 24 hours. Adult mites were then removed, and the number of eggs counted. The number of eggs per dish (Table 1) was adjusted by removing or adding eggs with a fine brush. To obtain larvae and protonymphs of T. urticae, 300 to 400 females of T. urticae were placed in each of two to three Petri dishes (100 × 42 mm; diameter × height), on a kidney bean leaf disc (70 mm diameter, upside down) placed upside down on water-saturated cotton. The T. urticae females were allowed to lay eggs for 48 hours, and then the females were removed. For three to seven days, daily observations were made to check the developmental stage of the T. urticae progeny. When mites were larvae or protonymphs, they were moved to the smaller (50 × 15 mm, diameter x height) test Petri dishes for experiments. The higher prey stages (deutonymphs and adult females of T. urticae) needed for the functional response experiments were collected directly from the rearing colony and moved to test Petri dishes. In the process of gathering prey mites, webbing of mites was excluded.
The starved adult females or males of A. eharai were individually placed in the test Petri dishes where the prey has been introduced. The A. eharai adults were allowed to attack prey for 24 hours, and the number of preys consumed were then counted. If the adult A. eharai died, that replicate was discarded.
Prey stage preference test (eggs vs. larvae of T. urticae)
For adult males and females of A. eharai we compared the numbers of prey consumed for eggs vs. larvae of T. urticae. This test was run under the same environmental conditions as the functional response experiment. The eggs and larvae of T. urticae were obtained by the same method as used in functional response experiments. Twenty-five eggs or larvae were transferred to the same kidney bean leaf disc (35 mm diameter, upside down) on water-saturated cotton in a test Petri dish (50 × 15 mm, diameter × height) with a fine brush. Adult females and males of A. eharai were obtained by the same method as used in functional response experiments, and they were placed individually in the test Petri dishes containing T. urticae eggs and larvae on a leaf disc. The predatory mites were allowed to feed on eggs and larvae of T. urticae for 24 hours, and then the numbers of eggs or larvae consumed were determined by counting the survivors. Eleven replications were run for females and ten for males.
Data analysis
The effects of prey density, life stage (egg vs. larva), and predator sex on the predation rate were analyzed by three-way analysis of variance (ANOVA) using PROC GLM in SAS (SAS Institute, Cary, NC) [34]. The effects of the prey density and predator sex on the predation of protonymphs, deutonymphs, and adult females were all analyzed by two-way ANOVA using PROC GLM in SAS (SAS Institute, Cary, NC) [34]. In the case of adult female predators, data on protonymph and adult female prey were analyzed using one-way ANOVA with PROC GLM in SAS (SAS Institute, Cary, NC) [34] due to lack of the male data in some density treatments. The numbers of consumed T. urticae eggs and larvae in the prey preference test were compared using a t-test with PROC TTEST in SAS (SAS Institute, Cary, NC) [34].
The functional response data were analyzed in two phases as per Juliano [35], using PROC GENMOD in SAS (SAS Institute, Cary, NC) [34]. Logistic regression was performed to assess the form of the functional response. The equation is:
(1)
where Ne is the number of prey consumed per predator, the initial number of prey is N0, and the likelihood of consumption is Ne/N0. The parameters include P0, P1, P2, and P3.
Using logistic regression, maximum likelihood estimates were obtained for parameters P0 to P3. To determine the functional response type, the logistic model parameter was estimated by the log-likelihood test. If P1 < 0, the functional response is type II. If P1 > 0 and P2 < 0, the functional response is type III [35]. In this study, all significant (P < 0.05) P1 values of A. eharai on different prey stages were negative. Thus, for type II functional responses, the random predator equation [36] was used to estimate functional response parameters such as an attack rate (α) and handling time (Th). The equation is:
(2)
where Ne is the number of preys consumed per predator during the test time T (24 h), N0 is the initial prey number. We compared the attack rate and handling time of A. eharai among the treatments at a 95% confidence interval.
Results
Predation rates of Amblyseius eharai on Tetranychus urticae life stages
Predation rates of A. eharai on T. urticae eggs and larvae were affected by prey density, prey stage, A. eharai sex, the prey density × stage interaction, and the prey density × A. eharai sex interaction (Table 2). However, predation rates were not significantly affected by the prey density × stage × A. eharai sex interaction (Table 2). Predation rates for A. eharai sexes on different densities of T. urticae larvae or eggs showed that as the initial prey density increased, the predation rate of A. eharai also increased (Table 3). The predation rate by female A. eharai was generally higher on larvae than on eggs of T. urticae (Table 3).
In the treatment of protonymph, deutonymph, and adult females of T. urticae to A. eharai, the predation rate was affected by prey density, sex of A. eharai, and their interaction (Table 2). In general, female predatory mites consumed more T. urticae than males (Table 4). Female A. eharai consumed significantly more prey at higher prey densities than at lower densities. However, such trend was not found in male A. eharai that rarely preyed on T. urticae adult female.
Prey preference of A. eharai between T. urticae eggs and larvae
The results of prey preference of A. eharai are presented in Fig 1. Both sexes of A. eharai consumed more T. urticae larvae than eggs (Female, T14.30 = 2.52, P = 0.024; Male, T18 = 2.16, P = 0.045). Female A. eharai consumed 22.7 larvae and 19.1 eggs of T. urticae during 24 hours. Male A. eharai consumed 11.8 larvae and 7.8 eggs of T. urticae during 24 hours.
*Means followed by the same letter in each sex are not significantly different at α = 0.05, t-test.
Functional response
In the treatment that showed significance of P1 value (Egg for females, P < 0.001; Egg for males, P = 0.003; Larva for females, P = 0.001; Larva for males, P < 0.001; Deutonymph for females, P = 0.015; Adult for males, P = 0.001), the P1 values for the relations between prey presented and prey consumed all were negative, indicating a type II functional response (Table 5, Fig 2). The attack rate was the highest in larval mite treatment fed on by female A. eharai (0.1087), and while it was lowest, also for larval mites as prey but fed on by male A. eharai (0.0186). Handling time was shortest for larval mites fed on by female A. eharai (0.1642) and longest for deutonymphs fed on by female A. eharai (1.0232). For protonymphs fed on by female A. eharai the relationship was not statistically significant (P = 0.108). For deutonymphs, the attack rate and handling time of female A. eharai were 0.0224 and 1.2105. Finally, the P1 values for male predators attacking deutonymphs and female predators attacking adult mites were both not significant (P = 0.116 and P = 0.181, respectively). Moreover, even though the P1 value for male predators attacking adult prey was significant (P = 0.014), this treatment was excluded in the comparison of parameters and model fitting because of its higher standard error of parameters.
Discussion
In this study, we explored the predation traits of adult A. eharai on the different developmental stages of T. urticae, and, as far as we know, this is the first study on the functional response of A. eharai. Amblyseius eharai showed a type II functional response against eggs, larvae, and deutonymphs of T. urticae. Although the type of functional response can vary among prey species or the predator’s physiological condition, type II is a common response in predatory mites [18–22, 37].
Also, predatory mites can be divided into four dietary groups, with most generalist predators belonging to types III or IV [2, 3], preying on various small-sized pests but also ingesting pollen. Which prey and pollen are most suitable for particular species of type III or IV phytoseiid mites vary depending on predator species [2, 3]. In our study, both sexes of A. eharai preferred T. urticae larvae over eggs. Blackwood et al. [38] reported that type III generalists either preferred T. urticae larvae over eggs, or showed no preference between the two. Amblyseius eharai is a type III generalist [2, 3], with a preference for T. urticae larvae. Interestingly, during our predation tests, we noticed variation in body shape of the predatory mites depending on prey life stage they had consumed (Fig 3). While starved females or one that had consumed eggs seemed to have a similar appearance, with a translucent body color and narrower body width (Fig 3A and 3B), females that had consumed larvae had a brownish color for their digestive tack and broader body (Fig 3C). Gut emptying could be used as a physiological indicator of phytoseiids [39, 40]. In Fig 3B, its gut is empty despite the food consumption, thus the predator probably has poor nutritional condition. Based on our results, we conclude that T. urticae eggs are a poor quality food for A. eharai adults.
(A) Starved female for 24 hours. (B) Female that consumed T. urticae eggs for 24 hours. (C) Female that consumed T. urticae larvae for 24 hours.
Webbing of spider mites is known to serve as a shield for them, and it can mediate the anti-predator behavior [41]. In this study, any webbing was discarded from the experimental arena, thus having potential of affecting the outcome of functional response analysis in the study.
In our study, A. eharai females had higher predation rates than males, for all T. urticae life stages. Moreover, A. eharai males usually failed to consume T. urticae adults. Like other phytoseiid mites, the body size of female A. eharai (length of dorsal shield: 0.352 mm) is greater than that of males (length of dorsal shield: 0.266 mm) [7]. Higher predation by females than by males and the male’s failure to prey on T. urticae adults may be a simple function of body size.
The attack rates and handling times of male and female A. eharai were compared with literature values of other phytoseiids at 95% confidence level (Table 6). The attack rate of male A. eharai on T. urticae eggs was significantly lower than other predatory mites. The attack rate of female A. eharai on T. urticae eggs was not significantly different from that of A. swirskii (at 25°C) or N. californicus (at 25 and 30°C). The handling time of the predatory mites species included here on T. urticae eggs was the shortest for N. longispinosus at 35°C (0.264) and for female A. eharai (0.181). Female A. eharai showed a higher maximum attack rate (T/Th) on T. urticae eggs (132.74) than other phytoseiids, even higher than N. longispinosus at 35°C (90.91).
For larval T. urticae as prey, female A. eharai showed the highest attack rate (0.109) and the shortest handling time (0.164) among the compared phytoseiid mites. The attack rate of male A. eharai (0.019) was significantly lower than that of female A. eharai; however, there was no significance in handling time between the two sexes of this species. The maximum attack rates of both sexes of A. eharai (Female, 146.16; Male, 102.43) were higher than those of other phytoseiid mites considered.
Comparing A. eharai to other phytoseiids studied previously, A. eharai seems to have a higher predation capacity on T. urticae than A. swirskii, N. californicus, N. womersleyi, or P. persimilis, under similar conditions. Also, it has been suggested that the predation capacity of A. eharai might increase with adult age [20, 21]. Overall, A. eharai seems to have good potential for commercial development as a new natural enemy for augmentative biological control of T. urticae in greenhouses in mid-latitude region. In the studies we compared, there were variations in experimental layout such as plant species, size of leaf disc, mite age, and temperature. The test temperature for the referenced studies on T. urticae eggs or larvae was 25°C for all species of phytoseiids, and additionally 30 and 35°C were also tested for N. californicus, N. longispinosus, and N. womersleyi on T. urticae eggs as a prey (Table 6). However, the test temperature for A. eharai was 26.6°C in this study. The temperature would be one of the significant factors affecting the predation ability of predatory mites, as shown by the results for N. californicus, N. longispinosus, and N. womersleyi [18, 22]. These species showed higher predation efficiency at high temperatures. Nevertheless, A. eharai showed higher predation ability on T. urticae.
For T. urticae eggs, even though the attack rate of A. eharai was lowered compared to other predatory mites we considered, the handling time of adult females of A. eharai was significantly shorter than for A. swirskii (at 25°C), N. californicus (at 25, 30, and 35°C), N. longispinosus (at 25 and 30°C), N. womersleyi (at 25, 30, and 35°C), P. persimilis (at 25°C), and N. californicus (at 20 and 25°C). Also, in adult males of A. eharai, handling time for eggs as prey was shorter than for N. womersleyi (at 25, 30, and 35°C) and P. persimilis (at 25°C) (Table 6). To some degree these shorter handling times may compensate for the lower attack rates of A. eharai. Thus, the maximum attack rate (T/Th) of female A. eharai for T. urticae eggs was the highest among the compared phytoseiid mites, even higher than that of N. californicus, N. longispinosus, and N. womersleyi at 30 or 35°C.
Morphological structure on plant leaves such as density of trichomes can be a cause of variation in predation capacity of predatory mites [42], with some exception such as study by Ahn et al. [18] who found no difference in functional response of N. californicus against T. urticae between two varieties of strawberry with different leaf trichome densities. Also, age of phytoseiids can affect predation. Fathipour et al. [20] reported that the predatory capability of A. swirskii on T. urticae egg increased until 17 days after egg hatch of the predator at 25°C. Such a trend was also reported in P. persimilis [21]. The predation capacity of P. persimilis on T. urticae egg increased until this phytoseiid was 30 days old (post egg hatch) at 25°C. The phytoseiid predator age in our test was three days old after molting to adult (equivalent to 8 to 9 days old after egg hatch). Further research might be needed to determine whether the maximum predation capacity, as affected by maturation of the A. eharai individuals, was reached more quickly than by the other phytoseiids or if the predatory capacity of A. eharai was simply always higher even from an early age.
A temperature-dependent intrinsic rate of increase (r) model was constructed to predict the adaptability of A. eharai to temperature changes using the data of Park and Lee [11]. Also, to compare with A. eharai, the r models of other phytoseiid mites were extracted from referenced studies or constructed in this study using the data in references. For A. eharai and N. californicus, the r values at each temperature were fitted against temperatures using the Briere 1 model [43] to construct temperature-dependent r models by the program TableCurve 2D (SYSTAT Software Inc., San Jose, CA) [44]. And, for A. swirskii, N. womersleyi, and N. longispinosus, we used models from the original studies. The lower and upper thresholds and optimal temperatures for the population growth of all species were estimated and compared using the models. There was an error in the original research paper for the model of N. longispinosus [15]. Thus, it was reconstructed in this study using the same equation (Lactin 1) [45].
The lower- and upper-threshold, and optimal temperatures of selected phytoseiid mites from the literature were calculated using the temperature-dependent r model (Table 7). Among the compared mites, the optimal temperature for r was the lowest in A. eharai (28.1°C) followed by N. californicus (29.0°C), A. swirskii (30.1°C), N. womersleyi (30.5°C), and N. longispinosus (32.0°C). Also, A. eharai had the lowest upper-threshold temperature (33.2°C) and the narrowest temperature range (12.5 to 33.2°C) for population growth.
Many studies have examined temperature-dependent models of the developmental rate in immature stages (and for total fecundity) of various insects and mites [11, 43, 45, 47–49]. These models can be used to predict and better understand the population dynamics of insects and mites. However, the mites have a very short life cycle; therefore, an r value that included information for both the immature and adult stages would be a better parameter to use in constructing a temperature-dependent model for mites. In comparing temperature-dependent r models among phytoseiids, A. eharai showed a slightly narrower temperature range for development (12.5 to 33.2°C) than did the other species considered (Table 7). Moreover, the optimal temperature of A. eharai (28.1°C, Table 7) was the lowest among the other mites we compared. The functional response studies on N. californicus, N. longispinosus, and N. womersleyi to T. urticae eggs were conducted at three temperatures (25, 30, and 35°C) (Table 6). Neoseiulus californicus showed a higher predation rate at 35°C, and this rate did not correspond with the optimal temperature (29.0°C) for N. californicus in the temperature-dependent r model (Tables 6 and 7). However, the higher predation rate of N. womersleyi (30°C) and N. longispinosus (35°C) seemed to match their optimal temperatures, 30.5 and 32.0°C, respectively (Tables 6 and 7). The optimal temperature for population development by A. eharai was 28.1°C. This means that this predator can reach faster to its maximum performance at the season of low temperature than the other phytoseiids because they had a lower optimal temperature. Generally, the mean temperature in greenhouses in mid-latitude region does not exceed 30°C, even in the summer season [50]. Thus, A. eharai would have a high potential as a biological control agent in the greenhouses.
Conclusion
In conclusion, A. eharai, a native species in Far East Asia and the dominant phytoseiid in crops in Korea [7, 16], has a high potential as a new biological control agent because of its higher predation capacity. Moreover, as this predatory mite has a relatively lower optimal temperature for population growth, it would confer an advantage for the use at the season of low temperature. Lastly, further studies on additional topics, such as this predator’s response to other pests or an outdoor field study, or a greenhouse crop study, are needed to advance our understanding of this species’ potential for commercial use.
Acknowledgments
We are grateful to Mrs. Jung Ran Choi for helping rear two-spotted spider mites and grow kidney bean plants.
References
- 1. Chant DA. Phytoseiid mites (Acarina: Phytoseiidae). Mem Ent Soc Can. 1959; 91(S12): 5–166.
- 2. McMurtry JA, de Moraes GJ, Sourassou NF. Revision of the lifestyles of phytoseiid mites (Acari: Phytoseiidae) and implications for biological control strategies. Sys Appl Acarol. 2013; 18(4): 297–320.
- 3. McMurtry JA, Croft BA. Life-style of phytoseiid mites and their roles in biological control. Annu Rev Entomol. 1997; 42: 291–321. pmid:15012316
- 4. Collier T, Van Steenwyk R. A critical evaluation of augmentative biological control. Biol Control. 2004; 31(2): 245–256.
- 5.
McMurtry JA, Sourassou NF, Demite PR. The Phytoseiidae (Acari: Mesostigmata) as biological control agents. In: Carrillo D, de Moraes GJ, Peña JE, editors. Prospects for biological control of plant feeding mites and other harmful organisms. Springer, Cham; 2015. pp. 133–149.
- 6. Lee MH, Zhang ZQ. Assessing the augmentation of Amblydromalus limonicus with the supplementation of pollen, thread, and substrates to combat greenhouse whitefly populations. Sci Rep. 2018; 8: 12189. pmid:30111848
- 7. Amitai S, Swirski E. A new species of Amblyseius (Acarina: Phytoseiidae) for the Far East. Isr J Entomol. 1981; 15: 59–66.
- 8. Waite GK, Gerson U. The predator guild associated with Aceria litchii (Acari: Eriophyidae) in Australia and China. Entomophaga. 1994; 39(3–4): 275–280.
- 9. Kakimoto K, Iguchi T, Inoue H, Kusigemati K. Predatory ability of Amblyseius eharai Amitai et Swirski (Acari: Phytoseiidae) on thrips. Kyushu Plant Prot Res. 2004; 50: 82–87.
- 10. Ji J, Zhang YX, Lin JZ, Chen X, Saito Y. Life histories of three predatory mites feeding upon Carpoglyphus lactis (Acari, Phytoseiidae; Carpoglyphidae). Syst Appl Acarol. 2013; 20(5): 491–496.
- 11. Park Y, Lee JH. Temperature-dependent development and oviposition models and life history characteristics of Amblyseius eharai (Amitai et Swirski) (Acari: Phytoseiidae) preying on Tetranychus urticae (Koch) (Acari: Tetranychidae). J Asia-Pac Entomol. 2020; 23(4): 869–878.
- 12. Kim T, Ahn J, Lee JH. Temperature-dependent developmental model of Neoseiulus californicus (McGregor) (Acari, Phytoseiidae). J Appl Entomol. 2009; 133: 284–291.
- 13. Lee HS, Gillespie DR. Life tables and development of Amblyseius swirskii (Acari: Phytoseiidae) at different temperatures. Exp Appl Acarol. 2011; 53: 17–27. pmid:20628894
- 14. Kim T, Ahn JJ, Lee JH. Age- and temperature-dependent oviposition model of Neoseiulus californicus (McGregor) (Acari: Phytoseiidae) with Tetranychus urticae as prey. J Appl Entomol. 2013; 137: 282–288.
- 15. Sugawara R, Ullah MS, Ho CC, Gökçe A, Chi H, Gotoh T. Temperature-dependent demography of two closely related predatory mites Neoseiulus womersleyi and N. longispinosus (Acari: Phytoseiidae). J Econ Entomol. 2017; 110(4): 1533–1546. pmid:28854663
- 16. Ryu MO, Lee WK, Kim TH. Habitats and abundances of Korean phytoseiid mites. Korean J Appl Entomol. 1997; 36(3): 224–230.
- 17. Cuthbert RN, Dick JT, Callaghan A, Dickey JW. Biological control agent selection under environmental change using functional responses, abundances and fecundities; the Relative Control Potential (RCP) metric. Biol Control. 2018; 121: 50–57.
- 18. Ahn JJ, Kim KW, Lee JH. Functional response of Neoseiulus californicus (Acari: Phytoseiidae) to Tetranychus urticae (Acari: Tetranychidae) on strawberry leaves. J Appl Entomol. 2010; 134(2): 98–104.
- 19. Song ZW, Zheng Y, Zhang BX, Li DS. Prey consumption and functional response of Neoseiulus californicus and Neoseiulus longispinosus (Acari: Phytoseiidae) on Tetranychus urticae and Tetranychus kanzawai (Acari: Tetranychidae). Syst Appl Acarol. 2016; 21(7): 936–946.
- 20. Fathipour Y, Karimi M, Farazmand A, Talebi A A. Age-specific functional response and predation rate of Amblyseius swirskii (Phytoseiidae) on two-spotted spider mite. Syst Appl Acarol. 2017; 22(2): 159–169.
- 21. Fathipour Y, Karimi M, Farazmand A, Ali AT. Age-specific functional response and predation capacity of Phytoseiulus persimilis (Phytoseiidae) on the two-spotted spider mite. Acarol. 2017; 58(1): 31–40.
- 22. Sugawara R, Ullah MS, Ho CC, Gotoh T. Impact of temperature-mediated functional responses of Neoseiulus womersleyi and N. longispinosus (Acari: Phytoseiidae) on Tetranychus urticae (Acari: Tetranychidae). Biol Control. 2018; 126: 26–35.
- 23. Fathipour Y, Maleknia B, Bagheri A, Soufbaf M, Reddy GV. Functional and numerical responses, mutual interference, and resource switching of Amblyseius swirskii on two-spotted spider mite. Biol Control. 2020; 146: 104266.
- 24. Park YL, Lee JH. Leaf cell and tissue damage of cucumber caused by twospotted spider mite (Acari: Tetranychidae). J Econ Entomol. 2002; 95(5): 952–957. pmid:12403421
- 25. Pascual-Ruiz S, Aguilar-Fenollosa E, Ibáñez-Gual V, Hurtado-Ruiz MA, Martínez-Ferrer MT, Jacas JA. Economic threshold for Tetranychus urticae (Acari: Tetranychidae) in clementine mandarins Citrus clementina. Exp Appl Acarol. 2014; 62(3): 337–362. pmid:24122147
- 26. Sousa VC, Zélé F, Rodrigues LR, Godinho DP, de la Masselière MC, Magalhães S. Rapid host-plant adaptation in the herbivorous spider mite Tetranychus urticae occurs at low cost. Curr Opin Insect Sci. 2019; 36: 82–89. pmid:31539789
- 27. van den Boom CEM, van Beek TA, Dicke M. Differences among plant species in acceptance by the spider mite Tetranychus urticae Koch. J Appl Entomol. 2003; 127(3): 177–183.
- 28. Marinosci C, Magalhaes S, Macke E, Navajas M, Carbonell D, Devaux C, et al. Effects of host plant on life-history traits in the polyphagous spider mite Tetranychus urticae. Ecol Evol. 2015; 5(15): 3151–3158. pmid:26356681
- 29. Dermauw W, Wybouw N, Rombauts S, Menten B, Vontas J, Grbić M, et al. A link between host plant adaptation and pesticide resistance in the polyphagous spider mite Tetranychus urticae. PNAS. 2013; 110(2): E113–E122. pmid:23248300
- 30. Kwon DH, Kang TJ, Kim YH, Lee SH. Phenotypic-and genotypic-resistance detection for adaptive resistance management in Tetranychus urticae Koch. PLoS ONE. 2015; 10(11): e0139934. pmid:26545209
- 31. Attia S, Grissa KL, Lognay G, Bitume E, Hance T, Mailleux AC. A review of the major biological approaches to control the worldwide pest Tetranychus urticae (Acari: Tetranychidae) with special reference to natural pesticides. J Pest Sci. 2013; 86(3): 361–386.
- 32. Rahman MA, Sarker S, Ham E, Lee JS, Lim UT. Development and fecundity of Orius minutus (Hemiptera: Anthocoridae) and O. laevigatus reared on Tetranychus urticae (Acari: Tetranychidae). J Econ Entomol. 2020; 113(4): 1735–1740. pmid:32333019
- 33.
Hajek AE, Eilenberg J. Natural enemies: an introduction to biological control. Cambridge, UK: Cambridge University Press; 2018.
- 34.
SAS Institute. SAS user’s guide; statistics, version 9.4ed. Cary. NC: SAS Institute; 2013.
- 35.
Juliano SA. Nonlinear curve fitting: predation and functional response curves. In: Scheiner SM, Gurevitch J, editors. Design and analysis of ecological experiments. Chapman & Hall, New York; 2001. pp. 178–196.
- 36. Rogers D. Random search and insect population models. J Anim Ecol. 1972; 41(2): 369–383.
- 37. Papanikolaou NE, Broufas GD, Papachristos DP, Pappas ML, Kyriakaki C, Samaras K, et al. On the mechanistic understanding of predator feeding behavior using the functional response concept. Ecosphere. 2020; 11(5): e03147.
- 38. Blackwood JS, Schausberger P, Croft BA. Prey-stage preference in generalist and specialist phytoseiid mites (Acari: Phytoseiidae) when offered Tetranychus urticae (Acari: Tetranychidae) eggs and larvae. Environ Entomol, 2001; 30(6): 1103–1111.
- 39.
Sabelis MW. The functional response of predatory mites to the density of two-spotted spider mites. In: Metz JAJ, Diekmann O, editors. Dynamics of physiologically structured populations. Springer, Berlin Heidelberg New York; 1986. pp. 298–328.
- 40. Dicke M, Sabelis MW, van den Berg H. Does prey preference change as a result of prey species being presented together? Analysis of prey selection by the predatory mite Typhlodromus pyri (Acarina: Phytoseiidae). Oecologia. 1989; 81(3): 302–309. pmid:28311180
- 41. Lemos F, Sarmento RA, Pallini A, Dias CR, Sabelis MW, Janssen A. Spider mite web mediates anti-predator behaviour. Exp Appl Acarol. 2010; 52(1): 1–10. pmid:20191311
- 42. Schmidt RA. Leaf structures affect predatory mites (Acari: Phytoseiidae) and biological control: a review. Exp Appl Acarol. 2014; 62(1): 1–17. pmid:23990040
- 43. Briere JF, Pracros P, Le Roux AY, Pierre JS. A novel rate model of temperature-dependent development for arthropods. Environ Entomol. 1999; 28(1): 22–29.
- 44.
SYSTAT Software Inc. TableCurve 2D; Automated curve fitting and equation discovery, version 5.01ed. San Jose. CA: SYSTAT Software Inc.; 2002.
- 45. Lactin DJ, Holliday NJ, Johnson DL, Craigen R. Improved rate model of temperature-dependent development by arthropods. Environ Entomol. 1995; 24(1): 68–75.
- 46. El Taj HF, Jung C. Effect of temperature on the life-history traits of Neoseiulus californicus (Acari: Phytoseiidae) fed on Panonychus ulmi. Exp Appl Acarol. 2012; 56: 247–260. pmid:22270114
- 47. Logan JA, Wollkind DJ, Hoyt SC, Tanigoshi LK. An analytic model for description of temperature dependent rate phenomena in arthropods. Environ Entomol. 1976; 5(6): 1133–1140.
- 48. Kim DS, Lee JH. Oviposition model of overwintered adult Tetranychus urticae (Acari: Tetranychidae) and mite phenology on the ground cover in apple orchards. Exp Appl Acarol. 2003; 31(3): 191–208.
- 49. Shi PJ, Reddy GV, Chen L, Ge F. Comparison of thermal performance equations in describing temperature-dependent developmental rates of insects: (I) empirical models. Ann Entomol Soc Am. 2016; 109(2): 211–215.
- 50.
Park Y. Eco-friendly pest management strategy for Frankliniella occidentalis and Bemisia tabaci in tomato smart greenhouse of Korea. Ph.D. Dissertation, Seoul National University. 2020. https://s-space.snu.ac.kr/bitstream/10371/169668/1/000000162178.pdf