Establishment of age- and sex-adjusted reference data for hand bone mass and investigation of hand bone loss in patients with rheumatoid arthritis treated in clinical practice: an observational study from the DANBIO registry and the Copenhagen Osteoarthritis Study

Background Rheumatoid arthritis is characterised by progressive joint destruction and loss of periarticular bone mass. Hand bone loss (HBL) has therefore been proposed as an outcome measure for treatment efficacy. A definition of increased HBL adjusted for age- and sex-related bone loss is lacking. In this study, we aimed to: 1) establish reference values for normal hand bone mass (bone mineral density measured by digital x-ray radiogrammetry (DXR-BMD)); and 2) examine whether HBL is normalised in rheumatoid arthritis patients during treatment with tumour necrosis factor alpha inhibitors (TNFI). Methods DXR-BMD was measured from hand x-rays in a reference cohort (1485 men/2541 women) without arthritis randomly selected from an urban Danish population. Sex- and age-related HBL/year was estimated. DXR-BMD was measured in rheumatoid arthritis patients (n = 350: at start of TNFI, and ~2 years after TNFI start), of which 135 patients had three x-rays (~2 years prior to TNFI, at start of TNFI, and ~2 years after TNFI start). Individual HBL/year prior to and during TNFI was calculated and compared to reference values. Results Estimated HBL/year varied strongly with age and sex. Compared to the reference values, 75 % of 135 patients had increased HBL prior to TNFI treatment and 59 % had increased HBL during TNFI treatment (p = 0.17, Chi-squared). In 38 % (38/101) of patients with increased HBL, HBL was normalised during TNFI treatment, whereas 47 % (16/34) of patients with normal HBL prior to TNFI had increased HBL during TNFI treatment. In the 350 patients, increased HBL during TNFI was associated with time-averaged 28-joint disease activity score (odds ratio 1.69 (95 % Confidence Interval 1.34-2.15)/unit increase, p < 0.001), and patients in time-averaged remission had lower HBL than patients without remission (0.0032 vs. 0.0058 g/cm2/year; p < 0.001, Mann-Whitney). Conclusions We established age- and sex-specific reference values for DXR-BMD in a large cohort without arthritis. HBL was increased in the majority of rheumatoid arthritis patients initiating TNFI in clinical practice, and only normalised in a minority during TNFI. Electronic supplementary material The online version of this article (doi:10.1186/s13075-016-0952-y) contains supplementary material, which is available to authorized users.


(Continued from previous page)
Results: Estimated HBL/year varied strongly with age and sex. Compared to the reference values, 75 % of 135 patients had increased HBL prior to TNFI treatment and 59 % had increased HBL during TNFI treatment (p = 0.17, Chi-squared). In 38 % (38/101) of patients with increased HBL, HBL was normalised during TNFI treatment, whereas 47 % (16/34) of patients with normal HBL prior to TNFI had increased HBL during TNFI treatment. In the 350 patients, increased HBL during TNFI was associated with time-averaged 28-joint disease activity score (odds ratio 1.69 (95 % Confidence Interval 1. 34-2.15)/unit increase, p < 0.001), and patients in time-averaged remission had lower HBL than patients without remission (0.0032 vs. 0.0058 g/cm 2 /year; p < 0.001, Mann-Whitney). Conclusions: We established age-and sex-specific reference values for DXR-BMD in a large cohort without arthritis. HBL was increased in the majority of rheumatoid arthritis patients initiating TNFI in clinical practice, and only normalised in a minority during TNFI.
Keywords: Rheumatoid arthritis, Bone mineral density, Anti-TNF Background Bone involvement in rheumatoid arthritis (RA) has three hallmarks: localised destruction of joint bone, periarticular bone loss and generalised bone loss. The two former features are visible on x-rays as erosions and periarticular osteoporosis and are suggested to be caused by the same pathophysiologic pathway [1,2]. Increased hand bone loss (HBL) measured by digital x-ray radiogrammetry (DXR) may precede the development of erosions, takes place throughout the disease course, and occurs in patients with no progression of radiographic joint damage [3][4][5]. Thus, HBL has been proposed as a more sensitive outcome measure for monitoring disease progression and treatment effect than the traditional evaluation of radiographs [5].
A common definition of increased HBL has not yet been established, which limits the clinical value of HBL and comparison between studies. Attempts to define increased HBL include HBL above the median in the study population, and the smallest detectable difference of HBL or cutoff values provided by the manufacturer based on the distribution of HBL in two early RA cohorts [6][7][8][9][10]. However, the normal HBL in healthy persons, shown to vary considerably with sex and age, has so far not been included in the definition of increased HBL [11,12].
Impact of tumour necrosis factor alpha inhibitors (TNFI) on HBL in RA patients has been investigated in the PREMIER and BEST randomised controlled trials (RCTs) in early RA, where patients treated with TNFI in combination with methotrexate (MTX) had a significantly lower HBL compared with MTX monotherapy [5,13]. No studies have analysed whether HBL is normalised during TNFI treatment.
The primary aims of the present study were: 1) to establish a large reference material for hand bone mass measured by DXR in order to estimate sex-and age-related HBL/year; and 2) to determine whether HBL in RA patients treated in clinical practice normalises during treatment with TNFI when switched from conventional synthetic disease-modifying anti-rheumatic drugs (csDMARDs) due to an unsatisfactory clinical response.
Secondarily, in RA patients treated in clinical practice, we aimed to: 3) investigate the association between HBL and radiographic progression during treatment with csDMARDs and TNFI; 4) investigate the association between HBL and inflammatory activity during TNFI; and 5) identify predictors of increased HBL during TNFI.

Reference cohort
The Copenhagen City Heart Study is a health survey of an adult urban Danish population selected by a random social security number algorithm. From 1991 to 1994, 10,135 participants were included and from this main cohort 1533 men and 2618 women without inflammatory joint disease had x-rays of the pelvis, knees, lumbar spine, hands and wrists performed once (Copenhagen Osteoarthritis Study) [14]. Age, sex and hand x-ray films were used in the present study.

Patients with RA
In the Danish DANBIO registry, disease outcomes are reported prospectively in RA patients treated in clinical practice. DANBIO covers >90 % of TNFI treated RA patients [15].
We included all RA patients in DANBIO who: 1) were TNFI naïve; 2) started treatment with adalimumab, etanercept or infliximab before 1 July 2007; and 3) had at least two relevant hand x-rays (baseline and follow-up).
A baseline x-ray had to precede the initiation of TNFI by less than 3 months, while the follow-up x-ray had to be obtained >6 months after the baseline x-ray. If available, a pre-baseline x-ray, which preceded both TNFI start and baseline x-ray with >6 months was collected. Five hundred and seventeen patients had three available x-rays, while 930 patients had two x-rays (baseline and follow-up) [16,17].

Ethical approval and patient consent
The individuals in the reference cohort signed an informed consent form. The Danish ethics committee for the City of Copenhagen and Frederiksberg approved the Copenhagen City Heart Study (100.2039/91).
The data from the patient cohort originates from the nationwide Danish DANBIO registry. DANBIO has been approved by The Danish Data Registry since the year 2000 (j. nr. 2007-58-0014 and j.nr. 2007-58-0006), and since October 2006 as a national quality registry by the National Board of Health (j. nr. 7-201-03-12/1). According to Danish law, informed consent and ethical approval were not required for the present study.

Radiographic assessment of patient radiographs
Patient x-rays were available and collected as x-ray films or digital images. To facilitate reading, x-ray films were scanned and converted to a digital format. Reading was performed by an experienced reader, blinded to patient identity and image sequence, according to the modified Sharp method [18].
Annual radiographic progression rates prior to baseline (csDMARD period) and during TNFI treatment (TNFI period) were calculated [16]. Radiographic progression was defined as: a) a change in total Sharp score (TSS) >0; and b) a change in TSS greater than the smallest detectable change (SDC) [19]. Intra-observer intra-class correlation coefficient (one-way random effects model) (ICC) for baseline TSS was 0.96 and for TSS change was 0.35 [20]. SDC was 3.9 TSS units/year [17,21].

DXR-BMD
DXR is a computerized version of the traditional radiogrammetry, previously described in detail [22]. Software automatically detects regions of interest in the three middle metacarpal bones and estimates the bone mineral density (DXR-BMD) based on cortical thickness. Shortterm coefficient of variation has been shown to be 0.28 % and smallest detectable difference 0.0046 g/cm [7,23].
X-ray films from the reference cohort were analysed with the Pronosco X-posure™ v. 2.0 equipment (Sectra, Sweden). DXR-BMD of both hands were obtained and the mean used for analyses (available in 1485/2541 men/ women as 48 men and 77 women were excluded due to missing or unsuitable x-rays). Patient x-rays in a digital format were analysed with the DXR-online system (Sectra, Sweden). The same algorithm was used for calculating DXR-BMD with the Pronosco X-posure™ and DXR-online systems [24].
As indicated in Fig. 1, DXR-BMD analysis or HBL calculation was not possible in a substantial number of patient x-rays due to methodological or disease-related issues. Both hands were measured if possible and the mean used for analyses. If a patient only had DXR-BMD of one hand at any one time point, all analyses were based on that hand (applicable in 38 patients). Three DXR-BMD measurements were available in 135 patients (csDMARD-to-TNFI cohort), while 350 patients had two DXR-BMD measurements (baseline and follow-up; TNFI cohort).

Clinical data
Health assessment questionnaire (HAQ) score and disease activity score in 28 joints based on three variables (DAS28) including C-reactive protein (CRP) were obtained from DANBIO at three visits. For the 135 patients included in the csDMARD-to-TNFI cohort, the three visits were selected to be: closest in time to the pre-baseline x-ray (pre-baseline), closest to the date of TNFI initiation (baseline) and closest in time to the follow-up x-ray (follow-up). For the 215 additional patients included in the TNFI cohort, the baseline and follow-up visits were selected in a similar manner, while the pre-baseline visit was the visit closest to 2 years prior to TNFI.
Patient files were reviewed and data on csDMARD and glucocorticoid treatment in the csDMARD and TNFI period registered.

Statistical analyses
All analyses were performed with R (version 2.15.3) [26]. Analyses were two-sided with a significance level of 0.05.

Reference cohort
Linear regression models for the relation between age and DXR-BMD were fitted for men and women separately. Model fits were compared with the Akaike information criterion (AIC) for non-nested models and analysis of variance (ANOVA) for nested models. Standard graphical tests of model assumptions were performed (plots inspected for linearity, homoscedasticity and normally distributed residuals). From the final models estimated mean annual change in DXR-BMD were calculated for all years of ages from 18 to 89 in both sexes. These estimates constitute reference values for normal HBL/year in the present study.

Patients with RA
HBL is presented as annual absolute (g/cm 2 ) and relative (%) change in DXR-BMD. Increased HBL in an individual patient was defined as a negative HBL/year exceeding the lower 95 % confidence interval (CI) of the normal HBL/year for the matching sex and year of age. For example, a female patient of 54 years would be said to have increased HBL if her HBL/year was lower than -0.0051 g/cm 2 (Additional file 1: supplementary table). HBL was compared between csDMARD and TNFI periods by non-parametric analyses due to a skewed distribution of HBL.
Univariate logistic and linear regression were used to analyse the association between inflammatory activity (assessed with time-averaged CRP, DAS28, 28SJC and 28TJC) and increased and absolute HBL, respectively.
Correlation between HBL and radiographic progression was analysed with Spearman's rho.
Possible predictors for increased HBL were analysed with univariate logistic regression, and significant variables (p < 0.10) included in a multiple logistic regression analysis with backward selection.

Reference cohort
Distribution of DXR-BMD in women and men in the reference cohort is presented in Fig. 2. In women, the best fitting model was DXR-BMD = 0.020 × age -0.00040 × age 2 + 0.0000021 × age 3 (R 2 = 0.54); in men, DXR-BMD = 0.0018 × age -0.0000371 × age 2 (R 2 = 0.27). Estimated mean annual changes in DXR-BMD (i.e. normal HBL/ year) in men and women derived from the models are presented in Additional file 1. Table 1 presents normal HBL/year averaged over 5-year age intervals (10-year intervals in the lowest and highest age groups). In men, the model estimated an increasing HBL/year from 35 years of age onwards reaching a maximum of -0.0047 g/cm 2 /-0.8 % per year at the age of 85. In women, an annual increase in DXR-BMD until 35 years of age was estimated, followed by a continuous HBL that accelerated between 55 and 70 years (greater than -0.0050/-1.0 %).

Patients with RA
A patient disposition is shown in Fig. 1, while demographic, clinical, treatment and radiographic characteristics of included RA patients are summarised in Table 2.
Patients with x-rays unsuitable for DXR-BMD had longer disease duration and more radiographic damage than patients included in the TNFI cohort, but had less functional disability. Other characteristics were similar between cohorts. In the csDMARD-to-TNFI cohort, the median (range) number of days from pre-baseline x-ray to baseline (TNFI initiation) was 607 (180-2989) days, from baseline to baseline x-ray 11 (90-866) days, and from baseline to follow-up x-ray 687 (198-1812) days. Fig. 1 Patient disposition. *Patients were excluded if: a) DXR-BMD could not be analysed due to technical problems (insufficient exposure of x-rays, insufficient positioning of hands) or disease-related factors (prostheses, severe bone damage); or b) DXR-BMD change (i.e. hand bone loss) could not be calculated due to change in acquisition modality between original x-rays (e.g. x-ray film at baseline and digital image at follow-up) or large change in digital image post-processing between x-rays (edge enhancement and change in x-ray resolution). csDMARD Conventional synthetic disease-modifying anti-rheumatic drug, DXR-BMD Bone mineral density estimated by digital x-ray radiogrammetry, TNFI Tumour necrosis factor alpha inhibitors HBL in the csDMARD-to-TNFI cohort In the 135 patients in the csDMARD-to-TNF cohort, prebaseline median (IQR) DXR-BMD was 0.545 (0.474-0.597) g/cm 2 decreasing to 0.516 (0.441-0.578) g/cm 2 at baseline (p < 0.001, Wilcoxon Signed Rank). At follow-up, DXR-BMD had further decreased to 0.504 (0.424-0.557) g/cm 2 (p < 0.001, Wilcoxon Signed Rank). In the csDMARD period, 101 (75 %) patients had increased HBL while increased HBL was found in 79 (59 %) patients in the TNFI period (p = 0.004, McNemars test). Thirty-eight patients had an increased HBL in the csDMARD period that normalised in the TNFI period, while 16 patients had a normal HBL in the DMARD period but increased HBL in the TNFI period. Of the latter 16 patients, 6 patients continued their first TNFI treatment throughout follow-up, while three and seven patients, respectively, withdrew from TNFI treatment or switched to a different biological drug. Corresponding numbers for the 38 patients with normalized HBL during TNFI were 26, 5 and 7.
HBL was significantly lower in the TNFI period (-0.0051 g/cm 2 /year and -1.15 %/year) compared to the csDMARD period (-0.0082 g/cm 2 /year and -1.55 %/year; p < 0.001 for both comparisons, Wilcoxon Signed Rank). Figure 3 illustrates a probability plot for HBL in the csDMARD and TNFI periods.
As was the case for the cut-off of 0, no relation between increased HBL and radiographic progression above the SDC was found in either the csDMARD or the TNFI period (Chi squared = 0.7 and 0.4, respectively). Correlations between TSS status scores and separate erosion score (ES) and joint space narrowing (JSN) scores with DXR-BMD were moderate, while correlations between TSS, JSN and ES changes with HBL were low ( Table 3).

Predictors of HBL during TNFI treatment
In univariate analyses of the TNFI cohort, high baseline DXR-BMD, longer disease duration, immunoglobulin M rheumatoid factor positivity and high DAS28 were associated with increased HBL, while age, HAQ score, CRP, sex, smoking status, type of TNFI, calendar year of treatment initiation and concomitant treatment with prednisolone and MTX were not. In the multivariable model, DXR-BMD (OR 1.005/mg increase (95 % CI 1.003-1.008)) and DAS28 (OR 1.43/unit increase (95% CI 1.15-1.81)) were independently associated with increased HBL. Fraction of explained variation (Nagelkerkes R 2 ) was 0.12. In a model including time-averaged DAS28 during TNFI treatment, baseline DXR-BMD and DAS28 remained independent predictors of increased HBL (R 2 increased to 0.20).

Discussion
In this paper we have established reference values for hand bone mass measured by DXR-BMD. Based on this reference material we estimated sex-and age-related mean HBL/year and used this to define whether HBL was increased or not in individual RA patients. Our main finding was that the majority of patients initiating TNFI treatment in clinical practice had increased HBL and that normalisation of HBL during TNFI treatment was achieved only in a minority of patients. Fig. 3 Cumulative probability plot illustrating HBL in each individual patient in the csDMARD period and TNFI period. The csDMARD period is the period between the pre-baseline and baseline x-rays, while the TNFI period is the period between baseline and follow-up x-rays (see Patients and methods for details). csDMARD Conventional synthetic disease-modifying anti-rheumatic drug, DXR-BMD Bone mineral density measured by digital x-ray radiogrammetry, TNFI Tumour necrosis factor inhibitors Table 3 Correlations between DXR-BMD and radiographic scores in the csDMARD-to-TNFI cohort The large reference cohort (n = 4026) allows precise estimation of sex-and age-related mean HBL/year. The reference cohort was randomly recruited and excluded individuals with inflammatory joint diseases. Individuals with other diseases affecting BMD (osteoporosis, kidney disease, etc.) were not excluded. Theoretically, this could result in overestimation of HBL in the reference cohort. However, a smaller Danish reference cohort excluded such individuals and no systematic differences between DXR-BMD in the two cohorts were found (data not shown) [27].
Based on the distribution of HBL in two early RA cohorts (EURIDISS and BARFOT), the manufacturer of DXR (Sectra, Sweden) defines moderately increased HBL as a negative change in DXR-BMD >0.25 mg/cm 2 / month and <2.5 mg/cm 2 /month (>0.003 g/cm 2 /year and <0.03 g/cm 2 /year) and highly increased HBL as negative change in DXR-BMD >2.5 mg/cm 2 /month (>0.03 g/cm 2 / year) irrespective of age and sex [6,7,10]. According to our reference values, normal HBL in women >50 years of age and men >70 years is >0.003 g/cm 2 /year (corresponding to moderately increased HBL according to the manufacturer), emphasizing the need for an improved definition of increased HBL. No data have demonstrated clinical relevance of the distinction between moderately and highly increased HBL, which is why we chose to simply dichotomise normal vs. increased HBL. The main limitation of our approach is the cross-sectional design of the reference cohort resulting in estimated mean values of HBL, not true HBL measured in individuals. This approach is an established method when no longitudinal studies are available, but should ideally be supplemented with longitudinal data [28]. Melton et al. compared rates of BMD loss by dual-energy x-ray absorptiometry (DEXA) estimated from cross-sectional baseline data with rates obtained from longitudinal assessment and found that cross-sectional data overestimated BMD loss in some skeletal sites (hip and spine) and underestimated the loss in others (radius and ulna) [29]. The background for these findings is unclear and the relevance for DXR is not known.
This study is the first to investigate HBL in established RA during csDMARD treatment and subsequent TNFI treatment. In the BEST and PREMIER trials, patients on MTX monotherapy had a relative HBL/year of -2.6 and -1.9 %, respectively, compared to -1.55 % in our cohort during csDMARD treatment [5,13]. The higher HBL in the RCT monotherapy arms may be explained by frequent use of combination treatment and intra-articular injections in our cohort and the RCT inclusion criteria selecting patients with early and aggressive disease.
As expected, we observed a decrease in HBL (from -1.55 %/year to -1.15 %/year) after initiation of TNFI treatment. The latter is considerably higher than HBL in the combination treatment arms in BEST (-0.6 and -0.9 %/year), but lower than the -1.63 %/year in the combination arm of PREMIER, probably reflecting stricter inclusion criteria in PREMIER than in BEST and our observational study. An HBL of -0.45 %/year was found in 215 patients with established RA receiving different treatments in the Oslo RA register [3]. This considerably lower HBL may be explained by moderate baseline disease activity (DAS28 = 4.0) in contrast to high disease activity in our cohort (DAS28 = 5.3).
Baseline DAS28 was an independent predictor of increased HBL in our cohort in accordance with previous studies documenting a predictive value of baseline DAS28 for future HBL [3,5,13,24]. Interestingly, baseline CRP was not predictive of increased HBL while time-averaged CRP during TNFI treatment was strongly associated with increased HBL. This suggests a key role of sustained inflammation in HBL supporting results from the BEST trial where differences in HBL between treatment groups disappeared after adjustment for changes in disease activity. A similar association is well established with regards to radiographic progression but has not previously been shown for HBL [17,30].
Previously, older age (PREMIER) and postmenopausal status (BEST) has been shown to be predictive of HBL in linear regression analyses [5,13]. In contrast, we found that older age was not associated with increased HBL. This difference stresses the importance of establishing reference data for normal HBL in different age groups and in both sexes.
In accordance with earlier studies, correlations between DXR-BMD and radiographic damage were moderate, and correlations between HBL and radiographic progression low [5-7, 22, 24]. Increased HBL was found in 110 patients with no radiographic progression during TNFI treatment, suggesting a higher sensitivity of HBL in detecting RA-related bone damage. However, we also identified 40 patients who progressed radiographically, but had normal HBL. This dissociation between the two types of bone damage was also seen in a study by Hoff et al. implying that the exact pathophysiological mechanisms may differ and that the traditional evaluation of xrays cannot be substituted by DXR or vice versa [7,31]. In our study, overall radiographic progression was low, reflected in a poor ICC for TSS change [32]. The SDC of 3.9 TSS units is comparable to SDCs reported in major RCTs of patients with established RA [33,34].
A definite limitation of our study is that many patient x-rays could not be analysed. This is a general weakness of DXR, limiting generalisability of our findings as patients with the most advanced and aggressive disease were excluded. A potential limitation of our study is the timing of x-rays, as the selection criteria theoretically allow for a HBL during TNFI treatment to be calculated based on x-rays obtained 90 days prior to baseline and 90 days after baseline. This could lead to an underestimation of the effectiveness of TNFI treatment. However, only two patients in the csDMARD-to-TNFI cohort had a baseline x-ray that preceded TNFI start in combination with a follow-up x-ray less than 300 days after TNFI start, and this is why we consider this bias minimal.
A potential limitation was the lack of information on the use of anti-resorptive treatment in the cohorts. Since previous studies have found no influence of alendronate on HBL, and the more potent anti-resorptive treatments were only marketed in the last part of our study period, we consider this limitation to be of minor significance [9,13].

Conclusions
In this study, we established a reference material for DXR-BMD in the general population and from this reference material we estimated normal sex-and age-related HBL/ year. Applying these reference values we found increased HBL in the majority of RA patients in a cohort from clinical practice initiating TNFI treatment. HBL was only normalised in a few patients during TNFI treatment, and more often in patients who were in time-averaged remission. Competing interests MØ received grants, personal fees and non-financial support from AbbVie; grants and personal fees from BMS; personal fees from Boehringer-Ingelheim; personal fees from Eli Lilly; personal fees and non-financial support from Janssen; grants, personal fees and non-financial support from Merck; personal fees and non-financial support from Pfizer; personal fees and non-financial support from Roche; grants, personal fees and non-financial support from UCB; personal fees from Celgene; personal fees from Sanofi; personal fees from Regeneron; and personal fees from Novartis. Outside the submitted work, MØ has no non-financial competing interests to declare. HL received personal fees from Berlin-Chemie Menarini; grants from Eli Lilly; grants from MSD; grants and personal fees from Roche; grants from Novartis; outside the submitted work, HL has no non-financial competing interests to declare. ORM has received research grants and/or consultancy/speaker fees from Abbott, BMS, MSD, Pfizer, Roche and UCB. Outside of the submitted work, ORM has no non-financial competing interests to declare. All other authors declare that they have no competing interests.
Authors' contributions LMØ participated in the design and conception of the study, contributed to x-ray and clinical data collection in the RA cohort, performed the statistical analyses, interpreted the results and drafted the manuscript. MØ participated in the design and conception of the study, interpreted results and revised the manuscript. TJ participated in the design and conception of the study, performed DXR analyses of the reference cohort x-rays, and helped with interpretation of results and revision of the manuscript. LH contributed to collection and analyses of x-rays in the reference cohort and revised the manuscript. PBM collected reference cohort x-rays, performed DXR analyses of the reference cohort and revised the manuscript. PB read all patient x-rays and helped to revise the manuscript. AT contributed substantially with x-ray collection in the RA cohort and helped to revise the manuscript. UT contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. WB contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. HL contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. UEP contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. AS contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. NG contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. AR contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. JE contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. GK contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. BG contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. ORM contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. DVJ contributed substantially with clinical and x-ray data collection in the RA cohort and helped to revise the manuscript. MLH participated in the design and conception of the study, interpretation of results and revised the manuscript. All authors read and approved the final manuscript.