Adenosine receptor expression in rheumatoid synovium: a basis for methotrexate action

Introduction Methotrexate (MTX) exerts at least part of its anti-inflammatory effects through adenosine receptors (ADOR). The aims of this study were to determine the expression of all four adenosine receptor genes (ADORA1, ADORA2A, ADORA2B, ADORA3 and ADORA3variant) in rheumatoid synovial tissue and any influence of MTX exposure on this expression. Furthermore, we investigated whether polymorphisms within ADORA3 were associated with response and/or adverse effects associated with MTX. Methods Adenosine receptor gene expression was undertaken using PCR in 20 rheumatoid arthritis (RA) synovial samples. A separate cohort of 225 RA patients receiving MTX was genotyped for SNPs in the ADORA3 receptor gene. Double immunofluorescence was used to identify cells expressing ADOR protein. Results All ADOR genes were expressed in all synovial samples. ADORA3 and A3variant were the dominant subtypes expressed irrespective of MTX therapy. Expression of ADORA2A and ADORA2B was increased in patients receiving MTX compared to those not receiving MTX. There was no association between the ADORA3 rs1544224 SNP and high and low disease activity or MTX-associated adverse effects. ADORA2B protein expression was most obvious in vascular endothelial cells whereas ADORA3 protein was more abundant and expressed by synovial fibroblasts. Conclusions We have shown that adenosine receptors are expressed in RA synovium. There is differential expression of receptors such that ADORA3 is expressed at significantly higher levels. This evidence demonstrates the potential for MTX to exert its anti-inflammatory effects at the primary site of pathology within the joints of patients with RA.


Introduction
Methotrexate (MTX) remains the first line drug for the treatment of rheumatoid arthritis (RA). The exact mechanism of action of MTX is complex. Following administration and absorption, serum MTX concentrations fall rapidly [1] as MTX is transported into a variety of cells including red blood cells (RBC) and synoviocytes. Intra-cellularly, glutamate moieties are added by folylpolyglutamate synthetase (FPGS) such that MTX is retained within cells as MTX polyglutamates (MTXGlu n ).
MTX polyglutamates inhibit a number of important intracellular enzymes in the folate pathway including dihydrofolate reductase and thymidylate synthase. Inhibition of 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) transformylase (ATIC) results in accumulation of AICAR and increased adenosine release into the circulation. Extracellular adenosine increases cAMP which inhibits production of proinflammatory cytokines including TNF-α, IFN-γ and IL-1β, which are important in the inflammatory process in RA.
Adenosine acts via four G-coupled adenosine receptors (ADOR); ADORA 1 , ADORA 2A , ADORA 2B and ADORA 3 . Activation of the ADORA 1 and ADORA 2B results in pro-inflammatory effects which include an increase in pro-inflammatory cytokine release from neutrophils and monocytes. In contrast, activation of the ADORA 2A and ADORA 3 results in anti-inflammatory effects with a reduction in IL-1β, TNF-α, IL-6, and decreased neutrophil superoxide production. In part, these opposing effects reflect patterns of ADOR subtype expression associated with distinct cell lineages and differentiation states, as well as linkage to divergent G protein mediated signaling pathways [2]. ADORA 3 is highly expressed in rat synovium, peripheral blood mononuclear cells (PBMC) and lymphocytes [3]. In murine studies, MTX has been shown to exert its anti-inflammatory effects via ADORA 2A and ADORA 3 [4]. In MTX-treated RA patients the ADORA 3 has been reported to be expressed at higher levels on PBMC compared to healthy controls [5]. Furthermore, expression of ADORA 2A , ADORA 2B and ADORA 3 appear to be regulated by inflammatory cytokines such as TNF-α [6][7][8].
In an experimental model using cultured human monocytes MTX, an ADORA 1 agonist and an ADORA 2 antagonist were shown to enhance the formation of multi-nucleated giant cells from monocytes. The authors suggest that within the joint MTX may increase adenosine concentrations to the extent that ADORA 2 are bound with resulting anti-inflammatory effects [9].
Whether the adenosine receptors are expressed in human joint synovial tissue, the primary lesion in RA is unknown. Furthermore, MTX through its anti-inflammatory effects may alter the expression of adenosine receptors within the synovium. We hypothesized that synovial tissue will express adenosine receptors and there will be a relationship between MTX exposure and adenosine receptor expression within the synovium. Furthermore, we hypothesized that polymorphisms within ADORA 3 may be associated with response and/ or adverse effects associated with MTX.

Methods
Ethical approval was obtained from the Multi-Regional Ethics Committee and the Upper South B Regional Ethics Committee, New Zealand. All patients gave written informed consent obtained according to the Declaration of Helsinki. Two cohorts of patients with RA as defined by the 1987 American Rheumatism Association criteria [10] were recruited. The first cohort (n = 20) were used for synovial tissue gene expression studies while the second cohort (n = 225) were used for the ADORA 3 genotyping. There was no overlap of the two patient cohorts.

Synovial tissue samples -PCR and QT-PCR
Synovial tissue was obtained from the first cohort of twenty RA patients undergoing joint surgery. Total RNA was extracted from the synovial tissue using Qiagen RNeasy mini kits. Either 0.5 or 1 μg of RNA was reverse transcribed at 42°C for 50 minutes using Superscript III (Life Technologies, Carlsbad, CA USA) and Oligo (dT) 12 -18 primers. Expression of ADORA 1 , A 2A , A 2B , A 3 and A 3variant were examined. Quantitative real-time PCR was undertaken using TaqMan gene expression assays (Applied Biosystems, Foster City, California, USA) for ADORA 1 (Assay ID: Hs00379752_m1), ADORA 2A (Hs00169123_m1) and ADORA 2B (Hs00386497_m1). We also measured three alternative transcripts produced from the ADORA 3 gene -the ADORA 3 assay (Hs00252933_m1) detects transcripts 1 and 3 while the ADORA 3variant assay (Hs00181232_m1) detects transcript variant 2 only. All values were normalized to GAPDH (Hs999999905_m1).

DNA sequencing and genotyping of variants in ADORA 3
A second cohort consisting of 234 RA patients receiving MTX, were recruited as part of an unrelated cross-sectional study [11]. DNA was extracted from 5 ml samples of peripheral blood using phenol-chloroform [12][13][14]. The promoter, 5' and 3' UTRs, and the open reading frame (ORF) of the ADORA 3 gene were amplified in five overlapping 680 bp to 750 bp fragments and subsequently sequenced in 15 patients with low disease activity (Disease Activity Score in 28 Joints (DAS28) < 3.2) and 15 patients with high disease activity (DAS28 ≥ 3.2). Briefly, each PCR was performed in a total volume of 20 μl containing 200 μM dNTPs, 2 mM MgCl 2 , 0.5 μM of the relevant forward and reverse primers [see Additional file 1, Table S1], 1 U of Platinum ® Taq DNA polymerase (Invitrogen, Carlsbad, CA, USA) and approximately 50 ng of DNA. Thermal cycling conditions comprised an initial denaturation of 94°C for two minutes, followed by 35 cycles of 94°C for 30 seconds, 60°C for 30 seconds, 72°C for one minute, and a final extension of 72°C for three minutes. Successful amplification was confirmed by analyzing 5 μl of each PCR by 1% agarose electrophoresis. Complementary DNA was purified from the remaining 15 μl of each PCR using PureLink™ PCR Purification Kit (Invitrogen) according to the manufacturer's instructions. Approximately 6 to 7 ng of purified PCR product was sequenced in both directions using BigDye ® Terminator Version 3.1 chemistry on an ABI 3730xl DNA Analyzer.
Genotyping of the 5'UTR SNP rs1544224 within ADORA 3 was conducted using a two-tube allele-specific PCR assay. Each reaction was performed in a total volume of 10 μl containing 200 μM dNTPs, 2 mM MgCl 2 , 0.5 μM of the control primers β2Mf and β2Mr, 0.5 μM of the common forward primer ADORA 3U5f [see Additional file 1, Table S1], 1 U of Platinum Taq DNA polymerase (5 U/μl) (Invitrogen) and approximately 50 ng of DNA. Reaction 1 contained 0.5 μM of primer rs1544224WT and reaction 2 contained 0.5 μM of primer rs1544224MU [see Additional file 1, Table S1]. Thermal cycling conditions comprised an initial denaturation of two minutes at 94°C, followed by 30 cycles of 30 seconds at 94°C, 30 seconds at 65°C, 30 seconds at 72°C, and a final extension of two minutes at 72°C. Sequence-validated positive controls for each genotype (rs1544224 C/C, C/T, T/T) were included in all genotyping runs. Reactions were visualized using 3% agarose gel electrophoresis. The primers β2Mf and β2Mr amplified a 567 bp region of the beta-2-microglobulin gene (β2M) which served as an internal control for amplification in samples where homozygosity for the SNP of interest resulted in allele-specific PCR product in only one reaction. The accuracy of the genotype was checked by repeat analysis of 5% of the samples and by sequencing PCR amplicons from a further 5% of the samples.

Immunohistology for ADOR protein expression
Cell-specific expression of ADOR protein was investigated using two-color immunofluorescence. Primary cell-specific monoclonal antibodies (mAbs) were used in combination with goat polyclonal antibody specific for ADORA 2B (1:50; Abcam, Cambridge, UK) as previously described [15]. For the detection of ADORA 2A or ADORA 3 , dual labeling with mAbs was utilized, based on the following outline. However, certain combinations of cell-and ADORA-specific mAbs required that the AlexaFluor-conjugated, detection antibodies were reversed. Briefly, tissue sections were first incubated overnight at 4°C with mAbs specific for human ADORA 2A (1:100; Lifespan Biosciences, Seattle, Washington, USA) or ADORA 3 (1:100; Abnova, Taipei, Taiwan), diluted in RPMI-1640 (GibcoBRL, Grand Island, New York, USA)/10% FCS. Sections were washed and bound mAb was detected by incubation with Alexa-Fluor568-conjugated goat anti-mouse IgG (Invitrogen) for 1.5 hours diluted in RPMI/10% FCS. Sections were then sequentially washed, incubated with 1:1 PBS/10% FCS: normal mouse serum (DAKO, Glostrup, Denmark) for 30 minutes, washed, then incubated with an excess of goat anti-mouse IgG(H+L) Fab fragments (Jackson Immuno Research, Baltimore, Philadelphia, USA) for one hour before a final wash. Primary cell-specific mAbs were applied for 1.5 hours, and the sections washed again before bound primary antibody was detected with AlexaFluor488-conjugated goat anti-mouse IgG (Invitrogen) for 1.5 hours diluted in RPMI/10% FCS containing 1.25 mg/ml Hoechst 33342 nuclear stain (Molecular Probes ® Invitrogen, Eugene, Oregon, USA).

Statistical analysis
The test for Hardy-Weinberg equilibrium (HWE) was performed using the internet-based program SHESIS [16]. Assessment of linkage disequilibrium between the two 5'UTR SNPs rs1544223 and rs1544224 was performed using Haploview 4.2 and HapMap Data Rel 28 Phase II + III (assembly dbSNP b126).
Data were analyzed using Graph-pad Prism 4. All data are presented as mean ± standard deviation (SD) unless otherwise stated. The statistical significance of the differences between the groups was determined by the Mann-Whitney U test and correlation was tested using the Spearman's test.

Results
Demographics for the two patient cohorts are shown in Table 1.

Expression of ADORA in synovial tissue
All five ADOR genes were expressed in all 20 synovial tissues. ADORA 3variant (A 3V ) and ADORA 3 were the dominant ADOR subtypes present. This was followed by ADORA 2B and ADORA 1 . ADORA 2A was the least expressed ( Figure 1). There was a significant correlation between ADORA 2A and ADORA 2B expression and between ADORA 3 and ADORA 3V expression ( Table 2).
Expression of ADORA 2A and ADORA 2B was significantly higher in those patients receiving MTX (n = 11) compared to those not receiving MTX (n = 9). However, MTX had no apparent influence on expression of the other ADOR genes ( Figure 2). There was no relationship between MTX dose and the amount of ADOR gene expression. Only three patients were receiving anti-TNF therapy without concomitant MTX. There was no significant difference in the expression of any ADOR genes between these three patients and the remaining patients (P > 0.05; data not shown). Similarly there was no significant difference in ADOR gene expression between those patients receiving prednisone (n = 6) and those not (n = 14) (P > 0.05; data not shown). There was no significant correlation between expression of any of the ADOR genes and ESR or serum CRP levels (data not shown).

Cells expressing ADOR genes in RA synovial tissue
Consistent with the gene expression there was widespread expression of ADOR A3 protein within rheumatoid synovium, co-localized to fibroblasts (Figure 3). Some, but not all infiltrating T cells stained weakly for ADOR A3 protein, as did interdigitating cells in the T cell areas of synovial lymphoid follicles. In contrast ADORA 2B protein was most obvious in vascular endothelial cells (ECs) that also expressed Factor VIII ( Figure 3). The ADORA 2B -positive ECs had a cuboidal morphology, consistent with activated endothelium. However ADORA 2B protein was also expressed by vascular ECs with a flattened morphology in osteoarthritis synovium (data not shown). ADOR 2B protein was also expressed by fibroblasts and by some aggregated T cells in rheumatoid synovia but in amounts considerably less than was observed in ECs. Expression of ADORA 2A protein was rare and was generally restricted to large isolated cells with a dendritic appearance. These latter cells have not been further identified. However ADOR 2A protein was also observed in a solitary germinal centerpositive lymphoid follicle, consistent with follicular dendritic cell expression.

Sequence variants in ADORA 3
A total of 15 SNPs were detected by sequencing ADORA 3 across 15 patients with low disease activity (DAS28 ≤ 3.2) and 15 patients with high disease activity (DAS28 > 3.2) [See Additional file 2, Figure S1]. Of these polymorphisms, two were located in the promoter, four in the 5'UTR, four in the ORF, and five in the 3'UTR. All SNPs identified have been previously reported in one or more populations. Of the SNPs within the ORF, two were synonymous (rs2789537 and rs2229155) and two were non-synonymous (rs35511654 and rs2800889). Although none of the variants were found exclusively in groups with high or low disease activity, the minor allele frequencies (MAF) of the 5'UTR SNPs rs1544223 and rs1544224 were greater in the high disease activity group (MAF high disease activity = 0.37 versus MAF low disease activity = 0.13, P = 0.037). As a consequence, SNP rs1544224 was genotyped in the complete RA cohort using allele-specific PCR. The second SNP rs1544223 was not genotyped as it was found to be in complete linkage disequilibrium (LD) with rs1544224 in the subset of patients sequenced and almost complete LD (r 2 = 0.974) in analysis of Utah residents with Northern and Western European ancestry (CEU) and Tuscan in Italy (TSI) from the HapMap Project.

Association between ADORA 3 rs1544224 SNP and disease activity, and MTX toxicity
Genotyping for rs1544224 was successful in 96.2% of patients (225/234) and no deviation from HWE was observed (P = 0.43). There was no association between ADORA 3 rs1544224 and any of the disease activity   (Table 3). There was no significant association between ADORA 3 rs1544224 genotype and low or high disease activity groups (P = 0.76). Patients homozygous for the minor allele were more likely to have 'hair loss' (P = 0.02) and 'forgetfulness' (P = 0.02), but there was no other difference in adverse effects between the genotypes [See Additional file 3, Table S2].

Discussion
MTX is an effective drug in the management of RA. Its ability to suppress inflammation results, at least in part, from its ability to increase extracellular concentrations of adenosine which stimulates adenosine A 2A and A 3 receptors [4]. The primary site of inflammation in RA is the synovium. Suppression of the inflammatory process that drives this inflammation in the synovium and formation of the tissue destructive pannus is critical to the control of inflammatory joint symptoms and prevention of joint erosion in RA.
Most attention with respect to the anti-inflammatory mechanism of action of MTX has focused on the folate pathway with less attention on the adenosine pathway. ADORA 3 has been reported to be expressed in PBMCs [6,17] and ADORA 1 and A 2 in synovial fluid cells from patients with RA [18] and isolated synoviocytes [19]. Both the adenosine A 2A and A 3 receptors, but not ADORA 1 or ADORA 2B , have been reported to be upregulated in lymphocyte and neutrophil membranes from patients with RA as compared to healthy controls [20]. However, there have been no previous studies examining expression of the adenosine receptors within synovial tissue from patients with RA. Herein we have demonstrated that the genes for all four adenosine receptors are expressed at varying levels in rheumatoid synovium.
ADORA 3 was expressed significantly more than the other ADOR genes. This may be the result of the inflammatory process occurring within the synovium, with previous reports suggesting that expression of ADORA 3 is up-regulated by TNF-α via activation of NFB [17,20]. ADORA 2A , through which MTX is also thought to exert its anti-inflammatory effects, was expressed in the lowest amounts of all the ADOR subtypes [7]. Several possibilities could explain the difference in expression of the various adenosine receptor subtypes. These include the cellular constituents of the tissue examined. Consistent with this possibility, ADORA 3 protein expression was more widespread in synovial tissue. In contrast ADORA 2B protein was prominent in vascular ECs but also present in fibroblasts and some T cells, albeit in much reduced amounts; ADORA 2A protein was restricted to cells with a dendritic appearance. Overall, the patterns of ADOR protein distribution were consistent with the hierarchical order we observed for expression of the corresponding ADOR genes. For all ADOR subtypes, we did not observe significant protein associated with the infiltrating inflammatory cells. We are not able to exclude that these patterns of protein expression and the differences in gene expression could also be affected by the inflammatory cytokine milieu within the synovium. However, we further considered the effect of anti-rheumatic therapy and the possible influence on ADORA expression.
Both MTX and anti-TNF therapy have been reported to have effects on adenosine receptor expression and affinity. MTX therapy results in up-regulation of the ADORA 2A and ADORA 3 , but not ADORA 1 or ADORA 2B expression in lymphocyte and neutrophil membranes from RA patients as compared to healthy controls [20]. The up-regulation of ADORA 3 expression on PBMCs has been observed after only ten  weeks of MTX therapy [5]. In addition, the affinity of the ADORA 2A and ADORA 3 was lower in patients receiving MTX compared to healthy controls [20]. In comparison, in RA patients receiving anti-TNF therapy the expression and affinity of ADORA 2A and ADORA 3 was similar to that observed in healthy controls [20]. This is consistent with an effect of TNF or the indirect effect of inflammation in general. In our cohort, expression of ADORA 2A and ADORA 2B but not ADORA 3 , was increased in those patients receiving MTX compared to those not receiving MTX. While we observed no effect of anti-TNF therapy, the number of patients was small. The presence of the adenosine receptors in the synovial membrane provides evidence that MTX has the potential to exert important antiinflammatory effects at the primary site of the inflammatory process in RA.
It is important to note that the specimens used in our study were from RA patients with late stage disease and these patients differed from the cohort used for the genotyping studies. There remains a possibility that changes in ADOR subtype expression related to MTX therapy could be different at earlier stages of disease. However, given that expression of ADORA 2A , ADORA 2B and ADORA 3 is regulated by inflammatory cytokines such as TNF-α [6][7][8] it would be reasonable to assume there may be similar findings. This will need to be confirmed in studies of patients with early RA.
Signalling through ADOR depends on a variety of factors including receptor expression, receptor sensitivity and extra-cellular concentration of ligand. Our data show greatest gene expression of ADOR 3 and the ADOR 3VAR in rheumatoid synovia and that expression of a designated ADOR subtype is not necessarily exclusive to a particular cell type. Furthermore, there are known variations in the kinetic properties of the different ADOR subtypes, with ADORA 1 , ADORA 2A and ADORA 3 several orders of magnitude more sensitive to adenosine than ADORA 2B [21]. It remains to be determined whether variation in expression compensates for the different kinetics attributed to the different ADOR subtypes and what consequences this has for individual synovial cell types.
The importance of ADORA 3 in the anti-inflammatory process has also been highlighted by the finding in activated PBMCs that addition of an adenosine agonist (IB-MECA, CF101) resulted in reduced ADORA 3 expression and down-regulation of TNF-α, thus breaking the autocrine loop and inhibiting the inflammatory process [17]. The clinical efficacy of an oral adenosine A 3 receptor agonist, CF101, in patients with RA has been examined in a small 12-week study. The level of expression of ADORA 3 in PBMCs at baseline correlated with ACR50 and ACR70 responses (P = 0.036) [22]. We have now demonstrated high levels of ADORA 3 expression in the joint tissues targeted by the inflammation. These data provide an encouraging line of evidence for further investigation of adenosine A 3 agonists as therapeutic agents in RA, particularly in patients with high levels of ADORA 3 expression.
MTX has also been suggested to have a beneficial effect on cardiovascular mortality in RA, an effect that is not observed with other disease modifying anti-rheumatic drugs (DMARDs) [23,24]. This effect appears to be mediated through ADORA 2A activation mediating reverse cholesterol transport and limiting foam cell formation [25]. One of the key challenges in the management of RA is rapid and effective disease control without provoking adverse effects. Currently there is no reliable means to predict which patients will respond to MTX, nor who will experience adverse effects. In this regard efforts have focused on the ability of polymorphisms of genes involved in the folate pathway to predict efficacy and toxicity of MTX. To date no SNP or combination of SNPs within the folate pathway has been identified that will reliably predict response or adverse effects associated with MTX.
Polymorphisms of adenosine monophosphate deaminase 1 (AMPD1), aminoimidazole carboxamide ribonucleotide transformylase (ATIC), and inosine triphosphate pyrophosphatase (ITPA), which are involved in the adenosine pathway, have been associated with good clinical response to MTX in patients with RA in some but not all studies [26,27]. Genetic variations in the ADOR genes have received less attention but have been reported to be associated with good clinical response to MTX [26] or have no effect [27]. Polymorphisms within the ADORA 2A gene have been associated with adverse gastrointestinal events (nausea, vomiting or diarrhea) associated with MTX in patients with RA [28]. Herein we have investigated polymorphisms within the ADORA 3 gene. We have shown that the SNP rs1544224 is not associated with disease activity in patients with RA receiving MTX. Whether this SNP can predict MTX efficacy will need to be formally examined in prospective clinical studies. The cross-sectional design of this study resulted in the cohort of patients being enriched for patients who tolerated MTX. Patients with more severe adverse effects such as hepatotoxicity or myelotoxicity would have discontinued MTX and not entered this study. Nevertheless, there are a number of other adverse effects associated with MTX for which patients elect not to discontinue therapy because the benefits of therapy outweigh the adverse effects. We have shown a weak association between the ADORA 3 rs1544224 SNP and forgetfulness and hair loss in our cohort. Confirmation of this will be required in other cohorts.