Factors contributing to the decline of leprosy in spain in the second half of the twentieth century
ABSTRACT
BACKGROUND: Leprosy is a chronic infectious disease that is considered to be declining, though it still remains prevalent in many parts of the world. A study was made to explore the health and socioeconomic factors that most influenced the trend of the disease in a typical Mediterranean country.RÉSUMÉ
CONTEXTE: La lèpre est une maladie chronique considérée comme étant sur le déclin, bien qu'elle reste encore prévalente dans plusieurs régions du monde. Le but de cette étude fut d'explorer dans un pays méditerranéen typique, les facteurs sanitaires et socio-économiques qui ont eu le plus d'impact sur l'évolution épidémiologique de la maladie.RESUMEN
BACKGROUND: Aunque la lepra es una enfermedad infecciosa crónica en decaimiento, aún permanece vigente en muchas partes del mundo. El presente, es un estudio realizado con el fin de explorar los factores de salud y sanitarios que tienen más influencia en la tendencia de la enfermedad en un típico país Mediterráneo.Leprosy is a chronic infectious disease produced by Mycobacterium leprae that results from an imbalance between the bacterium and the host immune response. With a history that goes back three thousand years, leprosy is one of the oldest diseases known.
For centuries, the care of leprosy patients lacked a scientific basis. Affected individuals were isolated from society and from their couples, and diseased children were separated from their families. Later, scientific evidence came to show that in some cases isolation is not particularly effective, and that no more than 5% of exposed individuals actually develop the disease.
In its third regional conference of leprosy, the World Health Organization (WHO) proposed elimination of the disease as a health problem in all countries by the year 2005. However, the number of new cases documented by the WHO in 2003 totaled 513,798, while the prevalence of the disease that same year was 457,792 cases. In other words, leprosy is still far from being controlled, and work remains to be done before the disease can be relegated to the history books. The majority of cases of leprosy are found in Southeast Asia. The present annual incidence of the disease worldwide is almost one new patient per minute, and is fundamentally attributable to developing countries like India-with 70% of all such cases. Some studies have even reported the incidence of the disease to be increasing in these parts of the world (25).
One of the principal factors relating to control of the disease, and which was already suspected many years ago, is the social and economical development of society. In this context, it has been shown that economical development is accompanied by important improvements in patient quality of life (QoL), and a decrease in mortality (3).
Although the incidence of leprosy in Spain has decreased greatly, thanks to factors such as economical development, and the prevalence of the disease has also decreased in recent years, a minimum number of fundamentally imported cases remain (7, 21), as a result of which it is difficult to speak of actual eradication of the disease.
The present study explores the influence of socioeconomic factors in the evolution of leprosy in Spain-a country illustrative of the Mediterranean setting of the disease in the last 50 years.
MATERIALS AND METHODS
Two leprosy information sources have been used: the Spanish State Leprosy Registry (Registro Estatal de Lepra), a database of the National Epidemiology Center (Centro Nacional de Epidemiología) (6), located in the Instituto de Salud Carlos III in Madrid, and which came into operation in 1991; and the Annual Statistical Bulletin (Anuario Estadístico) of the National Institute of Statistics (Instituto Nacional de Estadística, INE) (11) for prior data up until the year 1950.
Economic (national income) data and their provincial distribution were in turn obtained from the publications of the Service of Studies of the Banco de Bilbao Vizcaya (1) (presently BBVA). Demographic information was collected from the Natural Variations of the Spanish Population (Movimiento Natural de la Población Española) (12), and from the Spanish population census (13), both being publications of the INE. The latter document reflects the data corresponding to the population censuses conducted every decade.
The variables derived from these information sources were classified into four groups. Variables to be included were required to remain stable in terms of both the information sources and the way in which they were obtained during the 50-year period covered by the study.
The data were obtained at the provincial level, followed by grouping according to the different Spanish Autonomous Communities. The selected socio-demographic variables were: the proportion of rural population (i.e., the percentage of the population not living in the capital versus the total population), the employment rate with respect to the total population (per thousand), the emigration rate, the schooling rate (defined as the proportion of the population in the 614 years age range in school versus the total population), and the unemployment rate.
As economical variables, the gross domestic product (GDP) per inhabitant was recorded, along with the product per educational and health sectors. The magnitudes of the economical variables were determined for each year at constant price, using the GDP deflactor, and employing the 1986 prices as reference (expressed in euros). The variable corresponding to health care resources was the proportion of physicians per 1000 inhabitants, and the health variables were the infant mortality rate, the incidence of leprosy per 100,000 inhabitants, and the incidence of respiratory tuberculosis (TBC) per 1000. Base 100 index numbers were used for both the GDP and the incidence of leprosy and TBC. Retrospectively, the 5-year elasticity was calculated based on the ratio between the 5-year incidence variations and the economical variations corresponding to the same period.
The statistical analyses included, a complete study of the type of distribution, followed by simple and multiple regression, residues analysis with the Kolmogorov-Smirnov test for one sample, and correlation between typified residues. Likewise, the colinearity between variables was studied based on eigenvalues, the condition index and proportions of variance. Finally, fitting was carried out based on logarithmic transformation of the data and colinearity study via the tolerance, variance inflation factor (VIF), and also the study of residues-including histograms of typified residues and dispersion plots.
In practice, exact colinearity is rarely observed in studies of this kind, though socalled "near colinearity" is seen with some frequency (i.e., some variable is found to be "nearly" a linear combination of another variable/s). We have jointly applied the method proposed by Kleimbaum (18) and the criteria of Belsley (2). The latter author proposes utilization of the condition indices and proportion of variance decomposition to conduct the analysis of colinearity-establishing a value of 0.5 as the upper proportion threshold, as a result of which the conclusions are finally reached as follows: high condition indices (over 30) indicate the number of colinearities, and their magnitude reflects their relative importance. If a given component possesses a condition index of over 30, and two or more variables present a high proportion of variance in the latter, then these variables are considered to be colinear. Taking into account the long incubation period of leprosy (generally between 210 years) (26), the incidence rates were studied with the socioeconomic status of the population 110 years previously, i.e., at the time of contagion.
RESULTS
Regarding the mean annual incidence rates for leprosy, established for 5-year periods and 100,000 inhabitants during the period 19502000 (Fig. 1), and the interannual variations (Table 1), important differences are seen among the different regions. In effect, while most regions show reductions at a rate of 1% annually, some actually end up showing increases-such as the Balearic Islands, Navarra, and the Basque Country. Andalucía, in the south of the country, for many years exhibited the highest leprosy incidence for that period, though it is also the Spanish Autonomous Community with the greatest interannual reductions (almost 4.5%), followed by Galicia-with reductions of 3%-and the Canary Islands and Valencia (2%). Overall, Spain showed a decrease of 1.6% annually over this period of 50 years.
Fig. 1. Mean 5-yearly incidence rates for leprosy in Spain in the second half of the Twentieth Century. Note: Incidence rate = cases per 100,000.
In Table 2 the longitudinal study by regions presents the mean values from 1950 to 2000 for the variables considered in the study. Differences between regions are clearly seen, with important variability in both the economical values (with a nearly threefold difference between minimum and maximum), and particularly the incidence of tuberculosis (with an up to 5.5fold difference) and leprosy-the latter being the parameter showing the largest differences (up to 25-fold). In relation to this latter variable, the regions with the highest figures were Andalucía and Extremadura, along with Catalonia and the Canary Islands.
On the other hand, in relation to the 5-year means of the variables considered in the study by regions (Table 3), all Autonomous Communities were seen to exhibit a relatively regular evolution with gradual decreases or increases-the sole exception being emigration, which proved to be somewhat irregular.
In reference to 100 base index numbers of the Gross Domestic Product (GDP) and incidences for leprosy and tuberculosis in Spain (Table 4), a 7.23-fold economical increase was observed since the year 1950- this figure being similar to the reduction in the incidence of leprosy in the same period of time (7.4-fold). In turn, tuberculosis was seen to have decreased in incidence 3.7fold in the course of the period 19502000. On the other hand, the mean annual variation in GDP was seen to be 14%, while the mean variations for leprosy and tuberculosis were quite similar (1.66 and 1.57, respectively). The proportional variations taking economical growth as reference base revealed that leprosy varied 0.06 units for each proportional unit economical increment. The situation was similar in the case of tuberculosis (0.07).
Considering that economical status exerts an influence at the time of actual contagion, even though clinically manifest leprosy only appears some years later, we evaluated the data on economic growth corresponding to the preceding years (1, 2, 3, etc.) in relation to the incidence of Hansen’s Disease. To this effect we developed both simple and multiple regression models, using the annual information from all Spanish provinces. In this respect, Table 5 shows the results of multiple regression, fitted by Naperian logarithmic transformation. The explanatory contribution, measured in terms of R2, is seen to increase until a 6-year lag was reached; at this point 57% of the incidence of leprosy was explained. Later, with increasing delays, the values began to decrease- though in an irregular manner.
In the multiple regression analysis, based on the model involving a 6-year lag and using the annual and provincial data, we in turn developed a complete model (Table 6). The results revealed that the factors seen to influence the incidence of leprosy were: employment, the number of physicians, and economic growth, all with negative coefficients, while tuberculosis was associated with a positive coefficient. Nevertheless, the GDP was found to yield the most important coefficient, with a value of 0.5, followed by infant mortality and the physician rate.
In view of the similar behavior of tuberculosis as regards economic growth, a study was made of the interactions and resulting colinearity. In this context, tolerance and the variance inflation factor showed medium tolerance for the variables GDP, infant mortality and health education sector product, unemployment, and the number of physicians-while high tolerance was recorded for the rest.
DISCUSSION
The main problem posed by studies of this kind is that the data corresponding to the factors analyzed must offer continuity and reliability. This prevented us from examining a large number of variables which logically should have been included in the study.
The regions with the highest rate and most total cases of leprosy were Andalucía and Extremadura (both of which showed the lowest economic levels), followed by Catalonia and the Canary Islands-possibly due to the important number of immigrants found in these autonomous communities. Within Andalucía, the province of Córdoba clearly stands out-with a mean rate for the global study period of 0.6 new cases per 100,000 inhabitants per year-together with the province of Jaén (7) . In sum, a preferentially though somewhat irregular Mediterranean distribution is observed, in coincidence with the findings of other authors, in relation to the cases of leprosy (26).
On the other hand, although the study has been carried out with delimitations in the form of provinces and/or regions, the distribution of Hansen’s disease is known to be extremely irregular-with marked intra-regional variations in Spain and also elsewhere in the world (35).
Since 1981, the World Health Organization (WHO) has recommended multi-drug therapy for the management of leprosy, and has considered it possible to eliminate the disease as a public health problem worldwide. The leprosy-eliminating strategy proposed by the WHO has been based on the hypothesis that a reduction in the prevalence of the disease based on the administration of multi-drug regimens (dapsone, clofazimine and rifampicin) entails a low reservoir of cases in the community (fewer than 1/10,000 inhabitants)-as a result of which the transmission cycle of the disease cannot be perpetuated (25, 32). The concept underlying this strategy is the socalled "basic reproduction rate, R0" (the mean number of individuals directly infected by an infectious case). If R0 is >1, an epidemic may result, while if R0 is <1, the disease may disappear. In the case of leprosy, it has been estimated that R0 is <1 when the prevalence rate is less than 1/10,000-a situation that has been applicable in Spain for many years. In our study the observed annual reduction rate was 1.6%, while some sources in the literature conclude that countries with endemic disease have annual reduction rates of 4% (23).
However, certain data and contradictions speak against the optimism of the WHO and raise doubts as to the evolution of the disease-questioning the central argument of the elimination strategy whereby the only effective contagion source is the lepromatous patient, and a reduction in prevalence therefore necessarily implies a reduction in incidence. The principal unsettled questions refer to the evidence suggesting a slightly decreasing general incidence (though not in all countries); uncertainty over the true number of individuals with occult infection; a lack of knowledge of the magnitude of the problem posed by subclinical carriers of the disease (20); and the possible existence of animal reservoirs (5, 36). Many authors consider that the disease found in armadillos and monkeys is different from that seen in humans (27). Another important consideration is the present lack of an effective vaccine (10).
In this context, and in an attempt to account for the evolution of leprosy, sociosanitary and economic considerations have been postulated to be of greater relevance than multitherapy. In favor of this hypothesis would be the declining trend of the disease in China and Spain starting before implementation of the WHO multi-drug strategy, and the spontaneous disappearance of the disease in Norway (24), Japan (15)and Hawaii (37), as well as its continued persistence in India, Brazil (17) and Cuba despite application of the WHO program. In the case of Norway, the gradual decrease in leprosy incidence was accompanied by an increase in mean patient age at appearance of the disease, and a gradual prolongation of the incubation period (24).
Leprosy prevalence data were not used in this study, because they posed concordance problems with the incidence data. Such discordance was mainly due to the repetition of individuals included in the prevalence lists. Moreover, while the incidence rate was not affected as a result, the prevalence rates were falsely elevated due to persistence in the prevalence lists of cases that had likely died or healed.
Although the regression analysis cannot consider only the sum of the effects of variables, since the interference factors absorb part of the explanation or prediction of changes in the incidence of leprosy, the interactions among variables do not seem to be sufficiently intense to compromise interpretation of the results. The same applies to the colinearity among variables assessed by variance inflation factors and tolerance. The confounding variables and resulting attributable bias always pose a considerable risk and are a source of imprecision in ecological studies, due to the scant control afforded by the use of population instead of individual measures (29).
Although some authors have reported tuberculosis to be antagonistic to leprosy, attributing the almost total disappearance of the latter to the appearance of tuberculosis, it must be considered that this is an infectious disorder, and that the environmental conditions and factors of hygiene are therefore the bases for transmission. In this sense, socioeconomic status could be essential for creating a setting unfavorable for the transmission of both diseases, as observed in our study.
Previous authors have addressed the possible relationship between socioeconomic development and concrete health aspects, e.g., mortality in the classical study of Krishnan (19) in India, or Mckeown (22)in England; or neonatal mortality in the study of Shah and Abbey (33) in Baltimore (United States), among many others. However, a review of the literature has yielded almost no studies relating socioeconomic development to the evolution of leprosy incidence. In contrast, analyses have been made of the social and economic variations observed in leprosy patients, and the need to restore such parameters to their previous levels in order to help heal the patients, as in a published Indian study (29). Other parameters investigated include crowding, domestic conditions, family size (31), and even the inverse relation between years of schooling and the incidence of leprosy in a study from Malawi (28).
The present study made use of the gross domestic product (GDP), as recommended by some authors for conducting medium-and long-term studies (9). In this context, GDP was seen to be the factor exerting the greatest influence upon leprosy evolution, according to the multiple regression analysis with a 6-year lag period. It should be taken into account that economic development is the fundamental factor influencing family crowding conditions, schooling, and other factors such as nutritional status (33), likewise related to the incidence of leprosy. The family grouping of leprosy observed in some studies may also be a consequence of shared socioeconomic status (and therefore hygiene conditions) (16). Therefore some authors recommend considering this as a factor in chronic diseases (35).
The other two factors with negative and statistically significant coefficients were employment and the number of physicians, and the sum of their beta-coefficients was lower than that of GDP-i.e., they were seen to jointly exert much less influence upon leprosy incidence trends. The fourth significant factor was tuberculosis; in effect, since the latter is also a transmissible disease, it is also affected by social and economic development. Moreover, it should be taken into account that the results of the present study showed a similar interannual change for both transmissible diseases, and that Bacillus Calmette Guérin (BCG) vaccination produces a certain degree of immunity against leprosy as well (8, 38).
In relation to the lag period used to explain the incubation period of leprosy, the estimated period is 2.9 to 5.3 years in the case of tuberculoid leprosy, and 9.3 to 11.6 years in the case of the lepromatous form of the disease-though in our study a period of 6 years was seen to afford the best explanation for the observed statistical relationship. Other authors have reported periods of up to 30 years in veterans-a situation not examined in the present study.
It may have been advisable to separate the cases according to the type of leprosy, to establish the possible differences in the influence of socioeconomic factors according to the type of disease involved. However, the incidence rates would have been zero or very close to zero in the last years of the study period.
The declining tendency of leprosy in Spain was seen to be gradual since the 1950s, i.e., it started before the introduction of multi-drug therapy, and was more closely associated with the socioeconomic development of the country. Similar observations also apply to a number of other European countries (14), in addition to Norway as noted above (24). Our study illustrates the multicausal model of effects, whereby socioeconomic considerations can prove more important than the actual transmissibility conditions of leprosy.
REFERENCES
1. BANCO BILBAO VIZCAYA. Renta Nacional de España y su Distribución Provincial. Años 1950 a 1995. Servicio de publicaciones.
2. BELSLEY, D. Detecting and assessing the problems caused by multi-collinearity: A use of the singular-value decomposition. No. 66 in NBER Working Papers, National Bureau of Economic Research, Inc., 1974.
3. BHATTACHARYA, S. N., and SEHGAL, V. N. Leprosy elimination campaign and its impending fallout. Int. J. Dermatol. 39(9) (2000) 667669.
4. BJARNASON, O. The last lepra patient in Iceland. Nord. Medicinhist. Arsb. (1989) 197203.
5. BRUCE, S., SCHROEDER, T. L., ELLNER, K., RUBIN, H., WILLIAMS, T., and WOLF, J. Armadillo exposure and Hansen’s disease: an epidemiologic survey in southern Texas. J. Am. Acad. Dermatol. 43(2 Pt. 1) (2000) 223228.
6. CENTRO NACIONAL DE EPIDEMIOLOGÍA. Manual del registro estatal de lepra. Madrid. Instituto de Salud Carlos III. 1991.
7. DELGADO-RODRIGUEZ, M., RODRIGUEZ-CONTRERAS PELAYO, R., EXTREMERA-CASTILLO, F., SERRANO-ORTEGA, S., and GALVEZ-VARGAS, R. Epidemiologic aspects of leprosy in the province of Jaen. Rev. Clin. Esp. 185(2) (1989) 99103.
8. FINE, P. E. M. Variation in protection by BCG: implications of and for heterologous immunity. Lancet 346 (1995) 13391345.
9. GETZEN, T. E. Forecasting health expenditures: short, medium, and long term. J. Health Care Finance 26(3) (2000 Spring) 5672.
10. GUPTE, M. D., VALLISHAYEE, R. S., ANANTHARAMAN, D. S., NAGARAJU, B., SREEVATSA, BALASUBRAMANYAM S., DE BRITTO, R. L., ELANGO, N., UTHAYAKUMARAN, N., MAHALINGAM, V. N., LOURDUSAMY, G., RAMALINGAM, A., KANNAN, S., and AROKIASAMY, J. Comparative leprosy vaccine trial in south India. Indian J. Lepr. 70 (1998) 369388.
11. INSTITUTO NACIONAL DE ESTADÍSTICA. Anuario estadístico Ministerio de Economía. Madrid. Años 1950 hasta 2001.
12. INSTITUTO NACIONAL DE ESTADÍSTICA. Censo de la población española. Ministerio de Economía. Madrid. Años 1950 hasta 2000.
13. INSTITUTO NACIONAL DE ESTADÍSTICA. Movimiento natural de la población española Ministerio de Economía. Madrid. Años 1950 hasta 2001.
14. IRGENS, L. M., MELO-CAEIRO, F., and LECHAT, M. F. Leprosy in Portugal 194680: epidemiologic patterns observed during declining incidence rates. Lepr. Rev. 61(1) (1990) 3249.
15. ITO, I. The epidemiological situation in south east Asia. Lepr. Rev. 52 (1981) 4351.
16. JAIN, S., REDDY, R. G., OSMANI, S. N., LOCKWOOD, D. N., and SUNEETHA, S. Childhood leprosy in an urban clinic, Hyderabad, India: clinical presentation and the role of household contacts. Lepr. Rev. 73 (2002) 248253.
17. KLEINBAUM, D. G., KUPPER, LL., and MORGENSTERN, H. Epidemiologic research: Principles and quantitative methods. New York: Van Nostrand Reinhold. 1982.
18. KRISHNAN, P. Mortality decline in India. 1951 1961: developmental versus public health program hypothesis. Soc. Sc. Med. 9 (1975) 475479.
19. LECHAT, M. F. The source of infection: an unsolved issue. Indian J. Lepr. 72 (2000) 169173.
20. LOCKWOOD, D. N., and SUNEETHA, S. Leprosy: too complex a disease for a simple elimination paradigm. Bull. World Health Organ. 83 (2005) 230235.
21. LOPEZ-VELEZ, R., SAEZ-VAQUERO, T., BLANCOAREVALO, J. L., and GOMEZ-MAMPASO, E. Leprosy simulating other diseases. Rev. Clin. Esp. 199 (1999) 369372.
22. MCKEOWN, T., and BROWN, R. G. Medical evidence related to English populations change in the eighteenth century. Pop. Stud. 9 (1965) 119123.
23. MEIMA, A., GUPTE, M. D., VAN OORTMARSSEN, G. J., and HABBEMA, J. D. F. Trends in leprosy case detection rates. Int. J. Lepr. Other Mycobact. Dis. 65 (1997) 305319.
24. MEIMA, A., IRGENS, L. M., VAN OORTMARSSEN, G. J., RICHARDUS, J. H., and HABBEMA, J. D. Disappearance of leprosy from Norway: an exploration of critical factors using an epidemiological modelling approach. International Journal of Epidemiology 31 (2002) 9911000.
25. MEIMA, A., RICHARDUS, J. H., and HABBEMA, J. D. Trends in leprosy case detection worldwide since 1985. Lepr. Rev. 75 (2004) 1933.
26. NOORDEEN, S. K. Descriptive epidemiology of leprosy. In: Leprosy. Hasting eds. Edimburgh. Churchill Livingstone Publications, 1985.
27. NOORDEEN, S. K. Toward the elimination of leprosy, the challenges and opportunities. Int. J. Lepr. Other Mycobact. Dis. 66 (1998) 218221.
28. PONNIGHAUS, J. M., FINE, P. E., STERNE, J. A., MALEMA, S. S., BLISS, L., and WILSON, R. J. Extended schooling and good housing conditions are associated with reduced risk of leprosy in rural Malawi. Int. J. Lepr. Other Mycobact. Dis. 62 (1994 Sep) 345352.
29. RAO, V. P., RAO, I. R., and PALANDE, D. D. Socioeconomic rehabilitation programmes of LEPRA India-methodology, results and application of needs-based socio-economic evaluation. Lepr, Rev. 71 (2000) 466471.
30. ROTHMAN, K. J., and GREENLAND, S. Modern Epidemiology. Philadelphia. 2nd ed. Lippincott Williams & Wilkins, 1998.
31. SAIKAWA, K. The effect of rapid socio-economic development on the frequency of leprosy in a population. Lepr. Rev. 52 (1981) 167175.
32. SANSARRICQ, H. Epidémiologie de la lèpre. In: La lèpre. Chapitre 5. Paris. Edit Ellipses Universités francophones, 1995, pp. 5472.
33. SHAH, F., and ABBEY, H. Effects of some factors on neonatal and postneonatal mortality. Milbank Mem. Fund Q 49 (1971) 3357.
34. SKINSNES, O. K. Effect of malnutrition on leprosy. Int. J. Lepr. Other Mycobact. Dis. 44 (1976) 374375.
35. TERENCIO DE LAS AGUAS, J. La lepra, pasado, presente y futuro. Valencia. Ed Generalitat Valenciana, 1999, pp. 5775.
36. VALVERDE, C. R., CANFIELD, D., TARARA, R., ESTEVES, M. I., and GORMUS, B. J. Spontaneous leprosy in a wild-caught cynomolgus macaque. Int. J. Lepr. Other Mycobact. Dis. 66 (1998) 140148.
37. WORTH, R. M. The disappearance of leprosy in a semi-isolated population (Nihau Island, Hawai). Int. J. Lepr. Other Mycobact. Dis. 31 (1963) 3435.
38. ZODPEY, S. P., BANSOD, B. S., SHRIKHANDE, S. N., MALDHURE, B. R., and KULKARNI, S. W. Protective effect of Bacillus Calmette Guérin (BCG) against leprosy: a population-based case-control study in Nagpur, India. Lepr. Rev. 70 (1999) 287294.
1. Servicio de Medicina Preventiva, Consorcio Hospital General Universitario de Valencia.
2. Centro Salud Pública de Manises, Consellería de Sanidad de la Comunidad Valenciana.
3. Servicio de Dermatología, Consorcio Hospital General Universitario de Valencia.
4. Advisor to the International Leprosy Association (ILA).
Name and address to whom requests for reprints should be addressed:
José L. Alfonso MD, Ph.D.
Servicio de Medicina Preventiva, Consorcio Hospital General Universitario de Valencia, Departamento de Medicina Preventiva y Salud Pública, Universidad de Valencia
Blasco Ibáñez 15
46010 - Valencia (Spain)
e-mail: jose.l.alfonso@uv.es