Geographical distribution of human Schistosoma japonicum infection in The Philippines: tools to support disease control and further elimination

Ricardo J. Soares Magalhães, Maria Sonia Salamat, Lydia Leonardo, Darren J. Gray, Hélène Carabin, Kate Halton, Donald P. McManus, Gail M. Williams, Pilarita Rivera, Ofelia Saniel, Leda Hernandez, Laith Yakob, Stephen McGarvey, Archie Clements

Schistosoma japonicum infection is believed to be endemic in 28 of the 80 provinces of The Philippines and the most recent data on schistosomiasis prevalence have shown considerable variability between provinces. In order to increase the efficient allocation of parasitic disease control resources in the country, we aimed to describe the small-scale spatial variation in S. japonicum prevalence across The Philippines, quantify the role of the physical environment in driving the spatial variation of S. japonicum, and develop a predictive risk map of S. japonicum infection. Data on S. japonicum infection from 35,754 individuals across the country were geo-located at the barangay level and included in the analysis. The analysis was then stratified geographically for the regions of Luzon, the Visayas and Mindanao. Zero-inflated binomial Bayesian geostatistical models of S. japonicum prevalence were developed and diagnostic uncertainty was incorporated. Results of the analysis show that in the three regions, males and individuals aged ⩾ 20 years had significantly higher prevalence of S. japonicum compared with females and children < 5 years. The role of the environmental variables differed between regions of The Philippines. Schistosoma japonicum infection was widespread in the Visayas whereas it was much more focal in Luzon and Mindanao. This analysis revealed significant spatial variation in the prevalence of S. japonicum infection in The Philippines. This suggests that a spatially targeted approach to schistosomiasis interventions, including mass drug administration, is warranted. When financially possible, additional schistosomiasis surveys should be prioritized for areas identified to be at high risk but which were under-represented in our dataset.