Skip to main content

A cross-sectional study of asymptomatic Plasmodium falciparum infection burden and risk factors in general population children in 12 villages in northern Uganda

Abstract

Background

Plasmodium falciparum malaria is an important cause of morbidity in northern Uganda. This study was undertaken to assess village-, household-, and individual-level risk factors of asymptomatic falciparum malaria in children in 12 villages in northern Uganda.

Methods

Between 10/2011 and 02/2014, 1006 apparently healthy children under 16 years old were enrolled in 12 villages using a stratified, multi-stage, cluster survey design and assessed for P. falciparum malaria infection using the rapid diagnostic test (RDT) and thick film microscopy (TFM), and structured interviewer-administered questionnaires. Associations between weighted P. falciparum malaria prevalence (pfPR), based on RDT, and covariates were estimated as odds ratios and 95% confidence intervals (ORs, 95% CIs) using logistic models accounting for the survey design.

Results

Among 942 (93.5%) children successfully tested, pfPR was 52.4% by RDT and 32.7% by TFM. Overall pfPR was lower in villages where indoor residual insecticide spray (IRS) was, versus not, implemented (18.4% versus 75.2%, P < 0.0001). However, pfPR was heterogeneous both within IRS (10.6–34.8%) and non-IRS villages (63.6–86.2%). Elevated pfPR was associated with having a sibling who was RDT positive (OR 5.39, 95% CI 2.94–9.90, P = 0.0006) and reporting a fever at enrollment (aOR 4.80, 95% CI 1.94–11.9, P = 0.0094). Decreased pfPR was associated with living in an IRS village (adjusted OR 0.06, 95% CI 0.04–0.07, P < 0.0001), in a household with one (aOR 0.48, 95% CI 0.30–0.76) or more than one child below 5 years (aOR 0.23, 95% CI 0.12–0.44, Ptrend = 0.014), and reporting keeping a goat inside or near the house (aOR 0.42, 95% CI 0.29–0.62, P = 0.0021).

Conclusions

The results show high but heterogeneous pfPR in villages in northern Uganda, confirm significantly decreased pfPR associated with IRS implementation, and suggest significant associations with some household characteristics. Further research is needed to elucidate the factors influencing malaria heterogeneity in villages in Uganda.

Background

Plasmodium falciparum malaria remains an important cause of morbidity and mortality globally, particularly in children aged less than 5 years old in countries in sub-Saharan Africa [1, 2], including Uganda [3]. Scale-up of malaria control interventions since the mid-2000s has resulted in significant declines in malaria morbidity and mortality in many countries [4,5,6]. However, the changes in malaria morbidity and mortality are heterogenous in different countries, characterized by increases or slower declines in some countries [7], and heterogenous changes in the incidence patterns in different geographical areas of the same country [8]. However, malaria surveillance is mostly based on acute malaria mortality and morbidity at health facilities, which may not correlate perfectly with impact on asymptomatic malaria infections, which comprise ≥ 75% of all infections [9], and represent an important source of new infections in a population [7]. A better understanding of the epidemiology of asymptomatic malaria in high-burden countries may help to design interventions to reduce the local burden and clinical impact of malaria [10].

Continuing intense P. falciparum malaria transmission has been reported in northern Uganda [11], an area home to some of the world’s highest malaria transmission rates (annual entomological inoculation rates of 400–1500 infectious mosquito bites per person, per year [12]). Between 2008 and 2015, the Ugandan Government, with the support of international donor partners, implemented 4–6 monthly indoor residual spraying (IRS) cycles and distribution of long-lasting insecticide-impregnated bed nets (LLIN) to the general population in 10 districts in northern Uganda [13, 14], leading to strong suppression of malaria prevalence in the targeted districts [13, 14]. The dramatic declines in outpatient malaria cases in children aged below 5 years old led to optimism that prompted withdrawal of IRS in 2014. However, this decision was premature because it was followed by a rapid increase of malaria cases [15]. The rapid reversal of trends in malaria morbidity following withdrawal of IRS in northern Uganda highlights the likely role of undetected subclinical infections in contributing to new cases of malaria [7]. The patterns of asymptomatic falciparum infection in northern Uganda before or after IRS are not well known. Most malaria data are based on national surveys, which are not sufficiently granular to provide local insights [16], focus on participants in a restricted age range, usually children below 5 years old [14], or on clinical cases [17], which provide limited information about the asymptomatic cases that represent the main reservoir of infection that spawns new clinical cases. This study was undertaken to obtain village- household-, and individual-level information about risk factors of asymptomatic malaria infection in children in northern Uganda, which would complement data based on broad sampling at a national [14, 16] and regional levels [13]. However, the study had an additional motivation: P. falciparum malaria is thought to be important in endemic Burkitt’s lymphoma (eBL) [18]. Thus, the study also aimed to generate baseline data to enable precise evaluation of the role of malaria in eBL in Uganda.

Methods

Study setting

Between October 2011 and February 2014, a cross-sectional survey of children in 12 random villages in 22 districts Northwest and North-central regions of Uganda (Fig. 1). The children were enrolled as pilot population controls in a case–control study entitled the Epidemiology of Burrkitt Lymphoma in East African children and Minors (EMBLEM) study [13, 18], which is part of a larger study investigating risk factors of eBL in Uganda, Tanzania, and Kenya [19]. The target regions span from 621 to 2000 m above sea level and include areas covered by savannah, slow flowing rivers, swamps, lakes, which are separated into two regions (Northwest and North-central) by a gorge of the East African Rift valley. Both regions experience warm (20–30 °C) and wet climate (1250–1500 mm of rainfall per year [20] in two rainy seasons: April to June and September to December), which is conducive for mosquito breeding and perennial malaria transmission [12]. The malaria prevalence in these regions is high, but geographically heterogenous (prevalence ranging from 3.2 to 75.8% [13]). As noted in the introduction, IRS and LLIN were implemented in 10 districts in the study area between 2008 and 2015 [21], and universal LLIN has been implemented since 2017.

Fig. 1
figure 1

Map of Uganda showing the study area in northern Uganda, including district and census Enumeration area (EA) boundaries. The zoom out shows all-season geographical features, including lakes, rivers or streams. The 12-pilot population control (PPC) EAs are shown according to stratification characteristics (red: near surface water and urban, yellow: far from water and urban, sky blue: near surface water and urban, and azure blue: near water and rural; see “Methods”). The map also shows the 88 EAs that were included in the previous study by Maziarz et al. [13] as gray shaded areas to show their relative geographical distribution in relation to the 12 EAs included in the current study

Sampling design

The target population of children (0–< 16 years old) in northern Uganda was estimated to be 3,031,494 in 2015, based on the national population census data [20]. This population was sampled using a stratified multi-stage cluster sampling design (Fig. 2) [13]. In the first sampling stage, 12 census enumeration areas (EAs) were randomly selected from 4 strata, defined by ‘low population-density’, ‘high population-density’, ‘near water’, and ‘far from water’, from a list of all EAs obtained from the Uganda Bureau of Statistics (UBOS) [20]. These strata are expected to be associated with malaria transmission [22, 23]. Population density of an EA was categorized as low if it was below the median population count (n = 2683), and high if it was equal or higher than the median population count [13]. Proximity of an EA to water was defined as near surface water when the EA boundary was next to or within 500 m of an all-season swamp, river, or lake, based on distances estimated from national maps incorporating geographical information metadata; otherwise, it was defined as being far from water. Four EAs were sampled from the low-density near water stratum and the low density far from water stratum, and 2 EAs from the high density near water stratum and the high density far from water stratum. In the second stage, one village per EA was randomly selected and a household survey conducted, and eligible children in the household invited to participate. Children were eligible if they were aged 0–< 16 years old, were usual residents of the household, and were apparently healthy, i.e., not having symptoms requiring hospital care. Children were not eligible of they had a cancer diagnosis or symptoms requiring hospital treatment.

Fig. 2
figure 2

Flow chart showing the stratified, multi-stage cluster sampling design used to sample healthy children aged 0–15 years in 12 randomly selected villages in north-central and northwest regions of Uganda

Participant enrollment

An experienced field team consisting of laboratory technicians, interviewers, and a local guide, led by a local medical doctor, visited consecutive households in each village to enroll eligible children. For practical reasons, enrollment in each village continued until a sample size of at least 70 children was obtained. The sample size per village was based on a target sample size of 840–1080 children from 12 villages, which was needed to enable precise estimation of population of pfPR in healthy children in the 12 villages as well as in each village. This calculation assumed overall malaria prevalence of 50%, that each village would have at least 30 households with an average household size of 2–3 eligible children, yielding a sample size of 60–90 children per village. Enrollment in each village lasted 2–3 weeks per village, and between October 2011 to February 2014 for the 12 villages. The period of enrollment was prolonged because the study design stipulated that about half of the villages should be enrolled during the dry season and the other half during the wet season. This necessitated phased implementation to enable the single team would be able to mobilize and safely enroll children in all the villages, which were remote, located far away from one another and from the EMBLEM field laboratories (Fig. 1). Written informed permission for the child to take part was obtained from the child’s guardian (usually one or both parents), and assent from children aged 8 years or older. Structured questionnaires were administered to record each child’s age, sex and household characteristics, parental educational level, malaria prevention methods, including ownership and use of LLIN the night before interview and application of IRS in the house, and a history of outpatient or inpatient malaria treatment (≤ 6, 7–12, or ≥ 13 months ago).

Blood samples and malaria testing

Venous blood for research (10 mL) and clinical (4 mL) tests was drawn in EDTA tubes from each child. The research blood samples were transported in cold boxes to local EMBLEM field laboratories within 2 h from sampling, and centrifuged for 15 min at 1300g to separate plasma, buffy coat, and red cell fractions that were stored in barcoded cryovials at − 80 °C. Clinical samples were immediately tested for malaria parasites by experienced local technicians using thick film microscopy (TFM) to visualize asexual malaria parasite forms and commercial rapid diagnostic tests (RDT) (MALARIA DUAL kits, ICT Diagnostics, Muizenberg, Cape Town, South Africa) to detect malaria antigens [13]. TFM was performed on slides stained with 10% Giemsa solution for 10 min, and the visualized parasites were counted against 200 white blood cells (WBCs) and the results expressed as parasites/µL of blood, based on the measured WBC count. The RDT kits used in the current study detect the P. falciparum-specific malaria histidine-rich protein 2 (Pf-HRP2) and the pan-lactate dehydrogenase (pLDH) antigen and have a reported sensitivity and specificity of 92–100% in Uganda [24]. Because malaria antigens may appear in blood before asexual parasites and remain in blood for several weeks after asexual parasites have been cleared from blood, RDTs were used as the primary outcome measure to capture malaria infection defined by both visualized parasites and circulating antigens [25]. Although RDT tests can yield false-negative results, particularly for parasites that have deleted pfhrp2 or pfhrp3 genes that code for the RDT antigens [26], the frequency of such parasites is low (< 1.1) and was considered unlikely to distort the general malaria epidemiological patterns. Thin film smears were examined to identify Plasmodium species. Because ~ 98% of the parasites were P. falciparum, this species is assumed hereafter.

Data management

Questionnaires and result forms were reviewed for completeness and accuracy, and computerized by DataFax or computer data entry. DataFax uses intelligent character recognition built-in capabilities to perform data entry and reduce data entry errors. Logical consistency checks were performed on the electronic data, and data queries arising were corrected before creating analysis files.

Statistical analysis

The response rate was calculated as the proportion of enrolled children out of all eligible children in the selected households that were enrolled. The outcome variable for this analysis was P. falciparum prevalence (pfPR), based on RDT positivity, which captures infection defined by visualized asexual parasites and circulating falciparum antigens [25]. Although polymerase chain reaction (PCR) would provide a more accurate result for antigenemia. One child who reported a fever on their questionnaire and had a high parasite count on TFM (> 2500 parasites/uL) was considered to have clinical malaria excluded from analysis. Children who reported a fever on their questionnaire but were malaria negative (n = 6) or had low parasitemia (n = 11, parasite load ranging from scanty to 1200 parasites/uL) were considered to have incidental malaria and not excluded from analysis. The geometric mean parasite density (GMPD) per µL was assessed overall and by age group to confirm low-level parasitemia among children with visualized parasites, as would be expected in high-endemic areas [9].

The results were weighted by the inverse probability of being sampled into the study, and are weighted back to the general population of children aged 0–< 16 years old in Northwest and North-central Uganda in 2015 (projected to be 3,138,360). Weights were trimmed by replacing the value of the weights in the highest 3% of the weight distribution with the value of the weight at the 97th percentile to minimize the impact of outlier weights. Stratification and clustering at the village and household levels were accommodated in the variance computations of estimates. The reported results are weighted estimates, unless stated otherwise.

The analyses were performed using the “survey” package [27, 28] (v. 3.32-1) in R (version 3.4.3, r-project.org) following methods previously described [13]. Descriptive analyses of weighted pfPR are reported overall, by stratum, by village, and by demographical and geographical variables as in [13]. Associations with pfPR were assessed for: (a) household variables reflecting crowding (number of other children in the household), parental characteristics (mother’s income, parents’ occupation), malaria infection in the enrolled household members (any child with infection, infection in a child below 5 years old, infection in a younger sibling); (b) individual characteristics, including history of fever due to malaria or non-malaria conditions (at enrollment, in the past 6 months, in the past 12 months); and (c), use of herbs for treatment of skin or gum disease and lifetime number of admissions. When considering the number of other children in the household or the number of other children in the household below 5 years old, this variable was coded as 0 if there was no other child in the household or no child below 5 years in the household, otherwise it was coded as the number of other children or the number of children below 5 years in the household. Associations between pfPR and keeping different types of animals inside the house or nearby were evaluated because mosquito blood-feeding preferences (anthropophily/zoophily) could influence the risk of malaria positivity [29]. Unadjusted and adjusted odds ratios (ORs), standard errors [30] and Wald-type 95% CIs (95% CIs) [31] of association of pfPR with each variable were calculated. Forward stepwise logistic regression was used to construct adjusted models, using the covariates that in unadjusted models had associations with prevalence with P < 0.05. Covariates with several levels were coded with dummy variables for the categories in the univariate analysis and using trend coding in the adjusted analyses. As the objective of these analyses were descriptive and for hypothesis generation, statistical tests were not adjusted for multiple testing. A two-sided P < 0.05 was considered statistically significant.

Results

Demographic characteristics of study population

As shown in Fig. 2, 1007 (86.2%) of 1168 eligible children in 354 (79.5%) of 445 households randomly selected from 1121 eligible households in 12 villages were enrolled. These children represent a trimmed weighted sample of 2,829,988 children aged 0–15 years (Additional file 1: Table S1 and Additional file 2: Fig. S1). The demographic distribution was 45.4% being males and 54.6% females; 42.7% were children aged less than 5 years old, 34.4% aged 5–9 years old, and 23.9% aged 10–15 years old. The size and age distribution of the weighted population were comparable to the corresponding population in the study area in 2015, based on data from UBOS [20].

Malaria prevalence by rapid diagnostic test and validation by microscopy

Of 1007 enrolled children, 942 (93.5%) were successfully tested for malaria. One subject who reported a fever at enrollment and had parasite load > 2500 parasites/uL was excluded from analysis. In the remaining set, the weighted pfPR was 52.4% (95% CI 19.9–84.8%) by RDT read by experienced local technicians and 32.7% (9.3–56.2%) by TFM performed by experienced local technicians (Table 1). The sensitivity of RDT to detect asexual parasites was 97.5% (95% CI 95.8–99.7%) and the specificity was 69.7% (95% CI 43.4–96.1%), compared to TFM by experienced local technicians.

Table 1 Weighted Plasmodium falciparum results showing the validation of rapid diagnostic test to detect parasitaemia using thick film microscopy by experienced local technicians, among healthy children in northern Uganda

The overall and age-specific geometric mean parasite density

The overall weighted GMPD was low 733.3 parasites/µL. GMPD was significantly higher in children who were both RDT and TFM positive than those who were TFM-positive but RDT negative (768 versus 107 parasites/µL, P = 0.002). GMPD did not vary by sex (P = 0.167), and was not significantly different in children aged below 5 years versus older (955 versus 615 parasites/µL, Pdifference = 0.201). GMPDs showed modest variation across the villages (189–724 parasites/µL, except for two villages with an average of 1357 and 1807 parasites/µL).

Malaria prevalence in the 12 villages, by IRS status and season

The weighted pfPR in villages varied from 10.6% in village 19 in the North-central region to 86.2% in village 40 in the northwest region (Table 2). pfPR was lower in IRS villages and higher in non-IRS villages (18.4% versus 75.2%, P < 0.0001), but it varied both within IRS (10.6–34.8%) and non-IRS villages (63.6–86.2%). No overall differences were noted between wet versus dry seasons (66.9% versus 46.2%, P = 0.40). However, within IRS villages, pfPR was lower in villages enrolled in the dry season compared to those enrolled in the wet season (10.6–17.6% versus 20.1–34.8%). No differences in pfPR by season were observed within non-IRS villages (dry: 63.6–75.2% versus wet: 74.8–86.2%; Table 2).

Table 2 Weighted prevalence of asymptomatic Plasmodium falciparum malaria in 12 random villages in northern Uganda, by season and location in a district where indoor residual insecticide spraying (IRS) was implemented or not

pfPR and geometric mean parasite density, by population density strata, season, age

The stratum-specific weighted pfPR by season are shown in Fig. 3. Within the low-population density stratum, pfPR in the wet and dry season was similar (67.3% versus 70.8%, P = 0.802). Conversely, within the high-population density stratum, pfPR was higher in the wet season than in the dry season, but not statistically significant (62.1% versus 37.9%; Pheterogeneity = 0.589; Fig. 3a). GMPD did not vary by season in the low- or high-population density strata (Fig. 3b), but the GMPD tended to be higher in the dry season than the wet season in the high-density stratum.

Fig. 3
figure 3

Bar graphs showing weighted malaria weighted per cent Plasmodium falciparum parasite prevalence, based on the RDT, by season (wet or dry) in low- and high-population density villages (a) and the GMPD/µL among microscopy-positive children by visit season in low- and high-population density villages (b) among apparently healthy children enrolled in 12 random between November 2011 and February 2014 in north-central and northwest Uganda. Note: Orange shading is used for dry season months, while blue shading is used for wet season months. Wet and dry seasons are based on categorization by the Uganda Bureau of Statistics and generally correspond to ≥ 10 days/month for wet months and < 10 days/month for dry months. The unweighted number of participants in each group is shown

The age-specific pfPR was higher in non-IRS versus IRS villages, but slight differences in the actual peaks. pfPR peaked between 6 and 8 years old in IRS villages and slightly later between 9 and 11 years old in non-IRS villages (Fig. 4a). Among children with visualized parasites, the age-specific-GMPD was higher in non-IRS villages (Fig. 4b), in children up to 6–8 years of age and then decreased rapidly in older children. The GMPD in children in IRS villages was very low for all ages, although a spike was observed in children aged 6–8 years old (Fig. 4b), probably due to random fluctuation.

Fig. 4
figure 4

Bar graphs showing age-group patterns of weighted Plasmodium falciparum parasite prevalence, based on RDT, stratified by whether the village is in a district where indoor residual insecticide spraying (IRS) against mosquitoes was implemented (IRS district: red color) or not (non-IRS district: green color) (a) and the GMPD parasites/µL (b). Results are from apparently healthy, microscopy-positive children enrolled in 12 randomly selected villages between October 2011 and February 2014 in northern Uganda. The unweighted number of participants in each age-group is shown. In b, the open circle, or rectangles shows the GMPD parasites/µL, the lines show the 95% CIs of the GMPD. GMPD results were available on 328 with positive results of 339 (includes 1 subject who was RDT negative) tested by thick film microscopy; thick film negative subjects were not included (see Table 1)

Household- and individual-characteristics associations with weighted pfPR in children

The associations between pfPR and children’s demographical, malaria prevention, illness history, lifetime malaria treatment, environmental, and household characteristics are shown in Table 3. No differences in pfPR were noted between females and males (53.3% versus 51.4%; P = 0.768), by distance of the home from the source of drinking water being ≥ 1 or < 1 km (55.7% versus 49.4%; P = 0.723). pfPRwas associated with keeping goats inside or near the house (49.8% versus 60.5%, OR 0.65, 95% CI 0.48–0.87; P = 0.021).

Table 3 Plasmodium falciparum parasite prevalence (pfPR) among children 0–15 years old enrolled between October 2011 and Feb 2014 in 12 villages in north-central and northwest regions of Uganda and associations with sex and those characteristics that resulted in P < 0.05 in univariate logistic models

Compared to living in a non-IRS village, living in an IRS village was associated with decreased pfPR (18.5% versus 75.2%, OR 0.07, 95% CI 0.05–0.11). Consistent with this result, pfPR was decreased in children for whom their questionnaire reported that IRS was applied to their house, as compared to no application (23.7% versus 59.0%, OR 0.22, 95% CI 0.05–0.86). Consistent with previous results [13], sleeping under an LLIN on the night before enrollment, versus not, was not associated with decreased pfPR (64.6% versus 48.8%, P = 0.395). Compared to having no other child in the household, having one or more other children in the household, as was associated with decreased pfPR (Pheterogeneity < 0.0001; Table 3). The results were similar when the number of children below 5 years old was considered (Pheterogeneity = 0.0002).

Reporting a history of non-malaria related fever in the past 6 months, versus not, was unrelated to pfPR (60.6% versus 51.4%, P = 0.384). However, reporting four or more malaria-related fever episodes in the past 6 months, versus none, was significantly associated with elevated pfPR (OR 2.80, 95% CI 1.11–7.05, Pheterogeneity = 0.0003). Compared to children reporting a history of inpatient malaria treatment in the past 12 months, pfPR was lower in those reporting a history of inpatient malaria more than 12 months ago (OR 0.72) and those reporting no such a history ever (OR = 0.47, Ptrend = 0.026; Table 3). In contrast, compared to children reporting a history of outpatient malaria treatment in the past 12 months, the pfPR was elevated in those reporting a history of outpatient treatment of malaria more than 12 months ago (OR 1.83, 95% CI 1.15–2.91), but decreased in those reporting no such history ever (OR 0.57, 95% CI 0.41–0.80, Table 3).

A few children (n = 19) who reported a fever on their questionnaire were deemed to have incidental malaria at the time of enrollment. Compared to those who did not report a fever, reporting a fever was associated with pfPR (OR 5.38, 95% CI 2.01–14.39). pfPR was associated with having at least one sibling who was RDT positive (versus none: OR 4.24 95% CI 1.61–11.15) and with having a younger sibling who was RDT positive (versus none: OR 2.20 95% CI 1.40–3.46). Season, parental education and occupation, herbal treatment, and keeping other animals (other than goat) inside or near the house was not associated with pfPR (Table 4).

Table 4 Patterns of P. falciparum parasite prevalence according to environmental, parental, and household characteristics in children aged 0–15 years old enrolled in 12 villages in north-central and northwest regions of Uganda between October 2011 and Feb 2014

Only 6 variables of 14 evaluated using stepwise logistic regression remained significantly associated with pfPR in the multivariable model. Specifically, pfPR was inversely associated with living in an IRS village (adjusted OR 0.06, 95% CI 0.04–0.07), with having two or more other children in the household (versus none, Ptrend < 0.0001), or 2–4 children in the household aged below 5 years of age (versus none, Ptrend = 0.014), and keeping a goat inside or near the house (aOR 0.42, 95% CI 0.29–0.62). Conversely, pfPR was positively associated with reporting a fever at enrollment (aOR 5.02, 95% CI 1.95–12.92) and having at least one sibling who was RDT positive (aOR 5.39, 95% CI 2.94–9.90).

Discussion

The current analysis, conducted to generate recent village, household, and individual-level data about incidental malaria infections in children living in 12 villages in northern Uganda, confirms the high overall pfPR in the villages studied. Similar to findings based on a study of children in 100 villages [13], the study shows substantial heterogeneity on pfPR in different villages. The study also confirmed that IRS was major factor driving pfPR patterns, but noted substantial heterogeneity both within IRS and non-IRS villages. Further research is needed to elucidate factors influencing malaria heterogeneity in villages in Uganda, regardless of their IRS status. The pfPR patterns in non-IRS villages did not vary by season suggesting that the prevalence of incidental malaria is saturated in the dry and wet season due to perennial transmission and a high level of herd immunity [32]. Conversely, although not statistically significant, the pfPR in IRS villages increased two-fold from dry to wet season. This pattern suggests that the vector-host ratio in the IRS district villages increases during the wet season [33]. This pattern suggests that malaria prevalence increases in the wet season, presumably due to a corresponding increase in the vector-host ratio, and predicts that malaria prevalence will increase sharply during the wet season when IRS is withdrawn. Consistent with this, malaria epidemics were reported in northern Uganda following the withdrawal of IRS in 2014 [15], and the sharpest increases were noted during the wet seasons. The results presented here suggest that preemptive IRS before onset of the wet season or in areas where IRS has been withdrawn could offer some benefits to the population.

The use of insecticide-treated bed nets the night before interview was not associated with pfPR, confirming the results reported by other studies conducted in northern Uganda [13, 34] and elsewhere [35]. This finding may be pointing to the emergence of mosquito vector metabolic resistance to pyrethroids impregnated into bed nets [36] and a threat against the continued effectiveness of insecticide-treated bed nets in the fight against malaria [37, 38]. These results call for careful monitoring of malaria case activity at a district level and proactive correlation of uptake of universal distribution of insecticide-treated bed nets with clinical case activity, such as using electronic health management information systems (eHMIS) [39]. The results also point to a need for more research to investigate pyrethroid resistance and alternative chemicals if insecticide-treated bed nets are to be maintained as a key malaria control tool in malaria endemic countries [5, 6].

The interest in malaria patterns in northern Uganda had an addition motivation. Malaria is thought to be important in eBL [18], and a better understanding of its patterns in needed to improve the design of the EMBLEM study [40], and its goal to understand and control for geographical patterns of eBL [41]. The finding of mostly low-grade malaria parasitaemia in the BL-age children in this study is consistent with the hypothesis that children at risk of eBL generally have a high level of anti-parasite immunity [42]. If so, the baseline risk factors for pfPR (low grade malaria parasitaemia) in children with established immunity may be relevant for investigating and interpreting malaria-related risk factors for eBL, particularly in northern Uganda [43, 44]. For example, the total number of other children (or those below 5 years) living in a household was inversely associated with prevalence of low-grade malaria infection would be hypothesized to have an inverse relationship with eBL. Conversely, having a child with confirmed malaria infection in the household was associated with low-grade malaria infection, may confound associations with eBL. The finding that keeping goats inside or near the house was inversely associated with pfPR was unexpected, but interesting. This finding may be a clue about the role of zooprophylaxis or zoopotentiation [29], due to feeding preference of mosquito vectors for humans of animals [29], influencing malaria risk. Donnelly et al. [45] reported that keeping goats or chicken inside the house was significantly associated decreased risk of malaria in Kanungu District in Western Uganda. Bulterys et al. [46] reported that a goat: human ratio ≥ 1 was associated with reduced the risk of malaria infection in a study conducted in a village in Zambia. These findings suggest that keeping of certain animals inside or near the house may influence malaria risk, perhaps, by drawing vectors away from humans. Taken together, these results will enable hypothesis-driven research about factors associated with eBL risk in northern Uganda, and elsewhere.

The limitations of the study include reliance on cross-sectional data and lack of geospatial data to complement questionnaire data. The lack of recent entomological data limits the discussion about abundant vectors in northern Uganda, and specifically about feeding preference of the vectors prevalent in northern Uganda and whether their relative distribution has changed from before to after the implementation of malaria suppression [21]. Surveys conducted in the 1960s and 1970s suggested the abundant species as Anopheles gambiae and Anopheles funestus, which are antropophilic [47], while Anopheles arabiensis, which is feeds both on humans and animals was less abundant [48], but the current patterns are unknown.

Conclusions

The study confirmed high, but heterogenous asymptomatic pfPR in children in 12 villages in northern Uganda. The pfPR was strongly influenced by IRS, but showed substantial variation both within IRS and non-IRS villages. pfPR was inversely associated with living in an IRS village, having a greater number of other children, or children aged below 5 years of age in the household, and keeping a goat inside or near the house. Conversely, pfPR was positively associated with reporting a fever at enrollment and having at least one sibling who was RDT positive. Further research is needed to elucidate factors influencing malaria heterogeneity in villages in Uganda, and to determine the value of developing public health messages that stress household characteristics associated with pfPR.

Abbreviations

pfPR:

Plasmodium falciparum prevalence

pfhrp2 :

Plasmodium falciparum histidine-rich protein 2

pfhrp3 :

Plasmodium falciparum histidine-rich protein 3

GMPD:

geometric mean parasite density

OR:

odds ratio

CI:

confidence interval

AEIR:

annual entomological inoculation rate

BL:

Burkitt lymphoma

IRS:

indoor residual spraying

EA:

enumeration area

UBOS:

Uganda Bureau of Statistics

TFM:

thick film microscopy

WBC:

while blood cell

RDT:

rapid diagnostic test

pLDH:

pan lactate dehydrogenase

USHS:

Ugandan shillings

US$:

United States dollar

References

  1. Liu L, Oza S, Hogan D, Chu Y, Perin J, Zhu J, et al. Global, regional, and national causes of under-5 mortality in 2000–15: an updated systematic analysis with implications for the sustainable development goals. Lancet. 2016;388:3027–35.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Snow RW, Omumbo JA, Lowe B, Molyneux CS, Obiero JO, Palmer A, et al. Relation between severe malaria morbidity in children and level of Plasmodium falciparum transmission in Africa. Lancet. 1997;349:1650–4.

    Article  PubMed  CAS  Google Scholar 

  3. Talisuna AO, Noor AM, Okui AP, Snow RW. The past, present and future use of epidemiological intelligence to plan malaria vector control and parasite prevention in Uganda. Malar J. 2015;14:158.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Snow RW, Guerra CA, Noor AM, Myint HY, Hay SI. The global distribution of clinical episodes of Plasmodium falciparum malaria. Nature. 2005;434:214–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Hemingway J, Shretta R, Wells TN, Bell D, Djimde AA, Achee N, et al. Tools and strategies for malaria control and elimination: what do we need to achieve a grand convergence in malaria? PLoS Biol. 2016;14:e1002380.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  6. Rabinovich RN, Drakeley C, Djimde AA, Hall BF, Hay SI, Hemingway J, et al. malERA an updated research agenda for malaria elimination and eradication. PLoS Med. 2017;14:e1002456.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Snow RW, Sartorius B, Kyalo D, Maina J, Amratia P, Mundia CW, et al. The prevalence of Plasmodium falciparum in sub-Saharan Africa since 1900. Nature. 2017;550:515–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Jagannathan P, Muhindo MK, Kakuru A, Arinaitwe E, Greenhouse B, Tappero J, et al. Increasing incidence of malaria in children despite insecticide-treated bed nets and prompt anti-malarial therapy in Tororo, Uganda. Malar J. 2012;11:435.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Greenwood BM. Asymptomatic malaria infections–do they matter? Parasitol Today. 1987;3:206–14.

    Article  PubMed  CAS  Google Scholar 

  10. malERA Refresh Consultative Panel on Characterising the Reservoir, and Measuring Transmission. An updated research agenda for characterising the reservoir and measuring transmission in malaria elimination and eradication. PLoS Med. 2017;14:e1002452.

    Article  Google Scholar 

  11. Proietti C, Pettinato DD, Kanoi BN, Ntege E, Crisanti A, Riley EM, et al. Continuing intense malaria transmission in northern Uganda. Am J Trop Med Hyg. 2011;84:830–7.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Okello PE, Van Bortel W, Byaruhanga AM, Correwyn A, Roelants P, Talisuna A, et al. Variation in malaria transmission intensity in seven sites throughout Uganda. Am J Trop Med Hyg. 2006;75:219–25.

    Article  PubMed  Google Scholar 

  13. Maziarz M, Kinyera T, Otim I, Kagwa P, Nabalende H, Legason ID, et al. Age and geographic patterns of Plasmodium falciparum malaria infection in a representative sample of children living in burkitt lymphoma-endemic areas of northern Uganda. Malar J. 2017;16:124.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Ssempiira J, Nambuusi B, Kissa J, Agaba B, Makumbi F, Kasasa S, et al. Geostatistical modelling of malaria indicator survey data to assess the effects of interventions on the geographical distribution of malaria prevalence in children less than 5 years in Uganda. PLoS ONE. 2017;12:e0174948.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Okullo AE, Matovu JKB, Ario AR, Opigo J, Wanzira H, Oguttu DW, et al. Malaria incidence among children less than 5 years during and after cessation of indoor residual spraying in northern Uganda. Malar J. 2017;16:319.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Ssempiira J, Nambuusi B, Kissa J, Agaba B, Makumbi F, Kasasa S, et al. The contribution of malaria control interventions on spatio-temporal changes of parasitaemia risk in Uganda during 2009–2014. Parasit Vectors. 2017;10:450.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Okech B, Mujuzi G, Ogwal A, Shirai H, Horii T, Egwang TG. High titers of igg antibodies against Plasmodium falciparum serine repeat antigen 5 (sera5) are associated with protection against severe malaria in Ugandan children. Am J Trop Med Hyg. 2006;74:191–7.

    Article  PubMed  CAS  Google Scholar 

  18. Legason ID, Pfeiffer RM, Udquim KI, Bergen AW, Gouveia MH, Kirimunda S, et al. Evaluating the causal link between malaria infection and endemic Burkitt lymphoma in northern Uganda: a Mendelian randomization study. EBioMedicine. 2017;25:58–65.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Simbiri KO, Biddle J, Kinyera T, Were PA, Tenge C, Kawira E, et al. Burkitt lymphoma research in East Africa: highlights from the 9th African organization for research and training in cancer conference held in Durban, South Africa in 2013. Infect Agent Cancer. 2014;9:32.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Uganda Bureau of Statistics. 2002 Uganda population and housing census. Kampala: Government of Uganda; 2007.

    Google Scholar 

  21. Steinhardt LC, Yeka A, Nasr S, Wiegand RE, Rubahika D, Sserwanga A, et al. The effect of indoor residual spraying on malaria and anemia in a high-transmission area of northern Uganda. Am J Trop Med Hyg. 2013;88:855–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Clark TD, Greenhouse B, Njama-Meya D, Nzarubara B, Maiteki-Sebuguzi C, Staedke SG, et al. Factors determining the heterogeneity of malaria incidence in children in Kampala, Uganda. J Infect Dis. 2008;198:393–400.

    Article  PubMed  Google Scholar 

  23. Biggar RJ, Gardiner C, Lennette ET, Collins WE, Nkrumah FK, Henle W. Malaria, sex, and place of residence as factors in antibody response to Epstein–Barr virus in Ghana, West Africa. Lancet. 1981;2:115–8.

    Article  PubMed  CAS  Google Scholar 

  24. Hopkins H, Bebell L, Kambale W, Dokomajilar C, Rosenthal PJ, Dorsey G. Rapid diagnostic tests for malaria at sites of varying transmission intensity in Uganda. J Infect Dis. 2008;197:510–8.

    Article  PubMed  Google Scholar 

  25. Grandesso F, Nabasumba C, Nyehangane D, Page AL, Bastard M, De Smet M, et al. Performance and time to become negative after treatment of three malaria rapid diagnostic tests in low and high malaria transmission settings. Malar J. 2016;15:496.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Koita OA, Doumbo OK, Ouattara A, Tall LK, Konare A, Diakite M, et al. False-negative rapid diagnostic tests for malaria and deletion of the histidine-rich repeat region of the hrp2 gene. Am J Trop Med Hyg. 2012;86:194–8.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Korn EL, Graubard BI. Analysis of health surveys. New York: Wiley; 1999.

    Book  Google Scholar 

  28. Lumley T. Complex surveys: a guide to analysis using R. Hoboken: Wiley; 2010.

    Book  Google Scholar 

  29. Asale A, Duchateau L, Devleesschauwer B, Huisman G, Yewhalaw D. Zooprophylaxis as a control strategy for malaria caused by the vector Anopheles arabiensis (Diptera: Culicidae): a systematic review. Infect Dis Poverty. 2017;6:160.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Binder DA. On the variances of asymptotically normal estimators from complex surveys. Int Stat Rev. 1983;51:279–92.

    Article  Google Scholar 

  31. Rao JNK, Scott AJ. On Chi squared tests for multiway contigency tables with cell proportions estimated from survey data. Ann Stat. 1984;12:46–60.

    Article  Google Scholar 

  32. Crompton PD, Kayala MA, Traore B, Kayentao K, Ongoiba A, Weiss GE, et al. A prospective analysis of the ab response to Plasmodium falciparum before and after a malaria season by protein microarray. Proc Natl Acad Sci USA. 2010;107:6958–63.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Raouf S, Mpimbaza A, Kigozi R, Sserwanga A, Rubahika D, Katamba H, et al. Resurgence of malaria following discontinuation of indoor residual spraying of insecticide in an area of Uganda with previously high-transmission intensity. Clin Infect Dis. 2017;65:453–60.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Hamel MJ, Otieno P, Bayoh N, Kariuki S, Were V, Marwanga D, et al. The combination of indoor residual spraying and insecticide-treated nets provides added protection against malaria compared with insecticide-treated nets alone. Am J Trop Med Hyg. 2011;85:1080–6.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Katureebe A, Zinszer K, Arinaitwe E, Rek J, Kakande E, Charland K, et al. Measures of malaria burden after long-lasting insecticidal net distribution and indoor residual spraying at three sites in Uganda: a prospective observational study. PLoS Med. 2016;13:e1002167.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Weetman D, Wilding CS, Neafsey DE, Muller P, Ochomo E, Isaacs AT, et al. Candidate-gene based GWAS identifies reproducible DNA markers for metabolic pyrethroid resistance from standing genetic variation in east african Anopheles gambiae. Sci Rep. 2018;8:2920.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Glunt KD, Abilio AP, Bassat Q, Bulo H, Gilbert AE, Huijben S, et al. Long-lasting insecticidal nets no longer effectively kill the highly resistant Anopheles funestus of southern Mozambique. Malar J. 2015;14:298.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Anopheles gambiae 1000 Genomes Consortium, et al. Genetic diversity of the African malaria vector Anopheles gambiae. Nature. 2017;552:96–100.

    Google Scholar 

  39. Kiberu VM, Matovu JK, Makumbi F, Kyozira C, Mukooyo E, Wanyenze RK. Strengthening district-based health reporting through the district health management information software system: the Ugandan experience. BMC Med Inform Decis Making. 2014;14:40.

    Article  Google Scholar 

  40. Baik S, Mbaziira M, Williams M, Ogwang MD, Kinyera T, Emmanuel B, et al. A case-control study of Burkitt lymphoma in east Africa: are local health facilities an appropriate source of representative controls? Infect Agent Cancer. 2012;7:5.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Ogwang MD, Bhatia K, Biggar RJ, Mbulaiteye SM. Incidence and geographic distribution of endemic Burkitt lymphoma in northern Uganda revisited. Int J Cancer. 2008;123:2658–63.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Rodriguez-Barraquer I, Arinaitwe E, Jagannathan P, Boyle MJ, Tappero J, Muhindo M, et al. Quantifying heterogeneous malaria exposure and clinical protection in a cohort of Ugandan children. J Infect Dis. 2016;214:1072–80.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Johnston WT, Mutalima N, Sun D, Emmanuel B, Bhatia K, Aka P, et al. Relationship between plasmodium falciparum malaria prevalence, genetic diversity and endemic Burkitt lymphoma in Malawi. Sci Rep. 2014;4:3741.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  44. Emmanuel B, Kawira E, Ogwang MD, Wabinga H, Magatti J, Nkrumah F, et al. African Burkitt lymphoma: age-specific risk and correlations with malaria biomarkers. Am J Trop Med Hyg. 2011;84:397–401.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Donnelly B, Berrang-Ford L, Labbe J, Twesigomwe S, Lwasa S, Namanya DB, et al. Plasmodium falciparum malaria parasitaemia among indigenous Batwa and non-indigenous communities of Kanungu district, Uganda. Malar J. 2016;15:254.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  46. Bulterys PL, Chao A, Dalai SC, Zink MC, Dushimimana A, Katzenstein D, et al. Placental malaria and mother-to-child transmission of human immunodeficiency virus-1 in rural Rwanda. Am J Trop Med Hyg. 2011;85:202–6.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Onori E. Distribution of Plasmodium ovale in the eastern, western and northern regions of Uganda. Bull World Health Organ. 1967;37:665–8.

    PubMed  PubMed Central  CAS  Google Scholar 

  48. White GB. Anopheles gambiae complex and disease transmission in Africa. Trans R Soc Trop Med Hyg. 1974;68:278–301.

    Article  PubMed  CAS  Google Scholar 

Download references

Authors’ contributions

SMM, MDO, PK, SJR, KB, JJG, and RMP designed the study and provided administrative and fieldwork oversight. TK, IO, HN, and IDL conducted and monitored field work. MM and RMP conducted statistical analyses. MM drafted the manuscript. RJB and AOT interpreted data. All authors contributed to the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank the study population for their participation. We thank Ms. Janet Lawler-Heavner at Westat Inc, (Rockville, MD, USA) and Mr. Erisa Sunday at the African Field Epidemiology Network (Kampala, Uganda) for managing the study. We are grateful to Mr. Wilson Nyegenye at Uganda Bureau of Statistics (Kampala, Uganda) for training EMBLEM staff in field survey methods. We thank Ms. Laurie Buck, Dr. Carol Giffen, and Mr. Greg Rydzak at Information Management Services Inc. (Calverton, MD, USA) for coordinating data, and preparing data analysis files.

Competing interests

The authors declare that they have no competing interests.

Availability of data

The data and the code used for these analyses are available on reasonable request from the corresponding authors.

Ethics approval and consent to participate

The Uganda Virus Research Institute Research and Ethics Committee, the Uganda National Council of Science and Technology (HS-816) and the National Cancer Institute Special Studies Institutional Review Board reviewed and approved the study (10-C-N133). Written informed consent was obtained from the parents or guardians of the children and written informed assent was obtained from children aged 7 years or older prior to enrolment.

Funding

The study was funded by the Intramural Research Program of the Division of Cancer Epidemiology and Genetics, National Cancer Institute (NCI) (Contracts HHSN261201100063C and HHSN261201100007I) and, in part, by the Intramural Research Program, National Institute of Allergy and Infectious Diseases (SJR), National Institutes of Health, Department of Health and Human Services. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the US Government. The content of this manuscript is the sole responsibility of the authors.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sam M. Mbulaiteye.

Additional files

Additional file 1: Table S1.

Weighted distribution of characteristics of apparently healthy children aged 0–15 years enrolled between October 2011 and February 2014 in 12 villages in northern Uganda.

Additional file 2: Fig. S1.

Age distribution in single year age groups of the unweighted and weighted population of the children enrolled in 12 randomly selected villages in northern Uganda.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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 Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Maziarz, M., Nabalende, H., Otim, I. et al. A cross-sectional study of asymptomatic Plasmodium falciparum infection burden and risk factors in general population children in 12 villages in northern Uganda. Malar J 17, 240 (2018). https://doi.org/10.1186/s12936-018-2379-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12936-018-2379-1

Keywords