Introduction
Tetranychus urticae Koch (Acari: Tetranychidae) is a cosmopolitan and polyphagous pest that attacks many economic importance crops in Chile, such as vegetables, extensive cultures (cotton, corn, etc.), citrus, vine and other fruit and ornamental trees. Although different strategies for control of T. urticae are recommended worldwide, the use of natural enemies has been highlighted as a tool of first order for biological control of this pest, especially the use of predator mite Phytoseiulus persimilis Athias-Henriot (Acari: Phytoseiidae) (De Moraes, 2004). However, this practice has been less common that the use of chemical acaricides, which improperly used, can cause resurgence, resistance and other problems of pests due to the elimination of natural enemies (Vitelli & Sato, 2016).
Cydnodromus picanus, was described in 2000 on Citrus aurantium L. (Rutaceae) and reported only for Chile (Ragusa, 2000; De Moraes et al., 2004). Recent researches suggest C. picanus, could be considered as junior synonymy of Cydnodromus idaeus (Denmark and Muma) (Tixier et al., 2011). According to Ragusa et al. (2000), this species could be used in arid and warm agroecosystems, where is difficult to find phytoseiids that are adapted to these extreme conditions to control tetranychids. Based on classification of McMurtry & Croft (1997), C. picanus is a generalist predator belonging to type III lifestyle. In conditions of absence of prey can feed of pollen from different species (Ragusa et al. 2000), depredates on Oligonychus yothersi (McGregor) in Persea americana Mill. (Lauraceae) (Rioja & Vargas, 2009) and on Tetranychus cinnabarinus Boisduval in Dianthus caryophyllus Linn. (Caryophyllaceae) (Tello et al., 2009). The citrus are one of the more important crops in the zone of Pica, Chile (20º 29`12.4" S; 69º 19` 33.9" O), where C. picanus had been found associated to Panonychus citri (McGregor) (Ragusa et al., 2000).
The present research provides information about the predatory efficiency of C. picanus throughout functional response. The equations of the functional response are a key component of any mechanistic model, because they describe the interactions among consumers and its resources (Villemereuil & López-Sepulcre, 2011). This information is vital for the implementation of conservation programs for natural enemies and augmentation-type biological control (Ragusa et al., 2010). The effectivity of a predator is directly related to the type of its functional response. Predatory arthropods unfold one of the three typical functional responses, but these responses may vary with the crop phenology, the habitat heterogeneity, the predator age and other biotic and abiotic factors. In type I functional response, there is a lineal increase of predator attack rate respect to prey density (with a slope equal to the search efficiency), until reaching a point at which the maximum attack rate remains constant, even if the prey density increases. The response is estimated by a linear equation Na=α+βN0 ; where Na =number of prey eaten, N0 =prey density (number of prey offered), α and β are the y-intercept and the slope of prediction line, respectively. In type II functional response appears another parameter, in addition to the search efficiency, termed "handling time" (Th ) defined as the time that invests the predator in chasing, capturing, consuming and digesting a prey. In type II functional response, the number of prey consumed hyperbolically approaches to an asymptote in which the maximum attack rate is expressed (K=T/Th ). This type of response is estimated by a curved line function: Na =aN0T/(1+aN0Th), where Na =number of prey consumed, a=constant attack rate or instantaneous search rate, N0 =initial prey density, T=total available time and Th =handling time. The type III functional response is performed when the prey-consumed number approaches an asymptote as a sigmoid function (increases in the proportion of prey killed to an inflection point and then decreases similarly to type II response). The type III functional response may be derived using the model of Hassell (Hassell, 1978), where the attack rate (a) is a hyperbolic function of prey density: a=(d+hN0)/(1+cN0 ), substituting the value of the Holling equation, obtaining the type III model: Na=(dN0T+hN0 2T )/(1+cN0+dN0Th+bN0 2T h ), where Na =number of prey consumed, N0 =initial prey density, T=total time, Th =handling time and b, c, and d are constants.
Tello et al. (2009), proved the life parameters of C. picanus when fed with eggs, the same happened with survival, achieving 100% when fed with this stadium. For those experiments, eggs less than 20 h old were performed. In addition, the functional response of this phytoseiid was affected by the maturity of the egg offered as food.
The aim of this research was to evaluate the consumption rates of C. picanus (=C. idaeus) in relation to eggs density of T. urticae and to determine the effect of egg maturity on functional response parameters, under laboratory conditions.
Material and methods
Establishment of colonies
The bioassays were performed in the Regional Research Center INIA La Cruz (Instituto de Investigaciones Agropecuarias (INIA), La Cruz, Chile), at a temperature of 25±2°C, 50±2% of relative humidity (RH) and a photoperiod of 16:8 h (L:D). Cydnodromus picanus, was collected from trees of C. aurantium in Pica, Chile (20°29' S; 69°19' O) and bred in acrylic plates of 12x12 cm according to the methodology described by Swirski, Amitai & Dorzia (1970), in breeding rooms and fed only with pollen from Oxalis sp (Oxidales) to avoid habituation to prey. Tetranychus urticae obtained from breeding rooms located in INIA La Cruz. The colonies of this tetranychid had a time of about four months since they were collected in the field. Tetranychus urticae was reproduced on bean plants (Phaseolus vulgaris L.), that were infested with the mite when they had two true leaves. To maintain available plants, 100 pots weekly seeded with 20 bean seeds.
Consumption rate and functional response
The functional response of C. picanus on T. urticae eggs of different maturity time studied in separated bioassays. For each maturity time, one phytoseiid adult female, previously subjected to 24 h fasting, was confined in an artificial arena (black acrylic plate of 4 cm diameter) (Tello et al., 2009). These C. picanus females came from a cohort of eggs of the same age, using newly emerged females. Glue (Point Sticken Glue, Point Chile S.A.) was applied in the edges plates to prevent the mites escape. Therefore, T. urticae eggs of different maturity times (16, 42 and 65 h) were introduced as preys in the plates at densities of 5, 10, 20, 30, 40, 50, 60, 70, 80, 100 and 150 eggs per plate. T. urticae eggs were obtained by depositing on a bean leaf about 100 females to oviposit and at 16h were transferred into trial arenas, replicating this for eggs of 42 and 65 h. The number of T. urticae eggs consumed by the predator recorded at 24 h (experimental time). The observations were made using a stereomicroscope with 40X magnification (Zeiss Stemi 6V 6, Alemania).
Statistical analysis
The consumption among eggs with different maturity (16, 42 and 65 h) were compared. In addition, for each age, among different densities. Because consumption data is not normally distributed, even after applying transformations (e.g. logarithmic and square root). The non-parametric Kruskal-Wallis test with post-hoc comparison were performed to determine the type of functional response, data of each density were fitted to maximum likelihood logistic regression (CATMOD procedure, SAS Institute, Cary, NC, USA). (Equation 1).
Where Y, is the dichotomous dependent variable, representing a non-consumed (Y=0) or consumed (Y=1) prey at the end of the trial, Na is the number of prey consumed, No is the number of prey available; and P0, P1, P2 and P3 are the intercept, linear, quadratic and cubic coefficient, respectively. Negative or positive significant linear coefficients (e.g., P1 ) indicate a type II or III functional response, respectively. The slope of type II response decreases (e.g., P1 is negative), while that slope of type III increases (e.g., P1 is positive). After determining the type of functional response, its parameters (attack rate and handling time) are estimated through linearization of Holling' s disc equation depending on the type of functional response (Xiao & Fadamiro, 2010). Type II equation was used: a and Th parameters were calculated using a nonlinear regression by least squares based on Gauss-Newton method (SAS Institute, Cary, NC, USA). The a/Th value indicates the effectiveness of depredation. The maximum depredation rate (K) was calculated as T/Th . To perform comparisons between a and Th parameters, standard deviation was estimated with 95% confidence interval using Jackknife statistical technique available in SAS (SAS Institute, Cary, NC, USA). Later, the parameters were compared using non-parametric Kruskal-Wallis test (α=0.05) due to heterogeneity of variances and non-normally distributed residuals.
Results and discussion
Consumption rate
The average consumption rate of C. picanus did not show significant differences among the three types of eggs offered [KW=4.04; gl=2; P=0.131]. This indicates that the maturity time of eggs offered as prey have no influence in the preference by C. picanus. The predator consumed 54, 52 and 51% of initial egg density at 16, 42 and 65 h, respectively, not showing significant differences [KW =1.166; gl=2; P=0.556] (Table 1).
Equal letters in the same column indicate no significant differences according to Kruskal-Wallis test (P<0.05); Na=number of prey consumed; No=initial prey density; SE=standard error.
From a density of 50 eggs.predator-1, the age of eggs did not affect the predation by C. picanus; 5 (KW=0.32, gl=2, P=0,3679), 10 (KW=2.33, gl=2, P<0.1503), 20 KW=0.43, gl=2, P=0.7884), 30 (KW=4.76, gl=2, P<0.0913), 40 (KW=0.04, gl=2, P=0.9807) and 50 (KW=1.35, gl=2, P=0.5053). The Kruskal-Wallis test did not show significantly differences (P>0.05) for those prey densities (Table 2). Differences showed densities of 60 (KW=12.49, gl=2, P<0.05), 70 (KW=8.75, gl=2, P<0.05) and 100 (KW=8.54, gl=2; P<0.05) eggs.predator-1. The number of attacked eggs per capita was significantly affected by prey density for each egg maturity (16h: KW=54.87, gl=10, P< 0.0001), 42h: KW=43.23, gl=10, P<0.0001 and 65h: KW=46.31, gl=10, P<0.0001).
Functional response
The logistic regression delivered a significantly negative linear parameter (P1 <0) for the three classes of prey offered (eggs of 16, 42 and 65 h), suggesting that C. picanus performs a type II functional response for the three types of prey (Table 3). The functional response curves (Figure 1A), showed a high proportion of prey consumed at lower densities in the three kinds of prey. The proportion of eggs consumed by predator (Na/N0 ) declined with the increase of prey density (Figure 1B); this implies that there is an inverse density-dependence, being confirmed by logistic regression (Table 3) that, for all ages of eggs, presented a negative linear coefficient (P1 <0) and a positive quadratic coefficient (P2 >0), both significant (P<0.001; Table 3).
The comparison of functional response curves revealed significant differences among the three kinds of eggs offered (a: KW=11.94, gl=2, P<0.05; Th : KW=11.51, gl=2, P<0.05). The maximum number of prey consumed was significantly higher for eggs of maturity of 16 h (Table 4). The estimated parameters of functional response for C. picanus, for the three types of prey, showed that the attack rate was higher for eggs of 42 h and 65 h than for eggs of 16 h. On the other hand, C. picanus showed the lowest handling time (Th ) for eggs of 16 h (Table 4). The slightly higher a/Th values recorded for C. picanus on eggs of 42 and 65 h suggest that this predator was more effective on consuming that sort of T. urticae eggs.
Values with equal letters in the same column are significantly not different according to Kruskal-Wallis test (P<0.05); a=search or attack rate; Th =handling time (day); R2 =determination coefficient; a/Th =depredation efficiency; K=maximum depredation rate.
The average eggs consumption rate is similar to those obtained for adult females of this phytoseiid species. Ragusa et al. (2000) indicates a daily consumption of about 16 eggs of T. urticae at 26±1°C and 70±5% HR. Tello et al. (2009) registered the same consumption value (16 eggs.day-1) with T. cinnabarinus eggs at 29.44±1.47°C and 42.35±5.01% RH.
Cydnodromus picanus, showed a type II functional response for each maturity time of eggs. This type of functional response is registered for different phytoseiid especies, including: Amblyseius californicus (Mc-Gregor) (Gotoh, Nozawa & Yamaguchi, 2004); Neoseiulus longispinosus (Evans) (Rahman, Babu, Roobakkumar & Perumalsamy, 2012); Phytoseiulus persimilis Athias-Henriot, Galendromus occidentalis (Nesbitt) (Xiao & Fadamiro, 2010); Neoseiulus cucumeris (Oudemans) (Madadi et al., 2007). In general, the results obtained in this study for a and Th parameters of functional response for C. picanus are within the range of variation for different phytoseiid especies (Table 5).
The search (or attack) rate and the handling time are parameters used to determine the magnitude of the functional responses. This parameter could be affected by different factors, such as the predator velocity, the prey movements and the spent time in subduing each prey (Hassell, 1978). In this research, the handling time for eggs of 16 h was significantly lower than for eggs of 42 and 65 h, which could be explained for the presence of a higher amount of vitellus and the easiness to extract it. When eggs mature, the embryonic development makes diminish the amount of vitellus, making difficult its extraction and making to take longer time to predator in feeding. Cabrera, Donohuea, Khalil, Sonenshinec & Roea (2009), pointed that there is no available studies describing the characteristics of proteins in the yolk of mite eggs. According to these authors, in newly oviposited tetranychids eggs (<3 h), the embryo development would be minimal and the proteins would be more available for predators. Therefore, we can deduce that as time progress and egg maturity with it; the eggs of 42 and 65 h present higher embryo development, so the amount of vitellus must be minimal and thus of lower nutritional value.
Conclusion
This research compared the average consumption rates and the functional response of predator mite C. picanus fed with T. urticae eggs of different maturity times. The results show that adult females of C. picanus exhibita type II functional response and a proportion of prey consumed that ranged between 51% and 54% for the three types of prey offered. Studies of predation in the field, with different crops and preys (tetranychid species), are necessary to consider this predator as a potential biological control agent in arid zones of Chile.