SUMMARY
Background Direct and indirect COVID19-related mortality is uncertain. This study investigated COVID19-related and all-cause deaths among middle-aged and older adults during the first wave of COVID19 epidemic period, assessing mortality risks by pre-existing socio-demographic and medical underlying conditions.
Methods Population-based cohort study involving 79,083 individuals ≥50 years-old in Tarragona (Southern Catalonia, Spain). Baseline cohort characteristics (age/sex, comorbidities and medications/vaccinations history) were established at study start (01/03/2020) and main outcomes were COVID19-related deaths (occurred in patients diagnosed with the disease) and all-cause deaths occurred among cohort members between 01/03/2020-30/06/2020. Mortality risks were assessed by Cox regression analyses.
Results Cohort members were followed for 1,356,358 persons-weeks, occurring 576 all-cause deaths (124 COVID19-related). All-cause mortality rate was 42.5 deaths per 100,000 persons-week, being 22.8 in healthy/unrelated-COVID19 subjects, 236.4 in COVID19-excluded/PCR-negative subjects, 493.7 in COVID19-compatible/PCR-unperformed subjects and 4009.1 in COVID19-confirmed patients. In multivariable analyses, increasing age, sex male, nursing-home residence, cancer, neurologic, cardiac or liver disease, receiving diuretics, systemic corticosteroids, proton-pump inhibitors and benzodiazepines were associated with increased risk of all-cause mortality; conversely, receiving renin-angiotensin inhibitors and statins were associated with reduced risk. Age/years, sex male and nursing-home residence were strong predictors for COVID19-related mortality, but none comorbidity appeared significantly associated with an increased risk.
Conclusion Apart from direct COVID19-related deaths (which represented almost 22% of all-cause mortality), theoretically COVID19-excluded patients (PCR-negative) suffered considerable greater all-cause mortality than healthy/unrelated-COVID19 subjects, which could explain, in part, the large excess deaths observed across the COVID19 pandemic.
INTRODUCTION
A year ago, following an outbreak of pneumonia in Wuhan (China), a new coronavirus was detected. Said virus, later named Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), is the cause of the disease designated by the World Health Organization (WHO) as COVID19 in February 2020.[1] Currently, the pandemic caused by this virus continues to be a threat to public health worldwide.[2]
Available publications have reported different percentages of asymptomatic or mild cases of COVID19 but, even if they account for around 80%,[2,3] due to the high infectivity of the virus, the number of disease’s related deaths is very high.[2] The mortality rates reported are very heterogeneous, depending on the countries (with important differences between regions, even within the same country), and also on the availability / policy of conducting diagnostic tests.[2]. For example, in Lombardy (Italy) the case fatality rate ranged between 1.6% and 18.3%,[4] while the Chinese Centre for Disease Control and Prevention, in a cohort of 72,314 cases, estimated a 2.3% case-fatality rate (14.8% for patients ≥80 years)[5]. Most articles report in-hospital mortality or case fatality rates (mortality rate among confirmed cases); it must be taken into account that infection fatality rate (which consider all infected individuals) would be lower. On the other hand, many deaths from COVID19 occur in undiagnosed patients. Lastly, we have to keep in mind the “excess mortality” (deaths caused by other conditions, linked to a delay in care, overburden of the health system and socioeconomic determinants of health).[6,7]
In the case of Spain, the latest scientific-technical report of the Ministry of Health (dated November 12) informed a fatality rate of 11% in patients ≥50 years (reaching 19% in the group ≥70 years compared to 0.45% in aged 50-69 years).[8] However, the estimates that took into account the results of the seroprevalence study (done in May)[9] reported much lower rates: 0.04% in the younger (50-69 years) and 4.1% in those ≥70 years, with a global rate for ≥50 years of 1.5%.[8] According to the Statistics National Institute (INE), as compared with 2019, almost a 20% increased all-cause death occurred in Spain during 2020 (surely as consequence of COVID19 pandemic).[10]
Most reviews agree that age and gender are well-established risk factors of mortality in COVID19 patients,[11,12] but there is no clear consensus on the specific contribution of pre-existing conditions to the risk of mortality from the disease.[13] Knowing the relative contribution of factors related to the patient (including socio-demographic and medical aspects) to mortality is very important to identify the individuals with greater risk and guarantee a more efficient use of resources.
This study investigated COVID19-related and all-cause mortality among middle-aged and older adults in the region of Tarragona across the first epidemic wave (March-June 2020), assessing the possible association between previous conditions (demographic, comorbidities, chronic medications’ use) and risk of death (COVID19-related and/or any cause) among middle-aged and older adults (who support the greatest burden of severe disease).
METHODS
Design, setting and study population
This is part of a large population-based retrospective-prospective cohort study (initiated in March 2020 and called COVID19 TARRACO Cohort Study) involving 79,083 individuals 50 years or older in the region of Tarragona (Southern Catalonia, Spain).
The design, setting and study population have been extensively described elsewhere, in prior articles that evaluated the epidemiology of COVID-19 during the first 8-12 weeks of epidemic period in the study area.[14,15]
The 79,083 cohort members were all individuals ≥50 years-old affiliated to the 12 Primary Care Centres (PCCs) managed by the Institut Català de la Salut in the study area (Tarragonés, Alt Camp and Conca de Barberà counties) and represented approximately 75% of overall inhabitants in this age strata according to census data.[16] This report focuses on cohort members who died (by COVID-19 or any cause) in the study setting between March 1 and June 30, 2020. The study was approved by the Ethical Committee of the Institution (Ethics Committee IDIAP Jordi Gol, Barcelona, file 20/065-PCV) and was conducted in accordance with the general principles for observational studies.[17] The study was determined to be exempt for informed consent under the public health surveillance exception.
Data sources
Former CAPAMIS Research Database, an institutional clinical database used for earlier cohort studies in the region,[18] was updated and exploited as the prime data source in this COVID-19 epidemiological report. Summarily, this database collects inputs from the electronic system (e-CAP) of primary care clinical records, that includes administrative and clinical data, coded according to the International Classification of Diseases 10th Revision (ICD-10). It allowed us to identify socio-demographic aspects and previous conditions/comorbidities among cohort members in order to define their baseline profile at study onset (01/03/2020).
Two new electronic alerts, related to COVID-19’s laboratory results and ICD-10 codes for COVID-19 suspicion (B34.2, B97.29) were introduced to the e-CAP system at the beginning of the epidemic; later, both information sources were bonded to create the baseline research database utilised in this study.
Vital status was considered according to administrative data periodically updated in the e-CAP system. In addition, complementary data registered during emergency visits and/or hospital-stays, were used to complement vital status data in COVID-19 cases identified across study period. These data were transferred to a standardised form, after retrospective review of electronic records carried out by a research-trained group of family physicians and were afterwards linked with the prime research database containing baseline characteristics of 79,083 cohort members.
Outcome definitions
Primary outcome was death from any cause (COVID19-related or not) occurred among cohort members during the 4-month of study period (01/03/2020-30/06/2020). Laboratory-confirmed COVID-19 case was defined when a cohort member tested positive in SARS-COV-2 using Reverse Transcription-Polymerase Chain Reaction (RT-PCR) or serological tests, according institutional guidelines.[19] COVID19-related death was considered when the patient died from any cause within the first 30-days after the onset of the disease or at any time during hospital-stay. Death by COVID-19 (likely due to COVID-19) or death with COVID-19 (likely due to other concomitant causes) were differentiated according to family physician’s criteria after checking hospital and primary care clinical records in each of deceased COVID19 cases.
Covariates
Baseline covariates were age, sex, type of residence (community-dwelling or nursing-home), vaccinations’ history, previous comorbidities and chronic medications’ use; they were considered according to data registered in the e-CAP system.
Criteria used to define comorbidities and chronic medication use are extensively described in an appendix elsewhere.[15] Unrecorded comorbidities, medications or vaccinations were considered as absent.
In addition, at the end of the study period, according to their COVID-19 relation, cohort members were classified as: COVID19-unrelated/healthy, COVID19-excluded (PCR/serological test performed with negative result), COVID19-compatible (clinical suspicion alone, without any test performed) and COVID19-confirmed (positive PCR or serological test).
Statistical analysis
Descriptive statistics of socio-demographic and clinical variables include frequencies/ percentages for categorical variables; and means/ standard deviations (SD) for continuous variables. Chi-squared or Fisher test (categorical variables) and Student’s t test (continuous variables) were used for comparison of differences between groups.
Mortality Rates (MRs) for COVID19-related and all-cause deaths were calculated per 100,000 persons-week, considering in the denominator the sum of the persons-time contributed for each cohort member during the study period (17.4 weeks). The confidence intervals (CIs) for the MRs were calculated assuming a Poisson distribution for uncommon events.
Cox regression analyses were used to calculate unadjusted, age&sex-adjusted and multivariable-adjusted hazards ratios (HRs) and estimate the association between baseline conditions and the time to the first outcome (COVID19-related or all-cause death) occurred among cohort members throughout the study period (from 1st March to 30th June 2020). All exposure above mentioned covariates (i.e., age, sex, residence, comorbidities/underlying conditions, chronic medications’ use and vaccinations’ history) were considered for multivariable Cox models; the method to select a subset of covariates to include in final models was the purposeful selection.[20] Final models include both significant or confounder variables and all covariates judged clinically or epidemiologically relevant. We performed a main analysis including the total study cohort (N=79,083) and two subgroup analyses restricted to community-dwelling individuals (N=77,669) and nursing-home residents (N=1414), which are shown in supplementary tables. The analyses were performed using IBM SPSS Statistics for Windows, version 24 (IBM Corp., Armonk, N.Y., USA); in all of them, statistical significance was set at p<0.05 (two-tailed).
RESULTS
The study cohort included 79,083 persons (37,626 [47.6%] men and 41,407 [52.4%] women), with a mean age of 65.8 years (SD: 11.3). By age groups, 42,684 (54%) were aged 50-64 years and 10,386 (13.1%) were aged 80 years or older. A total of 77,669 (98.2%) were community-dwelling and 1414 (1.8%) nursing-home residents.
Cohort members were observed for a total of 1,356,358 persons-weeks (1,335,069 for community-dwelling and 21,289 for nursing-home residents). An amount of 4113 (5.2%) of the 79,083 cohort members was tested for SARS-COV-2 infection across study period; 3577 (87%) were negative and 536 (13%) positive (345 in community-dwelling and 191 in nursing-home residents). Two hundred and eighty-three cohort members (280 community-dwelling and 3 nursing-home residents) had a suspected COVID-19 without PCR testing (Table 1).
Across the four-month study period, 576 all-cause deaths were observed (413 in community-dwelling and 163 in nursing-home residents), of which 124 occurred among laboratory-confirmed COVID19 patients and were considered COVID19-related deaths (63 in community-dwelling and 61 in nursing-home residents). Of these 124 deceased patients with laboratory-confirmed COVID19, 112 (90.3%) died by (due to) COVID-19, while 12 (9.7%) died with COVID19 but likely due to other causes according to physician reviewers’ criteria after checking clinical records (ten cases in community dwelling and two cases in nursing-home residents).
All-cause mortality rate was 42.5 deaths per 100,000 persons-week (95% CI: 39.2-46.1) considering the total study cohort. By residence, it was 30.9 per 100,000 persons-week (95% CI: 28.0-34.1) in community-dwelling and 765.7 per 100,000 persons-week (95% CI: 653.9-896.6) in nursing-home residents.
COVID19-related mortality rate was 9.1 per 100,000 persons-week considering the total study cohort, being 4.7 per 100,000 persons-week among community-dwelling and 286.5 per 100,000 persons-week among nursing-home residents. Therefore, deaths from COVID19 represented a 21.5% (124/576) of the overall deaths observed in the total study cohort (15.3% [63/413] in community-dwelling vs 37.4% [61/163] in nursing-home residents).
Of the 576 all-cause deaths, 295 occurred among the 74,687 COVID19-unrelated/ healthy subjects (22.8 deaths per 100,000 persons-week), 150 among the 3577 COVID19-excluded/PCR-negative subjects (236.4 deaths per 100,000 persons-week), 7 among the 283 COVID19-suspected/PCR-unperformed (493.7 deaths per 100,000 persons-week) and 124 among the 536 COVID19-confirmed subjects (4009.1 deaths per 100,000 persons-week) (Table 2).
Table 3 shows univariate analyses reporting all-cause deaths (n=576) by age, sex and pre-existing conditions (comorbidities, chronic medications’ use and vaccinations’ history) in the total study cohort. Table 4 reports it for COVID19-related deaths (n=124). Tables 5 and 6 show Cox regression analyses assessing unadjusted, age&sex-adjusted and multivariable-adjusted risks of suffer all-cause death and COVID19-related death, respectively, in the total study cohort.
Supplementary tables S1-S4 show subgroup analyses restricted to community-dwelling individuals and supplementary tables S5-S8 show subgroup analyses restricted to nursing-home residents.
After multivariable-adjustments assessing risks of all-cause mortality in the total study cohort (see Table 5), increasing age/years (HR: 1.08; 95% CI: 1.07-1.09; p<0.001), male sex (HR: 1.47; 95% CI: 1.24-1.79; p<0.001), nursing-home residence (HR: 1.63; 95% CI: 1.28-2.09; p<0.001), neurological disease (HR: 1.48; 95% CI: 1.18-1.85; p=0.001), cancer (HR: 2.46; 95% CI: 2.03-2.98; p<0.001), cardiac disease (HR: 1.29; 95% CI: 1.06-1.56; p=0.009), liver disease (HR: 1.76; 95% CI: 1.14-2.73; p=0.012), receiving diuretics (HR: 1.55; 95% CI: 1.27-1.88; p<0.001), systemic corticosteroids (HR: 3.49; 95% CI: 2.59-4.70; p<0.001), proton-pump inhibitors (HR: 1.25; 95% CI: 1.03-1.51; p=0.026) and benzodiazepines (HR: 1.34; 95% CI: 1.11-1.61; p=0.003) were associated with a significant increased risk. Conversely, receiving Angiotensin Converting Enzyme inhibitors (ACEIs) (HR: 0.73; 95% CI: 0.60-0.88; p=0.001), Angiotensin II Receptor Blockers (ARBs) (HR: 0.59; 95% CI: 0.46-0.77; p<0.001) and statins (HR: 0.72; 95% CI: 0.59-0.88; p=0.002) were associated with a reduced risk. In subgroup analyses restricted to community-dwelling individuals, the results did not substantially change (see supplementary Table S3). In subgroup analysis restricted to nursing-home residents (see supplementary Table S7), history of tetanus vaccination was associated with a reduced risk (HR: 0.65; 95% CI: 0.44-0.97; p=0.036).
If we consider COVID19-related mortality (see Table 6), many comorbidities (cancer, neurological, renal, respiratory, cardiac disease, diabetes and hypertension) appeared associated with an increased risk in the crude analyses. However, after age&sex-adjustment only neurological, respiratory and cardiac disease were associated with increased risk, and none of them appeared significantly associated with an increased risk after multivariable-adjustment. In this analysis, increasing age/years (HR: 1.08; 95% CI: 1.06-1.10; p<0.001), male sex (HR: 1.82; 95% CI: 1.24-2.70; p=0.003) and nursing-home residence (HR: 12.56; 95% CI: 8.07-19.54; p<0.001) were the strongest predictors for COVID19-related mortality. Considering chronic medications’ use, while many medications were associated with an increased risk in crude analyses, none of them altered the risk of COVID19-related death in multivariable analysis (except for ARBs [HR: 0.52; 95% CI: 0.26-1.01; p=0.054] and oral antidiabetic drugs [HR: 2.14; 95% CI: 1.03-4.43; p=0.042], that were marginally associated with reduced and increased risk, respectively).
In subgroup analyses restricted to community-dwelling individuals evaluating risk of COVID19-related mortality, besides age/years and male sex, taking diuretics (HR: 1.89; 95% CI: 1.03-3.46; p=0.041) and systemic corticosteroids (HR: 3.10; 95% CI: 1.27-7.59; p=0.013) emerged associated with increased risk. None other medication, comorbidity and/or vaccination’s history altered the risk of COVID19-related death after multivariable-adjustments in these community-dwelling individuals (see supplementary Table S4).
In subgroup analyses restricted to nursing-home residents, apart from increasing age and male sex, only receiving antineoplastic agents (HR: 3.81; 95% CI: 1.00-14.44; p=0.049) and, surprisingly, influenza vaccination in prior autumn (HR: 2.34; 95% CI: 1.08-5.07; p=0.031) significantly altered the risk of COVID19-related death after multivariable adjustment (see supplementary Table S8).
DISCUSSION
Direct and indirect mortality due to COVID19 is controversial.[4,5,8,11,13,21,22] This study investigated all-cause and COVID19-related deaths among middle-aged and older adults in the region of Tarragona (Southern Catalonia, Spain) across the first wave of COVID19 epidemic period (March-June 2020), assessing mortality risks by demographic characteristics and underlying medical conditions. Of note, the study was conducted in a region with relatively low incidence of COVID19 as compared with other European and Spanish regions.[2,3,8]
As main findings, all-cause mortality rate across study period was 42.5 deaths per 100,000 persons-week, and COVID19-related deaths represented 21.5% of overall deaths. Considering specifically COVID19-related deaths (those occurred among patients with laboratory-confirmed COVID19), 90.3% died by COVID-19 and 9.7% died with COVID19 (i.e., the death may be attributed to other baseline or concomitant cause.
Increasing age, male sex, and nursing-home residence appear as the strongest predictors for both COVID19-related death and all-cause death. Other pre-existing conditions may also alter mortality risk, especially some comorbidities and chronic medications’ use, which emerged significantly related with all-cause mortality risk.
Importantly, apart from COVID19-related deaths (occurring among patients with laboratory-confirmed COVID19), theoretically COVID19-excluded patients (i.e., those persons with a negative PCR) suffered ten times greater all-cause mortality than unrelated-COVID19/healthy subjects. This fact suggests that, apart from baseline excess risk in PCR-negative subjects (who were older and had more comorbidities than healthy/COVID19-unrelated subjects), some of the observed excess mortality in PCR-negative subjects could be due to undiagnosed/undetected COVID19 despite PCR performed. In this sense, we note that the reliability of PCR testing depends on the quality of the nasopharyngeal samples collected, timing of collection and sensitivity of tests used,[19] which could explain, in part, the excess all-cause mortality observed in PCR-negative subjects in this study.
Regarding excess mortality, the Statistics National Institute (INE) reported an increase in deceases in 2020 (compared to 2019) of 18.95%; which corresponds in absolute numbers to 70,715 deaths, of which only 43,131 were attributed to COVID19.[10] The considerably larger all-cause mortality in PCR-negative subjects observed in the present study may reflect COVID19 underdiagnosis and could explain, in part, some excess all-cause deaths observed across the current COVID19 pandemic.
We note the importance of providing adjusted data instead of crude data. Thus, while crude all-cause mortality rate was almost 25-times greater in nursing-home (766 deaths per 100,000 persons-week) than in community-dwelling individuals (31 deaths per 100,000 persons-week), we found that after multivariable adjustments (by age/sex and pre-existing conditions) nursing-home residence only increased approximately 1.6 times the adjusted risk of all-cause mortality (HR: 1.63; 95% CI: 1.28-2.09) as compared with community-dwelling. In this same way, apart from increasing age and sex male, many pre-existing conditions appeared significantly associated with an increased risk of COVID19-related mortality in the crude analysis (and several in the age&sex-adjusted analysis), but none comorbidity emerged significantly associated with an increased risk after multivariable adjustment.
This highlights the importance of maximizing adjusted data assessing event’s risk in observational studies. In fact, this could explain distinct/opposite data reported in different studies evaluating relationships between some pre-existing comorbidities/conditions and susceptibility/risk of suffering COVID19 infection or death.[4,5,11,13, 21-26] To illustrate it, besides subgroup analyses, we reported all unadjusted, age&sex-adjusted and multivariable-adjusted results here, which is a major strength in this study.
Considering chronic medication, there is scarce data reporting the possible influence of previous use of these drugs on the risk of COVID19-related mortality. We did not observe any significant association in this sense but, interestingly, receiving statins and renin-angiotensin inhibitors were significantly associated with a lower risk of all-cause mortality in our study cohort.
Major strengths in this study were the large size and representativity of the cohort (that included more than 79,000 people, which represents approximately 75% of overall population over 50 years of age in the study area, and the use of survival analysis methods to accurately estimate both risks of all-cause and COVID19-related mortality by distinct important population subgroups (based in its relation with COVID19): COVID19-confirmed subjects, COVID19-suspected subjects (without PCR performed), COVID19-excluded subjects (negative PCR testing) and COVID19-unrelated/healthy subjects). The large size of the study cohort together with the adjustment of major possible confounding variables (e.g., age, type of residence, pre-existing underlying conditions/comorbidities and chronic medications’ use) in the multivariable analyses may provide an acceptable basis to assess mortality risk among the general population ≥50 years across the first wave of COVID19 epidemic period in our setting.
Major limitations in this study are related with its retrospective design and scarce availability of PCR tests during the first weeks of study period. Indeed, considering that PCR testing was not routinely performed for all clinically compatible/suspected COVID19 patients across the study period, the laboratory-confirmed COVID19 cases (and, consequently, the number of COVID19-related deaths) were likely underestimated. Then, all-cause mortality (which is not influenced by the frequency of PCR testing) may be a better measure of COVID19 pandemic impact. As another limitation, the study was conducted in a single geographical area and, logically, specific mortality data may not be directly extrapolated to other geographical regions with distinct epidemic conditions.
The authors recognise these inherent limitations but note that, opposite to many papers reporting only crude COVID19 data, the present study provides age&sex-adjusted and multivariable-adjusted data evaluating both all-cause and COVID19-related mortality risks. Importantly, the estimations may considerably vary depending on type of analyses/adjustments performed; therefore, we underline again the importance of maximizing adjustments. We did subgroup analysis (nursing-home/community-dwelling) and multivariable-adjustments, but, as in all observational studies, a residual confounding due to unmeasured factors (e.g., socio-economical, lifestyle, job-or healthcare-related factors) cannot be excluded.
CONCLUSIONS
In summary, in this large cohort study including 79,083 middle-aged and older adults followed across the first wave of COVID19 epidemic period in the region of Tarragona, we found that deaths from COVID19 represented 21.5% of all-cause mortality occurred across study period (March-June 2020). Increasing age, sex male, nursing-home residence, cancer, neurologic, cardiac or liver disease, receiving diuretics, systemic corticosteroids, proton-pump inhibitors and benzodiazepines were associated with an increased risk of all-cause mortality; conversely, receiving renin-angiotensin inhibitors and statins were associated with a reduced risk. Age/years, sex male and nursing-home residence were strong predictors for COVID19-related mortality, but no comorbidity was independently associated with increased risk.
Interestingly, apart from COVID19-related deaths (approximately 10% of them could be attributed to another baseline or concomitant cause), theoretically COVID19-excluded patients (PCR-negative) suffered considerable greater all-cause mortality than healthy/unrelated-COVID19 subjects. This data suggests COVID19 underdiagnosis and it could explain, in part, the unexpected excess deaths observed during the COVID19 pandemic.
Data Availability
These data have been obtained from the Catalonian Health Institute Information System. Interested authors might obtain CAPAMIS data (previous ethics and scientific approval by the ethics and clinical research committee of the Primary Care Research Institute Jordi Gol (IDIAP Jordi Gol)) addressing purposes to the Institution.
CONTRIBUTORS
AVC designed the study; AVC and ESG assessed outcomes and wrote the manuscript; ESG, CDC and MFP obtained data; ESG and AVR did statistical analyses and edited the manuscript; OOG revised the final version; AVC coordinated the study.
FUNDING
This study is supported by a grant from the Instituto de Salud Carlos III of the Spanish Health Ministry (file COV20/00852; call for the SARS-COV-2/COVID-19 disease, RDL 8/2020, March 17, 2020). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
COMPETING INTEREST
All authors, none declared.