A network-based transcriptomic landscape of HepG2 cells uncovering causal gene-cytotoxicity interactions underlying drug-induced liver injury

Drug-induced liver injury (DILI) remains the main reason for drug development attritions largely due to poor mechanistic understanding. Toxicogenomic to interrogate the mechanism of DILI has been broadly performed. Gene co-regulation network-based transcriptome analysis is a bioinformatics approach that potentially contributes to improve mechanistic interpretation of toxicogenomic data. Here we performed an extensive concentration time course response-toxicogenomic study in the HepG2 cell line exposed to 20 DILI compounds, 7 reference compounds for stress response pathways, and 10 agonists for cytokines and growth factor receptors. We performed whole transcriptome targeted RNA sequencing to more than 500 conditions to and applied weighted gene co-regulated network analysis (WGCNA) to the transcriptomics data followed by identification of gene co-regulated networks (modules) that were strongly modulated upon the exposure of DILI compounds. Preservation analysis on the module responses of HepG2 and PHH demonstrated highly preserved adaptive stress response gene co-regulated networks. We correlated gene co-regulated networks with cell death onset and causal relationships of 67 critical target genes of these modules with onset of cell death was evaluated using RNA interference screening. We identified GTPBP2, HSPA1B, IRF1, SIRT1 and TSC22D3 as essential modulators of DILI compound-induced cell death. These genes were also induced by DILI compounds in PHH. Altogether, we demonstrate the application of large transcriptome datasets combined with network-based analysis and biological validation to uncover the candidate determinants of DILI.