Skip to main content

Lower peripheral helper T cell levels in the synovium are associated with a better response to anti-TNF therapy in rheumatoid arthritis

Abstract

Background

The mechanisms by which only some rheumatoid arthritis (RA) patients respond favorably to TNF blockade are still poorly characterized. The goal of this study was to identify biological features that explain this differential response using a multilevel transcriptome analysis of the synovial membrane.

Methods

Synovial samples from 11 patients on anti-TNF therapy were obtained by arthroscopy at baseline and week 20. Analysis of the synovial transcriptome was performed at the gene, pathway, and cell-type levels. Newly characterized pathogenic cell types in RA, peripheral helper T cells (TPH), and CD34-THY1+ fibroblasts were estimated using a cell-type deconvolution approach. TPH association was validated using immunofluorescence. External validation was performed on an independent dataset.

Results

After multiple-test correction, 16 and 4 genes were differentially expressed at baseline and week 20, respectively. At the pathway level, 86 and 17 biological processes were significantly enriched at baseline and week 20, respectively. Longitudinal expression changes were associated with a drastic decrease of innate immune activity (P < 5e−30), and an activation of the bone and cartilage regeneration processes (P < 5e−10). Cell-type deconvolution revealed a significant association between low TPH cells at baseline and a better response (P = 0.026). Lower TPH cells were maintained in good responders up to week 20 (P = 0.032). Immunofluorescent analyses confirmed the accuracy of the cell-type estimation (r2 = 0.58, P = 0.005) and an association with response. TPH association with anti-TNF response was validated in an independent sample of RA patients (P = 0.0040).

Conclusions

A lower abundance in the synovial membrane of the pathogenic T cell type newly associated with RA, peripheral helper T lymphocyte, is associated with a good response to anti-TNF therapy. Major changes in the myeloid cell compartment were also observed in response to therapy. The results of this study could help develop more effective therapies aimed at treating the pathogenic mechanisms in RA that are currently not well targeted by anti-TNF agents.

Background

Anti-TNF therapies have been a major breakthrough in the management of rheumatoid arthritis (RA) [1]. To date, disease remission has become an attainable goal for many patients [2]. However, 30 to 40% of anti-TNF-treated patients do not show a significant clinical improvement [3]. The reasons behind this differential response are still elusive [3]. Clinical and basic studies have attempted to identify key factors, with very limited success [4]. Identifying the determinants of anti-TNF response would be of high value: these could be used to stratify patients, guide drug selection, and improve disease management, and could help to develop more efficacious therapies, either alone or in combination with anti-TNF agents.

The synovial membrane is the target tissue of chronic inflammation in RA [5], and its analysis should provide highly valuable information on the biological mechanisms underlying the differential response to TNF blocking agents. However, obtaining synovial biopsies is highly invasive and consequently, very few studies have been performed on this key issue with regard to drug efficacy. Using in-house two-color microarrays, van der Pouw Kraan [6] and Lindberg et al. [7] analyzed the baseline transcriptome of RA patients starting infliximab therapy, but neither study found differentially expressed genes. The higher technical variability associated with these early analysis platforms could have undermined the power to accurately identify transcript changes [8]. Using a single-color commercial array, Badot et al. analyzed the longitudinal change in expression of n = 8 adalimumab-treated patients [9] and found some expression changes mainly in genes related to cell division. However, differential expression was restricted to 14% of the analyzed transcripts, thereby limiting the global screening capacity to a minor fraction of the total biological variability.

The cellular composition of tissues and the extent of infiltration of immune cell types has proven to be predictive of disease subtypes of complex diseases like cancer [10]. Very recently, studies based on single-cell approaches in the synovial membrane have characterized new pathogenic subsets associated with RA. Using single-cell RNA-seq, Mizoguchi et al. [11] classified RA synovial fibroblasts into three different subsets, with the CD34-THY1+ subtype showing a significant expansion in comparison to osteoarthritis patients. Similarly, using mass-cytometry on T cells isolated from the synovium, Rao et al. [12] found an expanded population of PD-1hiCXCR5-CD4+T cells in RA. Due to their capacity to provide B cell help, this new T cell subset was described as the peripheral helper T cell (TPH). The relationship of these new pathogenic cell subsets with the response to anti-TNF therapy in RA has not yet been addressed.

In the present study, we have performed a multilevel analysis of the synovial membrane transcriptome to identify critical biological mechanisms associated with response to anti-TNF therapy in RA. Using a longitudinal cohort of patients starting anti-TNF treatment, we identify differentially expressed genes, pathways, and cell types associated with the clinical response at week 20. In the latter analysis, using a cell-deconvolution approach, we find that a lower TPH infiltration is associated with a good response to therapy, and these lower levels are maintained through time. A lower baseline presence of myeloid immune cells is also associated with a favorable response and drastic reduction of this cell compartment is the main biological change observed in responders. The results of this study provide a comprehensive look at the biological mechanisms associated with anti-TNF response in RA.

Methods

Patients and samples

Eleven patients (8 women, 3 men) fulfilling the American College of Rheumatology/European League Against Rheumatism 2010 criteria for RA [13], with a basal mean Disease Activity Score in 28 joints (DAS28) [14] of (5.3 [4.2–6.9]) (median and interquartile range [IQR]), non-responder to methotrexate therapy and who were prescribed an anti-TNF by their rheumatologist (n = 6 infliximab, n = 3 adalimumab, n = 2 etanercept), were enrolled in this study. At the time of enrolment, all patients were naïve to any biological therapy. Response to treatment was defined by the EULAR response criteria [15] after 20 weeks of therapy, aggregating moderate and good responders into a single responder group. During the 20-week treatment with anti-TNF, all patients followed the standard treatment pattern, with no adjustments. The study was approved by the local ethics committee, and all patients signed the informed consent. This study was conducted in accordance with the Declaration of Helsinki principles.

Synovial samples were obtained by guided arthroscopy using a 2.7-mm arthroscope (Storz, Tuttlingen, Germany) under local anesthesia. In all patients, 6–8 biopsies were taken from the suprapatellar pouch and medial and lateral gutters with a 3-mm grasping forceps before starting a TNF antagonist. All patients underwent a second arthroscopy with synovial sampling after 20 weeks of therapy. Synovial biopsies for mRNA analysis were immediately stored in the RNALater preserving agent (QIAGEN, USA) and frozen to − 80 °C until RNA extraction [16].

Microarray analysis

RNA was extracted from synovial biopsies using the RNA Mini Kit (Qiagen, USA), and the integrity was assessed using BioAnalyzer microfluidic gel analysis (Agilent, USA). All samples were of high quality (RNA Integrity Number > 8) and were subsequently analyzed using Sentrix whole genome Beadchips WG6 version 2 (Illumina, US). Briefly, after RNA isolation, biotin-labeled cRNA was prepared using the Ambion Illumina RNA amplification kit (Ambion, US) and Illumina TotalPrep RNA Amplification Kit (Ambion, US). Biotin-labeled cRNA (1.5 μg) was hybridized to WG6 v2 Beadchips and scanned on the 500x Illumina BeadStation. Data collection was performed using BeadStudio 3.1.1.0 software (Illumina, USA).

To perform an independent validation of the cell type associated with treatment response, publicly available transcriptomic datasets (microarray or RNA-seq) were searched in the Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/geo/). Using the keywords “synovial membrane” and “rheumatoid arthritis,” a total of 62 datasets were found that were performed on human samples. From these, after excluding in vitro culture studies and two-color microarrays, only two datasets, GSE57376 [17] and GSE15602 [9] on anti-TNF treatment were found. The former dataset was discarded since it was only based on 3 RA patients and had no information on clinical response. The latter dataset consisted of synovial membrane biopsies of n = 11 RA patients (3 Good, 5 Moderate and 3 None EULAR responders) were collected at week 12 of adalimumab treatment.

Gene-level and pathway-level association with anti-TNF response

After log2 transformation, raw gene expression data was quantile-normalized. Batch adjustment was performed using ComBat [18]. Differential gene expression was assessed using the t-test, and the resulting P values corrected for multiple testing using Bonferroni. The biological pathway association was tested using the Gene Ontology (GO) enrichment approach [19] implemented in the clusterProfiler R package. Genes with a nominally significant differential expression and a fold change > 1.2 were selected as input for this analysis; Bonferroni multiple test correction was applied to determine significant GO terms.

RT-QPCR analysis of associated genes

Real-time quantitative PCR analyses of the genes most significantly associated with treatment response and marker genes for PD-1hiCXCR5-CD4+ TPH cells were performed using Taqman assays (Life Technologies, USA). In the former group, FAM-labeled assays were used to measure the expression of PIK3CD (Hs00192399_m1) and CX3CL1 (Hs00171086_m1) genes, while TIMELESS (Hs01086960_g1) and CTLA4 (Hs01011591_m1) were used as TPH-associated markers. Ribosomal Protein L11 (VIC-labeled assay, Hs00831112_s1) was used as an endogenous control gene, and all assays were performed using the Taqman Universal PCR master mix (Life Technologies, USA). Gene expression analysis was performed using the recommended protocol with the Applied Biosystems HT7900 system (ABI, USA), and fold-change expression was calculated from the Ct values of the test and endogenous genes using the -2ΔΔCt method. The correlation between microarray gene expression values and RT-PCR values were calculated using Pearson’s product-moment correlation.

Deconvolution analysis of pathogenic cell subsets

In order to estimate the enrichment of relevant T cell and fibroblast subsets in the synovial membrane samples, we adapted the SPEC (Subset Prediction from Enrichment Correlation) method [20]. Briefly, this computational method uses a gene signature characteristic of specific cell types to estimate their contribution to the observed gene expression in a complex tissue like blood. The ranking of the cell-type signature genes within the global ranked list of genes is then used to compute the enrichment score with an adaptation of the Gene Set Enrichment Method [21]. The enrichment score corresponds to the maximum deviation from zero of a running sum statistic calculated by running down the ranked list of genes. This approach has proven to have a high accuracy in determining cell subsets within PBMC RNA analysis. Here, this deconvolution method was adapted to estimate fibroblast and T cell subtypes from synovial membrane transcriptomic data. To estimate the enrichment of the fibroblast types, we determined the top set of genes overexpressed in each subset. For this objective, we used the microarray gene expression data generated from the three isolated cell types by Mizoguchi et al. [11] and available at the GEO ID GSE107105 (supplementary Table 1). In the case of PD-1hiCXCR5-CD+ T cells, we used the set of n = 54 significantly overexpressed genes in TPH reported in the original discovery study [12] (supplementary Table 2).

Immunofluorescence analysis of TPH in RA synovial samples

Immunofluorescence analysis of the TPH lymphocytes was performed following the same methodology as described in the original discovery study [12]. Briefly, OCT frozen sections (7 μm) of the synovial membrane biopsies of the 11 patients treated with anti-TNF therapy and at both time points (baseline and week 20) were fixed with ice-cold acetone. Quenching of endogenous peroxidase was performed using Peroxidase-Blocking Solution (ref: S2023, Dako-Agilent). Blocking of unspecific unions was done using 5% goat normal serum (Life technologies, 16210064) and 2.5% bovine serum albumin (Sigma-Aldrich, 10735078001). Primary antibodies for PD-1 (NAT105,– ref.: 760-4895, Roche), CD4 (SP35, ref: 790-4423, Roche) and CXCR5 (51505, – ref: MAB190, R&D Systems) were used RTU for PD-1 and CD4 and 1:50 for CXCR5 to detect the presence of PD-1hiCXCR5-CD4+ cells. As secondary antibodies, Alexa Fluor 488 goat-anti mouse IgG1 (ref: A21121, ThermoFisher), Alexa Fluor Plus 647 goat anti-rabbit IgG (ref: A32733, ThermoFisher), and Biotinilated goat anti-mouse IgG2b (ref: ab98701, Abcam) plus Streptavidin Texas Red® (NEL721001EA, PerkinElmer) were used. Tissue samples were stained with DAPI (D9542, Sigma) and mounted with Fluorescence Mounting Medium (S3023, Dako-Agilent). The specificity of the primary antibodies was determined using rabbit and mouse IgG isotype controls (ref: ab27478 and ab37355, Abcam) and by the omission of the primary antibody.

Digitally scanned fluorescent images were acquired using a NanoZoomer-2.0 HT C9600 scanner (Hamamatsu, Photonics, France) with a ×20 objective and coupled with a mercury lamp unit L11600-05 and using NDP.scan2.5 software U10074-03 (Hamamatsu, Photonics, France). All images were visualized with the gamma correction set at 1.0 and the sharpen filter enabled in the image controls panel of the NDP.view 2 U12388-01 software (Hamamatsu, Photonics, France). Image analysis was performed using QuPath software [22]. Cell detection was performed using the cell detection algorithm based on the DAPI nuclei fluorescence and subsequently, cells were segmented between positive and negative according to the membranous labeling by the three different markers (Alexa 488 for PD-1, Texas Red® for CXCR5 and Alexa 647 for CD4). For each marker, isotype controls as well as negative samples (i.e., samples where the primary antibodies were omitted) were used to set the threshold values for each marker.

Results

From the 11 RA patients, 8 responded to anti-TNF treatment and 3 patients showed no significant clinical response at week 20 of therapy. The clinical and demographic data of the longitudinal patient cohort are detailed in Table 1.

Table 1 Clinical and demographical characteristics of the RA patient cohort

Single-gene and pathway association with response

At week 0, n = 16 genes were differentially expressed between responders and non-responders (Table 2), from which 11 were overexpressed in responders and 5 in non-responders. At week 20, four genes were differentially expressed after Bonferroni multiple test correction (Table 2), and one of these genes was overexpressed in responders (PDE4B, P = 1.09E-09) while the remaining 3 were overexpressed in non-responders. At the nominal level, 7 out of the 11 differential genes in the synovial membrane at week 0 continued to be differentially expressed between responders and non-responders at week 20 (P < 0.05, the same direction of change, Supplementary Table 3). The remaining 5 genes, while not reaching statistical significance, preserved the same direction of change as found at baseline (supplementary Figure 1). Conversely, none of the four differentially expressed genes at week 20 were differential at week 0 (Supplementary Table 3). This result suggests that the biological features that cause unresponsiveness to TNF blocking therapy prevail throughout anti-TNF treatment. Using RT-PCR, we validated mRNA expression of the gene most differentially expressed at baseline in responders (i.e., CX3CL1, correlation P value = 2.11e−11) and the gene most highly expressed in non-responders (PIK3CD, correlation P value = 2.42e−10) (Supplementary Figure 2).

Table 2 List of genes differentially expressed between anti-TNF responders and non-responders

At the pathway level, at week, 0 n = 16 GO terms were significantly enriched in responders, and n = 70 GO terms were enriched in non-responders (Table 3 and Supplementary Table 4). At week 20, n = 17 GO terms were significantly enriched, and all due to genes overexpressed in anti-TNF non-responders (Table 3 and Supplementary Table 5). From these, three gene ontologies were significantly enriched at baseline in non-responders: GO:0043299|leukocyte degranulation (adjusted Pwk0 = 9.8e−4, adjusted Pwk20 = 3.97e−11), GO:0002275|myeloid cell activation involved in immune response (adjusted Pwk0 = 3.1e−4, adjusted Pwk20 = 3.6e−10) and GO:0002444|myeloid leukocyte mediated immunity (adjusted Pwk0 = 6.98e−3, adjusted Pwk20 = 6.36e−10). This result suggests that the innate immune system composition in the synovial membrane conditions the response to anti-TNF therapy, with non-responders bearing a higher infiltrate of myeloid lineage cells (i.e., neutrophils and macrophages). Comparing the longitudinal changes in biological pathways in the responder group, we found a profound modification with n = 149 GO terms significantly enriched from week 0 to week 20 (Table 3 and Supplementary Table 6). Most of the changes are due to the inactivation of pathways (n = 123 GO terms, 82.6%), but also several biological processes are activated due to anti-TNF efficacy (n = 26 GO terms, 17.4%). In the former group, we found that the biological changes most associated in responders corresponds to the same three GO terms that are consistently different when comparing responders vs non-responders at the two time points: GO:0043299|leukocyte degranulation (adjusted P = 3.26e−32), GO:0002275|myeloid cell activation involved in immune response (adjusted P = 1.14e−31) and GO:0002444|myeloid leukocyte mediated immunity (adjusted P = 2.99e−30). This result corroborates that TNF efficacy is mainly mediated by the drastic inactivation of the innate immune system. Of relevance, the pathways that become activated due to TNF response include biological processes strongly related to the production of the specific components of the joint stromal tissue, including GO:0001501|skeletal system development (adjusted P = 2.11e−10), GO:0001503|ossification (adjusted P = 4.04e−8), GO:0001649|osteoblast differentiation (adjusted P = 2.61e−4), GO:0051216|cartilage development (adjusted P = 8−42e−4) and GO:0002063|chondrocyte development (adjusted P = 8.99e−4). This biological evidence shows that the effective clearance of inflammation is paralleled by the activation of joint tissue reconstruction programs.

Table 3 Top associated Gene Ontologies with the anti-TNF response at baseline and week 20 and with the longitudinal change

Cell type deconvolution of RA pathological cell subtypes

Using single-cell characterization technologies, recent studies have identified distinct cell subsets within RA synovium. In particular, two new subtypes peripheral helper CD4+ T cells and THY1+ synovial fibroblasts have been found to be key pathological types. Based on the genes that are the characteristic of each subtype and using a deconvolution approach, we estimated the relative proportion of each cell type. We found at baseline, response to anti-TNF therapy was associated with having lower numbers of TPH cells (P = 0.021, Fig. 1a). We did not find evidence that THY1+ fibroblasts were associated with the response to anti-TNF therapy. Instead, we observed that the baseline proportions of both CD34+THY1- and CD34-THY1- (i.e., double negative) fibroblasts were associated with a response at week 20 (P = 0.027 and P = 0.021, respectively, Supplementary Figure 3). At the time of response (week 20), we found that all three associated cell types (CD34+ and CD34-THY1- fibroblasts and TPH) showed the same trend of association as identified at baseline, although the difference was not statistically significant (P > 0.05, supplementary Figure 4).

Fig. 1
figure1

The TPH enrichment in the synovial membrane of RA patients at baseline therapy is associated with the clinical response. a PD-1hiCXCR5-CD4+ cell enrichment levels at baseline estimated using the cell-type deconvolution approach in our patient cohort. The estimated enrichment scores for TPH cells are significantly associated with the EULAR response at week 20 (P = 0.021). b PD-1hiCXCR5-CD4+ cell percentage at baseline quantified using immunofluorescence in our patient cohort. As seen in the cell-type deconvolution approach, lower TPH cells are associated with better responses to anti-TNF agents (P = 0.025). In both cases, moderate and good EULAR responders show distinctively lower proportions of this pathogenic T cell subset associated with RA

To perform technical validation of the observed gene expression associated with TPH cells, we performed RT-PCR analysis of TPH-marker genes CTLA4 and TIMELESS. We found a highly significant correlation between the microarray and the RT-PCR gene expression quantification (Pearson correlation coefficients r = 0.89, P = 3.6e−8 and r = 0.65 P = 0.0009 for CTLA4 and TIMELESS, respectively) (Supplementary Figure 5).

Immunofluorescence analysis of TPH cells

In order to corroborate the association between TPH enrichment and anti-TNF response, we performed immunofluorescence analysis on the synovial membrane biopsies from the same set of patients (Fig. 2). All samples were positive for all markers which presented a membranous pattern. From them, CXCR5 expression was minimal in comparison with PD-1 and CD4. Anti-TNF non-responder patients presented higher groups of PD-1hiCXCR5-CD4+ T cells (Fig. 1b, P = 0.025). We found a significant correlation between the proportions of PD-1hiCXCR5-CD4+ T cells estimated through our approach and the actual abundance of this pathogenic T cell subtype in the RA synovia (r2 = 0.58, P = 0.0051, Supplementary Figure 6).

Fig. 2
figure2

Representative immunofluorescence of TPH in synovial tissues of a good responder (a, c, e, g) and a non-responder (b, d, f, h) at week 20. c, d PD1-positive cells in green, e, f CXCR5-positive cells in orange, and g, h CD4-positive cells in red. a and b show the combined immunofluorescence image of a good and bad responder to anti-TNF therapy, respectively. White arrows indicate synovial membrane aggregations of CD4+PD1+ T cells

Having confirmed the association of TPH cells with the response to anti-TNF therapy, we applied the cell-type deconvolution method to an external microarray dataset of RA patients treated with anti-TNF therapy [9]. In the transcriptional data at week 12 of therapy (n = 3 Good, 5 Moderate, and 3 None), we found a highly similar TPH profile to that of patients at week 20 (P = 0.029 and P = 0.00019, respectively) (Fig. 3). In both cases, good responders show a markedly lower frequency of the pathogenic T cell subset (P = 0.00035 and P = 0.0040, t test comparing TPH cell levels in good vs moderate and none patients in our patient cohort and the external RA patient cohort, respectively).

Fig. 3
figure3

TPH enrichment analysis in the synovial membrane of the external dataset reproduces the patterns observed in our patient cohort. a PD-1hiCXCR5-CD4+ cell percentage quantified using immunofluorescence on synovial biopsies of our patient cohort at week 20 of anti-TNF therapy. b PD-1hiCXCR5-CD4+ cell enrichment levels estimated using the cell-type deconvolution approach on the external patient cohort at week 12 of anti-TNF therapy. In both cases, EULAR good responders show distinctively lower TPH cell levels than moderate or none responders (P = 3.5e−4 and P = 4e−3 for our patient cohort and the external RA patient cohort, respectively)

Discussion

The reasons behind the variable efficacy of anti-TNF agents in RA have been elusive. Using single-gene, pathway and cellular deconvolution analyses on longitudinal transcriptomic data from the synovial membrane of RA patients treated with anti-TNF therapy, we have found key biological determinants of the response to this drug. At the gene level, baseline overexpression of the chemokine CX3CL1 and downregulation of leukocyte signaling kinase PIK3CD were significantly associated with a better response. At the pathway level, anti-TNF non-responders had a higher activation of myeloid cell immune response processes at baseline. A drastic downregulation of these same innate immunological processes was strongly associated with the efficacy of the drug at week 20. At the cellular level, lower TPH cells were associated with a good response to anti-TNF agents, both at baseline and at week 20. The results from this study provide a multilevel description of the biological determinants of anti-TNF response in RA.

Recent studies based on single-cell analysis techniques have characterized new pathogenic cell subtypes associated with RA. Using single-cell mass cytometry, PD-1hiCXCR5-CD4+ T (TPH) cells were found to be significantly expanded in the synovial membrane of seropositive RA patients. Using the genes characterizing this T cell subset, we were able to estimate the presence of a significantly lower number of TPH cells in responders to anti-TNF therapy. This difference was clearly maintained up to week 20 in the good responder group of patients and validated in an independent cohort of patients. The transcriptional similarity between PD-1hiCXCR5-CD4+ T and follicular helper T (TFH) and the in vitro capacity to induce differentiation of memory B cells suggests a peripheral helper phenotype. Our results suggest that an important part of this TPH pathogenic activity in RA is associated with mechanisms that are independent of TNF. TPH may act by recruiting TFH and B cells, promoting local autoantibody production that may not be reflected in the serum, leading to pathogenic immune complex fixation. Furthermore, activated B cells can polarize into different subtypes to produce an array of proinflammatory cytokines other than TNF [23], including interferon-gamma, IL-6, and IL-17. The latter cytokine, for example, promotes the accumulation of neutrophils in the synovial fluid, where they produce large amounts of cartilage-degrading enzymes and reactive oxygen species [24]. This way, in patients with higher levels of TPH, the inflammation will have alternative mechanisms to progress despite the reduction of TNF levels by anti-TNF agents.

Our results also suggest that the stromal component of the synovial membrane participates in the favorable response to anti-TNF therapy. Single-cell analysis of the synovial membrane identified three different synovial fibroblast types [11]: CD34-THY1-, CD34-THY1+, and CD34+. CD34-THY1+ fibroblasts have been shown to be expanded when compared with osteoarthritis synovium and therefore have been identified as the principal pathogenic subset. In our cell type analysis, we did not find evidence for enrichment of CD34-THY1+ fibroblasts with anti-TNF response. Instead, we observed a significant increase of CD34+ and CD34-THY1- fibroblasts. As described previously [25], THY1 is also highly abundant in vascular cell types, and therefore, an immunofluorescence-based quantification of this fibroblast subtype was ruled out. Therefore, these results should be taken with caution. Nonetheless, functional evidence indicates that CD34+ fibroblasts have a leading role in the recruitment of monocytes into the inflamed synovial, producing large amounts of IL6, CXCL2, and CCL2 in response to TNF stimulation. In comparison, CD34-THY1-fibroblasts produce higher levels of matrix metalloproteinases MMP1 and MMP3. According to our results, both fibroblast types would cause more TNF-sensitive stromal cells, in which a reduction in local levels of this cytokine by therapy would be translated into a less proinflammatory and matrix-degrading activity and, ultimately, into a higher likelihood of clinical response. Recently, Croft et al. [25] have further refined the classification of synovial fibroblasts in RA by identifying fibroblast activation protein α (FAPα+) as a key stratification marker. Like in the previous study, THY1-fibroblasts were shown to localize in the synovial lining layer and specifically correlate with cartilage-damaging features. Instead, FAPα+ THY1+ fibroblasts were shown to correlate with markers of systemic and tissue inflammation. Studies integrating this refined fibroblast characterization into the association with the response to anti-TNF therapy in RA are therefore warranted.

At the single gene level, we found two genes, CX3CL1 and PIK3CD, showing distinctively high expression levels at baseline in responders and non-responders, respectively. From these, only CXC3L1 remains overexpressed in the synovial membrane of responders at week 20. CXC3L1, also known as fractalkine, is a chemokine that can be present as membrane-bound or in soluble forms and is highly expressed in endothelial cells activated by TNF [26]. Interaction of membrane-bound fractalkine with its receptor CX3CR1 mediates the capture of circulating immune cells and infiltration into the inflamed tissue. In particular, CX3CR1hi monocytes that patrol endothelium have been shown to rapidly accumulate in sites after activation of the endothelium [27]. Accordingly, our results are indicative of patients with vascular tissue richer in CX3CL1 expression are more likely to have a beneficial impact of TNF blocking.

PIK3CD is highly expressed in the synovial membrane of patients that will not respond to therapy. Together with PIK3CA and PIK3CB, PIK3CD is a class I phosphoinositide 3-kinase, which is a group of enzymes that phosphorylate inositol lipids and participates in different signaling processes [28]. While PIK3CA and PIK3CB are expressed in many different tissues, PIK3CD is limited to the cells of the immune system including T and B lymphocytes [29], neutrophils [30], and macrophages [31]. More recently, PIK3CD has also been found to be highly expressed in RA synovial fibroblasts [32]. In vitro studies have shown that TNF is a major inducer of PIK3CD expression in synovial fibroblasts but not in leukocytes [33]. In these previous studies, the expression of PIK3CD in RA synovial fibroblasts was shown to induce cell growth by sensitizing fibroblasts to platelet-derived growth factor (PDGF) [32]. Based on this evidence, non-responders would have a stromal component more prone to PDGF-derived hyperplasia. For this group of patients, a combination of drugs interfering with PDGF signaling might be an efficacious approach to rescue non-responder patients. Experimental and clinical evidence of the therapeutic potential of PDGF pathway interference in RA is clearly supported [34].

The pathway-level analysis of anti-TNF therapy response revealed an important contribution of the myeloid cell load. Our longitudinal data demonstrates that a high myeloid cell activity not only conditions for a poorer response, but also, that it is the biological process which is more highly altered by the drug’s efficacy. The existence of two patient groups in RA with different myeloid leukocyte activation is in high accordance with recent findings using a single-cell RNA-seq analysis. Analyzing the synovial membrane cell composition of a relatively large patient cohort of RA and osteoarthritis individuals, Zhang et al. [35] identified two subtypes of patients in RA, defined as leukocyte-poor RA and leukocyte-rich RA. Our study extends these findings into the therapeutic domain. Despite not showing distinguishable clinical activities at baseline, non-responder patients have a significantly higher myeloid activity, which is one of the top-ranking biological processes detected in the leukocyte-rich RA patients. It is important to note that, using a standard histological definition of synovial tissue (i.e., lymphoid/myeloid/fibroid), we found no association with a clinical response or with TPH proportions (data not shown). This suggests that more precise, molecular quantification techniques like transcriptomic analysis are needed to capture this patient variability. In this regard, the pathway analysis of genes associated with baseline with the response to anti-TNF therapy is in line with previous transcriptomic evidence [9], in which poor responders showed a predominance of pathways associated with cell cycle regulation. This common observation indicates that in those patients where there is very active cell division in the synovial membrane, the TNF blockade will likely be inefficient. Our analysis did not allow us to differentiate however if this higher cell cycle activation is due to the immune or stromal cell compartments (or both). Further exploring this biological feature could provide a simple means to classify RA patients that will likely be resistant to the therapeutic effect of anti-TNF agents.

While our results show an association between a new pathogenic cell type and the response to therapy, this does not imply that this information is sufficient to be used to predict treatment response in RA. Aspects like effect size and technical and biological variation can hinder the capacity of a biological marker like TPH to predict anti-TNF response. Examples of the disparity between association and prediction have been shown previously for example in genetic polymorphisms associated with the disease risk with insufficient disease prediction [36]. To fill this gap, specific studies will need to be carried out to evaluate the predictor performance of this new biomarker in extended patient cohorts under different technologies and varying clinical features. In this regard, finding less invasive biomarkers like blood analytes that correlate TPH abundance would therefore be a useful means to simplify patient recruitment and stratification and test its utility to predict response to anti-TNF therapy.

Conclusions

In the present study, we have performed a multilevel analysis of the RA synovial membrane transcriptome in relation to the anti-TNF therapy response. Using a cell-type deconvolution approach, we have found that a newly characterized pathogenic T cell subtype in RA, peripheral T helper cell, is associated with the lack of response to TNF blocking. Our longitudinal analysis shows that good response at week 20 is preserved by a sustained low TPH cell number. At the pathway level, response to anti-TNF treatment is strongly linked to an efficacious clearance of myeloid leukocyte activation. At the single-gene level, two genes previously associated with RA pathology, CX3CL1 and PIK3CD, were found to be markers for the anti-TNF response, implicating variability in the immune cell recruitment and fibroblastic growth capacities of the synovial membrane. The results from this study indicate new biological features of anti-TNF response in RA and could be useful starting points to allow for patient stratification and for the development of new therapeutic approaches, alone or in combination with TNF blockers, to achieve a favorable response for a larger number of patients.

Availability of data and materials

The datasets generated and analyzed during the current study are available at the NCBI Gene Expression Omnibus database (https://www.ncbi.nlm.nih.gov/gds) with accession numbers GSE140036, GSE47726, and GSE15602.

Abbreviations

EULAR:

European Leagues Against Rheumatism

GO:

Gene ontology

mRNA:

Messenger RNA

RA:

Rheumatoid arthritis

RT-PCR:

Real-time polymerase chain reaction

TNF:

Tumor necrosis factor-alpha

TPH :

Peripheral helper T cells

References

  1. 1.

    Lipsky PE, van der Heijde DM, St Clair EW, Furst DE, Breedveld FC, Kalden JR, Smolen JS, Weisman M, Emery P, Feldmann M, et al. Infliximab and methotrexate in the treatment of rheumatoid arthritis. Anti-tumor necrosis factor trial in rheumatoid arthritis with concomitant therapy study group. N Engl J Med. 2000;343(22):1594–602.

    CAS  Article  Google Scholar 

  2. 2.

    Felson DT, Smolen JS, Wells G, Zhang B, van Tuyl LH, Funovits J, Aletaha D, Allaart CF, Bathon J, Bombardieri S, et al. American College of Rheumatology/European League Against Rheumatism provisional definition of remission in rheumatoid arthritis for clinical trials. Arthritis Rheum. 2011;63(3):573–86.

    Article  Google Scholar 

  3. 3.

    Strand V, Kimberly R, Isaacs JD. Biologic therapies in rheumatology: lessons learned, future directions. Nat Rev Drug Discov. 2007;6(1):75–92.

    CAS  Article  Google Scholar 

  4. 4.

    Sieberts SK, Zhu F, García-García J, Stahl E, Pratap A, Pandey G, Pappas D, Aguilar D, Anton B, Bonet J. Crowdsourced assessment of common genetic contribution to predicting anti-TNF treatment response in rheumatoid arthritis. Nat Commun. 2016;7:12460.

    CAS  Article  Google Scholar 

  5. 5.

    Smolen JS, Aletaha D, McInnes IB. Rheumatoid arthritis. Lancet. 2016;388:2023–38.

  6. 6.

    van der Pouw Kraan TC, Wijbrandts CA, van Baarsen LG, Rustenburg F, Baggen JM, Verweij CL, Tak PP. Responsiveness to anti-tumour necrosis factor α therapy is related to pre-treatment tissue inflammation levels in rheumatoid arthritis patients. Ann Rheum Dis. 2008;67(4):563–6.

    Article  Google Scholar 

  7. 7.

    Lindberg J, Wijbrandts CA, Van Baarsen LG, Nader G, Klareskog L, Catrina A, Thurlings R, Vervoordeldonk M, Lundeberg J, Tak PP. The gene expression profile in the synovium as a predictor of the clinical response to infliximab treatment in rheumatoid arthritis. PLoS One. 2010;5(6):e11310.

    Article  Google Scholar 

  8. 8.

    Lowe R, Shirley N, Bleackley M, Dolan S, Shafee T. Transcriptomics technologies. PLoS Comp Biol. 2017;13(5):e1005457.

    Article  Google Scholar 

  9. 9.

    Badot V, Galant C, Toukap AN, Theate I, Maudoux A-L, Van den Eynde BJ, Durez P, Houssiau FA, Lauwerys BR. Gene expression profiling in the synovium identifies a predictive signature of absence of response to adalimumab therapy in rheumatoid arthritis. Arthrit Res Ther. 2009;11(2):R57.

    Article  Google Scholar 

  10. 10.

    Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD. The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med. 2015;21(8):938.

    CAS  Article  Google Scholar 

  11. 11.

    Mizoguchi F, Slowikowski K, Wei K, Marshall JL, Rao DA, Chang SK, Blazar PE. Functionally distinct disease-associated fibroblast subsets in rheumatoid arthritis. Nat Commun. 2018;9(1):1–11.

  12. 12.

    Rao DA, Gurish MF, Marshall JL, Slowikowski K, Fonseka CY, Liu Y, Donlin LT, Henderson LA, Wei K, Mizoguchi F. Pathologically expanded peripheral T helper cell subset drives B cells in rheumatoid arthritis. Nature. 2017;542(7639):110.

    CAS  Article  Google Scholar 

  13. 13.

    Villeneuve E, Nam J, Emery P. 2010 ACR-EULAR classification criteria for rheumatoid arthritis. Rev Bras Reumatol. 2010;50(5):481–3.

    Article  Google Scholar 

  14. 14.

    Prevoo ML, van 't Hof MA, Kuper HH, van Leeuwen MA, van de Putte LB, van Riel PL: Modified disease activity scores that include twenty-eight-joint counts. Development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum 1995, 38(1):44–48.

  15. 15.

    Fransen J, van Riel PL. The disease activity score and the EULAR response criteria. Rheum Dis Clin N Am. 2009;35(4):745–57 vii-viii.

    Article  Google Scholar 

  16. 16.

    van de Sande MG, Gerlag DM, Lodde BM, van Baarsen LG, Alivernini S, Codullo V, Felea I, Vieira-Sousa E, Fearon U, Reece R, et al. Evaluating antirheumatic treatments using synovial biopsy: a recommendation for standardisation to be used in clinical trials. Ann Rheum Dis. 2011;70(3):423–7.

    Article  Google Scholar 

  17. 17.

    Rosenberg A, Fan H, Chiu YG, Bolce R, Tabechian D, Barrett R, Moorehead S, Baribaud F, Liu H, Peffer N. Divergent gene activation in peripheral blood and tissues of patients with rheumatoid arthritis, psoriatic arthritis and psoriasis following infliximab therapy. PLoS One. 2014;9(10):e110657.

    Article  Google Scholar 

  18. 18.

    Johnson WE, Li C, Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2007;8(1):118–27.

    Article  Google Scholar 

  19. 19.

    Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G. GO:: TermFinder—open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.

    CAS  Article  Google Scholar 

  20. 20.

    Bolen CR, Uduman M, Kleinstein SH. Cell subset prediction for blood genomic studies. BMC Bioinformatics. 2011;12(1):258.

    Article  Google Scholar 

  21. 21.

    Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci. 2005;102(43):15545–50.

    CAS  Article  Google Scholar 

  22. 22.

    Bankhead P, Loughrey MB, Fernández JA, Dombrowski Y, McArt DG, Dunne PD, McQuaid S, Gray RT, Murray LJ, Coleman HG. QuPath: open source software for digital pathology image analysis. Sci Rep. 2017;7(1):16878.

    Article  Google Scholar 

  23. 23.

    Shen P, Fillatreau S. Antibody-independent functions of B cells: a focus on cytokines. Nat Rev Immunol. 2015;15(7):441–51.

    CAS  Article  Google Scholar 

  24. 24.

    Wright HL, Moots RJ, Edwards SW. The multifactorial role of neutrophils in rheumatoid arthritis. Nat Rev Rheumatol. 2014;10(10):593.

    CAS  Article  Google Scholar 

  25. 25.

    Croft AP, Campos J, Jansen K, Turner JD, Marshall J, Attar M, Savary L, Wehmeyer C, Naylor AJ, Kemble S. Distinct fibroblast subsets drive inflammation and damage in arthritis. Nature. 2019;570(7760):246–51.

    CAS  Article  Google Scholar 

  26. 26.

    Nanki T, Imai T, Kawai S. Fractalkine/CX3CL1 in rheumatoid arthritis. Mod Rheumatol. 2017;27(3):392–7.

    CAS  Article  Google Scholar 

  27. 27.

    Auffray C, Fogg D, Garfa M, Elain G, Join-Lambert O, Kayal S, Sarnacki S, Cumano A, Lauvau G, Geissmann F. Monitoring of blood vessels and tissues by a population of monocytes with patrolling behavior. Science. 2007;317(5838):666–70.

    CAS  Article  Google Scholar 

  28. 28.

    Okkenhaug K, Vanhaesebroeck B. PI3K in lymphocyte development, differentiation and activation. Nat Rev Immunol. 2003;3(4):317–30.

    CAS  Article  Google Scholar 

  29. 29.

    Okkenhaug K, Bilancio A, Farjot G, Priddle H, Sancho S, Peskett E, Pearce W, Meek SE, Salpekar A, Waterfield MD, et al. Impaired B and T cell antigen receptor signaling in p110delta PI 3-kinase mutant mice. Science. 2002;297(5583):1031–4.

    CAS  PubMed  Google Scholar 

  30. 30.

    Puri KD, Doggett TA, Douangpanya J, Hou Y, Tino WT, Wilson T, Graf T, Clayton E, Turner M, Hayflick JS, et al. Mechanisms and implications of phosphoinositide 3-kinase delta in promoting neutrophil trafficking into inflamed tissue. Blood. 2004;103(9):3448–56.

    CAS  Article  Google Scholar 

  31. 31.

    Sriskantharajah S, Hamblin N, Worsley S, Calver AR, Hessel EM, Amour A. Targeting phosphoinositide 3-kinase delta for the treatment of respiratory diseases. Ann N Y Acad Sci. 2013;1280:35–9.

    CAS  Article  Google Scholar 

  32. 32.

    Bartok B, Boyle DL, Liu Y, Ren P, Ball ST, Bugbee WD, Rommel C, Firestein GS. PI3 kinase delta is a key regulator of synoviocyte function in rheumatoid arthritis. Am J Pathol. 2012;180(5):1906–16.

    CAS  Article  Google Scholar 

  33. 33.

    Whitehead MA, Bombardieri M, Pitzalis C, Vanhaesebroeck B. Isoform-selective induction of human p110delta PI3K expression by TNFalpha: identification of a new and inducible PIK3CD promoter. Biochem J. 2012;443(3):857–67.

    CAS  Article  Google Scholar 

  34. 34.

    Charbonneau M, Lavoie RR, Lauzier A, Harper K, McDonald PP, Dubois CM. Platelet-derived growth factor receptor activation promotes the prodestructive invadosome-forming phenotype of synoviocytes from patients with rheumatoid arthritis. J Immunol. 2016;196(8):3264–75.

    CAS  Article  Google Scholar 

  35. 35.

    Zhang, F, Wei K, Slowikowski K, Fonseka CY, Rao DA, Kelly S, Watts GF. Defining inflammatory cell states in rheumatoid arthritis joint synovial tissues by integrating single-cell transcriptomics and mass cytometry. Nat Immunol. 2019;20(7):928–42.

  36. 36.

    Marigorta UM, Rodríguez JA, Gibson G, Navarro A. Replicability and prediction: lessons and challenges from GWAS. Trends Genet. 2018;34(7):504–17.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We would like to thank Prof. Joachim Schultze and Svenja Debey-Pascher (University of Bonn, Germany) for their technical assistance in the transcriptome analysis, as well as Dr. Neus Prats and Dr. Mònica Aguilera (Histopathology Facility, Institute for Research in Biomedicine, Spain) for their generous assistance in the synovial immunofluorescence analyses.

Funding

This study was supported by the Fondo de Investigación Sanitaria del Institudo de Salud Carlos III, grant numbers: P11/1890, PI040710, and PI12/01362. The funding source had no role in the collection, analysis or interpretation of the data, or in the writing of the manuscript.

Author information

Affiliations

Authors

Contributions

AJ: study design, methodology, statistical analysis, and writing; GA: resources, investigation, data curation, and writing; RC: resources and data curation; RS: resources and investigation; JR: resources and investigation; SM: resources, methodology, and writing; JDC: study design, methodology, and writing. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Antonio Julià.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Hospital Clinic de Barcelona ethics committee, and all the patients signed the informed consent. This study was conducted in accordance with the Declaration of Helsinki principles.

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1 : Supplementary Table 1

. Marker genes of the three synovial membrane fibroblast types. Supplementary Table 2. Significantly overexpressed genes in PD-1hi T cells according to Rao et al. Supplementary Table 3. Preservation of gene expression differences between baseline and response time point (week 20). Supplementary Figure 1. Preservation of differential expression changes associated to anti-TNF response in the synovium from baseline. Supplementary Figure 2. RTPCR replication of CX3CL1 and PIK3CD genes. Supplementary Figure 3. Association of the three fibroblast subsets at baseline with anti-TNF response. Supplementary Figure 4. Cell types associated with anti-TNF response: differences at week 20. Supplementary Figure 5. QRTPCR replication of TPH biomarker genes CTLA4 and TIMELESS. Supplementary Figure 6. Correlation between PD1hiCXCR5-CD4+ T cell abundance estimation using the cell deconvolution approach with the percentage quantified using immunofluorescence.

Additional file 2 :: Supplementary Table 4.

Significant GO terms of genes differentially expressed at week 0 between responders and non-responders to anti-TNF therapy. Supplementary Table 5. Significant GO terms of genes differentially expressed at week 20 between responders and non-responders to anti-TNF therapy. Supplementary Table 6. Significant GO terms of genes differentially expressed between week 0 and week 20 in the responder group of patients.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Julià, A., Ávila, G., Celis, R. et al. Lower peripheral helper T cell levels in the synovium are associated with a better response to anti-TNF therapy in rheumatoid arthritis. Arthritis Res Ther 22, 196 (2020). https://doi.org/10.1186/s13075-020-02287-9

Download citation

Keywords

  • Rheumatoid arthritis
  • Anti-TNF therapy
  • Synovial membrane
  • Clinical response
  • Deconvolution
  • Peripheral T helper
  • Gene expression
  • Immunofluorescence