Xiangyun Teng1†
Jiaqi Zhang1,2†Yaling Shi3†Yuntao Liu4†Yanqing Yang5Jinyong He1Shuhong Luo5Yile Huang1Yanxia Liu3Dongdong Liu2Yizhe Li1Shuangzhe Zhang5Ruo-Pan Huang5,6*Dawei Wang7*
Jianhua Xu1*- 1Department of Laboratory Medicine, Shunde Hospital of Guangzhou University of Chinese Medicine, Foshan, China
- 2Department of Laboratory Medicine, The Second Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, China
- 3Department of Laboratory Medicine, Guangzhou Eighth People’s Hospital, Guangzhou, China
- 4Emergency Department, The Second Affiliated Hospital of Guangzhou University of Chinese Medicine, Guangzhou, China
- 5Research and Development Department, RayBiotech, Inc., Guangzhou, China
- 6Raybiotech Center, RayBiotech, Inc., Norcross, GA, United States
- 7Department of Pulmonary and Critical Care Medicine, Shunde Hospital of Guangzhou University of Chinese Medicine, Foshan, China
To systematically explore potential biomarkers which can predict disease severity in COVID-19 patients and prevent the occurrence or development of severe COVID-19, the levels of 440 factors were analyzed in patients categorized according to COVID-19 disease severity; including asymptomatic, mild, moderate, severe, convalescent and healthy control groups. Factor candidates were validated by ELISA and functional relevance was uncovered by bioinformatics analysis. To identify potential biomarkers of occurrence or development of COVID-19, patient sera from three different severity groups (moderate, severe, and critical) at three time points (admission, remission, and discharge) and the expression levels of candidate biomarkers were measured. Eleven differential factors associated with disease severity were pinpointed from 440 factors across 111 patients of differing disease severity. The dynamic changes of GDF15 reflect the progression of the disease, while the other differential factors include TRAIL R1, IGFBP-1, IGFBP-4, VCAM-1, sFRP-3, FABP2, Transferrin, GDF15, IL-1F7, IL-5Rα, and CD200. Elevation of white blood cell count, neutrophil count, neutrophil-lymphocyte ratio (NLR), Alanine aminotransferase and Aspartate aminotransferase, low lymphocyte and eosinophil counts in the severe group were associated with the severity of COVID-19. GDF15 levels were observed to be associated with the severity of COVID-19 and the dynamic change of GDF15 levels was closely associated with the COVID-19 disease progression. Therefore, GDF15 might serve as an indicator of disease severity in COVID-19 patients.
Introduction
Coronavirus disease 2019 (COVID-19) resulting from severe acute respiratory syndrome coronary virus 2 (SARS-CoV-2) infection has been sweeping the globe at an unprecedented rate and is producing profound but still unfolding health and socio-economic impacts (1–4). Clinical experience to date has shown a high heterogeneity of COVID-19, with disease severity ranging from asymptomatic to mild to severe, resulting in a wide spectrum of clinical manifestations (cough, fever, myalgia, and malaise) that can eventually lead to acute respiratory distress syndrome (ARDS) or even cause death (5–7). However, it has been reported that many patients that became critically ill and died did not show severe clinical manifestations in the early stages of the disease, instead presenting with common symptoms such as mild fever, cough, or muscle soreness. These patients’ conditions suddenly deteriorated in the later stages of the disease or during recovery, where multiple organ failure, acute respiratory distress syndrome (ARDS), and death occurred within a very short period of time (2). Therefore, the research of biomarkers to monitor the COVID-19 process is particularly important.
Several studies have reported that cytokines are correlated with the progression of COVID-19. Cytokine synthesis, innate immune cell activation and adaptive immune response are key for rapidly targeting and combating the virus (8). However, uncontrolled excessive inflammation may lead to cytokine release syndrome (CRS), with impaired function causing tissue damage, immune cell apoptosis, and cytotoxicity (9). Cytokine storm is considered to be one of the major causes of ARDS and multiple-organ failure (10). Patients with severe COVID-19 have been reported to display unbalanced immune regulation (11, 12). Pro-inflammatory cytokines such as interleukin-1 (IL-1), IL-6, IL-8, and tumor necrosis factor alpha (TNF-α) have been found to be elevated in the sera of severely and critically ill COVID-19 patients, anti-inflammatory cytokines such as IL-10 have also been found in the sera of COVID-19 patients (13–15). Plasma IP-10, MCP-3 sTIM-3, HGF, IL-1α, IL-27, and PCT cytokines have also been reported to be related to the severity of the disease (16–19). The profile of cytokine gene expression by cDNA arrays and protein arrays are widely used in different infectious diseases to identify biomarkers (20–22).
To systematically explore potential severity-associated factors of COVID-19, we comprehensively analyzed the expression of 440 factors in asymptomatic, mild, moderate, severe, convalescent and healthy control groups (categorized according to COVID-19 disease severity). These factors encompass, cytokines, chemokines, growth factors, and secreted enzymes, etc., which relate to vascular growth, cell proliferation, cell migration, apoptosis, the inflammatory response as well as other processes. To further understand the dynamic changes of differentially expressed factors during the progression of COVID-19 in patients, sera were collected from patients at three time points (admission to the hospital, remission of symptoms, and discharge from the hospital) and the levels of the 440 factors were measured. To the best of our knowledge, this was the first study to exclusively characterize the factors associated with the severity and progression of COVID-19 using an antibody array profile. In order to provide a comprehensive dynamic characterization, longitudinal data from complete blood count, C-reaction protein (CRP), and liver function were also considered. We found eleven differentially expressed factors associated with disease severity, specifically TRAIL R1, IGFBP-1, IGFBP-4, VCAM-1, sFRP-3, FABP2, transferrin, GDF15, IL-1F7, IL-5Rα, and CD200. We further investigated GDF15 serum levels, and the results suggested that this cytokine may serve as an effective indicator of disease severity in COVID-19 patients.
Materials and Methods
Research Design and Patients
From 22 January 2020 to 13 May 2020, we carried out a retrospective study analyzing 111 patients diagnosed with COVID-19. These patients were diagnosed with real-time reverse transcriptase polymerase chain reaction (RT-PCR) oropharyngeal swab or nasopharyngeal swab tests upon admission to the Eighth People’s Hospital of Guangzhou. This hospital is the designated hospital for the treatment of COVID-19 in Guangzhou, and all COVID-19 patients in the current study were from this hospital.
To systematically explore potential severity associated factors, these 111 patients were categorized according to COVID-19 disease severity; including asymptomatic (n=14), mild (n=12), moderate (n=34), severe (n=18), and convalescent (n=33). Healthy controls (n = 20) came from the Shunde Hospital of Guangzhou University of Chinese Medicine. To further investigate the relative levels of factors at different points in the progression of the disease, we tracked 15 of the COVID-19 patients across three different severity groups (moderate n = 10, severe n = 4, and critical n = 1) and observed their statuses at three time points (admission to the hospital, remission of symptoms, and discharge from the hospital). This study was approved by the Ethics Committee of Shunde Hospital, Guangzhou University of Chinese Medicine (approval number: KY2020001).
Classification of Disease Severity
Following the Protocol on Prevention and Control of Novel Coronavirus Pneumonia (Edition 6) (23), asymptomatic cases are those that have a positive COVID-19 RT-PCR test or specific serum IgM antibodies, but do not know that they are infected or show clinically identifiable symptoms. In addition, based on the Guidelines for the Diagnosis and Treatment Protocol for Novel Coronavirus Pneumonia (Trial Version 7) (24), patients with mild cases displayed mild clinical symptoms with no radiological findings of pneumonia. Moderate cases are characterized by fever (temperature above 37.3°C) and respiratory symptoms with radiological findings of pneumonia. Severe cases are those that meet any of the following criteria: respiratory distress (≥ 30 breaths/min), oxygen saturation ≤ 93% at rest, or arterial partial pressure of oxygen (PaO2)/fraction of inspired oxygen (FiO2) ≤ 300 mmHg (l mmHg = 0.133 kPa). Cases with chest imaging showing obvious lesion progression within 24–48 hours of > 50% were also managed as severe cases. Convalescent patients are considered as those with a body temperature back to normal for more than 3 days, obviously improved respiratory symptoms, obvious absorption of inflammation on pulmonary imaging, and two consecutive negative nucleic acid tests with respiratory tract samples such as sputum or nasopharyngeal swabs (sampling interval being at least 24 hours). Those who met the above criteria were eligible for being discharged.
To investigate the dynamic progression of treatment, admission, remission, and discharge timepoints were observed. The admission timepoint includes admission to the hospital within 24 hours. The remission timepoint refers to when the patient’s lung CT images have improved, and symptoms have been alleviated after treatment. The discharge timepoint refers to having a negative nucleic acid (RT-PCR) test and meeting the criterion of Diagnosis and Treatment Protocol for Novel Coronavirus Pneumonia (Trial Version 7) for discharge.
Data Collection
The sex, age, epidemiological history, symptoms and other clinical characteristics of each patient were obtained from the electronic case system of the corresponding hospital. Routine laboratory parameters, including complete blood count, hemoglobin (Hb), C-reaction protein (CRP) and liver function data, were obtained from EDTA-anticoagulated whole blood. Sera were used for liver function, the Q440 factor antibody array, and enzyme-linked immunosorbent assay (ELISA) validation.
Profiling With the Antibody Array
Serum levels of 440 factors were examined using antibody array, the Human Cytokine Array Q440 (QAH-CAA-440; Raybiotech, Inc., Norcross, GA, USA) was used, and each sample was assayed four times. The array included cytokines, chemokines, growth factors, and secreted enzymes, The analysis was performed strictly following the manufacturer’s protocol. The array was dried at room temperature for 2 hours, and 50 μL of serum was then placed on the array and incubated overnight at 4°C. After washing the array, 80 μL of secondary antibody solution was added to the array and incubated for 2 hours at room temperature. The array was washed again, and 80 μL of detection solution was added to the array and incubated for 1 hour at room temperature. Finally, an InnoScan 310 microarray scanner (Innopsys, Carbonne, Midi-Pyrenees, France) was used to scan the microarray and extract the signal. RayBio Q Analyzer software (Raybiotech, Inc., Norcross, GA, USA) was used to convert these fluorescence data to concentration values.
ELISA Validation
To validate the results from the antibody array Q440, the expression levels of IL-10, IP-10, GDF15, eotaxin-2 and MCP-2 were measured using ELISA (Raybiotech, Inc., Norcross, GA, USA). The experiment was carried out according to the manufacturer’s instructions. Briefly, 100 μL of diluted serum was added into each well and incubated for 2.5 hours at room temperature. After washing the wells, 100 µL of solution containing the biotin-labeled antibody was added to each well and incubated for 1 hour at room temperature, followed by incubation of the resulting solution with a Horseradish Peroxidase (HRP)-streptavidin solution for 45 minutes at room temperature. The signal was developed using TMB One-Step Substrate Reagent which was incubated for 30 minutes at room temperature. Next, 50 μL of a stop solution was added and the absorbance at a wavelength of 450 nm was immediately recorded using an Elx800 microplate reader (BioTek Instruments, Inc., headquartered in Winooski, VT, USA). Analytical validation was performed using reference controls and biological replicates.
Statistical Analysis
Categorical variables (demographics and baseline characteristics of patients) were described as frequency rates and percentages, and continuous variables (laboratory parameters) were described using mean, median, and interquartile range (IQR) values. No imputation was made for missing data. Clinical characteristics of COVID-19 patients were analyzed using SPSS 23.0 (IBM Corporation, Armonk, NY, USA). The Mann-Whitney test or the Kruskal-Wallis H test were used to compare the differences among the groups, and the differences were considered statistically significant at P < 0.05.
The R software package limma from R/Bioconductor was used to analyze differentially expressed factors associated with disease severity in COVID-19 patients. This software was specifically used to compare the levels of factors between two different groups using moderated t-statistics. The results included (log2) fold changes and P-values for each factor. Next, FDR (the BH method) was used to calculate adjusted P-values (adj.P.Val). Differentially expressed factors were defined as those with adjusted P value (adj.P.Val) < 0.05, and a fold-change greater than 1.2 or less than 0.83 (absolute logFC > 0.263). A Venn diagram was made online to screen for and show the logical connections between the differentially expressed factors between groups (http://bioinformatics.psb.ugent.be/webtools/Venn/).
The R software package limma was also used to analyze the dynamic changes of differentially expressed factors during progression in COVID-19 patients, using a linear empirical Bayesian model. Intra-individual correlations were included in the analysis. Differential expression was determined as an adjusted P value (adj.P.Val) < 0.10. R software (version 3.6.3; R Foundation for Statistical Computing, Vienna, Austria) was used to statistically analyze all of the differential factor data.
The receiver operating characteristic (ROC) curve was constructed using GraphPad Prism® version 8.4.2 (GraphPad Software™, San Diego, CA, USA) to evaluate the diagnostic values of factors at different time points during COVID-19 disease progression. All corresponding figures were also constructed using this software.
Heat map analysis and correlation analysis (R package “gplots”), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis (R package “clusterProfiler”) and Gene Ontology (GO) enrichment analysis (R package “org.Hs. eg.db” and “clusterProfiler”) were executed using the open source program R (version 3.6.3).
Results
Demographics, Symptoms, and Comorbidities
Demographic and health condition data for the participants are summarized in Table 1. The median age of the non-severe group (asymptomatic, mild, and moderate) was 36 years (IQR, [24–49]). The median age of the asymptomatic group was 25 years (IQR, [23–45]), that of the mild group was 25 years (IQR, [22–41]), and that of the moderate group was 44 years (IQR, [28–58]). The median age of the severe group was 60 years (IQR, [55–72]), which was significantly increased compared with the non-severe group (P < 0.05). Further, 28% (17/60) of the patients in the non-severe group were male, whereas 56% (10/18) of the patients in the severe group were male, which was significantly increased compared with the non-severe groups (P < 0.05). During the course of the disease, coughing (13/18, 72%), fatigue as described by the patients (9/18, 50%), and fever (9/18, 50%) were the three most common three symptoms in the severe group, and the incidence of these symptoms were significantly increased in the severe group compared with in the non-severe groups (P < 0.05). Of the 111 patients with COVID-19, 77% (73/111) had one or more comorbidities, the most common of which were hypertension (32/111, 29%), liver cysts (20/111, 18%), and anemia (15/111, 14%). High blood pressure was the most common comorbidity in the severe group (9/18, 50%), followed by diabetes (4/18, 22%), and these occurrences significantly increased compared to occurrences in the non-severe groups (P < 0.05).
Table 1 Demographics and baseline characteristics of patients infected with SARS-CoV-2.
Laboratory Characteristics of COVID-19 Patients
The routine clinical laboratory characteristics of the participants are shown in Table 2 and Figure 1. In the severe group, the median white blood cell count was 7.74×109/L (IQR, [4.73–10.14]) and the median neutrophil count was 6.02×109/L (IQR, [3.61–8.70]), which was significantly increased compared to the median white blood cell count was 5.60×109/L (IQR, [4.69–6.86]) and the median neutrophil count was 3.72×109/L (IQR, [2.71–4.62]) in the non-severe group (P < 0.05). The median lymphocyte counts in the non-severe group and severe group were 1.42×109/L (IQR, [1.20–1.98]) and 1.08×109/L (IQR, [0.56–1.25]), respectively, whereas the median eosinophil counts in the non-severe group and severe group were 0.04×109/L (IQR, [0.01–0.09]) and 0.00×109/L (IQR, [0–0.015]), respectively. Both lymphocytes and eosinophils were significantly decreased in the severe group compared with the counts in the non-severe group (P < 0.001). We found a trend in the lymphocyte and eosinophil counts in the convalescent group to recovery. The lymphocyte count decreased with increasing disease severity in the mild, moderate and severe groups. The neutrophil-lymphocyte ratio (NLR) was 6.38×109/L (IQR, [3.15–16.55]) in the severe group, which was significantly increased compared with the non-severe group (2.59×109/L, IQR, [1.73-3.38], P < 0.001). Alanine aminotransferase (ALT) levels were 26.10 U/L (IQR, [18.75-49.08]) in the severe group and 15.25 U/L (IQR, [10.75–24.29]) in the non-severe group. Aspartate aminotransferase (AST) levels were 26.95 U/L (IQR, [19.48–51.73]) in the severe group and 16.25 (IQR, [13.15–21.25]) in the non-severe group. Both ALT and AST levels were significantly increased in the severe group compared with the non-severe group (P < 0.001). The C-reactive protein (CRP) expression in the severe group was 42.03 mg/L (IQR, [22.06–68.63]), and was significantly lower in the non-severe group (6.79 mg/L, IQR, [3.37–12.88], P < 0.001). Moreover, we found that the lymphocyte, eosinophil, and platelet counts continuously increased with disease progression and patient status (admission, remission, and discharge), as shown in Table 3 and Figure 2.
Table 2 Laboratory parameters of SARS-CoV-2-infected patients classified by disease severity.
Figure 1 Laboratory parameters from different COVID-19 severity groups (asymptomatic, mild, moderate, severe and convalescent) and healthy controls. Differential levels of routine complete blood count, hemoglobin (Hb), C-reaction protein (CRP) and liver function across different severity groups. White blood cell count, neutrophil count, neutrophil-lymphocyte ratio (NLR), alanine aminotransferase and aspartate aminotransferase levels were higher in the severe group compared with the non-severe groups (asymptomatic, mild, and moderate) (P < 0.05). Lymphocyte and eosinophil counts were significantly decreased in the severe group compared with the non-severe group (P < 0.001), with a trend to recovery in the convalescent group. Healthy control (n=20), asymptomatic group (n=14), mild group (n=12), moderate group (n=34), severe group (n=18), non-severe group (n=60), convalescent group (n=33). Significance is denoted as *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001. WBC, White blood cell count; NEU, Neutrophils count; LYM, Lymphocytes count; NLR, neutrophil-lymphocyte ratio; MON, Monocyte count; EOS, Eosinophil count; BAS, Basophils count; RBC, Erythrocyte count; Hb, Haemoglobin; PLT, Platelet count; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; CRP, C-reactive protein.
Table 3 Laboratory parameters of SARS-CoV-2-infected patients classified by disease/patient status progression.
Figure 2 Laboratory parameters with the progression of COVID-19/patient status (at three time points: admission, remission, and discharge). Different levels of routine complete blood count, hemoglobin (Hb), C-reaction protein (CRP) and liver function were found at different time-points of disease/patient status progression. Moreover, as patients improved, lymphocyte and eosinophil counts consistently increased, whereas NLR and AST consistently decreased. Admission group (n=14), remission group (n=14), discharge group (n=14). Significance is denoted as *P < 0.05, **P < 0.01 and ****P < 0.0001. WBC, White blood cell count; NEU, Neutrophils count; LYM, Lymphocytes count; NLR, neutrophil-lymphocyte ratio; MON, Monocyte count; EOS, Eosinophil count; BAS, Basophils count; RBC, Erythrocyte count; Hb, Haemoglobin; PLT, Platelet count; ALT, Alanine aminotransferase; AST, Aspartate aminotransferase; CRP, C-reactive protein.
Factor Profiling With Antibody Arrays
To investigate the factors associated with COVID-19, we first used an antibody array to screen 440 factors in sera from different disease severity groups (asymptomatic, mild, moderate, severe and convalescent) and with healthy controls. Differentially expressed factors from the 440 factors assessed in each of the different disease severity groups compared to the healthy control group were determined under the screening conditions of adj.P.Val < 0.05 and fold-change greater than 1.2 or less than 0.83 (absolute logFC > 0.263).There were 68 differential factors identified in the asymptomatic group, 65 differential factors identified in the mild group, 106 differential factors identified in the moderate group, 90 differential factors identified in the severe group, and 108 differential factors identified in the convalescent group. We further analyzed and compared taxon overlap among these intergroup differential factors using Venn diagrams (Figure 3A). Finally, a total of eleven differential factors related to COVID-19 disease severity were identified, including: TRAIL R1, IGFBP-1, IGFBP-4, VCAM-1, sFRP-3, FABP2, Transferrin, GDF15, IL-1F7 (also known as IL-37), IL-5Rα, and CD200. 9 of the 11 differentially expressed factors were significantly increased in COVID-19 patients compared with healthy controls: TRAIL R1, IGFBP-1, IGFBP-4, VCAM-1, sFRP-3, FABP2, Transferrin, GDF15, and CD200 (P < 0.01). 2 of the 11 differentially expressed factors: IL-1F7 and IL-5Rα were significantly decreased in COVID-19 patients compared with healthy controls (P < 0.01) (Figure 3B). In addition, we found a consistent increase in GDF15 levels with increasing disease severity among the asymptomatic, mild, moderate, and severe groups of COVID-19 patients. Moreover, GDF15 expression in the convalescent group was recovered to similar levels compared to the healthy control group, suggesting that GDF15 levels closely reflected the progression of the disease.
Figure 3 Eleven differentially expressed factors were highly associated with COVID-19 disease severity. (A) Venn diagram displaying the common differentially expressed cytokines in patients with COVID-19 of different severities. (B) Individual levels of the eleven differentially expressed cytokines in different COVID-19 severity groups. Healthy control (n=20), asymptomatic group (n=14), mild group (n=12), moderate group (n=34), severe group (n=18), convalescent group (n=33). Significance is denoted as *P < 0.05, **P < 0.01, ***P < 0.001, and ****P < 0.0001. HC, Healthy Control; ASY, Asymptomatic; MIL, Mild; MOD, Moderate; SEV, Severe; CON, Convalescent.
To investigate possible independent predictors for the progression of COVID-19, we further analyzed the eleven differentially expressed COVID-19-associated factors over the disease/treatment progression at 3 time points: admission, remission and discharge. Although most of the eleven factors showed dramatic differences in expression between the groups with different disease severity, the correlation between the magnitude of the changes and disease progression was not consistent.
To further validate GDF15 as an independent predictor of the severity of COVID-19, the correlation between GDF15 levels and disease/treatment stage (admission, remission, and discharge) was analyzed. In the moderate group, GDF15 levels increased mildly from admission to remission, indicating that the disease was progressing. However, GDF15 levels decreased markedly from remission to discharge. In the severe group, the levels of GDF15 gradually decreased with recovery progress across admission, remission, and discharge (Figure 4A). Taken together, the changes over time in GDF15 levels reflected the progression of the disease. The diagnostic value of GDF15 was evaluated using area under the curve (AUC) analysis by the receiver operating characteristic (ROC) curve. The AUC value of 0.89 (Figure 4B) indicated the serum GDF15 level to be an efficient diagnostic indicator of disease severity in COVID-19 patients.
Figure 4 GDF15 is an independent predictor of the progression of COVID-19. (A) GDF15 levels in critical (n=1), severe (n=4), and moderate (n=10) groups, each as a function of the progress of the infection/patient status. (B) The AUC value of 0.89 indicated the serum GDF15 level to be an efficient diagnostic indicator of disease severity in COVID-19 patients. (C) ELISA-determined GDF15 level as a function of disease severity. (D) ELISA-determined GDF15 levels compared to antibody array-determined levels showed a strong correlation (R = 0.87, P < 2.2e−16).
To validate the results of the antibody array profile, we used ELISA to measure the levels of IL-10, IP-10, GDF15, eotaxin-2 and MCP-2. We found the results of the antibody array to be well correlated with the ELISA results, indicating the data were reliable. The GDF15 expression as determined by ELISA showed a strong correlation with the expression as determined by antibody array (R = 0.87, P < 2.2e−16) (Figures 4C, D).
GO enrichment analysis included biological process (BP), cell component (CC), and molecular function (MF) terms. As displayed in Supplementary Figure 2A, CCs including GO: 0005576 (extracellular region) and GO: 0005615 (extracellular space) were significantly associated with differentiation. KEGG pathways with P < 0.05 are shown in Supplementary Figure 2B. Including hsa04061 (Viral protein interaction with cytokine and cytokine receptor) and hsa04060 (Cytokine-cytokine receptor interaction) were significantly associated with the identified differentially expressed factors.
Discussion
As the number of COVID-19 patients continues increasing dramatically worldwide, the treatment of severe cases has become a major challenge. Consequently, the early recognition of severe forms of COVID-19 is essential for achieving a timely triaging of patients.
In a cohort of 111 patients who tested positive for COVID-19, studied here, we found that both the average age and the proportion of males was increased in the severe group compared with the non-severe group, which is consistent with the general recognition of age being an important risk factor for developing severe COVID-19 (25). In our cohort, 65.77% of the patients had at least one comorbidity. High blood pressure was the most common comorbidity in the severe group, followed by diabetes, which is consistent with previous reports suggesting that SARS-CoV-2 is more likely to infect people with chronic comorbidities due to impaired immune function (26, 27). Most interestingly, other comorbidities including liver cysts, anemia, cardiovascular disease and high blood lipid levels were observed in our study, which, to the best of our knowledge, have not been previously reported.
Studies have reported a significant increase in white blood cell (WBC) count in patients with severe COVID-19 that is associated with poor prognosis. This phenomenon is consistent with the findings in our study. The increase in WBC counts in severely affected patients may have been driven by elevated neutrophil levels, as well as decreased levels of lymphocytes and eosinophils (28, 29). Several studies have reported that a decrease in lymphocytes and eosinophils was the most commonly found clinical characteristics of patients with COVID-19 (30). In the current study, the lymphocyte and eosinophil counts decreased gradually with increasing disease severity, but then increased at the convalescent stage. The observed increased NLR (i.e. the combination of a decrease in lymphocyte levels and an increase in neutrophil levels) is a hallmark of inflammation and infection throughout the body and has been studied as a predictor of bacterial infection (31). Our study found that the NLR was increased in COVID-19 patients, most significantly in the severe group, compared with healthy controls. This was consistent with the findings of other studies (32, 33), indicating a serious interference of the internal environment of the body by the potential critical conditions of these severe infections. CRP, a routinely measured inflammatory marker, that was also observed in the current study, has previously been shown to be elevated in most patients with COVID-19 (34). In our study, CRP was significantly elevated in the severe group, and appeared to correlate with disease severity.
The cytokines IL-6 and IL-8 have been previously suggested to be associated with the severity of COVID-19 (15). Here, we found the serum levels of IL-6 and IL-8 to have increased with increasing COVID-19 severity with COVID-19 disease progression. Moreover, we identified eleven factors: TNF-related apoptosis-inducing ligand (TRAIL)-R1, IGFBP-1, IGFBP-4, VCAM-1, sFRP-3, FABP2, transferrin, GDF15, IL-1F7, IL-5Rα, and CD200 differentially expressed among asymptomatic, mild, moderate, severe, and convalescent COVID-19 patients compared with healthy controls using an antibody array analysis of 440 factors.
TRAIL at high serum concentrations has been proposed to be an important mediator of adipose tissue inflammation and obesity-associated diseases (35, 36). Here, we found much higher levels of TRAIL R1 in the 8 patients with diabetes compared to non-diabetic patients, indicating that TRAIL R1 might contribute to diabetes-related COVID-19 disease severity. Insulin-like growth factor 1 (IGF-1) has been reported to play critical roles in diverse biological activities including cell metabolism, tissue regeneration, and survival (37). IGF-binding proteins (IGFBPs) bind both insulin-like growth factors (IGFs) I and II, which can circulate in the plasma in both glycosylated and non-glycosylated forms. Binding of any of these proteins to an IGF prolongs the half-life of the IGF and alters the interaction of the IGF with cell surface receptors (38). In this study, IGFBP-1 and IGFBP-4 had the highest expression levels in the severe group, suggesting a role in organ damage and consequent activation of the wound-healing process (39). In addition, the levels of IL-5Rα, FABP2, CD200 and transferrin consistently showed the highest levels in the moderate group, indicating that these factors may serve as sensitive biomarkers reflecting COVID-19 disease severity.
Growth differentiation factor-15 (GDF15) was previously reported to be associated with inflammation, and many biological processes and diseases, including cancer and obesity (40). The GDF15 gene is located on human chromosome 19p13.1-13.2. GDF15 is synthesized as pro-GDF15, which is then dimerized, cleaved and secreted into the circulation as a mature GDF15 dimer (41). Both the intracellular GDF15 and the circulating mature GDF15 are implicated in biological processes such as energy homeostasis and body weight regulation. GDF15 is a member of the transforming growth factor β (TGF-β) superfamily and is widely expressed in a variety of mammalian tissues, where it is highly regulated and often inducibly expressed under conditions associated with cellular stress (42). GDF15 is also well known to be a stress-induced factor released in response to tissue damage (43). Zimmers et al. reported that GDF15 expression is upregulated in mouse models of kidney and lung injury (44). Hsiao et al. have also shown a significant and rapid increase in GDF15 levels following injury to liver hepatocytes or bile duct epithelial cells (45). Serum levels of GDF15 are strongly associated with many diseases, and functional roles of GDF15 have been reported in cardiovascular, renal, cancer, and metabolic diseases (40). However, knowledge of its pathophysiological function at the molecular level is still limited and more studies are needed. The recent discovery that GFRAL is the endogenous receptor for GDF15 may provide additional insight into the molecular mechanism of GDF15 and its relationship to disease states (46). A recent experimental study has shown a novel function of GDF15 in promoting rhinovirus-associated lung inflammation in humans, which can lead to severe respiratory viral infections (47). A relationship between GDF15 and COVID-19 was reported recently. Myhre et al. showed that GDF15 concentrations were not only independently associated with risk but also used to identify key pathophysiological features of COVID-19, including hypoxemia in comparison of non-ICU and ICU (or death) patients. They also provided prognostic information superior to that provided by the established risk markers in COVID-19 (48). In another study, Notz et al. similarly reported an increase in the level of GDF15 during ICU treatment, suggesting increased tissue resilience (49). Remarkably, in the current study, we found GDF15 levels to be closely associated with disease progression of COVID-19.
To further evaluate the potential of GDF15 as a biomarker of disease severity, we longitudinally collected serum samples, specifically from three groups with different disease severities (moderate, severe, and critical) at three time points (admission, remission, and discharge). The average level of GDF15 increased with severity. Specifically, GDF15 levels were 305.930 ± 85.8pg/mL in the moderate group, 422.5 ± 80.9pg/mL in the severe group, and 621pg/mL in the critical group, with the critical group having a nearly 1.47-times increase compared with the severe group, and a 1.88-times increase compared to the moderate group. In the four patients with severe symptoms, the GDF15 levels showed a linear decrease from admission to the hospital to remission of symptoms to discharge from the hospital. In the ten patients with moderate symptoms, GDF15 increased slightly from admission to remission, but dropped dramatically from remission to discharge. Coincidentally, one critical patient had very high levels of GDF15 at admission, which decreased slightly at remission. Upon worsening of symptoms prior to death, we found that the GDF15 level increased. These observations provided evidence that GDF15 might serve as an indicator of disease severity in COVID-19 patients.
Our findings provide important knowledge related to the severity stratification of COVID-19. However, due to the improvement of the epidemic situation in China, the sample size of this study was limited, especially for the critically ill patients. Our findings should be validated in a larger, more diverse population and environment.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics Statement
The studies involving human participants were reviewed and approved by the ethics committees of Shunde Hospital of Guangzhou University of Chinese Medicine. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JHX, DWW, YLS, YTL, RPH designed the study; YLS, YXL provided samples; JYH, JQZ, YLH, DDL, YZL conducted experiments; YQY, SHL, SZZ analyzed data; XYT, JQZ participated in the manuscript writing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by the Science and Technology Innovation Project of Foshan Municipality (2020001000431), the National Key Research and Development Project (2020YFA0708001), the National Science and Technology Major Project of China (No. 2018ZX10732-401-003-012) and the Industry-University-Research Collaborative Innovation Special Project of Guangzhou (201802030001).
Conflict of Interest
Authors YQY, SHL, SZZ and RPH were employed by company RayBiotech Inc.
The remaining 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.
Acknowledgments
We thank Lin Li, Wei Hu, Wei Wen from Shunde Hospital of Guangzhou University of Chinese Medicine, Xing Chen, Mingkai Tan from Guangzhou Eighth People’s Hospital. We also thank Yunru Che, Wenli Yan from RayBiotech, Inc. for helping with the cytokine array experiment.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.662465/full#supplementary-material
References
1. Zhou P, Yang X-L, Wang X-G, Hu B, Zhang L, Zhang W, et al. A Pneumonia Outbreak Associated With a New Coronavirus of Probable Bat Origin. Nature (2020) 579(7798):270–3. doi: 10.1038/s41586-020-2012-7
2. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical Features of Patients Infected With 2019 Novel Coronavirus in Wuhan, China. Lancet (London England) (2020) 395(10223):497–506. doi: 10.1016/S0140-6736(20)30183-5
3. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A Novel Coronavirus From Patients With Pneumonia in China, 2019. New Engl J Med (2020) 382(8):727–33. doi: 10.1056/NEJMoa2001017
4. Ying T, Du L, Ju TW, Prabakaran P, Lau CC, Lu L, et al. Exceptionally Potent Neutralization of Middle East Respiratory Syndrome Coronavirus by Human Monoclonal Antibodies. J Virol (2014) 88(14):7796–805. doi: 10.1128/JVI.00912-14
5. Tang D, Comish P, Kang R. The Hallmarks of COVID-19 Disease. PloS Pathog (2020) 16(5):e1008536. doi: 10.1371/journal.ppat.1008536
6. Wang D, Hu B, Hu C, Zhu F, Liu X, Zhang J, et al. Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan, China. JAMA (2020) 323(11):1061–69. doi: 10.1001/jama.2020.1585
7. Wang W, Xu Y, Gao R, Lu R, Han K, Wu G, et al. Detection of SARS-CoV-2 in Different Types of Clinical Specimens. JAMA (2020) 323(18):1843–4. doi: 10.1001/jama.2020.3786
8. Catanzaro M, Fagiani F, Racchi M, Corsini E, Govoni S, Lanni C. Immune Response in COVID-19: Addressing a Pharmacological Challenge by Targeting Pathways Triggered by SARS-CoV-2. Signal Transduction Targeted Ther (2020) 5(1):84. doi: 10.1038/s41392-020-0191-1
9. Picchianti Diamanti A, Rosado MM, Pioli C, Sesti G, Laganà B. Cytokine Release Syndrome in COVID-19 Patients, A New Scenario for an Old Concern: The Fragile Balance Between Infections and Autoimmunity. Int J Mol Sci (2020) 21(9):3330. doi: 10.3390/ijms21093330
10. Ye Q, Wang B, Mao J. The Pathogenesis and Treatment of the `Cytokine Storm’ in COVID-19. J Infection (2020) 80(6):607–13. doi: 10.1016/j.jinf.2020.03.037
11. Song P, Li W, Xie J, Hou Y, You C. Cytokine Storm Induced by SARS-CoV-2. Clinica Chimica Acta; Int J Clin Chem (2020) 509:280–7. doi: 10.1016/j.cca.2020.06.017
12. Tay MZ, Poh CM, Rénia L, MacAry PA, Ng LFP. The Trinity of COVID-19: Immunity, Inflammation and Intervention. Nat Rev Immunol (2020) 20(6):363–74. doi: 10.1038/s41577-020-0311-8
13. Giamarellos-Bourboulis EJ, Netea MG, Rovina N, Akinosoglou K, Antoniadou A, Antonakos N, et al. Complex Immune Dysregulation in COVID-19 Patients With Severe Respiratory Failure. Cell Host Microbe (2020) 27(6):992–1000.e3. doi: 10.1016/j.chom.2020.04.009
14. Jamilloux Y, Henry T, Belot A, Viel S, Sève P. Should We Stimulate or Suppress Immune Responses in COVID-19? Cytokine and Anti-Cytokine Interventions. Autoimmun Rev (2020) 19(7):102567. doi: 10.1016/j.autrev.2020.102567
15. Del Valle DM, Kim-Schulze S, Huang H-H, Beckmann ND, Nirenberg S, Wang B, et al. An Inflammatory Cytokine Signature Predicts COVID-19 Severity and Survival. Nat Med (2020) 26(10):1636–43. doi: 10.1038/s41591-020-1051-9
16. Yang Y, Shen C, Li J, Yuan J, Wei J, Huang F, et al. Plasma IP-10 and MCP-3 Levels Are Highly Associated With Disease Severity and Predict the Progression of COVID-19. J Allergy Clin Immunol (2020) 146(1):119–27. doi: 10.1016/j.jaci.2020.04.027
17. Ueland T, Heggelund L, Lind A, Holten AR, Tonby K, Michelsen AE, et al. Elevated Plasma sTIM-3 Levels in Patients With Severe COVID-19. J Allergy Clin Immunol (2021) 147(1):92–8. doi: 10.1016/j.jaci.2020.09.007
18. Tamayo-Velasco Á, Martínez-Paz P, Peñarrubia-Ponce MJ, de la Fuente I, Pérez-González S, Fernández I, et al. HGF, IL-1α, and IL-27 Are Robust Biomarkers in Early Severity Stratification of COVID-19 Patients. J Clin Med (2021) 10(9):2017. doi: 10.3390/jcm10092017
19. Tang J, Lin J, Zhang E, Zhong M, Luo Y, Fu Y, et al. Serum IL-6 and Procalcitonin Are Two Promising Novel Biomarkers for Evaluating the Severity of COVID-19 Patients. Medicine (2021) 100(22):e26131. doi: 10.1097/MD.0000000000026131
20. Huang R, Zhang L, Gu Q, Zhou Y-H, Hao Y, Zhang K, et al. Profiles of Acute Cytokine and Antibody Responses in Patients Infected With Avian Influenza A H7n9. PloS One (2014) 9(7):e101788. doi: 10.1371/journal.pone.0101788
21. Cheung CY, Poon LLM, Lau AS, Luk W, Lau YL, Shortridge KF, et al. Induction of Proinflammatory Cytokines in Human Macrophages by Influenza A (H5N1) Viruses: A Mechanism for the Unusual Severity of Human Disease? Lancet (London England) (2002) 360(9348):1831–7. doi: 10.1016/S0140-6736(02)11772-7
22. Wiedemann A, Foucat E, Hocini H, Lefebvre C, Hejblum BP, Durand M, et al. Long-Lasting Severe Immune Dysfunction in Ebola Virus Disease Survivors. Nat Commun (2020) 11(1):3730. doi: 10.1038/s41467-020-17489-7
23. National Health Commission of the People’s Republic of China. Protocol on Prevention and Control of Novel Coronavirus Pneumonia, Edition 6. (2020).
24. National Health Commission of the People’s Republic of China. Diagnosis and Treatment Protocol for Novel Coronavirus Pneumonia (Trial Version 7). (2020).
25. Romero Starke K, Petereit-Haack G, Schubert M, Kämpf D, Schliebner A, Hegewald J, et al. The Age-Related Risk of Severe Outcomes Due to COVID-19 Infection: A Rapid Review, Meta-Analysis, and Meta-Regression. Int J Environ Res Public Health (2020) 17(16):5974. doi: 10.3390/ijerph17165974
26. Erener S. Diabetes, Infection Risk and COVID-19. Mol Metab (2020) 39:101044. doi: 10.1016/j.molmet.2020.101044
27. Guan W-j, Liang W-h, Zhao Y, Liang H-r, Chen Z-s, Li Y-m, et al. Comorbidity and its Impact on 1590 Patients With COVID-19 in China: A Nationwide Analysis. Eur Respir J (2020) 55(5):2000547. doi: 10.1183/13993003.00547-2020
28. Zhao Y, Nie H-X, Hu K, Wu X-J, Zhang Y-T, Wang M-M, et al. Abnormal Immunity of Non-Survivors With COVID-19: Predictors for Mortality. Infect Dis Poverty (2020) 9(1):108. doi: 10.1186/s40249-020-00723-1
29. Zhao K, Li R, Wu X, Zhao Y, Wang T, Zheng Z, et al. Clinical Features in 52 Patients With COVID-19 Who Have Increased Leukocyte Count: A Retrospective Analysis. Eur J Clin Microbiol Infect Dis Off Publ Eur Soc Clin Microbiol (2020) 39(12):2279–87. doi: 10.1007/s10096-020-03976-8
30. Sun S, Cai X, Wang H, He G, Lin Y, Lu B, et al. Abnormalities of Peripheral Blood System in Patients With COVID-19 in Wenzhou, China. Clinica Chimica Acta; Int J Clin Chem (2020) 507:174–80. doi: 10.1016/j.cca.2020.04.024
31. He L, Wang J, Wang F, Zhang L, Zhang L, Zhao W. Increased Neutrophil-to-Lymphocyte Ratio Predicts the Development of Post-Stroke Infections in Patients With Acute Ischemic Stroke. BMC Neurol (2020) 20(1):328. doi: 10.1186/s12883-020-01914-x
32. Zhang B, Zhou X, Zhu C, Song Y, Feng F, Qiu Y, et al. Immune Phenotyping Based on the Neutrophil-To-Lymphocyte Ratio and IgG Level Predicts Disease Severity and Outcome for Patients With COVID-19. Front Mol Biosci (2020) 7:157. doi: 10.3389/fmolb.2020.00157
33. Li X, Liu C, Mao Z, Xiao M, Wang L, Qi S, et al. Predictive Values of Neutrophil-to-Lymphocyte Ratio on Disease Severity and Mortality in COVID-19 Patients: A Systematic Review and Meta-Analysis. Crit Care (London England) (2020) 24(1):647. doi: 10.1186/s13054-020-03374-8
34. Zhang Y, Zheng L, Liu L, Zhao M, Xiao J, Zhao Q. Liver Impairment in COVID-19 Patients: A Retrospective Analysis of 115 Cases From a Single Centre in Wuhan City, China. Liver Int Off J Int Assoc Study Liver (2020) 40(9):2095–103. doi: 10.1111/liv.14455
35. Keuper M, Wernstedt Asterholm I, Scherer PE, Westhoff MA, Möller P, Debatin KM, et al. TRAIL (TNF-Related Apoptosis-Inducing Ligand) Regulates Adipocyte Metabolism by Caspase-Mediated Cleavage of PPARgamma. Cell Death Dis (2013) 4:e474. doi: 10.1038/cddis.2012.212
36. Ashley DT, O’Sullivan EP, Davenport C, Devlin N, Crowley RK, McCaffrey N, et al. Similar to Adiponectin, Serum Levels of Osteoprotegerin Are Associated With Obesity in Healthy Subjects. Metabolism: Clin Exp (2011) 60(7):994–1000. doi: 10.1016/j.metabol.2010.10.001
37. Vitale G, Pellegrino G, Vollery M, Hofland LJ. ROLE of IGF-1 System in the Modulation of Longevity: Controversies and New Insights From a Centenarians’ Perspective. Front Endocrinol (2019) 10:27. doi: 10.3389/fendo.2019.00027
38. Durai R, Davies M, Yang W, Yang S, Seifalian A, Goldspink G, et al. Biology of Insulin-Like Growth Factor Binding Protein-4 and its Role in Cancer (Review). Int J Oncol (2006) 28(6):1317–25. doi: 10.3892/ijo.28.6.1317
39. Fan X, Yin C, Wang J, Yang M, Ma H, Jin G, et al. Pre-Diagnostic Circulating Concentrations of Insulin-Like Growth Factor-1 and Risk 3 of COVID-19 Mortality: Results From UK Biobank. Eur J Epidemiol (2021) 36(3):311–8. doi: 10.1007/s10654-020-00709-1
40. Baek SJ, Eling T. Growth Differentiation Factor 15 (GDF15): A Survival Protein With Therapeutic Potential in Metabolic Diseases. Pharmacol Ther (2019) 198:46–58. doi: 10.1016/j.pharmthera.2019.02.008
41. Bauskin AR, Jiang L, Luo XW, Wu L, Brown DA, Breit SN. The TGF-Beta Superfamily Cytokine MIC-1/GDF15: Secretory Mechanisms Facilitate Creation of Latent Stromal Stores. J Interferon Cytokine Res Off J Int Soc Interferon Cytokine Res (2010) 30(6):389–97. doi: 10.1089/jir.2009.0052
42. Wischhusen J, Melero I, Fridman WH. Growth/Differentiation Factor-15 (GDF-15): From Biomarker to Novel Targetable Immune Checkpoint. Front Immunol (2020) 11:951. doi: 10.3389/fimmu.2020.00951
43. Emmerson PJ, Duffin KL, Chintharlapalli S, Wu X. GDF15 and Growth Control. Front Physiol (2018) 9:1712. doi: 10.3389/fphys.2018.01712
44. Zimmers TA, Jin X, Hsiao EC, McGrath SA, Esquela AF, Koniaris LG. Growth Differentiation Factor-15/Macrophage Inhibitory Cytokine-1 Induction After Kidney and Lung Injury. Shock (Augusta Ga.) (2005) 23(6):543–8. doi: 10.1097/01.shk.0000163393.55350.70
45. Hsiao EC, Koniaris LG, Zimmers-Koniaris T, Sebald SM, Huynh TV, Lee SJ. Characterization of Growth-Differentiation Factor 15, a Transforming Growth Factor Beta Superfamily Member Induced Following Liver Injury. Mol Cell Biol (2000) 20(10):3742–51. doi: 10.1128/.20.10.3742-3751.2000
46. Rochette L, Zeller M, Cottin Y, Vergely C. Insights Into Mechanisms of GDF15 and Receptor GFRAL: Therapeutic Targets. Trends Endocrinol Metabolism: TEM (2020) 31(12):939–51. doi: 10.1016/j.tem.2020.10.004
47. Wu Q, Jiang D, Schaefer NR, Harmacek L, O’Connor BP, Eling TE, et al. Overproduction of Growth Differentiation Factor 15 Promotes Human Rhinovirus Infection and Virus-Induced Inflammation in the Lung. Am J Physiol Lung Cell Mol Physiol (2018) 314(3):L514–27. doi: 10.1152/ajplung.00324.2017
48. Myhre PL, Prebensen C, Strand H, Røysland R, Jonassen CM, Rangberg A, et al. Growth Differentiation Factor 15 Provides Prognostic Information Superior to Established Cardiovascular and Inflammatory Biomarkers in Unselected Patients Hospitalized With COVID-19. Circulation (2020) 142(22):2128–37. doi: 10.1161/CIRCULATIONAHA.120.050360
Keywords: COVID-19, inflammatory factors, GDF15, biomarker, disease severity
Citation: Teng X, Zhang J, Shi Y, Liu Y, Yang Y, He J, Luo S, Huang Y, Liu Y, Liu D, Li Y, Zhang S, Huang R-P, Wang D and Xu J (2021) Comprehensive Profiling of Inflammatory Factors Revealed That Growth Differentiation Factor-15 Is an Indicator of Disease Severity in COVID-19 Patients. Front. Immunol. 12:662465. doi: 10.3389/fimmu.2021.662465
Received: 01 February 2021; Accepted: 24 June 2021;
Published: 15 July 2021.
Edited by:
Jörg Wischhusen, University Hospital Würzburg, GermanyReviewed by:
Florian Wedekink, University Hospital Würzburg, GermanyMichel Guy André Mittelbronn, National Health Laboratory, Luxembourg
Copyright © 2021 Teng, Zhang, Shi, Liu, Yang, He, Luo, Huang, Liu, Liu, Li, Zhang, Huang, Wang and Xu. 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.
*Correspondence: Ruo-Pan Huang, rhuang@raybiotech.com; Dawei Wang, david@gzucm.edu.cn; Jianhua Xu, jhxu1976@gzucm.edu.cn
†These authors have contributed equally to this work