key: cord-1006322-dbegm2h4
authors: Liu, Zeming; Li, Jinpeng; Chen, Danyang; Gao, Rongfen; Zeng, Wen; Chen, Sichao; Huang, Yihui; Huang, Jianglong; Long, Wei; Li, Man; Guo, Liang; Wang, Xinghuan; Wu, Xiaohui
title: Dynamic Interleukin-6 Level Changes as a Prognostic Indicator in Patients With COVID-19
date: 2020-07-17
journal: Front Pharmacol
DOI: 10.3389/fphar.2020.01093
sha: 0ceb60fcf5e46e1f48e02c44d0cc25f88d366e62
doc_id: 1006322
cord_uid: dbegm2h4

BACKGROUND: Interleukin-6 (IL-6), a proinflammatory cytokine, has been reported to be associated with disease severity and mortality in patients with coronavirus disease 2019 (COVID-19). Yet, dynamic changes in IL-6 levels and their prognostic value as an indicator of lung injury in COVID-19 patients have not been fully elucidated. OBJECTIVE: To validate whether IL-6 levels are associated with disease severity and mortality and to investigate whether dynamic changes in IL-6 levels might be a predictive factor for lung injury in COVID-19 patients. METHODS: This retrospective, single-center study included 728 adult COVID-19 patients and used data extracted from electronic medical records for analyses. RESULTS: The mortality rate was higher in the elevated IL-6 group than in the normal IL-6 group (0.16 vs 5%). Cox proportional hazards and logistic regression analyses for survival (adjusted hazard ratio, 10.39; 95% confidence interval [CI], 1.09–99.23; p = 0.042) and disease severity (adjusted odds ratio, 3.56; 95% CI, 2.06–6.19; p < 0.001) revealed similar trends. Curve-fitting analyses indicated that patient computed tomography (CT) scores peaked on days 22 and 24. An initial decline in IL-6 levels on day 16 was followed by resurgence to a peak, nearly in tandem with the CT scores. CONCLUSION: Increased IL-6 level may be an independent risk factor for disease severity and in-hospital mortality and dynamic IL-6 changes may serve as a potential predictor for lung injury in Chinese COVID-19 patients. These findings may guide future treatment of COVID-19 patients.

The outbreak of coronavirus disease , first detected in Wuhan, China, in December 2019, has been categorized as a pandemic since March 2020. As COVID-19 continues to spread worldwide, 5,722,859 cases and 356,279 deaths have been reported as of May 28, 2020.

The Leishenshan (Thunder God Mountain) Hospital, the largest makeshift hospital in China, played an essential role in the battle against COVID-19 in Wuhan, China. The hospital had the capacity to accommodate 1,600 patients and was well equipped with intensive care units and operating rooms as well as clinical laboratories providing relatively standardized and systematic medical information on COVID-19 patients.

Interleukin 6 (IL-6) and other components of the inflammatory cascade contribute to host defense against infections. However, excessive synthesis of IL-6 can lead to a severe acute systemic inflammatory response known as a "cytokine storm," which confers increased risks of vascular hyperpermeability, multiorgan failure, and eventually death (Meduri et al., 1995) . This pathophysiological process, which is supported by many studies, suggests that higher IL-6 levels may be an important predictor of COVID-19 severity (Cai et al., 2020; Chen et al., 2020) . However, the role of IL-6 as an independent predictor of severity and mortality in hospitalized COVID-19 patients has not been validated.

Acute respiratory distress syndrome (ARDS), one of the major factors contributing to mortality in COVID-19 patients, is similar to severe community-acquired pneumonia caused by other viruses or bacteria (Ranieri et al., 2012; Jose and Manuel, 2020) . The pathophysiological process of ARDS involves a cytokine storm caused by the overproduction of early response proinflammatory cytokines. Theoretically, monitoring the IL-6 level is crucial to understanding the pathophysiological process of ARDS; however, whether IL-6 levels can predict the severity of pulmonary injury has not been fully investigated in patients with Thus, in this study, we aimed to determine whether IL-6 levels are associated with disease severity and mortality in COVID-19 patients and investigate whether dynamic changes in IL-6 levels are predictive of pulmonary lesions caused by ARDS in hospitalized COVID-19 patients.

This retrospective study included 1,880 confirmed COVID-19 patients who were admitted to the Leishenshan Hospital in Wuhan, China, between February 9 and April 9, 2020. Exclusion criteria were patients: younger than 18, missing data on IL-6 levels, acute myocardial infarction during hospitalization, acute pancreatitis, or transferred to another hospital. The dataset included in the final analysis consisted of 728 cases ( Figure S1 ).

Detailed patient information, including demographic characteristics, history of comorbidities, symptoms, laboratory findings, computed tomography (CT) findings, length of hospital stay, oxygen support, illness severity on admission, and treatment, was collected from the electronic patient records stored in the hospital information management system and entered into a customized Excel form (Microsoft, Redmond, WA, USA).

The primary outcomes were patient survival and COVID-19 severity during hospitalization. Based on the 7 th Interim Guidance for the Diagnosis and Treatment of COVID-19 published by the National Health Commission of China, each patient was classified as mild, common, severe, or critical.

As chest CT has a high sensitivity for the diagnosis of COVID-19, studies have advocated for the use of CT images to evaluate the severity of acute respiratory syndrome (Wong et al., 2004; Ai et al., 2020; Lyu et al., 2020) . Furthermore, radiographic scores have been shown to play an important role in the diagnosis of severe acute respiratory syndrome (SARS) between November 2002 and August 2003 (Wong et al., 2004) . In our study, chest CT images were independently inspected by two radiologists who were blinded to the patients' clinical data. An semiquantitative scoring system, based on previous studies and the characteristics of COVID-19, was adopted for the assessment of pulmonary inflammation. Score 1 was generated based on CT imaging features such as ground-glass opacities (GGOs), reticulate or cord-like changes, consolidation, and pleural effusion. Each imaging feature received 1 point. While Score 1 comprised the sum of the points, Score 2 was generated based on the area of lung lobe involvement and graded from 0 to 4 as follows: no involvement, 0; <25% involvement, 1; 26-50% involvement, 2; 51-75% involvement, 3; and 76-100% involvement, 4. The overall score was the sum of Scores 1 and 2.

Patients were assigned to either the elevated or normal IL-6 group based on an IL-6 cutoff value of 7 pg/ml (obtained from medical records). Continuous variables are presented as medians and interquartile ranges (IQRs), whereas categorical variables are expressed as the number of patients (-percentages).

The Kruskal-Wallis or Mann-Whitney U test was used for intergroup comparisons between the elevated IL-6 and normal IL-6 groups with regard to the means of continuous variables. For categorical variables, intergroup comparisons were performed using the chi-squared test or Fisher's exact test.

Cox proportional hazards regression models were used to determine the effect of IL-6 levels on overall survival. We used a multivariate Cox proportional hazards regression model Abbreviations: ARDS, acute respiratory distress syndrome; AUC, area under the curve; CI, confidence interval; COVID-19, coronavirus disease 2019; CT, computed tomography; GGO, ground-glass opacity; HR, hazard ratio; Ig, immunoglobulin; IL-6, interleukin-6; INR, international normalized ratio; MERS-CoV, Middle East respiratory syndrome coronavirus; OR, odds ratio; ROC, receiver-operating characteristic; SARS, severe acute respiratory syndrome coronavirus; SARS-CoV-2, severe acute respiratory syndrome coronavirus 2. adjusted for variables such as age, history of cardiovascular disease, lymphocyte count, D-dimer level, and lactate dehydrogenase level. Survival analysis was conducted using Kaplan-Meier survival curves and the log-rank test or Mantel-Haenszel test. Moreover, the cumulative hazard function for COVID-19 progression in both groups was analyzed. The relationship between IL-6 levels and CT scores of lung injury was assessed using a curve-fitting analysis. Receiver-operating characteristic (ROC) curve analysis was used to evaluate the predictive ability for death and disease severity. All statistical analyses were performed using R version 3.6.1 (https://www.r-project.org/, R Core Team, Vienna, Austria) and EmpowerStats version 2.0-(http://www.empowerstats.com/cn/, X&Y Solutions, Inc., Boston, MA). A two-sided p-value <0.05 was considered statistically significant.

The study protocol was approved by the Research Ethics Committee of the Zhongnan Hospital of Wuhan University (approval no.: 2020074). The Research Ethics Committee waived the need for informed consent from the patients due to the urgent need for research insights in the context of this rapidly evolving infectious disease. All procedures undertaken in this study were in accordance with the ethical standards of the institutional and national research committee as well as the 1964 Declaration of Helsinki and its later amendments and other comparable ethical standards.

The baseline characteristics of 728 patients are listed in Table 1 . The median ages in the normal and elevated IL-6 groups were 57 and 68 years, respectively. A total of seven (0.96%) patients died. In the normal IL-6 group, the majority of patients (54.4%) were female. In contrast, there was a higher proportion (54.2%) of male patients in the elevated IL-6 group. Cardiovascular diseases were more prevalent in the elevated IL-6 group (44.17%) than in the normal IL-6 group (20.23%); similarly, neurological diseases were also more common in the elevated IL-6 than in the normal IL-6 group (11.76 vs. 4.93%). As shown in Table 2 , the levels of routine blood test parameters, including leukocyte count, neutrophil count, were significantly higher in the elevated IL-6 group than in the normal IL-6 group. Moreover, the lymphocyte, erythrocyte, and platelet counts and hemoglobin concentration were reduced in the elevated IL-6 group. With regard to liver function parameters, the direct bilirubin level was higher and accompanied by notable reductions in albumin and total protein levels in the elevated IL-6 group. An evaluation of renal function showed significantly elevated serum creatinine, urea, and cystatin C levels. Lactate dehydrogenase levels were increased in the elevated IL-6 group. Significant abnormalities in blood coagulation parameters, such as prolonged prothrombin time and activated partial thromboplastin time, higher International Normalized Ratio, and increased D-dimer levels, were noted in the elevated IL-6 group. Moreover, the proportion (84.62%) of patients with positivity for immunoglobulin G (IgG) was lower in the elevated IL-6 group than in the normal IL-6 group (94.03%). The clinical treatment and outcome analyses revealed that more patients in the elevated IL-6 group were likely to require antibiotics, antiviral drugs, anticoagulants, corticosteroids, and vitamin C therapy. Furthermore, in the normal IL-6 group, there was a higher use of antimalarial drugs (83.93%) and traditional Chinese medicine (88.82%). Patients in the elevated IL-6 group were more likely to be admitted to the intensive care unit during hospitalization (20.00%) than patients in the normal IL-6 group (0.82%) ( Table 3) . 

Levels Normal IL-6(n = 608), n (%) Elevated IL-6(n = 120), n (%) All patients(n = 728), n (%) P value 

The Kaplan-Meier curves and cumulative hazard function showed that patients with normal IL-6 levels had better survival (p < 0.05, Figures 1 and 2) . Further, the ROC analysis showed that the elevated IL-6 group had the largest areas under the ROC curve (AUCs) of 0.870 and 0.778 for disease severity and CT scores, respectively ( Figure 3 ).

The peak value of Score 1 for all patients was 2.35 on day 23 ( Figure 4A ); however, while the peak value for the elevated IL-6 group was 2.46 on day 23, the peak value in the normal IL-6 group was 2.24 on day 22 ( Figure 4B) . Nonetheless, the peak value of Score 1 was attained within the same period in both groups. The Score 2-fitted curve for all patients presented a linear line that indicated a decline over time ( Figure 4C ). In the normal IL-6 group, the fitted curve presented a more uniform declining linear line than that in the elevated IL-6 group ( Figure 4D ). For the overall score, all patients showed peak values on day 22 (4.60; Figure 4E ). Patients in the elevated IL-6 group showed the peak score on day 24 (4.79; Figure 4E ), whereas a fitted curve that presented a flattened, declining, linear line was seen for patients in the normal IL-6 group. The curve-fitting analysis was used to evaluate the changes in IL-6 levels: there was a slight decrease on day 16 and an increase to the peak level on day 26 ( Figure 5A) . A similar trend and an obvious dynamic change in IL-6 levels were observed in the elevated IL-6 subgroup ( Figure 5B ).

Using data from a unique makeshift hospital for COVID-19 patients, this results from this retrospective study showed that higher serum levels of IL-6 was an independent and reliable risk factor for COVID-19 patients and led to higher disease severity 

Levels Normal IL-6 (n = 608) n (%) Elevated IL-6 (n = 120) n (%) and mortality. Moreover, the trend fluctuation occurred earlier than the negative changes in lung CT scores, which indicated a potential predictive ability for lung injury. These findings provide clinical evidence in support of recently published guidance statements by Aziz et al. (2020) . Nevertheless, due to the programmatic, retrospective nature of this research, the inferences derived from the data require further validation from cohort studies and randomized controlled trials (RCTs). IL-6, a chemokine secreted by T-cells and macrophages to stimulate the immune response, is an important biomarker of inflammation, and elevated IL-6 levels have been demonstrated in inflammatory states induced by several pathological conditions (Cai et al., 2020) . In COVID-19 patients, overproduction of cytokines, such as IL-6, might activate coagulation pathways, with a resultant disruption of procoagulant-anticoagulant homeostasis, induction of disseminated intravascular coagulation, and multiorgan dysfunction or failure (Jose and Manuel, 2020) . In this study, we found that COVID-19 patients with elevated IL-6 levels had a higher likelihood of experiencing renal, hepatic, and lung injuries ( Table 2) .

Previous clinical studies have shown that older age, male sex, cardiovascular diseases, lower lymphocyte counts, higher Ddimer levels, and increased lactate dehydrogenase levels are significant risk factors for illness severity or mortality in COVID-19 patients Wang et al., 2020) . In this study, there were a larger proportion of patients with these in the elevated IL-6 group. Thus, patients with high IL-6 levels may have a worse prognosis. Research has demonstrated that a higher IL-6 level on admission can predict greater odds of complications and illness severity in COVID-19 patients ; however, mortality was not evaluated. Although a few studies have indicated that increases in IL-6 levels are potentially associated with an increased likelihood of mortality, they lacked statistical power as residual confounders were not considered or only indirectly validated through lower mortality after using an interleukin-6 receptor antagonist for blockading cytokine release in severe COVID-19 patients (Alattar et al., 2020; Aziz et al., 2020; Zhang et al., 2020) . A multivariate Cox regression model was developed in our study to adjust for confounders including older age, history of cardiovascular disease, lower lymphocyte count, higher D-dimer levels, and increased lactate dehydrogenase levels. These risk factors of mortality have been validated in previous research (Ji et al., 2020; Zhou et al., 2020) . Additionally, in our study, Kaplan-Meier and cumulative hazard analyses were conducted to verify the strong associations of IL-6 levels with disease severity and mortality. Thus, all analyses conducted in our study validated that elevated IL-6 levels are an independent risk factor for the severity and mortality of COVID-19.

Furthermore, we evaluated the potential predictive function of IL-6. COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) , which binds to the alveolar epithelial cells to subsequently activate innate and adaptive immune responses and induces the release of a large number of cytokines, including IL-6. Moreover, considering the inherent role of these proinflammatory factors, vascular permeability increases, followed by the influx of large amounts of fluid and a large number of blood cells into the alveoli, resulting in dyspnea or ARDS (Knudsen and Ochs, 2018; Leiva-Juarez et al., 2018) . The flux in IL-6 levels is associated with the pathological processes that mediate lung injuries, which was demonstrated by the similarity in peak timings of IL-6 levels and CT scores (Figure 4 ) in this study. Interestingly, IL-6 levels showed a slight decrease prior to the increase in CT scores in the fitting curve. This may be a result of SARS-CoV-2 eliciting attenuated innate immunity, leading to delayed proinflammatory cytokine induction, similar to the reaction observed in patients infected with Middle East respiratory syndrome coronavirus (MERS-CoV) (Lau et al., 2013) . Further, in the early phase of the cytokine storm, immune cells and cytokines may be clustered in the pulmonary lesions, resulting in a decrease in IL-6 levels in the peripheral blood. Consequently, we observed a mild decline occurred in IL-6 serum levels over the initial 10-20 days of the clinical course ( Figure 4) . In the management of SARS patients, dynamic changes in IL-6 levels have been shown to occur synchronously with the changes in radiographic scores (Chien et al., 2006) . Moreover, a similar association between IL-6 levels and CT scores was observed in COVID-19 patients over a period of 20-30 days ( Figure 5) . The decrease in IL-6 levels occurred earlier than the increase in CT scores that reflected the progressive deterioration of lung lesions ( Figure 5 ). The disparity in the two aforementioned scores was more obvious in the elevated IL-6 group. Thus, the decreasing trend of IL-6 levels might be predictive of the imminent worsening of lung injuries. In addition, the good predictive capability of IL-6 for disease severity and prognosis in COVID-19 patients, assessed by the ROC curves (AUC > 0.75; Figure 3 ), provided further evidence regarding its predictive function. Nonetheless, this study has several limitations. First, due to the programmatic and retrospective nature of this study, we could not measure serum IL-6 levels or assess CT images on the same day. Second, the size of the study sample was modest, and there was inadequate statistical power to detect whether there was a differential effect among the mild, severe, and critical COVID-19 cases. Third, as this was a retrospective study, data on all parameters were not available for all patients, and data regarding in-hospital medications may not have been fully recorded. Some patients were transferred from other hospitals and might not have received uniform treatment in the early stages of COVID-19. Fourth, testing for IL-6 levels was performed at different time intervals for each patient. Hence, diagnostic bias may have resulted from the increased number of tests performed in patients with high IL-6 levels. Fifth, we were unable to retrieve pre-hospital self-medication data from the in-hospital electronic records, considering the emergent circumstances of the COVID-19 pandemic. Sixth, median age was different between the elevated IL-6 group and the normal IL-6 group that may have caused bias in the analyses of different laboratory parameters, although we adjusted for age in our multivariate analysis. Finally, the use of some medications, such as immunomodulators used for suppressing an overactive cytokine response may have influenced the experimental results. There may have been unknown interactive effects between medications and IL-6.

The results from this study need to be interpreted considering other potential and residual confounders. It is likely that a higher IL-6 level is an independent risk factor for in-hospital severity and Frontiers in Pharmacology | www.frontiersin.org July 2020 | Volume 11 | Article 1093 mortality in Chinese COVID-19 patients, and it might be a potential predictor of lung injury in these patients. The findings of this study provide novel research insights that may guide the early therapeutic intervention for COVID-19 patients, such as dietary interventions (Messina et al., 2020) and medication with IL-6 receptor blockers . However, cohort studies, prospective studies, and RCTs are required to further validate the impact of these findings and the early therapeutic strategies.

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

The studies involving human participants were reviewed and approved by Research Ethics Committee of the Zhongnan Hospital of Wuhan University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

ZL, JL, DC, RG, and WZ designed the study, collected and analyzed the data, and wrote the manuscript. SC, YH, JH, WL, and ML collected and reviewed the clinical, laboratory, and radiological data. ZL, JL, and DC performed statistical analyses. RG and WZ reviewed, interpreted, and checked the clinical data.

LG, XWa, and XWu wrote the manuscript and provided valuable suggestions for the study design and data analysis.

LG, XWa, and XWu contributed equally, designed the project, edited the manuscript, and supervised the study. All authors contributed to the article and approved the submitted version.

We would like to thank all patients and their families. As Wuhan was hit hard by the pandemic not long ago, its people made huge sacrifices to push back the virus. We also would like to thank the doctors, nurses, and other medical team members from different places across China. In our most difficult days, they chose to stay with us at the Leishenshan Hospital, offered a helping hand, and helped us make it through.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar. 2020.01093/full#supplementary-material FIGURE S1 | Design schematic of the study.

Correlation of Chest CT and RT-PCR Testing in Coronavirus Disease 2019 (COVID-19) in China: A Report of 1014 Cases

Tocilizumab for the Treatment of Severe COVID-19

Elevated Interleukin-6 and Severe COVID-19: A Meta-Analysis

A B FIGURE 5 | Comparison of CT scores trend and IL-6 levels over time with curve fitting analyses. (A) All the COVID-19 patients

Interleukin-6 Level of COVID-19 Patients Frontiers in Pharmacology | www

COVID-19 in a designated infectious diseases hospital outside Hubei Province

Clinical and immunological features of severe and moderate coronavirus disease 2019

Temporal changes in cytokine/chemokine profiles and pulmonary involvement in severe acute respiratory syndrome

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China

Prediction for Progression Risk in Patients with COVID-19 Pneumonia: the CALL Score

COVID-19 cytokine storm: the interplay between inflammation and coagulation

The micromechanics of lung alveoli: structure and function of surfactant and tissue components

Delayed induction of proinflammatory cytokines and suppression of innate antiviral response by the novel Middle East respiratory syndrome coronavirus: implications for pathogenesis and treatment

Lung epithelial cells: therapeutically inducible effectors of antimicrobial defense

Risk factors for severity and mortality in adult COVID-19 inpatients in Wuhan

Prognostic value of interleukin-6, C-reactive protein, and procalcitonin in patients with COVID-19

The performance of chest CT in evaluating the clinical severity of COVID-19 pneumonia: identifying critical cases based on CT characteristics

Inflammatory cytokines in the BAL of patients with ARDS. Persistent elevation over time predicts poor outcome

Functional Role of Dietary Intervention to Improve the Outcome of COVID-19: A Hypothesis of Work

Acute respiratory distress syndrome: the Berlin Definition

15-day mortality and associated risk factors for hospitalized patients with COVID-19 in Wuhan, China: an ambispective observational cohort study

Severe acute respiratory syndrome: thin-section computed tomography features, temporal changes, and clinicoradiologic correlation during the convalescent period

The cytokine release syndrome (CRS) of severe COVID-19 and Interleukin-6 receptor (IL-6R) antagonist Tocilizumab may be the key to reduce the mortality

Clinical course and risk factors for mortality of adult inpatients with COVID-19 in Wuhan, China: a retrospective cohort study

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.Copyright © 2020 Liu, Li, Chen, Gao, Zeng, Chen, Huang, Huang, Long, Li, Guo, Wang and Wu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.