On the transmission dynamics of Buruli ulcer in Ghana: Insights through a mathematical model.

BACKGROUND: Mycobacterium ulcerans is know to cause the Buruli ulcer. The association between the ulcer and environmental exposure has been documented. However, the epidemiology of the ulcer is not well understood. A hypothesised transmission involves humans being bitten by the water bugs that prey on mollusks, snails and young fishes.
METHODS: In this paper, a model for the transmission of Mycobacterium ulcerans to humans in the presence of a preventive strategy is proposed and analysed. The model equilibria are determined and conditions for the existence of the equilibria established. The model analysis is carried out in terms of the reproduction number [Formula: see text]. The disease free equilibrium is found to be locally asymptotically stable for [Formula: see text] The model is fitted to data from Ghana.
RESULTS: The model is found to exhibit a backward bifurcation and the endemic equilibrium point is globally stable when [Formula: see text] Sensitivity analysis showed that the Buruli ulcer epidemic is highly influenced by the shedding and clearance rates of Mycobacterium ulcerans in the environment. The model is found to fit reasonably well to data from Ghana and projections on the future of the Buruli ulcer epidemic are also made.
CONCLUSIONS: The model reasonably fitted data from Ghana. The fitting process showed data that appeared to have reached a steady state and projections showed that the epidemic levels will remain the same for the projected time. The implications of the results to policy and future management of the disease are discussed.

Background

Buruli ulcer is caused by pathogenic bacterium where infection often leads to extensive destruction of skin and soft tissue through the formation of large ulcers usually on the legs or arms [28]. It is a devastating disease caused by Mycobacterium ulcerans. The ulcer is fast becoming a debilitating affliction in many countries [3]. It is named after a region called Buruli, near the Nile River in Uganda, where in 1961 the first large number of cases was reported. In Africa, close to 30,000 cases were reported between 2005 and 2010 [29]. Cote d’Ivoire, with the highest incidence, reported 2533 cases in 2010 [27]. This disease has dramatically emerged in several west African countries, such as Ghana, Cote d’Ivoire, Benin, and Togo in recent years [26].

The transmission mode of the ulcer is not well understood, however residence near an aquatic environment has been identified as a risk factor for the ulcer in Africa [6, 16, 25]. Transmission is thus likely to occur through contact with the environment [20]. Recent studies in West Africa have implicated aquatic bugs as transmission vectors for the ulcer [18, 24]. An attractive hypothesis for a possible mode of transmission to humans was proposed by Portaels et al. [22]: water-filtering hosts (fish, mollusks) concentrate the Mycobacterium ulcerans bacteria present in water or mud and discharge them again to this environment, where they are then ingested by aquatic predators such as beetles and water bugs. These insects, in turn, may transmit the disease to humans by biting [18]. Person to person transmission is less likely. Aquatic bugs are insects found throughout temperate and tropical environments with abundant freshwater. They prey, according to their size, on mollusks, snails, young fishes, and the adults and larvae of other insects that they capture with their raptorial front legs and bite with their rostrum. These insects can inflict painful bites on humans as well. In Ghana, where Buruli ulcer is endemic, the water bugs are present in swamps and rivers, where human activities such as farming, fishing, and bathing take place [18].

Research on Buruli ulcer has focused mainly on the socio-cultural aspects of the disease. The research recommends the need for Information, Education and Communication (IEC) intervention strategies, to encourage early case detection and treatment with the assumption that once people gain knowledge they will take the appropriate action to access treatment early [2]. IEC is defined as an approach which attempts to change or reinforce a set of behaviours to a targeted group regarding a problem. The IEC strategy is preventive in that it has a potential of enhancing control of the ulcer [5]. It is also important to note that Buruli ulcer is treatable with antibiotics. A combination of rifampin and streptomycin administered daily for 8 weeks has the potential to eliminate Mycobacterium ulcerans bacilli and promote healing without relapse.

Mathematical models have been used to model the transmission of many diseases globally. Many advances in the management of diseases have been born from mathematical modeling [11, 12, 14, 15]. Mathematical models can evaluate actual or potential control measures in the absence of experiments, see for instance [19]. To the best of our knowledge very few mathematical models have been formulated to analyse the transmission dynamics of Mycobacterium ulcerans. This could be largely due to the elusive epidemiology of the Buruli ulcer. Aidoo and Osei [3] proposed a mathematical model of the SIR-type in an endeavour to explain the transmission of Mycobacterium ulcerans and its dependence on arsenic. In this paper, we propose a model which takes into account the human population, water bugs as vectors and fish as potential reservoirs of Mycobacterium ulcerans following the transmission dynamics described in [8]. In addition we include the preventive control measures in a bid to capture the IEC strategy. Our main aim is to study the dynamics of the Buruli ulcer in the presence of a preventive control strategy, while emphasizing the role of the vector (water bugs) and fish and their interaction with the environment. The model is then validated using data from Ghana. This is crucial in informing policy and suggesting strategies for the control of the disease.

This paper is arranged as follows; in “Methods”, we formulate and establish the basic properties of the model. We also determine the steady states and analysed their stability. The results of this paper are given in “Results”. Parameter estimation, sensitivity analysis and the numerical results on the behavior of the model are also presented in this section. The paper is concluded in “Discussion”.

Methods

Model formulation

We consider a constant human population the vector population of water bugs and the fish population at any time t. The total human population is divided into three epidemiological subclasses of those that are susceptible the infected and the recovered who are still immune . Total population of vector (water bug) at any time t is divided into two subclasses to susceptible water bugs and those that are infectious and can transmit the Buruli ulcer to humans, The total population reservoir of small fish is also divided into two compartments of susceptible fish and infected fish We also consider the role of the environment by introducing a compartment U,  representing the density of Mycobacterium ulcerans in the environment. We make the following basic assumptions:

Mycobacterium ulcerans are transferred only from vector ( water bug) to the humans.

There is homogeneity of human, water bug and fish populations’ interactions.

Infected humans recover and are temporarily immune, but lose immunity.

Fish are preyed on by the water bugs.

Unlike some bacterial infections such as leprosy (caused by Mycobacterium leprae) and tuberculosis (caused by Mycobacterium tuberculosis), which are characterized by person-to-person contact transmission, it is hypothesized that Mycobacterium ulcerans is acquired through environmental contact and direct person-to-person transmission is rare [20].

Susceptible host (human population) can be infected through biting by an infectious vector (water bug). We represent the effective biting rate that an infectious vector has to susceptible host as and the incidence of new infections transmitted by water bugs is expressed by standard incidence rate One can interpret as a function of the biting frequency of the infected water bugs on humans, density of infectious water bugs per human, the probability that a bite will result in an infection and the efficacy of the IEC strategy. In particular we can set where is the efficacy of the IEC strategy, the number of water bugs per human host, the biting frequency (the biting rate of humans by a single water bug) and the probability that a bite by an infected vector to a susceptible human will produce an infection.

Susceptible water bugs are infected at a rate through predation of infected fish and representing other sources in the environment. Here differentiates the infectivity potential of the fish from that of the environment.

Assuming fish prey on infected water bugs, susceptible fish are infected at a rate through predation of infected fish and representing infection through the environment. Here is a modification parameter that models the relative infectivity of fish from that of the environment.

The vector population and the fish populations are assumed to be constant. The growth functions are respectively given by and where

It is important to note that other types of functions can be chosen as growth functions. In this work we however assume that the growth functions are linear.

There is a proposed hypothesis that environmental mycobacteria in the bottoms of swamps may be mechanically concentrated by small water-filtering organisms such as microphagous fish, snails, mosquito larvae, small crustaceans, and protozoa [8]. We assume that fish increase the environmental concentrations of Mycobacterium ulcerans at a rate Humans are are assumed not to shed any bacteria into the environment.

Aquatic bugs release bacteria into the environment at a rate

The model does not include a potential route of direct contact with the bacterium in water.

The birth rate of the human population is directly proportional to the size of the human population.

The recovery of infected individuals is assumed to occur both spontaneously and through treatment. Research has shown that localized lesions may spontaneously heal but, without treatment, most cases of Buruli ulcer result in physical deformities that often lead to physiological abnormalities and stigmas [4].

We now describe briefly, the transmission dynamics of Buruli ulcer:

New susceptibles enter the population at a rate of Buruli ulcer sufferers do not recover with permanent immunity, they loose immunity at a rate and become susceptible again. Susceptibles and infected through interaction with infected water bugs, with infection driven by water bugs biting susceptible humans. Once infected, individuals are allowed to recover either spontaneously or through antibiotic treatment at a rate In this model, the human population is assumed to be constant over the modeling time with the birth and death rates being equal. The compartment tracks the changes in the susceptible water bugs population that are recruited at a rate . The infection of water bugs is driven by two processes: their interaction infected fish and with the environment. The natural mortality of the water bugs occurs at a rate Similarly, the compartment tracks the changes in the susceptible fish population that are recruited at a rate . The infection of fish is also driven by two processes: their interaction infected water bugs and with the environment. Fish’s natural mortality rate is The growth of Mycobacterium ulcerans in the environment is driven by their shedding by infected water bugs and fish into the environment. They are assumed to die naturally at a rate The possible interrelations between humans, the water bug and fish are represented by the schematic diagram below (Fig. 1).

Fig. 1

Proposed transmission dynamics of the Buruli ulcer among humans, fish, water bugs and the environment (U)

The descriptions of the parameters that describe the flow rates between compartments are given in Table 1.

Table 1

Description of parameters used in the model

Symbol	Description
β_H	The effective contact rate between the vector and susceptible humans
β_V	The effective contact rate between fish and susceptible vectors
β_F	The effective contact rate between the susceptible fish and Mycobacterium ulcerans
γ	The recovery rate of infected humans
θ	The rate of loss of immunity of recovered humans
μ_H	Natural mortality rate/birth rate of the human population
μ_V	Natural mortality rate of the vector population
μ_F	Natural mortality rate of the fish population
r_V	The growth rate of the vector population
r_F	The growth rate of the fish population
K	The environmental carrying capacity of the bacteria population
σ_F	Rate of shedding of Mycobacterium ulcerans into the environment by fish
σ_V	Rate of shedding of Mycobacterium ulcerans into the environment by the water bugs
μ_E	Rate at which Mycobacterium ulcerans are cleared from the environment

The dynamics of the ulcer can be described by the following set of nonlinear differential equations:We assume that all the model parameters are positive and the initial conditions of the model system (1) are given byWe arbitrarily scale the time t by the quantity by letting and introduce the following dimensionless parameters;can be non dimensionalised bySo, system (1) settingThe forces of infection for humans, water bugs and fish are respectivelyGiven that the total number of bites made by the water bugs must equal the number of bites received by the humans, is a constant, see [9]. Similarly is constant and so is We also note that since and are constants, and are constants.

Given that and system (1) can be reduced to the following system of equations by conveniently maintaining the capitalised subscripts so that we can still respectively write and as and

Basic properties

Feasible region

Note that Through integration we obtain The feasible region (the region where the model makes biological sense) for the system (2) is in and is represented by the setwhere the basic properties of local existence, uniqueness and continuity of solutions are valid for the Lipschitzian system (2). The populations described in this model are assumed to be constant over the modelling time. The solutions of system (2) starting in remain in for all Thus , is positively invariant and it is sufficient to consider solutions in

Positivity of solutions

We desire to show that for any non-negative initial conditions of system (2), say the solutions remain non-negative for all We prove that all the state variables remain non-negative and the solutions of the system (2) with positive initial conditions will remain positive for all . We thus state the following lemma.

Lemma 1

Given that the initial conditions of system (2) are positive, the solutionsandare non-negative for all.

Proof

Assume thatThus and it follows directly from the first equation of the system (2) thatWe thus haveSince the exponential function is always positive and the solution will thus be always positive.

From the second equation of (2),Similarly, it can be shown that and for all and this completes the proof.

Steady states analysis

The disease free equilibrium

In this section, we solve for the equilibrium points by setting the right hand side of system (2) to zero. This direct calculation shows that system (2) always has a disease free equilibrium pointWe have the following result on the local stability of the disease free equilibrium.

Theorem 1

The disease free equilibriumwhenever it exists, is locally asymptotically stable ifand unstable otherwise.

The Jacobian matrix of system (2) at the equilibrium point is given byIt can be seen that the eigenvalues of are and the solution of the characteristic polynomialwhereforThe solutions of have negative real parts only if following the use of the Routh Hurwitz Criterion. We can thus conclude that the disease free equilibrium is locally asymptotically stable whenever

We note that is the model system (2)’s reproduction number and does not depend on the human population size. The model reproduction number is a sum of three terms. The terms and represent the contribution of fish and water bugs respectively to the infection dynamics. The term which is not very common in many epidemiological models, shows the combined contribution of the water bugs, fish and their shedding of Mycobacterium ulcerans into the environment. So, the infection is driven by the water bugs, fish and the density of the bacterium in the environment. The model reproduction number increases linearly with the shedding rates of the Mycobacterium ulcerans into the environment by fish and water bugs and the effective contact rates and . It decreases with increasing removal rates of the fish and Mycobacterium ulcerans. So the control of the ulcer depends largely on environmental management.

The endemic equilibrium

The endemic equilibrium is much more tedious to obtain. Given that from the first and second equations of system (2) we havewhere

The last equation of system (2) can be written asWe thus havewhere

From the third and fourth equations of system (2)we haveSubstituting (3) into (4) we obtain and the cubic equationwhereNote thatGiven thatthe turning points of equation (5) are given byThe discriminant of solutions (7) is . We now focus on the sign of the discriminant.

If , then has no real turning points, which implies that is a strictly monotonic function. The sign of is crucial in determining the monotonicity. Through completing the square, equation (6) can be written asClearly if , then . Since , then . Thus is a strictly monotone decreasing function. Note that . For the polynomial has no positive real roots for . However, if it has only one positive real root for and consequently only one endemic equilibrium.

If , then has only one real root with multiplicity two. This implies that and that . Thus the polynomial is a decreasing function. Given that the turning point is a point of inflexion for The polynomial has only one endemic equilibrium.

For , we consider two cases; and . If , then . This means that . Irrespective of the sign of , has two real positive and distinct roots. This implies that (5) has two positive turning points. If i.e then, has at least one positive real root, and hence at least one endemic equilibrium. On the other hand, if then, has at most two positive real roots when , and hence at most two endemic equilibria.

If , then , which implies that . For , has two real roots of opposite signs. Since for , then, has one positive root. For , has two negative real roots. Since for , then, has no positive real roots, and consequently no endemic equilibria.

Furthermore, we can use the Descartes’ Rule of Signs [7] to explore the existence of endemic equilibrium (or equilibria) for . We note the possible existence of backward bifurcation. The theorem below summarises the existence of endemic equilibria of the model system (2).

Theorem 2

The model system (2) has

a unique endemic equilibrium point if.

has two endemic equilibria forwhereis the critical threshold below which no endemic equilibrium exists.

Remark The evaluation of depends on the signs of and and the sign of the discriminant. The computations are algebraically involving and long and are not included here. Since the model system (2) possesses two endemic equilibria when , the model exhibits backward bifurcation for .

The consequence of the above remark is that bringing below unity is not sufficient to eradicate the disease. For eradication, must be brought below the critical value .

Global stability of the endemic equilibrium

Theorem 3

The endemic equilibrium pointof system (2), is globally asymptotically stable.

The global stability of the endemic equilibrium, can be determined by constructing a Lyapunov function such thatThe corresponding time derivative of is given byAt the endemic equilibrium, we have the following relationsEvaluating the components of the time derivative of the Lyapunov function using the relations (11) we haveLetSubstituting (13) into (12), we obtainwhereNext, we choose A, B, C and D so that none of the variable terms of are positive. It is important to group together the terms in that involve the same state variable terms, as well as grouping all of the constant terms together. So we can show that by expanding (15), writing out the constant term and the coefficients of the variable terms such as and so on. The only variable terms that appear with positive coefficients are x, y and z. We thus choose the Lyapunov coefficients so as to make the coefficients ofx, y and z equal to zero. We haveThe coefficients C and D can similarly be evaluated from the coefficients of y and z. Note that expressions such asemanating from the substitution of the coefficients into are less than or equal to zero by the arithmetic mean-geometric mean inequality. This implies that with equality only if

Therefore, and by the LaSalle’s Extension [17], it implies that the omega limit set of each solution lies in an invariant set contained in The only invariant set contained in is the singleton . This shows that each solution which intersects limits to the endemic equilibrium. This completes the proof.

Results

Parameter estimation

The biggest challenge in epidemic modeling is the estimation of parameters in the model validation process. In this section we endeavour to estimate some of the parameter values of system (2). The demographic parameters can be easily estimated from census population data. We begin by estimating the mortality rate We note that the average life expectancy of the human population in Ghana is 60 years [21]. This translates into per year or equivalently per day. Buruli ulcer is currently regarded as a vector borne disease. Recovery rates modelled by of vector borne diseases range from to 0.5 per day [23]. This translates to an average of between 0.00584 and 183 per year. The rate of loss of immunity for vector borne diseases range between 0 and per day[23]. The mortality rate of the water bugs is assumed to be 0.15 per day [3]. The rates per day can easily be transferred to yearly rates.

In this model we shall assume that we have more water bugs than humans so that Since the water bugs prey on the fish, a reasonable food chain structure leads to the assumption that we have more fish than water bugs hence and consequently If the water bug is assumed to interact more with the environment than fish then and The natural mortality of small fish in rivers is not well documented and data on the mortality of river fish in Ghana is not available. For the purpose of our simulations, we shall assume that per day. Given that we have We shall also assume that We summarise the parameters in the following Table 2.

Table 2

Parameter values used for the simulations and sensitivity analysis

Parameter	Value/range	Source
μ_h	4.5×10^-5	[21]
γ_h	1.6×10^-5-0.5	[23]
θ_h	0-1.1×10^-2	[23]
m_1,m_2	m_1<1,m_2>1	Estimated
m_3,m_4,m_5	(0,1)	Estimated
β_h,β_f,β_v	(0,1)	Estimated
η_v	(1,5)	Estimated
η_f	(0, 1)	Estimated
σ_f,σ_v	(0,1)	Estimated
μ_f	3×10^-3-7×10^-3	Estimated
μ_e	(0,1)	Estimated

Sensitivity analysis

Many of the parameters used in this paper are not determined experimentally. Therefore their accuracy is always in doubt. This can be overcome by observing responses of such parameters and their influence on the model variables through sensitivity and uncertainty analysis. In this subsection we present the sensitivity analysis of the model parameters to ascertain the degree to which the parameters affect the outputs of the model. We use the Partial Correlation Coefficients (PRCCs) analysis to determine the sensitivity of our model to each of the parameters used in the model. Through correlations, the association of the parameters and state variables can be established. In our case, we determine the correlation of our parameters and the state variable U. Alongside the PRCCs are the statistical significance test p-values for each of the parameters. If the magnitude of the PRCC value of a parameter is greater than 0.5 or less than −0.5 and the p-value less than 0.05, then the model is sensitive to the parameter. On the other hand, PRCC values close to or indicate that the parameter strongly influences the state variable output. The sign of a PRCC value indicates the qualitative relationship between the parameter and the output variable. A negative sign indicates that the parameter is inversely proportional to the outcome measure [10]. The parameters with negative PRCCs reduce the severity of Burili ulcer disease while those with positive PRCCs aggravate it. Using Latin Hypercube Sampling (LHS) scheme with 1000 simulations for each run, with U as the outcome variable. Our results show that the variable U is sensitive to the changes in the parameters and . The results are shown in Fig. 2.

Fig. 2

PRCC plots: The variable U largely depends on and . The bars pointing to the left indicate that U has an inverse dependence on the respective parameters. We observe that the parameters and aggravate the disease when they are increased while and reduce its severity when increased

The results from the PRCC analysis are summarized in Table 3. The significant parameters together with their PRCC values and p-values have been encircled.

Table 3

Outputs from PRCC analysis

In Fig. 3 the residuals for the ranked Latin Hypercube Sampling parameter values are plotted against the residuals for the ranked density of Mycobacterium ulcerans.The PRCC plots for parameters and show a strong linear correlation. The growth of Mycobacterium ulcerans increases as the number of infected fish that eventually shed bacteria into the environment increases. An increase in the parameters and leads a decrease in amount of bacteria in the environment.

Fig. 3

PRCC plots shows the PRCC plots for the parameters , , , and

Data and the fitting process

One of the most important steps in the model building chronology is model validation. We now focus on the data provided by the Ashanti Regional Disease Control Office for Buruli ulcer cases in Ghana per 10,0000 people. The data are given in the Table 4 below for the years 2003–2012.

Table 4

Data on Buruli ulcer cases in Ghana

Year	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012
Buruli ulcer cases	739	1159	1201	1096	1136	1300	1158	1428	1324	1292

Source [13]

We fit the model system (2) to the data of Buruli ulcer cases expressed as fractions. We use the least squares curve fit routine (lsqcurvefit) in Matlab with optimisation to estimate the parameter values. Many parameters are known to lie within limits. A few parameters such as the demographic parameters are known [13] and it is thus important to estimate the others. The process of estimating the parameters aims at finding the best concordance between computed and observed data. One tedious way to do it is by trial and error or by the use of software programs designed to find parameters that give the best fit. Here, the fitting process involves the use of the least squares-curve fitting method. Matlab code is used where unknown parameter values are given a lower and upper bound from which the set of parameter values that produce the best fit are obtained.

Figure 4 shows how system (2) fits to the available data on the incidence of the BU. The incidence solution curve shows a very reasonable fit to the data.

Fig. 4

Model fit to data. Model system (2) fitted to data of Burili ulcer cases in Ghana. The circles indicate the actual data and the solid line indicates the model fit to the data. The parameter values used for the fitting;

In planning for a long term response to the Buruli ulcer epidemic, it is important to have some reasonable projections to the epidemic. The fitting process allows us to envisage the Buruli ulcer epidemic in future. it is important to note that the projections are reasonably good over a short period of time since the current is evolving gradually based on the available data. We chose to project the epidemic beyond 5 years to 2017. Figure 5 show the projected Buruli ulcer epidemic.

Fig. 5

Projected model fit. Projection to fit in Fig. 4

Figures 6 and 7 show the changes in the prevalence of infected humans respectively when the shedding rate of Mycobacterium ulcerans in the environment and the removal rate of MU from the environment, are varied. Based on the sensitivity analysis, our model is very sensitive to the shedding rate of Mycobacterium ulcerans into the environment. Figure 6 shows that an increase in the shedding rate will lead increased human infections. We can actually quantify the related increases. For instance, if is increased from 0.51 to 0.52 on year 15, the percentage increase in the prevalence of human infections is 6 %. Minimising Mycobacterium ulceransin the environment is an important control measure that is, albeit impractical at the moment. We observe through our results that their decrease in the environment can lead to quantifiable changes in the prevalence of infected humans. Increasing leads to a decrease in the prevalence of infected humans.

Fig. 6

Prevalence of Buruli ulcer infection in humans. Shows prevalence humans when is varied

Fig. 7

Prevalence of Buruli ulcer in infected humans for different values of Shows prevalence the infected humans when is varied

Discussion

In this paper, a deterministic model on the dynamics of the Buruli ulcer in the presence of a preventive intervention strategy is presented. The model’s steady states are determined and their stabilities investigated in terms of the classic threshold In disease transmission modelling, it is well known that a classical necessary condition for disease eradication is that the basic reproductive number must be less than unity. The model has multiple endemic equilibria (in fact it exhibits a backward bifurcation). When a backward bifurcation occurs, endemic equilibria coexist with the disease free equilibrium for This means that getting the classic threshold less than 1, might not be sufficient to eliminate the disease. Thus the existence of backward bifurcation has important public health implications. This might explain why the disease has persisted in the human population over time. The endemic equilibrium is found to be globally stable if

The sensitivity analysis of model parameters showed some interesting results. These results suggest that efforts to remove Mycobacterium ulcerans and infected fish from the environment will greatly reduce the epidemic although the latter will be impracticable. This is because of the costs involved and the fact that many governments in affected areas operate on lean budgets.

The model is then fitted to data on the Buruli ulcer in Ghana. The model reasonably fits the data. The challenge in the fitting process was that the data appears to indicate that Buruli ulcer has reached a steady state. This then produced some parameter values that appeared unreasonable. Despite these challenges, the fit produced reasonable projections on the future of the ulcer. The model shows that in the near future, the number of cases will not change if everything remains the same. An important consideration that can be added to the model is the inclusion of probable policy shifts and the investigation of different scenarios on the progression of the epidemic as the policies change. Because not much of the disease is understood, parameter estimation was a daunting task. So we had to reasonably estimate some of the parameter using the hypothesis that Buruli ulcer is a vector borne disease. Due to the estimation of essential parameters sensitivity analysis was necessary and very important to determine how these parameters influence the model. The implications of varying some of the important epidemiological parameters such as the shedding rates were investigated. Important results were drawn through Figs. 6 and 7. The main result of this paper is that the management of Buruli ulcer depends mostly on the management of the environment.

Conclusions

This model can be improved by considering social interventions in the human population, modeled as functions and the inclusion of the different forms of treatment available as some individuals opt for traditional methods while others depend on the government health care system [1]. Social interventions include education, awareness, poverty reduction and provision of social services. While the mathematical representations of these interventions are insurmountable, they are vital to the dynamics of the disease and public health policy designs. Finally this model can be used to suggest the type of data that should be collected as research on the Buruli ulcer intensifies. The global burden of the disease and its epidemiology are not well understood, [28]. Clearly, gaps do exist in the nature and type of data available. Reports on the disease are often based on passive presentations of patients at health care facilities. As a result of the difficulties of accessing health care in affected areas, data on the disease is scanty.