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

Aromatase inhibitor-induced modulation of breast density: clinical and genetic effects

Background: Change in breast density may predict outcome of women receiving adjuvant hormone therapy for breast cancer. We performed a prospective clinical trial to evaluate the impact of inherited variants in genes involved in oestrogen metabolism and signalling on change in mammographic percent density (MPD) with aromatase inhibitor (AI) therapy. Methods: Postmenopausal women with breast cancer who were initiating adjuvant AI therapy were enrolled onto a multicentre, randomised clinical trial of exemestane vs letrozole, designed to identify associations between AI-induced change in MPD and single-nucleotide polymorphisms in candidate genes. Subjects underwent unilateral craniocaudal mammography before and following 24 months of treatment. Results: Of the 503 enrolled subjects, 259 had both paired mammograms at baseline and following 24 months of treatment and evaluable DNA. We observed a statistically significant decrease in mean MPD from 17.1 to 15.1% (P<0.001), more pronounced in women with baseline MPD ⩾20%. No AI-specific difference in change in MPD was identified. No significant associations between change in MPD and inherited genetic variants were observed. Conclusion: Subjects with higher baseline MPD had a greater average decrease in MPD with AI therapy. There does not appear to be a substantial effect of inherited variants in biologically selected candidate genes

Lack of association between oestrogen receptor polymorphisms and change in bone mineral density with tamoxifen therapy.

BACKGROUND: Tamoxifen, a selective oestrogen receptor (ER) modulator, increases bone mineral density (BMD) in postmenopausal women and decreases BMD in premenopausal women. We hypothesised that inherited variants in candidate genes involved in oestrogen signalling and tamoxifen metabolism might be associated with tamoxifen effects in bone. METHODS: A total of 297 women who were initiating tamoxifen therapy were enrolled in a prospective multicentre clinical trial. Lumbar spine and total hip BMD values were measured using dual-energy X-ray absorptiometry (DXA) at baseline and after 12 months of tamoxifen therapy. Single-nucleotide polymorphisms (SNPs) in ESR1, ESR2, and CYP2D6 were tested for associations in the context of menopausal status and previous chemotherapy, with a mean percentage change in BMD over 12 months. RESULTS: The percentage increase in BMD was greater in postmenopausal women and in those patients who had been treated with chemotherapy. No significant associations between tested SNPs and either baseline BMD or change in BMD with 1 year of tamoxifen therapy were detected. CONCLUSION: The evaluated SNPs in ESR and CYP2D6 do not seem to influence BMD in tamoxifen-treated subjects

Geographic clustering of elevated blood heavy metal levels in pregnant women

Abstract Background Cadmium (Cd), lead (Pb), mercury (Hg), and arsenic (As) exposure is ubiquitous and has been associated with higher risk of growth restriction and cardiometabolic and neurodevelopmental disorders. However, cost-efficient strategies to identify at-risk populations and potential sources of exposure to inform mitigation efforts are limited. The objective of this study was to describe the spatial distribution and identify factors associated with Cd, Pb, Hg, and As concentrations in peripheral blood of pregnant women. Methods Heavy metals were measured in whole peripheral blood of 310 pregnant women obtained at gestational age ~12 weeks. Prenatal residential addresses were geocoded and geospatial analysis (Getis-Ord Gi* statistics) was used to determine if elevated blood concentrations were geographically clustered. Logistic regression models were used to identify factors associated with elevated blood metal levels and cluster membership. Results Geospatial clusters for Cd and Pb were identified with high confidence (p-value for Gi* statistic <0.01). The Cd and Pb clusters comprised 10.5 and 9.2 % of Durham County residents, respectively. Medians and interquartile ranges of blood concentrations (μg/dL) for all participants were Cd 0.02 (0.01–0.04), Hg 0.03 (0.01–0.07), Pb 0.34 (0.16–0.83), and As 0.04 (0.04–0.05). In the Cd cluster, medians and interquartile ranges of blood concentrations (μg/dL) were Cd 0.06 (0.02–0.16), Hg 0.02 (0.00–0.05), Pb 0.54 (0.23–1.23), and As 0.05 (0.04–0.05). In the Pb cluster, medians and interquartile ranges of blood concentrations (μg/dL) were Cd 0.03 (0.02–0.15), Hg 0.01 (0.01–0.05), Pb 0.39 (0.24–0.74), and As 0.04 (0.04–0.05). Co-exposure with Pb and Cd was also clustered, the p-values for the Gi* statistic for Pb and Cd was <0.01. Cluster membership was associated with lower education levels and higher pre-pregnancy BMI. Conclusions Our data support that elevated blood concentrations of Cd and Pb are spatially clustered in this urban environment compared to the surrounding areas. Spatial analysis of metals concentrations in peripheral blood or urine obtained routinely during prenatal care can be useful in surveillance of heavy metal exposure