Abstract
Objective SARS-CoV-2 infection can lead to life-threatening clinical manifestations. Patients with cardiovascular disease (CVD) are at higher risk for severe courses of COVID-19. However, strategies to predict the course of SARS-CoV-2 infection in CVD patients at hospital admission are still missing. Here, we investigated whether the severity of SARS-CoV-2 infection can be predicted by analyzing the immunophenotype in the blood of CVD patients.
Approach and Results We prospectively analyzed the peripheral blood of 94 participants, including CVD patients with acute SARS-CoV-2 infection, uninfected CVD patients, and healthy donors using a 36-color spectral flow cytometry panel. Clinical assessment included blood sampling, echocardiography, and electrocardiography. Patients were classified by their ISARIC WHO 4C-Mortality-Score on the day of admission into three subgroups of an expected mild, moderate, or severe course of COVID-19. Unsupervised data analysis revealed 40 clusters corresponding to major circulating immune cell populations. This revealed little differences between healthy donors and CVD patients, whereas the distribution of the cell populations changed dramatically in SARS-CoV-2-infected CVD patients. The latter had more mature NK cells, activated monocyte subsets, central memory CD4+ T cells, and plasmablasts than uninfected CVD patients. In contrast, fewer dendritic cells, CD16+ monocytes, innate lymphoid cells, and CD8+ T cell subsets were detected in SARS-CoV-2-infected CVD patients. We identified an immune signature characterized by low frequencies of MAIT and intermediate effector CD8+ T cells in combination with a high frequency of NKT cells that is predictive for CVD patients with a severe course of SARS-CoV-2 infection on hospital admission.
Conclusion Acute SARS-CoV-2 infected CVD patients revealed marked changes in abundance and phenotype of several immune cell populations associated with COVID-19 severity. Our data indicate that intensified immunophenotype analyses can help identify patients at risk of severe COVID-19 at hospital admission, improving clinical outcomes through specific treatment.
Highlights
Patients with cardiovascular disease are at higher risk of severe courses of COVID-19 and may exhibit an altered immune response
Unsupervised data analysis revealed that patients with cardiovascular disease and SARS-CoV-2 infection showed significant changes in the abundance and the phenotype of various immune cell populations
We identified a disease-related immune signature in patients with cardiovascular disease and SARS-CoV-2 infection associated with the severity of COVID-19
Intensified immunophenotyping helps to identify cardiovascular patients at risk of a severe course of COVID-19 already at the early stages of the disease and might thereby improve clinical outcomes and specific treatment of COVID-19
1. Introduction
Coronavirus disease 2019 (COVID-19) refers to a broad spectrum of clinical manifestations caused by infection with the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2). The severity of the disease is related to risk factors such as age, sex, and pre-existing comorbidities1, 2 that correlate with the immune response during acute infection 1, 2. SARS-CoV-2 infection ranges from asymptomatic to fatal courses. In patients recovering from a non-severe infection, the immune system responded to SARS-CoV-2 infection with robust, broad-based, and transient regulatory features 3–6. In contrast, severe COVID-19 is characterized by hyperactivation of innate but also adaptive immune cells, an exuberant cytokine response, and high titers of SARS-CoV-2-specific antibodies 4, 5, 7, 8.
Common complications in hospitalized patients infected with SARS-CoV-2 include pneumonia, sepsis, acute respiratory distress syndrome (ARDS), and respiratory failure 9–11. The pathophysiology of SARS-CoV-2 is characterized by an early production of pro-inflammatory cytokines (e.g., tumor necrosis factor (TNF), IL-6, and IL-1β), often resulting in hyperinflammation 12. This so-called cytokine storm increases the risk of vascular hyperpermeability and, if persistent, multi-organ failure and eventual death 11.
Patients with cardiovascular disease (CVD) are critically susceptible to more severe courses of SARS-CoV-2 infection, which may lead to life-threatening complications like cardiac and pulmonary damage as the most common complications 13, 14. For example, an enhanced pro-inflammatory and pro-thrombotic immune response, characteristic of SARS-CoV-2 infected patients with pre-existing CVD, can trigger myocarditis or acute coronary syndrome with subsequent congestive heart failure 15–19. Furthermore, CVD patients are prone to ARDS and progressive respiratory failure and along with an increased risk of a pulmonary embolism due to infection-associated coagulopathy, which may cause not only acute right heart failure but also disseminated intravascular coagulation 19. These clinical scenarios explain the increased rate of unfavorable clinical outcomes like organ failure, admission to the intensive care unit (ICU) with rapidly progressive respiratory failure, and mortality in CVD patients with SARS-CoV-2 infection 14, 15, 17
CVD is characterized by alterations in inflammatory mediators like C-reactive protein (CRP), pro-inflammatory cytokines and chemokines, platelet and monocyte activation, and enhanced expression of adhesion molecules on endothelial and immune cells are critical players of atherogenesis and -progression 20–22. Up-regulation of the involved cytokines and chemokines and platelet activation recruits inflammatory cells like monocytes, macrophages, and dendritic cells to the arterial wall causing atherosclerotic lesions 21. In addition, pro-inflammatory adaptive immune cells like Th1 and Th17 cells are also critically involved in CVD, especially myocardial infarction, myocarditis, and heart failure 23–25. Due to a chronically enhanced pro-inflammatory immune response, patients with CVD are at risk for a severe course of COVID-19. However, strategies to identify these high-risk patients with CVD early in COVID-19 are still lacking. Thus, there is still an urgent clinical need to understand the complexity of the innate and adaptive immune system dysfunctions that underlie severe or even fatal COVID-19 in this subgroup of CVD patients to provide the best possible care upon hospitalization and improve clinical outcomes. In the present study, we aimed to identify the critical immune system components to predict a severe course of COVID-19 in CVD patients. Accordingly, we performed in-depth analyses of innate and adaptive immune cells and cytokines in the peripheral blood of CVD patients with and without acute SARS-CoV-2 infection. Our study revealed that a complex but specific immune signature is associated with the severity of COVID-19 in patients with CVD and can predict the course of the disease already upon first admission to the hospital.
2. Methods
Study design, participants, and assessment of clinical parameters
From March to April 2020, we prospectively studied a consecutive cohort of 94 participants at the Department of Cardiology and Angiology of the University Hospital Tübingen, Germany. Of these, 37 consecutive patients with pre-existing CVD and symptomatic, acute SARS-CoV-2 infection (CVD+SARS-CoV-2) were admitted to our emergency department. Twenty patients with pre-existing stable CVD without any infections were matched for the group of CVD+SARS-CoV-2 patients. 37 healthy donors (HD) served as controls (Table 1 and Tables S1-2). All patients underwent clinical and cardiac assessment, including echocardiography, electrocardiography, concomitant medication, comorbidities, and blood sampling for routine laboratory parameters within 12 hours of admission. SARS-CoV-2 infection was diagnosed by RNA detection from nasopharyngeal secretions with real-time reverse transcriptase polymerase chain reaction. Pre-existing CVD was defined as stable coronary artery disease, which had been determined by coronary angiography and diagnosed when there was luminal stenosis of one or more coronary vessels >25-50% diameter before hospital admission. Respiratory failure was defined by a Horovitz Index (HI) ≤ 200 mmHg. 23
The inclusion criteria for our study were confirmed CVD with or without SARS-CoV-2 infection and an age ≥18 years. Exclusion criteria were other viral or bacterial infections and malignancies. The study was approved by the local ethics committee (240/2018BO2) and complied with the Declaration of Helsinki and the Good Clinical Practice Guidelines on the approximation of the laws, regulations, and administrative provisions of the Member States relating to the implementation of good clinical practice in the conduct of clinical trials on medicinal products for human use. Written informed consent was obtained from each patient.
N terminal-pro-B-type natriuretic peptide (NT-pro-BNP, >300 ng/L), high sensitive troponin I (hs TNI, >37 ng/L), and C-reactive protein (CRP, >0.5 mg/dL) were classified as elevated laboratory markers of myocardial and inflammatory distress. Echocardiographic parameters included left and right ventricular function, right ventricular dilatation, presence of tricuspid valve regurgitation, and pericardial effusion according to current guidelines.24, 25
Isolation of peripheral blood mononuclear cells
Blood samples of patients with CVD (including SARS-CoV-2 infection) and HD were collected in CPDA monovettes and were processed within 4 hours after blood collection. Peripheral blood mononuclear cells (PBMCs) were isolated using SepMate tubes (Stem Cell Technologies) according to the manufacturer’s instructions by density gradient centrifugation at 1200 x g for 10 min at room temperature with 25 mL cell suspension (blood diluted with PBS 1:1) stacked on 15 mL Biocoll separation solution (Biochrom). The clear supernatant containing the PBMCs was decanted and washed three times with PBS+2% FCS. PBMCs were frozen as aliquots of 1-2×107 cells in RPMI1640 containing 20% FBS and 10% DMSO at −150°C until further use.
Flow Cytometry staining
PBMCs were thawed in a water bath at 37°C, followed by adding 10 ml RPMI (Sigma) containing 50 KU DNaseI (Merck). Cells were centrifuged at 400 x g for 5 minutes at room temperature. The supernatant was discarded, and the cells were resuspended in 5 ml RPMI containing 200 KU DNaseI and incubated for 20 minutes at 37°C. The cell numbers were determined by trypan blue exclusion using a Neubauer counting chamber. Cells were washed once with PBS, and 3×106 cells were stained with LIVE/DEAD Fixable Blue Dead Cell Stain Kit (Thermo Fisher Scientific) according to the manufacturer’s instructions to exclude dead cells (final dilution 1:1000). All following washing and incubation steps were performed with Cell Staining Buffer (BioLegend). Blocked Cells were incubated with Human True StainFcX (BioLegend) at room temperature for 10 minutes to prevent unspecific binding of antibodies to Fc receptors. After that, extracellular staining was performed for 1 hour at 4°C using the antibodies listed in Table S3, in a final volume of 100 µl containing 5 µl True-Stain Monocyte Blocker (BioLegend) to block non-specific binding of PE/Dazzle594, PE/Cy5, PE/Cy7, APC/Cy7, and APC/Fire750 to monocytes. Cells were washed twice with staining buffer, and at least 1.5×106 cells were acquired using an Aurora (Cytek Bioscience) with the SpectroFlow software. Unmixing was performed in SpectroFlow using cell- or bead-based single-stain controls and unstained cells for autofluorescence subtraction.
Flow Cytometry data analysis
We performed a classical gating strategy (Supplemental Figure 1) for quality control of the 36-color panel and unsupervised data analysis described below using OMIQ software (Dotmatics, Boston, MA, USA). First, the compensation and the scaling were set, followed by a subsampling to 1.5×106 events/sample. FlowAI was run (flow rate: second fraction 0.1; alpha 0.01; and dynamic range: both limits with negative value removal limit 1) including time, all fluorescent channels, and autofluorescence followed by gating and subsampling on 1×106 flowAI passed cells/group (HD, CVD or CVD+SARS-CoV2). The data were normalized using fdaNorm for all fluorescent channels, followed by adjusting the scaling where needed and gating on and subsampling of living 5×105 CD45+ cells per group. Subsequently, dimensionality reduction analysis was performed using Uniform Manifold Approximation and Projection (UMAP) to visualize the different sub-populations of the cells.25 UMAP settings were as follows: all files used, all fluorescent parameters were used except Live/Dead, CD45 and autofluorescence, Neighbors = 80, Minimum Distance = 0.7, Components = 2, Metric = Euclidean, Learning Rate = 1, Epochs = 200, Random Seed = 5733, Embedding Initialization = spectral. Following the UMAP26, 27 analysis, FlowSOM27 was run to cluster the data. FlowSOM settings were as follows: all files used, clustering features CD16, CD11c, CD56, CD8, CCR7, CD123, IgD, CD3, CD20, IgM, IgG, CD28, CD141, CD57, CD14, TCRγδ, CD25, CD4, CD27, CD1c, CD19, CD127, HLA-DR, CD38, umap_1, umap_2, xdim = 25, ydim = 25, # of training iterations = 10, Distance Metric = euclidian, consensus metaclustering with k = 100, Random Seed = 3919. A heatmap was generated with the metaclusters obtained from FlowSOM and clustered hierarchically using the samples from HD with the medians of all surface markers except for LIVE/DEAD, CD45, and autofluorescence with a euclidean row distance and a ward row leakage to assign the metaclusters to the cell populations.
Determination of plasma levels of cytokines/chemokines (LEGENDPlex)
Two LEGENDPlex (Inflammation Panel 1 and Pro-inflammatory Chemokine Panel; BioLegend, San Diego, California, USA) were performed to quantify the concentrations of several chemokines and cytokines in human plasma. Sixty-nine frozen plasma samples were analyzed, consisting of 27 HD, 15 CVD, and 27 SARS-CoV-2-infected CVD patients. The assays were performed according to the manufacturer’s manual. A FACS Lyric (BD Biosciences, Franklin Lakes, New Jersey, USA) was used for the measurement. Data analysis was performed with the LEGENDPlex Data Analysis Software (BioLegend, San Diego, California, USA).
Statistical analysis
We determined clinical and laboratory baseline characteristics in relation to measured immune cell phenotypes, marker expression, and clinical outcome. Continuous, not normally distributed variables are expressed as median with standard deviation and were compared using unpaired two-tailed Mann–Whitney U test for two-group comparison or Kruskal Wallis non-parametric test with Dunn’s post-test for three-group comparison. Categorical data are presented as total numbers and proportions and were analyzed by chi-squared test. Correlation analysis was performed by Spearman rank correlation coefficient r. To summarize, correlations of essential parameters in CVD patients with and without SARS-CoV-2 infection matrices were generated using Rstudio “Corrplot,” displaying correlations of cytokines with or without cell populations from the unsupervised data analysis. Spearman’s ρ is colored, and color intensity and size are plotted proportionally to correlation coefficients.
The “CombiROC” package was used in Rstudio to determine the optimal combination of cell populations to differentiate mild from severe or moderately SARS-CoV-2-infected CVD patients.28, 29 Orthogonal partial least-square discriminant analysis (OPLS-DA) was performed using the “MetaboAnalyst” package in RStudio.
Comparisons were considered statistically significant if the two-sided p-value was <0.05. Statistical analysis was performed using Prism Software Version 9.4.1 (GraphPad).
Results
Clinical characteristics of patients with cardiovascular disease and symptomatic acute SARS-CoV-2 infection
We prospectively studied a cohort of 94 participants from February to April 2020. This cohort consisted of 20 patients with pre-existing CVD and 37 with pre-existing CVD and symptomatic acute SARS-CoV-2 infection. 37 HD served as controls. The baseline characteristics and demographics of the overall cohort are given in Table 1. The population’s median age was 58 (IQR 42–74) years, and 45 (47.9%) patients were men. Detailed information on every subgroup of the cohort is given in Table S1. Out of 37 CVD+SARS-CoV-2 patients, 20 showed respiratory failure with HI ≤ 200 mmHg, while 11 (29.7%) were admitted to the ICU due to progressive respiratory, circulatory, or multi-organ failure. CVD+SARS-CoV-2 patients were further stratified by their ISARIC-WHO-4C-Mortality-Score into a group of expected mild, moderate, or severe course of COVID-19, patients’ characteristics stratified by ISARIC-WHO-4C-Mortality-Score are depicted in Table S2. Patients with expected mild COVID-19 were significantly younger (p=0.016) than those with an expected moderate or severe disease course. Interestingly, cardiovascular risk factors were equally distributed among the three groups, whereas laboratory parameters like lymphocyte count, CRP, and interleukin-(IL-)6 levels showed significant alterations.
Characterization of immune cell subsets in peripheral blood using a 36-color spectral flow cytometry panel
Frozen peripheral blood mononuclear cells (PBMCs) from the described patient cohort were stained with a 36-color antibody panel (Suppl. Figure 1) related to a previously published panel 30 and measured with a spectral flow cytometer (Aurora, Cytek) to characterize all common immune cell subsets (Fig. 1A). Data analysis was performed classically by manual gating (Suppl. Figure 2) as well as by unsupervised analysis using Uniform Manifold Approximation and Projection (UMAP)26 for dimensionality reduction followed by clustering using self-organized map (FlowSOM)27 (Figure 1C). Already in the UMAP visualization differences were seen between the concatenated files of the different groups of the cohort (Figure 1B). FlowSOM clustering was performed with 100 metaclusters (MCs) to identify cell populations with lower abundance, such as cDC1. 13 of these 100 MCs were excluded for further analysis as the percentage of cells in these MCs was less than 0.01. FlowSOM clustering in combination with UMAP plots displaying the expression of one lineage marker on the color-coded x-axis (Suppl. Figure 3) and a clustered heatmap showing the median expression of the analyzed markers (Figure 1D) was used to assign the obtained 87 clusters to distinct cell populations (Figure 1C). Finally, the 36-color staining followed by unsupervised data analysis revealed the discrimination of 40 manually assigned specific cell populations reflecting the manual gating strategy. These cell populations include B cell (lilac), CD4+ (orange) and CD8+ T cell (blue), NK cell (green), monocyte (pink), and dendritic cell subsets (cDCs, violet), innate lymphoid cells (ILCs, light green) and MAIT cells (turquoise) as well as basophils (yellow) and also neutrophils (purple).
SARS-CoV2-infected CVD patients showed significant differences in the distribution and the phenotype of immune cell populations compared to uninfected CVD patients
Next, we assessed the differences in the abundance of the 87 MCs obtained by FlowSOM (annotated to 40 different cell populations as described in Fig. 1) in PBMCs from HD and CVD patients with and without SARS-CoV2 infection (Suppl. Figure 4). Comparison of HD and uninfected CVD patients revealed a reduced frequency of clusters assigned to γδ T cells, MAIT cells, naïve CD8+ T cells, early-like effector memory CD4+ T cells (CD4+ early-like Tem), CD4-CD8- T cells, and naïve B cells in CVD patients. In contrast, clusters representing more innate inflammatory cells like CD14+ monocytes, neutrophils, mature NK cells, and CD8+ central memory T cells (CD8+ Tcm) were increased in frequency in CVD patients compared to HD (Suppl. Figure 5).
During symptomatic, acute SARS-CoV-2 infection in hospitalized CVD patients, we observed significant changes in their immune signature compared to uninfected CVD patients, particularly a greater proportion of mature NK cells (MC33), CD14+CD16+ monocytes (MC11, MC35), activated populations of CD14+CD45RA+ monocytes (MC27, MC50, MC63), CD4+ Tcm cells (MC97) expressing CD38+, CD8+ early Tem cells (MC60) expressing HLA-DR+CD38+PD-1+ as well as plasmablasts (MC31, MC39) (Figure 2B). In contrast, reduced proportions of ILCs (MC62), dendritic cell subsets cDC1 (MC14) and cDC2 (MC13), CD16+ monocytes (MC23, MC34), less activated/mature CD14+ monocytes (MC28, MC45, MC59), as well as CD4+CD8+ T cells (MC76), naïve, central memory, and early effector-memory CD8+ T cells (MC67, MC68, MC61) and CD4+ Tregs (MC100) were observed in SARS-CoV2-infected compared to uninfected CVD patients (Figure 2). To complement our unsupervised analysis, we applied manual gating (Suppl. Figure 2) for immune cell subsets contributing to SARS-CoV-2 infection. Similar to the unsupervised data analysis results, we found increased numbers of CD14+CD16+ intermediate monocytes and plasmablasts and, although not found in the unsupervised data analysis, increased numbers of CD4+ Temra cells in SARS-CoV2-infected CVD patients (Suppl. Figure 6). The numbers of ILCs, dendritic cell subsets pDCs, cDC1, and cDC2 (not statistically significant), CD16+ monocytes, as well as naïve, central memory, and early effector-memory CD4+ and CD8+ T cells and CD4+ Tregs, were reduced in SARS-CoV-2 infected compared to uninfected CVD patients. Moreover, fewer numbers of T helper 17 (Th17) and follicular T helper (Tfh) cells and γδ T cells were detected upon SARS-CoV-2 infection.
In addition to changes in the proportion of immune cells, we also observed phenotypical changes in SARS-CoV-2 infected CVD patients compared to uninfected CVD patients (Suppl. Figure 7). In particular, an overall activation of innate immune cells was observed. For example, the expression of CD38, CD95, and CCR5 was significantly increased on NK cell subsets of SARS-CoV-2 infected CVD patients. Plasmacytoid DCs, important in viral infections due to their type I interferon response 31, also showed a significantly more activated phenotype by increased expression of CCR7, CD38, and CD95. Similar effects were observed for basophils and innate lymphocyte cells (ILCs) with increased expression of CD38, CD95, and CXCR3 or CD45RA, CD28, and CD25, respectively.
In contrast, the cDC subsets from SARS-CoV-2 infected CVD patients showed less expression of HLA-DR and CD38 but increased expression of CCR5, CCR7, CD11b, CD45RA, and IgG (cDC2s) (Suppl. Figure 7), indicating cellular activation and migration, but also an association with a reduced antigen presentation capacity. Regarding the monocyte clusters, we observed an increased migratory and activated phenotype in SARS-CoV-2 infected CVD patients, determined by the high expression of CCR7 in all monocyte subsets. However, monocyte subsets from these patients revealed a less activated phenotype than their uninfected CVD counterparts, shown by significantly less expression of HLA-DR, CD38, CD11b (on CD14+ monocytes), and CD45RA. Most strikingly, we found an enhanced expression of the thrombin receptor CD141 on monocyte clusters from SARS-CoV-2 infected patients, possibly aiming toward impaired blood coagulation 32.
Concerning the adaptive immune cell populations, a more inhibitory phenotype was observed in SARS-CoV2-infected CVD patients compared to uninfected (Suppl. Figure 7). The only exception was the plasmablasts. Plasmablasts revealed reduced expression of the apoptotic marker CD95 combined with increased CD38 expression, indicating a more functional and long-living phenotype upon SARS-CoV2 infection. Similar to DCs, a phenotype with decreased activation was observed in the B-cell subsets of SARS-CoV2-infected CVD patients. This was evidenced by lower CD38, HLA-DR, CCR6, CD45RA, and IgD expression. CD8+ and most CD4+ T cell subsets showed an overall increase in CD38 expression, indicating activation upon SARS-CoV2 infection. However, the inhibitory molecule PD-1 on the cell surface of CD4+ and CD8+ naïve and Tcm cells, CD4+ Tregs, and HLA-DR+CD8+ activated T cells was increased, which was accompanied by less expression of the co-stimulatory molecule CD28 on CD8+ T cell subsets. CXCR3 expression was also reduced on several CD4+ and CD8+ T cell subsets of SARS-CoV-2 infected compared to uninfected CVD patients. The latter two observations hint towards an impaired migratory capacity and full T cell activation of CD4+ and CD8+ T cell subsets. In summary, our data show an infection-induced activation of innate immune cells in combination with reduced numbers and inhibitory phenotypes of DCs and adaptive immune cells in infected compared to uninfected CVD patients.
Chemokine and cytokine profiling showed significant differences between SARS-CoV-2 infected and uninfected CVD patients
In addition to spectral flow cytometry, 25 common cytokines and chemokines were analyzed in plasma samples from HD and CVD patients without and with SARS-CoV-2 infection, the latter being further subdivided into groups with expected mild and severe course of COVID-19 based on the ISARIC-WHO-4C-Mortality-Score. In line with other studies, the pro-inflammatory cytokines IL-6 and IL-18 were significantly increased in the plasma of SARS-CoV-2 infected compared to uninfected CVD patients and HD (Suppl. Figure 8A) 33, 34. Furthermore, IL-6 was also elevated in CVD patients compared to HD and significantly increased in severe compared to mild SARS-CoV2-infected CVD patients (Suppl. Figure 8A&B). Interestingly, TNF, IL-1b, IL-12p70, IL-23, and IL-33 were lower in SARS-CoV-2 infected CVD patients than uninfected CVD patients and HD, indicating impaired T helper cell differentiation in this subgroup. In contrast, the chemokines CCL2 and the CXCR3 ligands CXCL9, CXCL10, and CXCL11 were markedly elevated during SARS-CoV-2 infection. Differences between mildly and severely affected CVD patients were observed; in particular, a significant increase of IL-6 and IL-8, whereas CCL17, responsible for the recruitment of T cells35, was significantly less abundant in severely infected CVD patients (Suppl. Figure 8B). Correlation analysis revealed that predefined cell clusters and chemokine/cytokine profiling showed significant associations among the subgroup of SARS-CoV-2 infected CVD patients and CVD patients without any infections, as depicted in Suppl. Figure 8C.
An immune signature of SARS-CoV-2 infected CVD patients is associated with the severity of COVID-19
SARS-CoV-2 infected CVD patients were classified by their ISARIC WHO 4C-Mortality-Score on the day of admission into three subgroups of an expected mild (4C-Score <4), moderate (4C-Score 4-8), or severe (4C-Score ≥9) course of COVID-19. To determine an association between the immune signature and the severity of COVID-19, we analyzed the abundance and alterations of immune cell populations in the peripheral blood of CVD patients with expected mild, moderate, or severe course of COVID-19. The extent to which the immunophenotype within these three groups may differ was evaluated by spectral flow cytometry as described above.
Comparing the frequencies of cells in the above-defined 87 MC from patients with mild and severe or moderate and severe courses of COVID-19 revealed 19 significantly different MCs among the groups (Figure 3A, Suppl. Figure 9). These comprised reduced frequencies of neutrophils (MC8, MC16), MAIT cells (MC65, MC69)36, CD14+ monocytes expressing high levels of HLA-DR (MC12, MC28), and IgG+CD95+ B cells (MC25) in severe compared to mild SARS-CoV-2 infected CVD patients. In contrast, a greater proportion of mature NK cells (MC26, MC32), CD45RA+CD14+ monocytes (MC50), CD8+NKT cells (MC55), and several CD4+ and CD8+ T cell subsets, including CD4-CD8-HLA-DR+ T cells (MC49), naïve CD4+ T cells expressing CCR6 and CXCR3 (MC85), central memory (MC79), and effector subsets (MC51, MC87, MC74; Figure 3). Interestingly, the percentage of CD14+HLA-DR+ monocytes expressing high levels of HLA-DR (MC12, MC28), as well as that of CD8+ Tem (MC51) cells, was reduced in severe compared to moderate SARS-CoV-2 infected CVD patients (Figure 3). In addition to changes in the percentage of immune cell populations in the blood of SARS-CoV-2 infected CVD patients with different severity, we also found differences in the expression pattern of immune markers depending on disease severity (Suppl. Figure 10). For example, CD19 was significantly less expressed on MC that represent naïve and memory B cells (MC1, MC2, MC3, MC4, MC19, MC20, MC25, MC21) in SARS-CoV-2 infected CVD patients with a severe compared to mild and/or moderate course, indicating less B cell receptor signaling37. Similarly, the expression of the FcgRIII CD16 was reduced on neutrophils (MC7), CD14+ (MC12), CD14+CD16+ (MC11, MC18), CD16+ monocytes (MC37), mature NK cells (MC6, MC33), NKT cells (MC41) as well as on CD8+ Tem and Temra cells (MC51, MC52, MC40, MC46) and CD4+ Tcm cells (MC79) from CVD patients suffering from severe COVID-19 compared to those with a mild course, indicating less antibody-dependent cellular toxicity and thereby less killing of virus-infected cells (ADCC 38). Accordingly, CD14+CD16+ monocytes from severe SARS-CoV2-infected CVD patients express less HLA-DR on their surface (MC38). The chemokine receptor CXCR3, which is generally upregulated on activated NK and T cells39, had lower expression levels on immune cells from CVD patients with a severe course of COVID-19 namely on neutrophils (MC8), mature NK cells (MC32; MC33), MAIT cells (MC65) and several clusters representing CD4+ naïve, Tcm and Tem cells (MC85, MC93, MC73, MC79, MC95, MC86). CD127 (IL7Rα) was significantly less expressed by ILCs (MC62) and CD4+CD8+ T cells (MC76) from CVD patients with severe SARS-CoV-2 infection as compared to those with mild infection (Suppl. Figure 10). In contrast, the latter showed more CCR7 expression (MC76), indicating an altered function of these cells. Furthermore, CD16+ monocytes (MC23), naïve CD4+ and CD8+ T cells (MC84, MC93, MC67) from CVD patients with a severe course of infection, expressed higher levels of CD95 and CD161 (MC84, MC70), combined with lower levels of CD28 expression at least on naïve CD4+ T cells (MC85). CD16+ monocytes (MC23), naïve CD4+ T cells (MC92), and CD4-CD8-HLA-DR+ T cells (MC49) expressed higher levels of CD141 (Suppl. Figure 10).
In summary, severe SARS-CoV-2 infection resulted in fewer innate immune cells. Interestingly, although an increased proportion of NK cells, CD4+, and CD8+ T cell subsets were observed in CVD patients with severe COVID-19, their expression of functional markers like CD16 and CXCR3 was impaired, indicating an altered immune response in this subgroup.
Immune signature is predictive of severity and the course of SARS-CoV-2 infection in patients with pre-existing cardiovascular disease
Next, we looked for an immune signature that has the potential to predict the course of infection at hospital admission. We decided to use a combination of cell populations obtained through unsupervised data analysis that may provide a clear indication of disease progression through a single analysis method. To do this, we first compared the number of cells within different immune cell subpopulations in CVD patients with mild and severe courses as well as moderate and severe courses of SARS-CoV-2 infection (Figure 4A). Higher numbers of neutrophils, MAIT cells, and IgG+CD95+ B cells were found in mild-infected SARS-CoV-2 CVD patients. In contrast, the numbers of NKT cells, plasmablasts, CD4-CD8-HLA-DR+ T cells, and CD4+ Temra cells were higher in severely infected SARS-CoV-2 CVD patients (Figure 4A). When comparing moderately and severely infected patients, only the intermediate effector CD8+ T cells, CD8+ NKT cells, and cDC2 were higher in moderately SARS-CoV-2-infected CVD patients (Figure 4A&B).
After screening for the most relevant differentially abundant immune cell populations in the different patient groups, we used these to search for biomarker combinations to identify subjects at high risk of developing severe COVID-19. This search was performed using the recently published CombiROC package 20. Since moderately infected CVD patients may also develop a potentially dramatic course of infection and long-COVID, we also looked for a combination of cell populations that could distinguish this group from the cohort with a mild course. Thus, we used three settings for combined biomarker search: comparisons of (1) mild versus severe (Figure 4C), (2) mild versus moderately (Suppl. Figure 11A), and (3) mild versus moderately and severely (Figure 4C) infected CVD patients. Using CombiROC for the first comparison revealed that a combination of only three cell populations, namely MAIT cells, intermediate effector CD8+ T cells and, NKT cells predicts a severe course of SARS-CoV-2-infected CVD patients on hospital admission with a sensitivity and specificity of 1 (pink line) (Figure 4C & Table S4). The two other comparisons, (2) and (3), did not give quite as good results. In both cases, an area under the curve (AUC) of 0.98 was achieved at best, with a specificity of 1 and a sensitivity of 0.95 and 0.9267, respectively (Figure 4C & Tables S5-S6). However, a combination of 7 cell populations was necessary in both cases. The distinction between a moderate and a mild course of infection should be possible considering CD4+ T, CD8+ and CD8- NKT cells, CD8+HLA-DR+CD38+ T cells, IgG+CD95+ B cells, neutrophils, and plasmablasts. In the case of predicting a mild compared to a moderate or severe infection course if SARS-CoV-2 infection, a combination of CD4+ Temra, CD8+HLA-DR+CD38+ T cells, IgG+CD95+ B cells, neutrophils, plasmablasts, cDC1, and NKT cells is necessary.
Correlation analyses were performed between the patient groups’ most relevant differentially abundant immune cell populations (see Figure 4A) and the corresponding cytokine/chemokine plasma levels. Uninfected CVD patients were characterized mainly by negative correlations, e.g., cDC1s with IFN-a, IL-10, IL-12p70, TNF or CD8+HLA-DR+CD38+ T cells with CXCL10 (Suppl. Figure 11B). In contrast, SARS-CoV-2-infected CVD patients showed negative correlations of CD8+HLA-DR+CD38+ T cells with CCL5 or cDC2s with CXCL10, IFN-γ, and IL-10 as well as positive correlations like naïve B cells and plasmablasts with CXCL9, IL-10, and IL-18, or IFN-α with IgG+CD95+ B cells, naïve B cells, and CD8+ NKT cells (Suppl. Figure 11B). Mildly SARS-CoV-2-infected CVD patients showed positive CD8+ and CD8- NKT cells correlations with CCL2, CXCL10, CXCL9, IFN-α, IFN-γ, and IL-10, but negative correlations of cDC1 and cDC2 with CXCL10 and IFN-α (Suppl. Figure 11C). Significantly fewer and also less pronounced correlations were observed in patients with a moderate course of SARS-CoV2 infection. The positive correlations of IgG+CD95+ B cells and plasmablasts with IL-10, IL-12p70, IL-18, and, in the latter case, also with IL-6 should be noted. More correlations were generally observed again in the case of a severe SARS-CoV2 infection. Here, positive correlations of naïve B cells, plasmablasts, CD8+ NKT cells, and neutrophils with CCL2, IFN-α, IFN-γ, IL-17, and IL-6 dominated(Suppl. Figure 11C). Moreover, score plots of orthogonal partial least-squares discriminant analysis (OPLS-DA) were pereformed using immune cell populations and cytokine/chemokine plasma leverls from SARS-CoV-2 infected CVD patients (Figure 4D). Good between group variance of patients with mild (blue) compared to severe infection (red) points out a major impact of altered immune signatures on the clinical course of the disease.
In summary, using high-resolution flow cytometry, we identified an immune cell combination consisting of MAIT cells, intermediate effector CD8+ T cells, and NKT cells, stratifying CVD patients at high risk for severe SARS-CoV2 infection on the day of hospital admission.
Discussion
Our study identified a characteristic immune signature of patients with CVD and acute, symptomatic SARS-CoV-2 infection. In contrast to previous investigations, we focussed on patients with pre-existing CVD, which are at high risk of developing progressive cardiopulmonary failure due to the infection 13–15. To our knowledge, this is the first analysis to associate a specific immune signature of CVD patients with SARS-CoV-2 infection with the expected severity of COVID-19 during the disease. Here, we identified key immune system components critical to predicting a severe course of COVID-19 in our patient cohort through in-depth analyses of innate and adaptive immune cells 9, 10. Most other studies that have examined immune cell populations in the blood in SARS-CoV-2 infection have included an all-comers cohort of patients and have not stratified their analyses by concomitant comorbidities, such as diabetes, obesity, asthma, particularly the presence of CVD 4, 5, 8, 40–45. Furthermore, these studies focused on patients who either had SARS-CoV2 infection with a mild course or were hospitalized with a moderate-to-severe course. In addition, the previous analysis did not include the expected severity of COVID-19 through associations with established risk scores such as the ISARIC WHO 4C mortality score, another strength of our proposed study.
CVD is characterized by an augmented systemic inflammation triggering atherogenesis and - progression and is regulated by platelet and monocyte activation and the secretion of pro-inflammatory mediators, which play a crucial role in cell adhesion and migration 23. Consistent with and in addition to previous studies, we observed changes in the frequencies of innate immune cells like CD14+ monocytes, neutrophils, mature NK cells, and MAIT cells, as well as the adaptive immune cells CD4+CD8+, CD8+ naïve and central memory T cells, γδ T cells, early-like effector memory CD4+ T cells, and naïve B cell subsets in CVD patients compared to HD 21, 25, 45–49. Moreover, IL-6 levels were increased in CVD patients compared to HD. Thus, a chronic pro-inflammatory response in CVD patients might influence their immune signature compared to HD, leading to a more pronounced immune response to viral infections such as SARS-CoV-2.
We show that SARS-CoV-2-infected CVD patients show significant differences in immune cell populations’ distribution and phenotype compared to CVD patients without infections. Our data are consistent with other studies regarding, e.g., fewer non-classical monocytes in peripheral blood. These cells were shown to migrate to the site of organ injury, where they maintain the inflammatory processes within cardiac or pulmonary tissue 40, 42, 44, 46. As previously reported, we detected a more significant proportion of mature NK cells, CD14+CD16+ intermediate monocytes, activated populations of CD14+CD45RA+ monocytes, CD4+HLA-DR+CD38+ and CD8+ early Tem cells expressing HLA-DR+CD38+PD-1+ cells, as well as plasmablasts upon SARS-CoV-2 infection4, 5, 8, 40, 42, 43. In contrast, reduced numbers of ILCs, dendritic cell subsets cDC1 and cDC2, less activated/mature CD14+ monocytes, and CD8+ T cell subsets were observed in SARS-CoV-2 infected compared to uninfected CVD patients4, 5, 8, 41, 42, 47.
Furthermore, we describe phenotypical changes in SARS-CoV-2 infected CVD patients compared to uninfected CVD patients characterized by an overall higher expression of CD38 on plasmablasts, monocyte, NK cell, and T cell subsets 40, 42, 43 and lower expression of HLA-DR on DCs, CD14+ and CD14+CD16+ monocytes 4, 40, 42, 46, 48. Moreover, as described for NK cells 4, less CD16 was expressed on neutrophils, monocyte subsets, NKT cells, CD8+, and CD4+ T cell subsets from CVD patients with severe COVID-19 compared to those with a mild infection. These data point towards reduced antibody-dependent cellular toxicity, thereby less killing of virus-infected cells (ADCC 38). In contrast to Georg et al., we did not detect CD16+ activated T cells in severely infected SARS-CoV-2 CVD patients 49. This may be due to differences in the cohorts (all-comers versus CVD patients) and medication before blood collection (Partial medication such as antibiotics or steroids compared to no medication in our study).
Consistent with another study, we found increased expression of the thrombin receptor CD141 on monocyte clusters from SARS-CoV-2 infected CVD patients 42. CD141 reduces blood coagulation which, among other mechanisms, may be an additional explanation for thromboembolic complications, especially in CVD patients and COVID-19 32, 42, 50.
Plasmablasts and B cell subsets revealed a more inhibitory phenotype evidenced by reduced CD19 cell surface expression in SARS-CoV-2-infected CVD patients, whereas the latter was severity-dependent, indicating less B-cell receptor signaling8, 37.
In addition to other studies, we found that the remaining pDCs, regulating the type I interferon response 31, also showed a significantly more activated phenotype along with basophils and ILCs in SARS-CoV-2-infected CVD patients. pDCs and cDC subsets from SARS-CoV-2 infected CVD patients revealed a more migratory phenotype evidenced by their increased expression CCR7. In contrast, the chemokine receptor CXCR3, which is generally upregulated on activated NK and T cells39, was less expressed on neutrophils, mature NK cells, MAIT cells, and several CD4+ T cell subsets, implicating an impaired migratory capacity of these immune cells. In addition, the ligands for CXCR3, namely CXCL10, and CXCL11, were highly increased upon SARS-CoV-2 infection but also in several other viral infections and correlated with disease severity51. Thus one explanation for the reduced CXCR3 expression on many immune cells is downregulation upon ligand binding 52. Moreover, naïve CD4+ and CD8+ T cells from severely-infected CVD patients expressed higher levels of CD95 than those from mild-infected ones, indicating an impaired activation due to CD95-mediated inhibition of T-cell receptor signalling53. In summary, our data show an infection-induced activation of innate immune cells combined with reduced numbers and inhibitory phenotypes of DCs and adaptive immune cells in infected compared to uninfected CVD patients.
CVD patients with pre-existing cardiac and vascular dysfunction may benefit from intensified heart failure and anti-thrombotic therapy before progression to respiratory failure4. Since there is currently limited causal therapy for SARS-CoV-2 infection, early identification and treatment of prognostically relevant comorbidities are critical to prevent a fatal course triggered by the disease. We performed CombiROC analyses to predict the expected severity of COVID-19 at first patient contact based on a specific immune signature analyzed by flow cytometry. These analyses showed that a combination of only three immune cells detected in the blood at hospital admission (NKT cells, MAIT cells, and intermediate CD8+ effector T cells) is sufficient to distinguish CVD patients with mild from those with a severe course of infection.
Including characteristic findings in high-risk CVD patients, such as the specific immune signature, may improve the performance and generalizability of the 4C Mortality Score. This specific immune signature was independently associated with the severity of COVID-19 as determined by the ISARIC WHO 4C mortality score. In addition, there were no confounding factors in the OPLS-DA analysis, such as the comorbidities described above, suggesting that the immune signature may serve as an additional and objective tool for more intensive risk assessment of CVD patients with COVID-19 and increased risk of severe disease progression or infection-related disadvantages or complications, respectively. The ISARIC-WHO-4C mortality score is an easy-to-use, well-validated risk calculator for stratifying patients 54. However, the data collection of typical cardiovascular comorbidities such as hypertension, previous coronary artery disease, myocardial infarction, and stroke was not included in detail 55. Therefore, the risk assessment should be further improved for SARS-CoV-2 positive patients with pre-existing CVD, especially for very early stages of the viral infection. Here, we show that a pre-specified, objective analysis of the immune signature can identify CVD patients with an altered immune response at risk for severe courses of viral infections upon first admission or even at asymptomatic stages. These high-risk CVD patients could benefit from intensified monitoring and early anti-inflammatory treatment strategies, thereby preventing long-term ICU treatment or fatal outcomes. Another advantage of the immune signature is good discrimination between mild and severe courses of COVID-19 by an objective laboratory method that is not dependent on the sometimes hard-to-assess, sometimes poorly-defined, clinical parameters associated with unfavorable outcomes.
Nevertheless, the study has some limitations. Although we have studied a relatively large cohort of CVD patients with and without SARS-CoV-2 infection and HD, the number of patients with severe courses is limited. We analyzed circulating immune cells and their phenotype, more precise PBMCs. Therefore, we lose a large part of the neutrophil granulocytes in the analysis, shown to be crucial for the distinction between the fatal and non-fatal outcomes of the disease 40. And while knowledge of peripheral immune cells is essential for understanding pathogenic and protective immune responses to SARS-CoV-2 infection, it does not cover the immune response at the infection site.
In conclusion, our study revealed that a specific immune signature is associated with the severity of COVID-19 in patients with CVD and can predict the course of the disease already upon first admission to the hospital. The early determination of the immune signature might enable the treating clinicians to provide the best possible pharmacological and device-based care upon hospitalization and improve clinical outcomes.
Data Availability
All data produced in the present study are available upon reasonable request to the authors.
Sources of Funding
This work was supported by the German Research Foundation (DFG) – Project number 374031971–TRR 240 and by the Ministry of Science, Research and the Arts of the State of Baden-Württemberg (Long-COVID Funding). The funder had no role in the study design, data collection, analysis, interpretation, or manuscript writing. The corresponding author had full access to all data in the study and had final responsibility for the decision to submit for publication.
Author contributions
MG performed experimental analysis, data acquisition, and preparation of Figures. CP and SG performed the Legendplex analysis. MS and TH analyzed data using R. KALM, and MPG were involved in the conception and design of the study. KALM was responsible for the patient cohort and their classification and manuscript writing. SEA was responsible for the concept and study design, data analysis, data interpretation, drafting of the manuscript, preparation of Figures, and revising it critically for important intellectual content.
Conflict of interest
None of the authors has any conflict of interest to declare. Our study complies with the Declaration of Helsinki, our locally appointed ethics committee has approved the research protocol, and informed consent has been obtained from all participants.
Data Availability Statement
For original data, please contact Stella E. Autenrieth at stella.autenrieth{at}dkfz.de.
Acknowledgment
We thank the Core Facility Flow Cytometry of the Medical Faculty Tübingen for technical help.