key: cord-0916912-q0kmq9wx
authors: Liu, Jiao; Shen, Yanfei; Wen, Zhenliang; Xu, Qianghong; Wu, Zhixiong; Feng, Huibin; Li, Zhongyi; Dong, Xuan; Huang, Sisi; Guo, Jun; Zhang, Lidi; Chen, Yizhu; Li, Wenzhe; Zhu, Wei; Du, Hangxiang; Liu, Yongan; Wang, Tao; Chen, Limin; Teboul, Jean-Louis; Annane, Djillali; Chen, Dechang
title: Efficacy of Thymosin Alpha 1 in the Treatment of COVID-19: A Multicenter Cohort Study
date: 2021-08-02
journal: Front Immunol
DOI: 10.3389/fimmu.2021.673693
sha: f1fbdcf0a0a02ac6287784fb033df5ce1358bb8d
doc_id: 916912
cord_uid: q0kmq9wx

BACKGROUND: Thymosin alpha 1 (Tα1) is widely used to treat patients with COVID-19 in China; however, its efficacy remains unclear. This study aimed to explore the efficacy of Tα1 as a COVID-19 therapy. METHODS: We performed a multicenter cohort study in five tertiary hospitals in the Hubei province of China between December 2019 and March 2020. The patient non-recovery rate was used as the primary outcome. RESULTS: All crude outcomes, including non-recovery rate (65/306 vs. 290/1,976, p = 0.003), in-hospital mortality rate (62/306 vs. 271/1,976, p = 0.003), intubation rate (31/306 vs. 106/1,976, p = 0.001), acute respiratory distress syndrome (ARDS) incidence (104/306 vs. 499/1,976, p = 0.001), acute kidney injury (AKI) incidence (26/306 vs. 66/1,976, p < 0.001), and length of intensive care unit (ICU) stay (14.9 ± 12.7 vs. 8.7 ± 8.2 days, p < 0.001), were significantly higher in the Tα1 treatment group. After adjusting for confounding factors, Tα1 use was found to be significantly associated with a higher non-recovery rate than non-Tα1 use (OR 1.5, 95% CI 1.1–2.1, p = 0.028). An increased risk of non-recovery rate associated with Tα1 use was observed in the patient subgroups with maximum sequential organ failure assessment (SOFA) scores ≥2 (OR 2.0, 95%CI 1.4–2.9, p = 0.024), a record of ICU admission (OR 5.4, 95%CI 2.1–14.0, p < 0.001), and lower PaO2/FiO2 values (OR 1.9, 95%CI 1.1–3.4, p = 0.046). Furthermore, later initiation of Tα1 use was associated with a higher non-recovery rate. CONCLUSION: Tα1 use in COVID-19 patients was associated with an increased non-recovery rate, especially in those with greater disease severity.

Coronavirus 2019 disease (COVID- 19) represents an ongoing global threat to human health and has caused more than 2,400,000 deaths to date (JHU data-23/02/2021, https:// coronavirus.jhu.edu/). As the COVID-19 outbreak continues, novel and effective therapies are urgently needed.

Current evidence indicates that immune dysfunction in COVID-19 plays a critical role in disease severity and is associated with a poor prognosis. Several studies have reported that inflammatory markers such as C-reactive protein (CRP), ferritin, interleukin (IL)-6, IL-10, and tumor necrosis factor-a (TNF-a) were significantly elevated in severe COVID-19 cases. The resulting cytokine storm is thought to cause rapid COVID-19 progression and increased deaths (1) (2) (3) . Furthermore, Yu and colleagues (4) concluded that adaptive immune response dysregulation contributes to the cytokine storm and further leads to severe COVID-19. Dozens of studies have been performed to explore potential medical interventions to attenuate the harmful immune response, including statins (5), corticosteroids (6, 7) , and thymosin alpha 1 (Ta1) (8) (9) (10) (11) . However, the benefits of immunomodulatory treatment of patients with COVID-19 have remained controversial.

Ta1 is an endogenous polypeptide hormone secreted by thymic epithelial cells (12) . As an immunomodulatory therapy (12, 13) , Ta1 has been investigated in many diseases involving immune dysfunction [such as sepsis, cystic fibrosis, and hepatitis viral infection (14) (15) (16) ] and is associated with an improved patient outcome. The efficacy of Ta1 in COVID-19 has also been investigated in several studies (8) (9) (10) (11) . However, the conclusions remain inconsistent. For instance, in a retrospective study including 76 severe COVID-19 cases (8), Liu et al. reported that Ta1 use could increase T cell numbers and was associated with reduced mortality. However, other studies found that Ta1 use showed no effect on mortality in COVID-19 (9) or was even associated with increased mortality in patients with COVID-19 of all severity levels (10) . The reasons behind these inconsistent findings remain unclear. However, emerging evidence indicated that COVID-19 is a heterogeneous disease (17, 18) . Whether the role of Ta1 in COVID-19 is affected by different phenotypes remains unclear.

To this end, we performed a large multicenter cohort study to explore the efficacy of Ta1 in COVID-19 patients displaying a range of disease severities. Interaction analysis between Ta1 and predefined parameters was also performed to identify the heterogeneous effectiveness of Ta1 in different patient subgroups.

This was a multicenter cohort study conducted at five tertiary hospitals (in the Hubei province of China) that were designated 

Demographic characteristics, smoking and comorbidities data, and initial laboratory indices recorded within the first 24 h after hospital admission were extracted from each center through electronic medical records. The change in lymphocyte counts was defined as the difference between the initial lymphocyte count and the minimum lymphocyte count. Clinical classification, special drug combinations, and clinical outcomes data were also recorded. The sequential organ failure assessment (SOFA) scores (20) and acute physiology and chronic health evaluation II (APACHE II) scores (21) were calculated to assess the severity of illness using data obtained in the first 24 h after hospital admission. Ta1 use was defined as any record of subcutaneous/intramuscular Ta1 administration throughout the entire hospital stay. Both the duration and timing of Ta1 use were also extracted.

The detailed diagnosis of COVID-19 (19) was established as follows: 1. patients with an epidemiological history and chest imaging (computed tomography or radiography) that suggested viral pneumonia; 2. laboratory-confirmed SARS-COV-2 infection of throat swab, sputum, and/or lower respiratory tract samples by high-throughput sequencing or real-time reverse transcriptase polymerase chain reaction (RT-PCR); or 3. confirmed plasma positivity for specific antibody (IgM or/and IgG) against SARS-COV-2.

COVID-19 disease severity was classified according to the Chinese diagnosis and COVID-19 treatment guidelines (Trial Version 7) (19) . Briefly, patients were diagnosed with moderately severe COVID-19 if they had fever, respiratory symptoms, and chest imaging (computed tomography or radiography) suggesting viral pneumonia. Patients who met either conditions i-iii or iv-vi of the following criteria were diagnosed as having severe or critical disease, respectively: (i) respiratory rate ≥30 beats/min; (ii) resting stable oxygen saturation (SpO2) ≤93%; (iii) PaO2/FiO2 ≤300 mmHg (1 mmHg = 0.133 kPa); (iv) respiratory failure requiring mechanical ventilation; (v) shock; or (vi) multiple organ failure requiring ICU admission.

Diagnosis of acute respiratory distress syndrome (ARDS) was made according to the Berlin definition (22) . Acute kidney injury (AKI) was defined according to the "Kidney Disease: Improving Global Outcomes" document (23) .

Subgroup analysis was performed to interpret possible interactions between Ta1 use and disease severity. Data were stratified according to the median of the maximum SOFA score (≥/<2), PaO2/FiO2 at admission, admission to intensive care unit (ICU), application of mechanical ventilation, and development of ARDS. Non-recovery rate was the primary outcome at the time of data collection. Non-recovering patients were those who died as a result of COVID-19 or who were still in hospital, but in any one of severely deteriorated condition: a. after sufficient fluid resuscitation (≥30 ml/kg crystalloid solution), a dose of vasopressor therapy equal to or greater than an equivalent of norepinephrine 0.4 mg/kg/min for at least 24 h was hard to maintain mean arterial pressure (MAP) at 65 mmHg; b. PaO2/FiO2 was continuously lower than 100 mmHg for 24 h even with invasive mechanical ventilation. Secondary outcomes included in-hospital mortality, ICU admission, intubation rate, ARDS incidence, AKI incidence, duration of mechanical ventilation, length of ICU stay, and length of hospital stay. Sensitivity analysis was also performed using in-hospital mortality as the primary outcome.

Most continuous variables in the study had less than 5% missing data and were therefore replaced with their mean or median values. Variables with more than 20% missing data were not completed. For dichotomous variables, the missing value was replaced with the default value (zero).

Continuous variables were presented as mean ± standard deviation (SD) or median (interquartile range, IQR) according to the data distribution. The Student's t-test or Wilcoxon ranksum test were used as appropriate. Categorical data were compared using the Chi-squared test and were presented as percentages. To adjust for potential confounding factors, variables with p < 0.2 in the univariate comparison were included in the initial model. Stepwise regression was used to build the final logistic model. Multicollinearity was tested using the variance inflation factor (VIF) method, with a VIF ≥ 5 suggesting potential multicollinearity. To test the stability of these logistic models, bootstrap analysis (using 100 resamples) was performed.

The significant impact of maximum SOFA scores on the association between Ta1 use and the non-recovery rate was detected in the multivariable analysis. Thus, interactions between Ta1 use and disease severity indices were explored, and subgroup analysis was performed according to the following five parameters: median SOFA score (≥/<2), ICU admission, PaO2/FiO2 value, application of mechanical ventilation, and development of ARDS. For interpretation, the predicted marginal effect of Ta1 was estimated at different SOFA scores and PaO2/FiO2 values.

Propensity score matching (PSM) was applied to minimize the impact of confounding factors, such as biochemical indices and the disease severity score, which may lead to a biased outcome. The propensity score was assigned based on the probability that a patient would receive Ta1 therapy and estimated using a multivariable logistic regression model. A one-to-three nearest neighbor matching algorithm was applied with a caliper of 0.02. The following variables were selected to generate the propensity score: age, hypertension, chronic obstructive pulmonary disease (COPD), chronic liver disease, chronic cardiac disease, maximum SOFA score, ICU admission, serum hemoglobin and creatinine levels, and initial white blood cell, lymphocyte, and platelet counts. The bias in the means (or proportions) of covariates between two groups was examined using the standardized difference before and after PSM. Finally, to test bias owing to imbalances in unmeasured covariates in the PSM, sensitivity analyses were performed to quantify the degree of hidden bias (Gamma value) that would need to be present to invalidate our main conclusions. A two-tailed test was performed, and p < 0.05 was considered statistically significant. All statistical analyses were performed using STATA 14.0 (College Station, TX, USA).

Of the 2,411 patients admitted to the five participating tertiary hospitals, 73 were excluded because of unconfirmed COVID-19, 51 were excluded because of duplicate records, and 5 were excluded due to missing outcome data. Of the remaining 2,282 patients that were eligible to participate in our study, 306 received Ta1 therapy ( Figure 1 ). Clinical characteristics associated with COVID-19 patients in the Ta1 and non-Ta1 treatment groups are shown in Table 1 . There were no significant differences in terms of age or gender between the two groups. More patients received corticosteroids (107/306 vs. 382/1,976, p<0.001) and interferon (74/306 vs. 108/1,976, p<0.001) therapy in the Ta1 group than in the non-Ta1 group. Patients who received Ta1 therapy had almost the same disease severity according to the SOFA score, APACHE II score, initial serum creatinine and maximum lactate levels, although the patients in the Ta1 group included a slightly larger proportion of critically ill patients (59/306 vs. 240/1,976). Furthermore, although both the initial and minimal lymphocyte counts were significantly lower in the Ta1 group than in the non-Ta1 group (1.0 ± 0.5 vs. 1.1 ± 0.6 × 10^9/L, p = 0.007; 0.8 ± 0.5 vs. 0.9 ± 0.6 × 10^9/L, p < 0.001 for the Ta1 and non-Ta1 groups, respectively), the change in lymphocyte count was comparable between the groups (0.1 ± 0.3 vs. 0.1 ± 0.3 × 10^9/L, p = 0.337). However, as reported in Table 2 vs. 8.7 ± 8.2 days, p < 0.001), were significantly higher in the Ta1 group than in the non-Ta1 group. Information relating to the number of patients, Ta1 use, and the non-recovery rate for each center is provided in Supplementary Table S1 .

Potential confounding factors were adjusted for in the multivariable logistic models as shown in the Table 3 . The odds ratio (OR) for the non-recovery rate of Ta1 use was found to be significant (Model 2: OR 1.5, 95% confidence interval (CI) 1.1-2.1, p = 0.028) using bootstrap analysis with 1,000 resamples. However, after the addition of the maximum SOFA score to Model 3, the association between Ta1 use and the non-recovery rate became non-significant (OR 1.4, 95% CI 0.9-2.2, p = 0.118).

Interactions between Ta1 use and different disease severity indices were evaluated. There were significant interactions between Ta1 use and the maximum SOFA score (p = 0.024) and ICU admission (p < 0.001) ( Figure 2 ). Ta1 use was significantly associated with an increase in non-recovery rate in the subgroup with a maximum SOFA score ≥2 (OR 2.0, 95% CI 1.4-2.9) or ICU admission (OR 5.4, 95% CI 2.1-14.0), whereas this association was not significant in subgroups with a SOFA score <2 (OR 0.6, 95% CI 0.24-1.7) or without ICU admission (OR 1.1, 95% CI 0.7-1.7). There was no significant interaction between Ta1 use and PaO2/FiO2 values (≥300 mmHg) at admission, mechanical ventilation, or ARDS development. However, a trend towards higher ORs was observed in subgroups with low PaO2/FiO2 values (<300 mmHg), mechanical ventilation, and ARDS. Furthermore, significant interactions were also observed using maximum SOFA (p = 0.003) and PaO2/FiO2 values at admission (p = 0.016) as continuous variables. The predicted marginal effect of Ta1 on the non-recovery rate was estimated using different SOFA and PaO2/FiO2 values (Supplementary Figure S1) . The difference in the probability of non-recovery between the Ta1 and non-Ta1 groups correlated positively with the SOFA score or negatively with the PaO2/FiO2 value. In addition, aiming to evaluate whether the effect of Ta1 administration was influenced by other concurrent immunotherapies, we tested the potential interaction between Ta1, corticosteroids, and interferon and found no significant interaction. We also included both the corticosteroids and interferon factors in our multivariable model, and the overall results in each subgroup remained unchanged (Supplementary Table S2 ). The proportions of patients receiving Ta1, corticosteroids, and/or interferon are presented in a Venn diagram (Supplementary Figure S2) .

PSM was performed using a 1:3 algorithm; 306 cases from the Ta1 group and 720 cases from the non-Ta1 group were well matched (Supplementary Tables S3, 4) . The overall quality of the matched sample was assessed by comparing the standardized difference of included covariates before and after PSM. There was no significant difference between the matched groups relating to PaO2/FiO2, ratio of partial pressure of arterial oxygen to fraction of inspired oxygen; SOFA, sequential organ failure assessment; APACHE II, acute physiology and chronic health evaluation II. significantly affected by disease severity, and all the above conclusions were supported ( Table 4 ).

In the non-recovery group (n = 355), there are 18 patients still under treatment, but in a severely deteriorated condition during data extraction. For robustness, we performed sensitivity analysis under three patterns, and "death" was used as the dependent outcome (Supplementary Table S4 ). In pattern 1, these 18 patients were divided into the "alive" group. In pattern 2, these 18 patients were divided into the "death" group. In pattern 3, these 18 patients were excluded from the analysis. The results were stable in all these three patterns. In addition, interactions between thymosin a1 use and disease severity indexes on in-hospital mortality are also presented in Supplementary Table S5 .

Both the duration and timing of Ta1 were included in the multivariable logistic model in different forms ( Three logistic models were used to evaluate the association between non-recovery rate and thymosin a1 use. OR of thymosin a1 use was significantly associated with increased nonrecovery rate in models 1 and 2. However, after adjusting for SOFA score in model 3, the OR of thymosin a1 became non-significant. Bootstrap techniques (100 resamples) was used for calculating 95% CI in models 2 and 3. OR, odds ratio; aORs, adjusted odds ratios; PaO2/FiO2, ratio of partial pressure of arterial oxygen to fraction of inspired oxygen; SOFA, sequential organ failure assessment. median value; and Model 3: tertile analysis). Ta1 use at a later stage was significantly associated with a higher non-recovery rate than Ta1 use at an earlier stage (Models 1-3) . However, the duration of Ta1 use had a less significant effect on recovery rate. Inclusion of the maximum SOFA score as a confounding factor in the logistic model did not alter the results (Supplementary Table S6 ).

The current study has three major findings. First, Ta1 use in COVID-19 was associated with poor clinical outcomes. This finding was confirmed both in the multivariable logistic model and via PSM. Second, there were significant interactions between Ta1 use and disease severity indices, and the association between Ta1 use and the non-recovery rate was stronger in severe COVID-19 cases. Third, Ta1 use at a later stage was significantly associated with a higher non-recovery rate than Ta1 use at an earlier stage. However, the duration of Ta1 use had a less significant effect on recovery rate. Lymphopenia has been commonly reported in patients with severe COVID-19 (4, 24, 25) and is associated with a poor outcome (24) . Yu et al. (4) reported that a significant reduction in effector T cell numbers, accompanied by an increase in the frequency of effector B cells, was observed in COVID-19 patients with severe disease. This indicates that a decline in T cell numbers is likely the main cause of lymphopenia and may aggravate COVID-19 severity. Given the essential role of T cells in viral eradication, exploring medical interventions designed to boost T lymphocyte number and function may improve COVID-19 prognosis (26) . Previous studies have indicated that Ta1 can promote T cell development and proliferation, increase their number, and enhance their function (27, 28) . Yu's research (4) also demonstrated that Ta1 promoted the proliferation of effector T cells in vitro and relieved lymphopenia in COVID-19 patients. Matteucci et al. (29) found that genes associated with cytokine signaling and expression were upregulated in patients with COVID-19, and the ex vivo treatment with Ta1-mitigated cytokine expression, and inhibited lymphocyte activation. These findings provide the basis for the rationale use of Ta1 in COVID-19, but these findings need to be further supported by clinical studies. A recent cohort study showed that Ta1 can reduce mortality in patients with severe COVID-19 by reversing lymphocytopenia and restoring the function of exhausted T cells (8) . However, bias risk should be considered when interpreting this conclusion, due to the limited sample size used (only 76 patients with severe COVID-19 were included) and the presence of potential confounding factors. In another multicenter observation study, Wu et al. (11) also found that Ta1 significantly decreased 28-day mortality and attenuated acute lung injury in critical type COVID-19 patients. However, the definition of critical type is largely different from our study. In addition, the disease severity was significantly imbalanced between Ta1 and non-Ta1 groups, as patients in the Ta1 group had worse condition. Noteworthy, inconsistent findings were also As heterogeneous effects of thymosin a1 was observed in patients with different disease severity levels,PSM was performed both in all patients and in patients with SOFA ≥ 2. The following variables were selected to generate the propensity score: age, ICU admission, COPD, hypertension, chronic cardiac disease, the initial white blood cell, lymphocyte,and platelet count, initial hemoglobin and creatinine level, and maximum SOFA score during hospital admission. Gamma of sensitivity analysis was 1.2. SOFA, sequential organ failure assessment.; ARDS, acute respiratory distress syndrome; AKI, acute kidney injury; CRRT, continuous renal replacement therapy.

reported in other studies. Sun et al. (9) reported that in 771 patients with COVID-19, use of Ta1 was not associated with decreased mortality in critically ill COVID-19 patients. Wang et al. conducted a retrospective propensity score matched study including 317 COVID-19 patients, treatments with immunomodulatory therapies, including glucocorticoids, immunoglobulin, and thymosin, were significantly associated with a higher rate of COVID-19 death (10).

In the current study, Ta1 was used in 13.4% of COVID-19 patients. We noted that Ta1 use was associated with poor prognoses, such as higher non-recovery and in-hospital mortality rate and higher ARDS and AKI incidence. After adjusting for potential confounders, Ta1 use was still significantly associated with an increased non-recovery rate. However, the underlying mechanism remains unclear. One study (30) reported that Ta1 exhibited a dual effect on dendritic cells (DCs) upon different pathogen encounter in vitro. Ta1 was shown to promote DCs to secrete inflammatory cytokines, such as TNF-a, IL-6, and IL-8, in response to viral infection. In contrast, Ta1 dampened inflammation when the DCs were exposed to bacterial pathogens. Furthermore, as a result of the interaction analysis, we found that Ta1 use was associated with increased non-recovery rate mainly in patient subgroups with greater disease severity. Recent studies have found that COVID-19 follows a heterogeneous course of disease (17, 18) . Shi et al. (31) categorized COVID-19 cases into three stages based on inflammatory status: stage I, an asymptomatic incubation period; stage II, non-severe symptomatic period; and stage III, severe respiratory symptomatic stage. For patients in stage III, a hyper-inflammatory response characterized by the overexpression of inflammatory factors and cytokines played a key role in patient deterioration (32) . During this third stage, immunosuppressive therapy as opposed to immune enhancement therapy is needed (31) . However, studies have indicated that Ta1 plays an important role in the activation and maturation of DCs, which in turn increases the secretion of inflammatory factors, promotes the differentiation of T cell precursors into Th1, eventually exacerbating the inflammatory response to viral infection (28, 33) . In our study, the greatest disease severity indices, such as higher SOFA scores, lower PaO2/FiO2 values, or admission to ICU, to a certain degree represented a more severe inflammatory stage. This may be the potential mechanism explaining why adjuvant Ta1 use was associated with poor outcomes especially in severe COVID-19 cases.

In addition, we found that the timing of Ta1 treatment is also crucial as applying Ta1 at a later stage was significantly associated with increased non-recovery rate. Aiming for robustness, the timing of Ta1 was included in the multivariable logistic model in three different forms, and the findings remain stable. We inferred that the inflammatory response may be stronger in the later stage than in the early stage in COVID-19. Therefore, compared to early use, Ta1 use at a later stage may exacerbate the inflammatory response and lead to a poor prognosis. However, data about the inflammatory status at different stages were lacking. More studies are needed to validate our speculation. In addition, we also investigated the impact of Ta1 duration on prognosis in COVID-19. In the tertile analysis, we noted that compared to Tertile 1 (≤3 days), OR in Tertile 2 was significant while non-significant in Tertile 3. When we included the duration of Ta1 in other two forms, the results became nonsignificant. Therefore, there is a chance that this significant finding may be affected by different cut-off values, and the stability of this result should be verified.

To the best of our knowledge, this is the first large, multicenter observational study confirming the association between Ta1 therapy and increased non-recovery rate in severely ill COVID-19 patients. The main strengths of this study included the large sample size and the use of multiple statistical analyses methods to minimize confounding bias. However, several limitations need to be addressed. First, the use of Ta1 for the treatment of COVID-19 patients was mainly at the discretion of clinicians. Thus, the reasons for its use may have been quite varied, such as poor immune status, presence of hepatic disease, severe infection, or simply the clinician's preference. Second, we included most of the covariates, both in a multivariable model and PSM analysis. However, bias caused by potential confounding factors remains possible, such as the time of COVID-19 diagnosis. Third, owing to the retrospective nature of this study, the protocol for Ta1 administration was diverse among patients; for instance, the timing of Ta1 use was significantly associated with a poor outcome, whereas the duration was not. Thus, rigorously designed randomized controlled studies would be required to reach a more precise conclusion. Fourth, the causal relationship between Ta1 and a poor clinical outcome in the context of COVID-19 cannot be inferred from the current study. Although the results from the PSM supported our hypothesis, our work is observational. Therefore, future prospective studies are needed to verify our results.

The hyper-inflammatory response is a hallmark of COVID-19. Ta1 use in these patients is significantly associated with an increased non-recovery rate and a poor clinical outcome, especially in individuals with severe disease. As the COVID-19 outbreak is still ongoing, Ta1 use in this group of patients should be treated with caution.

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 the Ethics Committee of the Jin Yin-tan Hospital (KY-2020-03.01). Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

DC and XD designed the project. YC, LZ, JG, HD, YL,TW, LC, WL, and WZ carried out the data collection. YS analyzed the data and prepared the figures. JL, YS, ZLW, QX, ZXW, HF and ZL drafted the manuscript. JT and DA provided some writing suggestions. All authors contributed to the article and approved the submitted version. 

Clinical Predictors of Mortality Due to COVID-19 Based on an Analysis of Data of 150 Patients From Wuhan, China

Clinical Features of Patients Infected With 2019 Novel Coronavirus in Wuhan

Cytokine Release Syndrome in Severe COVID-19

Dysregulated Adaptive Immune Response Contributes to Severe COVID-19

Pros and Cons for Use of Statins in People With Coronavirus Disease-19 (COVID-19)

Dexamethasone in Hospitalized Patients With Covid-19

Corticosteroid Treatment in Severe COVID-19 Patients With Acute Respiratory Distress Syndrome

Thymosin Alpha 1 (Talpha1) Reduces the Mortality of Severe COVID-19 by Restoration of Lymphocytopenia and Reversion of Exhausted T Cells

The Effect of Thymosin Alpha1 on Mortality of Critical COVID-19 Patients: A Multicenter Retrospective Study

Risk Factors of Mortality and Contribution of Treatment in Patients Infected With COVID-19: A Retrospective Propensity Score Matched Study

Thymosin Alpha1 Therapy in Critically Ill Patients With COVID-19: A Multicenter Retrospective Cohort Study

Thymosin Alpha1: An Endogenous Regulator of Inflammation, Immunity, and Tolerance

Thymosin Alpha 1 Treatment for Patients With Sepsis

Immunotherapy for Hepatitis B in the Direct Acting Antiviral Era: Reevaluating the Thymosin Alpha1 Efficacy Trials in the Light of a Combination Therapy Approach

Thymosin Alpha1 Represents a Potential Potent Single-Molecule-Based Therapy for Cystic Fibrosis

The Efficacy of Thymosin Alpha 1 for Severe Sepsis (ETASS): A Multicenter, Single-Blind, Randomized and Controlled Trial

Cytokine Storm in COVID-19: Pathogenesis and Overview of Anti-Inflammatory Agents Used in Treatment

Contrast-Enhanced Ultrasound in Patients With COVID-19: Pneumonia, Acute Respiratory Distress Syndrome, or Something Else?

National Health Commission of the People's Republic of China. Diagnosis and Treatment of COVID-19 Guidelines (Trial Version 7)

The SOFA (Sepsis-Related Organ Failure Assessment) Score to Describe Organ Dysfunction/Failure. On Behalf of the Working Group on Sepsis-Related Problems of the European Society of Intensive Care Medicine

APACHE II: A Severity of Disease Classification System

Acute Respiratory Distress Syndrome: The Berlin Definition

Clinical Practice Guidelines for Acute Kidney Injury

Lymphopenia in Severe Coronavirus Disease-2019 (COVID-19): Systematic Review and Meta-Analysis

Patients With COVID-19 in 19 ICUs in Wuhan, China: A Cross-Sectional Study

Potential Interventions for Novel Coronavirus in China: A Systematic Review

Thymosin Alpha-1-Transformed Bifidobacterium Promotes T Cell Proliferation and Maturation in Mice by Oral Administration

Thymosin Alpha1 as a Stimulatory Agent of Innate Cell-Mediated Immune Response

Thymosin Alpha 1 Mitigates Cytokine Storm in Blood Cells From Coronavirus Disease

Dual Effect of Thymosin Alpha 1 on Human Monocyte-Derived Dendritic Cell In Vitro Stimulated With Viral and Bacterial Toll-Like Receptor Agonists

COVID-19 Infection: The Perspectives on Immune Responses

Clinical Characteristics of 138 Hospitalized Patients With 2019 Novel Coronavirus-Infected Pneumonia in Wuhan, China

Immune Modulation With Thymosin Alpha 1 Treatment

We would like to thank all the staff at the five selected hospitals, who were involved in COVID-19 patient care, as well as all the patients who took part in the study and their families.

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.

Conflict of Interest: 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.Publisher's Note: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.Copyright © 2021 Liu, Shen, Wen, Xu, Wu, Feng, Li, Dong, Huang, Guo, Zhang, Chen, Li, Zhu, Du, Liu, Wang, Chen, Teboul, Annane and Chen. This is an openaccess 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.