|
|
|
|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Departments of Pediatrics (M.A.K., E.T.C., R.N.H.) and Pharmacology and Toxicology (M.A.K., R.N.H.), Medical College of Wisconsin and Children's Research Institute (M.A.K., E.T.C., R.N.H.), Children's Hospital and Health System, Milwaukee, Wisconsin; and Departments of Medicinal Chemistry (A.E.R.) and Genome Sciences (M.J.R.), University of Washington, Seattle, Washington
Received for publication December 7, 2007.
Accepted for publication February 28, 2008.
| Abstract |
|---|
|
|
|---|
CYP2C9 is a major human cytochrome P450 enzyme that accounts for approximately 20% of the total cytochrome P450 protein content in adult human liver (Shimada et al., 1994
), and it is responsible for the metabolism of approximately 16% of clinically used drugs cleared by oxidative pathways (Williams et al., 2004
). Examples of small-molecular-weight therapeutics for which CYP2C9 is important for disposition include the anticoagulant warfarin; the antidiabetic agents tolbutamide and glipizide; the anticonvulsant phenytoin; the antihypertensive losartan; the antidepressant fluoxetine; and several nonsteroidal anti-inflammatory drugs, such as ibuprofen, diclofenac, and celecoxib (Klose et al., 1998
; Miners and Birkett, 1998
; Davies et al., 2000
). Over the past several years, multiple CYP2C9 single-nucleotide polymorphisms (SNPs) resulting in amino acid changes have been identified. More than 30 allelic variants are currently listed on the cytochrome P450 allele web site (http://www.cypalleles.ki.se); however, not all have been fully characterized. The most common CYP2C9 alleles include CYP2C9*1A (reference allele), CYP2C9*2 (three haplotype variants, all containing the g.3608C>T, R144C SNP), and CYP2C9*3 (two haplotype variants, all containing the g.42614A>C, I359L SNP). Less common are CYP2C9*4 (g.42615T>C, I359T) and CYP2C9*5 (g.42619C>G, D360E), and CYP2C9*6 (g.10601delA). The CYP2C9*2 and CYP2C9*3 variant alleles encode enzymes exhibiting significantly lower intrinsic clearance both in vivo and in vitro (for review, see Lee et al., 2002
).
As evidenced above, extensive research has been conducted to identify and characterize SNPs within the CYP2C9 coding region. However, much less is known about other variant sites, such as polymorphic promoter sequences, that also might contribute to observed interindividual differences in CYP2C9 expression. In a population study of Japanese epileptic patients, seven polymorphic sites were identified within the first 2000 bp upstream of the CYP2C9 transcription start site (Shintani et al., 2001
). Several of the identified SNPs were in linkage disequilibrium, resulting in six unique haplotypes. One of the haplotypes (g.-1912T>C, g.-1885C>, g.-1538G>A, g.-1189C>T, and g.-982G>A) exhibited a 60% reduction in in vitro promoter activity. Furthermore, an association between estimated phenytoin intrinsic clearance and promoter activity was observed with respect to the identified CYP2C9 upstream haplotypes (Shintani et al., 2001
). However, because the promoter SNPs defining this pattern also were in linkage disequilibrium with the SNP defining CYP2C9*3, it is unclear what impact these variants might have on CYP2C9 phenotype in vivo. Similar findings were reported in a population study of Japanese and non-Latino white subjects using warfarin (Takahashi et al., 2004
). However, when differences in S-warfarin CYP2C9 clearance were compared between Japanese and non-Latino white patients having reference alleles in both promoter (up to position -2100) and coding regions, the Japanese patients exhibited significantly greater intrinsic clearance than the non-Latino white patients. These data not only suggest the likelihood of significant interpopulation genetic differences but also the possible presence of additional important CYP2C9 promoter elements and variants upstream of position -2100. Such a conclusion is not only supported by the study of Takahashi et al. (2004
) but also has precedence based on our knowledge of other cytochrome P450 genes (Martinez-Jiménez et al., 2005
).
Given the data supporting the likely presence of important CYP2C9 regulatory polymorphisms, the evidence for interpopulation differences in both the presence and frequency of genetic variants, and the paucity of pharmacogenetic studies in the Hispanic population, the present study was designed to begin addressing this knowledge gap.
| Materials and Methods |
|---|
|
|
|---|
expression plasmid pCMVHNF-4
(Stoffel and Duncan, 1997
expression plasmid pBJ5HNF-1
(Kuo et al., 1990
(DCoH) pBJ5DCoH (Mendel et al., 1991Subjects. After consent, blood samples were collected from 193 women of self-reported Hispanic background who traced their Mexican ancestry back a minimum of two generations. Volunteers were recruited after admission to the labor and delivery units of Provena Saint Therese Medical Center or Victory Memorial Hospital (Waukegan, IL). DNA was extracted using the QIAmp DNA blood MIDI kit and stored at 4°C. This research protocol was approved by all involved Institutional Review Boards.
DNA Amplification for SNP Discovery. All PCR DNA amplification primers were designed using OLIGO version 6.45 (Molecular Biology Insights, Cascade, CO). Each primer pair was designed to have similar melting temperatures to facilitate high-throughput processing using a 96-well format. Primer sequences are provided in Supplement Table S1.
Common CYP2C9 genetic variants were identified by sequencing the approximate first 10 kbp of CYP2C9 5'-flanking sequence in each sample from the Coriell 24PDR. Templates were prepared by PCR DNA amplification using 20 to 25 ng of genomic DNA in a 20-µl reaction volume containing 0.2 mM each deoxyribonucleotide triphosphate, 0.5 µM each primer, and 1.25 units of Herculase high-fidelity polymerase (Stratagene). Cycles were as follows: denaturation at 92°C for 40 s and annealing and extension at 56°C to 62°C for 1 min for a total of 30 cycles. All amplification reactions included an initial 10-s hold at 94°C and a final 5-min hold at 75°C. After amplification, 4 µl of ExoSap-IT (United States Biochemical Corp.) was added, and the reaction was incubated at 37°C for 30 min to remove unincorporated deoxy- and dideoxyribonucleotide triphosphates and primers.
CYP2C9 Sequence Analysis. Sequence analysis was performed using 60 to 100 fmol of amplicon with 8 µl of Quick Start Sequencing reagent (Beckman Coulter) in a PTC-225 Peltier thermal cycler per manufacturer's recommendations (MJ Research, Watertown, MA). Analysis was performed on both DNA strands. SNPs identified on a single allele were verified by a repeat analysis of an independently generated amplicon. In all instances, the coordinates of the identified SNPs follow the recommendations of the HUGO nomenclature working group, wherein the "A" of the ATG start codon is assigned +1 and uses contig NT_030059 [GenBank] .12, build 36.1 as a reference. For purposes of clarity, the coordinates of all CYP2C9 DNA fragments also follow the same convention. Immediately after the sequencing reaction, unincorporated dideoxyribonucleotide triphosphates and primer were removed using CleanSeq magnetic beads and reagents (Agencourt) before analysis of samples by capillary electrophoresis in a CEQ8000 genetic analysis system (Beckman Coulter).
Genotyping. After reviewing all of the identified CYP2C9 SNPs, there were five unique clusters and eight remaining SNPs that fell outside of the clusters. Five PCR primer pairs were designed to amplify the DNA sequences containing the SNP clusters, resulting in amplicons ranging in size from 269 to 853 bp (Supplemental Table S1). Two CYP2C9 fragments containing the remaining eight sequence variants were amplified, resulting in amplicons of 1.6 and 2.3 kbp. DNA amplifications were performed as described above using 20 to 25 ng of genomic DNA as template. DNA amplification of sequences containing exons 3, 5, and 7 were generated for use as templates for genotyping of the previously identified SNP tags for the CYP2C9*2, CYP2C9*3, CYP2C9*4, CYP2C9*5, and CYP2C9*6 alleles (for review, see Lee et al., 2002
) (Supplemental Table S2).
The frequencies of previously identified CYP2C9 promoter variants (Veenstra et al., 2005
) as well as any new variants identified within the discovery phase of this project were determined within the Hispanic population using a combination of DNA sequencing (all SNPs that fell within the five unique cluster sets) and multiplexed single-base extension (SBE) reactions (Lindblad-Toh et al., 2000
) (eight remaining SNPs as well as the SNP tags for the CYP2C9*2, CYP2C9*3, CYP2C9*4, CYP2C9*5, and CYP2C9*6 alleles) (for review, see Lee et al., 2002
). DNA sequencing was performed as described above. For genotyping by SBE, two muliplex primer sets were designed, consisting of four primers in each set, using OLIGO version 6.45. The CYP2C9 -1537G>A variant was not identified in the initial SNP discovery. However, because of the previously reported linkage with the CYP2C9*3 haplotype (Veenstra et al., 2005
), the frequency of this SNP subsequently was determined in the Hispanic study population using SBE. SBE reactions were performed in a 20-µl reaction volume consisting of a mixture of all four dye-labeled dideoxynucleotide terminators and the supplied proprietary polymerase (Beckman Coulter), 20 fmol of template, and approximately 1 to 10 pmol of each primer. Individual primer concentrations were optimized to ensure adequate signal intensities for each primer within the multiplex. SBE reactions were as follows: denaturation at 96°C for 10 s, annealing at 50°C for 5 s, and extension at 72°C for 30 s for a total of 25 cycles. Unincorporated dideoxynucleotides were eliminated by adding one unit of shrimp alkaline phosphatase (United States Biochemical Corp.) and incubating for 30 min at 37°C. A 0.5-µl aliquot of the SBE reaction and 0.5 µl of SBE size standard 80 (Beckman Coulter) were mixed and analyzed by capillary electrophoresis on a CEQ8000 genetic analysis system (Beckman Coulter). As a quality control measure, 10% of the DNA samples were random selected and independently analyzed to confirm the original genotyping calls.
Plasmids and Cloning of CYP2C9 Promoter Region. A CYP2C9 genomic fragment containing positions -1515 to +184 was amplified from a single Hispanic DNA sample that, based on sequence analysis, was homozygous for the CYP2C9*1 allele. The amplification product was digested with the restriction enzymes HindIII (CYP2C9 position -1454) and SacI (CYP2C9 position 66) and subsequently cloned into these same sites within the multiple cloning site of the pBluescript II KS+ vector (Stratagene). Site-directed mutagenesis was performed to convert the adenosine residue at CYP2C9 position -1 to cytosine, thereby creating an NcoI site at the translation start codon. After mutagenesis, the complete CYP2C9 insert was sequenced using M13 forward, M13 reverse, and internal primers to verify the identity of the clone and to validate its fidelity relative to the reference sequence. The mutagenized amplicon was excised from the pBluescript II KS+ vector using the restriction enzymes NcoI and HindIII, and cloned into these same sites within the pGL3Basic vector, generating pRNH905.
Because of persistent amplification difficulties, an alternate approach was used to isolate a CYP2C9 genomic fragment containing positions -5910 to -1325. A human bacterial artificial chromosome clone containing the region of interest (CTD2343I24, chromosome 10, position 96,586,284 to position 96,697,248) was purchased from Open Biosystem's Clone Resource. The bacterial artificial chromosome clone DNA was isolated using a midi prep kit designed for large constructs (QIAGEN) and digested using EciI. The 16,589-bp fragment, representing CYP2C9 position -7870 to +8719, was gel purified and digested with HindIII and XmnI and subsequently cloned into the HindIII and SmaI sites of pRNH905, generating pRNH923.
A similar PCR approach as that described above for the proximal promoter fragment was used to amplify a CYP2C9 genomic fragment from position -10,108 to -5517. The resulting amplicon was treated with Taq polymerase (Invitrogen) to add a single deoxyadenosine (A) to the 3' end of the PCR product and subsequently cloned into pCR2.1 (Invitrogen). The relative orientation of the insert was checked by digestion with Bst1107I and SacI. A clone yielding product sizes of 4509 and 3990 bp was selected, and the 4509-bp fragment was subsequently cloned into these same sites in pRNH923, generating pRNH924. Thus, this final construct contained CYP2C9 sequences from position -10,008 to -1 (chromosomal coordinates 96678320 to 96688429) directing the expression of the luciferase gene.
Three CYP2C9 fragments were excised from pRNH924 and cloned into the pBluescript II KS+ vector to perform site-directed mutagenesis and to introduce the various haplotype sequences into the CYP2C9/luciferease reporter vector. pRNH971 was constructed by digesting pRNH924 with HindIII and SalI, and the resulting 3376-bp fragment containing the luciferase cassette and CYP2C9 position -1453 to -1 was cloned into the HindIII and SalI sites of the pBluescript II KS+ vector. pRNH972 was constructed by digesting pRNH924 with HindIII and XbaI, and the resulting 4379-bp fragment was cloned into these same sites of pBluescript II KS+. The third mutagenesis plasmid, containing CYP2C9 sequences from position -10,008 to -5833, was constructed by digesting pRNH924 with XbaI and SacI, and then cloning the resulting 4321-bp fragment into these same sites in pBluescript II KS+, generating pRNH973. After DNA sequence analysis of each plasmid, multiple rounds of site-directed mutagenesis were performed to introduce all combinations of SNPs present in each of the nine inferred CYP2C9 variant haplotypes as well as in the reference. The individual fragments were then reassembled within the pGL3Basic backbone, generating pRNH954 (haplotype 1), pRNH955 (haplotype 1A), pRNH974 (haplotype 1B), pRNH956 (haplotype 1C), pRNH957 (haplotype 1D), pRNH958 (haplotype 1E), pRNH959 (haplotype 2A), pRNH960 (haplotype 2B), pRNH962 (haplotype 3A), and pRNH961 (haplotype 3B).
Cell Culture and Transfection of HepG2 Cells. HepG2 human hepatoma cells were cultured in Eagle's minimal essential media supplemented with 10% fetal bovine serum, 50 U/ml penicillin, and 50 µg/ml streptomycin, and maintained in a humidified incubator at 37°C in an atmosphere of 5% CO2. Cultures were never allowed to grow beyond approximately 80% confluence. All experiments were performed with cells between passage 4 and 15. For transfection studies, 1.5 x 105 cells were seeded into 24-well culture dishes and 24 h later were transfected with 1.4 µg of Lipofectamine 2000, 0.5 µg of test luciferase reporter plasmid, and 50 ng each of pCMVβgal, pCMVHNF-4
, pBJ5HNF-1
, and pBJ5DCoH in Opti-MEM reduced serum medium. After incubation for 24 h at 37°C, transfection medium was replaced with normal growth medium, and cells were incubated for an additional 24 h. For induction experiments, transfection medium was replaced with Opti-MEM reduced serum medium supplemented with 10 µM rifampicin or vehicle (0.1% DMSO). The rifampicin concentration and the time of harvest post-treatment were based on results from initial optimization experiments in which rifampicin concentrations of 5, 10, and 25 µM were added for 10, 16, and 24 h (data not shown). Cell lysates were prepared and luciferase assays were performed according to the manufacturer's instructions. Data were normalized with respect to β-galactosidase activity to correct for transfection efficiency and expressed as relative luciferase activity. Data are reported as the mean ± S.D. of at least three determinations performed with two independently prepared luciferase reporter plasmids.
Sequence and Data Analysis. DNA sequence data were analyzed using DNAStar software (LaserGene, Madison, WI); minimal acceptable quality scores values were set at
12. Sequences were scanned for potential transcription factor recognition sequences using the Match Program and the TRANSFAC Professional version 11.3 database (BIOBASE Corporation, Wolfenbuettel, Germany). Both the Liver-Specific and Vertebrate NonRedundant Matrix Profiles were used with search criteria that minimized the identification of false positives. All SNPs identified were tested for deviations from Hardy-Weinberg equilibrium with the use of a
2 test. The frequencies of individual SNPs were compared using Fisher's exact test (GraphPad InStat version 3.05; GraphPad Software Inc., San Diego, CA). Based upon the observed frequencies of each SNP, haplotype analysis was inferred using PHASE version 2.1, with all parameters set at default values except that 5000 iterations were performed with a thinning value of 1 and a burn in of 1000 (Stephens and Donnelly, 2003
). Functional differences among the different haplotypes were assessed using transient expression assays and compared by one-way ANOVA with a Holm-Sidak post hoc test (SigmaStat version 3.11; Systat Software, Inc., Point Richmond, CA). Inferred haplotype frequencies between non-Latino whites and Hispanics of Mexican descent were compared using a Student's t test (SigmaStat version 3.11). In all instances, an
value of 0.05 was accepted as significant.
| Results |
|---|
|
|
|---|
|
CYP2C9 SNP Validation in a Hispanic Population of Mexican Descent. The PDR24 is an anonymous and blinded DNA panel from individuals representative of the United States population, but it provides no insight into the presence or frequency of specific genetic variants within defined ethnic or racial groups. Thus, further characterization of the identified SNPs was needed in the Hispanic population. To define CYP2C9 promoter region haplotypes, allelic frequencies of the 31 SNPs identified during the discovery phase, the g.-1537G>A SNP previously shown to be linked to the CYP2C9*3 allele in non-Latino whites (Veenstra et al., 2005
), as well as the SNPs defining the CYP2C9*2, CYP2C9*3, CYP2C9*4, CYP2C9*5, and CYP2C9*6 alleles were determined in 193 Hispanic Americans of Mexican descent by either sequence analysis (for SNP clusters) or multiplex SBE. All SNPs identified were in Hardy-Weinberg equilibrium except g.-1188T>C (
2 = 5.37) and g.-4302C>T(
2 = 6.19), both genotyped using multiplex SBE. To eliminate possible assay error for these two positions, all DNA samples with variant alleles at g.-1188T>C or g.-4302C>T were resequenced, and in all instances, the original genotype call was confirmed. Five of the nine novel SNPs identified in the 24PDR were observed in the Hispanic population (g.-8422A>G, g.-5146G>C, g.-5143A>C, g.5140A>T, and g.-4302C>T) (Table 2). The g.-4302C>T SNP was observed at a frequency of 10.9%, whereas all other novel SNPs were observed at frequencies <0.5%. The most common promoter variant in the Hispanic population was g.-1188T>C at an observed frequency of 20.7%, whereas the g.-8422A>G, g.-8416T>G, and g.-7336G>A SNPs were the least commonly observed variants at frequencies <1%. Two previously reported CYP2C9 structural variants, g.3608C>T(CYP2C9*2) and g.42614A>C(CYP2C9*3) (Rettie et al., 1994
; Haining et al., 1996
), also were observed in the Hispanic study population at 7.0 and 4.4%, respectively (Table 3). Other previously reported CYP2C9 exon variants, g.42615T>C (CYP2C9*4), g.42619C>G (CYP2C9*5), and g.10601delA (CYP2C9*6) (Imai et al., 2000
; Dickmann et al., 2001
; Kidd et al., 2001
), were not observed.
|
|
The CYP2C9 promoter SNP frequencies observed in the Hispanic population were compared with the previously reported frequencies in a non-Latino white cohort (Veenstra et al., 2005
) (Table 2). Five upstream genetic variants (g.-8422A>G, g.-5146G>C, g-5143A>C, g.-5140A>T, and g.-4302C>T) observed in the Hispanic population were absent in the non-Latino white population, one of which, g.-4302C>T, was not found because of the lack of ascertainment of that region in the later population (see GenBank accession no. AY702706
[GenBank]
; "region not scanned" because of repetitive sequences). In total, 21 SNPs were observed in both ethnic groups, five of which occurred at different frequencies. For example, the g.-620G>T, g.-1096A>G, g.-1188T>C, g.-2663delTG, and g.-3089G>A SNPs were observed in the Hispanic population at frequencies of 4.9, 7.0, 20.7, 7.3, and 10.9%, respectively, compared with 11.0, 12.0, 35.0, 17.0, and 17.0%, respectively, in the non-Latino white population (P < 0.05).
Using the determined minor allelic frequencies of upstream and structural variants, 34 CYP2C9 haplotypes in total were inferred using PHASE version 2.1, 10 of which occurred at a frequency greater than 1%. This latter group was compared with the previously reported haplotypes in a non-Latino white study group (Table 4) (Veenstra et al., 2005
). It is noteworthy that for the purposes of comparison, the SNPs reported previously in the non-Latino white cohort were reanalyzed, but only including the promoter and *2 and *3 SNPs. Haplotype 1 was designated as the reference sequence, and it was inferred to occur at a frequency not significantly different from that reported in the non-Latino white population (Table 4). Four CYP2C9 haplotypes were deduced in the Hispanic population, but not in the non-Latino white population. Haplotype 1A (g.-8553C>A) was present at 3.2%, 1B (g.-4302C>T) at 10.0%, 1E (g.-1188T>C) at 1.3%, and 2A (g.-3597A>G, g.-3360T>C, g.-1188T>C, g.-1096A>G, g.-485T>A, g.-484C>A, and g.3608C>T) at 2.1%. Five variant CYP2C9 haplotypes inferred in both the Hispanic and non-Latino white populations were present at different frequencies. For example, haplotype 1C (g.-3089G>A and g.-1188T>C) was deduced in the Hispanic population at a higher frequency than that observed in the non-Latino white population (i.e., 3.1 compared with 0.27%, respectively). In contrast, haplotype 1D was inferred to occur more frequently in the non-Latino white population, 17.2%, compared with the Hispanic population, 6.6%.
|
In addition to the five haplotypes (1A, 1B, 1C, 1D, and 1E) made up solely of SNPs within the CYP2C9 regulatory region, four inferred haplotypes (2A, 2B, 3A, and 3B) consisted of upstream SNPs in linkage disequilibrium with the previously described g.3608C>T (CYP2C9*2) and g.42614A>C (CYP2C9*3) structural allelic variants (Table 4). Of these four haplotype variants, haplotype 2A (g.-3579A>G, g.-3360T>C, g.-1188T>C, g.-1096A>G, g.-485T>A, g.-484C>A, and g.3608C>T) seemed to be unique to the Hispanic population. Haplotypes 2B, 3A, and 3B were present at different frequencies in the two population groups (Table 4).
Effect of Variant Promoter Haplotypes on Constitutive and Rifampicin-Induced CYP2C9 Promoter Activity. To explore the possible functional significance of the identified CYP2C9 variant haplotypes on constitutive promoter activity, we used site-directed mutagenesis to introduce nine common CYP2C9 haplotype variants into a reporter construct containing 10,008 bp of CYP2C9 upstream information driving the expression of the luciferase reporter gene. Negligible luciferase activity was observed when the CYP2C9/luciferase construct was transiently expressed alone (data not shown). Given the important roles HNF-1
, and HNF-4
(Chen et al., 2005
; Kawashima et al., 2006
) have in regulating CYP2C9 promoter activity, the initial experiments were repeated coexpressing both of these factors along with the CYP2C9/luciferase reporter construct. A substantial increase in promoter activity was observed in the presence of both HNF-1
and HNF-4
, but either factor alone had a minimal effect. As such, all subsequent transient expression experiments included expression vectors for both of these factors in the protocol. Of the nine variant haplotypes, 1B, 3A, and 3B exhibited 2.5-, 3.2-, and 2.6-fold decreased CYP2C9 promoter activity compared with the reference construct (Fig. 1). Comparing the SNPs constituting each of these hypomorphic promoter variants against the variants located within putative transcription factor binding sites (Table 1), only the g.-3360T>C SNP is located immediately 5' to a putative NF1/CTF core element (position -1). The T>C transition reduces the matrix match score from 0.747 to 0.722, and of the eight binding sites selected for defining this matrix, five hadaTat position -1 and none had a C.
|
|
| Discussion |
|---|
|
|
|---|
The observed frequency of two previously reported CYP2C9 structural variants, g.3608C>T (CYP2C9*2) and g.42614A>C(CYP2C9*3) were not significantly different between the Hispanics and non-Latino whites, nor were the determined frequencies different from those reported by LLerena et al. (2004
) for Hispanic Americans of Mexican descent. In contrast, there was a significant difference compared with the frequencies reported by Xie et al. (2002
). However, the ancestral background of the Hispanic population was not defined in the latter study. Thus, the discrepancy in these data is probably explained by a population of mixed ancestry.
Based upon the observed minor allele frequencies of 32 upstream and five structural variants, 10 common (>1%) CYP2C9 haplotypes were inferred. The observed haplotype 1 (reference) frequency in Hispanics was 60.53%, compared with 63.64% in non-Latino whites, consistent with a relatively simple haplotype structure. Four CYP2C9 haplotypes (1A, 1B, 1E, and 2A) were inferred in Hispanics, but not in non-Latino whites (Table 2).
The possible contribution of upstream CYP2C9 genetic variability to interindividual pharmacokinetic differences is more controversial. Variant effects on promoter activity in vitro have been observed, but they failed to correlate with differences in in vivo activity (Shintani et al., 2001
). Other studies have failed to show an independent effect of CYP2C9 upstream haplotype variants on mean warfarin clearance or dose (King et al., 2004
; Takahashi et al., 2004
; Veenstra et al., 2005
). However, the design of these studies precluded the ability to test the impact of the CYP2C9 upstream polymorphisms independently of the variants affecting CYP2C9 catalytic activity. In the current report, haplotype 1B, consisting of only the g.-4302C>T variant, exhibited a significant decrease in promoter activity in vitro. Combined with its relatively high frequency (i.e., 10.0%), these data would be consistent with reduced constitutive CYP2C9 expression in this population. It is interesting that the frequency of the g.-4302C>T SNP was not in Hardy-Weinberg equilibrium, suggesting a possible selective pressure for the presence of the variant allele, although this deviation may also be due to recent population admixture. Significant decreases in in vitro constitutive promoter activity also were observed with haplotypes 3A and 3B. Thus, the results of our study suggest that, in addition to the defective enzymatic function of CYP2C9.3, a decrease in basal CYP2C9 transcription also may contribute to the overall observed CYP2C9*3 phenotype.
In vitro studies using primary hepatocytes have shown that CYP2C9 mRNA, protein, and catalytic activity are all increased by drugs such as rifampicin, hyperforin, and phenobarbital through a PXR-dependent mechanism (Chen et al., 2004
). Rifampicin treatment also has been reported to enhance the clearance of CYP2C9 substrates, indicative of CYP2C9 induction in vivo (Williamson et al., 1998
; Niemi et al., 2001
). Cotransfection studies in HepG2 cells with a human PXR expression plasmid and CYP2C9 reporter constructs were used to investigate the possible functional effects of CYP2C9 haplotype variants on rifampicin-dependent enhancement of CYP2C9 promoter activity. With haplotype 1 (reference construct), addition of human PXR alone had no effect on basal CYP2C9 promoter activity, consistent with results reported by Ferguson et al. (2002
). These results, however, conflict with those of Chen et al. (2004
), who reported a 3.8-fold increase in CYP2C9 promoter activity upon addition of human PXR alone. When the cotransfected cells were treated with 10 µM rifampicin, a 3.2-fold increase in promoter activity was observed, similar to the approximate 3-fold increase reported by Chen et al. (2004
). Yet, Ferguson et al. (2002
) failed to see an effect with rifampicin. The discrepancies in these data may result from differences in the amount of CYP2C9 upstream sequences present in the reporter construct or from differences in transfection conditions. When cotransfection studies were performed with the CYP2C9 promoter variants, two of the inferred CYP2C9 haplotypes exhibited an altered induction profile; haplotype 1D resulted in a 1.6-fold reduction in inducibility, whereas haplotype 3B resulted in a 1.4-fold increase in inducibility. Given the frequency of these two haplotypes (6.6 and 2.4%, respectively) and the magnitude of this observed effect, it is possible these variants contribute to observed interindividual differences in CYP2C9 phenotype.
Previous studies have identified two CAR-responsive elements within the CYP2C9 promoter at positions -2898 and -1839 that bind human CAR in vitro and transactivate reporter constructs (Ferguson et al., 2002
). Furthermore, Gerbal-Chaloin et al. (2001
) described a 4-fold increase in CYP2C9 mRNA after treatment with phenobarbital, a known human CAR ligand. CAR-mediated phenobarbital induction of CYP2C9 promoter activity was not observed in the current study, although a 2-fold PXR-mediated phenobarbital enhancement of CYP2C9 promoter activity was observed. This latter observation supports the recent findings of Chen et al. (2004
) who determined that human PXR is responsible for the induction of human CYP2C9 by both rifampicin and phenobarbital.
To gain some insight into which SNP, or combination of SNPs, might be responsible for the observed altered activities, a careful comparison between all observed haplotypes was made. Reduced constitutive activity was observed with haplotypes 1B, 3A, and 3B. As indicated earlier, the g.-4302C>T SNP is unique to haplotype 1B and as such, it is assumed to be causative, yet it does not alter any known or putative transcription factor binding sites. Several SNPs are found in common between haplotypes 3A and 3B, but not in other haplotypes (g.-8897C>A, g.-7419A>G, g.-5813A>G, g.-5661C>A, g.-4877G>A, -1911T>C, -1885C>G, g.-1537G>A, and g.-981G>A). However, similar to the haplotype 1B g.-4302C>T SNP, none of these variants alter known or putative regulatory elements. Altered PXR-mediated inducibility was observed with both haplotypes 1D and 3B. In haplotype 1D (g.-3089G>A, g.-2663delTG, and g.-1188T>C), two other inferred haplotypes (1C and 1E) also contain g.-3089G>A, g.-1188T>C, or both, but they exhibit promoter activities no different from the reference construct. The remaining SNP in haplotype 1D, g.-2663delTG, is not present in any other inferred haplotype and as such, it may be causative for the reduced inducibility of this promoter variant. Haplotype 3B resulted in an increase in rifampicin-induced CYP2C9 promoter activity, whereas haplotype 3A exhibited no difference compared with the reference control. Yet, the only difference between these two haplotypes is the presence of the g.-1188T>C SNP in haplotype 3B. This observation suggests that the g.-1188T>C transition might be responsible for the observed difference in induction profiles. However, two other haplotypes (1C and 1E) also contain g.-1188T>C; yet, they exhibit no difference relative to the control construct. Furthermore, a recent study by Sandberg et al. (2004
) found that the g.-1188T>C variant did not affect gene expression in vitro. Finally, haplotype 1D also contains the g.-1188T>C SNP, and, in contrast to haplotype 3B, it was associated with decreased induction of CYP2C9 promoter activity. These findings suggest that g.-1188T>C cannot solely be responsible for the observed increase in activity, but rather it contributes to a combinatorial effect. The simplest explanation as to how individual or clusters of SNPs might function to alter gene regulation is by modifying or eliminating transcription factor binding. However, this does not seem to be the case for the variants identified herein, and it is clearly not the case for many regulatory polymorphisms. A survey of 247 known promoters and 647 haplotype variants by Buckland et al. (2005
) revealed that only 35% of the functional regulatory variants identified were localized within predicted transcription factor binding sites. Thus, many regulatory polymorphisms may affect gene expression through other yet to be identified sequence-specific mechanisms.
In summary, this is the first extensive study of CYP2C9 haplotype and in vitro functional analysis in a Hispanic population of Mexican descent. Interindividual and interethnic differences were observed in the incidence and frequency of CYP2C9 regulatory polymorphisms. Furthermore, based on in vitro assays, several of the inferred haplotypes are predicted to significantly decrease and/or alter PXR-mediated rifampicin-dependent CYP2C9 induction. These observations suggest that genetic variation within CYP2C9 regulatory sequences probably contributes to differences in CYP2C9 phenotype both within and among different populations. However, in vivo studies will be required to determine the relevance of these haplotype variants, if any, to clinical outcomes.
| Acknowledgements |
|---|
| Footnotes |
|---|
ABBREVIATIONS: P450, cytochrome P450; SNP, single-nucleotide polymorphism; bp, base pair(s); 24PDR, 24 sample set from the Coriell Polymorphism Discovery Resource; DMSO, dimethyl sulfoxide; NF1/CTF, nuclear factor 1/CCAAT transcription factor site; HNF, hepatic nuclear factor; DCoH, dimerization cofactor; PCR, polymerase chain reaction; kbp, kilobase pair(s); SBE, single-base extension; PXR, pregnane X receptor; CAR, constitutive androstane receptor; ANOVA, analysis of variance.
The online version of this article (available at http://molpharm.aspetjournals.org) contains supplemental material. ![]()
1 Current affiliation: Department of Experimental and Clinical Pharmacology, College of Pharmacy, University of Minnesota, Minneapolis, Minnesota. ![]()
Address correspondence to: Dr. Ronald N. Hines, Medical College of Wisconsin, TBRC/CRI/CPPT, 8701 Watertown Plank Rd., Milwaukee WI 53226. E-mail: rhines{at}mcw.edu
| References |
|---|
|
|
|---|
Blaisdell J, Jorge-Nebert LF, Coulter S, Ferguson SS, Lee SJ, Chanas B, Xi T, Mohrenweiser H, Ghanayem B, and Goldstein JA (2004) Discovery of new potentially defective alleles of human CYP2C9. Pharmacogenetics 14: 527-537.[CrossRef][Medline]
Buckland PR, Hoogendoorn B, Coleman SL, Guy CA, Smith SK, and O'Donovan MC (2005) Strong bias in the location of functional promoter polymorphisms. Hum Mutat 26: 214-223.[CrossRef][Medline]
Chen Y, Ferguson SS, Negishi M, and Goldstein JA (2004) Induction of human CYP2C9 by rifampicin, hyperforin, and phenobarbital is mediated by the pregnane X receptor. J Pharmacol Exp Ther 308: 495-501.
Chen Y, Kissling G, Negishi M, and Goldstein JA (2005) The nuclear receptors constitutive androstane receptor and pregnane X receptor cross-talk with hepatic nuclear factor 4
to synergistically activate the human CYP2C9 promoter. J Pharmacol Exp Ther 314: 1125-1133.
Davies NM, McLachlan AJ, Day RO, and Williams KM (2000) Clinical pharmacokinetics and pharmacodynamics of celecoxib: a selective cyclo-oxygenase-2 inhibitor. Clin Pharmacokinet 38: 225-242.[CrossRef][Medline]
Dickmann LJ, Rettie AE, Kneller MB, Kim RB, Wood AJJ, Stein CM, Wilkinson GR, and Schwarz UI (2001) Identification and functional characterization of a new CYP2C9 variant (CYP2C9*5) expressed among African Americans. Mol Pharmacol 60: 382-387.
Ferguson SS, LeCluyse EL, Negishi M, and Goldstein JA (2002) Regulation of human CYP2C9 by the constitutive androstane receptor: discovery of a new distal binding site. Mol Pharmacol 62: 737-746.
Gage BF, Eby C, Milligan PE, Banet GA, Duncan JR, and McLeod HL (2004) Use of pharmacogenetics and clinical factors to predict the maintenance dose of warfarin. Thromb Haemost 91: 87-94.[Medline]
Gerbal-Chaloin S, Pascussi J-M, Pichard-Garcia L, Daujat M, Waechter F, Fabre J-M, Carrère N, and Maurel P (2001) Induction of CYP2C genes in human hepatocytes in primary culture. Drug Metab Dispos 29: 242-251.
Haining RL, Hunter AP, Veronese ME, Trager WF, and Rettie AE (1996) Allelic variants of human cytochrome P450 2C9: baculovirus-mediated expression, purification, structural characterization, substrate stereoselectivity, and prochiral selectivity of the wild-type and I359L mutant forms. Arch Biochem Biophys 333: 447-458.[CrossRef][Medline]
Ibeanu GC and Goldstein JA (1995) Transcriptional regulation of human CYP2C genes: functional comparison of CYP2C9 and CYP2C18 promoter regions. Bio-chemistry 34: 8028-8036.[CrossRef][Medline]
Imai J, Ieiri I, Mamiya K, Miyahara S, Furuumi H, Nanba E, Yamane M, Fukumaki Y, Ninomiya H, Tashiro N, et al. (2000) Polymorphism of the cytrochrome P450 (CYP) 2C9 gene in Japanese epileptic patients: genetic analysis of the CYP2C9 locus. Pharmacogenetics 10: 85-89.[CrossRef][Medline]
Kamali F, Khan TI, King BP, Frearson R, Kesteven P, Wood P, Daly AK, and Wynne H (2004) Contribution of age, body size, and CYP2C9 genotype to anticoagulant response to warfarin. Clin Pharmacol Ther 75: 204-212.[CrossRef][Medline]
Kawashima S, Kobayashi K, Takama K, Higuchi T, Furihata T, Hosokawa M, and Chiba K (2006) Involvement of hepatocyte nuclear factor 4
in the different expression level between CYP2C9 and CYP2C19 in the human liver. Drug Metab Dispos 34: 1012-1018.
Kidd RS, Curry TB, Gallagher S, Edeki T, Blaisdell J, and Goldstein JA (2001) Identification of a null allele of CYP2C9 in an African-American exhibiting toxicity to phenytoin. Pharmacogenetics 11: 803-808.[CrossRef][Medline]
King BP, Khan TI, Aithal GP, Kamali F, and Daly AK (2004) Upstream and coding region CYP2C9 polymorphisms: correlation with warfarin dose and metabolism. Pharmacogenetics 14: 813-822.[CrossRef][Medline]
Klose TS, Ibeanu GC, Ghanayem BI, Pedersen LG, Li L, Hall SD, and Goldstein JA (1998) Identification of residues 286 and 289 as critical for conferring substrate specificity of human CYP2C9 for diclofenac and ibuprofen. Arch Biochem Biophys 357: 240-248.[CrossRef][Medline]
Kuo CJ, Conley PB, Hsieh CL, Francke U, and Crabtree GR (1990) Molecular cloning, functional expression, and chromosomal localization of mouse hepatocyte nuclear factor 1. Proc Natl Acad Sci U S A 87: 9839-9842.
Lee CR, Goldstein JA, and Pieper JA (2002) Cytochrome P450 2C9 polymorphisms: a comprehensive review of the in vitro and human data. Pharmacogenetics 12: 251-263.[CrossRef][Medline]
Lehmann JM, McKee DD, Watson MA, Wilson TM, Moore JT, and Kliewer SA (1998) The human orphan nuclear receptor PXR is activated by compounds that regulate CYP3A4 gene expression and cause drug interactions. J Clin Invest 102: 1016-1023.[Medline]
Lindblad-Toh K, Winchester E, Daly MJ, Wang DG, Hirschhorn JN, Laviolette J-P, Ardlie K, Reich DE, Robinson E, Sklar P, et al. (2000) Large-scale discovery and genotyping of single-nucleotide polymorphisms in the mouse. Nat Genet 24: 381-386.[CrossRef][Medline]
LLerena A, Dorado P, O'Kirwan F, Jepson R, Licinio J, and Wong ML (2004) Lower frequency of CYP2C9*2 in Mexican-Americans compared to Spaniards. Pharmacogenomics J 4: 403-406.[CrossRef][Medline]
Martinez-Jiménez CP, Gomez-Lechon MJ, Castell JV, and Jover R (2005) Transcriptional regulation of the human hepatic CYP3A4: identification of a new distal enhancer region responsive to CCAAT/enhancer-binding protein β isoforms (liver activating protein and liver inhibitory protein). Mol Pharmacol 67: 2088-2101.
Mendel DB, Khavari PA, Conley PB, Graves MK, Hansen LP, Admon A, and Crabtree GR (1991) Characterization of a cofactor that regulates dimerization of a mammalian homeodomain protein. Science 254: 1762-1767.
Miners JO and Birkett DJ (1998) Cytochrome P4502C9: an enzyme of major importance in human drug metabolism. Br J Clin Pharmacol 45: 525-538.[CrossRef][Medline]
Niemi M, Backman JT, Neuvonen M, Neuvonen PJ, and Kivisto KT (2001) Effects of rifampin on the pharmacokinetics and pharmacodynamics of glyburide and glipizide. Clin Pharmacol Ther 69: 400-406.[CrossRef][Medline]
Phillips KA, Veenstra DL, Oren E, Lee JK, and Sadee W (2001) Potential role of pharmacogenomics in reducing adverse drug reactions: a systematic review. JAMA 286: 2270-2279.
Rettie AE, Wienkers LC, Gonzalez FJ, Trager WF, and Korzekwa KR (1994) Impaired (S)-warfarin metabolism catalyzed by the R144C allelic variant of CYP2C9. Pharmacogenetics 4: 39-42.[Medline]
Sandberg M, Johansson I, Christensen M, Rane A, and Eliasson E (2004) The impact of CYP2C9 genetics and oral contraceptives on cytochrome P450 2C9 phenotype. Drug Metab Dispos 32: 484-489.
Schwarz UI (2003) Clinical relevance of genetic polymorphisms in the human CYP2C9 Gene. Eur J Clin Invest 33(Suppl 2): 23-30.[CrossRef][Medline]
Shimada T, Yamazaki H, Mimura M, Inui Y, and Guengerich FP (1994) Interindividual variations in human liver cytochrome P-450 enzymes involved in the oxidation of drugs, carcinogens and toxic chemicals: studies with liver microsomes of 30 Japanese and 30 Caucasians. J Pharmacol Exp Ther 270: 414-423.
Shintani M, Ieiri I, Inoue K, Mamiya K, Ninomiya H, Tashiro N, Higuchi S, and Otsubo K (2001) Genetic polymorphisms and functional characterization of the 5'-flanking region of the human CYP2C9 gene: in vitro and in vivo studies. Clin Pharmacol Ther 70: 175-182.[CrossRef][Medline]
Stephens M and Donnelly P (2003) A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet 73: 1162-1169.[CrossRef][Medline]
Stoffel M and Duncan SA (1997) The maturity-onset diabetes of the young (MODY1) transcription factor HNF4a regulates expression of genes required for glucose transport and metabolism. Proc Natl Acad Sci U S A 94: 13209-13214.
Takahashi H, Ieiri I, Wilkinson GR, Mayo G, Kashima T, Kimura S, Otsubo K, and Echizen H (2004) 5'-Flanking region polymorphisms of CYP2C9 and their relationship to S-warfarin metabolism in white and Japanese patients. Blood 103: 3055-3057.
Veenstra DL, Blough DK, Higashi MK, Farin FM, Srinouanprachan S, Rieder MJ, and Rettie AE (2005) CYP2C9 Haplotype structure in European American warfarin patients and association with clinical outcomes. Clin Pharmacol Ther 77: 353-364.[CrossRef][Medline]
Williams JA, Hyland R, Jones BC, Smith DA, Hurst S, Goosen TC, Peterkin V, Koup JR, and Ball SE (2004) Drug-drug interactions for UDP-glucuronosyltransferase substrates: a pharmacokinetic explanation for typically observed low exposure (AUCi/AUC) ratios. Drug Metab Dispos 32: 1201-1208.
Williamson KM, Patterson JH, McQueen RH, Adams KF Jr, and Pieper JA (1998) Effects of erythromycin or rifampin on losartan pharmacokinetics in healthy volunteers. Clin Pharmacol Ther 63: 316-323.[CrossRef][Medline]
Xie HG, Prasad HC, Kim RB, and Stein CM (2002) CYP2C9 Allelic variants: ethnic distribution and functional significance. Adv Drug Deliv Rev 54: 1257-1270.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||