Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Predictors of mortality among hospitalized COVID-19 patients and risk score formulation for prioritizing tertiary care—An experience from South India

  • Narendran Gopalan ,

    Contributed equally to this work with: Narendran Gopalan, Sumathi Senthil

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft

    gopalannaren@yahoo.co.in

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Sumathi Senthil ,

    Contributed equally to this work with: Narendran Gopalan, Sumathi Senthil

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Resources, Writing – review & editing

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Narmadha Lakshmi Prabakar ,

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Resources, Writing – review & editing

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Thirumaran Senguttuvan ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Adhin Bhaskar,

    Roles Data curation, Formal analysis, Methodology, Software, Validation, Writing – review & editing

    Affiliation Department of Statistics, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Muthukumaran Jagannathan,

    Roles Methodology, Resources, Writing – review & editing

    Affiliation Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Ravi Sivaraman,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation MDRU, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Jayalakshmi Ramasamy,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Ponnuraja Chinnaiyan,

    Roles Data curation, Validation, Writing – original draft

    Affiliation Department of Statistics, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Vijayalakshmi Arumugam,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of Microbiology, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Banumathy Getrude,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of Community Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Gautham Sakthivel,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Vignes Anand Srinivasalu,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Dhanalakshmi Rajendran,

    Roles Investigation, Resources, Writing – review & editing

    Affiliations Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India, Division of Epidemiology and Operational Research, ICMR-Vector Control Research Centre, Puducherry, India

  • Arunjith Nadukkandiyil,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Vaishnavi Ravi,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Sadiqa Nasreen Hifzour Rahamane,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Nirmal Athur Paramasivam,

    Roles Investigation, Resources, Writing – review & editing

    Affiliation Department of General Medicine, Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India

  • Tamizhselvan Manoharan,

    Roles Formal analysis, Validation, Writing – review & editing

    Affiliation Department of Statistics, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Maheshwari Theyagarajan,

    Roles Data curation, Resources, Writing – review & editing

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Vineet Kumar Chadha,

    Roles Formal analysis, Methodology, Validation, Writing – review & editing

    Affiliation Central Leprosy Teaching & Research Institute, Chengalpattu, Tamil Nadu, India

  • Mohan Natrajan,

    Roles Conceptualization, Writing – original draft

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Baskaran Dhanaraj,

    Roles Conceptualization, Writing – original draft

    Affiliation Department of Clinical Research, ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • Manoj Vasant Murhekar ,

    Roles Conceptualization, Methodology, Project administration, Writing – original draft, Writing – review & editing

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliations ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India, ICMR-National Institute of Epidemiology, Chennai, Tamil Nadu, India

  • Shanthi Malar Ramalingam ,

    Roles Methodology, Project administration, Supervision, Writing – original draft

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliations Government Chengalpattu Medical College & Hospital, Chengalpattu, Tamil Nadu, India, Government Kilpauk Medical College, Chennai, Tamil Nadu, India

  •  [ ... ],
  • Padmapriyadarsini Chandrasekaran

    Roles Project administration, Supervision, Writing – original draft

    ‡ These authors also contributed equally to this work. NG and SS are Joint First Authors. NLP and TS are Joint Second Authors. MVM, SMR and PC are Joint Senior Authors

    Affiliation ICMR-National Institute for Research in Tuberculosis (Formerly Tuberculosis Research Centre), Chennai, Tamil Nadu, India

  • [ view all ]
  • [ view less ]

Abstract

Background

We retrospectively data-mined the case records of Reverse Transcription Polymerase Chain Reaction (RT-PCR) confirmed COVID-19 patients hospitalized to a tertiary care centre to derive mortality predictors and formulate a risk score, for prioritizing admission.

Methods and findings

Data on clinical manifestations, comorbidities, vital signs, and basic lab investigations collected as part of routine medical management at admission to a COVID-19 tertiary care centre in Chengalpattu, South India between May and November 2020 were retrospectively analysed to ascertain predictors of mortality in the univariate analysis using their relative difference in distribution among ‘survivors’ and ‘non-survivors’. The regression coefficients of those factors remaining significant in the multivariable logistic regression were utilised for risk score formulation and validated in 1000 bootstrap datasets.

Among 746 COVID-19 patients hospitalised [487 “survivors” and 259 “non-survivors” (deaths)], there was a slight male predilection [62.5%, (466/746)], with a higher mortality rate observed among 40–70 years age group [59.1%, (441/746)] and highest among diabetic patients with elevated urea levels [65.4% (68/104)]. The adjusted odds ratios of factors [OR (95% CI)] significant in the multivariable logistic regression were SaO2<95%; 2.96 (1.71–5.18), Urea ≥50 mg/dl: 4.51 (2.59–7.97), Neutrophil-lymphocytic ratio (NLR) >3; 3.01 (1.61–5.83), Age ≥50 years;2.52 (1.45–4.43), Pulse Rate ≥100/min: 2.02 (1.19–3.47) and coexisting Diabetes Mellitus; 1.73 (1.02–2.95) with hypertension and gender not retaining their significance. The individual risk scores for SaO2<95–11, Urea ≥50 mg/dl-15, NLR >3–11, Age ≥50 years-9, Pulse Rate ≥100/min-7 and coexisting diabetes mellitus-6, acronymed collectively as ‘OUR-ARDs score’ showed that the sum of scores ≥ 25 predicted mortality with a sensitivity-90%, specificity-64% and AUC of 0.85.

Conclusions

The ‘OUR ARDs’ risk score, derived from easily assessable factors predicting mortality, offered a tangible solution for prioritizing admission to COVID-19 tertiary care centre, that enhanced patient care but without unduly straining the health system.

Introduction

The Coronavirus Disease 2019 (COVID-19) pandemic with its protean manifestation remains an unpredictable debacle; the spectrum of presentation varying from asymptomatic infection to a fulminant systemic inflammatory syndrome unleashed by the cytokine storm. Mortality has been 1.33%, in India, with 3,39,71,607 confirmed cases and 4,50,782 deaths, while Tamil Nadu, a southern state recorded 26,78,265 cases and 35,783 deaths, cumulative as of 11th October 2021 [1].

Risk factors for mortality in COVID-19 reported in various studies included advanced age [213], male gender [6,7,9,11,13] and comorbidities [3,5,6,916] like diabetes mellitus, obesity, systemic hypertension, renal diseases, coronary artery disease [12] and malignancy. Apart from manifestations such as fever [8], cough [8], haemoptysis [2], dyspnoea [2,6,8], fatigue [8], loss of consciousness [2,6], laboratory parameters such as elevation of Neutrophil-to-Lymphocyte Ratio (NLR) [2,15], and increased levels of creatinine [3], lactate dehydrogenase (LDH) [24,10,12], direct bilirubin [2] and alanine aminotransferase [3], which provide early clues to the severity of disease, an increased plasma level of biomarkers like D-dimer [35,12], C-Reactive Protein (CRP) [6,7], serum ferritin [3], Interleukin-6 (IL-6) [3], and procalcitonin (PCT) fortifies these findings[35]. With the efficacy of specific antiviral and targeted immunomodulatory therapy still remaining elusive, prediction of mortality and risk stratification offers a rational approach for health resources allocation.

Hence, we retrospectively analysed data from a designated COVID-19 tertiary care centre in Tamil Nadu, South India, to unearth the risk factors predisposing to severe COVID-19 infection and mortality to formulate a risk score based on symptoms, comorbidities, vital signs and simple laboratory investigations, that will help to stratify patients into a high-risk group mandating preferential admission to tertiary care without undue delay, while safely retaining the low-risk group either at a primary care centre or at home, rendering the best care available, both clinically and psychosocially, ensuring the best survival advantage in both the groups, while simultaneously reducing the strain to the health care system.

Methods

Data collection

De-identified data from case records of COVID-19 confirmed patients (by real-time Reverse Transcription Polymerase Chain Reaction tests obtained from nasal or oropharyngeal swabs), hospitalized to a public sector COVID-19 tertiary care centre between 1st May 2020 and 30th November 2020 were analysed; the period roughly coinciding with the ascending slope of the epidemic curve in this region. Case records with insufficient data or without blood investigations were excluded. Institutional Ethics Committee approvals were duly obtained along with waiver from getting individual informed consents as the de-identified data were generated during the process of routine medical management, with the first recorded information on clinical status including comorbidities, peripheral oxygen saturation (SaO2), other vital parameters along with additional information on therapy and oxygen supplementation providing the raw data for analysis. Data were transcribed into an electronic format by two clinicians. An independent statistician checked the accuracy of data and analysed the data using IBM SPSS Statistics for Windows (IBM Corp. Released 2017, Version 25.0. Armonk, NY) and R software version 4.0.4 (R Core Team 2021).

Categorisation of outcomes

Patients were broadly classified into ‘survivors’ (discharged or transferred to step-down COVID-19 care centres) and ‘non-survivors’ (deaths during hospitalisation) to derive prognostic factors determining mortality or recovery.

Variables considered

The univariate analysis included demographic details describing age, gender, days to admission from the onset of symptoms and pre-existing comorbidities. Multiple comorbidities were present in some of the patients, which is not surprising by clinical parlance. Vital signs first recorded on admission that included SaO2, Pulse Rate, and Blood Pressure were copied along with basic laboratory reports consisting of plasma sugar (random), complete blood count, serum electrolytes, liver, and renal function tests. The relative difference in the symptomatology of COVID-19 between survivors and non-survivors was also explored to see if they played a role in envisaging future clinical courses.

Biomarkers like D-Dimer, Ferritin, PCT, and CRP levels were available only in a subset of the cohort. The receiver operating characteristic (ROC) curves of the above biomarkers were constructed and the area under the curves (AUC) was calculated and superimposed on the ROC curves of components used in the risk score that includes NLR, blood urea and SaO2 curves to provide and perceive their relative diagnostic accuracy with respect to these standard biomarkers. Since all patients were treated uniformly as per prevailing National guidelines [17], we probed only into the type of steroid used in treatment, and to determine if any specific steroid had an undue therapeutic advantage in combating COVID-19 triggered cytokine storm. Details of complications and treatment in the cohort along with the requirement of Oxygen supplementation are provided in the S1 Table.

Statistical analysis and risk score construction

Univariate analysis of baseline variables was performed using the Chi-square test for proportions and Mann Whitney-U test for continuous variables, compared between survivors and non-survivors, to unearth candidate predictors among symptoms, comorbidities, vital signs, and basic lab investigations elucidated at admission, with the parameter being available in at least 20% of patients. The predictive model was built, applying the stepwise multivariable logistic regression, selecting those variables that were found significant at the 10% level in the univariate analysis. The ‘goodness of fit’ of the model was assessed using the Hosmer-Lemeshow test. The predictive capability was ascertained using the area under the curve of receiver operator characteristics (AUROC) and validated using the bootstrap method [18], fitted into 1000 bootstrap datasets. The risk score for each important parameter was arrived at, by multiplying the regression coefficient of the final model with ten and rounding it off to the nearest next integer; with the parameters remaining significant in at least 50% of the samples chosen from bootstrap data [19]. The capability of the risk score to predict mortality was once again endorsed by the AUROC value.

Results

Of 746 case records of hospitalised patients, with near-complete details on clinical presentation including past history, comorbid conditions, symptoms at onset along with all baseline investigations available for prognosticating risk factors for mortality, there was a slight male preponderance in the occurrence of COVID-19 (62.5%, 466/746). There were 487 survivors and 259 non-survivors. The days to admission from the onset of symptoms were similar [median (IQR) 3(2–5) days] among survivors and non-survivors. The age group between 40 and 70 years was affected the most (59.1%, 441/746).

The distribution of symptoms among patients was as follows; fever—62.6% (455/727), breathlessness—49.9% (363/727), cough—46.4% (337/727), and gastrointestinal symptoms- 10.9% (79/727) in the decreasing order of frequency. Asymptomatic patients accounted for 6.7% (49/727).

There was at least one pre-existing comorbidity in 65.2% (475/728) of the cohort investigated with the major ones being diabetes mellitus (41.1%, 299/728), hypertension (30.1%, 219/728), cardiovascular disease (10.3%, 75/728), and kidney disease (9.9%, 72/728). Other comorbidities like malignancy, obesity and airway disease had a meagre representation in our cohort (<5%), precluding us from analysing their role further. Univariate analysis of baseline characteristics of the patients segregated into ‘survivors’ and ‘non-survivors’ is provided in Table 1.

thumbnail
Table 1. Baseline characteristics of hospitalised patients categorised into ‘survivors’ and ‘non-survivors’.

https://doi.org/10.1371/journal.pone.0263471.t001

Symptoms such as anosmia [8.1% (39/483) vs 1.2% (3/244)], myalgia [24.2% (117/483) vs 6.1% (15/244)], and sore throat [14.5% (70/483) vs 6.1% (15/244)] favoured survival (p<0.001) while breathlessness [37.3% (180/483) vs 75% (183/244), p<0.001] and altered sensorium [1.7% (8/483) vs 7.8% (19/244), p<0.001] signalled a perilous outcome. There was not a single case of rhinorrhoea documented in our cohort.

Among hospitalised patients, survival was 84.6% (214/253) when there was no comorbidity, which reduced to 54.9% (261/475) when at least one comorbidity existed [p<0.001], which was further exacerbated by advancing age>50, [46.6% (149/320)]. There was a disproportionately higher incidence of mortality among diabetic patients compared to non-diabetic patients [48.5%, (145/299) vs 25.2%, (108/429)], p<0.001], the risk doubling with elevated urea levels compared to those with preserved renal function within the diabetic population [65.4%, (68/104) vs 27.3%, (33/121); p<0.001].

SaO2<90% at admission proved detrimental with a striking disparity perceived in the median (IQR) of SaO2 among survivors 98% (94–99%) compared to non-survivors 88% (76.5–93.5%). Tachycardia (PR ≥100/min) posed an ominous sign.

Multivariable analysis and model creation

The backward stepwise elimination logistic predictive model was fitted with clinically relevant candidate variables significant (p<0.1) in the univariate analysis and that could be feasible at a primary healthcare unit or screening centre. This included age, gender, SaO2, blood urea, pulse rate, diabetes, hypertension, and NLR. Fig 1 provides the influence of individual parameters on survival depicted as Kaplan-Meier survival curves. Gender and hypertension did not retain their significance in the multivariate logistic regression after adjustment for age and elevated urea and were subsequently excluded from the final model. The Hosmer-Lemeshow test (χ2 = 6.27, p = 0.617) established the goodness of fit of the model for this data. The AUC or the concordance index of 0.85 (95% CI: 0.81–0.89) indicated an adequate discriminating power, excluding overfitting bias in the model, which was further validated in the 1000 bootstrap datasets. The regression coefficients of parameters retained in the model were found to be significant in at least 70% of the bootstrap samples. The sign of the regression coefficient was also found to be consistent across the bootstrap samples. The multivariable logistic regression showing the adjusted odds ratio of the various parameters, with their regression coefficients used in deriving the risk score is provided in Table 2.

thumbnail
Fig 1. Kaplan-Meier analysis of the parameters significantly attributing to risk of mortality in the univariate analysis.

Kaplan-Meier mean survival estimates for the time to mortality censored at 40 days. The numbers given above reflect the number of individuals who were alive at that particular time point in days, with or without the risk factors evaluated. Significance was computed using the log-rank test. Definition of Abbreviations: SaO2-Peripheral Oxygen Saturation, NLR-Neutrophil-to-Lymphocyte Ratio, BPM-Beats per minute.

https://doi.org/10.1371/journal.pone.0263471.g001

thumbnail
Table 2. Multivariable logistic regression showing the adjusted odds ratio of predictors influencing mortality, among COVID-19 infected patients admitted to hospital.

https://doi.org/10.1371/journal.pone.0263471.t002

Construction and calculation of risk score—‘OUR ARDs’ score

A regression coefficient-based scoring system, acronymed ‘OUR-ARDs’ for easy recollection, was derived from the individual parameters that emerged significant for influencing mortality in the multivariate regression; the risk scores being: SaO2 <95% - 11, Urea ≥50 mg/dl -15, NLR>3–11, Age ≥50 years—9, Heart Rate ≥100 BPM—7 and history of Diabetes Mellitus—6. We elucidated that a sum of 25, as cut off, served as the critical value for predicting mortality with a sensitivity of 90% and specificity of 64%.

Even though biomarkers were available only in a subset of patients, the AUROC values were robust that compelled us to provide an illustrated comparison of the parameters utilised for risk scores formulation in combination with the Biomarkers to understand their relative precision is provided in Fig 2.

thumbnail
Fig 2. Illustrated depiction of ROC curves of important parameters used in score formulation along with biomarkers.

The relative positions in the ROC curve also provide the precision of parameters, compared to biomarkers, in predicting mortality. Definition of Abbreviations: ROC Curve—Receiver Operating Characteristic Curve, CRP- C-Reactive Protein, SaO2—Peripheral Oxygen Saturation, NLR—Neutrophil-to-Lymphocyte Ratio.

https://doi.org/10.1371/journal.pone.0263471.g002

On the therapeutic aspect of steroid selection, mortality rates were similar between usage of methylprednisolone (42.4%, 120/283) and dexamethasone (44.6%, 41/92), [p = 0.717].

Discussion

The ‘OUR-ARDs’ score, an anamnesis of the underlying pathophysiology of COVID-19, showed that a score of ≥25 denoted a higher risk of mortality and needed priority admission to a tertiary care centre for escalation of management. Utilising clinically relevant determinants of mortality coupled with simple and basic laboratory investigations, gleaned from admission records of hospitalised COVID-19 patients, the OUR-ARDs score provides early clues and guidance in triaging of patients to the hospital, with the dual advantage of neither compromising on patient safety nor straining the health system. We found no difference between survivors and non-survivors with respect to time to admission from symptom onset, similar to other studies, suggesting a robust health system was in place [3,5].

While the SMART-COP and SCAP scores were operating at a higher or later level of medical care, oriented towards the transfer of in-patients in hospital wards to intensive care, our risk score was more focussed at the entry-level of health care, triaging patients into those who require a referral to a higher centre and those who can safely be managed at home or at the primary health care facility, which would prove useful when there is an exponential increase in admissions avoiding the congestion of the health system [20,21].

A comparative table of some of the risk scores aimed at triaging of RT-PCR positive COVID-19 patients, reported from various parts of the globe along with their translational application is briefly described for comparison and clinical utility as a ready reckoner in Table 3. While the risk scores such as the 4C mortality score [6] and the CALL score [10] were stratified with multiple levels, our risk score had a cut off level to decide on admission like the MuLBSTA score [22]. The clinical components of our score were similar to that of Galloway et al. [7] except for hypoalbuminemia as an additional factor in the latter.

thumbnail
Table 3. A concise comparison of risk scores used for triaging COVID-19 patients and their potential utilisation in clinical practice and public health.

https://doi.org/10.1371/journal.pone.0263471.t003

Exploring the individual components of the risk score, levels of SaO2 at admission critically determined survival [23] apart from elevated urea and NLR. Tachycardia, (heart rate >100/min) had been a well-known prognostic factor in pneumonia, with COVID-19 being no exception. However, the cut-offs for heart rate varied among individual studies [24].

An intact functioning kidney is often a prerequisite for survival in COVID-19. Elevated Urea above 40 mg/dl was a strong predictor of mortality, frequented more commonly in the background of chronic kidney disease proving detrimental to survival. Renal impairment has been associated with delayed clearance of cytokines, along with elevated levels of von Willebrand factor (vWF), adding to both exaggerated inflammatory response and pro-coagulable state. Pre-existing vasodynamic insults caused by underlying diabetes mellitus or hypertension further exacerbated the condition [25,26]. Mudatsir et al. proclaimed that elevated levels of urea levels confronted at the time of COVID-19 were due to the direct involvement of renal cells by the Coronavirus rather than due to pre-existing renal disease which puts any patient with an increasing urea level into the high-risk category irrespective of the background renal status [27].

An elevated NLR, as in other viral infections, carried a poor prognosis, serving as a surrogate indicator for survival and the need for ventilator support [15,28]. In COVID-19, the increase in neutrophil number, after infiltration, generated greater quantities of cytokines, leading to immune fury and multi-systemic inflammation, ending fatally. The study by Zhou C et al. confirmed that profound lymphopenia was the strongest influencer of mortality and rightfully retained a higher score for severity in our study equivalent to blood oxygen saturation. (Table 2). Even among healthy young individuals without any comorbidities, a higher NLR and lymphopenia <900 cells/mm3 predicted a poor prognosis [29]. The reasons postulated for this corresponding decrease in lymphocyte number include a reduction of circulating lymphocytes, apoptotic destruction, decreased production in the spleen and thymus due to direct involvement by the virus, and the negative feedback signals received from elevated plasma levels of IL-6. The meta-analysis by Huang W et al. showed a profound decline in the entire repertoire of the lymphocyte family in COVID-19 including CD4+, CD8+, B-cells and NK-cells [30].

Age ≥50 years was an independent factor contributing to severity. Survivors were two decades younger than Non-Survivors in our study and this feature has been commonly reported across all studies (Table 3). Various studies established a linear relationship between the severity of COVID-19 and the increasing number of comorbidities in the same individual, supplemented by advancing age [6,31]. Our study elucidated that the diabetic patients with concomitant elevated renal parameters were the most likely to succumb to COVID-19. Hypertension and diabetes exacerbate many viral infections, with higher mortalities reported in dengue, Yellow Nile Fever and Zika Virus infection [32]. Our data did show that both isolated systemic hypertension and cardiovascular disease to be of higher significance in the univariate analysis (Table 1). However, hypertension did not retain its significance in the multivariable analysis, after adjusting for age and elevated urea. This aspect needs careful interpretation in lieu of possible overlap and interdependence between hypertension and other factors like age, diabetes, and the tendency for renal compromise or insufficiency, possibly mitigating the attributable role of hypertension individually in causing COVID-19 related mortality in our study.

The mortality rate was less than one per cent in the absence of any comorbidity while patients with solid organ transplants formed the other end of the spectrum having the greatest causality [33]. Underlying malignancies, associated cardiovascular diseases [34], diabetes, COPD, and hypertension have been reported as important determinants of survival in various other studies [31]. Incidentally, some of these did not have sufficient representation in our study population to draw a firm conclusion.

Venturing into the symptomatology of COVID-19, breathlessness and altered sensorium at admission were logically linked to a poorer prognosis in line with the findings of Myrstad et al. [35], while anosmia, myalgia, and sore throat offered a better chance of survival. Certain symptoms like fever, cough, abdominal pain and vomiting were non-committal in predicting survival similar to the report by Dependra et al. [15]. Intriguingly, rhinorrhoea among COVID19 patients was a rarity [29]. This categorisation of COVID-19 symptomatology however could be ‘wave-specific’ and needs to be interpreted in the light of region-specific association with other objective findings such as SaO2, NLR, etc. for better clinical application. A meta-analysis from Indonesia found that patients with breathlessness, anorexia, and fatigue had a worse prognosis [27].

Fig 2 depicts the AUROC of the components of the OUR-ARDs risk score superimposed with biomarkers estimated in a subset of our population, to bring out the comparative precision of the components of the risk score with the former. Based on clinical intuition, we hypothesize that those reaching risk scores above 19 but below 25 may be re-evaluated periodically using the same risk score for maximising patient safety or additionally correlated with biomarker estimation to precisely estimate the risk (S2 Table) as biomarkers strengthen the decision.

Therapeutically, we found no differences between the efficacies of intravenous steroids namely dexamethasone and methylprednisolone in preventing mortality. However, other studies have shown the latter to be associated with a lesser requirement of ventilator support and a shorter stay in hospital [36,37].

The strength of this study lies in the utilisation of basic parameters, easily available and readily assessable even in a peripheral hospital without the requirement of complex medical gadgets/tests for predicting mortality and formulation of a risk score from the statistically significant factors. We had avoided using symptoms in the multivariate analysis due to its subjective nature, inherent heterogeneity, with those symptoms crucially determining mortality presenting late during the course of the disease along with contradictory features like apparent “silent hypoxia”, taking the physician by surprise. The study has important limitations. Selection bias and missing data were inevitable with the pandemic still on a rampage. The internal validation was done using the bootstrap method that had provided reliable estimates for the predictive model. With the analysis dictated by the availability of data, we could confidently allude those who were not included in the analysis (data unavailable in hospital records) predominantly belonged to asymptomatic, incidental, or milder cases of COVID-19, who did not require much medical attention as such. Data on biomarkers were limited as priority was given to the higher risk group at the time when patients were exponentially pouring in. But the AUC values of the biomarkers were robust that we could illustrate the accuracy of the components of risk scores concerning crucial biomarkers. A recent review by Rod JE et al. showed only CRP and D-dimer to be the consistently performing predictors across studies [38]. Considering these limitations and that the risk score had to be applied at the entry into the health care system, where a physician should take a crucial decision whether to refer to tertiary care or retain the patients at a primary health care system, where biomarkers may not be readily available, we had not included them in the current model suggested. However, we do admit that the biomarkers D-Dimer and CRP could greatly enhance the accuracy of the risk score. Even though other comorbidities could potentially have a bearing on COVID-19 survival, it was difficult to ascertain their influence on COVID-19 severity due to the smaller numbers available in the current cohort study. For instance, a study reported that asthmatics had a better outcome compared to COPD patients [4]. Our data does not truly represent case fatality or true mortality rates, as data was extracted from a portion of patients with near-complete clinical history and baseline investigations. With the available data, our risk score may not be directly translatable to post-COVID-19 extra-pulmonary complications that may also demand admission. Regional, strain and ethnic variations in COVID-19 presentation may have to be accounted for, while interpreting the results.

Conclusion

With the compelling global travel and COVID-19 sparing no country, predictors of COVID-19 fatality and stratification of risk using a simplified risk score could provide a tangible solution to rationalise resource allocation for best results and efficient patient recovery, importantly without undue strain on the health system.

Supporting information

S1 Table. Details of complications and treatment in the cohort along with supplementation of O2.

a-Heparin includes both Low Molecular Weight Heparin (LMWH) and unfractionated heparin (UFH).

https://doi.org/10.1371/journal.pone.0263471.s001

(DOCX)

S2 Table. Summated OUR-ARDs risk scores and their corresponding sensitivity and specificity.

The OUR-ARDs risk score calculated using the 6 parameters namely, Peripheral Oxygen Saturation, Urea, Neutrophil-to-Lymphocyte Ratio, Age, Heart Rate and Diabetes Mellitus are provided, along with the sensitivity and specificity of each of the risk scores.

https://doi.org/10.1371/journal.pone.0263471.s002

(DOCX)

S1 STROBE statement. The STROBE checklist of items for this manuscript.

https://doi.org/10.1371/journal.pone.0263471.s003

(PDF)

Acknowledgments

We are grateful to the Ministry of Health & Family Welfare, Government of India; the Indian Council of Medical Research and the Department of Medicine & Family Welfare, Government of Tamil Nadu; for permitting the study and voluntary offering the services of the facilities to be utilised for a successful publication. Our special thanks to Prof. (Dr.) Balram Bhargava MD DM, Director General-ICMR & Secretary-Department of Health Research, Government of India, and Dr. J. Radhakrishnan IAS, Principal Health Secretary to the Government of Tamil Nadu for their readiness to share the data and provide insights and mould the study for maximising benefits to the community; Dr. R. Narayana Babu MD, Director of Medical Education (Tamil Nadu) for his help; and Prof. (Dr.) Sridhar Rathinam MD, Superintendent-Government Hospital of Thoracic Medicine, Tambaram Sanatorium; for his guidance and mentorship throughout the study period.

We also thank the former Directors of ICMR-NIRT, Dr. P.R. Narayanan PhD, and Dr. Srikanth Prasad Tripathy MD, the Institutional Ethics Committees of NIRT and GCMCH, the senior scientists of Department of Health Research, ICMR-NIRT and Dr Sekar Lakshmanan PhD, Former Principal Technical Officer, ICMR-NIRT for their critical inputs and valuable suggestions that had helped to fine-tune the study and Mrs. Devaki, Senior Administrative Officer, ICMR-NIRT and staff of ICMR-NIRT’s Transport Department for their logistic support.

We express our gratitude to the relentless work of our collaborators from GCMCH, Prof. (Dr.) Hariharan MS, MCh (Medical Superintendent), Prof. (Dr.) Gowthaman MD (Professor of Medicine), the Resident Medical Officer, Professors, Assistant Professors and Associate Professors of various departments, Senior & Junior Medical Residents of various departments, Nursing Staff, and Paramedical Staff of the Government Chengalpattu Medical College & Hospital whose untiring work has helped to save hundreds of precious lives during the raging pandemic.

We are grateful to Dr. Amarnath, Mr. Sakthi Vadivel and support staff of the departments of Biochemistry, Community Medicine, Radio-Diagnosis, and Medical Records Department of GCMCH who had spent extra hours to make it successful and to Messrs. Rajasekaran Subramani M.Sc., and Gokulraj Chandramohan M.A., of ICMR-NIRT for their zealous involvement in data collection. Our gratitude extends to the Director of CLT&RI, Chengalpattu and Dr. Vinoth Kumar of CLT&RI for their excellent hospitality. We acknowledge all the healthcare workers involved in the diagnosis and treatment of COVID-19 patients in GCMCH, Chengalpattu, Tamil Nadu.

References

  1. 1. COVID-19 India Dashboard #IndiaFightsCorona COVID-19 [Internet]. Mygov.in. 2020 [cited 2021 Oct 11]. Available from: https://www.mygov.in/covid-19.
  2. 2. Liang W, Liang H, Ou L, Chen B, Chen A, Li C, et al. Development and Validation of a Clinical Risk Score to Predict the Occurrence of Critical Illness in Hospitalized Patients With COVID-19. JAMA Internal Medicine. 2020;180(8):1081–9. pmid:32396163
  3. 3. Zhou F, Yu T, Du R, Fan G, Liu Y, Liu Z, et al. Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study. Lancet (London, England). 2020;395(10229):1054–62. pmid:32171076
  4. 4. Li X, Xu S, Yu M, Wang K, Tao Y, Zhou Y, et al. Risk factors for severity and mortality in adult COVID-19 inpatients in Wuhan. The Journal of allergy and clinical immunology. 2020;146(1):110–8. pmid:32294485
  5. 5. Shang Y, Liu T, Wei Y, Li J, Shao L, Liu M, et al. Scoring systems for predicting mortality for severe patients with COVID-19. EClinicalMedicine. 2020;24:100426. pmid:32766541
  6. 6. Knight SR, Ho A, Pius R, Buchan I, Carson G, Drake TM, et al. Risk stratification of patients admitted to hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: development and validation of the 4C Mortality Score. BMJ (Clinical research ed). 2020;370:m3339. pmid:32907855
  7. 7. Galloway JB, Norton S, Barker RD, Brookes A, Carey I, Clarke BD, et al. A clinical risk score to identify patients with COVID-19 at high risk of critical care admission or death: An observational cohort study. The Journal of infection. 2020;81(2):282–8. pmid:32479771
  8. 8. Salunke AA, Pathak SK, Dhanwate A, Warikoo V, Nandy K, Mendhe H, et al. A proposed ABCD scoring system for patient’s self assessment and at emergency department with symptoms of COVID-19. Diabetes & metabolic syndrome. 2020;14(5):1495–501. pmid:32795741
  9. 9. Joshi SR. Indian COVID-19 Risk Score, Comorbidities and Mortality. The Journal of the Association of Physicians of India. 2020;68(5):11–2. pmid:32610858
  10. 10. Ji Dong, Zhang Dawei, Xu Jing, et al. Prediction for Progression Risk in Patients With COVID-19-19 Pneumonia: The CALL Score, Clinical Infectious Diseases, Volume 71, Issue 6, 15 September 2020, Pages 1393–1399, pmid:32271369
  11. 11. Steinberg E, Balakrishna A, Habboushe J, Shawl A, Lee J. Calculated decisions: COVID-19 calculators during extreme resource-limited situations. Emergency medicine practice. 2020;22(4 Suppl):Cd1–cd5. pmid:32259420
  12. 12. Chauhan NK, Shadrach BJ, Garg MK, Bhatia P, Bhardwaj P, Gupta MK, et al. Predictors of Clinical Outcomes in Adult COVID-19 Patients Admitted to a Tertiary Care Hospital in India: an analytical cross-sectional study. Acta Biomed. 2021;92(3):e2021024. pmid:34212921
  13. 13. Ahmad S, Kumar P, Shekhar S, Saha R, Ranjan A, Pandey S. Epidemiological, Clinical, and Laboratory Predictors of In-Hospital Mortality Among COVID-19 Patients Admitted in a Tertiary COVID Dedicated Hospital, Northern India: A Retrospective Observational Study. J Prim Care Community Health. 2021;12:21501327211041486. pmid:34427136
  14. 14. Nakakubo S, Suzuki M, Kamada K, Yamashita Y, Nakamura J, Horii H, et al. Proposal of COVID-19 Clinical Risk Score for the management of suspected COVID-19 cases: a case control study. BMC infectious diseases. 2020;20(1):858. pmid:33208116
  15. 15. Rai D, Ranjan A, H A, Pandey S. Clinical and Laboratory Predictors of Mortality in COVID-19 Infection: A Retrospective Observational Study in a Tertiary Care Hospital of Eastern India. Cureus. 2021;13(9):e17660. pmid:34646702
  16. 16. Varghese GM, John R, Manesh A, Karthik R, Abraham OC. Clinical management of COVID-19. The Indian journal of medical research. 2020;151(5):401–10. pmid:32611911
  17. 17. Clinical Management Protocol for COVID-19. Gov.in. [cited 2021 Oct 11]. Available from: https://COVID-19.india.gov.in/document/clinical-management-protocol-for-COVID-19-3/.
  18. 18. Austin PC, Tu JV. Bootstrap Methods for Developing Predictive Models. The American Statistician. 2004;58(2):131–7.
  19. 19. Mehta HB, Mehta V, Girman CJ, Adhikari D, Johnson ML. Regression coefficient-based scoring system should be used to assign weights to the risk index. Journal of clinical epidemiology. 2016;79:22–8. pmid:27181564
  20. 20. Charles PG, Wolfe R, Whitby M, Fine MJ, Fuller AJ, Stirling R, et al. SMART-COP: a tool for predicting the need for intensive respiratory or vasopressor support in community-acquired pneumonia. Clin Infect Dis. 2008;47(3):375–84. pmid:18558884
  21. 21. Espana PP, Capelastegui A, Gorordo I, Esteban C, Oribe M, Ortega M, et al. Development and validation of a clinical prediction rule for severe community-acquired pneumonia. Am J Respir Crit Care Med. 2006;174(11):1249–56. pmid:16973986
  22. 22. Preetam M, Anurag A. MuLBSTA score in COVID-19 pneumonia and prediction of 14 day mortality risk: A study in an Indian cohort. J Family Med Prim Care 2021; 10:223–7. pmid:34017730
  23. 23. Li J, He X, Yuan Y, Zhang W, Li X, Zhang Y, et al. Meta-analysis investigating the relationship between clinical features, outcomes, and severity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pneumonia. American journal of infection control. 2021;49(1):82–9. pmid:32540370
  24. 24. Ahn JH, Choi EY. Expanded A-DROP Score: A New Scoring System for the Prediction of Mortality in Hospitalized Patients with Community-acquired Pneumonia. Scientific reports. 2018;8(1):14588. pmid:30275523
  25. 25. Bhargava A, Fukushima EA, Levine M, Zhao W, Tanveer F, Szpunar SM, et al. Predictors for Severe COVID-19 Infection. Clinical infectious diseases: an official publication of the Infectious Diseases Society of America. 2020;71(8):1962–8.
  26. 26. Khan IH, Savarimuthu S, Leung MST, Harky A. The need to manage the risk of thromboembolism in COVID-19 patients. Journal of vascular surgery. 2020;72(3):799–804. pmid:32417304
  27. 27. Mudatsir M, Fajar JK, Wulandari L, Soegiarto G, Ilmawan M, Purnamasari Y, et al. Predictors of COVID-19 severity: a systematic review and meta-analysis. F1000Research. 2020;9:1107. pmid:33163160
  28. 28. Garcia-Gordillo JA, Camiro-Zúñiga A, Aguilar-Soto M, Cuenca D, Cadena-Fernández A, Khouri LS, et al. COVID-IRS: A novel predictive score for risk of invasive mechanical ventilation in patients with COVID-19. PloS one. 2021;16(4):e0248357. pmid:33819261
  29. 29. Zhou C, Huang Z, Tan W, Li X, Yin W, Xiao Y, et al. Predictive factors of severe coronavirus disease 2019 in previously healthy young adults: a single-center, retrospective study. Respiratory research. 2020;21(1):157. pmid:32571410
  30. 30. Huang W, Berube J, McNamara M, Saksena S, Hartman M, Arshad T, et al. Lymphocyte Subset Counts in COVID-19 Patients: A Meta-Analysis. Cytometry Part A: the journal of the International Society for Analytical Cytology. 2020;97(8):772–6. pmid:32542842
  31. 31. Guan WJ, Liang WH, Zhao Y, Liang HR, Chen ZS, Li YM, et al. Comorbidity and its impact on 1590 patients with COVID-19 in China: a nationwide analysis. The European respiratory journal. 2020;55(5). pmid:32217650
  32. 32. Badawi A, Velummailum R, Ryoo SG, Senthinathan A, Yaghoubi S, Vasileva D, et al. Prevalence of chronic comorbidities in dengue fever and West Nile virus: A systematic review and meta-analysis. PloS one. 2018;13(7):e0200200. pmid:29990356
  33. 33. Pereira MR, Mohan S, Cohen DJ, Husain SA, Dube GK, Ratner LE, et al. COVID-19 in solid organ transplant recipients: Initial report from the US epicenter. American journal of transplantation: official journal of the American Society of Transplantation and the American Society of Transplant Surgeons. 2020;20(7):1800–8. pmid:32330343
  34. 34. Mullen B. COVID-19 Clinical Guidance for the Cardiovascular Care Team. American College of Cardiology. Acc.org. [cited 2021 Oct 11]. Available from: https://www.acc.org/~/media/665AFA1E710B4B3293138D14BE8D1213.pdf.
  35. 35. Myrstad M, Ihle-Hansen H, Tveita AA, Andersen EL, Nygård S, Tveit A, et al. National Early Warning Score 2 (NEWS2) on admission predicts severe disease and in-hospital mortality from Covid-19—a prospective cohort study. Scandinavian journal of trauma, resuscitation and emergency medicine. 2020;28(1):66. pmid:32660623
  36. 36. Ranjbar K, Moghadami M, Mirahmadizadeh A, Fallahi MJ, Khaloo V, Shahriarirad R, et al. Methylprednisolone or dexamethasone, which one is superior corticosteroid in the treatment of hospitalized COVID-19 patients: a triple-blinded randomized controlled trial. BMC infectious diseases. 2021;21(1):337. pmid:33838657
  37. 37. Ko JJ, Wu C, Mehta N, Wald-Dickler N, Yang W, Qiao R. A Comparison of Methylprednisolone and Dexamethasone in Intensive Care Patients With COVID-19. Journal of intensive care medicine. 2021;36(6):673–80. pmid:33632000
  38. 38. Rod JE, Oviedo-Trespalacios O, Cortes-Ramirez J. A brief-review of the risk factors for covid-19 severity. Rev Saude Publica. 2020;54:60. pmid:32491116