Skip to main content

Theory and Modern Applications

Modeling and analysis of taeniasis and cysticercosis transmission dynamics in humans, pigs and cattle

Abstract

Taeniasis and cysticercosis pose a significant challenge to food safety and public health. Cysticercosis reduces the market value for pigs and cattle by making pork and beef unsafe for consumption. In this paper, a mathematical model for the transmission dynamics of taeniasis and cysticercosis in humans, pigs and cattle is formulated and analyzed. The analysis shows that both the disease free equilibrium (DFE) and the endemic equilibrium (EE) exist. To study the dynamics of the diseases, we derived the basic reproduction number \(R_{0}\) by next generation matrix method. When \(R_{0}< 1\), the DFE is globally asymptotically stable whereas when \(R_{0} > 1\) the EE is globally asymptotically stable. The normalized forward sensitivity index was used to determine sensitive parameters to the diseases. Humans’ recruitment rate, probability of humans’ infection with taeniasis and the defecation rate of taenia eggs by humans with taeniasis are the most positive sensitive parameters to diseases’ transmission whereas the human natural death rate is the most negative sensitive parameter. However, it is biologically unethical and not practical to increase human natural mortality rate for disease control. In this case, other parameters with negative sensitivity indices such as death rate of taenia eggs and proportions of unconsumed infected beef and pork can be considered for disease control. Generally, to control the diseases, more efforts should be made directed to reducing the number of humans who have taeniasis and defecate in the open environment. Also meat inspection and indoor keeping of cattle and pigs should be emphasized.

1 Introduction

Taeniasis and cysticercosis are foodborne infections which are caused by the adult and larval form tapeworms, respectively [29]. Three tapeworm species that cause taeniasis in humans and cysticercosis in humans, pigs and cattle are Taenia solium, Taenia saginata and Taenia asiatica [36]. While T. asiatica is endemic only in Asia, other tapeworm species are distributed worldwide [35]. Taeniasis and cysticercosis are associated with poor hygiene and sanitation, open human defecation, free range farming system for pigs and cattle, and poverty [11, 35]. They affect mainly poor and marginalized communities who always have close contact with pigs and cattle, and depend on pork and beef as their sources of food [30, 36].

Humans acquire taeniasis when they consume raw or undercooked beef or pork that contains the tapeworm larval cysts, which develop into adult tapeworms in the human intestine [36]. When humans who are infected with taeniasis defecate in the fields, the tapeworm eggs that are passed out with human faeces spread over the environment by several means, such as water, winds, animal feet and insects, and contaminate the soil, fodder, pastures and water sources [8]. Domestic cattle acquire cysticercosis when they feed on the contaminated environment while pigs acquire cysticercosis through direct consumption of human faeces or indirectly when they ingest taenia eggs from the contaminated environment [29]. When taenia eggs are consumed by cattle and pigs, they hatch, penetrate the intestinal wall and reach the blood circulation system, where they spread throughout the body tissues and organs such as heart, diaphragm, kidney, lungs, liver and tongue, and develop into cysticerci [13, 35]. Human cysticercosis is a result of consuming T. solium eggs from the contaminated environment via contaminated water, fruits and vegetables, or by putting contaminated fingers in the mouth [7]. Once the eggs are ingested, they hatch in the small intestine, and develop into larvae which penetrate the intestinal wall and migrate to various parts of the body such as eyes, muscles, skin and the central nervous system through the blood circulatory system, where they form larval cysts [35]. When the cysts reach and infect the brain, they cause neurocysticercosis which is the most severe form of tapeworm infection of the central nervous system and is the major cause of epilepsy worldwide [36].

Taeniasis and cysticercosis are endemic in many developing countries of Latin America, Africa and Asia [29]. The diseases threaten people’s health and livelihood of subsistence farming communities, posing considerable challenges in food safety and reducing the market value of pigs and cattle by making pork and beef unsafe to eat [34, 35]. Globally, the diseases affect approximately 50 million people and nearly 50,000 people die annually due to cysticercosis [1]. In 2010, the World Health Organization (WHO) listed cysticercosis as a major neglected tropical disease and later in 2015, it was identified as a leading cause of deaths from foodborne diseases. In Tanzania, the disease has been reported to be a serious health problem whereby the parasite is spread in almost all regions [33]. Porcine cysticercosis has been reported to be highly endemic in southern, central and northern regions of Tanzania with the prevalence rate of 0.3–17.4%, 14.9% and 5.5–16.9%, respectively [21, 30]. Human infections have also been reported in the country [15, 23]. For example, in 2012 Tanzania had 17,853 cases and 212 deaths due to epilepsy while porcine cysticercosis cases were 18i and the economic burden for cysticercosis was estimated to be US$ 7.9 million [30].

The diagnosis of human taeniasis is done by examination of stool samples while human cysticercosis involves doing a biopsy of subcutaneous cysts, immunodiagnosis, radiography, computed tomography (CT) scan and magnetic resonance imaging (MRI). The diagnosis of neurocysticercosis requires both central nervous system imaging with CT brain scans or MRI and serological testing [35]. The diagnosis of cysticercosis in pigs and cattle involves meat inspection, serological tests and tongue inspection for pigs [34]. The treatment of human taeniasis and cysticercosis is through administration of prescribed medication of praziquantel, niclosamide, nitazoxanide or albendazole [35]. The treatment of neurocysticercosis may involve prolonged doses of albendazole and/or praziquantel with supporting therapy such as anti-epileptic drugs, corticosteroids, and possibly surgery for some cases [35]. Intervention strategies in pigs and cattle involve using vaccines such as S3Pvac and TSOL18 for pigs and TSA-18 and TSA-9 for cattle [18, 19], and anthelmintic treatment with flubendazole, fenbendazole, oxfendazole, praziquantel and nitazoxanide [34, 35].

The use of mathematical models plays an important role in studying the transmission dynamics of infectious diseases and is very useful for deciding on the appropriate disease control strategies. Over the past two decades, various mathematical models have been formulated and analyzed to study the dynamics and control of parasitic foodborne diseases such as cholera, brucellosis and echinococcus [6, 20, 24, 28, 32, 37, 38]. In particular, some statistical and deterministic models with some stochastic elements that have been formulated and analyzed to study the transmission dynamics and control of taeniasis and cysticercosis in humans and pigs can be found in Gonzalez et al. [12], Kyvsgaard et al. [16], Braae et al. [2], Winskill et al. [34], José et al. [14] and Sánchez-Torreset al. [27]. The study done by Winskill et al. [34] has assumed that human taeniasis occurs when susceptible humans interact with pigs that are infected with cysticercosis something that is practically wrong. Moreover, previous studies did not include the cattle population to study the transmission dynamics of the taeniasis and cysticercosis. To get new insights on the diseases transmission dynamics, we formulate and analyze the mathematical model for the transmission dynamics of taeniasis and cysticercosis by including cattle population and incorporating compartments of infected pork and beef. The cattle population is included due to its significance in accelerating diseases through feeding on a contaminated environment which later affects human beings. Also, the market value of cattle is reduced and thus affecting the livelihood of subsistence farming communities and leading to food safety problems.

This paper is organised as follows: the deterministic model and analysis are presented in Sect. 2, and numerical analysis in Sect. 3. A summary and a conclusion are presented in Sect. 4.

2 Model formulation and analysis

2.1 Model formulation

The basic model for the dynamics of taeniasis and cysticercosis is formulated by modifying the work by Winskill et al. [34] that studied intervention strategies against T. solium cysticercosis, by including cattle population and incorporating compartments of infected pork and beef as described by Kyvsgaard et al. [16]. Taeniasis and cysticercosis use humans as the definitive hosts, and pigs and cattle as the intermediate hosts [8].

The model divides human population into \(S_{H}, I_{HT}\) and \(I_{HC}\) classes that represent susceptible humans, humans infected with taeniasis and humans infected with cysticercosis, respectively. The pig population is divided into \(S_{P}\) and \(I_{P}\) classes that represent susceptible pigs and pigs that are infected with cysticercosis, respectively, while the cattle population is divided into \(S_{C}\) and \(I_{C}\), which represent susceptible cattle and cattle that are infected with cysticercosis, respectively. The compartments \(P_{I}\) and \(B_{I}\) represent infected pork and beef, respectively, and \(E_{T}\) represent the number of taenia eggs in the environment.

Susceptible humans are recruited through birth at per capita rate ψ and diminish through consumption of infected raw or insufficiently cooked pork or beef at rates \(\alpha _{p}\) and \(\alpha _{b}\), respectively. They also diminish by acquiring cysticercosis through consumption of T. solium eggs from the contaminated environment at a rate θ either in contaminated water, fruits, vegetables or by putting contaminated fingers in the mouth [7]. Humans with cysticercosis increase at a rate θ due to interaction of susceptible humans with T. solium eggs in the contaminated environment and they reduce due to disease induced death at rate \(\mu _{d}\). Humans with taeniasis increase at rates \(\alpha _{P}\) and \(\alpha _{b}\) when susceptible humans feed on infected pork and beef, respectively. The parameter \(\beta _{T}\) denotes the probability of susceptible humans acquiring taeniasis from an infected pork or beef, that is, the per capita rate at which susceptible individuals acquire infection. All human compartments suffer natural death rate \(\mu _{h}\). The number of taenia eggs in the environment grow as a result of open defecation by humans who are infected with taeniasis at a rate ν and decrease due to natural death at a rate \(\mu _{e}\).

Susceptible pigs are recruited at per capita rate Λ due to birth and reduce at a rate \(\gamma _{p}\) when they acquire cysticercosis from the contaminated environment. Similarly, susceptible cattle are recruited at per capita rate ϕ due to birth and diminish at a rate \(\gamma _{c}\) when they acquire cysticercosis from the contaminated environment. Both pigs and cattle are further slaughtered for consumption at rates ρ and σ, respectively. The infected pigs and cattle increase at rates \(\gamma _{p}\) and \(\gamma _{c}\), respectively, as a consequence of susceptible pigs and cattle feeding on contaminated environment. All pigs and cattle classes suffer natural death at rates \(\mu _{p}\) and \(\mu _{c}\), respectively. The infected pork increases when infected pigs are slaughtered at a rate ω and it decreases when consumed by humans at a rate \(\alpha _{p}\) whereas infected beef increases when infected cattle are slaughtered at a rate η and decreases when consumed by humans at a rate \(\alpha _{b}\). The parameters δ and ϵ are the proportions of infected pork and beef unconsumed by susceptible humans.

In model formulation we consider the free range farming system for both pigs and cattle populations, and we do not consider immigration. We assume that humans can be infected by either taeniasis or cysticercosis; the number of taenia eggs consumed by humans, pigs and cattle has negligible effect on the total number of eggs in the environment and that infected humans, pigs and cattle cannot recover naturally without treatment. We further assume that pigs and cattle do not suffer disease induced mortality, that they become carriers for their life and that the rates at which susceptible humans consume infected raw or undercooked pork or beef depend on the amount of infected pork or beef that is present. Humans, pigs and cattle contact rates with taenia eggs in the environment are assumed to be density dependent. The model for the transmission dynamics of taeniasis and cysticercosis in humans, pigs and cattle is summarized by using the flow diagram in Fig. 1.

Figure 1
figure 1

The model flow diagram

The state variables and parameters are summarized in Tables 1 and 2, respectively.

Table 1 Description of the state variables
Table 2 Parameters’ description and their values (unit: \(yr^{-1}\))

Some parameter values in this paper are assumed within realistic ranges due to the fact that only little has been done on this area, the diseases being common in rural areas where there is inadequate or no meat inspection and the treatment is not readily available [22]. The model that describes the transmission dynamics of taeniasis and cysticercosis is governed by the following system of differential equations:

$$\begin{aligned} & \frac{dS_{H}}{dt} = \psi - \beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})S_{H}- \theta S_{H} E_{T} -\mu _{h} S_{H}, \\ &\frac{dI_{HT}}{dt}=\beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})S_{H} -\mu _{h} I_{HT}, \\ &\frac{dI_{HC}}{dt}= \theta S_{H} E_{T} -(\mu _{h}+\mu _{d}) I_{HC}, \\ &\frac{dS_{P}}{dt}=\Lambda -\gamma _{p} S_{P} E_{T} -(\rho +\mu _{p}) S_{P}, \\ \begin{aligned} &\frac{dI_{P}}{dt}=\gamma _{p} S_{P} E_{T} -( \omega +\mu _{p}) I_{P}, \\ &\frac{dP_{I}}{dt}=\omega I_{P}-(\delta +\alpha _{P}) P_{I}, \end{aligned} \\ &\frac{dS_{C}}{dt}=\phi -\gamma _{c} S_{C} E_{T} -(\sigma +\mu _{c}) S_{C}, \\ &\frac{dI_{C}}{dt}=\gamma _{c} S_{C} E_{T} -( \eta +\mu _{c}) I_{C}, \\ &\frac{dB_{I}}{dt}=\eta I_{C}-(\epsilon +\alpha _{b}) B_{I}, \\ &\frac{dE_{T}}{dt}=\nu I_{HT}-\mu _{e} E_{T}, \end{aligned}$$
(1)

with initial conditions

$$\begin{aligned} &S_{H}(0)>0;\qquad I_{HT}(0)\geq 0;\qquad I_{HC}(0) \geq 0;\qquad S_{P}(0)>0;\qquad I_{P}(0) \geq 0;\\ & P_{I}(0)\geq 0;\qquad S_{C}(0)>0;\qquad I_{C}(0)\geq 0;\qquad B_{I}(0)\geq 0\quad \text{and}\quad E_{T}(0)\geq 0. \end{aligned}$$

2.2 Positivity of solutions and invariant region

For the model system (1) to be biologically and epidemiologically meaningful, we need to show that the model solutions are positive and bounded.

2.2.1 Positivity of solutions

From the first equation in the model system (1) for susceptible humans, we have

$$\begin{aligned} &\frac{dS_{H}}{dt} = \psi - \beta _{T}( \alpha _{p} P_{I}+ \alpha _{b} B_{I})S_{H}-\theta S_{H} E_{T} -\mu _{h} S_{H}, \\ &\frac{dS_{H}}{dt} \geq - \bigl( \beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})+\theta E_{T} +\mu _{h} \bigr)S_{H}, \\ &\frac{dS_{H}}{S_{H}} \geq - \bigl( \beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})+\theta E_{T} +\mu _{h} \bigr){dt}, \\ &S_{H}(t) \geq S_{H}(0){e}^{ \int _{0}^{t}- (\beta _{T}( \alpha _{p} P_{I}+ \alpha _{b} B_{I})+\theta E_{T} +\mu _{h} )\,ds} \geq 0, \quad\forall t \geq 0. \end{aligned}$$

Using the same approach, it can be shown that

$$\begin{aligned} &I_{HT}(t)\geq 0;\qquad I_{HC}(t)\geq 0;\qquad S_{P}(t)\geq 0;\qquad I_{P}(t)\geq 0;\qquad P_{I}(t) \geq 0;\\ & S_{C}(t)\geq 0;\qquad I_{C}(t)\geq 0;\qquad B_{I}(t)\geq 0;\qquad E_{T}(t)\geq 0,\quad \forall t \geq 0. \end{aligned}$$

Therefore, all solutions of the model system (1) are positive for all \(t \geq 0\).

2.2.2 Invariant region

To show that the model solutions are bounded, we let the total populations for humans, pigs and cattle be \(H=S_{H}+I_{HT}+I_{HC}\), \(T_{P}=S_{P}+I_{P}\) and \(T_{C}=S_{C}+I_{C}\), respectively. By considering the human population, we have

$$\begin{aligned} \begin{aligned}& \frac{dH}{dt}=\psi - \mu _{h} H-\mu _{d} I_{HC}, \\ &\frac{dH}{dt}+ \mu _{h} H\leq \psi. \end{aligned} \end{aligned}$$
(2)

Integrating throughout, we obtain

$$\begin{aligned} H(t)\leq \frac{\psi }{\mu _{h}}+ \biggl(H(0)- \frac{\psi }{\mu _{h}} \biggr){e}^{-\mu _{h}t}, \end{aligned}$$
(3)

where \(H(0)=S_{H}(0)+I_{HT}(0)+I_{HC}(0)\). Following the same procedures, it can be shown that the populations for pigs and cattle are given by

$$\begin{aligned} T_{P}(t)\leq \frac{\Lambda }{\mu _{p}}+ \biggl(T_{P}(0)- \frac{\Lambda }{\mu _{p}} \biggr){e}^{-\mu _{p}t} \quad\text{and}\quad T_{C}(t) \leq \frac{\phi }{\mu _{c}}+ \biggl(T_{C}(0)-\frac{\phi }{\mu _{c}} \biggr){e}^{-\mu _{c}t}, \end{aligned}$$
(4)

respectively, where \(T_{P}(0)=S_{P}(0)+I_{P}(0)\) and \(T_{C}(0)=S_{C}(0)+I_{C}(0)\). The analysis of the solutions (3) and (4) considers two cases: when \(H(0)>\frac{\psi }{\mu _{h}}\), \(T_{P}(0)>\frac{\Lambda }{\mu _{p}}\), \(T_{C}(0)>\frac{\phi }{\mu _{c}}\) and when \(H(0)<\frac{\psi }{\mu _{h}}\), \(T_{P}(0)<\frac{\Lambda }{\mu _{p}}\), \(T_{C}(0)<\frac{\phi }{\mu _{c}}\), which gives

$$\begin{aligned} \begin{aligned} &H(t)\leq \Phi _{t}=\operatorname{max} \biggl\lbrace \frac{\psi }{\mu _{h}},H(0) \biggr\rbrace , \\ &T_{P}(t)\leq \Pi _{t}=\operatorname{max} \biggl\lbrace \frac{\Lambda }{\mu _{p}},T_{P}(0) \biggr\rbrace , \\ &T_{C}(t)\leq \Psi _{t}=\operatorname{max} \biggl\lbrace \frac{\phi }{\mu _{c}},T_{C}(0) \biggr\rbrace . \end{aligned} \end{aligned}$$
(5)

Since \(H(t)=S_{H}+I_{HT}+I_{HC}\leq \Phi _{t}\), it follows that \(I_{HT}\leq \Phi _{t}\). We need to show that \(E_{T}\) is also bounded. By considering the last equation in model system (1), we have

$$\begin{aligned} \begin{aligned}& \frac{dE_{T}}{dt}= \nu I_{HT}-\mu _{e} E_{T}, \\ &\frac{dE_{T}}{dt}+\mu _{e} E_{T}= \nu I_{HT}, \\ &\frac{dE_{T}}{dt}+\mu _{e} E_{T}\leq \nu \Phi _{t}. \end{aligned} \end{aligned}$$
(6)

Integrating throughout, we obtain

$$\begin{aligned} E_{T}(t)\leq \frac{\nu \Phi _{t}}{\mu _{e}} + \biggl(E_{T}(0)- \frac{\nu \Phi _{t}}{\mu _{e}} \biggr){e}^{-\mu _{e}t}, \end{aligned}$$
(7)

from which we get

$$\begin{aligned} E_{T}(t)\leq \Gamma _{t}= \operatorname{max} \biggl\lbrace \frac{\nu \Phi _{t}}{\mu _{e}}, E_{T}(0) \biggr\rbrace . \end{aligned}$$

Using the same approach for the sixth and ninth equations in model system (1), we get

$$\begin{aligned} P_{I}(t)\leq \Theta _{t}= \operatorname{max} \biggl\lbrace \frac{\omega \Pi _{t} }{(\alpha _{b}+\delta )}, P_{I}(0) \biggr\rbrace \quad\text{and}\quad B_{I}(t)\leq \xi _{t}=\operatorname{max} \biggl\lbrace \frac{\eta \Psi _{t} }{(\epsilon +\alpha _{b})}, B_{I}(0) \biggr\rbrace . \end{aligned}$$

Therefore, all solutions of the model system (1) are positive invariant in the region

$$\begin{aligned} \Omega ={}& \bigl\{ (S_{H},I_{HT},I_{HC},S_{P},I_{P},P_{I},S_{C},I_{C},B_{I},E_{T}) \in \mathbb{R}_{+}^{10}: 0\leq H(t)\leq \Phi _{t}; 0\leq T_{p}(t) \leq \Pi _{t}; \\ &0\leq T_{C}(t) \leq \Psi _{t}; 0\leq E_{T}(t)\leq \Gamma _{t}; 0 \leq P_{I}(t)\leq \Theta _{t}; 0\leq B_{I}(t)\leq \xi _{t} \bigr\} . \end{aligned}$$

All solutions which start at the boundary of region Ω converge to this region. The model system (1) is biologically and epidemiologically meaningful and therefore we can consider the flow generated by the model for analysis.

2.3 Equilibrium states and reproduction number \(R_{0}\)

2.3.1 The disease free equilibrium (\(E^{0}\))

When there are no infections in humans, pigs and cattle populations, we obtain the disease free equilibrium \(E^{0}\) which is given by

$$\begin{aligned} E^{0}= \biggl(\frac{\psi }{\mu _{h}},0,0, \frac{\Lambda }{(\rho +\mu _{p})},0,0, \frac{\phi }{(\sigma +\mu _{c})},0,0,0 \biggr). \end{aligned}$$

2.3.2 The basic reproduction number \(R_{0}\)

The basic reproduction number \(R_{0}\) is the expected number of secondary infections that may occur as a result of introducing one infected individual in a fully susceptible population [9]. When \(R_{0}<1\), the disease clears whereas, when \(R_{0}>1\), the disease persists within the population. In computing \(R_{0}\), we adopt the next generation matrix method as used by Van den Driessche and Watmough [31]. Let \(\mathcal{F}_{i}\) be the new infections in compartment i and \(\mathcal{V}_{i}^{+}\) and \(\mathcal{V}_{i}^{-}\) be the transfer terms in and out of the compartment i, respectively, then the infected classes can be written as

$$\begin{aligned} \frac{dx_{i}}{dt}=\mathcal{F}_{i}(x)- \mathcal{V}_{i}^{+}(x)- \mathcal{V}_{i}^{-}(x). \end{aligned}$$

Using the next generation matrix method, we define \(\mathcal{F}_{i}\) and \(\mathcal{V}_{i}\) by

$$\begin{aligned} \mathcal{F}_{i}= \begin{pmatrix} \beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})S_{H} \\ \theta S_{H} E_{T} \\ \gamma _{p} S_{P} E_{T} \\ 0 \\ \gamma _{c} S_{C} E_{T} \\ 0 \\ 0 \end{pmatrix} ,\qquad \mathcal{V}_{i}= \begin{pmatrix} \mu _{h} I_{HT} \\ (\mu _{h}+\mu _{d}) I_{HC} \\ (\omega +\mu _{p}) I_{P} \\ -\omega I_{P}+(\delta +\alpha _{p}) P_{I} \\ (\eta +\mu _{c}) I_{C} \\ -\eta I_{C}+(\epsilon +\alpha _{b}) B_{I} \\ -\nu I_{HT}+\mu _{e} E_{T} \end{pmatrix} . \end{aligned}$$
(8)

The Jacobian of matrices \(\mathcal{F}_{i}\) and \(\mathcal{V}_{i}\) at the disease free equilibrium \(E^{0}\) are given by

$$\begin{aligned} F=\frac{\partial \mathcal{F}_{i}}{\partial x_{j}} \bigl(E^{0} \bigr), \qquad V= \frac{\partial \mathcal{V}_{i}}{\partial x_{j}} \bigl(E^{0} \bigr). \end{aligned}$$
(9)

The basic reproduction number \(R_{0}\) is given by

$$\begin{aligned} R_{0}=\rho \bigl(FV^{-1} \bigr). \end{aligned}$$
(10)

Using the definitions in (9), we have

$$\begin{aligned} &F= \begin{pmatrix} 0&0&0&\frac{\beta _{T}\alpha _{p}\psi }{\mu _{h}}&0&\frac{\beta _{T}\alpha _{b}\psi }{\mu _{h}}&0 \\ 0&0&0&0&0&0&\frac{\theta \psi }{\mu _{h}} \\ 0&0&0&0&0&0&\frac{\gamma _{p}\Lambda }{(\rho +\mu _{p})} \\ 0&0&0&0&0&0&0 \\ 0&0&0&0&0&0&\frac{\gamma _{c}\phi }{(\sigma +\mu _{c})} \\ 0&0&0&0&0&0&0 \\ 0&0&0&0&0&0&0 \end{pmatrix}\quad \text{and} \\ & V= \begin{pmatrix} \mu _{h}&0&0&0&0&0&0 \\ 0&(\mu _{h}+\mu _{d})&0&0&0&0&0 \\ 0&0&(\omega +\mu _{p})&0&0&0&0 \\ 0&0&-\omega &\delta +\alpha _{p}&0&0&0 \\ 0&0&0&0&\eta +\mu _{c}&0&0 \\ 0&0&0&0&-\eta &(\epsilon +\alpha _{b})&0 \\ -\nu &0&0&0&0&0&\mu _{e} \end{pmatrix}. \end{aligned}$$

The next generation matrix is given by

$$\begin{aligned} FV^{-1}= \begin{pmatrix} 0&0&\frac{\beta _{T}\omega \alpha _{p} S_{H}^{0}}{C_{1}}&\frac{\beta _{T}\alpha _{p} S_{H}^{0}}{(\delta +\alpha _{p})}&\frac{\beta _{T}\eta \alpha _{b} S_{H}^{0}}{C_{2}}&\frac{\beta _{T}\alpha _{b}S_{H}^{0}}{(\epsilon +\alpha _{b})}&0 \\ \frac{\nu \theta S_{H}^{0}}{\mu _{e}\mu _{h}}&0&0&0&0&0&\frac{\theta S_{H}^{0}}{\mu _{e}} \\ \frac{\nu \gamma _{p} S_{P}^{0}}{\mu _{e}\mu _{h}}&0&0&0&0&0&\frac{\gamma _{p}S_{P}^{0}}{\mu _{e}} \\ 0&0&0&0&0&0&0 \\ \frac{\nu \gamma _{c} S_{C}^{0}}{\mu _{e}\mu _{h}}&0&0&0&0&0&\frac{\gamma _{c}S_{C}^{0}}{\mu _{e}} \\ 0&0&0&0&0&0&0 \\ 0&0&0&0&0&0&0 \end{pmatrix}, \end{aligned}$$

where \(C_{1}=(\omega +\mu _{p})(\delta +\alpha _{p}), C_{2}=(\eta +\mu _{c})( \epsilon +\alpha _{b}), S_{H}^{0}=\frac{\psi }{\mu _{h}},S_{P}^{0}= \frac{\Lambda }{(\rho +\mu _{p})}\) and \(S_{C}^{0}=\frac{\phi }{(\sigma +\mu _{c})}\). Using the definition in (10), the basic reproduction number \(R_{0}\) is given by

$$\begin{aligned} R_{0}=\sqrt{R_{HP}+R_{HC}}, \end{aligned}$$
(11)

where

$$\begin{aligned} &R_{HP}= \frac{\beta _{T}\nu \alpha _{p}\gamma _{p}\omega \psi \Lambda }{\mu _{h}^{2}\mu _{e}(\omega +\mu _{p})(\alpha _{p}+\delta )(\mu _{p}+\rho )} \quad\text{and}\\ & R_{HC}= \frac{\beta _{T}\nu \alpha _{b}\gamma _{c}\eta \psi \phi }{\mu _{h}^{2}\mu _{e}(\eta +\mu _{c})(\alpha _{b}+\epsilon )(\mu _{c}+\sigma )}. \end{aligned}$$

\(R_{HP}\) is the partial reproduction number due to interaction of humans and pigs whereas \(R_{HC}\) is the partial reproduction number due to interaction of humans and cattle. To give the biological meaning for \(R_{0}\), we rewrite \(R_{HP}\) and \(R_{HC}\) in the forms

$$\begin{aligned} \begin{aligned}& R_{HP}=\beta _{T} \gamma _{p} \frac{\psi }{\mu _{h}} \frac{\nu }{\mu _{e}}\frac{1}{\mu _{h}} \frac{\Lambda }{(\mu _{p}+\rho )}\frac{\omega }{(\omega +\mu _{p})} \frac{\alpha _{p}}{(\alpha _{p}+\delta )}, \\ &R_{HC}=\beta _{T}\gamma _{c}\frac{\psi }{\mu _{h}} \frac{\nu }{\mu _{e}}\frac{1}{\mu _{h}} \frac{\phi }{(\mu _{c}+\sigma )}\frac{\eta }{(\eta +\mu _{c})} \frac{\alpha _{b}}{(\alpha _{b}+\epsilon )}. \end{aligned} \end{aligned}$$
(12)

The terms in (12) can be interpreted as follows: \(\nu /\mu _{e}\) is the density of taenia eggs released by humans with taeniasis, \(1/\mu _{h}\) is the human life expectancy and \(\beta _{T}\) is the probability of humans to be infected with taeniasis due to consumption of raw or insufficiently cooked pork or beef infected with tapeworm larval cysts. The terms \(\psi /\mu _{h}, \Lambda /(\mu _{p}+\rho )\) and \(\phi /(\mu _{c}+\sigma )\) are initial populations for susceptible humans, susceptible pigs and susceptible cattle, respectively; \(1/(\mu _{p}+\rho )\) and \(1/(\mu _{c}+\sigma )\) are the average times pigs and cattle spend in susceptible classes, respectively; \(1/(\omega +\mu _{p})\) and \(1/(\eta +\mu _{c})\) are the infectious periods of infected pigs and cattle, respectively; \(1/(\alpha _{p}+\delta )\) and \(1/(\alpha _{b}+\epsilon )\) are the average infectious period for infected pork and beef, respectively, whereas \(\omega /(\omega +\mu _{p})\) and \(\eta /(\eta +\mu _{c})\) are the proportions of infected pigs and cattle that are slaughtered for consumption, respectively. The terms \(\alpha _{p}/(\alpha _{p}+\delta )\) and \(\alpha _{b}/(\alpha _{b}+\epsilon ) \) are the proportions of infected pork and beef that are eaten by susceptible humans, respectively, whereas \(\gamma _{p}\) and \(\gamma _{c}\) are the rates at which T. solium eggs and T. sagnata eggs are consumed by pigs and cattle, respectively, from the contaminated environment.

2.3.3 Sensitivity analysis

To determine how sensitive the model parameters are to the diseases’ transmission, we adopt the normalized forward sensitivity index approach as used by Chitnis et al. [5]. If κ is a parameter in the basic reproduction number \(R_{0}\), then the sensitivity index of \(R_{0}\) with respect to κ is given by

$$\begin{aligned} \Gamma _{\kappa }^{R_{0}}=\frac{\partial R_{0}}{\partial \kappa } \times \frac{\kappa }{R_{0}}. \end{aligned}$$
(13)

Using Eq. (13) and parameter values in Table 2, we obtain sensitivity indices for each parameter as shown in Table 3. The positive sign of sensitivity index indicates that an increase or decrease of parameter value while keeping other parameters constant increases or decreases the basic reproduction number \(R_{0}\). The negative sign indicates that an increase or decrease of parameter value causes a decrease or increase in expected new average infection \(R_{0}\). For instance, \(\gamma _{p}=+0.1421\) means that an increase in \(\gamma _{p}\) by 10%, increases \(R_{0}\) by 1.421% and hence the disease transmission; and \(\delta =-0.1375\) means that an increase in δ by 20% causes a decrease in \(R_{0}\) by 2.75% and thus decrease the disease transmission.

Table 3 Sensitivity indices

The most positive sensitive parameters in the model are human’s recruitment (ψ), the probability of humans infection with taeniasis (\(\beta _{T}\)) and the defecation rate by humans who are infected with taeniasis (ν) whereas the most negative sensitive parameter is the natural mortality rate of humans (\(\mu _{h}\)). However, since it is unethical and not practical to increase human natural mortality, other parameters with negative sensitivity indices such as death rate of taenia eggs and proportions of unconsumed infected beef and pork can be considered for disease control.

The sensitivity indices for all parameters in the basic reproduction number \(R_{0}\) are plotted on a bar graph in Fig. 2.

Figure 2
figure 2

Graphical representation of sensitivity indices

2.3.4 The endemic equilibrium (\(E^{*}\))

The endemic equilibrium is a point at which taeniasis and cysticercosis exist in humans, pigs and cattle populations. It is obtained by setting the derivatives in model system (1) equal to zero and solving for state variables. The endemic equilibrium \(E^{*}=(S_{H}^{*},I_{HT}^{*},I_{HC}^{*},S_{P}^{*},I_{P}^{*},P_{I}^{*},S_{C}^{*},I_{C}^{*},B_{I}^{*},E_{T}^{*})\) is given by

$$\begin{aligned} &S_{H}^{*} = \frac{\psi }{ ( \beta _{T}F_{0} E_{T}^{*}+\theta E_{T}^{*} +\mu _{h} ) },\qquad I_{HT}^{*}=\frac{\beta _{T}F_{0}E_{T}^{*}S_{H}^{*}}{\mu _{h}}, \qquad I_{HC}^{*}= \frac{\theta S_{H}^{*} E_{T}^{*}}{(\mu _{h}+\mu _{d})}, \\ &S_{P}^{*}=\frac{\Lambda }{(\gamma _{p} E_{T}^{*} +\rho +\mu _{p})},\qquad I_{P}^{*}= \frac{\gamma _{p}\Lambda E_{T}^{*}}{(\gamma _{p} E_{T}^{*} +\rho +\mu _{p})(\omega +\mu _{p})}, \\ &P_{I}^{*}= \frac{\omega \gamma _{p}\Lambda E_{T}^{*}}{(\gamma _{p} E_{T}^{*} +\rho +\mu _{p})(\omega +\mu _{p})(\delta +\alpha _{p})},\qquad S_{C}^{*}= \frac{\phi }{(\gamma _{c} E_{T}^{*} +\sigma +\mu _{c})}, \\ &I_{C}^{*}= \frac{\gamma _{c}\phi E_{T}^{*}}{(\gamma _{c} E_{T}^{*} +\sigma +\mu _{c})(\eta +\mu _{c})}, \qquad B_{I}^{*}= \frac{\eta \gamma _{c}\phi E_{T}^{*}}{(\gamma _{c} E_{T}^{*} +\sigma +\mu _{c})(\eta +\mu _{c})(\epsilon +\alpha _{b})}, \end{aligned}$$

where

$$\begin{aligned} F_{0}=\frac{\alpha _{p} \omega \gamma _{p}\Lambda }{(\gamma _{p} E_{T}^{*} +\rho +\mu _{p})(\omega +\mu _{p})(\alpha _{p}+\delta )}+ \frac{\alpha _{b} \eta \gamma _{c}\phi }{(\gamma _{c} E_{T}^{*} +\sigma +\mu _{c})(\eta +\mu _{c})(\alpha _{b}+\epsilon )}. \end{aligned}$$

To obtain \(E_{T}^{*}\), we solve for the real roots of the polynomial:

$$\begin{aligned} a_{0} {E_{T}^{*}}^{3}+a_{1}{E_{T}^{*}}^{2}+a_{2} {E_{T}^{*}}+a_{3}=0, \end{aligned}$$
(14)

where

$$\begin{aligned} &a_{0}=1\textgreater 0, \\ &a_{1}=\frac{\alpha _{p}\omega \Lambda \beta _{T}}{\theta (\mu _{p}+\omega )(\alpha _{p}+\delta )}+ \frac{\alpha _{b}\eta \phi \beta _{T}}{\theta (\mu _{c}+\eta )(\alpha _{p}+\epsilon )}+ \frac{(\mu _{c}+\sigma )}{\gamma _{c}}+ \frac{(\mu _{p}+\rho )}{\gamma _{p}}+\frac{\mu _{h}}{\theta } \textgreater 0, \\ &a_{2}=\frac{(\rho +\mu _{p})(\sigma +\mu _{c})}{\gamma _{c} \gamma _{p}}+ \frac{\mu _{h}(\sigma +\mu _{c})}{\theta \gamma _{c}}+ \frac{\mu _{h}(\rho +\mu _{p})}{\theta \gamma _{p}}+ \frac{\alpha _{p}\omega \Lambda \beta _{T}(\sigma +\mu _{c})}{\theta \gamma _{c}(\mu _{p}+\omega )(\alpha _{p}+\delta )}\\ &\phantom{a_{2}=}{}+ \frac{\alpha _{b}\eta \phi \beta _{T}(\rho +\mu _{p})}{\theta \gamma _{p}(\mu _{c}+\eta )(\alpha _{p}+\epsilon )}- \biggl( \frac{\nu \psi \beta _{T}\alpha _{p}\omega \Lambda \gamma _{c}}{\mu _{h}\mu _{e}\gamma _{c}\theta (\mu _{p}+\omega )(\alpha _{p}+\delta )}+ \frac{\nu \psi \beta _{T}\alpha _{b}\eta \phi \gamma _{p}}{\mu _{h}\mu _{e}\gamma _{p}\theta (\mu _{c}+\eta )(\alpha _{p}+\epsilon )} \biggr), \\ &a_{3}=\frac{\mu _{h}(\sigma +\mu _{c})(\rho +\mu _{p})(1+R_{0})(1-R_{0})}{\theta \gamma _{c}\gamma _{p}}. \end{aligned}$$

To analyze the possible number of positive real roots of polynomial (14) when \(R_{0}< 1\) and \(R_{0}> 1\), we adopt the approach in Okosun et al. [25]. Using this approach, the number of possible real roots when \(R_{0}<1\) and \(R_{0}>1\) are summarized in Table 4.

Table 4 Number of possible positive real roots when \(R_{0}< 1\) and \(R_{0}> 1\)

Therefore the model system (1) has a unique endemic equilibrium when \(R_{0} > 1\) as shown in cases 2 and 4. Hence we state the following theorem.

Theorem 1

The model system (1) has a unique endemic equilibrium when the basic reproduction number \(R_{0}> 1\).

2.4 Stability analysis of equilibrium states

2.4.1 The global stability of the disease free equilibrium (\(E^{0}\))

Theorem 2

The disease free equilibrium (\(E^{0}\)) of the model system (1) is globally asymptotically stable when \(R_{0}<1\).

Proof

To analyze the global stability of the disease free equilibrium, we adopt the approach used in Castillo-Chavez et al. [4] and Dumont et al. [10]. Using this method, the system of differential equations (1) is written as

$$\begin{aligned} \begin{aligned} & \frac{dX_{r}}{dt}=B(X_{r}-X_{\mathrm{DFE}})+B_{1}X_{n}, \\ &\frac{dX_{n}}{dt}=B_{2}X_{n}, \end{aligned} \end{aligned}$$
(15)

where \(X_{r}\) and \(X_{n}\) are the non-transmitting and transmitting classes, respectively, \(X_{\mathrm{DFE}}\) is the disease free equilibrium, whereas \(B,B_{1}\) and \(B_{2}\) are the matrices to be computed. Here, we have

$$\begin{aligned} \begin{aligned} &X_{r}= (S_{H},S_{P},S_{C} )^{T}, \\ &X_{\mathrm{DFE}}= \biggl(\frac{\psi }{\mu _{h}},\frac{\Lambda }{\mu _{p}+\rho }, \frac{\phi }{\mu _{c}+\sigma } \biggr)^{T}. \end{aligned} \end{aligned}$$
(16)

Thus,

$$\begin{aligned} X_{r}-X_{\mathrm{DFE}}= \biggl(S_{H}- \frac{\psi }{\mu _{h}},S_{P}- \frac{\Lambda }{\mu _{p}+\rho },S_{C}- \frac{\phi }{\mu _{c}+\sigma } \biggr)^{T}. \end{aligned}$$
(17)

Therefore from Eq. (15), we have

$$\begin{aligned} \begin{pmatrix} \psi -\beta _{T}(\alpha _{p} P_{I}+\alpha _{b} B_{I})S_{H}-\theta S_{H} E_{T}-\mu _{h} S_{H} \\ \Lambda -\gamma _{p} S_{P} E_{T}-(\rho +\mu _{p})S_{P} \\ \phi -\gamma _{c} S_{C} E_{T}-(\sigma +\mu _{c})S_{P} \end{pmatrix} =B \begin{pmatrix} S_{H}-\frac{\psi }{\mu _{h}} \\ S_{P}-\frac{\Lambda }{\mu _{p}+\rho } \\ S_{C}-\frac{\phi }{\mu _{c}+\sigma } \end{pmatrix} +B_{1} \begin{pmatrix} I_{HT} \\ I_{HC} \\ I_{P} \\ P_{I} \\ I_{C} \\ B_{I} \\ E_{T} \end{pmatrix} \end{aligned}$$
(18)

and

$$\begin{aligned} \begin{pmatrix} \beta _{T}(\alpha _{p} P_{I}+\alpha _{b} B_{I})S_{H}-\mu _{h} I_{HT} \\ \theta S_{H} E_{T}-(\mu _{h}+\mu _{d})I_{HC} \\ \gamma _{p} S_{P} E_{T}-(\omega +\mu _{p})I_{P} \\ \omega I_{P}-(\delta +\alpha _{p})P_{I} \\ \gamma _{c} S_{C} E_{T}-(\eta +\mu _{c})I_{C} \\ \eta I_{C}-(\epsilon +\alpha _{b})B_{I} \\ \nu I_{HT}-\mu _{e} E_{T} \end{pmatrix} =B_{2} \begin{pmatrix} I_{HT} \\ I_{HC} \\ I_{P} \\ P_{I} \\ I_{C} \\ B_{I} \\ E_{T} \end{pmatrix}. \end{aligned}$$
(19)

From (18) and (19), the matrices B, \(B_{1}\) and \(B_{2}\) are given by

$$\begin{aligned} &B= \begin{pmatrix} -\mu _{h}&0&0 \\ 0&-(\rho +\mu _{p})&0 \\ 0&0&-(\sigma +\mu _{c}) \end{pmatrix}, \\ &B_{1}= \begin{pmatrix} 0&0&0&-\beta _{T}\alpha _{p}S_{H}&0&-\beta _{T}\alpha _{b}S_{H}&- \theta S_{H} \\ 0&0&0&0&0&0&\gamma _{p}S_{P} \\ 0&0&0&0&0&0&\gamma _{c}S_{C} \end{pmatrix}, \end{aligned}$$

and

$$\begin{aligned} B_{2}= \begin{pmatrix} -\mu _{h}&0&0&\beta _{T}\alpha _{p}S_{H}&0&\beta _{T}\alpha _{b}S_{H}&0 \\ 0&-(\mu _{h}+\mu _{d})&0&0&0&0&\theta S_{H} \\ 0&0&-(\omega +\mu _{p})&0&0&0&\gamma _{p}S_{P} \\ 0&0&\omega &-(\delta +\alpha _{p})&0&0&0 \\ 0&0&0&0&-(\eta +\mu _{c})&0&\gamma _{c}S_{C} \\ 0&0&0&0&\eta &-(\epsilon +\alpha _{b})&0 \\ \nu &0&0&0&0&0&-\mu _{e} \end{pmatrix}. \end{aligned}$$

It can be observed that matrix B has real and negative eigenvalues. Thus, the system

$$\begin{aligned} \frac{dX_{r}}{dt}=B(X_{r}-X_{\mathrm{DFE}})+B_{1}X_{n} \end{aligned}$$
(20)

is globally asymptotically stable at \(X_{\mathrm{DFE}}\).

To prove the stability of \(B_{2}\), we adopt the idea of stable Metzler matrix and apply the lemma in Dumont et al. [10]. A Metzler matrix is a matrix whose of-diagonal elements are non-negative, denoted by \(B_{2}(i,j)\geq 0\), for all \(i\neq j\). Thus, it can be observed that \(B_{2}\) is a Metzler matrix. □

Lemma 1

Let M be a square Metzler matrix written in block form:

$$\begin{aligned} M= \begin{pmatrix} P&Q \\ R&S \end{pmatrix}, \end{aligned}$$

where P and S are square matrices. M is Metzler stable if and only if matrices P and \(S-RP^{-1}Q\) are Metzler stable.

Proof

Comparing the Metzler matrix \(B_{2}\) with a square Metzler matrix M, the matrices P, Q, R and S are defined as

$$\begin{aligned} &P= \begin{pmatrix} -\mu _{h}&0&0 \\ 0&-(\mu _{h}+\mu _{d})&0 \\ 0&0&-(\omega +\mu _{p}) \end{pmatrix},\qquad Q= \begin{pmatrix} \beta _{T}\alpha _{p}S_{H}&0&\beta _{T}\alpha _{b}S_{H}&0 \\ 0&0&0&\theta S_{H} \\ 0&0&0&\gamma _{p}S_{P} \end{pmatrix}, \\ &R= \begin{pmatrix} 0&0&\omega \\ 0&0&0 \\ 0&0&0 \\ \nu &0&0 \end{pmatrix},\qquad S= \begin{pmatrix} -(\delta +\alpha _{p})&0&0&0 \\ 0&-(\eta +\mu _{c})&0&\gamma _{c}S_{C} \\ 0&\eta &-(\epsilon +\alpha _{b})&0 \\ 0&0&0&-\mu _{e} \end{pmatrix}. \end{aligned}$$

Clearly, P is a stable Metzler matrix. After some computations, we obtain

$$\begin{aligned} S-RP^{-1}Q= \begin{pmatrix} -(\delta +\alpha _{p})&0&0&\frac{\omega \gamma _{p} S_{P}}{(\omega +\mu _{p})} \\ 0&-(\eta +\mu _{c})&0&\gamma _{c}S_{C} \\ 0&\eta &-(\epsilon +\alpha _{b})&0 \\ \frac{\nu \beta _{T}\alpha _{p}S_{H}}{\mu _{h}}&0&\frac{\nu \beta _{T}\alpha _{b}S_{H}}{\mu _{h}}&- \mu _{e} \end{pmatrix}. \end{aligned}$$

Thus, \(S-RP^{-1}Q\) is Metzler stable if \(\operatorname{Det} (S-RP^{-1}Q)> 0\). That is,

$$\begin{aligned} \mathcal{W}(\delta +\alpha _{p})\mu _{e}- \frac{\nu \beta _{T} S_{H} (\omega \alpha _{p}\gamma _{p}\mathcal{W}S_{P}+\eta \alpha _{b}\gamma _{c}(\delta +\alpha _{b})(\omega +\mu _{p})S_{C} )}{\mu _{h}(\omega +\mu _{p})}> 0, \end{aligned}$$
(21)

where \(\mathcal{W}=(\epsilon +\alpha _{c})(\eta +\mu _{c})\). Substituting the values for \(S_{H}, S_{P}\) and \(S_{C}\) at disease free equilibrium into (21) and simplifying the expression, we obtain \(1-R_{0}^{2}> 0\), where \(R_{0}\) is given in (11). Therefore, the disease free equilibrium \(E^{0}\) is globally asymptotically stable when \(R_{0} < 1\). □

2.4.2 Global stability of the endemic equilibrium (\(E^{*}\))

Theorem 3

The endemic equilibrium (\(E^{*}\)) for the model system (1) is globally asymptotically stable when \(R_{0}>1\).

Proof

Since the endemic equilibrium exists if and only if \(R_{0} >1\), we adopt the approach in Osman et al. [26] to prove global stability of endemic equilibrium \(E^{*}\). Consider the non-linear Lyapunov function

$$\begin{aligned} \mathcal{L}={}&S_{H}^{*} \biggl( \frac{S_{H}}{S_{H}^{*}}-\ln { \frac{S_{H}}{S_{H}^{*}}} \biggr)+I_{HT}^{*} \biggl( \frac{I_{HT}}{I_{HT}^{*}}-\ln {\frac{I_{HT}}{I_{HT}^{*}}} \biggr)+I_{HC}^{*} \biggl(\frac{I_{HC}}{I_{HC}^{*}}-\ln {\frac{I_{HC}}{I_{HC}^{*}}} \biggr) \\ &{}+S_{P}^{*} \biggl(\frac{S_{P}}{S_{P}^{*}}-\ln { \frac{S_{P}}{S_{P}^{*}}} \biggr)+I_{P}^{*} \biggl( \frac{I_{P}}{I_{P}^{*}}-\ln {\frac{I_{P}}{I_{P}^{*}}} \biggr)+P_{I}^{*} \biggl(\frac{P_{I}}{P_{I}^{*}}-\ln {\frac{P_{I}}{P_{I}^{*}}} \biggr) +S_{C}^{*} \biggl(\frac{S_{C}}{S_{C}^{*}}-\ln { \frac{S_{C}}{S_{C}^{*}}} \biggr) \\ &{}+I_{C}^{*} \biggl( \frac{I_{C}}{I_{C}^{*}}-\ln {\frac{I_{C}}{I_{C}^{*}}} \biggr)+B_{I}^{*} \biggl(\frac{B_{I}}{B_{I}^{*}}-\ln {\frac{B_{I}}{B_{I}^{*}}} \biggr) +E_{T}^{*} \biggl(\frac{E_{T}}{E_{T}^{*}}-\ln { \frac{E_{T}}{E_{T}^{*}}} \biggr), \end{aligned}$$
(22)

The time derivative of the Lyapunov function \(\mathcal{L}\) is

$$\begin{aligned} \frac{d\mathcal{L}}{dt}={}& \biggl(1-\frac{S_{H}^{*}}{S_{H}} \biggr) \frac{dS_{H}}{dt}+ \biggl(1-\frac{I_{HT}^{*}}{I_{HT}} \biggr) \frac{dI_{HT}}{dt}+ \biggl(1-\frac{I_{HC}^{*}}{I_{HC}} \biggr) \frac{dI_{HC}}{dt}+ \biggl(1-\frac{S_{P}^{*}}{S_{P}} \biggr) \frac{dS_{P}}{dt} \\ &{}+ \biggl(1- \frac{I_{P}^{*}}{I_{P}} \biggr) \frac{dI_{P}}{dt}+ \biggl(1-\frac{P_{I}^{*}}{P_{I}} \biggr) \frac{dP_{I}}{dt}+ \biggl(1-\frac{S_{C}^{*}}{S_{C}} \biggr) \frac{dS_{C}}{dt}+ \biggl(1-\frac{I_{C}^{*}}{I_{C}} \biggr) \frac{dI_{C}}{dt} \\ &{}+ \biggl(1- \frac{B_{I}^{*}}{B_{I}} \biggr) \frac{dB_{I}}{dt}+ \biggl(1-\frac{E_{T}^{*}}{E_{T}} \biggr) \frac{dE_{T}}{dt}. \end{aligned}$$
(23)

Substituting the equations of model system (1) into Eq. (23), we have

$$\begin{aligned} \frac{d\mathcal{L}}{dt}={}& \biggl(1-\frac{S_{H}^{*}}{S_{H}} \biggr) ( \psi - \beta _{T}\alpha _{p} P_{I}S_{H}- \beta _{T}\alpha _{b} B_{I}S_{H}- \theta S_{H} E_{T} -\mu _{h} S_{H} ) \\ &{} + \biggl(1- \frac{I_{HT}^{*}}{I_{HT}} \biggr) \bigl(\beta _{T}(\alpha _{p} P_{I}+ \alpha _{b} B_{I})S_{H} -\mu _{h} I_{HT} \bigr) \\ &{}+ \biggl(1- \frac{I_{HC}^{*}}{I_{HC}} \biggr) \bigl(\theta S_{H} E_{T} -(\mu _{h}+ \mu _{d}) I_{HC} \bigr) \\ &{} + \biggl(1-\frac{S_{P}^{*}}{S_{P}} \biggr) \bigl(\Lambda -\gamma _{p} S_{P} E_{T} -( \rho +\mu _{p}) S_{P} \bigr)+ \biggl(1- \frac{I_{P}^{*}}{I_{P}} \biggr) \bigl(\gamma _{p} S_{P} E_{T} -(\omega +\mu _{p}) I_{P} \bigr) \\ &{} + \biggl(1- \frac{P_{I}^{*}}{P_{I}} \biggr) \bigl(\omega I_{P}-(\delta +\alpha _{p}) P_{I} \bigr)+ \biggl(1-\frac{S_{C}^{*}}{S_{C}} \biggr) \bigl( \phi -\gamma _{c} S_{C} E_{T} -( \sigma +\mu _{c}) S_{C} \bigr) \\ &{} + \biggl(1- \frac{I_{C}^{*}}{I_{C}} \biggr) \bigl(\gamma _{c} S_{C} E_{T} -(\eta +\mu _{c}) I_{C} \bigr)+ \biggl(1- \frac{B_{I}^{*}}{B_{I}} \biggr) \bigl( \eta I_{C}-(\epsilon + \alpha _{b}) B_{I} \bigr) \\ &{}+ \biggl(1-\frac{E_{T}^{*}}{E_{T}} \biggr) ( \nu I_{HT}- \mu _{e} E_{T} ). \end{aligned}$$
(24)

Equation (24) can also be written as

$$\begin{aligned} \frac{d\mathcal{L}}{dt}={}& \psi - \beta _{T}\alpha _{p} P_{I}S_{H}- \beta _{T}\alpha _{b} B_{I}S_{H}-\theta S_{H} E_{T} -\mu _{h} S_{H} \\ &{}- \frac{\psi S_{H}^{*}}{S_{H}}+ \beta _{T}\alpha _{p} P_{I}S_{H}^{*} +\beta _{T}\alpha _{b} B_{I}S_{H}^{*}+ \theta E_{T}S_{H}^{*} + \mu _{h} S_{H}^{*}+\beta _{T}\alpha _{p} P_{I}S_{H} \\ &{} + \beta _{T} \alpha _{b} B_{I}S_{H} -\mu _{h} I_{HT}- \frac{\beta _{T}\alpha _{p} P_{I}S_{H}I_{HT}^{*}}{I_{HT}} - \frac{ \beta _{T}\alpha _{b} B_{I}S_{H} I_{HT}^{*}}{I_{HT}} +\mu _{h} I_{HT}^{*}+\theta S_{H} E_{T} \\ &{} -(\mu _{h}+\mu _{d}) I_{HC}- \frac{\theta S_{H} E_{T}I_{HC}^{*}}{I_{HC}} +(\mu _{h}+\mu _{d}) I_{HC}^{*}+ \Lambda -\gamma _{p} S_{P} E_{T} -(\rho +\mu _{p}) S_{P} \\ &{}- \frac{\Lambda S_{P}^{*}}{S_{P}}+\gamma _{p} E_{T}S_{P}^{*} +( \rho +\mu _{p}) S_{P}^{*}+\gamma _{p} S_{P} E_{T} -(\omega +\mu _{p}) I_{P}- \frac{\gamma _{p} S_{P} E_{T}I_{P}^{*}}{I_{P}} \\ &{} +(\omega +\mu _{p}) I_{P}^{*} +\omega I_{P}-(\delta +\alpha _{p}) P_{I}- \frac{\omega I_{P} P_{I}^{*}}{P_{I}}+(\delta +\alpha _{p}) P_{I}^{*}+ \phi -\gamma _{c} S_{C} E_{T} \\ &{}-(\sigma +\mu _{c}) S_{C}- \frac{\phi S_{C}^{*}}{S_{C}}+\gamma _{c} E_{T}S_{C}^{*} +( \sigma +\mu _{c})S_{C}^{*}+\gamma _{c} S_{C} E_{T} -(\eta +\mu _{c}) I_{C} \\ &{}- \frac{\gamma _{c} S_{C} E_{T}I_{C}^{*}}{I_{C}}+(\eta +\mu _{c}) I_{C}^{*} +\eta I_{C}-(\epsilon +\alpha _{b}) B_{I}- \frac{\eta I_{C} B_{I}^{*}}{B_{I}}+(\epsilon +\alpha _{b}) B_{I}^{*}+ \nu I_{HT} \\ &{}-\mu _{e} E_{T}-\frac{\nu I_{HT}E_{T}^{*}}{E_{T}}+\mu _{e} E_{T}^{*}. \end{aligned}$$
(25)

Similarly, Eq. (25) can be written as

$$\begin{aligned} \frac{d\mathcal{L}}{dt}=\mathcal{J}-\mathcal{P}, \end{aligned}$$
(26)

where

$$\begin{aligned} \mathcal{J}={}& \psi + \beta _{T}\alpha _{p} P_{I}S_{H}^{*}+ \beta _{T}\alpha _{b} B_{I}S_{H}^{*}+\theta E_{T}S_{H}^{*}+\mu _{h} S_{H}^{*}+ \beta _{T}\alpha _{p} P_{I}S_{H} + \beta _{T}\alpha _{b} B_{I}S_{H}\\ &{}+ \mu _{h} I_{HT}^{*} +\theta S_{H} E_{T} +(\mu _{h}+\mu _{d}) I_{HC}^{*}+ \Lambda +\gamma _{p} E_{T}S_{P}^{*} +(\rho +\mu _{p}) S_{P}^{*}+ \gamma _{p} S_{P} E_{T} \\ &{}+(\omega +\mu _{p}) I_{P}^{*} +\omega I_{P}+( \delta +\alpha _{p}) P_{I}^{*}+\phi +\gamma _{c} E_{T}S_{C}^{*} +( \sigma +\mu _{c})S_{C}^{*}+\gamma _{c} S_{C} E_{T}\\ &{}+(\eta +\mu _{c}) I_{C}^{*}+ \eta I_{C}+(\epsilon +\alpha _{b}) B_{I}^{*}+\nu I_{HT}+\mu _{e} E_{T}^{*}, \\ \mathcal{P}={}& \beta _{T}\alpha _{p} P_{I}S_{H}+ \beta _{T}\alpha _{b} B_{I}S_{H}+ \theta S_{H} E_{T}+\mu _{h} S_{H}+\frac{\psi S_{H}^{*}}{S_{H}}+\mu _{h} I_{HT}+\frac{\beta _{T}\alpha _{p} P_{I}S_{H}I_{HT}^{*}}{I_{HT}} \\ &{} +\frac{ \beta _{T}\alpha _{b} B_{I}S_{H} I_{HT}^{*}}{I_{HT}}+( \mu _{h}+\mu _{d}) I_{HC}+ \frac{\theta S_{H} E_{T}I_{HC}^{*}}{I_{HC}}+\gamma _{p} S_{P} E_{T} +( \rho +\mu _{p}) S_{P}\\ &{}+\frac{\Lambda S_{P}^{*}}{S_{P}}+(\omega + \mu _{p}) I_{P}+\frac{\gamma _{p} S_{P} E_{T}I_{P}^{*}}{I_{P}}+ \frac{\omega I_{P} P_{I}^{*}}{P_{I}}+\gamma _{c} S_{C} E_{T}+( \sigma +\mu _{c}) S_{C}\\ &{}+\frac{\phi S_{C}^{*}}{S_{C}}+(\delta + \alpha _{p}) P_{I}+(\eta +\mu _{c}) I_{C}+ \frac{\gamma _{c} S_{C} E_{T}I_{C}^{*}}{I_{C}}+(\epsilon +\alpha _{b}) B_{I}+ \frac{\eta I_{C} B_{I}^{*}}{B_{I}}\\ &{}+\mu _{e} E_{T}+ \frac{\nu I_{HT}E_{T}^{*}}{E_{T}}. \end{aligned}$$

It can be seen from Eq. (26) that if \(\mathcal{J}<\mathcal{P}\) then \(\frac{d\mathcal{L}}{dt}< 0\) and if \(\Omega =\Omega ^{*}\) then \(\frac{d\mathcal{L}}{dt}=0\). Thus, the largest invariant set in Ω is the endemic equilibrium. Hence from the LaSalle invariant principle [17], we can conclude that, as \(t\to \infty \), the solution of the model system (1) approaches the endemic equilibrium when \(R_{0}>1\). Therefore, the endemic equilibrium is globally asymptotically stable in the invariant set Ω if \(\mathcal{J}<\mathcal{P}\). □

3 Numerical simulations

To understand well the dynamics of taeniasis and cysticercosis in humans, pigs and cattle, we simulate the model (1) using parameters from different literature and some are assumed as indicated in Table 2. To obtain initial conditions we consider a village with 5420 susceptible humans, 750 humans with taeniasis, 528 humans with cysticercosis, 1050 susceptible pigs, 620 infected pigs, 1250 susceptible cattle, 850 infected cattle and 1000 taenia eggs in the environment. We quantify pork and beef in terms of the number of infected pigs and cattle that are slaughtered for consumption. The Runge–Kutta order 4 numerical method is used to simulate the model in MATLAB software.

3.1 The fourth order Runge–Kutta method (RK4)

The 4th order Runge Kutta method (RK4) is a numerical method for solving a system of ordinary differential equations. The advantage of using RK4 over other numerical methods is that it is more accurate as it has high-order local truncation error \(O(h^{4})\) of the Taylor methods and eliminates the need of computing and evaluating the derivatives of a function. RK4 approximates the solution of the initial value system of the first order differential equation of the form

$$\begin{aligned} \frac{dy}{dt}=f \bigl(t,y(t) \bigr);\qquad y(t_{0})=y_{0}. \end{aligned}$$

The method uses the initial value of the function to start the algorithm:

$$\begin{aligned} &k_{1}=hf (t_{i},y_{i}, ),\qquad k_{2}=hf \biggl(t_{i}+\frac{h}{2},y_{i}+ \frac{1}{2}k_{1}, \biggr),\qquad k_{3}=hf \biggl(t_{i}+\frac{h}{2},y_{i}+ \frac{1}{2}k_{2}, \biggr),\\ & k_{4}=hf (t_{i+1},y_{i}+k_{3}, ),\qquad y_{i+1}=y_{i}+ \frac{1}{6} (k_{1}+2k_{2}+2k_{3}+k_{4} ), \end{aligned}$$

for each \(i=1,2,3,\ldots,N-1\), where \(h>0\) is the given step size and \(k_{1}, k_{2}, k_{3}, k_{4}\) are constants that are used to eliminate the need for successive nesting in the second variable of \(f(t,y)\) [3].

3.2 Model simulation

The infected cattle and pigs increase initially to their maximum in the first six months and later they decline and remain constant as illustrated in Fig. 3(a) and (b). The decline of infected cattle and pigs is in correspondence with the decline of taenia eggs in environment as shown in Fig. 3(d). This situation can happen if humans with taeniasis have toilets and hence do not shed taenia eggs in the environment. Susceptible cattle and pigs decrease to their lowest following infection by cysticercosis. However, they increase between the first and fourth year and thereafter remain constant.

Figure 3
figure 3

Dynamics of Humans, Pigs, Cattle and Taenia Eggs

In Fig. 3(c), susceptible humans decline rapidly in the first six months. However, at the end of the first year, susceptible humans increase gradually until the third year where they remain constant. Humans with cysticercosis increase rapidly in the first six months and then with gradual increase until the second year where they remain constant. These trends are caused by variation of taenia eggs in the environment as depicted by Fig. 3(d) and consumption of raw or undercooked infected pork and beef. On the other hand, humans who are infected with taeniasis increase with time in the first two years and then increases linearly with time. This can happen if humans have habit of consuming raw or insufficiently cooked infected beef and pork.

3.3 Effect of varying the most sensitive parameters

In this subsection we present numerical simulation by considering the most sensitive parameters to observe how they affect disease transmission in humans, pigs and cattle. From the sensitivity analysis, the most positive sensitive parameters are the human recruitment rate (ψ), probability of humans to be infected with taeniasis (\(\beta _{T}\)) and the defection rate of humans with taeniasis (ν) whereas the most negative sensitive parameter is the natural mortality rate of humans (\(\mu _{h}\)).

3.3.1 Effect of varying human recruitment rate (ψ)

The dynamics of taeniasis and cysticercosis shows that humans with taeniasis and cysticercosis, and infected cattle and pigs will increase in proportion to the human recruitment rate as illustrated in Fig. 4.

Figure 4
figure 4

Impact of varying human recruitment rate on infected populations

3.3.2 Effect of varying defecation rate (ν)

Infected pigs and cattle, and humans who are infected with cysticercosis increase with time as a result of increase in defecation rate as depicted by Fig. 5(a), (b), (d). A different trend can be observed for humans who are infected with taeniasis in Fig. 5(c) where there is a decrease in cases of humans who are infected when the human defecation rate is increased. This is due to the fact that a high defecation rate can result in an expel of a high number of taenia eggs and so the possibility of recovery.

Figure 5
figure 5

Impact of varying defecation rate on infected populations

3.3.3 Effect of varying probability of taeniasis infection \((\beta _{T})\)

In Fig. 6(a), (b), (c); the infected pigs, cattle and humans who are infected with taeniasis show an increase with time when the probability of taeniasis infection is increased. A different trend can be seen for humans who are infected with cysticercosis in Fig. 6(d), which shows the decrease in disease prevalence with increase in probability of taeniasis infection. This is because humans cysticercosis does not depend on the probability of human infection with taeniasis but on the rate at which humans consume T. solium eggs from the contaminated environment.

Figure 6
figure 6

Impact of varying probability of infection on infected populations

3.3.4 Effect of varying human natural death rate \((\mu _{h})\)

In Fig. 7, the results show that an increase in human natural mortality rate has negative impact in disease prevalence for all infected sub-populations. This indicates that human beings play an important role in transmission of the two diseases.

Figure 7
figure 7

Impact of varying human natural death rate on infected populations

4 Conclusion and recommendation

In this paper, a mathematical model for the transmission dynamics of taeniasis and cysticercosis in humans, pigs and cattle is presented and analyzed. The model is well posed since the model solutions are positive and bounded. The disease free and endemic equilibria exist and their stability are investigated. The next generation approach is used to compute the basic reproduction number \(R_{0}\). The analysis has shown that the disease free equilibrium is globally asymptotically stable when basic reproduction number \(R_{0}<1\) while the endemic equilibrium is globally asymptotically stable when \(R_{0}>1\). The effect of the most sensitive parameters in the diseases’ transmission dynamics was assessed. Numerical results indicate that increasing human recruitment rate \((\psi )\) leads to an increased disease prevalence whereas increasing human natural death rate \((\mu _{h})\) reduces disease prevalence in all populations. On the other hand, increasing probability of human infection with taeniasis \((\beta _{T})\) increases the number of infected pigs, cattle and humans with taeniasis and reduces the number of cases for humans with cysticercosis. Similarly, increasing open human defecation rate \((\nu )\) leads to an increased number of cases for infected pigs, cattle and humans with cysticercosis and decreases the number of cases for humans with taeniasis.

To control the diseases, the open human defecation rate should be reduced through use of toilets, especially in rural communities where pigs and cattle are kept under a free range system, treatment of infected individuals and proper cooking of pork and beef. To reduce the rate of transfer of infections from contaminated environment, we recommend pigs and cattle vaccination, indoor keeping of pigs and cattle as well as improvement in hygiene and sanitation. Moreover, infected pigs and cattle should be treated and meat inspection should be promoted to reduce the rate at which humans are infected with taeniasis.

Availability of data and materials

Most of data used in this research were found from different literature and some were assumed.

References

  1. Aung, A.K., Spelman, D.W.: Taenia solium taeniasis and cysticercosis in southeast Asia. Am. J. Trop. Med. Hyg. 94(5), 947–954 (2016)

    Article  Google Scholar 

  2. Braae, U.C., Devleesschauwer, B., Gabriël, S., Dorny, P., Speybroeck, N., Magnussen, P., Torgerson, P., Johansen, M.V.: Cystisim—an agent-based model for taenia solium transmission and control. PLoS Negl. Trop. Dis. 10(12), 1–17 (2016)

    Article  Google Scholar 

  3. Burden, R.L., Faires, J.D.: Numerical Analysis, Brooks. Cole, Belmont (1997)

    Google Scholar 

  4. Castillo-Chavez, C., Song, B.: Dynamical models of tuberculosis and their applications. Math. Biosci. Eng. 1(2), 361–404 (2004)

    Article  MathSciNet  Google Scholar 

  5. Chitnis, N., Hyman, J.M., Cushing, J.M.: Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull. Math. Biol. 70(5), 1272–1296 (2008)

    Article  MathSciNet  Google Scholar 

  6. Cui, J., Wu, Z., Zhou, X.: Mathematical analysis of a cholera model with vaccination. J. Appl. Math. 2014, 1–17 (2014)

    MathSciNet  MATH  Google Scholar 

  7. Del Brutto, O.H.: Human cysticercosis (taenia solium). Trop. Parasitol. 3(2), 100–103 (2013)

    Article  Google Scholar 

  8. Dermauw, V., Dorny, P., Braae, U.C., Devleesschauwer, B., Robertson, L.J., Saratsis, A., Thomas, L.F.: Epidemiology of taenia saginata taeniosis/cysticercosis: a systematic review of the distribution in southern and eastern Africa. Parasites Vectors 11(1), 1–12 (2018)

    Article  Google Scholar 

  9. Diekmann, O., Heesterbeek, J.A.P., Metz, J.A.: On the definition and the computation of the basic reproduction ratio r 0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 28(4), 365–382 (1990)

    Article  MathSciNet  Google Scholar 

  10. Dumont, Y., Chiroleu, F., Domerg, C.: On a temporal model for the Chikungunya disease: modeling, theory and numerics. Math. Biosci. 213(1), 80–91 (2008)

    Article  MathSciNet  Google Scholar 

  11. Flisser, A., Rodríguez-Canul, R., Willingham, A.L. III: Control of the taeniosis/cysticercosis complex: future developments. Vet. Parasitol. 139(4), 283–292 (2006)

    Article  Google Scholar 

  12. Gonzalez, A.E., Gilman, R.H., Garcia, H.H., Lopez, T.: Use of a simulation model to evaluate control programmes against taenia solium cysticercosis. In: Taenia Solium Cysticercosis: From Basic to Clinical Science, pp. 437–448. CAB International, Wallingford (2002)

    Chapter  Google Scholar 

  13. Ito, A., Wandra, T., Yamasaki, H., Nakao, M., Sako, Y., Nakaya, K., Margono, S.S., Suroso, T., Gauci, C., Lightowlers, M.W.: Cysticercosis/taeniasis in Asia and the Pacific. Vector-Borne Zoonot. Dis. 4(2), 95–107 (2004)

    Article  Google Scholar 

  14. José, M.V., Bobadilla, J.R., Sánchez-Torres, N.Y., Laclette, J.P.: Mathematical model of the life cycle of taenia-cysticercosis: transmission dynamics and chemotherapy (part 1). Theor. Biol. Med. Model. 15(1), 1–19 (2018)

    Article  Google Scholar 

  15. Kavishe, M.D., Mkupasi, E.M., Komba, E.V., Ngowi, H.A.: Prevalence and risk factors associated with porcine cysticercosis transmission in Babati District, Tanzania. Livestock Res. Rural Dev. 29(1), 1–16 (2017)

    Google Scholar 

  16. Kyvsgaard, N.C., Johansen, M.V., Carabin, H.: Simulating transmission and control of taenia solium infections using a reed-frost stochastic model. Int. J. Parasitol. 37(5), 547–558 (2007)

    Article  Google Scholar 

  17. LaSalle, J.P.: The Stability of Dynamical Systems, vol. 25. SIAM, Philadelphia (1976)

    Book  Google Scholar 

  18. Lightowlers, M.W.: Eradication of taenia solium cysticercosis: a role for vaccination of pigs. Int. J. Parasitol. 40(10), 1183–1192 (2010)

    Article  Google Scholar 

  19. Lightowlers, M.W., Rolfe, R., Gauci, C.G.: Taenia saginata: vaccination against cysticercosis in cattle with recombinant oncosphere antigens. Exp. Parasitol. 84(3), 330–338 (1996)

    Article  Google Scholar 

  20. Lolika, P.O., Mushayabasa, S.: Dynamics and stability analysis of a brucellosis model with two discrete delays. Discrete Dyn. Nat. Soc. 2018, 1–21 (2018)

    Article  MathSciNet  Google Scholar 

  21. Mkupasi, E., Kilemile, A., Mandike, O., Prosper, L., Ngowi, H.: Porcine cysticercosis—an emerging neglected food-borne parasitic zoonosis in urban settings in Tanzania: need for immediate control strategies. Tanzania Vet. J. 35(1), 112–117 (2018)

    Google Scholar 

  22. Mkupasi, E.M., Ngowi, H.A., Nonga, H.E.: Prevalence of extra-intestinal porcine helminth infections and assessment of sanitary conditions of pig slaughter slabs in Dar es Salaam city, Tanzania. Trop. Anim. Health Prod. 43(2), 417–423 (2011)

    Article  Google Scholar 

  23. Mwanjali, G., Kihamia, C., Kakoko, D.V.C., Lekule, F., Ngowi, H., Johansen, M.V., Thamsborg, S.M., Willingham, A.L. III: Prevalence and risk factors associated with human taenia solium infections in Mbozi District, Mbeya Region, Tanzania. PLoS Negl. Trop. Dis. 7(3), 1–8 (2013)

    Article  Google Scholar 

  24. Nyabadza, F., Aduamah, J.M., Mushanyu, J.: Modelling cholera transmission dynamics in the presence of limited resources. BMC Res. Notes 12(1), 1–8 (2019)

    Article  Google Scholar 

  25. Okosun, K.O., Mukamuri, M., Makinde, D.O.: Global stability analysis and control of leptospirosis. Open Math. 14(1), 567–585 (2016)

    Article  MathSciNet  Google Scholar 

  26. Osman, S., Otoo, D., Sebil, C.: Analysis of listeriosis transmission dynamics with optimal control. Appl. Math. 11(7), 712–737 (2020)

    Article  Google Scholar 

  27. Sánchez-Torres, N.Y., Bobadilla, J.R., Laclette, J.P., José, M.V.: How to eliminate taeniasis/cysticercosis: porcine vaccination and human chemotherapy (part 2). Theor. Biol. Med. Model. 16(1), 1–14 (2019)

    Article  Google Scholar 

  28. Sun, G.-Q., Xie, J.-H., Huang, S.-H., Jin, Z., Li, M.-T., Liu, L.: Transmission dynamics of cholera: mathematical modeling and control strategies. Commun. Nonlinear Sci. Numer. Simul. 45, 235–244 (2017)

    Article  MathSciNet  Google Scholar 

  29. Symeonidou, I., Arsenopoulos, K., Tzilves, D., Soba, B., Gabriël, S., Papadopoulos, E.: Human taeniasis/cysticercosis: a potentially emerging parasitic disease in Europe. Ann. Gastroenterol. 31(4), 406–412 (2018)

    Google Scholar 

  30. Trevisan, C., Devleesschauwer, B., Schmidt, V., Winkler, A.S., Harrison, W., Johansen, M.V.: The societal cost of taenia solium cysticercosis in Tanzania. Acta Trop. 165, 141–154 (2017)

    Article  Google Scholar 

  31. Van den Driessche, P., Watmough, J.: Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 180(1–2), 29–48 (2002)

    Article  MathSciNet  Google Scholar 

  32. Wang, K., Zhang, X., Jin, Z., Ma, H., Teng, Z., Wang, L.: Modeling and analysis of the transmission of echinococcosis with application to Xinjiang Uygur autonomous region of China. J. Theor. Biol. 333, 78–90 (2013)

    Article  MathSciNet  Google Scholar 

  33. WHO: World health organization: epilepsy due to pork tapeworm: a serious public health problem in the united republic of tanzania (2018). https://www.who.int/taeniasis/en/

  34. Winskill, P., Harrison, W.E., French, M.D., Dixon, M.A., Abela-Ridder, B., Basáñez, M.-G.: Assessing the impact of intervention strategies against taenia solium cysticercosis using the epicyst transmission model. Parasites Vectors 10(1), 1–14 (2017)

    Article  Google Scholar 

  35. World Health Organization: WHO/FAO/OIE Guidelines for the Surveillance, Prevention and Control of Taeniosis/Cysticercosis. World Organisation for Animal Health, Paris (2005)

    Google Scholar 

  36. World Health Organization: Taeniasis/cysticercosis fact sheet (2019). https://www.who.int/news-room/fact-sheets/detail/taeniasis-cysticercosis

  37. Wu, L., Song, B., Du, W., Lou, J.: Mathematical modelling and control of echinococcus in Qinghai province, China. Math. Biosci. Eng. 10(2), 425–444 (2013)

    Article  MathSciNet  Google Scholar 

  38. Yang, Y., Feng, Z., Xu, D., Sandland, G.J., Minchella, D.J.: Evolution of host resistance to parasite infection in the snail–schistosome–human system. J. Math. Biol. 65(2), 201–236 (2012)

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the Ministry of Education, Science and Technology (MoEST) in Tanzania for supporting PhD studies and Mkwawa University College of Education (MUCE) for giving study leave to J.A. Mwasunda.

Funding

The Ministry of Education, Science and Technology (MoEST) supports the corresponding author’s PhD studies at the Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania.

Author information

Authors and Affiliations

Authors

Contributions

JAM: Conceptualization, Model Formulation, Model Analysis and Drafting of the Manuscript; JII: Model Formulation and Supervision; DK: Supervision; and DK: Supervision. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Joshua A. Mwasunda.

Ethics declarations

Competing interests

The authors declare that they have no competing interests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mwasunda, J.A., Irunde, J.I., Kajunguri, D. et al. Modeling and analysis of taeniasis and cysticercosis transmission dynamics in humans, pigs and cattle. Adv Differ Equ 2021, 176 (2021). https://doi.org/10.1186/s13662-021-03341-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13662-021-03341-9

Keywords