key: cord-0966603-tu2o2k4s authors: Bermudi, Patricia Marques Moralejo; Lorenz, Camila; Aguiar, Breno Souza de; Failla, Marcelo Antunes; Barrozo, Ligia Vizeu; Neto, Francisco Chiaravalloti title: Spatiotemporal ecological study of COVID-19 mortality in the city of São Paulo, Brazil: Shifting of the high mortality risk from areas with the best to those with the worst socio-economic conditions date: 2020-12-02 journal: Travel Med Infect Dis DOI: 10.1016/j.tmaid.2020.101945 sha: 396c045611f2156762d1bdd8e7aa952d79409505 doc_id: 966603 cord_uid: tu2o2k4s BACKGROUND: Currently, Brazil is experiencing one of the fastest increasing coronavirus disease (COVID-19) mortality rates worldwide, with a minimum of 158,000 confirmed deaths presently. The city of São Paulo is particularly vulnerable because it is the most populated city in Brazil. Thus, this study aimed to analyse COVID-19 mortality in a spatiotemporal context in São Paulo, with respect to socio-economic levels. METHOD: We modelled the deaths using spatiotemporal architectures and Poisson probability distributions using a latent Gaussian Bayesian model approach. RESULTS: Both total deaths and confirmed deaths showed similar spatial patterns. Mortality was higher in men and increased with age. The most critical period regarding mortality occurred between the 20(th) and 23(rd) epidemiological weeks, followed by an apparent stabilisation of the epidemiological trend. The risk of death was greater in areas with the worst social conditions during the study period. However, this pattern was not uniform over time, since we identified a shift of high risk from the areas with the best socio-economic conditions to those with the worst conditions. CONCLUSIONS: Our study corroborated the relationship between COVID-19 mortality and socio-economic conditions, revealing the importance of geographic screening in the integration of better actions to face the pandemic. Mortality owing to coronavirus disease , caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is considered avoidable because some collective and individual measures can help prevent the infection, and appropriate health assistance greatly reduces the mortality risk. Nevertheless, mortality risk varies according to some individual and geographic risk factors, resulting in health inequity observed in several parts worldwide since the beginning of the pandemic [1, 2] . In Brazil, COVID-19 was first reported in the city of São Paulo on 25 February 2020 [3] . Until 3 August, 4 months after the first reported death, the disease had already resulted in 2,962,442 confirmed cases and 99,572 deaths in the country [4] . Currently, the United States of America and Brazil are the epicentres of the disease. Brazil is the fifth country worldwide in terms of surface area and population and is classified as having an upper-middle-income economy [5] . Almost 60% of the Brazilian population is concentrated in 6% of the large cities, among which São Paulo is the largest. The Brazilian Unified Health System guarantees healthcare for all citizens as well as for thousands of foreigners residing or passing through the country [6] . Despite this universal healthcare, geographic differences in mortality rates across different areas have been observed on a national and intra-urban scale. Spatial heterogeneity in population characteristics such as age, underlying health, household densities, partial lack of sanitation, socio-economic status, contact networks, and mobility patterns [7] has emerged as a potential propellant of the spatiotemporal spread of the disease. Since the beginning of the pandemic, mapping disease occurrence and spread has become a powerful tool to track the disease and establish measures to slow down the transmission of the infection both locally and globally [8] . Web-based Geographical Information Systems have allowed near real-time monitoring using map-centric dashboards [9] [10] [11] . Despite advances in the use of technology to reduce the impact of the pandemic, little is known about the spatiotemporal patterns of COVID-19 mortality, especially in intra-urban settings. Studying the spatiotemporal dynamics of mortality instead of the cases may help better evaluate inequity. All health disadvantages accumulated over decades of life, owing to any kind of deprivation, increase the risk of mortality due to COVID-19. The lack of a robust spatiotemporal analysis undermines the comprehension of mitigation strategies to potentialise disease-control efforts. Thus, this study aimed to unveil the spatiotemporal dynamics of COVID-19 mortality in a higher spatial resolution (in the city of São Paulo), considering the socio-economic context of the J o u r n a l P r e -p r o o f population. This approach may shed light on the urgent need for solid evidence on health inequities during the COVID-19 outbreak. This ecological study, based on COVID-19 secondary mortality data, was delineated in the city of São Paulo, state of São Paulo, Brazil ( Figure 1 ). In 2020, the estimated population of this city was 11,869,660 inhabitants, and the mean demographic density was 7,803 inhabitants/km² [12] . When using COVID-19 data, underreporting is always an issue that deserves attention, even for mortality. To minimise the effect of possible sub-notification of deaths, we analysed confirmed and suspected deaths due to COVID-19. Thus, mortality data Suspected deaths were coded as U04.9 (corresponding to severe acute respiratory syndrome). The places of residence of the COVID-19 deaths were geocoded with CEInfo/SMS-SP using its databases and Google Maps API geocoding script (that uses public places as the base map). The resulting geocoded addresses were validated by comparing the road or Zone Improvement Plan (ZIP code), whenever the record was allocated using the original ZIP code. Geocoded data were assigned to the 310 sample areas of the Brazilian Institute of Geography and Statistics, for which demographic and socio-economic census data are available [13] . We considered these areas as the spatial units in our models ( Figure 1 ). the e-SIC database. It was not necessary to submit this study to an ethics committee because we did not have access to personal data, such as names and addresses. The use of secondary data without personal identification and in a public domain dispenses the need for prior approval from the Ethics Committee on Research with Human Beings (as per Resolution No. 510/2016 of the National Health Council) [16] . We also used the data available in the TabNet and named this information as the "Tabnet database". The Tabnet database is an application, available at , provided by the Municipal Health Department of São Paulo. It was developed by DataSUS. This application provides free access (to any user) to population databases and to database information systems of SUS (such as the SIM), which is supplied by the Secretariat's Program (PRO-AIM). Through the Tabnet database application, it is possible to perform tabulations and to cross several variables of interest such as EW, sex, age group, and specific cause. The databases are updated periodically. Notably, the data sources of the Tabnet and e-SIC database are SIM. However, only the e-SIC database is subdivided according to sample area. To measure the socio-economic condition of the population, since individual-level data are not available in the mortality database, we used a socio-economic index elaborated for health research. The Socio-economic Index of the Geographic Context for Health Studies (GeoSES) [17] was developed using principal component analysis, starting with 41 variables. The index conceives the socio-economic condition by considering seven parameters based on the theoretical background [18, 19] : education, mobility, poverty, wealth, income, segregation, and deprivation of resources and services. The index was defined on three scales: national, Federative Unit, and intra-municipal. Figure 1 presents GeoSES for the sample areas of the city of São Paulo. It shows that the areas with the best socio-economic conditions (GeoSES equal to or close to 1) are located in the central part of the city and that the socio-economic conditions deteriorate towards the periphery, where they reach the worst levels (GeoSES equal to or close to -1). It has been shown to be useful in studies of mortality due to avoidable causes of deaths in individuals aged 5-74 years due to interventions at the Brazilian health system on a national scale as well as mortality due to circulatory system diseases in the city of São Paulo [17] . This index is publicly available at https://opendatasus.saude.gov.br/dataset/geoses. We used the information of confirmed (B34.2) and suspected (U04.9) deaths due to COVID-19 available for the entire city of São Paulo from the 11 th to 29 th EWs to calculate the weekly mortality rates of confirmed (B34.2), suspected (U04.9), and total (B34.2 + U04.9) COVID-19 deaths. This was similarly done using the e-SIC database from the 11 th to 12 th EW. We excluded (from the information e-SIC database) COVID-19 deaths that occurred in the 25 th EW because data from this week were incomplete (it was extracted on 18 June 2020, and it contained only part of the 25 th EW information). These rates were obtained by dividing the respective numbers of deaths in each week by the total population of the city and presented as death per 100,000 inhabitants per week. In this sequence, we obtained the mortality rates for confirmed, suspected, and total COVID-19 deaths according to sex and age for the entire period from the 11 th to 12 th EW using the Tabnet and e-SIC databases. The comparisons between the data from these two sources were useful to evaluate the completeness of the data we used for the spatial and spatiotemporal analysis. For the calculation of the mortality rates according to sex and age, we excluded the data without these information. As only one suspected COVID-19 death occurred during the 11 th EW, we restricted our spatial and spatiotemporal analysis from the 12 th to 24 th EW, and spatial or spatiotemporal architecture was considered in all models developed. We first modelled the confirmed and total COVID-19 deaths using spatiotemporal models involving only the intercept and random effects that account for spatial and temporal autocorrelation and the interaction between them. The spatial dependence was modelled considering the Besag-York-Mollié (BYM) model with two components representing the spatially structured and non-structured random effects [20, 21] . These two components were considered independent from one another and followed the parameterisation proposed previously [22] . The temporal dependence was modelled using a non-structured random effect and a structured random effect provided by a random walk autoregressive model of first order (RW1). The interaction between space and time was modelled considering spatial and temporal non-structured random effects [21] . The number of confirmed and total COVID-19 deaths per EW and per sample area was modelled using Poisson and zero-inflated Poisson probability distributions with a latent Gaussian Bayesian model approach. We considered the expected confirmed and total COVID-19 deaths for each EW and for each spatial unit as offsets in these models. The expected deaths were estimated with indirect standardisation, considering the age and sex structure of each sample area and the mortality rates for the entire study period and city. This enabled us to interpret the outcomes of our analysis as relative risks (RR) J o u r n a l P r e -p r o o f concerning the mortality rates for the entire study period and city. From these models, we obtained the temporal and spatiotemporal RR. Subsequently, we introduced the socioeconomic covariate (GeoSES) in these models and obtained the corresponding RR. Finally, we used a spatial approach to model the confirmed and total COVID-19 death per EW to evaluate the role of the socio-economic covariate in each EW. To this end, we considered spatial models with intercept, BYM spatial random effects, and GeoSES as a covariate. The expected COVID-19 deaths were obtained in a similar way for the spatiotemporal models. We considered the entire city's mortality rates for each EW, and as such, we could interpret the RR of the entire city's mortality rates for each EW. We performed our models in a Bayesian context using the integrated nested Laplace approximation (INLA) approach [23] . We selected the best models using the Deviance Information Criterion (DIC) so that the best-adjusted models were those with lower DIC values [21] . We used non-informative priors for the fixed effects and priors with penalised complexity for the precision parameters of the random effects [22] . We ran our models in the R environment [24] . We found 14,753 confirmed and suspected COVID-19 deaths in the Tabnet database (from EW 11 to 29) and 10,760 in the e-SIC database (from EW 11 to 25) in the city of São Paulo. We excluded 67 deaths from the e-SIC database because they were referent to the 25 th EW, which was not completed when the data were extracted. Accordingly, we were left with 10,693 deaths. Figure 2 shows the mortality rates for COVID-19 for both data sources considering the confirmed, suspected, and total deaths. The curves from the e-SIC database are similar to the curves from the Tabnet database. Both curves have the same data source with different extraction dates, and the differences among these curves in the 23 rd and 24 th EW are related to a delay in the notification of the COVID-19 deaths. Table 1 shows the numbers and mortality rates of confirmed, suspected, and total COVID-19 deaths in the different sex and age groups, obtained from the e-SIC and Tabnet databases (from EW 11 to 24). To build Table 1 , we excluded six deaths with ignored age and three deaths with ignored sex from the 10,693 deaths in the e-SIC database. We excluded six deaths with ignored age and four deaths with ignored sex from the 11,098 deaths in the Tabnet database. We can observe that the mortality was higher in men and that it increased with age, resulting in a mortality rate of 460.9 deaths per 100,000 J o u r n a l P r e -p r o o f inhabitants (in 14 weeks) for people aged ≥ 60 years. This pattern of increased mortality with increase in age was also observed after making a greater stratification among people aged 60 years or above, as shown in Supplementary Material 1. We present the results of the models with spatiotemporal architecture using only the intercept. Figure 3 shows the temporal RR from the 12 th to 24 th EW. Notable features include the maximum RR, which occurred during the 20 th EW for the total COVID-19 deaths and during the 23 rd EW for the confirmed ones during our study period. These results, even adjusted for the temporal autocorrelation, are similar to those presented in Figure 2 . Considering the data presented in Figures 2 and 3 , the apparent pattern of the temporal curves shows a tendency to stabilise. However, a new rise in mortality cannot be ruled out. Figures 4 and 5 show the posterior means of the spatiotemporal RR for the sample areas and EW, respectively, for confirmed and total COVID-19 deaths. Apart from the fact that the RR is greater for the total deaths than for the confirmed ones, the distribution of the RR is similar between them, and it follows the behaviour of the temporal RR. In the first two EWs, the sample areas presented lower values of RR that increased over time. However, this increase occurred with greater intensity in peripheral areas. We also used the spatiotemporal models, considering socio-economic variables (GeoSES). Table 2 shows the spatiotemporal RR and the 95% credibility intervals for GeoSES obtained for the models with confirmed and total COVID-19 deaths. In both models, it is noted that a high socio-economic level protected against the risk of mortality due to COVID-19 throughout the study period. Thus, an increase of one unit in the socioeconomic indicator represented a 25% reduction in the risk of mortality due to COVID-19 (for the model using only confirmed deaths) and a 33% reduction in the risk of dying for the model using total deaths of COVID-19. Moreover, the risk of mortality due to COVID-19 in the sample areas with the best socio-economic conditions (GeoSES close to 1), compared with that in areas with the worst conditions (GeoSES close to -1), was 50% lower for the model with confirmed deaths and 66% lower for the model with total deaths. Finally, we performed the spatial modelling of confirmed and total COVID-19 deaths in the different EWs separately, considering the socio-economic covariate. Figure 6 shows the RR and 95% credibility intervals for GeoSES, with respect to the confirmed and total COVID-19 deaths for each of the EW. We identified a shift in the pattern of the relationship between COVID-19 mortality and socio-economic status over time. The best socioeconomic level presented itself as a risk factor for COVID-19 deaths in the first two EWs in the city of São Paulo. From the 15 th EW, for total deaths, and from the 16 th EW, for confirmed and total deaths, the worst socio-economic condition became a risk factor. Even if some values were not significant, there was a continuous decrease in RR from the 12 th to 17 th EW, followed by its stabilisation. This is the first population-based study on the evolution of the spatiotemporal pattern of COVID-19 mortality in the intra-urban setting of the largest city of Brazil. Using two different datasets, analysing the confirmed and confirmed plus suspected deaths separately, we were able to evaluate how uncertainty would impact the association between the RR and the socio-economic condition. The robust EW models clearly showed when the high risk of death shifted from areas with the best to those with the worst socioeconomic conditions in the city. Our findings showed that the most critical period regarding mortality due to COVID-19 in the city of São Paulo occurred between the 20 th and 23 rd EW, followed by an apparent stabilisation of the temporal trend. However, it is not possible to predict a future scenario. As social distancing measures have been lifted in the city since the 25 th week [25] , this could increase the number of infected people and, consequently, the number of deaths. Although social distancing alone seems insufficient to contain COVID-19, many studies have frequently concluded that it remains a critical component of the outbreak control [26] . It is important to point out that both total deaths and confirmed deaths showed similar spatial patterns in our study, despite their differences. The suspected deaths, on the one hand, need to be treated with caution because they may not be COVID- 19-related and, on the other hand, could be considered as one of the strengths of the study. This is because there is a delay in the confirmation of suspected cases, and consequently, some of the suspected deaths are later confirmed to be COVID-19-related deaths. Furthermore, some of the suspected cases would be considered as confirmed COVID-19-related deaths if the case definition was not strict (as the required test is often not available or is not performed within the appropriate window time). From this perspective, the amplitudes of variation in rates and RR obtained from confirmed and total deaths could be considered as lower and upper limits (or vice versa) for the magnitudes of these measures. The elderly population represents one of the groups that are more prone to COVID-19 in the city of São Paulo. The risk of death was the highest among men aged > 70 years; this is similar to data from China and the United States [27] . Recently, Souza et al. [3] analysed the Brazilian population and found that most COVID-19 deaths were of men and that the most frequent comorbidities were cardiovascular disease and diabetes. Behavioural factors, especially social status, that may prejudice adherence to lockdown measures, have been shown to be potentially crucial in determining susceptibility to SARS-CoV-2 [28, 29] . This unequal death ratio in men may be interpreted considering many factors: the comparatively higher prevalence of comorbidities (i.e., hypertension, diabetes, cardiovascular disease, and chronic lung disease) [30] , higher risk behaviours J o u r n a l P r e -p r o o f (i.e., smoking and alcohol use), occupational exposure [31] , and sex differences in immune responses [32] . However, there may be other social and behavioural characteristics that favour women as reported in previous studies, which proposed that women are more likely than men to adopt hand hygiene practices [33] and to seek preventive care [34] . The spatial distribution of suspected and confirmed deaths due to COVID-19 in the city of São Paulo shows inequalities, with spatial dependence and positive correlation associated with socio-economic factors of the areas. This is remarkably similar to the results of Maciel et al. [35] . Our findings reveal that the socio-economic condition acts as a protective factor against the risk of mortality due to COVID-19. In the models considering only confirmed deaths and all deaths, an increase of one unit in the socio-economic indicator represented a 25% and 33% decrease in the risk of mortality, respectively. The first observation is that, when considering all deaths, the protective effect of the socioeconomic level is more evident, showing that there must be a higher incidence of suspected deaths in the less favoured areas than in the most favoured areas (in areas with better socio-economic level, confirmation tests for COVID-19 are probably more available). A study conducted by Souza et al. [3] reinforced this finding. They compared the spatial pattern of confirmed cases of COVID-19 and severe acute respiratory infection with unknown aetiology, with per capita income in the metropolitan region of São Paulo. They found that the COVID-19 cases were more associated with better levels than the latter. They pointed out that the degree of underreporting of COVID-19 cases would increase with a decrease in socio-economic status. Therefore, our results confirmed the association between COVID-19 and human development, pointing out the importance of geographic screening in locations with a potential for transmission of local infectious diseases, as a fundamental aspect to coordinate better actions against the pandemic [35] . The low socio-economic condition levels demonstrate not only the vulnerability of the population but also the difficulties health services face with respect to the diagnosis and treatment of the COVID-19. This is similar to the overall fragility expected from health services in Brazil [36] and Latin American countries facing the pandemic [37] . Living conditions may also be strongly influenced by low income in different ways, such as residence in poorer neighbourhoods and poorer housing conditions (particularly confined or overcrowded housing) [38] . The latter condition has been demonstrated to be associated with a greater risk of contagion from several other pathogens such as Helicobacter pylori [39] , tuberculosis [40] , or Epstein-Barr virus [41] . Regarding COVID- 19 , studies have shown that occupation is an explicit determinant of contagion and a J o u r n a l P r e -p r o o f secondary determinant of COVID-19 severity and deaths due to the association between occupational social class and comorbidities [38] . For example, workers such as cleaners, retail staff, teachers, or healthcare workers have a direct impact on COVID-19 incidence [42] . Individuals with underprivileged socio-economic conditions are more prone to be exposed to job stress, including burnout syndrome and unemployment, which may contribute to disrupted immune and inflammatory system responses [43, 44] as well as a higher risk of comorbidities with COVID-19 [45] . Until now, both debilitated immunity and the existence of comorbidities are recognised risk factors for COVID-19 severity [38] . We showed that the first cases of deaths occurred in the neighbourhoods with the best socio-economic conditions in the city of São Paulo. This may be because all the infected subjects had been abroad [46] . In the first two weeks, the best socio-economic conditions was presented as a risk factor. Then, there was a change in the spatial pattern from the fourth week onwards, when the worst socio-economic level became a risk factor. Similarly, Souza et al. [3] showed a higher risk of diagnosed COVID-19 cases in census tracts with higher per capita income in the São Paulo metropolitan region during the early phase of the COVID-19 epidemic. After these first cases in the richest areas, the virus started to circulate in the suburbs of the city, with high population density and worsened sanitary conditions [47] . This probably explains its fast transmission. The city of São Paulo is particularly vulnerable because it is the most populated city in the country, with approximately 12 million inhabitants [48] , and it is highly connected with other countries worldwide via its main airport, the São Paulo-Guarulhos International Airport, which is the largest in Brazil, with non-stop passenger flights to 103 destinations across 30 countries [37] . In this study, we used SIM instead of the SIVEP Gripe database, unlike other studies [3] . The recommendation of the State Health Department [49] to register the notification of death and the monitoring of mortality using SIM, in practice, leads to a timelapse between the event and the use of information. COVID-19 mortality data were improved in the city of São Paulo with PRO-AIM, using the Laboratory Environment Manager System (GAL) and Flu Surveillance Information System (Sivep Gripe). These characteristics, combined with the possibility of assessing home or unattended deaths, motivated the option of using SIM data (which considered the confirmed and suspected diagnoses in analysing mortality due to . It is necessary to emphasise that the recommendation to register confirmed and suspected deaths due to COVID-19 [49] (that motivated the change of the codes assigned by the health services since the beginning of J o u r n a l P r e -p r o o f the crisis, and which is applied by SMS-SP to monitor the evolution of events) has not been applied yet. Our study findings must be considered in the context of several assumptions and data limitations. We associated patients' addresses or postcodes to the area-based socioeconomic conditions using geolocalisation. This may provide some insight into the likelihood of exposure to health factors and COVID-19 risks. This approach is frequently used as representative of individual socio-economic conditions. Nevertheless, they are not a perfect picture of the individuals' conditions, and they could underestimate the magnitude of social disproportion related to individual social measures [50] . Rather, they are best employed along with individual-level variables to reflect geographical or aggregate-level risks [38] . We highlight that our spatial analysis is subject to methodological limitations caused by ecological fallacy and the modifiable areal unit problem. These constraints are intrinsic to any spatial analysis that uses aggregated data [51] . Despite these, our study still contributes to healthcare planning measures and to future precision studies focusing on the effects of social health factors on COVID-19 deaths. In addition, one of the strengths of our study was that it dealt with COVID-19 deaths instead of the cases (owing to the better accuracy and coverage of the data). When we consider only the cases, many asymptomatic cases may not be reported and this could hamper the conclusions. We used models with spatial and spatiotemporal architectures to investigate the patterns of confirmed and total (confirmed and suspected) COVID-19 deaths in the city of São Paulo. The obtained results, after considering both categories, showed differences regarding the magnitude of the rates and RR. However, there were no differences with respect to the derived conclusions. The risk of mortality due to COVID-19 was the highest between the 20 th and 23 rd EW, followed by an apparent stabilisation of the temporal trend. However, we did not rule out a possible future rise in mortality. A high socio-economic level was shown to protect against the risk of mortality due to COVID- 19 Posterior means of the relative risks and 95% credible interval for the socioeconomic covariate obtained with spatial models for confirmed and total COVID-19 deaths, according to each one of the epidemiologic weeks. City of São Paulo, 12 th to 24 th Epidemiology Week, 2020. J o u r n a l P r e -p r o o f COVID-19 and Inequity: a Comparative Spatial Analysis of New York City and Chicago Hot Spots Combating COVID-19: health equity matters Epidemiological and clinical characteristics of the COVID-19 epidemic in Brazil The World Bank DRG. World Development Indicators SUS 30 anos: o início, a caminhada e o rumo Mathematical model describing CoViD-19 in São Paulo, Brazil -evaluating isolation as control mechanism and forecasting epidemiological scenarios of release Geographical tracking and mapping of coronavirus disease COVID-19/severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) epidemic and associated events around the world: how 21st century GIS technologies are supporting the global fight against outbr Coronavirus COVID-19 Global Cases (dashboard) Novel coronavirus (COVID-19) situation (dashboard) 2020 Metodologia do Censo Demográfico SINAN -Calendário Epidemiológico Pedidos de informação protocolados à Prefeitura via sistema e-SIC 2020 Resolução n o 510 de 7 de abril de 2016 que trata das especificidades éticas das pesquisas nas ciências humanas e sociais e de outras que utilizam metodologias próprias dessas áreas GeoSES: A socioeconomic index for health and social research in Brazil Optimal Indicators of Socioeconomic Status for Health Research Measuring social class in US public health research: concepts, methodologies, and guidelines Spatial and spatio-temporal models with R-INLA Spatial and Spatio-temporal Bayesian Models with R-INLA Penalising Model Component Complexity: A Principled, Practical Approach to Constructing Priors Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations R: A Language and Environment for Statistical Computing Isolation Measures in the State of Sao Paulo Quarantine alone or in combination with other public health measures to control COVID-19: a rapid review The Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. Vital Surveillances: The Epidemiological Characteristics of an Outbreak of 2019 Novel Coronavirus Diseases (COVID-19) 2020 Sex Differences in Everyday Risk-Taking Behavior in Humans Non-white ethnicity, male sex, and higher body mass index, but not medications acting on the renin-angiotensin system are associated with coronavirus disease 2019 (COVID-19) hospitalisation: review of the first 669 cases from the UK biobank Sex Differences in Mortality From COVID-19 Pandemic COVID-19 Sex-Disaggregated Data Tracker 2020 Sex differences in immune responses Sex Differences in Public Restroom Handwashing Behavior Associated with Visual Behavior Prompts Gender differences in the utilization of health care services -PubMed Análise inicial da correlação espacial entre a incidência de COVID-19 e o desenvolvimento humano nos municípios do estado do Ceará no Brasil COVID-19) and healthcare-associated infections: Emerging and future challenges for public health in Brazil COVID-19 in Latin America: The implications of the first confirmed case in Brazil Importance of the COVID-19 outbreak onwards Relation between infection with Helicobacter pylori and living conditions in childhood: evidence for person to person transmission in early life Role of socio-economic factors in tuberculosis prevalence The role of the early social environment on Epstein Barr virus infection: A prospective observational design using the Millennium Cohort Study Occupational risks for COVID-19 infection Multi-cohort study identifies social determinants of systemic inflammation over the life course Work Stress as a Risk Factor for Cardiovascular Disease Pandemia do Coronavírus no Brasil: Impactos no Território Cearense Orientações para codificação das causas de morte no contexto da COVID-19 Using ecological socioeconomic position (SEP) measures to deal with sample bias introduced by incomplete individual-level measures: inequalities in breast cancer stage at diagnosis as an example Revisiting Robinson: The perils of individualistic and ecologic fallacy Supplementary Material 1. Number and mortality rates (per 100,000 inhabitants in fourteen-weeks) of suspected (severe acute respiratory syndrome -U04.9), confirmed (coronavirus infection disease -B34.2) and total (U04.9 + B34.2) COVID-19 deaths, according to e-SIC (São Paulo Electronic Information System) and Tabnet databases and the age groups among people over 60 years of age. City of São Paulo, 11 th to 24 th epidemiological weeks, 2020. Data source: Deaths: Mortality Information System of Mortality Information Improvement Program, from the Epidemiology and Information Coordination of the São Paulo Health Secretariat (SIM/PRO-AIM/CEInfo/SMS-SP) Population: Fundação Sistema Estadual de Análise de Dados Estatísticos (SEADE). Tabnet database was update on 7/23/2020 and e-SIC database was provided on 6/08/2020.