Short-term toxicogenomics as an alternative approach to chronic in vivo studies for derivation of points of departure: A case study in the rat with a triazole fungicide

Far from Their Origins: A Transcriptomic Investigation on How 2,4-Di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) Phenol Affects Rainbow Trout Alevins

Benzotriazole ultraviolet stabilizers (BUVSs) are a group of widely used chemicals added to a variety of consumer (e.g., plastics) and industrial (e.g., metal coating) goods. Although detected globally as an environmentally persistent pollutant, BUVSs have received relatively little toxicological attention and only recently have been acknowledged to affect development and the endocrine system in vivo. In our previous study, altered behavior, indicative of potential neurotoxicity, was observed among rainbow trout alevins (day 14 posthatching) that were microinjected as embryos with a single environmentally relevant dose of 2,4-di-tert-butyl-6-(5-chloro-2H-benzotriazol-2-yl) phenol (UV-327). In the present follow-up study, we performed whole-transcriptome profiling (RNA sequencing) of newly hatched alevins from the same batch. The primary aim was to identify biomarkers related to behavior and neurology. Dose-specifically, 1 to 176 differentially expressed genes (DEGs) were identified. In the group presenting altered behavior (273.4 ng g-1), 176 DEGs were identified, yet only a fraction was related to neurological functions, including water, calcium, and potassium homeostasis; acetylcholine transmission and signaling; as well insulin and energy metabolism. The second objective was to estimate the transcriptomic point of departure (tPOD) and assess if point estimate(s) are protective of altered behavior. A tPOD was established at 35 to 94 ng UV-327 g-1 egg, making this tPOD protective of behavioral alterations. Holistically, these transcriptomic alterations provide a foundation for future research on how BUVSs can influence rainbow trout alevin development, while providing support to the hypothesis that UV-327 can influence neurogenesis and subsequent behavioral endpoints. The exact structural and functional changes caused by embryonic exposure to UV-327 remain enigmatic and will require extensive investigation before being deciphered and understood toxicologically. Environ Toxicol Chem 2024;00:1-12. © 2024 The Author(s). Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Transcriptomic point of departure determination: a comparison of distribution-based and gene set-based approaches

A key step in assessing the potential human and environmental health risks of industrial and agricultural chemicals is to determine the toxicity point of departure (POD), which is the highest dose level that causes no adverse effect. Transcriptomic POD (tPOD) values have been suggested to accurately estimate toxicity POD values. One step in the most common approach for tPOD determination involves mapping genes to annotated gene sets, a process that might lead to substantial information loss particularly in species with poor gene annotation. Alternatively, methods that calculate tPOD values directly from the distribution of individual gene POD values omit this mapping step. Using rat transcriptome data for 79 molecules obtained from Open TG-GATEs (Toxicogenomics Project Genomics Assisted Toxicity Evaluation System), the hypothesis was tested that methods based on the distribution of all individual gene POD values will give a similar tPOD value to that obtained via the gene set-based method. Gene set-based tPOD values using four different gene set structures were compared to tPOD values from five different individual gene distribution methods. Results revealed a high tPOD concordance for all methods tested, especially for molecules with at least 300 dose-responsive probesets: for 90% of those molecules, the tPOD values from all methods were within 4-fold of each other. In addition, random gene sets based upon the structure of biological knowledge-derived gene sets produced tPOD values with a median absolute fold change of 1.3-1.4 when compared to the original biological knowledge-derived gene set counterparts, suggesting that little biological information is used in the gene set-based tPOD generation approach. These findings indicate using individual gene distributions to calculate a tPOD is a viable and parsimonious alternative to using gene sets. Importantly, individual gene distribution-based tPOD methods do not require knowledge of biological organization and can be applied to any species including those with poorly annotated gene sets.

Integrating high-throughput phenotypic profiling and transcriptomic analyses to predict the hepatosteatosis effects induced by per- and polyfluoroalkyl substances

Per- and polyfluoroalkyl substances (PFAS) is a large compound class (n > 12,000) that is extensively present in food, drinking water, and aquatic environments. Reduced serum triglycerides and hepatosteatosis appear to be the common phenotypes for different PFAS chemicals. However, the hepatosteatosis potential of most PFAS chemicals remains largely unknown. This study aims to investigate PFAS-induced hepatosteatosis using in vitro high-throughput phenotype profiling (HTPP) and high-throughput transcriptomic (HTTr) data. We quantified the in vitro hepatosteatosis effects and mitochondrial damage using high-content imaging, curated the transcriptomic data from the Gene Expression Omnibus (GEO) database, and then calculated the point of departure (POD) values for HTPP phenotypes or HTTr transcripts, using the Bayesian benchmark dose modeling approach. Our results indicated that PFAS compounds with fully saturated C-F bonds, sulfur- and nitrogen-containing functional groups, and a fluorinated carbon chain length greater than 8 have the potential to produce biological effects consistent with hepatosteatosis. PFAS primarily induced hepatosteatosis via disturbance in lipid transport and storage. The potency rankings of PFAS compounds are highly concordant among in vitro HTPP, HTTr, and in vivo hepatosteatosis phenotypes (ρ = 0.60-0.73). In conclusion, integrating the information from in vitro HTPP and HTTr analyses can accurately project in vivo hepatosteatosis effects induced by PFAS compounds.

Linking Mechanistic Effects of Pharmaceuticals and Personal Care Products to Ecologically Relevant Outcomes: A Decade of Progress

There are insufficient toxicity data to assess the ecological risks of many pharmaceuticals and personal care products (PPCPs). While data limitations are not uncommon for contaminants of environmental concern, PPCPs are somewhat unique in that an a priori understanding of their biological activities in conjunction with measurements of molecular, biochemical, or histological responses could provide a foundation for understanding mode(s) of action and predicting potential adverse apical effects. Over the past decade significant progress has been made in the development of new approach methodologies (NAMs) to efficiently quantify these types of endpoints using computational models and pathway-based in vitro and in vivo assays. The availability of open-access knowledgebases to curate biological response (including NAM) data and sophisticated bioinformatics tools to help interpret the information also has significantly increased. Finally, advances in the development and implementation of the adverse outcome pathway framework provide the critical conceptual underpinnings needed to translate NAM data into predictions of the ecologically relevant outcomes required by risk assessors and managers. The evolution and convergence of these various data streams, tools, and concepts provides the basis for a fundamental change in how ecological risks of PPCPs can be pragmatically assessed. Environ Toxicol Chem 2024;43:537-548. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Sub-chronic exposure to hexaconazole affects the lipid metabolism of rats through mTOR-PPAR-γ/SREBP1 signaling pathway mediated by oxidative stress

Hexaconazole (Hex) is a widely used and high frequency detected triazole fungicide in agricultural products and environment which may pose potential toxicity to the nontargeted organisms. Hex had been reported to affect lipid homeostasis while the mechanism was undefined. This study aims to explore the characteristic lipidomic profiles and clarify the underlying signaling pathways of Hex-induced lipid metabolism disorder in rat liver. The results showed that sub-chronic exposure to environmental related concentrations of Hex caused histopathological changes, oxidative stress, fat accumulation, lipid biochemical parameter increase in rats. Moreover, the untargeted lipidomic analysis showed that the levels of TAG, PC, and PE and the pathway of glycerophospholipid metabolism were heavily altered by Hex. We further analyzed the lipid metabolism related genes and proteins which revealed that Hex exposure increased amount of lipogenesis by activating oxidative stress-mediated mTOR-PPAR-γ/SREBP1 signaling pathways. The imbalance of lipid homeostasis induced by Hex exposure might further lead to obesity, cardiovascular diseases (CVDs), and hyperlipidemia. Our results provided systematic and comprehensive evidence for the mechanism of Hex-induced lipid metabolism disorder at environmental concentrations and supplied a certain basis for its health risks assessment.

Short-Term Transcriptomic Points of Departure Are Consistent with Chronic Points of Departure for Three Organophosphate Pesticides across Mouse and Fathead Minnow

New approach methods (NAMs) can reduce the need for chronic animal studies. Here, we apply benchmark dose (concentration) (BMD(C))-response modeling to transcriptomic changes in the liver of mice and in fathead minnow larvae after short-term exposures (7 days and 1 day, respectively) to several dose/concentrations of three organophosphate pesticides (OPPs): fenthion, methidathion, and parathion. The mouse liver transcriptional points of departure (TPODs) for fenthion, methidathion, and parathion were 0.009, 0.093, and 0.046 mg/Kg-bw/day, while the fathead minnow larva TPODs were 0.007, 0.115, and 0.046 mg/L, respectively. The TPODs were consistent across both species and reflected the relative potencies from traditional chronic toxicity studies with fenthion identified as the most potent. Moreover, the mouse liver TPODs were more sensitive than or within a 10-fold difference from the chronic apical points of departure (APODs) for mammals, while the fathead minnow larva TPODs were within an 18-fold difference from the chronic APODs for fish species. Short-term exposure to OPPs significantly impacted acetylcholinesterase mRNA abundance (FDR p-value <0.05, |fold change| ≥2) and canonical pathways (IPA, p-value <0.05) associated with organism death and neurological/immune dysfunctions, indicating the conservation of key events related to OPP toxicity. Together, these results build confidence in using short-term, molecular-based assays for the characterization of chemical toxicity and risk, thereby reducing reliance on chronic animal studies.

State of the Art of Genomic Technology in Toxicology: A Review

The rapid growth of genomics techniques has revolutionized and impacted, greatly and positively, the knowledge of toxicology, ushering it into a "new era": the era of genomic technology (GT). This great advance permits us to analyze the whole genome, to know the gene response to toxicants and environmental stressors, and to determine the specific profiles of gene expression, among many other approaches. The aim of this work was to compile and narrate the recent research on GT during the last 2 years (2020-2022). A literature search was managed using the PubMed and Medscape interfaces on the Medline database. Relevant articles published in peer-reviewed journals were retrieved and their main results and conclusions are mentioned briefly. It is quite important to form a multidisciplinary taskforce on GT with the aim of designing and implementing a comprehensive, collaborative, and a strategic work plan, prioritizing and assessing the most relevant diseases, so as to decrease human morbimortality due to exposure to environmental chemicals and stressors.

From vision toward best practices: Evaluating in vitro transcriptomic points of departure for application in risk assessment using a uniform workflow

The growing number of chemicals in the current consumer and industrial markets presents a major challenge for regulatory programs faced with the need to assess the potential risks they pose to human and ecological health. The increasing demand for hazard and risk assessment of chemicals currently exceeds the capacity to produce the toxicity data necessary for regulatory decision making, and the applied data is commonly generated using traditional approaches with animal models that have limited context in terms of human relevance. This scenario provides the opportunity to implement novel, more efficient strategies for risk assessment purposes. This study aims to increase confidence in the implementation of new approach methods in a risk assessment context by using a parallel analysis to identify data gaps in current experimental designs, reveal the limitations of common approaches deriving transcriptomic points of departure, and demonstrate the strengths in using high-throughput transcriptomics (HTTr) to derive practical endpoints. A uniform workflow was applied across six curated gene expression datasets from concentration-response studies containing 117 diverse chemicals, three cell types, and a range of exposure durations, to determine tPODs based on gene expression profiles. After benchmark concentration modeling, a range of approaches was used to determine consistent and reliable tPODs. High-throughput toxicokinetics were employed to translate in vitro tPODs (µM) to human-relevant administered equivalent doses (AEDs, mg/kg-bw/day). The tPODs from most chemicals had AEDs that were lower (i.e., more conservative) than apical PODs in the US EPA CompTox chemical dashboard, suggesting in vitro tPODs would be protective of potential effects on human health. An assessment of multiple data points for single chemicals revealed that longer exposure duration and varied cell culture systems (e.g., 3D vs. 2D) lead to a decreased tPOD value that indicated increased chemical potency. Seven chemicals were flagged as outliers when comparing the ratio of tPOD to traditional POD, thus indicating they require further assessment to better understand their hazard potential. Our findings build confidence in the use of tPODs but also reveal data gaps that must be addressed prior to their adoption to support risk assessment applications.

The contribution of larval zebrafish transcriptomics to chemical risk assessment

In Canada, the Canadian Environmental Protection Act (1999) requires human health and environmental risk assessments be conducted for new substances prior to their manufacture or import. While this toxicity data is historically obtained using rodents, in response to the international effort to eliminate animal testing, Health Canada is collaborating with the National Research Council (NRC) of Canada to develop a New Approach Method by refining existing NRC zebrafish models. The embryo/larval zebrafish model evaluates systemic (whole body) general toxicity which is currently unachievable with cell-based testing. The model is strengthened using behavioral, toxicokinetic and transcriptomic responses to assess non-visible indicators of toxicity following chemical exposure at sub-phenotypic concentrations. In this paper, the predictive power of zebrafish transcriptomics is demonstrated using two chemicals; Raloxifene and Resorcinol. Raloxifene exposure produced darkening of the liver and malformation of the nose/mandible, while Resorcinol exposure produced increased locomotor activity. Transcriptomic analysis correlated differentially expressed genes with the phenotypic effects and benchmark dose calculations determined that the transcriptomic Point of Departure (POD) occurred at subphenotypic concentrations. Correlating gene expression with apical (phenotypic) effects strengthens confidence in evaluation of chemical toxicity, thereby demonstrating the significant advancement that the larval zebrafish transcriptomics model represents in chemical risk assessment.

Towards achieving a modern science-based paradigm for agrochemical carcinogenicity assessment

The rodent cancer bioassay has been the standard approach to fulfill regulatory requirements for assessing human carcinogenic potential of agrochemicals, food additives, industrial chemicals, and pharmaceuticals. Decades of research have described the limitations of the rodent cancer bioassay leading to international initiatives to seek alternatives and establish approaches that modernize carcinogenicity assessment. Biologically relevant approaches can provide mechanistic information and increased efficiency for evaluating hazard and risk of chemical carcinogenicity to humans. The application of human-relevant mechanistic understanding to support new approaches to carcinogenicity assessment will be invaluable for regulatory decision-making. The present work outlines the challenges and opportunities that authorities should consider as they come together to build a roadmap that leads to global acceptance and incorporation of fit-for-purpose, scientifically defensible new approaches for human-relevant carcinogenicity assessment of agrochemicals.

Transcriptomic Points of Departure Calculated from Rainbow Trout Gill, Liver, and Gut Cell Lines Exposed to Methylmercury and Fluoxetine

Ethical and resource limitation concerns are pushing chemicals management to develop alternatives to animal testing strategies. The objective of our study was to determine whether transcriptomic point of departure (tPOD) values could be derived from studies that followed Organisation for Economic Co-operation and Development (OECD) Test No. 249 (rainbow trout gill cell line), as well as from studies on trout liver and gut cells. Gill, liver, and gut cell lines were exposed to methylmercury and fluoxetine. Concentrations causing 50% cytotoxicity (LC50) were derived, the whole transcriptome was sequenced, and gene tPOD and pathway benchmark dose (BMD) values were derived from transcriptomic dose-response analysis. Differences in LC50 and transcriptomic responses across the cell lines were noted. For methylmercury, the tPOD_mode values were 14.5, 20.5, and 17.8 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway (pathway BMDs in parentheses) was ferroptosis in the gill (3.1 ppb) and liver (3.5 ppb), and glutathione metabolism in the gut (6.6 ppb). For fluoxetine, the tPOD_mode values were 109.4, 108.4, and 97.4 ppb for the gill, liver, and gut cells, respectively. The most sensitive pathway was neurotrophin signaling in the gill (147 ppb) and dopaminergic signaling in the gut (86.3 ppb). For both chemicals, the gene tPOD and pathway BMD values were lower than cytotoxic concentrations in vitro, and within 10-fold below the in vivo LC50s. By bringing together transcriptomics and dose-response analysis with an OECD test method in three cell lines, the results help to establish an in vitro method yielding tPOD values that are hypothesized to be protective of in vivo concentrations associated with adverse outcomes, and also give insights into mechanisms of action. Environ Toxicol Chem 2022;41:1982-1992. © 2022 SETAC.

A microRNA or messenger RNA point of departure estimates an apical endpoint point of departure in a rat developmental toxicity model

Traditional developmental toxicity testing practice examines fetal apical endpoints to identify a point of departure (POD) for risk assessment. A potential new testing paradigm involves deriving a POD from a comprehensive analysis of molecular-level change. Here, the rat ketoconazole endocrine-mediated developmental toxicity model was used to test the hypothesis that maternal epigenomic (miRNA) and transcriptomic (mRNA) PODs are similar to fetal apical endpoint PODs. Sprague-Dawley rats were exposed from gestation day (GD) 6-21 to 0, 0.063, 0.2, 0.63, 2, 6.3, 20, or 40 mg/kg/day ketoconazole. Dam systemic, liver, and placenta PODs, along with GD 21 fetal resorption, body weight, and skeletal apical PODs were derived using BMDS software. GD 21 dam liver and placenta miRNA and mRNA PODs were obtained using three methods: a novel individual molecule POD accumulation method, a first mode method, and a gene set method. Dam apical POD values ranged from 2.0 to 38.6 mg/kg/day; the lowest value was for placenta histopathology. Fetal apical POD values were 10.9-20.3 mg/kg/day; the lowest value was for fetal resorption. Dam liver miRNA and mRNA POD values were 0.34-0.69 mg/kg/day, and placenta miRNA and mRNA POD values were 2.53-6.83 mg/kg/day. Epigenomic and transcriptomic POD values were similar across liver and placenta. Deriving a molecular POD from dam liver or placenta was protective of a fetal apical POD. These data support the conclusion that a molecular POD can be used to estimate, or be protective of, a developmental toxicity apical POD.

Rethinking chronic toxicity and carcinogenicity assessment for agrochemicals project (ReCAAP): A reporting framework to support a weight of evidence safety assessment without long-term rodent bioassays

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.

Using Transcriptomics and Metabolomics to Understand Species Differences in Sensitivity to Chlorpyrifos in Japanese Quail and Double-Crested Cormorant Embryos

Modern 21st-century toxicity testing makes use of omics technologies to address critical questions in toxicology and chemical management. Of interest are questions relating to chemical mechanisms of toxicity, differences in species sensitivity, and translation of molecular effects to observable apical endpoints. Our study addressed these questions by comparing apical outcomes and multiple omics responses in early-life stage exposure studies with Japanese quail (Coturnix japonica) and double-crested cormorant (Phalacrocorax auritus), representing a model and ecological species, respectively. Specifically, we investigated the dose-dependent response of apical outcomes as well as transcriptomics and metabolomics in the liver of each species exposed to chlorpyrifos, a widely used organophosphate pesticide. Our results revealed a clear pattern of dose-dependent disruption of gene expression and metabolic profiles in Japanese quail but not double-crested cormorant at similar chlorpyrifos exposure concentrations. The difference in sensitivity between species was likely due to higher metabolic transformation of chlorpyrifos in Japanese quail compared to double-crested cormorant. The most impacted biological pathways after chlorpyrifos exposure in Japanese quail included hepatic metabolism, oxidative stress, endocrine disruption (steroid and nonsteroid hormones), and metabolic disease (lipid and fatty acid metabolism). Importantly, we show consistent responses across biological scales, suggesting that significant disruption at the level of gene expression and metabolite profiles leads to observable apical responses at the organism level. Our study demonstrates the utility of evaluating effects at multiple biological levels of organization to understand how modern toxicity testing relates to outcomes of regulatory relevance, while also highlighting important, yet poorly understood, species differences in sensitivity to chemical exposure. Environ Toxicol Chem 2021;40:3019-3033. © 2021 SETAC.

In Silico Approaches In Carcinogenicity Hazard Assessment: Current Status and Future Needs

Historically, identifying carcinogens has relied primarily on tumor studies in rodents, which require enormous resources in both money and time. In silico models have been developed for predicting rodent carcinogens but have not yet found general regulatory acceptance, in part due to the lack of a generally accepted protocol for performing such an assessment as well as limitations in predictive performance and scope. There remains a need for additional, improved in silico carcinogenicity models, especially ones that are more human-relevant, for use in research and regulatory decision-making. As part of an international effort to develop in silico toxicological protocols, a consortium of toxicologists, computational scientists, and regulatory scientists across several industries and governmental agencies evaluated the extent to which in silico models exist for each of the recently defined 10 key characteristics (KCs) of carcinogens. This position paper summarizes the current status of in silico tools for the assessment of each KC and identifies the data gaps that need to be addressed before a comprehensive in silico carcinogenicity protocol can be developed for regulatory use.

Transcriptomic analyses of livers from mice exposed to 1,4-dioxane for up to 90 days to assess potential mode(s) of action underlying liver tumor development

1,4-Dioxane is a volatile organic compound with industrial and commercial applications as a solvent and in the manufacture of other chemicals. 1,4-Dioxane has been demonstrated to induce liver tumors in chronic rodent bioassays conducted at very high doses. The available evidence for 1,4-dioxane-induced liver tumors in rodents aligns with a threshold-dependent mode of action (MOA), with the underlying mechanism being less clear in the mouse than in rats. To gain a better understanding of the underlying molecular mechanisms related to liver tumor development in mice orally exposed to 1,4-dioxane, transcriptomics analysis was conducted on liver tissue collected from a 90-day drinking water study in female B6D2F₁/Crl mice (Lafranconi et al., 2020). Using tissue samples from female mice exposed to 1,4-dioxane in the drinking water at concentrations of 0, 40, 200, 600, 2,000 or 6,000 ppm for 7, 28, and 90 days, transcriptomic analyses demonstrate minimal treatment effects on global gene expression at concentrations below 600 ppm. At higher concentrations, genes involved in phase II metabolism and mitotic cell cycle checkpoints were significantly upregulated. There was an overall lack of enrichment of genes related to DNA damage response. The increase in mitotic signaling is most prevalent in the livers of mice exposed to 1,4-dioxane at the highest concentrations for 90 days. This finding aligns with phenotypic changes reported by Lafranconi et al. (2020) after 90-days of exposure to 6,000 ppm 1,4-dioxane in the same tissues. The transcriptomics analysis further supports overarching study findings demonstrating a non-mutagenic, threshold-based, mitogenic MOA for 1,4-dioxane-induced liver tumors.

A rat subchronic study transcriptional point of departure estimates a carcinogenicity study apical point of departure

Considerations of human relevance and animal use are driving research to identify new approaches to inform risk assessment of chemicals and replace guideline-based rodent carcinogenicity tests. Here, the hypothesis was tested across four agrochemicals that 1) a rat 90-day transcriptome-based BEPOD is protective of a rat carcinogenicity study and 2) a subchronic liver or kidney BEPOD would approximate a cancer bioassay apical POD derived from other organs and a rat subchronic BEPOD would approximate a mouse cancer bioassay apical POD. Using RNA sequencing and BMDExpress software, liver and/or kidney BEPOD values were generated in male rats exposed for 90 days to either Triclopyr Acid, Pronamide, Sulfoxaflor, or Fenpicoxamid. BEPOD values were compared to benchmark dose-derived apical POD values generated from rat 90-day and rodent carcinogenicity studies. Across all four agrochemicals, findings showed that a rat 90-day study BEPOD approximated the most sensitive apical POD (within 10-fold) generated from the 90-day rat study and long-term rodent carcinogenicity studies. This study supports the conclusion that a subchronic transcriptome-based BEPOD could be utilized to estimate an apical POD within a risk-based approach of chronic toxicity and carcinogenicity agrochemical assessment, abrogating the need for time- and resource-intensive rodent carcinogenicity studies and minimizing animal testing.