Next Article in Journal
Analysis of a State Degradation Model and Preventive Maintenance Strategies for Wind Turbine Generators Based on Stochastic Differential Equations
Next Article in Special Issue
Analyzing the Dynamics of a Periodic Typhoid Fever Transmission Model with Imperfect Vaccination
Previous Article in Journal
An Efficient Evolution-Based Technique for Moving Target Search with Unmanned Aircraft Vehicle: Analysis and Validation
Previous Article in Special Issue
Machine-Learning Approach for Risk Estimation and Risk Prediction of the Effect of Climate on Bovine Respiratory Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adapting Strategies for Effective Schistosomiasis Prevention: A Mathematical Modeling Approach

1
Department of Animal Ecology and Systematics, Justus Liebig University Giessen, Heinrich-Buff-Ring 26 (iFZ), 35392 Giessen, Germany
2
Department of Landscape Ecology and Resource Management, Justus Liebig University Giessen, Heinrich-Buff-Ring 26 (iFZ), 35392 Giessen, Germany
3
Bill and Joyce Cummings Institute of Global Health, University of Global Health Equity, Kigali Heights, Plot 772 KG 7 Ave., Kigali P.O. Box 6955, Rwanda
4
Centre for International Development and Environmental Research (ZEU), Justus Liebig University Giessen, Senckenbergstrasse 3, 35390 Giessen, Germany
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2609; https://doi.org/10.3390/math11122609
Submission received: 8 May 2023 / Revised: 30 May 2023 / Accepted: 5 June 2023 / Published: 7 June 2023

Abstract

:
One of the most deadly neglected tropical diseases known to man is schistosomiasis. Understanding how the disease spreads and evaluating the relevant control strategies are key steps in predicting its spread. We propose a mathematical model to evaluate the potential impact of four strategies: chemotherapy, awareness programs, the mechanical removal of snails and molluscicides, and the impact of a change in temperature on different molluscicide performances based on their half-lives and the length of time they persist in contact with target species. The results show that the recruitment rate of humans and the presence of cercaria and miracidia parasites are crucial factors in disease transmission. However, schistosomiasis can be entirely eradicated by combining all of the four strategies. In the face of climate change and molluscicide degradation, the results show that increasing the temperatures and the number of days a molluscicide persists in the environment before it completely degrades decreases the chemically induced mortality rate of snails while increasing the half-life of different molluscicides increases the death rate of snails. Therefore, eradicating schistosomiasis effectively necessitates a comprehensive integration of all preventative measures. Moreover, regions with different weather patterns and seasonal climates need strategies that have been adapted in terms of the appropriate molluscicide and time intervals for reapplication and effective schistosomiasis control.

Graphical Abstract

1. Introduction

Schistosomiasis is the second-most significant neglected tropical disease (NTD), a physically debilitating and persistent disease [1,2] that leads to severe morbidity and almost 12,000 deaths globally, of which at least 90% are from sub-Saharan Africa [2]. Human schistosomiasis infection caused by Trematoda worms depends on the availability of suitable freshwater intermediate host snails (IHs) and the final human host to be transmitted [3]. The main symptoms of infection include skin rash and itching, fever, cough, muscle pain, bloody urine, and growth retardation in children [2], while severe cases can lead to damage and failure of the liver, bladder, lungs, and intestines [1,2]. Unfortunately, unlike most NTDs, including lymphatic filariasis, leprosy, and leishmaniasis, there are currently no recommendations for intensive disease management for schistosomiasis [4], and current disease control strategies rely on mass drug administration (MDA) chemotherapy, which, at times, is combined with health education [5,6]. Although there has been a success in the implementation of chemotherapy, especially in Africa [7,8], this may have limited the implementation of new methods to interrupt parasite transmission through snail control, leading to a resurgence of the disease following drug treatment and education efforts [6,9,10]. Therefore, there is a need to recast the current and future schistosomiasis control strategies to focus on integrated measures and enhanced surveillance response to both changes in infection and snail abundance.
Various methods have been employed in the control of schistosomiasis. For instance, intermediate hosts (IHs) can be controlled through biological, chemical/molluscicides, or environmental/ecological methods [11,12,13]. Biological measures include the use of predators and competitor snails [12], while chemical control measures consider the use of molluscicides to reduce IH density and abundance [14,15]. The use of molluscicides and their effectiveness has been observed to depend on their proper concentration, half-life, technique and timing of application, and length of contact time with target species, as well as the temperature of the environment where the chemical is released [16]. In addition, environmental/ecological strategies are used to lower the density and abundance of IHs and the risks of snail-to-human transmission (referred to as the mechanical control approach in this study). The approach includes excavating deep channels that act as dry buriers for snails to spread to other water outlets, increasing the flow velocity in irrigation canals, and picking IHs out of the system. Other ways include adopting appropriate cultivation methods, such as shorter fallow times, modifying irrigation techniques, and regulating flooding [11,17].
With the predicted change in global temperature due to climate change, mathematical models have become valuable tools for exploring various infection and control scenarios as they provide both theoretical and practical insights into the epidemiology of infectious diseases. In understanding the transmission and control dynamics of schistosomiasis, different models have been used to create a range of transmission and control outcomes. These models explore dynamics, including free-living parasites, miracidia, and cercariae [18], chemotherapy treatment using praziquantel drugs and the removal of snails and cercariae [19], the impact of public literacy and snail control parameters [20], as well as health education and molluscicides [21].
Our study examines whether molluscicides and/or mechanical control for snail management, combined with mass drug administration and/or health education, can eliminate schistosomiasis. In addition, we postulate that the temperature rise, as predicted by climate change models [22,23,24], will strongly influence the use of chemicals/molluscicides, and their performance against the target species [25,26] and the general transmission dynamics of the disease. As such, it is important to understand how temperature affects molluscicide usage for the control of schistosomiasis. This can help in the identification and control of chemical-induced environmental hazards that could be harmful to even nontarget species and human health under various environmental conditions [27,28,29,30]. The current study will evaluate four controls: mass drug administration/chemotherapy, public literacy, mechanical measures, and chemicals/molluscicides to find the most effective way to eradicate schistosomiasis. In addition, the evaluation includes the effects of temperature rise on various molluscicide performances based on half-lives and contact durations with the targeted species.

2. Materials and Methods

2.1. Schistosomiasis Model Formulation

The dynamics of schistosomiasis [2,3] and stages of intervention strategies are the basis for model formulation (Figure 1). The deterministic model formulated and represented by ordinary differential equations (ODEs) (Equations (1)–(7)) is a modification of those presented by Abokwara and Madubueze [20] and Nur et al. [21], and it is based on the compartment diagram (Figure 1). A population is represented by state time variables that are associated with the various compartment models and change over time according to the interaction rules. t : the time variable for the sizes of susceptible humans; S h t : infected human; I h t : eggs released from infected humans; E h ( t ) : free-living miracidia; M f ( t ) : susceptible snail vector; S v ( t ) : infected snail vector; I v ( t ) : free-living cercariae C f t (Figure 1). The parameters of the model are based on the abiotic and biotic aspects of the Schistosoma cycle and transmission, namely the contact rates between the hosts and the free-living cercaria/miracidia population in water and the effective mortality rates.
The population of susceptible human hosts S h in the model system, Equation (1), increases by an exponential recruitment rate of Λ 1 e υ 1 x , where Λ 1 is the maximum per capita birth rate/immigration rate (recruitment rate) of human individuals, υ 1 is the natural mortality rate of humans, and x is the initial age of infection in children. Susceptible human hosts, S h , are infected only by contact with cercariae in an infested freshwater environment at a rate β 1 per contact. The saturation incidence is a Holling-type II function, which relates the rate of infection transmission and the inhibitory effect, where C o is the saturation coefficient that represents the transmissibility of infection due to the crowding effect of the cercariae, and ε is the limitation of the growth velocity (density) of the cercaria in contaminated freshwater. Thus, model system, Equation (2) represents the infected humans, I h . Moreover, due to schistosomiasis infection, infected people die at a rate of δ 1 . Thus, our model incorporates chemotherapy treatment for deworming schistosome worms in the infected human body and to stop or limit the release of parasite eggs through urine and feces. We assume that treated infected individuals recover at a rate of ψ and return to the susceptible class because there is no lasting immunity to schistosomiasis, and reinfection is inevitable [9,31]. In addition, public literacy strategy is considered and represented by the parameters α , α h ϵ [ 0 , 1 ] , where α is the proportion of people with knowledge/awareness about schistosomiasis, and α h is the effectiveness of acquired knowledge. This control method focuses on managing the human population by providing access to clean water and improving sanitation and hygiene (WASH) to lower the number of humans being infected [32].
When infected humans I h indiscriminately release parasite eggs through urine/feces into the environment, their population size increases logistically, as shown in Equation (3). Where ρ (in grams) is the portion of stool/urine per infected person due to open defecation or urination. θ h is the number of parasite eggs per gram of stool and/or urine, and K is the carrying capacity of the parasite eggs in the environment. The term 1 α α h represents the public literacy impact in reducing the portion of urine/feces ρ . It describes the impact of public literacy on the density of parasite eggs directly or indirectly deposited in freshwater and lowers the likelihood of human-to-snail transmission. The indirectly deposited parasite eggs find their way into a freshwater source. Each parasite egg E h hatches, releasing N E miracidia per parasite egg at a rate ω 1 under suitable conditions, or dies naturally at a rate υ 3 when there is no IHs to penetrate, as represented in miracidia population dynamics in Equation (4).
In the model system, Equation (5), the snail host population S v is recruited at a rate Λ 2 and is infected due to contact with miracidia at a saturation incidence given the saturation coefficient for miracidia infectivity M o and ε the limitation of the growth velocity of the miracidia. β 2 is the rate of miracidia-to-snail transmission per contact. Both susceptible S v and infected IH snails I v die naturally at a rate υ 2 , but the infected snail (in Equation (6) may also die at a rate of δ 2 due to infection host. Infected snails I v that survive release infectious cercariae C f at a rate of ω 2 , capable of invading and infecting humans, or dying naturally at a rate of u 5 in the absence of a human host, see Equation (7). Thus, we incorporated a snail management strategy through the mechanical measure, represented by parameter b , α v ϵ [ 0 , 1 ] , where b is the proportion of snail density removed and eliminated from the system, and α v is the effectiveness of such snail management practice. Thus, the impact of a mechanical measure 1 b α v reduces β 2 [20,21].
In addition, we assume that the application of molluscicides into the environment causes susceptible and infected IHs to die at a chemical-induced death rate of σ s ϵ [ 0 , 1 ] and σ I ϵ [ 0 , 1 ] , respectively. Furthermore, molluscicides also reduce viability and cause the deaths of infective cercaria and miracidia at a rate of σ m ϵ [ 0 , 1 ] and σ c ϵ [ 0 , 1 ] , respectively.
Thus, the model equations incorporating terms for control strategies highlighted in blue are given as follows:
S h ˙ t = Λ 1 e υ 1 x + ψ I h 1 α α h β 1 S h C f C o + ε C f υ 1 S h ,
I h ˙ t = 1 α α h β 1 S h C f C o + ε C f υ 1 + δ 1 + ψ I h ,
E h ˙ t = ρ θ h 1 α α h I h 1 E h K ω 1 + υ 3 E h ,
M f ˙ t = N E ω 1 E h υ 4 + σ m M f ,
S v ˙ t = Λ 2 1 b α v β 2 M f S v M o + ε M f υ 2 + σ s S v ,
I v ˙ t = 1 b α v β 2 M f S v M o + ε M f υ 2 + δ 2 + σ I I v ,
C f ˙ t = ω 2 I v υ 5 + σ c C f .

2.2. Temperature Control

Temperature is a crucial factor for the timing of molluscicide application [16] because molluscicide efficacy is temperature-dependent [33]. The temperature has a significant effect on how quickly the half-lives of molluscicide decrease in water [27]. This has a great influence on the chemical-induced death rates σ s , σ I , σ m , and σ c and can determine the rate of transmission of the disease. We examine and evaluate the effectiveness of temperature-dependent molluscicides on the death rates of targeted species. We derive the environmentally dependent chemical-induced deaths of species, where the effect of temperature can be reliably predicted using an Arrhenius equation [27,28,34]. The Arrhenius equation proposed and updated by the European Chemicals Agency, ECHA [29,30] to modify the half-life t 1 / 2 T at any temperature T of the environment is given as
t 1 / 2 T = t 1 / 2 T o e 0.08 T o T
where t 1 / 2 T o is the factory-predetermined half-life of the chemical at the experimental temperature T o . The preferred experimental temperature in most of the chemicals is T o = 20 ° C [27,35]. Furthermore, we assume that chemically induced mortalities ( σ s , σ I , σ m , σ c ) decrease exponentially with time t according to the following expressions:
σ s , I , m , c = σ s , I , m , c e k t
where σ s , σ I , σ m , and σ c are the maximum mortality rates on the first day ( t = 0 ) when the molluscicide is applied in the water. k , is the exponential decay constant of the mortality rate at degradation time t . The terms σ s , σ I , σ m and σ c were determined from the withdrawal terms ( d v + σ s ) , d v + δ 2 + σ I , υ 3 + σ m , and υ 4 + σ c present in the dynamics of IHs and parasite populations in the absence of humans (Equations (4)–(7)). A maximum value of one (1) is assumed for each term, the same approach as in Carvalho et al. [36]. Furthermore, we link k to the half-life ( t 1 / 2 ) of the chemical released by using the following equation:
t 1 / 2 = I n 2 / k
where 1 / k is the typical time the chemical persists in the environment and is, therefore, in contact with the HIs, the Schistosoma parasites, and other nontarget species. The results of substituting σ s , σ I , and k and substituting Equations (8) and (10) into Equation (9) are shown in Table 1.
Thus, chemical control guarantees the elimination of all three water-dwelling schistosomiasis agents [16]. Furthermore, based on Matthies and Beulke [27] and the European Commission [37], the chemical is nonpersistent (NP) in the environment when t 1 / 2 = 3.1 40 , persistent (P) when t 1 / 2 = 41 60 , and very persistent (VP) when t 1 / 2 = 61 232 . Based on this information, we discriminate the effects of different chemical performances based on their different half-lives t 1 / 2 , degradation k , and length of chemical exposure 1 / k on the death rate of the targeted species.
It is worth noting that in our model, we assume that chemotherapy treatment reduces disease morbidity among humans while molluscicide increases the mortality of IHs. Public literacy and mechanical control are management strategies for humans and snails, respectively, that directly and extrinsically prevent the likelihood of parasite-host contact for the transmission of schistosomiasis. Furthermore, there is no vertical transmission of the disease to humans nor immigration to the infected individuals. Infected snails are unable to reproduce as a side effect of infection and may die more frequently than susceptible snails.
This paper presents the model analysis first, showing that the model is epidemiologically meaningful, realistic, and of interest in a certain invariant region Ω (Appendix A). We then derive the reproduction number (Section 2.3) and show the steady states (Section 2.4).

2.3. Reproduction Number

In epidemiology, one of the most reliable indicators of infection risk is the basic reproduction number R 0 , which is the average number of new cases of infection caused by an infectious individual in a fully susceptible population [38]. We derived R 0 by rewriting the model system (Equations (1)–(7)), following Chavez et al. [39], into three components: the number of susceptible individuals ( U ), the infected schistosomiasis agents that cannot transmit the infection ( V ), and the infected schistosomiasis agents that transmit the disease ( W ), as follows:
d U d t = f ( U , V , W ) d V d t = g U , V , W d W d t = H U , V , W
where U = ( S h , S v ) , V = ( E h , I h , I V ) , and W = ( M f , C f ) . We let U 0 = ( U * , 0 , 0 ) represent disease-free at U 0 such that, g ~ ( U * , Z ) = ( g ~ 1 ( U * , W ) , g ~ 2 ( U * , W ) ), where g ~ 1 ( U * , W ) = ρ θ h β 1 Λ 1 1 α α h 2 e υ 1 x C f υ 1 C o ω 1 + υ 3 υ 1 + δ 1 + ψ + 1 α α h β 1 υ 1 υ 1 + δ 1 + ψ ω 1 + υ 3 + ρ θ h 1 α α h Λ 2 C f , and g ~ 2 ( U * , W ) = 1 b α v β 2 Λ 2 M f υ 2 υ 2 + δ 2 + σ I M o + ε M f . Suppose, A = D V h ( U * , g ~ ( U * , 0 ) ) , then A is given by
A = υ 4 + σ m ρ θ h N E ω 1 β 1 Λ 1 1 α α h 2 e υ 1 x υ 1 C o ω 1 + υ 3 υ 1 + δ 1 + ψ 1 b α v ω 2 β 2 Λ 2 υ 2 M o υ 2 + δ 2 + σ I υ 5 + σ c
Following the next-generation matrix, an approach used by van den Driessche and Watmough [40], and the concept of reproduction numbers by Diekmann et al. [38], the matrix A can be rewritten as A = F V , with F 0 (i.e., f i j 0 ) and V > 0 , a diagonal matrix, where
F = 0 ρ θ h N E ω 1 β 1 Λ 1 1 α α h 2 e υ 1 x υ 1 C o ω 1 + υ 3 υ 1 + δ 1 + ψ 1 b α v ω 2 β 2 Λ 2 υ 2 M o υ 2 + δ 2 + σ I 0 , V = υ 4 + σ m 0 0 υ 5 + σ p
The basic reproduction number is the spectral radius (dominant eigenvalue) of the matrix F V 1 given by
R 0 = ρ F V 1 = ρ θ h N E ω 1 ω 2 β 1 β 2 Λ 1 Λ 2 1 b α v 1 α α h 2 e υ 1 x υ 1 υ 2 M o C o υ 1 + δ 1 + ψ υ 2 + δ 2 + σ I ω 1 + υ 3 υ 4 + σ m υ 5 + σ c

2.4. Steady State

We show that the model system, Equations (1)–(7) have a disease-free state ( E 0 ) and an endemic steady state E 1 , where
E 0 = S h * , I h * , E h * , M f * S v * , I v * , C f * = Λ 1 e υ 1 x υ 1 , 0 , 0 , 0 , Λ 2 υ 2 , 0 , 0
It always exists in R + 0 7 provided R 0 < 1 .
E 1 = S h 1 , I h 1 , E h 1 , M f 1 S v 1 , I v 1 , C f 1 is expressed in terms of I v 1 and E h 1 by
S h 1 ( I v 1 ) = υ 1 + δ 1 + ψ Λ 1 υ 5 + σ c C o + ω 2 ε I v 1 υ 1 + δ 1 + ψ υ 5 + σ c C o + ω 2 υ 1 + δ 1 + ψ 1 α α h β 1 + ε ψ 1 α α h β 1 I v 1 , I h 1 ( I v 1 ) = β 1 ω 2 1 α α h υ 1 + δ 1 + ψ Λ 1 υ 5 + σ c C o + ω 2 ε I v 1 I v 1 υ 1 + δ 1 + ψ υ 5 + σ c C o + ω 2 I v 1 υ 1 + δ 1 + ψ υ 5 + σ c C o + ω 2 υ 1 + δ 1 + ψ 1 α α h β 1 + ε ψ 1 α α h β 1 I v 1 E h 1 I v 1 = β 1 ω 2 ρ θ h K 1 α α h υ 1 + δ 1 + ψ Λ 1 υ 5 + σ c C o + ω 2 ε I v 1 I v 1 υ 1 + δ 1 + ψ υ 5 + σ c ω 1 + υ 3 K C o + ( ω 1 + υ 3 K ω 2 + ρ θ h 1 α α h ω 2 ) I v 1 υ 1 + δ 1 + ψ υ 5 + σ c C o + ω 2 υ 1 + δ 1 + ψ 1 α α h β 1 + ε ψ 1 α α h β 1 I v 1 M f 1 I v 1 = β 1 ω 2 ρ θ h K N E ω 1 1 α α h υ 1 + δ 1 + ψ Λ 1 υ 5 + σ c C o + ω 2 ε I v 1 I v 1 υ 1 + δ 1 + ψ υ 4 + σ m υ 5 + σ c ω 1 + υ 3 K C o + ( ω 1 + υ 3 K ω 2 + ρ θ h 1 α α h ω 2 ) I v 1 υ 1 + δ 1 + ψ υ 5 + σ c C o + ω 2 υ 1 + δ 1 + ψ 1 α α h β 1 + ε ψ 1 α α h β 1 I v 1 S v 1 E h 1 I v 1 = Λ 2 M o + Λ 2 ε E h 1 I v 1 υ 2 + σ s M o + υ 2 + σ s ε + 1 b α v β 2 E h 1 I v 1 C f 1 I v 1 = ω 2 I v 1 ( υ 5 + σ c )
when I v ˙ t = 0, we solved I v 1 using Equation (6) of the model Equations (1)–(7), and this leads to the substitution of S v 1 E h 1 I v 1 and M f 1 I v 1 , which yields Equation (11) below.
( a 4 I v 1 4 + a 3 I v 1 3 + a 2 I v 1 2 + a 1 I v 1 + a 0 ) I v 1 = 0
where
a 0 = d 1 d 2 d 3 d 4 d 5 d 6 C 0 M 0 K Λ 1 K k 1 k 2 k 3 k 4 ω 2 d 1 d 6 M 0 Λ 1 · K d 1 d 2 d 3 ω 2 d 6 ε C 0 M 0 Λ 1 Λ 2 a 1 = d 5 ( d 1 d 2 d 3 d 4 d 5 d 6 C 0 M 0 K Λ 1 + M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε ω 2 M 0 C 0 d 6 d 4 d 3 d 1 K Λ 1 + ε Λ 2 K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ) K k 1 k 2 k 3 k 4 ω 2 d 1 d 6 M 0 Λ 1 K k 1 k 2 k 3 k 4 ω 2 d 1 ε C 0 Λ 1 · K d 1 d 2 d 3 ω 2 d 6 ε C 0 M 0 Λ 1 Λ 2 a 2 = d 5 ( d 1 d 2 d 3 d 4 d 5 d 6 C 0 M 0 K Λ 1 + M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε ω 2 M 0 C 0 d 6 d 4 d 3 d 1 K Λ 1   + ε Λ 2 K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 + M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ) ε M 0 C 0 d 1 d 6 Λ 1 Λ 2 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + M 0 C 0 d 1 d 2 2 d 3 2 K Λ 2 + ε Λ 2 K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ( K k 1 k 2 k 3 k 4 ω 2 d 1 d 6 M 0 Λ 1 + K k 1 k 2 k 3 k 4 ω 2 d 1 ε C 0 Λ 1 ) a 3 = d 5 ( M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε ω 2 M 0 C 0 d 6 d 4 d 3 d 1 K Λ 1 + ε Λ 2 K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 + M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ) K k 1 k 2 k 3 k 4 ω 2 d 1 ε C 0 Λ 1 · ( M 0 d 1 d 3 ε ω 2 Λ 1 Λ 2 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ) a 4 = d 5 ( M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε K k 1 k 2 k 3 ω 2 C 0 d 1 d 6 Λ 1 ) ( M 0 C 0 d 3 d 4 d 6 Λ 1 K d 1 d 2 ω 2 + k 1 k 2 ω 2 + ε ω 2 M 0 C 0 d 6 d 4 d 3 d 1 K Λ 1 ) k 1 = 1 α α h β 1 , k 2 = ρ θ h 1 α α h , k 3 = N E ω 1 , k 4 = 1 b α v β 2 , d 1 = υ 1 + δ 1 + ψ ,
d 2 = ω 1 + υ 3 , d 3 = υ 4 + σ m , d 4 = υ 2 + σ s , d 5 = υ 2 + δ 2 + σ I , d 6 = ( υ 5 + σ c ) . Equation (11) has one of the solutions I v 1 = 0 , which corresponds to a disease-free equilibrium. Other solutions can be determined according to the sign rule of Descartes [41], which states that depending on the change in the sign of the coefficients, a 0 , a 1 , a 2 , and a 3 , 3, 2, or 1 positive solutions exist; thus, the proof for existing conditions for E 1 for the model Equations (1)–(7) is complete.
In addition, the stability of the model system depends on the parameters used to build the model and the changes in these parameters affect R 0 , with numerous scientific and biological implications for the transmission and control of the disease. For instance, if we use β 1 = β * as the bifurcation parameter for model Equations (1)–(7), then R 0 = 1 and
β * = υ 1 υ 2 M o C o υ 1 + δ 1 + ψ υ 2 + δ 2 + σ I ω 1 + υ 3 υ 4 + σ m υ 5 + σ c ρ θ h N E ω 1 ω 2 β 2 Λ 1 Λ 2 1 b α v 1 α α h 2 e υ 1 x
is the crucial bifurcation value for cercaria-to-human transmission. When R 0 < 1 and β 1 < β * , for instance, only the asymptotically stable disease-free equilibrium point ( E 0 ) can exist, and R 0 provides conditions for disease extinction. The disease-free equilibrium is unstable, whereas the endemic equilibrium ( E 1 ) is asymptotically stable, and both equilibrium points exist when R 0 > 1 and β 1 > β * [40], and in this case, R 0 provides the conditions for the disease to persist. Thus, the asymptotic dynamic behavior of the infectious disease, which determines whether it will vanish or persist in the future, can be inferred from the steady states. Although we do not provide a stability analysis of the equilibrium points, we acknowledge that control strategies have a role to play in reducing R 0 and possibly eliminating schistosomiasis.

2.5. Parameter Data

For our study, we use data collected in the literature related to Bulinus and Biomphalaria (Table 2). In the absence of published data, parameter values are assumed or estimated based on expert knowledge using what is commonly known about vector and disease dynamics, as follows. The life expectancy of an adult Schistosoma worm within a human host is 3–10.5 years [42,43,44]; thus, δ 1 vary from 1 / ( 10.5 × 365 ) 0.000268 to 1 / ( 3 × 365 ) 0.000913 . In addition, Dabo et al. [45] establish schistosomiasis infection in preschool children aged 1–4 years through both passive and active exposure to infected water bodies. On average, a child`s initial infection age is 2 years [43]; thus, x = 2 × 365 = 730 days. Gryseels et al. [46] reported 2–10 weeks as the period a schistosomiasis patient takes to recover after chemotherapy treatment. Thus, we assume that the recovery period for a schistosomiasis patient varies between 14 and 70 days, which sets ψ between ψ = 1 / 14 0.0714 and ψ = 1 / 70 0.0143 individuals per day. The Schistosoma parasite egg can stay for 7 days, with a 0.14286 per egg death rate, and miracidia survive 12 days in freshwater ( υ 3 = 1 / 12 0.0833 ) [46,47]. Cercaria can only spend 10–40 h in freshwater [48], which sets υ 4 between υ 4 = 1 / ( 10 × 24 ) 0.00417 and υ 4 = 1 / ( 40 × 24 ) 0.00104 . Table 2 shows the numerical values (and potential ranges) of parameters set for the model system Equations (1)–(7) from the literature and estimations.
We performed a numerical simulation using the main R package ODE solver Version 1.10-4 for solving ordinary differential equations [53] in the R statistical environment version 4.0.3 [54]. The baseline values in Table 2 were used for all simulations and predictions. We conducted a sensitivity analysis of R 0 with respect to the baseline value using the Partial Rank Correlation Coefficients (PRCC) test to determine how robust our predictions are to changes in parameter values. In this way, we can determine quantitatively which key parameters can be targeted by control measures to reduce disease the most. For instance, R 0 increases when parameters with positive PRCC values are increased, increasing the likelihood of infection.

3. Results

The results show that the key parameters (bold in Table 3) that have the greatest influence on the transmission of the disease include the recruitment rate of human individuals Λ 1 , a portion of stool/urine per infected person ρ, number of miracidia per parasite egg N E , cercaria-to-human transmission per contact β1, rate of miracidia-to-snail transmission per contact β2, and rate of cercaria emergence from infected snails ω 2 .
The model system Equations (1)–(7) is solved numerically using the initial conditions S h 0 = 100,000 , S h 0 = 100,000 , I h 0 = 1 , E h 0 = 0 , M f 0 = 0 , S v 0 = 100,000 , I v 0 = 1 , and C f 0 = 0 , and the effects of different control strategies on the targetted human and snail populations are shown graphically in Figure 2. When a chemotherapy treatment of the infected humans is integrated with public literacy (two-tiered approach), the strategy is more successful in managing and reducing disease morbidity in the human population than when each strategy is carried out separately, i.e., the number of infected humans who recover, reduce, and join the susceptible humans, who grow in number (Figure 2a,b). Chemical/molluscicide control significantly increases IHs mortality, and when paired with a mechanical technique, this effect is essentially identical (Figure 2c,d).
The impact of individual control measures and various combinations of all control measures on disease transmission is determined by the value of R 0 , as shown in Table 4. The more successful a certain level of control is, the smaller the value of R 0 . We show that when a single control strategy is used, the chemical/molluscicides strategy ( σ I , M , C ) has a greater impact on reducing the value of R 0 followed by the chemotherapy treatment strategy ( ψ ), public literacy ( α α h ), and finally mechanical control ( b α v ). However, a single intervention cannot eradicate the disease, R 0 1 : see Table 4. The results of integrating the various strategies show that a combination of two-tiered control measures, treatment strategies ( ψ ), and molluscicides ( σ I , M , C ) have the greatest impact compared to any other pair of strategies (Table 4). Moreover, a combination of a three-tiered approach: treatment ( ψ ), molluscicides ( σ I , M , C ), and mechanical ( b α v ), as well as another three-tiered approach: public literacy ( α α h ) , molluscicides σ I , M , C , and mechanical ( b α v ) have the potential to lower R 0 below 1, but only if the effectiveness of each separate control is greater than 90% (Table 4). However, schistosomiasis can be eradicated, R 0 < 1 , if a four-tiered approach (treatment, public literacy, molluscicides, and mechanical) is applied, with each method being more than 70% effective when combined (Table 4).
When considering mortality rate functions (Table 1), we show that temperature, the length of chemical exposure, and the chemical’s half-life affect the molluscicide performance on the death rates of targeted species. In general, chemical-induced mortality rates for targeted species (host snails and free-living Schistosoma cercaria and miracidia) decline as the temperature increases. However, the decline is more rapid when the half-life of molluscicides decreases, i.e., shorter half-lived molluscicides (e.g., non-persistent and persistent molluscicides) lose their toxicity sooner, leading to lower mortality rates. For instance, Figure 3a shows how temperature control impacts the half-life of molluscicides and the consequent mortality rate of susceptible snails. In the same context, chemical-induced mortality rates decrease with an increase in the duration of molluscicide exposure. Molluscicides with a longer half-life (very persistent molluscicides) remain in the environment longer, resulting in a higher mortality rate compared to molluscicides with a shorter half-life (Figure 3b). Thus, increasing the half-life of molluscicides increases the chemically induced mortality rates, especially for cercariae σ c , and infected snails σ I (Figure 3c). As a result, increasing the temperature and days of chemical exposure decreases the mortality rate of snails (Figure 3d) associated with higher risks of disease transmission, thus increasing the R 0 value (Figure 3e).

4. Discussion

We formulate a mathematical model to assess the effectiveness of methods that target humans and those targeting IHs in the control of schistosomiasis. The results of the PRCC test suggest that control strategies focusing on the key parameters N E , β1, β2, and ω 2 in the transmission dynamics may be more effective in reducing R 0 and may result in disease-free conditions of R 0 < 1 , indicating disease control. Our results show that the chemotherapy treatment strategy (reducing N E ) and application of molluscicides (reducing β1, β2, and ω 2 ) are the most successful methods for reducing R 0 and disease morbidity among the human population and increasing snail mortality, respectively, when applied independently and compared to mechanical control (reducing ω 2 ) and public literacy (reducing β1 and β2) with limited impact. But each strategy cannot eradicate schistosomiasis individually, which complements the finding by Mangal et al. [52], King and Bertsch [9], and Zacharia et al. [10]. Thus, a single approach may reduce morbidity [2,55] but often without reducing the local and environmental transmission of parasites [6], which contributes to the development of the disease. In addition, chemotherapy has shown ineffectiveness against juvenile worms within the body. This is the reason why worms inside infected and treated children can continue to release eggs, observed in their feces three weeks after the drug has been administered [56,57,58]. This thus perpetuates local transmission, and residents who stay in endemic areas run the risk of re-infection [6,59]. The fragility of drug control alone is demonstrated by this.
On the other hand, mechanical control practices like the physical picking of snails are considered old-fashioned but prove to be a good complementary strategy to drugs, molluscicides, and public literacy [17]. However, its use may be limited, especially in large water bodies, such as lakes, rivers, and streams, or where a larger number of aquatic inhabitants must be controlled. In addition, public literacy advocates for significant behavior change in water use and contact. However, this ultimately requires the provision of safe and hygienic alternative water supplies, latrines, and washing facilities. When compared to using individual control strategies, integrated strategies are far more effective. Our study further shows that when a four-tiered control strategy is used, the results are much better than using one-, two- or three-tiered controls. Although a three-tiered strategy would be effective as well, it is practically difficult to achieve 90% effectiveness from each method. These findings support, to a certain extent, the findings of Mangal et al. [52], Thétiot-Laurent et al. [56], Lo et al. [58], and Zheng et al. [13], which suggest that the integrated intervention of chemotherapy and snail management has a greater impact on schistosomiasis control than a single method like chemotherapy treatment alone.
However, in such integrated efforts to control schistosomiasis, the time and season of application of molluscicides are crucial. This is relevant in providing a realistic assessment of the spread of schistosomiasis under the current and projected temperature increases in areas like sub-Saharan Africa, where the disease is most prevalent. Our results corroborate with the reports by Feng et al. [60], Fishel [25], Ziska [26], and Zheng et al. [13] that the temperature of the environment, chemical degradation, and half-life of the chemical are important factors in the chemical control of snails. Regardless of the chemical type used, non-persistent (NP), persistent (P), or very persistent (VP), the results show that the chemical strategy may be more effective in areas with low temperatures (15–25 °C). In the first two weeks (14 days) of application, the chemically induced mortality rate of snails is significantly higher, especially for cercariae and infected snails, and consequently, R 0 is less than at unity, indicating disease control. This corroborates the finding by Montanari et al. [33] that, at a standard temperature of 25 °C, molluscicides proved to be effective against adult snails; the snails reduced their food intake and stopped eating, whereas the same effects were significantly stronger at 18 °C. No effects on survival or feeding rate have been recorded at incubations of 30 °C. In regions with higher temperatures, an even longer treatment period using NP chemicals is associated with low rates of induced mortality and is, therefore, less effective in reducing snails, free-living parasite populations, and R 0 . The results also show that the molluscicides NP, P, and VP with increasing half-lives reduce species mortality by more than 50% only up to 2 weeks, 3 weeks, and 6 weeks, respectively, suggesting that the chemical should be reapplied regularly thereafter. Thus, the half-life of the molluscicide is directly proportional to the active length of the molluscicides in the water. This result supports the finding of Seligman et al. [61], who observed that the half-lives/chemical degradation of different chemicals ranges from 4 to 19 days and up to several months under specific conditions. For efficient snail management, integration efforts, and schistosomiasis control, different areas with varying weather patterns and seasonal variations may require different molluscicide-type and variable time intervals of reapplication.
Therefore, it is critical to quantify the use of chemical control agents in response to climatic conditions to determine alternative changes in the target intermediate host population and the best season (temperature) for control. This result complements the efforts to investigate molluscicidal action, the structure-activity relationship, and the potential mechanisms of several molluscicides [13,14,62] toward developing suitable chemicals.
Our findings, thus, demonstrate that effective schistosomiasis control programs necessitate a co-ordinated effort and comprehensive integration of preventative treatment, public literacy, and mechanical and molluscicide controls. Despite the promising results, we believe that controlling schistosomiasis with fewer than three different combinations of control strategies will be challenging and may not be worth the effort to eradicate schistosomiasis compared to when the four strategies are integrated and carried out simultaneously. However, achieving the schistosomiasis eradication goal may require a large investment of resources (high cost, which severely limits the use of the molluscicide on a big scale), time, and commitment, and also has associated ecological risks. For instance, chemical management has several negative effects, including the destruction of aquatic biodiversity due to toxicity [63] and the emergence of snail resistance to the chemicals are very likely. If chemicals are released into the environment, local people will not be able to use the water [13,64]. Several molluscicides have a half-life that ranges from 3.17 to 223 days [65], while others have a half-life of up to 2 years. This means that the presence of molluscicides in water would continue to pose a health risk to people for more than 30 years before it decays completely [61].
Furthermore, the development and use of molluscicide need to be blended with an increase in health education and an assessment of their levels of toxicity when ingested by humans and domestic animals through agricultural products and their effect on non-target aquatic fauna and flora [13]. However, our results suggest that chemicals with shorter half-lives, such as an extract from C. viminalis fruits, or natural insecticides, such as Nereistoxin, are effective in killing snails at low concentrations in a shorter amount of time [13,66]. They could provide an alternative to long-term exposures to common molluscicides with less effects on biodiversity. Furthermore, specific molluscicides, such as nicotinanilide molluscicides and silver (Ag) nanopowder (Ag-NPs) molluscicides [13], have been shown to exhibit strong molluscicidal effects with less toxicity to humans, animals, fish, and plants, can be applied in small waterbodies to manage snail populations despite their high cost and strong human skin irritant effects. The effect of change in temperature on the environment on molluscicide performance greatly affects the mortality rate of the targeted organism (IHs, cercariae, and miracidia) and general disease transmission. Therefore, better molluscicide performance should be achieved by developing harmless, low-toxicity, and temperature-dependent molluscicides, with low environmental impact to control the targeted species holistically.
Our study has some limitations given that the parameters that best reflect the biological representation and real-life situations of schistosomiasis disease and transmission had to be used based on the published record. Since neither the Biomphalaria-S. mansoni nor the Bulinus-S. haematobium infection system had a complete list of genus-specific model parameters, we used published parameter values for both infection systems. Due to the absence of specifically published data relating to the two systems, two parameter values were estimated, and two others were assumed based on expert knowledge. Furthermore, the literature data are not always consistent and are collected under different conditions. Nonetheless, the baseline values of such parameters provide a good first approximation of the parameter values, and the model can be used generally to develop a better understanding of the effects of interventions. This model can be improved to incorporate more features, such as the effect of temperature and molluscicides on the development of eggs into miracidia, snail recruitment, and biodiversity of nontarget species.

5. Conclusions

We conclude that the proposed model provides valuable mechanistic insights into how the impact of schistosomiasis and future control programs will depend on the integration of chemotherapy, public literacy, and mechanical and molluscicide types in different endemic regions. In addition, molluscicide performances in different regions with changing temperatures will vary depending on their type and half-life, degradation rate, and reapplication frequency.

Author Contributions

Conceptualization, Z.T.; Methodology, Z.T. and C.K.; Formal analysis, Z.T.; Funding acquisition, L.B. and C.A.; Supervision, C.A. and L.B.; Visualization, Z.T.; Writing—original draft, Z.T.; Writing—review & editing, Z.T., C.K., L.B. and C.A. All authors have read and agreed to the published version of the manuscript.

Funding

The German Academic Exchange Service (DAAD) supported Z.T with a Ph.D. scholarship (Grant No. 57507871). The Alexander von Humboldt Foundation supported C.K as a Georg Forster Research Fellow with a scholarship for his research visit and stay at Justus Liebig University Giessen. (Grant No. ZMB1217528GF-P). The funders had no role in study design, data collection or analysis, the decision to publish, or the preparation of the manuscript.

Data Availability Statement

All data generated or analyzed during this study are included in this published article.

Acknowledgments

It is with great gratitude that the authors acknowledge the German Academic Exchange Service (DAAD) for providing Z.T. with a Ph.D. scholarship, as well as Alexander von Humboldt Foundation for providing C.K. with a scholarship.

Conflicts of Interest

The authors declare that they do not have any competing interests.

Appendix A

Invariant region
In model system Equations (1)–(7), we prove that non-negative and bounded solutions exist in Ω, as follows:
Let L t = S h t , I h t , E h t , M f t , S v t , I v t , C f t T R + 7 be the solution set of the model equations Equations (1)–(7) from which L 0 0 are the non-negative initial conditions.
Lemma A1. 
L t is non-negative in  Ω R + 7  with  L 0 0  for  t 0 .
Proof. 
It is evident from model equation Equation (1) that
S h ˙ t 1 α α h β 1 C f C o + ε C f + υ 1 S h , where S h t S h 0 e 1 α α h β 1 C f C o + ε C f + υ 1 t , t 0 .
Thus, S h t 0 is non-negative and, in the same way, it can be shown that I h t 0 , E h t 0 , M f t 0 , S v t 0 , I v t 0 , and C f t 0 . Therefore, L t in the region Ω R + 7 is positively invariant. □
Lemma A2. 
L t  remains bounded in  Ω R + 7 .
Proof. 
Let the total human population S h t + I h t = N h t , and the total snail population S v t + I v t = N v t , such that N h ˙ t Λ 1 e υ 1 x υ 1 N h t and N v ˙ t Λ 2 υ 2 N v t and L t = N h t , E h t , M f t , N v t , C f t T . Thus, L t can be expressed analytically as follows:
L t = N h t Λ 1 e υ 1 x υ 1 Λ 1 e υ 1 x υ 1 N h 0 e υ 1 t E h t K K E h 0 e ρ θ h 1 α α h I h t M f t N E ω 1 E h υ 4 + σ m N E ω 1 E h υ 4 + σ m M f 0 e υ 4 + σ m t N v t Λ 2 υ 2 Λ 2 υ 2 N v 0 e υ 2 t C f t ω 2 I v υ 5 + σ c ω 2 I v υ 5 + σ c C f 0 e υ 5 + σ c t
Assuming N h 0 Λ 1 e υ 1 x υ 1 , and N v 0 Λ 2 υ 2 , we conclude that lim t S u p N h t = Λ 1 e υ 1 x υ 1 , and lim t S u p N v t = Λ 2 υ 2 , and therefore, both human ( S h t , I h t ) and snail ( S v t , I v t ) populations are biologically feasible. Similarly lim t S u p E h t , M f t , C f t T = K , N E ω 1 E h υ 4 + σ m , ω 2 I v υ 5 + σ c T , shows that the egg, miracidia, and cercaria populations are biologically feasible. Thus, L t with initial conditions remains bounded for all t 0 . □
As a result of Lemmas A1 and A2, our model is epidemiologically and mathematically well-posed.

References

  1. Adenowo, A.F.; Oyinloye, B.E.; Ogunyinka, B.I.; Kappo, A.P. Impact of human schistosomiasis in sub-Saharan Africa. Braz. J. Infect. Dis. 2015, 19, 196–205. [Google Scholar] [CrossRef] [Green Version]
  2. WHO, World Health Organization. Schistosomiasis. Available online: https://www.who.int/news-room/fact-sheets/detail/schistosomiasis (accessed on 3 February 2023).
  3. Wang, L.; Utzinger, J.; Zhou, X.-N. Schistosomiasis control: Experiences and lessons from China. Lancet 2008, 372, 1793–1795. [Google Scholar] [CrossRef]
  4. Rosenberg, M.; Utzinger, J.; Addiss, D.G. Preventive chemotherapy versus innovative and intensified disease management in neglected tropical diseases: A distinction whose shelf life has expired. PLoS Negl. Trop. Dis. 2016, 10, e0004521. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. WHO, World Health Organization. Elimination of Schistosomiasis. Sixty-Fifth World Health Assembly, WHA65. 21, Agenda item 13.11. 26 May 2012. Available online: http://apps.who.int/gb/ebwha/pdf_files/WHA65/A65_R21-en.pdf (accessed on 23 February 2023).
  6. French, M.D.; Churcher, T.S.; Gambhir, M.; Fenwick, A.; Webster, J.P.; Kabatereine, N.B.; Basáñez, M.-G. Observed reductions in Schistosoma mansoni transmission from large-scale administration of praziquantel in Uganda: A mathematical modelling study. PLoS Negl. Trop. Dis. 2010, 4, e897. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Stothard, J.R.; Chitsulo, L.; Kristensen, T.K.; Utzinger, J. Control of schistosomiasis in sub-Saharan Africa: Progress made, new opportunities and remaining challenges. Parasitology 2009, 136, 1665–1675. [Google Scholar] [CrossRef] [Green Version]
  8. Hotez, P.J.; Molyneux, D.H.; Fenwick, A.; Kumaresan, J.; Sachs, S.E.; Sachs, J.D.; Savioli, L. Control of neglected tropical diseases. N. Engl. J. Med. 2007, 357, 1018–1027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. King, C.H.; Bertsch, D. Historical perspective: Snail control to prevent schistosomiasis. PLoS Negl. Trop. Dis. 2015, 9, e0003657. [Google Scholar] [CrossRef] [Green Version]
  10. Zacharia, A.; Mushi, V.; Makene, T. A systematic review and meta-analysis on the rate of human schistosomiasis reinfection. PLoS ONE 2020, 15, e0243224. [Google Scholar] [CrossRef] [PubMed]
  11. Madsen, H. Biological methods for the control of freshwater snails. Parasitol. Today 1990, 6, 237–241. [Google Scholar] [CrossRef]
  12. Sokolow, S.H.; Huttinger, E.; Jouanard, N.; Hsieh, M.H.; Lafferty, K.D.; Kuris, A.M.; Riveau, G.; Senghor, S.; Thiam, C.; N’Diaye, A.; et al. Reduced transmission of human schistosomiasis after restoration of a native river prawn that preys on the snail intermediate host. Proc. Natl. Acad. Sci. USA 2015, 112, 9650–9655. [Google Scholar] [CrossRef] [Green Version]
  13. Zheng, L.; Deng, L.; Zhong, Y.; Wang, Y.; Guo, W.; Fan, X. Molluscicides against the snail-intermediate host of Schistosoma: A review. Parasitol. Res. 2021, 120, 3355–3393. [Google Scholar] [CrossRef] [PubMed]
  14. King, C.; Sutherland, L.J.; Bertsch, D. Systematic review and meta-analysis of the impact of chemical-based mollusciciding for control of Schistosoma mansoni and S. haematobium transmission. PLoS Negl. Trop. Dis. 2015, 9, e0004290. [Google Scholar] [CrossRef]
  15. Lo, N.C.; Gurarie, D.; Yoon, N.; Coulibaly, J.T.; Bendavid, E.; Andrews, J.R.; King, C.H. Impact and cost-effectiveness of snail control to achieve disease control targets for schistosomiasis. Proc. Natl. Acad. Sci. USA 2018, 115, E584–E591. [Google Scholar] [CrossRef] [Green Version]
  16. WHO, World Health Organization. Field Use of Molluscicides in Schistosomiasis Control Programmes: An Operational Manual for Programme Managers; Licence: CC BY-NC-SA 3.0 IGO; World Health Organization: Geneva, Switzerland, 2017. [Google Scholar]
  17. Sokolow, S.H.; Wood, C.L.; Jones, I.J.; Lafferty, K.D.; Kuris, A.M.; Hsieh, M.H.; De Leo, G.A. To reduce the global burden of human schistosomiasis, use ‘old fashioned’ snail control. Trends Parasitol. 2017, 34, 23–40. [Google Scholar] [CrossRef]
  18. Chiyaka, E.T.; Garira, W. Mathematical analysis of the transmission dynamics of schistosomiasis in the human-snail hosts. J. Biol. Syst. 2009, 17, 397–423. [Google Scholar] [CrossRef]
  19. Gao, S.; Liu, Y.; Luo, Y.; Xie, D. Control problems of a mathematical model for schistosomiasis transmission dynamics. Nonlinear Dyn. 2010, 63, 503–512. [Google Scholar] [CrossRef]
  20. Abokwara, A.; Madubueze, C.E. The Role of Non-pharmacological Interventions on the Dynamics of Schistosomiasis. J. Math. Fundam. Sci. 2021, 53, 243–260. [Google Scholar] [CrossRef]
  21. Nur, W.; Trisilowati; Suryanto, A.; Kusumawinahyu, W.M. Mathematical model of schistosomiasis with health education and molluscicide intervention. J. Phys. Conf. Ser. 2021, 1821, 012033. [Google Scholar] [CrossRef]
  22. Stensgaard, A.-S.; Booth, M.; Nikulin, G.; McCreesh, N. Combining process-based and correlative models improves predictions of climate change effects on Schistosoma mansoni transmission in eastern Africa. Geospat. Health 2016, 11, 406. [Google Scholar] [CrossRef]
  23. De Leo, G.A.; Stensgaard, A.-S.; Sokolow, S.H.; N’goran, E.K.; Chamberlin, A.J.; Yang, G.-J.; Utzinger, J. Schistosomiasis and climate change. BMJ 2020, 371, m4324. [Google Scholar] [CrossRef]
  24. Tabo, Z.; Neubauer, T.A.; Tumwebaze, I.; Stelbrink, B.; Breuer, L.; Hammoud, C.; Albrecht, C. Factors controlling the distribution of intermediate host snails of Schistosoma in crater lakes in Uganda: A machine learning approach. Front. Environ. Sci. 2022, 10, 871735. [Google Scholar] [CrossRef]
  25. Fishel, F.M. Storage Limitation Statements: Temperature–Herbicides: PI123/PI160, 4/2013. EDIS 2013, 4, 123–160. [Google Scholar]
  26. Ziska, L.H. Increasing minimum daily temperatures are associated with enhanced pesticide use in cultivated soybean along a latitudinal gradient in the mid-western United States. PLoS ONE 2014, 9, e98516. [Google Scholar] [CrossRef] [PubMed]
  27. Matthies, M.; Beulke, S. Considerations of temperature in the context of the persistence classification in the EU. Environ. Sci. Eur. 2017, 29, 15. [Google Scholar] [CrossRef] [Green Version]
  28. European Food Safety Authority (EFSA). Opinion on a request from EFSA related to the default Q10 value used to describe the temperature effect on transformation rates of pesticides in soil-Scientific Opinion of the Panel on Plant Protection Products and their Residues (PPR Panel). EFSA J. 2008, 6, 622. [Google Scholar] [CrossRef]
  29. ECHA. Guidance on Information Requirements and Chemical Safety Assessment; Chapter R.11: PBT/vPvB Assessment; Version 2; Technical Report ECHA-14-G-07-EN; European Chemicals Agency: Helsinki, Finland, 2014. [Google Scholar]
  30. ECHA. Guidance for Information Requirements and Chemical Safety Assessment; Chapter R.7b: Endpoint Specific Guidance; Draft Version 4.0 (Public); European Chemicals Agency: Helsinki, Finland, 2016. [Google Scholar]
  31. Ronoh, M.; Chirove, F.; Pedro, S.A.; Tchamga, M.S.S.; Madubueze, C.E.; Madubueze, S.C.; Addawe, J.; Mwamtobe, P.M.; Mbra, K.R. Modelling the spread of schistosomiasis in humans with environmental transmission. Appl. Math. Model. 2021, 95, 159–175. [Google Scholar] [CrossRef]
  32. Grimes, J.E.T.; Croll, D.; Harrison, W.E.; Utzinger, J.; Freeman, M.C.; Templeton, M.R. The roles of water, sanitation and hygiene in reducing schistosomiasis: A review. Parasites Vectors 2015, 8, 156. [Google Scholar] [CrossRef] [Green Version]
  33. Montanari, A.L.; Accorsi, A.; Nasi, M.; Malagoli, D. Effects of a nematode-based molluscicide on survival and antimicrobial peptide expression in Pomacea canaliculata. Invertebr. Surviv. J. 2019, 16, 37. [Google Scholar]
  34. Laidler, K.J.; Chen, H.; Ling, M.; Hencz, L.; Ling, H.Y.; Li, G.; Lin, Z.; Liu, G.; Zhang, S.; Abyazisani, M.; et al. The development of the Arrhenius equation. J. Chem. Educ. 1984, 61, 494. [Google Scholar] [CrossRef]
  35. Lewis, K.; Tzilivakis, J.; Green, A.; Warner, D. Pesticide Properties DataBase (PPDB). 2006. Available online: http://hdl.handle.net/2299/15375 (accessed on 14 February 2023).
  36. Carvalho, S.A.; da Silva, S.O.; Charret, I.D.C. Mathematical modeling of dengue epidemic: Control methods and vaccination strategies. Theory Biosci. 2019, 138, 223–239. [Google Scholar] [CrossRef] [Green Version]
  37. European Commission EC. Technical Guidance Document on Risk Assessment in Support of Commission Directive 93/67/EEC on Risk Assessment for New Notified Substances Part II. In Commission Regulation (EC) No 1488/94 on Risk Assessment for Existing Substances and of Directive 98/8/EC of the European Parliament and of the Council Concerning the Placing of Biocidal Products on the Market; Technical Report EUR 20418 EN/2; European Commission EC: Luxembourg, 2003. [Google Scholar]
  38. Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A.J. On the definition and the computation of the basic reproduction ratio R 0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 1990, 28, 365–382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Chavez, C.C.; Feng, Z.; Huang, W. On the computation of R0 and its role on global stability. In Mathematical Approaches for Emerging and Re-Emerging Infection Diseases: An Introduction; Springer: Berlin/Heidelberg, Germany, 2002; Volume IMA 125, pp. 229–250. [Google Scholar]
  40. Driessche, P.V.D.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2002, 180, 29–48. [Google Scholar] [CrossRef] [PubMed]
  41. Anderson, B.; Jackson, J.; Sitharam, M. Descartes’ rule of signs revisited. Am. Math. Mon. 1998, 105, 447. [Google Scholar] [CrossRef]
  42. Fulford, A.J.C.; Butterworth, A.E.; Ouma, J.H.; Sturrock, R.F. A statistical approach to schistosome population dynamics and estimation of the life-span of Schistosoma mansoni in man. Parasitology 1995, 110, 307–316. [Google Scholar] [CrossRef]
  43. Colley, D.G.; Bustinduy, A.L.; Secor, W.E.; King, C.H. Human schistosomiasis. Lancet 2014, 383, 2253–2264. [Google Scholar] [CrossRef]
  44. Anderson, R.; Turner, H.; Farrell, S.; Truscott, J. Studies of the transmission dynamics, mathematical model development and the control of schistosome parasites by mass drug administration in human communities. Adv. Parasitol. 2016, 94, 199–246. [Google Scholar] [CrossRef]
  45. Dabo, A.; Badawi, H.M.; Bary, B.; Doumbo, O.K. Urinary schistosomiasis among preschool-aged children in Sahelian rural communities in Mali. Parasites Vectors 2011, 4, 21. [Google Scholar] [CrossRef] [Green Version]
  46. Gryseels, B.; Polman, K.; Clerinx, J.; Kestens, L. Human schistosomiasis. Lancet 2006, 368, 1106–1118. [Google Scholar] [CrossRef]
  47. Michaels, R.M.; Prata, A. Evolution and characteristics of Schistosoma mansoni eggs laid in vitro. J. Parasitol. 1968, 54, 921. [Google Scholar] [CrossRef]
  48. Braun, L.; Grimes, J.E.T.; Templeton, M.R. The effectiveness of water treatment processes against schistosome cercariae: A systematic review. PLoS Negl. Trop. Dis. 2018, 12, e0006364. [Google Scholar] [CrossRef]
  49. Kanyi, E.; Afolabi, A.S.; Onyango, N.O. Mathematical modelling and analysis of transmission dynamics and control of schistosomiasis. J. Appl. Math. 2021, 2021, 1–20. [Google Scholar] [CrossRef]
  50. Feng, Z.; Eppert, A.; Milner, F.; Minchella, D. Estimation of parameters governing the transmission dynamics of schistosomes. Appl. Math. Lett. 2004, 17, 1105–1112. [Google Scholar] [CrossRef] [Green Version]
  51. Liang, S.; Spear, R.C.; Seto, E.; Hubbard, A.; Qiu, D. A multi-group model of Schistosoma japonicum transmission dynamics and control: Model calibration and control prediction. Trop. Med. Int. Health 2005, 10, 263–278. [Google Scholar] [CrossRef]
  52. Mangal, T.D.; 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. [Google Scholar] [CrossRef] [Green Version]
  53. Soetaert, K.; Petzoldt, T.; Setzer, R.W. Solving differential equations in R: Package deSolve. J. Stat. Softw. 2010, 33, 1–25. [Google Scholar] [CrossRef]
  54. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. 2016. Available online: http://www.R-project.org/ (accessed on 6 March 2021).
  55. Andrade, G.; Bertsch, D.J.; Gazzinelli, A.; King, C.H. Decline in infection-related morbidities following drug-mediated reductions in the intensity of Schistosoma infection: A systematic review and meta-analysis. PLoS Negl. Trop. Dis. 2017, 11, e0005372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Thétiot-Laurent, S.A.-L.; Boissier, J.; Robert, A.; Meunier, B. Schistosomiasis chemotherapy. Angew. Chem. Int. Ed. 2013, 52, 7936–7956. [Google Scholar] [CrossRef]
  57. Faust, C.L.; Osakunor, D.N.; Downs, J.A.; Kayuni, S.; Stothard, J.R.; Lamberton, P.H.; Reinhard-Rupp, J.; Rollinson, D. Schistosomiasis control: Leave no age group behind. Trends Parasitol. 2020, 36, 582–591. [Google Scholar] [CrossRef]
  58. Lo, N.C.; Bezerra, F.S.M.; Colley, D.G.; Fleming, F.M.; Homeida, M.; Kabatereine, N.; Kabole, F.M.; King, C.H.; Mafe, M.A.; Midzi, N.; et al. Review of 2022 WHO guidelines on the control and elimination of schistosomiasis. Lancet Infect. Dis. 2022, 22, e327–e335. [Google Scholar] [CrossRef]
  59. Satayathum, S.A.; King, C.; Muchiri, E.M.; Ouma, J.H.; Whalen, C.C. Factors affecting infection or reinfection with Schistosoma haematobium in coastal Kenya: Survival analysis during a nine-year, school-based treatment program. Am. J. Trop. Med. Hyg. 2006, 75, 83–92. [Google Scholar] [CrossRef] [Green Version]
  60. Feng, W. Progress of research on molluscicidal effect of nereistoxin pesticide. Chin. J. Schistosomiasis Control 2012, 19, 482. [Google Scholar]
  61. Seligman, P.F.; Maguire, R.J.; Lee, R.F.; Hinga, K.R.; Valkirs, A.O.; Stang, P.M. Persistence and fate of tributyltin in aquatic ecosystems. In Organotin; Springer: Dordrecht, The Netherlands, 1996; Volume 20, pp. 429–457. [Google Scholar] [CrossRef]
  62. Sokolow, S.H.; Wood, C.L.; Jones, I.J.; Swartz, S.J.; Lopez, M.; Hsieh, M.H.; Lafferty, K.D.; Kuris, A.M.; Rickards, C.; De Leo, G.A. Global assessment of schistosomiasis control over the past century shows targeting the snail intermediate host works best. PLoS Negl. Trop. Dis. 2016, 10, e0004794. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Sprecher, S.L.; Getsinger, K.D. Zebra Mussel Chemical Control Guide; ERDC/EL TR-00-1; U.S. Army Engineer Research and Development Center: Vicksburg, MS, USA, 2000; 114p. [Google Scholar]
  64. Zayed, K.M.; Guo, Y.-H.; Lv, S.; Zhang, Y.; Zhou, X.-N. Molluscicidal and antioxidant activities of silver nanoparticles on the multi-species of snail intermediate hosts of schistosomiasis. PLoS Negl. Trop. Dis. 2022, 16, e0010667. [Google Scholar] [CrossRef]
  65. Castle, G.D.; Mills, G.A.; Gravell, A.; Jones, L.; Townsend, I.; Cameron, D.G.; Fones, G.R. Review of the molluscicide metaldehyde in the environment. Environ. Sci. Water Res. Technol. 2017, 3, 415–428. [Google Scholar] [CrossRef] [Green Version]
  66. Gohar, A.A.; Maatooq, G.T.; Gadara, S.R.; Aboelmaaty, W.S.; El-Shazly, A.M. Molluscicidal activity of the methanol extract of Callistemon viminalis (Sol. ex Gaertner) G. Don ex Loudon fruits, bark and leaves against Biomphalaria alexandrina snails. Iran. J. Pharm. Res. 2014, 13, 505. [Google Scholar]
Figure 1. Transmission dynamics of schistosomiasis in (a) human definitive hosts (susceptible human S h ; infected human I h ), (b) schistosomiasis parasite forms (cercariae C f ; parasite eggs E h ; miracidia M f ), and (c) intermediate host snails (susceptible snail S v ; and infected snail I v ). The blue dotted arrows represent the interaction between free-living schistosomiasis parasites and the respective hosts. All parameters are as defined in Section 2.1, and the control parameters are shown in blue.
Figure 1. Transmission dynamics of schistosomiasis in (a) human definitive hosts (susceptible human S h ; infected human I h ), (b) schistosomiasis parasite forms (cercariae C f ; parasite eggs E h ; miracidia M f ), and (c) intermediate host snails (susceptible snail S v ; and infected snail I v ). The blue dotted arrows represent the interaction between free-living schistosomiasis parasites and the respective hosts. All parameters are as defined in Section 2.1, and the control parameters are shown in blue.
Mathematics 11 02609 g001
Figure 2. Comparisons between treatment and public literacy targeting the populations of (a) susceptible humans, (b) infected humans, and molluscicide (chemical) and the mechanical strategies targeting (c) susceptible snails and (d) infected snails of schistosomiasis in contrast with no control measures.
Figure 2. Comparisons between treatment and public literacy targeting the populations of (a) susceptible humans, (b) infected humans, and molluscicide (chemical) and the mechanical strategies targeting (c) susceptible snails and (d) infected snails of schistosomiasis in contrast with no control measures.
Mathematics 11 02609 g002
Figure 3. Various molluscicides are regulated by temperature according to their half-lives and durations of chemical exposure to snails. A molluscicide is classified as a non-permanent molluscicide (NP) with a half-life of t 1 / 2 = 21.6 , a persistent molluscicide (P) with a half-life of t 1 / 2 = 50.5 , or a very-permanent molluscicide (VP) with a half-life of t 1 / 2 = 146.5 . (a) Shows variations in the temperature of the environment and the duration of chemical exposure (b) for different chemicals affecting the death rate of targeted host snails. (c) Shows the effect of chemicals with different half-lives on the death rate of different target species (susceptible snails σ s , infected snails σ Ι , cercaria σ c , and miracidia σ m ). At the same time, the effect of combined chemical exposure and temperature on the death rate of host snails and the subsequent reproduction number is shown in (d,e).
Figure 3. Various molluscicides are regulated by temperature according to their half-lives and durations of chemical exposure to snails. A molluscicide is classified as a non-permanent molluscicide (NP) with a half-life of t 1 / 2 = 21.6 , a persistent molluscicide (P) with a half-life of t 1 / 2 = 50.5 , or a very-permanent molluscicide (VP) with a half-life of t 1 / 2 = 146.5 . (a) Shows variations in the temperature of the environment and the duration of chemical exposure (b) for different chemicals affecting the death rate of targeted host snails. (c) Shows the effect of chemicals with different half-lives on the death rate of different target species (susceptible snails σ s , infected snails σ Ι , cercaria σ c , and miracidia σ m ). At the same time, the effect of combined chemical exposure and temperature on the death rate of host snails and the subsequent reproduction number is shown in (d,e).
Mathematics 11 02609 g003
Table 1. The continuously exponentially decreasing mortality rates for intermediate hosts and Schistosoma parasite forms, as functions of the half-life of the chemical t 1 / 2 , the temperature T of the environment in which the chemical is released, and the duration t for which the targeted species are exposed to the chemicals.
Table 1. The continuously exponentially decreasing mortality rates for intermediate hosts and Schistosoma parasite forms, as functions of the half-life of the chemical t 1 / 2 , the temperature T of the environment in which the chemical is released, and the duration t for which the targeted species are exposed to the chemicals.
Intermediate HostsSchistosoma Parasite Forms
Susceptible snailsFree-living miracidia
σ S = ( 1 d v ) e ln 2 t 1 / 2 · e 0.08 20 T t σ M = ( 1 υ 3 ) e ln 2 t 1 / 2 e 0.08 20 T t
Infected snailsFree-living cercaria
σ I = 1 d v δ 2 e ln 2 t 1 / 2 e 0.08 20 T t σ C = ( 1 υ 4 ) e ln 2 t 1 / 2 e 0.08 20 T t
Table 2. Model parameters, their definition, baseline values, potential ranges, and sources.
Table 2. Model parameters, their definition, baseline values, potential ranges, and sources.
ParameterDefinitionBaseline ValueValues Range/DayReferences
Λ 1 Human recruitment rate 4127254–8000[18,49]
Λ 2 Snail recruitment rate 200200[18]
x Initial age of infection in children 730 d 730 d[33,46]
δ 1 The human death rate due to infection0.00390.0039[50]
υ 1 Natural death rate of human0.000040250.0000384–0.0000421[18,41]
ρProportion of stool/urine per person115 g70–160 g[51]
θ h Number of egg parasites in stool/urine262 g−110–513 g−1[51]
ω 1 Miracidia emergence rate 0.002320.00232[52]
υ 2 Natural death rate of IHs 0.011100.004–0.0182[52]
υ 3 Natural death rate of parasite eggs0.071930.001–0.14286[44,46,49]
υ 4 Natural death rate of miracidia0.491650.0833–0.9[46,49]
υ 5 Nautral death rate of cercaria0.0026050.00104–0.00417[52]
N E Number of miracidia released per egg500500[52]
β1Cercaria-human infection rate0.07500.028–0.122[52]
β2Miracidia-snail infection rate0.0012350.000127–0.615[52]
ω 2 Cercaria shedding rate2.62.6[19,49]
δ2Snail death rate due to infection0.0260.002–0.05[52]
K Parasite egg carrying capacity100,000100,000Estimated
CoSaturation coefficient for miracidia infectivity1,000,000 Estimated
MoSaturation coefficient for cercaria infectivity1,000,000 [19]
ε limitation of miracidia the growth velocity 0.250.2–0.3[18,19]
ψ , α , α h , b , α v Effective rates of control strategies0.50–1Varied
σ s , σ I , σ m , σ c Chemical-induced death rates0.50–1Varied
Table 3. The partial rank correlation coefficients for the baseline parameter values for R 0 (without control) as a response function, with the most significant parameters that influence the dynamics of the model highlighted in bold.
Table 3. The partial rank correlation coefficients for the baseline parameter values for R 0 (without control) as a response function, with the most significant parameters that influence the dynamics of the model highlighted in bold.
Parameter Λ 1 Λ 2 x δ 1 υ 1 ρ θ h ω 1 υ 2 υ 3
R 0 +0.8338+0.0451+0.0195−0.0089−0.7899+0.7725+0.7545+0.5504−0.8636−0.7537
Parameter υ 5 N E β1β2 ω 2 δ2CoMo ε υ 4
R 0 −0.5777+0.8186+0.7607+0.7928+0.8005−0.6731−0.7799−0.8032−0.0481−0.8303
Table 4. The variability of R 0 across rows for a single or different combination of control strategies, ψ , α , α h , σ I , σ M , σ C , b , α v ϵ 0 , 1 . The lower the R 0 , the more successful the strategy(s) are in controlling schistosomiasis.
Table 4. The variability of R 0 across rows for a single or different combination of control strategies, ψ , α , α h , σ I , σ M , σ C , b , α v ϵ 0 , 1 . The lower the R 0 , the more successful the strategy(s) are in controlling schistosomiasis.
Control   strategies   as   Functions   of   R 0 Effectiveness   of   Control   Strategies ; ψ , α , α h ,   σ I ,   σ M ,   σ C ,   b , and   α v
00.10.20.30.40.50.60.70.80.91
R 0 ( ψ ) 10,000194713901139988884808748700660643
R 0 ( α α h )10,000990096009100840075006400510036001900975
R 0 ( b α v )10,0009950979895399165866080007141600043593122
R 0 ( σ I , M , C )10,00029513886.961.646.837.230.625.722.020.5
R 0 ( ψ , α α h )10,00019281334103683066351738225212562.7
R 0 ( b α v , σ I , M , C )10,00029413682.956.540.529.821.815.49.66.401
R 0 ( ψ , σ I , M , C )10,00057.419.39.96.14.13.02.31.81.51.3
R 0 ( ψ , b α v )10,000193713621086905766646534420288201
R 0 ( α α h , σ I , M , C ) 10,00029213379.151.735.123.815.69.34.21.9
R 0 ( α α h , b α v )10,00098509406868176996495512036422160828304
R 0 ( ψ , α α h , σ I , M , C ) 10,00029213379.151.735.123.815.69.34.21.9
R 0 ( ψ , α α h , b α v ) 10,00098509406868176996495512036422160828304
R 0 ( ψ , σ I , M , C , b α v ) 10,00057.118.89.45.63.62.41.61.10.630.41
R 0 ( α α h , σ I , M , C , b α v ) 10,00029113075.547.430.419.111.15.51.80.62
R 0 ( ψ , α α h , σ I , M , C , b α v ) 10,00056.618.18.64.72.71.50.830.390.120.04
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tabo, Z.; Kalinda, C.; Breuer, L.; Albrecht, C. Adapting Strategies for Effective Schistosomiasis Prevention: A Mathematical Modeling Approach. Mathematics 2023, 11, 2609. https://doi.org/10.3390/math11122609

AMA Style

Tabo Z, Kalinda C, Breuer L, Albrecht C. Adapting Strategies for Effective Schistosomiasis Prevention: A Mathematical Modeling Approach. Mathematics. 2023; 11(12):2609. https://doi.org/10.3390/math11122609

Chicago/Turabian Style

Tabo, Zadoki, Chester Kalinda, Lutz Breuer, and Christian Albrecht. 2023. "Adapting Strategies for Effective Schistosomiasis Prevention: A Mathematical Modeling Approach" Mathematics 11, no. 12: 2609. https://doi.org/10.3390/math11122609

APA Style

Tabo, Z., Kalinda, C., Breuer, L., & Albrecht, C. (2023). Adapting Strategies for Effective Schistosomiasis Prevention: A Mathematical Modeling Approach. Mathematics, 11(12), 2609. https://doi.org/10.3390/math11122609

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop