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

Gut dysbiosis associated with worse disease activity and physical function in axial spondyloarthritis

Arthritis Research & Therapy volume 24, Article number: 42 (2022) Cite this article

Abstract

Background

Based on clinical and genetic associations, axial spondyloarthritis (axSpA) and inflammatory bowel disease (IBD) are suspected to have a linked pathogenesis. Gut dysbiosis, intrinsic to IBD, has also been observed in axSpA. It is, however, not established to what degree gut dysbiosis is associated with axSpA disease severity.

The objective of this study was to compare gut dysbiosis frequency between controls, non-radiographic axial spondyloarthritis (nr-axSpA), and ankylosing spondylitis (AS) patients and investigate whether gut dysbiosis is cross-sectionally associated with axSpA disease activity, physical function, mobility, or pain.

Methods

Gut dysbiosis was assessed by 16SrRNA analysis of feces from 44/88 nr-axSpA/AS patients (ASAS/mNY criteria) without inflammatory bowel disease (IBD) and 46 controls without IBD or rheumatic disease. The GA-map™ Dysbiosis Test was used, grading gut microbiota aberrations on a 1-5 scale, where ≥3 denotes dysbiosis. Proportions with dysbiosis were compared between the groups. Furthermore, standard axSpA measures of disease activity, function, mobility, and pain were compared between patients (nr-axSpA and AS combined) with and without dysbiosis, univariately, and adjusted for relevant confounders (ANCOVA).

Results

Gut dysbiosis was more frequent in AS than controls (36% versus 17%, p=0.023), while nr-axSpA (25% dysbiosis) did not differ significantly from either AS or controls. Univariately, most axSpA measures were significantly worse in patients with dysbiosis versus those without: ASDAS-CRP between-group difference 0.6 (95% CI 0.2–0.9); BASDAI 1.6 (0.8–2.4); evaluator’s global disease activity assessment (Likert scale 0–4) 0.3 (0.1–0.5), BASFI 1.5 (0.6–2.4), and VAS pain (cm) 1.3 (0.4–2.2). Differences remained significant after adjustment for demographics, lifestyle factors, treatments, gut inflammation (fecal calprotectin ≥50 mg/kg), and gut symptoms, except for VAS pain. BASMI and CRP were not associated with dysbiosis.

Conclusion

Gut dysbiosis, more frequent in AS patients than controls, is associated with worse axSpA disease activity and physical function, seemingly irrespective of both gut inflammation and treatments. This provides further evidence for an important link between disturbances in gastrointestinal homeostasis and axSpA.

Background

The link between spondyloarthritis (SpA) and inflammatory bowel disease (IBD) is well recognized. IBD patients are at increased risk of developing SpA, and the reverse is likewise true, with prevalence estimates of diagnosed IBD in axial SpA (axSpA) cohorts between 5 and 10% [1,2,3]. Several shared genetic risk loci are known [4, 5], and a common pathobiological pathway suggested [6, 7]. Additionally, approximately 50–60% of SpA patients display microscopic intestinal inflammation in biopsies of the ileum or colon, often reminiscent of Crohn’s disease [8,9,10].

IBD is associated with an altered intestinal bacteria composition [11], a condition known as intestinal or gut dysbiosis [12]. The intestinal inflammation characterizing IBD is thought to arise through a complex interplay between genetic, environmental, immunological, and gut microbial factors [13]. A variation of microbial aberrations has been reported in IBD, but a recurrent finding is a reduced species richness, particularly within the Firmicutes phylum [14,15,16]. Moreover, several studies have shown a decreased presence of Fecalibacterium prausnitzii, a species with protective, immune regulatory features, and increased abundance of Escherichia coli, possessing pathogenic properties [11, 17].

Recently, several studies have reported gut dysbiosis also in ankylosing spondylitis (AS) [18,19,20,21,22,23], whereas little is known about gut microbiota in non-radiographic axSpA (nr-axSpA) patients [24]. Furthermore, previous studies have linked intestinal inflammation, assessed either histologically [10, 25, 26], or estimated via elevated fecal calprotectin (F-calprotectin) levels [27, 28], to more active axSpA disease and worse prognosis. Yet, it remains largely unknown whether gut microbiota aberrations alone, independent of intestinal inflammation, is associated with more severe axSpA.

Methods

The objectives of this study were to compare the frequency and degree of gut dysbiosis between well-characterized nr-axSpA and AS patients without known IBD and controls and to investigate whether gut dysbiosis is cross-sectionally associated with worse axSpA disease activity, physical function, mobility, or pain.

Study population

Nr-axSpA and AS patients for the present study were enrolled from the population-based SPARTAKUS cohort of validated axSpA cases from southern Sweden. Comprehensive information about the SPARTAKUS cohort and its inclusion and diagnosis validation process has been published earlier [27]. Briefly, all patients from a defined geographical area of Skåne county, Sweden, with at least one outpatient visit to the Department of Rheumatology, Skåne University Hospital, with an ICD-10 diagnosis consistent with axial SpA (M45.9/M46.0/M46.1/M46.8/M46.9) during 2011–2014, were invited to enroll in the cohort. Among patients with a clinical diagnosis of undifferentiated SpA (M46.8/M46.9), only those reporting back pain ≥3 months with an onset before age 45 were eligible. Enrolled patients were extensively characterized through questionnaires, clinical examinations, blood/feces/urine sampling and imaging (X-ray and, if needed for the diagnosis validation/classification, also magnetic resonance imaging [MRI] of the sacroiliac [SI] joints), enabling classification into nr-axSpA or AS according to the Assessment of SpondyloArthritis international Society (ASAS) axSpA and modified New York criteria, respectively [29, 30]. For the present gut dysbiosis study, patients were consecutively enrolled from the SPARTAKUS cohort, whereby the first 150 subjects classified as nr-axSpA (n=50) and AS (n=100) without comorbid IBD, and who had provided a fecal sample (available for 88% of SPARTAKUS patients), were included. The study visits for these patients occurred 2015–2018.

A control group of 50 persons without rheumatic disease or IBD, and frequency-matched to the overall SPARTAKUS cohort by age and sex, was gathered among friends/colleagues/relatives of the authors. The controls followed a shortened protocol, including questionnaires and provision of blood and fecal samples (study visits during 2018).

Patients (n=18) and controls (n=4) reporting treatment with antibiotics within the last 3 months before they provided the fecal sample were excluded from the present study, with the final study population thus comprising 44 patients with nr-axSpA, 88 with AS and 46 controls.

Gut microbiota analysis

Gut microbiota was analyzed using the GA-map™ Dysbiosis Test (Genetic Analysis, Oslo, Norway), a validated method designed to identify and grade gut dysbiosis in fecal samples by means of a pre-determined target approach, assessing the presence and abundance of a selected panel of 48 bacteria markers at different taxonomic levels and translating this into a dysbiosis index [31]. In brief, the test employs DNA probes, based on the 16SrRNA sequence in seven variable regions (V3–V9), to measure the abundance of bacteria according to the strength of fluorescent signal detection (probe signal intensity). The DNA probes included in the GA-map™ Dysbiosis Test were selected based on published observations regarding gut microbiota aberrations in IBD and IBS (irritable bowel syndrome), and in all target ≥300 bacteria in the six phyla Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, Tenericutes, and Verrucomicrobia.

Based on bacterial composition and abundance by this sampling technique and applying a healthy Scandinavian control group as reference of normobiosis (normal microbiota composition), an algorithm was then developed to grade microbial aberrations from normobiosis, with a resulting Dysbiosis Index (DI) score of 1–5 [31]. A state of gut dysbiosis is considered present at DI ≥3 and severe dysbiosis at DI 5. A validation study of the GA-map™ Dysbiosis Test found dysbiosis (DI ≥3) in 16% of healthy volunteers and 74% of IBD patients [31]. Regarding the determination of gut dysbiosis, good agreement has also been demonstrated between the GA-map™ Dysbiosis Test and deep-sequencing by MiSeq Illumina sequencing technology [31].

The participants of the present study were instructed to collect their fecal sample as close in time to the study visit as possible and to store it in a −18°C freezer before delivery. After the visit, samples were stored frozen at −80°C until analysis (all samples were analyzed at the same time).

Outcomes

Mean DI score 1–5, as well as the frequency of dysbiosis (DI ≥3), were compared between axSpA patients (AS and nr-axSpA combined) and controls. Similar comparisons were also performed separately for nr-axSpA versus controls, AS versus controls, and for nr-axSpA versus AS.

Distributions of the following standard axSpA measures were assessed in relation to DI score 1–5 and presence of gut dysbiosis (DI ≥3), respectively: measures of disease activity (Ankylosing Spondylitis Disease Activity Score using CRP [ASDAS-CRP]; Bath Ankylosing Spondylitis Disease Activity Index [BASDAI]; Evaluator’s Global disease activity assessment [EvalGlobal, Likert scale 0–4, signifying remission/low/medium/high/maximal]; CRP), physical function (Bath Ankylosing Spondylitis Functional Index [BASFI]), mobility (Bath Ankylosing Spondylitis Metrology Index [BASMI]), and pain (patient’s visual analog scale for pain [VAS pain]).

While, as outlined above, the focus of the present study was on the DI as a global measure of gut microbiota aberrations, in order to also allow for some comparisons with prior findings, exploratory analyses comparing the probe signal intensities for the 48 bacterial markers included in the GA-map™ Dysbiosis Test between the axSpA patients (nr-axSpA and AS combined) and controls are presented in Additional file 1.

Statistics

Demographic, disease, and treatment characteristics were compared between groups by chi2-test, Fisher’s exact test or Mann-Whitney U test, as appropriate. Between-group comparisons of DI score 1–5 were performed univariately, as well as adjusted for age, sex, body mass index (BMI), and smoking (categorized as never smoker/quit ≥6 months ago/intermittent smoker or quit <6 months ago/everyday smoker) by analysis of covariance (ANCOVA). Due to the ordinal nature of the DI score 1-5, non-parametric bootstrapping with 10,000 iterations was used to calculate 95% confidence intervals (CI). BMI and smoking were adjusted for based on evidence that both are associated with intestinal microbiota alterations [32, 33]. The frequency of dysbiosis (DI ≥3) was compared between groups using chi2 analysis.

Comparisons of overall between-group differences in the various axSpA measures (regarding disease activity, function, mobility, and pain) between patients with different levels of DI score 1–5 were performed using one-way analysis of variance (ANOVA) or Kruskal-Wallis test, as appropriate. DI levels 4 and 5 were combined into one group, as there were too few cases with each score separately (n=3/n=5). Bootstrapped 95% CI were used for non-normally distributed measures (with normality assessed by Shapiro-Wilk test).

Comparisons of the axSpA measures (regarding disease activity, function, mobility, and pain) between patients with and without dysbiosis (DI ≥3) were performed univariately, as well as multivariately by ANCOVA, again applying bootstrapped 95% CI for non-normally distributed measures. The multivariate analyses were adjusted for age, sex, BMI, smoking (categorized as above), axSpA subtype (nr-axSpA/AS), HLA-B27 status, ongoing anti-TNF (tumor necrosis factor) therapy (yes/no), ASAS NSAID score for the preceding 3 months [34], gut inflammation as measured by F-calprotectin ≥50 mg/kg (yes/no; Calpro AS, Lysaker, Norway), and symptoms meeting the ROME III criteria for IBS (yes/no) [35]. Adjustment for axSpA subtype was considered relevant based on the numerically higher frequency of gut dysbiosis in AS, as compared to nr-axSpA, observed in the present study (see the “Results” section below) and the well-known differences in measures such as CRP and BASMI between the groups. Regarding treatments, anti-TNF agents are used to treat IBD and have been shown to affect gut microbiota in axSpA [20, 36], while NSAIDs can contribute to intestinal inflammation and also have potential to alter the gut microbiome [37]. Gut inflammation, as measured by F-calprotectin elevation, has been previously associated with both worse dysbiosis [19] and axSpA disease activity and function [27, 28] and was thus adjusted for to try to separate the effect of dysbiosis from that of inflammation. Finally, since a previously published study on the SPARTAKUS cohort found a high prevalence (30%) of IBS symptoms and that such symptoms were linked to comorbid fibromyalgia and worse patient-reported axSpA measures of disease activity, function, and pain [38], it was also deemed relevant to adjust for this, as gut dysbiosis is central to IBS [39].

To further ascertain that any observed associations between gut dysbiosis (DI ≥3) and axSpA measures (regarding disease activity, function, mobility, and pain) were linked to dysbiosis, and not to simultaneous intestinal inflammation or IBS, sensitivity analyses were performed excluding patients with elevated F-calprotectin (≥50 mg/kg) and IBS symptoms, respectively.

Welch’s t test with bootstrapped 95%CI (due to non-normal distributions; 10,000 iterations) was used for the exploratory comparisons of probe signal intensities for the 48 bacteria markers of the GA-map™ Dysbiosis Test between axSpA patients and controls.

Logarithmic values of ASAS NSAID score were used due to its skewed distribution. No imputations of missing data were performed. p<0.05 was considered statistically significant. SPSS, version 27 (IBM Corporation, NY, USA) was used for statistical analyses.

Results

Characteristics of the study population

Characteristics of the axSpA patients and controls are displayed in Table 1. Controls were well matched to the total patient group for age and gender. Compared to the AS group, nr-axSpA patients were on average younger and had shorter symptom duration, and a larger proportion were females. Ongoing or prior smoking was more common among AS patients than in the nr-axSpA or control groups. A numerically larger proportion of AS than nr-axSpA patients displayed gut inflammation (F-calprotectin ≥50 mg/kg), while presence of IBS symptoms was evenly distributed between the two axSpA groups. Similar proportions of nr-axSpA and AS patients received conventional synthetic and/or biologic DMARDs, with 41% anti-TNF use in both groups.

Table 1 Characteristics of the study population
Full size table

Regarding axSpA measures, BASMI values were significantly higher and CRP values numerically higher in the AS group, while other measures of disease activity, function, and pain were evenly distributed between the axSpA subtypes.

Gut dysbiosis in axSpA and controls

DI score 1–5 was significantly higher among axSpA patients than controls (unadjusted between-group difference 0.45 [95% CI 0.16 to 0.74], p=0.003; adjusted between-group difference 0.34 [0.04 to 0.64], p=0.028); Fig. 1A). Yet, the frequency of gut dysbiosis (DI ≥3) was only numerically more common in axSpA versus controls (33% versus 17%; p=0.050). Of the 8 controls with dysbiosis, none displayed more pronounced aberrations from normobiosis, as measured by a DI of 4 or 5, whereas in the axSpA group this was observed in 19% of the dysbiotic patients (n=8 out of 43 patients with dysbiosis).

Fig. 1
figure 1

Dysbiosis index distributions among axSpA patients and controls. A Distributions of dysbiosis index (DI) score 1–5 among axSpA patients (nr-axSpA and AS combined) and controls. Unadjusted between-group difference in DI score 1–5: 0.45 (95% CI 0.16 to 0.74), p=0.003. Adjusted (for age, sex, BMI, and smoking) between-group difference in DI score 1–5: 0.34 (0.04 to 0.64), p=0.028. B Distributions of DI score 1–5 among nr-axSpA patients, AS patients, and controls. Nr-axSpA versus controls: unadjusted between-group difference in DI score 1–5: 0.25 (−0.09 to 0.60), p=0.154; adjusted (as above): 0.17 (−0.19 to 0.54), p=0.409. AS versus controls: unadjusted: 0.55 (0.22 to 0.88), p<0.001; adjusted: 0.43 (0.07 to 0.78), p=0.019. AS versus nr-axSpA: unadjusted: 0.30 (−0.06 to 0.64), p=0.098; adjusted: 0.20 (−0.18 to 0.56), p=0.286. Number of subjects per group presented at the bottom of (or above) the bars. AS ankylosing spondylitis, axSpA axial spondyloarthritis, BMI body mass index, CI confidence interval, DI dysbiosis index, nr-axSpA non-radiographic axial spondyloarthritis

Full size image

In our predefined subgroup analysis, both DI score 1–5 (results in Fig. 1B) and presence of gut dysbiosis were significantly higher in the AS group compared to controls (DI ≥3: 36% versus 17%; p=0.023). Conversely, for nr-axSpA, neither DI score 1–5 (results in Fig. 1B), nor the frequency of gut dysbiosis differed significantly from controls (DI ≥3: 25% versus 17%; p=0.377), although in both cases the nr-axSpA point-estimates fell between those of the control and AS groups. Analogously, no significant between-group differences in DI score 1–5 (results in Fig. 1B) or presence of dysbiosis (p=0.189) were detected between the nr-axSpA and AS groups.

Results of the exploratory analyses comparing probe signal intensity for the 48 bacterial markers of the GA-map™ Dysbiosis Test between axSpA patients and controls are presented in Additional file 1.

AxSpA measures in relation to dysbiosis index levels

Among the axSpA patients (nr-axSpA and AS combined), higher DI score 1–5 (combining DI 4 and 5 into one group, due to too few cases) was associated with worse scores of most assessed axSpA measures (Fig. 2). Significant overall between-group differences across the DI score 1–5 categories were found for ASDAS-CRP, BASDAI, BASFI, and VAS pain, but not for EvalGlobal, CRP, or BASMI (Figs. 2 and 3; p=0.055 for EvalGlobal).

Fig. 2
figure 2

AxSpA measures in relation to dysbiosis index levels. Distributions of measures of disease activity, physical function, mobility, and pain among axSpA patients (nr-axSpA and AS combined) with different levels of dysbiosis index (DI) score 1–5. Number of analyzed subjects per group presented above the x-axes (DI 4 and 5 combined into one group due to too few cases). P values for overall between-group differences by one-way ANOVA (or Kruskal-Wallis test †) displayed in the graphs. *Bootstrapped 95% CI. ANOVA analysis of variance, AS ankylosing spondylitis, axSpA axial spondyloarthritis, ASDAS-CRP ankylosing spondylitis disease-activity score using C-reactive protein, BASDAI Bath ankylosing spondylitis disease activity index, BASFI Bath ankylosing spondylitis functional index, BASMI Bath ankylosing spondylitis metrology index, CI confidence interval, CRP C-reactive protein, DI dysbiosis index, nr-axSpA non-radiographic axial spondyloarthritis, VAS visual analog scale

Full size image
Fig. 3
figure 3

Evaluator’s global assessment of disease activity in relation to dysbiosis index levels. Distribution of Evaluator’s global assessment of disease activity (EvalGlobal; Likert scale 0–4, corresponding to Remission/Low/Medium/High/Maximal) among axSpA patients (nr-axSpA and AS combined) with different levels of dysbiosis index (DI) score 1–5. Number of patients in each group: DI 1 n=51, DI 2 n=33, DI 3 n=35, DI 4, or 5 n=8 (DI 4 and 5 combined into one group due to too few cases). No patients in any DI group were assessed as having high or maximal disease activity. P=0.055 for overall between-group difference by Kruskal-Wallis test. AS ankylosing spondylitis, axSpA axial spondyloarthritis, DI dysbiosis index, nr-axSpA non-radiographic axial spondyloarthritis

Full size image

Differences in axSpA measures in patients with versus without gut dysbiosis

When instead splitting the axSpA group (nr-axSpA and AS combined) according to presence of gut dysbiosis (DI <3 [n=89] versus ≥3 [n=43]), patients displaying dysbiosis had higher BMI and were more often ongoing or prior smokers (Table 2). HLA-B27 positivity was less frequent in the dysbiotic group, while gut inflammation (F-calprotectin ≥50 mg/kg) was numerically and IBS symptoms significantly more common in subjects with dysbiosis (Table 2).

Table 2 Characteristics of the axSpA patients (nr-axSpA and AS combined), stratified by the presence of gut dysbiosis
Full size table

Comparing the axSpA measures between patients with versus without gut dysbiosis, univariate analyses showed significantly higher ASDAS-CRP, BASDAI, BASFI, VAS pain, and EvalGlobal scores for patients with versus without dysbiosis (Fig. 4; Table 2). These between-group differences remained significant also after adjustment, except for VAS pain (p=0.066) (Fig. 4). Regarding CRP and BASMI, no differences were detected between patients with versus without dysbiosis, neither univariately nor in the adjusted analyses (Fig. 4).

Fig. 4
figure 4

Differences in axSpA measures in patients with versus without gut dysbiosis. Differences in measures of disease activity, physical function, mobility, and pain between axSpA patients (nr-axSpA and AS combined) with gut dysbiosis (DI ≥3) versus those without dysbiosis (DI <3). The data shown represent point-estimate differences (dots) with 95% CI (whiskers) and corresponding p values from unadjusted analyses, and after adjustment (*) for age, sex, BMI, smoking, axSpA subtype (nr-axSpA/AS), HLA-B27 status (positive/negative), ongoing anti-TNF therapy (yes/no), ASAS 3-month NSAID score, gut inflammation (F-calprotectin ≥50 mg/kg, yes/no), and IBS-symptoms (yes/no) (ANCOVA). † Bootstrapped 95% CI. ANCOVA analysis of covariance, AS ankylosing spondylitis, ASAS Assessment of SpondyloArthritis international Society, axSpA axial spondyloarthritis, ASDAS-CRP ankylosing spondylitis disease-activity score using C-reactive protein, BASDAI Bath ankylosing spondylitis disease activity index, BASFI Bath ankylosing spondylitis functional index, BASMI Bath ankylosing spondylitis metrology index, BMI body mass index, CI confidence interval, CRP C-reactive protein, DI dysbiosis index, EvalGlobal Evaluator’s global assessment of disease activity, F fecal, HLA human leukocyte antigen, IBS irritable bowel syndrome, nr-axSpA non-radiographic axial spondyloarthritis, NSAID non-steroidal anti-inflammatory drug, TNF tumor necrosis factor, VAS visual analog scale

Full size image

In the sensitivity analyses, the comparisons displayed in Fig. 4 were repeated after having respectively limited the study population to patients without gut inflammation (F-calprotectin <50 mg/kg; n=84) or to those not reporting IBS symptoms (n=88). In both cases, between-group differences for patients with versus without dysbiosis seen in the main analysis were confirmed or found even more distinct regarding BASDAI and VAS pain, and in relation to gut inflammation also for ASDAS-CRP and BASFI (see Additional file 1, Figures S2 and S3). However, regarding EvalGlobal, no significant difference remained between the groups after adjustment in the patient sample without gut inflammation (p=0.061; see Additional file 1, Figure S2). Also in the sample without IBS symptoms, no between-group differences were observed for EvalGlobal in either the univariate or adjusted analysis, nor for ASDAS-CRP or BASFI in the adjusted analysis (see Additional file 1, Figure S3).

Discussion

Main findings

In this population-based, cross-sectional study of well-characterized axSpA patients without IBD, we found gut dysbiosis to be significantly more frequent in AS (36%) than controls (17%). Likewise, the degree of gut microbiota aberration (as measured by DI score 1–5) was significantly higher in AS than controls. Conversely, neither dysbiosis frequency nor DI score 1–5 were significantly increased in nr-axSpA compared to controls, although numerically the nr-axSpA point-estimates fell between those of the AS patients and controls.

More importantly (and to the best of our knowledge not this clearly demonstrated before), in the overall axSpA group (nr-axSpA and AS combined), a state of gut dysbiosis was associated with significantly worse disease activity, physical function, and pain as measured by ASDAS-CRP, BASDAI, EvalGlobal, BASFI, and VAS pain. The associations remained also after adjustment for relevant confounders, including gut inflammation, anti-TNF, and NSAID treatment, except regarding VAS pain (where statistical significance was just barely lost). Restricting the analyses to include only patients without gut inflammation (as measured by elevated F-calprotectin ≥50 mg/kg) further demonstrated the associations to be present regardless of such inflammation. Thus, similar to prior observations regarding gut inflammation in axSpA [10, 25,26,27,28], but seemingly independent of such inflammation, gut dysbiosis may be a marker of more severe axSpA.

Previous research

Gastrointestinal pathology, including dysbiosis, and the pathogenesis of SpA

Gut dysbiosis is believed to play an important role in IBD pathogenesis through complex and incompletely understood interactions with genetic and environmental risk factors, gut barrier defects, and abnormal immune responses [13]. SpA is closely associated with IBD—in particular Crohn’s disease—both genetically and in terms of clinical comorbidity, and many of their shared genetic risk loci govern mucosal immune defenses [7]. Furthermore, several pathological, gastrointestinal processes and aberrations reminiscent of those in IBD have been observed also in SpA, including gut dysbiosis [18,19,20,21,22,23,24, 40], gut inflammation (present in 50–60% of SpA patients, although subclinical in most cases) [8,9,10], and signs of gut barrier defects [41]. Based on this, one central hypothesis for SpA pathogenesis implicates such IBD-like, gastrointestinal processes, including dysbiosis, as an upstream event, subsequently inducing an IL-23/IL-17-pathway driven inflammation of the spine and joints [6, 7, 42, 43]. The present results, linking gut dysbiosis to worse axSpA phenotype and more active disease appear to fit well with such a theory. This hypothesis remains unproven, however, and the close links between IBD/IBD-like pathology and SpA may also be due to their shared genetic predisposition, without a causal relationship [44]. Regardless of which, ongoing gut inflammation has been repeatedly associated with more severe SpA [10, 25,26,27,28], demonstrating an important interaction between the gut and the musculoskeletal apparatus in this disease.

Gut dysbiosis in axSpA

The current results are congruent with previous findings which have repeatedly associated AS with altered gut microbiota [18,19,20,21,22,23], whereas little prior information is available regarding nr-axSpA [24]. While compositional analyses of the gut microbiome have consistently been able to distinguish AS from healthy controls [18,19,20,21, 23], the specific microbiota alterations observed have been heterogeneous, and no specific changes at any level of the bacterial taxonomic hierarchy have so far been identified as typical for axSpA. Furthermore, several features reminiscent of the gut dysbiosis of IBD have been reported in AS [18,19,20,21, 40]. Compositional analyses, however, clearly distinguish the gut microbiome of AS from that of IBD [19, 23], showing that despite certain similarities, the type of dysbiosis in the two diseases is not the same.

As outlined above, 50–60% of SpA patients display microscopic gut inflammation [8,9,10], but whether there is a causal relationship between this and gut dysbiosis—and if so, in which direction—remains incompletely understood. The two processes are, however, intertwined, since the presence of inflammation has been associated with a significantly different gut microbiota composition compared to SpA patients without histological gut inflammation [24]. Furthermore, another study showed the gut microbiota composition to differ between AS patients with elevated versus normal F-calprotectin [19].

In the present study, gut microbiota aberrations were assessed by the validated, semi-quantitative DI of the GA-map™ Dysbiosis Test. Klingberg et al. previously applied the same method in a cohort of 150 AS patients, finding a substantially higher gut dysbiosis (DI ≥3) prevalence of 87% [19], as compared to the 36% in our AS sample. Of interest in light of our results, however, their cohort also displayed higher disease activity and worse physical function, with median BASDAI/BASFI scores of 3.2/2.3 versus 2.2/1.2 in our AS group, and thus represented a patient group with worse disease.

AxSpA disease severity in relation to gut dysbiosis

To our knowledge, only a few studies have previously investigated gut microbiota aberrations in SpA in relation to disease severity, most seeming to indicate an association with disease activity in line with our results. Tito et al. studied the intestinal microbiota in biopsies from 27 axSpA patients (AS and nr-axSpA) and found a significant positive correlation between the abundance of the genus Dialister and disease activity as measured by ASDAS (rs=0.62), and a similar trend for BASDAI (rs=0.53) [24]. In another study, Breban et al. observed an enrichment in Ruminococcus gnavus in SpA, which was significant in relation to controls only for cases with active disease (BASDAI ≥3) [40]. This enrichment in association with higher SpA disease activity was observed regardless of concomitant IBD, but the elevation was more pronounced in those with comorbid IBD, for whom there was also a significant positive correlation between Ruminucoccus gnavus abundance and BASDAI (rs=0.77). Finally, a recent abstract from a multinational study on gut microbiota in AS showed that the bacterial composition differed significantly between patients with different BASDAI levels (categorized as <2.5/2.5-5/5-7.5/>7.5), congruent with our findings [23].

Conversely, Klingberg et al. reported no associations between gut microbiota composition, assessed by the GA-map™ Dysbiosis Test, and ASDAS-CRP, BASDAI, BASFI, BASMI, CRP, or ESR (erythrocyte sedimentation rate) [19]. However, in contrast to our analyses, their study did not relate these measures to the DI, but performed compositional analysis, comparing the gut microbiota composition between patient groups with better or worse (dichotomized) levels of the different indices. Notably, as stated above, Klingberg et al. found dysbiosis (DI ≥3) in as many as 87% of their AS patients (n=130 versus n=20 without dysbiosis), a distribution which may have limited their ability to detect differences in disease measures based on microbiota composition.

Apart from excluding [19, 24] or separately analyzing [40] patients with comorbid IBD, none of the above-cited studies [19, 23, 24, 40] adjusted their analyses regarding the link between gut microbiota aberrations and disease severity for intestinal inflammation. Thus, some reported observations may have been driven by the known link between gut inflammation and more severe disease [10, 25,26,27,28], rather than by the dysbiosis per se. Here, we both adjusted for F-calprotectin and performed a sensitivity analysis excluding all patients with elevated F-calprotectin (≥50 mg/kg), aiming to more specifically assess a potential association with gut dysbiosis. In this context, it should be acknowledged that the ability of a F-calprotectin <50 mg/kg to rule out asymptomatic, microscopic (histological) gut inflammation, which is much more common in SpA than overt IBD [8, 27], has not been widely studied, and F-calprotectin may also be less sensitive for small-bowel than colonic inflammation [45]. Cypers et al., however, did report significantly elevated F-calprotectin in SpA patients displaying microscopic gut inflammation in the terminal ileum and/or colon, with an optimal F-calprotectin cut-off of 85 mg/kg for the detection of such cases [46]. Furthermore, the presence of macroscopic inflammation at capsule endoscopy and/or ileocolonoscopy in SpA patients with previously undiagnosed IBD has also been shown to be significantly associated with elevated F-calprotectin levels (>100 mg/kg) [47, 48]. Taken together, the current results seem to indicate that gut dysbiosis may be associated with worse axSpA disease activity and physical function, independently of intestinal inflammation.

Similarly, since IBS symptoms have previously been associated with worse patient-reported axSpA measures in the SPARTAKUS cohort [38] and that gut dysbiosis is intrinsic to IBS [39], the same adjustment and sensitivity analysis approach was also applied for IBS symptoms, again without relevant changes to our results. Anti-TNF therapy, used by 41% of our patients, has previously shown a potential to ameliorate gut dysbiosis in axSpA [20, 36] and was therefore adjusted for, without any considerable impact on the findings. Moreover, proton-pump inhibitor (PPI) use was considered as an additional adjustment factor due to its described association with gut dysbiosis [49], but was omitted since PPI-use during the last 3 months did not differ significantly between patients with versus without dysbiosis (DI ≥3) in our material (p=0.300 by chi2-test) and since no causal relationship between PPI use and the axSpA measures was suspected.

Strengths and limitations

This study provides new knowledge regarding how gut dysbiosis relates to axSpA phenotype and disease measures (regarding disease activity, function, mobility, and pain)—associations previously investigated by only a limited number of studies. Patients were consecutively enrolled from a population-based cohort study and were well characterized and classified according to the ASAS axSpA and modified New York criteria. The detailed protocol allowed for extensive adjustment for possible confounders, including BMI, smoking, axSpA subtype, HLA-B27 status, anti-TNF and NSAID therapy, intestinal inflammation, and IBS symptoms. By including F-calprotectin, we could—to a certain degree—disentangle a link between gut dysbiosis and worse axSpA disease activity and function from the impact of simultaneous intestinal inflammation.

Another strength is that, in addition to AS, our study also encompasses nr-axSpA, a disease subtype for which data regarding gut microbiota are scant. The relatively low number of nr-axSpA patients (n=44) means that we cannot rule out type 2 error from contributing to the lack of significant differences between this group and both AS and controls. Yet, the pattern by which the nr-axSpA point-estimates fall in between those of AS and controls is similar to prior observations for F-calprotectin [27]. The relatively long mean symptom duration of our nr-axSpA group should also be remembered when assessing the results, since it likely entails a lower average risk of future progression to AS than in a cohort of newly diagnosed nr-axSpA.

Regarding generalizability, only comorbid IBD and recent antibiotic use were exclusion criteria. Without restrictions regarding disease subtype (nr-axSpA or AS) or severity, treatment, or other comorbidities, we believe our results to be fairly generalizable to the wider axSpA population. However, in the absence of any concurrent ICD-10 codes for SpA (usually added in case of axial involvement), patients clinically diagnosed with psoriatic, IBD-associated, or reactive arthritis (by ICD-10 codes) were not included, although these conditions may sometimes have axial disease.

Regarding limitations, endoscopic examinations and biopsies were not part of the SPARTAKUS protocol, but could have provided valuable information regarding histological inflammation and mucosa-associated microbiota (fecal and mucosa-associated microbiota have been shown to differ [50, 51]). The use of deep-sequencing or shotgun metagenomics sequencing methods, rather than the pre-determined target approach of the GA-map™ Dysbiosis Test, would have provided a much more detailed characterization of the gut microbiota. Assessment of associations between specific bacteria and axSpA disease measures were, however, beyond the scope of the present study. For the determination of gut dysbiosis, the validated GA-map™ Dysbiosis Test has shown good agreement with deep-sequencing methodology [31], and by the DI, it provides a readily comprehensible, semiquantitative score of global dysbiosis, calibrated to healthy individuals as reference of normobiosis. Finally, information regarding probiotics use is missing, and the cross-sectional design did not allow for longitudinal evaluation of microbiota shifts or disease development.

Conclusions

The presence of gut dysbiosis is linked to worse axSpA phenotype and more active disease. While dysbiosis was more prevalent in AS than controls, no significant difference was observed between nr-axSpA and controls. In the overall axSpA group (nr-axSpA and AS combined), the presence of gut dysbiosis was independently associated with worse disease activity and physical function, seemingly irrespective of both gut inflammation and ongoing treatments. Further research is needed to assess the validity of these findings also in other axSpA cohorts, as well as to gain a deeper understanding of the connection between gut microbiota and axSpA, and examine whether a causal relationship exists. This would offer guidance as to whether targeting the intestine therapeutically, with regard to the luminal microbiota composition and/or mucosal inflammation, could be a relevant approach to try to interfere with axSpA development or mitigate the disease course. Nonetheless, the current results add to the steadily growing body of evidence for an important, potentially pathobiological, link between disturbances in gastrointestinal homeostasis and axSpA.

Availability of data and materials

The datasets used and analyzed for the current study are available from the corresponding author on reasonable request.

Abbreviations

ANCOVA:

Analysis of covariance

AS:

Ankylosing spondylitis

ASAS:

Assessment of SpondyloArthritis international Society

ASDAS-CRP:

Ankylosing spondylitis disease activity score using CRP

axSpA:

Axial spondyloarthritis

BASDAI:

Bath ankylosing spondylitis disease activity index

BASFI:

Bath ankylosing spondylitis functional index

BASMI:

Bath ankylosing spondylitis metrology index

bDMARD:

Biologic disease-modifying anti-rheumatic drug

BMI:

Body mass index

CI:

Confidence interval

CRP:

C-reactive protein

csDMARD:

Conventional synthetic disease-modifying anti-rheumatic drug

DI:

Dysbiosis index

DMARD:

Disease-modifying anti-rheumatic drug

ESR:

Erythrocyte sedimentation rate

EvalGlobal:

Evaluator’s global assessment of disease activity

HLA-B27:

Human leukocyte antigen B27

IBD:

Inflammatory bowel disease

IBS:

Irritable bowel syndrome

IL:

Interleukin

IQR:

Interquartile range

mNY:

Modified New York criteria

Nr-axSpA:

Non-radiographic axial spondyloarthritis

NSAID:

Non-steroidal anti-inflammatory drug

PPI:

Proton-pump inhibitor

SD:

Standard deviation

TNF:

Tumor necrosis factor

VAS pain:

Visual analog scale for pain

References

  1. Stolwijk C, van Tubergen A, Castillo-Ortiz JD, Boonen A. Prevalence of extra-articular manifestations in patients with ankylosing spondylitis: a systematic review and meta-analysis. Ann Rheum Dis. 2015;74:65–73.

    Article  PubMed  Google Scholar 

  2. van den Berg R, de Hooge M, van Gaalen F, Reijnierse M, Huizinga T, van der Heijde D. Percentage of patients with spondyloarthritis in patients referred because of chronic back pain and performance of classification criteria: experience from the Spondyloarthritis Caught Early (SPACE) cohort. Rheumatology (Oxford). 2013;52:1492–9.

    Article  Google Scholar 

  3. Dougados M, Etcheto A, Molto A, Alonso S, Bouvet S, Daures JP, et al. Clinical presentation of patients suffering from recent onset chronic inflammatory back pain suggestive of spondyloarthritis: The DESIR cohort. Joint Bone Spine. 2015;82:345–51.

    Article  PubMed  Google Scholar 

  4. Ellinghaus D, Jostins L, Spain SL, Cortes A, Bethune J, Han B, et al. Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci. Nat Genet. 2016;48:510–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Reveille JD. Genetics of spondyloarthritis--beyond the MHC. Nat Rev Rheumatol. 2012;8:296–304.

    Article  CAS  PubMed  Google Scholar 

  6. Asquith M, Elewaut D, Lin P, Rosenbaum JT. The role of the gut and microbes in the pathogenesis of spondyloarthritis. Best Pract Res Clin Rheumatol. 2014;28:687–702.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Choi RY, Asquith M, Rosenbaum JT. Fecal transplants in spondyloarthritis and uveitis: ready for a clinical trial? Curr Opin Rheumatol. 2018;30:303–9.

    Article  PubMed  Google Scholar 

  8. De Vos M, Cuvelier C, Mielants H, Veys E, Barbier F, Elewaut A. Ileocolonoscopy in seronegative spondylarthropathy. Gastroenterology. 1989;96:339–44.

    Article  PubMed  Google Scholar 

  9. Simenon G, Van Gossum A, Adler M, Rickaert F, Appelboom T. Macroscopic and microscopic gut lesions in seronegative spondyloarthropathies. J Rheumatol. 1990;17:1491–4.

    CAS  PubMed  Google Scholar 

  10. Van Praet L, Van den Bosch FE, Jacques P, Carron P, Jans L, Colman R, et al. Microscopic gut inflammation in axial spondyloarthritis: a multiparametric predictive model. Ann Rheum Dis. 2013;72:414–7.

    Article  PubMed  Google Scholar 

  11. Pittayanon R, Lau JT, Leontiadis GI, Tse F, Yuan Y, Surette M, et al. Differences in gut microbiota in patients with vs without inflammatory bowel diseases: a systematic review. Gastroenterology. 2020;158:930–46 e1.

    Article  PubMed  Google Scholar 

  12. Manasson J, Blank RB, Scher JU. The microbiome in rheumatology: where are we and where should we go? Ann Rheum Dis. 2020;79:727–33.

    Article  PubMed  Google Scholar 

  13. de Souza HS, Fiocchi C. Immunopathogenesis of IBD: current state of the art. Nat Rev Gastroenterol Hepatol. 2016;13:13–27.

    Article  PubMed  Google Scholar 

  14. Manichanh C, Rigottier-Gois L, Bonnaud E, Gloux K, Pelletier E, Frangeul L, et al. Reduced diversity of faecal microbiota in Crohn’s disease revealed by a metagenomic approach. Gut. 2006;55:205–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Dicksved J, Halfvarson J, Rosenquist M, Jarnerot G, Tysk C, Apajalahti J, et al. Molecular analysis of the gut microbiota of identical twins with Crohn's disease. ISME J. 2008;2:716–27.

    Article  CAS  PubMed  Google Scholar 

  16. Rajca S, Grondin V, Louis E, Vernier-Massouille G, Grimaud JC, Bouhnik Y, et al. Alterations in the intestinal microbiome (dysbiosis) as a predictor of relapse after infliximab withdrawal in Crohn’s disease. Inflamm Bowel Dis. 2014;20:978–86.

    PubMed  Google Scholar 

  17. Gevers D, Kugathasan S, Denson LA, Vazquez-Baeza Y, Van Treuren W, Ren B, et al. The treatment-naive microbiome in new-onset Crohn’s disease. Cell Host Microbe. 2014;15:382–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Costello ME, Ciccia F, Willner D, Warrington N, Robinson PC, Gardiner B, et al. Brief report: intestinal dysbiosis in ankylosing spondylitis. Arthritis Rheum. 2015;67:686–91.

    Article  Google Scholar 

  19. Klingberg E, Magnusson MK, Strid H, Deminger A, Stahl A, Sundin J, et al. A distinct gut microbiota composition in patients with ankylosing spondylitis is associated with increased levels of fecal calprotectin. Arthritis Res Ther. 2019;21:248.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Yin J, Sternes PR, Wang M, Song J, Morrison M, Li T, et al. Shotgun metagenomics reveals an enrichment of potentially cross-reactive bacterial epitopes in ankylosing spondylitis patients, as well as the effects of TNFi therapy upon microbiome composition. Ann Rheum Dis. 2020;79:132–40.

    Article  CAS  PubMed  Google Scholar 

  21. Zhang L, Han R, Zhang X, Fang G, Chen J, Li J, et al. Fecal microbiota in patients with ankylosing spondylitis: correlation with dietary factors and disease activity. Clin Chim Acta. 2019;497:189–96.

    Article  CAS  PubMed  Google Scholar 

  22. Wen C, Zheng Z, Shao T, Liu L, Xie Z, Le Chatelier E, et al. Quantitative metagenomics reveals unique gut microbiome biomarkers in ankylosing spondylitis. Genome Biol. 2017;18:142.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Sternes P, Brett L, Phipps J, Ciccia F, de Guzman E, Morrison M, et al. Gut microbiome changes are different between ankylosing spondylitis and inflammatory bowel disease, and correlate with disease activity in both diseases. Arthritis. Rheumatol. 2020;72(suppl 10). https://acrabstracts.org/abstract/gut-microbiome-changes-are-different-between-ankylosing-spondylitis-andinflammatory-bowel-disease-and-correlate-with-disease-activity-in-both-diseases/. Accessed 3 Feb 2022.

  24. Tito RY, Cypers H, Joossens M, Varkas G, Van Praet L, Glorieus E, et al. Brief report: dialister as a microbial marker of disease activity in spondyloarthritis. Arthritis Rheum. 2017;69:114–21.

    Article  CAS  Google Scholar 

  25. Van Praet L, Jans L, Carron P, Jacques P, Glorieus E, Colman R, et al. Degree of bone marrow oedema in sacroiliac joints of patients with axial spondyloarthritis is linked to gut inflammation and male sex: results from the GIANT cohort. Ann Rheum Dis. 2014;73:1186–9.

    Article  PubMed  Google Scholar 

  26. Mielants H, Veys EM, Cuvelier C, De Vos M, Goemaere S, De Clercq L, et al. The evolution of spondyloarthropathies in relation to gut histology. III. Relation between gut and joint. J Rheumatol. 1995;22:2279–84.

    CAS  PubMed  Google Scholar 

  27. Olofsson T, Lindqvist E, Mogard E, Andreasson K, Marsal J, Geijer M, et al. Elevated faecal calprotectin is linked to worse disease status in axial spondyloarthritis: results from the SPARTAKUS cohort. Rheumatology (Oxford). 2019;58:1176–87.

  28. Klingberg E, Strid H, Stahl A, Deminger A, Carlsten H, Ohman L, et al. A longitudinal study of fecal calprotectin and the development of inflammatory bowel disease in ankylosing spondylitis. Arthritis Res Ther. 2017;19:21.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Rudwaleit M, van der Heijde D, Landewe R, Listing J, Akkoc N, Brandt J, et al. The development of Assessment of SpondyloArthritis international Society classification criteria for axial spondyloarthritis (part II): validation and final selection. Ann Rheum Dis. 2009;68:777–83.

    Article  CAS  PubMed  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  Google Scholar 

  31. Casen C, Vebo HC, Sekelja M, Hegge FT, Karlsson MK, Ciemniejewska E, et al. Deviations in human gut microbiota: a novel diagnostic test for determining dysbiosis in patients with IBS or IBD. Aliment Pharmacol Ther. 2015;42:71–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Savin Z, Kivity S, Yonath H, Yehuda S. Smoking and the intestinal microbiome. Arch Microbiol. 2018;200:677–84.

    Article  CAS  PubMed  Google Scholar 

  33. Turnbaugh PJ, Ley RE, Mahowald MA, Magrini V, Mardis ER, Gordon JI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. 2006;444:1027–31.

    Article  PubMed  Google Scholar 

  34. Dougados M, Simon P, Braun J, Burgos-Vargas R, Maksymowych WP, Sieper J, et al. ASAS recommendations for collecting, analysing and reporting NSAID intake in clinical trials/epidemiological studies in axial spondyloarthritis. Ann Rheum Dis. 2011;70:249–51.

    Article  CAS  PubMed  Google Scholar 

  35. Drossman DA. The functional gastrointestinal disorders and the Rome III process. Gastroenterology. 2006;130:1377–90.

    Article  PubMed  Google Scholar 

  36. Bazin T, Hooks KB, Barnetche T, Truchetet ME, Enaud R, Richez C, et al. Microbiota composition may predict anti-Tnf alpha response in spondyloarthritis patients: an exploratory study. Sci Rep. 2018;8:5446.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Maseda D, Ricciotti E. NSAID-gut microbiota interactions. Front Pharmacol. 2020;11:1153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wallman JK, Mogard E, Marsal J, Andreasson K, Joud A, Geijer M, et al. Irritable bowel syndrome symptoms in axial spondyloarthritis more common than among healthy controls: is it an overlooked comorbidity? Ann Rheum Dis. 2020;79:159–61.

  39. Chong PP, Chin VK, Looi CY, Wong WF, Madhavan P, Yong VC. The microbiome and irritable bowel syndrome - a review on the pathophysiology, current research and future therapy. Front Microbiol. 2019;10:1136.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Breban M, Tap J, Leboime A, Said-Nahal R, Langella P, Chiocchia G, et al. Faecal microbiota study reveals specific dysbiosis in spondyloarthritis. Ann Rheum Dis. 2017;76:1614–22.

    Article  CAS  PubMed  Google Scholar 

  41. Ciccia F, Guggino G, Rizzo A, Alessandro R, Luchetti MM, Milling S, et al. Dysbiosis and zonulin upregulation alter gut epithelial and vascular barriers in patients with ankylosing spondylitis. Ann Rheum Dis. 2017;76:1123–32.

    Article  CAS  PubMed  Google Scholar 

  42. Asquith M, Rosenbaum JT. The interaction between host genetics and the microbiome in the pathogenesis of spondyloarthropathies. Curr Opin Rheumatol. 2016;28:405–12.

    Article  CAS  PubMed  Google Scholar 

  43. Babaie F, Hasankhani M, Mohammadi H, Safarzadeh E, Rezaiemanesh A, Salimi R, et al. The role of gut microbiota and IL-23/IL-17 pathway in ankylosing spondylitis immunopathogenesis: new insights and updates. Immunol Lett. 2018;196:52–62.

    Article  CAS  PubMed  Google Scholar 

  44. Gracey E, Dumas E, Yerushalmi M, Qaiyum Z, Inman RD, Elewaut D. The ties that bind: skin, gut and spondyloarthritis. Curr Opin Rheumatol. 2019;31:62–9.

    Article  PubMed  Google Scholar 

  45. Jung ES, Lee SP, Kae SH, Kim JH, Kim HS, Jang HJ. Diagnostic accuracy of fecal calprotectin for the detection of small bowel Crohn’s disease through capsule endoscopy: an updated meta-analysis and systematic review. Gut Liver. 2021;15:732–41.

    Article  PubMed  Google Scholar 

  46. Cypers H, Varkas G, Beeckman S, Debusschere K, Vogl T, Roth J, et al. Elevated calprotectin levels reveal bowel inflammation in spondyloarthritis. Ann Rheum Dis. 2016;75:1357–62.

    Article  CAS  PubMed  Google Scholar 

  47. Ostgard RD, Deleuran BW, Dam MY, Hansen IT, Jurik AG, Glerup H. Faecal calprotectin detects subclinical bowel inflammation and may predict treatment response in spondyloarthritis. Scand J Rheumatol. 2018;47:48–55.

  48. Kopylov U, Starr M, Watts C, Dionne S, Girardin M, Seidman EG. Detection of Crohn disease in patients with spondyloarthropathy: the SpACE capsule study. J Rheumatol. 2018;45:498–505.

    Article  CAS  PubMed  Google Scholar 

  49. Macke L, Schulz C, Koletzko L, Malfertheiner P. Systematic review: the effects of proton pump inhibitors on the microbiome of the digestive tract-evidence from next-generation sequencing studies. Aliment Pharmacol Ther. 2020;51:505–26.

    Article  PubMed  Google Scholar 

  50. Ringel Y, Maharshak N, Ringel-Kulka T, Wolber EA, Sartor RB, Carroll IM. High throughput sequencing reveals distinct microbial populations within the mucosal and luminal niches in healthy individuals. Gut Microbes. 2015;6:173–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Carstens A, Roos A, Andreasson A, Magnuson A, Agreus L, Halfvarson J, et al. Differential clustering of fecal and mucosa-associated microbiota in ‘healthy’ individuals. J Dig Dis. 2018;19:745–52.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We are indebted to all patients, controls, and staff involved in the SPARTAKUS study. A particular thanks to our research nurse Miriam Walsh Ingelström for study coordination.

Funding

Open access funding provided by Lund University. This study was supported by unrestricted grants from Skåne University Hospital, the Swedish Rheumatism Association, the Anna-Greta Crafoord Foundation, the Kock Foundation, The Österlund Foundation, the Nanna Svartz Foundation, and the Lundgren Foundation. Funding from the Hedlund Foundation and the Österlund Foundation contributed to financing JM’s research time. TO, JM, KA, and JKW received grants from the Swedish government to researchers in public health care (ALF). The sponsors had no role in study design, data collection, data analysis, data interpretation, or writing of the report.

Author information

Authors and Affiliations

  1. Section of Rheumatology, Department of Clinical Sciences Lund, Lund University, Kioskgatan 5, 22185, Lund, Sweden

    Jonas Sagard, Tor Olofsson, Elisabeth Mogard, Kristofer Andréasson, Lars Erik Kristensen, Elisabet Lindqvist & Johan K. Wallman

  2. Skåne University Hospital, Department of Rheumatology, Lund, Sweden

    Jonas Sagard, Tor Olofsson, Elisabeth Mogard, Kristofer Andréasson, Elisabet Lindqvist & Johan K. Wallman

  3. Department of Immunology, EMV, Lund University, Lund, Sweden

    Jan Marsal

  4. Department of Gastroenterology, Skåne University Hospital, Lund/Malmö, Sweden

    Jan Marsal

  5. Sahlgrenska Academy, Institute of Clinical Sciences, Department of Radiology, University of Gothenburg, Gothenburg, Sweden

    Mats Geijer

  6. Region Västra Götaland, Department of Radiology, Sahlgrenska University Hospital, Gothenburg, Sweden

    Mats Geijer

  7. Section of Radiology, Department of Clinical Sciences Lund, Lund University, Lund, Sweden

    Mats Geijer

  8. Parker Institute, Frederiksberg and Bispebjerg, Department of Rheumatology, Copenhagen University Hospital, Copenhagen, Denmark

    Lars Erik Kristensen

Authors
  1. Jonas Sagard
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Tor Olofsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Elisabeth Mogard
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Jan Marsal
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Kristofer Andréasson
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Mats Geijer
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Lars Erik Kristensen
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Elisabet Lindqvist
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Johan K. Wallman
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

JS participated in the study design, analysis, and interpretation of the data and draft and revision of the manuscript. JKW participated in the study design, acquisition of data, analysis and interpretation of data, and draft and revision of the manuscript. TO, EM, and EL participated in the study design, acquisition of data, interpretation of data, and revision of the manuscript (i.e., revising it critically for important intellectual content). KA, JM, MG, and LEK participated in the study design, interpretation of data, and revision of the manuscript (i.e., revising it critically for important intellectual content). The authors read and approved the final manuscript.

Corresponding author

Correspondence to Jonas Sagard.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for the SPARTAKUS study was granted by the Regional Ethics Committee in Lund, Sweden (Dnr. 2015/436 with amendment Dnr. 2018/238). Oral and written informed consent was acquired from all patients and controls before entry into the study, which was performed in accordance with the Helsinki Declaration.

Consent for publication

Not applicable.

Competing interests

JS: None declared. TO: Consultant of Eli Lilly, Merck Sharp & Dohme. EM: Consultant of Novartis. JM: Speakers bureau of AbbVie, Bristol-Myers Squibb, Ferring, Janssen Pharmaceuticals, Pfizer, Takeda, Tillotts; Consultant of AbbVie, Janssen Pharmaceuticals, Pfizer, Takeda; Grant/research support from AbbVie, Ferring, Pfizer, Takeda. KA: None declared. MG: Speakers bureau of AbbVie, Novartis, Pfizer, UCB Pharma. LEK: Speakers bureau of AbbVie, Amgen, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen Pharmaceuticals, Merck Sharp & Dohme, Novartis, Pfizer, UCB Pharma; Consultant of AbbVie, Amgen, Bristol-Myers Squibb, Celgene, Eli Lilly, Janssen Pharmaceuticals, Merck Sharp & Dohme, Novartis, Pfizer, UCB Pharma; EL: None declared. JKW: Consultant of AbbVie, Amgen, Celgene, Eli Lilly, Novartis.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1.

Exploratory analyses of bacterial abundance in axial spondyloarthritis patients versus controls, and sensitivity analyses regarding disease activity, physical function, mobility and pain in relation to gut dysbiosis status in the patient group. Exploratory analyses comparing probe signal intensity for the 48 bacterial markers of the GA-map™ Dysbiosis Test between the axial spondyloarthritis patients and controls (including Figure S1). Figure S2. Differences in measures of disease activity, physical function, mobility and pain between axial spondyloarthritis patients with gut dysbiosis versus those without dysbiosis, when only including patients without gut inflammation. Figure S3. Differences in measures of disease activity, physical function, mobility and pain between axial spondyloarthritis patients with gut dysbiosis versus those without dysbiosis, when only including patients without irritable bowel syndrome symptoms.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sagard, J., Olofsson, T., Mogard, E. et al. Gut dysbiosis associated with worse disease activity and physical function in axial spondyloarthritis. Arthritis Res Ther 24, 42 (2022). https://0-doi-org.brum.beds.ac.uk/10.1186/s13075-022-02733-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13075-022-02733-w

Keywords

Arthritis Research & Therapy

ISSN: 1478-6362

Contact us