ABSTRACT

Objective:

To analyze the spatio-temporal dynamics of COVID-19 in the Rio de Janeiro state within the nine health regions, between March 2020 and December 2022.

Methods:

The Poisson model with random effects was used to smooth and estimate the incidence of COVID-19 hospitalizations reported in the Influenza Epidemiological Surveillance Information System (SIVEP-Gripe) to verify the synchronicity of the epidemic in the state.

Results:

The COVID-19 epidemic in the state is characterized by the presence of seven peaks during the analyzed period corresponding to seven found. An asynchrony in hospitalizations was identified, varying according to the different virus variants in the nine health regions of the state. The incidence peaks of hospitalizations ranged from 1 to 12 cases per 100,000 inhabitants during the pandemic.

Conclusion:

This spatio-temporal analysis is applicable to other scenarios, enabling monitoring and decision-making for the control of epidemic diseases in different areas.

Keywords:
COVID-19; Spatio-temporal analysis; Hospitalization; Rio de Janeiro

INTRODUCTION

COVID-19 (coronavirus disease 2019) is an acute and potentially severe respiratory infection, of great importance for public health¹1. Hasoksuz M, Kiliç S, Saraç F. Coronaviruses and SARS-COV-2. Turk J Med Sci 2020; 50(Sl-1): 549-56. https://doi.org/10.3906/sag-2004-127
https://doi.org/10.3906/sag-2004-127... , which has been declared by the World Health Organization (WHO) as a public health emergency of international concern (PHEIC) for three years²2. Organização Pan-Americana da Saúde. Organização Mundial da Saúde. OMS declara fim da Emergência de Saúde Pública de importância internacional referente à COVID-19 [Internet]. 2023 [acessado em 05 jun. 2023]. Disponível em: Disponível em: https://www.paho.org/pt/noticias/5-5-2023-oms-declara-fim-da-emergencia-saude-publica-importancia-internacional-referente
https://www.paho.org/pt/noticias/5-5-202... . COVID-19 is caused by the SARS-CoV-2 virus¹1. Hasoksuz M, Kiliç S, Saraç F. Coronaviruses and SARS-COV-2. Turk J Med Sci 2020; 50(Sl-1): 549-56. https://doi.org/10.3906/sag-2004-127
https://doi.org/10.3906/sag-2004-127... and may range according to severity, including asymptomatic cases and mild clinical manifestations, moderate cases in the form of pneumonia, as well as severe and critical cases, such as Severe A cute Respiratory Illness(SARI ), leading to the need for hospitalization³3. Iser BPM, Silva I, Raymundo VT, Poleto MB, Schuelter-Trevisol F, Bobinski F. Definição de caso suspeito da COVID-19: uma revisão narrativa dos sinais e sintomas mais frequentes entre os casos confirmados. Epidemiol Serv Saúde 2020; 29(3): e2020233. https://doi.org/10.5123/S1679-49742020000300018
https://doi.org/10.5123/S1679-4974202000... . This situation generated a high number of cases in a short space of time around the world, increasing the challenge for health systems⁴4. Brasil. Ministério da Saúde. Secretaria de Vigilância em Saúde. Boletim Epidemiológico Especial. Doença pelo Coronavírus COVID-19 [Internet]. 2021 [acessado em 03 abr. 2023]. Disponível em: Disponível em: https://www.gov.br/saude/pt-br/centrais-de-conteudo/publicacoes/boletins/epidemiologicos/covid-19/2021/boletim_epidemiologico_covid_59.pdf/view
https://www.gov.br/saude/pt-br/centrais-... .

The SARS-CoV-2 virus has suffered mutations throughout the pandemic, which generated new variants, since it has extensively circulated in the population, causing many infections⁵5. Valverde R. O que são mutações, linhagens, cepas e variantes? [Internet]. 2021 [acessado em 05 jun. 2023]. Disponível em: Disponível em: https://portal.fiocruz.br/noticia/o-que-sao-mutacoes-linhagens-cepas-e-variantes-1
https://portal.fiocruz.br/noticia/o-que-... . The first identified variants were Gama (P.1), in December, 2020, in Manaus, capital of Amazonas, and Zeta (P.2), identified in October, 2020, in the Rio de Janeiro state. Respectively, they were classified as variants of concern (VOC) and interest (VOI), according to the WHO classification⁶6. Michelon CM. Principais variantes do SARS-CoV-2 notificadas no Brasil. Rev Bras Anal Clin 2021; 53(2): 109-16. https://doi.org/10.21877/2448-3877.202100961
https://doi.org/10.21877/2448-3877.20210... ^,77. Organização Pan-Americana da Saúde. Organização Mundial da Saúde. OMS anuncia nomenclaturas simples e fáceis de pronunciar para variantes de interesse e de preocupação do SARS-CoV-2 [Internet]. 2021 [acessado em 05 jun. 2023]. Disponível em: Disponível em: https://www.paho.org/pt/noticias/1-6-2021-oms-anuncia-nomenclaturas-simples-e-faceis-pronunciar-para-variantes-interesse-e
https://www.paho.org/pt/noticias/1-6-202... . The epidemic has spread fast and became worse with the disorganization of the health system, as well as the lack of efficient public policies for control⁸8. Coelho FC, Lana RM, Cruz OG, Villela DAM, Bastos LS, Piontti APY, et al. Assessing the spread of COVID-19 in Brazil: mobility, morbidity and social vulnerability. PLoS One 2020; 15(9): e0238214. https://doi.org/10.1371/journal.pone.0238214
https://doi.org/10.1371/journal.pone.023... .

The Rio de Janeiro statewas one of the main viral dissemination poles for other Brazilian locations. During the initial phase of the pandemic in the country, SARS-CoV-2 first spread locally and within state borders. The motion of the virus lineage in the state was more frequent than that in other states⁵5. Valverde R. O que são mutações, linhagens, cepas e variantes? [Internet]. 2021 [acessado em 05 jun. 2023]. Disponível em: Disponível em: https://portal.fiocruz.br/noticia/o-que-sao-mutacoes-linhagens-cepas-e-variantes-1
https://portal.fiocruz.br/noticia/o-que-... . Without a nationally coordinated strategy, states and cities varied in form, intensity and duration of implementation of pharmacological and non-pharmacological measures⁹9. Pereira AK, Oliveira MS, Sampaio TS. Heterogeneidades das políticas estaduais de distanciamento social diante da COVID-19: aspectos políticos e técnico-administrativos. Rev Adm Pública 2020; 54(4): 678-96. https://doi.org/10.1590/0034-761220200323
https://doi.org/10.1590/0034-76122020032... . Figure 1 shows, chronologically, the moment when variants entered the state.

Considering the importance of knowing about the dynamics of disease transmission in one of the main transmission poles in the country, the objective of this study was to characterize the synchronicity of spatiotemporal dynamics of the COVID-19 epidemic in the Rio de Janeiro state . The nine health regions of the state were chosen as units of analysis from March 2020 to December 2022, using hospitalization data from the Influenza Epidemiological Surveillance Information System (Sivep-Gripe). The analysis provided relevant information about the occurrence of an epidemic in the state, with potential to be used in other disease situations involving epidemic characteristics. The presented results are useful as a surveillance protocol for future conditions. Therefore, the analysis is an important tool to characterize epidemics in different Brazilian regions and around the world, increasing the knowledge in public health and promoting the preparation for the emergence of new variants of SARS-CoV-2 , and emergence or re-emergence of other communicable diseases.

METHODS

Study area

The study area contemplates the entire Rio de Janeiro state, which is located in the Southeast region, with an area of 43,781,588 km². The state has 92 cities distributed in nine health regions: Ilha Grande Bay, Baixada Litorânea, Center-South, Médio Paraíba, Metropolitan I, Metropolitan II, Northwest, North and Serrana (Figure 2).

Data

The data regarding cases of patients hospitalized due to COVID-19 used in this study were obtained in secondary bases from the Ministry of Health, via Sivep-Gripe, which includes information about SARI. These data are available in OpenDataSUS (https://opendatasus.saude.gov.br/).

Data referring to population size of the health regions came from aggregations of population projection numbers per city, carried out by the Laboratory of Estimations and Population Projections (LEPP) in the Postgraduate program in Demography (PPGDem), of the Department of Demography and Actuarial Sciences (DDCA) of Universidade Federal do Rio Grande do Norte (UFRN). The projections are available at (https://demografiaufrn.net/projecao-populacional/).

Case definition

A case of hospitalization due to SARI-COVID was defined by a hospitalized individual notified in Sivep-Gripe, who tested positive for SARS-CoV-2 in the RT-PCR molecular test and (or) reactive to the antigen test, and (or) presenting epidemiological clinical criteria.

Data Analysis

To verify the synchronicity of the epidemic in the Rio de Janeiro state , we used time series of SARI-COVID cases in the nine health regions of the state, on the ninth epidemiological week of 2020 (from February 23 to 29, 2020), until the 52^nd epidemiological week of 2022 (from December 25 to 31, 2022). The expected incidence per week and health region, (λ ^(t,r)), was estimated based on the following random effect Poisson model:

• t = 1, 2, …
• T = 159 weeks
• r = 1, 2, …, 9 health regions

In this formula, H ^(t,r) is the number of SARI-COVID cases on week t, and health region r, POP^(r), the population of region r, β^(r) is a fixed common effect (intercept) for regional r, and random effects θ^(t,r) follow the Gaussian model with time dependence of second order random walk, that is, the effect of week t on region r depends only on what happened in the two weeks before in that region, and assuming that all health regions are conditionally independent, that is, do not have a spatial structure.

The weekly expected incidence, (λ ^(t,r)), was estimated for all health regions under the Bayesian approach, using the Inla method (Integrated Nested Laplace Approximations), which is a statistics technique used to allow the approximation of a posteriori distribution of complex hierarchical models in a fast and precise manner¹⁰10. Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Series B Stat Methodol 2009; 71(2): 319-92. https://doi.org/10.1111/j.1467-9868.2008.00700.x
https://doi.org/10.1111/j.1467-9868.2008... .

Time was divided according to waves, in which each wave is defined by a continuous period of time, with the domain of one variant. One variant is described as dominant while being the most frequent one in sequencings deposited in a virus sequencing repository. The dates of beginning and end of each wave were defined based on sequencing data of the Rio de Janeiro state , compiled by the genomic network of Fundação Oswaldo Cruz (https://www.genomahcov.fiocruz.br/). Since we are evaluating the synchronicity of the epidemic, the division of time in waves is the same for all health regions, and we used data from the entire Rio de Janeiro state . The moment when the expected incidence reaches its highest value was calculated in each wave and region. This value is called peak. Synchronicity was assessed based on Metropolitan Region I (MRI) as reference, which contains the state capital. If the peak week of a specific region coincided with the weak peak of MRI, then, in this wave, the epidemics were considered as synchronous; otherwise, they were asynchronous. The analysis was carried out in R¹¹11. R Core Team. A language and environment for statistical computing. r foundation for statistical computing [Internet]. Vienna; 2022 [acessado em 04 jul. 2023]. Disponível em: Disponível em: https://www.scirp.org/reference/referencespapers?referenceid=3456808
https://www.scirp.org/reference/referenc... .

RESULTS

Until March 19, 2020, only 23 days after the first confirmed case in Brazil, the Rio de Janeiro state had already registered 64 confirmed cases, spread in five cities: Rio de Janeiro (55), Niterói (6), Barra Mansa (1), Miguel Pereira (1) and Guapimirim (1). Five of the nine health regions already identified the circulation of the virus: Metropolitan I, Metropolitan II, Médio Paraíba, Center-South and Serrana¹²11. Governo do Estado do Rio de Janeiro. Secretaria de Saúde. Primeiro caso do Novo Coronavírus é confirmado no Estado do Rio [Internet]. 2020 [acessado em 03 jul. 2023]. Disponível em: Disponível em: https://www.saude.rj.gov.br:443/noticias/2020/03/primeiro-caso-do-novo-coronavirus-e-confirmado-no-estado-do-rio
https://www.saude.rj.gov.br:443/noticias... . The COVID-19 epidemic in the state can be characterized with seven waves in the analyzed period. These waves correspond to the following variants: B.1.1.33, Zeta, Gama, Delta and Omicron (subvariants BA.1, BA.2 and BA.5).

In Figure 3, it is observed that, in the first wave, the nine health regions presented different moments of increase in hospitalizations. The predominant variable in the state was B.1.1.33. Metropolitan Region I was the one with the highest incidence, approximately five cases per 100 thousand inhabitants. Variant B.1.133 did not present a well-defined peak in the following regions: Baixada Litorânea, Center-South and Serrana, where the incidence of hospitalizations remained high during the entire period; it was not possible to define a peak. In the second wave, corresponding to the Zeta variant, and in the third wave, corresponding to the Gama variant, the highest hospitalization rates were in the following regions: Baixada Litorânea, Center-South and Serrana. In both moments, incidence rates reached more than 2.5 cases per 100 thousand inhabitants in these regions. The Zeta variant did not present a well-defined peak only in the Northwest region, whereas the Gama variant had a defined peak in all nine regions of the state. In the fourth wave, the prevalent variant in the state was Delta, which did not present a well-defined peak in the regions of Ilha Grande Bay, Baixada Litorânea, Northwest and North. In this context, the Serrana region had higher incidence of hospitalizations, approximately five cases per 100 thousand inhabitants. The Omicron variant, which circulated in the state from December 2021 until the end of the study (December, 2022), is the only variant that presents three waves (subvariants BA.1, BA.2 and BA.5) in the nine regions; the ones with highest incidence rates were: Médio Paraíba, Metropolitan II and Serrana, ranging from 1 to 12 cases per 100 thousand inhabitants when Omicron was dominant in the state.

Figure 4 shows the difference in time between the peak moment of each region and the peak moment of Metropolitan Region I (E). The regions in grey in the map represent places that did not have well-defined hospitalization waves. For variant B.1.1.33, this occurred in Baixada Litorânea, Center-South and Serrana. The regions in the extremes of the state, Northwest and Médio Paraíba were the ones that took three weeks or more to have their peaks, followed by the Ilha Grande Bay, in two weeks, and North, in one week. There was synchronicity between Metropolitan regions I and II.

The Zeta variant was the one with most synchronicity, being observed in five regions: Metropolitan I and II, Baixada Litorânea, Serrana and North. This one did not present a defined peak in the Northwest region; it took one week to present its peak in Ilha Grande Bay and three weeks or more in Center-South and Médio Paraíba. For the Gama variant, the peak in the Ilha Grande Bay region occurred one week before it did in Metropolitan region I. There was synchronicity between five regions: Metropolitan I, Metropolitan II, Baixada Litorânea, Center-South and Médio Paraíba. In the Serrana region, it took one week to occur, whereas in the North region it took two weeks; in the Northwest region, it took three weeks or more.

The Delta variant only presented synchronicity between Metropolitan I and Serrana regions. The Center-South region had its peak after one week; Metropolitan II and Médio Paraíba, after two weeks; and Northwest, North, Baixada Litorânea and Ilha Grande Bay did not have a defined peak. The Omicron variant BA.1 was present in all regions; however, its peak took two weeks to occur in the Northwest and North regions; one week, in Metropolitan II, Baixada Litorânea, Serrana, Center-South and Médio Paraíba. The peak in Ilha Grande Bay was synchronous with Metropolitan I. The peak of Omicron BA.2 subvariant was synchronous in regions of Ilha Grande Bay, Northwest, Metropolitan I and II, but did not have a well-defined peak in the Serrana region, whereas it happened one week before in the Center-South and Baixada Litorânea. A late peak, of one week, was observed in Médio Paraíba and North regions. The Omicron BA.5 variant presented synchronicity between Center-South, Metropolitan I, Metropolitan II, Serrana and North regions. It took one week for the peak to occur in the Ilha Grande Bay, Médio Paraíba and Baixada Litorânea, and three weeks or more in the Northwest region.

DISCUSSION

The main objective of this study was to carry out a spatiotemporal synchronicity analysis of the COVID-19 epidemic in the health regions of the Rio de Janeiro state , including the period from March 15, 2020 (EW 9) to December 31, 2022 (EW 52). The results of this spatiotemporal analysis revealed an introduction and asynchronous increase of confirmed cases of COVID-19 in the different health regions of the Rio de Janeiro state , showing seven waves throughout the study period.

The results pointed out to a first peak related to the B.1.1.33 variant, with the highest incidence of hospitalizations concentrated in Metropolitan Region I, which is the area that contains the state capital, the city of Rio de Janeiro. The fact that the city of Rio de Janeiro is an important entry point for national and international travelers, besides being a globally known tourism destination, may have contributed with the early introduction of the virus in the territory, making it the second city in the country to register confirmed cases of COVID-19¹³13. Candido DDS, Watts A, Abade L, Kraemer MUG, Pybus OG, Croda J, et al. Routes for COVID-19 importation in Brazil. J Travel Med 2020; 27(3): taaa042. https://doi.org/10.1093/jtm/taaa042
https://doi.org/10.1093/jtm/taaa042... .

The synchronicity analysis showed that, after the B.1.1.33 variant, there was a significant increase in cases not only in Metropolitan region I, but also in all other health regions. Initially, the propagation of the virus was influenced by the existing socioeconomic inequalities between the regions, and the lack of a coordinated, efficient and equalitarian response favored the asynchronous and generalized dissemination of the virus¹⁴14. Ahmed F, Ahmed N, Pissarides C, Stiglitz J. Why inequality could spread COVID-19. Lancet Public Health 2020; 5(5): e240. https://doi.org/10.1016/S2468-2667(20)30085-2
https://doi.org/10.1016/S2468-2667(20)30... . The vulnerable population, who depends on public transportation and on the public health system, had more difficulties to adjust to social distancing measures due to poor household conditions and overcrowded vehicles, health units and hospitals¹⁵15. Kopel J, Perisetti A, Roghani A, Aziz M, Gajendran M, Goyal H. Racial and gender-based differences in COVID-19. Front Public Health 2020; 8: 418. https://doi.org/10.3389/fpubh.2020.00418
https://doi.org/10.3389/fpubh.2020.00418... ^,1616. Souza CDF, Paiva JPS, Leal TC, Silva LF, Santos LG. Evolução espaçotemporal da letalidade por COVID-19 no Brasil, 2020. J Bras Pneumol 2020; 46(4): e20200208. https://dx.doi.org/10.36416/1806-3756/e20200208
https://doi.org/10.36416/1806-3756/e2020... .

The dissemination of the subsequent variant, Zeta, showed a pattern of higher incidence of hospitalizations in the following regions: Serrana, Center-South, Metropolitan II and Ilha Grande Bay, indicating an expansion of the virus to the countryside of the state, with higher peaks than in Metropolitan Region I. Besides, there was a more frequent movement of the virus lineage within the state in comparison with the movement between states, reinforcing the importance of control measures to prevent outbreaks¹³13. Candido DDS, Watts A, Abade L, Kraemer MUG, Pybus OG, Croda J, et al. Routes for COVID-19 importation in Brazil. J Travel Med 2020; 27(3): taaa042. https://doi.org/10.1093/jtm/taaa042
https://doi.org/10.1093/jtm/taaa042... . In that moment, control measures for the circulation of the virus had been reduced, which coincided with a second wave that began in November, 2020, in the electoral period of the country¹⁷17. Silva SJR, Pena L. Collapse of the public health system and the emergence of new variants during the second wave of the COVID-19 pandemic in Brazil. One Health 2021; 13: 100287. https://doi.org/10.1016/j.onehlt.2021.100287
https://doi.org/10.1016/j.onehlt.2021.10... . The cities softened restrictive measures in different moments, which facilitated the circulation of the virus in an asynchronous manner due to the mobility of the population, acting as a trigger for disease propagation¹⁸18. Kortessis N, Simon MW, Barfield M, Glass GE, Singer BH, Holt RD. The interplay of movement and spatiotemporal variation in transmission degrades pandemic control. Proc Natl Acad Sci U S A 2020; 117(48): 30104-6. https://doi.org/10.1073/pnas.2018286117
https://doi.org/10.1073/pnas.2018286117... .

The introduction of the vaccination campaign against COVID-19, in January 2021, with the Gama variant in the state, played a significant role in the reduction of hospitalizations, but the highest incidence rates of hospitalization remained in the countryside regions, due to the non-homogeneous distribution of vaccines among cities. However, the incidence of hospitalizations after this variant decreased due to the ongoing vaccination campaign, which reduced the number of hospitalized patients, relieving the load on health systems and allowing the flexibilization of social distancing measures¹⁹19. Kissler SM, Tedijanto C, Goldstein E, Grad YH, Lipsitch M. Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. Science 2020; 368(6493): 860-8. https://doi.org/10.1126/science.abb5793
https://doi.org/10.1126/science.abb5793... ^,2020. Sharpe HR, Gilbride C, Allen E, Belij-Rammerstorfer S, Bissett C, Ewer K, et al. The early landscape of coronavirus disease 2019 vaccine development in the UK and rest of the world. Immunology 2020; 160(3): 223-32. https://doi.org/10.1111/imm.13222
https://doi.org/10.1111/imm.13222... .

The fourth identified variant, Delta, presented lower hospitalization rates than the previous variants, since vaccination had begun in the previous months. The CoronaVac vaccine, the first one to be available in Brazil, especially in the first days of vaccination, proved to be relatively efficient against the most severe cases of this variant and its predecessors, reducing cases of hospitalization²¹21. Melo-González F, Soto JA, González LA, Fernández J, Duarte LF, Schultz BM, et al. Recognition of variants of concern by antibodies and T cells induced by a SARS-CoV-2 inactivated vaccine. Front Immunol 2021; 12: 747830. https://doi.org/10.3389/fimmu.2021.747830
https://doi.org/10.3389/fimmu.2021.74783... .

With the circulation of Omicron subvariants (BA.1, BA.2 e BA.5), there was an increase in cases in all health regions, even though the percentage of people who needed specialized care, such as hospitalization, was lower than in comparison to other variants. This variant presents several mutations with potential to increase transmissibility, providing resistance to therapy or partially escaping the immunity induced by infection or vaccine²²22. Schmidt F, Muecksch F, Weisblum Y, Silva J, Bednarski E, Cho A, et al. Plasma neutralization of the SARS-CoV-2 omicron variant. N Engl J Med 2022; 386(6): 599-601. https://doi.org/10.1056/NEJMc2119641
https://doi.org/10.1056/NEJMc2119641... . Mass vaccination stood out as an essential measure to prevent clinical complications caused by COVID-19. It is important to emphasize that this disease can be prevented from the point of view of clinical complications with high vaccine coverage²³23. Fundação Oswaldo Cruz. O avanço da variante Ômicron, a resposta das vacinas e o risco de desassistência [Internet]. Nota Técnica 24 de 10 de fevereiro de 2022 [acessado em 03 jul. 2023]. Disponível em: Disponível em: https://www.arca.fiocruz.br/bitstream/handle/icict/51252/nota_tecnica_24.pdf?sequence=2&isAllowed=y
https://www.arca.fiocruz.br/bitstream/ha... .

This study has limitations, even though hospitalizations are a fraction of cases, they are influenced by the vaccines. So, the use of hospitalizations to describe the dynamics of the epidemic in cases of COVID-19 in the Rio de Janeiro state is a study limitation. However, the data about hospitalizations are more robust to describe the dynamics of the COVID-19 epidemic than the data regarding cases, which face several biases, such as lack of an organized data system, unavailability of tests, lack of notification in self-test results, delay in notifications and quality of the filled out information.

The search for equality in vaccination in all cities of the state is recommended, aiming at the fast dissemination of the virus, preventing severe cases, preventing the appearance of new variants and reducing the demand for specialized health care. The knowledge of prevalent variables in each region allows for better planning of vaccine distribution and strategies of immunization, in order to increase vaccine coverage and reduce the dissemination of the virus. An additional and very relevant aspect is found in the supervision of the possible new ability of new variants to escape the immune response offered by the vaccines.

The results of this study have significant implications for the continuous monitoring of COVID-19 and the public health surveillance measures. When understanding the spatial and temporal dynamics of the epidemic, health administrators may address the available resources in a more efficient way, besides implementing preventive measures and punctual interventions in the most affected areas. Through this context, when implementing such monitoring in regions where the variant was previously established, and by identifying any signs of this occurrence, the warnings addressed to these areas where virus circulation is still less disseminated can be anticipated. The inclusion of these elements, when conducted in a concrete manner for surveillance, has the potential of considerably increasing a prompt response when facing this situation. This analysis can be applied in other epidemiological scenarios, allowing the monitoring of diseases with epidemic potential in different places, becoming an important tool to improve public health response, anticipating outbreaks and making decisions that are more based on the control of infectious diseases, thus contributing with the promotion of health and well-being of the population.

REFERENCES

Publication Dates

History