Identification of chemical modulators of the constitutive activated receptor (CAR) in a gene expression compendium

Case examples of an evaluation of the human relevance of the pyrethroids/pyrethrins-induced liver tumours in rodents based on the mode of action

Rodent carcinogenicity studies are useful for screening for human carcinogens but they are not perfect. Some modes of action (MOAs) lead to cancers in both experimental rodents and humans, but others that lead to cancers in rodents do not do so in humans. Therefore, analysing the MOAs by which chemicals produce tumours in rodents and determining the relevance of such tumour data for human risk are critical. Recently, experimental data were obtained as case examples of an evaluation of the human relevance of pyrethroid (metofluthrin and momfluorothrin)- and pyrethrins-induced liver tumours in rats based on MOA. The MOA analysis, based on the International Programme on Chemical Safety (IPCS) framework, concluded that experimental data strongly support that the postulated MOA for metofluthrin-, momfluorothrin- and pyrethrins-produced rat hepatocellular tumours is mediated by constitutive androstane receptor (CAR) activation. Since metofluthrin and momfluorothrin are close structural analogues, reproducible outcomes for both chemicals provide confidence in the MOA findings. Furthermore, cultured human hepatocyte studies and humanized chimeric mouse liver studies demonstrated species difference between human hepatocytes (refractory to the mitogenic effects of these compounds) and rat hepatocytes (sensitive to their mitogenic effects). These data strongly support the hypothesis that the CAR-mediated MOA for liver tumorigenesis is of low carcinogenic risk for humans. In this research, in addition to cultured human hepatocyte studies, the usefulness of the humanized chimeric liver mouse models was clearly demonstrated. These data substantially influenced decisions in regulatory toxicology. In this review I comprehensively discuss the human relevance of the CAR-mediated MOA for rodent liver tumorigenesis based on published information, including our recent molecular research on CAR-mediated MOA.

Human relevance of rodent liver tumour formation by constitutive androstane receptor (CAR) activators

A large number of nongenotoxic chemicals have been shown to increase the incidence of liver tumours in rats and/or mice by a mode of action (MOA) involving activation of the constitutive androstane receptor (CAR). Studies with the model CAR activator phenobarbital (PB) and its sodium salt (sodium phenobarbital; NaPB) have demonstrated that the key and associative events for rat and mouse liver tumour formation include CAR activation, increased hepatocyte replicative DNA synthesis (RDS), induction of cytochrome P450 CYP2B subfamily enzymes, liver hypertrophy, increased altered hepatic foci and hepatocellular adenomas/carcinomas. The key species difference between the rat and mouse compared to humans, is that human hepatocytes are refractory to the mitogenic effects of PB/NaPB and other CAR activators. While PB/NaPB and other CAR activators stimulate RDS in rat and mouse hepatocytes in both in vitro and in vivo studies, such compounds do not stimulate RDS in cultured human hepatocytes and in in vivo studies performed in chimeric mice with humanised livers. In terms of species differences in RDS, unlike the rat and mouse, humans are similar to other species such as the Syrian hamster and guinea pig in being nonresponsive to the mitogenic effects of CAR activators. Overall, the MOA for rat and mouse liver tumour formation by PB/NaPB and other CAR activators is considered qualitatively not plausible for humans. This conclusion is supported by data from a number of epidemiological studies, which demonstrate that chronic treatment with PB does not increase the incidence of liver cancer in humans.

From the Cover: Genomic Effects of Androstenedione and Sex-Specific Liver Cancer Susceptibility in Mice

Current strategies for predicting carcinogenic mode of action for nongenotoxic chemicals are based on identification of early key events in toxicity pathways. The goal of this study was to evaluate short-term key event indicators resulting from exposure to androstenedione (A4), an androgen receptor agonist and known liver carcinogen in mice. Liver cancer is more prevalent in men compared with women, but androgen-related pathways underlying this sex difference have not been clearly identified. Short-term hepatic effects of A4 were compared with reference agonists of the estrogen receptor (ethinyl estradiol, EE) and glucocorticoid receptor (prednisone, PRED). Male B6C3F1 mice were exposed for 7 or 28 days to A4, EE, or PRED. EE increased and PRED suppressed hepatocyte proliferation, while A4 had no detectable effects. In a microarray analysis, EE and PRED altered >3000 and >670 genes, respectively, in a dose-dependent manner, whereas A4 did not significantly alter any genes. Gene expression was subsequently examined in archival liver samples from male and female B6C3F1 mice exposed to A4 for 90 days. A4 altered more genes in females than males and did not alter expression of genes linked to activation of the mitogenic xenobiotic receptors AhR, CAR, and PPARα in either sex. A gene expression biomarker was used to show that in female mice, the high dose of A4 activated the growth hormone-regulated transcription factor STAT5b, which controls sexually dimorphic gene expression in the liver. These findings suggest that A4 induces subtle age-related effects on STAT5b signaling that may contribute to the higher risk of liver cancer in males compared with females.

Sex-Differential Responses of Tumor Promotion-Associated Genes and Dysregulation of Novel Long Noncoding RNAs in Constitutive Androstane Receptor-Activated Mouse Liver

Xenobiotic agonists of constitutive androstane receptor (CAR) induce many hepatic drug metabolizing enzymes, but following prolonged exposure, promote hepatocellular carcinoma, most notably in male mouse liver. Here, we used nuclear RNA-seq to characterize global changes in the mouse liver transcriptome following exposure to the CAR-specific agonist ligand 1,4-bis-[2-(3,5-dichloropyridyloxy)]benzene (TCPOBOP), including changes in novel long noncoding RNAs that may contribute to xenobiotic-induced pathophysiology. Protein-coding genes dysregulated by 3 h TCPOBOP exposure were strongly enriched in KEGG pathways of xenobiotic and drug metabolism, with stronger and more extensive gene responses observed in female than male liver. After 27 h TCPOBOP exposure, the number of responsive genes increased >8-fold in males, where the top enriched pathways and their upstream regulators expanded to include factors implicated in cell cycle dysregulation and hepatocellular carcinoma progression (cyclin-D1, oncogenes E2f, Yap, Rb, Myc, and proto-oncogenes β-catenin, FoxM1, FoxO1, all predicted to be activated by TCPOBOP in male but not female liver; and tumor suppressors p21 and p53, both predicted to be inhibited). Upstream regulators uniquely associated with 3 h TCPOBOP-exposed females include TNF/NFkB pathway members, which negatively regulate CAR-dependent proliferative responses and may contribute to the relative resistance of female liver to TCPOBOP-induced tumor promotion. These responses may be modified by the many long noncoding liver RNAs we show are dysregulated by TCPOBOP or pregnane-X-receptor agonist exposure, including lncRNAs proximal to CAR target genes Cyp2b10, Por, and Alas1. These data provide a comprehensive view of the CAR-regulated transcriptome and give insight into the mechanism of sex-biased susceptibility to CAR-dependent mouse liver tumorigenesis.

Epidermal Growth Factor Receptor Signaling Disruption by Endocrine and Metabolic Disrupting Chemicals

The purpose of this study is to identify an environmentally relevant shared receptor target for endocrine and metabolism disrupting chemical pollutants. A feature of the tested chemicals was that they induced Cyp2b10 in vivo implicating activation of the constitutive androstane receptor (CAR). Recent studies suggest that these compounds could be indirect CAR activators via epidermal growth factor receptor (EGFR) inhibition. Assays included a CAR activity reporter assay, EGF endocytosis assay, and EGFR phosphorylation assay. Docking simulations were used to identify putative binding sites for environmental chemicals on the EGFR. Whole-weight and lipid-adjusted serum mean pollutant exposures were determined using data from the National Health and Examination Survey (NHANES) and compared with the IC50 values determined in vitro. Chlordane, trans-nonachlor, PCB-126, PCB-153, and atrazine were the most potent EGFR inhibitors tested. PCB-126, PCB-153, and trans-nonachlor appeared to be competitive EGFR antagonists as they displaced bound EGF from EGFR. However, atrazine acted through a different mechanism and could be an EGFR tyrosine kinase inhibitor. EGFR inhibition relative effect potencies were determined for these compounds. In NHANES, serum concentrations of trans-nonachlor, PCB-126, and PCB-153 greatly exceeded their calculated IC50 values. A common mechanism of action through EGFR inhibition for three diverse classes of metabolic disrupting chemicals was characterized by measuring inhibition of EGFR phosphorylation and EGF-EGFR endocytosis. Based on NHANES data, EGFR inhibition may be an environmentally relevant mode of action for some PCBs, pesticides, and herbicides.

The embryonic developmental effect of sedaxane on zebrafish (Danio rerio)

The succinate dehydrogenase inhibitor (SDHI) fungicides have been extensively used in agriculture, and some of their potential ecological risks to aquatic organisms have been demonstrated recently. Sedaxane (SDX) is a broad spectrum SDHI fungicide. Despite being extensively used in environment, little is known about its potential developmental effect in zebrafish embryo. This study examined the effects of which SDX triggered in zebrafish through embryonic development assessments. Results show that SDX induced mortality, hatch delay and failure in zebrafish embryos, which were concentration dependent. In addition, several developmental abnormalities were observed at 2 mg/L and higher concentrations, including edema, microcephaly, body deformation, and swim bladder not fully inflated. SDX exposure influenced reactive oxygen species, malondialdehyde, peroxidase, glutathione S-transferase, superoxide dismutase and glutathione in live larvae, which indicated that oxidative stress was caused in zebrafish. Furthermore, SDX induced decrease of succinate dehydrogenase activity in zebrafish larvae. These results provide toxicological data of SDX on developing zebrafish embryo, which could be help for further understanding the potential risk on the environment.

DNA methylation as a genomic marker of exposure to chemical and environmental agents

Recent progress in interpreting comprehensive genetic and epigenetic profiles for human cellular states has contributed new insights into the developmental origins of disease, elucidated novel signalling pathways and enhanced drug discovery programs. A similar comprehensive approach to decoding the epigenetic readouts from chemical challenges in vivo would yield new paradigms for monitoring and assessing environmental exposure in model systems and humans.

Dual modulation of human hepatic zonation via canonical and non-canonical Wnt pathways

The hepatic lobule is divided into three zones along the portal-central vein axis. Hepatocytes within each zone exhibit a distinctive gene expression profile that coordinates their metabolic compartmentalization. The zone-dependent heterogeneity of hepatocytes has been hypothesized to result from the differential degree of exposure to oxygen, nutrition and gut-derived toxins. In addition, the gradient of Wnt signaling that increases towards the central vein seen in rodent models is believed to play a critical role in shaping zonation. Furthermore, hepatic zonation is coupled to the site of the homeostatic renewal of hepatocytes. Despite its critical role, the regulatory mechanisms that determine the distinctive features of zonation and its relevance to humans are not well understood. The present study first conducted a comprehensive zone-dependent transcriptome analysis of normal human liver using laser capture microdissection. Upstream pathway analysis revealed the signatures of host responses to gut-derived toxins in the periportal zone, while both the canonical Wnt pathway and the xenobiotic response pathway govern the perivenular zone. Furthermore, we found that the hypoxic environment of the perivenular zone promotes Wnt11 expression in hepatocytes, which then regulates unique gene expression via activation of the non-canonical Wnt pathway. In summary, our study reports the comprehensive zonation-dependent transcriptome of the normal human liver. Our analysis revealed that the LPS response pathway shapes the characteristics of periportal hepatocytes. By contrast, the perivenular zone is regulated by a combination of three distinct pathways: the xenobiotic response pathway, canonical Wnt signaling, and hypoxia-induced noncanonical Wnt signaling.

The PPARα-dependent rodent liver tumor response is not relevant to humans: addressing misconceptions

A number of industrial chemicals and therapeutic agents cause liver tumors in rats and mice by activating the nuclear receptor peroxisome proliferator-activated receptor α (PPARα). The molecular and cellular events by which PPARα activators induce rodent hepatocarcinogenesis have been extensively studied elucidating a number of consistent mechanistic changes linked to the increased incidence of liver neoplasms. The weight of evidence relevant to the hypothesized mode of action (MOA) for PPARα activator-induced rodent hepatocarcinogenesis is summarized here. Chemical-specific and mechanistic data support concordance of temporal and dose-response relationships for the key events associated with many PPARα activators. The key events (KE) identified in the MOA are PPARα activation (KE1), alteration in cell growth pathways (KE2), perturbation of hepatocyte growth and survival (KE3), and selective clonal expansion of preneoplastic foci cells (KE4), which leads to the apical event-increases in hepatocellular adenomas and carcinomas (KE5). In addition, a number of concurrent molecular and cellular events have been classified as modulating factors, because they potentially alter the ability of PPARα activators to increase rodent liver cancer while not being key events themselves. These modulating factors include increases in oxidative stress and activation of NF-kB. PPARα activators are unlikely to induce liver tumors in humans due to biological differences in the response of KEs downstream of PPARα activation. This conclusion is based on minimal or no effects observed on cell growth pathways and hepatocellular proliferation in human primary hepatocytes and absence of alteration in growth pathways, hepatocyte proliferation, and tumors in the livers of species (hamsters, guinea pigs and cynomolgus monkeys) that are more appropriate human surrogates than mice and rats at overlapping dose levels. Despite this overwhelming body of evidence and almost universal acceptance of the PPARα MOA and lack of human relevance, several reviews have selectively focused on specific studies that, as discussed, contradict the consensus opinion and suggest uncertainty. In the present review, we systematically address these most germane suggested weaknesses of the PPARα MOA.

Applying 'omics technologies in chemicals risk assessment: Report of an ECETOC workshop

Prevailing knowledge gaps in linking specific molecular changes to apical outcomes and methodological uncertainties in the generation, storage, processing, and interpretation of 'omics data limit the application of 'omics technologies in regulatory toxicology. Against this background, the European Centre for Ecotoxicology and Toxicology of Chemicals (ECETOC) convened a workshop Applying 'omics technologies in chemicals risk assessment that is reported herein. Ahead of the workshop, multi-expert teams drafted frameworks on best practices for (i) a Good-Laboratory Practice-like context for collecting, storing and curating 'omics data; (ii) the processing of 'omics data; and (iii) weight-of-evidence approaches for integrating 'omics data. The workshop participants confirmed the relevance of these Frameworks to facilitate the regulatory applicability and use of 'omics data, and the workshop discussions provided input for their further elaboration. Additionally, the key objective (iv) to establish approaches to connect 'omics perturbations to phenotypic alterations was addressed. Generally, it was considered promising to strive to link gene expression changes and pathway perturbations to the phenotype by mapping them to specific adverse outcome pathways. While further work is necessary before gene expression changes can be used to establish safe levels of substance exposure, the ECETOC workshop provided important incentives towards achieving this goal.

PPARα-independent transcriptional targets of perfluoroalkyl acids revealed by transcript profiling

Perfluoroalkyl acids (PFAAs) are ubiquitous and persistent environmental contaminants. Compounds such as perfluoroocanoic acid (PFOA), perfluorooctane sulfonate (PFOS), perfluorononanoic acid (PFNA), and perfluorohexane sulfonate (PFHxS) are readily found in the tissues of humans and wildlife. While PFOA and PFOS have been the subject of numerous studies since they were first described over a decade ago, less is known about the biological activity of PFHxS and PFNA. Most PFAAs are activators of peroxisome proliferator-activated receptor α (PPARα), although the biological effects of these compounds are likely mediated by other factors in addition to PPARα. To evaluate the effects of PFHxS and PFNA, male wild-type and Pparα-null mice were dosed by oral gavage with PFHxS (3 or 10mg/kg/day), PFNA (1 or 3mg/kg/day), or vehicle for 7days, and liver gene expression was evaluated by full-genome microarrays. Gene expression patterns were then compared to historical in-house data for PFOA and PFOS in addition to the experimental hypolipidemic agent, WY-14,643. While WY-14,643 altered most genes in a PPARα-dependent manner, approximately 11-24% of regulated genes in PFAA-treated mice were independent of PPARα. The possibility that PFAAs regulate gene expression through other molecular pathways was evaluated. Using data available through a microarray database, PFAA gene expression profiles were found to exhibit significant similarity to profiles from mouse tissues exposed to agonists of the constitutive activated receptor (CAR), estrogen receptor α (ERα), and PPARγ. Human PPARγ and ERα were activated by all four PFAAs in trans-activation assays from the ToxCast screening program. Predictive gene expression biomarkers showed that PFAAs activate CAR in both genotypes and cause feminization of the liver transcriptome through suppression of signal transducer and activator of transcription 5B (STAT5B). These results indicate that, in addition to activating PPARα as a primary target, PFAAs also have the potential to activate CAR, PPARγ, and ERα as well as suppress STAT5B.

RNA Sequencing Reveals Age and Species Differences of Constitutive Androstane Receptor-Targeted Drug-Processing Genes in the Liver

The constitutive androstane receptor (CAR/Nr1i3) is an important xenobiotic-sensing nuclear receptor that is highly expressed in the liver and is well known to have species differences. During development, age-specific activation of CAR may lead to modified pharmacokinetics and toxicokinetics of drugs and environmental chemicals, leading to higher risks for adverse drug reactions in newborns and children. The goal of this study was to systematically investigate the age- and species-specific regulation of various drug-processing genes (DPGs) after neonatal or adult CAR activation in the livers of wild-type, CAR-null, and humanized CAR transgenic mice. At either 5 or 60 days of age, the three genotypes of mice were administered a species-appropriate CAR ligand or vehicle once daily for 4 days (i.p.). The majority of DPGs were differentially regulated by age and/or CAR activation. Thirty-six DPGs were commonly upregulated by CAR activation regardless of age or species of CAR. Although the cumulative mRNAs of uptake transporters were not readily altered by CAR, the cumulative phase I and phase II enzymes as well as efflux transporters were all increased after CAR activation in both species. In general, mouse CAR activation produced comparable or even greater fold increases of many DPGs in newborns than in adults; conversely, humanized CAR activation produced weaker induction in newborns than in adults. Western blotting and enzyme activity assays confirmed the age and species specificities of selected CAR-targeted DPGs. In conclusion, this study systematically compared the effect of age and species of CAR proteins on the regulation of DPGs in the liver and demonstrated that the regulation of xenobiotic biotransformation by CAR is profoundly modified by age and species.

Perfluoroalkyl acids-induced liver steatosis: Effects on genes controlling lipid homeostasis

Persistent presence of perfluoroalkyl acids (PFAAs) in the environment is due to their extensive use in industrial and consumer products, and their slow decay. Biochemical tests in rodent demonstrated that these chemicals are potent modifiers of lipid metabolism and cause hepatocellular steatosis. However, the molecular mechanism of PFAAs interference with lipid metabolism remains to be elucidated. Currently, two major hypotheses are that PFAAs interfere with mitochondrial beta-oxidation of fatty acids and/or they affect the transcriptional activity of peroxisome proliferator-activated receptor α (PPARα) in liver. To determine the ability of structurally-diverse PFAAs to cause steatosis, as well as to understand the underlying molecular mechanisms, wild-type (WT) and PPARα-null mice were treated with perfluorooctanoic acid (PFOA), perfluorononanoic acid (PFNA), or perfluorohexane sulfonate (PFHxS), by oral gavage for 7days, and their effects were compared to that of PPARα agonist WY-14643 (WY), which does not cause steatosis. Increases in liver weight and cell size, and decreases in DNA content per mg of liver, were observed for all compounds in WT mice, and were also seen in PPARα-null mice for PFOA, PFNA, and PFHxS, but not for WY. In Oil Red O stained sections, WT liver showed increased lipid accumulation in all treatment groups, whereas in PPARα-null livers, accumulation was observed after PFNA and PFHxS treatment, adding to the burden of steatosis observed in control (untreated) PPARα-null mice. Liver triglyceride (TG) levels were elevated in WT mice by all PFAAs and in PPARα-null mice only by PFNA. In vitro β-oxidation of palmitoyl carnitine by isolated rat liver mitochondria was not inhibited by any of the 7 PFAAs tested. Likewise, neither PFOA nor PFOS inhibited palmitate oxidation by HepG2/C3A human liver cell cultures. Microarray analysis of livers from PFAAs-treated mice indicated that the PFAAs induce the expression of the lipid catabolism genes, as well as those involved in fatty acid and triglyceride synthesis, in WT mice and, to a lesser extent, in PPARα-null mice. These results indicate that most of the PFAAs increase liver TG load and promote steatosis in mice We hypothesize that PFAAs increase steatosis because the balance of fatty acid accumulation/synthesis and oxidation is disrupted to favor accumulation.

p21Cip1 plays a critical role in the physiological adaptation to fasting through activation of PPARα

Fasting is a physiological stress that elicits well-known metabolic adaptations, however, little is known about the role of stress-responsive tumor suppressors in fasting. Here, we have examined the expression of several tumor suppressors upon fasting in mice. Interestingly, p21 mRNA is uniquely induced in all the tissues tested, particularly in liver and muscle (>10 fold), and this upregulation is independent of p53. Remarkably, in contrast to wild-type mice, p21-null mice become severely morbid after prolonged fasting. The defective adaptation to fasting of p21-null mice is associated to elevated energy expenditure, accelerated depletion of fat stores, and premature activation of protein catabolism in the muscle. Analysis of the liver transcriptome and cell-based assays revealed that the absence of p21 partially impairs the transcriptional program of PPARα, a key regulator of fasting metabolism. Finally, treatment of p21-null mice with a PPARα agonist substantially protects them from their accelerated loss of fat upon fasting. We conclude that p21 plays a relevant role in fasting adaptation through the positive regulation of PPARα.

Development of an in vitro high content imaging assay for quantitative assessment of CAR-dependent mouse, rat, and human primary hepatocyte proliferation

Rodent liver tumors promoted by constitutive androstane receptor (CAR) activation are known to be mediated by key events that include CAR-dependent gene expression and hepatocellular proliferation. Here, an in vitro high content imaging based assay was developed for quantitative assessment of nascent DNA synthesis in primary hepatocyte cultures from mouse, rat, and human species. Detection of DNA synthesis was performed using direct DNA labeling with the nucleoside analog 5-ethynyl-2'-deoxyuridine (EdU). The assay was multiplexed to enable direct quantitation of DNA synthesis, cytotoxicity, and cell count endpoints. An optimized defined medium cocktail was developed to sensitize hepatocytes to cell cycle progression. The baseline EdU response to defined medium was greatest for mouse, followed by rat, and then human. Hepatocytes from all three species demonstrated CAR activation in response to the CAR agonists TCPOBOP, CITCO, and phenobarbital based on increased gene expression for Cyp2b isoforms. When evaluated for a proliferation phenotype, TCPOBOP and CITCO exhibited significant dose-dependent increases in frequency of EdU labeling in mouse and rat hepatocytes that was not observed in hepatocytes from three human donors. The observed species differences are consistent with CAR activators inducing a proliferative response in rodents, a key event in the liver tumor mode of action that is not observed in humans.

Regulation of Drug Disposition Gene Expression in Pregnant Mice with Car Receptor Activation

More than half of pregnant women use prescription medications in order to maintain both maternal and fetal health. The constitutive androstane receptor (Car) critically affects the disposition of chemicals by regulating the transcription of genes encoding metabolic enzymes and transporters. However, the effects of Car activation on chemical disposition during pregnancy are unclear. This study aims to determine the degree to which pregnancy alters the expression of drug metabolizing enzymes and transporters in response to the pharmacological activation of Car. To test this, pregnant C57BL/6 mice were administered IP doses of vehicle, or a potent Car agonist, TCPOBOP, on gestation days 14, 15 and 16. Hepatic mRNA and protein expression of Car target genes (phase I, II and transporters) were quantified on gestation day 17. Pregnancy-related changes, such as induction of Cyp2b10, Ugt1a1 and Sult1a1 and repression of Ugt1a6, Gsta1, Gsta2 and Mrp6, were observed. Interestingly, the induction of Cyp2b10, Gsta1, Gsta2 and Mrp2-4 mRNAs by TCPOBOP was attenuated in maternal livers suggesting that Car activation is impeded by the biochemical and/or physiological changes that occur during gestation. Taken together, these findings suggest that pregnancy and pharmacological activation of Car can differentially regulate the expression of drug metabolism and transport genes.

Application of the TGx-28.65 transcriptomic biomarker to classify genotoxic and non-genotoxic chemicals in human TK6 cells in the presence of rat liver S9

In vitro transcriptional signatures that predict toxicities can facilitate chemical screening. We previously developed a transcriptomic biomarker (known as TGx-28.65) for classifying agents as genotoxic (DNA damaging) and non-genotoxic in human lymphoblastoid TK6 cells. Because TK6 cells do not express cytochrome P450s, we confirmed accurate classification by the biomarker in cells co-exposed to 1% 5,6 benzoflavone/phenobarbital-induced rat liver S9 for metabolic activation. However, chemicals may require different types of S9 for activation. Here we investigated the response of TK6 cells to higher percentages of Aroclor-, benzoflavone/phenobarbital-, or ethanol-induced rat liver S9 to expand TGx-28.65 biomarker applicability. Transcriptional profiles were derived 3 to 4 hr following a 4 hr co-exposure of TK6 cells to test chemicals and S9. Preliminary studies established that 10% Aroclor- and 5% ethanol-induced S9 alone did not induce the TGx-28.65 biomarker genes. Seven genotoxic and two non-genotoxic chemicals (and concurrent solvent and positive controls) were then tested with one of the S9s (selected based on cell survival and micronucleus induction). Relative survival and micronucleus frequency was assessed by flow cytometry in cells 20 hr post-exposure. Genotoxic/non-genotoxic chemicals were accurately classified using the different S9s. One technical replicate of cells co-treated with dexamethasone and 10% Aroclor-induced S9 was falsely classified as genotoxic, suggesting caution in using high S9 concentrations. Even low concentrations of genotoxic chemicals (those not causing cytotoxicity) were correctly classified, demonstrating that TGx-28.65 is a sensitive biomarker of genotoxicity. A meta-analysis of datasets from 13 chemicals supports that different S9s can be used in TK6 cells, without impairing classification using the TGx-28.65 biomarker.

RNA-Seq reveals common and unique PXR- and CAR-target gene signatures in the mouse liver transcriptome

The pregnane X receptor (PXR) and constitutive androstane receptor (CAR) are well-known xenobiotic-sensing nuclear receptors with overlapping functions. However, there lacks a quantitative characterization to distinguish between the PXR and CAR target genes and signaling pathways in the liver. The present study performed a transcriptomic comparison of the PXR- and CAR-targets using RNA-Seq in livers of adult wild-type mice that were treated with the prototypical PXR ligand PCN (200mg/kg, i.p. once daily for 4days in corn oil) or the prototypical CAR ligand TCPOBOP (3mg/kg, i.p., once daily for 4days in corn oil). At the given doses, TCPOBOP differentially regulated many more genes (2125) than PCN (212), and 147 of the same genes were differentially regulated by both chemicals. As expected, the top pathways differentially regulated by both PCN and TCPOBOP were involved in xenobiotic metabolism, and they also up-regulated genes involved in retinoid metabolism, but down-regulated genes involved in inflammation and iron homeostasis. Regarding unique pathways, PXR activation appeared to overlap with the aryl hydrocarbon receptor signaling, whereas CAR activation appeared to overlap with the farnesoid X receptor signaling, acute-phase response, and mitochondrial dysfunction. The mRNAs of differentially regulated drug-processing genes (DPGs) partitioned into three patterns, namely TCPOBOP-induced, PCN-induced, as well as TCPOBOP-suppressed gene clusters. The cumulative mRNAs of the differentially regulated DPGs, phase-I and -II enzymes, as well as efflux transporters were all up-regulated by both PCN and TCPOBOPOP, whereas the cumulative mRNAs of the uptake transporters were down-regulated only by TCPOBOP. The absolute mRNA abundance in control and receptor-activated conditions was examined in each DPG category to predict the contribution of specific DPG genes in the PXR/CAR-mediated pharmacokinetic responses. The preferable differential regulation by TCPOBOP in the entire hepatic transcriptome correlated with a marked change in the expression of many DNA and histone epigenetic modifiers. In conclusion, the present study has revealed known and novel, as well as common and unique targets of PXR and CAR in mouse liver following pharmacological activation using their prototypical ligands. Results from this study will further support the role of these receptors in regulating the homeostasis of xenobiotic and intermediary metabolism in the liver, and aid in distinguishing between PXR and CAR signaling at various physiological and pathophysiological conditions. This article is part of a Special Issue entitled: Xenobiotic nuclear receptors: New Tricks for An Old Dog, edited by Dr. Wen Xie.

Disruption of STAT5b-Regulated Sexual Dimorphism of the Liver Transcriptome by Diverse Factors Is a Common Event

Signal transducer and activator of transcription 5b (STAT5b) is a growth hormone (GH)-activated transcription factor and a master regulator of sexually dimorphic gene expression in the liver. Disruption of the GH hypothalamo-pituitary-liver axis controlling STAT5b activation can lead to metabolic dysregulation, steatosis, and liver cancer. Computational approaches were developed to identify factors that disrupt STAT5b function in a mouse liver gene expression compendium. A biomarker comprised of 144 STAT5b-dependent genes was derived using comparisons between wild-type male and wild-type female mice and between STAT5b-null and wild-type mice. Correlations between the STAT5b biomarker gene set and a test set comprised of expression datasets (biosets) with known effects on STAT5b function were evaluated using a rank-based test (the Running Fisher algorithm). Using a similarity p-value ≤ 10(-4), the test achieved a balanced accuracy of 99% and 97% for detection of STAT5b activation or STAT5b suppression, respectively. The STAT5b biomarker gene set was then used to identify factors that activate (masculinize) or suppress (feminize) STAT5b function in an annotated mouse liver and primary hepatocyte gene expression compendium of ~1,850 datasets. Disruption of GH-regulated STAT5b is a common phenomenon in liver in vivo, with 5% and 29% of the male datasets, and 11% and 13% of the female datasets, associated with masculinization or feminization, respectively. As expected, liver STAT5b activation/masculinization occurred at puberty and suppression/feminization occurred during aging and in mutant mice with defects in GH signaling. A total of 70 genes were identified that have effects on STAT5b activation in genetic models in which the gene was inactivated or overexpressed. Other factors that affected liver STAT5b function were shown to include fasting, caloric restriction and infections. Together, these findings identify diverse factors that perturb the hypothalamo-pituitary-liver GH axis and disrupt GH-dependent STAT5b activation in mouse liver.

Chemical and Hormonal Effects on STAT5b-Dependent Sexual Dimorphism of the Liver Transcriptome

The growth hormone (GH)-activated transcription factor signal transducer and activator of transcription 5b (STAT5b) is a key regulator of sexually dimorphic gene expression in the liver. Suppression of hepatic STAT5b signaling is associated with lipid metabolic dysfunction leading to steatosis and liver cancer. In the companion publication, a STAT5b biomarker gene set was identified and used in a rank-based test to predict both increases and decreases in liver STAT5b activation status/function with high (≥ 97%) accuracy. Here, this computational approach was used to identify chemicals and hormones that activate (masculinize) or suppress (feminize) STAT5b function in a large, annotated mouse liver and primary hepatocyte gene expression compendium. Exposure to dihydrotestosterone and thyroid hormone caused liver masculinization, whereas glucocorticoids, fibroblast growth factor 15, and angiotensin II caused liver feminization. In mouse models of diabetes and obesity, liver feminization was consistently observed and was at least partially reversed by leptin or resveratrol exposure. Chemical-induced feminization of male mouse liver gene expression profiles was a relatively frequent phenomenon: of 156 gene expression biosets from chemically-treated male mice, 29% showed feminization of liver STAT5b function, while <1% showed masculinization. Most (93%) of the biosets that exhibited feminization of male liver were also associated with activation of one or more xenobiotic-responsive receptors, most commonly constitutive activated receptor (CAR) or peroxisome proliferator-activated receptor alpha (PPARα). Feminization was consistently associated with increased expression of peroxisome proliferator-activated receptor gamma (Pparg) but not other lipogenic transcription factors linked to steatosis. GH-activated STAT5b signaling in mouse liver is thus commonly altered by diverse chemicals, and provides a linkage between chemical exposure and dysregulated gene expression associated with adverse effects on the liver.

Moving Toward Integrating Gene Expression Profiling Into High-Throughput Testing: A Gene Expression Biomarker Accurately Predicts Estrogen Receptor α Modulation in a Microarray Compendium

Microarray profiling of chemical-induced effects is being increasingly used in medium- and high-throughput formats. Computational methods are described here to identify molecular targets from whole-genome microarray data using as an example the estrogen receptor α (ERα), often modulated by potential endocrine disrupting chemicals. ERα biomarker genes were identified by their consistent expression after exposure to 7 structurally diverse ERα agonists and 3 ERα antagonists in ERα-positive MCF-7 cells. Most of the biomarker genes were shown to be directly regulated by ERα as determined by ESR1 gene knockdown using siRNA as well as through chromatin immunoprecipitation coupled with DNA sequencing analysis of ERα-DNA interactions. The biomarker was evaluated as a predictive tool using the fold-change rank-based Running Fisher algorithm by comparison to annotated gene expression datasets from experiments using MCF-7 cells, including those evaluating the transcriptional effects of hormones and chemicals. Using 141 comparisons from chemical- and hormone-treated cells, the biomarker gave a balanced accuracy for prediction of ERα activation or suppression of 94% and 93%, respectively. The biomarker was able to correctly classify 18 out of 21 (86%) ER reference chemicals including "very weak" agonists. Importantly, the biomarker predictions accurately replicated predictions based on 18 in vitro high-throughput screening assays that queried different steps in ERα signaling. For 114 chemicals, the balanced accuracies were 95% and 98% for activation or suppression, respectively. These results demonstrate that the ERα gene expression biomarker can accurately identify ERα modulators in large collections of microarray data derived from MCF-7 cells.

Screening a mouse liver gene expression compendium identifies modulators of the aryl hydrocarbon receptor (AhR)

The aryl hydrocarbon receptor (AhR) is a ligand-activated transcription factor that mediates the biological and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), dioxin-like compounds (DLC) as well as some drugs and endogenous tryptophan metabolites. Short-term activation of AhR can lead to hepatocellular steatosis, and chronic activation can lead to liver cancer in mice and rats. Analytical approaches were developed to identify biosets in a genomic database in which AhR activity was altered. A set of 63 genes was identified (the AhR gene expression biomarker) that was dependent on AhR for regulation after exposure to TCDD or benzo[a]pyrene and includes the known AhR targets Cyp1a1 and Cyp1b1. A fold-change rank-based test (Running Fisher's test; p-value ≤ 10(-4)) was used to evaluate the similarity between the AhR biomarker and a test set of 37 and 41 biosets positive or negative, respectively for AhR activation. The test resulted in a balanced accuracy of 95%. The rank-based test was used to identify factors that activate or suppress AhR in an annotated mouse liver/mouse primary hepatocyte gene expression database of ∼ 1850 comparisons. In addition to the expected activation of AhR by TCDD and DLC, AhR was activated by AP20189 and phenformin. AhR was suppressed by phenobarbital and 1,4-Bis[2-(3,5-dichloropyridyloxy)] benzene (TCPOBOP) in a constitutive activated receptor (CAR)-dependent manner and pregnenolone-16α-carbonitrile in a pregnane X receptor (PXR)-dependent manner. Inactivation of individual genes in nullizygous models led to AhR activation (Pxr, Ghrhr, Taf10) or suppression (Ahr, Ilst6st, Hnf1a). This study describes a novel screening strategy for identifying factors in mouse liver that perturb AhR in a gene expression compendium.