Persistent gene expression and DNA methylation alterations linked to carcinogenic effects of dichloroacetic acid.

Background: Mechanistic understanding of transient exposures that lead to adverse health outcomes will enhance our ability to recognize biological signatures of disease. Here, we measured the transcriptomic and epigenomic alterations due to exposure to the metabolic reprogramming agent, dichloroacetic acid (DCA). Previously, we showed that exposure to DCA increased liver tumor incidence in B6C3F1 mice after continuous or early life exposures significantly over background level. Methods: Using archived formalin-fixed liver samples, we utilized modern methodologies to measure gene expression and DNA methylation levels to link to previously generated phenotypic measures. Gene expression was measured by targeted RNA sequencing (TempO-seq 1500+ toxicity panel: 2754 total genes) in liver samples collected from 10-, 32-, 57-, and 78-week old mice exposed to deionized water (controls), 3.5 g/L DCA continuously in drinking water ("Direct" group), or DCA for 10-, 32-, or 57-weeks followed by deionized water until sample collection ("Stop" groups). Genome-scaled alterations in DNA methylation were measured by Reduced Representation Bisulfite Sequencing (RRBS) in 78-week liver samples for control, Direct, 10-week Stop DCA exposed mice. Results: Transcriptomic changes were most robust with concurrent or adjacent timepoints after exposure was withdrawn. We observed a similar pattern with DNA methylation alterations where we noted attenuated differentially methylated regions (DMRs) in the 10-week Stop DCA exposure groups compared to the Direct group at 78-weeks. Gene pathway analysis indicated cellular effects linked to increased oxidative metabolism, a primary mechanism of action for DCA, closer to exposure windows especially early in life. Conversely, many gene signatures and pathways reversed patterns later in life and reflected more pro-tumorigenic patterns for both current and prior DCA exposures. DNA methylation patterns correlated to early gene pathway perturbations, such as cellular signaling, regulation and metabolism, suggesting persistence in the epigenome and possible regulatory effects. Conclusion: Liver metabolic reprogramming effects of DCA interacted with normal age mechanisms, increasing tumor burden with both continuous and prior DCA exposure in the male B6C3F1 rodent model.

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.

Vinyl chloride enhances high-fat diet-induced proteome alterations in the mouse pancreas related to metabolic dysfunction.

Alterations in physiological processes in pancreas have been associated with various metabolic dysfunctions and can result from environmental exposures, such as chemicals and diet. It was reported that environmental vinyl chloride (VC) exposure, a common industrial organochlorine and environmental pollutant, significantly exacerbated metabolic-related phenotypes in mice fed concurrently with high-fat diet (HFD) but not low-fat diet (LFD). However, little is known about the role of the pancreas in this interplay, especially at a proteomic level. The present study was undertaken to examine the protein responses to VC exposure in pancreas tissues of C57BL/6J mice fed LFD or HFD, with focus on the investigation of protein expression and/or phosphorylation levels of key protein biomarkers of carbohydrate, lipid, and energy metabolism, oxidative stress and detoxification, insulin secretion and regulation, cell growth, development, and communication, immunological responses and inflammation, and biomarkers of pancreatic diseases and cancers. We found that the protein alterations may indicate diet-mediated susceptibility in mouse pancreas induced by HFD to concurrent exposure of low levels of inhaled VC. These proteome biomarkers may lead to a better understanding of pancreas-mediated adaptive or adverse response and susceptibility to metabolic disease.

Circulating MicroRNAs, Polychlorinated Biphenyls, and Environmental Liver Disease in the Anniston Community Health Survey.

BACKGROUND: Polychlorinated biphenyl (PCB) exposures have been associated with liver injury in human cohorts, and steatohepatitis with liver necrosis in model systems. MicroRNAs (miRs) maintain cellular homeostasis and may regulate the response to environmental stress. OBJECTIVES: We tested the hypothesis that specific miRs are associated with liver disease and PCB exposures in a residential cohort. METHODS: Sixty-eight targeted hepatotoxicity miRs were measured in archived serum from 734 PCB-exposed participants in the cross-sectional Anniston Community Health Survey. Necrotic and other liver disease categories were defined by serum keratin 18 (K18) biomarkers. Associations were determined between exposure biomarkers (35 ortho-substituted PCB congeners) and disease biomarkers (highly expressed miRs or previously measured cytokines), and Ingenuity Pathway Analysis was performed. RESULTS: The necrotic liver disease category was associated with four up-regulated miRs (miR-99a-5p, miR-122-5p, miR-192-5p, and miR-320a) and five down-regulated miRs (let-7d-5p, miR-17-5p, miR-24-3p, miR-197-3p, and miR-221-3p). Twenty-two miRs were associated with the other liver disease category or with K18 measurements. Eleven miRs were associated with 24 PCBs, most commonly congeners with anti-estrogenic activities. Most of the exposure-associated miRs were associated with at least one serum hepatocyte death, pro-inflammatory cytokine or insulin resistance bioarker, or with both. Within each biomarker category, associations were strongest for the liver-specific miR-122-5p. Pathways of liver toxicity that were identified included inflammation/hepatitis, hyperplasia/hyperproliferation, cirrhosis, and hepatocellular carcinoma. Tumor protein p53 and tumor necrosis factor α were well integrated within the top identified networks. DISCUSSION: These results support the human hepatotoxicity of environmental PCB exposures while elucidating potential modes of PCB action. The MiR-derived liquid liver biopsy represents a promising new technique for environmental hepatology cohort studies. https://doi.org/10.1289/EHP9467.

Genomic comparisons between hepatocarcinogenic and non-hepatocarcinogenic organophosphate insecticides in the mouse liver.

Short-term biomarkers of toxicity have an increasingly important role in the screening and prioritization of new chemicals. In this study, we examined early indicators of liver toxicity for three reference organophosphate (OP) chemicals, which are among the most widely used insecticides in the world. The OP methidathion was previously shown to increase the incidence of liver toxicity, including hepatocellular tumors, in male mice. To provide insights into the adverse outcome pathway (AOP) that underlies these tumors, effects of methidathion in the male mouse liver were examined after 7 and 28 day exposures and compared to those of two other OPs that either do not increase (fenthion) or possibly suppress liver cancer (parathion) in mice. None of the chemicals caused increases in liver weight/body weight or histopathological changes in the liver. Parathion decreased liver cell proliferation after 7 and 28 days while the other chemicals had no effects. There was no evidence for hepatotoxicity in any of the treatment groups. Full-genome microarray analysis of the livers from the 7 and 28 day treatments demonstrated that methidathion and fenthion regulated a large number of overlapping genes, while parathion regulated a unique set of genes. Examination of cytochrome P450 enzyme activities and use of predictive gene expression biomarkers found no consistent evidence for activation of AhR, CAR, PXR, or PPARα. Parathion suppressed the male-specific gene expression pattern through STAT5b, similar to genetic and dietary conditions that decrease liver tumor incidence in mice. Overall, these findings indicate that methidathion causes liver cancer by a mechanism that does not involve common mechanisms of liver cancer induction.

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.

Transplacental arsenic exposure produced 5-methylcytosine methylation changes and aberrant microRNA expressions in livers of male fetal mice.

Arsenic is a known human carcinogen. Early-life exposure to inorganic arsenic induces tumors in humans and in C3H mice. We hypothesized that arsenic exposure in utero may induce epigenetic changes at the level of DNA methylation and miRNA alterations that could lead to greater postnatal susceptibility to cancer. To test this hypothesis, pregnant C3H mice were given sodium arsenite at doses known to cause liver cancer (42.5 and 85 ppm in the drinking water) from gestation day 8-19, and the livers from male fetal mice were collected for analysis. The antibody against 5-methylcytosine was used to perform chromatin-immunoprecipitation coupled with sequencing (ChIP-Seq) to determine genome-wide methylation alterations. In utero arsenic exposure produced global DNA hypomethylation and an array of gene-specific DNA methylation changes, including hypomethylation of Cyclin D1 and hypermethylation of Tp53. Illumina Correlation Engine analysis revealed 260 methylation alterations that would affect 143 microRNAs. MicroRNA array further revealed 140 aberrantly expressed miRNAs out of the 718 miRNAs. The increased expression of miR-205, miR-203, miR-215, miR-34a, and decreased expression of miR-217 were confirmed by qPCR. Comparison of the methylation changes to those of microarray analyses indicates little if any correspondence between gene methylation and gene expression. The increased expression of Xist, Prrc2, Krit1, Nish, and decreased expression of Prss2, Spp1, Col1a2, and Lox were confirmed by qPCR. In summary, in utero arsenic exposure induced global alterations in DNA methylation and aberrant miRNA expression that might contribute to adult adverse outcomes including liver cancer.

Activation of Nrf2 in the liver is associated with stress resistance mediated by suppression of the growth hormone-regulated STAT5b transcription factor.

The transcription factor Nrf2 (encoded by Nfe2l2) induces expression of numerous detoxifying and antioxidant genes in response to oxidative stress. The cytoplasmic protein Keap1 interacts with and represses Nrf2 function. Computational approaches were developed to identify factors that modulate Nrf2 in a mouse liver gene expression compendium. Forty-eight Nrf2 biomarker genes were identified using profiles from the livers of mice in which Nrf2 was activated genetically in Keap1-null mice or chemically by a potent activator of Nrf2 signaling. The rank-based Running Fisher statistical test was used to determine the correlation between the Nrf2 biomarker genes and a test set of 81 profiles with known Nrf2 activation status demonstrating a balanced accuracy of 96%. For a large number of factors examined in the compendium, we found consistent relationships between activation of Nrf2 and feminization of the liver transcriptome through suppression of the male-specific growth hormone (GH)-regulated transcription factor STAT5b. The livers of female mice exhibited higher Nrf2 activation than male mice in untreated or chemical-treated conditions. In male mice, Nrf2 was activated by treatment with ethinyl estradiol, whereas in female mice, Nrf2 was suppressed by treatment with testosterone. Nrf2 was activated in 5 models of disrupted GH signaling containing mutations in Pit1, Prop1, Ghrh, Ghrhr, and Ghr. Out of 59 chemical treatments that activated Nrf2, 36 exhibited STAT5b suppression in the male liver. The Nrf2-STAT5b coupling was absent in in vitro comparisons of chemical treatments. Treatment of male and female mice with 11 chemicals that induce oxidative stress led to activation of Nrf2 to greater extents in females than males. The enhanced basal and inducible levels of Nrf2 activation in females relative to males provides a molecular explanation for the greater resistance often seen in females vs. males to age-dependent diseases and chemical-induced toxicity.

Metabolic Disruption Early in Life is Associated With Latent Carcinogenic Activity of Dichloroacetic Acid in Mice.

Early-life environmental factors can influence later-life susceptibility to cancer. Recent evidence suggests that metabolic pathways may mediate this type of latency effect. Previously, we reported that short-term exposure to dichloroacetic acid (DCA) increased liver cancer in mice 84 weeks after exposure was stopped. Here, we evaluated time course dynamics for key events related to this effect. This study followed a stop-exposure design in which 28-day-old male B6C3F1 mice were given the following treatments in drinking water for up to 93 weeks: deionized water (dH2O, control); 3.5 g/l DCA continuously; or 3.5 g/l DCA for 4-52 weeks followed by dH2O. Effects were evaluated at eight interim time points. A short-term biomarker study was used to evaluate DCA effects at 6, 15, and 30 days. Liver tumor incidence was higher in all DCA treatment groups, including carcinomas in 82% of mice previously treated with DCA for only 4 weeks. Direct effects of DCA in the short-term study included decreased liver cell proliferation and marked mRNA changes related to mitochondrial dysfunction and altered cell metabolism. However, all observed short-term effects of DCA were ultimately reversible, and prior DCA treatment did not affect liver cell proliferation, apoptosis, necrosis, or DNA sequence variants with age. Key intermediate events resulting from transient DCA exposure do not fit classical cytotoxic, mitogenic, or genotoxic modes of action for carcinogenesis, suggesting a distinct mechanism associated with early-life metabolic disruption.

From the Cover: Genomic Effects of Androstenedione and Sex-Specific Liver Cancer Susceptibility in Mice.

Current strategies for predicting carcinogenic mode of action for nongenotoxic chemicals are based on identification of early key events in toxicity pathways. The goal of this study was to evaluate short-term key event indicators resulting from exposure to androstenedione (A4), an androgen receptor agonist and known liver carcinogen in mice. Liver cancer is more prevalent in men compared with women, but androgen-related pathways underlying this sex difference have not been clearly identified. Short-term hepatic effects of A4 were compared with reference agonists of the estrogen receptor (ethinyl estradiol, EE) and glucocorticoid receptor (prednisone, PRED). Male B6C3F1 mice were exposed for 7 or 28 days to A4, EE, or PRED. EE increased and PRED suppressed hepatocyte proliferation, while A4 had no detectable effects. In a microarray analysis, EE and PRED altered >3000 and >670 genes, respectively, in a dose-dependent manner, whereas A4 did not significantly alter any genes. Gene expression was subsequently examined in archival liver samples from male and female B6C3F1 mice exposed to A4 for 90 days. A4 altered more genes in females than males and did not alter expression of genes linked to activation of the mitogenic xenobiotic receptors AhR, CAR, and PPARα in either sex. A gene expression biomarker was used to show that in female mice, the high dose of A4 activated the growth hormone-regulated transcription factor STAT5b, which controls sexually dimorphic gene expression in the liver. These findings suggest that A4 induces subtle age-related effects on STAT5b signaling that may contribute to the higher risk of liver cancer in males compared with females.

MicroRNA Biomarkers of Toxicity in Biological Matrices.

Biomarker measurements that reliably correlate with tissue injury and that can be measured within accessible biofluids offer benefits in terms of cost, time, and convenience when assessing chemical and drug-induced toxicity in model systems or human cohorts. MicroRNAs (miRNAs) have emerged in recent years as a promising new class of biomarker for monitoring toxicity. Recent enthusiasm for miRNA biomarker research has been fueled by evidence that certain miRNAs are cell-type specific and are released during injury, thus raising the possibility of using biofluid-based miRNAs as a "liquid biopsy" that may be obtained by sampling extracellular fluids. As biomarkers, miRNAs demonstrate improved stability as compared with many protein markers and sequences are largely conserved across species, simplifying analytical techniques. Recent efforts have sought to identify miRNAs that are released into accessible biofluids following xenobiotic exposure, using compounds that target specific organs. Whereas still early in the discovery phase, miRNA biomarkers will have an increasingly important role in the assessment of adverse effects of both environmental chemicals and pharmaceutical drugs. Here, we review the current findings of biofluid-based miRNAs, as well as highlight technical challenges in assessing toxicologic pathology using these biomarkers.

Developing a Gene Biomarker at the Tipping Point of Adaptive and Adverse Responses in Human Bronchial Epithelial Cells.

Determining mechanism-based biomarkers that distinguish adaptive and adverse cellular processes is critical to understanding the health effects of environmental exposures. Shifting from in vivo, low-throughput toxicity studies to high-throughput screening (HTS) paradigms and risk assessment based on in vitro and in silico testing requires utilizing toxicity pathway information to distinguish adverse outcomes from recoverable adaptive events. Little work has focused on oxidative stresses in human airway for the purposes of predicting adverse responses. We hypothesize that early gene expression-mediated molecular changes could be used to delineate adaptive and adverse responses to environmentally-based perturbations. Here, we examined cellular responses of the tracheobronchial airway to zinc (Zn) exposure, a model oxidant. Airway derived BEAS-2B cells exposed to 2-10 µM Zn2+ elicited concentration- and time-dependent cytotoxicity. Normal, adaptive, and cytotoxic Zn2+ exposure conditions were determined with traditional apical endpoints, and differences in global gene expression around the tipping point of the responses were used to delineate underlying molecular mechanisms. Bioinformatic analyses of differentially expressed genes indicate early enrichment of stress signaling pathways, including those mediated by the transcription factors p53 and NRF2. After 4 h, 154 genes were differentially expressed (p < 0.01) between the adaptive and cytotoxic Zn2+ concentrations. Nearly 40% of the biomarker genes were related to the p53 signaling pathway with 30 genes identified as likely direct targets using a database of p53 ChIP-seq studies. Despite similar p53 activation profiles, these data revealed widespread dampening of p53 and NRF2-related genes as early as 4 h after exposure at higher, unrecoverable Zn2+ exposures. Thus, in our model early increased activation of stress response pathways indicated a recoverable adaptive event. Overall, this study highlights the importance of characterizing molecular mechanisms around the tipping point of adverse responses to better inform HTS paradigms.

A Polymorphic Antioxidant Response Element Links NRF2/sMAF Binding to Enhanced MAPT Expression and Reduced Risk of Parkinsonian Disorders.

The NRF2/sMAF protein complex regulates the oxidative stress response by occupying cis-acting enhancers containing an antioxidant response element (ARE). Integrating genome-wide maps of NRF2/sMAF occupancy with disease-susceptibility loci, we discovered eight polymorphic AREs linked to 14 highly ranked disease-risk SNPs in individuals of European ancestry. Among these SNPs was rs242561, located within a regulatory region of the MAPT gene (encoding microtubule-associated protein Tau). It was consistently occupied by NRF2/sMAF in multiple experiments and its strong-binding allele associated with higher mRNA levels in cell lines and human brain tissue. Induction of MAPT transcription by NRF2 was confirmed using a human neuroblastoma cell line and a Nrf2-deficient mouse model. Most importantly, rs242561 displayed complete linkage disequilibrium with a highly protective allele identified in multiple GWASs of progressive supranuclear palsy, Parkinson's disease, and corticobasal degeneration. These observations suggest a potential role for NRF2/sMAF in tauopathies and a possible role for NRF2 pathway activators in disease prevention.

Tipping the Balance: Hepatotoxicity and the 4 Apical Key Events of Hepatic Steatosis.

Hepatic steatosis is a condition were fat accumulates in the liver and it is associated with extra-hepatic diseases related to metabolic syndrome and systemic energy metabolism. If not reversed, steatosis can progress to steatohepatitis and irreversible stages of liver disease including fibrosis, cirrhosis, hepatocellular carcinoma, and death. From a public health standpoint, identifying chemical exposures that may be factors in steatosis etiology are important for preventing hepatotoxicity and liver disease progression. It is therefore important to identify the biological events that are key for steatosis pathology mediated by chemical exposure. In this review, we give a current overview of the complex biological cascades that can disrupt lipid homeostasis in hepatocytes in the context of 4 apical key events central to hepatic lipid retention: hepatic fatty acid (FA) uptake,de novoFA and lipid synthesis, FA oxidation, and lipid efflux. Our goal is to review these key cellular events and visually summarize them using a network for application in pathway-based toxicity testing. This effort provides a foundation to improve next-generation chemical screening efforts that may be used to prevent and ultimately reverse the growing incidence of fatty liver disease in our population.

Dose and Effect Thresholds for Early Key Events in a PPARα-Mediated Mode of Action.

Current strategies for predicting adverse health outcomes of environmental chemicals are centered on early key events in toxicity pathways. However, quantitative relationships between early molecular changes in a given pathway and later health effects are often poorly defined. The goal of this study was to evaluate short-term key event indicators using qualitative and quantitative methods in an established pathway of mouse liver tumorigenesis mediated by peroxisome proliferator-activated receptor alpha (PPARα). Male B6C3F1 mice were exposed for 7 days to di (2-ethylhexyl) phthalate (DEHP), di-n-octyl phthalate (DNOP), and n-butyl benzyl phthalate (BBP), which vary in PPARα activity and liver tumorigenicity. Each phthalate increased expression of select PPARα target genes at 7 days, while only DEHP significantly increased liver cell proliferation labeling index (LI). Transcriptional benchmark dose (BMDT) estimates for dose-related genomic markers stratified phthalates according to hypothetical tumorigenic potencies, unlike BMDs for non-genomic endpoints (relative liver weights or proliferation). The 7-day BMDT values for Acot1 as a surrogate measure for PPARα activation were 29, 370, and 676 mg/kg/day for DEHP, DNOP, and BBP, respectively, distinguishing DEHP (liver tumor BMD of 35 mg/kg/day) from non-tumorigenic DNOP and BBP. Effect thresholds were generated using linear regression of DEHP effects at 7 days and 2-year tumor incidence values to anchor early response molecular indicators and a later phenotypic outcome. Thresholds varied widely by marker, from 2-fold (Pdk4 and proliferation LI) to 30-fold (Acot1) induction to reach hypothetical tumorigenic expression levels. These findings highlight key issues in defining thresholds for biological adversity based on molecular changes.

Latent carcinogenicity of early-life exposure to dichloroacetic acid in mice.

Environmental exposures occurring early in life may have an important influence on cancer risk later in life. Here, we investigated carryover effects of dichloroacetic acid (DCA), a small molecule analog of pyruvate with metabolic programming properties, on age-related incidence of liver cancer. The study followed a stop-exposure/promotion design in which 4-week-old male and female B6C3F1 mice received the following treatments: deionized water alone (dH2O, control); dH2O with 0.06% phenobarbital (PB), a mouse liver tumor promoter; or DCA (1.0, 2.0 or 3.5g/l) for 10 weeks followed by dH2O or PB (n = 20-30/group/sex). Pathology and molecular assessments were performed at 98 weeks of age. In the absence of PB, early-life exposure to DCA increased the incidence and number of hepatocellular tumors in male and female mice compared with controls. Significant dose trends were observed in both sexes. At the high dose level, 10 weeks of prior DCA treatment induced comparable effects (≥85% tumor incidence and number) to those seen after continuous lifetime exposure. Prior DCA treatment did not enhance or inhibit the carcinogenic effects of PB, induce persistent liver cytotoxicity or preneoplastic changes on histopathology or alter DNA sequence variant profiles within liver tumors compared with controls. Distinct changes in liver messenger RNA and micro RNA profiles associated with prior DCA treatment were not apparent at 98 weeks. Our findings demonstrate that early-life exposure to DCA may be as carcinogenic as life-long exposures, potentially via epigenetic-mediated effects related to cellular metabolism.

Comparative time course profiles of phthalate stereoisomers in mice.

More efficient models are needed to assess potential carcinogenicity hazard of environmental chemicals based on early events in tumorigenesis. Here, we investigated time course profiles for key events in an established cancer mode of action. Using a case study approach, we evaluated two reference phthalates, di(2-ethylhexyl) phthalate (DEHP) and its stereoisomer di-n-octyl phthalate (DNOP), across the span of a two-year carcinogenicity bioassay. Male B6C3F1 mice received diets with no phthalate added (control), DEHP at 0.12, 0.60, or 1.20%, or DNOP at 0.10, 0.50, or 1.00% (n = 80-83/group) for up to 104 weeks with six interim evaluations starting at week 4. Mean phthalate doses were 139, 845, and 3147 mg/kg/day for DEHP and 113, 755, and 1281 mg/kg/day for DNOP groups, respectively. Incidence and number of hepatocellular tumors (adenoma and/or carcinoma) were greater at ≥ 60 weeks for all DEHP groups with time and dose trends, whereas DNOP had no significant effects. Key events supported a peroxisome proliferator-activated receptor alpha (PPARα) mode of action for DEHP, with secondary cytotoxicity at the high dose, whereas DNOP induced modest increases in PPARα activity without proliferative or cytotoxic effects. Threshold estimates for later tumorigenic effects were identified at week 4 for relative liver weight (+24%) and PPARα activity (+79%) relative to the control group. Benchmark doses (BMDs) for these measures at week 4 clearly distinguished DEHP and DNOP and showed strong concordance with values at later time points and tumorigenic BMDs. Other target sites included testis and kidney, which showed degenerative changes at higher doses of DEHP but not DNOP. Our results highlight marked differences in the chronic toxicity profiles of structurally similar phthalates and demonstrate quantitative relationships between early bioindicators and later tumor outcomes.

Novel hematopoietic target genes in the NRF2-mediated transcriptional pathway.

Nuclear factor- (erythroid-derived 2) like 2 (NFE2L2, NRF2) is a key transcriptional activator of the antioxidant response pathway and is closely related to erythroid transcription factor NFE2. Under oxidative stress, NRF2 heterodimerizes with small Maf proteins and binds cis-acting enhancer sequences found near oxidative stress response genes. Using the dietary isothiocyanate sulforaphane (SFN) to activate NRF2, chromatin immunoprecipitation sequencing (ChIP-seq) identified several hundred novel NRF2-mediated targets beyond its role in oxidative stress. Activated NRF2 bound the antioxidant response element (ARE) in promoters of several known and novel target genes involved in iron homeostasis and heme metabolism, including known targets FTL and FTH1, as well as novel binding in the globin locus control region. Five novel NRF2 target genes were chosen for followup: AMBP, ABCB6, FECH, HRG-1 (SLC48A1), and TBXAS1. SFN-induced gene expression in erythroid K562 and lymphoid cells were compared for each target gene. NRF2 silencing showed reduced expression in lymphoid, lung, and hepatic cells. Furthermore, stable knockdown of NRF2 negative regulator KEAP1 in K562 cells resulted in increased NQO1, AMBP, and TBXAS1 expression. NFE2 binding sites in K562 cells revealed similar binding profiles as lymphoid NRF2 sites in all potential NRF2 candidates supporting a role for NRF2 in heme metabolism and erythropoiesis.

Genetic variation and antioxidant response gene expression in the bronchial airway epithelium of smokers at risk for lung cancer.

Prior microarray studies of smokers at high risk for lung cancer have demonstrated that heterogeneity in bronchial airway epithelial cell gene expression response to smoking can serve as an early diagnostic biomarker for lung cancer. As a first step in applying functional genomic analysis to population studies, we have examined the relationship between gene expression variation and genetic variation in a central molecular pathway (NRF2-mediated antioxidant response) associated with smoking exposure and lung cancer. We assessed global gene expression in histologically normal airway epithelial cells obtained at bronchoscopy from smokers who developed lung cancer (SC, n = 20), smokers without lung cancer (SNC, n = 24), and never smokers (NS, n = 8). Functional enrichment analysis showed that the NRF2-mediated, antioxidant response element (ARE)-regulated genes, were significantly lower in SC, when compared with expression levels in SNC. Importantly, we found that the expression of MAFG (a binding partner of NRF2) was correlated with the expression of ARE genes, suggesting MAFG levels may limit target gene induction. Bioinformatically we identified single nucleotide polymorphisms (SNPs) in putative ARE genes and to test the impact of genetic variation, we genotyped these putative regulatory SNPs and other tag SNPs in selected NRF2 pathway genes. Sequencing MAFG locus, we identified 30 novel SNPs and two were associated with either gene expression or lung cancer status among smokers. This work demonstrates an analysis approach that integrates bioinformatics pathway and transcription factor binding site analysis with genotype, gene expression and disease status to identify SNPs that may be associated with individual differences in gene expression and/or cancer status in smokers. These polymorphisms might ultimately contribute to lung cancer risk via their effect on the airway gene expression response to tobacco-smoke exposure.

Probing the functional impact of sequence variation on p53-DNA interactions using a novel microsphere assay for protein-DNA binding with human cell extracts.

The p53 tumor suppressor regulates its target genes through sequence-specific binding to DNA response elements (REs). Although numerous p53 REs are established, the thousands more identified by bioinformatics are not easily subjected to comparative functional evaluation. To examine the relationship between RE sequence variation -- including polymorphisms -- and p53 binding, we have developed a multiplex format microsphere assay of protein-DNA binding (MAPD) for p53 in nuclear extracts. Using MAPD we measured sequence-specific p53 binding of doxorubicin-activated or transiently expressed p53 to REs from established p53 target genes and p53 consensus REs. To assess the sensitivity and scalability of the assay, we tested 16 variants of the p21 target sequence and a 62-multiplex set of single nucleotide (nt) variants of the p53 consensus sequence and found many changes in p53 binding that are not captured by current computational binding models. A group of eight single nucleotide polymorphisms (SNPs) was examined and binding profiles closely matched transactivation capability tested in luciferase constructs. The in vitro binding characteristics of p53 in nuclear extracts recapitulated the cellular in vivo transactivation capabilities for eight well-established human REs measured by luciferase assay. Using a set of 26 bona fide REs, we observed distinct binding patterns characteristic of transiently expressed wild type and mutant p53s. This microsphere assay system utilizes biologically meaningful cell extracts in a multiplexed, quantitative, in vitro format that provides a powerful experimental tool for elucidating the functional impact of sequence polymorphism and protein variation on protein/DNA binding in transcriptional networks.