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

Rewiring cis-2-butene-1,4-dial mediated urinary metabolomics fingerprints of short-term exposure to furan

Furan represents one of the dietary-sourced persistent organic pollutants and thermal processing contaminants. Given its widespread occurrence in food and various toxicological effects, accurately assessing furan exposure is essential for informing public health risks. Furan is metabolized to a reactive primary product, cis-2-butene-1,4-dial (BDA) upon absorption. Some of the resulting BDA-derived metabolites have been proposed as potential exposure biomarkers of furan. However, the lack of quantification for recognized and feasible furan biomarkers has hampered the development of internal exposure risk assessment of furan. In this study, we employed reliable non-targeted metabolomics techniques to uncover urinary furan metabolites and elucidate their chemical structures. We characterized 8 reported and 11 new furan metabolites derived from the binding of BDA with glutathione (GSH), biogenic amines, and/or amino acids in the urine of male rats subjected to varying doses of furan. Notably, a mono-GSH-BDA adduct named cyclic GSH-BDA emerged as a highly prospective specific biomarker of furan exposure, as determined by an ultrahigh-performance liquid chromatography-tandem mass spectrometry method. Cyclic GSH-BDA demonstrated a robust mass spectrometry ion response intensity and exhibited evident time- and dose response. Additionally, we conducted a comprehensive profiling of the kinetics of potential furan biomarkers over time to capture the metabolic dynamics of furan in vivo. Most urinary furan metabolites reached peak concentrations at either the first (3 h) or second (6 h) sampling time point and were largely eliminated within 36 h following furan treatment. The present study provides novel insights into furan metabolism and sheds light on the biomonitoring of furan exposure.

The role of microRNAs in acrylamide toxicity

The chemical compound known as Acrylamide (AA) is employed in different industries worldwide and is also found in thermal-processed food. AA has been acting as a reproductive toxicant, carcinogen, and neurotoxic in various animals, which may promote several toxic impacts in animal and human species. Up to now, various studies have focused on the harmful mechanisms and intervention actions of AA. However, the underlying mechanisms that AA and its toxic effects can exert have remained uncertain. MicroRNAs (miRNAs) are a class of short, non-coding RNAs that are able to act as epigenetic regulators. These molecules can regulate a wide range of cellular and molecular processes. In this regard, it has been shown that different chemical agents can dysregulate miRNAs. To determine the possible AA targets along with mechanisms of its toxicity, it is helpful to study the alteration in the profiles of miRNA regulation following AA intake. The current research aimed to evaluate the miRNAs' mediatory roles upon the AA's toxic potentials. This review study discussed the AA, which is made within the food matrix, the way it is consumed, and the potential impacts of AA on miRNAs and its association with different cancer types and degenerative diseases. The findings of this review paper indicated that AA might be capable of altering miRNA signatures in different tissues and exerting its carcinogen effects.

Comparative genomic analysis of Fischer F344 rat livers exposed for 90 days to 3-methylfuran or its parental compound furan

Furan is a naturally forming compound found in heat-processed foods such as coffee, canned meats, and jarred baby food. It is concurrently found with analogues including 2-methylfuran (2-MF) and 3-methylfuran (3-MF), and toxicity studies demonstrate all are potent liver toxins. Toxicity studies found 3-MF is more toxic than either furan, or 2-MF. The present analysis assesses the transcriptional response in liver samples taken from male Fischer (F344) rats exposed to furan or 3-MF from 0 to 2.0 and 0-1.0 mg/kg bw/day, respectively, for 90 days. Transcriptional analyses found decreased liver function and fatty acid metabolism are common responses to both furan and 3-MF exposure. Furan liver injury promotes a ductular reaction through Hippo and TGFB signalling, which combined with increased immune response results in ameliorating perturbed bile acid homeostasis in treated rats. Failure to activate these pathways in 3-MF exposed rats and decreased p53 activity leads to cholestasis, and increased toxicity. Finally, BMD analysis indicate many of the most sensitive pathways affected by furan and 3-MF exposure relate to metabolism - malate dehydrogenase and glucose metabolism with BMDLs of 0.03 and 0.01 mg/kg bw/day for furan and 3-MF exposure, respectively, which agrees with BMDLs previously reported for apical and microarray data.

High-Throughput Transcriptomics Differentiates Toxic versus Non-Toxic Chemical Exposures Using a Rat Liver Model

To address the challenge of limited throughput with traditional toxicity testing, a newly developed high-throughput transcriptomics (HTT) platform, together with a 5-day in vivo rat model, offers an alternative approach to estimate chemical exposures and provide reasonable estimates of toxicological endpoints. This study contains an HTT analysis of 18 environmental chemicals with known liver toxicity. They were evaluated using male Sprague Dawley rats exposed to various concentrations daily for five consecutive days via oral gavage, with data collected on the sixth day. Here, we further explored the 5-day rat model to identify potential gene signatures that can differentiate between toxic and non-toxic liver responses and provide us with a potential histopathological endpoint of chemical exposure. We identified a distinct gene expression pattern that differentiated non-hepatotoxic compounds from hepatotoxic compounds in a dose-dependent manner, and an analysis of the significantly altered common genes indicated that toxic chemicals predominantly upregulated most of the genes and several pathways in amino acid and lipid metabolism. Finally, our liver injury module analysis revealed that several liver-toxic compounds showed similarities in the key injury phenotypes of cellular inflammation and proliferation, indicating potential molecular initiating processes that may lead to a specific end-stage liver disease.

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.

Modulation of Disordered Bile Acid Homeostasis and Hepatic Tight Junctions Using Salidroside against Hepatocyte Apoptosis in Furan-Induced Mice

Furan, a processing-induced food contaminant, has attracted great attention due to its hepatotoxicity. To further investigate the underlying mechanism of salidroside (SAL) alleviating furan-induced liver damage, we divided Balb/c mice into the control group, the furan (8 mg/kg/day) group, and three groups of three different doses of SAL (10/20/40 mg/kg/day) in the current research. The shifted serum profile was observed through untargeted metabonomics, to which the bile acid metabolism was related, and the alleviating effect of SAL against furan-induced apoptosis was caused by the metabolism. Target bile acid quantification for the liver and serum showed that SAL positively regulated the homeostasis of bile acids disturbed by furan. Meanwhile, SAL significantly upregulated the synthesis genes of bile acids (Cyp7a1, Cyp7b1, Cyp8b1, and Cyp27a1) and the uptake transport genes (Ntcp and Oatps) and downregulated the efflux transport genes (Bsep, Ost-α, Ost-β, Mrp2, and Mrp4). Transmission electron microscopy of the bile canaliculi and tight junctions and the levels of tight junction marker proteins (ZO-1, occludin, and claudin-1) confirmed that the disruption of the hepatic tight junction was inhibited by SAL. The connection between the apoptosis- and tight junction-related proteins was observed through the construction of a protein-protein interaction network. SAL suppressed the furan-induced hepatocyte apoptosis evidenced by the detection of TUNEL and Bax, Bcl-2, and caspase-3 levels. Taken together, SAL alleviated furan-induced hepatocyte apoptosis via altering the disordered homeostasis of bile acids and hepatic tight junctions.

MiR-27a-5p regulates acrylamide-induced mitochondrial dysfunction and intrinsic apoptosis via targeting Btf3 in rats

Acrylamide (AA), a potential carcinogen, is commonly formed in foods rich in carbohydrates at high heat. It is known that AA-induced mitochondrial dysfunction is responsible for its toxicity. Previously we found AA exposure increased miR-27a-5p expression in livers of SD rats. Here, the regulation mechanism of miR-27a-5p in mitochondrial dysfunction was investigated in rat liver cell lines (IAR20) and SD rats. The results showed that the overexpressed miR-27a-5p contributes to modulating mitochondrial dysfunction and Btf3 is identified as its target gene. The knockdown of Btf3 increases the cleaved PARP1 level and the phosphorylation of ATM and p53, which results in mitochondria-dependent apoptosis. Therefore, the miR-27a-5p-Btf3-ATM-p53 axis might play a vital role in the promotion of AA-induced cell apoptosis through disrupting mitochondrial structure and function. This would provide a potential target for the assessment and intervention of AA toxicity.

Proteomic Analysis of Subchronic Furan Exposure in the Liver of Male Fischer F344 Rats

Furan is a volatile compound formed during the thermal processing of foods. Chronic exposure has been shown to cause cholangiocarcinoma and hepatocellular tumors in rodent models. We conducted a 90 day subchronic study in Fisher 344 rats exposed to various doses by gavage to determine the NOAEL. Previous reports have outlined changes in the liver using gross necropsy examination, histopathology, clinical biochemistry, hematology, immunohistochemistry, and toxicogenomics. The data revealed that males were more sensitive than females. The focus of this study was to evaluate the toxicoproteomic changes by 2-dimensional differential in gel electrophoresis followed by mass spectrometry analysis. To compliment previous studies, protein expression changes were evaluated of male animals after 90 days of exposure to doses of 0, 0.03, 0.5, and 8.0 mg/kg bw/d. Significant statistical treatment-related changes compared to the controls identified 45 protein spots containing 38 unique proteins. Proteins identified are implicated in metabolism, redox regulation, protein folding/proteolysis as well as structural and transport proteins. At lower doses, multiple cytoprotective pathways are activated to maintain a homeostasis but ultimately the loss of protein function and impairment of several pathways could lead to adverse health effects at higher doses of furan administration.

Protocatechuic acid protects against hepatorenal toxicities in rats exposed to Furan

Furan formed in processed food is hepatotoxic and likely carcinogenic in humans. We investigated protocatechuic acid (PCA) protective role in rats' hepatorenal function treated with furan. Rats were grouped and treated as follows: Control, PCA (50 mg/kg), furan alone (8 mg/kg), furan + PCA1 (25 + 8 mg/kg), and furan + PCA2 (50 + 8 mg/kg). Upon sacrifice, evaluation of hepatorenal function, oxidative stress status, reactive oxygen and nitrogen species (RONS), lipid peroxidation (LPO), myeloperoxidase (MPO) activity, among nitric oxide (NO) levels were performed. Cytokine levels (IL-10, IL-1ß, TNF-alpha), Caspase 3 and 9 activities, and histopathological examination were also assessed. We found that the final body and relative liver weights changed significantly (p < 0.05) in treated groups. Hepatic transaminases, urea, and creatinine increased (p < 0.05) in furan only treated group, and reduced in PCA co-treated groups. The furan-induced decrease in antioxidant status increased RONS, and LPO levels were alleviated (p < 0.05) by PCA co-treatment. Furthermore, furan-mediated increase in NO, IL-1ß, TNF-alpha levels, MPO, Cas-3, and 9 activities and suppressed IL-10 levels was reversed accordingly in rats' kidney and liver co-treated with PCA. The extent of furan-mediated hepatorenal lesions was lessened in PCA co-treated rats. Our findings suggest that PCA protects against oxido-inflammatory pathways, enhanced caspases 3 and 9 activations induced by furan in rat hepatorenal system.

A set of six Gene expression biomarkers and their thresholds identify rat liver tumorigens in short-term assays

Traditional methods for cancer risk assessment are retrospective, resource-intensive, and not feasible for the vast majority of environmental chemicals. In earlier studies, we used a set of six biomarkers to accurately identify liver tumorigens in transcript profiles derived from chemically-treated rats using either a Toxicological Priority Index (ToxPi) approach or using derived biomarker thresholds for cancer. The biomarkers consisting of 7-113 genes are used to predict the most common liver cancer molecular initiating events: genotoxicity, cytotoxicity and activation of the xenobiotic receptors AhR, CAR, ER, and PPARα. In the present study, we apply and evaluate the performance of these methods for cancer prediction in an independent rat liver study of 44 chemicals (6 h-7d exposures) examined by Affymetrix arrays. In the first approach, ToxPi ranking of biomarker scores consistently gave the highest scores to tumorigenic chemical-dose pairs; balanced accuracies for identification of liver tumorigenic chemicals were up to 89 %. The second approach used tumorigenic thresholds derived in the present study or from our earlier study that were set at the maximum value for chemical-dose exposures without detectable liver tumor outcomes. Using these thresholds, balanced accuracies were up to 90 %. Both approaches identified all tumorigenic chemicals. Almost all of the tumorigenic chemicals activated more than one MIE. We also compared biomarker responses between two types of profiling platforms (Affymetrix full-genome array, TempO-Seq 1500+ array containing ∼2600 genes) and found that the lack of the full set of biomarker genes on the 1500+ array resulted in decreased ability to identify chemicals that activate the MIEs. Overall, these results demonstrate that predictive approaches based on the 6 biomarkers could be used in short-term assays to identify chemicals and their doses that induce liver tumors, the most common endpoint in rodent bioassays.

Environmental epitranscriptomics

Chemical modifications of RNA molecules have gained increasing attention since evidence emerged for their substantive roles in a range of biological processes, such as the stability and translation of mRNA transcripts. More than 150 modifications have been identified in different organisms to date, collectively known as the 'epitranscriptome', with 6-methyladenosine (m⁶A), 5-methylcytidine (m⁵C), pseudouridine and N1-methyladenosine (m¹A) the most extensively investigated. Although we are just beginning to elucidate the roles of these modifications in cellular functions, there is already evidence for their dysregulation in diseases such as cancer and neurodevelopmental disorders. There is currently more limited knowledge regarding how environmental exposures affect the epitranscriptome and how this may mediate disease risk, but evidence is beginning to emerge. Here, we review the current evidence for the impact of environmental exposures such as benzo[a]pyrene, bisphenol A, pesticides, metals and nanoparticles upon RNA modifications and the expression of their 'writers' (methyl transferases), 'erasers' (demethylases) and 'readers'. We discuss future directions of the field and identify areas of particular promise and consider the technical challenges that are faced.

The crucial role of metabolic regulation in differential hepatotoxicity induced by furanoids in Dioscorea bulbifera

Diterpenoid lactones (DLs), a group of furan-containing compounds found in Dioscorea bulbifera L. (DB), have been reported to be associated with hepatotoxicity. Different hepatotoxicities of these DLs have been observed in vitro, but reasonable explanations for the differential hepatotoxicity have not been provided. Herein, the present study aimed to confirm the potential factors that contribute to varied hepatotoxicity of four representative DLs (diosbulbins A, B, C, F). In vitro toxic effects were evaluated in various cell models and the interactions between DLs and CYP3A4 at the atomic level were simulated by molecular docking. Results showed that DLs exhibited varied cytotoxicities, and that CYP3A4 played a modulatory role in this process. Moreover, structural variation may cause different affinities between DLs and CYP3A4, which was positively correlated with the observation of cytotoxicity. In addition, analysis of the glutathione (GSH) conjugates indicated that reactive intermediates were formed by metabolic oxidation that occurred on the furan moiety of DLs, whereas, GSH consumption analysis reflected the consistency between the reactive metabolites and the hepatotoxicity. Collectively, our findings illustrated that the metabolic regulation played a crucial role in generating the varied hepatotoxicity of DLs.

Markers involved in proinflammatory effects by environmental toxicants

Toxicological studies have identified polycyclic aromatic hydrocarbons (PAH) in human breast milk, smoked and barbequed food, although the largest contribution of PAH intake into the body are cereals and cereals products. The major effects attributable to PAH appeared to occur in the liver, lungs, the hematopoietic system, and the kidney. Nevertheless, more precise mechanisms by which PAH initiates its pathological features are not fully understood. In the present study, we evaluated levels of myeloperoxidase activity, its association with nitric oxide synthesis (NO), levels of uric acid (UA) in circulating blood and glucose in female rats exposed to environmental toxicants. A higher concentration of hydrogen peroxide activates myeloperoxidase, which acts as a leucocyte attractant, contributing to enhanced iNOS activity. In parallel, uric acid in addition to its pro-inflammatory effects aggravates insulin resistance and hyperglycemia, which worsens the process. Our findings suggest potential intermediate mechanisms involved in the inflammatory effects of PAH, which might give insight for the involvement of environmental toxicants not only in carcinogenesis but also in its association with acute cardiovascular disease and induction of multi-organ damage. The development of iNOS inhibitors might be beneficial in certain inflammatory disorders.

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.

Early microRNA indicators of PPARα pathway activation in the liver

MicroRNAs (miRNAs) are short non-coding RNA species that play key roles in post-transcriptional regulation of gene expression. MiRNAs also serve as a promising source of early biomarkers for different environmental exposures and health effects, although there is limited information linking miRNA changes to specific target pathways. In this study, we measured liver miRNAs in male B6C3F1 mice exposed to a known chemical activator of the peroxisome proliferator-activated receptor alpha (PPARα) pathway, di(2-ethylhexyl) phthalate (DEHP), for 7 and 28 days at concentrations of 0, 750, 1500, 3000, or 6000 ppm in feed. At the highest dose tested, DEHP altered 61 miRNAs after 7 days and 171 miRNAs after 28 days of exposure, with 48 overlapping miRNAs between timepoints. Analysis of these 48 common miRNAs indicated enrichment in PPARα–related targets and other pathways related to liver injury and cancer. Four of the 10 miRNAs exhibiting a clear dose trend were linked to the PPARα pathway: mmu-miRs-125a-5p, -182−5p, -20a−5p, and -378a−3p. mmu-miRs-182−5p and -378a−3p were subsequently measured using digital drop PCR across a dose range for DEHP and two related phthalates with weaker PPARα activity, di-n-octyl phthalate and n-butyl benzyl phthalate, following 7-day exposures. Analysis of mmu-miRs-182−5p and -378a−3p by transcriptional benchmark dose analysis correctly identified DEHP as having the greatest potency. However, benchmark dose estimates for DEHP based on these miRNAs (average 163; range 126−202 mg/kg-day) were higher on average than values for PPARα target genes (average 74; range 29−183 mg/kg-day). These findings identify putative miRNA biomarkers of PPARα pathway activity and suggest that early miRNA changes may be used to stratify chemical potency.

Mode of action-based risk assessment of genotoxic carcinogens

The risk assessment of chemical carcinogens is one major task in toxicology. Even though exposure has been mitigated effectively during the last decades, low levels of carcinogenic substances in food and at the workplace are still present and often not completely avoidable. The distinction between genotoxic and non-genotoxic carcinogens has traditionally been regarded as particularly relevant for risk assessment, with the assumption of the existence of no-effect concentrations (threshold levels) in case of the latter group. In contrast, genotoxic carcinogens, their metabolic precursors and DNA reactive metabolites are considered to represent risk factors at all concentrations since even one or a few DNA lesions may in principle result in mutations and, thus, increase tumour risk. Within the current document, an updated risk evaluation for genotoxic carcinogens is proposed, based on mechanistic knowledge regarding the substance (group) under investigation, and taking into account recent improvements in analytical techniques used to quantify DNA lesions and mutations as well as "omics" approaches. Furthermore, wherever possible and appropriate, special attention is given to the integration of background levels of the same or comparable DNA lesions. Within part A, fundamental considerations highlight the terms hazard and risk with respect to DNA reactivity of genotoxic agents, as compared to non-genotoxic agents. Also, current methodologies used in genetic toxicology as well as in dosimetry of exposure are described. Special focus is given on the elucidation of modes of action (MOA) and on the relation between DNA damage and cancer risk. Part B addresses specific examples of genotoxic carcinogens, including those humans are exposed to exogenously and endogenously, such as formaldehyde, acetaldehyde and the corresponding alcohols as well as some alkylating agents, ethylene oxide, and acrylamide, but also examples resulting from exogenous sources like aflatoxin B1, allylalkoxybenzenes, 2-amino-3,8-dimethylimidazo[4,5-f] quinoxaline (MeIQx), benzo[a]pyrene and pyrrolizidine alkaloids. Additionally, special attention is given to some carcinogenic metal compounds, which are considered indirect genotoxins, by accelerating mutagenicity via interactions with the cellular response to DNA damage even at low exposure conditions. Part C finally encompasses conclusions and perspectives, suggesting a refined strategy for the assessment of the carcinogenic risk associated with an exposure to genotoxic compounds and addressing research needs.

Protective effect of Spirulina platensis against physiological, ultrastructural and cell proliferation damage induced by furan in kidney and liver of rat

The modulatory role of the Spirulina platensis (SP) against furan-induced (FU) hepatic and renal damage was assessed in this study. For achieving this, sixty rats were distributed into six groups: control group, SP-administered group (300 mg/kg b.wt orally for 28 days), a FU-intoxicated group (16 mg/kg b.wt, orally, daily for 28 days), protective co-treated group SP/F (administered SP 300 mg/kg b.wt, one week before, and concurrently with FU intoxication), therapeutic co-treated group FU/SP (administered SP 300 mg/kg b.wt, one week after FU intoxication for 28 days) and protective/therapeutic co-treated group SP/FU/SP (administered SP one week before and after, concurrently with FU intoxication). Subsequently, the biochemical responses and the histology of hepatic and renal tissues in treated rats were assessed. The results indicated that FU intoxication induced a significant hepato- and nephropathy represented by the elevation in the values of tissue injury biomarkers and reduction in protein levels. Histologically, a wide range of morphological, cytotoxic, inflammatory, and vascular alterations as well as downregulation in the immunoexpression of the proliferating cell nuclear antigen (PCNA) and the proliferation-associated nuclear antigen (Ki-67) were induced by FU intoxication. Oral SP administration, particularly in the protective/therapeutic co-treated group, markedly supressed the serum levels of the tissue injury biomarkers, diminished the inflammatory response, restored the cytotoxic alterations, upregulated the immunoexpression of PCNA and Ki-67, and restored the perturbed morphology of the hepatic and renal tissues. In conclusion, the obtained data demonstrated that SP co-administration elicits both protective and therapeutic potential against the FU-induced hepato- and nephropathy.

Concentration-dependent transcriptome of zebrafish embryo for environmental chemical assessment

Mechanistic information is essential to screen and predict the adverse effects of a large number of chemicals during early-life exposure. Concentration-dependent omics can capture the extent of perturbations of biological pathways or processes and provide information on the mechanism of toxicity. However, the application of concentration-dependent transcriptome to assess the developmental toxicity of environmental chemicals is still limited. Here, twelve chemicals representing five different modes of action (MOAs) were tested by the concentration-dependent reduced zebrafish transcriptome approach (CRZT) in combination with a phenotype-based high content screen (PHCS). The responsiveness, sensitivity and mechanistic differentiation of CRZT were validated in comparison with PHCS. First, PHCS identified 10 chemicals with obvious embryotoxicity (LD₅₀ range: 2.11-70.68 μM), while the potencies of the biological pathways perturbed by 12 chemicals (POD_path20 range: 0.002-2.1 μM) were demonstrated by CRZT. Second, although the potency of the transcriptome perturbations was positively correlated with lethality (LD₅₀) (R² = 0.64, P-value < 0.05) for most tested chemicals, BbF was non-embryotoxic but was the most potent on the perturbance of biological pathways. Finally, the profiles of the perturbed biological processes and the transcriptome potency (POD_path20) captured by CRZT could effectively classify most chemicals corresponding to their known MOAs. In summary, CRZT could significantly improve testing the developmental toxicity of environmental chemicals.

A review on furan: Formation, analysis, occurrence, carcinogenicity, genotoxicity and reduction methods

Furan (C₄H₄O) is a volatile, heterocyclic and carcinogenic heterocyclic chemical compound occurring in a wide range of thermally processed foods. Several studies have been conducted to analyze the formation conditions, triggering furan formation via model systems. Furan can be encountered via various pathways including thermal degradation, oxidation of polyunsaturated fatty acids, thermal rearrangement of carbohydrates in the presence of amino acids, thermal degradation of certain amino acids. Furan has been proven to cause cancer in experimental animal models and classified as a possible human carcinogen by International agency for research on cancer based on sufficient evidences. Thus, different strategies should be developed to reduce furan contents in commercially available food stuffs while food processing. This review summarizes some current evidences of furan formation from different precursors, analytical methods for its detection, and its toxicity that might lead to carcinogenicity and genotoxicity with human risk assessment. In addition, furan occurrence in different thermally processed foods entailed by several recent studies as well as furan mitigation strategies during food processing have also been illustrated in this review.

Protective Role of Propolis on Low and High Dose Furan-induced Hepatotoxicity and Oxidative Stress in Rats

Introduction

The aim of this study was to evaluate potential protective effects of propolis on furan-induced hepatic damage by assessing the levels of malondialdehyde (MDA) and reduced glutathione (GSH), antioxidant enzyme activities, and histopathological changes in the liver.

Material and Methods

Albino Wistar rats were divided into six groups: a control, propolis-treated (100 mg/kg b.w./day), low-dose furan-treated (furan-L group; 2 mg/kg b.w./day), high-dose furan-treated (furan-H group; 16 mg/kg b.w./day), furan-L+propolis treated, and furan-H+propolis treated group. Propolis and furan were applied by gavage; propolis for 8 days, and furan for 20 days in furan-L groups and 10 days in furan-H groups.

Results

While MDA levels were elevated in furan-treated groups, levels of GSH and activities of antioxidant enzymes decreased (p < 0.001). The levels of MDA and GSH and activities of antioxidant enzymes were normal in the furan+propolis groups, especially in the furan-L+propolis group (p < 0.001). While the aspartate transaminase, alanine transaminase, alkaline phosphatase, and lactate pdehydrogenase activities were elevated in the furan-H treated group (p < 0.05 and p < 0.001), they were unchanged in the furan-L treated group. Histopathologically, several lesions were observed in the liver tissues of the furan-treated groups, especially in the higher-dose group. It was determined that these changes were milder in both of the furan+propolis groups.

Conclusion

The results indicate that propolis exhibits good hepatoprotective and antioxidant potential against furan-induced hepatocellular damage in rats.

Comparison of Furans Formation and Volatile Aldehydes Profiles of Four Different Vegetable Oils During Thermal Oxidation

The evolution of volatile aldehydes and the conversion of oxygenated ityβ-unsaturated aldehydes (OαβUAs) into furans were compared in four vegetable oils (soybean oil, olive oil [OVO], peanut oil [PO], and perilla oil [PAO]) thermally oxidized at temperatures of 150, 180, and 210 °C for 10 hr/day over a 3-day period. Results showed that 2 alkyl furans and 23 volatile aldehydes including 4 toxic OdβUAs were detected by GC-MS. The original fatty acid compositions of the oils played a key role in the type and concentration of those volatile compounds. 4-Hydroxy-2-hexenal (HHE) and ethyl furan were only detected in PAO with a high content of linolenic acid, while the greatest level of pentyl furan was detected in PO with abundant linoleic acid. Greater amounts of 4-hydroxy-(E)-2-nonenal (HNE) and 4-oxo-(E)-2-nonenal (ONE) were formed in the OVO with abundant oleic acid. The close relativity of HHE and ethyl furan was also demonstrated. With principal component analysis, these vegetable oils could be discriminated based on their fatty acids and volatile compounds. The loading plot confirmed that HHE and ethyl furan were derived from the linolenic acid oxidation and degradation. PRACTICAL APPLICATION: The chemometric results showed that the formation of the volatile components during heating in different vegetable oils has close correlation with the original fatty acids composition of vegetable oils. Our research has also confirmed the presence of toxic OɑβUAs in oils after heating. Considering that they are proven to generate lots of degenerative diseases, further studies are needed to establish the risk level of using certain oils in frying and seek effective methods to inhibit their formation.

Toxicogenomics - What added Value Do These Approaches Provide for Carcinogen Risk Assessment?

It is still a major challenge to protect humans at workplaces and in the environment. To cope with this task, it is a prerequisite to obtain detailed information on the extent of chemical perturbations of biological pathways, in particular, adaptive vs. adverse effects and the dose-response relationships. This knowledge serves as the basis for the classification of non-carcinogens and carcinogens and for further distinguishing carcinogens in genotoxic (DNA damaging) or non-genotoxic compounds. Basing on quantitative dose-response relationships, points of departures can be derived for chemical risk assessment. In recent years, new methods have shown their capability to support the established rodent models of carcinogenicity testing. In vitro high throughput screening assays assess more comprehensively cell response. In addition, omics technologies were applied to study the mode of action of chemicals whereby the term "toxicogenomics" comprises various technologies such as transcriptomics, epigenomics, or metabolomics. This review aims to summarize the current state of toxicogenomic approaches in risk science and to compare them with established ones. For example, measurement of global transcriptional changes generates meaningful information for toxicological risk assessment such as accurate classification of genotoxic/non-genotoxic carcinogens. Alteration in mRNA expression offers previously unknown insights in the mode of action and enables the definition of key events. Based on these, benchmark doses can be calculated for the transition from an adaptive to an adverse state. In short, this review assesses the potential and challenges of transcriptomics and addresses the impact of other omics technologies on risk assessment in terms of hazard identification and dose-response assessment.

PBDE-47 and PBDE mixture (DE-71) toxicities and liver transcriptomic changes at PND 22 after in utero/postnatal exposure in the rat

Pentabromodiphenyl ethers (PBDE) are found in human tissue, in household dust, and in the environment, and a particular concern is the potential for the induction of cancer pathways from these fat-soluble persistent organic pollutants. Only one PBDE cancer study has been conducted and that was for a PBDE mixture (DE-71). Because it is not feasible to test all PBDE congeners in the environment for cancer potential, it is important to develop a set of biological endpoints that can be used in short-term toxicity studies to predict disease outcome after long-term exposures. In this study, PBDE-47 was selected as the test PBDE congener to evaluate and compare toxicity to that of the carcinogenic PBDE mixture. The toxicities of PBDE-47 and the PBDE mixture were evaluated at PND 22 in Wistar Han rat (Crl: WI (Han)) pups after in utero/postnatal exposure (0, 0.1, 15, or 50 mg/kg; dams, GD6-21; pups, PND 12-PND 21; oral gavage daily dosing). By PND 22, PBDE-47 caused centrilobular hypertrophy and fatty change in liver, and reduced serum thyroxin (T4) levels; similar effects were also observed after PBDE mixture exposure. Transcriptomic changes in the liver included induction of cytochrome p450 transcripts and up-regulation of Nrf2 antioxidant pathway transcripts and ABC membrane transport transcripts. Decreases in other transport transcripts (ABCG5 & 8) provided a plausible mechanism for lipid accumulation, characterized by a treatment-related liver fatty change after PBDE-47 and PBDE mixture exposure. The benchmark dose calculation based on liver transcriptomic data was generally lower for PBDE-47 than for the PBDE mixture. The up-regulation of the Nrf2 antioxidant pathway and changes in metabolic transcripts after PBDE-47 and PBDE mixture exposure suggest that PBDE-47, like the PBDE mixture (NTP 2016, TR 589), could be a liver toxin/carcinogen after long-term exposure.

Differential susceptibility of kidneys and livers to proliferative processes and transcriptional level of the genes encoding desmin, vimentin, connexin 43, and nestin in rats exposed to furan

In this study, we aimed to assess the differential toxic impact, induced by furan exposure, on the liver and kidney tissues by estimating reactive oxygen species (ROS) level, total antioxidant capacity (TAC), oxidative damage, and the tissue injury markers in a male rat model. To explain such impacts, 20 rats were assigned into two groups: a control group, where rats were administered corn oil as a vehicle, and a furan-administered group, where furan was orally administered to rats at a dose of 16 mg/kg b wt/day (five days per week over eight weeks). The transcriptional levels of intermediate filament proteins (desmin, vimentin, nestin, and connexin 43) were assessed by using quantitative real-time polymerase chain reaction (PCR), and the cell proliferation markers (proliferating cell nuclear antigen [PCNA] and proliferation-associated nuclear antigen [Ki-67]) were recognized by immunohistochemical analysis. Furthermore, the ultrastructural changes of liver and kidney were monitored using electron microscopy. Our findings showed that furan exposure could induce hepatic and renal damage to different extents. Furan can increase the ROS content, oxidative damage indices, and liver tissue injury indices but not kidney injury indices. Furthermore, it decreases the TAC in the serum of exposed rats. In addition, furan exposure was associated with changes in the mRNA expression pattern of intermediate filament proteins in both kidney and liver tissues. Moreover, furan enhances the expression of PCNA and Ki-67 in the liver tissues but not in the kidney tissues. The ultrastructure evaluation revealed the incidence of glomerular podocyte degeneration and hepatocyte injury. These results conclusively demonstrate that the deleterious effects of furan are caused by promoting fibrosis and hepatocyte proliferation in liver tissues and triggering podocyte injury in the kidney tissues.

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

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

BIOLOGICAL SIGNIFICANCE

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

Toxicoproteomic assessment of liver responses to acute pyrrolizidine alkaloid intoxication in rats

A toxicoproteomic study was performed on liver of rats treated with retrorsine (RTS), a representative hepatotoxic pyrrolizidine alkaloid at a toxic dose (140 mg/kg) known to cause severe acute hepatotoxicity. By comparing current data with our previous findings in mild liver lesions of rats treated with a lower dose of RTS, seven proteins and three toxicity pathways of vascular endothelial cell death, which was further verified by observed sinusoidal endothelial cell losses, were found uniquely associated with retrorsine-induced hepatotoxicity. This toxicoproteomic study of acute pyrrolizidine alkaloid intoxication lays a foundation for future investigation to delineate molecular mechanisms of pyrrolizidine alkaloid-induced hepatotoxicity.

Exposure assessment of process-related contaminants in food by biomarker monitoring

Exposure assessment is a fundamental part of the risk assessment paradigm, but can often present a number of challenges and uncertainties. This is especially the case for process contaminants formed during the processing, e.g. heating of food, since they are in part highly reactive and/or volatile, thus making exposure assessment by analysing contents in food unreliable. New approaches are therefore required to accurately assess consumer exposure and thus better inform the risk assessment. Such novel approaches may include the use of biomarkers, physiologically based kinetic (PBK) modelling-facilitated reverse dosimetry, and/or duplicate diet studies. This review focuses on the state of the art with respect to the use of biomarkers of exposure for the process contaminants acrylamide, 3-MCPD esters, glycidyl esters, furan and acrolein. From the overview presented, it becomes clear that the field of assessing human exposure to process-related contaminants in food by biomarker monitoring is promising and strongly developing. The current state of the art as well as the existing data gaps and challenges for the future were defined. They include (1) using PBK modelling and duplicate diet studies to establish, preferably in humans, correlations between external exposure and biomarkers; (2) elucidation of the possible endogenous formation of the process-related contaminants and the resulting biomarker levels; (3) the influence of inter-individual variations and how to include that in the biomarker-based exposure predictions; (4) the correction for confounding factors; (5) the value of the different biomarkers in relation to exposure scenario's and risk assessment, and (6) the possibilities of novel methodologies. In spite of these challenges it can be concluded that biomarker-based exposure assessment provides a unique opportunity to more accurately assess consumer exposure to process-related contaminants in food and thus to better inform risk assessment.

Mechanism-based risk assessment strategy for drug-induced cholestasis using the transcriptional benchmark dose derived by toxicogenomics

Cholestasis is one of the major causes of drug-induced liver injury (DILI), which can result in withdrawal of approved drugs from the market. Early identification of cholestatic drugs is difficult due to the complex mechanisms involved. In order to develop a strategy for mechanism-based risk assessment of cholestatic drugs, we analyzed gene expression data obtained from the livers of rats that had been orally administered with 12 known cholestatic compounds repeatedly for 28 days at three dose levels. Qualitative analyses were performed using two statistical approaches (hierarchical clustering and principle component analysis), in addition to pathway analysis. The transcriptional benchmark dose (tBMD) and tBMD 95% lower limit (tBMDL) were used for quantitative analyses, which revealed three compound sub-groups that produced different types of differential gene expression; these groups of genes were mainly involved in inflammation, cholesterol biosynthesis, and oxidative stress. Furthermore, the tBMDL values for each test compound were in good agreement with the relevant no observed adverse effect level. These results indicate that our novel strategy for drug safety evaluation using mechanism-based classification and tBMDL would facilitate the application of toxicogenomics for risk assessment of cholestatic DILI.

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.

Unraveling the different toxic effect of flufenoxuron on the thyroid endocrine system of the Mongolia racerunner (Eremias Argus) at different stages

Flufenoxuron is a widely used pesticide to inhibit the synthesis of chitin during insect development and its effect on the growth of lizards has been little addressed. The hypothalamus-pituitary-thyroid (HPT) axis plays an important role on the development of lizards. In this study, the lizards at different development stages (proliferation and resting stages) were exposed to flufenoxuron for 21 days. The plasma thyroid hormone levels, thyroid gland histopathology and expression profiles of thyroid hormone receptors (trα, trβ), deiodinases (dio1, dio2), and transthyretin (ttr) genes were measured to evaluated the toxic effect of flufenoxuron on the HPT axis at different stages. The flufenoxuron exposure showed more seriously effect on the triiodothyronine (T3) level at resting phase than that at proliferation stage. The follicle epithelium cell height in the thyroid was only significantly increased when the exposed male lizards were at proliferation stage. The alteration of HPT axis-related genes expression was gender and tissue dependent after flufenoxuron treatment. The lizards exposed to flufenoxuron showed that the trα, trβ, dio1, dio2, and ttr genes in the female liver were more sensitive at the proliferation stage than that at the resting stage. In the male brain, the expressions of trα, trβ, dio1, and dio2 gene were significant decreased at proliferation stage while significant increased at resting stage after flufenoxuron exposure. Therefore, the thyroid endocrine system of lizards could be affected by the flufenoxuron exposure and the different development stage should also be considered when study the toxic effect of contaminants on the lizards.

Influence of ripeness and maceration of the grapes on levels of furan and carbonyl compounds in wine - Simultaneous quantitative determination and assessment of the exposure risk to these compounds

The validated method based on the use of headspace solid phase microextraction (HS-SPME) coupled with the comprehensive two-dimensional gas chromatography with time-of-flight mass spectrometric detection (GC×GC/TOFMS) proved to be appropriate for this first simultaneous quantitative determination of six toxic compounds (formaldehyde, acetaldehyde, ethyl carbamate, furan, furfural and acrolein) found in wines. Acetaldehyde and acrolein coeluted with other wine compounds, which indicated that difficulties could arise if only one-dimensional gas chromatography was used for the determination of these compounds. The advancement of the ripeness degree and increasing the grape maceration time seems to result in higher concentrations of toxic compounds. The exposure to furan, acrolein and ethyl carbamate through wine consumption may pose risks to consumer health, since calculated MOE values were lower than 10,000.

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.

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.

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.

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.