Research article Special Issues

Mathematical modeling of the influence of cultural practices on cholera infections in Cameroon


  • The Far North Region of Cameroon, a high risk cholera endemic region, has been experiencing serious and recurrent cholera outbreaks in recent years. Cholera outbreaks in this region are associated with cultural practices (traditional and religious beliefs). In this paper, we introduce a mathematical model of the influence of cultural practices on the dynamics of cholera in the Far North Region. Our model is an SEIR type model with a pathogen class and multiple susceptible classes based on traditional and religious beliefs. Using daily reported cholera cases from three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region from June 25, 2019 to August 16, 2019, we estimate parameter values of our model and use Akaike information criterion (AIC) to demonstrate that our model gives a good fit for our data on cholera cases. We use sensitivity analysis to study the impact of each model parameter on the threshold parameter (control reproduction number), Rc, and the number of model predicted cholera cases. Finally, we investigate the effect of cultural practices on the number of cholera cases in the region.

    Citation: Eric Che, Eric Numfor, Suzanne Lenhart, Abdul-Aziz Yakubu. Mathematical modeling of the influence of cultural practices on cholera infections in Cameroon[J]. Mathematical Biosciences and Engineering, 2021, 18(6): 8374-8391. doi: 10.3934/mbe.2021415

    Related Papers:

    [1] Chayu Yang, Jin Wang . A cholera transmission model incorporating the impact of medical resources. Mathematical Biosciences and Engineering, 2019, 16(5): 5226-5246. doi: 10.3934/mbe.2019261
    [2] Fred Brauer, Zhisheng Shuai, P. van den Driessche . Dynamics of an age-of-infection cholera model. Mathematical Biosciences and Engineering, 2013, 10(5&6): 1335-1349. doi: 10.3934/mbe.2013.10.1335
    [3] Jiazhe Lin, Rui Xu, Xiaohong Tian . Transmission dynamics of cholera with hyperinfectious and hypoinfectious vibrios: mathematical modelling and control strategies. Mathematical Biosciences and Engineering, 2019, 16(5): 4339-4358. doi: 10.3934/mbe.2019216
    [4] Jinliang Wang, Ran Zhang, Toshikazu Kuniya . A note on dynamics of an age-of-infection cholera model. Mathematical Biosciences and Engineering, 2016, 13(1): 227-247. doi: 10.3934/mbe.2016.13.227
    [5] Chenwei Song, Rui Xu, Ning Bai, Xiaohong Tian, Jiazhe Lin . Global dynamics and optimal control of a cholera transmission model with vaccination strategy and multiple pathways. Mathematical Biosciences and Engineering, 2020, 17(4): 4210-4224. doi: 10.3934/mbe.2020233
    [6] Conrad Ratchford, Jin Wang . Modeling cholera dynamics at multiple scales: environmental evolution, between-host transmission, and within-host interaction. Mathematical Biosciences and Engineering, 2019, 16(2): 782-812. doi: 10.3934/mbe.2019037
    [7] Shu Liao, Jin Wang . Stability analysis and application of a mathematical cholera model. Mathematical Biosciences and Engineering, 2011, 8(3): 733-752. doi: 10.3934/mbe.2011.8.733
    [8] Beryl Musundi . An immuno-epidemiological model linking between-host and within-host dynamics of cholera. Mathematical Biosciences and Engineering, 2023, 20(9): 16015-16032. doi: 10.3934/mbe.2023714
    [9] Conrad Ratchford, Jin Wang . Multi-scale modeling of cholera dynamics in a spatially heterogeneous environment. Mathematical Biosciences and Engineering, 2020, 17(2): 948-974. doi: 10.3934/mbe.2020051
    [10] Christopher Botelho, Jude Dzevela Kong, Mentor Ali Ber Lucien, Zhisheng Shuai, Hao Wang . A mathematical model for Vibrio-phage interactions. Mathematical Biosciences and Engineering, 2021, 18(3): 2688-2712. doi: 10.3934/mbe.2021137
  • The Far North Region of Cameroon, a high risk cholera endemic region, has been experiencing serious and recurrent cholera outbreaks in recent years. Cholera outbreaks in this region are associated with cultural practices (traditional and religious beliefs). In this paper, we introduce a mathematical model of the influence of cultural practices on the dynamics of cholera in the Far North Region. Our model is an SEIR type model with a pathogen class and multiple susceptible classes based on traditional and religious beliefs. Using daily reported cholera cases from three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region from June 25, 2019 to August 16, 2019, we estimate parameter values of our model and use Akaike information criterion (AIC) to demonstrate that our model gives a good fit for our data on cholera cases. We use sensitivity analysis to study the impact of each model parameter on the threshold parameter (control reproduction number), Rc, and the number of model predicted cholera cases. Finally, we investigate the effect of cultural practices on the number of cholera cases in the region.



    According to the World Health Organization (WHO), cholera is an acute intestinal infection caused by ingestion of food or water contaminated with the bacterium Vibrio cholerae [1]. In many parts of the world, cholera remains a significant threat to public health [2]. It continues to devastate impoverished populations with limited access to medication, clean water and proper sanitation amenities [3]. Annually, about 1.3–4.0 million cholera cases and 21,000–143,000 cholera-induced deaths are reported worldwide [1,4]. In recent years, there has been several reports of major cholera outbreaks. For example, in Yemen from 2016–2017 cholera outbreak led to more than one million cases. Major outbreaks were also reported in Haiti (2010–2012), Sierra Leone (2012), Ghana (2011), Nigeria (2010), Vietnam (2009), Zimbabwe (2008), and India (2007). Cameroon, a cholera endemic African country experienced large cholera outbreaks in 1991, 1996, 1998, 2004, 2010 and 2011 [5]. In 2017, 1,227,391 cases and 5654 cholera-related deaths were reported from 34 countries [6].

    In recent years, the Far North Region of Cameroon has been experiencing serious recurrent cholera outbreaks, and these outbreaks are linked to cultural practices (traditional and religious beliefs) of the people of this region. The impact of traditional and religious practices on the spread of infectious diseases is a known fact in the literature. For example, in 2015, A. Manguvo and B. Mafuvadze in their study demonstrated that the transmission of Ebola in West Africa is linked to traditional and religious practices [7]. In 2017, Ngwa et al. [8] carried out a field study that demonstrate that cultural practices in the Far North region influence cholera transmission. The high cholera attack rate in mountainous areas was attributed to cultural practices such as mountain burial. They also noted that many people in the Far North had limited scientific knowledge about the cholera disease and its transmission. Rich cultural practices are common in villages and urban areas of the Far North Region of Cameroon. For example, the tradition of communal eating from the same plate and drinking from the same cup is common place [9]. In addition to cultural practices, other factors that contribute to cholera cases include the lack of clean drinking water and poor health facilities.

    Transmission of cholera can be indirect, from a water source containing the bacterium to susceptible humans, or direct, from infectious humans to susceptible humans. Several continuous-time mathematical ODE models have been developed in an effort to gain a deeper understanding of cholera transmission dynamics, see for example [10,11,12,13,14] and the citations therein. In a recent paper, Che et al. introduced a continuous-time ODE low-high risk structured cholera model and used it to capture the annual reported cholera infections in Cameroon from 1987 to 2017 [15]. Furthermore, using the fitted risk structured cholera model, they studied the impact of three intervention strategies (vaccination, treatment and improved sanitation) on the number of cholera infections in Cameroon from 2004 to 2022.

    Education can also serve as a disease intervention strategy. For example, Al-arydah et al. in [16] formulated and analyzed a mathematical model that includes two essential and affordable control measures: water chlorination and education. Cholera education includes advising people with symptoms to seek medical care promptly, and improving sanitation and hygienic practices. In the Far North Region of Cameroon, even with education, some people still maintain their traditional and religious practices which are sometimes in conflict with proper cholera hygiene [8].

    In this paper, we formulate a mathematical model of the influence of cultural practices on cholera infections in the Far North Region of Cameroon. We compute the control reproduction number, Rc, and using cumulative reported cholera cases for June 25–August 16, 2019 from three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region (from the Cameroon Ministry of Public Health) [17], we estimate model parameter values and demonstrate that our model captures the dynamics of the cholera infections. Also, we use sensitivity analysis to study the impact of each model parameter on the threshold parameter, Rc, and the number of model predicted cholera cases.

    The rest of the paper is organized as follows: In Section 2, we formulate our mathematical model with different susceptible classes representing cultural practices and with vaccination, and compute the control reproduction number Rc. Also, we study the reduced model without cultural practices and vaccination. Using the cumulative cholera cases from three health districts in the Far North Region, we fit the parameters to the cumulative cases. Using the Akaike information criterion (AIC), we illustrate that in the Far North Region of Cameroon, our model with cultural practices and with vaccination (in comparison with the model without cultural practices and vaccination) gives a better fit for the cumulative cholera cases. In Section 3, we use sensitivity analysis to study the impact of each model parameter on the threshold parameter, Rc, and the number of model predicted cholera cases. The discussion of our results and conclusions are presented in Section 4.

    Here, we formulate a model that captures the influence of cultural practices on the dynamics of cholera in the Far North Region, a high risk cholera region [15] of Cameroon. Our model is based on the schematic diagram shown in Figure 1. The list of the model variables and their descriptions are in Table 1, while the parameters are in Table 3. In our model, S(t) represents susceptible humans, E(t) represents exposed humans, I(t) represents infectious humans and R(t) represents recovered individuals at time t0. Also, B(t) represents the concentration of pathogen in the contaminated water at time t. To investigate the impact of cultural practices on the dynamics of cholera infections in the Far North Region, we partition the susceptible human subpopulation S into three classes: those who are not careful in practicing safe cholera intervention measures and not willing to be vaccinated as a result of their cultural beliefs (Sr), those who are careful in practicing safe cholera intervention measures but are not willing to be vaccinated (Sc) and those who are careful in practicing safe cholera intervention measures and are willing to be vaccinated (Sv), so that S=Sr+Sc+Sv. The size of the total human population at time t is

    N(t)=Sr(t)+Sc(t)+Sv(t)+E(t)+I(t)+R(t).
    Figure 1.  Flow diagram for the cholera model. Solid lines represent flow between compartments, the straight dashed line represents the infected class shedding pathogen into the environment, and the curved dashed line represents the source of new infections resulting from susceptible individuals interacting with the pathogen. All other arrows represent natural deaths, birth or growth of pathogen.
    Table 1.  Model variables, descriptions and units.
    Variable Description Unit
    S Susceptible individuals individuals
    Sr Susceptible individuals who are not careful and not willing to be vaccinated individuals
    Sc Susceptible individuals who are careful but not willing to be vaccinated individuals
    Sv Susceptible individuals who are careful and willing to be vaccinated individuals
    E Exposed individuals individuals
    I Infectious individuals individuals
    R Recovered individuals individuals
    B Pathogen concentration in water environment cells ml1
    N Total population individuals

     | Show Table
    DownLoad: CSV

    At time t0, individuals in the Sr, Sc and Sv classes are recruited (by birth or immigration) at rates αβ, αcβ, αvβ respectively, where β is the constant recruitment rate and α,αc and αv are proportional constants in (0,1), and α+αc+αv=1. Individuals in the Sr, Sc and Sv classes become infected either by contact with infectious individuals (direct transmission) at rates ρr,ρc and ρv, respectively, or through contact with contaminated water (indirect transmission) at rates βr,βc and βv, respectively. We assume that individuals in the Sr class have higher direct transmission rates than those in the Sc and Sv classes, and that susceptible individuals in the Sv class are vaccinated at rate v. The vaccinated individuals progress to the recovered class. Following contact with infectious individuals or contaminated water, the susceptible individuals may become exposed to the infection. The time period from cholera pathogen exposure to the development of symptoms (incubation period) is relatively short, varying from about 12 hours to five days [1]. During the pathogen incubation period, exposed individuals are not infectious and do not shed pathogen into the water compartment. However, surviving exposed individuals may progress to the infectious class after the pathogen incubation period at rate σ. Infectious individuals can contaminate the water by shedding the cholera pathogen at rate ξ. An infectious individual can thus generate secondary infections in two ways: through direct contact with susceptible individuals, and by first shedding the pathogen into the water compartment, where it can eventually infect susceptible individuals. We assume that in Cameroon, the number of cholera infections from cholera infected corpses is very small and can thus be ignored. Cholera is a treatable disease. In Cameroon cholera is usually treated with oral rehydration salt solutions and antibiotics. In our model, infectious individuals can progress to the recovered class at rate γ. However, recovered individuals do not have permanent immunity, and can progress back to the Sr, Sc and Sv classes at rate τ(1(ηc+ηv)), τηc and τηv. The parameter δ represents natural death rate in each component and μ is the disease-related death rate. For the pathogen compartment, the concentration of cholera pathogen grows at rate g and decays at rate d. We let λ=dg. From the schematic diagram in Figure 1, with

    φ(t)=ρrI(t)+βrB(t),φc(t)=ρcI(t)+βcB(t),φv(t)=ρvI(t)+βvB(t),

    we obtain the following cholera model with cultural practices and vaccination:

    dSrdt=αβ(ρrI+βrB)Sr(δ+msc)Sr+mcsSc+τ(1(ηc+ηv))RdScdt=αcβ(ρcI+βcB)Sc(δ+mcv+mcs)Sc+mvcSv+mscSr+τηcRdSvdt=αvβ(ρvI+βvB)Sv(δ+mvc)Sv+mcvSc+τηvRvSvdEdt=(ρrI+βrB)Sr+(ρcI+βcB)Sc+(ρvI+βvB)Sv(σ+δ)EdIdt=σE(γ+μ+δ)IdRdt=γI(δ+τ)R+vSvdBdt=ξIλB, (2.1)

    with initial conditions:

    Sr(0)>0,Sc(0)>0,Sv(0)>0,E(0)0,I(0)0,R(0)0,B(0)0.

    From the B equation, when there is no shedding and λ<0, we see that the pathogen maintains itself in the environment. In our application, gd=0.33 [18]. Consequently, we assume throughout that λ=dg>0.

    To obtain the feasible region of solutions to model (2.1), we add the equations for Sr, Sc, Sv, E, I and R to obtain

    dNdt=βδNμI

    and noting that I(t)0 for all t, we get

    dNdtβδN.

    Thus

    lim suptN(t)βδ=N.

    Consequently, the B equation implies

    dBdtξNλB,

    and

    lim suptB(t)ξNλ=B.

    Assuming 0B(0)B, we have B(t)B for all t0. Then the feasible region of solutions of Model (2.1) is the compact set

    Ω={(Sr,Sc,Sv,E,I,R,B)R7+|0Sr,Sc,Sv,E,I,RN,0BB}.

    From the structure of our system of differential equations, with non-negative initial conditions the solutions can be shown to remain non-negative, and we have that the upper bounds in the feasible region are valid [19].

    When there are no cultural practices and no vaccination, and there is only one susceptible class S, Model (2.1) reduces to the following system of ordinary differential equations:

    dSdt=β(ρsI+βsB)S+τRδSdEdt=(ρsI+βsB)S(σ+δ)EdIdt=σE(γ+μ+δ)IdRdt=γI(δ+τ)RdBdt=ξIλB. (2.2)

    with initial conditions:

    S(0)=S0>0,E(0)=E00,I(0)=I00,R(0)=R00,B(0)=B00.

    In Model (2.1), ρs and βs are the direct and indirect transmission rates, respectively. The disease-free equilibrium (DFE) for model (2.2) is

    P0=(S0,E0,I0,R0,B0)=(β/δ,0,0,0,0).

    We use the next generation matrix (NGM) method to compute the basic reproduction number, R0 [13,20,21,22]. Assuming that λi in the B differential equation is not a new infection, then by the NGM method,

    R0=RI0+RB0, (2.3)

    where

    RI0=σρsS0(σ+δ)(γ+μ+δ)andRB0=ξσβsS0λ(σ+δ)(γ+μ+δ).

    The basic reproduction number, R0, is the sum of two terms, RI0 and RB0, with RI0 accounting for the infections from the infectious class, and RB0 accounting for the infections from the contaminated environment. Consequently, if R0<1 then the DFE, P0, is locally asymptotically stable, there is no cholera invasion, and the number of cholera infections eventually decrease. If R0>1, then P0 is unstable, cholera invades the Far North Region of Cameroon and the number of cholera infections increases. We note that if RI0>1 or RB0>1, then R0>1. Furthermore, it is possible for RI0<1, RB0<1 but R0>1.

    To compute R0 for the Far North Region of Cameroon, in Section 2.2, we use the daily reported cholera cases from the three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region from June 25 to August 16, 2019, from the Cameroon Ministry of Public Health [17] to estimate the remaining Model (2.2) parameter values. We will use model (2.2) to illustrate that in the Far North Region of Cameroon, R0>1 and cholera is an endemic infectious disease.

    Some model parameters are known from the literature and some are estimated from our data. The mean exposed period of cholera, 1/σ, is between 12 hours to 5 days [1]. In our simulation 1/σ=1.4 days, the estimated value for the median exposed period of toxigenic cholera [23]. For the immune period, 1/τ=103 days [24], for the net decay rate of pathogen, λ=0.33 day1, and for the shedding rate ξ=10 cells ml1 day1 individuals1 [18]. For the mean infectious period, 1/γ=5 days [14,25]. For the cholera-related death rate, μ=5.5329×104 per day [15]. We estimate the remaining model parameters using the multistart algorithm with the fmincon function in the MATLAB optimization toolbox. The fmincon function takes an iterative approach to solving an optimization problem that includes our parameter values as unknowns with the scalar value to be minimized being the relative error,

    Relative error=||YI(m)YI||2||YI||2, (2.4)

    where YI(m) is the vector of the cumulative number of daily infections given by the Model (2.2) with parameter vector m=(ρs,βs), and YI is the corresponding vector of the values of the cumulative daily reported cholera cases [17]. At time t0, the cumulative cases are computed by integrating the term σE in the I-equation of Model (2.2), from 0 to t, using the MATLAB solver ODE45. The total population of the three health districts in Far North Region as of August 17, 2019 is 89884 [17]. The time period of our simulation is 53 days (June 25, 2019 to August 16, 2019, ), so N(53)=89884. Here, t0 is in days and t=0 corresponds to June 25, 2019 and t=53 corresponds to August 16, 2019.

    We let

    β/δ=N=N(53)=89884,

    so that β=δ×89884. From Table 2, δ=9.1×103365 day1. Consequently, β=817944.4×103365 day1. We choose N(0) so that N(0)<N(53). So we let N(0)=89548. For the other initial conditions, we choose E(0)=0, I(0)=1, R(0)=0, S(0)=N(0)(E(0)+I(0)+R(0)). Similar initial conditions give similar results. As in [15], we take the total initial cholera bacteria concentration B(0)=1000. Table 2 gives estimates for the remaining Model (2.2) parameter values, ρs=2.1737×106 and βs=1.8694×108. Using these values and the known parameter values on Table 2, we obtain that R0=1.2279>1. That is, Model (2.1) predicts that cholera is endemic in the three health district Far North Region of Cameroon.

    Table 2.  Model without culture (Model (2.2)) parameters, descriptions and values.
    Parameter Description Value Source
    β constant recruitment rate 817944.4×103365 day1 estimated
    δ natural death rate 9.1×103365 day1 [26]
    ρs direct transmission rate between S and I 2.1737×106day1 estimated
    βs indirect transmission rate between S and I 1.8694×108 day1 estimated
    1/τ immune period 103 days [24]
    1/σ mean exposed period 1.4 days [23]
    1/γ mean infectious period 5 days [14,25]
    μ cholera-related death rate 5.5329×104 days [15]
    ξ shedding rate of pathogens by I 10 cells ml1day1individual1 [18]
    λ net decay rate of pathogen 0.33 day1 [18]

     | Show Table
    DownLoad: CSV

    In Figure 2, we illustrate the daily cumulative number of new cholera infections predicted by Model (2.2) in comparison to the actual cumulative daily cholera reported cases in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon from June 25 to August 16, 2019. From Figure 2, we see that with relative error 0.1623, Model (2.2) appears to predict an increasing trend in the cumulative number of daily cholera infections in the three health districts of the Far North Region of Cameroon, as in the actual reported cumulative cholera cases per day.

    Figure 2.  Scatterplot of the cumulative values of cholera disease cases in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon per day for June 25 to August 16, 2019 from [17] plotted against simulation output using estimated parameters.

    In order to compute the disease-free equilibrium (DFE) of the structured model with cultural practices, we introduce the following notation:

    A=δ+msc,C=τ(1(ηc+ηv)),D=δ+mcv+mcs,F=δ+mvc+v,J=δ+τ,K=αβ,L=αcβ,M=αvβ,P=JFτηvv,SN=JmscmcsMJmcvmscKAJDMAJmcvLSD=AJmcvmvc+msc(mcsP+mcvCv)+AmcvτηcvADP

    Thus, the DFE of model (2.1) is

    Pc=(Scr,Scc,Scv,Ec,Ic,Rc,Bc)=(Scr,Scc,Scv,0,0,Rc,0),

    where

    Scv=SNSD,Scr=KAmcsMAmcv+(mcsP+mcvCvAJmcv)Scv,Scc=(PJmcv)ScvMmcv,Rc=vJScv

    As in the model without culture, using the next generation matrix (NGM) method, we compute Rc, the control reproduction number. Assuming that λi is not a new infection,

    Rc=RIc+RBc, (2.5)

    where

    RIc=σ(ρrScr+ρcScc+ρvScv)(σ+δ)(γ+μ+δ)andRBc=σξ(βrScr+βcScc+βvScv)λ(σ+δ)(γ+μ+δ).

    Again, as in the Model without cultural practices, the control reproduction number, Rc, is the sum of two terms, RIc and RBc. RIc accounts for the infections from the infectious class, and RBc accounts for the infections from the contaminated environment.

    In Section 2.4, we compute Rc for the Far North Region of Cameroon using the daily reported cholera cases from the three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region from June 25, 2019 to August 16, 2019 [17], and estimate the remaining parameter values of model (2.1). Again, we will use model (2.1) to illustrate that in the Far North Region of Cameroon, Rc<1 and with vaccination, cholera can be eradicated. Moreover, we will use AIC to illustrate that model (2.1) better captures the impact of cultural practices on the dynamics of cholera infections in the three health districts (Kaélé, Kar Hay and Moutourwa) in the Far North Region, which is in agreement with the results of the field work in [8].

    As in model (2.2), we estimate the remaining model parameters using the multistart algorithm with the fmincon built-in function in the MATLAB optimization toolbox. In equation (2.4), YI(m) is the vector of the cumulative number of daily infections given by model (2.1) with parameter vector m, and YI is the corresponding vector of the values of the cumulative daily reported cholera cases [17]. Again, as in Model (2.2), at time t0, the cumulative cases are computed by integrating the term σE in the I-equation of Model (2.1), from 0 to t, using the MATLAB solver ODE45. Also, as in Model (2.2) N(0)=89548. For the other initial conditions, we choose the initial susceptible population so that Sr(0)>Sc(0)>Sv(0). Consequently, we start with Sc(0)=4500, Sv(0)=2100 E(0)=0, I(0)=1, R(0)=0, Sr(0)=N(0)(Sc(0)+Sv(0)+I(0)+E(0)+R(0)). Similar initial conditions give similar results. As in [15], we take the total initial cholera bacteria concentration B(0)=1000.

    In our simulations, to capture the fact that susceptible individuals in the Sr class (who because of their cultural beliefs are not careful in practicing safe cholera intervention measures and not willing to be vaccinated) have higher direct transmission rates than those in the Sc and Sv classes, we let ρc<ρr and ρv<ρr.

    Table 3 gives estimates for the remaining parameter values of model (2.1). Using these parameter values and equation (2.5), we estimate that Rc=0.0873<1. That is, model (2.1) predicts that with vaccination, cholera can be eradicated in the three health district Far North Region of Cameroon. In Figure 3, we illustrate the daily cumulative number of new cholera infections predicted by Model (2.1) in comparison to the actual cumulative daily cholera reported cases in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon for June 25 to August 16, 2019.

    Table 3.  Model with culture (Model (2.1)) parameters, descriptions and values.
    Parameter Description Value Source
    β constant recruitment rate 817944.4×103365 day1 day1 estimated
    δ natural death rate 9.1×103365 day1 [26]
    ρr direct transmission rate between Sr and I 3.0043×108day1 estimated
    ρc direct transmission rate between Sc and I 2.9990×108 day1 estimated
    ρv direct transmission rate between Sv and I 1.3225×109 day1 estimated
    βr indirect transmission rate between Sr and I 1.0023×1011 day1 estimated
    βc indirect transmission rate between Sc and I 1.0002×1010 day1 estimated
    βv indirect transmission rate between Sv and I 1.0×106 day1 estimated
    1/τ immune period 103 days [24]
    1/σ mean exposed period 1.4 days [23]
    1/γ mean infectious period 5 days [14,25]
    μ cholera-related death rate 5.5329×104 days [15]
    ξ shedding rate of pathogens by I 10 cells ml1day1individual1 [18]
    λ net decay rate of pathogen 0.33 day1 [18]
    msc transition rate from Sr to Sc 0.1538day1 estimated
    mcs transition rate from Sc to Sr 0.0001day1 estimated
    mcv transition rate from Sc to Sv 0.03 day1 estimated
    mvc transition rate from Sv to Sc 0.0001 day1 estimated
    ηc fraction of waning immunity from R to Sc 0.0523 estimated
    ηv fraction of waning immunity from R to Sv 0.70004 estimated
    α fraction of recruitment rate into Sr 0.9546 estimated
    αc fraction of recruitment rate into Sc 0.0321 estimated
    αv fraction of recruitment rate into Sv 0.0133 estimated
    v vaccination rate of Sv 0.1568 estimated

     | Show Table
    DownLoad: CSV
    Figure 3.  Scatterplot of the cumulative values of cholera disease cases in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon per day for June 25 to August 16, 2019 from [17] plotted against simulation output using estimated parameters. With the relative error 0.0676 and Rc=0.0873<1, model (2.1) is predicting an increasing trend in the daily number of cumulative cholera infections, like in the observed cumulative number of cholera infections.

    To compare our model without culture to the model with culture, in addition to their relative errors, we compute correct Akaike information criterion (AIC),

    AIC=nlog(ni=1(YI(i)YI(i))2n)+2k,

    where n is the number of days of the cholera infection, and for i=1,,n, YI(i) is of the actual cumulative number of daily reported cholera cases [17] and YI(i) is the cumulative number of daily infections given by our model, and k is the number of estimated parameters [27]. Table 4 summarizes the relative and AIC values for the model without culture and the model with culture.

    Table 4.  Comparative model fit for the two models.
    Model without cultureModel with culture
    n5353
    k216
    Relative error0.16230.0676
    AIC109.924797.6018

     | Show Table
    DownLoad: CSV

    From Table 4, we see that in comparison to the Model without culture, our Model with culture gives us a smaller relative error and AIC value.

    From Figures 2 and 3, we see that as in the actual reported cumulative cholera cases, both models without culture and with culture appear to predict an increasing trend in the cumulative number of daily cholera infections in the three health districts of the Far North Region of Cameroon. However, the model with culture (model (2.1)) gives us a better fit than the model without culture (model (2.2)), with a smaller relative error 0.0676 and smaller AIC = 97.6018. Thus, our model (2.1) captures the influence of cultural practices on the spread of cholera infections in the Far North Region as reported in [8].

    In this section, we use sensitivity analysis to study the impact of each model parameter on the threshold parameter, Rc, and the number of model predicted cholera cases. We also explore the effect of cultural practices on the spread of cholera by generating contour plots of the basic reproduction number (Rc) with respect to two of the parameters in Table 3.

    For a variable v and a parameter p, the relative sensitivity index to the parameter p is defined as

    Υvp=vp×pv.

    This index gives the proportional rate of change of v as p changes [28]. We compute the relative sensitivity indices for Rc with respect to each of the parameters of Table 3. Most of the relative sensitivity indices are complicated algebraic expressions with no obvious structure. As in [28], we evaluate the relative sensitivity indices at the baseline parameter values given in Table 3. Table 5 summarizes our results, with the parameters ordered from the most sensitive β to the least sensitive βr. From Table 5, ΥRcβ=+1.000. Thus, decreasing (or increasing) the most sensitive parameter β, the constant recruitment rate, by 10% decreases (or increases) Rc by 10.0000%. The least sensitive parameter is βr, the indirect transmission rate between the Sr and I classes. Because

    ΥRcβr=+2.7131×106,
    Table 5.  Sensitivity indices of Rc with respect to the parameters of Model (2.1), evaluated at the baseline parameter values given in Table 3.
    ParameterSensitivity Index
    β+1.0000
    ξ+9.9805×101
    λ9.9805×101
    βv+9.9788×101
    γ9.9712×101
    v9.9149×101
    δ9.7684×101
    τ+9.5790×101
    α+9.5459×101
    αc+3.2084×102
    ηv+2.2814×102
    αv+1.3321×102
    mcv+8.8344×103
    μ2.7585×103
    ρc+1.6370×103
    msc+1.4683×103
    ηc+2.8175×104
    ρr+2.6835×104
    βc+1.6544×104
    ρv+4.3551×105
    σ+3.4903×105
    mvc1.8273×105
    mcs5.9587×106
    βr+2.7131×106

     | Show Table
    DownLoad: CSV

    decreasing (or increasing) βr by 10% decreases (or increases) Rc by +2.7131×105%.

    We perform a global sensitivity analysis for our model using Latin Hypercube Sampling (LHS) to sample the parameter space and Partial Rank Correlation Coefficients (PRCC) to evaluate the sensitivity of the outcome variable, the total number of model predicted cholera cases, to uncertainty in the input variables [29,30,31]. For parameter intervals, we chose to go 50% above and below the values in Table 3. Uniform probability distributions were used for each parameter interval. From the recommendation in [32], we take N>4M/3 draws of the LHS design, where M is the number of input parameters and N is the number of LHS draws. In our case, M=24 and N=80. Table 6 gives the PRCCs and the p-values for each model parameter while Figure 4 show p-values of the sensitivity of the number of cholera cases to changes in parameters in Table 6 as computed by the Latin Hypercube Sampling-Partial Rank Correlation Coefficient (LHS-PRCC) index. Based on the p-values, we observe that PRCC values for βv, 0.5956; σ, 0.8981; ξ, 0.5172; λ, 0.9749 and v, 0.5769 are statistically significant. The parameters βv, σ, ξ and v show a significant positive correlation with the total number of cholera cases while λ shows a significant negative correlation with the total number of cholera cases. This is in agreement with our model predictions.

    Table 6.  Model parameters, corresponding PRCC and corresponding p-values resulting from the sensitivity analysis. Significant p-values (p<0.01) are bold-faced.
    ParameterPRCCP-value
    β0.17090.2037
    δ0.074180.5837
    ρr0.04690.7289
    ρc0.06890.6101
    ρv0.03710.7839
    βr0.01170.9311
    βc0.01840.8920
    βv0.59560.0000
    τ0.05940.6607
    σ0.89810.0000
    γ0.19730.1413
    μ0.15740.2422
    ξ0.51720.0000
    λ0.97490.0000
    msc0.19830.1392
    mcs0.03660.7867
    mcv0.11860.3794
    mvc0.18640.1649
    ηc0.25640.8499
    ηv0.02690.8426
    α0.13510.3163
    αc0.01590.9061
    αv0.00950.9438
    v0.57690.0000

     | Show Table
    DownLoad: CSV
    Figure 4.  Sensitivity of the number of cholera cases to changes in parameters in Table 3 as computed by the Latin Hypercube Sampling-Partial Rank Correlation Coefficient (LHS-PRCC) index.

    In this section, we explore the effects of cultural practices on the spread of cholera by generating a contour plot of the control reproduction number (Rc) with respect to the two parameters, mcs and mcv, in Table 3.

    A contour plot of the reproduction number, Rc, as a function of the transition rates mcs (from Sc to Sr) and mcv (from Sc to Sv) is depicted in Figure 5. Using the parameter values in Table 3, our contour plot indicates that an increase in the transition rate of susceptible individuals from being careful in practicing safe cholera intervention measures but unwilling to be vaccinated (Sc) to those who are not careful in practicing safe cholera intervention measures and not willing to be vaccinated (Sr), and a decrease in the transition rate of susceptible individuals from being careful in practicing safe cholera intervention measures but unwilling to be vaccinated (Sc) to those who are careful in practicing safe cholera intervention measures and are willing to be vaccinated (Sv) results in a control reproduction number of Rc=0.09<1. Even though there seem to be an increase in the control reproduction number, Rc is smaller than 0.1, which corroborates with our data.

    Figure 5.  Contour plot of Rc as a red function of mcs and mcv.

    We introduced a model for the dynamics of cholera in the Far North Region of Cameroon that captures cultural practices and vaccination of cholera in the Far North Region. Our model incorporates direct (human-to-human) and indirect (contaminated water-to-human) infection pathways. We use our model to capture the reported cholera cases in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon from June 25 to August 16, 2019. Using our model parameters we obtain that the control reproduction number Rc<1, and our model with cultural practices and vaccination predicts that cholera can be eradicated. The AIC shows that our model with cultural practices and vaccination appears to be a better fit to cumulative number of daily cholera infections in the three health districts of the Far North Region of Cameroon. Consequently, we use model (2.1) to study the effect of cultural practices on cholera cases [8].

    We perform sensitivity analysis to determine the impact of each model parameter on the threshold parameter, Rc, and on the total number of model predicted cholera cases in in the three health districts (Kaélé, Kar Hay and Moutourwa) of the Far North Region of Cameroon for June 25 to August 16, 2019. For Rc, the most sensitive parameter is β (constant recruitment rate), and the least sensitive parameter βr, the indirect transmission rate between the Sr and I classes. The statistically significant parameters to the number of model predicted cholera cases are βv, σ, ξ, λ and v. The parameters βv, σ, ξ and v show a significant positive correlation with the total number of cholera cases while λ shows a significant negative correlation with the total number of cholera cases. This is in agreement with our model predictions.

    A contour plot for the control reproduction number as a function of transition rates between different compartments of susceptible individuals suggests that as humans transition from a class of susceptible individuals who are careful in practicing safe cholera intervention measures but unwilling to be vaccinated (Sc) to those who are not careful in practicing safe cholera intervention measures and not willing to be vaccinated (Sr), the control reproduction number is below one. Thus, our work suggests that in the Far North Region of Cameroon, while implementing intervention schemes such as education and vaccination, public health officials should take into account cultural practices in the region. Since some of the cultural practices stem from religious beliefs of the populace, public health officials should consider working with local chiefs and religious priests when designing the intervention schemes for the control of cholera in the region.

    Our models do not take into account the impact of the cost of intervention strategies on the number of cholera infections. In addition to demonstrating that cultural practices in the Far North region influence cholera transmission, Ngwa et al. in [8] also noted that many people in the Far North had limited scientific knowledge about the cholera disease and its transmission. As more data becomes available, a study of a cholera model that captures education and cultural practices, and cost analysis of intervention strategies in the Far North Region of Cameroon will be an interesting extension of this work.

    The authors declare there is no conflict of interests.



    [1] WHO, Cholera, available from: https://www.who.int/immunization/diseases/cholera/en/.
    [2] WHO, Weekly Epidemiological Record, No 36, (2017), 92,521–536, available from: http://www.who.int/wer.
    [3] D. Posny, C. Modnak, J. Wang, A multigroup model for cholera dynamics and control, Int. J. Biomath., 9 (2016), 1650001. doi: 10.1142/S1793524516500017
    [4] M. Ali, A. R. Nelson, A. L. Lopez, D. A. Sack, Updated Global Burden of Cholera in Endemic Countries, PLoS Negl. Trop. Dis., 9 (2015), e0003832. doi: 10.1371/journal.pntd.0003832
    [5] UNICEF, Cholera epidemiology and response factsheet Cameroon, available from https://www.unicef.org/cholera/files/UNICEF-Factsheet-Cameroon-EN-FINAL.pdf.
    [6] WHO, Cholera, available from: https://www.who.int/news-room/fact-sheets/detail/cholera.
    [7] A. Manguvo, B. Mafuvadze, The impact of traditional and religious practices on the spread of Ebola in West Africa: time for a strategic shift, Pan Afr. Med. J., 22 (2015), 9.
    [8] M. C. Ngwa, A. Young, S. Liang, J. Blackburn, A. Mouhaman, J. G. Morris, Cultural influences behind cholera transmission in the Far North Region, Republic of Cameroon: a field experience and implications for operational level planning of interventions, Pan Afr. Med. J., 08 (2017), 311.
    [9] M. K. Nfor, Recurrent Cholera Outbreak in Far North Cameroon Highlights Development Gaps, Reliefweb, 2014, available from: https://reliefweb.int/report/cameroon/recurrent-cholera-outbreak-far-north-cameroon-highlights-development-gaps
    [10] K. R. Fister, H. Gaff, S. Lenhart, E. Numfor, E. Schaefer, J. Wang, Optimal Control of Vaccination in an Age-structured Cholera Model, In: Chowell G., Hyman J. (eds), Mathematical Modeling of Emerging and Re-emerging Infectious Diseases, Springer, Cham, (2016), 221–248.
    [11] M. R. Kelly, J. H. Tien, M. C. Eisenberg, S. Lenhart, The impact of spatial arrangements on the epidemic disease dynamics and intervention strategies, J. Biol. Dyn., 10 (2016), 222–249. doi: 10.1080/17513758.2016.1156172
    [12] J. H. Tien, D. J. D. Earn, Multiple transmission pathways and disease dynamics in a waterborne pathogen model, Bull. Math. Biol., 72 (2010), 1502–1533.
    [13] P. van den Driessche. Reproduction numbers of infectious disease models, Infect. Dis. Model., 2 (2017), 288–303.
    [14] J. Wang, C. Modnak, Modeling Cholera Dynamics with Controls, Canadian Appl. Math. Q., 19 (2011), 255–273.
    [15] E. Che, Y. Kang, A. Yakubu, Risk structured model of cholera infections in Cameroon, Math. Biosci., 320 (2020), 108303. doi: 10.1016/j.mbs.2019.108303
    [16] M. Al-arydah, A. Mwasa, J. M. Tchuenche, R. J. Smith, Modeling Cholera Disease with Education and Chlorination, J. Biol. Syst., 21 (2013), 1340007. doi: 10.1142/S021833901340007X
    [17] Ministry of Public Health Cameroon, Rapport de situation No10 Épidémie de choléra dans la région de l'Extrême-Nord, 17 août 2019. Available from https://reliefweb.int/report/cameroon/rapport-de-situation-n-10-pid-mie-de-chol-ra-dans-la-r-gion-de-l-extr-me-nord-17-ao.
    [18] C. T. Codeço, Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir, BMC Infect. Dis., 1 (2001), no. 1. doi: 10.1186/1471-2334-1-1
    [19] N. Siewe, A. A. Yakubu, A. R. Satoskar, A. Friedman. Immune response to infection by leishmania: A mathematical model, Math. Biosci., 276 (2016), 28–43. doi: 10.1016/j.mbs.2016.02.015
    [20] O. Diekmann, J. P. Heesterbeek, Mathematical epidemiology of infectious diseases, Wiley, (2000).
    [21] O. Diekmann, H. Heesterbeek, T. Britton, Mathematical tools for understanding infectious disease dynamics, Princeton University Press, (2012).
    [22] P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48. doi: 10.1016/S0025-5564(02)00108-6
    [23] A. S. Azman, K. E. Rudolph, D. A. T. Cummings, J. Lessler, The incubation period of cholera: A systematic review, J. Infect., 66 (2013), 432–438. doi: 10.1016/j.jinf.2012.11.013
    [24] M. Bani-Yaghoub, R. Gautam, Z. Shuai, P. van den Driessche, R. Ivanek, Reproduction numbers for infections with free-living pathogens growing in the environment, J. Biol. Dyn., 6 (2012), 923–940. doi: 10.1080/17513758.2012.693206
    [25] D. M. Hartley, J. G. Morris Jr., D. L. Smith, Hyperinfectivity: A critical element in the ability of V. cholerae to cause epidemics?, PLoS Med., 3 (2005), e7. doi: 10.1371/journal.pmed.0030007
    [26] Knoema, World Data Atlas, Cameroon. Available from: https://knoema.com/atlas/Cameroon.
    [27] K. L. Webster, I. F. Creed, M. D. Skowronski, Y. H. Kaheil, Comparison of the Performance of Statistical Models that Predict Soil Respiration from Forests, Soil Sci. Soc. Am. J., 73 (2009), 1157–1167. doi: 10.2136/sssaj2008.0310
    [28] N. Chitnis, J. M. Hyman, J. M. Cushing, Determining Important Parameters in the Spread of Malaria Through the Sensitivity Analysis of a Mathematical Model, Bull. Math. Biol., 70 (2008), 1272–1296. doi: 10.1007/s11538-008-9299-0
    [29] S. M. Blower, H. Dowlatabadi, Sensitivity and uncertainty analysis of complex models of disease transmission: an HIV model, as an example, Int. Stat. Rev., 62 (1994), 229–243. doi: 10.2307/1403510
    [30] C. Edholm, B. Levy, A. Abebe, T. Marijani, S. Le Fevre, S. Lenhart, et al., A Risk-Structured Mathematical Model of Buruli Ulcer Disease in Ghana. In Kaper H., Roberts F. (eds) Math. Planet Earth, Springer, Cham, 5 (2019), 109–128.
    [31] S. Marino, I. B. Hogue, C. J. Ray, D. E. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254 (2008), 178–196. doi: 10.1016/j.jtbi.2008.04.011
    [32] M. D. McKay, R. J. Beckman, W. J. Conover, A comparison of three methods for selecting values of input variables in the analysis of output from a computer code, Technometrics, 21 (1979), 239–245.
    [33] L. Allen, P. van den Driessche, The basic reproduction number in discrete-time epidemic models, J. Differ. Equ. Appl., 14 (2008), 1127–1147. doi: 10.1080/10236190802332308
  • This article has been cited by:

    1. Folashade B. Agusto, Eric Numfor, Karthik Srinivasan, Enahoro A. Iboi, Alexander Fulk, Jarron M. Saint Onge, A. Townsend Peterson, Impact of public sentiments on the transmission of COVID-19 across a geographical gradient, 2023, 11, 2167-8359, e14736, 10.7717/peerj.14736
    2. Jin Wang, Mathematical Models for Cholera Dynamics—A Review, 2022, 10, 2076-2607, 2358, 10.3390/microorganisms10122358
    3. Leul Mekonnen Anteneh, Bruno Enagnon Lokonon, Romain Glèlè Kakaï, Modelling techniques in cholera epidemiology: A systematic and critical review, 2024, 373, 00255564, 109210, 10.1016/j.mbs.2024.109210
    4. Makayla Preston, Alexandria Carter, Eric Numfor, Modeling the Effects of Media Awareness on SARS-CoV-2 Transmission in Georgia, 2024, 10, 2349-5103, 10.1007/s40819-024-01759-9
    5. Emily Howerton, Kyle Dahlin, Christina J. Edholm, Lindsey Fox, Margaret Reynolds, Brandon Hollingsworth, George Lytle, Melody Walker, Julie Blackwood, Suzanne Lenhart, The effect of governance structures on optimal control of two-patch epidemic models, 2023, 87, 0303-6812, 10.1007/s00285-023-02001-8
    6. Phoebe Amadi, George Lawı, Job Bonyo, A Metapopulation Model for Cholera with Variable Media Efficacy and Imperfect Vaccine, 2024, 7, 2636-8692, 20, 10.33187/jmsm.1289684
    7. Mark Z. Wang, Christina J. Edholm, Lihong Zhao, Using mathematical modeling to study the dynamics of Legionnaires' disease and consider management options, 2025, 22, 1551-0018, 1226, 10.3934/mbe.2025045
    8. Shuo Tian, Yaxin Zhou, Daqing Jiang, Dynamical behaviors of a stochastic SIRB model for cholera with general infection dynamics and Ornstein-Uhlenbeck process, 2025, 0924-090X, 10.1007/s11071-025-11227-8
    9. Murshed Ahmed Ovi, Andrei Afilipoaei, Hao Wang, Estimating Hidden Cholera Burden and Intervention Effectiveness, 2025, 87, 0092-8240, 10.1007/s11538-025-01460-y
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3611) PDF downloads(162) Cited by(9)

Figures and Tables

Figures(5)  /  Tables(6)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog