- Department of Mathematics, College of Natural and Computational Sciences, Debre Berhan University, Debre Berhan, Ethiopia
In this study, we have proposed and analyzed a new COVID-19 and syphilis co-infection mathematical model with 10 distinct classes of the human population (COVID-19 protected, syphilis protected, susceptible, COVID-19 infected, COVID-19 isolated with treatment, syphilis asymptomatic infected, syphilis symptomatic infected, syphilis treated, COVID-19 and syphilis co-infected, and COVID-19 and syphilis treated) that describes COVID-19 and syphilis co-dynamics. We have calculated all the disease-free and endemic equilibrium points of single infection and co-infection models. The basic reproduction numbers of COVID-19, syphilis, and COVID-19 and syphilis co-infection models were determined. The results of the model analyses show that the COVID-19 and syphilis co-infection spread is under control whenever its basic reproduction number is less than unity. Moreover, whenever the co-infection basic reproduction number is greater than unity, COVID-19 and syphilis co-infection propagates throughout the community. The numerical simulations performed by MATLAB code using the ode45 solver justified the qualitative results of the proposed model. Moreover, both the qualitative and numerical analysis findings of the study have shown that protections and treatments have fundamental effects on COVID-19 and syphilis co-dynamic disease transmission prevention and control in the community.
1. Introduction
Communicable diseases are illnesses caused by pathogenic microbial agents such as bacteria, viruses, fungi, and parasites, which affect human beings throughout the world [1]. The novel coronavirus (COVID-19) infection is a lethal disease that has been a major global public health concern. The COVID-19 pandemic has affected various animals mostly infecting millions of human beings in different nations throughout the world [2–6]. It has been spreading mainly through sneezing, individuals interacting with each other in a certain time frame, or through coughing [7]. Although different species of animals are thought to be the source of COVID-19 transmission, bats have been shown to be coronavirus hosts [8]. Many nations throughout the world have started to practice various prevention and control strategies such as lockdown approach, quarantine, isolation, and closing schools [3, 9].
Syphilis is a major sexually transmitted disease and has been affecting millions of individuals both in low- and high-income countries of the world [10]. It is a chronic systemic disease caused by Treponema pallidum bacterium which is mainly transmitted through sex, blood contact, and mother-to-child during birth [4, 10–16]. Diagnosis, treatment, and using a condom are the basic control mechanisms of syphilis spreading in the community [10]. If left untreated, syphilis progresses through four stages: primary, secondary, latent, and tertiary [17–19]. The first three infection stages can transmit the disease to other susceptible groups of individuals, the transmission can occur via sexual contact, and in most cases, the tertiary stage is not transmissible through sexual contact [19]. It can be a cause of different cardiovascular and neurological diseases [17]. Approximately 90% of new syphilis substantial morbidity and mortality data are recorded in low-income countries around the world [11, 16]. Co-infection is an infection of an individual with two or more microorganisms' species [20, 21]. COVID-19 is an opportunistic infection for people with a weak immune system who were already infected by acute and chronic infections such as pneumonia, TB, and HIV/AIDS.
Mathematical modeling approach research done by scholars using a deterministic method [10, 14], a stochastic method [7, 22], or a fractional order method [23–32] has made a great contribution to linking the scientific approach with real-world physical situations and also for the decision-making process for solving real-world problems [33]. Different scholars have formulated and analyzed mathematical models on COVID-19 transmission [7, 8, 22, 24–26, 29, 30, 34–37], syphilis transmission [10, 17–19, 23], and other infectious diseases transmission [20, 21, 27, 28, 33, 38–40]; however, no one has done analysis on COVID-19 and syphilis co-infection transmission dynamics.
Oshinubi et al. [41] proposed and analyzed a new age-dependent compartmental model for COVID-19 transmission. The qualitative analysis of the model includes the non-negativity and boundedness of the model solutions in a given region, and the existence, uniqueness, and stability of the model solutions. Using parameter estimation from three different nations Kuwait, France, and Cameroon, they carried out numerical simulations and have shown the fundamental role of vaccination on COVID-19 transmission. Babaei et al. [34] proposed and examined a model for novel coronavirus transmission with Caputo's fractional order approach. The finding of the study shows that quarantine has a very fundamental role to control transmission. Iboi et al. [17] formulated and analyzed a new multi-stage syphilis model to examine the role of transitory immunity loss in the spreading process. The analysis shows that the disease-free and unique endemic equilibrium points are globally asymptotically stable when the corresponding basic reproduction number is less than unity and greater than unity, respectively. The results show that high treatment rates in the primary and secondary stages have a positive effect on the remaining stages of infection. Nwankwo et al. [38] formulated a mathematical model to examine the interaction between HIV/AIDS and syphilis pathogens with syphilis treatment on the co-infection of syphilis and HIV/AIDS where treatment or HIV is not accessible. High treatment in the primary stage has a fundamental role in reducing both single infections and co-infections in the population. Teklu et al. [42] formulated a six-compartmental COVID-19 transmission model to examine the impacts of intervention measures. The results show that protection, treatment, and vaccinations are fundamental to minimizing infection in the population.
Because different scholars have been mainly concerned with studying COVID-19 and syphilis single infections, no one has studied syphilis and COVID-19 co-infection using a mathematical model approach. Therefore, in this study, we are interested in filling the gap by formulating and analyzing syphilis and COVID-19 model intervention strategies.
The remaining part of the article is organized as follows. Section 2 presents COVID-19 and syphilis co-infection model construction. Section 3 describes the qualitative model analysis. Section 4 presents the numerical and sensitivity analyses. Section 5 presents the discussions and conclusions of the whole research study.
2. Model construction
We have considered COVID-19 and syphilis co-infection by separating the four syphilis infection stages (primary, secondary, latent, and tertiary) into two, the asymptomatic and symptomatic groups, and we have divided the population N(t) into 10 mutually exclusive states, which are described in Table 1 as follows:
Assumptions and definitions of basic terms:
➢ Co-infectious humans do not transmit both infections simultaneously.
➢ COVID-19 infection is transmitted to susceptible individuals from Ic and Ics infectious groups at the transmission rate as follows:
➢ Syphilis infection is transmitted to susceptible individuals from Ias, Iss, and Ics infectious groups at the force of infection rate as follows:
Using variable and parameter definitions given in Tables 1, 2, respectively, the flowchart of the COVID-19 and syphilis co-infection model is represented in Figure 1.
Figure 1. Flowchart of the model (3) with forces of infections λC and λs as in (1) and (2), respectively.
Using the flowchart represented in Figure 1, the corresponding system of differential equations of the complete co-infection model (3) is written as follows:
2.1. Qualitative properties of the model (3)
System (3) represents the human population; we want to prove that all the solutions of the model are non-negative and bounded, respectively, in the following region:
Theorem 1: Let Pc(0) > 0, S(0) > 0, Ps(0), Ic(0) > 0, Qc(0) > 0, Ias(0) > 0, Iss(0) > 0, Ts(0) > 0, Ics(0) > 0, T(0) > 0 be the initial solutions of the system (3), then Pc(t), Ps(t), S(t), Ic(t), Qc(t), Ias(t), Iss(t), Ts(t), Tcs(t), and T(t) are positive in the region for any time t > 0.
Proof: Let τ = sup{t > 0:Pc(t) > 0, S(t) > 0, Ps(t), Ic(t) > 0, Qc(t) > 0, Ias(t) > 0, Iss(t) > 0, Ts(t) > 0, Ics(t) > 0, T(t) > 0}.
Since Pc(t), Ps(t), S(t), Ic(t), Qc(t), Ias(t), Iss(t), Ts(t), Tcs(t), and T(t) are continuous, and we deduce that τ > 0. If τ = +∞, then positivity holds, but, if 0 < τ < +∞, Pc(τ) = 0 or Ps(τ) = 0 or S (τ) = 0 or Ic(τ) = 0 or Qc(τ) = 0 or Ias(τ) = 0 Iss(τ) = 0 or Ts(τ) = 0 or Tcs(τ) = 0 or T(τ) = 0.
From model (3) first equation, we do have
After some calculations of integration, we got
From model (3) second equation, we have
After some calculations of integration, we have
From model (3) third equation, we have
After some calculations, we have
, where
, and by the definition of τ we have Pc(t) > 0, Ps(t), Ts(t) > 0, S(τ) > 0, so that S(τ) ≠ 0.
Similarly, by proving the remaining state variable, we have
Ic(τ) > 0, hence Ic(τ) ≠ 0, Qc(τ) > 0 hence Qc(τ) ≠ 0, Ias(τ) > 0 hence Ias(τ) ≠ 0, Iss(τ) > 0 hence Iss(τ) ≠ 0,Ts(τ) > 0 hence Ts(τ) ≠ 0, Tcs(τ) > 0 hence Tcs(τ) ≠ 0, and T(τ) > 0 hence T(τ) ≠ 0.
Thus, τ is not finite, and hence = +∞, which means all the model solutions are non-negative.
Theorem 2: The model feasible region Ω stated in (4) is bounded in .
Proof: The total human being of the model (3) is as follows:
Differentiating both sides gives the following result
After some steps, we have , and hence, the model solutions with positive initial solutions are bounded in Ω.
3. Model analysis in qualitative approach
The complete COVID-19 and syphilis co-infection model (3) depends on the results of the two sub-models analysis.
3.1. COVID-19 mono-infection model analysis
From the complete model (3), we have the COVID-19 mono-infection model taking values Ps = Ias = Iss = Ts = Ics = T = 0 as follows:
with N1(t) = Pc(t) + S(t) + Ic(t) + Qc(t) as a total population and λc = β2Ic.
3.1.1. COVID-19 infection-free equilibrium
The COVID-19 infection-free equilibrium of the model (5) at Ic = 0 is .
3.1.2. COVID-19 mono-infection reproduction number
This mono-infection model has one infectious class, Ic, and we can obtain basic reproduction numbers without a method of the next-generation matrix as follows:
3.1.3. COVID-19 incidence equilibrium point
The COVID-19 incidence equilibrium point of the system (5) is , where
Theorem 3: The COVID-19 mono-infection model (5) has a unique COVID-19 incidence (endemic) equilibrium point whenever.
Proof: Using equation (1), we have the following:
Then the non-zero value of after a simple simplification is as follows:
, if and only if .
Hence, the COVID-19 mono-infection model (5) has a unique incidence equilibrium point if and only if .
Theorem 4: COVID-19 infection-free equilibrium point of the model (5) is locally asymptotically stable if; otherwise, it is unstable.
Proof: The Jacobean matrix of the model (5) at the COVID-19 infection-free equilibrium point is
Further, the characteristics equation after a certain calculation gives us as follows:
Thus, each eigenvalue of the Jacobian matrix is negative if implies the COVID-19 infection-free equilibrium point is locally asymptotically stable whenever
Theorem 5: The COVID-19 infection-free equilibrium point denoted by of the COVID-19 mono-infection model is globally stable if ; otherwise, it is unstable.
Proof: Take the representative Lyapunov function l(Ic ) = aIc, where ,
Thus, , if , and the equality holds if , and hence the COVID-19 infection-free equilibrium point is globally asymptotically stable if
Theorem 6: The COVID-19 incidence denoted by of the COVID-19 mono-infection model (5) is locally asymptotically stable whenever; otherwise, it is unstable.
Proof: The Jacobean of the system (5) at
From the Jacobean matrix, the characteristics equation, after simplification, gives as follows:
Then we do have the eigenvalues λ1 = −μ < 0 or λ2 = −(β + μ) < 0 or
Hence, all the coefficients of the characteristics equations are positive when thus, the COVID-19 incidence equilibrium point has local asymptotic stability when .
3.2. Analysis of syphilis sub-model
The syphilis sub-model is obtained from the system (3) by making Pc = Ic = Qc = Ics = T = 0 and is as follows:
With N2(t) = Ps(t) + S(t) + Ias(t) + Iss(t) + Ts(t), and λs = β1(Ias + ϕ2Iss ).
3.2.1. Syphilis infection-free equilibrium
The syphilis infection-free equilibrium point of the model (6) was obtained by making Ias = Iss = 0 and is given by
3.2.2. Syphilis sub-model reproduction number
The syphilis sub-model (6) has two infectious classes, which are Ias, and Iss, then applying the next-generation matrix method stated in [43, 44] to obtain the basic reproduction number of the system (6) by computing FV−1 as follows:
and
Then, we applied Mathematica coding; we have
Thus, the reproduction number of the syphilis sub-model (6) is given by .
3.2.3. Syphilis incidence equilibrium point of the system (6)
Making the model (6) equation to zero, we have the syphilis incidence equilibrium point given by , where
Theorem 7: The syphilis incidence equilibrium point of syphilis in the model (6) is unique whenever
Proof: From the syphilis infection rate, we have
Hence, the syphilis sub-model (6) has a unique incidence equilibrium if .
Theorem 8: Syphilis infection-free equilibrium point of the model (3) has local asymptotic stability if; otherwise, it is unstable.
Proof: Jacobean of the model (6) at the syphilis infection-free equilibrium point is as follows:
From the Jacobian matrix, the characteristics equation after simplification is as follows:
Then the eigenvalues are λ1 = −(π + μ) < 0 or λ2 = −μ < 0 or λ3 = −(δγ + μ) < 0 or .
where,
Applying Routh–Hurwitz criteria stated in [33], each eigenvalue of the matrix is negative whenever; thus, the syphilis infection-free equilibrium point has local asymptotic stability if .
Theorem 9: Syphilis infection-free equilibrium point of the model (6) has global stability if ; otherwise, it is unstable.
Proof: Let the Lyapunov representative function be given as l(Ias, Iss ) = aIas + bIss, where , .
Hence, the syphilis-free equilibrium point is globally stable if .
3.3. COVID-19 and syphilis co-infection model analysis
3.3.1. The model (3) disease-free equilibrium
Making all the equations of (3) zero with Ic = Ias = Iss = Ics = 0, the disease-free equilibrium point of (3) is given by .
3.3.2. The model (3) reproduction number
The COVID-19 and syphilis co-infection model (3) reproduction number denoted by is calculated using next-generation matrix criteria, as stated in [44]. Since we have four infectious groups, those are Ic, Ias, Iss, and Ics, and we have
and
Then applying Mathematica, we have got
and
Then, the corresponding eigenvalues are λ1 = 0 or λ2 = 0 or or where , .
Therefore, the COVID-19–syphilis complete model (3) reproduction number denoted by is given by
3.3.3. Model (3) disease-free equilibrium local stability
Theorem 10: The full-model (3) disease-free equilibrium point has local asymptotic stability if otherwise, it is unstable.
Proof: The Jacobian of the COVID-19 and syphilis co-infection model (3) at E0 is as follows:
where a , , , , , , , , , , , , , , , , , , , , .
By using square block matrix properties, we rewrite the above determinant as follows:
From this, we do have
Then, the eigenvalue of the full model is as follows:
Therefore, the co-infection model disease-free equilibrium point has local asymptotic stability if
3.3.4. The full-model endemic equilibrium and stabilities
The COVID-19 and syphilis co-infection model endemic equilibrium point is denoted by
. The analysis of the COVID-19-only mono-infection system (5) and the syphilis-only sub-model (6) shows that there is no endemic equilibrium point whenever and , respectively, which means there is no endemic equilibrium point if for the co-infection model (3). In other words, the COVID-19 and syphilis co-infection disease-free equilibrium point have global stability if
The explicit calculation of the co-infection model endemic equilibrium in terms of model parameters is tedious analytically; however, the model (3) endemic equilibriums correspond to
1. if is the syphilis-free (COVID-19 persistence) equilibrium point.
The analysis of the equilibrium is similar to the endemic equilibrium in the model (5).
2. , if is the COVID-19-free (syphilis persistence) equilibrium point. The analysis of the equilibrium is similar to the endemic equilibrium in Equation (6).
3. is the COVID-19 and syphilis co-existence persistence equilibrium point. It exists when each component of is positive whenever for this case, we have shown its stability in the numerical simulation part given in Section 4.
4. Sensitivity analysis and numerical simulations
In this section, we carried out the sensitivity analysis to examine the most sensitive parameters in the disease spreading and numerical simulations to verify the qualitative results of the mathematical model (3). Particularly, some numerical verification is considered to illustrate the qualitative analysis and results of the preceding sections. Here, we have taken some parameter values from literature and assumed some of the parameter values that are not from real data since there is a lack of mathematical modeling analysis literature which have studied the COVID-19 and syphilis co-infection transmission dynamics in the community. The fundamental problem of numerical analysis of a mathematical model is how to estimate parameters. Randomly choosing the values of parameters in the model in plausible intervals followed by sensitivity to the parameters is possible partially to overcome the limitations of parameters [41].
Here, the numerical simulation is used for checking the behaviors of the full-model (3) solutions and the effects of parameters in the transmission as well as the controlling of COVID-19 infection, syphilis infection, and co-infection of COVID-19 and syphilis. For numerical simulation purposes, we have applied MATLAB ode45 code with parameter values given in Table 3.
4.1. Analysis of sensitivity
Definition: The syphilis and COVID-19 co-infection model (3) normalized forward sensitivity index for its variable reproduction number is denoted by its derivative depends on a parameter ζ is defined by SEI [20, 21, 42].
The syphilis and COVID-19 co-infection model sensitivity index values justify the significance of different parameters in the single infections and co-infection spreading in the community. The parameter which has the highest magnitude of the sensitivity index value compared to other parameter index values is the most sensitive. Here, we have calculated the sensitivity index values in terms of the basic reproduction number . Using parameter values stated in Table 3, the sensitivity index values of the model (3) are calculated in Tables 4, 5.
Using parameter values in Table 3, we have computed and biologically, it means that syphilis infection, COVID-19 infection, and their co-infection are spreading in the population. The sensitivity index values stated in Table 4 explain that the recruitment rate Λ and the COVID-19 spreading rate β2 have a high direct impact on the COVID-19 basic reproduction . That means the recruitment rate and the COVID-19 transmission rates are the most sensitive parameters where stakeholders can control the transmission rate by applying prevention and control measures. Similarly, the COVID-19 protection portion τ1 and the quarantine with treatment rate ρ also have an indirect impact on the COVID-19 reproduction number .
Sensitivity indices stated in Table 5 explain that the recruitment rate Λ and syphilis spreading rate β1 have a high direct impact on the syphilis basic reproduction . That means the recruitment rate and syphilis transmission rates are the most sensitive parameters where stakeholders can control the transmission rate by applying prevention and control measures. Similarly, the syphilis protection portion τ3 and syphilis treatment rate γ have a high indirect effect on the syphilis reproduction number.
4.2. Results of numerical simulations
4.2.1. Behaviors of solutions of model (3) whenever
In the numerical simulation given in Figure 2, we observed that all the COVID-19 and syphilis co-infection model (3) solutions converge toward the disease-free equilibrium point whenever and with β1 = 0.3 and β2 = 0.08, respectively. At the co-infection disease-free equilibrium point, each solution curve of the model converges to zero while the susceptible group increases and then becomes constant, implying that the disease-free equilibrium point of the COVID-19 and syphilis co-infection model has global asymptotic stability if Biologically it means the COVID-19 and syphilis co-infection diseases have been eradicated from the community through time whenever .
4.2.2. Behaviors of the model solutions whenever
Figure 3 shows that all the COVID-19 and syphilis co-infection model (3) solutions converge toward the endemic equilibrium point whenever and with β1 = 8 and β2 = 11, respectively. After 10 years, the full-model solutions converge to the endemic equilibrium, while the susceptible population decreases and then remains constant means the COVID-19 and syphilis co-infection model endemic equilibrium point has local asymptotic stability if Biologically, it means that COVID-19 and syphilis co-infection disease spreads throughout the community under consideration.
4.2.3. Effects of protection measures on reproduction numbers
The numerical simulation represented by Figure 4 shows that when we maximize the COVID-19 rate of protection τ1, the reproduction number decreases, implying that the COVID-19 spreading rate decreases. Its biological meaning is that whenever the COVID-19 rate of protection τ1 > 0.7 the reproduction number , that is, the COVID-19 infection will be eradicated throughout the community.
Here, the numerical simulation represented by Figure 5 shows that whenever we maximize the syphilis protection rate τ3, the syphilis reproduction number decreases, implying that the syphilis spreading rate decreases. Whenever τ3 > 0.686 then , biologically, it means the syphilis infection eradicate from the community.
4.2.4. Impact of treatment on co-infected population
The numerical simulation given in Figure 6 shows that whenever the combined treatment rate ε of the COVID-19 virus and syphilis microorganism Treponema pallidum bacterium co-infected individuals Ics increases, the number of co-infected individuals decreases; that is, whenever the value of ε increases from 0.3 to 0.8, then the co-infected group Ics going down.
The numerical simulation given in Figure 7 shows that if the treatment rate ρ of COVID-19 increases, then the number of infections in the population decreases; that is, whenever ρ value increases from 0.2 to 0.8 then the infected group Ic decreases.
5. Discussions and conclusion
In this study, we have formulated and analyzed a new deterministic mathematical model for gaining insight into the effects of protections and treatments on the transmission dynamics of COVID-19 and syphilis co-infection. Both the positivity and boundedness of the complete model solutions have been discussed to show that the model is both mathematically and biologically meaningful. COVID-19 infection-free equilibrium point, COVID-19 incidence equilibrium point, and local and global stabilities of COVID-19 infection-free and COVID-19 incidence equilibrium points have been examined. Syphilis infection-free equilibrium point, syphilis incidence equilibrium point, and local and global stabilities of syphilis-free and syphilis incidence equilibrium points have been carried out. Using data stated in Table 3, we have carried out and discussed both sensitivity and numerical analyses of the full COVID-19 and syphilis co-infection model. From the analytical and numerical results, we observed that the model disease-free equilibrium points have global asymptotic stability when the basic reproduction numbers are less than unity. Biologically, this means that diseases die out in the community, with the full-model solutions converging to their endemic equilibrium point whenever their basic reproduction number is greater than unity, the reproduction numbers of both the COVID-19 infection and syphilis infection sub-models decreasing when the corresponding protection and treatment rates are maximized, and the numbers of co-infected individuals decreasing when the co-infection treatment rate is increased.
Based on the findings of this study, we recommend public health stakeholders concentrate on increasing both the COVID-19 and syphilis protection rates, as well as the syphilis treatment rate, the COVID-19 isolation with treatment rate, and the co-infection treatment rate, in order to reduce and possibly eradicate syphilis and COVID-19 co-infection transmission in the community. Finally, since no other COVID-19 and syphilis mathematical modeling approach literature has been formulated and analyzed, this study is not exhaustive. Interested researchers can extend this study in different manners, such as including syphilis mother-to-child transmission, COVID-19 vaccination as a new compartment, two-strain COVID-19 co-infection with syphilis, age structure for both infections, the four infection stages of syphilis (primary, secondary, latent, and tertiary), optimal control approach, stochastic method, fractional order method, and applying appropriate real population data.
Data availability statement
All relevant data is contained within the article: The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. The patients/participants provided their written informed consent to participate in this study.
Author contributions
Both authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Acknowledgments
The authors thank Mr. Sitotaw Eshete for his great Wi-Fi contribution that we have used in the 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. Panovska-Griffiths J. Can mathematical modelling solve the current COVID-19 crisis? BMC Public Health. (2020) 20:1–3. doi: 10.1186/s12889-020-08671-z
2. Adiga A, Dubhashi D, Lewis B, Marathe M, Venkatramanan S, Vullikanti A. Mathematical models for COVID-19 pandemic: a comparative analysis. J Indian Inst Sci. (2020) 100:793–807. doi: 10.1007/s41745-020-00200-6
3. Musuuza JS, Watson L, Parmasad V, Putman-Buehler N, Christensen L, Safdar N. Prevalence and outcomes of co-infection and superinfection with SARS-CoV-2 and other pathogens: a systematic review and meta-analysis. PLoS ONE. (2021) 16:e0251170. doi: 10.1371/journal.pone.0251170
4. Kumar P, Erturk VS, Murillo-Arcila M. A new fractional mathematical modelling of COVID-19 with the availability of vaccine. Results Phys. (2021) 24:104213. doi: 10.1016/j.rinp.2021.104213
5. Omame A, Sene N, Nometa I, Nwakanma CI, Nwafor EU, Iheonu NO, et al. Analysis of COVID-19 and comorbidity co-infection model with optimal control. Optimal Control Appl Methods. (2021) 42:1568–90. doi: 10.1002/oca.2748
6. Tuite AR, Fisman DN, Greer AL. Mathematical modelling of COVID-19 transmission and mitigation strategies in the population of Ontario, Canada. CMAJ. (2020) 192:E497–505. doi: 10.1503/cmaj.200476
7. Boukanjime B, Caraballo T, El Fatini M, El Khalifi M. Dynamics of a stochastic coronavirus (COVID-19) epidemic model with Markovian switching. Chaos Solitons Fractals. (2020) 141:110361. doi: 10.1016/j.chaos.2020.110361
8. Din A, Khan A, Baleanu D. Stationary distribution and extinction of stochastic coronavirus (COVID-19) epidemic model. Chaos Solitons Fractals. (2020) 139:110036. doi: 10.1016/j.chaos.2020.110036
9. Abdela SG, Berhanu AB, Ferede LM, van Griensven J. Essential healthcare services in the face of COVID-19 prevention: experiences from a referral hospital in Ethiopia. Am J Trop Med Hyg. (2020) 103:1198. doi: 10.4269/ajtmh.20-0464
10. Gumel AB, Lubuma J, Sharomi O, Terefe YA. Mathematics of a sex-structured model for syphilis transmission dynamics. Math Methods Appl Sci. (2018) 41:8488–513. doi: 10.1002/mma.4734
11. Momoh AA, Bala Y, Washachi DJ, Déthié D. Mathematical analysis and optimal control interventions for sex structured syphilis model with three stages of infection and loss of immunity. Adv Diff Equ. (2021) 2021:1–26. doi: 10.1186/s13662-021-03432-7
12. Peeling RW, Mabey D, Kamb ML, Chen XS, Radolf JD. “Benzaken.” AS Syphilis. Nat Rev Dis Prim. (2017) 3:17073. doi: 10.1038/nrdp.2017.73
13. Peeling RW, Hook III EW. The pathogenesis of syphilis: the Great Mimicker, revisited. J Pathol. (2006) 208:224–32. doi: 10.1002/path.1903
14. Saad-Roy CM, Shuai Z, Van den Driessche P. A mathematical model of syphilis transmission in an MSM population. Math Biosci. (2016) 277:59–70. doi: 10.1016/j.mbs.2016.03.017
15. Tessema B, Yismaw G, Kassu A, Amsalu A, Mulu A, Emmrich F, et al. Seroprevalence of HIV, HBV, HCV and syphilis infections among blood donors at Gondar University Teaching Hospital, Northwest Ethiopia: declining trends over a period of five years. BMC Infect Dis. (2010) 10:1–7. doi: 10.1186/1471-2334-10-111
16. World Health Organization. WHO Guidelines for the Treatment of Treponema pallidum (Syphilis). Geneva: WHO (2016).
17. Iboi E, Okuonghae D. Population dynamics of a mathematical model for syphilis. Appl Math Model. (2016) 40:3573–90. doi: 10.1016/j.apm.2015.09.090
18. Milner FA, Zhao R. A new mathematical model of syphilis. Math Model Nat Phenom. (2010) 5:96–108. doi: 10.1051/mmnp/20105605
19. Tuite AR, Fisman DN, Mishra S. Screen more or screen more often? Using mathematical models to inform syphilis control strategies. BMC Public Health. (2013) 13:1–9. doi: 10.1186/1471-2458-13-606
20. Teklu SW, Rao KP. HIV/AIDS-Pneumonia codynamics model analysis with vaccination and treatment. Comput Math Methods Med. (2022) 2022:3105734. doi: 10.1155/2022/3105734
21. Teklu SW, Mekonnen TT. HIV/AIDS-pneumonia coinfection model with treatment at each infection stage: mathematical analysis and numerical simulation. J Appl Math. (2021) 2021:5444605. doi: 10.1155/2021/5444605
22. Danane J, Allali K, Hammouch Z, Nisar KS. Mathematical analysis and simulation of a stochastic COVID-19 Lévy jump model with isolation strategy. Results Phys. (2021) 23:103994. doi: 10.1016/j.rinp.2021.103994
23. Bonyah E, Chukwu CW, Juga ML. Modeling fractional order dynamics of syphilis via Mittag-Leffler law. medRxiv. (2021). doi: 10.3934/math.2021485
24. Nabi KN, Abboubakar H, Kumar P. Forecasting of COVID-19 pandemic: from integer derivatives to fractional derivatives. Chaos Solitons Fractals. (2020) 141:110283. doi: 10.1016/j.chaos.2020.110283
25. Nabi KN, Kumar P, Erturk VS. Projections and fractional dynamics of COVID-19 with optimal control strategies. Chaos Solitons Fractals. (2021) 145:110689. doi: 10.1016/j.chaos.2021.110689
26. Kumar P, Erturk VS, Murillo-Arcila M, Banerjee R, Manickam A. A case study of 2019-nCOV cases in Argentina with the real data based on daily cases from March 03, 2020 to March 29, 2021 using classical and fractional derivatives. Adv Diff Equ. (2021) 2021:1–21. doi: 10.1186/s13662-021-03499-2
27. Etemad S, Avci I, Kumar P, Baleanu D, Rezapour S. Some novel mathematical analysis on the fractal-fractional model of the AH1N1/09 virus and its generalized Caputo-type version. Chaos Solitons Fractals. (2022) 162:112511. doi: 10.1016/j.chaos.2022.112511
28. Zarin R, Khan A, Kumar P. Fractional-order dynamics of Chagas-HIV epidemic model with different fractional operators. AIMS Math. (2022) 7:18897–924. doi: 10.3934/math.20221041
29. Gao W, Veeresha P, Baskonus HM, Prakasha DG, Kumar P. A new study of unreported cases of 2019-nCOV epidemic outbreaks. Chaos Solitons Fractals. (2020) 138:109929. doi: 10.1016/j.chaos.2020.109929
30. Kumar P, Erturk VS, Abboubakar H, Nisar KS. Prediction studies of the epidemic peak of coronavirus disease in Brazil via new generalised Caputo type fractional derivatives. Alexandr Eng J. (2021) 60:3189–204. doi: 10.1016/j.aej.2021.01.032
31. Vellappandi M, Kumar P, Govindaraj V. Role of fractional derivatives in the mathematical modeling of the transmission of Chlamydia in the United States from 1989 to 2019. Nonlinear Dyn. (2022) 2022:1–15. doi: 10.1007/s11071-022-08073-3
32. Abbas S, Tyagi S, Kumar P, Ertürk VS, Momani S. Stability and bifurcation analysis of a fractional-order model of cell-to-cell spread of HIV-1 with a discrete time delay. Math Methods Appl Sci. (2022) 45:7081–95. doi: 10.1002/mma.8226
33. Teklu SW, Terefe BB. Mathematical modeling analysis on the dynamics of university students animosity towards mathematics with optimal control theory. Sci Rep. (2022) 12:1–19. doi: 10.1038/s41598-022-15376-3
34. Babaei A, Ahmadi M, Jafari H, Liya A. A mathematical model to examine the effect of quarantine on the spread of coronavirus. Chaos Solitons Fractals. (2021) 142:110418. doi: 10.1016/j.chaos.2020.110418
35. Basnarkov L. SEAIR. Epidemic spreading model of COVID-19. Chaos Solitons Fractals. (2021) 142:110394.
36. Mohamadou Y, Halidou A, Kapen PT. A review of mathematical modeling, artificial intelligence and datasets used in the study, prediction and management of COVID-19. Appl Intell. (2020) 50:3913–25. doi: 10.1007/s10489-020-01770-9
37. Ndaïrou F, Area I, Nieto JJ, Torres DF. Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan. Chaos Solitons Fractals. (2020) 135:109846. doi: 10.1016/j.chaos.2020.109846
38. Nwankwo A, Okuonghae D. Mathematical analysis of the transmission dynamics of HIV syphilis co-infection in the presence of treatment for syphilis. Bull Math Biol. (2018) 80:437–92. doi: 10.1007/s11538-017-0384-0
39. Teklu SW, Kotola BS. The impact of protection measures and treatment on pneumonia infection: a mathematical model analysis supported by numerical simulation. bioRxiv. (2022). doi: 10.1101/2022.02.21.481255
40. Vellappandi M, Kumar P, Govindaraj V. Role of vaccination, the release of competitor snails, chlorination of water, and treatment controls on the transmission of bovine schistosomiasis disease: a mathematical study. Phys Script. (2022) 97:074006. doi: 10.1088/1402-4896/ac7421
41. Oshinubi K, Buhamra SS, Al-Kandari NM, Waku J, Rachdi M, Demongeot J. Age dependent epidemic modeling of COVID-19 outbreak in Kuwait, France, and Cameroon. In: Send P, editor. Healthcare (Vol. 10, No. 3). MDPI (2022), p. 482.
42. Teklu SW. Mathematical analysis of the transmission dynamics of COVID-19 infection in the presence of intervention strategies. J Biol Dyn. (2022) 16:640–64. doi: 10.1080/17513758.2022.2111469
43. Castillo-Chavez C, Feng Z, Huang W. On the computation of ro and its role on. In: Castillo-Chavez PC, Blower S, Driessche P, Kirschner D, Yakubu A-A, editors. Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction. vol 1, 229 (2002). doi: 10.1007/978-1-4757-3667-0_13
44. Van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math Biosci. (2002) 180:29–48. doi: 10.1016/S0025-5564(02)00108-6
Keywords: syphilis, COVID-19, COVID-19 and syphilis co-infection, protection, numerical simulation
Citation: Teklu SW and Terefe BB (2023) COVID-19 and syphilis co-dynamic analysis using mathematical modeling approach. Front. Appl. Math. Stat. 8:1101029. doi: 10.3389/fams.2022.1101029
Received: 17 November 2022; Accepted: 30 December 2022;
Published: 01 February 2023.
Edited by:
Yasser Aboelkassem, University of Michigan–Flint, United StatesReviewed by:
Vedat Erturk, Ondokuz Mayıs University, TürkiyePushpendra Kumar, Central University of Punjab, India
Copyright © 2023 Teklu and Terefe. 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: Shewafera Wondimagegnhu Teklu, luelzedo2008@gmail.com; shewaferaw@dbu.edu.et