Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Measures of malaria transmission, infection, and disease in an area bordering two districts with and without sustained indoor residual spraying of insecticide in Uganda

  • Joaniter I. Nankabirwa ,

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    jnankabirwa@yahoo.co.uk.

    Affiliations Infectious Diseases Research Collaboration, Kampala, Uganda, Department of Medicine, Makerere University College of Health Sciences, Kampala, Uganda

  • Teun Bousema,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliations Department of Medical Microbiology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands, Department of Immunology and Infection, London School of Hygiene and Tropical Medicine, London, United Kingdom

  • Sara Lynn Blanken,

    Roles Data curation, Formal analysis, Writing – review & editing

    Affiliation Department of Medical Microbiology, Radboud University Nijmegen Medical Centre, Nijmegen, The Netherlands

  • John Rek,

    Roles Data curation, Methodology, Writing – review & editing

    Affiliation Infectious Diseases Research Collaboration, Kampala, Uganda

  • Emmanuel Arinaitwe,

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation Infectious Diseases Research Collaboration, Kampala, Uganda

  • Bryan Greenhouse,

    Roles Conceptualization, Methodology, Writing – review & editing

    Affiliation Department of Medicine, University of California San Francisco, San Francisco, CA, United States of America

  • Philip J. Rosenthal,

    Roles Methodology, Writing – review & editing

    Affiliation Department of Medicine, University of California San Francisco, San Francisco, CA, United States of America

  • Moses R. Kamya,

    Roles Conceptualization, Methodology, Supervision, Writing – review & editing

    Affiliations Infectious Diseases Research Collaboration, Kampala, Uganda, Department of Medicine, Makerere University College of Health Sciences, Kampala, Uganda

  • Sarah G. Staedke,

    Roles Data curation, Writing – review & editing

    Affiliation Department of Infectious Diseases, London School of Hygiene and Tropical Medicine, London, United Kingdom

  • Grant Dorsey

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Department of Medicine, University of California San Francisco, San Francisco, CA, United States of America

Abstract

Tororo District, in Eastern Uganda, experienced a dramatic decline in malaria burden starting in 2014 following the implementation of indoor residual spraying of insecticide (IRS) in the setting of repeated long-lasting insecticide treated nets (LLINs) distribution campaigns. However, in 2020 malaria began to resurge in Tororo following a change in the active ingredient used for IRS. In this study, epidemiological measures of malaria were compared shortly after the resurgence between two parishes in Tororo District (Kayoro and Osukuru) and one contiguous parish in Busia District (Buteba), where IRS has never been implemented. A cohort of 483 residents from 80 randomly selected households were followed from August 2020 to January 2021. Mosquitoes were collected every 2 weeks using CDC light traps in rooms where participants slept; parasitemia and gametoctyemia measured every 4 weeks by microscopy and PCR; and symptomatic malaria measured by passive surveillance. The annual entomological inoculation rate was significantly higher in Buteba (108.2 infective bites/person/year), compared to Osukuru (59.0, p = 0.001) and Kayoro (27.4, p<0.001). Overall, parasite prevalence was 19.5% by microscopy and 50.7% by PCR, with no significant differences between the three parishes. Among infected individuals, gametocyte prevalence by PCR was 45.5% and similar between sites. The incidence of malaria was significantly higher in Osukuru (2.46 episodes PPY) compared to Buteba (1.47, p = 0.005) and Kayoro (1.09, p<0.001). For participants over 15 years of age, the risk of symptomatic malaria if microscopic parasitemia was present was higher in Osukuru (relative risk [RR] = 2.99, p = 0.03) compared to Buteba. These findings highlight the complex relationships between measures of malaria transmission, infection, and disease, and the potential for excess disease burden, possibly due to waning immunity, in areas where vector control interventions begin to fail after a sustained period of highly effective control.

Introduction

Significant progress in malaria control has been realized in the last decade in Uganda [1]. This has been attributed to the scale up of effective control interventions including case management with artemisinin-based combination therapies, mass distribution of long-lasting insecticide treated nets (LLINs), intermittent preventive treatment in pregnancy (IPTp), and re-introduction of indoor residual spraying (IRS) in some districts. Together, these interventions have resulted in declines in measures of transmission intensity and parasite prevalence in children under 5 years of age [13], although this progress has not been unform across the country and has stalled in the last few years [4].

One dramatic example of progress is Nagongera sub-county, located in the central part of Tororo District in Eastern Uganda, a historically high transmission area. Comparing key malaria indicators measured from 2011–14, prior to IRS, and in 2017–19 after two national mass LLIN distribution campaigns and 5 years of sustained IRS, dramatic reductions in estimates were observed, with the entomological inoculation rate (EIR) falling from 238 to 0.43 infectious bites/person/year, prevalence of parasitemia measured by microscopy in children aged 0.5–10 years from 32% to 6.8%, and incidence of malaria from 3 to 0.05 episodes/person/year [2]. Although this dramatic decline in malaria burden following intense and sustained vector control suggests that elimination may be feasible in a previously hyperendemic area of Uganda, more recent data highlight the fragile nature of these gains. Beginning in 2020, a marked resurgence in malaria cases was documented in 5 districts of Uganda (including Tororo) following a change in IRS formulation, reaching pre-IRS levels within 1–2 years [5]. Although the causes of this resurgence are unclear, it is important to better understand the underlying epidemiology of malaria in areas where resurgence has occurred. In this study, we compared measures of malaria transmission, infection, and disease in a contiguous border area that included two parishes in Tororo District where control effects began to fail after 6 years of sustained IRS and one parish in Busia District, where IRS has never been implemented. Our underlying hypothesis was that there would be important difference in epidemiological measures of malaria between these two regions that could further our understanding of the epidemiological consequences of resurgent malaria in area where vector control interventions have diminished effectiveness.

Methods

Study design and setting

Between August 2020 and January 2021, a comprehensive cohort study and entomological surveillance were conducted in a contiguous area of Eastern Uganda that included two parishes in Tororo District (Osukuru and Kayoro parishes) and one parish in Busia District (Buteba parish) (Fig 1). To provide contemporaneous data from the interior of Tororo District (Nagongera Health Centre) and within the area of Tororo District where the cohort study was conducted (Osukuru Health Centre), routine health facility-based data on malaria morbidity was collected from January 2015 through June 2021.

thumbnail
Fig 1. Map of Uganda with an extract showing study sites.

https://doi.org/10.1371/journal.pone.0279464.g001

Before 2013, malaria control in Tororo District was limited to the distribution of LLINs through antenatal care services, promotion of intermittent preventive treatment during pregnancy, and malaria case management with artemether-lumefantrine. 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. Indoor residual spraying with the carbamate bendiocarb was initiated in December 2014–January 2015, with additional rounds administered in June–July 2015 and November–December 2015. In June–July 2016, the IRS formulation was changed to the organophosphate pirimiphos-methyl (Actellic), with repeated rounds in June–July 2017, June–July 2018, and March–April 2019. In March-April 2020, the formulation of IRS was changed to Fludora Fusion, a recently prequalified insecticide containing a mixture of clothianidin and deltamethrin [6]. IRS has also been sustained since December 2014 in Butaleja and Bugiri Districts, which border the northern side of Tororo District, but has not been implemented in other districts bordering Tororo or on the Kenyan side of the border with Tororo.

Busia District is a rural district bordering Tororo district to the south. Although no information on EIR is available for Busia District, data from health facility-based surveillance have documented very high-test positivity rates (TPRs) among patients presenting to outpatient facilities in this area [7]. As in Tororo, prior to 2013 malaria control in Busia District was limited to targeted distribution of LLINs through antenatal care services, promotion of intermittent preventive treatment during pregnancy, and malaria case management with AL. Universal distribution of free LLINs was conducted as part of repeated national campaigns in Busia District in 2013, 2017, and in December 2020. IRS has never been implemented in Busia District.

Routine health-facility based malaria surveillance data

Weekly aggregate data on the number of patients tested for malaria and the number of patients testing positive for malaria were collected retrospectively from the outpatient departments of Nagongera Health Centre IV and Osukuru Health Centre III. Temporal trends in malaria morbidity were generated by calculating the weekly test positivity rate, defined as the number of patients testing positive for malaria / the number of patients tested.

Screening and enrollment of cohort households

In May 2020, all 10,474 households in the parishes of Osukuru, Kayoro, and Buteba 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 (Fig 1). A household was defined as any single permanent or semi-permanent dwelling structure acting as the primary residence for a person or group of people that generally cook and eat together. Following enumeration, a cross-sectional survey was conducted in 300 randomly selected households to generate spatial estimates of parasite prevalence within the study area. Using the survey results, the study area was stratified into high, medium, and low transmission areas. In August 2020, households in the three transmission areas were randomly selected and screened for eligibility to join the cohort study. Households were enrolled in the cohort 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.

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 which included 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. 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 [8]. 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 January 10th, 2021 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.

Entomological surveillance

Mosquito collections were conducted every 2 weeks in all cohort study households. In each room where study participants slept, a miniature CDC light trap (Model 512; John W. Hock Company, Gainesville, Florida, USA) was positioned 1 m above the floor. Traps were set at 7 PM and collected at 7 AM the following morning. Female Anopheles were identified taxonomically to species level based on morphological criteria according to established taxonomic keys [9]. Every 2 weeks, up to 30 mosquitoes were randomly selected for PCR analysis to identify members of the Anopheles gambiae complex [10]. All female Anopheles mosquitoes were stored in dessicant and assessed for sporozoites using a standardized ELISA technique [11].

Laboratory evaluations

Thick blood smears were stained with 2% Giemsa for 30 minutes and evaluated for the presence of asexual and sexual (gametocytes) parasites. Parasite and gametocyte 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 [12].

Following automated extraction (Total Nucleic Acid Isolation Kit-High Performance; Roche Applied Science, Indianapolis, IN, USA), quantitative reverse transcriptase PCR (RT-qPCR) targeting male (PfMGET) and female (CCp4) mRNA transcripts was used to quantify gametocytes using 100 μL whole blood in RNA preservative (RNA protect Cell Reagent; Qiagen, Hilden, Germany) with a lower limit of detection of 0.01 gametocytes/ μL [13]. For high parasitemia samples (>1000 parasites/ μL), gametocyte densities were adjusted for background transcripts in asexual parasites [13]. Of all qPCR positive samples (n = 1912), the first qPCR positive follow-up sample from each participant (n = 387) was used for gametocyte quantification. Due to 6 extraction failures, 4 of the first qPCR positive follow-up visits were replaced by the second qPCR positive follow-up visit. This resulted in a total of 385 gametocyte measurements.

Ethics approval and consent to participate

Ethical approval was obtained from the Makerere University School of Medicine Research and Ethics Committee (REF 2019–134), the Uganda National Council of 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 study.

Statistical analysis

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) and R version 4.1.0 [14]. Baseline descriptive statistics included proportions for categorical variables and mean (SD) or median (range) values for continuous variables. Measures of transmission were based on entomological surveillance data. The human biting rate (HBR) was estimated using the total number of female Anopheles mosquitoes captured / number of CDC light trap collections. The sporozoite rate was calculated as the number of mosquitoes testing positive for sporozoites / the number of mosquitoes tested. The annual EIR was estimated using the product of the daily HBR and the sporozoite rate x 365 days/year.

Clinical metrics, including measures of infection and disease, were estimated and stratified by parish of residence and age categories. The prevalence of microscopic parasitemia/gametocytemia was calculated as the number of routine blood smears positive for asexual parasites or sexual parasites / total number of routine blood smears done. The prevalence of microscopic or sub-microscopic parasitemia was calculated as the number of routine samples positive by microscopy or qPCR / total number of routine assessments. The prevalence of gametocytes by RT-qPCR was calculated as the percentage of gametocyte positive samples among individuals at the time of their first follow-up visit when parasites were detected by qPCR. The incidence of malaria was calculated as the number of incident episodes of symptomatic malaria (defined as fever and a positive thick blood smear not proceeded by another episode of malaria in the prior 14 days) / person years of observation. To assess for measures of clinical immunity to malaria, individual observation time was divided into 28-day intervals (corresponding to the timing of routine clinic visits), and the interval risks of any symptomatic malaria, microscopic parasitemia, and symptomatic malaria if microscopic parasitemia was present were measured.

Pairwise comparisons between the parishes of entomological measures of transmission (HBR, sporozoite rate, and EIR) at the household level were made using Poisson regression models. Pairwise comparisons between the parishes of measures of infection (parasite and gametocyte prevalence using both microscopy and PCR) and interval risks of clinical outcomes were made using mixed effects generalized linear models with a Poisson family and adjustment for repeated measures from the same individual. Pairwise comparisons of malaria incidence between the parishes were made using negative binomial regression models. Assuming a p-value of < 0.05 (two-sided) was considered statistically significant.

Results

Test positivity rates at health facilities

To provide context for findings from the cohort study presented below, trends in malaria TPRs over a 6.5-year period following the implementation of IRS were compared between Nagongera Health Center, located in the interior of Tororo District, and Osukuru Health Center, located near the eastern border of Tororo District, where the cohort study was conducted. At both health centers, a gradual decline in TPR was recorded from 2015–2018, with peaks proceeding each round of IRS (Fig 2); in Osukuru mean weekly TPR decreased from 45.5% in 2015 to 26.8% in 2018, and in Nagongera mean weekly TPR decreased from 15.4% in 2015 to 6.5% in 2018. At both facilities, TPR remained relatively stable in 2019 and most of 2020, and then began to rise in late 2020; in Osukuru mean weekly TPR reached 39.1% in the first half of 2021, in Nagongera mean weekly TPR reached 24.4% in the first half of 2021. Throughout the 6.5-year period following the implementation of IRS in Tororo, TPR was consistently higher in Osukuru (mean weekly TPR = 35.4%) compared to Nagongera (mean weekly TPR = 15.0%).

thumbnail
Fig 2. Temporal changes in test positivity rates estimated from routine health facility-based data (January 2015 –June 2021) from the interior of Tororo District (Nagongera Health Centre IV) and within the border area of Tororo District where the cohort study was conducted (Osukuru Health Centre III).

(Pink bars signify rounds of mass LLIN distribution, yellow bars are rounds of IRS with bendiocarb, green bars are rounds of IRS with pirimiphos-methyl, and blue bars are rounds of IRS with Fludora Fusion).

https://doi.org/10.1371/journal.pone.0279464.g002

Characteristics of the households and cohort study participants

A total of 326 households were screened and 80 enrolled in the study (20 from Buteba Parish, 40 from Kayoro Parish, and 20 from Osukuru Parish). The primary reason for exclusion was having fewer than 2 household residents under 5 years of age (Fig 3). Characteristics of the study households and cohort participants are presented in Table 1. Houses from Osukuru were more likely to be constructed using traditional materials and houses from Buteba were more likely to be in the poorest wealth category. As expected, no households in Buteba, 90% of households in Kayoro, and all households in Osukuru reported receiving IRS in the last 12 months.

thumbnail
Table 1. Characteristics of households and cohort participants at enrolment.

https://doi.org/10.1371/journal.pone.0279464.t001

A mass campaign to distribute LLINs was conducted 3 months prior to cohort enrollment in Tororo district, and 3 months after enrollment in Busia. At the time of enrollment, the proportion of households that reported owning at least one LLIN was 45% in Kayoro and 65% in Osukuru (Tororo district), and only 20% in Buteba (Busia district). All cohort households were provided LLINs by the study team at the time of enrollment. All 483 household residents were enrolled in the cohort study (119 from Buteba, 247 from Kayoro, and 117 from Osukuru). Of the enrolled participants, 262 (54.2%) were female, and the median age at enrollment was 8.6 years (interquartile range 3.5–25.7), with a similar distribution of age and gender across the three parishes (Table 1).

Measures of transmission

Entomological measures of transmission stratified by parish are presented in Table 2. Overall, 19,563 female Anopheles mosquitoes were collected in 1,190 overnight CDC light trap collections (Table 2). Anopheles species composition differed between the two districts. In Buteba, where IRS was not implemented, the majority (85.4%) of mosquitoes collected were Anopheles gambiae s.s. In Kayoro and Osukuru, where IRS has been sustained for an extended period, less than 50% of mosquitoes were Anopheles gambiae s.s., with Anopheles arabiensis and Anopheles funestus much more common than in Buteba. The daily HBR was significantly higher in Buteba (23.6) compared to Kayoro (13.1 p<0.001), and Osukuru (17.4, p<0.001). The annual EIR was significantly higher in Buteba (108.2 infective bites/person/year) compared to Osukuru (59.0, p = 0.002) and Kayoro (27.4, p<0.001).

thumbnail
Table 2. Entomological metrics from bi-weekly CDC light trap collections stratified by parish.

https://doi.org/10.1371/journal.pone.0279464.t002

Measures of infection and disease

Routine visits were conducted at enrollment and then every 4 weeks, resulting in 2,523 scheduled visits over the study period (Table 3). Overall, parasite prevalence was 19.5% by microscopy and 50.7% by qPCR, with no significant differences across the three parishes. When stratified by age, patterns of parasite prevalence by microscopy and qPCR were similar across the 3 parishes (Fig 4). For all parishes combined, parasite prevalence by microscopy was higher among those 5–15 years of age (34.7%) compared to those < 5 years (14.6%, p<0.001) and those over 15 years (12.3%, p<0.001), with no significant difference between those < 5 years and those > 15 years (p = 0.29). Parasite prevalence by qPCR was also higher among those 5–15 years (66.3%) compared to those < 5 years (39.3%, p<0.001) and those over 15 years (50.3%, p<0.001), but also higher in those over 15 years compared to those < 5 years (p = 0.002). Overall, the prevalence of gametocytes was 2.0% by microscopy and 45.5% by RT-qPCR among those with parasites detected by qPCR, with no significant differences across the three parishes.

thumbnail
Fig 4. Prevalence of microscopic parasitemia and microscopic+ submicroscopic parasitemia as well as incidence of malaria across the three parishes stratified by age groups.

All estimates (points) and 95% confidence intervals (lines) were calculated using generalized estimating equations with a binomial family and exchangeable correlation structure to account for repeated measures within individuals.

https://doi.org/10.1371/journal.pone.0279464.g004

thumbnail
Table 3. Measures of infection and disease stratified by parish.

https://doi.org/10.1371/journal.pone.0279464.t003

A total of 280 episodes of symptomatic malaria were diagnosed over 181 person years of follow-up (Table 3); 277 episodes were uncomplicated and treated with AL and 3 episodes occurred in children with danger signs who were treated with IV artesunate, followed by oral AL (2 from Buteba and 1 from Kayoro). No cases met criteria for severe malaria. The incidence of malaria was significantly higher in Osukuru (2.46 episodes/person/per year) compared to Buteba (1.47, p = 0.005) and Kayoro (1.09, p<0.001). When comparing the incidence of malaria across parishes stratified by age, the incidence of malaria was more than twice as high in Oskuru compared to Buteba and Kayoro for children < 5 years of age and 5–15 years of age (Fig 4). As expected, the incidence of malaria was lower in participants aged 15 years or older compared to what was observed in children and there were no significant differences in estimates between the 3 parishes for those aged 15 years of older. When comparing Buteba to Kayoro, the incidence of malaria was significantly different only among children < 5 years (2.28 vs. 1.41, p = 0.03).

To assess for differences in clinical immunity, interval risks of symptomatic malaria, microscopic parasitemia, and the joint probability of symptomatic malaria if microscopic parasitemia present were compared between the parishes after stratification by age groups (Table 4). Using Buteba (the parish with the highest transmission intensity) as the reference group, the relative risks for symptomatic malaria increased with increasing age for both Kayoro and Osukuru, although pairwise comparisons were only statistically significant for Osukuru among the < 5 year and 5–15 year age groups. There were no significant differences in the interval risk of microscopic parasitemia between Buteba and Kayoro or Osukuru, although in Osukuru relative risks decreased with increasing age, the opposite of the trend seen with symptomatic malaria. For the combined risk of symptomatic malaria if microscopic parasitemia present (an indicator of less clinical immunity), compared to Buteba the relative risks increased with increasing age for both Kayoro and Osukuru, with those over 15 years of age in Osukuru being almost 3 times more likely to have symptomatic malaria if microscopic parasitemia present compared to Buteba (RR = 2.99, 95% CI 1.15–7.78, p = 0.03) (Table 4).

thumbnail
Table 4. Associations between parish and clinical outcomes stratified by age.

https://doi.org/10.1371/journal.pone.0279464.t004

Discussion

Like many other countries in sub-Saharan Africa, malaria prevention in Uganda has largely relied on vector control through large scale LLIN distribution campaigns and IRS in selected areas. This approach has been highly successful, with an 84% reduction in malaria cases documented in five districts of Uganda 4–5 years after IRS was first implemented in 2014. However, more recently malaria cases have markedly increased despite sustained IRS, exceeding pre-IRS levels by 39% in these same five districts [5]. Although causes of this resurgence are not clear, the timing corresponded to a change in active ingredient used for IRS. This study compared the epidemiology of malaria in a contiguous area that included two parishes in Tororo District early in the resurgence despite sustained IRS with one parish in Busia District where LLINs were IRS has never been implemented. Entomologic measures of transmission intensity were modestly lower in the two parishes from Tororo District compared to the parish in Busia District. In contrast, parasite and gametocyte prevalence were similar across the three parishes and malaria incidence was highest in one of the parishes in Tororo, where there was evidence of waning immunity against symptomatic malaria in older individuals.

Confirming recent observations from routine health facility data [5], the burden of malaria in the studied parishes of Tororo were much higher compared to those measured a few years earlier in Nagongera sub-county, located in the interior of Tororo District. Indeed, using the same methodology as this study, malaria metrics in a cohort of all ages in Nagongera had reached pre-elimination levels after 5 years of sustained IRS, with an annual EIR of < 1 infectious bite/person/year, parasite prevalence by microscopy < 2%, and an incidence of malaria of 0.04 episodes/person/year [2].

There are a several potential explanations as to why the burden of malaria was much higher a few years later in the two parishes from Tororo District included in this study that bordered areas were IRS has not been implemented. First, repeated rounds of IRS may not be as effective or have a shorter duration of effect in areas that border high transmission areas that have not received as intensive vector control measures. This is supported by the health facility-based data included in this study, which consistently showed higher TPRs from the facility located near the border where IRS has not been implemented compared to the facility located in the interior of Tororo District. Movement of people and mosquitos infected with malaria parasites across borders has been well documented as one of the major challenges to malaria control and elimination efforts especially between countries [15, 16]. Indeed, there are several examples from Southern Africa where IRS has been highly effective in reducing the burden of malaria, however, achieving elimination has been challenged by the importation of malaria across borders [17].

Second, is the possibility that the effectiveness of IRS in Tororo District may have declined recently, corresponding to the interval period from when the previous study in Nagongera sub-county and this study were conducted. This is also supported by the heath facility-based data included in this study which showed an increase in TPR since late 2019 at both facilities. The dramatic decline in malaria burden documented in Nagongera sub-county after 5 years of sustained IRS involved round of the carbamate bendiocarb every 6 months (2014–2015) followed by annual rounds with the organophosphate Actellic (2016–2019) [2]. However, in March 2020 there was a switch to Fludora Fusion, a novel insecticide combination containing clothianidin and deltamethrin that was recently added to the WHO’s list of pre-qualified insecticides for use in indoor residual spraying [6]. Fludora Fusion has been shown to be effective in reducing measures of transmission intensity, especially in areas where pyrethroid-resistance has been reported [18, 19]. However, data on Fludora Fusion have been limited to laboratory-based or experimental hut studies and none have evaluated clinically relevant outcomes in “real world” settings. In addition, the recent effectiveness of the IRS program may have been compromised by problems with implementation, especially in the context of the recent COVID-19 pandemic. Although, households included in this study reported high coverage levels during the last round of IRS, other problems may have occurred such as spray operators spending less time in homes due to fear of acquiring COVID. Indeed, disruptions of essential malaria services, including IRS, may have played a role in a surge in malaria cases documented during in the first half of 2020 in Zimbabwe [20].

Third, species composition and behavioral characteristics of mosquitoes may have played a role in the relatively high burden of malaria in the Kayoro and Osukuru parishes. Prior to the implementation of IRS, the predominant species collected in Nangongera was Anopheles gambiae s.s.. However, after 5 years of sustained IRS, there was a dramatic shift in species composition, with Anopheles arabiensis making up over 98% of mosquitoes collected [2, 3]. Compared to Nagongera, species composition was quite different in Kayoro and Osukuru parishes, with Anopheles arabiensis making up only a third of mosquitoes collected, and the remainder primarily Anopheles gambiae s.s or Anopheles funestes. Another contributing factor could have been a shift from indoor to outdoor biting, circumventing IRS targeting indoor biting mosquitoes, which has been reported in Tororo District [21].

Another interesting observation from this study was that, compared to Buteba parish in Busia district, where IRS has never been implemented, Osukuru parish in Tororo district had a higher incidence of malaria, despite lower transmission intensity. In addition, older parasitemic individuals in Osukuru were more likely to develop symptomatic malaria, suggesting less clinical immunity. It is well known that, in high transmission areas, older children and adults gradually develop immunity against malaria infection and disease, although this immunity is incomplete and not permanent [22]. One possible explanation for the findings of this study is that immunity waned among the population of Tororo District during the years when IRS was highly effective, evidenced by the more modest decrease in risk of developing symptomatic malaria in the setting of microscopic parasitemia seen with increasing age at these sites. If relatively high transmission intensity in Osukuru parish was a recent phenomenon, due to loss of IRS efficacy or other factors, this could have resulted in an unusually high burden of symptomatic malaria in a population that had become relatively nonimmune. An analogous scenario has been described in an area of northern Uganda, where the withdrawal of IRS after an extended period was associated with a remarkably high malaria burden, even exceeding pre-IRS levels in patients over 5 years of age [23].

This study had several limitations. First, data from the three study areas represented only a “snapshot” in time, and care should be taken when comparing data collected previously in a different area of Tororo District. Second, entomological measures of transmission intensity were generated using CDC light traps and therefore only surrogate measures of human biting behavior that did not capture outdoor biting mosquitoes. Third, in addition to impact of IRS, which was the primary focus of this study, differences in malaria metrics between the 3 parishes could have been due to a number of other factors that we were unable to measure. Fourth, given the labor-intensive nature of this comprehensive cohort study, the sample size was relatively small, limiting statistical power to assess associations and make precise estimates generalizable to the larger population.

Conclusion

In a border area where IRS has been sustained for 6 years, measures of malaria transmission, infection and disease were comparable to those in an adjacent area where IRS had not been implemented and much higher than recent estimates from the interior of the district receiving IRS. These observations highlight the importance of detailed surveillance studies to better understand the scope and durability of malaria control interventions targeting specific geographic areas. Further studies are warranted to continue to monitor the effectiveness of IRS in Uganda, the causes of the recent decline in effectiveness, and implications for sustaining gains and consideration of IRS exit strategies.

Availability of data and materials

Data from both cohort studies are available through a novel open-access clinical epidemiology database resource, ClinEpiDB [24].

Acknowledgments

We thank the study team and the Infectious Diseases Research Collaboration (IDRC) for administrative and technical support. We are grateful to the study participants who participated in this study and their families.

References

  1. 1. Ministry of Health. 2018–19 Uganda Malaria Indicator Survey. Acccessed from: https://dhsprogramcom/pubs/pdf/MIS34/MIS34pdf on 18 March 2022. 2020.
  2. 2. 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. pmid:32700666
  3. 3. Mawejje HD, Kilama M, Kigozi SP, Musiime AK, Kamya M, Lines J, et al. Impact of seasonality and malaria control interventions on Anopheles density and species composition from three areas of Uganda with differing malaria endemicity. Malaria Journal. 2021;20(1):138. pmid:33678166
  4. 4. World Health Organisation. World Malaria Report 2021. Accessed from: https://wwwwhoint/teams/global-malaria-programme/reports/world-malaria-report-2021 on 8 February 2022. 2021.
  5. 5. Epstein A, Maiteki-Ssebuguzi 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. medRxiv. 2022:2022.02.18.
  6. 6. World Health Organisation. Prequalification Team Vector Control Decision Document Fludora Fusion. Accessed from: https://wwwwhoint/pq-vector-control/prequalified-lists/FinalDecisionFludoraFusionpdf on 23 September 2021. 2018.
  7. 7. Okiring J, Epstein A, Namuganga JF, Kamya V, Sserwanga A, Kapisi J, et al. Relationships between test positivity rate, total laboratory confirmed cases of malaria, and malaria incidence in high burden settings of Uganda: an ecological analysis. Malar J. 2021;20(1):42. pmid:33441121
  8. 8. Ministry of Health Uganda. Uganda Clinical Guidelines 2016. Accessible at: http://library.health.go.ug/publications/guidelines/uganda-clinical-guidelines-2016 accessed on 23 July 20202016.
  9. 9. Gillies MT, Coetzee M. A supplement to the Anophelinae of Africa South of the Sahara. Johannesburg. South African Institute for Medical Research. 1987.
  10. 10. 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. pmid:8214283
  11. 11. Wirtz RA, Duncan JF, Njelesani EK, Schneider I, Brown AE, Oster CN, et al. ELISA method for detecting Plasmodium falciparum circumsporozoite antibody. Bulletin of the World Health Organization. 1989;67(5):535–42. pmid:2692869
  12. 12. 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 medicine. 2015;12(3):e1001788. pmid:25734259
  13. 13. Meerstein-Kessel L, Andolina C, Carrio E, Mahamar A, Sawa P, Diawara H, et al. A multiplex assay for the sensitive detection and quantification of male and female Plasmodium falciparum gametocytes. Malaria Journal. 2018;17(1):441. pmid:30497508
  14. 14. R Core Team. A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Accessed from: https://wwwR-projectorg/ on 19th March 2022. 2021.
  15. 15. Wangdi K, Gatton ML, Kelly GC, Clements AC. Cross-border malaria: a major obstacle for malaria elimination. Adv Parasitol. 2015;89:79–107. pmid:26003036
  16. 16. Recht J, Siqueira AM, Monteiro WM, Herrera SM, Herrera S, Lacerda MVG. Malaria in Brazil, Colombia, Peru and Venezuela: current challenges in malaria control and elimination. Malaria Journal. 2017;16(1):273. pmid:28676055
  17. 17. Khosa E, Kuonza LR, Kruger P, Maimela E. Towards the elimination of malaria in South Africa: a review of surveillance data in Mutale Municipality, Limpopo Province, 2005 to 2010. Malar J. 2013;12:7. pmid:23294805
  18. 18. Agossa FR, Padonou GG, Fassinou A, Odjo EM, Akuoko OK, Salako A, et al. Small-scale field evaluation of the efficacy and residual effect of Fludora(®) Fusion (mixture of clothianidin and deltamethrin) against susceptible and resistant Anopheles gambiae populations from Benin, West Africa. Malar J. 2018;17(1):484.
  19. 19. Fongnikin A, Houeto N, Agbevo A, Odjo A, Syme T, N’Guessan R, et al. Efficacy of Fludora® Fusion (a mixture of deltamethrin and clothianidin) for indoor residual spraying against pyrethroid-resistant malaria vectors: laboratory and experimental hut evaluation. Parasites & Vectors. 2020;13(1):466.
  20. 20. Gavi S, Tapera O, Mberikunashe J, Kanyangarara M. Malaria incidence and mortality in Zimbabwe during the COVID-19 pandemic: analysis of routine surveillance data. Malar J. 2021;20(1):233. pmid:34030711
  21. 21. 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. pmid:31881898
  22. 22. Rodriguez-Barraquer I, Arinaitwe E, Jagannathan P, Kamya MR, Rosenthal PJ, Rek J, et al. Quantification of anti-parasite and anti-disease immunity to malaria as a function of age and exposure. eLife. 2018;7. pmid:30044224
  23. 23. 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. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2017;65(3):453–60.
  24. 24. Ruhamyankaka E, Brunk BP, Dorsey G, Harb OS, Helb DA, Judkins J, et al. ClinEpiDB: an open-access clinical epidemiology database resource encouraging online exploration of complex studies. Gates Open Res. 2019;3:1661. pmid:32047873