# Anti-citrullinated protein antibody specificities and pulmonary fibrosis in relation to genetic loci in early rheumatoid arthritis

Mikael Brink, Lotta Ljung, Monika Hansson, Johan Rönnelid, Rickard Holmdahl, Karl Skriner, Guy Serre, Lars Klareskog, Solbritt Rantapää-Dahlqvist

► To cite this version:

Mikael Brink, Lotta Ljung, Monika Hansson, Johan Rönnelid, Rickard Holmdahl, et al.. Anti-citrullinated protein antibody specificities and pulmonary fibrosis in relation to genetic loci in early rheumatoid arthritis. Rheumatology, 2022, 61 (12), pp.4985-4990. 10.1093/rheumatology/keac280. hal-04220978

HAL Id: hal-04220978

https://hal.science/hal-04220978v1

Submitted on 28 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Concise report

Anti-citrullinated protein antibody specificities and pulmonary fibrosis in relation to genetic loci in early rheumatoid arthritis

Mikael Brink ID <fnref n="1" />, Lotta Ljung ID <fnref n="1" />, Monika Hansson², Johan Rönnelid ID ³, Rickard Holmdahl ID ⁴, Karl Skriner ID ⁵, Guy Serre⁶, Lars Klareskog² and Solbritt Rantapää-Dahlqvist ID <fnref n="1" />

Abstract

**Objectives.** Pulmonary manifestations in RA are common comorbidities, but the underlying mechanisms are largely unknown. The added value of a multiplex of ACPA and genetic risk markers was evaluated for the development of pulmonary fibrosis (PF) in an inception cohort.

**Methods.** A total of 1184 patients with early RA were consecutively included and followed prospectively from the index date until death or 31 December 2016. The presence of 21 ACPA fine specificities was analysed using a custom-made microarray chip (Thermo Fisher Scientific, Uppsala, Sweden). Three SNPs, previously found related to PF were evaluated, rs2609255 (FAM13A), rs111521887 (TOLLIP) and rs35705950 (MUC5B). ACPA and genetic data were available for 841 RA patients, of whom 50 developed radiologically defined PF.

**Results.** In unadjusted analyses, 11 ACPA specificities were associated with PF development. In multiple variable analyses, six ACPA specificities were associated with increased risk of PF: vimentin (Vim)60–75, fibrinogen (Fib)β62–78 (72), Fibα621–635, Bla26, collagen (C)II359–369 and F4-CIT-R (*P* < 0.01 to *P* < 0.05). The number of ACPA specificities was also related to PF development (*P* < 0.05 crude and adjusted models). In multiple variable models respectively adjusted for each of the SNPs, the number of ACPA specificities (*P* < 0.05 in all models), anti-Vim60–75 (*P* < 0.05, in all models), anti-Fibβ62–78 (72) (*P* < 0.001 to *P* < 0.05), anti-CII359–369 (*P* < 0.05 in all models) and anti-F4-CIT-R AQ4 (*P* < 0.01 to *P* < 0.05), anti-Fibα621–635 (*P* < 0.05 in one) and anti-Bla26 (*P* < 0.05 in two) were significantly associated with PF development.

**Conclusion.** The development of PF in an inception cohort of RA patients was associated with both presence of certain ACPA and the number of ACPA specificities and risk genes.

**Key words:** RA, pulmonary fibrosis, autoantibodies, ACPA specificities, genetic loci

Rheumatology key messages

* Development of pulmonary fibrosis in rheumatoid arthritis is associated with certain ACPA fine specificities.

* The number of positive ACPA specificities increased the risk for development of pulmonary fibrosis.

* Presence of the autoantibodies associated with pulmonary fibrosis was unrelated to the risk genes.

Introduction

RA is a chronic systemic inflammatory disease primarily
affecting the joints, but with a noticeable frequency of

severe extra-articular manifestations, especially lung dis-
ease [1]. Interstitial lung disease (ILD) with fibrosis in RA
shares several characteristics with idiopathic pulmonary
fibrosis (IPF), with a progressive course and poor survival

<footnote n="1">¹Department of Public Health and Clinical Medicine/Rheumatology, University Hospital, Umeå, ²Rheumatology Unit, Department of Medicine, Karolinska Institute at Karolinska University Hospital, Stockholm, ³Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, ⁴Medical Inflammation Research, Department of Medical Biochemistry and Biophysics, Karolinska Institute, Solna, Stockholm, Sweden, ⁵Department of Rheumatology and Clinical Immunology, Charite Universitätsmedizin, Berlin,</footnote>

<footnote>Germany and ⁶Institut Toulousain des Maladies Infectieuses et Inflammatoires, UMR 1291 Inserm, 5051 CNRS, Université de Toulouse, Toulouse, France</footnote>

<footnote>Submitted 14 February 2022; accepted 24 April 2022</footnote>

<footnote>Correspondence to: Solbritt Rantapää-Dahlqvist, Department of Public Health and Clinical Medicine, Rheumatology, University Hospital, SE-901 85 Umeå, Sweden. E-mail: solbritt.rantapaa.dahlqvist@umu.se</footnote>
[2]. The lungs have been proposed as a site of disease initiation with increased amounts of citrullinated peptides in lung tissue of smokers [3]. Thus, smoking in the presence of HLA-DRB1-shared epitope alleles may lead to the development of ACPA-positive RA [4]. However, infections and inflammation also increase citrullination of proteins, and citrullinated proteins have been found in lung tissue besides that from RA patients, in patients with idiopathic ILD [5, 6]. A broader ACPA repertoire was presented in RA-ILD suggesting a role for ACPA in the pathogenesis of ILD [7].

A strong genetic basis for IPF has been demonstrated in studies of familial aggregations [8] and in genome-wide association studies (GWASs) of the general population [9, 10] and of RA patients [11]. In a recent published study we found three SNPs, rs35705950 (MUC5B), rs111521887 (TOLLIP) and rs2609255 (FAM13A), to be significantly associated with pulmonary fibrosis (PF) in RA patients [12]. No significant relationship was found between anti-CCP2 antibodies and PF.

Here, we report the presence of ACPA fine specificities in relation to the development of PF and the three previously identified PF risk SNPs, analysed in our inception cohort of patients with RA followed prospectively within the catchment area of northern Sweden.

## Methods

The subjects and collection of clinical patient data and genetic analyses, as well as autoantibody detection, have been described previously [12, 13]. Briefly, an inception cohort of patients diagnosed with early RA (eRA) (symptomatic <12 months before diagnosis) according to the American Rheumatism Association classification criteria [14] were consecutively included in the study at the time of diagnosis (index date) between 1 January 1996 and 31 December 2016 at the five rheumatology clinics in northern Sweden. Clinical data, e.g. the 28-joint count and disease activity score (DAS28), were registered systematically and recorded at index date and at 6, 12, 18 and 24 months after diagnosis. Pharmacological treatments defined as corticosteroids and conventional synthetic DMARDs (csDMARDS; methotrexate, sulfasalazine, chloroquine, leflunomide, azathioprine, ciclosporin, mycophenolate mofetil and injectable gold salts) and biologic DMARDs (bDMARDs; abatacept, adalimumab, anakinra, etanercept, infliximab, rituximab, tocilizumab) were registered at the index date and continuously for the first 24 months. Smoking habits were registered as smoking ever vs non-smoking and as current smoker. Examinations and evaluations of the radiographs of the lungs were performed as previously described [12]. Data on high resolution CT (HRCT) examinations were collected during a 20 year period limiting the evaluation of the pulmonary manifestations to PF including reticular pattern, honeycombing or traction bronchiectasis of variable degree and with ground-glass in a few cases [15]. Thus, no further diagnostic evaluation was performed to identify and separate manifestations of ILD. Of the

original included patients (n = 1184) in 50 patients no information on radiologic pulmonary examinations could be retrieved; these patients were excluded, leaving 1134 patients for further analyses [12]. In the remaining eRA cohort, the patients had received a diagnosis of PF at inclusion or during follow-up, in 96% of the cases based on HRCT findings and in two on plain X-rays after a mean (s.d.) of 6.2 (4.7) years after onset of RA.

Of the included 1134 pulmonary evaluated eRA patients, plasma samples were analysed in 841 patients for IgG-specific ACPA fine specificities, using the delta value from analysis of both the citrullinated peptides and their arginine-containing counterpart using a custom-made microarray chip (Thermo Fisher Scientific, ImmunoDiagnostics, Uppsala, Sweden). Autoantibodies against the following peptides were analysed: α-enolase peptide 5-21 (CEP-1), collagen type II (CII359-369, F4-R-CIT, F4-CIT-CIT and F4-CIT-R, all in triple helical conformation [16, 17]), fibrinogen (Fib)α×36-50, Fibα563-583, Fibα580-600, Fibα621-635, Fibβ36-52, Fibβ60-74, Fibβ62-78(72), Fibβ62-78(74), filaggrin (Fil307-324), vimentin (Vim)2-17, Vim60-75 and hnRNP-A3 peptides (Bla26, Pept-1, Pept-5, PeptZ1 and PeptZ2) (presented in Supplementary Table S1, available at Rheumatology online). The methodology of analysing the antibodies is semi-quantitative, and thus dichotomized values of the antibodies are presented based on cut-off values set at the 98th percentile of 477 healthy controls; for more details see [13, 18]. At index date the analyses of RF used routine laboratory methods and those of anti-CCP2 antibodies were according to the manufacturer (EuroDiagnostica, Malmö, Sweden). All autoantibody analyses were done in plasma samples collected at the index date and preserved at -80°C until analysed. The antibody frequencies in the cohort with pulmonary examinations corresponded to the frequencies in the originally ACPA analysed cohort [13] (Supplementary Table S2, available at Rheumatology online).

GWAS genotyping of DNA samples from cases and controls was performed using the Global Screening Assay (GSA; Illumina, San Diego, CA, USA) to analyse 571151 genome-wide single-nucleotide polymorphisms (SNPs) at deCode Genetics (Reykjavík, Iceland). The three loci containing SNPs found to be significantly associated with PF development were used for further analyses: rs35705950 (MUC5B), rs111521887 (TOLLIP) and rs2609255 (FAM13A) [12]. HLA-shared epitope (HLA-SE) (0401/0404/0405/0408 or 0101) was identified from GWAS data using imputation from a European reference panel; for details see [12].

## Statistics

Statistical analysis was performed using SPSS software (v. 27.0, IBM Corp, Armonk, NY, USA) and R software version 4.1.2 (R Foundation for Statistical Computing, Vienna, Austria [19]). Descriptive data were summarized and presented as proportions with means and standard deviation. Frequencies were compared using the chi-square test, and for continuous data Student's t-test
was used. Associations between PF and possible predictors including genetic markers were analysed using logistic regression analysis and presented as the odds ratio (OR) with 95% confidence interval (CI). Frequency of ACPA positivity was calculated for each specificity and the number of ACPA specificities for each case. A crude P-value is presented for the simple model based on the same subsets of samples as the multiple variable models for which the adjusted P-value is presented. The significance level was set at *P* < 0.05. For missing data, imputations have been performed using linear/logistic regression based on age, sex, smoking, methotrexate use from index date (n = 33) and DAS28 at inclusion (n = 15). DAS28 was also summarized during the first 24 months to area under curve for DAS28 (DAS28-AUC24).

## Ethics

The study complies with the Declaration of Helsinki and the Regional Ethics Committees at Umeå University, Sweden approved this study (No. Dnr 2017-432-32M, 2019-02039) and the patients signed an informed consent to participate.

## Results

### Presence of ACPA fine specificities in eRA with pulmonary fibrosis

Of the 841 eRA patients included in our inception cohort and analysed for presence of ACPA fine specificities, 50 (5.9%) had PF at the time of diagnosis or developed PF during follow-up. Of the 21 analysed ACPA specificities, 11 antibodies were significantly associated with PF in unadjusted logistic regression analyses, namely anti-Fil307-324 (CCP1), anti-Vim60-75, anti-Vim2-17, anti-CEP1, anti-Fibβ62-78 (72), anti-Fibα621-635, anti-Bla26, anti-Pept5, anti-CII359-369, anti-F4-CIT-CIT and anti-F4-CIT-R antibodies, as well as the number of ACPA specificities, but not presence or level of anti-CCP2 (Table 1). In multiple variable analyses adjusted for sex, age at RA diagnosis, smoking ever status and presence of RF, six of the ACPA specificities remained significantly associated with increased risk of PF development, namely anti-Vim60-75 [OR 2.13 (95% CI 1.07, 4.50), *P* < 0.05], anti-Fibβ62-78 (72) [OR 2.43 (95% CI 1.20, 4.67), *P* < 0.01], anti-Fibα621-635 [OR 2.12 (95% CI 1.12, 3.85), *P* < 0.05], anti-Bla26 [OR1.92 (95% CI 1.05, 3.51), *P* < 0.05], anti-CII359-369 [OR 1.92 (95% CI 1.04, 3.56), *P* < 0.05] and anti-F4-CIT-R antibodies [OR 2.57 (95% CI 1.38, 4.79), *P* < 0.01] (data not shown). A higher number of ACPA specificities in an individual increased the risk for PF development, with the mentioned adjustments [OR 1.07 (95% CI 1.01, 1.13), *P* < 0.05]. The relationships between certain of the ACPA specificities and PF did not result from higher frequencies of these antibodies (see Supplementary Table S2, available at Rheumatology online).

### Presence of ACPA fine specificities and SNPs in eRA with pulmonary fibrosis

In subsequent models, the ACPA specificities that remained significantly associated with PF after adjustments were further analysed respectively adjusted for the three PF risk SNPs (rs35705950, rs111521887 and rs2609255) (Table 2). These SNPs were, as we have shown, significantly associated with PF in this patient cohort (Table 1, and ref. [12]). In the multiple variable models (adjusted for sex, age, ever smoker status, methotrexate at index date and RF status) and including each separately of the three SNPs, antibodies against Vim60-75, Fibβ62-78 (72) and F4-CIT-R and the number of ACPA specificities remained significantly associated with PF development (Table 2; *P* < 0.001 to *P* < 0.05). However, with the adjustments the association between anti-Fibα621-635 and PF was attenuated (Table 2).

No association was found between either of the three SNPs and positivity for any of the ACPA fine specificities (data not shown).

## Discussion

This exploratory study of a consecutively included inception cohort of patients with eRA is the first study to analyse both PF-related genetic risk loci and ACPA fine specificities for the PF development in RA patients. Associations were found between certain ACPA specificities and PF risk, irrespective of the risk genotype. The number of ACPA specificities was also associated with PF development, independent of presence of the selected risk genotypes. However, we did not find a significant association between PF and anti-CCP2 positivity or levels, but the findings of ACPA fine specificities are in line with a previous publication where the number of high levels of ACPAs was associated with the degree of PF [7]. Higher levels of all ACPA specificities analysed by Giles et. al. [7] showed significant association with ILD, including specificities to fibrinogen and vimentin, as shown in the present work although not using exactly the corresponding peptides. Further, Giles et al. presented an association between the level of antibodies against α-enolase and ILD, while no association of this antibody could be found with PF in our cohort (*P* = 0.067, *P*-value in unadjusted model). Higher frequency or higher level, respectively, of anti-filaggrin antibodies [Fil307-324 (CCP1)] and of anti-enolase 1A cyclic and anti-CEP1 antibodies showed significant associations in each of the two studies with the outcome measure, PF and ILD, respectively. In our study the significant association was lost after the mentioned adjustments (for anti-filaggrin antibodies *P* < 0.030 and 0.065 and for anti-CEP1 0.031 and 0.057 in unadjusted and adjusted models, respectively). An expanded ACPA repertoire was suggested by Giles et al. as a distinct feature of antibody associations with ILD, e.g. expanded ACPA specificities were most strongly associated with features of fibrosis, while anti-CCP2 level
<table><caption>TABLE 1 Demographic and clinical data of the included patients with early RA, with and without pulmonary fibrosis</caption><thead><tr><th>Variable</th><th>PF<br>(n = 50)</th><th>No PF<br>(n = 791)</th><th>OR (95% CI)</th><th>P-value</th></tr></thead><tbody><tr><td>Age at RA diagnosis, mean (s.d.), years</td><td>63.90 (10.35)</td><td>56.98 (13.85)</td><td>1.05 (1.02, 1.07)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Female, n (%)</td><td>29 (58.0)</td><td>548 (69.1)</td><td>0.61 (0.34, 1.10)</td><td>0.098</td></tr><tr><td>Ever-smoker, n (%)</td><td>38 (77.6)</td><td>505 (64.7)</td><td>1.90 (0.95, 3.77)</td><td>0.068</td></tr><tr><td>Current smoker, n (%)</td><td>11 (22.4)</td><td>167 (21.4)</td><td>0.97 (0.47, 1.83)</td><td>0.93</td></tr><tr><td>Anti-CCP2 antibody<sup>+</sup>, n (%)</td><td>39 (78.0)</td><td>535 (67.5)</td><td>1.69 (0.85, 3.35)</td><td>0.135</td></tr><tr><td>Anti-CCP2 level, mean (s.d.)</td><td>412.81 (523.80)</td><td>316.52 (624.42)</td><td>1.00 (1.00, 1.001)</td><td>0.337</td></tr><tr><td>RF<sup>+</sup>, n (%)</td><td>43 (86.0)</td><td>580 (73.1)</td><td>2.25 (0.996, 5.077)</td><td><strong>0.051</strong></td></tr><tr><td>HLA-shared epitope, n (%)</td><td>37 (74.0)</td><td>583 (73.7)</td><td>1.04 (0.85, 1.27)</td><td>0.72</td></tr><tr><td>DAS28 on index date, mean (s.d.)</td><td>4.76 (1.49)</td><td>4.69 (1.43)</td><td>1.03 (0.85, 1.24)</td><td>0.78</td></tr><tr><td>DAS28-AUC until PF or 24 months, mean (s.d.)</td><td>77.33 (33.17)</td><td>83.23 (22.18)</td><td>0.989 (0.977, 1.001)</td><td><strong>0.082</strong></td></tr><tr><td>DAS28-AUC<sub>24</sub>, mean (s.d.)</td><td>90.00 (20.05)</td><td>83.23 (22.18)</td><td>1.014 (1.001, 1.027)</td><td><strong>0.038</strong></td></tr><tr><td>Treatment at index date, n (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Glucocorticoids</td><td>28 (56.0)</td><td>463 (59.8)</td><td>0.85 (0.48, 1.52)</td><td>0.587</td></tr><tr><td>Methotrexate</td><td>35 (70.0)</td><td>630 (82.0)</td><td>0.51 (0.27, 0.97)</td><td><strong>0.039</strong></td></tr><tr><td>SNP</td><td></td><td></td><td></td><td></td></tr><tr><td>rs2609255 G, MAF</td><td>0.35</td><td>0.23</td><td>1.79 (1.17, 2.74)</td><td><strong>0.008</strong></td></tr><tr><td>rs111521887 G, MAF</td><td>0.32</td><td>0.21</td><td>1.78 (1.15, 2.75)</td><td><strong>0.009</strong></td></tr><tr><td>rs35705950 T, MAF</td><td>0.22</td><td>0.12</td><td>2.17 (1.25, 3.75)</td><td><strong>0.005</strong></td></tr><tr><td>Frequency of ACPA positivity (95%CI)</td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-Fil307-324 (CCP1)</td><td>72.0 (57.9, 82.8)</td><td>55.9 (52.3, 59.2)</td><td>2.03 (1.08, 3.82)</td><td><strong>0.03</strong></td></tr><tr><td>Anti-Vim60-75</td><td>74.0 (60.0, 84.4)</td><td>57.1 (53.4, 60.3)</td><td>2.14 (1.12, 4.08)</td><td><strong>0.02</strong></td></tr><tr><td>Anti-Vim2-17</td><td>34.0 (20.4, 47.6)</td><td>21.6 (18.7, 24.5)</td><td>1.87 (1.02, 3.43)</td><td><strong>0.04</strong></td></tr><tr><td>Anti-Fibβ62-78 (72)</td><td>28.0 (17.2, 42.1)</td><td>12.8 (10.6, 15.2)</td><td>2.66 (1.39, 5.10)</td><td><strong>0.002</strong></td></tr><tr><td>Anti-CEP1</td><td>70.0 (56.8, 83.2)</td><td>54.0 (50.5, 57.5)</td><td>1.99 (1.07, 3.79)</td><td><strong>0.03</strong></td></tr><tr><td>Anti-Fibα621-635</td><td>58.0 (43.9, 70.9)</td><td>38.9 (35.4, 42.2)</td><td>2.17 (1.21, 3.87)</td><td><strong>0.008</strong></td></tr><tr><td>Anti-Bla26</td><td>48.0 (34.5, 61.8)</td><td>32.4 (29.1, 35.6)</td><td>1.93 (1.09, 3.43)</td><td><strong>0.023</strong></td></tr><tr><td>Anti-Pept5</td><td>68.0 (54.6, 81.4)</td><td>52.7 (49.2, 56.2)</td><td>1.91 (1.04, 3.51)</td><td><strong>0.04</strong></td></tr><tr><td>Anti-CII359-369</td><td>46.0 (32.7, 59.9)</td><td>27.7 (24.6, 30.8)</td><td>2.23 (1.25, 3.96)</td><td><strong>0.006</strong></td></tr><tr><td>Anti-F4-CIT-CIT</td><td>55.1 (40.1, 67.3)</td><td>39.2 (35.7, 42.5)</td><td>1.90 (1.06, 3.40)</td><td><strong>0.028</strong></td></tr><tr><td>Anti-F4-CIT-R</td><td>46.9 (32.7, 59.9)</td><td>23.8 (20.9, 26.8)</td><td>2.83 (1.58, 5.08)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Number of ACPA positive specificities, mean (s.d.)</td><td>10.8 (6.1)</td><td>8.2 (6.5)</td><td>1.06 (1.01, 1.12)</td><td><strong>0.01</strong></td></tr></tbody></table>

P-values shown in bold indicate statistical significance. CI: confidence interval; Collagen type II: CII359-369, F4-CIT-CIT, F4-CIT-R; fibrinogen (Fib): Fibα621-635, Fibβ62-78 (72); filaggrin (Fil307-324); vimentin (Vim)60-75; and mutated proteins (Bla26). DAS28-AUC<sub>24</sub>: Area under curve for accumulated data for DAS28 during 24 months; MAF: minor allele frequency; OR: odds ratio; PF: pulmonary fibrosis.

indicated ILD outcome [7]. These results are in line with our findings as we have focused on PF separately and not the full pattern of ILD, and we found no association with anti-CCP2 antibodies and PF. The antibodies associated with PF in the present study were not among those of the highest frequency in the patient cohort. Three of the antibodies were directed against collagen type II, which is the main structural protein of hyaline cartilage. Another two were respectively directed against the α and β chain of fibrinogen. Presence of the antibodies associated with PF was unrelated to the risk genes and could represent a separate process ongoing in the lung parenchyma or just representing a para-phenomenon.

The strength of the current study is that the cohort of eRA patients was unselected and originated from a

homogeneous population of northern Sweden. Almost all (95%) individuals diagnosed with early RA within the catchment area of northern Sweden were willing to participate in the study. Furthermore, X-rays of the lungs of the patients were routinely performed at inclusion, providing baseline information about the lungs.

We have, however, identified some limitations of this study. The HRCT examinations in this study were not performed randomly or on all included patients, but for those with abnormalities on the plain X-rays or with clinical indications of development of defined symptoms. Further, the HRCT examinations were performed over a period of almost 20 years, and hence methodological improvements during this time could affect the results. Consequently, we refrained from further diagnostic evaluations of ILD besides PF.
<table><caption>TABLE 2 Association between RA with and without PF and ACPA specificities analysed in seven different models</caption><thead><tr><th rowspan="3">Variable</th><th colspan="6">Multiple variable regression analyses<sup>a</sup></th></tr><tr><th colspan="2">+rs2609255</th><th colspan="2">+rs111521887</th><th colspan="2">+rs35705950</th></tr><tr><th>OR (95% CI)</th><th>P-value</th><th>OR (95% CI)</th><th>P-value</th><th>OR (95% CI)</th><th>P-value</th></tr></thead><tbody><tr><td>Model 1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-Vim60-75</td><td>2.26 (1.09, 5.07)</td><td><strong>0.036</strong></td><td>2.28 (1.10, 5.14)</td><td><strong>0.035</strong></td><td>2.44 (1.06, 6.25)</td><td><strong>0.046</strong></td></tr><tr><td>SNP</td><td>1.74 (1.09, 2.74)</td><td><strong>0.018</strong></td><td>1.69 (1.06, 2.564)</td><td><strong>0.024</strong></td><td>2.08 (1.12, 3.70)</td><td><strong>0.015</strong></td></tr><tr><td>Model 2</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-Fibβ62-78 (72)</td><td>2.85 (1.40, 5.59)</td><td><strong>0.003</strong></td><td>2.72 (1.32, 5.33)</td><td><strong>0.005</strong></td><td>3.48 (1.58, 7.39)</td><td><strong>0.001</strong></td></tr><tr><td>SNP</td><td>1.74 (1.09, 2.75)</td><td><strong>0.019</strong></td><td>1.59 (1.01, 2.48)</td><td><strong>0.041</strong></td><td>2.25 (1.20, 4.07)</td><td><strong>0.008</strong></td></tr><tr><td>Model 3</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-Fibα621-635</td><td>1.92 (1.01, 3.72)</td><td><strong>0.047</strong></td><td>1.90 (1.00, 3.68)</td><td>0.052</td><td>1.95 (0.95, 4.11)</td><td>0.071</td></tr><tr><td>SNP</td><td>1.71 (1.07, 2.69)</td><td><strong>0.023</strong></td><td>1.63 (1.03, 2.55)</td><td><strong>0.034</strong></td><td>2.01 (1.08, 3.57)</td><td><strong>0.021</strong></td></tr><tr><td>Model 4</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-F4-CIT-R</td><td>2.46 (1.28, 4.74)</td><td><strong>0.007</strong></td><td>2.52 (1.31, 4.84)</td><td><strong>0.005</strong></td><td>2.34 (1.13, 4.85)</td><td><strong>0.021</strong></td></tr><tr><td>SNP</td><td>1.51 (0.93, 2.42)</td><td>0.087</td><td>1.53 (0.95, 2.41)</td><td><strong>0.074</strong></td><td>1.92 (1.02, 3.44)</td><td><strong>0.033</strong></td></tr><tr><td>Model 5</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-Bla26</td><td>1.89 (1.00, 3.57)</td><td><strong>0.049</strong></td><td>1.90 (1.01, 3.59)</td><td><strong>0.046</strong></td><td>1.80 (0.89, 3.67)</td><td><strong>0.10</strong></td></tr><tr><td>SNP</td><td>1.72 (1.07, 2.72)</td><td><strong>0.023</strong></td><td>1.65 (1.04, 2.58)</td><td><strong>0.029</strong></td><td>2.02 (1.10, 3.60)</td><td><strong>0.019</strong></td></tr><tr><td>Model 6</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Anti-CII359-369</td><td>2.05 (1.07, 3.94)</td><td><strong>0.030</strong></td><td>2.25 (1.17, 4.32)</td><td><strong>0.014</strong></td><td>2.10 (1.02, 4.32)</td><td><strong>0.042</strong></td></tr><tr><td>SNP</td><td>1.66 (1.04, 2.63)</td><td><strong>0.031</strong></td><td>1.72 (1.08, 2.69)</td><td><strong>0.019</strong></td><td>2.06 (1.12, 3.65)</td><td><strong>0.016</strong></td></tr><tr><td>Model 7</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Number of ACPA specificities</td><td>1.07 (1.01, 1.13)</td><td><strong>0.025</strong></td><td>1.07 (1.01, 1.13)</td><td><strong>0.021</strong></td><td>1.07 (1.01, 1.15)</td><td><strong>0.025</strong></td></tr><tr><td>SNP</td><td>1.60 (0.99, 2.55)</td><td>0.050</td><td>1.60 (0.99, 2.53)</td><td><strong>0.047</strong></td><td>1.95 (1.04, 3.50)</td><td><strong>0.029</strong></td></tr></tbody></table>

In this study we have shown that development of PF in patients with RA is associated to the number of ACPA specificities and to certain ACPA fine specificities independently of previously identified PF-related genetic risk loci.

## Acknowledgements

We would like to acknowledge the staff at the departments of rheumatology at the hospitals at Sunderbyn, Sundsvall Örnsköldsvik and Östersund and in particular Mrs Viktoria von Zweigbergk, Umeå for skilful management of the project. All authors were involved in drafting the article or revising it critically for important intellectual content, and all authors approved the final version to be published. Study conception and design: M.B., S.R.-D. Acquisition of data: M.B., S.R.-D., L.L., M.H., J.R., R.H., K.S., G.S., L.K.

**Funding:** This work was supported by grants from the Swedish Research Council (K2013-52X-20307-07-3 and Dnr : 2018-02551); King Gustaf V's 80-Year Fund; the Swedish Rheumatism Association, Region of Västerbotten; and Umeå University.

**Disclosure statement:** The authors have declared no conflicts of interest.

## Data availability statement

The data underlying this article cannot be shared publicly due to the privacy of individuals that participated in the study. The data will be shared on reasonable request to the corresponding author.

## Supplementary data

Supplementary data are available at *Rheumatology* online.

## References

1. Turesson C, O'Fallon WM, Crowson CS, Gabriel SE, Matteson EL. Extra-articular disease manifestations in rheumatoid arthritis: incidence trends and risk factors over 46 years. Ann Rheum Dis 2003;62:722-7.

2. Kelly CA, Saravanan V, Nisar M *et al.*; British Rheumatoid Interstitial Lung (BRILL) Network. Rheumatoid arthritis-related interstitial lung disease: associations, prognostic factors and physiological and radiological characteristics—a large multicentre UK study. Rheumatology 2014; 53:1676-82.

<footnote><sup>a</sup>Adjusted for; sex, age at onset, ever smoker status, methotrexate and RF status at index date. P-values shown in bold indicate statistical significance. CI: confidence interval; Collagen type II: F4-Cit-R; fibrinogen (Fib): Fibα621-635, Fibβ62-78 (72); vimentin (Vim)60-75. OR: odds ratio.</footnote>
3 Makrygiannakis D, Hermansson M, Ulfgren AK *et al.* Smoking increases peptidylarginine deiminase 2 enzyme expression in human lungs and increases citrullination in BAL cells. *Ann Rheum Dis* 2008;67:1488–92.

4 Klareskog L, Stolt P, Lundberg K *et al.* A new model for an etiology of rheumatoid arthritis: smoking may trigger HLA-DR (shared epitope)-restricted immune reactions to autoantigens modified by citrullination. *Arthritis Rheum* 2006;54:38–46.

5 Baka Z, György B, Géher P *et al.* Citrullination under physiological and pathological conditions. *Joint Bone Spine* 2012;79:431–6.

6 Bongartz T, Cantaert T, Atkins SR *et al.* Citrullination in extra-articular manifestations of rheumatoid arthritis. *Rheumatology* 2007;46:70–5.

7 Giles JT, Danoff SK, Sokolove J *et al.* Association of fine specificity and repertoire expansion of anticitrullinated peptide antibodies with rheumatoid arthritis associated interstitial lung disease. *Ann Rheum Dis* 2014;73: 1487–94.

8 Steele MP, Speer MC, Loyd JE *et al.* Clinical and pathologic features of familial interstitial pneumonia. *Am J Respir Crit Care Med* 2005;172:1146–52.

9 Fingerlin TE, Murphy E, Zhang W *et al.* Genome-wide association study identifies multiple susceptibility loci for pulmonary fibrosis. *Nat Genet* 2013;45:613–20.

10 Noth I, Zhang Y, Ma S-F *et al.* Genetic variants associated with idiopathic pulmonary fibrosis susceptibility and mortality: a genome-wide association study. *Lancet Respir Med* 2013;1:309–17.

11 Juge P-A, Lee JS, Ebstein E *et al.* MUC5B promoter variant and rheumatoid arthritis with interstitial lung disease. *N Engl J Med* 2018;379: 2209–19.

12 Jönsson E, Ljung L, Norrman E *et al.* Pulmonary fibrosis in relation to genetic loci in an inception cohort of patients with early rheumatoid arthritis from northern Sweden. *Rheumatology* (Oxford) 2021;61:943–52.

13 Boman A, Brink M, Lundquist A *et al.* Antibodies against citrullinated peptides are associated with clinical and radiological outcomes in patients with early rheumatoid arthritis: a prospective longitudinal inception cohort study. *RMD Open* 2019;5:e000946.

14 Arnett FC, Edworthy SM, Bloch DA *et al.* The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. *Arthritis Rheum* 1988;31:315–24.

15 Raghu G, Remy-Jardin M, Myers JL *et al.*; American Thoracic Society, European Respiratory Society, Japanese Respiratory Society, and Latin American Thoracic Society. Diagnosis of idiopathic pulmonary fibrosis. An official ATS/ERS/JRS/ALAT clinical practice guideline. *Am J Respir Crit Care Med* 2018;198:e44–68.

16 Burkhardt H, Sehnert B, Bockermann R *et al.* Humoral immune response to citrullinated collagen type II determinants in early rheumatoid arthritis. *Eur J Immunol* 2005;35:1643–52.

17 Haag S, Schneider N, Mason DE *et al.* Identification of new citrulline-specific autoantibodies, which bind to human arthritic cartilage, by mass spectrometric analysis of citrullinated type II collagen. *Arthritis Rheumatol* 2014;66:1440–9.

18 Hansson M, Mathsson L, Schlederer T *et al.* Validation of a multiplex chip-based assay for the detection of autoantibodies against citrullinated peptides. *Arthritis Res Ther* 2012;14:R201.

19 R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing, 2013.