Abstract
Objective While numerous studies have already compared the immune responses against SARS-CoV-2 in severely and mild-to-moderately ill COVID-19 patients, longitudinal trajectories are still scarce. We therefore set out to analyze serial blood samples from mild-to-moderately ill patients in order to define the immune landscapes for differently progressed disease stages.
Methods Twenty-two COVID-19 patients were subjected to consecutive venipuncture within seven days after diagnosis or admittance to hospital. Flow cytometry was performed to analyze peripheral blood immune cell compositions and their activation as were plasma levels of cytokines and SARS-CoV-2 specific immunoglobulins. Healthy donors served as controls.
Results Integrating the kinetics of plasmablasts and SARS-CoV-2 specific antibodies allowed for the definition of three disease stages of early COVID-19. The incubation phase was characterized by a sharp increase in pro-inflammatory monocytes and terminally differentiated cytotoxic T cells. The latter correlated significantly with elevated concentrations of IP-10. Early acute infection featured a peak in PD-1+ cytotoxic T cells, plasmablasts and increasing titers of virus specific antibodies. During late acute infection, immature neutrophils were enriched whereas all other parameters returned to baseline.
Conclusion Our findings will help to define landmarks that are indispensable for the refinement of new anti-viral and anti-inflammatory therapeutics, and may also inform clinicians to optimize treatment and prevent fatal outcome.
1 Introduction
The severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) is a highly infectious and rapidly transmittable β-coronavirus that led to the global pandemic of coronavirus disease 2019 (COVID-19). As of summer 2021, more than 180 million people were infected with almost 4 million casualties worldwide (1,2). Of great concern are new virus variants that emerge periodically and exhibit augmented pandemic potentials, can cause infection of individuals vaccinated against the original strain and re-infect patients previously recovered from COVID-19 (3,4).
COVID-19 is characterized by a diverse collection of isolated or combined symptoms ranging from mild to severe to life-threatening. While mild disease may proceed asymptomatic or show mere signs of a common cold, moderate and severe COVID-19 present with symptoms that range from fever, malaise and fatigue to neurological, dermatological, gastrointestinal and pulmonary manifestations (5). Pneumonia may lead to an acute respiratory distress syndrome (ARDS), necessitating oxygen supplementation and mechanical ventilation (6,7). Immense effort has been put into the pathophysiological and immunological profiling of the host reaction to the disease (8,9). Thus, a dysregulated host immune response was identified that causes a life-threatening cytokine storm and subsequent immune paralysis, leading to multiple organ failure and death. And even though mild-to-moderate COVID-19 does not seem to be associated with a fulminant immune activation, more and more data point to long lasting infection related sequelae including persisting nausea, fatigue and loss of smell and taste in up to 30% of patients who recovered from COVID-19 (10–12). Hence, even non-severe COVID-19 may pose a considerable hazard for global health and economy.
Although numerous studies have already addressed the differences between pathological profiles of severe and mild-to-moderate COVID-19, longitudinal trajectories within the circulating immune landscapes are still obscure (13,14). However, a rapid identification of distinct disease stages might in the future allow for refined strategies with immunomodulatory therapies that support the host’s immune response. These therapies may indeed gain importance as more and more variants may emerge that – despite efficient vaccination strategies – require the containment of viral loads after infection.
In order to rule out any bias from immunosuppressive therapy in severe COVID-19, we here concentrated on mild-to-moderate cases and analyzed the peripheral blood immune cell composition and their activation, plasma levels of cytokines and SARS-CoV-2 specific immunoglobulins. We then examined the complementary data sets for transient alterations in consecutive samples, indicative for differently progressed disease. Finally, we suggested distinct disease phases and characterized their immunological and serological layout.
2 Results
2.1 Mild-to-moderate COVID-19 induces minor changes in immune cell subpopulation distribution and activation
Twenty-eight participants were recruited from the local COVID-19 test center. Reasons for visiting were either return from high incidence regions, direct contact with SARS-CoV-2 infected individuals or experience of COVID-19 related symptoms, respectively. Eight of these participants were tested positive for SARS-CoV-2 however, did not require hospitalization due to mild symptoms. Fourteen inpatients were recruited from the division of tropical medicine and infectious diseases. All of them ran non-severe disease courses without the need for intensive care. Among all participants, two patients were available for three consecutive venipunctures, nine for two and eleven for a single one.
Table 1 lists the demographic data of patients and healthy donors showing an even distribution of sex and comparable age ranges.
In order to further our understanding of early immune cell responses in the peripheral blood during mild-to-moderate COVID-19, we performed 24-dimensional flow cytometry analyses and in a first approach compared healthy controls with SARS-CoV-2 infected patients from all time points. Our primary goal was to define whether subpopulations of neutrophilic granulocytes, monocytes and lymphocytes of the B or T lineage were specifically enriched during early infection. For that, we performed dimension reductions on our multiparametric data sets by using the embedding algorithm UMAP. Figure 1A shows the topological distributions of immune cell subpopulations based on different surface antigen expression patterns for healthy controls and COVID-19 patients, respectively. While the data imply ample variation for the abundance of CD19+ and CD177+ populations between both groups, the differences regarding CD3+ and CD14+ populations were less prominent.
(A) UMAP projection of flow cytometry data show the topological distribution of immune cell populations based on differentially expressed surface antigen patterns between healthy controls and COVID-19 patients. (B) Heatmaps of normalized protein phosphorylation data from sorted CD14+ monocytes (left) and CD19+ B lymphocytes (right) demonstrated contingent hierarchical clustering of healthy controls and COVID-19 patients.
We next sought to investigate whether SARS-CoV-2 infection induced the activation of particular signaling pathways in immune cells. For that, we sorted CD3+, CD11b+, CD14+ and CD19+ cells from peripheral blood samples and analyzed from lysates the extent of protein phosphorylation. Figure 1B exemplifies for CD14+ monocytes a hierarchical clustering that suggests increases of phosphorylated p38, ERK1/2, JNK, STAT3, STAT5 and CREB in COVID-19 patients. However, when analyzing individual signaling pathways in more detail, variances among groups turned out to be very high and obscured potential differences between patients and controls (Supplemental Table 1). Likewise, cluster analyses of phosphorylated signaling molecules in CD19+, CD3+ or CD11b+ cells did not reveal any COVID-19 specific activation of the respective pathways (Figure 1B & Supplemental Tables 2-4).
Our broad survey of surface expression- and intracellular activation data therefore hint at rather small-scale alterations among bulk immune cell responses during mild COVID-19 and hence call for a more sophisticated data-dependent stratification.
2.2 The fraction of plasmablasts in combination with antibody titers delineate infection phases
In order to explain the variation among the COVID-19 samples, we aimed to delineate a time line of early immune events following SARS-CoV-2 infection. However, neither the date of a positive test result for SARS-CoV-2 nor hospital admission are suitable reference points for the classification of different disease phases, as COVID-19 leads to a diverse collection of symptoms that may be delayed or may fail to develop altogether even in the presence of substantial viral loads (15–17). We therefore reviewed our multidimensional flow cytometry data for advanced clues that would allow for the assignment of the various samples to specific disease stages along this time line and for subsequent comparison of these disease stages with each other and to the controls.
As was recently described for COVID-19, enrichment of a distinct population of the B lineage cells may indicate a shift from the steady state immune constitution and point towards an early adaptive response (18–20). Indeed, we found a transient increase of CD19+CD45RA+CD27+CD38bright plasmablasts in consecutive samples from COVID-19 patients (Figure 2A). Accordingly, dynamics of early antibody producing cells and class-switching may translate into a time- or disease stage-dependent accumulation of immunoglobulins (Igs) of different isotypes. A very early stage post infection should be devoid of any antibodies and thereafter, IgM, IgA and IgG specific for viral proteins may appear almost simultaneously. In fact, it has been shown that IgM titers will increase for about two weeks after emergence before they plateau, while amounts of IgA and IgG are expected to increase over a longer period of time (20,21). Given these assumptions, we analyzed SARS-CoV-2 specific Igs in plasma samples and integrated the dynamics of peripheral IgD− plasmablasts in order to stratify our data set. We thus were able to define four different disease stages: (i) healthy controls that have no antibodies against SARS-CoV-2 and feature a mean of 1.1 ± 0.2% of plasmablasts among CD19+CD45RA+ lymphocytes, (ii) incubation phase with still negligible SARS-CoV-2 specific Igs and a mean of 2.9 ± 0.8% of plasmablasts, (iii) early acute infection characterized by elevated IgM, IgG and IgA titers along with an significantly increased proportion of plasmablasts (11.5 ± 2.6%) and thereafter (iv) late acute infection phase defined by a reduction in plasmablasts (4.4 ± 1.1%) together with stable IgM and minor increases in IgA and IgG titres (Figure 2B & supplemental Tables 5-6). Of note, we assumed that age would be a confounding variable influencing some of our data and therefore incorporated analyses of covariants (ANCOVA) in order to eliminate age from the statistical model (Figure 2B & supplemental Table 5). Indeed, we found that age specifically operated on null hypothesis rejection for the alteration of CD27+CD38dim early memory B-lymphocytes which were nonetheless altered between different disease stages when precluding the confounding covariable (supplemental Table 5).
(A) Representative pseudocolor plots of CD38 and CD27 expression data of CD19+CD45RA+ B lymphocytes showed the disease stage-dependent increase of plasmablast proportions (rectangular gate). (B) Top left: Quantitative data of plasmablast abundances at different disease stage. Top right, bottom left and bottom right: Semi-quantitative data of SARS-CoV-2 specific IgM, IgG and IgG titers in plasma samples. **p < 0.01; ***p < 0.001; Tukey-Kramer test for plasmablasts, Kruskal-Wallis test with Dunn’s correction for multiple comparisons for Antibody isotypes. MFI: Median fluorescence intensity. A450: Absorbance at 450 nm.
In summary, our data demonstrated that utilizing the dynamics of both, plasmablast trajectories and SARS-CoV-2 specific Ig titers allowed for the delineation of patient data into differently progressed disease phases.
2.3 Pro-inflammatory monocytes and CD16−CD177+ granulocytes follow different kinetics during mild-to-moderate COVID-19
Building on the fact that early phases of viral infection trigger innate immune responses, we sought to investigate the dynamics of monocyte and granulocyte subpopulations between the above postulated disease phases. Peripheral monocytes were defined via sideward scatter and expression of CD14. The additional analysis of CD16 expression allowed for the differentiation of classical (CD14+CD16−), pro-inflammatory (CD14+CD16+) and anti-viral monocytes (CD14+CD16bright). Figure 3 demonstrates for the contingent of pro-inflammatory monocytes a moderate but statistically not quite significant increase in the incubation phase when compared to healthy controls (p = 0.0766). However, upon progression into the acute infection phases, this population of monocytes was significantly reduced (Figure 3). Finally, we did not find any meaningful differences between the remaining monocyte subpopulations and the various phases of the disease (supplemental Table 7). For granulocytes we gated SSChi cells and then differentiated subpopulations by their expression patterns for CD11b, CD16 and CD177. While most of the granulocyte subpopulations were unaltered between all disease stages, we found that CD16−CD177+ cells were slightly increased during the incubation phase and early acute infection. However, during late acute infection, there was a significant increase in this population compared to healthy controls (Figure 4 & supplemental Table 8).
(A) Representative pseudocolor plots of CD14 and CD16 expression data of SSCmed monocytes demonstrated kinetics for the enrichment of pro-inflammatory monocytes (polygonal gate). (B) Quantitative data of pro-inflammatory monocyte proportions. *p < 0.05; **p < 0.01; Tukey-Kramer test.
(A) Representative pseudocolor plots of CD16 and CD177 expression data of SSChi leukocytes showed that CD16−CD177+ granulocytes specifically emerged at the late acute infection phase (Q3 gate). (B) Quantitative flow cytometry data indicated gradually increasing abundances of CD16−CD177+ granulocytes. **p < 0.01; Tukey-Kramer test.
In summary, the definition of early disease phases allowed for the detection of different kinetics among innate immune cells. While the proportion of pro-inflammatory monocytes increased rather early in the response to SARS-CoV2, CD16−CD177+ granulocytes were found to accumulate in the late acute infection stage.
2.4 Incubation and early acute infection phases are characterized by an increase in activated cytotoxic T cells and IP-10
In order to characterize in detail the cytotoxic T cell response during incubation as well as early and late acute infection phases, CD8+ cells were analyzed for their activation (CD25, CD38), exhaustion (PD-1) and differentiation status (CD45RA, CD27, CD127). While the fraction of naïve (CD45RA+CD27+) CD8-positive T cells was strictly age dependent, there were no further associations with any of the disease phases. However, the incubation phase was characterized by a significant relative increase in CD45RA+CD38+CD27− cytotoxic T cells (Figure 5A), indicating an activated and terminally differentiated status (22). This fraction of CD38+CD27− cytotoxic T cells gradually declined over the different phases of acute infection (Figure 5A-B). Following a different kinetic, CD45RA+CD38+CD27+ cytotoxic T cells, suggestive of TSCM (stem cell memory) started to increase during incubation yet reached a significant difference compared to healthy controls only during late acute infection (Figure 5A and supplemental Table 9). Likewise, CD45RA+CD127− CD27+ cells, indicative of TEMRA (effector memory re-expressing RA), started to increase during early acute infection and were significantly elevated during late acute infection (supplemental Table 9). Moreover, all disease phases featured an elevated expression of CD38 on CD8+ cells (Figure 5C & supplemental Table 10). Finally, we observed a significantly increased fraction of PD-1+ cells among cytotoxic T cells during the early acute infection phase. These PD-1+ cells were either CD27– or CD27+ and already in the decline during late acute infection (Figure 5D-E & supplemental Table 9). In contrast to the alterations among CD8+ subpopulations, we did not observe similarly prominent changes among CD4+ T helper cell differentiation/activation stages and the various disease phases (supplemental Table 11). Neither were there significant changes among CD25+CD127+ regulatory T cells nor NK cells (supplemental Tables 11-12).
(A) Representative pseudocolor plots of CD38 and CD27 expression data of CD8+ T cells demonstrated the accumulation of the CD38+CD27− subpopulation (Q3 gate) during the incubation and early acute infection phases. (B) Quantitative flow cytometry data of CD8+CD38+CD27− T cell proportions. (C) CD38 expression by CD8+ T cells at the different disease stages. (D) Representative plots of PD-1 and CD27 expression data showed an increase of PD-1+ T cells that are either CD27+ (Q2 gate) or CD27− (Q3 gate). (E) Quantitative data of CD8+PD-1+ T cell proportions showed their enrichment during the incubation and acute infection phases. *p < 0.05; **p < 0.01; ***p 0.001; Tukey-Kramer test.
Finally, we investigated the kinetics of 13 cytokines in an anti-virus response panel and found that IP-10 (CXCL10) was transiently elevated during both the incubation phase and early acute infection before returning to baseline in late acute infection (Figure 6 & supplemental Table 13). In support of mild-to-moderate COVID-19 courses in our cohort, we did not detect a cytokine storm for any of the analytes measured. We did though find slight increases of IL-6, IL-8, IL-12, TNFα, IFNγ, IFNα2, IFNλ1, IFNβ, GM-CSF and IL-10, mostly for the acute infection stages. However, none of these exhibited a meaningful effect size yet were subject to high variance (supplemental Table 13).
Quantitative data of IP-10 concentrations in plasma samples at different disease stages. *p < 0.05; **p < 0.01; Tukey-Kramer test.
In summary, while the increase of activated CD38+CD27– and activated or exhausted PD-1+ cytotoxic T cells was mostly restricted to the incubation phase and early acute infection respectively, we found the activation marker CD38 to be overexpressed during all disease phases.
2.5 Incubation and early acute infection are characterized by T lymphopenia unrelated to apoptosis
In order to also assess quantitative differences among the major leukocyte populations during the early disease phases, we determined absolute numbers of live granulocytes, lymphocytes and monocytes. The picture emerging revealed lymphopenia during incubation and early acute infection which slowly recovered during late acute infection (supplemental Table 14). This lymphopenia was based on a decrease in T lymphocytes and affected both, T helper cells and cytotoxic T cells (Figure 7A). Quantitative differences among B cell, granulocyte, and monocyte bulk populations did not reach statistical differences (supplemental Table 14). To rule out apoptosis as the reason for T lymphopenia, we calculated the respective proportions of apoptotic cells among live leukocytes. However, no significant variations emerged that paralleled the fluctuating numbers of live T cells nor were there any other populations affected (Figure 7B & supplemental Table 15). Apoptosis therefore could not be held responsible for the decrease of peripheral T lymphocytes during incubation and early acute infection.
(A) Total cell numbers of cytotoxic T cells (left) and T helper cells (right). (B) Quantitative data of apoptotic cell proportions among cytotoxic T cells (left) and T helper cells (right). *p < 0.05; **p < 0.01; ***p 0.001; Tukey-Kramer test.
2.6 IP-10 in patients correlates with different immune parameters than in healthy controls
Finally, we performed correlation analyses for all immune parameters assessed. As for the healthy controls, the only significant correlation existed between the fraction of pro-inflammatory monocytes (CD14+CD16bright) and serum concentrations of IP-10 (supplemental Table 16). In contrast, in the patients, serum IP-10 correlated significantly with CD38+CD27− terminally differentiated cytotoxic T cells (supplemental Table 17). Moreover, patient serum concentrations of all anti-S isotypes correlated significantly with each other and also with the fraction of immature neutrophilic granulocytes.
3 Discussion
By analyzing longitudinal trajectories of immune responses to SARS-CoV-2 in patients experiencing mild-to-moderate disease courses, we were able to define three different stages (Figure 8). The first one comprises the incubation phase, which for the patients manifested with symptoms that led to voluntary testing or – in case of elderly or multimorbid patients – to hospitalization. The immune response in this phase of disease progression is characterized by significant increases in pro-inflammatory monocytes (CD16bright CD14+), activated cytotoxic T cells (CD38+), and serum IP-10 paralleled by a significant decrease of peripheral T helper cell numbers (Figure 7). This incubation phase is followed by the early acute infection where peripheral plasmablasts (CD27+ CD38bright) are transiently yet significantly elevated as are exhausted cytotoxic T cells (PD-1+). IP-10 is still significantly upregulated while both, helper and cytotoxic T cells are significantly reduced. Virus-specific IgM, IgG and IgA antibodies are detectable and except for IgM still on the rise. Finally, during late acute infection, most values appear on their return to normal except for virus-specific antibodies and immature neutrophils which are now significantly elevated compared to the healthy controls (Figure 8).
We here show that the transient increase in peripheral plasmablasts in combination with seroconversion and relative increases in Ig subclasses serve as reference points to time early phases of COVID-19. While plasmablasts increase between days 1 and 9 after onset of symptoms and peak around day 9 (21), SARS-CoV-2 specific antibodies can be detected as early as 3-6 days post the onset of symptoms (19,20). Importantly, Ig subclasses in the serum did not appear to follow the classical IgM – IgG – IgA sequence of events as shown e.g. for the seroconversion following HIV infection (23,24). Instead, whether IgM or IgG appeared first seemed to depend on the antigen with anti-N and anti-S2 inducing a robust IgG and anti-S1 and anti-RBD inducing an IgM response first. IgA appeared simultaneously with either IgM or IgG however, while IgM on average plateaued during early acute infection - IgG and IgA were still on the rise during late acute infection. We here did not specifically assess SARS-CoV-2 neutralizing antibodies yet refer to a previous publication demonstrating that IgA dominates the early neutralizing response (21).
Taking the alteration of peripheral plasmablasts and virus-specific antibodies as a basis to define various phases of COVID-19 allowed for the detection of small scale alterations among bulk immune cell responses which otherwise were obscured. E.g. we here show that the incubation phase is characterized by a sharp increase in pro-inflammatory monocytes and indeed, this increase has previously been confirmed for acute Dengue and HI virus infections, alike (25,26). During Dengue virus infections, CD16brightCD14+ pro-inflammatory monocytes were suggested to take on an important role as they were shown to migrate to the draining lymph nodes where they encounter and activate antigen-specific memory B cells to differentiate into plasmablast that secrete IgG and IgM (25). Moreover, correlations with serum concentrations of IP-10 hinted at pro-inflammatory monocytes as the main producers of this chemokine (27). We can confirm this correlation for the healthy controls only. However, in the patients, IP-10 correlated significantly with terminally differentiated cytotoxic T cells. We therefore postulate that IP-10, which is held responsible for the cytokine storm in critically ill patients, is produced by cytotoxic T cells (28–30).
Immature neutrophils, characterized by a lack of CD16 expression, in our study continuously increased over time yet reached a significant difference compared to healthy controls only during late acute infection. Interestingly, a handful of longitudinal studies on hospitalized COVID-19 patients have shown that severity associated with an increase in neutrophil numbers and that these numbers were persistently elevated, even at clinical recovery (31). Even though we here did not assess severe cases, our findings also suggest an unresolved inflammation at a time when most immune parameters returned to normal.
What was conspicuously absent in our study was any major alteration among CD4+ T helper cell subpopulations including Tregs. In contrast, there were significant changes among cytotoxic subpopulations, among them an elevated CD38 expression, suggesting activation over all disease phases analyzed. Moreover, we observed a significant increase in terminally differentiated (CD27−) cytotoxic T cells during incubation which gradually declined during later phases. In chronic HIV infection this population was shown to be inefficient and to correlate with viral burden (32,33). Clearly, its role in COVID-19 awaits further elucidation. An opposing trend – increase during acute infection – was observed for both, the CD45RA+CD38+CD27+ and CD45RA+CD127−CD27+ populations, reminiscent of gradual differentiation into stem cell memory (TSCM) and effector memory (TEMRA), respectively (22,34). Finally there also was a significantly increased fraction of PD-1+ cytotoxic T cells during early acute infection. Even though PD-1 expression on SARS-CoV-2 specific CD8+ cells was previously associated with an activated phenotype (35), prolonged antigen stimulation and sustained expression of inhibitory receptors like PD-1 can ultimately lead to impaired T cell functionality and exhaustion (36).
As for significant alterations among the bulk populations of B cells, T cells, NK cells, monocytes, during the early phases of COVID-19, we only observed T-lymphopenia and this was more pronounced for the helper than the cytotoxic T cells. However, this drop in peripheral T cells was not paralleled by an increase in apoptotic cells suggesting that, rather than being eliminated, T cells possibly relocated to the sites of infection or to lymphatic tissues. Of note, we did not observe any significant alteration among NK cells, neither among the bulk nor among subpopulations. Likewise, there were no alterations among monocyte/NK or monocyte/granulocyte ratios which confirms previous findings about mild-to-moderate COVID-19 (31,37).
In summary, the present study defines a plasmablast-antibody landscape that defines different disease stages during mild-to-moderate COVID-19. These disease stages were further characterized by means of cellular parameters and identify the time frame where cytokinogenesis was most prominent. Upcoming research needs to focus on comparable disease stages during severe COVID-19 in order to allow for the identification of landmarks where resolution of viral infection segregates from progressing immune dysregulation and fatal outcomes. Moreover, as global efforts are underway to identify and develop new anti-viral and immunomodulatory therapeutics that reduce COVID-19 related hospitalization and deaths, need arises for the right timing (38). Indeed, early type I IFN treatment seemed to be critical for therapeutic success (39). However, prophylactic immunomodulation might not be feasible for most COVID-19 patients due to high costs and potential adverse events. Therefore, identifying immunological landmarks will enhance a therapeutic potential towards benefits rather than harm.
4 Methods
4.1 Study Population
COVID-19 patients were recruited either from the local COVID-19 test center or from the division of tropical medicine and infectious diseases. After an initial withdrawal of peripheral blood on the day of testing or admittance to the hospital (day 0), subsequent venipunctures followed on days 2 and 6. Individuals recruited at the test center who were tested negative, were considered healthy controls and they provided blood only once at the day of testing. The study was approved by the ethics committee of the University Medical Center of Rostock. It is filed under A 2020-0086 and written informed consent was obtained from all participants.
4.2 Flow Cytometry
For the analysis of surface expression markers, 100 µL of anticoagulated COVID-19 patient or healthy donor blood were used. In order to reduce unspecific antibody-conjugate binding, 10 µL fetal calf serum (Thermo Fisher, Waltham/MA, United States), 5 µL True-Stain Monocyte BlockerTM and 5 µL anti-Fc receptor TruStain FcXTM (Biolegend, San Diego/CA, United States) were added and incubated for 15 min on ice. The following amounts of antibody:fluorophore-combinations were used: 0.25 µg CD127:APC/R700 (clone HIL-7R-M21), 1 µg CD147:BV421 (TRA-1-85), 0.5 µg CD45RO:BV480 (UCHL1, BD Biosciences, Franklin Lakes/NJ, United States), 1 µg CD11b:PerCP/Cy5.5 (ICRF44), 0.8 µg CD11c:BV785 (3.9), 0.56 µg CD14:BV510 (63D3), 0.13 µg CD16:BV650 (3G8), 0.06 µg CD19:APC/Fire810 (HIB19), 0.13 CD20:SparkNIR685 (2H7), 0.5 µg CD27:BV605 (O323), 0.25 µg CD3:SparkBlue550 (SK7), 0.25 µg CD304:AlexFluor647 (12C2), 0.03 µg CD4:BV750 (SK3), 0.5 µg CD45RA:APC/Fire750 (HI100), 0.13 µg CD56:BV711 (5.1.H11), 0.13 µg CD8:BV570 (RPA-T8), 0.5 µg CD95:PE/Cy5 (DX2), 0.13 IgD:PE/Dazzle594 (IA6-2), 0.13 µg PD-1:APC (A17188B, Biolegend), 0.06 µg CD38:PerCP/eFluor710 (HB7, Thermo Fisher), 0.06 µg CD177:PE/Vio770 (REA258), 0.05 µg CD25:PE (REA570, Miltenyi, Bergisch Gladbach, Germany). Antibodies were incubated for 45 min on ice in the dark. In parallel, ApotrackerTM Green (Biolegend) was added according to the manufacturer’s instruction with an incubation for 30 min on ice in the dark. In order to lyse erythrocytes, 2.2 mL Fixative-Free Lysing Solution (Thermo Fisher) was added and incubated for 20 min at room temperature in the dark. Following this, 0.03 µg 4’,6-diamidino-2-phenylindole was added as a live/dead discriminator and incubated for 5 min. Finally, data acquisition was performed on the Cytek® Aurora flow cytometer running on the SpectroFlo Software version 2.2.0.3 (Cytek Biosciences, Fremont/CA, United States). Analysis of flow cytometry data was done using the FlowJo software version 10.7 (FlowJo, Ashland/OR, United States) with the manual gating strategy shown in supplemental Figure 1. Dimension reduction of down-sampled and concatenated data sets was performed using the FlowJo plugin for the algorithm “uniform manifold approximation and projection” (UMAP) (40).
4.3 Multipathway phosphoprotein analyses
T lymphocytes were obtained by processing 2 mL anticoagulated whole blood using CD3 MicroBeads on an automated magnetic activated cell sorting device (autoMACS, Miltenyi) according to the manufacturer’s instructions. The fraction enriched for CD3+ cells was centrifuged, washed with PBS and frozen at -80°C. Cells from the CD3− fraction were washed with Running Buffer (RB, Miltenyi) and erythrocytes were lysed by adding 5 mL Fixative-Free Lysing Solution (Thermo Fisher), mixing and incubating for 15 min at room temperature. The lysis was quenched by adding 10 mL RB and mixing. Subsequently, cells were pelleted, suspended in RB and B lymphocytes were isolated by utilizing CD19 MicroBeads (Miltenyi) on an autoMACS. The fraction enriched for B lymphocytes was then washed and the cell pellet was frozen thereafter as described above. For the isolation of monocytes we used 2 mL anticoagulated blood and CD14 MicroBeads (Miltenyi). The CD14+ fraction was processed as mentioned above. The CD14− fraction was submitted to the isolation of granulocytes by CD11b MicroBeads (Miltenyi) similar to the enrichment of CD19+ cells.
Proteins were obtained from sorted CD3+, CD19+, CD14+ and CD11b+ cells by adding 40 µL MILLIPLEX® Map Lysis Buffer containing 100-fold diluted Protease Inhibitor Cocktail Set III (Merck, Darmstadt, Germany) and mixing at 4°C and 700 rpm for at least 15 min. Samples were subsequently centrifuged at 10,000 x g for 10 min at 4°C. Supernatants were diluted fourfold and bulk protein concentrations were determined by using the PierceTM BCA Protein Assay Kit (Thermo Fisher) following the manufacturer’s guidelines. The absorbance was measured at 562 nm on the Infinite® 200 automated plate reader (Tecan, Männedorf, Switzerland). Subsequently, samples were processed using the MILLIPLEX® Multi-Pathway Cell Signaling Assay kits (Merck) for the semiquantitative analysis of either phosphorylated or unphosphorylated cAMP response element-binding protein (CREB), extracellular signal-regulated kinases (ERK)1/2, nuclear factor kappa-light-chain-enhancer of activated B cells (NF-κB), c-Jun N-terminal kinase (JNK), p38, p70S6 kinase, signal transducer and activator of transcription (STAT)3 and STAT5 to the manufacturer’s instructions. In brief, 1-25 µg of total protein was used. Samples where protein concentrations were insufficient were excluded from further analysis. For calibration purposes, unstimulated and stimulated cell lyophilisates were reconstituted to a protein concentration of 2 mg/mL. In detail, the kit provided unstimulated HeLa cells, HeLa cells stimulated with tumor necrosis factor (TNF)α/Calyculin A, MCF7 cells stimulated with insulin-like growth factor 1 and A431 cells stimulated with epidermal growth factor. Serial dilutions were performed and processed in parallel to the samples. Data acquisition was performed on the Luminex® 100/200TM flow cytometer with the acquisition software xPONENT® version 3.1.871.0 (Luminex, Austin/TX, United States).
Data analysis was performed by log10-transformation of raw median fluorescence intensity (MFI) values and annotating arbitrary units (AUs) for the respective standard dilution. In order to obtain AUs for samples, we utilized a linear regression model for the calibration standards onto which transformed sample MFI values were fitted. After back-transformation, the data for unphosphorylated proteins were normalized to the amount of whole protein from the BCA assay. AUs for phosphorylated proteins were normalized to AUs from their unphosphorylated counterparts. For heatmap visualization, the data was normalized to a 0 – 1 range (z value).
4.4 SARS-CoV-2 specific antibodies
Plasma samples were obtained by centrifugation of anticoagulated blood from healthy donors and COVID-19 patients, respectively. For the semi-quantitative analyses of IgM and IgG specific for the nucleocapsid (N) or the spike protein subunits S1, S2 and receptor binding domain (RBD), we used the respective MILLIPLEX® SARS-CoV-2 antigen panel kits (Merck) to the manufacturer’s instructions. In brief, plasma samples were diluted 100-fold in assay buffer and were incubated with antigen-conjugated beads for 2 h at room temperature with shaking at 600 rpm. After washing, PE-conjugated anti-human IgG/IgM detection antibodies were added to the samples followed by incubation for 1.5 h at room temperature and 600 rpm. Samples were then washed, immersed in sheath fluid (Merck) and analyzed on the Luminex® 100/200TM flow cytometer (Luminex). For the detection of trimeric spike protein specific IgA we used an enzyme-linked immunosorbent assay according to the manufacturer’s specifications (Euroimmun, Lübeck, Germany). Absorbances were detected at 450 nm (A450) on the Infinite® 200 automated plate reader (Tecan).
4.5 Cytokine analysis
For the determination of cytokine concentrations in plasma samples, the LEGENDplexTM Human Anti-Virus Response Panel (Biolegend) was used containing capture beads and detection antibodies for Interleukin (IL-)1β, IL-6, IL-8, IL-10, IL-12p70, Interferon (IFN)α, IFNβ, IFNλ1, IL-29, IFNλ2/3, IL-28, IFNγ, TNFα, interferon gamma-induced protein 10/C-X-C motif chemokine ligand 10 (IP10/CXCL10) and granulocyte-macrophage colony-stimulating factor (GM-CSF). The protocol was followed using the manufacturer’s instructions. Data acquisition was performed using the Cytek® Aurora (Cytek Biosciences).
4.6 Statistics
Data analysis and visualization were performed using R (version 3.5.1), InStat (GraphPad, San Diego/CA, United States) and SPSS (IBM, Armonk/NY, United States). Contingency tables were compared by the Fisher’s exact test. Normality of data sets was evaluated using the Kolmogorov-Smirnov test. Comparisons for multiple groups were done by analysis of variance (ANOVA) under the assumption for normally distributed sample data. Pairwise comparisons were conducted post-hoc with the Tukey-Kramer test. Non-normally distributed data were compared by the Kruskal-Wallis one-way analysis of variance with post-hoc pairwise comparisons using the Mann-Whitney-Wilcoxon (MWW) test and adjustment of p-values with the Holm-Bonferroni method. The Spearman rank method was used for correlation analyses. A p-value < 0.05 was considered statistically significant.
Data Availability
All data generated or analyzed during this study are included in this published article (and its supplementary information files).
5 Data Availability Statement
All data generated or analyzed during this study are included in this published article (and its supplementary information files).
6 Ethics Statement
The studies involving human participants were reviewed and approved by the ethics committee of the University Medical Center of Rostock (reference number A 2020-0086). Written informed consent was obtained from all study participants.
7 Author Contributions
MM, ECR and BM-H contributed to conception and design of the study. MM and JV acquired the data. BM-H and JV performed the statistical analyses and wrote the first draft of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version.
8 Funding
This study was financially supported by the Federal State of Mecklenburg-Western Pomerania via the Sondervermögen des MV Schutzfonds Säule Gesundheit.
11 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.
14 Supplementary Material
CD14+ Multi-Pathway Signaling
CD3+ Multi-Pathway Signaling
CD19+ Multi-Pathway Signaling
CD11b+ Multi-Pathway Signaling
Percentages of B cell Subpopulations
Serum Antibodies
Percentages of Monocyte Subpopulations
Percentages of Granulocyte Subpopulations
Percentages of CD8+ T Cell Subpopulations
Expression of CD38 among CD8+ T Cell Subpopulations
Percantages of CD4+ T Cell Subpopulations
Percantages of NK Cell Subpopulations
Serum Cytokines
Absolute Numbers of Live Leukocytes
Apoptotic among Live Leukocytes [%]
Correlations among Humoral and Cellular Immune Parameters in Healthy Controls
Correlations among Humoral and Cellular Immune Parameters in COVID-19 Patients
Hierarchical gating was performed with FlowJo Software. Representative dot plots were selected and combined from different samples to show all subpopulation mentioned in this study. Colored frames illustrate the main cell populations in peripheral blood (Granulocytes, Monocytes, B cells, CD4+ T cells and CD8+ T cells. Grey arrows indicate the hierarchical gating steps.
9 Acknowledgments
We thank all our participants for taking part in this study and providing blood samples. The authors want to particularly thank Wendy Bergmann (Core Facility for Cell Sorting and Cell Analysis, Rostock University Medical Center) for help during flow cytometry panel design and setup and Micha Löbermann (Department of Tropical Medicine and Infectious Diseases, Rostock University Medical Center) for helping recruit Participants.