- Open Access
Molecular insight into thiopurine resistance: transcriptomic signature in lymphoblastoid cell lines
Genome Medicinevolume 7, Article number: 37 (2015)
There has been considerable progress in the management of acute lymphoblastic leukemia (ALL) but further improvement is needed to increase long-term survival. The thiopurine agent 6-mercaptopurine (6-MP) used for ALL maintenance therapy has a key influence on clinical outcomes and relapse prevention. Genetic inheritance in thiopurine metabolism plays a major role in interindividual clinical response variability to thiopurines; however, most cases of thiopurine resistance remain unexplained.
We used lymphoblastoid cell lines (LCLs) from healthy donors, selected for their extreme thiopurine susceptibility. Thiopurine metabolism was characterized by the determination of TPMT and HPRT activity. We performed genome-wide expression profiling in resistant and sensitive cell lines with the goal of elucidating the mechanisms of thiopurine resistance.
We determined a higher TPMT activity (+44%; P = 0.024) in resistant compared to sensitive cell lines, although there was no difference in HPRT activity. We identified a 32-gene transcriptomic signature that predicts thiopurine resistance. This signature includes the GTPBP4 gene coding for a GTP-binding protein that interacts with p53. A comprehensive pathway analysis of the genes differentially expressed between resistant and sensitive cell lines indicated a role for cell cycle and DNA mismatch repair system in thiopurine resistance. It also revealed overexpression of the ATM/p53/p21 pathway, which is activated in response to DNA damage and induces cell cycle arrest in thiopurine resistant LCLs. Furthermore, overexpression of the p53 target gene TNFRSF10D or the negative cell cycle regulator CCNG2 induces cell cycle arrest and may also contribute to thiopurine resistance. ARHGDIA under-expression in resistant cell lines may constitute a novel molecular mechanism contributing to thiopurine resistance based on Rac1 inhibition induced apoptosis and in relation with thiopurine pharmacodynamics.
Our study provides new insights into the molecular mechanisms underlying thiopurine resistance and suggests a potential research focus for developing tailored medicine.
Approximately 6,000 patients are diagnosed with acute lymphoblastic leukemia (ALL) each year in the USA, including about two-thirds who are younger than 20 years, making ALL the most common malignancy in children and adolescents . Clinical outcomes of childhood ALL have improved considerably over time, and the overall 5-year event-free survival rate now reaches 80%, and even 90% when treatment strategies are selected based on the biological features of the leukemic cells and the pharmacodynamic and pharmacogenomic characteristics of the patient [1-3]. Nevertheless, this leaves about 20% of patients who experience relapses with far lower survival rates that decrease with each relapse . Furthermore, the outcomes are less favorable in adults, whose complete recovery rates rarely exceed 40% and who exhibit greater resistance to, and poorer tolerance of, therapeutic agents compared to children .
The treatment of ALL relies on combination chemotherapy. One of the cornerstone drugs for both intensification and maintenance therapy is the thiopurine agent 6-mercaptopurine (6-MP) . Effective maintenance therapy is essential for stabilizing the remission by suppressing the re-emergence of drug-resistant clones via a continuous reduction in the burden of residual leukemic cells. Recent studies have shown that, even when multiple drugs are used, the response to single drugs exerts a major influence on the relapse risk and that 6-MP is among the drugs with the greatest influence on remission maintenance . Therefore, elucidating the molecular basis of 6-MP resistance is crucial to relapse prediction, which allows optimization of the treatment strategy .
Genetic inheritance in thiopurine metabolism plays a major role in the interindividual variability that characterizes the clinical response to thiopurine agents [7,8]. Extensive pharmacogenetic studies have focused on the enzymes involved in thiopurine metabolism, such as thiopurine S-methyltransferase (TPMT), hypoxanthine phosphoribosyltransferase (HPRT), and inosine triphosphate pyrophosphatase (ITPA). These studies have shed light on the toxicity of thiopurine agents but have provided little information on thiopurine resistance [7-10]. To date, except rare cases of HPRT-deficiency or ultra-high TPMT activity, which impair the production of active thiopurine metabolites, most cases of thiopurine resistance remain unexplained [11-16].
Transcriptomic analysis is a powerful tool for characterizing susceptibility and resistance to drugs [17,18]. This approach can uncover previously unrecognized mechanisms of drug response and provides information on the associated biological pathways. Lymphoblastoid cell lines (LCLs) constitute a well-established pharmacogenomic model for genome-wide expression profiling [18,19]. Although there has been some debate about biological noise related to confounding factors, LCLs have been used to assess gene sets involved in responses to anticancer drugs such as bleomycin, gemcitabine, cytosine arabinoside, and 5-fluorouracil [20-24]. A gene set analysis of a vast panel of LCLs from different ethnic groups was conducted to assess associations between basal gene expression and thiopurine susceptibility . The 3′,5′-cyclic-AMP phosphodiesterase activity and the γ-aminobutyric acid catabolic process were found to be involved in the thiopurine response . However, this study was designed to evaluate gene expression profiles associated with a broad range of thiopurine susceptibility levels rather than with thiopurine resistance.
Here, we investigated the molecular basis of thiopurine resistance by exploring the whole-genome basal transcriptomic profiles involved in 6-MP resistant and sensitive phenotypes of LCLs originating from unrelated healthy individuals and selected by in vitro growth inhibition assays. We then used these transcriptomic profiles to identify genes predicting thiopurine resistance and relevant 6-MP metabolic pathways.
We screened 53 LCLs originating from consenting unrelated healthy adults and obtained via a collaboration program with the National Laboratory for the Genetics of Israeli Populations (NLGIP), Tel-Aviv University, Israel [25,26]. Six additional LCLs from male patients with the recessive genetic disease Lesch-Nyhan syndrome were obtained from the biobank of the Biochemistry and Molecular Biology Department of the Cochin University Hospital, Paris, France. They served as positive phenotypic controls for thiopurine resistance, as Lesch-Nyhan syndrome is characterized by HPRT deficiency. Cells were cultured as described elsewhere .
Ethical conduct of research
The authors state that they have obtained appropriate institutional review board approval for the collection of these samples in accordance with local legislation, or have followed the ethical principles outlined in the Declaration of Helsinki for experimental investigations. In addition, informed written consent has been obtained from the healthy donors and patients, or their guardians, involved in this study.
Roswell Park Memorial Institute medium, L-glutamine, and antibiotics for cell culture were purchased from Life Technologies (Carlsbad, CA, USA), fetal bovine serum from GE Healthcare (Little Chalfont, UK), and Falcon cell culture materials from Fisher Scientific (Waltham, MA, USA). Drugs (6-MP, azathioprine and 6-thioguanine [6-TG]) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Stock solutions (concentration, 5 mM) were prepared in 0.1 N sodium hydroxide and diluted in phosphate-buffered saline (PBS) to working solutions containing at least 25-fold the final tested concentrations.
Cell proliferation assay
Drug concentrations for assessing growth inhibition were 2 μM for 6-MP, 5 μM for azathioprine, and 0.5 μM for 6-TG. These concentrations were close to the mean half-maximal inhibitory concentration (IC50) and allowed optimal LCL classification . The cells were diluted to 200,000 /mL then incubated in Falcon 96-well plates (Waltham, MA, USA) in a volume of 200 μL (40,000 cells/well) for 3 days, with drug working solutions added as needed in three replicates, and six replicates for the controls (cells treated with 20 μL PBS). After 72 h, the tetrazolium derivative MTS reagent (CellTiter 96® AQueous One Solution Cell Proliferation Assay, Promega, Madison, WI, USA) was added (volume, 40 μL) to each well, including blank wells containing only PBS solution. After further incubation for 4 h, absorption at 490 nm was measured using a microplate reader spectrophotometer (Safir™ Tecan, Männedorf, Switzerland), which is ascribable to the living cells present in the medium. Growth inhibition relative to control was assessed for each cell line at least twice on two different batches from different cell vials thawed from the liquid nitrogen stock. Reproducibility of drug susceptibility for repeated thawing cycles was high, as previously reported .
Basal cell growth rate
Basal cell growth rates were estimated for each cell line during phenotyping experiments. After 72 h of the proliferation assay, cell density (N t ) was estimated in PBS control wells using MTS reagent, as described above, and calibration curve. The basal cell growth rate (r) was calculated using the following formula appropriate for the usual exponential kinetics of cell growth after a defined time (t): N t = N 0.2tr.
Nucleic acid extraction and quantification
Nucleic acids were extracted from cells incubated under optimal growth conditions with no added drugs. DNA and total RNA were extracted from cell pellets using QiAmp® DNA miniKit and miRNeasy® miniKit (Qiagen, Venlo, The Netherlands), respectively, according to the manufacturer’s instructions. DNA was quantified using an ND-1000 spectrophotometer (Nanodrop technologies, Wilmington, DE, USA). RNA quality and quantity were assessed using the 2100-Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
After validation of RNA quality (RIN score ≥8), 50 ng of total RNA was reverse-transcribed using the Ovation PicoSL WTA System V2 (NuGEN Technologies, West Cumbria, UK), according to the manufacturer’s instructions. Biotin-labelled cDNA was then hybridized to GeneChip® Human Gene 2.0ST microarrays (Affymetrix, Santa Clara, CA, USA) at 45°C for 17 h. The microarrays were washed on the fluidic station FS450 according to the specific manufacturer’s protocols and scanned using the GCS3000 7G (Affymetrix). The scans were then analyzed with Expression Console software (Affymetrix) to obtain raw data (.cel files) and metrics for quality controls. Examination of these quality-control metrics and of raw-data distribution showed no outlier samples.
Data were normalized using the Robust Multi-array Average (RMA) algorithm in R software with the custom chip description file (CDF) version 17.0.0 . Data are available on the NCBI Gene Expression Omnibus (GEO) via the accession number GSE61905 .
Differentially expressed gene enrichment analysis was carried out using the DAVID bioinformatics resources (NIH), based on gene ontology (GO) biological processes, and Ingenuity Pathways Analysis (Ingenuity® Systems, USA) [29-31]. To identify a transcriptomic signature predicting thiopurine resistance, we applied the ‘nearest shrunken centroids’ method using the Prediction Analysis of Microarrays (PAM) R package, which identifies predictive classifier genes .
Microfluidic-based RT-qPCR assay
Microfluidic-based quantitative PCR assay was performed to validate the differential microarray expression patterns of the molecular signature genes. Transcripts were quantified for 40 genes of interest and four reference genes (GUSB/GAPDH/RPL13A/B2M). The RT-qPCR assay is detailed in the Additional file 1: Data S1 and Table S2.
Intracellular ATP assay
Intracellular ATP was assayed on the day the cell proliferation assay was performed. Cell samples were kept on ice until the assay, as ATP is unstable. The assay was performed using the ATP Bioluminescence Assay Kit HS II (Roche, Germany) according to the manufacturer’s instructions. Briefly, luminescence measured in the microplate wells was related to intracellular ATP after the addition of luciferase and D-luciferin. Intracellular ATP concentrations were calculated using a calibration curve.
EBV and mtDNA copy number
Copy number variations (CNVs) of Epstein-Barr virus (EBV) and mitochondrial (mt) DNA were determined in the LCL DNA samples using TaqMan® CNV Assays (Life Technologies, Carlsbad, CA, USA). Custom TaqMan® assays were designed using Primer 3; primer and probe sequences are reported in Additional file 1: Table S3. The EBV CNV assay interrogated a 66-bp fragment at the DNA polymerase locus multiplexed with an 87-bp fragment assay from RNAse P (RPPH1 gene) as an internal reference, using the VIC® dye-labeled TAMRA™ probe and sequence-specific forward and reverse primers (Applied assay ID 4403328). The mtDNA copy number assay examined a 72-bp fragment at the ND2 locus multiplexed with a 90-bp fragment assay from NRF1 gene as an internal reference. Final concentrations for EBV primers, EBV probe, mtDNA primers, mtDNA probe, NRF1 primers and NRF1 probe were 90 nM, 250 nM, 30 nM, 250 nM, 900 nM, and 250 nM, respectively. RT-qPCR TaqMan® reactions were performed in 384-well plates (reaction volume, 10 μL) using 10 ng of DNA template (concentration, 5 ng/μL) and 5 μL of TaqMan® Genotyping Master Mix, according to the manufacturer’s protocol. Equal efficiency of amplification was observed for each assay in the multiplex reaction.
Gene CNVs were determined using the previously described 2-ΔΔCt method . The absolute EBV copy number was related to a calibrator DNA template from the Burkitt lymphoma-derived Namalwa cell line (ECACC, UK), which was determined by fluorescence in-situ hybridization to have integrated EBV copies in the diploid genome . The relative mtDNA copy number was related to the sample having the highest Ct (that is, the lowest copy number).
Intracellular enzymatic activities
TPMT activity was assessed using the previously described reverse-phase high-performance liquid chromatography (HPLC) method therefore adapted for lymphocyte pellets . The pellets were lysed by two freezing/thaw cycles in 200 μL PBS. Assay results are reported in pmol of 6-methylmercaptopurine formed per hour and per mg of total protein.
HPRT activity was measured as the rate of inosine monophosphate (IMP) production, using the Precice® kit (Novocib, Lyon, France). Briefly, as described elsewhere, the assay is based on the effect of IMP-dehydrogenase, which catalyzes the irreversible oxidation of IMP to xanthosine monophosphate while simultaneously reducing NAD to NADH2, whose production is monitored directly at 340 nm using a microplate reader spectrophotometer . Assay results are reported in nmol of IMP formed per hour and per mg of total protein.
Data are described as mean ± standard error (SEM). To compare 6-MP resistant and sensitive cell lines, we used the non-parametric Mann–Whitney t test. Correlations were assessed using the Spearman r s test.
Gene expression values were analyzed using unsupervised hierarchical clustering and principal component analysis (PCA). Then, to identify differentially expressed genes, we performed the parametric Student t test and computed fold-changes for the expression of each gene in thiopurine resistant versus sensitive cell lines. Genes associated with P values <0.01 were selected for functional bioinformatics analyses using DAVID and Ingenuity®. Statistical analyses were carried out using Partek® Genomics Suite™ (Partek Inc., St. Louis, MO, USA) and Prism 5.0 (GraphPad, San Diego, CA, USA).
LCL selection and characterization
The relative susceptibility to growth inhibition by thiopurine drugs was determined for 53 LCLs from healthy adults. Growth inhibition by 6-MP was closely correlated to growth inhibition by azathioprine and by 6-TG (r s = 0.95 and r s = 0.81, respectively; P <0.0001) (Additional file 2: Figure S1). Growth inhibition by 6-MP was about 18 times greater for the most sensitive than for the most resistant cell lines (10th to 90th percentile growth inhibition, 12.4% to 42.6%). We performed a genome-wide expression analysis of 11 cell lines with an extreme phenotype in terms of 6-MP susceptibility, selected among the 53 LCLs: five resistant and six sensitive, with mean growth inhibition by 6-MP of 11.8 ± 2.1% and 39.8 ± 2.3%, respectively (P = 0.008) (Figure 1). The LCLs from the six patients with Lesch-Nyhan syndrome were almost completely resistant to thiopurines (mean growth inhibition by 6-MP, 1 ± 2%) (Figure 1). We verified that these cell lines carrying an inherited HPRT deficiency exhibited undetectable HPRT activity (mean: 0.2 ± 2 nmol/h/mg protein).
Resistant cell lines showed a trend toward a lower basal growth rate per day compared to sensitive cell lines (0.3 ± 0.06 vs. 0.5 ± 0.04, respectively; P = 0.052). Moreover, TPMT activity was about 44% higher in resistant compared to sensitive cell lines (425 ± 20 vs. 295 ± 24 pmol/h/mg protein; P = 0.024). None of the other study parameters (EBV copy number, mtDNA copy number, intracellular ATP level, and HPRT activity) differed significantly between 6-MP resistant and sensitive cell lines (Table 1).
PCA graphically discriminated between the 6-MP resistant and sensitive 11 cell lines. The first three components explained 37.6% of the total variance (Additional file 3: Figure S2).
Of the 23,786 genes analyzed in the micro-array, 943 and 210 had different basal expression levels in the two groups at P values of <0.01 and <0.001, respectively. Of the 210 genes significantly different at the P < 0.001 level, 84 were upregulated and 126 downregulated in the resistant compared to the sensitive cell lines (Additional file 1: Table S4).
For identifying the most relevant genes that were differentially expressed in 6-MP resistant compared to sensitive cell lines, we performed PAM. This shrunken centroid algorithm identified 40 genomic loci that sorted the cell lines between resistant and sensitive, with a cross-validation error rate of 0.09 (Additional file 1: Table S5). To validate the expression signature determined from the micro-array, we performed RT-qPCR. Among the 40 genomic loci, seven were related to unknown transcripts, non-coding RNA, or unknown proteins and were consequently not selected for RT-qPCR. Furthermore, we were not able to amplify one gene (CNR1). The remaining 32 genes were validated for the transcriptomic signature using GUSB as the reference gene (r s = 0.87; P <0.0001) (Figure 2). These results were reproducible using the three other reference genes, GAPDH, RPL13A, and B2M (data not shown).
Comprehensive pathway analysis
The analysis of the GO biological processes involved in differential gene expression (P <0.01) led to 119 terms each containing at least five genes. After a Benjamini correction adjusted at a significance level (0.01/119), we found 10 GO terms, including 122 single genes, involved in the phenotype difference (Table 2; Additional file 1: Table S6). Half these terms were involved in cell proliferation functions and the other half chiefly in RNA processes. Ingenuity® predicted that nine upstream transcriptional regulators were activated (TP53, CD24, NUPR1, CDKN1A) or inhibited (FOXM1, FLI1, MYC, CSF2, CCDN1) in resistant cell lines (P <0.01). These regulators targeted 128 of the differentially expressed genes (Table 3). Furthermore 18 canonical pathways played a significant biological role in thiopurine resistance (P <0.01) (Figure 3, Additional file 1: Table S7 and Table S8).
Using a pharmacogenomic LCL-based model, we performed a comprehensive analysis of molecular resistance to thiopurines. To our knowledge, we identified for the first time a 32-gene transcriptomic signature predicting thiopurine resistance. Of the 32 genes, 22 were upregulated and 10 downregulated. Our transcriptomic analysis of untreated and phenotypically selected cell lines identified potential biomarkers for thiopurine resistance and suggested metabolic pathways that might constitute therapeutic targets for overcoming thiopurine resistance.
Our model based on LCL growth inhibition is of interest for studying thiopurine resistance, as thiopurines are used to target lymphoblasts in patients with ALL. Moreover, in autoimmune diseases, in which thiopurine therapy also plays a major role, lymphocytes are the target cells and lymphocyte apoptosis must be achieved to induce immunosuppression . Growth inhibition by 6-MP was closely correlated to growth inhibition by azathioprine and by 6-TG. A proof-of-concept study based on a panel of LCLs showed high goodness-of-fit values for linear regression plots comparing growth inhibition profiles of paired drugs acting via a shared pathway . Thus, cell lines resistant to 6-MP were also resistant to azathioprine and to 6-TG, indicating that our results on 6-MP resistance are likely to be representative of thiopurine drugs. In addition, LCLs from patients with Lesch-Nyhan syndrome, which lacked HPRT activity because of a recessive X-linked gene defect, were completely resistant to thiopurines. This finding reflects the inability of HPRT-deficient cells to bioactivate thiopurines into active cytotoxic metabolites, mainly 6-thioguanine nucleotides [11,13]. The resistance phenotype of HPRT-deficient LCLs constitutes an appropriate positive control for our in vitro model. Altogether, these findings support the validity of our LCL-based pharmacogenomics model for studying thiopurine resistance.
Targeted analysis of thiopurine metabolizing enzymes
The antiproliferative effects of the purine analogue 6-MP are due to metabolites that have three mechanisms of action: inhibition of de novo purine synthesis; cell cycle arrest due to metabolite incorporation into DNA; and, particularly in lymphocytes, increased apoptosis due to Rac1 inhibition . We first studied variations in the main thiopurine-metabolizing enzymes known to influence the pharmacological response to thiopurine drugs. We did not study xanthine oxidase, which is not present in lymphocytes. HPRT activity was not significantly different between resistant and sensitive cell lines. Variations in HPRT activity have been reported to be associated with thiopurine resistance and hematological toxicity [13,38]. However, HPRT activity does not vary widely within the general population and probably has little influence on clinical outcomes . By contrast, the considerable interindividual variations in TPMT activity in the general population, which are related to genetic polymorphisms, affect both the toxicity and the efficacy of thiopurine agents . Higher ALL remission rates have been reported in patients with a decreased TPMT activity . Interestingly, TPMT activity was slightly higher in the resistant than in the sensitive cell lines in our study. This finding confirms the association between drug metabolism and thiopurine susceptibility found in our LCLs, in keeping with data from treated patients . However, TPMT expression levels were not significantly different between 6-MP resistant and sensitive LCLs, probably because TPMT regulation is mainly post-translational [40,41].
Comprehensive pathway analysis
Transcriptomic analysis identified 210 genes that were significantly upregulated or downregulated in resistant cell lines. Our GO analysis of these genes showed that 10 GO terms were enriched in these genes, including five related with the cell cycle, most notably the M phase. Ingenuity® pathway analysis predicted CDKN1A activation in resistant cell lines. CDKN1A encodes a potent cyclin-dependent kinase (CDK) inhibitor, also called p21WAF1/CIP1, which binds to and inhibits the cyclin-CDK2 or -CDK4 complexes, preventing the phosphorylation of critical CDK substrates and blocking cell cycle progression . Thus, CDKN1A acts as a negative regulator of cell cycle progression at G1. More specifically, when located in the nucleus, p21WAF1/CIP1 controls the cell cycle and DNA replication, whereas cytoplasmic p21WAF1/CIP1 has been implicated in apoptosis inhibition . A study of human cancer cells showed that increased p21WAF1/CIP1 levels, related to phosphatidylinositol 3-kinase (PI3K) pathway inhibition, induced chemoresistance by causing a cell cycle delay . Moreover, resistance to another anticancer drug, taxol, has been reported in breast-cancer cells exhibiting upregulation of p21WAF1/CIP1 . A recent study also identified p21WAF1/CIP1 expression as a major factor in resistance to promising anticancer drugs acting within the cell cycle . The tumor-suppressor protein p53 tightly controls p21WAF1/CIP1, through which it mediates the p53-dependent cell cycle G1-phase arrest in response to a variety of stress stimuli. We found upregulation of p53 target genes and of ATM, a p53 upstream regulator, which is activated in response to DNA damage, to be related with thiopurine resistance. Moreover, p21WAF1/CIP1 mediates NUPR1-induced chemoresistance, and our analysis predicted NUPR1 activation in resistant cell lines .
Another p53 target gene is TNFRSF10D, which encodes tumor necrosis factor-related apoptosis-inducing ligand receptor 4 (TRAILR4) and whose overexpression was a component of the transcriptomic signature identified in our study. Increased TNFRSF10D expression has been found to be associated with chemoresistance . Taken in concert, these results suggest upregulation of the ATM/p53/p21 DNA damage response pathway in resistant cell lines, with resulting inhibition of the cyclin-CDK2 or -CDK4 complexes and cell cycle arrest (Additional file 4: Figure S3).
Furthermore, a study of the genome-scale protein-interaction profile of p53 showed that GTPBP4 was a p53 interactor involved in 60S ribosome biogenesis . This nucleolar GTP-binding protein, whose downregulation was a component of our transcriptomic signature in resistant cell lines, has been reported to activate p53 when silenced . Finally, in keeping with p53 activation, MYC inhibition, associated with cell cycle repression, was predicted in resistant cell lines .
Our canonical pathway analysis identified 18 significant processes that were differentially expressed between 6-MP resistant and sensitive cell lines, among which at least five were related to DNA repair in response to damage, including ‘Role of CHK proteins in cell cycle checkpoint control’ and ‘cell cycle: G2/M DNA damage checkpoint regulation’. Many anticancer drugs acting as antimetabolites require involvement of the DNA mismatch repair (MMR) system to exert their cellular responses . The primary function of the MMR system is to edit and repair DNA replication errors and DNA damage . Loss of MMR has been observed in a variety of human cancers and is associated with resistance to several anticancer agents such as etoposide, cisplatin, carboplatin, and 5-fluorouracil . Human cancer cell lines lacking the MMR system were resistant to high doses of 6-TG compared to MMR-proficient cell lines . A study established that cell cycle arrest in G2-M after thiopurine treatment was mediated by single-strand breaks in MMR-proficient cells . In keeping with these findings, a recent study showed that high expression of PKCζ, a protein kinase believed to stabilize the MMR protein MSH2, increased the response to thiopurine therapy in pediatric patients with ALL . A characteristic of resistant cell lines demonstrated by our transcriptomic analysis was a significantly decreased expression of RFC3 and POLDIP2, both known to interact with the MMR system and cell replication. Thus, MMR deficiency and, possibly, the expression levels of RFC3 and POLDIP2, may help to predict thiopurine resistance. Furthermore, a study based on MOLT-4 cell lines also suspected the role of alterations in the MMR system in the resistance phenotype to 6-MP . Moreover, an increase in induced mutations after 6-TG treatment has been reported in MMR-deficient cell lines . In a study of 228 children with ALL previously treated with anticancer agents including thiopurines, CCNG2 under-expression was a risk factor for treatment-related myeloid leukemia (t-ML) . CCNG2, a negative regulator of cell cycle progression independent from p53, is induced in cell cycle arrest in response to DNA damage . Alterations in CCNG2 expression may enhance cell cycle progression and contribute to failure of the cells to respond to DNA-damage stimuli that would otherwise promote exit from the cell cycle; subsequently, the proliferation of cells carrying misrepaired DNA may lead to leukemic transformation . Conversely, our resistant cell lines exhibited a high level of CCNG2 expression that might stop the cell cycle at the G1/S phase, preventing the MMR system from promoting thiopurine susceptibility and thereby contributing to thiopurine resistance.
Finally, we found under-expression of ARHGDIA, encoding for Rho-GDP dissociation inhibitor alpha (RhoGDIα) in resistant cell lines. RhoGDIα regulates and sequesters in cytoplasm inactive GDP-bound forms of RhoGTPase, including Rac1, a molecular target inhibited by thiopurine nucleotides . Thus, RhoGDIα prevents RhoGTPase from being recruited at the cellular membrane where it can be activated. Decreased ARHGDIA expression may therefore increase the amount of potentially active Rac1, preventing effective Rac1 inhibition by the thiopurine nucleotide 6-thioguanosine triphosphate (6-TGTP). RhoGDIα thinly regulates RhoGTPase activation, involved in cellular processes and contributing to tumor invasion and metastasis [57,58]. Moreover, loss of RhoGDIα has been previously associated with tamoxifene resistance . This mechanism, related to thiopurine pharmacodynamics, may contribute to thiopurine resistance by reducing 6-TGTP induced apoptosis via Rac1 inhibition in lymphocytes . It may represent an original biomarker of thiopurine resistance.
The transcriptomic content of LCLs includes many genes from diverse cellular pathways and has proven valuable for studying genome-wide individual differences in alternative mRNA splicing . The reliability of the association between genomic analysis results and drug response phenotypes in LCL-based models deserves discussion. In vitro biological noise may limit the usefulness of LCLs as a pharmacogenomics research tool . Several parameters measured in our study, including EBV and mtDNA copy numbers, intracellular ATP level, and basal cell growth rate, have been described as potential confounding factors influencing the LCL drug response phenotype [20,25,61,62]. These non-genetic variables unrelated to the genomic status of the cell line can alter growth inhibition. They are of major concern when using LCLs produced from different EBV strains and generated by different laboratories . In our study, however, all the LCLs came from a single biobank (NLGIP) and were generated by a single lab using the same stock of B-95 EBV-expressing marmoset cell line . Furthermore, our resistant and sensitive cell lines exhibited no significant differences regarding the EBV and mtDNA copy numbers or intracellular ATP levels measured during growth inhibition experiments. These parameters are thus unlikely to have affected our transcriptomic analysis. However, we observed a trend toward a lower basal cell growth rate in thiopurine-resistant compared to thiopurine-sensitive cell lines. Thus, thiopurine susceptibility may be, at least in part, associated with the basal cell growth rate. This finding has been previously observed in a mechanistic mathematical modelling of 6-MP resistance . As discussed above, some of the 210 genes identified in our study may be related to the cell growth rate and cell cycle. When Elion and Hitchings synthesized 6-MP as a drug for treating ALL, their goal was to selectively inhibit nucleic acid synthesis in rapidly dividing cells such as leukemic cells . Thus, our study unveils molecular pathways associated with the mechanisms of action of thiopurines such as cell cycle arrest upon incorporation of thioguanine nucleotides.
In conclusion, our study using an LCL-based model identified a transcriptomic signature of thiopurine resistance. We used a well-established pharmacogenomics approach involving transcriptomic profiling of basal mRNA in a cell-line model, taking advantage of the extra power afforded by analyzing extreme phenotype cell lines [64,65]. Thus, genome-wide transcriptomic analysis of LCLs coupled with drug susceptibility phenotyping can identify novel candidate genes and pathways that may help to explain individual response to thiopurine drugs (Figure 4). Our study provides new insights into the molecular mechanisms underlying thiopurine resistance suggesting potential research focus for developing tailored medicine.
Acute lymphoblastic leukemia
Database for annotation visualization and integrated discovery
Gene expression omnibus
Lymphoblastoid cell lines
National center for biotechnology information
National laboratory for the genetics of Israeli populations
Prediction analysis of microarrays
Principal component analysis
Robust multi-array average
Reverse transcriptase-quantitative polymerase chain reaction
Inaba H, Greaves M, Mullighan CG. Acute lymphoblastic leukaemia. Lancet. 2013;381:1943–55.
Pui C-H, Evans WE. Treatment of acute lymphoblastic leukemia. N Engl J Med. 2006;354:166–78.
Izraeli S, Shochat C, Tal N, Geron I. Towards precision medicine in childhood leukemia–insights from mutationally activated cytokine receptor pathways in acute lymphoblastic leukemia. Cancer Lett. 2014;352:15–20.
Ko RH, Ji L, Barnette P, Bostrom B, Hutchinson R, Raetz E, et al. Outcome of patients treated for relapsed or refractory acute lymphoblastic leukemia: a Therapeutic Advances in Childhood Leukemia Consortium study. J Clin Oncol Off J Am Soc Clin Oncol. 2010;28:648–54.
Beesley AH, Firth MJ, Anderson D, Samuels AL, Ford J, Kees UR. Drug-gene modeling in pediatric T-cell acute lymphoblastic leukemia highlights importance of 6-mercaptopurine for outcome. Cancer Res. 2013;73:2749–59.
Fotoohi AK, Coulthard SA, Albertioni F. Thiopurines: factors influencing toxicity and response. Biochem Pharmacol. 2010;79:1211–20.
Lennard L, Lilleyman JS, Van Loon J, Weinshilboum RM. Genetic variation in response to 6-mercaptopurine for childhood acute lymphoblastic leukaemia. Lancet. 1990;336:225–9.
Chouchana L, Narjoz C, Beaune P, Loriot M-A, Roblin X. Review article: the benefits of pharmacogenetics for improving thiopurine therapy in inflammatory bowel disease. Aliment Pharmacol Ther. 2012;35:15–36.
Gurwitz D, Rodríguez-Antona C, Payne K, Newman W, Gisbert JP, de Mesa EG, et al. Improving pharmacovigilance in Europe: TPMT genotyping and phenotyping in the UK and Spain. Eur J Hum Genet EJHG. 2009;17:991–8.
Farfan MJ, Salas C, Canales C, Silva F, Villarroel M, Kopp K, et al. Prevalence of TPMT and ITPA gene polymorphisms and effect on mercaptopurine dosage in Chilean children with acute lymphoblastic leukemia. BMC Cancer. 2014;14:299.
Rosman M, Lee MH, Creasey WA, Sartorelli AC. Mechanisms of resistance to 6-thiopurines in human leukemia. Cancer Res. 1974;34:1952–6.
Fotoohi AK, Lindqvist M, Peterson C, Albertioni F. Impaired transport as a mechanism of resistance to thiopurines in human T-lymphoblastic leukemia cells. Nucleosides Nucleotides Nucleic Acids. 2006;25:1039–44.
Pieters R, Huismans DR, Loonen AH, Peters GJ, Hählen K, van der Does-van den Berg A, et al. Hypoxanthine-guanine phosphoribosyl-transferase in childhood leukemia: relation with immunophenotype, in vitro drug resistance and clinical prognosis. Int J Cancer J Int. Cancer. 1992;51:213–7.
Van Egmond R, Chin P, Zhang M, Sies CW, Barclay ML. High TPMT enzyme activity does not explain drug resistance due to preferential 6-methylmercaptopurine production in patients on thiopurine treatment. Aliment Pharmacol Ther. 2012;35:1181–9.
Chouchana L, Roche D, Jian R, Beaune P, Loriot M-A. Poor response to thiopurine in inflammatory bowel disease: how to overcome therapeutic resistance? Clin Chem. 2013;59:1023–6.
Chouchana L, Narjoz C, Roche D, Golmard J-L, Pineau B, Chatellier G, et al. Interindividual variability in TPMT enzyme activity: 10 years of experience with thiopurine pharmacogenetics and therapeutic drug monitoring. Pharmacogenomics. 2014;15:745–57.
Morag A, Pasmanik-Chor M, Oron-Karni V, Rehavi M, Stingl JC, Gurwitz D. Genome-wide expression profiling of human lymphoblastoid cell lines identifies CHL1 as a putative SSRI antidepressant response biomarker. Pharmacogenomics. 2011;12:171–84.
Wheeler HE, Dolan ME. Lymphoblastoid cell lines in pharmacogenomic discovery and clinical translation. Pharmacogenomics. 2012;13:55–70.
Stark AL, Dolan ME. Lymphoblastoid cell lines in pharmacogenomics: how applicable are they to clinical outcomes? Pharmacogenomics. 2013;14:447–50.
Stark AL, Zhang W, Mi S, Duan S, O’Donnell PH, Huang RS, et al. Heritable and non-genetic factors as variables of pharmacologic phenotypes in lymphoblastoid cell lines. Pharmacogenomics J. 2010;10:505–12.
Cloos J, de Boer WPH, Snel MHJ, van den Ijssel P, Ylstra B, Leemans CR, et al. Microarray analysis of bleomycin-exposed lymphoblastoid cells for identifying cancer susceptibility genes. Mol Cancer Res MCR. 2006;4:71–7.
Fridley BL, Batzler A, Li L, Li F, Matimba A, Jenkins GD, et al. Gene set analysis of purine and pyrimidine antimetabolites cancer therapies. Pharmacogenet Genomics. 2011;21:701–12.
Wu T-Y, Fridley BL, Jenkins GD, Batzler A, Wang L, Weinshilboum RM. Mycophenolic acid response biomarkers: a cell line model system-based genome-wide screen. Int Immunopharmacol. 2011;11:1057–64.
Li L, Fridley B, Kalari K, Jenkins G, Batzler A, Safgren S, et al. Gemcitabine and cytosine arabinoside cytotoxicity: association with lymphoblastoid cell expression. Cancer Res. 2008;68:7050–8.
Morag A, Kirchheiner J, Rehavi M, Gurwitz D. Human lymphoblastoid cell line panels: novel tools for assessing shared drug pathways. Pharmacogenomics. 2010;11:327–40.
Gurwitz D, Kimchi O, Bonne-Tamir B. The Israeli DNA and Cell Line Collection: A Human Diversity Repository. In: Populations and genetics: legal and socio-ethical perspectives. Leiden: Martinus Nijhoff; 2003. p. 95–113.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.
Database for Annotation, Visualization and Integrated Discovery (DAVID). [http://david.abcc.ncifcrf.gov/].
Ingenuity® pathway analysis. [http://www.ingenuity.com/products/ipa].
Tibshirani R, Hastie T, Narasimhan B, Chu G. Diagnosis of multiple cancer types by shrunken centroids of gene expression. Proc Natl Acad Sci U S A. 2002;99:6567–72.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods San Diego Calif. 2001;25:402–8.
Lawrence JB, Villnave CA, Singer RH. Sensitive, high-resolution chromatin and chromosome mapping in situ: presence and orientation of two closely integrated copies of EBV in a lymphoma line. Cell. 1988;52:51–61.
Anglicheau D, Sanquer S, Loriot MA, Beaune P, Thervet E. Thiopurine methyltransferase activity: new conditions for reversed-phase high-performance liquid chromatographic assay without extraction and genotypic-phenotypic correlation. J Chromatogr B Analyt Technol Biomed Life Sci. 2002;773:119–27.
Balakireva LA, Godard NH. A continuous spectrophotometric assay for rapid measurement of hypoxanthine-guanine phosphoribosyltransferase (HPRT) activity in cell lysates. Nat Methods Appl Notes. 2011.
Tiede I, Fritz G, Strand S, Poppe D, Dvorsky R, Strand D, et al. CD28-dependent Rac1 activation is the molecular target of azathioprine in primary human CD4+ T lymphocytes. J Clin Invest. 2003;111:1133–45.
Ding L, Zhang F, Liu H, Gao X, Bi H, Wang X, et al. Hypoxanthine guanine phosphoribosyltransferase activity is related to 6-thioguanine nucleotide concentrations and thiopurine-induced leukopenia in the treatment of inflammatory bowel disease. Inflamm Bowel Dis. 2012;18:63–73.
Lennard L, Hale JP, Lilleyman JS. Red blood cell hypoxanthine phosphoribosyltransferase activity measured using 6-mercaptopurine as a substrate: a population study in children with acute lymphoblastic leukaemia. Br J Clin Pharmacol. 1993;36:277–84.
Tai HL, Krynetski EY, Schuetz EG, Yanishevski Y, Evans WE. Enhanced proteolysis of thiopurine S-methyltransferase (TPMT) encoded by mutant alleles in humans (TPMT*3A, TPMT*2): mechanisms for the genetic polymorphism of TPMT activity. Proc Natl Acad Sci U S A. 1997;94:6444–9.
Milek M, Smid A, Tamm R, Kuzelicki NK, Metspalu A, Mlinaric-Rascan I. Post-translational stabilization of thiopurine S-methyltransferase by S-adenosyl-L-methionine reveals regulation of TPMT*1 and *3C allozymes. Biochem Pharmacol. 2012;83:969–76.
Koster R, di Pietro A, Timmer-Bosscha H, Gibcus JH, van den Berg A, Suurmeijer AJ, et al. Cytoplasmic p21 expression levels determine cisplatin resistance in human testicular cancer. J Clin Invest. 2010;120:3594–605.
McDonald GT, Sullivan R, Paré GC, Graham CH. Inhibition of phosphatidylinositol 3-kinase promotes tumor cell resistance to chemotherapeutic agents via a mechanism involving delay in cell cycle progression. Exp Cell Res. 2010;316:3197–206.
Hawthorne VS, Huang W-C, Neal CL, Tseng L-M, Hung M-C, Yu D. ErbB2-mediated Src and signal transducer and activator of transcription 3 activation leads to transcriptional up-regulation of p21Cip1 and chemoresistance in breast cancer cells. Mol Cancer Res MCR. 2009;7:592–600.
Chen Y-P, Lin H-J, Chen J-S, Tsai M-Y, Hsieh H-P, Chang J-Y, et al. CDKN1A-mediated responsiveness of MLL-AF4-positive acute lymphoblastic leukemia to Aurora kinase-A inhibitors. Int J Cancer J Int Cancer. 2014;135:751–62.
Vincent AJ, Ren S, Harris LG, Devine DJ, Samant RS, Fodstad O, et al. Cytoplasmic translocation of p21 mediates NUPR1-induced chemoresistance: NUPR1 and p21 in chemoresistance. FEBS Lett. 2012;586:3429–34.
Liu X, Yue P, Khuri FR, Sun S-Y. Decoy receptor 2 (DcR2) is a p53 target gene and regulates chemosensitivity. Cancer Res. 2005;65:9169–75.
Lunardi A, Di Minin G, Provero P, Dal Ferro M, Carotti M, Del Sal G, et al. A genome-scale protein interaction profile of Drosophila p53 uncovers additional nodes of the human p53 network. Proc Natl Acad Sci U S A. 2010;107:6322–7.
Dang CV, O’Donnell KA, Zeller KI, Nguyen T, Osthus RC, Li F. The c-Myc target gene network. Semin Cancer Biol. 2006;16:253–64.
Fink D, Aebi S, Howell SB. The role of DNA mismatch repair in drug resistance. Clin Cancer Res Off J Am Assoc Cancer Res. 1998;4:1–6.
Glaab WE, Risinger JI, Umar A, Barrett JC, Kunkel TA, Tindall KR. Resistance to 6-thioguanine in mismatch repair-deficient human cancer cell lines correlates with an increase in induced mutations at the HPRT locus. Carcinogenesis. 1998;19:1931–7.
Yan T, Berry SE, Desai AB, Kinsella TJ. DNA mismatch repair (MMR) mediates 6-thioguanine genotoxicity by introducing single-strand breaks to signal a G2-M arrest in MMR-proficient RKO cells. Clin Cancer Res Off J Am Assoc Cancer Res. 2003;9:2327–34.
Hartsink-Segers SA, Beaudoin JJ, Luijendijk MWJ, Exalto C, Pieters R, Den Boer ML. PKCζ and PKMζ are overexpressed in TCF3-rearranged paediatric acute lymphoblastic leukaemia and are associated with increased thiopurine sensitivity. Leukemia. 2015;29:304–11.
Panetta JC, Evans WE, Cheok MH. Mechanistic mathematical modelling of mercaptopurine effects on cell cycle of human acute lymphoblastic leukaemia cells. Br J Cancer. 2006;94:93–100.
Bogni A, Cheng C, Liu W, Yang W, Pfeffer J, Mukatira S, et al. Genome-wide approach to identify risk factors for therapy-related myeloid leukemia. Leukemia. 2006;20:239–46.
Bennin DA, Don ASA, Brake T, McKenzie JL, Rosenbaum H, Ortiz L, et al. Cyclin G2 associates with protein phosphatase 2A catalytic and regulatory B’ subunits in active complexes and induces nuclear aberrations and a G1/S phase cell cycle arrest. J Biol Chem. 2002;277:27449–67.
Etienne-Manneville S, Hall A. Rho GTPases in cell biology. Nature. 2002;420:629–35.
Zhu Y, Tummala R, Liu C, Nadiminty N, Lou W, Evans CP, et al. RhoGDIα suppresses growth and survival of prostate cancer cells. Prostate. 2012;72:392–8.
Barone I, Brusco L, Gu G, Selever J, Beyer A, Covington KR, et al. Loss of Rho GDIα and resistance to tamoxifen via effects on estrogen receptor α. J Natl Cancer Inst. 2011;103:538–52.
Coulombe-Huntington J, Lam KCL, Dias C, Majewski J. Fine-scale variation and genetic determinants of alternative splicing across individuals. PLoS Genet. 2009;5:e1000766.
Choy E, Yelensky R, Bonakdar S, Plenge RM, Saxena R, De Jager PL, et al. Genetic analysis of human traits in vitro: drug response and gene expression in lymphoblastoid cell lines. PLoS Genet. 2008;4:e1000287.
Caliskan M, Cusanovich DA, Ober C, Gilad Y. The effects of EBV transformation on gene expression levels and methylation profiles. Hum Mol Genet. 2011;20:1643–52.
Elion GB. The purine path to chemotherapy. Science. 1989;244:41–7.
Gurwitz D. Expression profiling: a cost-effective biomarker discovery tool for the personal genome era. Genome Med. 2013;5:41.
Gurwitz D, McLeod HL. Genome-wide studies in pharmacogenomics: harnessing the power of extreme phenotypes. Pharmacogenomics. 2013;14:337–9.
We gratefully thank Dr. Sébastien Jacques and Dr. Franck Letourneur at the Cochin Institute (Paris, France) for performing the micro-array experiments.
We thank the anonymous donors of the NLGIP biobank at Tel Aviv University, Israel, whose altruism and trust in biomedical research have made this study possible.
This study was funded by the BioIntelligence collaborative program for the Institut National de la Recherche Médicale (INSERM).
The authors declare that they have no competing interests.
Study conception and design: LC, DG, MAL. Performed the experiments: LC, AAFR, CM. Contributed to materials and analysis tools: ICP, FD, DG. Analysis and interpretation of data: LC, FD, DG, MAL. Drafting of manuscript: LC, DG, MAL. Critical revision: ICP, PB. Final approval: LC, AAFR, FD, CM, ICP, PB, DG, MAL.
Data S1. Methodology of microfluidic-based RT-qPCR assay. Table S2. Primers used for the microfluidic-based quantitative RT-qPCR assays. Table S3. Primers used for the copy number variation assays of EBV and mitochondrial DNA in lymphoblastoid cell lines. Table S4. List of the 210 genes significantly different at the P <0.001 level in the resistant compared to the sensitive cell lines. Table S5. List of the 40 genomic loci identified using a PAM on micro-array data. Table S6. List of genes involved in gene ontology terms. Table S7. Top Ingenuity® canonical pathways enriched by genes that were significantly differentially expressed in resistant cell lines. Table S8. Genes involved in the top Ingenuity® canonical pathways.
Correlations between growth inhibitions by thiopurine drugs among 53 lymphoblastoid cell lines. (A) Correlation between growth inhibition by 6-MP (2 μM) and by azathioprine (5 μM) (r spearman = 0.95; P <0.0001). (B) Correlation between growth inhibition by 6-MP (2 μM) and by 6-TG (0.5 μM) (r spearman = 0.81; P <0.0001). Data of growth inhibition by 6-TG (0.5 μM) are available only for 43 cell lines.
Principal component analysis for the 11 lymphoblastoid cell lines used for the transcriptomic analysis. #xxxx, cell line ID; R, resistant cell line; S, sensitive cell line.
Functional downstream CDKN1A network.