Factors associated with high heterogeneity of malaria at fine spatial scale in the Western Kenyan highlands

Background The East African highlands are fringe regions between stable and unstable malaria transmission. What factors contribute to the heterogeneity of malaria exposure on different spatial scales within larger foci has not been extensively studied. In a comprehensive, community-based cross-sectional survey an attempt was made to identify factors that drive the macro- and micro epidemiology of malaria in a fringe region using parasitological and serological outcomes. Methods A large cross-sectional survey including 17,503 individuals was conducted across all age groups in a 100 km2 area in the Western Kenyan highlands of Rachuonyo South district. Households were geo-located and prevalence of malaria parasites and malaria-specific antibodies were determined by PCR and ELISA. Household and individual risk-factors were recorded. Geographical characteristics of the study area were digitally derived using high-resolution satellite images. Results Malaria antibody prevalence strongly related to altitude (1350–1600 m, p < 0.001). A strong negative association with increasing altitude and PCR parasite prevalence was found. Parasite carriage was detected at all altitudes and in all age groups; 93.2 % (2481/2663) of malaria infections were apparently asymptomatic. Malaria parasite prevalence was associated with age, bed net use, house construction features, altitude and topographical wetness index. Antibody prevalence was associated with all these factors and distance to the nearest water body. Conclusion Altitude was a major driver of malaria transmission in this study area, even across narrow altitude bands. The large proportion of asymptomatic parasite carriers at all altitudes and the age-dependent acquisition of malaria antibodies indicate stable malaria transmission; the strong correlation between current parasite carriage and serological markers of malaria exposure indicate temporal stability of spatially heterogeneous transmission.

. Individual malaria risk may also be associated with human genetic factors [7,8,15] or with behavioural factors [6][7][8]13] including those relating to occupation [16] and travel [17]. Variations in these factors over a small area can result in spatially heterogeneous transmission, resulting in foci or hotspots of malaria infection [2,18]. Targeting hotspots may be a highly efficacious approach for malaria control [1,19]; the operational feasibility of such targeted interventions depends on the stability of malaria hotspots in space and time [2,20] and the ability to readily detect them.
In Africa, highland fringe areas have traditionally been associated with unstable malaria transmission, epidemics and unpredictable disease patterns [21,22]. Over recent decades, however, it appears that this picture has been changing [23] with studies describing instances of relatively stable malaria transmission in the Kenyan highlands [24,25], characterized by age-dependent acquisition of anti-malarial immunity [24,26] and a substantial reservoir of asymptomatic malaria infections [27]. These studies were based on passively detected malaria cases [24,25] and active surveillance in children [26,27]; a more comprehensive, community-based, assessment of parasitological and serological outcomes in all age groups is needed to establish the macro-and micro-epidemiology of P. falciparum and to identify factors associated with exposure and infection so that more targeted and specific interventions can be locally deployed.

Study site and sampling
This study was conducted in highland fringe localities (1350-1600 m altitude) in Rachuonyo South District, western Kenya (Fig. 1). The main malaria vectors in the area are Anopheles funestus and Anopheles gambiae sensu lato (s.l.) [28]. Malaria transmission is seasonal, with two peaks in malaria cases reflecting the bimodal rainfall pattern, with the heaviest rainfall typically occurring between March and June and a smaller peak between October and November each year. The study procedures have been described in detail elsewhere as part of an online clinical trial protocol [29]. Briefly, a 5 × 20 km (100 km 2 ) area was selected and divided into 400 cells of 500 × 500 m that were further subdivided in four sub-cells of 250 × 250 m. All structures in the area were geo-located manually in ArcGIS [ArcGIS 9.2; Environmental Systems Research Institute, Redlands, CA, USA] using contemporaneous high-resolution satellite data [Quickbird; DigitalGlobe Services, Inc., Denver, CO, USA]. Where possible, a maximum of 16 compounds were chosen from each 500 × 500 m cell. The aim was to obtain measurements from ≥50 individuals per 500 × 500 m and sampling of individuals was guided by pre-defined age strata (≤5 years; 6-10 years; 11-15 years; 16-25 years and >25 years) to maximize the discriminative power of serological markers of exposure [30]. To ensure maximum geographical coverage, at least one compound was selected from each 250 × 250 m sub-cell while the number of compounds selected from each of the 250 × 250 m sub-cells was weighted by the density of structures in each sub-cell.
The survey was carried out in 2011, during what is considered to be the main malaria transmission season, between June and July. After initial consent during enumeration, participating compounds were visited and the name, gender, age, residency and travel history, use of insecticide-treated nets (ITNs), indoor residual spraying (IRS) in the past 12 months and sleeping times of each compound member were recorded. All compounds where at least one adult (>20 years) and one child (<15 years) were permanent residents (defined as sleeping regularly in the structure) qualified for enrolment during the survey. The axillary temperature of each compound member was measured by digital thermometer. For all febrile individuals a rapid diagnostic test [RDT; Paracheck ® , Orchid BiomedicalSystems, Goa, India] detecting P. falciparum-specific histidine rich protein-2 was performed. A single finger prick sample was taken for haemoglobin (Hb) measurement using a HemoCue photometer [HemoCue 201+, Angelholm, Sweden] and three droplets of blood transferred onto a filter paper [3MM Whatman, Maidstone, UK] for serum and DNA [31]. All individuals with an Hb ≤11 g/dL were given haematenics; individuals with an Hb ≤6 g/dL were accompanied to a nearby health centre for further care. Febrile individuals who were parasitaemic by RDT were given artemether-lumefantrine [AL, Coartem ® , Novartis, Switzerland]; women of child bearing age who were RDT positive were assessed for pregnancy and offered a pregnancy test if deemed appropriate. Febrile children below 6 months of age and pregnant women with malaria were referred to the nearest health facility for full assessment and treatment.

PCR and serology
A combined extraction of DNA and elution of antibodies was performed on the samples collected, as described elsewhere [31]. Antibodies against P. falciparum apical membrane antigen 1 (AMA-1) and merozoite surface protein 1 19 (MSP-1 19 ) were measured in all samples by ELISA [32,33]. Parasites were detected by nested PCR targeting the 18S rRNA gene [31,34]. For logistical reasons, PCR was performed on a subset of all available samples (12,912/16,381).

Geographical information
Altitude data for study compounds were derived from a high-resolution digital elevation model (DEM; ASTER GDEM). These DEM data were also used to derive a topographic wetness index (TWI) using the method of Cohen et al. [35]. Aggregated TWI estimates were derived for a 500 m circular window around each participating compound. Locations of rivers and streams were initially estimated using topographic modelling of DEM data and were later refined manually using Quickbird satellite data. The number of digitised structures within a 500 m circular window of each compound was used as a proxy for population density.

Statistical analysis
Broad patterns of transmission intensity were described by fitting age-seroprevalence curves [32,36] to samples collected from populations residing at different altitude bands 1350-1449 m, 1450-1499 m, 1500-1549 m and 1550-1641 m and quantifying parasite prevalence and antibody prevalence in 10 m altitude bands. For quantifying transmission intensity at a finer geographical scale, two individual level measures of transmission intensity were used: (1) combined antibody prevalence, i.e. seropositivity for AMA-1 and/or MSP-1 19 ; and (2) PCR-detected parasite prevalence. Correlations between both metrics were determined using a Chi square test and hypothesis testing with significance determined where p < 0.05 and odds ratios (OR) and corresponding 95 % confidence intervals were calculated. Potential factors associated with antibody prevalence or parasite prevalence were explored using multivariate logistic regression models accounting for correlations between observations from the same compound. In these models, an equal correlation model (exchangeable) was used to specify the within-household correlation structure. All univariate analyses were adjusted for clustering of observations from the same household and age but no other factors. Adjustment for age was performed because this was a very important determinant of both parasite prevalence and antibody prevalence. For multivariate models a forward selection method was used, using a p value of 0.05 to retain variables in the model. All analyses were performed using Stata [v. 13,StataCorp].

Ethical considerations
This study was approved by the ethical committees of the London School of Hygiene and Tropical Medicine (LSHTM 5721) and the Kenya Medical Research Institute (SSC 1802 & SSC2163). Approval was sought from district medical officers, local chiefs and communities. Individual informed consent was sought from all participants or guardians of those less than 18 years old by signature or a thumbprint accompanied with the signature of an independent witness. Assent was also sought from children above 13 years of age. As defined in the Kenya national guidelines, participants below 18 years of age who were pregnant, married, or a parent were considered "mature minors" and consented for themselves.

Characteristics of the study population
In total 17,503 individuals were sampled, coming from 3213 compounds across a 100 km 2 study area. The majority of individuals resided within a narrow altitude band of 1400-1550 m (92.5 %; 16,167/17,478 [37].

Altitude and malaria risk
The proportions of individuals with fever, clinical malaria (a positive RDT in combination with temperature ≥37.5 °C), parasites detected by PCR, malaria-specific antibody responses (Fig. 2) and anaemia (Hb <11 g/gL) were all negatively associated with increasing altitude (Table 1; p < 0.001 for all comparisons). The median age of individuals diagnosed with clinical malaria was 6 years (interquartile range 3-10 years), although 12.5 % (37/295) individuals with clinical malaria were >15 years of age (range 0-46 years); 93.2 % (2481/2663) of PCRdetected malaria infections were sub-clinical. The mean age of individuals with clinical malaria was not associated with altitude (p = 0.40). When altitude was categorized in bins of 10 m, there was a strong negative association between altitude and PCR parasite prevalence ( Fig. 3a; r = 0.92, p < 0.0001) and malaria antibody prevalence (r = 0.92, p < 0.0001). Age-seroprevalence curves were fitted for antibody responses to MSP-1 19 and/or AMA-1 for the four altitude bands and indicated a clear gradual decline of the seroconversion rate with increasing altitude (Fig. 3b). There are no apparent observed 'steps' in age-seroprevalence curves that can indicate changes in exposure due to intervention and/or age-associated behavioral changes patterns such as travel to malaria endemic regions outside the area of residence [38].
In univariate analysis, individual-level parasite prevalence and antibody prevalence was associated with age. When adjusting for age, the following factors were associated with parasite prevalence and antibody prevalence: the presence of eaves in the sleeping room, reported bed net use, construction of the sleeping room walls with mud, distance to water, altitude, the proportion of parasite positive individuals within a 500 m radius of the compound and the proportion of malaria antibody positive individuals within a 500 m radius of the compound ( Table 2). In a multivariate model that was constructed using forward selection of variables, nPCR parasite prevalence was negatively associated with bed net use, altitude and maximum TWI, and positively associated with mud walls, the presence of open eaves (Table 2). Similarly, malaria antibody prevalence was negatively associated with bed net use, distance to water and altitude, and positively associated with the presence of open eaves and minimum TWI (

Discussion
This study, conducted in the Kenyan highlands, shows the occurrence of ongoing stable malaria transmission across an altitudinal range of 1350-1600 m, characterized by marked spatial heterogeneity. The age-dependent acquisition of malaria antibodies, strong correlation antibody prevalence and current parasite prevalence, along with the considerable asymptomatic reservoir of P. falciparum infections in all age groups, suggests that malaria transmission is relatively stable in the study setting.
Much of the highlands of East Africa represent fringe regions between stable and unstable malaria transmission; seasonal and spatial patterns in malaria transmission are affected to some degree by annual variations in rainfall but primarily by ambient temperature [39]. The notion that malaria is largely absent in areas higher than 1500 m [40] has been challenged by findings of a large asymptomatic reservoir of malaria infections at altitudes [27] and an age-dependent acquisition of clinical immunity to malaria infections in highland communities [24]. In this study area area at 1350-1600 m above sea level, the special epidemiology of malaria infections was determined and markedly heterogeneous malaria transmission was observed. This is commonly observed in areas of low endemicity [2,20,25] and heterogeneity in clinical malaria cases has previously been reported in the Kenyan highlands [25]. This study adds detail to previous findings by describing the fine-scale spatial distribution of asymptomatic parasite carriage and immunological evidence of previous malaria exposure in a highland area. PCR-detectable P. falciparum infections were very common in this highland setting, apparently asymptomatic and negatively correlated with altitude. Parasite prevalence by PCR was 27. Only a small fraction of these P. falciparum infections resulted in fever at the time of sampling and apparently asymptomatic parasite carriage was prevalent at all altitudes [27] and in all age-groups [41]. Recently, a manuscript summarized the evidence on the clinical consequences of chronic low density infections, arguing that many infections are incorrectly classified as asymptomatic and have considerable health consequences in terms of anaemia, chronic inflammation, school performance and bacterial infections [42]. Since the cross-sectional design of this study does not allow to determine whether the nPCR detected infections had clinical implications for the study population, it cannot be concluded with certainty whether the detected infections were indeed asymptomatic. However, concurrent clinical symptoms were reported by a small minority of the examined population and infections were probably chronic in nature. Although individuals with limited previous exposure may harbor low density infections [41], the high prevalence of apparently asymptomatic parasite carriage, absence of travelling as obvious risk factor for malaria and the gradual age-and altitude-dependent acquisition of antibody responses to P. falciparum antigens at different altitudes suggests stable local malaria transmission in the area.
Whilst travel, a known risk factor for malaria in highland areas [43,44], was not significantly associated with individual risk of malaria infection or antibodies, several household factors such as the presence of open eaves and mud walls [6,9] were statistically significant predictors of malaria risk. This suggests that relatively simple household improvements may decrease malaria risk [45,46] in a region where the perceived and measured indoor exposure to malaria vectors is low [28]. Altitude, distance to . Symbols indicate parasite prevalence estimates for children below 5 years of age, 6-10 years, 11-15 years, 16-25 years and >26 years. Symbols are plotted at the median age for the different categories; for the highest age category parasite prevalence is plotted at 35 years water and the proportion of antibody or parasite positive individuals in the immediate vicinity of a compound were statistically significant and biologically plausible factors associated with malaria risk. The latter could suggest considerable between household transmission.

Conclusions
Evidence of relatively stable malaria transmission in this site at 1350-1600 m altitude was observed. Altitude was a major driver of malaria transmission in the study area, even across narrow altitude bands. Although malaria risk was spatially heterogeneous, the strong correlation between current parasite carriage and serological markers of malaria exposure and other established risk factors for malaria indicate temporal stability of geographical patterns in malaria exposure. The fine scale heterogeneity in this low-endemic setting may reflect a likely scenario for more endemic areas with active and effective control programmes as they reduce transmission to increasingly low levels. A priori knowledge of the factors that influence residual malaria in foci of low transmission is likely to further expedite control and elimination attempts.