key: cord-0793213-b4tyed8a
authors: Li, Juan; Wang, Lili; Liu, Chun; Wang, Zhengquan; Lin, Yi; Dong, Xiaoqi; Fan, Rui
title: Exploration of prognostic factors for critical COVID-19 patients using a nomogram model
date: 2021-04-14
journal: Sci Rep
DOI: 10.1038/s41598-021-87373-x
sha: 48f52868bbb1e0027517aab7e6a9afafaf721c9a
doc_id: 793213
cord_uid: b4tyed8a

The study aimed to explore the influencing factors on critical coronavirus disease 2019 (COVID-19) patients’ prognosis and to construct a nomogram model to predict the mortality risk. We retrospectively analyzed the demographic data and corresponding laboratory biomarkers of 102 critical COVID-19 patients with a residence time ≥ 24 h and divided patients into survival and death groups according to their prognosis. Multiple logistic regression analysis was performed to assess risk factors for critical COVID-19 patients and a nomogram was constructed based on the screened risk factors. Logistic regression analysis showed that advanced age, high peripheral white blood cell count (WBC), low lymphocyte count (L), low platelet count (PLT), and high-sensitivity C-reactive protein (hs-CRP) were associated with critical COVID-19 patients mortality risk (p < 0.05) and these were integrated into the nomogram model. Nomogram analysis showed that the total factor score ranged from 179 to 270 while the corresponding mortality risk ranged from 0.05 to 0.95. Findings from this study suggest advanced age, high WBC, high hs-CRP, low L, and low PLT are risk factors for death in critical COVID-19 patients. The Nomogram model is helpful for timely intervention to reduce mortality in critical COVID-19 patients.

Nomogram construction and validation. Critical COVID-19 patient mortality risk factors concluded from Model 4 were input into the software R. An estimation nomogram of critical COVID-19 patient mortality risk was then built. Each item in the model corresponded to a score, the sum of each item score became a total score. The point on the critical COVID-19 patient mortality risk axis corresponds to patient mortality risk. In the nomogram, the total scores ranged from 179 to 270, with corresponding mortality risks of 0.05 to 0.95. The higher the score, the higher the mortality risk (Fig. 1 ). In addition, the area under the receiver operating characteristic (ROC) curve in Table 3 was utilized to calculate the accuracy for different test indexes on the prediction of survival. The results showed that the area under the curve of the prediction model was the largest. Across all test indexes, hs-CRP had the largest area under the curve, which was the closest to the prediction model.

The C-index for this nomogram was 0.958 (95% CI 0.923-0.993), indicating a good discriminative ability. The calibration curve (Fig. 2) showed consistency between prediction probability and actual probability, which indicates the accuracy of the nomogram in estimating critical COVID-19 mortality risk.

The new coronavirus disease, COVID-19, is an acute respiratory infectious disease caused by the SARS-CoV-2 virus. As of January 22nd, 2021, there have been more than 9.8 million confirmed cases of COVID-19 globally, and the number of deaths had reached more than 2 million 8 . Investigations have found that only 80% of confirmed cases showed mild illness, but an admission survey of 166 British hospitals showed that around 33% of diagnosed patients required hospitalization and 45% of patients in intensive care died from the disease 9 . Table 1 . Comparability of clinical data and laboratory indicators between the death and survival group. Continuous variables with normal distribution were expressed as the mean ± standard deviation (SD), non-normal variables were expressed as the median (interquartile range (IQR)), and categorical data were expressed as number and percentage. The independent sample Student's t-test was used to compare the means of two continuous normally distributed variables. The means of two non-normally distributed variables were compared with the Mann-Whitney U test. The frequencies of categorical variables were compared by the χ 2 test.

Survival group (n = 50) Death group (n = 52) t/Z/χ 2 p www.nature.com/scientificreports/ Based on the comparison of various indicators of COVID-19 in patients in the survival group, our findings indicated that patients with advanced age, high WBC, low L, low PLT, and high hs-CRP were at high risk of mortality. A nomogram model including age, WBC, L, PLT, and CRP was constructed, which may be used for early and intuitive identification and evaluation of the prognosis of critical COVID-19 patients and provide a basis for personalized prevention and treatment.

COVID-19 infections appear in clusters, and are more likely to affect elderly males with comorbidities 10-12 . We analyzed patient's comorbidity and Model 2 showed that coronary heart disease was a risk factor for death in patients with severe COVID-19 (OR = 5.222). Zheng et al. 13 found that patients with cardiovascular disease have increased angiotensin converting enzyme 2 (ACE-2), compared to healthy individuals, while Goldstein et al. 14 have discussed that angiotensin receptor blocker's (ARB's) facilitate viral entry by increasing the expression Table 2 . Multivariate logistic regression analysis for risk of death in critical COVID-19 patients. Model 1: logistic regression analysis of age, sex, smoking habits, body temperature; Model 2: logistic regression analysis of the presence of hypertension, diabetes, and coronary heart disease; Model 3: logistic regression analysis of peripheral white blood cell count (WBC), lymphocyte count (L), platelet count (PLT), high-sensitivity C-reactive protein (hs-CRP), estimated glomerular filtration rate (eGFR), d-dimer (D-D), and troponin I (TnI); Model 4: logistic regression analysis of all items above. 15 . Changes in ACE-2 levels in the body can affect the pathogenicity of SARS-CoV-2. The high expression of ACE-2 enhances the ability of SARS-CoV-2 to infect the body, and the level of ACE-2 in patients with coronary heart disease is significantly higher than patients without. Similarly, the large percentage of patients with hypertension undertaking angiotensinconverting enzyme inhibitors (ACEIs) or ARBs is a remarkable characteristic among those presenting severe COVID-19 manifestations. However, Cao et al. have shown that the East Asian population has a higher expression than other countries 16, 17 . Only 16 of the 102 patients in this study had coronary heart disease, therefore the results of the survey were relatively limited. Whether coronary heart disease is associated with increased mortality in COVID-19 requires further study and there are also contradictory scenarios suggesting that the use of ACEIs/ARBs may exacerbate the deleterious conditions of the infection, therefore, this remains an area requiring further investigation. Many studies have considered high temperature to be related to the progression of COVID-19 13, 14, 18 . In this study, we found that body temperature in the death group was slightly lower than in the survival group. During the outbreak of COVID-19 in Wuhan, many patients were not able to be admitted into hospital. Many underwent cooling and other treatments at home at that time, so the body temperature measured after patients were admitted in this study may have been affected. As with the study by Chen et al., patients were also from Wuhan city, meaning there is a small number of fever patients included 19 . Model 1 showed that body temperature was different between the two groups in this study (OR = 0.513), but Model 4 did not find a difference in body temperature after adjusting for more factors, suggesting that the association between body temperature and prognosis may not be clear and requires further study.

Nomograms can show the probability of clinical events in statistical prediction models as simple scores through graphics, which can be used to predict prognosis by estimating clinical events and integrating significant prognostic factors in numerous diseases [20] [21] [22] . In the present study, we constructed a nomogram including age, WBC, L, PLT, and hs-CRP, which were screened through multivariate logistic regression analysis (Model 4). Individuals with higher total score from the nomogram had a greater risk of an unexpected prognosis. Through discriminatory capacity and calibration ability analysis, the results showed that this model was stable and had a high ability to predict prognosis of critical COVID-19 patients (C-index: 0.958, 95% CI 0.923-0.993). The calibration curve (Fig. 2) showed that the predicted probability and the actual probability were consistent with Compared to other studies which also constructed a nomogram for assessing survival in COVID-19 patients 19, 23, 24 , the indicators used in the current study are more basic and easier to obtain. However, there remains limitations. First, this was a retrospective study with a relatively limited number of patients, the causal relationship between the influencing factors and critical COVID-19 patients' prognosis was not certain. Second, although the confounding factors were adjusted by logistic regression, and the nomogram model was demonstrated to exhibit appropriate discriminatory capacity and calibration ability, the lack of symptomatic data and specific medications may affect the conclusions to some extent. Third, the participants were recruited from a single hospital, we therefore conducted an internal validation but lack externally validation.

Sample collection. This was a retrospective study with the samples being obtained from Tongji Hospital, Guanggu District, Wuhan City, Hubei, China from February 2020 to March 2020. A total of 102 critical COVID-19 patients were included and divided into the survival group and death group according to patient prognosis. The diagnosis and clinical classification of COVID-19 were based on the guidelines of the treatment of new coronavirus pneumonia (version 7) published by the National Health Commission of China 25 . The diagnosis was established on the basis of (1) epidemiological history, (2) fever and/or other respiratory symptoms, (3) presence or absence of imaging findings of novel coronavirus pneumonia, and (4) real time fluorescence RT-PCR for SARS-CoV-2 nucleic acid which yielded positive results. As for clinical classification, patients were deemed critical if they met at least one of the following criteria: (1) shortness of breath, with respiratory rate ≥ 30 breaths/min; (2) oxygen saturation (at resting state) ≤ 93%; (3) PaO 2 /FiO 2 ≤ 300 mmHg; (4) radiographical imaging showing lesion progression more than 50% within 24-48 h; or (5) respiratory failure, shock, or other organ failures. In addition, all critical patients were hospitalized for ≥ 24 h. This study was approved by the Institutional Ethics Committee of Union Hospital, Tongji Medical College, Huazhong University of Science and Technology (TJ-IRB20201014); the requirement for informed consent was exempted by the Ethics Committee. All methods were carried out in accordance with relevant guidelines and regulations.

Data collection and measurement. Demographic data including age, sex, smoking habits, body temperature, and comorbidity data including hypertension, diabetes, and coronary heart disease were collected. Fasting blood samples were extracted by venipuncture to collect clinical indicators in a core laboratory with a standard protocol. An ADVIA2400L automated biochemistry analyzer (Siemens AG, Munich, Germany) was used to measure peripheral WBC, L, PLT, hs-CRP, and D-D. A VITROS5600 (Johnson & Johnson, American) biochemistry analyzer was used to measure TnI. Clinical indicators were measured in a core laboratory with a standard protocol.

Statistical analysis. All statistical analyses were carried out using PASW Statistics 19.0 (SPSS, Inc., Somers, NY, USA) and R software version 3.6.1 (http:// www.R-proje ct. org). The Kolmogorov-Smirnov test was used to measure data normality. If the data satisfied normal distribution, mean ± standard deviation (SD) was shown and a T-test was used to find differences between two groups. Abnormally distributed data was shown as median (IQR), while a risk-sum test being used for comparison (Mann-Whitney U test). Categorical data were expressed as number and percentage and compared with the Chi-square test. Furthermore, logistic regression analysis by the stepwise method (Forward wald) was used. A nomogram was constructed based on the results of the logistic regression to predict the mortality risk using the nomogram function in the R software. Harrell's C-statistic was used to calculate C-index and to verify nomogram discrimination. Calibration curves were used to verify nomogram conformity. C-index and calibration curve were obtained by Bootstraps (1,000 multiple sampling). A two-sided p < 0.05 was considered statistically significant.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

The reproductive number of COVID-19 is higher compared to SARS coronavirus

Clinical features of patients infected with 2019 novel coronavirus in Wuhan

The novel coronavirus originating in Wuhan, China: Challenges for Global Health Governance

National Health Commission of the People's Republic of China

Clinical laboratory monitoring and analysis of COVID-19

New biomarkers in community-acquired pneumonia: Another step in improving outcome prediction

Analysis of clinical characteristics of 49 patients with coronavirus disease 2019 in Jiangxi

World Health Organization. Coronavirus disease (COVID-19) pandemic

Features of 20 133 UK patients in hospital with covid-19 using the ISARIC WHO Clinical Characterisation Protocol: Prospective observational cohort study

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: A descriptive study

Clinical characteristics of coronavirus disease 2019 in China

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

COVID-19 and the cardiovascular system

Are certain drugs associated with enhanced mortality in COVID-19?

Clinical and epidemiological features of 36 children with coronavirus disease 2019 (COVID-19) in Zhejiang, China: An observational cohort study

Comparative genetic analysis of the novel coronavirus (2019-nCoV/SARS-CoV-2) receptor ACE2 in different populations

Novel coronavirus (2019-nCoV): update on 3rd coronavirus outbreak of 21st century

Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation

Development and validation of a nomogram using on admission routine laboratory parameters to predict in-hospital survival of patients with COVID-19

A clinical nomogram construction method using genetic algorithm and naive Bayesian technique

How to build and interpret a nomogram for cancer prognosis

A nomogram to preoperatively predict 1-year disease-specific survival in resected pancreatic cancer following neoadjuvantchemoradiation therapy

Predictors for severe COVID-19 infection

Development and validation of a nomogram for assessing survival in patients with COVID-19 pneumonia

National Health Commission of China

This study was supported by the General Research Project of Zhejiang Provincial Department of Education (Y202043652).

L.L.W. and C.L. designed the research. R.F. and Y.L. reviewed the literatures. L.L.W. and Z.Q.W. collected the data. R.F. and Y.L. analyzed the data. L.L.W. and J.L. wrote the paper. J.L. and X.Q.D. modified the manuscript. All authors critically reviewed the manuscript for important intellectual content and approved the final version.

The authors declare no competing interests.

Correspondence and requests for materials should be addressed to R.F.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.