Global elimination of leprosy by 2020: are we on track?
© Blok et al. 2015
Received: 30 August 2015
Accepted: 3 October 2015
Published: 22 October 2015
Every year more than 200,000 new leprosy cases are registered globally. This number has been fairly stable over the past 8 years. WHO has set a target to interrupt the transmission of leprosy globally by 2020. The aim of this study is to investigate whether this target, interpreted as global elimination, is feasible given the current control strategy. We focus on the three most important endemic countries, India, Brazil and Indonesia, which together account for more than 80 % of all newly registered leprosy cases.
We used the existing individual-based model SIMCOLEP to predict future trends of leprosy incidence given the current control strategy in each country. SIMCOLEP simulates the spread of M. leprae in a population that is structured in households. Current control consists of passive and active case detection, and multidrug therapy (MDT). Predictions of leprosy incidence were made for each country as well as for one high-endemic region within each country: Chhattisgarh (India), Pará State (Brazil) and Madura (Indonesia). Data for model quantification came from: National Leprosy Elimination Program (India), SINAN database (Brazil), and Netherlands Leprosy Relief (Indonesia).
Our projections of future leprosy incidence all show a downward trend. In 2020, the country-level leprosy incidence has decreased to 6.2, 6.1 and 3.3 per 100,000 in India, Brazil and Indonesia, respectively, meeting the elimination target of less than 10 per 100,000. However, elimination may not be achieved in time for the high-endemic regions. The leprosy incidence in 2020 is predicted to be 16.2, 21.1 and 19.3 per 100,000 in Chhattisgarh, Pará and Madura, respectively, and the target may only be achieved in another 5 to 10 years.
Our predictions show that although country-level elimination is reached by 2020, leprosy is likely to remain a problem in the high-endemic regions (i.e. states, districts and provinces with multimillion populations), which account for most of the cases in a country.
Leprosy or Hansen disease is caused by an infection of Mycobacterium leprae, usually acquired through contact with an infected person. However, not everyone exposed to an infected contact will eventually develop the disease . Worldwide, more than 200,000 new leprosy cases are detected annually . This number has been fairly stable in the past 8 years, indicating ongoing transmission. In 2013, 14 countries reported more than 1000 new cases, of which three countries - India, Brazil and Indonesia - account for more than 80 % of all the cases in the world . The distribution of leprosy is becoming localized to a limited number of countries . Also within countries leprosy is found to be spatially unevenly distributed [3, 4].
The spread of leprosy in a population is highly dependent of the variation in susceptibility of individuals and the intensity of contact with other individuals [5, 6]. As a result leprosy patients are often found to be spatially clustered in regions, neighborhoods, families and households [7–9]. Contacts closest to the index patient, in particular household contacts, have the highest risk [5, 8].
Leprosy patients are treated with a combination of antibiotics known as multidrug therapy (MDT), which effectively heals the patient and thereby also reduces the infectivity in the community. However, MDT alone is not sufficient to prevent new cases. Therefore the main principles of leprosy control also include early detection of leprosy cases . Other methods of prevention such as immunoprophylaxis and chemoprophylaxis against M. leprae are not yet widely available .
Global elimination of leprosy has been a target for a long time. In 1991, the 44th World Health Assembly already adopted the target of elimination as a public health problem, defined as reducing the prevalence to less than 1 case per 10,000 population by the year 2000 . Although this target was met at a global level, elimination of leprosy as a public health problem has not been achieved in some endemic countries, in particular at a subnational level . Recently, WHO has formulated a roadmap for 17 neglected tropical diseases, including leprosy, to reduce their global impact. The targets for leprosy are (1) global interruption of transmission or elimination by 2020, and (2) reduction of grade-2 disabilities in newly detected cases to below 1 per million population at global level by 2020 .
The aim of this study is to investigate whether the first target is feasible given the current control strategy. To be able to assess this, we interpreted the first target of 'interruption of transmission' in terms of incidence reduction, contrary to the original WHO elimination target (as public health problem), which was defined in terms of prevelance reduction. Our operational definition of global elimination of leprosy as applied in this paper is: 'less than 10 new cases per 100,000 population. We focus on the three countries with the highest number of annual new cases: India, Brazil and Indonesia. Since leprosy is spatially unevenly distributed within countries, we will for comparison also make predictions for one high-endemic region within each country: Chhattisgarh (India), Pará State (Brazil) and Madura Island (Indonesia). We will use the existing individual-based model SIMCOLEP, which simulates the transmission and control of leprosy in a population structured by households [6, 13]. This model has been previously quantified to the leprosy situation in northwest Bangladesh with the aim of making future predictions of leprosy trends and testing the impact of various interventions .
In this study, we apply SIMCOLEP to make predictions of future leprosy incidence in India, Brazil, Indonesia and aforementioned regions until the year 2030. For each of these countries and regions, we will assess whether elimination, defined as less than 10 new cases per 100,000, will be met by 2020 given continuation of the current control strategies of MDT treatment and early detection.
The individual-based model SIMCOLEP simulates the spread of M. leprae in a population structured in households. Dynamics of the population are described by births, movements between households, and deaths. Births are determined by birth rates and newborns are placed in the household of their mother. Movements of individuals to another newly created or existing household occur after marriage, during adolescence or after becoming a widow(er). Deaths are determined by death rates .
In the model, transmission of M. leprae occurs after direct contact with an infectious individual. Two transmission processes are modelled separately: transmission in the general population and a within-household transmission. Infectivity is determined by the product of the contact rate and the probability of infection during a contact. Each transmission process has its own contact rate: general population, c pop , and within households, c hh . An infected individual will develop either paucibacillary leprosy (PB) or multibacillary leprosy (MB). Both types can be detected, treated and cured, but only MB leprosy is considered infectious in the model. The natural history of leprosy is modelled following the model of Meima et al. .
Leprosy can only be acquired by susceptible individuals. It is assumed that the majority of the population is never susceptible to leprosy . Since the mechanism that underlies susceptibility is still unknown and could not be identified in previous (modelling) work, susceptibility of an individual is randomly determined at birth. We assumed that 20 % of the population is susceptible (i.e. 80 % is never susceptible), based on model fitting results of a previous study . The type of leprosy (i.e. MB or PB) is randomly determined based on the distribution of the type of leprosy in each country. The MB proportion is about 48 %, 66 %, and 83 % in India, Brazil, and Indonesia, respectively.
Model fitting and data
Overview of data used to quantify the model
Census Indonesia (BPS)
1990, 2000, 2010
1993, 1997, 2000, 2007
1990, 2000, 2012
1991, 2000, 2010
1990, 2000, 2009
1992, 1998 , 2006, 2011, 2012, 2013
DHSa, Census India
1991, 1995, 2000, 2005, 2010
1971, 1980, 1990, 1991, 1997, 2003, 2007, 2012
1990, 2000, 2010
Distribution of household size
1993, 1999, 2006
1991, 1994, 1997, 2003, 2007
New case detection rate
- High-endemic region
After fitting the population and household structure in each country, we fitted the simulated leprosy trends to the observed trends in each country and region separately. Data used to fit leprosy trends include the new case detection rate and the MB proportion. Indian data were obtained from the National Elimination Program (NLEP) . Brazil leprosy data came from the Sistema de Informações de Agravos de Notificação (SINAN) database . The SINAN database is the Brazilian national database for communicable diseases. Finally, leprosy data for Indonesia were provided by the Netherlands Leprosy Relief (NLR) foundation (see Table 1).
In order to fit the leprosy trends, we also modeled leprosy control programs of the past decades in each country. The leprosy control program consisted of treatment with Dapsone until 1989 and MDT afterwards, and passive and active case detection. Changes in passive detection were expressed in terms of detection delays. These detection delays were fitted to match the trends of the new case detection rates. Based on literature, most recent detection delays were set to 2–3 years in our model . The estimated detection delays gradually improved since 1970 from 13 to 2 years, 18 to 3 years, and 10 to 2 years (gamma distributed) in India, Brazil and Indonesia, respectively. Active case detection, defined as examining contacts of the patient, was only included if data about coverage rates were available. In Indonesia active detection was included since 2010 with a coverage rate of 9–11 %. In Brazil the coverage rate was between 43 and 59 %, and started since 2003 . No active case detection was included in India . BCG vaccination in infants, which has a protective effect of 60 %, was also included . Coverage rates of BCG vaccination in infants, obtained from WHO, were used . We further assumed that the national leprosy control strategy was implemented consequently within regions.
Leprosy new case detection rates were fitted to the data by calibrating contact rates. In this study, we only calibrated the contact rate in the general population. Since no data was available about the prevalence and/or incidence of leprosy by household size and per household member, we fixed the contact rate within household to the optimal value of previous work (c hh = 0.98) . We assumed that contact rates within households would not differ between countries and regions. The contact rate in general population (c pop ) of each country and region was estimated separately by fitting the simulated new case detection rates to the observed new case detection rates from the data (see Table 1). To this end, the simulated new case detection rates (mean of 100 runs) were compared to the data by a log-likelihood function assuming a Poisson distribution. These likelihood ratios were fitted to a polynomial regression meta-model to obtain the optimal c pop value.
New simulations were performed with the optimal c pop values. A detailed description of the fitting procedure can be found elsewhere . A manual of the model, the model itself and the input files with country-specific quantifications are provided in Additional file 2 and Additional file 3.
Using the best fit for each country and high-endemic region, we made future predictions assuming continuation of current leprosy control programs: MDT treatment, passive detection and active detection (only Brazil and Indonesia). Each simulation run starts in the year 1000 with a modelled population size of 20,000 individuals that remains constant until 1850–1900. Afterwards the population increase with the annual growth rate to more than 125,000 individuals in 2030. All simulations continued to the year 2030 to predict future trends and to determine whether elimination could be reached by 2020. Simulation results were an average of 100 runs.
The new case detection rates in the regions Chhattisgarh, Pará State and Madura are 2–7 times higher than the national rates. Similarly to country-level predictions, the new case detection rates in the high-endemic regions will continue to drop over time. However, in these regions the elimination target will not be met in time. In 2020, the predicted new case detection rates are 16.2 (90 % CI: 9.4-24.7), 21.1 (90 % CI: 16.9-24.2), and 19.3 (90 % CI: 15.9-23.3) per 100,000 in Chhattisgarh, Pará State, and Madura, respectively. Elimination in these regions will not be established before 2026. The annual decrease between 2015 and 2030 is about 8.7 %, 8.4 %, and 7.2 % in Chhattisgarh, Pará State, and Madura, respectively. Predictions in terms of annual new cases are shown in Additional file 1: Figure S3.
The distribution of new leprosy cases by age group differs by country and is affected by control strategies. Improved detection will mainly increase the detection of older individuals (See Additional file 1: Figure S4). Also the distribution of new leprosy cases by household size differs by country. In Brazil most cases are found in bigger households, while in Indonesia most cases are found in smaller households (See Additional file 1: Figure S5).
We used the individual-based model SIMCOLEP to assess whether the WHO target of elimination of leprosy, defined as less than 10 new cases per 100,000 annually, would be met by 2020. We focused on three high-endemic countries, India, Brazil and Indonesia, which account for more than 80 % of all cases worldwide. Our study shows that elimination of leprosy in these countries has already been met or will be met in the very near future at country-level. However, in the high-endemic regions of these countries – Chhattisgarh, Pará State and Madura – elimination of leprosy will not be achieved by 2020 with the current leprosy control strategy. Regional new case detection rates are 2–7 times higher than at country-level.
Country-level predictions of the new case detection rates, which suggest elimination, provide a biased view of the leprosy situation, because these rates are masked by the large population size of each country. Focusing on high-endemic regions where this elimination target has not yet been achieved will give a more realistic representation of the actual leprosy situation in a country. It also reflects more accurately that the distribution of leprosy is becoming more localized [3, 4, 23].
Our results show a continuous decline in new case detection rates, demonstrating that the present leprosy control strategy, if consistently applied, will reduce transmission. Results also show that confidence bounds become smaller around future projections. This is caused by the downward trend of the NCDR given the control strategy and partially by the population growth. Predictions in terms of annual number of new cases illustrate that confidence bounds decrease less rapidly when the population growth is not considered (See Additional file 1: Figure S3 ). The annual rate of decline is greatest in India, followed by Brazil and Indonesia. Also within countries the speed of decline differs. This is caused by differences in intensity of contact, reflecting different living conditions based on varying socioeconomic and cultural circumstances. In high-endemic regions, which are generally less developed, the calibrated contact rate in the general population was higher than at country-level country rates, which reflects a nation-wide average.
The poorer fit of Madura in Indonesia might be explained by regions specific operational changes in control. Madura is an island in Indonesia and therefore relatively isolated. It has always been a high-endemic region. The poor fitting is probably the result of the focused intensified case finding activities in that particular area between 2003 and 2005 causing the NCDR to double in that region, which our model did not account for [24, 25].
The large number of undetected cases remains a threat to the elimination of leprosy generally. Leprosy is a slow disease with a long incubation time and long delays in detection [20, 26, 27]. These missing cases still contribute to the ongoing transmission. It has been estimated that over 4 million cases are missed between 2000 and 2020 worldwide . It has been argued that in India half of the cases have not been reported to meet the elimination targets of 2005 . Since our model heavily depends on the available case data, including their imperfections, the problem of missing cases is also inherent to our predictions. On the whole, the actual number of new leprosy cases is likely to be higher than presented in our predictions [20, 27, 28].
A concern of this study is that the predictions are made under the assumption that being susceptibility to leprosy was fully randomly determined. Previous work has shown that although this may be a valid assumption, other mechanisms, such as genetic inheritance of susceptibility, might also explain variation in susceptibility . Assuming another mechanism for susceptibility has proven to somewhat slow down the speed of decline of the new case detection trends . Our results can therefore be regarded as a best case scenario.
In our study elimination was defined as less than 10 new cases per 100,000, while elimination is officially defined as a prevalence rate of less than 10 per 100,000 population. This will, however, not alter our conclusions, because the prevalence rate of leprosy does not differ much from the new case detection rate.
Furthermore, our predictions are based on a continuation of the present leprosy control strategy, which also includes some contact tracing in Brazil and Indonesia. It assumes full adherence to the strategy. Any unexpected adverse events, such as famines and social upheaval, have not been considered. On the other hand, our study also does not account for possible additional leprosy control practices in regions, such as chemoprophylaxis for contacts of leprosy cases.
An important next step is to determine which interventions at population level would bring elimination forward and have the highest impact on future incidence of disease through the interruption of transmission. The focus should be on interventions targeting the contacts of leprosy patients to prevent new cases . Examples of such intervention are intensive contact tracing, administering chemoprophylaxis or immunoprophylaxis (e.g. BCG-like leprosy vaccine) and early diagnosis of leprosy using diagnostic tests for infection or tests that predict clinical disease.
Although it seems that country-level elimination is reached by 2020, leprosy still remains a problem in the high-endemic regions, which account for most of the cases in a country. These regions often have a multimillion population in which rapid progress of leprosy control, even if conducted optimally, will not be achieved soon. We therefore conclude that ongoing transmission of M. leprae will make global elimination of leprosy as a unlikely to occur by 2020 without further control measures.
We gratefully acknowledge funding of the NTD Modelling Consortium by Novartis Foundation. The views, opinions, assumptions and conclusions or any other information set out in this article are solely those of the authors and not of Novartis Foundation. Novartis Foundation was not involved in the establishment of this article and this article should not be attributed to Novartis Foundation or any person connected with Novartis Foundation. We are grateful to the Netherlands Leprosy Relief (NLR) foundation for providing the Indonesian leprosy data.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Fine PE. Leprosy: the epidemiology of a slow bacterium. Epidemiological Review. 1982;4:161–88.Google Scholar
- WHO. Global leprosy update, 2013; reducing disease burden. Wkly Epidemiol Rec. 2014;89(36):389–400.Google Scholar
- Penna ML, Wand-Del-Rey-de-Oliveira ML, Penna G. Spatial distribution of leprosy in the Amazon region of Brazil. Emerg Infect Dis. 2009;15(4):650–2. doi:https://doi.org/10.3201/eid1504.081378.PubMed CentralView ArticlePubMedGoogle Scholar
- Fischer E, Pahan D, Chowdhury S, Richardus J. The spatial distribution of leprosy cases during 15 years of a leprosy control program in Bangladesh: an observational study. BMC Infect Dis. 2008;8:126. doi:https://doi.org/10.1186/1471-2334-8-126.PubMed CentralView ArticlePubMedGoogle Scholar
- Moet FJ, Pahan D, Schuring RP, Oskam L, Richardus JH. Physical distance, genetic relationship, age, and leprosy classification are independent risk factors for leprosy in contacts of patients with leprosy. J Infect Dis. 2006;193(3):346–53. doi:https://doi.org/10.1086/499278.View ArticlePubMedGoogle Scholar
- Fischer EAJ, De Vlas SJ, Meima A, Habbema DF, Richardus JH. Different mechanisms for heterogeneity in leprosy susceptibility can explain disease clustering within households. PLoS One. 2010;5(11):e14061. doi:https://doi.org/10.1371/journal.pone.0014061.PubMed CentralView ArticlePubMedGoogle Scholar
- Bakker MI, May L, Hatta M, Kwenang A, Klatser PR, Oskam L, et al. Genetic, household and spatial clustering of leprosy on an island in Indonesia: a population-based study. BMC Med Genet. 2005;6:40. doi:https://doi.org/10.1186/1471-2350-6-40.PubMed CentralView ArticlePubMedGoogle Scholar
- Fine PE, Sterne JA, Ponnighaus JM, Bliss L, Saui J, Chihana A, et al. Household and dwelling contact as risk factors for leprosy in Northern Malawi. Am J Epidemiol. 1997;146(1):91–102.View ArticlePubMedGoogle Scholar
- Moet FJ, Meima A, Oskam L, Richardus JH. Risk factors for the development of clinical leprosy among contacts, and their relevance for targeted interventions. Lepr Rev. 2004;75(4):310–26.PubMedGoogle Scholar
- Richardus JH, Oskam L. Protecting people against leprosy: chemoprophylaxis and immunoprophylaxis. Clin Dermatol. 2015;33(1):19–25. doi:https://doi.org/10.1016/j.clindermatol.2014.07.009.View ArticlePubMedGoogle Scholar
- WHO. Elimination of leprosy: resolution of the 44th World Health Assembly (Resolution No. WHA 44.9). Geneva: WHO; 1991.Google Scholar
- WHO. Accelerating work to overcome the global impact of neglected tropical diseases - a roadmap for implementation. Geneva: World Health Organization; 2012.Google Scholar
- Blok DJ, de Vlas SJ, Fischer EA, Richardus JH. Mathematical modelling of leprosy and its control. Adv Parasitol. 2015;87:33–51. doi:https://doi.org/10.1016/bs.apar.2014.12.002.View ArticlePubMedGoogle Scholar
- Fischer EA, de Vlas SJ, Habbema JD, Richardus JH. The long-term effect of current and new interventions on the new case detection of leprosy: a modeling study. PLoS Negl Trop Dis. 2011;5(9):e1330. doi:https://doi.org/10.1371/journal.pntd.0001330.PubMed CentralView ArticlePubMedGoogle Scholar
- Meima A, Gupte MD, Van Oortmarssen GJ, Habbema JDF. SIMLEP: A simulation model for leprosy transmission and control. Int J Lepr Other Mycobact Dis. 1999;67(3):215–36.PubMedGoogle Scholar
- Programme) NNLE. NLEP - Progress report for the year 2013–14. New Delhi: Central Leprosy Division2014.Google Scholar
- Sistema de Informações de Agravos de Notificação (SINAN). Minstry of Health. http://dtr2004.saude.gov.br/sinanweb/index.php.
- Deps PD, Guedes BVS, Filho JB, Andreatta MK, Marcari RS, Rodrigues LC. Delay in the diagnosis of leprosy in the Metropolitan Region of Victoria, Brazil. Lepr Rev. 2006;77:41–7.PubMedGoogle Scholar
- Ministry Of Health. National Hansen’s Disease Control Programme Report 2009–2010: Ministry of Health 2011.Google Scholar
- Lockwood DN, Shetty V, Penna GO. Hazards of setting targets to eliminate disease: lessons from the leprosy elimination campaign. BMJ. 2014;348:g1136. doi:https://doi.org/10.1136/bmj.g1136.View ArticlePubMedGoogle Scholar
- Schuring RP, Richardus JH, Pahan D, Oskam L. Protective effect of the combination BCG vaccination and rifampicin prophylaxis in leprosy prevention. Vaccine. 2009;27(50):7125–8. doi:https://doi.org/10.1016/j.vaccine.2009.09.054.View ArticlePubMedGoogle Scholar
- WHO. Reported estimates of BCG coverage. World Health Organization. http://apps.who.int/immunization_monitoring/globalsummary/timeseries/tscoveragebcg.html. 2015.
- Brook CE, Beauclair R, Ngwenya O, Worden L, Ndeffo-Mbah M, Lietman TM et al. Spatial heterogeneity in projected leprosy trends in India Parasit Vectors. 2015. doi:https://doi.org/10.1186/s13071-015-1124-7.
- van Brakel WH. Technical Support Visit Indonesia no.2 (Surabaya and Madura). Amsterdam: Royal Tropical Institute; 2008.Google Scholar
- WHO. Global leprosy situation, 2005. Geneva: World Health Organization; 2005.Google Scholar
- WHO. WHO expert committee on leprosy. WHO Tech Rep Ser. 2012;968:1–61.Google Scholar
- Smith WC, van Brakel W, Gillis T, Saunderson P, Richardus JH. The missing millions: a threat to the elimination of leprosy. PLoS Negl Trop Dis. 2015;9(4):e0003658. doi:https://doi.org/10.1371/journal.pntd.0003658.PubMed CentralView ArticlePubMedGoogle Scholar
- Crump RE, Medley GF. Back-calculating the incidence of infection of leprosy in a Bayesian framework. Parasit Vectors. 2015. doi:https://doi.org/10.1186/s13071-015-1142-5.