Scabies is among the infestations almost forgotten due to its association with poor communities. We formulate a deterministic model to assess the possible impact vaccination will have on scabies control. The Descartes’s rule of signs is used to show the nature of the endemic equilibria. Analysis of the reproduction number and numerical simulations suggest that vaccination in addition to treatment will help greatly in reducing the spread of scabies infestation. This suggests there is a strong need for researchers to come up with a possible vaccine in that order to effectively control scabies especially among the disadvantaged communities.
1. Introduction
Scabies is a highly contagious skin infestation caused by skin infestation with the mite Sarcoptes scabiei and is spread by direct contact. Its common name is the “itch mite” derived from the severe pruritis it causes [1]. The female mite burrows the skin and lays its eggs which hatch into larvae within 2-3 days [2]. The larvae then burrow to reach the skin surface, then dig burrows where they pass through two further developmental stages (protonymphs and tritonymphs) before moulting into either male or female mites. As a result a result, S-shaped trails accompanied by small insect bites cover infested areas. It is only less than 1% of the eggs laid by female mites that develop into adult mites [3]. Following infection, the mite population increases rapidly up to 25 adult mites after 50 days and upto 500 mites by 100 days [3], and thereafter rapid decrease may be due to scratching as well as host immune's response. However, 10–12 mites are enough to cause a burden in a host [1].
It is estimated that, despite the availability of effective treatment, over 300 million people in the world are scabies infested [6]. Scabies is well described in the historical record with descriptions evident in ancient texts from China, India and the Middle East [7]. Scabies was also a scourge in London during the First World War upto 1920 [8]. Scabies can be treated with any of the following medications: permethrin (one of the most effective and most expensive), ivermectin (an oral medicine shown in clinical studies to be effective against scabies but not recommended for children), lindane (effective though there are concerns that it can cause neurotoxicity), benzyl benzoate, crotamiton, malathion, and sulphur preparations. Due to low cost, benzyl-benzoate and sulphur ointments are commonly used in developing countries. In general the prevalence of scabies is not well documented especially in Africa [2]. Scabies infestations like other infectious diseases are more common among the poor due to the overcrowded and unhygienic living conditions they are exposed to. However, the issues of associating poverty and scabies are complicated as even in the developed world scabies cases have also being recently reported [9]. In November 2007, the Public Health Service in Amsterdam was notified that the nursing home doctor in Amsterdam had been diagnosed with scabies and that another staff member who worked in the same nursing home had been diagnosed with scabies in September 2007 [9].
Currently, there is no vaccine against scabies; however, the observations that scabies infestation in humans is self limiting and that after the first infestation, significant protective immunity against re-infestation provide theoretical support for the development of a vaccine against scabies [10–12]. A vaccine against scabies mites would help improve the quality of life of poor people in the world who live in overcrowded and unhygenic conditions. Scabies being a neglected disease associated with the poor to the best knowledge of the authors no mathematical insight into its transmission dynamics have been attempted. This is possibly a first attempt to analyse the potential benefits a vaccine will have on scabies control. This work seeks to address whether the current treatment regime against scabies mites is enough or there is a need to come up with new vaccines.
The rest of this paper is organised as follows. In the next section, the model and its basic properties are presented. In Section 3, the equilibria states of the model are determined and detailed analysis of the model is carried out. The possible effects the intervention strategies have on scabies control are assessed through numerical simulations, and finally a discussion in the last section presented.
2. Model Formulation
The population is divided into four epidemiological classes according to an individual's disease status, namely, the unvaccinated susceptibles(Su), vaccinated susceptibles (Sv), infected individuals (I), and the recovered (R). Thus, the total human population is given by N(t)=Su(t)+Sv(t)+I(t)+R(t). The human population is recruited through birth at a rate Λ into the unvaccinated susceptible class. Unvaccinated susceptibles are vaccinated at a rate α to move into the vaccinated susceptible class (Sv). The vaccine wanes at a rate ρ for individuals in Sv to move back into Su class. Following prolonged contact and/or sexual contact with scabies infected person, the unvaccinated susceptibles are infected with scabies at a rate λ to move into the infected class (I) where λ=βI/N with β being the product of the probability of getting infected per contact with an infected case and the effective contact rate necessary for scabies infestation to occur. Even though currently there is no effective vaccine against scabies, we assume vaccinated susceptibles have a reduced rate of becoming infected; thus, vaccinated susceptibles are infected with scabies at a rate δλ with 0<δ<1 signifying the protective factor of the vaccine. The infected individuals are treated and move into the recovered class (R) at a rate γ. Recovery following treatment confers some partial immunity against reinfection [13], so individuals in R(t) class are reinfected at a rate σλ,σ∈(0,1). Individuals in all human subclasses experience natural death at a rate μ. On their own, scabies do not kill people, so there is no disease induced death in this model The model flow diagram is shown in Figure 1.
Structure of the model.
Based on the assumptions, the following system of differential equations describes the model:(1)Su′(t)=Λ+ρSv-λSu-(μ+α)Su,Sv′(t)=αSu-δλSv-(μ+ρ)Sv,I′(t)=λSu+δλSv+σλR-(μ+γ)I,R′(t)=γI-σλR-μR.
Adding all equations in model system (1), we have
(2)N′(t)=Λ-μN⟹N(t)=N(0)e-μt+Λμ(1-e-μt).
Thus, limt→∞N(t)=Λ/μ which is the limiting value of the population. Therefore without loss of generality, we can only consider solutions of model system (2) in the following positively invariant, that is, all solutions in biological space where all the variables are nonnegative:
(3)Ω={(Su,Sv,I,R)∣N≤Λμ}.
3. Disease-Free Equilibrium and Stability Analysis
The disease-free equilibrium of model system (1) is given by
(4)E0=(Su0,Sv0,I0,R0)=(Λ(μ+ρ)μ(μ+ρ+α),Λαμ(μ+ρ+α),0,0).
The reproduction number of model system (1) is given by
(5)ℛPT=β(μ+ρ+δα)(μ+γ)(μ+ρ+α).
Theorem 1 follows from [14, Theorem 2].
Theorem 1.
The disease-free equilibrium ℰ0 is locally asymptotically stable for ℛPT<1 and unstable otherwise.
3.1. Endemic Equilibrium
The endemic equilibrium of model system (1) is given by E*=(Su*,Sv*,I*,R*) with
(6)Su*=Λ(δλ*+μ+ρ)δλ*2+(δ(μ+α)+μ+ρ)λ*+μ(μ+α+ρ),Sv*=Λαδλ*2+(δ(μ+α)+μ+ρ)λ*+μ(μ+α+ρ),I*=(Λλ*(σλ*+μ)(δλ*+μ+ρ+δα))×(μ(σλ*+μ+γ)[δλ*2+(δ(μ+α)+μ+ρ)λ*+μ(μ+α+ρ)δλ*2])-1,I*=(Λλ*γ(δλ*+μ+ρ+δα))×(μ(σλ*+μ+γ)[δλ*2+(δ(μ+α)+μ+ρ)λ*+μ(μ+α+ρ)δλ*2])-1
in terms of the force of infection λ*. Substituting into the equation for the force of infection λ*=βI*/N* to obtain
(7)λ*h(λ*)=λ*(Aλ*3+Bλ*2+Cλ*+D)=0,
where λ*=0 corresponds to the disease-free equilibrium and h(λ*)=0 corresponds to the existence of endemic equilibria where
(8)A=σδμ(μ+γ)(μ+α+ρ),B=δμ(μ+α+ρ)+σ(μ+ρ+δ(μ+α-β))μ(μ+γ)(μ+α+ρ),C=μ+ρ+δαμ(μ+α+ρ)-σℛPTμ+σμ+γ+δ(μ+γ-β)(μ+γ)(μ+α+ρ),D=1-ℛPT.
Solving for λ* in h(λ*)=0, the roots of h(λ*)=0 are explored using the Descartes’s rule of signs. The various possibilities are tabulated in Table 1.
Number of possible positive roots of h(λ*)=0 for ℛPT>1 and ℛPT<1.
Cases
A
B
C
D
ℛPT
Number of sign changes
Number of possible real roots (endemic equilibrium)
1
+
+
+
+
ℛPT<1
0
0
+
+
+
−
ℛPT>1
1
1
2
+
−
+
+
ℛPT<1
2
0,2
+
−
+
−
ℛPT>1
3
1,3
3
+
−
−
+
ℛPT<1
2
0,2
+
−
−
−
ℛPT>1
1
1
4
+
+
−
+
ℛPT<1
2
0,2
+
+
−
−
ℛPT>1
1
1
The analysis of the results in Table 1 results in the following Theorem 2.
Theorem 2.
The scabies model (1)
has a unique endemic equilibrium if ℛPT>1 and whenever Cases 1, 2, 3, and 4 are satisfied;
could have more than one endemic equilibrium if ℛPT>1 and Case 2 is satisfied;
could have two endemic equilibria if ℛPT<1 and cases 2, 3, and 4 are satisfied.
3.2. Analysis of the Reproduction Number ℛPT
In the absence of both vaccination and treatment ℛPT is simply the basic reproduction number ℛ0 which is defined as the number of secondary scabies cases produced by one infectious individual during his/her entire infectious period in the absence of any intervention strategy and is given by
(9)ℛ0=βμ.
Now let us analyse what happens when their vaccination is only in the community. In that scenario, ℛPT becomes ℛP which is defined as the number of secondary scabies cases produced by one infected individual during his/her entire infectious period in presence of vaccination and is given by
(10)limγ→0ℛPT=ℛP=β(μ+ρ+δα)μ(μ+ρ+α)=μ+ρ+δαμ+ρ+αℛ0.
Subtracting ℛP from ℛ0, we have
(11)Δp=ℛP-ℛ0=(1-δ)αμ+ρ+αℛ0>0.
Also, partially differentiating ℛP with respect to α, one gets(12)∂ℛP∂α=-(μ+ρ)(1-δ)(μ+ρ+α)2ℛ0<0.
The fact that (11) is positive and (12) is negative suggests that the increase in the use of the vaccine will increase its benefit to the community. This result suggests more resources should be geared towards the development and deployment of this possible vaccine against scabies. Now we determine the critical vaccination rate necessary for the vaccine to have a positive impact on its control and this is done by expressing α as a subject of the formula in (10) to obtain
(13)αc=(μ+ρ)(1-ℛP/ℛ0)ℛP/ℛ0-δand exists forδ<ℛPℛ0<1.
Thus, as long as α>αc, then vaccination will aid in the control of scabies. However, for α<αc, vaccination levels will not be high enough to make a meaningful impact on scabies control. To effectively determine the effects vaccination has on scabies control, there is a need to know the critical vaccine waning rate which is found by making ρ the subject of formula in the equation of ℛP to obtain
(14)ρc=(μ+δα)(((μ+α)ℛP/(μ+δα)ℛ0)-1)1-ℛP/ℛ0which exists for μ+δαμ+α<ℛPℛ0<1.
Thus, vaccination will be able to control scabies provided ρ<ρc. If ρ>ρc, then the rate at which the vaccine wanes is very high that herd immunity cannot be achieved. Table 2 summarises the effects of vaccination and vaccine waning rates.
Effects of different levels of vaccination and vaccine waning rates.
α<αc
α>αc
ρ<ρc
As much as vaccine waning rates are low, the vaccination rates are not high enough to control scabies
Ideal situation where vaccination rates are high and vaccine waning rates are low. In this case, vaccination leads to the control of scabies.
ρ>ρc
Low vaccination rates with a vaccine which wanes too quickly cannot control scabies in a community
High vaccination rates with a vaccine which wanes too quickly have little effect on scabies control
Results from Table 2 show that an ideal situation for a vaccine to achieve herd immunity is when vaccination coverage rates are high and vaccine waning rates are low (α>αc,ρ<ρc). This result is further illustrated in Figure 2 which shows that high vaccination rates accompanied by low vaccine waning rates are ideal to make the reproduction number ℛP small. A reduction in the size of the reproduction number translates to a reduction in the number of new infections.
Effects of vaccination rate and vaccine waning rate on ℛP.
Since we are talking about vaccination and its potential impact of scabies control, there is a need to find the critical protective factor for a vaccine to be effective. This requires making δ the subject of the formulae in ℛPand doing so one gets
(15)δc=μ+ρα((μ+ρ+α)ℛP(μ+ρ)ℛ0-1)which exists for μ+ρμ+ρ+α<ℛPℛ0<1.
If δ<δc, then the vaccine will be effective enough to control scabies infection. However, when δ>δc, then vaccination will not contribute meaningfully to the reduction of scabies. Looking at the nature of ℛPT shows that it is a decreasing function of γ suggesting the benefits the community gets by increasing levels of treatment. Is vaccination necessary in a community where treatment is there? To answer this question, we partially differentiate ℛPT with respect to α and we obtain
(16)∂ℛPT∂α=-μℛ0(μ+ρ)(1-δ)(μ+γ)(μ+ρ+α)2<0.
The fact that (16) is negative suggests that even in communities where treatment of scabies is already available the community will also benefit scabies vaccination.
4. Numerical Simulations
In this section, we carry out detailed numerical simulations using MATLAB programming language to assess the possible effects intervention strategies will have on scabies control. The parameter values that we use for numerical simulations are in Table 3. In Table 3, CSOZ denotes Central Statistical Office of Zimbabwe. For vaccination and vaccine waning rates, we used estimates used in possible TB vaccines by Bhunu et al. [5], just for illustration purposes.
Model parameters.
Model parameter
Value (range)
Source
Λ
0.029 yr−1 × 10000
CSOZ
μ
0.02 yr-1
CSOZ
β
0.043 yr-1
Landwehr et al., 1998 [4]
δ
0.125 (0.0-1.0)
Assumed
σ
0.25 (0.0-1.0)
Assumed
α
0.8 yr-1
Bhunu et al., 2008 [5]
ρ
0.0002 yr-1
Bhunu et al., 2008 [5]
γ
0.75 yr-1
Assumed
Figure 3 shows the effect of vaccination and treatment on the actual scabies cases and on the reproduction numbers obtained for different scenarios. This in a way allows comparison of vaccination and treatment as intervention strategies to control scabies. Figure 3(a) is graphical representation of scabies infected population under different intervention strategies. It shows that in the absence of any intervention strategy the scabies infected population increases from the onset and reaches a state where it remains constant. In the presence of vaccination, only the number of scabies increases from the onset until it reaches its asymptotic state where it remains more or less constant. It is worth mentioning here that trend A is always greater than trend B apart from the starting point. This suggests that vaccination will partially reduce scabies infection. However, the situation is different when treatment is involved as the only intervention strategy, and the number of scabies case declines to asymptotically low levels, suggesting that treatment alone can control scabies infections. In the case where vaccination of susceptibles is coupled with treatment infectives, the scabies infected population falls off more rapidly to asymptotically low levels than when treatment is used as the only intervention strategy. Vaccination reduces the rate at which individuals are infected and treatment kills off scabies causing mites. Figure 3(b) illustrates relationship and behavior of ℛ0,ℛP,ℛT,ℛPT for varying rate of acquiring scabies infection. It shows that ℛ0>ℛP>ℛT>ℛPT for the given parameter values and when there is one scabies infected individual in a pool of exclusively susceptible individuals. This implies that vaccination coupled with treatment is the most effective way of combating scabies in a community. Comparing results in Figure 3(a) with those in Figure 3(b), we conclude that vaccination coupled with treatment is the best approach to control scabies as evidenced by the depletion of scabies infected cases in Figure 3(a) and a very small reproduction number in Figure 3(a). By looking at trends C and D (Figure 3(a)) and ℛTandℛPT (Figure 3(b)), we safely also conclude that treatment alone is also effective in scabies control.
Effects of vaccination and treatment. A denotes when there is no intervention, B denotes when vaccination is the only intervention strategy to fight scabies, C denotes when scabies treatment is the only intervention, and D denotes when vaccination and treatment are both used to control scabies. Parameter values used are in Table 3.
In Figure 4, effects of varying initial conditions when the vaccination induced reproduction is greater than unity are illustrated. It shows that the susceptibles and scabies infected converge to their corresponding endemic equilibrium states. It is worth mentioning that there is no recovered population in this case as once an individual has been infected, he/she remains so in the absence of treatment.
Simulations of model system (1) showing plots of the susceptibles and scabies infected as a function of time with various initial conditions when ℛP>1(ℛP=1.545). Parameter values used are in Table 3.
Figure 5 graphically shows what happens when ℛPT<1 for various initial conditions. It shows that all the sub-populations of people converge to their respective disease-free equilibrium states. This signifies that making the reproduction number less than unity is necessary for disease eradication.
Simulations of model system (1) showing plots of the susceptibles and scabies infected as a function of time for various initial conditions when ℛPT<1(ℛPT=0.875). Parameter values used are in Table 3.
5. Discussion
A mathematical model for the transmission dynamics of scabies has been proposed and investigated. The potential of a possible preexposure vaccine and treatment have been analysed and compared in an effort to understand the best possible scabies control strategies. The disease-free equilibrium is shown to be locally asymptotically stable when the corresponding reproduction number is less than unity. The possible endemic equilibria states of the model have been explored using the Descartes' rule of signs. Analysis of the reproduction number has shown that vaccination has some potential to control the scabies epidemic. Additionally, it has been shown that, in communities where scabies treatment is available, adding vaccination into that community will increase the potential to control scabies. It has been shown that to effectively control scabies through vaccination alone there should be high vaccination coverage rates with a vaccine which does not wane too quickly. Numerical simulations also support the analytical results obtained that is vaccination coupled with treatment is the best possible approach to successfully fight scabies in a given community. Furthermore, numerical simulations show that treatment alone is able to control scabies as noted by the depletion of scabies cases and low reproduction number when it is used as the only intervention strategy. This is good news considering that currently there is no preexposure vaccine against scabies. However, given the fact that the most effective medications for scabies are beyond the reach of many, a development and deployment of a “new” vaccine will help in controlling the spread of scabies especially among the disadvantaged communities.
McCarthyJ. S.KempD. J.WaltonS. F.CurrieB. J.Scabies: more than just an irritation2004809453823872-s2.0-324270852010.1136/pgmj.2003.014563TielschJ. M.BeecheA.Impact of ivermectin on illness and disability associated with onchocerciasis200494A45A562-s2.0-234265221610.1111/j.1365-3156.2004.01213.xMellanbyK.Scabies in 1976197797132402-s2.0-0347707160LandwehrD.KeitaS. M.PönnighausJ. M.TounkaraC.Epidemiologic aspects of scabies in Mali, Malawi, and Cambodia19983785885902-s2.0-003186726310.1046/j.1365-4362.1998.00367.xBhunuC. P.GariraW.MukandavireZ.MagombedzeG.Modelling the effects of pre-exposure and post-exposure vaccines in tuberculosis control200825436336492-s2.0-5134912384410.1016/j.jtbi.2008.06.023TaplinD.MeinkingT. L.ChenJ. A.SanchezR.Comparison of crotamiton 10% cream (Eurax) and permethrin 5% cream (Elimite) for the treatment of scabies in children19907167732-s2.0-0025101997ScabiesA. J. O.1984New York, NY, USASpringerGrayA. M. H.Treatment of scabies194153148van den HoekJ. A.van de WeerdJ. A.BaayenT. D.MolenaarP. M.SonderG. J.van OuwerkerkI. M.de VriesH. J.A persistent problem with scabies in and outside a nursing home in Amsterdam: indications for resistance to lindane and ivermectin20081348562-s2.0-60849129151MellanbyK.The development of symptoms, parasitic infection and immunity in human scabies194435197206ArlianL. G.RappC. M.MorganM. S.Resistance and immune response in scabies-infested hosts immunized with Dermatophagoides mites19955265395452-s2.0-0029011927ArlianL. G.MorganM. S.Vyszenski-MoherD. L.StemmerB. L.Sarcoptes scabiei: the circulating antibody responses and induced immunity to scabies199478137502-s2.0-002830542710.1006/expr.1994.1004ShrankA. B.AlexanderS. L.Scabies: another epidemic?196715416696712-s2.0-0014210051van den DriesscheP.WatmoughJ.Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission200218029482-s2.0-003684527410.1016/S0025-5564(02)00108-6