
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Regular Articles |
Department of Pediatrics, Kanazawa University School of Medicine, Kanazawa, Ishikawa 920-8641, Japan
| ABSTRACT |
|---|
|
|
|---|
T in the 3'-untranslated region of the human dihydrofolate reductase (DHFR) gene transcript, was identified in the study population of 37 patients with childhood leukemias/lymphomas and 83 healthy Japanese children. Frequencies of the DHFR 829C/C, 829C/T, and 829T/T genotypes were 83.8, 10.8, and 5.4%, respectively, in the cases and 74.7, 19.3, and 6.0% in the controls, showing no significant difference in genotype frequencies between the cases and controls. When determined by real-time quantitative reverse transcription-PCR analysis, the highest expression of the DHFR transcript was demonstrated in the samples with a DHFR 829T/T polymorphism (P < 0.001). Direct association of the presence of the SNP with methotrexate-related adverse events in each patient was not demonstrated in this limited analysis. These data suggest that the novel DHFR 829 polymorphism is associated with a positive role in gene expression and provide evidence of a functional SNP in the 3' regulatory region of the gene. | INTRODUCTION |
|---|
|
|
|---|
Recently developed technologies using high-throughput genotyping methods have accelerated SNP discovery by identifying the specific location in a gene at a frequency of >1% across a population (1, 2, 3, 4, 5) . It is necessary to build a comprehensive catalogue of SNPs in candidate genes in the human population to better understand the role of common genetic variants in common diseases. Recently, several studies characterizing SNPs on a large scale have been reported for this purpose (6, 7, 8, 9, 10) . In addition, one goal of pharmacogenetics must be to exploit all of the relevant SNP variations to improve the diagnosis and treatment of disease. In the near term, concerning adverse drug reactions, analyses of SNPs in the genes for the enzymes involved in the metabolic pathways of chemotherapeutic drugs may directly provide gene-based information that predicts drug response for the treatment of cancer patients.
DHFR (EC 1.5.1.3) catalyzes the reduction of dihydrofolate to THF, which rejoins the pool of active folate cofactors (11 , 12) . MTX acts as a tight-binding inhibitor of DHFR and is the most widely used antimetabolite in pediatric malignancies, being especially effective in the treatment of ALL, non-Hodgkins lymphoma, the histocytoses, and osteosarcoma (13) . MTX is also used in the treatment of rheumatoid arthritis, psoriasis, and autoimmune diseases, and for the prevention of graft-versus-host disease after transplantation. In actively proliferating tumor cells, inhibition of DHFR by MTX leads to accumulation of folates in the inactive dihydrofolate form, with a partial depletion of reduced folates. Similar to naturally occurring folates, which exist as polyglutamates within cells, intracellular MTX undergoes polyglutamation through the action of folylpolyglutamyl synthetase. MTX polyglutamates retained in cells subsequently inhibit DHFR and also enhance their inhibitory effects on folate-requiring enzymes of thymidylate and purine synthesis, thus contributing to MTX cytotoxicity (11 , 12) . In both experimental and clinical settings, increased levels of DHFR and a decreased DHFR binding affinity for MTX are common findings in cells manifesting an MTX-resistant phenotype, de novo or acquired (12 , 14) . The increased reductase activity may result from DHFR gene amplification, enhanced transcription, or efficient translation of DHFR messages. Several phenomena supporting these molecular alterations of the DHFR gene have been documented in MTX-resistant cells derived from clinical samples (12) . Thus, DHFR status in both target and nontarget cells has crucial roles in terms of sensitivity and resistance to MTX and could be involved in the occurrence of adverse reactions in individuals. In the present study, we therefore focused on identifying SNPs in the DHFR gene and evaluating the role of genetic variants in gene expression and drug response in a Japanese population.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Reverse Transcription-PCR and DNA Sequencing.
Total RNA was extracted from cryopreserved cell samples by a standard technique. The first-strand cDNA was generated with 2 µg of total RNA, random hexadeoxynucleotide primer (Takara, Tokyo, Japan), and RAV-2 reverse transcriptase (Takara). The oligonucleotide primers for the DHFR gene were 5'-CTGTCATGGTTGGTTCGC-3' (sense), with coordinates 38/55, and 5'-AAGCTTTTGGTATTTCCA-3' (antisense), with coordinates 962/979 [numbers correspond to the reference GenBank sequence (J100140)]. PCR amplification was performed in a DNA thermal cycler (GeneAmp PCR system 9600; Perkin-Elmer, Norwalk, CT) for 30 cycles (95°C for 30 s, 51°C for 30 s, and 72°C for 60 s). The PCR products were subsequently analyzed on 1.5% agarose gels (H14; Takara) stained with ethidium bromide. Gel-purified PCR products were sequenced directly in an automated sequencer (ABI model 310 genetic analyzer; Perkin-Elmer Applied Biosystems). Cycle-sequencing reactions were performed using ABI Prism Dye Terminator Sequencing Kits (Perkin-Elmer).
Sequence Analysis and Polymorphism Identification.
The ABI sequence software (DNA sequencing analysis, version 3.0; Perkin-Elmer Applied Biosystems) was used for sequence analysis. Sequencing traces in electropherograms were visually inspected, and variants were identified by comparison with the traces of the DHFR sequence relative to the reference GenBank sequence. Each variant position was confirmed by reamplification and resequencing of the variant site from the two strands. In addition, the DHFR 829C
T base pair substitution abolishes a TspR1 restriction site; thus, a PCR-RFLP assay was applied for the detection and confirmation of this variant. PCR products digested with 2 units of TspR1 restriction enzyme (New England BioLabs, Inc., Beverly, MA) were electrophoresed on a 2% agarose gel with ethidium bromide. The genotype according to DHFR 829 polymorphism was designated as 829C/C (prototype), 829C/T, or 829T/T.
Real-Time Quantitative PCR.
TaqMan probes were synthesized by Perkin-Elmer Japan Applied Biosystems (Chiba, Japan). The DHFR probe was 5'-6-carboxyfluorescein-CGTTCTTGCCGATGCCCATGTTCT-6-carboxy-tetramethylrhodamine-3' (nt 80103). The GAPDH probe used as an internal control was labeled at the 5' end with 27-dimethoxy-4,5-dichloro-6-carboxy-fluorescein and 6-carboxy-tetramethylrhodamine at the 3' end. The PCR primers for DHFR were designed as 5'-AAACTGCATCGTCGCTGTGTC-3' (sense; nt 5577) and 5'-ACCCATAATCACCAGATTCTGT-3' (antisense; nt 183204), and the primers for GAPDH were purchased from PE Biosystems.
Amplification reactions contained 5 µl of the first-strand cDNA reaction; 300 nM each PCR primer; 200 nM TaqMan probe; 5 µl of 10x PCR buffer A; 200 µM each of dATP, dCTP, and dGTP; 400 µM dUTP; 0.5 units of AmpErase uracil N-glycosylase; 5 mM MgCl2; and 1.25 units of AmpliTaq Gold DNA polymerase in a total volume of 50 µl. All components except the fluorogenic probes and amplification primers were supplied in a TaqMan PCR Core Reagent Kit (Perkin-Elmer). DNA amplifications were carried out in a 96-well reaction plate format in an ABI PRISM 7700 Sequence Detection System with a Sequence Detector V1.6.3 program (Perkin-Elmer Applied Biosystems). After a holding step at 50°C for 120 s to allow uracil N-glycosylase digestion in the reaction, the thermal cycling was initiated with a first denaturation step of 10 min at 95°C, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min.
To create a standard curve for relative quantification, a representative DNA sample without the SNP in the DHFR gene was chosen as a standard control, serially diluted in water (100, 10-1, and 10-2), and subjected to real-time quantitative PCR in triplicate. The dilution value (starting quantity) of the standard was plotted against the threshold cycle number (Ct) at which fluorescence first increased above background. The expression of the indicated gene in each sample was evaluated by this standard curve and shown as a value relative to the expression in the standard sample. These relative DHFR expression values were then corrected for the relative values obtained for GAPDH from the same DNA samples. Triplicate DHFR expression values were individually divided by the relative GAPDH values, and the mean was reported as DHFR/GAPDH normalized.
Statistical Analysis.
All data are presented as the mean ± SD or as a percentage. The Mann-Whitney test was used for differences between the gene expression values and a specific genotype. ORs and 95% CIs were used for estimating the risk of association with genotypes. Fishers exact probability tests were used to compute the risk estimates. All statistical analyses were performed using StatView (version 4.5; Abacus Concepts, Inc., Berkeley, CA). P < 0.05 was considered significant.
| RESULTS |
|---|
|
|
|---|
A substitution was observed in all alleles in this Japanese population and found identical to a candidate SNP proposed from the Cancer Genome Anatomy Project-Genetic Annotation Index. This 721A
T SNP may represent a population-specific polymorphism. A novel 829C
T variant (designated as DHFR 829) was identified in 34 of a total of 240 alleles. The representative sequence traces and a PCR-RFLP assay for validation were shown in Fig. 1, A and B
T was 14.2% in this Japanese population.
|
|
T SNP.
T SNP in the 3'-UTR, we quantitated the DHFR mRNA present in the samples from nonleukemic children as well as from leukemia/lymphoma patients. Fifteen randomly selected samples (10 from nonleukemic children and 5 from leukemia/lymphoma patients) with 829C/C, 12 random samples (8 from nonleukemic children and 4 from leukemia/lymphoma patients) with 829 C/T, and 10 (8 from nonleukemic children and 2 from leukemia/lymphoma patients) with 829T/T were analyzed for DHFR expression by real-time quantitative PCR analyses. The relative expression in the nonleukemic samples with C/C, C/T, and T/T was 0.22 ± 0.11 (median, 0.18), 0.40 ± 0.08 (0.41), and 2.54 ± 1.48 (2.53), respectively. The expression of DHFR transcript in the sample with 829T/T was significantly (P < 0.001) higher (
11-fold) than that with prototype (Fig. 2)
|
| DISCUSSION |
|---|
|
|
|---|
At least four pseudogenes derived from processed RNA molecules have been identified and located on separate chromosomes (19) . Despite highly homologous sequences, including the 3' portion of the DHFR pseudogenes to the DHFR gene (19) , a considerable number of single-nucleotide alterations existed around the SNP site, making it possible for us to discriminate between the functional DHFR transcripts and amplicons derived from the pseudogenes. The role of the 3'-UTR of mRNA is becoming clear and is seen to be as important as that of 5'-UTR in regulating gene expression. In addition to the key role of poly(A) tail, which protects the RNA chain from degradation by 3'-to-5' exonucleases, resulting in enhancement of translation, there are a number of motif sequences within the 3'-UTR that regulate mRNA stability, translational efficiency, or both (20) . Searching for the known regulatory UTR motifs, Day and Tuite (20) did not find any matches in the 3'-UTR sequences with or without the DHFR 829 polymorphism.
A comparison of the human and mouse DHFR 3'-UTR sequences revealed that only 100 nucleotides downstream from the terminator codon were conserved between the two species (18) . These structural features of the 3'-UTR of the human DHFR gene suggest that the downstream regions are less important for posttranscriptional gene regulation and are not part of the DHFR ancestral gene. However, the following evidence supports the possible involvement of the DHFR 829 polymorphism in the regulatory region within the 3'-UTR. Goldsmith and Cowan (21) constructed several DHFR minigenes with deletions, in which the first two polyadenylation sites and the DHFR 829 site turned out to be eliminated from the 3'-UTR. The level of DHFR induction and RNA transcripts by these constructs was reduced by 50% relative to the wild type, indicating the presence of regulatory elements in this region. Furthermore, studies using a human DHFR minigene construct supported the regulatory role of DHFR UTRs in gene expression during MTX exposure (22) . Recent studies on 3'-UTR-mediated regulation have focused on the interaction between trans-acting proteins and cis-acting mRNA sequences (20 , 23) . A specific RNA-protein interaction associated with the 3'-UTR regulatory elements would likely be involved in this enhancement, although the mechanism(s), e.g., increased mRNA stability, efficient translation, interaction of regulatory elements of 5'- and 3'-UTRs, or alteration of nucleo-cytoplasmic transport, of enhanced expression of DHFR mRNA associated with the DHFR 829 polymorphism need(s) to be clarified.
It would be anticipated that increased DHFR reduces MTX cytotoxicity in normal cells while conferring resistance in target cells. However, when we compared DHFR genotypes with MTX-related adverse events and patient outcomes, increased DHFR expression associated with the 829T/T genotype did not provide direct evidence of the clinical significance of the polymorphism (data not shown). This was possibly in part because MTX cytotoxicity is multifactorial. MTX-induced inhibition involves several folate-dependent enzymes in both de novo purine and thymidylate synthesis (12) . The coincidence and grades of MTX-related adverse reactions depend on the extent of the status of DHFR as well as the enzymes involved in the metabolic pathways of folates and their differential expression between normal and malignant cells. Therefore, the expression of these enzymes, including thymidylate synthetase, glycinamide ribonucleotide and aminoimidazole carboxamide ribonucleotide transformylases, folylpolyglutamyl synthetase, and folylpolyglutamate hydrolase, in both target and nontarget cells is the critical determinant of drug response. A high-density oligonucleotide array (DNA chip)-based analysis could be applied for this purpose, but it is still experimental in clinical practice. Profiling the combination of SNPs in particular reactions would help to circumvent unfavorable events and provide benefit for the treatment of patients. An additional, larger study with accumulation of clinical samples is required to find conclusive associations between the SNPs and drug-related adverse events.
A recent study by Skibola et al. (9) demonstrated that polymorphisms in the gene that encodes MTHFR, an enzyme that reduces methylene-THF to methyl-THF, decrease the risk of adult ALL. MTHFR is thought to have a protective role in leukemogenesis, with decreased MTHFR activity attributable to cSNPs consequently reducing the misincorporation of uracil into DNA, a process that is mediated by enhanced availability of methylene-THF in the DNA synthesis pathway (9) . Theoretically, MTX polyglutamates inhibit de novo thymidylate synthesis, resulting in enhanced uracil nucleotide misincorporation as a result of an increase in the dUTP pool (12) ; therefore, enhanced DHFR activity may anticipate a preventable function similar to reduced MTHFR activity. However, our results showed no evidence of risk association (either predictive or protective) between DHFR genotype and the disease. This is attributable in part to the difference in the metabolic fate of folates affected by the two enzymes (9 , 12) . Notably, the same MTHFR polymorphism failed to show evidence of a protective effect against AML (9) , indicating that cells may not be affected equally because of different biological features involved in the different leukemia types.
Efforts to analyze drug targets for polymorphisms related to drug treatment response are being made because gene-based information will make it possible to establish individual therapy as well as to improve the diagnosis of disease. Our data provide both evidence of the role of a genetic variant in the 3'-UTR in gene regulation and the framework to establish an additional therapeutic strategy based on genetically determined drug responsiveness in risk-stratified protocols for the treatment of childhood malignancies.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
1 This work was supported by Grants-in-Aid from the Ministry of Education, Science, Sports, and Culture of Japan (09670793, 11470170, and 11877128). ![]()
2 To whom requests for reprints should be addressed, at Department of Pediatrics, Kanazawa University School of Medicine, 13-1 Takara-machi, Kanazawa, Ishikawa 920-8641, Japan. Phone: 81-76-265-2314; Fax 81-76-262-1866; E-mail: saikawa@ped.m.kanazawa-u.ac.jp. ![]()
3 The abbreviations used are: SNP, single-nucleotide polymorphism; cSNP, SNP within the coding region; DHFR, dihydrofolate reductase; THF, tetrahydrofolate; MTX, methotrexate; ALL, acute lymphocytic leukemia; AML, acute myelogenous leukemia; nt, nucleotide, GAPDH, glyceraldehyde-3-phosphate dehydrogenase; UTR, untranslated region; OR, odds ratio; CI, confidence interval; MTHFR, methylenetetrahydrofolate reductase. ![]()
Received 12/27/00; revised 3/27/01; accepted 4/10/01.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. Fabbri, N. Valeri, and G. A. Calin MicroRNAs and genomic variations: from Proteus tricks to Prometheus gift Carcinogenesis, June 1, 2009; 30(6): 912 - 917. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. L. Hider, I. N. Bruce, and W. Thomson The pharmacogenetics of methotrexate Rheumatology, October 1, 2007; 46(10): 1520 - 1524. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. J. Mishra, R. Humeniuk, P. J. Mishra, G. S. A. Longo-Sorbello, D. Banerjee, and J. R. Bertino A miR-24 microRNA binding-site polymorphism in dihydrofolate reductase gene leads to methotrexate resistance PNAS, August 14, 2007; 104(33): 13513 - 13518. [Abstract] [Full Text] [PDF] |
||||
![]() |
U. Lim, S. S. Wang, P. Hartge, W. Cozen, L. E. Kelemen, S. Chanock, S. Davis, A. Blair, M. Schenk, N. Rothman, et al. Gene-nutrient interactions among determinants of folate and one-carbon metabolism on the risk of non-Hodgkin lymphoma: NCI-SEER Case-Control Study Blood, April 1, 2007; 109(7): 3050 - 3059. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. M. Troen, B. Mitchell, B. Sorensen, M. H. Wener, A. Johnston, B. Wood, J. Selhub, A. McTiernan, Y. Yasui, E. Oral, et al. Unmetabolized Folic Acid in Plasma Is Associated with Reduced Natural Killer Cell Cytotoxicity among Postmenopausal Women J. Nutr., January 1, 2006; 136(1): 189 - 194. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. L. Carroll, D. Bhojwani, D.-J. Min, E. Raetz, M. Relling, S. Davies, J. R. Downing, C. L. Willman, and J. C. Reed Pediatric Acute Lymphoblastic Leukemia Hematology, January 1, 2003; 2003(1): 102 - 131. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A. Chabner Cytotoxic Agents in the Era of Molecular Targets and Genomics Oncologist, August 1, 2002; 7(90003): 34 - 41. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |