key: cord-0886199-akgd14h0
authors: Undurraga, Eduardo A.; Chowell, Gerardo; Mizumoto, Kenji
title: COVID-19 case fatality risk by age and gender in a high testing setting in Latin America: Chile, March–August 2020
date: 2021-02-03
journal: Infect Dis Poverty
DOI: 10.1186/s40249-020-00785-1
sha: b9cc7f71d36c5537156c873b4e4118426d8648ea
doc_id: 886199
cord_uid: akgd14h0

BACKGROUND: Early severity estimates of coronavirus disease 2019 (COVID-19) are critically needed to assess the potential impact of the ongoing pandemic in different demographic groups. Here we estimate the real-time delay-adjusted case fatality rate across nine age groups by gender in Chile, the country with the highest testing rate for COVID-19 in Latin America. METHODS: We used a publicly available real-time daily series of age-stratified COVID-19 cases and deaths reported by the Ministry of Health in Chile from the beginning of the epidemic in March through August 31, 2020. We used a robust likelihood function and a delay distribution to estimate real-time delay-adjusted case-fatality risk and estimate model parameters using a Monte Carlo Markov Chain in a Bayesian framework. RESULTS: As of August 31, 2020, our estimates of the time-delay adjusted case fatality rate (CFR) for men and women are 4.16% [95% Credible Interval (CrI): 4.09–4.24%] and 3.26% (95% CrI: 3.19–3.34%), respectively, while the overall estimate is 3.72% (95% CrI: 3.67–3.78%). Seniors aged 80 years and over have an adjusted CFR of 56.82% (95% CrI: 55.25–58.34%) for men and 41.10% (95% CrI: 40.02–42.26%) for women. Results showed a peak in estimated CFR during the June peak of the epidemic. The peak possibly reflects insufficient laboratory capacity, as illustrated by high test positivity rates (33% positive 7-day average nationally in June), which may have resulted in lower reporting rates. CONCLUSIONS: Severity estimates from COVID-19 in Chile suggest that male seniors, especially among those aged ≥ 70 years, are being disproportionately affected by the pandemic, a finding consistent with other regions. The ongoing pandemic is imposing a high death toll in South America, and Chile has one of the highest reported mortality rates globally thus far. These real-time estimates may help inform public health officials' decisions in the region and underscore the need to implement more effective measures to ameliorate fatality. [Image: see text] [Image: see text]

The coronavirus disease 2019 (COVID-19) pandemic has strained or overwhelmed health systems across the world [1, 2] , with about 25.5 million COVID-19 cases and 850 thousand deaths as of August 31, 2020 [3, 4] . The first case in South America of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, the cause of COVID-19, was reported on February 25 in São Paulo, Brazil, a travel hub in the region [5] . A few weeks later, countries in the region had imposed major epidemic control measures, such as closed borders, restricted travel, closure of schools and universities, and enforced lockdowns [6, 7] . Despite these measures, the ongoing coronavirus pandemic has already imposed a high toll on most countries in South America, killing thousands in Brazil (121 100), Peru (28 800), Colombia (19 700) , Chile (15 700), Argentina (8700), Ecuador (6600), and Bolivia (5000), thus far [3] . In addition to alreadystrained healthcare systems, other factors have affected the dynamics of the pandemic in the region, including migration, sociopolitical crises, struggling economies, other infectious disease outbreaks, poor leadership, and challenges tied to the implementation of social distancing, hygiene, and lockdown strategies, due to inadequate water and sanitation infrastructure and precarious living conditions [6, [8] [9] [10] [11] [12] [13] [14] [15] .

The Ministry of Health reported the first COVID-19 case in Chile on March 3, 2020 [16] . The government activated a quick response effort, announcing restrictions on large gatherings on March 13 , and subsequently, closure of all daycares, schools, and universities (March 16), border controls, telework recommendations (March 18), closure of non-essential businesses (March 19) , national night curfews (March 22) , and a strategy of intermittent localized lockdowns in selected municipalities since March 26 [17] . Estimates for the early stage of the epidemic showed sustained transmission with an estimated reproduction number R 0 = 1.8 [95% confidence interval (CI): 1.6-1.9] [18] . Research suggests that containment measures at the start of the epidemic significantly slowed down the spread of the virus in Chile [18] . Early models estimated that a surge in the number of ill patients could overwhelm treatment capacity by June unless aggressive control measures were implemented [19] . Control measures prevented the epidemic's growth in some regions, particularly during the first months of the epidemic [20] [21] [22] .

Chile saw a rapid increase in reported COVID-19 cases during May, reaching a peak 7-day average incidence of about 6000 cases per day (31 per 100 000 population) by mid-June [16] . Health system capacity reached saturation levels by the end of June, with 95% intensive care unit occupancy in Greater Santiago and 89% at the national level [23, 24] . As of August 31, the Ministry of Health's surveillance system EPIVIGILA has recorded a total of 444 921 COVID-19 cases and 15 756 deaths [16, 25] . Most of these reported COVID-19 cases and deaths have been confirmed by reverse transcription-polymerase chain reaction (RT-PCR) assay. However, these reported cases and deaths also include some clinically diagnosed COVID-19, particularly during the peak of the epidemic when the healthcare system reached saturation [23] , and laboratories could not satisfy the demand for diagnostic tests. In June, 7-day average positivity rates of the RT-PCR tests reached 33% nationally. Evidence shows that the circulation of other respiratory viruses, such as influenza, has been negligible in Chile and elsewhere in the southern hemisphere [26] . With 130 total RT-PCR tests per 1000 population by August 31, Chile has tested at a higher rate than any other Latin American country [27, 28] .

During the first two months of the epidemic, the crude case fatality rate (CFR), the number of cumulative deaths over the number of cumulative cases, in Chile (1.2%) remained well below the global average (6.5%) [4, 16] . Because initial COVID-19 cases in Chile occurred among relatively low-risk groups (20-60 years of age), several infectious disease experts warned about authorities' excessive confidence over early success [29] . The CFR is a commonly used estimate of the severity of an epidemic [30] . It provides a reliable benchmark for public health officials to decide the intensity and duration of interventions to mitigate or suppress an epidemic [31] . Nevertheless, obtaining CFR estimates during an epidemic is challenging, as CFR is typically affected by right censoring and ascertainment bias [32] [33] [34] [35] [36] . Right-censoring occurs because of a time delay between the onset of symptoms and death, thus underestimating CFR. Underascertainment of cases occurs because mild or asymptomatic COVID-19 cases often go undetected by disease surveillance systems, which are not designed to detect all infections [33, 37, 38] . SARS-CoV-2 infection can result in a broad spectrum of clinical outcomes, including asymptomatic infection, mild symptoms, hospitalization, or death [39] [40] [41] .

Here we provide real-time estimates of adjusted agespecific CFR during the COVID-19 epidemic in Chile, March through August 2020, to gauge the severity of the SARS-CoV-2 epidemic. These CFR estimates could inform critical decisions by public health officials in Chile and Latin America.

We obtained daily cumulative numbers of reported COVID-19 cases and deaths stratified by age group and gender from March 3, 2020, through August 31, 2020, from the Ministry of Health's surveillance system EPIV-IGILA [16, 25] . The datasets used are available from the Ministry of Science's COVID-19 repository [25] . COVID-19 death definition follows the World Health Organization's updated guidelines for COVID-19 based on the International Statistical Classification of Diseases and Related Health Problems (ICD) standards [42] . Different groups had different starting times in the time series, which correspond to the day when the first death was reported. 

For the estimation of CFR in real-time, we employed the delay from hospitalization to death, h s . We assumed h s = H(s) -H(s-1) for s > 0 where H(s) is a cumulative density function of the delay from hospitalization to death and follows a gamma distribution with mean 10.1 days and standard deviation (SD) 5.4 days [43] . Let π a,ti be the time-delay adjusted case fatality ratio on reported day t i in area a, the likelihood function of the estimate π a,ti is where c a,t represents the number of new cases with reported day t in area a, and D a,ti is the cumulative number of deaths until reported day t i in an area a [44, 45] . Among the cumulative cases with reported day t in area a, D a,ti have died, and the remainder have survived the infection. The contribution of those who have died with biased death risk is shown in the middle parenthetical term. The contribution of survivors is presented in

c a,t −D a,t i the right parenthetical term. We assume that D a,ti is the result of the binomial sampling process with probability π a,ti . We estimated model parameters using a Monte Carlo Markov Chain (MCMC) method in a Bayesian framework. We evaluated the convergence of MCMC chains using the potential scale reduction statistic [46, 47] . Estimates and 95% credibility intervals (CrI) for these estimates are based on each parameter's posterior probability distribution and based on the samples drawn from the posterior distributions.

All statistical analyses were conducted in R version 3.6.1 (R Foundation for Statistical Computing, Vienna, Austria) using the 'rstan' package. This research was considered exempt from ethical review because it is limited to secondary data analysis based on publicly available data [16, 25] , which contains no information that can directly or indirectly identify an individual. All data used in this analysis are publicly available [25] .

The Fig. 1a , b). We found a similar pattern for women, except that most deaths were reported among women aged 80 years and older (44.4%).

The proportion of male cases is higher than 50% across age groups, except for 10-19 year-olds, 70-79 year-olds, and those aged 80 years and above (χ 2 test, P-value < 0.0001, Fig. 1c ). The proportion of male deaths is higher than 50% except for those aged 80 years and above (χ 2 test, P < 0.0001, Fig. 1d ). Figure 1e presents the cumulative morbidity ratio by gender and age group. The ratio of illness among males is significantly higher than among females across most age groups (proportion test, P < 0.000), except for 10-19 year-olds (proportion test, P < 0.00) and 20-29 year-olds (proportion test, P = 0.01). We found no statistically significant differences for 0-9 year-olds (proportion test, P = 0.19).

We found that mortality per 100 000 population directly caused by COVID-19 among males is higher than among females for all age groups 20 years and above (proportion test, P < 0.05, Fig. 1f ). We found no statistically significant differences for individuals aged 0-19 years.

The reported cumulative COVID-19 cases and deaths by age group and gender (A through J) over time are shown in Fig. 2 . The figure suggests cumulative cases of COVID-19 are growing faster than cumulative deaths. The growth curve for cumulative cases across all age groups increases exponentially after around day 30 (April 30, 2020) until around day 80 (June 19, 2020). Exponential growth in cumulative deaths for all age groups occurs approximately after day 45 (May 15, 2020).

The observed and model-based posterior estimates of the crude CFR of COVID-19 by age group (A-J) and time-delay adjusted CFR by age group (K-T) are shown in Fig. 3 . Black dots show crude case-fatality ratios, and light and dark indicate 95% and 50% credible intervals for posterior estimates, respectively. Overall, our model-based crude CFR fitted the observed data well. Crude CFR for all age groups (J) increased at an early stage of the epidemic, and once peaked around April 25 and was followed by an increase in CFR. There is an increasing CFR trend across most age groups, starting around May 10, 2020. Our model-based posterior estimates for the time-delay adjusted CFR are substantially higher than the crude observed CFR at an early stage. Both values are converging as the epidemic progresses.

The time-delay adjusted CFR results by gender and age group, as of August 31, 2020, are shown in Table 2 and Fig. 4 . Time-delay adjusted CFR for men and women are 4.16% (95% CrI: 4.09-4.24%) and 3.26% (95% CrI: 3.19-3.34%), respectively, while overall national estimate is 3.72% (95% CrI: 3.67-3.78%). Among men, senior citizens appear to be severely affected; the adjusted CFR is 10.16% (95% CrI: 9.76-10.60%) for men aged 60-69 years, 28.35% (95% CrI: 27.49-29.35%) for those aged 70-79 years, and 56.82% (95% CrI: 55.25-58.34%) for those 80 years old and above. We observe a similar pattern for women. The adjusted CFR is 6.44% (95% CrI: 6.10-6.81%) for women aged 60-69 years, 18.18% (95% CrI: 17.40-18.93%) for those aged 70-79 years, and 41.10% (95% CrI: 40.02-42.26%) for women aged 80 years old or more.

Few studies have estimated the time-delay adjusted CFR by age group for COVID-19 in Latin America [48] , a region with relatively little attention during the ongoing coronavirus pandemic. Consistent with other recent COVID-19 research [31, 33, 49] , our results show that the COVID-19 epidemic in Chile has disproportionately affected seniors, especially those aged ≥ 70 years. These results suggest that an aging population could exacerbate the fatality impact of COVID-19 [50] , similar to influenza and respiratory syncytial virus [51] , and consistent with data available from Italy [50] . The comparatively low CFR observed in Chile during the early stages of the epidemic [29] probably reflected the initial cases' age structure and socioeconomic status. During the first weeks of the pandemic, COVID-19 cases occurred among relatively young age groups, with most transmission occurring among individuals between 20 and 60 years of age, and in highincome communities with better access to healthcare and lower prevalence of risk factors for severe COVID-19 [13, 14, 16, 21, 52, 53] .

Our latest adjusted CFR estimates for those aged 80 and over reach values as high as 56.82% (95%CrI: 55.25-58.34) for men and 41.10% (95%CrI: 40.02-42.26%) for women. Adjusted CFR for men aged 80 years of age and older is 2.0-fold higher than our estimates for men aged 70-79 years and 5.6-fold higher than men aged 60-69.

Adjusted CFR for women 80 years of age and older is 2.3-fold higher than our estimates for women aged 70-79 years and 6.4-fold higher than women aged 60-69. The large majority (68.4%) of COVID-19 cases occurred in Greater Santiago during the study period. As seen for the overall population starting in April until about day 75 (June 14), the upward trend in the crude CFR suggests that viral transmission was spreading to more vulnerable populations. Interestingly, throughout March, most transmissions occurred among relatively young, betteroff populations. In Greater Santiago, localized lockdowns were put in place starting March 28 through mid-April in seven municipalities, six of which are among the richest in Chile [54] . Viral transmission moved towards lowerincome municipalities during April [16, 21] . Among these economically disadvantaged populations, social distancing measures are harder to comply with due to a higher proportion of the population participating in the informal economy, greater population density and overcrowding, more deficient sanitary infrastructure, and lower quality healthcare [14, 21, 55] . Estimates of human mobility during the city-wide lockdown based on cellphone data across municipalities show that the lockdown was less effective in lower-income municipalities [14, 56] . The higher the rate of multidimensional poverty (based on education, health, labor, housing, and social capital), the less effective was the lockdown in reducing mobility, with up to a three-fold difference in mobility reductions between rich and poor municipalities in Santiago. The disproportionate impact that the COVID-19 pandemic has had on lower-income or more vulnerable populations has also been documented elsewhere [5, 57, 58] . An upward trend in the crude CFR could also have resulted from an increasing number of unreported cases due to a saturated testing capacity [55] . Chile's testing capabilities were greatly expanded during the epidemic, partly as a coordinated effort from the Ministry of Science to include testing from public and private laboratories, including universities that converted research labs to PCR testing labs. Processing capacity grew from a few hundred tests per day to up to 35 000 per day (cumulative of about 130 tests per 1000 population) so far [25] . Despite this massive increase in testing capabilities, the 7-day average test positivity (positive tests/total tests) reached 33% nationally and 48% in Greater Santiago during the mid-June peak of the epidemic. This dramatic increase in positive test rates suggests there may be a substantial proportion of under-diagnosed COVID-19 cases in Chile, at least during the peak of transmission in Greater Santiago [59] .

The health system capacity also reached saturation levels during the pandemic's peak by the end of June, with 95% intensive care unit occupancy in Greater Santiago and 89% at the national level [23, 24] . A major disease outbreak, such as COVID-19, can strain the healthcare system and hospital resources [19] . This excess demand may potentially degrade treatment quality or lead to sub-optimal treatment decisions, such as diversion from hospitalization or early discharge, and cause delays in laboratory work and reporting, as has been observed in other public health emergencies [60] [61] [62] . However, our data do not allow us to test this hypothesis.

The downward trend in the adjusted CFR at the early stage of the study period may have been influenced by reporting delays. The observed differences in our crude and adjusted CFR estimates are directly due to the timedelay, which we assume fixed during the epidemic.

The small proportion of men (39.2%) among COVID-19 cases in people aged 80 and over is probably attributable to the relatively small male population size for that age group; men represent only 37% of the population > 80 in Chile [63] . Life expectancy at birth in Chile is 77.2 years for men and 82.1 years for women. Life expectancy at age 60 is 21.8 years for men and 24.9 years for women. China and the United States have reported higher mortality among men [64, 65] ; our data provide the opportunity to examine the CFR by gender and age. Our study has limitations. First, our CFR estimates are probably affected by under ascertainment, as has been estimated elsewhere [33, 37, 38] . Under ascertainment may have pushed our estimates upwards if the surveillance system captured the most severe cases. This bias may be particularly problematic during the peak of the epidemic from mid-May through the end of June. Infectious diseases with a substantial share of asymptomatic or mild infections, such as COVID-19, may be more accurately characterized by infection-fatality risk (deaths / infected people). However, those data are not yet available in Chile. Second, we used laboratory-confirmed and probable COVID-19 cases in our analysis, which may result in an overestimation of COVID-19 cases and thus may have pushed our CFR estimates downwards. The demand for diagnostic tests exceeded laboratory capacity in Chile during the peak of the epidemic, as suggested by high rates of positive tests [16] . Test positivity reached a 7-day average of 33% nationally and 48% at the peak of the epidemic in Greater Santiago, suggesting that a substantial proportion of cases were left without testing. In our favor, the circulation of other respiratory viruses has been of no importance during the COVID-19 pandemic in Chile and elsewhere [26] in the southern hemisphere.

Using real-time epidemiological data from a high COVID-19 testing setting in Latin America, we found that the COVID-19 epidemic in Chile has disproportionately affected seniors, especially those aged ≥ 70 years, suggesting an aging population could exacerbate the death toll brought by the emergence of SARS-CoV-2. COVID-19 is already imposing a high death toll in Latin America. These real-time estimates may help inform public health officials' decisions in the region and underscore the need to implement more effective measures to ameliorate mortality, especially among those facing the highest risk of death. 

The impact of COVID-19 and strategies for mitigation and suppression in low-and middle-income countries

Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe

An interactive web-based dashboard to track COVID-19 in real time

Coronavirus Disease (COVID-19) Outbreak Geneva: WHO

Epidemiological and clinical characteristics of the COVID-19 epidemic in Brazil

COVID-19 in Latin America

COVID-19 response in Latin America

The unfolding migrant crisis in Latin America

Global behaviors and perceptions at the onset of the COVID-19 pandemic

COVID-19 and dengue, co-epidemics in Ecuador and other countries in Latin America: pushing strained health care systems over the edge

South America prepares for the impact of COVID-19

The Economy in the Time of CoVID-19 Washington DC: World Bank

Epidemiology, causes, clinical manifestation and diagnosis, prevention and control of coronavirus disease (COVID-19) during the early outbreak period: a scoping review

CoVID-19 has exposed how 'the other half' (still) lives

Political and institutional perils of Brazil's COVID-19 crisis

Cifras Oficiales COVID-19 Santiago: Ministry of Health

Plan de acción por coronavirus Santiago: Ministerio de Salud, Gobierno de Chile

Transmission dynamics and control of COVID-19 in Chile

A dynamic modeling tool for estimating healthcare demand from the COVID19 epidemic and evaluating population-wide interventions

Effectiveness of localized lockdowns in the SARS-CoV-2 pandemic

All things equal? Heterogeneity in policy effectiveness against COVID-19 spread in Chile

Impact of small-area lockdowns for the control of the COVID-19 pandemic

Epidemic trends, public health response and health system capacity: the Chilean experience in four months of the COVID-19 pandemic

Sociedad Chilena de Medicina Intensiva SOCHIMI. Encuesta nacional sobre ocupación de unidades críticas durante contingencia CoVID-19

Universidad Finis Terrae

Decreased influenza activity during the COVID-19 pandemic-United States

Chile charts new path with rolling lockdowns, immunity cards

Coronavirus pandemic (COVID-19)

Chile doctors fear complacency over Covid-19 after initial successes. The Guardian. 2020

Novel framework for assessing epidemiologic effects of influenza epidemics and pandemics

Estimating risk for death from 2019 novel coronavirus disease, China

Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (COVID-19)

Estimation of SARS-CoV-2 mortality during the early stages of an epidemic: a modeling study in Hubei, China, and six regions in Europe

Methods for estimating the case fatality ratio for a novel, emerging infectious disease

Nonparametric estimation of the case fatality ratio with competing risks data: an application to Severe Acute Respiratory Syndrome (SARS)

Potential biases in estimating absolute and relative case-fatality risks during outbreaks

Real-time estimation of the risk of death from novel coronavirus (COVID-19) infection: inference using exported cases

Association of public health interventions with the epidemiology of the COVID-19 outbreak in Wuhan

Clinical characteristics of coronavirus disease in China

Clinical features and short-term outcomes of 18 patients with corona virus disease 2019 in intensive care unit

Estimating clinical severity of COVID-19 from the transmission dynamics in Wuhan, China

International guidelines for certification and classification (coding) of COVID-19 as cause of death Geneva: WHO. 2020. • fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year submit your research ? Choose BMC

Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship

Early epidemiological assessment of the virulence of emerging infectious diseases: a case study of an influenza pandemic

Dynamics of the pneumonic plague epidemic in Madagascar

Markov chain Monte Carlo: stochastic simulation for Bayesian inference

Inference from iterative simulation using multiple sequences

Risk of death by age and gender from CoVID-19 in Peru

The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19)-China

Demographic science aids in understanding the spread and fatality rates of COVID-19

The burden of influenza and respiratory syncytial viruses in Japan, 2006-2014: a region-and agespecific excess mortality study

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

Does comorbidity increase the risk of patients with COVID-19: evidence from meta-analysis

Sistema Nacional de Información Municipal. Datos Municipales Santiago. 2020

Desiguales: Orígenes, cambios y desafíos de la brecha social en Chile: Programa de las Naciones Unidas para el Desarrollo (PNUD)

Informe sobre la evolución de la epidemia de Covid-19 en Chile, 30 de mayo 2020 Santiago de Chile

Evidence mounts on the disproportionate effect of COVID-19 on ethnic minorities

The disproportionate impact of COVID-19 on racial and ethnic minorities in the United States

The effect of COVID-19 on the economy: evidence from an early adopter of localized lockdowns

Effects of response to 2014-2015 Ebola outbreak on deaths from malaria, HIV/AIDS, and tuberculosis, West Africa

The changing epidemiology of malaria elimination: new strategies for new challenges

Health-care worker mortality and the legacy of the Ebola epidemic

Estimaciones y proyecciones de la población de Chile

Presenting characteristics, comorbidities, and outcomes among 5700 patients hospitalized with COVID-19 in the New York City area

The Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19)-China

Not applicable.

GC and KM conceived the initial study idea. KM implemented statistical analysis. GC, KM, EU contributed to the study implementation, interpretation of results, and writing of the manuscript. All authors read and approved the final manuscript. ] to EU. The study sponsors had no role in the study design, collection, analysis, and interpretation of data; in the writing of the manuscript; and in the decision to submit the manuscript for publication.

The datasets used and analyzed during the current study are available from Base de Datos COVID-19 repository, http://www.minci encia .gob.cl/covid 19.

This research was considered exempt from ethical review because it is limited to secondary data analysis based on publicly available data [16, 25] , which contains no information that can directly or indirectly identify an individual.

Not applicable.

The authors (EU, GC, KM) declare that they have no competing interests.