Protective reactive thymus hyperplasia in COVID-19 acute respiratory distress syndrome

Abstract

Background

Patients with COVID-19 (COVID) may develop acute respiratory distress syndrome with or without sepsis, coagulopathy and visceral damage. While chest CT scans are routinely performed in the initial assessment of patients with severe pulmonary forms, thymus involvement and reactivation have not been investigated so far.

Methods

In this observational study, we systematically scored the enlargement of the thymus and the lung involvement, using CT scans, in all adult patients admitted to the ICU for COVID or any other cause (control group) at one centre between March and April 2020. Initial biological investigations included nasal detection of SARS-CoV-2 ribonucleic acid by polymerase chain reaction (PCR). In a subgroup of 24 patients with different degrees of pulmonary involvement and thymus hypertrophy, plasma cytokine concentrations were measured and the export of mature T cells from the thymus was estimated simultaneously by PCR quantification of T cell receptor excision circles (TRECs).

Results

Eighty-seven patients were studied: 50 COVID patients and 37 controls. Non-atrophic or enlarged thymus was more commonly observed in COVID patients than in controls (66% vs. 24%, p < 0.0001). Thymus enlargement in COVID patients was associated with more extensive lung injury score on CT scans (4 [3–5] vs. 2 [1.5–4], p = 0.01), but a lower mortality rate (8.6% vs. 41.2%, p < 0.001). Other factors associated with mortality were age, lymphopaenia, high CRP and co-morbidities. COVID patients had higher concentrations of IL-7 (6.00 [3.72–9.25] vs. 2.17 [1.76–4.4] pg/mL; p = 0.04) and higher thymic production of new lymphocytes (sj/βTREC ratio = 2.88 [1.98–4.51] vs. 0.23 [0.15–0.60]; p = 0.004). Thymic production was also correlated with the CT scan thymic score (r = 0.38, p = 0.03) and inversely correlated with the number of lymphocytes (r = 0.56, p = 0.007).

Conclusion

In COVID patients, thymus enlargement was frequent and associated with increased T lymphocyte production, which appears to be a beneficial adaptation to virus-induced lymphopaenia. The lack of thymic activity/reactivation in older SARS-CoV-2 infected patients could contribute to a worse prognosis.

Introduction

Less than one year after the beginning of the coronavirus 2019 pandemic (COVID) due to the virus identified as SARS-CoV-2, an increasing number of studies has markedly improved our knowledge of its epidemiology and illustrated the large spectrum of clinical consequences of the infection. Although a large proportion of infected people remains asymptomatic or develops a more or less severe flu-like syndrome, some patients can develop severe forms and require hospitalisation [1, 2]. In France, the infection fatality rate ranges from 0.001% in individuals under 20 years of age to 8.3% in patients over 80 years old [3].

Several studies have identified clinical and biological risk factors for death, which include age, gender, co-morbidities (such as obesity, diabetes, hypertension), D-dimer concentration, low lymphocyte counts and high C-reactive protein [2, 4]. Recent investigations have highlighted the role of increased pro-inflammatory cytokines (cytokine storm) [5, 6], impaired type I interferon responses [7], profoundly altered T cell phenotypes [8] and functional exhaustion of antiviral lymphocytes [9, 10] in the severity of COVID. These findings, similar to those reported in previous pathogenic human coronavirus epidemics due to SARS-CoV-1 and Middle East respiratory syndrome coronavirus (MERS-CoV), reflect the shift from a protective regulated inflammatory response against the virus to pathogenic dysregulated inflammation [11].

In this context, during the initial chest CT scan assessment of adult patients admitted to our intensive care unit (ICU) for COVID, we noticed a previously unreported marked thymus enlargement in many individuals. Although the thymus continues to generate new T lymphocytes in adults, it undergoes progressive physiological involution with age [12, 13]. Enlargement of the thymus region, especially in the elderly, is mostly observed in autoimmune conditions and tumours, or in response to profound lymphopaenia [14], whether or not associated with sepsis. Such lymphopaenia is frequently observed in patients with severe COVID [15]. However, thymic volume is not directly related to the extent of thymic function, which consists of producing and exporting mature T cells into the blood. Accordingly, the most accurate measurement of thymic activity in humans relies on the quantification of T cell receptor excision circles (TRECs), which are by-products of chromosomal rearrangements that occur at T cell receptor loci (TCRA and TCRB, leading to sjTREC and βTRECs, respectively) during T cell development in the thymus. The ratio of both types of TRECs (sj/βTREC), quantified in circulating lymphocytes, is considered to be a surrogate marker of the extent of thymic production of new T cells [16, 17].

The goal of this study was to quantify the enlargement of the thymus in all chest CT scans performed in patients admitted to our intensive care unit in March and April 2020, and to compare the CT scan scores of COVID patients with those of patients hospitalised for other reasons. The clinical condition of the patients, their thymic function and outcomes were then examined according to these scores.

Methods

Patients

This observational study was conducted in all 87 adult patients (50 COVID and 37 controls) hospitalised in the intensive care unit (ICU) of the Clinique Ambroise Paré (Neuilly, France) from the beginning of March to the end of April 2020. During this time period, this ICU, which under normal circumstances principally treats patients admitted for cardiovascular diseases, contributed to the national French program aimed at enhancing the hospitalisation capacity in the ICU for severely affected SARS-CoV-2-infected patients.

According to French national regulations, systematic written informed consent was requested from all patients upon admission, for the potential anonymous use in clinical research studies of clinical and paraclinical data obtained during hospitalisation.

Thoracic CT scans

All patients had CT scans at admission. Thoracic CT scans were performed for the initial evaluation of pulmonary conditions upon admission. Specific settings were used to examine the thymus area. All images were independently reviewed and classified by two radiologists. In case of disagreement, they reviewed the scans together and were asked to give a final common classification. The classification of pulmonary parenchymal images and of thymus aspects is shown in Tables 1 and 2 as well as in Fig. 1.

Table 1 CT scan classification of COVID-associated pneumopathy

Table 2 CT scan classification of thymus aspects

Fig. 1

Chest CT scan images. Representative chest CT scan images of the thymic area of patients hospitalised for COVID, illustrating the classification of thymus aspects depicted in Table 2. Arrows indicate the different types of thymus hyperplasia

Laboratory procedures

Initial clinical laboratory investigations included a complete blood count, serum biochemical tests including liver and kidney function, creatine kinase, lactate dehydrogenase, electrolytes, fibrinogen, C reactive protein (CRP) and a coagulation profile. SARS-CoV-2 infection was confirmed by nasal viral RNA detection using reverse transcriptase polymerase chain reaction (RT-PCR) following the procedure recommended by the manufacturer and/or clinical history associated with the notion of plausible contamination by a proven infected contact.

Plasma concentrations of 23 cytokines (IL-1β, IL-4, IL-6, IL-7, IL-8, IL-10, IL-12p70, IL-15, IL-17A, IL-17E, IL-17F, IL-21, IL-22, IL-23, IL-33, TSLP, GM-CSF, TNFα, MDC, MIP-1α, MIP-1β, MIP-3β and SDF-1) were measured on the same day in 24 hospitalised patients, including 18 COVID patients and 6 control individuals with different degrees of pulmonary involvement and thymus hypertrophy, using a multiplex ELISA (U-plex technique; Meso Scale Diagnostics, Rockville, Maryland) according to the manufacturer’s instructions.

The capacity of the thymus to export mature T cells into the bloodstream was estimated in the same patients through PCR quantification of two types of T cell receptor excision circles (sjTRECs and βTRECs) in peripheral blood mononuclear cells and by the calculation of the sj/βTREC ratio, a good estimate of the thymic output ratio, as described previously [16,17,18].

Statistical analysis

Descriptive statistics were used to present baseline characteristics for the total population and for the various groups. The normality of continuous data was assessed using the Shapiro–Wilk test. Continuous variables are presented as mean ± SD or median [25–75% interquartile range], according to their distribution, and categorical variables are presented as the number of patients in each category and the corresponding percentages. Missing values concerned four patients for leukocyte counts, two patients for haemoglobin levels, one patient for CRP, and 34 patients for D-dimer and fibrinogen. They were not replaced and were considered as empty cells in all calculations. Quantitative variables were compared using Student’s t test or the Mann–Whitney U test according to the normality of the distribution. Categorical variables were compared using Pearson tests. Kruskal–Wallis tests were used to compare three groups. Bivariate (correlation coefficient) analyses were performed to determine the relationships between the variables. The Bonferroni correction was used for multiple comparison adjustments. JMP software was used for the statistical analysis, and a value of p < 0.05 was considered to be statistically significant.

Results

Among the 87 patients enrolled in the study, 37 were considered as non-infected by SARS-CoV-2 (control group). Fifty patients had confirmed COVID (COVID group). The control group (non-COVID) was composed of patients admitted during the same period of time for other pathologies and who had a CT scan: medical cardiology (16), cardiac surgery (11), urology (5) and other (5). They had no clinical or radiological signs of COVID, and their PCR test was negative. A descriptive analysis of the two groups at the time of the first CT scan is summarised in Table 3.

Table 3 Descriptive analysis of the investigated patients

Thymus enlargement characterises COVID patients

The thymus aspects were classified into three groups (Table 2). Seventeen patients displayed fatty atrophy of the thymus (the most common in senior adults) and were classified as subgroup 0. Fourteen patients with a homogeneous non-fatty thymus, commonly observed in young adults, or a fatty thymus area associated with micronodules, were classified as subgroup A. Nineteen patients showing various levels of thymus hyperplasia or a pseudo-tumoral thymus represented subgroup B. Thymus persistence (subgroup A) or enlargement (subgroup B) was more frequent in the COVID group than in the control group (66% vs. 24%, p < 0.0001), except in patients over 80 years of age (Fig. 2, p < 0.05).

Fig. 2

Thymic enlargement in COVID and control groups according to age categories. Thymic scores were determined through CT scan imaging in COVID (n = 50, red symbols) and COVID (n = 37, blue symbols) patients and displayed as a function of age categories. *p < 0.05; **p < 0.001, NS, not significant

COVID patients with a persistent or enlarged thymus were younger than patients with normal thymic fatty atrophy (75.4 ± 12.6 years in subgroup 0 as compared to 59.3 ± 12.7 years and 55.3 ± 16.2 years in subgroups A and B, respectively; p < 0.0001, Table 3), and displayed, on average, more severe pulmonary involvement (4 [3,4,5] vs. 2 [1.5–4]; p = 0.01, Fig. 3). They were also less frequently hypertensive (p < 0.05), developed less frequent renal failure (p < 0.001) and had a lower mortality rate (7.1% and 5.2% vs. 41.2%, p < 0.001; Table 4 and Fig. 4). The sex ratio, lymphopaenia and D-dimer concentration were not different among the subgroups, even in younger patients.

Fig. 3

Lung damage in COVID patients, according to thymus enlargement group. Lung CT scan score is presented for each patient, according to individual thymus CT scan. Each dot represents an individual patient. Median values with the 25–75% interquartile ranges are shown; **p = 0.01 between groups

Fig. 4

Survival as a function of thymic CT scan score. Thymus CT scan score is presented for each patient in survivor and non-survivor groups. Statistical differences between groups are shown on top; ***p < 0.0001

Table 4 Mortality in the COVID group, univariate analysis

Enhanced thymic function in COVID patients

As expected, sjTREC values decreased as a function of age in both COVID patients and controls. sjTREC values were slightly higher in COVID individuals and correlated with the CT scan thymic score (p = 0.02, n = 24). In COVID patients, the sj/βTREC ratio was higher than in controls, particularly in older individuals, consistent with enhanced thymic production of new lymphocytes (Fig. 5). The sj/βTREC ratio was also correlated with the CT scan thymic score.

Fig. 5

Enhanced thymic function in COVID patients. The sj/βTREC ratio, which reflects the extent of thymic output, is represented as a function of age (left) or as a function of CT scan thymic score (right). Statistical differences between groups (ANCOVA using age as a concomitant variable), nonparametric correlation Pearson’s r and associated p value (one tailed) are indicated

Increased IL-7 plasma concentrations in COVID patients

Among the cytokines implicated in the immune response, IL-6, IL-10 and MIP-1α concentrations were correlated with COVID intensity as measured by chest CT scan (r = 0.47, r = 0.39 and r = 0.49; p = 0.024, p = 0.056 and 0.019, respectively; Fig. 6). Interestingly, IL-6 and IL-10 were not significantly elevated in COVID patients than in controls. In contrast, we found significantly higher concentrations of IL-7 in the plasma of COVID patients compared to uninfected controls (Fig. 7). Surprisingly, contrary to previous observations in lymphopaenic patients, IL-7 plasma levels were directly proportional to blood lymphocyte counts in both groups (Fig. 8).

Fig. 6

Plasma concentrations of three cytokines as a function of chest CT scan score. Plasma concentrations of IL-10, IL-6 and MIP-1α were measured in patients hospitalised for COVID (n = 18) and are presented as a function of chest CT scan scores. Nonparametric correlation Pearson’s r and associated p value (one tail) are indicated

Fig. 7

IL-7 concentrations in the plasma of control and COVID patients. IL-7 plasma concentration was measured in patients hospitalised for COVID (n = 18, red symbols) and controls (n = 6, blue symbols). Statistical significance of observed differences between groups are shown on top (Mann–Whitney U test)

Fig. 8

IL-7 concentrations correlated with blood lymphocyte counts. IL-7 plasma concentration was measured in patients hospitalised for COVID (n = 18, red symbols) and controls (n = 6, blue symbols) and presented as a function of blood lymphocyte counts. Nonparametric correlation Pearson’s r and associated p value (one tail) are indicated

Discussion

We found that 66% of COVID patients hospitalised in our ICU displayed a surprisingly large thymus. This thymus enlargement was associated with more severe pulmonary involvement but less mortality (8.6% vs. 41.2% in the control group). These findings were not mentioned in previous reports on features associated with COVID. The emergency context in which scans of patients hospitalised in intensive care for severe forms of COVID pneumonia were performed probably contributed to underestimating the changes occurring in the thymic area. Indeed, analysis of the CT scan images of the pulmonary parenchyma requires specific settings (for example, adjustment of grey levels), which are not suitable for mediastinum analysis, causing a lack of contrast between the various mediastinal tissues.

In adults, hyperplasia of the thymus is unusual and can be observed in a limited number of pathological conditions: auto-immune diseases such as myasthenia gravis [19] and Graves-Basedow disease [20], after high-dose chemotherapy associated or not with autologous stem cell transplantation [21, 22], in lymphopaenic patients infected by human immunodeficiency virus (HIV) with maintained naïve T cell counts [23] or after antiretroviral therapy [24]. In the setting of severe T cell depletion, caused by HIV infection or cytoreductive transplant or chemotherapy regimens, thymus hyperplasia is critical for the restoration of peripheral T cell populations [22].

In patients with severe forms of COVID, severe lymphopaenia is commonly observed [15]. In our study, it was associated with overall reduced survival and inversely correlated with the intra-thymic proliferation of T cell precursors. Together, these data indicate that the increase in thymic mass observed in COVID patients is a beneficial adaptation to virus-induced lymphopaenia. Associated with increased thymic production and at least partly triggered by enhanced IL-7 levels, this adaptation appears to decline in patients older than 80 years of age, possibly contributing to the higher mortality observed in more senior patients. Severe COVID has been shown to be associated with massive infiltration of macrophages in the lungs, while recovery is dependent on T cell responses to the virus [25]. A similar progressive reduction in T cell precursor proliferation with age was previously reported in healthy HIV individuals [16] and in a systematic study of 1,000 healthy individuals showing that age and sex strongly affect thymic function [26]. Moreover, while potential genetic factors are currently under investigation in COVID patients who develop a severe form of the disease despite the absence of risk factors, it should be noted that thymopoiesis is genetically determined in both mice and healthy human adults [26, 27].

We also confirmed that the concentrations of IL-6, IL10 and MIP-1α, which are considered to be key factors of pulmonary involvement, were effectively correlated with the degree of lung involvement in COVID patients. Surprisingly, despite the already demonstrated relationship with disease intensity, IL-6 was not significantly enhanced in the COVID patient group as compared to control individuals. However, it should be noted that most of our COVID patients were not in a very severe condition, and the only patient with a very high IL-6 plasma concentration died. Moreover, our control patients cannot be considered as ‘healthy’ since they were all hospitalised in the ICU; the condition of these patients may have caused some lung inflammation, independently of the original pathological condition [28]. In contrast, we observed a significant increase in the concentration of IL-7 in the plasma of COVID patients compared to uninfected controls. A similar observation has already been reported in patients with severe COVID [29]. Unlike previous observations showing higher plasma IL-7 levels in lymphopaenic patients, plasma IL-7 levels were directly proportional to the number of blood lymphocytes in COVID patients. This observation indicates that, during acute SARS-CoV-2 infection, high plasma levels of IL-7 were not due to its low consumption by T cells, as previously suggested [30], but most likely originated from active overproduction, as previously evidenced in acute SIV infection in rhesus macaques [31]. Interestingly, IL-7 is a critical cytokine for the survival, proliferation and differentiation of immature thymocytes [32] and is a cofactor for rearrangement of the T cell receptor beta chain during early T cell development [33]. It is thus possible that the enhanced production of IL-7 in COVID patients participates in the observed thymus enlargement and increased thymic function. Unlike thymic production, the actual size of the thymus in the CT scans of patients belonging to subgroups A and B of the thymus classification was not inversely correlated with lymphopaenia. Nodular and/or minimal areas of hyperplastic thymus might be sufficient to restore an adequate number of peripheral lymphocytes.

The observation that the pulmonary involvement of COVID patients quantified by the CT scan score was significantly higher in patients with an enlarged thymus appears counter-intuitive based on the above data and hypotheses. A plausible explanation of this apparent paradox is that, in patients with an ‘activated productive’ thymus, the influx of immune cells into the infected lungs is different from that which occurs in patients with an ‘unreactive’ thymus. This hypothesis is supported by a recent study in which bronchoalveolar lavage fluid cells from COVID patients were characterised using single-cell RNA sequencing. Monocyte-derived inflammatory macrophages were found in severe COVID pneumonia, contrasting with clonal expansion of CD8 + effector T cells in mild cases, indicating a role for CD8 + T cells and the adaptive immune response in the clearance of SARS-CoV-2 [34].

This observational study, conducted in an emergency context, has some limitations. TRECs could only be evaluated in patients who were still hospitalised; some immunological investigations in the blood or bronchoalveolar lavage fluid had not been anticipated, and patient follow-up was limited in time. Prospective studies will be necessary to further characterise the role of thymic function in the control of SARS-CoV-2 infection and the involved molecular and cellular mechanisms.

Conclusion

In response to SARS-CoV-2 infection, thymic reactivation is frequent and seems to be a good prognostic factor because it indicates the activation of compensatory mechanisms in the context of lymphopaenia, contributing to an efficient adaptive immune response in the lungs. A CT scan examination of the thymic area with appropriate settings is recommended in all COVID patients with pulmonary involvement.