key: cord-0944550-tiqkunae authors: Zhao, Sheng-Chao; Yu, Xian-Qiang; Lai, Xue-Feng; Duan, Rui; Guo, De-Liang; Zhu, Qian title: Dose-response relationship between risk factors and incidence of COVID-19 in 325 hospitalized patients: A multicenter retrospective cohort study date: 2022-04-06 journal: World J Clin Cases DOI: 10.12998/wjcc.v10.i10.3047 sha: 00fe95c1170dcf29b3b425d54261af2eb85a21fa doc_id: 944550 cord_uid: tiqkunae BACKGROUND: The epidemiological and clinical characteristics of coronavirus disease 2019 (COVID-19) patients have been widely reported, but the assessment of dose-response relationships and risk factors for mortality and severe cases and clinical outcomes remain unclear. AIM: To determine the dose-response relationship between risk factors and incidence of COVID-19. METHODS: In this retrospective, multicenter cohort study, we included patients with confirmed COVID-19 infection who had been discharged or had died by February 6, 2020. We used multivariable logistic regression and Cox proportional hazard models to determine the dose-response relationship between risk factors and incidence of COVID-19. RESULTS: It clarified that increasing risk of in-hospital death were associated with older age (HR: 1.04, 95%CI: 1.01-1.09), higher lactate dehydrogenase [HR: 1.04, 95% confidence interval (CI): 1.01-1.10], C-reactive protein (HR: 1.10, 95%CI: 1.01-1.23), and procalcitonin (natural log-transformed HR: 1.88, 95%CI: 1.22-2.88), and D-dimer greater than 1 μg/mL at admission (natural log transformed HR: 1.63, 95%CI: 1.03-2.58) by multivariable regression. D-dimer and procalcitonin were logarithmically correlated with COVID-19 mortality risk, while there was a linear dose-response correlation between age, lactate dehydrogenase, D-dimer and procalcitonin, independent of established risk factors. CONCLUSION: Higher lactate dehydrogenase, D-dimer, and procalcitonin levels were independently associated with a dose-response increased risk of COVID-19 mortality. In December 2019, many cases of unknown viral pneumonia were reported [1] . A novel coronavirus, capable of infecting humans, was detected in January 2020 [2, 3] and the disease caused was termed coronavirus disease 2019 (COVID-19) by WHO [2] . As of March 19, 2020 > 200000 Laboratory-confirmed cases had been documented globally[1, [4] [5] [6] [7] . With the increasing awareness of COVID-19 pneumonia, a variety of diagnostic protocols and guidelines have evolved to guide clinical practice[1, 5, [8] [9] [10] . Many patients in some case series which had been published, were hospitalized at the time of reporting. Studies in patients who were not discharged may have misclassified outcomes due to patients developing severe disease or dying during subsequent hospitalization. Consequently, It might be inaccurate and unreliable to estimate of risk factors for severe illness and death in these early case series. Furthermore, although several large studies have reported risk factors for mortality and severe disease in COVID-19 patients, studies that systematically explored the potential associations were limited. Thus, the association of risk factors with COVID-19 outcomes remained unknown. Therefore, we detail all laboratory-confirmed COVID-19 patients admitted to two designated hospitals as of February 2020, along with clear clinical outcomes (death or discharge). The purpose of this study was to investigate risk factors for death in hospital and to clarify hospitalization characteristics of COVID-19 patients. This retrospective cohort study included two cohorts of adult inpatients (≥ 20 years old) from two designated hospitals. All patients who were diagnosed with COVID-19 according to the WHO interim guideline were screened [2] , and those who died or were discharged by February 6, 2020 were included in the study. We extracted demographic, clinical, laboratory, treatment, and outcome data from the hospital electronic medical records using a standardized data collection form modified from the version of the WHO/International Severe Acute Respiratory and Emerging Infection Consortium. The primary endpoint was in-hospital death occurring beyond 24 h but within 28 d and composite severe cases referred to the admission to intensive care unit, intubation, or death during hospitalization. Patients, who had normalized temperature for over 3 days, relief of clinical symptoms, substantial improvement in the imaging of both lungs and throat-swab samples negative twice for at least 24h apart, were allowed to be discharged. All data were collected by two physicians, double-checked independently, and verified by a third researcher. The computed tomography (CT) demonstrations were described according to the internationally standard nomenclature defined by the Fleischner Society and peerreviewed literature on viral pneumonia, using the terms including ground glass opacity (GGO), crazypaving pattern, and consolidation [11, 12] . A semi-quantitative scoring system was used to quantitatively estimate the pulmonary involvement of all these abnormalities on the basis of the area involved [13] . The total CT score was the sum of the individual lobar scores and ranged from 0 (no involvement) to 25 (maximum involvement). The distribution of lung abnormalities was recorded as predominantly subpleural (involving mainly the peripheral one-third of the lungs), random (without predilection for subpleural or central regions), or diffuse (continuous involvement without respect to lung segments) [14] . This study was approved by the Ethics of Committees of Zhongnan Hospital of Wuhan University, and in accordance with the Helsinki Declaration. Written informed consent was obtained from all patients before examination. The anonymous data was collected and analyzed to optimize clinical decision and treatment. Continuous and categorical variables were presented as median and n (%), respectively. We used the Mann-Whitney U test, χ² test, or Fisher's exact test to compare differences between survivors and nonsurvivors where appropriate. To explore the risk factors associated with in-hospital death, multivariable logistic regression models and the Cox proportional hazards model was used to determine the independent factors, which were based on the variables selected by a univariate analysis. To generate the Receiver operating characteristic (ROC) curves, patients were classified as survivor or non-survivors and CT total score of different stages excluded patients who were lost to follow-up. We compared patients' characteristics between the two hospitals and used a generalized linear model to adjust for possible differences in patients' characteristics and treatment between the two study centers. Statistical tests and P values were two-sided. Differences were considered significant with a value of P < 0.05. All statistical analyses were carried out using the SAS software (version 9.4), unless otherwise indicated. We assessed potential dose-response associations of incident COVID-19 mortality and severe cases risk by restricted cubic splines logistic and Cox regression using 3 knots at 25 th , 50 th , and 75 th percentiles of the corresponding risk factors with the median value of the above risk factors as the reference group [15] . The basic characteristics of the patients are shown in Table 1 . Medical workers accounted for 5.8% (19/325), and those who had a history of contact with wildlife accounted for 1.2%. The median incubation period was 6 d (interquartile range, 2-15 d). The median age was 45 years (interquartile range, 34-61 years). Female accounted for 57.8%. 77.5% of patients had fever on admission and 85.5% had fever during hospitalization. The second most common symptoms were cough (63.7%) and fatigue (48.0%), but nausea or vomiting (7.7%) and difficulty breathing (4.6%) were uncommon. In the total population, 21.2 % have at least one co-existing condition (e.g., hypertension and diabetes). At admission, the severity of COVID-19 was classified as not severe 265 cases and severe 60 cases. Patients with severe disease had a median age of 16 years older than those without severe disease, and any comorbidities were more common (66.7% vs26.4%), but exposure histories were similar. On admission, lymphocytopenia, thrombocytopenia and leukopenia were present in 61.8%, 19.4% and 28.6%, respectively. Most patients (69.5%) had increased C-reactive protein (CRP) levels. Laboratory abnormalities, including lymphocytopenia and leukopenia, were more pronounced in critically ill patients than in non-critically ill patients (Table 2) . 1 Quantitative data were presented as mean ± SD (minimum-maximum), while the counting data were presented as count (percentage of the total). IQR: Inter quartile range. All patients underwent computed tomography scans at the time of admission, and 97.8% revealed abnormal results. The most common patterns on chest CT were GGO (61.1%) and bilateral patchy shadowing (84.7%). No CT abnormality was found in seven of 308 (2.2%) patients who survived and in none of 17 patients who died. GGO, crazy-paving pattern and consolidation were the most frequent CT findings in mild COVID-19 pneumonia (Supplementary Figure 1) . Most patients (279/325), the total CT score increasedabout10 d after the onset of symptoms, and then gradually decreased (Table 3 , Supplementary Figure 2 ). There were statistically significant differences between the bilateral lower lobe CT scores at stage 1 and the corresponding upper/middle lobe CT scores (left lower lobe vs left upper lobe: 1 ± 1 vs 0 ± 1, P < 0.001, P = 0.004) ( Table 3 ). According to the degree of lung involvement and the quartile of patients 0-26 days after onset, there were six stages starting from the onset of symptoms ( Table 4 , Supplementary Figure 3 ). Overall, subpleural lesions were more common than changes in central lung disease. Bilateral lung involvement occurred in most patients during the course of the disease (Supplementary Figure 4) . ROC curve analysis showed that the area under the curve (AUC) of stage 5 disease was higher than either of stage, and the combined AUC for stages 2 and 5 was highest among all stages ( Figure 1 ). After univariate analysis, patients with diabetes or hypertension had a higher chance of death in hospital (Tables 1 and 2) . Age, sex, leukocytosis, and elevated glucose level, lactate dehydrogenase (LDH), high-sensitivity C-reactive protein (CRP), D-dimer (DD), total cholesterol, triglyceride, creatinine, and procalcitonin (PCT) were associated with death or severe illness. Older age [hazard ratio (HR):1.04, 95% confidence interval (CI): 1.01-1.09], higher LDH (HR: 1.04, 95% 95%CI:1.01-1.10), higher CRP (HR:1.10, 95%CI: 1.01-1.23), and elevated PCT (log n transformed HR: 1.88, 95%CI: 1.22-2.88), and DD > 1 μg/mL at admission (log n transformed HR: 1.63, 95%CI: 1.03-2.58) were associated with increasing odds of in-hospital death (Table 5) . Furthermore, DD and PCT were loglinearly correlated with COVID-19 mortality risk, while there were linear dose-response correlations between age, LDH, DD and PCT. In particular, It was evident that the dose-response association of LDH and PCT occurred in severe patients (all P for overall association < 0.05). The dose-response relationship between LDH and PCT was more obvious in severe patients in the meantime (all P for interaction < 0.05) (Figure 2 , Tables 5 and 6) . ROC curve analysis indicated that the combined AUC for age, sex, high-sensitivity CRP, DD, LDH and PCT (0.947) was higher than that of any one of these variables alone (Figure 1 ). These results show that combination of age, sex, high-sensitivity CRP, DD, LDH and PCT was more precise in predicting clinical outcome than single factors alone. Consistent with most studies[1,2,16], we found that the clinical features of COVID-19 were similar to those of SARS. Fever, cough, gastrointestinal symptoms were rare[17]. Lymphocytopenia was common, a finding that was consistent with two recent reports[1, 16] . We found that the fatality rate (5.2%) was lower than recently reported [1, 16] . This may be due to differences in sample size and case inclusion criteria. Our findings were higher than the national official statistics, which showed a mortality rate of 3.9% among 81003 cases of COVID-19 as of March 13, 2020. In this study, patients underwent multiple lung CT scans (≥ 3 times), providing reliable dynamic radiographic pattern data. During the first 2 wk, the number and severity of abnormal lesions on chest CT increased. Subsequently, there was a short plateau phase and a gradual decrease in abnormalities. There were six stages of lung involvement in patients who have recovered from COVID-19, which could be more accurately evaluate the time course of lung changes, compared with the previous 4 stages [18] . Combined profiling of stages 2 and 5provides a more precise clinical outcome prediction than conventional stages 1-4 classification [18] , suggesting a novel valuable prognostic indicator for COVID-19 patients after antiviral therapy. Our retrospective cohort study demonstrated several risk factors for death in patients who were hospitalized with COVID-19. Particularly, older age, LDH > 285 U/L, creatinine > 111 ng/mL, PCT > 0.05 ng/mL, and DD > 1 μg/mL on admission were associated with higher odds of in-hospital death. Previously, older age, DD > 1 μg/mL and sequential organ failure assessment (SOFA) score (including creatinine level) have been reported as important independent predictors of mortality in COVID-19 [4] , which is in accordance with our current study. The most plausible explanation included an agedependent defect in T-cell and B-cell function and excess type 2 cytokines, which predispose to ischemia Table 3 The computed tomography score of the pulmonary involvement in four stages Total CT score of the pulmonary involvement 2 ± 4 (0-18) 5 ± 5 (0-22) 7 ± 7 (0-22) 7 ± 7 (0-25) 5 ± 7 (0-24) 4 ± 6 (0-25) < 0.0001 1 Number of involved lobes 22 ± 2 (0-5) 3 ± 2 (1-5) 4 ± 2 (1-5) 3 ± 2 (1-5) 3 ± 2 (1-5) 4 ± 2 (1-5) < 0.0001 1 CT score in each lobe < 0.0001 1 Left upper lobe 0 ± 1 (0-3) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-4) 1 ± 1 (0-5) Left lower lobe 1 ± 1 (0-5) 1 ± 2 (0-5) 2 ± 2 (0-5) 2 ± 1 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) Right upper lobe 0 ± 1 (0-3) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) Right middle lobe 0 ± 1 (0-3) 1 ± 1 (0-5) 1 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 1 (0-5) 0 ± 1 (0-5) Right lower lobe 1 ± 2 (0-12) 2 ± 1 (0-5) 2 ± 2 (0-5) 1 ± 2 (0-5) 1 ± 1 (0-5) 1 ± 1 (0-5) Quantitative data were presented as mean ± SD (minimum-maximum). Mann-Whitney U test showed significant difference between Stage-1 and Stage-2-6 (P < 0.05); Wilcoxon test showed significant difference between the left lower lobe and the left upper lobe (P < 0.05); Wilcoxon test showed significant difference between the right lower lobe and the right upper/middle lobe (P < 0.05). 1 P < 0.05 by t test. CT: Computed tomography. and thrombosis, potentially leading to poor outcome [4, [19] [20] [21] [22] . SOFA score is a good diagnostic marker for renal function, and reflects the state and degree of multiorgan dysfunction [20, 23 ]. In the current study, higher PCT and LDH levels were independently associated with prognosis of COVID-19. Additionally, we found that most patients had lower white blood cell count, and no bacterial pathogens were detected. Viral infections is one of the cause of sepsis syndrome, despite that bacterial infections are used to be the primary cause of sepsis, PCT, as an inflammatory indicator, could better stratify the degree of infection. The level of LDH is important in assessing the risk of cardiac and liver dysfunction, which has great significance for both patient isolation decision-making and guidance around the length of antiviral treatment. Effective antiviral therapy may improve the outcome of COVID-19 in spite of that we did not observe a reduction in viral shedding time after antiviral therapy in the current study. However, further research is needed to investigate the pathogenesis of sepsis in COVID- 19 . We showed that CT stage is a powerful indicator in the evaluation of COVID-19 prognosis. We characterized specific factors-prognostic factors model (PFM), age, high-sensitivity CRP, DD, LDH and PCT as a valuable independent prognostic tool of COVID-19 from CT stage. Predictive value of PFM was comparable to that of CT stage. Thus, these results consistently point to the notion that high PFM and CT stage are pivotal factors in evaluating COVID-19, but further research is needed to investigate the prognostic value. Table 4 Distribution and frequency of the major of lung lesions on computed tomography in different stages defined by the time of onset of symptoms The counting data were presented as count (percentage of the total). GGO: Ground glass opacity. Hazard ratios are reported per 1-SD increase in loge D-dimer (DD) and Procalcitonin levels; 1-SD higher loge DD and Procalcitonin was approximately equivalent to 2-fold higher DD and Procalcitonin levels. Model 1: adjusted age, sex. Model 2: adjusted model 1 plus smoking status, history of hypertension, diabetes, cancers, cardiac disease, and chronic pulmonary disease, systolic blood pressure, fasting blood glucose, total cholesterol, triglyceride, white blood cell count, C-reactive protein, creatinine, DD, LDH, procalcitonin. CRP: Creactive protein; DD: D-dimer; LDH: Lactate dehydrogenase; NLT: Natural log-transformation; CI: Confidence interval. However, no published works were found about the dose-response relationship between mortality and severe illness in adult patients with COVID-19. In recent studies, the relationship of prognostic factors with risk of COVID-19 incidence has not been reported. Of note, we found that higher LDH, DD April 6, 2022 Volume 10 Issue 10 and PCT levels were independently associated with a dose-response increased mortality risk in patients with COVID-19. Notably, the dose-response relationship between LDH and PCT levels and incidence of COVID-19 was seen in survivors and patients with severe illness. To our knowledge, this is the first study to demonstrate that the higher risk of COVID-19 incidence associated with LDH and PCT levels provides evidence of the dose-response relationship. Several potential mechanisms might explain the association between LDH, DD and PCT levels and COVID-19 [4, 19, 20, 22, 23] . Although the underlying pathophysiological mechanisms are unclear, it is possible that the presence of COVID-19 risk factors could cover up the effect of LDH, DD and PCT on the risk of COVID-19 among high-risk persons and leave the pernicious effects prominent in relatively healthy adults. Further studies should be performed, which is the key for the development of specific inhibitors targeting COVID-19. There are some limitations to our study. First, contact histories and laboratory testing records for some cases were incomplete. Second, we could only estimate the incubation period in patients who have recorded information. Uncertainty about the exact date (recall bias) might have inevitably influenced our assessment. Third, since our study did not include patients with mild illness who did not seek medical attention, the case fatality rate would likely have been lower in real-world situations. Meanwhile, during the beginning of the pandemic, we a little about COVID-19, so the treatment regimens have been improving. Also, due to limited medical resources, older patients and patients with serious symptoms may have been preferentially admitted, and this may have resulted in bias. Fourth, Figure 2 Adjusted hazard ratios (solid lines) and 95% confidence interval (dashed lines) for coronavirus disease 2019 mortality from restricted cubic splines in a multivariate-adjusted Cox proportional hazard model. The model was adjusted for age, gender, smoking status, history of hypertension, diabetes, cancers, cardiac disease, and chronic pulmonary disease, systolic blood pressure, fasting blood glucose, total cholesterol, triglyceride, white blood cell count, CRP, creatinine, DD, LDH, procalcitonin. CRP: C-reactive protein; LDH: Lactate dehydrogenase; DD: D-dimer. data generation was clinically driven and not systematic. Lastly, this was a retrospective study. To our knowledge, this is the largest retrospective cohort study of COVID-19 patients who have experienced clear results and systematically explored almost all potential risk factors associated with mortality and severe illness. Six stages of lung involvement could be more accurately defined to evaluate the prognosis of COVID-19. The combination of PFM and six stages could provide the rationale for testing novel coronavirus management to improve outcomes. We found that older age, higher LDH and creatinine, and elevated PCT and DD at admission were risk factors for death of patients with COVID-19. These findings suggested that higher LDH, DD and PCT levels were independently associated with increased risk of COVID-19 incidence. Higher LDH, DD and PCT levels were independently associated with a dose-response increased risk of COVID-19 mortality. April 6, 2022 Volume 10 Issue 10 Dose-response assessments and risk factors for mortality, severe cases and clinical outcomes for coronavirus disease 2019 (COVID- 19) have not been well described. To screen for dose-response relationships between risk factors and incidence of COVID-19. To explore risk factors of in-hospital death and describe the clinical course of symptoms, viral shedding, and temporal changes of laboratory findings during hospitalization. This retrospective cohort study included two cohorts of adult inpatients from two designated hospitals. Multivariate logistic regression and Cox proportional risk models were used to determine the doseresponse relationship between risk factors and the incidence of COVID-19. D-dimer and procalcitonin were log-linear correlated with the risk of death from COVID-19, while there was a linear dose-response relationship between age, LDH, D-dimer and procalcitonin, independent of identified risk factors. High lactate dehydrogenase, D-dimer and procalcitonin levels were independently associated with an increased dose-response risk of death from COVID-19. This study provides ideas and basis for prospective observation of dose-response relationships between risk factors and incidence of COVID-19. Author contributions: Zhao SC and Yu XQ made equal contributions to the article; Yu XQ and Zhu Q had the idea for and designed the study and had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis; Zhao SC and Lai XF drafted the paper; Zhao J and Guo DL did the analysis, and all authors critically revised the manuscript for important intellectual content and gave final approval for the version to be published; all authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Institutional review board statement: The study was reviewed and approved by the Zhongnan Hospital of Wuhan Informed consent statement: All study participants or their legal guardian provided informed written consent about personal and medical data collection prior to study enrolment. There has no conflict of interest of this study. The datasets used or analysed during the current study are available from the corresponding author on reasonable request. The authors have read the STROBE Statement-checklist of items, and the manuscript was prepared and revised according to the STROBE Statement-checklist of items. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia China Novel Coronavirus Investigating and Research Team. A Novel Coronavirus from Patients with Pneumonia in China Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study Clinical characteristics and co-infections of 354 hospitalized patients with COVID-19 in Wuhan, China: a retrospective cohort study Genomic characterization of a novel SARS-CoV-2 Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan, China The positive impact of lockdown in Wuhan on containing the COVID-19 outbreak in China National Clinical Research Center for Child Health and Disorders and Pediatric Committee of Medical Association of Chinese People's Liberation Army. A contingency plan for the management of the 2019 novel coronavirus outbreak in neonatal intensive care units Therapeutic and triage strategies for 2019 novel coronavirus disease in fever clinics Imaging of pulmonary viral pneumonia Fleischner Society: glossary of terms for thoracic imaging Pulmonary function and COVID-19 Thin-Section Computed Tomography Manifestations During Convalescence and Long-Term Follow-Up of Patients with Severe Acute Respiratory Syndrome (SARS) Dose-response analyses using restricted cubic spline functions in public health research Epidemiological and clinical features of pediatric COVID-19 Primary to tertiary COVID-19 transmission in a hospital -A cluster outbreak analysis COVID-19 pneumonia: the great radiological mimicker The immunopathogenesis of sepsis in elderly patients Sepsis and Septic Shock Strategies Beneficial Effects of Vaccination on Cardiovascular Events: Myocardial Infarction, Stroke, Heart Failure Cardiac complications in patients with community-acquired pneumonia: incidence, timing, risk factors, and association with short-term mortality Wang H, Kang X, Shi Y, Bai ZH, Lv JH, Sun JL, Pei HH. SOFA score is superior to APACHE-II score in predicting the prognosis of critically ill patients with acute kidney injury undergoing continuous renal replacement therapy.