Article Text

Download PDFPDF

Original research
Association of body mass index on disease activity in axial spondyloarthritis: systematic review and meta-analysis
  1. Jean W Liew1,
  2. Irvin J Huang1,
  3. Diana N Louden2,
  4. Namrata Singh1 and
  5. Lianne S Gensler3
  1. 1Medicine, Division of Rheumatology, University of Washington, Seattle, Washington, USA
  2. 2University Libraries, University of Washington, Seattle, Washington, USA
  3. 3Division of Rheumatology, University of California San Francisco, San Francisco, California, USA
  1. Correspondence to Jean W Liew; liew.jw{at}gmail.com

Abstract

Objectives In axial spondyloarthritis (axSpA), higher body mass index (BMI) is associated with worse outcomes including response to biologics. Further clarity is needed on whether BMI is associated with disease activity overall, independent of treatment response. We performed a systematic review and meta-analysis to assess the association between BMI and disease activity as reported by the Bath Ankylosing Spondylitis Disease Activity Index (BASDAI) or Ankylosing Spondylitis Disease Activity Score (ASDAS) in axSpA.

Methods We systematically searched for studies evaluating BMI and disease activity as the exposure and outcome of interest, respectively, in axSpA. Using random effects models, we estimated summary standardised mean differences (SMDs) and 95% CIs of BASDAI or ASDAS, comparing obese (BMI>30 kg/m2) or overweight/obese (BMI>25 kg/m2) individuals to those with normal BMI (18.5–24.9 kg/m2).

Results Twelve studies were included in the meta-analysis. Among all studies reporting the BASDAI at baseline, the pooled SMD of the BASDAI for those with an obese or overweight/obese BMI compared to a normal BMI was 0.38 (95% CI 0.21 to 0.55, I2 =75.2%), indicating a significant association of higher BMI with higher BASDAI score. The pooled SMD of the ASDAS for those with an obese or overweight/obese BMI compared to a normal BMI was 0.40 (95% CI 0.27 to 0.54, I2=0%). Findings were robust across subgroup analyses.

Conclusion These results demonstrate an association between an overweight/obese BMI and higher disease activity in studies of axSpA. Future longitudinal studies of BMI and disease activity should assess how this association changes over time.

  • Ankylosing spondylitis
  • Disease activity
  • Epidemiology
  • Outcomes research
  • Adalimumab
  • Rheumatoid arthritis
  • Autoimmune diseases
  • Spondyloarthritis
  • Arthritis
http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

INTRODUCTION

Axial spondyloarthritis (axSpA) is a chronic inflammatory arthritis that affects the spine and sacroiliac joints and can be separated into ankylosing spondylitis (AS), which is also known as radiographic axSpA, and nonradiographic axSpA (nr-axSpA). Cardiovascular (CV) disease is the leading cause of death worldwide, and multiple population-based studies have demonstrated increased CV events and CV-related mortality in axSpA.1–5 Traditional CV risk factors in the general population include diabetes mellitus, hypertension, dyslipidemia, tobacco smoking and obesity.

Although obesity is an important and modifiable CV risk factor, there have been limited studies addressing the impact of higher body mass index (BMI), or of being overweight or obese, on clinical outcomes in axSpA. In a US-based registry of patients with AS and psoriatic arthritis (PsA), obesity was a significant predictor of tumour necrosis factor inhibitor (TNFi) switching or discontinuation.6 Prior systematic reviews and meta-analyses have assessed the effect of BMI on TNFi response in multiple inflammatory diseases with the finding that higher BMI was associated with increased odds of an inadequate response to TNFi treatment in individuals with axSpA.7 8 However, these studies have been limited by the inclusion of only patients initiating biological therapy.

A knowledge gap remains regarding whether higher BMI is associated with higher disease activity overall in axSpA, exclusive of treatment response to biologics. The aim of this study is to conduct a systematic literature review and meta-analysis of observational studies and randomised clinical trials (RCTs) in individuals with AS or nr-axSpA assessing the association of BMI on disease activity as measured by the Bath AS Disease Activity Index (BASDAI) or AS Disease Activity Score (ASDAS).

METHODS

Literature search and study selection

We conducted this systematic review according to an established protocol based on the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement and reporting recommendations from the Meta-analysis of Observational Studies in Epidemiology (MOOSE) guidelines.9 10 The PRISMA checklist is included as a Supplemental file.

With input from the clinical investigators, an experienced medical librarian (DNL) performed an online database search for studies examining the association of BMI with disease activity among individuals with AS or nr-axSpA in PubMed, Embase, the Cochrane Central Register of Controlled Trials and ClinicalTrials.gov from the date of database inception to 15 December 2019. For the study population and the exposure of interest, we employed the following terms in PubMed: (‘Spondylitis, Ankylosing’[Mesh] OR ‘ankylosing spondylitis’ OR ‘axial spondyloarthritis’ OR ‘ankylosing spondylarthritis’ OR ‘axial spondylitis’ OR axSpA) AND (‘Body Mass Index’[Mesh] OR ‘Overweight’[Mesh] OR ‘Obesity’[Mesh] OR ‘body mass index’ OR BMI OR overweight OR obese OR obesity). Comparable combinations of subject headings and keywords were used in the other databases. Full search strategies are included in Supplementary Materials. We included RCTs and observational studies including cohort, case–control and cross-sectional studies. Conference abstracts were also included. Reference lists from retrieved articles and existing reviews were manually searched for additional studies. The titles and abstracts of each citation retrieved from these searches were reviewed based on the following inclusion criteria:

  1. The main study population was AS, nr-axSpA or all axSpA;

  2. The primary exposure of interest was BMI;

  3. The primary outcome of interest was a validated measure of disease activity, either the BASDAI, ASDAS or a treatment response measure based on the BASDAI or ASDAS.

Studies not fulfilling these inclusion criteria were excluded. Duplicates and non-English-language articles and abstracts were excluded. In the case of multiple publications reporting the results for the same study subjects, we included the data from the most recent report.

Exposure and outcome of interest

The exposure of interest was BMI reported as either a continuous or categorical variable. Categorical BMI was either dichotomised as normal and overweight/obese, or by separate levels for normal, overweight and obese BMI. For Caucasians, the WHO defines obese as a BMI ≥30 kg/m2, overweight as 25.0–29.9 kg/m2, normal as 18.5–24.9 kg/m2 and underweight as <18.5 kg/m2.11

The two most commonly used measures for axSpA in clinical practice are the BASDAI and the ASDAS. The BASDAI was developed in 1994 and comprises six questions addressing five major symptoms in AS: fatigue, spinal pain, peripheral joint pain and swelling, localised tenderness and morning stiffness.12 The ASDAS includes three questions from BASDAI, as well as the patient global assessments, and laboratory measures (either the C reactive protein (CRP) or the erythrocyte sedimentation rate).13 For the main analysis, we considered the BASDAI or ASDAS, as a continuous or categorical variable, measured at baseline as the primary outcome. In the sub-analysis of cohort studies, we considered the outcome of TNFi response as a measure of change in the BASDAI or ASDAS.

Data extraction and quality assessment

The full text of eligible articles was reviewed by two independent reviewers (JWL and IJH). Data from each study were abstracted based on a structured data abstraction form. The decision on whether to include the study in the final qualitative and quantitative analyses was based on an agreement between the two reviewers with disagreements resolved by a third reviewer (NS).

We used the Newcastle–Ottawa scale (NOS) for the assessment of study quality in observational studies and the Cochrane risk of bias tool for RCTs. With the NOS, each study is judged on the selection of study groups, the comparability of study groups and the ascertainment of the exposure or outcome and given a score between 0 and 2.14 A modification of the NOS was used for cross-sectional studies.15 16 The quality score is based on a scale out of 10 for cross-sectional studies and 9 for cohort studies, with a higher score indicating better study quality.

Data synthesis and analysis

We extracted the mean and SD for BASDAI or ASDAS values at baseline, stratified by categories of BMI, from all included studies. If the study only reported the median and IQR, we calculated the mean and SD using the formulae from Hozo et al.17 We calculated the standardised mean difference (SMD) using Hedges’ g18 in either the BASDAI or ASDAS comparing those with either overweight/obese or obese BMI to those with normal BMI. If BMI was a three-level categorical measure, we compared the highest, or obese, category with the normal BMI category as the referent. If BMI was a binary measure, we compared the higher category (overweight/obese) with the lower category (normal). The SMD represents the difference between the weighted mean and the SD of the disease activity measure in individuals with higher versus lower BMI categories.19 Due to study heterogeneity, we used the random effects model described by DerSimonian and Laird to calculate the pooled SMD and 95% CIs across studies.20 With this method, study weights depend partly on the size of individual studies and partly on the amount of variation in the results of studies.

We additionally performed sensitivity analyses to assess the robustness of our findings, in which we restricted analyses to the following more homogeneous subgroups:

  1. Cross-sectional studies;

  2. Studies in which the study population consisted of AS patients only;

  3. Studies in which BMI was modelled as a three-level categorical variable (obese, overweight, normal);

  4. Studies performed in European populations.

For cohort studies, we prespecified a plan to extract the raw data, where available, for the crude ORs (OR) of TNFi response comparing a higher BMI to a lower BMI. In studies where these data were not available, we planned to extract the adjusted OR and 95% CIs.21 All primary outcome measures of TNFi response at follow-up would then be combined using a random effects model to calculate the pooled OR for the failure of TNFi response.

We used Cochran’s Q and the I2 index to examine heterogeneity across studies. Cochran’s Q is a test for heterogeneity where the null hypothesis is that all studies are evaluating the same effect. However, this test has low power if there are few included studies.22 23 The I2 index describes the proportion of variation across studies that is due to heterogeneity rather than random chance. An index of 25% indicates low heterogeneity, 50% moderate and 100% high.24

Publication bias for meta-analyses of 10 or more studies was visually assessed using funnel plots with pseudo-95% CIs, which indicate regions where 95% of studies are expected to lie if all studies are estimating the same effect.25 We used the trim-and-fill method to quantify the impact of potential publication bias.26 We performed Egger’s test for meta-analyses of 10 or more studies to additionally evaluate publication bias.

Analyses were significant at the alpha level of 0.05. Analyses were conducted in R version 3.6.027 using the packages metafor,28 esc29 and MAd.30

RESULTS

We identified 961 titles from our database search. Of these, 223 were duplicates, which we excluded. From the remaining 738 studies, we reviewed the full text of 20 studies. Of the full-text articles, we excluded three for not having the primary exposure or outcome of interest, three for not reporting disease activity stratified by BMI at baseline and one for incomplete reporting of results. We included 13 studies in our qualitative analysis and 12 in our quantitative meta-analysis. No RCTs met our inclusion criteria. One RCT underwent full-text review but was excluded because disease activity was not stratified by BMI at baseline.31 The PRISMA flow diagram is shown in figure 1.

Figure 1

PRISMA flow chart for study inclusion.

PRISMA, Preferred Reporting Items for Systematic Reviews and Meta-Analyses.

Characteristics of included studies and risk of bias

Tables 1 and 2 provide details of included studies and participants. We included eight cross-sectional studies and five cohort studies in the systematic review. All of the cohort studies assessed the association of BMI and response to TNFi. Six studies included only AS patients (by the modified New York criteria) and seven included axSpA patients (by the Assessment in SpondyloArthritis International Society [ASAS] criteria). Study sample sizes ranged from 46 to 684, with a total of 4054 subjects in the meta-analysis. The proportion of male patients ranged from 22% to 93%. Eleven studies were conducted in Europe; one study was conducted in Asia and one in the Middle East. The risk of bias assessment is shown in tables 3 and 4. Based on the NOS, the mean score was 7.8 for cross-sectional studies. All included cohort studies had a score of 8 out of 10. The main source of bias was selection bias.

Table 1

Characteristics of included cross-sectional studies (n=8)

Table 2

Characteristics of included cohort studies (n=5)

Table 3

Risk of bias assessment for included cross-sectional studies

Table 4

Risk of bias assessment for included cohort studies

Relationship between BMI and disease activity at baseline in cross-sectional and cohort studies

There were 12 studies reporting the BASDAI at baseline and five reporting the ASDAS at baseline. The cohort study of 190 Romanian AS patients by Ancuta et al was included in the qualitative analysis but not the meta-analysis due to incomplete data.44 Among all studies reporting the BASDAI at baseline (12 studies, 3864 subjects), the pooled SMD of the BASDAI for those with an obese or overweight/obese BMI compared to a normal BMI was 0.38 (95% CI 0.21 to 0.55; I2=75.2%, Q-test p<0.001), indicating a significant association of higher BMI with higher BASDAI score (figure 2). Among all studies reporting the ASDAS at baseline (5 studies, 1469 subjects), the pooled SMD of the ASDAS for those with an obese or overweight/obese BMI compared to a normal BMI was 0.40 (95% CI 0.27 to 0.54; I2=0%, Q-test p=0.70) (figure 3).

Figure 2

Forest plot for the standardised mean difference in the BASDAI comparing obese or overweight/obese BMI to normal BMI in cross-sectional and cohort studies.

Negative SMDs indicate the association of higher BMI with lower BASDAI; positive SMDs indicate: the association of higher BMI with higher BASDAI.

BASDAI, Bath Ankylosing Spondylitis Disease Activity Index; BMI, body mass index; SMD, standardised mean difference.

Figure 3

Forest plot for the standardised mean difference in the ASDAS comparing obese or overweight/obese BMI to normal BMI in cross-sectional and cohort studies.

Negative SMDs indicate the association of higher BMI with lower ASDAS; positive SMDs indicate the association of higher BMI with higher ASDAS.

ASDAS, Ankylosing Spondylitis Disease Activity Score; BMI, body mass index; SMD, standarised mean difference.

In the sub-analysis restricted to the cross-sectional studies, compared to those with a normal BMI, those with an obese or overweight/obese BMI had a higher disease activity measure, which was significant for both BASDAI (eight studies, 2858 subjects) and ASDAS (three studies, 788 subjects) (pooled SMD 0.38, 95% CI 0.23 to 0.52, I2=56.9%, Q-test p<0.001 for BASDAI; pooled SMD 0.42, 95% CI 0.25 to 0.59, I2 0%, Q-test p=0.76 for ASDAS). Five studies (587 subjects) limited to an AS population reported the BASDAI (pooled SMD 0.46, 95% CI −0.10 to 1.02, I2=89.8%, Q-test p<0.001) and only two reported the ASDAS (meta-analysis not performed). Studies using categorical (three-level) BMI as the exposure included seven studies (2344 subjects) that reported the BASDAI (pooled SMD 0.36, 95% CI 0.10 to 0.62, I2=79.0%, Q-test p<0.001), and three studies (1142 subjects) that reported the ASDAS (pooled SMD 0.42, 95% CI 0.26 to 0.59, I2=0%, Q-test p=0.38). There were 10 studies (3500 subjects) conducted in Europe with BASDAI was the outcome measure (pooled SMD 0.33, 95% CI 0.15 to 0.50, I2=73.4%, Q-test p<0.001). All of the studies reporting ASDAS were performed in Europe, so an additional sub-analysis was not necessary.

Relationship between BMI and TNFi response at follow-up in cohort studies

Five cohort studies were included in the qualitative analysis that assessed as a primary outcome of TNFi response at follow-up comparing BMI categories. In the four studies that reported the results of a multivariable analysis of the association of BMI and TNFi response, all found that a higher BMI was significantly associated with the lower odds of a TNFi response. Two of these studies reported the BASDAI50 response (50% or greater improvement in the BASDAI) as the primary outcome, one reported the ASAS40 response (40% improvement according to ASAS criteria) and one reported BASDAI≤4. We were only able to extract the raw data for two cohort studies and were unable to verify the subgroup comparisons for adjusted ORs in the remainder of studies. Thus, we were unable to pool ORs for the failure of TNFi response in the included cohort studies.

Publication bias

Funnel plots of the SMD plotted against the SE are shown for studies reporting the BASDAI (figure 4). There is asymmetry on visual inspection, which can be attributed to reporting bias for small studies with null or negative results. The trim-and-fill method did not identify any potentially missing studies. However, this procedure does not address other sources of funnel plot asymmetry, such as the heterogeneity of study methods or random chance.45 Egger’s test failed to reject the null hypothesis of symmetry in the funnel plots (p=0.50); however, this test suffers from low power with fewer studies.46 We did not perform funnel plots or Egger’s test for studies reporting the ASDAS as there were too few studies.

Figure 4

Funnel plot assessing heterogeneity in included studies reporting the BASDAI at baseline.

BASDAI, Bath Ankylosing Spondylitis Disease Activity Index.

DISCUSSION

In our main analyses, higher BMI was significantly associated with higher disease activity as measured by both BASDAI and ASDAS. Results were robust across most sensitivity analyses restricted to prespecified study characteristics. These findings suggest that a higher BMI, potentially through increased adiposity, may contribute to the disease burden in axSpA.

Our findings are similar to prior systematic reviews and meta-analyses in axSpA evaluating the impact of obesity or higher BMI on various clinical outcome measures. Lee et al performed a systematic review of obesity and disease activity outcomes across multiple rheumatic diseases.47 In axSpA, they found one study reporting a neutral association and six with a positive association. These findings were limited by the inclusion of mostly cross-sectional studies from a restricted geographical area, as most were performed in Europe. Given the limitations of their review findings, the authors elected to not pool results in a meta-analysis. Singh et al performed a meta-analysis of the impact of BMI and TNFi response in observational cohort studies and RCTs, and found that in axSpA, higher BMI was associated with increased odds of an inadequate response to TNFi treatment (six studies; pooled OR 3.36, 95% CI 1.33 to 8.51, I2=81%)).7 Shan and Zhang performed a similar meta-analysis, with similar findings, and included four studies examining the association of BMI with TNFi response in axSpA.8 Our qualitative analysis of cohort studies is in accordance with these findings. Due to differences in our prespecified inclusion criteria and different study question, we did not update these meta-analyses. Compared to these prior studies, which had limited generalisability beyond those patients initiating biological therapy, we addressed the association of BMI and disease activity among a broader population of axSpA patients.

The impact of obesity on imaging measures has also been explored. Bakirci et al performed a systematic review on the association of BMI and imaging-defined inflammation and damage in SpA including PsA.48 In four studies, higher BMI was associated with new syndesmophyte formation. In one study, higher BMI was also associated with a higher structural damage score by the modified Stoke Ankylosing Spondylitis Spinal Score. There were no studies using MRI that met their inclusion criteria. This review was limited by the small number of studies included. Limited inference can be made regarding the association of BMI with active inflammation as reflected by disease activity.

The association of high BMI with disease activity has also been studied in PsA, with similar findings as in axSpA. In prospective cohorts of PsA patients on TNFi therapy, Di Minno et al have shown that obesity was associated with a lower probability of achieving or maintaining minimal disease activity and that a weight loss of ≥5% in overweight or obese individuals predicts treatment response.49 50 Obesity was also associated with worse response to TNFi in large Scandinavian PsA registries.51

The biological mechanisms that are believed to underlie obesity as a chronic inflammatory state involve the production of pro-inflammatory cytokines by adipose tissue.52 Higher BMI and higher fat mass have been associated with chronic pain in multiple populations.5355 Obesity may be related to disease activity independent of inflammation, such as through mechanical loading and stress.56

Some limitations of our systematic review and meta-analysis should be considered. The primary limitation is heterogeneity, in both study design and methodology. For example, we have included studies with different populations of interest (AS vs nr-axSpA) and studies that have used different measures of the exposure and outcome, as well as measurement of the outcome at different times relative to treatment initiation. To address this, we performed sensitivity analyses restricted to studies with similar characteristics, which supported the robustness of our results. Because the ASDAS calculation includes the CRP level13 and elevations in CRP are correlated with higher BMI in the general population,57 this may limit the inference from our findings. However, the similar associations seen with the BASDAI, which does not depend upon the CRP, further support our findings. Four of our included studies provided comparisons of the baseline CRP between obese and normal BMI individuals34 35 37 40; only one study found a statistically significant difference.34 Most of our studies were cross-sectional, which limits inference as data were drawn from one point in time rather than derived over a period of follow-up. Cohort studies were limited to those in which all patients were starting biological therapy with TNFi; selection bias may be at play as these studies may only include patients with more severe and active disease. Generalisability is limited as most of these studies were conducted in Europe.

CONCLUSION

We conducted a systematic literature review and meta-analysis that included observational studies in individuals with AS or axSpA assessing the association of BMI on disease activity. In our main meta-analysis, higher BMI was significantly associated with higher disease activity as measured by both the BASDAI and the ASDAS. Results were similar across sensitivity analyses. Future work should involve longitudinal studies of BMI and disease activity to assess how this association changes over time. It is also important to parse out what components of BMI and disease activity are important in this relationship, and particularly those factors that are modifiable and can be targeted by specific interventions.

Key messages

What is already known about this subject

  • Obesity is associated with a poorer response to biological therapy in axial spondyloarthritis (axSpA); higher body mass index (BMI) may be associated with higher disease activity overall.

What does this study add

  • This systematic review and meta-analysis demonstrates that higher BMI is significantly associated with higher disease activity in patients with axial spondyloarthritis (axSpA).

How might this impact on clinical practice or future developments

  • Interventions on weight loss among those with overweight or obese BMIs should be considered in addition to directed therapies for axSpA.

REFERENCES

Footnotes

  • Twitter Jean Liew @rheum_cat and Irvin Huang @Irvin_Huang

  • Contributors JWL substantially contributed to the conception and design of the work, and the acquisition, analysis and interpretation of data; involved in drafting the work and revising it critically for important intellectual content; provided final approval of the version to be published. IJH substantially contributed to the acquisition, analysis and interpretation of data for the work; involved in revising the work critically for important intellectual content; provided final approval of the version to be published. DNL substantially contributed to the conception and design of the work, and the acquisition of data for the work; involved in revising the work critically for important intellectual content; provided final approval of the version to be published. NS substantially contributed to the conception and design of the work, and the analysis and interpretation of data; involved in revising the work critically for important intellectual content; provided final approval of the version to be published. LSG substantially contributed to the conception and design of the work, and the interpretation of data; involved in revising the work critically for important intellectual content; provided final approval of the version to be published.

  • Funding JWL: NIH T32 training grant (T32AR007108) and funding from the Assessment of Spondyloarthritis Society (ASAS) and the Spondyloarthritis Research and Treatment Network (SPARTAN); NS: Rheumatology Research Foundation and American Heart Association; LSG: Spondylitis Association of America and the Russel Engelman Rheumatology Research Center at the University of California, San Francisco.

  • Competing interests LSG: Consulting for AbbVie, GSK, Eli Lilly, Novartis, Pfizer and UCB (under $10 000); Grant/research support from UCB, Novartis, Pfizer;

    All others: None.

  • Patient consent for publication Not required.

  • Data sharing statement Data are available upon reasonable request.

  • Provenance and peer review Not commissioned; externally peer reviewed.