Incremental impact on malaria incidence following indoor residual spraying in a highly endemic area with high standard ITN access in Mozambique: results from a cluster‐randomized study

Attaining the goal of reducing the global malaria burden is threatened by recent setbacks in maintaining the effectiveness of vector control interventions partly due to the emergence of pyrethroid resistant vectors. One potential strategy to address these setbacks could be combining indoor residual spraying (IRS) with non-pyrethroids and standard insecticide-treated nets (ITNs). This study aimed to provide evidence on the incremental epidemiological benefit of using third-generation IRS product in a highly endemic area with high ITN ownership. A cluster-randomized, open-label, parallel-arms, superiority trial was conducted in the Mopeia district in Zambezia, Mozambique from 2016 to 2018. The district had received mass distribution of alphacypermethrin ITNs two years before the trial and again mid-way. 86 clusters were defined, stratified and randomized to receive or not receive IRS with pirimiphos-methyl (Actellic®300 CS). Efficacy of adding IRS was assessed through malaria incidence in a cohort of children under five followed prospectively for two years, enhanced passive surveillance at health facilities and by community health workers, and yearly cross-sectional surveys at the peak of the transmission season. A total of 1536 children were enrolled in the cohort. Children in the IRS arm experienced 4,801 cases (incidence rate of 3,532 per 10,000 children-month at risk) versus 5,758 cases in the no-IRS arm (incidence rate of 4,297 per 10,000 children-month at risk), resulting in a crude risk reduction of 18% and an incidence risk ratio of 0.82 (95% CI 0.79–0.86, p-value < 0.001). Facility and community passive surveillance showed a malaria incidence of 278 per 10,000 person-month in the IRS group (43,974 cases over 22 months) versus 358 (95% CI 355–360) per 10,000 person-month at risk in the no-IRS group (58,030 cases over 22 months), resulting in an incidence rate ratio of 0.65 (95% CI 0.60–0.71, p < 0.001). In the 2018 survey, prevalence in children under five in the IRS arm was significantly lower than in the no-IRS arm (OR 0.54, 95% CI, 0.31–0.92, p = 0.0241). In a highly endemic area with high ITN access and emerging pyrethroid resistance, adding IRS with pirimiphos-methyl resulted in significant additional protection for children under five years of age. Trial registration: ClinicalTrials.gov identifier NCT02910934, registered 22 September 2016, https://clinicaltrials.gov/ct2/show/NCT02910934?term=NCT02910934&draw=2&rank=1.


Background
There has been remarkable success in the global fight against malaria since 2000. During the period 2000-2015, coordinated malaria control efforts helped reduce worldwide malaria mortality rates in all ages by 47%, averting an estimated 4.3 million malaria deaths [1]. This progress was particularly impressive in Africa, where infection prevalence was halved and clinical cases reduced by 40%, averting an estimated 663 million cases, during the same time period [2]. Most of this progress (81% of the cases averted) in Africa can be attributed to the successful scale-up of malaria vector control with conventional pyrethroid insecticidetreated nets (ITNs) and indoor residual spraying (IRS) [2].
Despite this overall success, recent trends indicate that maintaining high intervention coverage is challenging and that the number of malaria cases has increased slightly, but consistently, every year since 2016, this is mainly driven by a few high-burden countries [3,4]. This interrupted progress has put the malaria fight at a crossroads [3] and threatens attaining the disease burden reduction targets set forth by the World Health Organization (WHO) in the Global Technical Strategy for Malaria 2016-2030 (GTS) [5]. Further complicating the picture for vector control [5,6] is knowledge that the continued effectiveness of currently available tools is threatened by the spread of insecticide resistance (especially pyrethroid resistance) in key vector populations [7,8]. Indeed, extensive modelling conducted in preparation of the GTS suggests that innovative approaches are needed to get back on track to achieving the proposed goals [9]. These needs include better access to prevention and treatment interventions, better distribution systems, better tools with non-pyrethroid insecticides, and optimized combinations of available tools.
The efficacy of ITNs to reduce malaria incidence, prevalence, and even all-cause child mortality has been well established [10,11], and so this intervention has become the main malaria vector control method worldwide with an estimated 72% of households at risk in sub-Saharan Africa owning at least one ITN in 2017, as compared with 47% in 2010 [4]. Usage has also steadily increased; it is estimated that 50% of the population at risk in sub-Saharan Africa, including 61% of children under 5 years and 61% of pregnant women, slept under an ITN in 2017 [4].
The impact and cost-effectiveness of IRS as a malaria control intervention has also been clearly established by historical and programme documentation [12,13]. One major challenge has been the increased cost of IRS with new insecticides or third generation IRS products (3GIRS). This increased cost of IRS products was associated with a reduction in IRS coverage throughout sub-Saharan Africa. Globally, the proportion of the population at risk protected by IRS was 5% in 2010 but declined to 3% in 2017 as countries identified insecticide resistance and new effective 3GIRS insecticides were more expensive [4]. In sub-Saharan Africa, IRS coverage experienced a marked decline from 10.1% (80 million people protected) in 2010 to 5.4% (51 million people protected) in 2016 before rising again to 6.6% (64 million people protected) in 2017 [4]. This and other intervention coverage gaps, as well as a funding plateau, have been identified as important contributors to the stall in progress seen in 2017 and 2018 [4,14].
The use of 3GIRS, with longer residual activity, in addition to high ITN coverage is one approach that could improve vector control and enhance disease burden reduction in some situations. The current evidence for this potential benefit is mixed. Although modelling suggests additional incremental impact and observational studies suggest added value for IRS in addition to ITNs [15][16][17][18], experimental hut studies [19][20][21], non-randomized [22], and cluster-randomized trials [23][24][25][26][27] show variable impact that is highly dependent on transmission intensity, vector bionomics, ITN coverage, insecticide resistance profiles, implementation strategies, and other factors. A recent metanalysis on the combined used of IRS and long-lasting insecticidal nets (LLINs) concluded that care is needed when using the limited available evidence for policy decisions [28]. Regarding cost-effectiveness, there are important logistical costs associated with IRS, and while, a 2011 systematic review of IRS found it was cost-effective in low income setting [29], only one trial has explicitly evaluated the cost-effectiveness of the combined approach, and did so in the unique context of a low-burden region of Ethiopia [30]. Moreover, the added value and cost-effectiveness of IRS in addition to ITNs in the context of intense transmission areas are critical questions for the new "High burden to high impact" strategy developed by the WHO and RBM [31].
Although implementation is often sub-national, at least 35 countries in Africa already recommend combining ITNs and IRS [4], and the latter is often deployed in areas targeted by mass ITN distribution campaigns. Combining IRS and ITNs can result in different insecticides in the same area [32]. Indeed, the WHO guidelines for vector control [33] suggests that combined deployment can be used as part of an insecticide resistance management strategy, but specifically cautions against introducing a second intervention to compensate for deficiencies in the implementation of the first.
Robust data are needed to guide decisions about prioritizing and combining vector control strategies in the context of different transmission dynamics, changing insecticide resistance patterns, and limited funds [34]. To help address this information gap in a highintensity transmission setting with evidence of emerging pyrethroid resistance, a cluster-randomized trial was conducted, assessing the impact of IRS with a microencapsulated formulation of the organophosphate insecticide pirimiphos-methyl (PM) on malaria transmission, compared to no IRS. Both arms received ITNs in accordance with the national distribution campaigns, which resulted in high ITN access in the IRS and no-IRS arms. This study explores the epidemiological outcomes of malaria incidence and prevalence over two years.

Methods
The overall study concept, setting, and methods of this open-label, controlled, parallel-arm, superiority trial have been previously published [35].

Study setting
The study occurred in rural Mopeia District (population 162,000) [36] in the Zambezia Province of Mozambique during 2016-2018. Zambezia is highly endemic for malaria, with parasite prevalence exceeding 60% and significant direct and indirect costs associated with the disease in some recent assessments [37,38]. The main vectors were Anopheles funestus and Anopheles gambiae sensu lato (s.l.) and data from neighbouring districts showed pyrethroid resistance in Anopheles gambiae s.l. [39].
Access to ITNs was relatively high at baseline in 2016, as Mopeia District received 175,000 pyrethroid ITNs in a mass distribution campaign in 2013 (which represented more than one ITN per habitant) and benefits from routine distribution in antenatal clinics. Mopeia received IRS (with DDT and then pyrethroids) from 2007-2011 and in 2014 [40]. In Mozambique, ITN coverage is sustained through routine distribution at antenatal care clinics.
In 2017, the NMCP conducted a mass ITN distribution campaign with alphacypermethrin-treated ITNs in Mopeia. Ownership among all ages in Mopeia was 54% during the 2017 cross-sectional study and 95% in 2018. Net use in Zambezia among households with at least one ITN was 89% in the 2018 Malaria Indicator Survey. The standard of care at public health facilities (testing of all fevers with a rapid diagnostic test (RDT) or microscopy and provision of treatment with artemisinin-based combination therapy to all positive cases) and from community health workers remained unaltered beyond study efforts to prevent stock outs of malaria commodities. There were 30 community health workers providing passive testing, treatment and reporting in Mopeia throughout the study period.

Intervention
Given expected impact at community level, IRS with PM was implemented only in the IRS-assigned clusters by the President's Malaria Initiative Africa Indoor Residual Spraying (PMI AIRS) project from October-November in both 2016 and 2017 (Fig. 1). Spraying was conducted according to PMI AIRS standard operating procedures, including community and household consent.

Study design
The study employed a two-arm, cluster randomized, controlled study design. A household and population enumeration was conducted from June-July 2016, which identified 139,286 total residents (26,320 under five years old) living in 21,328 households distributed across 194 villages. Cluster limits were delineated using expanded village borders through Voronoi polygons, with villages not reaching the minimum population for inclusion combined with the nearest neighbouring village to form a single cluster [27].
The primary research question was: In an area with high malaria endemicity and high ITN access, what is the incremental benefit of IRS with PM on reducing malaria transmission in a cohort of children under five years of age? The primary outcome was malaria infection incidence in an active cohort of children under five years of age at community level. Secondary outcomes included: (1) passively reported confirmed case incidence in all ages through the national health system, including health facilities and community health workers and, (2) malaria prevalence in all ages from annual cross-sectional surveys near the peak of the transmission season (April-May).

Randomization and masking
The 168 clusters were stratified into three groups according to the number of households (< 69 = small; 69-125 = medium; > 125 = large), and randomized 1:1 into one of the two arms, IRS and no-IRS, by drawing lots during a public community-engagement ceremony.

Entomological surveillance
The standard PMI AIRS Mozambique vector surveillance methods and study-specific sampling strategies have been previously described [35,41]. In short, vector densities were monitored monthly in a subset of ten sentinel study villages (selected based on preliminary mosquito density surveys, ease of access and safety) five IRS and five no-IRS villages. In each sentinel village, overnight CDC light trap collections were conducted at eight houses for three consecutive nights every month. At one additional house per village, paired indoor-outdoor human landing collections were conducted overnight on the same three consecutive nights. Subsequent molecular analyses of the collected specimens (species confirmation, Plasmodium spp. infection rates, and appropriate pyrethroid resistance marker frequencies) have been reported in a separate publication [42].
Standard WHO cone wall bioassay tests were performed at a subset of five randomly selected households in each of three villages in Mopeia to assess initial spray quality and estimate the residual efficacy of PM [43]. Larval collections and subsequent insecticide resistance profiling of An. gambiae s.l. (2017 and 2018) and An. funestus s.l. (2018) using the WHO tube test bioassay also followed standard PMI AIRS Mozambique methods [41,43].

Primary outcome measures
Active cohort 86 total clusters (43 IRS, 43 no-IRS) were selected for participation in the active cohort component of the study. Eligible households were selected from the core zones of each cluster using a fried-egg design [44] with a 1-km buffer zone at the margins of each cluster, effectively leaving a buffer of at least 2 km between spraydiscordant core zones [35]. No buffers were included between clusters that had been randomized to the same study arm. Malaria infection incidence at community level was determined by enrolling a cohort of children under five years of age under parental informed consent (18 children per cluster, 774 per study arm). These children were visited monthly by a trained field worker that administered a short questionnaire to the caregiver and performed an HRP2-based rapid diagnostic test (RDT). Every child with a positive RDT received treatment with artemether-lumefantrine (AL) according to Mozambique National Malaria Control Programme guidelines at baseline and in every subsequent visit. Person-time at risk was reduced by ten days after each treatment to account for the prophylactic effect of lumefantrine [45].

Passive case detection
The incidence of confirmed malaria cases (defined as fever, either reported or measured plus a positive RDT) that sought care in the public health system in Mopeia was measured using an enhanced passive surveillance approach: a study worker was placed in each of the 13 health facilities in the district to assure the quality of malaria case recording and to register the village origin of every case by village study-code.

Cross-sectional surveys
A cross-sectional survey was conducted in April-May in 2017 and again in 2018 to assess malaria prevalence in all-ages and to gather behavioural information as well data on costing, and health care expenditure.

Statistical considerations
For the active cohort, 42 clusters of 12 children per arm had 80% power at a 5% significance level to detect a reduction in baseline incidence of 30% (from estimated 700/1,000 [46] children-years to 490/1,000 childrenyears), using a robust K of 0.5. The number of clusters per arm was 43 and the number of children per cluster was 18 at enrolment to account for potential sample loss. Power and sample size calculations were conducted using the Hayes and Bennett formula [47]. The sampling strategy for each cross-sectional survey (770 individuals, half under five years of age) aimed for 5% precision to measure an estimated prevalence of 50% in a population of 128,000.
Primary analysis was done on intention-to-treat, assuming that all individuals living in an IRS cluster received IRS in their household. The effect of IRS was estimated using negative binomial regression models with the generalized estimating equations (GEE) approach. This effect was adjusted for the variables identified as potential confounders in univariate models; the interaction term between IRS and ITNs was included in the multivariate analysis. Sensitivity analyses and additional per protocol analysis adjustments were done considering ITN ownership and usage, household socioeconomic status, and cluster size (as defined by number of households). The analysis was performed using Stata Statistical Software (StataCorp 2017).

Ethical reviews and registration
All procedures were reviewed and approved by PATH's Research Ethics Committee, CISM's IRB, and the National Ethics Committee of Mozambique as well as the PMI Operational Research Committee. This study was reviewed by the Centers for Disease Control and Prevention (CDC) and determined to be human subjects research with non-engagement by CDC staff. The trial was registered at clinicaltrials.gov with the identifier NCT02910934.

Total population and study flow
The study enumeration and enrolment process are depicted in Fig. 2.

IRS quality, acceptance and ITN distribution
The IRS campaigns were well accepted in Mopeia, with 16,500 structures (83%) sprayed of the 19,992 target in 2016 and 16,936 structures (85%) of the 19,950 target sprayed in 2017. These targets were based on PMI AIRSled structure enumeration which was conducted prespray each year. Standard WHO cone bioassays using susceptible An. gambiae sensu stricto (s.s.) in houses from Mopeia district indicated that both IRS campaigns were of high quality, with all houses tested in both years showing 100% mortality within 48 h of spraying. Additionally, PM was efficacious for at least three months in 2017 [41]. In 2018, results again showed residual efficacy for a minimum of three months on all wall surface types tested, though on mud walls efficacy was of longer duration and lasted for at least four months [43].
In June-July 2017, all villages in Mopeia received 120,765 ITNs in the context of the mass distribution campaign. There were no reported stock-outs of RDTs or anti-malarials reported at health facilities in Mopeia during the study period. Following the 2017 campaign, the four-month time point measurement showed 82% mosquito mortality on mud walls of Cero village and a five-month measurement in neighbouring Mocuba and Morrumbala districts showed 95% mortality.

Active cohort detection
Baseline characteristics were calculated using the active cohort first measurement to ensure comparability between clusters as the passive case recording did not delineate village of origin prior to the study period. A total of 1,536 children under five years of age (765 the no-IRS arm and 771 in the IRS arm), were enrolled in the active cohort from the 86 clusters (43 ITNs-only and 43 ITNs + IRS). The distribution of cluster size was equal in both groups, with 14 small, 14 medium and 15 large clusters per arm.
The baseline characteristics of the cohort are shown in Table 1. There were no major differences in terms of distance to the nearest health facility from the cluster´s centroid, ITN ownership, basic socioeconomic characteristics, age, or gender of the children enrolled. More than 60% of the children had a positive RDT at enrolment.
The comparison of factors potentially associated with a positive RDT between both groups at baseline is presented in the Supplementary Materials. Specifically, there were slight associations between malaria test positivity and living in a medium or large cluster, having a household sibling also testing positive, younger age, longer distances to the nearest health facility, and a history of fever in the last 48 h (Additional file 1: Table S1).
ITN ownership data collected at baseline and after the mass distribution campaign showed consistency across both arms and a large increase from 61-63% in January 2017 to 90% ownership of at least one ITN by the end of the trial (Additional file 1: Fig. S1).
The children in the IRS arm experienced a significantly lower malaria infection incidence throughout the study. There were 4,801 cases in the IRS arm (incidence rate of 3,532 per 10,000 children-month at risk) versus 5,758 cases in the no-IRS arm (incidence rate of 4,297 per 10,000 children-month at risk). The crude risk reduction was 18% and the incidence risk ratio (IRR) was 0.82 (95% CI: 0.79, 0.86, p-value < 0.001) ( Table 2 and Fig. 3).
Using these data, the IRS campaign in Mopeia averted between 15,697 and 21,651 malaria infections in the 12,670 children under five years of age living in the IRS clusters from January 2017 to October 2018. A sensitivity analysis was conducted, adjusting to account for residual  Fig. 2).
The coefficient of variation (k) between clusters for RDT positive test results was calculated to be 0.336 using GEE. Given the crude incidence reduction of 18%, from 4,297 infections per 10,000 children-month (5.1 cases per child-year) in the no-IRS arm to 3,532 per 10,000 children-month (4.2 cases per child-year) in the IRS arm, the study had 74% power for its primary outcome with 43 clusters of 18 children per arm. Univariate GEE with negative binomial models were used to explore which potentially confounding factors should be included in the multivariate model. These results are shown in Table 3.
Variables identified as having a significant influence on the IRR in the univariate analysis were included in a multivariate model using GEE. Table 4 shows the corrected IRR associated with IRS alone (i.e. no ITN owned), ITN use the night before, combined IRS + ITN use, having a sibling who tested positive, cluster size or distance to the nearest HF. The combined effect of sleeping under an ITN the night before in a cluster that received the IRS intervention was significantly greater than the effect of either intervention used alone: the adjusted IRR for the interaction term was 0.62 (95% CI 0.57 -0.67; p < 0.001) corresponding to an incidence reduction of 38% (95% CI 33%-43%). The incidence reduction associated with IRS alone was 19% (95% CI 13-26%) and 23% (95% CI 18-28%) with ITN use alone (Table 4). There was a 21% risk increase in children with at least one other sibling in the cohort that tested positive. There was small but statistically significant reduction in the IRR in larger clusters and a higher risk of malaria in clusters with longer distances to health facilities (Table 4).
Sensitivity analyses were performed including only data after the ITN distribution campaign or adjusting the reference category for the IRR and no major changes in these results were noted (Additional file 1: Table S2).

Passive case detection at HF
There was a total of 188 distinct villages coded during the district-wide pre-study enumeration, 81 that received IRS and 107 that did not, as per randomization which excluded a few villages that were not accessible for logistical or instability reasons. The total enumerated population was 138,685, of which 18.8% (26,097) were under the age of 5 years. Slightly less than half the total enumerated population lived in IRS villages. (68,725 = 49.6%).
There were 380,727 total visits to health facilities and to community health workers in Mopeia recorded during the study period; of these, 365,741 (96%) had a village code corresponding to a spray status with the rest corresponding to patients from outside the district boundaries, patients unwilling to disclose their home address, From no-IRS villages, a total of 58,030 RDT-confirmed malaria cases were recorded over 22 months, resulting in a crude all-ages case incidence rate of 361 per 10,000 person-months at risk. Almost half (48.7%) of these confirmed cases were in children under five years of age, an age-specific case incidence rate of 916 per 10,000 childmonths at risk in this population. There were significantly fewer confirmed cases of malaria recorded from IRS villages: 43,973 total cases, resulting in a crude all-ages case incidence rate of 278 per 10,000 person-months at risk. Of the cases from IRS villages, 45.5% were in children under five years of age, an age-specific case incidence rate of 687 per 10,000 child-months at risk.
The malaria incidence at health facility in the overall population was 358 (95% CI: 355-360) per 10,000 person-month at risk in the no-IRS group (58,030 cases over 22 months) and 278 per 10,000 person-month in the IRS group (43,974 cases over 22 months), resulting in an incidence rate ratio of 0.65 (95% CI 0.60-0.71, p < 0.001). The number of averted cases was estimated to be between 15,697 and 21,651. Monthly case incidence in both arms for the overall population and children under five years of age are shown in Table 5 and Fig. 4.
The crude incidence was adjusted using a negative binomial regression model with variables identified via univariate regression. These results confirmed the lower IRR in larger clusters (IRR: 0.98 per every 100 population increase 95% CI 0.98-0.99 p < 0.0001) and also revealed an increased risk of malaria detection in clusters with shorter linear distance to a health facility; people living closer to a health facility received an RDT with higher frequency (IRR: 0.68 per every 5-km increase in distance 95% CI 0.65-0.71, p < 0.0001). Data are presented in Supplementary Tables 3 and 4.

Cross sectionals
A total of 822 participants were surveyed in 2017 and 805 in 2018. Both samples were balanced in terms of age, gender, ITN ownership, and other relevant factors (Table 6). In the 2017 survey, conducted before the mass distribution of ITNs, there was no significant difference in prevalence at the peak of the transmission season between both study arms, even when correcting by age of the participant or ITN ownership (Tables 7 and 8). In the 2018 survey, conducted ten months after ITN distribution, prevalence in children under five years of age in the IRS arm was significantly lower than in the no-IRS arm (OR 0.54, 95% CI 0.31-0.92, p = 0.0241). The incremental protective effect was particularly marked among ITN owners compared to those with no ITNs (Table 8).

Entomological characterization
A full analysis of the entomological impact of the IRS campaigns will be presented in a complementary manuscript [42]. In terms of characterizing the underlying vector bionomics at the sentinel sites, more than 90% of all anophelines collected (23, [41,42]. The WHO tube test results from 2015 showed that pyrethroid resistance was evident in An. gambiae s.l. populations from the nearby districts of Mocuba (52% mortality against deltamethrin/40% against lambda cyhalothrin) and Morrumbala (34% mortality against deltamethrin/33% against lambda cyhalothrin) [39], although data from Mopeia district in 2017 showed that An. gambiae s.l. was 100% susceptible to both alphacypermethrin and to PM. Anopheles funestus s.l. from Mopeia were tested in 2018 and were 100% susceptible to PM and DDT, but showed signs of emerging resistance to alphacypermethrin (85% mortality), deltamethrin (88% mortality), and bendiocarb (89% mortality) [42,43].

Discussion
The IRS campaigns of 2016 and 2017 made positive contributions to malaria control in this high transmission district of Mozambique, as evident by: (1) reduced infection incidence in IRS clusters relative to no-IRS clusters, even in the presence of high ITN ownership; (2) reduced confirmed clinical case incidence at public health clinics; and among those detected by community health workers, and (3) reduced odds of malaria infection in the population under five years of age during the 2018 prevalence survey.
Malaria policy makers and implementers face difficult decisions regarding the best available tools and their optimal deployment. Prior cluster-randomized trials have provided differing results, suggesting that the added value of the combination of ITNs and IRS is variable and likely. Fig. 3 Cohort incidence by spray status (a); cohort cumulative incidence by spray status (b); and spray IRR (with 95% confidence interval) at cohort level (c). IRS campaigns highlighted in blue and mass ITN distribution highlighted in grey, ACT treatment correction of time at risk: 10 days dependent on local factors like transmission intensity, vector bionomics, insecticide resistance profiles, and implementation strategy [23][24][25][26][27].
This study generated robust evidence to help support those making policy and implementation decisions about the use of IRS with a non-pyrethroid insecticide in communities with high rates of malaria transmission, high ITN ownership of standard pyrethroid-only ITNs, and evidence of emerging pyrethroid resistance in the local vector populations.
This study found significant added malaria protection by adding IRS with PM to a policy of universal coverage with a pyrethroid-only ITN in Mopeia. This was quantified as 18% protective efficacy when considering new P. falciparum infections detected in the incidence cohort, and around 28% protective efficacy when considering confirmed cases reporting to the public health system for treatment. These suggest that when resources are available, combining these two interventions will reduce malaria incidence. Some of the reasons contributing to this incremental impact include the indoor insecticide-mosaic created by combining pyrethroid ITNs and organophosphate IRS and the benefit obtained at the household level from at least one insecticide present independently of compliance with ITN use. The adjusted analysis performed with the active cohort data confirmed that the interaction of IRS and ITNs leads to the greater incidence reduction, namely 38%. While this study provides valuable evidence on the combination of interventions, it also highlights the need to generate evidence on the value of IRS in combination with ITNs with piperonyl butoxide and non-pyrethroid ITNs to inform programmatic decision-making.
The reduced odds ratio of malaria infection observed in the population under five years of age during the 2018 prevalence survey is particularly interesting as it aligns with the protective effect of IRS also observed in the   (2) even as prevalence in the overfive population fell from 40 to 25% during the same time.
Additionally, the combination of IRS and ITNs appeared to significantly reduce the odds of malaria infection in the under-five population by almost 50% during the five months after the second spray campaign and following the mass distribution campaign. These interesting trends highlight how complex the relationship between malaria infection incidence, malaria clinical case incidence, and malaria infection prevalence can be, particularly in very highly endemic areas with year-round transmission. These results are in contrast from those of previous cluster-randomized trials conducted in lower transmission settings in which no benefit with a combined IRS and ITN approach showed no added benefit [23,24,27], but in concordance with the positive results seen in higher transmission settings [25,32]. This raises the question of whether the incremental impact maybe dependent on the transmission level.
This study employed a robust cluster-randomized design, with a multiplicity of outcome measures, a large sample size, and close community engagement, but there are some important limitations. Children in the active cohort were subject to screen and treatment every visit, resulting in early detection of infections followed by prompt treatment, which also provided temporary prophylaxis. Given that malaria infections were more common in children in the no-IRS arm, they received proportionally more treatments (and potential prophylactic benefit) which may have reduced the difference between arms, resulting in an underestimation of the true added benefit of the combined approach as adjustments done at analysis cannot fully correct for this prophylaxis. Additionally, the net reduction in malaria incidence was lower than originally expected for sample size calculations; this was, however, partly compensated by the coefficient of variation, which, once retrospectively calculated from empirical data, was lower than assumed. Another potential source of bias could be falsepositive RDT results. The RDTs used are based on the HRP2 antigen, which can persist for several weeks after treatment potentially inflating estimates of incidence in the active cohort and adding uncertainty around the number of true infections [48]. This effect would, however, occur equally in both groups. A sensitivity analysis adjusting incidence rates by censoring positive RDT results from a second consecutive household visit showed no major difference with the main findings presented here ( Supplementary Fig. 2). Both study arms also benefitted from the study team efforts to avoid stock-outs which could have contributed to lower the incidence at cohort and health facility level in both arms. Despite these potential biases, it is reassuring to note the consistency among all outcome measures with active cohort, passive surveillance and cross sectionals.
Malaria remains a challenge that will require multiple preventive as well as therapeutic intervention strategies, and many of these tools will need to be used in combination to maximize impact and reach elimination goals. Understanding when and where to combine vector control strategies requires locally relevant data to ensure that resources are invested wisely, particularly in the context of the "High burden for high impact" strategy. This study demonstrated added value for IRS with a non-pyrethroid active ingredient in the context of high coverage with standard (pyrethroid-only) ITNs, as well as good access to malaria case management commodities. This supports consideration of co-investment strategy (IRS and ITNs) in areas such as Zambezia, where transmission is high and the local primary vector species, An. funestus s.s., shows moderate levels of pyrethroid resistance. The cost-effectiveness of this combined approach has been analysed in the context of this trial resulting in a separate manuscript.

Conclusion
In 2017, Mozambique had 5% of the global share of malaria cases [4]. The results of this trial suggest consideration for the combined deployment of non-pyrethroid IRS with ITNs in areas of high transmission and emerging pyrethroid resistance. This strategy could prove especially valuable in the context of an overall increase in malaria burden and strategy put in place in an attempt to get back on track to achieving the 2030 goals as outlined in the WHO Global Technical Strategy [31]. Monthly population incidence at health facilities by spray status (a); cumulative population incidence at health facilities by spray status (b); and monthly incidence rate ratio (c). IRS campaigns highlighted in blue and mass ITN distribution highlighted in grey