ORIGINAL RESEARCH article

Front. Immunol., 24 November 2023

Sec. Microbial Immunology

Volume 14 - 2023 | https://doi.org/10.3389/fimmu.2023.1260283

High-dimensional phenotyping of the peripheral immune response in community-acquired pneumonia

  • 1. Center for Experimental and Molecular Medicine (CEMM), Amsterdam UMC location University of Amsterdam, Amsterdam, Netherlands

  • 2. Department of Molecular Cell Biology and Immunology, Amsterdam UMC location Vrije Universiteit Amsterdam, Amsterdam, Netherlands

  • 3. Department of Internal Medicine, Flevo Hospital, Almere, Netherlands

  • 4. Department of Internal Medicine, BovenIJ Hospital, Amsterdam, Netherlands

  • 5. Application Department, Cytek Biosciences, Inc., Fremont, CA, United States

  • 6. Division of Infectious Diseases, Amsterdam UMC location University of Amsterdam, Amsterdam, Netherlands

  • 7. Infectious Disease, Leiden Universitair Medisch Centrum, Leiden, Netherlands

Abstract

Background:

Community-acquired pneumonia (CAP) represents a major health burden worldwide. Dysregulation of the immune response plays an important role in adverse outcomes in patients with CAP.

Methods:

We analyzed peripheral blood mononuclear cells by 36-color spectral flow cytometry in adult patients hospitalized for CAP (n=40), matched control subjects (n=31), and patients hospitalized for COVID-19 (n=35).

Results:

We identified 86 immune cell metaclusters, 19 of which (22.1%) were differentially abundant in patients with CAP versus matched controls. The most notable differences involved classical monocyte metaclusters, which were more abundant in CAP and displayed phenotypic alterations reminiscent of immunosuppression, increased susceptibility to apoptosis, and enhanced expression of chemokine receptors. Expression profiles on classical monocytes, driven by CCR7 and CXCR5, divided patients with CAP into two clusters with a distinct inflammatory response and disease course. The peripheral immune response in patients with CAP was highly similar to that in patients with COVID-19, but increased CCR7 expression on classical monocytes was only present in CAP.

Conclusion:

CAP is associated with profound cellular changes in blood that mainly relate to classical monocytes and largely overlap with the immune response detected in COVID-19.

Introduction

Community-acquired pneumonia (CAP) is a preeminent driver of hospitalization, short-term mortality and long-term decline (). CAP is initiated by a pathogen, but symptoms and life-threatening complications are primarily caused by the collateral damage of dysregulated host defense mechanisms aimed at eradicating the causative microorganism (, , ). While in vitro and animal studies have generated extensive knowledge of the phenotype and function of immune cells involved in pneumonia pathophysiology, studies that report in-depth cellular phenotyping in patients with CAP remain scarce. In stark contrast, numerous deep immunological profiling studies have been performed in patients with CAP caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) – the pathogen responsible for coronavirus disease (COVID)-19 – relaying intricate details on immune phenotypes and functional states at single cell level (). This knowledge has contributed to various immunotherapeutic trials, and the successful application of some of these therapies in clinical practice (). No such immunotherapies are considered standard of care in patients with other forms of CAP.

We recently reported on the immune response in different forms of CAP based on single-cell RNA sequencing and an extensive set of host response plasma biomarkers (, ), and we described differences in lung and blood immunophenotypes in critically ill patients with COVID-19 (). Here, we use high-dimensional spectral flow cytometry to provide a comprehensive analysis of the peripheral blood immune system compartment at the single-cell protein level in patients with CAP, using age- and sex-matched control subjects without acute infection and patients with COVID-19 as comparator groups.

Materials and methods

Study design and participants

This project was part of the ELDER-BIOME study (ClinicalTrials.gov Identifier: NCT02928367) (, ). Patients for this project were consecutively enrolled between October 2018 and March 2020 for CAP and in April and May 2020 for COVID-19 (prior to the introduction of the Alpha variant, the availability of vaccines, and dexamethasone becoming standard of care). Inclusion criteria were as reported previously (), with additional exclusion criteria for comorbidities associated with an immunocompromised state (Supplemental Methods). Age- and sex-matched control subjects without signs of acute infection were recruited from the outpatient clinic. The study was approved by the medical ethical committee of the Amsterdam University Medical Centers. All patients or their representatives provided written informed consent.

Sampling

Blood samples were obtained from patients within 48 hours of admission to a general hospital ward and processed within 4 hours after collection. Peripheral blood mononuclear cells (PBMCs) were isolated from heparin-anticoagulated blood by density gradient centrifugation using Ficoll-Paque plus medium (GE Healthcare Life sciences, Eindhoven, The Netherlands) and cryopreserved in liquid nitrogen until further use, as described (). Plasma was obtained from EDTA-anticoagulated blood by centrifugation and stored at -80 °C until use.

Spectral flow cytometry – staining, data acquisition and analysis

We performed 36-color spectral flow cytometry based on the Optimized Multicolor Immunofluorescence Panel (OMIP) 069 (). Details on experimental procedures and antibodies used can be found in the Supplemental Methods and Supplementary Table 1. Live single CD45+ cells (Supplementary Figure 1) were clustered using FlowSOM () (version 1.18) into 600 SOM clusters. The optimal number of metaclusters (MCs) was determined as described in the Supplemental Methods and Supplementary Figure 2A. MCs were phenotyped by comparison to manual gating strategies as described (). MCs with the same immunological phenotype were given numbers, and represent differential (activation) states of these phenotypes.

Protein biomarker assays

Plasma protein biomarkers were measured using Luminex (R&D systems, Minneapolis, Minnesota, USA) or cytometric bead array (BioLegend, San Diego, California, USA) as described ().

Statistical analysis

We initially determined which MCs were differentially abundant between patient groups in an untargeted manner (i.e. the volcano plot) by performing Welch’s t-test on log2-transformed MC proportions (of total cells), and applied Benjamini-Hochberg (BH) adjustment to the resulting P-values. In subsequent analyses, we expressed the difference in log2-transformed MC proportions between patient groups and controls by Hedges’ g, a commonly used effect size measure that incorporates both the differences in means and a shared measure of variance (). An overview of the study design, including analysis levels and comparisons, is depicted in Figure 1A. Further details can be found in the Supplemental Methods.

Figure 1

Results

Patient characteristics

We enrolled 40 consecutive patients with CAP within 48 hours after hospital admission, and 31 age- and sex-matched control subjects without signs of acute infection (Table 1). Patients and controls were well-matched, although patients more often had chronic obstructive pulmonary disease. Vital signs, disease severity scores and clinical outcomes reflected an overall low to moderate severity of disease; all patients survived up to day 28.

Table 1

CAP
(n = 40)
Controls
(n = 31)
P-value
DEMOGRAPHICS
 Age, years68.8 (15.4)64.3 (15.6)0.23
 Sex, male24 (60.0)20 (64.5)0.81
COMORBIDITIES
 Chronic obstructive pulmonary disease15 (37.5)2 (6.5)0.002
 Asthma6 (15.0)1 (3.2)0.13
 Hypertension16 (40.0)15 (48.4)0.63
 Diabetes mellitus, type 24 (10.0)7 (22.6)0.19
 Chronic kidney disease3 (7.5)3 (9.7)>0.99
LABORATORY TESTS*
 Platelets, x109 cells/L273.9 (119.7)
 Leukocytes, x109 cells/L13.2 [9.4, 17.0]
 Neutrophils, x109 cells/L11.5 [7.2, 14.9]
 Lymphocytes, x109 cells/L0.80 [0.60, 1.28]
 Neutrophil-to-lymphocyte ratio9.9 [7.3, 15.4]
 C-reactive protein, mg/L150.8 [68.2, 264.5]
VITAL SIGNS AND DISEASE SEVERITY*
 Mean arterial pressure, mmHg95 (14)
 Respiratory rate, bpm22 [19, 28]
 Temperature, °C37.8 (1.3)
 Modified Early Warning Score3 [2, 4]
 Pneumonia Severity Index4 [2, 4]
 CURB-651 [1, 2]
 qSOFA1 [0, 1]
CLINICAL COURSE AND OUTCOMES
 Symptoms to admission, days3 [2, 5]
 Pathogen identified16 (40.0)
 ICU stay (at any point during admission)3 (7.5)
 Hospital length of stay, days5 [3, 9]
 Time to clinical stability†† or discharge, days5 [2, 8]
 28-day mortality0 (0)

Baseline characteristics and outcomes.

CAP, community-acquired pneumonia; bpm, breaths/beats per minute; CURB-65, confusion, blood urea nitrogen, respiratory rate, blood pressure, age 65 or older; qSOFA, quick sequential organ failure assessment score.

Normally distributed continuous data are displayed as mean (standard deviation) and compared using Welch’s t¬-test; non-normally distributed continuous data are displayed as median [interquartile range]; categorical data are displayed as count (percentage) and compared using Fisher’s exact test.

* Measured upon presentation to the emergency ward.

† Pathogens identified: influenza A (n = 4), Pseudomonas aeruginosa (n = 4), Streptococcus pneumoniae (n = 3), Moraxella catarrhalis (n = 2), Staphylococcus aureus (n = 2), Haemophilus influenzae (n = 1), rhinovirus (n = 1), respiratory syncytial virus (n = 1), and Aspergillus spp. (n = 1).

Co-infections in three patients: (1) Pseudomonas aeruginosa and Staphylococcus aureus, (2) Streptococcus pneumoniae and influenza A, and (3) Streptococcus pneumoniae and respiratory syncytial virus.

†† Defined as the modified Halm’s criteria (): temperature ≤37,2°C, heart rate ≤100 bpm, systolic blood pressure ≤90 mmHg, respiratory rate ≤ 24 bpm, and oxygen saturation ≥90% for the entire day.

Differential abundance of monocyte/dendritic cell and unconventional T cell metaclusters characterize the immune response in CAP

We identified 86 metaclusters (MCs) in PBMCs by 36-color spectral flow cytometry followed by unsupervised clustering using FlowSOM () (Supplementary Figures 2, 3). We classified these MCs into seven parent lineages: monocytes and dendritic cells (DCs), natural killer (NK) cells, B cells, CD4 T cells, CD8 T cells, unconventional T cells, and granulocytes (Figures 1B–E and Supplementary Figure 4A). The proportions of these lineages were similar between patients and controls, except for unconventional T cells, which were low in abundance in both groups, yet relatively less abundant in CAP patients (median for CAP, 1.8% [IQR, 1.2-2.6], median for control, 2.2% [IQR, 1.6-3.9]; P = 0.043; Figure 1E and Supplementary Figure 4B). We next compared the abundance of all 86 MCs between patients and controls. This untargeted analysis resulted in 19 (22.1%) significantly differentially abundant MCs (BH-adjusted P <0.05) between groups (Figure 1F). MCs within the monocyte and DC lineage featured prominently in this comparison (7/19 [36.8%] of total MCs significantly different), as did unconventional T cell MCs (5/19 [26.3%]). MCs 50 (neutrophils) and 54 (basophils) are likely low density granulocytes that commonly contaminate PBMC fractions of patients with inflammatory conditions (). These cells comprised only a minor fraction of total cells and will not be discussed further. As alterations in immune phenotypes have been reported when comparing patients with COPD with healthy controls (), we repeated the comparison between patients with CAP and control subjects after excluding patients (n = 15) and controls (n = 2) with COPD. The effect sizes in this comparison correlated almost perfectly to the effect sizes of the whole cohort (Pearson’s r = 0.97, P <0.0001), indicating that the higher prevalence of COPD in patients with CAP did not confound the comparison of circulating immune cell phenotypes between CAP and control subjects (Supplementary Figure 5). Thus, substantial differences emerged in the relative abundance of immune cell MCs between patients with CAP and controls, driven primarily by monocytes, DCs and unconventional T cells.

Monocytes and DCs express markers of immune suppression

We delineated 17 monocyte/DC MCs (Figure 1D), divided over 8 cell subsets (Figure 2A). The four monocyte/DC MCs that were significantly more abundant in patients with CAP compared with controls were three classical monocytes MCs and one CD14dim monocyte MC, whereas all less abundant MCs within this lineage were DC MCs (Figure 2B; Supplementary Figure 6A for the relative abundance of all monocyte and DC MCs). We examined the functional profile of these MCs by assessing the expression of seven surface markers relevant for immune cell activation and suppression: CD11b, CD11c, CD38, CD95/Fas, human leukocyte antigen – DR isotype (HLA-DR), programmed death 1 (PD-1), and PD-ligand 1 (PD-L1; Supplementary Figure 6B for all markers). First, we compared the expression of these markers on the significantly more abundant classical monocyte MCs (MCs 64, 65 and 76) with the remaining non-significant classical monocyte MCs (MCs 63, 69, 70, 71, 79, and 83; Figure 2C). MCs 64, 65, and 76 all showed increased expression of PD-L1, and to a lesser extent PD-1 and CD38. Expression of HLA-DR was evidently decreased on MC 65, but increased on MC 76. CD11c was decreased on MCs 64 and 65, but slightly increased on MC 76. Concurrent expression of HLA-DR and CD11c on MC 76 may indicate that these cells are pro-inflammatory, whereas the concurrent downregulation of these molecules on MC 65 likely indicates an immunosuppressive phenotype. The differential expression patterns on the significantly more abundant classical monocyte MCs (Figure 2C) may reflect both immune suppression (e.g. reduced HLA-DR and increased PD-L1) and inflammatory activity (e.g. increased CD38 expression), which are both key features of the immunopathology in critically ill patients with sepsis ().

Figure 2

We next contrasted expression patterns of these markers across all monocyte/DC MCs. For this, we compared the expression (as median fluorescence intensity [MFI] per subject) between CAP and control cells within each MC (Figure 2D) and each subset (Figure 2E and Supplementary Figure 7). These analyses revealed that cells from patients with CAP exhibited markedly decreased expression of HLA-DR – even in MC 76, which had an overall high expression of HLA-DR compared with other classical monocyte MCs (Figure 2C) – and CD11c, and increased expression of PD-1 and PD-L1. To a lesser extent, CD11b was reduced, and CD38 and CD95/Fas increased. We then explored whether the expression of these markers corresponds to the relative abundance of monocyte/DC MCs (Supplementary Figure 8). We found linear relationships between the overall expression of PD-L1 and HLA-DR (the MFI for all cells together) and the relative abundance of each MC (quantified as the Hedges’ g effect size compared with control); for example, the more abundant a monocyte/DC MC was in CAP, the higher the expression of PD-L1 (Pearson’s r = 0.82; P <0.0001), and the lower the expression of HLA-DR (r = -0.57; P = 0.018; Figure 2F). Taken together, these results indicate that, when compared with control subjects, monocytes and DCs from patients with CAP display variable signs of inflammation-induced activation, yet clear subset-wide signs of immune suppression (increased PD-L1 and PD-1, lower HLA-DR) and susceptibility to (Fas ligand-mediated) apoptosis (increased CD95/Fas).

Universal upregulation of CCR7 and CXCR5 on monocytes and DCs

Considering their importance in immune cell migration, we next examined the expression of chemokines receptors. The three more abundant classical monocyte MCs all showed increased expression of all measured chemokine receptors (CCR5, CCR6, CCR7, CXCR3, and CXCR5), when compared with the classical monocyte MCs that were not differentially abundant (Figure 3A). However, when comparing the MFI for these chemokine receptors between CAP and controls for each subject within each MC (Figure 3B) and each subset (Figure 3C and Supplementary Figure 9), we found a near-universal upregulation of the chemokine receptors CXCR3, CCR7, and CXCR5 across most monocyte/DC MCs and cell subsets in patients with CAP. In addition, we found strong linear relationships between the relative abundance of monocyte/DC MCs and the overall expression of CCR7 (r = 0.78; P = 0.0002) and CXCR5 (r = 0.80; P = 0.0001; Figure 3D); such relationships were not present for the other chemokine receptors (Supplementary Figure 8). Expression of CCR7 and CXCR5 positively correlated with PD-1, PD-L1 and CD95/Fas, and negatively correlated with CD11c and to a lesser extent HLA-DR (Figure 3E), and thus appear related to the immunosuppressed and apoptosis-susceptible monocyte/DC phenotype described in Figure 2.

Figure 3

Classical monocyte phenotypes are associated with systemic inflammation and clinical outcomes

We next explored whether K-means clustering – an unsupervised machine learning algorithm – could use patterns of surface marker expression on classical monocyte MCs to distinguish subgroups of patients with CAP. The expression of CD11c, CD11b, CD38, CD95/Fas, PD-1, PD-L1, HLA-DR, CCR5, CCR6, CCR7, CXCR3, and CXCR5 on classical monocyte MCs separated patients in two distinct groups (KM-cluster 1, n = 21, and KM-cluster 2, n = 19; Figure 4A and Supplementary Figure 10A). Separation of the two KM-clusters was driven primarily by expression of CCR7 and CXCR5 (Supplementary Figure 10B). Consistently, when we compared the overall classical monocyte phenotype between the two KM-clusters, KM-cluster 2 was predominantly characterized by increased expression of chemokine receptors, except for CCR6 (Figure 4B and Supplementary Figure 10C). Remarkably, KM clustering using only CCR7 or CXCR5 resulted in almost the exact same cluster assignment for patients (Supplementary Table 2), emphasizing the key role of CCR7 and CXCR5 in separating these two groups of patients with CAP.

Figure 4

): temperature ≤37,2°C, heart rate ≤100 bpm, systolic blood pressure ≤90 mmHg, respiratory rate ≤ 24 bpm, and oxygen saturation ≥90% for the entire day. The hazard ratio (HR) for KM-cluster 2 (compared with KM-cluster 1) from the Cox proportional hazards model is displayed in the graph. There were no competing risks for reaching clinical stability or discharge in this cohort (i.e. no patients were transferred or died prior to the event taking place). HLA-DR, human leukocyte antigen – DR isotype; IL(-RA), interleukin (receptor antagonist); PD-(L)1, programmed death (ligand) 1.

Patients in KM-cluster 2 had a longer time from start of symptoms to hospital admission (median, 4 vs 2 days; P = 0.03), and higher admission C-reactive protein (median, 235 vs 74 mg/L; P = 0.03; Supplementary Table 3). We next compared host response plasma biomarkers between these clusters of patients. While limited in statistical power due to sample size, the overall pattern suggested an exaggerated inflammatory state in KM-cluster 2 – consistent with the elevated C-reactive protein at hospital admission – with significantly higher IL-6 and IL-1RA levels (Figure 4C and Supplementary Figure 10D). Patients in KM-cluster 2 had an increase, albeit not statistically significant, in time to clinical stability or discharge (median, 6 vs 3 days; hazard ratio 0.56; 95% confidence interval 0.29-1.07; P = 0.08; Figure 4D). Together, this exploratory analysis provides preliminary evidence that classical monocyte phenotypes – in particular high expression of CCR7 and CXCR5 – are associated with the inflammatory status and disease course in patients with CAP.

Other differentially abundant metaclusters in patients with CAP

Besides monocytes and DCs, other differentially abundant MCs – in the untargeted comparison of patients with CAP versus controls (Figure 1F) – included five unconventional T cell MCs: three less abundant TCRγδ T cell MCs (MCs 18, 24, 47), one less abundant double-positive T cell MC (MC 10), and one more abundant double-negative T cell MC (MC 17) with an activated phenotype (high HLA-DR, CD38 and CD95 expression; Figure 5A). Double-negative T cells are heterogeneous and can be either TCRαβ+ or TCRγδ+ and may be both immunosuppressive and proinflammatory (, ). While data on double-negative T cells in pneumonia are limited, increased circulating numbers have been reported to correlate with disease severity in pediatric pneumonia () – specifically the CD3low subset, which may correspond to the relatively low CD3 expression in MC 17 – and also in COVID-19 ().

Figure 5

For the remaining differentially abundant MCs (Figure 1F), we focused on those increased in patients with CAP. These included a CD4 central memory T cell cluster with an activated and possibly exhausted phenotype (high CD38, HLA-DR, PD-1, and CD95; Figure 5B); a plasmablast MC (MC 84, Figure 5C); and a B cell MC (MC 38) with high CXCR3 and high IgG, possibly reflecting an activated, class-switched pre-plasmablast/plasma-cell memory B cell (Figure 5D) (). These results highlight that cells throughout the circulating immune system are altered in patients with CAP, and may point towards a hitherto underexplored role for unconventional T cells in CAP pathophysiology.

Immune features largely overlap between patients with CAP and COVID-19

We next assessed to what degree peripheral immune features in patients hospitalized for CAP overlap with those in 35 patients hospitalized for COVID-19, with – similar to CAP – a low to moderate severity of disease (Supplementary Table 4). The proportions of lineages between COVID-19 and controls were similar (Supplementary Figure 11A). Differences in the relative abundance of MCs between COVID-19 patients and controls were largely in line with the COVID-19 peripheral immune signature reported in literature () (Supplementary Figure 11B). Specifically, for the innate immune system, this included an increase in classical monocytes (MCs 64 and 65) and neutrophils (MC 50, likely low density), with a concurrent decrease in non-classical and intermediate monocytes (MCs 72 and 73) and DCs (MCs 66 and 80); for the adaptive immune system this included an increase in plasmablasts (MC 86) and activated effector CD8 T cells (MC 33).

We assessed whether the immune features documented in patients with CAP were also present in patients with COVID-19. Interestingly, the majority of MCs with a biologically relevant difference compared with controls (defined as Hedges’ g ≥0.5 versus controls in CAP or COVID-19) changed in the same direction for CAP and COVID-19 (26/30 [86.7%]), indicating a large overlap in peripheral immune response (Figure 6A). Indeed, the proportional changes across all MCs strongly correlated between patients with CAP and patients with COVID-19 (Pearson’s r = 0.65; P <0.001). Directly comparing CAP and COVID-19 in an untargeted manner resulted in a modest number (8/86 [9.3%]) of differentially abundant MCs (with BH-adjusted P-values mostly very close to 0.05), including higher frequencies of several CD4, CD8 and TCRγδ T cell MCs (MCs 12, 14, 18, 19, 21, and 24) in COVID-19 and higher frequencies of DCs (MC 66) and non-classical monocytes (MC 73) in CAP (Figure 6B). Together, these data show that the proportional change of circulating immune cell frequencies is highly similar between patients with CAP and patients with COVID-19.

Figure 6

Finally, we investigated whether key findings in monocytes and DCs in CAP were also present in COVID-19. Comparing COVID-19 and controls at the level of MFI for each subject within each MC and each subset revealed expression patterns highly comparable to CAP, such as decreased HLA-DR and CD11c, and increased PD-L1, PD-1, CD95/Fas, CD38 and CXCR5 (Figures 6C, D). However, in clear contrast with CAP, classical monocytes of patients with COVID-19 did not express more CCR7 than controls (Figures 6C, D). In summary, these comparisons show that the immune response – at the level of circulating immune cell frequencies and monocyte/DC surface marker expression – is highly similar between CAP and COVID-19, with the notable exception of CCR7 expression on classical monocytes.

Discussion

In this study we aimed to characterize the peripheral immune response in patients with CAP at single-cell protein levels using 36-color spectral flow cytometry. Compared with matched controls without signs of acute infection, the most profound alterations in patients with CAP were in the frequency and phenotype of monocyte and DC MCs. Specifically, we demonstrate both an increased frequency in distinct classical monocyte MCs, as well as phenotypic alterations across virtually all classical monocyte MCs reminiscent of immunosuppression (low HLA-DR, high PD-1 and PD-L1) and clear upregulation of chemokine receptors, most notably CCR7 and CXCR5. In an exploratory analysis, expression patterns on classical monocytes, driven by CCR7 and CXCR5, separated patients into two KM-clusters with diverging inflammatory status and disease course. The peripheral immune response in patients with COVID-19 was highly comparable to patients with CAP, yet lacked the increased expression of CCR7 on classical monocytes.

Immunosuppression is a key feature of the immunopathology of sepsis (, ). Reduced expression of HLA-DR (important for antigen presentation) and increased expression of PD-L1 (an inhibitory immune checkpoint protein) on monocytes are used as surrogates for an immunosuppressive state and, in septic shock, have been linked to reduced cytokine secretion upon stimulation, secondary infections, and mortality (). In line with previous studies by our group (, , ), these indicators of immunosuppression were already present at hospital admission in this cohort of patients with CAP of only moderate disease severity. It remains to be established whether signs of immunosuppression in patients hospitalized with moderately severe CAP relate to clinical outcomes such as new infections (, ), an association that has been made in patients with septic shock (, ). Our understanding of the association between immunological profiles and clinical consequences of pneumonia would be enhanced by future longitudinal studies documenting immune and clinical readouts before, during, and after pneumonia, although such investigations are logistically and methodologically challenging.

We report a remarkable upregulation in the expression of CCR7 and CXCR5 across virtually all circulating classical monocyte MCs. Additionally, separation of KM-clusters 1 and 2 was driven by these chemokine receptors, and high expression of CCR7 and CXCR5 on classical monocytes in KM-cluster 2 was associated with systemic inflammation and a clear trend towards a longer time to clinical stability. Moreover, enhanced CCR7 expression on classical monocytes was restricted to CAP and not present in COVID-19. CCR7 is best known for its role in the migration of DCs and T cells to lymph nodes for antigen presentation, and CXCR5 for migration of B and T cells to lymphoid follicles (, ). In agreement with our results, exposure of human monocytes to toll-like receptor agonists, β-adrenergic agents and supernatants of activated platelets, increased CCR7 and CXCR5 expression (, ). Particularly monocyte CCR7 may impact the host response during bacterial infection. CCR7 is considered a marker for pro-inflammatory M1 macrophages (), and binding of CCR7 on monocytes/macrophages by its ligands CCL19 or CCL21 may potentiate the secretion of pro-inflammatory mediators such as tumor necrosis factor and IL-8 (, ). Patients with sepsis showed increased CCR7 expression on circulating monocytes (, ), and a bioinformatics analysis of whole blood leukocyte RNA expression in datasets of patients with CAP-induced sepsis implicated CCR7 as a gene driving sepsis development (). In mice with Pseudomonas pneumonia CCR7 deficiency resulted in a stronger proinflammatory response and a more efficient clearance of bacteria from the lungs (). Collectively, these results identify monocyte CCR7 as a potentially important player in the host response during pneumonia.

This study is the first to comprehensively characterize PBMC immunophenotypes at single cell resolution in patients with CAP. We used two matched control groups, including patients with COVID-19, which allowed identification of common and distinctive immune features. Limitations include that measurements were performed at a single time point and in a single compartment (i.e. blood), and the lack of functional testing of cell subsets/MCs. Furthermore, while our panel was extensive, it does not capture all rare cell subsets. Although absence of the confounding influence of dexamethasone treatment in our comparison of CAP with COVID-19 can be considered a strength, results obtained in the latter group may be less generalizable to patients currently hospitalized with COVID-19. The K-means clustering analysis was aimed at linking immune phenotypes to clinical characteristics, but should be considered exploratory and requires validation in a larger cohort.

We provide a comprehensive immunological map of PBMCs in patients with moderately severe CAP at single-cell protein level, revealing monocyte MCs with immunosuppressive features and enhanced chemokine receptor expression, and disclosing increased classical monocyte CCR7 expression as a feature unique to CAP when compared with COVID-19. Our results may offer an entry point toward developing novel immunomodulatory treatments of patients with CAP.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by the medical ethical committee of the Amsterdam University Medical Centers. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

TR: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing. ARS: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing. JV: Conceptualization, Data curation, Formal Analysis, Methodology, Software, Visualization, Writing – review & editing. MB: Formal Analysis, Writing – review & editing. RD: Investigation, Writing – review & editing. DF: Investigation, Writing – review & editing. AP: Data curation, Formal Analysis, Investigation, Software, Writing – review & editing. WW: Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing. AS: Conceptualization, Formal Analysis, Investigation, Supervision, Writing – review & editing. JG: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Visualization, Writing – review & editing. TP: Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. TR is supported by research program “NACTAR,” project “MDR-phage” (grant number 16447) which is financed by the Dutch research council (NWO). TP is supported by grants from the Innovative Medicines Initiative (IMI) and the European Union (COMBACTE-CARE and COMBACTE-MAGNET).

Acknowledgments

We are grateful to Dr. A. Davison (Cytek Biosciences) for discussions and continuous support on data quality control of our spectral flow cytometry dataset.

Conflict of interest

Author AP was employed by the company Cytek Biosciences, Inc.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1260283/full#supplementary-material

References

  • 1

    QuintonLJWalkeyAJMizgerdJP. Integrative physiology of pneumonia. Physiol Rev (2018) 98:1417–64. doi: 10.1152/physrev.00032.2017

  • 2

    SalahHMMinhasAMKKhanMSPandeyAMichosEDMentzRJet al. Causes of hospitalization in the USA between 2005 and 2018. Eur Heart J Open (2021) 1. doi: 10.1093/ehjopen/oeab001

  • 3

    AbbafatiCAbbasKMAbbasi-KangevariMAbd-AllahFAbdelalimAAbdollahiMet al. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet (2020) 396:1204–22. doi: 10.1016/S0140-6736(20)30925-9

  • 4

    TorresACillonizCNiedermanMSMenéndezRChalmersJDWunderinkRGet al. Pneumonia. Nat Rev Dis Primers (2021) 7:128. doi: 10.1038/s41572-021-00259-0

  • 5

    ReijndersTDYSchuurmanARvan der PollT. The immune response to respiratory viruses: from start to memory. Semin Respir Crit Care Med (2021) 42:759–70. doi: 10.1055/S-0041-1736459

  • 6

    TianYCarppLNMillerHERZagerMNewellEWGottardoR. Single-cell immunology of SARS-CoV-2 infection. Nat Biotechnol (2022) 40:3041. doi: 10.1038/s41587-021-01131-y

  • 7

    van de VeerdonkFLGiamarellos-BourboulisEPickkersPDerdeLLeavisHvan CrevelRet al. A guide to immunotherapy for COVID-19. Nat Med (2022) 28:3950. doi: 10.1038/s41591-021-01643-9

  • 8

    SchuurmanARReijndersTDYSarisAMoralIRSchinkelMde BrabanderJet al. Integrated single-cell analysis unveils diverging immune features of covid-19, influenza, and other community-acquired pneumonia. Elife (2021) 10. doi: 10.7554/ELIFE.69661

  • 9

    SchuurmanARReijndersTDYvan EngelenTSRLeopoldVde BrabanderJvan LingeCet al. The host response in different aetiologies of community-acquired pneumonia. EBioMedicine (2022) 81:104082. doi: 10.1016/J.EBIOM.2022.104082

  • 10

    SarisAReijndersTDYNossentEJSchuurmanARVerhoeffJVan AstenSet al. Distinct cellular immune profiles in the airways and blood of critically ill patients with COVID-19. Thorax (2021) 76:1010–9. doi: 10.1136/THORAXJNL-2020-216256

  • 11

    ParkLMLanniganJJaimesMC. OMIP-069: forty-color full spectrum flow cytometry panel for deep immunophenotyping of major cell subsets in human peripheral blood. Cytometry Part A (2020) 97:1044–51. doi: 10.1002/CYTO.A.24213

  • 12

    van GassenSCallebautBvan HeldenMJLambrechtBNDemeesterPDhaeneTet al. FlowSOM: Using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A (2015) 87:636–45. doi: 10.1002/CYTO.A.22625

  • 13

    HedgesLV. Distribution theory for glass’s estimator of effect size and related estimators. J Educ Stat (1981) 6:107. doi: 10.2307/1164588

  • 14

    HalmEAFineMJMarrieTJColeyCMKapoorWNObroskyDSet al. Time to clinical stability in patients hospitalized with community-acquired pneumonia: implications for practice guidelines. JAMA (1998) 279:1452–7. doi: 10.1001/JAMA.279.18.1452

  • 15

    SchenzJObermaierMUhleSWeigandMAUhleF. Low-density granulocyte contamination from peripheral blood mononuclear cells of patients with sepsis and how to remove it – A technical report. Front Immunol (2021) 12:684119/BIBTEX. doi: 10.3389/FIMMU.2021.684119/BIBTEX

  • 16

    BhatTAPanzicaLKalathilSGThanavalaY. Immune dysfunction in patients with chronic obstructive pulmonary disease. Ann Am Thorac Soc (2015) 12:S169–75. doi: 10.1513/AnnalsATS.201503-126AW

  • 17

    van der PollTShankar-HariMWiersingaWJ. The immunology of sepsis. Immunity (2021) 54:2450–64. doi: 10.1016/j.immuni.2021.10.012

  • 18

    WuZZhengYShengJHanYYangYPanHet al. CD3+CD4-CD8- (Double-negative) T cells in inflammation, immune disorders and cancer. Front Immunol (2022) 13:816005. doi: 10.3389/fimmu.2022.816005

  • 19

    VoelklSGaryRMackensenA. Characterization of the immunoregulatory function of human TCR-αβ+ CD4- CD8- double-negative T cells. Eur J Immunol (2011) 41:739–48. doi: 10.1002/eji.201040982

  • 20

    WangYLuWLiASunZWangL. Elevated CD3low double negative T lymphocyte is associated with pneumonia and its severity in pediatric patients. PeerJ (2018) 2018:117. doi: 10.7717/peerj.6114

  • 21

    Ekşioğlu-DemiralpEAlanSSiliUBakanDOcakİYürekliRet al. Peripheral innate and adaptive immune cells during COVID-19: Functional neutrophils, pro-inflammatory monocytes, and half-dead lymphocytes. Cytometry B Clin Cytom (2022) 102:153–67. doi: 10.1002/cyto.b.22042

  • 22

    MuehlinghausGCiglianoLHuehnSPeddinghausALeyendeckersHHauserAEet al. Regulation of CXCR3 and CXCR4 expression during terminal differentiation of memory B cells into plasma cells. Blood (2005) 105:3965–71. doi: 10.1182/BLOOD-2004-08-2992

  • 23

    ArunachalamPSWimmersFMokCKPPereraRAPMScottMHaganTet al. Systems biological assessment of immunity to mild versus severe COVID-19 infection in humans. Science (2020) 369:1210–20. doi: 10.1126/science.abc6261

  • 24

    MathewDGilesJRBaxterAEOldridgeDAGreenplateARWuJEet al. Deep immune profiling of COVID-19 patients reveals distinct immunotypes with therapeutic implications. Science (2020) 369:129. doi: 10.1126/science.abc8511

  • 25

    WilkAJRustagiAZhaoNQRoqueJMartínez-ColónGJMcKechnieJLet al. A single-cell atlas of the peripheral immune response in patients with severe COVID-19. Nat Med (2020) 26:1070–6. doi: 10.1038/s41591-020-0944-y

  • 26

    AhernDJAiZAinsworthMAllanCAllcockAAngusBet al. A blood atlas of COVID-19 defines hallmarks of disease severity and specificit. Cell (2022) 185. doi: 10.1016/J.CELL.2022.01.012

  • 27

    Schulte-SchreppingJReuschNPaclikDBaßlerKSchlickeiserSZhangBet al. Severe COVID-19 is marked by a dysregulated myeloid cell compartment. Cell (2020) 182:14191440.e23. doi: 10.1016/J.CELL.2020.08.001

  • 28

    JouanYGuillonAGonzalezLPerezYBoisseauCEhrmannSet al. Phenotypical and functional alteration of unconventional T cells in severe COVID-19 patients. J Exp Med (2020) 217. doi: 10.1084/JEM.20200872

  • 29

    Kuri-CervantesLPampenaMBMengWRosenfeldAMIttnerCAGWeismanARet al. Comprehensive mapping of immune perturbations associated with severe COVID-19. Sci Immunol (2020) 5. doi: 10.1126/SCIIMMUNOL.ABD7114

  • 30

    AgratiCSacchiABordoniVCiminiENotariSGrassiGet al. Expansion of myeloid-derived suppressor cells in patients with severe coronavirus disease (COVID-19). Cell Death Differ (2020) 27:3196–207. doi: 10.1038/S41418-020-0572-6

  • 31

    LaingAGLorencAdel Molino del BarrioIDasAFishMMoninLet al. A dynamic COVID-19 immune signature includes associations with poor prognosis. Nat Med (2020) 26:1623–35. doi: 10.1038/S41591-020-1038-6

  • 32

    VenetFMonneretG. Advances in the understanding and treatment of sepsis-induced immunosuppression. Nat Rev Nephrol (2018) 14:121–37. doi: 10.1038/nrneph.2017.165

  • 33

    TorresLKPickkersPvan der PollT. Sepsis-induced immunosuppression. Annu Rev Physiol (2022) 84:157–81. doi: 10.1146/ANNUREV-PHYSIOL-061121-040214

  • 34

    ZhangYLiJLouJZhouYBoLZhuJet al. Upregulation of programmed death-1 on T cells and programmed death ligand-1 on monocytes in septic shock patients. Crit Care (2011) 15:19. doi: 10.1186/cc10059

  • 35

    GuignantCLepapeAHuangXKheroufHDenisLPoitevinFet al. Programmed death-1 levels correlate with increased mortality, nosocomial infection and immune dysfunctions in septic shock patients. Crit Care (2011) 15:111. doi: 10.1186/cc10112

  • 36

    ShaoRFangYYuHZhaoLJiangZLiCS. Monocyte programmed death ligand-1 expression after 3-4 days of sepsis is associated with risk stratification and mortality in septic patients: A prospective cohort study. Crit Care (2016) 20:110. doi: 10.1186/s13054-016-1301-x

  • 37

    MonneretGVenetF. Sepsis-induced immune alterations monitoring by flow cytometry as a promising tool for individualized therapy. Cytometry B Clin Cytom (2016) 90:376–86. doi: 10.1002/cyto.b.21270

  • 38

    BrandsXHaakBWKlarenbeekAMButlerJUhelFQinWet al. An epigenetic and transcriptomic signature of immune tolerance in human monocytes through multi-omics integration. Genome Med (2021) 13:117. doi: 10.1186/s13073-021-00948-1

  • 39

    BrandsXHaakBWKlarenbeekAMOttoNAFaberDRLutterRet al. Concurrent immune suppression and hyperinflammation in patients with community-acquired pneumonia. Front Immunol (2020) 11:796/PDF. doi: 10.3389/FIMMU.2020.00796/PDF

  • 40

    CharoIFRansohoffRM. The many roles of chemokines and chemokine receptors in inflammation. New Engl J Med (2006) 354:610–21. doi: 10.1056/nejmra052723

  • 41

    ShiCPamerEG. Monocyte recruitment during infection and inflammation. Nat Rev Immunol (2011) 11:762–74. doi: 10.1038/nri3070

  • 42

    HalvorsenBSmedbakkenLMMichelsenAESkjellandMBjerkeliVSagenELet al. Activated platelets promote increased monocyte expression of CXCR5 through prostaglandin E2-related mechanisms and enhance the anti-inflammatory effects of CXCL13. Atherosclerosis (2014) 234:352–9. doi: 10.1016/j.atherosclerosis.2014.03.021

  • 43

    ParadisABernierSDumaisN. TLR4 induces CCR7-dependent monocytes transmigration through the blood-brain barrier. J Neuroimmunol (2016) 295–296:12–7. doi: 10.1016/j.jneuroim.2016.03.019

  • 44

    JaguinMHoulbertNFardelOLecureurV. Polarization profiles of human M-CSF-generated macrophages and comparison of M1-markers in classically activated macrophages from GM-CSF and M-CSF origin. Cell Immunol (2013) 281:5161. doi: 10.1016/J.CELLIMM.2013.01.010

  • 45

    DamåsJKSmithCØieEFevangBHalvorsenBWæhreTet al. Enhanced expression of the homeostatic chemokines CCL19 and CCL21 in clinical and experimental atherosclerosis: Possible pathogenic role in plaque destabilization. Arterioscler Thromb Vasc Biol (2007) 27:614–20. doi: 10.1161/01.ATV.0000255581.38523.7c

  • 46

    FevangBYndestadADamåsJKHalvorsenBHolmAMBeiskeKet al. Chemokines and common variable immunodeficiency; Possible contribution of CCL19, CCL21 and CCR7 to immune dysregulation. Clin Exp Immunol (2009) 158:237–45. doi: 10.1111/j.1365-2249.2009.04013.x

  • 47

    FaivreVLukaszewiczACAlvesACharronDPayenDHaziotA. Accelerated in vitro differentiation of blood monocytes into dendritic cells in human sepsis. Clin Exp Immunol (2007) 147:426–39. doi: 10.1111/J.1365-2249.2006.03287.X

  • 48

    Flores-MejíaLACabrera-RiveraGLFerat-OsorioEMancilla-HerreraITorres-RosasRBoscó-GarateIBet al. Function is dissociated from activation-related immunophenotype on phagocytes from patients with SIRS/sepsis syndrome. Shock (2019) 52:E68–75. doi: 10.1097/SHK.0000000000001314

  • 49

    LiangGLiJPuSHeZ. Screening of sepsis biomarkers based on bioinformatics data analysis. J Healthc Eng (2022) 2022. doi: 10.1155/2022/6788569

  • 50

    EppertBLMotzGTWorthamBWFluryJLBorchersMT. CCR7 deficiency leads to leukocyte activation and increased clearance in response to pulmonary Pseudomonas aeruginosa infection. Infect Immun (2010) 78:2099–107. doi: 10.1128/IAI.00962-09

Summary

Keywords

pneumonia, immunophenotyping, host response, COVID-19, spectral flow cytometry, monocytes, immunosuppression, chemokine receptors

Citation

Reijnders TDY, Schuurman AR, Verhoeff J, van den Braber M, Douma RA, Faber DR, Paul AGA, Wiersinga WJ, Saris A, Garcia Vallejo JJ and van der Poll T (2023) High-dimensional phenotyping of the peripheral immune response in community-acquired pneumonia. Front. Immunol. 14:1260283. doi: 10.3389/fimmu.2023.1260283

Received

17 July 2023

Accepted

30 October 2023

Published

24 November 2023

Volume

14 - 2023

Edited by

Patricia Fitzgerald-Bocarsly, Rutgers, The State University of New Jersey, United States

Reviewed by

Hong Yien Tan, Xiamen University, Malaysia; Katrina Traber, Boston University, United States

Updates

Copyright

*Correspondence: Tom D. Y. Reijnders,

†These authors share first authorship

‡These authors share senior authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics