Compendium of gene expression profiles comprising a baseline model of the human liver drug metabolism transcriptome

Determinants of gene expression in the human liver: Impact of aging and sex on xenobiotic metabolism

There is a need to characterize the potential susceptibility of older adults to toxicity from environmental chemical exposures. Liver xenobiotic metabolizing enzymes (XMEs) play important roles in detoxifying and eliminating xenobiotics. We examined global gene expression in the livers of young (21-45 years) and old (69+ years) men and women. Differentially expressed genes (DEG) were identified using two-way ANOVA (p ≤ 0.05). We identified 1437 and 1670 DEGs between young and old groups in men and women, respectively. Only a minor number of the total number of genes overlapped (146 genes). Aging increased or decreased pathways involved in inflammation and intermediary metabolism, respectively. Aging led to numerous changes in the expression of XME genes or genes known to control their expression (~90 genes). Out of 10 cytochrome P450s activities examined, there were increased activities of CYP1A2 and CYP2C9 enzymes in the old groups. We also identified sex-dependent genes that were more numerous in the young group (1065) than in the old group (202) and included changes in XMEs. These studies indicate that the livers from aging humans when compared to younger adults exhibit changes in XMEs that may lead to differences in the metabolism of xenobiotics.

Deregulation of the Genes that Are Involved in Drug Absorption, Distribution, Metabolism, and Excretion in Hepatocellular Carcinoma

Genes involved in drug absorption, distribution, metabolism, and excretion (ADME) are called ADME genes. Currently, 298 genes that encode phase I and II drug metabolizing enzymes, transporters, and modifiers are designated as ADME genes by the PharmaADME Consortium. ADME genes are highly expressed in the liver and their levels can be influenced by liver diseases such as hepatocellular carcinoma (HCC). In this study, we obtained RNA-sequencing and microRNA (miRNA)-sequencing data from 371 HCC patients via The Cancer Genome Atlas liver hepatocellular carcinoma project and performed ADME gene-targeted differential gene expression analysis and expression correlation analysis. Two hundred thirty-three of the 298 ADME genes (78%) were expressed in HCC. Of these genes, almost one-quarter (58 genes) were significantly downregulated, while only 6% (15) were upregulated in HCC relative to healthy liver. Moreover, one-half (14/28) of the core ADME genes (CYP1A2, CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, CYP2E1, CYP3A4, NAT1, NAT2, UGT2B7, SLC22A1, SLCO1B1, and SLCO1B3) were downregulated. In addition, about one-half of the core ADME genes were positively correlated with each other and were also positively (AHR, ARNT, HNF4A, PXR, CAR, PPARA, and RXRA) or negatively (PPARD and PPARG) correlated with transcription factors known as ADME modifiers. Finally, we show that most miRNAs known to regulate core ADME genes are upregulated in HCC. Collectively, these data reveal 1) an extensive transcription factor-mediated ADME coexpression network in the liver that efficiently coordinates the metabolism and elimination of endogenous and exogenous compounds; and 2) a widespread deregulation of this network in HCC, most likely due to deregulation of both transcriptional and post-transcriptional (miRNA) pathways.

Mathematical Models in the Description of Pregnane X Receptor (PXR)-Regulated Cytochrome P450 Enzyme Induction

The pregnane X receptor (PXR) is a drug/xenobiotic-activated transcription factor of crucial importance for major cytochrome P450 xenobiotic-metabolizing enzymes (CYP) expression and regulation in the liver and the intestine. One of the major target genes regulated by PXR is the cytochrome P450 enzyme (CYP3A4), which is the most important human drug-metabolizing enzyme. In addition, PXR is supposed to be involved both in basal and/or inducible expression of many other CYPs, such as CYP2B6, CYP2C8, 2C9 and 2C19, CYP3A5, CYP3A7, and CYP2A6. Interestingly, the dynamics of PXR-mediated target genes regulation has not been systematically studied and we have only a few mechanistic mathematical and biologically based models describing gene expression dynamics after PXR activation in cellular models. Furthermore, few indirect mathematical PKPD models for prediction of CYP3A metabolic activity in vivo have been built based on compartmental models with respect to drug–drug interactions or hormonal crosstalk. Importantly, several negative feedback loops have been described in PXR regulation. Although current mathematical models propose these adaptive mechanisms, a comprehensive mathematical model based on sufficient experimental data is still missing. In the current review, we summarize and compare these models and address some issues that should be considered for the improvement of PXR-mediated gene regulation modelling as well as for our better understanding of the quantitative and spatial dynamics of CYPs expression.

mRNA transfection retrofits cell-based assays with xenobiotic metabolism

The US EPA's ToxCast program is designed to assess chemical perturbations of molecular and cellular endpoints using a variety of high-throughput screening (HTS) assays. However, existing HTS assays have limited or no xenobiotic metabolism which could lead to false positive (chemical is detoxified in vivo) as well as false negative results (chemical is bioactivated in vivo) and thus potential mischaracterization of chemical hazard. To address this challenge, the ten most prevalent human liver cytochrome P450 (CYP) enzymes were introduced into a human cell line (HEK293T) with low endogenous metabolic capacity. The CYP enzymes were introduced via transfection of modified mRNAs as either singlets or as a mixture in relative proportions as expressed in human liver. Initial experiments using luminogenic substrates demonstrate that CYP enzyme activities are significantly increased when co-transfected with an mRNA encoding a CYP accessory protein, P450 oxidoreductase (POR). Transfected HEK293T cells demonstrate the ability to produce predicted metabolites following treatment with well-studied CYP substrates for at least 18 h post-treatment. As a demonstration of how this method can be used to retrofit existing HTS assays, a proof-of-concept screen for cytotoxicity in HEK293T cells was conducted using 56 test compounds. The results demonstrate that the xenobiotic metabolism conferred by transfection of CYP-encoding mRNAs shifts the dose-response relationship for some of the tested chemicals such as aflatoxin B1 (bioactivation) and fenazaquin (detoxification). Overall, transfection of CYP-encoding mRNAs is an effective and portable solution for retrofitting existing cell-based HTS assays with metabolic competence.

Two persistent organic pollutants which act through different xenosensors (alpha-endosulfan and 2,3,7,8 tetrachlorodibenzo-p-dioxin) interact in a mixture and downregulate multiple genes involved in human hepatocyte lipid and glucose metabolism

Individuals, typically, are exposed to mixtures of environmental xenobiotics affecting multiple organs and acting through different xenosensors and pathways in species and cell-type specific manners. 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and α-endosulfan are Persistent Organic Pollutants (POPs) and endocrine disruptors which act through different xenosensors and accumulate in the liver. Our objective in this HEALS study was to investigate the effects of the mixture of these POPs on gene expression in a human-derived hepatocyte cell line, HepaRG. We found that, in spite of having largely uncorrelated effects, TCDD and α-endosulfan, when mixed, alter the expression of genes. The combined effects of the mixture of the POPs significantly altered the expression of 100 genes (42 up- and 58 down-regulated) whereas the same concentration of either POP alone did not alter significantly the expression of these genes. For 32 other genes, selective inhibitory crosstalk between TCDD and α-endosulfan was observed. One of the POPs inhibited the effect, on gene expression, of the other in the mixture although, when used alone, that POP did not affect expression. The expression of another 82 genes was significantly altered (up- or down-regulated) by a single POP. The addition of the second POP either increased, in the same direction, the effect on gene expression or had no further effect. At low concentrations (0.2 nM TCDD and 1 μM α-endosulfan), the POPs still had significant effects and the levels of expression of the corresponding proteins were found to be affected for some genes. Particularly striking was the 80-90% inhibition, by the mixture, of the expression of a number of genes of several hepatic intermediary metabolic pathways (glycerolipid metabolism, FXR/RXR activation, glycolysis/gluconeogenesis, retinoid and bile acid biosynthesis), whereas each pollutant alone had only a moderate effect.

Human PXR modulates hepatotoxicity associated with rifampicin and isoniazid co-therapy

Rifampicin and isoniazid co–therapy frequently causes liver injury in humans. A pregnane X receptor–humanized mouse model revealed that rifampicin and isoniazid co–treatment causes accumulation of protoporphyrin IX, an endogenous hepatotoxin, in the liver via a pregnane X receptor–mediated alteration of heme biosynthesis pathway. These results provide novel insight into the mechanism of rifampicin and isoniazid–induced liver injury that may be applied to clinical management of the hepatotoxicity associated with tuberculosis chemotherapy.

Sources of variance in baseline gene expression in the rodent liver

The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variation due to individual, environmental, and technical factors. Analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies has yielded useful information on baseline fluctuations in liver gene expression in the rodent. Here, studies which highlight contributions of different factors to gene expression variability in the rodent liver are discussed including a large meta-analysis of rat liver, which identified genes that vary in control animals in the absence of chemical treatment. Genes and their pathways that are the most and least variable were identified in a number of these studies. Life stage, fasting, sex, diet, circadian rhythm and liver lobe source can profoundly influence gene expression in the liver. Recognition of biological and technical factors that contribute to variability of background gene expression can help the investigator in the design of an experiment that maximizes sensitivity and reduces the influence of confounders that may lead to misinterpretation of genomic changes. The factors that contribute to variability in liver gene expression in rodents are likely analogous to those contributing to human interindividual variability in drug response and chemical toxicity. Identification of batteries of genes that are altered in a variety of background conditions could be used to predict responses to drugs and chemicals in appropriate models of the human liver.

Systematic genetic and genomic analysis of cytochrome P450 enzyme activities in human liver

Liver cytochrome P450s (P450s) play critical roles in drug metabolism, toxicology, and metabolic processes. Despite rapid progress in the understanding of these enzymes, a systematic investigation of the full spectrum of functionality of individual P450s, the interrelationship or networks connecting them, and the genetic control of each gene/enzyme is lacking. To this end, we genotyped, expression-profiled, and measured P450 activities of 466 human liver samples and applied a systems biology approach via the integration of genetics, gene expression, and enzyme activity measurements. We found that most P450s were positively correlated among themselves and were highly correlated with known regulators as well as thousands of other genes enriched for pathways relevant to the metabolism of drugs, fatty acids, amino acids, and steroids. Genome-wide association analyses between genetic polymorphisms and P450 expression or enzyme activities revealed sets of SNPs associated with P450 traits, and suggested the existence of both cis-regulation of P450 expression (especially for CYP2D6) and more complex trans-regulation of P450 activity. Several novel SNPs associated with CYP2D6 expression and enzyme activity were validated in an independent human cohort. By constructing a weighted coexpression network and a Bayesian regulatory network, we defined the human liver transcriptional network structure, uncovered subnetworks representative of the P450 regulatory system, and identified novel candidate regulatory genes, namely, EHHADH, SLC10A1, and AKR1D1. The P450 subnetworks were then validated using gene signatures responsive to ligands of known P450 regulators in mouse and rat. This systematic survey provides a comprehensive view of the functionality, genetic control, and interactions of P450s.

Development of a quantitative model of pregnane X receptor (PXR) mediated xenobiotic metabolizing enzyme induction

The pregnane X receptor plays an integral role in the regulation of hepatic metabolism. It has been shown to regulate CYP3A4, which is the most abundant cytochrome P450 in the human liver. With its large and flexible ligand-binding domain, PXR can be activated by an enormous range of relatively small, hydrophobic, exogenous compounds. Upon activation, PXR partners with the retinoid X receptor (RXR) to form a heterodimer. The newly formed heterodimer binds to an appropriate DNA response element, causing increased transcription. This leads to an induction in the level of CYP3A4. These mechanistic steps are included into a biologically-based mathematical model. The quantitative model predicts fold level inductions of CYP3A4 mRNA and protein in response to PXR activation. Model parameter values have been taken from literature when appropriate. Unknown parameter values are estimated by optimizing the model results to published in vivo and in vitro data sets. A sensitivity analysis is performed to evaluate the model structure and identify future data needs which would be critical to revising the model.

Hepatic gene expression in morbidly obese women: implications for disease susceptibility

The objective of this study was to determine the molecular bases of disordered hepatic function and disease susceptibility in obesity. We compared global gene expression in liver biopsies from morbidly obese (MO) women undergoing gastric bypass (GBP) surgery with that of women undergoing ventral hernia repair who had experienced massive weight loss (MWL) following prior GBP. Metabolic and hormonal profiles were examined in MO vs. MWL groups. Additionally, we analyzed individual profiles of hepatic gene expression in liver biopsy specimens obtained from MO and MWL subjects. All patients underwent preoperative metabolic profiling. RNAs were extracted from wedge biopsies of livers from MO and MWL subjects, and analysis of mRNA expression was carried out using Affymetrix HG-U133A microarray gene chips. Genes exhibiting greater than twofold differential expression between MO and MWL subjects were organized according to gene ontology and hierarchical clustering, and expression of key genes exhibiting differential regulation was quantified by real-time-polymerase chain reaction (RT-PCR). We discovered 154 genes to be differentially expressed in livers of MWL and MO subjects. A total of 28 candidate disease susceptibility genes were identified that encoded proteins regulating lipid and energy homeostasis (PLIN, ENO3, ELOVL2, APOF, LEPR, IGFBP1, DDIT4), signal transduction (MAP2K6, SOCS-2), postinflammatory tissue repair (HLA-DQB1, SPP1, P4HA1, LUM), bile acid transport (SULT2A, ABCB11), and metabolism of xenobiotics (GSTT2, CYP1A1). Using gene expression profiling, we have identified novel candidate disease susceptibility genes whose expression is altered in livers of MO subjects. The significance of altered expression of these genes to obesity-related disease is discussed.

Activation of CAR and PXR by Dietary, Environmental and Occupational Chemicals Alters Drug Metabolism, Intermediary Metabolism, and Cell Proliferation

The constitutive androstane receptor (CAR) and the pregnane × receptor (PXR) are activated by a variety of endogenous and exogenous ligands, such as steroid hormones, bile acids, pharmaceuticals, and environmental, dietary, and occupational chemicals. In turn, they induce phase I-III detoxification enzymes and transporters that help eliminate these chemicals. Because many of the chemicals that activate CAR and PXR are environmentally-relevant (dietary and anthropogenic), studies need to address whether these chemicals or mixtures of these chemicals may increase the susceptibility to adverse drug interactions. In addition, CAR and PXR are involved in hepatic proliferation, intermediary metabolism, and protection from cholestasis. Therefore, activation of CAR and PXR may have a wide variety of implications for personalized medicine through physiological effects on metabolism and cell proliferation; some beneficial and others adverse. Identifying the chemicals that activate these promiscuous nuclear receptors and understanding how these chemicals may act in concert will help us predict adverse drug reactions (ADRs), predict cholestasis and steatosis, and regulate intermediary metabolism. This review summarizes the available data on CAR and PXR, including the environmental chemicals that activate these receptors, the genes they control, and the physiological processes that are perturbed or depend on CAR and PXR action. This knowledge contributes to a foundation that will be necessary to discern interindividual differences in the downstream biological pathways regulated by these key nuclear receptors.

Toxicogenomic biomarkers for liver toxicity

Toxicogenomics (TGx) is a widely used technique in the preclinical stage of drug development to investigate the molecular mechanisms of toxicity. A number of candidate TGx biomarkers have now been identified and are utilized for both assessing and predicting toxicities. Further accumulation of novel TGx biomarkers will lead to more efficient, appropriate and cost effective drug risk assessment, reinforcing the paradigm of the conventional toxicology system with a more profound understanding of the molecular mechanisms of drug-induced toxicity. In this paper, we overview some practical strategies as well as obstacles for identifying and utilizing TGx biomarkers based on microarray analysis. Since clinical hepatotoxicity is one of the major causes of drug development attrition, the liver has been the best documented target organ for TGx studies to date, and we therefore focused on information from liver TGx studies. In this review, we summarize the current resources in the literature in regard to TGx studies of the liver, from which toxicologists could extract potential TGx biomarker gene sets for better hepatotoxicity risk assessment.

Sexually dimorphic regulation and induction of P450s by the constitutive androstane receptor (CAR)

The constitutive androstane receptor (CAR) is a xenosensing nuclear receptor and regulator of cytochrome P450s (CYPs). However, the role of CAR as a basal regulator of CYP expression nor its role in sexually dimorphic responses have been thoroughly studied. We investigated basal regulation and sexually dimorphic regulation and induction by the potent CAR activator TCPOBOP and the moderate CAR activator Nonylphenol (NP). NP is an environmental estrogen and one of the most commonly found environmental toxicants in Europe and the United States. Previous studies have demonstrated that NP induces several CYPs in a sexually dimorphic manner, however the role of CAR in regulating NP-mediated sexually dimorphic P450 expression and induction has not been elucidated. Therefore, wild-type and CAR-null male and female mice were treated with honey as a carrier, NP, or TCPOBOP and CYP expression monitored by QPCR and Western blotting. CAR basally regulates the expression of Cyp2c29, Cyp2b13, and potentially Cyp2b10 as demonstrated by QPCR. Furthermore, we observed a shift in the testosterone 6alpha/15alpha-hydroxylase ratio in untreated CAR-null female mice to the male pattern, which indicates an alteration in androgen status and suggests a role for androgens as CAR inverse agonists. Xenobiotic-treatments with NP and TCPOBOP induced Cyp2b10, Cyp2c29, and Cyp3a11 in a CAR-mediated fashion; however NP only induced these CYPs in females and TCPOBOP induced these CYPs in both males and females. Interestingly, Cyp2a4, was only induced in wild-type male mice by TCPOBOP suggesting Cyp2a4 induction is not sensitive to CAR-mediated induction in females. Overall, TCPOBOP and NP show similar CYP induction profiles in females, but widely different profiles in males potentially related to lower sensitivity of males to either indirect or moderate CAR activators such as NP. In summary, CAR regulates the basal and chemically inducible expression of several sexually dimorphic xenobiotic metabolizing P450s in a manner that varies depending on the ligand.

Sources of variation in baseline gene expression levels from toxicogenomics study control animals across multiple laboratories

Background

The use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variance due to individual, environmental, and technical factors. Meta-analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies could yield useful information on baseline fluctuations in gene expression, although control animal data has not been available on a scale and in a form best served for data-mining.

Results

A dataset of control animal microarray expression data was assembled by a working group of the Health and Environmental Sciences Institute's Technical Committee on the Application of Genomics in Mechanism Based Risk Assessment in order to provide a public resource for assessments of variability in baseline gene expression. Data from over 500 Affymetrix microarrays from control rat liver and kidney were collected from 16 different institutions. Thirty-five biological and technical factors were obtained for each animal, describing a wide range of study characteristics, and a subset were evaluated in detail for their contribution to total variability using multivariate statistical and graphical techniques.

Conclusion

The study factors that emerged as key sources of variability included gender, organ section, strain, and fasting state. These and other study factors were identified as key descriptors that should be included in the minimal information about a toxicogenomics study needed for interpretation of results by an independent source. Genes that are the most and least variable, gender-selective, or altered by fasting were also identified and functionally categorized. Better characterization of gene expression variability in control animals will aid in the design of toxicogenomics studies and in the interpretation of their results.

Expression profiling of interindividual variability following xenobiotic exposures in primary human hepatocyte cultures

To examine the magnitude of human variability across the entire transcriptome after chemical challenge, we profiled gene expression responses to three different prototypic chemical inducers in primary human hepatocyte cultures from ten independent donors. Correlation between basal expression in any two hepatocyte donors ranged from r(2) values of 0.967 to 0.857, and chemical treatment tended to negatively impact correlation between donors. Including anticipated target genes, 10,812, 8373, and 7847 genes were changed in at least one donor by Aroclor 1254 (A1254), di(2-ethylhexyl) phthalate (DEHP), and phenobarbital (PB), respectively. A subset of these gene targets (n=41) were altered with a high level of reproducibility in at least 9 donors, gene responses that correlated well with literature-reported mechanism of action. Filtering responses to the level of gene subsets clarified the biological impact associated with the respective chemical effectors, in lieu of substantial interindividual variation among donor responses. In these respects, the use of hierarchical clustering methods successfully grouped seven of the ten donors into chemical-specific rather than donor-specific clusters. However, at the whole-genome level, the magnitude of conserved gene expression changes among donors was surprisingly small, with fewer than 50% of the gene responses altered by a single chemical conserved in more than one donor. The use of higher level descriptors, such as those defined by the PANTHER classification system, may enable more consistent categorization of gene expression changes across individuals, as increased reproducibility was identified using this method.

The liver pharmacological and xenobiotic gene response repertoire

We have used a supervised classification approach to systematically mine a large microarray database derived from livers of compound-treated rats. Thirty-four distinct signatures (classifiers) for pharmacological and toxicological end points can be identified. Just 200 genes are sufficient to classify these end points. Signatures were enriched in xenobiotic and immune response genes and contain un-annotated genes, indicating that not all key genes in the liver xenobiotic responses have been characterized. Many signatures with equal classification capabilities but with no gene in common can be derived for the same phenotypic end point. The analysis of the union of all genes present in these signatures can reveal the underlying biology of that end point as illustrated here using liver fibrosis signatures. Our approach using the whole genome and a diverse set of compounds allows a comprehensive view of most pharmacological and toxicological questions and is applicable to other situations such as disease and development.

Expression of constitutive androstane receptor, hepatic nuclear factor 4 alpha, and P450 oxidoreductase genes determines interindividual variability in basal expression and activity of a broad scope of xenobiotic metabolism genes in the human liver

Identification of genetic variation predictive of clearance rate of a wide variety of prescription drugs could lead to cost-effective personalized medicine. Here we identify regulatory genes whose variable expression level among individuals may have widespread effects upon clearance rate of a variety of drugs. Twenty liver samples with variable CYP3A activity were profiled for expression level and activity of xenobiotic metabolism genes as well as genes involved in the regulation thereof. Regulatory genes whose expression level accounted for the highest degree of collinearity among expression levels of xenobiotic metabolism genes were identified as possible master regulators of drug clearance rate. Significant linear correlations (p < 0.05) were identified among mRNA levels of CYP2A6, CYP2B6, CYP2C8, CYP2C9, CYP2C19, MRP2, OATP2, P450 oxidoreductase (POR), and UDP-glucuronosyltranferase 1A1, suggesting that these xenobiotic metabolism genes are coregulated at the transcriptional level. Using partial regression analysis, constitutive androstane receptor (CAR) and hepatic nuclear factor 4 alpha (HNF4 alpha) were identified as the nuclear receptors whose expression levels are most strongly associated with expression of coregulated xenobiotic metabolism genes. POR expression level, which is also associated with CAR and HNF4 alpha expression level, was found to be strongly associated with the activity of many cytochromes P450. Thus, interindividual variation in the expression level of CAR, HNF4 alpha, and POR probably determines variation in expression and activity of a broad scope of xenobiotic metabolism genes and, accordingly, clearance rate of a variety of xenobiotics. Identification of polymorphisms in these candidate master regulator genes that account for their variable expression among individuals may yield readily detectable biomarkers that could serve as predictors of xenobiotic clearance rate.

Microarray-based compendium of hepatic gene expression profiles for prototypical ADME gene-inducing compounds in rats and mice in vivo

To examine species-specific aspects of the induction of absorption, distribution, metabolism and excretion (ADME)-related genes, we used 25 000 gene oligonucleotide microarrays to construct a rodent gene-response compendium that compared hepatic gene expression profiles and developed consensus aryl hydrocarbon receptor (AhR), constitutive androstane receptor (CAR) and pregnane X-receptor (PXR) ligand signatures relevant to drug clearance. Twenty-six inducer compounds were chosen from the literature. Rats and mice received one of six dose levels (log2 dose escalation, 32-fold dose range) of each compound daily for 3 days. Animals were necropsied 6-9 h after the last dose, and tissues were collected for RNA analysis. Hepatic gene expression profiles were obtained using Rosetta Resolver expression analysis system, and ADME-related genes were extracted. Cross-talk among nuclear receptors or hepatoxicity at high dose levels resulted in large signatures (usually >1000 genes at p < 0.01) for most compounds. After ADME gene transcript enrichment, agglomerative clustering separated AhR ligands from CAR/PXR ligands, but it was difficult to distinguish CAR from PXR ligands. Consensus signatures were derived from groups of AhR, CAR and PXR ligands; and cross-talk among responding genes was determined. Many compounds had distinct log dose-response profiles, and relative potencies for ligands were established. Robust responses by CYP1A1, CYP2B10 (CAR responsive in mice) and CYP2B15 (CAR responsive in rats) and CYP3A1 (PXR responsive in rats) were used to benchmark the relative potency of different ligands and to determine the relative selectivity for AhR, CAR or PXR. By using a compendium of gene expression profiles, we defined species-specific induction patterns across the ADME transcriptome.