A Network Toxicology Approach for Mechanistic Modelling of Nanomaterial Hazard and Adverse Outcomes.

Hazard assessment is the first step in evaluating the potential adverse effects of chemicals. Traditionally, toxicological assessment has focused on the exposure, overlooking the impact of the exposed system on the observed toxicity. However, systems toxicology emphasizes how system properties significantly contribute to the observed response. Hence, systems theory states that interactions store more information than individual elements, leading to the adoption of network based models to represent complex systems in many fields of life sciences. Here, they develop a network-based approach to characterize toxicological responses in the context of a biological system, inferring biological system specific networks. They directly link molecular alterations to the adverse outcome pathway (AOP) framework, establishing direct connections between omics data and toxicologically relevant phenotypic events. They apply this framework to a dataset including 31 engineered nanomaterials with different physicochemical properties in two different in vitro and one in vivo models and demonstrate how the biological system is the driving force of the observed response. This work highlights the potential of network-based methods to significantly improve their understanding of toxicological mechanisms from a systems biology perspective and provides relevant considerations and future data-driven approaches for the hazard assessment of nanomaterials and other advanced materials.

Advancing chemical safety assessment through an omics-based characterization of the test system-chemical interaction.

Assessing chemical safety is essential to evaluate the potential risks of chemical exposure to human health and the environment. Traditional methods relying on animal testing are being replaced by 3R (reduction, refinement, and replacement) principle-based alternatives, mainly depending on in vitro test methods and the Adverse Outcome Pathway framework. However, these approaches often focus on the properties of the compound, missing the broader chemical-biological interaction perspective. Currently, the lack of comprehensive molecular characterization of the in vitro test system results in limited real-world representation and contextualization of the toxicological effect under study. Leveraging omics data strengthens the understanding of the responses of different biological systems, emphasizing holistic chemical-biological interactions when developing in vitro methods. Here, we discuss the relevance of meticulous test system characterization on two safety assessment relevant scenarios and how omics-based, data-driven approaches can improve the future generation of alternative methods.

Expanding adverse outcome pathways towards one health models for nanosafety.

The ever-growing production of nano-enabled products has generated the need for dedicated risk assessment strategies that ensure safety for humans and the environment. Transdisciplinary approaches are needed to support the development of new technologies while respecting environmental limits, as also highlighted by the EU Green Deal Chemicals Strategy for Sustainability and its safe and sustainable by design (SSbD) framework. The One Health concept offers a holistic multiscale approach for the assessment of nanosafety. However, toxicology is not yet capable of explaining the interaction between chemicals and biological systems at the multiscale level and in the context of the One Health framework. Furthermore, there is a disconnect between chemical safety assessment, epidemiology, and other fields of biology that, if unified, would enable the adoption of the One Health model. The development of mechanistic toxicology and the generation of omics data has provided important biological knowledge of the response of individual biological systems to nanomaterials (NMs). On the other hand, epigenetic data have the potential to inform on interspecies mechanisms of adaptation. These data types, however, need to be linked to concepts that support their intuitive interpretation. Adverse Outcome Pathways (AOPs) represent an evolving framework to anchor existing knowledge to chemical risk assessment. In this perspective, we discuss the possibility of integrating multi-level toxicogenomics data, including toxicoepigenetic insights, into the AOP framework. We anticipate that this new direction of toxicogenomics can support the development of One Health models applicable to groups of chemicals and to multiple species in the tree of life.

The in vitro immunomodulatory effect of multi-walled carbon nanotubes by multilayer analysis.

The study of multi-walled carbon nanotube (MWCNT) induced immunotoxicity is crucial for determining hazards posed to human health. MWCNT exposure most commonly occurs via the airways, where macrophages are first line responders. Here we exploit an in vitro assay, measuring dose-dependent secretion of a wide panel of cytokines, as a measure of immunotoxicity following the non-lethal, multi-dose exposure (IC5, IC10 and IC20) to 7 MWCNTs with different intrinsic properties. We find that a tangled structure, and small aspect ratio are key properties predicting MWCNT induced immunotoxicity, mediated predominantly by IL1B cytokine secretion. To assess the mechanism of action giving rise to MWCNT immunotoxicity, transcriptomics analysis was linked to cytokine secretion in a multilayer model established through correlation analysis across exposure concentrations. This reinforced the finding that tangled MWCNTs have greater immunomodulatory potency, displaying enrichment of immune system, signal transduction and pattern recognition associated pathways. Together our results further elucidate how structure, length and aspect ratio, critical intrinsic properties of MWCNTs, are tied to immunotoxicity.

Identification of a Gene Signature Predicting (Nano)Particle-Induced Adverse Lung Outcome in Rats.

Nanoparticles are extensively used in industrial products or as food additives. However, despite their contribution to improving our quality of life, concerns have been raised regarding their potential impact on occupational and public health. To speed up research assessing nanoparticle-related hazards, this study was undertaken to identify early markers of harmful effects on the lungs. Female Sprague Dawley rats were either exposed to crystalline silica DQ-12 with instillation, or to titanium dioxide P25, carbon black Printex-90, or multi-walled carbon nanotube Mitsui-7 with nose-only inhalation. Tissues were collected at three post-exposure time points to assess short- and long-term effects. All particles induced lung inflammation. Histopathological and biochemical analyses revealed phospholipid accumulation, lipoproteinosis, and interstitial thickening with collagen deposition after exposure to DQ-12. Exposure to the highest dose of Printex-90 and Mitsui-7, but not P25, induced some phospholipid accumulation. Comparable histopathological changes were observed following exposure to P25, Printex-90, and Mitsui-7. Comparison of overall gene expression profiles identified 15 potential early markers of adverse lung outcomes induced by spherical particles. With Mitsui-7, a distinct gene expression signature was observed, suggesting that carbon nanotubes trigger different toxicity mechanisms to spherical particles.

ESPERANTO: a GLP-field sEmi-SuPERvised toxicogenomics metadAta curatioN TOol.

SUMMARY: Biological data repositories are an invaluable source of publicly available research evidence. Unfortunately, the lack of convergence of the scientific community on a common metadata annotation strategy has resulted in large amounts of data with low FAIRness (Findable, Accessible, Interoperable and Reusable). The possibility of generating high-quality insights from their integration relies on data curation, which is typically an error-prone process while also being expensive in terms of time and human labour. Here, we present ESPERANTO, an innovative framework that enables a standardized semi-supervised harmonization and integration of toxicogenomics metadata and increases their FAIRness in a Good Laboratory Practice-compliant fashion. The harmonization across metadata is guaranteed with the definition of an ad hoc vocabulary. The tool interface is designed to support the user in metadata harmonization in a user-friendly manner, regardless of the background and the type of expertise. AVAILABILITY AND IMPLEMENTATION: ESPERANTO and its user manual are freely available for academic purposes at https://github.com/fhaive/esperanto. The input and the results showcased in Supplementary File S1 are available at the same link.

A curated gene and biological system annotation of adverse outcome pathways related to human health.

Adverse outcome pathways (AOPs) are emerging as a central framework in modern toxicology and other fields in biomedicine. They serve as an extension of pathway-based concepts by depicting biological mechanisms as causally linked sequences of key events (KEs) from a molecular initiating event (MIE) to an adverse outcome. AOPs guide the use and development of new approach methodologies (NAMs) aimed at reducing animal experimentation. While AOPs model the systemic mechanisms at various levels of biological organisation, toxicogenomics provides the means to study the molecular mechanisms of chemical exposures. Systematic integration of these two concepts would improve the application of AOP-based knowledge while also supporting the interpretation of complex omics data. Hence, we established this link through rigorous curation of molecular annotations for the KEs of human relevant AOPs. We further expanded and consolidated the annotations of the biological context of KEs. These curated annotations pave the way to embed AOPs in molecular data interpretation, facilitating the emergence of new knowledge in biomedicine.

Toxicogenomics Data for Chemical Safety Assessment and Development of New Approach Methodologies: An Adverse Outcome Pathway-Based Approach.

Mechanistic toxicology provides a powerful approach to inform on the safety of chemicals and the development of safe-by-design compounds. Although toxicogenomics supports mechanistic evaluation of chemical exposures, its implementation into the regulatory framework is hindered by uncertainties in the analysis and interpretation of such data. The use of mechanistic evidence through the adverse outcome pathway (AOP) concept is promoted for the development of new approach methodologies (NAMs) that can reduce animal experimentation. However, to unleash the full potential of AOPs and build confidence into toxicogenomics, robust associations between AOPs and patterns of molecular alteration need to be established. Systematic curation of molecular events to AOPs will create the much-needed link between toxicogenomics and systemic mechanisms depicted by the AOPs. This, in turn, will introduce novel ways of benefitting from the AOPs, including predictive models and targeted assays, while also reducing the need for multiple testing strategies. Hence, a multi-step strategy to annotate AOPs is developed, and the resulting associations are applied to successfully highlight relevant adverse outcomes for chemical exposures with strong in vitro and in vivo convergence, supporting chemical grouping and other data-driven approaches. Finally, a panel of AOP-derived in vitro biomarkers for pulmonary fibrosis (PF) is identified and experimentally validated.

Biomarkers of nanomaterials hazard from multi-layer data.

There is an urgent need to apply effective, data-driven approaches to reliably predict engineered nanomaterial (ENM) toxicity. Here we introduce a predictive computational framework based on the molecular and phenotypic effects of a large panel of ENMs across multiple in vitro and in vivo models. Our methodology allows for the grouping of ENMs based on multi-omics approaches combined with robust toxicity tests. Importantly, we identify mRNA-based toxicity markers and extensively replicate them in multiple independent datasets. We find that models based on combinations of omics-derived features and material intrinsic properties display significantly improved predictive accuracy as compared to physicochemical properties alone.

Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration.

The molecular effects of exposures to engineered nanomaterials (ENMs) are still largely unknown. In classical inhalation toxicology, cell composition of bronchoalveolar lavage (BAL) is a toxicity indicator at the lung tissue level that can aid in interpreting pulmonary histological changes. Toxicogenomic approaches help characterize the mechanism of action (MOA) of ENMs by investigating the differentially expressed genes (DEG). However, dissecting which molecular mechanisms and events are directly induced by the exposure is not straightforward. It is now generally accepted that direct effects follow a monotonic dose-dependent pattern. Here, we applied an integrated modeling approach to study the MOA of four ENMs by retrieving the DEGs that also show a dynamic dose-dependent profile (dddtMOA). We further combined the information of the dddtMOA with the dose dependency of four immune cell populations derived from BAL counts. The dddtMOA analysis highlighted the specific adaptation pattern to each ENM. Furthermore, it revealed the distinct effect of the ENM physicochemical properties on the induced immune response. Finally, we report three genes dose-dependent in all the exposures and correlated with immune deregulation in the lung. The characterization of dddtMOA for ENM exposures, both for apical endpoints and molecular responses, can further promote toxicogenomic approaches in a regulatory context.

Nextcast: A software suite to analyse and model toxicogenomics data.

The recent advancements in toxicogenomics have led to the availability of large omics data sets, representing the starting point for studying the exposure mechanism of action and identifying candidate biomarkers for toxicity prediction. The current lack of standard methods in data generation and analysis hampers the full exploitation of toxicogenomics-based evidence in regulatory risk assessment. Moreover, the pipelines for the preprocessing and downstream analyses of toxicogenomic data sets can be quite challenging to implement. During the years, we have developed a number of software packages to address specific questions related to multiple steps of toxicogenomics data analysis and modelling. In this review we present the Nextcast software collection and discuss how its individual tools can be combined into efficient pipelines to answer specific biological questions. Nextcast components are of great support to the scientific community for analysing and interpreting large data sets for the toxicity evaluation of compounds in an unbiased, straightforward, and reliable manner. The Nextcast software suite is available at: ( https://github.com/fhaive/nextcast).

Microarray Data Preprocessing: From Experimental Design to Differential Analysis.

DNA microarray data preprocessing is of utmost importance in the analytical path starting from the experimental design and leading to a reliable biological interpretation. In fact, when all relevant aspects regarding the experimental plan have been considered, the following steps from data quality check to differential analysis will lead to robust, trustworthy results. In this chapter, all the relevant aspects and considerations about microarray preprocessing will be discussed. Preprocessing steps are organized in an orderly manner, from experimental design to quality check and batch effect removal, including the most common visualization methods. Furthermore, we will discuss data representation and differential testing methods with a focus on the most common microarray technologies, such as gene expression and DNA methylation.

Manually curated transcriptomics data collection for toxicogenomic assessment of engineered nanomaterials.

Toxicogenomics (TGx) approaches are increasingly applied to gain insight into the possible toxicity mechanisms of engineered nanomaterials (ENMs). Omics data can be valuable to elucidate the mechanism of action of chemicals and to develop predictive models in toxicology. While vast amounts of transcriptomics data from ENM exposures have already been accumulated, a unified, easily accessible and reusable collection of transcriptomics data for ENMs is currently lacking. In an attempt to improve the FAIRness of already existing transcriptomics data for ENMs, we curated a collection of homogenized transcriptomics data from human, mouse and rat ENM exposures in vitro and in vivo including the physicochemical characteristics of the ENMs used in each study.

TinderMIX: Time-dose integrated modelling of toxicogenomics data.

BACKGROUND: Omics technologies have been widely applied in toxicology studies to investigate the effects of different substances on exposed biological systems. A classical toxicogenomic study consists in testing the effects of a compound at different dose levels and different time points. The main challenge consists in identifying the gene alteration patterns that are correlated to doses and time points. The majority of existing methods for toxicogenomics data analysis allow the study of the molecular alteration after the exposure (or treatment) at each time point individually. However, this kind of analysis cannot identify dynamic (time-dependent) events of dose responsiveness. RESULTS: We propose TinderMIX, an approach that simultaneously models the effects of time and dose on the transcriptome to investigate the course of molecular alterations exerted in response to the exposure. Starting from gene log fold-change, TinderMIX fits different integrated time and dose models to each gene, selects the optimal one, and computes its time and dose effect map; then a user-selected threshold is applied to identify the responsive area on each map and verify whether the gene shows a dynamic (time-dependent) and dose-dependent response; eventually, responsive genes are labelled according to the integrated time and dose point of departure. CONCLUSIONS: To showcase the TinderMIX method, we analysed 2 drugs from the Open TG-GATEs dataset, namely, cyclosporin A and thioacetamide. We first identified the dynamic dose-dependent mechanism of action of each drug and compared them. Our analysis highlights that different time- and dose-integrated point of departure recapitulates the toxicity potential of the compounds as well as their dynamic dose-dependent mechanism of action.

NanoSolveIT Project: Driving nanoinformatics research to develop innovative and integrated tools for in silico nanosafety assessment.

Nanotechnology has enabled the discovery of a multitude of novel materials exhibiting unique physicochemical (PChem) properties compared to their bulk analogues. These properties have led to a rapidly increasing range of commercial applications; this, however, may come at a cost, if an association to long-term health and environmental risks is discovered or even just perceived. Many nanomaterials (NMs) have not yet had their potential adverse biological effects fully assessed, due to costs and time constraints associated with the experimental assessment, frequently involving animals. Here, the available NM libraries are analyzed for their suitability for integration with novel nanoinformatics approaches and for the development of NM specific Integrated Approaches to Testing and Assessment (IATA) for human and environmental risk assessment, all within the NanoSolveIT cloud-platform. These established and well-characterized NM libraries (e.g. NanoMILE, NanoSolutions, NANoREG, NanoFASE, caLIBRAte, NanoTEST and the Nanomaterial Registry (>2000 NMs)) contain physicochemical characterization data as well as data for several relevant biological endpoints, assessed in part using harmonized Organisation for Economic Co-operation and Development (OECD) methods and test guidelines. Integration of such extensive NM information sources with the latest nanoinformatics methods will allow NanoSolveIT to model the relationships between NM structure (morphology), properties and their adverse effects and to predict the effects of other NMs for which less data is available. The project specifically addresses the needs of regulatory agencies and industry to effectively and rapidly evaluate the exposure, NM hazard and risk from nanomaterials and nano-enabled products, enabling implementation of computational 'safe-by-design' approaches to facilitate NM commercialization.

BMDx: a graphical Shiny application to perform Benchmark Dose analysis for transcriptomics data.

MOTIVATION: The analysis of dose-dependent effects on the gene expression is gaining attention in the field of toxicogenomics. Currently available computational methods are usually limited to specific omics platforms or biological annotations and are able to analyse only one experiment at a time. RESULTS: We developed the software BMDx with a graphical user interface for the Benchmark Dose (BMD) analysis of transcriptomics data. We implemented an approach based on the fitting of multiple models and the selection of the optimal model based on the Akaike Information Criterion. The BMDx tool takes as an input a gene expression matrix and a phenotype table, computes the BMD, its related values, and IC50/EC50 estimations. It reports interactive tables and plots that the user can investigate for further details of the fitting, dose effects and functional enrichment. BMDx allows a fast and convenient comparison of the BMD values of a transcriptomics experiment at different time points and an effortless way to interpret the results. Furthermore, BMDx allows to analyse and to compare multiple experiments at once. AVAILABILITY AND IMPLEMENTATION: BMDx is implemented as an R/Shiny software and is available at https://github.com/Greco-Lab/BMDx/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

FunMappOne: a tool to hierarchically organize and visually navigate functional gene annotations in multiple experiments.

BACKGROUND: Functional annotation of genes is an essential step in omics data analysis. Multiple databases and methods are currently available to summarize the functions of sets of genes into higher level representations, such as ontologies and molecular pathways. Annotating results from omics experiments into functional categories is essential not only to understand the underlying regulatory dynamics but also to compare multiple experimental conditions at a higher level of abstraction. Several tools are already available to the community to represent and compare functional profiles of omics experiments. However, when the number of experiments and/or enriched functional terms is high, it becomes difficult to interpret the results even when graphically represented. Therefore, there is currently a need for interactive and user-friendly tools to graphically navigate and further summarize annotations in order to facilitate results interpretation also when the dimensionality is high. RESULTS: We developed an approach that exploits the intrinsic hierarchical structure of several functional annotations to summarize the results obtained through enrichment analyses to higher levels of interpretation and to map gene related information at each summarized level. We built a user-friendly graphical interface that allows to visualize the functional annotations of one or multiple experiments at once. The tool is implemented as a R-Shiny application called FunMappOne and is available at https://github.com/grecolab/FunMappOne . CONCLUSION: FunMappOne is a R-shiny graphical tool that takes in input multiple lists of human or mouse genes, optionally along with their related modification magnitudes, computes the enriched annotations from Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, or Reactome databases, and reports interactive maps of functional terms and pathways organized in rational groups. FunMappOne allows a fast and convenient comparison of multiple experiments and an easy way to interpret results.