ToxicoDB: an integrated database to mine and visualize large-scale toxicogenomic datasets

A strategy to detect metabolic changes induced by exposure to chemicals from large sets of condition-specific metabolic models computed with enumeration techniques

BACKGROUND

The growing abundance of in vitro omics data, coupled with the necessity to reduce animal testing in the safety assessment of chemical compounds and even eliminate it in the evaluation of cosmetics, highlights the need for adequate computational methodologies. Data from omics technologies allow the exploration of a wide range of biological processes, therefore providing a better understanding of mechanisms of action (MoA) related to chemical exposure in biological systems. However, the analysis of these large datasets remains difficult due to the complexity of modulations spanning multiple biological processes.

RESULTS

To address this, we propose a strategy to reduce information overload by computing, based on transcriptomics data, a comprehensive metabolic sub-network reflecting the metabolic impact of a chemical. The proposed strategy integrates transcriptomic data to a genome scale metabolic network through enumeration of condition-specific metabolic models hence translating transcriptomics data into reaction activity probabilities. Based on these results, a graph algorithm is applied to retrieve user readable sub-networks reflecting the possible metabolic MoA (mMoA) of chemicals. This strategy has been implemented as a three-step workflow. The first step consists in building cell condition-specific models reflecting the metabolic impact of each exposure condition while taking into account the diversity of possible optimal solutions with a partial enumeration algorithm. In a second step, we address the challenge of analyzing thousands of enumerated condition-specific networks by computing differentially activated reactions (DARs) between the two sets of enumerated possible condition-specific models. Finally, in the third step, DARs are grouped into clusters of functionally interconnected metabolic reactions, representing possible mMoA, using the distance-based clustering and subnetwork extraction method. The first part of the workflow was exemplified on eight molecules selected for their known human hepatotoxic outcomes associated with specific MoAs well described in the literature and for which we retrieved primary human hepatocytes transcriptomic data in Open TG-GATEs. Then, we further applied this strategy to more precisely model and visualize associated mMoA for two of these eight molecules (amiodarone and valproic acid). The approach proved to go beyond gene-based analysis by identifying mMoA when few genes are significantly differentially expressed (2 differentially expressed genes (DEGs) for amiodarone), bringing additional information from the network topology, or when very large number of genes were differentially expressed (5709 DEGs for valproic acid). In both cases, the results of our strategy well fitted evidence from the literature regarding known MoA. Beyond these confirmations, the workflow highlighted potential other unexplored mMoA.

CONCLUSION

The proposed strategy allows toxicology experts to decipher which part of cellular metabolism is expected to be affected by the exposition to a given chemical. The approach originality resides in the combination of different metabolic modelling approaches (constraint based and graph modelling). The application to two model molecules shows the strong potential of the approach for interpretation and visual mining of complex omics in vitro data. The presented strategy is freely available as a python module ( https://pypi.org/project/manamodeller/ ) and jupyter notebooks ( https://github.com/LouisonF/MANA ).

Application of transfer learning to predict drug-induced human in vivo gene expression changes using rat in vitro and in vivo data

The liver is the primary site for the metabolism and detoxification of many compounds, including pharmaceuticals. Consequently, it is also the primary location for many adverse reactions. As the liver is not readily accessible for sampling in humans; rodent or cell line models are often used to evaluate potential toxic effects of a novel compound or candidate drug. However, relating the results of animal and in vitro studies to relevant clinical outcomes for the human in vivo situation still proves challenging. In this study, we incorporate principles of transfer learning within a deep artificial neural network allowing us to leverage the relative abundance of rat in vitro and in vivo exposure data from the Open TG-GATEs data set to train a model to predict the expected pattern of human in vivo gene expression following an exposure given measured human in vitro gene expression. We show that domain adaptation has been successfully achieved, with the rat and human in vitro data no longer being separable in the common latent space generated by the network. The network produces physiologically plausible predictions of human in vivo gene expression pattern following an exposure to a previously unseen compound. Moreover, we show the integration of the human in vitro data in the training of the domain adaptation network significantly improves the temporal accuracy of the predicted rat in vivo gene expression pattern following an exposure to a previously unseen compound. In this way, we demonstrate the improvements in prediction accuracy that can be achieved by combining data from distinct domains.

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.

Metabolic modeling of sex-specific tissue predicts mechanisms of differences in toxicological responses

Male subjects in animal and human studies are disproportionately used for toxicological testing. This discrepancy is evidenced in clinical medicine where females are more likely than males to experience liver-related adverse events in response to xenobiotics. While previous work has shown gene expression differences between the sexes, there is a lack of systems-level approaches to understand the direct clinical impact effect of these differences. Here, we integrate gene expression data with metabolic network models to characterize the impact of transcriptional changes of metabolic genes in the context of sex differences and drug treatment. We used Tasks Inferred from Differential Expression (TIDEs), a reaction-centric approach to analyzing differences in gene expression, to discover that androgen, ether lipid, glucocorticoid, tryptophan, and xenobiotic metabolism have more activity in the male liver, and serotonin, melatonin, pentose, glucuronate, and vitamin A metabolism have more activity in the female liver. When TIDEs is used to compare expression differences in treated and untreated hepatocytes, we see little response in those sex-altered subsystems, and the largest differences are in subsystems related to lipid metabolism. Finally, using sex-specific transcriptomic data, we create individual and averaged male and female liver models and find differences in the import of bile acids and salts. This result suggests that the sexually dimorphic behavior of the liver may be caused by differences in enterohepatic recirculation, and we suggest an investigation into sex-specific microbiome composition as an avenue of further research.

Author Summary

Male-bias in clinical testing of drugs has led to a disproportionate number of hepatotoxic events in women. Previous works use gene-by-gene differences in biological sex to explain this discrepancy, but there is little focus on the systematic interactions of these differences. To this end, we use a combination of gene expression data and metabolic modeling to compare metabolic activity between the male and female liver and treated and untreated hepatocytes. We find several subsystems with differential activity in each sex; however, when comparing these subsystems with those pathways altered by hepatotoxic agents, we find little overlap. To explore these differences on a reaction-by-reaction basis, we use the same sex-specific transcriptomic data to contextualize the previously published Human1 human cell metabolic model. In these models we find a difference in flux for the import of bile acids and salts, suggesting a potential difference in enterohepatic circulation. These findings can help guide future drug design, toxicological testing, and sex-specific research to better account for the entire human population.

TOXRIC: a comprehensive database of toxicological data and benchmarks

The toxic effects of compounds on environment, humans, and other organisms have been a major focus of many research areas, including drug discovery and ecological research. Identifying the potential toxicity in the early stage of compound/drug discovery is critical. The rapid development of computational methods for evaluating various toxicity categories has increased the need for comprehensive and system-level collection of toxicological data, associated attributes, and benchmarks. To contribute toward this goal, we proposed TOXRIC (https://toxric.bioinforai.tech/), a database with comprehensive toxicological data, standardized attribute data, practical benchmarks, informative visualization of molecular representations, and an intuitive function interface. The data stored in TOXRIC contains 113 372 compounds, 13 toxicity categories, 1474 toxicity endpoints covering in vivo/in vitro endpoints and 39 feature types, covering structural, target, transcriptome, metabolic data, and other descriptors. All the curated datasets of endpoints and features can be retrieved, downloaded and directly used as output or input to Machine Learning (ML)-based prediction models. In addition to serving as a data repository, TOXRIC also provides visualization of benchmarks and molecular representations for all endpoint datasets. Based on these results, researchers can better understand and select optimal feature types, molecular representations, and baseline algorithms for each endpoint prediction task. We believe that the rich information on compound toxicology, ML-ready datasets, benchmarks and molecular representation distribution can greatly facilitate toxicological investigations, interpretation of toxicological mechanisms, compound/drug discovery and the development of computational methods.

Machine learning and artificial intelligence in toxicological sciences

Machine learning and artificial intelligence approaches have revolutionized multiple disciplines, including toxicology. This review summarizes representative recent applications of machine learning and artificial intelligence approaches in different areas of toxicology, including physiologically based pharmacokinetic (PBPK) modeling, quantitative structure-activity relationship modeling for toxicity prediction, adverse outcome pathway analysis, high-throughput screening, toxicogenomics, big data and toxicological databases. By leveraging machine learning and artificial intelligence approaches, now it is possible to develop PBPK models for hundreds of chemicals efficiently, to create in silico models to predict toxicity for a large number of chemicals with similar accuracies compared to in vivo animal experiments, and to analyze a large amount of different types of data (toxicogenomics, high-content image data, etc.) to generate new insights into toxicity mechanisms rapidly, which was impossible by manual approaches in the past. To continue advancing the field of toxicological sciences, several challenges should be considered: (1) not all machine learning models are equally useful for a particular type of toxicology data, and thus it is important to test different methods to determine the optimal approach; (2) current toxicity prediction is mainly on bioactivity classification (yes/no), so additional studies are needed to predict the intensity of effect or dose-response relationship; (3) as more data become available, it is crucial to perform rigorous data quality check and develop infrastructure to store, share, analyze, evaluate, and manage big data; and (4) it is important to convert machine learning models to user-friendly interfaces to facilitate their applications by both computational and bench scientists.

Orchestrating and sharing large multimodal data for transparent and reproducible research

Reproducibility is essential to open science, as there is limited relevance for findings that can not be reproduced by independent research groups, regardless of its validity. It is therefore crucial for scientists to describe their experiments in sufficient detail so they can be reproduced, scrutinized, challenged, and built upon. However, the intrinsic complexity and continuous growth of biomedical data makes it increasingly difficult to process, analyze, and share with the community in a FAIR (findable, accessible, interoperable, and reusable) manner. To overcome these issues, we created a cloud-based platform called ORCESTRA ( orcestra.ca ), which provides a flexible framework for the reproducible processing of multimodal biomedical data. It enables processing of clinical, genomic and perturbation profiles of cancer samples through automated processing pipelines that are user-customizable. ORCESTRA creates integrated and fully documented data objects with persistent identifiers (DOI) and manages multiple dataset versions, which can be shared for future studies.

Mexican Ganoderma Lucidum Extracts Decrease Lipogenesis Modulating Transcriptional Metabolic Networks and Gut Microbiota in C57BL/6 Mice Fed with a High-Cholesterol Diet

Prevention of hyperlipidemia and associated diseases is a health priority. Natural products, such as the medicinal mushroom Ganoderma lucidum (Gl), have demonstrated hypocholesterolemic, prebiotic and antidiabetic properties. However, the underlying transcriptomic mechanisms by which Gl exerts bioactivities are not completely understood. We report a comprehensive hepatic and renal transcriptome profiling of C57BL/6 mice under the consumption of a high-cholesterol diet and two standardized Gl extracts obtained from basidiocarps cultivated on conventional substrate (Gl-1) or substrate containing acetylsalicylic acid (ASA; Gl-2). We showed that Gl extracts modulate relevant metabolic pathways involving the restriction of lipid biosynthesis and the enrichment of lipid degradation and secretion. The Gl-2 extract exerts a major modulation over gene expression programs showing the highest similarity with simvastatin druggable-target-genes and these are enriched more in processes related to human obesity alterations in the liver. We further show a subset of Gl-modulated genes correlated with Lactobacillus enrichment and the reduction of circulating cholesterol-derived fats. Moreover, Gl extracts induce a significant decrease of macrophage lipid storage, which occurs concomitantly with the down-modulation of Fasn and Elovl6. Collectively, this evidence suggests a new link between Gl hypocholesterolemic and prebiotic activity, revealing thereby that standardized Mexican Gl extracts are a novel transcriptome modulator to prevent metabolic disorders associated with hypercholesterolemia.