# African Trypanosomiasis Dynamics: Modelling the Effects of Treatment, Education, and Vector Trapping

Authors: Yustina A. Liana; Nyimvua Shaban; Goodluck Mlay; Anitha Phibert
Journal: International Journal of Mathematics and Mathematical Sciences
 (2020)

---

## Abstract

African trypanosomiasis is a vector-borne disease that is mainly transmitted by infected tsetse flies. A deterministic model of tsetse fly vector, human, and cattle hosts is formulated and analyzed to gain insights into the disease dynamics. The roles of public health education, treatment, and tsetse fly traps are studied. The effective reproduction number, a threshold used to determine whether the disease persists or dies out in the population, is determined. The sensitivity analysis of the model parameters is performed to determine their relationship with the effective reproduction number. The results show that the tsetse fly biting rate is the most sensitive parameter to the effective reproduction number. Furthermore, the model’s numerical simulation shows that a combination of all three interventions has the most significant impact on the control of African trypanosomiasis. Thus, we recommend that these control measures be put concurrently in endemic areas for effective control of the disease transmission.


---

## Body

## 1. Introduction

African trypanosomiasis is a disease caused by microscopic parasites of the speciesTrypanosoma brucei, and it is transmitted through bites of infected tsetse flies of the genus Glossina which are most common in woodland and savannah areas of sub-Saharan Africa. Trypanosoma brucei infects both humans and animals, and if not treated early, it can lead to death [1]. Human African trypanosomiasis (HAT) is caused by Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense while African animal trypanosomiasis (AAT) is mainly caused by Trypanosoma brucei vivax, Trypanosoma brucei congolense, and Trypanosoma brucei rhodesiense [2]. Even though both humans and animals are infected by Trypanosoma brucei, cattle are mostly infected because of tsetse flies’ feeding preferences. Both male and female tsetse flies can transmit Trypanosoma and depend only on hosts’ blood to survive or for all their nutritional needs, unlike other vector-borne diseases such as malaria, where only a female mosquito can feed on blood and can transmit the disease [3].The disease has affected at least 37 countries in sub-Saharan Africa, threatening the lives of millions of people in rural areas. Around 10 million square kilometers in sub-Saharan Africa have been affected byTrypanosoma brucei species [3]. On average, 70,000 cases of HAT are reported each year in sub-Saharan Africa, and more than 1 million cattle die every year due to trypanosomiasis and cause an economic loss of between $2 and $4.5 billion annually [4, 5]. It is claimed that the high mortality rate of livestock can lead to low production of meat and milk up to 50% in sub-Saharan Africa every year [6].So far, much work has been done to investigate the dynamics of African trypanosomiasis. For example, Moore et al. [4] developed a model on the effect of climate change on African trypanosomiasis dynamics. The result from their study predicted that, by the year 2090, about 46–77 million people would be exposed to trypanosomiasis disease. Otieno et al. [7] studied the dynamics of trypanosomiasis in a cattle population by including the wild animals as an alternative feeding source for tsetse flies. The results obtained from the study show that wild animals accelerate the disease in the cattle population. The sensitivity analysis also revealed that the vector biting rate and the vector survival rate are parameters with the greatest influence on the disease spread. These results indicate that the control strategies that target to decrease the contact between the vector and cattle populations would be the best way to eliminate the disease in the population. Kajunguri et al. [8] developed a model to control tsetse flies using insecticide-treated cattle in a multihost population. The results showed that the treatment of both infected humans and cattle combined with insecticide-treated cattle effectively decreases trypanosomes’ prevalence. Meisner et al. [9] researched the role of trypanocide treatment on cattle. The result showed that as the coverage of treated cattle with trypanocide increases, the disease prevalence decreases in both humans and cattle. Ndondo et al. [10] analyzed gambiense sleeping sickness on both human and cattle by including tsetse fly growth from its larval stage to the adult stage. The study’s findings show that human African trypanosomiasis cannot persist in the human population in the absence of cattle.Despite the numerous studies conducted to control African trypanosomiasis, the disease remains a major health threat to both human livelihood and livestock production and affects economic development in Africa. African trypanosomiasis has caused around 500 million farmers in rural Africa villages to live under food shortage and poverty [11]. The African trypanosomiasis model developed by Ndondo et al. [10] considered treatment as the only control strategy while leaving out public health education and tsetse fly traps. Therefore, this paper aimed to extend Ndondo et al.’s [10] work by incorporating public health education and tsetse fly traps to address the question “How does tsetse-fly traps, public health education, and treatment of both human and cattle affect the dynamics of African trypanosomiasis?”. The rest of this paper is structured as follows: in Section 2, the model is formulated, and in Section 3, the model analysis is carried out. In Section 4, sensitivity analysis and its interpretations are made, and in Section 5, we have numerical simulation while concluding remarks are covered in Section 6.


## 2. The Model

The model considered here consists of the submodels of human, cattle, and tsetse fly populations. For the control of trypanosomiasis disease infection, we consider three interventions, namely, public health education, treatment, and trapping of tsetse flies. Human population at timet is subdivided into five subpopulations of uneducated susceptible Sut, educated susceptible Set, exposed Eht, infected Iht, and recovered individuals Rht. The total human population Nht is thus given by Nht = Sut + Set + Eht + Iht + Rht. Human individuals are recruited into the population at a constant rate of Λh. We assume that education strategy is implemented at the rate θ only to susceptible people to make them aware of how to protect themselves from tsetse fly biting (public health education is given on the importance of clearing the environment, wearing long-sleeved clothes, and using repellents). Uneducated susceptible individuals can acquire infection and move to exposed class through the bite of infectious tsetse flies at a rate of λh. Due to the education campaign, it is assumed that only a small fraction of educated individuals move to the exposed class at a rate of 1−ϵλh, where ϵ is the efficacy of the education campaign. Both educated susceptible and uneducated susceptible may also leave their respective classes through natural death at a rate of μh. Exposed humans become infectious at a rate αh, and the infectious humans leave the infected class through natural death, disease-induced death, or recovery at the rates μh,σh, and βh, respectively. Individuals in the recovered class may leave the compartment by either natural death at a rate of μh or through losing temporary immunity and move to the susceptible class at a rate of ψh. We further assume that the infected human recovers through treatment, implying that no human individual recovers naturally, and also, this study assumes that all humans are born susceptible. The model for the human population takes the following form:(1)dSudt=Λh+ψhRh−λhSu−μh+θSu,dSedt=θSu−1−ϵλhSe−μhSe,dEhdt=λhSu+1−ϵλhSe−μh+αhEh,dIhdt=αhEh−μh+σh+βhIh,dRhdt=βhIh−μh+ψhRh,where λh is the force of infection given by λh = (1 − ρ)bkIv/Nh, with k as the probability that the infectious vector infects a susceptible individual, b is the tsetse fly blood-feeding rate per day, and ρ is the proportion of tsetse fly feeding on cattle, and the complimentary (1 − ρ) is the proportion of tsetse fly feeding on a human per day. It is assumed that a tsetse fly can only become infected at its first blood meal and remains so throughout its lifespan [12].The cattle population at timet is subdivided into four compartments of susceptible Sct, exposed Ect, infected Ict, and recovered Rct. The total cattle population denoted by Nct is given as Nct=Sct+Ect+Ict+Rct. At any moment in time, it is assumed that cattle are recruited into the population at a constant rate of Λc. Cattle leave the susceptible class through natural death at a rate of μc or by getting infected and joining the exposed class at λc. The exposed cattle become infectious and move to the infected class at a rate of αc. The infectious cattle leave the infected class through natural death, disease-induced death, or recovery at the rates μc, σc, and βc, respectively. We also assume that cattle acquire temporary immunity, and recovered cattle may leave the recovered class either by natural death at the rate μc or through the waning of temporary immunity and move to susceptible class at a rate ψc. Furthermore, it is assumed in this study that all cattle are born susceptible, and there are no cattle that can experience natural recovery from the disease.The model for cattle population takes the following form:(2)dScdt=Λc+ψcRc−λcSc−μcSc,dEcdt=λcSc−μc+αcEc,dIcdt=αcEc−μc+σc+βcIc,dRcdt=βcIc−μc+ψcRc,where λc is the force of infection given by λc = bdρIv/Nc and d is the probability that an infectious tsetse fly infects cattle.The tsetse fly vector population at timet is divided into three compartments of susceptible Svt, exposed Evt, and infectious Ivt. Therefore, the total vector population is Nvt=Svt+Evt+Ivt. It is assumed that tsetse flies are recruited through birth at a rate of Λv. We assume that the tsetse fly population, regardless of its status, dies naturally at a rate of μv or by being trapped at a rate of ω. The tsetse flies move to exposed class after acquiring the infection at a rate of λv, and also, the infected flies progress to being infectious at a rate of αv. In this study, all tsetse flies are assumed to be born susceptible, and once a tsetse fly becomes infected, it is assumed to remain in that condition throughout its lifespan.The model for vector population takes the following form:(3)dSvdt=∧v−λvSv−μv+ωSv,dEvdt=λvSv−μv+ω+αvEv,dIvdt=αvEv−μv+ωIv,where λv=1−ρbgIh/Nh+ρbzIc/Nc is the force of infection in the vector population.The tsetse fly becomes infected either by biting infected humans or by biting infected cattle. Therefore,g is a probability that a susceptible tsetse fly becomes infected after biting an infected human host and z is the probability that a susceptible tsetse fly becomes infected after biting infected cattle. Tables 1 and 2 summarize the definitions of all state variables and the associated parameters, respectively.Table 1
Descriptions of variables of model (4).
VariableDescriptionsSuSusceptible uneducated humanSeSusceptible educated humanEhExposed humanIhInfected humanRhRecovered humanScSusceptible cattleλhForce of infection for humanλcForce of infection for cattleEcExposed cattleIcInfected cattleRcRecovered cattleSvSusceptible vectorEvExposed vectorIvInfected vectorλvForce of infection for vectorTable 2
Description of parameters of model (4).
ParameterDescriptionsΛhHuman birth rateμhHuman natural death rateϵEducation efficacyψhImmunity waning rate for the recovered humanψcImmunity waning rate for the recovered cattleαhHuman incubation rateαcCattle incubation rateαvVector incubation rateβhHuman recovery rate/treatment rateβcCattle recovery rate/treatment rateθRate at which humans become aware of the diseasezProbability that a vector becomes infected after biting infected cattleσhHuman death rate due to diseaseσcCattle death rate due to diseaseΛcCattle birth rateμcCattle natural death rateμvVector natural death rateΛvVector recruitment rateωVector death rate due to trappingdProbability that infectious vector infects cattlebVector biting rate per daykProbability that infectious vector infects humangProbability that a vector becomes infected by biting infected humanρProportion of vector feeding on cattleFrom the model diagram in Figure1 and the relevant assumptions, we get a full model that describes the dynamics of African trypanosomiasis disease:(4)dEhdt=λhSu+1−ϵλhSe−μh+αhEh,dIhdt=αhEh−μh+σh+βhIh,dRhdt=βhIh−μh+ψhRh,dScdt=Λc+ψcRc−λcSc−μcSc,dEcdt=λcSc−μc+αcEc,dIcdt=αcEc−μc+σc+βcIc,dRcdt=βcIc−μc+ψcRc,dSvdt=Λv−λvSv−μv+ωSv,dEvdt=λvSv−μv+ω+αvEv,dIvdt=αvEv−μv+ωIv.Figure 1
Flow diagram for the dynamics of African trypanosomiasis.The corresponding initial conditions areSu0≥0,Se0≥0,Rh0≥0,Sc0≥0,Eh0≥0,Ih0≥0,Ec0≥0,Ic0≥0, and Rc0≥0,Sv0≥0,Ev0≥0,Iv0≥0.


## 3. Model Analysis

Since model (4) monitors population of humans, cattle, and tsetse flies, all its associated variables are assumed to be nonnegative ∀t≥0. We need to show that model (4) is mathematically and epidemiologically well defined by proving that all the state variables are always positive ∀t≥0.Theorem 1.
Let the feasible region for the three populations beΩ=Ωh∪Ωc∪Ωv∈ℝ+5×ℝ+4×ℝ+3, where Ωh=Su,Se,Eh,Ih,Rh∈ℝ+5,Su+Se+Eh+Ih+Rh=Nh≤Λh/μh, Ωc=Sc,Ec,Ic,Rc,Sc+Ec+Ic+Rc=Nc≤Λc/μc, and Ωv=Sv,Ev,Iv,Sv+Ev+Iv=Nv≤Λv/μv. It is sufficient to consider the solutions in Ω since it is positively invariant and attracting with respect to model (4).Proof.
To prove the feasible region, we use model (4) and compute the total population of humans, cattle, and tsetse flies.
For human population,(5)dNhdt=Λh−μhNh−σhIh+.
For cattle population,(6)dNcdt=Λc−μcNc−σcIc≤Λc−μcNc.
For tsetse fly population,(7)dNvdt=Λv−μv+ωNv≤Λv−μvNv.
Then, solving forNh,Nc, and Nv from equations (5)–(7), we obtain the following inequalities as in [13] (see “theorem”):(8)Nht≤Λhμh−Λh−μhNh0μhe−μht⟶Nht≤Λhμhast⟶∞,(9)Nct≤Λcμc−Λc−μcNc0μce−μct⟶Nct≤Λcμcast⟶∞,(10)Nvt≤Λvμv−Λv−μvNv0μve−μvt⟶Nvt≤Λvμvast⟶∞.
Therefore, from equations (8)–(10), we see that the solutions for human, cattle, and tsetse fly populations enter the following invariant regions:(11)Ωh=Su,Se,Eh,Ih,Rh∈ℝ+5,Su+Se+Eh+Ih+Rh=Nh≤Λhμh,Ωc=Sc,Ec,Ic,Rc,Sc+Ec+Ic+Rc=Nc≤Λcμc,andΩv=Sv,Ev,Iv,Sv+Ev+Iv=Nv≤Λvμv.
The results imply that the region is bounded, well posed, and biologically meaningful as it attracts all solutions inΩ.


### 3.1. Local Stability of Disease-Free Equilibrium

Model (4) has a unique disease-free equilibrium (DFE) which was obtained by setting the right-hand sides of the equations to zero. The DFE is given by(12)E0=Su∗,Se∗,Eh∗,Ih∗,Rh∗,Sc∗,Ec∗,Ic∗,Rc∗,Sv∗,Ev∗,Iv∗=Λhθ+μh,θμhΛhθ+μh,0,0,0,Λcμc,0,0,0,Λvμv+ω,0,0.Using the next generation method as in Van den Driessche and Watmough [14], the associated matrices F for new infection terms and V for the remaining transition terms are evaluated at E0 and given by(13)F=00000bkθB1B2ΛhA1Λh/A1+θΛh/A1μhμh+bkB2ΛhA1Λh/A1+θΛh/A1μh00000000000bdρ0000000bgB2ΛVA8Λh/A1+θΛh/A1μh0bzρΛVμcA8Λc00000000,V=A200000−αhA3000000A500000−αcA6000000A900000−αvA8,where A1=μh+θ, A2=μh+αh, A3=μh+σh+βh, A4=μh+ψh, A5=μc+αc, A6=μc+σc+βc, A7=μc+ψc, A8=μv+ω, A9=μv+ω+αv, B1=1−ϵ, and B2=1−ρ.It follows that the effective or control reproduction number denoted byRe is the spectral radius of the next generation matrix ρFV−1 given as(14)Re=b2αvΛvA1A5A6B22gkΛcαhμhB1θ+μh+A2A3dρ2zαcμcΛhμh+θ2A2A3A5A6A82A9ΛcΛhμh+θ2,=b2αvαcdρ2Λvμczμv+ω2Λcμc+αcA6A9+b2Λvαvαh1−ρ2gkμhμh+1−ϵθμv+ω2Λhμh+αhμh+θA3A9,=Rec+Reh,where Rec represents the effective reproduction number for cattle and Reh is the effective reproduction number for humans. The reproduction number is used to determine whether the disease will die out or persist. That is, the disease persists in the community if Re>1 and it dies out if Re<1. The following theorem summarizes the result.Theorem 2.
The DFE of the model system (4) is locally asymptotically stable if Re<1 and unstable if Re>1.Biologically, it implies that African trypanosomiasis can be eliminated from the community provided that the initial size of the subpopulations of the model (4) is in the basin of attraction of disease-free equilibrium point E0 when Re<1. In other words, the introduction of infectious individuals into a population of susceptible individuals does not induce an epidemic outbreak. On the contrary, if Re>1, it implies that African trypanosomiasis disease will persist in the population.


### 3.2. Global Stability of the Disease-Free Equilibrium Point

Here, we apply Castillo-Chavez et al. [15] approach to analyze the global stability of the disease-free equilibrium point of the model (4).Theorem 3.
The disease-free equilibrium point is globally asymptotically stable ifRe<1 and unstable if Re>1.Proof.
The model (4) can be written in the following format:(15)dYndt=A1Yn−YEo+A2Yi,dYidt=A3Yi,where Yn is a vector of nontransmitting variables, YEo is Yn at disease-free equilibrium point E0, and Yi is a vector consisting of infectious variables.
From model (4), we have Yn = Su,Se,Rh,Sc,Rc,SvT and Yi = Eh,Ih,Ec,Ic,Ev,IvT from which we get(16)Yn−YEo=SuSeRhScRcSv−Λhθ+μhθμhΛhθ+μh0Λcμc0Λvμv+ω=Su−Λhθ+μhSe−θμhΛhθ+μhRhSc−ΛcμcRcSv−Λvμv+ω.
The disease-free equilibrium is globally asymptotically stable if the matrixA1 has real negative eigenvalues and A3 is a Metzler matrix (i.e., the off diagonal elements of A3 are nonnegative, which means A1i,j≥0 for all indices of i≠j). That is,(17)A1=∂Yn∂Su,Se,Rh,Sc,Rc,Sv=−θ+μh0ψh000θ−μh000000−μh+ψh000000−μcψc00000−μc+ψc000000−μv+ω,with eigenvalues λ1=−μh, λ2=−θ+μh, λ3=−μh+ψh, λ4=−μc, λ5=−μc+ψc, and λ6=−μv+ω. Also,(18)A2=∂Yn∂Eh,Ih,Ec,Ic,Ev,Iv=00000−1−ρbkSuNh00000−1−ϵbkSeNh0βh000000000−bdρScNc000βc000−1−ρbgSvNh0−ρbzSvNc0−0,A3=∂Yi∂Eh,Ih,Ec,Ic,Ev,Iv=−P100001−ρbkSu+1−ϵSeNhαh−P2000000−P300bdρScNc00αc−P40001−ρbgSvNh0ρbzSvNc−P500000αv−P6,where P1=μh+αh, P2=μh+σh+βh, P3=μc+αc, P4=μc+σc+βc, P5=μv+ω+αv, and P6=μv+ω. It can be seen that matrix A1 has all eigenvalues which are real and negative and matrix A3 is the Metzler matrix as its off diagonal elements are positive. Thus, the system(19)dYndt=A1Yn−YEo+A2Yiis globally asymptotically stable at disease-free equilibrium point.


## 3.1. Local Stability of Disease-Free Equilibrium

Model (4) has a unique disease-free equilibrium (DFE) which was obtained by setting the right-hand sides of the equations to zero. The DFE is given by(12)E0=Su∗,Se∗,Eh∗,Ih∗,Rh∗,Sc∗,Ec∗,Ic∗,Rc∗,Sv∗,Ev∗,Iv∗=Λhθ+μh,θμhΛhθ+μh,0,0,0,Λcμc,0,0,0,Λvμv+ω,0,0.Using the next generation method as in Van den Driessche and Watmough [14], the associated matrices F for new infection terms and V for the remaining transition terms are evaluated at E0 and given by(13)F=00000bkθB1B2ΛhA1Λh/A1+θΛh/A1μhμh+bkB2ΛhA1Λh/A1+θΛh/A1μh00000000000bdρ0000000bgB2ΛVA8Λh/A1+θΛh/A1μh0bzρΛVμcA8Λc00000000,V=A200000−αhA3000000A500000−αcA6000000A900000−αvA8,where A1=μh+θ, A2=μh+αh, A3=μh+σh+βh, A4=μh+ψh, A5=μc+αc, A6=μc+σc+βc, A7=μc+ψc, A8=μv+ω, A9=μv+ω+αv, B1=1−ϵ, and B2=1−ρ.It follows that the effective or control reproduction number denoted byRe is the spectral radius of the next generation matrix ρFV−1 given as(14)Re=b2αvΛvA1A5A6B22gkΛcαhμhB1θ+μh+A2A3dρ2zαcμcΛhμh+θ2A2A3A5A6A82A9ΛcΛhμh+θ2,=b2αvαcdρ2Λvμczμv+ω2Λcμc+αcA6A9+b2Λvαvαh1−ρ2gkμhμh+1−ϵθμv+ω2Λhμh+αhμh+θA3A9,=Rec+Reh,where Rec represents the effective reproduction number for cattle and Reh is the effective reproduction number for humans. The reproduction number is used to determine whether the disease will die out or persist. That is, the disease persists in the community if Re>1 and it dies out if Re<1. The following theorem summarizes the result.Theorem 2.
The DFE of the model system (4) is locally asymptotically stable if Re<1 and unstable if Re>1.Biologically, it implies that African trypanosomiasis can be eliminated from the community provided that the initial size of the subpopulations of the model (4) is in the basin of attraction of disease-free equilibrium point E0 when Re<1. In other words, the introduction of infectious individuals into a population of susceptible individuals does not induce an epidemic outbreak. On the contrary, if Re>1, it implies that African trypanosomiasis disease will persist in the population.


## 3.2. Global Stability of the Disease-Free Equilibrium Point

Here, we apply Castillo-Chavez et al. [15] approach to analyze the global stability of the disease-free equilibrium point of the model (4).Theorem 3.
The disease-free equilibrium point is globally asymptotically stable ifRe<1 and unstable if Re>1.Proof.
The model (4) can be written in the following format:(15)dYndt=A1Yn−YEo+A2Yi,dYidt=A3Yi,where Yn is a vector of nontransmitting variables, YEo is Yn at disease-free equilibrium point E0, and Yi is a vector consisting of infectious variables.
From model (4), we have Yn = Su,Se,Rh,Sc,Rc,SvT and Yi = Eh,Ih,Ec,Ic,Ev,IvT from which we get(16)Yn−YEo=SuSeRhScRcSv−Λhθ+μhθμhΛhθ+μh0Λcμc0Λvμv+ω=Su−Λhθ+μhSe−θμhΛhθ+μhRhSc−ΛcμcRcSv−Λvμv+ω.
The disease-free equilibrium is globally asymptotically stable if the matrixA1 has real negative eigenvalues and A3 is a Metzler matrix (i.e., the off diagonal elements of A3 are nonnegative, which means A1i,j≥0 for all indices of i≠j). That is,(17)A1=∂Yn∂Su,Se,Rh,Sc,Rc,Sv=−θ+μh0ψh000θ−μh000000−μh+ψh000000−μcψc00000−μc+ψc000000−μv+ω,with eigenvalues λ1=−μh, λ2=−θ+μh, λ3=−μh+ψh, λ4=−μc, λ5=−μc+ψc, and λ6=−μv+ω. Also,(18)A2=∂Yn∂Eh,Ih,Ec,Ic,Ev,Iv=00000−1−ρbkSuNh00000−1−ϵbkSeNh0βh000000000−bdρScNc000βc000−1−ρbgSvNh0−ρbzSvNc0−0,A3=∂Yi∂Eh,Ih,Ec,Ic,Ev,Iv=−P100001−ρbkSu+1−ϵSeNhαh−P2000000−P300bdρScNc00αc−P40001−ρbgSvNh0ρbzSvNc−P500000αv−P6,where P1=μh+αh, P2=μh+σh+βh, P3=μc+αc, P4=μc+σc+βc, P5=μv+ω+αv, and P6=μv+ω. It can be seen that matrix A1 has all eigenvalues which are real and negative and matrix A3 is the Metzler matrix as its off diagonal elements are positive. Thus, the system(19)dYndt=A1Yn−YEo+A2Yiis globally asymptotically stable at disease-free equilibrium point.


## 4. Sensitivity Analysis

The sensitivity analysis is carried out to determine the parameters that have a higher impact on the effective reproduction number. To reduce disease transmission, the parameters that have a higher impact on the effective reproduction number should be targeted for control purposes. Analytically, the sensitivity index ofRe is calculated by using the normalized forward sensitivity index defined as ΥpRe = ∂Re/∂p×p/Re, where p stands for any parameter in effective reproduction number Re [16]. For example, the sensitivity index of Re corresponding to the parameter αv is given as ΥαvRe = ∂Re/∂αv×αv/Re=+0.3182. Other indices are calculated using similar approach, and the results are summarized in Table 3 and in Figure 2.Table 3
Sensitivity indices ofRe using parameter values in Table 4.
ParameterIndex valueParameterIndex valueΛh−0.0036βc−0.2775θ−0.0638αv+0.3182μh−0.0042μv−0.5649ϵ−0.034ω−0.7532αh+0.000192Λv+0.5βh−0.025b+1σh−0.0011k+0.0036Λc−0.4964g+0.0036μc−0.4748z+0.4964σc−0.2006d+0.4964αc+0.0033ρ+0.9761Table 4
Epidemiological data for model (4).
ParameterValues (days)SourceΛh27.5Kajunguri et al. [8]θ0.0002Misra et al. [17]μh0.000046Assumedϵ0.6Assumedψh0.02Gervas et al. [18]αh0.083Ndondo et al. [10]βh0.009Moore et al. [4]σh0.004Gervas et al. [18]k0.62Ndondo et al. [10]Λc22Otieno et al. [7]μc0.00055Kajunguri et al. [8]σc0.006Otieno et al. [7]αc0.083Ndondo et al. [10]βc0.0083Moore et al. [4]ψc0.013Ndondo et al. [10]d0.62Ndondo et al. [10]ρ0.7Ndondo et al. [10]μv0.03Kajunguri et al. [8]ω0.04AssumedΛv120Assumedb0.33Meisner et al. [9]g0.01Ndondo et al. [10]z0.1Assumedαv0.04Meisner et al. [9]Figure 2
Sensitivity analysis ofRe using results in Table 3.


### 4.1. Interpretation of the Sensitivity Indices

From Table3, the parameters αh,αc,Λv,b,k,g,z,d,ρ, and αv have positive indices, indicating that increasing one of these parameters while keeping others constant increases the effective reproduction number, hence increasing the possibility of the disease outbreak. On the contrary, the parameters μh,μc,Λh,θ,ϵ,βh,σh,Λc,σc,βc,μv, and ω have negative indices, implying that increasing one of these parameters and keeping others constant decreases the effective reproduction number, hence reducing the disease burden among human, cattle, and vector population. From these results, the most sensitive parameters are the tsetse fly biting rate b, the proportion of tsetse fly biting on cattle ρ, the rate at which the tsetse flies die due to trapping ω, and tsetse fly natural mortality rate μv, followed by the tsetse fly recruitment rate Λv.Therefore, increasing tsetse fly death rate and reducing tsetse fly biting rate through public health education campaigns on the importance of wearing long-sleeved clothes, clearing the bushes, and the use of repellent solutions to avoid vector-host contact rate would have a higher positive impact in controlling trypanosomiasis transmission in a community.


## 4.1. Interpretation of the Sensitivity Indices

From Table3, the parameters αh,αc,Λv,b,k,g,z,d,ρ, and αv have positive indices, indicating that increasing one of these parameters while keeping others constant increases the effective reproduction number, hence increasing the possibility of the disease outbreak. On the contrary, the parameters μh,μc,Λh,θ,ϵ,βh,σh,Λc,σc,βc,μv, and ω have negative indices, implying that increasing one of these parameters and keeping others constant decreases the effective reproduction number, hence reducing the disease burden among human, cattle, and vector population. From these results, the most sensitive parameters are the tsetse fly biting rate b, the proportion of tsetse fly biting on cattle ρ, the rate at which the tsetse flies die due to trapping ω, and tsetse fly natural mortality rate μv, followed by the tsetse fly recruitment rate Λv.Therefore, increasing tsetse fly death rate and reducing tsetse fly biting rate through public health education campaigns on the importance of wearing long-sleeved clothes, clearing the bushes, and the use of repellent solutions to avoid vector-host contact rate would have a higher positive impact in controlling trypanosomiasis transmission in a community.


## 5. Numerical Simulation

In this section, we simulate model (4) using parameter values shown in Table 4. The Matlab ODE45 solver is used to simulate the model system (4). The initial conditions of the state variables are given as follows: Su=200, Se=180, Eh=150, Ih=100, Rh=80, Sc=200, Ec=170, Ic=140, Rc=100, Sv=3000, Ev=800, and Iv=500. The initial conditions of the state variables are arbitrarily chosen to illustrate specific behaviour of the model (4).


### 5.1. Effects of Interventions on Infected Human and Cattle Population

Figure3 illustrates the impact of different combinations of interventions (human treatment, public health education, and tsetse fly trapping) on the dynamics of African trypanosomiasis in the human population. Combining all three interventions tends to diminish disease transmission in a community faster than treating infected humans only. It is also observed from Figure 3 that treating the infected humans while increasing people’s awareness about the disease has a greater impact than using treatment alone. This implies that apart from using other control measures like public health education and human treatment, there must be an effort to eliminate the tsetse fly vector to eradicate the disease from a community.Figure 3
Impact of various interventions on infected human population whenθ=0.5.Figure4 shows that treatment alone cannot control the disease in the cattle populations. Therefore, to control the disease, the tsetse fly traps must be used as a complementary intervention to cattle treatment.Figure 4
Impact of interventions on infected cattle population.


### 5.2. Effect of Varying Some Parameter Values

Figure5 illustrates the effects of varying human and cattle treatment rates while fixing rates of tsetse fly trapping and public health education constant. As the rate of treating the infected human and cattle increases, the number of infected humans and cattle reduced with time, as shown in Figures 5(a) and 5(b). It indicates that treating the infected humans and cattle has a significant impact on controlling the African trypanosomiasis disease as it reduces the number of sick individuals in the population. Since tsetse flies depend only on blood for their survival, reducing the number of infected humans and cattle by treating them reduces the probability of tsetse flies biting the infected human/cattle, hence reducing the disease’s spread.Figure 5
Variation in treatment rate on infected human and cattle population.
(a)(b)Figure6 shows the effect of varying tsetse fly trapping rates on infected human and cattle populations while keeping human/cattle treatment and public health education rates constant. It is clear from Figures 6(a) and 6(b) that increasing the trapping tsetse flies’ rates has a great contribution toward African trypanosomiasis elimination. It reduces the number of infected humans and cattle in the population. The reason behind this is that tsetse fly trapping reduces the density of tsetse flies in the population. As the number of vectors is reduced in the population, the contact rate of vector to human/cattle is also reduced.Figure 6
Variation in trapping rate on infected human and cattle population.
(a)(b)Figure7 shows the effect of varying public health education rates while constantly fixing trapping and human/cattle treatment rates. As more individuals become aware of the disease and use protective measures such as insect repellent solution, wearing long-sleeved clothes, and clearing their environment, the number of infected humans is diminished in the population with time. This scenario shows that if public health education campaigns are properly conducted, particularly in endemic areas, there is a plausibility of reducing disease transmission, as it reduces the vector human-contact rate.Figure 7
Variation in public health education on infected human population.From Figure8(a), we observe that as the rates of public health education increase, the number of susceptible uneducated individuals decreases as well. This implies that more people have become aware of the disease and join the susceptible educated population. On the contrary, as public health education rates increase, the number of susceptible educated individuals increases, as shown in Figure 8(b). This indicates that more people have become aware of the disease and avoid contact with the tsetse fly by wearing long-sleeved clothes, using insect repellent, and clearing their surroundings.Figure 8
Variation in public health education rate on susceptible human population.
(a)(b)Figure9 indicates the effect of varying the efficacy of public health education on the susceptible educated human population. It is clear that as the efficacy of public health education increases, the number of susceptible educated human population increases with time. It means that increasing the public health education efficacy rates reduces the number of people contacting tsetse flies and hence reducing the number of people joining the exposed class.Figure 9
Variation in efficacy of public health education on susceptible human population.


### 5.3. Effects of Interventions on Effective Reproduction Number

When only one intervention is used to control African trypanosomiasis in a population, it is observed that the disease will not be cleared out since the effective reproduction number is greater than the unit, as shown in Figure10. The tsetse fly trapping seems to be the best control as its effective reproduction number is lower than treatment and public health education. This observation indicates that, for the African trypanosomiasis to be eliminated in the population, a combination of interventions should be considered rather than the application of a single control method.Figure 10
Effects of single intervention on effective reproduction number with respect to vector biting rate.From Figure11, we see that the effective reproduction number increases as the vector biting rates increases. When only two interventions are used, treating the infectious human/cattle and using tsetse fly trapping seem to be effective at reducing the threshold, Re compared to a combination of public health education and tsetse fly trapping, as well as treatment and public health education. The effective reproduction number seems to be very high when only a combination of treatment and public health education is in place. This is because treatment is only applied to sick human/cattle, preventing new cases from occurring. Public health education helps prevent new cases from occurring in a human population, but it does not prevent new cases from occurring in the cattle populations. It can also be observed from Figure 11 that a combination of all three interventions (treatment, public health education, and tsetse fly trapping) is the best strategy for eliminating African trypanosomiasis as it adequately reduces the effective reproduction number than when only two interventions are used.Figure 11
Effects of combining interventions on effective reproduction number with respect to vector biting rate.


## 5.1. Effects of Interventions on Infected Human and Cattle Population

Figure3 illustrates the impact of different combinations of interventions (human treatment, public health education, and tsetse fly trapping) on the dynamics of African trypanosomiasis in the human population. Combining all three interventions tends to diminish disease transmission in a community faster than treating infected humans only. It is also observed from Figure 3 that treating the infected humans while increasing people’s awareness about the disease has a greater impact than using treatment alone. This implies that apart from using other control measures like public health education and human treatment, there must be an effort to eliminate the tsetse fly vector to eradicate the disease from a community.Figure 3
Impact of various interventions on infected human population whenθ=0.5.Figure4 shows that treatment alone cannot control the disease in the cattle populations. Therefore, to control the disease, the tsetse fly traps must be used as a complementary intervention to cattle treatment.Figure 4
Impact of interventions on infected cattle population.


## 5.2. Effect of Varying Some Parameter Values

Figure5 illustrates the effects of varying human and cattle treatment rates while fixing rates of tsetse fly trapping and public health education constant. As the rate of treating the infected human and cattle increases, the number of infected humans and cattle reduced with time, as shown in Figures 5(a) and 5(b). It indicates that treating the infected humans and cattle has a significant impact on controlling the African trypanosomiasis disease as it reduces the number of sick individuals in the population. Since tsetse flies depend only on blood for their survival, reducing the number of infected humans and cattle by treating them reduces the probability of tsetse flies biting the infected human/cattle, hence reducing the disease’s spread.Figure 5
Variation in treatment rate on infected human and cattle population.
(a)(b)Figure6 shows the effect of varying tsetse fly trapping rates on infected human and cattle populations while keeping human/cattle treatment and public health education rates constant. It is clear from Figures 6(a) and 6(b) that increasing the trapping tsetse flies’ rates has a great contribution toward African trypanosomiasis elimination. It reduces the number of infected humans and cattle in the population. The reason behind this is that tsetse fly trapping reduces the density of tsetse flies in the population. As the number of vectors is reduced in the population, the contact rate of vector to human/cattle is also reduced.Figure 6
Variation in trapping rate on infected human and cattle population.
(a)(b)Figure7 shows the effect of varying public health education rates while constantly fixing trapping and human/cattle treatment rates. As more individuals become aware of the disease and use protective measures such as insect repellent solution, wearing long-sleeved clothes, and clearing their environment, the number of infected humans is diminished in the population with time. This scenario shows that if public health education campaigns are properly conducted, particularly in endemic areas, there is a plausibility of reducing disease transmission, as it reduces the vector human-contact rate.Figure 7
Variation in public health education on infected human population.From Figure8(a), we observe that as the rates of public health education increase, the number of susceptible uneducated individuals decreases as well. This implies that more people have become aware of the disease and join the susceptible educated population. On the contrary, as public health education rates increase, the number of susceptible educated individuals increases, as shown in Figure 8(b). This indicates that more people have become aware of the disease and avoid contact with the tsetse fly by wearing long-sleeved clothes, using insect repellent, and clearing their surroundings.Figure 8
Variation in public health education rate on susceptible human population.
(a)(b)Figure9 indicates the effect of varying the efficacy of public health education on the susceptible educated human population. It is clear that as the efficacy of public health education increases, the number of susceptible educated human population increases with time. It means that increasing the public health education efficacy rates reduces the number of people contacting tsetse flies and hence reducing the number of people joining the exposed class.Figure 9
Variation in efficacy of public health education on susceptible human population.


## 5.3. Effects of Interventions on Effective Reproduction Number

When only one intervention is used to control African trypanosomiasis in a population, it is observed that the disease will not be cleared out since the effective reproduction number is greater than the unit, as shown in Figure10. The tsetse fly trapping seems to be the best control as its effective reproduction number is lower than treatment and public health education. This observation indicates that, for the African trypanosomiasis to be eliminated in the population, a combination of interventions should be considered rather than the application of a single control method.Figure 10
Effects of single intervention on effective reproduction number with respect to vector biting rate.From Figure11, we see that the effective reproduction number increases as the vector biting rates increases. When only two interventions are used, treating the infectious human/cattle and using tsetse fly trapping seem to be effective at reducing the threshold, Re compared to a combination of public health education and tsetse fly trapping, as well as treatment and public health education. The effective reproduction number seems to be very high when only a combination of treatment and public health education is in place. This is because treatment is only applied to sick human/cattle, preventing new cases from occurring. Public health education helps prevent new cases from occurring in a human population, but it does not prevent new cases from occurring in the cattle populations. It can also be observed from Figure 11 that a combination of all three interventions (treatment, public health education, and tsetse fly trapping) is the best strategy for eliminating African trypanosomiasis as it adequately reduces the effective reproduction number than when only two interventions are used.Figure 11
Effects of combining interventions on effective reproduction number with respect to vector biting rate.


## 6. Conclusion

In this study, we formulated and analyzed the African trypanosomiasis model with interventions. The model consists of three interventions, namely, public health education for humans, trapping for tsetse flies, and treatment for humans and cattle. Human, cattle, and vector populations were subdivided into different classes concerning their disease status. The invariant set was derived, and the model’s solution was found to be biologically and mathematically meaningful by using the theory of differential equations. We computed the thresholdRe and used it to discuss the local and global stability of the equilibria points. The disease-free equilibrium point was established, and by using the effective reproduction number, its stability was also investigated. The disease-free equilibrium point was locally asymptotically stable when the reproduction number is less than one and unstable when the reproduction number is greater than one. By applying the Metzler stability theory, the disease-free equilibrium point was globally asymptotically stable when the reproduction number is less than one. The sensitivity analysis shows that the control measures based on public health education, tsetse fly trapping, and human and cattle treatment have negative values because increasing them reduces disease transmission in a community. The numerical simulations showed that the combination of all three interventions (treatment, public health education, and trapping) considerably cleared out the population’s disease burden faster than using only two interventions. Furthermore, the numerical results showed that, with an increase in public health education rates against African trypanosomiasis disease, the number of susceptible educated individuals increases gradually. We recommend that, to keep the disease under control, public health education campaigns through seminars and media like television, radio, magazines, and mobile networks should be spread, especially in the rural areas, to make people aware of the disease. We also recommend that the community, especially those living in endemic areas, should be encouraged to use tsetse traps as they are cheap. They have shown a great impact on tsetse control and trypanosomiasis elimination. The model presented in this study is not exhaustive. Therefore, the assumptions made during model formulation can be relaxed to incorporate the aspect of infected immigrants, age structure, climatic change, and the cost-effectiveness of the control strategies.



---

*Source: 3690472-2020-11-26.xml*