Introduction

Spodoptera litura (Fabricius) (Lepidoptera: Noctuidae) is a cosmopolitan and obnoxious pest of more than 100 crops, ornamentals, and vegetable species1. A single mother can lay over 2,000 eggs and huge populations are present when conditions are favourable2. Larvae are gregarious leaf eaters, and uncontrolled feeding can induce 100 percent losses in some sensitive crops. Pesticides are primary tools that are used to control this pest, but inappropriate use has several associated problems. It has led to the development of pesticide resistance against routinely used organophosphates, carbamates, pyrethroids, indoxacarb, abamectin, emamectin benzoate, and chlorantraniliprole3. As a result of resistance development, it becomes difficult to control the pest and outbreaks occur due to control failures4. Further challenges include potential risks to biodiversity, humans, and the environment. Environmental-friendly strategies such as biological control as replacements of conventionally used pesticides are thus need of time5,6,7.

Biological control is one of the finest means of suppressing pest problems with reduced reliance on synthetic insecticides8. Many factors that significantly determine biocontrol fate include species identity, resource quality/quantity, and other ecological considerations9. Temperature is primarily important for biological systems through its close link with species metabolic rates and fitness parameters i.e., survival, development, reproduction, handling time, and searching efficiency of predators10. Since insects are ectotherms, they rely greatly on temperature for their growth, reproduction, maturation and trophic relations11. Insects vary greatly in thermal sensitivity levels, hunting behaviour, and resource-specificity. Assessment of the effect of temperature at individual species level12 is typically desired to improve our understanding about temperature regulation of trophic cascades and population dynamics13. This will pave a way towards the development of effective biocontrol strategies under altered climates11. The use of age-stage, two-sex life table approach, developed by Chi and Liu (1985) and Chi (1988), represents an effective way of studying fitness of the individuals under influence of stressors14.

Polyphagous predators can be mass-reared and produced with artificial diets, and alternative resource can increase their numbers in the field15. The family Coccinellidae includes a globally widespread group of predatory beetles that are polyphagous16 and indicate preference for resource including transitions across kingdoms (Animalia, Fungi, and Plantae) and trophic levels (carnivorous and herbivorous)17. Many coccinellids are important predators of hemipterans, including aphids10,18, whitefly19, mealybugs20, scale insects21 and immature lepidopterans22,23. Resource quality, thermal surroundings, mobility24, nutritional status25, as well as density and size of the prey26 are critical determinants of predator development, biology16 and predation parameters12,27,28,29. The biological control effectiveness of the predator can thus be evaluated by combining development, survival rates, and reproductive potential with age-specific predation rates30.

The multicolored Asian ladybird beetle, Harmonia axyridis (Pallas) (Coleoptera: Coccinellidae), native to Asia, is a globally widespread predator that feeds on many resources, including 77 insect taxa on 85 plant species in 35 families15,31. This predator is highly voracious, multi-voltine, greatly dispersed32 and has ability to survive complex tropical and sub-tropical environments around the world31. H. axyridis has been commercially reared and extensively deployed as a potential biocontrol agent in many biocontrol programs worldwide. Its feeding is reported on hemipterans and lepidopterans, including immature stages of Danaus plexippus (Lepidoptera: Nymphalidae)22, Anagasta kuehniella (Lepidoptera: Pyralidae)33, Sitotroga cerealella (Lepidoptera: Gelechiidae)34, Spodoptera litura35 and Spodoptera frugiperda (Lepidoptera: Noctuidae)36.

Life table of this predator has been developed on various insects and temperatures10,33,34. The population parameters and feeding rate of this predator have been shown to depend on the temperature and resource30. Using laboratory experiments, we showed that foraging by H. axyridis on the eggs and larvae of S. litura (eggs and larvae) is temperature-derived23,35. Predator (larvae and adults) increases its foraging with warmer temperatures and decreases with colder temperatures. What remains unexplored is the impact of temperature-influenced feeding on subsequent predator development and predation rate. The present study answers this question by feeding and rearing H. axyridis on the eggs of S. litura at four constant temperatures by hypothesizing that predation and life history parameters of H. axyridis on S. litura eggs will be different across four temperatures tested.

Materials and methods

Mass rearing of insects

Harmonia axyridis population culture was established from a laboratory reared stock at Key Laboratory of Huazhong Agricultural University (HZAU), Wuhan, China, during February 2019. The adult insects were grouped in transparent plastic boxes (12 cm × 16 cm × 7 cm) with moistened cotton swabs, and reared on ad libitum diet of Acyrthosiphon pisum nymphs Harris (Hemiptera: Aphididae). A. pisum was reared on faba bean Vicia faba L. plants under conditions of 22 ± 1 °C, 65 ± 5% RH and 16:8 (L:D) h photoperiod. The seeds of faba bean were purchased online and grown inside cages. A. pisum was then introduced to these faba beans to start an aphid laboratory culture. Eggs of H. axyridis were removed daily and kept in 9 cm Petri dishes on moistened tissue paper. Hatchlings were transferred to separate Petri dishes to avoid cannibalism. These Petri dishes had folded tissue paper, A. pisum nymphs as food and moist cotton as water source. All the insects were kept under controlled laboratory conditions at 70 ± 5% RH, 24 ± 1 °C, and 16:8 (L:D) h photoperiod. All larval, pupal, and adult stages were maintained under same settings as above. The plant materials used were obtained with prior permission, and the present study complies with relevant guidelines and legislation.

The pupae of S. litura were purchased online from Henan Jiyuan Baiyun Industrial Co., Ltd, China, in September 2018 and placed in a transparent glass jar (20 cm × 12 cm × 8 cm) and incubated inside the growth chamber (Jiangnan Instruments RXZ-430, Ningbo, China) at settings of 26 ± 1 °C and 70 ± 5% RH under a 16:8 h photoperiod (L: D) to obtain adult populations. Hatchlings from S. litura eggs were fed with an artificial diet and cultured until the reach of final instar stage inside a mesh covered transparent circular glass jars (1 L). The artificial diet prepared followed Saljoqi, et al.37. The artificial diet (semi-solid) consisted of yeast powder (24 g), ascorbic acid (2.35 g), methyl-4-hydroxy benzoate (1.5 g), distilled water 550 mL, kidney bean flour (150 g), agar (8.4 g), sorbic acid (0.75 g), formaldehyde solution (1 mL) and streptomycin (0.75 g). The final caterpillars were isolated and subsequently cultured until maturity inside similar jars with 20% glucose solution as adult’s diet and clean pieces of paper covering jar’s bottom as oviposition substratum. The cotton wool strips (2 cm wide, 7–12 cm long) were placed in jars as a suitable oviposition substrate in order to obtain new oviposition. Oviposition was checked daily. Hatchlings were reared together in a separate Petri dish. Same age larvae were reared together using the same artificial diet as above to maintain population homogeneity.

Life table and predation trials

The experimental setup used here was a Petri dish (6 cm), and one leaf disc of cotton (Gossypium hirsutum L. cv. Varamin) placed on a 0.5 cm thick layer of 1.5% agar to fix the leaf disk. Paired H. axyridis adults were separated from the stock culture and reared with S. litura eggs in transparent plastic cups (9 cm in diameter and 7 cm in height) covered with fine nylon netting for ventilation under conditions of 70 ± 5% RH, 24 ± 1 °C, and 16:8 (L:D) h photophase for two consecutive generations to avoid interbreeding. Hundred H. axyridis eggs laid within 24 h were collected and transferred to Petri dish for the life table trials. Hatchlings were transferred to new Petri dishes and reared on sufficient eggs until pupation. The duration and viability of the individuals were recorded daily. Instar changes were evaluated based on the exuvia presence. The experiment was performed at four different temperatures (i.e., 15, 20, 25 and 30 °C), 70 ± 5% RH, and 16:8 (Light: Dark) h photoperiod. The temperature regimes tested here represent the thermal conditions that H. axyridis encounters in different protected plantations and field crops in tropical and temperate regions38. Experimental cohorts were generated as fifty individuals per temperature. Survival rate and time to hatch were recorded for each stage in each temperature treatment. The number of surviving larvae and the total developmental period from larva to pupa were also recorded. Predation was recorded for both larval and adult stages of the predator. After eclosion, the predator larvae were fed with S. litura eggs until pupation. In addition to this, male and female adults that emerged from the pupae on the same day from a given temperature were numbered and coupled in the insect boxes. Adult pairs were monitored to ensure that mating occurred, and time allowed for mating event was ~ 10 h a day, and no prey was made available during mating events throughout the trial duration. One or a few matings are reported sufficient for coccinellids females to fertilize all eggs. No eggs provision during mating events was meant to avoid any potential disturbance during this activity. After mating, the pairs were separated and shifted cautiously to their respective Petri dishes with eggs. The life span of an adult male is shorter than female, hence extra pairs of insects were raised in parallel to supply supplementary males for mating depending on the need at each temperature, simultaneously. The number of eggs laid per female including the cannibalized ones and adult mortality was noted twice daily until all female adults died. About 400–500 eggs of S. litura were supplied daily to the larvae and adults on cotton leaves. The counted numbers of S. litura eggs were fed to each stage of the predator. Predation data were recorded daily as numbers of eggsdamaged or consumed. Petri dishes and boxes were cleaned daily with 75% ethanol.

Data analysis

The data were analyzed by using TWOSEX-MSChart39 according to age-stage, two-sex life table theory14,40. The age stage-specific survival rate (sxj), age-specific survival rate (lx), age-specific fecundity (mx), the net maternity (lxmx), age-stage life expectancy (exj), age-stage reproductive value (vxj), population and reproduction parameters were calculated accordingly. All oviposited eggs per female were counted to calculate the fecundity.

The age-stage survival rate sxj was calculated as:

$$s_{xj} = \frac{{n_{xj} }}{{n_{01} }} \,$$
(1)

where n01 is the number of individuals used at the beginning of the life table study and nxj is the number of individuals surviving to age x and stage j41.

The age-specific survival rate (lx) was calculated as:

$$l_{x} = \sum\limits_{j = 1}^{\beta } {s_{xj} } \,$$
(2)

and age-specific fecundity (mx) were obtained from the following formula, where β is the number of stage sxj is the age-stage specific survival rate, i.e. the possibility that a new insect (newly hatched individual) will live or exist to age x and stage j.

$$m_{x} = \frac{{\sum\limits_{j = 1}^{\beta } {s_{xj} f_{xj} } }}{{\sum\limits_{j = 1}^{\beta } {s_{xj} } }}$$
(3)

The net reproductive rate (R0) is defined as the total number of offspring that an insect or individual can produce during its lifetime and was calculated as:

$$R_{{0}} = \sum\limits_{x = 0}^{\infty } {l_{x} m_{x} }$$
(4)

The intrinsic rate of increase (r) was then predicted iteratively following the Euler–Lotka equation with age indexed from 042 by the following formula,

$$\sum\limits_{x = 0}^{\infty } {e^{{ - r\left( {x + 1} \right)}} l_{x} m_{x} } = 1$$
(5)

The finite rate of increase (λ) was calculated as:

$$\lambda = e^{r}$$
(6)

The mean generation time (T) is the time duration that a population needs to increase to \({\text{R}}_{{0}}\)-fold of its size at the stable age-stage distribution, and was calculated as:

$$T = \frac{{\ln R_{0} }}{r}$$
(7)

Age-stage life expectancy (exj) is the duration that an insect of x age and j stage is predicted to live43 as:

$$e_{xj} = \sum\limits_{i = x}^{\infty } {\sum\limits_{y = j}^{\beta } {s^{\prime}_{iy} } }$$
(8)

where \(s^{\prime}_{iy}\) is the probability that an individual of age x and stage j will survive to age i and stage y and it is calculated by assuming sxj = 1.

The age-stage reproductive value (vxj) was defined as the contribution of individuals of age x and stage j to the future population and is calculated according to Tuan et al.44,45.

$$v_{xj} = \frac{{e^{{r\left( {x + 1} \right)}} }}{{s_{xj} }}\sum\limits_{i = x}^{\infty } {e^{{ - r\left( {i + 1} \right)}} } \sum\limits_{y = j}^{\beta } {s^{\prime}_{iy} f_{iy} }$$
(9)

Raw data on daily predation rates were analyzed according to46 by CONSUME-MSChart47. The age-specific predation rate (kx, the number of prey consumed by the surviving individuals at age x), age-stage specific predation rate (cxj, the mean predation rate of individuals at age x and stage j), age-specific net predation rate (qx), net predation rate (C0), transformation rate (Qp), finite predation rate (ω) and stable predation rate (ψ) were calculated.

The age-stage specific predation rate cxj was calculated as:

$$c_{xj} = \frac{{\sum\limits_{i = 1}^{{n_{xj} }} {d_{xj,i} } }}{{n_{xj} }}$$
(10)

where dxj,i is the calculated predation rate for the ith individual at age x and stage j.

The age specific predation rate (kx) was calculated as:

$$k_{x} = \frac{{\sum\limits_{j = 1}^{\beta } {s_{xj} c_{xj} } }}{{\sum\limits_{j = 1}^{\beta } {s_{xj} } }}$$
(11)

The age specific net predation rate (qx) by considering survival rate was calculated as:

$$q_{x} = l_{x} k_{x} = \sum\limits_{j = 1}^{\beta } {s_{xj} c_{xj} }$$
(12)

The net predation rate (C0) is defined as the number of prey consumed by an average individual predator during its lifetime. It includes all individuals that died in the preadult stages and those that survived to the adult stage. It was calculated as:

$$C_{0} = \sum\limits_{x = 0}^{\infty } {\sum\limits_{j = 1}^{\beta } {s_{xj} c_{xj} } } = \sum\limits_{x = 0}^{\infty } {l_{x} k_{x} }$$
(13)

The transformation rate (Qp) is the number of prey needed for a predator to produce a single offspring. It was calculated as:

$$Q_{p} = \frac{{C_{0} }}{{R_{0} }}$$
(14)

The finite predation rate (ω) was calculated according to Yu, et al.30 as:

$$\omega = \lambda \psi = \lambda \sum\limits_{x = 0}^{\infty } {\sum\limits_{j = 1}^{\beta } {a_{xj} c_{xj} } }$$
(15)

where axj is the proportion of individuals belonging to age x and stage j in a stable age-stage distribution, ψ is the stable predation rate.

Variances and standard errors (SE) of all parameters were determined by using 100,000 bootstraps replications48,47,50. We used paired bootstrap tests to compare differences among the temperatures at 5% significance level based on the confidence interval of difference51. All graphs were plotted with SigmaPlot (ver. 12, Systat Software, Palo Alto, CA).

Population and predation projection

The population and predation projection began with 10 pairs of Harmonia axyridis at all temperatures and it was projected for 90 d to obtain the total population size assuming a scenario without suppression by biotic and/or abiotic factors and predation potential with unlimited prey, i.e., S. litura eggs.

The total population size at time t was calculated as:

$$N\left( t \right) = \sum\limits_{j = 1}^{\beta } {\sum\limits_{x = 0}^{\infty } {n_{xj, \, t} } }$$
(16)

where nxj,t is the number of individuals of age x and stage j at time t52.

The predation potential at time t was calculated following Huang, et al.52:

$$P\left( t \right) = \sum\limits_{j = 1}^{\beta } {\sum\limits_{x = 0}^{\infty } {c_{xj} n_{xj, \, t} } }$$
(17)

To reveal the variability of both projections, we sorted the 100,000 bootstrap results of the finite rate of increase (λ) to find the 2.5th and 97.5th percentiles, i.e., the 2,500th and 97,500th sorted bootstrap samples. We then used the bootstrap life table samples that generated the 2.5th and 97.5th percentiles of the finite rate of increase (λ) to project the population to represent the confidence interval of the projected populations52. Projections were made based on the method described by41,53 and via TIMING- MSChart computer program54.

Results

Development time of each stage of H. axyridis

Increasing temperature from 15 to 30 °C significantly (P < 0.05) reduced development times for eggs, larvae, pupae and adult stages of the H. axyridis (Table 1). When temperature increased from 15 to 30 °C, the developmental time decreased from 7.08 to 2.24 days for first instar, 5.33 to 1.64 days for second instar, 5.19 to 2.11 days for third instar, 8.85 to 3.21 days for fourth instar, 88.32 to 37.61 days for adult males, and 105.94 to 46.65 days for adult females. Increasing temperature from 15 to 30 °C significantly reduced the duration of the total juvenile stages. Male and female longevities were far greater at colder temperature (15 °C) than at warmer temperatures (25–30 °C). The total preadult stage was 15.39 d at 30 °C and 44.26 d at 15 °C.

Table 1 Duration of each stage and preadult survival rate (Mean ± SE) of Harmonia axyridis on Spodoptera litura eggs at four different temperatures.

The sxj of H. axyridis fed with S. litura eggs at four constant temperatures is presented in Supplementary Fig. S1. Results showed an increase in temperature reduced the survival rates of new hatchlings. The survival was of 161 and 171 days at 15 °C, 124 and 131 d at 20 °C, 95 and 112 d at 25 °C, and 61 and 74 d at 30 °C (Supplementary Fig. S1d), with survival peaks of 0.56 and 0.36 at 15 °C, 0.44 and 0.48 at 20 °C, 0.44 and 0.38 at 25 °C, and 0.36 and 0.40, at 30 °C for male and female adults, respectively.

Reproduction and population parameters of H. axyridis

The temperature effects were assessed on the estimated intrinsic rate of increase (r), finite rate of increase (λ), net reproductive rate (R0), and mean generation time (T) for H. axyridis feeding on S. litura eggs by applying the bootstrap technique with 100,000 resamplings (Table 2). Population parameters i.e. r, λ, R0, and T of H. axyridis on S. litura eggs under influence of temperature are presented in Table 2. The r (the intrinsic rate of increase) and λ (finite rate of increase) highest values were 0.1378 d−1 and 1.1477 d−1, respectively, at 30 °C. The shortest mean generation time (35.55 d) was obtained at 30 °C. The lowest r and λ values (0.0662 d−1 and 1.0685 d−1) and the longest generation time (79.84 d) were observed at 15 °C. The net reproductive rate, R0, was significantly higher (237.96) at 20 °C when compared with other temperatures.

Table 2 Reproduction and population parameters (Mean ± SE) of Harmonia axyridis fed on Spodoptera litura eggs at four different temperatures.

The mean fecundity and oviposition days decreased with warming temperatures. Fecundity was highest at 15 °C with 550.06 eggs per female, which decreased as the temperature increased i.e., 495.75 eggs at 20 °C, 484.74 eggs at 25 °C, and 334.80 eggs at 30 °C. The oviposition days were significantly shortest at 30 °C with 16.65 d and longest at 15 °C with 26.94 d. Temperature impacts were obvious on the adult pre-oviposition period (APOP) and total pre-oviposition period (TPOP) (Table 2). APOP of females lasted shortest at 30 °C (i.e., 8.85 days) and longest at 20 °C (i.e., 12.88 d). TPOP was shortest at 30 °C (i.e., 24.25 d) and colder temperatures prolonged it from 56.17 d at 15 °C to 45.46 d at 20 °C, and 30.63 d at 25 °C.

Figure 1 shows lx (age-specific survival rate), mx (total fecundity of population), and lxmx (net maternity) of H. axyridis population fed with S. litura eggs. The effects of temperature were observed on all the mentioned parameters. The lx is described as the simple form of the sxj curves. The age-specific survival rate of H. axyridis population fed with S. litura eggs at 15 °C decreased to 0.92 on the 47th d when the first fecundity occurred. At 20 °C, the first fecundity was obtained on the 42nd d with 0.92 survival rates. The fecundity of the population was seen early at 25 °C on the 27th d with 0.82 lx and 22nd d with survival rate decreasing to 0.76 at 30 °C. The lxmx curves depended on lx and mx values. The highest lxmx peaks in H. axyridis populations were 5.24 (on 70th d), 6.18 (on 78th d), 5.44 (on 58th d) and 7.32 (on 39th) (Fig. 1).

Figure 1
figure 1

Age-specific survival rate (lx), age-specific fecundity (mx), and age-specific maternity (lxmx) of the Harmonia axyridis feeding on Spodoptera litura eggs at different temperatures.

The exj (age-stage life expectancy) of H. axyridis fed with S. litura eggs at four constant temperatures are plotted in Supplementary Fig. S2. The exj estimates the length of duration of the newly hatched individual. Results showed temperature impacts exj curve of H. axyridis, but with varied degrees, i.e. highest life expectancy values for females and males were 111.61 d and 93.32 d at 15 °C, 84.83 d and 75.96 d at 20 °C, 72.21 d and 63.46 d at 25 °C, 72.21 d and 63.46 d at 30 °C (Supplementary Fig. S2). Supplementary Fig. S3 presents vxj (age-stage-specific reproduction) of H. axyridis fed with S. litura eggs at four constant temperatures. The vxj tells the contribution every individual will likely have in the next progeny. The curves show that temperature change also changes the age stage-specific reproduction. The highest vxj value of H. axyridis female fed with S. litura eggs was 102.16 eggs/d (on 61st d) at 15 °C which changed to 95.53 eggs/day (on 61st d) at 20 °C, 92.54 eggs/d (on 42nd d) at 25 °C and 85.11eggs/d (on 34th d) at 30 °C (Supplementary Fig. S3).

Predation rate of H. axyridis

Daily predation of 1st-4th instar increased at all tested temperatures, exponentially. Daily mean prey consumption during the pre-adult stages was significantly different between temperatures as 22.57, 39.87, 44.86 and 52.46 eggs/larva at 15, 20, 25 and 30 °C, respectively (Table 3). Female predators consumed significantly more than did males within the same temperature and total daily predation of adults varied significantly with 79.58, 130.50, 172.02 and 198.67 eggs/adults at 15, 20, 25 and 30 °C, respectively (Table 3). The age-stage specific predation rate (cxj) of H. axyridis which means the mean predation rate of individuals at age x and stage j was presented in Supplementary Fig. S4. While the long survival durations increase total predations in lower temperatures i.e. 15 and 20 °C, daily predation through the stages was generally higher at 25 and 30 °C. The net predation rate (C0) was greatest i.e., 10,466.28 and 10,139.38 at 20 and 25 °C than those at 15 and 30 °C with 7935.54 and 7126.36 eggs/individual, respectively (Table 4). There was no significant difference in transformation rate (Qp) at all tested temperatures. The transformation rate (Qp) was highest at 25 °C, which indicates that H. axyridis requires 55.05 S. litura eggs for the production of each egg at 20 °C. Stable and finite predation rates differed among all tested temperatures. The highest values for both parameters were obtained at 30 °C, whereas lowest values were at 15 °C.

Table 3 Stage daily predation (Dj) (preys/individual) (Mean ± SE) of Harmonia axyridis fed on Spodoptera litura eggs at four different temperatures.
Table 4 Predation rates (Mean ± SE) of Harmonia axyridis fed on Spodoptera litura eggs at four different temperatures.

Age-specific predation rate (kx), age-specific net predation rate (qx) cumulative net predation rate (C0) for H. axyridis feeding on S. litura eggs at four different temperatures were presented in Fig. 2. The kx gives the average number of S. litura eggs eaten by H. axyridis of age x, while qx gives the weighted number of aphids eaten by H. axyridis of age x. As H. axyridis feed with S. litura eggs only in the larval and adult period, the increase in kx and qx was linked to survival rates until the adult period with a stepwise reduction beginning from the adult period. There were two gaps in the curves for kx and qx at all tested temperatures due to non-predator i.e., egg and pupa stages of the H. axyridis. Cumulative predation rates were 7935 (on 171st d), 10,466 (on 131st d), 10,139 (on 113rd d) and 7126 (on 75th d) at 15, 20, 25 and 30 °C, respectively (Fig. 2).

Figure 2
figure 2

Age-specific survival rate (lx), age-specific predation rate (kx), age-specific net predation rate (qx) cumulative net predation rate (C0) of Harmonia axyridis feeding on Spodoptera litura eggs at different temperatures.

Population and predation projection

Both the population and predation projection for 90 days began with 10 H. axyridis pairs at all temperatures (Fig. 3). The variability of both population growth and predation potential was assessed by utilizing the life tables that generated the 0.025th and 0.975th percentiles of λ of 100,000 bootstrap results. The total population size of H. axyridis was highest at 30 °C, projected to exceed 8.9 × 106 individuals, whereas the population at 15 °C was estimated at 4.3 × 104 individuals. The total population at 25 and 30 °C increased significantly than those at 15 and 20 °C. At lower temperatures (i.e., 15 and 20 °C), the curves of different stages tended to be straight lines (Fig. 3). Total consumption of S. litura eggs was highest and projected to reach 5.24 × 108, while egg consumption was approximately 4.9 × 105 at 15 °C after 90 days (Fig. 3).

Figure 3
figure 3

Population and predation projection over a 90-day period of Harmonia axyridis feeding on Spodoptera litura eggs at different temperatures (the original cohort and the cohorts constructed based on the 2.5 and 97.5% percentiles of λ, finite rate of increase are all in log base 10).

Discussion

Coccinellid species are generalist predators and the type of prey they consume shows great impacts on their fitness and predation. The amount and quality of food is very important because it influences directly biological aspects of predators55. When ingested food is scarce or inferior quality, the development time usually increases and the reproductive rate, i.e., oviposition, fecundity and fertility decreases56. Given the relevance of coccinellids for biological control, much attention has been paid to documenting feeding habits within this group. Food that provides better development and reproduction is considered essential for coccinellids. Many aphids, psyllids and mealybugs are important preys. Many other energy sources as preys that prolong survival are characterized as alternative57. Coccinellidae are shown to reproduce even more on non-aphid preys in some studies, which demonstrates the importance of alternative resource58,59. Eggs of lepidopterans were the best diet for H. axridis because shorter duration and higher survival of the larval stage of this predator was evident on the eggs34,60. Abdel‐Salam and Abdel‐Baky34 estimated the R0 of H. axyridis when feeding on frozen and fresh eggs of the Sitotrega cerealella, and showed moth eggs to be a better diet based on findings obtained by applying the female age-specific life table and adult age. This shows that scales and chorion of eggs were not barriers for larval as well as adults feeding, as reported for H. axyridis with S. litura eggs35. The current study is the first reporting on the survivorship, development, reproduction, and longevity of H. axyridis on S. litura eggs at different temperatures. The current study reveals evidence of temperature influence on the development and predation of H. axyridis. We also used computer projections based on bootstrap percentile confidence intervals to better understand the variability of population growth and predation parameters.

In our findings, the development and reproduction of H. axyridis on S. litura eggs were temperature-dependent. The total larval duration was the longest at 15 °C, intermediate at 20 and 25 °C, and the shortest at 30 °C. The development periods of the egg, preadult stages, as well as total development period of the adult stages (Females: Males), were significantly different in length with the four different temperatures tested. The shortest developmental periods for the preadult and adult stages of H. axyridis occurred at 25 °C and 30 °C. The shorter development duration of preadult and adult H. axyridis at higher temperatures (i.e., 30 °C) is most likely due to increased attack and metabolic rate in response to increased food consumption and energy requirements, allowing the predator to rapidly advance to the next stage. The mean longevities of male and female adults also prolonged at 15 °C and reduced with warming temperatures, in accordance with De Oliveira Ramos, et al.60, who under similar thermal conditions reared H. axyridis on different prey, Anagasta kuehniella (Zeller) (Lepidoptera: Pyralidae). de Castro-Guedes, et al.33 reported that female H. axyridis could survive on average 91.13 ± 25.53 d on the eggs of A. kuehniella at 25 °C. Abdel‐Salam and Abdel‐Baky34 reported female H. axyridis to survive on average 62.2 ± 2.70 d on S. cerealella eggs at 24 °C. The fecundity of female H. axyridis fed on different lepidopterans eggs also changed with temperature. For example, the fecundity of H. axyridis increased from 806.00 ± 78.4 to 1.685.2 ± 209.6 per female as the temperature increased from 18 to 30 °C. Similarly, the fecundity was recorded 715.3 ± 33.62 per female, when H. axyridis was allowed to feed on S. cerealella eggs at 24 °C. Accelerated metabolic rate, which leads to increased food consumption to satisfy energy requirements may has allowed the predator complete its development rapidly and advance to the next stage, relatively in a short period of time under influence of warmer temperatures.

Insect population dynamics rely heavily on fecundity61, was the highest at 15 °C (i.e., 550.06 eggs) in our findings, suggesting S. litura to be better for the development of H. axyridis at lower temperature than higher temperatures. The findings support H. axyridis can lay more eggs at optimum temperature, which is in consistence with the previous research27. Castro, et al.27 reported that H. axyridis fed with C. atlantica showed the greatest fecundity at 15 °C (i.e., 805.7 ± 127.3) than when fed on the same prey at 20 and 25 °C (i.e., 608.5 ± 113.7 and 614 ± 129.2, respectively). In our experiment, the fecundity comparatively reduced at higher temperatures but not the predation that increased with warmer temperatures. Higher temperatures can reduce reproductive fitness of H. axyridis, and other researchers have reported similar findings62, which means H. axyridis survives thermally complex and diverse environments with such a trade off. This research was conducted using total fecundity without the exclusion of unhatched eggs, however it was showed by Mou, et al.63 that the life table and predation parameters differ depending on whether the egg viability factor is included or excluded. As coccinellid's eggs are not always fertile, and as that fertility rate varies with temperature and predator diet, with resulting consequences on biological procedures27, other investigations should look at fertility as well as fecundity to gain a better understanding of this predator's biological performance. The duration of adult pre-oviposition days including APOP and TPOP was shortened with higher temperatures. These results are in agreement with those of other researchers such as Castro, et al.27, who reported the highest and the lowest APOP and TPOP for H. axyridis when it was fed with C. atlantica at 15 °C (6.10 and 82.40 d) and 25 °C (5.80 and 76.87 d), respectively.

We found that predation of H. axyridis was also temperature-dependent and increased with warmer temperature. Warming is expected to expedite predation because it boosts the predator's capacity to grab and handle prey in response to increased metabolic rate that can result from increased development and fostered growth12. In comparison to first three instars, the fourth larval instar and female H. axyridis consumed more S. litura eggs and showed the maximum predation capacity at these stages. Other researchers have reported that the predation rate of preadult stages of H. axyridis generally increased with larval sizes12,35. Females consumed more eggs than did males because females need to oviposit. The sexes differed in our findings in a similar way. The female age-specific fecundity correlates with the trend of predation rate. Similar results have also been reported by another researcher who showed fourth instar of Hippodamia variegata (Coleoptera: Coccinellidae) to consume significantly more Aphis fabae Scopoli (Hemiptera: Aphididae) than other stages of the predator at 23 °C64.

This study concludes that H. axyridis development depends on the temperature and predation varies with predator stage, temperature, and resource being consumed. On the eggs of S. litura, longevity , reproductive rate and mean generation time of H. axyridis were much improved at 15 and 25 °C than at 25 and 30 °C. Lower temperatures can thus be favorable for mass rearing the predator. Warmer temperatures (30 °C) can lead to an accelerated development and growth of this predator, which means an increased predation rate. These findings may help with devising effective biocontrol programs against S. litura using H. axyridis as the predator, probably by identifying potential predator stage and optimal release time. This study further provides useful information about the future consequences of warming on this predator development and predation parameters, which can help attune bio-based IPM programs according to changing climates11,12. Hance, et al.65 reported that divergences between thermal optimums of predators and preys may promote disturbance of geographical or temporal synchronization, increasing the risk of pest outbreaks. Temperature limits pest distribution; however, warming has resulted in more suitable disciplines for pests, e.g., insect range expansions. Our experiments used a thermal maximum of 30 °C, and temperatures in Asia can reach as high as > 40 °C for short periods of time at midday during hot summer66. Presumably, S. litura in this region is locally adapted to short bouts of high heat and can withstand it, although additional research is needed to confirm that host eggs exposed to high heat are suitable for predator development. Although the present study was performed under laboratory conditions, it may still provide information that can complement field studies. Future field experiments under controlled conditions are recommended to provide a much realistic view of the effectiveness of H. axyridis in the field against this serious pest.