Content » Vol 48, Issue 5

Original report

Explanatory factors and predictors of fatigue in persons with rheumatoid arthritis: A longitudinal study

Caroline Feldthusen, PhD1,2, Anna Grimby-Ekman, PhD3, Helena Forsblad-d’Elia, PhD4, Lennart Jacobsson, PhD1 and Kaisa Mannerkorpi, PhD2,5

From the 1Department of Rheumatology and Inflammation Research, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, 2The University of Gothenburg Centre for Person-Centred Care (GPCC), 3Department of Public Health and Community Medicine, Sahlgrenska Academy and University Hospital, University of Gothenburg, 4Department of Public Health and Clinical Medicine, Rheumatology, Umeå University, Umeå and 5Section of Physiotherapy, Institute of Neuroscience and Physiology, Sahlgrenska Academy, University of Gothenburg, Gothenburg, Sweden

OBJECTIVE: To investigate the impact of disease-related aspects on long-term variations in fatigue in persons with rheumatoid arthritis.

DESIGN: Observational longitudinal study.

METHODS: Sixty-five persons with rheumatoid arthritis, age range 20–65 years, were invited to a clinical examination at 4 time-points during the 4 seasons. Outcome measures were: general fatigue rated on visual analogue scale (0–100) and aspects of fatigue assessed by the Bristol Rheumatoid Arthritis Fatigue Multidimensional Questionnaire. Disease-related variables were: disease activity (erythrocyte sedimentation rate), pain threshold (pressure algometer), physical capacity (six-minute walk test), pain (visual analogue scale (0–100)), depressive mood (Hospital Anxiety and Depression scale, depression subscale), personal factors (age, sex, body mass index) and season. Multivariable regression analysis, linear mixed effects models were applied.

RESULTS: The strongest explanatory factors for all fatigue outcomes, when recorded at the same time-point as fatigue, were pain threshold and depressive mood. Self-reported pain was an explanatory factor for physical aspects of fatigue and body mass index contributed to explaining the consequences of fatigue on everyday living. For predicting later fatigue pain threshold and depressive mood were the strongest predictors.

CONCLUSION: Pain threshold and depressive mood were the most important factors for fatigue in persons with rheumatoid arthritis.

Key words: fatigue; rheumatoid arthritis; longitudinal study; outcome assessment.

J Rehabil Med 2016; 48: 469–476

Correspondence address: Caroline Feldthusen, Department of Rheumatology and Inflammation Research, Box 480, SE-405 30 Gothenburg, Sweden. E-mail: caroline.feldthusen@vgregion.se

Accepted Feb 25, 2016; Epub ahead of print Apr 15, 2016

INTRODUCTION

Pharmacological treatment of rheumatoid arthritis (RA) has improved substantially over the last decade, contributing to better control of inflammation, reduced joint damage and improved maintenance of function (1). However, despite adequate pharmacological treatment, many persons with RA report that they experience fatigue. This has been noticed increasingly and has become established as one of the most important issues in rheumatology, from the perspective of the patients (2). After pain, fatigue is the most prominent symptom of the disease, and for many patients it has greater impact on daily life than pain (3).

The prevalence of fatigue in persons with RA varies in different investigations, from 42% to 80%, and the prevalence of severe fatigue is approximately 50% (4). Persons with RA describe fatigue as multidimensional with physical, cognitive and emotional components (3, 5).

Although fatigue is commonly reported in RA, little is known about its causes and consequences (6). A conceptual model for understanding fatigue in RA proposes that fatigue is influenced by an interaction between disease processes, thoughts, feelings and behaviours, and personal life issues (7). A subsequent review suggested that factors contributing to fatigue would be found among disease-related aspects, such as pain, inflammatory activity, physical functioning, cognitive/emotional functioning and social/environmental aspects (6). However, none of these variables show profound and stable relationships with fatigue across studies. For example, 5 of 6 studies reporting clinical characteristics of inflammatory activity (e.g. erythrocyte sedimentation rate (ESR), C-reactive protein (CRP)) showed no association with fatigue (6). The strongest associations with fatigue were found with pain, disability/physical functioning and depression/depressive mood (6).

Pain is a core symptom in RA and many patients continue to have pain when inflammation has declined (8). At joint site, pain threshold, i.e. the degree of pressure eliciting pain, is found to be lower in persons with RA than in healthy controls, and sleep problems are found to be related to pain and pain threshold (9, 10).

Most studies investigating associations with fatigue have been cross-sectional and have applied a global rating scale, such as a one-dimensional visual analogue scale (VAS), for measuring fatigue (6). Due to the multidimensional nature of fatigue, several aspects, such as physical and mental aspects, are recommended for evaluating fatigue (6). It appears that measures of one-dimensional or multiple aspects of fatigue are useful in different contexts (11).Therefore, longitudinal studies that measure both one-dimensional and multiple aspects of fatigue (6) regularly over time (4) are necessary to yield more knowledge about fatigue. To the best of our knowledge no such study has been published.

The aim of this study was to investigate how disease-related aspects contribute to the variation in general fatigue, as well as in multiple aspects of fatigue, in persons with RA at 4 time-points during the 4 seasons, using multivariable regression analysis.

METHODS

Participants

Participants were identified by the hospital administrative register at the rheumatology clinic, Sahlgrenska University Hospital, Gothenburg, Sweden.

Inclusion criteria were: diagnosis of RA according to the International Classification of Diseases 10 (diagnosis codes M05 and M06) (12); of working age (20–65 years); disease duration > 3 years; and stable pharmacological treatment > 3 months prior to study entry concerning disease modifying anti-rheumatic drugs (DMARDs) including biological DMARDS and glucocorticosteroids. Exclusion criteria were: other severe somatic or psychiatric diseases or not having the capacity to communicate effectively in Swedish.

Persons with RA, with a diagnosis code, fulfilling the age criterion and with a visit to the clinic in the last 1.5 years were recruited. Of these, 250 persons were randomly selected and inclusion and exclusion criteria were checked from the medical records. If eligible, a letter with an invitation to participate in the study was sent to the persons (n = 140). A total of 72 persons (51%) answered the invitation letter and, of those, 65 persons entered the study (7 out of 72 declined to participate because of: other health problems n = 3, no time n = 2, family reasons n = 1, to much effort n = 1). Three persons, 65 years of age at the time of recruitment had turned 66 years when entering the study. Ninety-one percent (n = 61) of subjects participated in > 50% of the measurements.

Measures

Demographic data, including age, sex, disease duration, medication, work status and body mass index (BMI), were obtained by a standardized interview.

Clinical examination included:

• Disease activity, assessed by disease activity score, DAS-28 (13) based on the number of tender and swollen joints (0–28), patient global assessment of general health (VAS 0–100 mm) and ESR.

• Pain threshold, assessed using a pressure algometer (Somedic Sales AB, Hörby, Sweden) with 1 cm2 probe area and a pressure increase of approximately 50 kPa/s (14). Pain thresholds were measured in kPa and assessed bilaterally in the upper and lower limb (thumb nail base; m. trepezius; m. supraspinatus; knee m. vastus medialis; and m. gluteus). A mean value of the 10 locations assessed was calculated.

• Physical capacity assessed by the six-minute walk test (6MWT), measuring the walking distance, in m, during 6 min (15).

Self-reported questionnaires included:

• General fatigue rated on a one-dimensional VAS (0–100 mm) during the last week, with anchors “no fatigue” and “worst imaginable fatigue”.

• Multiple aspects of fatigue rated on the Bristol Rheumatoid Arthritis Fatigue–Multi-dimensional Questionnaire (BRAF-MDQ) (Swedish version) (16, 17). The summa-score “Total” measures general fatigue (0–70) and the 4 sub-scores are “Physical” (0–22; a measure of physical fatigue), “Living” (0–21; describing sequelae due to the unpredictability of fatigue), “Cognition” (0–15; describing the cognitive effects of fatigue) and “Emotion” (0–12; describing the effects of fatigue on emotions and mood). A higher score denotes more severe fatigue.

• Self-reported pain rated on VAS (0–100 mm) during last week, with anchors “no pain” and “worst imaginable pain”.

• Depressive mood rated on the Hospital Anxiety and Depression scale, depression subscale HADS-D (0–21). A score ≥ 11 indicates probable depression and a score ≥ 8 indicates possible depression (18).

• Functional status (disability) rated on the Swedish version of Health Assessment Questionnaire (HAQ) range 0–3. A score ≤ 1.0 is regarded as indicating a low level of disability (19).

Procedure

All participants were invited to a clinical examination at 4 time-points during the course of the study, approximately 3 months apart, to cover the 4 seasons. At these time-points demographic data, disease activity, inflammation parameters, pain thresholds and physical capacity were assessed and self-reported questionnaires to measure the level of fatigue and other health aspects in RA were administered.

The primary objective of this study was to analyse whether disease activity, physical capacity, pain, pain thresholds, depressive mood, age, sex, body mass index (BMI) and season, explain variations in fatigue over time. Variables recorded at the same time-point as the fatigue scorings of the participants were used as explanatory factors of fatigue. Variables measured at the previous time-point, approximately 3 months previously, were used as predictors of fatigue.

The study was approved by the Regional Ethical Review Board in Gothenburg. Written and oral information were given to all participants and written consent was obtained from all participants.

Trial registration. ClinicalTrials.gov Identifier NCT01697202.

Statistical analysis

Descriptive data are presented as mean ± standard deviation (SD), median (max; min) or number (n) and percentage for demographic data.

Correlation analyses using Pearson’s correlation coefficient were made for the independent variables to test for risk of collinearity (r ≥ 0.7). Due to risk of collinearity, disease activity measured with DAS-28 was omitted from the regression analysis in favour of ESR.

Multivariable regression analysis, the linear mixed effects model, was used to analyse the outcome fatigue, with the 4 time-points included as an explanatory variable (factor). The regression model was specified to have a random intercept to handle the repeated measurements over time (20). In the linear mixed effect model all measured values are used and missing values do not lead to exclusion of any measured data. The selected fixed variables were: ESR, VAS pain, Algometer, 6MWT and HADS-D. Also the personal factors Age, BMI, and Sex were included in the analyses as well as Season. Sex and Season were included as factors, Sex was defined as female and male, and Season was defined as winter (December, January, February), spring (March, April, May), summer (June, July, August) and autumn (September, October, November). The other variables were included as continuous variables.

Two analyses were tested: (i) explanatory factors for fatigue with variables and outcome recorded at the same time-point; (ii) predictors of fatigue with variables recorded at the time-point prior to the outcome (approximately 3 months). In both analyses, longitudinal data with measures at repeated time-points were analysed according to the model building procedure described below.

To identify important variables a process of model building steps were performed similar to the steps proposed earlier for logistic regression (21). As a first step, to determine which variables to include in the multivariable regression analysis full model, univariate analyses were performed with each independent variable one at a time. Variables with p < 0.25 using F-test were included in the multivariable regression analysis full model. In the multivariable analysis full model variables with p ≥ 0.25 using t-test were excluded, resulting in a reduced model. To determine which variables to include in a best-fitting final model the reduced model was compared with the full model and the excluded variables from the full model and the variables with high p-value in the univariate regression analyses were reinserted one by one in the reduced model and tested once again. This model building procedure was performed separately for all the outcome variables (VAS fatigue, BRAF-MDQ total score and the 4 subscores).

Estimates with 95% confidence intervals, are presented for the variables included in the final multivariable models.

The variance explained by the final models is presented as 2 ratios; 1 for the residual variance and 1 for the between-individual variance (22). The explained residual variance, denoted R2resid, is the ratio between the difference between the residual variance in the final model and the model only including random intercept, divided by the residual variance in the model only including random intercept. The explained between-individual variance, denoted R2between-ind, is the difference between the between-individual variance in the final model and the model only including random intercept, divided by the between-individual variance in the model only including random intercept.

Statistical analyses was performed using SPSS version 15.0 (SPSS Inc., IBM, Chicago USA).

RESULTS

Descriptive statistics of the baseline characteristics of the 65 persons participating in the study are shown in Table I. The mean age and standard deviation (SD) of the participants was 54 years (SD 9.9) , with disease duration of 15 years (SD 9.6), and 74% were women. Baseline fatigue ranged from no fatigue to maximum fatigue (see Table II). Medication at baseline is shown in Table I. During the course of the study 6 participants increased the dosage of, or added, a DMARD, 4 participants changed to another DMARD, and 10 participants decreased the dosage or discontinued a DMARD. Six participants received an intra-articular injection of glucocorticosteroids and 5 participants received short-term use of oral glucocorticosteroids during the investigation. The number of missing values ranges from 0 to 10 in different analyses. No trends were found in missing values in terms of age and symptom severity.

Table I. Demographic data of the 65 persons with rheumatoid arthritis at baseline

Variables

Mean (SD)

Median (min;max)

n (%)

Age, years

53.7 (9.9)

56 (23;66)

 

Sex, female

 

 

48 (74)

Disease duration, years

15.3 (9.6)

12 (4;45)

 

BMI, kg/m2

26.7 (5.2)

26 (19;42)

 

DAS 28

3.7 (1.4)

3.8 (0.8;6.9)

 

Tender joints 0–28

6.9 (6.5)

5 (0;27)

 

Swollen joints 0–28

3.7 (3.5)

3 (0;13)

 

Erythrocyte sedimentation rate, mm

10.8 (10.2)

7 (2;53)

 

General health, VAS 0–100, mm

36.4 (22.7)

32 (1;83)

 

HAQ 0–3

0.6 (0.6)

0.6 (0;2.4)

 

Medication

 

 

 

 

No DMARD

 

 

8 (12)

 

Conventional synthetic DMARD

 

57 (88)

 

Biological DMARD

 

 

15 (23)

 

NSAID

 

 

43 (66)

 

Painkiller

 

 

25 (38)

 

Antidepressive drug

 

 

2 (3)

Work status

 

 

 

 

Working or studying (full-time or part-time)

 

41 (63)

 

Unemployed

 

 

7 (11)

 

Retired

 

 

6 (9)

 

Disability benefits (full-time or part-time)

 

24 (37)

 

Parent’s allowance

 

2 (3)

VAS: visual analogue scale; BMI: body mass index; DAS: Disease Activity Score; HAQ: Health Assessment Questionnaire; DMARD: disease-modifying anti rheumatic drug; NSAID: non-steroidal anti-inflammatory drugs.

Table II. Descriptive data of the outcome variable of fatigue, assessed with VAS fatigue and BRAF- MDQ and key explanatory variables, in the 65 persons with RA at the 4 clinical examinations

 

1

 

2

 

3

 

4

Mean (SD)

Median (min;max)

 

Mean (SD)

Median (min;max)

 

Mean (SD)

Median (min;max)

 

Mean (SD)

Median (min;max)

Outcome variables

 

 

 

 

 

 

 

 

 

 

 

VAS fatigue (0–100)

47.1 (27.1)

49 (0;100)

 

49.8 (26.4)

48 (2;98)

 

50.6 (28.0)

53 (0;98)

 

46.3 (27.9)

44 (1;94)

BRAF-MDQ

 

 

 

 

 

 

 

 

 

 

 

Total (0–70)

29.0 (14.8)

29 (0;65)

 

27.0 (12.8)

26 (0;56)

 

26.2 (13.8)

26 (0;51)

 

24.3 (14.2)

27 (0;56)

Physical (0–22)

13.9 (5.1)

15 (0;21)

 

13.8 (4.7)

15 (0;22)

 

13.6 (5.1)

15 (0;20)

 

12.5 (5.6)

14 (0;21)

Living (0–21)

6.1 (4.5)

6 (0;21)

 

5.4 (4.2)

4 (0;16)

 

5.0 (4.0)

5 (0;14)

 

4.6 (4.0)

4 (0;16)

Cognitive (0–15)

5.3 (4.0)

5 (0;15)

 

4.6 (3.6)

4 (0;15)

 

4.6 (3.5)

5 (0;13)

 

4.3 (3.5)

4 (0;14)

Emotion (0–12)

3.6 (3.1)

3 (0;12)

 

3.2 (2.8)

3 (0;12)

 

2.9 (3.1)

2 (0;11)

 

2.9 (3.1)

3 (0;11)

Explanatory variables

 

 

 

 

 

 

 

 

 

 

 

ESR (mm)

10.8 (10.2)

7 (2;53)

 

11.3 (10.1)

8 (2;48)

 

11.5 (11.1)

8 (2;56)

 

11.6 (10.1)

10 (1;50)

VAS Pain (0–100)

38.0 (25.5)

34 (0;100)

 

34.4 (24.8)

30 (1;100)

 

37.2 (25.2)

35 (0;84)

 

36.5 (27.7)

32 (0;97)

Algometer, kPa

301.2 (119.0))

281 (95;650)

 

301.0 (114.8)

272 (153;593)

 

304.8 (121.9)

278 (85;662)

 

328.4 (130.9)

283 (103;663)

6MWT, m

535.2 (98.1)

538 (232;775)

 

544.9 (84.8)

544 (273;704)

 

537.1 (86.8)

539 (302;727)

 

537.6 (87.0)

531 (316;699)

HADS-D (7–21)

4.8 (3.8)

4 (0;16)

 

4.4 (3.6)

3 (0;15)

 

4.2 (3.5)

4 (0;14)

 

4.2 (3.5)

3 (0;14)

VAS: visual analogue scale; BRAF-MDQ: Bristol Rheumatoid Arthritis Fatigue Multidimensional Questionnaire; ESR: erythrocyte sedimentation rate; 6MWT: six-minute walk test; HADS-D: Hospital Anxiety and Depression Scale – Depression subscale; SD: standard deviation.

 

Analyses of explanatory factors of fatigue recorded at the same time-point as the outcome

The result of univariate analyses to select which explanatory factors (p < 0.25) to include in the multivariable regression analysis full model are shown in Table III.

Table III. Explanatory factors: univariate regression analyses to determine which explanatory factors (p < 0.25) to include in the multivariable analysis full model

Explanatory factors

VAS

Fatigue

BRAF-MDQ

Total

Physical

Living

Cognition

Emotion

ESR

p = 0.328

p = 0.880

p = 0.790

p = 0.696

p = 0.491

p = 0.969

VAS Pain

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.033

Algometer

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.000

6MWT

p = 0.000

p = 0.000

p = 0.000

p = 0.001

p = 0.031

p = 0.005

HADS-D

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.000

p = 0.000

Age

p = 0.566

p = 0.473

p = 0.567

p = 0.727

p = 0.273

p = 0.560

Sex

p = 0.063

p = 0.303

p = 0.112

p = 0.450

p = 0.505

p = 0.506

BMI

p = 0.758

p = 0.079

p = 0.177

p = 0.012

p = 0.526

p = 0.081

Season

p = 0.085

p = 0.523

p = 0.171

p = 0.490

p = 0.957

p = 0.842

Significant values are shown in bold.

VAS: visual analogue scale; BRAF-MDQ: Bristol Rheumatoid Arthritis Fatigue – Multidimensional Questionnaire; ESR: erythrocyte sedimentation rate; 6MWT: six-minute walk test; HADS-D: Hospital Anxiety and Depression Scale – Depression subscale; BMI: body mass index.

Multivariable analyses

The final models are shown in Table IV:

Table IV. Explanatory factors: multivariable analysis of explanatory factors of fatigue. The final models for each fatigue outcome are presented

Explanatory factors

 

Estimate

95% CI

VAS Fatigue (n = 65, 60, 58, 56)

Intercept

 

33.594

20.321;46.868

VAS pain

 

0.451

0.337;0.565

Algometer

 

–0.0417

–0.072;–0.012

HADS-D

 

1.457

0.510;2.405

Season

 

 

 

Winter (n = 60)

 

6.68

0.996;12.366

Spring (n = 58)

 

6.61

1.018;12.506

Summer (n = 48)

 

6.60

0.307;12.885

Autumn (n = 73)

 

 

 

BRAF-MDQ Total score (n = 65, 59, 58, 55)

Intercept

 

15.290

4.071;26.509

VAS pain

 

0.0722

0.0262;0.118

Algometer

 

–0.0382

–0.0521;–0.0243

HADS-D

 

1.969

1.533;2.405

BMI

 

0.458

0.0546;0.861

BRAF-MDQ Physical subscale (n = 65, 59, 58, 55)

Intercept

 

9.889

5.305;14.472

VAS pain

 

0.0342

0.0150;0.0534

Algometer

 

–0.0145

–0.0202;–0.00872

HADS-D

 

0.513

0. 334;0.692

BMI

 

0.152

–0.0109;0.315

Season

 

 

 

Winter (n = 60)

 

0.930

0.00797;1.853

Spring (n = 57)

 

0.721

–0.167;1.609

Summer (n = 47)

 

0.711

–0.264;1.686

Autumn (n = 73)

 

 

 

BRAF-MDQ Living subscale (n = 65, 60, 58, 55)

Intercept

 

–0.455

–4.401;3.491

VAS pain

 

0.0177

0.00190;0.0334

Algometer

 

–0.0106

–0.0154;–0.00584

HADS-D

 

0.506

0. 356;0.655

BMI

 

0.238

0.0959;0.380

BRAF-MDQ Cognitive subscale (n = 65, 60, 58, 55)

Intercept

 

4.089

2.476;5.702

VAS pain

 

0.0139

–0.000248;0.0281

Algometer

 

–0.00724

–0.0111;–0.00336

HADS-D

 

0.549

0. 425;0.673

BRAF-MDQ Emotion subscale (n = 65, 60, 58, 55)

Intercept

 

0.733

–1.992;3.458

Algometer

 

–0.00624

–0.00954;–0.00293

HADS-D

 

0.435

0.329;0.540

BMI

 

0.095

–0.00427;0.194

VAS: visual analogue scale; BRAF-MDQ: Bristol Rheumatoid Arthritis Fatigue-Multidimensional Questionnaire; HADS-D: Hospital Anxiety and Depression Scale–Depression subscale; BMI: body mass index.

 

• VAS fatigue: VAS pain, Algometer, 6MWT, HADS-D, Sex and Season were included in the full multivariable model (see Table III). The final model included VAS pain, Algo­meter, HADS-D and Season. VAS pain (p < 0.001), Algometer (p < 0.01) and HADS-D (p < 0.05) were statistically significant.

• BRAF-MDQ Total: VAS pain, Algometer, 6MWT, HADS-D, and BMI were included in the full multivariable model (see Table III). The final model included VAS pain, Algo­meter, HADS-D, and BMI. Algometer (p < 0.001), HADS-D (p < 0.001), VAS pain (p < 0.01) and BMI (p < 0.05) were statistically significant.

• BRAF-MDQ Physical: VAS pain, Algometer, 6MWT, HADS-D, BMI, Sex and Season were included in the full multivariable model (see Table III). The final model included VAS pain, Algometer, HADS-D, BMI and Season. Algometer (p < 0.001), HADS-D (p < 0.001) and VAS pain (p < 0.01) were statistically significant.

• BRAF-MDQ Living: VAS pain, Algometer, 6MWT, HADS-D and BMI were included in the full multivariable model (see Table III). The final model included VAS pain, Algometer, HADS-D and BMI. Algometer (p < 0.001), HADS-D (p < 0.001), BMI (p < 0.01) and VAS pain (p < 0.05) were statistically significant.

• BRAF-MDQ Cognitive: VAS pain, Algometer, 6MWT and HADS-D were included in the full multivariable model (see Table III). The final model included VAS pain, Algo­meter and HADS-D. Algometer (p < 0.001) and HADS-D (p < 0.001) were statistically significant.

• BRAF-MDQ Emotion: VAS pain, Algometer, 6MWT, HADS-D and BMI were included in the full multivariable model (see Table III). The final model included Algometer, HADS-D and BMI. Algometer (p < 0.001) and HADS-D (p < 0.001) were statistically significant.

In the explanatory models the explained residual variances R2resid were 0.11–0.30, with the 2 highest values (0.30 and 0.25) for BRAF-MDQ Total and BRAF-MDQ Physical. The explained between-individual variances R2between-ind were 0.46–0.64.

Analysis of predictors of fatigue with variables recorded at the time-point prior to the outcome (approximately 3 months)

The results from univariate analyses to select which predictors (p < 0.25) to include in the multivariable regression analysis full model are shown in Table V.

Table V. Predictions: univariate regression analyses to determine which predictors (p < 0.25) to include in the multivariable analysis full model

Predictors

VAS Fatigue

BRAF-MDQ

Total

Physical

Living

Cognition

Emotion

ESR

p = 0.593

p = 0.721

p = 0.682

p = 0.257

p = 0.993

p = 0.507

VAS pain

p = 0.273

p = 0.058

p = 0.083

p = 0.060

p = 0.087

p = 0.096

Algometer

p = 0.000

p = 0.001

p = 0.001

p = 0.016

p = 0.002

p = 0.002

6MWT

p = 0.126

p = 0.483

p = 0.237

p = 0.488

p = 0.407

p = 0.848

HADS-D

p = 0.003

p = 0.000

p = 0.006

p = 0.000

p = 0.000

p = 0.000

Age

p = 0.547

p = 0.358

p = 0.492

p = 0.679

p = 0.144

p = 0.463

Sex

p = 0.078

p = 0.114

p = 0.131

p = 0.104

p = 0.279

p = 0.159

BMI

p = 0.970

p = 0.243

p = 0.257

p = 0.124

p = 0.868

p = 0.288

Season

p = 0.157

p = 0.052

p = 0.145

p = 0.122

p = 0.139

p = 0.265

Significant values are shown in bold.

VAS: visual analogue scale; BRAF-MDQ: Bristol Rheumatoid Arthritis Fatigue-Multidimensional Questionnaire; ESR: erythrocyte sedimentation rate; 6MWT: six-minute walk test; HADS-D: Hospital Anxiety and Depression Scale-Depression subscale; BMI: body mass index.

Multivariable analyses

The final models are shown in Table VI:

Table VI. Predictors: multivariable analyses of predictors of fatigue. The final models for each fatigue outcome are presented

Predictors

 

Estimate

95% CI

VAS Fatigue (n = 60, 56, 57)

Intercept

 

63.820

49.636;78.004

Algometer

 

–0.0760

–0.114;–0.0380

HADS-D

 

1.869

0.596;3.141

BRAF-MDQ Total score (n = 59, 56, 56)

Intercept

 

11.014

–3.377;25.406

VAS pain

 

0.0450

–0.0261;0.116

Algometer

 

–0.0261

–0.0449;–0.00737

HADS-D

 

1.528

0.980;2.0758

BMI

 

0.412

–0.00976;0.834

Sex

 

 

 

Women (n = 6,135)

 

4.389

–1.265;10.0422

Men (n = 636)

 

 

 

BRAF-MDQ Physical (n = 59, 56, 56)

Intercept

 

11.376

6.211;16.540

Algometer

 

–0.0138

–0.0210;–0.00660

HADS-D

 

0.335

0.0981;0.571

BMI*

 

0.179

–0.00823;0.367

BRAF-MDQ Living (n = 60, 56, 56)

Intercept

 

–1.222

–6.130;3.686

Algometer

 

–0.00484

–0.0109;0.00117

HADS-D

 

0.430

0.249;0.611

BMI

 

0.153

0.00885;0.298

Sex

 

 

 

Women (n = 6136)

 

2.075

0.133;4.016

Men (n = 636)

 

 

 

BRAF-MDQ Cognitive (n = 60, 56, 56)

Intercept

 

4.846

3.249;6.443

Algometer

 

–0.00757

–0.0119;–0.00328

HADS-D

 

0.445

0.302;0.589

BRAF-MDQ Emotion (n = 60, 56, 56)

Intercept

 

0.673

–2.939;4.284

Algometer

 

–0.00547

–0.00983;–0.00110

HADS-D

 

0.273

0.141;0.406

BMI*

 

0.0734

–0.0346;0.181

Sex

 

 

 

Women (n = 6136)

 

1.038

–0.378;2.454

Men (n = 636)

 

 

 

*Did not enter the model according to univariate analysis, but was tested in the reduced model and was included.

VAS: visual analogue scale; BRAF-MDQ: Bristol Rheumatoid Arthritis Fatigue – Multidimensional Questionnaire; HADS-D: Hospital Anxiety and Depression Scale – Depression subscale; BMI: body mass index.

 

• VAS fatigue: Algometer, 6MWT, HADS-D, Sex and Season were included in the full multivariable model (see Table V). The final model included Algometer and HADS-D. Both Algometer (p < 0.001) and HADS-D (p < 0.01) were statistically significant.

• BRAF-MDQ Total: VAS pain, Algometer, HADS-D, BMI, Sex and Season were included in the full multivariable model (see Table V). The final model included VAS pain, Algometer, HADS-D, BMI and Sex. HADS-D (p < 0.001) and Algometer (p < 0.01) were statistically significant.

• BRAF-MDQ Physical: VAS pain, Algometer, 6MWT, HADS-D, Sex and Season were included in the full multivariable model (see Table V). The final model included Algometer, HADS-D and BMI. Algometer (p < 0.001) and HADS-D (p < 0.01) were statistically significant.

• BRAF-MDQ Living: VAS pain, Algometer, HADS-D, BMI, Sex and Season were included in the full multivariable model (see Table V). The final model included Algometer, HADS-D, BMI and Sex. HADS-D (p < 0.001), BMI (p < 0.05) and Sex (p < 0.05) were statistically significant.

• BRAF-MDQ Cognitive: VAS pain, Algometer, HADS-D, Age and Season were included in the full multivariable model (see Table V). The final model included Algometer and HADS-D. Both HADS-D (p < 0.001) and Algometer (p < 0.01) were statistically significant.

• BRAF-MDQ Emotion: VAS pain, Algometer, HADS-D and Sex were included in the full multivariable model (see Table V). The final model included Algometer, HADS-D, BMI and Sex. HADS-D (p < 0.001) and Algometer (p < 0.05) were statistically significant.

In the predictive models the explained residual variance R2resid ranged from –0.07 to –0.23 and the explained between-individual variances R2 between-ind were 0.44–0.69. The negative values of R2resid are explained by that the most negative values for R2resid were connected to the highest positive values on R2 between-ind.

To investigate whether the results were influenced by baseline fatigue additional analyses of the predictive approach were performed with baseline fatigue inserted as a variable to predict fatigue. For each fatigue outcome measure the baseline fatigue ratings were added to the respective final model, i.e. the outcome VAS fatigue was controlled by baseline VAS fatigue in the final model. The results showed that HADS-D, Algometer and VAS pain remained as predictors of fatigue in 5 of 6 analyses. The exception was HADS-D for BRAF-MDQ Physical and Algometer for BRAF-MDQ Living. BMI and sex did not remain to be included in any of the models. Baseline fatigue significantly predicted later fatigue in all fatigue outcomes (see Table SI1).

At baseline there was a significant correlation (p < 0.001) between all fatigue variables (r = 0.45–0.77). The lowest correlation was seen between VAS fatigue and BRAF-MDQ Emotion (r = 0.45) and the highest correlation was seen between VAS fatigue and BRAF-MDQ Physical (r = 0.77).

DISCUSSION

This longitudinal study of fatigue, recorded at 4 time-points during 4 seasons in persons with RA, showed that pain threshold and depressive mood were the strongest factors explaining the fatigue, when measured at the same time-point as fatigue, and to predict fatigue 3 months later. Since this was found for general fatigue as well as for the other aspects of fatigue studied, pain threshold and depressive mood seem to relate to overall fatigue. Additional analysis found that fatigue at baseline was a predictor for later fatigue, and previous interviews in persons with RA indicate that fatigue itself uses up energy and thereby contributes to persistence of the fatigue (5). However, when controlling for baseline fatigue pain threshold and depressive mood remained important variables in predicting fatigue.

Depression has, in previous cross-sectional studies, been found to be one of the factors with strong association with fatigue in RA (6), but the results of studies analysing depression as a potential predictor of later fatigue are inconsistent. When comparing baseline fatigue with fatigue 1–2 years later, in one study depression was shown to be a statistically significant predictor of fatigue (23), while in 2 other studies depression failed to predict fatigue (24, 25). This implies that depressive mood is an important factor to take into account in the care of persons with RA.

Decreased pain thresholds have been found to be more common in persons with RA compared with a healthy control group, which may be due to both peripheral and central pain sensitization (26). In this study pain threshold was clearly related to all investigated aspects of fatigue in the multivariable regression models, both as an explanatory factor and as a predictor of fatigue.

Sleep problems are common in persons with RA and have been found to be associated with fatigue in these subjects (27). Previous studies show that poor sleep quality is associated with increased pain sensitivity (9, 10) as well as with depression (27), implying that the relationship between decreased pain thresholds, depression and fatigue may be mediated by poor sleep. However, due to cross-sectional study designs in the previous studies, the causal relationships between sleep quality, pain thresholds, depressive mood and fatigue are unclear and investigations with a longitudinal approach are necessary.

Interestingly, we found the ratings of pain on VAS to be an explanatory factor of fatigue when recorded at the same time-point, but not to predict fatigue 3 months later. This is in line with previous results concluding that fatigue and pain have an association, but not that one precedes or causes the other (28). In this study we found VAS pain to strongly explain the physical aspects of fatigue (i.e. BRAF-MDQ Physical and Living) suggesting that pain intensity may influence the energy and the ability to cope, physically and socially. Previous studies show a weak association between pain threshold and VAS pain in persons with rheumatic diseases and chronic pain (29, 30), while stronger association has been observed between pain threshold and depression (30), both of which relate to the central nervous system. This highlights the multifaceted nature of pain in RA, which can be influenced by peripheral and central pain sensitization beyond inflammation (10, 26). Our results confirm this complex relationship between pain threshold, depressive mood and fatigue.

Physical activity has been found to be inversely associated with fatigue (31, 32) and the results of exercise interventions in persons with RA show positive effects on fatigue (33) as well as on pain (33, 34), depression (33) and sleep quality (35). A large international study in persons with RA found significantly higher fatigue levels among physically inactive persons compared with those who were physically active (36). This implies that physical activity and exercise interventions may be a useful tool in planning interventions to diminish fatigue. Furthermore, physical activity is assumed to result in positive effects on pain and depressive mood, factors found to be of importance for fatigue in this study.

Physical capacity measured with 6MWT showed an ambiguous association with fatigue. In the univariate analysis when physical capacity was recorded at the same time-point as fatigue, a significant association was seen between 6MWT and all fatigue outcomes. However, in the multiple regression models, when covariability was accounted for, 6MWT did not remain significant. Physical capacity measured with 6MWT was neither found to be a predictor of fatigue 3 months later in the univariate or multivariable analysis. According to these results there does not seem to be a causal relationship between physical capacity measured with 6MWT and the variation in fatigue. One possible explanation of the association in the univariate analyses, when recorded at same time-point as fatigue, may be that fatigue influences the performance on the 6MWT and not the other way around (32).

BMI appears to be of potential importance when studying fatigue, given that higher BMI was associated with higher ratings of fatigue. A significant impact of BMI was seen in BRAF-MDQ Living, reflecting the impact of fatigue on the ability to perform activities in everyday life. In RA the inflammatory process has been found to lead to a degradation of lean tissue, especially muscle mass. In combination with inactive lifestyle persons with RA are predisposed to reduced muscle mass and accumulation of body fat, leading to increased adiposity despite stable weight (37, 38). Therefore, BMI is questioned as a reliable measure of overweight and obesity among persons with RA; therefore future studies of the association of body composition and fatigue in RA are of interest.

Seasonal variations in symptoms are often described by persons with RA, but clinical studies have so far failed to show such association (39). In this study season was, in the explanatory multivariable model, found to be a factor of potential importance for general fatigue and the physical aspects of fatigue.

In this study no associations were found between ESR and fatigue. This is in line with the results found in 5 of 6 studies, which showed no association between ESR or CRP and fatigue (6). However, other inflammatory biomarkers might be associated with fatigue (40), which warrants further investigation.

Age or sex did not have an impact on fatigue in our study. However, the analysis was limited to the age criteria of 20–65 years, with a low number of participants in the younger ages and although the distribution of women and men was in concordance with a general RA population, the number of men in this study was limited (n = 17). Previous studies have also failed to show a relationship between fatigue and age, and the relationship between sex and fatigue are showing inconsistent associations (6).

In the current study we applied a multidimensional instrument to rate fatigue. Other researchers suggest that fatigue can be assessed as a unidimensional construct (41). In addition, our analyses show that the different aspects of fatigue correlate with each other and with the one-dimensional VAS. However, it appears that some aspects of fatigue may be easier to treat than others, as in a recent treatment study we found most impact on physical fatigue from an intervention comprising physical activity (42).

The variables in the final models were able to explain a considerable ratio of the variation in fatigue, both at the individual level and at the group level. In the predictive models the result was inconsistent. The variables in the final models seemed to mainly handle the between-individual, rather than the residual variation in fatigue. This indicates that the models were able to predict the level of fatigue rather than change in fatigue. The goodness-of-fit measure used here is R2. It has the useful property of providing an absolute value for the goodness-of-fit of a model, and the amount of variance explained can be of biological interest. Information criteria, e.g. Akaike Information Criteria, cannot provide this information. One problem, however, of R2 for mixed-effects models is that it can be defined in several ways. Despite the problems of the R2 chosen here (e.g. decreased or negative R2 values in larger models), common for several definitions of R2 for mixed-effects models, the choice here is based on their relatively intuitive interpretation and simplicity.

Strengths and weaknesses

Strengths of this study were the longitudinal design and the use of multivariable analysis method, enabling us to draw conclusions about possible explanatory factors and predictors of fatigue.

Limitations of this study were the use of parametric methods with ordinal data which limits the possibility to draw conclusions about the magnitude of the estimates. Also, one should be cautious in generalizing the results due to small sample size. In this study, all results of all models are presented, but only the main patterns in the results are interpreted.

Conclusion

In conclusion, pain threshold and depressive mood appear to be important factors for general as well as physical, cognitive and emotional aspects of fatigue in persons with rheumatoid arthritis. Pain intensity appears to be an important factor to discuss when a patient reports physical fatigue.

ACKNOWLEDGEMENTS

This work was supported by the Gothenburg Centre for Person-Centred Care (www.gpcc.gu.se), the Swedish Research Council and the ALF at Sahlgrenska University Hospital. The authors also would like to thank Sahlgrenska University Hospital/Rheumatology where the study was conducted.

14498.png

1http://www.medicaljournals.se/jrm/content/?doi=10.2340/16501977-2090

 

REFERENCES

1. Scott DL, Wolfe F, Huizinga TW. Rheumatoid arthritis. Lancet 2010; 376: 1094–1108.

2. Kirwan JR, Minnock P, Adebajo A, Bresnihan B, Choy E, de Wit M, et al. Patient perspective: fatigue as a recommended patient centered outcome measure in rheumatoid arthritis. J Rheumatol 2007; 34: 1174–1177.

3. Repping-Wuts H, Uitterhoeve R, van Riel P, van Achterberg T. Fatigue as experienced by patients with rheumatoid arthritis (RA): a qualitative study. Int J Nurs Stud 2008; 45: 995–1002.

4. Repping-Wuts H, Fransen J, van Achterberg T, Bleijenberg G, van Riel P. Persistent severe fatigue in patients with rheumatoid arthritis. J Clin Nurs 2007; 16: 377–383.

5. Feldthusen C, Bjork M, Forsblad-d’Elia H, Mannerkorpi K. Perception, consequences, communication, and strategies for handling fatigue in persons with rheumatoid arthritis of working age-a focus group study. Clin Rheumatol 2013; 32: 557–566.

6. Nikolaus S, Bode C, Taal E, van de Laar MA. Fatigue and factors related to fatigue in rheumatoid arthritis: a systematic review. Arthritis Care Res (Hoboken) 2013; 65: 1128–1146.

7. Hewlett S, Chalder T, Choy E, Cramp F, Davis B, Dures E, et al. Fatigue in rheumatoid arthritis: time for a conceptual model. Rheumatology (Oxford) 2011; 50: 1004–1006.

8. Lee YC, Frits ML, Iannaccone CK, Weinblatt ME, Shadick NA, Williams DA, et al. Subgrouping of patients with rheumatoid arthritis based on pain, fatigue, inflammation, and psychosocial factors. Arthritis Rheum 2014; 66: 2006–2014.

9. Lee YC, Lu B, Edwards RR, Wasan AD, Nassikas NJ, Clauw DJ, et al. The role of sleep problems in central pain processing in rheumatoid arthritis. Arthritis Rheum 2013; 65: 59–68.

10. Lee YC, Chibnik LB, Lu B, Wasan AD, Edwards RR, Fossel AH, et al. The relationship between disease activity, sleep, psychiatric distress and pain sensitivity in rheumatoid arthritis: a cross-sectional study. Arthritis Res Ther 2009; 11: R160.

11. Ericsson A, Bremell T, Mannerkorpi K. Usefulness of multiple dimensions of fatigue in fibromyalgia. J Rehabil Med 2013; 45: 685–693.

12. Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum 1988; 31: 315–324.

13. Fuchs HA. The use of the disease activity score in the analysis of clinical trials in rheumatoid arthritis. J Rheumatol 1993; 20: 1863–1866.

14. Kosek E, Ekholm J, Hansson P. Sensory dysfunction in fibromyalgia patients with implications for pathogenic mechanisms. Pain 1996; 68: 375–383.

15. Mannerkorpi K, Svantesson U, Carlsson J, Ekdahl C. Tests of functional limitations in fibromyalgia syndrome: a reliability study. Arthritis Care Res 1999; 12: 193–199.

16. Nicklin J, Cramp F, Kirwan J, Greenwood R, Urban M, Hewlett S. Measuring fatigue in rheumatoid arthritis: a cross-sectional study to evaluate the Bristol Rheumatoid Arthritis Fatigue Multi-Dimensional questionnaire, visual analog scales, and numerical rating scales. Arthritis Care Res (Hoboken) 2010; 62: 1559–1568.

17. Dures EK, Hewlett SE, Cramp FA, Greenwood R, Nicklin JK, Urban M, et al. Reliability and sensitivity to change of the Bristol Rheumatoid Arthritis Fatigue scales. Rheumatology (Oxford) 2013; 52: 1832–1839.

18. Herrmann C. International experiences with the Hospital Anxiety and Depression Scale--a review of validation data and clinical results. J Psychosom Res 1997; 42: 17–41.

19. Fries JF, Spitz P, Kraines RG, Holman HR. Measurement of patient outcome in arthritis. Arthritis Rheum 1980; 23: 137–145.

20. Vittinghoff E, Glidden DV, Shiboski SC, McCulloch CE. Regression methods in biostatistics: linear, logistic, survival, and repeated measures. Second ed. Boston: Springer Science & Business Media; 2012.

21. Bursac Z, Gauss CH, Williams DK, Hosmer DW. Purposeful selection of variables in logistic regression. Source Code Biol Med 2008; 3: 17.

22. Nakagawa S, Schielzeth H, O’Hara RB. A general and simple method for obtaining R2 from generalized linear mixed effects models. Methods Ecol Evolut 2013; 4: 133–142.

23. Brekke M, Hjortdahl P, Kvien TK. Self-efficacy and health status in rheumatoid arthritis: a two-year longitudinal observational study. Rheumatology (Oxford) 2001; 40: 387–392.

24. Mancuso CA, Rincon M, Sayles W, Paget SA. Psychosocial variables and fatigue: a longitudinal study comparing individuals with rheumatoid arthritis and healthy controls. J Rheumatol 2006; 33: 1496–1502.

25. Treharne GJ, Lyons AC, Hale ED, Goodchild CE, Booth DA, Kitas GD. Predictors of fatigue over 1 year among people with rheumatoid arthritis. Psychol Health Med 2008; 13: 494–504.

26. Meeus M, Vervisch S, De Clerck LS, Moorkens G, Hans G, Nijs J. Central sensitization in patients with rheumatoid arthritis: a systematic literature review. Semin Arthritis Rheum 2012; 41: 556–567.

27. Sariyildiz MA, Batmaz I, Bozkurt M, Bez Y, Cetincakmak MG, Yazmalar L, et al. Sleep quality in rheumatoid arthritis: relationship between the disease severity, depression, functional status and the quality of life. J Clin Med Res 2014; 6: 44–52.

28. van Dartel SA, Repping-Wuts JW, van Hoogmoed D, Bleijenberg G, van Riel PL, Fransen J. Association between fatigue and pain in rheumatoid arthritis: does pain precede fatigue or does fatigue precede pain? Arthritis Care Res (Hoboken) 2013; 65: 862–869.

29. Laursen BS, Bajaj P, Olesen AS, Delmar C, Arendt-Nielsen L. Health related quality of life and quantitative pain measurement in females with chronic non-malignant pain. Eur J Pain 2005; 9: 267–275.

30. Bagnato G, De Andres I, Sorbara S, Verduci E, Corallo G, Ferrera A, et al. Pain threshold and intensity in rheumatic patients: correlations with the Hamilton Depression Rating scale. Clin Rheumatol 2015; 34: 555–561.

31. Rongen-van Dartel SA, Repping-Wuts H, van Hoogmoed D, Knoop H, Bleijenberg G, van Riel PL, et al. Relationship between objectively assessed physical activity and fatigue in patients with rheumatoid arthritis: inverse correlation of activity and fatigue. Arthritis Care Res (Hoboken) 2014; 66: 852–860.

32. Henchoz Y, Bastardot F, Guessous I, Theler JM, Dudler J, Vollenweider P, et al. Physical activity and energy expenditure in rheumatoid arthritis patients and matched controls. Rheumatology (Oxford) 2012; 51: 1500–1507.

33. Neuberger GB, Aaronson LS, Gajewski B, Embretson SE, Cagle PE, Loudon JK, et al. Predictors of exercise and effects of exercise on symptoms, function, aerobic fitness, and disease outcomes of rheumatoid arthritis. Arthritis Rheum 2007; 57: 943–952.

34. Conn VS, Hafdahl AR, Minor MA, Nielsen PJ. Physical activity interventions among adults with arthritis: meta-analysis of outcomes. Semin Arthritis Rheum 2008; 37: 307–316.

35. Durcan L, Wilson F, Cunnane G. The effect of exercise on sleep and fatigue in rheumatoid arthritis: a randomized controlled study. J Rheumatol 2014; 41: 1966–1973.

36. Sokka T, Hakkinen A, Kautiainen H, Maillefert JF, Toloza S, Mork Hansen T, et al. Physical inactivity in patients with rheumatoid arthritis: data from twenty-one countries in a cross-sectional, international study. Arthritis Rheum 2008; 59: 42–50.

37. Stavropoulos-Kalinoglou A, Metsios GS, Koutedakis Y, Kitas GD. Obesity in rheumatoid arthritis. Rheumatology (Oxford) 2011; 50: 450–462.

38. Summers GD, Deighton CM, Rennie MJ, Booth AH. Rheumatoid cachexia: a clinical perspective. Rheumatology (Oxford) 2008; 47: 1124–1131.

39. Hawley DJ, Wolfe F, Lue FA, Moldofsky H. Seasonal symptom severity in patients with rheumatic diseases: a study of 1,424 patients. J Rheumatol 2001; 28: 1900–1909.

40. Davis MC, Zautra AJ, Younger J, Motivala SJ, Attrep J, Irwin MR. Chronic stress and regulation of cellular markers of inflammation in rheumatoid arthritis: implications for fatigue. Brain Behav Immun 2008; 22: 24–32. 

41. Cella D, Lai JS, Stone A. Self-reported fatigue: one dimension or more? Lessons from the Functional Assessment of Chronic Illness Therapy–Fatigue (FACIT-F) questionnaire. Support Care Cancer 2011; 19: 1441–1450. 

42. Feldthusen C, Dean E, Forsblad-d’Elia H, Mannerkorpi K. Effects of person-centered physical therapy on fatigue-related variables in persons with rheumatoid arthritis: a randomized controlled trial. Arch Phys Med Rehabil 2016; 97: 26–36.

Supplementary content
Table SI

Comments

Do you want to comment on this paper? The comments will show up here and if appropriate the comments will also separately be forwarded to the authors. You need to login/create an account to comment on articles. Click here to login/create an account.