Skip to main content
Download PDF
Download ePub
  • Research article
  • Open access
  • Published:

Genetic polymorphisms of enzyme proteins and transporters related to methotrexate response and pharmacokinetics in a Japanese population

Journal of Pharmaceutical Health Care and Sciences volume 2, Article number: 35 (2016) Cite this article

Abstract

Background

Methotrexate (MTX) is currently the anchor drug widely used worldwide in the treatment of rheumatoid arthritis (RA). However, the therapeutic response to MTX has been shown to vary widely among individuals, genders and ethnic groups. The reason for this has been not clarified but it is considered to be partially due to several mechanisms in the cellular pathway of MTX including single-nucleotide polymorphisms (SNPs). The purpose of this study was to investigate the allelic frequencies in different ethnic and/or population groups in the 10 polymorphisms of enzyme proteins and transporters related to the MTX response and pharmacokinetics including MTHFR, TYMS, RFC1, FPGS, GGH, ABCB1, ABCC2 and ABCG2 in unrelated healthy Japanese adults and patients with RA.

Methods

Ten polymorphisms, methylenetetrahydrofolate reductase (MTHFR) 1298, thymidylate synthase (TYMS) 3'-UTR, reduced folate carrier 1 (RFC1) 80 and−43, folypolyglutamyl synthase (FPGS) 1994, γ-glutamyl hydrolase (GGH) 452 and−401, the ABC transporters (ABCB1 3435, ABCC2 IVS23 + 56, ABCG2 914) of enzyme proteins and transporters related to MTX response and pharmacokinetics in 299 unrelated healthy Japanese adults and 159 Japanese patients with RA were investigated to clarify their contributions to individual variations in response and safety to MTX and establish personalized MTX therapy. SNPs were evaluated using real-time polymerase chain reaction (PCR).

Results

Comparison of allelic frequencies in our study with other ethnic/population groups of healthy adults and RA patients showed significant differences in 10 polymorphisms among healthy adults and 7 among RA patients. Allelic frequencies of MTHFR 1298 C, FPGS 1994A and ABCB1 3435 T were lower in Japanese than in Caucasian populations and those of ABCC2 IVS23 + 56 C and ABCG2 914A were higher in Japanese than in Caucasian/European populations in both healthy adults and RA patients. Allelic frequencies of MTHFR 1298 C, GGH−401 T, ABCB1 3435 T, and ABCG2 914A were higher in healthy Japanese adults than in an African population, and those of RFC1 80A, RFC1−43C and ABCC2 IVS23 + 56 C in healthy Japanese adults were lower than in Africans. However, no significant differences were seen in the distribution of allelic frequencies between healthy Japanese adults and RA patients.

Conclusion

The variations in allelic frequencies in different ethnic and/or population groups in healthy adults and RA patients may contribute to individual variations in MTX response and toxicity.

Background

Rheumatoid arthritis (RA) is a systemic autoimmune inflammatory disease and its pathogenesis remains unclear. Although no curative therapy for RA has been established, the therapeutic goal is to delay symptom progression using disease-modifying anti-rheumatic drugs, gold preparations and biologics and achieve pain relief using nonsteroidal anti-inflammatory drugs (DMARDs). Among the therapeutic options, methotrexate (MTX) is currently the anchor drug widely used worldwide in the treatment of RA because it is inexpensive, effective and safe. In the 2014 Japanese therapeutic guidelines for RA treatment published by the Japan College of Rheumatology (JCR) [1], MTX was also cited as a first-line drug. However, the therapeutic response to MTX has been shown to vary widely among individuals, genders and ethnic groups [2, 3]. The reason for this has been not clarified but it is considered to be partially due to several mechanisms in the cellular pathway of MTX including single-nucleotide polymorphisms (SNPs) in transporters, glutamation, the folate pathway and adenosine pathway [2] such as reduced folate carrier 1 (RFC1), folypolyglutamyl synthase (FPGS), γ-glutamyl hydrolase (GGH), the ABC transporters ABCB1, ABCC2 and ABCG2, methylenetetrahydrofolate reductase (MTHFR), and thymidylate synthase (TYMS). We previously reported that the distribution of MTHFR C677T between black and Japanese populations, of TYMS 5'-UTR alleles between Caucasian or black and Japanese populations, and of TYMS 3'-UTR alleles between Caucasian and Japanese populations showed significant differences, as well as gender differences in TYMS 3'-UTR allelic frequency in Japanese [3].

Up to the present, a clear genetic difference has not been confirmed between healthy individuals and RA patients which would indicate susceptibility to RA, as such genetic variations do for many other diseases. There are few data on the possible genetic differences in RA, and therefore, it is important to clarify them to obtain useful information from the viewpoint of human genetics and/or population genetics, leading to the establishment of personalized medicine for those susceptible to the development of RA.

The purpose of this study was to investigate the allelic frequencies in different ethnic and/or population groups in the 10 polymorphisms of enzyme proteins and transporters related to MTX response and pharmacokinetics including MTHFR, TYMS, RFC1, FPGS, GGH, ABCB1, ABCC2 and ABCG2 in unrelated healthy Japanese adults and patients with RA.

Methods

Volunteers

A total of 299 unrelated healthy Japanese adult volunteers (148 men and 151 women) were recruited from Sumida Hospital (Tokyo) and 159 unrelated Japanese adult patients with RA (40 men and 119 women) admitted to PS Clinic (Fukuoka, Kyushu) were enrolled in this study. All were older than 20 years of age. MTX was administered to 99 (62.3%) of 159 RA patients, while 43 (27%) were receiving DMARDs, and 55 (34.6%) biologics as combination drugs. The daily dosage of MTX was 4 mg (3 patients), 6 mg (11 patients), 8 mg (75 patients), 10 mg (7 patients), and 12 mg (3 patients).

DNA extraction

For genetic analysis, a 5-mL peripheral blood sample was obtained from each study participant using the standard venipuncture technique. The samples were placed in tubes containing ethylenediaminetetraacetic disodium salt (EDTA 2Na) and stored at−20 °C until DNA extraction. DNA was extracted using the Wizard Genomic DNA Purification Kit (Promega Corporation, WI, USA) following the manufacturer’s instructions. Total genomic DNA was quantified and its purity and integrity were analyzed using the NanoDrop 1000 spectrophotometer v3.7 (Thermo Scientific, Wilmington, DE, USA).

Allele genotyping

Genotyping of alleles of MTHFR 1298A > C (rs1801131), TYMS 3'-UTR (−6/+6) (rs16430), RFC1 80G > A (rs1051266), RFC1−43 T > C (rs1131596), FPGS 1994G > A (rs10106), GGH 452C > T (rs11545078), GGH 401C > T (rs3758149), ABCB1 3435C > T (rs1045642), ABCC2 IVS23 + 56 T > C (rs4148396), and ABCG2 914C > A (rs2231142) was performed using the TaqMan SNP Genotyping Assay from Applied Biosystems (Foster City, CA, USA) with fluorogenic binding probes. PCR amplification with the real-time PCR method was performed in 25 μL of reaction mixture including genomic DNA 20 ng (60 ng for FPGS 1994G > A), 0.63 μL of 40 × TaqMan SNP Genotyping Assay Mix (0.32 μL for RFC1 80G > A), and TaqMan Universal PCR Master Mix 12.5 μL. The PCR reaction conditions were: initial denaturation for 10 min at 95 °C; 40 cycles at 92 °C/15 s; and annealing and extension for 1 min at 60 °C (55 °C for FPGS 1994G > A). For genotyping of alleles of ABCB1 3435C > T and ABCG2 914C > A, the PCR reaction conditions were: initial denaturation for 10 min at 95 °C; 50 cycles at 92 °C/15 s; and annealing and extension for 1 min at 60 °C. The PCR system used was the StepOnePlus real-time PCR system (Applied Biosystems).

Assessment of ethnic or population and gender differences

We identified eligible comparative studies by searching MEDLINE (1966–May 2016), Embase (1974–May 2016), the Cochrane Library (Issue 1 of 12, January 2016), Japana Centra Revuo Medicina (Ichushi) (1981–May 2016), and the database of the HapMap Project (accessed in May 2016) [4] on the assumption that it contained healthy adult data. The search terms were MTHFR, methylenetetrahydrofolate reductase, TS, TYMS, thymidylate synthase, dihydrofolate reductase, RFC1, reduced folate carrier 1, Solute carrier family 19 (folate transporter) member 1, SLC19A1, folylpolyglutamate synthase, FPGS, gamma-glutamyl hydrolase, GGH, ATP-binding Cassette Sub-family B Member 1, ABCB1, ATP-binding cassette sub-family C member 2, ABCC2, ATP-binding cassette transporter G2 and ABCG2, in combination with genotype and polymorphism. All languages were included. Additionally, a manual search of reference listings from all of the articles retrieved from the electronic databases was performed. Selection criteria were case–control, cross-sectional or prospective cohort studies, which assessed any allelic frequencies of MTHFR 1298 A > C, TYMS 3'-UTR 6-bp deletion, RFC1 80 G > A, RFC−43 C > T, FPGS 1994 G > A, GGH 452 C > T, GGH−401C > T, ABCB 3435 C > T, ABCC2 IVS23 + 56 T > C, and ABCG2 914 C > A in healthy adults and RA patients and in which ethnic group was clearly identified. In case–control studies, only the allelic frequencies of the healthy adult control group were extracted. If we could not identify the ethnic group in the literature, population groups were assumed to be from the country of residence of the lead author.

Statistical analysis

Comparisons of ethnic and gender differences in the distribution of allelic frequencies among genotypes, and/or gender differences and tests for Hardy-Weinberg equilibrium were carried out using the chi-square test. A p-value of <0.05 was considered to represent a significant difference in all statistical analyses. All statistical analyses were performed using JMP 12.0.1 software (SAS Institute Inc., Cary, NC, USA).

Results

Allelic frequencies in healthy adult populations and RA patients

Among the 298 studies retrieved from the electronic databases and those retrieved from their references by manual search, a total of 14,000 healthy adults from 16 studies [4–19] and 4284 RA patients from 25 studies [5, 6, 9, 14, 20–40], involving HapMap Projects [4] that described Japanese, Asian, European and African groups, were included in the comparative analysis (Tables 1 and 2). The number of individuals and published genotype frequency data for each ethnic/population group and comparisons with our Japanese allelic frequency results are shown in Tables 1 and 2.

Table 1 Comparison of distribution of MTX-related enzyme gene and transporter polymorphisms among ethnic or population group in healthy adults
Full size table
Table 2 Comparison of distribution of MTX-related enzyme gene and transporter polymorphisms among ethnic or population group in the patients with rheumatoid arthritis
Full size table

The distributions of the RFC1 80G > A allele in African [4, 11] in healthy adults and Caucasian [24, 34–36] in RA patients, RFC1−43 T > C allele in Caucasian [4] in healthy adults, ABCB1 3435 C > T allele in Caucasian [4] in healthy adults, and ABCC2 IVS23 + 56 T > C allele in Asian [4] and African [4] healthy adults were not in agreement with Hardy-Weinberg equilibrium. The distributions of the MTHFR 1298 A > C allele in Caucasian [23–31] in RA patients were also not in agreement with Hardy-Weinberg equilibrium.

Statistically significant differences were found in the distributions of the MTHFR 1298A > C, RFC1 80 G > A, RFC−43 T > C, FPGS 1994 G > A, GGH−401 C > T, ABCB1 3435 C > T, ABCC2 IVS23 + 56 T > C, and ABCG 914 C > A alleles in healthy adults between our Japanese population and Caucasians and Africans. Statistically significant differences were found in the distributions of the RFC1 80 G > A and ABCC2 IVS23 + 56 T > C alleles in healthy adults between our Japanese population and Asians (Table 1).

Statistically significant differences were found in the distributions of the MTHFR 1298A > C, TYMS 3'-UTR−6 > +6, RFC1 80 G > A, FPGS 1994 G > A, GGH 452 C > T, ABCB1 3435 C > T, ABCC2 IVS23 + 56 T > C, and ABCG 914 C > A alleles in RA patients between our Japanese population and Caucasians (Table 2).

Allelic frequencies and gender differences in frequencies in healthy Japanese adults

The distribution of MTHFR 1298A > C, and TYMS 3'-UTR (−6/+6), RFC1 80G > A, RFC1−43 T > C, FPGS 1994G > A, GGH 452C > T, GGH−401C > T, ABCB1 3435C > T, ABCC2 IVS 23 + 56 T > C, and ABCG2 914C > A polymorphisms in 299 healthy Japanese adults are summarized in Table 3. The proportions of genotypes at each site for their transporter carrier proteins and enzymes were generally in agreement with Hardy-Weinberg equilibrium. Therefore, the proportions at each site followed Mendelian principles. Only three genotype combination patterns were found for RFC1 80G > A and RFC1−43 T > C, consisting of 80G/G,−43 T/T genotype, 80G/A,−43 T/C genotype, and 80A/A,−43C/C genotype, respectively and a strong linkage disequilibrium was observed in healthy Japanese adults.

Table 3 Distribution of MTX-related enzyme gene and transporter polymorphisms in healthy Japanese adults and patients with rheumatoid arthritis
Full size table

Gender differences in the allelic frequencies are also summarized in Table 3. Allelic frequencies of RFC1 80G > A and RFC1−43 T > C were 61% in men and 52% in women for the A and C alleles and those of ABCC2 IVS23 + 56 T > C were 80% in men and 69% in women for the C allele. These differences in distribution were statistically significant (p = 0.0226 for RFC1 80G > A and RFC1−43 T > C, p = 0.0009 for ABCC2 IVS23 + 56 T > C, chi-square test). No gender differences in the distribution of other genotypes were observed (Table 3), although a tendency for a difference by gender was noted in the distribution of FPGS 1994 G > A (p = 0.0775, chi-square test) and GGH 452 C > T (p = 0.097, chi-square test).

Allelic frequencies and gender differences in frequencies in Japanese RA patients

The distributions of MTHFR 1298A > C, and TYMS 3'-UTR (−6/+6), RFC1 80G > A, RFC1−43 T > C, FPGS 1994G > A, GGH 452C > T, GGH−401C > T, ABCB1 3435C > T, ABCC2 IVS 23 + 56 T > C, and ABCG2 914C > A polymorphisms in 159 Japanese RA patients are summarized in Table 3. The proportions of genotypes at each site for their transporter carrier proteins and enzymes were generally in agreement with Hardy-Weinberg equilibrium, following Mendelian principles. Only three genotype combination patterns occurred for RFC1 80G > A and RFC1−43 T > C, consisting of the 80G/G,−43 T/T genotype, 80G/A,−43 T/C genotype, and 80A/A,−43C/C genotype, respectively and a strong linkage disequilibrium was observed in Japanese RA patients. Allelic frequencies of TYMS 3'-UTR−6 > +6 were 38% in men and 32% in women for the +6 allele, and a significant gender difference in TYMS 3'-UTR−6 > +6 (p = 0.0064, chi-square test) was seen, but the proportions of genotypes in men were not in agreement with Hardy-Weinberg equilibrium (p = 0.003, chi-square test). Gender differences in the distribution of other genotypes were not observed (Table 3), although a tendency toward gender differences was seen in the distribution of MTHFR 1298 A > C (p = 0.0736, chi-square test).

Comparison of allelic frequencies between healthy Japanese adults and RA patients

The distributions of MTHFR 1298A > C, TYMS 3'-UTR (−6/+6), RFC1 80G > A, RFC1−43 T > C, FPGS 1994G > A, GGH 452C > T, GGH−401C > T, ABCB1 3435C > T, ABCC2 IVS 23 + 56 T > C, and ABCG2 914C > A polymorphisms in 299 healthy Japanese healthy adults and 159 Japanese RA patients are summarized in Table 3. There were no statistically significant differences in allelic frequencies between the healthy and RA groups.

Discussion

We comprehensively investigated the allelic frequencies in the gene polymorphisms of enzymes and transporter proteins including MTHFR, TYMS, RFC1, FPGS, GGH, ABCB1, ABCC2 and ABCG2, which affect MTX pharmacokinetics and the therapeutic response to it, between a healthy Japanese population and Japanese RA patients to obtain the fundamental data for the personalized MTX therapy for RA.

Characteristics of MTX-related enzyme gene and transporter polymorphisms are summarized in Table 4. In the comparison of allelic frequencies in ethnic and/or population groups in healthy adults and RA patients, the frequency of the MTHFR1298 C allele in our Japanese study was significantly lower than in Caucasian [4, 7, 8], but higher than in Africans [4]. The MTHFR gene is associated with the generation of 5-methyl THF, the MTHFR 1298A > C polymorphism decreases MTHFR activity, and is thereby associated with MTX efficacy [41, 42] and toxicity [26, 43]. In terms of the effect of each SNP on MTX efficacy and toxicity, it is considered to be the greatest in Caucasians [4, 7, 8], followed by Japanese and Africans [4]. However, currently there are few data comparing Japanese with Africans in relation to the efficacy and toxicity of MTX treatment. As such one example, Hughes et al. [26] reported that there was an association between scores of MTX toxicity and the rs4846051 C allele, that is, a higher mean toxicity score among African-Americans than among Caucasians, and haplotypes containing this allele in African-Americans, but not in Caucasians.

Table 4 Summarized characteristics of MTX-related enzyme gene and transporter polymorphisms
Full size table

The frequency of RFC1 80 A and RFC1−43C alleles in the present Japanese study was found to be higher than in Caucasian healthy adults [4, 7, 8, 11] and RA patients [24, 34–36]. The present study found it to be lower in Japanese than in African healthy adults [4]. RA patients with the RFC1 80 A/A genotype had increased MTX-PG concentrations in RBCs [9, 41], and the mean MTX-PG concentration in RBCs in RA patients with that genotype was reported to be 3.4-fold greater than in other genotypes [44]. Those patients showed a good clinical response to MTX treatment [35]. On the other hand, the−43 T > C change decreased the expression of RFC1 protein in patients with RA [34]. The directions in MTX-PG influx between RFC1 80 G > A and−43 T > C are completely opposite. To clarify the contributions of these two genotypes to MTX-PG concentrations in RBCs, it will be necessary to measure the concentrations in RFC1 80G > A and RFC1−43 T > C RA patients.

The frequency of the FPGS1994 A allele in our Japanese study was significantly lower than reported in Caucasian [4] and Africans [4]. FPGS 1994 A > G polymorphism reported to have no effect on MTX efficacy or toxicity [45], although there was an association between FPGS mRNA expression in peripheral blood mononuclear cells and a poor response to MTX in RA patients with that genotype [46]. Therefore, it is unknown whether the FPGS 1994 A > G polymorphism is associated with the inter-individual difference in MTX efficacy or toxicity.

For the frequency of the GGH 452 T allele, our Japanese study found it to be significantly lower than reported in Caucasian healthy individuals [4] and was the almost same as in Caucasian RA patients [27, 37, 38]. The GGH−401 T allelic frequency in our study was also lower than in HapMap-JPT [4] and in Caucasians [4], but it was higher than in Africans [4]. The GGH gene is involved in reversing polyglutamation by removing glutamate moieties, and a specific functional SNP (452C > T) in the human GGH gene is associated with lower catalytic activity and greater accumulation of long-chain MTX-PGs in leukemia cells of patients treated with high-dose MTX [47]. The−401C > T mutation resulting in greater GGH promotor activity increases the hydrolytic activity of MTX-PGs and decreases polyglutamation in the−401TT genotype compared with the−401CC or CT genotype [44]. Thus, MTX might be more effective in Caucasian than in Japanese RA patients at the same dose, when we speculate on the intracellular MTX-PG concentration.

The frequency of the ABCB1 3435 T allele in our study was lower than reported in Caucasian [4] but higher than in African healthy adults [4]. Conversely, the ABCC2 IVS23 + 56 C allelic frequency was higher in our study than reported in Caucasians [4] but lower than in Africans [4]. We also found a higher frequency of the ABCG2 914 A allele in our study than reported in Caucasians [4] and Africans [4]. The ABCB1 3435C > T mutation resulting in decreased stability and expression of mRNA reduces the activity of efflux transporters, may affect P-glycoprotein function and MTX efflux from cells [48], and is associated with MTX efficacy. It was reported that the length of time before it became necessary to reduce the MTX dose or discontinue administration due to the occurrence of toxicity was 2 months in patients with the T/T genotype, 23 months in those with the T/C genotype, and 29 months in those with the C/C genotype, and the time period was significantly correlated with genotype [40]. Therefore, it is considered that the ABCC2 IVS23 + 56 T > C mutation may affect ABCC2 enzyme activity involved in MTX efflux from cells [40]. The ABCG2 914C > A mutation was reported to increase MTX-PG1 and MTX-PG2 concentrations in RA patients [27]. It is considered that MTX efficacy is the highest in Caucasians, followed by Japanese and Africans, from the viewpoint of ABCB1 3435C > T; conversely, it is the highest in Caucasians, followed by Japanese and Africans, from the viewpoint of ABCC2 IVS 23 + 56 T > C; and the highest in Japanese, followed by Caucasians and then Africans, from the viewpoint of ABCG2 914 C > A. However, we need to investigate the net contribution in MTX-PG efflux among three ABC transporters polymorphism.

There were no previous reports of RFC1−43 T > C, FPGS 1994 G > A, GGH−401 C > T, and ABCC2 IVS23 + 56 T > C in a healthy Japanese population except for HapMap-JPT and RFC1−43 T > C, FPGS 1994 G > A, and ABCC2 IVS23 + 56 T > C in Japanese patients with RA for the evaluation of allelic frequency. Therefore, this is the first comprehensive report on these genetic polymorphisms in healthy Japanese adults and Japanese RA patients.

When comparing healthy Japanese adults with Japanese RA patients, there were no significant differences in allelic frequencies of the MTHFR 1298A > C, TYMS 3'-UTR (−6/+6), RFC1 80G > A, RFC1−43 T > C, FPGS 1994G > A, GGH 452C > T, GGH−401C > T, ABCB1 3435C > T, ABCC2 IVS 23 + 56 T > C, and ABCG2 914C > A polymorphisms in 299 healthy Japanese adults and 159 Japanese RA patients. Therefore, we confirmed that the allelic frequencies of those SNPs in healthy Japanese adults were almost the same as in Japanese patients with RA in this study.

The allelic frequencies of SNPs in Japanese RA patients in this study were in agreement with Hardy-Weinberg equilibrium and were considered representative of the Japanese population as a whole. Chatzikyriakidou et al. [34] reported a strong linkage disequilibrium between RFC1 80 G > A and RFC1−43 T > C genotypes for only three patterns in the combination of the RFC1 80 G > A with RFC1−43 T > C genotype, that is, 80G/G and−43 T/T, 80G/A and−43 T/C, and 80A/A and−43C/C in Greek patients with RA. In agreement with the results of Chatzikyriakidou et al. [34], we also confirmed that the three patterns in the combination of the RFC1 80 G > A with RFC1−43 T > C genotype had a strong linkage disequilibrium between RFC1 80 G > A and RFC1−43 T > C genotypes in healthy Japanese adults and RA patients.

Gender differences in allelic frequencies of RFC1 80G > A, RFC1−43 T > C, and ABCC2 IVS23 + 56 T > C in healthy Japanese adults and TYMS 3'-UTR−6 > +6 in Japanese RA patients were observed. Meanwhile, Kameda et al. [48] reported that the tender joint count, swollen joint count and decrease in serum C-reactive protein (CRP) levels as indices of MTX treatment response in Japanese patients with RA were better in men than in women and gender differences in the MTX therapeutic response were suggested in Japanese RA patients. Currently, however, we cannot clearly explain the reason for the gender differences in these allelic frequencies between healthy Japanese adults and Japanese RA patients, but Hardy-Weinberg disequilibrium in the distribution of those SNPs in either men or women may be one factor responsible for the differences.

This study had several limitations: 1) There were only a small number of reports in the literature allowing comparisons of ethnic and/or population groups, because most did not clearly describe the ethnicity/race and/or population groups studied. These were assumed based on the country of residence of the lead author, and therefore might have misclassified ethnicity/race and/or population groups. 2) We collected allelic frequency data from the HapMap Project for comparisons with our Japanese study results in cases when sufficient data could not be retrieved from the literature search. Although the HapMap data contained the same rs numbers, they included different ss numbers. 3) Several papers used for comparison with our results did not report the rs numbers for ethnic/race and/or population groups. Therefore, our comparisons may not be completely accurate in using the same rs numbers of SNPs. 4) All SNPs studied here could not be compared among ethnic/race and/or population groups due to a lack of data. This study may also include mixed ethnicities/races in a population group that we treated as a single national population, because that was assumed based on the lead author’s country of residence. 5) The distributions of several SNPs in the ethnic/population group in healthy adults and RA patients were also not in agreement with Hardy-Weinberg equilibrium. 6) We did not evaluate allelic frequencies using combinations of SNPs such as diplotypes or haplotypes in our Japanese population. 7) There are few studies of MTX pharmacokinetics among ethnic/population groups. We also did not investigate MTX-PGs pharmacokinetics among ethnic/population groups by comparing the genotypes of MTX-related SNPs to clarify interethnic variations in MTX pharmacokinetics and therapeutic response. This is the greatest limitation in this study.

Conclusion

This study identified allelic frequencies in ethnic and/or population groups compared with healthy Japanese and RA patients. The differences in frequencies may be variables in the interethnic variations in MTX response. Further study is needed to confirm the association of these genetic/phenotypic factors and clinical outcomes and whether these SNPs may be useful in determining personalized MTX therapy for RA.

Abbreviations

3'-UTR:

3'-untranslated region

ABCB1:

ATP binding cassette subfamily B member 1

ABCC2:

ATP-binding cassette sub-family C member 2

ABCG2:

ATP-binding cassette sub-family G member 2

FPGS:

Folypolyglutamyl synthase

GGH:

γ-glutamyl hydrolase

IVS:

Intervening sequence

MTHFR:

Methylenetetrahydrofolate reductase

RFC1:

Reduced folate carrier 1

TYMS:

Thymidylate synthase

References

  1. Japan College of Rheumatology. Guideline for the management of rheumatoid arthritis. Osaka: Medical Review Co., Ltd; 2014.

  2. Ranganathan P. An update on methotrexate pharmacogenetics in rheumatoid arthritis. Pharmacogenomics. 2008;9:439–51.

    Article  CAS  PubMed  Google Scholar 

  3. Inoue S, Hashiguchi M, Chiyoda T, Sunami Y, Tanaka T, Mochizuki M. Pharmacogenetic study of methylenetetrahydrofolate reductase and thymidylate synthase in Japanese and assessment of ethnic and gender differences. Pharmacogenomics. 2007;8:41–7.

    Article  CAS  PubMed  Google Scholar 

  4. National Center for Biotechnology Information SNP database (https://www.ncbi.nlm.nih.gov/snp). Accessed 1 May 2016.

  5. Fukino K, Kawashima T, Suzuki M, Ueno K. Methylenetetrahydrofolate reductase and reduced folate carrier-1 genotypes and methotrexate serum concentrations in patients with rheumatoid arthritis. J Toxicol Sci. 2007;32(4):449–52.

    Article  CAS  PubMed  Google Scholar 

  6. Kumagai K, Hiyama K, Oyama T, Maeda H, Kohno N. Polymorphisms in the thymidylate synthase and methylenetetrahydrofolate reductase genes and sensitivity to the low-dose methotrexate therapy in patients with rheumatoid arthritis. Int J Mol Med. 2003;11(5):593–600.

    CAS  PubMed  Google Scholar 

  7. Herrlinger KR, Cummings JR, Barnardo MC, Schwab M, Ahmad T, Jewell DP. The pharmacogenetics of methotrexate in inflammatory bowel disease. Pharmacogenet Genomics. 2005;15(10):705–11.

    Article  CAS  PubMed  Google Scholar 

  8. Kotsopoulos J, Zhang WW, Zhang S, et al. Polymorphisms in folate metabolizing enzymes and transport proteins and the risk of breast cancer. Breast Cancer Res Treat. 2008;112(3):585–93.

    Article  CAS  PubMed  Google Scholar 

  9. Takatori R, Takahashi KA, Tokunaga D, et al. ABCB1 C3435T polymorphism influences methotrexate sensitivity in rheumatoid arthritis patients. Clin Exp Rheumatol. 2006;24(5):546–54.

    CAS  PubMed  Google Scholar 

  10. Hiraoka M. Folate intake, serum folate, serum total homocysteine levels and methylenetetrahydrofolate reductase C677T polymorphism in young Japanese women. J Nutr Sci Vitaminol. 2004;50(4):238–45.

    Article  CAS  PubMed  Google Scholar 

  11. Summers CM, Mitchell LE, Stanislawska-Sachadyn A, et al. Genetic and lifestyle variables associated with homocysteine concentrations and the distribution of folate derivatives in healthy premenopausal women. Birth Defects Res A Clin Mol Teratol. 2010;88(8):679–88.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Oppeneer SJ, Ross JA, Koh WP, Yuan JM, Robien K. Genetic variation in folylpolyglutamate synthase and gamma-glutamyl hydrolase and plasma homocysteine levels in the Singapore Chinese Health Study. Mol Genet Metab. 2012;105(1):73–8.

    Article  CAS  PubMed  Google Scholar 

  13. Hayashi H, Fujimaki C, Inoue K, Suzuki T, Itoh K. Genetic polymorphism of C452T (T127I) in human gamma-glutamyl hydrolase in a Japanese population. Biol Pharm Bull. 2007;30(4):839–41.

    Article  CAS  PubMed  Google Scholar 

  14. Hamada A, Kato T, Saito H, Mori S. Pharmacogenetic study to develop the personalized medicine of methotrexate monotherapy in patients with rheumatoid arthritis. Recent Adv Clin Pharmacol. 2010;31:154–62.

    Google Scholar 

  15. Komoto C, Nakamura T, Sakaeda T, et al. MDR1 haplotype frequencies in Japanese and Caucasian, and in Japanese patients with colorectal cancer and esophageal cancer. Drug Metab Pharmacokinet. 2006;21(2):126–32.

    Article  CAS  PubMed  Google Scholar 

  16. Honda M, Ogura Y, Toyoda W, et al. Multiple regression analysis of pharmacogenetic variability of carvedilol disposition in 54 healthy Japanese volunteers. Biol Pharm Bull. 2006;29(4):772–28.

    Article  CAS  PubMed  Google Scholar 

  17. Sakaeda T, Nakamura T, Horinouchi M, et al. MDR1 genotype-related pharmacokinetics of digoxin after single oral administration in healthy Japanese subjects. Pharm Res. 2001;18(10):1400–4.

    Article  CAS  PubMed  Google Scholar 

  18. Tanabe M, Ieiri I, Nagata N, et al. Expression of P-glycoprotein in human placenta: relation to genetic polymorphism of the multidrug resistance (MDR)-1 gene. J Pharmacol Exp Ther. 2001;297(3):1137–43.

    CAS  PubMed  Google Scholar 

  19. Yamagishi K, Tanigawa T, Kitamura A, Köttgen A, Folsom AR, Iso H, CIRCS Investigators. The rs2231142 variant of the ABCG2 gene is associated with uric acid levels and gout among Japanese people. Rheumatology. 2010;49(8):1461–5.

    Article  CAS  PubMed  Google Scholar 

  20. Kato T, Hamada A, Mori S, Saito H. Genetic polymorphisms in metabolic and cellular transport pathway of methotrexate impact clinical outcome of methotrexate monotherapy in Japanese patients with rheumatoid arthritis. Drug Metab Pharmacokinet. 2012;27(2):192–9.

    Article  CAS  PubMed  Google Scholar 

  21. Hayashi H, Fujimaki C, Daimon T, Tsuboi S, Matsuyama T, Itoh K. Genetic polymorphisms in folate pathway enzymes as a possible marker for predicting the outcome of methotrexate therapy in Japanese patients with rheumatoid arthritis. J Clin Pharmacol Ther. 2009;34(3):355–61.

    Article  CAS  Google Scholar 

  22. Xiao H, Xu J, Zhou X, et al. Associations between the genetic polymorphisms of MTHFR and outcomes of methotrexate treatment in rheumatoid arthritis. Clin Exp Rheumatol. 2010;28(5):728–33.

    CAS  PubMed  Google Scholar 

  23. Palomino-Morales R, Gonzalez-Juanatey C, Vazquez-Rodriguez TR, et al. A1298C polymorphism in the MTHFR gene predisposes to cardiovascular risk in rheumatoid arthritis. Arthritis Res Ther. 2010;12(2):R71.

    Article  PubMed  PubMed Central  Google Scholar 

  24. James HM, Gillis D, Hissaria P, et al. Common polymorphisms in the folate pathway predict efficacy of combination regimens containing methotrexate and sulfasalazine in early rheumatoid arthritis. J Rheumatol. 2008;35(4):562–71.

    CAS  PubMed  Google Scholar 

  25. Kurzawski M, Pawlik A, Safranow K, Herczynska M, Drozdzik M. 677C > T and 1298A > C MTHFR polymorphisms affect methotrexate treatment outcome in rheumatoid arthritis. Pharmacogenomics. 2007;8(11):1551–9.

    Article  CAS  PubMed  Google Scholar 

  26. Hughes LB, Beasley TM, Patel H, et al. Racial or ethnic differences in allele frequencies of single-nucleotide polymorphisms in the methylenetetrahydrofolate reductase gene and their influence on response to methotrexate in rheumatoid arthritis. Ann Rheum Dis. 2006;65(9):1213–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Stamp LK, Chapman PT, O’Donnell JL, et al. Polymorphisms within the folate pathway predict folate concentrations but are not associated with disease activity in rheumatoid arthritis patients on methotrexate. Pharmacogenet Genomics. 2010;20(6):367–76.

    Article  CAS  PubMed  Google Scholar 

  28. Pawlik A, Kurzawski M, Gawronska-Szklarz B, et al. The effect of 677C > T and 1298A > C MTHFR polymorphisms on sulfasalazine treatment outcome in rheumatoid arthritis. Braz J Med Biol Res. 2009;42(7):660–4.

    Article  CAS  PubMed  Google Scholar 

  29. Lee YC, Cui J, Costenbader KH, Shadick NA, Weinblatt ME, Karlson EW. Investigation of candidate polymorphisms and disease activity in rheumatoid arthritis patients on methotrexate. Rheumatology. 2009;48(6):613–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Pawlik A, Kurzawski M, Górnik W, Dabrowska-Zamojcin E, Droździk M. 677C > T and 1298A > C MTHFR polymorphisms affect arechin treatment outcome in rheumatoid arthritis. Pharmacol Rep. 2007;59(6):721–6.

    CAS  PubMed  Google Scholar 

  31. Berkun Y, Levartovsky D, Rubinow A, et al. Methotrexate related adverse effects in patients with rheumatoid arthritis are associated with the A1298C polymorphism of the MTHFR gene. Ann Rheum Dis. 2004;63(10):1227–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Hayashi H, Tazoe Y, Tsuboi S, et al. A single nucleotide polymorphism of reduced folate carrier 1 predicts methotrexate efficacy in Japanese patients with rheumatoid arthritis. Drug Metab Pharmacokinet. 2013;28(2):164–8.

    Article  CAS  PubMed  Google Scholar 

  33. Ando Y, Shimada H, Matsumoto N, et al. Role of methotrexate polyglutamation and reduced folate carrier 1 (RFC1) gene polymorphisms in clinical assessment indexes. Drug Metab Pharmacokinet. 2013;28(5):442–5.

    Article  CAS  PubMed  Google Scholar 

  34. Chatzikyriakidou A, Georgiou I, Voulgari PV, Papadopoulos CG, Tzavaras T, Drosos AA. Transcription regulatory polymorphism−43 T > C in the 5'-flanking region of SLC19A1 gene could affect rheumatoid arthritis patient response to methotrexate therapy. Rheumatol Int. 2007;27(11):1057–61.

    Article  CAS  PubMed  Google Scholar 

  35. Drozdzik M, Rudas T, Pawlik A, Gornik W, Kurzawski M, Herczynska M. Reduced folate carrier-1 80G > A polymorphism affects methotrexate treatment outcome in rheumatoid arthritis. Pharmacogenomics J. 2007;7(6):404–7.

    Article  CAS  PubMed  Google Scholar 

  36. Lima A, Bernardes M, Azevedo R, et al. SLC19A1, SLC46A1 and SLCO1B1 polymorphisms as predictors of methotrexate-related toxicity in Portuguese rheumatoid arthritis patients. Toxicol Sci. 2014;142(1):196–209.

    Article  CAS  PubMed  Google Scholar 

  37. Wessels JA, van der Kooij SM, le Cessie S, et al. Pharmacogenetics Collaborative Research Group. A clinical pharmacogenetic model to predict the efficacy of methotrexate monotherapy in recent-onset rheumatoid arthritis. Arthritis Rheum. 2007;56(6):1765–75.

    Article  CAS  PubMed  Google Scholar 

  38. Jekic B, Lukovic L, Bunjevacki V, et al. Association of the TYMS 3G/3G genotype with poor response and GGH 354GG genotype with the bone marrow toxicity of the methotrexate in RA patients. Eur J Clin Pharmacol. 2013;69(3):377–83.

    Article  CAS  PubMed  Google Scholar 

  39. Drozdzik M, Rudas T, Pawlik A, et al. The effect of 3435C > T MDR1 gene polymorphism on rheumatoid arthritis treatment with disease-modifying antirheumatic drugs. Eur J Clin Pharmacol. 2006;62(11):933–7.

    Article  CAS  PubMed  Google Scholar 

  40. Ranganathan P, Culverhouse R, Marsh S, et al. Methotrexate (MTX) pathway gene polymorphisms and their effects on MTX toxicity in Caucasian and African American patients with rheumatoid arthritis. J Rheumatol. 2008;35(4):572–9.

    CAS  PubMed  Google Scholar 

  41. Wessels JA, de Vries-Bouwstra JK, Heijmans BT, et al. Efficacy and toxicity of methotrexate in early rheumatoid arthritis are associated with single-nucleotide polymorphisms in genes coding for folate pathway enzymes. Arthritis Rheum. 2006;54(4):1087–95.

    Article  CAS  PubMed  Google Scholar 

  42. Taniguchi A, Urano W, Tanaka E, et al. Validation of the associations between single nucleotide polymorphisms or haplotypes and responses to disease-modifying antirheumatic drugs in patients with rheumatoid arthritis: a proposal for prospective pharmacogenomic study in clinical practice. Pharmacogenet Genomics. 2007;17(6):383–90.

    Article  CAS  PubMed  Google Scholar 

  43. Dervieux T, Greenstein N, Kremer J. Pharmacogenomic and metabolic biomarkers in the folate pathway and their association with methotrexate effects during dosage escalation in rheumatoid arthritis. Arthritis Rheum. 2006;54(10):3095–103.

    Article  CAS  PubMed  Google Scholar 

  44. Dervieux T, Kremer J, Lein DO, et al. Contribution of common polymorphisms in reduced folate carrier and gamma-glutamylhydrolase to methotrexate polyglutamate levels in patients with rheumatoid arthritis. Pharmacogenetics. 2004;14(11):733–9.

    Article  CAS  PubMed  Google Scholar 

  45. van der Straaten RJ, Wessels JA, de Vries-Bouwstra JK, et al. Exploratory analysis of four polymorphisms in human GGH and FPGS genes and their effect in methotrexate-treated rheumatoid arthritis patients. Pharmacogenomics. 2007;8(2):141–50.

    Article  PubMed  Google Scholar 

  46. Stranzl T, Wolf J, Leeb BF, Smolen JS, Pirker R, Filipits M. Expression of folylpolyglutamyl synthetase predicts poor response to methotrexate therapy in patients with rheumatoid arthritis. Clin Exp Rheumatol. 2003;21(1):27–32.

    CAS  PubMed  Google Scholar 

  47. Cheng Q, Wu B, Kager L, et al. A substrate specific functional polymorphism of human gamma-glutamyl hydrolase alters catalytic activity and methotrexate polyglutamate accumulation in acute lymphoblastic leukaemia cells. Pharmacogenetics. 2004;14(8):557–67.

    Article  CAS  PubMed  Google Scholar 

  48. Kameda H, Amano K, Sekiguchi N, et al. Factors predicting the response to low-dose methotrexate therapy in patients with rheumatoid arthritis: a better response in male patients. Mod Rheumatol. 2004;14(6):442–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We are grateful to Dr. Anthoula Chatzikyriakidou, Genetics Unit, Department of Obstetrics and Gynecology, Ioannina University School of Medicine, Ioannina, Greece, for useful information on the strict linkage disequilibrium in RFC-1 80G > A and−43 T > C polymorphisms.

Funding

This study was supported in part by Grants for Encouragement of Research from the Japan Research Foundation for Clinical Pharmacology.

Availability of data and materials

Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study. Please contact the authors for data requests.

Authors’ contribution

MH, KM, JH, OT and MS carried out the genotype determination. TT, TT, MS, TC, JH and SI carried out the clinical study. MH, JH, KM and MS performed the statistical analyses. MH, KM, TT, SI, MS, JM and MM participated in study design and coordination and drafted the manuscript. MH, MS and MM completed the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

The study protocol was approved by the Ethics Committee of the Clinical Pharmacology Center, Sumida Hospital (B-28) and PS Clinic (C-72), and written informed consent was given by all volunteers and patients prior to enrollment.

Author information

Authors and Affiliations

  1. Division for Evaluation and Analysis of Drug Information, Faculty of Pharmacy, Keio University, 1-5-30 Shibakoen, Minato-ku, Tokyo, 105-8512, Japan

    Masayuki Hashiguchi, Jun Hakamata, Kumika Miyawaki, Junya Maruyama & Mayumi Mochizuki

  2. Department of Hygienic Chemistry, Faculty of Pharmacy, Keio University, 1-5-30 Shibakoen, Minato-ku, Tokyo, 105-8512, Japan

    Mikiko Shimizu

  3. PS Clinic, LTA Clinical Pharmacology Center, 6-18 Tenyacho, Hakata-ku, Fukuoka, 812-0025, Japan

    Tomomi Tsuru & Midori Suzaki

  4. Sumida Hospital, LTA Clinical Pharmacology Center, 1-29-1 Honjo, Sumida-ku, Tokyo, 130-0004, Japan

    Takanori Tanaka, Takeshi Chiyoda & Jiro Hiratsuka

  5. Division of Research, BioMedical Laboratory, Kitasato University Kitasato Institute Hospital, 5-9-1 Shirokane, Minato-ku, Tokyo, 108-8642, Japan

    Osamu Takeuchi

  6. Kyushu Clinical Pharmacology Research Clinic, LTA Clinical Pharmacology Center, 6-18 Tenyacho, Hakata-ku, Fukuoka, 812-0025, Japan

    Shin Irie

Authors
  1. Masayuki Hashiguchi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mikiko Shimizu
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jun Hakamata
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Tomomi Tsuru
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Takanori Tanaka
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Midori Suzaki
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Kumika Miyawaki
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Takeshi Chiyoda
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Osamu Takeuchi
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Jiro Hiratsuka
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Shin Irie
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Junya Maruyama
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Mayumi Mochizuki
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Masayuki Hashiguchi.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hashiguchi, M., Shimizu, M., Hakamata, J. et al. Genetic polymorphisms of enzyme proteins and transporters related to methotrexate response and pharmacokinetics in a Japanese population. J Pharm Health Care Sci 2, 35 (2016). https://doi.org/10.1186/s40780-016-0069-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40780-016-0069-0

Keywords

Journal of Pharmaceutical Health Care and Sciences

ISSN: 2055-0294

Contact us