Susan D. Hester

Dose-Response Modeling of Early Molecular and Cellular Key Events in the CAR-Mediated Hepatocarcinogenesis Pathway

Low-dose extrapolation and dose-related transitions are paramount in the ongoing debate regarding the quantification of cancer risks for nongenotoxic carcinogens. Phenobarbital (PB) is a prototypical nongenotoxic carcinogen that activates the constitutive androstane receptor (CAR) resulting in rodent liver tumors. In this study, male and female CD-1 mice administered dietary PB at 0, 0.15, 1.5, 15, 75, or 150 mg/kg-day for 2 or 7 days to characterize multiple apical and molecular endpoints below, at (∼75 mg/kg-day), and above the carcinogenic dose level of PB and examine these responses using benchmark dose modeling. Linear toxicokinetics were observed for all doses. Increased liver weight, hepatocellular hypertrophy, and mitotic figures were seen at 75 and 150 mg/kg-day. CAR activation, based on Cyp2b qPCR and pentoxyresorufin dealkylase activity, occurred at doses ≥ 1.5 mg/kg-day. The no-observable transcriptional effect level for global gene expression was 15 mg/kg-day. At 2 days, several xenobiotic metabolism and cell protective pathways were activated at lower doses and to a greater degree in females. However, hepatocellular proliferation, quantified by bromodeoxyuridine immunohistochemistry, was the most sensitive indicator of PB exposure with female mice more sensitive than males, contrary to sex-specific differences in sensitivity to hepatocarcinogenesis. Taken together, the identification of low-dose cellular and molecular transitions in the subtumorigenic dose range aids the understanding of early key events in CAR-mediated hepatocarcinogenesis.

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

BackgroundThe use of gene expression profiling in both clinical and laboratory settings would be enhanced by better characterization of variance due to individual, environmental, and technical factors. Meta-analysis of microarray data from untreated or vehicle-treated animals within the control arm of toxicogenomics studies could yield useful information on baseline fluctuations in gene expression, although control animal data has not been available on a scale and in a form best served for data-mining.ResultsA dataset of control animal microarray expression data was assembled by a working group of the Health and Environmental Sciences Institutes Technical Committee on the Application of Genomics in Mechanism Based Risk Assessment in order to provide a public resource for assessments of variability in baseline gene expression. Data from over 500 Affymetrix microarrays from control rat liver and kidney were collected from 16 different institutions. Thirty-five biological and technical factors were obtained for each animal, describing a wide range of study characteristics, and a subset were evaluated in detail for their contribution to total variability using multivariate statistical and graphical techniques.ConclusionThe study factors that emerged as key sources of variability included gender, organ section, strain, and fasting state. These and other study factors were identified as key descriptors that should be included in the minimal information about a toxicogenomics study needed for interpretation of results by an independent source. Genes that are the most and least variable, gender-selective, or altered by fasting were also identified and functionally categorized. Better characterization of gene expression variability in control animals will aid in the design of toxicogenomics studies and in the interpretation of their results.

Follow-Up Actions from Positive Results of In Vitro Genetic Toxicity Testing

Appropriate follow‐up actions and decisions are needed when evaluating and interpreting clear positive results obtained in the in vitro assays used in the initial genotoxicity screening battery (i.e., the battery of tests generally required by regulatory authorities) to assist in overall risk‐based decision making concerning the potential effects of human exposure to the agent under test. Over the past few years, the International Life Sciences Institute (ILSI) Health and Environmental Sciences Institute (HESI) Project Committee on the Relevance and Follow‐up of Positive Results in In Vitro Genetic Toxicity (IVGT) Testing developed a decision process flow chart to be applied in case of clear positive results in vitro. It provides for a variety of different possibilities and allows flexibility in choosing follow‐up action(s), depending on the results obtained in the initial battery of assays and available information. The intent of the Review Subgroup was not to provide a prescriptive testing strategy, but rather to reinforce the concept of weighing the totality of the evidence. The Review Subgroup of the IVGT committee highlighted the importance of properly analyzing the existing data, and considering potential confounding factors (e.g., possible interactions with the test systems, presence of impurities, irrelevant metabolism), and chemical modes of action when analyzing and interpreting positive results in the in vitro genotoxicity assays and determining appropriate follow‐up testing. The Review Subgroup also examined the characteristics, strengths, and limitations of each of the existing in vitro and in vivo genotoxicity assays to determine their usefulness in any follow‐up testing. Environ. Mol. Mutagen., 2011.

Discrimination of Tumorigenic Triazole Conazoles from Phenobarbital by Transcriptional Analyses of Mouse Liver Gene Expression

Conazoles are fungicides used to control fungal growth in environmental settings and to treat humans with fungal infections. Mouse hepatotumorigenic conazoles display many of the same hepatic toxicologic responses as the mouse liver carcinogen phenobarbital (PB): constitutive androstane receptor (CAR) activation, hypertrophy, Cyp2b induction, and increased cell proliferation. The goal of this study was to apply transcriptional analyses to hepatic tissues from mice exposed to PB, propiconazole (Pro) or triadimefon (Tri) at tumorigenic exposure levels to reveal similarities and differences in response among these treatments. Mice were administered diets containing PB (850 ppm), Pro (2500 ppm), or Tri (1800 ppm) for 4 and 30 days. Targeted transcriptomic analyses were conducted at the gene level examining differentially expressed genes (DEGs), and subsets of DEGs: cell cycle genes, and transcription factors. Analyses were also conducted on function, pathway and network levels examining Ingenuity Pathway Analysis Tox Lists and Canonical Pathways, and Gene-Go MetaCore dynamic networks and their central hubs. Genes expressed by PB or the two conazoles were also compared with those genes associated with human hepatocellular cancer. The results from these analyses indicated greater differences between PB and the two conazoles than similarities. Significant commonalities between the two conazole treatments were also noted. We posit that the transcriptional profiles of tissues exposed to toxic chemicals inherently contain their mechanisms of toxicity. We conclude that although PB and these 2 conazoles induce mouse liver tumors and exhibit similar toxicological responses, their transcriptional profiles are significantly different and thus their mechanisms of tumorigenic action are likely to differ.

The Hepatocarcinogenic Conazoles: Cyproconazole, Epoxiconazole, and Propiconazole Induce a Common Set of Toxicological and Transcriptional Responses

Conazoles are fungicides used as agricultural pesticides and pharmaceutical products. We investigated whether a common core of toxicological and transcriptional responses underlies the observed carcinogenic effects of three conazoles: cyproconazole, epoxiconazole, and propiconazole. In studies where mice were fed diets of these conazoles for 30 days, we found a common set of toxicological effects altered by these conazoles: hepatomegaly, hepatocellular hypertrophy, decreased serum cholesterol, decreased hepatic levels of all-trans-retinoic acid, and increased hepatic cell proliferation. Microarray-based transcriptional analysis revealed 330 significantly altered probe sets common to these conazoles, many of which showed strong dose responses for cytochrome P450, glutathione S-transferase, and oxidative stress genes. More detailed analyses identified a subset of 80 altered genes common to the three conazoles that were associated with cancer. Pathways associated with these genes included xenobiotic metabolism, oxidative stress, cell signaling, and cell proliferation. A common TGFα-centric pathway was identified within the 80-gene set, which, in combination with the toxicological and other transcriptomic findings, provides a more refined toxicity profile for these carcinogenic conazoles.

The Cellular and Genomic Response of an Immortalized Microglia Cell Line (BV2) to Concentrated Ambient Particulate Matter

Ambient particulate matter (PM) damages pulmonary tissue through oxidative stress (OS) pathways. Several reports indicate that the brain is another affected target of PM exposure. Since microglia (brain macrophages) are critical to OS-mediated neurodegeneration, the cellular and genomic response of immortalized mouse microglia (BV2) was examined in response to fine (≤ 2.5 μ m) concentrated ambient particles (CAPs) collected from Tuxedo, NY. Samples of CAPs were labeled as high potency (HP) or low potency (LP) depending on their stimulation of nuclear factor (NF)-κ B activity in human bronchial epithelial cells. Compositional analysis of these samples, performed during their original collection, indicated a strong correlation between HP CAPs and and the presence of nickel and vanadium (). Exposure of the BV2 microglia to LP CAPs reduced intracellular levels of ATP (≥ 250 μ g/ml) and depolarized mitochondrial membranes (≥ 6 μ g/ml) within 15 min of exposure. HP and LP CAPs (≥ 25 μ g/ml) differentially affected the endogenous scavengers, glutathione and nonprotein sulfhydryl in BV2 microglia after 1.5 h of exposure. Both HP and LP CAPs stimulated the release of proinflammatory cytokines tumor necrosis factor (TNF) α and interleukin (IL)-6 after 6 h of exposures. Microarray analysis of BV2 microglia exposed to either HP or LP CAPs (75 μ g/ml, 4 h) identified 3200 (HP CAPs) and 160 (LP CAPs) differentially expressed (up- and downregulated) genes relative to media controls. Of the 3200 genes significantly affected by HP CAPs, the most prominent upregulated gene probes related to inflammatory pathways associated with Toll-like receptor signaling, MAPK signaling, T- and B-cell receptor signaling, apoptosis, and various proinflammatory cytokines and their receptors. LP CAPs significantly affected 160 genes that related to pathways associated with cellular maintenance and division, cell cycling and nuclear events. These data suggest that HP CAPs, which contained higher levels of nickel and vanadium than LP CAPs, appear to be more inflammatory and selectively upregulated the expression of inflammatory and innate immunity pathways in BV2 microglia.

Increased Transcription of Immune and Metabolic Pathways in Naïve and Allergic Mice Exposed to Diesel Exhaust

Diesel exhaust (DE) has been shown to enhance allergic sensitization in animals following high-dose instillation or chronic inhalation exposure scenarios. The purpose of this study was to determine if short-term exposures to diluted DE enhance allergic immune responses to antigen, and identify possible mechanisms using microarray technology. BALB/c mice were exposed to filtered air or diluted DE to yield particle concentrations of 500 or 2000 mug/m(3) 4 h/day on days 0-4. Mice were immunized intranasally with ovalbumin (OVA) antigen or saline on days 0-2, challenged on day 18 with OVA or saline, and all mice were challenged with OVA on day 28. Mice were necropsied either 4 h after the last DE exposure on day 4, or 18, 48, and 96 h after the last challenge. Immunological endpoints included OVA-specific serum IgE, biochemical and cellular profiles of bronchoalveolar lavage (BAL), and cytokine production in the BAL. OVA-immunized mice exposed to both concentrations of DE had increased eosinophils, neutrophils, lymphocytes, and interleukin-6 (high dose only) post-challenge compared with OVA control, whereas DE/saline exposure yielded increases in neutrophils at the high dose only. Transcriptional microarray analysis 4 h after the last DE exposure demonstrated distinct gene expression profiles for the high-dose DE/OVA and DE/saline groups. DE/OVA induced oxidative stress and metabolism pathways, whereas DE in the absence of immunization modulated cell cycle control, growth and differentiation, G-proteins, and cell adhesion pathways. This study shows for the first time early changes in gene expression induced by the combination of DE inhalation and mucosal immunization, which resulted in stronger development of allergic eosinophilia.

Gene expression of normal human epidermal keratinocytes modulated by trivalent arsenicals

Chronic exposure to inorganic arsenic (iAs) is associated with the development of benign and malignant human skin lesions including nonmelanoma skin cancers. The precise arsenical form(s) responsible for this carcinogenic effect are unknown, although trivalent inorganic arsenic (iAsIII) and two of its toxic metabolites, monomethylarsonous acid (MMAIII) and methylarsinous acid (DMAIII), are attractive candidates. In an effort to better understand and compare their toxic effects in the skin, we compared the global gene expression profiles of normal human epidermal keratinocytes (NHEKs) exposed to varying noncytotoxic/slightly cytotoxic concentrations of iAsIII, MMAIII, and DMAIII for 24u2009h. Exposure to each arsenical treatment group exhibited a dose effect in the number of altered genes and the magnitude of expression change in NHEKs. The most significant gene expression changes associated with iAsIII and MMAIII exposure were consistent with several key events believed to be important to As‐driven skin carcinogenesis, namely induction of oxidative stress, increased transcript levels of keratinocyte growth factors, and modulation of MAPK and NF‐κB pathways. At both comparable arsenical concentrations and comparable NHEK toxicity, greater potential carcinogenic effects were observed in MMAIII‐exposed NHEKs than those exposed to iAsIII, including involvement of more proinflammatory signals and increased transcript levels of more growth factor genes. In contrast, none of these above‐mentioned transcriptional trends were among the most significantly altered functions in the DMAIII treatment group. This study suggests the relative capacity of each of the tested arsenicals to drive suspected key events in As‐mediated skin carcinogenesis is MMAIIIu2009>u2009iAsIII with little contribution from DMAIII.

Concordance of Transcriptional and Apical Benchmark Dose Levels for Conazole-Induced Liver Effects in Mice

The ability to anchor chemical class-based gene expression changes to phenotypic lesions and to describe these changes as a function of dose and time informs mode-of-action determinations and improves quantitative risk assessments. Previous global expression profiling identified a 330-probe cluster differentially expressed and commonly responsive to 3 hepatotumorigenic conazoles (cyproconazole, epoxiconazole, and propiconazole) at 30 days. Extended to 2 more conazoles (triadimefon and myclobutanil), the present assessment encompasses 4 tumorigenic and 1 nontumorigenic conazole. Transcriptional benchmark dose levels (BMDL(T)) were estimated for a subset of the cluster with dose-responsive behavior and a ≥ 5-fold increase or decrease in signal intensity at the highest dose. These genes primarily encompassed CAR/RXR activation, P450 metabolism, liver hypertrophy- glutathione depletion, LPS/IL-1-mediated inhibition of RXR, and NRF2-mediated oxidative stress pathways. Median BMDL(T) estimates from the subset were concordant (within a factor of 2.4) with apical benchmark doses (BMDL(A)) for increased liver weight at 30 days for the 5 conazoles. The 30-day median BMDL(T) estimates were within one-half order of magnitude of the chronic BMDLA for hepatocellular tumors. Potency differences seen in the dose-responsive transcription of certain phase II metabolism, bile acid detoxification, and lipid oxidation genes mirrored each conazoles tumorigenic potency. The 30-day BMDL(T) corresponded to tumorigenic potency on a milligram per kilogram day basis with cyproconazole > epoxiconazole > propiconazole > triadimefon > myclobutanil (nontumorigenic). These results support the utility of measuring short-term gene expression changes to inform quantitative risk assessments from long-term exposures.

Transcriptional Profile of Diuron-Induced Toxicity on the Urinary Bladder of Male Wistar Rats to Inform Mode of Action

Diuron (3-(3,4-dichlorophenyl)-1,1-dimethylurea) is a substituted urea herbicide that induces rat urinary bladder urothelial tumors at high dietary levels (2500 ppm). The specific mode of action and molecular alterations triggered by diuron, however, have not been clarified. The present study evaluated the dose-dependent effects of mucosal alterations and transcriptional changes in the urinary bladder of rats exposed to diuron. Six-week-old male Wistar rats were treated with 0, 60, 125, 1250, and 2500 ppm of diuron in the diet for 20 weeks. Histologic examination showed urothelial hyperplasia present in rats treated with either 1250 or 2500 ppm of diuron but not 60 or 125 ppm. Comprehensive gene expression analyses of urothelial cell RNA were conducted using Affymetrix microarrays. The numbers of differentially expressed transcripts between each treatment group and control increased with diuron dose. Based on similar histology and gene expression responses, the treatment groups were regrouped into a high-dose (1250 and 2500 ppm) and low-dose group (60 and 125 ppm). These data suggest that persistent exposure to high dietary concentrations of diuron induces oxidative stress, increases cellular metabolism, and enhances cell death that is associated with sustained urothelial hyperplasia.