Skip to main content

Association of cytokine and matrix metalloproteinase profiles with disease activity and function in ankylosing spondylitis



The pathology of ankylosing spondylitis (AS) suggests that certain cytokines and matrix metalloproteinases (MMPs) might provide useful markers of disease activity. Serum levels of some cytokines and MMPs have been found to be elevated in active disease, but there is a general lack of information about biomarker profiles in AS and how these are related to disease activity and function. The purpose of this study was to investigate whether clinical measures of disease activity and function in AS are associated with particular profiles of circulating cytokines and MMPs.


Measurement of 30 cytokines, five MMPs and four tissue inhibitors of metalloproteinases was carried out using Luminex® technology on a well-characterised population of AS patients (n = 157). The relationship between biomarker levels and measures of disease activity (Bath ankylosing spondylitis disease activity index (BASDAI)), function (Bath ankylosing spondylitis functional index) and global health (Bath ankylosing spondylitis global health) was investigated. Principal component analysis was used to reduce the large number of biomarkers to a smaller set of independent components, which were investigated for their association with clinical measures. Further analyses were carried out using hierarchical clustering, multiple regression or multivariate logistic regression.


Principal component analysis identified eight clusters consisting of various combinations of cytokines and MMPs. The strongest association with the BASDAI was found with a component consisting of MMP-8, MMP-9, hepatocyte growth factor and CXCL8, and was independent of C-reactive protein levels. This component was also associated with current smoking. Hierarchical clustering revealed two distinct patient clusters that could be separated on the basis of MMP levels. The high MMP cluster was associated with increased C-reactive protein, the BASDAI and the Bath ankylosing spondylitis functional index.


A profile consisting of high levels of MMP-8, MMP-9, hepatocyte growth factor and CXCL8 is associated with increased disease activity in AS. High MMP levels are also associated with smoking and worse function in AS.


Ankylosing spondylitis (AS) is a chronic inflammatory disease of the spine, three times more common in men than in women. It is characterised by sacroiliitis, accompanied by inflammation of the entheses and the spine [1]. Involvement of peripheral joints occurs in about 30% of patients. Extraskeletal manifestations include acute anterior uveitis, inflammatory bowel symptoms and renal, pulmonary and cardiac involvement. Although genetic factors are important in disease development, there is no clear way of predicting which patients will suffer severe disease [2].

Standard measures of acute phase response (erythrocyte sedimentation rate/C-reactive protein (CRP)) are not generally felt to be reliable indicators of activity of spinal disease [3]: these measures correlate weakly with disease activity in AS, do not fully reflect the disease process, and their value in clinical trials is limited. Neither measure appears to be clearly superior in terms of validity [3].

Although the serum levels of some cytokines and matrix metalloproteinases (MMPs) have been found to be elevated in active disease, there is a general lack of information about biomarker profiles in AS and how these are related to disease activity. Individual cytokines/cytokine receptors that have been associated with disease activity in AS include IL-6, transforming growth factor beta-1, vascular endothelial growth factor (VEGF), macrophage colony-stimulating factor and soluble IL-2 receptor [49]. Reductions in the circulating levels of transforming growth factor beta-1, IL-6, VEGF and macrophage colony-stimulating factor have been seen in patients treated with TNF inhibitors [1014], although not all studies have shown changes after anti-TNF treatment [9, 15].

Recent studies have suggested that serum MMP-3 is a useful marker of disease activity in AS, particularly in patients with peripheral synovitis [9, 16, 17]. A number of studies have shown that MMP-3 levels are reduced in response to TNF antagonists, although its usefulness as a marker of response is debatable [9, 13, 16, 1822]. A recent study has suggested that MMP-3 levels were not useful for monitoring and predicting response to etanercept in terms of disease activity and functional assessments [22]. However, higher levels of MMP-3 have been shown to be predictive of greater radiographic progression in AS, especially in patients with pre-existing damage [23].

Clinical disease activity and particular disease features in AS will possibly be associated with discrete cytokine profiles, and combinations of cytokines and associated markers (for example, MMPs) may be more informative than individual markers. With the availability of technology for measuring many different biomarkers simultaneously from small serum samples, we have carried out an investigation of the relationship between a large panel of cytokines/MMPs/tissue inhibitors of metalloproteinases (TIMPs) and disease activity in a well-characterised population of patients with established AS.

Materials and methods

Patients with AS (n = 180) were recruited from six secondary care rheumatology centres in the UK (Bath, North Bristol, Cannock, Durham, Romford and Stoke-on-Trent). Patients invited to take part in the study were over 21 years of age, had AS according to the modified New York Criteria of 1984 [24], and provided written informed consent according to the declaration of Helsinki. The Trent Research Ethics Committee and the six site-specific National Health Service trusts approved the multicentre study. The majority of patients had been on a stable dose of nonsteroidal anti-inflammatory drugs for at least the previous 3 months. Patients on anti-TNF therapy, systemic steroid use in the previous 3 months or bisphosphonates in the previous 12 months were excluded from the study. Other exclusions included current use of methotrexate, sulfasalazine, raloxifene, calcitonin, phenylbutazone or hormone replacement therapy.

The principal aim of this study was to evaluate the effect of bisphosphonate (alendronate) on global health (Bath ankylosing spondylitis global health (BAS-G)) in AS (Bisphosphonates in Ankylosing Spondylitis trial). Secondary aims were to determine whether there were changes in disease activity, function and bone status, and whether changes in biomarker levels were different between patients treated with or without alendronate. A preliminary report has been published, which indicated that alendronate had no significant effect on changes in CRP, cytokine and MMP levels [25]. A further aim of the study was to investigate the relationship between circulating biomarker levels and measures of disease activity and function in AS patients on standard nonsteroidal anti-inflammatory drug therapy. Data collected from the baseline visit are reported here.

Baseline assessment included questionnaires to assess disease activity (Bath ankylosing spondylitis disease activity index (BASDAI)), function (Bath ankylosing spondylitis functional index (BASFI)) and well-being (BAS-G). The BASDAI is based on six questions related to fatigue, spinal pain, peripheral arthritis, enthesitis and morning stiffness (both severity and duration). The BASFI is a set of 10 questions that consider activities related to functional anatomy and the patients' ability to cope with everyday life. The BAS-G consists of two questions that assess the effect of the disease on the patient's well-being. A 10-cm visual analogue scale is used to answer questions to the BASDAI, BASFI and BAS-G. All measures are scored between 0 and 10, with higher values indicating worse disease activity, function or well-being.

Blood samples were taken for measurement of cytokines, MMPs and TIMPs, as well as a full blood count, urea and electrolyte levels, serum calcium and CRP. Information was also collected on smoking status, smoking duration, average number of cigarettes smoked per day and age of smoking cessation. Pack years were calculated (1 pack year = 20 cigarettes/day for 1 year) to provide a quantitative measure of smoking history. The effects of smoking intensity were assessed by categorising participants according to pack-year history, as in previous studies in AS [26]: category 1, 0 pack years; category 2, 1 to 15 pack years; category 3, 16 to 30 pack years; and category 4, > 30 pack years.

Measurement of cytokines, MMPs and TIMPs

Sera were separated from bloods collected in plain Becton Dickinson Vacutainer® tubes (Becton Dickinson, Oxford, Oxfordshire, UK) at study entry. All sera were stored at -70°C until required. Measurement of the various cytokines, MMPs and TIMPs was performed using multiplex, bead-based (Luminex®) assays on a Bio-Plex™ 200 suspension array system (Bio-Rad Laboratories, Hemel Hempstead, Hertfordshire, UK). The levels of 30 cytokines (Human cytokine 30-plex panel; Life Technologies, Paisley, UK), five MMPs and four TIMPS (Fluorokine multi-analyte MMP and TIMP kits; R&D Systems Europe, Abingdon, UK) were measured in separate multiplex assays according to the manufacturers' instructions. The following biomarkers were measured: cytokines and cytokine receptors - IL-1β, IL-1 receptor antagonist, IL-2, IL-2R, IL-4, IL-5, IL-6, IL-7, IL-10, IL-12p70, IL-13, IL-15, IL-17, IFNα, IFNγ, TNFα, epidermal growth factor, basic fibroblast growth factor, VEGF, hepatocyte growth factor (HGF), granulocyte colony-stimulating factor and granulocyte-macrophage colony-stimulating factor; chemokines - CXCL8, CXCL10, Eotaxin, macrophage inflammatory protein-1α (CCL3), macrophage inflammatory protein-1β (CCL4), monokine induced by gamma interferon (CXCL9), monocyte chemotactic protein-1(CCL2) and regulated upon activation, normal T-cell expressed, and secreted (CCL5); MMPs - MMP-1, MMP-2, MMP-3, MMP-8 and MMP-9; and TIMPs - TIMP-1, TIMP-2, TIMP-3 and TIMP-4. High and low control samples were used in each assay. To test the reproducibility of MMP measurements we also measured the same samples using standard ELISA kits for MMP-3 and MMP-8 (R&D Systems).

Statistical analysis

All data were tested for normality and the appropriate parametric or nonparametric tests were selected. Continuous data were expressed as mean ± standard deviation or median (interquartile range), as appropriate. Univariate correlations between biomarker levels and disease measures were carried out using Spearman's correlation. A multivariate variable selection procedure using the algorithm of McHenry [27] was initially used to select cytokine, MMP or TIMP variables that showed the strongest association with each clinical assessment. Multiple regression analysis or multivariate logistic regression analysis were used to investigate the association between clinical measures and biomarker levels while adjusting for other possible confounders. Where appropriate, data transformation to normality (log or square-root transformation) was carried out before analysis. When cytokine levels were below the level of detection, we carried out imputation of the lowest standard for that particular cytokine [28].

Principal component analysis

Principal component analysis (PCA) is an exploratory technique that reduces the dimensionality of a large number of variables to a smaller set of uncorrelated independent components, and allows the identification of combinations of variables that best explain the differences between observations. PCA allows the identification of patterns within the variables and expresses them in a way that highlights the similarities and differences between them. Principal components (PCs) were extracted using varimax rotation, with the factor selection based on an eigenvalue cutoff of 1.0. The PCs identified were used in multivariate analysis to look for associations with clinical measures.

Hierarchical cluster analysis

Biomarker levels were first converted to log2 and expressed relative to the normalised mean value. These measurements were used to generate heat maps using Genesis software (version 1.7.2; Alexander Sturn, Institute for Genomics and Bioinformatics, Graz University of Technology, Graz, Austria). The Genesis program uses a hierarchical clustering method that enables groups of variables with similar expression levels to be clustered together, as well as grouping together patient samples with similar expression patterns.

Statistical analyses were carried out using Number Cruncher Statistical Software package for Windows (Number Cruncher Statistical System 2000; NCSS Statistical Software, Kaysville, UT, USA). The significance level was set at P = 0.05.


Characteristics of the ankylosing spondylitis patients

Table 1 displays the demographic and clinical characteristics of the 180 AS subjects from the six centres involved in the study. Serum samples were available on 157 patients at baseline. No difference was found in clinical disease measures between patients with and without serum samples available. Information on smoking status (never, past, current) was available for all patients investigated, although the pack-year history was only obtained for 68/96 (70.8%) of patients who had smoked.

Table 1 Demographic and clinical characteristics of the ankylosing spondylitis patients at baseline

Correlation between serum biomarkers and clinical disease measures at baseline

Variables showing significant correlations (Spearman) are shown in Additional File 1. Baseline CRP levels correlated significantly with baseline levels of IL-6 (r = 0.25, P = 0.004), MMP-1 (r = 0.17, P = 0.04), MMP-2 (r = 0.17, P = 0.04), MMP-3 (r = 0.24, P = 0.002), MMP-8 (r = 0.23, P = 0.003) and MMP-9 (r = 0.23, P = 0.003). Baseline CRP was correlated with the BASFI (r = 0.19, P = 0.01) but not the BASDAI or BAS-G, and no correlation was found between IL-6 levels and the BASDAI, BASFI or BAS-G. Significant correlations with all three clinical assessments were found, however, for HGF and MMP-8 levels (r ≥ 0.20, P ≥ 0.02). For the BASDAI, significant correlations were also found with CXCL10 (r = 0.17, P = 0.04) and MMP-9 levels (r = 0.24, P = 0.002). MMP-9 levels were also correlated with the BASFI (r = 0.17, P = 0.03), while MMP-1 levels were correlated with both the BASFI and BAS-G (r = 0.20, P = 0.01 and r = 0.19, P = 0.01, respectively).

The reproducibility of MMP-3 and MMP-8 correlations with clinical measures was tested by also measuring these MMPs using standard ELISA kits (R&D Systems). The results were similar to those achieved using the Luminex® system. MMP-3 levels were correlated with CRP levels but not with the BASDAI, BASFI or BAS-G. MMP-8 levels were significantly correlated with CRP and all three clinical assessments (Additional file 2).

Multivariate analyses

A multivariate variable selection procedure was initially used to select cytokine, MMP or TIMP variables that showed the strongest association with each clinical assessment. To further determine which variables were independently associated with clinical measures we then carried out multiple regression analyses for each clinical assessment (dependent variable) in models that were adjusted for age, sex and disease duration, and included all significant variables found in the preliminary analyses.

Associations with C-reactive protein

MMP-2, MMP-3 and MMP-8 levels were found to be independently associated with the CRP level (log transformed) at baseline (Table 2). MMP-3 and MMP-8 were positively associated while MMP-2 showed a negative association. No association of IL-6 or any other cytokine, MMP or TIMP was found in this model. Adjustment for age, sex and disease duration made little or no difference to the association of CRP levels with MMP levels. A similar model of association was found when MMP-3 and MMP-8 levels measured by ELISA replaced those obtained by Luminex® analysis (data not shown).

Table 2 Multiple regression analysis showing variables associated with CRP levels in ankylosing spondylitis patients at baseline

Associations with the BASDAI

Female sex and MMP-8 levels provided the best model for association with the BASDAI (square root transformed) at baseline (Table 3, Model 1), although a very similar model was found with female sex and MMP-9 levels as independent variables (Additional file 3). MMP-8 and MMP-9 levels were not independently associated since both lost significance in models containing both MMPs as independent variables. CRP was independently associated in a model adjusted for female sex only, but lost significance in models containing MMP-8 or MMP-9 levels. No other MMP, TIMP or cytokine were associated in these models. MMP-8 levels measured by ELISA were also significantly associated with the BASDAI in a model that also contained female sex (data not shown).

Table 3 Multiple regression models showing variables associated with the BASDAI, BASFI and BAS-G at baseline

Associations with the BASFI

The strongest association with the BASFI (square root transformed) was found in a model containing age and CRP levels as independent variables (Table 3, Model 2). There were no independent associations of any cytokine, MMP or TIMP with the BASFI in models that contained CRP.

Associations with the BAS-G

As for the BASDAI, female sex and MMP-8 levels provided the strongest independent association with the BAS-G (Table 3, Model 3). In this case, MMP-8 levels remained independently associated in models that also contained MMP-9 as a variable (Additional file 3). A similar model was obtained when MMP-8 levels measured by ELISA replaced those obtained by Luminex® analysis (data not shown). No association was found with CRP levels, although this approached significance (P = 0.07) in models that did not include MMP-8.

Principal component analysis

We carried out PCA to investigate patterns of cytokines and MMPs in the data. This allowed the identification of factors containing particular biomarker profiles, which were examined for their association with clinical assessments in AS. Eight PCs were identified, explaining 78.3% of the total variance. These PCs are shown in Table 4 along with their respective factor loadings after varimax rotation.

Table 4 Principal component analysis of serum biomarkers in ankylosing spondylitis patients at baseline

The only PC showing an association with any clinical assessments was PC3, consisting of MMP-8, MMP-9, CXCL8 and HGF (Table 4). This component was associated with both CRP (log transformed) and the BASDAI (square root transformed). The latter association was independent of CRP, and remained significant after adjustment for age, sex and disease duration (Table 5). Interestingly, PC3 was also significantly higher in current smokers than in past smokers (0.455 vs. 0.089, P = 0.004) or in patients who had never smoked (0.455 vs. -0.350, P = 0.01). There was also a significant trend (P = 0.0005) of increasing PC3 level with increasing pack-year category (Table 6). Analysis of the individual markers of the PC3 profile demonstrated that MMP-8 and MMP-9 were both higher in current smokers than in past or nonsmokers (Table 7). MMP-8 levels also showed an increasing trend (P < 0.0001) with increasing pack-year category, although no trend was seen with MMP-9 levels (Additional file 4).

Table 5 Multiple regression analysis showing baseline association of PC3 (MMP-8, MMP-9, HGF, CXCL8) with the BASDAI
Table 6 Association between principal component 3 and pack-year category
Table 7 Levels of individual biomarkers comprising principal component 3 stratified by smoking status

Hierarchical clustering analysis

Hierarchical cluster analysis was carried out on cytokines, MMPs and TIMPS separately and in combination. Heat maps generated from these analyses were used to identify patient groups that clustered together according to the similarity or dissimilarity of their profiles.

Discrimination between patient groups was best achieved by hierarchical cluster analysis of the five MMPs alone. This resulted in clustering of two major groups of patients based on relatively low or high MMP levels (Figure 1). The high MMP group was characterised by significantly higher levels of MMP-1 (20% higher), MMP-3 (43% higher), MMP-8 (222% higher) and MMP-9 (128% higher) (Additional file 5). There was no difference in MMP-2 levels between the two groups.

Figure 1
figure 1

Hierarchical cluster analysis of serum matrix metalloproteinase levels in patients with ankylosing spondylitis. Results are displayed as a heat map and dendrogram in which the relative levels of matrix metalloproteinases (MMPs) are represented by shades of yellow/red (high) and blue (low). Each row represents the MMP profile for an individual patient, each of which is represented by a number on the vertical axis. Each column represents a different MMP (MMP-1, MMP-2, MMP-3, MMP-8 and MMP-9, left to right). C1, high MMP cluster; C2, low MMP cluster.

The high MMP group of patients had significantly higher CRP levels (median 12.0 vs. 6.0, P = 0.001), higher BASDAI scores (median, 4.67 vs. 3.27, P = 0.006) and higher BASFI scores (median, 3.83 vs. 2.94, P = 0.03), although the association with the BAS-G (median 5.0 vs. 3.6, P = 0.1) was not significant. The level of PC3 was also significantly higher in the high MMP group (1.125 vs. -0.511, P < 0.0001). In addition, patients in the high MMP group were more likely to be current smokers than those in the low MMP group (36.5% vs. 17.9%, odds ratio 2.60, 95% CI 1.26 to 5.38, P = 0.01). An increasing frequency of the high MMP group was also found with increasing pack-year category (P trend = 0.001) (Table 8).

Table 8 Relationship between pack-year category and frequency of AS patients in low and high MMP clusters


In this study we analysed the serum levels of a large panel of cytokines, chemokines, MMPs and TIMPs for their association with measures of disease activity and function in AS, using univariate and multivariate methods. We believe this is the largest study to date to examine such a large panel of biomarkers in AS. Our data indicate that serum levels of MMPs in particular show the strongest association with disease activity, as measured by the BASDAI. In contrast to some previous studies, we found that the strongest association was with MMP-8 (or MMP-9) levels, rather than MMP-3 levels. Importantly, we obtained similar findings using the Luminex® system and ELISA methods. It is worth noting that previously reported relationships between MMP-3 levels and the BASDAI have often been based on univariate correlations rather than multivariate regression analysis, and correlations have not been consistently found in all studies [9, 14, 16, 17, 23]. This may reflect differences in the patient groups since it has also been suggested that MMP-3 levels are higher in patients with peripheral arthritis than those with axial disease only [14, 16, 17], although this has not always been observed [23]. The same may apply to IL-6, which has also been shown to be higher in patients with peripheral synovitis [5] but is not associated with the BASDAI in all studies [14]. One of the limitations of the present study was the absence of information on peripheral joint disease in the patients studied.

The correlations of MMP-3 and IL-6 with disease activity in some studies may be explained in part by the association of these molecules with CRP levels. Like previous studies, we found correlations of MMP-3 and IL-6 with CRP levels - although we also found that MMP-2 and MMP-8 levels were independently associated in a multivariate model containing MMP-3. IL-6 was not associated in this model. Interestingly the association with MMP-2 was negative, suggesting an anti-inflammatory role for this MMP, which is consistent with other studies in inflammatory arthritis [29, 30].

The associations of MMP-8 and MMP-9 levels with the BASDAI were not independent of each other, and probably reflect involvement of the same pathway in the release of these particular MMPs. MMP-8 (neutrophil collagenase) is primarily produced by activated neutrophils, while MMP-9 is released mainly by neutrophils and macrophages.

In the PCA, MMP-8 and MMP-9 are also associated with CXCL8 and HGF (in the PC3 profile), both of which have also been associated with neutrophil activation [31, 32]. The association of the BASDAI with this component may thus reflect a possible relationship between disease activity and neutrophil activation in AS. However, the biomarkers making up the PC3 profile are also associated with aspects of angiogenesis [3335], although this may be linked also to an association with neutrophil activation. For example, neutrophil-derived MMPs can promote neoangiogenesis through release and activation of angiogenesis promoters (for example, VEGF-A) from the extracellular matrix [33]. Other studies in AS have shown a correlation between VEGF levels and disease activity [8, 13, 14], although we failed to find an association in this study. Further work is needed to determine with which particular aspects of the disease process the markers associated with PC3 are most closely associated. The association of PC3 with the BASDAI independent of CRP suggests that it is not merely a surrogate for systemic inflammation.

It is noteworthy that PC3 was also associated with current smoking, and that a significant trend was seen in relation to pack-year history. This appears to be largely due to associations with MMP-8 and MMP-9 levels, which have been shown in other studies to be increased in smokers [36, 37]. Hierarchical clustering analysis also revealed an association between patients with high MMP levels and current smoking, as well as a quantitative relationship with pack-year history. All of these data add weight to the idea that the PC3 profile may reflect neutrophil activation, which is characteristically found in smokers [38, 39]. A number of studies have indicated that smoking is associated with increased disease activity, worse functional outcome and poorer quality of life in patients with AS [26, 4043]. Smoking has also been associated with radiographic spinal progression in early axial spondyloarthritis independently of baseline radiographic damage and elevated acute phase reactants [44]. These associations may be explained in part by the increased levels of MMPs (particularly MMP-8 and MMP-9) in smokers. This association is suggested by the current study since the association of the BASDAI and BASFI with the pack-year category disappears when MMP-8 or MMP-9 levels are included in multivariate regression models along with pack-year history (data not shown). Smoking may thus exacerbate the production of neutrophil activation markers, but there also appears to be an association of these markers with disease activity independent of smoking.

The association of disease activity with markers of neutrophil and/or macrophage activation is interesting in light of studies suggesting that the innate immune pathway might be more important in axial spondyloarthritis than the adaptive immune response [4548]. Infiltration of the synovium in axial spondyloarthritis by macrophage subsets and neutrophils has been shown to reflect global disease activity [47], while more recently the role of IL-17 in axial spondyloarthritis has been linked to increased numbers of IL-17+ neutrophils and macrophage subsets in the subchondral bone marrow of affected facet joints [48]. Other cell types such as CD3+ T cells, mast cells, B cells and natural killer cells provided little or no source of IL-17. The hypothesis that the innate immune pathway, mostly mediated through neutrophils, might be of greater relevance in AS inflammation than the T-helper type 17 cell-mediated adaptive immune response was therefore suggested. Such a suggestion is in line with earlier studies that indicated significant alterations in neutrophil function in patients with AS [4954]. These include alterations in chemotaxis, phagocytosis and superoxide radical anion generation. Priming of neutrophils has been suggested to be a probable causative factor in the onset of AS, and increased production of reactive oxygen species and MMPs from these cells may lead to tissue damage in AS [53]. Increased levels of advanced oxidation protein products, a novel oxidative stress marker of protein, have been demonstrated in AS patients, and provide evidence of oxidative stress mediated by neutrophil myeloperoxidase-hypochlorous acid in these patients [54]. The cause of increased neutrophil activation in AS is unclear - although smoking is likely to have an exacerbating effect, and cessation of smoking may therefore provide one area for limited control of disease.


We have shown, using several methods of analysis, that higher disease activity and worse function in patients with established AS are associated with increased serum levels of MMPs (particularly MMP-8 and MMP-9) and cytokines/chemokines (HGF, CXCL8) associated with neutrophil activation and/or angiogenesis. Serum MMP-8 and/or MMP-9 levels are more strongly associated with disease activity than MMP-3, and biomarker profiles containing high MMP-8/MMP-9 levels are associated with smoking in AS.

Further work is needed to investigate the role of these particular molecules in determining radiographic outcome in AS, and their potential use as markers of response to therapy.



ankylosing spondylitis


Bath ankylosing spondylitis disease activity index


Bath ankylosing spondylitis functional index


Bath ankylosing spondylitis global health


C-reactive protein


enzyme-linked immunosorbent assay


hepatocyte growth factor






matrix metalloproteinase


principal component


principal component analysis


tissue inhibitor of metalloproteinase


tumour necrosis factor


vascular endothelial growth factor.


  1. Gran JT, Husby G: Ankylosing spondylitis: prevalence and demography. Rheumatology. Edited by: Klippel JH, Dieppe PA. 1998, London: Mosby, 15.1-15.6. 2

    Google Scholar 

  2. Hamersma J, Cardon LR, Bradbury L, Brophy S, van-der Horst-Bruinsma I, Calin A, Brown MA: Is disease severity in ankylosing spondylitis genetically determined?. Arthritis Rheum. 2001, 44: 1396-1400. 10.1002/1529-0131(200106)44:6<1396::AID-ART233>3.0.CO;2-A.

    Article  CAS  PubMed  Google Scholar 

  3. Ruof J, Stucki G: Validity aspects of erythrocyte sedimentation rate and C-reactive protein in ankylosing spondylitis: a literature review. J Rheumatol. 1999, 26: 966-970.

    CAS  PubMed  Google Scholar 

  4. Gratacós J, Collado A, Fillela X, Sanmarti R, Canete J, Llena J, Molina R, Ballesta A, Muñoz-Gómez J: Serum cytokines (IL-6, TNF-alpha, IFN-gamma) in ankylosing spondylitis: a close correlation between serum IL-6 and disease activity and severity. Br J Rheumatol. 1994, 33: 927-931. 10.1093/rheumatology/33.10.927.

    Article  PubMed  Google Scholar 

  5. Claudepierre P, Rymer J-C, Authier F-J, Allanore Y, Larget-Piet B, Gherardi R, Chevalier X: A relationship between TGF-β1 or IL-6 plasma levels and clinical features of spondyloarthropathies. Br J Rheumatol. 1997, 36: 400-401. 10.1093/rheumatology/36.3.400.

    Article  CAS  PubMed  Google Scholar 

  6. Bal A, Unlu E, Bahar G, Aydog E, Eksioglu E, Yorgancioglu R: Comparison of serum IL-1 beta, IL-2R, IL-6, and TNF-alpha levels with disease activity parameters in ankylosing spondylitis. Clin Rheumatol. 2007, 26: 211-215.

    Article  CAS  PubMed  Google Scholar 

  7. Howe HS, Cheung PL, Kong KO, Badsha H, Thong BY, Leong KP, Koh ET, Lian TY, Cheng YK, Lam S, Teo D, Lau TC, Leung BP: Transforming growth factor beta-1 and gene polymorphisms in oriental ankylosing spondylitis. Rheumatology. 2005, 44: 51-54. 10.1093/rheumatology/keh426.

    Article  CAS  PubMed  Google Scholar 

  8. Drouart M, Saas P, Billot M, Cedoz JP, Tiberghien P, Wendling D, Toussirot E: High serum vascular endothelial growth factor correlates with disease activity of spondyloarthropathies. Clin Exp Immunol. 2003, 132: 158-162. 10.1046/j.1365-2249.2003.02101.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Yang C, Gu J, Rihl M, Baeten D, Huang F, Zhao M, Zhang H, Maksymowych WP, De Keyser F, Veys EM, Yu DTY: Serum levels of matrix metalloproteinase 3 and macrophage colony stimulating factor 1 correlate with disease activity in ankylosing spondylitis. Arthritis Rheum (Arthritis Care Res). 2004, 51: 691-699. 10.1002/art.20696.

    Article  CAS  Google Scholar 

  10. Woo JH, Lee HJ, Sung IH, Kim TH: Changes of clinical response and bone biochemical markers in patients with ankylosing spondylitis taking etanercept. J Rheumatol. 2007, 34: 1753-1759.

    CAS  PubMed  Google Scholar 

  11. Appel H, Janssen L, Listing J, Heydrich R, Rudwaleit M, Sieper J: Serum levels of biomarkers of bone and cartilage destruction and new bone formation in different cohorts of patients with axial spondyloarthritis with and without tumor necrosis factor-alpha blocker treatment. Arthritis Res Ther. 2008, 10: R125-10.1186/ar2537.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Visvanathan S, Wagner C, Marini JC, Baker D, Gathany T, Han J, van der Heijde D, Braun J: Inflammatory biomarkers, disease activity and spinal measures in patients with ankylosing spondylitis after treatment with infliximab. Ann Rheum Dis. 2008, 67: 511-517.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Pedersen SJ, Hetland ML, Sorensen IJ, Ostergaard M, Nielsen H, Johansen JS: Circulating levels of interleukin-6, vascular endothelial growth factor, YKL-40, matrix metalloproteinase-3, and total aggrecan in spondyloarthritis patients during 3 years of treatment with TNFα inhibitors. Clin Rheumatol. 2010, 29: 1301-1309. 10.1007/s10067-010-1528-x.

    Article  PubMed  Google Scholar 

  14. Pedersen SJ, Sørensen IJ, Garnero P, Johansen JS, Madsen OR, Tvede N, Hansen MS, Thamsborg G, Smedegaard L, Majgaard O, Loft AG, Erlendsson J, Asmussen K, Jurik AG, Møller J, Hasselquist M, Mikkelsen D, Skjødt T, Lambert R, Hansen A, Østergaard M: ASDAS, BASDAI and different treatment responses and their relation to biomarkers of inflammation, cartilage and bone turnover in patients with axial spondyloarthritis treated with with TNFα inhibitors. Ann Rheum Dis. 2011, 70: 1375-1381. 10.1136/ard.2010.138883.

    Article  CAS  PubMed  Google Scholar 

  15. Keller C, Webb A, Davis J: Cytokines in the seronegative spondyloarthropathies and their modification by TNF blockade: a brief report and literature review. Ann Rheum Dis. 2003, 62: 128-132.

    Article  Google Scholar 

  16. Vandooren B, Kruithof E, Yu DTY, Rihl M, Gu J, De Rycke L, Van den Bosch F, Veys EM, De Keyser F, Baeten D: Involvement of matrix metalloproteinases and their inhibitors in peripheral synovitis and down-regualtion by tumor necrosis factor α blockade in spondyloarthropathy. Arthritis Rheum. 2004, 50: 2942-2953. 10.1002/art.20477.

    Article  CAS  PubMed  Google Scholar 

  17. Chen C-H, Lin K-C, Yu DTY, Yang C, Huang F, Chen H-A, Liang T-H, Liao H-T, Tsai C-Y, Wei JCC, Chou C-T: Serum matrix metalloproteinases and tissue inhibitors of metalloproteinases in ankylosing spondylitis: MMP-3 is a reproducibly sensitive and specific biomarker of disease activity. Rheumatology. 2006, 45: 414-420. 10.1093/rheumatology/kei208.

    Article  CAS  PubMed  Google Scholar 

  18. Maksymowych WP, Poole RA, Hiebert L, Webb A, Ionescu M, Lobanok T, King L, Davis JC: Etanercept exerts beneficial effects on articular cartilage biomarkers of degradation and turnover in patients with ankylosing spondylitis. J Rheumatol. 2005, 32: 1911-1917.

    CAS  PubMed  Google Scholar 

  19. Maksymowych WP, Rahman P, Shojania K, Olsynski WP, Thomson GT, Ballal S, Wong RL, Inman RD, M03-606 Study Group, Bell MJ, Keystone E, Martin WO, Beaulieu A, Choquette D, Chardonnet L, He Y, Nissen M, Lush L: Beneficial effects of adalimumab on biomarkers reflecting structural damage in patients with ankylosing spondylitis. J Rheumatol. 2008, 35: 2030-2037.

    CAS  PubMed  Google Scholar 

  20. Wendling D, Cedoz J-P, Racador E: Serum levels of MMP-3 and cathepsin K in patients with anklosing spondyltis: effect of TNFα antagonist therapy. Joint Bone Spine. 2008, 75: 559-562. 10.1016/j.jbspin.2008.01.026.

    Article  CAS  PubMed  Google Scholar 

  21. Romero-Sanchez C, Robinson WH, Tomooka BH, Londono J, Valle-Onate R, Huang F, Deng X, Zhang L, Yang C, Yu DT: Identification of acute phase reactants and cytokines useful for monitoring infliximab therapy in ankylosing spondylitis. Clin Rheumatol. 2008, 27: 1429-1435. 10.1007/s10067-008-0941-x.

    Article  PubMed  Google Scholar 

  22. Arends S, van der Veer E, Groen H, Houtman PM, Jansen TL, Leijsma MK, Bijzet J, Limburg PC, Kallenberg CG, Spoorenberg A, Brouwer E: Serum MMP-3 level as a biomarker for monitoring and predicting response to etanercept treatment in ankylosing spondylitis. J Rheumatol. 2011, 38: 1644-1650. 10.3899/jrheum.101128.

    Article  CAS  PubMed  Google Scholar 

  23. Maksymowych WP, Landewe R, Conner-Spady B, Dougados M, Mielants H, van der Tempel H, Poole RA, Wang N, van der Heijde D: Serum matrix metalloproteinase 3 is an independent predictor of structural damage progression in patients with ankylosing spondylitis. Arthritis Rheum. 2007, 56: 1846-1853. 10.1002/art.22589.

    Article  CAS  PubMed  Google Scholar 

  24. Van der Linden S, Valkenburg HA, Cats A: Evaluation of diagnostic criteria for ankylosing spondylitis. Arthritis Rheum. 1984, 27: 361-368. 10.1002/art.1780270401.

    Article  CAS  PubMed  Google Scholar 

  25. Coates L, Bhalla A, Creamer P, Hailwood S, Mattey D, Packham J, Taylor G, BIAS Study Group: The effect of alendronate on ankylosing spondylitis; the results of the bisphosphonates in ankylosing spondylitis trial (BIAS). Ann Rheum Dis. 2011, 70 (Suppl 3): 128.Abstract.-

    Google Scholar 

  26. Mattey DL, Dawson SR, Healey EL, Packham JC: Relationship between smoking and patient reported measures of disease outcome in ankylosing spondylitis. J Rheumatol. 2011, 38: 2608-2615. 10.3899/jrheum.110641.

    Article  PubMed  Google Scholar 

  27. McHenry CE: Computation of a best subset in multivariate analysis. Appl Stat. 1978, 27: 291-296. 10.2307/2347164.

    Article  Google Scholar 

  28. Little RJA, Rubin DB: Statistical Analysis with Missing Data. 2002, Wiley-Interscience; New York, 2

    Google Scholar 

  29. Itoh T, Matsuda H, Tanioka M, Kuwabara K, Itohara S, Suzuki R: The role of matrix metalloproteinase-2 and matrix metalloproteinase-9 in antibody induced arthritis. J Immunol. 2002, 169: 2643-2647.

    Article  CAS  PubMed  Google Scholar 

  30. Rodriguez-Lopez J, Perez-Pampin E, Gomez-Reino JJ, Gonzalez A: Regulatory polymorphisms in extracellular matrix protease genes and susceptibility to rheumatoid arthritis: a case-control study. Arthritis Res Ther. 2006, 8: R1-10.1186/ar1849.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Scapini P, Lapinet-Vera JA, Gasperini S, Calzetti F, Bazzoni F, Cassatella MA: The neutrophil as a cellular source of chemokines. Immunol Rev. 2000, 177: 195-203. 10.1034/j.1600-065X.2000.17706.x.

    Article  CAS  PubMed  Google Scholar 

  32. Matsushima A, Ogura H, Koh T, Fujita K, Yoshiya K, Sumi Y, Hosotsubo H, Kuwagata Y, Tanaka H, Shimazu T, Sugimoto H: Hepatocyte growth factor in polymorphonuclear leukocytes is increased in patients with systemic inflammatory response syndrome. J Trauma. 2004, 56: 259-264. 10.1097/01.TA.0000111752.60500.DA.

    Article  CAS  PubMed  Google Scholar 

  33. van Hinsbergh VW, Engelse MA, Quax PH: Pericellular proteases in angiogenesis and vasculogenesis. Arterioscler Thromb Vasc Biol. 2006, 26: 716-728. 10.1161/01.ATV.0000209518.58252.17.

    Article  CAS  PubMed  Google Scholar 

  34. Rosenkilde MM, Schwartz TW: The chemokine system - a major regulator of angiogenesis in health and disease. APMIS. 2004, 112: 481-495. 10.1111/j.1600-0463.2004.apm11207-0808.x.

    Article  CAS  PubMed  Google Scholar 

  35. Sengupta S, Gherardi E, Sellers LA, Wood JM, Sasisekharan R, Fan TP: Hepatocyte growth factor/scatter factor can induce angiogenesis independently of vascular endothelial growth factor. Arterioscler Thromb Vasc Biol. 2003, 23: 69-75. 10.1161/01.ATV.0000048701.86621.D0.

    Article  CAS  PubMed  Google Scholar 

  36. Aquilante CL, Beitelshees AL, Zineh I: Correlates of serum matrix metalloproteinase-8 (MMP-8) concentrations in nondiabetic subjects without cardiovascular disease. Clin Chim Acta. 2007, 379: 48-52. 10.1016/j.cca.2006.12.006.

    Article  CAS  PubMed  Google Scholar 

  37. Garvin P, Nilsson L, Carstensen J, Jonasson L, Kristenson M: Circulating matrix metalloproteinase-9 is associated with cardiovascular risk factors in a middle aged normal population. PloS One. 2008, 3: e1774-10.1371/journal.pone.0001774. doi: 10.1371/journal.pone.0001774

    Article  PubMed Central  PubMed  Google Scholar 

  38. Blann AD, Kirkpatrick U, Devine C, Naser S, McCollum CN: The influence of acute smoking on leucocytes, platelets and the endothelium. Atherosclerosis. 1998, 141: 133-139. 10.1016/S0021-9150(98)00163-4.

    Article  CAS  PubMed  Google Scholar 

  39. Gustaffson A, Asman B, Bergstrom K: Cigarette smoking as an aggravating factor in inflammatory tissue-destructive diseases. Increase in tumor-necrosis factor-alpha priming of peripheral neutrophils measured as generation of oxygen radicals. Int J Clin Lab Res. 2000, 30: 187-190. 10.1007/s005990070005.

    Article  Google Scholar 

  40. Averns HL, Oxtoby J, Taylor HG, Jones PW, Dziedzic K, Dawes PT: Smoking and outcome in ankylosing spondylitis. Scand J Rheumatol. 1996, 25: 138-142. 10.3109/03009749609080003.

    Article  CAS  PubMed  Google Scholar 

  41. Ward MM: Predictors of the progression of functional disability in patients with ankylosing spondylitis. J Rheumatol. 2002, 29: 1420-1425.

    PubMed  Google Scholar 

  42. Reed MD, Dharmage S, Boers A, Martin BJ, Buchanan RR, Schachna L: Ankylosing spondylitis: an Australian experience. Intern Med J. 2008, 38: 321-327. 10.1111/j.1445-5994.2007.01471.x.

    Article  CAS  PubMed  Google Scholar 

  43. Chung HY, Machado P, van der Heijde D, D'Agostino MA, Dougados M: Smokers in early axial spondyloarthritis have earlier disease onset, more disease activity, inflammation and damage, and poorer function and health-related quality of life: results from the DESIR cohort. Ann Rheum Dis. 2012, 71: 809-816. 10.1136/annrheumdis-2011-200180.

    Article  PubMed  Google Scholar 

  44. Poddubnyy D, Haibel H, Listing J, Marker-Hermann E, Zeidler H, Braun J, Sieper J, Rudwaleit M: Baseline radiographic damage, elevated acute phase reactants and cigarette smoking status predict radiographic progression in the spine in early axial spondyloarthritis. Arthritis Rheum. 2012, 64: 1388-1398. 10.1002/art.33465.

    Article  PubMed  Google Scholar 

  45. Baeten D, Van Damme N, Van den Bosch F, Kruithof E, De Vos M, Mielants H, Veys EM, De Keyser F: Impaired Th1 cytokine production in spondyloarthropathy is restored by anti-TNFα. Ann Rheum Dis. 2001, 60: 750-755. 10.1136/ard.60.8.750.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Baeten D, Kruithof E, De Rycke L, Boots AM, Mielants H, Veys EM, de Keyser F: Infiltration of the synovial membrane with macrophage subsets and polymorphonuclear cells reflects global disease activity in spondyloarthropathy. Arthritis Res Ther. 2005, 7: R359-R369. 10.1186/ar1501.

    Article  PubMed Central  PubMed  Google Scholar 

  47. McConagle D, McDermott MF: A proposed classification of the immunological diseases. PloS Med. 2006, 3: e297-10.1371/journal.pmed.0030297.

    Article  Google Scholar 

  48. Appel H, Maier R, Wu P, Scheer R, Hempfing A, Kayser R, Thiel A, Radbruch A, Loddenkemper C, Sieper J: Analysis of IL-17+ cells in facet joints of patients with spondyloarthritis suggests that the innate immune pathway might be of greater relevance than the Th17-mediated adaptive immune response. Arthritis Res Ther. 2011, 13: R95-10.1186/ar3370.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Repo H, Leirisalo-Repo M, Koivuranta-Vaara P: Exaggerated inflammatory responsiveness plays a part in the pathogenesis of HLA-B27 linked disease - hypothesis. Ann Clin Res. 1984, 16: 47-50.

    CAS  PubMed  Google Scholar 

  50. Pease CT, Fennell M, Brewerton DA: Polymorphonuclear leucocyte motility in men with ankylosing spondylitis. Ann Rheum Dis. 1989, 48: 35-41. 10.1136/ard.48.1.35.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Wendling D, Didier JM, Vuitton DA: The phagocyte oxidative metabolism function in ankylosing spondylitis. Rheumatol Int. 1991, 11: 187-189. 10.1007/BF00332560.

    Article  CAS  PubMed  Google Scholar 

  52. Biasi D, Carletto A, Caramaschi P, Bellavite P, Andrioli G, Caraffi M, Pacor ML, Bambara LM: Neutrophil functions, spondyloarthropathies and HLA-B27: a study of 43 patients. Clin Exp Rheumatol. 1995, 13: 623-627.

    CAS  PubMed  Google Scholar 

  53. Ho K-J, Chen P-Q, Chang C-Y, Lu F-J: The oxidative metabolism of circulating phagocytes in ankylosing spondyltis: determination by whole blood chemiluminescence. Ann Rheum Dis. 2000, 59: 338-341. 10.1136/ard.59.5.338.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Yazici C, Kose K, Calis M, Kuzuguden S, Kirnap M: Protein oxidation status in patients with ankylosing spondylitis. Rheumatology. 2004, 43: 1235-1239. 10.1093/rheumatology/keh317.

    Article  CAS  PubMed  Google Scholar 

Download references


Bisphosphonates in Ankylosing Spondylitis study group: Dr D Mulherin, Prof. K Chakravarty, Dr P Lewis, and Ms M Griffin; British Society for Rheumatology; and Arthritis Research UK. The authors are grateful for help from Miss Ashley Bhalla in the data collection and entry. The work was supported by the British Society for Rheumatology and the Arthritis Research Campaign (now Arthritis Research UK).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Derek L Mattey.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DLM and NBN carried out the biomarker measurements. DLM and GJT carried out the statistical analysis. JCP, LC, PC, SH and AKB participated in the design of the study, and recruitment of patients. DLM and JCP conceived the study, participated in its design and coordination, carried out analysis and interpretation of data, and drafted the final manuscript. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1 presenting correlations between clinical measures and biomarkers in ankylosing spondylitis patients at baseline. (PDF 18 KB)


Additional file 2: Table S2 presenting correlations between clinical measures and MMP-3 and MMP-8 levels measured by ELISA in ankylosing spondylitis patients at baseline. (PDF 8 KB)


Additional file 3: Table S3 presenting alternative multiple regression models showing variables associated with the BASDAI and BAS-G at baseline. (PDF 8 KB)


Additional file 4: Table S4 presenting levels of individual biomarkers comprising principal component 3 stratified by pack-year category. (PDF 47 KB)


Additional file 5: Table S5 presenting a comparison of MMP levels in patient clusters with low or high MMP levels as selected by hierarchical cluster analysis. (PDF 45 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

Reprints and permissions

About this article

Cite this article

Mattey, D.L., Packham, J.C., Nixon, N.B. et al. Association of cytokine and matrix metalloproteinase profiles with disease activity and function in ankylosing spondylitis . Arthritis Res Ther 14, R127 (2012).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: