Pre-silencing of genes involved in the electron transport chain (ETC) pathway is associated with responsiveness to abatacept in rheumatoid arthritis
Arthritis Research & Therapy volume 19, Article number: 109 (2017)
In the current context of personalized medicine, one of the major challenges in the management of rheumatoid arthritis (RA) is to identify biomarkers that predict drug responsiveness. From the European APPRAISE trial, our main objective was to identify a gene expression profile associated with responsiveness to abatacept (ABA) + methotrexate (MTX) and to understand the involvement of this signature in the pathophysiology of RA.
Whole human genome microarrays (4 × 44 K) were performed from a first subset of 36 patients with RA. Data validation by quantitative reverse-transcription (qRT)-PCR was performed from a second independent subset of 32 patients with RA. Gene Ontology and WikiPathways database allowed us to highlight the specific biological mechanisms involved in predicting response to ABA/MTX.
From the first subset of 36 patients with RA, a combination including 87 transcripts allowed almost perfect separation between responders and non-responders to ABA/MTX. Next, the second subset of patients 32 with RA allowed validation by qRT-PCR of a minimal signature with only four genes. This latter signature categorized 81% of patients with RA with 75% sensitivity, 85% specificity and 85% negative predictive value. This combination showed a significant enrichment of genes involved in electron transport chain (ETC) pathways. Seven transcripts from ETC pathways (NDUFA6, NDUFA4, UQCRQ, ATP5J, COX7A2, COX7B, COX6A1) were significantly downregulated in responders versus non-responders to ABA/MTX. Moreover, dysregulation of these genes was independent of inflammation and was specific to ABA response.
Pre-silencing of ETC genes is associated with future response to ABA/MTX and might be a crucial key to susceptibility to ABA response.
Recent advances in our understanding of the pathophysiology of rheumatoid arthritis (RA) have led to the development of new biologic treatments designed to act against a precise therapeutic target: these include tumor necrosis factor-alpha (TNFα), interleukin-1-receptor antagonist (IL1-Ra, anakinra), cytotoxic T-lymphocyte-associated protein 4 (CTLA4-Ig or abatacept (ABA)), CD20 expressed on B cells (rituximab) and IL-6 receptor (tocilizumab) . All these biologic agents have proven efficacy in stopping joint inflammation and structural damage in association with the anchor drug methotrexate (MTX) . No one molecule has proven its clinical superiority over others in terms of efficacy [2,3,4,5,6]. In addition, no response to these treatments is obtained in approximately 30% of patients with RA, and the response to all these medications is highly variable from one patient to another. Nevertheless, prescription of biologics agents remains highly empirical [7,8,9]. Presently, we are still unable to predict the clinical efficacy of these treatments in a given patient because of the heterogeneity of RA and the subgroups of patients susceptible to responding better to one molecule than to another. Given the increasing number of available molecules, it is crucial to identify predictive markers in order to optimize drug prescription only to those patients susceptible to responding and thus avoiding side effects.
In the current context of personalized medicine, large-scale analysis of gene expression to predict drug response is a relevant and original approach, which has already shown its utility in cancer or kidney transplantation [10,11,12]. Gene expression profiling is clearly a powerful method for the identification of biomarkers and the development of personalized medicine in the field of rheumatology [13,14,15,16]. This approach has already allowed us to identify and validate two gene combinations able to predict response to infliximab or anakinra although in small cohorts of patients with RA [17, 18]. The proof of concept of gene expression profiles or signatures as a predictor of good response to drugs in RA was further confirmed by other teams using other biologic agents, but in cohorts that included very few patients with RA [19,20,21]. Thus, we conducted this present ancillary study as a follow-up to the APPRAISE trial [22, 23]. Our main objective was to identify and validate a gene expression profile associated with a good response to ABA/MTX, and to understand the involvement of this signature in RA pathophysiology.
Patients from the APPRAISE trial
A total of 68 patients with RA from the APPRAISE trial were enrolled in this ancillary study [22, 23]. The APPRAISE study assessed the capability of the composite power of Doppler and greyscale ultrasound score to measure the early effect and time course of response to treatment with ABA in biologic-naïve patients with active RA despite MTX therapy. APPRAISE (NCT00767325) initially including 104 patients with RA was a 24-week, Phase IIIb, open-label, multicenter, single-arm study conducted at 21 sites across Europe (Denmark, France, Germany, Hungary, Italy, Norway, Spain and the UK) [22, 23]. Eligible patients were ≥18 years of age, had American College of Rheumatology (ACR)-defined RA according to the 1987 classification criteria for at least 6 months , and had been on MTX (≥15 mg/week) for at least 3 months prior to baseline, with a stable MTX dose for at least 28 days before baseline (except in cases of intolerance to MTX). Patients were required to have active disease, defined by a baseline disease activity score in 28 joints (calculated with C-reactive protein (CRP)) (DAS28(CRP)) >3.2 or tender and ≥6 swollen joint counts, and CRP above the upper limit of normal. All patients received intravenous (IV) infusions of ABA at a weight-titered dose of 10 mg/kg at baseline (day 1), and at weeks 2, 4, 8, 12, 16, 20 and 24, in addition to stable doses of concomitant MTX (≥15 mg/week). MTX dose increases were not permitted, and dose decreases were allowed only in cases of intolerance. Oral corticosteroid use (stable dose of ≤10 mg prednisone/day) was permitted during the study. For this study, 5 ml of whole blood was collected in PAXgene RNA tube (PreAnalytiX, Qiagen) just before the first infusion and 6 months later and was stored at –80 °C until use.
Clinical evaluation and response to ABA/MTX
Several clinical characteristics were collected at baseline and 6 months later: age, gender, disease duration and MTX and corticosteroid doses. Disease activity was evaluated at all assessment visits (baseline, weeks 1, 2, 4, 6, 8, 12, 16, 20 and 24) using the DAS28(CRP) calculated from 28 tender joints, 28 swollen joints, CRP and patient global assessment (visual analog scale (VAS); 0–10 scale).
The response to ABA/MTX was evaluated at 6 months using European League Against Rheumatism (EULAR) response criteria based on DAS28(CRP) . Since we were looking for biomarkers associated with good response to ABA/MTX, patients were categorized according to their EULAR response as responders (R) (n = 36) or non-responders (NR) (comprising moderate responders (n = 25) and no responders (n = 7)) .
Total RNAs from whole blood were extracted with PAXgene blood RNA kit according to the manufacturer’s recommendations (Qiagen PreAnalytiX GmbH, Courtaboeuf, France) and stored at –80 °C until use. Total RNA from 10 healthy donors (5 women and 5 men) was pooled and used as an internal standard reference (control pool). The quality and quantity of isolated mRNAs were assessed using the 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and the Nanodrop device (Thermo Scientific, Wilmington, USA). Only RNA samples with a minimal RNA integrity number of 7 were used for subsequent experiments.
Whole human genomic DNA microarrays were used to analyze two-colored gene expression profiling (4 × 44 K Whole Human Genome, Agilent Technologies, Les Ulis, France). Each RNA sample from patients with RA was labeled by Cyanine-5 and co-hybridized with a Cyanine-3 labeled RNA control pool according to the manufacturer’s instructions (Low Input QuickAmp Labeling Kit, Agilent Technologies, Les Ulis, France). Briefly, 100 ng of RNAs were labeled with cyanine-5 CTP (patients with RA) or cyanine-3 CTP (control pool). After hybridization reaction using a hybridization kit (Agilent Technologies, Les Ulis, France) co-hybridization was performed at 65 °C for 17 hours. After wash steps, the microarrays were scanned with a 5-μM pixel size using the DNA Microarray Scanner GB (Agilent Technologies, Les Ulis, France). Image analysis and extraction of raw and normalized signal intensities (lowess) were performed using Feature Extraction Software 10.5.1.1 (Agilent Technologies). The data were in agreement with the guidelines for minimum information about a microarray experiment and were deposited in the database of the National Center for Biotechnology Information Gene Expression Omnibus (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE68215). The data are accessible [GEO:GSE68215]. Non-uniform spots and saturated spots or spots with intensities below the background were not taken into account. Only spots that passed these quality controls on 100% of arrays were selected for further analysis. Hierarchical clustering was performed using the Pearson coefficient metric and complete linkage to build the transcripts and sample dendrograms.
Quantitative reverse transcription-PCR (qRT-PCR)
cDNA was synthesized from 1-μg RNA samples using random primers and M-MLV enzyme (Invitrogen™, Carlsbad, USA). qRT-PCR was performed using a Lightcycler as instructed by the manufacturer (Roche™, Meylan, France). qRT-PCR reactions were performed for each sample in duplicate using SYBR-Green (Roche™, Meylan, France) and values were normalized using the geometric mean of three control genes (18S, ACTB, B2M) defined by the geNorm algorithm . Sequences of primers (Eurogentec™, Fremont, USA) used for qRT-PCR were: BLOC1S1 forward, 5’-AAGCAGACAGGCCAGTGGAT-3’; BLOC1S1 reverse, 5’-CAGTGCAGTGGCAATGGTG-3’; RNASE3 forward, 5’-CAGGAGCCACAGCTCAGAGA-3’; RNASE3 reverse, 5’-GAGCCCTCCACACCCATAAG-3’; COX6A1 forward, 5’-CCACTTCCAACTGGCTACGA-3’; COX6A1 reverse, 5’-AAGCAAAGGGATGGGAGACC-3’; PTRH2 forward, 5’-GCTGTTGGAGTTGCTTGTGG-3’; PTRH2 reverse, 5’-AGGCTGAAACAGCAGCATGA-3’; 18S forward, 5’-GTGGAGCGATTTGTCTGGTT-3’; 18S reverse, 5’-CGCTGAGCCAGTCAGTGTAG-3’; ACTB forward, 5’-CTGGAACGGTGAAGGTGACA-3’; ACTB reverse, 5’-AAGGGACTTCCTGTAACAATGCA-3’; B2M forward, 5’-TGCTGTCTCCATGTTTGATGTATCT-3’; B2M reverse, 5’-TCTCTGCTCCCCACCTCTAAGT-3’.
Statistical and functional analysis
Comparisons of clinical and biological data between R and NR were performed at baseline using Student’s t test for continuous variables. Comparisons of R or NR before and after treatment were performed using the paired t test. Identification of clinical parameters able to predict good response to ABA/MTX was performed in two different multivariate analyses: (1) a logistic regression model with variable selection using Bayesian information criterion (BIC) and (2) linear discriminant analysis (LDA), which estimates a coefficient for each variable.
Data from transcriptomic analysis were analyzed using GeneSpring GX V.13.0 (Agilent Technologies, Les Ulis, France). The normality of log 2 ratio of gene expression was evaluated using the Shapiro–Wilk statistical test. The unpaired Student’s t test (p value < 0.05), with the Benjamini–Hochberg correction to check the false discovery rate (FDR), was used to determine the statistical significance of differences in gene expression levels between R and NR. Gene Ontology (GO) analysis was used to investigate the biological processes, molecular function or cellular localization enriched in the transcripts list, showing a significant fluctuation in gene expression between R and NR. The p value was computed by standard hypergeometric distribution. The GeneSpring Single Experiment Analysis (SEA) bio-informatics tool was used for computational analysis to identify potential curated canonical pathways with setting parameters (reactome and GenMAPP for pathway source), which are enriched in the differentially expressed transcripts list, using the WikiPathways database (http://www.wikipathways.org/index.php/Pathway:WP111). The significance of the association between the genes and the pathways was measured by Fisher’s exact test.
Characteristics of patients with RA and their response to ABA/MTX
Of the 104 patients with RA included in the original APPRAISE trial, clinical and biological data for subsequent analysis were available for 91 patients: 68 of these patients were recruited to this present ancillary study based on the quality of the RNA samples (Fig. 1). These 68 patients with RA were split into two subsets at random: subset 1 (n = 36) to identify clinical or biological (including transcripts) markers associated with response to ABA/MTX and subset 2 (n = 32) to validate a gene expression profile able to predict good response to ABA/MTX (Fig. 1). After 6 months of treatment, patients with RA were categorized according to their EULAR response as either responders (R: n = 17 and 19, respectively in subsets 1 and 2) or non-responders (NR: n = 19 and 13, respectively in subsets 1 and 2). Table 1 provides demographic and clinical information for these 68 patients with RA at baseline and after 6 months of treatment.
Clinical and biological variables were collected before treatment and 6 months later. Values are expressed as mean ± standard error of the mean unless stated otherwise. R responders, NR non-responders, F female, M male, VAS visual analog scale, DAS Disease Activity Score in 28 joints, CRP C-reactive protein. Response to abatacept/methotrexate was assessed by the DAS28 calculated with CRP at 6 months of treatment. Patients were categorized as indicated in “Methods”. P values were determined by Student’s independent samples t test or paired t test as appropriate: *p < 0.05; **p < 0.01 for comparisons between responders and non-responders at baseline; $$$ p < 0.001 for comparison of values at baseline and 6 months in responders; ££ p < 0.01; £££ p < 0.001 for comparison of values at baseline and 6 months in non-responders. All other comparisons were not significant.
The baseline characteristics of patients with RA from subsets 1 and 2 were comparable for all variables, suggesting absence of bias in randomization (Table 1). Whatever the subset, tender joint count (TJC), swollen joint count, global assessment of disease measured by the patient with the VAS, CRP and DAS28(CRP) improved significantly after 6 months of treatment in both R and in NR. However, the range of improvement was significantly higher in R than in NR for each parameter except the global VAS and CRP (Table 1 and Additional file 1). In each subset, there were statistical differences between R and NR at baseline in clinical and biological variables. Indeed, disease duration was longer in NR than in R (p < 0.01) (Table 1). DAS28, patient assessment of disease and CRP were higher in NR than in R (p < 0.05). Moreover, the sex ratio was significantly different between R and NR (p = 0.047). There was better response in patients with RA of short disease duration and moderate disease activity. These observations raise the question of using clinical and biological variables as predictors of response.
Clinical and biological variables are not associated with response to ABA/MTX
Clinical and/or biological variables as predictors of response were identified by two different multivariate statistical approaches, logistic regression and linear discriminant analysis. Logistic regression based on more pertinent variables as MTX, TJC, disease duration and CRP (identified in subset 1) did not allow an accurate prediction of response to ABA/MTX in subset 2 (sensitivity = 63%; specificity = 46%) (Fig. 2a). Similarly, linear discriminant analysis, which balanced biological and clinical parameters, did not allow accurate prediction of response to ABA/MTX in subset 2 (sensitivity = 68%; specificity = 46%) (Fig. 2b). As a result of these shortcomings, we used a microarray approach to identify biomarkers associated with response to ABA/MTX.
Gene expression profiling associated with response to ABA/MTX
With subset 1, the cRNA from 17 R and 19 NR were co-hybridized with a cRNA internal reference from 10 healthy subjects on whole human genome microarrays, by two-color technology. After elimination of spikes and flagged probes, 19,620 probes were detected in all samples. From these 19,620 transcripts, we identified 87 transcripts with statistically significant variation between R and NR according to the t test with correction for multiple testing (FDR, Benjamini–Hochberg correction) (Fig. 3a). These transcripts are listed in Additional file 2 (69 transcripts are referenced with a Ref Seq accession number while 18 probes were unknown). Since the sex ratio in R and in NR was significantly different, this maldistribution between male and female patients did not introduce bias into our study (data not shown) Finally, we performed hierarchical clustering of the 36 patients from subset 1 based on the levels of the 87 transcripts indicated previously, resulting in almost perfect separation of R and NR into two major clusters. Indeed, only one R patient was misclassified in the patient’s dendrogram (Fig. 3b).
We wanted to confirm that a combination of the aforementioned transcript levels was associated with response to ABA/MTX. For this purpose, we used linear discriminant analysis to select a minimal combination of transcripts from “identification” subset 1 that were able to correctly classify patients with RA from “validation” subset 2. Among these 87 mRNA, four transcripts (BLOC1S1, RNASE3, COX6A1, PTRH2) were retained by linear discriminant analysis although they did not have the smallest p value (Table 2). As shown in Fig. 4, hierarchical clustering of the 32 patients from subset 2, based on the levels of these four transcripts measured by qRT-PCR, resulted in two major clusters of R versus NR, with 26 well-categorized patients. This procedure identified six misclassified patients and indicated that this set of four transcripts provides 75% sensitivity, 85% specificity, and 85% negative predictive value and 75% positive predictive value for identification of future R and NR to ABA/MTX. To ascertain if this signature was specific to ABA, we used this combination of 87 mRNA associated with drug response in an independent cohort of patients with RA, with the same level of disease activity, treated by TNFα blocking agents given as etanercept or adalimumab (Additional file 3), both associated with MTX. This combination was unable to correctly classify R and NR to TNFα blocking agents whatever the molecular structure (fusion protein or monoclonal antibody directed against TNFα), suggesting specificity of this signature for ABA/MTX (Additional file 3). This combination of 87 transcripts also raised the question of understanding the biological significance of these genes in response to ABA/MTX in the pathophysiology of RA.
Pre-silencing of electron transport chain pathway associated with response to ABA/MTX
Different approaches were used to understand the involvement of this 87 mRNA signature in biological processes. These 87 mRNA were submitted to GO analysis leading us to discover enrichment (from 6 to 49 times more) of 18 GO classes in these 87 mRNA compared to the whole human genome (Table 2). Among these 87 transcripts, 11 were associated with these 18 GO classes (Additional file 4).
All these 18 GO classes, based on the 87 mRNA, were relative to the mitochondrial respiratory chain located in the inner membrane of mitochondria (corrected p value <0.005). For instance, 6/87 mRNA of genes linked to the electron transport chain (ETC; GO: 0022900) were included in our signature while 101 genes from the whole human genome are known to be relative to this same GO class (corrected p value = 0.001). The fold enrichment in this GO term is × 16 (Table 2). Moreover, these 87 mRNA were submitted to WikiPathways using the SEA tool from GeneSpring software leading us to identify enrichment of four signaling pathways: Electron Transport Chain WP111_41171 (version of 1 March 2011; p value = 1.6.10-9), Oxidative Phosphorylation WP623_45305 (version of 7 October 2011; p value = 2.5.10-4), Proteasome Degradation WP183_45274 (version of 7 October 2011; p value = 0.008) and TSH signaling pathway WP2032_44635 (version of 22 September 2011 ; p value = 0.008). GO analysis and SEA were in agreement as the ETC pathway localized in mitochondria was found by using these two tools.
As shown in Fig. 5a, the ETC pathway includes 104 proteins split up into five complexes of ETC, embedded in the inner membrane of mitochondria. Among the 87 mRNA of our signature, 7 transcripts from the ETC pathway (NDUFA6, NDUFA4, UQCRQ, ATP5J, COX7A2, COX7B, COX6A1) covering four out of five complexes of this pathway, were significantly (p value <0.05) downregulated in 17 R versus 19 NR to ABA/MTX (Fig. 5a and b). This differential gene expression profile including these seven genes suggested reduced expression of the mitochondrial respiratory chain pathway in R before ABA administration (Fig. 5b). In addition, taking into account gene expression profiling from the whole ETC pathway (n = 104 genes), principal component analysis correctly separated R from NR before treatment, even if all the transcripts from this pathway were not significantly dysregulated between R and NR (Fig. 5c). Nevertheless, this ETC combination including these 104 transcripts did not better separate R and NR (73% specificity) than did the minimal combination of four transcripts (85% specificity) (Fig. 4).
The dysregulation of these transcripts could be linked to the more elevated CRP rate observed in NR compared to R at baseline (Table 1). Correlations between each significant transcript level from the ETC pathway (NDUFA6, NDUFA4, UQCRQ, ATP5J, COX7A2, COX7B, COX6A1) and CRP or DAS28(CRP) were performed to discard a CRP or disease activity effect on the gene expression levels from this signature linked to a supposedly more pronounced mitochondrial activity in the inflammatory state. A weak positive coefficient of correlation (r) and coefficient of determination (r 2) between CRP (0.25 < r < 0.49; 0.06 < r 2 < 0.24) or DAS28(CRP) (0.32 < r < 0.43; 0.10 < r 2 < 0.18) and gene expression level was observed for the majority of genes except ATP5J. In subset 2, weak positive correlation between CRP or DAS28(CRP) and gene expression data was also found for NDUFA4, COX6A1, COX7A2 and COX7B (data not shown). All these data led us to discard the impact of CRP or disease activity on the dysregulation of the ETC pathway between R and NR at baseline (Additional file 5).
Finally, we measured gene expression fluctuation in an independent transcriptomic analysis between baseline and 6 months of treatment by ABA/MTX in R and NR. Only five genes (COX7B, COX11, UQCRC2, NDUFU3, NDUFS1) among the 104 genes involved in the ETC pathway were significantly upregulated in R after 6 months of treatment, while these genes were invariant in NR to ABA/MTX (data not shown). This observation suggests that ABA increased gene expression levels of some genes from the ETC pathway only in R and not in NR independently of disease activity, which decreased whatever the response to ABA/MTX. Overall, these results suggest that these gene levels were not related to disease but rather they are genes of susceptibility to ABA and that silent expression of the ETC signature is associated more with good response to ABA/MTX.
In this ancillary study, we identified 87 transcripts with relative abundances, which were able to separate R and NR to ABA/MTX in 36 patients with RA before treatment (Fig. 3 and Additional file 2). Next, for the first time we validated a minimal combination of four transcripts associated with response to ABA/MTX in an independent subset of 32 patients with RA (Fig. 4). According to GO and WikiPathways analysis, we identified enrichment of 7 genes among the aforementioned 87 belonging to the ETC pathway, suggesting a specific signature of response (Fig. 5). Indeed, patients with RA with pre-silencing of the ETC pathway before treatment were those most susceptible to responding to ABA/MTX 6 months later.
In the current context of personalized medicine for the management of RA, the optimization of drug prescription is crucial. Hence, drug prescription should be tailored exclusively only to patients susceptible to responding to the drug. The identification of biomarkers to predict of drug response is paramount. To date, there are only two biomarkers associated with good response to ABA/MTX in RA, anti-CCP positivity and low baseline number of CD8+CD28– T cells [27,28,29,30]. Furthermore, there is no clinical and/or biological biomarker used in routine practice able to predict response to ABA/MTX prior to treatment initiation. Indeed, in our study, two statistical approaches to multivariate analysis were unable to identify predictive clinical or biological parameters of response to ABA/MTX, even if disease duration, CRP, global patient assessment of disease and DAS28(CRP) were significantly different between R and NR at baseline (Table 1 and Fig. 2). It would be interesting to establish a relationship between response to ABA and anti-citrullinated protein antibodies (ACPA) and/or rheumatoid factor status, given that ACPA positivity is associated with better response to ABA in RA. Unfortunately, we were unable to include ACPA or rheumatoid factor status in our analysis because these data were not collected during the APPRAISE trial and were not available. The difficulties we faced in predicting drug response with clinical and/or biological variables led us to use a more global functional genomic approach without a priori having to identify a signature associated with good response to ABA/MTX.
The transcriptomic approach based on the whole genome allowed us to identify signatures able to separate R from NR to several drugs (infliximab, tocilizumab, rituximab) used in RA [17,18,19,20,21, 31]. For ABA, one study previously measured the type-I interferon-regulated transcripts from peripheral blood mononuclear cells in patients with RA, independently of response to ABA . Another study identified gene sets associated with remission according to clinical disease activity index (CDAI) in ABA-treated patients with RA . In our study, even if there were differences at baseline in disease duration, CRP, patient assessment of disease and DAS28(CRP) between R and NR, which might have influenced gene expression, we identified and validated a gene signature associated with response to ABA/MTX. This signature included 87 genes variously involved in the ETC, proteasome, interferon and RNA processes, etc. This combination of 87 genes is also specific to ABA as it was not able to predict TNFα blocking agent response. Next, we validated this signature with an independent subset of 32 patients with RA by means of a minimal combination of four genes with expression levels measured by qRT-PCR, which are more easily usable in routine practice than microarrays. This signature is the optimized combination of genes associated with drug response with good accuracy, because it correctly predicted the future response in 81% of patients with RA (75% sensitivity, 85% specificity, 85% negative predictive value and 75% positive predictive value) prior to treatment (Fig. 4). Each gene taken separately was unable to predict response to ABA/MTX (data not shown) but all genes together were associated with response to ABA/MTX with good accuracy. Nevertheless, further validation studies in independent cohorts are essential before considering this signature as a predictive biomarker and of use in clinical practice.
Of the four genes in the minimal combination, the RNAseIII gene codes for the RNAseIII enzyme that specifically cleaves double-stranded RNA and is involved in the processing of ribosomal RNA precursors of some mRNAs . Biogenesis of lysosomal organelles complex-1, subunit 1 (BLOC1S1) codes for the protein BLOC1S1, also known as GCN5L1, and is an essential component of the mitochondrial acetyltransferase machinery and modulates mitochondrial respiration via acetylation of ETC proteins . COX6A1 codes for the mitochondrial protein cytochrome c oxidase subunit 6A1 (COX6A1) located in complex IV. This is the last enzyme in the mitochondrial ETC that drives ATP synthesis . PTRH2 codes for the peptidyl-tRNA hydrolase 2, which is a mitochondrial protein released from mitochondria to the cytoplasm during apoptosis.
Out of these 4 transcripts, 3 (BLOC1S1, COX6A1 and PTRH2) and 13 probes out of 87 code for proteins located in mitochondria. In addition, some of them are involved in the ETC pathway, suggesting implication of mitochondrial metabolism in response to ABA/MTX. Moreover, GO analysis and SEA are in agreement as these two analyses point to the involvement of ETC in the response to ABA/MTX. Interestingly, we found that seven transcripts with levels that are significantly lower in R than in NR at baseline were similarly regulated in patients in remission and patients not in remission defined by clinical disease activity index (CDAI) in a recent study . The ETC is a series of five complexes anchored to the inner membrane of mitochondria that transfers electrons via redox reactions, which drives ATP synthesis, generating reactive oxygen species (ROS) and subsequent oxidative stress . Redox balance in mitochondria is a critical component in T cell activation and proliferation . The production of ROS by the ETC complex III leads to production of large amounts of ATP to enhance activity of proliferating T cells after TCR cross-linking [39,40,41,42]. In our study, ETC genes were downregulated in R compared to NR, suggesting potentially low ROS production and potentially less T cell activation before treatment. Also, a previous study suggested it would be interesting to determine whether ROS production in T cells might be a predictor of clinical response to ABA . Our results suggest the possible involvement of ROS in the pathophysiological processes of ABA response, but further functional studies are necessary to confirm this hypothesis.
This 87 mRNA signature also included RASSF5, which was significantly upregulated in R compared to NR. RASSF5, also known as RAPL, is the effector of Rap1, which plays a central role in T cell response through TCR and co-stimulation signals. Indeed, a model was proposed in which inactivation of Rap1 plays a central role in establishing oxidative stress and can influence T cell response in RA [43, 44]. These data suggest less oxidative stress in future responders to ABA while NR present with high expression of genes from the ETC pathway, showing oxidative stress. A reduction in the expression level of ETC genes seems to increase the sensitivity of patients with RA to ABA/MTX. This model has already been highlighted in esophageal adenocarcinoma and colorectal cancer treated by chemotherapy [45, 46]. Indeed, ATP5J and COX7A2 included in our combination were also found to be downregulated and associated with response to chemotherapy, respectively in colorectal cancer and esophageal adenocarcinoma [45, 46]. As in cancer, pre-treatment conditions targeting the mitochondrial metabolism might be a determinant of susceptibility to ABA/MTX .
After 6 months of treatment with ABA, we showed that five genes (COX7B, COX11, UQCRC2, NDUFU3, NDUFS1) involved in the ETC pathway were significantly upregulated in R, while these genes were invariant in NR to ABA/MTX. So, whereas ABA does not seem to affect the oxidative state in NR, it seems to modulate oxidative stress in R, as genes from the ETC pathway were upregulated under ABA. This drug restores the expression of genes involved in redox balance.
In summary, we identified and validated a gene expression signature associated with good response to ABA/MTX. This signature involves genes from the ETC pathway. The pre-silencing of these genes could be a specific determinant of susceptibility to MTX + ABA. The downregulation of these genes reflected the redox imbalance observed in R. The intracellular redox balance is crucial for the antigenic response of T cells. After 6 months of treatment, ABA significantly upregulated some ETC genes in R probably leading to a slight increase in ROS, restoring the redox balance and improving T cell response. Further studies, including the investigation of ROS levels in R and NR to ABA/MTX, are necessary for large-scale validation of this signature before evaluating its pertinence in clinical practice.
Anti-citrullinated protein antibodies
American College of Rheumatology
ATP Synthase H+ transporting Mitochondrial Fo Complex
Biogenesis of lysosomal organelles complex-1 subunit 1
Clinical Disease Activity Index
Cytochrome C oxidase subunit 1
Cytochrome c oxidase subunit 6A1
Cytochrome C oxidase subunit VIIa polypeptide 2
Cytochrome C oxidase subunit VIIb
Complementary ribonucleic acid
Disease Activity Score
Electron transport chain
European League Against Rheumatism
False discovery rate
General control of amino acid synthesis 5 like 1
Interleukin 1 receptor antagonist
Messenger ribonucleic acid
NADH Dehydrogenase ubiquinone 1 alpha subcomplex 4
NADH dehydrogenase ubiquinone 1 alpha subcomplex 6
NADH ubiquinone oxidoreductase core subunit S1
NADH dehydrogenase-ubiquinone Fe-S subunit 3
Power Doppler grayscale ultrasound
Proteasome subunit alpha type 6
Proteasome subunit beta type-6
Peptidyl-tRNA hydrolase 2
Quantitative reverse-transcription polymerase chain reaction
Ras-related protein 1
Ras association domain family member 5
Reactive oxygen species
Single Experiment Analysis
tender joint count
Tumor necrosis factor alpha
Ubiquitin like 5
Ubiquinol cytochrome C reductase complex core protein 2
Ubiquinol-cytochrome C reductase complex III subunit VII
Visual analog scale
Smolen JS, Aletaha D. Rheumatoid arthritis therapy reappraisal: strategies, opportunities and challenges. Nat Rev Rheumatol. 2015;11(5):276–89.
Nam JL, Ramiro S, Gaujoux-Viala C, Takase K, Leon-Garcia M, Emery P, et al. Efficacy of biological disease-modifying antirheumatic drugs: a systematic literature review informing the 2013 update of the EULAR recommendations for the management of rheumatoid arthritis. Ann Rheum Dis. 2014;73(3):516–28.
Orme ME, Macgilchrist KS, Mitchell S, Spurden D, Bird A. Systematic review and network meta-analysis of combination and monotherapy treatments in disease-modifying antirheumatic drug-experienced patients with rheumatoid arthritis: analysis of American College of Rheumatology criteria scores 20, 50, and 70. Biologics. 2012;6:429–64.
Ramiro S, Gaujoux-Viala C, Nam JL, Smolen JS, Buch M, Gossec L, et al. Safety of synthetic and biological DMARDs: a systematic literature review informing the 2013 update of the EULAR recommendations for management of rheumatoid arthritis. Ann Rheum Dis. 2014;73(3):529–35.
Singh JA, Christensen R, Wells GA, Suarez-Almazor ME, Buchbinder R, Lopez-Olivo MA, et al. A network meta-analysis of randomized controlled trials of biologics for rheumatoid arthritis: a Cochrane overview. CMAJ. 2009;181(11):787–96.
Singh JA, Wells GA, Christensen R, Tanjong Ghogomu E, Maxwell L, Macdonald JK, et al. Adverse effects of biologics: a network meta-analysis and Cochrane overview. Cochrane Database Syst Rev. 2011;2:CD008794.
Gaujoux-Viala C, Gossec L, Cantagrel A, Dougados M, Fautrel B, Mariette X, et al. Recommendations of the French Society for Rheumatology for managing rheumatoid arthritis. Joint Bone Spine. 2014;81(4):287–97.
Singh JA, Furst DE, Bharat A, Curtis JR, Kavanaugh AF, Kremer JM, et al. 2012 update of the 2008 American College of Rheumatology recommendations for the use of disease-modifying antirheumatic drugs and biologic agents in the treatment of rheumatoid arthritis. Arthritis Care Res (Hoboken). 2012;64(5):625–39.
Smolen JS, Breedveld FC, Burmester GR, Bykerk V, Dougados M, Emery P, et al. Treating rheumatoid arthritis to target: 2014 update of the recommendations of an international task force. Ann Rheum Dis. 2015;75(1):3–15.
Lenz G, Wright G, Dave SS, Xiao W, Powell J, Zhao H, et al. Stromal gene signatures in large-B-cell lymphomas. N Engl J Med. 2008;359(22):2313–23.
Suthanthiran M, Schwartz JE, Ding R, Abecassis M, Dadhania D, Samstein B, et al. Urinary-cell mRNA profile and acute cellular rejection in kidney allografts. N Engl J Med. 2013;369(1):20–31.
Markopoulos C, van de Velde C, Zarca D, Ozmen V, Masetti R. Clinical evidence supporting genomic tests in early breast cancer: do all genomic tests provide the same information? Eur J Surg Oncol. 2017;43:909–20.
Burska AN, Roget K, Blits M, Soto Gomez L, van de Loo F, Hazelwood LD, et al. Gene expression analysis in RA: towards personalized medicine. Pharmacogenomics J. 2014;14(2):93–106.
de Jong TD, Vosslamber S, Verweij CL. Moving towards personalized medicine in rheumatoid arthritis. Arthritis Res Ther. 2014;16(3):110.
Haupl T, Stuhlmuller B, Grutzkau A, Radbruch A, Burmester GR. Does gene expression analysis inform us in rheumatoid arthritis? Ann Rheum Dis. 2010;69(1):i37–42.
Folkersen L, Brynedal B, Diaz-Gallo LM, Ramskold D, Shchetynsky K, Westerlind H, et al. Integration of known DNA, RNA and protein biomarkers provides prediction of anti-TNF response in rheumatoid arthritis: results from the COMBINE study. Mol Med. 2016;22:322–8.
Bansard C, Lequerre T, Derambure C, Vittecoq O, Hiron M, Daragon A, et al. Gene profiling predicts rheumatoid arthritis responsiveness to IL-1Ra (anakinra). Rheumatology (Oxford). 2010;50(2):283–92.
Lequerre T, Gauthier-Jauneau AC, Bansard C, Derambure C, Hiron M, Vittecoq O, et al. Gene profiling in white blood cells predicts infliximab responsiveness in rheumatoid arthritis. Arthritis Res Ther. 2006;8(4):R105.
Krintel SB, Dehlendorff C, Hetland ML, Horslev-Petersen K, Andersen KK, Junker P, et al. Prediction of treatment response to adalimumab: a double-blind placebo-controlled study of circulating microRNA in patients with early rheumatoid arthritis. Pharmacogenomics J. 2015;16(2):141–6.
Sanayama Y, Ikeda K, Saito Y, Kagami S, Yamagata M, Furuta S, et al. Prediction of therapeutic responses to tocilizumab in patients with rheumatoid arthritis: biomarkers identified by analysis of gene expression in peripheral blood mononuclear cells using genome-wide DNA microarray. Arthritis Rheum. 2014;66(6):1421–31.
Sellam J, Marion-Thore S, Dumont F, Jacques S, Garchon HJ, Rouanet S, et al. Use of whole-blood transcriptomic profiling to highlight several pathophysiologic pathways associated with response to rituximab in patients with rheumatoid arthritis: data from a randomized, controlled, open-label trial. Arthritis Rheum. 2014;66(8):2015–25.
D'Agostino MA, Boers M, Wakefield RJ, Berner Hammer H, Vittecoq O, Filippou G, et al. Exploring a new ultrasound score as a clinical predictive tool in patients with rheumatoid arthritis starting abatacept: results from the APPRAISE study. RMD Open. 2016;2(1):e000237.
D'Agostino MA, Wakefield RJ, Berner-Hammer H, Vittecoq O, Filippou G, Balint P, et al. Value of ultrasonography as a marker of early response to abatacept in patients with rheumatoid arthritis and an inadequate response to methotrexate: results from the APPRAISE study. Ann Rheum Dis. 2016;75(10):1763–9.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.
van der Heijde DM, van’t Hof MA, Van Riel PL, Theunisse LA, Lubberts EW, van Leeuwen MA, et al. Judging disease activity in clinical practice in rheumatoid arthritis: first step in the development of a disease activity score. Ann Rheum Dis. 1990;49(11):916–20.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3(7):RESEARCH0034.
Gottenberg JE, Neto D, Gomez-Reino J, Iannone F, Lie E, Canhão H, et al. Positivity for rheumatoid factor and anti-cyclic citrullinated peptide is associated with a better drug retention of abatacept: data from a pan European analysis of RA registries. Ann Rheum Dis. 2014;73(supp 2):502–3.
Gottenberg JE, Ravaud P, Cantagrel A, Combe B, Flipo RM, Schaeverbeke T, et al. Positivity for anti-cyclic citrullinated peptide is associated with a better response to abatacept: data from the ‘Orencia and Rheumatoid Arthritis’ registry. Ann Rheum Dis. 2012;71(11):1815–9.
Nusslein HG, Alten R, Galeazzi M, Lorenz HM, Nurmohamed MT, Bensen WG, et al. Prognostic factors for abatacept retention in patients who received at least one prior biologic agent: an interim analysis from the observational, prospective ACTION study. BMC Musculoskelet Disord. 2015;16:176.
Scarsi M, Ziglioli T, Airo P. Baseline numbers of circulating CD28-negative T cells may predict clinical response to abatacept in patients with rheumatoid arthritis. J Rheumatol. 2011;38(10):2105–11.
Toonen EJ, Gilissen C, Franke B, Kievit W, Eijsbouts AM, den Broeder AA, et al. Validation study of existing gene expression signatures for anti-TNF treatment in patients with rheumatoid arthritis. PLoS One. 2012;7(3):e33199.
Serikawa KA, Jacobsen S, Lundsgaard D, Fox BA, Hummelshoj L, Poulsen LK, et al. Detection of gene expression signatures related to underlying disease and treatment in rheumatoid arthritis patients. Mod Rheumatol. 2012;23(4):729–40.
Nakamura S, Suzuki K, Iijima H, Hata Y, Lim CR, Ishizawa Y, et al. Identification of baseline gene expression signatures predicting therapeutic responses to three biologic agents in rheumatoid arthritis: a retrospective observational study. Arthritis Res Ther. 2016;18:159.
Nicholson AW. Ribonuclease III, mechanisms of double-stranded RNA cleavage. Wiley Interdiscip Rev RNA. 2013;5(1):31–48.
Scott I, Webster BR, Li JH, Sack MN. Identification of a molecular component of the mitochondrial acetyltransferase programme: a novel role for GCN5L1. Biochem J. 2012;443(3):655–61.
Fornuskova D, Stiburek L, Wenchich L, Vinsova K, Hansikova H, Zeman J. Novel insights into the assembly and function of human nuclear-encoded cytochrome c oxidase subunits 4, 5a, 6a, 7a and 7b. Biochem J. 2010;428(3):363–74.
Chang HY, Huang TC, Chen NN, Huang HC, Juan HF. Combination therapy targeting ectopic ATP synthase and 26S proteasome induces ER stress in breast cancer cells. Cell Death Dis. 2014;5:e1540.
Leavy O. T cells: Mitochondria and T cell activation. Nat Rev Immunol. 2013;13(4):224.
Sena LA, Li S, Jairaman A, Prakriya M, Ezponda T, Hildeman DA, et al. Mitochondria are required for antigen-specific T cell activation through reactive oxygen species signaling. Immunity. 2013;38(2):225–36.
Devadas S, Zaritskaya L, Rhee SG, Oberley L, Williams MS. Discrete generation of superoxide and hydrogen peroxide by T cell receptor stimulation: selective regulation of mitogen-activated protein kinase activation and fas ligand expression. J Exp Med. 2002;195(1):59–70.
Ledderose C, Bao Y, Lidicky M, Zipperle J, Li L, Strasser K, et al. Mitochondria are gate-keepers of T cell function by producing the ATP that drives purinergic signaling. J Biol Chem. 2014;289(37):25936–45.
Turrens JF. Mitochondrial formation of reactive oxygen species. J Physiol. 2003;552(2):335–44.
Remans PH, Wijbrandts CA, Sanders ME, Toes RE, Breedveld FC, Tak PP, et al. CTLA-4IG suppresses reactive oxygen species by preventing synovial adherent cell-induced inactivation of Rap1, a Ras family GTPASE mediator of oxidative stress in rheumatoid arthritis T cells. Arthritis Rheum. 2006;54(10):3135–43.
Remans PH, Gringhuis SI, van Laar JM, Sanders ME, der Papendrecht-van der Voort EA, Zwartkruis FJ, et al. Rap1 signaling is required for suppression of Ras-generated reactive oxygen species and protection against oxidative stress in T lymphocytes. J Immunol. 2004;173(2):920–31.
Aichler M, Elsner M, Ludyga N, Feuchtinger A, Zangen V, Maier SK, et al. Clinical response to chemotherapy in oesophageal adenocarcinoma patients is linked to defects in mitochondria. J Pathol. 2013;230(4):410–9.
Zhu H, Chen L, Zhou W, Huang Z, Hu J, Dai S, et al. Over-expression of the ATP5J gene correlates with cell migration and 5-fluorouracil sensitivity in colorectal cancer. PLoS One. 2013;8(10):e76846.
Ni Chonghaile T, Sarosiek KA, Vo TT, Ryan JA, Tammareddi A, Moore Vdel G, et al. Pretreatment mitochondrial priming correlates with clinical response to cytotoxic chemotherapy. Science. 2011;334(6059):1129–33.
The authors would like to thank Manuela Le Bars and Corine Gaillez from Bristol-Myers Squibb for designing the study and providing study logistics. They also wish to thank the principal investigators of the APPRAISE study. APPRAISE principal investigators are Silvano Adami, Vivi Bakkenheim, Hilde Berner Hammer, Stefano Bombardieri, Maria-Antonietta D'Agostino, Paul Emery, Liana Euller-Ziegler, Gianfranco Ferraccioli, Maurizio Galeazzi, Philippe Gaudin, Walter Grassi, Annamaria Iagnocco, Herbert Kellner, Thierry Lequerré, Ingrid Möller, Esperanza Naredo, Mikkel Østergaard, Fredeswinda Romero, Istvan Szombati, Lene Terslev, Jacqueline Uson, Esther Vicente, Olivier Vittecoq and Richard Wakefield. The authors are grateful to Nikki Sabourin-Gibbs, Rouen University Hospital, for her help in editing the manuscript.
The APPRAISE trial was funded by Bristol-Myers Squibb and this ancillary study was funded by Bristol-Myers Squibb, which was involved in the design of the study and in the collection of data and samples.
Availability of data and materials
The complete datasets are available in the Gene Expression Omnibus database [GEO:GSE68215] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE68215). All materials were generated in Inserm unit 905 and they are available upon request in accordance with institutional regulations and policies.
Professor Thierry Lequerré had full access to all of the data in the study and takes responsibility for the integrity of the data and accuracy of data analysis. Design of the APPRAISE study: MAD, OV. Microarray Study conception: CD, OV, TL. Acquisition of data: CD, MH, GDT. Analysis and interpretation of data: CD, C.B., NV, OV, TL. Drafting of the manuscript: CD, PM, OV, TL. All authors were involved in drafting the article or revising it critically for important intellectual content and all authors have given final approval for the submitted version.
TL and OV received a grant from Bristol-Myers Squibb. CD, OV, TL, Bristol-Myers Squibb and Inserm Transfert have a patent application in conjunction with the University of Rouen.
Consent for publication
Ethics approval and consent to participate
The study was approved by the Institutional Review Board/Independent Ethics Committee (IRB/IEC) and local ethics committees, and was conducted in accordance with the ethical principles underlying the European Union Directive 2001/20/EC and the United States Code of Federal Regulations on Good Clinical Practice, as defined by the International Conference on Harmonisation. All patients provided written informed consent.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Variation in clinical and biological parameters between baseline and 6 months according to response (PDF 273 kb)
Transcripts differentially expressed in ABA/MTX responders (R) vs non-responders (NR) (PDF 305 kb)
Ineffectiveness of 87 mRNA associated with response to ABA/MTX to predict response to adalimumab or etanercept (PDF 49 kb)
Transcripts associated with Gene Ontology enrichment analysis (PDF 225 kb)
Weak correlation between inflammatory state or disease activity and gene expression level from the ETC pathway at baseline (PDF 201 kb)
About this article
Cite this article
Derambure, C., Dzangue-Tchoupou, G., Berard, C. et al. Pre-silencing of genes involved in the electron transport chain (ETC) pathway is associated with responsiveness to abatacept in rheumatoid arthritis. Arthritis Res Ther 19, 109 (2017). https://doi.org/10.1186/s13075-017-1319-8