Introduction

While a great deal of information regarding the evolution of patients with COVID-19 infection [1,2,3,4,5,6,7,8,9,10] has been recorded in hospitals over the world since the onset of the novel coronavirus pandemic, only a few data have been published up to date regarding the course of the disease in patients admitted in long-term care facilities [11].

Studies by McMichael T.M. and Aron MM [11, 12] from a long-term care facility suggest that nursing home populations (e.g., older adults often with underlying chronic medical conditions) could be especially vulnerable to SARS-CoV-2 infection and at higher risk for COVID-19-associated morbidity and mortality.

In the current study, we aimed to provide information about the clinical evolution of older patients infected by COVID-19 admitted in nursing homes and to identify prognostic mortality factors in such a population.

The study was approved by the Research Ethics Committee of the Andorran Healthcare Services.

Methods

Study design

The study was performed at a public nursing home from March 15 to June 5, 2020, identifying all the confirmed COVID-19 patients. This reference nursing home was transformed into a COVID-19 center for the admission of institutionalized seniors considered probable or confirmed cases of COVID-19 from the same center or other nursing homes. To this end, a joint action strategy was settled with the other nursing homes in the country, establishing a referral circuit to support them in the event of an outbreak of COVID-19 disease in their institution. Within 72 h of detecting the first COVID case, 35 residents were moved to a hotel to prevent further spread following medical criteria (residents with a low level of functional and cognitive dependence were transferred). In the different serological and PCR controls carried out, only three residents housed in the hotel had COVID-19 infection, and the mortality in this population was null.

Data collection

The following data were recorded for all patients: age, sex, date of admission, length of stay, the origin of the patient, dementia stage, Barthel index, Charlson comorbidity index, previous flu vaccination, clinical presentation, laboratory results, treatment, Reverse Transcription Polymerase Chain Reaction (RT-PCR) for SARS-CoV-2, hospital referral, mortality rate, and case fatality rate.

Definitions

Four categories of the patient’s origin were defined: the reference nursing home, other nursing homes, home, or hotel. Patients from home belonged to a program of day care or home support for frail patients deployed by the same long-term care center.

Dementia was measured by the Global Deterioration Scale (GDS), developed by Dr. Reisberg [13], which consists of seven stages. 1: no cognitive decline, 2: very mild cognitive decline, 3: mild cognitive decline, 4: moderate cognitive decline, 5: moderately severe cognitive decline, 6: moderately severe decline, 7: very severe cognitive decline.

Barthel index [14] is an ordinal scale of functional capacity used to measure performance in daily living activities, with values ranging from 0 (totally dependent) to 100 (totally independent). Proposed guidelines for interpreting Barthel index are: 0–15: totally dependent, 20–35: very dependent, 40–55: partially dependent, 60–75: minimally dependent, 80–100: totally independent [15].

Underlying diseases were considered as the presence of comorbid illness with the age adjusted Charlson comorbidity index [16], which predicts 10-year survival in patients with multiple comorbidities.

Patients were grouped into three clinical categories: asymptomatic, respiratory symptoms (rhinitis, pharyngitis, cough, expectoration, and dyspnea) and digestive symptoms (diarrhea). Fever was defined as an axillary temperature of at least 37.5 °C.

An analytical control was carried out, at least, on the admission and the discharge, recording the presence of lymphocytopenia, anemia, and or thrombocytopenia. Inflammation markers such as C-reactive protein (CRP), Ferritin, Dimer D, Lactate dehydrogenase (LDH), and troponin were also collected. Lymphocytopenia was defined as a total lymphocyte count of less than 1.0 × 109/L (1000/μL) and thrombocytopenia as a platelet count of less than 150 × 103 per μL. The diagnosis of anemia in men was based on a hemoglobin of less than 13–14 g/dL; in women, less than 12–13 g/dL.

As for the treatment, five categories were defined: patients who received hydroxychloroquine and azithromycin, hydroxychloroquine only, hydroxychloroquine plus another antibiotic, beta-lactam, or quinolone antibiotics, and no treatment.

We also recorded the percentage of patients with two consecutive (48 h apart) negative molecular detection results for SARS-CoV-2 from a nasopharyngeal swab after 14 days of treatment with hydroxychloroquine and azithromycin.

Patients whose symptoms resolved and who had two consecutive (48 h apart) RT-PCR for SARS-CoV-2 negative were considered successfully treated and cured [17].

Four categories of discharge destination were defined: reference nursing home in non-COVID areas, other nursing homes, home, or death.

The mortality rate associated with COVID-19 in the center was defined as the death rate in the reference nursing home population.

The case fatality rate in the reference nursing center is the proportion of deaths from a COVID-19 disease compared to the total number of people diagnosed with the disease admitted in the center.

Statistical analysis

Statistical analyses have been carried out to identify risk factors for mortality in patients with COVID 19. Bivariate tests have been performed between the qualitative and quantitative variables and the successfully treated variable COVID (yes/no). A chi-squared test was used for categorical variables and the t test/Mann–Whitney test for quantitative variables. Additional bivariate analyses were performed to compare subject characteristics between treatment groups, to identify possible confounders. Variables with a p value lower than 0.05 were included in a multivariate logistic regression model [18] to identify independent predictors of mortality. The final model was obtained after removing all non-statistically significant variables (backward selection procedure). Furthermore, a second multivariate regression model was adjusted, including only baseline variables to obtain a score to predict mortality independently of the effect of pharmacological treatment. ROC curves have been obtained to evaluate the fit of the models, and the area under the curve (AUC) has been calculated. The best cut-off point has been determined, and sensitivity and specified values have been obtained. Statistical analyses were performed with the SAS system version 9.4 (SAS Institute Inc., Cary, North Carolina, USA). The statistical significance level was set at 0.05.

Data availability

The datasets generated and/or analysed for this study during the current study are available from the corresponding author upon request on reasonable request.

Results

The study population comprised 100 confirmed COVID-19 cases with a mean age of 85 (IQR 65–103) years. Sixty-two percent were female, and the average length of stay was 22 days. Fifty-two patients were initially from the reference nursing home, 36 came from other nursing homes, nine from home, and three comings back from the hotel.

According to the global deterioration scale for the assessment of primary degenerative dementia, 73% of patients presented dementia, and 91.6% of dead had severe dementia (GDS6–GDS7). Seventy-six percent of patients had some level of functional dependence (Fig. 1). Clinically, 57% had respiratory symptoms, and 39% overcome the infection without any symptoms. Blood tests revealed anemia in 36 cases and lymphocytopenia in 38 patients.

Fig. 1
figure 1

Functional capacity and percentage of COVID-19 patients cured

Eighty-three percent received pharmacological treatment, mostly with hydroxychloroquine and azithromycin (70%) (Fig. 2). Only five patients had diarrhea as a side effect related to hydroxychloroquine. Cardiac monitoring was performed by electrocardiogram, and no rhythm changes were observed with this treatment in any patient.

Fig. 2
figure 2

Type of treatment and percentage of COVID-19 patients cured. H + A hydroxychloroquine and azithromycin. H hydroxychloroquine. No treatment includes Others: Beta-lactam or Quinolone antibiotics

After 14 days, 12 patients had 2 negative molecular detection results for SARS-CoV-2 from a nasopharyngeal swab. Seventy-six patients were considered cured, and 24 died, of which 4 happened in the reference hospital. Thirty-five patients were discharged and sent back to non-COVID areas of Cedre’s nursing home, 31 to other nursing homes, and 9 to home.

The mortality rate associated with COVID-19 only for the original residents of the reference nursing center was 14% (18 deaths/125 residents). The overall case fatality rate at the reference nursing home was 20% (20 deaths/100 COVID-19 residents).

In the bivariate analysis, the following factors showed a significantly greater risk of mortality among COVID-19 patients: male-gender, fever, respiratory symptoms, pharmacological treatment, type of treatment, serum therapy, oxygen therapy, dementia, Barthel index, lymphocytopenia, LDH, and D-dimer (Table 1).

Table 1 Risk factors associated with COVID-19 mortality on bivariate analysis

No statistically significant differences were observed in Charlson’s index, treatment started within 24 h, anemia, low platelet count, ferritin, CRP, troponin levels, or previous flu vaccination.

Table 2 shows the results of the bivariate analysis results between treatment groups. Only ‘treatment started within 24 h’ was found statistically imbalanced between groups, but we did not find any relation between this variable and mortality. Temperature, serum therapy and D-dimer presented some differences between groups, although not statistically significant.

Table 2 Bivariate analysis with treatment

In the multivariate regression analysis, the independent risk factors associated with a higher COVID-19 related mortality were: male gender, type of treatment, Barthel index, and lymphocytopenia (Table 3). Although statistically significant differences were observed between LDH, D-dimer and mortality values in the bivariate analysis, (as their levels increase, the risk of mortality increases) they were not included in the multivariate logistic regression due to the high number (20%) of missings in these two parameters. A second multivariate analysis was performed including ‘treatment started within 24 h’, temperature and serum therapy, and similar results as in the initial multivariate analysis were obtained.

Table 3 Risk factors associated with COVID-19 mortality on multivariate analysis

To obtain a score that allows predicting COVID-19 progression from baseline, a new model was adjusted, excluding pharmacological treatment.

To calculate a prognostic score, the following formula was used:

Score = 1.4 × (Sex = “Male”) – 0.04 × Barthel + 0.9 × (Lymphopenia = “Yes”)

The higher the score, the higher was the probability of dying. We obtained an AUC of 0.85 (Fig. 3). The cut-off point that best ranks between patients who are successfully treated and those who are not cured is 0.2, meaning all those patients with a score higher than 0.2 are most at risk of dying, with a specificity of 79.7% and a sensitivity of 80.0%.

Fig. 3
figure 3

Prognostic COVID-19 mortality score-ROC curve for model

Discussion

With the outbreak of COVID-19, the Andorran Healthcare Service acted quickly and proactively to control the pandemic in nursing homes. The early reaction transforming the reference nursing home into an intermediate healthcare facility and moving to a hotel 35 healthy residents from this center with a low level of functional and cognitive dependence, probably avoided excess mortality since none of the people transferred to the hotel died. One hundred COVID-19 older patients were admitted to the reference intermediate care nursing home, and they received the same care as they would get in hospital wards, 76% of cases were cured. The case fatality rate in those remaining in the nursing home was lower (20%) than that reported in other settings (34%) [11, 12] despite being mostly patients with high functional dependence. Seven patients were transferred to a tertiary care hospital for presenting medical criteria for mechanical respiratory ventilation.

As far as we know, this is the first reported study describing how a nursing home was transformed into an intermediate care facility to face the COVID-19 outbreak and analyzing which prognostic factors could predict infection-related mortality in this segment of the older population, despite the deployment of a significant amount of medical resources.

Several factors related to gender, respiratory symptoms, supportive and specific therapy, cognitive and functional deterioration, and inflammatory and immune factors were associated with mortality in the bivariate analysis. In the multivariate logistic regression analysis, only male gender, Barthel index, lymphocytopenia, and hydroxychloroquine plus azithromycin were identified as independent factors for mortality. Variables that presented more differences between treatment groups were included in a second multivariate model to evaluate potential confounding factors, and similar results were obtained.

Many reports have outlined that, although rates of infection are similar among men and women, men are at higher risk of death from the COVID-19 infection [3, 19]. The explanation for that phenomenon is unclear, although previous studies have suggested that women are less susceptible to viral infection, possibly because of the protection of X chromosome and sex hormones, which play an essential role in innate and adaptive immunity [20].

The institutionalized older persons present a cognitive and functional decline that may worsen with acute illness such as COVID-19 [21]. This study reveals a direct relationship between decreased functional capacity measured by the Barthel index and mortality. Functional status in older people has been identified as a prognostic mortality factor in respiratory infections [22]. The current study is the first to prove a direct relationship between the Barthel index and mortality in COVID-19 infected patients.

Contrary to what was expected, we did not find any relationship between clinical comorbidity and mortality in these patients, a factor that is usually present in studies evaluating mortality for other causes in older people [23,24,25], and it was also relevant in previous reports of COVID-19 infection in China [5, 6, 26]. Moreover, flu vaccination did not have any impact on mortality either.

The other independent risk factor detected was lymphocytopenia, a common finding in COVID-19 infected patients that has been previously related to mortality [4, 27, 28]. This could be due to a direct effect of the virus or the consequence of the cytokine-mediated inflammatory cascade inducing lymphocytes migration [27]. In any case, T-lymphocytes deficiency or dysregulation may reflect oversized immune reactions that can contribute to disease severity and mortality [29].

Regarding the treatment, hydroxychloroquine and azithromycin were prescribed for 5 days as the initial treatment for all patients. Before and during this initial treatment period, electrocardiograms were performed to monitor patient’s tolerance to the treatment. When this pharmacological combination was contraindicated based on electrocardiographic alterations, either hydroxychloroquine alone or beta-lactams were prescribed instead.

It is worth highlighting that the multivariate logistic regression analysis identified hydroxychloroquine plus azithromycin treatment as an independent factor favoring survival compared with no treatment or other treatments. At present, contradictory results have been published regarding the effectiveness of this treatment in COVID-19-infected patients and the associated risk of adverse events using those medications [30,31,32,33,34]. It might be possible that this treatment could only work when used early and in a particular group of patients with similar specific characteristics as in the present study. In any case, only randomized clinical trials could definitively clear these uncertainties.

Using only those independent risk factors identified in the multivariate analysis, which could be easily collected at onset of the disease infection, we have elaborated a score that, beyond the potential beneficial effect of pharmacological treatment, may be able to predict mortality with reasonable specificity and sensitivity.

This study has some limitations. First, it is a retrospective study, and some missing laboratory data precluded to include LDH and D-dimer parameters in the multivariate regression analysis. Second, the advanced age of all patients included has probably been prevented from identifying this parameter as a prognostic factor for mortality. Third, the high prevalence of dementia and functional dependence among the patients included may have generated results that might not be precisely reproducible in other cohorts of older patients with better functional capabilities. Finally, it would be interesting to evaluate the presence of interactions between variables for which a bigger sample size would be required.

Conclusion

In summary, this study describes the feasibility of transforming a nursing home into an intermediate healthcare facility to face an epidemic outbreak of COVID-19 in this setting and limiting referrals to a tertiary care hospital. To take care of COVID-19-infected older person, the independent prognostic factors identified in the present study can be of enormous help to adjust the adequate healthcare resource provision and to limitate therapeutic efforts in case of new outbreaks of the pandemic.