Transcriptomic re-analyses of human hepatocyte spheroids treated with PFAS reveals chain length and dose-dependent modes of action.

To identify pathway perturbations and examine biological modes of action (MOAs) for various perfluoroalkyl substances, we re-analyzed published in vitro gene expression studies from human primary liver spheroids. With treatment times ranging from 10 to 14 days, shorter-chain PFAS (those with 6 or fewer fluorinated carbon atoms in the alkyl chain) showed enrichment for pathways of fatty acid metabolism and fatty acid beta-oxidation with upregulated genes. Longer-chain PFAS compounds, specifically PFOS (perfluorooctane sulfonate), PFDS (perfluorodecane sulfonate), and higher doses of PFOA (perfluorooctanoic acid), had enrichment for pathways involved in steroid metabolism, fatty acid metabolism, and biological oxidation for downregulated genes. Although PFNA (perfluorononanoic acid), PFDA (perfluorodecanoic acid), and PFUnDA (perfluoroundecanoic acid) were more toxic and could only be examined after a 1-day treatment, all three had enrichment patterns similar to those observed with PFOS. With PFOA there were dose-dependent changes in pathway enrichment, shifting from upregulation of fatty acid metabolism and downregulation of steroid metabolism to downregulation of both at higher doses. The response to PFHpS (perfluoroheptanesulfonic acid) was similar to the PFOA pattern at the lower treatment dose. Based on results of transcription factor binding sites analyses, we propose that downregulation of pathways of lipid metabolism by longer chain PFAS may be due to inhibitory interactions of PPARD on genes controlled by PPARA and PPARG. In conclusion, our transcriptomic analysis indicates that the biological MOAs of PFAS compounds differ according to chain length and dose, and that risk assessments for PFAS should consider these differences in biological MOAs when evaluating mixtures of these compounds.

Transcriptomic analysis of AHR wildtype and Knock-out rat livers supports TCDD's role in AHR/ARNT-mediated circadian disruption and hepatotoxicity.

Single, high doses of TCDD in rats are known to cause wasting, a progressive loss of 30 to 50% body weight and death within several weeks. To identify pathway perturbations at or near doses causing wasting, we examined differentially gene expression (DGE) and pathway enrichment in centrilobular (CL) and periportal (PP) regions of female rat livers following 6 dose levels of TCDD - 0, 3, 22, 100, 300, and 1000 ng/kg/day, 5 days/week for 4 weeks. At the higher doses, rats lost weight, had increased liver/body weight ratios and nearly complete cessation of liver cell proliferation, signs consistent with wasting. DGE curves were left shifted for the CL versus the PP regions. Canonical Phase I and Phase II genes were maximally increased at lower doses and remained elevated at all doses. At lower doses, ≤ 22 ng/kg/day in the CL and ≤ 100 ng/kg/day, upregulated genes showed transcription factor (TF) enrichment for AHR and ARNT. At the mid- and high-dose doses, there was a large number of downregulated genes and pathway enrichment for DEGs which showed downregulation of many cellular metabolism processes including those for steroids, fatty acid metabolism, pyruvate metabolism and citric acid cycle. There was significant TF enrichment of the hi-dose downregulated genes for RXR, ESR1, LXR, PPARalpha. At the highest dose, there was also pathway enrichment with upregulated genes for extracellular matrix organization, collagen formation, hemostasis and innate immune system. TCDD demonstrates most of its effects through binding the aryl hydrocarbon receptor (AHR) while the downregulation of metabolism genes at higher TCDD doses is known to be independent of AHR binding to DREs. Based on our results with DEG, we provide a hypothesis for wasting in which high doses of TCDD shift circadian processes away from the resting state, leading to greatly reduced synthesis of steroids and complex lipids needed for cell growth, and producing gene expression signals consistent with an epithelial-to-mesenchymal transition in hepatocytes.

Disinfection Byproducts of Haloacetaldehydes Disrupt Hepatic Lipid Metabolism and Induce Lipotoxicity in High-Fat Culture Conditions.

Unhealthy lifestyles, obesity, and environmental pollutants are strongly correlated with the development of nonalcoholic fatty liver disease (NAFLD). Haloacetaldehyde-associated disinfection byproducts (HAL-DBPs) at various multiples of concentrations found in finished drinking water together with high-fat (HF) were examined to gauge their mixed effects on hepatic lipid metabolism. Using new alternative methods (NAMs), studying effects in human cells in vitro for risk assessment, we investigated the combined effects of HF and HAL-DBPs on hepatic lipid metabolism and lipotoxicity in immortalized LO-2 human hepatocytes. Coexposure of HAL-DBPs at various multiples of environmental exposure levels with HF increased the levels of triglycerides, interfered with de novo lipogenesis, enhanced fatty acid oxidation, and inhibited the secretion of very low-density lipoproteins. Lipid accumulation caused by the coexposure of HAL-DBPs and HF also resulted in more severe lipotoxicity in these cells. Our results using an in vitro NAM-based method provide novel insights into metabolic reprogramming in hepatocytes due to coexposure of HF and HAL-DBPs and strongly suggest that the risk of NAFLD in sensitive populations due to HAL-DBPs and poor lifestyle deserves further investigation both with laboratory and epidemiological tools. We also discuss how results from our studies could be used in health risk assessments for HAL-DBPs.

Comparison of Rat Hepatocyte 2D-Monocultures and Hepatocytes Non-Parenchymal Cell Co-Cultures for Assessing Chemical Toxicity.

Liver responses are the most common endpoints used as the basis for setting exposure standards. Liver hepatocytes play a vital role in biotransformation of xenobiotics, but non-parenchymal cells (NPCs) in the liver are also involved in certain liver responses. Development of in vitro systems that more faithfully capture liver responses to reduce reliance on animals is a major focus of New Approach Methodology (NAMs). Since rodent regulatory studies are frequently the sole source safety assessment data, mode-of-action data, and used for risk assessments, in vitro rodent models that reflect in vivo responses need to be developed to reduce reliance on animal models. In the work presented in this paper, we developed a 2-D hepatocyte monoculture and 2-D liver cell co-culture system using rat liver cells. These models were assessed for conditions for short-term stability of the cultures and phenotypic and transcriptomic responses of 2 prototypic hepatotoxicants compounds - acetaminophen and phenobarbital. The optimized multi-cellular 2-D culture required use of freshly prepared hepatocytes and NPCs from a single rat, a 3:1 ratio of hepatocytes to NPCs and growth medium using 50% Complete Williams E medium (WEM) and 50% Endothelial Cell Medium (ECM). The transcriptomic responses of the 2 model systems to PB were compared to previous studies from TG-Gates on the gene expression changes in intact rats and the co-culture model responses were more representative of the in vivo responses. Transcriptomic read-outs promise to move beyond conventional phenotypic evaluations with these in vitro NAMs and provide insights about modes of action.

NRF2-ARE signaling is responsive to haloacetonitrile-induced oxidative stress in human keratinocytes.

Humans are exposed to disinfection by-products through oral, inhalation, and dermal routes, during bathing and swimming, potentially causing skin lesions, asthma, and bladder cancer. Nuclear factor E2-related factor 2 (NRF2) is a master regulator of the adaptive antioxidant response via the antioxidant reaction elements (ARE) orchestrating the transcription of a large group of antioxidant and detoxification genes. Here we used an immortalized human keratinocyte model HaCaT cells to investigate NRF2-ARE as a responder and protector in the acute cytotoxicity of seven haloacetonitriles (HANs), including chloroacetonitrile (CAN), bromoacetonitrile (BAN), iodoacetonitrile (IAN), bromochloroacetonitrile (BCAN), dichloroacetonitrile (DCAN), dibromoacetonitrile (DBAN), and trichloroacetonitrile (TCAN) found in drinking water and swimming pools. The rank order of cytotoxicity among the HANs tested was IAN ≈ BAN Ëƒ DBAN Ëƒ BCAN ˃ CAN Ëƒ TCAN Ëƒ DCAN based on their LC50. The HANs induced intracellular reactive oxygen species accumulation and activated cellular antioxidant responses in concentration- and time-dependent fashions, showing elevated NRF2 protein levels and ARE activity, induction of antioxidant genes, and increased glutathione levels. Additionally, knockdown of NRF2 by lentiviral shRNAs sensitized the HaCaT cells to HANs-induced cytotoxicity, emphasizing a protective role of NRF2 against the cytotoxicity of HANs. These results indicate that HANs cause oxidative stress and activate NRF2-ARE-mediated antioxidant response, which in turn protects the cells from HANs-induced cytotoxicity, highlighting that NRF2-ARE activity could be a sensitive indicator to identify and characterize the oxidative stress induced by HANs and other environmental pollutants.

Integration of paraquat pharmacokinetic data across species using PBPK modelling.

Paraquat dichloride (PQ) is a non-selective herbicide which has been the subject of numerous toxicology studies over more than 50 years. This paper describes the development of a physiologically-based pharmacokinetic (PBPK) model of PQ kinetics for the rat, mouse and dog, firstly to aid the interpretation of studies in which no kinetic measurements were made, and secondly to enable the future extension of the model to humans. Existing pharmacokinetic data were used to develop a model for the rat and mouse. Simulations with this preliminary model were then used to identify key data gaps and to design a new blood binding study to reduce uncertainty in critical aspects of the model. The new data provided evidence to support the model structure, and its predictive performance was then assessed against dog and rat datasets not used in model development. The PQ-specific model parameters are the same for all three species, with only the physiological parameters varying between species. This consistency across species provides a strong basis for extrapolation to other species, as demonstrated here for the dog. The model enables a wide range of PQ data to be linked together to provide a broad understanding of PQ pharmacokinetics in rodents and the dog, showing that the key aspects of PQ kinetics in these species are understood and adequately encapsulated within the model.

Using quantitative modeling tools to assess pharmacokinetic bias in epidemiological studies showing associations between biomarkers and health outcomes at low exposures.

Biomarkers of exposure can be measured at lower and lower levels due to advances in analytical chemistry. Using these sensitive methods, some epidemiology studies report associations between biomarkers and health outcomes at biomarker levels much below those associated with effects in animal studies. While some of these low exposure associations may arise from increased sensitivity of humans compared with animals or from species-specific responses, toxicology studies with drugs, commodity chemicals and consumer products have not generally indicated significantly greater sensitivity of humans compared with test animals for most health outcomes. In some cases, these associations may be indicative of pharmacokinetic (PK) bias, i.e., a situation where a confounding factor or the health outcome itself alters pharmacokinetic processes affecting biomarker levels. Quantitative assessment of PK bias combines PK modeling and statistical methods describing outcomes across large numbers of individuals in simulated populations. Here, we first provide background on the types of PK models that can be used for assessing biomarker levels in human population and then outline a process for considering PK bias in studies intended to assess associations between biomarkers and health outcomes at low levels of exposure. After providing this background, we work through published examples where these PK methods have been applied with several chemicals/chemical classes - polychlorinated biphenyls (PCBs), perfluoroalkyl substances (PFAS), polybrominated biphenyl ethers (PBDE) and phthalates - to assess the possibility of PK bias. Studies of the health effects of low levels of exposure will be improved by developing some confidence that PK bias did not play significant roles in the observed associations.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole.

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

A toxicogenomic approach for the risk assessment of the food contaminant acetamide.

Acetamide (CAS 60-35-5) is detected in common foods. Chronic rodent bioassays led to its classification as a group 2B possible human carcinogen due to the induction of liver tumors in rats. We used a toxicogenomics approach in Wistar rats gavaged daily for 7 or 28 days at doses of 300 to 1500 mg/kg/day (mkd) to determine a point of departure (POD) and investigate its mode of action (MoA). Ki67 labeling was increased at doses ≥750 mkd up to 3.3-fold representing the most sensitive apical endpoint. Differential gene expression analysis by RNA-Seq identified 1110 and 1814 differentially expressed genes in male and female rats, respectively, following 28 days of treatment. Down-regulated genes were associated with lipid metabolism while up-regulated genes included cell signaling, immune response, and cell cycle functions. Benchmark dose (BMD) modeling of the Ki67 labeling index determined the BMD10 lower confidence limit (BMDL10) as 190 mkd. Transcriptional BMD modeling revealed excellent concordance between transcriptional POD and apical endpoints. Collectively, these results indicate that acetamide is most likely acting through a mitogenic MoA, though specific key initiating molecular events could not be elucidated. A POD value of 190 mkd determined for cell proliferation is suggested for risk assessment purposes.

The role of fit-for-purpose assays within tiered testing approaches: A case study evaluating prioritized estrogen-active compounds in an in vitro human uterotrophic assay.

Chemical risk assessment relies on toxicity tests that require significant numbers of animals, time and costs. For the >30,000 chemicals in commerce, the current scale of animal testing is insufficient to address chemical safety concerns as regulatory and product stewardship considerations evolve to require more comprehensive understanding of potential biological effects, conditions of use, and associated exposures. We demonstrate the use of a multi-level new approach methodology (NAMs) strategy for hazard- and risk-based prioritization to reduce animal testing. A Level 1/2 chemical prioritization based on estrogen receptor (ER) activity and metabolic activation using ToxCast data was used to select 112 chemicals for testing in a Level 3 human uterine cell estrogen response assay (IKA assay). The Level 3 data were coupled with quantitative in vitro to in vivo extrapolation (Q-IVIVE) to support bioactivity determination (as a surrogate for hazard) in a tissue-specific context. Assay AC50s and Q-IVIVE were used to estimate human equivalent doses (HEDs), and HEDs were compared to rodent uterotrophic assay in vivo-derived points of departure (PODs). For substances active both in vitro and in vivo, IKA assay-derived HEDs were lower or equivalent to in vivo PODs for 19/23 compounds (83%). Activity exposure relationships were calculated, and the IKA assay was as or more protective of human health than the rodent uterotrophic assay for all IKA-positive compounds. This study demonstrates the utility of biologically relevant fit-for-purpose assays and supports the use of a multi-level strategy for chemical risk assessment.

Gene co-regulation and co-expression in the aryl hydrocarbon receptor-mediated transcriptional regulatory network in the mouse liver.

Four decades after its discovery, the aryl hydrocarbon receptor (AHR), a ligand-inducible transcription factor (TF) activated by the persistent environmental contaminant 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD), remains an enigmatic molecule with a controversial endogenous role. Here, we have assembled a global map of the AHR gene regulatory network in female C57BL/6 mice orally gavaged with 30 µg/kg of TCDD from a combination of previously published gene expression and genome-wide TF-binding data sets. Using Kohonen self-organizing maps and subspace clustering, we show that genes co-regulated by common upstream TFs in the AHR network exhibit a pattern of co-expression. Directly bound, indirectly bound, and non-genomic AHR target genes exhibit distinct expression patterns, with the directly bound targets associated with highest median expression. Interestingly, among the directly bound AHR target genes, the expression level increases with the number of AHR-binding sites in the proximal promoter regions. Finally, we show that co-regulated genes in the AHR network activate distinct groups of downstream biological processes. Although the specific findings described here are restricted to hepatic effects under short-term TCDD exposure, this work describes a generalizable approach to the reconstruction and analysis of transcriptional regulatory cascades underlying cellular stress response, revealing network hierarchy and the nature of information flow from the initial signaling events to phenotypic outcomes. Such reconstructed networks can form the basis of a new generation of quantitative adverse outcome pathways.

Updating physiologically based pharmacokinetic models for manganese by incorporating rapid association/dissociation processes in tissues.

Previously, we developed a series of physiologically based pharmacokinetic (PBPK) models for manganese (Mn) in which saturable tissue binding and dose-dependent increases in biliary excretion captured key aspects of Mn homeostasis biology. These models reproduced the non-linear behavior of Mn kinetics in different tissues, accounting for dose-dependent changes in Mn kinetics. The original model construct had relatively slow association and dissociation rate constants for Mn binding in tissues. In this updated model, both rates of entry into tissue and the interaction of Mn with binding sites are rapid, and the step limiting Mn accumulation is the saturation of tissue binding sites. This binding reflects general cellular requirements for Mn with high affinity but rapid exchange between bound and free forms, which we captured using a dissociation constant (KD) of ~ 0.5 µM across tissues while maintaining different maximum binding capacities in each tissue. Variability in the binding capacities accounted for different background levels of Mn in particular tissues. This alternative structure successfully described Mn kinetics in tissues in adult rats exposed to Mn either in their diet or by inhalation, indicating that both the original and the present models capture the dose-dependent and tissue-specific kinetic behavior of Mn in adult rats. Although the published models that emphasize the role of smaller tissue binding rate constants in non-linear behaviors capture all relevant dose-dependent kinetic behaviors of this metal, increasing biological relevance of the model structure and parameters should provide greater confidence in applying the Mn PBPK models to risk assessment.

Incorporation of in vitro metabolism data and physiologically based pharmacokinetic modeling in a risk assessment for chloroprene.

Objective: To develop a physiologically based pharmacokinetic (PBPK) model for chloroprene in the mouse, rat and human, relying only on in vitro data to estimate tissue metabolism rates and partitioning, and to apply the model to calculate an inhalation unit risk (IUR) for chloroprene.Materials and methods: Female B6C3F1 mice were the most sensitive species/gender for lung tumors in the 2-year bioassay conducted with chloroprene. The PBPK model included tissue metabolism rate constants for chloroprene estimated from results of in vitro gas uptake studies using liver and lung microsomes. To assess the validity of the PBPK model, a 6-hr, nose-only chloroprene inhalation study was conducted with female B6C3F1 mice in which both chloroprene blood concentrations and ventilation rates were measured. The PBPK model was then used to predict dose measures - amounts of chloroprene metabolized in lungs per unit time - in mice and humans.Results: The mouse PBPK model accurately predicted in vivo pharmacokinetic data from the 6-hr, nose-only chloroprene inhalation study. The PBPK model was used to conduct a cancer risk assessment based on metabolism of chloroprene to reactive epoxides in the lung, the target tissue in mice. The IUR was over100-fold lower than the IUR from the EPA Integrated Risk Information System (IRIS), which was based on inhaled chloroprene concentration. The different result from the PBPK model risk assessment arises from use of the more relevant tissue dose metric, amount metabolized, rather than inhaled concentrationDiscussion and conclusions: The revised chloroprene PBPK model is based on the best available science, including new test animal in vivo validation, updated literature review and a Markov-Chain Monte Carlo analysis to assess parameter uncertainty. Relying on both mouse and human metabolism data also provides an important advancement in the use of quantitative in vitro to in vivo extrapolation (QIVIVE). Inclusion of the best available science is especially important when deriving a toxicity value based on species extrapolation for the potential carcinogenicity of a reactive metabolite.

Considerations for refining the risk assessment process for formaldehyde: Results from an interdisciplinary workshop.

Anticipating the need to evaluate and integrate scientific evidence to inform new risk assessments or to update existing risk assessments, the Formaldehyde Panel of the American Chemistry Council (ACC), in collaboration with the University of North Carolina, convened a workshop: "Understanding Potential Human Health Cancer Risk - From Data Integration to Risk Evaluation" in October 2017. Twenty-four (24) invited-experts participated with expertise in epidemiology, toxicology, science integration and risk evaluation. Including members of the organizing committee, there were 29 participants. The meeting included eleven presentations encompassing an introduction and three sessions: (1) "integrating the formaldehyde science on nasal/nasopharyngeal carcinogenicity and potential for causality"; (2) "integrating the formaldehyde science on lymphohematopoietic cancer and potential for causality; and, (3) "formaldehyde research-data suitable for risk assessment". Here we describe key points from the presentations on epidemiology, toxicology and mechanistic studies that should inform decisions about the potential carcinogenicity of formaldehyde in humans and the discussions about approaches for structuring an integrated, comprehensive risk assessment for formaldehyde. We also note challenges expected when attempting to reconcile divergent results observed from research conducted within and across different scientific disciplines - especially toxicology and epidemiology - and in integrating diverse, multi-disciplinary mechanistic evidence.

Adipocyte-specific deficiency of Nfe2l1 disrupts plasticity of white adipose tissues and metabolic homeostasis in mice.

The maintenance of healthy adipose tissues is essential for efficient regulation of energy homeostasis. Nuclear factor-erythroid 2-related factor 1 (NFE2L1, also known as Nrf1), a CNC-bZIP protein, is a master regulator of the cellular adaptive response to stresses. To investigate the role of NFE2L1 in adipocytes, we bred a line of mice with adipocyte-specific Nfe2l1 knockout (Nfe2l1(f)-KO), and found that Nfe2l1(f)-KO mice exhibited a dramatically reduced subcutaneous adipose tissue (SAT) mass, insulin resistance, adipocyte hypertrophy, and severe adipose inflammation. Mechanistic studies revealed that Nfe2l1 deficiency may disturb the expression of lipolytic genes in adipocytes, leading to adipocyte hypertrophy followed by inflammation, pyroptosis, and insulin resistance. Our findings reveal a novel role for NFE2L1 in regulating adipose tissue plasticity and energy homeostasis.

Physiologically-based pharmacokinetic modeling suggests similar bioavailability of Mn from diet and drinking water.

Due to concerns for enhanced absorption of manganese (Mn) from drinking water compared to diet, bioavailability of Mn from drinking water remains a major data gap in understanding Mn kinetics. In this study, PBPK models for adult rats and humans were updated with a drinking water exposure route and were used to assess the homeostatic control of Mn uptake, excretion and tissue kinetics between the two different ingestion modes. Drinking water model parameters were estimated from tissue kinetic data from a drinking water study in rats. The published study included a 10 ppm-Mn diet with additional Mn added to drinking water to give a total ingested Mn dose equivalent to that from a 200 ppm diet. The 200 ppm diet and equivalent mixed drinking water/diet exposures provided Mn concentrations for brain (striatum, olfactory bulb, and cerebellum), liver and bone after 7 and 61 days of Mn exposure. Modeling of these data sets indicated that (1) the oral Mn bioavailability is similar for diet or drinking water and (2) homeostatic control of gut uptake of Mn occurs with either drinking water or dietary ingestion. This updated description for absorption and distribution of Mn from gut was added to a human Mn-PBPK model to simulate Mn exposure from multiple routes of exposure (i.e. dietary intake, drinking water, and inhalation). This increases the utility of the Mn PBPK model by allowing for the simulation of multiple Mn exposure scenarios, including variable daily food and drinking water exposures in a human population.

Based on an analysis of mode of action, styrene-induced mouse lung tumors are not a human cancer concern.

Based on 13 chronic studies, styrene exposure causes lung tumors in mice, but no tumor increases in other organs in mice or rats. Extensive research into the mode of action demonstrates the key events and human relevance. Key events are: metabolism of styrene by CYP2F2 in mouse lung club cells to ring-oxidized metabolites; changes in gene expression for metabolism of lipids and lipoproteins, cell cycle and mitotic M-M/G1 phases; cytotoxicity and mitogenesis in club cells; and progression to preneoplastic/neoplastic lesions in lung. Although styrene-7,8-oxide (SO) is a common genotoxic styrene metabolite in in vitro studies, the data clearly demonstrate that SO is not the proximate toxicant and that styrene does not induce a genotoxic mode of action. Based on complete attenuation of styrene short-term and chronic toxicity in CYP2F2 knockout mice and similar attenuation in CYP2F1 (humanized) transgenic mice, limited metabolism of styrene in human lung by CYP2F1, 2 + orders of magnitude lower SO levels in human lung compared to mouse lung, and lack of styrene-related increase in lung cancer in humans, styrene does not present a risk of cancer to humans.

Strain-related differences in mouse lung gene expression over a two-year period of inhalation exposure to styrene: Relevance to human risk assessment.

Both CD-1 and C57BL/6 wildtype (C57BL/6-WT) mice show equivalent short-term lung toxicity from exposures to styrene, while long-term tumor responses are greater in CD-1 mice. We analyzed lung gene expression from styrene exposures lasting from 1-day to 2-years in male mice from these two strains, including a Cyp2f2(-/-) knockout (C57BL/6-KO) and a Cyp2F1/2A13/2B6 transgenic mouse (C57BL/6-TG). With short term exposures (1-day to 1-week), CD-1 and C57BL/6-WT mice had thousands of differentially expressed genes (DEGs), consistent with changes in pathways for cell proliferation, cellular lipid metabolism, DNA-replication and inflammation. C57BL/6-WT mice responded within a single day; CD-1 mice required several days of exposure. The numbers of exposure related DEGs were greatly reduced at longer times (4-weeks to 2-years) with enrichment only for biological oxidations in C57BL/6-WT and metabolism of lipids and lipoproteins in CD-1. Gene expression results indicate a non-genotoxic, mouse specific mode of action for short-term styrene responses related to activation of nuclear receptor signaling and cell proliferation. Greater tumor susceptibility in CD-1 mice correlated with the presence of the Pas1 loci, differential Cytochrome P450 gene expression, down-regulation of Nr4a, and greater inflammatory pathway activation. Very few exposure-related responses occurred at any time in C57BL/6-KO or -TG mice indicating that neither the short term nor long term responses of styrene in mice are relevant endpoints for assessing human risks.

Assessing molecular initiating events (MIEs), key events (KEs) and modulating factors (MFs) for styrene responses in mouse lungs using whole genome gene expression profiling following 1-day and multi-week exposures.

Styrene increased lung tumors in mice at chronic inhalation exposures of 20ppm and greater. MIEs, KEs and MFs were examined using gene expression in three strains of male mice (the parental C57BL/6 strain, a CYP2F2(-/-) knock out and a CYP2F2(-/-) transgenic containing human CYP2F1, 2A13 and 2B6). Exposures were for 1-day and 1, 4 and 26weeks. After 1-day exposures at 1, 5, 10, 20, 40 and 120ppm significant increases in differentially expressed genes (DEGs) occurred only in parental strain lungs where there was already an increase in DEGs at 5ppm and then many thousands of DEGs by 120ppm. Enrichment for 1-day and 1-week exposures included cell cycle, mitotic M-M/G1 phases, DNA-synthesis and metabolism of lipids and lipoproteins pathways. The numbers of DEGs decreased steadily over time with no DEGs meeting both statistical significance and fold-change criteria at 26weeks. At 4 and 26weeks, some key transcription factors (TFs) - Nr1d1, Nr1d2, Dbp, Tef, Hlf, Per3, Per2 and Bhlhe40 - were upregulated (|FC|>1.5), while others - Npas, Arntl, Nfil3, Nr4a1, Nr4a2, and Nr4a3 - were down-regulated. At all times, consistent changes in gene expression only occurred in the parental strain. Our results support a MIE for styrene of direct mitogenicity from mouse-specific CYP2F2-mediated metabolites activating Nr4a signaling. Longer-term MFs include down-regulation of Nr4a genes and shifts in both circadian clock TFs and other TFs, linking circadian clock to cellular metabolism. We found no gene expression changes indicative of cytotoxicity or activation of p53-mediated DNA-damage pathways.

A tissue dose-based comparative exposure assessment of manganese using physiologically based pharmacokinetic modeling-The importance of homeostatic control for an essential metal.

A physiologically-based pharmacokinetic (PBPK) model (Schroeter et al., 2011) was applied to simulate target tissue manganese (Mn) concentrations following occupational and environmental exposures. These estimates of target tissue Mn concentrations were compared to determine margins of safety (MOS) and to evaluate the biological relevance of applying safety factors to derive acceptable Mn air concentrations. Mn blood concentrations measured in occupational studies permitted verification of the human PBPK models, increasing confidence in the resulting estimates. Mn exposure was determined based on measured ambient air Mn concentrations and dietary data in Canada and the United States (US). Incorporating dietary and inhalation exposures into the models indicated that increases in target tissue concentrations above endogenous levels only begin to occur when humans are exposed to levels of Mn in ambient air (i.e. >10µg/m3) that are far higher than those currently measured in Canada or the US. A MOS greater than three orders of magnitude was observed, indicating that current Mn air concentrations are far below concentrations that would be required to produce the target tissue Mn concentrations associated with subclinical neurological effects. This application of PBPK modeling for an essential element clearly demonstrates that the conventional application of default factors to "convert" an occupational exposure to an equivalent continuous environmental exposure, followed by the application of safety factors, is not appropriate in the case of Mn. PBPK modeling demonstrates that the relationship between ambient Mn exposures and dose-to-target tissue is not linear due to normal tissue background levels and homeostatic controls.