Model-Based Geostatistical Mapping of the Prevalence of Onchocerca volvulus in West Africa

Background The initial endemicity (pre-control prevalence) of onchocerciasis has been shown to be an important determinant of the feasibility of elimination by mass ivermectin distribution. We present the first geostatistical map of microfilarial prevalence in the former Onchocerciasis Control Programme in West Africa (OCP) before commencement of antivectorial and antiparasitic interventions. Methods and Findings Pre-control microfilarial prevalence data from 737 villages across the 11 constituent countries in the OCP epidemiological database were used as ground-truth data. These 737 data points, plus a set of statistically selected environmental covariates, were used in a Bayesian model-based geostatistical (B-MBG) approach to generate a continuous surface (at pixel resolution of 5 km x 5km) of microfilarial prevalence in West Africa prior to the commencement of the OCP. Uncertainty in model predictions was measured using a suite of validation statistics, performed on bootstrap samples of held-out validation data. The mean Pearson’s correlation between observed and estimated prevalence at validation locations was 0.693; the mean prediction error (average difference between observed and estimated values) was 0.77%, and the mean absolute prediction error (average magnitude of difference between observed and estimated values) was 12.2%. Within OCP boundaries, 17.8 million people were deemed to have been at risk, 7.55 million to have been infected, and mean microfilarial prevalence to have been 45% (range: 2–90%) in 1975. Conclusions and Significance This is the first map of initial onchocerciasis prevalence in West Africa using B-MBG. Important environmental predictors of infection prevalence were identified and used in a model out-performing those without spatial random effects or environmental covariates. Results may be compared with recent epidemiological mapping efforts to find areas of persisting transmission. These methods may be extended to areas where data are sparse, and may be used to help inform the feasibility of elimination with current and novel tools.

[1]  M. Boussinesq,et al.  Effect of single-dose ivermectin on Onchocerca volvulus: a systematic review and meta-analysis. , 2008, The Lancet. Infectious diseases.

[2]  J. Remmé,et al.  Distribution and severity of onchocerciasis in southern Benin, Ghana and Togo. , 1992, Acta tropica.

[3]  M. Basáñez,et al.  Modelling the impact of ivermectin on River Blindness and its burden of morbidity and mortality in African Savannah: EpiOncho projections , 2014, Parasites & Vectors.

[4]  P. Ngoumou,et al.  A rapid mapping technique for the prevalence and distribution of onchocerciasis: a Cameroon case study. , 1994, Annals of tropical medicine and parasitology.

[5]  Moussa S. Sanfo,et al.  Validation of a Remote Sensing Model to Identify Simulium damnosum s.l. Breeding Sites in Sub-Saharan Africa , 2013, PLoS neglected tropical diseases.

[6]  Wilson,et al.  Deforestation and the spatio-temporal distribution of savannah and forest members of the Simulium damnosum complex in southern Ghana and south-western Togo. , 2002, Transactions of the Royal Society of Tropical Medicine and Hygiene.

[7]  Damien Sulla-Menashe,et al.  MODIS Collection 5 global land cover: Algorithm refinements and characterization of new datasets , 2010 .

[8]  M. Basáñez,et al.  Reaching the London Declaration on Neglected Tropical Diseases Goals for Onchocerciasis: An Economic Evaluation of Increasing the Frequency of Ivermectin Treatment in Africa , 2014, Clinical infectious diseases : an official publication of the Infectious Diseases Society of America.

[9]  P. Diggle,et al.  The geographic distribution of onchocerciasis in the 20 participating countries of the African Programme for Onchocerciasis Control: (2) pre-control endemicity levels and estimated number infected , 2014, Parasites & Vectors.

[10]  María-Gloria Basáñez,et al.  Prediction of community prevalence of human onchocerciasis in the Amazonian onchocerciasis focus: Bayesian approach. , 2003, Bulletin of the World Health Organization.

[11]  E. Pebesma,et al.  Classes and Methods for Spatial Data , 2015 .

[12]  H. Duerr,et al.  High infection rates at low transmission potentials in West African onchocerciasis. , 2006, International journal for parasitology.

[13]  M. Thomson,et al.  The transmission of onchocerciasis at a forest village in Sierra Leone. I. Simulium damnosum s.l. biting densities and infection with Onchocerca volvulus at five representative sites. , 1990, Annals of Tropical Medicine and Parasitology.

[14]  R. Feachem.  Disease and Mortality in Sub-Saharan Africa , 1991 .

[15]  I. Tada,et al.  Onchocerciasis in San Vicente Pacaya, Guatemala. , 1979, The American journal of tropical medicine and hygiene.

[16]  R. Baker,et al.  Onchocerciasis distribution and severity in five West African countries. , 1991, Bulletin of the World Health Organization.

[17]  J. Maguire,et al.  Forest onchocerciasis in Cameroon: its distribution and implications for selection of communities for control programmes. , 1997, Annals of tropical medicine and parasitology.

[18]  J. Hervouët,et al.  Les niveaux d'endémicité dans l'onchocercose , 1979 .

[19]  Rob J Hyndman,et al.  Another look at measures of forecast accuracy , 2006 .

[20]  S. Wanji,et al.  Integrated rapid mapping of onchocerciasis and loiasis in the Democratic Republic of Congo: impact on control strategies. , 2011, Acta tropica.

[21]  Guo-Jing Yang,et al.  A Research Agenda for Helminth Diseases of Humans: Modelling for Control and Elimination , 2012, PLoS neglected tropical diseases.

[22]  D. Boakye,et al.  Prevalence and intensity of Onchocerca volvulus infection and efficacy of ivermectin in endemic communities in Ghana: a two-phase epidemiological study , 2007, The Lancet.

[23]  James P. Verdin,et al.  A topological system for delineation and codification of the Earth's river basins , 1999 .

[24]  M. Little,et al.  Density-Dependent Mortality of the Human Host in Onchocerciasis: Relationships between Microfilarial Load and Excess Mortality , 2012, PLoS neglected tropical diseases.

[25]  P. Enyong,et al.  Rapid epidemiological mapping of onchocerciasis (REMO): its application by the African Programme for Onchocerciasis Control (APOC) , 2002, Annals of tropical medicine and parasitology.

[26]  Edwin W. Pak,et al.  An extended AVHRR 8‐km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data , 2005 .

[27]  S. Cousens,et al.  Impact of long-term treatment of onchocerciasis with ivermectin in Kaduna State, Nigeria: first evidence of the potential for elimination in the operational area of the African Programme for Onchocerciasis Control , 2012, Parasites & Vectors.

[28]  Wolfgang Schramm,et al.  Team , 2018, Spaces of Intensity.

[29]  David L. Smith,et al.  A World Malaria Map: Plasmodium falciparum Endemicity in 2007 , 2009, PLoS medicine.

[30]  B. Thylefors,et al.  Ocular onchocerciasis and intensity of infection in the community. II. West African rainforest foci of the vector Simulium yahense. , 1989, Tropical medicine and parasitology : official organ of Deutsche Tropenmedizinische Gesellschaft and of Deutsche Gesellschaft fur Technische Zusammenarbeit.

[31]  B. Kirkwood,et al.  Variations in the prevalence and intensity of microfilarial infections by age, sex, place and time in the area of the Onchocerciasis Control Programme. , 1983, Transactions of the Royal Society of Tropical Medicine and Hygiene.

[32]  A. Prost,et al.  Essai de normalisation de la méthodologie des enquêtes clinico-parasitologiques sur l'onchocercose en Afrique de l'Ouest , 1978 .

[33]  A. Prost,et al.  [Epidemiologic status of onchocerciasis]. , 1979, Bulletin of the World Health Organization.

[34]  P. Wenk,et al.  Studies on the dynamics of transmission of onchocerciasis in a Sudan-savanna area of North Cameroon I. Prevailing Simulium vectors, their biting rates and age-composition at different distances from their breeding sites. , 1987, Annals of tropical medicine and parasitology.

[35]  F. White.  The vegetation of Africa : a descriptive memoir to accompany the Unesco/AETFAT/UNSO vegetation map of Africa , 1985 .

[36]  M. Thomson,et al.  Satellite mapping of Loa loa prevalence in relation to ivermectin use in west and central Africa , 2000, The Lancet.

[37]  M. Boussinesq,et al.  Country-wide rapid epidemiological mapping of onchocerciasis (REMO) in Cameroon. , 1997, Annals of tropical medicine and parasitology.

[38]  J. Remmé,et al.  A force-of-infection model for onchocerciasis and its applications in the epidemiological evaluation of the Onchocerciasis Control Programme in the Volta River basin area. , 1986, Bulletin of the World Health Organization.

[39]  A. Gorbalenya,et al.  Topley and Wilson's Microbiology and Microbial Infections , 2005 .

[40]  D. Boakye,et al.  Onchocerciasis Control: Vision for the Future from a Ghanian perspective , 2009, Parasites & Vectors.

[41]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[42]  A. Ghani,et al.  Estimating Air Temperature and Its Influence on Malaria Transmission across Africa , 2013, PloS one.

[43]  J. B. Davies,et al.  Onchocerciasis transmission levels and Simulium damnosum complex biting activity at riverside and rice field sites in Sierra Leone , 1988, Medical and veterinary entomology.

[44]  M. Boussinesq,et al.  Severe adverse reaction risks during mass treatment with ivermectin in loiasis-endemic areas. , 1996, Parasitology today.

[45]  Penelope Vounatsou,et al.  Bayesian Geostatistical Modeling of Leishmaniasis Incidence in Brazil , 2013, PLoS neglected tropical diseases.

[46]  M. Tanner,et al.  Soil-transmitted helminth infection in South America: a systematic review and geostatistical meta-analysis. , 2013, The Lancet. Infectious diseases.

[47]  M. Murdoch,et al.  A clinical classification and grading system of the cutaneous changes in onchocerciasis , 1993, The British journal of dermatology.

[48]  D. Boakye,et al.  Onchocerciasis Transmission in Ghana: Persistence under Different Control Strategies and the Role of the Simuliid Vectors , 2015, PLoS neglected tropical diseases.

[49]  R. Tibshirani,et al.  Improvements on Cross-Validation: The 632+ Bootstrap Method , 1997 .

[50]  Edzer J. Pebesma,et al.  Applied Spatial Data Analysis with R - Second Edition , 2008, Use R!.

[51]  B. Boatin.  The Current State of the Onchocerciasis Control Programme in West Africa , 2003, Tropical doctor.

[52]  María-Gloria Basáñez,et al.  River Blindness: A Success Story under Threat? , 2006, PLoS medicine.

[53]  L. Konaté,et al.  Feasibility of Onchocerciasis Elimination with Ivermectin Treatment in Endemic Foci in Africa: First Evidence from Studies in Mali and Senegal , 2009, PLoS neglected tropical diseases.

[54]  J. Remmé,et al.  Detailed epidemiological mapping of three onchocerciasis foci in West Africa. , 1991, Acta tropica.

[55]  P J Diggle,et al.  Spatial modelling and the prediction of Loa loa risk: decision making under uncertainty , 2007, Annals of tropical medicine and parasitology.

[56]  B. Duke,et al.  Studies on factors influencing the transmission of onchocerciasis. VII. A comparison of the Onchocerca volvulus transmission potentials of Simulium damnosum populations in four Cameroon rain-forest villages and the pattern of onchocerciasis associated therewith. , 1972, Annals of tropical medicine and parasitology.

[57]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[58]  R. Baker,et al.  Ocular onchocerciasis and intensity of infection in the community. III. West African rainforest foci of the vector Simulium sanctipauli. , 1990, Tropical medicine and parasitology : official organ of Deutsche Tropenmedizinische Gesellschaft and of Deutsche Gesellschaft fur Technische Zusammenarbeit.

[59]  Young Eun Kim,et al.  Control, Elimination, and Eradication of River Blindness: Scenarios, Timelines, and Ivermectin Treatment Needs in Africa , 2015, PLoS neglected tropical diseases.

[60]  M. Basáñez,et al.  Association between microfilarial load and excess mortality in onchocerciasis: an epidemiological study , 2004, The Lancet.

[61]  C. Vajime,et al.  The distribution of the simulium damnosum complex in West Africa with particular reference to the onchocerciasis control programme area. , 1978, Tropenmedizin und Parasitologie.

[62]  Peter C Austin,et al.  Bootstrap Methods for Developing Predictive Models , 2004 .

[63]  L. Coffeng,et al.  Onchocerciasis: The Pre-control Association between Prevalence of Palpable Nodules and Skin Microfilariae , 2013, PLoS neglected tropical diseases.

[64]  Dean B. Gesch,et al.  TECHNIQUES FOR DEVELOPMENT OF GLOBAL I-KILOMETER DIGITAL ELEVATION MODELS , 2003 .

[65]  R. Frentzel-Beyme.  The geographical distribution of Onchocerca volvulus infection in Liberia. , 1975, Tropenmedizin und Parasitologie.

[66]  D. Boakye,et al.  Potential effects of warmer worms and vectors on onchocerciasis transmission in West Africa , 2015, Philosophical Transactions of the Royal Society B: Biological Sciences.

[67]  E. Samba,et al.  The Onchocerciasis Control Programme in West Africa: An Example of Effective Public Health Management , 1994 .

[68]  B. Boatin.  The Onchocerciasis Control Programme in West Africa (OCP) , 2008, Annals of tropical medicine and parasitology.

[69]  J. Habbema,et al.  Can ivermectin mass treatments eliminate onchocerciasis in Africa? , 2002, Bulletin of the World Health Organization.

[70]  L. Coffeng,et al.  African Programme for Onchocerciasis Control 1995–2015: Model-Estimated Health Impact and Cost , 2013, PLoS neglected tropical diseases.

[71]  Tim Burt,et al.  Topographic controls of soil moisture distributions , 1985 .

[72]  Anand P. Patil,et al.  Bayesian geostatistics in health cartography: the perspective of malaria. , 2011, Trends in parasitology.

[73]  M. Grillet,et al.  Landscape Epidemiology of Human Onchocerciasis in Southern Venezuela , 2011 .

[74]  D. Boakye,et al.  Onchocerciasis transmission in Ghana: biting and parous rates of host-seeking sibling species of the Simulium damnosum complex , 2014, Parasites & Vectors.

[75]  Peter M. Atkinson,et al.  Bayesian Geostatistical Analysis and Prediction of Rhodesian Human African Trypanosomiasis , 2010, PLoS neglected tropical diseases.

[76]  A. N. Strahler.  DYNAMIC BASIS OF GEOMORPHOLOGY , 1952 .

[77]  Molly E. Brown,et al.  EMD CORRECTION OF ORBITAL DRIFT ARTIFACTS IN SATELLITE DATA STREAM , 2010 .

[78]  B. Boatin,et al.  Twenty-two years of blackfly control in the onchocerciasis control programme in West Africa. , 1997, Parasitology today.

[79]  Danielle K. Golon.  The Land Processes Distributed Active Archive Center (LP DAAC) , 2016 .

[80]  J. Walsh,et al.  Studies on the reinvasion of the Onchocerciasis Control Programme in the Volta River Basin by Simulium damnosum s.I. with emphasis on the south-western areas. , 1979, Tropenmedizin und Parasitologie.

[81]  K. Beven,et al.  A physically based, variable contributing area model of basin hydrology , 1979 .

[82]  J. L. Parra,et al.  Very high resolution interpolated climate surfaces for global land areas , 2005 .

[83]  M. Grillet,et al.  Onchocerciasis in the Amazonian focus of southern Venezuela: altitude and blackfly species composition as predictors of endemicity to select communities for ivermectin control programmes. , 1998, Transactions of the Royal Society of Tropical Medicine and Hygiene.

[84]  E. LESTER SMITH,et al.  AND OTHERS , 2005 .

[85]  Peter M Atkinson,et al.  Developing geostatistical space-time models to predict outpatient treatment burdens from incomplete national data. , 2008, Geographical analysis.

[86]  Archie C A Clements,et al.  The applications of model-based geostatistics in helminth epidemiology and control. , 2011, Advances in parasitology.

[87]  B. Thylefors,et al.  Ocular onchocerciasis and intensity of infection in the community. I. West African savanna. , 1989, Tropical medicine and parasitology : official organ of Deutsche Tropenmedizinische Gesellschaft and of Deutsche Gesellschaft fur Technische Zusammenarbeit.

[88]  D. Skole,et al.  Land Use and Land Cover Change , 2014 .

[89]  P. Diggle,et al.  Model‐based geostatistics , 2007 .

[90]  C. Back,et al.  Sibling species distributions of the Simulium damnosum complex in the West African Onchocerciasis Control Programme area during the decade 1984–93, following intensive larviciding since 1974 , 1998, Medical and veterinary entomology.

[91]  M. Boussinesq,et al.  Serious reactions after mass treatment of onchocerciasis with ivermectin in an area endemic for Loa loa infection , 1997, The Lancet.

[92]  Simon Brooker,et al.  Spatial heterogeneity of parasite co-infection: Determinants and geostatistical prediction at regional scales , 2009, International journal for parasitology.

[93]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[94]  Burkina Faso,et al.  Eliminating onchocerciasis after 14 years of vector control: a proved strategy. , 2001, The Journal of infectious diseases.

[95]  David L. Smith,et al.  A new world malaria map: Plasmodium falciparum endemicity in 2010 , 2011, Malaria Journal.

[96]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[97]  M. Little,et al.  Transmission intensity and the patterns of Onchocerca volvulus infection in human communities. , 2002, The American journal of tropical medicine and hygiene.

[98]  T. C. Haas,et al.  Model-based geostatistics. Discussion. Authors' reply , 1998 .

[99]  S. Brooker.  Spatial epidemiology of human schistosomiasis in Africa: risk models, transmission dynamics and control , 2007, Transactions of the Royal Society of Tropical Medicine and Hygiene.

[100]  L. Konaté,et al.  Proof-of-Principle of Onchocerciasis Elimination with Ivermectin Treatment in Endemic Foci in Africa: Final Results of a Study in Mali and Senegal , 2012, PLoS neglected tropical diseases.

[101]  Brian D. Ripley,et al.  geoRglm: A Package for Generalised Linear Spatial Models , 2002 .

[102]  P. Diggle,et al.  The Geographic Distribution of Loa loa in Africa: Results of Large-Scale Implementation of the Rapid Assessment Procedure for Loiasis (RAPLOA) , 2011, PLoS neglected tropical diseases.

[103]  Hannah C. Slater,et al.  Mapping, Bayesian Geostatistical Analysis and Spatial Prediction of Lymphatic Filariasis Prevalence in Africa , 2013, PloS one.

[104]  M. Basáñez,et al.  Uncertainty Surrounding Projections of the Long-Term Impact of Ivermectin Treatment on Human Onchocerciasis , 2013, PLoS neglected tropical diseases.

[105]  M. Basáñez,et al.  The potential impact of moxidectin on onchocerciasis elimination in Africa: an economic evaluation based on the Phase II clinical trial data , 2015, Parasites & Vectors.