_{1}

Improvements in sanitation and the provision of clean drinking water led to the elimination of typhoid fever from developed countries in the beginning of the 20th century. However, Salmonella typhi and paratyphi remain a major source of morbidity and mortality in many developing countries today. The dynamics of typhoid transmission are poorly understood. In this study, we develop a novel mathematical model that captures the role of both human to human interaction and human to environment interaction in the transmission dynamics of typhoid fever. Our results have shown the feasible impact of different methods of typhoid control, including vaccination, improved treatment strategies, and investment in clean water and sanitation.

Typhoid fever is a communicable disease, found only in man and occurs due to systemic infection mainly by Salmonella typhi organism [

Our main objective in the present paper is to develop an SIIcR-B (susceptible, symptomatic infectious, asymptomatic infectious, recovered, bacteria concentration) model of typhoid fever with vaccination, treatment and water sanitation as control strategies that has not been investigated in prior studies. Our major assumption is existence of both direct transmission of typhoid from infected individuals to susceptible and indirect transmission of bacteria from the environment to the susceptible, the other assumptions; all susceptible individuals are equally likely to be infected by infectious individuals in case of contact, and we also assume direct transmission of typhoid from infected to susceptible individuals and that there is a constant recruitment rate to the susceptible population. Furthermore, we assume that the rate of transmission for carriers is greater than that of symptomatic infectious individuals.

We first develop a more realistic model for typhoid. The model subdivides the human population of interest into four compartments: susceptible humans (S), infected humans (I), carrier humans (I_{c}), and recovered humans (R). Previous models of typhoid dynamics [

Individuals in the class, I, can recover from typhoid at the rate, η 1 . The carrier individuals can also either progress to carrier class, I c , at rate α or recover from typhoid, but with a significantly slow rate, η 2 . Infected individuals in both infectious state and carrier state excrete bacteria into the environment. However, the rate of excretion by the infectious group, ϵ 1 , is significantly higher than that by the carrier group, ϵ 2 . Note that despite low excretion of bacteria by the carrier group, because of its extremely long duration without showing any sickness the carrier group plays an important role on infection dynamics of typhoid. Growth curves of organisms are often described well with the 2 logistic models [

From the assumptions, descriptions and the compartmental diagram in

d S d t = Λ + ϕ R − ( λ + μ h + θ ) S (1)

d I d t = p λ S + α I c − ( μ h + d 1 + η 1 + ϵ 1 ) I (2)

d I c d t = ( 1 − p ) λ S − ( α + μ h + d 2 + η 2 + ϵ 2 ) I c (3)

d B d t = r ( 1 − B K ) B + ϵ 1 I + ϵ 2 I c − μ b B (4)

d R d t = η 1 I + η 2 I c + θ S − ( μ h + ϕ ) R (5)

d N d t = Λ − μ h N − d 1 I − d 2 I c (6)

where

λ = β I + γ I c + τ B K + B

Equation (1) describes the dynamics of susceptible in the community of size N. The death rate of the susceptible individuals is represented by, is the rate of recruitment into susceptible class, is the rate of exposure to contaminated food and water, is the probability of susceptible catching typhoid fever, is the susceptible and is the concentration of Salmonella typhi bacteria in food and water. Equations (2) and (3) describes the dynamics of infected people in the community, their number increases as susceptible become infected and decreases as the infected recovers or die from the disease or natural death. Measures to limit the spread of the disease, such as hygiene and total sanitation reduce the amount of Salmonella typhi bacteria in the environment. Equation (4) describes the dynamics of pathogenic Salmonella typhi bacteria in environment, comprising the contaminated food or water consumed by people and unhygienic handling of typhoid fever patients and their waste products. Equation (5) describes the dynamics of effect of treatment or lack of treatment to the population of infected people.

We show that if the system starts with non-negative initial conditions

( S 0 , I 0 , I c 0 , R 0 , B 0 ) , the solutions/trajectories of (1)-(5) will remain non-negative for all t ∈ ( 0, ∞ ) . This is an ideal condition to check since the model monitors human population and the pathogen concentration in the aquatic environment. We thus have the following theorem: Theorem 1. Given that the initial conditions of the system (1)-(5) ( S 0 > 0 , I 0 > 0 , I c 0 > 0 , R 0 > 0 , B 0 > 0 ) , the resulting solutions ( S ( t ) , I ( t ) , I c ( t ) , R ( t ) , B ( t ) ) are all non-negative for all t ∈ ( 0, ∞ ) .

Proof. To show positivity of solution, it is enough to show that each of the tra-jectories of system (1)-(5) is non negative for all t > 0 . From Equation (1), the differential inequality describing the evolution of the susceptible population over time is given by

d S d t ≥ Λ − ( λ + μ h + θ ) S (7)

The resulting differential inequality can be solved by separation of variables. Since at t = 0 , S ( 0 ) = S 0 , then the complete solution to the differential inequality for the susceptible population is given by

S ( t ) ≥ Λ μ h + θ + ( S 0 − Λ μ h + θ ) e − ( μ h + θ ) t (8)

lim n → ∞ inf S ( t ) ≥ 0

Using the same principle, the rest of the phase space variables as t approaches infinity can be shown to satisfy

I ( t ) ≥ I 0 e − ( μ h + d 1 + η 1 + ϵ 1 ) t ,

I c ( t ) ≥ I 0 e − ( μ h + d 2 + η 2 + ϵ 2 ) t

R ( t ) ≥ R 0 e − ( μ h + ϕ ) t

From which the limit inf of the corresponding state variables can be shown to be non-negative. Using the equation describing the evolution of the pathogen concentration, we have a differential inequality given by

d B d t + ( μ b − r ) B ≥ − B 2 K (9)

Equation (9) is a Bernoulli type of equation. It is solved by substitution i.e. B = y − 1 to obtain

B ( t ) ≥ − e − ( μ b − r ) t A K ( μ b − r ) (10)

The model can be separated into two parts which include, the human population T H and the concentration of the pathogen in the aquatic environment T B such that T H = { ( S ( t ) , I ( t ) , I c ( t ) , R ( t ) ) ∈ ℝ + 4 : S + I + I c + R = N } and

T B = { B ( t ) ∈ ℝ + 1 } respectively. From Equation (1) the differential inequality of the susceptible population is given by

d S d t = Λ + ϕ R − ( λ + μ h + θ ) S

d S d t + ( μ h + θ ) S ≤ Λ + ϕ R (11)

Using a suitable integrating factor, I . F = e − μ t , the differential inequality (11) can be solved to obtain

S ( t ) ≤ Λ μ h + e − μ t ∫ 0 t ω R ( x ) e − μ x d x (12)

Following the theorem of differential inequality by Birkhoff and Rota [

lim t → ∞ Sup S ( t ) ≤ Λ μ h

Therefore, the state variable describing the evolution of the susceptible population is less or equal to the ratio of the recruitment rate and the natural mortality rate. We note also that the total population is given as

N = S + I + I c + R . If we take the time derivative of N i.e. d N d t and substitute

the Equations (1)-(4) into the resulting expression we obtain

d N d t = Λ − μ h N − d 1 I − d 2 I c (13)

The solution (12) can be obtained by separating variables and integrating both sides with respect to the corresponding variable. This result into

ln | Λ − μ h N | ≥ − μ h t + c

where c is a constant of integration. If we exponentiate both sides and assume that the initial total population is N 0 , the solution becomes

N ≤ Λ μ h − ( Λ μ h − N 0 ) e − μ h t (14)

therefore,

lim t → ∞ Sup N ( t ) ≤ Λ μ h

Since N is the sum of all state space variables, then each of the individual

state variables is less or equal to Λ μ h . Using Equation (5), we assume that the

growth rate of the pathogen in linear at a constant rate r . We therefore obtain a differential inequality

d B d t + ( μ b − r ) B ≤ ϵ 1 I + ϵ 2 I c (15)

d B d t + ( μ b − r ) B ≤ ( ϵ 1 + ϵ 2 ) Λ μ h (16)

The solution to this equation can be obtained by using a suitable integrating factor to obtain

B ( t ) ≤ ( ϵ 1 + ϵ 2 ) Λ μ h ( μ b − r ) + A e ( μ b − r ) t (17)

therefore,

lim t → ∞ Sup B ( t ) ≤ ( ϵ 1 + ϵ 2 ) Λ μ h ( μ b − r )

The domain of biological significance of the system (1)-(5) is

T = [ S , I , I c , R , B ( t ) ≥ 0 : S + I + I c + R ≤ Λ μ h , B ( t ) ≤ ( ϵ 1 + ϵ 2 ) Λ μ h ( μ b − r ) ] . (18)

The domain T is positively invariant under the flow induced by the system (1)-(5). Therefore, the system (1)-(5) is biologically meaningful and it is feasible to analyze the model in the domain T .

The basic reproduction number denoted by R 0 is the average number of secondary infections caused by an infectious individual during his or her entire period of infectiousness Diekmann et al. [

F V − 1 = [ ∂ F i ( E 0 ) ∂ x j ] [ ∂ V i ( E 0 ) ∂ x j ] − 1 (19)

where F i is the rate of appearance of new infection in compartment i , V i is the transfer of infections from one compartment i to another and E 0 is the disease-free equilibrium. From system equation of the system (1)-(4), we re-write the equations with infectious classes, I , I c and B . This leads to the system

d I d t = p λ S + α I c − ( μ h + d 1 + η 1 + ϵ 1 ) I

d I c d t = ( 1 − p ) λ S − ( α + μ h + d 2 + η 2 + ϵ 2 ) I c

d B d t = r ( 1 − B K ) B + ϵ 1 I + ϵ 2 I c − μ b B

Jacobian at diseases free point ( E 0 ) is computed and found to be

J ( E 0 ) = [ p τ B S K + B + p β I S + p γ I c S ( 1 − p ) τ B S K + B + ( 1 − p ) β I S + ( 1 − p ) γ I c S 0 ] (20)

from which we obtain:

F i = [ p τ B S K + B + p β I S + p γ I c S ( 1 − p ) τ B S K + B + ( 1 − p ) β I S + ( 1 − p ) γ I c S 0 ] (21)

V i = [ ( μ h + d 1 + η 1 + ϵ 1 ) I ( μ h + d 2 + η 2 + ϵ 2 ) I c μ b B − ϵ 1 I ϵ 2 I c + r ( 1 − B K ) B ] (22)

Partial differentiation of F i and V i with respect to I , I c and B gives

F = [ p β S 0 p γ S 0 p τ S 0 K ( 1 − p ) β S 0 ( 1 − p ) β S 0 ( 1 − p ) τ S 0 K 0 0 0 ] (23)

V = [ ( μ h + d 1 + η 1 + ϵ 1 ) 0 0 0 ( μ h + d 2 + η 2 + ϵ 2 ) 0 − ϵ 1 − ϵ 2 μ b − r ] (24)

V − 1 was computed and found to be:

V − 1 = 1 a 1 a 2 ( μ b − r ) [ a 2 ( μ b − r ) α ( μ b − r ) 0 0 a 1 ( μ b − r ) 0 ϵ 1 a 2 a ϵ 2 a 1 ]

Lastly F V − 1 was calculated and results is:

The eigenvalues of F V − 1 was calculated as follows:

let

Z 1 = p β S 0 a 2 ( u b − r ) + p τ S 0 ϵ 1 a 2 K

Z 2 = p β S 0 α ( u b − r ) + p γ a 1 S 0 ( u b − r ) + p τ S 0 f K

Z 3 = p τ S 0 K

Z 4 = ( 1 − p ) β S 0 a 2 ( u b − r ) + ( 1 − p ) τ S 0 ϵ 1 a 2 K

Z 5 = ( 1 − p ) γ a 1 S 0 ( u b − r ) + ( 1 − p ) a 1 τ S 0 f K

Z 6 = ( 1 − p ) a 1 τ S 0 K

We get

| λ − Z 1 Z 2 Z 3 Z 4 λ − Z 5 Z 6 0 0 λ | = 0

λ [ ( λ − Z 5 ) ( λ − Z 1 ) − Z 2 Z 5 ] = 0

λ ( λ 2 − ( Z 1 + Z 5 ) λ + ( Z 1 Z 5 − Z 2 Z 4 ) ) = 0

it can be shown that Z 1 Z 5 − Z 2 Z 4 = 0

λ 2 ( λ − ( Z 1 + Z 5 ) ) = 0 (25)

λ = 0 , 0 , Z 1 + Z 5 (26)

thus the maximum eigenvalues is

Z 1 + Z 5

which gives

R e = Z 1 + Z 5

therefore

R e = p β S 0 a 2 ( μ b − r ) + a 1 ( 1 − p ) γ S 0 ( μ b − r ) + p τ S 0 ϵ 1 a 2 K + ( 1 − p ) τ S 0 a 1 ϵ 2 K (27)

= S 0 a 1 a 2 ( μ b − r ) ( ( μ b − r ) [ p a 2 β + a 1 ( 1 − p ) γ ] + τ K [ a 2 p ϵ 1 + ( 1 − p ) a 1 ϵ 2 ] ) (28)

= S 0 p a 1 ( μ b − r ) τ K ( β ( μ b − r ) + τ ϵ 1 K ) + ( 1 − p ) S 0 a 2 ( μ b − r ) (29)

= ( S 0 p a 1 β + S 0 p τ a 1 ( μ b − r ) K ϵ 1 ) + ( ( 1 − p ) S 0 a 2 γ + τ K ( 1 − p ) S 0 γ a 2 ( μ b − r ) ϵ 2 ) (30)

This can be written as

R e = ( R 01 + R 03 ) + ( R 02 + R 04 )

where

R 01 = ( S 0 p a 1 ) (31)

R 02 = S 0 ( 1 − p ) γ a 2 (32)

R 03 = ( S 0 p τ a 1 ( μ b − r ) K ) ϵ 1 (33)

R 04 = S 0 ( 1 − p ) τ a 2 ( μ b − r ) K ϵ 2 (34)

S 0 = ( μ h + ϕ ) Λ μ h ( μ h + θ + ϕ ) (35)

a 1 = μ h + d 1 + η 1 + ϵ 1 (36)

a 2 = μ h + d 2 + η 2 + ϵ 2 (37)

Thus R_{e} is the effective reproduction number (basic reproduction number

with controls). The terms 1 a 1 and 1 a 2 indicate the maximum time an individual

is expected to stay in compartments I and I c respectively. The reproduction number consists of four terms which characterize the contribution from the different pathways to new infections with typhoid.

In this study our objective is to model transmission dynamics of typhoid fever via the direct and indirect paths, we want to analyze what happens to the system when control measures like vaccination, treatment and water sanitation either effected or not. Furthermore we want to know the role of carriers and symptomatic individuals to this dynamical system. Various graphical representations have been generated with the help of MATLAB which will support our analytical results. Since, most of the parameters were not readily available; it was found convenient to pick from other sources and unavailable data were estimated. In order to perform simulations, baseline values of parameters from

It can be seen from Figures 2(a)-(d) that sanitation, vaccination, treatment of both symptomatic and asymptomatically infected individuals respectively, are observed to reduce the severity of the disease if such parameters are increased.

Parameter | Value | Description | Source |
---|---|---|---|

10^{6} | Constant human recruitment rate | [ | |

0.4/year | Mortality rate for bacteria, including phage degradation | Estimated | |

0.167/year | Natural human mortality rate | Estimated | |

0.15/year | Disease induced death rate | [ | |

0.6/year | Disease induced death rate | Estimated | |

0.02/year | Effective contact rate between individuals (contact sufficient) | Estimated | |

0.7/year | Per capita contact rate for humans and contaminated water | Estimated | |

0.04/year | Recovery rate of infectious humans | [ | |

0.05/year | Recovery rate of infectious humans | Estimated | |

0.5/year | Bacteria shed rate into the water supply by infectious human | Estimated | |

0.4/year | Bacteria shed rate into the water supply by infectious human | Estimated | |

0.33/year | Per capita rate at which recovered humans are susceptible | Estimated | |

0.4/year | Per capita rate at which susceptible humans are vaccinated | [ | |

0.01/year | (Maximum) per capita growth rate for S. typhi bacteria | [ | |

100/year | Carrying capacity for S. typhi | Estimated | |

0.8/year | Proportion of infected individuals who are symptomatic | Estimated |

Recovery of infected individuals has a twofold benefit in the fight against the infection; 1) it leads to reduction the likelihood getting new infections through direct person to person transmission, and 2) a negligible amount of the pathogen would be shed into the aquatic reservoir greatly reducing the infection risk, most especially for the community that may have direct contact with a potentially contaminated water source. In addition, recovering individuals acquire some immunity to the disease which only wanes over a reasonable period of time hence reducing the susceptible. On the other hand, it can be seen from Figures 3(a)-(e) that increased discharge of fresh bacteria into the aquatic environment by either the symptomatic or asymptomatic individuals, increased probability of contracting bacteria from the environment. Likewise, high personal contact rate with either the symptomatic or asymptomatic individuals increase the risk of contracting typhoid fever. The worst case scenario in case of epidemic outbreak may be experienced when the level of hygiene is poor; maximizing person-to- person contact rate and when there is no accesses to clean water; which maximizes disease transmission through contact with the contaminated reservoir. If immunization of the susceptible population alone does not bring about typhoid elimination, then measures to reduce per case or per carrier infectivity, such as improved sanitation or hand washing with soap, might be considered instead

of or in conjunction with vaccination [

In this paper, a deterministic model which incorporates person-to-person contact rate and person-environment was presented and analyzed. Important mathematical features of the models such as the threshold for the epidemic, steady states, positivity and boundedness of solutions as well as the region of biological significance were determined. The model was shown to have a disease free equilibrium which is both locally and globally asymptotically stable when the reproduction number is less than unity. This disease free equilibrium is unstable when the disease threshold is greater than unity. The model has a unique endemic equilibrium for R e > 1 . In general, vaccination coverage, therapeutic treatment and water sanitation have been shown to play a pivotal role of diminishing the outbreak when they are encouraged in the community. On the other hand, it has been found beneficial to minimize contact with typhoid patients, avoid spoiling water sources with feces, and people should use Latrines. We acknowledge the fact that this work may have shortfalls as follows. The model does not take into account education campaigns. However, education is recommended as a preventive measure for typhoid. Proposed improvements of the model include consideration of a combination of hygiene, vaccination, education campaigns and biological control with salmonella specific bacteriophage.

The author would like to thank the editors of this journal for their inputs which has helped to improve this paper.

Edward, S. (2017) A Deterministic Mathematical Model for Direct and Indirect Transmission Dyna- mics of Typhoid Fever. Open Access Li- brary Journal, 4: e3493. https://doi.org/10.4236/oalib.1103493