Article Text

Download PDFPDF

Original research
Added value of multiple autoantibody testing for predicting progression to inflammatory arthritis in at-risk individuals
  1. Frederique Ponchel1,
  2. Laurence Duquenne1,
  3. Xuanxiao Xie1,
  4. Diane Corscadden1,
  5. Farag Shuweihdi2,
  6. K Mankia1,
  7. L A Trouw3 and
  8. Paul Emery1,4
  1. 1Leeds Institute of Rheumatic and Musculoskeletal Medicine, University of Leeds, Leeds, UK
  2. 2Leeds Institute of Health Sciences, Faculty of Medicine, University of Leeds, Leeds, UK
  3. 3Department of Immunology, Leiden University Medical Center, Leiden, Netherlands
  4. 4Leeds NIHR Biomedical Research Centre, Leeds Teaching Hospitals Trust, Leeds, UK
  1. Correspondence to Dr Frederique Ponchel; mmefp{at}leeds.ac.uk

Abstract

Background Predicting progression to clinical arthritis in individuals at-risk of developing rheumatoid arthritis is a prerequisite to developing stratification groups for prevention strategies. Selecting accurate predictive criteria is the critical step to define the population at-risk. While positivity for anti-citrullinated protein antibodies (ACPA) remains the main recruitment biomarker, positivity for other autoantibodies (AutoAbs) identified before the onset of symptoms, may provide additional predictive accuracy for stratification.

Objective To perform a multiple AutoAbs analysis for both the prediction and the time of progression to inflammatory arthritis (IA).

Methods 392 individuals were recruited based on a new musculoskeletal complaint and positivity for ACPA or rheumatoid factor (RF). ELISAs were performed for ACPA, RF, anti-nuclear Ab, anti-carbamylated protein (anti-CarP) and anti-collagen AutoAbs. Logistic and COX regression were used for analysis.

Results Progression to IA was observed in 125/392 (32%) of cases, of which 78 progressed within 12 months. The AutoAbs ACPA, RF, anti-CarP were individually associated with progression (p<0.0001) and improved prediction when combined with demographic/clinical data (Accuracy >77%; area under the curve (AUC) >0.789), compared with prediction using only demographic/clinical data (72.9%, AUC=0.760). Multiple AutoAbs testing provided added value, with +6.4% accuracy for number of positive AutoAbs (AUC=0.852); +5.4% accuracy for AutoAbs levels (ACPA/anti-CarP, AUC=0.832); and +6.2% accuracy for risk-groups based on high/low levels (ACPA/RF/anti-CarP, AUC=0.837). Time to imminent progression was best predicted using ACPA/anti-CarP levels (AUC=0.779), while the number of positive AutoAbs was/status/risk were as good (AUC=0.778).

Conclusion We confirm added value of multiple AutoAbs testing for identifying progressors to clinical disease, allowing more specific stratification for intervention studies.

  • Arthritis, Rheumatoid
  • Anti-Citrullinated Protein Antibodies
  • Rheumatoid Factor

Data availability statement

No data are available.

https://creativecommons.org/licenses/by/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.

Statistics from Altmetric.com

Request Permissions

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

What is already known about this subject?

  • Rheumatoid arthritis is now recognised as a continuum of disease with a preclinical phase identified by the presence of Autoantibodies highlighting a breach of tolerance. Both rheumatoid factor and anti-citrullinated protein antibodies have been associated with progression to inflammatory arthritis from this at-risk, preclinical stage.

What does this study add?

  • Other autoantibodies have been described but have not yet been assessed for their additional value in improving prediction models. Whether combination of autoantibodies may have added value also remains to be investigated.

How might this impact on clinical practice or further developments?

  • Being able to predict the highest risk of progression with better accuracy and/or the imminence of it, would allow the design of clinical trials testing interventions aiming at preventing progression to disease, notably by shortening the follow-up time needed to observe progression.

Introduction

Rheumatoid arthritis (RA) is a chronic disease with substantial impact on the lives of people worldwide, with >1 million patients having severe disability from the disease in Europe alone. Despite the recent progress achieved with ‘treat-to-target’ therapeutic strategies1 and the earlier access to treatment enabled by the European Alliance of Associations for Rheumatology (EULAR)-2010 diagnostic criteria,2 the outcome of RA is still a major concern as well as a financial burden on health services, patients and society with considerable socioeconomic cost.3–5

Over a decade ago, research highlighted a preclinical phase of the disease with the recognition of an inflammatory arthritis continuum (IAC),6 moving the knowledge gap to the earliest stages of disease progression. Although the exact pathogenesis of RA remains unclear, autoimmune processes are believed to play a role and recent research has proposed a series of events underpinning the IAC.6 A breach of tolerance needs to occur (ie, development of autoimmunity) and this assumption has largely been verified by the presence of specific autoantibodies (AutoAb), years before the development of symptoms,7–11 occurring in association with an major histocompatitbility complex/T-cell signalling genetic background.12 13 Systemic autoimmunity is not sufficient to initiate clinical arthritis and a second series of events (the second hit) needs to happen, for which the cells and the possible triggers are not yet known, although several hypotheses have been proposed (implicating cytokines, netosis, pain, osteoclast and/or T-cell or Th17 cells).14–18 This understanding offers an unprecedented opportunity to study ‘at-risk’ individuals and ultimately intervene with the ambition of preventing the progression to arthritis by means of life-style modification or pharmacotherapy.

Several risk stratification models have been developed including demographic, genetic, clinical, imaging, serological and immunological biomarkers with variable predictive value.19–26 Reported predictive biomarkers however have low sensitivity for the prediction of progression to disease stages while specificities are relatively high (>75%). Other studies have reported potential biomarkers, while not modelling for the actual outcome.8 10 11 27–30

The presence of disease specific AutoAbs (notably anti citrullinated protein antibodies, ACPA) prior to disease development remains the best-known risk-factor,7 8 although it is also the main biomarker used to recruit most at-risk cohorts, (including arthralgia and first degree relationship to RA patients). Different ACPA tests also showed different results (positivity and levels), based on the repertoire of peptides used in second and third generation kits.31 Other AutoAbs have been detected years before the onset of inflammatory symptoms including anti-carbamylated protein (anti-CarP), rheumatoid factor (RF) (of the IgM class, further referred to as RF), anti-native and anti-glycated collagen-II AutoAbs.9–11

In this study, we aimed to investigate the added value of testing for the presence of multiple AutoAbs for the prediction of progression to clinical synovitis, in a cohort of at-risk individuals selected on the basis of positivity for ACPA and/or RF associated with a new onset of non-specific musculoskeletal symptoms. Several AutoAbs were explored, notably using in-house ELISA for anti-CarP and anti-collagen AutoAbs. Our data suggest added-value for ACPA, RF, anti-CarP and anti-collagen-II AutoAbs when tested individually either as positive/negative status, continuous levels or dichotomised as high-risk/low-risk groups. Detecting positivity of multiple AutoAbs improved predictive accuracy, which may facilitate more precise selection of study populations as well as determination of high risk of progression in a short time frame to prioritise intervention.

Methods

Subjects and samples

Serum samples were obtained from participants attending an at-risk research clinic at Chapel Allerton Hospital in Leeds between 2008 and 2018.

Briefly, individuals aged ≥18 years, with a new non-specific musculoskeletal joint symptom presenting to their primary care physician or other health professional were referred as previously described.22 23 A new musculoskeletal symptom was defined as pain in any joint and/or musculoskeletal symptom (including but not limited to rotator cuff tendonitis, subacromial bursitis, carpal tunnel syndrome, tendonitis, or epicondylitis), which was not previously reported. Individuals were then recruited to attend the research clinic if they tested positive for ACPA or RF as tested by National Health Service (NHS) services as described below. Exclusion criteria were previous diagnosis of an inflammatory condition, exposure to DMARDs and presence of a swollen joint at first visit. Participants were followed 3 monthly for 1 year and then yearly, until progression to inflammatory arthritis (IA) developed (ie, clinical swelling of at least one joint), evaluated by a rheumatologist/rheumatology research nurse. Progression could occur at any time, and patients had open access to immediate appointment if experiencing a new symptoms. Classification for RA criteria were recorded at study visit only and progression to RA was observed although later during follow-up.

The dataset used in this report was frozen in late 2019 and do not include visits, samples or data collected over the pandemic. We only included participants up to the point in time for which they had visit data. The strategy for data analysis was to use (1) the exact time of progression for participants who progressed at any time during follow-up; (2) for non-progressors, we only included those having at least 12-month follow-up (up to 10 years), while follow-up duration was calculated from baseline until their last available appointment.

A serum sample was collected in clotting tube at inclusion in the cohort, spun for 10 min, no shorter than 30 min and no longer than 1 hour after collection, aliquoted and stored at −30°C.

Autoantibody status and levels

In the absence of absolute standards for AutoAbs assays, arbitrary units (AU) or optical densities (OD) values were used to describe levels observed. Hospital services were used to initially assess positivity (yes/no) for ACPA and/or RF.

RF levels (IgM) were measured by nephrology and considered positive above a cut-off of 20 AU/mL. Anti-nuclear Ab (ANA) were detected as part of our routine assessment and detected by immunofluorescence (homogeneous, nucleolar or specked staining) and reported as negative/positive. Levels of ACPA and RF measured in 15 HCs were all below cut-off and ANA staining were all negative.

ACPA IgG levels were then measured using a seconnd-generation research ELISA (CCP-2, immunoCap Phadia). Cut-off for positivity was set according to manufacturers at 10 AU/ml. Patients with borderline results were retested at 3 months follow-up (mostly confirming negative).

In-house ELISA was developed for IgG autoantibodies for Anti-CarP and anti-native and anti-glycated collagen. In addition, IgG-RF was measured using a commercial ELISA (EuroImmune, Germany) as per manufacturer’s instruction (cut-off 20 RU/mL).

Anti-CarP Abs were measured in the Leiden University Medical Centre (Leiden, The Netherland) as previously described,10 using carbamylated foetal calf serum (10 µg/mL) as antigen. Positivity for Anti-CarP Abs in HC (n=175, median age 51, range 27–65) was previously described32 33 and a cut-off set at 235 UA/mL.

In-house ELISA for the detection of anti-native collagen (CI and CII) and glycated collagen (GLY-CI and GLY-CII) was performed as previously described.11 33 Briefly, collagen CI (Cellsystems) and CII (MD Biosciences) were chemically modified by Maillard reaction in a single batch to generate post-translationally modified GLY-CI/GLY-CII. ELISA plates were coated with 10 µg/mL of GLY-CI/GLY‐CII or native-CI/CII. ELISA OD values were recorded following normalisation to BSA and GLY-BSA (10 µg/mL). Ranges of results in HC were previously established33 (n=98, median age 51, range 27–65) and set at 1.8 and 3 OD/mL for native CI/CII, and 2 and 2.4 for GLY-CI/GLY-CII respectively.

Enough serum was available from only for 152 and 71 individuals for RF-IgG and the collagen related AutoAb testing respectively. ANA were tested in 290 participants. We had missing data for a few samples for the main three autoAbs tested. This affected mainly RF with 11 missing data (status and levels). For ACPA and anti-CarP, this affected 29 (levels only) and 20 cases, respectively (status and level).

Statistical analysis

Frequencies were compared using Pearson’s χ2 test. Numerical variables were not normally distributed (Kolmogorov-Smirnov test) and compared using non-parametric Mann‐Whitney U tests. A value of p<0.05 was considered statistically significant. Where needed, corrections for multiple testing were applied.

Missing data frequencies ranged from none to 31/451 (6.9%) for both clinical and AutoAbs data. Multiple data imputation was performed using five cycles. Logistic regressions were used to model the added value of multiple AutoAb testing for the prediction of progression, using a stepwise Forward approach. A Cox regression was also used to model time to progression. Data were presented using GraphPad Prism V.8 and analysed using SPSS V.26 and /or R V.4.1.2 packages.

Results

Population description

Individuals at-risk of RA, followed since 2008 and up to December 2019, was included in our AutoAbs analysis, when progression to IA could be established over at least 12 months of follow-up, excluding recently recruited patients with less than 12 months of follow-up data, withdrawal by choice. Patients with whom we lost contact (over 12 year) were included to the last appointment available. This allowed us to include 392 participants. Participants had a median age of 51 years (IRQ 18 years, range 18–82) and 70% were women. Participant’s demographic and clinical characteristics (including missing data) according to their progression status are shown in table 1.

Table 1

Overall cohort clinical description (n=392)

We observed progression to clinical IA in 125/392 (32%). The median time of progression was 12 months (IQR 21 months and range 1–113) and the follow-up in non- progressors had a median time of 38 months (IQR 38 and range 12–144). Progression over years is illustrated in online supplemental figure 1. Age was not different between progressors and non-progressors but female gender tended to be more frequent in progressors (p=0.108). A patient-reported case of first degree arthritis in the family was also not associated with progression. Previously reported risk factors from this cohort22 were associated with progression, including smoking history, early morning stiffness, erythrocyte sedimentation rate (ESR), C reactive protein and painful joint counts (all p<0.001). Comorbidities such as diabetes were not associated with progression.

There was randomness of missing data at low frequencies (maximum 5.1% missing data excluding RF-IgG, ANA and anti-collagen, see below), therefore, missing data were imputed for both clinical as well as for AutoAbs. Unadjusted OR for each variable were calculated from the pooled 5 cycles of imputation and were not significantly different from the original dataset.

AutoAbs determination

Because of limited amount of serum, all AutoAbs could not be tested in all samples (missing tests detailed in table 2), notably for the anti-collagen AutoAbs and RF-IgG, while ANA were not performed routinely in 36% of participants, resulting in 7 AutoAbs tested.

Table 2

AutoAbs individual predictive performance (n=392)

Most participants were ACPA+ (206/392, 52%), accounting for individual with initial borderline levels, found negative on retesting, and 61 were ACPA+only. 167 participants were RF+ (figure 1A, median 65 UA/mL, IQR 140 UA/mL) of which, 39 participants were RF+only. RF-IgG and ANA were positive in 50% and 20% of samples tested respectively. Positivity for anti-CarP was observed in 126/372 (33.9%) of participants, 10 being anti-CarP+only. Anti-native-CI and CII AutoAbs were detected at low frequencies (9.9% and 29.6%, respectively) while anti-GLY-CII (50.7%) was more frequent but not anti-GLY-CI (24%).

Figure 1

Positivity (A), as number of patients and (B) levels of AutoAb in at-risk individuals. (A) Bars represent the number of individual tested positive (black) or negative (white) for each individual. X2 tests were used individually to assess associations with progression (***p<0.0001, #p<0.100). ACPA n=363, RF n=381, anti-CarP n=372, RF-IgG n=152, ANA n=290, all anti-collagen autoAbs n=71. (B) Violin plots represent the distribution of AutoAb levels observed. Solid lines across the plot indicate the positivity cut-off for each test. Medians and quartiles of distribution are indicated by dotted and dashed lines within violin plot respectively. MWU tests comparing levels were used individually to assess associations with progression (*p<0.05, ***p<0.001, #p<0.100). ACPA, anticitrullinated protein antibody; ANA, anti-nuclear Ab; RF, rheumatoid factor.

Using positivity status (figure 1A), ACPA, RF and anti-CarP were highly associated with progression (table 2A, all p<0.0001). Although the assays could only be performed in a small number of samples (n=71), anti-GLY-CII and anti-native-CII showed a trend for association with progression (p=0.075 and p=0.111, respectively). RF-IgG, ANA and Anti-CI AutoAbs were not associated with progression.

Using continuous levels of AutoAbs (figure 1B), there was significantly higher levels of ACPA (irrespective of the test used), RF and anti-CarP AutoAbs in progressors (table 2B, p<0.0001), while levels were significant higher for anti-native CII (p=0.048) but not for anti-GLY-CII.

Results presented below account for the cohort of 392 patients following imputations for ACPA, RF and anti-CarP, while RF-Ig and ANA AutoAbs were not pursued further due to lack of association with progression and anti-collagen AutoAbs as it could only be tested in 71 patients.

Prediction of progression: individual AutoAbs

Altogether, ACPA, RF, anti-CarP and anti-CII AutoAbs were taken forward, while RF-IgG, ANA, and CI AutoAbs were no longer included in the analysis as not associated with progression. The rate of progression in ACPA+, RF+ and anti-CarP+individual are illustrated in figure 2A, confirming clear separation between status.

Figure 2

Modelling the predictive value of autoAbs in at-risk individuals. (A) Individual time of progression with respect to positivity for the three main autoAbs tested (red positive, blue negative). (B) Logistic regression: AUC for prediction models using clinical data and multiple autoAbs based on status, levels, high and low risk groups or autoAbs count. (C, D) Cox regression: Survival and AUC for prediction models using clinical data and multiple autoAbs based on status, levels, risk groups or autoAbs count. (B–D) Black clinical data only; green clinical data+status; red clinical data+levels; blue clinical data+risk group; purple clinical data+autoAbs count. ACPA, anticitrullinated protein antibody; AUC, area under the curve; RF, rheumatoid factor.

ACPA positivity performed best, the OR being 10.95 (table 2A), with best sensitivity (86%) and 90% negative predictive value (NPV) notably compared with RF (81%) and anti-CarP (76%). RF and anti-CarP had high OR>4.5 compared with other variable (table 1).

An analysis performed on continuous AutoAbs levels (table 2B) confirmed the predictive performance of the different AutoAbs with area under the curve (AUC) values ranging from 0.702 for RF to 0.770 for ACPA. This allowed for levels to be dichotomised for a more practical high versus a low risk of progression, using a cut-off set at 80% specificity for progression. Individual predictive performances for each AutoAbs were calculated for high vs low risk (table 2C). ORs were not better than when using AutoAbs status but AUC were improved for anti-CarP (0.675) while but not for ACPA or RF. All positive predictive value (PPV) were relatively similar (range 54%–60%), while NPV were best for CCP-res assays (81%).

AUC is individually presented for ACPA/RF/anti-CarP in online supplemental figure 2 for status, levels and risk groups.

Prediction of progression: added-value of Individual autoAbs

To demonstrate the possible added value of AutoAbs testing, we first established which demographic and clinical variables were predictive of progression (from eight parameters routinely recorded). A regression model built from the clinical data retained 4-variables (table 1, AUC=0.760 and displayed in online supplemental figure 3), including a genetic risk factor (shared epitope SE), an environmental risk factor (smoking), pain (tender joint count in 78 joints, TJC78) and an inflammation marker (ESR) which reproduced previously reported results22 and allowed 72.9% of individual’s outcome to be predicted correctly (Accuracy) and all four variable being highly predictive (all p<0.0001). This 4-variables model was considered the reference for further comparison.

We evaluated the added value of AutoAbs individually (table 3 and online supplemental figure 3) for status (A), levels (B) and risk categories (C). All models retained the four original clinical variables with each AutoAb, suggesting that they all had individual added value, improving accuracy (+4.1% to +4.4%) over the reference model and better AUC (+0.029 to +0.076). Overall, the AutoAb status with the best prediction accuracy were ACPA or anti-Carp (7.3%), best specificity-PPV was with anti-Carp (90% and 70%) although the AUC=0.789 was less good than ACPA (0.836) or RF (0.800).

Table 3

Added-value of individual autoAbs (n=392)

For continuous levels (table 3B), ACPA and anti-CarP but not RF, demonstrated added value over the reference model. The best performing model was with anti-CarP with+5.5% accuracy but the best AUC was for ACPA (0.824).

Using AutoAbs risk categories (table 3C), the best models were similar for ACPA or anti-Carp (+5.4% accuracy) but the best AUC was for ACPA (0.832). All model had similar specificity (anti-CarP being the best with 92%) while ACPA had highest sensitivity (57%).

Added value of multiple autoantibody testing

We used combination ACPA with RF and anti-CarP to evaluate the added value of multiple AutoAbs testing. A progression rate is presented in online supplemental figure 4 with respect to combined positivity for the three autoAbs.

First, Status were combined in a single variable, counting the number of positive AutoAbs. Data suggested significant association of higher number of positive AutoAbs with progression (p<0.0001) with triple negativity (123/392 (28% of the cohort) poorly associated with progression (5/123 (4.0%)), 1 AutoAb+ (117/392 (30% of cohort)) showing increased progression rate (28/117 (23%)), 2 AutoAbs+ (83/392 (21% of cohort)) a further increase (44/83 (53%)) and triple positivity (69/392 (17.6% of the cohort)) the highest progression rate (48/63 (76%). The number of positive AutoAbs had an individual OR=3.422 and AUC=0.813.

The AutoAbs count was then analysed with the variables of the reference model to address its added-value (table 4A, figure 2B). The combination provided a clear improvement on accuracy (+6.4%) and AUC=0.852). To account more for the importance of each AutoAbs to the prediction, models were then constructed using status, levels and risk groups. Using status, the model retained both ACPA and RF but removed anti-CarP. It increased slightly the prediction (77.8% accuracy) of the models with only ACPA or anti-Carp (table 3, 77.3%), while showing an increased AUC=0.855, but not better specificity over anti-Carp-alone (90%). There was added accuracy (+5.4%) of multiple AutoAbs testing using continuous levels keeping ACPA and anti-CarP above the model using anti-CarP levels alone (77.7%), as well as improved AUC=0.832. Using risk categories, the model retained again ACPA and anti-CarP, and showed good accuracy (79.1%) although slightly bellow AutoAb counts (−0.2%) but above autoAb status (+1.3%).

Table 4

Added value of multiple autoAb testing: logistic and COX regression T(est of proportional hazards (PH) models overtime were valid (p>0.05) for all predictors)

We then performed an analysis of time to progression using a COX regression (table 4B). We chose to look at imminent progression within 12 months which did not include patients who progressed late (n=344, as illustrated in online supplemental figure 1). The proportional hazards assumption of Cox model was checked and found satisfied (p>0.05). The reference model selected the different variables (SE, ESR EMS and TJC78) and had an AUC=0.725 (table 1). The COX models did not retain smoking but showed clear increase in AUCs (table 4B, figure 2C). The status and risk category models kept ACPA and RF and were as good as counting autoAbs model (all same AUC=0.778). The best improvement in AUC was observed for AutoAbs levels (AUC=0.779, figure 2D) with ACPA and anti-CarP, suggesting value for utilising levels for predicting time of imminent progression.

Discussion

Multiple AutoAbs have been associated with RA, including triple positivity for ACPA/RF/anti-CarP34 35 and more recently with its pre-clinical stages.7–11 Here, we show the individual and combined value of AutoAbs used for RA classification for the prediction of progression to clinical disease. In this large cohort of 392 individuals selected on the basis of having a new non-specific musculoskeletal complain and being ACPA+and/or RF+, progression to IA was associated with ACPA, RF, anti-CarP and trends for anti-collagen-II Abs (likely due to small numbers), while not with anti-collagen-I and ANA. These associations were individually predictive improving the performance of the prediction based solely on demographic and clinical characteristics.22 Levels (although not for RF) consistently provided higher specificity and better AUC, and PPV, while sensitivity was more variable and NPV were similar. Here, the combination of AutoAbs levels with demographic and clinical data showed improvement over an individual AutoAb model, while dichotomising levels into high/low risk groups was better than using positive status.

Our data aligned with a meta-analysis of 12 studies,36 showing evidence of the value of triple positivity, ACPA, RF and anti-CarP, in identifying individuals at-risk from the healthy population (however not including demographic and clinical data in the models). We further demonstrate that triple positivity as well as analysis by status, levels or risk groups all suggested clear added value of using anti-CarP AutoAb for a better prediction of progression overall as well as its timing.

Participants included in our cohort based on NHS-ACPA positivity showed a high rate of false positive compared with research ELISA test (also previously observed).31It may be beneficial to use a higher cut-off (at 10 AU/mL for example) for recruitment of individual with non-specific musculoskeletal symptoms. This would allow >75% of false positive cases, not to be referred for clinical follow-up. On the other hand, these individuals still present with a new MSK pain complaint and may actually offer an interesting ACPA-negative research control group. This had notably allowed for RF+only participants (n=39, 8 progressed) and anti-CarP+only (n=10, 1 progressed) participants to be identified (see online supplemental figure 2 for illustration). Altogether, 123 participants had a triple negative status and five progressed (1.2%) compared with 26/110 (23.6%) in the group with at least one auto-Ab positive (p<0.0001). Importantly for the 186 ACPA-negative participant 14 of the 17 that progressed were either RF+ (n=9) or anti-CarP+ (n=1) or both (n=4), leaving three progressors currently not identified by an autoantibody (although one was positive for anti-Gly-CII+).

Currently, many biomarker have been proposed in individual at-risk of RA. From a single blood test, our study has confirmed the previously reported better individual predictive value of ACPA (using second or third generation tests,37 38 while still demonstrating that combinations of autoantibodies are more informative that each autoantibody alone. In addition to serological testing, other biomarkers such as high resolution imaging (ultrasound and MRI)39–41 or T-cell subsets23 24 or the combination of some, may still provide increased predictive value. The impact of certain biomarkers may also be greatest in those with ‘imminent’ arthritis, given most individuals develop subclinical joint and tendon inflammation prior to the onset of clinical joint swelling41 while other may provide additional value for understanding pathogenesis.

Limitations to the study, despite a large number of participants was the uniqueness of recruitment criteria, which currently limits any replication. There are different recognised populations considered at-risk of RA, notably as most recently redefined by a EULAR taskforce.42 These include seropositive arthralgia, clinically suspected arthralgia (CSA), first degree relatives of RA patients and ACPA+individuals with non-specific MSK symptoms. In Leeds, we chose to recruit the later, as at-risk individuals often initially present to primary care clinicians while CSA requires specialist assessment (ie, secondary care). Furthermore, there is a much higher overall rate of progression in these individual than reported in CSA, while ACPA+positive RA is usually more severe compared with seronegative disease.43 44 As criteria are needed to define an at-risk population and to balance the specificity of the recruitment with the number of cases needed to develop models. The overall cohort tested over 9000 potential participants since 2008, while about >750 were included being ACPA+in line with reported prevalence of ACPA-positivity.45 The populations and models developed for ACPA+arthralgia22 46 are therefore not directly transferable to cohort of patients with CSA (as defined by EULAR taskforce42 47 48), or in first degree relatives of patients with IA,49–51 which are the other main characteristics used for selecting at risk individuals. In addition, progression rate in these different groups are different and if 32% for IA in our cohort, progression to RA was only observed following from the development of IA in less than 25% of cases.

In conclusion, our data confirm the value of multiple AutoAbs testing for three particular antigens in preclinical RA, while suggesting that others may also have value (anti-collagen-II AutoAbs notably) and allowing for the exclusion of some (ANA, RF-IgG, anti-Collagen-I but not II). A baseline assessment of multiple AutoAbs may, therefore, inform a follow-up strategy for individual complaining of non-specific musculoskeletal symptoms. A research strategy based on an observational design may direct participants for 3 monthly vs annual review, using an overall prediction of progression (logistic models) based on counting AutoAbs or using risk groups. Alternatively, a strategy aiming at modifying the risk of progression (ie, intervention) may be better tested using a stratification based on time to progression (Cox regression) reducing the need for long-term follow-up.

Data availability statement

No data are available.

Ethics statements

Patient consent for publication

Ethics approval

The at-risk of RA study (coordinated Programme to Prevent Arthritis) was approved under references REC-06/Q1205/169. Healthy controls (HC) were also included (REC-09/H1307/98). All participants provided written informed consent for their blood samples to be used in research.

Acknowledgments

XX received support from a joint PhD scholarship between the University of Leeds and the China Research Council (No. 201708330243). LAT received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 724517).

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors Design of the study: FP and PE. Data acquisition, data analysis: FP, LD, XX, FS, LAT and DC. Drafting of the manuscript: FP, LAT, KM and LD. Review/editing of the manuscript: all authors.

  • Disclaimer The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health

  • Competing interests LAT is mentioned as an inventor on a patent describing a method to detect anti-CarP antibodies.

  • Patient and public involvement statement Patient and Public were involved in the study design, recruitment and samples collection. Biomarker research was included in the study design but the lab work was not directly involving pPPI.

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

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.