A total of 828 axSpA patients were included in the analysis, comprising 568 (68.6%) males and 260 (31.4%) females. Table 1 summarizes the demographic and clinical data of the sample, comparing men and women. Regarding demographic data, women had a later median age at symptom onset (29 years, IQR 19-34 vs. 27 years, IQR 21-38; p=0.001) and age at diagnosis (37 years, IQR 30–47 vs. 35 years, IQR 27–44; p=0.002), but shorter disease duration (16 years, IQR 9–26 vs. 18 years, IQR 10–28; p=0.041) when compared to men. Despite these differences, diagnostic delay was similar between sexes (median 5 years, IQR 1–13 for both; p=0.975). Although HLA-B27 positivity was more frequent among males (77.6% vs. 69.2%; p=0.009; OR=1.33, 95%CI 1.07–1.64), family history of SpA was more frequently reported by females (26.2% vs. 17.2%; p=0.002; OR=1.71, 95%CI 1.21–2.40).
While axial symptoms presented similar frequencies in both sexes, active peripheral involvement was more common in females (25.3% vs. 19.1%; p=0.040; OR=1.43, 95%CI 1.01–2.02). Emphasizing this finding, upper limb arthritis was significantly more reported by females (21.6% vs. 14.8%; p=0.010; OR=1.58, 95% CI 1.11–2.25). Concerning extra-musculoskeletal manifestations, no statistically significant differences were found between the sexes. Regarding comorbidities, psychological distress was more prevalent among females (43.8% vs. 20.9%; p<0.001; OR=2.95 95%CI 2.16-4.01) while chronic kidney disease was more frequently reported by males (3.0% vs. 0.7%; p=0.028; OR=0.22, 95%CI 0.05–0.96) (table 1).
The prevalence of arterial hypertension, diabetes and obesity was similar between sexes. Psychological distress was more prevalent among females (43.8% vs. 20.9%; p<0.001; OR=2.95 95%CI 2.16-4.01) and chronic kidney disease was more frequently reported among males (3.0% vs. 0.7%; p=0.028; OR=0.22, 95%CI 0.05–0.96).
The overall proportion of physically active individuals did not differ significantly between sexes. However, a greater proportion of men reported engaging in physical exercise with higher frequency, specifically 4–7 times per week (45.7% vs. 29.9%). Regarding employment status, men were significantly more actively employed (48.9% vs. 32.1%; p<0.001; OR=0.27, 95% CI 0.17–0.41), as well as were more frequently retired due to disability (63.6% vs. 28.5%; p<0.001; OR=2.06, 95% CI 1.63–2.60) among those who were not currently working (table 1).
There were no significant differences in prior treatment patterns. Regarding current treatments, sulfasalazine was more commonly used by men (16.0% vs. 10.7%; p=0.029 ; OR=0.62, 95%CI 0.40–0.95), and oral corticosteroid use was more frequent among women (4.5% vs. 2.1%; p=0.042; OR=2.18, 95%CI 1.01–4.70). Table 1 summarizes the demographic and clinical data of the sample, comparing men and women (table 1).
Sex-stratified analysis of continuous clinical variables is presented in Table 2. Women had a later median age at symptom onset (29 years, IQR 21–38) compared to men (27 years, IQR 19–34) (p=0.001) and were also diagnosed at an older age (37 years, IQR 30–47 vs. 35 years, IQR 27–44; p=0.002). Despite these differences, diagnostic delay was similar between sexes (median 5 years, IQR 1–13 for both; p=0.975). There were no significant differences in CRP levels, enthesitis count, mSASSS, ASAS NSAID intake score, swollen joint count (SJC), or tender joint count (TJC). However, VAS PGA (5.0, IQR 3.5–7.0 vs. 5.0, IQR 2.0–7.0; p<0.001), pain (VAS 5.0, IQR 3.0–7.0 vs. 4.0, IQR 2.0–7.0; p<0.001), and VAS MGA (4.0, IQR 2.0–6.0 vs. 3.0, IQR 1.0–5.0; p=0.006) were all higher in women.
Sex-related differences in key clinical outcomes are illustrated in Figure 1. Female patients had significantly higher disease activity scores, both by BASDAI (median 4.5 vs. 3.2; p<0.0001) and by ASDAS-CRP (median 2.2 vs. 1.9; p<0.0001). Functional impairment, assessed by BASFI (median 5.0 vs. 4.0; p=0.0012), and quality of life, measured by ASQoL (median 9.0 vs. 7.0; p<0.0001), were also worse in women. In contrast, spinal mobility assessed by BASMI was more restricted in men (median 4.0 vs. 3.5; p=0.0023).
Table 1. Comparative analysis between male (n=568) and female (n=260) patients with axSpA across demographic, clinical and disease-related variables
|
|
Male (n=568)
|
Female (n=260)
|
|
|
|
|
n
|
%
|
n
|
%
|
p
|
OR (95%CI)
|
|
Skin color
|
|
|
|
|
|
|
|
White
|
279
|
41.9
|
121
|
41.6
|
0.928
|
0.98 (0.74-1.30)
|
|
Non White
|
387
|
58.1
|
170
|
58.4
|
|
Active at work
|
301
|
48.9
|
137
|
49.8
|
0.809
|
0.27 (0.17-0.41)
|
|
Retired due to illness
|
196
|
63.6
|
43
|
32.1
|
<0.001
|
2.06 (1.63-2.60)
|
|
Comorbidities
|
|
|
|
|
|
|
|
Hypertension
|
238
|
35.7
|
103
|
35.4
|
0.932
|
1.00 (0.82-1.23)
|
|
Diabetes
|
95
|
14.2
|
37
|
12.7
|
0.528
|
0.87 (0.58-1.31)
|
|
Obesity
|
61
|
9.1
|
25
|
8.6
|
0.783
|
0.93 (0.57-1.52)
|
|
Chronic renal diseases
|
20
|
3.0
|
2
|
0.7
|
0.028
|
0.22 (0.05-0.96)
|
|
Coronary heart disease
|
9
|
1.3
|
5
|
1.7
|
0.662
|
1.27 (0.42-3.84)
|
|
Psychological distress
|
125
|
20.9
|
121
|
43.8
|
<0.001
|
2.95 (2.16-4.01)
|
|
No comorbidities
|
299
|
44.8
|
102
|
35.1
|
0.005
|
0.66 (0.49-0.88)
|
|
Smoking
|
78
|
12.2
|
19
|
6.6
|
0.010
|
1.64 (1.08-2.49)
|
|
Alcoholism
|
41
|
6.5
|
4
|
1.4
|
<0.001
|
3.63 (1.41-9.30)
|
|
Physically active
|
270
|
43.8
|
118
|
41.8
|
0.577
|
0.92 (0.69-1.22)
|
|
1-3 times/week
|
146
|
54.3
|
82
|
70.1
|
0.004
|
0.50 (0.31-0.80)
|
|
4-7 times/week
|
123
|
45.7
|
35
|
29.9
|
|
|
SpA family history
|
106
|
17.2
|
73
|
26.2
|
0.002
|
1.71 (1.21-2.40)
|
|
HLA-B27 positive
|
441
|
77.6
|
180
|
69.2
|
0.009
|
1.33 (1.07-1.64)
|
|
First symptoms
|
|
|
|
|
|
|
|
Cervical pain
|
236
|
35.4
|
92
|
31.6
|
0.258
|
1.12 (0.91-1.38)
|
|
Hip pain
|
151
|
22.6
|
62
|
21.3
|
0.648
|
0.92 (0.66-1.29)
|
|
Dactylitis
|
29
|
4.3
|
12
|
4.1
|
0.875
|
1.04 (0.64-1.69)
|
|
Buttock pain
|
257
|
38.5
|
120
|
41.2
|
0.430
|
1.12 (0.84-1.48)
|
|
Enthesitis
|
166
|
24.9
|
64
|
22.0
|
0.335
|
0.85 (0.61-1.18)
|
|
Low back pain
|
583
|
87.4
|
247
|
84.9
|
0.291
|
0.80 (0.54-1.19)
|
|
Peripheral arthritis
|
278
|
41.7
|
125
|
43.0
|
0.713
|
1.05 (
|
|
EMM
|
|
|
|
|
|
|
|
Uveitis
|
203
|
30.4
|
98
|
33.7
|
0.320
|
1.16 (0.86-1.57)
|
|
IBD
|
16
|
2.4
|
12
|
4.1
|
0.145
|
1.75 (0.81-3.74)
|
|
Psoriasis
|
9
|
1.3
|
9
|
3.1
|
0.068
|
2.33 (0.91-5.9)
|
|
Active peripheral involvement
|
111
|
19.1
|
68
|
25.3
|
0.040
|
1.43 (1.01-2.02)
|
|
Previous treatment
|
|
|
|
|
|
|
|
NSAID
|
287
|
43.0
|
123
|
42.3
|
0.827
|
0.96 (0.73-1.28)
|
|
Methotrexate
|
168
|
25.2
|
72
|
24.7
|
0.884
|
0.97 (0.71-1.34)
|
|
TNFi
|
236
|
35.4
|
115
|
39.5
|
0.222
|
1.19 (0.89-1.58)
|
|
SSZ
|
210
|
31.5
|
90
|
30.9
|
0.864
|
0.97 (0.72-1.31)
|
|
Oral corticosteroid
|
55
|
8.2
|
28
|
9;6
|
0.486
|
1.18 (0.73-1.91)
|
|
IL17i
|
15
|
2.2
|
10
|
3.4
|
0.289
|
1.54 (0.68-3.48)
|
|
Current treatment
|
|
|
|
|
|
|
|
TNFi
|
418
|
62.7
|
171
|
58.8
|
0.253
|
0.84 (0.64-1.12)
|
|
IL17i
|
54
|
8.1
|
30
|
10.3
|
0.265
|
1.30 (0.81-2.08)
|
|
Methotrexate (oral)
|
82
|
12.3
|
35
|
12.0
|
0.908
|
0.97 (0.64-1.48)
|
|
SSZ
|
107
|
16.0
|
31
|
10.7
|
0.029
|
0.62 (0.40-0.95)
|
|
Oral corticosteroid
|
14
|
2.1
|
13
|
4.5
|
0.042
|
2.18 (1.01-4.70)
|
Categorical variables are shown as absolute frequencies and percentages. P-values indicate the significance of the difference between sexes, calculated using the chi-square test. Odds Ratios (OR) with 95% Confidence Intervals (CI) were calculated to represent the likelihood of occurrence in females compared to males. Statistical significance was considered when p < 0.05.
OR – Odds Ratio; 95% CI – 95% Confidence Interval; SD – Standard Deviation; IBD – Inflammatory Bowel Disease; NSAID: Nonsteroidal Anti-Inflammatory Drug; TNFi: Tumor Necrosis Factor inhibitor; SSZ: sulfasalazine IL-17i: Interleukin-17 inhibitor
Table 2. Sex-based comparison of clinical and functional measures in patients with axSpA
|
Variable
|
Male Median
(P25–P75)
|
Female Median (P25–P75)
|
p-value
|
|
Age
|
49 (38-57)
|
50 (40-57)
|
0.351
|
|
Age at symptoms onset
|
27 (19-34)
|
29 (21-38)
|
0.001
|
|
Age at diagnosis
|
35 (27-44)
|
37 (30-47)
|
0.002
|
|
Diagnosis delay
|
5 (1-13)
|
5 (1-13)
|
0.975
|
|
Disease duration
|
18 (10-28)
|
16 (9-26)
|
0.041
|
|
Enthesitis count
|
4 (2-9)
|
4 (2–11)
|
0.315
|
|
mSASSS
|
14 (5-35)
|
14 (4–23)
|
0.244
|
|
ASAS NSAID intake score
|
23.4 (9.5–53.6)
|
19.0 (7.10–50.0)
|
0.204
|
|
VAS PGA
|
5.0 (2.0–7.0)
|
5.0 (3.5–7.0)
|
<0.001
|
|
VAS pain
|
4.0 (2.0–7.0)
|
5.0 (3.0–7.0)
|
<0.001
|
|
VAS MGA
|
3.0 (1.0–5.0)
|
4.0 (2.0–6.0)
|
0.006
|
|
SJC
|
0.0 (0.0–1.0)
|
0.00 (0.0–1.0)
|
0.760
|
|
TJC
|
3.0 (2.0–6.0)
|
4.00 (2.0-13.5)
|
0.241
|
Data are presented as median (IQR: P25–P75). Comparisons between groups were conducted using the Mann Whitney U test for independent samples (95%CI). VAS: visual analogic scale; PGA: patient global assessment; MGA: medical global assessment; TJC: tender joint count; SJC: swollen joint count ; NSAID score: Nonsteroidal Anti-Inflammatory Drug score.
Sex-stratified bivariate analyses showed that musculoskeletal manifestations, particularly shoulder and hip pain and respective functional limitations, reflecting root-joint involvement, were significantly associated with worse BASDAI, ASDAS-CRP, BASFI, BASMI, and ASQoL scores in both men and women. Among male patients, additional factors consistently associated with worse outcomes across multiple domains included physical inactivity, not currently working, active peripheral arthritis, and comorbidities such as hypertension, diabetes, smoking, and obesity. In female patients, psychological distress was significantly associated with higher BASFI and ASQoL scores, while obesity and reduced physical activity were also linked to worse ASQoL. Detailed results are presented in the supplementary material.
In the multivariable binomial logistic regression model, male sex was significantly associated with HLA-B27 positivity (OR=1.66; 95%CI 1.13-2.42; p=0.009), smoking (OR=2.55; 95%CI 1.34-4.88; p=0.005), and hip limitation (OR=1.85; 95%CI 1.26-2.71; p=0.002). In contrast, female patients showed higher rates of psychological distress (OR=0.38; 95%CI 0.26–0.54; p<0.001) and reported higher disease activity scores as measured by BASDAI (OR=0.82 per unit increase; 95%CI 0.76–0.88; p<0.001) (Figure 2).
Table 3 presents the results of multivariable linear regression models for key outcomes in axSpA, stratified by sex. For each outcome (BASDAI, ASDAS, BASFI, BASMI, and ASQOL), the variables retained in the final models differed between men and women, highlighting distinct patterns of association. Model fit, as expressed by adjusted R², ranged from 0.398 to 0.710 across models, with similar performance in both sexes overall. However, the specific predictors and the strength of their associations varied, suggesting sex-specific pathways influencing disease outcomes in axSpA.
Table 3 : Sex-specific multivariable linear regression models for key Axial Spondyloarthritis outcomes
|
|
Male
|
Female
|
|
|
β
|
CI95%
Lower
|
CI95%
Upper
|
P-value
|
|
β
|
CI95%
Lower
|
CI95%
Upper
|
P-value
|
|
BASDAI
|
|
Shoulder pain
|
0.556
|
0.207
|
0.904
|
0.002
|
Shoulder pain
|
0.096
|
0.016
|
1.024
|
0.043
|
|
|
|
|
|
|
Hip pain
|
0.137
|
0.229
|
1.192
|
0.004
|
|
Active peripheral
|
0.370
|
0.008
|
0.733
|
0.045
|
|
|
|
|
|
|
BASMI
|
-0.241
|
-0.321
|
-0.161
|
<0.001
|
BASMI
|
-0.106
|
-0.276
|
-0.004
|
0.044
|
|
BASFI
|
0.392
|
0.310
|
0.474
|
<0.001
|
BASFI
|
0.483
|
0.307
|
0.564
|
<0.001
|
|
ASQOL
|
0.170
|
0.131
|
0.209
|
<0.001
|
|
|
|
|
|
|
R=0.787; R²=0.619; Adjusted R² =0.608
|
R=0.777;R²=0.604; Adjusted R²=0.595
|
|
|
ASDAS
|
|
Shoulder pain
|
0.235
|
0.032
|
0.438
|
0.023
|
|
|
|
|
|
|
Hip pain
|
0.217
|
0.039
|
0.395
|
0.017
|
Hip pain
|
0.129
|
0.059
|
0.435
|
0.010
|
|
Active enthesitis
|
0.178
|
0.013
|
0.343
|
0.035
|
|
|
|
|
|
|
BASMI
|
-0.081
|
-0.121
|
-0.042
|
<0.001
|
BASMI
|
-0.130
|
-0.115
|
-0.009
|
0.022
|
|
BASFI
|
0.144
|
0.105
|
0.184
|
<0.001
|
BASFI
|
0.473
|
0.105
|
0.205
|
<0.001
|
|
ASQOL
|
0.060
|
0.041
|
0.078
|
<0.001
|
ASQOL
|
0.336
|
0.035
|
0.082
|
<0.001
|
|
R=0.725; R²=0.525; Adjusted R² =0.511
|
R=0.766;R²=0.587; Adjusted R²=0.574
|
|
|
BASFI
|
|
Hip pain
|
0.053
|
-0.002
|
0.639
|
0.052
|
|
|
|
|
|
|
Current working
|
-0.105
|
-0.934
|
-0.294
|
<0.001
|
|
|
|
|
|
|
Physically active
|
-0.051
|
-0.599
|
-0.003
|
0.048
|
|
|
|
|
|
|
ASDAS
|
0.271
|
0.61
|
1.011
|
<0.001
|
ASDAS
|
1.007
|
0.659
|
1.320
|
<0.001
|
|
BASMI
|
0.353
|
0.402
|
0.546
|
<0.001
|
BASMI
|
0.531
|
0.419
|
0.643
|
<0.001
|
|
ASQOL
|
0.399
|
0.181
|
0.257
|
<0.001
|
ASQOL
|
0.214
|
0.159
|
0.269
|
<0.001
|
|
R=0.45; R²=0.713; Adjusted R² =0.710
|
R=0.844;R²=0.712; Adjusted R²=0.708
|
|
|
BASMI
|
|
Hypertension
|
0.166
|
0.429
|
1.06
|
<0.001
|
|
|
|
|
|
|
Smoking
|
0.159
|
0.607
|
1.528
|
<0.001
|
|
|
|
|
|
|
Radiographic SpA
|
0.120
|
0.467
|
1.739
|
<0.001
|
|
|
|
|
|
|
|
|
|
|
|
Current working
|
-1.233
|
-1.64
|
-0.825
|
<0.001
|
|
|
|
|
|
|
Shoulder pain
|
-0.954
|
-1.436
|
-0.472
|
<0.001
|
|
|
|
|
|
|
Shoulder limitation
|
0.712
|
0.095
|
1.33
|
0.024
|
|
Hip limitation
|
0.141
|
0.306
|
0.94
|
<0.001
|
|
|
|
|
|
|
ASDAS
|
-0.199
|
-0.656
|
-0.234
|
<0.001
|
|
|
|
|
|
|
BASFI
|
0.685
|
0.43
|
0.588
|
<0.001
|
BASFI
|
0.319
|
0.244
|
0.394
|
<0.001
|
|
ASQOL
|
-0.185
|
-0.118
|
-0.033
|
<0.001
|
|
|
|
|
|
|
R=0.669; R²=0.447 Adjusted R² =0.439
|
R=0.639; R²=0.408; Adjusted R² =0.398
|
|
|
ASQOL
|
|
Psychological distress
|
0.138
|
1.043
|
2.556
|
<0.001
|
|
|
|
|
|
|
Current working
|
-0.070
|
-1.439
|
-0.055
|
0.035
|
|
|
|
|
|
|
ASDAS
|
0.299
|
1.210
|
2.062
|
<0.001
|
ASDAS
|
2.288
|
1.394
|
3.182
|
<0.001
|
|
BASFI
|
0.569
|
0.863
|
1.207
|
<0.001
|
BASFI
|
0.830
|
0.507
|
1.154
|
<0.001
|
|
BASMI
|
-0.128
|
-0.492
|
-0.133
|
<0.001
|
|
|
|
|
|
|
R=0.780; R²=0.609; Adjusted R² =0.605
|
R = 0.815. R² = 0.664. and adjusted R² = 0.619
|
|
Results of multivariable linear regression analyses for five outcome variables (BASDAI, ASDAS, BASFI, BASMI, and ASQOL) in male and female patients with axial spondyloarthritis. For each predictor retained in the final model, the standardized regression coefficient (β), 95% confidence interval (CI), and p-value are shown. Only statistically significant predictors (p < 0.05) are included. Coefficients reflect the unique contribution of each variable to the outcome after adjustment for other predictors in the model. Model fit is expressed using R, R², and adjusted R².