Introduction

Rheumatoid arthritis (RA), the most common autoimmune disease, is thought to be a complex disease in which a combination of risk alleles from different susceptibility genes predisposes to development of the disease, following exposure to as yet unknown environmental factors. Several genome scans have suggested multiple RA loci [1], and recent case-control association studies have suggested new RA genes [2]. However, only HLADRB1 alleles have been both linked to and associated with RA, fulfilling the criteria for a fully demonstrated genetic factor [3].

The HLA-DRB1 locus accounts for approximately one-third of genetically determined susceptibility to the disease [4]. The identification of further RA susceptibility loci, both in candidate gene approaches and genome-wide linkage studies, has been hindered in the past by difficulties replicating such results in different study populations.

The role of the carnitine system in cell metabolism is mainly known in the mitochondria, where the interaction between fatty acid and glucose metabolism is fundamental for cell energy production [5]. Adequate carnitine levels are required for normal lipid metabolism and are important for energy metabolism [6]. An important component of the carnitine system is the plasma membrane carnitine transporters, also called organic cation transporters, i.e. OCTN1 and OCTN2 encoded by the SLC22A4 and SLC22A5 genes, respectively. Both genes map to the cytokine gene cluster on chromosome 5q31 and show 88% homology and 77% identity in their sequences. Although OCTN1 and OCTN2 are considered to be carnitine transporters, only OCTN2 is a high-affinity human carnitine transporter, while the carnitine transport activity of OCTN1 is very low [7].

In fact, a recent study has reported that the main substrate of this transporter is ergothioneine, an intracellular antioxidant with metal ion affinity that is transported one hundred times more efficiently than carnitine [8]. OCTN2 is widely expressed in many adult tissues, among them the pancreas, and it participates, at least in part, in proton/organic cation antiport at the renal apical plasma membrane level [9].

Previous reports have revealed the association of some polymorphisms within SLC22A4 and SLC22A5 genes with autoimmune complex diseases, namely RA and Crohn’s disease [10, 11]. The aim of this study was to investigate the association between SLC22A5 polymorphism – 207C>G (rs 2631367) and RA.

Material and methods

Patients

The study was carried out on 404 patients (322 women, 82 men, mean age 57.9 ±11.7 years) diagnosed with RA according to the criteria of the American College of Rheumatology (ACR) [12]. Patients were recruited from the outpatient and inpatient population of the Department of Rheumatology, County Hospital in Szczecin, Poland. All subjects were Caucasians from the Pomerania region of Poland. Subjects enrolled in the study underwent routine biochemical blood analysis and, when required, assays for anticardiolipin antibodies, antinuclear antibodies and immunological complexes. X-rays of the chest, hands and feet (erosive or non-erosive RA) were obtained in all patients, and, when necessary, radiographs of other joints. These were interpreted by two different expert radiologists.

The evaluation of subjects included physical examination with particular focus on the pattern of joint involvement, the presence of extra-articular features (such as vasculitis, anemia, sicca syndrome, amyloidosis, organ involvement), and laboratory features such as rheumatoid factor (RF). Amyloidosis was diagnosed by histomorphology (skin and bowel or duodenum biopsy), vasculitis by histomorphology (skin biopsy) and angiogram. The patients were treated with low doses of methotrexate (MTX) and glucocorticosteroids.

The control group consisted of 560 healthy subjects recruited from the Pomerania region (452 women and 108 men, mean age 59.9 ±12.6 years).

The study was approved by the local ethics committee and written informed consent was obtained from all subjects.

Evaluation of treatment efficacy

Good responders were defined as patients who were receiving MTX and had a DAS28 of ≤ 2.5 at 6 months of therapy (patients with remission of disease symptoms). Poor responders were defined as patients who were receiving MTX and had a DAS28 of > 2.5.

Genotyping

DNA was extracted from 200 μl of whole blood samples using a GeneMATRIX Quick Blood DNA Purification Kit (EURx, Poland). A SNP within the SLC22A5 gene – 207C>G (rs 2631367) was genotyped using pre-validated TaqMan genotyping assays (Life Technologies, USA). Fluorescence data were captured using a 7500 Fast Real-Time PCR System (Applied Biosystems, USA).

Statistical analysis

χ2 or Fisher’s exact test was used to compare genotype and allele frequencies between the study groups; p < 0.05 was considered statistically significant. The age at onset was compared between genotypes using the Kruskal-Wallis test.

Results

The distribution of SLC22A5 genotypes and alleles in RA patients did not differ significantly from that in healthy controls (Table I). As shown in Table II, there was no significant correlation between age at time of disease diagnosis and the SLC22A5 genotypes.

Table I

SLC22A5 genotype and allele distribution in rheumatoid arthritis patients and healthy controls

ParametersRA patientsControl group
n%n%
GG13232.715928.4
GC20350.228851.4
CC6917.111320.2
G allele46757.860654.1
C allele34142.251445.9

[i] RA – rheumatoid arthritis; no significant differences were noted (χ2 test).

Table II

. Association between SLC22A5 genotypes and clinical rheumatoid arthritis features

GenotypeAge of disease diagnosisRFErosive diseaseExtra-articular manifestations
n%n%n%
GG (n =132)46.9 ±12.89068.211587.12216.7
GC (n = 203)46.9 ±13.514470.916380.3199.4
CC (n = 69)47.7 ±13.04768.15884.157.2
p-valueaNSNSNS0.022b

RF – rheumatoid factor; NS – not significant

a χ2 test for nominal variables and Kruskal-Wallis test for age of disease diagnosis

b Fisher’s exact test for comparison of GG homozygotes vs. C allele carriers (CG + CC)

OR = 2.06 (95% CI: 1.11–3.85, p = 0.022).

Rheumatoid factor was diagnosed in 68.2% of subjects with SLC22A5 GG genotype, 70.9% with GC and 68.1% with CC. These differences were not significant (Table II).

Erosive RA was diagnosed in 87.1% of GG homozygous patients, in 80.3% of patients with the GC genotype and 84.1% of those with CC. These differences were not significant (Table II).

Extra-articular manifestations were diagnosed in 16.7% of SLC22A5 GG homozygous patients, in 9.4% with the GC genotype and in 7.2% of homozygous CC patients. The frequency of extra-articular manifestations was two-fold greater in homozygous GG patients as compared with carriers of the C allele (GG vs. GC + CC), OR = 2.06 (95% CI: 1.11–3.85, p = 0.022).

The efficacy of RA therapy with MTX is presented in Table III. Under MTX therapy remission of RA symptoms was achieved in 48.5% of patients with GG genotype, in 49.3% of patients with GC genotype and in 46.3% of patients with CC genotype. The differences were statistically non-significant.

Table III

Association between SLC22A5 gene polymorphism and response to treatment of rheumatoid arthritis patients with methotrexate

SLC22A5 genotypeGood responders n = 196Poor responders n = 208p-valuea
n%n%
GG (n = 132)6448.56851.5GG + GC vs. CCNS
GC (n = 203)10049.310350.7GG vs. GC + CCNS
CC (n = 69)3246.33753.7GG vs. CCNS

a Fisher’s exact test.

Discussion

In the present study we examined the association between SLC22A5 polymorphism and RA. The SLC22A5 genotypes did not correlate significantly with RA susceptibility, age at disease diagnosis, erosive or seropositive disease, or response to treatment with MTX, but were associated with extra-articular manifestations.

Organic cation transporters are among a large family of solute carrier transporters that number more than 200 in humans. Three organic cation transporters, SLC22A1 (OCT1), SLC22A2 (OCT2) and SLC22A3 (OCR3), were initially isolated and thereafter two more, SLC22A4 (OCTN1) and SLC22A5 (OCTN2), were isolated as a new family. The first two, SLC22A1 and SLC22A2, seem to transport organic cations in the renal basolateral membrane in a potentially dependent fashion. SLC22A5 is present in various tissues including kidney, skeletal muscle, heart, and placenta [13]. In the kidney SLC22A5 is expressed at the apical membrane of the proximal tubular epithelial cells. Expression was enhanced by tumor necrosis factor (TNF-α). A murine homologue of SLC22Q4 was found to be expressed in inflamed synovial tissue of collagen-induced arthritis [10].

A mouse strain with a point mutation of the mouse counterpart of human SLC22A5 exhibited an abnormal distribution of carnitine. This point mutation was identified as being involved in a familial carnitine metabolic disorder [14]. SLC22A5 is a physiological transporter of carnitine, and an in vitro analysis revealed that SLC22A5 transports tetraethylammonium (TEA) and carnitine. SLC22A4 and SLC22A5 are 76% homologous at the amino acid level, and an in vitro investigation of SLC22A4 indicated that it transports carnitine as well as TEA. However, the efficiency of carnitine transportation by SLC22A4 was far lower than that of SLC22A5.

The whole genome linkage study identified IBD5 as an inflammatory bowel disease (IBD)-linked locus [15]. Then the region was further closely evaluated with a dense SNP map that identified haplotypes consisting of two SNPs in SLC22A4 and SLC22A5. One of the SNPs substitutes 503 L to F with non-conservative effects on the tertiary structure of SLC22A4 and the other disrupts the heat shock element in the 5′-UTR of SLC22A5. Pharmacological assays revealed that the polymorphic amino acid substitution of SLC22A4 affected the transporting function of the molecule with several changes in Vmax and Km for some potential transport compounds. The allelic difference in the 5′UTR SNP in SLC22A5 was observed in the binding region of nuclear factors and affected in vitro transcription assays [11]. No conclusion was reached on whether one or both of the genes with the SNPs were responsible for the disease susceptibility.

So far the association between SLC22A5 and RA has not been widely investigated. Previous studies regarding the involvement of the SLC22A4 and SLC22A5 polymorphisms in RA pathogenesis are inconsistent.

Tokuhiro et al. [10] revealed a significant association between RA and the organic cation transporter gene SLC22A4. The authors showed that the expression of SLC22A4 was specific to hematological and immunological tissues, and that SLC22A4 was also highly expressed in the inflammatory joints of mice with collagen-induced arthritis. A SNP affects the transcriptional efficiency of SLC22A4 in vitro, owing to an allelic difference in affinity to runt-related transcription factor 1 (RUNX1), a transcriptional regulator in the hematopoietic system. A SNP in RUNX1 was also strongly associated with RA. These data indicate that the regulation of SLC22A4 expression by RUNX1 is associated with susceptibility to RA, which may represent an example of an epistatic effect of two genes on this disorder.

However, Barton et al. [16] found no evidence for an association between RA and either the SNP or the haplotype previously reported to be associated with RA in a Japanese population by Tokuhiro et al. [10]. These authors suggested that functional polymorphisms of the OCTN gene locus that have previously been associated with RA in Japanese subjects were not found to be associated with RA in a UK population. The findings do not provide support for a major role of these genes in the etiology of RA in this Caucasian population.

Similarly, Martinez et al. [17] studied whether SLC22A4/SLC22A5 haplotypes are relevant for RA predisposition in a Spanish population. The SLC22A4 and RUNX1 polymorphisms described previously did not show a significant role in RA susceptibility in the above population.

The results of a study by Komlósi et al. [18] suggested no influence of SLC22A4 C6607T genotypic variants on the circulating carnitine ester profile in patients with RA. Frigeni et al. [19] examined the functional properties of the OCTN1 and OCTN2 transporters encoded by SLC22A4 and SLC22A5 genes respectively. The results indicated that the OCTN1 transporter is tolerant of amino acid substitutions in the SLC22A4 gene, whereas the amino acid substitution in the SLC22A5 gene causes reduced transporter activity. Previous studies have shown that the activity of SLC22A5 transport protein in exosomes involved in RA pathogenesis is regulated by the pro-inflammatory cytokine INF-γ [20, 21].

Conclusions

The results of the present study indicated that the SLC22A5 polymorphism does not influence RA susceptibility or the response to methotrexate treatment in studied RA patients from a Polish population, but may be associated with the development of extra-articular manifestations of RA.