Metabolomics Simultaneously Derives Benchmark Dose Estimates and Discovers Metabolic Biotransformations in a Rat Bioassay.

Benchmark dose (BMD) modeling estimates the dose of a chemical that causes a perturbation from baseline. Transcriptional BMDs have been shown to be relatively consistent with apical end point BMDs, opening the door to using molecular BMDs to derive human health-based guidance values for chemical exposure. Metabolomics measures the responses of small-molecule endogenous metabolites to chemical exposure, complementing transcriptomics by characterizing downstream molecular phenotypes that are more closely associated with apical end points. The aim of this study was to apply BMD modeling to in vivo metabolomics data, to compare metabolic BMDs to both transcriptional and apical end point BMDs. This builds upon our previous application of transcriptomics and BMD modeling to a 5-day rat study of triphenyl phosphate (TPhP), applying metabolomics to the same archived tissues. Specifically, liver from rats exposed to five doses of TPhP was investigated using liquid chromatography-mass spectrometry and 1H nuclear magnetic resonance spectroscopy-based metabolomics. Following the application of BMDExpress2 software, 2903 endogenous metabolic features yielded viable dose-response models, confirming a perturbation to the liver metabolome. Metabolic BMD estimates were similarly sensitive to transcriptional BMDs, and more sensitive than both clinical chemistry and apical end point BMDs. Pathway analysis of the multiomics data sets revealed a major effect of TPhP exposure on cholesterol (and downstream) pathways, consistent with clinical chemistry measurements. Additionally, the transcriptomics data indicated that TPhP activated xenobiotic metabolism pathways, which was confirmed by using the underexploited capability of metabolomics to detect xenobiotic-related compounds. Eleven biotransformation products of TPhP were discovered, and their levels were highly correlated with multiple xenobiotic metabolism genes. This work provides a case study showing how metabolomics and transcriptomics can estimate mechanistically anchored points-of-departure. Furthermore, the study demonstrates how metabolomics can also discover biotransformation products, which could be of value within a regulatory setting, for example, as an enhancement of OECD Test Guideline 417 (toxicokinetics).

Multi-omics bioactivity profile-based chemical grouping and read-across: a case study with Daphnia magna and azo dyes.

Grouping/read-across is widely used for predicting the toxicity of data-poor target substance(s) using data-rich source substance(s). While the chemical industry and the regulators recognise its benefits, registration dossiers are often rejected due to weak analogue/category justifications based largely on the structural similarity of source and target substances. Here we demonstrate how multi-omics measurements can improve confidence in grouping via a statistical assessment of the similarity of molecular effects. Six azo dyes provided a pool of potential source substances to predict long-term toxicity to aquatic invertebrates (Daphnia magna) for the dye Disperse Yellow 3 (DY3) as the target substance. First, we assessed the structural similarities of the dyes, generating a grouping hypothesis with DY3 and two Sudan dyes within one group. Daphnia magna were exposed acutely to equi-effective doses of all seven dyes (each at 3 doses and 3 time points), transcriptomics and metabolomics data were generated from 760 samples. Multi-omics bioactivity profile-based grouping uniquely revealed that Sudan 1 (S1) is the most suitable analogue for read-across to DY3. Mapping ToxPrint structural fingerprints of the dyes onto the bioactivity profile-based grouping indicated an aromatic alcohol moiety could be responsible for this bioactivity similarity. The long-term reproductive toxicity to aquatic invertebrates of DY3 was predicted from S1 (21-day NOEC, 40 µg/L). This prediction was confirmed experimentally by measuring the toxicity of DY3 in D. magna. While limitations of this 'omics approach are identified, the study illustrates an effective statistical approach for building chemical groups.

Energy Starvation in Daphnia magna from Exposure to a Lithium Cobalt Oxide Nanomaterial.

Growing evidence across organisms points to altered energy metabolism as an adverse outcome of metal oxide nanomaterial toxicity, with a mechanism of toxicity potentially related to the redox chemistry of processes involved in energy production. Despite this evidence, the significance of this mechanism has gone unrecognized in nanotoxicology due to the field's focus on oxidative stress as a universal─but nonspecific─nanotoxicity mechanism. To further explore metabolic impacts, we determined lithium cobalt oxide's (LCO's) effects on these pathways in the model organism Daphnia magna through global gene-expression analysis using RNA-Seq and untargeted metabolomics by direct-injection mass spectrometry. Our results show that a sublethal 1 mg/L 48 h exposure of D. magna to LCO nanosheets causes significant impacts on metabolic pathways versus untreated controls, while exposure to ions released over 48 h does not. Specifically, transcriptomic analysis using DAVID indicated significant enrichment (Benjamini-adjusted p ≤0.0.5) in LCO-exposed animals for changes in pathways involved in the cellular response to starvation (25 genes), mitochondrial function (70 genes), ATP-binding (70 genes), oxidative phosphorylation (53 genes), NADH dehydrogenase activity (12 genes), and protein biosynthesis (40 genes). Metabolomic analysis using MetaboAnalyst indicated significant enrichment (γ-adjusted p <0.1) for changes in amino acid metabolism (19 metabolites) and starch, sucrose, and galactose metabolism (7 metabolites). Overlap of significantly impacted pathways by RNA-Seq and metabolomics suggests amino acid breakdown and increased sugar import for energy production. Results indicate that LCO-exposed Daphnia respond to energy starvation by altering metabolic pathways, both at the gene expression and metabolite levels. These results support altered energy production as a sensitive nanotoxicity adverse outcome for LCO exposure and suggest negative impacts on energy metabolism as an important avenue for future studies of nanotoxicity, including for other biological systems and for metal oxide nanomaterials more broadly.

Knowledge-Driven Approaches to Create the MTox700+ Metabolite Panel for Predicting Toxicity.

Endogenous metabolite levels describe the molecular phenotype that is most downstream from chemical exposure. Consequently, quantitative changes in metabolite levels have the potential to predict mode-of-action and adversity, with regulatory toxicology predicated on the latter. However, toxicity-related metabolic biomarker resources remain highly fragmented and incomplete. Although development of the S1500+ gene biomarker panel has accelerated the application of transcriptomics to toxicology, a similar initiative for metabolic biomarkers is lacking. Our aim was to define a publicly available metabolic biomarker panel, equivalent to S1500+, capable of predicting pathway perturbations and/or adverse outcomes. We conducted a systematic review of multiple toxicological resources, yielding 189 proposed metabolic biomarkers from existing assays (BASF, Bowes-44, and Tox21), 342 biomarkers from databases (Adverse Outcome Pathway Wiki, Comparative Toxicogenomics Database, QIAGEN Ingenuity Pathway Analysis, and Toxin and Toxin-Target Database), and 435 biomarkers from the literature. Evidence mapping across all 8 resources generated a panel of 722 metabolic biomarkers for toxicology (MTox700+), of which 462 (64%) are associated with molecular pathways and 575 (80%) with adverse outcomes. Comparing MTox700+ and S1500+ revealed that 418 (58%) metabolic biomarkers associate with pathways shared across both panels, with further metabolites mapping to unique pathways. Metabolite reference standards are commercially available for 646 (90%) of the panel metabolites, and assays exist for 578 (80%) of these biomarkers. This study has generated a publicly available metabolic biomarker panel for toxicology, which through its future laboratory deployment, is intended to help build foundational knowledge to support the generation of molecular mechanistic data for chemical hazard assessment.

Toxicity by descent: A comparative approach for chemical hazard assessment.

Toxicology is traditionally divided between human and eco-toxicology. In the shared pursuit of environmental health, this separation does not account for discoveries made in the comparative studies of animal genomes. Here, we provide evidence on the feasibility of understanding the health impact of chemicals on all animals, including ecological keystone species and humans, based on a significant number of conserved genes and their functional associations to health-related outcomes across much of animal diversity. We test four conditions to understand the value of comparative genomics data to inform mechanism-based human and environmental hazard assessment: (1) genes that are most fundamental for health evolved early during animal evolution; (2) the molecular functions of pathways are better conserved among distantly related species than the individual genes that are members of these pathways; (3) the most conserved pathways among animals are those that cause adverse health outcomes when disrupted; (4) gene sets that serve as molecular signatures of biological processes or disease-states are largely enriched by evolutionarily conserved genes across the animal phylogeny. The concept of homology is applied in a comparative analysis of gene families and pathways among invertebrate and vertebrate species compared with humans. Results show that over 70% of gene families associated with disease are shared among the greatest variety of animal species through evolution. Pathway conservation between invertebrates and humans is based on the degree of conservation within vertebrates and the number of interacting genes within the human network. Human gene sets that already serve as biomarkers are enriched by evolutionarily conserved genes across the animal phylogeny. By implementing a comparative method for chemical hazard assessment, human and eco-toxicology converge towards a more holistic and mechanistic understanding of toxicity disrupting biological processes that are important for health and shared among animals (including humans).