Objective Identify factors associated with presence and extension of spinal and sacroiliac joints (SIJ)–MRI lesions suggestive of axial spondyloarthritis (axSpA) in a population-based cohort (Study of Health in Pomerania) aged <45 years.
Methods Spinal (sagittal T1/T2) and SIJ (semicoronal STIR sequences) MRIs were evaluated by two trained blinded readers. The presence (yes/no) and extension (Berlin MRI Score) of bone marrow oedema (BME) were captured. Degenerative spinal lesions were excluded and discrepancies resolved by consensus. Cross-sectional associations between clinical factors and presence/extension of BME were analysed by logistic/negative binomial regression. Record linkage of claims data was applied to identify participants with axSpA.
Results MRIs of 793 volunteers were evaluated. The presence of SIJ–BME (odds ratio) was strongly associated delivery during the last year (4.47, 1.49–13.41). For SIJ–BME extension, associations (incidence rate ratios, 95% CI) were found for delivery ((during last year) 4.52, 1.48–13.84), human leucocyte antigen (HLA)-B27+ (2.32, 1.30–4.14), body mass index (25–30 vs <25 kg/m²; 1.86 (1.19–2.89)) and back pain ((last 3 months) 1.55, 1.04–2.31), while for spinal BME, associations were found for age per decade (1.46, 1.13–1.90) and physically demanding work (1.46, 1.06–2.00). Record linkage was available for 694 (87.5%) participants and 9/694 (1.3%) had a record of axSpA (ICD M45.09).
Conclusion These population-based data support the hypothesis of mechanic strain contributing to BME in the general population aged <45 years and the role of HLA-B27+ as a severity rather than a susceptibility factor for SIJ–BME.
- magnetic resonance imaging
Statistics from Altmetric.com
What is already known about this subject?
There is a relatively high frequency of inflammatory and fatty spinal/inflammatory sacroiliac joints (SIJ) and spinal MRI lesions suggestive of has been suspected.
The reasons for these false-positive signals are still insufficiently investigated.
What does this study add?
Human leucocyte antigen (HLA)-B27+, delivery during the last year in female adults and presence of back pain in the last 3 months seem to be the most important predictors for the extent of bone marrow oedema (BME) in the SIJ, while age and physically demanding work seem to be the most relevant predictors for the extent of BME in the spine.
These data also support the role of HLA-B27 as a severity rather than a susceptibility factor for the occurrence of BME in the SIJ.
How might this impact on clinical practice or future developments?
There is a different association between the occurrence and extension of BME in the SIJ in persons with high body mass index.
The associations presented here are relevant for the interpretation of spinal MRIs of young people with back pain in blue collar jobs and their possible referral to rheumatologists for exclusion or confirmation of spondyloarthritis.
Axial spondyloarthritis (axSpA) is a chronic rheumatic disease that is characterised by inflammatory back pain and several other musculoskeletal and extramusculoskeletal disease manifestations and comorbidities.1 The classification criteria published in 2009 by the Assessment in Spondyloarthritis International Society (ASAS) have set the scene for a differentiation between the classical ankylosing spondylitis (AS), or radiographic SpA, and the non-radiographic form of SpA, based on the presence or absence of definite radiographic changes in the sacroiliac joints (SIJ).2 Conventional radiography and MRI as imaging methods and human leucocyte antigen (HLA)-B27 as a laboratory test play an important role for diagnosis and classification of axSpA.3 Recently, ASAS has published updates on the recommendations of how to interprete MRIs of the SIJ and the spine in patients with axSpA.4 5 Yet, recent studies have shown that not only axSpA can be related to such changes in the axial skeleton and that differential diagnosis may sometimes be difficult.6 Among others, we have recently demonstrated that there is a large proportion of false-positive MRIs suggestive of inflammatory activity as seen in axSpA, as well as frequent other MRI changes such as fat lesions, in the axial skeleton in the normal population.7 8 Furthermore, these changes may also occur in situations with physical demanding activities9 or after pregnancy.10
The pathogenesis of axSpA is not fully understood to date but it seems to be largely genetically determined.11 HLA-B27 is itself responsible for about 20% of the total genetic risk11 and HLA-B27+ individuals carry a 10–20-fold increased risk of developing SpA.12 In the last two decades, genetic studies have provided major insights into this topic, by identifying susceptibility alleles including >100 established loci which contribute roughly to 10% of the heritability of the disease, over and above the major effect of HLA-B27.11 Although the association of AS with HLA-B27 was already reported in 1973,12 it took more than 30 years to establish its role for the classification,2 diagnosis13 and the referral14 of patients with axSpA. The clinical relevance of HLA-B27 has recently also been demonstrated in patients with psoriatic arthritis15 and chronic inflammatory bowel disease.16
Using data from participants of a large population-based study <45 years,17 18 we conducted a cross-sectional analysis to identify factors associated with MRI changes that are known to occur frequently in patients with axSpA in the SIJ and the spine.
The population-based project ‘Study of Health in Pomerania’ (SHIP) comprises the two separate cohorts SHIP and SHIP-TREND that are sampled in a north-eastern region of Germany and followed up every 4–6 years. SHIP is part of the Community Medicine Research net of the University Medicine of Greifswald, Germany, and conducts an extensive clinical examination programme in each follow-up of the cohorts17
Within the second follow-up of SHIP-2 and the baseline examination of SHIP-TREND-0, whole-body MRI was performed. For the present study, MRIs of the 793 volunteers being <45 years at the day of the MRI examination and who had complete MRI sets (both spine and SIJ) were included. The selection of volunteers from the SHIP project has been described in more detail elsewhere.17
MRI and reading of images
MRI was performed by the Department of Diagnostic Radiology and Neuroradiology in at one study site using one MRI device (Magnetom Avanto, Siemens Medical Systems, Erlangen, Germany) and under the same standardised protocol.18 For the entire spine T1 (TR 6.760 ms/TE 120 ms, flip angle 150°, slice thickness 4 mm, scan time 2:42 min) and T2 (TR 37.600 ms/TE 1.060 ms, flip angle 180°, slice thickness 4 mm, scan time 2:04 min) MRI sequences were available in a sagittal and for the SIJ respective STIR (TR 48.910 ms/TE 670 ms, flip angle 180°, slice thickness 5 mm) sequences in a semicoronal orientation were available. All images were blinded for additional participant information such as age and gender. Two readers, who first completed a training session of reference images including patients with axSpA, evaluated the MRI independently in a paired fashion to assess bone marrow oedema (BME) (defined as hyperintense signal on T2-weighted and hypointense signal on T1-weighted images in the spine or as hyperintense signal on STIR sequences in the SIJs). The ASAS definitions were used for defining lesions as ‘positive’.4 19 Spinal lesions with pathologic changes involving the vertebral endplate or being accompanied by abnormalities of the intervertebral disc (obvious dehydration, protrusion or prolapse) were considered as degenerative and were not counted. In case of disagreement for a lesion being present or not present, both readers adjudicated together in order to reach consensus.
In addition to a the binary approach of lesions being present or absent, quantification of the extension of BME lesions was performed based on the definitions of the Berlin SIJ and spine MRI Score.20 Briefly, this score captures BME in a scale of 0–3. Both for the SIJ and the spine, a score of 0 means no BME, while a score of 1, 2 and 3 mean 0% to ≤33%, >33% to ≤66% and >66% of an SIJ quadrant or and 0% to ≤25%, >25% to ≤50% and >50% of a discovertebral unit, respectively. The final scores used for analysis were calculated based on consent of both readers in cases of disagreement.
Collection of clinical information
Clinical information was collected for age (in years), sex, smoking habit (current smoker, previous smoker, no smoker), mean spinal back pain level in the last 3 months prior to the MRI examination (on a numerical rating scale (NRS) rated 0–10), high-sensitivity C reactive protein (hsCRP, in mg/dL), HLA-B27 status (positive or negative), body mass index (BMI) (categorised according to the definitions of WHO, WHO to normal or underweight, overweight and obese), blue-collar or white-collar job, and whether the (female) patients had given birth within the last 12 months prior to the MRI examination.
Linkage of claims data
Claims data were available from the regional Association of Statutory Health Insurance including ICD-10 codes between 2002 and 2018. These data were used to obtain all ICD-10 codes of AS M45.xx including the quarterly period after the SHIP examination has been conducted.
As the claims data and those of the SHIP study have no common key for linkage, we applied record linkage as proposed in Vatsalan et al 21 based on: surname, name, date of birth and sex of participants with linkage consent. The names were normalised to upper-case letters and by removal of special characters, the indexing of candidate pairs was blocked via birth date. After comparison of candidate pairs, respective claims data were available for 694 of 769 (90.2%) participants who consented (769 of 793, 97.0%) in record linkage.
Descriptive measures (mean, median) are shown with SD, minimum and maximum. Frequencies and percentages are provided for categorical data (N). For all variables, the amount of missing values is reported.
The outcomes of affected SIJ quadrants or spinal segments are count data and were, therefore, modelled using negative binomial regression. The decision for a negative binomial regression model versus a Poisson model was based on likelihood ratio test in all outcome categories (SIJ, spine) and in the complete data. Effect estimates of the regression models were exponentiated resulting in incidence rate ratios (IRRs).22 Associations for the presence of BME was modelled using multiple logistic regression (BME yes/no).
Missing values were accounted for by multiple imputations using the R-package mice.23 Given the low number of missing values in each variable, m=10 imputations were considered sufficient. The imputation model comprised all variables of the analysis model. Imputations for categorical variables were conducted using logistic or generalised logistic regression; imputations of numerical variables were conducted using predictive means. All results from the regression models are pooled over all imputations according to Rubins’ rule.24 The variance increase due to missing data was highest for HLA-B27 with 3.3%. We found no relevant difference between observations with complete data on those with incomplete data (data not shown).
In a sensitivity analysis we added information on known ICD-10 codes (M45.xx) to the regression models to examine changes and the robustness of results.
The interpretation of results is conducted in compliance with recommendations to avoid dichotomisation into statistically significant and insignificant results.25 The use of p values is restricted to model diagnostics; all effects are described using effect estimates and 95% CIs.
A total of 793 MRIs of volunteers were evaluated by the two readers, with mean age 37.3±6.3 years, 49.4% male, 67 (8.4%) HLA-B27+, 53 (6.7%) with hsCRP >0.5 mg/dL, 451 (56.9%) with back pain in the last 3 months, (228 (28.8%) with back pain ≥4/10 on an NRS) were available. A physically demanding job was reported by 283 individuals (35.7%), 511 (64.3%) reported to work at a desktop. More than half (n=436) of the participants (55%) had a BMI >25 kg/m2, and 247 were current smokers (31.1%). Delivery in the last year before the MRI examination was reported by 16 females (5%). Average physical activity >1 hour/day was reported by 653 subjects (82.4%). All characteristics including the numbers of missing values are shown in the table 1.
For the extension of BME on SIJ-MRIs, the largest IRR (95% CI) were found for delivery in the last year: 4.52 (1.48 to 13.84), HLA-B27+ : 2.32 (1.30 to 4.14) and BMI 25–30 versus <25 kg/m2: 1.86 (1.19 to 2.89). For the extension of BME on spinal MRIs, IRRs were overall lower in size with the largest effects found for age per decade increase: 1.46 (1.13 to 1.90) and physically demanding work: 1.46 (1.06 to 2.00), see figure 1.
Regarding the presence of positive lesions on MRI only for the occurrence of delivery during the last year a strong association with SIJ–BME was found in females, with an OR of 4.47 (95% CI: 1.49 to 13.41) (figure 1).
Finally, in the comparison between SIJ and spinal MRIs, participants with back pain in the last 3 months (62.5% vs 56.9%) had more often SIJ BME than spinal BME, while spinal BME was more frequent than SIJ BME in participants working at a desktop (61.5% vs 54.4%), while smokers (66.9% vs 63.8%) and participants with back pain in the last 3 months (62.5% vs 56.9%) had more often SIJ BME than spinal BME, respectively, (data not shown).
In a subset of participants of the SHIP cohorts, information on ICD codes of all participants who consented to record linkage was available from the claims data: in 9/694 (1.3%) the ICD-10 code M45.09 (AS without specified location) was documented prior to the SHIP examination. However, an affirmation of the ‘real’ diagnoses of AS or axSpA by an expert was not obtained, since this was not part of the study protocol.
Results from sensitivity analyses did not reveal relevant differences to the results presented above. The largest change was found regarding the IRR for HLA-B27+ in SIJ–BME which decreased from 2.32 (1.30; 4.14) to 2.21 (1.24; 3.95).
In this population-based study, HLA-B27+, delivery during the last year in female participants and presence of back pain in the last 3 months were the most important predictors for the extent of BME in the SIJ, while age and physically demanding work were the most relevant predictors for the extent of BME in the spine.
These data support the hypothesis that mechanic strain contributes to BME found in the general population in subjects not diagnosed with SpA. In addition, HLA-B27+ was associated with the extent of BME but only contributes to a minor extent as a susceptibility factor for BME in the SIJ, while no association with spinal MRI changes was detected. Although the link between mechanically induced inflammation, genes and the immune system still needs to be elucidated, our data suggest, for the first time, that the recently supported26–28 hypothesis coming from studies with patients also applies on a general population level.
Our finding of a clear association of BME in the SIJ after delivery within a relatively short period of time prior to the performance of the MRI is based on a small number of events (n=16). However, the result is clinically plausible, looking at the reported prevalence of MRI29 and CT30 findings of osteitis condensans ilii (OCI) in young patients with symptoms of low back pain. In a recent Belgian study,10 27 (77%) out of 35 women developed sacroiliac BME immediately postpartum, with 60% fulfilling the ASAS definition of a positive MRI in axSpA. Importantly, after 6 months, BME was still present in 46% of the participants. Whether and how OCI, a condition with a rather obvious mechanic pathogenesis, is associated with axSpA is unclear at present. In a recent study on OCI from Berlin, sacroiliac BME occurred frequently and equally often in both groups but only 7.4% of patients with OCI developed erosions in the SIJ even though about 35% of women with OCI were HLA-B27+.31
The association of back pain with the extension of BME in the SIJ but not the spine, even though the latter was more frequent, confirms its pathologic relevance related to the development of SpA, including the size of such lesions.4 Nevertheless, it needs to be taken into account that the participants in this cohort study were not patients complaining about chronic or inflammatory back pain but were individuals from the general population. Thus, it may be expected that about 6% of the participants in this study had inflammatory back pain,32 and based on a population prevalence of axSpA of 0.5%–1%, about four to eight patients may also have had axSpA.33 Since it was not planned in this cohort study to directly diagnose patients, we cannot provide data on classification. Nevertheless, results obtained from record linkage of claims data showed that there were no established cases of axSpA in this study population, even though nine participants had been given the ICD code M45.09 somewhere in outpatient care at least once. This result is in line with the clinical experience that axSpA cases will be identified late or not at all.34 These nine subjects were more likely to be HLA-B27+ and to have MRI changes in the SIJ (data not shown); this is consistent with this interpretation. It is also worth mentioning that the successful use of non-steroidal anti-rheumatic drugs (NSAIDs), which is not only part of the classification criteria for axSpA but sometimes also used as a diagnostic tool, was only present in a few patients and not associated with the presence of BME in this study (data not shown).
Furthermore, the prevalence of back pain, known to increase with age35 during the first life decades,36 was confirmed for subjects below the age of 45 in our study. The result that participants working at a desktop had more often BME in the spine than in the SIJ may be methodologically explained since the differentiation of the severity of the physical workload (eg, blue or white collar) was not very sharp in this study. However, the finding that the extent of spinal BME was significantly associated with physically demanding work makes sense from the clinical point of view and is consistent with previous work.37 Taken together, these results confirm that BME does occur in the younger general population with no established diagnosis of axSpA. Since physically demanding work is known to be associated with chronic back pain,38 these associations are relevant for the interpretation of spinal MRIs of young people with back pain in blue collar jobs. Other assessments of physical activity showed smaller associations with the occurrence or extent of BME in our study. However, we did not assess the intensity of physical activity in more detail. In another study, physical inactivity over 14 years was shown to be associated with disc degeneration in the thoracic and the lumbar spine.39 Similarly, physically demanding work may also have an effect of patients with axSpA. This is not necessarily related to inflammation.40
In addition, while the association of back pain with BMI and fat mass is well established,41 our data suggest that there is also a different association between the occurrence and extension of BME in the SIJ in persons with high BMI (OR vs IRR, figure 1). This also seems to point in the direction of mechanical strain as a cause of more severe BME and is likely to stimulate research in this area.
Finally, only no or small associations between the occurrence of BME with any other of the investigated factors, such as smoking, hsCRP or female sex were found in our cohort. All of those factors are special interest in the overall topic of SpA-related MRI findings, since they represent established risk factors for low back pain42 and also have been found to have an impact on the radiographic course of axSpA.43
In conclusion, in this cross-sectional analysis of MRIs of the SIJ and spine of participants of a population-based study, we identified predictors relevant to the presence of BME suggestive of axSpA. Our findings support the hypothesis of a mechanic origin of BME in the general population aged <45 years. On the other hand, HLA-B27 was identified to be a severity rather than a susceptibility factor especially for the presence of BME in the SIJ. Longitudinal analyses may be able to demonstrate causal relationships.
Study of Health in Pomerania is part of the Community Medicine Research Net of the University Medicine Greifswald, which is supported by the Federal State of Mecklenburg-West Pomerania, Germany, as well as by a joint grant from Siemens Healthineers, Erlangen.
COS and JB are joint senior authors.
Handling editor Josef S Smolen
XB and AR contributed equally.
Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.
Competing interests None declared.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Patient consent for publication Not required.
Ethics approval The study was approved by the local ethics committee of the University Medicine Greifswald.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data are available upon reasonable request. All data relevant to the study are included in the article or uploaded as supplementary information. All analyses are described in the manuscript. For additional information, please contact the corresponding author.
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.