Skip to main content

Age and geographic patterns of Plasmodium falciparum malaria infection in a representative sample of children living in Burkitt lymphoma-endemic areas of northern Uganda



Falciparum malaria is an important risk factor for African Burkitt lymphoma (BL), but few studies have evaluated malaria patterns in healthy BL-age children in populations where both diseases are endemic. To obtain accurate current data, patterns of asymptomatic malaria were investigated in northern Uganda, where BL is endemic.


Between 2011 and 2015, 1150 apparently healthy children under 15 years old were sampled from 100 villages in northern Uganda using a stratified, multi-stage, cluster survey design. Falciparum malaria prevalence (pfPR) was assessed by questionnaire, rapid diagnostic test (RDT) and thick film microscopy (TFM). Weighted pfPR and unadjusted and adjusted associations of prevalence with covariates were calculated using logistic models and survey methods.


Based on 1143 children successfully tested, weighted pfPR was 54.8% by RDT and 43.4% by TFM. RDT sensitivity and specificity were 97.5 and 77.8%, respectively, as compared to TFM, because RDT detect malaria antigens, which persist in peripheral blood after clinical malaria, thus results based on RDT are reported. Weighted pfPR increased from 40% in children aged under 2 years to 61.8% in children aged 6–8 years (odds ratio 2.42, 95% confidence interval (CI) 1.26–4.65), then fell slightly to 49% in those aged 12–15 years. Geometric mean parasite density was 1805.5 parasites/µL (95% CI 1344.6–2424.3) among TFM-positive participants, and it was higher in children aged <5 years at 5092.9/µL (95% CI 2892.7–8966.8) and lower in those aged ≥10 years at 983.8/µL (95% CI 472.7–2047.4; P = 0.001). Weighted pfPR was lower in children residing in sub-regions employing indoor residual spraying (IRS) than in those residing in non-IRS sub-regions (32.8 versus 65.7%; OR 0.26, 95% CI 0.14, 0.46). However, pfPR varied both within IRS (3.2–55.3%) and non-IRS sub-regions (29.8–75.8%; Pheterogeneity <0.001). pfPR was inversely correlated with a child’s mother’s income (P = 0.011) and positively correlated with being enrolled in the wet season (P = 0.076), but sex was irrelevant.


The study observed high but geographically and demographically heterogenous patterns of asymptomatic malaria prevalence among children living in northern Uganda. These results provide important baseline data that will enable precise evaluation of associations between malaria and BL.


Falciparum malaria is an important cause of childhood morbidity and mortality in young children in Equatorial Africa [1, 2]. It is also an important risk factor for African Burkitt lymphoma (BL) [3], a B cell tumour that is the most common cancer in children in many countries in Equatorial Africa [4], including Uganda. The association of falciparum malaria with BL is based on the geographical correlation of high incidence of severe malaria and high incidence of BL [57], and on the significant positive or inverse association of high titres of particular malaria antibodies in children with BL [810]. The epidemiological evidence is supported by evidence from mouse studies, which have shown that infection of mice with Plasmodium chabaudi increases the risk of genetic abnormalities in their B cells, including those similar to the characteristic abnormalities seen in BL involving the translocation of exons in the MYC gene on chromosome 8 into the vicinity of enhancer or promoter regions for immunoglobulin genes on chromosomes 14, 2, or 22 [11, 12].

Whether clinical malaria or the significantly more prevalent asymptomatic parasitaemia similarly affect BL risk is unknown [3]. The overlap between severe malaria and BL by age is imperfect [1, 2] because severe malaria occurs most commonly in young children of around 2 years old [13], whereas BL peaks in children around 6–9 years old [14, 15]. Since BL is an extremely fast-growing tumour [16], the imperfect overlap between the peak ages for severe malaria and BL does not support the idea that severe malaria is the antigenic insult that triggers BL onset. This reasoning is supported by the casual observations that a history of severe malaria is not a commonly elicited feature in children with BL suggesting that older children in the age range of BL are protected against severe malaria by acquired immunity [17, 18] in agreement with findings that the probability of clinical malaria falls by 6% per year of age and by 2% per episode of clinical malaria [19]. Thus, the timing of peak BL risk is consistent with recurrent mild or asymptomatic malaria parasitaemia being the malaria insults that triggers the onset of the genetic abnormalities that lead to BL.

Clarifying the malaria phenotype that is relevant for BL aetiology is an important question that remains unanswered because research efforts are hampered by the lack of reliable data about malaria experience among children of a similar age and sex distribution as BL cases living in BL-endemic areas. To address this knowledge gap, detailed epidemiological data and biological samples were collected from a representative population-based sample of healthy children from northern Uganda with a similar age and sex distribution as typical BL cases. This area has historically experienced holoendemic malaria [2022] (weighted Plasmodium falciparum parasite prevalence (pfPR) over 75%) [22] and high incidence of BL [7]. Baseline data from field-collected data were analysed to obtain a more granular picture of the geographical and demographical patterns of asymptomatic malaria infection in this region. This would pave the way for follow-up studies conducted using molecular and immunological methods to evaluate the malaria response in patients with and without BL and to perform more precise adjustment of the geographical and demographical co-factors.


The aim, design and setting of the study

To obtain geographically and demographically representative data about malaria in northern Uganda, between January 2011 and April 2015, apparently healthy children (0–15 years of age) from the northwest and north-central regions of Uganda were enrolled as part of the EMBLEM study of BL [23]. These two regions were selected because they have historically experienced holoendemic malaria transmission [20] and they experience a correspondingly high, albeit geographically variable, BL endemicity [7]. The annual entomological inoculation rate (AEIR) varies from 397 in  the Arua District in the northwest region to 1586, the highest number recorded worldwide, in  the Apac District in the north-central region [24]. Between 2009 and 2012, the Government of Uganda responded to this high AEIR by implementing indoor residual spraying (IRS) of insecticides and distributing insecticide-treated bed nets to pregnant women and children under 5 years old in 10 districts in the north-central region, hereafter called the IRS sub-region [21]. The non-sprayed area is called the non-IRS sub-region.

Sampling design

A stratified multi-stage cluster sampling design was used to select a regionally representative, population-based sample of healthy children ages 0–15 years old from was selected 100 villages in the study region (Fig. 1). To select the villages, 100 census enumeration areas (EAs) were selected from a complete list of EAs in the study region obtained from the Uganda Bureau of Statistics (UBOS) [25]. This EA sample frame was stratified on ‘low-density’ versus ‘high-density’ population and ‘near water’ versus ‘far from water’ before sampling. Since the stratifying variables are thought to influence mosquito ecology and through that malaria transmission [26, 27], they were considered domains of special interest when estimating geographical patterns of malaria. Low-density EA stratum, considered a surrogate for rural areas, was defined as including those EAs having a population count of children younger than 15 years fewer than the EA mean population count (n = 2683) [25]; otherwise the EAs were categorized as belonging to a high-density EA stratum, considered a surrogate for urban areas. Near water EA stratum was defined as including those EAs with a boundary next to or within 500 m of an all-season surface water body (a swamp, river, lake), based on distances estimated from national maps incorporating geographical information metadata; otherwise the EAs were categorized as belonging to a far from water EA stratum. Because 90% of the Ugandan population resides in rural (or low-density) areas, the high-density EA stratum were oversampled by selecting low-density: high-density EAs in a ratio of 2:1. Thus, 68 EAs were sampled from low-density strata and 32 EAs from high-density with equal probability in each stratum, maintaining a ratio of 1:1 for near water versus far from water EAs within each stratum.

Fig. 1

Flow chart showing the sampling of healthy children aged 0–15 years from 100 villages in north-central and northwest regions of Uganda between January 2011 and April 2015

Because the EAs are large geographical areas that were logistically challenging for research staff to accurately survey, they were segmented into four to 17 units per EA using existing local administrative units called Local Council 1 (LC-1), also known as villages and that have distinct non-overlapping boundaries, as segments. One segment was randomly selected with equal probability per EA. Thus, one segment or village represented the entire EA. In the next stage, 22–25 households per segment were selected with equal probability from a household list constructed by study staff. In the third stage, four to 25 children were selected from a list of all children in the selected households in that segment, based on residing in the study region for at least 4 months and according to the predetermined age and sex frequency distribution of typical BL cases in the region [7].

Participant enrolment

Experienced fieldwork teams visited participants in their homes and invited them to participate in the study after obtaining informed consent. Structured questionnaires were used to collect participant information about the demographics and household characteristics, educational level of a child’s parents, a child’s exposure to malaria suppression methods (use of mosquito bed nets the night before interview and use of IRS in the house), and history of outpatient or inpatient malaria treatment (up to 6 months ago, 7–12 months ago, or 13 or more months ago). A venous blood sample for immediate malaria testing and for storage for future studies was obtained from the participant child using EDTA tubes.

Malaria testing

Malaria infection was diagnosed by experienced laboratory technicians using thick film microscopy (TFM) to identify asexual malaria parasite forms in thick film slides stained with 10% Giemsa solution for 10 min. Malaria parasites were counted against 200 white blood cells (WBCs) and standardized to parasites/µL of blood based on the measured WBC/µL. Thin film smears were examined to identify Plasmodia species. Since ~98% were P. falciparum, this species is assumed hereafter. The technicians also used antigen–antibody capture rapid diagnostic tests (RDT) (MALARIA DUAL kits, ICT Diagnostics, Muizenberg, Cape Town, South Africa) to diagnose malaria infection. These kits detect the P. falciparum-specific malaria histidine-rich protein 2 (Pf-HRP2) and the pan-lactate dehydrogenase (pLDH) antigen shared by other Plasmodia that parasitize humans. The sensitivity and specificity of the kits for malaria in clinical samples in Uganda was reported to be 92–100% [28], which is adequate for the objectives of assessing epidemiological patterns of malaria in the current study. As only 2.1% of the children studied had fever at the time of enrolment, this number was considered negligible and unlikely to impact the overall results so these children were not excluded. The children with fever and positive RDT or TFM were treated for possible clinical malaria with artemisinin-based combination therapy (ACT) as first line drug for treatment of mild clinical malaria following national treatment guidelines. Parents of children without symptoms were advised to seek treatment from their local health centre if a new fever developed.

Data management

Questionnaires were edited in the field and computerized using DataFax and the data submitted via a secure share-portal to Information Management Services, Inc. (IMS, Calverton, MD, USA) to generate analysis files.

Statistical analysis

Descriptive analyses of pfPR, based on RDT, by geographical and demographical characteristics of participants are reported. To obtain regionally representative results, prevalence estimates were calculated as weighted probabilities of selecting the children in the study sample using standard survey statistical packages [29, 30] implemented in R (version 3.2.3). Details of the calculation of sampling weights are described in Additional file 1: Appendix 1 and 2. Thus, the reported estimates are representative for all the children in this region and the variance estimation took the weights into account and accounted for the clustering of the sample of children in the segments (villages). Unadjusted and adjusted odds ratios (ORs), standard errors [31] and Wald-type 95% CIs (95% CIs) [32] of association of pfPR with each variable were calculated. The adjusted models included only those variables with P < 0.10. Because parasite load in asymptomatic children is partly a function of cumulative acquired immunity to malaria, the overall geometric mean parasite density (GMPD) per µL were calculated and compared by sex and age using the Student’s t test statistic. As the objective of these analyses was to describe patterns to obtain granular baseline data for hypothesis generation, the results were not adjusted for multiple testing. Thus, a two-sided P < 0.05 was considered as evidence of statistical significance and values P < 0.10 as suggesting a trend. Analyses of malaria patterns were repeated separately using pfPR defined based on only TFM and on both TFM and RDT and similar results were obtained.


Characteristics of study population

Overall, 2484 households were selected in 100 villages and 2467 of them participated. From these households, 1167 children were invited to participate and 1150 of them were enrolled, representing a weighted sample of 585,762 children in the region. Consistent with the design, to select children with a typical age and gender distribution typical of BL, 44.6% children participating were 6–10 years old and 52.1% were males (Table 1). Consistent with the predominantly rural population of northern Uganda, 67.6% of the children lived in low-density villages and 89.3% lived in villages near water. By region, 60.8% were from the north-central region and 39.2% from the northwest region (the location of these regions is shown in Fig. 2a). With respect to malaria suppression, 32.9% lived in the sub-region employing IRS and 67.1% in the non-IRS sub-region (as defined in Fig. 2b). To further explore malaria variation in sub-regions, the study area was sub-divided into seven sub-regions based on boundaries that correspond roughly to tribal areas (as shown in Fig. 2c). Two of these sub-regions, which contribute to the majority of BL cases in this region [7], contributed 50% of the study population. Some 28.9% of the children reportedly used a mosquito bed net on the night before interview, but 31.3% lived in houses using IRS, which is similar to the per cent living in an IRS sub-region. Regarding socio-economic status, 55.3% of mothers of participants had fewer than 5 years of formal education and 55.7% had a monthly median income below 30,000 Ugandan shillings (equivalent to ~US$10). Regarding history of malaria, 63.5% of the participants reported a history of outpatient malaria treatment, with 55.6% reporting such treatment in the past 12 months. Conversely, only 36.6% reported at least one inpatient treatment for malaria, defined as requiring some observation at a clinic or hospital, with 12.6% reporting such treatment in the past 12 months.

Table 1 Weighted distribution of characteristics of apparently healthy children aged 0–15 years enrolled between January 2011 and April 2015 from the north-central and northwest regions of Uganda
Fig. 2

Map of north-central and northwest regions in Uganda. Showing weighted Plasmodium falciparum parasite prevalence based on the rapid diagnostic test in a representative sample of 1142 healthy children 0–15 years old enrolled between January 2011 and April 2015. Panel a pfPR patterns in the two north-central and northwest regions; Panel b pfPR patterns in two sub-regions employing or not IRS; Panel c pfPR patterns in seven sub-regions. A inset map is included to show the location of north-central and northwest regions in relation to Kampala, the capital city of Uganda, and the location of collaborating hospitals in the study area. It also shows major Ugandan water features (rivers and lakes)

Weighted malaria prevalence and density

Malaria testing was successful for 1143 children; the seven children who failed were excluded from further analysis (Table 2). The weighted pfPR based on RDT was 54.8% (95% CI 47.3–62.4%) and it was 43.4% (36.1–50.7%) based on TFM. As expected, the sensitivity of RDT was 97.5% and the specificity was 77.8%, compared to TFM. RDT is more sensitive than TFM because it detects malaria antigens, which can be detected during early infection and for up to 35–42 days after treatment of symptomatic cases due to persistence of antigens in peripheral blood [33]. Using positivity to both RDT and TFM to define current parasitaemia, the weighted pfPR of current parasitaemia was 42.3%. Conversely, using only RDT positivity to define malaria antigenaemia alone (due to early or persistent antigen after infection has resolved), the weighted pfPR of malaria antigenaemia alone was 12.6%. The weighted prevalence of negative RDT among TFM-positive subjects (including two with parasitaemia identified only on thin film) was 1.1%. These cases (TFM-positive and RDT-negative) had lower GMPD of 178.5 parasites/µL (95% CI 115.0–277.0) than those with of current infection (RDT and TFM positive) with GMPD of 1912.3 parasites/µL (95% CI 1423.0–2570.0; P < 0.001).

Table 2 Weighted Plasmodium falciparum parasite prevalence according to the rapid diagnostic test and thick film microscopy among healthy children in northern Uganda

Geographical and seasonal variation in weighted pfPR

Similar to prior studies from this region [20, 34, 35], there was no difference in the weighted pfPR estimate for the north-central and northwest regions (Fig. 2a). However, these similar malaria patterns masked variation in sub-regional analysis. First, pfPR was a substantially lower in IRS than non-IRS sub-regions (32.8 versus 65.7%; Fig. 2b). However, the apparent homogeneity in IRS and non-IRS sub-regions masked further, substantial variation in sub-regions in the IRS (3.2–55.3%) and the non-IRS sub-regions (55.5–75.8%; Fig. 2c). This additional variation was not explained by patterns in other covariates (Table 3). As expected, the weighted pfPR (Fig. 3a) and the GMPD (Fig. 3b) showed strong seasonal effects, both being higher during the wet season than the dry seasons, but this pattern was more striking in the high-density than in the low-density villages.

Table 3 Plasmodium falciparum parasite prevalence among apparently healthy children 0–15 years old enrolled between January 2011 and April 2015 from the north-central and northwest regions of Uganda
Fig. 3

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

Weighted pfPR and geometric mean parasite density by children’s age

The weighted pfPR increased from 40.0% in children 1–2 years old to 61.8% in those 6–8 years old and decreased slightly to 48.9% in children 12–15 years old (Fig. 4a, b). Using logistic regression to model the weighted pfPR across five age categories that were previously defined in Emmanuel et al. (under 2, 3–5, 6–8, 9–11, 12+ years) [14], the age-specific ORs of pfPR were heterogeneous across these age categories. Compared to children under 2 years old, the ORs for pfPR peaked at 2.42 (95% CI 1.26–4.65) in children 6–8 years old and then decreased slightly to 1.43 (95% CI 0.62–3.31) in children 12–15 years old. There was a significant linear relationship between age and the weighted pfPR (Ptrend = 0.006) and a significant non-linear effect of age on pfPR (Ptrend = 0.004). These contrasting trends are consistent with the idea that the weighted pfPR increases with age as children become increasingly mobile with age and are more likely to be exposed to malaria infection. Conversely, increasing age is associated with increasing cumulative acquired malaria immunity and the chance that older children are more likely to successful resist new infection. Interestingly, these age trends were masked when pfPR across was modeled across three age categories (Table 3, Pheterogeneity = 0.185), highlighting the importance of using finer categories of age to capture rapid changes in the trends. Because BL epidemiology is characterized by male and rural predominance, the age-group specific patterns of pfPR were examined by sex and in high- and low-density villages, but no significant differences in age-specific patterns were noted by sex in the low-density (Fig. 5a) and high-density (Fig. 5b) villages.

Fig. 4

Bar graphs showing age-group patterns of weighted per cent Plasmodium falciparum parasite prevalence. Based on the RDT (a) and the GMPD/µL (b) among apparently healthy, microscopy-positive children enrolled between January 2011 and April 2015 from north-central and northwest Uganda. The unweighted number of participants in each age-group is shown. In b, the solid circle indicates the GMPD, the lines indicate the 95% CIs of the GMPD, except for the under 2 years old age group, where the upper boundary is beyond the plotted points. GMPD results were available on 436 of 441 microscopy positive subjects (see Table 2)

Fig. 5

Bar graphs showing the weighted proportion of children (males and females) positive for Plasmodium falciparum parasite prevalence. Among apparently healthy children enrolled between January 2011 and April 2015 from north-central and northwest Uganda from low population density villages (a) and high population density villages (b). The left y-axis shows the proportion of children by pfPR results. The right y-axis shows the re-weighted population counts of the children represented by the sample in each stratum. The legend and colouring in bars shows pfPR categories based on both the RDT and TFM (see "Methods")

Consistent with their asymptomatic status, GMPD among the 436 of 441 TFM-positive participants with quantitative data was comparatively low at 1805.5 parasites/µL (95% CI 1344.6–2424.3). However, GMPD was significantly and inversely correlated with age. GMPD was higher in children aged <5 versus ≥10 years old (5092.9 parasites/µL, 95% CI 2892.7–8966.8 versus 983.8 parasites/µL, 95% CI 472.7–2047.4; P = 0.001), consistent with the idea that immunity against malaria parasitaemia increases with age among highly exposed children.

Association of pfPR with children’s characteristics

In unadjusted analyses, pfPR was positively associated with enrolment during the wet versus dry season (OR 2.35, 95% CI 1.26–4.40), and inversely associated with living in a sub-region where IRS was implemented versus living in a sub-region where IRS was not implemented (OR 0.26, 95% CI 0.14–0.46). Similar results were obtained when exposure to IRS was assessed using self-reported questionnaire information about whether a child lived or not in a house where IRS was applied in the past year (OR 0.28, 95% CI 0.15–0.53), confirming the validity of the ecological association between IRS and pfPR. The weighted pfPR was inversely associated with a child’s mother’s income ≥30,000 versus <30,000 USHS (OR 0.68, 95% CI 0.47–1.00) and positively associated with a child reporting outpatient treatment for malaria in the past 12 months versus no treatment (OR 1.73, 95% CI 1.1.0–2.72). However, sex, proximity of village to water, population density, region, and mother’s education level, sleeping under a mosquito bed net the previous night, or history of inpatient malaria treatment were unrelated to weighted pfPR (all P > 0.10). In multivariable analyses, including only those variables with P < 0.10 in the unadjusted analyses, the weighted pfPR remained significantly heterogeneous across the seven geographical sub-regions (P = 0.001). In addition, the weighted pfPR was inversely associated with a child’s mother income (OR 0.57, 95% CI 0.37–0.87) and with living in an IRS sub-region (OR 0.37, 95% CI 0.14–1.03), but positively associated with being enrolled during the wet season calendar months (OR 1.81, 95% CI 0.95–3.46). Similar patterns were observed when these analyses were repeated using pfPR defined based on only TFM or combined positivity TFM and RDT results.


The current study is the largest to present accurate granular baseline data about geographical and demographical patterns of asymptomatic malaria infection in healthy children living in a BL-endemic region. As expected, overall malaria prevalence in northern Uganda was high, consistent with the results from previous surveys from that region [20, 34, 35]. In agreement with the earlier studies, the weighted pfPR in the two neighboring regions studied was comparable. However, analysis of sub-regional patterns revealed 50% lower pfPR in sub-regions where IRS was employed compared to those where it was not applied, highlighting IRS as a major determinant of geographical patterns of malaria. However, substantial, significant geographical variation in pfPR was observed within IRS and non-IRS sub-regions, underscoring the contribution from other poorly understood geographical co-factors. This geographical variation, which had been masked by the overall regional-specific pfPR, included the entire spectrum of endemicity ranging from hypoendemic (1–10%), mesoendemic (pfPR 10–50%), hyper-endemic (pfPR 50–74%) to holoendemic (pfPR ≥75%) both within IRS and non-IRS sub-regions (Fig. 2c). The geographical variation in pfPR patterns observed here is likely correlated with the underlying variation in mosquito vector population relative to human hosts—the so called vector-to-host ratio [36], the prevailing vector species, e.g., Anopheles gambiae and Anopheles funestus, and their flying and feeding habits [37, 38]. Since recurrent immunological challenge from parasitaemia or malaria antigenaemia is the malaria exposure of interest in the study of BL, knowledge of these geographical patterns in pfPR can enable more precise evaluation of associations between malaria and BL through matching or adjustment on geography as a proxy for malaria pressure.

In addition to the environmental factors noted above, the geographical variation in the weighted pfPR might be due to genetic variation in prevalent genetic mutations like the sickle cell disease polymorphism that has been linked to malaria resistance [39]. Interestingly, the sub-region observed to have the highest pfPR in this study (sub-region 7; Fig. 1c) was also reported to have the highest prevalence of the sickle cell trait in northern Uganda (18–24%) [40]. This correlation is consistent with published literature suggesting positive selection of the sickle cell trait in areas of high malaria prevalence [41, 42]. This positive correlation could be due to a higher number of secondary malaria infections per infected person (also called the reproductive number of infection) among subjects with the sickle cell trait compared to those without the trait [43]. Specifically, because sickle cell trait carriers typically experience low GMPD infections [44] and low GMPD infections are more likely to persist and support higher rates of gametocyte production and lead to higher transmission rates in populations with a high frequency of the sickle cell trait [45]. Interestingly, Gouagna and colleagues [46] recently reported that parasites also may have adapted to be more efficiently transmitted to mosquitoes by persons with genetic variants in human haemoglobin gene as compared to non-carriers, thereby presenting evidence of parasite adaptation to host genetic mutations. This parasite adaptation to the human genetic mutations could potentially lead to stronger correlations between pfPR and the frequency of sickle cell trait in some areas.

Both pfPR and the GMPD were noted to vary with age, but in different ways. The weighted pfPR increased with age and peaked at ages 6–8 years, then declined slightly. Conversely, GMPD peaked quickly in children aged under 2, 4–6 years before peak pfPR age, then declined rapidly. These contrasting pfPR/GMPD patterns suggest that children are exposed to malaria infection at all ages and infection is usually successfully established until they reach 6–8 years, when the risk of new infection decreases. The GMPD patterns, which decrease rapidly with age, suggest that high parasite counts are experienced mostly at the young ages when children still lack immunity, but infections at older ages occur in the context of immunity to malaria, which suppresses parasitaemia [47]. Given the interest in malaria infection as a risk factor for BL, the observation that GMPD peaks about 4–8 years earlier than BL [14] argues against high GMPD being the insult that triggers the onset of such a rapidly growing tumor [16]. Rather, the occurrence of BL in older children, who are more likely to have lower parasitaemia (Fig. 4b), suggests that recurrent low-count afebrile parasitaemia, probably due to novel allelic variants that cause break-through infections in children with established immunity [48, 49], may be the underlying malaria antigenic insult that triggers the onset of this rapidly growing tumor [50]. If so, examination of immune responses to malaria or allelic variants in patients with BL and well-matched controls could reveal both evidence of strong exposure to malaria, strong protective immune response against malaria, and infection with multi-clonal parasites in children with BL.

The observation of an inverse correlation between pfPR and the use of IRS agrees with results from national malaria indicator surveys conducted in Uganda [21]. Similar to those surveys, the current study also found that mosquito bed nets, which are widely recommended and distributed [51], appear to have had a negligible impact on geographical or individual pfPR patterns in this region. This lack of impact of mosquito bed nets on pfPR in this study, which agrees with reports by others in the region [52], may, in part, be due to misclassification of responses about bed net use or to their inconsistent and/or improper use. However, mosquito bed nets only suppress mosquito populations around bed areas, which may not lead to community-wide suppression of the vector-host ratio like IRS. More worrying, the lack of impact could be an indicator of emerging pyrethroid resistance in this region as has been reported in some studies [53]. If so, the lack of impact of mosquito bed nets on geographical patterns of pfPR should be evaluated more carefully to understand its basis and to exclude more serious reasons such as emerging pyrethroid resistance in regional mosquito populations.

The strengths of the current study include selecting a regionally representative, population-based sample of children to generate accurate regional and sub-regional statistics of pfPR in a BL-endemic area. The use of RDT, which is more sensitive than TFM and provides information about population prevalence of malaria infection over a three- to four-week period [33] is a strength. The finding of similar patterns of pfPR when using results defined by TFM or both TFM and RDT suggests that the patterns reported here are likely valid. Only 1.1% of participants had false negative RDT results, which is similar to prior reports [54]. This percentage is negligible for descriptive studies of geographical and demographical patterns, but the molecular reasons for false negative results, which may be due to deletions of pfhrp2 or pfhrp3 genes that code for the RDT antigens [54], need careful evaluation to preserve the value of these tests for clinical and epidemiological studies. The limitations of this study include relying on recall for information about history of outpatient or inpatient malaria treatment, which is subject to recall bias. Review of health records at health centers was considered to strengthen the results, but the quality of records at health centres in rural Uganda was found to be unreliable. Another limitation is the classification of wet and dry season based on average rain patterns during the study period, hence there might be some misclassification, which would attenuate differences.


The study confirmed high but geographically and demographically heterogeneous malaria prevalence in BL-age children residing in northern Uganda. The finding that the peak age for pfPR and BL are comparable, while GMPD peaks 4–8 years earlier, raises the hypothesis that recurrent asymptomatic malaria exposure in children with established malarial immunity may be the underlying malaria antigenic insult associated with the onset of BL. Understanding the geographical patterns of malaria in this BL-endemic area will pave the way for a more granular examination of the malaria/BL associations in follow-up studies using molecular and immunological methods to measure the malaria response in patients with and without BL.



Plasmodium falciparum prevalence

pfhrp2 :

Plasmodium falciparum histidine-rich protein 2

pfhrp3 :

Plasmodium falciparum histidine-rich protein 3


geometric mean parasite density


odds ratio


confidence interval


annual entomological inoculation rate


Burkitt lymhoma


indoor residual spraying


enumeration area


Uganda Bureau of Statistics


thick film microscopy


while blood cell


rapid diagnostic test


pan lactate dehydrogenase


Ugandan shillings


United States dollar


  1. 1.

    Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526:207–11.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  2. 2.

    Snow RW, Korenromp EL, Gouws E. Pediatric mortality in Africa: Plasmodium falciparum malaria as a cause or risk? Am J Trop Med Hyg. 2004;71:16–24.

    PubMed  Google Scholar 

  3. 3.

    Bouvard V, Baan RA, Grosse Y, Lauby-Secretan B, El Ghissassi F, Benbrahim-Tallaa L, et al. Carcinogenicity of malaria and of some polyomaviruses. Lancet Oncol. 2012;13:339–40.

    Article  PubMed  Google Scholar 

  4. 4.

    Parkin DM, Sitas F, Chirenje M, Stein L, Abratt R, Wabinga H. Part I: cancer in indigenous Africans—burden, distribution, and trends. Lancet Oncol. 2008;9:683–92.

    Article  PubMed  Google Scholar 

  5. 5.

    Morrow RH Jr. Epidemiological evidence for the role of falciparum malaria in the pathogenesis of Burkitt’s lymphoma. IARC Sci Publ. 1985;60:177–86.

    Google Scholar 

  6. 6.

    Rainey JJ, Omenah D, Sumba PO, Moormann AM, Rochford R, Wilson ML. Spatial clustering of endemic Burkitt’s lymphoma in high-risk regions of Kenya. Int J Cancer. 2007;120:121–7.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  8. 8.

    Carpenter LM, Newton R, Casabonne D, Ziegler J, Mbulaiteye S, Mbidde E, et al. Antibodies against malaria and Epstein-Barr virus in childhood Burkitt lymphoma: a case–control study in Uganda. Int J Cancer. 2008;122:1319–23.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Aka P, Vila MC, Jariwala A, Nkrumah F, Emmanuel B, Yagi M, et al. Endemic Burkitt lymphoma is associated with strength and diversity of Plasmodium falciparum malaria stage-specific antigen antibody response. Blood. 2013;122:629–35.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Mutalima N, Molyneux E, Jaffe H, Kamiza S, Borgstein E, Mkandawire N, et al. Associations between Burkitt lymphoma among children in Malawi and infection with HIV, EBV and malaria: results from a case–control study. PLoS ONE. 2008;3:e2505.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Manolov G, Manolova Y. Marker band in one chromosome 14 from Burkitt lymphomas. Nature. 1972;237:33–4.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Robbiani DF, Deroubaix S, Feldhahn N, Oliveira TY, Callen E, Wang Q, et al. Plasmodium infection promotes genomic instability and AID-dependent B cell lymphoma. Cell. 2015;162:727–37.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    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.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    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 

  15. 15.

    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  Google Scholar 

  16. 16.

    Ziegler JL, Bluming AZ, Fass L, Morrow RH Jr, Iversen OH. Burkitt’s lymphoma: cell kinetics, treatment and immunology. Bibl Haematol. 1973;39:1046–52.

    CAS  PubMed  Google Scholar 

  17. 17.

    Goncalves BP, Huang CY, Morrison R, Holte S, Kabyemela E, Prevots DR, et al. Parasite burden and severity of malaria in Tanzanian children. N Engl J Med. 2014;370:1799–808.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Cham GK, Turner L, Lusingu J, Vestergaard L, Mmbando BP, Kurtis JD, et al. Sequential, ordered acquisition of antibodies to Plasmodium falciparum erythrocyte membrane protein 1 domains. J Immunol. 2009;183:3356–63.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    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  Google Scholar 

  20. 20.

    Yeka A, Lameyre V, Afizi K, Fredrick M, Lukwago R, Kamya MR, et al. Efficacy and safety of fixed-dose artesunate-amodiaquine vs. artemether-lumefantrine for repeated treatment of uncomplicated malaria in Ugandan children. PLoS ONE. 2014;9:e113311.

    Article  PubMed  PubMed Central  Google Scholar 

  21. 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.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Metselaar D, Van Thiel PH. Classification of malaria. Trop Geogr Med. 1959;11:157–61.

    Google Scholar 

  23. 23.

    Simbiri KO, Biddle J, Kinyera T, Were PA, Tenge C, Kawira E, et al. Burkitt lymphoma research in East Africa: highlights from the 9(th) 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 

  24. 24.

    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.

    PubMed  Google Scholar 

  25. 25.

    2002 Uganda Population and Housing Census. Kampala: Government of Uganda (Uganda Bureau of Statistics); 2007.

  26. 26.

    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 

  27. 27.

    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.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    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 

  29. 29.

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

    Google Scholar 

  30. 30.

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

    Google Scholar 

  31. 31.

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

    Article  Google Scholar 

  32. 32.

    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 

  33. 33.

    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 

  34. 34.

    Prugger C, Engl M, Ogwang M, Ploner F, Ploner M, Gluderer D, et al. Malariological baseline survey and in vitro antimalarial drug resistance in Gulu district, Northern Uganda. Wien Klin Wochenschr. 2008;120:63–8.

    Article  PubMed  Google Scholar 

  35. 35.

    Kafuko GW, Burkitt DP. Burkitt’s lymphoma and malaria. Int J Cancer. 1970;6:1–9.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Zhou G, Munga S, Minakawa N, Githeko AK, Yan G. Spatial relationship between adult malaria vector abundance and environmental factors in western Kenya highlands. Am J Trop Med Hyg. 2007;77:29–35.

    PubMed  Google Scholar 

  37. 37.

    Lefevre T, Gouagna LC, Dabire KR, Elguero E, Fontenille D, Renaud F, et al. Beyond nature and nurture: phenotypic plasticity in blood-feeding behavior of Anopheles gambiae s.s. when humans are not readily accessible. Am J Trop Med Hyg. 2009;81:1023–9.

    Article  PubMed  Google Scholar 

  38. 38.

    Githeko AK, Adungo NI, Karanja DM, Hawley WA, Vulule JM, Seroney IK, et al. Some observations on the biting behavior of Anopheles gambiae s.s., Anopheles arabiensis, and Anopheles funestus and their implications for malaria control. Exp Parasitol. 1996;82:306–15.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Allison AC. The distribution of the sickle-cell trait in East Africa and elsewhere, and its apparent relationship to the incidence of subtertian malaria. Trans R Soc Trop Med Hyg. 1954;48:312–8.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Ndeezi G, Kiyaga C, Hernandez AG, Munube D, Howard TA, Ssewanyana I, et al. Burden of sickle cell trait and disease in the Uganda Sickle Surveillance Study (US3): a cross-sectional study. Lancet Glob Health. 2016;4:e195–200.

    Article  PubMed  Google Scholar 

  41. 41.

    Aidoo M, Terlouw DJ, Kolczak MS, McElroy PD, ter Kuile FO, Kariuki S, et al. Protective effects of the sickle cell gene against malaria morbidity and mortality. Lancet. 2002;359:1311–2.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Elguero E, Delicat-Loembet LM, Rougeron V, Arnathau C, Roche B, Becquart P, et al. Malaria continues to select for sickle cell trait in Central Africa. Proc Natl Acad Sci USA. 2015;112:7051–4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Smith DL, McKenzie FE. Statics and dynamics of malaria infection in Anopheles mosquitoes. Malar J. 2004;3:13.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Komba AN, Makani J, Sadarangani M, Ajala-Agbo T, Berkley JA, Newton CR, et al. Malaria as a cause of morbidity and mortality in children with homozygous sickle cell disease on the coast of Kenya. Clin Infect Dis. 2009;49:216–22.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Price R, Nosten F, Simpson JA, Luxemburger C, Phaipun L, ter Kuile F, et al. Risk factors for gametocyte carriage in uncomplicated falciparum malaria. Am J Trop Med Hyg. 1999;60:1019–23.

    CAS  PubMed  Google Scholar 

  46. 46.

    Gouagna LC, Bancone G, Yao F, Yameogo B, Dabire KR, Costantini C, et al. Genetic variation in human HBB is associated with Plasmodium falciparum transmission. Nat Genet. 2010;42:328–31.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Teo A, Feng G, Brown GV, Beeson JG, Rogerson SJ. Functional antibodies and protection against blood-stage malaria. Trends Parasitol. 2016;32:887–98.

    CAS  Article  PubMed  Google Scholar 

  48. 48.

    Newbold CI, Schryer M, Boyle DB, McBride JS, McLean A, Wilson RJ, et al. A possible molecular basis for strain specific immunity to malaria. Mol Biochem Parasitol. 1984;11:337–47.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Smith T, Beck HP, Kitua A, Mwankusye S, Felger I, Fraser-Hurt N, et al. Age dependence of the multiplicity of Plasmodium falciparum infections and of other malariological indices in an area of high endemicity. Trans R Soc Trop Med Hyg. 1999;93(Suppl 1):15–20.

    Article  PubMed  Google Scholar 

  50. 50.

    Rochford R, Cannon MJ, Moormann AM. Endemic Burkitt’s lymphoma: a polymicrobial disease? Nat Rev Microbiol. 2005;3:182–7.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Snow RW, Lengeler C, de Savigny D, Cattani J. Insecticide-treated bed nets in control of malaria in Africa. Lancet. 1995;345:1056–7.

    CAS  Article  PubMed  Google Scholar 

  52. 52.

    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 

  53. 53.

    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 

  54. 54.

    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.

    CAS  Article  PubMed  PubMed Central  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, IDL, SK, RJB, AWB, BE, and MJ conducted and monitored field work. AWB provided technical advice. MS, MM and RMP conducted statistical analyses. MM drafted the manuscript. All authors contributed to the manuscript. All authors read and approved the final manuscript.

Authors’ information

SMM is a medical graduate from Makerere University, Kampala, in Uganda with training in epidemiology and biostatistics from the University of Cambridge. His career started in Uganda as an HIV and cancer doctor, but he transferred to the Infections and Immunoepidemiology Branch (IIB) in the Division of Cancer Epidemiology and Genetics (DCEG) at the national Cancer Institute (NCI), where he is a senior investigator. His research interests are Burkitt lymphoma and Kaposi sarcoma and he conducts epidemiological studies to characterize risk factors for these conditions. Both malignancies are endemic in Africa and their risk is substantially increased in the setting of HIV/AIDS. He is currently the principal investigator for Epidemiology of Burkitt Lymphoma in East-African Children and Minors (EMBLEM) study, a population-based case–control study being conducted in 6 regions in three countries of East Africa (Uganda, Kenya, and Tanzania) to investigate infectious and genetic risk factors of BL. This study aims to create a rich granular data and sample repository  for discovery  of risk factors for BL in East Africa.


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, Mr. Greg Rydzak, and Mr. Jeremy Lyman at Information Management Services Inc. (Calverton, MD, USA) for coordinating data, preparing data analysis files and drawing the maps. We thank Mr. David Check and Ms. Adrienne Rolls of Biostatistics Branch at NCI (Bethesda, MD, USA) for their help formatting the Tables. We thank Mr. Graham Kalton at Westat Inc, (Rockville, MD, USA) for comments on the manuscript.

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.


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.

Author information



Corresponding authors

Correspondence to Ruth M. Pfeiffer or Sam M. Mbulaiteye.

Additional file

Appendices I and II.

Additional file 1.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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 ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Maziarz, M., Kinyera, T., Otim, I. 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 16, 124 (2017).

Download citation


  • Burkitt lymphoma
  • Africa
  • Plasmodium falciparum
  • Malaria
  • Epidemiology
  • Non-Hodgkin lymphoma
  • Uganda