Search

Identification of chemicals that affect hormone-regulated systems will help to predict endocrine disruption. In our previous study, a 46 gene biomarker was found to be an accurate predictor of estrogen receptor (ER) α modulation in chemically treated MCF-7 cells. Here, potential ERα modulators were identified using the biomarker by screening a microarray compendium consisting of ∼1600 gene expression comparisons representing exposure to ∼1200 chemicals. A total of ∼170 chemicals were identified as potential ERα modulators. In the Connectivity Map 2.0 collection, 75 and 39 chemicals were predicted to activate or suppress ERα, and they included 12 and six known ERα agonists and antagonists/selective ERα modulators, respectively. Nineteen and eight of the total number were also identified as active in an ERα transactivation assay carried out in an MCF-7-derived cell line used to screen the Tox21 10K chemical library in agonist or antagonist modes, respectively. Chemicals predicted to modulate ERα in MCF-7 cells were examined further using global and targeted gene expression in wild-type and ERα-null cells, transactivation assays, and cell-free ERα coregulator interaction assays. Environmental chemicals classified as weak and very weak agonists were confirmed to activate ERα including apigenin, kaempferol, and oxybenzone. Novel activators included digoxin, nabumetone, ivermectin, and six progestins. Novel suppressors included emetine, mifepristone, niclosamide, and proscillaridin. Our strategy will be useful to identify environmentally relevant ERα modulators in future high-throughput transcriptomic screens.

High-throughput transcriptomics (HTTr) has the potential to support efforts to reduce or replace some animal tests. In past studies, we described a computational approach utilizing a gene expression biomarker consisting of 46 genes to predict estrogen receptor (ER) activity after chemical exposure in ER-positive human breast cancer cells including the MCF-7 cell line. We hypothesized that the biomarker model could identify ER activities of chemicals examined by Endocrine Disruptor Screening Program (EDSP) Tier 1 screening assays in which transcript profiles of the same chemicals were examined in MCF-7 cells. For the 62 chemicals examined including 5 chemicals examined in this study using RNA-Seq, the ER biomarker model accuracy was 1) 97% for in vitro reference chemicals, 2) 76-85% for guideline uterotrophic assays, and 3) 87-88% for guideline and nonguideline uterotrophic assays. For the same chemicals, these accuracies were similar or slightly better than those of the ToxCast ER model based on 18 in vitro assays. The performance of the ER biomarker model indicates that HTTr interpreted using the ER biomarker correctly identifies active and inactive ER reference chemicals. As part of the HTTr screening program the approach could rapidly identify chemicals with potential ER bioactivities for additional screening and testing.

Rodent cancer bioassays have been long-required studies for regulatory assessment of human cancer hazard and risk. These studies use hundreds of animals, are resource intensive, and certain aspects of these studies have limited human relevance. The past 10 years have seen an exponential growth of new technologies with the potential to effectively evaluate human cancer hazard and risk while reducing, refining, or replacing animal use. To streamline and facilitate uptake of new technologies, a workgroup comprised of scientists from government, academia, non-governmental organizations, and industry stakeholders developed a framework for waiver rationales of rodent cancer bioassays for consideration in agrochemical safety assessment. The workgroup used an iterative approach, incorporating regulatory agency feedback, and identifying critical information to be considered in a risk assessment-based weight of evidence determination of the need for rodent cancer bioassays. The reporting framework described herein was developed to support a chronic toxicity and carcinogenicity study waiver rationale, which includes information on use pattern(s), exposure scenario(s), pesticidal mode-of-action, physicochemical properties, metabolism, toxicokinetics, toxicological data including mechanistic data, and chemical read-across from similar registered pesticides. The framework could also be applied to endpoints other than chronic toxicity and carcinogenicity, and for chemicals other than agrochemicals.

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.

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.

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.

Endocrine disrupting chemicals raise a lot of interest and concern regarding their risk for human health and the environment. They represent a broad variety of natural and synthetic chemicals with different levels of endocrine activity evaluation. In particular, for high production volume chemicals, new methods are required to enable the evaluation of the vast number of chemicals for their potential to alter the endocrine system and prioritize them for deeper characterization. The ToxCast program from the US EPA provides data from high throughput screening assays to develop computational tools aimed at rapid in vitro bioactivity screening and prioritization. Using publicly available data (ToxCast and ToxRef databases), we evaluate whether in vitro assay evaluations could predict in vivo outcomes observed in rat long-term studies for more than 400 chemicals. We focus on effects observed in three endocrine and two sex accessory organs and 42 in vitro assays related to pathways associated with endocrine related toxicity. First, using simple statistical correlation we demonstrate that there is no mutual linear correlation between the selected in vitro assays and any in vivo outcome, with balanced accuracies around 50% for each assay-outcome pair. Then, by applying machine learning to investigate potential non-linear correlations, we show that the combination of different in vitro assays is not correlated with the long-term in vivo effects and cannot help to predict them since balanced accuracies are also around 50%. Moreover, the prediction based on in vitro assays is not better than the one based on classical QSAR methods. This study highlights that the selected in vitro assays do not provide information about in vivo outcomes observed in endocrine and associated organs in long-term rat in vivo studies and stresses the need for the development of in vitro assays that reflect the compounds’ pharmacokinetic properties.