<!DOCTYPE article
PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD with MathML3 v1.2 20190208//EN" "JATS-archivearticle1-mathml3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article"><?properties open_access?><front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id><journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id><journal-title-group><journal-title>Scientific Reports</journal-title></journal-title-group><issn pub-type="epub">2045-2322</issn><publisher><publisher-name>Nature Publishing Group UK</publisher-name><publisher-loc>London</publisher-loc></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">32161302</article-id><article-id pub-id-type="pmc">7066154</article-id><article-id pub-id-type="publisher-id">61350</article-id><article-id pub-id-type="doi">10.1038/s41598-020-61350-2</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Impact of indoor residual spraying with pirimiphos-methyl (Actellic 300CS) on entomological indicators of transmission and malaria case burden in Migori County, western Kenya</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-5038-604X</contrib-id><name><surname>Abong&#x02019;o</surname><given-names>Bernard</given-names></name><address><email>abongoben@gmail.com</email></address><xref ref-type="aff" rid="Aff1">1</xref><xref ref-type="aff" rid="Aff2">2</xref><xref ref-type="aff" rid="Aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Gimnig</surname><given-names>John E.</given-names></name><xref ref-type="aff" rid="Aff4">4</xref></contrib><contrib contrib-type="author"><name><surname>Torr</surname><given-names>Stephen J.</given-names></name><xref ref-type="aff" rid="Aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Longman</surname><given-names>Bradley</given-names></name><xref ref-type="aff" rid="Aff1">1</xref></contrib><contrib contrib-type="author"><name><surname>Omoke</surname><given-names>Diana</given-names></name><xref ref-type="aff" rid="Aff1">1</xref></contrib><contrib contrib-type="author"><name><surname>Muchoki</surname><given-names>Margaret</given-names></name><xref ref-type="aff" rid="Aff1">1</xref></contrib><contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0003-3663-5617</contrib-id><name><surname>ter Kuile</surname><given-names>Feiko</given-names></name><xref ref-type="aff" rid="Aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Ochomo</surname><given-names>Eric</given-names></name><xref ref-type="aff" rid="Aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Munga</surname><given-names>Stephen</given-names></name><xref ref-type="aff" rid="Aff3">3</xref></contrib><contrib contrib-type="author"><name><surname>Samuels</surname><given-names>Aaron M.</given-names></name><xref ref-type="aff" rid="Aff4">4</xref></contrib><contrib contrib-type="author"><name><surname>Njagi</surname><given-names>Kiambo</given-names></name><xref ref-type="aff" rid="Aff5">5</xref></contrib><contrib contrib-type="author"><name><surname>Maas</surname><given-names>James</given-names></name><xref ref-type="aff" rid="Aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Perry</surname><given-names>Robert T.</given-names></name><xref ref-type="aff" rid="Aff6">6</xref></contrib><contrib contrib-type="author"><name><surname>Fornadel</surname><given-names>Christen</given-names></name><xref ref-type="aff" rid="Aff7">7</xref></contrib><contrib contrib-type="author"><name><surname>Donnelly</surname><given-names>Martin J.</given-names></name><xref ref-type="aff" rid="Aff2">2</xref></contrib><contrib contrib-type="author"><name><surname>Oxborough</surname><given-names>Richard M.</given-names></name><xref ref-type="aff" rid="Aff8">8</xref></contrib><aff id="Aff1"><label>1</label>Abt Associates, PMI VectorLink Project, White House, Milimani, Ojijo Oteko Road, P.O. Box 895-40123, Kisumu, Kenya </aff><aff id="Aff2"><label>2</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1936 9764</institution-id><institution-id institution-id-type="GRID">grid.48004.38</institution-id><institution>Liverpool School of Tropical Medicine, Pembroke Place, </institution></institution-wrap>Liverpool, L3 5QA UK </aff><aff id="Aff3"><label>3</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0155 5938</institution-id><institution-id institution-id-type="GRID">grid.33058.3d</institution-id><institution>Centre for Global Health Research, Kenya Medical Research Institute, P.O. Box 1578, </institution></institution-wrap>Kisumu, Kenya </aff><aff id="Aff4"><label>4</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 0540 3132</institution-id><institution-id institution-id-type="GRID">grid.467642.5</institution-id><institution>Division of Parasitic Diseases and Malaria, Center for Global Health, Centers for Disease Control and Prevention, </institution></institution-wrap>Atlanta, GA 30333 USA </aff><aff id="Aff5"><label>5</label>Kenya National Malaria Control Programme (NMCP), Ministry of Health, PO Box 19982, Kenyatta National Hospital, Nairobi, 00202 Kenya </aff><aff id="Aff6"><label>6</label>The United States Presidents Malaria Initiative (PMI), US Embassy Nairobi, United Nations Avenue, Nairobi, Kenya </aff><aff id="Aff7"><label>7</label><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 1955 0561</institution-id><institution-id institution-id-type="GRID">grid.420285.9</institution-id><institution>The United States Presidents Malaria Initiative (PMI), US Agency for International Development, </institution></institution-wrap>Washington, DC USA </aff><aff id="Aff8"><label>8</label>PMI VectorLink Project, Abt Associates 6130 Executive Blv, Rockville, MD 20852 USA </aff></contrib-group><pub-date pub-type="epub"><day>11</day><month>3</month><year>2020</year></pub-date><pub-date pub-type="pmc-release"><day>11</day><month>3</month><year>2020</year></pub-date><pub-date pub-type="collection"><year>2020</year></pub-date><volume>10</volume><elocation-id>4518</elocation-id><history><date date-type="received"><day>8</day><month>9</month><year>2019</year></date><date date-type="accepted"><day>25</day><month>2</month><year>2020</year></date></history><permissions><copyright-statement>&#x000a9; The Author(s) 2020</copyright-statement><license license-type="OpenAccess"><license-p><bold>Open Access</bold> This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article&#x02019;s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article&#x02019;s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>.</license-p></license></permissions><abstract id="Abs1"><p id="Par1">Indoor residual spraying (IRS) of insecticides is a major vector control strategy for malaria prevention. We evaluated the impact of a single round of IRS with the organophosphate, pirimiphos-methyl (Actellic 300CS), on entomological and parasitological parameters of malaria in Migori County, western Kenya in 2017, in an area where primary vectors are resistant to pyrethroids but susceptible to the IRS compound. Entomological monitoring was conducted by indoor CDC light trap, pyrethrum spray catches (PSC) and human landing collection (HLC) before and after IRS. The residual effect of the insecticide was assessed monthly by exposing susceptible <italic>An. gambiae</italic> s.s. Kisumu strain to sprayed surfaces in cone assays and measuring mortality at 24&#x02009;hours. Malaria case burden data were extracted from laboratory records of four health facilities within the sprayed area and two adjacent unsprayed areas. IRS was associated with reductions in <italic>An. funestus</italic> numbers in the intervention areas compared to non-intervention areas by 88% with light traps (risk ratio [RR] 0.12, 95% CI 0.07&#x02013;0.21, p&#x02009;&#x0003c;&#x02009;0.001) and 93% with PSC collections (RR&#x02009;=&#x02009;0.07, 0.03&#x02013;0.17, p&#x02009;&#x0003c;&#x02009;0.001). The corresponding reductions in the numbers of <italic>An. arabiensis</italic> collected by PSC were 69% in the intervention compared to the non-intervention areas (RR&#x02009;=&#x02009;0.31, 0.14&#x02013;0.68, p&#x02009;=&#x02009;0.006), but there was no significant difference with light traps (RR&#x02009;=&#x02009;0.45, 0.21&#x02013;0.96, p&#x02009;=&#x02009;0.05). Before IRS, <italic>An. funestus</italic> accounted for over 80% of <italic>Anopheles</italic> mosquitoes collected by light trap and PSC in all sites. After IRS, <italic>An. arabiensis</italic> accounted for 86% of <italic>Anopheles</italic> collected by PSC and 66% by CDC light trap in the sprayed sites while the proportion in non-intervention sites remained unchanged. No sporozoite infections were detected in intervention areas after IRS and biting rates by <italic>An. funestus</italic> were reduced to near zero. <italic>Anopheles funestus</italic> and <italic>An. arabiensis</italic> were fully susceptible to pirimiphos-methyl and resistant to pyrethroids. The residual effect of Actellic 300CS lasted ten months on mud and concrete walls. Malaria case counts among febrile patients within IRS areas was lower post- compared to pre-IRS by 44%, 65% and 47% in Rongo, Uriri and Nyatike health facilities respectively. A single application of IRS with Actellic 300CS in Migori County provided ten months protection and resulted in the near elimination of the primary malaria vector <italic>An. funestus</italic> and a corresponding reduction of malaria case count among out-patients. The impact was less on <italic>An. arabiensis</italic>, most likely due to their exophilic nature.</p></abstract><kwd-group kwd-group-type="npg-subject"><title>Subject terms</title><kwd>Population dynamics</kwd><kwd>Entomology</kwd><kwd>Epidemiology</kwd></kwd-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><funding-group><award-group><funding-source><institution>Funder: USAID President&#x00026;apos;s Malaria Initiative (PMI) Grant Reference Number: Task Order #1, AID-OAA-TO-17-00027</institution></funding-source></award-group></funding-group><custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name><meta-value>&#x000a9; The Author(s) 2020</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="Sec1" sec-type="introduction"><title>Introduction</title><p id="Par2">Over the last two decades, malaria control has been scaled up throughout sub-Saharan Africa with an emphasis on the distribution of long-lasting insecticidal nets (LLINs), targeted application of indoor residual spraying (IRS), and improved diagnostics and case management. As a result, the burden of malaria has declined substantially with a 40% reduction in incidence and a 50% reduction in prevalence between 2000 and 2015. While LLINs contributed an estimated 68% of the decline in malaria prevalence, IRS was responsible for 13%<sup><xref ref-type="bibr" rid="CR1">1</xref></sup>.</p><p id="Par3">The efficacy of insecticide-treated nets was demonstrated in a series of cluster randomized, controlled trials<sup><xref ref-type="bibr" rid="CR2">2</xref></sup>. Formal randomized controlled trials of IRS have also demonstrated the efficacy of IRS<sup><xref ref-type="bibr" rid="CR3">3</xref>,<xref ref-type="bibr" rid="CR4">4</xref></sup>. Furthermore, there is a long history of programmatic implementation of IRS in many settings of the world which resulted in reduced malaria burden and even elimination in some settings<sup><xref ref-type="bibr" rid="CR2">2</xref></sup>. The use of both LLINs and IRS for malaria control has a direct impact on mosquito bionomics. LLINs and IRS have multiple effects on mosquito populations which may result in reduced malaria transmission including: reduced indoor <italic>Anopheles</italic> densities<sup><xref ref-type="bibr" rid="CR5">5</xref>&#x02013;<xref ref-type="bibr" rid="CR7">7</xref></sup>, shifts in vector species composition<sup><xref ref-type="bibr" rid="CR8">8</xref></sup>, changes in the time and location of mosquito biting<sup><xref ref-type="bibr" rid="CR9">9</xref>&#x02013;<xref ref-type="bibr" rid="CR11">11</xref></sup>, and changes in host selection<sup><xref ref-type="bibr" rid="CR12">12</xref></sup>, and increases in early exophily<sup><xref ref-type="bibr" rid="CR13">13</xref></sup>.</p><p id="Par4">In western Kenya, vector control included universal coverage of LLINs through periodic mass campaigns and routine distribution to high-risk groups as well as IRS in specifically targeted areas. The first mass LLIN distribution occurred in 2006 and targeted children &#x0003c;5 years of age. Additional distributions aiming for universal coverage occurred in 2011 and 2014 leading to 54% of households in the lake endemic zone having one LLIN for every two residents<sup><xref ref-type="bibr" rid="CR14">14</xref></sup>. The region also bears the highest malaria burden nationally<sup><xref ref-type="bibr" rid="CR14">14</xref>&#x02013;<xref ref-type="bibr" rid="CR16">16</xref></sup>. Implementation challenges facing LLINs include incomplete coverage<sup><xref ref-type="bibr" rid="CR14">14</xref>,<xref ref-type="bibr" rid="CR17">17</xref>&#x02013;<xref ref-type="bibr" rid="CR19">19</xref></sup>, widespread pyrethroid resistance<sup><xref ref-type="bibr" rid="CR20">20</xref>&#x02013;<xref ref-type="bibr" rid="CR23">23</xref></sup> and possible changes in&#x000a0;vector behaviour<sup><xref ref-type="bibr" rid="CR9">9</xref></sup>.</p><p id="Par5">IRS in western Kenya was based exclusively on pyrethroids until 2012<sup><xref ref-type="bibr" rid="CR24">24</xref>&#x02013;<xref ref-type="bibr" rid="CR26">26</xref></sup>. However, spraying was interrupted between 2013 and 2017 due to widespread pyrethroid resistance in local malaria vector populations and the lack of a registered, non-pyrethroid insecticide in the country. In response to widespread pyrethroid resistance, the Kenyan National Malaria Control Programme (NMCP) developed an insecticide resistance management strategy involving the rotation of different non-pyrethroid classes of insecticides used in IRS every two years in endemic and epidemic-prone areas where 80% or more households own one or more LLIN<sup><xref ref-type="bibr" rid="CR27">27</xref></sup>. This is in accordance with global insecticide resistance management strategy aimed at delaying the rise and spread of insecticide resistance to new classes of insecticide while preserving pyrethroids for use in bednets<sup><xref ref-type="bibr" rid="CR28">28</xref></sup>. In 2017, IRS was re-introduced using a microencapsulated formulation of pirimiphos-methyl (Actellic 300CS&#x000a0;(capsule suspension)). The insecticide has been reported to be effective against pyrethroid-resistant <italic>Anopheles</italic> mosquitoes<sup><xref ref-type="bibr" rid="CR29">29</xref>&#x02013;<xref ref-type="bibr" rid="CR31">31</xref></sup> and has a relatively long residual effect on sprayed wall surfaces, of up to twelve months<sup><xref ref-type="bibr" rid="CR29">29</xref>,<xref ref-type="bibr" rid="CR30">30</xref>,<xref ref-type="bibr" rid="CR32">32</xref></sup>.</p><p id="Par6">Between 2008 and 2012, IRS campaigns in western Kenya were conducted using pyrethroid based insecticides. These campaigns resulted in significantly lower parasitaemia and cases of clinical malaria in IRS districts compared to non-IRS districts<sup><xref ref-type="bibr" rid="CR26">26</xref></sup>. The protective efficacy of ITN plus IRS compared with ITN only in the same area was measured at 62%, whereas, 60% fewer anophelines were observed in the sprayed districts compared to unsprayed districts post-IRS<sup><xref ref-type="bibr" rid="CR25">25</xref></sup>. <italic>Anopheles funestus</italic> accounted for 1.7% of total anophelines collected in the sprayed areas before IRS and 2.7% after IRS<sup><xref ref-type="bibr" rid="CR25">25</xref></sup>, an indication that the species was not affected by spraying. Following reports of increasing insecticide resistance, the malaria vector control policy changed to the use of non-pyrethroid insecticides for IRS. IRS was reintroduced in 2017 with an organophosphate, Actellic 300CS. We sought to determine efficacy of the Actellic&#x000a0;CS based IRS campaign on entomological and epidemiological indices against a background of moderate to high coverage of pyrethroid LLINs<sup><xref ref-type="bibr" rid="CR14">14</xref></sup> and extensive pyrethroid resistance.</p></sec><sec id="Sec2"><title>Methods</title><sec id="Sec3"><title>Study sites</title><p id="Par7">Entomological monitoring was conducted in 12 villages in Migori (&#x02212;1.0667S; 34.4667E) and Homa Bay (&#x02212;0.5396S; 34.4565E) counties from July 2016 to February 2018. Six IRS intervention sub-counties were in Migori County, and six control sub-counties were in neighbouring Homa Bay County (n&#x02009;=&#x02009;4) and unsprayed areas of Migori County (n&#x02009;=&#x02009;2) (Fig.&#x000a0;<xref rid="Fig1" ref-type="fig">1</xref>). The residents in the study area are mainly of the Luo ethnic group and are subsistence farmers with a few growing cash crops such as sugar cane and tobacco. Residents mostly live in small houses, clustered into family social units called compounds. The region has bimodal peaks of rainfall with the long rains between April and June and short rains in October and November. The Lake Victoria region of western Kenya is malaria endemic; the most recent Malaria Indicator Survey in 2015 documented a malaria prevalence of 27% by microscopy. The last mass net campaign in the region was conducted in April 2017 with DawaPlus 2.0, a deltamethrin long-lasting (coated) insecticidal net distributed in this area. Though 87% of households own at least one LLIN and 60% own more than one LLIN, only 54% of households have an adequate number of nets, defined as one LLIN for every two residents. Among residents of western Kenya, 68.4% reported to have slept under any net the night before the survey while 66.9% had slept under LLINs<sup><xref ref-type="bibr" rid="CR14">14</xref></sup>. <italic>Anopheles funestus</italic>, <italic>An. arabiensis</italic> and <italic>An. gambiae</italic> s.s. are the main malaria vectors in the region<sup><xref ref-type="bibr" rid="CR33">33</xref>,<xref ref-type="bibr" rid="CR34">34</xref></sup>. Insecticide resistance among major malaria vectors in Kenya has been reported for pyrethroids, carbamates, and organochlorines<sup><xref ref-type="bibr" rid="CR35">35</xref></sup>. The reported resistance mechanisms include both target site mutations and increased activity of enzymes involved in metabolic detoxification<sup><xref ref-type="bibr" rid="CR21">21</xref>&#x02013;<xref ref-type="bibr" rid="CR23">23</xref></sup>.<fig id="Fig1"><label>Figure 1</label><caption><p>Map of Kenya, showing study sites in western Kenya with the names of sub-counties. Yellow shading represents non-intervention sites and red dots represent sampled houses. The green shading is the intervention site with the blue dots representing sampled houses (The map was created on ArcGIS 10.2.2).</p></caption><graphic xlink:href="41598_2020_61350_Fig1_HTML" id="d29e633"/></fig></p><sec id="Sec4"><title>IRS campaign</title><p id="Par8">IRS was conducted in February-March 2017. A total of 212,029 houses in Migori County were sprayed, representing coverage of 97.7% of houses sprayed against houses found. The campaign covered a population of 906,388 people, including 16,932 pregnant women and 127,157 children below five years of age<sup><xref ref-type="bibr" rid="CR36">36</xref></sup>.</p></sec><sec id="Sec5"><title>Vector surveillance</title><p id="Par9">Vector surveillance was conducted in the twelve villages from July 2016 to February 2018. Houses were randomly selected in each village every month for mosquito collections by PSC and indoor CDC light trap (CDC-LT). Household information including roof type, wall type, open or closed eaves, the presence of nets, number of people that slept under a net the previous night and those that did not, and the presence of cattle were collected on a tablet computer. The mosquito density for each method were expressed as the mean number of mosquitoes caught per&#x000a0;house per collection visit.</p><p id="Par10">Indoor-resting mosquitoes were collected between 07:00 and 11:00 by PSC in five houses per site per month. PSCs were done by laying white sheets on the floor and over the furniture within the house. Two collectors, one inside the house and another outside, sprayed around the eaves with 0.025% pyrethrum emulsifiable concentrate mixed with 0.1% piperonyl butoxide (supplied by the Kenya Pyrethrum Board) in kerosene. The collector inside the house then sprayed the roof space. The house was closed for 10&#x02013;15&#x02009;minutes after which knocked-down mosquitoes were collected from the sheets and transferred to the laboratory in scintillation vials containing 70% ethanol.</p><p id="Par11">Indoor host-seeking mosquitoes were collected by CDC-LT in 10 houses per site once per month. A single 12-volt CDC-LT was hung in each house in the sleeping area, approximately 1.5 meters from the floor, adjacent to an occupied bed net owned by a member of the household. The traps were run from 18:00 to 07:00 the following morning. The trapped mosquitoes were transferred into paper cups and transported to the laboratory for further analysis.</p><p id="Par12">Human landing catches (HLC) were used to assess biting time and location (indoor vs outdoor) of the local vector population before and after spraying. HLC was done during the short rains pre-IRS in November 2016, and after the long rains in June 2017. Collections were performed at six sites used for routine surveillance, two in non-IRS areas and four IRS areas. In each site, five houses were randomly selected, and collections were performed for five consecutive nights in each house once before and after IRS.</p><p id="Par13">During HLC, one volunteer sat outside within 5 meters from the house, and another sat inside the house in the living room. Collectors kept their trousers folded to knee length and aspirated any mosquitoes landing on their lower legs. Each house had a team of six collectors, each working in pairs during one of three six-hour shifts running from 17:00 to 11:00 the next morning. Collections were performed for 45&#x02009;minutes, and the collectors rested for 15&#x02009;minutes&#x02009;in each collection hour. The collectors recorded the location of members of the household observed at the end of each hour as either outdoor, in the living room, or in the bedroom. Collected mosquitoes were separated by time and location of collection and sustained on 10% sugar solution before being transported to the laboratory for analysis. Estimation of exposure of individuals to bites by <italic>An. funestus</italic> was performed using models previously described by Seyoum <italic>et al</italic>.<sup><xref ref-type="bibr" rid="CR37">37</xref></sup>.</p></sec><sec id="Sec6"><title>Persistence of insecticidal activity on sprayed walls</title><p id="Par14">To assess the persistence of insecticidal activity on sprayed walls following IRS, WHO cone bioassays<sup><xref ref-type="bibr" rid="CR38">38</xref></sup> were conducted each month using laboratory-reared, 2&#x02013;5 day old, non-blood fed susceptible colony of <italic>An. gambiae</italic> s.s. Kisumu strain. Mosquitoes were exposed in 10 randomly selected sprayed houses, seven with mud walls and three with cement walls, in each of four sub-counties in Migori county. Exposures were performed with ten mosquitoes per cone each month in the same houses at three different heights (0.5&#x02009;m, 1&#x02009;m, and 1.5&#x02009;m) from the floor for 30&#x02009;minutes, on three different walls of the living room of each sprayed house. A control cone with ten mosquitoes was set on an unsprayed plywood board outside of each sprayed house in a shaded area close to the house. Temperature and relative humidity were recorded at every house where mosquitoes were exposed.</p></sec><sec id="Sec7"><title>Insecticide resistance monitoring</title><p id="Par15">WHO insecticide susceptibility tests were performed in Rongo, Nyatike, Awendo and Uriri sub-counties in Migori County (IRS sites) and Homa Bay and Ndhiwa sub-counties in Homa Bay County (no IRS). Larval stages of <italic>An. gambiae</italic> s.l. were collected from Homa Bay, Ndhiwa, Rongo, Nyatike and Awendo sub-counties before IRS and in Homa Bay, Ndhiwa and Rongo after IRS. Larval samples were collected mainly from borrow pits, tire tracks, farmlands and trenches. The collected larvae were raised to three-day-old adults before testing. Adult <italic>An. funestus</italic> were also collected by hand aspiration inside houses for insecticide resistance tests, as larvae were difficult to find. Collections were performed in Homa Bay and Ndhiwa, Rongo, Awendo and Uriri sub-counties before IRS. However, after IRS few adult mosquitoes were in the sprayed sites, so no <italic>An. funestus</italic> s.l. were available for testing from these areas. CDC bottle intensity assay was used to assess intensity of resistance for <italic>An. arabiensis</italic> from Awendo, Ndhiwa and Homa Bay counties. Three-day-old adults raised from field-collected larvae were exposed in bottle coated with doses of 1X, 2X, 5X and 10X of permethrin. Survival frequency of <italic>An. arabiensis</italic> to the different doses was determined at 30&#x02009;minutes exposure period.</p><p id="Par16">Insecticide resistance status was assessed using the WHO diagnostic concentrations of deltamethrin (0.05%), permethrin (0.75%), pirimiphos-methyl (0.25%) and alpha-cypermethrin (0.05%). All papers were prepared by the WHO collaborating centre, Universiti Sains Malaysia. The WHO bioassay was done using 2- to 5-day-old <italic>An. gambiae</italic> s.l. emerging from collected larvae or by direct exposure of field-collected adult <italic>An. funestus</italic> since these were difficult to collect as larvae. At least 100 mosquitoes (four replicates of 25) of each species were exposed to each insecticide per sub-county. The samples were then transferred to a holding tube, provided with cotton wool soaked in 10% sugar solution and held for 24&#x02009;hours. Mortality was scored 24&#x02009;hours after exposure.</p></sec><sec id="Sec8"><title>Mosquito species identification, sporozoite infection and blood meal identification</title><p id="Par17">All <italic>Anopheles</italic> collected were identified morphologically to species using the keys of Gillies and DeMeillon or Gillies and Coetzee<sup><xref ref-type="bibr" rid="CR39">39</xref>,<xref ref-type="bibr" rid="CR40">40</xref></sup>. The physiological status was determined by observation of the abdomen to classify female mosquitoes as either blood-fed, gravid, half gravid or unfed. Female mosquitoes were dissected into three parts for various procedures: heads and thoraces were used for determination of <italic>Plasmodium falciparum</italic> sporozoite infection by enzyme-linked immunosorbent assay (ELISA) using the MR4 Methods in <italic>Anopheles</italic> Research adapted from Wirtz <italic>et al</italic>.<sup><xref ref-type="bibr" rid="CR41">41</xref>,<xref ref-type="bibr" rid="CR42">42</xref></sup>; the abdomens of blood-fed females were used to determine the source of mosquito blood meals by targeting <italic>cytochrome b</italic> protein using a multiplexed PCR protocol<sup><xref ref-type="bibr" rid="CR43">43</xref></sup>, with slight modifications. The legs and wings were used in PCR analyses to identify to species level members of the <italic>An. gambiae</italic> species complex and <italic>Anopheles funestus</italic> group<sup><xref ref-type="bibr" rid="CR44">44</xref></sup>. All mosquitoes morphologically identified as <italic>An. gambiae</italic> s.l. and 20% of randomly selected <italic>An. funestus</italic> s.l., were analyzed by PCR each month. This approach was done due to the greater number of <italic>An. funestus</italic> collected and based on previous studies in the area showing that <italic>An. gambiae</italic> and <italic>An. arabiensis</italic> are found in sympatry, while <italic>An. funestus</italic> s.s. was the only member of the species group routinely collected<sup><xref ref-type="bibr" rid="CR33">33</xref>,<xref ref-type="bibr" rid="CR34">34</xref></sup>.</p></sec><sec id="Sec9"><title>Health facility surveillance</title><p id="Par18">Health facility laboratory data were collected from Rongo, Uriri, and Macalder sub-county hospitals within Migori County (IRS) and Marindi health centre and Ndhiwa sub-County hospital in Homa Bay County (No IRS). The facilities were chosen based on proximity to entomological surveillance sites, availability of health records and catchment area as falling within either IRS or non-IRS area. Febrile cases were tested by health facility staff using light microscopy as part of routine health care and data were recorded in registers provided by the Kenya Ministry of Health. Data were abstracted from laboratory registers of the selected health facilities for the period from January 2015 until June 2018. Each page of the register was photographed using a smartphone camera, and the photographs converted to PDF files using CamScanner-Phone PDF creator, (INTSIG Information Co., Ltd). To ensure confidentiality, the column containing the patient&#x02019;s name was covered when taking the photograph. The PDF copies were then printed and filed.</p></sec><sec id="Sec10"><title>Data management and analysis</title><p id="Par19">Field entomological data collection used Open Data Kit software (ODK) run on tablets with an interface designed to limit data entry errors. Data entry screens used drop-down menus and automatic data checks to reduce errors. Each house sampled received a unique code and a study number. Individual mosquitoes were placed in Eppendorf tubes labeled with pre-printed barcodes and linked to the field data by house code and a unique study number. Results of additional testing, including sporozoite ELISA, species identification by PCR and blood meal analysis, were linked to an&#x000a0;individual mosquito by the unique barcode label. Individual patient records included date of testing, age, gender, village, clinical diagnosis, test performed, and test results from scanned copies of health facility registers. Data were entered into a Microsoft Access database.</p><p id="Par20">Data analysis was performed using R statistical software version 3.4.1 or SAS version 9.4. The risk ratio (RR) was used to assess the statistical significance of differences in mosquito densities pre and post IRS, between intervention and non-intervention sites. Data were fitted using Generalized Linear Mixed Effects Statistical Models (GLMMs). Since the data were over-dispersed, we used the package Generalized Linear Mixed Models using Template Model Builder (glmmTMB) or PROC GLIMMIX, to fit negative binomial distribution models for the analysis of mosquito numbers. The numbers of female <italic>Anopheles</italic> mosquitoes were assessed as a function of the period of collection (before or after IRS) and intervention status (sprayed or non-sprayed) as a fixed effect, while village was treated as a random effect. To obtain the risk ratios (RR) and confidence intervals, we exponentiated the model coefficients. Models were adjusted for reported net use, the presence of open eaves, and the presence of cattle on the compound. A test of interaction was performed to compare differences in estimates of mosquito numbers between the period of mosquito collection and intervention status<sup><xref ref-type="bibr" rid="CR45">45</xref></sup>. Conditional estimates of the change in mosquito densities pre- and post-IRS conditional on the IRS or non-IRS County were generated. To analyse <italic>Anopheles</italic> species proportions pre-and post-IRS in intervention and non-intervention areas for each trapping method, a binomial GLM model was used. The model was also used to analyse sporozoite rates (proportion of sporozoite ELISA tests that are positive) and human biting rates between intervention and non-intervention sites, before and after IRS and proportions of the types of mosquito host blood meals.</p><p id="Par21">To detect changes in numbers of malaria cases before and after IRS within each health facility Auto-Regressive Integrated Moving Average (ARIMA) analysis was performed. Data from each facility was analysed using the &#x0201c;Time Series Analysis&#x0201d; (TSA)<sup><xref ref-type="bibr" rid="CR46">46</xref></sup> and &#x0201c;Alternative Time Series Analyses&#x0201d; (aTSA)<sup><xref ref-type="bibr" rid="CR47">47</xref>,<xref ref-type="bibr" rid="CR48">48</xref></sup> packages in R. The ARIMA model was derived by observation of the autocorrelation and partial-autocorrelation functions to determine the most parsimonious solution of the &#x0201c;order&#x0201d; (p), &#x0201c;differencing&#x0201d; (d), and &#x0201c;moving-average&#x0201d; (q) parameter values. The model was then regressed on the absence (prior to) or presence of IRS in the village to estimate the value of the number of positive malaria cases prior to, and during the period of IRS.</p></sec><sec id="Sec11"><title>Ethical considerations</title><p id="Par22">The study was approved by the Kenya Medical Research Institute/ Scientific and Ethics Review Unit (KEMRI/SERU), number 2776 and by CDC through a reliance agreement with KEMRI/SERU (CDC IRB 6728). Individuals participating in HLC gave informed consent. They were screened for malaria before the start of the study and treated if positive. Collectors were placed on mefloquine malaria prophylaxis, (Mephaquin, Acino Pharma AG, Switzerland) one week before collections began, with repeat doses once every week through the collection period, until four weeks after collections ended. During routine mosquito collections, verbal consent was sought from the household head to use CDC-LT and PSC in their compound. All methods were performed in accordance with relevant guidelines and regulations.</p></sec></sec></sec><sec id="Sec12" sec-type="results"><title>Results</title><sec id="Sec13"><title>Vector species composition and seasonality</title><p id="Par23">The mean number of <italic>An. funestus</italic> and <italic>An. gambiae</italic>. s.l.&#x000a0;found in indoor CDC-LT and PSCs are presented by IRS status and period (pre- or post-IRS) in Table&#x000a0;<xref rid="Tab1" ref-type="table">1</xref>. The number of each species of mosquito collected by the two different methods was compared using negative binomial regression models incorporating IRS status, period, an interaction between IRS status and period, net use, the presence of open eaves and the presence of cattle on the compound (Supplemental Tables&#x000a0;<xref rid="MOESM1" ref-type="media">1</xref>&#x02013;<xref rid="MOESM1" ref-type="media">4</xref>). For all models except for the <italic>An. gambiae</italic> s.l. collected by PSC, the interaction term was statistically significant indicating a differential effect of period based upon the IRS status. Conditional estimates of the effect of period controlling for IRS status are provided in Table&#x000a0;<xref rid="Tab1" ref-type="table">1</xref>.<table-wrap id="Tab1"><label>Table 1</label><caption><p>Comparison of mean numbers of <italic>An. funestus</italic> and <italic>An. arabiensis</italic> collected indoors by CDC-LTs and PSCs pre- and post-IRS in intervention and non-intervention areas. Risk ratios of post- versus pre-IRS periods conditional on intervention status are also provided for each species and collection method. Models include terms for IRS status, pre/post spray period and an interaction term (Table 1). See Supplemental Tables 1&#x02013;4 for full models.</p></caption><table frame="hsides" rules="groups"><thead><tr><th><italic>Anopheles</italic> Species</th><th>Collection Method</th><th>IRS Status</th><th>Level</th><th>Mean</th><th>Risk Ratio</th><th>Lower CL</th><th>Upper CL</th><th>P-value</th></tr></thead><tbody><tr><td rowspan="8"><italic>Anopheles funestus</italic></td><td rowspan="4">Light trap</td><td rowspan="2">IRS</td><td>Post-Spray</td><td>0.05</td><td>0.12</td><td>0.07</td><td>0.19</td><td><bold>&#x0003c;0.001</bold></td></tr><tr><td>Pre-Spray</td><td>0.45</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="2">Non-IRS</td><td>Post-Spray</td><td>0.88</td><td>0.98</td><td>0.69</td><td>1.38</td><td>0.899</td></tr><tr><td>Pre-Spray</td><td>0.92</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="4">PSC</td><td rowspan="2">IRS</td><td>Post-Spray</td><td>0.04</td><td>0.04</td><td>0.02</td><td>0.07</td><td><bold>&#x0003c;0.001</bold></td></tr><tr><td>Pre-Spray</td><td>0.99</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="2">Non-IRS</td><td>Post-Spray</td><td>1.05</td><td>0.64</td><td>0.41</td><td>1.00</td><td>0.052</td></tr><tr><td>Pre-Spray</td><td>2.05</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="8"><italic>Anopheles arabiensis</italic></td><td rowspan="4">Light trap</td><td rowspan="2">IRS</td><td>Post-Spray</td><td>0.19</td><td>1.39</td><td>0.78</td><td>2.47</td><td>0.266</td></tr><tr><td>Pre-Spray</td><td>0.10</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="2">Non-IRS</td><td>Post-Spray</td><td>0.21</td><td>3.06</td><td>1.59</td><td>5.92</td><td><bold>0.001</bold></td></tr><tr><td>Pre-Spray</td><td>0.05</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="4">PSC</td><td rowspan="2">IRS</td><td>Post-Spray</td><td>0.24</td><td>0.60</td><td>0.33</td><td>1.09</td><td>0.093</td></tr><tr><td>Pre-Spray</td><td>0.52</td><td>Ref</td><td/><td/><td/></tr><tr><td rowspan="2">Non-IRS</td><td>Post-Spray</td><td>0.41</td><td>1.64</td><td>0.87</td><td>3.09</td><td>0.123</td></tr><tr><td>Pre-Spray</td><td>0.27</td><td>Ref</td><td/><td/><td/></tr></tbody></table></table-wrap></p><p id="Par24">The number of <italic>An. funestus</italic> collected in light traps in intervention sites were significantly lower in the post-IRS compared to pre-IRS period (RR&#x02009;=&#x02009;0.12, 95% CI: 0.07&#x02013;0.19, P&#x02009;&#x0003c;&#x02009;0.001). No significant difference in the mean number of <italic>An. funestus</italic> was observed in the non-intervention sites between pre- and post-IRS (IRR&#x02009;=&#x02009;0.98, 95% CI: 0.69&#x02013;1.38, p&#x02009;=&#x02009;0.899). A statistically significant difference-of-differences between the period of mosquito collection and intervention status was observed based on the statistically significant interaction term (RR&#x02009;=&#x02009;0.12, 95% CI: 0.07&#x02013;0.21) (Supplemental Table&#x000a0;<xref rid="MOESM1" ref-type="media">1</xref>). From PSC collections, significantly fewer numbers of <italic>An. funestus</italic> were observed in both IRS and non-IRS sites in the post-IRS period compared to pre-IRS period (RR&#x02009;=&#x02009;0.04, 95% CI: 0.02&#x02013;0.07, p&#x02009;&#x0003c;&#x02009;0.001). The number of <italic>An. funestus</italic> in the non-IRS area also declined but the conditional difference between pre-IRS and post-IRS was not statistically significant (RR&#x02009;=&#x02009;0.64, 95% CI: 0.41&#x02013;1.00, p&#x02009;=&#x02009;0.052). A statistically significant difference-of-differences was observed between period of mosquito collection and intervention status post-IRS indicating a stronger decline in the IRS sites compared to the non-IRS sites (RR&#x02009;=&#x02009;0.06, 95% CI: 0.03&#x02013;0.13, p&#x02009;&#x0003c;&#x02009;0.001) (Supplemental Table&#x000a0;<xref rid="MOESM1" ref-type="media">2</xref>).</p><p id="Par25">The mean numbers of <italic>An. arabiensis</italic> collected in indoor CDC-LTs in both intervention and non-intervention sites increased in the post-IRS compared to pre-IRS period with a statistically significant increase in the non-IRS sites (IRS sites: RR&#x02009;=&#x02009;1.39, 95% CI: 0.78&#x02013;2.47, p&#x02009;=&#x02009;0.266; non-IRS sites: RR&#x02009;=&#x02009;3.06, 95% CI: 1.59&#x02013;5.92, p&#x02009;=&#x02009;0.001). The conditional estimates are provided in Table&#x000a0;<xref rid="Tab1" ref-type="table">1</xref> although the interaction term was not significant (RR&#x02009;=&#x02009;0.45, 95% CI: 0.2&#x02013;1.01, p&#x02009;=&#x02009;0.052) (Table&#x000a0;<xref rid="Tab1" ref-type="table">1</xref>) indicating the increase was not statistically greater in the non-IRS sites than the IRS sites (Supplemental Table&#x000a0;<xref rid="MOESM1" ref-type="media">3</xref>).</p><p id="Par26">The mean numbers of <italic>An. arabiensis</italic> collected by PSC in the intervention sites were not significantly different in the post-IRS compared to pre-IRS period (RR&#x02009;=&#x02009;0.60, 95% CI: 0.33&#x02013;1.09, p&#x02009;=&#x02009;0.093). For the non-IRS areas, the number of <italic>An. arabiensis</italic> collected by PSC increased although not significantly (RR&#x02009;=&#x02009;1.64, 95% CI: 0.87&#x02013;3.09, p&#x02009;=&#x02009;0.123). Although no significant difference in the mean numbers of <italic>An. arabiensis</italic> was observed pre- and post-IRS in either the IRS or the non-IRS sites, a statistically significant difference-of-differences was observed between the period of mosquito collection and intervention status indicating a significant difference between the IRS and non-IRS areas after IRS implementation (RR&#x02009;=&#x02009;0.36, 95% CI: 0.16&#x02013;0.82, p&#x02009;=&#x02009;0.015) (Supplemental Table&#x000a0;<xref rid="MOESM1" ref-type="media">4</xref>).</p><p id="Par27">The mean number of <italic>An. funestus</italic> collected by CDC-LT and PSC in both intervention and non-intervention areas varied by month, with the highest numbers collected during the short rainy season before IRS (Nov&#x02013;Dec 2016) and during the long rainy season in the unsprayed area (March&#x02013;June 2017) (Fig.&#x000a0;<xref rid="Fig2" ref-type="fig">2</xref>). After IRS, the mean numbers collected by both CDC-LT and PSC in the intervention areas remained low, with no seasonal variation throughout the study period. The mean number of <italic>An. arabiensis</italic> collected by either method was lower compared to <italic>An. funestus</italic> with little monthly variation before and after IRS. No clear difference was observed in the seasonality of <italic>An. arabiensis</italic> before and after IRS (Fig.&#x000a0;<xref rid="Fig2" ref-type="fig">2</xref>).<fig id="Fig2"><label>Figure 2</label><caption><p>Mean number of observed female <italic>Anopheles</italic> (means&#x02009;&#x000b1;&#x02009;std error) per trap-night per month in indoor CDC light trap and PSC before and after IRS in sprayed and unsprayed areas. The grey shaded area is the&#x000a0;period post-spray with residual efficacy above 80%. The primary scale shows <italic>Anopheles</italic> density while the secondary scale shows rainfall in milliliters. Sixty trap-nights for CDC-LT and 30 for PSC per study arm per month. Standard errors were calculated by dividing the standard deviation by the square root of the sample size, a function was developed in R that generated the standard error around each mean.</p></caption><graphic xlink:href="41598_2020_61350_Fig2_HTML" id="d29e1295"/></fig></p><p id="Par28">Differences in <italic>Anopheles</italic> species composition were observed in intervention areas before and after IRS. In both CDC-LT and PSC collections, <italic>An. funestus</italic> comprised over 80% of the total <italic>Anopheles</italic> collected in both intervention and non-intervention sites before IRS. While <italic>An. funestus</italic> remained predominant in non-intervention sites after IRS, <italic>An. arabiensis</italic> formed the bulk of all <italic>Anopheles</italic> collected in the intervention sites after IRS (Fig.&#x000a0;<xref rid="Fig3" ref-type="fig">3</xref>).<fig id="Fig3"><label>Figure 3</label><caption><p>Proportions of <italic>Anopheles</italic> species collected by CDC-LT and PSC before and after IRS in sprayed and unsprayed regions.</p></caption><graphic xlink:href="41598_2020_61350_Fig3_HTML" id="d29e1331"/></fig></p><sec id="Sec14"><title>Insecticide decay rate and insecticide resistance monitoring</title><p id="Par29">Mortality rates of susceptible <italic>An. gambiae</italic> s.s females were over 80% up to 10 months post-IRS across all study sites&#x000a0;on both mud and concrete walls as measured by WHO cone assay (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">1</xref>). Using WHO tube bioassays both <italic>An. funestus</italic> (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">2A</xref>) and <italic>An. arabiensis</italic> (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">2B</xref>) before IRS and <italic>An. arabiensis</italic> after IRS (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">2C</xref>) were fully susceptible to pirimiphos-methyl and bendiocarb but resistant to the pyrethroids, deltamethrin, permethrin, and alpha-cypermethrin. Survival frequencies of <italic>An. arabiensis</italic> to different doses of permethrin varied between sites (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">3</xref>).</p></sec><sec id="Sec15"><title>Factors affecting <italic>Anopheles</italic> numbers</title><p id="Par30">Table&#x000a0;<xref rid="Tab2" ref-type="table">2</xref> presents data showing modelled estimates of the effect of net use, open eaves, and presence of cattle in the compound on the indoor occurrence of <italic>An. funestus</italic> and <italic>An. arabiensis</italic> in sprayed and unsprayed houses, measured by CDC-LT and PSC collections. For <italic>An. funestus</italic>, significantly fewer were collected by light traps in houses with completely closed eaves (RR&#x02009;=&#x02009;0.68, 95% CI: 0.48&#x02013;0.96, p&#x02009;=&#x02009;0.030) while significantly more were collected from houses where cattle were kept on the compound (RR&#x02009;=&#x02009;1.62, 95%CI: 1.22&#x02013;2.13, p&#x02009;=&#x02009;0.001). No other comparisons were statistically significant. By PSC, there were again significantly more <italic>An. funestus</italic> collected in households where cattle were kept on the compound (RR&#x02009;=&#x02009;1.63, 95% CI: 1.12&#x02013;2.35, p&#x02009;=&#x02009;0.010). There were significantly more <italic>An. funestus</italic> in houses where some but not all residents used a net the previous night compared to houses where no one used a net (RR&#x02009;=&#x02009;2.02, 95% CI: 1.13&#x02013;3.59, p&#x02009;=&#x02009;0.017). No other comparisons were statistically significant.<table-wrap id="Tab2"><label>Table 2</label><caption><p>Model estimates comparing the mean number of indoor <italic>An. funestus</italic> and <italic>An. arabiensis</italic> collected, by collection type, eave type, net use, and presence of cattle in intervention and non-intervention areas. Models include terms for IRS status, pre/post spray period and an interaction term (Table&#x000a0;<xref rid="Tab1" ref-type="table">1</xref>). See Supplemental Tables&#x000a0;<xref rid="MOESM1" ref-type="media">1&#x02013;4</xref> for full models.</p></caption><table frame="hsides" rules="groups"><thead><tr><th>Species</th><th>Collection Method</th><th>Parameter</th><th>Level</th><th>Proportion of Houses</th><th>Risk Ratio</th><th>Lower CL</th><th>Upper CL</th><th>t-value</th><th>P-value</th></tr></thead><tbody><tr><td rowspan="16"><italic>Anopheles funestus</italic></td><td rowspan="8">Light Trap</td><td rowspan="3">Net Use</td><td>All under net</td><td>0.47</td><td>1.11</td><td>0.77</td><td>1.6</td><td>0.576</td><td>0.565</td></tr><tr><td>Some under net</td><td>0.09</td><td>1.2</td><td>0.77</td><td>1.86</td><td>0.816</td><td>0.415</td></tr><tr><td>None under net</td><td>0.44</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="3">Eaves</td><td>Closed</td><td>0.20</td><td>0.68</td><td>0.48</td><td>0.96</td><td>&#x02212;2.174</td><td><bold>0.030</bold></td></tr><tr><td>Partially open</td><td>0.10</td><td>0.84</td><td>0.56</td><td>1.27</td><td>&#x02212;0.817</td><td>0.414</td></tr><tr><td>Open</td><td>0.70</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="2">Cattle</td><td>Yes</td><td>0.71</td><td>1.62</td><td>1.22</td><td>2.13</td><td>3.395</td><td><bold>0.001</bold></td></tr><tr><td>No</td><td>0.29</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="8">PSC</td><td rowspan="3">Net Use</td><td>All under net</td><td>0.43</td><td>0.96</td><td>0.61</td><td>1.5</td><td>&#x02212;0.187</td><td>0.852</td></tr><tr><td>Some under net</td><td>0.08</td><td>2.02</td><td>1.13</td><td>3.59</td><td>2.383</td><td><bold>0.017</bold></td></tr><tr><td>None under net</td><td>0.49</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="3">Eaves</td><td>Closed</td><td>0.18</td><td>0.8</td><td>0.5</td><td>1.3</td><td>&#x02212;0.889</td><td>0.374</td></tr><tr><td>Partially open</td><td>0.12</td><td>1.08</td><td>0.64</td><td>1.83</td><td>0.291</td><td>0.771</td></tr><tr><td>Open</td><td>0.70</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="2">Cattle</td><td>Yes</td><td>0.71</td><td>1.63</td><td>1.12</td><td>2.35</td><td>2.583</td><td><bold>0.010</bold></td></tr><tr><td>No</td><td>0.29</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="16"><italic>Anopheles arabiensis</italic></td><td rowspan="8">Light Trap</td><td rowspan="3">Net Use</td><td>All under net</td><td>0.47</td><td>1.95</td><td>0.99</td><td>3.84</td><td>1.94</td><td>0.052</td></tr><tr><td>Some under net</td><td>0.09</td><td>2.17</td><td>1.02</td><td>4.62</td><td>2.008</td><td><bold>0.045</bold></td></tr><tr><td>None under net</td><td>0.44</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="3">Eaves</td><td>Closed</td><td>0.20</td><td>0.57</td><td>0.33</td><td>0.96</td><td>&#x02212;2.131</td><td><bold>0.033</bold></td></tr><tr><td>Partially open</td><td>0.09</td><td>0.78</td><td>0.43</td><td>1.42</td><td>&#x02212;0.814</td><td>0.416</td></tr><tr><td>Open</td><td>0.70</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="2">Cattle</td><td>Yes</td><td>0.71</td><td>1.33</td><td>0.89</td><td>1.98</td><td>1.383</td><td>0.167</td></tr><tr><td>No</td><td>0.29</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="8">PSC</td><td rowspan="3">Net Use</td><td>All under net</td><td>0.43</td><td>1.61</td><td>0.9</td><td>2.87</td><td>1.594</td><td>0.111</td></tr><tr><td>Some under net</td><td>0.08</td><td>1.85</td><td>0.89</td><td>3.84</td><td>1.652</td><td>0.099</td></tr><tr><td>None under net</td><td>0.49</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="3">Eaves</td><td>Closed</td><td>0.18</td><td>0.34</td><td>0.18</td><td>0.67</td><td>&#x02212;3.165</td><td><bold>0.002</bold></td></tr><tr><td>Partially open</td><td>0.11</td><td>0.60</td><td>0.32</td><td>1.13</td><td>&#x02212;1.583</td><td>0.114</td></tr><tr><td>Open</td><td>0.70</td><td>Ref</td><td/><td/><td/><td/></tr><tr><td rowspan="2">Cattle</td><td>Yes</td><td>0.71</td><td>1.53</td><td>1.00</td><td>2.34</td><td>1.944</td><td>0.052</td></tr><tr><td>No</td><td>0.29</td><td>Ref</td><td/><td/><td/><td/></tr></tbody></table></table-wrap></p><p id="Par31">From light trap collections, closed eaves were associated with significantly lower numbers of <italic>An. arabiensis</italic> (RR&#x02009;=&#x02009;0.57, 95% CI: 0.33&#x02013;0.96, p&#x02009;=&#x02009;0.033) while significantly more <italic>An. arabiensis</italic> were collected in houses where some but not all residents of the household used a net the night before compared to houses where no one used a net (RR&#x02009;=&#x02009;2.17, 95% CI: 1.02&#x02013;4.62, p&#x02009;=&#x02009;0.045). The number of <italic>An. arabiensis</italic> collected by PSC also was significantly lower in houses with closed eaves compared to those with open eaves (RR&#x02009;=&#x02009;0.34, 95% CI: 0.18&#x02013;0.67, p&#x02009;=&#x02009;0.002). No other comparisons were statistically significant.</p></sec><sec id="Sec16"><title>Sporozoite infection rates</title><p id="Par32">Sporozoite infection rates in <italic>Anopheles</italic> mosquitoes were determined in intervention and non-intervention sites before and after IRS. Before IRS, 4.8% (48/1,000) of <italic>An. funestus</italic> were sporozoite positive in non-intervention sites compared to 2.2% (10/447) in the intervention sites whereas&#x000a0;for <italic>An. arabiensis</italic>, sporozoite positivity rate was 2.8% (10/357) in non-intervention sites and 1.5% (3/192) in the intervention sites before IRS. Sporozoite infection rates for both species combined were not significantly different between intervention and non-intervention sites before IRS (RR&#x02009;=&#x02009;0.27, 95% CI: 0.06&#x02013;1.26, P&#x02009;=&#x02009;0.09). After IRS, sporozoite infections were detected only in the non-intervention sites, where 3.5% (40/1,132) of <italic>An. funestus</italic> and 3.3% (22/643) of <italic>An. arabiensis</italic> were positive. (Fig.&#x000a0;<xref rid="Fig4" ref-type="fig">4</xref>).<fig id="Fig4"><label>Figure 4</label><caption><p>Sporozoite rates in <italic>An. funestus</italic> and <italic>An. arabiensis</italic> in sprayed and unsprayed areas, pre- and post-IRS.</p></caption><graphic xlink:href="41598_2020_61350_Fig4_HTML" id="d29e2181"/></fig></p></sec><sec id="Sec17"><title>Vector biting behaviour</title><p id="Par33">We estimated the exposure of humans to the risk of mosquito bites based on their observed behaviour and time and location of <italic>An. funestus</italic> biting. The numbers of <italic>An. arabiensis</italic> were insufficient to be included in the analysis. Over 70% of people within the study area were observed to be outdoors at 17:00, the beginning of mosquito collection. The number of people outdoors declined steadily over time with an increase in the number of individuals observed indoors, either in the living room (indoors not asleep) or in the bedroom (indoors and in bed). Over 90% of the people were observed to be indoors and in bed between 23:00 and 05:00 (Supplemental Fig.&#x000a0;<xref rid="MOESM1" ref-type="media">4</xref>). In both intervention and non-intervention sites, before IRS, exposure to <italic>An. funestus</italic> was estimated to occur mostly, although not exclusively, indoors, late at night when people were asleep (Fig.&#x000a0;<xref rid="Fig5" ref-type="fig">5a,b</xref>). In the post-IRS period, no change in the estimated exposure to bites by <italic>An. funestus</italic> was observed in the non-intervention sites (Fig.&#x000a0;<xref rid="Fig5" ref-type="fig">5c</xref>). However, in the intervention sites, the risk of exposure to mosquito bites was nearly zero post-IRS (Fig.&#x000a0;<xref rid="Fig5" ref-type="fig">5d</xref>). The relative proportion of bites by <italic>An. funestus</italic> increased both indoors and outdoors at dawn (05:00 am&#x02013;08:00 am), corresponding to the time when most individuals woke up. Low levels of biting continued until 11:00 am when collection ceased.<fig id="Fig5"><label>Figure 5</label><caption><p>Profiles of biting by <italic>An. funestus</italic> experienced by the human population in intervention and non-intervention sites before and after IRS. The black area represents biting that occurs outdoors, the dark red represents biting that occurs indoors when people are away from their bed nets and the blue represents biting that occurs while people are asleep.</p></caption><graphic xlink:href="41598_2020_61350_Fig5_HTML" id="d29e2226"/></fig></p></sec><sec id="Sec18"><title>Blood meal type</title><p id="Par34">Blood meal analysis using mosquitoes collected by PSC was conducted on 236 fed <italic>Anopheles</italic> mosquitoes, 151 <italic>An. funestus</italic> and 85 <italic>An. arabiensis</italic>. <italic>An. funestus</italic> fed mostly on humans 52.3% (79/151), followed by cattle 40.4% (61/151), goat 3.3% (5/151), pig 0.7% (1/151) and mixed-blood meals 3.3% (5/151, 2 human/cow, 2 human/goat and 1 human/pig). <italic>An. arabiensis</italic> had fed mostly on cattle blood 70.6% (60/85), followed by pig 12.9% (11/85), human 9.4% (8/85), goat 4.7% (4/85) and mixed-blood meal, human/goat 1.9% (1/85) (Fig.&#x000a0;<xref rid="Fig6" ref-type="fig">6</xref>).<fig id="Fig6"><label>Figure 6</label><caption><p>Comparison of mammalian host blood meal type between <italic>An</italic>. <italic>funestus</italic> and <italic>An. arabiensis</italic> (Numbers tested; <italic>An. funestus</italic>- 61 cow, 5 goats, 79 humans, 1 pig and 5 mixed and <italic>An. arabiensis</italic> &#x02013; 60 cows, 4 goats, 8, humans, 11 pigs and 3 mixed).</p></caption><graphic xlink:href="41598_2020_61350_Fig6_HTML" id="d29e2275"/></fig></p></sec><sec id="Sec19"><title>Malaria case count</title><p id="Par35">A total of 137,972 laboratory test records from patients attending the out-patient departments were extracted from the five health facilities. For the two-year period before IRS (January 2015&#x02013;February 2017), malaria test positive proportions were similar at 33.2% (18,036/54,404) in intervention and 33.3% (12,920/38,835) in non-intervention sites respectively. For the post-IRS period (March 2017&#x02013;May 2018), the test positivity rates were 30.4% (6,347/20,882) in the non-intervention sites and 20.6% (4905/23,851) in the intervention sites.</p><p id="Par36">ARIMA analysis of malaria case counts for health facilities within IRS areas showed a reduction in malaria cases at the facilities post-IRS. Estimated mean monthly malaria cases in Rongo sub-county hospital dropped by 44% from 323 cases per month before IRS to 178 after IRS [mean difference&#x02009;=&#x02009;&#x02212;142; 95% CI: &#x02212;236 to &#x02212;48; P&#x02009;=&#x02009;0.003]. A similar reduction in mean monthly malaria cases with&#x000a0;a 65.0% drop from 301 before IRS to 78 cases after IRS [mean difference&#x02009;=&#x02009;&#x02212;196; 95% CI: &#x02212;345 to &#x02212;47; P&#x02009;=&#x02009;0.01] was observed in Uriri sub-county hospital. In Nyatike sub-county hospital, the mean monthly malaria cases dropped by 47.4% from 118 cases before IRS to 72 after IRS [mean difference&#x02009;=&#x02009;&#x02212;56; 95% CI: &#x02212;123 to 11; P&#x02009;=&#x02009;0.1]. For the two health facilities within non-IRS sites, no significant changes in malaria case counts were observed post-IRS, [Ndhiwa hospital: mean difference&#x02009;=&#x02009;&#x02212;82; 95%CI: &#x02212;230 to 65; P&#x02009;=&#x02009;0.3; Marindi hospital: mean difference = 9.3; 95%CI: &#x02212;132 to 151; = 0.9]. A plot of positive malaria cases over time, before and after IRS, showed a decline in the number of cases detected at facilities within sprayed areas compared to those in unsprayed regions (Fig.&#x000a0;<xref rid="Fig7" ref-type="fig">7</xref>).<fig id="Fig7"><label>Figure 7</label><caption><p>Health facility laboratory test positivity rates among febrile out-patients in Homa Bay (non-intervention) and Migori (Intervention) covering two years pre-IRS and over one-year post-IRS.</p></caption><graphic xlink:href="41598_2020_61350_Fig7_HTML" id="d29e2294"/></fig></p></sec></sec></sec><sec id="Sec20" sec-type="discussion"><title>Discussion</title><p id="Par37">Our results demonstrate a significant reduction in <italic>An. funestus</italic> indoor resting densities, biting rates and sporozoite infections, as well as a decline in malaria test positivity rates and case-count at health facilities after one round of Actellic 300CS IRS in Migori County, western Kenya. Human biting rates and sporozoite infections in <italic>Anopheles</italic> mosquitoes are the most direct entomological measures of malaria infection risk. We observed moderate biting and sporozoite rates in both intervention and non-intervention sites before IRS and the unsprayed sites after IRS. However, after IRS, <italic>An. funestus</italic> biting rates were nearly zero, and no sporozoite infections were detected post-IRS. Susceptibility tests confirmed that the major vector species, <italic>An. arabiensis</italic> and <italic>An. funestus</italic>, were both resistant to pyrethroid insecticides but were susceptible to pirimiphos-methyl, the active ingredient in Actellic 300CS.</p><p id="Par38">Similar reductions in <italic>An. funestus</italic> populations to near elimination were observed in the Asembo Bay area of western Kenya following the scale-up of pyrethroid-treated nets<sup><xref ref-type="bibr" rid="CR5">5</xref></sup> although <italic>An. funestus</italic> later returned as the primary malaria vector in the region presumably due to the development of pyrethroid resistance<sup><xref ref-type="bibr" rid="CR34">34</xref></sup>. The complete elimination of <italic>An. funestus</italic> following effective IRS campaigns have been reported in South Africa, Mauritius and the Pare/Taveta area of Tanzania and Kenya<sup><xref ref-type="bibr" rid="CR49">49</xref>,<xref ref-type="bibr" rid="CR50">50</xref></sup>. <italic>An. funestus</italic> is particularly sensitive to effective insecticides used in IRS, probably due to the highly endophilic and anthropophilic nature of the species<sup><xref ref-type="bibr" rid="CR51">51</xref>&#x02013;<xref ref-type="bibr" rid="CR53">53</xref></sup>, leading to high levels of exposure to treated surfaces.</p><p id="Par39">In contrast, <italic>An. arabiensis</italic> indoor resting densities, human-biting rates, and sporozoite infection rates all reduced only marginally in sprayed areas post-IRS. With the decline in <italic>An. funestus</italic>, <italic>An</italic>. <italic>arabiensis</italic> became the predominant vector species in the sprayed areas (although the densities of <italic>An. arabiensis</italic> did not increase). IRS has a limited impact on the population of <italic>An. arabiensis</italic>, despite full susceptibility to pirimiphos-methyl in WHO susceptibility tests. This lesser impact of IRS on <italic>An. arabiensis</italic> is therefore unlikely to be due to insecticide resistance but may be attributable to the behaviour of this species. Blood-meal host analysis showed that <italic>An. arabiensis</italic> fed more frequently on cattle than humans, unlike <italic>An. funestus</italic> that fed more frequently on humans. This finding is supported by a previous study in western Kenya that reported <italic>An. arabiensis</italic> fed predominantly on cattle (65% of blood meals on cattle; 22% mixed bovine/human; 13% human)<sup><xref ref-type="bibr" rid="CR8">8</xref></sup>. Furthermore, results from a deterministic model, developed using data from Kilombero, Tanzania, suggested that <italic>An. arabiensis</italic> fed outdoors on both humans and cattle and rapidly exited houses without fatal exposure to insecticidal nets or IRS<sup><xref ref-type="bibr" rid="CR54">54</xref></sup>. A recent trap evaluation in western Kenya observed over sevenfold more <italic>An. arabiensis</italic> collected by cow odour compared to human odour outdoor<sup><xref ref-type="bibr" rid="CR55">55</xref></sup>. Therefore, it is likely that a significant population of <italic>An. arabiensis</italic> rests predominantly outdoors and feeds primarily on cattle, but occasionally bites humans and transmits malaria, albeit less efficiently than the more anthropophilic vector <italic>An. funestus</italic>. Therefore, it is possible that the population of <italic>An. arabiensis</italic> collected indoors by light traps and PSC represents only a proportion of a larger outdoor population. These factors may explain the lesser impact of IRS on <italic>An. arabiensis</italic>.</p><p id="Par40">IRS with Actellic 300CS had a prolonged residual activity of at least ten months post-IRS, as measured by wall bioassays. As spraying was conducted in February, the insecticide provided protection throughout the periods of highest malaria transmission during the long (April&#x02013;June) and short (October&#x02013;November) rainy seasons. This long residual efficacy of Actellic 300CS makes the insecticide particularly useful in providing all-year-round protection with just one spray round each year. Similar prolonged residual activity of Actellic 300CS and control of pyrethroid-resistant mosquitoes have been reported in other countries<sup><xref ref-type="bibr" rid="CR32">32</xref></sup>. As this is the first report of the efficacy of IRS with a non-pyrethroid insecticide in Kenya in the wake of pyrethroid resistance, strict adherence to the insecticide resistance management strategy<sup><xref ref-type="bibr" rid="CR27">27</xref>,<xref ref-type="bibr" rid="CR28">28</xref></sup> is needed. The Kenya National Malaria Control Program (NMCP) has a responsibility to fast track registration of additional non-pyrethroid classes of insecticides to enable rotation of insecticide in IRS for effective resistance management.</p><p id="Par41">This study presents three novel findings that are crucial for malaria control and elimination in Kenya. Our study demonstrated that a single spray-round with Actellic 300CS is effective for year-round control of malaria in the region. Well-timed IRS, just before the peak malaria transmission, suppressed the local malaria vector populations through both the long and short rain periods, hence reducing disease transmission even during the peak transmission seasons that coincide with the long and short rain seasons. Second, the impact of IRS on <italic>An. arabiensis</italic> was less than that on other vectors, making it the dominant malaria vector sustaining residual transmission. While IRS was effective in controlling <italic>An. funestus, An. arabiensis</italic> would require alternative control approaches for malaria elimination. Finally, the study demonstrates late morning biting by <italic>An. funestus</italic>, a phenomenon not before investigated in the region. This presents a new challenge facing vector control given that the tools currently deployed for vector control are limited to the household and mostly only to the times when people are asleep (LLINs).</p><p id="Par42">Biting by <italic>An. funestus</italic> in the intervention and non-intervention areas before IRS and non-intervention areas after IRS occurred mostly indoors late at night corresponding to the period when most people were indoors and in bed. Late night, indoor biting by <italic>An. funestus</italic> has been previously reported, dating back to the pre-bed net era. For instance, in 1975, 94% of <italic>An. funestus</italic> were observed to bite after midnight, with another a peak in the hours before dawn in Kano plain of western Kenya<sup><xref ref-type="bibr" rid="CR56">56</xref></sup>. Similar late night, indoor biting was observed in the same study area in 1996<sup><xref ref-type="bibr" rid="CR53">53</xref></sup> and more recently the vector species have been reported to persistently bite indoors, late in the night despite high coverage in insecticide-treated nets<sup><xref ref-type="bibr" rid="CR52">52</xref></sup>. While all these studies observed a biting peak at dawn, the collections ceased at 07:00&#x02009;hours. With extended collections to monitor <italic>An. funestus</italic> biting in the morning, we observed biting until 11:00&#x02009;hours. Similar findings of day-biting <italic>An. funestus</italic> in the presence of LLINs have been recently reported in Senegal<sup><xref ref-type="bibr" rid="CR9">9</xref></sup> and Benin<sup><xref ref-type="bibr" rid="CR57">57</xref></sup>. This observed extended biting behaviour in <italic>An. funestus</italic> not previously investigated may potentially undermine the effectiveness of LLINs as people may be exposed to mosquito bites while away from the protection of their bed nets. However, one round of IRS with Actellic 300CS substantially reduced the number of <italic>An. funestus</italic> collected and it was not possible to detect biting either indoors or outdoors in the sprayed areas post-IRS.</p><p id="Par43">In an analysis of risk factors associated with the indoor occurrence of mosquitoes, fewer <italic>An. funestus</italic> and <italic>An. arabiensis</italic> were collected by light traps in houses with closed eaves. Similar results were observed for <italic>An. arabiensis</italic> collected by PSC. Open eaves are known to be the main route for indoor entry of <italic>Anopheles</italic> mosquitoes<sup><xref ref-type="bibr" rid="CR11">11</xref>,<xref ref-type="bibr" rid="CR58">58</xref></sup> and blocking them has been demonstrated to be effective in preventing <italic>Anopheles</italic> house entry<sup><xref ref-type="bibr" rid="CR59">59</xref></sup>. Closing eave spaces or deploying vector control tools in these spaces may present an additional intervention to the current vector control tool kit for reducing the indoor occurrence of mosquitoes in addition to IRS and LLINs. However, new tools targeting outdoor <italic>Anopheles</italic> populations, in particular, <italic>An. arabiensis</italic>, are needed to control vector species that are less well controlled by interventions that target the indoor environment such as IRS and LLINs.</p><p id="Par44">Reductions in malaria cases at the health facilities within sprayed areas post-IRS provided further evidence of the impact of a single round of IRS on malaria transmission. Health facility-based surveys of malaria cases in febrile patients have been useful as part of rapid analysis of changes in local malaria epidemiology<sup><xref ref-type="bibr" rid="CR60">60</xref>&#x02013;<xref ref-type="bibr" rid="CR63">63</xref></sup>. Malaria infection is highly correlated with febrile cases reported at the health facilities<sup><xref ref-type="bibr" rid="CR62">62</xref></sup>. Furthermore, a systematic review of febrile illness over 20 years in sub-Saharan Africa reported a dramatic reduction in the proportion of fevers associated with <italic>Plasmodium falciparum</italic> malaria<sup><xref ref-type="bibr" rid="CR64">64</xref></sup>. Consequently, reductions in malaria cases likely contribute considerably to the reductions in febrile illnesses presenting at health facilities. The use of routine Health Management Information System (HMIS) data to evaluate malaria control interventions<sup><xref ref-type="bibr" rid="CR63">63</xref></sup>, however, suffers from incompleteness in reporting and variation in the utilization of the health system<sup><xref ref-type="bibr" rid="CR65">65</xref></sup>. We extracted data from the primary records, health facility laboratory registers and observed a reduction in confirmed malaria cases in health facilities within sprayed areas post-IRS, with no change in the number of cases detected at the facilities in control regions. Reduction in <italic>An. funestus</italic> densities, sporozoite rates and man-biting rates coupled with reduced malaria cases following one round IRS with pirimiphos-methyl provide compelling evidence of the effectiveness of IRS in malaria transmission reduction when implemented with an effective insecticide to which mosquito populations are susceptible.</p><p id="Par45">This study compared entomological indices and malaria test positivity rates between one district receiving IRS in addition to LLINs and the neighboring districts receiving only LLINs. This approach was driven by the fact that it was an effectiveness study rather than a highly controlled study to assess efficacy. However, the study does raise questions about the optimal approach to evaluate vector control interventions such as IRS or larval source management which target population-wide suppression of malaria vectors and where the scale of the study may affect the magnitude with which the results can be generalized. While cluster randomized trials might be ideal for the evaluation of these tools, the flight range of mosquitoes should be accounted for as migration from outside clusters and/or spill-over effects from intervention clusters into control clusters may bias results towards the null<sup><xref ref-type="bibr" rid="CR66">66</xref></sup>. However, larger clusters, in addition to being logistically unfeasible, may result in clusters that are spatially distant and less similar which may also bias the results. In this study, the non-random selection of areas for IRS, the low number of clusters and their spatial distance are a potential limitation of the study. However, the before and after design coupled with multiple endpoints indicating similar trends in mosquito densities and malaria-related outcomes strongly suggest that the results are robust.</p></sec><sec id="Sec21" sec-type="conclusion"><title>Conclusion</title><p id="Par46">IRS with pirimiphos-methyl was highly effective for the control of indoor biting and indoor resting, pyrethroid-resistant <italic>An. funestus</italic> and resulted in substantially reduced numbers of this primary vector species coupled with reduced malaria cases. Due to the long residual effect of pirimiphos-methyl, it was possible to achieve year-round protection with a single round of IRS. Sustaining these gains is a priority for the Kenya NMCP and development partners and IRS should continue to be implemented to sustain the impact on <italic>An. funestus</italic>. However, there was less of an impact of spraying on <italic>An. arabiensis</italic> populations, likely due to their exophilic nature. Additional control measures are needed to control outdoor biting and resting <italic>An. arabiensis</italic>.</p></sec><sec sec-type="supplementary-material"><title>Supplementary information</title><sec id="Sec22"><p>
<supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="41598_2020_61350_MOESM1_ESM.pdf"><caption><p>Supplementary Information.</p></caption></media></supplementary-material>
</p></sec></sec></body><back><fn-group><fn><p><bold>Publisher&#x02019;s note</bold> Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p></fn></fn-group><sec><title>Supplementary information</title><p>is available for this paper at 10.1038/s41598-020-61350-2.</p></sec><ack><title>Acknowledgements</title><p>We would like to thank the residents of Migori and Homa Bay Counties from whose houses we conducted mosquito sampling. We also thank the County Health Management Teams in Migori and Homa Bay counties and health facility management for Uriri, Rongo, Sori Karungu and Nytike, Ndhiwa and Marindi health facilities for access to the hospital laboratory health records. Kate Zinszer provided invaluable help with the ARIMA analyses. Lastly, we thank the PMI-VectorLink Kenya project staff, Alexia Mshambala for participation in health facility surveillance and entomology staff who conducted mosquito collection and sample analysis, especially, Benjamin Oloo, Tobias Odongo, Amos Wabwile, Celestine Wekesa, and Mercy Nduta. The study was supported by the United States President&#x02019;s Malaria Initiative through the United States Agency for International Development Africa Indoor Residual Spraying Project. The findings and conclusions in this report do not necessarily reflect the official position of the U.S. Department of Health and Human Services or the U.S. Centers for Disease Control and Prevention. The use of trade names is for identification only and does not imply endorsement by the U.S. Department of Health and Human Services or the U.S. Centers for Disease Control and Prevention.</p></ack><notes notes-type="author-contribution"><title>Author contributions</title><p>B.A., J.E.G., B.L. and R.M.O. wrote the protocol and designed the study. B.A., D.O., M.M. and R.M.O., participated in data collection and sample analysis. J.E.G., S.J.T., F.K., A.M.S., M.J.D., R.T.P., S.M., K.N. and R.M.O., provided technical oversight during the study. B.A., J.E.G., M.M., F.K., J.M. and M.J.D., analyzed the data and prepared manuscript figures. B.A., J.E.G. M.J.D. and R.M.O. wrote the manuscript. All authors reviewed and approved the final manuscript.</p></notes><notes notes-type="data-availability"><title>Data availability</title><p>All data generated or analyzed during this study are included in this published article and its Supplementary Information files.</p></notes><notes notes-type="COI-statement"><title>Competing interests</title><p id="Par48">The authors declare no competing interests.</p></notes><ref-list id="Bib1"><title>References</title><ref id="CR1"><label>1.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bhatt</surname><given-names>S</given-names></name><etal/></person-group><article-title>The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015</article-title><source>Nature</source><year>2015</year><volume>526</volume><fpage>207</fpage><lpage>211</lpage><pub-id pub-id-type="doi">10.1038/nature15535</pub-id><pub-id pub-id-type="pmid">26375008</pub-id></element-citation></ref><ref id="CR2"><label>2.</label><mixed-citation publication-type="other">Pluess, B., Tanser, F. C., Lengeler, C. &#x00026; Sharp, B. L. Indoor residual spraying for preventing malaria. <italic>Cochrane Database Syst Rev</italic>, CD006657, 10.1002/14651858.CD006657.pub2 (2010).</mixed-citation></ref><ref id="CR3"><label>3.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Protopopoff</surname><given-names>N</given-names></name><etal/></person-group><article-title>Combination of Insecticide Treated Nets and Indoor Residual Spraying in Northern Tanzania Provides Additional Reduction in Vector Population Density and Malaria Transmission Rates Compared to Insecticide Treated Nets Alone: A Randomised Control Trial</article-title><source>PloS one</source><year>2015</year><volume>10</volume><fpage>e0142671</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0142671</pub-id><pub-id pub-id-type="pmid">26569492</pub-id></element-citation></ref><ref id="CR4"><label>4.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Protopopoff</surname><given-names>N</given-names></name><etal/></person-group><article-title>Effectiveness of a long-lasting piperonyl butoxide-treated insecticidal net and indoor residual spray interventions, separately and together, against malaria transmitted by pyrethroid-resistant mosquitoes: a cluster, randomised controlled, two-by-two factorial design trial</article-title><source>The Lancet</source><year>2018</year><volume>391</volume><fpage>1577</fpage><lpage>1588</lpage><pub-id pub-id-type="doi">10.1016/s0140-6736(18)30427-6</pub-id></element-citation></ref><ref id="CR5"><label>5.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gimnig</surname><given-names>JE</given-names></name><etal/></person-group><article-title>Effect of permethrin-treated bed nets on the spatial distribution of malaria vectors in western Kenya</article-title><source>American Journal of Tropical Medicine and Hygiene</source><year>2003</year><volume>68</volume><fpage>115</fpage><lpage>120</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2003.68.115</pub-id></element-citation></ref><ref id="CR6"><label>6.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gimnig</surname><given-names>JE</given-names></name><etal/></person-group><article-title>Impact of permethrin-treated bed nets on entomologic indices in an area of intense year-round malaria transmission</article-title><source>American Journal of Tropical Medicine and Hygiene</source><year>2003</year><volume>68</volume><fpage>16</fpage><lpage>22</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2003.68.16</pub-id><pub-id pub-id-type="pmid">12749481</pub-id></element-citation></ref><ref id="CR7"><label>7.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mbogo</surname><given-names>CN</given-names></name><name><surname>Baya</surname><given-names>NM</given-names></name><name><surname>Ofulla</surname><given-names>AV</given-names></name><name><surname>Githure</surname><given-names>JI</given-names></name><name><surname>Snow</surname><given-names>RW</given-names></name></person-group><article-title>The impact of permethrin-impregnated bednets on malaria vectors of the Kenyan coast</article-title><source>Medical Veterinary Entomology</source><year>1996</year><volume>10</volume><fpage>251</fpage><lpage>259</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2915.1996.tb00739.x</pub-id><pub-id pub-id-type="pmid">8887336</pub-id></element-citation></ref><ref id="CR8"><label>8.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bayoh</surname><given-names>MN</given-names></name><etal/></person-group><article-title>Anopheles gambiae: historical population decline associated with regional distribution of insecticide-treated bed nets in western Nyanza Province, Kenya</article-title><source>Malaria journal</source><year>2010</year><volume>9</volume><fpage>62</fpage><pub-id pub-id-type="doi">10.1186/1475-2875-9-62</pub-id><pub-id pub-id-type="pmid">20187956</pub-id></element-citation></ref><ref id="CR9"><label>9.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sougoufara</surname><given-names>S</given-names></name><etal/></person-group><article-title>Biting by Anopheles funestus in broad daylight after use of long-lasting insecticidal nets: a new challenge to malaria elimination</article-title><source>Malaria journal</source><year>2014</year><volume>13</volume><fpage>125</fpage><pub-id pub-id-type="doi">10.1186/1475-2875-13-125</pub-id><pub-id pub-id-type="pmid">24678587</pub-id></element-citation></ref><ref id="CR10"><label>10.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ototo</surname><given-names>EN</given-names></name><etal/></person-group><article-title>Surveillance of malaria vector population density and biting behaviour in western Kenya</article-title><source>Malaria journal</source><year>2015</year><volume>14</volume><fpage>244</fpage><pub-id pub-id-type="doi">10.1186/s12936-015-0763-7</pub-id><pub-id pub-id-type="pmid">26082138</pub-id></element-citation></ref><ref id="CR11"><label>11.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Takken</surname><given-names>W</given-names></name></person-group><article-title>Do insecticide-treated bednets have an effect on malaria vectors?</article-title><source>Tropical Medicine &#x00026; International Health</source><year>2002</year><volume>7</volume><fpage>1022</fpage><lpage>1030</lpage><pub-id pub-id-type="doi">10.1046/j.1365-3156.2002.00983.x</pub-id><pub-id pub-id-type="pmid">12460393</pub-id></element-citation></ref><ref id="CR12"><label>12.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ndenga</surname><given-names>BA</given-names></name><etal/></person-group><article-title>Malaria vectors and their blood-meal sources in an area of high bed net ownership in the western Kenya highlands</article-title><source>Malaria journal</source><year>2016</year><volume>15</volume><fpage>76</fpage><pub-id pub-id-type="doi">10.1186/s12936-016-1115-y</pub-id><pub-id pub-id-type="pmid">26857915</pub-id></element-citation></ref><ref id="CR13"><label>13.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mathenge</surname><given-names>EM</given-names></name><etal/></person-group><article-title>Effect of permethrin-impregnated nets on exiting behavior, blood feeding success, and time of feeding of malaria mosquitoes (Diptera: Culicidae) in western Kenya</article-title><source>J Med Entomol</source><year>2001</year><volume>38</volume><fpage>531</fpage><lpage>536</lpage><pub-id pub-id-type="doi">10.1603/0022-2585-38.4.531</pub-id><pub-id pub-id-type="pmid">11476333</pub-id></element-citation></ref><ref id="CR14"><label>14.</label><mixed-citation publication-type="other">KMIS. Kenya Malaria Indicator Survey (National Malaria Control Programe, Kenya Ministry of Heath, Nairobi, Kenya, 2015).</mixed-citation></ref><ref id="CR15"><label>15.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Noor</surname><given-names>AM</given-names></name><etal/></person-group><article-title>The risks of malaria infection in Kenya in 2009</article-title><source>BMC Infectious Diseases</source><year>2009</year><volume>9</volume><fpage>180</fpage><pub-id pub-id-type="doi">10.1186/1471-2334-9-180</pub-id><pub-id pub-id-type="pmid">19930552</pub-id></element-citation></ref><ref id="CR16"><label>16.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wesolowski</surname><given-names>A</given-names></name><etal/></person-group><article-title>Quantifying the impact of human mobility on malaria</article-title><source>Science</source><year>2012</year><volume>338</volume><fpage>267</fpage><lpage>270</lpage><pub-id pub-id-type="doi">10.1126/science.1223467</pub-id><pub-id pub-id-type="pmid">23066082</pub-id></element-citation></ref><ref id="CR17"><label>17.</label><mixed-citation publication-type="other">Nyunt, M. H. <italic>et al</italic>. Challenges in universal coverage and utilization of insecticide-treated bed nets in migrant plantation workers in Myanmar. <italic>Malaria journal</italic><bold>13</bold> (2014).</mixed-citation></ref><ref id="CR18"><label>18.</label><mixed-citation publication-type="other">Atieli, H. E. <italic>et al</italic>. Insecticide-treated net (ITN) ownership, usage, and malaria transmission in the highlands of western Kenya. <italic>Parasites &#x00026; vectors</italic><bold>4</bold> (2011).</mixed-citation></ref><ref id="CR19"><label>19.</label><mixed-citation publication-type="other">WHO. World Malaria Report. (WHO, Geneva, 2018).</mixed-citation></ref><ref id="CR20"><label>20.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ochomo</surname><given-names>EO</given-names></name><etal/></person-group><article-title>The efficacy of long-lasting nets with declining physical integrity may be compromised in areas with high levels of pyrethroid resistance</article-title><source>Malaria journal</source><year>2013</year><volume>12</volume><fpage>368</fpage><pub-id pub-id-type="doi">10.1186/1475-2875-12-368</pub-id><pub-id pub-id-type="pmid">24156715</pub-id></element-citation></ref><ref id="CR21"><label>21.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ochomo</surname><given-names>E</given-names></name><etal/></person-group><article-title>Presence of the knockdown resistance mutation, Vgsc-1014F in Anopheles gambiae and An. arabiensis in western Kenya</article-title><source>Parasites &#x00026; vectors</source><year>2015</year><volume>8</volume><fpage>616</fpage><pub-id pub-id-type="doi">10.1186/s13071-015-1223-5</pub-id><pub-id pub-id-type="pmid">26626424</pub-id></element-citation></ref><ref id="CR22"><label>22.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ochomo</surname><given-names>E</given-names></name><etal/></person-group><article-title>Pyrethroid resistance in Anopheles gambiae s.s. and Anopheles arabiensis in western Kenya: phenotypic, metabolic and target site characterizations of three populations</article-title><source>Medical and veterinary entomology</source><year>2013</year><volume>27</volume><fpage>156</fpage><lpage>164</lpage><pub-id pub-id-type="doi">10.1111/j.1365-2915.2012.01039.x</pub-id><pub-id pub-id-type="pmid">22861380</pub-id></element-citation></ref><ref id="CR23"><label>23.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ochomo</surname><given-names>E</given-names></name><etal/></person-group><article-title>Pyrethroid susceptibility of malaria vectors in four Districts of western Kenya</article-title><source>Parasit Vectors.</source><year>2014</year><volume>4</volume><fpage>310</fpage><pub-id pub-id-type="doi">10.1186/1756-3305-7-310</pub-id></element-citation></ref><ref id="CR24"><label>24.</label><mixed-citation publication-type="other">KMIS. Kenya Malaria Indicator Survey. (Division of Malaria Control, Ministry of Public Health and Sanitation, Nairobi, Kenya, 2010).</mixed-citation></ref><ref id="CR25"><label>25.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hamel</surname><given-names>MJ</given-names></name><etal/></person-group><article-title>The combination of indoor residual spraying and insecticide-treated nets provides added protection against malaria compared with insecticide-treated nets alone</article-title><source>American Journal of Tropical Medicine and Hygiene</source><year>2011</year><volume>85</volume><fpage>1080</fpage><lpage>1086</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2011.10-0684</pub-id><pub-id pub-id-type="pmid">22144448</pub-id></element-citation></ref><ref id="CR26"><label>26.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gimnig</surname><given-names>JE</given-names></name><etal/></person-group><article-title>The Effect of Indoor Residual Spraying on the Prevalence of Malaria Parasite Infection, Clinical Malaria and Anemia in an Area of Perennial Transmission and Moderate Coverage of Insecticide Treated Nets in Western Kenya</article-title><source>PloS one</source><year>2016</year><volume>11</volume><fpage>e0145282</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0145282</pub-id><pub-id pub-id-type="pmid">26731524</pub-id></element-citation></ref><ref id="CR27"><label>27.</label><mixed-citation publication-type="other">NMCP. (ed Ministry of Health) (Nairobi, Kenya, 2015).</mixed-citation></ref><ref id="CR28"><label>28.</label><mixed-citation publication-type="other">WHO. Global Plan For Insecticide Resistance Management In Malaria Vectors. (World Health Organization, Geneva, Switzerland, 2012).</mixed-citation></ref><ref id="CR29"><label>29.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rowland</surname><given-names>Mark</given-names></name><name><surname>Boko</surname><given-names>Pelagie</given-names></name><name><surname>Odjo</surname><given-names>Abibatou</given-names></name><name><surname>Asidi</surname><given-names>Alex</given-names></name><name><surname>Akogbeto</surname><given-names>Martin</given-names></name><name><surname>N&#x02019;Guessan</surname><given-names>Raphael</given-names></name></person-group><article-title>A New Long-Lasting Indoor Residual Formulation of the Organophosphate Insecticide Pirimiphos Methyl for Prolonged Control of Pyrethroid-Resistant Mosquitoes: An Experimental Hut Trial in Benin</article-title><source>PLoS ONE</source><year>2013</year><volume>8</volume><issue>7</issue><fpage>e69516</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0069516</pub-id><pub-id pub-id-type="pmid">23936033</pub-id></element-citation></ref><ref id="CR30"><label>30.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mashauri</surname><given-names>FM</given-names></name><etal/></person-group><article-title>Indoor residual spraying with micro-encapsulated pirimiphos-methyl (Actellic(R) 300CS) against malaria vectors in the Lake Victoria basin, Tanzania</article-title><source>PloS one</source><year>2017</year><volume>12</volume><fpage>e0176982</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0176982</pub-id><pub-id pub-id-type="pmid">28489935</pub-id></element-citation></ref><ref id="CR31"><label>31.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chanda</surname><given-names>E</given-names></name><etal/></person-group><article-title>Preventing malaria transmission by indoor residual spraying in Malawi: grappling with the challenge of uncertain sustainability</article-title><source>Malaria journal</source><year>2015</year><volume>14</volume><fpage>254</fpage><pub-id pub-id-type="doi">10.1186/s12936-015-0759-3</pub-id><pub-id pub-id-type="pmid">26104657</pub-id></element-citation></ref><ref id="CR32"><label>32.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Haji</surname><given-names>KA</given-names></name><etal/></person-group><article-title>Efficacy, persistence and vector susceptibility to pirimiphos-methyl (Actellic 300CS) insecticide for indoor residual spraying in Zanzibar</article-title><source>Parasites &#x00026; vectors</source><year>2015</year><volume>8</volume><fpage>628</fpage><pub-id pub-id-type="doi">10.1186/s13071-015-1239-x</pub-id><pub-id pub-id-type="pmid">26652708</pub-id></element-citation></ref><ref id="CR33"><label>33.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Okara</surname><given-names>RM</given-names></name><etal/></person-group><article-title>Distribution of the main malaria vectors in Kenya</article-title><source>Malaria journal</source><year>2010</year><volume>9</volume><fpage>69</fpage><pub-id pub-id-type="doi">10.1186/1475-2875-9-69</pub-id><pub-id pub-id-type="pmid">20202199</pub-id></element-citation></ref><ref id="CR34"><label>34.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McCann</surname><given-names>RS</given-names></name><etal/></person-group><article-title>Reemergence of Anopheles funestus as a vector of Plasmodium falciparum in western Kenya after long-term implementation of insecticide-treated bed nets</article-title><source>American Journal of Tropical Medicine and Hygiene</source><year>2014</year><volume>90</volume><fpage>597</fpage><lpage>604</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.13-0614</pub-id><pub-id pub-id-type="pmid">24470562</pub-id></element-citation></ref><ref id="CR35"><label>35.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ondeto</surname><given-names>BM</given-names></name><etal/></person-group><article-title>Current status of insecticide resistance among malaria vectors in Kenya</article-title><source>Parasites &#x00026; vectors</source><year>2017</year><volume>10</volume><fpage>429</fpage><pub-id pub-id-type="doi">10.1186/s13071-017-2361-8</pub-id><pub-id pub-id-type="pmid">28927428</pub-id></element-citation></ref><ref id="CR36"><label>36.</label><mixed-citation publication-type="other">AIRS. Kenya End of Spray Report. (Abt Associates Inc, Bethesda, MD, May 2017).</mixed-citation></ref><ref id="CR37"><label>37.</label><mixed-citation publication-type="other">Seyoum, A. <italic>et al</italic>. Human exposure to anopheline mosquitoes occurs primarily indoors, even for users of insecticide-treated nets in Luangwa Valley, South-east Zambia. <italic>Parasites &#x00026; vectors</italic><bold>5</bold> (2012).</mixed-citation></ref><ref id="CR38"><label>38.</label><mixed-citation publication-type="other">WHO. Guidelines for testing mosquito adulticides for indoor residual spraying and treatment of mosquito nets (Worl Health Organization, Geneva, Switzerland, 2006).</mixed-citation></ref><ref id="CR39"><label>39.</label><mixed-citation publication-type="other">Gillies, M. T. &#x00026; Coetzee, M. A supplement to the Anophelinae of Africa South of the Sahara (Afrotropical region). (Johannesburg: South Africa Medical Resaerch Institute, 1987).</mixed-citation></ref><ref id="CR40"><label>40.</label><mixed-citation publication-type="other">Gillies, M. T. &#x00026; De Meillon, B. The Anophelinae of Africa, South of the Sahara (Ethopian zoogeographical region). 1968:1954 (Johannesburg: South African Institute for Madical Research; 1968).</mixed-citation></ref><ref id="CR41"><label>41.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wirtz</surname><given-names>RA</given-names></name><name><surname>Burkot</surname><given-names>TR</given-names></name><name><surname>Graves</surname><given-names>PM</given-names></name><name><surname>Andre</surname><given-names>RG</given-names></name></person-group><article-title>Field evaluation of enzyme-linked immunosorbent assays for Plasmodiumfalciparum and Plasmodium vivax sporozoites in mosquitoes (Diptera: Culicidae) from Papua New Guinea</article-title><source>Journal of Medical Entomology</source><year>1987</year><volume>24</volume><fpage>433</fpage><lpage>437</lpage><pub-id pub-id-type="doi">10.1093/jmedent/24.4.433</pub-id><pub-id pub-id-type="pmid">3305949</pub-id></element-citation></ref><ref id="CR42"><label>42.</label><mixed-citation publication-type="other">MR4. <italic>Methods in Anopheles Research</italic>. Vol. Chapter 5 267&#x02013;269 (Mark Q. Benedict CDC, 2015).</mixed-citation></ref><ref id="CR43"><label>43.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kent</surname><given-names>RJ</given-names></name><name><surname>Norris</surname><given-names>DE</given-names></name></person-group><article-title>Identification of mammalian blood meals in mosquitoes by a multiplexed polymerase chain reaction targeting cytochrome B</article-title><source>The American journal of tropical medicine and hygiene</source><year>2005</year><volume>73</volume><fpage>336</fpage><lpage>342</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2005.73.336</pub-id><pub-id pub-id-type="pmid">16103600</pub-id></element-citation></ref><ref id="CR44"><label>44.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Scott</surname><given-names>JA</given-names></name><name><surname>Brogdon</surname><given-names>WG</given-names></name><name><surname>Collins</surname><given-names>FH</given-names></name></person-group><article-title>Identification of single specimens of the <italic>Anopheles gambiae</italic> complex by the polymerase chain reaction</article-title><source>American Journal of Tropical Medicine and Hygiene</source><year>1993</year><volume>49</volume><fpage>520</fpage><lpage>529</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.1993.49.520</pub-id><pub-id pub-id-type="pmid">8214283</pub-id></element-citation></ref><ref id="CR45"><label>45.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Altman</surname><given-names>DG</given-names></name><name><surname>Bland</surname><given-names>JM</given-names></name></person-group><article-title>Interaction revisited: the difference between two estimates</article-title><source>BMJ</source><year>2003</year><volume>326</volume><fpage>219</fpage><pub-id pub-id-type="doi">10.1136/bmj.326.7382.219</pub-id><pub-id pub-id-type="pmid">12543843</pub-id></element-citation></ref><ref id="CR46"><label>46.</label><mixed-citation publication-type="other">Team, R. C. (ed. R Foundation for Statistical Computing) (Vienna, Austria, 2018).</mixed-citation></ref><ref id="CR47"><label>47.</label><mixed-citation publication-type="other">Chan, K.-S. &#x00026; Ripley, B. TSA: Time Series Analysis. <italic>R package version 1.2</italic> (2018).</mixed-citation></ref><ref id="CR48"><label>48.</label><mixed-citation publication-type="other">Qiu, D. aTSA: Alternative Time Series Analysis. <italic>R package version 3.1.2</italic> (2015).</mixed-citation></ref><ref id="CR49"><label>49.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smith</surname><given-names>A</given-names></name><name><surname>Draper</surname><given-names>CC</given-names></name></person-group><article-title>Malaria in the Taveta area of Kenya and Tanganyika. Part I. Epidemiology</article-title><source>East Africa Medical Journal</source><year>1959</year><volume>36</volume><fpage>99</fpage></element-citation></ref><ref id="CR50"><label>50.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bruce-Chwatt</surname><given-names>LJ</given-names></name><name><surname>Gockel</surname><given-names>CW</given-names></name></person-group><article-title>A study of blood feeding patterns of <italic>Anopheles</italic> mosquitoes through precipitin test</article-title><source>Bull. Wld. Hlth. Org</source><year>1960</year><volume>22</volume><fpage>685</fpage></element-citation></ref><ref id="CR51"><label>51.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Githeko</surname><given-names>AK</given-names></name><name><surname>Service</surname><given-names>MW</given-names></name><name><surname>Mbogo</surname><given-names>CM</given-names></name><name><surname>Atieli</surname><given-names>FK</given-names></name><name><surname>Juma</surname><given-names>FO</given-names></name></person-group><article-title>Origin of blood meals in indoor and outdoor resting malaria vectors in western Kenya</article-title><source>Acta tropica</source><year>1994</year><volume>58</volume><fpage>307</fpage><lpage>316</lpage><pub-id pub-id-type="doi">10.1016/0001-706X(94)90024-8</pub-id><pub-id pub-id-type="pmid">7709869</pub-id></element-citation></ref><ref id="CR52"><label>52.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bayoh</surname><given-names>MN</given-names></name><etal/></person-group><article-title>Persistently high estimates of late night, indoor exposure to malaria vectors despite high coverage of insecticide treated nets</article-title><source>Parasites &#x00026; vectors</source><year>2014</year><volume>7</volume><fpage>380</fpage><pub-id pub-id-type="doi">10.1186/1756-3305-7-380</pub-id><pub-id pub-id-type="pmid">25141761</pub-id></element-citation></ref><ref id="CR53"><label>53.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Githeko</surname><given-names>AK</given-names></name><etal/></person-group><article-title>Some observations on the biting behavior of Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus and their implications for malaria control</article-title><source>Exp Parasitol</source><year>1996</year><volume>82</volume><fpage>306</fpage><lpage>315</lpage><pub-id pub-id-type="doi">10.1006/expr.1996.0038</pub-id><pub-id pub-id-type="pmid">8631382</pub-id></element-citation></ref><ref id="CR54"><label>54.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Killeen</surname><given-names>GF</given-names></name><name><surname>Govella</surname><given-names>NJ</given-names></name><name><surname>Lwetoijera</surname><given-names>DW</given-names></name><name><surname>Okumu</surname><given-names>FO</given-names></name></person-group><article-title>Most outdoor malaria transmission by behaviourally-resistant Anopheles arabiensis is mediated by mosquitoes that have previously been inside houses</article-title><source>Malaria journal</source><year>2016</year><volume>15</volume><fpage>225</fpage><pub-id pub-id-type="doi">10.1186/s12936-016-1280-z</pub-id><pub-id pub-id-type="pmid">27093890</pub-id></element-citation></ref><ref id="CR55"><label>55.</label><mixed-citation publication-type="other">Abong&#x02019;o, B. <italic>et al</italic>. Host Decoy Trap (HDT) with cattle odour is highly effective for collection of exophagic malaria vectors. <italic>Parasites &#x00026; vectors</italic><bold>11</bold>, 10.1186/s13071-018-3099-7 (2018).</mixed-citation></ref><ref id="CR56"><label>56.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chandler</surname><given-names>JA</given-names></name><name><surname>Highton</surname><given-names>RB</given-names></name><name><surname>Hill</surname><given-names>MN</given-names></name></person-group><article-title>Mosquitoes of the Kano Plain, Kenya. I. Results of indoor collections in irrigated and nonirrigated areas using human bait and light traps</article-title><source>J Med Entomol</source><year>1975</year><volume>12</volume><fpage>504</fpage><lpage>510</lpage><pub-id pub-id-type="doi">10.1093/jmedent/12.5.504</pub-id><pub-id pub-id-type="pmid">4620</pub-id></element-citation></ref><ref id="CR57"><label>57.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Moiroux</surname><given-names>N</given-names></name><etal/></person-group><article-title>Human exposure to early morning Anopheles funestus biting behavior and personal protection provided by long-lasting insecticidal nets</article-title><source>PloS one</source><year>2014</year><volume>9</volume><fpage>e104967</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0104967</pub-id><pub-id pub-id-type="pmid">25115830</pub-id></element-citation></ref><ref id="CR58"><label>58.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lindsay</surname><given-names>SW</given-names></name><name><surname>Snow</surname><given-names>RW</given-names></name></person-group><article-title>The trouble with eaves; house entry by vectors of malaria</article-title><source>Trans R Soc Trop Med Hyg</source><year>1988</year><volume>82</volume><fpage>645</fpage><lpage>646</lpage><pub-id pub-id-type="doi">10.1016/0035-9203(88)90546-9</pub-id><pub-id pub-id-type="pmid">3256125</pub-id></element-citation></ref><ref id="CR59"><label>59.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lindsay</surname><given-names>SW</given-names></name><name><surname>Emerson</surname><given-names>MP</given-names></name><name><surname>D.</surname><given-names>CJ</given-names></name></person-group><article-title>Reducing malari a by mosquito-proofing houses</article-title><source>Trends in Parasitology</source><year>2002</year><volume>18</volume><fpage>510</fpage><lpage>514</lpage><pub-id pub-id-type="doi">10.1016/S1471-4922(02)02382-6</pub-id><pub-id pub-id-type="pmid">12473368</pub-id></element-citation></ref><ref id="CR60"><label>60.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Oduro</surname><given-names>AR</given-names></name><etal/></person-group><article-title>Monitoring malaria using health facility based surveys: challenges and limitations</article-title><source>BMC Public Health</source><year>2016</year><volume>16</volume><fpage>354</fpage><pub-id pub-id-type="doi">10.1186/s12889-016-2858-7</pub-id><pub-id pub-id-type="pmid">27102913</pub-id></element-citation></ref><ref id="CR61"><label>61.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Githinji</surname><given-names>S</given-names></name><etal/></person-group><article-title>A national health facility survey of malaria infection among febrile patients in Kenya, 2014</article-title><source>Malaria journal</source><year>2016</year><volume>15</volume><fpage>591</fpage><pub-id pub-id-type="doi">10.1186/s12936-016-1638-2</pub-id><pub-id pub-id-type="pmid">27931229</pub-id></element-citation></ref><ref id="CR62"><label>62.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Macedo de Oliveira</surname><given-names>A</given-names></name><etal/></person-group><article-title>Prevalence of malaria among patients attending public health facilities in Maputo City, Mozambique</article-title><source>The American journal of tropical medicine and hygiene</source><year>2011</year><volume>85</volume><fpage>1002</fpage><lpage>1007</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2011.11-0365</pub-id><pub-id pub-id-type="pmid">22144434</pub-id></element-citation></ref><ref id="CR63"><label>63.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ashton</surname><given-names>Ruth A.</given-names></name><name><surname>Bennett</surname><given-names>Adam</given-names></name><name><surname>Al-Mafazy</surname><given-names>Abdul-Wahid</given-names></name><name><surname>Abass</surname><given-names>Ali K.</given-names></name><name><surname>Msellem</surname><given-names>Mwinyi I.</given-names></name><name><surname>McElroy</surname><given-names>Peter</given-names></name><name><surname>Kachur</surname><given-names>S. Patrick</given-names></name><name><surname>Ali</surname><given-names>Abdullah S.</given-names></name><name><surname>Yukich</surname><given-names>Joshua</given-names></name><name><surname>Eisele</surname><given-names>Thomas P.</given-names></name><name><surname>Bhattarai</surname><given-names>Achuyt</given-names></name></person-group><article-title>Use of Routine Health Information System Data to Evaluate Impact of Malaria Control Interventions in Zanzibar, Tanzania from 2000 to 2015</article-title><source>EClinicalMedicine</source><year>2019</year><volume>12</volume><fpage>11</fpage><lpage>19</lpage><pub-id pub-id-type="doi">10.1016/j.eclinm.2019.05.011</pub-id><pub-id pub-id-type="pmid">31388659</pub-id></element-citation></ref><ref id="CR64"><label>64.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>D&#x02019;Acremont</surname><given-names>V</given-names></name><name><surname>Lengeler</surname><given-names>C</given-names></name><name><surname>Genton</surname><given-names>B</given-names></name></person-group><article-title>Reduction in the proportion of fevers associated with Plasmodium falciparum parasitaemia in Africa: a systematic review</article-title><source>Malaria journal</source><year>2010</year><volume>9</volume><fpage>240</fpage><pub-id pub-id-type="doi">10.1186/1475-2875-9-240</pub-id><pub-id pub-id-type="pmid">20727214</pub-id></element-citation></ref><ref id="CR65"><label>65.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cibulskis</surname><given-names>RE</given-names></name><etal/></person-group><article-title>Estimating trends in the burden of malaria at country level</article-title><source>The American journal of tropical medicine and hygiene</source><year>2007</year><volume>77</volume><fpage>133</fpage><lpage>137</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2007.77.133</pub-id></element-citation></ref><ref id="CR66"><label>66.</label><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hawley</surname><given-names>WA</given-names></name><etal/></person-group><article-title>Community-wide effects of permethrin-treated bed nets on child mortality and malaria morbidity in western Kenya</article-title><source>The American journal of tropical medicine and hygiene</source><year>2003</year><volume>68</volume><fpage>121</fpage><lpage>127</lpage><pub-id pub-id-type="doi">10.4269/ajtmh.2003.68.121</pub-id><pub-id pub-id-type="pmid">12749495</pub-id></element-citation></ref></ref-list></back></article>