Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation.

Gene transcription profiles across tissues are largely defined by the activity of regulatory elements, most of which correspond to regions of accessible chromatin. Regulatory element activity is in turn modulated by genetic variation, resulting in variable transcription rates across individuals. The interplay of these factors, however, is poorly understood. Here we characterize expression and chromatin state dynamics across three tissues-liver, lung, and kidney-in 47 strains of the Collaborative Cross (CC) mouse population, examining the regulation of these dynamics by expression quantitative trait loci (eQTL) and chromatin QTL (cQTL). QTL whose allelic effects were consistent across tissues were detected for 1,101 genes and 133 chromatin regions. Also detected were eQTL and cQTL whose allelic effects differed across tissues, including local-eQTL for Pik3c2g detected in all three tissues but with distinct allelic effects. Leveraging overlapping measurements of gene expression and chromatin accessibility on the same mice from multiple tissues, we used mediation analysis to identify chromatin and gene expression intermediates of eQTL effects. Based on QTL and mediation analyses over multiple tissues, we propose a causal model for the distal genetic regulation of Akr1e1, a gene involved in glycogen metabolism, through the zinc finger transcription factor Zfp985 and chromatin intermediates. This analysis demonstrates the complexity of transcriptional and chromatin dynamics and their regulation over multiple tissues, as well as the value of the CC and related genetic resource populations for identifying specific regulatory mechanisms within cells and tissues.

A systematic approach to evaluate plausible modes of actions for mouse lung tumors in mice exposed to 4-methylimidozole.

The National Toxicology Program (NTP) reported that chronic dietary exposure to 4-methylimidazole (4-MeI) increased the incidence of lung adenomas/carcinomas beyond the normally high spontaneous rate in B6C3F1 mice. To examine plausible modes of action (MoAs) for mouse lung tumors (MLTs) upon exposure to high levels of 4-MeI, and their relevance in assessing human risk, a systematic approach was used to identify and evaluate mechanistic data (in vitro and in vivo) in the primary and secondary literature, along with high-throughput screening assay data. Study quality, relevance, and activity of mechanistic data identified across the evidence-base were organized according to key characteristics of carcinogens (KCCs) to identify potential key events in known or novel MLT MoAs. Integration of these evidence streams provided confirmation that 4-MeI lacks genotoxic and cytotoxic activity with some evidence to support a lack of mitogenic activity. Further evaluation of contextual and chemical-specific characteristics of 4-MeI was consequently undertaken. Due to lack of genotoxicity, along with transcriptomic and histopathological lung changes up to 28 and 90 days of exposure, the collective evidence suggests MLTs observed following exposure to high levels of 4-MeI develop at a late stage in the mouse chronic bioassay, albeit the exact MoA remains unclear.

Assessment of the biochemical pathways for acetaminophen toxicity: Implications for its carcinogenic hazard potential.

In 2019 California's Office of Environmental Health Hazard Assessment (OEHHA) initiated a review of the carcinogenic hazard potential of acetaminophen. In parallel with this review, herein we evaluated the mechanistic data related to the steps and timing of cellular events following therapeutic recommended (≤4 g/day) and higher doses of acetaminophen that may cause hepatotoxicity to evaluate whether these changes indicate that acetaminophen is a carcinogenic hazard. At therapeutic recommended doses, acetaminophen forms limited amounts of N-acetyl-p-benzoquinone-imine (NAPQI) without adverse cellular effects. Following overdoses of acetaminophen, there is potential for more extensive formation of NAPQI and depletion of glutathione, which may result in mitochondrial dysfunction and DNA damage, but only at doses that result in cell death - thus making it implausible for acetaminophen to induce the kind of stable, genetic damage in the nucleus indicative of a genotoxic or carcinogenic hazard in humans. The collective data demonstrate a lack of a plausible mechanism related to carcinogenicity and are consistent with rodent cancer bioassays, epidemiological results reviewed in companion manuscripts in this issue, as well as conclusions of multiple international health authorities.

Assessment of the Mode of Action Underlying the Effects of GenX in Mouse Liver and Implications for Assessing Human Health Risks.

GenX is an alternative to environmentally persistent long-chain perfluoroalkyl and polyfluoroalkyl substances. Mice exposed to GenX exhibit liver hypertrophy, elevated peroxisomal enzyme activity, and other apical endpoints consistent with peroxisome proliferators. To investigate the potential role of peroxisome proliferator-activated receptor alpha (PPARα) activation in mice, and other molecular signals potentially related to observed liver changes, RNA sequencing was conducted on paraffin-embedded liver sections from a 90-day subchronic toxicity study of GenX conducted in mice. Differentially expressed genes were identified for each treatment group, and gene set enrichment analysis was conducted using gene sets that represent biological processes and known canonical pathways. Peroxisome signaling and fatty acid metabolism were among the most significantly enriched gene sets in both sexes at 0.5 and 5 mg/kg GenX; no pathways were enriched at 0.1 mg/kg. Gene sets specific to the PPARα subtype were significantly enriched. These findings were phenotypically anchored to histopathological changes in the same tissue blocks: hypertrophy, mitoses, and apoptosis. In vitro PPARα transactivation assays indicated that GenX activates mouse PPARα. These results indicate that the liver changes observed in GenX-treated mice occur via a mode of action (MOA) involving PPARα, an important finding for human health risk assessment as this MOA has limited relevance to humans.

Population-Based Analysis of DNA Damage and Epigenetic Effects of 1,3-Butadiene in the Mouse.

Metabolism of 1,3-butadiene, a known human and rodent carcinogen, results in formation of reactive epoxides, a key event in its carcinogenicity. Although mice exposed to 1,3-butadiene present DNA adducts in all tested tissues, carcinogenicity is limited to liver, lung, and lymphoid tissues. Previous studies demonstrated that strain- and tissue-specific epigenetic effects in response to 1,3-butadiene exposure may influence susceptibly to DNA damage and serve as a potential mechanism of tissue-specific carcinogenicity. This study aimed to investigate interindividual variability in the effects of 1,3-butadiene using a population-based mouse model. Male mice from 20 Collaborative Cross strains were exposed to 0 or 635 ppm 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. We evaluated DNA damage and epigenetic effects in target (lung and liver) and nontarget (kidney) tissues of 1,3-butadiene-induced carcinogenesis. DNA damage was assessed by measuring N-7-(2,3,4-trihydroxybut-1-yl)-guanine (THB-Gua) adducts. To investigate global histone modification alterations, we evaluated the trimethylation and acetylation of histones H3 and H4 across tissues. Changes in global cytosine DNA methylation were evaluated from the levels of methylation of LINE-1 and SINE B1 retrotransposons. We quantified the degree of variation across strains, deriving a chemical-specific human variability factor to address population variability in carcinogenic risk, which is largely ignored in current cancer risk assessment practice. Quantitative trait locus mapping identified four candidate genes related to chromatin remodeling whose variation was associated with interstrain susceptibility. Overall, this study uses 1,3-butadiene to demonstrate how the Collaborative Cross mouse population can be used to identify the mechanisms for and quantify the degree of interindividual variability in tissue-specific effects that are relevant to chemically induced carcinogenesis.

Comparison of Gene Expression Responses in the Small Intestine of Mice Following Exposure to 3 Carcinogens Using the S1500+ Gene Set Informs a Potential Common Adverse Outcome Pathway.

Carcinogenesis of the small intestine is rare in humans and rodents. Oral exposure to hexavalent chromium (Cr(VI)) and the fungicides captan and folpet induce intestinal carcinogenesis in mice. Previously (Toxicol Pathol. 330:48-52), we showed that B6C3F1 mice exposed to carcinogenic concentrations of Cr(VI), captan, or folpet for 28 days exhibited similar histopathological responses including villus enterocyte cytotoxicity and regenerative crypt epithelial hyperplasia. Herein, we analyze transcriptomic responses from formalin-fixed, paraffin-embedded duodenal sections from the aforementioned study. TempO-Seq technology and the S1500+ gene set were used to analyze transcription responses. Transcriptional responses were similar between all 3 agents; gene-level comparison identified 126/546 (23%) differentially expressed genes altered in the same direction, with a total of 25 upregulated pathways. These changes were related to cellular metabolism, stress, inflammatory/immune cell response, and cell proliferation, including upregulation in hypoxia inducible factor 1 (HIF-1) and activator protein 1 (AP1) signaling pathways, which have also been shown to be related to intestinal injury and angiogenesis/carcinogenesis. The similar molecular-, cellular-, and tissue-level changes induced by these 3 carcinogens can be informative for the development of an adverse outcome pathway for intestinal cancer.

Sex-specific differences in genotoxic and epigenetic effects of 1,3-butadiene among mouse tissues.

Exposure to environmental chemicals has been shown to have an impact on the epigenome. One example is a known human carcinogen 1,3-butadiene which acts primarily by a genotoxic mechanism, but also disrupts the chromatin structure by altering patterns of cytosine DNA methylation and histone modifications. Sex-specific differences in 1,3-butadiene-induced genotoxicity and carcinogenicity are well established; however, it remains unknown whether 1,3-butadiene-associated epigenetic alterations are also sex dependent. Therefore, we tested the hypothesis that inhalational exposure to 1,3-butadiene will result in sex-specific epigenetic alterations. DNA damage and epigenetic effects of 1,3-butadiene were evaluated in liver, lung, and kidney tissues of male and female mice of two inbred strains (C57BL/6J and CAST/EiJ). Mice were exposed to 0 or 425 ppm of 1,3-butadiene by inhalation (6 h/day, 5 days/week) for 2 weeks. Strain- and tissue-specific differences in 1,3-butadiene-induced DNA adducts and crosslinks were detected in the liver, lung and kidney; however, significant sex-specific differences in DNA damage were observed in the lung of C57BL/6J mice only. In addition, we assessed expression of the DNA repair genes and observed a marked upregulation of Mgmt in the kidney in female C57BL/6J mice. Sex-specific epigenetic effects of 1,3-butadiene exposure were evident in alterations of cytosine DNA methylation and histone modifications in the liver and lung in both strains. Specifically, we observed a loss of cytosine DNA methylation in the liver and lung of male and female 1,3-butadiene-exposed C57BL/6J mice, whereas hypermethylation was found in the liver and lung in 1,3-butadiene-exposed female CAST/EiJ mice. Our findings suggest that strain- and sex-specific effects of 1,3-butadiene on the epigenome may contribute to the known differences in cancer susceptibility.

Tissue- and strain-specific effects of a genotoxic carcinogen 1,3-butadiene on chromatin and transcription.

Epigenetic effects of environmental chemicals are under intense investigation to fill existing knowledge gaps between environmental/occupational exposures and adverse health outcomes. Chromatin accessibility is one prominent mechanism of epigenetic control of transcription, and understanding of the chemical effects on both could inform the causal role of epigenetic alterations in disease mechanisms. In this study, we hypothesized that baseline variability in chromatin organization and transcription profiles among various tissues and mouse strains influence the outcome of exposure to the DNA damaging chemical 1,3-butadiene. To test this hypothesis, we evaluated DNA damage along with comprehensive quantification of RNA transcripts (RNA-seq), identification of accessible chromatin (ATAC-seq), and characterization of regions with histone modifications associated with active transcription (ChIP-seq for acetylation at histone 3 lysine 27, H3K27ac). We collected these data in the lung, liver, and kidney of mice from two genetically divergent strains, C57BL/6J and CAST/EiJ, that were exposed to clean air or to 1,3-butadiene (~600 ppm) for 2 weeks. We found that tissue effects dominate differences in both gene expression and chromatin states, followed by strain effects. At baseline, xenobiotic metabolism was consistently more active in CAST/EiJ, while immune system pathways were more active in C57BL/6J across tissues. Surprisingly, even though all three tissues in both strains harbored butadiene-induced DNA damage, little transcriptional effect of butadiene was observed in liver and kidney. Toxicologically relevant effects of butadiene in the lung were on the pathways of xenobiotic metabolism and inflammation. We also found that variability in chromatin accessibility across individuals (i.e., strains) only partially explains the variability in transcription. This study showed that variation in the basal states of epigenome and transcriptome may be useful indicators for individuals or tissues susceptible to genotoxic environmental chemicals.

An Adverse Outcome Pathway (AOP) for forestomach tumors induced by non-genotoxic initiating events.

The utility of rodent forestomach tumor data for hazard and risk assessment has been examined for decades because humans do not have a forestomach, and these tumors occur by varying modes of action (MOAs). We have used the MOA for ethyl acrylate (EA) to develop an Adverse Outcome Pathway (AOP) for forestomach tumors caused by non-genotoxic initiating events. These tumors occur secondary to site of contact induced epithelial cytotoxicity and regenerative repair-driven proliferation. For EA, the critical initiating event (IE) is epithelial cytotoxicity, and supporting key events (KEs) at the cellular and tissue level are increased cell proliferation (KE1) resulting in sustained hyperplasia (KE2), with the adverse outcome of forestomach papillomas and carcinomas. For EA, a pre-molecular initiating event (pre-MIE) of sustained glutathione depletion is probable. Supporting data from butylated hydroxyanisole (BHA) are also reviewed. Although there may be some variability in the pre-MIEs and IEs for BHA and EA, they share the same KEs, and evidence for BHA confers support for the AOP. Evolved Bradford Hill considerations of biological plausibility, essentiality, and empirical support were evaluated per OECD guidance. Although an MIE is not specifically described, overall confidence in the AOP is high due to well-developed and accepted evidence streams, and the AOP can be used for regulatory applications including hazard identification and risk assessment for chemicals that act by this AOP.

Assessment of the mode of action underlying development of forestomach tumors in rodents following oral exposure to ethyl acrylate and relevance to humans.

Chronic repeated gavage dosing of high concentrations of ethyl acrylate (EA) causes forestomach tumors in rats and mice. For two decades, there has been general consensus that these tumors are unique to rodents because of: i) lack of carcinogenicity in other organs, ii) specificity to the forestomach (an organ unique to rodents which humans do not possess), iii) lack of carcinogenicity by other routes of exposure, and iv) obvious site of contact toxicity at carcinogenic doses. In 1986, EA was classified as possibly carcinogenic to humans by the International Agency for Research on Cancer (IARC). However, by applying a MOA analyses and human relevance framework assessment, the weight-of-evidence supports a cytotoxic MOA with the following key events: i) bolus delivery of EA to forestomach lumen and subsequent absorption, ii) cytotoxicity likely due to saturation of enzymatic detoxification, iii) chronic regenerative hyperplasia, and iv) spontaneous mutation due to increased cell replication and cell population. Clonal expansion of initiated cells thus results in late onset tumorigenesis. The key events in this 'wound and healing' MOA provide high confidence in the MOA as assessed by evolved Bradford-Hill Criteria. The weight-of-evidence supported by the proposed MOA, combined with a unique tissue that does not exist in humans, indicates that EA is highly unlikely to pose a human cancer hazard.

Benchmark Dose Modeling Estimates of the Concentrations of Inorganic Arsenic That Induce Changes to the Neonatal Transcriptome, Proteome, and Epigenome in a Pregnancy Cohort.

Prenatal inorganic arsenic (iAs) exposure influences the expression of critical genes and proteins associated with adverse outcomes in newborns, in part through epigenetic mediators. The doses at which these genomic and epigenomic changes occur have yet to be evaluated in the context of dose-response modeling. The goal of the present study was to estimate iAs doses that correspond to changes in transcriptomic, proteomic, epigenomic, and integrated multi-omic signatures in human cord blood through benchmark dose (BMD) modeling. Genome-wide DNA methylation, microRNA expression, mRNA expression, and protein expression levels in cord blood were modeled against total urinary arsenic (U-tAs) levels from pregnant women exposed to varying levels of iAs. Dose-response relationships were modeled in BMDExpress, and BMDs representing 10% response levels were estimated. Overall, DNA methylation changes were estimated to occur at lower exposure concentrations in comparison to other molecular endpoints. Multi-omic module eigengenes were derived through weighted gene co-expression network analysis, representing co-modulated signatures across transcriptomic, proteomic, and epigenomic profiles. One module eigengene was associated with decreased gestational age occurring alongside increased iAs exposure. Genes/proteins within this module eigengene showed enrichment for organismal development, including potassium voltage-gated channel subfamily Q member 1 (KCNQ1), an imprinted gene showing differential methylation and expression in response to iAs. Modeling of this prioritized multi-omic module eigengene resulted in a BMD(BMDL) of 58(45) µg/L U-tAs, which was estimated to correspond to drinking water arsenic concentrations of 51(40) µg/L. Results are in line with epidemiological evidence supporting effects of prenatal iAs occurring at levels <100 µg As/L urine. Together, findings present a variety of BMD measures to estimate doses at which prenatal iAs exposure influences neonatal outcome-relevant transcriptomic, proteomic, and epigenomic profiles.

Differentially expressed MicroRNAs provide mechanistic insight into fibrosis-associated liver carcinogenesis in mice.

Hepatocellular carcinoma (HCC) is one of the most prevalent human cancers, with a rising incidence worldwide. The molecular mechanisms associated with the development of HCC are complex and include multiple interconnected molecular alterations with mounting evidence indicating an important role of microRNAs (miRNAs) in the pathogenesis of HCC. In humans, the development of HCC is commonly associated with liver cirrhosis. To study fibrosis-associated liver carcinogenesis, we used a mouse model designed to emulate the development of HCC in cirrhotic liver. Specifically, we were interested in evaluating the role of miRNAs in the molecular pathogenesis of liver carcinogenesis in male B6C3F1/J mice treated with N-nitrosodiethylamine (DEN) or carbon tetrachloride (CCl4 ) alone or a combination of DEN and CCl4 and characterized by a differential tumor incidence that increased in the following order: DEN

A mouse model of alcoholic liver fibrosis-associated acute kidney injury identifies key molecular pathways.

Clinical data strongly indicate that acute kidney injury (AKI) is a critical complication in alcoholic hepatitis, an acute-on-chronic form of liver failure in patients with advanced alcoholic fibrosis. Development of targeted therapies for AKI in this setting is hampered by the lack of an animal model. To enable research into molecular drivers and novel therapies for fibrosis- and alcohol-associated AKI, we aimed to combine carbon tetrachloride (CCl4)-induced fibrosis with chronic intra-gastric alcohol feeding. Male C57BL/6J mice were administered a low dose of CCl4 (0.2ml/kg 2× week/6weeks) followed by alcohol intragastrically (up to 25g/kg/day for 3weeks) and with continued CCl4. We observed that combined treatment with CCl4 and alcohol resulted in severe liver injury, more pronounced than using each treatment alone. Importantly, severe kidney injury was evident only in the combined treatment group. This mouse model reproduced distinct pathological features consistent with AKI in human alcoholic hepatitis. Transcriptomic analysis of kidneys revealed profound effects in the combined treatment group, with enrichment for damage-associated pathways, such as apoptosis, inflammation, immune-response and hypoxia. Interestingly, Havcr1 and Lcn2, biomarkers of AKI, were markedly up-regulated. Overall, this study established a novel mouse model of fibrosis- and alcohol-associated AKI and identified key mechanistic pathways.

Genetic and epigenetic changes in fibrosis-associated hepatocarcinogenesis in mice.

Hepatocellular carcinoma (HCC) is one of the most prevalent cancers and is rising in incidence worldwide. The molecular mechanisms leading to the development of HCC are complex and include both genetic and epigenetic events. To determine the relative contribution of these alterations in liver tumorigenesis, we evaluated epigenetic modifications at both global and gene specific levels, as well as the mutational profile of genes commonly altered in liver tumors. A mouse model of fibrosis-associated liver cancer that was designed to emulate cirrhotic liver, a prevailing disease state observed in most humans with HCC, was used. Tumor and nontumor liver samples from B6C3F1 mice treated with N-nitrosodiethylamine (DEN; a single ip injection of 1 mg/kg at 14 days of age) and carbon tetrachloride (CCl4; 0.2 ml/kg, 2 times/week ip starting at 8 weeks of age for 14 weeks), as well as corresponding vehicle control animals, were analyzed for genetic and epigenetic alterations. H-ras, Ctnnb1 and Hnf1α genes were not mutated in tumors in mice treated with DEN+CCl4 . In contrast, the increased tumor incidence in mice treated with DEN+CCl4 was associated with marked epigenetic changes in liver tumors and nontumor liver tissue, including demethylation of genomic DNA and repetitive elements, a decrease in histone 3 lysine 9 trimethylation (H3K9me3) and promoter hypermethylation and functional downregulation of Riz1, a histone lysine methyltransferase tumor suppressor gene. Additionally, the reduction in H3K9me3 was accompanied by increased expression of long interspersed nucleotide elements 1 and short interspersed nucleotide elements B2, which is an indication of genomic instability. In summary, our results suggest that epigenetic events, rather than mutations in known cancer-related genes, play a prominent role in increased incidence of liver tumors in this mouse model of fibrosis-associated liver cancer.

Enabling novel paradigms: a biological questions-based approach to human chemical hazard and drug safety assessment.

Throughput needs, costs of time and resources, and concerns about the use of animals in hazard and safety assessment studies are fueling a growing interest in adopting new approach methodologies for use in product development and risk assessment. However, current efforts to define "next-generation risk assessment" vary considerably across commercial and regulatory sectors, and an a priori definition of the biological scope of data needed to assess hazards is generally lacking. We propose that the absence of clearly defined questions that can be answered during hazard assessment is the primary barrier to the generation of a paradigm flexible enough to be used across varying product development and approval decision contexts. Herein, we propose a biological questions-based approach (BQBA) for hazard and safety assessment to facilitate fit-for-purpose method selection and more efficient evidence-based decision-making. The key pillars of this novel approach are bioavailability, bioactivity, adversity, and susceptibility. This BQBA is compared with current hazard approaches and is applied in scenarios of varying pathobiological understanding and/or regulatory testing requirements. To further define the paradigm and key questions that allow better prediction and characterization of human health hazard, a multidisciplinary collaboration among stakeholder groups should be initiated.

Chemistry domain of applicability evaluation against existing estrogen receptor high-throughput assay-based activity models.

Introduction: The U. S. Environmental Protection Agency's Endocrine Disruptor Screening Program (EDSP) Tier 1 assays are used to screen for potential endocrine system-disrupting chemicals. A model integrating data from 16 high-throughput screening assays to predict estrogen receptor (ER) agonism has been proposed as an alternative to some low-throughput Tier 1 assays. Later work demonstrated that as few as four assays could replicate the ER agonism predictions from the full model with 98% sensitivity and 92% specificity. The current study utilized chemical clustering to illustrate the coverage of the EDSP Universe of Chemicals (UoC) tested in the existing ER pathway models and to investigate the utility of chemical clustering to evaluate the screening approach using an existing 4-assay model as a test case. Although the full original assay battery is no longer available, the demonstrated contribution of chemical clustering is broadly applicable to assay sets, chemical inventories, and models, and the data analysis used can also be applied to future evaluation of minimal assay models for consideration in screening. Methods: Chemical structures were collected for 6,947 substances via the CompTox Chemicals Dashboard from the over 10,000 UoC and grouped based on structural similarity, generating 826 chemical clusters. Of the 1,812 substances run in the original ER model, 1,730 substances had a single, clearly defined structure. The ER model chemicals with a clearly defined structure that were not present in the EDSP UoC were assigned to chemical clusters using a k-nearest neighbors approach, resulting in 557 EDSP UoC clusters containing at least one ER model chemical. Results and Discussion: Performance of an existing 4-assay model in comparison with the existing full ER agonist model was analyzed as related to chemical clustering. This was a case study, and a similar analysis can be performed with any subset model in which the same chemicals (or subset of chemicals) are screened. Of the 365 clusters containing >1 ER model chemical, 321 did not have any chemicals predicted to be agonists by the full ER agonist model. The best 4-assay subset ER agonist model disagreed with the full ER agonist model by predicting agonist activity for 122 chemicals from 91 of the 321 clusters. There were 44 clusters with at least two chemicals and at least one agonist based upon the full ER agonist model, which allowed accuracy predictions on a per-cluster basis. The accuracy of the best 4-assay subset ER agonist model ranged from 50% to 100% across these 44 clusters, with 32 clusters having accuracy ≥90%. Overall, the best 4-assay subset ER agonist model resulted in 122 false-positive and only 2 false-negative predictions compared with the full ER agonist model. Most false positives (89) were active in only two of the four assays, whereas all but 11 true positive chemicals were active in at least three assays. False positive chemicals also tended to have lower area under the curve (AUC) values, with 110 out of 122 false positives having an AUC value below 0.214, which is lower than 75% of the positives as predicted by the full ER agonist model. Many false positives demonstrated borderline activity. The median AUC value for the 122 false positives from the best 4-assay subset ER agonist model was 0.138, whereas the threshold for an active prediction is 0.1. Conclusion: Our results show that the existing 4-assay model performs well across a range of structurally diverse chemicals. Although this is a descriptive analysis of previous results, several concepts can be applied to any screening model used in the future. First, the clustering of the chemicals provides a means of ensuring that future screening evaluations consider the broad chemical space represented by the EDSP UoC. The clusters can also assist in prioritizing future chemicals for screening in specific clusters based on the activity of known chemicals in those clusters. The clustering approach can be useful in providing a framework to evaluate which portions of the EDSP UoC chemical space are reliably covered by in silico and in vitro approaches and where predictions from either method alone or both methods combined are most reliable. The lessons learned from this case study can be easily applied to future evaluations of model applicability and screening to evaluate future datasets.

Optimizing androgen receptor prioritization using high-throughput assay-based activity models.

Introduction: Computational models using data from high-throughput screening assays have promise for prioritizing and screening chemicals for testing under the U.S. Environmental Protection Agency's Endocrine Disruptor Screening Program (EDSP). The purpose of this work was to demonstrate a data processing method for the determination of optimal minimal assay batteries from a larger comprehensive model, to provide a uniform method of evaluating the performance of future minimal assay batteries compared with the androgen receptor (AR) pathway model, and to incorporate chemical cluster analysis into this evaluation. Although several of the assays in the AR pathway model are no longer available through the original vendor, this approach could be used for future evaluations of minimal assay models for prioritization and screening. Methods: We compared two previously published models and found that an expanded 14-assay model had higher sensitivity for antagonists, whereas the original 11-assay model had slightly higher sensitivity for agonists. We then investigated subsets of assays in the original AR pathway model to optimize overall testing strategies that minimize cost while maintaining sensitivity across a broad chemical space. Results and Discussion: Evaluation of the critical assays across subset models derived from the 14-assay model identified three critical assays for predicting antagonism and two critical assays for predicting agonism. A minimum of nine assays is required for predicting agonism and antagonism with high sensitivity (95%). However, testing workflows guided by chemical structure-based clusters can reduce the average number of assays needed per chemical by basing the assays selected for testing on the likelihood of a chemical being an AR agonist, according to its structure. Our results show that a multi-stage testing workflow can provide 95% sensitivity while requiring only 48% of the resources required for running all assays from the original full models. The resources can be reduced further by incorporating in silico activity predictions. Conclusion: This work illustrates a data-driven approach that incorporates chemical clustering and simultaneous consideration of antagonism and agonism mechanisms to more efficiently screen chemicals. This case study provides a proof of concept for prioritization and screening strategies that can be utilized in future analyses to minimize the overall number of assays needed for predicting AR activity, which will maximize the number of chemicals that can be tested and allow data-driven prioritization of chemicals for further screening under the EDSP.

Updated systematic assessment of human, animal and mechanistic evidence demonstrates lack of human carcinogenicity with consumption of aspartame.

Aspartame has been studied extensively and evaluated for its safety in foods and beverages yet concerns for its potential carcinogenicity have persisted, driven primarily by animal studies conducted at the Ramazzini Institute (RI). To address this controversy, an updated systematic review of available human, animal, and mechanistic data was conducted leveraging critical assessment tools to consider the quality and reliability of data. The evidence base includes 12 animal studies and >40 epidemiological studies reviewed by the World Health Organization which collectively demonstrate a lack of carcinogenic effect. Assessment of >1360 mechanistic endpoints, including many guideline-based genotoxicity studies, demonstrate a lack of activity associated with endpoints grouped to key characteristics of carcinogens. Other non-specific mechanistic data (e.g., mixed findings of oxidative stress across study models, tissues, and species) do not provide evidence of a biologically plausible carcinogenic pathway associated with aspartame. Taken together, available evidence supports that aspartame consumption is not carcinogenic in humans and that the inconsistent findings of the RI studies may be explained by flaws in study design and conduct (despite additional analyses to address study limitations), as acknowledged by authoritative bodies.

Crypt and Villus Transcriptomic Responses in Mouse Small Intestine Following Oral Exposure to Hexavalent Chromium.

Oral exposure to hexavalent chromium (Cr(VI)) induces tumors in the mouse duodenum. Previous microarray-based transcriptomic analyses of homogenized mouse duodenal tissue have demonstrated Cr(VI)-induced alterations in various cellular pathways and processes. However, X-ray fluorescence microscopy indicates that chromium localizes primarily to the duodenal villi following exposure to Cr(VI), suggesting that previous transcriptomic analyses of homogenized tissue provide an incomplete picture of transcriptomic responses in the duodenum. Herein, transcriptomic analyses were conducted separately on crypt and villus tissue from formalin-fixed paraffin-embedded transverse duodenal sections from the same study in which microarray-based analyses were previously conducted. A total of 28 groups (7 doses × 2 timepoints × 2 tissue compartments) were analyzed for differential gene expression, dose-response, and gene set enrichment. Tissue compartment isolation was confirmed by differences in expression of typical markers of crypt and villus compartments. Fewer than 21 genes were altered in the crypt compartment of mice exposed to 0.1-5 ppm Cr(VI) for 7 or 90 days, which increased to hundreds or thousands of genes at ≥20 ppm Cr(VI). Consistent with histological evidence for crypt proliferation, a significant, dose-dependent increase in genes that regulate mitotic cell cycle was prominent in the crypt, while subtle in the villus, when compared with samples from time-matched controls. Minimal transcriptomic evidence of DNA damage response in either the crypts or the villi is consistent with published in vivo genotoxicity data. These results are also discussed in the context of modes of action that have been proposed for Cr(VI)-induced small intestine tumors in mice.

Epigenetic alterations induced by genotoxic occupational and environmental human chemical carcinogens: An update of a systematic literature review.

Epigenetic alterations, such as changes in DNA methylation, histones/chromatin structure, nucleosome positioning, and expression of non-coding RNAs, are recognized among key characteristics of carcinogens; they may occur independently or concomitantly with genotoxic effects. While data on genotoxicity are collected through standardized guideline tests, data collected on epigenetic effects is far less uniform. In 2016, we conducted a systematic review of published studies of genotoxic carcinogens that reported epigenetic endpoints to better understand the evidence for epigenetic alterations of human carcinogens, and the potential association with genotoxic endpoints. Since then, the number of studies of epigenetic effects of chemicals has nearly doubled. This review stands as an update on epigenetic alterations induced by occupational and environmental human carcinogens that were previously and recently classified as Group 1 by the International Agency for Research on Cancer. We found that the evidence of epigenetic effects remains uneven across agents. Studies of DNA methylation are most abundant, while reports concerning effects on non-coding RNA have increased over the past 5 years. By contrast, mechanistic toxicology studies of histone modifications and chromatin state alterations remain few. We found that most publications of epigenetic effects of carcinogens were studies in exposed humans or human cells. Studies in rodents represent the second most common species used for epigenetic studies in toxicology, in vivo exposures being the most predominant. Future studies should incorporate dose- and time-dependent study designs and also investigate the persistence of effects following cessation of exposure, considering the dynamic nature of most epigenetic alterations.