Hepatic transcriptional dose-response analysis of male and female Fischer rats exposed to hexabromocyclododecane

Derivation of Transcriptomics-Based Points of Departure for 20 Per- or Polyfluoroalkyl Substances Using a Larval Fathead Minnow (Pimephales promelas) Reduced Transcriptome Assay

Traditional toxicity testing has been unable to keep pace with the introduction of new chemicals into commerce. Consequently, there are limited or no toxicity data for many chemicals to which fish and wildlife may be exposed. Per- and polyfluoroalkyl substances (PFAS) are emblematic of this issue in that ecological hazards of most PFAS remain uncharacterized. The present study employed a high-throughput assay to identify the concentration at which 20 PFAS, with diverse properties, elicited a concerted gene expression response (termed a transcriptomics-based point of departure [tPOD]) in larval fathead minnows (Pimephales promelas; 5-6 days postfertilization) exposed for 24 h. Based on a reduced transcriptome approach that measured whole-body expression of 1832 genes, the median tPOD for the 20 PFAS tested was 10 µM. Longer-chain carboxylic acids (12-13 C-F); an eight-C-F dialcohol, N-alkyl sulfonamide; and telomer sulfonic acid were among the most potent PFAS, eliciting gene expression responses at concentrations <1 µM. With a few exceptions, larval fathead minnow tPODs were concordant with those based on whole-transcriptome response in human cell lines. However, larval fathead minnow tPODs were often greater than those for Daphnia magna exposed to the same PFAS. The tPODs overlapped concentrations at which other sublethal effects have been reported in fish (available for 10 PFAS). Nonetheless, fathead minnow tPODs were orders of magnitude higher than aqueous PFAS concentrations detected in tributaries of the North American Great Lakes, suggesting a substantial margin of safety. Overall, results broadly support the use of a fathead minnow larval transcriptomics assay to derive screening-level potency estimates for use in ecological risk-based prioritization. Environ Toxicol Chem 2024;00:1-16. © 2024 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

The brominated flame retardant tetrabromobisphenol A-bis(2,3-dibromo-2-methylpropyl ether) as well as hexabromocyclododecane lead to lipid disorders in mice

The brominated flame retardant tetrabromobisphenol A-bis(2,3-dibromo-2-methylpropyl ether) (TBBPA-DBMPE) is a recommended substitute for hexabromocyclododecane (HBCD), a banned persistent organic pollutant, yet its potential toxicities remains largely unexplored. Here, we investigated the effects of a long-term exposure to TBBPA-DBMPE at nominal doses of 50 and 1000 μg/kg/d on lipid homeostasis in CD-1 mice, in comparison with 50 μg/kg/d HBCD as a positive control. Male pups received chemical treatments through maternal administration via drinking water from postnatal day 0-21, followed by direct administration through drinking water after weaning. On the 23rd week after treatment, the oral lipid tolerance test revealed that low-dose TBBPA-DBMPE as well as HBCD affected lipid tolerance, although the fasting serum triglyceride (TG) levels were not altered. When chemical treatment was extended to the 32nd week, TBBPA-DBMPE-treated animals displayed adipocyte hypertrophy in both white adipose tissue (eWAT) and brown adipose tissue (BAT) and hepatic steatosis, which was largely consistent with the effects of HBCD. These findings indicate that like HBCD, TBBPA-DBMPE led to increased lipid load in mice. Interestingly, we also observed intestinal histological changes, coupled with increased expression of lipid absorption-related genes in both HBCD and TBBPA-DBMPE treatments, suggesting increased lipid absorption. This was supported by in vitro findings that both HBCD and TBBPA-DBMPE promoted lipid accumulation in IEC-6 cells under the stress of oleic acid for 6 h, implying that altered lipid absorption by the intestine may partly contributed to increased lipid load in mice. Overall, the effects of 50 μg/kg/d TBBPA-DBMPE in terms of some parameters were comparable with 50 μg/kg/d HBCD, suggesting that TBBPA-DBMPE may not be an ideal substitute of HBCD.

Persistent immune injury induced by short-term decabromodiphenyl ether (BDE-209) exposure to female middle-aged Balb/c mice

Decabromodiphenyl ether (BDE-209), widely used in various industries for its excellent flame-retardant performance, could be enriched in humans and is closely associated with immune impairment. In addition, immune system is gradually declined and becoming more sensitive to environmental pollutants in the ageing process. Therefore, the immunotoxicity of BDE-209 (4, 40, and 400 mg/kg/day) to middle-aged mice and its recovery and susceptibility was first to be comprehensively investigated in this study. The results showed that BDE-209 exposure could lead to oxidative injury to immune organs (spleen, thymus, and liver), impair humoral (immunoglobulins), cellular (lymphopoiesis), and non-specific immunity, and disturb the expressions of the genes related to Th1/Th2 balance (T helper cells) in the middle-aged mice. In addition, Integrated Biomarker Response (IBR) indicated that BDE-209-induced immune impairment was challenging to self-regulated, and even exacerbated after 21 days of recovery and oxidative injury in immune organs could be the main reason. Furthermore, factorial analysis showed that middle-aged mice exposed to BDE-209 suffered from greater immune impairment than adult mice, and the immune impairment in aged mice is more difficult to be self-repaired than that in adult mice. It can be seen that the aged tend to suffer from BDE-209-induced persistent immune impairment and health threats.

Adverse outcome pathway for pregnane X receptor-induced hypercholesterolemia

Pharmaceuticals and environmental contaminants contribute to hypercholesterolemia. Several chemicals known to cause hypercholesterolemia, activate pregnane X receptor (PXR). PXR is a nuclear receptor, classically identified as a sensor of chemical environment and regulator of detoxification processes. Later, PXR activation has been shown to disrupt metabolic functions such as lipid metabolism and recent findings have shown PXR activation to promote hypercholesterolemia through multiple mechanisms. Hypercholesterolemia is a major causative risk factor for atherosclerosis and greatly promotes global health burden. Metabolic disruption by PXR activating chemicals leading to hypercholesterolemia represents a novel toxicity pathway of concern and requires further attention. Therefore, we constructed an adverse outcome pathway (AOP) by collecting the available knowledge considering the molecular mechanisms for PXR-mediated hypercholesterolemia. AOPs are tools of modern toxicology for systematizing mechanistic knowledge to assist health risk assessment of chemicals. AOPs are formalized and structured linear concepts describing a link between molecular initiating event (MIE) and adverse outcome (AO). MIE and AO are connected via key events (KE) through key event relationships (KER). We present a plausible route of how PXR activation (MIE) leads to hypercholesterolemia (AO) through direct regulation of cholesterol synthesis and via activation of sterol regulatory element binding protein 2-pathway.

From vision toward best practices: Evaluating in vitro transcriptomic points of departure for application in risk assessment using a uniform workflow

The growing number of chemicals in the current consumer and industrial markets presents a major challenge for regulatory programs faced with the need to assess the potential risks they pose to human and ecological health. The increasing demand for hazard and risk assessment of chemicals currently exceeds the capacity to produce the toxicity data necessary for regulatory decision making, and the applied data is commonly generated using traditional approaches with animal models that have limited context in terms of human relevance. This scenario provides the opportunity to implement novel, more efficient strategies for risk assessment purposes. This study aims to increase confidence in the implementation of new approach methods in a risk assessment context by using a parallel analysis to identify data gaps in current experimental designs, reveal the limitations of common approaches deriving transcriptomic points of departure, and demonstrate the strengths in using high-throughput transcriptomics (HTTr) to derive practical endpoints. A uniform workflow was applied across six curated gene expression datasets from concentration-response studies containing 117 diverse chemicals, three cell types, and a range of exposure durations, to determine tPODs based on gene expression profiles. After benchmark concentration modeling, a range of approaches was used to determine consistent and reliable tPODs. High-throughput toxicokinetics were employed to translate in vitro tPODs (µM) to human-relevant administered equivalent doses (AEDs, mg/kg-bw/day). The tPODs from most chemicals had AEDs that were lower (i.e., more conservative) than apical PODs in the US EPA CompTox chemical dashboard, suggesting in vitro tPODs would be protective of potential effects on human health. An assessment of multiple data points for single chemicals revealed that longer exposure duration and varied cell culture systems (e.g., 3D vs. 2D) lead to a decreased tPOD value that indicated increased chemical potency. Seven chemicals were flagged as outliers when comparing the ratio of tPOD to traditional POD, thus indicating they require further assessment to better understand their hazard potential. Our findings build confidence in the use of tPODs but also reveal data gaps that must be addressed prior to their adoption to support risk assessment applications.

Occurrence and Health Effects of Hexabromocyclododecane: An Updated Review

Hexabromocyclododecane (HBCD) is a non-aromatic compound belonging to the bromine flame retardant family and is a known persistent organic pollutant (POP). This compound accumulates easily in the environment and has a high half-life in water. With a variety of uses, the HBCD is found in house dust, electronics, insulation, and construction. There are several isomers and the most studied are α-, β-, and γ-HBCD. Initially used as a substitute for other flame retardants, the polybrominated diphenyl ethers (PBDEs), the discovery of its role as a POP made HBCD use and manufacturing restricted in Europe and other countries. The adverse effects on the environment and human health have been piling, either as a result from its accumulation or considering its power as an endocrine disruptor (ED). Furthermore, it has also been proven that it has detrimental effects on the neuronal system, endocrine system, cardiovascular system, liver, and the reproductive system. HBCD has also been linked to cytokine production, DNA damage, increased cell apoptosis, increased oxidative stress, and reactive oxygen species (ROS) production. Therefore, this review aims to compile the most recent studies regarding the negative effects of this compound on the environment and human health, describing the possible mechanisms by which this compound acts and its possible toxic effects.

Considerations for application of benchmark dose modeling in radiation research: workshop highlights

BACKGROUND

Exposure to different forms of ionizing radiation occurs in diverse occupational, medical, and environmental settings. Improving the accuracy of the estimated health risks associated with exposure is therefore, essential for protecting the public, particularly as it relates to chronic low dose exposures. A key aspect to understanding health risks is precise and accurate modeling of the dose-response relationship. Toward this vision, benchmark dose (BMD) modeling may be a suitable approach for consideration in the radiation field. BMD modeling is already extensively used for chemical hazard assessments and is considered statistically preferable to identifying low and no observed adverse effects levels. BMD modeling involves fitting mathematical models to dose-response data for a relevant biological endpoint and identifying a point of departure (the BMD, or its lower bound). Recent examples in chemical toxicology show that when applied to molecular endpoints (e.g. genotoxic and transcriptional endpoints), BMDs correlate to points of departure for more apical endpoints such as phenotypic changes (e.g. adverse effects) of interest to regulatory decisions. This use of BMD modeling may be valuable to explore in the radiation field, specifically in combination with adverse outcome pathways, and may facilitate better interpretation of relevant in vivo and in vitro dose-response data. To advance this application, a workshop was organized on June 3rd, 2022, in Ottawa, Ontario that brought together BMD experts in chemical toxicology and the radiation scientific community of researchers, regulators, and policy-makers. The workshop's objective was to introduce radiation scientists to BMD modeling and its practical application using case examples from the chemical toxicity field and demonstrate the BMDExpress software using a radiation dataset. Discussions focused on the BMD approach, the importance of experimental design, regulatory applications, its use in supporting the development of adverse outcome pathways, and specific radiation-relevant examples.

CONCLUSIONS

Although further deliberations are needed to advance the use of BMD modeling in the radiation field, these initial discussions and partnerships highlight some key steps to guide future undertakings related to new experimental work.

The effects of hexabromocyclododecane on the transcriptome and hepatic enzyme activity in three human HepaRG-based models

The disruption of thyroid hormone homeostasis by hexabromocyclododecane (HBCD) in rodents is hypothesized to be due to HBCD increasing the hepatic clearance of thyroxine (T4). The extent to which these effects are relevant to humans is unclear. To evaluate HBCD effects on humans, the activation of key hepatic nuclear receptors and the consequent disruption of thyroid hormone homeostasis were studied in different human hepatic cell models. The hepatoma cell line, HepaRG, cultured as two-dimensional (2D), sandwich (SW) and spheroid (3D) cultures, and primary human hepatocytes (PHH) cultured as sandwich were exposed to 1 and 10 µM HBCD and characterized for their transcriptome changes. Pathway enrichment analysis showed that 3D models, followed by SW, had a stronger transcriptome response to HBCD, which is explained by the higher expression of hepatic nuclear receptors but also greater accumulation of HBCD measured inside cells in these models. The Pregnane X receptor pathway is one of the pathways most upregulated across the three hepatic models, followed by the constitutive androstane receptor and general hepatic nuclear receptors pathways. Lipid metabolism pathways had a downregulation tendency in all exposures and in both PHH and the three cultivation modes of HepaRG. The activity of enzymes related to PXR/CAR induction and T4 metabolism were evaluated in the three different types of HepaRG cultures exposed to HBCD for 48 h. Reference inducers, rifampicin and PCB-153 did affect 2D and SW HepaRG cultures' enzymatic activity but not 3D. HBCD did not induce the activity of any of the studied enzymes in any of the cell models and culture methods. This study illustrates that for nuclear receptor-mediated T4 disruption, transcriptome changes might not be indicative of an actual adverse effect. Clarification of the reasons for the lack of translation is essential to evaluate new chemicals' potential to be thyroid hormone disruptors by altering thyroid hormone metabolism.

Microbiota and Transcriptomic Effects of an Essential Oil Blend and Its Delivery Route Compared to an Antibiotic Growth Promoter in Broiler Chickens

This study evaluated the effect of the delivery of a commercial essential oil blend containing the phytonutrients star anise, cinnamon, rosemary, and thyme oil (via different routes) on broiler chickens’ ileal and ceca microbiota and liver transcriptome compared to an antibiotic growth promoter. Eggs were incubated and allocated into three groups: non-injected, in ovo saline, and in ovo essential oil. On day 18 of incubation, 0.2 mL of essential oil in saline (dilution ratio of 2:1) or saline alone was injected into the amnion. At hatch, chicks were assigned to post-hatch treatment combinations: (A) a negative control (corn-wheat-soybean diet), (B) in-feed antibiotics, (C) in-water essential oil (250 mL/1000 L of drinking water), (D) in ovo saline, (E) in ovo essential oil, and (F) in ovo essential oil plus in-water essential oil in eight replicate cages (six birds/cage) and raised for 28 days. On days 21 and 28, one and two birds per cage were slaughtered, respectively, to collect gut content and liver tissues for further analysis. Alpha and beta diversity differed significantly between ileal and ceca samples but not between treatment groups. In-feed antibiotic treatment significantly increased the proportion of specific bacteria in the family Lachnospiraceae while reducing the proportion of bacteria in the genus Christensenellaceae in the ceca, compared to other treatments. Sex-controlled differential expression of genes related to cell signaling and tight junctions were recorded. This study provides data that could guide the use of these feed additives and a foundation for further research.

Latent Variables Capture Pathway-Level Points of Departure in High-Throughput Toxicogenomic Data

Estimation of points of departure (PoDs) from high-throughput transcriptomic data (HTTr) represents a key step in the development of next-generation risk assessment (NGRA). Current approaches mainly rely on single key gene targets, which are constrained by the information currently available in the knowledge base and make interpretation challenging as scientists need to interpret PoDs for thousands of genes or hundreds of pathways. In this work, we aimed to address these issues by developing a computational workflow to investigate the pathway concentration–response relationships in a way that is not fully constrained by known biology and also facilitates interpretation. We employed the Pathway-Level Information ExtractoR (PLIER) to identify latent variables (LVs) describing biological activity and then investigated in vitro LVs’ concentration–response relationships using the ToxCast pipeline. We applied this methodology to a published transcriptomic concentration–response data set for 44 chemicals in MCF-7 cells and showed that our workflow can capture known biological activity and discriminate between estrogenic and antiestrogenic compounds as well as activity not aligning with the existing knowledge base, which may be relevant in a risk assessment scenario. Moreover, we were able to identify the known estrogen activity in compounds that are not well-established ER agonists/antagonists supporting the use of the workflow in read-across. Next, we transferred its application to chemical compounds tested in HepG2, HepaRG, and MCF-7 cells and showed that PoD estimates are in strong agreement with those estimated using a recently developed Bayesian approach (cor = 0.89) and in weak agreement with those estimated using a well-established approach such as BMDExpress2 (cor = 0.57). These results demonstrate the effectiveness of using PLIER in a concentration–response scenario to investigate pathway activity in a way that is not fully constrained by the knowledge base and to ease the biological interpretation and support the development of an NGRA framework with the ability to improve current risk assessment strategies for chemicals using new approach methodologies.

Integration of Toxicogenomics and Physiologically Based Pharmacokinetic Modeling in Human Health Risk Assessment of Perfluorooctane Sulfonate

Toxicogenomics and physiologically based pharmacokinetic (PBPK) models are useful approaches in chemical risk assessment, but the methodology to incorporate toxicogenomic data into a PBPK model to inform risk assessment remains to be developed. This study aimed to develop a probabilistic human health risk assessment approach by integrating toxicogenomic dose-response data and PBPK modeling using perfluorooctane sulfonate (PFOS) as a case study. Based on the available human in vitro and mouse in vivo toxicogenomic data, we identified the differentially expressed genes (DEGs) at each exposure paradigm/duration. Kyoto Encyclopedia of Genes and Genomes and disease ontology enrichment analyses were conducted on the DEGs to identify significantly enriched pathways and diseases. The dose-response data of DEGs were analyzed using the Bayesian benchmark dose (BMD) method. Using a previously published PBPK model, the gene BMDs were converted to human equivalent doses (HEDs), which were summarized to pathway and disease HEDs and then extrapolated to reference doses (RfDs) by considering an uncertainty factor of 30 for mouse in vivo data and 10 for human in vitro data. The results suggested that the median RfDs at different exposure paradigms were similar to the 2016 U.S. Environmental Protection Agency's recommended RfD, while the RfDs for the most sensitive pathways and diseases were closer to the recent European Food Safety Authority's guidance values. In conclusion, genomic dose-response data and PBPK modeling can be integrated to become a useful alternative approach in risk assessment of environmental chemicals. This approach considers multiple endpoints, provides toxicity mechanistic insights, and does not rely on apical toxicity endpoints.

R-ODAF: Omics data analysis framework for regulatory application

Despite the widespread use of transcriptomics technologies in toxicology research, acceptance of the data by regulatory agencies to support the hazard assessment is still limited. Fundamental issues contributing to this are the lack of reproducibility in transcriptomics data analysis arising from variance in the methods used to generate data and differences in the data processing. While research applications are flexible in the way the data are generated and interpreted, this is not the case for regulatory applications where an unambiguous answer, possibly later subject to legal scrutiny, is required. A reference analysis framework would give greater credibility to the data and allow the practitioners to justify their use of an alternative bioinformatic process by referring to a standard. In this publication, we propose a method called omics data analysis framework for regulatory application (R-ODAF), which has been built as a user-friendly pipeline to analyze raw transcriptomics data from microarray and next-generation sequencing. In the R-ODAF, we also propose additional statistical steps to remove the number of false positives obtained from standard data analysis pipelines for RNA-sequencing. We illustrate the added value of R-ODAF, compared to a standard workflow, using a typical toxicogenomics dataset of hepatocytes exposed to paracetamol.

Benchmarking of a Bayesian single cell RNAseq differential gene expression test for dose-response study designs

The application of single-cell RNA sequencing (scRNAseq) for the evaluation of chemicals, drugs, and food contaminants presents the opportunity to consider cellular heterogeneity in pharmacological and toxicological responses. Current differential gene expression analysis (DGEA) methods focus primarily on two group comparisons, not multi-group dose-response study designs used in safety assessments. To benchmark DGEA methods for dose-response scRNAseq experiments, we proposed a multiplicity corrected Bayesian testing approach and compare it against 8 other methods including two frequentist fit-for-purpose tests using simulated and experimental data. Our Bayesian test method outperformed all other tests for a broad range of accuracy metrics including control of false positive error rates. Most notable, the fit-for-purpose and standard multiple group DGEA methods were superior to the two group scRNAseq methods for dose-response study designs. Collectively, our benchmarking of DGEA methods demonstrates the importance in considering study design when determining the most appropriate test methods.

Determination of Hexabromocyclododecane in Expanded Polystyrene and Extruded Polystyrene Foam by Gas Chromatography-Mass Spectrometry

A gas chromatography-mass spectrometry (GC/MS) method for the determination of hexabromocyclododecane (HBCD) in expanded polystyrene and extruded polystyrene foam (EPS/XPS) was developed. The EPS/XPS samples were ultrasonically extracted with acetone and the extracts were purified by filtration through a microporous membrane (0.22 μm) and solid-phase extraction. The samples were analyzed using a GC/MS using the selected ion monitoring mode. The ions 157, 319 and 401 were selected as the qualitative ions, while ion 239 was chosen as the quantitative ion. An HBCD standard working solution with a concentration range of 1.0–50.0 mg/L showed good linearity. The detection limit of HBCD was 0.5 mg/kg, meeting the LPC limit (<100 or 1000 mg/kg). Six laboratories were selected to verify the accuracy of the method, and 10 samples were tested. The interlaboratory relative standard deviation range was 3.68–9.80%. This method could play an important role in controlling HBCD contamination in EPS/XPS.

Transcriptomic Analysis of the Differential Nephrotoxicity of Diverse Brominated Flame Retardants in Rat and Human Renal Cells

Brominated flame retardants (BFRs) are environmentally persistent, are detected in humans, and some have been banned due to their potential toxicity. BFRs are developmental neurotoxicants and endocrine disruptors; however, few studies have explored their potential nephrotoxicity. We addressed this gap in the literature by determining the toxicity of three different BFRs (tetrabromobisphenol A (TBBPA), hexabromocyclododecane (HBCD), and tetrabromodiphenyl ether (BDE-47)) in rat (NRK 52E) and human (HK-2 and RPTEC) tubular epithelial cells. All compounds induced time- and concentration-dependent toxicity based on decreases in MTT staining and changes in cell and nuclear morphology. The toxicity of BFRs was chemical- and cell-dependent, and human cells were more susceptible to all three BFRs based on IC₅₀s after 48 h exposure. BFRs also had chemical- and cell-dependent effects on apoptosis as measured by increases in annexin V and PI staining. The molecular mechanisms mediating this toxicity were investigated using RNA sequencing. Principal components analysis supported the hypothesis that BFRs induce different transcriptional changes in rat and human cells. Furthermore, BFRs only shared nine differentially expressed genes in rat cells and five in human cells. Gene set enrichment analysis demonstrated chemical- and cell-dependent effects; however, some commonalities were also observed. Namely, gene sets associated with extracellular matrix turnover, the coagulation cascade, and the SNS-related adrenal cortex response were enriched across all cell lines and BFR treatments. Taken together, these data support the hypothesis that BFRs induce differential toxicity in rat and human renal cell lines that is mediated by differential changes in gene expression.

Paving the way for application of next generation risk assessment to safety decision-making for cosmetic ingredients

Next generation risk assessment (NGRA) is an exposure-led, hypothesis-driven approach that has the potential to support animal-free safety decision-making. However, significant effort is needed to develop and test the in vitro and in silico (computational) approaches that underpin NGRA to enable confident application in a regulatory context. A workshop was held in Montreal in 2019 to discuss where effort needs to be focussed and to agree on the steps needed to ensure safety decisions made on cosmetic ingredients are robust and protective. Workshop participants explored whether NGRA for cosmetic ingredients can be protective of human health, and reviewed examples of NGRA for cosmetic ingredients. From the limited examples available, it is clear that NGRA is still in its infancy, and further case studies are needed to determine whether safety decisions are sufficiently protective and not overly conservative. Seven areas were identified to help progress application of NGRA, including further investments in case studies that elaborate on scenarios frequently encountered by industry and regulators, including those where a ‘high risk’ conclusion would be expected. These will provide confidence that the tools and approaches can reliably discern differing levels of risk. Furthermore, frameworks to guide performance and reporting should be developed.

Flame Retardants-Mediated Interferon Signaling in the Pathogenesis of Nonalcoholic Fatty Liver Disease

Nonalcoholic fatty liver disease (NAFLD) is a growing concern worldwide, affecting 25% of the global population. NAFLD is a multifactorial disease with a broad spectrum of pathology includes steatosis, which gradually progresses to a more severe condition such as nonalcoholic steatohepatitis (NASH), fibrosis, cirrhosis, and eventually leads to hepatic cancer. Several risk factors, including exposure to environmental toxicants, are involved in the development and progression of NAFLD. Environmental factors may promote the development and progression of NAFLD by various biological alterations, including mitochondrial dysfunction, reactive oxygen species production, nuclear receptors dysregulation, and interference in inflammatory and immune-mediated signaling. Moreover, environmental contaminants can influence immune responses by impairing the immune system’s components and, ultimately, disease susceptibility. Flame retardants (FRs) are anthropogenic chemicals or mixtures that are being used to inhibit or delay the spread of fire. FRs have been employed in several household and outdoor products; therefore, human exposure is unavoidable. In this review, we summarized the potential mechanisms of FRs-associated immune and inflammatory signaling and their possible contribution to the development and progression of NAFLD, with an emphasis on FRs-mediated interferon signaling. Knowledge gaps are identified, and emerging pharmacotherapeutic molecules targeting the immune and inflammatory signaling for NAFLD are also discussed.

Update of the risk assessment of hexabromocyclododecanes (HBCDDs) in food

The European Commission asked EFSA to update its 2011 risk assessment on hexabromocyclododecanes (HBCDDs) in food. HBCDDs, predominantly mixtures of the stereoisomers α-, β- and γ-HBCDD, were widely used additive flame retardants. Concern has been raised because of the occurrence of HBCDDs in the environment, food and in humans. Main targets for toxicity are neurodevelopment, the liver, thyroid hormone homeostasis and the reproductive and immune systems. The CONTAM Panel concluded that the neurodevelopmental effects on behaviour in mice can be considered the critical effects. Based on effects on spontaneous behaviour in mice, the Panel identified a lowest observed adverse effect level (LOAEL) of 0.9 mg/kg body weight (bw) as the Reference Point, corresponding to a body burden of 0.75 mg/kg bw. The chronic intake that would lead to the same body burden in humans was calculated to be 2.35 μg/kg bw per day. The derivation of a health-based guidance value (HBGV) was not considered appropriate. Instead, the margin of exposure (MOE) approach was applied to assess possible health concerns. Over 6,000 analytical results for HBCDDs in food were used to estimate the exposure across dietary surveys and age groups of the European population. The most important contributors to the chronic dietary LB exposure to HBCDDs were fish meat, eggs, livestock meat and poultry. The CONTAM Panel concluded that the resulting MOE values support the conclusion that current dietary exposure to HBCDDs across European countries does not raise a health concern. An exception is breastfed infants with high milk consumption, for which the lowest MOE values may raise a health concern.

Di-(2-ethylhexyl) phthalate-induced hepatotoxicity exacerbated type 2 diabetes mellitus (T2DM) in female pubertal T2DM mice

Di-(2-ethylhexyl) phthalate (DEHP), one of the most common plasticizers, is closely associated with a high prevalence of pubertal type 2 diabetes mellitus (T2DM). Numerous studies have indicated that DEHP-induced metabolic toxicity exhibits sex differences. In this study, the sex differences in the effect of DEHP on pubertal T2DM (P-T2DM) mice, the susceptibility of female P-T2DM mice to DEHP-induced metabolic toxicity, and the underlying mechanisms were investigated. DEHP exposure exacerbated metabolic disorders in female P-T2DM mice. Factorial analysis showed that female P-T2DM mice were more sensitive to DEHP exposure than female normal mice and male P-T2DM mice. It was determined by integrated biomarker response results that female P-T2DM mice had higher risks of developing T2DM, metabolic disorders, cardiovascular events and hepatotoxicity than male P-T2DM mice. Moreover, hepatic transcriptome analysis emphasized the effects of DEHP on the expression of oxidative injury- and metabolic function-related genes. Western blotting indicated that DEHP activated Jun-N-terminal kinase (JNK) and impaired insulin sensitivity in the liver, which were the main causes of DEHP-exacerbated metabolic abnormalities in P-T2DM mice. Our study revealed that compared with normal mice and male P-T2DM mice, female P-T2DM mice tend to suffer from increased DEHP-induced metabolic toxicity, which was primarily attributed to hepatotoxicity.

Probing the Mechanism of Hepatotoxicity of Hexabromocyclododecanes through Toxicological Network Analysis

The prediction and mechanism analysis of hepatotoxicity of contaminants, because of their various phenotypes and complex mechanisms, is still a key problem in environmental research. We applied a toxicological network analysis method to predict the hepatotoxicity of three hexabromocyclododecane (HBCD) diastereoisomers (α-HBCD, β-HBCD, and γ-HBCD) and explore their potential mechanisms. First, we collected the hepatotoxicity related genes and found that those genes were significantly localized in the human interactome. Therefore, these genes form a disease module of hepatotoxicity. We also collected targets of α-, β-, and γ-HBCD and found that their targets overlap with the hepatotoxicity disease module. Then, we trained a model to predict hepatotoxicity of three HBCD diastereoisomers based on the relationship between the hepatotoxicity disease module and targets of compounds. We found that 593 genes were significantly located in the hepatotoxicity disease module (Z = 11.9, p < 0.001) involved in oxidative stress, cellular immunity, and proliferation, and the accuracy of hepatotoxicity prediction of HBCD was 0.7095 ± 0.0193 and the recall score was 0.8355 ± 0.0352. HBCD mainly affects the core disease module genes to mediate the adenosine monophosphate-activated kinase, p38MAPK, PI3K/Akt, and TNFα pathways to regulate the immune reaction and inflammation. HBCD also induces the secretion of IL6 and STAT3 to lead hepatotoxicity by regulating NR3C1. This approach is transferable to other toxicity research studies of environmental pollutants.

A Rat Liver Transcriptomic Point of Departure Predicts a Prospective Liver or Non-liver Apical Point of Departure

Identifying a toxicity point of departure (POD) is a required step in human health risk characterization of crop protection molecules, and this POD has historically been derived from apical endpoints across a battery of animal-based toxicology studies. Using rat transcriptome and apical data for 79 molecules obtained from Open TG-GATES (Toxicogenomics Project-Genomics Assisted Toxicity Evaluation System) (632 datasets), the hypothesis was tested that a short-term exposure, transcriptome-based liver biological effect POD (BEPOD) could estimate a longer-term exposure "systemic" apical endpoint POD. Apical endpoints considered were body weight, clinical observation, kidney weight and histopathology and liver weight and histopathology. A BMDExpress algorithm using Gene Ontology Biological Process gene sets was optimized to derive a liver BEPOD most predictive of a systemic apical POD. Liver BEPODs were stable from 3 h to 29 days of exposure; the median fold difference of the 29-day BEPOD to BEPODs from earlier time points was approximately 1 (range: 0.7-1.1). Strong positive correlation (Pearson R = 0.86) and predictive accuracy (root mean square difference = 0.41) were observed between a concurrent (29 days) liver BEPOD and the systemic apical POD. Similar Pearson R and root mean square difference values were observed for comparisons between a 29-day systemic apical POD and liver BEPODs derived from 3 h to 15 days of exposure. These data across 79 molecules suggest that a longer-term exposure study apical POD from liver and non-liver compartments can be estimated using a liver BEPOD derived from an acute or subacute exposure study.

Application of transcriptomic and microRNA profiling in the evaluation of potential liver carcinogens

Hepatocarcinogens are agents that increase the incidence of liver cancer in exposed animals or humans. It is now established that carcinogenic exposures have a widespread impact on the transcriptome, inducing both adaptive and adverse changes in the activities of genes and pathways. Chemical hepatocarcinogens have also been shown to affect expression of microRNA (miRNA), the evolutionarily conserved noncoding RNA that regulates gene expression posttranscriptionally. Considerable effort has been invested into examining the involvement of mRNA in chemical hepatocarcinogenesis and their potential usage for the classification and prediction of new chemical entities. For miRNA, there has been an increasing number of studies reported over the past decade, although not to the same degree as for transcriptomic studies. Current data suggest that it is unlikely that any gene or miRNA signature associated with short-term carcinogen exposure can replace the rodent bioassay. In this review, we discuss the application of transcriptomic and miRNA profiles to increase mechanistic understanding of chemical carcinogens and to aid in their classification.

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.