Skip to main content
  • Research article
  • Open access
  • Published:

An explorative study on deep profiling of peripheral leukocytes to identify predictors for responsiveness to anti-tumour necrosis factor alpha therapies in ankylosing spondylitis: natural killer cells in focus



Therapeutic targeting of tumour necrosis factor (TNF)-α is highly effective in ankylosing spondylitis (AS) patients. However, since one-third of anti-TNF-treated AS patients do not show an adequate clinical response there is an urgent need for new biomarkers that would aid clinicians in their decision-making to select appropriate therapeutic options. Thus, the aim of this explorative study was to identify cell-based biomarkers in peripheral blood that could be used for a pre-treatment stratification of AS patients.


A high-dimensional, multi-parametric flow cytometric approach was applied to identify baseline predictors in 31 AS patients before treatment with the TNF blockers adalimumab (TNF-neutralisation) and etanercept (soluble TNF receptor).


As the major result, the frequencies of natural killer (NK) cells, and in particular CD8-positive (CD8+) NK cell subsets, were most predictive for therapeutic outcome in AS patients. While an inverse correlation between classical CD56+/CD16+ NK cells and reduction of disease activity was observed, the CD8+ NK cell subset behaved in the opposite direction. At baseline, responders showed significantly increased frequencies of CD8+ NK cells compared with non-responders.


This is the first study demonstrating that the composition of the NK cell compartment has predictive power for prediction of therapeutic outcome for anti-TNF-α blockers, and we identified CD8+ NK cells as a potential new player in the TNF-α-driven chronic inflammatory immune response of AS.


Ankylosing spondylitis (AS) is a multifactorial chronic inflammatory rheumatic disease belonging to the group of rheumatic diseases known as spondyloarthritis (SpA) which primarily affects the axial skeleton [1]. AS has a prevalence of about 1.43 million in the European population [2] with an onset in adolescence [3] and a two-times higher occurrence in men than women [4]. The pathogenesis of AS is still obscure; it is assumed that AS is mainly caused by both genetic factors, which implies the expression of the major histocompatibility complex (MHC) I antigen human leukocyte antigen B27 (HLA-B27) [5], and also by environmental factors such as enterobacterial antigens [6]. To alleviate the axial symptoms of AS patients, non-steroidal anti-inflammatory drugs (NSAIDs) are delivered as first-line therapy. Anti-tumour necrosis factor (TNF) blocking therapy is applied only in patients with constantly high disease activity who are non-responsive to conventional NSAID treatment [7, 8]. At present, there are five anti-TNF-α agents approved for the treatment of AS: infliximab [9], a monoclonal chimeric antibody; etanercept, a soluble human TNF receptor (sTNFR)2 fusion protein [10]; adalimumab, a humanised monoclonal antibody [11]; golimumab, a fully human monoclonal antibody [12]; and certoluzimab, a Fab fragment of a humanised monoclonal antibody [13]. Most of these biologics are also successfully administered in rheumatoid arthritis (RA), psoriasis (Pso), juvenile inflammatory arthritis (JIA), and inflammatory bowel disease (IBD). The biological functions of TNF-α are mediated by binding to the membrane receptors TNFR1 (p55) or TNFR2 (p75). While TNFR1 is ubiquitously expressed in all lymphoid and myeloid immune cells and body cells [14], the expression of TNFR2 is mainly restricted to T cells [15] and natural killer (NK) cells [16]. In addition, TNFR2 can be found to be expressed in endothelial and mesenchymal cells, cells of the central nervous system (CNS), oligodendrocytes and cardiac myocytes [17], and a few other cell types [18]. According to their functions, TNFR1 is primarily associated with pro-apoptotic processes, while TNFR2 is responsible for processes ensuring survival of cells [19].

Although targeting of TNF-α is very effective in AS, around one-third of treated patients show only a poor response that can be partly attributed to the development of anti-drug antibodies (ADAb) resulting in reduced bioavailability [20]. The most likely precondition for swapping to another anti-TNF agent is partial or entire failure of effectiveness along with side effects [21]. With respect to these adverse reactions and the high costs of anti-TNF agents leading to high economic burden for the health care systems, it is desirable to stratify patients according to treatment predictors prior to biological therapy. Various demographic and clinical parameters such as high baseline disease activity, short disease duration, young age, male sex, and presence of HLA-B27 have been shown to correlate with adequate clinical short- and long-term response [22,23,24,25,26]. In addition, modern imaging techniques, such as magnetic resonance imaging (MRI), are used to correlate bone and tissue destruction with treatment response [27]. However, these techniques are time consuming and expensive when used as a standard pre-treatment assessment. Another level of therapy response prediction was investigated when responsiveness to anti-TNF agents was related to the presence of different TNF-α genotypes. It was reported that patients with TNF-α–308 G/G, –857 C/C, or –1031 T/T genotypes showed a better response to anti-TNF agents than patients without these polymorphisms [28, 29]. Apart from these observations, there are no reliable predictive biomarkers for anti-TNF responsiveness in AS. Using a multi-parametric flow cytometric approach, we aimed to identify cell-based biomarkers in the peripheral blood of AS patients that are able to predict a successful therapeutic response to TNF inhibitors before starting therapy. As a result, we found that a low pre-treatment frequency of a CD8-expressing subpopulation of NK cells is associated with a lack of therapeutic response.



Ethics statement

The study was performed in accordance with the 1964 Declaration of Helsinki and approved by the Charité University Medicine ethics committee I of Charité Campus Mitte. All patients provided written informed consent to participate in the study. Furthermore, we declare that this manuscript contains no information or images that could lead to identification of a study participant.

A total of 31 AS patients (22 male, 9 female) of whom 81% were positive for HLA-B27, recruited from the rheumatology outpatient clinics of the Charité, and 10 healthy controls (HC; 7 male, 3 female) participated in the study. The patients had an average age of 38 ± 10.2 years and the HC averaged 34 ± 10.7 years. All patients fulfilling the modified New York criteria [30] and who were eligible for anti-TNF inhibitor treatment because of persistently high disease activity (Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) > 4) despite treatment with NSAIDs or who were unable to take NSAIDs due to contraindications were included in the study (Table 1). Disease activity was assessed according to the BASDAI index consisting of a score range from 0 (no symptoms) to 10 (high disease activity). The mean baseline BASDAI prior to TNF inhibitor therapy was 6.2 ± 1.3 (Table 1).

Table 1 Demographic and disease characteristics of AS patients treated with ETN and ADA, respectively

Prior to the start of TNF inhibitor therapy, 10 ml heparinised blood was taken to perform flow cytometric analysis. Fifteen patients were treated with etanercept (Enbrel; Amgen, and Pfizer) and 16 patients with adalimumab (Humira; AbbVie Inc.). The BASDAI score was obtained at baseline and at follow-up visits [31]. The response to treatment was assessed between 1 and 6 months after the start of therapy and defined as a 50% BASDAI reduction (BASDAI50 response) relative to baseline BASDAI (Additional file 1: Table S1).

Blood sample preparation, antibody staining, and flow cytometry measurement

Blood sample preparation and antibody staining procedures were as described previously [32]. Cells obtained from the blood of patients prior to treatment were stained for 50 different surface antigens in a seven-colour staining combined to 10 tubes (Table 2). After staining, cells were fixed with 1% paraformaldehyde and analysed within 24 h. We did not include a live/dead cell staining, but cell debris, erythrocytes, and thrombocytes were excluded according to their SSC/FSC characteristics.

Table 2 Staining matrix showing antibodies and their corresponding fluorochrome conjugates measured in ten separate staining tubes

Data acquisition was accomplished with a FACSCanto™ II Flow Cytometer (BD Biosciences, USA) with an average cell count of one million cells per sample. To warrant reproducibility and to survey the instruments’ performance, a BD™ Cytometer Setup and Tracking Beads were regularly used before each measurement. In addition, we have always monitored the quality of antibody staining directly after data acquisition by monitoring each individual fluorescence channel used for each particular staining tube. For this 20,000 randomly selected events were plotted. Samples which did not pass this quality check were excluded from further analysis.

Data analysis and statistical analysis

Two different software tools were applied to analyse the complex datasets generated by this unsupervised flow cytometry approach, which is based on both manual and automatic bioinformatic strategies identify potential candidate phenotypes. In the first approach, the relevant statistics such as mean fluorescence intensities (MFIs) and absolute cell numbers of manually analysed data were transferred as a comma-separated value (CSV) file format to an Access database as shown previously [32]. This primary data analysis including colour compensation and gate setting was performed by FACSDIVA v6.0 software (BD Biosciences, USA). The second approach utilised the automated classification algorithm immunoClust, which processes uncompensated raw data and therefore excludes any operator-dependent gating or compensation artefacts [33]. Population clustering and comparative meta-clustering of immunoClust assume finite mixture models and use Expectation Maximisation (EM)-iterations with an integrated classification likelihood (ICL) criterion to stabilise the number of reasonable clusters. For meta-clustering, a probability measure on Gaussian distributions was invented, which is based on the Bhattacharyya Coefficients. Meta-clusters were manually annotated and classified.

Linear regression and receiver operating characteristic (ROC) analysis performed with Prism 5 (GraphPad Software, Inc.) was used to elucidate associations between candidate markers and clinical parameters. For statistical data analysis, the Welch corrected t test was used where p values < 0.05 were determined as statistically significant.


Patient baseline characteristics and their clinical responses

The study design encompassed 31 AS patients with high disease activity indicated by a baseline BASDAI of 6.2 ± 1.3 before treatment with adalimumab (ADA; n = 16) or etanercept (ETN; n = 15). The patient demographic and baseline clinical characteristics are summarised in Table 1 and showed no significant differences between the ADA- and ETN-treated patients.

The average BASDAI assessment date was 3.1 ± 1.4 month after the start of treatment. After 1 to 6 months of treatment, the relative reduction of disease activity assessed by BASDAI was 51.1 ± 31.7% for ETN-treated patients and 44.3 ± 31.6% for patients who received ADA. According to the BASDAI response criteria, five patients in the ETN group and seven patients in the ADA group failed to respond (Additional file 1: Table S1).

Clinical parameters available at baseline, such as disease duration, baseline BASDAI, C-reactive peptide (CRP) levels, and erythrocyte sedimentation rate (ESR), as well as the expression of HLA-B27, allowed no discrimination between patients who would respond to TNF-α blockers and those who would fail.

For data analysis, we have applied an unsupervised, automated cell clustering approach, immunoClust [33], to identify potential immunophenotypic parameters that enable classification of AS patients into responders (R) and non-responders (NR) prior to anti-TNF treatment. Since the freshly obtained patient blood was immediately processed we have not performed dead cell staining, but we gated for live cells according to cell size and granularity defined by forward scatter (FSC) and sideward scatter (SSC) characteristics, respectively.

The presented two-dimensional clustering approach with patients clustered in the vertical direction and immunophenotypic parameters clustered in the horizontal direction gives an overview of all leukocyte subsets including certain activation markers that are differentially expressed in R and NR. For the identification of significant parameters by the immunoClust algorithm, FCS files (Flow Cytometry Standard file format) of uncompensated raw data were used and finally disclosed 36 parameters when all patient samples were considered. Analysing ETN- and ADA-treated patients separately revealed 21 and 27 parameters, respectively (Fig. 1).

Fig. 1
figure 1

Two-dimensional hierarchical cluster analyses of immunophenotypic parameters that are differentially expressed in R and NR prior to treatment when considering all patients (a), etanercept (ETN) alone (b), and adalimumab (ADA) alone (c). Columns represent individual patients and the colour code below the vertical dendrogram indicates responders (Resp.; blue) or non-responders (Non R.; red) and the TNF-blocker used (either ETN (cyan) or ADA (pink)). Rows of the cluster diagram represent immunophenotypic parameters with p values < 0.1. The magnitude of parameter expression is colour coded with red for a relatively increased and blue for a relatively decreased expression. The colour code for the horizontal dendrogram indicates the expression in a particular cell type, such as natural killer (NK) cells (cyan), B cells (green), T cells (raspberry-red), monocytes (mo; orange), granulocytes (gr; blue), and CD3-negative lymphocytes (ly CD3–; white). In total, one million cells were acquired per sample to ensure that even rare cell populations with frequencies around 0.1% could be reliable detected

Although using all these parameters did not allow an error-free classification of R and NR, all samples were grouped into two main clusters which were enriched for R and NR, respectively (Fig. 1a). Surprisingly, more than 50% of the discriminating parameters could be clearly assigned to NK cell subsets if all patients were analysed together (Fig. 1a). For further analysis of NK cell-related subsets, and knowing that ETN and ADA have different modes of action to neutralise the effect of TNF-α, we continued to investigate both treatment groups separately to identify therapy-specific response signatures. Using this approach, the majority of parameters that significantly discriminate between NR and R in the ETN group (Fig. 1b) and ADA group (Fig. 1c) were related to the NK cell compartment. The best classification of R and NR was achieved in the group of ETN-treated patients. Here, only two of 10 R were grouped as NR and one single of five NR was grouped as R.

Validation of classical NK cells and CD8-positive NK cells as potential immunological biomarkers for an anti-TNF-α therapy prediction

Since both the percentage distribution of classical NK cells in general and the ratio of CD8-positive and CD8-negative NK cells in particular appeared to be the most promising predictors for an anti-TNF therapy response, we looked at the CD8 receptor expression on NK cells in more detail. In Fig. 2 the general gating strategy to define classical CD56dimCD16+ NK cells is shown. In Additional file 2 (Figure S1), backgating for monocytes and NK cells demonstrates that, despite labelling multiple antigens by the same fluorochrome (CD14/CD56/IgD labelled to PE), unravelling of these complex stainings with respect to NK cells, monocytes, T cells, and B cells is possible.

Fig. 2
figure 2

Representative gating strategy to identify CD8-positive and CD8-negative NK cell subsets. To increase parameter diversity detectable in a single seven-colour staining setup we have combined specifically expressed cell lineage markers that are labelled by the same fluorochrome, such as CD14-PE/CD56-PE/IgD-PE. Its cell-specific expression can be deconvoluted by sequential gating as exemplified by a: (left panel) gating on small cells (G1) with low granularity according to FSC/SSC characteristics allows exclusion of monocytes (CD14) and granulocytes (CD16); (middle panel) exclusion of B lymphocytes by CD19 and T lymphocytes by CD3 (G2); and (right panel) CD19/CD3-double negative cells were analysed for CD56 versus CD16 allowing us to identify NK cells by the co-expression of CD56 and CD16 (G3). b Classical NK cells were further analysed according to their expression of the CD8 receptor in five exemplary responders (R; top row) and non-responders (NR; bottom row)

At first, a lymphocyte scatter gate (G1) was set (Fig. 2a). Next T and B cells were excluded out of the lymphocyte population (G2) to determine the percentage of CD56dimCD16+ NK cells (G3). Subsequently, the CD8-positive NK cells were quantified by sub-gating as shown exemplarily for five responders (R2, R9, R16, R12, and R18) and five non-responders (NR1, NR7, NR8, NR10, and NR12) (Fig. 2b).

Figure 3 shows the frequencies of CD8-positive NK cells for healthy controls (n = 10), anti-TNF R (n = 19), and anti-TNF NR (n = 12). A significantly higher frequency of CD8 receptor-bearing NK cells was observed in the R group compared with NR (Fig. 3). Comparing the percentage of CD8-positive NK cells to HC, AS patients who will not respond showed significantly lower frequencies of CD8-expressing NK cells (Fig. 3). Individual values for frequencies and absolute counts of NK cells and CD8-positive subsets are given in Additional file 1 (Table S1).

Fig. 3
figure 3

Frequencies of CD8-positive natural killer (NK) cells in responders (R) and non-responders (NR) assessed before treatment for all patients or for adalimumab (ADA) and etanercept (ETN) separately. Significance was determined by Welch’s corrected t test. HC healthy controls

Next, we performed Spearman’s rank correlation and linear regression analyses which showed a significant inverse correlation of frequencies of classical NK cells and an improved therapy outcome if all 31 patients were included (Fig. 4a; p = 0.01, r2 = 0.19).

Fig. 4
figure 4

Linear regression and Spearman’s rank correlation analyses showing relative reduction of Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) values within 1 to 6 months after treatment with ADA and ETN (n = 31), associated with changes in frequencies of classical natural killer (NK) cells (a) and CD8-positive NK cells (b). ROC curves for all patients demonstrating the value of classical (c) and CD8-expressing NK cells (d) as prognostic markers according to increasing area under the curve (AUC) values

Surprisingly, a positive, linear correlation was found when we focused on the analysis of CD8- expressing NK cells (Fig. 4b). Here, the frequency of CD8-positive cells related to the cells of the total NK cells (G3) clearly correlated with a successful therapy response. This correlation was statistically significant when all patients were considered (p = 0.002, r2 = 0.29).

Furthermore, reverse regression analyses, displayed by ROC curves, were performed to verify the quality of classical and especially CD8-positive NK cells as suitable cellular biomarkers for predicting an anti-TNF therapy outcome (Fig. 4c, d). If all samples were included for this analysis, the frequency of classical NK cells was the most promising parameter as a baseline predictor (Fig. 4c; p = 0.008, area under the curve (AUC) = 0.79). Slightly improved values were obtained if the frequency of CD8-positive NK cells was used (Fig. 4d; p = 0.007, AUC = 0.79). Therefore, our data implicate that the appearance of CD8-positive NK cells is a robust biomarker for the prediction of an anti-TNF-α response at baseline.


To our knowledge, this is the first study aiming to identify cellular biomarkers in peripheral blood to stratify AS patients upfront with regard to subsequent responsiveness to anti-TNF-α therapy. To date, there are no quantifiable laboratory parameters available that can be used for a personalised response prediction [34]. Although higher values of CRP, ESR, and serum amyloid A (SAA), together with the presence of HLA-B27, have been reported to be useful baseline predictors for a successful anti-TNF-α therapy response in AS, the robustness, sensitivity, and specificity fail if applied to individual patients [35, 36]. Even our data, such as the age of patients, disease duration, CRP levels, or blood sedimentation, did not allow any prediction of future responsiveness.

Other potential new biomarkers described, such as single nucleotide polymorphisms (SNPs) and activity of endoplasmic reticulum aminopeptidase (ERAP)-1 [37], serum levels of matrix metalloproteinase (MMP)-3 [38], or vascular endothelial growth factor (VEGF) [39], have not yet reached a standard in the clinical diagnostic routine of AS. Generally, it is challenging to correlate quantifiable parameters to the disease activity BASDAI score, which is calculated on the basis of a patient’s subjective assessment of well-being and therefore is only of limited value when used as an absolute variable.

In our explorative study presented here we used an integrated multi-parametric flow cytometry and a new unsupervised data clustering approach to identify possible cellular biomarkers that are qualitatively or quantitatively different in responders and non-responders. A similar approach has been successfully used to classify active AS patients from healthy donors according to specific phenotypical changes in blood [32]. A comprehensive overview about immunophenotypical changes described so far in different autoimmune diseases is given by Alegria et al. [40].

Principally, at first view, an increased frequency of the major CD56dimCD16+ NK cell subset was indicative for a weak therapeutic anti-TNF response. This finding is in line with other reports showing that innate immunity and particularly NK cells may play a central role in the pathogenesis of various autoimmune diseases both in a protective and pathogenic manner. Their frequency and functionality were investigated in chronic inflammation, such as rheumatoid arthritis [41], multiple sclerosis [42], psoriasis [43], and systemic lupus erythematosus [44]. Despite some conflicting results, an overall decrease and a cytolytic impairment of circulating NK cells could be observed. In AS, the functionality of NK cells is described to be likewise impaired [45] but, in contrast to other autoimmune diseases, an increased NK cell number is reported [46, 47]. We could confirm this observation if comparing the non-responder group with normal donors (Fig. 3a); however, the NK cell frequency in the responder cohort was similar to the number of healthy individuals. The strong association between NK cells and the aetiology of AS is underlined by the capability of the disease-dependent elevated killer cell immunoglobulin-like receptors (KIRs) to recognise the HLA-B27 antigen, which is expressed in 80–90% of AS patients [48, 49]. Concomitant with a pathogenetic upregulated frequency of KIR3DL1+ NK cells in AS, interferon (IFN)γ production is correspondingly diminished [50]. Moreover, activated KIR3DL2+ NK cells are increased in SpA and may play a pathogenic role.

An in-depth analysis of the NK cell compartment revealed that the frequency of classical NK cells expressing the CD8 antigen showed a significant and positive correlation with anti-TNF responsiveness. It is known that about 40% of NK cells variably express CD8 in an α/α homodimeric form whereas the CD8-positive subset is described to exhibit enhanced cytotoxic features as compared with its CD8-negative counterpart [51, 52]. We could not detect significant age- or sex-related differences with respect to the frequency of CD8-positive NK cells either in the group of healthy controls or in that of AS patients. Thus, our data imply that NK cells expressing the CD8αα homodimer are directly or indirectly involved in the immunosuppressive effect exerted by anti-TNF-α blockers.

If CD8αα-expressing NK cells are directly involved it can be assumed that their increased cytotoxic activity and/or their diminished behaviour regarding cytotoxicity-induced apoptosis are responsible for the improved responsiveness of TNF blockers. Alternatively, it can be hypothesised that an engagement of CD8αα receptors on NK cells to molecules of the HLA-I family can cause an increased secretion of the pro-inflammatory cytokines TNF-α and IFNγ [53], which in turn promote both TNF receptor (TNFR) synthesis and its proteolysis to a soluble form [54]. By this mechanism, endogenously synthesised sTNFRs [55, 56] can co-operatively support the action of therapeutic TNF inhibitors [57, 58]. These non-signalling ‘decoy’ receptors are still competent for binding TNF and thus may function as a natural TNF antagonist comparable to ETN [59]. Since we have included adalimumab- and etanercept-treated patients in our study, it was interesting to know if differences in the prediction of responsiveness were detectable when either TNF-α was neutralised by ADA or scavenged by the sTNFR2 fusion protein (ETN). Although group size reduction was responsible for a less statistical power of prediction analysis, we could ascertain a better correlation with respect to ETN-treated patients compared with ADA-treated patients. Therefore, these findings are encouraging for validation in new independent cohorts of appropriate group sizes. If this result could be validated it would indicate diverse and more complex modes of action going beyond the mere neutralising effect of TNF-α blockers.

To elucidate possible differences in the TNFR expression in CD8-positive and CD8-negative NK cell subsets isolated from healthy individuals, we performed global gene expression analyses but could not detect any differential expression of TNFR1, TNFR2, or TNF-α (data not shown). Comparing the expression of TNFR1 and TNFR2 in healthy individuals, it was obvious that higher expression levels were detectable for TNFR2 (data not shown). Unfortunately, thus far we have had no opportunity for analysing cells isolated from AS patients to test if TNFRs or TNF-α were differentially expressed under chronic inflammatory conditions. Thus, according to our data, it can be postulated that CD8-positive NK cells are obviously capable of amplifying the neutralising effect of TNF-blockers, but it remains unclear why this amplification is preferentially observed in ETN- and not in ADA-treated patients.


Although our findings are promising, a further validation of CD8-positive NK cells as a potential biomarker for TNF responsiveness is necessary in an independent cohort of AS and other rheumatic and gastrointestinal diseases where anti-TNF blockers are successfully administered. Nevertheless, our study is a first proof of the concept that cellular response signatures can be identified in peripheral blood by an extensive immunophenotyping approach. Even though the exact mechanism of how CD8-positive NK cells and the therapeutic effect of TNF-α blockers are interrelated in AS is only poorly understood thus far, monitoring these cells by flow cytometry offers an interesting new diagnostic option with respect to the challenges of an individualised therapy concept, at least in the field of chronic inflammatory rheumatic diseases.





Ankylosing spondylitis


Bath Ankylosing Spondylitis Disease Activity Index


C-reactive protein


Endoplasmic reticulum aminopeptidase


Erythrocyte sedimentation rate




Healthy controls


Inflammatory bowel disease


Juvenile inflammatory arthritis


Killer cell immunoglobulin-like receptor


Mean fluorescence intensity


Matrix metalloproteinase


Natural killer




Non-steroidal anti-inflammatory drug






Rheumatoid arthritis


Serum amyloid A


Single nucleotide polymorphism




Soluble human tumour necrosis factor receptor


Tumour necrosis factor


Vascular endothelial growth factor


  1. Braun J, Sieper J. Ankylosing spondylitis. Lancet. 2007;369:1379–90.

    Article  PubMed  Google Scholar 

  2. Dean LE, et al. Global prevalence of ankylosing spondylitis. Rheumatology (Oxford). 2014;53:650–7.

    Article  Google Scholar 

  3. Feldtkeller E, Khan M, van der Heijde D, van der Linden S, Braun J. Age at disease onset and diagnosis delay in HLA-B27 negative vs. positive patients with ankylosing spondylitis. Rheumatol Int. 2003;23:61–6.

    PubMed  Google Scholar 

  4. Hill HF, Hill AG, Bodmer JG. Clinical diagnosis of ankylosing spondylitis in women and relation to presence of HLA-B27. Ann Rheum Dis. 1976;35:267–70.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  5. Brewerton DA, et al. Ankylosing spondylitis and HL-A 27. Lancet. 1973;1:904–7.

    Article  PubMed  CAS  Google Scholar 

  6. Trull A, Ebringer A, Panayi G, Ebringer R, James DC. HLA-B27 and the immune response to enterobacterial antigens in ankylosing spondylitis. Clin Exp Immunol. 1984;55:74–80.

    PubMed  PubMed Central  CAS  Google Scholar 

  7. Coates LC, Marzo-Ortega H, Bennett AN, Emery P. Anti-TNF therapy in ankylosing spondylitis: insights for the clinician. Ther Adv Musculoskelet Dis. 2010;2:37–43.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Sieper J, Braun J, Rudwaleit M, Boonen A, Zink A. Ankylosing spondylitis: an overview. Ann Rheum Dis. 2002;61(Suppl 3):iii8–18.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Braun J, et al. Efficacy and safety of infliximab in patients with ankylosing spondylitis over a two-year period. Arthritis Rheum. 2008;59:1270–8.

    Article  PubMed  Google Scholar 

  10. Dijkmans B, et al. Etanercept in the long-term treatment of patients with ankylosing spondylitis. J Rheumatol. 2009;36:1256–64.

    Article  PubMed  CAS  Google Scholar 

  11. van der Heijde D, et al. Adalimumab effectiveness for the treatment of ankylosing spondylitis is maintained for up to 2 years: long-term results from the ATLAS trial. Ann Rheum Dis. 2009;68:922–9.

    Article  PubMed  CAS  Google Scholar 

  12. Deodhar A, et al. Golimumab administered subcutaneously every 4 weeks in ankylosing spondylitis: 5-year results of the GO-RAISE study. Ann Rheum Dis. 2015;74:757–61.

    Article  PubMed  CAS  Google Scholar 

  13. Landewe R, et al. Efficacy of certolizumab pegol on signs and symptoms of axial spondyloarthritis including ankylosing spondylitis: 24-week results of a double-blind randomised placebo-controlled phase 3 study. Ann Rheum Dis. 2014;73:39–47.

    Article  PubMed  CAS  Google Scholar 

  14. Aggarwal BB, Gupta SC, Kim JH. Historical perspectives on tumor necrosis factor and its superfamily: 25 years later, a golden journey. Blood. 2012;119:651–65.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Ware C, et al. Tumor necrosis factor (TNF) receptor expression in T lymphocytes. Differential regulation of the type I TNF receptor during activation of resting and effector T cells. J Immunol. 1991;147:4229–38.

    PubMed  CAS  Google Scholar 

  16. Mason A, et al. Regulation of NK cells through the 80-kDa TNFR (CD120b). J Leukoc Biol. 1995;58:249–55.

    Article  PubMed  CAS  Google Scholar 

  17. Yang L, Lindholm K, Konishi Y, Li R, Shen Y. Target depletion of distinct tumor necrosis factor receptor subtypes reveals hippocampal neuron death and survival through different signal transduction pathways. J Neurosci. 2002;22:3025–32.

    Article  PubMed  CAS  Google Scholar 

  18. Faustman DL, Davis M. TNF receptor 2 and disease: autoimmunity and regenerative medicine. Front Immunol. 2013;4(478):1–8.

    CAS  Google Scholar 

  19. Naude PJ, den Boer JA, Luiten PG, Eisel UL. Tumor necrosis factor receptor cross-talk. FEBS J. 2011;278:888–98.

    Article  PubMed  CAS  Google Scholar 

  20. Vincent FB, et al. Antidrug antibodies (ADAb) to tumour necrosis factor (TNF)-specific neutralising agents in chronic inflammatory diseases: a real issue, a clinical perspective. Ann Rheum Dis. 2013;72:165–78.

    Article  PubMed  CAS  Google Scholar 

  21. Haberhauer G, Strehblow C, Fasching P. Observational study of switching anti-TNF agents in ankylosing spondylitis and psoriatic arthritis versus rheumatoid arthritis. Wien Med Wochenschr. 2010;160:220–4.

    Article  PubMed  Google Scholar 

  22. Arends S, et al. Baseline predictors of response and discontinuation of tumor necrosis factor-alpha blocking therapy in ankylosing spondylitis: a prospective longitudinal observational cohort study. Arthritis Res Ther. 2011;13:R94.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  23. Arends S, van der Veer E, Kallenberg CG, Brouwer E, Spoorenberg A. Baseline predictors of response to TNF-alpha blocking therapy in ankylosing spondylitis. Curr Opin Rheumatol. 2012;24:290–8.

    Article  PubMed  CAS  Google Scholar 

  24. Lord PA, et al. Predictors of response to anti-TNF therapy in ankylosing spondylitis: results from the British Society for Rheumatology biologics register. Rheumatology (Oxford). 2010;49:563–70.

    Article  CAS  Google Scholar 

  25. Maneiro JR, Souto A, Salgado E, Mera A, Gomez-Reino JJ. Predictors of response to TNF antagonists in patients with ankylosing spondylitis and psoriatic arthritis: systematic review and meta-analysis. RMD Open. 2015;1:e000017.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Rudwaleit M, Listing J, Brandt J, Braun J, Sieper J. Prediction of a major clinical response (BASDAI 50) to tumour necrosis factor alpha blockers in ankylosing spondylitis. Ann Rheum Dis. 2004;63:665–70.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  27. Rudwaleit M, et al. MRI in predicting a major clinical response to anti-tumour necrosis factor treatment in ankylosing spondylitis. Ann Rheum Dis. 2008;67:1276–81.

    Article  PubMed  CAS  Google Scholar 

  28. Seitz M, Wirthmuller U, Moller B, Villiger PM. The −308 tumour necrosis factor-alpha gene polymorphism predicts therapeutic response to TNFalpha-blockers in rheumatoid arthritis and spondyloarthritis patients. Rheumatology (Oxford). 2007;46:93–6.

    Article  CAS  Google Scholar 

  29. Tong Q, et al. TNF-alpha −857 and −1031 polymorphisms predict good therapeutic response to TNF-alpha blockers in Chinese Han patients with ankylosing spondylitis. Pharmacogenomics. 2012;13:1459–67.

    Article  PubMed  CAS  Google Scholar 

  30. van der Linden S, Valkenburg HA, Cats A. Evaluation of diagnostic criteria for ankylosing spondylitis. A proposal for modification of the New York criteria. Arthritis Rheum. 1984;27:361–8.

    Article  PubMed  CAS  Google Scholar 

  31. Garrett S, et al. A new approach to defining disease status in ankylosing spondylitis: the Bath Ankylosing Spondylitis Disease Activity Index. J Rheumatol. 1994;21:2286–91.

    PubMed  CAS  Google Scholar 

  32. Steinbrich-Zollner M, et al. From transcriptome to cytome: integrating cytometric profiling, multivariate cluster, and prediction analyses for a phenotypical classification of inflammatory diseases. Cytometry A. 2008;73:333–40.

    Article  PubMed  Google Scholar 

  33. Sörensen T, Baumgart S, Durek P, Grützkau A, Häupl T. ImmunoClust—an automated analysis pipeline for the identification of immunophenotypic signatures in high-dimensional cytometric datasets. Cytometry A. 2015;87:603–15.

    Article  PubMed  Google Scholar 

  34. Reveille J. Biomarkers for diagnosis, monitoring of progression, and treatment responses in ankylosing spondylitis and axial spondyloarthritis. Clin Rheumatol. 2015;34:1009–18.

    Article  PubMed  PubMed Central  Google Scholar 

  35. de Vries MK, et al. Erythrocyte sedimentation rate, C-reactive protein level, and serum amyloid a protein for patient selection and monitoring of anti-tumor necrosis factor treatment in ankylosing spondylitis. Arthritis Rheum. 2009;61:1484–90.

    Article  PubMed  CAS  Google Scholar 

  36. Luc M, et al. C-reactive protein predicts tumor necrosis factor-alpha blocker retention rate in axial ankylosing spondylitis. J Rheumatol. 2007;34:2078–81.

    PubMed  CAS  Google Scholar 

  37. Evans D, et al. Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility. Nat Genet. 2011;43:761–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Arends S, et al. Serum MMP-3 level as a biomarker for monitoring and predicting response to etanercept treatment in ankylosing spondylitis. J Rheumatol. 2011;38:1644–50.

    Article  PubMed  CAS  Google Scholar 

  39. Poddubnyy D, et al. Elevated serum level of the vascular endothelial growth factor predicts radiographic spinal progression in patients with axial spondyloarthritis. Ann Rheum Dis. 2014;73:2137–43.

    Article  PubMed  CAS  Google Scholar 

  40. Carvajal Alegria G, Gazeau P, Hillion S, Daien CI, Cornec DYK. Could lymphocyte profiling be useful to diagnose systemic autoimmune diseases? Clinic Rev Allerg Immunol. 2017;53:219–36.

    Article  CAS  Google Scholar 

  41. Aramaki T, et al. A significantly impaired natural killer cell activity due to a low activity on a per-cell basis in rheumatoid arthritis. Mod Rheumatol. 2009;19:245–52.

    Article  PubMed  Google Scholar 

  42. Vranes Z, Poljakovic Z, Marusic M. Natural killer cell number and activity in multiple sclerosis. J Neurol Sci. 1989;94:115–23.

    Article  PubMed  CAS  Google Scholar 

  43. Cameron AL, Kirby B, Griffiths CE. Circulating natural killer cells in psoriasis. Br J Dermatol. 2003;149:160–4.

    Article  PubMed  CAS  Google Scholar 

  44. Ewan PW, Barrett HM, Pusey CD. Defective natural killer (NK) and killer (K) cell function in systemic lupus erythematosus. J Clin Lab Immunol. 1983;10:71–6.

    PubMed  CAS  Google Scholar 

  45. Conigliaro P, Scrivo R, Valesini G, Perricone R. Emerging role for NK cells in the pathogenesis of inflammatory arthropathies. Autoimmun Rev. 2011;10:577–81.

    Article  PubMed  CAS  Google Scholar 

  46. Lopez-Botet M, Moretta L, Strominger J. NK-cell receptors and recognition of MHC class I molecules. Immunol Today. 1996;17:212–4.

    Article  PubMed  CAS  Google Scholar 

  47. Mousavi T, et al. Phenotypic study of natural killer cell subsets in ankylosing spondylitis patients. Iran J Allergy Asthma Immunol. 2009;8:193–8.

    PubMed  CAS  Google Scholar 

  48. Azuz-Lieberman N, et al. The involvement of NK cells in ankylosing spondylitis. Int Immunol. 2005;17:837–45.

    Article  PubMed  CAS  Google Scholar 

  49. Chan AT, Kollnberger SD, Wedderburn LR, Bowness P. Expansion and enhanced survival of natural killer cells expressing the killer immunoglobulin-like receptor KIR3DL2 in spondylarthritis. Arthritis Rheum. 2005;52:3586–95.

    Article  PubMed  CAS  Google Scholar 

  50. Scrivo R, Morrone S, Spadaro A, Santoni A, Valesini G. Evaluation of degranulation and cytokine production in natural killer cells from spondyloarthritis patients at single-cell level. Cytometry B Clin Cytom. 2011;80:22–7.

    Article  PubMed  Google Scholar 

  51. Addison EG, et al. Ligation of CD8alpha on human natural killer cells prevents activation-induced apoptosis and enhances cytolytic activity. Immunology. 2005;116:354–61.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  52. Srour EF, Leemhuis T, Jenski L, Redmond R, Jansen J. Cytolytic activity of human natural killer cell subpopulations isolated by four-color immunofluorescence flow cytometric cell sorting. Cytometry. 1990;11:442–6.

    Article  PubMed  CAS  Google Scholar 

  53. Gibbings D, Befus AD. CD4 and CD8: an inside-out coreceptor model for innate immune cells. J Leukoc Biol. 2009;86:251–9.

    Article  PubMed  CAS  Google Scholar 

  54. Fischer R, Kontermann RE, Maier O. Targeting sTNF/TNFR1 signaling as a new therapeutic strategy. Antibodies. 2015;4:48–70.

    Article  CAS  Google Scholar 

  55. Solomaon KA, Pesti N, Wu G, Newton RC. Cutting edge: a dominant negative form of TNF-alpha converting enzyme inhibits proTNF and TNFRII secretion. J Immunol. 1999;163:4105–8.

    Google Scholar 

  56. Wallach D, et al. Soluble and cell surface receptors for tumor necrosis factor. Agents Actions Suppl. 1991;35:51–7.

    PubMed  CAS  Google Scholar 

  57. Meusch U, et al. In vitro response pattern of monocytes after tmTNF reverse signaling predicts response to anti-TNF therapy in rheumatoid arthritis. J Transl Med. 2015;13:256.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Pinckard JK, Sheehan KC, Arthur CD, Schreiber RD. Constitutive shedding of both p55 and p75 murine TNF receptors in vivo. J Immunol. 1997;158:3869–73.

    PubMed  CAS  Google Scholar 

  59. Higuchi M, Aggarwal BB. TNF induces internalization of the p60 receptor and shedding of the p80 receptor. J Immunol. 1994;152:3550–8.

    PubMed  CAS  Google Scholar 

Download references


This study was funded by the following grants: IMI JU-funded Project BeTheCure (grant no. 115142–2), the German Research Foundation (DFG) (grant no. SFB 650 TPZ06), the German Federal Ministry of Education and Research (BMBF) within the framework of the e:Med research and funding concept sysINFLAME (grant no. 01ZX1306B), and Leibniz ScienceCampus Chronic Inflammation ( These funders had no involvement in the study design, in the collection, analysis, and interpretation of the data, in the writing of the report, or in the decision to submit the paper for publication.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations



US-W performed experiments, analysed data, and wrote the manuscript. TS, JRG, and TH analysed data. MS-Z, HH, and PW performed experiments. AR and JS designed the experimental strategy. US and DP characterised the patients and determined clinical inclusion criteria. AG designed the experimental strategy, analysed data, and wrote the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Andreas Grützkau.

Ethics declarations

Ethics approval and consent to participate

The study was performed in accordance with the 1964 Declaration of Helsinki and approved by the “Charité University Medicine ethics committee I of Charité Campus Mitte”, and all patients provided written informed consent to participate in the study. Furthermore, we declare that this manuscript contains no information or images that could lead to identification of a study participant.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Individual data on treatment, responsiveness, NK cell counts per μl blood, and frequencies related to total leukocytes. Moreover, absolute numbers and frequencies of CD8-positive NK cells in relation to all CD56/CD16-double positive NK cells are given. n.d. not determined. (DOCX 25 kb)

Additional file 2:

Figure S1. Clarification of the gating strategy for NK cells shown in Fig. 2. By backgating, monocytes can be clearly separated from NK cells (Gate A, PBMC; Gate B, CD3/CD19-double negative PBMC; Gate C, monocytes; and Gate D, NK cells). (DOC 87 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schulte-Wrede, U., Sörensen, T., Grün, J.R. et al. An explorative study on deep profiling of peripheral leukocytes to identify predictors for responsiveness to anti-tumour necrosis factor alpha therapies in ankylosing spondylitis: natural killer cells in focus. Arthritis Res Ther 20, 191 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: