Next Article in Journal
Transcriptomic Profiling Identifies Differentially Expressed Genes in Palbociclib-Resistant ER+ MCF7 Breast Cancer Cells
Previous Article in Journal
Base Editing: The Ever Expanding Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) Tool Kit for Precise Genome Editing in Plants
Previous Article in Special Issue
Transcriptome Analysis Identifies LINC00152 as a Biomarker of Early Relapse and Mortality in Acute Lymphoblastic Leukemia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment

1
Laboratory for Molecular Biomedicine, Institute of Molecular Genetics and Genetic Engineering, University of Belgrade, Vojvode Stepe 444a, 11 000 Belgrade, Serbia
2
Department of Hematology and Oncology, University Children’s Hospital, University of Belgrade, Tirsova 10, 11 000 Belgrade, Serbia
3
Medical Faculty, University of Belgrade, Doktora Subotica 8, 11 000 Belgrade, Serbia
*
Author to whom correspondence should be addressed.
Genes 2020, 11(4), 468; https://doi.org/10.3390/genes11040468
Submission received: 2 March 2020 / Revised: 20 April 2020 / Accepted: 22 April 2020 / Published: 24 April 2020
(This article belongs to the Special Issue Genetics of Childhood Acute Lymphoblastic Leukemia)

Abstract

:
Methotrexate (MTX) is one of the staples of pediatric acute lymphoblastic leukemia (ALL) treatment. MTX targets the folate metabolic pathway (FMP). Abnormal function of the enzymes in FMP, due to genetic aberrations, leads to adverse drug reactions. The aim of this study was to investigate variants in pharmacogenes involved in FMP and their association with MTX pharmacokinetics (MTX elimination profile) and toxicity in the consolidation therapy phase of pediatric ALL patients. Eleven variants in the thymidylate synthetase (TYMS), methylenetetrahydrofolate reductase (MTHFR), dihydrofolate reductase (DHFR), SLC19A1 and SLCO1B genes were analyzed in 148 patients, using PCR- and sequencing-based methodology. For the Serbian and European control groups, data on allele frequency distribution were extracted from in-house and public databases. Our results show that the A allele of SLC19A1 c.80 variant contributes to slow MTX elimination. Additionally, the AA genotype of the same variant is a predictor of MTX-related hepatotoxicity. Patients homozygous for TYMS 6bp deletion were more likely to experience gastrointestinal toxicity. No allele frequency dissimilarity was found for the analyzed variants between Serbian and European populations. Statistical modelling did not show a joint effect of analyzed variants. Our results indicate that SLC19A1 c.80 variant and TYMS 6bp deletion are the most promising pharmacogenomic markers of MTX response in pediatric ALL patients.

1. Introduction

Acute lymphoblastic leukemia (ALL) is a disease characterized by abnormal immature B or T lymphocytes—lymphoblasts. ALL accounts for 25% of childhood cancers. Childhood ALL is treated according to standardized protocols, such as the Berlin-Frankfurt-Munster protocol (BFM) [1]. Despite the existence of efficient treatment protocols for childhood ALL, 75% of patients experience therapy related adverse effects, and 1–3% of patients die from toxicities of therapy [2].
Methotrexate (MTX) is one of the staples of standardized protocols of pediatric ALL treatment. MTX is used in every phase of the BFM protocol: in the remission induction (intrathecal administration), the consolidation and the maintenance therapy phase [3]. In the consolidation phase, patients receive high dose of MTX and after each dose, MTX pharmacokinetics is evaluated [3]. High plasma MTX concentration is associated with MTX therapy related toxicities. It is not expected that improvement of disease outcome and reduction of adverse effects and toxicities will be achieved by eliminating MTX from therapy protocols, but by finding a way to calibrate the dosage of MTX.
MTX is a folate antagonist, and its main therapeutic effect is an inhibition of DNA synthesis. Folate metabolic pathway (FMP) is affected by MTX in a well-defined way (Figure 1). MTX and its polyglutamate forms inhibit key enzymes in folate metabolism, dihydrofolate reductase (DHFR) and thymidylate synthetase (TYMS), by the mechanism of competitive inhibition [4]. The DHFR enzyme replenishes the reserves of the reduced folate form by converting the folic acid and dihydrofolate into tetrahydrofolate (THF). There are several bioactive forms of reduced folate, including 5,10-methylen-THF and 5-methyl-THF, used for thymidylate synthesis (catalyzed by TYMS) and methionine synthesis, for which methylenetetrahydrofolate reductase (MTHFR) enzyme is needed. Methionine is a precursor of S-adenosyl methionine (SAM), which is used for the reactions of methylation, including the methylation of DNA. A lack of THF related to MTX therapy inhibits the numerous processes of DNA synthesis and methylation, thus leading to cell death (therapeutic effect of MTX). Moreover, pharmacokinetics of MTX and other folates is of the essence when it comes to dosing, in order to prevent adverse drug reactions (ADR) [5]. Membrane transporters of solute carrier (SLC) superfamily are important for MTX distribution and clearance. SLC19A1 transporter is responsible for MTX entering the cell, while SLCO1B1 transporter is involved in MTX elimination via bile [6].
Abnormal function of any of the enzymes and transporters in the FMP, due to genetic aberrations, can influence MTX drug efficacy and toxicity and lead to ADR [7]. Pharmacogenomics (PGx) investigates an influence of the individual genetic signature on the response of a patient to treatment [8]. Therefore, pharmacogenomics of MTX needs to be studied in order to determine if it can be implemented in treatment protocols, with the aim to administer the right dose to each pediatric ALL, based on the variants present in the pharmacogenes [9].
The most extensively studied pharmacogenes associated with patients’ response to MTX are DHFR, TYMS, MTHFR, SLC19A1 and SLCO1B1. Numerous variants in these pharmacogenes have been found to be associated with different prognosis and probability of toxic events after MTX administration. Several studies found polymorphisms in DHFR, MTHFR, TYMS, SLC19A1 and SLCO1B1 to be significantly associated with event-free survival (EFS) of ALL patients [10,11,12,13,14]. The toxic events can be gastrointestinal, hepatic, neurological, hematologic, etc. There have been numerous studies in favor of classifying some of variants as biomarkers of MTX induced toxicity [12,13,15,16,17,18,19].
Several studies investigated MTX plasma concentration as a surrogate marker of MTX toxicities in association with variants in MTX pharmacogenes, without taking into account clinical signs of MTX related toxicity [16,18]. Other studies have dealt with laboratory and clinical features of patients to establish severity of MTX toxicity, the majority of which included limited number of ALL patients. Meta-analyses were carried out in order to address limited number of patients in those studies. Meta-analyses which included variants of MTHFR, TYMS and SLC19A1 pharmacogenes [20,21,22,23] could not definitely disregard or establish analyzed variants as MTX pharmacogenomic markers of ALL treatment.
In order to tackle these issues, we aimed to investigated genetic variants in pharmacogenes involved in the FMP, DHFR, TYMS, MTHFR, SLC19A1 and SLCO1B1, and their association with MTX pharmacokinetics and toxicity in the consolidation phase of treatment of pediatric ALL patients in the context of the Serbian population, making this study the first of its kind.

2. Materials and Methods

2.1. Subjects

This study enrolled childhood ALL patients diagnosed between 2003 and 2012 at University children’s hospital, University of Belgrade, Serbia. Patients were treated according to the ALL IC-BFM 2002 or the ALL IC-BFM 2009 protocol (Ethic Code Number 017/6-990/67). Only patients who reached consolidation phase were included in the study. The study was approved by the University children’s hospital Ethic Committee, and was carried out in accordance with the Declaration of Helsinki. The control group comprised of 238 individuals of Serbian origin. Informed consent was obtained from the parents or guardians of each patient, as well as from the control group subjects.

2.2. Pharmacokinetics of MTX in the Consolidation Therapy Phase

In the consolidation phase, patients receive 4 doses of MTX, either 2 g/m2—medium dose MTX (MD-MTX) or 5 g/m2—high dose MTX (HD-MTX) [3]. High risk (HR) patients received HD-MTX, while the others received either MD-MTX (mM protocol) or HD-MTX (M protocol), depending on immunophenotype, risk class and randomization. The estimated risk of each patient (standard, intermediate or high risk) was assessed according to age, white blood cell count, early treatment response and genetic aberrations profile, following BFM protocol guidelines.
MTX plasma concentration was measured only in patients who received HD-MTX, and this data is used to discriminate between normal and slow elimination MTX profile, the latter associated with higher rate of ADR. MTX plasma level that exceed either 1 or 0.4 μmol/L after 36h and 48h, respectively, is indicative of slow MTX elimination, according to the BFM protocol.

2.3. MTX Toxicity in the Consolidation Therapy Phase

Multiple organ toxicity has been evaluated after each MD- and HD-MTX dose, which included: oral mucositis (OM) and gastrointestinal tract (GIT), liver, renal, skin and neurological toxicities, grading from 0 (no toxicity) to 4 (severe toxicity). Toxicity was estimated through clinical examination, laboratory and radiological findings according to Common Toxicity Criteria of the National Cancer Institute (Common Terminology Criteria for Adverse Events, version 5.0), as presented in Table S1.

2.4. Genetic Variants Detection

In this study, variants in TYMS, MTHFR, DHFR, SLC19A1 and SLCO1B were analyzed (Table 1). Variants in TYMS, MTHFR, DHFR and SLC19A1 were detected, as previously described [11]. The region of the SLCO1B1 gene containing the rs4149056 variant was amplified in a total reaction volume of 30 μL. The reaction mix contained 10 pmol of each primer (forward: CAGCCATGAGGAACTATGAGTCC; reverse: CAGAGATCCCAGGGTAAAGCC), 20–100 ng of genomic DNA, 0.5 mM of each dNTP (Fermentas), 1 × PCR reaction buffer, 2.75 mM MgCl2 and 1 U DNA polymerase (FastGene® Taq DNA Polymerase, NIPPON Genetics EUROPE GmbH, Dueren, Germany). The temperature profile of the PCR reactions was starting with the initial activation of DNA polymerase at 95 °C for 5 min, followed by 35 cycles of 30 s denaturation at 95 °C, 30 s annealing at 57 °C and 30 s elongation at 72 °C, ending with a final extension period of 10 min at 72 °C. The PCR fragments were visualized on 2.0% agarose gel, and subsequently analyzed on an ABI PRISM® 3130 DNA analyzer (Applied Biosystems, Foster City, CA, USA).
Variant detection for the control group of Serbian origin was carried out, as described above, on 105 subjects. For 133 control group individuals, their clinical exomes were sequenced employing TruSight One Sequencing Panel on Miseq platform, and the relevant information was extracted using VariantStudio software (version 3, Illumina, San Diego, CA, USA).
For the European control group, data were extracted from 1000 genome project [24] via Ensembl database, which contains the genome data of 503 individuals of European descent. From the Ensembl database, variant frequencies were extracted of all genes of interest, except TYMS. Frequencies of TYMS variants in the European population were estimated from British, American, Italian and Portuguese Caucasians [25,26,27,28,29].

2.5. Statistical Analysis

All detected variants were tested for the Hardy–Weinberg equilibrium using an exact test. Haplotype phases were determined using maximum likelihood algorithm in Arlequin software (version 3.5.1.3, University of Bern, Bern, Switzerland). Differences in allele frequencies between the Serbian and European control groups were tested using chi square test. Genetic variants in transporter genes (SLC19A1 and SLCO1B1) were tested for association with number of low MTX clearance episodes using liner regression. The association of variants in all analyzed genes were tested for association with toxicity using logistic regression. All probabilities were adjusted for demographic and clinical characteristics, as well as MTX dose. Probability values less than 0.05 were considered statistically significant. Multiple comparison correction of all probabilities was applied using false discovery rate (FDR) method. Association analyses were performed using SPSS (version 21, IBM, Armonk, NY, USA) and R (version 3.4.3, R Foundation for Statistical Computing, Vienna, Austria) software.

3. Results

3.1. Population PGx

A total of 9 single nucleotide variants were analyzed in the promoter of DHFR and the coding region of MTHFR, SLC19A1 and SLCO1B1 genes. One variable number of tandem repeats (VNTR) and one indel, both in UTRs of TYMS gene, were also included in analyses. Variant frequencies of childhood ALL patients, as well as control groups of Serbian and European origin, are presented in Table 1. The Hardy–Weinberg equilibrium testing showed that the TYMS VNTR variant (rs34743033) in the ALL group and MTHFR c.677C>T in the control group were not in the equilibrium. However, after FDR correction for multiple testing, none of the analyzed variants showed statistically significant violation of Hardy–Weinberg equilibrium. Next, we tested for differences in allele frequency distribution between the Serbian and European control groups, and we did not find any statistically significant dissimilarity (Table 1).

3.2. Patients

The childhood ALL cohort comprised of 148 patients age 0.9–17.6 years, median 5.5 years, with 94 boys (60.8%). One hundred and eleven patients (75.0%) were treated according to the ALL IC-BFM 2002 protocol, and the others were treated in accordance with the current ALL IC-BFM 2009 protocol. MTX treatment during the consolidation phase did not differ between these two protocols. There were 32 standard risk, 96 intermediate risk and 20 high risk patients. In consolidation, 102, 28 and 18 patients were treated according to the mM, M and HR protocol, respectively. Two patients classified in HR group were treated according to mM and M protocol, because ADR were reported in the previous phase of the therapy.

3.3. MTX Kinetics during Consolidation Therapy Phase

MTX plasma level data in the consolidation phase was available for 41 of 46 patients, who received HD-MTX protocol (28 according to M and 18 according to HR protocol). MTX plasma concentration 36 h and 48 h after each of four HD-MTX doses was used to determine MTX elimination profile (normal or low). Data were available to evaluate MTX elimination profile after 128 HD MTX doses (on average 3.2 per patient). For each patient, MTX elimination profile after each HD-MTX was used to determine number of MTX low clearance episodes (from 0 to 4). In patients for whom MTX elimination profile was not complete after each of the four HD-MTX doses, missing data were estimated taking into account available data for that patient. (For instance, if two normal and one low MTX elimination profile was noted, the number of MTX low clearance episodes would be estimated to be 1/3*4 = 1.33). The number of MTX low clearance episodes was correlated against variants in transporter genes (SLC19A1 and SLCO1B1). Our results show that with each A allele of SLC19A1 c.80 variant, the odds for having one additional slow MTX elimination episode increased by around 40% (linear regression, p = 0.009, β = 0.42 adjusted for age, gender and protocol in consolidation) (Figure 2). The result remained statistically significant after FDR correction. The SLCO1B1 c.521 variant was not associated with MTX elimination (p = 0.85). In order to examine joint effect of SLC19A1 and SLCO1B1 variants in MTX elimination, we included both variants in the liner regression model. Results showed that SLC19A1 remained significant predictor of the MTX elimination profile (p = 0.014, β = 0.42), while SLCO1B1 variant did not contribute to the model (p = 0.94, β = 0.012).

3.4. MTX Toxicity during Consolidation Therapy Phase

After each MD or HD-MTX dose, toxicity was evaluated according to the Common Toxicity Criteria of the National Cancer Institute criteria. HR patients were excluded from toxicity analyses, because they receive aggressive therapy comprising of multiple chemotherapeutics, thus, attributing observed toxicity to MTX would be challenging. Therefore, for toxicity evaluation, only patients who were treated according to mM and M protocol are included (n = 130). For each non-HR patient, after each MD and HD-MTX dose, toxicity of multiple organs was evaluated, and data regarding the highest toxicity grade per organ is presented in Table 2.
Next, we evaluated genetic markers in the folate pathway as predictive markers of toxicity attributed to MD-MTX and HD-MTX therapy. Patients were divided into low (grade 0–1) and high toxicity groups (grade 2–4), and genetic markers were tested for association using a recessive genetic model for variants in TYMS, MTHFR, SLC19A1 and SLCO1B1 genes. Five DHFR variants located less than 400 bp apart and in strong mutual linkage disequilibrium were analyzed as haplotypes. To evaluate the association of a haplotype with a certain type of toxicity, patients with the haplotype of interest were compared to other patients who did not carry that haplotype. The association analyses with genetic markers were carried out for oral mucositis, hepatotoxicity and GIT toxicity, because a considerable number of patients were affected with these kinds of ADR (Table 2). The results showed that non-HR patients with an AA genotype of SLC19A1 c.80G>A variant were almost 10 times more likely to experience grade 2 or higher hepatotoxicity than the non-HR patients with GG or GA genotype (p = 0.011, logistic regression, adjusted for age, gender and MTX dose). Non-HR patients homozygous for TYMS 6bp deletion were more than four times more likely to experience grade 2 or higher GIT toxicity than the non-HR patients with other genotypes (p = 0.020, logistic regression, adjusted for age, gender and MTX dose). Other associations were not statistically significant (Table 3 and Table S2). When probabilities were corrected for multiple comparisons using FDR, no association remained statistically significant.
Next, we analyzed selected genetic variants in the multivariate logistic models with an aim to examine their joint effect on MTX toxicity. We did not find any significant contribution of variants when they were analyzed together (data not shown), other than that already identified in the previous analysis (Table 3). We concluded that these variants do not interact in prediction of MTX toxicity.

4. Discussion

Research efforts in pediatric ALL are focused not only towards etiopathology of the disease and new therapeutic options comprehension, but also towards finding ways to avoid adverse side effects of therapy, therapy toxicity, in a way to efficiently predict them [30]. Our study has contributed by validation of genetic variants, candidates for PGx markers associated with MTX pharmacokinetics and toxicity in pediatric ALL patients in Serbia. Our study validated the SLC19A1 c.80 variant as being associated with delayed methotrexate clearance and the TYMS 6bp deletion variant as a risk of developing gastrointestinal toxicity in pediatric ALL patients.
The SLC19A1 gene is located on chromosome 21q22.3, and codes for the SLC19A1 transporter [31]. The main route of cellular folate intake is by active transport using folate transporters, including SLC19A1. The most represented folate form in circulation, 5-methyl-THF, as well as MTX, enter the cell efficiently via SLC19A1 [31,32]. SLC19A1 c.80 G>A is the most common and the most studied SLC19A1 variant [33]. This variant is located in exon 2 of the SLC19A1 gene, and affects amino acid change p.His27Arg. As the variant position is in the region coding for transmembrane domain of the SLC19A1 protein, it could disrupt the structure itself, influencing the kinetics of the SLC19A1 substrate. Several variants identified in the region coding for the transmembrane domain of SLC19A1 may influence the capacity of SLC19A1 to transport MTX, preventing the transporter from situating itself in a proper position in the cell membrane [34,35,36,37,38,39,40]. An effect of the SLC19A1 c.80 variant to reduce the SLC19A1 MTX transport was investigated in transport-impaired K562 cells transfected with both His27 and Arg27 transporters [41]. Two-fold differences were detected between the Arg27 and His27 transporters in the interaction affinity of different anti-folates. The results suggested that there were at most minor functional differences between Arg27 and His27 transporters in terms of substrate affinities and transport efficiencies. SLC19A1 c.80 G>A variant reduces the SLC19A1 transfer mediated efficacy of both natural folates and MTX, however, the functional significance and mechanism have to be yet fully elucidated [30,42].
Because of the very important physiologic function of the SLC19A1 transporter, its reduced activity could influence the onset of certain diseases associated with low intracellular folate intake such as malignant, cardiovascular and neurological diseases, as well as modulate the efficacy and toxicity of MTX. Studies examining the association between the genetic variant SLC19A1 c.80A in children with ALL and disease outcome have yielded conflicting results. SLC19A1 c.80 AA and GA genotypes are associated with ALL relapse and poor survival [43]. One study reported that of all the genetic variations associated with the folate pathway, the presence of the SLC19A1 c.80AA genotype conveys the largest risk for developing ALL [44], while another study showed the opposite results: a protective effect of SLC19A1 c.80 AA genotype [45]. A study found an association between SLC19A1 c.80 AA genotype and a higher probability of staying in remission, compared with GG or GA variants, and a lower likelihood of relapse among 500 children with ALL, as well as its effect on greater bone marrow toxicity [46]. Furthermore, the relationship between a greater number of copies of the gene present in ALL patients with trisomy 21 or hyperdiploid ALL, and a better outcome or greater MTX toxicity, was shown in the same study, possibly due to greater intracellular polyglutamated MTX in this setting.
Our results in pediatric ALL patients receiving HD-MTX doses have shown that the possibility for having one additional low clearance episode increases by 40% with each A allele of SLC19A1.c80 variant. This finding suggests that SLC19A1.c80 variant is a strong candidate for MTX kinetics PGx marker.
Additionally, our results demonstrated that the SLC19A1 c.80G>A variant could be a PGx marker for hepatotoxicity in pediatric ALL patients. According to a meta-analysis, it was demonstrated that G80A, a variant which affects the activity of SLC19A1, is not a good marker of MTX-related toxicity in pediatric ALL [33]. Another study concluded that the same variant is associated with increased MTX toxicity [47]. This is another example of the inconsistency of the assessment of relevance of PGx markers for MTX toxicity.
We have also investigated the association of MTX kinetics during the consolidation with variants in SLCO1B1, another transporter encoding gene. MTX clearance was associated with numerous variants in SLCO1B1 [16]. The functional variant c.521 in SLCO1B1 gene has been suggested to play an important role in methotrexate clearance and related toxicities in patients with pediatric ALL in several studies [15,17,18,48]. However, some studies did not replicate those results [49,50,51].
Our results, pointing out that the SLCO1B1 c.521 variant was not associated with MTX elimination, do not support the proposed role of the c.521 variant in decreased MTX clearance and the development of related toxicities. We hypothesize that the lack of consistency across the studies could be due to the fact that the pharmacokinetics of methotrexate can be influenced by both non-genetic and genetic factors, since it has been shown that all SLCO1B1 variants account for 10.7% of the population variability in MTX clearance [17]. Among the non-genetic influences, the differences in treatment protocols, estimations of methotrexate toxicity and homogeneity of the studied populations could be the main sources of differences across the studies.
TYMS encodes the enzyme essential for DNA replication and repair, converting dUMP (uridine monophosphate) into dTMP (thymidine monophosphate). As the function of TYMS enzyme in fast-dividing cells is of crucial importance, the TYMS enzyme is the target in the treatment of various types of malignant and inflammatory diseases. TYMS enzyme inhibitors include pyrimidine analogues, fluoropyrimidines and folate analogs, including MTX [52,53]. The most studied genetic variations in the TYMS gene are located in the 5′UTR and 3′UTR regions, and influence mRNA stability, localization and translational efficiency [54]. A minisatellite region, usually containing three or two repeats, each 28 base pairs (bp) long, is located in the 5′UTR region of the TYMS gene. The allele with two repeats (2R) has been shown to cause decreased translation of TYMS relative to the three-repeat (3R) allele [55], but the exact mechanism is not fully known. It has been reported that patients with rheumatoid arthritis and TYMS 3R3R genotype show a weaker response to MTX therapy [56]. The 3′UTR region of the TYMS gene contains and 6bp-long indel. Alleles with deletion exhibit decreased mRNA stability and result in lower levels of thymidylate synthetase, due to higher degradation susceptibility [57]. Polyglutamate forms of MTX bind directly to the TYMS enzyme and thus inhibit it [58]. This mechanism is related to the observed finding of increased gastrointestinal toxicity. Intestinal mucositis, bleeding and peptic ulcers are well-known gastrointestinal adverse effects of MTX. The pathogenic mechanism underlying this mucositis remains unknown, but immune response and increased blood concentration could be considered as the causes [59]. Polyglutamate forms of MTX bind directly to the TYMS enzyme, and thus inhibit it. If the cell is naturally enriched in the enzyme, due to the presence of the TYMS allele causing increased gene expression, less toxicity but also less effectiveness of therapy is expected [60]. Conversely, patients harboring the TYMS gene variants causing decreased gene expression, display a higher likelihood of adverse drug reactions [61]. We believe that our finding, increased gastrointestinal toxicity in ALL patient carriers of TYMS 6bp deletion, is a result of the interaction of MTX and TYMS, whose low expression level enables more effective action of MTX, leading both to a better therapeutic effect and more prominent toxicity.
Our results have demonstrated that among children with ALL, homozygous carriers TYMS 6bp deletion, were more than four times more likely to experience grade 2 or higher GIT toxicity than the non-HR patients with other genotypes. In addition, non-HR patients with the AA genotype of the SLC19A1 c.80G>A variant were almost 10 times more likely to experience grade 2 or higher hepatotoxicity than the non-HR patients with GG or GA genotype. It was already shown that SLC19A1 c.80A allele was associated with gastrointestinal MTX toxicity in children with ALL [62]. Our findings are in line with the published literature, suggesting that the TYMS and SLC19A1 variants could be explored as pharmacogenetics variants.
We have tried to highlight our results in population-specific genetic background. It has been previously shown that the allele frequencies of several pharmacogenes are significantly different when the Serbian population is compared with other European populations [63]. However, for the variants in pharmacogenes analyzed in this study, no statistically significant dissimilarity was found. Therefore, the population-specific genetic background of the patients included in this study did not contribute to particular findings reported here. The fact that pharmacogenetics variants mostly are not cause of the disease is furthermore seen from the variant frequencies around 0.5 for TYMS rs34743033 and SLC19A1 rs1051266. The majority of variants analyzed in our study had a variant frequency of around 30% (TYMS rs34489327, MTHFR rs1801133 and rs1801131, DHFR rs442767, rs1643641, rs1650695 and rs1650696). These relatively high numbers could indicate the importance of the relevant pharmacogenetics markers in European populations, including the Serbian population.
Some limitations should be considered when interpreting the results of this study. They include relatively small number of pediatric ALL patients—only 130 for toxicity study and 41 for methotrexate clearance follow-up. The candidate gene approach we used could not reveal new potential pharmacogenes but only assess already proposed ones. The study did not directly assess gene expression nor protein levels in the FMP, as these pharmacotranscritomics markers are recognized as important factors in methotrexate clearance and toxicity [11,57,64]. However, our findings contribute to the knowledge of PGx of MTX. Most previous studies have analyzed the association of particular genetic variants in FMP with MTX plasma levels or MTX resistance. The findings of our study that the SLC19A1c.80 variant had significant role in methotrexate clearance and that TYMS 6bp variant could be considered a valid pharmacogenetics marker of MTX gastrointestinal toxicity, which represents an important piece in a complex MTX PGx puzzle.
MTX toxicity has been recognized as potentially life-threatening event during chemotherapy for a while. International guides that identify and grade MTX toxicity [65], already integrated into modern treatment protocols for ALL in children [66], still need evidence-based revision and continuous upgrade. It is demonstrated again that the pharmacogenomics of MTX is a complex research field, since the MTX/folate metabolic pathway is not a simple one [7,58,67]. A number of polymorphic genes encode enzymes and transporters involved in different pathway branches of folate metabolism. Moreover, gene–gene interactions in folate-related genes were observed, and could contribute to complex interdependencies between the pathway participants [68]. Additionally, there are gene-nutrient interactions and the evident modifying role of dietary components on the sophisticated folate machinery [69].
It appears that genome wide approach has not been successful in defining valid PGx markers due to limited statistical power. A candidate gene approach is incapable of identifying new genes and genetic variants that could be used as drug targets or predictive markers [70], because they cannot include all the relevant aspects of MTX metabolism. Nevertheless, continuous technology progress provides us with more genomic and clinical data that area likely to help researchers understand and predict adverse reactions to therapy, as well as contributing to better outcomes of the disease. Molecular profiling of patients using high-throughput technology, as well as the implementation of artificial intelligence through advanced bioinformatics tools and predictive machine learning algorithms, would enable us to create reliable genetic panels for routine clinical pharmacogenomics testing, to be used before administering different cytotoxic drugs, combinations of drugs or drugs used in different phases of therapeutic protocols for pediatric ALL patients [2]. It is a long way ahead to reach personalized treatment in everyday clinical practice, but pediatric ALL is one of the most encouraging examples that it is achievable.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/4/468/s1, Table S1: Criteria for grading drug related toxicity according to Common Toxicity Criteria of the National Cancer Institute. Table S2: Childhood ALL patients (n = 130) experiencing low (grade 0–1) or high (grade 2–4) toxicity divided by presence of specific genotype/haplotype and protocol regiment during the consolidation phase of BFM guided therapy.

Author Contributions

Conceptualization: N.K. (Nikola Kotur), J.L., N.K. (Nada Krstovski), B.Z. and S.P; Formal Laboratory Analysis- Investigation N.K., B.R., V.G. and J.L.; Data Curation: N.K. (Nikola Kotur), B.S., J.L., L.D., G.M., N.K. (Nada Krstovski), D.J., B.Z. and S.P.; Writing—Original Draft Preparation, N.K. (Nikola Kotur), B.S., V.G., B.Z., J.L. and S.P.; Writing—Review & Editing, N.K. (Nikola Kotur), B.S., V.G., B.R., J.L., N.K. (Nada Krstovski), L.D., G.M., D.J., B.Z. and S.P.; Supervision: S.P. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This work was supported by Ministry of Education, Science and Technological Development Republic of Serbia, EB: 451-03-68/2020-14/ 200042.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schrappe, M.; Reiter, A.; Ludwig, W.-D.; Harbott, J.; Zimmermann, M.; Hiddemann, W.; Niemeyer, C.; Henze, G.; Feldges, A.; Zintl, F.; et al. Improved outcome in childhood acute lymphoblastic leukemia despite reduced use of anthracyclines and cranial radiotherapy: Results of trial ALL-BFM 90. German-Austrian-Swiss ALL-BFM Study Group. Blood 2000, 95, 3310–3322. [Google Scholar] [PubMed]
  2. Pavlovic, S.; Kotur, N.; Stankovic, B.; Zukic, B.; Gasic, V.; Dokmanovic, L. Pharmacogenomic and Pharmacotranscriptomic Profiling of Childhood Acute Lymphoblastic Leukemia: Paving the Way to Personalized Treatment. Genes 2019, 10, 191. [Google Scholar] [CrossRef] [Green Version]
  3. ALL IC-BFM 2009. A Randomized Trial of the I-BFM-SG for the Management of Childhood non-B Acute Lymphoblastic Leukemia Final Version of Therapy Protocol from August-14-2009. Available online: http://www.bialaczka.org/wp-content/uploads/2016/10/ALLIC_BFM_2009.pdf (accessed on 28 February 2020).
  4. Farber, S.; Diamond, L.K.; Mercer, R.D.; Sylvester, R.F.; Wolff, J.A. Temporary Remissions in Acute Leukemia in Children Produced by Folic Acid Antagonist, 4-Aminopteroyl-Glutamic Acid (Aminopterin). N. Engl. J. Med. 1948, 238, 787–793. [Google Scholar] [CrossRef] [PubMed]
  5. Lopez-Lopez, E.; Martin-Guerrero, I.; Ballesteros, J.; Piñán, M.; Garcia-Miguel, P.; Navajas, A.; Garcia-Orad, A. Polymorphisms of the SLCO1B1 Gene Predict Methotrexate-Related Toxicity in Childhood Acute Lymphoblastic Leukemia. Pediatr. Blood Cancer 2011, 57, 612–619. [Google Scholar] [CrossRef] [PubMed]
  6. Lopez-Lopez, E.; Gutierrez-Camino, A.; Bilbao-Aldaiturriaga, N.; Pombar-Gomez, M.; Martin-Guerrero, I.; Garcia-Orad, A. Pharmacogenetics of childhood acute lymphoblastic leukemia. Pharmacogenomics 2014, 15, 1383–1398. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Assaraf, Y. The role of multidrug resistance efflux transporters in antifolate resistance and folate homeostasis. Drug Resist. Updat. Rev. Comment. Antimicrob. Anticancer Chemother. 2006, 9, 227–246. [Google Scholar] [CrossRef]
  8. Wang, L.; Weinshilboum, R. Thiopurine S-methyltransferase pharmacogenetics: Insights, challenges and future directions. Oncogene 2006, 25, 1629–1638. [Google Scholar] [CrossRef] [Green Version]
  9. Rudin, S.; Marable, M.; Huang, R. The Promise of Pharmacogenomics in Reducing Toxicity During Acute Lymphoblastic Leukemia Maintenance Treatment. Genom. Proteom. Bioinform. 2017, 15, 82–93. [Google Scholar] [CrossRef]
  10. Ceppi, F.; Gagné, V.; Douyon, L.; Quintin, C.; Colombini, A.; Parasole, R.; Buldini, B.; Basso, G.; Conter, V.; Cazzaniga, G.; et al. DNA variants in DHFR gene and response to treatment in children with childhood B ALL: Revisited in AIEOP-BFM protocol. Pharmacogenomics 2018, 19, 105–112. [Google Scholar] [CrossRef]
  11. Dulucq, S.; St-Onge, G.; Gagne, V.; Ansari, M.; Sinnett, D.; Labuda, D.; Moghrabi, A.; Krajinovic, M. DNA variants in the dihydrofolate reductase gene and outcome in childhood ALL. Blood 2008, 111, 3692–3700. [Google Scholar] [CrossRef]
  12. Krajinovic, M.; Costea, I.; Primeau, M.; Dulucq, S.; Moghrabi, A. Combining several polymorphisms of thymidylate synthase gene for pharmacogenetic analysis. Pharmacogenom. J. 2005, 5, 374–380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Liu, S.-G.; Gao, C.; Zhang, R.-D.; Zhao, X.-X.; Cui, L.; Li, W.-J.; Chen, Z.-P.; Yue, Z.-X.; Zhang, Y.-Y.; Wu, M.-Y.; et al. Polymorphisms in methotrexate transporters and their relationship to plasma methotrexate levels, toxicity of high-dose methotrexate, and outcome of pediatric acute lymphoblastic leukemia. Oncotarget 2017, 8, 37761–37772. [Google Scholar] [CrossRef]
  14. Sepe, D.; McWilliams, T.; Chen, J.; Kershenbaum, A.; Zhao, H.; La, M.; Devidas, M.; Lange, B.; Rebbeck, T.; Aplenc, R. Germline genetic variation and treatment response on CCG-1891. Pediatr. Blood Cancer 2012, 58, 695–700. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Radtke, S.; Zolk, O.; Renner, B.; Paulides, M.; Zimmermann, M.; Moricke, A.; Stanulla, M.; Schrappe, M.; Langer, T. Germline genetic variations in methotrexate candidate genes are associated with pharmacokinetics, toxicity, and outcome in childhood acute lymphoblastic leukemia. Blood 2013, 121, 5145–5153. [Google Scholar] [CrossRef] [Green Version]
  16. Ramsey, L.; Panetta, J.; Smith, C.; Yang, W.; Fan, Y.; Winick, N.; Martin, P.; Cheng, C.; Devidas, M.; Pui, C.-H.; et al. Genome-wide study of methotrexate clearance replicates SLCO1B1. Blood 2012, 121. [Google Scholar] [CrossRef] [PubMed]
  17. Ramsey, L.B.; Bruun, G.H.; Yang, W.; Trevino, L.R.; Vattathil, S.; Scheet, P.; Cheng, C.; Rosner, G.L.; Giacomini, K.M.; Fan, Y.; et al. Rare versus common variants in pharmacogenetics: SLCO1B1 variation and methotrexate disposition. Genome Res. 2012, 22, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Trevino, L.R.; Shimasaki, N.; Yang, W.; Panetta, J.C.; Cheng, C.; Pei, D.; Chan, D.; Sparreboom, A.; Giacomini, K.M.; Pui, C.H.; et al. Germline genetic variation in an organic anion transporter polypeptide associated with methotrexate pharmacokinetics and clinical effects. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2009, 27, 5972–5978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Zaïr, Z.M.; Singer, D.R.J. Influx transporter variants as predictors of cancer chemotherapy-induced toxicity: Systematic review and meta-analysis. Pharmacogenomics 2016, 17, 1189–1205. [Google Scholar] [CrossRef]
  20. He, H.R.; Liu, P.; He, G.H.; Dong, W.H.; Wang, M.Y.; Dong, Y.L.; Lu, J. Association between reduced folate carrier G80A polymorphism and methotrexate toxicity in childhood acute lymphoblastic leukemia: A meta-analysis. Leuk. Lymphoma 2014, 55, 2793–2800. [Google Scholar] [CrossRef]
  21. Lopez-Lopez, E.; Martin-Guerrero, I.; Ballesteros, J.; Garcia-Orad, A. A systematic review and meta-analysis of MTHFR polymorphisms in methotrexate toxicity prediction in pediatric acute lymphoblastic leukemia. Pharmacogenom. J. 2013, 13, 498–506. [Google Scholar] [CrossRef] [Green Version]
  22. Oosterom, N.; Berrevoets, M.; den Hoed, M.A.H.; Zolk, O.; Hoerning, S.; Pluijm, S.M.F.; Pieters, R.; de Jonge, R.; Tissing, W.J.E.; van den Heuvel-Eibrink, M.M.; et al. The role of genetic polymorphisms in the thymidylate synthase (TYMS) gene in methotrexate-induced oral mucositis in children with acute lymphoblastic leukemia. Pharmacogenet. Genom. 2018, 28, 223–229. [Google Scholar] [CrossRef] [PubMed]
  23. Zhu, C.; Liu, Y.W.; Wang, S.Z.; Li, X.L.; Nie, X.L.; Yu, X.T.; Zhao, L.B.; Wang, X.L. Associations between the C677T and A1298C polymorphisms of MTHFR and the toxicity of methotrexate in childhood malignancies: A meta-analysis. Pharmacogenom. J. 2018, 18, 450–459. [Google Scholar] [CrossRef] [PubMed]
  24. 1000 Genomes Project Consortium; Auton, A.; Brooks, L.D.; Durbin, R.M.; Garrison, E.P.; Kang, H.M.; Korbel, j.; Marchini, J.L.; McCarthy, S.; McVean, G.A.; et al. A global reference for human genetic variation. Nature 2015, 526, 68–74. [Google Scholar] [PubMed] [Green Version]
  25. Carmona, B.; Guerreiro, C.; Cravo, M.; Nobre-Leitao, C.; Brito, M. 5′ and 3′ UTR thymidylate synthase polymorphisms modulate the risk of colorectal cancer independently of the intake of methyl group donors. Mol. Med. Rep. 2008, 1, 747–752. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Graziano, F.; Ruzzo, A.; Loupakis, F.; Santini, D.; Catalano, V.; Canestrari, E.; Maltese, P.; Bisonni, R.; Fornaro, L.; Baldi, G.; et al. Liver-only metastatic colorectal cancer patients and thymidylate synthase polymorphisms for predicting response to 5-fluorouracil-based chemotherapy. Br. J. Cancer 2008, 99, 716–721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Hubner, R.; Liu, J.-F.; Sellick, G.; Logan, R.; Houlston, R.; Muir, K. Thymidylate synthase polymorphisms, folate and B-vitamin intake, and risk of colorectal adenoma. Br. J. Cancer 2007, 97, 1449–1456. [Google Scholar] [CrossRef] [Green Version]
  28. Shi, Q.; Zhang, Z.; Neumann, A.; Li, G.; Spitz, M.; Wei, Q. Case–control analysis of thymidylate synthase polymorphism and risk of lung cancer. Carcinogenesis 2005, 26, 649–656. [Google Scholar] [CrossRef]
  29. Zhang, Z.; Shi, Q.; Sturgis, E.; Spitz, M.; Hong, W.; Wei, Q. Thymidylate synthase 5′- and 3′-untranslated region polymorphisms associated with risk and progression of squamous cell carcinoma of the head and neck. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2005, 10, 7903–7910. [Google Scholar] [CrossRef] [Green Version]
  30. Cheok, M.; Evans, W. Acute lymphoblastic leukemia: A model for the pharmacogenomics of cancer therapy. Nat. Rev. Cancer 2006, 6, 117–129. [Google Scholar] [CrossRef]
  31. Hou, Z.; Matherly, L. Biology of the major facilitative folate transporters SLC19A1 and SLC46A1. Curr. Top. Membr. 2014, 73, 175–204. [Google Scholar] [CrossRef] [Green Version]
  32. Matherly, L.; Hou, Z.; Deng, Y. Human reduced folate carrier: Translation of basic biology to cancer etiology and therapy. Cancer Metastasis Rev. 2007, 26, 111–128. [Google Scholar] [CrossRef] [PubMed]
  33. Rady, P.L.; Szucs, S.; Matalon, R.K.; Grady, J.; Hudnall, S.D.; Kellner, L.H.; Nitowsky, H. Genetic Polymorphism (G80A) of Reduced Folate Carrier Gene in Ethnic Populations. Mol. Genet. Metab. 2001, 73, 285–286. [Google Scholar] [CrossRef] [PubMed]
  34. Drori, S.; Jansen, G.; Mauritz, R.; Peters, G.; Assaraf, Y. Clustering of Mutations in the First Transmembrane Domain of the Human Reduced Folate Carrier in GW1843U89-resistant Leukemia Cells with Impaired Antifolate Transport and Augmented Folate Uptake. J. Biol. Chem. 2000, 275, 30855–30863. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Gifford, A.; Haber, M.; Witt, T.; Whetstine, J.; Taub, J.; Matherly, L.; Norris, M. Role of the E45K-reduced folate carrier gene mutation in methotrexate resistance in human leukemia cells. Leuk. Off. J. Leuk. Soc. Am. Leuk. Res. Fund U.K 2003, 16, 2379–2387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Jansen, G.; Mauritz, R.; Drori, S.; Sprecher, H.; Kathmann, I.; Bunni, M.; Priest, D.G.; Noordhuis, P.; Schornagel, J.H.; Pinedo, H.M.; et al. A Structurally Altered Human Reduced Folate Carrier with Increased Folic Acid Transport Mediates a Novel Mechanism of Antifolate Resistance. J. Biol. Chem. 1998, 273, 30189–30198. [Google Scholar] [CrossRef] [Green Version]
  37. Tse, A.; Brigle, K.; Taylor, S.M.; Moran, R.G. Mutations in the Reduced Folate Carrier Gene Which Confer Dominant Resistance to 5,10-Dideazatetrahydrofolate. J. Biol. Chem. 1998, 273, 25953–25960. [Google Scholar] [CrossRef] [Green Version]
  38. Zhao, R.; Assaraf, Y.; Goldman, I. A Mutated Murine Reduced Folate Carrier (RFC1) with Increased Affinity for Folic Acid, Decreased Affinity for Methotrexate, and an Obligatory Anion Requirement for Transport Function. J. Biol. Chem. 1998, 273, 19065–19071. [Google Scholar] [CrossRef] [Green Version]
  39. Zhao, R.; Assaraf, Y.; Goldman, I.D. A reduced folate carrier mutation produces substrate-dependent alterations in carrier mobility in murine leukemia cells and methotrexate resistance with conservation of growth in 5-formyltetrahydrofolate. J. Biol. Chem. 1998, 273, 7873–7879. [Google Scholar] [CrossRef] [Green Version]
  40. Zhao, R.; Sharina, I.; Goldman, I. Pattern of Mutations that Results in Loss of Reduced Folate Carrier Function under Antifolate Selective Pressure Augmented by Chemical Mutagenesis. Mol. Pharmacol. 1999, 56, 68–76. [Google Scholar] [CrossRef] [Green Version]
  41. Whetstine, J.R.; Gifford, A.J.; Witt, T.; Liu, X.Y.; Flatley, R.M.; Norris, M.; Haber, M.; Taub, J.W.; Ravindranath, Y.; Matherly, L.H. Single Nucleotide Polymorphisms in the Human Reduced Folate Carrier. Clin. Cancer Res. 2001, 7, 3416–3422. [Google Scholar]
  42. Galbiatti-Dias, A.L.; Ruiz, M.; Pinto, D.; Raposo, L.; Maníglia, J.; Pavarino, É.; Goloni-Bertollo, E. A80G polymorphism of reduced folate carrier 1 (RFC1) gene and head and neck squamous cell carcinoma etiology in Brazilian population. Mol. Biol. Rep. 2011, 38, 1071–1078. [Google Scholar] [CrossRef]
  43. Leyva-Vazquez, M.; Organista-Nava, J.; Gómez-Gómez, Y.; Contreras-Quiroz, A.; Flores-Alfaro, E.; Illades-Aguiar, B. Polymorphism G80A in the Reduced Folate Carrier Gene and its Relationship to Survival and Risk of Relapse in Acute Lymphoblastic Leukemia. J. Investig. Med. 2012, 60, 1064–1067. [Google Scholar] [CrossRef]
  44. Jonge, R.; Tissing, W.; Hooijberg, J.; Jansen, G.; Kaspers, G.; Lindemans, J.; Peters, G.; Pieters, R. Polymorphisms in folate-related genes and risk of pediatric acute lymphoblastic leukemia. Blood 2008, 113, 2284–2289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Silva, R.; Fontes, A.; Silva, K.; Sant’Ana, T.; Ramos, F.; Marques-Salles, T.; Pombo-de-Oliveira, M.; Muniz, M. Polymorphisms Involved in Folate Metabolism Pathways and the Risk of the Development of Childhood Acute Leukemia. Genet. Test. Mol. Biomark. 2013, 17, 147–152. [Google Scholar] [CrossRef] [PubMed]
  46. Gregers, J.; Christensen, I.; Dalhoff, K.; Lausen, B.; Schroeder, H.; Rosthøj, S.; Carlsen, N.; Schmiegelow, K.; Peterson, C. The association of reduced folate carrier 80G>A polymorphism to outcome in childhood acute lymphoblastic leukemia interacts with chromosome 21 copy number. Blood 2010, 115, 4671–4677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Giletti, A.; Esperon, P. Genetic markers in methotrexate treatments. Pharmacogenom. J. 2018, 18, 689–703. [Google Scholar] [CrossRef] [PubMed]
  48. Csordas, K.; Lautner-Csorba, O.; Semsei, A.; Harnos, A.; Md, M.; Erdélyi, D.; Eipel, O.; Szalai, C.; Kovacs, G. Associations of novel genetic variations in the folate-related and ARID5B genes with the pharmacokinetics and toxicity of high-dose methotrexate in paediatric acute lymphoblastic leukaemia. Br. J. Haematol. 2014, 166, 410–420. [Google Scholar] [CrossRef] [Green Version]
  49. Fukushima, H.; Fukushima, T.; Sakai, A.; Suzuki, R.; Nakajima-Yamaguchi, R.; Kobayashi, C.; Iwabuchi, A.; Saito, M.; Yoshimi, A.; Nakao, T.; et al. Polymorphisms of MTHFR Associated with Higher Relapse/Death Ratio and Delayed Weekly MTX Administration in Pediatric Lymphoid Malignancies. Leuk. Res. Treat. 2013, 2013, 238528. [Google Scholar] [CrossRef]
  50. Giletti, A.; Vital, M.; Lorenzo, M.; Cardozo, P.; Borelli, G.; Gabus, R.; Martínez, L.; Díaz, L.; Assar, R.; Rodriguez, M.N.; et al. Methotrexate pharmacogenetics in Uruguayan adults with hematological malignant diseases. Eur. J. Pharm. Sci. 2017, 109, 480–485. [Google Scholar] [CrossRef]
  51. Lopez-Lopez, E.; Ballesteros, J.; Piñán, M.; Toledo, J.; García de Andoin, N.; Garcia-Miguel, P.; Navajas, A.; Garcia-Orad, A. Polymorphisms in the methotrexate transport pathway: A new tool for MTX plasma level prediction in pediatric acute lymphoblastic leukemia. Pharmacogenet. Genom. 2013, 23, 53–61. [Google Scholar] [CrossRef]
  52. Lima, A.; Azevedo, R.; Sousa, H.; Seabra, V.; Medeiros, R. Current approaches for TYMS polymorphisms and their importance in molecular epidemiology and pharmacogenetics. Pharmacogenomics 2013, 14, 1337–1351. [Google Scholar] [CrossRef] [PubMed]
  53. Lurje, G.; Manegold, P.; Ning, Y.; Pohl, A.; Zhang, W.; Lenz, H.-J. Thymidylate synthase gene variations: Predictive and prognostic markers. Mol. Cancer Ther. 2009, 8, 1000–1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lynch, M.; Scofield, D.; Hong, X. The Evolution of Transcription-Initiation Sites. Mol. Biol. Evol. 2005, 22, 1137–1146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kawakami, K.; Salonga, D.; Park, J.M.; Danenberg, K.D.; Uetake, H.; Brabender, J.; Omura, K.; Watanabe, G.; Danenberg, P.V. Different Lengths of a Polymorphic Repeat Sequence in the Thymidylate Synthase Gene Affect Translational Efficiency but Not Its Gene Expression. Clin. Cancer Res. 2001, 7, 4096–4101. [Google Scholar] [PubMed]
  56. Lima, A.; Seabra, V.; Bernardes, M.; Azevedo, R.; Sousa, H.; Medeiros, R. Role of Key TYMS Polymorphisms on Methotrexate Therapeutic Outcome in Portuguese Rheumatoid Arthritis Patients. PLoS ONE 2014, 9, e108165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Pullmann, R.; Abdelmohsen, K.; Lal, A.; Martindale, J.L.; Ladner, R.D.; Gorospe, M. Differential Stability of Thymidylate Synthase 3′-Untranslated Region Polymorphic Variants Regulated by AUF1. J. Biol. Chem. 2006, 281, 23456–23463. [Google Scholar]
  58. Allegra, C.J.; Chabner, B.A.; Drake, J.C.; Lutz, R.; Rodbard, D.; Jolivet, J. Enhanced inhibition of thymidylate synthase by methotrexate polyglutamates. J. Biol. Chem. 1985, 260, 9720–9726. [Google Scholar]
  59. Tsukada, T.; Nakano, T.; Miyata, T.; Sasaki, S. Life-Threatening Gastrointestinal Mucosal Necrosis during Methotrexate Treatment for Rheumatoid Arthritis. Case Rep. Gastroenterol. 2013, 7, 470–475. [Google Scholar] [CrossRef]
  60. Pullarkat, S.; Stoehlmacher, J.; Ghaderi, V.; Xiong, Y.P.; Ingles, S.; Sherrod, A.; Warren, R.; Tsao-Wei, D.; Groshen, S.; Lenz, H. Thymidylate synthase gene polymorphism determines response and toxicity of 5-FU chemotherapy. Pharmacogenom. J. 2001, 1, 65–70. [Google Scholar] [CrossRef] [Green Version]
  61. Wang, B.; Walsh, S.; Saif, W. Pancytopenia and Severe Gastrointestinal Toxicities Associated with 5-Fluorouracil in a Patient with Thymidylate Synthase (TYMS) Polymorphism. Cureus 2016, 8, e798. [Google Scholar] [CrossRef] [Green Version]
  62. Kishi, S.; Cheng, C.; French, D.; Pei, D.; Das, S.; Cook, E.; Hijiya, N.; Rizzari, C.; Rosner, G.; Frudakis, T.; et al. Ancestry and pharmacogenetics of antileukemic dug toxicity. Blood 2007, 109, 4151–4157. [Google Scholar] [CrossRef] [PubMed]
  63. Mizzi, C.; Dalabira, E.; Kumuthini, J.; Dzimiri, N.; Balogh, I.; Basak, A.N.; Böhm, R.; Borg, J.; Borgiani, P.; Božina, N.; et al. A European Spectrum of Pharmacogenomic Biomarkers: Implications for Clinical Pharmacogenomics. PLoS ONE 2016, 11, e0162866. [Google Scholar] [CrossRef] [PubMed]
  64. Wang, S.-M.; Sun, L.-L.; Zeng, W.-X.; Wu, W.-S.; Zhang, G.-L. Effects of a microRNA binding site polymorphism in SLC19A1 on methotrexate concentrations in Chinese children with acute lymphoblastic leukemia. Med. Oncol. (Northwood Lond. Engl.) 2014, 31, 62. [Google Scholar] [CrossRef] [PubMed]
  65. Nacional Cancer Institute. Cancer Therapy Evaluation Program: Common Toxicity Criteria Manual; National Cancer Institute: Bethesda, MD, USA, 1999; pp. 1–29. [Google Scholar]
  66. BFM International Trials. Available online: http://www.bfm-international.org/trials.php (accessed on 28 February 2020).
  67. Neidhart, M. DNA Methylation and Complex Human Disease, 1st ed.; Academic Press: Cambridge, MA, USA, 2016. [Google Scholar]
  68. Relton, C.; Wilding, C.; Pearce, M.; Laffling, A.J.; Jonas, P.A.; Lynch, S.A.; Tawn, E.; Burn, J. Gene-gene interaction in folate-related genes and risk of neural tube defects in a UK population. J. Med. Genet. 2004, 41, 256–260. [Google Scholar] [CrossRef] [PubMed]
  69. Zinck, J.; MacFarlane, A. Approaches for the Identification of Genetic Modifiers of Nutrient Dependent Phenotypes: Examples from Folate. Front. Nutr. 2014, 1, 8. [Google Scholar] [CrossRef] [Green Version]
  70. Amos, W.; Driscoll, E.; Hoffman, J. Candidate genes versus genome-wide associations: Which are better for detecting genetic susceptibility to infectious disease? Proc. Biol. Sci. R. Soc. 2010, 278, 1183–1188. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Folate metabolic pathway disrupted by methotrexate (MTX). MTX—methotrexate; pg-MTX—polyglutamated MTX; THF—etrahydrofolate; DHF—dihydrofolate; SAM—S-adenosyl methionine; DHFR—dihydrofolate reductase; TYMS—thymidylate synthetase; MTHFR—methylenetetrahydrofolate reductase; SLC19A1—Solute Carrier Family 19 Member 1; SLCO1B1—Solute Carrier Organic Anion Transporter Family Member 1B1.
Figure 1. Folate metabolic pathway disrupted by methotrexate (MTX). MTX—methotrexate; pg-MTX—polyglutamated MTX; THF—etrahydrofolate; DHF—dihydrofolate; SAM—S-adenosyl methionine; DHFR—dihydrofolate reductase; TYMS—thymidylate synthetase; MTHFR—methylenetetrahydrofolate reductase; SLC19A1—Solute Carrier Family 19 Member 1; SLCO1B1—Solute Carrier Organic Anion Transporter Family Member 1B1.
Genes 11 00468 g001
Figure 2. Number of MTX low clearance episodes in SLC19A1 c.80G>A genotype carriers.
Figure 2. Number of MTX low clearance episodes in SLC19A1 c.80G>A genotype carriers.
Genes 11 00468 g002
Table 1. Folate pathway genetic variant frequencies of Serbian acute lymphoblastic leukemia (ALL) patients, the Serbian control group and the European descent control group. For the European control group, data were extracted from 1000 genome project for variants in all genes of interest except thymidylate synthetase (TYMS). Frequencies of TYMS variants in European population was estimated from British, American, Italian and Portuguese Caucasians [25,26,27,28,29].
Table 1. Folate pathway genetic variant frequencies of Serbian acute lymphoblastic leukemia (ALL) patients, the Serbian control group and the European descent control group. For the European control group, data were extracted from 1000 genome project for variants in all genes of interest except thymidylate synthetase (TYMS). Frequencies of TYMS variants in European population was estimated from British, American, Italian and Portuguese Caucasians [25,26,27,28,29].
Genetic VariantsGenotypeSerbian Childhood ALL PatientsSerbian Control GroupEuropean Caucasian Populationp Value b
N%HWMAF (%)N%HWMAF (%)N%MAF (%)
TYMSa rs34743033
28bp VNTR
3R/3R4933.10.04745.73331.70.8142.898129.946.30.39
3R/2R6141.25149.0156647.7
2R/2R3725.01918.373522.4
3R/4R10.7NANA11.0NANANANANANA
TYMS rs34489327
6bp indel
Ins/Ins7047.30.5733.14745.20.9233.2153246.332.80.92
Ins/Del6141.24543.3138641.8
Del/Del1711.51211.539411.9
MTHFR rs1801133
c.677 C>T
C/C7047.30.5732.111347.50.01333.820440.636.50.32
C/T6141.28937.423145.9
T/T1711.53615.16813.5
MTHFR rs1801131
c.1298 A>C
A/A7148.60.25231.411548.30.76130.923947.531.30.86
A/C5739.09941.621342.3
C/C1812.32410.15110.1
DHFR rs442767
-680 C>A
C/C6344.7131.44543.30.8634.623246.131.90.45
C/A6344.74644.222143.9
A/A1510.61312.5509.9
DHFR rs1643641
-675 A>G
A/A7452.50.3127.45552.90.3126.028656.925.20.82
A/G5337.64442.318035.8
G/G149.954.8377.4
DHFR rs1650695
-556 T>C
T/T7351.80.5427.45552.90.5226.428556.725.30.74
T/C5539.04341.318136.0
C/C139.265.8377.4
DHFR rs1650696
-464 A>T
A/A7251.10.6827.75552.90.5226.428656.925.20.72
A/T5639.74341.318035.8
T/T139.265.8377.4
DHFR rs408626
-317 A>G
A/A4934.80.6028.23735.60.9940.415029.845.00.22
A/G7150.45048.125350.3
G/G2114.91716.310019.9
SLC19A1 rs1051266
c.80 G>A
G/G3926.90.4049.07230.30.8065.715831.445.10.92
G/A6746.211648.723646.9
A/A3926.95021.010921.7
SLCO1B1 rs4149056
c.521 T>C
T/T10071.90.0615.510074.6113.435169.816.10.28
T/C3223.03223.914228.2
C/C75.021.5102.0
a Rare 4R allele of TYMS variant rs34743033 was disregarded in all analyses. b Probability obtained after chi square testing for differences in allele frequencies between Serbian and European control groups; HW—Hardy–Weinberg statistics. The number represents p value for HW equilibrium testing; NA—not available/not applicable; MAF—minor allele frequency. R—repeats. N—number of subjects
Table 2. Number of ALL patients experiencing different grade of toxicity during consolidation phase.
Table 2. Number of ALL patients experiencing different grade of toxicity during consolidation phase.
mM Protocol (MD-MTX), n = 102
Toxicity Grade
M Protocol (HD-MTX), n = 28
Toxicity Grade
0123401234
GIT toxicity6117222-15391-
Oral mucositis751962-1945--
Liver toxicity811641-214111
Nephrotoxicity993---2341--
Skin toxicity102----26-2--
Neurotoxicity9633--2611--
GIT—gastrointestinal tract.
Table 3. Genetic variation in folate pathway genes as predictive markers of GIT toxicity, hepatotoxicity and oral mucositis (Toxicity grade 0-1 versus 2-4). To assess the association of genetic variants with medium to severe toxicity, logistic regression was used, and all probabilities were adjusted for age, gender and MTX dose. The recessive genic model was applied for variants in TYMS, methylenetetrahydrofolate reductase (MTHFR), SLC19A1 and SLCO1B1 genes, with the reference group comprising common homozygotes and heterozygotes. Patients with denoted dihydrofolate reductase (DHFR) haplotypes were compared against the patients who carried all other haplotypes, the latter group being the reference group for statistical analyses.
Table 3. Genetic variation in folate pathway genes as predictive markers of GIT toxicity, hepatotoxicity and oral mucositis (Toxicity grade 0-1 versus 2-4). To assess the association of genetic variants with medium to severe toxicity, logistic regression was used, and all probabilities were adjusted for age, gender and MTX dose. The recessive genic model was applied for variants in TYMS, methylenetetrahydrofolate reductase (MTHFR), SLC19A1 and SLCO1B1 genes, with the reference group comprising common homozygotes and heterozygotes. Patients with denoted dihydrofolate reductase (DHFR) haplotypes were compared against the patients who carried all other haplotypes, the latter group being the reference group for statistical analyses.
GeneVariantdbSNPAllelesMAFOral MucositisGIT toxicityHepatotoxicity
OR95% CIpOR95% CIpOR95% CIp
MTHFRc.677 C>Trs1801133C/T0.321.540.29–8.070.610.970.28–3.410.97NANA1.00
MTHFRc.1298 A>Crs1801131A/C0.321.170.21–6.520.860.360.07–1.730.203.370.55–20.700.19
TYMS28bp repeatsrs347430332R–4R2R: 0.46 4R: 0.0030.430.09–2.130.300.750.29–1.920.551.570.34–7.320.56
TYMS6bp indel rs34489327Ins/Del0.322.770.60–12.730.194.171.25–13.910.020NANA1.00
DHFR hap1promoter region (−680–−317)rs442767–rs1643641– rs1650695–rs1650696– rs408626A-A-T-A-G0.330.370.09–1.600.181.020.44–2.370.965.830.65–52.610.12
DHFR hap2C-A-T-A-A0.308.070.97–67.50.0541.420.62–3.260.411.180.24–5.650.84
DHFR hap3C-A-T-A-G0.070.730.08–6.880.780.650.16–2.610.553.080.50–18.900.22
DHFR hap4C-G-C-T-A0.271.150.30–4.410.841.190.52–2.730.680.490.09–2.700.41
SLC19A1c.80 G>Ara1051266G/A0.501.090.26–4.570.901.380.56–3.400.499.701.70–55.780.011
SLCO1B1c.521 T>Crs4149056T/C0.16NANA1.001.650.27–10.020.59NANA1.00
GIT—gastrointestinal tract; MAF—minor allele frequency; OR—odds ratio; CI—confidence interval; NA—not applicable because OR could not be estimated.

Share and Cite

MDPI and ACS Style

Kotur, N.; Lazic, J.; Ristivojevic, B.; Stankovic, B.; Gasic, V.; Dokmanovic, L.; Krstovski, N.; Milosevic, G.; Janic, D.; Zukic, B.; et al. Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment. Genes 2020, 11, 468. https://doi.org/10.3390/genes11040468

AMA Style

Kotur N, Lazic J, Ristivojevic B, Stankovic B, Gasic V, Dokmanovic L, Krstovski N, Milosevic G, Janic D, Zukic B, et al. Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment. Genes. 2020; 11(4):468. https://doi.org/10.3390/genes11040468

Chicago/Turabian Style

Kotur, Nikola, Jelena Lazic, Bojan Ristivojevic, Biljana Stankovic, Vladimir Gasic, Lidija Dokmanovic, Nada Krstovski, Goran Milosevic, Dragana Janic, Branka Zukic, and et al. 2020. "Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment" Genes 11, no. 4: 468. https://doi.org/10.3390/genes11040468

APA Style

Kotur, N., Lazic, J., Ristivojevic, B., Stankovic, B., Gasic, V., Dokmanovic, L., Krstovski, N., Milosevic, G., Janic, D., Zukic, B., & Pavlovic, S. (2020). Pharmacogenomic Markers of Methotrexate Response in the Consolidation Phase of Pediatric Acute Lymphoblastic Leukemia Treatment. Genes, 11(4), 468. https://doi.org/10.3390/genes11040468

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop