Case study on the utility of hepatic global gene expression profiling in the risk assessment of the carcinogen furan

Epigenotoxicity: Decoding the epigenetic imprints of genotoxic agents and their implications for regulatory genetic toxicology

Regulatory genetic toxicology focuses on DNA damage and subsequent gene mutations. However, genotoxic agents can also affect epigenetic marks, and incorporation of epigenetic data into the regulatory framework may thus enhance the accuracy of risk assessment. Additionally, epigenetic alterations may identify non-genotoxic carcinogens that are not captured with the current battery of tests. Epigenetic alterations could also explain long-term consequences and potential transgenerational effects in the absence of DNA mutations. Therefore, at the 2022 International Workshops on Genotoxicity Testing (IWGT) in Ottawa (Ontario, Canada), an expert workgroup explored whether including epigenetic endpoints would improve regulatory genetic toxicology. Here we summarize the presentations and the discussions on technical advancements in assessing epigenetics, how the assessment of epigenetics can enhance strategies for detecting genotoxic and non-genotoxic carcinogens and the correlation between epigenetic alterations with other relevant apical endpoints.

Adverse outcome pathway exploration of furan-induced liver fibrosis in rats: Genotoxicity pathway or oxidative stress pathway through CYP2E1 activation?

Furan is a widespread endogenous contaminant in heat-processed foods that can accumulate rapidly in the food chain and has been widely detected in foods, such as wheat, bread, coffee, canned meat products, and baby food. Dietary exposure to this chemical may bring health risk. Furan is classified as a possible category 2B human carcinogen by the International Agency for Research on Cancer, with the liver as its primary target organ. Hepatic fibrosis is the most important nontumoral harmful effect of furan and also an important event in the carcinogenesis of furan. Although the specific mechanism of furan-induced liver fibrosis is still unclear, it may involve oxidative stress and genetic toxicity, in which the activation of cytochrome P450 2E1 (CYP2E1) may be the key event. Thus, we conducted a study using an integrating multi-endpoint genotoxicity platform in 120-day in vivo subchronic toxicity test in rats. Results showed that the rats with activated CYP2E1 exhibited DNA double-strand breaks in D4, gene mutations in D60, and increased expression of reactive oxygen species and nuclear factor erythroid 2-related factor 2 in D120. Necrosis, apoptosis, hepatic stellate cell activation, and fibrosis also occurred in the liver, suggesting that furan can independently affect liver fibrosis through oxidative stress and genotoxicity pathways. Point of Departure (PoD) was obtained by benchmark-dose (BMD) method to establish health-based guidance values. The human equivalent dose of PoD derived from BMDL05 was 2.26 μg/kg bw/d. The findings laid a foundation for the safety evaluation and risk assessment of furan and provided data for the further construction and improvement of the adverse outcome pathway network in liver fibrosis.

Early microRNA responses in rodent liver mediated by furan exposure establish dose thresholds for later adverse outcomes

To reduce reliance on long-term in vivo studies, short-term data linking early molecular-based measurements to later adverse health effects is needed. Although transcriptional-based benchmark dose (BMDT) modeling has been used to estimate potencies and stratify chemicals based on potential to induce later-life effects, dose-responsive epigenetic alterations have not been routinely considered. Here, we evaluated the utility of microRNA (miRNA) profiling in mouse liver and blood, as well as in mouse primary hepatocytes in vitro, to indicate mechanisms of liver perturbation due to short-term exposure of the known rodent liver hepatotoxicant and carcinogen, furan. Benchmark dose modeling of miRNA measurements (BMDmiR) were compared to the referent transcriptional (BMDT) and apical (BMDA) estimates. These analyses indicate a robust dose response for 34 miRNAs to furan and involvement of p53-linked pathways in furan-mediated hepatotoxicity, supporting mRNA and apical measurements. Liver-sourced miRNAs were also altered in the blood and primary hepatocytes. Overall, these results indicate mechanistic involvement of miRNA in furan carcinogenicity and provide evidence of their potential utility as accessible biomarkers of exposure and disease.

Toxicity, genotoxicity, and carcinogenicity of 2-methylfuran in a 90-day comprehensive toxicity study in gpt delta rats

2-Methylfuran (2-MF) exists naturally in foods and is used as a flavoring agent. Furan, the core structure of 2-MF, possesses hepatocarcinogenicity in rodents. Accumulation of toxicological information on furan derivatives is needed to elucidate their carcinogenic mode of action. In the current study, we examined the comprehensive toxicological studies of 2-MF using gpt delta rats. 2-MF was intragastrically administered to groups of 10 male and 10 female Sprague-Dawley gpt delta rats at a dose of 0, 1.2, 6, or 30 mg/kg/day for 13 weeks. Effects of 2-MF on the hepatobiliary system including an increase in serum alkaline phosphatase were observed in the 6 and 30 mg/kg groups, and cholangiofibrosis was found in the 30 mg/kg group. The no observed adverse effect level was set at 1.2 mg/kg/day for both sexes and 1.14 mg/kg/day was determined as the benchmark dose low. The acceptable daily intake was calculated to be 11.4 μg/kg/day. Increases in the number and areas of glutathione S-transferase placental form-positive foci in the 30 mg/kg group were apparent, suggesting the hepatocarcinogenicity of 2-MF in rats. By contrast, the lack of increase in in vivo mutagenicity in the liver implied that 2-MF hepatocarcinogenesis may not involve genotoxic mechanisms.

Evaluating the Influences of Confounding Variables on Benchmark Dose using a Case Study in the Field of Ionizing Radiation

Purpose A vast amount of data regarding the effects of radiation stressors on transcriptional changes has been produced over the past few decades. These data have shown remarkable consistency across platforms and experimental design, enabling increased understanding of early molecular effects of radiation exposure. However, the value of transcriptomic data in the context of risk assessment is not clear and represents a gap that is worthy of further consideration. Recently, benchmark dose (BMD) modeling has shown promise in correlating a transcriptional point of departure (POD) to that derived using phenotypic outcomes relevant to human health risk assessment. Although frequently applied in chemical toxicity evaluation, our group has recently demonstrated application within the field of radiation research. This approach allows the possibility to quantitatively compare radiation-induced gene and pathway alterations across various datasets using BMD values and derive meaningful biological effects. However, before BMD modeling can confidently be used, an understanding of the impact of confounding variables on BMD outputs is needed. Methods: To this end, BMD modeling was applied to a publicly available microarray dataset (Gene Expression Omnibus #GSE23515) that used peripheral blood ex-vivo gamma-irradiated at 0.82 Gy/min, at doses of 0, 0.1, 0.5 or 2 Gy, and assessed 6 hours post-exposure. The dataset comprised six female smokers (F-S), six female non-smokers (F-NS), six male smokers (M-S), and six male non-smokers (M-NS). Results: A combined total of 412 genes were fit to models and the BMD distribution was noted to be bi-modal across the four groups. A total of 74, 41, 62 and 62 genes were unique to the F-NS, M-NS, F-S and M-S groups. Sixty-two BMD modeled genes and nine pathways were common across all four groups. There were no differential sensitivity of responses in the robust common genes and pathways. Conclusion: For radiation-responsive genes and pathways common across the study groups, the BMD distribution of transcriptional activity was unaltered by sex and smoking status. Although further validation of the data is needed, these initial findings suggest BMD values for radiation relevant genes and pathways are robust and could be explored further in future studies.

Glycidamide and cis-2-butene-1,4-dial (BDA) as potential carcinogens and promoters of liver cancer - An in vitro study

Acrylamide and furan are environmental and food contaminants that are metabolized by cytochrome P450 2E1 (CYP2E1), giving rise to glycidamide and cis-2-butene-1,4-dial (BDA) metabolites, respectively. Both glycidamide and BDA are electrophilic species that react with nucleophilic groups, being able to introduce mutations in DNA and perform epigenetic remodeling. However, whereas these carcinogens are primarily metabolized in the liver, the carcinogenic potential of acrylamide and furan in this organ is still controversial, based on findings from experimental animal studies. With the ultimate goal of providing further insights into this issue, we explored in vitro, using a hepatocyte cell line and a hepatocellular carcinoma cell line, the putative effect of these metabolites as carcinogens and cancer promoters. Molecular alterations were investigated in cells that survive glycidamide and BDA toxicity. We observed that those cells express CD133 stemness marker, present a high proliferative capacity and display an adjusted expression profile of genes encoding enzymes involved in oxidative stress control, such as GCL-C, GSTP1, GSTA3 and CAT. These molecular changes seem to be underlined, at least in part, by epigenetic remodeling involving histone deacetylases (HDACs). Although more studies are needed, here we present more insights towards the carcinogenic capacity of glycidamide and BDA and also point out their effect in favoring hepatocellular carcinoma progression.

A microRNA or messenger RNA point of departure estimates an apical endpoint point of departure in a rat developmental toxicity model

Traditional developmental toxicity testing practice examines fetal apical endpoints to identify a point of departure (POD) for risk assessment. A potential new testing paradigm involves deriving a POD from a comprehensive analysis of molecular-level change. Here, the rat ketoconazole endocrine-mediated developmental toxicity model was used to test the hypothesis that maternal epigenomic (miRNA) and transcriptomic (mRNA) PODs are similar to fetal apical endpoint PODs. Sprague-Dawley rats were exposed from gestation day (GD) 6-21 to 0, 0.063, 0.2, 0.63, 2, 6.3, 20, or 40 mg/kg/day ketoconazole. Dam systemic, liver, and placenta PODs, along with GD 21 fetal resorption, body weight, and skeletal apical PODs were derived using BMDS software. GD 21 dam liver and placenta miRNA and mRNA PODs were obtained using three methods: a novel individual molecule POD accumulation method, a first mode method, and a gene set method. Dam apical POD values ranged from 2.0 to 38.6 mg/kg/day; the lowest value was for placenta histopathology. Fetal apical POD values were 10.9-20.3 mg/kg/day; the lowest value was for fetal resorption. Dam liver miRNA and mRNA POD values were 0.34-0.69 mg/kg/day, and placenta miRNA and mRNA POD values were 2.53-6.83 mg/kg/day. Epigenomic and transcriptomic POD values were similar across liver and placenta. Deriving a molecular POD from dam liver or placenta was protective of a fetal apical POD. These data support the conclusion that a molecular POD can be used to estimate, or be protective of, a developmental toxicity apical POD.

Benchmark Dose Modeling Approaches for Volatile Organic Chemicals using a Novel Air-Liquid Interface In Vitro Exposure System

Inhalation is the most relevant route of volatile organic chemical (VOC) exposure; however, due to unique challenges posed by their chemical properties and poor solubility in aqueous solutions, in vitro chemical safety testing is predominantly performed using direct application dosing/submerged exposures. To address the difficulties in screening toxic effects of VOCs, our cell culture exposure system permits cells to be exposed to multiple concentrations at air-liquid interface (ALI) in a 24-well format. ALI exposure methods permit direct chemical-to-cell interaction with the test article at physiological conditions. In the present study, BEAS-2B and primary normal human bronchial epithelial cells (pHBEC) are used to assess gene expression, cytotoxicity, and cell viability responses to a variety of volatile chemicals including acrolein, formaldehyde, 1,3-butadiene, acetaldehyde, 1-bromopropane, carbon tetrachloride, dichloromethane, and trichloroethylene. BEAS-2B cells were exposed to all the test agents, while pHBECs were only exposed to the latter four listed above. The VOC concentrations tested elicited only slight cell viability changes in both cell types. Gene expression changes were analyzed using benchmark dose (BMD) modeling. The BMD for the most sensitive gene set was within one order of magnitude of the threshold-limit value reported by the American Conference of Governmental Industrial Hygienists, and the most sensitive gene sets impacted by exposure correlate to known adverse health effects recorded in epidemiologic and in vivo exposure studies. Overall, our study outlines a novel in vitro approach for evaluating molecular-based points-of-departure in human airway epithelial cell exposure to volatile chemicals.

Salidroside alleviates liver inflammation in furan-induced mice by regulating oxidative stress and endoplasmic reticulum stress

Furan is a genotoxic and carcinogenic toxicant formed during the food thermal processing. Our previous studies confirmed that salidroside (SAL) displayed excellent protective effects against furan-induced hepatotoxicity and inflammation, whereas the underlying mechanism was still unclear. In the current study, Balb/c mice were divided to the control group (CON), the furan model group (FUR8, 8 mg/kg BW furan for 30 days) and SAL intervention groups (SAL10/20/40, 8 mg/kg BW furan for 30 days + 10/20/40 mg/kg BW SAL from day 16 to day 30). The alleviative effects and the mechanisms of SAL against furan-induced liver inflammation in mice were investigated through oxidative stress (OS) and endoplasmic reticulum stress (ERS). Liver metabonomics data, molecular docking and Western-blotting results implied that SAL suppressed the activity and the high expression of hepatic CYP2E1, and alleviated liver OS induced by furan. Levels of key markers (GRP78, CHOP and Caspase-12) of ERS and proteins in IRE1α pathway of the UPR branch increased by furan were prominently reduced after SAL treatment. Levels of phosphorylated proteins JNK, ERK, p38, IKKα/β, IκB and p65 in MAPK and NF-κB pathways were also suppressed by SAL. We further confirmed that SAL inhibited furan-induced inflammation by reducing the levels of NLRP3, ASC, Cleaved Caspase-1 and IL-1β and decreasing the production of pro-inflammatory cytokines. Our results shed light into the alleviating mechanisms behind furan-induced liver inflammation, and suggested that SAL inhibited OS, ERS and related MAPK and NF-κB pathways and therefore inhibited the NLRP3 inflammasome activation, which may be its potential mechanism of alleviating liver inflammation.

Meta-analysis of transcriptomic datasets using benchmark dose modeling shows value in supporting radiation risk assessment

PURPOSE

Benchmark dose (BMD) modeling is used to determine the dose of a stressor at which a predefined increase in any biological effect above background occurs (e.g. 10% increase from control values). BMD analytical tools have the capacity to model transcriptional dose-response data to derive BMDs for genes, pathways and gene ontologies. We recently demonstrated the value of this approach to support various areas of radiation research using predominately 'in-house' generated datasets.

MATERIALS AND METHODS

As a continuation of this work, transcriptomic studies of relevance to ionizing radiation were retrieved through the Gene Expression Omnibus (GEO). The datasets were compiled and filtered, then analyzed using BMDExpress. The objective was to determine the reproducibility of BMD values in relation to pathways and genes across different exposure scenarios and compare to those derived using cytogenetic endpoints. A number of graphic visualization approaches were used to determine if BMD outputs could be correlated to parameters such as dose-rate, radiation quality and cell type.

RESULTS

Curated studies were diverse and derived from experiments with varied design and intent. Despite this, common genes and pathways were identified with low and high dose thresholds. The higher BMD values were associated with immune response and cell death, while transcripts with lower BMD values were generally related to the classic DNA damage response/repair processes, centered on TP53 signaling. Analysis of datasets with relatively similar dose-ranges under comparable experimental conditions showed a bi-modal distribution with a high degree of consistency in BMD values across shared genes and pathways, particularly for those below the 25th percentile of total distribution by dose. The median BMD values were noted to be approximately 0.5 Gy for genes/pathways that comprised mode 1. Furthermore, transcriptional BMD values derived from a subset of genes using in vivo and in vitro datasets were in accord to those using cytogenetic endpoints.

CONCLUSION

Overall, the results from this work highlight the value of the BMD methodology to derive meaningful outputs that are consistent across different models, provided the studies are conducted using a similar dose-range.

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.

Transcriptomics in Toxicogenomics, Part I: Experimental Design, Technologies, Publicly Available Data, and Regulatory Aspects

The starting point of successful hazard assessment is the generation of unbiased and trustworthy data. Conventional toxicity testing deals with extensive observations of phenotypic endpoints in vivo and complementing in vitro models. The increasing development of novel materials and chemical compounds dictates the need for a better understanding of the molecular changes occurring in exposed biological systems. Transcriptomics enables the exploration of organisms' responses to environmental, chemical, and physical agents by observing the molecular alterations in more detail. Toxicogenomics integrates classical toxicology with omics assays, thus allowing the characterization of the mechanism of action (MOA) of chemical compounds, novel small molecules, and engineered nanomaterials (ENMs). Lack of standardization in data generation and analysis currently hampers the full exploitation of toxicogenomics-based evidence in risk assessment. To fill this gap, TGx methods need to take into account appropriate experimental design and possible pitfalls in the transcriptomic analyses as well as data generation and sharing that adhere to the FAIR (Findable, Accessible, Interoperable, and Reusable) principles. In this review, we summarize the recent advancements in the design and analysis of DNA microarray, RNA sequencing (RNA-Seq), and single-cell RNA-Seq (scRNA-Seq) data. We provide guidelines on exposure time, dose and complex endpoint selection, sample quality considerations and sample randomization. Furthermore, we summarize publicly available data resources and highlight applications of TGx data to understand and predict chemical toxicity potential. Additionally, we discuss the efforts to implement TGx into regulatory decision making to promote alternative methods for risk assessment and to support the 3R (reduction, refinement, and replacement) concept. This review is the first part of a three-article series on Transcriptomics in Toxicogenomics. These initial considerations on Experimental Design, Technologies, Publicly Available Data, Regulatory Aspects, are the starting point for further rigorous and reliable data preprocessing and modeling, described in the second and third part of the review series.

Toxicological effects of furan on the reproductive system of male rats: An "in vitro" and "in vivo"-based endocrinological and spermatogonial study

Furan is a colorless toxic chemical produced in various food items during heat processing and in chemical industries. Both in vitro and in vivo studies have reported that it induces oxidative stress and endocrine disruption; however, limited data are available regarding the effects of furan on the reproduction of mammals. In the present study, an in vitro experiment was designed to investigate the direct effects of furan exposure on oxidative stress and testosterone concentration in rat testicular tissue. Furan not only generated high oxidative stress but also decreased antioxidant enzyme activity in the testicular tissue. On the basis of in vitro study results, an in vivo sub-chronic exposure study was performed. Male rats were orally exposed to different concentrations of furan (0, 5, 10, 20, and 40 mg kg^-1). An increase (P < 0.05) of reactive oxygen species (ROS) and of the lipid profile (cholesterol, triglycerides, and LDL) in higher dose treatment groups of furan was observed, while total protein content and antioxidant enzyme activity were considerably decreased after furan exposure. Also, plasma and intratesticular testosterone concentrations were reduced in high-dose treatment groups. Sperm parameters such as sperm viability, sperm count, and sperm motility showed a decrease (P < 0.05) in a dose-dependent manner. Histopathological findings revealed significant alterations in testis and epididymis tissues. These results confirm that furan can induce toxic effects on the reproductive system of male rats.

Evaluation of Polymer Nanoformulations in Hepatoma Therapy by Established Rodent Models

Hepatoma is one of the most severe malignancies usually with poor prognosis, and many patients are insensitive to the existing therapeutic agents, including the drugs for chemotherapy and molecular targeted therapy. Currently, researchers are committed to developing the advanced formulations with controlled drug delivery to improve the efficacy of hepatoma therapy. Numerous inoculated, induced, and genetically engineered hepatoma rodent models are now available for formulation screening. However, animal models of hepatoma cannot accurately represent human hepatoma in terms of histological characteristics, metastatic pathways, and post-treatment responses. Therefore, advanced animal hepatoma models with comparable pathogenesis and pathological features are in urgent need in the further studies. Moreover, the development of nanomedicines has renewed hope for chemotherapy and molecular targeted therapy of advanced hepatoma. As one kind of advanced formulations, the polymer-based nanoformulated drugs have many advantages over the traditional ones, such as improved tumor selectivity and treatment efficacy, and reduced systemic side effects. In this article, the construction of rodent hepatoma model and much information about the current development of polymer nanomedicines were reviewed in order to provide a basis for the development of advanced formulations with clinical therapeutic potential for hepatoma.

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

Hexabromocyclododecane (HBCD) is a brominated flame retardant found in the environment and human tissues. The toxicological effects of HBCD exposure are not clearly understood. We employed whole-genome RNA-sequencing on liver samples from male and female Fischer rats exposed to 0, 250, 1250, and 5000 mg technical mixture of HBCD/kg diet for 28 days to gain further insight into HBCD toxicity. HBCD altered 428 and 250 gene transcripts in males and females, respectively, which were involved in metabolism of xenobiotics, oxidative stress, immune response, metabolism of glucose and lipids, circadian regulation, cell cycle, fibrotic activity, and hormonal balance. Signature analysis supported that HBCD operates through the constitutive androstane and pregnane X receptors. The median transcriptomic benchmark dose (BMD) for the lowest statistically significant pathway was within 1.5-fold of the BMD for increased liver weight, while the BMD for the lowest pathway with at least three modeled genes (minimum 5% of pathway) was similar to the lowest apical endpoint BMD. The results show how transcriptional analyses can inform mechanisms underlying chemical toxicity and the doses at which potentially adverse effects occur. This experiment is part of a larger study exploring the use of toxicogenomics and high-throughput screening for human health risk assessment.

Furan-induced cardiotoxicity in diabetic rats and protective role of lycopene

The objective of this current study is to search the impacts of furan and lycopene on a diabetic rat's heart. Diabetes increases prevalence with a sedentary lifestyle and obesity. Furan, a carcinogen, was detected in foods that had undergone thermal treatment. Two groups were formed, such as: the control and diabetic groups. Diabetic groups fall into four allocated groups: control, lycopene, furan, and furan + lycopene. The changes in malondialdehide levels, antioxidant enzymes' activities, and histopathology of the heart were pointed out. We observed that the changes in both the MDA level and the antioxidant enzymes' activities were due to diabetes. Furan-induced increment in MDA levels, but GST, CAT, SOD, and GPx activities showed a decrease. Lycopene influenced these changes positively. In terms of the pathological studies, changes were also observed in diabetic rats. The histological damage in the diabetic furan groups was found to be very serious. Lycopene was protective against furan that caused histopathological changes. Diabetes and furan resulted into heart damages in rats and lycopene showed preventive effects. But this was not a complete protection. PRACTICAL APPLICATIONS: Furan by product of chlorinated chemical compounds production and chemical processes including combustion, are virtually ubiquitous in the environment. Because of this, furan occur potential risk for human health. Furan is known to exist at higher levels in fatty foods, such as full-fat milk and dairy products, meat and eggs. They can accumulate in the fatty tissues of animals and humans. Diabetes increase in the prevalence of a sedentary lifestyle and obesity. Diabetes with furan-induced cardiotoxicity. Foods should not be heated over and over again, especially diabetic heart patients should stay away from this situation.

The application of transcriptional benchmark dose modeling for deriving thresholds of effects associated with solar-simulated ultraviolet radiation exposure

Considerable data has been generated to elucidate the transcriptional response of cells to ultraviolet radiation (UVR) exposure providing a mechanistic understanding of UVR-induced cellular responses. However, using these data to support standards development has been challenging. In this study, we apply benchmark dose (BMD) modeling of transcriptional data to derive thresholds of gene responsiveness following exposure to solar-simulated UVR. Human epidermal keratinocytes were exposed to three doses (10, 20, 150 kJ/m² ) of solar simulated UVR and assessed for gene expression changes 6 and 24 hr postexposure. The dose-response curves for genes with p-fit values (≥ 0.1) were used to derive BMD values for genes and pathways. Gene BMDs were bi-modally distributed, with a peak at ∼16 kJ/m² and ∼108 kJ/m² UVR exposure. Genes/pathways within Mode 1 were involved in cell signaling and DNA damage response, while genes/pathways in the higher Mode 2 were associated with immune response and cancer development. The median value of each Mode coincides with the current human exposure limits for UVR and for the minimal erythemal dose, respectively. Such concordance implies that the use of transcriptional BMD data may represent a promising new approach for deriving thresholds of actinic effects. Environ. Mol. Mutagen. 59:502-515, 2018. © 2018 The Authors Environmental and Molecular Mutagenesis published by Wiley Periodicals, Inc. on behalf of Environmental Mutagen Society.

Editor's Highlight: Application of Gene Set Enrichment Analysis for Identification of Chemically Induced, Biologically Relevant Transcriptomic Networks and Potential Utilization in Human Health Risk Assessment

The rate of new chemical development in commerce combined with a paucity of toxicity data for legacy chemicals presents a unique challenge for human health risk assessment. There is a clear need to develop new technologies and incorporate novel data streams to more efficiently inform derivation of toxicity values. One avenue of exploitation lies in the field of transcriptomics and the application of gene expression analysis to characterize biological responses to chemical exposures. In this context, gene set enrichment analysis (GSEA) was employed to evaluate tissue-specific, dose-response gene expression data generated following exposure to multiple chemicals for various durations. Patterns of transcriptional enrichment were evident across time and with increasing dose, and coordinated enrichment plausibly linked to the etiology of the biological responses was observed. GSEA was able to capture both transient and sustained transcriptional enrichment events facilitating differentiation between adaptive versus longer term molecular responses. When combined with benchmark dose (BMD) modeling of gene expression data from key drivers of biological enrichment, GSEA facilitated characterization of dose ranges required for enrichment of biologically relevant molecular signaling pathways, and promoted comparison of the activation dose ranges required for individual pathways. Median transcriptional BMD values were calculated for the most sensitive enriched pathway as well as the overall median BMD value for key gene members of significantly enriched pathways, and both were observed to be good estimates of the most sensitive apical endpoint BMD value. Together, these efforts support the application of GSEA to qualitative and quantitative human health risk assessment.

Risks for public health related to the presence of furan and methylfurans in food

The European Commission asked EFSA for a scientific evaluation on the risk to human health of the presence of furan and methylfurans (2-methylfuran, 3-methylfuran and 2,5-dimethylfuran) in food. They are formed in foods during thermal processing and can co-occur. Furans are produced from several precursors such as ascorbic acid, amino acids, carbohydrates, unsaturated fatty acids and carotenoids, and are found in a variety of foods including coffee and canned and jarred foods. Regarding furan occurrence, 17,056 analytical results were used in the evaluation. No occurrence data were received on methylfurans. The highest exposures to furan were estimated for infants, mainly from ready-to-eat meals. Grains and grain-based products contribute most for toddlers, other children and adolescents. In adults, elderly and very elderly, coffee is the main contributor to dietary exposure. Furan is absorbed from the gastrointestinal tract and is found in highest amounts in the liver. It has a short half-life and is metabolised by cytochrome P450 2E1 (CYP2E1) to the reactive metabolite, cis-but-2-ene-1,4-dialdehyde (BDA). BDA can bind covalently to amino acids, proteins and DNA. Furan is hepatotoxic in rats and mice with cholangiofibrosis in rats and hepatocellular adenomas/carcinomas in mice being the most prominent effects. There is limited evidence of chromosomal damage in vivo and a lack of understanding of the underlying mechanism. Clear evidence for indirect mechanisms involved in carcinogenesis include oxidative stress, gene expression alterations, epigenetic changes, inflammation and increased cell proliferation. The CONTAM Panel used a margin of exposure (MOE) approach for the risk characterisation using as a reference point a benchmark dose lower confidence limit for a benchmark response of 10% of 0.064 mg/kg body weight (bw) per day for the incidence of cholangiofibrosis in the rat. The calculated MOEs indicate a health concern. This conclusion was supported by the calculated MOEs for the neoplastic effects.

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

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

Integration of the TGx-28.65 genomic biomarker with the flow cytometry micronucleus test to assess the genotoxicity of disperse orange and 1,2,4-benzenetriol in human TK6 cells

In vitro gene expression signatures to predict toxicological responses can provide mechanistic context for regulatory testing. We previously developed the TGx-28.65 genomic biomarker from a database of gene expression profiles derived from human TK6 cells exposed to 28 well-known compounds. The biomarker comprises 65 genes that can classify chemicals as DNA damaging or non-DNA damaging. In this study, we applied the TGx-28.65 genomic biomarker in parallel with the in vitro micronucleus (MN) assay to determine if two chemicals of regulatory interest at Health Canada, disperse orange (DO: the orange azo dye 3-[[4-[(4-Nitrophenyl)azo]phenyl] benzylamino]propanenitrile) and 1,2,4-benzenetriol (BT: a metabolite of benzene) are genotoxic or non-genotoxic. Both chemicals caused dose-dependent declines in relative survival and increases in apoptosis. A strong significant increase in MN induction was observed for all concentrations of BT; the top two concentrations of DO also caused a statistically significant increase in MN, but these increases were <2-fold above controls. TGx-28.65 analysis classified BT as genotoxic at all three concentrations and DO as genotoxic at the mid and high concentrations. Thus, although DO only caused a small increase in MN, this response was sufficient to induce a cellular DNA damage response. Benchmark dose modeling confirmed that BT is much more potent than DO. The results strongly suggest that follow-up work is required to assess whether DO and BT are also genotoxic in vivo. This is particularly important for DO, which may require metabolic activation by bacterial gut flora to fully induce its genotoxic potential. Our previously published data and this proof of concept study suggest that the TGx-28.65 genomic biomarker has the potential to add significant value to existing approaches used to assess genotoxicity.

High-Throughput Screening Data Interpretation in the Context of In Vivo Transcriptomic Responses to Oral Cr(VI) Exposure

The toxicity of hexavalent chromium [Cr(VI)] in drinking water has been studied extensively, and available in vivo and in vitro studies provide a robust dataset for application of advanced toxicological tools to inform the mode of action (MOA). This study aimed to contribute to the understanding of Cr(VI) MOA by evaluating high-throughput screening (HTS) data and other in vitro data relevant to Cr(VI), and comparing these findings to robust in vivo data, including transcriptomic profiles in target tissues. Evaluation of Tox21 HTS data for Cr(VI) identified 11 active assay endpoints relevant to the Ten Key Characteristics of Carcinogens (TKCCs) that have been proposed by other investigators. Four of these endpoints were related to TP53 (tumor protein 53) activation mapping to genotoxicity (KCC#2), and four were related to cell death/proliferation (KCC#10). HTS results were consistent with other in vitro data from the Comparative Toxicogenomics Database. In vitro responses were compared to in vivo transcriptomic responses in the most sensitive target tissue, the duodenum, of mice exposed to ≤ 180 ppm Cr(VI) for 7 and 90 days. Pathways that were altered both in vitro and in vivo included those relevant to cell death/proliferation. In contrast, pathways relevant to p53/DNA damage were identified in vitro but not in vivo. Benchmark dose modeling and phenotypic anchoring of in vivo transcriptomic responses strengthened the finding that Cr(VI) causes cell stress/injury followed by proliferation in the mouse duodenum at high doses. These findings contribute to the body of evidence supporting a non-mutagenic MOA for Cr(VI)-induced intestinal cancer.

Recommended approaches in the application of toxicogenomics to derive points of departure for chemical risk assessment

There is increasing interest in the use of quantitative transcriptomic data to determine benchmark dose (BMD) and estimate a point of departure (POD) for human health risk assessment. Although studies have shown that transcriptional PODs correlate with those derived from apical endpoint changes, there is no consensus on the process used to derive a transcriptional POD. Specifically, the subsets of informative genes that produce BMDs that best approximate the doses at which adverse apical effects occur have not been defined. To determine the best way to select predictive groups of genes, we used published microarray data from dose–response studies on six chemicals in rats exposed orally for 5, 14, 28, and 90 days. We evaluated eight approaches for selecting genes for POD derivation and three previously proposed approaches (the lowest pathway BMD, and the mean and median BMD of all genes). The relationship between transcriptional BMDs derived using these 11 approaches and PODs derived from apical data that might be used in chemical risk assessment was examined. Transcriptional BMD values for all 11 approaches were remarkably aligned with corresponding apical PODs, with the vast majority of toxicogenomics PODs being within tenfold of those derived from apical endpoints. We identified at least four approaches that produce BMDs that are effective estimates of apical PODs across multiple sampling time points. Our results support that a variety of approaches can be used to derive reproducible transcriptional PODs that are consistent with PODs produced from traditional methods for chemical risk assessment.

A microchip for integrated single-cell genotoxicity assay

With the development of large-scale biologic databases, precision medicine is becoming a frontier in biomedical research. As a main focus of precision medicine study, cancer has been widely accepted as a disease born out of inherited genetic variations or accumulating genomic damage. At the single-cell level, microfluidics or lab-on-a-chip technology for cancer study is an emerging tool for improving risk assessment, diagnostic categories and therapeutic strategies. This work presents a multi-layer microchip for single-cell gene expression profiling. Treated by three drug reagents (i.e. methyl methanesulfonate, docetaxel and colchicine) with varied concentrations and time lengths, individual human breast cancer cells (MCF-7) are then lysed on-chip, and the released mRNA templates are captured and reversely transcribed into cDNA on microbead surface. Three genes (GAPDH, CDKN1A, AURKA) are amplified and quantified simultaneously through triplex real-time polymerase chain reactions (qPCR). Readout per run is set to be eighteen, and can be further improved following same approach. The microchip is able to integrate all steps of single-cell gene expression profiling, and provide precision study of drug induced genotoxicity with reduced reagents consumption per reaction and instrumental cost.

Furan-induced hepatotoxic and hematologic changes in diabetic rats: the protective role of lycopene

Furan forms as a result of thermal treatment of food and induces harmful effects on organisms. In our work, lycopene, furan, and a combination of the two were given to diabetic male rats for 28 days. Hematological changes, total protein and cholesterol, triglyceride, and albumin levels, alanine aminotransferase, aspartate aminotransferase, lactate dehydrogenase, and alkaline phosphatase activities of the serum, malondialdehyde levels, glutathione peroxidase, catalase, glutathione-S-transferase, superoxide dismutase activities, DNA damage in liver tissues and hepatic histopathological alterations were compared to a control group. There were significant changes in the liver function tests, DNA damage, activities of antioxidant enzymes, and malondialdehyde levels between diabetic control and non-diabetic control groups, between diabetic control and diabetic lycopene groups, and also between diabetic furan and diabetic control groups. In diabetic lycopene and diabetic furan + lycopene treated groups we designated the preventive effects of lycopene against diabetes and furan, however, on the analysed parameters only. In spite of some pathological alterations designated in diabetic furan treated group’s liver, fewer pathological alterations were observed in furan+lycopene treated groups at the end of week 4. Consequently, lycopene significantly reduced furan- and diabetes-induced toxicity in rat liver.

Transcriptional benchmark dose modeling: Exploring how advances in chemical risk assessment may be applied to the radiation field

Recent advances in "-omics" technologies have simplified capacity to concurrently assess expression profiles of thousands of targets in a cellular system. However, compilation and analysis of "omics" data in support of human health protection remains a challenge. Benchmark dose (BMD) modeling is currently being employed in chemical risk assessment to estimate acceptable levels of exposure. Although typically applied to conventional endpoints, newer software has enabled this application to be extended to transcriptomic datasets. BMD analytical tools now have the capacity to model transcriptional dose-response data to derive meaningful BMD values for genes, pathways and gene ontologies. In this report, radiation data obtained from the Gene Expression Omnibus (GEO) were analyzed to generate BMD values for transcriptional responses. The datasets comprised microarray analyses of human blood gamma-irradiated ex vivo (0-20 Gy) and human-derived cell lines exposed to alpha particle radiation (0.5-1.5 Gy). The distributions of BMDs for statistically significant genes and pathways in response to radiation exposure were examined and compared across studies. BMD modeling could identify pathway/gene sensitivities across wide radiation dose ranges, experimental conditions (time-points, cell types) and radiation qualities. BMD analysis offered a new approach to examine transcriptional data. The results were shown to provide information on transcriptional thresholds of effects to support refined risk assessments for low dose ionizing radiation exposures, derive gene-based values for relative biological effectiveness and identify pathways involved in radiation sensitivities across cell types which may extend to applications a clinical setting. Environ. Mol. Mutagen. 57:589-604, 2016. © 2016 Wiley Periodicals, Inc.

Toxicogenomics Case Study: Furan

Development of pragmatic methodologies for human health risk assessment is required to address current regulatory challenges. We applied three toxicogenomic approaches—quantitative, predictive, and mechanistic—to a case study in mice exposed for 3 weeks to the hepatocarcinogen furan. We modeled the dose response of a variety of transcriptional endpoints and found that they produced benchmark doses similar to the furan-dependent cancer benchmark doses. Meta-analyses showed strong similarity between furan-dependent gene expression changes and those associated with several hepatic pathologies. Molecular pathways facilitated the development of a molecular mode of action for furan-induced hepatocellular carcinogenicity. Finally, we compared transcriptomic profiles derived from formalin-fixed and paraffin-embedded (FFPE) samples with those from high-quality frozen samples to evaluate whether archival samples are a viable option for toxicogenomic studies. The advantage of using FFPE tissues is that they are very well characterized (phenotypically); the disadvantage is that formalin degrades biomacromolecules, including RNA. We found that FFPE samples can be used for toxicogenomics using a ribo-depletion RNA-seq protocol. Our case study demonstrates the utility of toxicogenomics data to human health risk assessment, the potential of archival FFPE tissue samples, and identifies viable strategies toward the reduction of animal usage in chemical testing.

Furan-induced transcriptomic and gene-specific DNA methylation changes in the livers of Fischer 344 rats in a 2-year carcinogenicity study

Furan is a significant food contaminant and a potent hepatotoxicant and rodent liver carcinogen. The carcinogenic effect of furan has been attributed to genotoxic and non-genotoxic, including epigenetic, changes in the liver; however, the mechanisms of the furan-induced liver tumorigenicity are still unclear. The goal of the present study was to investigate the role of transcriptomic and epigenetic events in the development of hepatic lesions in Fischer (F344) rats induced by furan treatment in a classic 2-year rodent tumorigenicity bioassay. High-throughput whole-genome transcriptomic analysis demonstrated distinct alterations in gene expression in liver lesions induced in male F344 rats treated with 0.92 or 2.0 mg furan/kg body weight (bw)/day for 104 weeks. Compared to normal liver tissue, 1336 and 1541 genes were found to be differentially expressed in liver lesions in rats treated with 0.92 and 2.0 mg furan/kg bw/day, respectively, among which 1001 transcripts were differentially expressed at both doses. Pairing transcriptomic and next-generation bisulfite sequencing analyses of the common differentially expressed genes identified 42 CpG island-containing genes in which the methylation level was correlated inversely with gene expression. Forty-eight percent of these genes (20 genes, including Areg, Jag1, and Foxe1) that exhibited the most significant methylation and gene expression changes were involved in key pathways associated with different aspects of liver pathology. Our findings illustrate that gene-specific DNA methylation changes have functional consequences and may be an important component of furan hepatotoxicity and hepatocarcinogenicity.

Transcriptional profiling of the mouse hippocampus supports an NMDAR-mediated neurotoxic mode of action for benzo[a]pyrene

Benzo[a]pyrene (BaP) is a genotoxic carcinogen and a neurotoxicant. The neurotoxicity of BaP is proposed to arise from either genotoxicity leading to neuronal cell death, or perturbed expression of N-methyl-d-aspartate receptor (NMDAR) subunits. To explore these hypotheses, we profiled hippocampal gene expression of adult male Muta(™) Mouse administered 0, 1, 35, or 70 mg BaP/kg bw per day by oral gavage for 3 days. Transcriptional profiles were examined by RNA-sequencing (RNA-seq), DNA microarrays, and real-time quantitative reverse transcription polymerase chain reaction (RT-PCR). BaP-DNA adducts in the cerebellum were quantified by (32) P-post-labeling to measure genotoxicity. RNA-seq revealed altered expression of 0, 260, and 219 genes (P-value < 0.05, fold-change ≥ ± 1.5) following exposure to the low, medium, and high doses, respectively; 54 genes were confirmed by microarrays. Microarray and RT-PCR analysis showed increased expression of NMDAR subunits Grina and Grin2a. In contrast, no effects on DNA-damage response genes were observed despite comparable BaP-DNA adduct levels in the cerebellum and in the lungs and livers of mice at similar BaP doses in previous studies. The results suggest that DNA-damage response does not play a major role in BaP-induced adult neurotoxicity. Meta-analysis revealed that BaP-induced transcriptional profiles are highly correlated with those from the hippocampus of transgenic mice exhibiting similar neurotoxicity outcomes to BaP-exposed mice and rats (i.e., defects in learning and memory). Overall, we suggest that BaP-induced neurotoxicity is more likely to be a consequence of NMDAR perturbation than genotoxicity, and identify other important genes potentially mediating this adverse outcome. Environ. Mol. Mutagen. 57:350-363, 2016. © 2016 Her Majesty the Queen in Right of Canada. Environmental and Molecular Mutagenesis © 2016 Environmental Mutagen Society.

Protective role of garlic oil against oxidative damage induced by furan exposure from weaning through adulthood in adult rat testis

Furan is produced in a wide variety of heat-treated foods via thermal degradation. Furan contamination is found to be relatively high in processed baby foods, cereal products, fruits juices, and canned vegetables. Several studies have demonstrated that furan is a potent hepatotoxin and hepatocarcinogen in rodents. However, few studies have investigated the toxic effects of furan in the testis. In addition, the exact mechanism(s) by which furan exerts toxicity in the testis has not been fully elucidated. In this study, we investigated the potential of furan exposure from weaning through adulthood to induce oxidative stress in adult rat testis, as well as the potential of garlic oil (GO) to ameliorate the induced toxicity. Our results reveal that furan administration significantly reduced serum testosterone levels and increased the levels of malondialdehyde (MDA); furthermore, furan administration decreased significantly the enzymatic activity of testicular antioxidants, including glutathione (GSH), superoxide dismutase (SOD), and catalase (CAT) and induced histopathological alterations in the testis. GO co-administration ameliorated the reduction in testosterone levels and dramatically attenuated the furan-induced oxidative and histopathological changes. In addition, Go significantly down-regulated the increased caspase-3 and cytochrome P450 2E1 (CYP2E1) expression in the furan-treated testis. To the best of our knowledge, this study is the first to demonstrate the furan-induced oxidative changes in the adult rat testis and the protective role of GO to ameliorate these changes through its antioxidant effects and its ability to inhibit CYP2E1 production.

Lack of genotoxic mechanisms in early-stage furan-induced hepatocellular tumorigenesis in gpt delta rats

Furan has been used as an intermediate in the chemical-manufacturing industry and has been shown to contaminate various foods. Although furan induces hepatocellular tumors in rodents, equivocal results from in vitro and in vivo mutagenicity tests have caused controversy regarding the involvement of genotoxic mechanisms in furan-induced carcinogenesis. In the present study, to elucidate the possible mechanisms underlying furan-induced hepatocarcinogenesis, a comprehensive medium-term analysis was conducted using gpt delta rats treated with furan at carcinogenic doses for 13 weeks. In the liver, the frequencies of gpt and Spi^- mutants derived mainly from point and deletion mutations, respectively, were not changed, and there were no furan-specific gpt mutations in furan-treated rats. In contrast, the number and area of glutathione S-transferase placental form (GST-P)- positive foci were significantly increased in the high-dose group. Also, the ratio of PCNA-positive hepatocytes was significantly elevated in the same group, as supported by significant increases in cyclin d1 and cyclin e1 mRNA levels. Thus, it is highly probable that cell proliferation, but not genotoxic mechanisms, contribute to the development of GST-P foci in furan-treated rats. Based on the close relationship between GST-P and neoplastic hepatocytes, these data allowed us to hypothesize that cell proliferation following signal transduction other than the mitogen-activated protein kinase (MAPK)/ERK pathway may play a crucial role in early-stage furan-induced hepatocarcinogenesis. Copyright © 2016 John Wiley & Sons, Ltd.

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

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

Nano-risk Science: application of toxicogenomics in an adverse outcome pathway framework for risk assessment of multi-walled carbon nanotubes

BACKGROUND

A diverse class of engineered nanomaterials (ENMs) exhibiting a wide array of physical-chemical properties that are associated with toxicological effects in experimental animals is in commercial use. However, an integrated framework for human health risk assessment (HHRA) of ENMs has yet to be established. Rodent 2-year cancer bioassays, clinical chemistry, and histopathological endpoints are still considered the 'gold standard' for detecting substance-induced toxicity in animal models. However, the use of data derived from alternative toxicological tools, such as genome-wide expression profiling and in vitro high-throughput assays, are gaining acceptance by the regulatory community for hazard identification and for understanding the underlying mode-of-action. Here, we conducted a case study to evaluate the application of global gene expression data in deriving pathway-based points of departure (PODs) for multi-walled carbon nanotube (MWCNT)-induced lung fibrosis, a non-cancer endpoint of regulatory importance.

METHODS

Gene expression profiles from the lungs of mice exposed to three individual MWCNTs with different physical-chemical properties were used within the framework of an adverse outcome pathway (AOP) for lung fibrosis to identify key biological events linking MWCNT exposure to lung fibrosis. Significantly perturbed pathways were categorized along the key events described in the AOP. Benchmark doses (BMDs) were calculated for each perturbed pathway and were used to derive transcriptional BMDs for each MWCNT.

RESULTS

Similar biological pathways were perturbed by the different MWCNT types across the doses and post-exposure time points studied. The pathway BMD values showed a time-dependent trend, with lower BMDs for pathways perturbed at the earlier post-exposure time points (24 h, 3d). The transcriptional BMDs were compared to the apical BMDs derived by the National Institute for Occupational Safety and Health (NIOSH) using alveolar septal thickness and fibrotic lesions endpoints. We found that regardless of the type of MWCNT, the BMD values for pathways associated with fibrosis were 14.0-30.4 μg/mouse, which are comparable to the BMDs derived by NIOSH for MWCNT-induced lung fibrotic lesions (21.0-27.1 μg/mouse).

CONCLUSIONS

The results demonstrate that transcriptomic data can be used to as an effective mechanism-based method to derive acceptable levels of exposure to nanomaterials in product development when epidemiological data are unavailable.

Hepatic miRNA profiles and thyroid hormone homeostasis in rats exposed to dietary potassium perfluorooctanesulfonate (PFOS)

Perfluorooctanesulfonate (PFOS) has been widely used in a variety of industrial and commercial applications as a surfactant and stain repellent. PFOS causes liver damage (including liver tumors) in experimental animals, primarily via interaction with PPARα and CAR/PXR. We investigated the involvement of microRNAs (miRNAs) in PFOS-induced hepatotoxicity, and mechanisms involved in abnormal thyroid hormone (TH) homeostasis, in the livers of adult male rats exposed in feed to 50mg PFOS/kg diet for 28 days. PFOS-treated rats exhibited expected histopathological and clinical chemistry changes, and global gene expression changes consistent with the involvement of PPARα and CAR/PXR. Thirty-eight miRNAs were significantly altered. Three members of the miR-200 family were the most increased, while miR-122-5p and miR-21-5p were the most decreased, in PFOS-treated rats. Expression of the miR-23b-3p/27b-3p/24-3p cluster also decreased in PFOS-treated animals. Pathway analysis of miRNAs and associated gene expression changes suggests involvement of epithelial to mesenchymal transition (EMT), which is a primary process of tumor cell motility and cancer metastasis. Our analysis also revealed transcripts that may mediate PFOS-induced effects on TH homeostasis including: activation of the CAR/PXR pathway, phase II/III enzymes, and deiodinase. These changes are consistent with low serum TH due to enhanced metabolic clearance of TH. However, most TH hepatic target genes were not altered in a manner consistent with reduced TH signaling, suggesting that PFOS exposure did not induce functional hypothyroidism. Collectively, the study suggests an important role for miRNAs in PFOS-induced hepatotoxicity and provides insight into the effects of PFOS on TH homeostasis.

Transcriptional Profiling of Dibenzo[def,p]chrysene-induced Spleen Atrophy Provides Mechanistic Insights into its Immunotoxicity in MutaMouse

Dibenzo[def,p]chrysene (DBC) is the most carcinogenic polycyclic aromatic hydrocarbon (PAH) examined to date. We investigated the immunotoxicity of DBC, manifested as spleen atrophy, following acute exposure of adult MutaMouse males by oral gavage. Mice were exposed to 0, 2.0, 6.2, or 20.0 mg DBC /kg-bw per day, for 3 days. Genotoxic endpoints (DBC-DNA adducts and lacZ mutant frequency in spleen and bone marrow, and red blood cell micronucleus frequency) and global gene expression changes were measured. All of the genotoxicity measures increased in a dose-dependent manner in spleen and bone marrow. Gene expression analysis showed that DBC activates p53 signaling pathways related to cellular growth and proliferation, which was evident even at the low dose. Strikingly, the expression profiles of DBC exposed mouse spleens were highly inversely correlated with the expression profiles of the only published toxicogenomics dataset of enlarged mouse spleen. This analysis suggested a central role for Bnip3l, a pro-apoptotic protein involved in negative regulation of erythroid maturation. RT-PCR confirmed expression changes in several genes related to apoptosis, iron metabolism, and aryl hydrocarbon receptor signaling that are regulated in the opposite direction during spleen atrophy versus benzo[a]pyrene-mediated splenomegaly. In addition, benchmark dose modeling of toxicogenomics data yielded toxicity estimates that are very close to traditional toxicity endpoints. This work illustrates the power of toxicogenomics to reveal rich mechanistic information for immunotoxic compounds and its ability to provide information that is quantitatively similar to that derived from standard toxicity methods in health risk assessment.

Application of biclustering of gene expression data and gene set enrichment analysis methods to identify potentially disease causing nanomaterials

BACKGROUND

The presence of diverse types of nanomaterials (NMs) in commerce is growing at an exponential pace. As a result, human exposure to these materials in the environment is inevitable, necessitating the need for rapid and reliable toxicity testing methods to accurately assess the potential hazards associated with NMs. In this study, we applied biclustering and gene set enrichment analysis methods to derive essential features of altered lung transcriptome following exposure to NMs that are associated with lung-specific diseases. Several datasets from public microarray repositories describing pulmonary diseases in mouse models following exposure to a variety of substances were examined and functionally related biclusters of genes showing similar expression profiles were identified. The identified biclusters were then used to conduct a gene set enrichment analysis on pulmonary gene expression profiles derived from mice exposed to nano-titanium dioxide (nano-TiO2), carbon black (CB) or carbon nanotubes (CNTs) to determine the disease significance of these data-driven gene sets.

RESULTS

Biclusters representing inflammation (chemokine activity), DNA binding, cell cycle, apoptosis, reactive oxygen species (ROS) and fibrosis processes were identified. All of the NM studies were significant with respect to the bicluster related to chemokine activity (DAVID; FDR p-value = 0.032). The bicluster related to pulmonary fibrosis was enriched in studies where toxicity induced by CNT and CB studies was investigated, suggesting the potential for these materials to induce lung fibrosis. The pro-fibrogenic potential of CNTs is well established. Although CB has not been shown to induce fibrosis, it induces stronger inflammatory, oxidative stress and DNA damage responses than nano-TiO2 particles.

CONCLUSION

The results of the analysis correctly identified all NMs to be inflammogenic and only CB and CNTs as potentially fibrogenic. In addition to identifying several previously defined, functionally relevant gene sets, the present study also identified two novel genes sets: a gene set associated with pulmonary fibrosis and a gene set associated with ROS, underlining the advantage of using a data-driven approach to identify novel, functionally related gene sets. The results can be used in future gene set enrichment analysis studies involving NMs or as features for clustering and classifying NMs of diverse properties.

BMDExpress Data Viewer - a visualization tool to analyze BMDExpress datasets

Regulatory agencies increasingly apply benchmark dose (BMD) modeling to determine points of departure for risk assessment. BMDExpress applies BMD modeling to transcriptomic datasets to identify transcriptional BMDs. However, graphing and analytical capabilities within BMDExpress are limited, and the analysis of output files is challenging. We developed a web‐based application, BMDExpress Data Viewer (http://apps.sciome.com:8082/BMDX_Viewer/), for visualizing and graphing BMDExpress output files. The application consists of “Summary Visualization” and “Dataset Exploratory” tools. Through analysis of transcriptomic datasets of the toxicants furan and 4,4′‐methylenebis(N,N‐dimethyl)benzenamine, we demonstrate that the “Summary Visualization Tools” can be used to examine distributions of gene and pathway BMD values, and to derive a potential point of departure value based on summary statistics. By applying filters on enrichment P‐values and minimum number of significant genes, the “Functional Enrichment Analysis” tool enables the user to select biological processes or pathways that are selectively perturbed by chemical exposure and identify the related BMD. The “Multiple Dataset Comparison” tool enables comparison of gene and pathway BMD values across multiple experiments (e.g., across timepoints or tissues). The “BMDL‐BMD Range Plotter” tool facilitates the observation of BMD trends across biological processes or pathways. Through our case studies, we demonstrate that BMDExpress Data Viewer is a useful tool to visualize, explore and analyze BMDExpress output files. Visualizing the data in this manner enables rapid assessment of data quality, model fit, doses of peak activity, most sensitive pathway perturbations and other metrics that will be useful in applying toxicogenomics in risk assessment. © 2015 Her Majesty the Queen in Right of Canada. Journal of Applied Toxicology published by John Wiley & Sons, Ltd.

We developed BMDExpress Data Viewer, which contains two collections of tools, “Summary Visualization Tools” and “Dataset Exploratory Tools,” to visualize and analyze BMDExpress output files. Through two case studies, we demonstrate the capabilities of graphically examining transcriptomic dose–response datasets in a risk assessment context by comparing and observing trends in transcriptomic benchmark doses (BMDs) for biological processes and pathways. Our results illustrate that BMDExpress Data Viewer is a useful tool to visualize, explore and analyze BMDExpress output files.

Integrating toxicogenomics into human health risk assessment: lessons learned from the benzo[a]pyrene case study

The use of short-term toxicogenomic tests to predict cancer (or other health effects) offers considerable advantages relative to traditional toxicity testing methods. The advantages include increased throughput, increased mechanistic data, and significantly reduced costs. However, precisely how toxicogenomics data can be used to support human health risk assessment (RA) is unclear. In a companion paper ( Moffat et al. 2014 ), we present a case study evaluating the utility of toxicogenomics in the RA of benzo[a]pyrene (BaP), a known human carcinogen. The case study is meant as a proof-of-principle exercise using a well-established mode of action (MOA) that impacts multiple tissues, which should provide a best case example. We found that toxicogenomics provided rich mechanistic data applicable to hazard identification, dose-response analysis, and quantitative RA of BaP. Based on this work, here we share some useful lessons for both research and RA, and outline our perspective on how toxicogenomics can benefit RA in the short- and long-term. Specifically, we focus on (1) obtaining biologically relevant data that are readily suitable for establishing an MOA for toxicants, (2) examining the human relevance of an MOA from animal testing, and (3) proposing appropriate quantitative values for RA. We describe our envisioned strategy on how toxicogenomics can become a tool in RA, especially when anchored to other short-term toxicity tests (apical endpoints) to increase confidence in the proposed MOA, and emphasize the need for additional studies on other MOAs to define the best practices in the application of toxicogenomics in RA.

Comparison of toxicogenomics and traditional approaches to inform mode of action and points of departure in human health risk assessment of benzo[a]pyrene in drinking water

Toxicogenomics is proposed to be a useful tool in human health risk assessment. However, a systematic comparison of traditional risk assessment approaches with those applying toxicogenomics has never been done. We conducted a case study to evaluate the utility of toxicogenomics in the risk assessment of benzo[a]pyrene (BaP), a well-studied carcinogen, for drinking water exposures. Our study was intended to compare methodologies, not to evaluate drinking water safety. We compared traditional (RA1), genomics-informed (RA2) and genomics-only (RA3) approaches. RA2 and RA3 applied toxicogenomics data from human cell cultures and mice exposed to BaP to determine if these data could provide insight into BaP's mode of action (MOA) and derive tissue-specific points of departure (POD). Our global gene expression analysis supported that BaP is genotoxic in mice and allowed the development of a detailed MOA. Toxicogenomics analysis in human lymphoblastoid TK6 cells demonstrated a high degree of consistency in perturbed pathways with animal tissues. Quantitatively, the PODs for traditional and transcriptional approaches were similar (liver 1.2 vs. 1.0 mg/kg-bw/day; lungs 0.8 vs. 3.7 mg/kg-bw/day; forestomach 0.5 vs. 7.4 mg/kg-bw/day). RA3, which applied toxicogenomics in the absence of apical toxicology data, demonstrates that this approach provides useful information in data-poor situations. Overall, our study supports the use of toxicogenomics as a relatively fast and cost-effective tool for hazard identification, preliminary evaluation of potential carcinogens, and carcinogenic potency, in addition to identifying current limitations and practical questions for future work.

Mining the Archives: A Cross-Platform Analysis of Gene Expression Profiles in Archival Formalin-Fixed Paraffin-Embedded Tissues

Formalin-fixed paraffin-embedded (FFPE) tissue samples represent a potentially invaluable resource for transcriptomic research. However, use of FFPE samples in genomic studies has been limited by technical challenges resulting from nucleic acid degradation. Here we evaluated gene expression profiles derived from fresh-frozen (FRO) and FFPE mouse liver tissues preserved in formalin for different amounts of time using 2 DNA microarray protocols and 2 whole-transcriptome sequencing (RNA-seq) library preparation methodologies. The ribo-depletion protocol outperformed the other methods by having the highest correlations of differentially expressed genes (DEGs), and best overlap of pathways, between FRO and FFPE groups. The effect of sample time in formalin (18 h or 3 weeks) on gene expression profiles indicated that test article treatment, not preservation method, was the main driver of gene expression profiles. Meta- and pathway analyses indicated that biological responses were generally consistent for 18 h and 3 week FFPE samples compared with FRO samples. However, clear erosion of signal intensity with time in formalin was evident, and DEG numbers differed by platform and preservation method. Lastly, we investigated the effect of time in paraffin on genomic profiles. Ribo-depletion RNA-seq analysis of 8-, 19-, and 26-year-old control blocks resulted in comparable quality metrics, including expected distributions of mapped reads to exonic, untranslated region, intronic, and ribosomal fractions of the transcriptome. Overall, our results indicate that FFPE samples are appropriate for use in genomic studies in which frozen samples are not available, and that ribo-depletion RNA-seq is the preferred method for this type of analysis in archival and long-aged FFPE samples.

Impact of Genomics Platform and Statistical Filtering on Transcriptional Benchmark Doses (BMD) and Multiple Approaches for Selection of Chemical Point of Departure (PoD)

Many regulatory agencies are exploring ways to integrate toxicogenomic data into their chemical risk assessments. The major challenge lies in determining how to distill the complex data produced by high-content, multi-dose gene expression studies into quantitative information. It has been proposed that benchmark dose (BMD) values derived from toxicogenomics data be used as point of departure (PoD) values in chemical risk assessments. However, there is limited information regarding which genomics platforms are most suitable and how to select appropriate PoD values. In this study, we compared BMD values modeled from RNA sequencing-, microarray-, and qPCR-derived gene expression data from a single study, and explored multiple approaches for selecting a single PoD from these data. The strategies evaluated include several that do not require prior mechanistic knowledge of the compound for selection of the PoD, thus providing approaches for assessing data-poor chemicals. We used RNA extracted from the livers of female mice exposed to non-carcinogenic (0, 2 mg/kg/day, mkd) and carcinogenic (4, 8 mkd) doses of furan for 21 days. We show that transcriptional BMD values were consistent across technologies and highly predictive of the two-year cancer bioassay-based PoD. We also demonstrate that filtering data based on statistically significant changes in gene expression prior to BMD modeling creates more conservative BMD values. Taken together, this case study on mice exposed to furan demonstrates that high-content toxicogenomics studies produce robust data for BMD modelling that are minimally affected by inter-technology variability and highly predictive of cancer-based PoD doses.

Toxicogenomic assessment of liver responses following subchronic exposure to furan in Fischer F344 rats

Furan is a widely used industrial chemical and a contaminant in heated foods. Chronic furan exposure causes cholangiocarcinoma and hepatocellular tumors in rats at doses of 2 mg/kg bw/day or greater, with gender differences in frequency and severity. The hepatic transcriptional alterations induced by low doses of furan (doses below those previously tested for induction of liver tumors) and the potential mechanisms underlying gender differences are largely unexplored. We used DNA microarrays to examine the global hepatic mRNA and microRNA transcriptional profiles of male and female rats exposed to 0, 0.03, 0.12, 0.5 or 2 mg/kg bw/day furan over 90 days. Marked gender differences in gene expression responses to furan were observed, with many more altered genes in exposed males than females, confirming the increased sensitivity of males even at the low doses. Pathway analysis supported that key events in furan-induced liver tumors in males include gene expression changes related to oxidative stress, apoptosis and inflammatory response, while pathway changes in females were consistent with primarily adaptive responses. Pathway benchmark doses (BMDs) were estimated and compared to relevant apical endpoints. Transcriptional pathway BMDs could only be examined in males. These median BMDs ranged from 0.08 to 1.43 mg/kg bw/day and approximated those derived from traditional histopathology. MiR-34a (a P53 target) was the only microRNA significantly increased at the 2 mg/kg bw/day, providing evidence to support the importance of apoptosis and cell proliferation in furan hepatotoxicity. Overall, this study demonstrates the use of transcriptional profiling to discern mode of action and mechanisms involved in gender differences.