Background
Human schistosomiasis, caused by five major species of the trematode
Schistosoma, is endemic in 78 countries, and more than 800 million people worldwide live in areas endemic for schistosomiasis [
1
]. According to the road map for neglected tropical diseases (NTDs) 2021-2030 issued by the World Health Organization (WHO), all endemic countries are expected to achieve the elimination of schistosomiasis as a public health problem and 25 countries should achieve transmission interruption by 2030 [
2
]. Asian schistosomiasis is endemic in five countries, i.e., Cambodia, the People’s Republic of China (P.R. China), Indonesia, Laos, and the Philippines. According to the WHO’s NTDs roadmap 2021-2030, Asian schistosomiasis will be eliminated by 2030.
Schistosomiasis is a parasitic disease with a life cycle that requires a snail intermediate host [
3
]. For this reason, snail control is an important means of control. The 65th World Health Assembly (WHA) in 2012 proposed taking advantage of programmes outside the immediate health area, such as general improvement of the environment, to accelerate the interruption of schistosomiasis transmission by elimination of the snail intermediate host [
4
], a message further underscored at the 70th WHA in 2017 [
5
]. In response, the WHO has refocused on snail control, a strategy that was largely abandoned in 1985 [
6
]and has recently published guidelines for laboratory and field testing of molluscicides in the control of schistosomiasis [
7
].
P.R. China was once one of the most affected countries with an estimated 11.8 million cases of schistosomiasis in the mid-1950s [
8
]. Meanwhile, more than 95% of endemic counties have achieved transmission interruption (TI) or elimination in P.R. China [
9
]. Two hundred and one out of 453 endemic counties had become free of
Oncomelania spp. snails and the nationwide snail-infested areas declined from 8.2 billion m
2 in the mid-1950s to 3.6 billion m
2 in 2020 [
9
,
10
]. The simultaneous decline in schistosomiasis prevalence and snail-infested areas implies that snail control has played an essential role.
Snail control in P.R. China can be stratified into chemical control using molluscicides, environmental modification as well as physical removal [
11
]. Although the application of molluscicides was able to reduce the snail density, and hence, the risk of schistosomiasis transmission, snail elimination was rarely achieved by molluscicides alone [
7
,
8
]. Of note,
Oncomelania hupensis are amphibious. The female snails lay eggs into wet soil that is washed into freshwater bodies, often during flooding events, where they hatch. Schistosomiasis transmission is therefore sensitive to hydrological events. Water conservancy projects against flooding can change the hydrological conditions through attenuation of fluctuation of water levels through the use of artificial reservoirs, and by reducing runoff downstream through the construction of dams. Both of these measures have a devastating impact on snail populations [
12
,
13
,
14
]. Development of agricultural irrigation systems is also playing an important role in governing the distribution of
O. hupensis, as artificial barriers were added to original water networks [
15
]. Other environmental modifications include replacing old irrigation ditches with new cemented canals; construction of water storages for fishing and irrigation, change of land use, among others [
11
].
The Yangtze River basin, which is the major endemic region of schistosomiasis, is home to many water conservancy developments [
16
]. In 1949-1957, embankments were strengthened along the Yangtze River to deter flooding and construction of flood diversion and storage systems were undertaken. Embankment to a large extent blocked the spread of
O. hupensis from the river and lakes to irrigated region, which facilitated the elimination of snails in the latter. Over the next 20 years, from 1958 to 1977, many small reservoirs and irrigation systems abruptly emerged to ensure the water supply for agriculture, accompanied by land reclamation. The snails disappeared soon after dam constructions [
13
,
14
], which was in contrast to the impacts of dam construction on most other
Schistosoma species [
1
]. The fastest elimination of snail habitats occurred in the 1960s and 1970s [
17
], which coincided with the development of agriculture and water conservancy [
18
,
19
].
In this paper, we present results from a nationwide snail survey conducted between March 2016 and December 2017. We mapped the spatiotemporal dynamics of
O. hupensis distribution and explored the relationship between snail control and schistosomiasis elimination at a national scale for the first time.
Methods
Study area and habitat definitions
All counties (
n= 453) historically known to be endemic for schistosomiasis in 12 provincial-level administrative divisions (PLADs) of P.R. China were included in the study (
Figure S1
). A habitat is defined as a relatively independent physical environment infested with
O. hupensis, with adjacent habitats separated by obvious barriers, or their ecological features are essentially different. A habitat can be an artificial man-made environment (e.g., irrigation ditch or channel), or naturally formed (e.g., marshland). Since many habitats remained from the 1950s when the national schistosomiasis control programme was launched, we accepted conventional definitions for the purpose of maintaining the integrity of data. We defined an extinct habitat when
O. hupensis were not found in two or more consecutive years. Otherwise, the habitat was considered as extant.
A habitat is not necessarily fully infested by
O. hupensis. To consider this aspect, we introduced the term "snail-infested range" (SIR) to indicate the actual distribution of snails in the habitat. Hence, a SIR can be less or equal to the overall area. In addition, since the actual distribution of
O. hupensis in a habitat may change or shift from one year to another, the term "accumulated snail-infected range" (ASR) was employed to indicate the maximum distribution range in the habitat by overlapping annual actual distribution ranges from the time that
O. hupensis had been discovered in the habitat.
Data collection and analysis
The national protocol of snail survey was prepared by the National Health Commission and the National Institute of Parasitic Diseases (NIPD) at Chinese Center for Disease Control and Prevention (China CDC). National training sessions were held with the principal investigators at the province level and they in turn trained the staff from local official centres for disease control and prevention (local CDCs) or institutes of schistosomiasis control (ISCs). The habitat survey was performed by skilled staff. The principal investigators sampled and visited 20% of the counties for quality control, (i) to check data completeness in database, and (ii) to confirm the status (extant or extinct) of habitats.
All habitats documented by annual records were registered in the first step of this investigation. The annual records could date back the 1950s and kept in the local agencies of schistosomiasis control. Each habitat was coded by a unique 13-digit identification number that included 2 digits each for each spatial step (PLAD, city, county, township and community) plus a 3-digit serial number. The information pertaining to the habitats were extracted from annual records, including the setting they came from (i.e. marshland, water network or mountainous/hilly landscape), habitat type (e.g., ditch, pond, marshland, etc.), habitat size (m
2), year for first discovery and year of elimination of
O. hupensis, initial SIR when snails were first discovered (m
2), ASR (m
2) and extant SIR (m
2). Data were entered into a Microsoft Excel (Microsoft Corp., Redmond, WA, USA) spreadsheet, doublechecked at the county level and pooled into a Microsoft Access (Microsoft Corp., Redmond, WA, USA) database at the national level. Statistical analysis was performed by SPSS version 19.0 (IBM Corp., Armonk, NY, USA).
The spatial habitat data were obtained by global positioning system (GPS) and organized by a geographical information system (GIS). Briefly, small habitats were drawn up by tracing with hand-held GPS in WGS1984. For larger habitats, the coordinates of four or more key boundary points were obtained by GPS with habitat shapes determined in GIS based on remotely sensed images (
Figure S2
). Each habitat shape file was referred to by its unique 13-digit ID. Habitats destroyed by a land use change approach were mapped according to previous annual records by senior staff. The attribute and shape data of habitat were linked by their unique ID in ArcGIS version 10.1 (ESRI, Redlands, CA, USA).
The occurrence of
O. hupensis was confirmed by standardized sampling procedure, namely systematic sampling and environmental sampling. The description of the methods has been published elsewhere [
20
].
The elevation data covering the study area were obtained from a digital elevation model (DEM; spatial resolution: 90 × 90 m) of Shuttle Radar Topography Mission (SRTM) (http://srtm.csi.cgiar.org/srtmdata/). A grid with slope data was produced based on DEM and then a merged spatial shape without steep depressions across the study area was created in ArcGIS 10.1. The habitats could then be plotted in three dimensions (3-D) by ASR, elevation and slope in Origin version 9.1 (Origin-Lab Corp., Northampton, MA, USA). We explored habitat clustering at the elevation level. We first divided the elevation range into subgroups every 10 m from 0 m. The frequency of habitats in each subgroup was calculated. A cluster was consist of 10 or more consecutive subgroups with a frequency of more than 250 habitats.
We also noted information on TI in the endemic counties, defined as absence of infection in humans, animals or
O. hupensis for at least 5 years [
21
]. The ratio of eliminated ASR was defined as the proportion of the eliminated snail-infested areas in the ASR by the year when a county achieved TI criteria or by 2016. The relation between TI and the ratio of eliminated ASR was analysed. Comparison between TI counties and non-TI counties using non-parametric test (Mann-Whitney) with a significance level of 0.05 was performed in SPSS.
Discussion
In P.R. China, a major feature of successful schistosomiasis control is the shrinking of snail-infested areas. Indeed, the reduction in ASR by over 99% in five provinces has helped to maintain the TI status for more than 20 years. Our results also suggest that the endemic counties with smaller ASR and/or a larger ratio of eliminated SIR are more likely to achieve TI. Such conclusions are corroborated by data from Japan, where
S. japonicum was endemic, but the disease was eliminated in 1976 [
5
]. Of note,
Oncomelania spp. snails once occurred in six separate areas in Japan. Snail control served as the major approach to interrupt the transmission of schistosomiasis. Only few foci of the Japanese snail intermediate host,
O. h. nosophora, remain extant, and these lie in two different endemic areas (Kofu Basin and Obitsu) [
22
], which maintains the necessity for continued attention to schistosomiasis control in the country.
In Indonesia, the overall prevalence of schistosomiasis in humans, rats and snails was 2.6, 8.6 and 2.4% between 2008 and 2011, respectively [
23
]. Although more than half of old foci of
O. h. lindoensis disappeared in some areas, several new foci were observed in conventional endemic valleys (Lindu and Napu) and even new endemic areas (Bada valley) was recently identified [
23
,
24
]. The national baseline survey in the Philippines between 2005 and 2008 indicated that the prevalence of
S. japonicum ranged from 0.1 to 6.3% in 23 provinces [
25
]. Although snail control is included as an additional measure in schistosomiasis control in the two countries, no pronounced effect was observed [
26
,
27
]. WHO is reinforcing snail control as part of its strategic approach to achieve the target of eliminating schistosomiasis as a public health problem by 2030 [
2
,
6
]. Snail control can be implemented where there is high prevalence.
The achievement of schistosomiasis control in P.R. China is attributed to intersectoral collaboration. Although application of chemical molluscicides conducted by the public health sector reduce snail densities, and hence, effectively prevent infections, elimination of
Oncomelania spp. snails by means of environmental modification is more effective in a long term [
11
]. However, schistosomiasis control programmes for sectors other than public health were terminated in recent years due to very low prevalence [
28
], which may be a potential trigger of resurgence of
O. hupensis. Hence, collaboration across sectors (i.e., water conservancy, agriculture, land source management, forestry and public health), coordinated by the central and local governments, plays an essential role in the national snail control programme leading to reduced schistosomiasis transmission [
29
].
The snail habitats in the marshlands are the most difficult ones to deal with. The huge areas make it impossible to eliminate
Oncomelania spp. snails by conventional measures (e.g., application of molluscicides, land reclamation or small water conservancy projects). However, recent evidence indicates that the construction of the Three Gorges dam contributed to decline of snail density in the marshlands. Both snail density and snail-infested areas in previous habitats are declining in the whole middle-lower reach of the Yangtze River [
30
,
31
,
32
]. Nevertheless, the lower water level has led to emerging new marshland ranges with the potential of becoming areas of new
O. hupensis habitats [
33
]. Long-term effects are being monitored. On the other hand, infections in wild animals in these areas of marshlands are not taken into account by the national schistosomiasis control programme, which is making more challenges in reaching the elimination targets in P.R. China [
34
,
35
,
36
]. Boatmen or fishermen there have frequent contact with water in the marshland habitats, which thus might become the bridge between natural transmission in the marshlands and other sectors of the population [
37
]. Taken together, the success of schistosomiasis elimination in P.R. China depends on the effectiveness of snail control in the marshlands. Hence, the Chinese Government should update the strategy for schistosomiasis control in middle-lower reach of the Yangtze River.
Although the present nationwide survey involved all agencies of schistosomiasis control, and quality control guidelines were adhered to, there are still some shortcomings. First, we did not collect the year of resurgence for those habitats in which snails were found to have reappeared. Resurgence of snails in historically endemic areas occurred frequently, which challenges the consolidation of schistosomiasis transmission interruption [
38
]. In order to respond to the resurgence, the counties where schistosomiasis was endemic historically should maintain the capacity of surveillance and response by routine trainings. Second, missing habitats and data were inevitable since some habitats were eliminated very early and the information was not complete. There are 8199 habitats with unknown year of discovery, accounting for 1.1% ASR. Third, we followed the conventional system for classification of habitat landscape type, as used throughout much of China’s snail control history. However, some habitats around lakes might be classified into water network rather than marshland because the original marshland had been transformed into agricultural land. Such habitats are more similar to those in plain water network, as is commonly observed in the eastern part of P.R. China. Finally, we used a constant ASR of snail habitat to replace the annual real snail-infested area as the indicator for temporal pattern analysis, which might lead to an overestimation of the annual nationwide ASR of extant habitats.