Acessibilidade / Reportar erro

Spatial diffusion of Zika fever epidemics in the Municipality of Salvador-Bahia, Brazil, in 2015-2016: does Zika fever have the same spread pattern as Dengue and Chikungunya fever epidemics?

Abstract

INTRODUCTION

The recent emergence and rapid spread of Zika and Chikungunya fevers in Brazil, occurring simultaneously to a Dengue fever epidemic, together represent major challenges to public health authorities. This study aimed to identify and compare the 2015-2016 spatial diffusion pattern of Zika, Chikungunya, and Dengue epidemics in Salvador-Bahia.

METHODS

We used two study designs comprising a cross-sectional-to-point pattern and an ecological analysis of lattice data. Residential addresses involving notified cases were geocoded. We used four spatial diffusion analysis techniques: (i) visual inspection of the sequential kernel and choropleth map, (ii) spatial correlogram analysis, (iii) spatial local autocorrelation (LISA) changes analysis and, (iv) nearest neighbor index (NNI) modeling.

RESULTS

Kernel and choropleth maps indicated that arboviruses spread to neighboring areas near the first reported cases and occupied these new areas, suggesting a diffusion expansion pattern. A greater case density occurred in central and western areas. In 2015 and 2016, the NNI best-fit model had an S-curve compatible with an expansion pattern for Zika (R2 = 0.94; 0.95), Chikungunya (R2 = 0.99; 0.98) and Dengue (R2 = 0.93; 0.99) epidemics, respectively. Spatial correlograms indicated a decline in spatial lag autocorrelations for the three diseases (expansion pattern). Significant LISA changes suggested different diffusion patterns, although a small number of changes were detected.

CONCLUSIONS

These findings indicate diffusion expansion, a unique spatial diffusion pattern of Zika, Chikungunya, and Dengue epidemics in Salvador-Bahia, namely. Knowing how and where arboviruses spread in Salvador-Bahia can help improve subsequent specific epidemic control interventions.

Keywords:
Zika; Chikungunya; Dengue; Arboviruses; Spatial diffusion

INTRODUCTION

In recent years, the emergence of arboviruses in different territories and countries of the Americas, including Brazil, has been observed11. Lima-Camara TN. Emerging arboviruses and public health challenges in Brazil. Rev Saúde Pública [Internet]. 2016;50:1-7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27355468
http://www.ncbi.nlm.nih.gov/pubmed/27355...
. The most recent arbovirus identified in Brazil was Zika fever in 2015, which triggered an epidemic of undetermined exanthematous disease in the northeast region22. Campos G, Bandeira A, Sardi S. Zika Virus Outbreak, Bahia Brazil. Emerg Infect Dis. 2015;21(10):1881..

The disease presented in most individuals with a mild and benign evolution; however, a concentration of Guillain-Barré syndrome cases generated hypotheses that the Zika virus infection was associated with neurological and autoimmune complications in Brazil, as had also been suspected in French Polynesia in 201333. Cao-Lormeau V-M, Blake A, Mons S, Lastère S, Roche C, Vanhomwegen J, et al. Guillain-Barré Syndrome outbreak associated with Zika virus infection in French Polynesia: a case-control study. Lancet [Internet]. 2016;387(10027):1531-9. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0140673616005626
http://linkinghub.elsevier.com/retrieve/...
,44. Cardoso CW, Paploski IAD, Kikuti M, Rodrigues MS, Silva MMO, Campos GS, et al. Outbreak of exanthematous illness associated with Zika, Chikungunya, and Dengue viruses, Salvador, Brazil. Emerg Infect Dis . 2015;21(12):2274-6.. Moreover, a considerable increase in the number of reported cases of microcephaly in newborns provided evidence supporting a connection between Zika fever and congenital neurological malformations55. Martines RB, Bhatnagar J, Keating MK, Silva-Flannery L, Muehlenbachs A, Gary J, et al. Notes from the field: evidence of Zika virus infection in brain and placental tissues from two congenitally infected newborns and two fetal losses--Brazil, 2015. MMWR Morb Mortal Wkly Rep. 2016;65(6):159-60..

Prior to the sudden rise in Zika fever cases, Chikungunya fever emerged in Brazil in 201466. Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde - Relatório de Situação: Brasília. 2014. V. 45, N. 22. 6 p.. Along with its severe articular pain characteristic, severe manifestations such as neurological, cardiac, renal, and ocular complications have been reported to be associated with the disease77. Rajapakse S, Rodrigo C, Rajapakse A. Atypical manifestations of chikungunya infection. Trans R Soc Trop Med Hyg. 2010;104(2):89-96..

These two arboviruses emerged amidst an endemic Dengue fever circulation present since 1986, involving four viral serotypes and considered one of the most prevalent infectious diseases. A triple epidemic in 2016 resulted in approximately 2 million reported cases88. Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde - Relatório de Situação: Brasília. 2017. V. 48, N. 7. 10 p.,99. Osanai CH, Travassos da Rosa APA, Tang AT, do Amaral RS, Passos AD, Tauil PL. Surto De Dengue Em Boa Vista, Roraima. Nota Previa. Rev Inst Med Trop Sao Paulo. 1983;25(1):53-4..

These arboviruses have several characteristics in common; transmission occurs through the same vector, namely, the Aedes aegypti mosquito, which is abundant throughout Brazil, and these arboviruses are derived from very close families that have similarities in signs and symptoms, making diagnosis difficult1010. Gubler DJ. Dengue/dengue haemorrhagic fever: history and current status. Novartis Found Symp [Internet]. 2016;277:3-16. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17319151
http://www.ncbi.nlm.nih.gov/pubmed/17319...
,1111. Ioos S, Mallet HP, Leparc Goffart I, Gauthier V, Cardoso T, Herida M. Current Zika virus epidemiology and recent epidemics. Med Mal Infect [Internet]. 2014;44(7):302-7. Available from: http://dx.doi.org/10.1016/j.medmal.2014.04.008
http://dx.doi.org/10.1016/j.medmal.2014....
. The vector control measures and the main strategy for disease control were inefficient and favored the propagation of arboviruses in Brazil1212. Yakob L, Walker T. Zika virus outbreak in the Americas: The need for novel mosquito control methods. Lancet Glob Heal [Internet]. 2016;4(3):e148-9. Available from: http://dx.doi.org/10.1016/S2214-109X(16)00048-6
https://doi.org/http://dx.doi.org/10.101...
.

Given this situation, an analysis of the geographical processes of spatial diffusion can be used for monitoring epidemics and planning interventions, since it enables an understanding of how and where diseases propagate through space and in the course of time1313. Cromley EK, McLafferty S. Analyzing the risk and spread of infectious diseases. In: GIS and Public Health. 2 ed. 2012. p. 234-62.. In 2014, Sant-Julien stated that the process of spatial diffusion is the action, or the result of an action of propagating a phenomenon homogeneously in a system, whatever the force that drives the dispersion of this phenomenon1414. Sant-Julien T. The process of spatial diffusion and modeling change. Model Spat Anal [Internet]. 2010;127-57. Available from: http://doi.wiley.com/10.1002/9780470612255.ch5
http://doi.wiley.com/10.1002/97804706122...
. The visualization of spatial patterns in data and the description of these patterns aims to predict changes and uses this information to guide the formulation of policies1515. Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, Clements A. Spatial analysis in epidemiology. Oxford University Press. 2008. 1-14 p..

In the most commonly used classification, it is stated that spatial diffusion can occur through four phenomena. The first is expansion, in which the phenomenon propagates from one place to the other, often intensifying in the region of origin. The second phenomenon is relocation, when the phenomenon propagates to new localities but when moving, leaves the areas where it originated. The third is contagion, dependent on direct contact and, therefore, strongly influenced according to distance (individuals in nearby regions have a much higher probability of contact than individuals in distant regions). Lastly, a hierarchical phenomenon is characterized by the occurrence of diffusion respecting an ordered sequence of classes or places, as from a large metropolitan center to remote villages1616. Cliff AD, Haggett P, Ord JK VG. Spatial diffusion. Cambridge University Press; 1981..

A very limited number of studies have sought to specifically analyze the spatial diffusion of arboviruses, with an emphasis on investigating Dengue fever in the state of Bahia, Brazil1717. Barreto FR, Teixeira MG, Costa MDCN, Carvalho MS, Barreto ML. Spread pattern of the first dengue epidemic in the city of Salvador, Brazil. BMC Public Health [Internet]. 2008;8(1):51. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-8-51
http://bmcpublichealth.biomedcentral.com...
,1818. Morato DG, Barreto FR, Braga JU, Natividade MS, Costa M da CN, Morato V, et al. The spatiotemporal trajectory of a dengue epidemic in a medium-sized city. Mem Inst Oswaldo Cruz [Internet]. 2015;110(4):528-33. Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0074-02762015000400528&lng=en&nrm=iso&tlng=en
http://www.scielo.br/scielo.php?script=s...
. This study aimed to identify and compare the 2015-2016 spatial diffusion pattern of Zika, Chikungunya, and Dengue epidemics in Salvador-Bahia.

METHODS

Two study designs were used, a sectional study for the point spatial data, and an ecologic study for the area data. The study area comprised the municipality of Salvador, the capital of state of Bahia which is the third most populous city in Brazil with an estimated population of 2,938,092 inhabitants in 20161919. IBGE. População do município de Salvador em 2016 [Internet]. Instituto brasileiro de geografia e estatística. 2016. Available from: http://cidades.ibge.gov.br/xtras/perfil.php?codmun=292740
http://cidades.ibge.gov.br/xtras/perfil....
and 163 neighborhoods (Figure 1)2020. SEDUR. Plano Diretor de Desenvolvimento Urbano do Município de Salvador [Internet]. Secretaria de desenvolvimento e urbanismo. 2016. Available from: http://www.sucom.ba.gov.br/wp-content/uploads/2016/07/LEI-n.-9.069-PDDU-2016.pdf
http://www.sucom.ba.gov.br/wp-content/up...
.

FIGURE 1:
Location of the municipality of Salvador in the state of Bahia and its territorial division.

In the data analysis, three island neighborhoods namely Ilha dos Frades, Ilha de Maré, and Ilha de Bom Jesus dos Passos were excluded, because we considered that, through discontinuity with the territory, the processes that act for spatial diffusion of arboviruses could have differed. Centro Administrativo da Bahia (CAB) and Aeroporto were also excluded since they are not residential neighborhoods, making 5 of 163 neighborhoods excluded from the study.

The study period covered the epidemic waves that occurred in the years 2015 and 2016. For Zika fever, the first wave comprised epidemiological weeks (EW) 15/2015 to 52/2015, and for Chikungunya and Dengue fever, the first wave comprised the EWs from 01/2015 to 52/2015. For the three arboviruses, we considered the second wave comprised the EWs from 01/2016 to 52/2016.

Confirmed cases were used, as notified in the Sistema de Informação de Agravos de Notificação-SINAN of the Brazilian Ministry of Health. Residential addresses for each case were georeferenced to generate punctual spatial data using three automated geocoding application programming interfaces to obtain a greater number of georeferenced cases: Google Maps, Bing Maps, and Open Street Maps.

For aggregated data, we opted for smoothed incidence rates (empirical Bayesian) because of the large variation of gross rates that would have strongly influenced the analyses. The number of people living in the neighborhoods of Salvador in 2015 was obtained from estimates according to the Brazilian Institute of Geography and Statistics (IBGE), based on a 2010 census, and a linear interpolation was performed to estimate the population for the year 20161919. IBGE. População do município de Salvador em 2016 [Internet]. Instituto brasileiro de geografia e estatística. 2016. Available from: http://cidades.ibge.gov.br/xtras/perfil.php?codmun=292740
http://cidades.ibge.gov.br/xtras/perfil....
,2121. STATACORP. Stata Statistical Software: Release 13.0. 2015.. Cartographic digital meshes of the municipalities and the neighborhoods were obtained from the IBGE.

An analysis of spatial diffusion patterns was performed using various techniques that have been previously applied in other studies, namely visual inspection of the sequential kernel and choropleth map, spatial correlogram analysis, local spatial autocorrelation (LISA) changes analysis, and the nearest neighbor index (NNI) modeling.

For visual inspection of sequential maps, 30 choropleth maps and 30 kernel maps were generated with data from every three EWs of Zika fever. For Chikungunya and Dengue fevers, 35 thematic and 35 kernel maps were generated.

Kernel estimation is a method for smoothing point events through inserting a continuous surface over each point and requires a bandwidth of the region of influence. A smoothed surface is displayed on the territory representing the levels of case intensity2222. Andrade ALSS de, Silva SA e, Martelli CMT, Oliveira RM de, Morais Neto OL de, Siqueira Júnior JB, et al. Population-based surveillance of pediatric pneumonia: use of spatial analysis in an urban area of Central Brazil. Cad saude publica / Minist da Saude, Fund Oswaldo Cruz, Esc Nac Saude Publica. 2004;20(2):411-21.. Kernel maps were developed with a bandwidth of 2500 m, using QGIS software2323. Team Qgis Development. QGIS Geographic Information System [Internet]. Open Source Geospatial Foundation Project. 2016. Available from: http://www.qgis.org
http://www.qgis.org...
.

Choropleth maps display levels of the same attribute per area at each time range for sequential demonstration of the distribution of the disease. The maps were constructed using five classes of incidence rates using QGIS software2323. Team Qgis Development. QGIS Geographic Information System [Internet]. Open Source Geospatial Foundation Project. 2016. Available from: http://www.qgis.org
http://www.qgis.org...
.

In the modeling of the NNI, the ratio of the nearest neighbor, or R, is calculated by dividing the mean of the distances of the nearest points in a point pattern by the average of the distances of a random distribution of the same number of points in the same area. R is calculated cumulatively whenever a new case is diagnosed. The Rs are presented as a function of time and modeled with a set of regression curves to find the one that best fits the distribution. The results were interpreted according to Lee et al. (2014)2424. Lee J, Lay J-G, Chin WCB, Chi Y-L, Hsueh Y-H. An experiment to model spatial diffusion process with nearest neighbor analysis and regression estimation. Int J Appl Geospatial Res [Internet]. 2014;5(1):1-15. Available from: http://services.igi-global.com/resolvedoi/resolve.aspx?doi=10.4018/ijagr.2014010101
http://services.igi-global.com/resolvedo...
. Their approach indicates that the processes of diffusion through contagion and expansion are better adjusted using an inverse curve, and the hierarchical and relocation processes are better modeled using a cubic curve.

In our study, linear, logarithmic, inverse, quadratic, cubic, power, compound, S-curve, logistic, growth, and exponential regression curves were constructed using the “curvefit” function of Stata software2121. STATACORP. Stata Statistical Software: Release 13.0. 2015.. These curves were compared to the annual R curves plotted against time so that the best fit was identified with the coefficient of R².

For the analysis of the spatial correlograms, the values ​​corresponding to Moran's I (global autocorrelation indicators) were calculated using an adjacency matrix (determining the spatial autocorrelation) and plotted against the spatial lags referring to the distance intervals. The correlograms were developed for each arbovirus annually, using “R statistical software”2525. R (Core Team). R: A language and environment for statistical computing. [Internet]. R Foundation for Statistical Computing. Vienna, Austria; 2016. Available from: https://www.r-project.org/
https://www.r-project.org/...
in conjunction with QGIS software2323. Team Qgis Development. QGIS Geographic Information System [Internet]. Open Source Geospatial Foundation Project. 2016. Available from: http://www.qgis.org
http://www.qgis.org...
.

The incidence of each arbovirus was considered autocorrelated in space using definitions proposed by Lam et al. (1996)2626. Lam NS, Fan M. Spatial-temporal spread of the AlDS epidemic , regions of the United States. 1996;28(2)., in which a propagation by contagion model is considered when the correlogram shows a tendency to autocorrelation in a soft decline, pointing to distance as a principal factor in determining the similarity of incidence rates. A hierarchical diffusion is suggested when the trajectory decreases and then rises, forming a V-shaped curve, indicating that distant areas show similarity in rates and that the diffusion then occurs in jumps.

To apply the technique of change in local autocorrelation, the Cohen and Tita (1999) study was used as a reference2727. Cohen J, Tita G. Diffusion in homicide : exploring a general method for detecting spatial diffusion processes. J Quant Criminol [Internet]. 1999;15(4). Available from: https://www.jstor.org/stable/23366752?seq=1#page_scan_tab_contents
https://www.jstor.org/stable/23366752?se...
. LISA indicators consist of the decomposition of Moran's I2828. Anselin L. Local Indicators of Spatial Association-LISA. 1995;27(2)..

The values ​​of a variable for a location and its neighboring areas provide the most basic representation of local spatial associations. Each Local-Neighbor (LN) pair consists of standardized levels of a variable in the local spatial unit L and in neighboring spatial units N. Each element of the pair is low (L) or high (H) relative to the local and neighbor value distributions in all observations2727. Cohen J, Tita G. Diffusion in homicide : exploring a general method for detecting spatial diffusion processes. J Quant Criminol [Internet]. 1999;15(4). Available from: https://www.jstor.org/stable/23366752?seq=1#page_scan_tab_contents
https://www.jstor.org/stable/23366752?se...
.

In a scatterplot, pairs in which both local and neighbor values ​​are above their respective averages fall in the upper right quadrant (HH), or into the lower left quadrant (LL). When they differ, pairs fall into the HL or LH quadrant.

For a dynamic view of the process, Cohen and Tita (1999)2727. Cohen J, Tita G. Diffusion in homicide : exploring a general method for detecting spatial diffusion processes. J Quant Criminol [Internet]. 1999;15(4). Available from: https://www.jstor.org/stable/23366752?seq=1#page_scan_tab_contents
https://www.jstor.org/stable/23366752?se...
suggested analyzing changes in local-neighbor pair levels over time, looking for evidence of diffusion that involves the propagation of high rates (or low rates) to other spatial units. The authors claimed that the types of pair changes in successive periods were compatible with a different type of diffusion, and allowed for a determination of the mechanisms behind the change corresponded to each type of diffusion.

LISA values were calculated using GeoDa software2929. Anselin L. Geoda: User’s guide. Spatial Analysis Laboratory (SAL). Department of Agricultural and Consumer Economics, University of Illinois, Urbana-Champaign, IL; 2017., based on the incidence rates of each arbovirus, with reference to periods of three EWs. Analyzed annually, changes in the predominant local autocorrelation levels determined the spatial diffusion type for each arbovirus.

RESULTS

The distribution curve of Zika fever cases showed two epidemic waves. The first occurred between EW 15 and EW 41 in 2015, with an epidemic peak at EW 27 (112 cases). The second less significant epidemic, occurred between EW 01 and EW 22 in 2016, with a higher concentration at EW 9 (43 cases) (Figure 2).

FIGURE 2:
Distribution curves of reported cases of Zika fever, Chikungunya fever and Dengue in Salvador, 2015-2016.

The distribution curve of Chikungunya fever cases shows an epidemic wave in 2015, with a higher case concentration between EW 26 and EW 38, and an epidemic peak at EW 33 (76 cases). Waves of lesser intensity were identified throughout 2016.

The distribution curve of Dengue fever cases shows a large epidemic wave that covered the entire year in 2015, with a higher concentration of cases between EW 14 and EW 39. During 2016, there was a regular level of reported cases.

The majority of reported cases involving the three arboviruses lacked confirmatory diagnostic criteria (undetermined), and diagnosis was presumably based on clinical-epidemiological criteria, as laboratory confirmation is usually performed at the beginning of epidemics to confirm autochthony.

The georeferencing of Zika fever cases in Salvador resulted in 1,780 (93%) georeferenced addresses for the 1,914 cases confirmed between 2015 and 2016. These cases were organized according to the date of first symptoms, and 32 cases were discarded due to the date of first symptoms field having been incorrectly completed. In total, 1,748 cases were considered in the statistical analysis and, of these, there were 937 cases in 2015 and 811 in 2016. Of 2,120 Chikungunya fever cases, 1,974 (93%) were georeferenced, 38 were discarded, and we used 1,210 cases from the year 2015 and 726 from the year 2016 in the statistical analysis. Of 9,302 cases of confirmed Dengue fever, 8,938 (96%) were georeferenced, 301 discarded, and we used 5,449 cases for the year 2015, and 3,188 for the year 2016.

Zika, Chikungunya, and Dengue fever cases were recorded in approximately 95% of neighborhoods analyzed in the municipality of Salvador. The analysis of kernel maps (Figure 3) and choropleth maps suggested that the spatial diffusion process of Zika fever occurred through expansion, as did Dengue and Chikungunya fevers. Some areas of the municipality, mainly in the west/southwest region, showed an emergence of cases that then intensified and reached nearby regions.

FIGURE 3:
Kernel maps for analysis of the spread of Zika fever in Salvador, EW 15-44/2015 (the maps of the last EW were hidden due their similarities).

Regarding Zika fever, the northeast of the municipality also presented a greater concentration of cases in 2015. In Chikungunya and Dengue fevers, the central region was highlighted. In 2016, the diseases appeared to have occurred with less intensity, with a greater concentration of cases observed in the central region.

For Zika fever, 144 and 311 R-values ​​were generated for the years 2015 and 2016, respectively. For Chikungunya fever, 285 and 303 R-values ​​were generated for the years 2015 and 2016, respectively. For Dengue fever, 357 and 414 R-values were generated for the years 2015 and 2016, respectively.

The regression curve that best fit the R curves was the S-curve (Figure 4), with R2 0.94;0.95 (Zika fever), R2 0.99;0.98 (Chikungunya fever), and R2 0.93;0.99 (Dengue fever) in 2015 and 2016, respectively; therefore, showing a process of diffusion through expansion.

FIGURE 4:
Regression curves of better fit to Nearest Neighbor Index (NNI) curves of Zika fever, Chikungunya and Dengue in Salvador, 2015-2016.

The correlograms of the three arboviruses show a smooth decline in spatial autocorrelation as the distance between spatial lags increases (Figure 5). This decline shows that expansion diffusion was the dominant factor for the spread of Zika, Chikungunya, and Dengue fevers in the municipality of Salvador. That is, the municipality presented areas with high incidence rates and neighboring areas that had lower incidence rates, indicating that these diseases progressively spread from one region to their neighboring regions.

FIGURE 5:
Spatial correlograms of Zika fever, Chikungunya and Dengue in Salvador, 2015 - 2016.

In the analysis of LISA changes, LISA values ​​were generated with reference to periods of three EWs, for each of the 160 neighborhoods within the Salvador municipality. Analyzing the changes in these values, 19 of 1,280 LISA changes were significant for Zika fever in 2015 and 34 of 2,400 LISA changes were significant for Zika fever in 2016. The local autocorrelation changes characteristic of expansion diffusion[14] were predominant in 2015 and hierarchical diffusion[20] in 2016. Relocation diffusion was considered likely here, along with hierarchical diffusion, due to similarity in the characteristics of the processes that occurred in jumps.

For Chikungunya fever, 42 out of 2,400 LISA changes were significant in 2015 and 19 out of 2,080 LISA changes were significant in 2016. In two years, the type of spatial diffusion was mixed because the number of changes related to the expansion and hierarchical diffusion were equal in 2015[21;21], and similar in 2016[10;19].

For Dengue fever in 2015, 62 out of 2,560 LISA changes were significant and 79 out of 2560 changes were significant in 2016. The number of changes characterizing the hierarchical diffusion process[34] surpassed the expansion[28] in 2015 but was very close in 2016[39;40] with diffusion mixed in this year.

DISCUSSION

The analysis of the spatial diffusion of Zika fever, Chikungunya fever, and Dengue fever in Salvador demonstrated that the propagation of the three arboviruses occurred through expansion in the two years of study.

In 2015, Zika and Chikungunya fevers broke out in Salvador and, as expected, due to population susceptibility and abundance of the main vector, Aedes aegypti, the arboviruses caused an explosive number of cases.

Although Dengue fever had already established in this period in Salvador, it showed a significant increase in notifications, which was observed throughout Brazil. Experts attribute this increase to the reemergence of the DENV1 serotype in 2013 and to a reduced DENV4 circulation, taking into account that the presence of the four serotypes of the Dengue virus in Brazil often trigger epidemics of the disease after alternating or replacing the predominant serotype3030. Teixeira MG, Siqueira JB, Ferreira GLC, Bricks L, Joint G. Epidemiological trends of dengue disease in Brazil (2000-2010): a systematic literature search and analysis. PLoS Negl Trop Dis. 2013;7(12).,3131. Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde. Plano de contingência nacional para epidemias de dengue. Brasília. 2015.. The antigenic diversity of the Dengue virus is a contributing factor to the continuous occurrence of the disease, since multiple sequential infections may occur due to the lack of cross-immunity between virus serotypes3232. Gubler DJ. Epidemic dengue / dengue hemorrhagic fever as a public health , social and economic problem in the 21st century. TRENDS Microbiol. 2002;10(2):100-3..

Our study demonstrated that the three diseases emerged and intensified in some areas of Salvador, and initially expanded to nearby sites through the expansion diffusion process. This finding is compatible with the dynamics of infectious diseases, where transmission is more likely to occur among individuals who are closely connected in space and time3030. Teixeira MG, Siqueira JB, Ferreira GLC, Bricks L, Joint G. Epidemiological trends of dengue disease in Brazil (2000-2010): a systematic literature search and analysis. PLoS Negl Trop Dis. 2013;7(12)..

For the year 2016, our results indicated that the three arboviruses had spread through expansion, although with less intensity. The progressive decrease in the number of new cases was consistent with the dynamics of the arboviruses in the second year of propagation because, although the number of new cases continued to intensify in places with greater population density, a decrease in the susceptible population in some areas led to a lower occurrence of diseases in the greater part of Salvador.

The diffusion through expansion model was identified by Barreto et al. (2008)1717. Barreto FR, Teixeira MG, Costa MDCN, Carvalho MS, Barreto ML. Spread pattern of the first dengue epidemic in the city of Salvador, Brazil. BMC Public Health [Internet]. 2008;8(1):51. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-8-51
http://bmcpublichealth.biomedcentral.com...
in relation to the first Dengue epidemic in Salvador in 1995, during the emergency period of this arbovirus. Lee et al. (2004)2424. Lee J, Lay J-G, Chin WCB, Chi Y-L, Hsueh Y-H. An experiment to model spatial diffusion process with nearest neighbor analysis and regression estimation. Int J Appl Geospatial Res [Internet]. 2014;5(1):1-15. Available from: http://services.igi-global.com/resolvedoi/resolve.aspx?doi=10.4018/ijagr.2014010101
http://services.igi-global.com/resolvedo...
identified the same pattern concerning the spatial diffusion of Dengue fever between 2003 and 2008 in Taiwan.

Morato et al. (2015)1818. Morato DG, Barreto FR, Braga JU, Natividade MS, Costa M da CN, Morato V, et al. The spatiotemporal trajectory of a dengue epidemic in a medium-sized city. Mem Inst Oswaldo Cruz [Internet]. 2015;110(4):528-33. Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0074-02762015000400528&lng=en&nrm=iso&tlng=en
http://www.scielo.br/scielo.php?script=s...
identified a different pattern in their study on the spatial diffusion of Dengue fever in the municipality of Jequié, Bahia, in 2009, where they showed that the disease epidemic occurred initially through expansion and later through relocation.

In addition to the identification of the spatial diffusion processes, analysis of kernel maps indicated that the three diseases had a main focus of transmission in an area that extended from the west/southwest region of the municipality, where administrative districts are concentrated, to the central region, with both regions having higher population density and social interaction networks. This area sustained a high number of cases throughout almost the entire study period, possibly due to the higher concentration of susceptible cases compared to other areas.

Areas with high population density and economic development that can act as transmission focuses reinforce the pattern of diffusion through expansion as a characteristic of communicable diseases, reaching a greater number of people who are interacting in close quarters1515. Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, Clements A. Spatial analysis in epidemiology. Oxford University Press. 2008. 1-14 p..

These areas are compatible with the epicenter of the first Dengue fever epidemic, also occurring in Salvador, as reported by Barreto et al. (2008)1717. Barreto FR, Teixeira MG, Costa MDCN, Carvalho MS, Barreto ML. Spread pattern of the first dengue epidemic in the city of Salvador, Brazil. BMC Public Health [Internet]. 2008;8(1):51. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-8-51
http://bmcpublichealth.biomedcentral.com...
. This similarity suggests that, if interventions had been directed to those areas of previous Dengue fever epidemics, the expansion of arboviruses and the recent triple epidemic could have been minimized.

Among the visualization methods, the kernel maps provided a better appreciation of the spatial diffusion than the choropleth maps, because they exhibited areas with a higher density of cases where the risk of arbovirus transmission may have been higher due to prevalent cases.

Choropleth maps were more difficult to interpret because of excessive variation in incidence rates, a phenomenon less frequent in large populations and more important when small numbers of cases are related to small populations. The color scale has been reported to be also very sensitive to variations, which hinders interpretation1515. Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, Clements A. Spatial analysis in epidemiology. Oxford University Press. 2008. 1-14 p.,3333. Gatrell AC, Bailey TC. Interactive spatial data analysis in medical geography. Soc Sci Med. 1996;42(6):843-55.. Apart from this, the diffusion processes were evaluated in conjunction with other techniques used.

Despite a high degree of subjectivity, the use of visual presentation methods combined with statistical techniques can make quantitative results easier to interpret.

Among the statistical methods used, the LISA changes analysis technique produced a different result for the three arboviruses compared to the other techniques. However, the LISA result should be interpreted with caution because of the small proportion of statistically significant changes found in our study (<3%).

The spatial correlograms presented similar configurations. The NNI analysis technique appeared to be the most robust method of analysis as it was based on a statistical approach that used daily information, was modeled using regression curves, and could be compared to other analyses previously simulated by Lee et al. (2008)2424. Lee J, Lay J-G, Chin WCB, Chi Y-L, Hsueh Y-H. An experiment to model spatial diffusion process with nearest neighbor analysis and regression estimation. Int J Appl Geospatial Res [Internet]. 2014;5(1):1-15. Available from: http://services.igi-global.com/resolvedoi/resolve.aspx?doi=10.4018/ijagr.2014010101
http://services.igi-global.com/resolvedo...
.

The main limitation of this study is that the surveillance data may have contained a poor diagnostic classification in reports, possibly due to the rapid emergence of arboviruses and similarity in the symptoms when defining the diagnoses according to clinical epidemiological criteria. Simultaneous circulation of the three arboviruses was a factor that may have compromised the investigation of some of the diseases and led to more intense investigation of others. Since most cases were asymptomatic, we considered that a substantially lower number of cases might have been reported. It is likely that diagnostic cross-classification errors occurred during the simultaneous triple epidemic, but it is difficult to evaluate the magnitude of such bias. It might be that the cross-diagnosis brings the specific spatial diffusion patterns close to that observed for Dengue. This would obscure specific characteristics that might happen for Zika and Chikungunya.

Despite this limitation, the combination of statistical and visualization techniques appeared to allow the analysis of spatial diffusion. This information should allow the implementation of more effective control programs by directing prevention and control strategies to specific areas and periods, as recommended by Cliff et al. (1981)1616. Cliff AD, Haggett P, Ord JK VG. Spatial diffusion. Cambridge University Press; 1981..

Our study of the spatial diffusion of Zika, Chikungunya, and Dengue arboviruses in in Salvador-Bahia, in the period of triple epidemics between 2015 and 2016, enabled an identification of the propagation process of these arboviruses, indicating that the diffusion pattern of the Zika fever epidemics was the same as that of the Dengue and Chikungunya fever epidemics. During the two years of study, these diseases spread through diffusion expansion. This knowledge on how and where arboviruses spread in Salvador in a triple epidemic situation provides important information for improving epidemiological surveillance, both in terms of monitoring cases and in terms of targeting interventions for disease control.

REFERENCES

  • 1
    Lima-Camara TN. Emerging arboviruses and public health challenges in Brazil. Rev Saúde Pública [Internet]. 2016;50:1-7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27355468
    » http://www.ncbi.nlm.nih.gov/pubmed/27355468
  • 2
    Campos G, Bandeira A, Sardi S. Zika Virus Outbreak, Bahia Brazil. Emerg Infect Dis. 2015;21(10):1881.
  • 3
    Cao-Lormeau V-M, Blake A, Mons S, Lastère S, Roche C, Vanhomwegen J, et al. Guillain-Barré Syndrome outbreak associated with Zika virus infection in French Polynesia: a case-control study. Lancet [Internet]. 2016;387(10027):1531-9. Available from: http://linkinghub.elsevier.com/retrieve/pii/S0140673616005626
    » http://linkinghub.elsevier.com/retrieve/pii/S0140673616005626
  • 4
    Cardoso CW, Paploski IAD, Kikuti M, Rodrigues MS, Silva MMO, Campos GS, et al. Outbreak of exanthematous illness associated with Zika, Chikungunya, and Dengue viruses, Salvador, Brazil. Emerg Infect Dis . 2015;21(12):2274-6.
  • 5
    Martines RB, Bhatnagar J, Keating MK, Silva-Flannery L, Muehlenbachs A, Gary J, et al. Notes from the field: evidence of Zika virus infection in brain and placental tissues from two congenitally infected newborns and two fetal losses--Brazil, 2015. MMWR Morb Mortal Wkly Rep. 2016;65(6):159-60.
  • 6
    Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde - Relatório de Situação: Brasília. 2014. V. 45, N. 22. 6 p.
  • 7
    Rajapakse S, Rodrigo C, Rajapakse A. Atypical manifestations of chikungunya infection. Trans R Soc Trop Med Hyg. 2010;104(2):89-96.
  • 8
    Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde - Relatório de Situação: Brasília. 2017. V. 48, N. 7. 10 p.
  • 9
    Osanai CH, Travassos da Rosa APA, Tang AT, do Amaral RS, Passos AD, Tauil PL. Surto De Dengue Em Boa Vista, Roraima. Nota Previa. Rev Inst Med Trop Sao Paulo. 1983;25(1):53-4.
  • 10
    Gubler DJ. Dengue/dengue haemorrhagic fever: history and current status. Novartis Found Symp [Internet]. 2016;277:3-16. Available from: http://www.ncbi.nlm.nih.gov/pubmed/17319151
    » http://www.ncbi.nlm.nih.gov/pubmed/17319151
  • 11
    Ioos S, Mallet HP, Leparc Goffart I, Gauthier V, Cardoso T, Herida M. Current Zika virus epidemiology and recent epidemics. Med Mal Infect [Internet]. 2014;44(7):302-7. Available from: http://dx.doi.org/10.1016/j.medmal.2014.04.008
    » http://dx.doi.org/10.1016/j.medmal.2014.04.008
  • 12
    Yakob L, Walker T. Zika virus outbreak in the Americas: The need for novel mosquito control methods. Lancet Glob Heal [Internet]. 2016;4(3):e148-9. Available from: http://dx.doi.org/10.1016/S2214-109X(16)00048-6
    » https://doi.org/http://dx.doi.org/10.1016/S2214-109X(16)00048-6
  • 13
    Cromley EK, McLafferty S. Analyzing the risk and spread of infectious diseases. In: GIS and Public Health. 2 ed. 2012. p. 234-62.
  • 14
    Sant-Julien T. The process of spatial diffusion and modeling change. Model Spat Anal [Internet]. 2010;127-57. Available from: http://doi.wiley.com/10.1002/9780470612255.ch5
    » http://doi.wiley.com/10.1002/9780470612255.ch5
  • 15
    Pfeiffer D, Robinson T, Stevenson M, Stevens K, Rogers D, Clements A. Spatial analysis in epidemiology. Oxford University Press. 2008. 1-14 p.
  • 16
    Cliff AD, Haggett P, Ord JK VG. Spatial diffusion. Cambridge University Press; 1981.
  • 17
    Barreto FR, Teixeira MG, Costa MDCN, Carvalho MS, Barreto ML. Spread pattern of the first dengue epidemic in the city of Salvador, Brazil. BMC Public Health [Internet]. 2008;8(1):51. Available from: http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-8-51
    » http://bmcpublichealth.biomedcentral.com/articles/10.1186/1471-2458-8-51
  • 18
    Morato DG, Barreto FR, Braga JU, Natividade MS, Costa M da CN, Morato V, et al. The spatiotemporal trajectory of a dengue epidemic in a medium-sized city. Mem Inst Oswaldo Cruz [Internet]. 2015;110(4):528-33. Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0074-02762015000400528&lng=en&nrm=iso&tlng=en
    » http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0074-02762015000400528&lng=en&nrm=iso&tlng=en
  • 19
    IBGE. População do município de Salvador em 2016 [Internet]. Instituto brasileiro de geografia e estatística. 2016. Available from: http://cidades.ibge.gov.br/xtras/perfil.php?codmun=292740
    » http://cidades.ibge.gov.br/xtras/perfil.php?codmun=292740
  • 20
    SEDUR. Plano Diretor de Desenvolvimento Urbano do Município de Salvador [Internet]. Secretaria de desenvolvimento e urbanismo. 2016. Available from: http://www.sucom.ba.gov.br/wp-content/uploads/2016/07/LEI-n.-9.069-PDDU-2016.pdf
    » http://www.sucom.ba.gov.br/wp-content/uploads/2016/07/LEI-n.-9.069-PDDU-2016.pdf
  • 21
    STATACORP. Stata Statistical Software: Release 13.0. 2015.
  • 22
    Andrade ALSS de, Silva SA e, Martelli CMT, Oliveira RM de, Morais Neto OL de, Siqueira Júnior JB, et al. Population-based surveillance of pediatric pneumonia: use of spatial analysis in an urban area of Central Brazil. Cad saude publica / Minist da Saude, Fund Oswaldo Cruz, Esc Nac Saude Publica. 2004;20(2):411-21.
  • 23
    Team Qgis Development. QGIS Geographic Information System [Internet]. Open Source Geospatial Foundation Project. 2016. Available from: http://www.qgis.org
    » http://www.qgis.org
  • 24
    Lee J, Lay J-G, Chin WCB, Chi Y-L, Hsueh Y-H. An experiment to model spatial diffusion process with nearest neighbor analysis and regression estimation. Int J Appl Geospatial Res [Internet]. 2014;5(1):1-15. Available from: http://services.igi-global.com/resolvedoi/resolve.aspx?doi=10.4018/ijagr.2014010101
    » http://services.igi-global.com/resolvedoi/resolve.aspx?doi=10.4018/ijagr.2014010101
  • 25
    R (Core Team). R: A language and environment for statistical computing. [Internet]. R Foundation for Statistical Computing. Vienna, Austria; 2016. Available from: https://www.r-project.org/
    » https://www.r-project.org/
  • 26
    Lam NS, Fan M. Spatial-temporal spread of the AlDS epidemic , regions of the United States. 1996;28(2).
  • 27
    Cohen J, Tita G. Diffusion in homicide : exploring a general method for detecting spatial diffusion processes. J Quant Criminol [Internet]. 1999;15(4). Available from: https://www.jstor.org/stable/23366752?seq=1#page_scan_tab_contents
    » https://www.jstor.org/stable/23366752?seq=1#page_scan_tab_contents
  • 28
    Anselin L. Local Indicators of Spatial Association-LISA. 1995;27(2).
  • 29
    Anselin L. Geoda: User’s guide. Spatial Analysis Laboratory (SAL). Department of Agricultural and Consumer Economics, University of Illinois, Urbana-Champaign, IL; 2017.
  • 30
    Teixeira MG, Siqueira JB, Ferreira GLC, Bricks L, Joint G. Epidemiological trends of dengue disease in Brazil (2000-2010): a systematic literature search and analysis. PLoS Negl Trop Dis. 2013;7(12).
  • 31
    Ministério da Saúde (MS). Secretaria de Vigilância em Saúde. Sistema Nacional de Vigilância em Saúde. Plano de contingência nacional para epidemias de dengue. Brasília. 2015.
  • 32
    Gubler DJ. Epidemic dengue / dengue hemorrhagic fever as a public health , social and economic problem in the 21st century. TRENDS Microbiol. 2002;10(2):100-3.
  • 33
    Gatrell AC, Bailey TC. Interactive spatial data analysis in medical geography. Soc Sci Med. 1996;42(6):843-55.
  • Financial Support: This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES).

Publication Dates

  • Publication in this collection
    03 Apr 2020
  • Date of issue
    2020

History

  • Received
    20 Dec 2019
  • Accepted
    17 Feb 2020
Sociedade Brasileira de Medicina Tropical - SBMT Sociedade Brasileira de Medicina Tropical - SBMT, Núcleo de Medicina Tropical – UnB, Sala 43C – 70904-970, E-mails: rsbmt@uftm.edu.br | artes.rsbmt@gmail.com | sbmt@sbmt.org.br , WhatsApp: SBMT (61) 9.9192-6496, WhatsApp: RSBMT (34) 9.9996-5807 - Brasília - DF - Brazil
E-mail: rsbmt@uftm.edu.br