Dramatic resurgence of malaria after 7 years of intensive vector control interventions in Eastern Uganda ======================================================================================================== * Moses R. Kamya * Joaniter I. Nankabirwa * Emmanuel Arinaitwe * John Rek * Maato Zedi * Catherine Maiteki-Sebuguzi * Jimmy Opigo * Sarah G. Staedke * Ambrose Oruni * Martin J. Donnelly * Bryan Greenhouse * Jessica Briggs * Paul J. Krezanoski * Teun Bousema * Philip J. Rosenthal * Peter Olwoch * Prasanna Jagannathan * Isabel Rodriguez-Barraquer * Grant Dorsey ## Abstract **Background** Tororo District, Uganda experienced a dramatic decrease in malaria burden from 2015-19 following 5 years of indoor residual spraying (IRS) with carbamate (Bendiocarb) and then organophosphate (Actellic) insecticides. However, a marked resurgence occurred in 2020, which coincided with a change to a clothianidin-based IRS formulations (Fludora Fusion/SumiShield). To quantify the magnitude of the resurgence, investigate causes, and evaluate the impact of a shift back to IRS with Actellic in 2023, we assessed changes in malaria metrics in regions within and near Tororo District. **Methods** Malaria surveillance data from Nagongera Health Center, Tororo District was included from 2011-2023. In addition, a cohort of 667 residents from 84 houses was followed from August 2020 through September 2023 from an area bordering Tororo and neighboring Busia District, where IRS has never been implemented. Cohort participants underwent passive surveillance for clinical malaria and active surveillance for parasitemia every 28 days. Mosquitoes were collected in cohort households every 2 weeks using CDC light traps. Female *Anopheles* were speciated and tested for sporozoites and phenotypic insecticide resistance. Temporal comparisons of malaria metrics were stratified by geographic regions. **Findings** At Nagongera Health Center average monthly malaria cases varied from 419 prior to implementation of IRS; to 56 after 5 years of IRS with Bendiocarb and Actellic; to 1591 after the change in IRS to Fludora Fusion/SumiShield; to 155 after a change back to Actellic. Among cohort participants living away from the border in Tororo, malaria incidence increased over 8-fold (0.36 vs. 2.97 episodes per person year, p<0.0001) and parasite prevalence increased over 4-fold (17% vs. 70%, p<0.0001) from 2021 to 2022 when Fludora Fusion/SumiShield was used. Incidence decreased almost 5-fold (2.97 vs. 0.70, p<0.0001) and prevalence decreased by 39% (70% vs. 43%, p<0.0001) after shifting back to Actellic. There was a similar pattern among those living near the border in Tororo, with increased incidence between 2021 and 2022 (0.93 vs. 2.40, p<0.0001) followed by a decrease after the change to Actellic (2.40 vs. 1.33, p<0.001). Among residents of Busia, malaria incidence did not change significantly over the 3 years of observation. Malaria resurgence in Tororo was temporally correlated with the replacement of *An. gambiae s.s.* by *An. funestus* as the primary vector, with a marked decrease in the density of *An. funestus* following the shift back to IRS with Actellic. In Busia, *An. gambiae s.s.* remained the primary vector throughout the observation period. Sporozoite rates were approximately 50% higher among *An. funestus* compared to the other common malaria vectors. Insecticide resistance phenotyping of *An. funestus* revealed high tolerance to clothianidin, but full susceptibility to Actellic. **Conclusions** A dramatic resurgence of malaria in Tororo was temporally associated with a change to clothianidin-based IRS formulations and emergence of *An. funestus* as the predominant vector. Malaria decreased after a shift back to IRS with Actellic. This study highlights the ability of malaria vectors to rapidly circumvent control efforts and the importance of high-quality surveillance systems to assess the impact of malaria control interventions and generate timely, actionable data. Key words * malaria * effective malaria control * resurgence * long-lasting insecticidal nets (LLINs) * indoor residual spraying of insecticide (IRS) * clothianidin * pirimiphos-methyl ## Introduction Between 2000 and 2015 it was estimated the incidence of malaria in Sub-Saharan Africa decreased by 40%, with the scale up of vector control interventions responsible for the majority of cases averted.(1) However, since 2015 progress has stalled and even reversed course in some of the highest burden countries of Africa.(2) Indeed, the global burden of malaria has become increasingly concentrated in Africa, which accounted for over 94% of cases and 96% of deaths in 2022.(2) Turning the tide on malaria will require a better understanding of the root causes of stalled progress, better use of local data to inform policy decision making, and a more flexible and targeted approach to control interventions.(3) Uganda is emblematic of other high burden African countries, where progress in reducing the burden of malaria has been slow and difficult to sustain despite the scale up of proven control interventions. Uganda was the first country to implement universal distribution of free long-lasting insecticidal nets (LLINs) starting in 2014, with repeated campaigns every 3-4 years. Uganda also has one of the largest indoor residual spraying of insecticide (IRS) programs, focusing on selected high transmission districts using different formulations of insecticides rotated every few years. The success of this intensive approach to vector control has been well documented in Tororo District, a historically high transmission area of southeastern Uganda, where our group has been conducting comprehensive cohort-based malaria surveillance studies since 2011.(4) IRS was first implemented in Tororo District in December 2014, initially using a carbamate (Bendiocarb), then switching to an organophosphate (Actellic) in 2016. Comparing key malaria indicators measured from 2011-14, prior to IRS, with indicators in 2017-19, after two rounds of universal LLIN distribution and 5 years of sustained IRS, we documented a 500-fold decrease in malaria transmission intensity, a 60-fold decrease in the incidence of symptomatic malaria, and a 5-fold decrease in parasite prevalence among children 0.5-10 years of age.(5) In addition, over this time frame there was a marked shift in the predominant vector species from *An. gambiae s.s.* to *An. arabiensis*. However, following a change to clothianidin-based formulations of IRS in 2019-20, a marked resurgence of malaria cases was documented using health facility-based data from 5 districts of Uganda (including Tororo), reaching pre-IRS levels within 1-2 years.(6) To better quantify the magnitude of the malaria resurgence in Tororo and to investigate potential causes, we compared temporal changes in malaria incidence, prevalence, and entomological measures between September 2020 and September 2023 in a cohort of 667 residents living in two areas within Tororo District and in neighboring Busia District, where IRS has never been implemented. In addition, we assessed changes after the formulation of IRS in Tororo District was shifted back to an organophosphate in March 2023. ## Methods ### Study setting and population level vector control interventions Our team has conducted a series of comprehensive cohort-based malaria surveillance studies (referred to as PRISM) in southeastern Uganda starting in October 2011 (Figure 1). Briefly, the “PRISM 1” study was conducted from October 2011 – September 2017 and the “PRISM 2” study conducted from October 2017 – October 2019 in Nagongera subcounty, Tororo District.(5) The “PRISM Border Cohort” study (the focus of this report) was conducted from August 2020 – September 2023 in both Tororo District and neighboring Busia District. Tororo is historically a high malaria transmission district with an estimated entomological inoculation rate (EIR) of 310 infective bites per person per year in 2011-12.(7) Prior to 2013, vector control in Tororo District was limited to the distribution of LLINs through antenatal care services. In November 2013, universal distribution of free LLINs was conducted as part of a national campaign, and similar campaigns were repeated in May 2017 and June 2020. All LLIN distribution campaigns in Tororo District utilized standard pyrethroid LLINs. IRS with the carbamate (Bendiocarb) was introduced for the first time in Tororo in December 2014–January 2015, with additional rounds administered in June–July 2015 and November–December 2015. In June–July 2016, the IRS compound was changed to an organophosphate (Actellic), with repeated rounds in June–July 2017, June–July 2018, and March–April 2019. The IRS program changed to clothianidin-based formulations in 2020, with Fludora Fusion (clothianidin + deltamethrin) administered in March 2020 and March 2021, and SumiShield (clothianidin alone) administered in March 2022. In March 2023, a decision was made to shift back to IRS with Actellic in response to the resurgence of malaria. Busia is a high malaria transmission district bordering Tororo District to the south. As in Tororo, prior to 2013 vector control in Busia District was limited to targeted distribution of LLINs through antenatal care services. Universal distribution of free LLINs were conducted in Busia District in May 2013, May 2017, and December 2020. The first two LLINs distribution campaigns in Busia District utilized standard pyrethroid LLINs and the third campaign utilized LLINs containing deltamethrin plus piperonyl butoxide (PermaNet 3.0). IRS has never been implemented in Busia District. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/16/2024.03.15.24304352/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/F1) Figure 1: Map of eastern Uganda with an extract showing study sites. ### Routine health facility-based malaria surveillance data To provide context, we illustrate continuous trends in disease burden from October 2011 – September 2023 at Nagongera Health Center IV, Tororo District (Figure 2). This health center is part of a national health facility-based malaria surveillance network that provides high-quality data to monitor geographic and temporal trends in malaria burden and assess the impact of population level control interventions.(6) Briefly, routine individual level data are collected from all patients who present to the outpatient departments, including the number of laboratory confirmed cases of malaria. Since November 2018, after mapping of the region, we have been able to estimate the incidence of malaria from a target area around Nagongera Health Center IV, defined as the number of cases of laboratory-confirmed malaria diagnosed at the health center among patients residing in the target area, per unit time, divided by the total population of the target area. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/16/2024.03.15.24304352/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/F2) Figure 2: Monthly trends in laboratory confirmed cases of malaria and malaria incidence from Nagongera Health Center. PY = person years. ### Screening and enrollment of households into the cohort study Screening and enrollment of households and study participants was previously described.(8) Briefly, in May 2020, all households in a contiguous study area including two sub-counties in Tororo District and one sub-county in Busia District were enumerated and mapped using handheld global positioning systems (Garmin e-Trex 10 GPS unit, Garmin International Inc., Olathe, KS) to provide a sampling frame for recruitment of households into the cohort study. Households were enrolled from three distinct geographic regions (Figure 1). In August 2020, 20 households from Busia District located 1.9-3.2 km from the border with Tororo District and 30 households from Tororo district located 0.7-3.5 km from the border (referred to as Tororo, near border) were enrolled. Preliminary data revealed that malaria burden was surprisingly similar between houses in Busia District and those in Tororo near the border.(8) To capture more of a gradient in transmission intensity, an additional 30 households from Tororo District located 5.5-10.8 km from the border (referred to as Tororo, away from border) were enrolled in January 2021. Households were randomly selected from our enumeration list for screening and enrolled if they met the following criteria: 1) having at least two members aged 5 years or younger; 2) no more than 7 permanent residents currently residing; 3) no plans to move from the study catchment area in the next 2 years; and 4) willingness to participate in entomological surveillance studies. In March 2022, 4 households (3 from Busia and 1 from Tororo away from the border) were enrolled to replace households that had been withdrawn from the cohort study. ### Screening, enrollment and follow-up of cohort participants All permanent residents from enrolled households were screened and enrolled in the cohort study if they met the following criteria: 1) the selected household was considered their primary residence; 2) agreement to come to the study clinic for any febrile illness and scheduled routine visits; 3) agreement to avoid antimalarial medications outside the study; and 4) provision of written informed consent (for parent or guardian in the case of children). The cohort was dynamic, such that over the course of the study any permanent residents that joined a household were screened for enrollment. At enrollment, a baseline evaluation was conducted including a detailed medical history, focused physical examination, and blood collection by venipuncture for hemoglobin measurement, thick blood smear, and storage for future molecular studies. A household survey was conducted to collect information on characteristics of the household and LLIN ownership; all household members were provided access to an LLIN after the survey. A wealth index was generated for each household using principal components analysis based on common assets and categorized into tertiles. Cohort study participants were encouraged to come to a dedicated study clinic open 7 days per week for all their medical care. Routine visits were conducted every 4 weeks and included a standardized evaluation and collection of blood by finger prick/heel stick (if *<* 6 months of age) or venipuncture (if aged 6 months and older) for thick blood smear, hemoglobin measurement (every 12 weeks), and storage for future molecular studies. Study participants found to have a fever (tympanic temperature *>* 38.0°C) or history of fever in the previous 24 hours at the time of any clinic visit had a thick blood smear read immediately. If the thick blood smear was positive by light microscopy, the patient was diagnosed with malaria and managed according to national guidelines. Study participants who missed their scheduled routine visits were visited at home and requested to come to the study clinic as soon as possible. All enrolled participants were followed through September 30th, 2023 unless they were prematurely withdrawn. Participants were withdrawn if they: 1) moved out of the cohort household; 2) were unable to be located for *>* 4 months; 3) withdrew informed consent; or 4) were unable to comply with the study schedule and procedures. ### Laboratory evaluations Thick blood smears were stained with 2% Giemsa for 30 minutes and evaluated for the presence of asexual parasites. Parasite densities were calculated by counting the number of asexual parasites per 200 leukocytes (or per 500, if the count was less than 10 parasites per 200 leukocytes), assuming a leukocyte count of 8,000/μL. A thick blood smear was considered negative if examination of 100 high power fields revealed no asexual parasites. For quality control, all slides were read by a second microscopist, and a third reviewer settled any discrepant readings. Quantitative PCR (qPCR) was performed at the time of enrollment, at each routine visit (every 4 weeks), and when malaria was diagnosed. At each of these time points, DNA was extracted from approximately 200 μL of whole blood using Qiagen spin columns, and extraction products were tested for the presence and quantity of *P. falciparum* DNA using a highly sensitive qPCR assay targeting the multicopy conserved var gene acidic terminal sequence with a lower limit of detection of 50 parasite/mL.(9) ### Entomological surveillance Mosquito collections were conducted every 2 weeks in all rooms of study houses where cohort study participants slept using CDC light traps positioned 1 m above the floor (Model 512; John W. Hock Company, Gainesville, Florida, USA). Traps were set at 7 PM and contents collected at 7 AM the following morning. Entomology technicians assessed whether cohort study participants reported sleeping under an LLIN each morning when trap contents were collected. All female *Anopheles* were enumerated and identified taxonomically to species level based on morphological criteria according to established taxonomic keys.(10) Every 2 weeks, up to 30 mosquitoes identified as from the *An. gambiae s.l.* complex from each of the 3 study sites were randomly selected for PCR analysis to distinguish *An. gambiae s.s.* from *An. arabiensis*.(11) All female Anopheles mosquitoes were stored in desiccant and up to 50 mosquitoes per CDC light trap collection were assessed for sporozoites using ELISA.(12) ### Insecticide resistance phenotyping Insecticide resistance testing was performed using field caught mosquitoes collected from neighboring Mayuge district to compare the response of both *An. funestus s.l.* and *An. Gambiae s.l.* to clothianidin and pirimiphos-methyl. Blood fed indoor-resting mosquitoes were collected in May 2023 using a Prokopack electric aspirator, kept for 3-4 days, and then made to lay eggs by forced-egg laying.(13) The F1 progeny were reared and then used for bioassays using CDC bottles for clothianidin and WHO tubes for pirimiphos-methyl at the standard diagnostic dose as per WHO protocol.(14) Mosquitoes were exposed at different times to assess the time response curve for both species and then kept for 5 days post exposure to clothianidin and 1 day post exposure to pirimiphos-methyl to assess delayed mortality. ### Data analyses All data were collected using standardized case record forms and double-entered using Microsoft Access (Microsoft Corporation, Redmond, Washington, USA). Analyses were performed using Stata, version 14 (Stata Corporation, College Station, Texas, USA). The primary objective of the study was to compare temporal trends in key malaria metrics in the context of population level vector control interventions stratified by 3 geographic regions: 1) Tororo, away from border, 2) Tororo, near border, and 3) Busia (Figure 1). The incidence of symptomatic malaria was defined as the number of incident cases of malaria (fever plus a positive thick blood smear by microscopy) divided by the person time of observation. Episodes of malaria occurring within 10 days of a prior episode were not considered incident events. The prevalence of microscopic parasitemia was defined as the proportion of routine visits conducted every 4 weeks in which asexual parasites were detected by microscopy. When estimating the prevalence of microscopic or sub-microscopic parasitemia, samples that were positive by qPCR but negative by microscopy were added to the numerator. Vector density was defined as the number of female *Anopheles* collected (stratified by species) divided by the number CDC light trap collections. The number of *An. gambiae s.s.* and *An. arabiensis* collected were estimated by multiplying the observed proportions using the subset tested with species PCR by the total number of *An. gambiae s.l.* collected stratified by each 2-week collection period and geographic region. Selected comparisons of measures of malaria incidence were made using mixed effects generalized linear models with a negative binomial regression family, with person time of observation included as an offset, and random effects at the level of the individual. Selected comparisons of measures of parasitemia at the time of routine visits were made using mixed effects generalized linear models with a Poisson family and random effects at the level of the individual. Selected comparisons of measures of vector density were made using negative binomial regression models with the number of CDC light trap collections included as an offset. For all comparisons, a p-value < 0.05 was considered statistically significant. ### Ethics approval and informed consent Ethical approval was obtained from the Makerere University School of Medicine Research and Ethics Committee (REF 2019–134), the Uganda National Council for Science and Technology (HS 2700), the London School of Hygiene & Tropical Medicine Ethics Committee (17777), and the University of California, San Francisco Committee on Human Research (257790). Written informed consent was obtained for all participants prior to enrolment into the cohort study. ## Results ### Temporal changes in malaria burden from Nagongera Health Center Data on temporal changes in malaria burden from Nagongera Health Center (located in the interior of Tororo District) covering from approximately 3 years before until 9 years after the initiation of IRS are presented in Figure 2. From October 2011 – November 2014, prior to the implementation of IRS, an average of 419 laboratory confirmed cases were diagnosed each month. After IRS was implemented, first with Bendiocarb (administered every 6 months) followed by Actellic (administered annually), there was a sharp decline in malaria burden, reaching an average of 56 cases per month and an incidence of 25 episodes per 1000 person years in the surrounding community from March 2019 – February 2020. After switching to annual rounds of IRS with clothianidin-based formulations in March 2020, malaria burden gradually began to increase. This increase accelerated after the 2nd round of Fludora Fusion in March 2021. From March 2022 – February 2023, following a round of SumiShield (the 3rd and final round of clothianidin-based IRS) malaria burden reached an average of 1591 cases per month and an incidence of 587 episodes per 1000 person years. Malaria burden declined sharply after switching back to IRS with Actellic in March 2023, reaching 155 cases and an incidence of 71 episodes per 1000 person years in September 2023. ### Characteristics of cohort households and study participants A cohort study including 3 distinct geographic regions in Tororo and Busia districts was designed to further characterize the resurgence of malaria and explore potential causes. Characteristics of the households and cohort participants at the time of enrollment are presented in Table 1. Household characteristics unrelated to vector control interventions were similar across the three geographic regions, with the exception of household wealth, as a higher proportion of houses in Busia were in the poorest category compared to houses in Tororo. The proportion of houses that reported owning at least one LLIN at enrollment ranged from 30% in Busia to 74% in Tororo, away from border. Adequate LLIN coverage (defined as 1 LLIN per 2 persons) was 19% and similarly low across the three geographic regions. All cohort households were provided LLINs by the study team following enrollment and 90% of household residents reported sleeping under an LLIN at the time of home assessments conducted every 2 weeks over the course of follow-up, with no significant differences between the three geographic regions. As expected, no households in Busia reported receiving IRS 12 months prior to enrollment or during the course of the study. Household coverage for the 4 annual rounds of IRS conducted in Tororo district from 2020 (assessed at enrollment) through 2023 was 90-100% within Tororo, near border and 80-100% within Tororo, away from border. View this table: [Table 1.](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/T1) Table 1. Characteristics of households and cohort participants at enrolment ### Temporal changes in the incidence of malaria Monthly trends over the entire observation period in the incidence of malaria among cohort participants of all ages are presented in Figure 3A stratified by geographic region. To allow for direct comparisons over time and in relationship to rounds of IRS, malaria incidence (stratified by geographic region and age categories) was quantified over 6-month time periods from April 2021 through September 2023 (Table 2). Among residents of all ages from Tororo, away from border, the incidence of clinical malaria remained below 0.5 episodes per person year (PY) for each of the first 7 months of observation from February-August 2021 before there was a sharp increase, reaching over 4 episodes per PY in March 2022. The incidence remained above 2.2 episodes per PY for all months from April 2022 until Actellic was sprayed in March 2023, after which there was a sharp decline, reaching 0.2 episodes per PY in September 2023 (Figure 3A). Comparing 6-month time periods from April–September following annual rounds of IRS, the incidence of malaria increased over 8-fold between 2021 and 2022 (0.36 vs. 2.97, p<0.0001) corresponding to the second and third rounds of IRS with clothianidin-based formulations, then decreased almost 5-fold between 2022 and 2023 (2.97 vs. 0.70, p<0.0001) following the shift back to Actellic (Table 2). Among residents from Tororo, near border, the incidence of clinical malaria increased over 2.5-fold between April-September of 2021 and 2022 (0.93 vs. 2.40, p<0.0001), then decreased almost 2-fold between 2022 and 2023 (2.40 vs. 1.33, p<0.0001). Among residents of Busia, the incidence of malaria between April-September did not change significantly from 2021 to 2022 to 2023. Across the 3 geographic regions, the incidence of malaria was between 3 and 4 times higher in children compared to adults, and the temporal patterns described above were similar across different age categories (Table 2). Despite the high burden of malaria in this study, only 4 of 2,266 episodes (0.2%) met WHO criteria for severe malaria: 2 cases of severe anemia (Hb < 5.0 gm/dL) occurred in children with confirmed sickle cell disease, one of whom died; 1 case of severe anemia occurred in a 1 year old child who recovered; and 1 case of jaundice occurred in a 3 year old child who recovered. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/16/2024.03.15.24304352/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/F3) Figure 3: Monthly trends in the incidence of clinical malaria, prevalence of microscopic parasitemia, and prevalence of microscopic or sub-microscopic parasitemia stratified by geographic regions. Vertical bars indicate rounds of IRS (Tororo District only): yellow = Fludora Fusion, purple = SumiShield, and orange = Actellic. View this table: [Table 2.](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/T2) Table 2. Temporal changes in measures of malaria burden stratified by age categories and study site ### Temporal changes in parasite prevalence Temporal changes in parasite prevalence mirrored those seen for incidence. Among residents of all ages from Tororo, away from border, the prevalence of microscopic parasitemia remained below 6% and the prevalence of microscopic or sub-microscopic parasitemia remained below 18% for each of the first 7 months of observation (February-August 2021) before there was a sharp increase, reaching over 34% and 78%, respectively, in the year after the final round of clothianidin-based IRS was administered (Figure 3B and 3C). After IRS was changed back to Actellic in March 2023, there was a sharp decline in these measures, with parasitemia reaching 13% and 30%, respectively. Comparing 6-month time periods from April–September following annual rounds of IRS, the prevalence of microscopic parasitemia increased over 5-fold (5% vs. 29%, p<0.0001) and the prevalence of microscopic or sub-microscopic parasitemia increased over 4-fold (17% vs 70%, p<0.0001) between April-September of 2021 and April-September of 2022. After spraying with Actellic, the prevalence of microscopic parasitemia decreased by 30% (29% vs. 20%, p<0.0001) and the prevalence of microscopic or sub-microscopic parasitemia decreased by 39% (70% vs 43%, p<0.0001). Among residents from Tororo, near border, the prevalence of microscopic parasitemia increased by 62% (24% vs. 40%, p<0.0001) and the prevalence of microscopic or sub-microscopic parasitemia increased by 47% (52% vs. 76%, p<0.0001) between April-September of 2021 and April-September of 2022. After spraying with Actellic, the prevalence of microscopic parasitemia decreased by 13% (40% vs. 35%, p=0.04) and the prevalence of microscopic or sub-microscopic parasitemia decreased by 17% (76% vs. 63%, p<0.0001). Among residents from Busia, between April-September of 2021 and April-September of 2022 there were modest increases in the prevalence of microscopic parasitemia (18% vs. 24%, p=0.03) and microscopic or sub-microscopic parasitemia (46% vs. 56%, p=0.03), and, in contrast to the areas in Tororo receiving IRS with Actellic, remained similarly high in April-September of 2023. Across the 3 geographic regions, the prevalence of parasitemia was highest among children 5-15 years of age and the temporal patterns described above were similar across different age categories. ### Temporal changes in entomological measures Among houses from Tororo, away from border, through December 2021 the predominant species was *An. arabiensis* (71% of female *Anopheles* collected) followed by *An. funestus* (14%) and *An. gambiae s.s.* (12%). After December 2021, *An. funestus* became the predominant species, a trend that continued until after IRS was shifted back to Actellic in March 2023, when the density of *An. funestus* decreased sharply (Figure 4A). Between mid-2021 and mid-2022, the density of *An. funestus* increased over 4-fold (0.2 vs 1.0 mosquitoes per collection, p<0.0001); this density then decreased over 10-fold (1.0 vs 0.08, p<0.0001) between mid-2022 and mid-2023. During these same 6-month time periods, the densities of *An. gambiae s.s.* and *An. arabiensis* decreased between 2021 and 2022, with no significant change from 2022 to 2023 (Table 3). Among houses from Tororo, near border, vector densities were initially more evenly distributed across the 3 major *Anopheles* species. After December 2021, *An. funestus* became the predominant species until after the round of IRS with Actellic in March 2023, although the decrease in *An. funestus* was not as sustained as in houses away from the border (Figure 4B, Table 3). In Busia, *An. gambiae s.s.* made up 73% of all female *Anopheles* collected and was the predominant species throughout the observation period, with the exception of June-August each year, when *An. funestus* was predominant (Figure 4C, Table 3). ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/16/2024.03.15.24304352/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/F4) Figure 4: Monthly trends in the vector densities of common *Anopheles* species stratified by geographic regions. Vertical bars indicate rounds of IRS: yellow = Fludora Fusion, purple = SumiShield, and orange = Actellic. View this table: [Table 3.](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/T3) Table 3. Temporal changes in vector density stratified by Anopheles species and study site Overall, sporozoite rates were approximately 50% higher among *An. funestus* (1.75%) compared to *An. arabiensis* (1.16%) and *An. gambiae s.s.* (1.12%), with sporozoite rates highest for *An. funestus* across all 3 geographic regions (Table 3). However, the study lacked precision to evaluate temporal trends in the proportion of mosquitoes infected with sporozoites. ### Insecticide resistance phenotyping Time response curves for *An. funestus s.l.* exposed to clothianidin and pirimiphos-methyl (active ingredient in Actellic) showed high resistance to clothianidin, but full susceptibility to pirimiphos-methyl (Figure 5). *An. gambiae s.l.* showed susceptibility to clothianidin even at reduced exposure times but mild tolerance to pirimiphos-methyl when exposure time was decreased. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/03/16/2024.03.15.24304352/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2024/03/16/2024.03.15.24304352/F5) Figure 5: Mortality *of An. gambiae s.l.* and *An. funestus* exposed to clothianidin and pirimiphos-methyl for varying durations. ## Discussion The combination of repeated universal LLIN distribution campaigns and sustained IRS has been associated with marked reductions in the burden of malaria at several historically holoendemic areas of Uganda, including Tororo District.(15) However, in 2020-21 malaria began to resurge in these areas following a change in IRS from an organophosphate to clothianidin-based formulations.(6) In this study, we quantified this resurgence and investigated potential causes in a cohort of residents of Tororo District, where IRS had been sustained since 2015, and neighboring Busia District, where IRS has never been implemented. In Tororo, resurgence of malaria increased approximately 18 months after the change to annual rounds of clothianidin-based IRS formulations, reaching a malaria incidence in children of over 4 episodes as compared to 0.5 episodes per year before the change, with over 75% of the population infected with malaria parasites. This resurgence was temporally correlated the emergence of *An. funestus* as the predominant malaria vector. In addition, *An. funestus* had a higher rate of infection with sporozoites compared to the other common malaria vectors and it exhibited high tolerance to clothianidin. Following a shift back to IRS with Actellic, the burden of malaria and the density of *An. funestus* declined sharply in Tororo. In contrast, in Busia district the burden of malaria remained high and relatively stable throughout this period, with *An. gambiae s.s.* the primary vector. This study and others provide several lines of evidence to support a causal link between the change from an organophosphate (Actellic) to clothianidin-based (Fludora Fusion/SumiShield) IRS formulations and the marked resurgence of malaria. First, the timing of the resurgence corresponded to a shift in active ingredient from pirimiphos-methyl to clothianidin-based formulations similar to what has been documented in 4 other districts of Uganda.(6) Second, the timing of the resurgence also corresponded with the emergence of *An. funestus*, which was found to have high tolerance to clothianidin and relatively high sporozoite rates, which likely contributed to a marked increase in transmission intensity. Third, the burden of malaria and density of *An. funestus* were sharply reduced when the formulation of IRS was shifted back to Actellic, a formulation for which *An. funestus* was found to be fully susceptible. Finally, just a few kilometers away in the neighboring district of Busia, malaria burden remained relatively unchanged, and the predominant vector was *An. gambiae s.s.*, consistent with historical reports from high transmission areas of Uganda.(16) Clothianidin is a neonicotinoid insecticide, one of 4 chemical classes prequalified for use in IRS by the WHO. The two most widely used clothianidin-based IRS formulations have been SumiShieldTM (active ingredient clothianidin alone) and Fludora FusionTM (active ingredients clothianidin combined with deltamethrin, a pyrethroid), which were prequalified by the WHO in 2017 and 2018, respectively. In preparation for and following the subsequent rollout of clothianidin-based IRS formulations in sub-Saharan Africa, several field-based studies demonstrated broad susceptibility of major Anopheles vector species to clothianidin, with residual activity lasting up to 18 months.(17) In a study of wild pyrethroid resistant *An. Gambiae s.l.* carried out in experimental huts in Benin from 2013-14, clothianidin caused high overall mortality rates across a range of wall surfaces.(18) In studies of insectary reared and wild collected mosquitoes conducted from 2016-17 across Africa, *An. gambiae s.l.* from 14 countries and *An. funestus* from 2 countries were widely susceptible to clothianidin.(19) Subsequent studies of field populations of common *Anopheles* species from Kenya, Benin, and Cote d’Ivoire also reported high mortality rates following exposure to clothianidin.(20, 21, 22, 23) It should be noted that clothianidin is slower acting than other commonly used insecticide classes, requiring a longer holding period to fully assess post-exposure mortality.(18, 24, 25) Despite these encouraging reports, some recent data have suggested the existence of reduced susceptibility to clothianidin among African vectors, as reported in this study. In a rural area of Cameroon where neonicotinoids have been used for crop protection, field caught *An. gambiae* showed resistance to clothianidin.(26) In a study of laboratory susceptible and resistant laboratory strains using a variety of wall surfaces, the potency of clothianidin was lower against *An. funestus* compared to *An. gambiae*.(17) In a study from Cameroon, Malawi, Ghana, and Uganda, field collected *An. funestus* were broadly susceptible to clothianidin, however, possible cross-resistance was detected in mosquitoes with genetic mutations associated with metabolic resistance.(27) Importantly, in the only other study of the clinical impact of clothianidin-based IRS formulations, from Northern Zambia (where *An. funestus* is the primary vector) in 2019-20, IRS with Fludora Fusion had no impact on parasite prevalence, demonstrating, as in the current study, a lack of clinically relevant efficacy of the insecticide.(28) One of the most striking findings in this and previous studies conducted in Tororo District was the dramatic shift in *Anopheles* species predominance following the implementation of different formulations of IRS. Prior to the implementation of IRS, the primary vector was *An. gambiae s.s.* (74% of female *Anopheles* collected), followed by *An. arabiensis* (22%) and *An. funestus* (4%). After 5 years of IRS with Bendiocarb followed by Actellic, vector densities declined dramatically with *An. arabiensis* making up 99% of female *Anopheles* collected.(29) Entomologic surveillance was not conducted in the area from November 2019 – August 2020, corresponding to the period when IRS with Fludora Fusion was first implemented. However, just prior to the second round of Fludora Fusion, *An. arabiensis* remained the primary vector in houses away from the border in Tororo with a mix of *An. arabiensis* and *An. gambiae s.s.* predominating in houses near the border in Tororo and *An. gambiae s.s.* predominating in in houses on the Busia side of the border. Although low levels of *An. funestus* were detected at baseline in all 3 geographic regions in this study, it was only just before the second round of clothianidin-based IRS that this species emerged as the predominant vector in Tororo, corresponding to the dramatic resurgence of malaria. Interestingly, in houses away from the border in Tororo, *An. arabiensis* again became the predominant species after IRS was shifted back to Actellic, consistent with the trend previously seen when Actellic was used. *An. funestus* has emerged as an increasingly important malaria vector in many parts of sub-Saharan Africa. In a systematic review of studies from east and southern Africa, the primary vector shifted from *An. gambiae s.l.* between 2000 and 2010 to *An. funestus* between 2011 and 2021.(30) In addition, several studies have reported higher infection rates among *An. funestus* relative to other *Anopheles* species.(30, 31, 32) *An. funestus* has rapidly developed insecticide resistance in many parts of Africa,(33, 34, 35, 36, 37) is highly anthropophilic, can survive longer than other vectors,(38) and may have the ability to take multiple blood meals within each gonotrophic cycle,(39) all factors that may have contributed to the increasing importance of this vector following the scale up vector control throughout Africa. Our prospective cohort study design allowed for a detailed description of the dynamics and magnitude of clinically relevant measures of the malaria resurgence. The rapid increase in the incidence of symptomatic malaria over a 6-9 month period was similar to historically high transmission areas of Uganda where IRS had been discontinued,(15, 40) suggesting that after the 2nd round of Fludora Fusion there was little or no impact, even in a setting of high reported utilization of LLINs. In addition, the peak incidence of malaria during the resurgence exceeded that in children enrolled in another cohort study prior to the implementation of IRS in Tororo District, an observation that was also documented using health facility-based data covering the period from 2011 through 2023 presented in Figure 1 of this study. Not only did the resurgence affect children, but malaria incidence also reached relatively high levels in adults. These findings suggest that the resurgence may have been compounded by a relative loss of naturally acquired immunity to malaria following an extended period when malaria had been well controlled. Despite this suggestion of loss of immunity, there was no evidence from our cohort that there was a “excess” in the risk of severe forms of malaria, highlighting the importance of prompt and effective treatment available 7 days a week at our study clinic. However, a nearby hospital based study conducted during the malaria resurgence reported severe malaria affecting older children with unusual presentations including high risks of prostration, jaundice, severe anemia, and black water fever.(41) This study was not without limitations. First, we used an observational study design, with measures of impact based on comparisons of temporal trends between districts; a cluster randomized controlled trial would have better compared trends, but available resources did not allow this design. Second, there may have been other contributing factors to the resurgence that we did not fully explore, including the dynamics of host immunity. Third, we may have underestimated the overall burden of disease in the community during periods of resurgence, since our cohorts receive prompt treatment and close follow up. Fourth, entomological data were limited to CDC light traps, which do not provide direct measures of biting rates and information on the timing and location of mosquito-to-human transmission. Fifth, measures of insecticide susceptibility did not utilize mosquitoes collected in the study districts; although they were from a nearby district, susceptibility may have differed from that in Tororo District. Despite these limitations, this study benefited from the availability of prior comprehensive surveillance data collected from the area over an extended period, resulting in important scientific findings with significant policy implications. ## Conclusion Our findings show that resurgence of malaria in Tororo was temporally and geographically correlated with a change in IRS formulation from Actellic to Fludora Fusion and the unique emergence of *An. funestus* as the predominant vector. This unprecedented increase in malaria burden despite sustained and intense vector control combining repeated rounds of universal LLIN distribution and IRS underscores the fragile nature of malaria control in Africa. Rotating IRS formulations has been promoted as a way to limit the selection of insecticide resistance, but in this case such a rotation was associated with a dramatic malaria resurgence in Uganda. Future changes in insecticides should take into consideration on-going surveillance that is comprehensive and includes clinical and entomological metrics, including vector species composition and insecticide resistance patterns. Surveillance data should also be timely and actionable, allowing for well-coordinated responses by policy makers. Indeed, the WHO Global Malaria Programme now encourages the use of multiple sources of data to better understand malaria risk at the sub-national level and to target interventions as part of its “High Burden High Impact” strategy.(42) ## Availability of data and materials Data from the cohort study and health facility-based surveillance will be made available through an established open-access clinical epidemiology database resource, ClinEpiDB. [https://clinepidb.org/ce/app](https://clinepidb.org/ce/app) ## Competing interests None ## Funding This project was funded by the US National Institutes of Health as part of the International Centers of Excellence in Malaria Research (ICEMR) program (U19AI089674). The funders played no role in the design of the study; in the collection, analyses, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication. ## Author’s contributions **Conceptualization:** Joaniter I. Nankabirwa, Emmanuel Arinaitwe, Moses R. Kamya, Grant Dorsey. **Data curation:** John Rek, Maato Zedi, Ambrose Oruni, Peter Olwoch, Jessica Briggs, Paul Krezanoski, Bryan Greenhouse, Grant Dorsey. **Formal analysis:** Grant Dorsey. **Funding acquisition:** Moses R. Kamya, Joaniter I. Nankabirwa, Grant Dorsey. **Methodology:** Joaniter I. Nankabirwa, John Rek, Bryan Greenhouse, Ambrose Oruni, Moses R. Kamya, Grant Dorsey. **Supervision:** Joaniter I. Nankabirwa, Emmanuel Arinaitwe, Moses R. Kamya, Grant Dorsey. **Writing – original draft:** Moses R. Kamya, Joaniter I. Nankabirwa, Emmanuel Arinaitwe, Ambrose Oruni, Grant Dorsey. **Writing – review & editing:** all authors ## Acknowledgements We thank all the study team members for successfully conducting the PRISM studies over the years and the Tororo and Busia district administrations for their support. We are grateful to the study participants who participated in this study and their families. * Received March 15, 2024. * Revision received March 15, 2024. * Accepted March 16, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526(7572):207-11. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature15535&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26375008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 2. 2.World Health Organization. World malaria report 2023. Geneva: World Health Organization; 2023. Licence: CC BY-NC-SA 3.0 IGO. 3. 3.World Health Organization. Global technical strategy for malaria 2016–2030, 2021 update. Geneva: World Health Organization; 2021. Licence: CC BY-NC-SA 3.0 IGO. 4. 4.Nankabirwa JI, Rek J, Arinaitwe E, Namuganga JF, Nsobya SL, Asua V, et al. East Africa International Center of Excellence for Malaria Research: Summary of Key Research Findings. Am J Trop Med Hyg. 2022;107(4_Suppl):21–32. 5. 5.Nankabirwa JI, Arinaitwe E, Rek J, Kilama M, Kizza T, Staedke SG, et al. Malaria Transmission, Infection, and Disease following Sustained Indoor Residual Spraying of Insecticide in Tororo, Uganda. Am J Trop Med Hyg. 2020;103(4):1525–33. 6. 6.Epstein A, Maiteki-Sebuguzi C, Namuganga JF, Nankabirwa JI, Gonahasa S, Opigo J, et al. Resurgence of malaria in Uganda despite sustained indoor residual spraying and repeated long lasting insecticidal net distributions. PLoS Global Public Health. 2022:doi:10.1371/journal.pgph.0000676. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgph.0000676&link_type=DOI) 7. 7.Kamya MR, Arinaitwe E, Wanzira H, Katureebe A, Barusya C, Kigozi SP, et al. Malaria transmission, infection, and disease at three sites with varied transmission intensity in Uganda: implications for malaria control. Am J Trop Med Hyg. 2015;92(5):903–12. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI5Mi81LzkwMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzE2LzIwMjQuMDMuMTUuMjQzMDQzNTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 8. 8.Nankabirwa JI, Bousema T, Blanken SL, Rek J, Arinaitwe E, Greenhouse B, et al. Measures of malaria transmission, infection, and disease in an area bordering two districts with and without sustained indoor residual spraying of insecticide in Uganda. PLoS ONE. 2022;17(12):e0279464. 9. 9.Hofmann N, Mwingira F, Shekalaghe S, Robinson LJ, Mueller I, Felger I. Ultra-sensitive detection of Plasmodium falciparum by amplification of multi-copy subtelomeric targets. PLoS Med. 2015;12(3):e1001788. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001788&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25734259&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 10. 10.Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa South of the Sahara. Johannesburg: The South African Institute for Medical Research; 1987. 11. 11.Scott JA, Brogdon WG, Collins FH. Identification of single specimens of the Anopheles gambiae complex by the polymerase chain reaction. Am J Trop Med Hyg. 1993;49(4):520–9. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI0OS80LzUyMCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzAzLzE2LzIwMjQuMDMuMTUuMjQzMDQzNTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 12. 12.Wirtz RA, Duncan JF, Njelesani EK, Schneider I, Brown AE, Oster CN, et al. ELISA method for detecting Plasmodium falciparum circumsporozoite antibody. Bull World Health Organ. 1989;67(5):535–42. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2692869&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989CF53600009&link_type=ISI) 13. 13.Morgan JC, Irving H, Okedi LM, Steven A, Wondji CS. Pyrethroid resistance in an Anopheles funestus population from Uganda. PLoS ONE. 2010;5(7):e11872. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0011872&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20686697&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 14. 14.World Health Organization. Manual for monitoring insecticide resistance in mosquito vectors and selecting appropriate interventions. Geneva: World Health Organization; 2022. Licence: CC BY-NC-SA 3.0 IGO. 15. 15.Namuganga JF, Epstein A, Nankabirwa JI, Mpimbaza A, Kiggundu M, Sserwanga A, et al. The impact of stopping and starting indoor residual spraying on malaria burden in Uganda. Nat Commun. 2021;12(1):2635. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-22896-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 16. 16.Yeka A, Gasasira A, Mpimbaza A, Achan J, Nankabirwa J, Nsobya S, et al. Malaria in Uganda: challenges to control on the long road to elimination: I. Epidemiology and current control efforts. Acta Trop. 2012;121(3):184–95. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actatropica.2011.03.004&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21420377&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000302668500005&link_type=ISI) 17. 17.Lees RS, Praulins G, Lissenden N, South A, Carson J, Brown F, et al. The Residual Efficacy of SumiShield 50WG and K-Othrine((R)) WG250 IRS Formulations Applied to Different Building Materials against Anopheles and Aedes Mosquitoes. Insects. 2022;13(2). 18. 18.Ngufor C, Fongnikin A, Rowland M, N’Guessan R. Indoor residual spraying with a mixture of clothianidin (a neonicotinoid insecticide) and deltamethrin provides improved control and long residual activity against pyrethroid resistant Anopheles gambiae sl in Southern Benin. PLoS ONE. 2017;12(12):e0189575. 19. 19.Oxborough RM, Seyoum A, Yihdego Y, Dabire R, Gnanguenon V, Wat’senga F, et al. Susceptibility testing of Anopheles malaria vectors with the neonicotinoid insecticide clothianidin; results from 16 African countries, in preparation for indoor residual spraying with new insecticide formulations. Malar J. 2019;18(1):264. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-019-2888-6&link_type=DOI) 20. 20.Agumba S, Gimnig JE, Ogonda L, Ombok M, Kosgei J, Munga S, et al. Diagnostic dose determination and efficacy of chlorfenapyr and clothianidin insecticides against Anopheles malaria vector populations of western Kenya. Malar J. 2019;18(1):243. 21. 21.Fongnikin A, Houeto N, Agbevo A, Odjo A, Syme T, N’Guessan R, et al. Efficacy of Fludora(R) Fusion (a mixture of deltamethrin and clothianidin) for indoor residual spraying against pyrethroid-resistant malaria vectors: laboratory and experimental hut evaluation. Parasit Vectors. 2020;13(1):466. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-020-04341-6&link_type=DOI) 22. 22.Kouassi BL, Edi C, Ouattara AF, Ekra AK, Bellai LG, Gouamene J, et al. Entomological monitoring data driving decision-making for appropriate and sustainable malaria vector control in Cote d’Ivoire. Malar J. 2023;22(1):14. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-023-04439-z&link_type=DOI) 23. 23.Odjo EM, Salako AS, Padonou GG, Yovogan B, Adoha CJ, Adjottin B, et al. What can be learned from the residual efficacy of three formulations of insecticides (pirimiphos-methyl, clothianidin and deltamethrin mixture, and clothianidin alone) in large-scale in community trial in North Benin, West Africa? Malar J. 2023;22(1):150. 24. 24.Agossa FR, Padonou GG, Koukpo CZ, Zola-Sahossi J, Azondekon R, Akuoko OK, et al. Efficacy of a novel mode of action of an indoor residual spraying product, SumiShield(R) 50WG against susceptible and resistant populations of Anopheles gambiae (s.l.) in Benin, West Africa. Parasit Vectors. 2018;11(1):293. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-018-2869-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29747684&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 25. 25.Marti-Soler H, Maquina M, Opiyo M, Alafo C, Sherrard-Smith E, Malheia A, et al. Effect of wall type, delayed mortality and mosquito age on the residual efficacy of a clothianidin-based indoor residual spray formulation (SumiShield 50WG) in southern Mozambique. PLoS ONE. 2021;16(8):e0248604. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0248604&link_type=DOI) 26. 26.Fouet C, Ashu FA, Ambadiang MM, Tchapga W, Wondji CS, Kamdem C. Clothianidin-resistant Anopheles gambiae adult mosquitoes from Yaounde, Cameroon, display reduced susceptibility to SumiShield(R) 50WG, a neonicotinoid formulation for indoor residual spraying. BMC Infect Dis. 2024;24(1):133. 27. 27.Assatse T, Tchouakui M, Mugenzi L, Menze B, Nguiffo-Nguete D, Tchapga W, et al. Anopheles funestus Populations across Africa Are Broadly Susceptible to Neonicotinoids but with Signals of Possible Cross-Resistance from the GSTe2 Gene. Trop Med Infect Dis. 2023;8(5). 28. 28.Ferriss E, Chaponda M, Muleba M, Kabuya JB, Lupiya JS, Riley C, et al. The Impact of Household and Community Indoor Residual Spray Coverage with Fludora Fusion in a High Malaria Transmission Setting in Northern Zambia. Am J Trop Med Hyg. 2023;109(2):248–57. 29. 29.Musiime AK, Smith DL, Kilama M, Rek J, Arinaitwe E, Nankabirwa JI, et al. Impact of vector control interventions on malaria transmission intensity, outdoor vector biting rates and Anopheles mosquito species composition in Tororo, Uganda. Malar J. 2019;18(1):445. 30. 30.Msugupakulya BJ, Urio NH, Jumanne M, Ngowo HS, Selvaraj P, Okumu FO, et al. Changes in contributions of different Anopheles vector species to malaria transmission in east and southern Africa from 2000 to 2022. Parasit Vectors. 2023;16(1):408. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-023-06019-1&link_type=DOI) 31. 31.Akoton R, Sovegnon PM, Djihinto OY, Medjigbodo AA, Agonhossou R, Saizonou HM, et al. Vectorial competence, insecticide resistance in Anopheles funestus and operational implications for malaria vector control strategies in Benin Republic. Malar J. 2023;22(1):385. 32. 32.Aschale Y, Getachew A, Yewhalaw D, De Cristofaro A, Sciarretta A, Atenafu G. Systematic review of sporozoite infection rate of Anopheles mosquitoes in Ethiopia, 2001-2021. Parasit Vectors. 2023;16(1):437. 33. 33.Choi KS, Christian R, Nardini L, Wood OR, Agubuzo E, Muleba M, et al. Insecticide resistance and role in malaria transmission of Anopheles funestus populations from Zambia and Zimbabwe. Parasit Vectors. 2014;7:464. 34. 34.Djouaka R, Riveron JM, Yessoufou A, Tchigossou G, Akoton R, Irving H, et al. Multiple insecticide resistance in an infected population of the malaria vector Anopheles funestus in Benin. Parasit Vectors. 2016;9:453. 35. 35.Menze BD, Riveron JM, Ibrahim SS, Irving H, Antonio-Nkondjio C, Awono-Ambene PH, et al. Multiple Insecticide Resistance in the Malaria Vector Anopheles funestus from Northern Cameroon Is Mediated by Metabolic Resistance Alongside Potential Target Site Insensitivity Mutations. PLoS ONE. 2016;11(10):e0163261. 36. 36.Mzilahowa T, Chiumia M, Mbewe RB, Uzalili VT, Luka-Banda M, Kutengule A, et al. Increasing insecticide resistance in Anopheles funestus and Anopheles arabiensis in Malawi, 2011-2015. Malar J. 2016;15(1):563. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27876046&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 37. 37.Pinda PG, Eichenberger C, Ngowo HS, Msaky DS, Abbasi S, Kihonda J, et al. Comparative assessment of insecticide resistance phenotypes in two major malaria vectors, Anopheles funestus and Anopheles arabiensis in south-eastern Tanzania. Malar J. 2020;19(1):408. 38. 38.Takken W, Verhulst NO. Host preferences of blood-feeding mosquitoes. Annu Rev Entomol. 2013;58:433–53. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-ento-120811-153618&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23020619&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316406000023&link_type=ISI) 39. 39.Omondi S, Kosgei J, Musula G, Muchoki M, Abong’o B, Agumba S, et al. Late morning biting behaviour of Anopheles funestus is a risk factor for transmission in schools in Siaya, western Kenya. Malar J. 2023;22(1):366. 40. 40.Raouf S, Mpimbaza A, Kigozi R, Sserwanga A, Rubahika D, Katamba H, et al. Resurgence of Malaria Following Discontinuation of Indoor Residual Spraying of Insecticide in an Area of Uganda With Previously High-Transmission Intensity. Clin Infect Dis. 2017;65(3):453–60. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/cix251&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F03%2F16%2F2024.03.15.24304352.atom) 41. 41.Namayanja C, Eregu EEI, Ongodia P, Okalebo CB, Okiror W, Okello F, et al. Unusual clinical spectra of childhood severe malaria during malaria epidemic in eastern Uganda: a prospective study. Malar J. 2023;22(1):169. 42. 42.World Health Organization. Malaria surveillance, monitoring & evaluation: a reference manual. Geneva: World Health Organization; 2018. Licence: CC BY-NC-SA 3.0 IGO..