Mapping the baseline prevalence of lymphatic filariasis across Nigeria
Obiora A. Eneanya1,
Claudio Fronterre2,
Ifeoma Anagbogu3,
Chukwu Okoronkwo3,
Tini Garske1,
Jorge Cano2 na1 &
Christl A. Donnelly1,4 na1
Parasites & Vectors volume 12, Article number: 440 (2019) Cite this article
The baseline endemicity profile of lymphatic filariasis (LF) is a key benchmark for planning control programmes, monitoring their impact on transmission and assessing the feasibility of achieving elimination. Presented in this work is the modelled serological and parasitological prevalence of LF prior to the scale-up of mass drug administration (MDA) in Nigeria using a machine learning based approach.
LF prevalence data generated by the Nigeria Lymphatic Filariasis Control Programme during country-wide mapping surveys conducted between 2000 and 2013 were used to build the models. The dataset comprised of 1103 community-level surveys based on the detection of filarial antigenemia using rapid immunochromatographic card tests (ICT) and 184 prevalence surveys testing for the presence of microfilaria (Mf) in blood. Using a suite of climate and environmental continuous gridded variables and compiled site-level prevalence data, a quantile regression forest (QRF) model was fitted for both antigenemia and microfilaraemia LF prevalence. Model predictions were projected across a continuous 5 × 5 km gridded map of Nigeria. The number of individuals potentially infected by LF prior to MDA interventions was subsequently estimated.
Maps presented predict a heterogeneous distribution of LF antigenemia and microfilaraemia in Nigeria. The North-Central, North-West, and South-East regions displayed the highest predicted LF seroprevalence, whereas predicted Mf prevalence was highest in the southern regions. Overall, 8.7 million and 3.3 million infections were predicted for ICT and Mf, respectively.
QRF is a machine learning-based algorithm capable of handling high-dimensional data and fitting complex relationships between response and predictor variables. Our models provide a benchmark through which the progress of ongoing LF control efforts can be monitored.
LF is thought to be endemic in large parts of Nigeria [1]. Although endemicity mapping is mostly completed at district level nationwide [2], some districts in the North-East remain unmapped mainly due to security issues. For mapped areas, there is a need to further understand intra-district heterogeneity in prevalence. Modelling has shown that the success of control programmes to interrupt LF transmission highly relied upon the intensity of transmission prior to the scale-up of MDA interventions [3]. Therefore, knowing the spatial heterogeneity in the intensity of infection throughout implementation areas, as opposed to simply endemic/non-endemic classification, would enable control programmes to identify areas which may require enhanced interventions as they approach the endgame in the elimination pathway. Furthermore, producing maps based on a modelling approach serves as a tool to validate endemicity maps that are already in use for control programmes and potentially identify any discrepancies in endemicity classifications.
Control of LF is largely based on MDA interventions, whereby entire endemic populations are treated with repeated rounds of antifilarial medications. Prior to the implementation of this interventions it was necessary to determine the endemicity status of the geographical areas to be treated [4]. Mapping surveys relied upon the detection of circulating filarial antigens in blood samples of adults from selected communities using a rapid ICT card [5]. However, prior to the existence of ICT tests, LF diagnosis was based on the detection of circulating filarial worm, Mf, by microscopic examination of thick blood smears [6]. These tools are key to LF control programmes because they inform decisions regarding endemicity status, allow for monitoring control interventions and ultimately provide the necessary evidence of interruption of infection transmission.
Although Rapid Assessment of the Geographical Distribution of Bancroftian Filariasis (RAGFIL) surveys have been used to generate initial estimates of the burden of LF [7, 8], these estimates have been improved upon by the use of geostatistics. Geostatistical modelling has enabled the prediction of infection prevalence in unsampled locations across large geographical areas using a suite of potential disease drivers such as climate, environmental and demographical data [9,10,11,12,13,14,15]. Models that take into account the spatial structure of the infection distribution are commonly used both in frequentist [16] and Bayesian [17] modelling frameworks for prevalence mapping. The importance of accounting for spatial effects in prevalence models has been previously explored [18,19,20], and methods for handling spatially correlated data have been suggested [16, 17]. Generalised linear models with spatially correlated random effects, otherwise known as generalised linear spatial models (GLSMs), are widely used to fit binomial data with spatial structure [18]. Despite this being the method of choice for modelling prevalence data for a variety of diseases including schistosomiasis [21], LF [10] and malaria [22], a major drawback of GLSMs is their limitations to handle high-dimensional, non-linear and collinear predictors and response datasets [23]. Machine learning based algorithms have proven to be powerful tools to handle complex relationships between continuous and binary data and independent covariates [24, 25], but little is known about their performance when it comes to modelling binomial data obtained through randomised surveys [26].
In this study, a model was trained based on baseline prevalence data collected through mapping surveys conducted across Nigeria and a suite of environmental and demographic data using a machine learning algorithm, Quantile Regression Forest (QRF). Then, the trained model was used to predict the prevalence and related uncertainty for unsampled locations based on the selected predictors.
Lymphatic filariasis data
Community-level prevalence data (both ICT and Mf) collected during nationwide mapping surveys conducted by the Nigeria Lymphatic Filariasis Control Programme from 2000–2013 was used in this analysis. Also, historical data, mostly surveys based on parasitological diagnosis (Mf detection), publicly available and assembled by the Global Atlas of Helminth Infections [27] were included. All surveys were conducted prior to the implementation of MDA interventions. Up to two communities were surveyed by local government area (LGA) during nationwide mapping. The LGA corresponds to the second administrative level for Nigeria and is considered the geographical area for the implementation of control interventions, thus called the implementation unit (IU).
Within each IU at least one sample village was randomly selected for survey and a buffer zone of at least 50 km separated pairs of sample villages. LF endemicity was estimated by testing for filarial antigenemia in peripheral blood using rapid ICT. For sample villages recording ≥ 1% prevalence, the entire IU is considered endemic for LF. IUs which have frequent reports of hydrocele and lymphedema are strongly suspected to have high endemicity for LF. Such villages are thus identified as sentinel sites for evaluating control programmes and, in addition to ICT card test, LF burden is estimated by Mf load in venous blood collected during night-time surveys. The timing of blood collection coincides with the appearance of Mf in blood, known as nocturnal periodicity. The sampling methodology for the LF surveys has been described in greater detail elsewhere [28]. The mapping methodology described above, however, does not apply to urban LF infections mainly due to the differences in LF vectors and living conditions between urban and rural areas [29].
For this analysis, we had 1103 ICT and 184 Mf prevalence estimates (Fig. 1). These were all pre-intervention data testing 142,881 survey participants, 15 years or older. In building the model, we utilized both ICT and Mf observed prevalence, creating a new covariate 'diagnostic type' (an indicator variable) to assess the effects of diagnostic method on the spatial distribution of LF. Mean observed prevalence for ICT and Mf were 9.5% and 5.0%, respectively.
Location of study sites in Nigeria showing the prevalence of lymphatic filariasis. a Immunochromatographic test (ICT). b Microfilaria (Mf)
The spatial correlation in the observed LF prevalence dataset was measured using a variogram analysis. The variogram gives a measure of the variability between pairs of datapoints [18]. This is important as it was used to assess the degree of spatial autocorrelation that remains on the residuals of the fitted models. Figure 1 shows the distribution of survey locations for ICT and Mf respectively.
Climate and environmental predictors
Continuous gridded maps of climate, topography, vegetation and land use for Nigeria were obtained from different sources. All variables considered had biologically plausible association to LF prevalence. In total, 17 environmental variables were considered in the construction of this model. Precipitation and temperature variables processed from the WorldClim database were all long-term (1950–2000) averages of data from weather stations distributed across the world [30]. A detailed description of these predictors and their processing is provided in our previous work [28]. All input grids raster covariates were resampled to a common spatial resolution of 1 × 1 km2 using the nearest-neighbour algorithm [31].
Quantile regression forest models
QRF is an extension of Random Forest (RF) and is adept at handling large, complex datasets [32]. RF algorithm is an ensemble learning method for classification and regression based on the construction of regression decision trees. This machine learning-based algorithm has proven to outperform other approaches under similar modelling problems [33, 34]. Briefly, trees are grown through recursive partitioning into binary splits from a primary root node which contains all data. For each split, a random subset of predictor variables (approx. one third) is used to grow new root nodes. Each partition contains a different random bootstrapped sample (approx. two thirds) of the dataset. Using bootstrapped samples avoids the problem of overfitting in RF models [35]. This process is repeated until a terminal node is reached, and the average of all the trees is used to make predictions. The response variables not selected during binary node splits, known as the 'out-of-bag' cases, are used to evaluate the predictive accuracy of the model and generating estimating the variable importance.
In RF models (for regression), only the mean of the bootstrapped response variable is considered when splitting/growing trees and for quantile determination [34]. All other features of the response variable of possible interest are neglected. QRF was thus developed to consider all the values in the response variable for splitting and quantile determination [32, 36, 37]. Therefore, QRF enables the estimation of any quantile from the entire posterior conditional distribution for a modelled outcome. Accounting for all the features of the response variable is thought to give a more complete picture of the dataset and resulting predictions [38]. The RF and QRF algorithm uses bagging to randomly resample the training dataset (with replacement of original data) and builds a forest of trees, whereas boosted regression trees use boosting to randomly resample training dataset (without replacement) and builds a sequence of trees with each added tree focussing on poorly fitted nodes.
In this work, an RF model is initially fitted to tune parameters for use in the QRF. Here a 10-fold internal cross-validation was performed and repeated 5 times on empirical logit-transformed infection prevalence and set of predictors. The empirical logit-transformation of infection prevalence was weighted by number of individuals examined using the following formula:
$${\tilde{Y}_{i} = \log \left( {\frac{{Y_{i} + \frac{1}{2}}}{{m_{i} - Y_{i} + \frac{1}{2}}}} \right)}{:} \quad i = 1, \ldots , n$$
where Yi is the observed number of people infected at location i and \(m_{i}\) is the number of people examined.
This tuning process informs an optimum number of predictor variables to be considered at each node split
Using the optimal number of predictors yielded above, a QRF model was then constructed. For the QRF model, data were partitioned into two, with a random subset of 25% of the complete dataset retained for model validation and the remaining 75% used to train the model. The mean, median and prediction intervals estimates were obtained and projected over a continuous geographical space at a spatial resolution of 5 × 5 km. The RF and QRF models were implemented using the randomForest [39] and quantregForest [32] packages, respectively, in R (v.3.3.2) [40]. Predictive maps were exported into ArcGIS v10.3 for preparing map layouts [41].
Model evaluation was performed using the validation dataset based on the root mean squared error (RMSE) and R-squared scores (R2). Variable importance was represented by percentage increase in mean square error (%IncMSE). The %IncMSE is estimated with out-of-bag cross-validation as a result of a variable being permuted (values randomly shuffled). The difference between the calculated mean square errors is then averaged over all trees and then normalised by the standard deviation of the differences [42]. If a predictor is important in the model, then assigning other values for that predictor randomly should have a negative influence on prediction, thereby resulting in a higher %IncMSE value. To check for the existence of spatial structure on the data not captured by the predictors after implementing the QRF model, an empirical variogram was calculated based on the residuals of the QRF model. The Pearson's correlation coefficient was calculated between pairs of observed and predicted ICT and Mf values. Results were presented as 95% prediction intervals and P-values.
From the final predicted prevalence maps, the number of people infected with LF was calculated. By overlaying the prevalence predictions on 2010 gridded population density estimates obtained from the WorldPop Africa repository [43], on a cell-by-cell basis, estimates for the infected population in each cell was generated. Population density data available for Nigeria were for the years 2006, 2010, 2015 and 2020. As the survey data used for building the models spanned from 2000–2013, we estimated the population infected based on population density estimates of 2010. All cells were summed up to get estimates for mainland Nigeria. This analysis was calculated using the Zonal Statistics function within the Spatial Analyst Tool in ArcGIS v10.3 [41].
Analysis for the models was performed using a total of 1287 site-level infection prevalence surveys for ICT (1103 surveys) and Mf (184 surveys) tests respectively as shown in Fig. 1.
Variogram analysis
The results of the variogram analysis (Fig. 2) indicate that there is significant spatial correlation in the observed ICT prevalence data. The range of spatial correlation is about 250 km after which the points start flattening out, an indication of the limit of spatial correlation between datapoints. Conversely, for Mf prevalence, there is limited evidence of spatial correlation, even at shorter distances.
Variogram plot showing the spatial correlation in observed LF data. a Immunochromatographic test (ICT). b Microfilaria (Mf). The empirical variogram is represented by black dots; the theoretical variogram is represented by a solid line
Variable importance plot
Figure 3 shows the variable importance plot of the QRF model trained using LF prevalence data. Here, %IncMSE shows that diagnostic type, precipitation in the driest and wettest quarter, distance to permanent water bodies and land surface temperature were the 5 most important predictors for constructing our model.
Variable importance for the trained model analysing both immunochromatographic test and microfilaria data simultaneously
Predicted ICT and Mf prevalence
Predicted prevalence estimates were projected on the map of mainland Nigeria based on a suite of climate and environmental predictors at a spatial resolution of 5 × 5 km. RMSE and R-square values for the model were 1.24 and 0.40, respectively.
The maps presented in Fig. 4 is the predicted median and upper and lower bounds of ICT prevalence of LF in Nigeria. The median ICT map suggests a high prevalence of LF primarily in three out of six regions in Nigeria, namely North-West, North-Central and South-East of Nigeria. Also, transmission appears to be high in the South-West state of Ekiti.
Median and 95% prediction interval of the immunochromatographic test (ICT) prevalence prediction
The median Mf maps in Fig. 5 shows a distinct spatial pattern of LF, predicting a higher Mf prevalence in much of the southern region of Nigeria and along the course of Niger and Benue rivers. All regions (except the North-West) had predicted prevalence exceeding 1% for Mf.
Median and 95% prediction interval of the microfilaria (Mf) prevalence prediction
Validation of predictive models
The variogram fitted on the residuals demonstrates that the trained QRF model, despite this not being a spatially explicit model, was able to capture, through some of the spatially varying predictors, the spatial structure in the observed ICT and Mf prevalence (Fig. 6). When exploring the correlation between observed and predicted ICT and Mf prevalence, there was a significant positive correlation: Pearson's coefficient of 0.63 (95% CI: 0.57–0.67) and 0.51 (95% CI: 0.37–0.60) between observed and fitted values for ICT and Mf prevalence respectively.
Empirical variogram of the residuals resulting from trained Quantile Regression Forest model
Also, cross-validation of the predicted LF prevalence using a subsample of 25% of the observed data was performed (Fig. 7). Predictive intervals (shadow area) have been centred and observed prevalence for held-out subsample plot on it. 75.2% of the surveys fall within the prediction intervals. In Additional file 1: Figure S1, predicted prevalence values are plotted against observed prevalence. Those plotted in red had observations outside of the 95% prediction intervals while those plotted in blue had observations inside the 95% prediction intervals. Many of the observations plotted in red had 0% observed prevalence.
Cross-validation of the predicted lymphatic filariasis prevalence using a subsample of 25% of the observed data
Estimating population infected with lymphatic filariasis
The mean human population infected with LF is estimated to be 8.7 million and 3.3 million for ICT and Mf respectively (Table 1). This amounts to a national prevalence of 5.3% for ICT and 2.0% for Mf. Total national population for Nigeria for 2010 was derived from gridded population density estimates from the WorldPop repository [43, 44].
Table 1 Estimated number of people infected with lymphatic filariasis prior to MDA using 2010 population estimates
Maps produced in this analysis are intended to estimate the prevalence of LF in unsampled locations, highlight intra-district heterogeneity of infection, and estimate the population infected with LF. These would help guide programme activities for a more focussed intervention. Here, our results are particularly helpful in classifying LGAs in Borno State which are yet to be mapped under the national mapping survey.
The predicted prevalence levels of LF ICT and Mf presented here demonstrate that LF distribution in Nigeria is largely ubiquitous. For ICT, prevalence estimates were highest in the North-Central (8.2%), North-West (7.8%), and South-East (7.1%) regions, and lowest in the South-South (2.5%) region. In contrast, prevalence estimates for Mf were generally higher in the southern than in the northern regions. Overall, predicted mean national prevalence was 5.3% and 2.0% for ICT and Mf, respectively. As MDA for LF is ongoing in Nigeria, robust estimates of baseline prevalence are important for evaluating the efficiency of control efforts.
Models provided state-level LF prevalence predictions for Nigeria, delineating within-region heterogeneities in infection prevalence. The ICT prevalence estimates generated from this analysis are in keeping with earlier work describing the environmental suitability of LF in Nigeria (Additional file 2: Figure S2) [28]. Both maps also correspond well with the distribution of Anopheles spp., which is the chief vector for LF transmission in Nigeria [29] and endemicity maps of the national LF control programme [2]. The climate and environmental variables that contributed the most to predicting the prevalence of LF in unsampled locations were precipitation, land surface temperature, and distance to permanent water bodies. The influence of these variables on the mosquito vectors and their effect in driving the transmission of LF has been discussed in our previously published work [28]. Machine learning methods have been widely used to model distribution of various parasitic diseases, both for species classification [45] or for regression analysis [34]. Their ability to handle non-linear associations between response and predictor variables, control for interactions among predictor variables and handle large complex datasets is a major advantage of these methods [34].
In this work, the prevalence of LF in Nigeria was modelled using the QRF algorithm which is an extension of the RF. The models were constructed by combining infection prevalence data from both ICT and Mf diagnostic types; however, final maps were projected according to diagnostic type. The main reason for the uneven availability of ICT and Mf data is that the ICT-based survey is considered the method of choice for the mapping of LF prior to intervention, and has been used for this purpose since 2000. In contrast, Mf surveys were conducted only in areas suspected to be highly endemic for LF, using lymphedema and hydrocele cases as an indicator for high endemicity. This left large portions of the country without Mf survey points, with some states with no survey points at all, while others were very sparsely surveyed. A visual observation of the Mf survey plot shows more dense clustering of the surveys in the southern parts of the country. This selective sampling of sites may have biased the Mf survey locations as hard-to-reach and more rural areas are more likely to be ignored. Also, as blood testing is performed at night (between 10 pm and 2 am) to coincide with the nocturnal periodicity of the parasite in blood, this cumbersome approach may be a contributing factor to a biased and selective survey.
The ICT and Mf prevalence predictions presented in this work are distinctly different. It is well known, however, that estimates for ICT are generally higher than Mf estimates even in surveys conducted in similar locations [5]. Understanding this contrasting prevalence values is an ongoing challenge in LF research and previous works have attempted to model the relationship between ICT and Mf prevalence [46, 47]. Irvine et al. [46] demonstrated that ICT and Mf prevalence is a consequence of the distribution of adult worms and the subsequent microfilariae production, although it is suggested that ICT prevalence is relatively uninformative in providing estimates of the infective pool [46]. This is mainly due to the therapeutic action of MDA and the nonlinear relationship between adult worm burden and Mf output [46]. Mf prevalence provides a more accurate estimate of microfilariae worm load and thus, a good tool for measuring infectious pool within LF endemic communities.
Furthermore, the treatment regimen used for MDA is known to be more effective against the microfilariae and less so for the adult worms [48]. Therefore, following treatment, ICT prevalence tend to decline more slowly due to the continued presence of the adult worms and their production of the filarial antigen which is still detected in blood using ICT cards [48]. Mf prevalence declines much more quickly due to the stronger microfilaricidal effects of treatment. In addition, in Nigeria, where there are vast areas co-endemic for both LF and onchocerciasis [49], and of which MDA for onchcocerciasis, using ivermectin, pre-dates the survey data used in this study. Although this study models the baseline prevalence of LF, it does not account for the therapeutic effects of onchocerciasis treatment in LF co-endemic areas [50].
Although ICT and Mf prevalence maps are contrasting, this is not to say that one is better or that predictions from the other should be disregarded. Overall, the main determinant factor in the QRF model for the LF prevalence maps is the diagnostic type (Fig. 3), indicating that this accounts for most of the variability in distribution. Understandably, ICT and Mf prevalence maps appear quite different from each other. However, as mentioned earlier, the main goal of the ICT surveys is to better understand the general geographical extent of infection transmission [7], while the Mf surveys, due to the rapid microfilaricidal effects of treatment, are useful for providing a more precise numerical measure with which to evaluate control interventions and track control progress.
The predictive accuracy of models is usually tested by exploring the ability of the model to correctly predict on an independent dataset [51]. As there was no independent dataset to test the model on, predictive accuracy was evaluated by calibrating the QRF model on a random sample of 75%, and then predicting on the held-out 25%. Predictive accuracy was measured by the value of the R2, which is the percentage of variation explained by the climate and environmental covariates included in the model. With an R2 value of 40%, more than half of the variation in the model predictions is not explained by factors included in our analysis. In this work, effects of the malaria control programmes (particularly bednet usage) and its influence on LF prevalence in areas of co-endemicity were not accounted for. This is perhaps a factor to consider going forward, as malaria and LF endemic areas largely overlap [52] and bednets used for malaria prevention have been widespread in Nigeria [53] and have been demonstrated to be protective against LF [54]. Further, with the increase in gross domestic product of Nigeria by more that 10-fold (from $46.4 billion in 2000 to $514.96 billion in 2013) [55], it is believed that general living conditions are improved, providing better protection against mosquito vectors [56]. A larger proportion of the population living above the poverty line and better awareness of the aetiology of LF and malaria may also have led to personal protection measures from the mosquito vectors.
Although the random forest algorithm is growing in popularity for use for spatial predictions, it fails to account for residual spatial correlations in observations [57], however, the inclusion of corresponding geographical coordinates of the survey dataset as a predictor could address this issue. The existence of spatial autocorrelation on the cross-validation residuals is an indication of suboptimal model predictions [57]. Our results (Fig. 6) suggest no evidence of spatial autocorrelation in the cross-validation residuals.
After more than five rounds of MDA, the North-Central states of Plateau and Nasarawa have demonstrated evidence of interruption of transmission and in 2017 transmission assessment survey commenced [58]. For large portions of the country that MDA is currently ongoing, findings from this work will aid re-assessment of programme activities. For instance, in the allocation of preventive chemotherapy and making sure that number of treatments offered are enough to achieve the stipulated population and programme coverage. Furthermore, there is a risk of resurgence in areas where transmission has been interrupted mainly due to the continued presence of mosquito vectors and within-country human migration. Areas previously identified as highly endemic for LF will be key in monitoring prevalence levels going forward. Additionally, periodic entomological examination of mosquito vectors for the presence of the filarial antigen (xenomonitoring) is an effective tool to determine whether the parasite is still present in populations where transmission had been interrupted [59]. This should be used in combination with Mf and ICT surveys, though caution must be exercised when evaluating control measures by testing with ICT because filarial antigenemia is still detected in blood samples long after MDA has been completed [6].
Maps are presented in this work have relatively wide prediction intervals. Given these wide intervals, predictions are in keeping with previous knowledge of LF endemicity in Nigeria [60]. However, there are several ways to narrow intervals in QRF models. It has been suggested that decreasing the spatial resolution reduces uncertainty [38]. Lower resolution maps may be more useful for predicting infections at a larger geographical scale (for instance, continent-wide or on a global scale) as is intended to give estimates for regions rather than smaller areas within a country. Here maps have been projected at a 5 × 5 km spatial resolution in order to better delineate intra-state prevalence levels. Maps produced at lower resolution will be a lot coarser and risk losing their intended purpose of estimating prevalence levels at the lowest administrative level. However, mapped estimates of the mean prevalence estimates for ICT closely mirror maps describing the ecological niche of LF in Nigeria (Additional file 2: Figure S2) [28] and it is believed these prevalence estimates are a fair picture of the LF distribution for Nigeria.
The human population infected with LF in Nigeria was previously estimated to be 13.53 million [10]. This estimate, however, was derived from modelling 27 infection data points of which the most recent survey was conducted in the year 1990 using population estimates of the year 2000. These surveys were mostly conducted by individual researchers in their region of interest with no inter-survey standardisation. The geographical extent and standardised protocol of the survey dataset used for the present analysis coupled with a robust modelling method provides the most comprehensive picture of LF endemicity in Nigeria to date. We however did not account for age structure in our population estimates; therefore, figures presented in this work are likely to be higher than the actual infected population.
Mathematical models have been used to assess the impact of intervention (mainly MDA and vector control) on the LF transmission [61]. The three most commonly used models for evaluating LF interventions are the population-based model, EPIFIL [62], and individual-based models TRANSFIL and LYMFASIM [63, 64]. These models have been trialled on data from LF endemic communities in India, Papua New Guinea and Kenya [61]. One key parameter in these models is baseline LF prevalence levels. These prevalence levels, however, were obtained from surveys conducted by in-country programmes which were sometimes patchy and did not have complete geographical coverage of the area of interest. Further, these data follow the structure of RAGFIL mapping, where entire districts are classed as either endemic or non-endemic according to results from a single survey point within the district. Estimates from our work, however, provide a comprehensive picture of the baseline LF burden for Nigeria. This will be invaluable data to parameterise these models in different settings and to assess the spatial heterogeneity of control efforts.
Modelling the baseline endemicity of infection should ideally be performed prior to scaling up control programmes. Good knowledge of the extent of disease burden is also useful for raising awareness and serves as framework for advocacy for community/institutional engagement. Since the LF control programme is already ongoing in Nigeria, these model estimates provide a basis with which to evaluate control efforts and encourage more coordination towards reaching the elimination targets. Prevalence estimates provided may also serve as a proxy for estimating the burden LF morbidity (lymphedema and hydrocele) [65] in line with planning morbidity management and disability prevention programmes.
The datasets used for this work are available in the Expanded Special Project for Elimination of Neglected Tropical Diseases (ESPEN) repository http://espen.afro.who.int/countries/nigeria. predictor datasets and R code are available upon reasonable request to the corresponding author.
GLSM:
generalised linear spatial model
ICT:
immunochromatographic test
IU:
implementation unit
LF:
LGA:
local government area
MDA:
microfilaria
QRF:
quantile regression forest
RAGFIL:
Rapid Assessment of the Geographical Distribution of Bancroftian Filariasis
Federal Ministry of Health Nigeria. Neglected Tropical Diseases Nigeria Multi-Year Master Plan 2015–2020. Abuja: Federal Ministry of Health Nigeria; 2016.
Expanded special project for elimination of neglected tropical diseases W-A. ESPEN - NTD portal Brazzaville, Republic of Congo; 2017. http://espen.afro.who.int/countries/nigeria. Accessed Jan 2019.
Michael E, Malecela-Lazaro MN, Kabali C, Snow LC, Kazura JW. Mathematical models and lymphatic filariasis control: endpoints and optimal interventions. Trends Parasitol. 2006;22:226–33.
World Health Organization. Monitoring and Epidemiological Assessment of Mass Drug Administration for the Global Programme to Eliminate Lymphatic Filariasis (GPELF). A Manual for National Elimination Programmes WHO/HTM/NTD/PCT/2011.4. Geneva: World Health Organization; 2011.
Weil GJ, Lammie PJ, Weiss N. The ICT filariasis test: a rapid-format antigen test for diagnosis of bancroftian filariasis. Parasitol Today. 1997;13:401–4.
Weil GJ, Ramzy RMR. Diagnostic tools for filariasis elimination programs. Trends Parasitol. 2007;23:78–82.
Gyapong JO, Remme JHF. The use of grid sampling methodology for rapid assessment of the distribution of bancroftian filariasis. Trans Roy Soc Trop Med Hyg. 2001;95:681–6.
World Health Organization. Operational guidelines for rapid mapping of bancroftian filariasis in Africa (WHO/CDS/CPE/CEE/2000.9). Geneva: World Health Organization; 2000.
O'Hanlon SJ, Slater HC, Cheke RA, Boatin BA, Coffeng LE, Pion SDS, et al. Model-based geostatistical mapping of the prevalence of Onchocerca volvulus in West Africa. PLoS Negl Trop Dis. 2016;10:e0004328.
Slater H, Michael E. Mapping, Bayesian geostatistical analysis and spatial prediction of lymphatic filariasis prevalence in Africa. PloS One. 2013;8:e71574.
Moraga P, Cano J, Baggaley RF, Gyapong JO, Njenga SM, Nikolay B, et al. Modelling the distribution and transmission intensity of lymphatic filariasis in sub-Saharan Africa prior to scaling up interventions: integrated use of geostatistical and mathematical modelling. Parasit Vectors. 2015;8:221–37.
Adigun AB, Gajere EN, Oresanya O, Penelope V. Malaria risk in Nigeria: Bayesian geostatistical modelling of 2010 malaria indicator survey data. Malar J. 2015;14:77–85.
Ekpo UF, Hurlimann E, Schur N, Oluwole AS, Abe EM, Mafe MA, et al. Mapping and prediction of schistosomiasis in Nigeria using compiled survey data and Bayesian geospatial modelling. Geospatial Health. 2013;7:355–66.
Deribe K, Cano J, Newport MJ, Golding N, Pullan RL, Sime H, et al. Mapping and modelling the geographical distribution and environmental limits of podoconiosis in Ethiopia. PLoS Negl Trop Dis. 2015;9:e0003946.
Deribe K, Cano J, Giorgi E, Pigott D, Golding N, Pullan R, et al. Estimating the number of cases of podoconiosis in Ethiopia using geostatistical methods Wellcome Open Res. 2017;2:78.
Carsten FD, Jana MM, Miguel BA, Roger B, Janine B, Gudrun C, et al. Methods to account for spatial autocorrelation in the analysis of species distributional data: a review. Ecography. 2007;30:609–28.
Lunn DJ, Thomas A, Best N, Spiegelhalter D. WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comput. 2000;10:325–37.
Diggle PJ, Tawn JA, Moyeed RA. Model-based geostatistics. J R Stat Soc Ser C. 1998;47:299–350.
Paterson S, Lello J. Mixed models: getting the best use of parasitological data. Trends Parasitol. 2003;19:370–5.
Raso G, Matthys B, N'Goran EK, Tanner M, Vounatsou P, Utzinger J. Spatial risk prediction and mapping of Schistosoma mansoni infections among schoolchildren living in western Côte d'Ivoire. Parasitology. 2005;131:97–108.
Clements AC, Lwambo N, Blair L, Nyandindi U, Kaatano G, Kinung'hi S, et al. Bayesian spatial analysis and disease mapping: tools to enhance planning and implementation of a schistosomiasis control programme in Tanzania. Trop Med Int Health. 2006;11:490–503.
Kleinschmidt I, Bagayoko M, Clarke GPY, Craig M, Le Sueur D. A spatial statistical approach to malaria mapping. Int J Epidemiol. 2000;29:355–61.
Cade Brian S, Noon Barry R. A gentle introduction to quantile regression for ecologists. Front Ecol Environ. 2003;1:412–20.
Singh A, Thakur N, Sharma A, editors. A review of supervised machine learning algorithms. In: 3rd International conference on computing for sustainable global development (INDIACom), 16–18 March 2016, New Delhi. India; 2016.
Dey A. Machine learning algorithms: a review. Int J Comput Sci Inf Technol. 2016;7:1174–9.
Guo P, Liu T, Zhang Q, Wang L, Xiao J, Zhang Q, et al. Developing a dengue forecast model using machine learning: a case study in China. PLoS Negl Trop Dis. 2017;11:e0005973.
Cano J, Rebollo MP, Golding N, Pullan RL, Crellen T, Soler A, et al. The global distribution and transmission limits of lymphatic filariasis: past and present. Parasit Vectors. 2014;7:101–21.
Eneanya OA, Cano J, Dorigatti I, Anagbogu I, Okoronkwo C, Garske T, et al. Environmental suitability for lymphatic filariasis in Nigeria. Parasit Vectors. 2018;11:513–26.
Okorie PN, McKenzie FE, Ademowo OG, Bockarie M, Kelly-Hope L. Nigeria Anopheles vector database: an overview of 100 years' research. PLoS One. 2011;6:e28347.
WorldClim. Free climate data for ecological modeling and GIS. 2017. http://worldclim.com/. Accessed Jan 2019.
Yates D, Gangopadhyay S, Rajagopalan B, Strzepek K. A technique for generating regional climate scenarios using a nearest-neighbor algorithm. Water Resour Res. 2003;39:1199–214.
Meinshausen N. Quantile regression forests. J Mach Learn Res. 2006;7:983–99.
Breiman L. Random forests. Mach Learn. 2001;45:5–32.
Roy M-H, Larocque D. Robustness of random forests for regression. J Nonparametr Stat. 2012;24:1–14.
Swatantran A, Dubayah R, Goetz S, Hofton M, Betts MG, Sun M, et al. Mapping migratory bird prevalence using remote sensing data fusion. PLoS One. 2012;7:e28922.
Bandreddy N, Devabhaktuni V, Alam M, and Kumar A. Estimation of unmeasured radon concentrations in Ohio using quantile regression forest. ProQuest Dissertations and Thesis. Gurgaon, India: ProQuest Dissertations Publishing; 2014.
Vaysse K, Lagacherie P. Using quantile regression forest to estimate uncertainty of digital soil mapping products. Geoderma. 2017;291:55–64.
Kampichler C, Sierdsema H. On the usefulness of prediction intervals for local species distribution model forecasts. Ecol Inf. 2017;47:67–72.
Liaw A, Wiener M. Classification and Regression by randomForest. R News. 2002;2:18–22.
R Developement Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.
Esri - ArcGIS 10.3. Geographic information system software. 2017.
Forkuor G, Hounkpatin OKL, Welp G, Thiel M. High resolution mapping of soil properties using remote sensing variables in south-western Burkina Faso: a comparison of machine learning and multiple linear regression models. PLoS One. 2017;12:e0170478.
Tatem AJ. WorldPop, open data for spatial demography. Sci Data. 2017;4:170004.
Linard C, Gilbert M, Snow RW, Noor AM, Tatem AJ. Population distribution, settlement patterns and accessibility across Africa in 2010. PLoS One. 2012;7:e31743.
Marmion M, Parviainen M, Luoto M, Heikkinen RK, Thuiller W. Evaluation of consensus methods in predictive species distribution modelling. Divers Distrib. 2009;15:59–69.
Irvine MA, Njenga SM, Gunawardena S, Njeri Wamae C, Cano J, Brooker SJ, et al. Understanding the relationship between prevalence of microfilariae and antigenemia using a model of lymphatic filariasis infection. Tran R Soc Trop Med Hyg. 2016;110:118–24.
Cano J, Moraga P, Nikolaya B, Rebollob MP, Okorie PN, Davies E, et al. An investigation of the disparity in estimates of microfilaraemia and antigenemia in lymphatic filariasis surveys. Tran R Soc Trop Med Hyg. 2015;109:529–31.
Taylor MJ, Hoerauf A, Bockarie M. Lymphatic filariasis and onchocerciasis. Lancet. 2010;376:1175–85.
Okorie PN, Ademowo GO, Saka Y, Davies E, Okoronkwo C, Bockarie MJ, et al. Lymphatic filariasis in Nigeria; micro-stratification overlap mapping (MOM) as a prerequisite for cost-effective resource utilization in control and surveillance. PLoS Negl Trop Dis. 2013;7:e2416.
Brant TA, Okorie PN, Ogunmola O, Ojeyode NB, Fatunade SB, Davies E, et al. Integrated risk mapping and landscape characterisation of lymphatic filariasis and loiasis in South West Nigeria. Parasit Epidemiol Control. 2018;3:21–35.
Allouche O, Tsoar A, Kadmon R. Assessing the accuracy of species distribution models: prevalence, kappa and the true skill statistic (TSS). J Appl Ecol. 2006;43:1223–32.
Adigun AB, Gajere EN, Oresanya O, Vounatsou Penelope. Malaria risk in Nigeria: Bayesian geostatistical modelling of 2010 malaria indicator survey data. Malar J. 2015;14:156.
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.
Kelly-Hope LA, Molyneux DH, Bockarie MJ. Can malaria vector control accelerate the interruption of lymphatic filariasis transmission in Africa; capturing a window of opportunity? Parasit Vectors. 2013;6:52–64.
The World Bank Group. World Bank Open Data Washington DC, USA; 2018. https://data.worldbank.org/country/nigeria?view=chart.
Kaindoa EW, Finda M, Kiplagat J, Mkandawile G, Nyoni A, Coetzee M, et al. Housing gaps, mosquitoes and public viewpoints: a mixed methods assessment of relationships between house characteristics, malaria vector biting risk and community perspectives in rural Tanzania. Malar J. 2018;17:298–314.
Hengl T, Nussbaum M, Wright MN, Heuvelink GBM, Gräler B. Random forest as a generic framework for predictive modeling of spatial and spatio-temporal variables. PeerJ. 2018;6:e5518.
The Cater Centre. Two states in Nigeria eliminate disfiguring parasitic disease lymphatic filariasis as public health problem; 2017. https://www.cartercenter.org/news/pr/nigeria-101317.html. Accessed Jan 2019.
Irish SR, Al-Amin HM, Paulin HN, Mahmood ASMS, Khan RK, Muraduzzaman AKM, et al. Molecular xenomonitoring for Wuchereria bancrofti in Culex quinquefasciatus in two districts in Bangladesh supports transmission assessment survey findings. PLoS Negl Trop Dis. 2018;12:e0006574.
Expanded Special Projects for Elimination of Neglected Tropical Diseases (ESPEN). Status of Lymphatic filariasis MDA (2005–2016)—Nigeria WHO; 2018. http://espen.afro.who.int/system/files/content/maps/WHO_LF_IU_MDA_TC_trend_NG.pdf.
Smith ME, Singh BK, Irvine MA, Stolk WA, Subramanian S, Hollingsworth TD, et al. Predicting lymphatic filariasis transmission and elimination dynamics using a multi-model ensemble framework. Epidemics. 2017;18:16–28.
Norman RA, Chan MS, Srividya A, Pani SP, Ramaiah KD, Vanamail P, et al. EPIFIL: the development of an age-structured model for describing the transmission dynamics and control of lymphatic filariasis. Epidemiol Infect. 2000;124:529–41.
Irvine MA, Hollingsworth TD. Making transmission models accessible to end-users: the example of TRANSFIL. PLoS Negl Trop Dis. 2017;11:e0005206.
Stolk WA, De Vlas SJ, Borsboom GJ, Habbema JD. LYMFASIM, a simulation model for predicting the impact of lymphatic filariasis control: quantification for African villages. Parasitology. 2008;135(Special Issue 13):1583–98.
Gyapong JO, Webber RH, Morris J, Bennett S. Prevalence of hydrocele as a rapid diagnostic index for lymphatic filariasis. Tran R Soc Trop Med Hyg. 1998;92:40–3.
We are grateful to the Ministry of Health staff in Nigeria who conducted data collection and to all participants who gave their time generously during the surveys.
The Commonwealth Scholarship Commission (NGCA-2015-051) provided funding for the doctoral studentship of OAE. JC and CF are funded by grants from the Bill and Melinda Gates Foundation (OPP 1033751). We acknowledge joint Centre funding from the UK Medical Research Council and Department for International Development – grant reference: MR/R015600/1. The funders played no role in the study design, data collection, formal analysis, interpretation of results and decision to publish this study.
Jorge Cano and Christl A. Donnelly contributed equally to this work
MRC Centre for Global Infectious Disease Analysis, Department of Infectious Disease Epidemiology, Imperial College London, London, UK
Obiora A. Eneanya
, Tini Garske
& Christl A. Donnelly
Faculty of Infectious and Tropical Diseases, London School of Hygiene and Tropical Medicine, London, UK
Claudio Fronterre
& Jorge Cano
Federal Ministry of Health, Abuja, Nigeria
Ifeoma Anagbogu
& Chukwu Okoronkwo
Department of Statistics, University of Oxford, Oxford, UK
Christl A. Donnelly
Search for Obiora A. Eneanya in:
Search for Claudio Fronterre in:
Search for Ifeoma Anagbogu in:
Search for Chukwu Okoronkwo in:
Search for Tini Garske in:
Search for Jorge Cano in:
Search for Christl A. Donnelly in:
EOA, CAD, TG and JC conceived and designed the study. EOA carried out formal analysis and wrote the first draft of manuscript. JC, CAD and CF advised on spatial and machine learning models. IA and CO contributed survey data. All authors read and approved the final manuscript.
Correspondence to Obiora A. Eneanya.
The process of obtaining ethical approvals, informed consent, and arranging logistical procedures for field surveys were handled in-country by the Nigeria Ministry of Health, with technical support by the WHO.
Additional file 1: Figure S1. Cross-validation of the predicted lymphatic filariasis prevalence using a subsample of 25% of the observed data. Predicted prevalence values are plotted against observed prevalence. Those plotted in red had observations outside of the 95% prediction intervals while those plotted in blue had observations inside the 95% prediction intervals.
Additional file 2: Figure S2. a Mean predicted prevalence of lymphatic filariasis in Nigeria. b Predicted environmental suitability of lymphatic filariasis in Nigeria [28].
Eneanya, O.A., Fronterre, C., Anagbogu, I. et al. Mapping the baseline prevalence of lymphatic filariasis across Nigeria. Parasites Vectors 12, 440 (2019) doi:10.1186/s13071-019-3682-6
Antigenaemia
Microfilaraemia
The LCNTDR Collection: Advances in scientific research for NTD control