Explore the vast range of titles available.

Steatosis, characterized by excessive fat accumulation in the liver, is a significant precursor to chronic liver disease and hepatocarcinoma. This condition is influenced by multiple contributing factors such as obesity, alcohol consumption, and exposure to chemicals or drugs. Systems biology approaches including transcriptomics and metabolomics can aid in grouping chemicals according to their mode of action. In this study, we analyze transcriptomic and metabolomic data from primary human and transformed hepatocytes, respectively, to differentiate between steatotic and non-steatotic chemicals. Rather than assessing each steatotic compound individually, we pooled several steatotic chemicals in order to minimize compound-specific noise and better identify features associated with the underlying process of steatosis. Differential gene expression analysis revealed established mechanisms involved in steatosis, consistent with the recently updated adverse outcome pathway. Likewise, metabolomic data enabled clear discrimination between steatotic and non-steatotic chemicals. These findings highlight the potential of omics technologies to support chemical grouping based on insights into the molecular mechanisms that drive steatosis development.

Adaptive stress response pathways play a key role in the switch between adaptation and adversity, and are important in drug-induced liver injury. Previously, we have established an HepG2 fluorescent protein reporter platform to monitor adaptive stress response activation following drug treatment. HepG2 cells are often used in high-throughput primary toxicity screening, but metabolizing capacity in these cells is low and repeated dose toxicity testing inherently difficult. Here, we applied our bacterial artificial chromosome-based GFP reporter cell lines representing Nrf2 activation (Srxn1-GFP and NQO1-GFP), unfolded protein response (BiP-GFP and Chop-GFP), and DNA damage response (p21-GFP and Btg2-GFP) as long-term differentiated 3D liver-like spheroid cultures. All HepG2 GFP reporter lines differentiated into 3D spheroids similar to wild-type HepG2 cells. We systematically optimized the automated imaging and quantification of GFP reporter activity in individual spheroids using high-throughput confocal microscopy with a reference set of DILI compounds that activate these three stress response pathways at the transcriptional level in primary human hepatocytes. A panel of 33 compounds with established DILI liability was further tested in these six 3D GFP reporters in single 48 h treatment or 6 day daily repeated treatment. Strongest stress response activation was observed after 6-day repeated treatment, with the BiP and Srxn1-GFP reporters being most responsive and identified particular severe-DILI-onset compounds. Compounds that showed no GFP reporter activation in two-dimensional (2D) monolayer demonstrated GFP reporter stress response activation in 3D spheroids. Our data indicate that the application of BAC-GFP HepG2 cellular stress reporters in differentiated 3D spheroids is a promising strategy for mechanism-based identification of compounds with liability for DILI.

Thirteen pharmaceutical companies have shared and integrated preclinical and clinical data for creating computational resources that enhance translational drug safety assessment.Thirteen pharmaceutical companies have shared and integrated preclinical and clinical data for creating computational resources that enhance translational drug safety assessment.

Mechanism-based risk assessment is urged to advance and fully permeate into current safety assessment practices, possibly at early phases of drug safety testing. Toxicogenomics is a promising source of mechanisms-revealing data, but interpretative analysis tools specific for the testing systems (e.g. hepatocytes) are lacking. In this study, we present the TXG-MAPr webtool (available at https://txg-mapr.eu/WGCNA_PHH/TGGATEs_PHH/), an R-Shiny-based implementation of weighted gene co-expression network analysis (WGCNA) obtained from the Primary Human Hepatocytes (PHH) TG-GATEs dataset. The 398 gene co-expression networks (modules) were annotated with functional information (pathway enrichment, transcription factor) to reveal their mechanistic interpretation. Several well-known stress response pathways were captured in the modules, were perturbed by specific stressors and showed preservation in rat systems (rat primary hepatocytes and rat in vivo liver), with the exception of DNA damage and oxidative stress responses. A subset of 87 well-annotated and preserved modules was used to evaluate mechanisms of toxicity of endoplasmic reticulum (ER) stress and oxidative stress inducers, including cyclosporine A, tunicamycin and acetaminophen. In addition, module responses can be calculated from external datasets obtained with different hepatocyte cells and platforms, including targeted RNA-seq data, therefore, imputing biological responses from a limited gene set. As another application, donors’ sensitivity towards tunicamycin was investigated with the TXG-MAPr, identifying higher basal level of intrinsic immune response in donors with pre-existing liver pathology. In conclusion, we demonstrated that gene co-expression analysis coupled to an interactive visualization environment, the TXG-MAPr, is a promising approach to achieve mechanistic relevant, cross-species and cross-platform evaluation of toxicogenomic data.

Mitochondrial perturbation is a key event in chemical-induced organ toxicities that is incompletely understood. Here, we studied how electron transport chain (ETC) complex I, II, or III (CI, CII and CIII) inhibitors affect mitochondrial functionality, stress response activation, and cell viability using a combination of high-content imaging and TempO-Seq in HepG2 hepatocyte cells. CI and CIII inhibitors perturbed mitochondrial membrane potential (MMP) and mitochondrial and cellular ATP levels in a concentration- and time-dependent fashion and, under conditions preventing a switch to glycolysis attenuated cell viability, whereas CII inhibitors had no effect. TempO-Seq analysis of changes in mRNA expression pointed to a shared cellular response to CI and CIII inhibition. First, to define specific ETC inhibition responses, a gene set responsive toward ETC inhibition (and not to genotoxic, oxidative, or endoplasmic reticulum stress) was identified using targeted TempO-Seq in HepG2. Silencing of one of these genes, NOS3, exacerbated the impact of CI and CIII inhibitors on cell viability, indicating its functional implication in cellular responses to mitochondrial stress. Then by monitoring dynamic responses to ETC inhibition using a HepG2 GFP reporter panel for different classes of stress response pathways and applying pathway and gene network analysis to TempO-Seq data, we looked for downstream cellular events of ETC inhibition and identified the amino acid response (AAR) as being triggered in HepG2 by ETC inhibition. Through in silico approaches we provide evidence indicating that a similar AAR is associated with exposure to mitochondrial toxicants in primary human hepatocytes. Altogether, we (i) unravel quantitative, time- and concentration-resolved cellular responses to mitochondrial perturbation, (ii) identify a gene set associated with adaptation to exposure to active ETC inhibitors, and (iii) show that ER stress and an AAR accompany ETC inhibition in HepG2 and primary hepatocytes.

Background Drug-induced liver injury (DILI) is an adverse reaction caused by the intake of drugs of common use that produces liver damage. The impact of DILI is estimated to affect around 20 in 100,000 inhabitants worldwide each year. Despite being one of the main causes of liver failure, the pathophysiology and mechanisms of DILI are poorly understood. In the present study, we developed an ensemble learning approach based on different features (CMap gene expression, chemical structures, drug targets) to predict drugs that might cause DILI and gain a better understanding of the mechanisms linked to the adverse reaction. Results We searched for gene signatures in CMap gene expression data by using two approaches: phenotype-gene associations data from DisGeNET, and a non-parametric test comparing gene expression of DILI-Concern and No-DILI-Concern drugs (as per DILIrank definitions). The average accuracy of the classifiers in both approaches was 69%. We used chemical structures as features, obtaining an accuracy of 65%. The combination of both types of features produced an accuracy around 63%, but improved the independent hold-out test up to 67%. The use of drug-target associations as feature obtained the best accuracy (70%) in the independent hold-out test. Conclusions When using CMap gene expression data, searching for a specific gene signature among the landmark genes improves the quality of the classifiers, but it is still limited by the intrinsic noise of the dataset. When using chemical structures as a feature, the structural diversity of the known DILI-causing drugs hampers the prediction, which is a similar problem as for the use of gene expression information. The combination of both features did not improve the quality of the classifiers but increased the robustness as shown on independent hold-out tests. The use of drug-target associations as feature improved the prediction, specially the specificity, and the results were comparable to previous research studies.

The transcription factor NRF2, governed by its repressor KEAP1, protects cells against oxidative stress. There is interest in modelling the NRF2 response to improve the prediction of clinical toxicities such as drug-induced liver injury (DILI). However, very little is known about the makeup of the NRF2 transcriptional network and its response to chemical perturbation in primary human hepatocytes (PHH), which are often used as a translational model for investigating DILI. Here, microarray analysis identified 108 transcripts (including several putative novel NRF2-regulated genes) that were both downregulated by siRNA targeting NRF2 and upregulated by siRNA targeting KEAP1 in PHH. Applying weighted gene co-expression network analysis (WGCNA) to transcriptomic data from the Open TG-GATES toxicogenomics repository (representing PHH exposed to 158 compounds) revealed four co-expressed gene sets or ‘modules’ enriched for these and other NRF2-associated genes. By classifying the 158 TG-GATES compounds based on published evidence, and employing the four modules as network perturbation metrics, we found that the activation of NRF2 is a very good indicator of the intrinsic biochemical reactivity of a compound (i.e. its propensity to cause direct chemical stress), with relatively high sensitivity, specificity, accuracy and positive/negative predictive values. We also found that NRF2 activation has lower sensitivity for the prediction of clinical DILI risk, although relatively high specificity and positive predictive values indicate that false positive detection rates are likely to be low in this setting. Underpinned by our comprehensive analysis, activation of the NRF2 network is one of several mechanism-based components that can be incorporated into holistic systems toxicology models to improve mechanistic understanding and preclinical prediction of DILI in man.

eTRANSAFE is a research project funded within the Innovative Medicines Initiative (IMI), which aims at developing integrated databases and computational tools (the eTRANSAFE ToxHub) that support the translational safety assessment of new drugs by using legacy data provided by the pharmaceutical companies that participate in the project. The project objectives include the development of databases containing preclinical and clinical data, computational systems for translational analysis including tools for data query, analysis and visualization, as well as computational models to explain and predict drug safety events.

Toxicogenomic data represent a valuable source of biological information at molecular and cellular level to understand unanticipated organ toxicities. Weighted gene co-expression networks analysis can reduce the complexity of gene-level transcriptomic data to a set of biological response-networks useful for providing insights into mechanisms of drug-induced adverse outcomes. In this study, we have built co-regulated gene networks (modules) from the TG-GATEs and DrugMatrix rat kidney datasets consisting of time- and dose-response data for 180 compounds, including nephrotoxicants. Data from the 347 modules were incorporated into the rat kidney TXG-MAPr web tool, a user-friendly interface that enables visualization and analysis of module perturbations, quantified by a module eigengene score (EGS) for each treatment condition. Several modules annotated for cellular stress, renal injury and inflammation were statistically associated with concurrent renal pathologies, including modules that contain both well-known and novel renal biomarker genes. In addition, many rat kidney modules contain well annotated, robust gene networks that are preserved across transcriptome datasets, suggesting that these biological networks translate to other (drug-induced) kidney injury cases. Moreover, preservation analysis of human kidney transcriptomic data provided a quantitative metric to assess the likelihood that rat kidney modules, and the associated biological interpretation, translate from non-clinical species to human. In conclusion, the rat kidney TXG-MAPr enables uploading and analysis of kidney gene expression data in the context of rat kidney co-expression networks, which could identify possible safety liabilities and/or mechanisms that can lead to adversity for chemical or drug candidates.Competing Interest StatementJSR reports funding from GSK, Pfizer and Sanofi and fees/honoraria from Travere Therapeutics, Stadapharm, Astex, Owkin, Pfizer and Grunenthal. PT is a Sanofi employee and may hold shares and/or stock options in the company. All the other authors have declared no competing interests.Footnotes* Abstract and introduction slightly updated; Methods section updated by more clearly describing the module association with pathology, including a resource table and including a link to the R-script to run the WGCNA and module preservation; Discussion section updated by including comparison of co-expression networks with literature and to make some conclusions more clear; Figure 5 revised; Included a new Figure S10; Author affiliations updated; Updated some supplementary figures and tables.* https://txg-mapr.eu/login/