Introduction

Coronavirus disease 2019 (COVID-19) is a global outbreak of a new respiratory infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). According to the last WHO report as of 18 August 2021, there have been over 207 million confirmed cases of COVID-19 worldwide, including over four million deaths [1].

The presence of comorbidities among hospitalised COVID-19 patients [2, 3], particularly in the elderly, are linked to acute respiratory syndrome, poor clinical outcomes and increased mortality rate [4, 5]. These patients often present factors that may increase the risk of bone fragility, such as systemic inflammation, older age, glucocorticoid treatment and immobilisation. Vertebral fractures (VFs) are the most frequent fragility fracture caused by osteoporosis and commonly occur without a recognisable trauma [6]. Among women over 50 years old, the incidence of VFs increases with age following the same trend as the osteoporosis prevalence, with 25% of women older than 80 having experienced at least one VF [7]. VFs are associated with reduced pulmonary function [8,9,10], disability [11] and an over 60% increased risk of mortality [12, 13]. Recently, Di Filippo et al. retrospectively assessed the VFs prevalence and clinical impact in COVID-19 patients in a tertiary health care hospital in Italy. In their study, including 114 patients with available lateral chest X-rays, the authors found that 36% of COVID-19 patients had a thoracic VF [14]. Fracture severity was associated with COVID-19 mortality, implying a potential role of bone fragility in disease severity. Moreover, Kottlors et al. found that low bone mineral density (BMD) measured by computer tomography (CT) scan was a risk factor for intensive care unit (ICU) admission in a cohort of 58 COVID-19 patients, although BMD analysis did not result in a prognostic advantage over simply considering age [15]. These data suggest that VF could be even more prevalent than expected based on large cohort studies, where VF prevalence is between 8 and 25% [16]. However, it is unknown whether VF prevalence is higher compared to a background population of non-COVID-19 subjects. It is also unknown whether the effect of VF on mortality is greater in COVID-19 compared to those without COVID-19. Therefore, the aim of this study was to evaluate VF and their effect on survival in COVID-19 and non-COVID-19 patients. To this aim, we analysed a large cohort of patients consecutively admitted to the emergency department (ED) for clinical suspicion of symptomatic SARS-CoV-2 infection. We hypothesised that COVID-19 patients with VFs have lower survival compared to those without COVID-19.

Methods

Study design and population

This was a single centre retrospective cohort study including patients admitted to the ED of Trauma Center Public Hospital Bufalini, Cesena, Italy, during the first pandemic wave, who were subjected to CT scan for detection of pulmonary COVID-19 when clinical suspicion was present based on symptoms and clinical findings. The first case of COVID-19 (patient 0) at this Hospital was registered on February 26, 2020. This report includes patients admitted to the ED until April 28, 2020 for clinical suspicion of SARS-CoV-2 infection. During this time frame, CT scans were performed on a total of 501 patients [age 64.36 (18.71) years, 277 males], together with clinical evaluation and real-time reverse-transcriptase polymerase chain reaction (RT-PCR) from a nasal and/or throat swab. A diagnosis of confirmed COVID-19 was based on a SARS-CoV-2 positive RT-PCR together with signs, symptoms and radiological findings suggestive of COVID-19 pneumonia. Patients evaluated in the ED for clinical suspicion of SARS-CoV-2 infection who tested negative by RT-PCR nasopharyngeal swab served as control group (non-COVID-19). The primary exposure and outcome measures were prevalent VFs at hospital admission and death at 30 and 120 days from hospital admission, respectively. For all patients, we obtained age, gender and admission to ICU during hospital stay. In addition, the following clinical and biochemical features were collected in COVID-19 patients: leucocyte count, C-reactive protein (CRP), serum creatinine and estimated glomerular filtration rate (using the CKD-EPI equation), plasma lactic dehydrogenase, hospital length of stay. Comorbidities were identified in medical history by patients or caregivers and included: history of cardiovascular disease (any disease of the cardiovascular system including hypertension), diabetes and chronic obstructive pulmonary disease. Visceral and subcutaneous adipose tissues were assessed as previously described [3].

Vertebral fractures and bone density assessment

VFs were assessed using sagittal image reconstruction from chest CT scanner (Philips Diamond Select Brillance CT 64-slice) standard acquired covering all thorax volume from D1 to a plane transverse to L2. To define osteoporotic VFs, we used the visual semiquantitative method proposed by Genant et al. [17] by two radiologists in consensus (S.B. and M.V.). Vertebral deformities were defined as follows: (1) wedge deformity: marked reduction of the vertebral anterior height without evidence of bone discontinuity; anterior height reduced by ≥4 mm compared with posterior height; (2) endplate deformity: central height of the body reduced by ≥4 mm with respect to the posterior and anterior walls; (3) compression deformity: all three heights reduced by ≥4 mm compared with adjacent vertebrae. Spinal trabecular bone density was assessed on CT scans using Hounsfield unit (HU) quantification [18]. HUs were measured at the sagittal cross sections of the trabecular regions of L1 vertebral body, avoiding cortex value. The vertebral body was divided into three axial segments and HUs were calculated by placing a circular region of interest over an area of trabecular bone on the vertebral body. D12 or any cranial vertebral body HU value were used when L1 was subjected to VF. CT-attenuation values ≤110 HU (high specificity cut-off) were used for definition of low bone density [19, 20].

Statistical analysis

Patients’ characteristics were described using means and standard deviations or medians and interquartile ranges, as appropriate, and percentages. Student’s t test or the Wilcoxon rank-sum test were used as appropriate to compare variables between groups. When variables were not normally distributed, we used a logarithmic transformation. Categorical variables between groups were compared by the chi-square test. The risk of mortality was estimated using the Kaplan–Meier method with log‐rank test for groups comparison. All statistical tests were two-tailed. Cox regression models were used to estimate the hazard ratios (HR) of death by VFs. Statistical analyses were performed in R Statistical Software 3.3 (R Foundation for Statistical Computing, Wien, Austria).

Results

Clinical features of COVID-19 and non-COVID-19 patients

Two hundred thirty-nine patients had a confirmed diagnosis of COVID-19. Of those, 188 patients (78.7%) had CT findings suggestive of COVID-19 pneumonia. Fifty-eight (24.2%) COVID-19 patients required ICU admission (ICU-COVID-19 group). The non-COVID-19 group consisted of 262 patients (Table 1). There were no differences among COVID-19 and controls in terms of age [mean (standard deviation): 63.7 (18.1) vs. 64.9 (19.3) years, respectively; p = 0.468], gender or body mass index, although a male preponderance was observed in COVID-19 patients (58.2% vs. 52.7%, respectively, p = 0.217) (Table 1). The 30-day death rate was 15.8% (38/239) and 12.2% (32/262) in the COVID-19 and non-COVID-19 groups, respectively (p = 0.234). At the longest follow-up (120 days), number of deaths increased slightly further up to 21.8% (52/239) and 17.6% (46/262) in the two groups, respectively (p = 0.236).

Table 1 Clinical and biochemical features of the studied population

Prevalence of vertebral fractures

Prevalence of VFs was similar between non-COVID-19 and COVID-19 groups [50/262 (19%) vs. 53/239 (22.2%); p = 0.458] (Table 1). A similar prevalence was found also in COVID-19 patients with CT findings of pneumonia [36/188 (19.1%)]. VFs were not associated with ICU admission (Table 2). Multiple VFs tended to be more frequent in ICU-COVID-19 than in nICU-COVID-19, but this difference was not significant (57% vs. 38% of total fracture number, respectively; p = 0.34). Consistently, spinal trabecular density was lower in fractured patients compared to those without fractures, in both COVID-19 [71.41(38.56) vs. 134.3(57.12) HU; p < 0.001; Table 2] and controls [77.03(36.99) vs. 131.4(53.64) HU; p = 0.008]. Using a high specificity cut-off (HU ≤ 110) osteoporosis was detected in 45% COVID-19 vs. 47% control subjects (p = ns). A more sensitive cut-off (HU ≤ 135) [19, 20] identified osteoporosis in 62 COVID-19 vs. 64% control subjects (p = ns). Table 2 shows differences in clinical and biochemical features of COVID-19 patients according to VFs. Patients with VFs were older than those without VFs [76.4(10.8) vs. 60.3(18); p < 0.0001] and more commonly diagnosed with cardiovascular comorbidities. Gender, BMI, TC-derived adiposity measures (visceral and subcutaneous adipose tissue), inflammatory markers (CRP, leucocyte counts) and length of in-hospital stay were similar between patients with and without VFs. Although diabetes tended to be more prevalent in the VFs group, this difference was not significant (30.2 vs. 17.3%, p = 0.065).

Table 2 Clinical and biochemical features of COVID-19 subjects with or without vertebral fractures (VFs)

Mortality risk in COVID-19 and non-COVID-19 patients according to VFs

Mortality was associated with VFs in both COVID-19 and non-COVID-19 control groups. COVID-19 patients without VF had a 12.8% (95% CI: 7.9–17.5%) death probability within 30 days, which increased progressively with the number of VFs (20.7%, 95% CI 4.5–34.1%, in those with a single VF, and 34.8%, 95% CI: 12.1–51.6%, in those with multiple VFs) (Fig. 1). Non-COVID-19 patients without VF had a 6.6% (95% CI: 3.2–9.9%) death probability within 30 days from ED evaluation, which increased over five fold in fractured patients regardless of fracture number (40%, 95% CI: 14.2–58.0% in patients with a single VF; 33.3%, 95% CI: 14.1–48.2% in those with multiple VFs) (Fig. 2). Table 3 shows mortality HR in both COVID-19 and non-COVID-19 groups according to VFs. In the unadjusted models, the risk for 30 days mortality was higher in COVID-19 patients, at 30 and 120 days, but statistical significance was reached only in those with multiple VFs (30-day HR 3.03, 95% CI 1.36–6.75; and 120-day HR 2.91, 95% CI 1.43–5.91, respectively, Table 3). Risk of mortality was significantly higher in non-COVID-19 patients either with one (HR 7.46, 95% CI 3.12–17.8, p < 0.001) or multiple VFs (HR 6.2, 95% CI 2.75–13.98) compared to patients without VFs at 30 and 120 days (Table 3). After adjustment for age, sex and spinal trabecular bone density, risk of death remained associated with VFs in non-COVID-19 (30-day HR for a single VF: 3.94, 95% CI 1.57–9.87; 30-day HR for multiple VFs 3.17, 95% CI 1.30–7.74), but not in the COVID-19 group (30-day HR in patients with a single VF: 0.82, 95% CI 0.31–2.17; 30-day HR in patients with multiple VFs 1.58, 95% CI 0.45–5.57). A similar trend was shown in both groups at day 120 (Table 3). When restricting analysis to patients with COVID-19 pneumonia detected by CT, fractures were not significantly associated with mortality at either day 30 or day 120 (Table 4).

Fig. 1
figure 1

Survival probability in subjects with COVID-19 according to vertebral fractures. Survival probability was lower in subjects with fractures and decreased faster in subjects with multiple vertebral fractures compared to those without fractures

Fig. 2
figure 2

Survival probability in subjects without COVID-19 according to vertebral fractures. Survival probability decreased faster in subjects with vertebral fractures compared with those without fractures, regardless the number of fractures

Table 3 Risk of death at days 30 and 120 in both COVID-19 and non-COVID-19 groups according to vertebral fractures
Table 4 Risk of death in subjects with COVID-19 pneumonia according to vertebral fractures

Discussion

In this study, we found that 22% of COVID-19 patients have at least one CT-detectable VF. VFs were associated with mortality although this effect was not significant in COVID-19 patients after adjusting for age and sex. By contrast, VFs strongly increased mortality risk in non-COVID-19 patients in the short term after hospitalisation.

In agreement with our data, Di Filippo et al. did not find a significant association between VFs prevalence and mortality, although survival was lower in those with severe VFs in comparison with moderate or mild VFs [14]. Compared to Di Filippo et al., we observed a lower prevalence of fractures in COVID-19 patients (22.2%), which was similar to the non-COVID-19 control group (19%). Our finding is in line with data from the European population [16] where 18–26% of women older than 50 years, and 8–23% of men aged 50–75 years had a morphometric fracture. Differences between our data and those by De Filippo may be explained by study design, patient’s selection and different clinical conditions.

Our findings are also consistent with previous literature that identified an association between VFs and increased mortality [13, 21]. In a prospective Australian cohort study mortality rates were higher for men than women during the 18-year follow-up period and with greater risk soon after the fracture event [21]. Another prospective cohort study has shown a 1.6 relative risk for mortality in post-menopausal women followed for an average of 2.9 years [12]. Besides, a retrospective study of Swedish elderly women followed for an average of 5 years following a VF found that risk of mortality was increased immediately after the fracture event and decreased by 16% per annum [13]. However, such studies focused mainly on long-term mortality risk after a VF. To the best of our knowledge, our study is the first to assess the association between prevalent VFs and short-term mortality. In our non-COVID-19 cohort, VF determined a spike in mortality risk (up to seven times within 30 days from admission), which tended to plateau over the total observation timeframe of 6 months. We may speculate that VF may exacerbate mortality risk in presence of weakened clinical condition. Therefore, our data expand previous findings highlighting the importance of VFs as strong predictors of short-term mortality in hospitalised patients.

The biological explanation on the difference in risk of mortality related to VFs among COVID-19 and non-COVID-19 patients still needs to be clarified. We may speculate that COVID-19 patients carry a background mortality risk higher than that simply conferred by VFs. In this context, age and gender have been consistently associated with COVID-19 severity and fatality [22,23,24] as they are among the main factors shaping immune response and susceptibility to infection [25, 26]. Thus, while gender and age are main drivers of COVID-19 fatality, our data do not support the role of VFs as independent cause of death in this population. In contrast, given the potential heterogeneity of the non-COVID-19 group, we were unable to account for potential confounding factors other than age and gender. Lack of clinical information on the non-COVID-19 patients does not allow us to determine a causal relationship between VFs and mortality. However, soon after the spread of the COVID-19 pandemic and social lockdown restrictions, health care systems have been re-organised with detrimental effects on the care of non-COVID-19 patients [27]. Accordingly, a 30% reduction in admissions for acute cardiovascular events [28,29,30,31,32,33,34,35,36,37] and an increased rate of deaths not related to COVID-19 during pandemic have been described [36,37,38,39,40,41,42]. We may speculate that by refraining from seeking medical attention, some patients reached the hospital later and in worse clinical condition during the pandemic. VFs, despite the barriers imposed by the pandemic [43], should not be neglected from healthcare providers even during hospitalisation. Increasing the risk of future fragility fractures and mortality, VFs should be carefully evaluated and treated, as indicated by recent ASBMR guidelines [44].

Our study has limitations and strengths. In the non-COVID-19 patients we were not able to obtain clinical information related to admission in ED, cause of death or possible clinical complications. Despite these limitations, both controls and COVID-19 patients were selected based on the same criterion (clinical suspicion for COVID-19), from the same geographic area and timeframe of the first COVID-19 pandemic wave. The total lumbar spine scan was not available and therefore VFs prevalence might be underestimated. We did not have information on a previous diagnosis of osteoporosis, recency of VF or anti-osteoporotic treatments, but differences between COVID-19 and the non-COVID-19 patients were unlikely. Indeed, main factors associated with osteoporosis such as age and CT-derived bone density were similar between the two groups. Furthermore, our data may not be generalised to other COVID-19 cohorts as the clinical algorithm implemented for ICU admission and hospital/ICU capacity may vary according to different clinical settings. Strengths of our approach include: the use of CT scans which have higher accuracy over plain X-rays; the systematic assessment of VFs in a large cohort of unselected COVID-19 patients consecutively admitted to the ED, which minimised selection biases; the inclusion of a large control group of non-COVID-19 patients from same geographic area and timeframe as COVID-19 patients.

In conclusion, VF did not have an independent effect on mortality in COVID-19 patients. They strongly increased mortality risk in non-COVID-19 patients. However, VF remain often neglected and associated with worse outcome. We advocate prompt recognition and treatment with vitamin D and anti-osteoporosis drugs as soon as fragility fractures are identified in both COVID-19 and non-COVID-19 patients [45].