ABSTRACT
Toxicogenomics measures molecular features, such as transcripts, proteins, metabolites and epigenomic modifications, to understand and predict the toxicological effects of environmental and pharmaceutical exposures. Transcriptomics has become an integral tool in contemporary toxicology research owing to innovations in gene expression profiling that can provide mechanistic and quantitative information at scale. These data can be used to predict toxicological hazards through the use of transcriptomic biomarkers, network inference analyses, pattern-matching approaches and artificial intelligence. Furthermore, emerging approaches, such as high-throughput dose-response modelling, can leverage toxicogenomic data for human health protection even in the absence of predicting specific hazards. Finally, single-cell transcriptomics and multi-omics provide detailed insights into toxicological mechanisms. Here, we review the progress since the inception of toxicogenomics in applying transcriptomics towards toxicology testing and highlight advances that are transforming risk assessment.
ABSTRACT
New approach methodologies (NAMs) aim to accelerate the pace of chemical risk assessment while simultaneously reducing cost and dependency on animal studies. High Throughput Transcriptomics (HTTr) is an emerging NAM in the field of chemical hazard evaluation for establishing in vitro points-of-departure and providing mechanistic insight. In the current study, 1201 test chemicals were screened for bioactivity at eight concentrations using a 24-h exposure duration in the human- derived U-2 OS osteosarcoma cell line with HTTr. Assay reproducibility was assessed using three reference chemicals that were screened on every assay plate. The resulting transcriptomics data were analyzed by aggregating signal from genes into signature scores using gene set enrichment analysis, followed by concentration-response modeling of signatures scores. Signature scores were used to predict putative mechanisms of action, and to identify biological pathway altering concentrations (BPACs). BPACs were consistent across replicates for each reference chemical, with replicate BPAC standard deviations as low as 5.6 × 10-3 µM, demonstrating the internal reproducibility of HTTr-derived potency estimates. BPACs of test chemicals showed modest agreement (R2 = 0.55) with existing phenotype altering concentrations from high throughput phenotypic profiling using Cell Painting of the same chemicals in the same cell line. Altogether, this HTTr based chemical screen contributes to an accumulating pool of publicly available transcriptomic data relevant for chemical hazard evaluation and reinforces the utility of cell based molecular profiling methods in estimating chemical potency and predicting mechanism of action across a diverse set of chemicals.
ABSTRACT
High-throughput transcriptomics (HTTr) uses gene expression profiling to characterize the biological activity of chemicals in in vitro cell-based test systems. As an extension of a previous study testing 44 chemicals, HTTr was used to screen an additional 1751 unique chemicals from the EPA's ToxCast collection in MCF7 cells using eight concentrations and an exposure duration of 6 hours. We hypothesized that concentration-response modeling of signature scores could be used to identify putative molecular targets and cluster chemicals with similar bioactivity. Clustering and enrichment analyses were conducted based on signature catalog annotations and ToxPrint chemotypes to facilitate molecular target prediction and grouping of chemicals with similar bioactivity profiles. Enrichment analysis based on signature catalog annotation identified known mechanisms-of-action (MeOAs) associated with well-studied chemicals and generated putative MeOAs for other active chemicals. Chemicals with predicted MeOAs included those targeting estrogen receptor (ER), glucocorticoid receptor (GR), retinoic acid receptor (RAR), the NRF2/KEAP/ARE pathway, AP-1 activation and others. Using reference chemicals for ER modulation, the study demonstrated that HTTr in MCF7 cells was able to stratify chemicals in terms of agonist potency, distinguish ER agonists from antagonists, and cluster chemicals with similar activities as predicted by the ToxCast ER Pathway model. Uniform manifold approximation and projection (UMAP) embedding of signature-level results identified novel ER modulators with no ToxCast ER Pathway model predictions. Finally, UMAP combined with ToxPrint chemotype enrichment was used to explore the biological activity of structurally-related chemicals. The study demonstrates that HTTr can be used to inform chemical risk assessment by determining in vitro points-of-departure, predicting chemicals' molecular mechanism(s)-of-action (MeOA) and grouping chemicals with similar bioactivity profiles.
ABSTRACT
Utilization of in vitro (cellular) techniques, like Cell Painting and transcriptomics, could provide powerful tools for agrochemical candidate sorting and selection in the discovery process. However, using these models generates challenges translating in vitro concentrations to the corresponding in vivo exposures. Physiologically based pharmacokinetic (PBPK) modeling provides a framework for quantitative in vitro to in vivo extrapolation (IVIVE). We tested whether in vivo (rat liver) transcriptomic and apical points of departure (PODs) could be accurately predicted from in vitro (rat hepatocyte or human HepaRG) transcriptomic PODs or HepaRG Cell Painting PODs using PBPK modeling. We compared two PBPK models, the ADMET predictor and the httk R package, and found httk to predict the in vivo PODs more accurately. Our findings suggest that a rat liver apical and transcriptomic POD can be estimated utilizing a combination of in vitro transcriptome-based PODs coupled with PBPK modeling for IVIVE. Thus, high content in vitro data can be translated with modest accuracy to in vivo models of ultimate regulatory importance to help select agrochemical analogs in early stage discovery program.
Subject(s)
Agrochemicals , Animals , Rats , Humans , Agrochemicals/pharmacokinetics , Agrochemicals/toxicity , Hepatocytes/metabolism , Liver/metabolism , Models, Biological , Male , Transcriptome , Cell Line , Risk AssessmentABSTRACT
Single, high doses of TCDD in rats are known to cause wasting, a progressive loss of 30 to 50% body weight and death within several weeks. To identify pathway perturbations at or near doses causing wasting, we examined differentially gene expression (DGE) and pathway enrichment in centrilobular (CL) and periportal (PP) regions of female rat livers following 6 dose levels of TCDD - 0, 3, 22, 100, 300, and 1000 ng/kg/day, 5 days/week for 4 weeks. At the higher doses, rats lost weight, had increased liver/body weight ratios and nearly complete cessation of liver cell proliferation, signs consistent with wasting. DGE curves were left shifted for the CL versus the PP regions. Canonical Phase I and Phase II genes were maximally increased at lower doses and remained elevated at all doses. At lower doses, ≤ 22 ng/kg/day in the CL and ≤ 100 ng/kg/day, upregulated genes showed transcription factor (TF) enrichment for AHR and ARNT. At the mid- and high-dose doses, there was a large number of downregulated genes and pathway enrichment for DEGs which showed downregulation of many cellular metabolism processes including those for steroids, fatty acid metabolism, pyruvate metabolism and citric acid cycle. There was significant TF enrichment of the hi-dose downregulated genes for RXR, ESR1, LXR, PPARalpha. At the highest dose, there was also pathway enrichment with upregulated genes for extracellular matrix organization, collagen formation, hemostasis and innate immune system. TCDD demonstrates most of its effects through binding the aryl hydrocarbon receptor (AHR) while the downregulation of metabolism genes at higher TCDD doses is known to be independent of AHR binding to DREs. Based on our results with DEG, we provide a hypothesis for wasting in which high doses of TCDD shift circadian processes away from the resting state, leading to greatly reduced synthesis of steroids and complex lipids needed for cell growth, and producing gene expression signals consistent with an epithelial-to-mesenchymal transition in hepatocytes.
Subject(s)
Aryl Hydrocarbon Receptor Nuclear Translocator , Liver , Polychlorinated Dibenzodioxins , Receptors, Aryl Hydrocarbon , Animals , Receptors, Aryl Hydrocarbon/genetics , Receptors, Aryl Hydrocarbon/metabolism , Female , Liver/drug effects , Liver/metabolism , Liver/pathology , Polychlorinated Dibenzodioxins/toxicity , Rats , Aryl Hydrocarbon Receptor Nuclear Translocator/genetics , Aryl Hydrocarbon Receptor Nuclear Translocator/metabolism , Circadian Rhythm/drug effects , Circadian Rhythm/genetics , Basic Helix-Loop-Helix Transcription Factors/genetics , Basic Helix-Loop-Helix Transcription Factors/metabolism , Chemical and Drug Induced Liver Injury/genetics , Chemical and Drug Induced Liver Injury/metabolism , Chemical and Drug Induced Liver Injury/pathology , Transcriptome/drug effects , Gene Expression Profiling/methods , Rats, Sprague-Dawley , Dose-Response Relationship, DrugABSTRACT
Per- and polyfluoroalkyl substances (PFAS) are widely used, and their fluorinated state contributes to unique uses and stability but also long half-lives in the environment and humans. PFAS have been shown to be toxic, leading to immunosuppression, cancer, and other adverse health outcomes. Only a small fraction of the PFAS in commerce have been evaluated for toxicity using in vivo tests, which leads to a need to prioritize which compounds to examine further. Here, we demonstrate a prioritization approach that combines human biomonitoring data (blood concentrations) with bioactivity data (concentrations at which bioactivity is observed in vitro) for 31 PFAS. The in vitro data are taken from a battery of cell-based assays, mostly run on human cells. The result is a Bioactive Concentration to Blood Concentration Ratio (BCBCR), similar to a margin of exposure (MoE). Chemicals with low BCBCR values could then be prioritized for further risk assessment. Using this method, two of the PFAS, PFOA (Perfluorooctanoic Acid) and PFOS (Perfluorooctane Sulfonic Acid), have BCBCR values < 1 for some populations. An additional 9 PFAS have BCBCR values < 100 for some populations. This study shows a promising approach to screening level risk assessments of compounds such as PFAS that are long-lived in humans and other species.
ABSTRACT
Per- and polyfluoroalkyl substances (PFAS) represent a large group of contaminants of concern based on their widespread use, environmental persistence, and potential toxicity. Many traditional models for estimating toxicity, bioaccumulation, and other toxicological properties are not well suited for PFAS. Consequently, there is a need to generate hazard information for PFAS in an efficient and cost-effective manner. In the present study, Daphnia magna were exposed to multiple concentrations of 22 different PFAS for 24 h in a 96-well plate format. Following exposure, whole-body RNA was extracted and extracts, each representing five exposed individuals, were subjected to RNA sequencing. Following analytical measurements to verify PFAS exposure concentrations and quality control on processed cDNA libraries for sequencing, concentration-response modeling was applied to the data sets for 18 of the tested compounds, and the concentration at which a concerted molecular response occurred (transcriptomic point of departure; tPOD) was calculated. The tPODs, based on measured concentrations of PFAS, generally ranged from 0.03 to 0.58 µM (9.9-350 µg/L; interquartile range). In most cases, these concentrations were two orders of magnitude lower than similarly calculated tPODs for human cell lines exposed to PFAS. They were also lower than apical effect concentrations reported for seven PFAS for which some crustacean or invertebrate toxicity data were available, although there were a few exceptions. Despite being lower than most other available hazard benchmarks, D. magna tPODs were, on average, four orders of magnitude greater than the maximum aqueous concentrations of PFAS measured in Great Lakes tributaries. Overall, this high-throughput transcriptomics assay with D. magna holds promise as a component of a tiered hazard evaluation strategy employing new approach methodologies. Environ Toxicol Chem 2024;00:1-16. © 2024 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.
ABSTRACT
The presence of numerous chemical contaminants from industrial, agricultural, and pharmaceutical sources in water supplies poses a potential risk to human and ecological health. Current chemical analyses suffer from limitations, including chemical coverage and high cost, and broad-coverage in vitro assays such as transcriptomics may further improve water quality monitoring by assessing a large range of possible effects. Here, we used high-throughput transcriptomics to assess the activity induced by field-derived water extracts in MCF7 breast carcinoma cells. Wastewater and surface water extracts induced the largest changes in expression among cell proliferation-related genes and neurological, estrogenic, and antibiotic pathways, whereas drinking and reclaimed water extracts that underwent advanced treatment showed substantially reduced bioactivity on both gene and pathway levels. Importantly, reclaimed water extracts induced fewer changes in gene expression than laboratory blanks, which reinforces previous conclusions based on targeted assays and improves confidence in bioassay-based monitoring of water quality.
Subject(s)
Water Pollutants, Chemical , Water Purification , Humans , Environmental Monitoring , Water Pollutants, Chemical/analysis , Water Quality , Gene Expression Profiling , Biological AssayABSTRACT
Multiple new approach methods (NAMs) are being developed to rapidly screen large numbers of chemicals to aid in hazard evaluation and risk assessments. High-throughput transcriptomics (HTTr) in human cell lines has been proposed as a first-tier screening approach for determining the types of bioactivity a chemical can cause (activation of specific targets vs. generalized cell stress) and for calculating transcriptional points of departure (tPODs) based on changes in gene expression. In the present study, we examine a range of computational methods to calculate tPODs from HTTr data, using six data sets in which MCF7 cells cultured in two different media formulations were treated with a panel of 44 chemicals for 3 different exposure durations (6, 12, 24 hr). The tPOD calculation methods use data at the level of individual genes and gene set signatures, and compare data processed using the ToxCast Pipeline 2 (tcplfit2), BMDExpress and PLIER (Pathway Level Information ExtractoR). Methods were evaluated by comparing to in vitro PODs from a validated set of high-throughput screening (HTS) assays for a set of estrogenic compounds. Key findings include: (1) for a given chemical and set of experimental conditions, tPODs calculated by different methods can vary by several orders of magnitude; (2) tPODs are at least as sensitive to computational methods as to experimental conditions; (3) in comparison to an external reference set of PODs, some methods give generally higher values, principally PLIER and BMDExpress; and (4) the tPODs from HTTr in this one cell type are mostly higher than the overall PODs from a broad battery of targeted in vitro ToxCast assays, reflecting the need to test chemicals in multiple cell types and readout technologies for in vitro hazard screening.
Subject(s)
Gene Expression Profiling , Transcriptome , Humans , High-Throughput Screening Assays/methods , Estrogens , Cell Line , Risk Assessment/methodsABSTRACT
The predominantly animal-centric approach of chemical safety assessment has increasingly come under pressure. Society is questioning overall performance, sustainability, continued relevance for human health risk assessment and ethics of this system, demanding a change of paradigm. At the same time, the scientific toolbox used for risk assessment is continuously enriched by the development of "New Approach Methodologies" (NAMs). While this term does not define the age or the state of readiness of the innovation, it covers a wide range of methods, including quantitative structure-activity relationship (QSAR) predictions, high-throughput screening (HTS) bioassays, omics applications, cell cultures, organoids, microphysiological systems (MPS), machine learning models and artificial intelligence (AI). In addition to promising faster and more efficient toxicity testing, NAMs have the potential to fundamentally transform today's regulatory work by allowing more human-relevant decision-making in terms of both hazard and exposure assessment. Yet, several obstacles hamper a broader application of NAMs in current regulatory risk assessment. Constraints in addressing repeated-dose toxicity, with particular reference to the chronic toxicity, and hesitance from relevant stakeholders, are major challenges for the implementation of NAMs in a broader context. Moreover, issues regarding predictivity, reproducibility and quantification need to be addressed and regulatory and legislative frameworks need to be adapted to NAMs. The conceptual perspective presented here has its focus on hazard assessment and is grounded on the main findings and conclusions from a symposium and workshop held in Berlin in November 2021. It intends to provide further insights into how NAMs can be gradually integrated into chemical risk assessment aimed at protection of human health, until eventually the current paradigm is replaced by an animal-free "Next Generation Risk Assessment" (NGRA).
Subject(s)
Artificial Intelligence , Toxicity Tests , Humans , Reproducibility of Results , Toxicity Tests/methods , Risk Assessment/methodsABSTRACT
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.
ABSTRACT
In a joint effort involving scientists from academia, industry and regulatory agencies, ECETOC's activities in Omics have led to conceptual proposals for: (1) A framework that assures data quality for reporting and inclusion of Omics data in regulatory assessments; and (2) an approach to robustly quantify these data, prior to interpretation for regulatory use. In continuation of these activities this workshop explored and identified areas of need to facilitate robust interpretation of such data in the context of deriving points of departure (POD) for risk assessment and determining an adverse change from normal variation. ECETOC was amongst the first to systematically explore the application of Omics methods, now incorporated into the group of methods known as New Approach Methodologies (NAMs), to regulatory toxicology. This support has been in the form of both projects (primarily with CEFIC/LRI) and workshops. Outputs have led to projects included in the workplan of the Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) group of the Organisation for Economic Co-operation and Development (OECD) and to the drafting of OECD Guidance Documents for Omics data reporting, with potentially more to follow on data transformation and interpretation. The current workshop was the last in a series of technical methods development workshops, with a sub-focus on the derivation of a POD from Omics data. Workshop presentations demonstrated that Omics data developed within robust frameworks for both scientific data generation and analysis can be used to derive a POD. The issue of noise in the data was discussed as an important consideration for identifying robust Omics changes and deriving a POD. Such variability or "noise" can comprise technical or biological variation within a dataset and should clearly be distinguished from homeostatic responses. Adverse outcome pathways (AOPs) were considered a useful framework on which to assemble Omics methods, and a number of case examples were presented in illustration of this point. What is apparent is that high dimension data will always be subject to varying processing pipelines and hence interpretation, depending on the context they are used in. Yet, they can provide valuable input for regulatory toxicology, with the pre-condition being robust methods for the collection and processing of data together with a comprehensive description how the data were interpreted, and conclusions reached.
Subject(s)
Adverse Outcome Pathways , Genomics , Genomics/methods , Risk Assessment , Toxicogenetics , Research DesignABSTRACT
'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to â¼100 µM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.
Subject(s)
Biological Assay , High-Throughput Screening Assays , Humans , Risk Assessment/methods , High-Throughput Screening Assays/methods , Cells, Cultured , Biological Assay/methodsABSTRACT
Screening new compounds for potential bioactivities against cellular targets is vital for drug discovery and chemical safety. Transcriptomics offers an efficient approach for assessing global gene expression changes, but interpreting chemical mechanisms from these data is often challenging. Connectivity mapping is a potential data-driven avenue for linking chemicals to mechanisms based on the observation that many biological processes are associated with unique gene expression signatures (gene signatures). However, mining the effects of a chemical on gene signatures for biological mechanisms is challenging because transcriptomic data contain thousands of noisy genes. New connectivity mapping approaches seeking to distinguish signal from noise continue to be developed, spurred by the promise of discovering chemical mechanisms, new drugs, and disease targets from burgeoning transcriptomic data. Here, we analyze these approaches in terms of diverse transcriptomic technologies, public databases, gene signatures, pattern-matching algorithms, and statistical evaluation criteria. To navigate the complexity of connectivity mapping, we propose a harmonized scheme to coherently organize and compare published workflows. We first standardize concepts underlying transcriptomic profiles and gene signatures based on various transcriptomic technologies such as microarrays, RNA-Seq, and L1000 and discuss the widely used data sources such as Gene Expression Omnibus, ArrayExpress, and MSigDB. Next, we generalize connectivity mapping as a pattern-matching task for finding similarity between a query (e.g., transcriptomic profile for new chemical) and a reference (e.g., gene signature of known target). Published pattern-matching approaches fall into two main categories: vector-based use metrics like correlation, Jaccard index, etc., and aggregation-based use parametric and nonparametric statistics (e.g., gene set enrichment analysis). The statistical methods for evaluating the performance of different approaches are described, along with comparisons reported in the literature on benchmark transcriptomic data sets. Lastly, we review connectivity mapping applications in toxicology and offer guidance on evaluating chemical-induced toxicity with concentration-response transcriptomic data. In addition to serving as a high-level guide and tutorial for understanding and implementing connectivity mapping workflows, we hope this review will stimulate new algorithms for evaluating chemical safety and drug discovery using transcriptomic data.
Subject(s)
Gene Expression Profiling , Transcriptome , Gene Expression Profiling/methods , Workflow , Databases, Factual , Drug DiscoveryABSTRACT
To date, approximately 200 chemicals have been tested in US Environmental Protection Agency (EPA) or Organization for Economic Co-operation and Development (OECD) developmental neurotoxicity (DNT) guideline studies, leaving thousands of chemicals without traditional animal information on DNT hazard potential. To address this data gap, a battery of in vitro DNT new approach methodologies (NAMs) has been proposed. Evaluation of the performance of this battery will increase the confidence in its use to determine DNT chemical hazards. One approach to evaluate DNT NAM performance is to use a set of chemicals to evaluate sensitivity and specificity. Since a list of chemicals with potential evidence of in vivo DNT has been established, this study aims to develop a curated list of "negative" chemicals for inclusion in a "DNT NAM evaluation set". A workflow, including a literature search followed by an expert-driven literature review, was used to systematically screen 39 chemicals for lack of DNT effect. Expert panel members evaluated the scientific robustness of relevant studies to inform chemical categorizations. Following review, the panel discussed each chemical and made categorical determinations of "Favorable", "Not Favorable", or "Indeterminate" reflecting acceptance, lack of suitability, or uncertainty given specific limitations and considerations, respectively. The panel determined that 10, 22, and 7 chemicals met the criteria for "Favorable", "Not Favorable", and "Indeterminate", for use as negatives in a DNT NAM evaluation set. Ultimately, this approach not only supports DNT NAM performance evaluation but also highlights challenges in identifying large numbers of negative DNT chemicals.
Subject(s)
Neurotoxicity Syndromes , Toxicity Tests , Animals , Neurotoxicity Syndromes/etiology , Research Design , Toxicity Tests/methods , United States , United States Environmental Protection AgencyABSTRACT
The United States Environmental Protection Agency has proposed a tiered testing strategy for chemical hazard evaluation based on new approach methods (NAMs). The first tier includes in vitro profiling assays applicable to many (human) cell types, such as high-throughput transcriptomics (HTTr) and high-throughput phenotypic profiling (HTPP). The goals of this study were to: (1) harmonize the seeding density of U-2 OS human osteosarcoma cells for use in both assays; (2) compare HTTr- versus HTPP-derived potency estimates for 11 mechanistically diverse chemicals; (3) identify candidate reference chemicals for monitoring assay performance in future screens; and (4) characterize the transcriptional and phenotypic changes in detail for all-trans retinoic acid (ATRA) as a model compound known for its adverse effects on osteoblast differentiation. The results of this evaluation showed that (1) HTPP conducted at low (400 cells/well) and high (3000 cells/well) seeding densities yielded comparable potency estimates and similar phenotypic profiles for the tested chemicals; (2) HTPP and HTTr resulted in comparable potency estimates for changes in cellular morphology and gene expression, respectively; (3) three test chemicals (etoposide, ATRA, dexamethasone) produced concentration-dependent effects on cellular morphology and gene expression that were consistent with known modes-of-action, demonstrating their suitability for use as reference chemicals for monitoring assay performance; and (4) ATRA produced phenotypic changes that were highly similar to other retinoic acid receptor activators (AM580, arotinoid acid) and some retinoid X receptor activators (bexarotene, methoprene acid). This phenotype was observed concurrently with autoregulation of the RARB gene. Both effects were prevented by pre-treating U-2 OS cells with pharmacological antagonists of their respective receptors. Thus, the observed phenotype could be considered characteristic of retinoic acid pathway activation in U-2 OS cells. These findings lay the groundwork for combinatorial screening of chemicals using HTTr and HTPP to generate complementary information for the first tier of a NAM-based chemical hazard evaluation strategy.
Subject(s)
Bone Neoplasms , Tretinoin , Humans , Phenotype , RNA-Seq , Receptors, Retinoic Acid/genetics , Tretinoin/pharmacology , United StatesABSTRACT
Inhalation is the most relevant route of volatile organic chemical (VOC) exposure; however, due to unique challenges posed by their chemical properties and poor solubility in aqueous solutions, in vitro chemical safety testing is predominantly performed using direct application dosing/submerged exposures. To address the difficulties in screening toxic effects of VOCs, our cell culture exposure system permits cells to be exposed to multiple concentrations at air-liquid interface (ALI) in a 24-well format. ALI exposure methods permit direct chemical-to-cell interaction with the test article at physiological conditions. In the present study, BEAS-2B and primary normal human bronchial epithelial cells (pHBEC) are used to assess gene expression, cytotoxicity, and cell viability responses to a variety of volatile chemicals including acrolein, formaldehyde, 1,3-butadiene, acetaldehyde, 1-bromopropane, carbon tetrachloride, dichloromethane, and trichloroethylene. BEAS-2B cells were exposed to all the test agents, whereas pHBECs were only exposed to the latter 4 listed above. The VOC concentrations tested elicited only slight cell viability changes in both cell types. Gene expression changes were analyzed using benchmark dose (BMD) modeling. The BMD for the most sensitive gene set was within one order of magnitude of the threshold-limit value reported by the American Conference of Governmental Industrial Hygienists, and the most sensitive gene sets impacted by exposure correlate to known adverse health effects recorded in epidemiologic and in vivo exposure studies. Overall, our study outlines a novel in vitro approach for evaluating molecular-based points-of-departure in human airway epithelial cell exposure to volatile chemicals.
Subject(s)
Air Pollutants , Volatile Organic Compounds , Acetaldehyde , Benchmarking , Formaldehyde , Humans , Volatile Organic Compounds/analysisABSTRACT
In vivo developmental neurotoxicity (DNT) testing is resource intensive and lacks information on cellular processes affected by chemicals. To address this, DNT new approach methodologies (NAMs) are being evaluated, including: the microelectrode array neuronal network formation assay; and high-content imaging to evaluate proliferation, apoptosis, neurite outgrowth, and synaptogenesis. This work addresses 3 hypotheses: (1) a broad screening battery provides a sensitive marker of DNT bioactivity; (2) selective bioactivity (occurring at noncytotoxic concentrations) may indicate functional processes disrupted; and, (3) a subset of endpoints may optimally classify chemicals with in vivo evidence for DNT. The dataset was comprised of 92 chemicals screened in all 57 assay endpoints sourced from publicly available data, including a set of DNT NAM evaluation chemicals with putative positives (53) and negatives (13). The DNT NAM battery provides a sensitive marker of DNT bioactivity, particularly in cytotoxicity and network connectivity parameters. Hierarchical clustering suggested potency (including cytotoxicity) was important for classifying positive chemicals with high sensitivity (93%) but failed to distinguish patterns of disrupted functional processes. In contrast, clustering of selective values revealed informative patterns of differential activity but demonstrated lower sensitivity (74%). The false negatives were associated with several limitations, such as the maximal concentration tested or gaps in the biology captured by the current battery. This work demonstrates that this multi-dimensional assay suite provides a sensitive biomarker for DNT bioactivity, with selective activity providing possible insight into specific functional processes affected by chemical exposure and a basis for further research.
Subject(s)
Neurotoxicity Syndromes , Toxicity Tests , Humans , Neurogenesis , Neuronal Outgrowth , Neurons , Neurotoxicity Syndromes/etiology , Toxicity Tests/methodsABSTRACT
Omics methodologies are widely used in toxicological research to understand modes and mechanisms of toxicity. Increasingly, these methodologies are being applied to questions of regulatory interest such as molecular point-of-departure derivation and chemical grouping/read-across. Despite its value, widespread regulatory acceptance of omics data has not yet occurred. Barriers to the routine application of omics data in regulatory decision making have been: 1) lack of transparency for data processing methods used to convert raw data into an interpretable list of observations; and 2) lack of standardization in reporting to ensure that omics data, associated metadata and the methodologies used to generate results are available for review by stakeholders, including regulators. Thus, in 2017, the Organisation for Economic Co-operation and Development (OECD) Extended Advisory Group on Molecular Screening and Toxicogenomics (EAGMST) launched a project to develop guidance for the reporting of omics data aimed at fostering further regulatory use. Here, we report on the ongoing development of the first formal reporting framework describing the processing and analysis of both transcriptomic and metabolomic data for regulatory toxicology. We introduce the modular structure, content, harmonization and strategy for trialling this reporting framework prior to its publication by the OECD.
Subject(s)
Metabolomics/standards , Organisation for Economic Co-Operation and Development/standards , Toxicogenetics/standards , Toxicology/standards , Transcriptome/physiology , Documentation/standards , HumansABSTRACT
Using in vitro data to estimate point of departure (POD) values is an essential component of new approach methodologies (NAMs)-based chemical risk assessments. In this case study, we evaluated a NAM for hepatotoxicity based on rat primary hepatocytes, high-content imaging (HCI), and toxicokinetic modeling. First, we treated rat primary hepatocytes with 10 concentrations (0.2-100 µM) of 51 chemicals that produced hepatotoxicity in repeat-dose subchronic and chronic exposures. Second, we used HCI to measure endoplasmic reticulum stress, mitochondrial function, lysosomal mass, steatosis, apoptosis, DNA texture, nuclear size, and cell number at 24, 48, and 72 h and calculated concentrations at 50% maximal activity (AC50). Third, we estimated administered equivalent doses (AEDs) from AC50 values using toxicokinetic modeling. AEDs using physiologically based toxicokinetic models were 4.1-fold (SD 6.3) and 8.1-fold (SD 15.5) lower than subchronic and chronic lowest observed adverse effect levels (LOAELs), respectively. In contrast, AEDs from ToxCast and Tox21 assays were 89.8-fold (SD 149.5) and 168-fold (SD 323.7) lower than subchronic and chronic LOAELs. Individual HCI endpoints also estimated AEDs for specific hepatic lesions that were lower than in vivo PODs. Lastly, AEDs were similar for different in vitro exposure durations, but steady-state toxicokinetic models produced 7.6-fold lower estimates than dynamic physiologically based ones. Our findings suggest that NAMs from diverse cell types provide conservative estimates of PODs. In contrast, NAMs based on the same species and cell type as the adverse outcome may produce estimates closer to the traditional in vivo PODs.