Predicting the effects of climate change on Schistosoma mansoni transmission in eastern Africa
© McCreesh et al.; licensee BioMed Central. 2014
Received: 24 October 2014
Accepted: 19 December 2014
Published: 6 January 2015
Survival and fitness attributes of free-living and sporocyst schistosome life-stages and their intermediate host snails are sensitive to water temperature. Climate change may alter the geographical distribution of schistosomiasis by affecting the suitability of freshwater bodies for hosting parasite and snail populations.
We have developed an agent-based model of the temperature-sensitive stages of the Schistosoma mansoni and intermediate host snail lifecycles. The model was run using low, moderate and high warming climate projections over eastern Africa. For each climate projection, eight model scenarios were used to determine the sensitivity of predictions to different relationships between air and water temperature, and different snail mortality rates. Maps were produced showing predicted changes in risk as a result of increasing temperatures over the next 20 and 50 years.
Baseline model output compared to prevalence data indicates suitable temperatures are necessary but not sufficient for both S. mansoni transmission and high infection prevalences. All else being equal, infection risk may increase by up to 20% over most of eastern Africa over the next 20 and 50 years. Increases may be higher in Rwanda, Burundi, south-west Kenya and eastern Zambia, and S. mansoni may become newly endemic in some areas. Results for 20-year projections are robust to changes in simulated intermediate host snail habitat conditions. There is greater uncertainty about the effects of different habitats on changes in risk in 50 years’ time.
Temperatures are likely to become suitable for increased S. mansoni transmission over much of eastern Africa. This may reduce the impact of control and elimination programmes. S. mansoni may also spread to new areas outside existing control programmes. We call for increased surveillance in areas defined as potentially suitable for emergent transmission.
The transmission potential of most neglected tropical diseases in a particular location is partly dependent on abiotic factors affecting either free-living life stages and/or those which occur in poikilothermic organisms such as snails and mosquitoes. As an exemplar, both the schistosome parasite and its intermediate host snails are very sensitive to water temperature [1-4]. Increasing temperatures in freshwater bodies in sub-tropical and tropical areas may therefore alter the geographic distribution of schistosomiasis. There is some empirical evidence that this may be occurring already in Uganda, with transmission occurring at altitudes previously considered too cold [5,6].
Climate change projections show increasing temperatures across Africa , where the majority of people infected with schistosome parasites are located . What remains unclear is how this phenomenon might affect the transmission potential of schistosomiasis in different locations, given the non-linear relationship between water temperature and schistosome transmission .
Perhaps because of the lack of research into this issue, the implications of climate change for schistosomiasis control and elimination have been largely ignored, and were not mentioned in the World Health Organization’s 2012 ‘Roadmap to implementation’ . Only a few mathematical models have explored the effects of temperature on schistosomiasis transmission [9-15], and while some incorporated projected increases in mean temperature, no previous models have been run using climate projections. Although many geostatistical models of schistosomiasis prevalence incorporate temperature (e.g. [16-18]) (and, in one case, climate change ), limitations of the empirical data and the complex, non-linear relationship between water temperature and infection risk  mean that they are unlikely to capture the full contribution of temperature. In this study, we advance the field by producing, for the first time, high-resolution maps of eastern Africa highlighting areas where temperatures may become suitable for increased or decreased transmission, and where Schistosoma mansoni may spread to new areas.
The main output of the model is ‘infection risk’, a measure of the mean number of cercariae in the model each hour, adjusted for the fact that their probability of causing infection decreases as the amount of time since they were shed from the snail increases. It can be viewed as an estimate of the suitability of temperatures for schistosome transmission in an area. The non-temperature dependent stages of the schistosome lifecycle – the stages in humans – are not simulated, and the rate of miracidium introduction into the model is constant.
The model has a time step of one hour, allowing the cercariae and miracidia, which have lifespans of the order of hours, to be accurately represented. Hourly temperatures are simulated using a sine wave, with maximum temperatures reached at 3 to 4 pm, and minimum temperatures at 3 to 4 am .
To allow schistosomes to become established in the model in new locations, a snail egg is introduced into the model each hour with a probability of 0.00012 (which gives an average rate of one snail egg per model year). As non-temperature dependent egg mortality is simulated using reduced egg production rates in the model, this is equivalent to a ‘real life’ egg introduction rate of 10 eggs per year.
An ensemble of regional climate simulations over Africa was used to provide projected daily maximum and minimum temperature data for the eastern Africa study region. The ensemble consists of three members of the Rossby Centre Regional Climate Model (RCM) – RCA4 , driven by a coupled atmosphere ocean general circulation model (AOGCM) – EC-EARTH . Three Representative Concentration Pathway (RCP) scenarios were used – RCP2.6, RCP4.5 and RCP8.5 – which represent low, moderate and high levels of warming respectively. All three regional simulations were made within the African branch of the Coordinated Regional Downscaling Experiment (CORDEX), and cover the whole African continent at about a 50 km (0.44°) resolution. A smaller sub-domain in eastern Africa was selected for the study and has size of 1470 (35 x 42) grid boxes.
The mean absolute increases in temperature over the study area over the next 20 years (2026–2035 relative to 2006–2015) were 0.34°C, 0.42°C, and 0.66°C in the low, moderate, and high warming projections, respectively. Over the next 50 years (2056–2065 relative to 2006–2035) the mean increases were 0.53°C, 1.2°C, and 1.8°C (Additional file 1: Figure S1). To prevent confusion, the different climate projections or scenarios are referred to as ‘projections’ throughout the paper, and the term ‘scenarios’ reserved for the model scenarios described below.
daily minimum and maximum water temperatures equal to daily minimum and maximum air temperatures;
minimum and maximum water temperatures 2°C higher;
minimum and maximum water temperatures 2°C lower; and
minimum water temperatures 2°C higher, and maximum water temperatures 2°C lower.
These were each simulated using two sets of snail mortality rates: rates estimated from experimental and field data , and double these rates. Therefore, in total, eight model scenarios were simulated for each of three climate projections.
The model was initially run using the first two years of the climate data (2006–2008) to allow snail numbers to reach plausible initial values. Following this, the model was run using the full climate data between 2006 and 2065. The model was run separately for each of the 1470 grid squares. For each location, the model was run 20 times and the results averaged. Outputs were averaged over each of three 10-year time periods: baseline (2006–2015), 20 years in the future (2026–2035), and 50 years in the future (2056–2065). Kernel interpolation using a fifth order polynomial function was used to produce smooth risk maps for each scenario and time period.
Comparison with empirical data
Geo-referenced data on the prevalence of S. mansoni in human populations in the model output region were extracted from the open-access Global Neglected Tropical Disease (GNTD) database . The database brings together historical and contemporary survey data on schistosomiasis and other neglected tropical diseases, and the data have been widely used in statistical models of disease prevalence (e.g. [18,19]). For grid squares containing more than one data point, the un-weighted mean of the prevalence estimates was calculated. The mean prevalence for each grid square was then plotted against the mean (across scenarios) model output ‘infection risk’ for the same square, and the area under the receiver operating characteristic curve (AUC) calculated for the ability of the model to predict prevalences of above 0%, 10%, 20% and 50%. The analysis was then repeated, restricting the prevalence data to estimates from children, collected from 2004 onwards, and where at least 10 children were tested (‘selected data’).
To enable easy interpretation of model results, for each climate projection and future time period median changes in risk (across model scenarios) were calculated to give a central estimate of the magnitude of changes that may occur. Means were not used as the predicted changes in each scenario were highly skewed in some locations, due to very high relative risks in some areas and scenarios where baseline risks were very low.
The level of agreement between scenarios in the direction of change in risk was also investigated. For each scenario, climate projection, and future time period, areas were given the value of +1 if risk was predicted to increase from its baseline value, −1 if risk was predicted to decrease, and zero if it was predicted to stay the same. Risk was considered to stay the same if it changed by less than ±10%. The values were then summed over all eight scenarios for each climate projection and future time period to give an indication of the overall predicted direction of change. For areas where both an increase and a decrease in risk were predicted by one or more scenarios, the number of scenarios that disagreed with the overall predicted direction of change was calculated.
Additional analyses were conducted on each of the model scenarios to explore the possibility of schistosome transmission becoming established at new sites. For these analyses, areas were given the value of +1 if risk was predicted to increase from below to above a cut-off between the two time periods, and zero otherwise. The analysis was conducted for all possible cut-offs between 1% and 99% of the maximum risk in the scenario in any time period, increasing in increments of 1% of the maximum risk. The cut-offs were split into three equal-sized groups (1-33%, 34-66%, and 67-99% of maximum risk), and the proportion of cut-offs crossed was calculated for each location, for each group, time period and climate projection.
Comparison with empirical data
Comparison of model output ‘infection risk’ at baseline (2006–2015) with empirical prevalence estimates
All data (N = 279)
Selected data (N = 100)
Number above cut-off (%)
AUC* (95% CI)
Number above cut-off (%)
AUC* (95% CI)
When all data were used, the AUC’s were 0.6 (0.52-0.67), 0.62 (0.55-0.69), 0.63 (0.56-0.71), and 0.57 (0.46-0.69) for prevalence cut-offs of above 0%, 10%, 20% and 50% respectively. When selected data were used, the AUC’s were 0.56 (0.45-0.68), 0.62 (0.47-0.77), 0.68 (0.50-0.86), and 0.76 (0.39-1.00) (Table 1). Figure 2b illustrates the mean model output at baseline divided by the maximum model output at any location, plotted against the prevalence data. The mean prevalence is higher than the adjusted model output in 9/279 (3.2%) grid squares when all data were used, and 1/100 (1.0%) grid squares when selected data were used.
Increase in risk
There is widespread agreement between scenarios that infection risk may decrease by more than 50% over the next 20 and 50 years in parts of north and east Kenya, southern South Sudan, and eastern People’s Democratic Republic of Congo. The size of the area over which reductions may occur is larger with higher levels of warming, and in 50 years’ time.
Newly endemic areas and new foci of transmission
The results suggest that, all else being equal, S. mansoni infection risk may increase across much of eastern Africa as temperatures increase over the next few decades. In most areas, the predicted increases are less than 20%. Temperature driven increases in risk may be much larger in Rwanda, Burundi, south-west Kenya, and eastern Zambia. Conversely, infection risk may decrease by more than 50% in parts of north and east Kenya, southern South Sudan, and eastern People’s Democratic Republic of Congo. The results also highlight areas where schistosome transmission may occur at new sites.
The results predict changes in risk that are attributable to increasing temperatures only. Other climatic changes, such as changes in patterns of rainfall, flooding and droughts, will also have an impact on future schistosomiasis prevalence . Changes in patterns of rainfall may be particularly important in determining the seasonality of transmission in many areas, and therefore results are not presented by season. In addition to climatic changes, non-climatic changes will play a large role in determining future prevalence and infection intensity. These include changes in land use, population growth and mobility, water contact behaviour and sanitation infrastructure. Potential modifying factors of risk also include mass treatment control programmes, and snail control efforts. For these reasons, the model results should not be taken as predictive of future schistosomiasis prevalence as such, but instead as indicative of areas where temperature changes are likely to influence schistosomiasis transmission potential.
A novel feature of the model presented here is consideration of the relationship between air and water temperature in different types of water body and in different seasons. Compared to minimum and maximum air temperatures, minimum and maximum water temperatures may be both higher [21,26], both lower , or may be higher and lower respectively [26,27]. Furthermore, snails may be capable of escaping above optimum water temperatures by moving to deeper water , where water temperatures may be up to 2°C cooler , or by burrowing in mud . For these reasons, a number of scenarios were simulated, with different relationships between air temperature and water temperature. Scenarios with high and low snail mortality rates explored the effect of differences between different types of habitat further. With the exception of the moderate and high warming projections in 50 years’ time, there is little disagreement between scenarios in the direction of change in risk, suggesting that the model results are robust to variations between snail habitats. The eight scenarios also act as a form of sensitivity analysis, as shifting temperatures and keeping snail and parasite temperature preferences constant is equivalent to keep the temperatures constant but altering model parameterisation.
Validation of the model is challenging as the model produces an estimate of the contribution of temperature to infection risk, whereas infection prevalence is the standard measure of population infection risk. Estimating prevalence was outside of the scope of this study due to a lack of available information (at high spatial resolution) on most of the underlying processes and phenomena that drive infection prevalence, including the presence or absence of suitable snail habitats, water contact behaviour, and sanitation infrastructure and usage. In addition, prevalence is a poor indicator of infection risk, and data on infection intensity are rarely routinely collected . Existing prevalence data are also of varying quality, having been collected for a wide range of different purposes, at different times, in different age groups, using different sampling methods, and using different tests with widely varying sensitivities . Areas where temperatures are considered to be unsuitable for schistosomiasis may be greatly under-sampled, and these are the areas where the predictive ability of the model is likely to be highest. These limitations both impede model validation, and reduce the ability of geostatistical models to accurately capture the effects of temperature.
Despite the limitations of the empirical data, a comparison of model output at baseline with prevalence data gives AUC’s of around 0.6 (Table 1), which suggests that temperature does play an important role in determining empirical prevalence (and that the model is able to capture that role), but that other factors are also important. Figure 2 illustrates clearly that zero or low prevalences of schistosomiasis can occur in any area, regardless of temperatures, but that higher prevalences only occur in areas where temperatures are suitable for higher levels of transmission. In other words, suitable temperatures are necessary but not sufficient for both schistosome transmission and for high prevalences of schistosomiasis.
Snail eggs are introduced into the model at an average rate of one a year. This is necessary to allow simulated snail populations to become established in areas where conditions become newly suitable. In real snail populations, this may occur through movement of snails from interconnected water bodies or during flooding , or through live snails being transported short distances on objects such as fishing nets. These events are all extremely stochastic, and the ‘newly endemic’ model results should therefore be taken to indicate areas that may become newly suitable for the establishment of schistosome transmission only, and not as a definite guide to the spread of schistosomiasis. The exception to this is cooler areas where snails are found but schistosome transmission does not currently occur, where new transmission foci could quickly become established as temperatures become suitable.
Taking regional climate simulations made by only one regional climate model, downscaling one global circulation model (under different RCP scenarios), is not enough to assess the full range of uncertainties in future climate projections. However, the main focus of the study is to provide a range of plausible future changes in schistosomiasis risk, suitable for public health planning purposes, and to explore the effect of uncertainties associated with the disease model and its parameterisation. Using a wider range of different climate projections (from different RCMs downscaling different AOGCMs) is unlikely to have any real effect on the public health consequences of our findings. As the main focus of the study was infectious disease modelling, there was also no bias correction applied to the RCM data.
To maximise the potential utility of our results for public health purposes, we used short (10-year) time slices for the analysis, namely: 2006–2015 (baseline), 2026–2035 and 2056–2065. Using 10-year time periods can lead to errors if natural multi-decadal variability is attributed to the global warming signal. We do not believe that this was the case for our study however, as climate projections for daily maximum and minimum temperatures in eastern Africa averaged over 10-years were robust and consistent across the three RCP trajectories (Additional file 1: Figure S1).
The model was parameterised using data from B. pfeifferi snails. This is the most widespread intermediate host snail species in sub-Saharan Africa . Other species of intermediate host snail are found at many transmission sites however, and previous modelling work has shown that the species of snail used in parameterising a model can have a large effect on the relationship between water temperature and infection risk . Simulating different relationships between air temperature and water temperature goes some way towards testing the sensitivity of the model to the choice of snail species, as the higher and lower simulated temperatures could also be viewed as alterations in minimum, maximum, and optimum temperatures for snail survival and reproduction. The widespread agreement between scenarios therefore suggests that the model is not very sensitive to small differences in temperature preferences between snail species. Nevertheless, experimental work with the different species of intermediate host snail found in eastern Africa would increase confidence in the model results by allowing the model to be explicitly parameterised to other species of snail.
A practical application of the model presented here is associated with planning the long-term control of schistosomiasis. Current efforts rely on periodic mass administration of medicines within specific areas. As climate change begins to exert an influence on the African environment, the situation on the ground is likely to change. Mitigation and/or adaptation measures to climate change will require greater surveillance efforts to capture those changes in advance of them becoming significantly deleterious. A rapid mapping of snail populations within the areas identified as being at risk of new endemicity would pinpoint the communities at highest risk, allowing targeted monitoring of human populations.
Climate change may have a large effect on both the distribution and intensity of S. mansoni infection over coming decades. Temperatures are predicted to become suitable for increased transmission over much of eastern Africa over the next 20 years. This may lead to increased prevalences and intensities of infection in some areas, and is likely to reduce the effects of control and elimination programmes. In some areas, particularly in Rwanda, Burundi, south-west Kenya and eastern Zambia, increases in infection risk may be large. Schistosomiasis may spread to new areas, outside the current range of control programmes. Increased surveillance of these areas would enable education and control programmes to be promptly implemented in newly endemic areas, minimising disease morbidity.
The research leading to these results received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 266327. The funders had no role in the design of the experiments; the collection, analysis, and interpretation of data; in the writing of the manuscript; or in the decision to submit the manuscript for publication.
- Foster R: The effect of temperature on the development of Schistosoma mansoni Sambon 1907 in the intermediate host. J Trop Med Hyg 1964, 67:289–292.PubMedGoogle Scholar
- Appleton CC: Review of literature on abiotic factors influencing the distribution and life cycle of bilharzia intermediate host snails. Malacol Rev 1978, 11:1–25.Google Scholar
- Pfluger W: Experimental epidemiology of schistosomiasis. I. The prepatent period and cercarial production of Schistosoma mansoni in Biomphalaria snails at various constant temperatures. Z Parasitenkd 1980, 63:159–169.PubMedView ArticleGoogle Scholar
- McCreesh N, Arinaitwe M, Arineitwe W, Tukahebwa EM, Booth M: Effect of water temperature and population density on the population dynamics of Schistosoma mansoni intermediate host snails. Parasit Vectors 2014, 7:503.PubMed CentralPubMedView ArticleGoogle Scholar
- Kabatereine NB, Brooker S, Tukahebwa EM, Kazibwe F, Onapa AW: Epidemiology and geography of Schistosoma mansoni in Uganda: implications for planning control. Trop Med Int Health 2004, 9:372–380.PubMedView ArticleGoogle Scholar
- Rubaihayo J, Moghusu E, Clouds P, Abaasa A: Schistosomiasis transmission at high altitude crater lakes in Western Uganda. BMC Infect Dis 2008, 8:110.View ArticleGoogle Scholar
- Stocker TF, Qin D, Plattner G-K, Tignor M, Allen SK, Boschung J, Nauels A, Xia Y, Bex V, Midgley PM: Climate change 2013: The physical science basis. In Intergovernmental Panel on Climate Change, Working Group I Contribution to the IPCC Fifth Assessment Report (AR5). New York: Cambridge Univ Press; 2013.Google Scholar
- World Health Organization: Accelerating work to overcome the global impact of neglected tropical diseases: a roadmap for implementation. Geneva, Switzerland: World Health Organization; 2012.Google Scholar
- McCreesh N, Booth M: The effect of increasing water temperatures on Schistosoma mansoni transmission and Biomphalaria pfeifferi population dynamics: An agent-based modelling study. PLoS One 2014, 9:e101462.PubMed CentralPubMedView ArticleGoogle Scholar
- Martens WJM, Jetten TH, Rotmans J, Niessen LW: Climate change and vector-borne diseases: a global modelling perspective. Glob Environ Chang 1995, 5:195–209.View ArticleGoogle Scholar
- Martens WJM, Jetten TH, Focks DA: Sensitivity of malaria, schistosomiasis and dengue to global warming. Clim Change 1997, 35:145–156.View ArticleGoogle Scholar
- Mangal TD, Paterson S, Fenton A: Predicting the impact of long-term temperature changes on the epidemiology and control of schistosomiasis: a mechanistic model. PLoS One 2008, 3:e1438.PubMed CentralPubMedView ArticleGoogle Scholar
- Zhou XN, Yang GJ, Yang K, Wang XH, Hong QB, Sun LP, Malone JB, Kristensen TK, Bergquist NR, Utzinger J: Potential impact of climate change on schistosomiasis transmission in China. Am J Trop Med Hyg 2008, 78:188–194.PubMedGoogle Scholar
- Mas-Coma S, Valero MA, Bargues MD: Climate change effects on trematodiases, with emphasis on zoonotic fascioliasis and schistosomiasis. Vet Parasitol 2009, 163:264–280.PubMedView ArticleGoogle Scholar
- McCreesh N, Booth M: The effect of simulating different intermediate host snail species on the link between water temperature and schistosomiasis risk. PLoS One 2014, 9:e87892.PubMed CentralPubMedView ArticleGoogle Scholar
- Stensgaard A, Jorgensen A, Kabatereine N, Malone J, Kristensen T: Modeling the distribution of Schistosoma mansoni and host snails in Uganda using satellite sensor data and Geographical information systems. Parassitologia 2005, 47:115–125.PubMedGoogle Scholar
- Brooker S: Spatial epidemiology of human schistosomiasis in Africa: risk models, transmission dynamics and control. Trans R Soc Trop Med Hyg 2007, 101:1–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Schur N, Hürlimann E, Stensgaard A-S, Chimfwembe K, Mushinge G, Simoonga C, Kabatereine NB, Kristensen TK, Utzinger J, Vounatsou P: Spatially explicit Schistosoma infection risk in eastern Africa using Bayesian geostatistical modelling. Acta Trop 2013, 128:365–377.PubMedView ArticleGoogle Scholar
- Stensgaard A-S, Utzinger J, Vounatsou P, Hürlimann E, Schur N, Saarnak CFL, Simoonga C, Mubita P, Kabatereine NB, Tchuenté L-AT, Rahbek C, Kristensen TK: Large-scale determinants of intestinal schistosomiasis and intermediate host snail distribution across Africa: Does climate matter? Acta Trop 2013, 128:378–390.PubMedView ArticleGoogle Scholar
- Wilensky U: Netlogo. In Center for Connected Learning and Computer-Based Modeling. Evanston, IL: Northwestern University; 1999.Google Scholar
- Paaijmans KP, Jacobs AFG, Takken W, Heusinkveld BG, Githeko AK, Dicke M, Holtslag AAM: Observations and model estimates of diurnal water temperature dynamics in mosquito breeding sites in western Kenya. Hydrol Process 2008, 22:4789–4801.View ArticleGoogle Scholar
- Nikulin G, Jones C, Giorgi F, Asrar G, Büchner M, Cerezo-Mota R, Christensen OB, Déqué M, Fernandez J, Hänsler A: Precipitation climatology in an ensemble of CORDEX-Africa regional climate simulations. J Clim 2012, 25:6057–6078.View ArticleGoogle Scholar
- Hazeleger W, Severijns C, Semmler T, Stefanescu S, Yang S, Wang X, Wyser K, Dutra E, Baldasano JM, Bintanja R: EC-Earth: a seamless earth-system prediction approach in action. Bull Am Meteorol Soc 2010, 91:1357–1363.View ArticleGoogle Scholar
- Hürlimann E, Schur N, Boutsika K, Stensgaard A-S, de Himpsl ML, Ziegelbauer K, Laizer N, Camenzind L, Di Pasquale A, Ekpo UF, Simoonga C, Mushinge G, Saarnak CFL, Utzinger J, Kristensen TK, Vounatsou P: Toward an open-access global database for mapping, control, and surveillance of neglected tropical diseases. PLoS Negl Trop Dis 2011, 5:e1404.PubMed CentralPubMedView ArticleGoogle Scholar
- McCreesh N, Booth M: Challenges in predicting the effects of climate change on Schistosoma mansoni and Schistosoma haematobium transmission potential. Trends Parasitol 2013, 29:548–555.PubMedView ArticleGoogle Scholar
- Jacobs AG, Heusinkveld B, Kraai A, Paaijmans K: Diurnal temperature fluctuations in an artificial small shallow water body. Int J Biometeorol 2008, 52:271–280.PubMed CentralPubMedView ArticleGoogle Scholar
- Adebisi A: The physico-chemical hydrology of a tropical seasonal river-upper Ogun river. Hydrobiologia 1981, 79:157–165.View ArticleGoogle Scholar
- Beadle L: Prolonged stratification and deoxygenation in tropical lakes. I. Crater Lake Nkugute, Uganda, compared with Lakes Bunyoni and Edward. Limnol Oceanogr 1966, 11:152–163.View ArticleGoogle Scholar
- Anderson R, Truscott J, Hollingsworth TD: The coverage and frequency of mass drug administration required to eliminate persistent transmission of soil-transmitted helminths. Phil Trans R Soc B 2014, 369:20130435.PubMed CentralPubMedView ArticleGoogle Scholar
- Booth M, Vounatsou P, N’goran E, Tanner M, Utzinger J: The influence of sampling effort and the performance of the Kato-Katz technique in diagnosing Schistosoma mansoni and hookworm co-infections in rural Côte d’Ivoire. Parasitology 2003, 127:525–531.PubMedView ArticleGoogle Scholar
- Wu XH, Zhang SQ, Xu XJ, Huang YX, Steinmann P, Utzinger J, Wang TP, Xu J, Zheng J, Zhou XN: Effect of floods on the transmission of schistosomiasis in the Yangtze River valley, People’s Republic of China. Parasitol Int 2008, 57:271–276.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.