Bumblebees. Experiments were conducted between April and May 2019. Twenty-five commercial colonies of B. terrestris were purchased from Koppert Agricultural Co., Ltd. (Beijing, China). Each colony contained about 200 workers, a brood at all developmental stages, and a laying queen. The bumblebees were reared on a diet that included pollen and nectar and were provided by the company in an incubator with continuous darkness, at a temperature of 25 ± 1 ℃ and a relative humidity of 60 ± 10%.
Insecticides. Chlorpyrifos (CAS No. 2921-88-2, 96% Technical material) was supplied by the Hunan Research Institute of Chemical Industry (Hunan, China). Imidacloprid (CAS No. 138261-41-3, 96% TC) was supplied by Shandong Zhongnong United Biological Technology Co., Ltd. (Shandong, China). Thiamethoxam (CAS No. 153719-23-4, 97% TC) was obtained from the Hailier Pesticides and Chemicals Group (Shandong, China). Each insecticide was dissolved in dimethyl sulfoxide and diluted in a 50% (w/w) sugar solution as described by Yue et al.34. Each stock solution was diluted to six test concentrations by using a calibrated micropipette and volumetric flasks.
Toxicity assessment. The acute oral toxicity of the insecticides to the worker was tested according to the method recommended by the OECD35(Organization for Economic Co-operation and Development). Briefly, one leg of the same size workers was clamped gently with a forceps, and the bees were quickly transferred to a thermostat-controlled wooden box (dimensions 12 cm ´ 8 cm ´ 8 cm; Fig. 1). Fifteen workers were placed in each wooden box in the dark at room temperature (25 ± 1 ℃) and a relative humidity of 60 ± 10% with a sufficient amount of noncontaminated 50% sugar solution (w/w). The bees were left alone for at least 8 hours for adaptation. The experiment was conducted when the mortality rate of bumblebees in the wooden box did not exceed 10%. Then 300 μL of the 50% sugar solution was either contaminated with a insecticide or fed uncontaminated to the worker bumblebees via a 5 mL syringe with the tip removed (Fig. 2) for 6 hours, followed by 2 hours of starvation. The sugar solution was immediately replaced with a sufficient amount of uncontaminated sugar solution once the 300 μL sugar solution had been consumed over the 6 hours. The mass of each test solution was weighed and recorded before and after each feeding. Six groups of workers (a total of 90 bees) from the same colony were treated with a sugar solution containing six different concentrations of insecticides. Triplicate experiments were performed with 270 worker bees for each treatment. In total, seven treatments were performed: (1) chlorpyrifos, (2) imidacloprid, (3) thiamethoxam, (4) chlorpyrifos + imidacloprid, (5) chlorpyrifos + thiamethoxam, (6) imidacloprid + thiamethoxam, (7) chlorpyrifos + imidacloprid + thiamethoxam. The total number of bees used for the experiment was 1,890.
All binary and ternary insecticide toxicities were administered as described by Liu et al.36. Stock solutions of chlorpyrifos, imidacloprid, and thiamethoxam were prepared as described above and used in three binary combinations (chlorpyrifos + imidacloprid; imidacloprid + thiamethoxam; chlorpyrifos + thiamethoxam) and a ternary combination (chlorpyrifos + imidacloprid + thiamethoxam). The constant combination ratios were chlorpyrifos:imidacloprid, 0568:0.310; imidacloprid:thiamethoxam, 0.31:0.438; chlorpyrifos:thiamethoxam, 0.568:0.438; and chlorpyrifos:imidacloprid:thiamethoxam, 0.568:0.310:0.438 based on the individual median lethal dose (LD50) toxicity such that the effects of the individual insecticides within the combination would be approximately equal. In addition, the mixed insecticides were diluted to six concentrations. All six concentrations were tested simultaneously to minimize experimental variations.
Data analysis. A preliminary experiment suggested that evaporation of the sugar solution in the syringe did not significantly affect the mass change (a loss of about 0.001 g). Therefore, the consumption of the sugar solution could be inferred from the differences before and after insecticide exposure. The mixtures were then converted from concentrations into doses in micrograms of the active ingredient per worker. The LD50 values were calculated by probit analysis using POLO-PC software37.
The individual and combined toxic effects of insecticides on bumblebees were assessed using the median-effect equation described by Liu et al.36 and Chou and Talalay37:
fa/fu = (D/Dm)m (1)
where D is the dose of an insecticide, Dm is the dose for a 50% effect, fa is the mortality influenced by D (percentage of mortality), fu is the survival rate uninfluenced by D (percentage of survival, fu = 1 - fa), and m is the coefficient determining the shape of the dose–effect relationship.
By rearranging equation (1), we can obtain the following equations:
fa = 1/[1 + (Dm/D)m] (2)
D = Dm[fa/(1 - fa)]1/m (3)
Therefore, if we know the values for m and Dm, we can easily assess the effect (fa) for any given dose (D) in equation (2). In the same way, the dose (D) can easily be calculated by the effect (fa) given in equation (3). In addition, if we take the logarithm of both sides of equation (1) and assume that x = log(D) and y = log(fa/fu), we can obtain the following middle-effect diagram:
log(fa/fu) = mlog(D) - mlog(Dm) (4)
In the median-effect plot in equation (4), we can easily determine the Dm, where m for the Dm means the antilog of the x-intercept and m is the slope. Here, m > 1, m = 1, and m < 1 signify sigmoidal, hyperbolic, and flat sigmoidal dose–effect curves, respectively. In addition, the linear correlation coefficient (r) of the median-effect plot can reveal how the data conform to the median-effect plot, where r = 1 shows excellent conformity.
Therefore, we can easily calculate the combination index (CI) values by using the CI equation for a combination of n insecticides, which is given as
(5)
where (CI)x is the combination index for n insecticides at x% effect (fa); (Dx)1-n is the sum of the dose of n insecticides causing x% effect (fa) in combination; [D]j/ is the proportionality of the dose of n individual insecticides causing x% effect (fa) in combination; (Dm)j{(fax)j/[1 - (fax)j]1/mj} is the dose of individual insecticides causing x% effect (fa); and fax is the fractional effect (fa) at x% effect (fa), where CI > 1, CI < 1, and CI = 1 indicate an antagonistic, synergistic, and an additive effect, respectively.
The computer program CompuSyn39 was used to calculate the parameters including the dose–response curve parameters, CI values, fa–CI plot representing CI versus fa, the fraction influenced by a specific dose, and the polygonogram representation describing the antagonistic, additive, or synergistic effect of the insecticide combination.