1. Introduction
Rheumatoid arthritis (RA) is a chronic autoimmune inflammatory disease, which may lead to articular bone and cartilage destruction, disability and reduced life expectancy [
1,
2]. RA develops upon aberrant activation of the immune system mainly due to failure of self-tolerance mechanisms [
3]. T helper (Th) lymphocytes have been demonstrated to play a central role in disease pathogenesis and progression through the release of cytokines and chemokines [
4,
5]. Both Th1 and Th17 cells have been proved to contribute to RA inflammatory responses. On the other hand, CD4
+ T regulatory cells (FoxP3
+ Tregs) as well as myeloid-derived suppressor cells (MDSCs) that exert an immunoregulatory role and re-establish homeostasis are defective during active RA [
6]. Restoration of immune regulation has been shown in RA patients upon effective treatment.
T cell proliferation and acquisition of effector functions require intracellular signals elicited by both the T cell receptor a co-receptor, CD28 [
7]. CD28 delivers co-stimulatory signals upon engagement of its ligands, B7-1 (CD80) or B7-2 (CD86). Cytotoxic T lymphocyte antigen 4 (CTLA4), binds to the same ligands as CD28, with higher affinity and restricts T cell activation [
8]. CTLA4 expression is induced in all T cells transiently after T cell receptor activation, to regulate T cell proliferation and maintenance of tolerance [
9]. Abatacept (CTLA4-Ig) is a recombinant fusion protein comprising the extracellular domain of human CTLA4 and a fragment of the Fc domain of human IgG1, which has been modified to prevent complement fixation [
10]. Abatacept, like CTLA4, competes with CD28 for CD80 and CD86 binding on DCs, and thereby can be used to selectively modulate T-cell activation [
11,
12]. Abatacept has been approved for the treatment of RA, based on the results of an extensive clinical development program assessing its effectiveness and safety in different RA populations [
13,
14,
15].
Biologic therapies (bDMARDs) have revolutionized the treatment of RA patients. Nevertheless, data from registries and cohort studies have shown that in clinical practice, 50-60% of RA patients starting a bDMARD will stop treatment due to inefficacy or toxicity in the long-term, while clinical response (remission or low disease activity) is unpredictable and is achieved by 20-40% [
16,
17]. Development of predictive markers of response for everyday clinical practice will be an important step in optimizing the treatment of autoimmune diseases. Available data for predictors of abatacept response are rather limited [
18].
In this study we propose that a detailed immunological profile early at treatment initiation with abatacept, could differentiate responders vs non-responders. Moreover, sera proteomics analysis revealed a signature of proteins associated with clinical responses. Immunological studies of peripheral blood (PB) and sera were performed at baseline, on 3rd and 6th month after treatment initiation further characterized the biological effects of abatacept. Based on the above data we propose a signature which could predict response to abatacept.
2. Materials and Methods
2.1. Study approval
This is an observational, prospective, single-center study of RA patients, starting treatment with abatacept, due to residual disease activity. All patients were recruited by the Rheumatology and Clinical Immunology Clinic at the University Hospital in Heraklion, Crete. Decision regarding treatment scheme was made by the treating rheumatologist and according to the guidelines for RA treatment of the Hellenic Society of Rheumatology and EULAR guidelines. The study was approved by the Ethics Committee of the University Hospital of Heraklion (protocol number:3601, 17/06/2015, Heraklion) and informed consent was obtained from all individuals prior to sample collection.
2.2. Human subjects
Peripheral blood (PB) samples were obtained from RA individuals before the initiation of abatacept treatment and at 3 and 6 months upon treatment. Subjects with RA fulfilled the 2010 EULAR/ACR classification criteria, at Rheumatology and Clinical Immunology Clinic at the University Hospital of Heraklion [
19]. At the time of sampling all patients had active disease state (DAS28≥3.2) and were naïve for biological therapy. Patients were followed clinically every 3 months and all immunological studies were performed at the indicated time points.
2.3. Peripheral blood mononuclear cell isolation
Heparinized blood (20ml) was collected from individuals with RA. Peripheral blood mononuclear cells (PBMCs) were isolated on Histopaque-1077 (Sigma-Aldrich) density gradient. Briefly, blood was diluted 1:1 with PBS and layered over Histopaque medium. Tubes were centrifuged at 400g for 30 min at room temperature (RT) to obtain PBMC layer. PBMC layer was collected and cells were washed with PBS. Erythrocytes were eliminated by hypotonic lysis (1ml ddH2O for 35 seconds and 1ml 1.8% NaCl). Then, samples were processed for staining. Additionally, blood in appropriate serum separator tubes was collected from individuals with RA. Tubes were centrifuged at 2500g for 15 min at RT to obtain serum.
2.4. Antibodies
Fluorescence-conjugated antibodies against human with the following specificities were used to detect antigens by flow cytometry: CD4 (OKT4), IL-17 (BL168), CD127 (A019D5), CD25 (BC96), HLADR (L243), CD3 (OKT3), CD1α (HI149), CD14 (HCD14), CD15 (SSEA-1), CD33 (WM53) from Biolegend; IFN-γ (B27) from BD Pharmingen; and FoxP3 (PCH101) from eBioscience.
2.5. Flow cytometry
Single-cell suspensions were stained with fluorescence-conjugated antibodies and analyzed by flow cytometry. In brief, 3–5 × 106 PBMCs were washed with FACS buffer (5% FBS in PBS), and cells were then stained with the indicated antibodies for 20 min at 4°C. Stained cells were washed twice with FACS buffer, resuspended and acquired on a FACS Calibur (BD Biosciences) and analyzed using the FlowJo software (Tree Star).
2.6. Cellular Cytokine Expression Assays
Cytokine production was assessed through intracellular staining. In brief, to study Th cell populations (Th1, Th17), PBMCs were resuspended into a concentration of 3 × 106 cells/mL in RPMI 1640 supplemented with 10% FBS, 100IU/ml penicillin, and 100μg/ml streptomycin cell culture media in 96-well flat-bottom plates and stimulated with PMA (50 ng/ml) and ionomycin (2 µg/ml) for 6 hours in the presence of Brefeldin A (10 µg/ml) in a 5% CO2 atmosphere in a 37oC incubator. Stimulation was stopped by washing the cells with cold PBS. Then, cells resuspended in FACS buffer and stained with anti-CD4. For the intracellular staining, CD4-stained cells were fixed with formaldehyde solution (4%) and permeabilized using saponin (10% w/v) and stained with conjugated antibodies against human IFN-γ and IL-17. For FoxP3 intracellular staining, CD4-stained cells were permeabilized using the FoxP3 Staining Buffer Set (eBioscience), according to the manufacturer’s instructions. Finally, the cells were washed with FACS buffer before their flow cytometric analysis.
2.7. Sample preparation and Proteomics analysis
The protein concentration of plasma samples was measured by the Bradford assay. An appropriate volume containing 200 μg total protein per sample was processed with the filter aided sample preparation (FASP) protocol [
20]. Briefly, proteins were reduced with DTE (0.1 M), alkylated with iodoacetamide (0.05 M), and digested overnight by trypsin in 50 mM NH4HCO3 pH 8.5. The peptides originating from tryptic digestion were lyophilized and kept at −80 °C. The peptides were further cleaned-up using the Sp3 protocol for peptides. Samples were analyzed on a liquid chromatography tandem mass spectrometry (LC-MS/MS) setup consisting of a Dionex Ultimate 3000 nanoRSLC coupled in line with a Thermo Q Exactive HF-X Orbitrap mass spectrometer. Peptidic samples were directly injected and separated on an 25 cm-long analytical C18 column (PepSep, 1.9μm3 beads, 75 µm ID) using an one-hour long run, starting with a gradient of 7% Buffer B (0.1% Formic acid in 80% Acetonitrile) to 35% for 40 min and followed by an increase to 45% in 5 min and a second increase to 99% in 0.5 min and then kept constant for equilibration for 14.5min. A full MS was acquired in profile mode using a Q Exactive HF-X Hybrid Quadrupole-Orbitrap mass spectrometer, operating in the scan range of 375- 1400 m/z using 120K resolving power with an AGC of 3x 106 and maximum IT of 60 ms followed by data independent acquisition method using 8 windows (a total of 39 loop counts) each with 15K resolving power with an AGC of 3x 105 and max IT of 22ms and normalized collision energy (NCE) of 26.
Orbitrap raw data were analyzed in DIA-NN 1.8.1 (Data-Independent Acquisition by Neural Networks) through searching against the Human Proteome (downloaded from Uniprot 20583 proteins entries, 8/11/2022) using the library free mode of the software, allowing up to two tryptic missed cleavages and a maximum of three variable modifications/peptide. A spectral library was created from the DIA runs and used to re-analyze them (double search mode). DIA-NN search was used with oxidation of methionine residues and acetylation of the protein N-termini set as variable modifications and carbamidomethylation of cysteine residues as fixed modification. The match between runs feature was used for all analyses and the output (precursor) was filtered at 0.01 FDR and finally the protein inference was performed on the level of genes using only proteotypic peptides.
The Perseus software (version1.6.15.0) was used for the statistical analysis of the raw data. The baseline proteomic data set of the cohort was grouped based on active swollen and tender joints count upon 6 months abatacept therapy. The groupings were responders and non-responders. The dataset was Log2 transformed and filtered based on a threshold of at least 70% valid values in at least one of the three groups. The groupings were statistically evaluated using an ANOVA test. The final list of significantly altered proteins contained ten proteins with P-values less than 0.05.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD046112 [
21].
2.8. Statistics
Unpaired t test, one-way ANOVA or Wilcoxon matched-pairs test were applied in all experimental settings. All data were analyzed using GraphPad Prism v8 software. Differences were considered statistically significant at P-value < 0.05.
4. Discussion
Although several trials have proven the clinical effectiveness of CTLA4-Ig, only 40-50% of patients attain desired clinical responses. Responses to abatacept as well as to most of the biologic agents used to treat RA are unpredictable, since no clinically reliable prognostic markers are available to assist individual agent selections. In this study, we investigated whether the PB immunological profiling and serum proteome of RA patients may be used as a biomarker to predict clinical responses to abatacept. We found a strong association of high levels of Th1, MDSCs and DCs with a better response to abatacept treatment at 6 months. Moreover, baseline levels of 10 out of 303 proteins in peripheral serum showed differential expression according to clinical responses at 6 months. Interestingly a composite index based on the above described cellular and protein markers showed a high discriminative ability to predict response to abatacept.
Several animal and human studies have shown that both Th1 and Th17 effector T cells contribute to RA pathogenesis. Earlier studies have shown that IFN-γ producing CD4
+ cells are present in RA synovium [
25,
26], while more recent synovium RNA-seq analysis data support the presence of IFN-γ in early RA patients [
27]. Mostly animal studies strongly support the contribution of Th17 cells in synovial inflammation [
28,
29,
30]. On the contrary, FoxP3
+ Tregs detected in RA tissues have compromised function, highlighting that the imbalance between pathogenic effector T cells and regulatory T cells is an essential factor in the development of RA [
6]. Given that CTLA4 is expressed and regulates functionally the above effector and regulatory T cells, we herein questioned whether these subpopulations may affect response to abatacept. We found that responders had higher numbers of CD4
+IFN-γ
+ and a trend for higher CD4
+IL-17
+ T cells at baseline compared to patients resistant to abatacept. These findings indirectly support that the therapeutic effect of abatacept in RA patients in PB, could be partially attributed to a direct inhibitory effect on these effector T cells. Although several studies have assessed for soluble biomarkers to predict response to abatacept, only a few assessed whether immune cells in PB may be associated to clinical responses.
Microarray analysis of PB of RA patients before abatacept therapy revealed that high expression of type I IFN-related genes was associated with a clinical response to abatacept [
31]. In our study, we did not assess for type-I IFN levels, however we showed that IFN-γ-producing CD4
+ T cells are strongly associated with clinical response to abatacept in RA patients. Inamo
et al have shown that elevated levels of Th17 and Treg cells could be associated with clinical responses to abatacept in seropositive early RA patients [
32]. However, in our cohort seropositivity was not a co-factor – together with CD4
+IFN-γ
+ levels – which affected clinical responses. Differences in results in the above studies could be attributed to clinical and genetic differences between patients studied, stage of RA (early vs established) and methods applied to investigate peripheral immune cells.
Additionally, we found that CTLA4-Ig significantly reduces PB pathogenic CD4
+ T cells (
Figure 2b), nevertheless the effect was more significant and irrespective of clinical response for Th17 cells (
Figure 3f), while the reduction in CD4
+IFN-γ
+ was significant only for clinical responders (
Figure 3e). This finding is in accordance with data showing that abatacept decreased circulating CD8
+CD28
- T cells and other effector T cells [
33], as well as CD25
+ Tregs [
32]. Nevertheless, the mechanisms that abatacept mediates these effects have not been investigated.
MDSCs are a heterogeneous population of myeloid cells which have a regulatory role through increase of regulatory cells and inhibition of pathogenic T cells [
34]. Although some reports have shown that MDSCs can aggravate inflammatory arthritis in mice [
35,
36,
37], in most of the studies a regulatory role of MDCSs has been found. MDSCs play crucial role in the suppression of collagen-induced arthritis (CIA) mouse model, by inhibiting CD4
+ T cell proliferation, differentiation into Th17 cells and cytokine production [
38].
In vitro studies have indicated that MDSCs regulate Th17/Treg cells, decrease Th1 and Th17 cells while Tregs are increased via IL-10, ameliorating inflammatory arthritis
in vivo [
39,
40]. Furthermore, it has been found that RA patients with high disease activity have increased frequency of PB MDSCs and are also detectable in the synovial tissue of patients [
41,
42,
43]. However, the impact of abatacept on MDSCs as well as any association between the proportion of MDSCs and therapeutic response are unknown. In this study, although we found no depleting effect in MDSCs population after 6 months of abatacept (
Figure 2b, 4b), interestingly the baseline levels of MDSCs were positively associated with clinical responses at 6 months (
Figure 4a). It has been shown that MDSCs express low levels of CD80 [
37], but whether any effect of CTLA4-Ig on these cells could be mediated through reverse signaling in this mostly suppressive population is not known.
DCs are in the crossroad of immune responses through their crucial role in the antigen presentation process, and their dual role for both T cell activation or T cell inhibition. RA DCs have an increased capability to recruit macrophages, neutrophils and monocytes due to their enhanced secretion of chemokines, CXCL8 and CCL3, leading to exacerbated inflammation [
44]. Moreover, DCs exist in high concentrations in RA synovial joint tissues and secrete increased amounts of the pro-inflammatory cytokines IL-12, IL-23, inducing the generation of Th17 cells [
45]. Synovial DCs, also secrete high levels of chemokines CCL17, CXCL9 and CXCL10, attracting effector T cells which increase local inflammation, underling their role in the initiation of joint inflammation [
46]. There are clinical trials with
ex-vivo manipulated tolerogenic DCs in patients with RA, however they have limited efficacy [
47,
48]. Regarding abatacept, it is well studied that CTLA4-Ig induces its immunoregulatory function through binding to co-stimulatory molecules on DCs mediating a reverse signaling effect [
49,
50]. Our findings suggest that a higher frequency of DCs at abatacept initiation is positively correlated with response to abatacept therapy (
Figure 4c). This finding further supports that the anti-inflammatory effect in RA patients may at least partially mediated through the above-described reverse signaling on DCs. Interestingly, our finding that DCs’ number is stable during the 6 months of therapy (
Figure 2b, 4e), could support the concept of a continuous immunomodulatory effect through CTLA4-Ig on DCs which could be of clinical significance.
Herein, we performed serum proteomic analysis before abatacept treatment and we detected proteins considered to contribute to the inflammatory response in RA. Among the 303 proteins assessed, we concluded with 10 differentially expressed regarding the response to therapy after 6 months, and the majority of them are known to contribute to the pathogenesis of disease.
Interestingly, most of the proteins positively correlated to response have a pro-inflammatory role. Thus, CSF1R which binds to Colony-stimulating factor-1 (CSF-1) and promotes the proliferation and differentiation of myeloid cells, has been found to be expressed in synovial tissue of RA and PsA patients [
51]. Pharmacological inhibitors of CSF1R reduce the inflammatory activation of RA synovial tissue and severity of experimental arthritis [
52]. ADIPOQ which is a collagen-like protein shown to have a pro-inflammatory role in RA through the secretion of inflammatory mediators was found to be positively correlated to response [
53]. Interestingly, high levels of adiponectin are detected in serum and synovial fluid of RA patients [
54]. TNXB is a component of synovial membrane and immunofluorescence labeling in synovium from RA patients shows significant expression, suggesting that inflammatory mediators may increase TNXB production [
55]. Finally, CTBS is a lysosomal glycosidase and participates in degradation of glycoproteins, was also increased in serum of responders. A previous report revealed increased levels of CTBS in RA patients after tocilizumab treatment [
56].
Furthermore, 5 of the proteins detected in our analysis were shown to have a negative correlation with response to abatacept. ApoC3 is a regulator of lipoproteins and serves as a link between insulin resistance (IR) and beta-cell dysfunction that are present in RA patients. Notably, high ApoC3 serum levels have a positive and significant association with disease activity in RA patients [
57]. IGHV2-70D and IGHD are immunoglobulin heavy chain (IGH) proteins that influence B cell receptor altering response to infections. A variety of IGH proteins, as well as polymorphisms in IGH loci are associated with autoimmune diseases, like RA and systemic lupus erythematosus (SLE) [
58,
59]. Regarding properdin complement system, it has been found that neutralization of properdin has a protective role in development arthritis in mice models [
60]. Herein, we found that Complement Factor Properdin (CFP) is high in serum of non-responders to abatacept.
Although abatacept is an effective treatment for RA only 30-40% responds to therapy and individual responses are unpredictable. Research has been focused in the area of biomarkers to predict response in biological agents at the individual level, but reliable and clinically applicable biomarkers are not available. Herein we formulated a “composite index” based on the expression levels of the cellular signature and proteins found to be associated with clinical responses. The index was found to have a high discriminative value to predict response (AUC=0.93, 95% CI 0.83 to 1.00), with 90% sensitivity and 88.24% specificity for clinical responders (
Figure S1). Nevertheless, we were not able to confirm these results in a replication cohort neither to assess its specificity towards another biologic agent. Should these data be confirmed in a larger cohort would be of clinical value to assist individual treatment choices and further optimize RA treatment.
Author Contributions
Conceptualization, P.V. and P.S.; methodology, P.G., and G.P.; software, M.S.; validation, P.G.; formal analysis, P.G., and A.B.; investigation, P.G. and G.P.; resources, A.R., N.A., E.K., I.F., and G.B.; data curation, P.G., and J.Z.; writing—original draft preparation, P.G., and P.S.; writing—review and editing, M.S., P.V. and P.S.; visualization, P.G.; supervision, P.V. and P.S.; project administration, P.S.; funding acquisition, P.S. All authors have read and agreed to the published version of the manuscript.
Scheme 1.
and 6 months of CTLA4-Ig therapy using flow cytometry and proteomic analysis (Created with BioRender.com).
Scheme 1.
and 6 months of CTLA4-Ig therapy using flow cytometry and proteomic analysis (Created with BioRender.com).
Figure 1.
Proteomics analysis of RA patients sera at baseline. (a) Heatmap of the 303 proteins that are identified at the serum from the proteomic analysis. (b) Bubble plot of enrichment analysis representing enriched pathways in the set of identified proteins. The color of the bubbles represents the statistical significance.
Figure 1.
Proteomics analysis of RA patients sera at baseline. (a) Heatmap of the 303 proteins that are identified at the serum from the proteomic analysis. (b) Bubble plot of enrichment analysis representing enriched pathways in the set of identified proteins. The color of the bubbles represents the statistical significance.
Figure 2.
(a) Gating strategy for the identification of PB IFN-γ and IL-17 expressing CD4+ T cell (Th1, Th17, respectively), CD4+CD127-CD25+FoxP3+ T cells (FoxP3+ T cells), CD14+HLADRint/-CD15+CD33+ cells (MDSCs) and CD3-HLADR+CD1α+ cells (DCs). (b) Frequencies of above-mentioned cell populations in the PB of patients at 0 and 6 months post treatment (**P-value=0.0087, ***P-value= 0.0004 and ***P-value= 0.0002, respectively). Patients n=27. Statistical significance was obtained by Wilcoxon matched-pairs test.
Figure 2.
(a) Gating strategy for the identification of PB IFN-γ and IL-17 expressing CD4+ T cell (Th1, Th17, respectively), CD4+CD127-CD25+FoxP3+ T cells (FoxP3+ T cells), CD14+HLADRint/-CD15+CD33+ cells (MDSCs) and CD3-HLADR+CD1α+ cells (DCs). (b) Frequencies of above-mentioned cell populations in the PB of patients at 0 and 6 months post treatment (**P-value=0.0087, ***P-value= 0.0004 and ***P-value= 0.0002, respectively). Patients n=27. Statistical significance was obtained by Wilcoxon matched-pairs test.
Figure 3.
Evaluation of Th1, Th17 and FoxP3+ T cells frequencies in the PB of RA patients in response to abatacept. Frequencies of (a) CD4+ T cells (b) Th1, (c) Th17 and (d) FoxP3+ T cells at baseline (0 months) of responders and non-responders at 6 months of abatacept therapy (*P-value=0.0224, **P-value=0.009). (e) Th1 (f) Th17 and (g) FoxP3+ T cells percentages in the PB of patients at 0, 3 and 6 months post treatment (*P-value=0.0277 and **P-value= 0.0042). Classification of patients (responders or non-responders) based on the “Swollen joints”, 6 months after abatacept. RA patients with “Swollen joints”≤2 belong to low disease activity group (responders). Responders n=9, non-responders n=18. Statistical significance was obtained by un-paired t test and one-way ANOVA.
Figure 3.
Evaluation of Th1, Th17 and FoxP3+ T cells frequencies in the PB of RA patients in response to abatacept. Frequencies of (a) CD4+ T cells (b) Th1, (c) Th17 and (d) FoxP3+ T cells at baseline (0 months) of responders and non-responders at 6 months of abatacept therapy (*P-value=0.0224, **P-value=0.009). (e) Th1 (f) Th17 and (g) FoxP3+ T cells percentages in the PB of patients at 0, 3 and 6 months post treatment (*P-value=0.0277 and **P-value= 0.0042). Classification of patients (responders or non-responders) based on the “Swollen joints”, 6 months after abatacept. RA patients with “Swollen joints”≤2 belong to low disease activity group (responders). Responders n=9, non-responders n=18. Statistical significance was obtained by un-paired t test and one-way ANOVA.
Figure 4.
Evaluation of myeloid cell frequencies in the PB of RA patients in response to abatacept. Frequencies of MDSCs in the PB of patients at (a) baseline, (b) 3 and 6 months post treatment (*P-value=0.0322). Responders n=9, non-responders n=18. (c) Frequencies of DCs at baseline (0 months) and (d) 6 months of responders and non-responders (*P-value=0.0339 and 0.0438, respectively). (e) Frequencies of DCs in the PB of patients at 0, 3 and 6 months post treatment. Responders n=8, non-responders n=18. Statistical significance was obtained by unpaired t test.
Figure 4.
Evaluation of myeloid cell frequencies in the PB of RA patients in response to abatacept. Frequencies of MDSCs in the PB of patients at (a) baseline, (b) 3 and 6 months post treatment (*P-value=0.0322). Responders n=9, non-responders n=18. (c) Frequencies of DCs at baseline (0 months) and (d) 6 months of responders and non-responders (*P-value=0.0339 and 0.0438, respectively). (e) Frequencies of DCs in the PB of patients at 0, 3 and 6 months post treatment. Responders n=8, non-responders n=18. Statistical significance was obtained by unpaired t test.
Figure 5.
Proteomics analysis of sera, at baseline, of responders vs non- responders to abatacept. (a) Volcano plots of the differential expression analysis comparing responders vs non-responders. Upregulated proteins are denoted by red and downregulated by blue. Proteins not reaching our significance threshold (P-value <0.05) depicted in gray. (b) Enrichment analysis of Kegg pathways for differentially expressed proteins of responders vs non-responders in 303 proteins background, sorted by P-value ranking. Color gradient correlates with statistical significance. Pathways not reaching significance threshold (P-value <0.05) are shown in gray.
Figure 5.
Proteomics analysis of sera, at baseline, of responders vs non- responders to abatacept. (a) Volcano plots of the differential expression analysis comparing responders vs non-responders. Upregulated proteins are denoted by red and downregulated by blue. Proteins not reaching our significance threshold (P-value <0.05) depicted in gray. (b) Enrichment analysis of Kegg pathways for differentially expressed proteins of responders vs non-responders in 303 proteins background, sorted by P-value ranking. Color gradient correlates with statistical significance. Pathways not reaching significance threshold (P-value <0.05) are shown in gray.
Table 1.
Clinical characteristics of patients.
Table 1.
Clinical characteristics of patients.
|
All |
Responders |
Non-responders |
Patients (n) |
29 |
10 |
19 |
Gender (% females) |
25 (86.2%) |
7 (70.0%) |
18 (94.7%) |
Age,Median (min,max;IQR)
|
64 (32,77;17) |
58 (32,71;20) |
67 (40,77;18) |
Disease duration,Median (min,max;IQR)
|
23 (3,240;45) |
32.5 (8,120;50) |
22 (3,240;40) |
R.F. positive (%) |
8 (27.6%) |
2 (20.0%) |
6 (31.6%) |
Anti-CCP positive (%) |
15 (51.7%) |
5 (50.0%) |
10 (52.6%) |
Methotrexate (%) |
26 (89.65%) |
9 (90%) |
17 (89.47%) |
Steroids (%) |
10 (34.48%) |
4 (40%) |
6 (31.57%) |
Disease characteristics before abatacept therapy (Baseline) |
ESR,Median (min,max;IQR)
|
26 (4,77;34) |
25.5 (4,69;34) |
27 (4,77;40) |
CRP,Median (min,max;IQR)
|
0.33 (0.24,6.03;0.53) |
0,52 (0.31,3.61;0.70) |
0.33 (0.24,6.03;0.53) |
Swollen 28,Median (min,max;IQR)
|
6 (2,18;5) |
9 (3,15;5) |
6 (2,18;4) |
Tender 28,Median (min,max;IQR)
|
7 (0,23;6) |
9 (2,23;8) |
6 (0,20;6) |
VAS global,Median (min,max;IQR)
|
70 (40,100;20) |
70 (40,100;60) |
60 (40,100;20) |
DAS28-ESR,Median (min,max;IQR)
|
5.48 (3.42,6.85; 1.90) |
5.91 (3.42,6.79;1.89) |
4.63 (3.94,6.85;1.89) |
Disease characteristics after abatacept therapy (6 months) |
ESR,Median (min,max;IQR)
|
28 (6,68,21) |
24.5 (8,68;22) |
31 (6,59;18) |
CRP,Median (min,max;IQR)
|
0.36 (0.10,3.01;0.27) |
0.36 (0.10,1.45;0.18) |
0.36 (0.30,3.01;0.30) |
Swollen 28,Median (min,max;IQR)
|
4 (0,17;6) |
1.5 (0,2;2) |
6 (3,17;8) |
Tender 28,Median (min,max;IQR)
|
2 (0,27;6) |
1 (0,4;2) |
3 (0,27;9) |
VAS global,Median (min,max;IQR)
|
30 (0,100;35) |
30 (0,60;23) |
40 (10,100;50) |
DAS28-ESR,Median (min,max;IQR)
|
4.20 (2.22,8.04;2.13) |
2.93 (2.22,4.61;1.80) |
4.67 (2.85,8.04;2.05) |
ΔDAS28,Median (min,max;IQR)
|
1.30 (-1.84, 4.30; 2.61) |
2.42 (0.16, 4.30;2.48) |
0.48 (-1.84,3.44;2.76) |
All values are medians (IQR) unless otherwise specified. R.F.: Rheumatoid Factor; anti-CCP: Anti-cyclic Citrullinated Peptide; ESR: Erythrocyte Sedimentation Rate; CRP: C-reactive protein; VAS: Visual Analogue Scale; DAS28: Disease activity score using 28 joints |
Table 2.
Baseline association of cell populations with clinical characteristics.
Table 2.
Baseline association of cell populations with clinical characteristics.
|
TH1 % |
TH17 % |
FoxP3 % |
MDSCs % |
DCs % |
|
Mean (SD) |
Mean (SD) |
Mean (SD) |
Mean (SD) |
Mean (SD) |
Seropositive |
|
|
|
|
|
No |
27.9 (13.4) |
2.32 (1.43) |
4.58 (4.40) |
0.25 (0.22) |
1.18 (1.19) |
Yes |
18.8 (15.8) |
1.62 (1.08) |
3.05 (2.43) |
0.77 (0.20) |
2.07 (2.21) |
P-value |
0.120 |
0.167 |
0.273 |
0.050 |
0.224 |
ESR (P-value) |
-0.237 (0.233) |
0.028 (0.890) |
-0.190 (0.343) |
-0.143 (0.476) |
-0.202 (0.321) |
CRP (P-value) |
0.083 (0.681) |
-0.087 (0.271) |
-0.125 (0.535) |
0.261 (0.188) |
0.081 (0.692) |
Swollen 28 (P-value) |
0.346 (0.077) |
-0.220 (0.271) |
-0.051 (0.799) |
-0.151 (0.452) |
-0.124 (0.547) |
Tender 28 (P-value) |
0.386* (0.047) |
0.173 (0.389) |
0.440* (0.022) |
-0.348 (0.075) |
-0.133 (0.518) |
VAS global (P-value) |
-0.102 (0.614) |
0.097 (0.631) |
0.085 (0.674) |
-0.274 (0.167) |
-0.271 (0.180) |
DAS28-ESR (P-value) |
0.184 (0.359) |
0.195 (0.329) |
0.272 (0.169) |
-0.363 (0.063) |
-0.204 (0.317) |
Independent samples T-test (comparison of means) was performed on the upper part of the table and Spearman's correlation coefficient on the lower part. *Significant associations |
Table 3.
Description of the 10 differentially expressed proteins at baseline in responders compared to non-responders.
Table 3.
Description of the 10 differentially expressed proteins at baseline in responders compared to non-responders.
Protein Name |
Ratio Responders/ non-Responders |
P-value |
Function |
CSF1R (Macrophage colony-stimulating factor 1 receptor) |
2,21997 |
1,4471 |
Biological regulation; cell differentiation |
GC (Vitamin D-binding protein) |
2,19503 |
1,42407 |
Lipid metabolic process |
ADIPOQ (Adiponectin) |
2,11482 |
1,35079 |
Adiponectin-mediated signaling pathway |
TNXB (Tenascin-X) |
2,09943 |
1,33687 |
Actin cytoskeleton organization |
CTBS (Chitobiase, Di-N-Acetyl-) |
2,93148 |
2,14787 |
Amine metabolic process |
APOC3 (Apolipoprotein C-III) |
-2,60697 |
1,81868 |
Acylglycerol metabolic process |
IGHV2-70D (Immunoglobulin heavy variable 2-70D) |
-2,20792 |
1,43596 |
Immune response |
IGHD (Immunoglobulin heavy constant delta) |
-2,11844 |
1,35407 |
Immune response |
BTD (Biotinidase) |
-2,02902 |
1,27382 |
Metabolic process |
CFP (Complement factor properdin) |
-2,01735 |
1,26348 |
Activation of immune response |