* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Download MODELING THE EFFECTS OF CARRIERS ON TRANSMISSION
Human cytomegalovirus wikipedia , lookup
West Nile fever wikipedia , lookup
Middle East respiratory syndrome wikipedia , lookup
Neonatal infection wikipedia , lookup
Onchocerciasis wikipedia , lookup
Neglected tropical diseases wikipedia , lookup
Trichinosis wikipedia , lookup
Chagas disease wikipedia , lookup
Sexually transmitted infection wikipedia , lookup
Marburg virus disease wikipedia , lookup
Oesophagostomum wikipedia , lookup
Hospital-acquired infection wikipedia , lookup
African trypanosomiasis wikipedia , lookup
Schistosomiasis wikipedia , lookup
Coccidioidomycosis wikipedia , lookup
Eradication of infectious diseases wikipedia , lookup
Leptospirosis wikipedia , lookup
Infectious mononucleosis wikipedia , lookup
Hepatitis C wikipedia , lookup
MATHEMATICAL BIOSCIENCES AND ENGINEERING Volume 8, Number 3, July 2011 doi:10.3934/mbe.2011.8.711 pp. 711–722 MODELING THE EFFECTS OF CARRIERS ON TRANSMISSION DYNAMICS OF INFECTIOUS DISEASES Darja Kalajdzievska Department of Mathematics University of Manitoba, Winnipeg, Manitoba R3T 2N2 Canada Michael Y. Li Department of Mathematical and Statistical Sciences University of Alberta, Edmonton, Alberta, T6G 2G1 (Communicated by Yang Kuang) Abstract. An S-Ic -I-R epidemic model is investigated for infectious diseases that can be transmitted through carriers, infected individuals who are contagious but do not show any disease symptoms. Mathematical analysis is carried out that completely determines the global dynamics of the model. The impacts of disease carriers on the transmission dynamics are discussed through the basic reproduction number and through numerical simulations. 1. Introduction. For certain infectious diseases, there are individuals who are able to transmit their illness but do not exhibit any symptoms. These individuals are called “carriers” and they play an important role in the transmission of the disease. There are two types of carriers. Genetic carriers carry the illness on their recessive genes. They can only pass on their disease to their children and are not contagious. The focus of our study is on infectious disease carriers. These individuals are asymptomatic and are likely unaware of their conditions, and therefore are more likely to infect others. An infectious disease that produces long-term asymptotic carriers is the Typhoid fever caused by the bacteria Salmonella Typhi. Typhoid fever reached public notoriety at the beginning of the 20th century with the cases of “Mr. N the milker” in England and Typhoid Mary in the US. These individuals infected hundreds of people over the decades while they worked in the food production industry and private homes. Even today, Typhoid fever infects 21 million people and kills 200,000 worldwide every year. Asymptomatic carriers are believed to play an essential role in the evolution and global transmission of Typhi, and their presence greatly hinders the eradication of Typhoid fever using treatment and vaccination [13]. Another major infectious disease that causes long-term asymptomatic carriage is hepatitis B, a liver disease caused by the HBV virus of the Hepadnavirus family. Most people infected with HBV recover completely and develop a lifelong immunity to the virus. However, about 5-10% of adults will develop chronic HBV infection, 2000 Mathematics Subject Classification. Primary: 92D40, 92D25; Secondary: 34D20. Key words and phrases. Disease carriers, epidemic models, global stability, Lyapunov functions. 711 712 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI and 15-25% of these will develop liver disease. Hepatitis B’s symptoms include jaundice, abdominal pain, nausea, fatigue and joint pain. About 30% of people with the disease do not show any of these symptoms. A major public-health challenge in the control of hepatitis B infection in many countries is the existence of a large pool of chronic carriers who are responsible for transmitting most of the new infections. Infections of other pathogens are also know to produce asymptomatic carriers. The Epstein-Barr Virus (EBV) of the herpes family is one of the most common viruses in humans. EBV infection commonly causes infectious mononucleosis, also known as glandular fever. Most people infected with EBV are asymptomatic, as it remains dormant in those who have had it for the rest of their lives in the cells of the throat and the immune system. Clostridium difficile is a bacterium that causes Clostridium difficile-associated diseases (CDAD). CDAD remains the most common cause of acute hospital-acquired diarrhea, responsible for more than 300,000 cases of diarrhea annually in acute-care facilities in the United States. Asymptomatic carriage rates of up to 30% have been reported in long-term care facilities. It is believed that carriers are responsible for transmission and large outbreaks of CDAD in Europe and North America [12]. Despite their public health significance, the effects of carriers on the transmission dynamics of the disease have not received adequate research attention in the mathematical modeling literature. One of the earlier attempts was Kemper [7], in which a general mathematical model that incorporates disease carriers was developed and analyzed. Medley et al. [10] used a mathematical model for hepatitis B with carriers to discussed the effects of HBV vaccination. Several other studies using large-scale computational models with carriers are specifically aimed at hepatitis B and other diseases [15, 1, 2, 14, 11]. In the present paper, we propose a general mathematical model for infectious diseases with asymptomatic carriers to investigate the effects of carriers on the transmission dynamics. We have derived the basic reproduction number R0 and show that the global dynamics of the model are completely determined by the values of R0 . Since R0 explicitly involves parameters related to disease carriage, we are able to discuss the impact of disease carriage on R0 . We have also carried out numerical simulations of the model using parameter values that are pertinent to hepatitis B infection, and investigated the effects of carriers on the HBV transmission dynamics. Mathematically, our proof of the global stability of the unique endemic equilibrium when R0 > 1 nontrivially utilizes the method of global Lyapunov functions. The model derivation is given in the next section. The basic reproduction number is derived and discussed in Section 3. Global stability of the disease-free and endemic equilibrium is proved in Section 4 and 5, respectively. In Section 6, effects of carriers on the transmission of chronic hepatitis B infection are discussed. 2. A general epidemic model with asymptomatic carriers. 2.1. Model formulation. We formulate an S-Ic -I-R compartmental model where S, Ic , I, and R represent the susceptible, carrier, symptomatically infectious or infectious for short, and removed classes, respectively. A susceptible individual can be infected through direct contact with an infectious individual or a carrier. A newly infected individual can become a carrier with probability p, or shows disease symptoms with probability 1 − p. We assume that the rate of transmission β for carriers is higher than the rate γ of symptomatically infected individuals due to the fact that they are more likely to be unaware of their condition, and therefore continue CARRIERS’ EFFECTS ON INFECTIOUS DISEASES 713 with their regular behaviours. Carriers may become symptomatic at a rate α. For infections such as HBV for which carriage can remain life-long, α can be regarded as rate of diagnosis. We assume a constant influx b of susceptibles, and let d1 , d2 , d3 , d4 denote the death rates of those in the susceptible, carrier, infectious, and removed classes, respectively. Here d1 , d4 can be considered as natural death rates, while d2 and d3 may include both natural and disease-related death. We incorporate a simple vaccination strategy in which a fraction θ of the susceptible population is vaccinated and is fully protected by the vaccine. Symptomatically infected individual recover with rate π, and we assume that recovered individuals are permanently immune. Parameters in the model are summarized and explained in Table 1, and the model is depicted in the transfer diagram in Figure 1. We assume that all parameters in the model are nonnegative and that b > 0, di > 0, i = 1, 2, 3, 4. Table 1. Parameters in the Model b: d1 , d4 : d2 , d3 : β: γ: α: π: p: θ: Rate of influx of susceptibles Natural death rates Death rates for Ic and I compartments, respectively, including both natural and disease-caused death Transmission coefficient for the carrier compartment Ic Transmission coefficient for the symptomatically infected compartment I Rate at which carriers develop symptoms Rate of recovery Probability of a newly infected individual is asymptomatic Vaccination rate 6 d2 Ic 6 d1 S b - S 6 Ic * d4 R pS(βI c + γI) HH H Hc H (1 − p)S(βI + γI) HH j αIc ? I πI * d3 I ? θS Figure 1. Transfer diagram of model (1). R 6 714 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI Based on our assumptions and the transfer diagram, we can derive the following system of differential equations that govern our model. S 0 = b − d1 S − S(βIc + γI) − θS Ic0 = pS(βIc + γI) − (d2 + α)Ic I 0 = (1 − p)S(βIc + γI) − (d3 + π)I + αIc (1) R0 = πI + θS − d4 R. We note that disease carriage is different from disease latency in that individuals in the carrier state are infectious while those in the latent period are not. Our model (1) is thus different from the traditional SEIR models that incorporate disease latency. In the special case when β = 0, the Ic class is not infectious and can be considered as latent, and our model becomes a modified SEIR model in which new infections can be either latent or infectious. While both Ic and I are infectious, model (1) is different from the differential infectivity models considered in [5], since new infections from Ic or I may enter either compartment with certain probability. Our model is more general than the carrier model in [7] in that we incorporate demography and disease-caused death. We also allow carriers to become symptomatic over their life time. Our model is different from the carrier model in [10] in that we allow new infections to be either symptomatic or asymptomatic with certain probabilities. 2.2. Feasible region and equilibria. From (1) we have that S 0 ≤ b − (d1 + θ)S, and thus lim supt→∞ S(t) ≤ (d1b+θ) along each solution. Also from (1) we see that ¯ N 0 = b − d1 S − d2 Ic − d3 I − d4 R ≤ b − dN, ¯ The equation for where d¯ = min{d1 , d2 , d3 , d4 }. Therefore lim supt→∞ N (t) ≤ b/d. R can be omitted in our analysis as R does not appear in the other equations. This shows that the model can be studied in the feasible region Γ = { (S, Ic , I) ∈ R3 : S ≤ b/(d1 + θ), S + Ic + I ≤ b/d¯}. + It can be verified that Γ is positively invariant with respect to (1). Once the dynamics of (S, IC , I) are understood, those of R can then be determined from the equation R0 = πI + θS − d4 R. The first step in our analysis is to find equilibria (S ∗ , Ic∗ , I ∗ ) from equations 0 = b − d1 S ∗ − S ∗ (βIc∗ + γI ∗ ) − θS ∗ , 0 = pS ∗ (βIc∗ + γI ∗ ) − (d2 + α)Ic∗ , (2) 0 = (1 − p)S ∗ (βIc∗ + γI ∗ ) − (d3 + π)I ∗ + αIc∗ . Model (1) always has a disease-free equilibrium P0 = ( d1b+θ , 0, 0). An endemic equilibrium P ∗ = (S ∗ , Ic∗ , I ∗ ) satisfies S ∗ , Ic∗ , I ∗ > 0. From the equilibrium equations we can show that a unique P ∗ exists with S∗ = (d3 + π)(d2 + α) . pd3 β + (d2 + α)γ + p(πβ − d2 γ) For P ∗ to exist in the feasible region Γ, it is necessary and sufficient that 0 < S ∗ ≤ b b d1 +θ , or equivalently, (d1 +θ)S ∗ ≥ 1. Define R0 = b b pd3 β + (d2 + α)γ + p(πβ − d2 γ) 1 = . ∗ S d1 + θ d1 + θ (d3 + π)(d2 + α) (3) CARRIERS’ EFFECTS ON INFECTIOUS DISEASES 715 Then R0 is a threshold parameter that determines the number of equilibria. We will show in Section 3 that R0 is the basic reproduction number. Proposition 1. If R0 ≤ 1 then P0 is the only equilibrium in Γ; if R0 > 1, then there are two equilibria, P0 and a unique endemic equilibrium P ∗ . 3. The basic reproduction number. Rewrite R0 in (3) as h 1 α 1 i b 1 +p β + γ . (4) R0 = (1 − p) γ d3 + π d2 + α d2 + α d3 + π d1 + θ In the following, we show that R0 is the basic reproduction number, namely, it represents the average number of secondary infections caused by a single infective in an entirely susceptible population during its entire infectious period. When a single infective is introduced into the population, with probability 1−p it is a non-carrier, hence makes γ effective contacts per unit time. This is multiplied by the average infectious period d31+π for non-carriers; with probability p the infective is a carrier, and hence makes β effective contacts per unit time during the average period d21+α it remains a carrier. This number should be augmented by the number of infections γ d31+π caused by this infective after it becomes a non-carrier, with probability d2α+α to survive the carrier stage. Therefore, the expression in the big square brackets in (4) is the per capita average number of secondary infections. This number multiplied by the number of susceptibles at the disease-free equilibrium, b d1 +θ , gives R0 . The carriers in our system can have a great effect on R0 . The parameters β, α, and p are all related to the carrier class and all appear in the basic reproductive number. It is straightforward from (3) that R0 increases as β increases. This agrees with the intuition that higher transmissibility increases the basic reproduction number. To see the effect of p on R0 we note h ∂R0 γ β γ α i b = − + + ∂p d3 + π d2 + α d3 + π d2 + α d1 + θ i b(d2 + α) h d2 = β− γ , d1 + θ d3 + π and thus ∂R0 ∂p > 0 if d2 γ. (5) d3 + π We see that a greater probability to develop carriage will increase the basic reproduction number under the condition (5). We can also analyze the effect of diagnosis rate α on R0 . Straightforward computation gives h i b β d2 ∂R0 γ =p − + ∂α (d2 + α)2 d3 + π (d2 + α)2 d1 + θ h i d2 bp =− β− γ , 2 (d2 + α) (d1 + θ) d3 + π β> 0 and thus ∂R ∂α > 0 if the same condition (5) holds. From these analysis we see that parameter p and α have opposite effects on R0 : while a higher probability p of carriage increases R0 , a higher diagnosis rate α of 716 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI carriage decreases R0 . The latter aspect can be a very useful control strategy and will be further explored in Section 6 through numerical simulations. Biologically, condition (5) only requires that the transmissibility β of carriers is not too small compared to that of the symptomatically infected. This is likely to hold for many diseases with carriers since carriers can unknowingly infect many people. 4. Stability of the disease-free equilibrium. To examine the local stability of the disease-free equilibrium P0 we evaluate the Jacobian matrix at P0 = ( d1b+θ , 0, 0) −d1 − θ −β( d1b+θ ) −γ( d1b+θ ) b b . 0 pβ( d1 +θ ) − (d2 + α) pγ( d1 +θ ) J(P0 ) = b b 0 (1 − p)β( d1 +θ ) + α (1 − p)γ( d1 +θ ) − (d3 + π) We have the following stability result that shows R0 is a sharp threshold. Proposition 2. P0 is locally asymptotically stable if R0 < 1 and is unstable if R0 > 1. Proof. One eigenvalue of J(P0 ) is λ1 = −(d1 + θ) < 0. The other two eigenvalues λ2 , λ3 are eigenvalues of the 2 × 2 matrix pβ( d1b+θ ) − (d2 + α) pγ( d1b+θ ) A= . (1 − p)β( d1b+θ ) + α (1 − p)γ( d1b+θ ) − (d3 + π) We want to show, when R0 < 1, that the Routh-Hurwitz conditions hold, namely, tr(A) < 0 and det(A) > 0. Simple calculations show that h pβ b i h (1 − p)γ b i d1 +θ d1 +θ tr(A) = (d2 + α) − 1 + (d3 + π) −1 . d2 + α d3 + π Using our assumption that R0 = pβ d1b+θ <1 and −1<0 and d2 + α Therefore pβ d1b+θ pβ b d1 +θ [ (d2 +α) + pαγ (d2 +α)(d3 +π) (1 − p)γ d1b+θ d3 + π + (1−p)γ (d3 +π) ] < 1 we have < 1. (1 − p)γ d1b+θ − 1 < 0. d2 + α d3 + π This shows that tr(A) < 0. Now we calculate h ih i b b det(A) = pβ − (d2 + α) (1 − p)γ − (d3 + π) d1 + θ d1 + θ h i b b − (1 − p)β + α pγ d1 + θ d1 + θ b b = (d2 + α)(d3 + π) − (d3 + π)pβ − d2 (1 − p)γ d1 + θ d1 + θ = (d2 + α)(d3 + π)[1 − R0 ]. Therefore, det(A) > 0 if and only if R0 < 1. This proves the proposition. Theorem 4.1. P0 is globally asymptotically stable in the feasible region Γ if R0 ≤ 1. Proof. To prove the global asymptotic stability of P0 we use the method of Lyapunov functions. Define h β i γα γ L= + Ic + I. d2 + α (d3 + π)(d2 + α) d3 + π CARRIERS’ EFFECTS ON INFECTIOUS DISEASES 717 Then i γ dL h β γα I0 + = + I0 dt d2 + α (d3 + π)(d2 + α) c d3 + π h pβ pγα (1 − p)γ i = S(βIc + γI) − (βIc + γI) + + d2 + α (d3 + π)(d2 + α) d3 + π hd + θ i 1 = R0 S − 1 (βIc + γI). b Using S ≤ b d1 +θ we know dL ≤ (R0 − 1)(βIc + γI) ≤ 0. dt dL b So dL dt ≤ 0 if R0 ≤ 1. Furthermore, dt = 0 ⇔ Ic = I = 0 or R0 = 1 and S = d1 +θ . dL Therefore the largest invariant set in the closure Γ̄ of Γ where dt = 0 is the singleton {P0 }. By LaSalle’s Invariance Principle [9], P0 is globally asymptotically stable in Γ, completing the proof. 5. Stability of the endemic equilibrium P ∗ . Theorem 5.1. If R0 > 1, then P ∗ is globally asymptotically stable with respect to the interior of Γ. Proof. To study the global stability of the endemic equilibrium, we make use of a Lyapunov function V of form V (S, Ic , I) = x1 (S − S ∗ ln S) + x2 (Ic − Ic∗ ln Ic ) + x3 (I − I ∗ ln I), (6) where x1 , x2 , x3 > 0 are constants to be specified. Note that V has a global minimum at P ∗ = (S ∗ , Ic∗ , I ∗ ) and V (S, Ic , I) − V (P ∗ ) is positive definite. We show that suitable constants x1 , x2 , x3 can be chosen such that the Lyapunov derivative of V is negative definite with respect to P ∗ . Direct calculation and applying the identity b = d1 S ∗ + θS ∗ + βIc∗ S ∗ + γI ∗ S ∗ lead to dV S∗ 0 I∗ I∗ =x1 (S 0 − S ) + x2 (Ic0 − c Ic0 ) + x3 (I 0 − I 0 ) dt S Ic I h i S∗ =x1 b − (d1 + θ)S − (βIc + γI)S − b + (d1 + θ)S ∗ + (βIc + γI)S ∗ S h βIc SIc∗ γISIc∗ + x2 (1 − p)(βIc + γI)S − (d2 + α)Ic − (1 − p) − (1 − p) Ic Ic i h ∗ βI SI c + (d2 + α)Ic∗ + x3 p(βIc + γI)S − αIc − (d3 + π)I − p I i γISI ∗ αIc I ∗ −p − + (d3 + π)I ∗ I I h S i S∗ − ∗) = x1 (d1 + θ)S ∗ (2 − S S h i ∗ ∗ ∗ + x1 (βIc + γI )S + x2 (d2 + α)Ic∗ + x3 (d3 + π)I ∗ h (βI ∗ + γI ∗ )S ∗2 γISIc∗ c − x1 + x2 (1 − p)βSIc∗ + x2 (1 − p) S Ic ∗i βIc SI ∗ αI I c + x3 pγSI ∗ − x3 . + x3 p I I 718 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI Positive constants x1 , x2 , and x3 are chosen as x1 = 1, x2 = (d3 + π)βS ∗ + γαS ∗ , (d2 + α)(d3 + π) x3 = γS ∗ . (d3 + π) (7) It can be verified that they satisfy relations −x1 + x2 (1 − p) + x3 p = 0, x1 γS ∗ − x3 (d3 + π) = 0, (8) ∗ x1 βS − x2 (d2 + α) + x3 α = 0. We re-group terms in dV dt such that dV dt = V1 + V2 + V3 , where S S∗ − ), ∗ S S V2 = x1 (βS ∗ Ic∗ + γI ∗ S ∗ ) + x2 (d2 + α)Ic∗ + x3 (d3 + π)I ∗ , V1 = (d1 + θ)S ∗ (2 − x1 (βIc∗ S ∗ + γI ∗ S ∗ )S ∗ − x2 (1 − p)βSIc∗ − x3 pγSI ∗ S x2 (1 − p)γISIc∗ x3 αIc I ∗ x3 pβSIc I ∗ − − . − Ic I I V3 = − We see that V1 ≤ 0 from the inequality x + x1 ≥ 2 for all x > 0, and that V1 = 0 if and only if S = S ∗ . We are left to show that V2 + V3 ≤ 0. We begin by examining V2 . Using the values for x1 , x2 , and x3 in (7), relations in (8), and the equilibrium relation (pd2 + α)Ic∗ = (1 − p)(d3 + π)I ∗ , (9) we can rewrite V2 as V2 = 2px3 γI ∗ S ∗ + 2(1 − p)x2 βIc∗ S ∗ + 4px3 βIc∗ S ∗ + 3(1 − p)α ∗ ∗ γI S . (pd2 + α) (10) Similarly, we can rewrite V3 as h (1 − p)x2 βS ∗2 Ic∗ i V3 = − x2 (1 − p)βIc∗ S − S h ∗ ∗2 i px γI S 3 + − x3 pγSI ∗ − S h x2 (1 − p)γISIc∗ x3 αIc I ∗ (1 − p)x2 γI ∗ S ∗2 i + −y − −y Ic I S h x2 (1 − p)γISIc∗ x3 pβIc SI ∗ (1 − p)x2 γI ∗ S ∗2 + − (1 − y) − − (1 − y) Ic I S px3 βS ∗2 Ic∗ i − , S where (1 − p)α (1 − p)pβS ∗ y= , 1−y = . (pd2 + α)(1 − p)x2 (pd2 + α)(1 − p)x2 Write V3 = Va + Vb + Vc + Vd , with each term representing the expression enclosed in a pair of big square brackets. We will estimate each term in V3 by applying the inequality a1 + a2 + · · · + an ≥ (a1 · a2 · · · an )1/n , for ai > 0. n CARRIERS’ EFFECTS ON INFECTIOUS DISEASES 719 We obtain (1 − p)x2 βS ∗2 Ic∗ Va = −x2 (1 − p)βIc∗ S − S p ≤ −2 (x2 (1 − p))2 (βIc∗ S ∗ )2 = −2(1 − p)x2 βIc∗ S ∗ , (11) and Vb = −x3 pγSI ∗ − p px3 γI ∗ S ∗2 ≤ −2 (x3 p)2 (γI ∗ S ∗ )2 = −2px3 γI ∗ S ∗ . S (12) Similarly, x2 (1 − p)γISIc∗ x3 αIc I ∗ (1 − p)x2 γI ∗ S ∗2 − −y Ic I S 1 3(1 − p)α ∗ ∗ ≤ −3[(x2 (1 − p))2 x3 αIc∗ y 2 (γI ∗ S ∗ )2 ] 3 = − γI S , (pd2 + α) Vc = −y (13) and (1 − p)x2 γI ∗ S ∗2 x2 (1 − p)γISIc∗ x3 pβIc SI ∗ − (1 − y) − Ic I S px3 βS ∗2 Ic∗ (14) − S 1 ≤ −4[(x2 (1 − p))2 (px3 )2 (1 − y)2 γI ∗ S ∗ 2βIc∗ S ∗2 γβI ∗ Ic∗ ] 4 = −4px3 βIc∗ S ∗ . Vd = −(1 − y) Therefore, (11) - (14) imply V3 ≤ −2(1 − p)x2 βIc∗ S ∗ − 2px3 γI ∗ S ∗ − 4px3 βIc∗ S ∗ − 3(1 − p)α ∗ ∗ γI S . (pd2 + α) (15) It follows from (10) and (15) that V2 + V3 ≤ 0 and thus dV dt ≤ 0. Furthermore, dV = 0 if and only if V = 0 and V + V = 0. Using (10) (15), we can show that 1 2 3 dt dV dV ∗ ∗ ∗ dt = 0 ⇔ (S, Ic , I) = (S , Ic , I ), and thus dt is negative definite with respect to P ∗ . The global stability of P ∗ follows from the classical stability theorem of Lyapunov. We remark that the form of Lyapunov function in (6) was motivated by those used in Korobeinikov and Maini [8], Guo [3], and Guo and Li [4]. 6. Impact of carriers on the transmission dynamics: chronic hepatitis B infection. To further illustrate the impact of disease carriers on the transmission dynamics, we use chronic hepatitis B infection as an example. Hepatitis B is a liver disease caused by the HBV virus. It is transmitted through sexual contact, the sharing of infected needles, or from mother to infant. Chronic HBV infection is more common: children infected with HBV rarely develop acute illness and up to 90% of infected children become chronically infected; adults infected with HBV usually recover from acute illness, but 5-10% will become chronically infected. About 30% of people infected with HBV do not show any symptoms. These people are the asymptomatic carriers. According to WHO statistics, about 2 billion people worldwide have been infected with HBV and about 350 million live with chronic infection. An estimated 600,000 persons die each year due to the acute or chronic consequences of hepatitis B. Safe and effective hepatitis B vaccines became available in 1982. Integration of the HBV vaccines into childhood immunization programs since 1991 has produced a great decline in the amount of children infected. In many countries where 8% to 15% of children used to become chronically infected 720 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI with HBV, vaccination has reduced the rate of chronic infection to less than 1% among immunized children [18, 16, 17]. Our model (1) can be used as a crude approximation for the transmission dynamics of chronic hepatitis B infection among an adult population. The compartment S contains individuals who are susceptible to HBV infection, compartment I contains individuals who are chronically infected with HBV and are symptomatic or have been tested and are aware of their condition, and compartment Ic contains individuals who are asymptomatic carriers of HBV and have no knowledge that they are infected. Recall that in this model, α represents the rate of diagnosis - the rate at which people carrying the disease are made aware of their infection, either through testing or through appearance of symptoms. Based on epidemiological data from WHO, CDC (US), and PHAC (Canada) [18, 16, 17], we have estimated the the values of our model parameters as follows: b = 90, 000, d1 = 1/80, d2 = d1 + 0.004, d3 = d2 , β = 1.5γ, π = 0.75. (16) We carry out numerical simulations of our model (1) in a hypothetical population of size 200, 000. We will vary key parameters to investigate the impact of asymptomatic carriers and HBV vaccinations. 1000 1000 800 800 Ic 600 600 200 0 Ic 400 400 200 I 0 10 20 30 (a) 40 50 t 0 I 0 10 20 30 40 t 50 (b) Figure 2. Simulation results showing the impact of testing and diagnosis of carriers. In (a), diagnosis rate α = 0.1. In (b), α = 0.5. Other parameter values are the same as in (16) In the first set of simulations, we fix a vaccination coverage rate at 70% where the vaccine has a success rate of 85%. This means that θ = 0.85 × 0.7 = 0.585. We vary the parameter α to see the effects of diagnosis rate at which carriers move into the infected class. We see in Figure 2 that, if only 1% of chronic carriers become aware of their disease, the number of symptomatically infected individuals decreases significantly, but the number of carriers is still high. This is not a desirable result as it is the carriers that are responsible for most of the new infections. If we increase α from 1% to 5%, a more dramatic change occurs in the disease dynamics: the number of carriers shows a much greater decline while the number of symptomatically infected remains low. This demonstrates that testing and diagnosis of carriers can be an effective control measure in high HBV prevalence countries. In the second set of simulations, we will fix α = 0.01 and vary θ to see the effect of increasing the vaccination rate. If we set θ = 0.1, we see in Figure 3 that though the number of symptomatically infected reduces rapidly, the number of carriers remains CARRIERS’ EFFECTS ON INFECTIOUS DISEASES 721 high. Increasing θ to 0.6 only slightly alters the disease dynamics; the number of carriers only shows a moderate decline. 1000 1000 Ic 800 800 600 600 400 400 200 0 200 I 0 10 Ic 20 30 40 50 t 0 I 0 (a) 10 20 30 40 50 t (b) Figure 3. Simulation results showing the impact of vaccination. In (a), vaccination rate θ = 0.1. In (b), θ = 0.6. Other parameter values are the same as in (16) Our model simulations demonstrate the challenges of chronic HBV infection: the existence of a large number of carriers who are infectious but show no symptoms. Because carriers do not show symptoms, they will not be part of any treatment program. Comparing our simulation results in Figures 2 and 3, we conclude that, in high HBV prevalence countries, testing and increasing awareness of carriers will have a much greater impact on the disease burden than increasing vaccination rates. While this conclusion may have practical implications for the control of chronic HBV infections, more realistic models that are specific for HBV infection and more detailed data need to be employed to further explore its significance in future studies. Acknowledgments. This work was based on the MSc Thesis of Darja Kalajdzievska at the University of Alberta [6]. DK acknowledges support of a NSERC PGSM Scholarship. MYL’s research is supported in part by grants from the Natural Sciences and Engineering Research Council of Canada (NSERC) and Canada Foundation for Innovation (CFI). Both authors acknowledge the financial support from the NCE-MITACS. The authors wish to thank Dr. Hongbin Guo for his help on the choice of Lyapunov function in the proof of Theorem 5.1. REFERENCES [1] M. Ghosh, P. Chandra, P. Sinha and J. B. Shukla, Modelling the spread of carrier-dependent infectious diseases with environmental effect, Appl. Math. Comput., 152 (2004), 385–402. [2] S. Goldstein, F. Zhou, S. C. Hadler, B. P. Bell, E. E. Mast and H. S. Margolis, A mathematical model to estimate global hepatits B disease burden and vaccination impact, Int. J. Epidemiol., 34 (2005), 1329–1339. [3] H. Guo, Global dynamics of a mathematical model of tuberculosis, Canadian Appl. Math. Quart., 13 (2005), 313–323. [4] H. Guo and M. Y. Li, Global dynamics of a staged progression model for infectious diseases, Math. Biosci. Eng., 3 (2006), 513–525. [5] J. M. Hyman and J. Li, Differential susceptibility and infectivity epidemic models, Math. Biosci. Eng., 3 (2006), 89–100. 722 DARJA KALAJDZIEVSKA AND MICHAEL Y. LI [6] D. Kalajdzievska, “Modeling the Effects of Carriers on the Transmission Dynamics of Infectious Diseases,” M.Sc. thesis, University of Alberta, 2006. [7] J. T. Kemper, The effects of asymptotic attacks on the spread of infectious disease: A deterministic model, Bull. Math. Bio., 40 (1978), 707–718. [8] A. Korobeinikov and P. K. Maini, A Lyaponov function and global properties for SIR and SEIR epedimiological models with nonlinear incidence, Math. Biosci. Eng., 1 (2004), 57–60. [9] J. P. LaSalle, “The Stability of Dynamical Systems,” Regional Conference Series in Applied Mathematics, SIAM, Philadelphia, 1976. [10] G. F. Medley, N. A. Lindop, W. J. Edmunds and D. J. Nokes, Hepatitis-B virus edemicity: Heterogeneity, catastrophic dynamics and control, Nat. Med., 7 (2001), 617–624. [11] R. Naresh, S. Pandey and A. K. Misra, Analysis of a vaccination model for carrier dependent infectious diseases with environmental effects, Nonlinear Analysis: Modelling and Control, 13 (2008), 331–350. [12] M. M. Riggs, A. K. Sethi, T. F. Zabarsky, E. C. Eckstein, R. L. Jump and C. J. Donskey, Asymptomatic carriers are a potential source for transmission of epidemic and nonepidemic Clostridium difficile strains among long-term care facility residents, Clin. Infect. Dis., 45 (2007), 992–998. [13] P. Roumagnac, et al., Evolutionary history of Salmonella typhi, Science, 314 (2006), 1301– 1304. [14] C. L. Trotter, N. J. Gay and W. J. Edmunds, Dynamic models of meningococcal carriage, disease, and the impact of serogroup C conjugate vaccination, Am. J. Epidemiol., 162 (2005), 89–100. [15] S. Zhao, Z. Xu and Y. Lu, A mathematical model of hepatitis B virus transmission and its application for vaccination strategy in China, Int. J. Epidemiol., 29 (2000), 744–752. [16] “The ABCs of Hepatitis,” Center for Disease Control and Prevention (CDC), 2009. Available from: http://www.cdc.gov/hepatitis/Resources/Professionals/PDFs/ABCTable_BW.pdf. [17] “Viral Hepatitis and Emerging Bloodborne Pathogens in Canada,” CCDR, 27S3, Public Health Agency of Canada (PHAC), 2001. [18] WHO, “Fact Sheet on Hepatitis B,” 2008. Available from: http://www.who.int/mediacentre/ factsheets/fs204/en/index.html. Received September 13, 2010; Accepted March 8, 2011. E-mail address: [email protected] E-mail address: [email protected]