![]() |
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Division of Medicinal Safety Science (N.K., K.K., M.T., R.H.), Division of Environmental Chemistry (H.J., T.T-K.), and Project Team for Pharmacogenetics (Y.S., M.S., J.S.), National Institute of Health Sciences, Tokyo, Japan
(Received August 10, 2004; accepted November 23, 2004)
| Abstract |
|---|
|
|
|---|
The wide interindividual variability in SN-38G formation in hepatic tissues is known and has been shown to correlate with UGT1A1 genetic polymorphisms (Iyer et al., 1999
). To date, a wide variety of allelic polymorphisms of human UGT1A1 have been reported, most of which lead to a partial or complete deficiency of enzyme activity and are associated with diseases, such as the Criegler-Najjar (CN) syndrome types I and II or Gilbert's syndrome, with hyperbilirubinemia and jaundice as main symptoms (http://som.flinders.edu.au/FUSA/ClinPharm/UGT/1A1aalleles.html; Mackenzie et al., 1997
; Tukey and Strassburg, 2000
). To understand these syndromes and their possible genetic causes, polymorphisms in the TATA box have been the most extensively studied. An inverse relationship between the number of TA repeats (five to eight) and transcriptional activity has been reported (Beutler et al., 1998
). It has been reported that serum bilirubin levels in individuals with seven TA repeats are higher than the levels in individuals without a seventh repeat (Bosma et al., 1995
; Raijmakers et al., 2000
). Patients suffering from Gilbert's syndrome, CN syndrome, or prolonged neonatal jaundice also often carry seven or eight TA repeats (Iolascon et al., 1999
; Burchell et al., 2000
). The UGT1A1 promoter polymorphisms are also considered a determinant of CPT-11 disposition and toxicity (Iyer et al., 1999
; Ando et al., 2000
; Iyer et al., 2002
). Other variations often associated with neonatal hyperbilirubinemia, CN syndromes, Gilbert's disease, or irinotecan toxicity are 211G>A (G71R), 1456T>G (Y486D), and 686C>A (P229Q) (Ando et al., 2000
; Kraemer and Klinker, 2002
; Yamamoto et al., 2002
) (In this report, the A of the translational start codon of the cDNA is designated 1.) Recently, the possibility of a synergistic effect of 3'-untranslated region variations with 3279T>G on the bilirubin level and SN-38 pharmacokinetics was also proposed (Sai et al., 2004
). The single nucleotide polymorphism (SNP) 3279T>G is located in the phenobarbital-responsive enhancer module, which is involved in activation of UGT1A1 transcription by the constitutive androstane receptor (Sugatani et al., 2002
).
A recent pharmacogenetic study has suggested that there are advantages in the use of haplotypes rather than individual SNPs to investigate the association between genotypes and phenotypes (Judson et al., 2000
). Racial variability in haplotype frequencies of UGT1A1 may contribute to an interpretation of the racial diversity of pharmacokinetics and toxicity of drugs metabolized by the enzyme. Recently, Sai et al. (2004
) revealed the haplotype structure of UGT1A1 in Japanese cancer patients using 3279T>G, TA repeats in the TATA box, 211G>A, 686C>A, and three SNPs (1813C>T, 1941C>G, and 2042C>G) in the 3'-untranslated region in exon 5 as markers. Therefore, we studied in this report the differences in the haplotype frequencies of UGT1A1 among African-American, Caucasian, and Japanese populations using genomic DNA obtained from 150 individuals in each population, assessing the genotypes of these markers with dideoxy sequencing and pyrosequencing methods previously reported by Saeki et al. (2003
). The function of a novel SNP 686C>T (P229L) that was found in an African-American individual during this study was also investigated by utilizing a heterologous expression system with COS-1 cells.
| Materials and Methods |
|---|
|
|
|---|
Peripheral Blood Samples. One hundred fifty peripheral blood samples were obtained from healthy Japanese volunteers, 100 of which were kindly provided by Professor Ichiro Ieiri from the Tottori University Hospital with permission of the ethics committees of the Tottori University Faculty of Medicine. Written informed consent was obtained from all participants. Peripheral blood samples of healthy Caucasian and African-American volunteers (150 each) were purchased from the Tennessee Blood Service Corporation (Memphis, TN). The ethics committee of the National Institute of Health Sciences approved this study.
UGT1A1 Genotyping by Pyrosequencing and Dideoxy Sequencing. Genotyping of UGT1A1 using Pyrosequence methods was performed as previously reported (Saeki et al., 2003
). Briefly, genomic DNA (1015 ng) was amplified by Ex-Taq (1 U; Takara Shuzo, Ootsu, Japan) with a specific primer pair in which one of the primers was biotinylated. The primers for amplification and sequencing used for detection of polymorphisms 3279T>G, the TA repeat, 211G>A (G71R), and 686C>A/T (P229Q/L) have been shown previously (Saeki et al., 2003
). Genotyping of 1941C>G for the Japanese samples was performed using an amplification primer pair (biotin-ATTTGAATATG-TATCGTGCCC and CATTCATTCATTTCACCTACACT) and a sequencing primer (CAGTAGGGGCAGC). Generation of the single-stranded fragment and annealing of the sequencing primers have been described previously (Saeki et al., 2003
). Genotypes were determined using the PSQ 96MA (Biostage AB, Uppsala, Sweden) and the PSQ 96 SNP reagent set (Biostage AB). Genotyping of 1813C>T, 1941C>G, and 2042C>G for African-American and Caucasian samples as well as the detection of a genotype 686C>T in the African-American individual were done by dideoxy sequencing as reported previously (Saeki et al., 2002
; Sai et al., 2004
).
Haplotype Analysis. Genotyping was successfully determined for 150 Japanese, 147 Caucasians, and 148 African-Americans. Diplotypes (combinations of haplotypes) in both blocks 1 and 2 were inferred separately for each population by PHASE version 2.0 (Stephens et al., 2001
; Stephens et al., 2003).
Expression and Enzyme Assay of Wild-Type and Variant UGT1A1s. Expression of wild-type and variant UGT1A1s in COS-1 cells was carried out as described previously (Jinno et al., 2003b
) with a minor modification; attB-flanked UGT1A1 cDNA was amplified by the two-step attB adaptor polymerase chain reaction (PCR) (Jinno et al., 2003a
) using pcDNA3.1-UGT1A1/WT (Jinno et al., 2003b
) as a template and gene-specific primers, 5'-AAAAAGCAGGCTGCAAAGGCGCCATGGCTGT-3' and 5'-AGAAA-GCTGGGTCTCAATGGGTCTTGGATTTGTGGG-3'. The resulting PCR fragment was cloned into the pDONR201 vector (Invitrogen, Carlsbad, CA). The 686C>T mutation was introduced into the wild-type UGT1A1 cDNA clone in pDONR201, using a QuickChange multisite-directed mutagenesis kit (Stratagene, La Jolla, CA) with the 5'-phosphorylated oligonucleotide primer 5'-phospho-GCGACGTGGTTTATTCCCTGTATGCAACCCTTGCCTC-3' (the nucleotide base changed is underlined). Subcloning of each UGT1A1 fragment from pDONR201 into pcDNA-DEST40, a mammalian expression vector, was performed by the Gateway LR reaction.
Expression levels of UGT1A1 mRNA and protein were determined as described previously (Jinno et al., 2003b
). SN-38 glucuronidation activity of the wild-type and P229L variant UGT1A1 was assayed according to the method of Hanioka et al. (2001
).
Data Analysis. The
2 test in the SAS Preclinical Package (SAS Institute Japan Ltd., Tokyo, Japan) was applied to s x r contingency tables for comparison of genotype frequencies and analysis of haplotype distribution patterns among the three ethnic groups. In vitro Michaelis-Menten kinetic parameters were estimated by nonlinear regression analysis using Prism version 3.0 (GraphPad Software Inc., San Diego, CA). The average values of the kinetic parameters were calculated using results from three independent preparations. The t test supplied in Prism version 3.0 was applied to the comparison of the average values of protein expression and mRNA levels between wild-type and variant UGT1A1 at a significance level of 0.05.
| Results |
|---|
|
|
|---|
2 test). Regarding the variation in the TATA box, the allele TA6 was most common in the Japanese population, whereas frequencies of TA6 and TA7 (UGT1A1*28) were equally frequent in African-Americans. The allele frequencies of TA5 (UGT1A1*36) or TA8 (UGT1A1*37) in African-Americans were about 5%, but no carriers with TA5 or TA8 were found in the Japanese. The allele distribution patterns for the wild-type and three variants in the TATA box of the three ethnic groups were significantly different (p < 0.0001 by the
2 test). The variant 211G>A (G71R) (UGT1A1*6) was frequently found in the Japanese population, but not in the Caucasians population (only two heterozygous subject), and none were found in African-Americans. The variant 686C>A (P229Q) (UGT1A1*27) was very rare in all the ethnic groups, and only one heterozygous carrier was found in the Japanese group. A novel variation in the same position as the variant 686C>A, but with a different nucleotide change (686C>T), was found in an African-American. This novel variant led to an amino acid change from proline to leucine (P229L). Pyrograms for 686C/C, 686C/A, and 686 C/T are compared in Fig. 1A and the nucleotide change of the novel variant was confirmed by the dideoxy sequencing method (Fig. 1B). This is the first report detailing allele frequencies of three SNPs in the 3'-untranslated region in exon 5, 1813C>T, 1941C>G and 2042C>G in African-Americans and Caucasians. The frequencies of 1813C>T and 2042C>G in the Japanese were not determined in this study, because these alleles were reported to be in complete association with 1941C>G in a Japanese population without exception (Sai et al., 2004
2 test).
|
|
Differences in Haplotype Frequencies. Sai et al. (2004
) previously reported that UGT1A1 could be divided into two blocks, with the transcription-regulating and promoter regions and exon 1 in block 1, and exons 2 to 5 in block 2, according to linkage disequilibrium analysis. Therefore, haplotype/diplotype analysis was performed using four marker variations in block 1 (3279T>G, 211G>A, and 686C>A, and the TATA box), and three marker variations in block 2 (1813C>T, 1941C>G, and 2042C>G). The diplotype configurations (combinations of haplotypes) were estimated for each subject by PHASE software (Stephens et al., 2001
; Stephens et al., 2003).
Concerning block 1, the diplotype configurations were inferred with greater than 0.99 certainties for 145 Japanese, 144 Caucasians, and 147 African-Americans. The haplotypes identified are summarized in Table 2 along with their frequencies for each population, where the 14 subjects with ambiguous diplotypes were excluded. We followed the haplotype nomenclature previously reported (Sai et al., 2004
). One Japanese subject carrying heterozygous 686C>A also carried homozygous TA7. Thus, the association of variation 686C>A with TA7 was also confirmed in this study, as previously reported by Huang et al. (2000
) and Sai et al. (2004
). A new haplotype, designated *6d, having variations in both positions 3279 and 211, was identified that was not found in the previous study (Sai et al., 2004
). We could not determine whether the novel variation 686T and the *28 allele were on the same chromosome, because the subject with 686T carried heterozygous TA repeats, TA6/TA7.
|
In block 1, haplotype *1a was predominant in the Japanese population (0.61 in frequency). In contrast, its frequencies in Caucasians and African-Americans were 0.45 and 0.15, respectively. Major haplotypes in African-Americans were *28b and *60a. Haplotype frequencies of *1a and *28b in Caucasians were comparable. Thus, the haplotype distribution patterns in block 1 for the individual populations were significantly different (p < 0.0001 by the
2 test).
In block 2, the diplotype configurations using the three SNPs were inferred with certainties greater than 0.98 for all Caucasians and African-Americans (150 each). The haplotypes and their frequencies are summarized in Table 3. Four novel haplotypes, *IC, *ID, *IE, and *IF, were identified in the Caucasian and African-American populations, although *IE was not found in the latter. These haplotypes were not found in the Japanese, as shown in a previous study (Sai et al., 2004
). Although the three SNPS were reported to be completely associated with each other in the Japanese (Sai et al., 2004
), they were not always linked with each other in the rest of the populations. However, it is noteworthy that 1941G allele was always associated with 1813T and 2042G alleles except for the haplotype *IE, found only in two Caucasians. Haplotype *IA was predominant for all ethnic groups. The second major haplotype was *IB in both the Japanese and Caucasians. However, the frequencies of *IB and *IC were similar in African-Americans (0.183 and 0.163, respectively). The haplotype *IC was also found in Caucasians (0.06). Thus, the haplotype distribution patterns in block 2 for the three populations were also significantly different (p < 0.0001 by the
2 test).
|
Expression and SN-38 Glucuronidation of a Novel Variant P229L Compared with Wild-Type UGT1A1. The relative expression level of the novel variant UGT1A1/P229L in the membrane fraction of COS-1 cells was determined by Western blotting using a polyclonal anti-human UGT1A antibody (Fig. 2). Because two bands, one at the same position as the wild-type and the other at a higher molecular weight, were detected for P229L, the total chemiluminescence of the two bands was used to calculate the protein expression level for P229L. The protein expression level of P229L was approximately 60% of the wild type, and this difference was statistically significant by the t test (p = 0.044). The UGT1A1 mRNA expression level was then measured by real-time reverse transcription-PCR (RT-PCR) using SYBR Green. As shown in Fig. 3, however, no significant difference in the mRNA level was detected among the COS-1 cells transfected with the expression plasmids carrying wild-type and P229L UGT1A1 cDNAs.
|
|
The glucuronidation activity of SN-38 by P229L expressed in COS-1 cells was compared with that of the wild-type under 11 substrate concentrations ranging from 2.5 to 150 µM. The representative curves of the Michaelis-Menten kinetics are shown in Fig. 4, and the estimated apparent kinetic parameters (Km, Vmax, and Vmax/Km) are summarized in Table 4. Vmax values normalized to the expression levels are also shown. Wild-type UGT1A1 catalyzed SN-38 glucuronidation with an average apparent Km value of 8.67 µM, whereas P229L catalyzed SN-38 glucuronidation with a Km value of 37.6 µM. The average Vmax values were 71.2 and 5.27 pmol/min/mg membrane protein for the wild-type and P229L, respectively. The average intrinsic clearances of SN-38 by glucuronidation (Vmax/Km) normalized to expressed UGT1A1 protein levels were 8.26 and 0.24 µl/min/mg protein for the wild-type and P229L, respectively.
|
|
| Discussion |
|---|
|
|
|---|
Innocenti et al. (2002
) reported on the haplotype frequencies of UGT1A1 in African-Americans and Caucasians using 3279T>G, the TA repeat, and other SNPs in the transcription-regulating region. Although they did not use the marker SNPs 211G>A and 686C>A (likely due to their rarity in both populations), frequencies of the *28, *36, *37, and *60 haplotype groups in African-Americans can be calculated according to their data as 0.35, 0.04, 0.12, and 0.33, and the calculated frequencies in Caucasian as 0.36, 0.01, 0.01, and 0.09. These frequencies were also comparable to our data. We found two individuals carrying 211G>A in the Caucasian group; however, we could not infer haplotypes for these individuals with high certainty. The observed haplotype *28 group frequency in our study of 0.100 (*28b plus *28c) was slightly lower than the reported values of 0.131 in Japanese patients (Sai et al., 2004
) and 0.127 in Koreans (Ki et al., 2003
). However, the previously reported values are included in the 95% confidence interval of our data. The observed haplotype *6 group frequency (*6a plus *6b) was also comparable to the previously reported data in Japanese and Koreans. A strong association of TA7 with 3279G was suggested for Japanese cancer patients (Sai et al., 2004
). This association was also observed in all populations measured in this study. Other variations, TA5 and TA8 in the TATA box, also had linkage with 3279G.
Recent studies indicate that 3'-untranslated regions may regulate mRNA stability (Day and Tuite, 1998
; Conne et al., 2000
). Exon 1 is unique for each member of the UGT1A subfamily, whereas exons 2 to 5 are common to all members of the subfamily. The variants in the 3'-untranslated region of UGT1A1 in exon 5, therefore, could have various effects on all enzymes of the UGT1A subfamily. However, until recently, their physiological effects have remained unknown. Acuna et al. (2003) reported a protective effect of the variant 1941C>G on liver transaminase levels caused by tolcapone, for which glucuronidation (mainly by UGT1A9) is one metabolic pathway (Acuna et al., 2003). The SNP 1941C>G is a marker for haplotype *IB in block 2 as shown in Table 3. Sai et al. (2004
) reported the *IB-dependent decreasing trend of the AUC ratio (SN-38G/SN-38) and an increase in serum total bilirubin levels. Highly differential haplotype distributions in both blocks 1 and 2 observed in different ethnic groups may cause different toxicity profiles during therapy using drugs, such as CPT-11, that are metabolized by UGT1A1.
A novel SNP 686C>T was found in an African-American sample in this study, leading to an amino acid change, P229L, which is different from the known SNP 686C>A (P229Q). Studies using a larger sample size may be necessary to elucidate an accurate frequency. The known SNP 686C>A (P229Q) was reported to be related to Gilbert's syndrome and to have very low bilirubin glucuronidation activity (Koiwai et al., 1995
). However, Jinno et al. (2003a
,b
) reported that protein expression and mRNA levels of P229Q expressed in COS-1 cells were comparable to those of the wild type, and the decrease in its SN-38 glucuronidation was marginal. Since the association of P229Q with TA7 has been suggested previously (Huang et al., 2000
; Sai et al., 2004
) and confirmed in this study, hyperbilirubinemia observed in Japanese and Taiwanese patients carrying the P229Q variant can be mainly attributed to the TA7 variation. On the other hand, the novel variant P229L expressed in COS-1 cells was found to have extremely low SN-38 glucuronidation efficacy in this study (less than 3% of that for the wild type). The low glucuronidation activity of P229L was not caused by its low transcription levels, because a significant difference in the mRNA level was not observed between the wild type and variant. The low UGT1A1/P229L expression level in COS-1 cells suggests that this enzyme is unstable compared with the wild-type enzyme. Moreover, a band shifted to a higher molecular weight than the wild-type UGT1A1 was detected by a rabbit anti-human UGT1A antibody in a Western blot. Such additional bands were also observed in Western blots of variant UGT1A1 having extremely low enzyme activities, such as F83L and Y486D in our group (Jinno et al., in press). The band shift to a higher molecular weight in SDS-polyacrylamide gel electrophoresis is considered to be the result of products with different post-translational modifications such as phosphorylation or glycosylation, which may occur due to conformational changes caused by amino acid substitution. This may also contribute to the extremely low catalytic activity of the variant enzyme compared with that of the wild type as shown in Table 4. Although the clinical significance of 229L remains to be confirmed, this novel SNP may be a potential cause of diseases such as hyperbilirubinemia or jaundice and may cause severe adverse effects after administration of CPT-11.
| Acknowledgments |
|---|
| Footnotes |
|---|
Article, publication date, and citation information can be found at http://dmd.aspetjournals.org.
ABBREVIATIONS: SN-38, 7-ethyl-10-hydroxycamptothecin; CPT-11, 7-ethyl-10-[4-(1-piperidino)-1-piperidino]carbonyloxycamptothecin; SN-38G, 7-ethyl-10-hydroxycamptothecin glucuronide; CN, Crigler-Najjar; SNP, single nucleotide polymorphism; PCR, polymerase chain reaction; RT-PCR, reverse transcription-polymerase chain reaction.
Address correspondence to: Dr. Nahoko Kaniwa, Division of Medicinal Safety Science, National Institute of Health Sciences, 1-18-1 Kamiyoga, Setagaya-ku, Tokyo 158-8501, Japan. E-mail: nkaniwa{at}nihs.go.jp
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. L. Hong, D. Huo, H.-J. Kim, Q. Niu, D. L. Fackenthal, S. A. Cummings, E. M. John, D. W. West, A. S. Whittemore, S. Das, et al. UDP-Glucuronosyltransferase 1A1 Gene Polymorphisms and Total Bilirubin Levels in an Ethnically Diverse Cohort of Women Drug Metab. Dispos., August 1, 2007; 35(8): 1254 - 1261. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Saeki, Y. Saito, K. Sai, K. Maekawa, N. Kaniwa, J.-i. Sawada, M. Kawamoto, A. Saito, and N. Kamatani A Combinatorial Haplotype of the UDP-Glucuronosyltransferase 1A1 Gene (#60-#IB) Increases Total Bilirubin Concentrations in Japanese Volunteers Clin. Chem., February 1, 2007; 53(2): 356 - 358. [Full Text] [PDF] |
||||
![]() |
N. Risch Dissecting Racial and Ethnic Differences N. Engl. J. Med., January 26, 2006; 354(4): 408 - 411. [Full Text] [PDF] |
||||
![]() |
W. Sadee and Z. Dai Pharmacogenetics/genomics and personalized medicine Hum. Mol. Genet., October 15, 2005; 14(suppl_2): R207 - R214. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||