A General Model of Zoonosis Where the Diseases Can Only Be Transmitted from Animals to Humans ()
1. Introduction
Infectious diseases seriously endanger human health and the development of community economy all over the world [1] [2] [3], historically, infectious diseases have caused millions of deaths. SARS attracted worldwide attention in 2003, more than 8000 infected people have been found worldwide, resulting in more than 900 deaths [4]; Dengue fever in Thailand has reached 5284 people, and the disease is prevalent in Ecuador, Paraguay, Australia and other countries to varying degrees by 2006 [5]. Thus, researching and controlling the dynamic process of infections is very significant to society.
Among numerous infectious diseases, zoonosis is an important factor affecting human economic development and population mortality. There are about 200 confirmed zoonoses in the world, such as avian influenza, rabies, malaria, brucellosis and so on. From 1996 to 2015, 30,300 human rabies cases were reported [6]. Every year, about 50,000 people are newly infected with brucellosis [7]. What’s more, there are about 300 to 500 million new cases because of malaria annual, among which maybe one million to three million deaths, most of them are children [8]. Malaria kills a child every 40 seconds, costing the world more than 2000 young lives a day [8] [9].
Many scholars have studied infectious disease models. At the beginning Kermack and Mckendrick opened the door of infectious disease research. In 1927 and 1932, they separately established the classical SIR [10] and SIS [11] model of plague. In 2002 van Driessche and Watmough [12] provided a method of establishing Lyapunov function which furnishes strong theoretical support for the latter study of the dynamic process of the infections. On this basis, more and more scholars established different models for different infectious diseases [13] [14] [15]. However, a basic model on zoonosis has not been established at now. Many models consider the vaccination [12] [16], but there are only a few people get vaccinated in poor area, even there is not an effective vaccination for some newly discovered zoonoses. Some models consider human transmission to human or animal [17] [18], in fact there are some zoonosis that can only be transmitted from animals to humans such as rabies, brucellosis, tapeworm disease and so on. What’s more, few animals are treated after they get sick in the real life. This paper studies on this kind of zoonosis.
In this article, a general model of the interaction between infected animals and humans are established. Referencing to classical model [10] [11] we generalize a SIS model which consider the cure rate those who are infected on human population and a SI model on animal population. The organizational structure of the article is shown below. In Section 2 a mathematical expression of the model is constructed. The third section studies the global dynamics of disease from the perspective of disease-free equilibrium and endemic equilibrium. Section 4 makes numerical simulation to explain the consistency of the theory and numerical results, after that the sensitive analysis is performed. Finally, the conclusion is shown in Section 5.
2. Model Formulation
We consider both animal and human, then classify each of them into tow subclasses: susceptible and infectious. Defining
as the susceptible and infectious human,
as the susceptible and infectious animals respectively. If a susceptible human individual is contact with infected animals, this human individual may be infected. If this one received effective treatment, he can change back to a susceptible individual, but if he isn’t cure in time, he is in the risk of death from infection. The model proposed in this paper is a system composed of four ordinary differential equations:
(2.1)
The dynamic transmissions for the model of infectious of animal and humans are demonstrated in Figure 1.
Figure 1 illustrate that the people and animal population import new susceptible individuals by the birth of baby with the number of
and
. Under the influence of infected animals, the susceptible individuals will transfer to infected one with the rate of
and
. In the people population, the infected individual can return back to a susceptible individual with the rate
by medical treatment. The infected individuals will die and then be removed from the population for the reason of nature factor with probability of
and disease factor with probability of
. All the parameters with subscript 1 define the people population property while the parameters with subscript 2 define the animal population property. Based on the reality, we can know there is not a parameter is negative for human population.
3. Dynamic Analysis
Following, we express the derivative of t with
. Adding the first two equations, we can get
(3.1)
There is
(3.2)
Similarly
(3.3)
Figure 1. Transmission diagram of zoonosis among animals and humans.
Accordingly, we can conclude an invariant set for the model
(3.4)
3.1. The Basic Reproduction Number
Obviously, setting the right-hand side of the four differential equations to zero we can find the equilibrium. When set
The model (1) has a unique disease-free equilibrium
:
(3.5)
when
,we can find a unique equilibrium
:
(3.6)
where
(3.7)
Based on this, the knowledge of next generation matrix and basic reproduction number are proposed [10], we now consider the follow auxiliary model
(3.8)
It can be written as
(3.9)
Following the recipe from [10] we have
(3.10)
is the next generation matrix and
is the spectral of
. We call
as the basic reproduction number
(3.11)
3.2. Dynamic Analysis on Disease-Free Equilibrium
Theorem 1. If
,the disease-free equilibrium
of the model is locally asymptotic stability, or else,
is instability.
Proof:
The disease-free equilibrium Jacobian matrix of model (1) is
(3.12)
The eigenvalues of
can be obtained from the above equation as follows:
(3.13)
when
,we can know
,then every eigenvalue of the model has the negative real parts,
is locally asymptotic stability. Otherwise,
is unstable and there will be exist diseases case.
Theorem 2. If
,then the disease-free equilibrium
is globally asymptotic stability.
Proof:
The Jacobian matrix of model (2) on the disease-free equilibrium
is
(3.14)
It is easy to see that when
,there is
,where
Since the invariant set of the model (1) is
,choose positive and small enough
,
,and
,such for all
,there is
Consider the following auxiliary model
(3.15)
Similarly, the Jacobian matrix of model (3) on the disease-free equilibrium
is
(3.16)
Since
are infinitely close to zero,
is continuous, on account of
,It’s easy to derive
for positive and small enough
. By the comparison principle [19], it can be concluded that
and
,as
. By the theory of [20], which defined that
is called asymptotically autonomous―with limit equation
if
,
,locally uniformly in
.
Thus, there exist
and
as
. Accordingly, when
,the disease-free equilibrium of the model
is global asymptotic stability.
3.3. Dynamic Analysis on Endemic Equilibrium
We have calculated the endemic equilibrium
in the above article. Next, the global exponentially stable of
will be shown.
Theorem 3. If
,then the endemic equilibrium
of the model is globally asymptotic stability.
Proof:
On the endemic equilibrium
,we derive the model (1) satisfied the follow equation:
(3.17)
Thus, the model (1) can be converted into
(3.18)
Define
(3.19)
It is easy to get
(3.20)
The corresponding differential equation is
(3.21)
We construct the Lyapunov function as
(3.22)
Therefore, the derivative of L is as follow
(3.23)
It is easy to see,
for all
. In addition, if and only if
,
. That is to say,
if and only if
,
,
,
. Thus, the endemic equilibrium
of the model is globally asymptotic stability in the interior of
.
4. Numerical Simulation
In this section, we discuss two situations in which animal population size is larger or smaller than population size. For the purpose of proving the universality of the model, we analyze the global stability by establishing the simulation data including animals and setting the parameters and initial values of the population size. After that, we give the sensitive analysis to see what is the most important factor to infect the existence or extinction of disease.
4.1. Numerical Results
Because it is not specific to infectious diseases, we use dimensionless parameter to conduct the numerical experiments. All the parameters values and initial values are displayed in Table 1, when people size is larger than animal size.
Since the transmission rate of animal to animal is bigger than animal to human and disease mortality rate is bigger than nature mortality rate, for
,we set
and
,this moment
the corresponding results are shown in Figure 2(a). For
,we set
and
,this moment
the corresponding results will be displayed in Figure 2(b).
We can conclude from Figure 2(a) that when
,the number of infected animal individuals is stable at 0, that is to say the disease-free equilibrium
is globally stable. Otherwise, from Figure 2(b) when
,the result shows that the quantity of infected animal individuals is finally stable around 10.9, in the meantime, we calculate
,it illustrates that the endemic equilibrium
is globally stable when
.
Next, we analysis the situation of the people size is smaller than animal size. Values of parameters and initial values are shown in Table 2.
For
,we set
and
,now we calculate
Table 1. Values of parameters and initial values when people size is larger than animal size.
Table 2. Values of parameters and initial values when people size is smaller than animal size.
Figure 2. Numerical results of the infected animal individuals when people size is larger than animal size.
and the corresponding results are shown in Figure 3(a). For
,
and
,now
and the corresponding results are shown in Figure 3(b).
Similarly, Figure 3(a) we know when
,the disease-free equilibrium
is globally stable. When
,the number of infected animal individuals is stable around 545.5, and we calculate
,it illustrates that the endemic equilibrium
is globally stable when
.
Figure 3. Numerical results of the infected animal individuals when people size is smaller than animal size.
4.2. Sensitive Analysis
With the view of investigating the role of each parameter on
,we use the following formula:
where
is the parameters of
,and
shows the sensitivity of
to
. The greater the absolute value of
,the greater the effect of
on
. It is easy to see
and
are equal to 1,
and
are as follows:
Obviously,
,
. It follows that increase the mortality of animal can promoting the extinction of disease. However, massive killing of animals is not advocated for those animals like dogs, sheep. So we can find ways to reduce the infection rate and birth rate.
5. Conclusion
In this article, we construct a general model of zoonoses transmitted from animals to humans. By dynamic analysis we can know that when
,the disease-free equilibrium
is globally stable. When
,
is globally stable. Theory tells us that infectious diseases will not unrestrictedly grow, but if not take action to solve the epidemic of infectious diseases, it will exist all the time and then continuously make harm to human society. We can take measures to control the birth rate of animals to reduce the incidence of infectious diseases.