key: cord-0006378-gbui4lgl authors: Wen, Buyu; Wang, Jianpeng; Teng, Zhidong title: A discrete-time analog for coupled within-host and between-host dynamics in environmentally driven infectious disease date: 2018-02-26 journal: Adv Differ Equ DOI: 10.1186/s13662-018-1522-1 sha: dd253ec6a68c90e59479c7d2fdd71d61fdb68d8b doc_id: 6378 cord_uid: gbui4lgl In this paper, we establish a discrete-time analog for coupled within-host and between-host systems for an environmentally driven infectious disease with fast and slow two time scales by using the non-standard finite difference scheme. The system is divided into a fast time system and a slow time system by using the idea of limit equations. For the fast system, the positivity and boundedness of the solutions, the basic reproduction number and the existence for infection-free and unique virus infectious equilibria are obtained, and the threshold conditions on the local stability of equilibria are established. In the slow system, except for the positivity and boundedness of the solutions, the existence for disease-free, unique endemic and two endemic equilibria are obtained, and the sufficient conditions on the local stability for disease-free and unique endemic equilibria are established. To return to the coupling system, the local stability for the virus- and disease-free equilibrium, and virus infectious but disease-free equilibrium is established. The numerical examples show that an endemic equilibrium is locally asymptotically stable and the other one is unstable when there are two endemic equilibria. As is well known, viruses have caused abundant types of epidemic and occur almost everywhere on Earth, infecting humans, animals, plants, and so on. There are a large number of diseases, for example: influenza, hepatitis, HIV, AIDS, SARS, Ebola, MERS, etc., which are caused by viruses. Therefore, it is important to study viral infection, which can supply theory evidence for controlling diseases breaking out. In recent years, many authors have established and investigated the various kinds of viral infection dynamical systems which are described by differential equations and difference equations. Many important and valuable results were established and successfully applied to viral infections in practice. See, for example, and the references therein. In [1] [2] [3] [4] , the authors proposed a coupled within-host and between-host continuoustime dynamical system: where S, I, E, T, T * and V denote the numbers of susceptible and infectious individuals, the level of environment contamination, the densities of healthy cells and infected cells, and the parasite load, respectively. In system (1), the parameter A and denote the recruitment rate of susceptible and healthy cells, respectively. β is the infection rate of hosts in a contaminated environment. μ is the natural mortality rate of host. θ is the rate of contamination. γ is the clearance rate. k is the infection rate of cells. m and d denote the natural and infection-induced mortality rates of infected cells, respectively. p is the parasite reproduction rate by an infected cell. c is the within-host mortality rate of parasites. α is the induced mortality rate of host. It is assumed that the rate of environment contamination is proportional to the number of infected hosts and the parasite load V within a host, which has the form θ VI. The function g(E) denotes the rate at which an average host is inoculated. In [1] , for system (1) the authors introduced a slow time variable τ = εt, where 0 < ε 1. In this case, t as a fast time variable. The authors further considered the parameters associated with the dynamics at the population level to be small based on the assumption that the between-host dynamics occur on a slower time scale than that of parasite-cell dynamics within the host. Let = ε˜ , β = εβ, μ = εμ, α = εα, θ = εθ and γ = εγ . Then, under the faster time variable t and the slower time variable τ , system (1) can be written as the two singular perturbation systems, see systems (6) and (7) given in [1] , respectively. Using the techniques from the singular perturbation theory in [23] , the authors in [1] analyzed system (1) by analyzing the corresponding fast and slow dynamics, and the fast and slow dynamics can be analyzed using the fast and slow time subsystems, see systems (8) and (9) given in [1] , respectively. Here, we see that the concepts of fast and slow time systems were introduced in [23] early. We easily see that system (1) also can be described by using fast and slow time variables t and τ in the following form: (2) From this system, by the transformations t = τ ε for the last three equations of system (2) and τ = εt for the first three equations of system (2), we can easily obtain two similar singular perturbation systems and two similar fast and slow time subsystems as systems (6), (7) and systems (8) , (9) proposed in [1] . In recent years, more and more attention was paid on the discrete-time epidemic models. The reasons are as follows. Firstly, because the statistic data about infectious disease is collected by day, week, month, or year, so it is more direct, more convenient and more accurate to describe the epidemic by using discrete-time models than continuous-time models. Secondly, it is very difficult to solve a nonlinear differential equation with a given initial condition to obtain the exact solution. Thus, for many practical requirements, such as numerical calculation, it is often necessary to discretize a continuous model into the corresponding discrete model. Therefore, we see that the discrete-time analog is also alike important for studying coupled system (1) . At the present time, there are various discretization methods to discretize a continuous model, including the standard methods, such as the Euler method, the Runge-Kutta method, and some other standard finite difference schemes, and the non-standard finite difference scheme, which is originally developed by Mickens (see [24] [25] [26] ). In this paper, we propose a discrete-time analog for above continuous-time system (2) by using discretization method of Micken's non-standard finite difference scheme, the model is given as follows: and where system (3) denotes the slow dynamics with slow time t, and system (4) denotes the fast dynamics with fast time s. However, in slow system (3) there is a fast time term V (s+1), and in fast system (4) there is a slow time term g(E(t + 1)). Therefore, systems (3) and (4) form a coupled system. It is clear that in order to study the dynamical properties of coupled systems (3)- (4) we can firstly analyze the fast and slow two subsystems which are determined by two time scales t and s. In other words, we can treat the within-host subsystem (4) as the fast system and the between-host subsystem (3) as the slow system. For fast time s and slow time t, we may assume that there exists a certain relation between s and t. For example, it may be assumed that there is a large enough integer K such that s = Kt in slow system (3) and t = [ s K ] in fast system (4), where [ s K ] denotes the maximum integer which is not more than s K . Thus, coupled systems (3)-(4) will acquire the following form: However, since there are terms V (Kt + 1) and g(E([ s K ] + 1)) in system (5), it is very difficult to readily investigate system (5) theoretically. Therefore, in this paper we firstly separate the coupled systems (3)-(4) into a fast system and a slow system by using the idea of limit systems. For fast system (4), we assume that the environmental contamination E keeps a constant owing to the faster time scale is enough quick. For slow system (3), we can assume that V (s) steadies to an equilibriumV (E). Thus, coupled systems (3)-(4) are separated into the following two isolated subsystems: and ⎧ ⎪ ⎪ ⎨ ⎪ ⎪ ⎩ S(t + 1) -S(t) = A -βE(t)S(t + 1) -μS(t + 1), whereV (E) is given in Sect. 3. In this paper, for fast system (6) we will investigate the dynamical behaviors, including the positivity, boundedness, basic reproduction number, the existence of equilibria and the local stability of equilibria by using the discretization method. For slow system (7), we will investigate the dynamical properties, including the positivity, boundedness, the existence of disease-free equilibrium, only a unique endemic equilibrium, and two endemic equilibria, and the local asymptotic stability for the disease-free and endemic equilibria. Next, we will investigate the dynamical behaviors for the coupled systems (3)-(4) basing on the research results obtained for the fast and slow subsystems. We will establish some criteria on the local asymptotic stability for the infection-and disease-free equilibrium, virus infectious but disease-free equilibrium and the endemic equilibrium. Furthermore, for the special cases which there is a unique endemic equilibrium, and two endemic equilibria in coupled systems (3)-(4), by means of the numerical examples, we will indicate that the unique endemic equilibrium may be locally asymptotically stable, and an endemic equilibrium may be locally asymptotically stable but the other one may be unstable. This paper is organized as follows. In Sects. 2 and 3, fast system (6) and slow system (7) are discussed. Some criteria on the positivity, boundedness, existence of equilibria and local asymptotic stability are stated and proved. In Sect. 4, coupled systems (3)-(4) is discussed. Some criteria on the existence of equilibria and local asymptotic stability are stated and proved. In Sect. 5, the numerical examples are given. Lastly, a discussion is presented in Sect. 6. We firstly introduce the following lemmas on the quadratic and cubic polynomial equations which are given in [27] . Lemma 2 All roots λ of the cubic equation f (λ) = λ 3 + a 2 λ 2 + a 1 λ + a 0 = 0 satisfy |λ| < 1 if and only if the following conditions are satisfied: For coupled systems (3)-(4), function g(E) is assumed to satisfy the following basic assumption. (H) g(E) is defined for all 0 ≤ E ≤ 1 and is continuously differentiable, which satisfies g(0) = 0, g(E) ≥ 0, g (E) > 0 and g (E) ≤ 0 for all 0 ≤ E ≤ 1. From the biological background of system (6) , it is assumed that any solution (T(s), T * (s), V (s)) satisfies the following initial value: Firstly, on the positivity and boundedness of the solutions and the existence of nonnegative equilibria for system (6) we have the following results. Proof System (6) is equivalent to the following: If (8) is satisfied, then from system (9) it follows that (T(1), T * (1), V (1)) exists and is positive. Hence, by induction, we see that (T(s), T * (s), V (s)) exists and is positive for all s ≥ 0. From the first equation of system (9), we have Hence, Particularly, when T(0) ≤ T 0 , we also have T(s) ≤ T 0 for all s > 0. From the second equation of system (9), we have By (10), we can obtain Particularly, we also can prove that when T(0) ≤ T 0 and T * (0) ≤T * , then T * (s) ≤T * for all s > 0. From the third equation of system (9) , it follows that . By (10) and (11), we can obtain Similarly, we also can prove that V (s) ≤V for all s > 0 if T(0) ≤ T 0 , T * (0) ≤T * and V (0) ≤V . Therefore, solution (T(s), T * (s), V (s)) with initial value (8) is ultimately bounded. This completes the proof. Remark 1 Define a set as follows: Then from the proof of Lemma 3 we see that is a positive invariable and globally attractive set for system (6) . We define the baseline within-host reproduction number as follows: . Lemma 4 Let E = 0, then system (6) always has infection-free equilibrium B 0 (T 0 , 0, 0), and when R f > 1, system (6) has a unique infectious equilibrium B * (T,T * ,V ). Proof Let E = 0. It is obvious that system (5) has a unique infection-free equilibrium B 0 (T 0 , 0, 0). For infectious equilibrium B * (T,T * ,V ), we have Hence, This shows that, when R f > 1, an infectious equilibrium B * (T,T * ,V ) exists and is unique. This completes the proof. Let E > 0 in system (6) . If (T(E),T * (E),V (E)) is a nonnegative equilibrium of system (6), then we have We further havȇ andT(E) satisfies the following equation: where Since Eq. (13) has always two positive real solutions given by the following: Note that a 1 (E) = g (E) m+d pm > 0, and Owing to a 2 > 0, we have from (15) for any E ≥ 0,T + (E) > 0 andT -(E) < 0. In (14) , when E = 0, by calculating we obtain Therefore, andT SinceT + (E) > 0 andT + (0) ≥ T 0 from (16), we obtainT + (E) > T 0 for all E > 0. But, from the first equation of (12), we haveT + (E) ≤ T 0 . This leads to a contradiction. SinceT -(E) < 0 andT -(0) ≤ T 0 from (17), we obtainT -(E) 0, system (6) has a unique positive equilibrium B 1 (T(E),T * (E),V (E)), wherȇ By calculating, we also havȇ Furthermore, when E → 0 + , from (14) and (17), by calculating we can obtain Summarizing the above discussions, we finally get the following result. Lemma 5 Let E > 0, then fast system (6) always has a unique infected equilibrium Next, we discuss the stability of the infection-free equilibrium and infectious equilibrium for system (6) . We have the following theorems. which is equivalent to the following: The characteristic equation of system (18) is We have eigenvalues λ 1 = 1 1+m , and λ 2 and λ 3 satisfying the equation f (λ) = 0. It is clear that f (0) < 1, f (-1) > 0 and When R f < 1, then f (1) > 0. By Lemma 1, it follows |λ 2 | < 1 and |λ 3 | < 1. Therefore, equilibrium B 0 is locally asymptotically stable. When R f > 1, then f (1) < 0. Since lim λ→∞ f (λ) = +∞, we see that f (λ) = 0 has a root λ 3 ∈ (1, +∞). This implies that equilibrium B 0 is unstable. This completes the proof. Theorem 2 Let E = 0 in system (6) . If R f > 1, then infectious equilibrium B * is locally asymptotically stable. The proof of Theorem 2 will be given in Theorem 3 as the special case with E = 0. We hence omit it here. Proof We will prove Theorem 3 in the case E ≥ 0, and when E = 0 we assume that R f > 1 and In fact, from the third equation of (12) and . This shows that Therefore,T(E) ≤ c(m+d) pk . For convenience, letT(E) =T,T * (E) =T * andV (E) =V . The linearization system of system (6) at equilibrium B 1 is The characteristic equation of system (20) is Obviously, |a 3 | < 1. According to (19) , we obtain Further, we have we have and according to (19) , we further have Therefore, by Lemma 2 all roots λ of the equation g(λ) = 0 satisfy |λ| < 1. Thus, equilibrium B 1 is local asymptotically stable. Particularly, when E = 0 we also see that equilibrium B * is local asymptotically stable. This completes the proof. Now, we consider slow system (7) . We assume that fast system (6) has steadied at the equilibriumB(T(E),T * (E),V (E)), whereB(T(E),T * (E),V (E)) is defined as follows: It is clear that fast system (6) is locally asymptotically stable in equilibriumB, that is, when V (E(t + 1)) in slow system (7). From the biological background of system (7), we assume that any solution (S(t), I(t), E(t)) of system (7) satisfies the initial value Firstly, on the positivity and boundedness of the solutions and the existence of nonnegative equilibria for slow system (7) we have the following lemmas. The solution (S(t), I(t), E(t)) of system (7) with initial value (21) is positive for all t ≥ 0 and ultimately bounded. Furthermore, 0 ≤ E(t) ≤ 1 for all t ≥ 0. Proof We know that system (7) is equivalent to the following form: 1+γ +θI(t)V (E(t+1)) . When t = 0, we prove that (S(1), I(1), E(1)) is positive. In fact, according to the first equation of system (22) (1)) . Let x = E(1), then the above equality becomes f (x) = 0, where When E > 0 we havê Hence, According to assumption (H), (15) andT(x) =T -(x), we know thatV (x) > 0 and where a 1 (x) = m+d pm g (x) < 0. This showsT (x) > 0, which leads toV (x) < 0. Therefore, f (x) > 0, this implies that f (x) = 0 has at most one solution. Hence, there is a unique x > 0 such that f (x) = 0. Therefore, (S(1), I(1), E(1)) exists uniquely and is positive. Furthermore, from 0 ≤ E(0) ≤ 1 we also have 0 < E(1) < 1. When t = 1, by a similar argument to above, we can prove that (S(2), I(2), E(2)) exists uniquely and is positive. Owing to 0 < E(1) < 1, we also have 0 < E(2) < 1. Using induction, for any t ≥ 0, we know that (S(t), I(t), E(t)) exists uniquely and is positive. Furthermore, we finally have 0 < E(t) < 1 for all t ≥ 0. Hence, lim sup t→∞ S(t) ≤ A μ . From the second equation of system (22), we have It follows that lim sup t→∞ I(t) ≤ βA μ(μ+α) . This completes the proof. Noticing the slow system (7) and the quick system (6) is linked by the terms V (s) and g (E) . Then if R f > 1, there are the steadied infectious equilibrium in the quick system (6) . In this case, we give the basic reproduction number R s for slow system (7): Obviously, we see that if R f < 1 then R s < 0, if R f = 1 then R s = 0 and if R f > 1 then R s > 0. Furthermore, when R s ≥ 1 then we must have R f > 1. We denote the functions as follows: Based on the reproduction number R s , we have the following lemma. (7) always has a disease-free equilibrium P 0 ( A μ , 0, 0). Proof It is obvious that system (7) has disease-free equilibrium P 0 ( A μ , 0, 0). The endemic equilibrium P * (S,Ī,Ē) satisfies equation Hence, we havē which is equivalent to H(Ē) = 0. By calculating, we have Furthermore, by calculating we see that, when 0 < E ≤ 1, Hence, H (E) < 0 for all 0 < E ≤ 1. This shows that H(E) is as above a convex function. If condition (a) holds, then from H(0) = 0 and H M > 0, we easily see that H(E) = 0 has a unique positive rootĒ. Hence, endemic equilibrium P * (S,Ī,Ē) exists and is unique. If condition (b) holds, then from H(0) > 0, it follows that H(E) = 0 has a unique positive rootĒ, and hence endemic equilibrium P * (S,Ī,Ē) also exists and is unique. Assume that condition (c) holds, then owing to H(0) < 0 and H M > 0, H(E) = 0 has only two positive roots. Hence, system (7) has only two endemic equilibria P 1 and P 2 . Lastly, we prove that system (7) has only disease-free equilibrium P 0 ( A μ , 0, 0) if one of the conditions (d) and (e) holds. In fact, when H M < 0 we see that H(E) = 0 has no root. When R s = 1, then H(0) = 0. Therefore, by H M = 0 there is only disease-free equilibrium P 0 ( A μ , 0, 0). This completes the proof. Remark 2 From the proof of Lemma 7, we further see that when system (7) has a unique endemic equilibrium P * (S,Ī,Ē), then since H(E) is a above convex function for 0 ≤ E ≤ 1, we also have H (Ē) ≤ 0. Next, we discuss the stability of the disease-free equilibrium and endemic equilibrium for system (7) . We have the following theorems. (a) If R s ≤ 1, then disease-free equilibrium P 0 of system (7) is locally asymptotically stable. Proof The linearization part of system (22) is The characteristic equation of system (23) is . When R f ≤ 1, we haveV (0) = 0. Hence f (λ) = 0 has roots λ 1 = 1 1+μ+α and λ 2 = 1 1+γ . This implies that disease-free equilibrium E 0 is locally asymptotically stable. When R f > 1, we haveV (0) = V * . Owing to R s ≤ 1, by calculating we obtain We see that by Lemma 1 two roots λ 1 and λ 2 of f (λ) = 0 satisfy |λ 1 | < 1 and |λ 2 | < 1. Therefore, disease-free equilibrium P 0 is locally asymptotically stable. When R s > 1, we have and lim λ→∞ f (λ) = ∞, there is aλ > 1 such that f (λ) = 0. Therefore, P 0 is unstable. In order to discuss the stability of unique endemic equilibrium P * (S,Ī,Ē) of system (7), we need to introduce the following assumption. (A) (ϕ + 1)(n + 1)(w -1)M + (1 + 2r + ϕ)N > 0, where M = ϕw + ϕn + nw + nwϕ + r(1 + μ)ϕ -2, N = wϕ + nϕ + nw + 2(ϕ + n + w) + 2rϕ, ϕ = μ + α, w = γ -θĪF (Ē) and n = μ + βĒ. Theorem 5 Assume that (A) holds and one of conditions (a) and (b) in Lemma 7 holds. Then unique endemic equilibrium P * (S,Ī,Ē) of system (7) is locally asymptotically stable. Proof The linearization part of system (22) is Here, by calculating we have The characteristic equation of system (24) is where a 1 = -1 1 + μ + βĒ , Obviously, |a 3 | < 1. By calculating, we can obtain Further, we have (-1) 3 f (-1) = 1 + 1 1 + μ + βĒ Now, we prove |b 0 | + b 2 > 0 and |b 0 |b 2 > 0, where b 0 = 1a 2 3 and b 2 = a 1 a 3a 2 . Since by calculating we can obtain From assumption (A) we obtain |b 0 | + b 2 > 0. By Lemma 2, all roots λ of the equation f (λ) = 0 satisfy |λ| < 1. Therefore, equilibrium P * locally asymptotically stable. It is difficult to discuss the local stability for the case of two positive equilibria P 1 and P 2 in condition (c) of Lemma 7 by using the linearization method. However, we can give the following conjecture. Assume that R s < 1 and H M > 0. Let P 1 (S 1 ,Ī 1 ,Ē 1 ) and P 2 (S 2 ,Ī 2 ,Ē 2 ) be two positive equilibria of slow system (7) withĒ 1 <Ē 2 . Then P 2 is locally asymptotically stable, and P 1 is unstable. Now, we return to coupled systems (3)-(4). IfD(S,Ĩ,Ẽ,T,T * ,Ṽ ) is the equilibrium of coupled systems (3)-(4), then we have From Lemmas 4, 5 and 7, we have the following result. (1) Coupled system (3)-(4) always has a disease-free and infection-free equilibrium D 0 ( A μ , 0, 0, T 0 , 0, 0). On the stability of equilibriumD(S,Ĩ,Ẽ,T,T * ,Ṽ ) of coupled system (3) and (4), we have the following definition. (1)D is said to be stable, if for any constant > 0, there is a δ = δ( ) > 0 such that, for any initial point (S 0 , I 0 , E 0 , T 0 , T * 0 , V 0 ) at time s = 0 and t = 0 satisfying |S 0 -S| < δ, |I 0 -Ĩ| < δ, |E 0 -Ẽ| < δ, |T 0 -T| < δ, |T * 0 -T * | < δ, and |V 0 -Ṽ | < δ, one has |S(t) -S| < δ, |I(t) -Ĩ| < δ, |E(t) -Ẽ| < δ, |T(s) -T| < δ, |T * (s) -T * | < δ, and |V (s) -Ṽ | < δ, for all t ≥ 0 and s ≥ 0. (2)D is said to be locally asymptotically stable, ifD is stable and there is a constant δ > 0 such that, for any solution (S(t), I(t), E(t), T(s), T * (s), V (s)) with initial point (S 0 , I 0 , E 0 , T 0 , T * 0 , V 0 ) at time s = 0 and t = 0 satisfying |S 0 -S| < δ, |I 0 -Ĩ| < δ, |E 0 -Ẽ| < δ, |T 0 -T| < δ, |T * 0 -T * | < δ, and |V 0 -Ṽ | < δ, one has Furthermore, by applying the theory of limit equations, from Theorems 1, 2 and 4, we have the following result. (1) If R f < 1 and R s ≤ 1, then equilibrium D 0 ( A μ , 0, 0, T 0 , 0, 0) is locally asymptotically stable, and if R f > 1, then D 0 is unstable. (2) If R f > 1 and R s ≤ 1, then equilibrium D 1 ( A μ , 0, 0,T,T * ,V ) is locally asymptotically stable, and if R s > 1, then D 1 is unstable. Proof In fact, the linearization systems of coupled systems (3)-(4) at equilibrium D 0 is It is easy to see that when R f < 1 and R s ≤ 1, then, by conclusion (a) of Theorem 1 and conclusion (a) of Theorem 4 and from the first three equation of (25), we know (X(t), Y (t), Z(t)) → (0, 0, 0) as t → ∞ and further from the last three equations of (25), we also have (U(s), V (s), W (s)) → (0, 0, 0) as s → ∞. Therefore, D 0 is locally asymptotically stable. When R f > 1, then, by conclusion (b) of Theorem 1, we see that equilibrium (0, 0, 0) of the last three equations of (25) is unstable. In addition, when R s > 1 we also have R f > 1. Therefore, D 0 is unstable. The linearization system of coupled system (3)-(4) at equilibrium D 1 is It is clear that when R s ≤ 1, by conclusion (a) of Theorem 4, from first three equation of (26), we know (X(t), Y (t), Z(t)) → (0, 0, 0) as t → ∞. By Theorem 2 and from last three equations of (26), we further have (U(s), V (s), W (s)) → (0, 0, 0) as s → ∞. Therefore, D 1 is locally asymptotically stable. When R s > 1, by conclusion (b) of Theorem 4, we see that equilibrium (0, 0, 0) of the first three equations of (26) is unstable. Therefore, D 1 is unstable. This completes the proof. However, to establish the criteria of stability for endemic equilibrium D * and two positive equilibria D 2 and D 3 is very difficult. We here only give the following conjectures. Assume that R s < 1 and H M > 0. Then D 3 is locally asymptotically stable and D 2 is unstable. In the following section, we will give a numerical example to show that Conjectures 2 and 3 may be right. In this section, we give the numerical examples to discuss assumption (A) and the stability for two endemic equilibria P 1 and P 2 of slow system (7) and for two equilibria D 2 and D 3 of coupled systems (3)-(4). We assume that s = Kt and t = [ s K ] in coupled systems (3)-(4). For convenience, we choose K = 366 and function g(E) = wE. The parasite reproduction rate of infected cell p is chosen as a free parameter. The rest of the parameters in coupled systems (3)-(4) are chosen in Table 1 . By calculating, we obtain (ϕ + 1)(n + 1)(w -1)M + (1 + 2r + ϕ)N = 0.000026 > 0. Therefore, assumption (A) is satisfied. We first take the parasite reproduction rate of infected cell p = 1200. By calculating, we see that the basic reproduction numbers R f . = 1.33 > 1 and R s . = 5 > 1. We see that system The numerical simulations given in Fig. 1 show that equilibrium P * and D * is locally asymptotically stable only when R s > 1. We next take the parasite reproduction rate of infected cell p = 950. By calculating, we see that the basic reproduction numbers R f . = 1.0556 > 1 and R s . = 0.8333 < 1. Furthermore, we also have H M = 1891.39 > 0. Hence, slow system (7) has two endemic equilibria P 1 (9503.94, 248.02, 0.0348) and P 2 (6900.24, 1549.87, 0.2995) and coupled systems (3) Fig. 2 show that equilibria P 2 and D 3 are locally asymptotically stable and equilibria P 1 and D 2 are unstable. Therefore, Conjecture 1, Conjecture 2 and Conjecture 3 may be right. In this paper we studied a discrete coupled within-host and between-host models (3)-(4) in environmentally driven infectious disease obeying Micken's non-standard finite difference scheme. Since there are two fast and slow time scales in the model, and the fast time scale is sufficiently quicker than the slow time scale, the model is separated into a fast system (6) and a slow system (7) . The basic properties for fast system (6) , including the existence of infection-free equilibrium B 0 , infected equilibrium B * (when E = 0) and infected equilibrium B 1 (when E > 0), the positivity and ultimate boundedness of the solutions with positive initial values, are established. Under assumption (H), the local stability of equilibria for system (6) is completely determined by basic reproduction number R f . That is, when E = 0 in system (6), if R f < 1 then B 0 is locally asymptotically stable, and if R f > 1 then B 0 is unstable and B * is locally asymptotically stable. When E > 0 in system (6), then infectious equilibrium B 1 exists always and also is locally asymptotically stable. For slow system (7), the basic properties on the existence of disease-free equilibrium P 0 , unique endemic equilibrium P * and two positive equilibria P 1 and P 2 , and the positivity and ultimate boundedness of the solutions with positive initial values are established. The sufficient conditions on the local stability of disease-free equilibria P 0 and unique endemic equilibria P * are established by virtue of basic reproduction number R s , the quantity H M and condition (A). However, it is very difficult to discuss the local stability of two endemic equilibria P 1 and P 2 . Here we only show the local stability of P 1 and P 2 by the numerical examples in Sect. 5 . We see that assumption (A) is a pure mathematical condition. It is only used in the proofs of theorems on the local stability of endemic equilibria P * to obtain |b 0 | + b 2 > 0 (see the proof of Theorem 5). Generally, we expect that the local stability of equilibria of slow system (7) can be determined only by basic reproduction number R s . Therefore, an open problem is whether condition (A) can be thrown off in Theorem 5. Furthermore, we also do not obtain the global asymptotic stability of equilibria for system (7) . The reason is that the construction of Lyapunov function is very difficult. For whole coupled systems (3)-(4), the basic properties on the existence of infectionand disease-free equilibrium D 0 , viral infection and disease-free equilibrium D 1 , unique endemic equilibrium D * and two endemic equilibria D 2 and D 3 , and the local stability of equilibria D 0 and D 1 are established, respectively. However, it is difficult to discuss the local stability for unique endemic equilibrium D * , and two endemic equilibria D 2 and D 3 . Here, we only show the local stability of D * , D 2 and D 3 by the numerical examples in Sect. 5. Comparing the results established in this paper with the results obtained in [1, 3] , we see that the dynamical properties of equilibria for discrete-time model (3)-(4) and continuous-time model (2) (see Theorems 1-3 in [1] ) in fast time and slow time subsystems, respectively, are very oncoming. This shows that discrete-time model (3)-(4), as a discrete-time analog of continuous-time model (2), is fairly appropriate. Particularly, we can use model (3)-(4) to calculate the numerical approximative solution of model (2) in a neighborhood of equilibrium. In addition, in this paper we further investigate the dynamical properties for whole coupled systems (3)-(4), such as the existence of equilibrium and the local stability of equilibrium. Coupled within-host and between-host dynamics and evolution of virulence A mathematical model for coupling within-host and between-host dynamics in an environmentally infectious disease Emerging disease dynamics in a model coupling within-host and between-host systems A model for coupling within-host and between-host dynamics in an infectious disease Modeling Infectious Diseases in Humans and Animals Virus Dynamics: Mathematical Principles of Immunology and Virology Modelling within-host HIV-1 dynamics and the evolution of drug resistance: trade-offs between viral enzyme function and drug susceptibility Linking within-and between-host dynamics in the evolutionary epidemiology of infectious diseases Evaluating the importance of within-and between-host selection pressures on the evolution of chronic pathogens. Theor Epidemic dynamics and host immune response: a nested approach Evolution of virulence: interdependence, constraints, and selection using nested models. Theor An immuno-epidemiological model with threshold delay: a study of the effects of multiple exposures to a pathogen Mathematical analysis of a virus dynamics model with general incidence rate and cure rate Global stability of a five-dimensional model with immune responses and delay Global properties of basic virus dynamics models Global dynamics of an in-host viral model with intracellular delay Global dynamics of a virus dynamical model with general incidence rate and cure rate Modeling the adaptive immune response in HBV infection Stability and Hopf bifurcation for a virus infection model with delayed humoral immunity response A model of HIV-1 infection with two time delays: mathematical analysis and comparison with patient data Viral dynamics model with CTL immune response incorporating antiretroviral therapy Global properties of a discrete viral infection model with general incidence rate Dynamical Systems. Lecture Notes in Mathematics A note on an NSFD scheme for a mathematical model of respiratory virus transmission Application of Nonstandard Finite Difference Scheme Dynamics consistency: a fundamental principle for constructing nonstandard finite difference scheme for differential equation Mathematical Models in Population Biology and Epidemiology We are very grateful to the reviewers for their very helpful comments and careful reading of our manuscript. This research was supported by the National Natural Science Foundation of China (Grant Nos. 11771373, 11401512, 11261056) . A the recruitment rate of individuals 4Ref. [5] β the infection rate of hosts in a contamination 0.0006 Ref. [5] μ the natural mortality rate of host 0.0004 Ref. [5] α the induced mortality rate of host 0.0004 Ref. [5] the rate which an average host is inoculated g(E) = 4 × 10 5 E Refs. [5, 6, 23] θ the rate of contamination 1.5 × 10 -10 Ref. [5] γ clearance rate 0.015 Ref. [ No potential conflict of interest was reported by the authors. The main idea of this paper was proposed by ZT. BW prepared the manuscript initially and performed all the steps of the proofs in this research. JW performed the numerical examples and simulations. All authors read and approved the final manuscript. Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Received: 26 August 2017 Accepted: 9 February 2018