Study site
The field study was conducted in agricultural districts of Ulanga and Kilombero in south-eastern Tanzania, where the main malaria vectors include Anopheles funestus and Anopheles arabiensis [52,54]. Four villages were covered, namely Minepa (V1) 8.23°S, 36.75°E; altitude 268m) and Lupiro (V2) (8.41°S, 36.81°E; altitude = 389m) in Ulanga district and Kisawasawa (V3) (7.91°S, 36.82°E; altitude = 728m) and Njage (V4) (8.3°S, 36.14°E; altitude = 519m) in Kilombero district. The area experiences 15 - 35°C daily temperatures, and 1300 - 3600mm annual precipitation [55]. Economic activities include crop-farming, fishing, and livestock keeping [16,56]. The main malaria vector control method in the area is ITNs [57].
Qualitative assessment of the use of agricultural pesticides
A qualitative assessment was conducted in the four study villages, by carrying out focus group discussions (FGDs) to identify the common agricultural pesticides used in each village, followed by direct observations of pesticide use in the farms.
A discussion guide was prepared for the FGDs to capture the following key areas of discussion: a) farming practices, b) types of chemicals used in the farm, c) farming seasons and crop types, d) sources of agricultural pesticides and e) management and disposal of the pesticides. An FGD was conducted in each village, bringing together eight local community representatives recruited upon signing the informed consent. These representatives were community leaders, practicing subsistence farming, and included both males and female adults. The discussions were facilitated by three researchers from Ifakara Health Institute, who were all knowledgeable on agricultural practices and malaria control interventions in the region and had experience conducting qualitative research in rural communities. Before starting the discussions, facilitators presented a brief overview of the study topic and the reason for conducting the FGDs. Each session lasted about 60-75 minutes.
In addition, direct observations were made in the farms in the study villages to record signs of pesticide use, disposal, and methods of application. Where possible, pictorial evidence was gathered to complement the direct observations.
Mosquito collections and rearing of larvae
Field mosquitoes were collected in their larval stages during the dry season months of July and September 2021. The study villages were inspected for presence of An. arabiensis habitats in or near the rice farms before the actual mosquito collections were done. The standard 350mL dippers were used to collect the larvae; and the collections were transported to the VectorSphere Mosquito Laboratory facility at Ifakara Health Institute for further experiments. In the insectary, the larvae were kept in 5L rearing dishes, under controlled conditions of 25-27°C and 80% relative humidity. The larvae were fed on Tetramin® fish food and the emergent adults were supplied with a 10% sugar solution. As described below, An. arabiensis females from a laboratory colony maintained under the same conditions since 2009, were also used in some comparative tests.
Assessment of insecticide susceptibility of the field-collected mosquitoes
Emergent adult female mosquitoes aged 3–5-day old and not previously blood-fed were used according to the WHO insecticide susceptibility test procedures [58]. A total of 120 mosquitoes per test were used over four replicates, each requiring 20–25 mosquitoes. The bioassays were done at 25°C ± 2 °C and 80% ± 10% relative humidity.
The candidate insecticides and doses included: 0.75% permethrin (pyrethroids class I), 0.05% lambda-cyhalothrin (pyrethroid class II), 0.25% pirimiphos-methyl (organophosphate), 4% DDT (organochloride) and 0.1% bendiocarb (carbamate). The mosquitoes were exposed for one hour, and knockdown recorded at intervals of 10, 15, 20, 30, 40, 50, 60 minutes. The mosquitoes were then transferred to clean holding tubes, supplied with 10% sugar solution, and observed to assess mortality after 24 hours. For comparison, similar tests were conducted using the laboratory-reared mosquitoes.
The synergist, 4% Piperonyl Butoxide (PBO), was used to assess possible involvement of mixed-function oxidases in metabolic resistance [58]. In these tests, the mosquitoes were pre-exposed to PBO or control and then to the candidate insecticides. Four groups of 20-25 mosquitoes each were used, and the groups were exposed for one hour to either PBO alone, PBO followed by respective candidate insecticide for one hour, the candidate insecticide alone, or the control papers. The mosquitoes were monitored for 24-hour mortality and the tests replicated three times.
Assessment of fitness parameters of field-collected mosquitoes
Field collected larvae from the four study villages (V1-V4) that were brought to the insectary and kept till they pupated, after which they were collected in cups and placed in standard 15 x 15 x 15cm cages. Upon adult emergence, the mosquitoes were fed on 10% glucose solution via soaked cotton pads. Mosquitoes were provided their first blood meal on the third day after emergence, to give them time to mate. Fully fed mosquitoes were moved into individual cups with a damp filter paper at the base of the cup to stimulate oviposition conditions. Eggs laid by each mosquito were counted under a stereomicroscope as a measure of fecundity.
Wing sizes were also measured and used as a proxy for body size. Adult female mosquitoes emerging from the field larvae collections, were anaesthetized at -10°C for seven minutes and a single wing was removed from either the left or right side of the mosquito. Distilled water was used to fix the wings onto the slides and the length from the apical notch to the auxiliary margins was measured using a micrometer ruler under a stereo microscope.
Exposure of mosquitoes to sub-lethal doses of agricultural pesticides
These tests were done to simulate field exposures to common agrochemicals and assess effects on emergent adults. Groups of 3rd instar An. arabiensis larvae from the laboratory colonies were exposed to a range of pesticide concentrations (Table 1), to observe the effect of the active treatments and determine the doses below which there were still substantial emergent adults for further experimentation. For this experiment, three insecticides from three pesticide classes (pyrethroids, organophosphates and carbamates) were selected, so as to be representative of the commonest pesticide classes used by communities in the study villages (as determined by the FGDs and the field observations).
In each test, 25 larvae were introduced into 1.2L basins containing a 1L of water, into which the candidate pesticides were introduced. The pesticide concentrations ranged from 1´10-4 to 1´10-8 g/L (for lambda cyhalothrin and pyrimiphos methyl) and 5´10-4 to 3.5´10-3 g/L (bendiocarb). Three replicates were completed for each concentration (treatment) and control, all under the same conditions. All larvae exposed and unexposed were fed on Tetramin ® fish food and monitored every 12 hours for mortality until 120hours.
A lethal concentration (LC) of 15% mortality was selected for subsequent experiments, since at this level there was sufficient adult emergence despite significant pesticide exposure. Fresh batches of fourth-instar larvae were therefore exposed to the LC15 pesticides concentrations in the 1.2L trays for 48hours, during which the larvae were fed Tetramin® fish food. Pupae from these trays were transferred to cups containing only water and placed in small cages (15 × 15 × 15cm) for emergence.
Assessment of insecticide susceptibility and fitness parameters of the laboratory mosquitoes
Once the sub-lethal dose was determined, the 3rd to 4th instar larvae were introduced and maintained until pupation. The pupae collected from each treated dish (exposed mosquitoes) and from concurrent controls (unexposed mosquitoes) were counted and placed in cages for emergence. Three-to-five-days old female adult mosquitoes were then used to carry out insecticide susceptibility bioassays as described above for field-collected mosquitoes. Fecundity and wing length were also assessed in the same way as described above for field-collected mosquitoes.
Data analysis
The quantitative data was analysed using the open-source software, R programming version 4.0.5 [59]. In tests using mosquitoes pre-exposed to sub-lethal aquatic doses of pesticides, Probit analysis was done using the ‘ecotox ’ library [60] to determine appropriate sub-lethal concentrations to be used in subsequent experiments. To assess the insecticide resistance profiles, data analysis was done according to the WHO susceptibility test guidelines. The resistance or susceptibility status was defined based on the WHO criteria, where 98–100% mortality indicates susceptibility; 90–97% mortality requires further confirmation of possible resistance, and less than 90% mortality indicates resistance [58]. The resistance/susceptibility graphs were plotted using ggplot2 package [61]. Knockdown times were calculated using the PoLo Plus software [62] using log-probit analysis. For fecundity and wing length, the dabestr package in R was used to generate a two-group estimation (Gardner-Altman) plots [63]. The estimation plots were used to display distributions of residual mean differences in the number of eggs laid and wing lengths in the different experiments. The results were presented in summary graphs or tables.
For the qualitative assessment, audio recordings of the focus group discussions and key informant interviews were transcribed and translated from Swahili to English. Any notes written during the discussion were added to the written transcripts. The written transcripts were reviewed and analysed on Microsoft word. Analyses were conducted separately for different villages. The FGD guide and objective of the study, were used to develop deductive codes. Inductive codes were derived from detailed studying of the written transcripts. Once the coding was completed, codes were grouped, and emerging patterns were used to identify themes.