Introduction: STAT5A and STAT5B are important transcription factors that play a key role in regulation of several important physiological processes including proliferation, survival, mediation of responses to cytokines and in regulating gender differences in drug response genes such as the hepatic cytochrome P450s (CYPs) that are responsible for a large majority of drug metabolism reactions in the human body. STAT5A and STAT5b have a high degree of sequence homology and have been reported to have largely similar functions. Recent studies have, however, indicated that they can also often have distinct and unique roles in regulating gene expression.Objective: In this study, we evaluated the association of STAT5A and STAT5B mRNA expression levels with those of several key hepatic cytochrome P450s (CYPs) and hepatic transcription factors (TFs) and evaluated the potential roles of STAT5A and 5b in mediating gender differences in these CYPs and TFs.Methods: Expression profiling for major hepatic CYP isoforms and transcription factors was performed using RNA sequencing (RNA-seq) in 102 human liver samples (57 female, 45 male). Real time PCR gene expression data for selected CYP5 and TFs was available on a subset of 50 human liver samples (25 female, 25 male) and was used to validate the RNA-seq findings.Results: While STAT5A demonstrated significant negative correlation with expression levels of multiple hepatic transcription factors (including NR1I2 and HNF4A) and DMEs such as CYP3A4 and CYP2C19, STAT5B expression was observed to demonstrate positive associations with several CYP5 and TFs analyzed. As STAT5A and STAT5B have been shown to be important in regulation of gender differences in CYPs, we also analyzed STAT5A and 5b associations with CYPs and TFs separately in males and females and observed gender dependent differential associations of STATs with several CYPs and TFs. Results from the real time PCR validation largely supported our RNA-seq findings.Conclusions: Using both RNA sequencing and real time PCR, we examined the association of STAT5A and STAT5B mRNA expression with CYP and TF gene expression. While STAT5A demonstrated significant negative correlations with expression levels of multiple hepatic TFs (including NR1I2 and HNF4 alpha) and CYPs (eg. CYP3A4, CYP2C19), STAT5B expression was observed to demonstrate positive association with most of the CYPs/TFs analyzed suggesting that STAT5A and STAT5b have potentially different and distinct roles in regulating expression of hepatic drug response genes. Further studies are needed to elucidate the potential roles of STAT5A and 5b in regulation of CYPs/TFs and the potential implications of these findings. (C) 2016 Published by Elsevier Inc.
microRNA (miRNA) mediated regulation of gene expression has emerged as a significant mechanism contributing to variation in gene expression. In this study, we evaluated the potential role of miRNAs in regulating expression of hepatic cytochromes P450 and their transcriptional regulatory genes. We screened the Targetscan database for high scoring miRNA binding site predictions in selected hepatic DMEs and transcription factors. Expression profiling for candidate miRNAs (n=22) and their target genes (n=20) was performed in 50 human liver samples (25 female, 25 male). Significant negative correlations were observed between expression levels of several CYPs/hepatic transcription factors and the hepatic miRNAs studied. Interestingly, hepatic miR-34a demonstrated significant negative correlation with expression levels of multiple hepatic transcription factors (including NR1I2 and HNF4α) and DMEs (CYP3A4, CYP2C19). miR-34a expression was also significantly higher in males than in females in congruence with previous observations of higher CYP3A4 expression in females versus males. A mediation analysis revealed that miR-34a was involved in significant mediation of the association observed between CYP2C19 and several hepatic transcription factors (HNF4α, NR1I2). miR-34a may thus play a key regulatory role and be a key contributory factor to the inter-individual variability observed in expression of key drug metabolizing genes in humans.
The main objective of this study was to investigate the feasibility of using PharmGKB, a pharmacogenomic database, as a source of training data in combination with text of MEDLINE abstracts for a text mining approach to identification of potential gene targets for pathway-driven pharmacogenomics research. We used the manually curated relations between drugs and genes in PharmGKB database to train a support vector machine predictive model and applied this model prospectively to MEDLINE abstracts. The gene targets suggested by this approach were subsequently manually reviewed. Our quantitative analysis showed that a support vector machine classifiers trained on MEDLINE abstracts with single words (unigrams) used as features and PharmGKB relations used for supervision, achieve an overall sensitivity of 85% and specificity of 69%. The subsequent qualitative analysis showed that gene targets "suggested" by the automatic classifier were not anticipated by expert reviewers but were subsequently found to be relevant to the three drugs that were investigated: carbamazepine, lamivudine and zidovudine. Our results show that this approach is not only feasible but may also find new gene targets not identifiable by other methods thus making it a valuable tool for pathway-driven pharmacogenomics research.
The CYP3A locus encodes hepatic enzymes that metabolize many clinically used drugs. However, there is marked interindividual variability in enzyme expression and clearance of drugs metabolized by these enzymes. We utilized comparative genomics and computational prediction of transcriptional factor binding sites to evaluate regions within CYP3A that were most likely to contribute to this variation. We then used a haplotype tagging single-nucleotide polymorphisms (htSNPs) approach to evaluate the entire locus with the fewest number of maximally informative SNPs. We investigated the association between these htSNPs and in vivo CYP3A enzyme activity using a single-point IV midazolam clearance assay. We found associations between the midazolam phenotype and age, diagnosis of hypertension and one htSNP (141689) located upstream of CYP3A4 . 141689 lies near the xenobiotic responsive enhancer module (XREM) regulatory region of CYP3A4 . Cell-based studies show increased transcriptional activation with the minor allele at 141689, in agreement with the in vivo association study findings. This study marks the first systematic evaluation of coding and noncoding variation that may contribute to CYP3A phenotypic variability.
The purpose of this study was to evaluate the affinity of docetaxel for 14 transporter proteins and assess the functional significance of 17 variants in five genes involved in drug elimination. among the transfected models investigated, OATP1B3 (SLCO1B3) was identified as the most efficient influx transporter for docetaxel. none of the observed genotypes (SLCO1B3, ABCB1, and ABCC2) was related with docetaxel clearance in 92 white patients (P > 0.17). however, the simultaneous presence of the CYP3A4*1B and CYP3A5*1A alleles was associated with a 64% increase in docetaxel clearance (P = 0.0015), independent of both sex and CYP3A activity (as determined using the erythromycin breath test). This haplotype was also associated with increased midazolam clearance in another population (P = 0.0198). an analysis of the CYP3A locus among CEPH-HapMap samples revealed that CYP3A4*1B is present exclusively among a subset of CYP3A5 expressors. Therefore, future studies should first stratify the population on the basis of CYP3A5 genotype and then compare CYP3A activity between individuals with and without the CYP3A4*1B allele.
Objectives. Variant cytochrome P450 (CYP) 3A4 alleles cannot explain human variation in CYP3A4 expression. This study investigated whether common single-nucleotide polymorphisms (SNPs) in multidrug resistance 1 (MDRI), encoding P-glycoprotein, or the pregnane X receptor (PXR) were associated with basal or inducible CYP3A4 expression.Methods: MDRI G2677T and C3435T SNPs and a PXR 6-base pair (bp) deletion were genotyped in the deoxyribonucleic acid from 144 human livers in 3 cohorts each phenotyped for basal or rifampin (INN, rifampicin)-inducible hepatic CYP3A4 expression (or both) and in 57 human small bowel biopsy specimens from 3 cohorts each phenotyped for either basal or rifampin-induced CYP3A4 expression (or both).Results: Hepatic CYP3A4 expression/function was significantly higher in persons homozygous for the MDRI 2677T (Ser893) allele compared with persons homozygous for 2677G (Ala893) in all 3 hepatic cohorts. For example, homozygous MDR1 2677 TT livers had higher midazolam hydroxylase activity than homozygous 2677 GG livers (1831 +/- 1336 pmol (.) min(-1) (.) mg protein(-1) versus 1060 +/- 552 pmol (.) min(-1.) mg protein(-1), P =.03). In 2 of the 3 groups the association was observed in men but not in women. For example, homozygous MDR1 2677 TT male hepatocytcs had significantly higher testosterone 6 beta-hydroxylase activity compared with homozygous 2677 GG livers (0.120 +/- 0.06 pmol (.) min(-1) (.) mg protein(-1) versus 0.069 +/- 0.04 pmol (.) min(-1) (.) mg protein(-1), P =.0002). Conversely, rifampin induction of testosterone 6 beta-hydroxylation in primary human hepatocytes was significantly higher in persons homozygous for 2677G (12.0 +/- 5.7-fold) compared with MDR1 homozygous TT carriers (7.3 +/- 4.6-fold) (P =.01). Suggestive evidence for higher CYP3A4 expression in MDR1 2677T carriers was also observed in human intestines. CYP3A4 expression was also related to a 6-bp deletion in PXR in 2 of the liver cohorts. Two-factor ANOVA analysis revealed a significant interaction between the MDRI 2677 SNP and the PXR 6-bp deletion influencing CYP3A4 expression (P =.007).Conclusions: Individuals homozygous for the MDRI 2677T allele show enhanced constitutive CYP3A4 expression in the liver and intestine, as compared with those homozygous for the MDRI 2677G allele,
Previous studies have found that, compared with Whites, Hispanic donor livers had elevated expression of CYP2 enzymes, gene products regulated by the constitutive androstane receptor (CAR). The objectives of the current study were to determine (1) the CAR activation signature in human liver (2) whether other drug detoxification (absorption, distribution, metabolism and excretion (ADME)) genes were differentially expressed in Hispanic versus White livers, and (3) the extent of overlap in the CAR and Hispanic liver transcriptomes. The CAR transcriptome (ADME genes differentially expressed following phenobarbital versus vehicle treatment of human hepatocytes) and the Hispanic liver transcriptome (ADME genes differentially expressed in Hispanic versus White livers) were identified using Affymetrix oligonucleotide arrays. Quantitative real-time polymerase chain reaction (PCR) was used to verify candidate genes in a larger sample size. Comparison of the CAR and Hispanic liver ADME transcriptomes revealed a significant association between the gene changes. Sixty-four per cent of the ADME genes induced more than twofold by phenobarbital were also induced in Hispanics, and 14% of the ADME genes repressed more than twofold by phenobarbital were repressed in Hispanics. In conclusion, compared with Whites, Hispanic donor livers have increased expression of many genes that are transcriptionally regulated by CAR. This result has practical implications to the drug treatment of Hispanic patients.
Objective: Saquinavir, a widely prescribed human immunodeficiency virus I protease inhibitor, has a low and variable oral bioavailability that has been attributed to extensive first-pass extraction mediated by hepatic or intestinal cytochrome P450 (CYP) 3A4 and intestinal P-glycoprotein (P-gp). The polymorphic CYP3A5 has also been shown to influence the saquinavir metabolite/parent urinary ratio, suggesting a role for CYP3A5.Methods: Twenty healthy subjects received a single oral dose of saquinavir (600 mg) with water (control) and, on a separate occasion, with Seville orange juice (a selective intestinal CYP3A4/5 inhibitor). Hepatic CYP3A4 activity was evaluated by use of the erythromycin breath test. Duodenal biopsy specimens were used to assess relative intestinal CYP3A4 and CYP3A5 protein contents. Relative P-gp content was also assessed in the biopsy specimens and in lymphocytes. Genetic polymorphisms in MDR1 (in exon 21 and 26), CYP3A5 (*1 and *3), and CYP3A4*1B were identified by direct sequencing. Saquinavir plasma concentrations were measured by tandem liquid chromatography-mass spectrometry. Pharmacokinetic parameter estimates (maximum concentration, time to reach maximum concentration, area under the concentration-time curve, apparent oral clearance [CL/F]) were computed by standard noncompartmental methods. Stepwise multiple regression analysis was used to identify the hepatic or intestinal variables that predicted variation in saquinavir pharmacokinetic measures.Results: Baseline saquinavir CL/F was not correlated with liver CYP3A4 activity (the erythromycin breath test result), intestinal CYP3A4 content, or intestinal P-gp content (r(2) = 0.08, 0.08, and 0.007, respectively; P > .2). MDR1 genotype and lymphocyte P-gp content were also not predictive. Among the 6 subjects expressing intestinal CYP3A5, the mean saquinavir CL/F was almost twice as high as for the 14 nonexpressors (36.7 L/h [95% confidence interval (CI), 18.7-54.6 L/h] and 19.3 L/h [95% CI, 11.2-27.4 L/h], respectively; P = .03). However, among the 6 CYP3A5 expressors, there was an unexpected negative correlation between CL/F and intestinal CYP3A5 content (r(2) = 0.58, P = .05). Seville orange juice decreased the mean CL/F in all 20 subjects from 24.5 L/h (95% CI, 16.7-32.3 L/h) to 14.7 L/h (95% CI, 8.4-20.6 L/h) (P = .05). The effect size did not appear to be influenced by CYP3A5 expression.Conclusions. The CYP3A5*1 genotype is associated with increased saquinavir CL/F. This does not appear to reflect intestinal CYP3A5 expression and presumably reflects the contribution of hepatic CYP3A5. The interaction with Seville orange juice in subjects not expressing CYP3A5 supports a role for intestinal CYP3A4. However, the modest nature of the interaction, combined with the inability to detect a correlation between CL/F and CYP3A4 enterocyte content, supports our recent in vitro work suggesting a smaller contribution of intestinal CYP3A4 than has been assumed.
The defense mechanisms responsible for protecting the body from endogenous toxins are also involved in the metabolism of drugs and are composed of phase I and phase II drug metabolizing enzymes, as well as drug transporters. Numerous drugs and chemicals have been shown to modulate the expression of the genes involved in these three drug-detoxifying processes. Induction of these genes contributes to both auto-induction of drug clearance and to drug-drug interactions in combination therapies. The orphan nuclear receptors PXR (pregnane X receptor) and CAR (Constitutive androstane receptor) are xenosensors that mediate drug-induced changes by increasing transcription of genes that are involved in drug clearance and disposition. Co-administration of drugs, one of which is a nuclear receptor agonist or antagonist, can either lead to altered clearance of the second drug and severe toxicity, or a loss of therapeutic efficacy or an imbalance in physiological substrate concentrations, providing a novel molecular mechanism for drug-drug interactions. Thus, genetic variability in these nuclear receptors will contribute to human variation in the magnitude of clinically significant drug-drug interactions. This review describes common PXR and CAR genetic variants that have been identified to date in the human population and the functional consequence of these variant alleles. In addition, alternatively spliced variants of PXR and CAR that may also contribute to individual variability as well as tissue specific expression of these receptors are also described. Identification of PXR and CAR genetic variants and alternatively spliced mRNAs may ultimately allow predictions of interindividual differences in target gene induction and drug-drug interactions.
The constitutive androstane receptor (CAR) NR1I3 is a transcription factor that upon activation by xenobiotics induces transcription of drug-metabolizing and drug transporter genes. Our goal was to identify whether alternative splicing of CAR makes an important contribution to the generation of novel CAR proteins. The wild-type CAR mRNA (CAR.1) and splice variants (SVs) were amplified from human liver cDNAs and in a panel of cDNAs from 36 human tissues, using exon 1- and 3'-untranslated region primers, cloned and sequenced. Twenty-two unique hCAR splice variants (CAR-SVs) containing different combinations of splicing (deletion of exons 2, 4, 5, 7, partial deletion of exon 9, or insertion of 12 or 15 base pairs from introns 6 or 7) were identified. CAR mRNAs were expressed in small intestine, kidney, testis, adrenal, and brain caudate nucleus. Intestine expressed only CAR.1 mRNA, whereas spleen, heart, and prostate expressed only CAR-SVs. In vitro transcription and translation of CAR-SVs lacking exon 2 (missing ATG start site) generated CAR proteins that differed in M(r) from CAR.1. These CAR-SVs used a translation initiation site in exon 1, generating CAR with a unique amino-terminal sequence. Transcripts lacking part of exon 9 altered the CAR reading frame generating CAR proteins with a unique carboxy-terminal end. CAR-SVs demonstrated compromised binding to CYP2B6 NR elements and transcriptional activation of a CYP2B6 luciferase reporter. The expression of CAR in additional human cell types increases the range of tissue specific transcriptional responses regulated by this receptor, suggesting additional biological roles for CAR and CAR-SV proteins in these tissues.
CYP2B6 metabolizes many drugs, and its expression varies greatly. CYP2B6 genotype-phenotype associations were determined using human livers that were biochemically phenotyped for CYP2B6 ( mRNA, protein, and CYP2B6 activity), and genotyped for CYP2B6 coding and 5'-flanking regions. CYP2B6 expression differed significantly between sexes. Females had higher amounts of CYP2B6 mRNA ( 3.9-fold, P < 0.001), protein ( 1.7-fold, P < 0.009), and activity (1.6-fold, P < 0.05) than did male subjects. Furthermore, 7.1% of females and 20% of males were poor CYP2B6 metabolizers. Striking differences among different ethnic groups were observed: CYP2B6 activity was 3.6- and 5.0-fold higher in Hispanic females than in Caucasian ( P < 0.022) or African-American females ( P < 0.038). Ten single nucleotide polymorphisms ( SNPs) in the CYP2B6 promoter and seven in the coding region were found, including a newly identified 13072A > G substitution that resulted in an Lys139Glu change. Many CYP2B6 splice variants (SV) were observed, and the most common variant lacked exons 4 to 6. A nonsynonymous SNP in exon 4 (15631G > T), which disrupted an exonic splicing enhancer, and a SNP 15582C > T in an intron-3 branch site were correlated with this SV. The extent to which CYP2B6 variation was a predictor of CYP2B6 activity varied according to sex and ethnicity. The 1459C > T SNP, which resulted in the Arg487Cys substitution, was associated with the lowest level of CYP2B6 activity in livers of females. The intron-3 15582C > T SNP ( in significant linkage disequilibrium with a SNP in a putative hepatic nuclear factor 4 (HNF4) binding site) was correlated with lower CYP2B6 expression in females. In conclusion, we found several common SNPs that are associated with polymorphic CYP2B6 expression.