FormalPara Key Summary Points

Why carry out this study?

Classical acute phase reactants do not always reflect the synovitis in the patient and may be influenced by other factors.

The search for new biomarkers of disease activity has been a growing field of research in recent years.

Extensive research has reported that the innate immune system, especially the involvement of neutrophils, plays a central role in the initiation and perpetuation of rheumatoid arthritis (RA).

The objective of this study is to determine the performance of neutrophil activation markers, including neutrophil extracellular traps (NETs) and calprotectin, as biomarkers of disease activity using clinical and joint ultrasound (US) parameters in patients with established RA.

We also analyse the relationship between NETs and various types of therapies as well as their association with autoimmunity.

What was learned from this study?

While NET formation induced by neutrophils may play a role in RA pathogenesis, our study raises questions about the utility of NET remnants in peripheral circulation as a biomarker for inflammatory activity.

This study strongly supports the usefulness of calprotectin as a biomarker of inflammatory activity in patients with RA.

This study represents an important contribution to the growing body of research exploring the association among NETs, disease activity and autoantibodies in RA.

Introduction

Rheumatoid arthritis (RA) is a chronic, systemic autoimmune disease characterized by synovial inflammation and cartilage and joint destruction [1]. Over the past few decades, the prognosis of this disease has improved significantly thanks to the adoption of a treat-to-target approach, early diagnosis and treatment, and the utilization of more effective antirheumatic drugs [2, 3]. Currently, close monitoring of the inflammatory disease activity in RA is strongly recommended using a therapeutic goal of achieving remission or low disease activity [2, 4]. In addition to clinical parameters, traditional acute phase reactants (APRs) such as the erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP) serum concentration are commonly employed as laboratory markers to assess disease activity in RA. These parameters are also integrated into composite disease activity indices like the Simplified Disease Activity Index (SDAI) or the Disease Activity Score 28 (DAS28) for longitudinal evaluation [5]. However, while commonly used, APRs may not always reflect the degree of active synovitis in patients with RA accurately.

Extensive research has demonstrated the significant involvement of adaptive immunity in the pathogenesis of RA. Systemic autoimmunity is found in most patients with RA, with the presence of characteristic autoantibodies such as rheumatoid factor (RF) or anti-citrullinated peptide/protein antibodies (ACPAs). However, mounting evidence indicates that the innate immune system, especially the involvement of neutrophils, plays a central role in the initiation and perpetuation of RA, both directly through its effects on the synovium and indirectly via its own inflammatory response and by modulating that of the adaptive immune system [6,7,8]. In recent years, increasing attention has been paid to the role of neutrophils and their products in RA. Calprotectin (S100A9/S100A8), a protein member of the S100 family which is secreted by activate neutrophils, seems to play a significant role in the inflammatory cascade in RA. Calprotectin is a significant proinflammatory factor of innate immunity and acts as an endogenous damage-associated molecular pattern (DAMP) through the activation of TLR4 [9]. Moreover, when secreted into the synovium, calprotectin subsequently enters peripheral blood circulation. Levels of calprotectin in serum or plasma have been found to correlate closely with inflammatory disease activity in patients with RA and therefore provide a more accurate reflection of the patient's clinical status than APRs [10,11,12,13].

In 2004, a mechanism of action was reported for neutrophils through the formation of neutrophil extracellular traps (NETs) [14]. NETs are extracellular structures that contain neutrophil DNA, histones and proteins and which are ejected into the extracellular medium in the form of a network. NETs provide an alternative defence mechanism and may play a crucial role in regulating the immune response and maintaining bodily homeostasis. They are formed through a process called NETosis, which can be triggered by various stimuli, including both infectious and sterile factors such as autoantibodies or immunocomplexes [15].

The role of NETs in RA has garnered significant interest in recent years. These inflammatory structures are believed to have both local and systemic effects, potentially contributing to the initiation and perpetuation of the disease through autoantigen exposure [16,17,18,19,20,21]. Furthermore, NETs have been identified as a source of specific autoantigens in RA, including ACPAs and anti-carbamylated protein (anti-CarP) antibodies [22,23,24,25,26]. Patients with RA demonstrate elevated NET formation compared to the general population, revealing NETs to be a promising potential biomarker for the disease [27]. However, the precise role of NETs in disease activity remains uncertain as studies evaluating their utility as a biomarker have produced conflicting results [17, 25, 26, 28].

The objective of this study is to determine the performance of neutrophil activation markers, including NETs and calprotectin, as biomarkers of disease activity using clinical and joint ultrasound (US) parameters in patients with established RA. The study also aims to examine the relationship between NETs and various types of therapies as well as their association with serum autoantibodies.

Methods

Design and Study Population

We performed an observational cross-sectional study. Patients with RA (ACR/EULAR 2010 criteria) [29] from our Arthritis Unit who had been undergoing treatment with biological disease-modifying antirheumatic drugs (DMARDs), tumour necrosis factor inhibitors (anti-TNF), monoclonal antibodies against IL-6 receptors (anti-IL6r) or Janus kinase inhibitors (JAKi) for at least 3 months were consecutively included, irrespective of their disease activity status, previous use of DMARDs (including biological therapies or JAKi) or concomitant treatment (such as with methotrexate or others). Patients who displayed signs of active infection or other clinical conditions that could impact the determination of neutrophil activation markers (calprotectin or NETs) or APR measurements, as determined by the researcher concerned, were excluded. Demographic data, disease duration, autoantibody status (RF, ACPAs and anti-CarP), radiological data (presence/absence of erosive disease) and information on previous biological therapy and concurrent treatments were collected.

Two control groups were included to assess NET levels. The first consisted of patients with RA [29] from our Arthritis Unit with high disease activity, defined as those with at least four inflamed joints, regardless of treatment or disease duration. The second control group was comprised of healthy subjects.

Measurement of Clinical Disease Activity

All patients received a clinical evaluation that included the counting of 28 swollen and tender joints (28SJC and 28TJC) as well as physician and patient self-reported global assessments (PhGA and PGA respectively) using visual analogue scales ranging from 0 to 10 along a centimetre scale. Subsequently, composite disease activity indices such as DAS28, SDAI and CDAI (Clinical Disease Activity Index) were calculated. Additionally, patients were asked to complete two questionnaires: the Health Assessment Questionnaire (HAQ) and Routine Assessment of Patient Index Data 3 (RAPID3).

Assessment of Biomarkers

Acute Phase Reactants

Blood samples were collected concurrently with the clinical assessment. ESR was measured using the Westergren method (normal value [NV] < 20 mm/h) and high-sensitivity C-reactive protein (hsCRP) plasma concentration using an immunoturbidimetric method measured with Siemens Atellica® Solution (lowest limit of detection: 0.02 mg/dl; NV < 0.4 mg/dl).

Plasma Calprotectin

Calprotectin was measured using an enzyme-linked immunosorbent assay (ELISA) test kit (CALPROLAB ALP [CALPRO], Menarini Diagnósticos S.A.) following the manufacturer's instructions. Briefly, 100 µl aliquots of each standard, control and diluted 1:20 sample in duplicate wells were incubated at room temperature for 40 min; three washings were performed, 100 µl of the conjugated enzyme was added and the plates were incubated at room temperature for another 40 min. After three further washes and the addition of the enzyme substrate, optical density values at 405 nm were determined using an ELISA reader. To reduce variation in calprotectin determinations, the whole procedure was performed in a Triturus Autoanalyzer; the coefficients of variation were 5% within and 13% between assays. This technique had already been used by our group for previous studies[11].

Assessment of NETs

NETs were indirectly determined by measuring their remnants in plasma. Neutrophil elastase-DNA (NE-DNA) and histone-DNA (H3-DNA) complex levels were determined using a tailor-made ELISA[30, 31]. A 96-well plate was coated with 100 µl of capture antibodies, anti-elastase (Millipore cat. no. 481001) and anti-histone H3 citrulline (Abcam cat. no. ab5103) at 2.5 µg/ml in sterile PBS. The ELISA plates were incubated overnight at 4 °C without agitation to allow the binding of capture antibodies. The next day, the wells were washed three times with PBS-Tween 20 (200 μl/well) and blocked with PBS-1% BSA (200 μl/well) at room temperature for 1 h. Again, three repeated washes with PBS-Tween 20 (200 μl/well) were performed. Diluted plasma (1:100, with PBS-1% BSA) was added to the wells (100 μl/well) and again incubated overnight at 4 °C but with agitation (100–150 rpm). Three washes with PBS-Tween 20 (200 μl/well) were performed to remove unbound NETs and other components. Mouse anti-dsDNA (Millipore cat no. MAB030) diluted 1:100 with PBS-1% BSA was added to the wells (100 μl/well) and incubated for 1 h, at room temperature, with agitation (100–150 rpm). The plates were washed three more times with PBS-Tween 20 (200 μl/well). Goat anti-mouse, HRP conjugate (Millipore cat no. AP127P) diluted 1:10,000 was added (100 μl/well) and incubated for 1 h at room temperature. The plates were washed five times with PBS-Tween 20 (200 μl/well). Then, TMB peroxidase substrate was added to each well (50 μl/well) and the plate was incubated in the dark for between 10 and 15 min. The reaction was stopped by the addition of 1 M HCl to each well (50 μl/well). Finally, absorbance was measured at 450 nm using a microplate photometer, and an OD index (OD each well/OD blank) was calculated for each well. The assay was performed in duplicate and mean values were calculated.

Assessment of Autoantibodies (RF, ACPAs and Anti-CarP)

Autoantibody status was measured in sera collected at enrolment in the study.

IgM RF was determined using a chemiluminescence assay (CIA) (QUANTA Flash®, Inova Diagnostics, CA). The assay's analytical measurement range (AMR) is from 0.3 to 490.0 IU/ml and the cut-off recommended by the manufacturer is < 5 IU/ml.

ACPA levels were determined by a tailor-made ELISA using a chimeric fibrin/filaggrin citrullinated peptide (anti-CFFCP1) as an antigen [32] biotinylated at its terminal amino group with a biotin derivative containing two polyethylene glycol (PEG) chains bound to neutravidin derivatized microtitre plates. This derivatization strategy allows for better antigenic exposure of the peptide in the ELISA plates.

Anti-CarP antibodies were also determined using a tailor-made ELISA test using carbamylated foetal calf serum (FCS-CarP) as the antigen and the non-homocitrullinated version of the protein as the control for the homocitrulline anti-FCS-CarP detected [33].

A positive cut-off value for the tests was defined as ≥ 11.5 AU/ml and ≥ 173.5 AU/ml for ACPAs (anti-CFFCP1) and anti-CarP antibodies (anti-FCS-CarP), respectively. A result was only considered positive and specific for citrulline and/or homocitrulline when the UA/ml values were higher than the respective cut-offs and the difference in OD values between native and post-translationally modified antigens was > 0.1. Further details of the ELISA techniques are provided in the Supplementary Material.

Imaging Biomarkers: Ultrasound Score

Sonographic assessments were carried out using high-sensitivity US equipment (MyLab9®; Esaote, Genoa, Italy), a longitudinal probe with a 10–14 MHz frequency range and a pulse repetition frequency of between 800 and 900 Hz. The joint musculoskeletal US findings were defined based on the Outcome Measures in Rheumatoid Arthritis Clinical Trials (OMERACT) definitions [34].

In accordance with the EULAR guidelines [35], synovial hypertrophy (SH) and intra-articular power Doppler (PD) signalling of 11 joints and tendons on each hand, including the proximal interphalangeal joints, metacarpophalangeal joints and wrists, were evaluated by a single experienced sonographer (A.P.). The sonographer was blinded to the results of the clinical joint examination. SH and PD signals were graded using a four-grade semi-quantitative scoring system (0 = no, 1 = mild, 2 = moderate and 3 = severe) according to the methodology of Szkudlarek et al. [36]. The highest grade of SH and PD signalling detected during the scans was selected as representative of each joint, respectively.

The PD score (sum of PD scores in all joints, range: 0–66), SH score (sum of SH scores in all joints, range: 0–66) and global score (sum of the PD and SH scores, range: 0–132) were calculated by adding up the scores for elementary lesions in each joint. To ensure strict criteria for defining US synovitis, only patients with SH grade ≥ 2 plus PD signal ≥ 1 were considered to have active synovitis [37]. This US assessment method had been previously employed in prior studies conducted by our group [11]. Clinically symptomatic joints were also evaluated to assess whether US synovitis was present. However, these joints were not used to calculate the sonographic score (PD score, SH score and global US score).

Statistical Analysis

Continuous data were presented as mean and standard deviation (SD) while categorical variables were given as absolute frequency with percentages. Groups were compared using parametric or nonparametric tests according to the distribution of the variables.

Correlation analysis (Spearman’s correlation coefficient) was used to assess the association between the different activity indices (clinical and US activity) and the neutrophil activation markers (NET and calprotectin plasma levels). To specifically investigate the relationship between calprotectin and the US parameters (including the US total score, SH score, and PD score), we utilized a linear regression model. This regression analysis was adjusted to account for potential confounding factors, with the selection of these variables guided solely by clinical criteria. Furthermore, similar models were constructed to assess the association among ESR, hsCRP and the US parameters. We also analysed the correlation between autoantibody levels and NETs. The performance of NETs (NE-DNA complexes and H3-DNA complexes) and calprotectin in the diagnosis of US synovitis was analysed using receiver-operating characteristic curves (ROC) with US synovitis yes/no (yes: SH grade ⩾ 2 plus PD signal ⩾ 1) as the gold standard. The ROC curves made it possible to identify the best cut-off point in terms of sensitivity, while specificity enabled us to calculate the area under the curve (AUC) as a measure of the overall discriminative power. The performance of ESR and hsCRP was also studied, and the three ROC curves were compared to determine which of the three parameters had the highest discriminatory power for the diagnosis of US synovitis. The test for the equality of the AUC uses an algorithm suggested by DeLong, DeLong and Clarke-Pearson (1988) [38]. Additionally, a logistic regression model was conducted to investigate the correlation between calprotectin and the presence of US synovitis (categorized as "yes" when SH grade ⩾ 2 plus PD signal ⩾ 1 were present). Similar to the linear regression model, these logistic regression models were adjusted for the same confounding variables.

Analysis was carried out using SPSS software (version 27.0, IBM Inc., Chicago, IL, USA) and Stata version 13.1 (StataCorp LP 4905 Lakeway Drive College Station, TX, USA).

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Clinical Research Ethics Committee of the Hospital Clinic of Barcelona (reg. HCB20210783). Written informed consent was obtained from all patients before enrolment in the study.

Results

Demographic, Clinical and Therapeutic Characteristics

The study included 114 patients: 90.4% women; mean age 56 ± 11 years; mean disease duration 16 ± 20 years. Most patients (86%) were seropositive for either RF or ACPAs, 67.5% had erosive disease and 31.6% had extra-articular manifestations of RA. Fifty-six patients (49%) were receiving treatment with anti-IL6r (49 tocilizumab and 7 sarilumab), 28 (25%) were receiving JAKi (19 baricitinib and 9 tofacitinib) and 30 (26%) were receiving anti-TNF treatment (21 etanercept and 9 adalimumab). The mean duration of these treatments was 56 ± 58 months. Nearly half of the patients were receiving concomitant conventional DMARD treatment (methotrexate or leflunomide) (49.1%), and 44.7% were also receiving glucocorticoid treatment (mean prednisone equivalent dose 3.3 ± 2.8 mg/day) (Table 1).

Table 1 Demographic characteristics and treatment of patients with RA

As for disease activity, the mean CDAI was 12 ± 10; 64.3% of the patients were in remission and/or exhibited low disease activity according to CDAI criteria (Table 2). Regarding laboratory parameters, the mean hsCRP was 0.4 ± 1.3 mg/dl, the ESR was 14 ± 18 mm/h and the mean plasma calprotectin concentration was 0.99 ± 1.5 μg/ml. When comparing the parameters between the groups, we found that the classical APRs (both hsCRP and ESR) were significantly lower in the group of patients treated with anti-IL6r. However, we did not find statistically significant differences in plasma calprotectin levels between the groups (Table 2).

Table 2 Clinical, ultrasound and laboratory characteristics of patients with RA

Ninety-nine patients were studied with US, of which 66 (67%) were active according to US criteria (SH score ≥ 2 and PD score ≥ 1). The mean scores for SH, PD and total US activity were 6.6 ± 6.4, 4.1 ± 5.5 and 11 ± 12, respectively.

When differentiating by treatment group, we found no differences between the patients in terms of the different clinical and US parameters of disease activity, except for DAS28, which was lower in the group of patients treated with anti-IL6r (Table 2).

Two control groups were included to compare NET levels. The active control group consisted of 15 patients: 93% female; mean age 54 ± 12 years. The mean disease duration was 13 ± 15 years. All these patients were seropositive, 67% had erosive disease, and 13% had extra-articular manifestations of RA. Eleven patients (73%) were not receiving targeted therapy, neither biological DMARDs nor JAKi. Four (27%) were receiving treatment with anti-IL6r, and three (20%) with rituximab. At the time of inclusion, 73% were receiving concomitant glucocorticoid therapy (mean prednisone equivalent dose 7.0 ± 7.8 mg/day). The second control group consisted of 30 healthy subjects (mean age 47 ± 10 years, 73% women).

NET Plasma Levels in RA According to the Antirheumatic Drug Group

We analysed the levels of circulating NET remnants, measured as DNA-NE complexes and DNA-H3 complexes. Plasma levels of neither NET remnant showed statistically significant differences between the treatment groups (anti-Il6r, JAKi and anti-TNF) (Table 3).

Table 3 NET remnant levels in the main group vs control groups (active patients and healthy subjects)

Since the assessment of NETs is a novel and non-standardized technique, we included a group of healthy subjects and a group of patients with highly active RA, as previously defined, as negative and positive controls, respectively. Even though the active patient group exhibited higher levels of NETs than the main group and the healthy subject group, this difference did not reach statistical significance (Table 3).

Association Between NETs and Plasma Calprotectin with Clinical and Ultrasonographic Disease Activity

We analysed the correlation between neutrophil activation markers (calprotectin and NET levels) and clinical and US activity parameters as well as APRs.

We found that NET levels did not show a correlation with any of the clinical and US activity parameters analysed, either as a single group (Table 4) or when considering the different treatment groups separately (Table 5).

Table 4 Correlation of NETs and plasma calprotectin with clinical and ultrasonographic disease activity
Table 5 Correlation of NETs and plasma calprotectin with clinical and ultrasonographic disease activity by treatment group

Plasma calprotectin showed a moderate correlation with clinical parameters (SJC rho = 0.49, PhGA rho = 0.39, CDAI rho = 0.30; P < 0.001) and US parameters (rho > 0.50; P < 0.001). This correlation was stronger than that observed with APR (Table 4). To further examine the relationship between plasma calprotectin levels and US parameters, we conducted a linear regression model. This analysis was adjusted to account for potential confounding factors such as gender, smoking status, disease duration, use of glucocorticosteroids and drug therapy.

Our findings revealed a positive association between plasma calprotectin and SH and PD scores. Specifically, each increase in plasma calprotectin units corresponded to an increase of 2.14 units in HS score (95% confidence interval: 0.71–3.58; P value = 0.004) and 1.66 units in PD score (95% confidence interval: 0.35–2.7; P value = 0.013). However, no significant association was observed between plasma calprotectin levels and the total US score. Similar models were applied to ESR and hsCRP but failed to demonstrate a significant association (refer to Supplementary Material for details).

When we analysed the correlations by treatment groups, we found that calprotectin showed a moderate or strong correlation with US parameters, unlike classical APR; this is especially relevant for the group of patients treated with IL6 inhibitors (Table 5).

When analysing the discriminatory capacity by the AUC of both the neutrophil biomarkers of the presence of active US synovitis, calprotectin demonstrated an excellent capacity (Fig. 1) but plasma NETs did not (Figs. 2 and 3). This discriminatory capacity of plasma calprotectin was reflected by an AUC of 0.812 (95% confidence interval [CI]: 0.729–0.894), in contrast to the lower discriminatory capacity of hsCRP and ESR, with AUCs of 0.598 (0.531–0.665) and 0.560 (0.444–0.676), respectively (Fig. 1a). When we analysed the treatment groups separately, we found that these differences were maintained (Fig. 1b, c and e).

Fig. 1
figure 1

ROC curves of plasma calprotectin, hsCRP and ESR vs ultrasound synovitis according to type of treatment. a Main group (n = 114); b anti-Il6r (n = 56); c JAKi (n = 28); d anti-TNF (n = 30). Anti-Il6r monoclonal antibodies against IL-6 receptors; JAKi JAK inhibitors; anti-TNF tumour necrosis factor inhibitors; hsCRP high sensitivity C-reactive protein; ESR erythrocyte sedimentation rate; ROC receiver-operating characteristic

Fig. 2
figure 2

ROC curves of NE-DNA complexes, hsCRP and ESR vs ultrasound synovitis according to type of treatment. a Main group (n = 114); b anti-Il6r (n = 56); c JAKi (n = 28); d anti-TNF (n = 30). Anti-Il6r monoclonal antibodies against interleukin-6 receptors; JAKi JAK inhibitors; anti-TNF tumour necrosis factor inhibitors; NE-DNA complex neutrophil elastase-DNA complexes (NETs); hsCRP high sensitivity C-reactive protein; ESR erythrocyte sedimentation rate; ROC receiver-operating characteristic

Fig. 3
figure 3

ROC curves of H3-DNA complexes, hsCRP and ESR vs ultrasound synovitis according to type of treatment. a Main group (n = 114); b anti-Il6r (n = 56); c JAKi (n = 28); d anti-TNF (n = 30). Anti-Il6r monoclonal antibodies against interleukin-6 receptors; JAKi Janus kinase inhibitors; anti-TNF tumour necrosis factor inhibitors; H3-DNA complex histone-DNA complexes (NETs); hsCRP high-sensitivity C-reactive protein; ESR erythrocyte sedimentation rate; ROC receiver-operating characteristic

Furthermore, a logistic regression model was employed to examine the association between calprotectin and the presence of US synovitis (categorized as "yes" when SH grade ⩾ 2 plus PD signal ⩾ 1 was present). Using a calprotectin level of 0.4 mcg/dl as the cut-off point, our analysis indicated that patients with plasma calprotectin values > 0.4 mcg/dl faced an 11-fold higher risk of being identified as active via ultrasound compared to those with lower values (OR 11.31; 95% CI 3.18–40.17; P < 0.001) (refer to Supplementary Material for details).

Correlation Results Between NETs and Autoantibody Levels

We analysed potential correlations between NET levels and autoantibodies (ACPAs, anti-CarP antibodies and RF). However, our results showed no correlations, even when autoantibody status was categorized into different quartiles (data not shown). We also did not find any differences in plasma NET levels when we classified patients as seropositive or seronegative for each of the autoantibodies analysed (Table 6).

Table 6 NET titres in patients according to autoantibody positivity

Correlation Results Between Calprotectin and ACPA Status

We have analysed whether the correlation of calprotectin with clinical and ultrasound parameters of disease activity differs based on ACPA status (ACPA positive vs ACPA negative).

We have observed that patients who are ACPA positive maintain a strong correlation between plasma calprotectin levels and indices of clinical disease activity. However, in the ACPA-negative patient group, this correlation diminishes. On the other hand, when analysing the correlation with US parameters, we did find any differences between the two subgroups (Table 7).

Table 7 Correlation of plasma calprotectin with clinical and ultrasonographic disease activity, by ACPA status

Discussion

One of the main objectives of the present study was to assess the utility of plasma NET remnants as a parameter for evaluating clinical disease activity in patients with established RA. However, our findings demonstrated that plasma NET levels did not reflect the inflammatory disease status in our patients, regardless of the specific targeted therapy group analysed. Furthermore, we investigated the potential relationship between NET remnants and US-detected synovitis, but no significant association was observed. These results are consistent with previous studies that also failed to find a significant association between peripheral blood levels of NETs and disease activity [17, 39, 40]. However, some studies have shown an association with clinical disease activity, but those studies are heterogeneous, included limited numbers of patients and the associations described are weak and only have a few of the activity parameters analysed [21, 25, 26, 28, 41, 42].

The lack of correlation between NETs and inflammatory disease activity in our cohort may be attributed to the inclusion of patients undergoing biological or JAKi treatment with low disease activity. Nevertheless, even in our control group consisting of patients with high disease activity, NET levels were only slightly elevated and there were no significant differences compared to the main group. These findings indicate very strongly that NET remnants cannot be considered reliable biomarkers of disease activity in patients with RA, particularly those receiving targeted treatment. It is important to note, however, that the observed low levels of NETs in our patients may be influenced by the effect of biologics or JAKi, as previous research has reported reductions in NET formation following different treatments with biological drugs [28, 41].

The involvement of NETs in the pathogenesis of RA and their potential as a antigenic source and autoantibodies has been suggested [16, 22, 43]. Therefore, in this study, we aimed to determine the association between plasma NETs and the characteristic autoantibodies seen in RA. We simultaneously measured the levels of RF, ACPAs and anti-CarP antibodies in the same blood samples used for analysing NET levels. We did not find any correlation between the titres of autoantibodies and the levels of NETs in peripheral blood. Low levels of NETs together with the effect of biological therapy may affect the possible relationship between antibody levels and NET remnants in our cohort. One study found an association between ACPA titres and serum NET remnants (myeloperoxidase-DNA and NE-DNA complexes) in patients with RA, but only in those patients with extremely high ACPA levels [24]. The discrepancies observed between those findings and our study may be attributed to factors such as geographic origin, disease duration, methodological issues in NET determination or differences in drug therapies among the study populations.

The lack of association of plasma NET remnants and autoantibody levels in our study does not rule out a role for NETs in autoantibody production in RA as a consequence of higher antigenic exposition. Several in vitro studies corroborate the role of NETs in autoantibody production and the idea that autoantibodies stimulate neutrophils to form NETs [16, 22, 43]. In vitro studies have also demonstrated that neutrophils from patients with RA exhibit increased spontaneous NET formation, which correlates with levels of ACPAs [17, 22]. One major factor that predisposes patients to exacerbated NETosis is the proinflammatory environment present in RA-affected tissues. This inflammatory milieu, generated by cells such as synovial fibroblasts, dendritic cells and macrophages releasing proinflammatory cytokines like TNF-α, IL-1β and interferon-gamma (IFN-γ), can activate neutrophils and promote their migration to inflamed joints. Once there, the neutrophils are continuously stimulated by the proinflammatory environment, triggering the release of NETs. Exposure of neutrophils to sera from patients with RA, especially those with elevated ACPA and RF levels, as well as the inflammatory cytokines IL-17A and TNF, induces NET formation in isolated RA neutrophils, which are more susceptible to NET formation than those from healthy individuals [22]. NETs are also an important source of carbamylated proteins in patients with RA [44]. Anti-CarP antibodies are frequently present in RA and are associated with increased joint destruction, mortality and interstitial lung disease [33, 45]. These antibodies can activate osteoclasts and contribute to bone resorption, thus suggesting a causal relationship with bone damage in RA [44]. Moreover, the specific types of NETs involved in RA pathogenesis warrant further investigation. While some studies have identified specific NET components, such as citrullinated proteins, others have highlighted the presence of histones or antimicrobial peptides. Understanding the heterogeneity of NETs and their contribution to disease processes is crucial if we are to gain deeper insight into the complex mechanisms of RA.

Overall, the induction of NETs by autoantibodies, the association between blood NET levels and ACPA titres and the variations observed in different experimental settings together emphasize the intricate nature of the involvement of NETs in RA pathogenesis. Further research is needed to elucidate the specific subtypes and functional roles of NETs in the context of disease progression and treatment response.

Additionally, we examined a well-known neutrophil activity marker: calprotectin. Our results confirm that calprotectin is a good biomarker of disease activity in patients with RA, showing a correlation with activity measured by clinical and US parameters. It is a promising biomarker with potential use in monitoring patients, especially those using treatments that modify the levels of classical APR, such as anti-IL6r [46, 47] or JAKi [46, 47]. These results are in line with the existing literature. Previous studies have indicated that this protein reflects a patient's clinical status more accurately than APRs [10,11,12,13], including clinical remission or low disease activity. Furthermore, calprotectin serves as a prognostic biomarker for radiographic progression [48], a marker of therapeutic response to specific targeted DMARDs [20], and it has proved useful in predicting disease flares in seemingly controlled patients [49]. A recent publication by Sejersen et al. explored the relationship between calprotectin levels and autoantibody-characterized subgroups [50]. They described association between circulating calprotectin and inflammation in patients who are ACPA positive but not in patients with RA who are ACPA negative. In line with this research, our analysis of the cohort, stratified based on autoantibody status, revealed a notable correlation between calprotectin and clinical parameters in patients with RA who are ACPA positive but not in patients who are ACPA negative. However, interestingly, we did not observe any significant differences in ultrasound (US) parameters. It is important to note that the sample size in the seronegative subgroup in our cohort was small.

The present study has some limitations. First, the methodology used to assess NETs relied on an ELISA that detects proteins expressed in NETs but which are not exclusive to these structures. This indirect approach may have introduced variability and potential bias into the results. Another limitation of our study was its cross-sectional design, which restricted our ability to establish causality or observe changes in NET levels over time. Additionally, the patient cohort in our study consisted of individuals with long-standing RA who were undergoing targeted therapies. These treatments may influence the levels of NETs, potentially leading to lower levels than those observed in patients who were naïve to treatment. It is essential to acknowledge the potential confounding effect of medications on NET levels and to consider studying patients at different disease stages to better understand the dynamics of NETs in RA.

Despite these limitations, this study constitutes a significant contribution to the expanding body of research investigating the interplay among NETs, disease activity and autoantibodies in RA. However, based on the results of this study, we cannot advocate for the routine clinical measurement of NETs to monitor disease activity. Conversely, our findings solidify the notion that calprotectin serves as an exceptionally valuable biomarker for assessing inflammatory activity (both clinical and ultrasonographic). It plays a pivotal role, particularly in evaluating patients undergoing treatments that inhibit the production of classical RFA, such as rIL6 antagonists or JAKi.

Conclusion

NET formation induced by neutrophils may play a potential role in RA, but our results question the utility of NET remnants in peripheral circulation as a biomarker of inflammatory activity in this disease. No association between NETS and clinical and US disease activity or with autoantibodies was observed. By contrast, our study clearly confirms the usefulness of blood calprotectin, a protein derived from neutrophil activation, as a biomarker of inflammatory activity and local synovitis in patients with RA, thereby confirming the results of previous studies.