 Original Article
 Open Access
 Published:
Taenia solium taeniasis and cysticercosis: extinction or outbreak
Animal Diseases volume 3, Article number: 10 (2023)
Abstract
Taenia solium taeniasis and cysticercosis are neglected zoonotic diseases that affect human health and economies of developing countries. In this work, we formulate and analyze deterministic and continuous time Markov chain (CTMC) stochastic models to determine parameters that drive Taenia solium taeniasis and cysticercosis and the likelihood of their extinction. The basic reproduction number \(R_0\) is computed by the next generation matrix approach, sensitivity index of each parameter in \(R_0\) is derived by the normalized forward sensitivity index and the likelihood of diseases’ extinction is computed by the multitype branching process. The analysis shows that humans with Taenia solium taeniasis, infectious pork and Taenia solium eggs in the environment play an important role in the transmission of Taenia solium taeniasis and cysticercosis, and the model exhibits forward bifurcation at \(R_0=1\). This implies that \(R_0<1\) is a sufficient condition to eliminate Taenia solium taeniasis and cysticercosis. For CTMC model, analysis shows that the probability of Taenia solium taeniasis and cysticercosis extinction is high if the diseases emerge from humans with Taenia solium cysticercosis and there is an outbreak if the diseases emerge from either humans with Taenia solium taeniasis or infectious pork or Taenia solium eggs in the environment. To control Taenia solium taeniasis and cysticercosis, the intervention strategies should focus on improving hygiene and sanitation for reducing shedding rate of Taenia solium eggs in the environment, inspection of pork for reducing the rate of acquiring Taenia solium taeniasis and spraying of insecticides for killing Taenia solium eggs in the environment.
Introduction
Taenia solium taeniasis and cysticercosis are neglected diseases that impose setback to health and economic sectors in developing countries (Mwanjali et al., 2013; Trevisan et al., 2018). The diseases are caused by Taenia solium tapeworm which infects both humans and pigs (Ngowi et al., 2007). The adult stage of Taenia solium tapeworm causes intestinal infection in humans called taeniasis and larval stage infects tissues in both humans and pigs causing cysticercosis (Gebrie and Engdaw, 2015; Mwanjali et al., 2013; Mwidunda et al., 2015; Schmidt et al., 2019). The diseases are mainly associated with poor sanitation and hygiene, pork consumption and free range pig keeping (Mwanjali et al., 2013; Trevisan et al., 2017). Developed countries are also vulnerable to the diseases due to global business, international tourism, and increased migration and refugees (Fernanda et al., 2010; José et al., 2018; Schantz et al., 1998; Sorvillo et al., 2011). Approximately 2.5 million people globally carry adult Taenia solium and at least 20 million are infected with cysticerci of Taenia solium (Pawlowski et al., 2005).
In the dynamics of Taenia solium taeniasis and cysticercosis, pig is an intermediate host and human being is a definitive host (Gebrie and Engdaw, 2015; Maridadi et al., 2011; Trevisan et al., 2017). Interaction of human beings and pigs in the environment maintains the disease life cycle (Trevisan et al., 2017). Humans acquire Taenia solium taeniasis when they eat raw or prebaked pork contaminated with cysticerci (Dixon et al., 2019; Gebrie and Engdaw, 2015; Kyvsgaard et al., 2007; Mwanjali et al., 2013; Mwidunda et al., 2015; Schmidt et al., 2019; Shonyela et al., 2017; Yanagida et al., 2012). In the human intestine, cysticerci develop into adult tapeworms Taenia solium (Sorvillo et al., 2007) whose eggbearing gravid proglottids pass into the environment during defecation (Canadian Paediatric Society, 2019; Gebrie and Engdaw, 2015) and live in the environment for a number of months (Gebrie and Engdaw, 2015). Adult tapeworm carrier can pass millions of eggs on a daily basis either freely in faeces (Schmidt et al., 2019) or as intact segments containing about 250, 000 eggs each (Gebrie and Engdaw, 2015). Pigs contract Taenia solium cysticercosis through direct consumption of human faeces containing Taenia solium eggs or through feeding on the contaminated environment (Mwidunda et al., 2015; Ngowi et al., 2007; Schmidt et al., 2019). The eggs develop oncospheres embryos which enter the intestinal wall and later disseminate to various tissues via bloodstream to form cysticerci (Garcia and Del Brutto, 2000; Sorvillo et al., 2011).
Human beings acquire Taenia solium cysticercosis through ingesting Taenia solium eggs in contaminated food or water (Gebrie and Engdaw, 2015; Mwanjali et al., 2013; Sorvillo et al., 2011; Trevisan et al., 2017; WHO, 2019). The ingested Taenia solium eggs hatch larvae in the intestine which later migrate to various tissues to form cysticerci (Canadian Paediatric Society, 2019). Taenia solium cysticercosis becomes critical when cysticerci migrate to the central nervous system. Infection of central nervous system by cysticerci leads to neurocysticercosis (Del Brutto et al., 2014; García et al., 2003; Mwidunda et al., 2015) which causes neurologic disease (Garcia et al. 2012; Sorvillo et al., 2011; White, 2000) with symptoms such as severe headache, epilepsy, seizure, sightlessness, mental interruption and even fatality (García et al., 2003; Garcia et al., 2016; José et al., 2018; Phiri et al., 2003; Trevisan et al., 2017). Neurocysticercosis is estimated to cause 50, 000 deaths annually (Boa et al., 2006; Garcia and Del Brutto, 2000). In Tanzania, 17, 853 epilepsy cases and 212 deaths due to neurocysticercosis were reported in 2012. It is estimated that 183, 927 pigs (which is \(11.7\%\) of total pigs in Tanzania) were affected with cysticercosis and their market value depreciated to half (Trevisan et al., 2017) or \(60\%\) of uninfected pig (Boa et al., 2006).
Though the impact of Taenia solium taeniasis and cysticercosis is felt among the communities in the developing countries, public knowledge on modes of transmission in many endemic areas is poor (Gebrie and Engdaw, 2015; Mwidunda et al., 2015; Ngowi et al., 2007; Trevisan et al., 2017). This is the main factor that makes human beings to take part in diseases’ transmission (Mwidunda et al., 2015; Ngowi et al., 2007). Lack of knowledge for Taenia solium taeniasis and cysticercosis has led communities in developing world especially Africa to associate epilepsy (which is caused by neurocysticercosis) with witchcraft or demonic possession (Winkler et al., 2010). Studies conducted in Tanzania have shown that the communities’ knowledge gap on the diseases and the modes of transmission is the problem that facilitate diseases’ persistence. Ngowi et al. (2007) have shown that, communities’ poor knowledge on modes of transmission and poor pig management practices are factors which influence people to behaviour that facilitate diseases’ transmission. Other studies also found communities’ knowledge gap (Maridadi et al., 2011; Mwidunda et al., 2015; Shonyela et al., 2017).
To bridge the knowledge gap, deterministic and continuous time Markov chain (CTMC) stochastic models are employed to study the transmission dynamics of Taenia solium taeniasis and cysticercosis. CTMC stochastic models are important tools for studying the dynamics of infectious diseases and determining the likelihood of diseases’ extinction or outbreak. Deterministic models use the basic reproduction number \(R_0\) to predict the disease dynamics and ignore uncertainties that can influence the disease dynamics. In deterministic approach, the disease clears when the basic reproduction number \(R_0<1\) and persists when \(R_0>1\). While deterministic modeling approach ignores uncertainties that influence the disease dynamics, its continuous time Markov chain modeling counterpart takes into account the uncertainties that are inherent to disease transmission and it has the likelihood for disease extinction. Unlike deterministic models, stochastic models can predict extinction of the disease even if the stochastic threshold is greater than unity. Most of mathematical models that are developed to study Taenia solium taeniasis and cysticercosis are deterministic in nature and ignore the influence of uncertainties that are inherent to the diseases’ transmission.
In this work, CTMC models are formulated and analyzed to study dynamics of Taenia solium taeniasis and cysticercosis. The multitype branching process is used to determine the likelihood for diseases’ extinction. The organization of this article is as follows: we begin with deterministic model and its analysis followed by the stochastic model. Numerical analysis is presented before conclusion and recommendations.
Methods and Results
Deterministic model
Development and analysis
Taenia solium taeniasis and cysticercosis model considers human beings, pigs and Taenia solium eggs in the environment. Human population is divided into susceptible \(S_{H}\), humans with Taenia solium taeniasis \(I_{T}\), and humans with Taenia solium cysticercosis \(I_{C}\). Pigs’ population is also divided into susceptible pigs \(S_{P}\) and infected pigs \(I_{P}\). The infectious pork and Taenia solium eggs are represented by M and \(P_{a}\) respectively.
Susceptible humans \(S_{H}\) are replenished at a rate \(\Lambda _{H}\) and contract Taenia solium taeniasis and cysticercosis at respective rates
Humans with Taenia solium taeniasis \(I_{T}\) increase as susceptible individuals acquire Taenia solium taeniasis at a rate \(\lambda _T\) and humans with Taenia solium cysticercosis \(I_{C}\) increase when susceptible humans eat food or drink water which is contaminated with Taenia solium eggs at a rate \(\lambda _C\). Humans with Taenia solium cysticercosis \(I_{C}\) suffer death due to disease at a rate \(\alpha _C\) (Fernanda et al., 2010; Sorvillo et al., 2007) and all human classes suffer death due natural causes at a rate \(\mu _{H}\).
Susceptible pigs \(S_{P}\) are recruited at a constant rate \(\Lambda _{P}\) and diminish by contracting Taenia solium cysticercosis at a rate
and when they are harvested at a rate \(\alpha _P\). Pigs which are Taenia solium cysticercosis carriers \(I_{P}\) increase at a rate \(\Lambda _{P}\) and decrease as they are slaughtered at a rate \(\gamma _P\). All pig classes are assumed to suffer natural mortality at a rate \(\mu _{P}\). Infectious pork M increase following slaughtering of infected pigs at rate \(\gamma _P\). However, it decreases following consumption by susceptible humans at a rate \(\alpha _H\) and decays at a rate \(\alpha _m\) (\(\alpha _m\) is the proportion of infectious pork that is not consumed by human beings). Taenia solium eggs \(P_{a}\) are replenished when humans with Taenia solium taeniasis defecate in the open environment at a rate \(\gamma _E\), and diminish due to natural mortality at a rate \(\mu _E\).
The model assumes that: the number of Taenia solium eggs which are ingested to cause Taenia solium cysticercosis in humans and pigs is negligible thus has no significant impact on their total number in the environment (Winskill et al., 2017). For model symplicity, auto infection and individuals who have both Taenia solium taeniasis and cysticercosis are not considered. Transmission of Taenia solium taeniasis and cysticercosis is density dependent. Susceptible humans who acquire both Taenia solium taeniasis and cysticercosis are not considered. Susceptible pigs are harvested at a constant rate and all pigs are assumed to be kept under free range system. Table 1 described model parameters and interactions of variables are described in Fig. 1.
Using parameter values in Table 1 and state variables in Fig. 1, the model system for Taenia solium taeniasis and cysticercosis is given by:
Positivity and boundedness of solutions
The positivity of model solutions and their bounds are proved by the following theorem:
Theorem 1
For nonnegative initial conditions, the system (3) has nonnegative solutions which are positively invariant and globally attracting in the region
where \(N_H\) and \(N_P\) denote the total populations for humans and pigs respectively.
Proof
To prove Theorem 1, we rewrite the system (3) in the form
where \(y=(y_1, y_2, y_3,y_4,y_5, y_6,y_7)=(S_{H}, I_{C}, I_{T}, S_{P}, I_{P}, M, P_{a})\) and \(h(y)=\left( h_1(y), h_2(y),\cdots , h_7(y)\right)\) represent the rate of change for the corresponding state variables. In this form, the system (5) is;;
We can see that for all \(i=1,\cdots , 7\), \(h_i(y)\ge 0\) if \(y\in [0,\infty )^7\) and \(y_i=0\). Since humans, pigs and Taenia solium eggs populations are nonnegative, then the rate of change for each variable in system (5) is locally Lipschitz in \(\Pi\) and system (3) has a unique solution.
Now using total humans and pigs populations, we have:
Since \(\frac{dN_H}{dt}\le \Lambda _{H}\mu _{H}N_H\) and \(\frac{dN_P}{dt}\le \Lambda _{P}(\mu _P+\alpha _P)N_P\), then \(\frac{dN_H}{dt}\le 0\) if \(N_H(t)\ge \frac{\Lambda _{H}}{\mu _{H}}\) and \(\frac{dN_P}{dt}\le 0\) if \(N_P(t)\ge \frac{\Lambda _{P}}{\alpha _P+\mu _P}\) . By standard comparison theorem, it is clear that: \(N_H(t)\le \frac{\Lambda _{H}}{\mu _{H}}\) when \(N_H(0)\le \frac{\Lambda _{H}}{\mu _{H}}\) and \(N_P(t)\le \frac{\Lambda _{P}}{\alpha _P+\mu _P}\) when \(N_P(0)\le \frac{\Lambda _{P}}{\alpha _P+\mu _P}\).
Considering infectious pork and Taenia solium eggs, we can show that
\(M(t)\le \frac{\gamma _P\Lambda _{P}}{\alpha _P+\mu _P}\) when \(M(0)\le \frac{\gamma _P\Lambda _{P}}{\alpha _P+\mu _P}\) and \(P_{a}(t)\le \frac{\gamma _E\Lambda _{H}}{\mu _{H}}\) when \(P_{a}(0)\le \frac{\gamma _E\Lambda _{H}}{\mu _{H}}\).
Thus the region:
is positive invariant and the solutions of the model system (3) enter the region in finite time and the flow generated can be considered for analysis.
Equilibrium states and basic reproduction number \(R_0\)
In the absence of Taenia solium taeniasis and cysticercosis, the system has disease free equilibrium:
To determine whether Taenia solium taeniasis and cysticercosis clear or persist, the basic reproduction number \(R_0\) that measures the average new infections is computed by next generation matrix approach (Van den Driessche and Watmough, 2002; Winskill et al., 2017). Taenia solium taeniasis and cysticercosis clear when \(R_0<1\) and persist when \(R_0>1\) (Diekmann et al., 1990; Van den Driessche and Watmough, 2002). When \(R_0<1\), the average new infections are less than one thus Taenia solium taeniasis and cysticercosis die out (Van den Driessche and Watmough, 2002) while when \(R_0>1\), average new infections are greater than one therefore the diseases persist.
If we define the vectors for new infections and transition terms from model (3) by \(F_i\) and \(V_i\) respectively, then the basic reproduction number \(R_0\) is given as the maximum eigenvalue of the next generation matrix \(FV^{1}\) (Van den Driessche and Watmough, 2002), that is:
where
From equation (9), the basic reproduction number \(R_0\) is given by:
The basic reproduction number \(R_0\) depends on the proportion at which infected pigs are slaughtered, the proportion of infectious pork that decays, the rate at which Taenia solium eggs are shed in the environment, the rates of acquiring Taenia solium taeniasis by humans and Taenia solium cysticercosis by pigs, the life expectancy for human, pig and Taenia solium eggs, human initial population and the rate at which pigs are recruited. The relative impact of each parameter in the basic reproduction number \(R_0\) is determined by sensitivity analysis.
Sensitivity analysis
The relative impact of each parameter in the basic reproduction number \(R_0\) is determined by normalized forward sensitivity index. If \(\beta _i\) is a parameter in the basic reproduction number \(R_0\), then its relative impact is measured by sensitivity index given by:
Using Eq. (11), the sensitivity indices for parameters in the basic reproduction number \(R_0\) are summarized in Table 2. Analysis shows that, the average new infections will increase if pigs are recruited and acquire Taenia solium cysticercosis, humans acquire Taenia solium taeniasis, and when Taenia solium eggs are shed in the environment. However, the average new infections will diminish if natural mortality for pigs and Taenia solium eggs increases, a large proportion of infectious pork decays and susceptible pigs are highly harvested.
Analysis of equilibrium states
We obtain the equilibrium states by setting the right handside of model system (3) equal to zero and solve for the state variables. The model system (3) is therefore rewritten as:
Finding the equilibrium states, we solve for \(S_{H}^*\), \(I_{C}^*\), \(S_{P}^*\), \(I_{P}^*\), \(M^*\) and \(P_{a}^*\) in terms of \(I_{T}^*\) to obtain:
Substituting \(S_{H}^*\) and \(M^*\) in Eq. (14), we obtain the polynomial in \(I_{T}^{*}\) given by
where
The polynomial (19) has two solutions which are \(I_{T}^{*}=0\) and \(AI_{T}^{*2}+BI_{T}^{*}+C=0\). The solution \(I_{T}^{*}=0\) represents the disease free equilibrium which is stable when \(R_0<1\) and \(AI_{T}^{*2}+BI_{T}^{*}+C=0\) represents endemic equilibrium which is stable when \(R_0>1\).
The number of roots for \(AI_{T}^{*2}+BI_{T}^{*}+C=0\) depends on the signs of A, B and C. Using the general formula, the roots that represent endemic equilibria are given by
The model system (3) has two endemic equilibria if \(B^2\ge 4AC\) and \(R_0>1\). The endemic point \(I_{T}^{*}=\frac{B\sqrt{B^24AC}}{2A}\) is not stable as we do not have negative infected pigs while the endemic point \(I_{T}^{*}=\frac{B+\sqrt{B^24AC}}{2A}\) is stable when \(R_0>1\). The model system (3) undergoes forward bifurcation at \(R_0=1\) as illustrated in Fig. 2. Therefore \(R_0<1\) is a sufficient condition to eliminate Taenia solium taeniasis and cysticercosis.
Stochastic model
Disease dynamics are exposed to influences that cannot explicitly modeled (Ditlevsen and Samson, 2010). Modeling of such dynamics needs consideration of uncertainties which are inherent to diseases’ transmission. CTMC stochastic models capture well uncertainties that are inherent to diseases’ transmission by considering the movement of individuals from one epidemiological class to another as discrete (Lloyd et al., 2007; Maliyoni et al., 2017). In a discrete movement of individuals from one epidemiological class to another, there is a possibility of infectious individuals to die or recover before an outbreak occurs (Lahodny and Allen, 2013; Lloyd et al., 2007; Maliyoni, 2021; Maliyoni et al., 2017). This fact is not captured by deterministic models which rely on the basic reproduction number \(R_0\) to conclude the extinction or persistence of the disease. The inability to take into account the possibility of an infectious individual to die or recover from the disease before an outbreak occurs shows that results from deterministic models need to be complemented with results from stochastic models. In this research, we develop CTMC model and apply multitype branching process to compute the probability of diseases’ extinction.
Continuous Time Markov Chain (CTMC) Model
Formulation of a CTMC model uses the notations for the state variables and parameters that were used in the deterministic model (3). The time is assumed to be continuous and the random variables \(S_{H}\), \(I_{C}\), \(I_{T}\), \(S_{P}\), \(I_{P}\), M, \(P_{a}\) are discrete with a finite space
where N is a positive integer that represents the maximum size of the population in the finite space (Allen, 2017).
Let
be the associated random vector for all discretevalued random variables. The model is assumed to be time homogeneous and satisfies Markov chain property that the future state of the process at \(t+\bigtriangleup t\) depends on the current state t, and time between events is exponentially distributed with parameter (Maliyoni, 2021; Maliyoni et al., 2017)
The time interval \(\bigtriangleup t\) is sufficiently small (Allen, 2017) so that the transition probabilities of the process are within the interval [0, 1]. We further assume that at most one event takes place during the time interval \(\bigtriangleup t\) (Maliyoni et al., 2017). The state transition rates for the continuous time Markov chain model are defined in Table 3. The positive, negative and zero values denote increase by one, decrease by one and no change respectively.
Multitype branching process and stochastic threshold
The multitype branching process is used to approximate nonlinear CTMC model near the disease free equilibrium (Maliyoni et al., 2017). The theory of the multitype branching process is applied to infected classes \(I_{C}\), \(I_{T}\), \(I_{P}\), M and \(P_{a}\). \(S_{H}\) and \(S_{P}\) are assumed to be at disease free equilibrium i.e \(S_{H}(0)\approx N_H(0)\) and \(S_{P}(0)\approx N_P(0)\), where \(N_H(0)\) and \(N_P(0)\) are initial populations for humans and pigs respectively.
The infection of susceptible humans or pigs is referred to as a birth of an offspring and we say that, the infectious individual of type i, \(I_i\) has given birth to an infectious individual of type j, \(I_j\). The number of offspring produced by an individual of type i is independent of the number of offspring produced by individuals of type j, \(j\ne i\), and infectious individuals of type i have the same offspring probability generating function (pgf). With linearity of multitype branching process near disease free equilibrium and timehomogeneous with independent births and deaths, we can define offspring pgf for birth and death of the infectious humans and pigs, and use them to compute probability of disease extinction or major outbreak.
Let \(X_{ji}\) denote the offspring random variable of type j produced by infectious individuals of type i, then the offspring pgf associated with birth of infectious individuals of type i or death of infectious individuals \(I_j\) provided the process begins with one infectious individual of type j, \(I_j=1\), \(I_i=0\) for \(i\ne j\) (Dorman et al., 2004; Pénisson, 2010) is;
where
is the probability that one individual of type i gives birth to individual \(r_j\). Using Eq. (24), we have \(n\times n\) non negative expectation matrix \(\textbf{M}=[m_{ji}]\), where \(m_{ji}\) is the expected number of infectious offspring of type j generated by infectious individual of type i (Maliyoni et al., 2017) defined by
The probability of disease extinction or outbreak can be determined by finding the size of the spectral radius \(\rho (\textbf{M})\) of the matrix \(\textbf{M}\). When \(\rho (\textbf{M})<1\) then the probability of ultimate disease extinction is one, that is
We apply multitype branching process to define the offspring probability generating functions for all infected classes. If there is one human with Taenia solium cysticercosis \(I_{C}(0) = 1\) and there are no infectives in other clasess (\(I_{T}(0) = 0, I_{P} (0) = 0, M(0) = 0, P_{a}(0) = 0\)), then there is no chance for the diseases to spread. Therefore, the initially infected human can only die of Taenia solium cysticercosis or die naturally. Thus the offspring probability generating function is
If there is one human with Taenia solium taeniasis \(I_{T}(0)=1\) and there are no infectives for other clasess \((I_{C}(0) = 0, I_{P} (0) = 0, M(0) = 0, P_{a}(0) = 0)\) then the offspring probability generating function is given by
The term \(\gamma _E /(\gamma _E+\mu _{H})\) is the probability that the initial human that is infected with Taenia solium taeniasis sheds Taenia solium eggs in the environment at a rate \(\gamma _E\) and does not die thus resulting into one human that is infected with Taenia solium taeniasis and one Taenia solium egg in the environment whereas the term \(\mu _{H}/(\gamma _E+\mu _{H})\) is the probability that the initial human that is infected with Taenia solium taeniasis dies before shedding Taenia solium eggs in the environment thus resulting into zero human with Taenia solium taeniasis and zero Taenia solium eggs in the environment.
When there is one infected pig \(I_{P}(0)=1\) and there are no infectives for other clasess (\(I_{C}(0) = 0, I_{T} (0) = 0, M(0) = 0, P_{a}(0) = 0\)) then the offspring probability generating function is
The term \(\gamma _P /(\gamma _P+\mu _P)\) defines the probability that the initial infected pig is slaughtered resulting into no infected pig and the presence of infectious pork whereas \(\mu _P/(\gamma _P+\mu _P)\) is the probability that, the initial pig with Taenia solium cysticercosis dies before it is slaughtered thus resulting into zero infected pig and infectious pork.
If there is an infectious pork \(M(0)=1\) and there are no infectives for other classes (\(I_{C}(0) = 0, I_{T}(0) = 0, I_{P}(0) = 0, P_{a}(0) = 0\)) then the offspring probability generating function is
The term \(\alpha _H\beta _T /(\alpha _H\beta _T+\alpha _m)\) is the probability that a proportion of infectious pork is consumed by human resulting into one human with Taenia solium taeniasis and the presence of infectious pork whereas the term \(\alpha _m /(\alpha _H\beta _T+\alpha _m)\) is the probability that the initial infectious pork decays resulting into no infectious pork and human with Taenia solium taeniasis.
If initially there is one Taenia solium egg in the environment \(P_{a}(0)=1\) and there are no infectives for other clasess (\(I_{C}(0) = 0, I_{T}(0) = 0, I_{P}(0) = 0, M(0) = 0\)) then the offspring probability generating function is
The term \(r\beta _C / (r\beta _C+c_P\beta _P+\mu _E)\) is the probability that the initial Taenia solium egg from the environment infect a susceptible human resulting into one human with Taenia solium cysticercosis and Taenia solium egg. The term \(c_P\beta _P /(r\beta _C+c_P\beta _P+\mu _E)\) is the probability that the initial Taenia solium egg from the environment infect a susceptible pig resulting into one infected pig and Taenia solium egg whereas \(\mu _E /(r\beta _C+c_P\beta _P+\mu _E)\) is the probability that the initial Taenia solium egg dies before it infects a human or pig resulting into zero infected human and pig.
The expectation matrix \(\mathbb {M}\) is given by
The stochastic threshold is given by \(\rho (\textbf{M})\). The disease goes to entinction if \(\rho (\textbf{M})<1\). This fact is related to the basic reproduction \(R_0\) that when \(R_0<1\) the disease dies out in the population. While in deterministic models the disease persists when \(R_0>1\), in CTMC stochastic models there are chances for the disease to die or persist even if \(\rho (\textbf{M})>1\) depending on the initial number of infectives in the population. Therefore if \(\rho (\textbf{M}) > 1\), there exists a fixed point \((q_1,q_2,q_3,q_4,q_5)\in (0,1)^5\) of the offspring generating functions (27)(31) which is used to calculate the likelihood of disease extinction. The fixed point is obtained by setting \(f_i(q_1,q_2,q_3,q_4,q_5) = q_i\) for \(i= 1, 2, 3, 4,5\). That is:
Thus, the likelihood of disease outbreak for type i infectives is
where in for \(n=1,2,3,4,5\) are the initial values for human with cysticercosis, human with taeniasis, infected pigs, infectious pork and Taenia solium eggs that were exposed in a susceptible population respectively.
Numerical Analysis
In this section, we study the dynamics of Taenia solium taeniasis and cysticercosis using parameter values in Table 1. Most of the parameter values are assumed due to the fact that, Taenia solium taeniasis and cysticercosis are neglected diseases which are not much studied and therefore, have no recorded data. We simulate the model system (3) and compare the results with continuous time Markov chain stochastic model.
The effects of Taenia solium taeniasis and cysticercosis in humans and pigs appear after five years. Susceptible humans decrease as they contract Taenia solium taeniasis and cysticercosis whereas humans with Taenia solium taeniasis and cysticercosis increase as shown in Fig. 3A. The steady state is attained after 20 years. The same trend is observed in pig population. The effects of Taenia solium cysticercosis in human and pig populations are observed in the first 20 years where steady state is attained as illustrated in Fig. 3B. The number of Taenia solium eggs in the environment increases when humans with Taenia solium taeniasis defecate in the open spaces while infectious pork increases when infected pigs are slaughtered as in Fig. 4.
Using 10000 sample paths, we concurrently simulate the deterministic and continuous time Markov chain models to compare their results. Susceptible humans and pigs are assumed to be at a disease free equilibrium.
Figures 4, 5 and 6 show that deterministic and stochastic results do not deviate significantly from each other and the number of Taenia solium taeniasis and cysticercosis cases are at the peak in the first 20 and 15 years respectively. From 15 years, susceptible humans and pigs decrease rapidly as depicted in Fig. 7. Generally, the results show that as Taenia solium taeniasis and cysticercosis invade human population, endemic level is attained in the first twenty years and stabilize thereafter. The solutions for both deterministic and CTMC stochastic models predict persistence of Taenia solium taeniasis and cysticercosis as illustrated in Figs. 4, 5 and 6. As long as there are infected pigs and Taenia solium eggs in the environment, Taenia solium taeniasis and cysticercosis will persist if interventions are not administered. To eradicate the diseases, efforts should focus on treating infected humans and pigs, pork inspection and improving hygiene and sanitation.
Probability of Taenia solium taeniasis and cysticercosis Extinction or Outbreak
The fixed point of multitype branching process is applied to compute the probability for Taenia solium taeniasis and cysticercosis extinction \(\textbf{P}_{\textbf{0}}\) when initial conditions are varied. Using a proportion of 10000 sample paths for CTMC stochastic model, we compute numerical probability for Taenia solium taeniasis and cysticercosis extinction \(\textbf{P}_{\textbf{a}}\) and compare \(\textbf{P}_{\textbf{0}}\) and \(\textbf{P}_{\textbf{a}}\) results. The results from Table 4 show that \(\textbf{P}_{\textbf{a}}\) is a good estimate of \(\textbf{P}_{\textbf{0}}\). The fixed point due to the offspring probability generating function in \((0,~1)^5\) is \((q_1,q_2,q_3,q_4,q_5)=(1,0.0420,0.5729,0.2507,0.0814)\). This reveals that, the probability for Taenia solium taeniasis and cysticercosis extinction is high if the diseases emerge from humans with Taenia solium cysticercosis or infected pigs. However, there is major outbreak for Taenia solium taeniasis and cysticercosis if the diseases emerge from either humans with taeniasis or infectious pork or Taenia solium eggs in the environment as shown in Fig. 8.
Conclusion
Taenia Solium taeniasis and cysticercosis are neglected diseases which affect poor communities in developing countries. Deterministic and CTMC stochastic models were formulated to study the dynamics of the diseases. The next generation matrix approach was used to compute the basic reproduction number \(R_0\) and forward normalized sensitivity index is applied in deriving sensitivity indices to determine parameters that drive the diseases. The multitype branching process was applied to a CTMC stochastic model to determine the probability of diseases’ extinction.
The analysis shows that the model undergoes forward bifurcation at \(R_0=1\), implying that \(R_0<1\) is a sufficient condition to eliminate Taenia Solium taeniasis and cysticercosis. Thus the disease free equilibrium is globally asymptotically stable when \(R_0<1\) and endemic equilibrium is globally asymptotically stable when \(R_0>1\). Sensitivity results show that humans with taeniasis, infectious pork and Taenia solium eggs in the environment play an importan role in the persistence of Taenia solium taeniasis and cysticercosis. However, the natural mortality for pigs and Taenia solium eggs in the environment, and the rate at which infectious pork decays decrease the average new infections for Taenia solium taeniasis and cysticercosis.
The results for CTMC stochstic model reveal that there is a high likelihood of Taenia solium taeniasis and cysticercosis extinction if they emerge from either humans with Taenia solium cysticercosis or infected pigs. However, there is a major outbreak if the diseases emerge from humans with Taenia solium taeniasis or infectious pork or Taenia solium eggs in the environment. Simulation shows that the results for deterministic and CTMC stochastic models do not deviate from each other. Both models show that, humans with Taenia solium taeniasis, infectious pork and Taenia solium eggs in the environment play an important role in the persistence of Taenia solium taeniasis and cysticercosis. To effectively control Taenia solium taeniasis and cysticercosis, intervention strategies should focus on treatment of infected humans and pigs, improving hygiene and sanitation, conduct pork inspection and spray of insecticides for killing Taenia solium eggs in the environment.
Availability of data and materials
This work used data from the literature to simulate the model.
Abbreviations
 CTMC:

Continuous Time Markov Chain
 pgf:

probability generation function
References
Allen, L.J. 2010. An introduction to stochastic processes with applications to biology. CRC Press.
Boa, M., E. Mahundi, A. Kassuku, and A. Willingham III, and N.C. Kyvsgaard. 2006. Epidemiological survey of swine cysticercosis using antemortem and postmortem examination tests in the southern highlands of tanzania. Veterinary Parasitology 139 (1–3): 249–255.
Canadian Paediatric Society. 2019. Cysticercosis and Taeniasis  Medical Conditions  Caring for Kids New to Canada. https://www.kidsnewtocanada.ca/conditions/cysticercosis. Accessed 14 Mar 2020.
Del Brutto, O.H., H.H. García, et al. 2014. Cysticercosis of the human nervous system. Springer.
Diekmann, O., J.A.P. Heesterbeek, and J.A. Metz. 1990. On the definition and the computation of the basic reproduction ratio r 0 in models for infectious diseases in heterogeneous populations. Journal of Mathematical Biology 28 (4): 365–382.
Ditlevsen, S., and A. Samson. 2013. Introduction to stochastic models in biology. Stochastic Biomathematical Models: with Applications to Neuronal Modeling, 3–35.
Dixon, M.A., U.C. Braae, P. Winskill, M. Walker, B. Devleesschauwer, S. Gabriel, and M.G. Basanez. 2019. Strategies for tackling taenia solium taeniosis/cysticercosis: A systematic review and comparison of transmission models, including an assessment of the wider taeniidae family transmission models. PLoS Neglected Tropical Diseases 13 (4): 0007301.
Dorman, K.S., J.S. Sinsheimer, and K. Lange. 2004. In the garden of branching processes. SIAM Review 46 (2): 202–229.
Fernanda, M., S. Celeste, M. Liliana, F. Victoria, R. Dominique, and T. Alcides. 2010. Neurocysticercosis causing sudden death: a case report. Asian Pacific Journal of Tropical Medicine 3 (8): 669–670.
Garcia, H.H., and O.H. Del Brutto. 2000. Taenia solium cysticercosis. Infectious disease clinics of North America 14 (1): 97–119.
García, H.H., A.E. Gonzalez, C.A. Evans, R.H. Gilman, C. W. G. in Peru, et al. 2003. Taenia solium cysticercosis. Lancet 362 (9383): 547–556.
Garcia, H.H., A.E. Gonzalez, V.C. Tsang, S.E. O’Neal, F. LlanosZavalaga, G. Gonzalvez, J. Romero, S. Rodriguez, L.M. Moyano, V. Ayvar, et al. 2016. Elimination of Taenia solium transmission in northern peru. New England Journal of Medicine 374 (24): 2335–2344.
Garcia, H.H., S. Rodriguez, R.H. Gilman, A.E. Gonzalez, V.C. Tsang, and C. W. G. in Peru. 2012. Neurocysticercosis: is serology useful in the absence of brain imaging? Tropical Medicine & International Health 17 (8): 1014–1018.
Gebrie, M., and T.A. Engdaw. 2015. Review on taeniasis and its zoonotic importance. European Journal of Applied Sciences 7 (14): 182–191.
José, M.V., J.R. Bobadilla, N.Y. SánchezTorres, and J.P. Laclette. 2018. Mathematical model of the life cycle of taeniacysticercosis: transmission dynamics and chemotherapy (part 1). Theoretical Biology and Medical Modelling 15 (1): 18.
Kyvsgaard, N.C., M.V. Johansen, and H. Carabin. 2007. Simulating transmission and control of Taenia solium infections using a reedfrost stochastic model. International Journal for Parasitology 37 (5): 547–558.
Lahodny, G.E., and L.J. Allen. 2013. Probability of a disease outbreak in stochastic multipatch epidemic models. Bulletin of mathematical biology 75 (7): 1157–1180.
Lloyd, A.L., J. Zhang, and A.M. Root. 2007. Stochasticity and heterogeneity in hostvector models. Journal of the Royal Society Interface 4 (16): 851–863.
Maliyoni, M. 2021. Probability of disease extinction or outbreak in a stochastic epidemic model for west nile virus dynamics in birds. Acta Biotheoretica 69 (2): 91–116.
Maliyoni, M., F. Chirove, H.D. Gaff, and K.S. Govinder. 2017. A stochastic tickborne disease model: Exploring the probability of pathogen persistence. Bulletin of mathematical biology 79 (9): 1999–2021.
Maridadi, A., J. Lwelamira, and F. Simime. 2011. Knowledge and practices related to t. solium cysticercosistaeniasis among smallholder farmers in selected villages in kilolo district in iringa region in southern highlands of tanzania. International Journal of Animal and Veterinary Advances 3(3):196–201.
Mwanjali, G., C. Kihamia, D. V. C. Kakoko, F. Lekule, H. Ngowi, M. V. Johansen, S. M. Thamsborg, and A. L. Willingham III. 2013. Prevalence and risk factors associated with human Taenia solium infections in mbozi district, mbeya region, tanzania. PLoS Neglected Tropical Diseases 7 (3): e2102.
Mwasunda, J.A., J.I. Irunde, D. Kajunguri, and D. Kuznetsov. 2021. Modeling and analysis of taeniasis and cysticercosis transmission dynamics in humans, pigs and cattle. Advances in Difference Equations 2021 (1): 1–23.
Mwasunda, J.A., J.I. Irunde, D. Kajunguri, and D. Kuznetsov. 2022. Outbreak or extinction of bovine cysticercosis and human taeniasis: A stochastic modelling approach. Applied Mathematical Modelling 106: 73–85.
Mwidunda, S. A., H. Carabin, W. B. Matuja, A. S. Winkler, and H. A. Ngowi. 2015. A school based cluster randomised health education intervention trial for improving knowledge and attitudes related to Taenia solium cysticercosis and taeniasis in mbulu district, northern tanzania. PLoS ONE 10 (2): e0118541.
Ngowi, H., E. Tolma, A. Kassuku, M. Mlozi, J. Mlangwa, H. Carabin, and A. Willingham Iii. 2007. Using the precede model to plan a health promotion strategy for control of Taenia solium infections in northern tanzania. International Journal of Health Promotion and Education 45 (2): 41–48.
Pawlowski, Z., J. Allan, and E. Sarti. 2005. Control of Taenia solium taeniasis/cysticercosis: from research towards implementation. International Journal for Parasitology 35 (11–12): 1221–1232.
Pénisson, S. 2010. Conditional limit theorems for multitype branching processes and illustration in epidemiological risk analysis.
Phiri, I.K., H. Ngowi, S. Afonso, E. Matenga, M. Boa, S. Mukaratirwa, S. Githigia, M. Saimo, C. Sikasunge, N. Maingi, et al. 2003. The emergence of Taenia solium cysticercosis in eastern and southern africa as a serious agricultural problem and public health risk. Acta Tropica 87 (1): 13–23.
Schantz, P. M., P. P. Wilkins, and V. C. Tsang. 1998. Immigrants, imaging, and immunoblots: the emergence of neurocysticercosis as a significant public health problem. In Emerging infections 2, 213–242. American Society of Microbiology.
Schmidt, V., M.C. O’Hara, B. Ngowi, K.H. Herbinger, J. Noh, P. P. Wilkins, V. Richter, C. Kositz, W. Matuja, and A. S. Winkler. 2019. Taenia solium cysticercosis and taeniasis in urban settings: Epidemiological evidence from a healthcenter based study among people with epilepsy in dar es salaam, tanzania. PLoS Neglected Tropical Diseases 13 (12): e0007751.
Shonyela, S.M., E.M. Mkupasi, S.C. Sikalizyo, E.M. Kabemba, H.A. Ngowi, and I. Phiri. 2017. An epidemiological survey of porcine cysticercosis in nyasa district, ruvuma region, tanzania. Parasite Epidemiology and Control 2 (4): 35–41.
Sorvillo, F.J., C. DeGiorgio, and S.H. Waterman. 2007. Deaths from cysticercosis, united states. Emerging Infectious Diseases 13 (2): 230.
Sorvillo, F., P. Wilkins, S. Shafir, and M. Eberhard. 2011. Public health implications of cysticercosis acquired in the united states. Emerging Infectious Diseases 17 (1): 1.
Trevisan, C., B. Devleesschauwer, V. Schmidt, A.S. Winkler, W. Harrison, and M.V. Johansen. 2017. The societal cost of taenia solium cysticercosis in tanzania. Acta tropica 165: 141–154.
Trevisan, C., B. Devleesschauwer, N. Praet, A. Pondja, Y.A. Assane, P. Dorny, S.M. Thamsborg, P. Magnussen, and M.V. Johansen. 2018. Assessment of the societal cost of taenia solium in angónia district, mozambique. BMC Infectious Diseases 18 (1): 1–11.
Van den Driessche, P., and J. Watmough. 2002. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Mathematical biosciences 180 (1–2): 29–48.
White, A.C., Jr. 2000. Neurocysticercosis: updates on epidemiology, pathogenesis, diagnosis, and management. Annual Review of Medicine 51 (1): 187–206.
WHO. 2019. Taeniasis/cysticercosis. https://www.who.int/newsroom/factsheets/detail/taeniasiscysticercosis. Accessed 14 Mar 2020.
Winkler, A.S., M. Mayer, S. Schnaitmann, M. Ombay, B. Mathias, E. Schmutzhard, and L. JilekAall. 2010. Belief systems of epilepsy and attitudes toward people living with epilepsy in a rural community of northern tanzania. Epilepsy & Behavior 19 (4): 596–601.
Winskill, P., W.E. Harrison, M.D. French, M.A. Dixon, B. AbelaRidder, and M.G. Basáñez. 2017. Assessing the impact of intervention strategies against Taenia solium cysticercosis using the epicyst transmission model. Parasites & vectors 10 (1): 73.
Yanagida, T., Y. Sako, M. Nakao, K. Nakaya, and A. Ito. 2012. Taeniasis and cysticercosis due to Taenia solium in japan. Parasites & Vectors 5 (1): 18.
Acknowledgements
Authors thank Dr. J.A. Mwasunda and Mr. Mussa Amos for proof reading this work.
Funding
There is no external funding for this work.
Author information
Authors and Affiliations
Contributions
All authors have contributed equally in this work. The author(s) read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Authors have agreed to send this work for publication.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Irunde, J.I., Luhanda, F.B. Taenia solium taeniasis and cysticercosis: extinction or outbreak. Animal Diseases 3, 10 (2023). https://doi.org/10.1186/s44149023000711
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s44149023000711
Keywords
 Taeniasis
 Cysticercosis
 Basic reproduction number
 Sensitivity analysis
 Multitype branching process
 Stochastic model