Pharmacological activation of constitutive androstane receptor induces female-specific modulation of hepatic metabolism

Background & Aims

The constitutive androstane receptor (CAR) is a nuclear receptor that binds diverse xenobiotics and whose activation leads to the modulation of the expression of target genes involved in xenobiotic detoxification and energy metabolism. Although CAR hepatic activity is considered to be higher in women than in men, its sex-dependent response to an acute pharmacological activation has seldom been investigated.

Methods

The hepatic transcriptome, plasma markers, and hepatic metabolome, were analysed in Car+/+ and Car-/- male and female mice treated either with the CAR-specific agonist 1,4-bis[2-(3,5-dichloropyridyloxy)]benzene (TCPOBOP) or with vehicle.

Results

Although 90% of TCPOBOP-sensitive genes were modulated in a sex-independent manner, the remaining 10% showed almost exclusive female liver specificity. These female-specific CAR-sensitive genes were mainly involved in xenobiotic metabolism, inflammation, and extracellular matrix organisation. CAR activation also induced higher hepatic oxidative stress and hepatocyte cytolysis in females than in males. Hepatic expression of flavin monooxygenase 3 (Fmo3) was almost abolished and was associated with a decrease in hepatic trimethylamine-N-oxide (TMAO) concentration in TCPOBOP-treated females. In line with a potential role in the control of TMAO homeostasis, CAR activation decreased platelet hyper-responsiveness in female mice supplemented with dietary choline.

Conclusions

More than 10% of CAR-sensitive genes are sex-specific and influence hepatic and systemic responses such as platelet aggregation. CAR activation may be an important mechanism of sexually-dimorphic drug-induced liver injury.

Impact and implications

CAR is activated by many drugs and pollutants. Its pharmacological activation had a stronger impact on hepatic gene expression and metabolism in females than in males, and had a specific impact on liver toxicity and trimethylamine metabolism. Sexual dimorphism should be considered when testing and/or prescribing xenobiotics known to activate CAR.

An integrated assessment of the 1,4-dioxane cancer mode of action and threshold response in rodents

1,4-Dioxane is an environmental contaminant that has been shown to cause cancer in rodents after chronic high dose exposures. We reviewed and integrated information from recently published studies to update our understanding of the cancer mode of action of 1,4-dioxane. Tumor development in rodents from exposure to high doses of 1,4-dioxane is preceded by pre-neoplastic events including increased hepatic genomic signaling activity related to mitogenesis, elevation of Cyp2E1 activity and oxidative stress leading to genotoxicity and cytotoxicity. These events are followed by regenerative repair and proliferation and eventual development of tumors. Importantly, these events occur at doses that exceed the metabolic clearance of absorbed 1,4-dioxane in rats and mice resulting in elevated systemic levels of parent 1,4-dioxane. Consistent with previous reviews, we found no evidence of direct mutagenicity from exposure to 1,4-dioxane. We also found no evidence of CAR/PXR, AhR or PPARα activation resulting from exposure to 1,4-dioxane. This integrated assessment supports a cancer mode of action that is dependent on exceeding the metabolic clearance of absorbed 1,4-dioxane, direct mitogenesis, elevation of Cyp2E1 activity and oxidative stress leading to genotoxicity and cytotoxicity followed by sustained proliferation driven by regenerative repair and progression of heritable lesions to tumor development.

A tiered testing strategy based on in vitro phenotypic and transcriptomic data for selecting representative petroleum UVCBs for toxicity evaluation in vivo

Hazard evaluation of substances of "unknown or variable composition, complex reaction products and biological materials" (UVCBs) remains a major challenge in regulatory science because their chemical composition is difficult to ascertain. Petroleum substances are representative UVCBs and human cell-based data have been previously used to substantiate their groupings for regulatory submissions. We hypothesized that a combination of phenotypic and transcriptomic data could be integrated to make decisions as to selection of group-representative worst-case petroleum UVCBs for subsequent toxicity evaluation in vivo. We used data obtained from 141 substances from 16 manufacturing categories previously tested in 6 human cell types (induced pluripotent stem cell [iPSC]-derived hepatocytes, cardiomyocytes, neurons, and endothelial cells, and MCF7 and A375 cell lines). Benchmark doses for gene-substance combinations were calculated, and both transcriptomic and phenotype-derived points of departure (PODs) were obtained. Correlation analysis and machine learning were used to assess associations between phenotypic and transcriptional PODs and to determine the most informative cell types and assays, thus representing a cost-effective integrated testing strategy. We found that 2 cell types-iPSC-derived-hepatocytes and -cardiomyocytes-contributed the most informative and protective PODs and may be used to inform selection of representative petroleum UVCBs for further toxicity evaluation in vivo. Overall, although the use of new approach methodologies to prioritize UVCBs has not been widely adopted, our study proposes a tiered testing strategy based on iPSC-derived hepatocytes and cardiomyocytes to inform selection of representative worst-case petroleum UVCBs from each manufacturing category for further toxicity evaluation in vivo.

Bridging Sex-Specific Differences in the CAR-Mediated Hepatocarcinogenesis of Nitrapyrin Using Molecular and Apical Endpoints

Nitrapyrin, a nitrification inhibitor, produces liver tumors in B6C3F1 mice. In a 2-year oncogenicity study, increased incidence of mice with hepatocellular tumors was observed following exposure to 125 (females only) or 250 mg/kg/day (males and females) nitrapyrin in the diet. Previous data was generated in male mice to support a mode-of-action (MoA) characterized by constitutive androstane receptor (CAR) nuclear receptor (NR) activation, increased hepatocellular proliferation, and subsequent hepatocellular foci and tumor formation. Uncertainty as to the relevance of this MoA for females remained given the increased sensitivity to tumor formation in female mice. A targeted MoA study was conducted to evaluate CAR activation and hepatic responses in female mice treated with the female carcinogenic dose of nitrapyrin for 4 days. Nitrapyrin induced a treatment-related increase in hepatocellular hypertrophy and hepatocellular proliferation. Nitrapyrin also induced a dose-related increase in the Cyp2b10/CAR-associated transcript and liver weights. Nitrapyrin-induced liver weights and Cyp2b10 gene expression for both males and females were compared to data generated from three other established CAR activators; methyl isobutyl ketone, phenobarbital, and sulfoxaflor. The response observed in female mice following exposure to nitrapyrin was within range of the degree of change observed in mice following exposure to tumorigenic doses of other CAR activators. Consistent with the liver MoA in male mice, these data support a CAR-mediated mode of action for nitrapyrin-induced liver tumors in female mice, with the understanding that a focused approach minimizing animal use can bridge male and female datasets when sex-specific carcinogenic differences are observed.

Rethinking chronic toxicity and carcinogenicity assessment for agrochemicals project (ReCAAP): A reporting framework to support a weight of evidence safety assessment without long-term rodent bioassays

Rodent cancer bioassays have been long-required studies for regulatory assessment of human cancer hazard and risk. These studies use hundreds of animals, are resource intensive, and certain aspects of these studies have limited human relevance. The past 10 years have seen an exponential growth of new technologies with the potential to effectively evaluate human cancer hazard and risk while reducing, refining, or replacing animal use. To streamline and facilitate uptake of new technologies, a workgroup comprised of scientists from government, academia, non-governmental organizations, and industry stakeholders developed a framework for waiver rationales of rodent cancer bioassays for consideration in agrochemical safety assessment. The workgroup used an iterative approach, incorporating regulatory agency feedback, and identifying critical information to be considered in a risk assessment-based weight of evidence determination of the need for rodent cancer bioassays. The reporting framework described herein was developed to support a chronic toxicity and carcinogenicity study waiver rationale, which includes information on use pattern(s), exposure scenario(s), pesticidal mode-of-action, physicochemical properties, metabolism, toxicokinetics, toxicological data including mechanistic data, and chemical read-across from similar registered pesticides. The framework could also be applied to endpoints other than chronic toxicity and carcinogenicity, and for chemicals other than agrochemicals.

A 90-day drinking water study in mice to characterize early events in the cancer mode of action of 1,4-dioxane

Studies demonstrate that with sufficient dose and duration, 1,4-dioxane (1,4-DX) induces liver tumors in laboratory rodent models. The available evidence aligns with a threshold-dependent, tumor promotion mode of action (MOA). The MOA and key events (KE) in rats are well developed but less so in the mouse. Therefore, we conducted a 90-day drinking water study in female mice to evaluate early KE at 7, 28, and 90 days. Female B6D2F1/Crl mice consumed drinking water containing 0, 40, 200, 600, 2000 or 6000 ppm 1,4-DX. 1,4-DX was detected in blood at 90-days of exposure to 6000 ppm, but not in the other exposure groups, indicating a metabolic clearance threshold between 2000 and 6000. Early events identified in this study include glycogen-like vacuolization, centrilobular hypertrophy, centrilobular GST-P staining, apoptosis, and pan-lobular increase in cell proliferation observed after 90-days of exposure to 6000 ppm 1,4-DX. There was minimal evidence of hepatotoxicity over the duration of this study. These findings demonstrate a previously unreported direct mitogenic response following exposures exceeding the metabolic clearance threshold of 1,4-DX. Collectively, the information generated in this study supports a threshold MOA for the development of liver tumors in mice after exposure to 1,4-DX.

Direct formalin fixation induces widespread transcriptomic effects in archival tissue samples

Sequencing technologies now provide unprecedented access to genomic information in archival formalin-fixed paraffin-embedded (FFPE) tissue samples. However, little is known about artifacts induced during formalin fixation, which could bias results. Here we evaluated global changes in RNA-sequencing profiles between matched frozen and FFPE samples. RNA-sequencing was performed on liver samples collected from mice treated with a reference chemical (phenobarbital) or vehicle control for 7 days. Each sample was divided into four parts: (1) fresh-frozen, (2) direct-fixed in formalin for 18 h, (3) frozen then formalin-fixed, and (4) frozen then ethanol-fixed and paraffin-embedded (n = 6/group/condition). Direct fixation resulted in 2,946 differentially expressed genes (DEGs) vs. fresh-frozen, 98% of which were down-regulated. Freezing prior to formalin fixation had ≥ 95% fewer DEGs vs. direct fixation, indicating that most formalin-derived transcriptional effects in the liver occurred during fixation. This finding was supported by retrospective studies of paired frozen and FFPE samples, which identified consistent enrichment in oxidative stress, mitochondrial dysfunction, and transcription initiation pathways with direct fixation. Notably, direct formalin fixation in the parent study did not significantly impact response profiles resulting from chemical exposure. These results advance our understanding of FFPE samples as a resource for genomic research.

Chemical Activation of the Constitutive Androstane Receptor Leads to Activation of Oxidant-Induced Nrf2

Exposure to environmentally relevant chemicals that activate the xenobiotic receptors aryl hydrocarbon receptor (AhR), constitutive androstane receptor (CAR), and peroxisome proliferator-activated receptor alpha (PPARα) in rodent test systems often leads to increases in oxidative stress (OS) that contributes to liver cancer induction. We hypothesized that activation of the oxidant-induced transcription factor Nrf2 could be used as a surrogate endpoint for increases in OS. We examined the relationships between activation of xenobiotic receptors and Nrf2 using previously characterized gene expression biomarkers that accurately predict modulation. Using a correlation approach (Running Fisher Test), the biomarkers were compared with microarray profiles in a mouse liver gene expression compendium. Out of the 163 chemicals examined, 47% from 53 studies activated Nrf2. We found consistent coupling between CAR and Nrf2 activation. Out of the 41 chemicals from 32 studies that activated CAR, 90% also activated Nrf2. CAR was activated earlier and at lower doses than Nrf2, indicating CAR activation preceded Nrf2 activation. Nrf2 activation by 2 CAR activators was abolished in CAR-null mice. We hypothesized that Nrf2 is activated by reactive oxygen species from the increased activity of enzymes encoded by Cyp2b family members. However, Nrf2 was similarly activated in the livers of both TCPOBOP-treated wild-type and Cyp2b9/10/13-null mice. This study provides evidence that Nrf2 activation (1) often occurs after exposure to xenobiotic chemicals, (2) is tightly linked to activation of CAR, and (3) does not require induction of 3 Cyp2b genes secondary to CAR activation.

Demodifying RNA for Transcriptomic Analyses of Archival Formalin-Fixed Paraffin-Embedded Samples

Archival formalin-fixed paraffin-embedded (FFPE) tissue samples offer a vast but largely untapped resource for genomic research. The primary technical issues limiting use of FFPE samples are RNA yield and quality. In this study, we evaluated methods to demodify RNA highly fragmented and crosslinked by formalin fixation. Primary endpoints were RNA recovery, RNA-sequencing quality metrics, and transcriptional responses to a reference chemical (phenobarbital, PB). Frozen mouse liver samples from control and PB groups (n = 6/group) were divided and preserved for 3 months as follows: frozen (FR); 70% ethanol (OH); 10% buffered formalin for 18 h followed by ethanol (18F); or 10% buffered formalin (3F). Samples from OH, 18F, and 3F groups were processed to FFPE blocks and sectioned for RNA isolation. Additional sections from 3F received the following demodification protocols to mitigate RNA damage: short heated incubation with Tris-Acetate-EDTA buffer; overnight heated incubation with an organocatalyst using 2 different isolation kits; or overnight heated incubation without organocatalyst. Ribo-depleted, stranded, total RNA libraries were built and sequenced using the Illumina HiSeq 2500 platform. Overnight incubation (± organocatalyst) increased RNA yield >3-fold and RNA integrity numbers and fragment analysis values by > 1.5- and >3.0-fold, respectively, versus 3F. Postsequencing metrics also showed reduced bias in gene coverage and deletion rates for overnight incubation groups. All demodification groups had increased overlap for differentially expressed genes (77%-84%) and enriched pathways (91%-97%) with FR, with the highest overlap in the organocatalyst groups. These results demonstrate simple changes in RNA isolation methods that can enhance genomic analyses of FFPE samples.

Pilot studies evaluating the nongenotoxic rodent carcinogens phenobarbital and clofibrate in the rat Pig-a assay

The Pig-a assay is an emerging and promising in vivo method to determine mutagenic potential of chemicals. Since its development in 2008, remarkable progress has been made in harmonizing and characterizing the test procedures, primarily using known mutagenic chemicals. The purpose of the present study was to evaluate specificity of the Pig-a assay using two nongenotoxic and well-characterized rodent liver carcinogens, phenobarbital and clofibrate, in male F344/DuCrl rats. Daily oral administration of phenobarbital or clofibrate at established hepatotoxic doses for 28 days resulted in substantial hepatic alterations, however, did not increase the frequency of Pig-a mutation markers (RET^CD59- and RBC^CD59- ) compared to vehicle control or pre-exposure (Day -5) mutant frequencies. These results are consistent with the existing literature on the nonmutagenic mode of action (MoA) of phenobarbital and clofibrate liver tumors. The present study contributes to the limited, but expanding evidence on the specificity of the Pig-a assay and further for the investigations of carcinogenic MoAs, i.e., mutagenic or nonmutagenic potential of chemicals. Environ. Mol. Mutagen. 60:42-46, 2019. © 2018 Wiley Periodicals, Inc.

Effects on the hepatic transcriptome of chicken embryos in ovo exposed to phenobarbital

This work aimed at evaluating the toxic effects of in ovo exposure to phenobarbital (PB) and unveiling the mode of action by transcriptome analysis in the embryonic liver of a model avian species, chicken (Gallus gallus). Embryos were initially treated with saline or 1 μg PB /g egg at Hamburger Hamilton Stage (HHS) 1 (1st day), followed by 20 days of incubation to HHS 46. At 21st day, chicks that pipped successfully were euthanized and dissected for assessing the PB caused effects on phenotypes and the liver transcriptome in both genders. In the PB treatment group, a 7% attenuation in tarsus length was found in females. While no adverse phenotypic effect on the liver somatic index (LSI) was observed, PB caused significant changes in the expressions of 52 genes in males and 516 genes in females (False Discovery Rate < 0.2, p value < 0.05, and absolute fold change > 2). PB exposure modulated the genes primarily enriched in the biological pathways of the cancer, cardiac development, immune response, lipid metabolism, and skeletal development in both genders, and altered expressions of genes related to the cellular process and neural development in females. However, mRNA expressions of chicken xenobiotic receptor (CXR)-mediated CYP genes were not induced in the PB treatment groups, regardless of males and females. On the contrary, PB exposure repressed the mRNA expressions of CYP2AC2 in males and CYP2R1, CYP3A37, and CYP8B1 in females. Although transcription factors (TFs) including SREBF1 and COUP-TFII were predicted to be commonly activated in both genders, some TFs were activated in a gender-dependent manner, such as PPARa in males and BRCA1 and IRF9 in females. Taken together, our results provided an insight into the mode of action of PB on the chicken embryos.

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.

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.

Reprint of: Environmental toxicology and omics: A question of sex

Molecular initiating events and downstream transcriptional/proteomic responses provide valuable information for adverse outcome pathways, which can be used predict the effects of chemicals on physiological systems. There has been a paucity of research that addresses sex-specific expression profiling in toxicology and due to cost, time, and logistic considerations, sex as a variable has not been widely considered. In response to this deficiency, federal agencies in the United States, Canada, and Europe have highlighted the importance of including sex as a variable in scientific investigations. Using case studies from both aquatic and mammalian toxicology, we report that there can be less than ~20-25% consensus in how the transcriptome and proteome of each sex responds to chemicals. Chemicals that have been shown to elicit sex-specific responses in the transcriptome or proteome include pharmaceuticals, anti-fouling agents, anticorrosive agents, and fungicides, among others. Sex-specific responses in the transcriptome and proteome are not isolated to whole animals, as investigations demonstrate that primary cell cultures isolated from each sex responds differently to toxicants. This signifies that sex is important, even in cell lines. Sex has significant implications for predictive toxicology, and both male and female data are required to improve robustness of adverse outcome pathways.

BIOLOGICAL SIGNIFICANCE

Clinical toxicology recognizes that sex is an important variable, as pharmacokinetics (ADME; absorption, distribution, metabolism, and excretion) can differ between females and males. However, few studies in toxicology have explored the implication of sex in relation to the transcriptome and proteome of whole organisms. High-throughput molecular approaches are becoming more frequently applied in toxicity screens (e.g. pre-clinical experiments, fish embryos, cell lines, synthetic tissues) and such data are expected to build upon reporter-based cell assays (e.g. receptor activation, enzyme inhibition) used in toxicant screening programs (i.e. Tox21, ToxCast, REACH). Thus, computational models can more accurately predict the diversity of adverse effects that can occur from chemical exposure within the biological system. Our studies and those synthesized from the literature suggest that the transcriptome and proteome of females and males respond quite differentially to chemicals. This has significant implications for predicting adverse effects in one sex when using molecular data generated in the other sex. While molecular initiating events are not expected to differ dramatically between females and males (i.e. an estrogen binds estrogen receptors in both sexes), it is important to acknowledge that the downstream transcriptomic and proteomic responses can differ based upon the presence/absence of co-regulators and inherent sex-specific variability in regulation of transcriptional and translational machinery. Transcriptomic and proteomic studies also reveal that cell processes affected by chemicals can differ due to sex, and this can undoubtedly lead to sex-specific physiological responses.

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.

Dose-response analysis of epigenetic, metabolic, and apical endpoints after short-term exposure to experimental hepatotoxicants

Identification of sensitive and novel biomarkers or endpoints associated with toxicity and carcinogenesis is of a high priority. There is increasing interest in the incorporation of epigenetic and metabolic biomarkers to complement apical data; however, a number of questions, including the tissue specificity, dose-response patterns, early detection of those endpoints, and the added value need to be addressed. In this study, we investigated the dose-response relationship between apical, epigenetic, and metabolomics endpoints following short-term exposure to experimental hepatotoxicants, clofibrate (CF) and phenobarbital (PB). Male F344 rats were exposed to PB (0, 5, 25, and 100 mg/kg/day) or CF (0, 10, 50, and 250 mg/kg/day) for seven days. Exposure to PB or CF resulted in dose-dependent increases in relative liver weights, hepatocellular hypertrophy and proliferation, and increases in Cyp2b1 and Cyp4a1 transcripts. These changes were associated with altered histone modifications within the regulatory units of cytochrome genes, LINE-1 DNA hypomethylation, and altered microRNA profiles. Metabolomics data indicated alterations in the metabolism of bile acids. This study provides the first comprehensive analysis of the apical, epigenetic and metabolic alterations, and suggests that the latter two occur within or near the dose response curve of apical endpoint alterations following exposure to experimental hepatotoxicants.

Sodium bisulfite pyrosequencing revealed that developmental exposure to environmental contaminant mixtures does not affect DNA methylation of DNA repeats in Sprague-Dawley rats

Hypomethylation of DNA repeats has been linked to diseases and cancer predisposition. Human studies suggest that higher blood concentrations of environmental contaminants (EC) correlate with levels of hypomethylation of DNA repeats in blood. The objective of this study was to examine the effect of in utero and/or lactational exposure to EC on the methylation of DNA repeats (LINE-1 and identifier element) in Sprague-Dawley rat pups at birth, at postnatal day (PND) 21, and in adulthood (PND78-86). From gestation day 0 to PND20, dams were exposed to a mixture "M" of polychlorinated biphenyls (PCB), pesticides, and methylmercury (MeHg), at 0.5 or 1 mg/kg/d (0.5M and M). At birth, some control (C) and M litters were cross-fostered to create the following in utero/postnatal exposure groups: C/C, M/C, C/M, M/M. Additional dams received 1.8 ng/kg/d of a mixture of aryl-hydrocarbon receptor (AhR) agonists (non-ortho-PCB, PC-dibenzodioxins, and PC-dibenzofurans) without or with 0.5M (0.5MAhR). Measurements of EC residue levels confirmed differences in their accumulation across treatments, age, and tissues. Although induction of hepatic detoxification enzyme activities (cytochrome P-450) demonstrated biological effects of treatments, the assessment of methylation in DNA repeats by sodium bisulfite pyrosequencing of liver, spleen, and thymus samples revealed no marked treatment-related effects but significant tissue- and age-related methylation differences. Further studies are required to determine whether absence of significant observable treatment effects on methylation of DNA repeats in the rat relate to tissue, strain, or species differences.

Methyl isobutyl ketone-induced hepatocellular carcinogenesis in B6C3F1 mice: A constitutive androstane receptor (CAR)-mediated mode of action

In a National Toxicology Program (NTP) chronic inhalation study with methyl isobutyl ketone (MIBK), increases in hepatocellular adenomas and hepatocellular adenomas and carcinomas (combined) were observed in male and female B6C3F₁ mice at 1800 ppm. A DNA reactive Mode-of-Action (MOA) for this liver tumor response is not supported by the evidence as MIBK and its major metabolites lack genotoxicity in both in vitro and in vivo studies. Constitutive androstane receptor (CAR) nuclear receptor-mediated activation has been hypothesized as the MOA for MIBK-induced mouse liver tumorigenesis. To further investigate the MOA for MIBK-induced murine liver tumors, male and female B6C3F1, C57BL/6, and CAR/PXR Knockout (KO) mice were exposed to either 0 or 1800 ppm MIBK for 6 h/day, 5 days/week for a total of 10 days. On day 1, mice were implanted with osmotic mini-pumps containing 5-Bromo-2-deoxyuridine (BrdU) 1 h following exposure and humanely euthanized 1-3 h following the final exposure. B6C3F₁ and C57BL/6 mice had statistically significant increases in liver weights compared to controls that corresponded with hepatocellular hypertrophy and increased mitotic figures. Hepatocellular proliferation data indicated induction of S-phase DNA synthesis in B6C3F₁ and C57BL/6 mice exposed to 1800 ppm MIBK compared to control, and no increase was observed in MIBK exposed CAR/PXR KO mice. Liver gene expression changes indicated a maximally-induced Cyp2b10 (CAR-associated) transcript and a slight increase in Cyp3a11(PXR-associated) transcript in B6C3F₁ and C57BL/6 mice exposed to 1800 ppm MIBK compared to controls, but not in Cyp1a1 (AhR-associated) or Cyp4a10 (PPAR-α-associated) transcripts. CAR/PXR KO mice exposed to 1800 ppm MIBK showed no evidence of activation of AhR, CAR, PXR or PPAR-α nuclear receptors via their associated transcripts. MIBK induced hepatic effects are consistent with a phenobarbital-like MOA where the initiating events are activation of the CAR and PXR nuclear receptors and resultant hepatocellular proliferation leading to rodent liver tumors.

Editor's Highlight: Dose-Response Analysis of RNA-Seq Profiles in Archival Formalin-Fixed Paraffin-Embedded Samples

Use of archival resources has been limited to date by inconsistent methods for genomic profiling of degraded RNA from formalin-fixed paraffin-embedded (FFPE) samples. RNA-sequencing offers a promising way to address this problem. Here, we evaluated transcriptomic dose responses using RNA-sequencing in paired FFPE and frozen (FROZ) samples from 2 archival studies in mice, one <2 years old and the other >20 years old. Experimental treatments included 3 different doses of di(2-ethylhexyl)phthalate or dichloroacetic acid for the recently archived and older studies, respectively. Total RNA was ribo-depleted and sequenced using the Illumina HiSeq platform. In the recently archived study, FFPE samples had 35% lower total counts compared to FROZ samples but high concordance in fold-change values of differentially expressed genes (DEGs) (r²  =0.99), highly enriched pathways (90% overlap with FROZ), and benchmark dose estimates for preselected target genes (<5% difference vs FROZ). In contrast, older FFPE samples had markedly lower total counts (3% of FROZ) and poor concordance in global DEGs and pathways. However, counts from FFPE and FROZ samples still positively correlated (r² = 0.84 across all transcripts) and showed comparable dose responses for more highly expressed target genes. These findings highlight potential applications and issues in using RNA-sequencing data from FFPE samples. Recently archived FFPE samples were highly similar to FROZ samples in sequencing quality metrics, DEG profiles, and dose-response parameters, while further methods development is needed for older lower-quality FFPE samples. This work should help advance the use of archival resources in chemical safety and translational science.

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.

Dose and Effect Thresholds for Early Key Events in a PPARα-Mediated Mode of Action

Current strategies for predicting adverse health outcomes of environmental chemicals are centered on early key events in toxicity pathways. However, quantitative relationships between early molecular changes in a given pathway and later health effects are often poorly defined. The goal of this study was to evaluate short-term key event indicators using qualitative and quantitative methods in an established pathway of mouse liver tumorigenesis mediated by peroxisome proliferator-activated receptor alpha (PPARα). Male B6C3F1 mice were exposed for 7 days to di (2-ethylhexyl) phthalate (DEHP), di-n-octyl phthalate (DNOP), and n-butyl benzyl phthalate (BBP), which vary in PPARα activity and liver tumorigenicity. Each phthalate increased expression of select PPARα target genes at 7 days, while only DEHP significantly increased liver cell proliferation labeling index (LI). Transcriptional benchmark dose (BMDT) estimates for dose-related genomic markers stratified phthalates according to hypothetical tumorigenic potencies, unlike BMDs for non-genomic endpoints (relative liver weights or proliferation). The 7-day BMDT values for Acot1 as a surrogate measure for PPARα activation were 29, 370, and 676 mg/kg/day for DEHP, DNOP, and BBP, respectively, distinguishing DEHP (liver tumor BMD of 35 mg/kg/day) from non-tumorigenic DNOP and BBP. Effect thresholds were generated using linear regression of DEHP effects at 7 days and 2-year tumor incidence values to anchor early response molecular indicators and a later phenotypic outcome. Thresholds varied widely by marker, from 2-fold (Pdk4 and proliferation LI) to 30-fold (Acot1) induction to reach hypothetical tumorigenic expression levels. These findings highlight key issues in defining thresholds for biological adversity based on molecular changes.

Assessing the carcinogenic potential of low-dose exposures to chemical mixtures in the environment: the challenge ahead

Low-dose exposures to common environmental chemicals that are deemed safe individually may be combining to instigate carcinogenesis, thereby contributing to the incidence of cancer. This risk may be overlooked by current regulatory practices and needs to be vigorously investigated.

Lifestyle factors are responsible for a considerable portion of cancer incidence worldwide, but credible estimates from the World Health Organization and the International Agency for Research on Cancer (IARC) suggest that the fraction of cancers attributable to toxic environmental exposures is between 7% and 19%. To explore the hypothesis that low-dose exposures to mixtures of chemicals in the environment may be combining to contribute to environmental carcinogenesis, we reviewed 11 hallmark phenotypes of cancer, multiple priority target sites for disruption in each area and prototypical chemical disruptors for all targets, this included dose-response characterizations, evidence of low-dose effects and cross-hallmark effects for all targets and chemicals. In total, 85 examples of chemicals were reviewed for actions on key pathways/mechanisms related to carcinogenesis. Only 15% (13/85) were found to have evidence of a dose-response threshold, whereas 59% (50/85) exerted low-dose effects. No dose-response information was found for the remaining 26% (22/85). Our analysis suggests that the cumulative effects of individual (non-carcinogenic) chemicals acting on different pathways, and a variety of related systems, organs, tissues and cells could plausibly conspire to produce carcinogenic synergies. Additional basic research on carcinogenesis and research focused on low-dose effects of chemical mixtures needs to be rigorously pursued before the merits of this hypothesis can be further advanced. However, the structure of the World Health Organization International Programme on Chemical Safety ‘Mode of Action’ framework should be revisited as it has inherent weaknesses that are not fully aligned with our current understanding of cancer biology.

Mode of action analysis for pesticide-induced rodent liver tumours involving activation of the constitutive androstane receptor: relevance to human cancer risk

A number of non-genotoxic chemicals, including some pesticides, have been shown to increase the incidence of liver tumours in rats and/or mice. Frameworks for analysing the modes of action (MOAs) by which chemicals produce liver tumours in rodents and the relevance of such tumour data for human risk assessment have now been established. One common MOA for rodent liver tumour formation by non-genotoxic chemicals involves activation of the constitutive androstane receptor (CAR). Key and associative events for a CAR-activation MOA include receptor activation, liver hypertrophy, induction of cytochrome P450 enzyme activities, increased replicative DNA synthesis, altered hepatic foci and liver tumours. While some effects of rodent CAR activators can be observed in human liver, a major species difference is that, unlike rodents, CAR activators do not increase replicative DNA synthesis in human hepatocytes. The CAR-activation MOA for rodent liver tumour formation is thus not plausible for humans, and hence such compounds do not pose a hepatocarcinogenic hazard for humans.

Risk assessments for chronic exposure of children and prospective parents to ethylbenzene (CAS No. 100-41-4)

Potential chronic health risks for children and prospective parents exposed to ethylbenzene were evaluated in response to the Voluntary Children's Chemical Evaluation Program. Ethylbenzene exposure was found to be predominately via inhalation with recent data demonstrating continuing decreases in releases and both outdoor and indoor concentrations over the past several decades. The proportion of ethylbenzene in ambient air that is attributable to the ethylbenzene/styrene chain of commerce appears to be relatively very small, less than 0.1% based on recent relative emission estimates. Toxicity reference values were derived from the available data, with physiologically based pharmacokinetic models and benchmark dose methods used to assess dose-response relationships. An inhalation non-cancer reference concentration or RfC of 0.3 parts per million (ppm) was derived based on ototoxicity. Similarly, an oral non-cancer reference dose or RfD of 0.5 mg/kg body weight/day was derived based on liver effects. For the cancer assessment, emphasis was placed upon mode of action information. Three of four rodent tumor types were determined not to be relevant to human health. A cancer reference value of 0.48 ppm was derived based on mouse lung tumors. The risk characterization for ethylbenzene indicated that even the most highly exposed children and prospective parents are not at risk for non-cancer or cancer effects of ethylbenzene.

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

The nuclear receptor family member constitutive activated receptor (CAR) is activated by structurally diverse drugs and environmentally-relevant chemicals leading to transcriptional regulation of genes involved in xenobiotic metabolism and transport. Chronic activation of CAR increases liver cancer incidence in rodents, whereas suppression of CAR can lead to steatosis and insulin insensitivity. Here, analytical methods were developed to screen for chemical treatments in a gene expression compendium that lead to alteration of CAR activity. A gene expression biomarker signature of 83 CAR-dependent genes was identified using microarray profiles from the livers of wild-type and CAR-null mice after exposure to three structurally-diverse CAR activators (CITCO, phenobarbital, TCPOBOP). A rank-based algorithm (Running Fisher’s algorithm (p-value ≤ 10^-4)) was used to evaluate the similarity between the CAR biomarker signature and a test set of 28 and 32 comparisons positive or negative, respectively, for CAR activation; the test resulted in a balanced accuracy of 97%. The biomarker signature was used to identify chemicals that activate or suppress CAR in an annotated mouse liver/primary hepatocyte gene expression database of ~1850 comparisons. CAR was activated by 1) activators of the aryl hydrocarbon receptor (AhR) in wild-type but not AhR-null mice, 2) pregnane X receptor (PXR) activators in wild-type and to lesser extents in PXR-null mice, and 3) activators of PPARα in wild-type and PPARα-null mice. CAR was consistently activated by five conazole fungicides and four perfluorinated compounds. Comparison of effects in wild-type and CAR-null mice showed that the fungicide propiconazole increased liver weight and hepatocyte proliferation in a CAR-dependent manner, whereas the perfluorinated compound perfluorooctanoic acid (PFOA) increased these endpoints in a CAR-independent manner. A number of compounds suppressed CAR coincident with increases in markers of inflammation including acetaminophen, concanavalin A, lipopolysaccharide, and 300 nm silica particles. In conclusion, we have shown that a CAR biomarker signature coupled with a rank-based similarity method accurately predicts CAR activation. This analytical approach, when applied to a gene expression compendium, increased the universe of known chemicals that directly or indirectly activate CAR, highlighting the promiscuous nature of CAR activation and signaling through activation of other xenobiotic-activated receptors.

Latent carcinogenicity of early-life exposure to dichloroacetic acid in mice

Environmental exposures occurring early in life may have an important influence on cancer risk later in life. Here, we investigated carryover effects of dichloroacetic acid (DCA), a small molecule analog of pyruvate with metabolic programming properties, on age-related incidence of liver cancer. The study followed a stop-exposure/promotion design in which 4-week-old male and female B6C3F1 mice received the following treatments: deionized water alone (dH2O, control); dH2O with 0.06% phenobarbital (PB), a mouse liver tumor promoter; or DCA (1.0, 2.0 or 3.5g/l) for 10 weeks followed by dH2O or PB (n = 20-30/group/sex). Pathology and molecular assessments were performed at 98 weeks of age. In the absence of PB, early-life exposure to DCA increased the incidence and number of hepatocellular tumors in male and female mice compared with controls. Significant dose trends were observed in both sexes. At the high dose level, 10 weeks of prior DCA treatment induced comparable effects (≥85% tumor incidence and number) to those seen after continuous lifetime exposure. Prior DCA treatment did not enhance or inhibit the carcinogenic effects of PB, induce persistent liver cytotoxicity or preneoplastic changes on histopathology or alter DNA sequence variant profiles within liver tumors compared with controls. Distinct changes in liver messenger RNA and micro RNA profiles associated with prior DCA treatment were not apparent at 98 weeks. Our findings demonstrate that early-life exposure to DCA may be as carcinogenic as life-long exposures, potentially via epigenetic-mediated effects related to cellular metabolism.

Phenobarbital-mediated tumor promotion in transgenic mice with humanized CAR and PXR

The nuclear receptors CAR (constitutive androstane receptor) and possibly PXR (pregnane X receptor) mediate the hepatic effects of phenobarbital (PB) and similar-acting compounds. Although PB is a potent nongenotoxic tumor promoter in rodent liver, epidemiological data from epilepsy patients treated with phenobarbital do not show a specific role of PB in human liver cancer risk. That points to species differences in the susceptibility to tumor promotion by PB, which might be attributed to divergent functions of the PB receptors CAR and PXR in mice and humans. In the present study, male transgenic mice expressing human CAR and PXR were used to detect possible differences between wild-type (WT) and humanized mice in their response to CAR activation in a tumor initiation/promotion experiment with a single injection of the tumor initiator N-nitrosodiethylamine preceding chronic PB treatment for 10 months. Analysis of liver tumor burden revealed that PB strongly promoted the outgrowth of hepatocellular adenoma driven by activated β-catenin in WT mice, whereas the tumor-promoting effect of PB was much less pronounced in the humanized group. In conclusion, the present findings demonstrate that human CAR and PXR support tumor promotion by PB in mouse liver, but to a significantly lesser extent than the WT murine receptors.