Spatial analysis of Chikungunya fever incidence and the associated socioeconomic, demographic, and vector infestation factors in municipalities of Pernambuco, Brazil, 2015–2021

Maísa Aguiar-Santos Liana Gabriele da Cruz Mendes Diogenes Ferreira dos Passos Tamyris Gomes da Silva Santos Rosanny Holanda Freitas Benevides Lins Ana Cristina Pedrosa do Monte About the authors

ABSTRACT

Objective:

To identify the spatial patterns of chikungunya fever (CHIKF) and the associated socioeconomic, demographic, and vector infestation factors in the 1st Health Region of Pernambuco (1st HRP).

Methods:

This ecological study used a spatial analysis of Mean Incidence Rates (MIR) of probable cases of CHIKF reported among residents of the 19 municipalities of the 1st HRP, in 2015–2021. The univariate and bivariate global Moran indexes (I) were estimated. From the significant associations (p<0.05), clusters were identified using the local Moran index and maps.

Results:

A predominance of the largest CHIKF rates was identified in the east. However, there was a heterogeneous distribution of rates across municipalities, which may have contributed to the absence of spatial autocorrelation of CHIKF (I=0.03; p=0.294) in univariate I. The bivariate I revealed a positive spatial correlation between CHIKF and the Municipal Human Development Index (MHDI) (I=0.245; p=0.038), but with a cluster of cities with low incidences and low MHDI in the west. There was no spatial correlation between CHIKF and the other variables analyzed: population density, Gini index, social vulnerability index, and building infestation index for Aedes aegypti.

Conclusions:

The results suggest that only the MHDI influenced the occurrence of CHIKF in the 1st HRP, so that municipalities in the west demonstrated spatial dependence between lower values of MHDI and MIR. However, this spatial correlation may have occurred due to possible underreporting in the area. These findings can assist in the (re)orientation of resources for surveillance and health care services.

Keywords:
Arbovirus infections; Aedes aegypti; Chikungunya virus; Spatial analysis; Socioeconomic factors; Public health surveillance

INTRODUCTION

Arboviruses transmitted by the Aedes aegypti mosquito are one of the greatest obstacles to Public Health worldwide, including Chikungunya fever (CHIKF), which reached 137,025 reported cases in the Americas in 202111 Pan American Health Organization. Epidemiological update for Dengue, Chikungunya and Zika in 2021 [Internet]. Washington: PAHO; 2021 [cited on Jul. 1, 2022]. Available at: https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/annual-arbovirus-bulletin-2021.html
https://www3.paho.org/data/index.php/en/...
. Of these records, approximately 97% occurred in Brazil11 Pan American Health Organization. Epidemiological update for Dengue, Chikungunya and Zika in 2021 [Internet]. Washington: PAHO; 2021 [cited on Jul. 1, 2022]. Available at: https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/annual-arbovirus-bulletin-2021.html
https://www3.paho.org/data/index.php/en/...
, where the Pernambuco state stood out for the highest incidence rate of CHIKF (329 probable cases/100,000 inhabitants)22 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos pelo mosquito Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 51, 2021 [Internet]. Boletim Epidemiologico 2021; 52(48): 1-6 [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/boletins-epidemiologicos/edicoes/2021/boletim-epidemiologico-vol-52-no-48.pdf/view
https://www.gov.br/saude/pt-br/centrais-...
.

Autochthonous cases of CHIKF in Pernambuco were reported for the first time in 201533 Pernambuco. Secretaria Estadual de Saúde. Circulação do Vírus de Chikungunya (CHIKV) no estado de Pernambuco - orientações aos serviços de saúde. Nota Técnica 38/15 [Internet]. 2015 [cited on Feb. 3, 2022]. Available at: https://www.cievspe.com/_files/ugd/3293a8_fe1aec3817e44caa89a4eb0c851f1f43.pdf
https://www.cievspe.com/_files/ugd/3293a...
. After that, the state experienced alternating between endemic and epidemic periods, with considerable increases in the years 2016, 2020, and 202122 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos pelo mosquito Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 51, 2021 [Internet]. Boletim Epidemiologico 2021; 52(48): 1-6 [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/boletins-epidemiologicos/edicoes/2021/boletim-epidemiologico-vol-52-no-48.pdf/view
https://www.gov.br/saude/pt-br/centrais-...
,44 Silva NM, Teixeira RAG, Cardoso CG, Siqueira Junior JB, Coelho GE, Oliveira ESF. Vigilância de chikungunya no Brasil: desafios no contexto da Saúde Pública. Epidemiol Serv Saúde 2018; 27(3): e2017127. https://doi.org/10.5123/S1679-49742018000300003
https://doi.org/10.5123/S1679-4974201800...
,55 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos por Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 53, 2020 [Internet]. Boletim Epidemiologico 2021; 52(3): 1-15. [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/epidemiologicos/edicoes/2021/boletim_epidemiologico_svs_3.pdf/view
https://www.gov.br/saude/pt-br/centrais-...
. Of the 12 health regions into which the state is divided66 Pernambuco. I Gerência Regional de Saúde. Secretaria de Saúde. Mapa de Saúde da I Região de Saúde/I Gerência Regional de Saúde [Internet]. Pernambuco: Secretaria de Saúde; 2021 [cited on Feb. 1, 2022]. Available at: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/mapa_de_saude_2020_i_regiao_de_saude.pdf
http://portal.saude.pe.gov.br/sites/port...
, the 1st Health Region of Pernambuco (1st HRP) was responsible for 81.2% of all probable cases notified in the previous year (2021)22 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos pelo mosquito Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 51, 2021 [Internet]. Boletim Epidemiologico 2021; 52(48): 1-6 [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/boletins-epidemiologicos/edicoes/2021/boletim-epidemiologico-vol-52-no-48.pdf/view
https://www.gov.br/saude/pt-br/centrais-...
, which contributed to registering the highest incidence rate (600.5/100,000 residents) among regions.

Characterized by intense arthralgia and possible atypical manifestations (such as neurological manifestations) that can persist for years77 Brasil. Ministério da Saúde. Vigilância em Saúde no Brasil 2003|2019: da criação da Secretaria de Vigilância em Saúde aos dias atuais [Internet]. Boletim Epidemiologico. 2019;50(n. esp.):1-154. [cited on Feb. 5, 2022]. Available at: https://ameci.org.br/wp-content/uploads/2019/09/boletim-especial-21ago19-web.pdf
https://ameci.org.br/wp-content/uploads/...
, CHIKF has the potential to cause overload of health services and work absenteeism77 Brasil. Ministério da Saúde. Vigilância em Saúde no Brasil 2003|2019: da criação da Secretaria de Vigilância em Saúde aos dias atuais [Internet]. Boletim Epidemiologico. 2019;50(n. esp.):1-154. [cited on Feb. 5, 2022]. Available at: https://ameci.org.br/wp-content/uploads/2019/09/boletim-especial-21ago19-web.pdf
https://ameci.org.br/wp-content/uploads/...
. Both the magnitude and dissemination of this disease are influenced by the presence and abundance of the A. aegypti vector88 Costa SSB, Branco MRFC, Vasconcelos VV, Queiroz RCS, Araujo AS, Câmara APB, et al. Autoregressive spatial modeling of possible cases of dengue, chikungunya, and Zika in the capital of Northeastern Brazil. Rev Soc Bras Med Trop 2021; 54: e0223. https://doi.org/10.1590/0037-8682-0223-2021
https://doi.org/10.1590/0037-8682-0223-2...
, whose proliferation depends on socioeconomic, demographic, and environmental factors88 Costa SSB, Branco MRFC, Vasconcelos VV, Queiroz RCS, Araujo AS, Câmara APB, et al. Autoregressive spatial modeling of possible cases of dengue, chikungunya, and Zika in the capital of Northeastern Brazil. Rev Soc Bras Med Trop 2021; 54: e0223. https://doi.org/10.1590/0037-8682-0223-2021
https://doi.org/10.1590/0037-8682-0223-2...
1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
. Studies show that urban areas with precarious sanitation conditions and high population density provide ideal breeding grounds for these vectors, favoring an increase in the incidence of arboviruses1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
,1111 Brasil. Ministério da Saúde. Lixo é o principal criadouro do mosquito da dengue nas regiões Norte, Centro-Oeste e Sul [Internet]. Brasília: Ministério da Saúde; 2014. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/noticia/lixo-e-o-principal-criadouro-do-mosquito-da-dengue-nas-regioes-norte-centro-oeste-e-sul#:~:text=Segundo%20dados%20do%20Levantamento%20de,50%2C1%25%20dos%20criadouros.
https://www.unasus.gov.br/noticia/lixo-e...
. Spaces like this can be common in Pernambuco, a state that also stands out for having the highest poverty rate1212 Instituto Mobilidade e Desenvolvimento Social. Atualização dos resultados de pobreza: 2021. [Internet] 2022. [cited on Jul. 1, 2022]. Available at: https://imdsbrasil.org/atualizacao-dos-resultados-de-pobreza-2021
https://imdsbrasil.org/atualizacao-dos-r...
and for occupying 19th place in the ranking of the Human Development Index (HDI) in the country1313 Instituto Brasileiro de Geografia e Estatística. Pernambuco. Panorama [Internet]. 2021 [cited on Feb. 10, 2022]. Available at: https://cidades.ibge.gov.br/brasil/pe/panorama
https://cidades.ibge.gov.br/brasil/pe/pa...
. However, there are contradictory findings in the literature1414 Falcão Neto PAO, Branco MRFC, Costa SSB, Câmara APB, Marques TMNF, Araujo AS, et al. Análise espacial da taxa de detecção de casos suspeitos de síndrome congênita pelo vírus Zika, Maranhão, 2015 a 2018. Rev Bras Epidemiol 2022; 25: e220002. https://doi.org/10.1590/1980-549720220002
https://doi.org/10.1590/1980-54972022000...
1717 Teixeira TRA, Cruz OG. Spatial modeling of dengue and socio-environmental indicators in the city of Rio de Janeiro, Brazil. Cad Saúde Pública 2011; 27(3): 591-602. https://doi.org/10.1590/S0102-311X2011000300019
https://doi.org/10.1590/S0102-311X201100...
that indicate both the presence of associations between high incidence of arboviruses and better socioeconomic conditions and the absence of any association.

In view of the above, it is essential to know which determinants interfere in the health-disease process of a territory, so that there is guidance for effective control and prevention actions. That said, spatial analysis techniques become extremely useful instruments, as they have been used to investigate the spatial dynamics of diseases such as arboviruses in the social context in different locations, addressing possible associated factors and detecting vulnerable areas88 Costa SSB, Branco MRFC, Vasconcelos VV, Queiroz RCS, Araujo AS, Câmara APB, et al. Autoregressive spatial modeling of possible cases of dengue, chikungunya, and Zika in the capital of Northeastern Brazil. Rev Soc Bras Med Trop 2021; 54: e0223. https://doi.org/10.1590/0037-8682-0223-2021
https://doi.org/10.1590/0037-8682-0223-2...
,99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
,1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
1818 Lippi CA, Stewart-Ibarra AM, Romero M, Lowe R, Mahon R, Van Meerbeeck CJ, et al. Spatiotemporal tools for emerging and endemic disease hotspots in small areas: an analysis of dengue and chikungunya in Barbados, 2013–2016. Am J Trop Med Hyg 2020; 103(1): 149-56. https://doi.org/10.4269/ajtmh.19-0919
https://doi.org/10.4269/ajtmh.19-0919...
. However, there are gaps in the knowledge of these factors in the context of CHIKF in Pernambuco.

Therefore, the present study aimed to identify the spatial patterns of CHIKF and the associated socioeconomic, demographic, and vector infestation factors, in the municipalities of the 1st HRP, in the period from 2015 to 2021. Thus, valuable directions can be indicated for Surveillance in Health.

METHODS

This is an ecological study with spatial analysis referring to the incidence rates of probable cases of CHIKF in residents of the 1st HRP, from 2015 to 2021, as described by Costa et al.1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
. The units of spatial analysis used were municipalities.

Study area

Pernambuco is divided into 12 health regions, including the 1st Region, which concentrates most of the state's resident population (43%; 4,259,679) and has a population density of 1,444 inhabitants/km2.66 Pernambuco. I Gerência Regional de Saúde. Secretaria de Saúde. Mapa de Saúde da I Região de Saúde/I Gerência Regional de Saúde [Internet]. Pernambuco: Secretaria de Saúde; 2021 [cited on Feb. 1, 2022]. Available at: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/mapa_de_saude_2020_i_regiao_de_saude.pdf
http://portal.saude.pe.gov.br/sites/port...
With a nucleus in the capital Recife, the 1st Region is organized into three health micro-regions and has 19 municipalities, in addition to the State District of Fernando de Noronha (Figure 1)66 Pernambuco. I Gerência Regional de Saúde. Secretaria de Saúde. Mapa de Saúde da I Região de Saúde/I Gerência Regional de Saúde [Internet]. Pernambuco: Secretaria de Saúde; 2021 [cited on Feb. 1, 2022]. Available at: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/mapa_de_saude_2020_i_regiao_de_saude.pdf
http://portal.saude.pe.gov.br/sites/port...
.

Figure 1
Map referring to the division of municipalities in the 1st Health Region of Pernambuco, Brazil, into microregions.

Definition of probable cases and data collection

Probable cases of CHIKF correspond to the total number of reported cases, excluding discarded cases that had a negative laboratory diagnosis or that were diagnosed for other diseases, according to the Ministry of Health1919 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de dengue, febre de chikungunya e doença aguda pelo vírus Zika até a Semana Epidemiológica 34 de 2018. Boletim Epidemiologico 2018; 49(2): 1 [cited on Jan. 12, 2023]. Available at: https://repositorio.observatoriodocuidado.fiocruz.br/bitstream/handle/handle/1591/2018-043.pdf?sequence=1&isAllowed=y
https://repositorio.observatoriodocuidad...
. Data were obtained from the Notifiable Diseases Information System (Sistema de Informação de Agravos de Notificação – SINAN), considering the municipality of residence and epidemiological year of onset of symptoms, made available, in April 2022, by the Epidemiological Surveillance of the 1st Regional Health Management of the state of Pernambuco (Iᵃ Gerência Regional de Saúde do estado de Pernambuco – 1st GERES/SES-PE).

Socioeconomic and demographic variables

The socioeconomic and demographic variables included in the spatial analysis were: Gini Index and Municipal Human Development Index (MHDI), extracted from Atlas Brasil2020 Atlas do Desenvolvimento Humano no Brasil. Índice de Desenvolvimento Humano Municipal [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#idhm-all
http://www.atlasbrasil.org.br/perfil/mun...
,2121 Atlas do Desenvolvimento Humano no Brasil. Índice de Gini [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#sec-renda
http://www.atlasbrasil.org.br/perfil/mun...
; Social Vulnerability Index (SVI), obtained by the Institute of Applied Economic Research (Instituto de Pesquisa Econômica Aplicada – IPEA)2222 Instituto de Pesquisa Econômica Aplicada. Atlas da vulnerabilidade social nos municípios brasileiros [Internet]. Brasília: IPEA; 2015 [cited on May 2, 2022]. Available at: http://ivs.ipea.gov.br/images/publicacoes/Ivs/publicacao_atlas_ivs.pdf
http://ivs.ipea.gov.br/images/publicacoe...
; and population density (inhabitants/km2), provided by the Brazilian Institute of Geography and Statistics (Instituto Brasileiro de Geografia e Estatística – IBGE)2323 Instituto Brasileiro de Geografia e Estatística. Estimativas de população. População residente estimada [Internet] [cited on Mar. 1, 2022]. Available at: https://sidra.ibge.gov.br/tabela/6579
https://sidra.ibge.gov.br/tabela/6579...
. The data were from the last Census, carried out in 2010.

All socioeconomic indices used range from zero to one. The MHDI indicates the degree of human development in the municipality related to three subindices — longevity, education, and income —, classified as: very low (<0.500), low (0.500–0.599), medium (0.600–0.699), high (0.700–0.799), and very high (0.800–1)2020 Atlas do Desenvolvimento Humano no Brasil. Índice de Desenvolvimento Humano Municipal [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#idhm-all
http://www.atlasbrasil.org.br/perfil/mun...
. The Gini index highlights the degree of inequality in the per capita household income distribution; the closer to one, the greater the social inequality: low (<0.550), medium (0.550–0.699), high (0.700–0.799), and very high (0.800–1)2121 Atlas do Desenvolvimento Humano no Brasil. Índice de Gini [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#sec-renda
http://www.atlasbrasil.org.br/perfil/mun...
. SVI expresses the indication of exclusion based on the subindices urban infrastructure, human capital and income and work, being categorized as: very low (<0.200), low (0.201–0.300), medium (0.301–0.400), high (0.401– 0.500), and very high (0.501–1)2222 Instituto de Pesquisa Econômica Aplicada. Atlas da vulnerabilidade social nos municípios brasileiros [Internet]. Brasília: IPEA; 2015 [cited on May 2, 2022]. Available at: http://ivs.ipea.gov.br/images/publicacoes/Ivs/publicacao_atlas_ivs.pdf
http://ivs.ipea.gov.br/images/publicacoe...
.

Aedes aegypti infestation

The vector infestation of the A. aegypti mosquito was estimated from the Building Infestation Index (BII), which assesses the risk situation of arbovirus transmission by measuring the population level of the vector when calculating the percentage of properties with the presence of larvae of this Culicidae within the total number of surveyed properties2424 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância das Doenças Transmissíveis. Levantamento rápido de índices para Aedes Aegypti (LIRAa) para vigilância entomológica do Aedes aegypti no Brasil: metodologia para avaliação dos índices de Breteau e Predial e tipo de recipientes [Internet]. Brasília: Ministério da Saúde; 2013 [cited on Feb. 11, 2022]. Available at: https://bvsms.saude.gov.br/bvs/publicacoes/manual_liraa_2013.pdf
https://bvsms.saude.gov.br/bvs/publicaco...
. The calculation is performed using the Rapid Index Survey for A. aegypti (Levantamento Rápido de Índices para A. aegypti – LIRAa) and the BII is classified as: satisfactory (<1%), alert (1–3.9%), and at risk (>3.9%)2424 Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância das Doenças Transmissíveis. Levantamento rápido de índices para Aedes Aegypti (LIRAa) para vigilância entomológica do Aedes aegypti no Brasil: metodologia para avaliação dos índices de Breteau e Predial e tipo de recipientes [Internet]. Brasília: Ministério da Saúde; 2013 [cited on Feb. 11, 2022]. Available at: https://bvsms.saude.gov.br/bvs/publicacoes/manual_liraa_2013.pdf
https://bvsms.saude.gov.br/bvs/publicaco...
. These data were provided by the Environmental Surveillance of 1st GERES/SES-PE.

Calculation of incidence rates

The outcome of interest in this study was the mean CHIKF incidence rate per municipality, calculated using the ratio between the average of the total number of probable CHIKF cases in each municipality over the 7-year period (2015 to 2021) and the estimated population, multiplied by 100,000 inhabitants, according to Costa et al.1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
. The estimated population for 2018, equivalent to half of the period studied, was obtained from IBGE2323 Instituto Brasileiro de Geografia e Estatística. Estimativas de população. População residente estimada [Internet] [cited on Mar. 1, 2022]. Available at: https://sidra.ibge.gov.br/tabela/6579
https://sidra.ibge.gov.br/tabela/6579...
.

Spatial analysis

All municipalities were included in the spatial analysis, except for the Fernando de Noronha archipelago, due to different environmental conditions and the difficult visualization of spatial patterns on maps based on their distance. In order to establish the spatial autocorrelation linked only to the average municipal incidence rate of CHIKF, the univariate Global Moran Index was measured. The bivariate Global Moran Index was calculated in order to determine the spatial correlation between this rate and the possible associated factors (population density, Gini Index, HDI, SVI, and BII).

The Moran Global Index evaluates the pattern of spatial distribution and provides a general measure of the degree of spatial association present in the dataset, considering geographic adjacency2525 Druck S, Carvalho MS, Câmara G, Monteiro AMV. Análise espacial de dados geográficos. Brasília: Embrapa; 2004.. It can vary from -1 to +1: indices close to +1 mean that neighboring areas present similar behavior (positive spatial relationship) and indices close to -1 portray that municipalities with low values tend to be close to neighbors with high values and vice versa (negative spatial association). When the index is close to zero, the absence of spatial dependence is revealed due to the random spatial distribution of the data2525 Druck S, Carvalho MS, Câmara G, Monteiro AMV. Análise espacial de dados geográficos. Brasília: Embrapa; 2004..

From the significant global spatial associations, it was possible to evidence them locally through the graphical representation (LisaMap) of the Local Moran Index and its calculation referring to a specific value for each area2626 Anselin L. GeoDa: local spatial autocorrelation. LISA and Local Moran [Internet]. 2020 [cited on Apr. 10, 2022]. Available at: https://geodacenter.github.io/workbook/6a_local_auto/lab6a.html
https://geodacenter.github.io/workbook/6...
. In these representations, clusters of municipalities with significant spatial dependencies are observed. When a municipality has a positive local spatial association with neighbors, it means that both have similar scenarios (high-high or low-low relationship for the value(s) of the variable(s)), analyzing the clusters. In contrast, when municipalities have different conditions (high-low or low-high ratio), outliers are analyzed, which indicate points of negative spatial association2626 Anselin L. GeoDa: local spatial autocorrelation. LISA and Local Moran [Internet]. 2020 [cited on Apr. 10, 2022]. Available at: https://geodacenter.github.io/workbook/6a_local_auto/lab6a.html
https://geodacenter.github.io/workbook/6...
.

The indices were estimated using the GeoDa software, version 1.20.0.8. The “first-order” Queen neighborhood matrix was used, involving data from neighbors that touch at least one geographic side of each municipality; and 999 permutations in the pseudosignificance tests, considering a significance level of p<0.05 for spatial relationships. Subsequently, the maps were made using the Quantum Geographic Information (QGis) software, version 3.22.6.

RESULTS

Between 2015 and 2021, 53,269 probable cases (31,961 confirmed) of CHIKF were reported in residents of the municipalities that make up the Iᵃ RSP, for an average incidence rate of 180.8 cases/100,000 inhabitants over 7 years.

Table 1 shows the number of probable cases by year of onset of symptoms, according to demographic characterization, and it is possible to note that there was a greater involvement of female patients, aged 20 to 39 years and of brown race/color.

Table 1
Demographic description of probable cases of Chikungunya fever by year of symptom onset, 1st Health Region of Pernambuco (including Fernando de Noronha), Brazil, 2015–2021. Percentages include discarded/blank data.

As for the annual incidence rates, 2021 (600.5), 2016 (414.6), and 2020 (104.9) stood out, in contrast to the years 2015 (62.9), 2019 (35.3), 2017 (23.2), and 2018 (16.7), which showed the lowest values.

The distribution of municipal incidence rates showed that 10 of the 19 municipalities had high (92.57–168.38) and very high (>168.30) rates: Ilha de Itamaracá (369.92), Cabo de Santo Agostinho (277.41), Recife (262.93), Camaragibe (196.51), Olinda (171.14), Ipojuca (165.92), Moreno (142.44), Abreu e Lima (135.37), and Paulista (108.12). Thus, there was a predominance on the east coast, covering, in sequence, microregions I and III; and a municipality in microregion II, located in the west (Figures 1 and 2).

Figure 2
Map with municipalities categorized according to the parameters of the incidence rate of probable cases of Chikungunya fever, 1st Health Region of Pernambuco, Brazil, 2015–2021.

The Global Moran Index (I) identified a positive spatial autocorrelation for the CHIKF incidence rate, but not significant (I=0.03; p=0.294). However, in the bivariate analyses, there was a significant spatial correlation only between the MHDI and the incidence rate of CHIKF (I=0.245; p=0.038), making the MHDI the only spatially associated factor. The maps elaborated from the analysis of the Local Moran Index showed the nucleus of a “low-low” cluster in Pombos, represented by the dark blue color, portraying that the municipality had low incidence, being surrounded by neighbors that had low values of MHDI (Figure 3). It is noteworthy that population density, the Gini Index, SVI, and BII were non-significant variables and, therefore, factors not associated with CHIKF (Table 2).

Figure 3
Map of the bivariate local spatial analysis (LisaMap) between the incidence rate of probable cases of Chikungunya fever and the Municipal Human Development Index, with its spatial associations, 1st Health Region of Pernambuco, Brazil, 2015–2021.
Table 2
Univariate analysis of the incidence rate of probable cases of Chikungunya fever and bivariate analyses between this rate and demographic, socioeconomic, and environmental variables, 1st Health Region of Pernambuco, Brazil, 2015–2021.

DISCUSSION

This study is the first record to explore a spatial analysis between the incidence of CHIKF and demographic, socioeconomic, and vector infestation indicators in areas of the state of Pernambuco. The univariate analysis identified that the municipalities of 1st HRP did not show spatial autocorrelation for the incidence rates of probable CHIKF cases. In the bivariate analyses, the MHDI revealed a spatial correlation with the disease, showing that municipalities in the west showed spatial dependence between lower MHDI values and lower incidence rates of this arbovirus.

There was a similarity between the demographic characterization of dengue in Pernambuco (2015–2020)2727 Lima Filho CA, Lima AES, Arcanjo RMG, Silva DL, Jesus GF, Albuquerque AOBC, et al. Perfil epidemiológico dos casos de dengue no estado de Pernambuco, Brasil. Research, Society and Development 2022; 11(2): e36711225891. http://dx.doi.org/10.33448/rsd-v11i2.25891
http://dx.doi.org/10.33448/rsd-v11i2.258...
and that evidenced for CHIK in the 1st Region. In the entire study area, CHIKF was present, with a predominance of incidence rates considered high and very high in the Metropolitan Region of Recife (MRR) (microregions I and III); and low and medium in microregion II. However, the spatial distribution of these rates was heterogeneous in the 19 municipalities of the 1st Region, which may have determined the random spatial pattern observed by the Global Moran Index, as found in investigations carried out in the state of Maranhão1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
and in Barbados1818 Lippi CA, Stewart-Ibarra AM, Romero M, Lowe R, Mahon R, Van Meerbeeck CJ, et al. Spatiotemporal tools for emerging and endemic disease hotspots in small areas: an analysis of dengue and chikungunya in Barbados, 2013–2016. Am J Trop Med Hyg 2020; 103(1): 149-56. https://doi.org/10.4269/ajtmh.19-0919
https://doi.org/10.4269/ajtmh.19-0919...
, an island country in the Caribbean.

In contrast, studies conducted in the cities of Rio de Janeiro99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
and São Luís88 Costa SSB, Branco MRFC, Vasconcelos VV, Queiroz RCS, Araujo AS, Câmara APB, et al. Autoregressive spatial modeling of possible cases of dengue, chikungunya, and Zika in the capital of Northeastern Brazil. Rev Soc Bras Med Trop 2021; 54: e0223. https://doi.org/10.1590/0037-8682-0223-2021
https://doi.org/10.1590/0037-8682-0223-2...
demonstrated positive and significant global spatial autocorrelation for the incidence rates of suspected and probable cases of this arbovirus, respectively, indicating that adjacent neighborhoods have similar rates. Freitas et al.2828 Freitas JR, Nascimento GILA, Ferreira DSA, Santiago EJP, Moreira GR, Silva ASA, et al. Análise espaço-temporal da incidência de febre Chikungunya no estado de Pernambuco. Research, Society and Development 2020; 9(9): e288997114. http://dx.doi.org/10.33448/rsd-v9i9.7114
http://dx.doi.org/10.33448/rsd-v9i9.7114...
showed spatial dependence in municipalities with high incidences of CHIKF in the MRR, Pernambuco, in 2018. The authors, however, did not analyze the possible explanatory factors for the distribution of CHIKF in the state.

In the present spatial analysis, it was found that the incidence rate of CHIKF was spatially directly associated with the MHDI. Consistent with this result, Costa et al.1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
, when investigating CHIKF and Zika between 2015 and 2016, indicated a positive spatial association between the incidence of arboviruses transmitted by A. aegypti and the HDI in municipalities in Maranhão. These researchers found relationships between better MHDI values and high incidence rates, suggesting that municipalities with higher MHDI are generally the most populous and have greater access to health services, influencing the quantification of notifications.

Diverging, the current study noted spatial dependence in four municipalities with low incidence rates and low MHDI values in the west, more specifically in microregion II. However, these low incidences may be related to underreporting. With the lowest population densities and lowest levels of development in the entire region, micro-region II probably did not have the infrastructure that would allow greater access to health services1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
or technical-operational conditions to detect, notify, and investigate cases2929 Brasil. Ministério da Saúde. Características dos Indicadores. Fichas de Qualificação. Taxa de incidência de dengue – D.2.3 – 2012 [Internet]. [cited on Jun. 30, 2022]. Available at: http://fichas.ripsa.org.br/2012/d-2-3/?l=pt_BR
http://fichas.ripsa.org.br/2012/d-2-3/?l...
, thus being able to have cooperated for the lowest number of reported cases. However, the quality of these services was not measured.

The other socioeconomic variables analyzed, the SVI and the Gini Index (social inequality), did not present spatial associations with the incidence rate of CHIKF. Studies with conflicting results about the socioeconomic parameter and the triad of arboviruses (dengue, Zika, and CHIKF) have been recorded: while some indicate the absence of spatial dependence related to the Gini Index and the incidence of these diseases (for dengue1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
), others point to positive (for suspected cases of dengue15) or inverse correlations (for CHIKF1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
and Zika1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
).

The differences between the associations found in different studies may be related to the type of spatial unit used, indicating that secondary data used by the same type of grouping supposedly present similar results1717 Teixeira TRA, Cruz OG. Spatial modeling of dengue and socio-environmental indicators in the city of Rio de Janeiro, Brazil. Cad Saúde Pública 2011; 27(3): 591-602. https://doi.org/10.1590/S0102-311X2011000300019
https://doi.org/10.1590/S0102-311X201100...
. Another issue is that socioeconomic levels are difficult to measure due to the diverse and complex variables involved for aggregate data1717 Teixeira TRA, Cruz OG. Spatial modeling of dengue and socio-environmental indicators in the city of Rio de Janeiro, Brazil. Cad Saúde Pública 2011; 27(3): 591-602. https://doi.org/10.1590/S0102-311X2011000300019
https://doi.org/10.1590/S0102-311X201100...
. Therefore, deeper investigations are needed on the relationship between socioeconomic aspects and CHIKF and other arboviruses, as well as vector infestation.

With regard to demographic data, the spatial dependence analysis did not identify a correlation between population density and the incidence rate of CHIKF. Similar to this result, a survey carried out between 2000 and 2018, in Brazilian municipalities, identified that population density was not spatially associated with the incidence rates of dengue3030 Azevedo TS, Lorenz C, Chiaravalloti-Neto F. Spatiotemporal evolution of dengue outbreaks in Brazil. Trans R Soc Trop Med Hyg 2020; 114(8): 593-602. https://doi.org/10.1093/trstmh/traa030
https://doi.org/10.1093/trstmh/traa030...
— an arbovirus transmitted by the same vector as CHIKF. On the other hand, studies suggest that more populated environments favor the spatial diffusion of infectious diseases, especially those transmitted by the A. aegypti mosquito99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
,1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
,3131 Skalinski LM, Dalvi APR, Natividade M, Braga JU, Costa MCN, Santana EB, et al. The triple epidemics of arboviruses in Feira de Santana, Brazilian Northeast: epidemiological characteristics and diffusion patterns. Epidemics. 2022; 38: 100541. https://doi.org/10.1016/j.epidem.2022.100541
https://doi.org/10.1016/j.epidem.2022.10...
, due to the supply of a high number of individuals susceptible to new infections and the adaptation of the vector to perilous and intradomiciliary environments99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
.

Areas with a high number of inhabitants per square kilometer are usually the result of disorderly urban growth1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
. Together with this, the process of integration of territories with common interests from a “core city”, the so-called metropolization, has been causing several effects to Brazilian cities, such as socioeconomic advances that act, for example, in the MHDI; but also includes social inequality1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
.

The MRR was the most unequal metropolis in Brazil in the fourth quarter of 20203232 Salata AR, Ribeiro MG. Boletim Desigualdade nas Metrópoles [Internet] 2021. Porto Alegre: Observatório das Metrópoles. Pontifícia Universidade Católica do Rio Grande do Sul; 2021. [cited on Jul. 1, 2022]. Available at: https://www.pucrs.br/datasocial/wp-content/uploads/sites/300/2022/07/boletim_desigualdade_nas_metropoles_06.pdf
https://www.pucrs.br/datasocial/wp-conte...
and represented the majority (63%) of the municipalities in the 1st Region that did not have a Municipal Sanitation Plan or that had a plan under development in 20173333 Instituto Brasileiro de Geografia e Estatística. Base de dados – Munic 2017 – saneamento básico [Internet]. [cited on Jun. 28, 2022]. Available at: https://censoagro2017.ibge.gov.br/agencia-detalhe-de-midia.html?view=mediaibge&catid=2103&id=2341
https://censoagro2017.ibge.gov.br/agenci...
. This may have contributed for the high incidence rates of CHIKF found in the area, although there was no spatial association of this rate in the 1st Region.

Precarious sanitation conditions, such as inadequate water supply, sewage, and garbage collection, are also consequences of the intense disorderly urbanization1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
and are part of the components in the calculation of the SVI2222 Instituto de Pesquisa Econômica Aplicada. Atlas da vulnerabilidade social nos municípios brasileiros [Internet]. Brasília: IPEA; 2015 [cited on May 2, 2022]. Available at: http://ivs.ipea.gov.br/images/publicacoes/Ivs/publicacao_atlas_ivs.pdf
http://ivs.ipea.gov.br/images/publicacoe...
. Such characteristics create favorable scenarios for the proliferation of the A. aegypti vector, due to the opportunity of standing water reservoirs for the laying of its eggs77 Brasil. Ministério da Saúde. Vigilância em Saúde no Brasil 2003|2019: da criação da Secretaria de Vigilância em Saúde aos dias atuais [Internet]. Boletim Epidemiologico. 2019;50(n. esp.):1-154. [cited on Feb. 5, 2022]. Available at: https://ameci.org.br/wp-content/uploads/2019/09/boletim-especial-21ago19-web.pdf
https://ameci.org.br/wp-content/uploads/...
,1111 Brasil. Ministério da Saúde. Lixo é o principal criadouro do mosquito da dengue nas regiões Norte, Centro-Oeste e Sul [Internet]. Brasília: Ministério da Saúde; 2014. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/noticia/lixo-e-o-principal-criadouro-do-mosquito-da-dengue-nas-regioes-norte-centro-oeste-e-sul#:~:text=Segundo%20dados%20do%20Levantamento%20de,50%2C1%25%20dos%20criadouros.
https://www.unasus.gov.br/noticia/lixo-e...
. However, a study carried out in the capital of Rio de Janeiro99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
showed that the incidence rate of suspected cases of CHIKF showed a positive spatial correlation with variables contrary to these components, indicating that high incidences were spatially related to high percentages of households connected to the public water supply network and adequate sanitation.

With regard to A. aegypti infestation, the BII did not obtain a spatial correlation with the CHIKF incidence rate, which is consistent with the analysis conducted in Maranhão1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
during 2015 and 2016, despite being a contradictory finding. The LIRAa was not performed in some cycles of seven municipalities, mainly in 2020 and 2021, which influenced the measurement of the annual average of the BII and may have contributed to the failure to identify a correlation between this index and the disease incidence rate, in addition to the possibility of deficiency in the execution of the LIRAa also in other years. However, it is possible that the BII does not reflect the reality of the risk of spreading arboviruses in the municipalities, given that the index does not measure the density of A. aegypti vector adults, who are responsible for transmitting the virus3434 Ribeiro MS, Ferreira DF, Azevedo RC, Santos GBG, Medronho RA. Índices larvais de Aedes aegypti e incidência de dengue: um estudo ecológico no Estado do Rio de Janeiro, Brasil. Cad Saúde Pública 2021; 37(7): e00263320. https://doi.org/10.1590/0102-311X00263320
https://doi.org/10.1590/0102-311X0026332...
. New studies are needed to investigate this association and, perhaps, new methodological analyses to predict the risk of transmission of arboviruses, with appreciation and integration of epidemiological and entomological surveillance.

When evaluating notifications of arboviruses in Brazil during the COVID-19 pandemic, Lisboa et al.3535 Lisboa TR, Serafim IBM, Serafim JCM, Ramos AC, Nascimento RM, Roner MNB. Relação entre incidência de casos de arboviroses e a pandemia da COVID-19. Rev Interdiscip Ciênc Apl 2022; 6(10): 31-6. https://doi.org/10.18226/25253824.v6.n10.04
https://doi.org/10.18226/25253824.v6.n10...
found that, from 2019 to 2020 (first pandemic year), there was a reduction in CHIKF cases. However, in the 1st Region, it was possible to notice that CHIKF cases have increased continuously since 2019, especially in the pandemic: compared to the past years, 2020 and 2021 showed, respectively, substantial increases of approximately 298 and 576% in the records. This increase may have been due to the consequences of measures to prevent and control COVID-193636 Brasil. Ministério da Saúde. Nota Informativa n° 8/2020. Recomendações aos Agentes de Combate a Endemias (ACE) para adequação das ações de vigilância e controle de zoonoses frente à atual situação epidemiológica referente ao Coronavírus (COVID-19) [Internet]. Brasília: Ministério da Saúde; 2020. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/especial/covid19/pdf/106
https://www.unasus.gov.br/especial/covid...
, such as home isolation and the reduction in routine home visits by Endemic Disease Control Agents, which are essential for the control of arboviruses3535 Lisboa TR, Serafim IBM, Serafim JCM, Ramos AC, Nascimento RM, Roner MNB. Relação entre incidência de casos de arboviroses e a pandemia da COVID-19. Rev Interdiscip Ciênc Apl 2022; 6(10): 31-6. https://doi.org/10.18226/25253824.v6.n10.04
https://doi.org/10.18226/25253824.v6.n10...
. Such actions made it possible, for example, to increase peri- and indoor exposure to A. aegypti and to increase waste production, which, added to inadequate disposal, may have favored the proliferation of the CHIKF vector1111 Brasil. Ministério da Saúde. Lixo é o principal criadouro do mosquito da dengue nas regiões Norte, Centro-Oeste e Sul [Internet]. Brasília: Ministério da Saúde; 2014. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/noticia/lixo-e-o-principal-criadouro-do-mosquito-da-dengue-nas-regioes-norte-centro-oeste-e-sul#:~:text=Segundo%20dados%20do%20Levantamento%20de,50%2C1%25%20dos%20criadouros.
https://www.unasus.gov.br/noticia/lixo-e...
. It should be noted that, even so, it is possible that there is underreporting due, also, to similar clinical manifestations in the early stages of COVID-193737 Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med 2020; 382(18): 1708-20. https://doi.org/10.1056/NEJMoa2002032
https://doi.org/10.1056/NEJMoa2002032...
and arboviruses3838 Rodriguez-Morales AJ, Villamil-Gómez WE, Franco-Paredes C. The arboviral burden of disease caused by co-circulation and co-infection of dengue, chikungunya and Zika in the Americas. Travel Med Infect Dis 2016; 14(3): 177-9. https://doi.org/10.1016/j.tmaid.2016.05.004
https://doi.org/10.1016/j.tmaid.2016.05....
, meaning that the increase in cases verified may be even greater.

In addition to the possibility of an ecological fallacy, which can occur in ecological studies, this research had some limitations. The use of aggregated secondary data to describe large regions such as municipalities may not have captured local variations, thus suggesting the use of smaller spatial analysis units, such as neighborhoods or census tracts99 Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
https://doi.org/10.1017/S095026882100180...
,1717 Teixeira TRA, Cruz OG. Spatial modeling of dengue and socio-environmental indicators in the city of Rio de Janeiro, Brazil. Cad Saúde Pública 2011; 27(3): 591-602. https://doi.org/10.1590/S0102-311X2011000300019
https://doi.org/10.1590/S0102-311X201100...
. There may also have been inconsistencies or delays in completing, quantifying, and processing information in the CHIKF notification forms and in the subsequent entry of data1616 Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
https://doi.org/10.1590/S1678-9946201860...
. Another limiting factor is related to the use of data from the 2010 Census for socioeconomic indicators, given that they are the most recent available. The outdated data, together with an alleged underreporting of the disease, may not be showing the real situation of the municipalities.

In view of the above, it is important to train health professionals about notifications and active and continuous investigation of CHIKF cases. Furthermore, the investigation of other socioeconomic, demographic, and environmental variables, as well as climatic aspects1010 Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
https://doi.org/10.1590/1413-81232020251...
, can show the reality of each location and better explain the spatial dynamics of CHIKF, taking as examples: productivity of vector breeding sites, temperature, precipitation, and percentage coverage of the Family Health Strategy. In addition, studies that analyze the simultaneous circulation of arboviruses are necessary in regions of Pernambuco.

Despite the limitations, the study had the following advantages: knowledge of the profile of CHIKF in the 1st HRP since the first autochthonous case registered in the state, showing that high incidence rates were concentrated in the MRR (east), while low rates were concentrated in microregion II (west), where there was possibly underreporting; and the identification of aspects that may be interfering with the incidence rates of CHIKF in the spatial scope of the largest health region in Pernambuco.

Although the only possible explanatory factor for the distribution pattern of CHIKF cases was the MHDI, covering low incidences, it should be noted that socioeconomic, demographic, and environmental variables are important to explain part of the phenomenon and to serve as subsidies to (re)formulation of effective public policy planning and the strengthening of the Unified Health System (Sistema Único de Saúde – SUS). The development of public policies aimed at improving indicators would make urban spaces less unequal, better planned and healthier, as well as promoting a better quality of life for residents. In this way, it is expected that there will be a contribution to reducing the transmission of this arbovirus, as well as other diseases, improving collective health.

  • Funding: none.

REFERENCES

  • 1
    Pan American Health Organization. Epidemiological update for Dengue, Chikungunya and Zika in 2021 [Internet]. Washington: PAHO; 2021 [cited on Jul. 1, 2022]. Available at: https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/annual-arbovirus-bulletin-2021.html
    » https://www3.paho.org/data/index.php/en/mnu-topics/indicadores-dengue-en/annual-arbovirus-bulletin-2021.html
  • 2
    Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos pelo mosquito Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 51, 2021 [Internet]. Boletim Epidemiologico 2021; 52(48): 1-6 [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/boletins-epidemiologicos/edicoes/2021/boletim-epidemiologico-vol-52-no-48.pdf/view
    » https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/boletins-epidemiologicos/edicoes/2021/boletim-epidemiologico-vol-52-no-48.pdf/view
  • 3
    Pernambuco. Secretaria Estadual de Saúde. Circulação do Vírus de Chikungunya (CHIKV) no estado de Pernambuco - orientações aos serviços de saúde. Nota Técnica 38/15 [Internet]. 2015 [cited on Feb. 3, 2022]. Available at: https://www.cievspe.com/_files/ugd/3293a8_fe1aec3817e44caa89a4eb0c851f1f43.pdf
    » https://www.cievspe.com/_files/ugd/3293a8_fe1aec3817e44caa89a4eb0c851f1f43.pdf
  • 4
    Silva NM, Teixeira RAG, Cardoso CG, Siqueira Junior JB, Coelho GE, Oliveira ESF. Vigilância de chikungunya no Brasil: desafios no contexto da Saúde Pública. Epidemiol Serv Saúde 2018; 27(3): e2017127. https://doi.org/10.5123/S1679-49742018000300003
    » https://doi.org/10.5123/S1679-49742018000300003
  • 5
    Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de arboviroses urbanas causados por vírus transmitidos por Aedes (dengue, chikungunya e zika), semanas epidemiológicas 1 a 53, 2020 [Internet]. Boletim Epidemiologico 2021; 52(3): 1-15. [cited on Feb. 5, 2022]. Available at: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/epidemiologicos/edicoes/2021/boletim_epidemiologico_svs_3.pdf/view
    » https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/epidemiologicos/edicoes/2021/boletim_epidemiologico_svs_3.pdf/view
  • 6
    Pernambuco. I Gerência Regional de Saúde. Secretaria de Saúde. Mapa de Saúde da I Região de Saúde/I Gerência Regional de Saúde [Internet]. Pernambuco: Secretaria de Saúde; 2021 [cited on Feb. 1, 2022]. Available at: http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/mapa_de_saude_2020_i_regiao_de_saude.pdf
    » http://portal.saude.pe.gov.br/sites/portal.saude.pe.gov.br/files/mapa_de_saude_2020_i_regiao_de_saude.pdf
  • 7
    Brasil. Ministério da Saúde. Vigilância em Saúde no Brasil 2003|2019: da criação da Secretaria de Vigilância em Saúde aos dias atuais [Internet]. Boletim Epidemiologico. 2019;50(n. esp.):1-154. [cited on Feb. 5, 2022]. Available at: https://ameci.org.br/wp-content/uploads/2019/09/boletim-especial-21ago19-web.pdf
    » https://ameci.org.br/wp-content/uploads/2019/09/boletim-especial-21ago19-web.pdf
  • 8
    Costa SSB, Branco MRFC, Vasconcelos VV, Queiroz RCS, Araujo AS, Câmara APB, et al. Autoregressive spatial modeling of possible cases of dengue, chikungunya, and Zika in the capital of Northeastern Brazil. Rev Soc Bras Med Trop 2021; 54: e0223. https://doi.org/10.1590/0037-8682-0223-2021
    » https://doi.org/10.1590/0037-8682-0223-2021
  • 9
    Queiroz ERS, Medronho RA. Spatial analysis of the incidence of Dengue, Zika and Chikungunya and socioeconomic determinants in the city of Rio de Janeiro, Brazil. Epidemiol Infect 2021; 149: e188. https://doi.org/10.1017/S0950268821001801
    » https://doi.org/10.1017/S0950268821001801
  • 10
    Almeida LS, Cota AL, Rodrigues DF. Saneamento, arboviroses e determinantes ambientais: impactos na saúde urbana. Ciên Saúde Coletiva 2020; 25(10): 3857-68. https://doi.org/10.1590/1413-812320202510.30712018
    » https://doi.org/10.1590/1413-812320202510.30712018
  • 11
    Brasil. Ministério da Saúde. Lixo é o principal criadouro do mosquito da dengue nas regiões Norte, Centro-Oeste e Sul [Internet]. Brasília: Ministério da Saúde; 2014. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/noticia/lixo-e-o-principal-criadouro-do-mosquito-da-dengue-nas-regioes-norte-centro-oeste-e-sul#:~:text=Segundo%20dados%20do%20Levantamento%20de,50%2C1%25%20dos%20criadouros
    » https://www.unasus.gov.br/noticia/lixo-e-o-principal-criadouro-do-mosquito-da-dengue-nas-regioes-norte-centro-oeste-e-sul#:~:text=Segundo%20dados%20do%20Levantamento%20de,50%2C1%25%20dos%20criadouros
  • 12
    Instituto Mobilidade e Desenvolvimento Social. Atualização dos resultados de pobreza: 2021. [Internet] 2022. [cited on Jul. 1, 2022]. Available at: https://imdsbrasil.org/atualizacao-dos-resultados-de-pobreza-2021
    » https://imdsbrasil.org/atualizacao-dos-resultados-de-pobreza-2021
  • 13
    Instituto Brasileiro de Geografia e Estatística. Pernambuco. Panorama [Internet]. 2021 [cited on Feb. 10, 2022]. Available at: https://cidades.ibge.gov.br/brasil/pe/panorama
    » https://cidades.ibge.gov.br/brasil/pe/panorama
  • 14
    Falcão Neto PAO, Branco MRFC, Costa SSB, Câmara APB, Marques TMNF, Araujo AS, et al. Análise espacial da taxa de detecção de casos suspeitos de síndrome congênita pelo vírus Zika, Maranhão, 2015 a 2018. Rev Bras Epidemiol 2022; 25: e220002. https://doi.org/10.1590/1980-549720220002
    » https://doi.org/10.1590/1980-549720220002
  • 15
    Silva ETC, Olinda RA, Pachá AS, Costa AO, Brito AL, Pedraza DF. Análise espacial da distribuição dos casos de dengue e sua relação com fatores socioambientais no estado da Paraíba, Brasil, 2007-2016. Saúde Debate 2020; 44(125): 465-77. https://doi.org/10.1590/0103-1104202012514
    » https://doi.org/10.1590/0103-1104202012514
  • 16
    Costa SS, Branco MRFC, Aquino Junior J, Rodrigues ZMR, Queiroz RCS, Araujo AS, et al. Spatial analysis of probable cases of dengue fever, chikungunya fever and zika virus infections in Maranhão State, Brazil. Rev Inst Med Trop São Paulo 2018; 60: e62. https://doi.org/10.1590/S1678-9946201860062
    » https://doi.org/10.1590/S1678-9946201860062
  • 17
    Teixeira TRA, Cruz OG. Spatial modeling of dengue and socio-environmental indicators in the city of Rio de Janeiro, Brazil. Cad Saúde Pública 2011; 27(3): 591-602. https://doi.org/10.1590/S0102-311X2011000300019
    » https://doi.org/10.1590/S0102-311X2011000300019
  • 18
    Lippi CA, Stewart-Ibarra AM, Romero M, Lowe R, Mahon R, Van Meerbeeck CJ, et al. Spatiotemporal tools for emerging and endemic disease hotspots in small areas: an analysis of dengue and chikungunya in Barbados, 2013–2016. Am J Trop Med Hyg 2020; 103(1): 149-56. https://doi.org/10.4269/ajtmh.19-0919
    » https://doi.org/10.4269/ajtmh.19-0919
  • 19
    Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Monitoramento dos casos de dengue, febre de chikungunya e doença aguda pelo vírus Zika até a Semana Epidemiológica 34 de 2018. Boletim Epidemiologico 2018; 49(2): 1 [cited on Jan. 12, 2023]. Available at: https://repositorio.observatoriodocuidado.fiocruz.br/bitstream/handle/handle/1591/2018-043.pdf?sequence=1&isAllowed=y
    » https://repositorio.observatoriodocuidado.fiocruz.br/bitstream/handle/handle/1591/2018-043.pdf?sequence=1&isAllowed=y
  • 20
    Atlas do Desenvolvimento Humano no Brasil. Índice de Desenvolvimento Humano Municipal [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#idhm-all
    » http://www.atlasbrasil.org.br/perfil/municipio/261160#idhm-all
  • 21
    Atlas do Desenvolvimento Humano no Brasil. Índice de Gini [Internet]. [cited on May 2, 2022]. Available at: http://www.atlasbrasil.org.br/perfil/municipio/261160#sec-renda
    » http://www.atlasbrasil.org.br/perfil/municipio/261160#sec-renda
  • 22
    Instituto de Pesquisa Econômica Aplicada. Atlas da vulnerabilidade social nos municípios brasileiros [Internet]. Brasília: IPEA; 2015 [cited on May 2, 2022]. Available at: http://ivs.ipea.gov.br/images/publicacoes/Ivs/publicacao_atlas_ivs.pdf
    » http://ivs.ipea.gov.br/images/publicacoes/Ivs/publicacao_atlas_ivs.pdf
  • 23
    Instituto Brasileiro de Geografia e Estatística. Estimativas de população. População residente estimada [Internet] [cited on Mar. 1, 2022]. Available at: https://sidra.ibge.gov.br/tabela/6579
    » https://sidra.ibge.gov.br/tabela/6579
  • 24
    Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Departamento de Vigilância das Doenças Transmissíveis. Levantamento rápido de índices para Aedes Aegypti (LIRAa) para vigilância entomológica do Aedes aegypti no Brasil: metodologia para avaliação dos índices de Breteau e Predial e tipo de recipientes [Internet]. Brasília: Ministério da Saúde; 2013 [cited on Feb. 11, 2022]. Available at: https://bvsms.saude.gov.br/bvs/publicacoes/manual_liraa_2013.pdf
    » https://bvsms.saude.gov.br/bvs/publicacoes/manual_liraa_2013.pdf
  • 25
    Druck S, Carvalho MS, Câmara G, Monteiro AMV. Análise espacial de dados geográficos. Brasília: Embrapa; 2004.
  • 26
    Anselin L. GeoDa: local spatial autocorrelation. LISA and Local Moran [Internet]. 2020 [cited on Apr. 10, 2022]. Available at: https://geodacenter.github.io/workbook/6a_local_auto/lab6a.html
    » https://geodacenter.github.io/workbook/6a_local_auto/lab6a.html
  • 27
    Lima Filho CA, Lima AES, Arcanjo RMG, Silva DL, Jesus GF, Albuquerque AOBC, et al. Perfil epidemiológico dos casos de dengue no estado de Pernambuco, Brasil. Research, Society and Development 2022; 11(2): e36711225891. http://dx.doi.org/10.33448/rsd-v11i2.25891
    » http://dx.doi.org/10.33448/rsd-v11i2.25891
  • 28
    Freitas JR, Nascimento GILA, Ferreira DSA, Santiago EJP, Moreira GR, Silva ASA, et al. Análise espaço-temporal da incidência de febre Chikungunya no estado de Pernambuco. Research, Society and Development 2020; 9(9): e288997114. http://dx.doi.org/10.33448/rsd-v9i9.7114
    » http://dx.doi.org/10.33448/rsd-v9i9.7114
  • 29
    Brasil. Ministério da Saúde. Características dos Indicadores. Fichas de Qualificação. Taxa de incidência de dengue – D.2.3 – 2012 [Internet]. [cited on Jun. 30, 2022]. Available at: http://fichas.ripsa.org.br/2012/d-2-3/?l=pt_BR
    » http://fichas.ripsa.org.br/2012/d-2-3/?l=pt_BR
  • 30
    Azevedo TS, Lorenz C, Chiaravalloti-Neto F. Spatiotemporal evolution of dengue outbreaks in Brazil. Trans R Soc Trop Med Hyg 2020; 114(8): 593-602. https://doi.org/10.1093/trstmh/traa030
    » https://doi.org/10.1093/trstmh/traa030
  • 31
    Skalinski LM, Dalvi APR, Natividade M, Braga JU, Costa MCN, Santana EB, et al. The triple epidemics of arboviruses in Feira de Santana, Brazilian Northeast: epidemiological characteristics and diffusion patterns. Epidemics. 2022; 38: 100541. https://doi.org/10.1016/j.epidem.2022.100541
    » https://doi.org/10.1016/j.epidem.2022.100541
  • 32
    Salata AR, Ribeiro MG. Boletim Desigualdade nas Metrópoles [Internet] 2021. Porto Alegre: Observatório das Metrópoles. Pontifícia Universidade Católica do Rio Grande do Sul; 2021. [cited on Jul. 1, 2022]. Available at: https://www.pucrs.br/datasocial/wp-content/uploads/sites/300/2022/07/boletim_desigualdade_nas_metropoles_06.pdf
    » https://www.pucrs.br/datasocial/wp-content/uploads/sites/300/2022/07/boletim_desigualdade_nas_metropoles_06.pdf
  • 33
    Instituto Brasileiro de Geografia e Estatística. Base de dados – Munic 2017 – saneamento básico [Internet]. [cited on Jun. 28, 2022]. Available at: https://censoagro2017.ibge.gov.br/agencia-detalhe-de-midia.html?view=mediaibge&catid=2103&id=2341
    » https://censoagro2017.ibge.gov.br/agencia-detalhe-de-midia.html?view=mediaibge&catid=2103&id=2341
  • 34
    Ribeiro MS, Ferreira DF, Azevedo RC, Santos GBG, Medronho RA. Índices larvais de Aedes aegypti e incidência de dengue: um estudo ecológico no Estado do Rio de Janeiro, Brasil. Cad Saúde Pública 2021; 37(7): e00263320. https://doi.org/10.1590/0102-311X00263320
    » https://doi.org/10.1590/0102-311X00263320
  • 35
    Lisboa TR, Serafim IBM, Serafim JCM, Ramos AC, Nascimento RM, Roner MNB. Relação entre incidência de casos de arboviroses e a pandemia da COVID-19. Rev Interdiscip Ciênc Apl 2022; 6(10): 31-6. https://doi.org/10.18226/25253824.v6.n10.04
    » https://doi.org/10.18226/25253824.v6.n10.04
  • 36
    Brasil. Ministério da Saúde. Nota Informativa n° 8/2020. Recomendações aos Agentes de Combate a Endemias (ACE) para adequação das ações de vigilância e controle de zoonoses frente à atual situação epidemiológica referente ao Coronavírus (COVID-19) [Internet]. Brasília: Ministério da Saúde; 2020. [cited on Oct. 28, 2021]. Available at: https://www.unasus.gov.br/especial/covid19/pdf/106
    » https://www.unasus.gov.br/especial/covid19/pdf/106
  • 37
    Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med 2020; 382(18): 1708-20. https://doi.org/10.1056/NEJMoa2002032
    » https://doi.org/10.1056/NEJMoa2002032
  • 38
    Rodriguez-Morales AJ, Villamil-Gómez WE, Franco-Paredes C. The arboviral burden of disease caused by co-circulation and co-infection of dengue, chikungunya and Zika in the Americas. Travel Med Infect Dis 2016; 14(3): 177-9. https://doi.org/10.1016/j.tmaid.2016.05.004
    » https://doi.org/10.1016/j.tmaid.2016.05.004

Publication Dates

  • Publication in this collection
    20 Feb 2023
  • Date of issue
    2023

History

  • Received
    11 Aug 2022
  • Reviewed
    04 Nov 2022
  • Accepted
    18 Nov 2022
Associação Brasileira de Pós -Graduação em Saúde Coletiva São Paulo - SP - Brazil
E-mail: revbrepi@usp.br