Article Text

Original research
Characteristics, outcomes and risk factors for mortality of 522 167 patients hospitalised with COVID-19 in Brazil: a retrospective cohort study
  1. Marcia C Castro1,
  2. Susie Gurzenda1,
  3. Eduardo Marques Macário2,
  4. Giovanny Vinícius A França2
  1. 1Global Health and Population, Harvard TH Chan School of Public Health, Boston, Massachusetts, USA
  2. 2Secretariat of Health Surveillance, Brazilian Ministry of Health, Brasilia, Brazil
  1. Correspondence to Marcia C Castro; mcastro{at}hsph.harvard.edu

Abstract

Objective To provide a comprehensive description of demographic, clinical and radiographic characteristics; treatment and case outcomes; and risk factors associated with in-hospital death of patients hospitalised with COVID-19 in Brazil.

Design Retrospective cohort study of hospitalised patients diagnosed with COVID-19.

Setting Data from all hospitals across Brazil.

Participants 522 167 hospitalised patients in Brazil by 14 December 2020 with severe acute respiratory illness, and a confirmed diagnosis for COVID-19.

Primary and secondary outcome measures Prevalence of symptoms and comorbidities was compared by clinical outcomes and intensive care unit (ICU) admission status. Survival was assessed using Kaplan Meier survival estimates. Risk factors associated with in-hospital death were evaluated with multivariable Cox proportional hazards regression.

Results Of the 522 167 patients included in this study, 56.7% were discharged, 0.002% died of other causes, 30.7% died of causes associated with COVID-19 and 10.2% remained hospitalised. The median age of patients was 61 years (IQR, 47–73), and of non-survivors 71 years (IQR, 60–80); 292 570 patients (56.0%) were men. At least one comorbidity was present in 64.5% of patients and in 76.8% of non-survivors. From illness onset, the median times to hospital and ICU admission were 6 days (IQR, 3–9) and 7 days (IQR, 3–10), respectively; 15 days (IQR, 9–24) to death and 15 days (IQR, 11–20) to hospital discharge. Risk factors for in-hospital death included old age, Black/Brown ethnoracial self-classification, ICU admission, being male, living in the North and Northeast regions and various comorbidities. Age had the highest HRs of 5.51 (95% CI: 4.91 to 6.18) for patients≥80, compared with those ≤20.

Conclusions Characteristics of patients and risk factors for in-hospital mortality highlight inequities of COVID-19 outcomes in Brazil. As the pandemic continues to unfold, targeted policies that address those inequities are needed to mitigate the unequal burden of COVID-19.

  • public health
  • COVID-19
  • public health

Data availability statement

Data are available in a public, open access repository. The SIVEP-Gripe data set is publicly available on the Ministry of Health’s DATASUS website (https://opendatasus.saude.gov.br/nl/dataset).

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Strengths and limitations of this study

  • The strength of this study is that it leverages Brazil’s established national Influenza Epidemiological Surveillance Information System data to present comprehensive characteristics, clinical course and risk factors for COVID-19 in-hospital deaths across Brazil.

  • Administrative records lack details available in hospital medical records and may have accuracy and completeness problems.

  • We did not have access to laboratory results other than COVID-19 tests (eg, complete blood count) that would allow for a better characterisation of the clinical course of the disease.

  • COVID-19 deaths at home likely follow a different clinical course than deaths in the hospital and are not included in this analysis.

  • In-hospital deaths due to COVID-19 are likely under-reported and are limited by the testing protocol and capacity of each hospital.

Introduction

On 11 March 2020, the WHO declared COVID-19 as a pandemic. Caused by the novel coronavirus SARS-CoV-2, it emerged in China and quickly spread across the country and beyond. As of 16 February 2021, it was present in 223 countries and territories, with 108 822 960 confirmed cases and 2 403 641 confirmed deaths.1 Brazil recorded the first confirmed COVID-19 case on 26 February 2020 and the first death on 12 March, both in São Paulo State. In 24 days, the disease had spread to all Federal Units. As of 16 February 2021, 9 834 513 cases (9% of worldwide cases) and 239 245 deaths (10% of worldwide deaths) had been reported in Brazil, the second-highest in the world, behind only the USA. These numbers are underestimated since most mild cases are not being tested and thus are not likely to be reported, and some deaths may be reported with ill-defined causes, or not reported at all.

Brazil has a comprehensive health information system,2 with the systematic collection of births, deaths, hospitalisations and diseases of mandatory notification, among others. However, a complete and linked registry of records combining data from ambulatory and inpatient care, laboratory and radiologic results and outcome of the disease is not available. Therefore, there is limited information on the course of the disease for every case reported in Brazil.

Currently, the most detailed data available in Brazil refer to hospitalisations due to severe acute respiratory illness (SARI). Here, we use these data to provide a comprehensive description of demographic, clinical and radiographic characteristics, treatment, case outcome and risk factors associated with in-hospital death of patients hospitalised with SARI with a confirmed diagnosis for COVID-19, as of 14 December 2020. We analyse the largest retrospective number of cases (n=522 167) and we assess whether the Brazilian case is comparable to patterns previously described for other countries.

Methods

Data sources

We used deidentified records from the Influenza Epidemiological Surveillance Information System (Sistema de Informação de Vigilância Epidemiológica da Gripe, SIVEP-Gripe, in Portuguese), an information system of the Ministry of Health that captures all notifications of SARI hospitalisations in both public and private hospitals. The system is updated daily, and every 2 weeks, a new data set is made publicly available (https://opendatasus.saude.gov.br/nl/dataset). Here, we analysed records as of 14 December 2020 (n=1 029 684 notifications), after 15 419 duplicate records were removed by the Ministry of Health. Each record has data on patient’s age, sex, place of residence and of hospitalisation, ethnoracial self-classification,3 pregnancy status, comorbidities and symptoms; drug treatment; radiologic test results; and dates of illness onset, hospitalisation, ICU admission and outcome (death, release, still hospitalised). We considered only records of patients hospitalised with a confirmed diagnosis for COVID-19 (n=522 167). Diagnosis followed the Ministry of Health guidelines.4

Statistical analysis

Characteristics of inpatients were summarised in three groups: demographic, clinical and radiographic and treatment and outcomes. Medians and IQR ranges were used to describe continuous variables, and counts and percentages to describe categorical variables. Differences between inpatients that needed and did not need ICU admission and those that survived and did not survive were assessed by Whitney U, χ2 or Fisher’s exact test, as appropriate. No data imputation was performed for missing data (see online supplemental table 1 for information on data completeness).

Survival curves of inpatients at 60 days of hospitalisation by age, sex, ethnoracial self-classification, region and ICU admission were estimated using the Kaplan-Meier estimator and compared with the log-rank test. Factors associated with inpatient death were identified by univariable and multivariable logistic regression (excluding from the analysis those that remained hospitalised). Considering time to death as the outcome, HRs were estimated using Cox proportional-hazards models. Based on previous studies5–7 and on our available information, covariates included in both logistic and Cox models were age (0–19, 20–39, 40–59, 60–69, 70–79 and 80 or more years), sex, ethnoracial self-classification (White, Black/Brown, other, not reported), region (North—where Amazonia is located, Northeast, South, Southeast—where the cities of São Paulo and Rio de Janeiro are located and Center-West), comorbidities (diabetes, asthma, chronic liver disease, chronic neurological disease, chronic lung disease, immunodeficiency and chronic kidney disease), obesity and ICU admission. The variable ethnoracial self-classification was missing in 23.1% of the records, and we added those as a separate category (not reported). Distances between municipalities of residence and hospitalisation were calculated in ArcMap, V.10.6 (ESRI, Redlands, CA, USA). All analyses were performed in Stata, V.15.1 (Stata Corp., College Station, TX, USA), and R V.4.0.0 (RStudio Team, Boston, MA, USA).

Patient and public involvement

Our analysis used administrative records, and thus study participants were not involved in the design of the study. Public involvement was achieved through collaboration with the Ministry of Health, with whom we defined the research questions to fill in knowledge gaps and inform decision-making. Results were discussed and shared with the Ministry, and their wide dissemination with public health officials, researchers and through the media will reach the broader public.

Results

As of 14 December 2020, 522 167 patients had been hospitalised with confirmed COVID-19 since the beginning of the epidemic in Brazil. Of those, 296 002 (56.7%) were discharged, 1004 (0.002%) died of other causes, 160 495 (30.7%) died of causes associated with COVID-19, 53 503 (10.2%) remained hospitalised. Clinical outcome was unknown for 11 126 (2.1%) patients (table 1). The cumulative curve of hospital admissions (online supplemental figure 1) shows the fast increase in severe cases that required hospitalisation, following the steep increase in COVID-19 transmission in Brazil since the end of March. The median age of patients was 61 years (IQR, 47–73), and much higher for non-survivors, 71 years (IQR, 60–80), as shown by the age distribution in figure 1A,B. Patients aged 60 years or more represented 50.1% of hospitalisations, 59.0% of ICU admissions and 74.0% of deaths associated with COVID-19. Patients were mostly males (56.0%) and from the Southeast region (49.3%). Among females, 2.5% were pregnant or puerperal at the time of hospitalisation, and 7.5% of those died in the hospital. A total of 172 473 (33.0%) patients with median age of 65 years (IQR, 52–76) needed ICU admission. Of all hospitalisations, 37.7% of the patients were White, and 37.9% were Black/Brown. Among survivors, 38.8% were White, while among non-survivors 41.7% were Black/Brown. About 25% of the patients travelled a median of 32.0 km (IQR, 18.6–64.1) to be hospitalised in a municipality different from where they reside (table 1).

Table 1

Demographic characteristics of patients

Figure 1

Age distribution of patients, density curves of length of time from hospital admission to death and survival curve 60 days after hospital admission. (A) Age distribution of patients hospitalised. (B) Age distribution of in-hospital deaths. (C) Density of number of days from hospital admission to death up to 60 days after hospital admission, detailed by ICU admission status. (D) Survival curve estimated with Kaplan Meier and considering 60 days from hospital admission. ICU, intensive care unit.

Comorbidities were observed in 64.5% of the patients, 74.6% of those who needed ICU admission, 76.8% of non-survivors and 58.5% of the survivors and those whose death was not associated with COVID-19. With the exception of asthma, all comorbidities had a higher prevalence among non-survivors (compared with all patients). The most common comorbidities were chronic cardiovascular disease (34.5% of patients and 43.5% of non-survivors) and diabetes (25.7% of patients and 33.0% of non-survivors). Obesity was reported in 7.4% of the patients and 10.5% of those who needed ICU admission. The most common symptoms were fever, cough, shortness of breath, low oxygen saturation and respiratory distress symptoms (table 2).

Table 2

Clinical and radiographic characteristics of patients

The median time from illness onset to hospital admission was 6 days (IQR, 3–9), slightly shorter among non-survivors, 5 days (IQR, 2–8). Mechanical ventilation was needed by 62.2% of all patients, and by 75.6% of those who died. Invasive ventilation was more common in the ICU (44.0%). Oseltamivir, an antiviral medication, was the most common drug used during treatment (15.8% overall, and 17.6% among those in ICU), and the median time from illness onset to treatment was 5 days (IQR, 3–8). Of the patients who needed ICU, 51.8% died from causes associated with COVID-19, and 19.0% remained hospitalised after ICU discharge for 5 days (IQR, 2–10). The median time from illness onset to ICU admission was 7 days (IQR, 3–10), and the medium length of ICU stay was 8 days (IQR, 3–15) for all patients, 9 days (IQR, 4–16) for the deceased. Among the 160 495 patients who died of causes associated with COVID-19 by 14 December, the median time from illness onset to death was 15 days (IQR, 9–24) (table 3). Medium length of hospital stay was 8 days (IQR, 4–17), but longer for those who needed ICU admission, 12 days (IQR, 6–22). The density of time from hospital admission to death is positively skewed, more so for those who did not get admitted to the ICU (figure 1C).

Table 3

Treatment and outcomes of patients

Kaplan Meier curves (figure 1D and online supplemental figure 2) for a period of up to 60 days after hospital admission showed that survival curves were significantly different by age, region, sex, ethnoracial self-classification and ICU admission.

Univariable logistic analysis indicated that the odds of in-hospital death progressively increased with age, and were higher for patients who were male, non-white, from the North and Northeast regions, needed ICU care, were obese and had diabetes and other comorbidities (table 4). The multivariable analysis included 168 936 records (65 670 non-survivors) that had no missing data for covariates. The odds of in-hospital death for males are 1.23 times that of females, and for those in the North and Northeast regions were, respectively, 1.83 and 1.48 times that of patients in the Southeast. The Cox proportional-hazards model included 176 559 records (64 809 non-survivors). Variables associated with in-hospital death were age, sex, ethnoracial self-classification, region, ICU care and various comorbidities. Age, however, had the highest hazard ratios, ranging from 1.67 (95% CI: 1.49 to 1.89) for those aged 20–39 to 5.51 (95% CI: 4.91 to 6.18) for those 80 or older, compared with patients younger than 20 years (table 4).

Table 4

ORs and HRs for death among hospitalised patients with confirmed COVID-19

Discussion

This study described demographic, clinical and radiographic characteristics, treatment, case outcome and risk factors associated with in-hospital death of 522 167 patients hospitalised with confirmed COVID-19 in Brazil. Results show that 56.7% were discharged, 0.002% died of other causes, 30.7% died of causes associated with COVID-19 and 10.2% remained in the hospital as of 14 December. Patients were mostly older than 40 years, predominantly from the Southeast region, with about one-fourth needing to travel to a different municipality for hospitalisation. At least one comorbidity was present in 64.5% of patients and in 76.8% of the non-survivors. From illness onset, the median time to hospital and ICU admission was 6 and 7 days, respectively; 15 days to death (17 to those admitted to the ICU) and 15 days to hospital discharge (18 days to those admitted to the ICU). Risk factors for in-hospital death were older age, being male, of Black/Brown ethnoracial self-classification, living in the North or Northeast regions, with a history of ICU admission and various co-morbidities.

Our results can be analysed in light of findings from other countries. The use of mechanical ventilation was higher in Brazil (62.2% among all patients, 75.6% of the non-survivors) compared with patterns described for China (ranging from 17.2% to 38.7% of patients),8–10 and Germany (17% of patients),11 but lower than Italy (81.7% of all patients).12 While ICU mortality in Italy was 26%,12 in Brazil, it was 51.8%. In Brazil, 33.0% of hospitalised patients were admitted to the ICU, against 16% in Italy,13 and 19% in France.14 In-hospital death was observed in 18.1% of patients in France,14 22% in Germany11 and 30.7% in Brazil. The time from illness onset to hospitalisation in China9 was 11 days (IQR, 8–14), but much shorter in Brazil, 6 days (IQR, 3–9). The length of hospital stay, however, was about 4 days longer in China.9 15 These comparisons need to be taken with caution. First, studies from China had a smaller sample size, and regional variability is very large, as reported for France.14 In Brazil, for example, in-hospital death varied from 25.9% in the South region to 36.4% in the Northeast, and ICU mortality from 48.0% in the Southeast to 66.5% in the North. The time from illness to hospitalisation was also longer in the North and Center-West regions, 7 days (IQR, 4–10).

Our results confirm previous findings regarding symptoms and comorbidities. Hypertension, a very common comorbidity in China9 could not be measured from our data, but over one-third of the adult population in Brazil has that condition.16 In Brazil, 35.5% of the patients reported no comorbidities, while in New York City, this number was 6.1%,17 and in China, it was 52%.9 Diabetes was reported in 19% of patients in China,9 33.8% in New York City17 and 25.7% in Brazil. Part of these differences reflects the disease burden in each locality, but also the lack of standardised data collection (eg, conditions systematically collected in one country and only reported in the ‘other’ category in another country).

The observed associations of age, sex, obesity and diabetes with in-hospital death corroborate previous findings.5 6 14 The higher risk among non-white patients was previously reported in Brazil, the UK and the USA.18–22 In Brazil, this reflects structural inequalities that made large fractions of the population more vulnerable to COVID-19 (eg, people living in areas with precarious infrastructure, overcrowded households, regions with low supply of physicians and hospital services and who depend on informal labour).23 24 Those inequalities are also captured by a higher HR in the North and Northeast regions, where Brazil consistently reports worse socioeconomic indicators.25 Currently, the North region has the lowest rates of hospital beds, ICU beds and physicians per person.26 Indeed, the region had the worst indicators in terms of mortality and time to hospitalisation, and patients who were hospitalised in a municipality different from the one where they live had to travel 122.0 km (IQR, 58.3–258.6), while those in the Southeast travelled 22.3 km (IQR, 16.1–36.3). Hospitalisation in a different municipality occurs when the place of residence has no hospitals, has no available hospital beds or when the closest hospital is actually outside the municipality of residence. In Brazil, the size of municipalities varies widely: 23% have 5000 residents or less, and 5% have 100 000 or more residents. Of the 5570 municipalities, 37% and 75% have no hospitals and no ICU care, respectively. A regionalisation process guarantees that all the population has access to hospital care.27 However, when hospitals reach capacity, as was observed in several cities in Brazil in late April and May, municipalities without hospitals and ICU units are unable to provide proper care, which may have contributed to higher COVID-19 mortality. Therefore, risk factors for in-hospital mortality due to COVID-19 expose local and structural inequalities.

This study has some limitations. First, we used administrative records captured in structured surveillance forms. Those lack details available in hospital medical records and may have accuracy and completeness problems. In addition, it limits the types of comorbidities and symptoms reported, as those listed under the ‘other’ category were not systematically collected from all cases (eg, loss of taste and smell). Second, we did not have access to laboratory results other than COVID-19 tests (eg, complete blood count). While this does not change any of our results, they would allow for a better characterisation of the clinical course of the disease. Third, 23.1% of patients did not report information on ethnoracial self-classification, 10.8% did not have information on ICU admission and for 2.1%, there was no information on clinical outcome. This is not uncommon in the analysis of administrative records.22 Here, we report all records and included an additional category (ethnoracial self-classification not reported) in the risk factor analysis.

Despite these limitations, this study provides a comprehensive description of characteristics, outcomes and risk factors for mortality of patients hospitalised with COVID-19 in Brazil, and the largest cohort of patients so far analysed (n=5 22 167). Results shed light on commonalities and differences between Brazil and other countries affected by COVID-19, and highlight inequalities in disease outcomes. Most importantly, our results could be used to inform coordinated actions to target those who currently bear the highest morbidity and mortality burden. Brazil has a network of almost 270 000 community health workers that reach out to about 70% of the Brazilian population. These agents could actively identify vulnerable people who face higher risk of mortality, could act as agents of information to sensitise the population and boost adherence to control measures (eg, use of masks) and could continue to deliver community-based primary care services that have been, for the most part, interrupted by the pandemic. These agents will also be important to support the delivery of vaccination to the most vulnerable. Leveraging and strengthening the existing network of primary healthcare is paramount to contain the sustained and unequal burden of COVID-19 in Brazil.

Data availability statement

Data are available in a public, open access repository. The SIVEP-Gripe data set is publicly available on the Ministry of Health’s DATASUS website (https://opendatasus.saude.gov.br/nl/dataset).

Acknowledgments

We would like to thank Nicholas Arisco, MS, for technical assistance.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors MCC and GVAdF conceived the study, were responsible for data analysis and interpretation. MCC wrote the manuscript. GVAdF and EMM acquired the data. GVAdF and SG were responsible for data curation. SG ran the statistical models and contributed to writing. All authors edited and approved the final version of the manuscript.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.