Impact Factor 6.429

The 5th most cited journal in Immunology

Original Research ARTICLE

Front. Immunol., 12 September 2017 | https://doi.org/10.3389/fimmu.2017.01130

Mass Cytometry Identifies Distinct Lung CD4+ T Cell Patterns in Löfgren’s Syndrome and Non-Löfgren’s Syndrome Sarcoidosis

imageYlva Kaiser1*, imageTadepally Lakshmikanth2,3, imageYang Chen2,3, imageJaromir Mikes2,3, imageAnders Eklund1, imagePetter Brodin2,3,4, imageAdnane Achour2,3 and imageJohan Grunewald1
  • 1Respiratory Medicine Unit, Department of Medicine, Solna and Center for Molecular Medicine, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden
  • 2Science for Life Laboratory, Department of Medicine, Karolinska Institutet, Stockholm, Sweden
  • 3Department of Infectious Diseases, Karolinska University Hospital, Stockholm, Sweden
  • 4Department of Neonatology, Karolinska University Hospital, Stockholm, Sweden

Sarcoidosis is a granulomatous disorder of unknown etiology, characterized by accumulation of activated CD4+ T cells in the lungs. Disease phenotypes Löfgren’s syndrome (LS) and “non-LS” differ in terms of clinical manifestations, genetic background, HLA association, and prognosis, but the underlying inflammatory mechanisms largely remain unknown. Bronchoalveolar lavage fluid cells from four HLA-DRB1*03+ LS and four HLA-DRB1*03 non-LS patients were analyzed by mass cytometry, using a panel of 33 unique markers. Differentially regulated CD4+ T cell populations were identified using the Citrus algorithm, and t-stochastic neighborhood embedding was applied for dimensionality reduction and single-cell data visualization. We identified 19 individual CD4+ T cell clusters differing significantly in abundance between LS and non-LS patients. Seven clusters more frequent in LS patients were characterized by significantly higher expression of regulatory receptors CTLA-4, PD-1, and ICOS, along with low expression of adhesion marker CD44. In contrast, 12 clusters primarily found in non-LS displayed elevated expression of activation and effector markers HLA-DR, CD127, CD39, as well as CD44. Hierarchical clustering further indicated functional heterogeneity and diverse origins of T cell receptor Vα2.3/Vβ22-restricted cells in LS. Finally, a near-complete overlap of CD8 and Ki-67 expression suggested larger influence of CD8+ T cell activity on sarcoid inflammation than previously appreciated. In this study, we provide detailed characterization of pulmonary T cells and immunological parameters that define separate disease pathways in LS and non-LS. With direct association to clinical parameters, such as granuloma persistence, resolution, or chronic inflammation, these results provide a valuable foundation for further exploration and potential clinical application.

Introduction

Sarcoidosis is a multisystem granulomatous disorder of unknown etiology, characterized by accumulation of activated CD4+ T cells in the lungs (1). Löfgren’s syndrome (LS) is distinguished by an acute disease onset with fever, bilateral hilar lymphadenopathy (BHL), erythema nodosum, and/or ankle arthritis, as well as a good prognosis and typically self-limiting disease course (2). In contrast, “non-LS” sarcoidosis constitutes a more heterogeneous patient group that is usually characterized by an insidious onset and disease progression, with substantial risk of developing chronic disease and pulmonary fibrosis. LS is associated with carriage of the HLA-DRB1*03 allele (3), which also strongly influences spontaneous disease resolution (4). The presence of clonal T cell receptor (TCR)-restricted cells in HLA-DRB1*03+ patients further implicates specific antigen recognition in the lungs (5, 6).

Detailed genetic profiling indicates that LS and non-LS share much fewer features than previously acknowledged (7), and arguments have been made for LS to be considered a separate disease entity. However, understanding of pulmonary immune pathways underlying the observed clinical heterogeneity remains incomplete and, until recently, methodology for adequately addressing this question has been lacking.

In this study, we, therefore, subjected bronchoalveolar lavage fluid (BALF) cells from four HLA-DRB1*03+ LS and four HLA-DRB1*03 non-LS patients (Table 1), which are believed to represent two clinical opposites, to mass cytometry, utilizing a T cell panel of 33 unique markers. Unbiased data analysis was performed using two different clustering algorithms to identify differentially regulated cell populations. Specific attention was paid to the recently identified subpopulation preferentially expressing TCR variable segments Vα2.3 and Vβ22 in HLA-DRB1*03+ LS patients (8). With the aim of identifying novel cell populations and pathways that associate with disease resolution or progression, this study demonstrates the sensitivity of mass cytometry, its ability to detect significant differences between patient groups despite limited sample sizes, and the benefit of its application in complex diseases. Albeit in a small study population, we here present novel data on immunological markers that differentiate between sarcoidosis subtypes at the molecular level, and that can be targeted by future exploration in larger patient cohorts.

TABLE 1
www.frontiersin.org

Table 1. Clinical characteristics of sarcoidosis patients.

Materials and Methods

Study Subjects, Bronchoscopy, and Bronchoalveolar Lavage (BAL)

Bronchoscopy with BAL was performed as previously described (9). Eight newly diagnosed sarcoidosis patients (one female) with a median age of 43 years (Table 1) were included in the study. All patients were HLA-typed and diagnosed with sarcoidosis according to criteria established by the World Association of Sarcoidosis and Other Granulomatous Disorders (WASOG) (10). Specifically, these included typical clinical and radiographic manifestations, findings at bronchoscopy with BAL including an elevated CD4/CD8 ratio and, if required, positive biopsies, as well as exclusion of other diagnoses. Four patients were diagnosed with LS, defined as an acute onset, chest radiographic findings of BHL alone (n = 2), or BHL with pulmonary infiltrates (n = 2), and erythema nodosum and/or bilateral ankle arthritis. Four patients presented with non-LS sarcoidosis, with insidious onset, fatigue, dry cough, and parenchymal changes. In addition, BALF cells from one healthy volunteer and peripheral blood mononuclear cells (PBMCs) from three healthy donors were used as staining controls. Informed consent was obtained from all subjects and ethical approval granted from the Stockholm County Regional Ethical Committee (approval numbers: 2005/1031-31/2, 2009/20-32, and 2011/35-32).

HLA Typing

Genomic DNA was extracted from whole blood samples of patients and healthy volunteers. HLA-DRB1 and DRB3 alleles were subsequently determined by the PCR-sequence-specific primer technique (Olerup SSP-DR Low Resolution Kit, Saltsjöbaden, Sweden) as previously described (11).

Assessment of TCR Expression by Flow Cytometry

Freshly isolated BALF cells were stained ex vivo using the following antibodies: CD3-Pacific Blue, clone UCHT1 (BD Pharmingen, San Diego, CA, USA), CD4-APC-H7, clone SK3 (BD Biosciences, San Jose, CA, USA), Vα2.3-FITC, clone F1 (Thermo Scientific, Rockford, IL, USA), and Vβ22-PE, clone IMMU 546 (Beckman Coulter Immunotech, Marseille, France). Live/Dead Fixable Aqua Dead Cell Stain Kit (Life Technologies, Eugene, OR, USA) was used for assessment of cell viability. Cells were sequentially gated on lymphocytes (based on FSC vs. SSC), single cells (based on FSC-A vs. FSC-H), viable cells (defined as Aqua negatively stained cells), CD3+, and CD4+ cells. The CD4+ gate was set as threshold gate for acquisition, with a minimum of 15,000 events being collected. Flow cytometry was run on a BD FACSVerse (Beckton Dickinson, San Jose, CA, USA) and results were analyzed using FlowJo X (TreeStar, Ashland, OR, USA) software.

Mass Cytometry

Cryopreserved BALF cells obtained from eight sarcoidosis patients (LS, n = 4; non-LS, n = 4) were thawed in RPMI medium supplemented with fetal bovine serum (FBS), penicillin–streptomycin, and benzonase (Sigma-Aldrich, St. Louis, MO, USA). BALF cells from one healthy individual and PBMCs from three healthy donors were included as staining controls but not subjected to bioinformatic analysis. Briefly, for live–dead cell distinction, cells were stained with 2.5 µM Cisplatin (Fluidigm, South San Francisco, CA, USA) in RPMI without serum for 5 min at RT and quenched with RPMI containing FBS. An Agilent robotic platform was used for staining cells, as this offers more uniform staining and less variability among samples. Cells were then re-suspended in CyFACS buffer (PBS with 0.1% BSA, 0.05% sodium azide, and 2 mM EDTA) and counted; approximately 1–2 million live cells were used for staining in a 96-well round-bottom plate. Next, cells were incubated for 30 min at 4°C with a 30 µl cocktail of metal-conjugated antibodies targeting surface antigens in a T cell panel (Table 2). Following wash with CyFACS buffer and overnight fixation using 1% formaldehyde made in PBS (Polysciences Inc., Warminster, PA, USA), cells were permeabilized using intracellular fixation and permeabilization buffer set (eBiosciences Inc., San Diego, CA, USA) according to the manufacturer’s recommendations and stained with 30 µl of intracellular antibody cocktail (Ki-67) for 60 min at RT. Cells were washed and fixed in 1% formaldehyde at 4°C until acquisition. On the day of acquisition (within a week after staining), cells were stained with DNA intercalator (0.125 µM Iridium-191/193 or MaxPar® Intercalator-Ir; Fluidigm) in 1% formaldehyde made in PBS for 20 min at RT. After multiple washes with CyFACS, PBS, and MilliQ water, cells were filtered through a 35-µm nylon mesh and diluted to 500,000 cells/ml. Cells were acquired at a rate of 300–500 cells/s using a CyTOF2 mass cytometer (Fluidigm), CyTOF software v.6.0.626 with noise reduction, a lower convolution threshold of 200, event length limits of 10–150 pushes, a sigma value of 3, and a flow rate of 0.045 ml/min.

TABLE 2
www.frontiersin.org

Table 2. Mass cytometry staining panel.

Antibodies and Reagents

Purified antibodies were obtained in protein-free buffer and subsequently coupled to isotopically purified metals using MaxPar® antibody conjugation kits (Fluidigm) (12), as per the manufacturer’s protocol. After determining protein concentration by measurement of absorbance at 280 nm, metal-labeled antibodies were diluted in Candor PBS Antibody Stabilization solution (Candor Bioscience, Wangen, Germany) for long-term storage at 4°C. Antibodies used in this study are listed in Table 2.

Dimensionality Reduction, Cluster Analysis, and Visualization

Sarcoid BALF samples were sequentially gated on single cells, viable (DNA+Cisplatin) cells, and CD3+CD4+CD8 T cells. Based on the collective expression of the remaining 28 T cell-specific markers in the panel (Table 2), the Citrus algorithm (13) was used to identify differently abundant T cell populations in LS and non-LS samples, and clusters were visualized by ACCENSE (14). Citrus was run using the original R implementation (v. citrus_0.08) in R v.3.3.1 (2016-06-21) with a random down-sampling to 1,500 manually gated CD4+ T-cells from each file. The following markers were used for cluster analysis: CD45, CD57, CCR6, CD5, CCR5, CD4, CD8a, CD31, ICOS, TCRαβ, CD3e, CCR4, Vα2.3, CXCR3, Vβ22, CD28, CD161, Ki-67, HLA-DR, CD44, CD127, CD27, CD38, CD45RA, CTLA-4, PD-1, CD39, and CXCR5. The collective expression across all markers was used to define cell populations. The “pamr” model was applied for calculating differential regulation, using a relaxed FDR-constrained model for exploratory purposes (FDR < 0.25). Cells from all samples were combined and subjected to hierarchical clustering. Descriptive features of identified clusters were calculated on a per-sample basis to train a regularized regression model predictive of the two patient groups. Classification models were constructed using the nearest shrunken centroid and lasso-regularized logistic regression methods, both of which build a series of predictive models using automatically selected informative subsets of supplied regressors. Internal cross-validation was used to evaluate model fit and select an appropriate threshold for the final model.

Dimensionality reduction of unclustered data was performed using the t-stochastic neighborhood embedding (t-SNE) algorithm implemented in the Cytofkit library (1517), supplied by Bioconductor v.3.4 (18, 19) and run in RStudio v.1.0.44 (2016-11-01) (20). The markers used for clustering were CD45, CD57, CCR6, CD5, CCR5, CD4, CD8a, CD31, ICOS, TCRαβ, CD3e, CCR4, Vα2.3, CXCR3, Vβ22, CD28, CD161, Ki-67, HLA-DR, CD44, CD127, CD27, CD38, CD45RA, CTLA-4, PD-1, CD39, CXCR5, CD19, and CD11c. A fixed number of 5,000 cells were sampled without replacement from each file and combined for analysis. Resulting t-SNE plots were subsequently filtered by marker expression and patient group to visualize differences between LS and non-LS.

Complementary analysis and validation of marker expression in BALF CD4+ T cells from LS and non-LS patients, respectively, as well as from a single healthy individual, were performed using FlowJo X (TreeStar) software, based on the above mentioned gating strategy [i.e., single cells, viable (DNA+Cisplatin) cells, and CD3+CD4+CD8 T cells]. Univariate statistical analysis of marker expression in LS vs. non-LS CD4+ T cells was performed using the non-parametric Mann–Whitney U test in GraphPad Prism v.5.02 software (GraphPad Software, Inc., La Jolla, CA, USA.). p < 0.05 was considered significant.

Results

Abundance of Individual CD4+ T Cell Clusters Differs Significantly between LS and Non-LS

Within the BALF CD4+ T cell population, we identified 19 clusters of cells whose abundance differed significantly between LS and non-LS patients. Seven of these were more abundant in LS, while 12 primarily appeared in non-LS (Figure 1A). Hierarchical clustering showed these 19 populations to assemble in three distinct groups, where clusters in group 1 were dominated by cells from LS patients and defined by a significant reduction in adhesion molecule CD44 compared to background cells (i.e., all CD4+ T cells). In group 2, clusters were distinguished by effector cell phenotypes, and in the case of LS-dominated populations, expression of TCR variable chains Vα2.3 and/or Vβ22. The sole cluster in group 3 showed an elevation in CD57 expression and was significantly reduced in LS patients (Figures 1A,B).

FIGURE 1
www.frontiersin.org

Figure 1. Significant differences in CD4+ T cell patterns between Löfgren’s syndrome (LS) and non-LS. (A) Citrus network tree visualizing the hierarchical relationship between identified bronchoalveolar lavage fluid CD4+ T cell populations in LS (n = 4) and non-LS (n = 4). Circle size reflects number of cells within a given cluster. Clusters differing significantly in abundance between the two conditions are divided into three main groups (highlighted). Populations more abundant in LS are indicated in blue, while those more abundant in non-LS are marked in red. (B) Citrus-generated box plots for four representative and differentially regulated populations in LS and non-LS, illustrative of all three cluster groups shown in panel (A). Cluster abundances are shown on a logarithmic scale for clusters 24990, 24987, 24992, and 24961, respectively, annotated according to panel (A). All differences in abundance are significant at FDR < 0.25.

Regulatory Elements are More Pronounced in LS CD4+ T Cells

Most notably, the regulatory profile of CD4+ T cells was markedly enhanced in LS patients. Citrus analysis and subsequent t-SNE visualization revealed the expression of immune-inhibitory receptors CTLA-4 and PD-1, as well as co-stimulator ICOS to be significantly elevated in LS CD4+ T cells compared to background CD4+ T cells. Most prominently, a twofold to threefold increase in PD-1 expression was observed in LS patients compared to background, while such a shift was virtually absent in non-LS patients (Figures 2A,B). Manual gating of cells using FlowJo revealed higher frequencies of CTLA-4+ (4.9%) and PD-1+ (35.7%) CD4+ T cells in BALF from a healthy individual compared to sarcoidosis patients of both phenotypes. Notably, however, CD4+ T cells from LS patients displayed higher median expression of both receptors (CTLA-4+ 1.8%, PD-1+ 27.8%) compared to non-LS cells (CTLA-4+ 1.4%, PD-1+ 14.6%), indicating an overall downregulation of these receptors in disease, but also that regulatory capacity is retained to a greater extent in LS lungs (Figure 2C; Table 3).

FIGURE 2
www.frontiersin.org

Figure 2. Regulatory elements dominate in Löfgren’s syndrome (LS) CD4+ T cells. (A) t-stochastic neighborhood embedding (t-SNE) plots visualizing distribution and intensity of CTLA-4 and PD-1 expression in LS and non-LS CD4+ T cells, respectively, following clustering of all samples combined and subsequent separation by condition. (B) Citrus-generated histograms showing significant changes in marker distributions in individual clusters (red) compared to background total CD4+ T cells (blue). Representative LS- and non-LS-specific populations are exemplified by clusters 24987 and 24988, respectively. (C) Representative FlowJo contour plots manually gated for single, live (DNA+ Cisplatin), CD3+, CD4+, CD8 T cells in LS and non-LS patients, respectively, showing expression of CTLA-4 and PD-1 in the two patient groups, as well as in a bronchoalveolar lavage fluid sample from a healthy control (HC).

TABLE 3
www.frontiersin.org

Table 3. Summary of CTLA-4 and PD-1 expression in bronchoalveolar lavage fluid (BALF) CD4+ T cells.

Non-LS CD4+ T Cells Consistently Exhibit Effector Profiles

In clusters more abundant in non-LS patients, HLA-DR, CD127, and CD39 were all significantly elevated compared to background CD4+ T cells, consistently displaying a twofold, or in some cases nearly threefold, increase (Figure 3A). In addition, two clusters (24961 and 24966) showed a significant increase in CD57 expression that was not detected in LS-dominated clusters (Figure 3B).

FIGURE 3
www.frontiersin.org

Figure 3. Non-Löfgren’s syndrome (non-LS) CD4+ T cells exhibit a pronounced effector profile. Enhanced influence of HLA-DR, CD127, and CD39 in non-LS CD4+ T cells as visualized by Citrus-generated histograms. A representative cluster (24971) abundant in non-LS patients is shown in panel (A), exemplifying the significantly higher expression of these effector markers in non-LS. Similarly, cluster 24966 illustrates elevated expression of CD57 in one of two clusters prominent in non-LS (B).

Adhesion Marker CD44 is Reduced in LS Compared to Non-LS CD4+ T Cells

Expression of the adhesion molecule CD44, believed to influence granuloma formation and persistence, was significantly reduced in LS CD4+ T cells with an expression onefold to twofold lower than background (Figure 4A). This reduction was consistent across all LS clusters but most prominent in group 1 cells (Figure 1A). Conversely, non-LS cells showed a shift in CD44 expression of similar magnitude, but opposite direction (Figure 4B). The differentially regulated marker expressions that characterize LS and non-LS are summarized in Table 4.

FIGURE 4
www.frontiersin.org

Figure 4. Adhesion marker CD44 is consistently reduced in Löfgren’s syndrome (LS). Expression of CD44 was significantly reduced in clusters more abundant in LS, as exemplified by cluster 24975 (A), compared to non-LS (B) (cluster 24988).

TABLE 4
www.frontiersin.org

Table 4. Summary of differentially regulated marker expression.

TCR Vα2.3+CD4+ T Cells in LS are Functionally Heterogeneous

By t-SNE visualization, a small cluster defined by high expression of TCR segment Vα2.3 was identified in LS patients (Figure 5A). Interestingly, this cluster concomitantly expressed chemokine receptors CXCR3, CCR4, and CCR6 (Figures 5B–D), as well as high levels of activation marker HLA-DR (Figure 5E), CD11c, associated with T cell differentiation and migratory potential (Figure 5F), and CD31, involved in leukocyte adhesion and a marker of recent thymic emigrant cells (Figure 5G). Expression of regulatory markers, e.g., CTLA-4 (Figure 5H), varied with TCR expression. Citrus comparison of Vα2.3+Vβ22 (Figure 5I) and Vα2.3+Vβ22+ (Figure 5J) clusters (24987 and 24983, respectively) revealed abundance of both populations to be significantly increased in LS compared to non-LS patients, but the distribution of receptors associated with cell activation and effector functions (CD127, CD28, CD161) and immune regulation (PD-1, CTLA-4, ICOS) molecules, respectively, to differ between the two Vα2.3+ clusters. Judging by the network tree presented in Figure 1A, these two clusters also appear to have arisen from different origins.

FIGURE 5
www.frontiersin.org

Figure 5. T cell receptor (TCR) Vα2.3+ cells in Löfgren’s syndrome (LS) constitute a functionally heterogeneous population. t-stochastic neighborhood embedding (t-SNE) plots filtered for expression of Vα2.3 (A), CXCR3 (B), CCR4 (C), CCR6 (D), HLA-DR (E), CD11c (F), CD31 (G), and CTLA-4 (H) in LS bronchoalveolar lavage fluid cells. The main Vα2.3+ cluster is highlighted in red. Citrus comparison of Vα2.3+Vβ22 (I) and Vα2.3+Vβ22+ (J) clusters (24987 and 24983, respectively) reveals distribution of effector (e.g., CD127, CD28) and regulatory markers (e.g., PD-1, ICOS) to vary with TCR expression.

Overlapping Expression of CD8 and Ki-67 in Both LS and Non-LS

Previous studies and the established diagnostic criterion of an elevated BALF CD4/CD8 ratio have hitherto assumed CD4+ T cells to be the driving force of disease in sarcoidosis. In both patient groups, however, expression of Ki-67, representing actively proliferating cells, and CD8 coincided remarkably well (Figures 6A,B), indicating an active role of CD8+ T cells in sarcoid inflammation. A significant increase in Ki-67 expression could also be observed for group 2 CD4+ T cell clusters (Figure 1A), i.e., those primarily exhibiting an effector profile, but much less pronounced than for CD8+ T cells.

FIGURE 6
www.frontiersin.org

Figure 6. Near-complete overlap of CD8 and Ki-67 expression. t-stochastic neighborhood embedding (t-SNE) plots revealing an almost perfect overlap in expression of CD8 (A) and Ki-67 (B), the latter being a marker of actively proliferating cells, in both Löfgren’s syndrome (LS) and non-LS patients.

Discussion

While marked differences in clinical presentation and outcome have long defined LS and non-LS, comprehensive genetic disparities were only recently demonstrated (7), and the immunological mechanisms driving these two forms of sarcoidosis are still incompletely understood. To address key characteristics of CD4+ T cells in LS and non-LS in an unbiased fashion, we subjected BALF cells from two genetically and clinically distinct patient groups to mass cytometry and cluster analysis, identifying cell populations of significantly different abundance in the two disease phenotypes. On the individual marker level, significant changes were most notably observed for inhibitory receptors CTLA-4 and PD-1, adhesion marker CD44, effector and activation markers HLA-DR, CD127 and CD39, as well as a range of markers associated with cells expressing TCR variable segments Vα2.3 and/or Vβ22. These findings connect immunological profiles to established clinical characteristics of disease, and the identified novel markers may allow for both prognostic evaluation and more precise differentiation between sarcoidosis subtypes. Despite sample size being limited by methodological requirements and patient availability, these statistically robust differences further argue for LS and non-LS to be considered as immunologically distinct conditions.

Believed to influence tolerance induction and long-term tolerance maintenance, respectively, inhibitory receptors CTLA-4and PD-1 act to constrain development of autoimmunity. Notably, anti-CTLA-4 and anti-PD-1 treatments have the potential to initiate or exacerbate sarcoid-like inflammation (21, 22), reinforcing the notion of altered immune activation and possibly an autoimmune component in sarcoidosis. Broos et al. also previously found reduced CTLA-4 expression on sarcoidosis T cells compared to healthy individuals, suggestive of a shift in balance between immune activation and regulation that results in exacerbated T cell activity and, at length, chronic inflammation, and tissue damage (23). In this study, a similar reduction in CTLA-4 and PD-1 expression could be observed in sarcoidosis compared to a healthy individual. Among sarcoidosis patients, however, expression of both receptors was primarily found in LS CD4+ T cells, indicating a higher degree of immune self-restriction compared to non-LS patients. In addition, ICOS expression was significantly increased in some, though not all, LS clusters compared to background CD4+ T cells. Elevated ICOS expression has previously been demonstrated in LS (24) and could associate with the higher IL-10 production seen in these patients (25, 26). From a clinical perspective, LS patients in general and HLA-DRB1*03+ individuals in particular have a favorable prognosis, and in our opinion, it is reasonable to believe that superior regulatory capacity of CD4+ T cells in LS compared to non-LS influences the characteristic resolving disease course. Accordingly, it is interesting to speculate that stimulation of regulatory elements could have therapeutic potential in non-LS patients prior to progression to a chronic disease state.

In contrast, compared to background BALF CD4+ T cells from non-LS patients displayed an elevated expression of activation markers HLA-DR, CD127, and CD39 consistent across all clusters. Polymorphisms in the CD127 (IL-7R) gene are known risk factors for sarcoidosis development (27), and the IL-7 receptor is crucial during T cell activation, homeostasis, and differentiation. Importantly, it is differentially expressed on effector and regulatory T cells (Tregs), with the latter being defined by low levels of this receptor. The lack of CD127 in LS-dominated clusters, combined with the expression of CTLA-4 and PD-1, suggests immune regulation to be maintained in LS CD4+ T cells, while the balance is shifted in favor of effector functions in non-LS CD4+ T cells. The role of CD39 in this context is less established, but CD39+CD4+ T cells have been proposed as inducers of T cell proliferation and cytokine production (28). Moreover, a recent study showed reduced suppressive function of CD39+ compared to CD39 Tregs (29), implicating dysfunctional regulatory mechanisms particularly in the pathogenesis of non-LS. Notably, CD39+ cells have also been found in sarcoid granulomas (30).

While expression of CD57 is most commonly associated with NK cells, its presence has been documented in T cells with a more advanced differentiation state, reduced replicative capacity but retained cytokine production, particularly of IFNγ (3133), and possibly cytotoxic functions (34). In the present study, two clusters, both of which more abundant in non-LS patients, demonstrated elevated CD57 expression compared to background CD4+ T cells. Higher IFNγ production in non-LS patients has been demonstrated previously (25, 26), and could possibly be at least partly attributed to such terminally differentiated cells with cytotoxic potential.

The adhesion molecule CD44 has been detected on lymphocytes in granulomas of both Crohn’s disease and pulmonary sarcoidosis, and suggested to affect their homing and activation, as well as macrophage differentiation (35). Interestingly, in the present study, CD44 was consistently reduced in LS CD4+ T cell clusters compared to background, while the opposite was observed in non-LS clusters. As T cells are essential for maintenance of granulomatous structures, reduced adhesive capacity of T cells may well play a role in spontaneous resolution of granulomas, thereby influencing disease resolution in LS. Of note, any clinically relevant differences between LS and non-LS granulomas have yet to be investigated, but the prospect that alternate granulomatous development might ultimately affect disease outcome is worth considering. Recently, involvement of the mTORC1 pathway was proposed to promote disease progression in sarcoidosis (36). Intriguingly, mTORC1 complexes can be inhibited by CD44 blockade (37), suggesting a dual role for CD44 in granuloma persistence, acting both during cell–cell adhesion and signal transduction.

The recently noted clonality of Vα2.3+Vβ22+CD4+ T cells (8) indicates them to be a driving force in HLA-DRB1*03+-mediated antigen recognition. The current study further suggests that cells expressing Vα2.3 can initiate and/or act within different immune pathways. Moreover, expression of Vα2.3 together with Vβ22 appear to associate with effector-like markers, while Vα2.3+ populations with other Vβ segment pairings preferentially expressed markers associated with immune regulation, such as PD-1 and ICOS. As hierarchical clustering suggests these clusters to arise from separate origins, it is plausible that several oligoclonal Vα2.3+ populations expand in the lungs in response to different epitopes derived from the same antigen, and that these cells serve discrete purposes during antigen recognition and clearance. Previous observations of both an effector phenotype of Vα2.3+ cells (5, 6, 8, 38), as well as reduced IFNγ and Th1 responses in LS patients (25, 26), could possibly be at least partly attributed to divergent TCRβ chain preferences. Expression of multiple chemokine receptors (CXCR3, CCR4, and CCR6) suggests that Vα2.3+CD4+ T cells might follow different migratory patterns compared to Vα2.3CD4+ T cells, possibly allowing for more efficient interaction with other immune cells. The concurrent high expression of HLA-DR and CD31 might also signify their recent activation and ongoing homing to the site of antigen presentation. Activated, antigen-specific T cells expressing CD11c have been found to associate with both effector and regulatory profiles (3941), and in this study, we show that CD11c+Vα2.3+ cells can express either co-stimulatory or co-inhibitory molecules. Alternatively, this functional discrepancy may reflect a transitional state of activated T cells following antigen recognition, implying that multiple stages of differentiation occur throughout the disease process. A gradual shift toward a more regulatory profile may thus be an inherent characteristic of LS lung CD4+ T cells, ultimately contributing to a better prognosis.

In both patient groups, the near-complete overlap of Ki-67 and CD8 expression conceivably points to a more pronounced activation of CD8+ T cells in sarcoidosis than previously appreciated. It is intriguing to speculate that lack of regulatory capacity in non-LS CD4+ T cells may result in unrestrained proliferation and increased cytotoxic activity of CD8+ T cells, which could promote tissue damage and chronic inflammation. In line with this hypothesis, a stronger cytotoxic response in non-LS CD8+ T cells was recently observed by Brighenti et al. (in preparation). Accordingly, tighter regulation of the CD8+ response by CTLA-4- and PD-1-expressing CD4+ T cells in LS patients could contribute to spontaneous resolution and absence of permanent tissue scarring.

Despite the differences observed, the low sample number is an obvious limitation to the study, as is the inability to perform validation experiments in the same cohort. Mass cytometry is a powerful method for detecting robust differences in small sample sizes, thus making it suitable for unbiased discovery studies. However, the cell number requirement is higher than for flow cytometry, which is a restraint when using clinical samples, and due to sample vaporization, no additional sorting or functional characterization experiments can be performed on the same cells. Moreover, established techniques such as flow cytometry or western blot have lower sensitivity and are, therefore, not ideal means of validation of mass cytometry data. Nevertheless, the current study demonstrates a novel approach to the study of CD4+ T cells in the lung and lays a foundation for future studies, ideally comprising larger groups of patients. The deliberately high FDR threshold for Citrus analysis enables identification of subpopulations that can be further explored in more stringent and comprehensive studies. Of note, several differentially regulated clusters in this study are derived from the same “parent” cluster and, therefore, exhibit similarities in marker expression, thus strengthening the accuracy of identified differences at the individual marker level despite a higher FDR. Future analyses should also aim to expand the current T cell panel to include markers of interest when investigating tissue-resident T cells, such as CD69, or for detailed delineation of Tregs, e.g., CD73 and CD25. Targeting of additional cytokine receptors and a broad range of transcription factors could also contribute to a more comprehensive subclassification of T cells and their role in different disease phenotypes.

Conclusion

Together, these mass cytometry-based findings contribute new information on lung T cell differentiation and activation in two forms of sarcoidosis that may in time be considered as separate diseases. In LS, a stronger regulatory profile along with decreased adhesive capacity of CD4+ T cells, offer a stark contrast to the non-LS effector profile and a direct link to characteristic clinical symptoms. The demonstrated ability to detect significant differences at the molecular level, despite low sample numbers and limited availability of cells from the lungs, could be employed to further distinguish between clinical subtypes, e.g., based on presence of extra-pulmonary symptoms.

In conclusion, mass cytometry provides a powerful and sensitive means of understanding cellular mechanisms of disease that could prove especially valuable in heterogeneous conditions such as non-LS sarcoidosis, and markers that differentiate disease phenotypes could potentially be exploited for future therapeutic purposes in non-resolving disease.

Ethics Statement

This study was carried out in accordance with the recommendations of Karolinska Institutet, Karolinska University Hospital and the Stockholm County Regional Ethical Committee with written informed consent from all subjects in accordance with the Declaration of Helsinki. The protocol was approved by the Stockholm County Regional Ethical Committee.

Data Reposition

Mass cytometry data in the form of FCS files can be accessed from FlowRepository (https://flowrepository.org/) via accession number FR-FCM-ZY9W.

Author Contributions

Conception and study design: YK, AA, and JG. Sample collection: YK and AE. Mass cytometry panel development, optimization, and implementation: TL, JM, and PB. Data pre-processing and cluster analysis (Citrus, ACCENSE): YC and PB. Data visualization (t-SNE): YK. Interpretation and assessment of clinical parameters: YK, AA, and JG. Manuscript preparation: YK, AA, and JG. Critical reading and intellectual assessment of manuscript: YK, TL, AE, PB, AA, and JG. All authors read and approved the final manuscript.

Conflict of Interest Statement

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.

Acknowledgments

The authors thank research nurses Gunnel de Forest, Margitha Dahl, Heléne Blomqvist and biomedical analyst Benita Dahlberg (Respiratory Medicine Unit, Department of Medicine, Solna, Karolinska University Hospital) for skillful assistance in sample preparation and processing. Thanks also to Michael Hagemann-Jensen (Respiratory Medicine Unit, Department of Medicine, Solna, Karolinska Institutet) for valuable input regarding RStudio-based data analysis and interpretation.

Funding

Financial support in the form of research grants was received from The Swedish Heart-Lung Foundation (Hjart-Lungfonden) 20130321, The Swedish Research Council (Vetenskapsradet) K2013-57X-14182-12-3, The Mats Kleberg Foundation, The King Oscar II Jubilee Foundation, The Swedish Order of Freemasons’ King Gustav V and Queen Victoria Foundation, and Karolinska Institutet 3-1818/2013. Support was also provided through the regional agreement on medical training and clinical research (ALF) between Stockholm County Council and Karolinska Institutet.

References

1. Iannuzzi MC, Rybicki BA, Teirstein AS. Sarcoidosis. N Engl J Med (2007) 357(21):2153–65. doi:10.1056/NEJMra071714

CrossRef Full Text | Google Scholar

2. Lofgren S, Lundback H. The bilateral hilar lymphoma syndrome; a study of the relation to tuberculosis and sarcoidosis in 212 cases. Acta Med Scand (1952) 142(4):265–73. doi:10.1111/j.0954-6820.1952.tb13865.x

CrossRef Full Text | Google Scholar

3. Berlin M, Fogdell-Hahn A, Olerup O, Eklund A, Grunewald J. HLA-DR predicts the prognosis in Scandinavian patients with pulmonary sarcoidosis. Am J Respir Crit Care Med (1997) 156(5):1601–5. doi:10.1164/ajrccm.156.5.9704069

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Grunewald J, Eklund A. Lofgren’s syndrome: human leukocyte antigen strongly influences the disease course. Am J Respir Crit Care Med (2009) 179(4):307–12. doi:10.1164/rccm.200807-1082OC

CrossRef Full Text | Google Scholar

5. Grunewald J, Berlin M, Olerup O, Eklund A. Lung T-helper cells expressing T-cell receptor AV2S3 associate with clinical features of pulmonary sarcoidosis. Am J Respir Crit Care Med (2000) 161(3 Pt 1):814–8. doi:10.1164/ajrccm.161.3.9906001

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Planck A, Eklund A, Grunewald J. Markers of activity in clinically recovered human leukocyte antigen-DR17-positive sarcoidosis patients. Eur Respir J (2003) 21(1):52–7. doi:10.1183/09031936.03.00059103

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Rivera NV, Ronninger M, Shchetynsky K, Franke A, Nöthen MM, Müller-Quernheim J, et al. High-density genetic mapping identifies new susceptibility variants in sarcoidosis phenotypes and shows genomic-driven phenotypic differences. Am J Respir Crit Care Med (2016) 193(9):1008–22. doi:10.1164/rccm.201507-1372OC

CrossRef Full Text | Google Scholar

8. Grunewald J, Kaiser Y, Ostadkarampour M, Rivera NV, Vezzi F, Lotstedt B, et al. T-cell receptor-HLA-DRB1 associations suggest specific antigens in pulmonary sarcoidosis. Eur Respir J (2016) 47(3):898–909. doi:10.1183/13993003.01209-2015

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Olsen HH, Grunewald J, Tornling G, Skold CM, Eklund A. Bronchoalveolar lavage results are independent of season, age, gender and collection site. PLoS One (2012) 7(8):e43644. doi:10.1371/journal.pone.0043644

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Costabel U, Hunninghake GW. ATS/ERS/WASOG statement on sarcoidosis. Sarcoidosis Statement Committee. American Thoracic Society. European Respiratory Society. World Association for Sarcoidosis and Other Granulomatous Disorders. Eur Respir J (1999) 14(4):735–7. doi:10.1034/j.1399-3003.1999.14d02.x

CrossRef Full Text | Google Scholar

11. Olerup O, Zetterquist H. HLA-DR typing by PCR amplification with sequence-specific primers (PCR-SSP) in 2 hours: an alternative to serological DR typing in clinical practice including donor-recipient matching in cadaveric transplantation. Tissue Antigens (1992) 39(5):225–35. doi:10.1111/j.1399-0039.1992.tb01940.x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Lou X, Zhang G, Herrera I, Kinach R, Ornatsky O, Baranov V, et al. Polymer-based elemental tags for sensitive bioassays. Angew Chem Int Ed Engl (2007) 46(32):6111–4. doi:10.1002/anie.200700796

CrossRef Full Text | Google Scholar

13. Bruggner RV, Bodenmiller B, Dill DL, Tibshirani RJ, Nolan GP. Automated identification of stratifying signatures in cellular subpopulations. Proc Natl Acad Sci U S A (2014) 111(26):E2770–7. doi:10.1073/pnas.1408792111

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Shekhar K, Brodin P, Davis MM, Chakraborty AK. Automatic classification of cellular expression by nonlinear stochastic embedding (ACCENSE). Proc Natl Acad Sci U S A (2014) 111(1):202–7. doi:10.1073/pnas.1321405111

CrossRef Full Text | Google Scholar

15. Becher B, Schlitzer A, Chen J, Mair F, Sumatoh HR, Teng KW, et al. High-dimensional analysis of the murine myeloid cell system. Nat Immunol (2014) 15(12):1181–9. doi:10.1038/ni.3006

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Chen H, Lau MC, Wong MT, Newell EW, Poidinger M, Chen J. Cytofkit: a bioconductor package for an integrated mass cytometry data analysis pipeline. PLoS Comput Biol (2016) 12(9):e1005112. doi:10.1371/journal.pcbi.1005112

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Wong MT, Chen J, Narayanan S, Lin W, Anicete R, Kiaang HT, et al. Mapping the diversity of follicular helper T cells in human blood and tonsils using high-dimensional mass cytometry analysis. Cell Rep (2015) 11(11):1822–33. doi:10.1016/j.celrep.2015.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol (2004) 5(10):R80. doi:10.1186/gb-2004-5-10-r80

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Huber W, Carey VJ, Gentleman R, Anders S, Carlson M, Carvalho BS, et al. Orchestrating high-throughput genomic analysis with bioconductor. Nat Methods (2015) 12(2):115–21. doi:10.1038/nmeth.3252

PubMed Abstract | CrossRef Full Text | Google Scholar

20. RStudio Team. RStudio: Integrated Development for R. (2015). Available from: http://www.rstudio.com/

Google Scholar

21. Danlos FX, Pages C, Baroudjian B, Vercellino L, Battistella M, Mimoun M, et al. Nivolumab-induced sarcoid-like granulomatous reaction in a patient with advanced melanoma. Chest (2016) 149(5):e133–6. doi:10.1016/j.chest.2015.10.082

CrossRef Full Text | Google Scholar

22. Eckert A, Schoeffler A, Dalle S, Phan A, Kiakouama L, Thomas L. Anti-CTLA4 monoclonal antibody induced sarcoidosis in a metastatic melanoma patient. Dermatology (2009) 218(1):69–70. doi:10.1159/000161122

CrossRef Full Text | Google Scholar

23. Broos CE, van Nimwegen M, In’t Veen JC, Hoogsteden HC, Hendriks RW, van den Blink B, et al. Decreased cytotoxic T-lymphocyte antigen 4 expression on regulatory T cells and Th17 cells in sarcoidosis: double trouble? Am J Respir Crit Care Med (2015) 192(6):763–5. doi:10.1164/rccm.201503-0635LE

CrossRef Full Text | Google Scholar

24. Sakthivel P, Grunewald J, Eklund A, Bruder D, Wahlstrom J. Pulmonary sarcoidosis is associated with high-level inducible co-stimulator (ICOS) expression on lung regulatory T cells – possible implications for the ICOS/ICOS-ligand axis in disease course and resolution. Clin Exp Immunol (2016) 183(2):294–306. doi:10.1111/cei.12715

CrossRef Full Text | Google Scholar

25. Idali F, Wiken M, Wahlstrom J, Mellstedt H, Eklund A, Rabbani H, et al. Reduced Th1 response in the lungs of HLA-DRB1*0301 patients with pulmonary sarcoidosis. Eur Respir J (2006) 27(3):451–9. doi:10.1183/09031936.06.00067105

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Kaiser Y, Lepzien R, Kullberg S, Eklund A, Smed-Sorensen A, Grunewald J. Expanded lung T-bet+RORgammaT+ CD4+ T-cells in sarcoidosis patients with a favourable disease phenotype. Eur Respir J (2016) 48(2):484–94. doi:10.1183/13993003.00092-2016

CrossRef Full Text | Google Scholar

27. Heron M, Grutters JC, van Moorsel CH, Ruven HJ, Huizinga TW, van der Helm-van Mil AH, et al. Variation in IL7R predisposes to sarcoid inflammation. Genes Immun (2009) 10(7):647–53. doi:10.1038/gene.2009.55

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ndhlovu LC, Leal FE, Eccles-James IG, Jha AR, Lanteri M, Norris PJ, et al. A novel human CD4+ T-cell inducer subset with potent immunostimulatory properties. Eur J Immunol (2010) 40(1):134–41. doi:10.1002/eji.200939258

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Cook L, Munier CML, Seddiki N, van Bockel D, Ontiveros N, Hardy MY, et al. Circulating gluten-specific FOXP3+CD39+ regulatory T cells have impaired suppressive function in patients with celiac disease. J Allergy Clin Immunol (2017). doi:10.1016/j.jaci.2017.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Fujimura T, Kambayashi Y, Aiba S. Expression of CD39/Entpd1 on granuloma-composing cells and induction of Foxp3-positive regulatory T cells in sarcoidosis. Clin Exp Dermatol (2013) 38(8):883–9. doi:10.1111/ced.12094

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Batista MD, Tincati C, Milush JM, Ho EL, Ndhlovu LC, York VA, et al. CD57 expression and cytokine production by T cells in lesional and unaffected skin from patients with psoriasis. PLoS One (2013) 8(2):e52144. doi:10.1371/journal.pone.0052144

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Fernandez S, French MA, Price P. Immunosenescent CD57+CD4+ T-cells accumulate and contribute to interferon-gamma responses in HIV patients responding stably to ART. Dis Markers (2011) 31(6):337–42. doi:10.1155/2011/217860

CrossRef Full Text | Google Scholar

33. Palmer BE, Blyveis N, Fontenot AP, Wilson CC. Functional and phenotypic characterization of CD57+CD4+ T cells and their association with HIV-1-induced T cell dysfunction. J Immunol (2005) 175(12):8415–23. doi:10.4049/jimmunol.175.12.8415

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Appay V, Zaunders JJ, Papagno L, Sutton J, Jaramillo A, Waters A, et al. Characterization of CD4(+) CTLs ex vivo. J Immunol (2002) 168(11):5954–8. doi:10.4049/jimmunol.168.11.5954

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ikehata A, Tomichi N. CD44v6 expression in granuloma of Crohn’s disease. Am J Gastroenterol (2000) 95(11):3325–6. doi:10.1111/j.1572-0241.2000.03324.x

CrossRef Full Text | Google Scholar

36. Linke M, Pham HT, Katholnig K, Schnoller T, Miller A, Demel F, et al. Chronic signaling via the metabolic checkpoint kinase mTORC1 induces macrophage granuloma formation and marks sarcoidosis progression. Nat Immunol (2017) 18(3):293–302. doi:10.1038/ni.3655

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Gadhoum SZ, Madhoun NY, Abuelela AF, Merzaban JS. Anti-CD44 antibodies inhibit both mTORC1 and mTORC2: a new rationale supporting CD44-induced AML differentiation therapy. Leukemia (2016) 30(12):2397–401. doi:10.1038/leu.2016.221

CrossRef Full Text | Google Scholar

38. Katchar K, Wahlstrom J, Eklund A, Grunewald J. Highly activated T-cell receptor AV2S3(+) CD4(+) lung T-cell expansions in pulmonary sarcoidosis. Am J Respir Crit Care Med (2001) 163(7):1540–5. doi:10.1164/ajrccm.163.7.2005028

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Chen Z, Han Y, Gu Y, Liu Y, Jiang Z, Zhang M, et al. CD11c(high)CD8+ regulatory T cell feedback inhibits CD4 T cell immune response via Fas ligand-Fas pathway. J Immunol (2013) 190(12):6145–54. doi:10.4049/jimmunol.1300060

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Beyer M, Wang H, Peters N, Doths S, Koerner-Rettberg C, Openshaw PJ, et al. The beta2 integrin CD11c distinguishes a subset of cytotoxic pulmonary T cells with potent antiviral effects in vitro and in vivo. Respir Res (2005) 6:70. doi:10.1186/1465-9921-6-70

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Qualai J, Li LX, Cantero J, Tarrats A, Fernandez MA, Sumoy L, et al. Expression of CD11c is associated with unconventional activated T cell subsets with high migratory potential. PLoS One (2016) 11(4):e0154253. doi:10.1371/journal.pone.0154253

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: sarcoidosis, mass cytometry, bronchoalveolar lavage, CD4+ T cells, Löfgren’s syndrome, granuloma, disease phenotypes

Citation: Kaiser Y, Lakshmikanth T, Chen Y, Mikes J, Eklund A, Brodin P, Achour A and Grunewald J (2017) Mass Cytometry Identifies Distinct Lung CD4+ T Cell Patterns in Löfgren’s Syndrome and Non-Löfgren’s Syndrome Sarcoidosis. Front. Immunol. 8:1130. doi: 10.3389/fimmu.2017.01130

Received: 12 June 2017; Accepted: 28 August 2017;
Published: 12 September 2017

Edited by:

Kjetil Taskén, University of Oslo, Norway

Reviewed by:

Xing Chang, Shanghai Institutes for Biological Sciences (CAS), China
Hu Zeng, Mayo Clinic Minnesota, United States

Copyright: © 2017 Kaiser, Lakshmikanth, Chen, Mikes, Eklund, Brodin, Achour and Grunewald. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ylva Kaiser, ylva.kaiser@ki.se