Assessment of control strategies against Clonorchis sinensis infection based on a multi-group dynamic transmission model

Clonorchiasis is one of the most important food-borne trematodiases affecting millions of people. Strategies were recommended by different organizations and control programmes were implemented but mostly in short-time periods. It’s important to assess the long-term benefits and sustainability of possible control strategies on morbidity control of the disease. We developed a multi-group transmission model to describe the dynamics of C. sinensis transmission among different groups of people with different raw-fish-consumption behaviors, based on which, a full model with interventions was proposed and three common control measures (i.e., preventive chemotherapy, information, education, and communication (IEC) and environmental modification) and their possible combinations were considered. Under a typical setting of C. sinensis transmission, we simulated interventions according to different strategies and with a series of values of intervention parameters. We found that combinations of measures were much beneficial than those singly applied; higher coverages of measures had better effects; and strategies targeted on whole population performed better than that on at-risk population with raw-fish-consumption behaviors. The strategy recommended by the government of Guangdong Province, China shows good and sustainable effects, under which, the infection control (with human prevalence <5%) could be achieved within 7.84 years (95% CI: 5.78–12.16 years) in our study setting (with original observed prevalence 33.67%). Several sustainable strategies were provided, which could lead to infection control within 10 years. This study makes the effort to quantitatively assess the long-term effects of possible control strategies against C. sinensis infection under a typical transmission setting, with application of a multi-group dynamic transmission model. The proposed model is easily facilitated with other transmission settings and the simulation outputs provide useful information to support the decision-making of control strategies on clonorchiasis.


Introduction
Clonorchiasis is one of the most important food-borne trematodiases, caused by infection with Clonorchis sinensis [1]. About 15 million people are infected with C. sinensis worldwide, predominantly located in East Asia and over 85% of infected cases happened in China [2]. Chronic infection of C. sinensis can lead to complications in liver and biliary systems. Furthermore, C. sinensis is classified as a defined carcinogenic to humans [2,3]. The global disease burden was estimated to 522,863 disability-adjusted life years (DALYs) by WHO in 2010 [4]. Therefore, control of C. sinensis infection becomes one of the major health challenges in endemic society.
The life cycle of C. sinensis involves snails as the first intermediate hosts, fish or shrimp as the second, and humans or other piscivorous mammals as the definitive hosts [5]. The most common way of human getting infected with metacercariae of C. sinensis is eating the raw or insufficiently cooked infected fish/shrimp, the frequency and intensity of which can directly influence the probability of getting infection [6]. Moreover, people can get infected without reporting raw-fish consumption, probably due to poor hygiene habits, e.g., eating food contaminated by metacercariae of C. sinensis through tableware or chopping board previously used to prepare raw fish without effective cleaning [6,7]. Different behaviors of human in rawfish consumption may show different transmission characteristics, which is important particularly when studying the dynamics of the transmission process.
The major control measures of clonorchiasis includes preventive chemotherapy, information, education, and communication (IEC), environmental modification and the possible combinations of the above [8]. Chemotherapy is the most important measure for the management of clonorchiasis, the targeted population of which can be all residents (i.e., mass drug administration, MDA) or selective population (e.g., infected individuals and at-risk groups with raw-fish-eating behavior) [8,9]. Repeated MDA or chemotherapy with selective population showed good morbidity control effects during short periods (e.g., 2 or 3 years) [9]. IEC combined with chemotherapy is suggested to enhance the control sustainability [10]. And environmental modification, an additional way of control measure by removing unimproved toilets built adjacent to fish ponds, can prevent the contamination of water by faeces and decrease the number of infected snails and fish [11].
Effective and sustainable control strategies on clonorchiasis are important, assisting local governments and policy makers to implement the interventions. WHO recommended a control strategy of clonorchiasis that advocates praziquantel administration for all residents every year in high endemic areas (with prevalence�20%) and for all residents every two years or individuals regularly eating raw fish every year in moderate endemic areas (with prevalence <20%) [12]. The Chinese government set a goal in 2015 to reduce prevalence of clonorchiasis by 30% or more in areas with prevalence higher than 10% by 2020. Intervention strategy was developed, including chemotherapy every year for at-risk population with a coverage by 80%, improvement of sanitation toilets by 80% and increasing the awareness of key knowledge on clonorchiasis control among students by 95% [13]. The government of Guangdong Province, one of the most endemic provinces in China, set a more ambitious goal to reduce the prevalence by 40% in areas with prevalence higher than 20%. Stronger interventions were recommended, including coverage of chemotherapy by 80% on at-risk population every year, coverage of sanitation toilets by 90% and improvement of people's healthy behaviors by 90% [14]. However, studies were absence on evaluation of the long-term (>10 years) effects of the above strategies. Several control programmes have been implemented in endemic areas, focusing on either chemotherapy only or interventions combined with several measures [9,[15][16][17]. But the study periods were short, mostly less than 5 years, thus it's difficult to assess the longterm effects and the sustainability. There is an urgent need to explore long-term benefits and effectiveness of different potential strategies on morbidity control of clonorchiasis within a quantitative framework.
Mathematical models play important roles for understanding the transmission dynamics of parasitic diseases [18]. Furthermore, they test hypotheses about likely dynamics and epidemiology under proposed strategies, thus support control and elimination programmes [19]. Four papers have been published on mathematical models depicting the transmission dynamics of C. sinensis infection. Dai and colleagues proposed a deterministic SI compartmental model describing the transmission of C. sinensis infection between human, snail, and fish. The authors also verified the globally asymptotically stable of the model in both disease-free and endemic equilibria [20]. Yuan and colleagues applied a deterministic SEIR and a SIR model to describe the spread of clonorchiasis in Foshan and Guangzhou, China, respectively [21,22]. Compared to Dai's model, the models developed by Yuan included recovered population as well as different status (e.g., eggs and cercariae) of the parasite. Gao and colleagues proposed a SIRS model to evaluate the effects of snail control and health education on clonorchiasis control in Foshan, China [23]. These papers analyzed relationships between model parameters and the basic production number (R 0 ), giving insights to possible directions of intervention strategies. However, no further indicators (e.g., time to achieve transmission control or extinction, percentages of prevalence decreased through time), which are more applicable for selection of specific strategies, were used for assessment. Furthermore, the models assumed homogeneity of human population, yet people with different behaviors of raw-fish consumption show different transmission patterns, which should be considered when the transmission dynamics is depicted. And control strategies focus on population with raw fish consumption instead of whole population cannot be assessed under these models. Interestingly, Bürli and colleges [24,25] developed deterministic models of Opisthorchis viverrini transmission dynamics to evaluate the long-term effectiveness of different interventions (e.g., MDA, IEC, and environmental modification), by including the parameters of interventions inside the models and simulating different scenarios. O. viverrini, known as the other major species of liver fluke, have very similar life cycle as C. sinensis [26].
In this paper, with extension of Dai's transmission model and reference of Bürli's intervention model [20,24], we developed a multi-group deterministic SIS transmission model to: (1) describing the transmission dynamics among groups of population with different raw-fishconsumption behaviors and (2) quantitatively assessing the long-term effectiveness of different intervention strategies on control of C. sinensis infection through simulating scenarios. As a practical example, the model was applied in Fusha Town, a typical clonorchiasis endemic area with raw-fish-eating culture in Guangdong province, China. The model parameters were estimated by fitting the model with the epidemiological survey data collected in the town.

Method overview
We developed a basic model to describe the transmission of C. sinensis, with consideration of different groups of population with different raw-fish-consumption behaviors. Using the survey data in Fusha Town, Guangdong Province, we estimated the unknown model parameters under Bayesian melding approach. Validation and sensitivity analysis were undertaken afterwards. A full model with interventions were further developed by adapting the basic model with specific intervention parameters, and simulations were carried out to evaluate the longterm effects of different control strategies. In particular, we assessed: (1) the current strategies recommended by WHO, the Chinese government and the government of Guangdong Province; and (2) the three commonly used measures (i.e., preventive chemotherapy, IEC and environmental modification) when applied singly or combinedly, with different coverages, durations and frequencies of chemotherapy and improvement rates of IEC and environmental modification. The resulted sustainable intervention strategies were listed for practical application. A flow diagram of the methodology is show in Fig 1.

The basic model
As individuals previous infected with C. sinensis show no apparent protective immunity after recovery from treatment [27] and re-infection is common in endemic areas [9,28], we proposed a deterministic Susceptible-Infected-Susceptible (SIS) model to describe the transmission of C. sinensis among human-snail-fish populations, where infected human recovered from infection will instantly return to the susceptible state. According to different frequencies of raw-fish consumption, we classified the human population into four groups, that is seldom, moderately (less than 5 times per year), often (5-10 times per year) and very often (more than 10 times per year) eating of raw fish. We assumed each group of people get infection under different transmission rates. Fig 2A describes the model structure graphically, and Tables 1 and 2 list the descriptions of state variables and model parameters, respectively. We assumed each host has two states (i.e., susceptible S and infectious I) and the transmission dynamics were modeled by a system of twelve ordinary differential equations. More detailed description of the model is shown in S1 Text.
The basic reproduction number R 0 (i.e., the average number of new cases generated by a typical infected individual in an otherwise susceptible population) was derived by the next generation matrix approach (with details in S3 Text). The expression for R 0 is as following:

Data
As a practical example, we applied the model in Fusha Town (coordinates: 22˚40'2.64''N, 113˚20'58.2''E, area: 37 km 2 , and population: 59,593 in the year 2017), located at the northern part of Zhongshan City, Guangdong Province (S1 Fig), where C. sinensis infection is highly endemic [5]. Situated among the Pearl River Delta's vast water network, the town is crisscrossed by rivers and the fishing industry is flourishing. Local residents like eating raw fish and have poor hygiene habits of using the same chopping blocks for both raw and cooked food [33]. All these geographical, environmental and socioeconomic factors make the town a typical place for the disease transmission. The prevalence data of human was obtained from a cross-sectional study done between October 2011 and June 2012 in the town, the detailed result of which is shown in Table 3 [33]. As no prevalence data of fish is available for Fusha town, we used the data from a cross-sectional survey between October 2008 to September 2009 in Zhongshan City instead, where the town belongs. The survey shows the prevalence of fish infected with cercariae of C. sinensis is 20.97% (39/186) [34]. Even though no survey data of snails is available for either Fusha Town or Zhongshan City, series of studies reported the prevalence of infected snail between 0.1% and 3.7% in different endemic areas of Guangdong province, with study period between 2004 to 2009 [34,35]. This interval was set as prior information for the following parameter estimation.

Parameter estimation
In absence of time-series longitudinal data, we assumed the system of basic model on the endemic equilibrium point under the observed situation. In order to capture the local transmission dynamics, we applied Bayesian melding approach to estimate the unknown parameters based on the basic model and the survey data [36]. More details were described in the S4 Text.
Validation was undertaken by calculating the modeling efficiency statistic (EF) [37] and the coverage of observed data within 95% credible intervals of posterior estimations. Here, O is the mean of the observed values, and E i is the estimated value from the best set of estimated parameters. Chi-square test was used to see whether there is significant difference between observed and estimated prevalence [37,38]. Local sensitivity index [25] and partial rank correlation coefficients (PRCC) [39] were calculated to assess the influence of different parameters on R 0 , locally and globally, respectively.

Full model with interventions
We considered common types of interventions against C. sinensis infection, that is preventive chemotherapy, IEC, environmental modification, and possible combinations of the above measures. In particular, two kinds of improvement through IEC were considered: (1) improvement of hygiene habits and (2) changing the behavior of raw-fish-consumption. The improvements through IEC and environmental modification were assumed continuous, that is, once it happens, it will last through time in the same population. Intervention parameters were set as following: γ 2,i the recovery rate of infected humans through preventive chemotherapy, C m,i the coverage of chemotherapy, T i the time interval of treatment, C e1,i the improvement rate of hygiene habits, C e2,i the rate of stopping raw-fish-eating behavior and C d the coverage of sanitation toilets. i = 1,2,3,4 represent human groups who seldom, moderately, often and very often consume raw fish, respectively. A full model is show in Fig 2B, which describes transmission dynamics of C. sinensis infection under different control interventions. More details are shown in S2 Text. The approach for calculation of control reproduction number (R c ) in the system with interventions is similar as that of R 0 . With the assumption that chemotherapy is distributed continuously [24], R c is expressed as following: ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffiffi

Simulation of interventions
A baseline scenario was initiated at the endemic equilibrium state according to the situation of Fusha Town, which presents a typical high endemic setting. To access the current recommended strategies, we simulated the corresponding interventions recommended by WHO, the The increased transmission rate from an infected fish to a susceptible human who moderately consumes raw fish, compared to that who seldom 3 The increased transmission rate from an infected fish to a susceptible human who often consumes raw fish, compared to that who seldom The increased transmission rate from an infected fish to a susceptible human who very often consumes raw fish, compared to that who seldom Chinese government and the government of Guangdong Province, respectively, with an intervention duration of 5 years and a further post-intervention period of 45 years. The values set for the corresponding intervention parameters are listed in Table 4. In addition, to explore the effects of intervention durations, coverage rates and frequencies of chemotherapy, improvement rates of the two kinds of healthy behaviors by IEC, improvement rates of environmental modification, and possible combinations of the three measures, we set a series of values for intervention parameters and simulated the corresponding scenarios. Indicators including control reproduction number R c , prevalence and its reduced rates in 5, 10 and 15 years from the beginning of intervention, were calculated to evaluate the corresponding interventions. We considered the intervention sustainable if R c <1 and the prevalence in the post-intervention period remain decreasing. In addition, with reference to the schistosomiasis control criteria in China [40], we defined areas with clonorchiasis prevalence less than 5% as infection control areas, and less than 1% as transmission control. We calculated the years from the beginning of intervention to infection control and transmission control status, respectively, for each scenario. We assumed the drug efficacy h to be 92% according to the median of the corresponding studies [16,[41][42][43][44][45].

Parameter estimation
The estimated results of unknown parameters are shown in Table 2. The transmission rates were much higher for people who regularly eat raw fish than that for people who seldom. Particularly, the rates were 7.5 (95% CI: 6.0-10.4), 25

The three common measures when applied singly
Chemotherapy, when applied singly, were able to reduce human prevalence fast during the intervention, but it was not sustainable (S3 Table). Of noted, it performed better on whole population than at-risk population with raw-fish-consumption behaviors under the same frequencies, coverage rates and duration time of intervention (S4 Fig and S3 Table). Higher frequencies of treatment, higher coverage rates and longer duration of intervention could result in larger reduction of human prevalence at the end of intervention. However, even with twice treatment each year on whole population under 100% coverage and 20 years' duration of intervention, the prevalence would still rebound after intervention stopped. Under the assumption of continuation of the improvements, the IEC measure was sustainable when applied singly, if both kinds of improvement rates reached 93.29% (95% CI: 91.87%-94.06%) on whole population (S5 Fig and S4 Table). The effects increased when the improvement rates rose. To be noted, IEC only focus on improvement of hygiene habits of people or changing the behavior of raw-fish-consumption showed less effective than that focus on both. For example, even people achieve a 100% improvement of hygiene habits, the prevalence would remain at high level (22.08%, 95% CI: 19.50%-24.46%) after 20 years from the beginning of intervention, if no improvement was undertaken in changing people's behavior of raw-fish-consumption. Similarly, the prevalence was not able to decrease below 9.33% (95% CI: 7.07%-14.16%) within 50 years without improvement of hygiene habits, even when 100% of at-risk population with raw-fish-consumption behaviors would not consume raw fish any more. We found the effects of IEC slow, as even with a 90% improvement of both kinds of healthy behaviors, the study area needed more than 33.87 years (95% CI: 24.65-50.27 years) to reach infection control, and was not able to reach transmission control within 50 years. Similar

Measures when applied combinedly
Generally, interventions combined with two or three measures performed much better than that only applied singly, and higher coverages resulted in larger reductions of prevalence (S6 and S7 Tables and Fig 5). Table 6 lists several selected sustainable combinations, under which the study area could achieve infection control within 10 years from the beginning of intervention. With a high coverage (i.e., 90%) of chemotherapy on whole population (once per year for five years), the improvement of people's healthy behaviors should reach 94%, or the improvement of people's healthy behaviors and the coverage of sanitation toilets should both reach 70% to achieve the goal. If chemotherapy only focus on at-risk population with raw-fish-consumption behaviors, a higher improvement (85%) of both people's healthy behaviors and coverage of sanitation toilets should be applied. Under a high coverage (90%) of sanitation toilets  and high improvement (90%) of people's healthy behaviors, chemotherapy (once per year for five years) with 20% coverage on whole population or 40% on at-risk groups with raw-fishconsumption behaviors could result in transmission control within 10 years. In addition, a coverage of chemotherapy targeted on whole population by 60% combined with improvements of people's healthy behaviors and environmental modification by 70%, or a coverage of chemotherapy on at-risk population by 80% combined with improvements of people's healthy behaviors and environmental modification by 85%, could also achieve the goal.

Advantages of the methodology
In this study, we developed a multi-group transmission model to describe the dynamics of C. sinensis transmission among people with different raw-fish-consumption behaviors. Based on the above model, we proposed a comprehensive model with interventions and simulated the long-term dynamics under various control strategies. The intervention model also allows us to investigate the control effectiveness among different raw-fish-consumption human population and with different improvements of IEC. To our knowledge, it is the first time that a multi-group model was proposed for transmission dynamics of C. sinensis infection. Values of many parameter could not be obtained from existing studies. Data-driven approaches, integrating survey data into models to estimate parameters that capture the local transmission dynamics [46], are most suitable approaches of parameter estimation for our study setting. Bayesian melding approach is one of the most common data-driven approaches [36,47,48]. It has important advantages compared to the maximum likelihood approach, including the ability to produce probability distributions of the estimated parameters and model outputs which can be used for inference and projection, and the good computational efficiency when a large number of parameters are estimated simultaneously [49]. The performance of our model estimation seems good, as the outputs were able to capture the 100% of the observed data within 95% CI, the EF score was close to 1, and the Chi-square test show there was no significant difference between observed and estimated prevalence.

Three recommended strategies
Based on a typical setting of C. sinensis transmission (Fusha Town), we simulated scenarios according to strategies recommended by WHO, the Chinese government and the government of Guangdong Province, with intervention duration of five years. On one hand, the WHO strategy, mainly focusing on chemotherapy of whole population, could decrease the human prevalence rapidly but were not sustainable (Fig 4), even with an extension of intervention duration to 30 years (S4 Fig), which may due to the high reinfection of the disease [9,28]. On the other hand, strategies recommended by the Chinese government and the government of Guangdong Province could reach their setting goals within 5 years of intervention (Fig 4). However, the former strategy, combined with chemotherapy and environmental modification, showed slower effects and human prevalence won't reduce to the transmission control level within 50 years. This effect might be slightly underestimated, as we ignored one component of the strategy (i.e., increasing the awareness of key knowledge on clonorchiasis control among students through IEC by 95%). Students are not among the high-risk groups of transmission [50], and it's difficult to assess how much improvement of healthy behaviors the increase of knowledge will result in. The latter strategy, even though showed the best control effects, is quite strong and ambitious, as a 90% coverage should be reached for the improvements of both environmental modification and healthy behaviors, and an 80% coverage for chemotherapy on at-risk population, which may be hard to achieve in practice, especially the one for  � Results were expressed based on the best set of parameter estimates and their 95% CI; the duration of chemotherapy was 5 years; C d , C e and C m indicate the coverage of sanitation toilets, the improvement rate of both kinds of healthy behaviors by IEC, and the coverage of chemotherapy on targeted population, respectively; chemotherapy is applied once per year for five years; R c −Y 1% have the same meanings as that in Table 5. # The targeted population of chemotherapy; "at-risk" means at-risk population with raw-fish-consumption behaviors. https://doi.org/10.1371/journal.pntd.0008152.t006 improvements of people's behaviors [51]. From a practical implementation perspective, we might overestimate this strategy.

Measures when applied singly or combinedly
In order to find more applicable strategies, we set a series of values for intervention parameters and assessed the corresponding simulated results. We found that when applied singly, IEC or environmental modification showed slow and moderate control effects but kept decreasing the prevalence of C. sinensis infection within 50 years, while chemotherapy indicated much quicker and stronger effects but prevalence was easy to rebound after intervention stopped. Only with very strong chemotherapy (e.g., 5 times per year for 33 years, or 2 times per year for 70 years) on whole population could lead to sustainable transmission control (S4 Fig). Besides, IEC only focus on improvement of hygiene habits of people or changing the behavior of rawfish-consumption showed less effective than that focus on both. The higher coverages of measures resulted in better effects, which is similar as previous studies on other parasitic diseases [52]. Combinations of measures were much beneficial than those singly applied; and particularly, chemotherapy combined with both IEC and environmental modification showed much better effects than that combined with only one measure. Chemotherapy could reduce the number of infected people but could not prevent reinfections or new infections [9], while IEC helps to reduce the transmission rate from infected fish to susceptible human [25], and environmental modification decreases the chance of the transmission from infected human to snails [25]. Thus, all of the three measures implemented together will effectively control the transmission, which is strongly recommended for morbidity control of the disease. We also listed several sustainable combined strategies leading to infection control within 10 years in our current setting, with chemotherapy on either whole or at-risk population under coverages from low to high (Table 6). To be noted, durations, levels of coverage and kinds of targeted population (e.g., whole or selective) may influence the compliance of chemotherapy [9,17,53], thus both simulated effects and compliance in practice should be considered before selection of strategies.

Insights from sensitivity analysis
Sensitivity analysis, characterizing the response of model outputs to parameter variation, shows strength and direction of the influence of different parameters on the outcomes of interest, thus provides insights on possible intervention strategies [54]. Results from sensitivity analysis (Fig 3) should be considered together with simulated outputs of transmission models and actual practical challenges, for finding better interventions. In our basic transmission model, the death rates (μ h , μ s and μ f ), the transmission rates (β h,1 , β s and β f ), the recruitments of fish and snail (λ s and λ f ), and the basic recovery rate γ 1 showed relatively strong influence on R 0 . The three common measures, that is chemotherapy, IEC and environmental modification, could influence recovery rate, β h,1 and β f , respectively, which is quite reasonable to consider for interventions. Application of molluscicides may decrease λ s and increase μ s . However, due to the wide distribution of snails, their low prevalence of C. sinensis infection, and the potential toxic effects of molluscicides to fish, this measure is not recommended in endemic areas [8,55]. β f could be decreased by developing vaccines targeting the intermediate hosts, which is still under investigation [56]. Interestingly, parameters related to the group of people who very often consume raw fish have much stronger influence on R 0 compared to that of other human groups, but chemotherapy or IEC focusing only on this group showed much softer effects (S3 and S4 Tables), which may be explained by the small proportion of population in the group. On one hand, people who very often consume raw fish play an important role in transmission and need to pay more special attention to. On the other hand, measures should be taken on wider range of people instead of just on this particular group for a better morbidity control of the disease.

Limitations
This study is subject to some limitations. First of all, we made simplifying assumptions to several potential heterogeneities. We assumed IEC and environmental modification equally efficacious and continuous over the simulation period; however, the effects may change along with social and economic development, while it's difficult to assess how those effects might change. We ignored the impact of seasonality on the transmission, whilst temperature and rainfall, showing close relationships with seasonality, may influence the egg hatching, activities and distributions of intermediate hosts [57], thus affect the transmission of the disease. As our studied area is relatively independent from surroundings and the population are stable (S1 Table), we assumed the transmission of the disease was influenced very little by nearby areas and did not considered the interactions of disease transmission between areas. However, the distance between Fusha Town and Zhangshan City is quite close (around 29 km) and migration of population may happen, which could have an impact on the transmission. We didn't consider diversity of compliance with preventive chemotherapy on different groups of people, however, people who seldom consume raw fish may show less compliant than those very often consume. We assumed the basic transmission rate β h,1 influenced mainly by improvement of hygiene habits and ignored the effect of changing raw-fish-eating behavior of people. The latter may decrease the chance of tableware or chopping board being contaminated by metacercariae of C. sinensis as fewer people would prepare the raw-fish-dish, thus may reduce β h,1 . We may underestimate the effect of this improvement. Due to unavailability of the corresponding data, we didn't address the potential heterogeneities mentioned above to avoid large uncertainties, as the performance of models with increased complexity is critically conditional on an increase of available data for parameterization [58]. Besides, we assumed human the major definitive host and ignored the contributions of reservoir hosts (e.g., dogs and cats), which may lead to slightly over-estimation of effects of the studied strategies, as they mainly focus on human. We made the parsimonious assumption that the system was at an endemic equilibrium and fitted and validated the basic model using the cross-sectional survey data, based on which parameters pertaining to the transmission process were estimated. However, we could not validate the intervention model due to lack of longitudinal data on different intervention scenarios. Preventive chemotherapy recommended by WHO is praziquantel at a dosage of 40mg/kg in one single day [45], with drug efficacy of 87% [43], while control programs in China applied different modes, commonly with praziquantel at a dosage of 75mg/kg/day for 2 days [16] or albendazole 0.8g/day for 4 days [44], the corresponding efficacies of which were around 98%-100% or 83%-92%, respectively [41,42]. Different administration modes of chemotherapy may result in different efficacies and for convenience, we used the median efficacy (92%) of the above studies. To be noted, the simulation results may have minor changes with varied drug efficacies (e.g., the simulated prevalence after five years of preventive chemotherapy under the WHO strategy would increase or decrease by around 5% if drug efficacy was set to 87% or 98%, respectively). In addition, we assessed different intervention strategies under conditions on Fusha Town, a typical setting of C. sinensis transmission in Guangdong Province, China, the results of which can be generalized to other endemic areas with similar epidemic and under similar socioeconomic and environmental conditions. However, for areas far away the generalization should be taken very cautiously. Nevertheless, the multi-group model can be facilitated with epidemiological data from those areas and produce assessment results fit to their local dynamics. Besides, this study didn't consider economic factors for implementation of the interventions. We are currently extending this work to cost-effective analysis of different strategies and developing models adjusting spatial disparities to provide more comprehensive and generalized assessments.

Conclusions
In summary, under the typical transmission setting we applied, our simulated results suggested that combinations of measures were much beneficial; higher coverages of measures had better effects; and strategies targeted on whole population performed better than that on at-risk population with raw-fish-consumption behaviors. The strategy of Guangdong Province performed best among the three recommended ones. Sustainable strategies, such as combined measures of preventive chemotherapy (once per year for five years) with 60% coverage on whole population, IEC and environmental modification with improvement rates of 70%, or combined measures of the above preventive chemotherapy with 80% coverage on at-risk population, IEC and environmental modification with improvement rates of 85%, could achieve infection control within 10 years. In conclusion, this study makes the effort to quantitatively assess the longterm effects of possible control strategies against C. sinensis infection. The model we proposed is facilitated to other transmission settings and the simulation outputs provide useful information that supports the selection of control strategies on clonorchiasis by researchers and decision makers.