A Penalized Mixture Model Approach in Genotype/Phenotype Association Analysis for Quantitative Phenotypes

A mixture normal model has been developed to partition genotypes in predicting quantitative phenotypes. Its estimation and inference are performed through an EM algorithm. This approach can conduct simultaneous genotype clustering and hypothesis testing. It is a valuable method for predicting the distribution of quantitative phenotypes among multi-locus genotypes across genes or within a gene. This mixture model’s performance is evaluated in data analyses for two pharmacogenetics studies. In one example, thirty five CYP2D6 genotypes were partitioned into three groups to predict pharmacokinetics of a breast cancer drug, Tamoxifen, a CYP2D6 substrate (p-value = 0.04). In a second example, seventeen CYP2B6 genotypes were categorized into three clusters to predict CYP2B6 protein expression (p-value = 0.002). The biological validities of both partitions are examined using established function of CYP2D6 and CYP2B6 alleles. In both examples, we observed genotypes clustered in the same group to have high functional similarities. The power and recovery rate of the true partition for the mixture model approach are investigated in statistical simulation studies, where it outperforms another published method

Aromatase inhibitors augment nociceptive behaviors in rats and enhance the excitability of sensory neurons

Although aromatase inhibitors (AIs) are commonly used therapies for breast cancer, their use is limited because they produce arthralgia in a large number of patients. To determine whether AIs produce hypersensitivity in animal models of pain, we examined the effects of the AI, letrozole, on mechanical, thermal, and chemical sensitivity in rats. In ovariectomized (OVX) rats, administering a single dose of 1 or 5mg/kg letrozole significantly reduced mechanical paw withdrawal thresholds, without altering thermal sensitivity. Repeated injection of 5mg/kg letrozole in male rats produced mechanical, but not thermal, hypersensitivity that extinguished when drug dosing was stopped. A single dose of 5mg/kg letrozole or daily dosing of letrozole or exemestane in male rats also augmented flinching behavior induced by intraplantar injection of 1000nmol of adenosine 5'-triphosphate (ATP). To determine whether sensitization of sensory neurons contributed to AI-induced hypersensitivity, we evaluated the excitability of neurons isolated from dorsal root ganglia of male rats chronically treated with letrozole. Both small and medium-diameter sensory neurons isolated from letrozole-treated rats were more excitable, as reflected by increased action potential firing in response to a ramp of depolarizing current, a lower resting membrane potential, and a lower rheobase. However, systemic letrozole treatment did not augment the stimulus-evoked release of the neuropeptide calcitonin gene-related peptide (CGRP) from spinal cord slices, suggesting that the enhanced nociceptive responses were not secondary to an increase in peptide release from sensory endings in the spinal cord. These results provide the first evidence that AIs modulate the excitability of sensory neurons, which may be a primary mechanism for the effect of these drugs to augment pain behaviors in rats

Population Pharmacokinetic Modeling To Estimate the Contributions of Genetic and Nongenetic Factors to Efavirenz Disposition

Efavirenz pharmacokinetics is characterized by large between-subject variability, which determines both therapeutic response and adverse effects. Some of the variability in efavirenz pharmacokinetics has been attributed to genetic variability in cytochrome P450 genes that alter efavirenz metabolism, such as CYP2B6 and CYP2A6. While the effects of additional patient factors have been studied, such as sex, weight, and body mass index, the extent to which they contribute to variability in efavirenz exposure is inconsistently reported. The aim of this analysis was to develop a pharmacometric model to quantify the contribution of genetic and nongenetic factors to efavirenz pharmacokinetics. A population-based pharmacokinetic model was developed using 1,132 plasma efavirenz concentrations obtained from 73 HIV-seronegative volunteers administered a single oral dose of 600 mg efavirenz. A two-compartment structural model with absorption occurring by zero- and first-order processes described the data. Allometric scaling adequately described the relationship between fat-free mass and apparent oral clearance, as well as fat mass and apparent peripheral volume of distribution. Inclusion of fat-free mass and fat mass in the model mechanistically accounted for correlation between these disposition parameters and sex, weight, and body mass index. Apparent oral clearance of efavirenz was reduced by 25% and 51% in subjects predicted to have intermediate and slow CYP2B6 metabolizer status, respectively. The final pharmacokinetic model accounting for fat-free mass, fat mass, and CYP2B6 metabolizer status was consistent with known mechanisms of efavirenz disposition, efavirenz physiochemical properties, and pharmacokinetic theory. (This study has been registered at ClinicalTrials.gov under identifier NCT00668395.

Rifampin modulation of xeno- and endobiotic conjugating enzyme mRNA expression and associated microRNAs in human hepatocytes

Rifampin is a pleiotropic inducer of multiple drug metabolizing enzymes and transporters. This work utilized a global approach to evaluate rifampin effects on conjugating enzyme gene expression with relevance to human xeno- and endo-biotic metabolism. Primary human hepatocytes from 7 subjects were treated with rifampin (10 μmol/L, 24 hours). Standard methods for RNA-seq library construction, EZBead preparation, and NextGen sequencing were used to measure UDP-glucuronosyl transferase UGT, sulfonyltransferase SULT, N acetyltransferase NAT, and glutathione-S-transferase GST mRNA expression compared to vehicle control (0.01% MeOH). Rifampin-induced (>1.25-fold) mRNA expression of 13 clinically important phase II drug metabolizing genes and repressed (>1.25-fold) the expression of 3 genes (P < .05). Rifampin-induced miRNA expression changes correlated with mRNA changes and miRNAs were identified that may modulate conjugating enzyme expression. NAT2 gene expression was most strongly repressed (1.3-fold) by rifampin while UGT1A4 and UGT1A1 genes were most strongly induced (7.9- and 4.8-fold, respectively). Physiologically based pharmacokinetic modeling (PBPK) was used to simulate the clinical consequences of rifampin induction of CYP3A4- and UGT1A4-mediated midazolam metabolism. Simulations evaluating isolated UGT1A4 induction predicted increased midazolam N-glucuronide exposure (~4-fold) with minimal reductions in parent midazolam exposure (~10%). Simulations accounting for simultaneous induction of both CYP3A4 and UGT1A4 predicted a ~10-fold decrease in parent midazolam exposure with only a ~2-fold decrease in midazolam N-glucuronide metabolite exposure. These data reveal differential effects of rifampin on the human conjugating enzyme transcriptome and potential associations with miRNAs that form the basis for future mechanistic studies to elucidate the interplay of conjugating enzyme regulatory elements

Efavirenz inhibits the human ether-a-go-go related current (hERG) and induces QT interval prolongation in CYP2B6*6*6 allele carriers

Background Efavirenz (EFV) has been associated with torsade de pointes despite marginal QT interval lengthening. Since EFV is metabolized by the cytochrome P450 (CYP) 2B6 enzyme, we hypothesized that EFV would lengthen the rate-corrected QT (QTcF) interval in carriers of the CYP2B6*6 decreased functional allele. Objective The primary objective of this study was to evaluate EFV-associated QT interval changes with regard to CYP2B6 genotype and to explore mechanisms of QT interval lengthening. Methods EFV was administered to healthy volunteers (n=57) as a single 600 mg dose followed by multiple doses to steady-state. Subjects were genotyped for known CYP2B6 alleles and ECGs and EFV plasma concentrations were obtained serially. Whole-cell, voltage-clamp experiments were performed on cells stably expressing hERG and exposed to EFV in the presence and absence of CYP2B6 expression. Results EFV demonstrated a gene-dose effect and exceeded the FDA criteria for QTcF interval prolongation in CYP2B6*6/*6 carriers. The largest mean time-matched differences ΔΔQTcF were observed at 6 hrs (14 ms; 95% CI [1; 27]), 12 hrs (18 ms; 95% CI [−4; 40] and 18 hrs (6 ms; 95% CI [−1; 14]) in the CYP2B6*6/*6 genotype. EFV concentrations exceeding 0.4 µg/mL significantly inhibited outward hERG tail currents (P<0.05). Conclusions This study demonstrates that homozygous carriers of CYP2B6*6 allele may be at increased risk for EFV-induced QTcF interval prolongation via inhibition of hERG

Effects of exemestane and letrozole therapy on plasma concentrations of estrogens in a randomized trial of postmenopausal women with breast cancer

PURPOSE: Inter-individual differences in estrogen concentrations during treatment with aromatase inhibitors (AIs) may contribute to therapeutic response and toxicity. The aim of this study was to determine plasma concentrations of estradiol (E2), estrone (E1), and estrone sulfate (E1S) in a large cohort of AI-treated breast cancer patients. METHODS: In a randomized, multicenter trial of postmenopausal women with early-stage breast cancer starting treatment with letrozole (n = 241) or exemestane (n = 228), plasma estrogen concentrations at baseline and after 3 months were quantitated using a sensitive mass spectrometry-based assay. Concentrations and suppression below the lower limit of quantification (LLOQ) were compared between estrogens and between drugs. RESULTS: The ranges of baseline estrogen concentrations were <LLOQ-361 pg/mL for E2, <LLOQ-190 pg/mL for E1, and 8.3-4060 pg/mL for E1S. For E2, the frequency of suppression below the LLOQ was not statistically significantly different between AIs (exemestane: 89.0%, letrozole: 86.9%, p = 0.51). However, patients on letrozole were more likely to achieve suppression below the LLOQ of both E1 (exemestane: 80.1%, letrozole: 90.1%, p = 0.005) and E1S (exemestane: 17.4%, letrozole: 54.9%, p = 4.34e-15). After 3 months of AI therapy, the ranges of estrogen concentrations were <LLOQ-63.8 pg/mL, <LLOQ-36.7 pg/mL, and <LLOQ-1090 pg/mL for E2, E1, and E1S, respectively. During treatment, 16 patients had an increased concentration compared to the baseline concentration of at least one estrogen. CONCLUSIONS: Letrozole had greater suppression of plasma E1 and E1S than exemestane, though the response was highly variable among patients. Additional research is required to examine the clinical relevance of differential estrogen suppression

Efavirenz inhibits the human ether-a-go-go related current (hERG) and induces QT interval prolongation in CYP2B6*6*6 allele carriers

Background Efavirenz (EFV) has been associated with torsade de pointes despite marginal QT interval lengthening. Since EFV is metabolized by the cytochrome P450 (CYP) 2B6 enzyme, we hypothesized that EFV would lengthen the rate-corrected QT (QTcF) interval in carriers of the CYP2B6*6 decreased functional allele. Objective The primary objective of this study was to evaluate EFV-associated QT interval changes with regard to CYP2B6 genotype and to explore mechanisms of QT interval lengthening. Methods EFV was administered to healthy volunteers (n=57) as a single 600 mg dose followed by multiple doses to steady-state. Subjects were genotyped for known CYP2B6 alleles and ECGs and EFV plasma concentrations were obtained serially. Whole-cell, voltage-clamp experiments were performed on cells stably expressing hERG and exposed to EFV in the presence and absence of CYP2B6 expression. Results EFV demonstrated a gene-dose effect and exceeded the FDA criteria for QTcF interval prolongation in CYP2B6*6/*6 carriers. The largest mean time-matched differences ΔΔQTcF were observed at 6 hrs (14 ms; 95% CI [1; 27]), 12 hrs (18 ms; 95% CI [−4; 40] and 18 hrs (6 ms; 95% CI [−1; 14]) in the CYP2B6*6/*6 genotype. EFV concentrations exceeding 0.4 µg/mL significantly inhibited outward hERG tail currents (P<0.05). Conclusions This study demonstrates that homozygous carriers of CYP2B6*6 allele may be at increased risk for EFV-induced QTcF interval prolongation via inhibition of hERG

Estrogen Receptor Genotypes Influence Hot Flash Prevalence and Composite Score Before and After Tamoxifen Therapy

PURPOSE: Hot flashes are common and frequently lead to drug discontinuation among women prescribed tamoxifen. We determined whether genetic polymorphisms in estrogen receptors (ESRs) alpha and beta (ESR1 and ESR2, respectively) are associated with tamoxifen-induced hot flashes. PATIENTS AND METHODS: We determined ESR1 PvuII and XbaI and ESR2-02 genotypes in 297 women who were initiating tamoxifen. One-week hot flash diaries were collected to calculate a hot flash score (frequency x severity) before and 1, 4, 8, and 12 months after starting tamoxifen. RESULTS: Approximately 80% of 297 participants reported hot flashes before or during the first year of tamoxifen. After 4 months of tamoxifen, premenopausal women who did not receive adjuvant chemotherapy had a four-fold increase in hot flash score (from 5.9 to 23.6; P = .003) compared with a 1.17-fold increase (from 19.6 to 23; P = .34) in those who received chemotherapy. In premenopausal women, increased number of ESR1 PvuII and XbaI CG alleles was associated with higher baseline hot flash scores compared with those who had other haplotypes (P = .0026). At 4 months, postmenopausal women with ESR1 PvuII CC and ESR2-02 GG genotypes had 4.6 times increases in hot flash scores than other postmenopausal women (56 v 12; P = .0007). Women who had the ESR2-02 AA genotype were significantly less likely to experience tamoxifen-induced hot flashes than women who carried at least one ESR-02 G allele (hazard ratio, 0.26; 95% CI, 0.10 to 0.63; P = .001). CONCLUSION: Knowledge of menopausal status, prior chemotherapy, and ESR genotype may help predict which women are most likely to suffer hot flashes during tamoxifen treatment

Composite Functional Genetic and Comedication CYP2D6 Activity Score in Predicting Tamoxifen Drug Exposure Among Breast Cancer Patients

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/97292/1/0091270009359182.pd