Search

High-throughput transcriptomics (HTTr) is increasingly being used to identify molecular targets of chemicals that can be linked to adverse outcomes. Cell proliferation (CP) is an important key event in chemical carcinogenesis. Here, we describe the construction and characterization of a gene expression biomarker that is predictive of the CP status in human and rodent tissues. The biomarker was constructed from 30 genes known to be increased in expression in prostate cancers relative to surrounding tissues and in cycling human MCF-7 cells after estrogen receptor (ER) agonist exposure. Using a large compendium of gene expression profiles to test utility, the biomarker could identify increases in CP in (i) 308 out of 367 tumor vs. normal surrounding tissue comparisons from 6 human organs, (ii) MCF-7 cells after activation of ER, (iii) after partial hepatectomy in mice and rats, and (iv) the livers of mice and rats after exposure to nongenotoxic hepatocarcinogens. The biomarker identified suppression of CP (i) under conditions of p53 activation by DNA damaging agents in human cells, (ii) in human A549 lung cells exposed to therapeutic anticancer kinase inhibitors (dasatinib, nilotnib), and (iii) in the mouse liver when comparing high levels of CP at birth to the low background levels in the adult. The responses using the biomarker were similar to those observed using conventional markers of CP including PCNA, Ki67, and BrdU labeling. The CP biomarker will be a useful tool for interpretation of HTTr data streams to identify CP status after exposure to chemicals in human cells or in rodent tissues., (Published by Oxford University Press on behalf of the Society of Toxicology 2024.)

There is growing recognition across broad sectors of the scientific community that use of genomic biomarkers has the potential to reduce the need for conventional rodent carcinogenicity studies of industrial chemicals, agrochemicals, and pharmaceuticals through a weight-of-evidence approach. These biomarkers fall into 2 major categories: (1) sets of gene transcripts that can identify distinct tumorigenic mechanisms of action; and (2) cancer driver gene mutations indicative of rapidly expanding growth-advantaged clonal cell populations. This call-to-action article describes a collaborative approach launched to develop and qualify biomarker gene expression panels that measure widely accepted molecular pathways linked to tumorigenesis and their activation levels to predict tumorigenic doses of chemicals from short-term exposures. Growing evidence suggests that application of such biomarker panels in short-term exposure rodent studies can identify both tumorigenic hazard and tumorigenic activation levels for chemical-induced carcinogenicity. In the future, this approach will be expanded to include methodologies examining mutations in key cancer driver gene mutation hotspots as biomarkers of both genotoxic and nongenotoxic chemical tumor risk. Analytical, technical, and biological validation studies of these complementary genomic tools are being undertaken by multisector and multidisciplinary collaborative teams within the Health and Environmental Sciences Institute. Success from these efforts will facilitate the transition from current heavy reliance on conventional 2-year rodent carcinogenicity studies to more rapid animal- and resource-sparing approaches for mechanism-based carcinogenicity evaluation supporting internal and regulatory decision-making., (© The Author(s) 2022. Published by Oxford University Press on behalf of the Society of Toxicology.)

Environmental exposure to metals is known to cause a number of human toxicities including cancer. Metal-responsive transcription factor 1 (MTF-1) is an important component of metal regulation systems in mammalian cells. Here, we describe a novel method to identify chemicals that activate MTF-1 based on microarray profiling data. MTF-1 biomarker genes were identified that exhibited consistent, robust expression across 10 microarray comparisons examining the effects of metals (zinc, nickel, lead, arsenic, mercury, and silver) on gene expression in human cells. A subset of the resulting 81 biomarker genes was shown to be altered by knockdown of the MTF1 gene including metallothionein family members and a zinc transporter. The ability to correctly identify treatment conditions that activate MTF-1 was determined by comparing the biomarker to microarray comparisons from cells exposed to reference metal activators of MTF-1 using the rank-based Running Fisher algorithm. The balanced accuracy for prediction was 93%. The biomarker was then used to identify organic chemicals that activate MTF-1 from a compendium of 11 725 human gene expression comparisons representing 2582 chemicals. There were 700 chemicals identified that included those known to interact with cellular metals, such as clioquinol and disulfiram, as well as a set of novel chemicals. All nine of the novel chemicals selected for validation were confirmed to activate MTF-1 biomarker genes in MCF-7 cells and to lesser extents in MTF1-null cells by qPCR and targeted RNA-Seq. Overall, our work demonstrates that the biomarker for MTF-1 coupled with the Running Fisher test is a reliable strategy to identify novel chemical modulators of metal homeostasis using gene expression profiling.

Traditional methods for cancer risk assessment are resource-intensive, retrospective, and not feasible for the vast majority of environmental chemicals. In this study, we investigated whether quantitative genomic data from short-term studies may be used to set protective thresholds for potential tumorigenic effects. We hypothesized that gene expression biomarkers measuring activation of the key early events in established pathways for rodent liver cancer exhibit cross-chemical thresholds for tumorigenesis predictive for liver cancer risk. We defined biomarker thresholds for 6 major liver cancer pathways using training sets of chemicals with short-term genomic data (3-29 days of exposure) from the TG-GATES (n = 77 chemicals) and DrugMatrix (n = 86 chemicals) databases and then tested these thresholds within and between datasets. The 6 pathway biomarkers represented genotoxicity, cytotoxicity, and activation of xenobiotic, steroid, and lipid receptors (aryl hydrocarbon receptor, constitutive activated receptor, estrogen receptor, and peroxisome proliferator-activated receptor α). Thresholds were calculated as the maximum values derived from exposures without detectable liver tumor outcomes. We identified clear response values that were consistent across training and test sets. Thresholds derived from the TG-GATES training set were highly predictive (97%) in a test set of independent chemicals, whereas thresholds derived from the DrugMatrix study were 96%-97% predictive for the TG-GATES study. Threshold values derived from an abridged gene list (2/biomarker) also exhibited high predictive accuracy (91%-94%). These findings support the idea that early genomic changes can be used to establish threshold estimates or "molecular tipping points" that are predictive of later-life health outcomes., (Published by Oxford University Press on behalf of the Society of Toxicology 2020. This work is written by US Government employees and is in the public domain in the US.)

Chemical-induced liver cancer occurs in rodents through well-characterized adverse outcome pathways. We hypothesized that measurement of the 6 most common molecular initiating events (MIEs) in liver cancer adverse outcome pathways in short-term assays using only gene expression will allow early identification of chemicals and their associated doses that are likely to be tumorigenic in the liver in 2-year bioassays. We tested this hypothesis using transcript data from a rat liver microarray compendium consisting of 2013 comparisons of 146 chemicals administered at doses with previously established effects on rat liver tumor induction. Five MIEs were measured using previously characterized gene expression biomarkers composed of gene sets predictive for genotoxicity and activation of 1 or more xenobiotic receptors (aryl hydrocarbon receptor, constitutive activated receptor, estrogen receptor, and peroxisome proliferator-activated receptor α). Because chronic injury can be important in tumorigenesis, we also developed a biomarker for cytotoxicity that had a 96% balanced accuracy. Characterization of the genes in each biomarker set using the unsupervised TXG-MAP network model demonstrated that the genes were associated with distinct functional coexpression modules. Using the Toxicological Priority Index to rank chemicals based on their ability to activate the MIEs showed that chemicals administered at tumorigenic doses clearly gave the highest ranked scores. Balanced accuracies using thresholds derived from either TG-GATES or DrugMatrix data sets to predict tumorigenicity in independent sets of chemicals were up to 93%. These results show that a MIE-directed approach using only gene expression biomarkers could be used in short-term assays to identify chemicals and their doses that cause tumors., (Published by Oxford University Press on behalf of the Society of Toxicology 2020. This work is written by US Government employees and is in the public domain in the US.)

Exposure to environmentally relevant chemicals that activate the xenobiotic receptors aryl hydrocarbon receptor (AhR), constitutive androstane receptor (CAR), and peroxisome proliferator-activated receptor alpha (PPARα) in rodent test systems often leads to increases in oxidative stress (OS) that contributes to liver cancer induction. We hypothesized that activation of the oxidant-induced transcription factor Nrf2 could be used as a surrogate endpoint for increases in OS. We examined the relationships between activation of xenobiotic receptors and Nrf2 using previously characterized gene expression biomarkers that accurately predict modulation. Using a correlation approach (Running Fisher Test), the biomarkers were compared with microarray profiles in a mouse liver gene expression compendium. Out of the 163 chemicals examined, 47% from 53 studies activated Nrf2. We found consistent coupling between CAR and Nrf2 activation. Out of the 41 chemicals from 32 studies that activated CAR, 90% also activated Nrf2. CAR was activated earlier and at lower doses than Nrf2, indicating CAR activation preceded Nrf2 activation. Nrf2 activation by 2 CAR activators was abolished in CAR-null mice. We hypothesized that Nrf2 is activated by reactive oxygen species from the increased activity of enzymes encoded by Cyp2b family members. However, Nrf2 was similarly activated in the livers of both TCPOBOP-treated wild-type and Cyp2b9/10/13-null mice. This study provides evidence that Nrf2 activation (1) often occurs after exposure to xenobiotic chemicals, (2) is tightly linked to activation of CAR, and (3) does not require induction of 3 Cyp2b genes secondary to CAR activation.

High-throughput transcriptomic (HTTr) technologies are increasingly being used to screen environmental chemicals in vitro to identify molecular targets and provide mechanistic context for regulatory testing. Here, we describe the development and validation of a novel gene expression biomarker to identify androgen receptor (AR)-modulating chemicals using a pattern matching method. Androgen receptor biomarker genes were identified by their consistent expression after exposure to 4 AR agonists and 4 AR antagonists and included only those genes that were regulated by AR. The 51 gene biomarker was evaluated as a predictive tool using the fold-change, rank-based Running Fisher algorithm. Using 158 comparisons from cells treated with 95 chemicals, the biomarker gave balanced accuracies for prediction of AR activation or AR suppression of 97% or 98%, respectively. The biomarker correctly classified 16 out of the 17 AR reference antagonists including those that are "weak" and "very weak". Predictions based on microarray profiles from AR-positive LAPC-4 cells treated with 28 chemicals in antagonist mode were compared with those from an AR pathway model which used 11 in vitro HT assays. The balanced accuracy for suppression was 93%. Using our approach, we identified conditions in which AR was modulated in a large collection of microarray profiles from prostate cancer cell lines including (1) constitutively active mutants or knockdown of AR, (2) decreases in availability of androgens by castration or removal from media, and (3) exposure to chemical modulators that work through indirect mechanisms including suppression of AR expression. These results demonstrate that the AR gene expression biomarker could be a useful tool in HTTr to identify AR modulators.

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., (Published by Oxford University Press on behalf of the Society of Toxicology 2017. This work is written by US Government employees and is in the public domain in the United States.)

Plasticizers with estrogenic activity, such as bisphenol A (BPA), have potential adverse health effects in humans. Due to mounting evidence of these health effects, BPA is being phased out and replaced by other bisphenol variants in "BPA-free" products. We have compared estrogenic activity of BPA with 6 bisphenol analogues [bisphenol S (BPS); bisphenol F (BPF); bisphenol AP (BPAP); bisphenol AF (BPAF); bisphenol Z (BPZ); bisphenol B (BPB)] in 3 human breast cancer cell lines. Estrogenicity was assessed (10-11-10-4 M) by cell growth in an estrogen receptor (ER)-mediated cell proliferation assay, and by the induction of estrogen response element-mediated transcription in a luciferase assay. BPAF was the most potent bisphenol, followed by BPB > BPZ ∼ BPA > BPF ∼ BPAP > BPS. The addition of ICI 182,780 antagonized the activation of ERs. Data mining of ToxCast high-throughput screening assays confirm our results but also show divergence in the sensitivities of the assays. Gene expression profiles were determined in MCF-7 cells by microarray analysis. The comparison of transcriptome profile alterations resulting from BPA alternatives with an ERα gene expression biomarker further indicates that all BPA alternatives act as ERα agonists in MCF-7 cells. These results were confirmed by Illumina-based RNA sequencing. In conclusion, BPA alternatives are not necessarily less estrogenic than BPA in human breast cancer cells. BPAF, BPB, and BPZ were more estrogenic than BPA. These findings point to the importance of better understanding the risk of adverse effects from exposure to BPA alternatives, including hormone-dependent breast cancer., (© The Author 2017. Published by Oxford University Press on behalf of the Society of Toxicology.)

Recent international efforts have led to proposals for modified carcinogenicity testing paradigms based on data from shorter-term studies. The main goal of the current study was to evaluate the negative predictive value (NPV) of short-term toxicity indicators on carcinogenicity study outcomes and cancer classifications for chemicals previously reviewed by the U.S. Environmental Protection Agency (EPA). Pathology data were analyzed from over 900 acceptable 2-sex guideline subchronic (3-month) and carcinogenicity studies in the U.S. EPA Toxicity Reference Database. Chemical cancer classifications were obtained from annual reports of the U.S. EPA Office of Pesticide Programs. Histopathologic risk signals and evidence of hormonal perturbation in subchronic rat studies provided 56% NPV for any tumor outcome in the rat or mouse and 75% NPV for cancer classifications not requiring quantitative risk assessment (qRA). In comparison, lack of activity in a battery of 35 in vitro cytotoxicity assays from the U.S. EPA ToxCast library provided 49% NPV for any tumor outcome and 80% NPV for cancer classifications not requiring qRA. These findings support the idea that the absence of short-term bioactivity may provide useful information for prioritizing chemicals based on potential carcinogenic risk. Additional data streams are needed to further refine these models., (Published by Oxford University Press on behalf of the Society of Toxicology 2016. This work is written by US Government employees and is in the public domain in the US.)

Microarray profiling of chemical-induced effects is being increasingly used in medium- and high-throughput formats. Computational methods are described here to identify molecular targets from whole-genome microarray data using as an example the estrogen receptor α (ERα), often modulated by potential endocrine disrupting chemicals. ERα biomarker genes were identified by their consistent expression after exposure to 7 structurally diverse ERα agonists and 3 ERα antagonists in ERα-positive MCF-7 cells. Most of the biomarker genes were shown to be directly regulated by ERα as determined by ESR1 gene knockdown using siRNA as well as through chromatin immunoprecipitation coupled with DNA sequencing analysis of ERα-DNA interactions. The biomarker was evaluated as a predictive tool using the fold-change rank-based Running Fisher algorithm by comparison to annotated gene expression datasets from experiments using MCF-7 cells, including those evaluating the transcriptional effects of hormones and chemicals. Using 141 comparisons from chemical- and hormone-treated cells, the biomarker gave a balanced accuracy for prediction of ERα activation or suppression of 94% and 93%, respectively. The biomarker was able to correctly classify 18 out of 21 (86%) ER reference chemicals including "very weak" agonists. Importantly, the biomarker predictions accurately replicated predictions based on 18 in vitro high-throughput screening assays that queried different steps in ERα signaling. For 114 chemicals, the balanced accuracies were 95% and 98% for activation or suppression, respectively. These results demonstrate that the ERα gene expression biomarker can accurately identify ERα modulators in large collections of microarray data derived from MCF-7 cells., (Published by Oxford University Press on behalf of the Society of Toxicology 2016. This work is written by US Government employees and is in the public domain in the US.)

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., (Published by Oxford University Press on behalf of the Society of Toxicology 2015. This work is written by US Government employees and is in the public domain in the US.)

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.

Peroxisome proliferator chemicals (PPC) are thought to mediate their effects in rodents on hepatocyte growth and liver cancer through the nuclear receptor peroxisome proliferator-activated receptor (PPAR) alpha. Recent studies indicate that the plasticizer di-(2-ethylhexyl) phthalate (DEHP) increased the incidence of liver tumors in PPARalpha-null mice. We hypothesized that some PPC, including DEHP, induce transcriptional changes independent of PPARalpha but dependent on other nuclear receptors, including the constitutive-activated receptor (CAR) that mediates phenobarbital (PB) effects on hepatocyte growth and liver tumor induction. To determine the potential role of CAR in mediating effects of PPC, a meta-analysis was performed on transcript profiles from published studies in which rats and mice were exposed to PPC and compared the profiles to those produced by exposure to PB. Valproic acid, clofibrate, and DEHP in rat liver and DEHP in mouse liver induced genes, including Cyp2b family members that are known to be regulated by CAR. Examination of transcript changes by Affymetrix ST 1.0 arrays and reverse transcription-PCR in the livers of DEHP-treated wild-type, PPARalpha-null, and CAR-null mice demonstrated that (1) most (approximately 94%) of the transcriptional changes induced by DEHP were PPARalpha-dependent, (2) many PPARalpha-independent genes overlapped with those regulated by PB, (3) induction of genes Cyp2b10, Cyp3a11, and metallothionine-1 by DEHP was CAR dependent but PPARalpha-independent, and (4) induction of a number of genes (Cyp8b1, Gstm4, and Gstm7) was independent of both CAR and PPARalpha. Our results indicate that exposure to PPARalpha activators including DEHP leads to activation of multiple nuclear receptors in the rodent liver.

In order to gain insight into the effects of aging on susceptibility to environmental toxins, we characterized the expression of xenobiotic metabolizing enzymes (XMEs) from the livers of male F344 and Brown Norway (BN) rats across the adult lifespan. Using full-genome Affymetrix arrays, principal component analysis showed a clear age-dependent separation between young and old animals in both rat strains. Out of 1135 or 1435 genes altered between the old and young groups in the F344 or BN rats, 7 or 3% were XMEs and included members of the phase I, II, and III classes of genes. There was a 20 or 32% overlap in the gene expression profile between the two strains for F344 or BN, respectively. Lipid, ergosterol, alcohol, and fatty acid metabolism genes were also altered with age in both strains. Some of the genes altered by age exhibited a gender-dependent expression pattern in young adult rats, suggesting an increasingly feminized pattern of gene expression with age in male rats. To examine transcriptional responses across lifespan after challenge with a xenobiotic compound, BN rats were exposed to toluene by oral gavage. Toluene exposure decreased the expression of glutathione synthetase, and dramatically increased the number of phase III genes being downregulated. The expression of CYP2B2 and glutathione-S-transferase decreased with age but increased in all age groups after toluene exposure. Decreased ability to detoxify and transport chemicals out of the body with age could result in increased susceptibility to some classes of chemicals in the aging population.

A number of perfluorinated alkyl acids including perfluorooctanoic acid (PFOA) elicit effects similar to peroxisome proliferator chemicals (PPC) in mouse and rat liver. There is strong evidence that PPC cause many of their effects linked to liver cancer through the nuclear receptor peroxisome proliferator-activated receptor alpha (PPAR alpha). To determine the role of PPAR alpha in mediating PFOA transcriptional events, we compared the transcript profiles of the livers of wild-type or PPAR alpha-null mice exposed to PFOA or the PPAR alpha agonist WY-14,643 (WY). After 7 days of exposure, 85% or 99.7% of the genes altered by PFOA or WY exposure, respectively were dependent on PPAR alpha. The PPAR alpha-independent genes regulated by PFOA included those involved in lipid homeostasis and xenobiotic metabolism. Many of the lipid homeostasis genes including acyl-CoA oxidase (Acox1) were also regulated by WY in a PPAR alpha-dependent manner. The increased expression of these genes in PPAR alpha-null mice may be partly due to increases in PPAR gamma expression upon PFOA exposure. Many of the identified xenobiotic metabolism genes are known to be under control of the nuclear receptor CAR (constitutive activated/androstane receptor) and the transcription factor Nrf2 (nuclear factor erythroid 2-related factor 2). There was excellent correlation between the transcript profile of PPAR alpha-independent PFOA genes and those of activators of CAR including phenobarbital and 1,4-bis[2-(3,5-dichloropyridyloxy)] benzene (TCPOBOP) but not those regulated by the Nrf2 activator, dithiol-3-thione. These results indicate that PFOA alters most genes in wild-type mouse liver through PPAR alpha, but that a subset of genes are regulated by CAR and possibly PPAR gamma in the PPAR alpha-null mouse.

Data from toxicology and toxicogenomics studies are valuable, and can be combined for meta-analysis using public data repositories such as Chemical Effects in Biological Systems Knowledgebase, ArrayExpress, and Gene Expression Omnibus. In order to fully utilize the data for secondary analysis, it is necessary to have a description of the study and good annotation of the accompanying data. This study annotation permits sophisticated cross-study comparison and analysis, and allows data from comparable subjects to be identified and fully understood. The Minimal Information About a Microarray Experiment Standard was proposed to permit deposition and sharing of microarray data. We propose the first step toward an analogous standard for a toxicogenomics/toxicology study, by describing a checklist of information that best practices would suggest be included with the study data. When the information in this checklist is deposited together with the study data, the checklist information helps the public explore the study data in context of time, or identify data from similarly treated subjects, and also explore/identify potential sources of experimental variability. The proposed checklist summarizes useful information to include when sharing study data for publication, deposition into a database, or electronic exchange with collaborators. It is not a description of how to carry out an experiment, but a definition of how to describe an experiment. It is anticipated that once a toxicology checklist is accepted and put into use, then toxicology databases can be configured to require and output these fields, making it straightforward to annotate data for interpretation by others.

Toxicogenomic analysis of five environmental chemicals was performed to investigate the ability of genomics to predict toxicity, categorize chemicals, and elucidate mechanisms of toxicity. Three triazole antifungals (myclobutanil, propiconazole, and triadimefon) and two perfluorinated chemicals [perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS)] were administered daily via oral gavage for one, three, or five consecutive days to male Sprague-Dawley rats at single doses of 300, 300, 175, 20, or 10 mg/kg/day, respectively. Clinical chemistry, hematology, and histopathology were measured at all time points. Gene expression profiling of livers from three rats per treatment group at all time points was performed on the CodeLink Uniset Rat I Expression array. Data were analyzed in the context of a large reference toxicogenomic database containing gene expression profiles for over 630 chemicals. Genomic signatures predicting hepatomegaly and hepatic injury preceded those results for all five chemicals, and further analysis segregated chemicals into two distinct classes. The triazoles caused similar gene expression changes as other azole antifungals, particularly the induction of pregnane X receptor (PXR)-regulated xenobiotic metabolism and oxidative stress genes. In contrast, PFOA and PFOS exhibited peroxisome proliferator-activated receptor alpha agonist-like effects on genes associated with fatty acid homeostasis. PFOA and PFOS also resulted in downregulation of cholesterol biosynthesis genes, matching an in vivo decrease in serum cholesterol, and perturbation of thyroid hormone metabolism genes matched by serum thyroid hormone depletion in vivo. The concordance of in vivo observations and gene expression findings demonstrated the ability of genomics to accurately categorize chemicals, identify toxic mechanisms of action, and predict subsequent pathological responses.

Chronic exposure to peroxisome proliferators (PPs) leads to increased incidence of liver tumors in rodents. Liver tumor induction is thought to require increased hepatocyte proliferation and suppression of apoptosis. Transcript profiling showed increased expression of proapoptotic genes and decreased expression of antiapoptotic genes in the livers of mice exposed to the PP WY-14,643 (WY). We tested the hypothesis that prior exposure to WY would increase susceptibility to apoptosis inducers such as Jo2, an antibody which activates the Fas (Apo-1/CD95) death pathway. When compared to their untreated counterparts, wild-type mice pretreated with WY exhibited increased caspase-3 activation and hepatocyte apoptosis following challenge with Jo2. Livers from WY-treated peroxisome proliferator-activated receptor alpha (PPARalpha)-null mice were resistant to the effects of Jo2. In the absence of Jo2 and detectable apoptosis, wild-type mice treated with WY exhibited increases in the activated form of caspase-9. As caspase-9 is a component of the apoptosome, we examined the expression of upstream effectors of apoptosome activity including members of the Bcl-2 family. The levels of the antiapoptotic Mcl-1 transcript and protein were significantly decreased by PPs. PPARalpha-null mice were also resistant to another treatment (concanavalin A) that induces hepatocyte apoptosis. These results (1) indicate that PPARalpha activation increases sensitivity of the liver to apoptosis and (2) identify a mechanism by which PPARalpha could serve as a pharmacological target in diseases where apoptosis is a contributing feature.

Fumonisin mycotoxins occur worldwide in corn and corn-based foods. Fumonisin B1 (FB1) is a rodent liver carcinogen and suspected human carcinogen. It inhibits ceramide synthase and increases tissue sphinganine (Sa) and sphingosine (So) concentrations. Events linking disruption of sphingolipid metabolism and fumonisin toxicity are not fully understood; however, Sa and So were shown to bind mouse recombinant peroxisome proliferator-activated receptor alpha (PPARalpha) in vitro. To investigate the role of PPARalpha in fumonisin hepatotoxicity in vivo, wild-type (WT) and PPARalpha-null mice were fed control diets or diets containing 300 ppm FB1, Fusarium verticillioides culture material (CM) providing 300 ppm FB1, or 500 ppm of the peroxisome proliferator WY-14,643 (WY) for 1 week. WY-fed WT mice exhibited hepatomegaly, an effect not found in WY-fed PPARalpha-null mice, and WY did not change liver sphingoid base concentrations in either strain. Hepatotoxicity found in FB1- and CM-fed WT and PPARalpha-null mice was similar, qualitatively different from that found in WY-treated animals, and characterized by increased Sa concentration, apoptosis, and cell proliferation. Transcript profiling using oligonucleotide arrays showed that CM and FB1 elicited similar expression patterns of genes involved in cell proliferation, signal transduction, and glutathione metabolism that were different from that altered by WY. Real-time RT-PCR analysis of gene expression demonstrated PPARalpha-dependence of lipid metabolism gene expression in WY-treated mice, whereas PPARalpha-independent alterations of genes in lipid metabolism, and other categories, were found in CM- and FB1-fed mice. Together, these findings demonstrate that FB1- and CM-induced hepatotoxicity in mice does not require PPARalpha.

Many phthalate ester plasticizers are classified as peroxisome proliferators (PP), a large group of industrial and pharmaceutical chemicals. Like PP, exposure to some phthalates increases hepatocyte peroxisome and cellular proliferation, as well as the incidence of hepatocellular adenomas in mice and rats. Most effects of PP are mediated by three nuclear receptors called peroxisome proliferator-activated receptors (PPARalpha,beta,gamma). An obligate role for PPARalpha in PP-induced events leading to liver cancer is well-established. Exposure of rats in utero or in the neonate to a subset of phthalate esters causes profound, sometimes irreversible malformations in the male reproductive tract. We review here the data that supports or discounts roles for PPARs in phthalate-induced testis toxicity including (1) toxic effects of phthalates on the male reproductive tract, (2) expression of PPARs in the testis, (3) activation of PPARs by phthalates, (4) role of PPARalpha in testis toxicity, (5) gene targets of phthalates involved in steroid biosynthesis and catabolism, and (6) interactions between PPARs and other nuclear receptors that play roles in testis development and homeostasis. Critical research needs are identified that will help determine the significance of PPARs in phthalate-induced effects in the rat male reproductive tract and the relevance of toxicity to humans.

In this study, we hypothesized that many of the reported effects of phthalate esters and other peroxisome proliferators (PPs) in the testis are mediated by members of the PP- activated receptor (PPAR) family of transcription factors through alterations in proteins involved in steroidogenesis. Exposure of Leydig cells to PPs prevented cholesterol transport into the mitochondria after hormonal stimulation and inhibited steroid synthesis, without altering total cell protein synthesis or mitochondrial and DNA integrity. PPs also reduced the levels of the cholesterol-binding protein peripheral-type benzodiazepine receptor (PBR) because of a direct transcriptional inhibition of PBR gene expression in MA-10 Leydig cells. MA-10 cells contain mRNAs for PPARalpha and PPARbeta/delta, but not for PPARgamma. In vivo treatment of mice with PPs resulted in the reduction of both testis PBR mRNA and circulating testosterone levels, in agreement with the proposed role of PBR in steroidogenesis. By contrast, liver PBR mRNA levels were increased, in agreement with the proposed role of PBR in cell growth/tumor formation in nonsteroidogenic tissues. However, PPs did not inhibit testosterone production and testis PBR expression in PPARalpha-null mice. These results suggest that the antiandrogenic effect of PPs is mediated by a PPARalpha-dependent inhibition of Leydig cell PBR gene expression.

The purpose of this study was to investigate whether activation of the nuclear receptor PPARalpha is needed for protection from acetaminophen (APAP) hepatotoxicity produced by repeated administration of the peroxisome proliferator clofibrate (CFB). Female wild-type and PPARalpha-null mice received corn oil vehicle or 500 mg CFB/kg, ip, daily for 10 days. They were then fasted overnight (18 h) and either killed at 4 or 24 h after challenge with 400 mg APAP/kg. Controls received 50% propylene glycol vehicle only. In this model of CFB hepatoprotection, liver injury was assessed by measuring plasma sorbitol dehydrogenase activity and by histopathology at 24 h after APAP challenge. Significant hepatocellular necrosis was evident in both corn oil-pretreated PPARalpha-null and wild-type mice at 24 h after APAP challenge. In agreement with previous studies, CFB-pretreated wild-type mice showed marked protection against APAP toxicity. In contrast, CFB did not provide protection against APAP hepatotoxicity in the PPARalpha-null mice. Similarly, at 4 h after APAP challenge, hepatic glutathione depletion and selective arylation of cytosolic proteins were reduced significantly in CFB-pretreated wild-type mice, but not in PPARalpha-null mice. The lack of changes in APAP binding and NPSH depletion in CFB-pretreated, PPARalpha-null mice is consistent with the presence of significant liver injury at 24 h in this treatment group. These findings demonstrate that the protection against APAP hepatotoxicity by peroxisome proliferator treatment is mediated by the activation of PPARalpha.

The rat, mouse and human estrogen receptor (ER) exists as two subtypes, ER alpha and ER beta, which differ in the C-terminal ligand-binding domain and in the N-terminal transactivation domain. In this study, we investigated the estrogenic activity of environmental chemicals and phytoestrogens in competition binding assays with ER alpha or ER beta protein, and in a transient gene expression assay using cells in which an acute estrogenic response is created by cotransfecting cultures with recombinant human ER alpha or ER beta complementary DNA (cDNA) in the presence of an estrogen-dependent reporter plasmid. Saturation ligand-binding analysis of human ER alpha and ER beta protein revealed a single binding component for [3H]-17beta-estradiol (E2) with high affinity [dissociation constant (Kd) = 0.05 - 0.1 nM]. All environmental estrogenic chemicals [polychlorinated hydroxybiphenyls, dichlorodiphenyltrichloroethane (DDT) and derivatives, alkylphenols, bisphenol A, methoxychlor and chlordecone] compete with E2 for binding to both ER subtypes with a similar preference and degree. In most instances the relative binding affinities (RBA) are at least 1000-fold lower than that of E2. Some phytoestrogens such as coumestrol, genistein, apigenin, naringenin, and kaempferol compete stronger with E2 for binding to ER beta than to ER alpha. Estrogenic chemicals, as for instance nonylphenol, bisphenol A, o, p'-DDT and 2',4',6'-trichloro-4-biphenylol stimulate the transcriptional activity of ER alpha and ER beta at concentrations of 100-1000 nM. Phytoestrogens, including genistein, coumestrol and zearalenone stimulate the transcriptional activity of both ER subtypes at concentrations of 1-10 nM. The ranking of the estrogenic potency of phytoestrogens for both ER subtypes in the transactivation assay is different; that is, E2 >> zearalenone = coumestrol > genistein > daidzein > apigenin = phloretin > biochanin A = kaempferol = naringenin > formononetin = ipriflavone = quercetin = chrysin for ER alpha and E2 >> genistein = coumestrol > zearalenone > daidzein > biochanin A = apigenin = kaempferol = naringenin > phloretin = quercetin = ipriflavone = formononetin = chrysin for ER beta. Antiestrogenic activity of the phytoestrogens could not be detected, except for zearalenone which is a full agonist for ER alpha and a mixed agonist-antagonist for ER beta. In summary, while the estrogenic potency of industrial-derived estrogenic chemicals is very limited, the estrogenic potency of phytoestrogens is significant, especially for ER beta, and they may trigger many of the biological responses that are evoked by the physiological estrogens.