New approach methods to assess developmental and adult neurotoxicity for regulatory use: a PARC work package 5 project.

In the European regulatory context, rodent in vivo studies are the predominant source of neurotoxicity information. Although they form a cornerstone of neurotoxicological assessments, they are costly and the topic of ethical debate. While the public expects chemicals and products to be safe for the developing and mature nervous systems, considerable numbers of chemicals in commerce have not, or only to a limited extent, been assessed for their potential to cause neurotoxicity. As such, there is a societal push toward the replacement of animal models with in vitro or alternative methods. New approach methods (NAMs) can contribute to the regulatory knowledge base, increase chemical safety, and modernize chemical hazard and risk assessment. Provided they reach an acceptable level of regulatory relevance and reliability, NAMs may be considered as replacements for specific in vivo studies. The European Partnership for the Assessment of Risks from Chemicals (PARC) addresses challenges to the development and implementation of NAMs in chemical risk assessment. In collaboration with regulatory agencies, Project 5.2.1e (Neurotoxicity) aims to develop and evaluate NAMs for developmental neurotoxicity (DNT) and adult neurotoxicity (ANT) and to understand the applicability domain of specific NAMs for the detection of endocrine disruption and epigenetic perturbation. To speed up assay time and reduce costs, we identify early indicators of later-onset effects. Ultimately, we will assemble second-generation developmental neurotoxicity and first-generation adult neurotoxicity test batteries, both of which aim to provide regulatory hazard and risk assessors and industry stakeholders with robust, speedy, lower-cost, and informative next-generation hazard and risk assessment tools.

Evaluation of Physicochemical Properties of Ipsapirone Derivatives Based on Chromatographic and Chemometric Approaches.

Drug discovery is a challenging process, with many compounds failing to progress due to unmet pharmacokinetic criteria. Lipophilicity is an important physicochemical parameter that affects various pharmacokinetic processes, including absorption, metabolism, and excretion. This study evaluated the lipophilic properties of a library of ipsapirone derivatives that were previously synthesized to affect dopamine and serotonin receptors. Lipophilicity indices were determined using computational and chromatographic approaches. In addition, the affinity to human serum albumin (HSA) and phospholipids was assessed using biomimetic chromatography protocols. Quantitative Structure-Retention Relationship (QSRR) methodologies were used to determine the impact of theoretical descriptors on experimentally determined properties. A multiple linear regression (MLR) model was calculated to identify the most important features, and genetic algorithms (GAs) were used to assist in the selection of features. The resultant models showed commendable predictive accuracy, minimal error, and good concordance correlation coefficient values of 0.876, 0.149, and 0.930 for the validation group, respectively.

Environmental impact of PFAS: Filling data gaps using theoretical quantum chemistry and QSPR modeling.

Per- and polyfluorinated alkyl substances (PFAS), known for their widespread environmental presence and slow degradation, pose significant concerns. Of the approximately 10,000 known PFAS, only a few have undergone comprehensive testing, resulting in limited experimental data. In this study, we employed a combination of physics-based methods and data-driven models to address gaps in PFAS bioaccumulation potential. Using the COnductor-like Screening MOdel for Realistic Solvents (COSMO-RS) method, we predicted n-octanol/water partition coefficients (logKOW), crucial for PFAS bioaccumulation. Our developed Quantitative Structure-Property Relationship (QSPR) model exhibited high accuracy (R2 = 0.95, RMSEC = 0.75) and strong predictive ability (Q2LOO = 0.93, RMSECV = 0.83). Leveraging the extensive NORMAN, we predicted logKOW for over 4,000 compounds, identifying 244 outliers out of 4519. Further categorizing the database into eight Organisation for Economic Co-operation and Development (OECD) categories, we confirmed fluorine atoms role in enhanced bioaccumulation. Utilizing predicted logKOW, water solubility logSW, and vapor pressure logVP values, we calculated additional physicochemical properties that are responsible for the transport and dispersion of PFAS in the environment. Parameters such as Henry's Law (kH), air-water partition coefficient (KAW), octanol-air coefficient (KOA), and soil adsorption coefficient (KOC) exhibited favorable correlations with literature data (R2 > 0.66). Our study successfully filled data gaps, contributing to the understanding of ubiquitous PFAS in the environment and estimating missing physicochemical data for these compounds.

Toward Nano-Specific In Silico NAMs: How to Adjust Nano-QSAR to the Recent Advancements of Nanotoxicology?

The rapid development of engineered nanomaterials (ENMs) causes humans to become increasingly exposed to them. Therefore, a better understanding of the health impact of ENMs is highly demanded. Considering the 3Rs (Replacement, Reduction, and Refinement) principle, in vitro and computational methods are excellent alternatives for testing on animals. Among computational methods, nano-quantitative structure-activity relationship (nano-QSAR), which links the physicochemical and structural properties of EMNs with biological activities, is one of the leading method. The nature of toxicological experiments has evolved over the last decades; currently, one experiment can provide thousands of measurements of the organism's functioning at the molecular level. At the same time, the capacity of the in vitro systems to mimic the human organism is also improving significantly. Hence, the authors would like to discuss whether the nano-QSAR approach follows modern toxicological studies and takes full advantage of the opportunities offered by modern toxicological platforms. Challenges and possibilities for improving data integration are underlined narratively, including the need for a consensus built between the in vitro and the QSAR domains.

A bibliometric analysis of the recent achievements in pulmonary safety of nanoparticles.

Assessing research activity is an important step for planning future initiatives oriented toward filling the remaining gaps in a field. Therefore, the objective of the current study was to review recently published research on pulmonary toxicity caused by nanomaterials. However, here, instead of reviewing possible toxic effects and discussing their mode of action, the goal was to establish trends considering for example examined so far nanomaterials or used testing strategies. A total of 2316 related articles retrieved from the three most cited databases (PubMed Scopus, Web of Science), selected based on the title and abstract requirements, were used as the source of the review. Based on the bibliometric analysis, the nano-meter metal oxides, and carbon-based nanotubes were identified as the most frequently studied nanomaterials, while quantum dots, which might induce possible harmful effects, were not considered so far. The majority of testing of pulmonary safety is based on in vitro studies with observed growth of the contribution of novel testing strategies, such as 3D lung model, air-liquid interface system, or omic analysis.

The Effects of Two Kinds of Dietary Interventions on Serum Metabolic Profiles in Haemodialysis Patients.

The goal of this study was to evaluate the effects of two kinds of 24-week dietary interventions in haemodialysis patients, a traditional nutritional intervention without a meal before dialysis (HG1) and implementation of a nutritional intervention with a meal served just before dialysis (HG2), in terms of analysing the differences in the serum metabolic profiles and finding biomarkers of dietary efficacy. These studies were performed in two homogenous groups of patients (n = 35 in both groups). Among the metabolites with the highest statistical significance between HG1 and HG2 after the end of the study, 21 substances were putatively annotated, which had potential significance in both of the most relevant metabolic pathways and those related to diet. After the 24 weeks of the dietary intervention, the main differences between the metabolomic profiles in the HG2 vs. HG1 groups were related to the higher signal intensities from amino acid metabolites: indole-3-carboxaldehyde, 5-(hydroxymethyl-2-furoyl)glycine, homocitrulline, 4-(glutamylamino)butanoate, tryptophol, gamma-glutamylthreonine, and isovalerylglycine. These metabolites are intermediates in the metabolic pathways of the necessary amino acids (Trp, Tyr, Phe, Leu, Ile, Val, Liz, and amino acids of the urea cycle) and are also diet-related intermediates (4-guanidinobutanoic acid, indole-3-carboxyaldehyde, homocitrulline, and isovalerylglycine).

Core, Coating, or Corona? The Importance of Considering Protein Coronas in nano-QSPR Modeling of Zeta Potential.

To control stability in a biological medium, several factors affecting the zeta potential (ζ) of nanoparticles (NPs) must be considered, including complex interactions between the nanostructure and the composition of the protein corona (PC). Effective in silico methods (based on machine learning and quantitative structure-property relationship (QSPR) models) could help predict and characterize the relationship between the physicochemical properties of NP and the formation of PC and biological outcomes in the medium at an early stage of the experiment. However, the models currently developed are limited to simple descriptors that do not represent the complex interactions between the core, the coating, and their PC fingerprints. To be useful, the models developed should be described as a function of both the structural properties determined by the core and coating of the NPs and the biological medium determined by the formation of the protein corona. We have developed a set of complex descriptors that describe the quantitative relationship between the value of the zeta potential (ζ), core, the coating of NPs, and their PC fingerprints (the so-called nano-QSPR model). The nano-QSPR model was developed based on a genetic algorithm using a partial least-squares regression method (GA-PLS), which is characterized by high external predictive power (Q2EXT = 0.89). The GA-PLS model was developed using descriptors that describe (i) the core structure (determined by 7 different types of polymer-based NMs in the range of 20 different sizes), (ii) the coating structure with 7 different functional groups, and (iii) 80 different types of protein compositions adsorbed on the surface of the NPs. The presented study answers the question of how complex interactions between the corona and NP determine the zeta potential (ζ) of NP in a given medium. Moreover, our current study is a proof-of-concept that the zeta potential of NPs modeled on the original structure depends not only on the NPs themselves but also on the structure and properties determined by the NP core and coating, as well as the biological medium determined by the formation of the protein corona. On the basis of these results, our studies will be useful in determining the stability and mechanism of cell uptake, toxicity, and ability to predict the zeta potential of compounds not yet tested.

Representing and describing nanomaterials in predictive nanoinformatics.

Engineered nanomaterials (ENMs) enable new and enhanced products and devices in which matter can be controlled at a near-atomic scale (in the range of 1 to 100 nm). However, the unique nanoscale properties that make ENMs attractive may result in as yet poorly known risks to human health and the environment. Thus, new ENMs should be designed in line with the idea of safe-and-sustainable-by-design (SSbD). The biological activity of ENMs is closely related to their physicochemical characteristics, changes in these characteristics may therefore cause changes in the ENMs activity. In this sense, a set of physicochemical characteristics (for example, chemical composition, crystal structure, size, shape, surface structure) creates a unique 'representation' of a given ENM. The usability of these characteristics or nanomaterial descriptors (nanodescriptors) in nanoinformatics methods such as quantitative structure-activity/property relationship (QSAR/QSPR) models, provides exciting opportunities to optimize ENMs at the design stage by improving their functionality and minimizing unforeseen health/environmental hazards. A computational screening of possible versions of novel ENMs would return optimal nanostructures and manage ('design out') hazardous features at the earliest possible manufacturing step. Safe adoption of ENMs on a vast scale will depend on the successful integration of the entire bulk of nanodescriptors extracted experimentally with data from theoretical and computational models. This Review discusses directions for developing appropriate nanomaterial representations and related nanodescriptors to enhance the reliability of computational modelling utilized in designing safer and more sustainable ENMs.

Predicting electrophoretic mobility of TiO2, ZnO, and CeO2 nanoparticles in natural waters: The importance of environment descriptors in nanoinformatics models.

Natural and engineered nanoparticles (NPs) entering the environment are influenced by many physicochemical processes and show various behavior in different systems (e.g., natural waters showing different characteristics). Determining the physicochemical characteristics and predicting the behavior of nanoparticles ending up in the natural aquatic environment are key aspects of their risk assessment. Here, we show that the quantitative structure-property relationship modeling method used in nanoinformatics (nano-QSPR) can be successfully applied to predict environmental fate-relevant properties (electrophoretic mobility) of TiO2, ZnO, and CeO2 nanoparticles. However, in contrast to the previous works, we postulate to use, in parallel: (i) the nanoparticles' structure descriptors (S-descriptors) and (ii) the environment descriptors (E-descriptors) as the input variables. Thus, the method should be abbreviated more precisely as nano-QSEPR ("E" stands for the "environment"). As a proof-of-the-concept, we have developed a group of models (including MLR, GA-PLS, PCR, and Meta-Consensus models) with high predictive capabilities (QEXT2 = 0.931 for the GA-PLS model), where the S-descriptors are represented by the core-shell model descriptor and the E-descriptors - by different ambient water features (including ions concentration and the ionic strength). The newly proposed nano-QSEPR modeling scheme can be efficiently used to design safe and sustainable nanomaterials.

Improving lithium dose prediction using population pharmacokinetics and pharmacogenomics: a cohort genome-wide association study in Sweden.

BACKGROUND: Lithium is the most effective treatment for bipolar disorder, resulting in strong suicide prevention effects. The therapeutic range of lithium, however, is narrow and treatment initiation requires individual titration to address inter-individual variability. We aimed to improve lithium dose prediction using clinical and genomic data. METHODS: We performed a population pharmacokinetic study followed by a genome-wide association study (GWAS), including two clinical Swedish cohorts. Participants in cohort 1 were from specialised outpatient clinics at Huddinge Hospital, in Stockholm, Sweden, and participants in cohort 2 were identified using the Swedish National Quality Registry for Bipolar disorder (BipoläR). Patients who received a lithium dose corresponding to at least one tablet of lithium sulphate (6 mmol) per day and had clinically relevant plasma concentrations of lithium were included in the study. Data on age, sex, bodyweight, height, creatinine concentration, estimated glomerular filtration rate (eGFR), lithium preparation, number of tablets of lithium per day, serum lithium concentration, and medications affecting kidney function (C09 antihypertensives, C03 [except C03D] sodium-retaining diuretics, and non-steroidal anti-inflammatory drugs) were obtained retrospectively for several timepoints when possible from electronic health records, BipoläR, and the Swedish prescription registry. The median time between timepoints was 1·07 years for cohort 1 and 1·09 years for cohort 2. The primary outcome of interest was the natural logarithm of total body clearance for lithium (CLLi) associated with the clinical variables. The residual effects after accounting for age and sex, representing the individual-level effects (CLLi,age/sex), were used as the dependent variable in a GWAS. FINDINGS: 2357 patients who were administered lithium (1423 women [60·4%] and 934 men [39·6%]; mean age 53·6 years [range 17-89], mainly of European descent) were included and 5627 data points were obtained. Age (variance explained [R2]: R2cohort1=0·41 and R2cohort2=0·31; both p<0·0001), sex (R2cohort1=0·0063 [p=0·045] and R2cohort2=0·026 [p<0·0001]), eGFR (R2cohort1=0·38 and R2cohort2=0·20; both p<0·0001), comedication with diuretics (R2cohort1=0·0058 [p=0·014] and R2cohort2=0·0026 [p<0·0001]), and agents acting on the renin-aldosterone-angiotensin system (R2cohort1=0·028 and R2cohort2=0·015; both p<0·0001) were clinical predictors of CLLi. Notably, an association between CLLi and serum lithium was observed, with a lower CLLi being associated with higher serum lithium (R2cohort1=0·13 and R2cohort2=0·15; both p<0·0001). In a GWAS of CLLi,age/sex, one locus was associated with a change in CLLi (rs583503; ß=-0·053 [95% CI -0·071 to -0·034]; p<0·00000005). We also found enrichment of the associations with genes expressed in the medulla (p=0·0014, corrected FDR=0·04) and cortex of the kidney (p=0·0015, corrected FDR=0·04), as well as associations with polygenic risk scores for eGFR (p value threshold: 0·05, p=0·01), body-mass index (p value threshold: 0·05, p=0·00025), and blood urea nitrogen (p value threshold: 0·001, p=0·00043). The model based on six clinical predictors explained 61·4% of the variance in CLLi in cohort 1 and 49·8% in cohort 2. Adding genetic markers did not lead to major improvement of the models: within the subsample of genotyped individuals, the variance explained only increased from 59·32% to 59·36% in cohort 1 and from 49·21% to 50·03% in cohort 2 when including rs583503 and the four first principal components. INTERPRETATION: Our model predictors could be used clinically to better guide lithium dosage, shortening the time to reach therapeutic concentrations, thus improving care. Identification of the first genomic locus and PRS to be associated with CLLi introduces the opportunity of individualised medicine in lithium treatment. FUNDING: Stanley Medical Research Institute, Swedish Research Council, Swedish Foundation for Strategic Research, Swedish Brain Foundation, Swedish Research Council, Söderström-Königska Foundation, Bror Gadelius Minnesfond, Swedish Mental Health Fund, Karolinska Institutet and Hospital.

HBM4EU Chromates Study: Urinary Metabolomics Study of Workers Exposed to Hexavalent Chromium.

Exposure to hexavalent chromium Cr(VI) may occur in several occupational activities, placing workers in many industries at risk for potential related health outcomes. Untargeted metabolomics was applied to investigate changes in metabolic pathways in response to Cr(VI) exposure. We obtained our data from a study population of 220 male workers with exposure to Cr(VI) and 102 male controls from Belgium, Finland, Poland, Portugal and the Netherlands within the HBM4EU Chromates Study. Urinary metabolite profiles were determined using liquid chromatography mass spectrometry, and differences between post-shift exposed workers and controls were analyzed using principal component analysis. Based on the first two principal components, we observed clustering by industrial chromate application, such as welding, chrome plating, and surface treatment, distinct from controls and not explained by smoking status or alcohol use. The changes in the abundancy of excreted metabolites observed in workers reflect fatty acid and monoamine neurotransmitter metabolism, oxidative modifications of amino acid residues, the excessive formation of abnormal amino acid metabolites and changes in steroid and thyrotropin-releasing hormones. The observed responses could also have resulted from work-related factors other than Cr(VI). Further targeted metabolomics studies are needed to better understand the observed modifications and further explore the suitability of urinary metabolites as early indicators of adverse effects associated with exposure to Cr(VI).

In vitro to in vivo extrapolation to support the development of the next generation risk assessment (NGRA) strategy for nanomaterials.

There is growing interest in developing novel strategies to support assessment of human health risks due to chemicals. Regulatory and decision-making agencies have recommended that non-animal-based alternatives should be applied whenever possible instead of experimentation on living animals. These alternative methods are beneficial because they are ethical, inexpensive, and rapid. Herein, we review recent activities aimed at developing in vitro to in vivo extrapolation (IVIVE) models as a part of the Next Generation Risk Assessment (NGRA) of nanomaterials. In this context, we show the adverse outcome pathway (AOP)-based methodology for the identification of mechanistically relevant events serving as biomarkers for the targeted selection of in vitro assays. Considered events need to be biologically plausible, regulatory relevant, and crucial for the examination of occurrence of adverse outcomes. The promising advantages of using high-throughout-based omics data are highlighted. Furthermore, the application of 3D in vitro models and nano genome atlases to study nanoparticle toxicity is briefly summarized. Additionally, the challenges related to the extrapolation of in vitro doses into in vivo-relevant responses are presented. We also discuss the limitations of models applied thus far to study the fate of chemicals in the human body, which exist due to the lack of available knowledge regarding transformations of nanomaterials occurring in biological systems.

AOP173 key event associated pathway predictor - online application for the prediction of benchmark dose lower bound (BMDLs) of a transcriptomic pathway involved in MWCNTs-induced lung fibrosis.

Nano-QSAR model allows for prediction of the toxicity of materials that have not been experimentally tested before by linking the nano-related structural properties with the biological responses induced by nanomaterials. Prediction of adverse effects caused by substances without having to perform time- and cost-consuming experiments makes QSAR models promising tools for supporting risk assessment. However, very often, newly developed nano-QSAR models are not used in practice due to the complexity of their algorithms, the necessity to have experience in chemoinformatics, and their poor accessibility. In this perspective, the aim of this paper is to encourage developers of the QSAR models to take the effort to prepare user-friendly applications based on predictive models. This would make the developed models accessible to a wider community, and, in effect, promote their further application by regulators and decision-makers. Here, we describe a web-based application that enables to predict the transcriptomic pathway-level response perturbated in the lungs of mice exposed to multiwalled carbon nanotubes. The developed application is freely available at http://aop173-event1.nanoqsar-aop.com/apps/aop_app. It requires only two types of input information related to analyzed nanotubes (their length and diameter) to assess the doses that initiate the inflammation process that may lead to lung fibrosis.

Fusobacterium nucleatum - Friend or foe?

Fusobacterium nucleatum (F. nucleatum) is one of the most abundant Gram-negative anaerobic bacteria, part of the gut, and oral commensal flora, generally found in human dental plaque. Its presence could be associated with various human diseases, including, e.g., periodontal, angina, lung and gynecological abscesses. This bacteria can enter the blood circulation as a result of periodontal infection. It was proven that F. nucleatum migrates from its primary site of colonization in the oral cavity to other parts of the body. It could cause numerous diseases, including cancers. On the other hand, it was shown that Fusobacterium produces significant amounts of butyric acid, which is a great source of energy for colonocytes (anti-inflammatory cells). Therefore, it is very interesting to get to know the two faces of F. nucleatum.

Comparative studies on vitamin B1 deficiency inwholeblood of chronically haemodialysed patients: chromatographic, fluorimetricandPCAstudy.

The levels of thiamine diphosphate (ThDP), the most active biologically form of vitamin B1, were assessed in whole blood oflong-term haemodialysed patients (n = 50), by applying chromatographic methods based on RP-HPLC technique with isocratic elution and fluorescence detection. The target analyte, thiochrome diphosphate (ThODP), was obtained by pre-column derivatization of vitamin B1 contained in blood samples, applying deproteination with trichloroacetic acid, following by oxidation with alkaline solution of potassium ferricyanide(III) and stabilization with DTT before assays. A simple and sensitive assay was developed, and the results were referenced to the commercially available test. Steady-state and time-resolved studies on emissive properties of ThODP enabled optimization of the proposed assay. The F-Snedecor test shown no statistically significant differences between both approaches. Assessed parameters of the proposed assay, such as linearity, precision, sensitivity, and recovery, were satisfactory if compared to the reference one. The LOQ value for ThDP in whole blood of studied group of patients was of 0.5 ng/mL and the recovery of88%. The results disclosed high individual variabilities in the interdialytic deficiencies of ThDP among the patients - ranged from afew percent to values close to 100%. A comprehensive clinical data, characterizing patients under study, were processed together, and analysed by employing achemometric discriminative tool, the Principal Components Analysis,to find interdependences among clinical data characterizing patients. The three Principal Components were disclosed, that in sum explained almost 50% of the observed variability of the clinical data set. Among the clinical parameters involved in PCs were dialyzer membrane and type, duration as well as levels of creatinine, haemoglobin, and red blood cells in patients' whole blood.

A strategy towards the generation of testable adverse outcome pathways for nanomaterials.

Manufactured nanomaterials (NMs) are increasingly used in a wide range of industrial applications leading to a constant increase in the market size of nano-enabled products. The increased production and use of NMs are raising concerns among different stakeholder groups with regard to their effects on human and environmental health. Currently, nanosafety hazard assessment is still widely performed using in vivo (animal) models, however the development of robust and reg­ulatory relevant strategies is required to prioritize and/or reduce animal testing. An adverse outcome pathway (AOP) is a structured representation of biological events that start from a molecular initiating event (MIE) leading to an adverse outcome (AO) through a series of key events (KEs). The AOP framework offers great advancement to risk assessment and regulatory safety assessments. While AOPs for chemicals have been more frequently reported, the AOP collection for NMs is limited. By using existing AOPs, we aimed to generate simple and testable strategies to predict if a given NM has the potential to induce a MIE leading to an AO through a series of KEs. Firstly, we identified potential MIEs or initial KEs reported for NMs in the literature. Then, we searched the identified MIE or initial KEs as keywords in the AOP-Wiki to find associated AOPs. Finally, using two case studies, we demonstrate how in vitro strategies can be used to test the identified MIE/KEs.

Transcriptomics-Based and AOP-Informed Structure-Activity Relationships to Predict Pulmonary Pathology Induced by Multiwalled Carbon Nanotubes.

This study presents a novel strategy that employs quantitative structure-activity relationship models for nanomaterials (Nano-QSAR) for predicting transcriptomic pathway level response using lung tissue inflammation, an essential key event (KEs) in the existing adverse outcome pathway (AOP) for lung fibrosis, as a model response. Transcriptomic profiles of mouse lungs exposed to ten different multiwalled carbon nanotubes (MWCNTs) are analyzed using statistical and bioinformatics tools. Three pathways "agranulocyte adhesion and diapedesis," "granulocyte adhesion and diapedesis," and "acute phase signaling," that (1) are commonly perturbed across the MWCNTs panel, (2) show dose response (Benchmark dose, BMDs), and (3) are anchored to the KEs identified in the lung fibrosis AOP, are considered in modelling. The three pathways are associated with tissue inflammation. The results show that the aspect ratio (κ) of MWCNTs is directly correlated with the pathway BMDs. The study establishes a methodology for QSAR construction based on canonical pathways and proposes a MWCNTs grouping strategy based on the κ-values of the specific pathway associated genes. Finally, the study shows how the AOP framework can help guide QSAR modelling efforts; conversely, the outcome of the QSAR modelling can aid in refining certain aspects of the AOP in question (here, lung fibrosis).

Transcriptomics in Toxicogenomics, Part II: Preprocessing and Differential Expression Analysis for High Quality Data.

Preprocessing of transcriptomics data plays a pivotal role in the development of toxicogenomics-driven tools for chemical toxicity assessment. The generation and exploitation of large volumes of molecular profiles, following an appropriate experimental design, allows the employment of toxicogenomics (TGx) approaches for a thorough characterisation of the mechanism of action (MOA) of different compounds. To date, a plethora of data preprocessing methodologies have been suggested. However, in most cases, building the optimal analytical workflow is not straightforward. A careful selection of the right tools must be carried out, since it will affect the downstream analyses and modelling approaches. Transcriptomics data preprocessing spans across multiple steps such as quality check, filtering, normalization, batch effect detection and correction. Currently, there is a lack of standard guidelines for data preprocessing in the TGx field. Defining the optimal tools and procedures to be employed in the transcriptomics data preprocessing will lead to the generation of homogeneous and unbiased data, allowing the development of more reliable, robust and accurate predictive models. In this review, we outline methods for the preprocessing of three main transcriptomic technologies including microarray, bulk RNA-Sequencing (RNA-Seq), and single cell RNA-Sequencing (scRNA-Seq). Moreover, we discuss the most common methods for the identification of differentially expressed genes and to perform a functional enrichment analysis. This review is the second part of a three-article series on Transcriptomics in Toxicogenomics.

Transcriptomics in Toxicogenomics, Part III: Data Modelling for Risk Assessment.

Transcriptomics data are relevant to address a number of challenges in Toxicogenomics (TGx). After careful planning of exposure conditions and data preprocessing, the TGx data can be used in predictive toxicology, where more advanced modelling techniques are applied. The large volume of molecular profiles produced by omics-based technologies allows the development and application of artificial intelligence (AI) methods in TGx. Indeed, the publicly available omics datasets are constantly increasing together with a plethora of different methods that are made available to facilitate their analysis, interpretation and the generation of accurate and stable predictive models. In this review, we present the state-of-the-art of data modelling applied to transcriptomics data in TGx. We show how the benchmark dose (BMD) analysis can be applied to TGx data. We review read across and adverse outcome pathways (AOP) modelling methodologies. We discuss how network-based approaches can be successfully employed to clarify the mechanism of action (MOA) or specific biomarkers of exposure. We also describe the main AI methodologies applied to TGx data to create predictive classification and regression models and we address current challenges. Finally, we present a short description of deep learning (DL) and data integration methodologies applied in these contexts. Modelling of TGx data represents a valuable tool for more accurate chemical safety assessment. This review is the third part of a three-article series on Transcriptomics in Toxicogenomics.

Transcriptomics in Toxicogenomics, Part I: Experimental Design, Technologies, Publicly Available Data, and Regulatory Aspects.

The starting point of successful hazard assessment is the generation of unbiased and trustworthy data. Conventional toxicity testing deals with extensive observations of phenotypic endpoints in vivo and complementing in vitro models. The increasing development of novel materials and chemical compounds dictates the need for a better understanding of the molecular changes occurring in exposed biological systems. Transcriptomics enables the exploration of organisms' responses to environmental, chemical, and physical agents by observing the molecular alterations in more detail. Toxicogenomics integrates classical toxicology with omics assays, thus allowing the characterization of the mechanism of action (MOA) of chemical compounds, novel small molecules, and engineered nanomaterials (ENMs). Lack of standardization in data generation and analysis currently hampers the full exploitation of toxicogenomics-based evidence in risk assessment. To fill this gap, TGx methods need to take into account appropriate experimental design and possible pitfalls in the transcriptomic analyses as well as data generation and sharing that adhere to the FAIR (Findable, Accessible, Interoperable, and Reusable) principles. In this review, we summarize the recent advancements in the design and analysis of DNA microarray, RNA sequencing (RNA-Seq), and single-cell RNA-Seq (scRNA-Seq) data. We provide guidelines on exposure time, dose and complex endpoint selection, sample quality considerations and sample randomization. Furthermore, we summarize publicly available data resources and highlight applications of TGx data to understand and predict chemical toxicity potential. Additionally, we discuss the efforts to implement TGx into regulatory decision making to promote alternative methods for risk assessment and to support the 3R (reduction, refinement, and replacement) concept. This review is the first part of a three-article series on Transcriptomics in Toxicogenomics. These initial considerations on Experimental Design, Technologies, Publicly Available Data, Regulatory Aspects, are the starting point for further rigorous and reliable data preprocessing and modeling, described in the second and third part of the review series.