Plasmodium falciparum malaria parasitaemia among indigenous Batwa and non-indigenous communities of Kanungu district, Uganda

Background The indigenous Batwa of southwestern Uganda are among the most highly impoverished populations in Uganda, yet there is negligible research on the prevalence of malaria in this population. Plasmodium falciparum malaria parasitaemia prevalence was estimated in an indigenous Batwa and a non-indigenous neighbouring population, and an exploration of modifiable risk factors was carried out to identify potential entry points for intervention. Additionally, evidence of zooprophylaxis was assessed, hypothesizing that livestock ownership may play a role in malaria risk. Methods Two cross-sectional surveys of Batwa and non-Batwa communities were carried out in Kanungu District, Uganda in July 2013 and April 2014 based on a census of adult Batwa and a two-stage systematic random sample of adult non-Batwa in ten Local Councils where Batwa settlements are located. A community-based questionnaire and antigen rapid diagnostic test for P. falciparum were carried out in the cross-sectional health surveys. A multivariable logistic regression model was built to identify risk factors associated with positive malaria diagnostic test. A subset analysis of livestock owners tested for zooprophylaxis. Results Batwa experienced higher prevalence of malaria parasitaemia than non-Batwa (9.35 versus 4.45 %, respectively) with over twice the odds of infection (OR 2.21, 95 % CI 1.23–3.98). Extreme poverty (OR 1.96, 95 % CI 0.98–3.94) and having an iron sheet roof (OR 2.54, 95 % CI 0.96–6.72) increased the odds of infection in both Batwa and non-Batwa. Controlling for ethnicity, wealth, and bed net ownership, keeping animals inside the home at night decreased the odds of parasitaemia among livestock owners (OR 0.29, 95 % CI 0.09–0.94). Conclusion A health disparity exists between indigenous Batwa and non-indigenous community members with Batwa having higher prevalence of malaria relative to non-Batwa. Poverty was associated with increased odds of malaria infection for both groups. Findings suggest that open eaves and gaps in housing materials associated with iron sheet roofing represent a modifiable risk factor for malaria, and may facilitate mosquito house entry; larger sample sizes will be required to confirm this finding. Evidence for possible zooprophylaxis was observed among livestock owners in this population for those who sheltered animals inside the home at night.


Background
There are renewed calls for malaria eradication with a focus on Africa [1,2]. Malaria mortality rates are decreasing in many populations, with global incidence having fallen by approximately 37 % since 2000 to 214 million new cases in 2015 [3,4]. Despite these gains, malaria remains a major global disease burden, with approximately 438,000 deaths annually [3,5].
The elimination of malaria from many western nations has been attributed primarily to social and economic development allowing for screening of windows and doors, destruction of vector breeding sites, and rapid diagnosis and treatment [6]. The feasibility of Plasmodium falciparum malaria elimination in most of sub-Saharan Africa is low, with Uganda being among the countries with lowest feasibility [7]. Sub-Saharan African countries are disproportionately affected by malaria [4] due to the presence of highly competent mosquito vectors, widespread poverty, limited infrastructure, and overburdened health systems [8,9]. Those living in extreme poverty are most vulnerable to infectious diseases, yet within-country disparities are often ignored [10,11]. African indigenous populations in particular have consistently poorer health outcomes than their nonindigenous counterparts [12]. Social determinants of indigenous health include, but are not limited to, poverty, discrimination, limited access to health care, and loss of traditional lands [13]. Indigenous and ethnic minority populations outside of Uganda experience higher rates of malaria, which have been attributed to relative impoverishment, marginalization, and geographic remoteness [14][15][16][17][18]. In some cases, genetic variations have been identified as drivers of ethnic differences in malaria parasitaemia and immunological response [19].
Risk factors for malaria can be conceptualized as nonmodifiable and modifiable. Non-modifiable factors such as age and sex have been inconsistently associated with higher malaria risk. In endemic areas, children under 5 years of age have high risk of malaria due to their immunological naïveté [10,20,21]. Beyond this, the relationship between age and malaria in adults is less clear [22]. It is known that regular exposure to malaria results in a functional immunity to the disease which quickly wanes in the absence of exposure. Excluding immunologically-suppressed pregnant women, sex-related variations in malaria risk are generally linked to gender-role and occupational exposures [23,24]. Women in roles as household water-collectors may spend more time near mosquito breeding sites. Men with forest-related jobs may spend their time in mosquito-dense areas [10,23].
Modifiable risk factors include environmental conditions and human behaviour. Local vector ecology, including the locations of swamps, forests, rice paddies where vectors breed, and the proximity of these sites to human habitations may bring humans into frequent contact with mosquitoes [25][26][27]. Housing conditions may also affect transmission. Open eaves and windows, for example, may permit mosquito entry into sleeping quarters [25,28,29]. Education and wealth are known to be protective against malaria. Understanding malaria transmission and prevention may result in behaviours such as staying indoors during peak vector activity, and having access to preventive strategies that may reduce exposure of humans to vectors [27,30]. Many of these risk factors are mediated by poverty whereby access to building materials and bed nets may be dependent upon income [31][32][33]. Poverty has been described as a key modifiable determinant of malaria burden [6].
Livestock are routinely included in analyses of risk factors for malaria infection, yet their role in malaria transmission is not well understood. Livestock are dead-end hosts for human-infectious Plasmodium parasites and may reduce human malaria risk by drawing vectors away from humans (zooprophylaxis) [34][35][36]. In some cases, however, livestock act as additional blood meal sources, and they may alter vector longevity and population density to increase human malaria risk (zoopotentiation) [37,38]. Empirical research exploring the association between livestock and malaria risk has been complicated by the confounding role of wealth [39]. Livestock symbolize social standing, provide food and services, and act as an asset to be sold in times of financial need [40][41][42][43][44]. Both livestock and wealth are generally associated with lower malaria prevalence. Given efforts to include zooprophylaxis in integrated vector control programs, there have been calls for further research to understand the impact of livestock on malaria transmission [39].
The Ugandan Batwa are an indigenous population with life expectancy and child mortality rates significantly worse than national averages [45]. There remain considerable gaps regarding our understanding of Batwa health; to our knowledge only two peer-reviewed studies currently report the prevalence and risk factors of health outcomes for Batwa in Uganda [23,46]. While literature on Batwa livelihoods is currently limited, it is thought that Batwa engage in livestock livelihoods at a much lower rate than non-Batwa due to financial restrictions, in addition to their historical hunter-gatherer culture. Further, Namanya identified consistent differences across risk factors for malaria, including housing and education between Batwa and non-Batwa [47].
The aims of this paper were (1) to estimate malaria prevalence in indigenous Batwa of Kanungu District and compare this with their non-indigenous neighbours, (2) explore modifiable risk factors for malaria parasitaemia in order to identify potential entry points for intervention to reduce malaria prevalence among Batwa and non-Batwa and, (3) to test the hypothesis of zooprophylaxis in a sub-set of livestock owners from both Batwa and non-Batwa populations.

Study population
Kanungu District is located in southwestern Uganda bordering the Democratic Republic of Congo to the west. It contains part of the Bwindi Impenetrable National Park (Bwindi) [48] and has a population of approximately 250,000, 80 % of whom live in rural settlements [49]. The majority of the population are of Bakiga ethnicity [45,50]. The non-Batwa (non-indigenous populations, including Bakiga) rely mainly on subsistence farming consisting of cash-and food-cropping, and small-scale livestock holdings [48,50], but tourism centred around gorilla trekking in Bwindi also provides local employment [45]. The indigenous Batwa population of Kanungu District is approximately 900 [51]. When Bwindi received its national park designation in 1991, the Batwa were no longer permitted to enter the park. Thus, they lost their ancestral lands and livelihood as traditional hunter-gatherers. Since then, the Batwa have been transitioning to settled agriculture [45]. Kanungu District is remote and has limited infrastructure and service delivery [50]. It receives rainfall throughout the year, with two dryer seasons from December to February and from June to July. Malaria in Kanungu is characterized by low transmission intensity and low endemicity [8,23,52,53].

Study design and sample
Two cross-sectional, in-person surveys were administered in the 10 Local Councils (LCs, the smallest government administration unit in Uganda) containing the 10 Batwa settlements in July 2013 and April 2014 (Fig. 1). These 10 LCs are located within the subcounties of Kayantorogo, Kayonza, Kirima, Mpungu and Butogota Town Council. Considering the small total population size of the Kanungu Batwa population, a census of all adult Batwa present (>18 years) was attempted. A two-step proportional systematic random-sample of adult non-Batwa was conducted, representing approximately 40 % of non-Batwa households in each LC. Sampling frames (a list of households in the LC) were provided by Batwa settlement and non-Batwa Local Council 1 (LC1) chairpersons. The sampling strategy consisted of the following: a number between 1 and 10 was randomly selected as the first household to be selected from the sampling frame. The total number of households was multiplied by 0.4 to obtain the sampling interval required to obtain a 40 % sample of the LC. Every nth household was then selected from the sampling frame. For the second stage, the LC1 chairperson conducted a census of all adults living in the selected households from stage 1. The LC1 chairperson selected the adult to be sampled out of a hat based on the census of adults in those households. The same individuals were interviewed in July and April.

Data collection
Plasmodium falciparum is the primary cause of malaria in this region [8]. Parasitaemia (outcome variable) was classified as a positive result on CareStart ™ malaria rapid diagnostic test (RDT) detecting histidine-rich protein 2 (HRP2) of P. falciparum (Access Bio, Inc. NJ, USA) at either the July or April data collections [54]. The RDT has a sensitivity of 98 % and specificity of 97.5 % [55]. Having a positive test result at either collection period was defined as "RDT ever" while a negative outcome at both collection periods was defined as "RDT never".
At each survey (July and April), participants were interviewed in Rukiga (the local language for Batwa and non-Batwa) by trained local survey administrators. The survey consisted of three parts: (1) an individual-level health questionnaire, (2) a household questionnaire and (3) a malaria RDT administered by a community health nurse. Part one (individual questionnaire) collected information pertaining to demographics and personal information, health and risk behaviour, and livelihoods. Part two (household questionnaire) collected information on household structure, housing materials, assets, water, and sanitation. For Batwa, each individual completed parts one and three, with the household head completing part two. The individual questionnaires were later matched to the household questionnaires. For non-Batwa, each individual completed all components. Data were collected in hard copy and entered into Microsoft Access (Microsoft Corp, USA).

Ethics
This study was designed using a community-based participatory research approach and was approved by the Research Ethics Boards at McGill University and the University of Guelph. The study design is consistent with the Canadian Tri-Council's Policies and requirements of the Ethical Conduct of Research Involving Human Subjects. Informed consent was obtained verbally from all participants. This study is a part of the Indigenous Health Adaptation to Climate Change (IHACC) project, an international initiative with parallel field sites in the Canadian Arctic and Peruvian Amazon [56].

Malaria parasitaemia prevalence and risk factors
Among Batwa and non-Batwa a total of 649 questionnaires were carried out in July 2013 and 546 surveys were completed in April 2014. All questionnaires from July and April were pooled to create a single study population. This was done to: (1) increase the sample size, (2) achieve a prevalence reflective across seasons, and (3) account for the potential sporadic nature of infection and risk factors using only one point estimate in time. For individuals who were present for both the July and April surveys, and who were RDT negative at both times (n = 709), the risk factor data from the July surveys were used as baseline.
If an individual was present in July and April and was RDT positive at one time (n = 49), data from the time at which they tested positive were used. If they were present at both July and April and positive at both times (n = 1), the July data were used. The period prevalence was calculated by dividing the total number of positive cases by the total number of participants tested.
To identify risk factors associated with parasitaemia, a series of univariable logistic regression models were constructed. Risk factors of interest were selected based on theorized relevance in the literature and data availability. Covariates included: age, sex, knowledge of malaria transmission (as reflected by knowing that one should avoid mosquito bites to prevent malaria), roofing material, bed net ownership, and wealth (Table 1). To measure wealth, a household-level asset-based wealth index was created using principal component analysis [46,57,58]. The variables included in the index were cell phone ownership, radio ownership, bicycle ownership, having a source of electricity, land ownership, receiving remittances, toilet type, and having access to hand washing facilities. The resultant component was categorized as a dichotomous variable representing the 50 % least poor  Livestock have been associated with increases in malaria prevalence through zoopotentiation [25,38], and decreases in prevalence through zooprophylaxis [76,77,83] and wealth [26,33] Control  [24,28,29] Animals kept inside at night Question: If you own animals, do your animals come into the house during the night?
Dichotomous dummy variable (yes/no) Zooprophylaxis/zoopotentiation may be determined by relative proximity of animals to human sleeping quarters [29,[86][87][88] relative to the 50 % poorest individuals. Livestock were not included in the wealth index to allow the evaluation of their specific role in malaria transmission. Collinearity between risk factor variables was assessed using a Spearman's rank correlation coefficient with a cut-off of 60 %. Where collinearity was identified, the variable with the strongest association with the outcome variable was retained. A multivariable logistic regression model was built, and a final, reduced multivariable model based on variables significant in either the univariable or multivariable model. A mixed-effects logistic regression model with random intercept was tested, but there was no evidence of clustering at the LC-level. Data were analysed using Stata v.11 (Stata Corp., USA).

Zooprophylaxis
A subset analysis of livestock owners assessed the impact of relative location of livestock to humans. The subset was achieved by selecting from the survey dataset all Batwa and non-Batwa individuals who had a positive outcome on the "any livestock" variable. Among those who owned any livestock, the risk factor of interest was whether livestock were housed inside or outside of the home at night. This analysis controlled for bed net use and non-livestock wealth, which have been proposed as potential confounding and effect modifying variables, respectively [39]. Sensitivity analyses were used to assess the effect of housing specific livestock species (goats and chickens) inside the house on parasitaemia. Although age and sex were not significant in the reduced multivariable model, sensitivity analysis also tested their potential effect on the final model.

Batwa and non-Batwa demographics
A total of 759 individual questionnaires were completed in July 2013 and April 2014, of which 59.10 % (449) of participants were non-Batwa and 40.90 % (310) Batwa (Table 2). Females were overrepresented in the study population compared to males, slightly more so among non-Batwa than among Batwa. Men were more likely than women to be employed or away finding work and therefore unavailable to participate.

Malaria parasitaemia prevalence and risk factors
Of the 759 questionnaires completed in July 2013 and April 2014, 50 participants tested positive for P. falciparum antigen on RDT. One participant tested positive at both data collection events. Only the July data were considered for this participant, resulting in a total of 49 positive cases. Of the positive cases, 29 ( There was no association between LC and parasitaemia, nor was there a significant difference between models that controlled for clustering at the LC-level and those that did not. There was minimal difference in the fit of the model with and without random intercepts. The model results were not sensitive to survey date, and showed consistent results on the removal of ethnicity and/or wealth.

Zooprophylaxis
Amongst the 339 Batwa and non-Batwa livestock owners within the total survey population of 759 respondents, 17 individuals had a positive RDT outcome. Of the 339 livestock owners, 99 (29.2 %) were Batwa and 240 (70.8 %) were non-Batwa. Keeping animals inside the house at night was significantly protective against malaria (OR 0.29, CI 0.09-0.94) in adjusted analyses ( Table 4). The sample size was insufficient to disaggregate by livestock species (chickens and goats). Sensitivity analyses indicated consistency of results across species, as well as no effect on the results when adding age and/or sex in the model.

Discussion
The period prevalence of P. falciparum parasitaemia for July 2013 and April 2014 for adult Batwa was 9.35 %, which is higher than the 4.1 % rate previously reported for Batwa over the age of 5 years [23]. A bed net distribution was carried out in November 2012, which could explain the low prevalence rate found previously in this population in January 2013 [59]. These findings are consistent with parasite prevalence estimates for Kanungu District as being relatively low (10-40 % prevalence when including children) compared to other parts of Uganda [8,60]. Ethnicity was associated with malaria parasitaemia, with Batwa at higher risk. This is an important result given that wealth and other theorized risk factors for parasitaemia were controlled for; there remains a residual effect of ethnicity beyond these factors. There are marked differences in livelihoods between the two populations that may be driving these results. Batwa are currently undergoing a transition as they adapt to life outside of the forest. This transition is reflected in the low levels of education and livestock ownership relative to non-Batwa. There may also be unmeasured genetic differences Table 3 Univariable and multivariable explanatory    between these ethnic groups that influence malaria prevalence [61][62][63][64]. For instance, ethnic differences in immune responses [19,65] and lower susceptibility to malaria infection among ethnic minority Fulani tribes relative to sympatric ethnic groups in West Africa have been attributed to genetic differences at various loci [66,67]. It might be expected, given their higher prevalence of infection, that Batwa mount a weaker immune response relative to non-Batwa in the face of exposure. Immunological studies would be required to test this hypothesis, but it is unclear how the results of such studies would lead to intervention strategies [6]. These results provide tentative indication that there may be important modifiable risks for malaria infection among Batwa and non-Batwa. Iron sheet roofing may play a role in malaria risk. House construction is an important risk factor for malaria infection [68][69][70] and thatched houses have previously been related to higher malaria prevalence of inhabitants in other parts of sub-Saharan Africa [28] due to the propensity of mosquitoes to rest indoors [25]. However, a systematic review and meta-analysis on the effect of house construction on malaria found that modern roof materials (iron sheets, tiles) were not consistently associated with decreased odds of infection [71]. Thatched or wooden roofs may confer some protection to individuals relative to corrugated iron sheets. It is possible that the open eaves of iron sheet roofs facilitate mosquito house entry and thereby increased malaria risk. Open eaves and gaps in housing materials are frequently associated with higher rates of parasitaemia [24,29,71]. However, this result may be confounded by wealth. Among non-Batwa, iron sheet roofs are a reflection of wealth since they must be purchased at a high cost compared to thatch or banana fibre, which may be obtained at no cost from the surrounding environment. In contrast, iron sheet roofing is purchased for Batwa by a local NGO with priority given to the most impoverished families as identified by the community. A Ugandan study carried out in part in Kanungu District found that closed eaves and modern house construction were associated with significantly decreased human biting rates and incidence rates of malaria in children when controlling for age, sex, and household wealth [72]. Further evaluation of the impact of house construction on malaria risk among Batwa may help to inform and improve these community development strategies.
The analysis of livestock owners suggested that when controlling for wealth and bed net use, keeping animals inside at night reduced the odds of malaria infection. Previous studies suggest that keeping animals indoors at night increases malaria risk where mosquitoes are zoophilic [39]. Entomological studies in the region suggest that Anopheles funestus and Anopheles gambiae sensu lato are the predominant vector species [60]. Anopheles funestus tends towards anthropophily; however, An. gambiae s.l. consists of seven morphologically indistinguishable species [73] of which the most important vectors may be more zoophilic, as with Anopheles arabiensis, or anthropophilic, as with An. gambiae sensu stricto. This finding might suggest that the predominant local species are anthropophilic and, in turn, that keeping animals indoors at night results in reduced malaria risk among livestock owners, however further entomological evaluation of local vector populations is required. The mechanism through which livestock infer protection against malaria remains poorly understood and their role in malaria transmission has been much contested [32,35,37,38,74,75]. Livestock may draw mosquitoes away from humans, reducing their exposure to malaria [76,77], or may provide an abundance of blood-meals, increasing vector density and longevity [32,38,74]. The zooprophylatic effect is, however, complicated by the relationship between livestock and wealth [32,39]. It is well recognized that livestock contribute significantly to household economies for the rural poor, including the most marginalized [40,41,43]. Wealthy community members in Kanungu are those with the greatest livestock and land holdings. The findings of this study are consistent with others showing that poverty is positively associated with malaria prevalence [31,78,79].
The cross-sectional nature of this study limits the ability to infer causal relationships between the risk factors and malaria infection. Unemployed people who were available to be surveyed during data collection events were overrepresented in the sample. This may have reduced variation within the population and perhaps led to an underestimation of the importance of employment for malaria infection. Similarly, men were more likely to be working away, leading to an overrepresentation of women. Given that the role of sex in malaria risk remains unclear, it is difficult to predict the direction of this effect. The small sample size and number of cases limited statistical power to detect effect sizes. Sample size also prevented the stratification of uni-and multi-variable analyses of parasitemia based on ethnicity, although sensitivity analysis suggested results were consistent between ethnic groups. Similarly, sample size prevented ethnicity-based stratification for the zooprophylaxis analysis. As a result, the identification of risk factors, and recommendations for malaria control, apply to the survey population as a whole. In highly vulnerable remote indigenous populations, sample sizes are typically small and demographics frequently differ from non-indigenous populations; lack of statistical power must be balanced with the importance of prioritizing research within vulnerable and at-risk sub-populations.

Conclusions
Indigenous Batwa in Kanungu District, Uganda experience a two-fold increase in malaria risk compared to their non-indigenous neighbours. This inequitable burden mirrors health disparities experienced by indigenous peoples worldwide. Tentative support for the role of housing construction and wealth in accounting for some of the risk differential was found. Investigation of roof construction and related vector entry may be a prudent follow-up measure given these results. High baseline poverty among both populations, but particularly among Batwa, will remain a major determinant of health inequity and transmission risk. The influence of livestock may be related to asset-based wealth and suggests that keeping livestock indoors at night may play a role in reducing human exposure to malaria in this setting. This evidence may support livestock-based interventions here as both a poverty reduction strategy as well as a component of malaria vector control.