The geographical range of schistosomiasis is affected by the ecology of schistosome parasites and their obligate host snails, including their response to temperature. Previous models predicted schistosomiasis' thermal optimum at 21.7 °C, which is not compatible with the temperature in sub-Saharan Africa (SSA) regions where schistosomiasis is hyperendemic. We performed an extensive literature search for empirical data on the effect of temperature on physiological and epidemiological parameters regulating the free-living stages of S. mansoni and S. haematobium and their obligate host snails, i.e., Biomphalaria spp. and Bulinus spp., respectively. We derived nonlinear thermal responses fitted on these data to parameterize a mechanistic, process-based model of schistosomiasis. We then re-cast the basic reproduction number and the prevalence of schistosome infection as functions of temperature. We found that the thermal optima for transmission of S. mansoni and S. haematobium range between 23.1-27.3 °C and 23.6-27.9 °C (95 % CI) respectively. We also found that the thermal optimum shifts toward higher temperatures as the human water contact rate increases with temperature. Our findings align with an extensive dataset of schistosomiasis prevalence in SSA. The refined nonlinear thermal-response model developed here suggests a more suitable current climate and a greater risk of increased transmission with future warming for more than half of the schistosomiasis suitable regions with mean annual temperature below the thermal optimum.
Keywords: Basic reproduction number; Climate Change; Mechanistic models; Schistosomiasis; Thermal response.