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.

Hepatic genotoxicity and toxicogenomic responses in Muta™Mouse males treated with dibenz[a,h]anthracene.

Dibenz[a,h]anthracene (DB[a,h]A) is a polycyclic aromatic hydrocarbon that is a by-product of combustion and a potent carcinogen. Few studies have investigated the effects of DB[a,h]A on mRNA and microRNA expression to dissect the mechanisms involved in carcinogenesis. In this study, mature male mice (Muta(™)Mouse) were exposed to 6.25, 12.5 and 25mg/kg/day DB[a,h]A by oral gavage for 28 consecutive days. Results were compared with mice similarly exposed to benzo[a]pyrene (B[a]P) in our previous work. Liver DNA adduct levels and lacZ mutant frequency increased dose dependently for both chemicals. Aryl hydrocarbon receptor (AhR) potency was greater for DB[a,h]A than B[a]P using the chemical-activated luciferase expression assay. Microarray analysis revealed 19 up-regulated and 22 down-regulated genes (false discovery rate-adjusted P ≤ 0.05; fold change ≥ 1.5) following treatment with 6.25 mg/kg/day DB[a,h]A. Thirteen transcripts were up-regulated and 32 down-regulated in the 12.5mg/kg/day group. The 25mg/kg/day dose had major effects on mRNA expression with 135 up-regulated and 104 down-regulated genes. Overall, perturbations were greater for DB[a,h]A than for B[a]P; in vitro chemical-activated luciferase expression supports that this may be driven by the AhR. Many of the DB[a,h]A-affected genes are implicated in cancer and are essential in vital biological functions including circadian rhythm, glucose metabolism, lipid metabolism, immune response, cell cycle and apoptosis. Although a number of functional groups were similarly affected by B[a]P and DB[a,h]A, in general the responses generated by each chemical were quite distinct. Commonalities included a DNA damage response leading to induction of cell cycle arrest and apoptosis in both Tp53-dependent and Tp53-independent manners. MicroRNA expression was identical for both chemicals, with only miR-34a showing a dose-dependent increase in treated mice.

Increased transcript levels and kinetic function of pyruvate kinase during severe dehydration in aestivating African clawed frogs, Xenopus laevis.

The African clawed frog, Xenopus laevis, can withstand extremely arid conditions through aestivation, resulting in dehydration and urea accumulation. Aestivating X. laevis reduce their metabolic rate, and rely on anaerobic glycolysis to meet reduced ATP demands. The present study investigated how severe dehydration affected the transcript levels, kinetic profile, and phosphorylation state of the key glycolytic enzyme pyruvate kinase (PK) in the liver and skeletal muscle of X. laevis. Compared to control frogs, severely dehydrated frogs showed an increase in the transcript abundance of both liver and muscle isoforms of PK. While the kinetics of muscle PK did not differ between dehydrated and control frogs, PK from the liver of dehydrated frogs had a lower Km for phosphoenolpyruvate (PEP) (38%), a lower Ka for fructose-1,6-bisphosphate (F1,6P2) (32%), and a greater activation of PK via F1,6P2 (1.56-fold). PK from dehydrated frogs also had a lower phosphorylation-state (25%) in comparison to the enzyme from control frogs in the liver. Experimental manipulation of the phosphorylation-state of liver PK taken from control frogs by endogenous protein phosphatases resulted in decreased phosphorylation, and a similar kinetic profile as seen in dehydrated frogs. The physiological consequence of dehydration-induced PK modification appears to adjust PK function to remain active during a metabolically depressed state. This study provides evidence for the maintenance of PK activity through elevated mRNA levels and a dephosphorylation event which activates frog liver PK in the dehydrated state in order to facilitate the production of ATP via anaerobic glycolysis.

The regulation of heat shock proteins in response to dehydration in Xenopus laevis.

African clawed frogs (Xenopus laevis) endure bouts of severe drought in their natural habitats and survive the loss of approximately 30% of total body water due to dehydration. To investigate molecular mechanisms employed by X. laevis during periods of dehydration, the heat shock protein response, a vital component of the cytoprotective stress response, was characterized. Using western immunoblotting and multiplex technology, the protein levels of HSP27, HSP40, HSP60, HSP70, HSC70, and HSP90 were quantified in the liver, skeletal muscle, kidney, lung, and testes from control frogs and those that underwent medium or high dehydration (~16 or ~30% loss of total body water). Dehydration increased HSP27 (1.45-1.65-fold) in the kidneys and lungs, and HSP40 (1.39-2.50-fold) in the liver, testes, and skeletal muscle. HSP60 decreased in response to dehydration (0.43-0.64 of control) in the kidneys and lungs. HSP70 increased in the liver, lungs, and testes (1.39-1.70-fold) during dehydration, but had a dynamic response in the kidneys (levels increased 1.57-fold with medium dehydration, but decreased to 0.56 of control during high dehydration). HSC70 increased in the liver and kidneys (1.20-1.36-fold), but decreased in skeletal muscle (0.27-0.55 of control) during dehydration. Lastly, HSP90 was reduced in the kidney, lung, and skeletal muscle (0.39-0.69 of control) in response to dehydration, but rose in the testes (1.30-fold). Overall, the results suggest a dynamic tissue-specific heat shock protein response to whole body dehydration in X. laevis.

Hepatic mRNA, microRNA, and miR-34a-target responses in mice after 28 days exposure to doses of benzo(a)pyrene that elicit DNA damage and mutation.

Benzo(a)pyrene (BaP) is a mutagenic carcinogen that is ubiquitous in our environment. To better understand the toxic effects of BaP and to explore the relationship between toxicity and toxicogenomics profiles, we assessed global mRNA and microRNA (miRNA) expression in Muta™Mouse. Adult male mice were exposed by oral gavage to 25, 50, and 75 mg/kg/day BaP for 28 days. Liver tissue was collected 3 days following the last treatment. Initially, we established that exposure to BaP led to the formation of hepatic DNA adducts and mutations in the lacZ transgene of the Muta™Mouse. We then analyzed hepatic gene expression profiles. Microarray analysis of liver samples revealed 134 differentially expressed transcripts (adjusted P < 0.05; fold changes > 1.5). The mRNAs most affected were involved in xenobiotic metabolism, immune response, and the downstream targets of p53. In this study, we found a significant 2.0 and 3.6-fold increase following exposure to 50 and 75 mg/kg/day BaP, respectively, relative to controls for miR-34a. This miRNA is involved in p53 response. No other significant changes in miRNAs were observed. The protein levels of five experimentally confirmed miR-34a targets were examined, and no major down-regulation was present. The results suggest that liver miRNAs are largely unresponsive to BaP doses that cause both DNA adducts and mutations. In summary, the validated miRNA and mRNA expression profiles following 28 day BaP exposure reflect a DNA damage response and effects on the cell cycle, consistent with the observed increases in DNA adducts and mutations.