Explore the vast range of titles available.

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.

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.

Nephrotoxicity caused by drug or chemical exposure involves complex mechanisms as well as a temporal integration of injury and repair responses in different nephron segments. Distinct cellular transcriptional programs regulate the time-dependent tissue injury and regeneration responses. Whole kidney transcriptome analysis cannot dissect the complex spatio-temporal injury and regeneration responses in the different nephron segments. Here, we used laser capture microdissection of formalin-fixed paraffin embedded sections followed by whole genome targeted RNA-sequencing-TempO-Seq and co-expression gene-network (module) analysis to determine the spatial–temporal responses in rat kidney glomeruli (GM), cortical proximal tubules (CPT) and outer-medulla proximal tubules (OMPT) comparison with whole kidney, after a single dose of the nephrotoxicant cisplatin. We demonstrate that cisplatin induced early onset of DNA damage in both CPT and OMPT, but not GM. Sustained DNA damage response was strongest in OMPT coinciding with OMPT specific inflammatory signaling, actin cytoskeletal remodeling and increased glycolytic metabolism with suppression of mitochondrial activity. Later responses reflected regeneration-related cell cycle pathway activation and ribosomal biogenesis in the injured OMPT regions. Activation of modules containing kidney injury biomarkers was strongest in OMPT, with OMPT Clu expression highly correlating with urinary clusterin biomarker measurements compared the correlation of Kim1. Our findings also showed that whole kidney responses were less sensitive than OMPT. In conclusion, our LCM-TempO-Seq method reveals a detailed spatial mechanistic understanding of renal injury/regeneration after nephrotoxicant exposure and identifies the most representative mechanism-based nephron segment specific renal injury biomarkers. Graphical Abstract Highlights • Different nephron segments exhibit distinct transcriptomic perturbation with different degrees of sensitivity. • Sustained activation of DNA damage responses upon cisplatin exposure is linked to progressive outcomes of injured nephron regions. • Mechanistic kidney injury biomarkers such as urinary clusterin outperform conventional biomarkers in reflecting the condition of the damaged nephron segments.

Renal tubular epithelial cells are exposed to mechanical forces due to fluid flow shear stress within the lumen of the nephron. These cells respond by activation of mechano-sensors located at the plasma membrane or the primary cilium, having crucial roles in maintenance of cellular homeostasis and signaling. In this paper, we applied fluid shear stress to study TGF-β signaling in renal epithelial cells with and without expression of the Pkd1 -gene, encoding a mechano-sensor mutated in polycystic kidney disease. TGF-β signaling modulates cell proliferation, differentiation, apoptosis, and fibrotic deposition, cellular programs that are altered in renal cystic epithelia. SMAD2/3-mediated signaling was activated by fluid flow, both in wild-type and Pkd1 −/− cells. This was characterized by phosphorylation and nuclear accumulation of p-SMAD2/3, as well as altered expression of downstream target genes and epithelial-to-mesenchymal transition markers. This response was still present after cilia ablation. An inhibitor of upstream type-I-receptors, ALK4/ALK5/ALK7, as well as TGF-β-neutralizing antibodies effectively blocked SMAD2/3 activity. In contrast, an activin-ligand trap was ineffective, indicating that increased autocrine TGF-β signaling is involved. To study potential involvement of MAPK/ERK signaling, cells were treated with a MEK1/2 inhibitor. Surprisingly, fluid flow-induced expression of most SMAD2/3 targets was further enhanced upon MEK inhibition. We conclude that fluid shear stress induces autocrine TGF-β/ALK5-induced target gene expression in renal epithelial cells, which is partially restrained by MEK1/2-mediated signaling.

Next generation risk assessment of chemicals revolves around the use of mechanistic information without animal experimentation. In this regard, toxicogenomics has proven to be a useful tool to elucidate the underlying mechanisms of adverse effects of xenobiotics. In the present study, two widely used human in vitro hepatocyte culture systems, namely primary human hepatocytes (PHH) and human hepatoma HepaRG cells, were exposed to liver toxicants known to induce liver cholestasis, steatosis or necrosis. Benchmark concentration-response modelling was applied to transcriptomics gene co-expression networks (modules) in order to derive benchmark concentrations (BMCs) and to gain mechanistic insight into the hepatotoxic effects. BMCs derived by concentration-response modelling of gene co-expression modules recapitulated concentration-response modelling of individual genes. Although PHH and HepaRG cells showed overlap in deregulated genes and modules by the liver toxicants, PHH demonstrated a higher responsiveness, based on the lower BMCs of co-regulated gene modules. Such BMCs can be used as transcriptomics point of departure (tPOD) for assessing module-associated cellular (stress) pathways/processes. This approach identified clear tPODs of around maximum systemic concentration (Cmax) levels for the tested drugs, while for cosmetics ingredients the BMCs were 10-100 fold higher than the estimated plasma concentrations. This approach could serve next generation risk assessment practice to identify early responsive modules at low BMCs, that could be linked to key events in liver adverse outcome pathways. In turn, this can assist in delineating potential hazards of new test chemicals using in vitro systems and used in a risk assessment when BMCs are paired with chemical exposure assessment.Competing Interest StatementThe authors have declared no competing interest.Footnotes* Updated introduction paragraph on gene co-expression modules and the TXG-MAP tool; extended the method description of the transcriptomics data analysis steps; added PHH cytotoxicity data; added connection of CSA and cholestasis and VPA and steatosis; included or explained often used abbreviations; changed the discussion on BHT metabolism by CYP2B1 for rat and CYP2B6 for human; clarified the discussion on the potential application of module level BMC to determine tPOD and conection to KEs in AOPs; Supplementary figures and tables revised and renumbered when needed.* https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-12668* https://www.ebi.ac.uk/biostudies/arrayexpress/studies/E-MTAB-12677

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/

Nephrotoxicity caused by drug or chemical exposure involves different mechanisms and nephron segments as well as a complex temporal integration of injury and repair responses. Distinct cellular transcriptional programs regulate the time-dependent tissue injury and regeneration responses. Whole kidney transcriptome analysis cannot dissect the complex the nephron segment spatio-temporal injury and regeneration responses. Here, we used laser capture microdissection of formalin-fixed paraffin embedded sections followed by whole genome targeted RNA-sequencing-TempO-Seq and co-expression gene-network (module) analysis to determine the spatial-temporal responses in rat kidney glomeruli (GM), cortical proximal tubules (CPT) and outer-medulla proximal tubules (OMPT) comparison with whole kidney, after a single dose of the nephrotoxicant cisplatin. We demonstrate that cisplatin induced early onset of DNA damage in both CPT and OMPT, but not GM. Sustained DNA damage response was strongest in OMPT coinciding with OMPT specific inflammatory signaling, actin cytoskeletal remodeling and increased glycolytic metabolism coincident with suppression of mitochondrial activity. Later responses reflected regeneration-related cell cycle pathway activation and ribosomal biogenesis in the injured OMPT regions. Activation of modules containing kidney injury biomarkers was strongest in the OMPT, with OMPT Clu expression best correlating with urinary clusterin biomarker measurements compared the correlation of Kim1. Our findings also showed that whole kidney responses were less sensitive than OMPT. In conclusion, our LCM-TempO-Seq method reveals a detailed spatial mechanistic understanding of renal injury/regeneration after nephrotoxicant exposure and identifies the most representative mechanism-based nephron segment specific renal injury biomarkers.Competing Interest StatementThe authors have declared no competing interest.

To minimise unexpected toxicities in early phase clinical studies of new drugs, it is vital to understand fundamental similarities and differences between preclinical test species and humans. We have used physiologically-based pharmacokinetic modelling to identify doses of the model hepatotoxin acetaminophen yielding similar hepatic burdens of the reactive metabolite N-acetyl-p-benzoquinoneimine in mice and rats, to enable comparison of tissue adaptive responses under conditions of equivalent chemical insult. Mice exhibited a greater degree of liver injury than rats, despite the equivalent hepatic NAPQI burden. Transcriptomic and proteomic analyses highlighted the stronger activation of stress response pathways (including the Nrf2 oxidative stress response and autophagy) in the livers of rats. Components of these pathways were also found to be expressed at a higher basal level in the livers of rats compared with both mice and humans. Our findings exemplify a systems approach to understanding differential species sensitivity to hepatotoxicity, and have important implications for species selection and human translation in the safety testing of new drug candidates.

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 comprehensive and mechanisms-revealing data, but analysis tools to interpret mechanisms of toxicity and 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 networks (WGCNA) obtained from the Primary Human Hepatocytes (PHH) TG-GATEs dataset. 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, are perturbed by specific stressors and show preserved in rat systems (rat primary hepatocytes and rat in vivo liver), highlighting stress responses that translate across species/testing systems. The TXG-MAPr tool was successfully applied to investigate the mechanism of toxicity of TG-GATEs compounds and using external datasets obtained from different hepatocyte cells and microarray platforms. Additionally, we suggest that module responses can be calculated from targeted RNA-seq data therefore imputing biological responses from a limited gene. By analyzing 50 different PHH donors’ responses to a common stressor, tunicamycin, we were able to suggest modules associated with donor’s traits, e.g. pre-existing disease state, therefore connected to donors’ variability. 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.