Skip to main content

ORIGINAL RESEARCH article

Front. Appl. Math. Stat., 01 June 2023
Sec. Mathematical Biology
This article is part of the Research Topic Using Mathematical Models to Understand, Assess, and Mitigate Vector-Borne Diseases View all 7 articles

Mathematical modeling for the control of fly-borne mastitis disease in cattle

  • Department of Pure and Applied Mathematics, Ladoke Akintola University of Technology, Ogbomoso, Nigeria

Several diseases cause losses in cattle farming, especially in the dairy industry among which mastitis disease (Bovine mastitis) is the leading cause of health and economic damages globally as it results in animals' ill health and reduced quality and quantity of milk produced by infected cows. Some mathematical studies have been conducted that focused on mastitis transmission from one udder-quarter to another in an infected cow, even though clinical studies established the cow–cow and flies–cow transmissions. The present study, therefore, proposed a mathematical model for the control of mastitis disease in cattle in the presence of flies as vectors. The formulated model was shown to have non-negative solutions in feasible regions for both cattle and flies populations. Furthermore, the model has a stable disease-free equilibrium if the sum of the effective reproduction numbers for cattle–cattle and fly–cattle transmissions (ch and ch) is less than unity, and there is a possibility of multiple endemic equilibria if otherwise. The numerical results indicated that the infectious populations can be reduced if the rates of the control parameters are increased, thereby curtailing or eradicating mastitis in the cattle population.

1. Introduction

The word mastitis is derived from two Greek words: masto (meaning mammary gland) and itis (meaning inflammation) [1]. Thus, mastitis is referred to as the inflammation of the mammary gland. Mastitis disease is caused by bacterial infection, trauma, or injury to the udder. It is the most common and expensive disease that affects dairy cattle throughout the world. Mastitis is caused by several different bacteria or other microorganisms (fungi, yeasts, and possibly viruses) that invade the udder, multiply and produce harmful substances there and then result in inflammation [2]. The most common microorganisms, causing mastitis, can be divided into two groups: (i) contagious bacteria (e.g., Staphylococcus aureus, Mycoplasma spp., and Streptococcus agalactiae) and (ii) environmental bacteria (e.g., Streptococcus uberis, Streptococcus dysgalactiae subsp. dysgalactiae, coliforms including Escherichia coli and Klebsiella spp.) [3].

Mastitis disease has well-recognized detrimental effects on animals' wellbeing and dairy farm profitability [4]. Mastitis can, directly, affect the mammary gland of the cow, leading to a significant reduction of the quantity and quality of the milk and, therefore, to reduced value of the production [5]. Despite decades of implementation of control strategies, mastitis continues to be one of the most significant and economically challenging problems of dairy cows [6, 7]. It reduces the productivity of the cows and the quality of milk causing enormous losses for breeders and, consequently, to the national income of the country [2]. According to the estimation of the National Mastitis Council (USA), mastitis costs more than $2 billion to dairy farmers annually [3, 7].

Mastitis is transmitted by direct contact between an infected and a susceptible animal which becomes infected by inhaling droplets disseminated in the environment. It can also be transmitted through an infected vector (such as a horn fly and, possibly, a housefly) [810]. The presentation of mastitis in a cow can be either a clinical or subclinical infection. In the subclinical type, the infection is asymptomatic and there are no visible changes in the appearance of the milk and/or the udder, milk production decreases by 10–20% with undesirable effect on its constituents and nutritional value rendering it of low quality and unfit for processing [11]. Subclinical (or acute) mastitis is the most common and economically harmful infection. Clinical mastitis, on the other hand, is characterized by heat, pain, swelling, and redness of the udder, along with reduced and an abnormal nature of milk yield. It is usually accompanied by a mild fever and animal depression. The affected quarter is sensitive to touch and painful to the animal. If acute mastitis is not attended, the inflammatory process persists for long, it gets converted into chronic mastitis which renders the milk-secreting tissue unable to produce any more milk. These changes are generally incurable and permanent [2]. The control of subclinical mastitis is more important than simply treating clinical cases because of the following reasons: The cows that have subclinical mastitis are reservoirs of organisms that lead to infection of other cows; and most clinical cases start as subclinical; thus, controlling subclinical mastitis is the best way to reduce the clinical cases [5]. Therefore, routine physical examination of the udder and diagnostic screening tests for early detection of mastitis and proper treatment of affected animals are of paramount importance to minimize losses due to subclinical and clinical mastitis [2].

Cows suffering from mastitis may recover spontaneously, but usually, drug therapy is required to maintain productivity [2]. Treatment of mastitis disease is achieved via systemic and intramammary administration of antibiotics, supportive fluid, and anti-inflammatory therapy [4]. Strict hygiene has to be maintained in cattle bedding to prevent transmission of infection during milking via contaminated milk, hands of the milker, and udder cloths (in the case of the milking machine) [2]. Culling of cows suffering from recurrent clinical mastitis (i.e., chronic mastitis) [2], and vector control via the application of insecticides (such as methoprene, synthetic pyrethroids, organophosphates, and abamectins) [9] are also good mastitis disease control measures.

Some mathematical studies have been conducted that focused on mastitis transmission from one udder-quarter to another in an infected cow [8, 1214], even though clinical studies established the cow–cow and flies–cow transmissions [4, 5, 7, 10]. Therefore, the present study proposed an epidemic mathematical model for the control of mastitis disease in cattle in the presence of flies as vectors. The goals were to analyze the dynamic transmission and to investigate the impacts of various control measures available for vector-borne mastitis. The rest of this study is presented with the following organization: Section 2 provided the detail of the model's description and a domain in which it is feasible epidemiologically and well-posed mathematically. Section 3 demonstrated the analyses of the equilibrium points, which include the fundamental reproduction number derivation and analyses of local and global stabilities of the equilibrium points and the sensitivity analysis. Section 4 contained the model's numerical simulations and graphical illustrations and a discussion of the results. Section 5 included the conclusion of the study.

2. Mathematical model

To study the control of fly-borne mastitis disease in cattle using a mathematical model, a mathematical model is hereby formulated using a system of ordinary differential equations as follows. Furthermore, the model was shown to be well-posed by establishing the positivity of solutions and the region where the model is epidemiologically feasible.

2.1. Model formulation

The cattle population and fly population are the interacting host and vector populations, respectively. The cattle population is subdivided into the subpopulations of susceptible [Sc(t)], exposed [E(t)], subclinical infective [A(t)], clinical infective [Ic(t)], chronic infective [C(t)], and recovered cattle [R(t)], so that the total cattle population at any time, t, is given as:

Nc(t)=Sc(t)+E(t)+A(t)+Ic(t)+C(t)+R(t).    (1)

The vector (flies) population is subdivided into the subpopulations of susceptible [Sv(t)] and infected flies [Iv(t)], and so that the total vector population at any time, t, is given as:

Nv(t)=Sv(t)+Iv(t).    (2)

The model takes into consideration cattle–cattle transmission, flies–cattle transmission, and cattle–flies transmission with transmission coefficients β, bp, and bq, respectively. Thus, the forces of infection for both populations are λc = [β(A(t) + aIc(t)) + bpIv(t)] and λv = bq(A(t) + Ic(t)). The interaction between the cattle and fly populations is, therefore, governed by the following proposed model, depicted by the schematic diagram in Figure 1.

FIGURE 1
www.frontiersin.org

Figure 1. Schematic diagram showing the dynamics of mastitis in cattle and fly populations.

2.2. Model's assumptions

Other assumptions on which the model's formulation is based are as follows:

  i. Both populations are assumed not to be constant since births, immigrations, emigrations, and deaths occur in the populations.

 ii. The model assumes a homogeneous mixing of individuals in both populations where all individuals have an equal likelihood of contracting mastitis if they come into effective contact with infectious individuals, and that transmission of the infection occurs with a mass action incidence.

iii. The probability of survival till the infectious state for cattle exposed to mastitis is less than or equal to unity. As a result, the exposed class of cattle is included in the model for the cattle population but ignored for the fly population since flies generally have a short lifecycle.

iv. The chronic infective cows are assumed not to transmit mastitis because they are weak and, therefore, are culled from the herd at a rate, d.

 v. Mastitis-induced deaths of cows only occur at the chronic stage of the infection at a rate δ. The flies also experience induced deaths, at a rate ϕ, due to the application of insecticides such as methoprene, synthetic pyrethroids, organophosphates, and abamectins [9].

vi. Individual cattle who recover from mastitis, spontaneously or as a result of treatment via systemic and intramammary administration of antibiotics, supportive fluid, and anti-inflammatory therapy [4], lose their immunity, and become susceptible since mastitis is not known to confer permanent immunity.

From the schematic diagram in Figure 1 together with the set of the model's assumptions, the model is formulated as a system of ordinary differential equations, using the mass action incidence rate, as given in Equation (3).

dSc(t)dt=Λc-[β(A(t)+aIc(t))+bpIv(t)]Sc(t)-μcSc(t)+κR(t)dE(t)dt=[β(A(t)+aIc(t))+bpIv(t)]Sc(t)-(μc+ρ)E(t)dA(t)dt=ερE(t)-(μc+α+γ)A(t)dIc(t)dt=(1-ε)ρE(t)+αA(t)-(μc+θ+γ+τ)Ic(t)dC(t)dt=θIc(t)-(μc+δ+d)C(t)dR(t)dt=γA(t)+(γ+τ)Ic(t)-(μc+κ)R(t)dSv(t)dt=Λv-bq(A(t)+Ic(t))Sv(t)-(μv+ϕ)Sv(t)dIv(t)dt=bq(A(t)+Ic(t))Sv(t)-(μv+ϕ)Iv(t);}    (3)

subject to the initial conditions:

Sc(0)=Sc0,E(0)=E0,A(0)=A0,Ic(0)=Ic0,C(0)=C0,  R(0)=R0,Sv(0)=Sv0,Iv(0)=Iv0.    (4)

The feasible region for the model is Ω = Ω1Ω2, where

Ω1={(Sc(t),E(t), A(t),Ic(t), C(t), R(t))ϵ+6:Sc(t) + E(t) + A(t) + Ic(t) + C(t) + R(t) = Nc(t)} and Ω2={(Sv(t), Iv(t))ϵ+2:Sv(t)+Iv(t)=Nc(t)} are the region of feasibility for cattle and fly populations, respectively.

The parameters of the model are described in Table 1.

TABLE 1
www.frontiersin.org

Table 1. Description of the model's parameters.

2.3. Positivity and boundedness of solutions of the model

Theorem 1: Suppose the initial conditions (4) of the state variables are non-negative in the region Ω = Ω1Ω2, then the solutions set {Sc(t), E(t), A(t), Ic(t), C(t), R(t), Nc(t)} ∈ Ω1 and {Sv(t), Iv(t), Nv(t)} ∈ Ω2 are non-negative in the region Ω for all time t ≥ 0.

Proof: It follows from the first equation of the model (3) that

dSc(t)dt-[μc+β(A(t)+aIc(t))+bpIv(t)]Sc(t).

Separating the variables and integrating both sides gives

     ln Sc(t)-[μc+β(A(t)+aIc(t))+bpIv(t)]dt     Sc(t)A1e-[μct+(β(A(t)+aIc(t))+bpIv(t))dt].             S(0)=S0A1=S0.      Sc(t)S0e-[μct+(β(A(t)+aIc(t))+bpIv(t))dt]Thus Sc(t)0,  provided that  S00.

Similarly, it can be shown that E(t), A(t), Ic(t), C(t), R(t), Sv(t) and Iv(t) are non-negative. Hence, all the state variables are non-negative for all time t ≥ 0, whenever the initial conditions are non-negative.

Theorem 2: Every solution in the region Ω = Ω1Ω2 are positively invariant with respect to the mastitis model (3) in cattle and fly populations.

Proof: Differentiating (1) and then substituting the first six equations of the model (3) gives,

dNc(t)dt=Λc-μcNc(t)-(δ+d)C(t)dNc(t)dtΛc-μcNc(t)dNc(t)dt+μcNc(t)ΛcN(0)=N0A2=N0-Λcμc.

Solving, using the integration factor, I.F= eμcdt , gives,

Nc(t)e-μct[Λceμctdt]Nc(t)e-μct[Λcμceμct+A2]Nc(t)Λcμc+A2e-μct.N(0)=N0A2=N0-Λcμc.Nc(t)Λcμc+(N0-Λcμc)e-μct.

As t0,Nc(t)N0 and As t,Nc(t)Λcμc

Thus  0Nc(t)Λcμc.

Moreover, differentiating (2) and substituting the last two equations of the model (3) and then following the above procedures give 0Nv(t)Λv(μv+ϕ).

These results imply that the mastitis model (3) has non-negative and bounded solutions in the region Ω = Ω1Ω2, ∀t ≥ 0. Hence, the proposed mastitis model is both epidemiologically feasible and mathematically well-posed, and so, it suffices to consider the dynamics of the model in the feasible region Ω.

3. Analysis of the model

In this section, the mastitis model (3) is analyzed as follows: the equilibrium points and the effective reproduction number were obtained. Moreover, the local and global stabilities for both equilibria were investigated, and a sensitivity analysis was carried out.

3.1. Equilibria and effective reproduction number

Equilibrium points are steady-state solutions of the model satisfying

dSc(t)dt=dE(t)dt=...=dIv(t)dt=0.    (5)

3.1.1. Mastitis-free equilibrium

In the absence of mastitis, i.e., when E(t) = A(t) = Ic(t) = C(t) = Iv(t) = 0, then both cattle and fly populations are in a state of mastitis-free equilibrium, E0, obtained by solving the model equation (3) subject to (5) and E(t) = A(t) = Ic(t) = C(t) = Iv(t) = 0 as follows:

E0=(Λcμc,0,0,0,0,0,Λv(μc+ϕ),0).    (6)

3.1.2. Effective reproduction number

The effective reproduction number, ℜc, is the actual average of secondary cases per primary case observed in a population with an infectious disease in the presence of control measures. Using the next generation operator, described by [15, 16],

c=ρ(FV-1)    (7)

where F is the Jacobian of the transmission matrix, obtained at E0 as

F=(0βΛcμcaβΛcμc0bpΛcμc0000000000000000bqΛv(μv+ϕ)bqΛv(μv+ϕ)00),    (8)

and V is the Jacobian of the transition matrix, obtained at E0 as

V=(μc+ρ0000-ερμc+α+γ000-(1-ε)ρ-αμc+θ+γ+τ0000-θμc+δ+d00000μv+ϕ),    (9)

Therefore, by (7), the effective reproduction number is the leading eigenvalue of the matrix FV−1, computed as

c=βΛc2μck1[ερk2+ak3(εραk2+(1-ε)ρ)]     +{βΛc2μck1[ερk2+ak3(εραk2+(1-ε)ρ)]}2+b2pqΛcΛvμc(μv+ϕ)k1k5[ερk2+1k3(εραk2+(1-ε)ρ)],    (10)

where k1 = μc + ρ, k2 = μc + α + γ, k3 = μc + θ + γ + τ, k4 = μc + δ+d, and k5 = μv + ϕ.

This can be summarized as

c=12ch+(12ch)2+cv,    (11)

where ch=βΛcμck1[ερk2+ak3(εραk2+(1-ε)ρ)] and cv=b2pqΛcΛvμc(μv+ϕ)k1k5[ερk2+1k3(εραk2+(1-ε)ρ)] represent the reproduction numbers for cattle–cattle transmission and fly–cattle transmission, respectively.

3.1.3. Mastitis-endemic equilibrium

In the presence of mastitis disease, i.e., when E(t) ≠ A(t) ≠ Ic(t) ≠ C(t) ≠ Iv(t) ≠ 0, then both cattle and fly populations are said to be in a state of mastitis-endemic equilibrium, Ee, obtained by solving (3) subject to (5) and E(t) ≠ A(t) ≠ Ic(t) ≠ C(t) ≠ Iv(t) ≠ 0 as

Ee=(Sc*(t),E*(t),A*(t),Ic*(t),C*(t),R*(t),Sv*(t),Iv*(t)),    (12)

where Sc*(t)=1μc{Λc+κk6[εργk2+γ+τk3(εραk2+(1-ε)ρ)-k1] E*(t)};

A*(t)=ερk2E*(t);Ic*(t)=1k3(εραk2+(1-ε)ρ)E*(t);C*(t)=θk3k4(εραk2+(1-ε)ρ)E*(t);R*(t)=1k6[εργk2+γ+τk3(εραk2+(1-ε)ρ)]E*(t);Sv*(t)=Λvk5-Iv*(t);Iv*(t)=Λv[ερk2+1k3(εραk2+(1-ε)ρ)]E*(t)1+ερk2+1k3(εραk2+(1-ε)ρ)E*(t);

and E*(t) is obtained as the roots of the quadratic equation

AE*2(t)+BE*(t)+C=0,    (13)

or equivalently, using the approach of [17]

A1(cv-1)E*2(t)+[k5+B1(ch-1)+B2(cv-1)]E*(t)                                      +k5[(ch+cv)-1]=0,    (14)

with positive constants A1, B1, and B2 defined as follows:

A1=κbqβk3k6[ερk2+1k3(εραk2+(1-ε)ρ)][εργk2      +γ+τk3(εραk2+(1-ε)ρ)];B1=bq[ερk2+1k3(εραk2+(1-ε)ρ)]; andB2=1k6[εργk2+γ+τk3(εραk2+(1-ε)ρ)];

so that the coefficients of (13) are given as

A=A1(cv-1),B=[k5+B1(ch-1)+B2(cv-1)],andC=k5{(ch+cv)-1}.

We will now explore different cases of the mastitis-endemic equilibrium with respect to conditions placed on A1,B1,B2,ch,cv and (ch+cv) as follows:

  i. If ch<1,cv<1 and (ch+cv)>1; then there exists a unique mastitis-endemic equilibrium;

 ii. If ch<1,cv<1,(ch+cv)=1 and ck5(B1+B2)<1, where c is a constant; then the model has two endemic equilibria, such that E*(t)=0 or E*(t)=k5B1[1-ck5(B1+B2) ];

iii. If ch<1,cv<1,(ch+cv)<1 and ck5(B1+B2)<1, where c is a constant; then the model has two endemic equilibria, one of which might coincide with the mastitis-free equilibrium;

iv. If ch>1,cv>1, and (ch+cv)>1; then there exists no endemic equilibria for the mastitis model;

 v. If ch>1,cv>1, and (ch+cv)>1; then the model has two distinct endemic equilibria.

3.2. Local asymptotic stability of the mastitis-free equilibrium

Theorem 3: If the sum (ch+cv)<1, and the quantities βΛc[ερ+(1-ε)ρ]μc[k1k2+k3k5+(k1+k2)(k3+k5)]<1 and βΛc(k3+k5)[ερ+(ερα(k3+k5)+(1-ε)ρ)]+b2pqΛcΛvερμc[k1k2(k3+k5)+(k1+k2)k3k5]<1, then the mastitis-free equilibrium is locally asymptotically stable. Otherwise, it is a saddle point that is unstable.

Proof: The stability of the mastitis-free equilibrium depends on the Jacobian matrix of the mastitis model (3), evaluated at E0 as

J|E0=(-(μc+λ)0-βΛcμc-aβΛcμc0κ0-bpΛcμc0-(k1+λ)βΛcμcaβΛcμc000bpΛcμc0ερ-(k2+λ)000000(1-ε)ρα-(k3+λ)0000000θ-(k4+λ)00000γγ+τ0-(k6+λ)0000-bqΛv(μv+ϕ)-bqΛv(μv+ϕ)00-(k5+λ)000bqΛv(μv+ϕ)bqΛv(μv+ϕ)000-(k5+λ)),    (15)

This has eigenvalues λ = −μc, λ = −k4, λ = −k5, λ = −k6; others being the roots of the quartic equation (16) as follows:

λ4+(k1+k2+k3+k5)λ3     +[k1k2+k3k5+(k1+k2)(k3+k5)][1-R1]λ2     +[k1k2(k3+k5)+(k1+k2)k3k5][1-R2]λ     +k1k2k3k5[1-(ch+cv)]=0.,    (16)

where

R1=βΛc[ερ+(1-ε)ρ]μc[k1k2+k3k5+(k1+k2)(k3+k5)],R2=βΛc(k3+k5)[ερ+(ερα(k3+k5)+(1-ε)ρ)]+b2pqΛcΛvερμc[k1k2(k3+k5)+(k1+k2)k3k5].

Now, if (ch+cv)<1, and the quantities R1 and R2 are both less than unity, then the quartic equation (16) has negative roots or complex roots with negative real parts. Hence, the mastitis-free equilibrium is locally asymptotically stable. Otherwise, it is unstable.

3.3. Local asymptotic stability of the mastitis-endemic equilibria

In Section 3.1.3, we established that the mastitis model has some cases of the existence of endemic equilibrium. Of interest are the cases with unique (one) and two endemic equilibria. When two endemic equilibria exist (especially in cases (ii) and (iii) of Section 3.1.3), then this might cause the mastitis model to exhibit a subcritical (backward) bifurcation in which case a stable endemic equilibrium co-exists with stable disease-free equilibrium for ℜc < 1. When this happens, then, we have the occurrence of bi-stability [18, 19] as a result of backward bifurcation and the condition ℜc < 1 is a necessary but not a sufficient condition to eradicate mastitis in both cattle and the fly populations [20]. The condition for the existence of this backward bifurcation is described in the result of Theorem 4. When only one equilibrium exists, it is unique and the stabilities of this unique endemic equilibrium are established in Theorem 4 and Theorem 6.

Theorem 4: If ℜc ≥ 1, then the mastitis-endemic equilibrium is locally asymptotically stable. Otherwise, it is unstable.

Proof: The center manifold theory approach, described by Castillo–Chavez and Song [21, 22], is employed to establish the local asymptotic stability of the mastitis-endemic equilibrium around ℜc = 1.

From (11), we have,

c=12ch+(12ch)2+cv.

Let β = β*, where β* is the chosen bifurcation parameter at ℜc = 1.

         1=12ch+(12ch)2+cv  1=ch+cv  1=β*Λcμck1[ερk2+ak3(εραk2+(1-ε)ρ)]           +b2pqΛcΛvμc(μv+ϕ)k1k5[ερk2+1k3(εραk2+(1-ε)ρ)]β*=μck1Λc[ερk2+ak3(εραk2+(1-ε)ρ)]-b2pqΛv(μv+ϕ)k1k5.,

so that the mastitis-free equilibrium is locally asymptotically stable if β<β*,i.e. if c<1, and unstable whenever β>β*,i.e. if c>1.

Now, let Sc(t) = x1, E(t) = x2, A(t) = x3, Ic(t) = x4, C(t) = x5, R(t) = x6, Sv(t) = x7, Iv(t) = x8, then the model in terms of these new state variables becomes

dx1dt=f1=Λc-[β(x3+ax4)+bpx8]x1-μcx1+κx6dx2dt=f2=[β(x3+ax4)+bpx8]x1-(μc+ρ)x2dx3dt=f3=ερx2-(μc+α+γ)x3dx4dt=f4=(1-ε)ρx2+αx3-(μc+θ+γ+τ)x4dx5dt=f5=θx4-(μc+δ+d)x5dx6dt=f6=γx3+(γ+τ)x4-(μc+κ)x6dx7dt=f7=Λv-bq(x3+x4)x7-(μv+ϕ)x7dx8dt=f8=bq(x3+x4)x7-(μv+ϕ)x8.}    (17)

The Jacobian matrix of the model (17) at the mastitis-free equilibrium, E0, when β=β*, such that c=1, is obtained as in (15) with the calculated eigenvalues. If ℜc = 1, this implies that ch+cv=1. Then, the quartic equation (16) has a zero eigenvalue λ = 0; with other eigenvalues being negative.

Applying the center manifold theory, let w=(w1,w2,w3,...,w8)T be the right eigenvector associated with the zero eigenvalues, λ = 0, where

w1=-w2μc{β*Λcμc[ερk2+ak3(εραk2+(1-ε)ρ)]      -κk6[εργk2+θ(γ+τ)k3k4(εραk2+(1-ε)ρ)]      +b2pqΛcΛvμc(μv+ϕ)k1k5[ερk2+1k3(εραk2+(1-ε)ρ)]};w2=w2>0 (afree eigenvector);w3=ερk2w2;w4=ak3(εραk2+(1-ε)ρ)w2;w5=θk3k4(εραk2+(1-ε)ρ)w2;w6=w2k6[εργk2+θ(γ+τ)k3k4(εραk2+(1-ε)ρ)];w7=-bqΛv(μv+ϕ)k5[ερk2+1k3(εραk2+(1-ε)ρ)]w2;w8=bqΛv(μv+ϕ)k5[ερk2+1k3(εραk2+(1-ε)ρ)]w2.

Similarly, the left eigenvector associated with the zero eigenvalues, λ = 0, is given by v = (v1, v2, v3, ..., v8), such that v. w = 1, where

v1=v5=v6=v7=0;v2=v2>0 (a free eigenvector);v3=Λcμck2(1-αk3)[β*+b2pqΛv(μv+ϕ)k5]v2;v4=Λcμck3[aβ*+b2pqΛv(μv+ϕ)k5]v2;v8=bpΛcμck5v2.

Computation of a and b

The coefficients a and b as defined by Castillo–Chavez and Song [21] are a=k,i,j=18vkwiwj2fk(E0,β*)xixj  and b=k,i=18vkwi2fk(E0,β*)xiβ*; and are algebraically computed as follows, considering only the non-zero components fk: k = 2, 3, 4, 8. With these definitions, the following was obtained:

a=2v2w1(w3β*+w4aβ*+w8bp)+2v8w7bq(w3+w4).

On substitution of necessary quantities and simplifying, the above reduces to

a=-v2w22k1Λc(ch+cv)2(1-Q);    (18)

where Q=κΛck1k6[εργk2+θ(γ+τ)k3k4(εραk2+(1-ε)ρ)] + cv2k1(ch+cv)2[ερk2+1k3(εραk2+(1-ε)ρ)].

Similarly, the coefficient b is obtained as

b=v2x1(w3+aw4)  =v2w2Λcμc[ερk2+1k3(εραk2+(1-ε)ρ)].    (19)

Thus b > 0, since v2, w2 > 0, and all parameters are positive.

Now, from (18), if the quantity Q < 1, then a < 1, b > 0, and the model exhibit a forward bifurcation around ℜc = 1; when β* < 0, then there exists a positive stable equilibrium. This implies that the disease-free is globally stable while the endemic equilibrium is locally asymptotically stable when ℜc > 1, but close to 1. This scenario is depicted in Figure 2A, and in this case, the condition ℜc < 1 is sufficient to eradicate mastitis disease from the populations.

FIGURE 2
www.frontiersin.org

Figure 2. Bifurcation curve at the threshold value, ℜc = 1: (A) forward bifurcation plot and (B) backward bifurcation plot.

If, however, Q > 1, then a > 1, b > 0, and so the model exhibits a backward bifurcation at ℜc = 1; when β* < 0, then x = 0 is unstable and there exists a negative and locally asymptotically stable endemic equilibrium. This implies that the disease-free is not globally stable while the endemic equilibrium is locally asymptotically stable, whenever ℜc > 1. This scenario is shown in Figure 2B, and in this case, the condition ℜc < 1 is not sufficient to eradicate mastitis in both populations. For mastitis disease to be curtailed or eradicated from the population, the condition c*<c<1 must hold.

3.4. Global asymptotic stability of the equilibria

In this section, the global asymptotic stability of the mastitis-free and mastitis-endemic equilibrium are investigated. It should be noted that whenever a stable endemic equilibrium co-exists with the disease-free equilibrium (for the existence of backward bifurcation), then the mastitis-free equilibrium cannot be globally stable and the endemic equilibrium also might not be globally stable [18, 23, 24]. Given this, the global stabilities are considered for the case with unique endemic equilibrium or when Q < 1 in Section 3.3 above and the model exhibits a forward bifurcation.

3.4.1. Global asymptotic stability of the mastitis-free equilibrium

To show the global stability of the mastitis-free equilibria, the two conditions for the global stability of disease-free as contained in [2527] for ℜc < 1 were established.

We begin by re-writing the model equation (3) in the form

dXdt=F(X,Y) and dYdt=G(X,Y),G(X,0)=0;

where X = (Sc(t), R(t), Sc(t)) and Y = (E(t), A(t), Ic(t), C(t), Iv(t)) are the non-infectious and the infectious states of the system; F(X, Y) and G(X, Y) are the right-hand sides of X. and Y., respectively, with E(t) = A(t) = Ic(t) = C(t) = Iv(t) = 0.

Now, from the model (3), the reduced system dXdt=F(X,0) is given as

dSc(t)dt=Λc-μcSc(t)+κR(t)dR(t)dt=-(μc+κ)R(t)dSv(t)dt=Λv-(μv+ϕ)Sv(t);    (20)

Let X*=(Sc*(t),R*(t),Sc*(t))=(Λcμc,0,Λv(μc+ϕ)) be the equilibrium of the reduced system (20) then the following result holds.

Theorem 5: The mastitis-free equilibrium, E0, of the model (3) is globally asymptotically stable if ℜc < 1; provided that the conditions (H1) and (H2) below are satisfied:

H1: For dXdt=F(X,0),X* is globally asymptotically stable;

H2: G(X, Y) = AY − Ĝ(X, Y), Ĝ(X, Y) ≥ 0 for (X, Y) ∈ Ω.

where A=DYG(X*,0) is a Metzler matrix (the off-diagonal elements of A are non-negative) and Ω is the region where the model is biologically feasible.

Proof: First, we show that X* is globally asymptotically stable by solving each equation in the reduced system (20) and then taking the limit as t → ∞.

Thus, solving the second equation in (20) subject to the initial condition R(0) = R0 gives

R(t)=R0e-(μc+κ)t.    (21)

As t → ∞, R(t) → 0, and independent of R0.

Moreover, from the first equation in (20), we have that

                          dSc(t)dt+μcSc(t)=Λc+κR(t)dSc(t)dt+μcSc(t)=Λc+κR0e-(μc+κ)t.

Solving this with its integrating factor, eμct, and subject to the initial condition Sc(0) = Sc0 gives

Sc(t)=Λcμc(1-e-μct)+Sc0e-μct+R0(e-(μc+κ)t-e-μct).    (22)

As t,Sc(t)Λcμc, and independent of Sc0.

From the last equation in (20), we have,

dSv(t)dt+(μv+ϕ)Sv(t)=Λv.

Solving this with its integrating factor, e(μv+ϕ)t, and subject to the initial condition Sv(0) = Sv0 gives

Sv(t)=Λvμv+ϕ+(Sv0-Λvμv+ϕ)e-(μv+ϕ)t.    (23)

As t,Sv(t)Λvμv+ϕ, and independent of Sv0.

This implies that the asymptotic dynamics of the system are independent of the initial conditions, i.e., every solution of the model with initial conditions in Ω approaches X*=E0, as t. Hence, the equilibrium X* is globally asymptotically stable for dXdt=F(X,0).

Next, it is left to show H2: G(X, Y) = AY − Ĝ(X, Y), Ĝ(X, Y) ≥ 0 for (X, Y) ∈ Ω.

From the model (3),

G(X,Y)=([β(A(t)+aIc(t))+bpIv(t)]Sc(t)-(μc+ρ)E(t)ερE(t)-(μc+α+γ)A(t)(1-ε)ρE(t)+αA(t)-(μc+θ+γ+τ)Ic(t)θIc(t)-(μc+δ+d)C(t)bq(A(t)+Ic(t))Sv(t)-(μv+ϕ)Iv(t));

and

A=DYG(X*,0)  =(-(μc+ρ)βΛcμcaβΛcμc0bpΛcμcερ-(μc+α+γ)000(1-ε)ρα-(μc+θ+γ+τ)0000θ-(μc+δ+d)00bqΛv(μv+ϕ)bqΛv(μv+ϕ)0-(μc+ϕ)),AY=(-(μc+ρ)βΛcμcaβΛcμc0bpΛcμcερ-(μc+α+γ)000(1-ε)ρα-(μc+θ+γ+τ)0000θ-(μc+δ+d)00bqΛv(μv+ϕ)bqΛv(μv+ϕ)0-(μc+ϕ))(E(t)A(t)Ic(t)C(t)Iv(t))   =(-(μc+ρ)E(t)+βΛcμcA(t)+aβΛcμcIc(t)+bpΛcμcIv(t)ερE(t)+-(μc+α+γ)A(t)(1-ε)ρE(t)+αA(t)-(μc+θ+γ+τ)Ic(t)θIc(t)-(μc+δ+d)C(t)bqΛv(μv+ϕ)A(t)+bqΛv(μv+ϕ)Ic(t)-(μc+ϕ)Iv(t)).

Then G(X, Y) = AY − Ĝ(X, Y) implies Ĝ(X, Y) = AYG(X, Y), so that

Ĝ(X,Y)=([β(A(t)+aIc(t))+bpIv(t)](Λcμc-Sc(t))000bq(A(t)+Ic(t))(bqΛv(μv+ϕ)-Sv(t)))          =(Ĝ1(X,Y)Ĝ2(X,Y)Ĝ3(X,Y)Ĝ4(X,Y)Ĝ5(X,Y)).    (24)

Obviously, Ĝ2(X, Y) = Ĝ3(X, Y) = Ĝ4(X, Y) = 0. Moreover, Ĝ1(X, Y) ≥ 0 and Ĝ5(X, Y) ≥ 0 since Sc(t)Λcμc and Sv(t)Λvμv+ϕ at any time t, as shown in H1.

Now, it has been established that: (i) For dXdt=F(X,0),X* is globally asymptotically stable; and

(ii) G(X, Y) = AY − Ĝ(X, Y), Ĝ(X, Y) ≥ 0 for (X, Y) ∈ Ω.

Hence, the mastitis-free equilibrium, E0, is globally asymptotically stable whenever ℜc < 1, and unstable otherwise.

3.4.2. Global asymptotic stability of the mastitis-endemic equilibrium

Theorem 6: If ℜc > 1, then the mastitis-endemic equilibrium, whenever it is unique, is globally asymptotically stable in Ω, provided that IcIc*;CC*;RR*;Sc*EScE*;Sv*Iv SvIv*.

Proof: The global stability of the mastitis-endemic equilibrium is analyzed by constructing a suitable Lyapunov function V(t) of Goh–Volterra type, following the approach of [2630], such that

V(t)=(Sc-Sc*-Sc*ln ScSc*)+(E-E*-E*ln EE*)       +K(A-A*-A*ln AA*)+(Ic-Ic*-Ic*ln IcIc*)       +(C-C*-C*(t)ln CC*)       +κR*γA*+(γ+τ)Ic*(R-R*-R*ln RR*)       +(Sv-Sv*-Sv*ln SvSv*)+(Iv-Iv*-Iv*ln IvIv*),    (25)

where K=[β(A*+aIc*)+bpIc*ερE*]Sc*.

The time derivative of V(t) along the solutions of the system (3) is obtained as

dV(t)dt=(1-Sc*Sc)dScdt+(1-E*E)dEdt+K(1-A*A)dAdt            +(1-Ic*Ic)dIcdt+(1-C*C)dCdt            +κR*γA*+(γ+τ)Ic*(1-R*R)dRdt            +(1-Sv*Sv)dSvdt+(1-Iv*Iv)dIvdt.

Putting the appropriate equations of the system (3) into the derivative of V(t) gives

dV(t)dt-=(1-Sc*Sc)(Λc-[β(A+aIc)+bpIv]Sc-μcSc+κR)                 +(1-E*E)([β(A+aIc)+bpIv]Sc-(μc+ρ)E)                 +K(1-A*A)(ερE-(μc+α+γ)A)                 +(1-Ic*Ic)((1-ε)ρE+αA-(μc+θ+γ+τ)Ic)                 +(1-C*C)(θIc-(μc+δ+d)C)                 +κR*γA*+(γ+τ)Ic*(1-R*R)(γA+(γ+τ)Ic-(μc+κ)R)                 +(1-Sv*Sv)(Λv-bq(A+Ic)Sv-(μv+ϕ)Sv)                 +(1-Iv*Iv)(bq(A+Ic)Sv-(μv+ϕ)Iv).    (26)

At the mastitis-endemic equilibria of the model, the following relationships hold from the model (3):

                          Λc=[β(A*+aIc*)+bpIv*]Sc*+μcSc*-κR*;               (μc+ρ)=[β(A*+aIc*)+bpIv*]Sc*E*(t);       (μc+α+γ)=ερE*A*;(μc+θ+γ+τ)=(1-ε)ρE*+αA*Ic*;       (μc+δ+d)=θIc*C*;               (μc+κ)=γA*+(γ+τ)Ic*R*;                          Λv=bq(A*+Ic*)Sv*+(μv+ϕ)Sv*;               (μv+ϕ)=bq(A*+Ic*)Sv*Iv*.

Using these relations in (26) yields

dV(t)dt=(1-Sc*Sc){([β(A*+aIc*)+bpIv*]Sc*-[β(A+aIc)+bpIv(t)]Sc)+μc(Sc*-Sc)+κ(R-R*)}+(1-E*E)([β(A+aIc)+bpIv]Sc-[β(A*+aIc*)+bpIv*]Sc*EE*)+Sc*[β(A*+aIc*)+bpIc*ερE*](1-A*A)(ερE-ερE*AA*)+(1-Ic*Ic)([(1-ε)ρE+αA]-[(1-ε)ρE*+αA*]IcIc*)+(1-C*C)(θIc-θIc*CC*)+κR*γA*+(γ+τ)Ic*(1-R*R)([γA+(γ+τ)Ic]-[γA*+(γ+τ)Ic*]RR*)+(1-Sv*Sv)(bq(A*+Ic*)Sv*-bq(A+Ic(t))Sv+(μv+ϕ)[Sv*-Sv])+(1-Iv*Iv)(bq(A+Ic(t))Sv-bq(A*+Ic*)Sv*IvIv*).

Expanding and simplifying this gives

dV(t)dt=Sc*[β(A*+aIc*)+bpIv*]{3-Sc*Sc-AA*-A*EAE*}+Sc[β(A+aIc)+bpIv(t)]{1-ScE*Sc*E}+μc(Sc*-Sc)(1-Sc*Sc)+κ(R-R*)(1-Sc*Sc)+[(1-ε)ρE+αA](1-Ic*Ic)+[(1-ε)ρE*+αA*](1-IcIc*)+θIc*{1-CC*-C*IcCIc*+IcIc*}+κR*{1-RR*-γA+(γ+τ)IcγA*+(γ+τ)Ic*+[γA+(γ+τ)IcγA*+(γ+τ)Ic*]R*R}+bqSv*(A*+Ic*){2-Sv*Sv-IvIv*}+bqSv(A+Ic(t)){1-SvIv*Sv*Iv}+Sv*(μv+ϕ){2-SvSv*-Sv*Sv}.

This simplifies to

dV(t)dt=Sc*[β(A*+aIc*)+bpIv*]{3-[Sc*Sc+AA*+A*EAE*-β(A+aIc)+bpIvβ(A*+aIc*)+bpIv*(1-ScE*Sc*E)]}+bqSv*(A*+Ic*){2-[Sv*Sv+IvIv*-bqSv(A+Ic(t))bqSv*(A*+Ic*)(1-SvIv*Sv*Iv)]}+μcSc*{2-[Sc*Sc+ScSc*]}+Sv*(μv+ϕ){2-[SvSv*-Sv*Sv]}+[(1-ε)ρE*+αA*]{1-[IcIc*-(1-ε)ρE+αA(1-ε)ρE*+αA*(1-Ic*Ic)]}+κR{1-[Sc*Sc+R*R+Sc*R*ScR]}+κR*{1-[RR*-γA+(γ+τ)IcγA*+(γ+τ)Ic*(1-R*R)]}+θIc*{1-[CC*-IcIc*(1-C*C)]}.}    (27)

Now, by this theorem's hypothesis and the arithmetic–geometric mean inequality (i.e., the fact that the arithmetic mean is greater than or equals the geometric mean) [18, 31], the following inequalities are true:

Sc*Sc+AA*+A*EAE*-β(A+aIc)+bpIvβ(A*+aIc*)+bpIv*(1-ScE*Sc*E)3;Sv*Sv+IvIv*-bqSv(A+Ic(t))bqSv*(A*+Ic*)(1-SvIv*Sv*Iv)2;Sc*Sc+ScSc*2;SvSv*-Sv*Sv2;IcIc*-(1-ε)ρE+αA(1-ε)ρE*+αA*(1-Ic*Ic)1;Sc*Sc+R*R+Sc*R*ScR1;RR*-γA+(γ+τ)IcγA*+(γ+τ)Ic*(1-R*R)1;CC*-IcIc*(1-C*C)1.}    (28)

It follows from (27) that dV(t)dt0, and that dV(t)dt=0 iff Sc=Sc*;E=E*;A=A*;Ic=Ic*C=C*;R=R*;Sv=Sv*;Iv=Iv*. Thus, the largest positively invariant set in the feasible region Ω is the unique endemic equilibrium. Hence, by LaSalle's invariance principle [32], the endemic equilibrium, whenever it is unique, is globally asymptotically stable in Ω if ℜc > 1. Otherwise, it is unstable.

3.5. Sensitivity analysis of the reproduction number

In this section, the robustness of the mastitis model to its parameter values was explored by performing the sensitivity analysis of the effective reproduction number, ℜc, of the model. By doing this, we have information on those parameters which significantly impact the mastitis model with respect to the effective reproduction number, ℜc. To achieve this, the approach of the normalized forward sensitivity index of a variable to a parameter was employed as used in [31, 3335]. By this approach, the normalized forward sensitivity index of ℜc that depends differentially on a parameter p is defined as:

χpc=cp×pc.    (29)

With this definition, we derive an analytical expression for the sensitivity of ℜc with respect to each component parameter of ℜc. For instance, the sensitivity index of ℜc with respect to the fly's rate of biting, b, is obtained as

χbc=cb×bc=+0.9367111925.    (30)

Similarly, the values for the sensitivity index of ℜc with respect to the other parameters for the given baseline parameter values are displayed in Table 2 and Figure 3.

TABLE 2
www.frontiersin.org

Table 2. Sensitivity indices of ℜc in terms of the model's parameter baseline values.

FIGURE 3
www.frontiersin.org

Figure 3. Graphical representation of the sensitivity indices of ℜc.

It can be read from the above index table in Table 2 and from Figure 3, that the most sensitive parameter is the fly's biting rate (b). Other parameters such as the spontaneous recovery rate (γ), rate of administration of antibiotics (τ), and rate of insecticide application (ϕ), among others, are also sensitive to the mastitis effective reproduction number, ℜc. To interpret these, χbc=+ 0.9367111925 implies that increasing (or decreasing) b by 10% increases (or decreases) ℜc by 9.367%; while χϕc=-0.3857046086 means that increasing (or decreasing) ϕ by 10% decreases (or increases) ℜc by 3.857%. The interpretation of the sensitivity indices of other parameters follows a similar manner as for b and ϕ.

4. Numerical simulations and discussion

This section presents the numerical simulation results for the mastitis model (3). The simulations were carried out with the aid of Maple 18.0 software to show the solution of the model equation, the global stability of the mastitis-free equilibrium, and the effects of sensitive parameters on mastitis transmission and control.

The parameter values used for the numerical simulations were sourced from existing literature, and where unavailable, values were assumed and fixed within a reasonable range as given in Table 1. The initial conditions are: Sc(0) = 500, E(0) = 100, A(0) = 30, Ic(0) = 50, C(0) = 20, R(0) = 10, for cattle population; and Sv(0) = 750, Iv(0) = 100, for fly population.

The results of the numerical simulations are given in Figures 411 to illustrate the system's behavior for different values of the mastitis model's parameters.

FIGURE 4
www.frontiersin.org

Figure 4. Plot of the global stability of the equilibria with various initial conditions: (A, C, E) are the global stability plots for the Mastitis-free equilibrium; (B, D, F) are the global stability plots for the Mastitis-endemic equilibrium.

The plots in Figure 4 illustrated the global stabilities of the mastitis-free and mastitis-endemic equilibria, E0 and Ee, and they agree with the results of the global stability analyses given in Theorems 5 and 6. These results imply that regardless of the initial value of the infective, mastitis can be wiped out from both cattle and fly populations, whenever ℜc < 1, since the plots in Figures 4A, C, E revealed that the solutions converge at the mastitis-free equilibrium points for subclinical and clinical infective cattle population and infected fly population. It should be noted that this result may not hold for the cases where the mastitis-endemic equilibria co-exist with the mastitis-free equilibrium at ℜc < 1, because at this point, the mastitis model exhibits backward bifurcation as shown in Figure 2B. However, the reverse is the case when ℜc > 1, since the plots in Figures 4B, D, F showed that the infective populations converge at the endemic equilibrium above the x-axis.

Figure 5 is the plots of all the solutions for both populations against time for the given parameter values in Table 1. The plot in Figure 5A indicates that the susceptible cattle declines as a result of the infection. This consequently increases the population of exposed cattle. The recovered cattle population increases due to the presence of various control measures. In Figure 5B, the plot shows the behavior of the fly population with time. The plot revealed that both susceptible and infected flies increase with time and the susceptible later decreases just above 3,000 flies.

FIGURE 5
www.frontiersin.org

Figure 5. Plot of the (A) cattle and (B) fly populations with time at the given model's parameters.

The effects of the fly's biting rate (b) on the dynamics of mastitis disease were shown in Figures 6AF. It was revealed by the plot in Figures 6A, E that increase in the fly's rate of biting (b) decreases the populations of the susceptible cattle and fly. This, in turn, caused an increase in the population of the subclinical infective, clinical infective, and chronic infective cattle, and the infected flies, as shown by the plots in Figures 6BD, F, respectively. In Figure 7, the effects of the cattle transmission rate, β, were investigated. The plots in the Figure showed that an increase in the rate of transmission via contact with the infectious cattle declines the population of susceptible cattle, as can be seen in the plot in Figure 7A, while the increase in this rate caused an increase in the populations of the subclinical infective, clinical infective, and chronic infective cattle as depicted by the plots in Figures 7BD, respectively.

FIGURE 6
www.frontiersin.org

Figure 6. Plot of the effect of fly's biting rate “b” on the (A) susceptible cattle, (B) subclinical infective cattle, (C) clinical infective cattle, (D) chronic infective cattle, (E) susceptible fly, and (F) infected fly populations.

FIGURE 7
www.frontiersin.org

Figure 7. Plot of the effect of cattle's transmission rate “β” on the (A) susceptible cattle, (B) subclinical infective cattle, (C) clinical infective cattle, and (D) chronic infective cattle populations.

From Figure 8, which showed the effect of the rate of spontaneous recovery, γ, on the cattle populations, it was observed that an increased rate of spontaneous recovery, γ, significantly reduced the infectious cattle, including the chronic infective ones, as shown by the plots in Figures 8AC. Consequent to this reduction, the recovered cattle population was increased. This is depicted by the plot in Figure 8D. Furthermore, we investigated the effect of boosting the recovery from mastitis disease via the administration of pharmaceutical/non-pharmaceutical therapy (such as the use of antibiotics), the rate of which is denoted by τ. It was then observed that the use of this control measure has a significant effect on decreasing mastitis disease, as can be seen from the plots in Figures 9A, B, that increase in τ decreased the populations of the clinical infective and chronic infective cattle but increased the population of the recovered cattle as in Figure 9C. This control measure also, interestingly, decreased the population of the infected flies, as shown in Figure 9D, thereby reducing the number of infected flies available to cause mastitis disease whenever biting occurs.

FIGURE 8
www.frontiersin.org

Figure 8. Plot of the effect of spontaneous recovery rate “γ” on the (A) subclinical infective cattle, (B) clinical infective cattle, (C) chronic infective cattle, and (D) recovered cattle populations.

FIGURE 9
www.frontiersin.org

Figure 9. Plot of the effect of rate of antibiotics administration “τ” on the (A) clinical infective cattle, (B) chronic infective cattle, (C) recovered cattle, and (D) infected fly populations.

The plot in Figure 10 showed the effect of the rate of culling “τ” on the dynamics of mastitis disease. Culling is a control measure that involved the total separation of infected, and weak cattle from the herd, to prevent the further spread of the disease. The result of this control measure as shown in Figure 10 revealed that an increase in the rate of culling will decrease the population of the chronically infective cattle and so they do not partake in the transmission of mastitis disease. This is one of the assumptions on which the proposed mastitis model in this study is based.

FIGURE 10
www.frontiersin.org

Figure 10. Plot of the effect of rate of culling “d” on the dynamics of mastitis disease.

This study also simulated the effect of the rate of application of insecticides “ϕ” on the populations of flies that transmit mastitis disease, and the results are plotted in Figure 11. These plots showed that the more the use of insecticides against such flies, the lesser the population of the flies, irrespective of their infection status. This implies that both susceptible and infected fly's populations are reduced whenever insecticides are applied as revealed by the plots in Figures 11A, B, respectively.

FIGURE 11
www.frontiersin.org

Figure 11. Plot of the effect of rate of application of insecticides “ϕ” on the (A) susceptible fly and (B) infected fly populations.

5. Conclusion

In this study, we proposed and formulated an epidemic mathematical model for the dynamics of mastitis disease in cattle (host) and fly (vector) populations with standard incidence rates. The intervention strategy focused on the treatment of infected cattle, culling of chronic infective cattle, and application of insecticides on the flies that bite cattle. The proposed model was shown to be epidemiologically feasible and mathematically well-posed by establishing that the solutions set is non-negative and bounded in a feasible region. The existence and stability of both disease-free and endemic equilibria were obtained, and these depended on the effective reproduction number, ℜc. If this threshold value, ℜc < 1, mastitis disease is under control, but if ℜc > 1, then mastitis will persist in the population.

We performed a sensitivity analysis on the effective reproduction number, ℜc, and the result revealed that the most sensitive parameter was the biting rate of flies. This can be significantly reduced by the use of insecticides, such as methoprene, synthetic pyrethroids, organophosphates, and abamectins [9] against the flies, as revealed by the results of the simulations. Other parameters such as the spontaneous recovery rate (γ), rate of administration of antibiotics (τ), and rate of insecticide application (ϕ), among others, are also sensitive to the mastitis effective reproduction number, ℜc. Therefore, mastitis control measures should be based on these parameters, among others, so that the spread of mastitis disease could be reduced in both cattle and fly populations.

The numerical simulations performed revealed that parameters like the biting rate of flies and cattle transmission rate both reduced healthy cows and increased mastitis infection in both populations. However, control measures such as infected cattle treatment (τ), culling of chronic infective cattle (d), and the use of insecticides against flies (ϕ), all have the effect of decreasing the spread of mastitis disease and decreasing the number of infective cattle and increasing healthy ones. It is, therefore, recommended that these control measures be adopted in the livestock and dairy industries, as this will increase the number of healthy cows available, and also improve the quantity and quality of milk produced for human consumption. Then meat and milk economy can also be greatly enhanced. Future study would be based on determining which of the control measures would yield an optimal result for controlling the disease. Control measures such as media advertisement and community awareness [3638] would be explored, among others, in the optimal control analysis of the mastitis model.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

MA contributed to the conception of the study, model formulation of the study, methodology and analysis, numerical simulations, and write-up of sections of the manuscript. TO contributed to the model formulation of the study, methodology and analysis, and numerical simulations. Both authors contributed to the manuscript revision, read, and approved the final version of the submitted manuscript.

Acknowledgments

The authors are thankful to the editor and the reviewers for their comments and suggestions that helped in improving this study.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Kumar P, Ojasvita Deora A, Sharma H, Sharma S, Mittal D, Bhanot V, et al. Bovine mastitis: a review. Middle East J Sci Res. (2020) 28:497–507. doi: 10.5829/idosi.mejsr.2020.497.507

CrossRef Full Text

2. Tewari A. Bovine mastitis: an important diary cattle disease, technical report. Dairyman. (2014) 2014:62–5.

Google Scholar

3. National Mastitis Council (NMC). Current Concept of Bovine mastitis. 5th ed. Verona: National Mastitis Council (2011).

Google Scholar

4. Ruegg PL. A 100-year review: mastitis detection, management and prevention. J Diary Sci. (2017) 100:10381–97. doi: 10.3168/jds.2017-13023

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Wattiaux MA. Mastitis: The Disease and its Transmission. Diary essentials, Babcock Institute for International Dairy Research and Development. Wisconsin: University of Wisconsin-Madison (2011).

6. Ruegg PL. New perspectives in udder health management. Vet Clin Food Anim. (2012) 28:149–63. doi: 10.1016/j.cvfa.2012.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Arsenopoulos K, Triantafillous E, Filioussis G, Papadopoulos E. Fly repellency using deltamethrin may reduce intramammary infections of dairy cows under intensive management. Comp Immunol Microbiol Infect Dis. (2018) 61:16–23. doi: 10.1016/j.cimid.2018.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Aligaz AA, Munganga JMW. Mathematical modelling of the transmission dynamics of contagious Bovine pleuropneumonia with vaccination and antibiotic treatment. J Appl Math. (2018) 2019:1–10. doi: 10.1155/2019/2490313

CrossRef Full Text | Google Scholar

9. Arnold M. Preventing Summer Mastitis in Heifers Begins with Horn Fly Control UK College of Agriculture, Food Environment, Cooperative Extension Service. (2019). Available online at: https://entomology.ca.uky.edu/files/recs0/ent12-dairy.pdf (accessed May 1, 2022).

10. Yeruham I, Braverman Y, Shipgel NY, Chizov-Ginzburg A, Saran A, Winkler M. Mastitis in dairy caused by Corynebacterium pseudotuberculosis and the feasibility of transmission by housefly I. Vet Quart. (1996) 18:87–9. doi: 10.1080/01652176.1996.9694623

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Abdel-Rady A, Sayed M. Epidemiological studies on subclinical mastitis in dairy cows in assiut governorate. Veterinary World. (2009) 2:373–80. doi: 10.5455/vetworld.2009.373-380

PubMed Abstract | CrossRef Full Text

12. Rachah A, Dalen G, Norstebo H, Reksen O, Barlow JW. Deterministic modelling of the transmission dynamics of intramammary infections. J Phys Conf Ser. (2018) 1132:012053. doi: 10.1088/1742-6596/1132/1/012053

CrossRef Full Text | Google Scholar

13. Dalen G, Rachah A, Norstebo H, Schukken YH, Gröhn YT, Barlow JW, et al. Transmission dynamics of intramammary infections caused by corynebacterium species. J Dairy Sci. (2018) 101:472. doi: 10.3168/jds.2017-13162

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Oltenacu PA, Natzke RP. Mathematical modelling of the mastitis infection process. J Dairy Sci. (1975) 15:515–21. doi: 10.3168/jds.S0022-0302(76)84233-6

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Van-den-Driessche P, Watmough J. Reproduction number and sub-threshold endemic equilibria for computational models of diseases transmission. Math Biosci. (2002) 180:29–48. doi: 10.1016/s0025-5574(02)00108-6

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and the computation of the basic reproduction ratio, Ro, in models for infectious diseases in heterogeneous populations. J Math Biol. (1990) 28:365. doi: 10.1007/BF00178324

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Ouifki R, Banasia J. Epidemiological models with quadratic equation for endemic equilibrium – a bifurcation atlas. Math Model Appl Sci. (2020) 43:10413–29. doi: 10.1002/mma.6389

CrossRef Full Text | Google Scholar

18. Martcheva M. An Introduction to Mathematical Epidemiology. New York, NY: Springer (2015). doi: 10.1007/978-1-4899-7612-3

CrossRef Full Text | Google Scholar

19. Arino J, Milliken E. Bistability in deterministic and stochastic sliar-type models with imperfect and waning vaccine protection. J Math Biol. (2022) 84:61. doi: 10.1007/s00285-022-01765-9

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Piazza N, Wang H. Bifurcation and sensitivity of immunity duration in an epidemic model. Int J Num Anal Modell Ser B. (2013) 4:179–202.

Google Scholar

21. Castillo-Chavez C, Song B. Dynamical models of tuberculosis and their applications. Math Biosci Eng. (2004) 1:361–404. doi: 10.3934/mbe.2004.1.361

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Buonomo B, Lacitignola D. On the backward bifurcation of a vaccination model with nonlinear incidence. Nonlinear Anal Modell Control. (2011) 16:30–46. doi: 10.15388/NA.16.1.14113

CrossRef Full Text | Google Scholar

23. Madubueze CE, Kimbir AR, Aboiyar T. Global stability of ebola disease model with contact tracing and quarantine. Applic Appl Math Int J. (2018) 13:382–403.

Google Scholar

24. Bonyah E, Khan MA, Islam S. A theoretical Model for Zika Virus Transmission. PLoS ONE. (2017) 12:e0185540. doi: 10.1371/journal.pone.0185540

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Castillo-Chavez C, Feng Z, Huang W. On the computation of R[[sb]]0[[/s]] and its role on global stability. Math Approach Emerg Reemerg Infect Dis Introd. (2002) 1:229–54. doi: 10.1007/978-1-4757-3667-0_13

CrossRef Full Text | Google Scholar

26. Adeyemi MO, Oluyo TO, Oladejo JK. Modelling the transmission and control dynamics of coronavirus disease with social distancing and contact tracing. Int J Innovat Sci Res Technol. (2020) 5:948–64.

PubMed Abstract | Google Scholar

27. Oladejo JK, Adeyemi MO, Oluyo TO. The impact of vaccine on the dynamical spread of Covid-19 pandemic. Int J Creat Res Thoughts. (2020) 8:4034–48.

Google Scholar

28. DeLeon CV. Construction of Lyapunov Functions for Classics SIS, SIR and SIRS Epidemic Model with Variable Population Size. Unite Academia de Mathematicas, Universidadn Autonoma de Guerrere, Mexico Facultad de Estudios Superiores Zarogoza, UNAM, Mexico (2009).

Google Scholar

29. Okuonghae D, Gumel AB, Ikhimwin BO, Iboi E. Mathematical assessment of the role of early latent infections and targeted control strategies of syphilis transmission dynamics. Acta Biotheor. (2018) 67:47–84. doi: 10.1007/s10441-018-9336-9

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Obabiyi OS, Olaniyi S. Global stability analysis of malaria transmission dynamics with vigilant compartment. Electron J Differ Equ. (2019) 2019:1–10.

Google Scholar

31. Akanni JO, Akinpelu FO, Olaniyi S, Oladipo AT, Ogunsola AW. Modelling financial crime population dynamics: optimal control and cost-effectiveness analysis. Int J Dyn Control. (2019) 8:531–44. doi: 10.1007/s40435-019-00572-3

CrossRef Full Text | Google Scholar

32. LaSalle JP. The Stability of Dynamical systems, Regional conference Series in Applied Mathematics. Philadelphia, PA: SIAM (1976).

Google Scholar

33. Oluyo TO, Adeyemi MO. Sensitivity analysis of Zika epidemic model. Int J Sci Eng Res. (2018) 9:587–601.

34. Chitnis NR, Hyman JM, Chusing JM. Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull Math Biol. (2008) 70:1272–96. doi: 10.1007/s11538-008-9299-0

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Onuorah MO, Akinwande NI, Nasir MO, Ojo MS. Sensitivity analysis of lassa fever model. Int J Math Stat. (2016) 4:30–49.

Google Scholar

36. Ghosh I, Tiwari PK, Samanta S, Elmojtaba IM, Al-Salti N, Chattopadyay J. A simple SI-type model for HIV/AIDS with media and self-imposed psychological fear. Math Biosci. (2018) 306:160–9. doi: 10.1016/j.mbs.2018.09.014

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Tiwari PK, Rai RK, Khajanchi S, Gupta RK, Misra AK. Dynamics of corona virus pandemic: effects of community awareness and global information campaigns. Eur Phys J Plus. (2021) 136:994. doi: 10.1140/epjp/s13360-021-01997-6

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Rai RK, Khajanchi S, Tiwari PK, Venturino E, Misra AK. Impact of social media advertisement on the transmission dynamics of COVID-19 pandemic in India. J Appl Math Comput. (2022) 68:19–44. doi: 10.1007/s12190-021-01507-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: mastitis, clinical and subclinical infective, treatment, culling, effective reproduction number, multiple equilibria, stability, simulations

Citation: Adeyemi MO and Oluyo TO (2023) Mathematical modeling for the control of fly-borne mastitis disease in cattle. Front. Appl. Math. Stat. 9:1171157. doi: 10.3389/fams.2023.1171157

Received: 21 February 2023; Accepted: 24 April 2023;
Published: 01 June 2023.

Edited by:

Meksianis Ndii, University of Nusa Cendana, Indonesia

Reviewed by:

Pankaj Tiwari, University of Kalyani, India
Hasan S. Panigoro, Universitas Negeri Gorontalo, Indonesia
Asep K. Supriatna, Padjadjaran University, Indonesia

Copyright © 2023 Adeyemi and Oluyo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Moses Olayemi Adeyemi, moadeyemi64@lautech.edu.ng

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.