The Mathematical Model for Streptococcus suis Infection in Pig-Human Population with Humidity Effect

Inthira Chaiya，Kamonchat Trachoo，Kamsing Nonlaopon and Din Prathumwan，＊

1Department of Mathematics，Faculty of Science，Mahasarakham University，Mahasarakham，44150，Thailand

2Department of Mathematics，Faculty of Science，Khon Kaen University，Khon Kaen，40002，Thailand

Abstract: In this paper， we developed a mathematical model for Streptococcus suis， which is an epidemic by considering the moisture that affects the infection.The disease is caused by Streptococcus suis infection found in pigs which can be transmitted to humans.The patients of Streptococcus suis were generally found in adults males and the elderly who contacted pigs or who ate uncooked pork.In human cases， the infection can cause a severe illness and death.This disease has an impact to the financial losses in the swine industry.In the development of models for this disease，we have divided the population into 7 related groups which are susceptible pig compartment，infected pig compartment，quarantined pig compartment，recovered pig compartment， susceptible human compartment， infected human compartment，and recovered human compartment.After that， we use this model and a quarantine strategy to analyze the spread of the infection.In addition， the basic reproduction number R0 is determined by using the next-generation matrix which can analyze the stability of the model.The numerical simulations of the proposed model are illustrated to confirm the results from theorems.The results showed that there is an effect from moisture to the disease transmission.When the moisture increases the disease infection also increases.

Keywords: Streptococcus suis; mathematical model; stability analysis;infectious disease;reproduction number

1 Introduction

Nowadays，one of the major problems in the world is zoonotic pathogen such as Rift Valley fever，SARS，pandemic influenza H1N1 2009，Yellow fever，Avian Influenza(H5N1)and(H7N9)，West Nile virus[1]，the Middle East respiratory syndrome coronavirus(MERSCoV)[2]，and novel coronavirus(COVID-19) [3].Streptococcus suisis an infectious agent which is widely found in pigs around the world.This pathogen is generally found in the upper respiratory tract，genital，and alimentary tracts of the pigs[4，5].Streptococcus suisis a Gram-positive bacterium and a facultative anaerobic bacterium which is a primary agent of sepsis and meningitis in pigs[6，7].In pig farms，the disease can transmit from pig to pig rapidly.The reported rate of mortality in pigs fromStreptococcus suisis about 20%[8].It can classify at least 35 serotypes，for examples，serotypes 1 to 31，33，and 1/2[9，10].Serotypes 1 to 9 and 14 are most commonly found in pigs around the world[11，12].

In 1968，the first human who was infected withStreptococcus suiswas found in Denmark[4，5].This situation confirms thatStreptococcus suiscan transmit from pig to human.Serotype 2 is reported that it is frequently found in human infections[13，14]，but in some cases，the infections are caused by other serotypes.Kerdsin et al.reported that serotypes 2 and 14 are involved in human cases that found in Thailand[15，16].Dutkiewicz et al.presented that disease found in human caused byS.suisserotypes 2，4，5，14，16，21 and 24[12].S.suisserotype 9 which is the most common found in pigs in Europe was found the first case infection in human in May 2013 locating in northern of Thailand[17].In addition，Hatrongjit et al.reported the first case of serotype 31 in human[18]in the central region of Thailand.The common symptoms of this disease in human are fever，headache，meningitis，septicemia，arthritis，pneumonia，and hearing loss[19].

To prevent the spread of the disease， the mathematical models become powerful tools for describing the dynamics of the disease.They can forecast the future behaviors of the disease with many assumptions.The solutions of the model can be simulated satisfying the given parameters from the hypothesis.Recently，many mathematical models have been used to describe the behavior of infectious disease[20–24].Rahman et al.proposed the nonlinear SEITR fractional order model of tuberculosis disease by using Atangana-Baleanu derivative which consists of 5 population groups[25].Arfan et al.[26]presented the fractional order model to predict the dynamic of tumor with drug intervention.The nonlinear fractional order mathematical model for predicting the dynamic of COVID-19 was created by Shah et al.[27].Almuqrin et al.proposed the fractional model for forecasting the transmission of Ebola virus in bats population by using Atangana-Baleanu fractional derivative[28].Srivastava et al.improved the mathematical model for studying the dynamic of diabetes patients[29].Alzaid et al.[30] studied the dynamic of HIV-1 infection by using the mathematical model and showed that the mathematical model is efficient.The fractional order model of HIV with source term was proposed by Shah et al.to consider the behavior of CD4+T-cells which depending on the concentration of the viral load[31].The symmetry and asymmetry concepts can be linked to the epidemic model[32–34].

A number of research studied the disease fromStreptococcus suison epidemiology and medicine.However，a few pieces of the research proposed and studied the mathematical model forStreptococcus suis.Shen et al.[35] proposed the SIQRW model to explore the outbreaks of S.Suis.Giang et al.proposed the stochastic model and SEI model to predict the behavior of the disease and fitted the model parameters with collected data [9].The proposed models consider the transmission on pig population only.However，to our knowledge，no research considering the disease transmission between pigs and humans forStreptococcus suis.

In this paper，we propose a mathematical model to describe the disease transmission ofStreptococcus suisin humans and pigs.The novelty of this work is the improved mathematical model from one species to two species population.We classify the population of pigs and humans into 4 classes and 3 classes，respectively.This work contains the analysis of the model and numerical simulations to study the dynamic behavior of the disease.We also find the basic reproduction numbers to explain the rapid transmission and study the stability of the solution of the model.

This paper is constructed as follows; In Section 2， we present the mathematical model which expresses the Streptococcus suis transmission.The analysis of the proposed model is shown in Section 3 which includes a basic reproduction number，equilibrium points，and their stability.The numerical simulations are given in Section 4 and follow by a conclusion in Section 5.

2 Model Formulation

In this section， we describe the model formulation ofStreptococcus suisdisease transmission in humans and pigs.Based on the classical epidemiology model，we propose a new generalized model ofStreptococcus suisinfection which is SIQR-SIR model.There are some facts about the environmental factors such as temperature and relative humidity[36–38].So，we consider the effect of the moisture in the air for disease transmission in the pig farm.

The investigated population is divided into two subpopulations.These are pig population and human population.Then，both subpopulations of pigs and human are separated into four classes and three classes， respectively.These are a pig susceptible class (Sp)， pig infectious class (Ip)， pig isolated class(Qp)，pig recovery class(Rp)，human susceptible class(Sh)，human infectious class(Ih)，and human recovery class(Rp).We assume that the total population at timetisN(t).We have that

Sp(t)+Ip(t)+Qp(t)+Rp(t)=Np(t)，

Sh(t)+Ih(t)+Rh(t)=Nh(t)，

and

Sp(t)+Ip(t)+Qp(t)+Rp(t)+Sh(t)+Ih(t)+Rh(t)=Np(t)+Nh(t)=N(t)，

whereNp(t)is the total population of pigs，andNh(t)is the total population of human at timet.

We can write the dynamic equations forNp(t)andNh(t)as the following.

The new infection in pigs can be described .whereβ1is the transmission coefficient per unit of time per pig in the susceptible class contact with infectious class[39].The new infection in human can be described byβ2ShIp+β3ShIhwhereβ2is the transmission coefficient per unit of time per person in the susceptible class in contact with pig infectious class，β3is the transmission coefficient per unit of time per person in the susceptible class contact with the infectious class.In our model，we investigate the relative humidity for considering the infection of this disease in pigs.Then，the new infection in pigs can be expressed byβ1MSpIpinstead ofβ1SpIp，whereMis the moisture in the air(relative humidity).In addition，we assume that the disease cannot be transmitted from human to pigs.Throughout this paper，we assume that all parameters are positive constants.

The model of disease transmission byStreptococcus suiscan be represented by the system of differential equations as follows:

wherebis the pig removal rate，ais the pig death rate induced by the disease，δis the rate from infectious class to isolated class in pigs，eis the transition rate from isolated class to recovery class，μis the human natural death rate，αis the death rate induced by disease，andγis the transition rate from infectious class to recovery class.The model form is valid only ifa≤bandα≤μ.

Eq.(1a) expresses the rate of change of the population of susceptible pig.The number of susceptible pigs increases only via birthN1described in the first term on the right-hand side.On the other hand， the population of them decreases due to the removal rate of pig(b)， transmission ofStreptococcus suisfrom infected pig to susceptible pig(β1)and moisture effect(M).

Eq.(1b)expresses the rate of change of the population of the infected pig.On the right-hand side，the first term represents the number of infected pigs increases due to susceptible pig becoming infected pig with the rateβ1Mfollowed by the negative effects of disease death rate(a)，removal rate(b)，and isolated rate(δ).

Eq.(1c)expresses the rate of change of the population of the isolated pig.On the right-hand side，the first term represents the increase in the population of isolated pig due to the isolation of infected pig with the rate Δ followed by the negative effects of disease death rate(a)， removal rate(b)， and recovered rate(∈).

Eq.(1d) expresses the rate of change of the population of the recovered pig.On the right-hand side，the first term represents the increase in the population of the recovered pig due to treatment with the rate∈，and the last term represents the removal rate(b)of pigs.

Eq.(1e) expresses the rate of change of the population of susceptible human.The number of susceptible pigs increases only via birthN2described in the first term on the right-hand side.On the other hand，the population of them decreases due to the transmission ofStreptococcus suisfrom infected pig to susceptible human(β2)and infected human to susceptible human(β3).The last term represents the natural death rate(μ)of human.

Eq.(1f) expresses the rate of change of the population of the infected human.On the righthand side， the first term and second term represent the number of the infected human increases due to susceptible human becoming infected human by contact with the infected pigs(β2) and the infected human(β3)， respectively.The remaining terms are the negative effects of disease death rate(α)，recovered rate(γ)，and natural death rate(μ).

Eq.(1g) expresses the rate of change of the population of recovered human.On the right-hand side，the first term represents the increase in the population of recovered human due to treatment with the rate(γ)，and the last term represents the natural death rate(μ)of human

LetSp(0) =Sp0，Ip(0) =Ip0，Qp(0) =Qp0，Rp(0) =Rp0，Sh(0) =Sh0，Ih(0) =Ih0，Rh(0) =Rh0.A flowchart of the SIQR-SIR model of pigs and human which describe by the system(1)is shown in Fig 1.

Figure 1:Flowchart of SIQR-SIR(Susceptible-Infectious-Isolation-Recovery-Susceptible-Infectious-Recovery)model for human and pigs

The values of population classesSp(t)，Ip(t)，Qp(t)，Rp(t)，Sh(t)，Ih(t)，Rh(t) at timetis nonnegative numbers.

Therefore，the positive invariant of the system(1)is

3 Model Analysis

3.1 Equilibria

The equilibria are obtained by setting all equations of the system(1)to be zero.

By solving the system(4)，We get three equilibrium points:

i) Disease-free equilibrium

ii) Pig disease free equilibrium

whereA≡α+γ+μ.

Note that，E2exists ifβ3N2-μA ＞0.

iii)The endemic equilibrium

where

withX≡β3N2-A(β2I*p+μ)andB≡a+b+δ.

Note that，E*exists ifβ1MN1-bB ＞0.

3.2 Basic Reproduction Number

The basic reproduction number is the expected number of secondary cases produced by a single infection in a completely susceptible population.To compute the basic reproduction number and to study the local stability of the equilibrium，we use the next generation matrix method[40，41].We define x′=f(x)-v(x)，where x = [Ip，Qp，Ih]Twheref(x)is the matrix of new infection terms，andv(x)is the matrix of transfer terms into compartment and out of compartment as:

The corresponding Jacobian matrices are

Then，we obtain

We obtain three eigenvalues ofFV-1as:

Hence，the spectral radii are，and

3.3 The Stability of Disease-Free Equilibrium(E1)

Theorem 1.The disease-free equilibrium point(E1)is locally asymptotically stable if R01＜1 andR02＜1.

Proof.The model system(1)has the Jacobian matrix at the pointE1as:

The corresponding eigenvalues ofJ(E1)are

Note thatλ6is negative ifR02＜1 andλ7is negative ifR01＜1.

Therefore，the disease-free equilibrium(E1)of the model system(1)is locally asymptotically stable due to all negative real part eigenvalues ifR01＜1 andR02＜1.□

3.4 The Stability of Pig Disease Free Equilibrium(E2)

Theorem 2.The pig disease free equilibrium(E2) is locally asymptotically stable ifR01＞1 andR02＜1.

Proof.The pig disease-free equilibrium

exists and is positive ifR01＞1.

The model system(1)has the Jacobian matrix at the pointE2as:

The corresponding eigenvalues ofJ(E2)are

Note thatλ5is negative ifR02＜1 whileλ6andλ7have negative real part ifR01＞1.

Therefore，the pig disease free equilibrium(E2)of the system(1)is locally asymptotically stable ifR01＞1 andR02＜1.□

3.5 The Stability of Endemic Equilibrium(E*)

Theorem 3.The endemic equilibrium(E*)is locally asymptotically stable if R01＜1 andR02＞1.

Proof.The endemic equilibrium

exists and is positive ifR02＞1.

The model system(1)has the Jacobian matrix at the pointE*as:

where

The corresponding eigenvalues ofJ(E*)are

where

and

Next，we consider

and

We obtain thatλ4andλ5has negative real part ifR02＞1.

Then，let us consider

Then，we get

IfR01＜1，I*pandI*hare positive i.e.，R02＞1，then we obtain

Thus，the eigenvaluesλ5andλ6have a negative real part ifR01＜1 andR02＞1.

Therefore，the endemic equilibrium(E*)of the model system(1)is locally asymptotically stable ifR01＜1 andR02＞1.□

We can draw the bifurcation diagram to describe the basic reproduction numbers and stability as shown in Fig.2.

Figure 2:The bifurcation region

4 Numerical Examples and Discussion

The numerical results of the system (1) are computed by using MATLAB with the given initial values:

Sp(0)=4057，Ip(0)=1000，Qp(0)=0，Rp(0)=0，Sh(0)=50000，Ih(0)=1000，Rh(0)=0.

The numerical results of the system(1)with the parameter values as shown in Tab.1.Note that the parameters and initial values are obtained from data in[35]and[42].

Table 1:Parameter values of the system(1)

The solution trajectories tend to the disease-free equilibrium(E1)which satisfy Theorem 1 with the remaining parameter valuesμ=0.9，α=0.9，γ=0.9，M=0.9，a=0.9，β2=0.1，and Δ =0.9 as shown in Fig.3.The calculated reproduction numbers of this case areR01= 0.9568＜1 andR02=0.5625＜1.

Figure 3:The simulation results of the system(1)，(a)The plots of asymptotic population density of pig susceptible population(Sp)，pig infectious population(Ip)，pig isolated population(Qp)，and pig recovery population (Rp) with respect to time t tends to the equilibrium (E1)， and (b) The plots of asymptotic population density of human susceptible population (Sh)， human infectious population(Ih)，and human recovery population(Rh)with respect to time t tends to the equilibrium(E1)

The results showed that the number of pig infectious cases dramatically decreased in the first quarter.After that， the decreasing of pig infectious cases slowly decreased.Then， it tended to the equilibrium value.Fig.3a， the number of pig susceptible population increased after infectious cases were recovered.In addition， the isolated population increased when the number of infectious cases increased to control the spread of the disease.Then， it decreased after the number of infectious cases decreased.In human case (Fig.3b)， the number of infectious populations increased in the first quarter，then the number increased and approached to the equilibrium number.Therefore，the parameters of this case provided that the disease died out.

The solution trajectories tend to the disease-free equilibrium(E2)which satisfy Theorem 2 with the remaining parameter valuesμ= 0.6，α= 0.5，γ= 0.4，M= 0.3，a= 0.7，β2= 0.1，and Δ = 0.8 as shown in Fig.4.The calculated reproduction numbers of this case areR01= 2.5833＞1andR02=0.2125＜1.

The results showed that the number of pig infectious cases dramatically decreased in the first quarter.After that， the decreasing of pig infectious cases slowly decreased.Then， it tended to the equilibrium value.The number of pig susceptible population increased after infectious cases decreased.In human case (Fig.4b)， the number of infectious populations dramatically increased in the first quarter， then it decreased and approached to the equilibrium number.Therefore， the result of this case indicated that the disease in pig died out while it still appeared in human.

Figure 4:The simulation results of the system(1)，(a)The plots of asymptotic population density of pig susceptible population (Sp)，pig infectious population (Ip)， pig isolated population (Qp)， and pig recovery population (Rp)with respect to time t tends to the equilibrium (E2)， and (b) The plots of asymptotic population density of human susceptible population (Sh)， human infectious population(Ih)，and human recovery population(Rh)with respect to time t tends to the equilibrium(E2)

The numerical results of the system(1)with the remaining parameter valuesμ=0.5，α=0.6，γ=0.1，M= 0.9，a= 0.3，β2= 0.01，and Δ = 0.1 withR01= 0.9568＜1 andR02= 1.2473＞1.The solution trajectories tend to the endemic equilibrium(E*) which satisfy Theorem 3 as shown in Fig.5.The number of pig susceptible population decreased in the first five quarters then it increased and approached to the equilibrium.The number of pig infectious population increase in the first quarter then it tends to equilibrium after five quarters.In human case，the number of human infectious population increased sharply in the first two quarters.After that，it decreased and tend to equilibrium value.The number of recovery population increased in the first quarter then it decreased and approached to equilibrium value.The result of this case showed that the disease still appeared in both human and pigs.

To consider the effect of moisture， the numerical results of the system (1) with the remaining parameter valuesμ= 0.5，α= 0.6，γ= 0.1，a= 0.01，β2= 0.01， and Δ = 0.01.The solutions trajectories are computed with various values of relative humidityMas 0.1， 0.3， 0.5， 0.7，and 0.9 is shown in Fig.6.

In Fig.6， the result indicates that when the relative humidity(M) increases the pig infectious population density also increases.Particularly， the trajectories of pig infectious population decrease and approach to constant values whenM= 0.1， 0.3，and 0.5.On the other hand，the number of pig infectious population increase and go to equilibrium value whenM= 0.7 and 0.9.Moreover，whenM=0.1 and 0.3 the trajectories approach zero.This means that the disease will die out when there is less moisture in the air.Therefore，the result confirms that there is an effect of the moisture in the air to the infection of Streptococcussuisas mentioned in the literature.

Figure 5:The simulation results of the system(1)，(a)The plots of asymptotic population density of pig susceptible population(Sp)，pig infectious population(Ip)，pig isolated population(Qp)，and pig recovery population (Rp) with respect to time t tends to the equilibrium (E*)， and (b) The plots of asymptotic population density of human susceptible population (Sh)， human infectious population(Ih)，and human recovery population(Rh)with respect to time t tends to the equilibrium(E*)

Figure 6:The simulation results of the pig infectious population (Ip) for various values of relative humidity M with respect to time t

5 Conclusions

We have proposed the mathematical model for predicting the disease transmission ofStreptococcus suisbetween pigs and humans，which is developed from the SIQR model that considered only in the pig population.The pigs-human epidemiology model ofStreptococcus suistransmission which is SIQR-SIR model is investigated.The model combined with susceptible(Sp)，infection(Ip)，isolation(Qp)，and recovery(Rp)for pig population and susceptible(Sh)，infectious(Ih)，and recovery(Rh)for human population.Moreover， we studied the effect of moisture in the air on disease transmission.The equilibria of the model were analyzed.Then， the next-generation matrix was used to find the basic reproduction numbersR01，andR02.We obtained the conditions of the equilibrium pointsE1，E2，andE*which are locally asymptotically stable in Theorem 1，Theorem 2，and Theorem 3，respectively.The results show that ifR01＜1 andR02＜1，the equilibrium pointsE1is locally asymptotically stable whileE2is locally asymptotically stable ifR01＞1 andR02＜1.Finally，the endemic equilibrium(E*)is locally asymptotically stable ifR01＜1 andR02＞1.The bifurcation diagram is also presented.The numerical examples confirmed the results of Theorems.The results indicated that moisture affects disease transmission.When the air has higher moisture the transmission of the pathogen is also high.Moreover，the model can be improved by considering infection of newborn pigs from infectious adult female pigs.

Acknowledgement:This research was financially supported by the young researcher development project of Khon Kaen University.Appreciation is extended to Faculty of Science， Khon Kaen University，Khon Kaen，Thailand.

Funding Statement:This research was financially supported by the young researcher development project of Khon Kaen University.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.