Changes in malaria burden and transmission in sentinel sites after the roll-out of long-lasting insecticidal nets in Papua New Guinea

Background Papua New Guinea exhibits a complex malaria epidemiology due to diversity in malaria parasites, mosquito vectors, human hosts, and their natural environment. Heterogeneities in transmission and burden of malaria at various scales are likely to affect the success of malaria control interventions, and vice-versa. This manuscript assesses changes in malaria prevalence, incidence and transmission in sentinel sites following the first national distribution of long-lasting insecticidal nets (LLINs). Methods Before and after the distribution of LLINs, data collection in six purposively selected sentinel sites included clinical surveillance in the local health facility, household surveys and entomological surveys. Not all activities were carried out in all sites. Mosquitoes were collected by human landing catches. Diagnosis of malaria infection in humans was done by rapid diagnostic test, light microscopy and PCR for species confirmation. Results Following the roll-out of LLINs, the average monthly malaria incidence rate dropped from 13/1,000 population to 2/1,000 (incidence rate ratio = 0.12; 95 % CI: 0.09–0.17; P < 0.001). The average population prevalence of malaria decreased from 15.7 % pre-LLIN to 4.8 % post-LLIN (adjusted odds ratio = 0.26; 95 % CI: 0.20–0.33; P < 0.001). In general, reductions in incidence and prevalence were more pronounced in infections with P. falciparum than with P. vivax. Additional morbidity indicators (anaemia, splenomegaly, self-reported fever) showed a decreasing trend in most sites. Mean Anopheles man biting rates decreased from 83 bites/person/night pre-LLIN to 31 post-LLIN (P = 0.008). Anopheles species composition differed between sites but everywhere diversity was lower post-LLIN. In two sites, post-LLIN P. vivax infections in anophelines had decreased but P. falciparum infections had increased despite the opposite observation in humans. Conclusions LLIN distribution had distinct effects on P. falciparum and P. vivax. Higher resilience of P. vivax may be attributed to relapses from hypnozoites and other biological characteristics favouring the transmission of P. vivax. The effect on vector species composition varied by location which is likely to impact on the effectiveness of LLINs. In-depth and longer-term epidemiological and entomological investigations are required to understand when and where residual transmission occurs and whether observed changes are sustained. Electronic supplementary material The online version of this article (doi:10.1186/s13071-016-1635-x) contains supplementary material, which is available to authorized users.


Background
Papua New Guinea (PNG) exhibits a diverse and complex malaria epidemiology [1]. Four Plasmodium species (P. falciparum, P. vivax, P. ovale and P. malariae) are endemic and more than ten Anopheles species filling different ecological niches have been incriminated in the parasite's transmission [2]. Historical malaria control measures in the 1960s and 1970s included indoor residual spraying of insecticides, mass drug administration as well as environmental management in certain areas resulting in substantial initial reductions in malaria in many locations [3,4]. The cessation of the spraying program in the 1980s coinciding with the decentralization of responsibilities for malaria control as well as emerging resistance of the malaria parasites to commonly used drugs [5,6] led to a subsequent resurgence, particularly of P. falciparum malaria, across most parts of PNG [7].
In 2004, the Government of Papua New Guinea reintensified its malaria control efforts with the financial support of a round 3 grant from the Global Fund to Fight AIDS, Tuberculosis and Malaria. The central component of this program was the first country-wide free distribution of long-lasting insecticidal mosquito nets (LLIN) [8]. Insecticide-treated nets had been shown to reduce the incidence and prevalence of falciparummalaria in children in PNG as early as 1985 [9]. Yet, while these results were later reflected in policy documents [10], little effort was made in practice to scale up mosquito net use. Coverage with mosquito nets, particularly insecticide-treated nets, therefore remained patchy and low in most parts of the country until 2004 [11][12][13]. The Global Fund supported program that subsequently facilitated a single round of LLIN distribution resulted in an increase in ownership and use of bednets, particularly LLINs. A national household survey after the distribution found 80 % of all households owning a mosquito net of any type and 65 % owning a LLIN, while 33 % of people reported to sleep under a LLIN [8]. In six sentinel surveillance sites established by the Papua New Guinea Institute of Medical Research (PNGIMR), average LLIN ownership increased from 9 % before the distribution to 89 % thereafter and LLIN use rose from 6 to 55 % [8].
This manuscript assesses malaria morbidity and transmission indicators in sentinel sites before and after the LLIN distribution. The study was part of the evaluation of the Global Fund round 3 malaria grant.

Study sites
Sentinel site locations were selected purposively in 2008 in places which were yet to be covered with the LLIN campaign. The sites were located in the Momase and Highlands regions of the main island of PNG ( Fig. 1) and consisted of three to four randomly selected villages in the catchment area of a sentinel health facility. Five sites were located in Momase: Dreikikir in the hills on the edges of the Sepik River basin (altitude 200-400 m), Finschhafen on the Morobe coastline of the Huon peninsula (altitude 0-50 m), Sausi in the Ramu River valley

Data collection
Data were collected prior to the LLIN campaign (10/ 2008-08/2009) and one year later after LLINs had been distributed (10/2009-08/2010). Data collection included clinical surveillance in the local health centre in three sites (Dreikikir, Mumeng and Sausi) over a period of approximately two months pre-and post-LLIN, coinciding roughly with the main malaria transmission season. A household survey in two to three randomly selected villages in the health centre's catchment area was carried out in five sites (Finschhafen, Mumeng, Sausi, Tabibuga and Yapsie) and entomological surveys were conducted in five sites (Dreikikir, Finschhafen, Mumeng, Sausi and Yapsie). For operational reasons not all activities could be carried out in all sites (Additional file 1: Table S1).

Clinical surveillance
For the clinical surveillance, all patients attending the health facility were screened for a history of fever within the previous three days. A research nurse collected a finger-prick blood sample from all consenting fever patients and malaria was diagnosed by rapid diagnostic test (RDT; ICT Malaria Combo, ICT Diagnostics, South Africa). From the same blood sample, thick and thin blood films were prepared on one slide and haemoglobin (Hb) levels were measured using a portable analyser (Hemo-Cue Hb 201+, HemoCue AB, Sweden). Anaemia and severe anaemia were defined in accordance with World Health Organization definitions applying age groupspecific cut-offs and altitude adjustment [14]. A dry blood spot was prepared on Whatman 3MM filter paper (GE Healthcare), whenever possible, and stored in individual plastic zip-bags with desiccant silica gel. Axillary temperature was measured with an electronic thermometer, and in patients between 2 and 9 years of age the spleen was palpated and graded according to Hackett [15]. Demographic details of the patient were recorded in a paper form. The subsequent clinical assessment, final diagnosis and treatment were then taken over by the health facility's clinician following standard procedures and results were recorded on the patient's study form.

Household survey
For the household survey, three (in the case of Yapsie four) villages were randomly selected from the sentinel health facility's catchment area. In each village, 30 to 35 households were randomly sampled from a list of all households. A structured questionnaire about the coverage with malaria control interventions was administered to the heads of sampled households. A finger-prick blood sample was collected from household members above five months of age for preparing a microscopy slide, a dry blood spot on filter paper, and for measuring Hb levels as described above. Symptomatic individuals were diagnosed on the spot using a malaria RDT and positive cases were treated according to standard treatment guidelines [16]. Axillary temperature was measured with an electronic thermometer. Details of the household survey methodology have been described in more details elsewhere [8,17].

Entomological surveys
Entomological surveys were conducted in two villages per site. Mosquitoes were collected by outdoor human landing catch from 6 households per village prior to the LLIN distribution and again 12 months later. The number of person-nights per collection ranged from 16 to 48. Collectors worked in pairs with one member collecting from 18:00 to 24:00 h and the second collecting from 24:00 to 06:00 h. Mosquitoes were stored per collection hour and identified to morphological species before storage on silica gel. Lysates from whole mosquitoes were screened for P. falciparum, P. vivax 210 and P. vivax 247 circumsporozoite proteins by enzyme-linked immunosorbent assay [18]. DNA was extracted from a portion of the lysate using DNEasy blood and tissue kit (QIAGEN, Maryland, USA). Mosquitoes that were morphologically identified as members of the An. punctulatus group were confirmed to species by polymerase chain reaction (PCR) restriction length polymorphism of the ITS2 region [2].

Laboratory procedures
Thin blood smears were first fixed with methanol and thick and thin smears stained with Giemsa and read by light microscopy independently by two PNGIMR microscopists. Discordant reads were confirmed by a senior microscopist. The number of parasites was counted until reaching 200 white blood cells and a slide was declared negative only after reading a minimum of 200 thick film fields. Molecular diagnosis was used to complement missing second reads and to disambiguate discordant species read results (clinical surveillance n = 127; household survey n = 87). DNA was extracted from filter papers using the DNeasy Blood and Tissue Kit (Qiagen, Valencia, USA) or the Favorgen 96-Well Genomic DNA Kit (Favorgen Biotech Corp., Taiwan) following the manufacturer's protocols. The molecular assay was a semiquantitative post-PCR, ligase detection reaction/ fluorescent microsphere assay described in more detail elsewhere [19][20][21].

Data analysis
Data were double-entered into a Microsoft FoxPro or Microsoft Access database at PNGIMR Goroka and analysed with Stata (StataCorp, USA) software. Differences in binary and categorical variables between study sites were assessed by Chi-square test and logistic regression, differences in continuous variables by t-test and linear regression and a non-parametric test for median age. Multivariable logistic and linear regression models assessing changes in malaria prevalence, parasite density and morbidity indicators were adjusted for age group and study site, anaemia furthermore for sex. In the clinical surveillance samples, primary diagnosis of malaria was based on the result of the RDT as performed in routine clinical practice while species-results were based on light microscopic/PCR diagnosis, which comprised only a sub-set of patients in the case of Mumeng post-LLIN. For the calculation of incidence rates, a stable population denominator was used for both years in the absence of available village-level growth rates.
The parasite species composition in the febrile patient sample changed in Dreikikir (but not in the other sites) from a clear dominance of P. falciparum over P. vivax (67.9 vs 12.9 %, Fisher's exact test, P < 0.001) to equal proportions of the two species (7.9 vs 8.6 %, P = 0.88). Pre-LLIN, P. falciparum infection prevalence peaked in fever patients aged 5-9 years, and P. vivax in 1-4 year olds. Post-LLIN, the P. falciparum peak had shifted to later age groups, while for P. vivax infections a clear peak could not be identified due to the low number of positive patients. Changes in test-positivity and species composition by site are shown in Fig. 2.

Malaria prevalence and morbidity in the general population
Across five sentinel sites, the average prevalence of malaria by light microscopy in the general population of the health facility catchment areas decreased from 15.7 % pre-LLIN to 4.8 % post-LLIN (AOR = 0.26; 95 % CI: 0.20-0.33; P < 0.001). The largest reduction was observed in Yapsie, the site with the highest pre-LLIN prevalence of both P. falciparum (18.9 %) and P. vivax (11.8 %). Across all sites, the reductions were more pronounced in infections with P. falciparum than with P. vivax ( Table 2). The changes in P. vivax infections were only statistically significant in two sites (Finschhafen and Yapsie) ( Table 2). In three sites (Finschhafen, Tabibuga and Yapsie), there was a change from P. falciparum dominance to P. vivax dominance, which reached statistical significance (P = 0.005) only in Tabibuga (Fig. 4). There was no statistically significant difference in the observed reductions between age groups or sex (Mantel Haenszel test of homogeneity of odds ratios, all P > 0.1). Axillary temperature, self-reported fever and anaemia were assessed as additional morbidity indicators in the general population in four sites (Finschhafen, Mumeng, Sausi, Tabibuga). Acute fever, defined as axillary temperature of > 37.5°C, was measured in 2.  Table S4).  Prior to LLIN roll-out, the majority of infected mosquitoes in three sites (Finschhafen, Mumeng and Sausi) carried P. vivax sporozoites, while in two sites (Dreikikir and Yapsie) P. falciparum infections were dominant. Sporozoite prevalence was lowest in Sausi, where the MBR was highest. Post-LLIN, no infected mosquitoes were found in Dreikikir, Mumeng and Yapsie. In the remaining sites, pre-LLIN infection of mosquitoes and monthly entomological inoculation rate (EIR) were higher for P. vivax than for P. falciparum. Post-LLIN, P. vivax infections had decreased, but P. falciparum infections had increased in both sites (Table 3).

Discussion
The roll-out of LLINs in 2009 was followed by substantial reductions in key indicators of malaria morbidity and transmission in six sites in PNG within less than Fig. 6 a Mean nightly Anopheles man biting rate (MBR; bites/person/night) before (grey) and after (white) LLIN distribution. b Anopheles species composition before and after LLIN distribution. Error bars show 95 % confidence intervals. *P < 0.05; **P < 0.001 one year. Average LLIN use in these sites had increased from 6 to 55 %, and 89 % of households in the sites owned at least one LLIN after the distribution [8]. The LLINs were new and efficacious [22] and there was no indication of resistance to pyrethroids in PNG, even though data at that time was limited [23]. The prevalence of Plasmodium spp. infection in the general population dropped significantly over the same period mainly due to a 76 % reduction in infections with P. falciparum (AOR = 0.23, P < 0.001). The prevalence of blood-stage P. vivax parasitaemia decreased less prominently (-60 %) reaching statistical significance in two of five study sites (Finschhafen and Yapsie). A significant reduction was also observed in the incidence of clinical malaria cases in three sites, primarily as a result of an 81 % drop in the proportion of P. falciparum-infected febrile patients (AOR = 0.1, P < 0.001). The proportion of cases infected with P. vivax decreased less across all age groups (-46 %), and statistically significantly only in children under 5 years of age. In the Dreikikir site, no significant decrease in the proportion of P. vivax cases was detected in any age group. These findings confirm results from a cohort study in young children conducted near Dreikikir which found that LLIN use was associated with a reduction in both P. falciparum infections and clinical episodes, while it only reduced P. vivax infections but not clinical episodes [24].
A relative increase in P. vivax over P. falciparum following implementation of malaria control measures confirms previous observations both from PNG [7,9,20,25] and elsewhere [26,27]. The shift in species composition was observed in both asymptomatic community members and febrile patients although not to the same degree in all sites. The differential impact of interventions on the two parasites species is likely a result of the parasite's [28] and potentially the vector's biology [29].
The ability of P. vivax to relapse from long-lasting liver stages and the higher probability of gametocytaemia are most likely key factors. A recent study confirmed that relapses from P. vivax hypnozoites contribute 80 % of the burden of malaria infection and clinical episodes in PNG [30]. In tropical strains, such as the South-West Pacific Chesson strains, relapses occur rapidly and frequently [31]. Although the majority of hypnozoites may therefore activate within 12 months of an initial infection, hypnozoites can survive in the human liver for up to 3-5 years. As a consequence, P. vivax blood-stage infections detected post-LLIN roll-out are not only due to new infections from mosquito bites but also due to relapses from hypnozoites established both pre-and post-LLIN [32]. Gametocytaemia was previously found to be driven by asexual blood stage parasitaemia, with a 10fold increase in parasite density leading to a 1.8-fold and 3.3-fold increase in the odds of carrying P. falciparum and P. vivax gametocytes, respectively [33]. Reducing the chance of P. vivax transmission would therefore require a stronger reduction of P. vivax parasitaemia. Also, P. vivax gametocytes appear early in an infection, when asexual densities are still very low [34] and almost all P. vivax infections are thus thought to be gametocyte positive [35] hence increasing the probability of transmission. Previous studies in PNG further found evidence of earlier biting Anopheles mosquitoes being more likely infected with P. vivax than P. falciparum [29,36], thus increasing the relative probability of P. vivax transmission before people retire to bed. This is further aggravated by the shorter extrinsic incubation period for P. vivax. As LLINs reduce the biting density as well as the lifespan of mosquitoes, an equal reduction in lifespan is more likely to interrupt P. falciparum transmission [37]. As a consequence of the above, particularly indoor vector control interventions are likely to have less impact on P. vivax than P. falciparum gametocytaemia and hence on transmission. However, once the pre-LLIN hypnozoite reservoir is exhausted, P. vivax prevalence can be expected to drop in line with reductions in the P. vivax EIR.
The entomological results indicate a significant impact of LLINs on transmission but a varying effect by species and study site. In the site with the largest number of mosquitoes pre-and post-LLIN (Sausi, An. farauti four in both rounds), only the P. vivax EIR was found to decrease after LLIN distribution, while the P. falciparum EIR increased. The same was observed in Finschhafen, yet with a small number of mosquitoes post-LLIN (n = Further in-depth entomological and epidemiological investigations over longer periods of time are needed to clarify the relationship between infections in mosquitoes and in humans. Sub-microscopic infections in humans may explain part of the differences in parasite species detected in humans and mosquitoes [30,40]. A better understanding is urgently required of when and where (at both macroand micro-levels) human-mosquito contacts result in malaria transmission in order to target both vector control and vector surveillance. It remains unclear how well current entomological monitoring practices capture actual transmission in space and time. A longer-term follow-up will be required to measure the full impact of LLINs on P. vivax prevalence, incidence, and transmission.
Few additional morbidity indicators were assessed in this study, all of which showed a reduction in frequency (even though not all reached statistical significance) suggesting generally positive health developments. Splenomegaly, a syndrome closely associated with chronic malaria in PNG [41,42], was less frequent among malaria patients in two sites following LLIN distribution. The trend in anaemia differed between sites. This study confirmed previous findings that anaemia is highly prevalent in some parts of PNG [17,43,44] with pre-LLIN prevalence of mild to moderate anaemia between 59 and 77 %. The prevalence of anaemia decreased in the general population (and in fever patients in Mumeng) post-LLIN. No statistically significant decrease was found in Sausi, interestingly the site in which parasite density in malaria patients remained higher than in the two other sites with clinical surveillance. However, the aetiology of anaemia is known to be multi-factorial and a causal link between LLIN and anaemia is likely to be confounded by a number of other contributing factors. Manning et al. [44] for example identified parvovirus B19 infection, P. falciparum infection, vitamin A deficiency, wasting, and incomplete vaccination as primary risk factors for severe anaemia in PNG. As expected, the vast majority of infections were asymptomatic both pre-and post-LLIN.
This study was not without limitations. Due to financial and operational constraints not all data collection components (clinical surveillance, household survey, entomology) could be implemented in all study sites, thus limiting the extent to which trends in these indicators can be compared. Household survey and entomology survey were cross sectional in nature and the clinical surveillance limited to approximately three months each round. While this approach effectively disregards seasonality, pre-LLIN and post-LLIN assessments were implemented during the same time of the year in an effort to maximize comparability. Differences in weather patterns potentially impacting on entomological indicators could not be reliably investigated due to a lack of detailed weather records (U.S. National Centers for Environmental Information; http://www.ncdc.noaa.gov/), but rainfall data from Madang airport presented by Reimer et al. [38] indicates no abnormalities in the seasonal fluctuations over the study period. However, as the comparison of different entomology datasets from Dreikikir shows, small differences in data collection may alter results and care should therefore be taken when making inferences particularly when measurements are based on small numbers of samples, in this specific case number of mosquitoes. The design of this study does not allow establishing a causal link between the LLIN and the observed health effects. However, to the best of our knowledge no other important health interventions were implemented in the study areas over the course of the two data collection rounds making it highly probable that much of the observed reductions in malaria morbidity and transmission is linked to this single intervention.

Conclusions
The distribution of LLIN had distinct effects on P. falciparum and P. vivax. Higher resilience of P. vivax may be attributed to relapses from hypnozoites and other biological characteristics favouring the transmission of P. vivax. The effect on vector species composition varied by location which is likely to impact on the effectiveness of LLINs. The relationship between prevalence, incidence and transmission is a function of several factors, including many that could not be investigated in the frame of this study. The findings of this analysis however clearly show the complexity of these relationships in