Geostatistical analysis of active human cysticercosis: Results of a large-scale study in 60 villages in Burkina Faso

Cysticercosis is a neglected tropical disease caused by the larval stage of the zoonotic tapeworm (Taenia solium). While there is a clear spatial component in the occurrence of the par-asite, no geostatistical analysis of active human cysticercosis has been conducted yet, nor has such an analysis been conducted for Sub-Saharan Africa, albeit relevant for guiding prevention and control strategies. The goal of this study was to conduct a geostatistical analysis of active human cysticercosis, using data from the baseline cross-sectional component of a large-scale study in 60 villages in Burkina Faso. The outcome was the prevalence of active human cysticercosis (hCC), determined using the B158/B60 Ag-ELISA, while various environmental variables linked with the transmission and spread of the disease were explored as potential explanatory variables for the spatial distribution of T. solium. A gener-alized linear geostatistical model (GLGM) was run, and prediction maps were generated. Analyses were conducted using data generated at two levels: individual participant data and grouped village data. The best model was selected using a backward variable selection procedure and models were compared using likelihood ratio testing. The best individual-level GLGM included precipitation (increasing values were associated with an increased odds of positive test result), distance to the nearest river (decreased odds) and night land tempera-ture (decreased odds) as predictors for active hCC, whereas the village-level GLGM only retained precipitation and distance to the nearest river. The range of spatial correlation was estimated at 45.0 [95%CI: 34.3; 57.8] meters and 28.2 [95%CI: 14.0; 56.2] km for the indi-vidual-and village-level datasets, respectively. Individual-and village-level GLGM unrav-elled large areas with active hCC predicted prevalence estimates of at least 4% in the south-east, the extreme south, and north-west of the study area, while patches of prevalence estimates below 2% were seen in the north and west. More research designed to analyse the spatial characteristics of hCC is needed with sampling strategies ensuring appropriate char-acterisation of spatial variability, and incorporating the uncertainty linked to the measure-ment of outcome and environmental variables in the geostatistical analysis. Trial registration: ClinicalTrials.gov; NCT0309339.