Temporal evolution and spatial distribution of leprosy in a municipality with low endemicity in São Paulo state, Brazil

Antônio Carlos Vieira Ramos José Francisco Martoreli Júnior Thaís Zamboni Berra Yan Mathias Alves Tatiana Pestana Barbosa Alessandro Rolim Scholze Ivaneliza Simionato de Assis Pedro Fredemir Palha Dulce Gomes Ricardo Alexandre Arcêncio About the authors

Abstract

Objective

To analyze the spatial and temporal distribution of leprosy in a scenario of low endemicity in the state of São Paulo, Brazil.

Methods

Ecological study with leprosy cases in Ribeirão Preto, between 2006 to 2016. The temporal trend of leprosy detection was verified through the decomposition of time series and identified areas of high and low occurrence of the disease using the Getis-Ord Gi* technique.

Results

There were 890 cases, and the detection rate showed an increasing trend in the period from 2011 to 2015, with an average growth of 1% per month. Areas of high occurrence of the disease were identified in the northern region of the city (99% and 95% confidence).

Conclusion

The temporal analysis showed that the rate of detection of leprosy presented an increasing trend, and the spatial analysis showed that the region of the municipality with the highest occurrence of the disease is characterized by presenting the greatest social inequalities.

Keywords:
Leprosy; Epidemiology; Spatial Analysis; Time Series Studies; Ecological Studies

Study contributions

Main results

The leprosy detection rate showed a rising trend between 2006 and 2016 in Ribeirão Preto, SP, whereby the region of the municipality identified as having the greatest occurrence of the disease was also characterized as a region with marked social inequality.

Implications for services

The results can contribute to planning and implementing strategic health surveillance actions to control leprosy in the municipality of Ribeirão Preto, SP, especially in the region where the disease occurs most.

Perspectives

Epidemiological investigations of leprosy with the aim of identifying areas or territories with greater occurrence of and vulnerability to the disease should be encouraged, in particular scenarios in which leprosy endemicity is considered to be low.

Introduction

According to World Health Organization (WHO) data, 202,185 new leprosy cases were detected worldwide in 2019, and of these some 80% were reported in India, Brazil and Indonesia.11. World Health Organization. Weekly epidemiological record: global leprosy (hansen disease) update, 2019: time to step-up prevention initiatives. Genebra: World Health Organization; 2020 [update 2020 sep 3; cited 2021 set 7]. Available from: https://www.who.int/publications/i/item/who-wer9536
https://www.who.int/publications/i/item/...
In the same period, Brazil recorded 27,863 new cases, being the country with the second highest number of cases globally.11. World Health Organization. Weekly epidemiological record: global leprosy (hansen disease) update, 2019: time to step-up prevention initiatives. Genebra: World Health Organization; 2020 [update 2020 sep 3; cited 2021 set 7]. Available from: https://www.who.int/publications/i/item/who-wer9536
https://www.who.int/publications/i/item/...

Despite the high levels of leprosy endemicity in Brazil, there have been significant improvements in the control of the disease in recent decades, especially with effect from the year 2000, following the putting in place of the international policies on leprosy elimination, which caused a reduction in the new case detection rate, from 26.2 per 100,000 inhabitants in 2001, to 11.2 cases per 100,000 inhabitants in 2019.22. Ministério da Saúde (BR). Departamento de Informática do SUS. Acompanhamento dos dados de hanseníase - Brasil [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 set 7]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswbr.def
http://tabnet.datasus.gov.br/cgi/tabcgi....
,33. Ribeiro MDA, Silva JCA, Oliveira SB. Estudo epidemiológico da hanseníase no Brasil: reflexão sobre as metas de eliminação. Rev Panam Salud Publica. 2018;42:1-7. doi: 10.26633/RPSP.2018.42
https://doi.org/10.26633/RPSP.2018.42...

Although there has been a reduction in the occurrence of the disease nationwide, it is noteworthy that its distribution is heterogeneous over the country, with differences being found between macro-regions, states and municipalities.44. Freitas LRS, Duarte EC, Garcia LP. Leprosy in Brazil and its association with characteristics of municipalities: ecological study, 2009-2011. Trop Med Int Health. 2014;19(10):1216-25. doi: 10.1111/tmi.12362
https://doi.org/10.1111/tmi.12362...
In 2019, the highest leprosy detection rates were found in the Midwest (51.8 per 100,000 inhab.), North (38.8 per 100,000 inhab.) and Northeast (26.6 per 100,000 inhab.) regions, while the lowest rates were reported by the Southeast (5.6 per 100,000 inhab.) and Southern (3.8 per 100,000 inhab.) regions of the country.22. Ministério da Saúde (BR). Departamento de Informática do SUS. Acompanhamento dos dados de hanseníase - Brasil [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 set 7]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswbr.def
http://tabnet.datasus.gov.br/cgi/tabcgi....

Despite their having the lowest detection rates, the data for the Southern and Southeast regions should be interpreted with caution, since low case detection does not mean evidence of low transmission.55. Oliveira KS, Souza J, Campos RB, Zilly A, Silva-Sobrinho RA. Avaliação dos indicadores epidemiológicos e operacionais para a hanseníase em municípios prioritários no estado do Paraná, 2001 a 2010. Epidemiol Serv Saude. 2015;24(3):507-16. doi: 10.5123/S1679-49742015000300016
https://doi.org/10.5123/S1679-4974201500...
Studies conducted in municipalities with low or medium leprosy endemicity have described substantial increases in detection rates after intensive surveillance and active searching for cases, revealing that low endemicity sites can have ‘hidden endemicity’ and active Mycobacterium leprae transmission hotspots.66. Bernardes Filho F, Paula NA, Leite MN, Abi-Rached TLC, Vernal S, Silva MB, et al. Evidence of hidden leprosy in a supposedly low endemic area of Brazil. Mem Inst Oswaldo Cruz. 2017;112(12):822-8. doi: 10.1590/0074-02760170173
https://doi.org/10.1590/0074-02760170173...
,77. Frade MAC, Paula NA, Gomes CM, Vernal S, Bernardes Filho F, Lugão HB, et al. Unexpectedly high leprosy seroprevalence detected using a random surveillance strategy in midwestern Brazil: a comparison of ELISA and a rapid diagnostic test. PLoS Negl Trop Dis. 2017;11(2):e0005375. doi: 10.1371/journal.pntd.0005375
https://doi.org/10.1371/journal.pntd.000...

In endemic countries, investigations focused on smaller units, such as municipalities, can provide a more realistic representation of the leprosy situation, reflecting the distribution of the disease more accurately.88. Blok DJ, De Vlas SJ, Richardus JH. Global elimination of leprosy by 2020: are we on track?. Parasit Vectors. 2015;8:548. doi: 10.1186/s13071-015-1143-4
https://doi.org/10.1186/s13071-015-1143-...
Another relevant aspect is that most scientific investigations of leprosy are carried out in regions where its endemicity is high, so that regions with lower detection rates, where endemicity is stable or falling, are little studied. Areas with low detection rates need to be investigated with the same attention as hyperendemic areas, as they can become neglected spaces that contribute to the maintenance of the transmission chain.99. Ribeiro GC, Barreto JG, Bueno IC, Vasconcelos BF, Lana FCF. Prevalence and spatial distribution of mycobacterium leprae infection in a medium endemicity municipality. Rev Rene. 2019;20:e39497. doi: 10.15253/2175-6783.20192039497
https://doi.org/10.15253/2175-6783.20192...

In view of these gaps in knowledge, the objective of this study was to analyze the spatial and temporal distribution of leprosy in a scenario of low endemicity in the state of São Paulo, Brazil.

Methods

This was an ecological study of new leprosy cases reported on the Notifiable Health Conditions Information System (SINAN) in Ribeirão Preto, SP, from 2006 to 2016, using the municipality’s 988 census tracts as units of ecological analysis.1010. Instituto Brasileiro de Geografia e Estatística. Base de informações do Censo Demográfico 2010: resultados do universo agregados por setor censitário - documentação do arquivo [Internet]. Brasília: Instituto Brasileiro de Geografia e Estatística; 2011 [citado 2021 set 7]. Disponível em: https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
https://www.ibge.gov.br/estatisticas/dow...

Ribeirão Preto has a population of 604,682 inhabitants.1111. Instituto Brasileiro de Geografia e Estatística. Cidades e estados - Ribeirão Preto [Internet].Brasília: Instituto Brasileiro de Geografia e Estatística; 2021 [citado 2021 set 7]. Disponível em: https://www.ibge.gov.br/cidades-e-estados/sp/ribeirao-preto.html
https://www.ibge.gov.br/cidades-e-estado...
Its public health system is divided into five Health Districts - Central, Eastern, Northern, Western and Southern -, comprising 49 Primary Health Care establishments, of which five are primary district health centers, 18 are Family Health Strategy centers and 26 are ordinary primary health care centers.1212. Secretaria Municipal de Saúde (Ribeirão Preto). Relação das unidades de saúde [Internet]. Ribeirão Preto: Secretaria Municipal de Saúde; 2021 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/portal/saude/relacao-unidades-saude
https://www.ribeiraopreto.sp.gov.br/port...
Figure 1 shows the location of Ribeirão Preto.

Figure 1
Map showing location of Ribeirão Preto, state of São Paulo, Brazil, 2016

In 2016, leprosy prevalence in Ribeirão Preto was 0.8 case per 10,000 inhab. while the detection rate was 7.7 new cases per 100,000 inhab.1313. Ministério da Saúde (BR). Departamento de Informática do SUS (DATASUS). Acompanhamento dos dados de hanseníase - São Paulo [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 set 7]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswsp.def
http://tabnet.datasus.gov.br/cgi/tabcgi....
According to parameters established by the Ministry of Health for the detection of new cases, the municipality presents low endemicity, considering the prevalence values and average strength of morbidity.1414. Ministério da Saúde (BR). Secretaria de Vigilância em Saúde. Departamento de Vigilância das Doenças Transmissíveis. Diretrizes para vigilância, atenção e eliminação da hanseníase como problema de saúde pública [Internet]. Brasília: Ministério da Saúde; 2016 [citado 2021 set 7]. Disponível em: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/diretrizes_para_._eliminacao_hanseniase_-_manual_-_3fev16_isbn_nucom_final_2.pdf
http://portal.saude.pe.gov.br/sites/port...

The data used in the study were retrieved from the SINAN database and included clinical and sociodemographic information about reported cases, such as operational classification (paucibacillary; multibacillary), grade of physical disability (GD) at diagnosis (G0D; G1D; G2D;), date of case notification, date of birth, sex (male; female), education (no education; elementary school; high school; higher education) and home address. Regarding GD, its classification varies according to the neural impairment of the affected individual: G0D corresponds to cases without any type of functional disability; G1D corresponds to those with loss of protective sensation; and G2D corresponds to individuals who, in addition to loss of protective sensation, present complications that include trophic ulcers, claw deformity, bone resorption in hands and/or feet and various eye lesions.1515. Amaral EP, Lana FCF. Análise espacial da hanseníase na microrregião de Almenara, MG, Brasil. Rev Bras Enferm. 2008;61(Esp):701-7. doi: 10.1590/S0034-71672008000700008
https://doi.org/10.1590/S0034-7167200800...

Access to the SINAN database was provided by the Epidemiological Survey Division of the Ribeirão Preto City Health Department.

A descriptive analysis of the clinical and sociodemographic variables of the leprosy cases was performed using RStudio software version 3.5.2.

With regard to the temporal analysis, the time series was built according to month of notification. The detection rates were calculated taking the total number of leprosy cases notified in the month as the numerator, while the denominator was the population of the municipality according to the 2010 Demographic Census and the estimated populations for each intercensal year (between the periods 2006-2009 and 2011-2016), multiplied by 100,000 inhabitants.1616. Ministério da Saúde (BR). Departamento de Informática do SUS. População Residente - Estimativas para o TCU - São Paulo [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 dez 10]. Disponível em: http://tabnet.datasus.gov.br/cgi/deftohtm.exe?ibge/cnv/poptsp.def
http://tabnet.datasus.gov.br/cgi/deftoht...

The detection rate time series was then decomposed using the time series decomposition method originally called Seasonal-Trend Decomposition using LOESS (STL).1717. Cleveland RB, Cleveland WS, McRae JE, Terpenning I. STL: a seasonal-trend decomposition procedure based on Loess [Internet]. J Off Stat. 1990 [cited 2021 set 7];6(1):3-33 Available from: https://www.wessa.net/download/stl.pdf
https://www.wessa.net/download/stl.pdf...
Assuming an additive form of decomposition, the leprosy detection rate in month t (Y t) is written according to the following formula:

Yt=St+Tt+Rt

where: St is the seasonal component; Tt is the trend component; and Rt is the residual or noise component.

Trend refers to the direction in which a time series develops over a given time interval, and may follow a rising, falling or stationary pattern.Seasonality is defined as identical patterns that a time series may follow that repeat regularly at fixed periods of time. “Noise” represents the fluctuations observed during the period covered by the series, generally irregular and random, perceptible only when the other components of the time series are removed.1818. Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. doi: 10.5123/S1679-49742015000300024
https://doi.org/10.5123/S1679-4974201500...

We chose to select only the trend component from all the components of the time series, in order to characterize the behavior of the leprosy detection rate over time. After decomposing the time series, we calculated Average Monthly Percentage Change (AMPC),1919. Bras AL, Gomes D, Filipe PA, Sousa B, Nunes C. Trends, seasonality and forecasts of pulmonary tuberculosis in Portugal. Int J Tuberc Lung Dis. 2014;18(10):1202-10. doi: 10.5588/ijtld.14.0158
https://doi.org/10.5588/ijtld.14.0158...
as the summary measure of the trend component over time, where n was the total number of months of study, expressed as follows:

AMPC=i=2nRatemonth iRatemonth i-1-1n-1×100

With regard to the temporal analyses, we used Rstudio version 3.5.2 to build the time series distribution graph and the detection rate trend.

In the spatial analysis, the leprosy cases were initially georeferenced based on the latitude and longitude information for their home addresses. Google Earth® was applied to identify latitude and longitude, and ArcGIS® 10.8 was used for georeferencing, building a file in shapefile format containing the location of each case. In this step, we excluded cases with no address or incomplete address information (street name or number unavailable), those who lived in the rural area of Ribeirão Preto and/or cases notified in Ribeirão Preto but who were not resident there.

After georeferencing the cases, the number of cases per unit of analysis was computed, enabling the calculation of the leprosy detection rate by census tract, by taking the number of cases per census tract as the numerator and the total population of each census tract as the denominator, and multiplying by 100,000 inhabitants.

We performed the analysis of the spatial clustering of leprosy cases using the Getis-Ord Gi* technique, based on the detection rates. This technique indicates local spatial association, taking the values for each census tract in the municipality based on a neighborhood matrix.2020. Getis A, Ord JK. The analysis of Spatial Association by use of distance statistics. Geogr Anal. 1992;24(3):189-206. doi: 10.1111/j.1538-4632.1992.tb00261.x
https://doi.org/10.1111/j.1538-4632.1992...
Getis-Ord Gi* is obtained using the following formula:

Gi*=j=1nwijxjj=1nxj

where: Gi* is the spatial autocorrelation statistic of an event i over n events (total number of components); xj characterizes the magnitude of variable x in j events over all n; and w is the spatial weight between components i and j. Inference as to the significance of Gi* is based on a standardized distribution:

Z(Gi*)=j=inwi, jxj-Χ-j=1nwi, jSnj=1nwi,j2-(j=1nwi,j)2n-1

where: X- is the mean and S is the variance:

Χ-=j=1nxjn

S=j=1nxj2n-(Χ-)2

Interpretation of this statistic is performed based on the Z-score and significance level values (α). A positive Z value with statistical evidence indicates a spatial cluster with greater occurrence of the event (hot spot), while a negative Z value with a statistical evidence indicates a spatial cluster with lesser occurrence of the event (cold spot). We used 90%, 95% and 99% confidence levels. We tested spatial association of leprosy in four different time periods (2006-2009; 2010-2013; 2014-2016; and the entire study period, 2006-2016), with the aim of verifying hot spot and cold spot distribution over time. We used ArcGIS® 10.8 to analyze spatial clustering and to prepare maps.

The study was submitted to the Research Ethics Committee of the Universidade de São Paulo Nursing School in Ribeirão Preto and was approved as per Certificate of Submission for Ethical Appraisal No. 44637215.0.0000.5393, issued on March 23rd 2018.

Results

A total of 890 leprosy cases were reported from 2006 to 2016, with a predominance of multibacillary cases (732; 82.2%) and presence of G1D at the time of diagnosis (347; 39.0%). There was a higher frequency of cases in the 30-59 age group (494; 55.5%), followed by those aged 60 or older (234; 26.3%), as well as among males (542; 60.9%) and those with elementary school education (476; 53.5%). Table 1 shows the descriptive statistics of the cases.

Table 1
Characteristics of leprosy cases in Ribeirão Preto, state of São Paulo, Brazil, 2006-2016

Figure 2 shows the leprosy detection rate time series and trend. From 2006 to 2010 the trend was stationary, followed by a rising trend as of 2011, which lasted until 2015, when the rate began to show a slight falling trend until the end of the period studied. Despite the falling trend after 2015, the detection rate values were higher at the end of the time series than at the beginning. The AMPC value indicated an average increase of 1.0% per month for the series.

Figure 2
Time series and temporal trend of leprosy detection rates, Ribeirão Preto, state of São Paulo, Brazil, 2006-2016

Of the total number of notified cases, 206 (23.1%) were excluded because their home address was outside Ribeirão Preto (they were notified in Ribeirão Preto but did not live there) and 29 (3.3%) because they lived in the rural area of the city. This left a total of 655 cases (73.6%) to be standardized for georeferencing. Of this total, 47 (7.2%) had a blank and/or incomplete address, so that the final number of cases to be georeferenced was 608 cases (92.8%) (Figure 3).

Figure 3
Leprosy cases georeferenced according to census tracts, Ribeirão Preto, state of São Paulo, Brazil, 2006-2016

Figure 4 shows the results of the leprosy detection rate local spatial association analysis using the Getis-Ord Gi* technique, which allowed the identification of areas of high case occurrence (hot spots) and low case occurrence (cold spots) and case distribution over time.

Figure 4
Clusters of high and low leprosy detection rate values, distributed according to census tracts, Ribeirão Preto, state of São Paulo, Brazil, 2006-2016

In the period from 2006 to 2016 (Figure 4A), the clusters with higher occurrence of leprosy were found in the north of Ribeirão Preto and included census tracts located in the Northern and Central Health Districts and a small part of the Eastern Health District. The clusters with lower occurrence were identified in the southern census tracts of the city, located in part of the area of the Central, Southern and Eastern Health Districts. Although some hot spots and cold spots are found in the same Health Districts, the areas of higher and lower occurrence are well delimited and geographically opposite to each other.

In the analysis of the spatial-temporal distribution of leprosy, the periods from 2006 to 2009 (Figure 4B) and 2010 to 2013 (Figure 4C) presented similar distribution, with clusters of higher occurrences of leprosy in the north of the city; Northern and Central Health Districts and a small part of the Eastern Health District.

In the period from 2014 to 2016 (Figure 4D), similarly to the other periods, clusters with high occurrence were found in the census tracts located in the Central and Northern Health Districts. However, unlike the other years, hot spots were also identified in census tracts of the Western Health District.

Discussion

Of the total leprosy cases reported in Ribeirão Preto, SP, the highest percentages corresponded to multibacillary cases, with a higher proportion of G1D at the time of diagnosis. Despite the predominance of G1D, it is noteworthy that 16.4% of total cases had G2D, and together G1D and G2D accounted for about 50% of notified cases. The predominance of multibacillary cases together with GD suggests the occurrence of active transmission of the disease, and the high proportion of G2D among new cases may be an indication of failures in case access to diagnosis and follow-up in Ribeirão Preto.2121. Santos AR, Ignotti E. Prevenção de incapacidade física por hanseníase no Brasil: análise histórica. Cien Saude Colet. 2020;25(10):3731-44. doi: 10.1590/1413-812320202510.30262018
https://doi.org/10.1590/1413-81232020251...

The 30-59 age group was the most affected by leprosy. This is an age range that includes a large portion of the economically active population, among whom the disease can affect job performance, result in limitations and time off work due to the disabilities caused by leprosy and, as a consequence, lead to workers facing sick leave, early retirement and decreased quality of life.2222. Monteiro LD, Alencar CHM, Barbosa JC, Braga KP, Castro MD, Heukelbach J. Incapacidades físicas em pessoas acometidas pela hanseníase no período pós-alta da poliquimioterapia em um municípiono norte do Brasil. Cad Saude Publica. 2013;29(5):909-20. doi: 10.1590/S0102-311X2013000500009
https://doi.org/10.1590/S0102-311X201300...

Besides the 30-59 age group, high proportions of cases were found in those aged 60 and over, this being the second most affected age group. The elderly constitute an important group for the control of leprosy, since the disabilities resulting from the disease are more frequent in this age group, which together with the aging process and other comorbidities, contribute to greater vulnerability and loss of autonomy in this group age.2323. Oliveira JSS, Reis ALM, Margalho LP, Lopes GL, Silva AR, Moraes NS, et al. Leprosy in elderly people and the profile of a retrospective cohort in an endemic region of the Brazilian Amazon. PLoS Negl Trop Dis. 2019;13(9):e0007709. doi: 10.1371/journal.pntd.0007709
https://doi.org/10.1371/journal.pntd.000...

Surveillance actions aimed at early diagnosis and prevention of disabilities should be targeted at all ages and reinforced in the most affected groups, in order to reduce transmission of the disease and prevent disabilities in the population as a whole.

In Ribeirão Preto, most of the cases reported were male, a result that is in keeping with the literature.2424. Marciano LHSC, Belone AFF, Rosa PS, Coelho NMB, Ghidella CC, Nardi SMT, et al. Epidemiological and geographical characterization of leprosy in a Brazilian hyperendemic municipality. Cad Saude Publica. 2018;34(8):e00197216. doi: 10.1590/0102-311X00197216
https://doi.org/10.1590/0102-311X0019721...
Occurrence of leprosy is more frequent in males, who are also more vulnerable when compared to females.2525. Nery JS, Ramond A, Pescarini JM, Alves A, Strina A, Ichihara MY, et al. Socioeconomic determinants of leprosy new case detection in the 100 Million Brazilian cohort: a population-based linkage study. Lancet Glob Heal. 2019;7(9):e1226-36. doi: 10.1016/S2214-109X(19)30260-8
https://doi.org/10.1016/S2214-109X(19)30...
The higher occurrence of the disease in men may be related to a greater exposure to M. leprae due to individual vulnerabilities, such as low demand and/or late demand for health services when compared to women, due to difficulties in accessing health services, which increase the risk of developing disabilities.2626. Nobre ML, Illarramendi X, Dupnik KM, Hacker MA, Nery JAC, Jerônimo SMB, et al. Multibacillary leprosy by population groups in Brazil: lessons from an observational study. PLoS Negl Trop Dis. 2017;11(2): e0005364. doi: 10.1371/journal.pntd.0005364
https://doi.org/10.1371/journal.pntd.000...
The finding that males were more affected by leprosy in Ribeirão Preto contributes to the planning of health practices that minimize the risk of leprosy transmission in this population, such as developing public campaigns.

Regarding education levels, most of the notified cases had only elementary schooling, with an average of less than ten years of study. The association between education and leprosy has been studied and has revealed evidence of the following association: low levels or absence of education are related to increased transmission of the disease, and may contribute to an up to twofold increase in its incidence and, consequently, to an increase in physical disabilities.2525. Nery JS, Ramond A, Pescarini JM, Alves A, Strina A, Ichihara MY, et al. Socioeconomic determinants of leprosy new case detection in the 100 Million Brazilian cohort: a population-based linkage study. Lancet Glob Heal. 2019;7(9):e1226-36. doi: 10.1016/S2214-109X(19)30260-8
https://doi.org/10.1016/S2214-109X(19)30...
,2727. Leano HAM, Araújo KMFA, Bueno IC, Niitsuma ENA, Lana FCF. Socioeconomic factors related to leprosy: an integrative literature review. Rev Bras Enferm. 2019;72(5):1405-15. doi: 10.1590/0034-7167-2017-0651
https://doi.org/10.1590/0034-7167-2017-0...
Populations with little or no schooling also have greater difficulty in accessing health services and less understanding of the meaning of health actions.2828. Grantz KH, Chabaari W, Samuel RK, Gershom B, Blum L, Worden L, et al. Spatial distribution of leprosy in India: an ecological study. Infect Dis Poverty. 2018;7(1):20. doi: 10.1186/s40249-018-0402-y
https://doi.org/10.1186/s40249-018-0402-...

The time series analysis of leprosy identified a rising trend in the detection rate from 2006 to 2016. Although a slight falling trend was seen at the end of the study period, it is noteworthy that the detection rate values were higher at the end of the series than at the beginning. During the study period, leprosy diagnosis and treatment actions were mainly carried out in municipal Reference Centers for leprosy care, so they were not fully decentralized to Primary Health Care. Despite this characteristic, professional training and active searching for cases are regularly carried out in the municipality, especially with effect from 2011 and 2012, coinciding with the sharp increase in the detection rate trend curve.2929. Neves LAS, Frade MAC, Silva ECS, Yano TK, Ferraz AN, Coelho WS, et al. Hanseníase: parcerias institucionais para promover a busca ativa e a qualificação da assistência. In: Anais do II Conferência Internacional em Epidemiologia [Internet]. 2012 nov 12-14; São Paulo, Brasil. São Paulo: Secretaria de Estado da Saúde; 2012 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/ssaude/pdf/hanseniase_parcerias_institucionais_promover.pdf
https://www.ribeiraopreto.sp.gov.br/ssau...

While the present study did not directly evaluate the impact of decentralized care and active case-finding actions, following professional training and community outreach measures, studies carried out in low-endemicity leprosy settings identified a significant increase in the number of reported cases, with a clear impact on the increase in detection rates, evidencing that such low-endemicity settings may harbor “hidden leprosy endemicity”.66. Bernardes Filho F, Paula NA, Leite MN, Abi-Rached TLC, Vernal S, Silva MB, et al. Evidence of hidden leprosy in a supposedly low endemic area of Brazil. Mem Inst Oswaldo Cruz. 2017;112(12):822-8. doi: 10.1590/0074-02760170173
https://doi.org/10.1590/0074-02760170173...
,77. Frade MAC, Paula NA, Gomes CM, Vernal S, Bernardes Filho F, Lugão HB, et al. Unexpectedly high leprosy seroprevalence detected using a random surveillance strategy in midwestern Brazil: a comparison of ELISA and a rapid diagnostic test. PLoS Negl Trop Dis. 2017;11(2):e0005375. doi: 10.1371/journal.pntd.0005375
https://doi.org/10.1371/journal.pntd.000...

The analysis of the time series indicated that the leprosy detection rate showed a rising trend over the study period as a whole, providing evidence that, even in a scenario of low endemicity, in which occurrence of the disease presents stability, there may be active transmission of M. leprae. This finding reinforces the importance of the need for constant actions involving active searching for and detection of new cases.

By using the Getis-Ord Gi* technique it was possible to identify areas of high and low occurrence of leprosy in Ribeirão Preto, over the period of time observed, showing that the disease is heterogeneously distributed in the city. In all analyzed periods, the areas of higher occurrence were located in the northern region of Ribeirão Preto, while the areas of low risk were located in the southern region of the city (identified here only on the map of the complete period, from 2006 to 2016). There was an exception in the final period of investigation, from 2014 to 2016, when hot spots were present in the Western Health District. Some Health Districts shared both hot spots and cold spots (Central and Eastern), although these are in opposite areas of these regions.

The Northern Health District, which had the highest occurrence of clusters of cases in all periods, is composed of the neighborhoods with the lowest social indicators among the five Health Districts, with the highest percentage of people with no income or income below one minimum wage, the lowest municipal human development index (HDI-M), the lowest school attendance rates and the highest number of crowded substandard living conditions (slums) in the municipality.3030. Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
https://www.ribeiraopreto.sp.gov.br/file...

The Central Health District, composed of both high-risk (to the north) and low-risk (to the south) clusters, is made up of the city’s oldest neighborhoods, with a high percentage of residents over 60 years old.3030. Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
https://www.ribeiraopreto.sp.gov.br/file...
The Central Health District has the main referral center for leprosy cases, although it has low Family Health Strategy coverage.3030. Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
https://www.ribeiraopreto.sp.gov.br/file...
Although the Central Health District has both hot spots and cold spots, they are geographically opposite to each other, so that its hot spots have similar characteristics to those of the Northern Health District, while its cold spots are similar in their characteristics to those of the Eastern and Southern Health Districts.

The Western Health District, which had hot spots in the period 2014-2016, is characterized by having a complex health service network, with a greater number of Primary Health Care centers. Regarding social indicators, in this district there is a predominance of households with families whose income is between one and five minimum wages and crowded substandard living conditions.3030. Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
https://www.ribeiraopreto.sp.gov.br/file...
The hot spots found in Western Health District in the final years of the study may reflect increased case detection in this region, thus influencing spatial distribution of the disease.

The high-risk clusters in the Eastern Health District comprise few census tracts and border the Central Health District tracts, so they may share similar characteristics with this region.

In the southern region of the city, besides the Central Health District cold spots, cold spots are also found in a small number of census tracts in the Eastern and Southern Health Districts; the latter Health District did not have hot spots, and is therefore the region in which occurrence of leprosy is lowest. The Eastern and Southern Health Districts are characterized by better socioeconomic indicators, better HDI-M, a population with higher levels of education and a low proportion of people with no income. Another characteristic of these Health Districts is that they have the highest concentration of neighborhoods where real estate speculation is more active.3030. Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
https://www.ribeiraopreto.sp.gov.br/file...

In Ribeirão Preto, analysis of the spatial association of the occurrence of cases and their distribution over time reinforces the fact that leprosy is a social disease: the areas with the highest concentration are located predominantly in the northern region of the city of Ribeirão Preto, historically characterized by great social inequality. These results can support the formulation of intervention strategies in these territories, especially actions involving active searching for cases, which should impact indicators and contribute to reducing transmission of the disease.

The main limitation of this study refers to its ecological design, which does not allow its results to be generalized to the individual level. Another limitation of the study is the use of secondary data, which are sometimes unavailable or incomplete. Finally, the time series rates were calculated based on population estimates, which may not fully reflect the size of the population in the period.

In conclusion, the time series analysis in Ribeirão Preto identified a rising trend in the leprosy detection rate during the study period, raising the hypothesis that, even in a scenario of low endemicity, existence of high transmission of the disease remains. The spatial-temporal analysis made clear that the disease is heterogeneous in Ribeirão Preto, with well-defined behavior over time, and that the north of the city is the region with the greatest social inequalities and areas with the highest rates of the disease.

The results of this investigation can inform health surveillance actions aimed at leprosy control. The understanding of aspects related to the profile, time trend and spatial distribution of cases in the municipality of Ribeirão Preto, São Paulo, can serve as a guide to where its main leprosy control actions should be focused.

References

  • 1
    World Health Organization. Weekly epidemiological record: global leprosy (hansen disease) update, 2019: time to step-up prevention initiatives. Genebra: World Health Organization; 2020 [update 2020 sep 3; cited 2021 set 7]. Available from: https://www.who.int/publications/i/item/who-wer9536
    » https://www.who.int/publications/i/item/who-wer9536
  • 2
    Ministério da Saúde (BR). Departamento de Informática do SUS. Acompanhamento dos dados de hanseníase - Brasil [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 set 7]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswbr.def
    » http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswbr.def
  • 3
    Ribeiro MDA, Silva JCA, Oliveira SB. Estudo epidemiológico da hanseníase no Brasil: reflexão sobre as metas de eliminação. Rev Panam Salud Publica. 2018;42:1-7. doi: 10.26633/RPSP.2018.42
    » https://doi.org/10.26633/RPSP.2018.42
  • 4
    Freitas LRS, Duarte EC, Garcia LP. Leprosy in Brazil and its association with characteristics of municipalities: ecological study, 2009-2011. Trop Med Int Health. 2014;19(10):1216-25. doi: 10.1111/tmi.12362
    » https://doi.org/10.1111/tmi.12362
  • 5
    Oliveira KS, Souza J, Campos RB, Zilly A, Silva-Sobrinho RA. Avaliação dos indicadores epidemiológicos e operacionais para a hanseníase em municípios prioritários no estado do Paraná, 2001 a 2010. Epidemiol Serv Saude. 2015;24(3):507-16. doi: 10.5123/S1679-49742015000300016
    » https://doi.org/10.5123/S1679-49742015000300016
  • 6
    Bernardes Filho F, Paula NA, Leite MN, Abi-Rached TLC, Vernal S, Silva MB, et al. Evidence of hidden leprosy in a supposedly low endemic area of Brazil. Mem Inst Oswaldo Cruz. 2017;112(12):822-8. doi: 10.1590/0074-02760170173
    » https://doi.org/10.1590/0074-02760170173
  • 7
    Frade MAC, Paula NA, Gomes CM, Vernal S, Bernardes Filho F, Lugão HB, et al. Unexpectedly high leprosy seroprevalence detected using a random surveillance strategy in midwestern Brazil: a comparison of ELISA and a rapid diagnostic test. PLoS Negl Trop Dis. 2017;11(2):e0005375. doi: 10.1371/journal.pntd.0005375
    » https://doi.org/10.1371/journal.pntd.0005375
  • 8
    Blok DJ, De Vlas SJ, Richardus JH. Global elimination of leprosy by 2020: are we on track?. Parasit Vectors. 2015;8:548. doi: 10.1186/s13071-015-1143-4
    » https://doi.org/10.1186/s13071-015-1143-4
  • 9
    Ribeiro GC, Barreto JG, Bueno IC, Vasconcelos BF, Lana FCF. Prevalence and spatial distribution of mycobacterium leprae infection in a medium endemicity municipality. Rev Rene. 2019;20:e39497. doi: 10.15253/2175-6783.20192039497
    » https://doi.org/10.15253/2175-6783.20192039497
  • 10
    Instituto Brasileiro de Geografia e Estatística. Base de informações do Censo Demográfico 2010: resultados do universo agregados por setor censitário - documentação do arquivo [Internet]. Brasília: Instituto Brasileiro de Geografia e Estatística; 2011 [citado 2021 set 7]. Disponível em: https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
    » https://www.ibge.gov.br/estatisticas/downloads-estatisticas.html
  • 11
    Instituto Brasileiro de Geografia e Estatística. Cidades e estados - Ribeirão Preto [Internet].Brasília: Instituto Brasileiro de Geografia e Estatística; 2021 [citado 2021 set 7]. Disponível em: https://www.ibge.gov.br/cidades-e-estados/sp/ribeirao-preto.html
    » https://www.ibge.gov.br/cidades-e-estados/sp/ribeirao-preto.html
  • 12
    Secretaria Municipal de Saúde (Ribeirão Preto). Relação das unidades de saúde [Internet]. Ribeirão Preto: Secretaria Municipal de Saúde; 2021 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/portal/saude/relacao-unidades-saude
    » https://www.ribeiraopreto.sp.gov.br/portal/saude/relacao-unidades-saude
  • 13
    Ministério da Saúde (BR). Departamento de Informática do SUS (DATASUS). Acompanhamento dos dados de hanseníase - São Paulo [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 set 7]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswsp.def
    » http://tabnet.datasus.gov.br/cgi/tabcgi.exe?sinannet/cnv/hanswsp.def
  • 14
    Ministério da Saúde (BR). Secretaria de Vigilância em Saúde. Departamento de Vigilância das Doenças Transmissíveis. Diretrizes para vigilância, atenção e eliminação da hanseníase como problema de saúde pública [Internet]. Brasília: Ministério da Saúde; 2016 [citado 2021 set 7]. Disponível em: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/diretrizes_para_._eliminacao_hanseniase_-_manual_-_3fev16_isbn_nucom_final_2.pdf
    » http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/diretrizes_para_._eliminacao_hanseniase_-_manual_-_3fev16_isbn_nucom_final_2.pdf
  • 15
    Amaral EP, Lana FCF. Análise espacial da hanseníase na microrregião de Almenara, MG, Brasil. Rev Bras Enferm. 2008;61(Esp):701-7. doi: 10.1590/S0034-71672008000700008
    » https://doi.org/10.1590/S0034-71672008000700008
  • 16
    Ministério da Saúde (BR). Departamento de Informática do SUS. População Residente - Estimativas para o TCU - São Paulo [Internet]. Brasília: Ministério da Saúde; 2021 [citado 2021 dez 10]. Disponível em: http://tabnet.datasus.gov.br/cgi/deftohtm.exe?ibge/cnv/poptsp.def
    » http://tabnet.datasus.gov.br/cgi/deftohtm.exe?ibge/cnv/poptsp.def
  • 17
    Cleveland RB, Cleveland WS, McRae JE, Terpenning I. STL: a seasonal-trend decomposition procedure based on Loess [Internet]. J Off Stat. 1990 [cited 2021 set 7];6(1):3-33 Available from: https://www.wessa.net/download/stl.pdf
    » https://www.wessa.net/download/stl.pdf
  • 18
    Antunes JLF, Cardoso MRA. Uso da análise de séries temporais em estudos epidemiológicos. Epidemiol Serv Saude. 2015;24(3):565-76. doi: 10.5123/S1679-49742015000300024
    » https://doi.org/10.5123/S1679-49742015000300024
  • 19
    Bras AL, Gomes D, Filipe PA, Sousa B, Nunes C. Trends, seasonality and forecasts of pulmonary tuberculosis in Portugal. Int J Tuberc Lung Dis. 2014;18(10):1202-10. doi: 10.5588/ijtld.14.0158
    » https://doi.org/10.5588/ijtld.14.0158
  • 20
    Getis A, Ord JK. The analysis of Spatial Association by use of distance statistics. Geogr Anal. 1992;24(3):189-206. doi: 10.1111/j.1538-4632.1992.tb00261.x
    » https://doi.org/10.1111/j.1538-4632.1992.tb00261.x
  • 21
    Santos AR, Ignotti E. Prevenção de incapacidade física por hanseníase no Brasil: análise histórica. Cien Saude Colet. 2020;25(10):3731-44. doi: 10.1590/1413-812320202510.30262018
    » https://doi.org/10.1590/1413-812320202510.30262018
  • 22
    Monteiro LD, Alencar CHM, Barbosa JC, Braga KP, Castro MD, Heukelbach J. Incapacidades físicas em pessoas acometidas pela hanseníase no período pós-alta da poliquimioterapia em um municípiono norte do Brasil. Cad Saude Publica. 2013;29(5):909-20. doi: 10.1590/S0102-311X2013000500009
    » https://doi.org/10.1590/S0102-311X2013000500009
  • 23
    Oliveira JSS, Reis ALM, Margalho LP, Lopes GL, Silva AR, Moraes NS, et al. Leprosy in elderly people and the profile of a retrospective cohort in an endemic region of the Brazilian Amazon. PLoS Negl Trop Dis. 2019;13(9):e0007709. doi: 10.1371/journal.pntd.0007709
    » https://doi.org/10.1371/journal.pntd.0007709
  • 24
    Marciano LHSC, Belone AFF, Rosa PS, Coelho NMB, Ghidella CC, Nardi SMT, et al. Epidemiological and geographical characterization of leprosy in a Brazilian hyperendemic municipality. Cad Saude Publica. 2018;34(8):e00197216. doi: 10.1590/0102-311X00197216
    » https://doi.org/10.1590/0102-311X00197216
  • 25
    Nery JS, Ramond A, Pescarini JM, Alves A, Strina A, Ichihara MY, et al. Socioeconomic determinants of leprosy new case detection in the 100 Million Brazilian cohort: a population-based linkage study. Lancet Glob Heal. 2019;7(9):e1226-36. doi: 10.1016/S2214-109X(19)30260-8
    » https://doi.org/10.1016/S2214-109X(19)30260-8
  • 26
    Nobre ML, Illarramendi X, Dupnik KM, Hacker MA, Nery JAC, Jerônimo SMB, et al. Multibacillary leprosy by population groups in Brazil: lessons from an observational study. PLoS Negl Trop Dis. 2017;11(2): e0005364. doi: 10.1371/journal.pntd.0005364
    » https://doi.org/10.1371/journal.pntd.0005364
  • 27
    Leano HAM, Araújo KMFA, Bueno IC, Niitsuma ENA, Lana FCF. Socioeconomic factors related to leprosy: an integrative literature review. Rev Bras Enferm. 2019;72(5):1405-15. doi: 10.1590/0034-7167-2017-0651
    » https://doi.org/10.1590/0034-7167-2017-0651
  • 28
    Grantz KH, Chabaari W, Samuel RK, Gershom B, Blum L, Worden L, et al. Spatial distribution of leprosy in India: an ecological study. Infect Dis Poverty. 2018;7(1):20. doi: 10.1186/s40249-018-0402-y
    » https://doi.org/10.1186/s40249-018-0402-y
  • 29
    Neves LAS, Frade MAC, Silva ECS, Yano TK, Ferraz AN, Coelho WS, et al. Hanseníase: parcerias institucionais para promover a busca ativa e a qualificação da assistência. In: Anais do II Conferência Internacional em Epidemiologia [Internet]. 2012 nov 12-14; São Paulo, Brasil. São Paulo: Secretaria de Estado da Saúde; 2012 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/ssaude/pdf/hanseniase_parcerias_institucionais_promover.pdf
    » https://www.ribeiraopreto.sp.gov.br/ssaude/pdf/hanseniase_parcerias_institucionais_promover.pdf
  • 30
    Secretaria Municipal de Saúde (Ribeirão Preto). Departamento de Vigilância em Saúde e Planejamento. Fatores de risco relacionados à saúde da população residente na zona urbana de Ribeirão Preto (SP) 2008-2011 [Internet]. Ribeirão Preto: Secretaria Municipal da Saúde; 2011 [citado 2021 set 7]. Disponível em: https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf
    » https://www.ribeiraopreto.sp.gov.br/files/ssaude/pdf/fatores-risco.pdf

  • Associated academic work

    Work derived from the Ph.D. project of author Antônio Carlos Vieira Ramos, entitled ‘Markers of social inequality and their relationship to leprosy: a geo-epidemiological approach’, begun in 2018. The author is enrolled at the ‘Nursing in Public Health’ Postgraduate Program at the Universidade de São Paulo Nursing School in Ribeirão Preto.
  • Financial

    This work received support from the Coordination for the Improvement of Higher Education Personnel/Ministry of Education - Funding Code 001.

Publication Dates

  • Publication in this collection
    20 Apr 2022
  • Date of issue
    2022

History

  • Received
    04 Oct 2021
  • Accepted
    05 Jan 2022
Secretaria de Vigilância em Saúde - Ministério da Saúde do Brasil Brasília - Distrito Federal - Brazil
E-mail: leilapgarcia@gmail.com