Computational and Mathematical Methods in Medicine
Volume 2018 (2018), Article ID 7383170, 14 pages
https://doi.org/10.1155/2018/7383170

Research Article

Avian Influenza A (H7N9) Model Based on Poultry Transport Network in China

Juping Zhang,1,2 Wenjun Jing,1,2 Wenyi Zhang,³ and Zhen Jin1,2

¹Complex Systems Research Center, Shanxi University, Taiyuan, Shanxi 030006, China
²Shanxi Key Laboratory of Mathematical Techniques and Big Data Analysis on Disease Control and Prevention, Shanxi University, Taiyuan, Shanxi 030006, China
³Institute of Disease Control and Prevention of PLA, Beijing 100071, China

Correspondence should be addressed to Zhen Jin; jinzhn@263.net

Received 1 June 2018; Accepted 27 September 2018; Published 4 November 2018

Academic Editor: Konstantin Blyuss

Copyright © 2018 Juping Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In order to analyze the spread of avian influenza A (H7N9), we construct an avian influenza transmission model from poultry (including poultry farm, backyard poultry farm, live-poultry wholesale market, and wet market) to human according to poultry transport network. We obtain the threshold value for the prevalence of avian influenza A (H7N9) and also give the existence and number of the boundary equilibria and endemic equilibria in different conditions. We can see that poultry transport network plays an important role in controlling avian influenza A (H7N9). Finally, numerical simulations are presented to illustrate the effects of poultry in different places on avian influenza. In order to reduce human infections in China, our results suggest that closing the retail live-poultry market or preventing the poultry of backyard poultry farm into the live-poultry market is feasible in a suitable condition.

1. Introduction

Avian influenza A (H7N9) is a subtype of influenza viruses that have been detected in birds in the past. Until 2013 outbreak in China, no human infections with H7N9 viruses had ever been reported. But from March 31 to August 31, 2013, 134 cases had been reported in mainland China, resulting in 45 deaths [1]. However, the virus came back in November 2013 again. Afterwards the disease came back in November every year. In fact, the second outbreak occurred from November 2013 to May 2014. The third outbreak occurred from November 2014 to June 2015. The fourth outbreak occurred from November 2015 to June 2016. And the fifth outbreak occurred from September 2016 to May 2017 (NHFPC [1]). The disease causes a high death rate. In China, from March 2013 to May 2017, H7N9 has resulted in 1263 human cases including 459 deaths with a death rate of nearly 37%. In China, from September 2016 to May 2017, provinces with human cases are shown as Figure 1. H7N9 virus does not induce clinical signs in poultry and is classified as a low pathogenicity avian influenza virus (LPAIV) [2]. However, the virus can infect humans and most of the reported cases of human H7N9 infection have resulted in severe respiratory illness [3].

Figure 1: Provinces with avion influenza A (H7N9) from September 2016 to May 2017.

Jones et al. [4] demonstrated that interspecies transmission of H7N9 virus occurs readily between society finches and bobwhite quail but only sporadically between finches and chickens, and transmission occurs through shared water. Pantin-Jackwood et al. [3] showed that quail and chickens are susceptible to infection, shed large amounts of virus, and are likely important in the spread of the virus to humans, and it is therefore conceivable that passerine birds may serve as vectors for transmission of H7N9 virus to domestic poultry [4]. Zhang et al. [5] concluded that migrant birds are the original infection source. Many authors investigated the epidemic model which describes the transmission of avian influenza among birds and humans [8–15]. Liu et al. [16] constructed two avian influenza bird-to-human transmission models with different growth laws of the avian population, one with logistic growth and the other with Allee effect, and analyzed their dynamical behavior. Lin et al. [17] developed three different SIRS models to fit the observed human cases between March 2013 and July 2015 in China and found that environmental transmission via viral shedding of infected chickens had contributed to the spread of H7N9 human cases in China. Chen and Wen [18] took into account gene mutation in poultry. Guo et al. [19] proposed and analyzed an SE-SEIS avian-human influenza model. Mu and Yang [20] analyzed an SI-SEIR avian-human influenza model with latent period and nonlinear recovery rate. Gourley et al. [21] analyzed the patchy model for the spatiotemporal distribution of a migratory bird species. Bourouiba et al. [22] investigated the role of migratory birds in the spread of H5N1 avian influenza among birds by considering a system of delay differential equations for the numbers of birds on patches, where the delays represent the flight times between patches. In China, in 2013, to control the outbreak, local authorities of the provinces and municipalities, such as Jiangsu, Shanghai, and Zhejiang, temporarily closed the retail live-poultry markets which proved to be an effective control measure. Data indicate that the novel avian influenza A (H7N9) virus was most likely transmitted from the secondary wholesale market to the retail live-poultry market and then to humans [6, 7]. How is avian influenza A (H7N9) transmitted from live-poultry to human in China? In order to reveal the fact, the global network model of avian influenza A (H7N9) is constructed based on poultry transport network. The relationship between the global system and subsystem is analyzed. The corresponding risk indices are obtained. We study the impact of subsystems on the risk index of the global system. When the disease occurs, it can provide theoretical guidance for the global and local transport of poultry.

In this paper, we construct an avian influenza A (H7N9) transmission model from live poultry (including poultry farm, backyard poultry farm, live-poultry wholesale market, and wet market) to human for the heterogenous environments which affect the spread of H7N9. The remaining part of this paper is organized as follows: in Section 2, we first establish the model based on poultry transport network. We derive the threshold value of the model. In Sections 3 and 4, we discuss the different boundary and endemic equilibrium in the different thresholds. Section 5 gives the effect of different transmission rate on H7N9 by numerical simulation. Finally, concluding remarks are made in Section 6.

2. Model Based on Poultry Transport Network

The avian population is classified into poultry farm, backyard poultry farm, live-poultry wholesale market, and wet market (the retail live-poultry market). According to the present situation in China, the backyard poultry feeding is regarded as a large node, which is considered to be connected with all other nodes (except poultry farm) in network. The relationship diagram of poultry transport and contacts between human and poultry are described in Figure 2. Let , , and be the total number of poultry in ith poultry farm, jth live-poultry wholesale market, and kth wet market at time t, respectively, where , , and are classified into two subclasses: susceptible and infective, denoted by and , and , and and , respectively. Suppose there are L poultry farms, M live-poultry wholesales, and K wet markets, namely, . And they are independent of each other. Let be the total number of human at time t. The human population is classified into three subclasses: susceptible, infective, and recovered, denoted by , , and , respectively. All new recruitments of human population and avian population are susceptible. The avian influenza virus is not contagious from an infective human to a susceptible human. It is only contagious from an infective avian to a susceptible avian and a susceptible human. An infected avian keeps in the state of disease and cannot recover, but an infected human can recover, and the recovered human has permanent immunity. We neglect death rates of the poultry individuals during the transport process. The detailed description of dynamical transmission of H7N9 avian influenza is described in the following flowchart (Figure 3).

Figure 2: A possible network of H7N9 avian influenza.

Figure 3: Detailed transfer diagram on the dynamical transmission of H7N9 avian influenza.

The corresponding dynamical model can be seen in the following equation:

The interpretations of parameters of system (1) are described in Table 1. The parameters in system (1) are all nonnegative constants.

Table 1: Parameters of system (1).

Parameter
Interpretation

All new recruitments of the avian in ith poultry farm

All new recruitments of the avian in backyard poultry farm

The natural death rate (including slaughter) of the avian in different places

The disease-related death rate of the infected avian in different places

The transmission rate from infective avian to susceptible avian in different places

The transport rate of individuals from ith poultry farm to jth live-poultry wholesale market

The transport rate of individuals from jth live-poultry wholesale market to kth wet market

The transport rate of individuals from backyard poultry farm to jth live-poultry wholesale market

The transport rate of individuals from backyard poultry farm to kth wet market

All new recruitments of the human

The natural death rate of the human

The transmission rate from the infective avian in kth wet market to the susceptible human

The transmission rate from the infective avian in backyard farm to the susceptible human

The disease-related death rate of the infected human

The recovery rate of the infective human

The variation of the number of poultry in ith poultry farm is

and thus,

Similarly, the variation of the number of poultry in backyard poultry farm is

and thus,

The variation of the number of poultry in jth live-poultry wholesale market is

and thus,

The variation of the number of poultry in kth wet market is

and thus

The variation of the number of human is

and thus,

For convenience, we denote the positive solution of system (1) by .

Let then G is a positively invariant for system (1).

In order to find the disease-free equilibrium of system (1), we consider

System (12) has the unique positive equilibrium , where , , , , and . Thus, is the disease-free equilibrium of system (1).

According to the concepts of the next generation matrix and reproduction number presented in [23, 24], we define

where

Set , where represents the spectral radius of the matrix. Then, is called the reproduction number for system (1), where

If , then system (1) has the disease-free equilibrium , and is locally asymptotically stable.

Remark 1. If we do not consider backyard poultry farm, then system (1) becomes

A similar analysis is available for the above system.

3. Analysis of Subsystems of System (1)

Consider the poultry of the poultry farm subsystem, given by the first two equations of system (1), as follows

Let the right-hand side of system (17) equals to zero; when , we obtain

If , system (17) has the positive equilibrium . If , system (17) has only the disease-free equilibrium .

Remark 2.

(1)If , then each farm has the positive equilibrium.(2)If , then some of the poultry farms have the positive equilibrium, and the others have only the disease-free equilibrium.

Consider the poultry of the backyard poultry farm subsystem, given by the third and fourth equations of system (1), as follows

Let the right-hand side of system (19) equals to zero; when , we obtain

If , system (19) has the positive equilibrium . If , system (19) has only the disease-free equilibrium .
Consider the poultry of the live-poultry wholesale market subsystem, given by the fifth and sixth equations of system (1), as follows

Let the right-hand side of system (21) equals to zero; when , we can divide it into two cases.
If and , then we have

If , then system (21) has the positive equilibrium .
If or , then we obtain

where

Because , the solutions of the above equation are

If , system (21) has two positive equilibria (, ) and (, ). If and , system (21) has one positive equilibrium (, ). If , system (21) has no positive equilibrium.
Consider the poultry of the wet market (the retail live-poultry market) subsystem, given by the seventh and eighth equations of system (1), as follows:

Let the right-hand side of system (26) equals to zero, when , we can divide it into two cases.
If and , then we have

If , then system (26) has the positive equilibrium .
If or , then we have

where

Because , the solutions of the above equation are

If , system (26) has two positive equilibria (, ) and (, ). If and , system (26) has one positive equilibrium (, ). If , system (26) has no positive equilibrium.
Consider the human subsystem, given by the last three equations of system (1), as follows:

Since the first two equations of system (31) are independent of the variable , we only need to analyze the first two equations of system (31). Let the right-hand side of system (31) equals to zero, when , if or , then we have

4. Analysis of the Full System (1)

We analyze the following equivalent system:

For the sake of discussion, without loss of generality, we assume that a node has at least one link with the nodes in the next layer. So we have the following cases.

Case 1. If , system (33) has only the disease-free equilibrium . Namely, when all poultry has no avian influenza, human will not be infected with avian influenza.

Case 2. If , , and , system (33) has the boundary equilibrium

This shows that avian influenza A (H7N9) virus is most likely transmitted from the retail live-poultry market to humans when poultry has no disease in other types of farms.

Case 3. If , , and , system (33) has the boundary equilibrium as described next.
If and , system (33) has one boundary equilibrium:

If , system (33) has two boundary equilibria:

This shows that avian influenza A (H7N9) virus is most likely transmitted from the secondary wholesale market to the retail live-poultry market and then to humans [6, 7]. And there may be two boundary equilibria.

Case 4. If and , system (33) has the boundary equilibrium as described next.
If , , , and , system (33) has one boundary equilibrium:

If , , and , system (33) has two boundary equilibria:

If , , and , system (33) has two boundary equilibria:

If and , system (33) has four boundary equilibria:

When the poultry of poultry farms has no avian influenza, and the poultry of backyard poultry farm has avian influenza, we can obtain four cases. In four cases, human is most likely transmitted from the backyard poultry farm to the secondary wholesale market then to the retail live-poultry market, and finally to humans, or direct transmission from backyard poultry to humans.

Case 5. If and , system (33) has the boundary equilibrium as described next.
If , , , and , system (33) has one boundary equilibrium:

If , , and , system (33) has two boundary equilibria:

If , , and , system (33) has two boundary equilibria:

If and , system (33) has four boundary equilibria

When the poultry of poultry farms has avian influenza, and the poultry of backyard poultry farm has no avian influenza, we can obtain four cases. In four cases, human is most likely transmitted from the poultry farm to the secondary wholesale market, then to the retail live-poultry market, and finally to humans.

Case 6. If and , system (33) has the positive equilibrium as described next.
If , , , and , system (33) has one positive equilibrium:

If , , and , system (33) has two positive equilibria:

If , , and , system (33) has two positive equilibria:

If and , system (33) has four positive equilibria:

When the poultry of poultry farms has avian influenza, and the poultry of backyard poultry farm has avian influenza, we can obtain four cases. In four cases, human is most likely transmitted from the poultry farm and backyard poultry farm to the secondary wholesale market, then to the retail live-poultry market, and finally to humans, or direct transmission from backyard poultry to humans.

Remark 3. If we assume that there is an edge between each node of the upper layer and each node of the next layer, that is, each node of the upper layer transport poultry to each node of the next layer in the network, when , , or , according to the actual situation, it can be calculated and analyzed by a similar method. Hence, we omit them here.

5. Numerical Simulations

In this section, we first use , , and submodel to simulate. The course of the infected human is typically 1–4 weeks, and we assume that it is 2.5 weeks on average. Thus, the recovery rate of the infective human is . The disease-related death rate of the infected human is . The disease-induced death rates of poultry are assumed to be and . We assume that human can survive 70 years, and the poultry can survive 2 months in the farm, 1 week in wholesale market, 1 month in wet market, and 8 months in backyard farm, respectively. These rates also referred to removal due to slaughtering. Hence, these rates referred to removal due to slaughtering and the natural death. We take the parameter values as , , , and , respectively.

We estimate that the number of susceptible poultry population is between and , the number of infective poultry population is between 0 and 1000 in farm, the number of susceptible poultry population is between and , the number of infective poultry population is between 0 and 500 in live-poultry wholesale market, the number of susceptible poultry population is between and , the number of infective poultry population is between 0 and 100 in wet market, and the number of susceptible human population is between and in the region. So, we choose the initial values as , , , , and .

The difficulty in parameter estimations is that there is no scientifically or officially reported data of live-poultry transportation in China. The values of , , , and used in simulations may be estimated based on living habits of people of regions, the density of human population, and so on. Now, we assume that the transport rates of the backyard poultry are the same to each node, namely, and , where . Let and We assumed the replenishment rate to be 2 months which is the mean lifetime of farm poultry. Let , , , , and , respectively.

The transmission rates from the infective poultry in kth wet market to the susceptible human are , . The transmission rate from the infective poultry in backyard farm to the susceptible human is .

The transmission rates from infective poultry to susceptible poultry in different places are , , , and , respectively. Then, . Solution is asymptotically stable and converges to the disease-free equilibrium in Figure 4.

Figure 4: Solution is asymptotically stable and converges to the disease-free state value.

The transmission rates from infective poultry to susceptible poultry in different places are , , and , respectively. These parameters are fixed. is varied. Let , , and . From Figure 5(a), we can see that the beginning is almost the same, but the later is different. Therefore, the transmission rate has a small impact in the earliest stages but has an important impact on the late disease.

(a)

(b)

(c)

(d)

Figure 5: The plots display the changes of with varying.

The transmission rates from infective poultry to susceptible poultry in different places are , , and , respectively. These parameters are fixed. is varied. Let , , and . This only affects the number of infected humans, whereas it has no effect on the arrival time of the peak (Figure 5(b)). Therefore, preventing poultry of backyard poultry farm into the live-poultry market is feasible in a suitable condition.

The transmission rates from infective poultry to susceptible poultry in different places are , , and , respectively. These parameters are fixed. is varied. Let , , and . The bigger the , the more the human infected (Figure 5(c)). The secondary wholesale market plays an amplifier role.

The transmission rates from infective poultry to susceptible poultry in different places are , , and , respectively. These parameters are fixed. is varied. Let , , and . The impact is relatively small (Figure 5(d)).

6. Conclusion

In this paper, we construct the avian influenza transmission model from poultry (including poultry farm, backyard poultry farm, live-poultry wholesale market, and wet market) to human. We obtain the threshold value for the prevalence of avian influenza and the number of the boundary equilibria and endemic equilibria in different conditions. Numerical simulations show the effects of different transmission rates of different layer on the infected human. And, we can obtain the following cases:

(1)The poultry of poultry farm, backyard poultry farm, and poultry wholesale market have no avian influenza, but there is a possible outbreak of avian influenza in wet market (the retail live-poultry market), and avian influenza A (H7N9) virus is most likely transmitted from the retail live-poultry market to humans.(2)The poultry of poultry farm and backyard poultry farm has no avian influenza, but there is a possible outbreak of avian influenza in poultry wholesale market, and then avian influenza A (H7N9) virus is most likely transmitted from the poultry wholesale market to the retail live-poultry market and finally to humans.(3)The poultry of poultry farm has avian influenza, and the poultry of backyard poultry farm has no avian influenza, but there is a possible outbreak of avian influenza in poultry farm, and then avian influenza A (H7N9) virus is most likely transmitted from poultry market to poultry wholesale market, then to the retail live-poultry market, and finally to humans.(4)The poultry of poultry farm has no avian influenza, and the poultry of backyard poultry farm has avian influenza, but there is a possible outbreak of avian influenza in backyard poultry farm, and then avian influenza A (H7N9) virus is most likely transmitted from backyard poultry farm to poultry wholesale market, then to the retail live-poultry market, and finally to humans, or direct transmission from backyard poultry to humans.(5)The poultry of poultry farm and backyard poultry farm has avian influenza, but there is a possible outbreak of avian influenza in poultry farm and backyard poultry farm, and then avian influenza A (H7N9) virus is most likely transmitted from poultry farm and backyard poultry farm to poultry wholesale market, then to the retail live-poultry market, and finally to humans, or direct transmission from backyard poultry to humans.

Hence, the poultry of some nodes on network has avian influenza, and then all edges connected to the node should be cut off. It has a great inhibitory on preventing the spread of disease. So, the network of poultry transportation plays an important role in controlling avian influenza A (H7N9). Moreover, we find that there may have been avian influenza A (H7N9) among humans when there is avian influenza A (H7N9) in the retail live-poultry market, so closing the live-poultry market can reduce the spread of disease to humans at a certain time. In addition, we find that there may have been avian influenza A (H7N9) among humans when there is avian influenza A (H7N9) in the backyard poultry farm. But the spread of backyard poultry to human is quit complex. It can be either direct infection or indirect infection. In China, there are many backyard poultry, so there are still some difficulties in the prevention and control of avian influenza A (H7N9).

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partly supported by the National Sciences Foundation of China (11331009, 11314091, and 11501339).

References

National Health and Family Planning Commission of China (NHFPC), “National notifiable disease situation,” http://en.nhfpc.gov.cn/diseases.html, in English, http://www.nhfpc.gov.cn/zhuzhan/yqxx/lists.shtml, in Chinese.

Q. Li, L. Zhou, M. Zhou et al., “Epidemiology of human infections with avian influenza a (H7N9) virus in China,” New England Journal of Medicine, vol. 370, pp. 520–532, 2014.

M. J. Pantin-Jackwood, P. J. Miller, E. Spackman et al., “Role of poultry in the spread of novel H7N9 influenza virus in China,” Journal of Virology, vol. 88, no. 10, pp. 5381–5390, 2014.

J. C. Jones, S. Sonnberg, R. J. Webby, and R. G. Webster, “Influenza a (H7N9) virus transmission between finches and poultry,” Emerging Infectious Diseases, vol. 21, no. 4, pp. 619–628, 2015.

J. Zhang, Z. Jin, G. Q. Sun, X. D. Sun, Y. M. Wang, and B. Huang, “Determination of original infection source of H7N9 avian influenza by dynamical model,” Scientific Reports, vol. 4, p. 4846, 2014.

Y. Chen, W. Liang, S. Yang et al., “Human infections with the emerging avian influenza a H7N9 virus from wet market poultry: clinical analysis and characterisation of viral genome,” The Lancet, vol. 381, no. 9881, pp. 1916–1925, 2013.

C. Bao, L. Cui, M. Zhou, L. Hong, G. F. Gao, and H. Wang, “Live-animal markets and influenza A (H7N9) virus infection,” New England Journal of Medicine, vol. 368, no. 24, pp. 2337–2339, 2013.

S. Iwami, Y. Takeuchi, and X. N. Liu, “Avian-human influenza epidemic model,” Mathematical Biosciences, vol. 207, no. 1, pp. 1–25, 2007.

K. I. Kim, Z. G. Lin, and L. Zhang, “Avian-human influenza epidemic model with diffusion,” Nonlinear Analysis: Real World Applications, vol. 11, no. 1, pp. 313–322, 2010.

Y.-H. Hsieh, J. H. Wu, J. Fang, Y. Yang, and J. Lou, “Quantification of bird-to-bird and bird-to-human infections during 2013 novel H7N9 avian influenza outbreak in China,” PLoS One, vol. 9, no. 12, Article ID e111834, 2014.

J. C. Jones, S. Sonnberg, Z. A. Kocer et al., “Possible role of songbirds and parakeets in transmission of influenza a (H7N9) virus to humans,” Emerging Infectious Diseases, vol. 20, no. 3, pp. 380–385, 2014.

X. Ma and W. Wang, “A discrete model of avian influenza with seasonal reproduction and transmission,” Journal of Biological Dynamics, vol. 4, no. 3, pp. 296–314, 2010.

X. S. Wang and J. Wu, “Periodic systems of delay differential equations and avian influenza dynamics,” Journal of Mathematical Sciences, vol. 201, no. 5, pp. 693–704, 2014.

N. K. Vaidya and L. M. Wahl, “Avian influenza dynamics under periodic environmental conditions,” SIAM Journal on Applied Mathematics, vol. 75, no. 2, pp. 443–467, 2015.

S. Liu, L. Pang, S. Ruan, and X. Zhang, “Global dynamics of avian influenza epidemic models with psychological effect,” Computational and Mathematical Methods in Medicine, vol. 2015, Article ID 913726, 12 pages, 2015.

S. Liu, S. Ruan, and X. Zhang, “Nonlinear dynamics of avian influenza epidemic models,” Mathematical Biosciences, vol. 283, pp. 118–135, 2017.

Q. Y. Lin, Z. G. Lin, A. P. Y. Chiu, and D. He, “Seasonality of influenza a (H7N9) virus in China fitting simple epidemic models to human cases,” PLoS One, vol. 11, no. 3, Article ID e0151333, 2016.

Y. Chen and Y. Wen, “Global dynamic analysis of a H7N9 avian-human influenza model in an outbreak region,” Journal of Theoretical Biology, vol. 367, pp. 180–188, 2015.

S. M. Guo, J. Wang, M. Ghosh, and X. Z. Li, “Analysis of avian influenza a (H7N9) model based on the low pathogenicity in poultry,” Journal of Biological Systems, vol. 25, no. 2, pp. 1–16, 2017.

R. Mu and Y. Yang, “Global dynamics of an avian influenza A (H7N9) epidemic model with latent period and nonlinear recovery rate,” Computational and Mathematical Methods in Medicine, vol. 2018, Article ID 7321694, 11 pages, 2018.

S. A. Gourley, R. Liu, and J. Wu, “Spatiotemporal distributions of migratory birds: patchy models with delay,” SIAM Journal on Applied Dynamical Systems, vol. 9, no. 2, pp. 589–610, 2010.

L. Bourouiba, S. A. Gourley, R. Liu, and J. Wu, “The interaction of migratory birds and domestic poultry and its role in sustaining avian influenza,” SIAM Journal on Applied Mathematics, vol. 71, no. 2, pp. 487–516, 2011.

O. Diekmann, J. A. P. Heesterbeek, and J. A. J. Metz, “On the definition and the computation of the basic reproduction ratio in the models for infectious disease in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 1-2, pp. 29–48, 2002.