Mapping Interindividual Variability of Toxicodynamics Using High-Throughput Transcriptomics and Primary Human Hepatocytes from Fifty Donors.

BACKGROUND: Understanding the variability across the human population with respect to toxicodynamic responses after exposure to chemicals, such as environmental toxicants or drugs, is essential to define safety factors for risk assessment to protect the entire population. Activation of cellular stress response pathways are early adverse outcome pathway (AOP) key events of chemical-induced toxicity and would elucidate the estimation of population variability of toxicodynamic responses. OBJECTIVES: We aimed to map the variability in cellular stress response activation in a large panel of primary human hepatocyte (PHH) donors to aid in the quantification of toxicodynamic interindividual variability to derive safety uncertainty factors. METHODS: High-throughput transcriptomics of over 8,000 samples in total was performed covering a panel of 50 individual PHH donors upon 8 to 24 h exposure to broad concentration ranges of four different toxicological relevant stimuli: tunicamycin for the unfolded protein response (UPR), diethyl maleate for the oxidative stress response (OSR), cisplatin for the DNA damage response (DDR), and tumor necrosis factor alpha (TNFα) for NF-κB signaling. Using a population mixed-effect framework, the distribution of benchmark concentrations (BMCs) and maximum fold change were modeled to evaluate the influence of PHH donor panel size on the correct estimation of interindividual variability for the various stimuli. RESULTS: Transcriptome mapping allowed the investigation of the interindividual variability in concentration-dependent stress response activation, where the average of BMCs had a maximum difference of 864-, 13-, 13-, and 259-fold between different PHHs for UPR, OSR, DDR, and NF-κB signaling-related genes, respectively. Population modeling revealed that small PHH panel sizes systematically underestimated the variance and gave low probabilities in estimating the correct human population variance. Estimated toxicodynamic variability factors of stress response activation in PHHs based on this dataset ranged between 1.6 and 6.3. DISCUSSION: Overall, by combining high-throughput transcriptomics and population modeling, improved understanding of interindividual variability in chemical-induced activation of toxicity relevant stress pathways across the human population using a large panel of plated cryopreserved PHHs was established, thereby contributing toward increasing the confidence of in vitro-based prediction of adverse responses, in particular hepatotoxicity. https://doi.org/10.1289/EHP11891.

Assessing the transferability and reproducibility of 3D in vitro liver models from primary human multi-cellular microtissues to cell-line based HepG2 spheroids.

To reduce, replace, and refine in vivo testing, there is increasing emphasis on the development of more physiologically relevant in vitro test systems to improve the reliability of non-animal-based methods for hazard assessment. When developing new approach methodologies, it is important to standardize the protocols and demonstrate the methods can be reproduced by multiple laboratories. The aim of this study was to assess the transferability and reproducibility of two advanced in vitro liver models, the Primary Human multicellular microtissue liver model (PHH) and the 3D HepG2 Spheroid Model, for nanomaterial (NM) and chemical hazard assessment purposes. The PHH model inter-laboratory trial showed strong consistency across the testing sites. All laboratories evaluated cytokine release and cytotoxicity following exposure to titanium dioxide (TiO2) and zinc oxide (ZnO) nanoparticles. No significant difference was observed in cytotoxicity or IL-8 release for the test materials. The data were reproducible with all three laboratories with control readouts within a similar range. The PHH model ZnO induced the greatest cytotoxicity response at 50.0 µg/mL and a dose-dependent increase in IL-8 release. For the 3D HepG2 spheroid model, all test sites were able to construct the model and demonstrated good concordance in IL-8 cytokine release and genotoxicity data. This trial demonstrates the successful transfer of new approach methodologies across multiple laboratories, with good reproducibility for several hazard endpoints.

Model-based translation of DNA damage signaling dynamics across cell types.

Interindividual variability in DNA damage response (DDR) dynamics may evoke differences in susceptibility to cancer. However, pathway dynamics are often studied in cell lines as alternative to primary cells, disregarding variability. To compare DDR dynamics in the cell line HepG2 with primary human hepatocytes (PHHs), we developed a HepG2-based computational model that describes the dynamics of DDR regulator p53 and targets MDM2, p21 and BTG2. We used this model to generate simulations of virtual PHHs and compared the results to those for PHH donor samples. Correlations between baseline p53 and p21 or BTG2 mRNA expression in the absence and presence of DNA damage for HepG2-derived virtual samples matched the moderately positive correlations observed for 50 PHH donor samples, but not the negative correlations between p53 and its inhibitor MDM2. Model parameter manipulation that affected p53 or MDM2 dynamics was not sufficient to accurately explain the negative correlation between these genes. Thus, extrapolation from HepG2 to PHH can be done for some DDR elements, yet our analysis also reveals a knowledge gap within p53 pathway regulation, which makes such extrapolation inaccurate for the regulator MDM2. This illustrates the relevance of studying pathway dynamics in addition to gene expression comparisons to allow reliable translation of cellular responses from cell lines to primary cells. Overall, with our approach we show that dynamical modeling can be used to improve our understanding of the sources of interindividual variability of pathway dynamics.

Towards an advanced testing strategy for genotoxicity using image-based 2D and 3D HepG2 DNA damage response fluorescent protein reporters.

In vitro assessment of mutagenicity is an essential component in the chemical risk assessment. Given the diverse modes of action by which chemicals can induce DNA damage, it is essential that these in vitro assays are carefully evaluated for their possibilities and limitations. In this study, we used a fluorescent protein HepG2 reporter test system in combination with high content imaging. To measure induction of the DNA damage response (DDR), we used three different green fluorescent protein reporters for p53 pathway activation. These allowed for accurate quantification of p53, p21 and BTG2 (BTG anti-proliferation factor 2) protein expression and cell viability parameters at a single cell or spheroid resolution. The reporter lines were cultured as 2D monolayers and as 3D spheroids. Furthermore, liver maturity and cytochrome P450 enzyme expression were increased by culturing in an amino acid-rich (AAGLY) medium. We found that culture conditions that support a sustained proliferative state (2D culturing with normal DMEM medium) give superior sensitivity when genotoxic compounds are tested that do not require metabolisation and of which the mutagenic mode of action is dependent on replication. For compounds, which are metabolically converted to mutagenic metabolites, more differentiated HepG2 DDR reporters (e.g. 3D cultures) showed a higher sensitivity. This study stratifies how different culture methods of HepG2 DDR reporter cells can influence the sensitivity towards diverse genotoxicants and how this provides opportunities for a tiered genotoxicity testing strategy.

Recapitulating Cholangiopathy-Associated Necroptotic Cell Death In Vitro Using Human Cholangiocyte Organoids.

BACKGROUND & AIMS: Liver and bile duct diseases often are associated with extensive cell death of cholangiocytes. Necroptosis represents a common mode of programmed cell death in cholangiopathy, however, detailed mechanistic knowledge is limited owing to the lack of appropriate in vitro models. To address this void, we investigated whether human intrahepatic cholangiocyte organoids (ICOs) can recapitulate cholangiopathy-associated necroptosis and whether this model can be used for drug screening. METHODS: We evaluated the clinical relevance of necroptosis in end-stage liver diseases and liver transplantation by immunohistochemistry. Cholangiopathy-associated programmed cell death was evoked in ICOs derived from healthy donors or patients with primary sclerosing cholangitis or alcoholic liver diseases by the various stimuli. RESULTS: The expression of key necroptosis mediators, receptor-interacting protein 3 and phosphorylated mixed lineage kinase domain-like, in cholangiocytes during end-stage liver diseases was confirmed. The phosphorylated mixed lineage kinase domain-like expression was etiology-dependent. Gene expression analysis confirmed that primary cholangiocytes are more prone to necroptosis compared with primary hepatocytes. Both apoptosis and necroptosis could be specifically evoked using tumor necrosis factor α and second mitochondrial-derived activator of caspases mimetic, with or without caspase inhibition in healthy and patient-derived ICOs. Necroptosis also was induced by ethanol metabolites or human bile in ICOs from donors and patients. The organoid cultures further uncovered interdonor variable and species-specific drug responses. Dabrafenib was identified as a potent necroptosis inhibitor and showed a protective effect against ethanol metabolite toxicity. CONCLUSIONS: Human ICOs recapitulate cholangiopathy-associated necroptosis and represent a useful in vitro platform for the study of biliary cytotoxicity and preclinical drug evaluation.

The human hepatocyte TXG-MAPr: gene co-expression network modules to support mechanism-based risk 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.

Fluorescent tagging of endogenous Heme oxygenase-1 in human induced pluripotent stem cells for high content imaging of oxidative stress in various differentiated lineages.

Tagging of endogenous stress response genes can provide valuable in vitro models for chemical safety assessment. Here, we present the generation and application of a fluorescent human induced pluripotent stem cell (hiPSC) reporter line for Heme oxygenase-1 (HMOX1), which is considered a sensitive and reliable biomarker for the oxidative stress response. CRISPR/Cas9 technology was used to insert an enhanced green fluorescent protein (eGFP) at the C-terminal end of the endogenous HMOX1 gene. Individual clones were selected and extensively characterized to confirm precise editing and retained stem cell properties. Bardoxolone-methyl (CDDO-Me) induced oxidative stress caused similarly increased expression of both the wild-type and eGFP-tagged HMOX1 at the mRNA and protein level. Fluorescently tagged hiPSC-derived proximal tubule-like, hepatocyte-like, cardiomyocyte-like and neuron-like progenies were treated with CDDO-Me (5.62-1000 nM) or diethyl maleate (5.62-1000 µM) for 24 h and 72 h. Multi-lineage oxidative stress responses were assessed through transcriptomics analysis, and HMOX1-eGFP reporter expression was carefully monitored using live-cell confocal imaging. We found that eGFP intensity increased in a dose-dependent manner with dynamics varying amongst lineages and stressors. Point of departure modelling further captured the specific lineage sensitivities towards oxidative stress. We anticipate that the newly developed HMOX1 hiPSC reporter will become a valuable tool in understanding and quantifying critical target organ cell-specific oxidative stress responses induced by (newly developed) chemical entities.

Integration of temporal single cell cellular stress response activity with logic-ODE modeling reveals activation of ATF4-CHOP axis as a critical predictor of drug-induced liver injury.

Drug-induced liver injury (DILI) is the most prevalent adversity encountered in drug development and clinical settings leading to urgent needs to understand the underlying mechanisms. In this study, we have systematically investigated the dynamics of the activation of cellular stress response pathways and cell death outcomes upon exposure of a panel of liver toxicants using live cell imaging of fluorescent reporter cell lines. We established a comprehensive temporal dynamic response profile of a large set of BAC-GFP HepG2 cell lines representing the following components of stress signaling: i) unfolded protein response (UPR) [ATF4, XBP1, BIP and CHOP]; ii) oxidative stress [NRF2, SRXN1, HMOX1]; iii) DNA damage [P53, P21, BTG2, MDM2]; and iv) NF-κB pathway [A20, ICAM1]. We quantified the single cell GFP expression as a surrogate for endogenous protein expression using live cell imaging over > 60 h upon exposure to 14 DILI compounds at multiple concentrations. Using logic-based ordinary differential equation (Logic-ODE), we modelled the dynamic profiles of the different stress responses and extracted specific descriptors potentially predicting the progressive outcomes. We identified the activation of ATF4-CHOP axis of the UPR as the key pathway showing the highest correlation with cell death upon DILI compound perturbation. Knocking down main components of the UPR provided partial protection from compound-induced cytotoxicity, indicating a complex interplay among UPR components as well as other stress pathways. Our results suggest that a systematic analysis of the temporal dynamics of ATF4-CHOP axis activation can support the identification of DILI risk for new candidate drugs.

Systematic transcriptome-based comparison of cellular adaptive stress response activation networks in hepatic stem cell-derived progeny and primary human hepatocytes.

Various adaptive cellular stress response pathways are critical in the pathophysiology of liver disease and drug-induced liver injury. Human-induced pluripotent stem cell (hiPSC)-derived hepatocyte-like cells (HLCs) provide a promising tool to study cellular stress response pathways, but in this context there is limited insight on how HLCs compare to other in vitro liver models. Here, we systematically compared the transcriptomic profiles upon chemical activation in HLCs, hiPSC, primary human hepatocytes (PHH) and HepG2 liver cancer cells. We used targeted RNA-sequencing to map concentration transcriptional response using benchmark concentration modeling for the various stress responses in the different test systems. We found that HLCs are very sensitive towards oxidative stress and inflammation conditions as corresponding genes were activated at over 3 fold lower concentrations of the corresponding pathway inducing compounds as compared to PHH. PHH were the most sensitive model when studying UPR related effects. Due to the non-proliferative nature of PHH and HLCs, these do not pose a good/sensitive model to pick up DNA damage responses, while hiPSC and HepG2 were more sensitive in these conditions. We envision that this study contributes to a better understanding on how HLCs can contribute to the assessment of cell physiological stress response activation to predict hepatotoxic events.

In Vitro Three-Dimensional Liver Models for Nanomaterial DNA Damage Assessment.

Whilst the liver possesses the ability to repair and restore sections of damaged tissue following acute injury, prolonged exposure to engineered nanomaterials (ENM) may induce repetitive injury leading to chronic liver disease. Screening ENM cytotoxicity using 3D liver models has recently been performed, but a significant challenge has been the application of such in vitro models for evaluating ENM associated genotoxicity; a vital component of regulatory human health risk assessment. This review considers the benefits, limitations, and adaptations of specific in vitro approaches to assess DNA damage in the liver, whilst identifying critical advancements required to support a multitude of biochemical endpoints, focusing on nano(geno)toxicology (e.g., secondary genotoxicity, DNA damage, and repair following prolonged or repeated exposures).

Utility of Extrapolating Human S1500+ Genes to the Whole Transcriptome: Tunicamycin Case Study.

The TempO-Seq S1500+ platform(s), now available for human, mouse, rat, and zebrafish, measures a discrete number of genes that are representative of biological and pathway co-regulation across the entire genome in a given species. While measurement of these genes alone provides a direct assessment of gene expression activity, extrapolating expression values to the whole transcriptome (~26 000 genes in humans) can estimate measurements of non-measured genes of interest and increases the power of pathway analysis algorithms by using a larger background gene expression space. Here, we use data from primary hepatocytes of 54 donors that were treated with the endoplasmic reticulum (ER) stress inducer tunicamycin and then measured on the human S1500+ platform containing ~3000 representative genes. Measurements for the S1500+ genes were then used to extrapolate expression values for the remaining human transcriptome. As a case study of the improved downstream analysis achieved by extrapolation, the "measured only" and "whole transcriptome" (measured + extrapolated) gene sets were compared. Extrapolation increased the number of significant genes by 49%, bringing to the forefront many that are known to be associated with tunicamycin exposure. The extrapolation procedure also correctly identified established tunicamycin-related functional pathways reflected by coordinated changes in interrelated genes while maintaining the sample variability observed from the "measured only" genes. Extrapolation improved the gene- and pathway-level biological interpretations for a variety of downstream applications, including differential expression analysis, gene set enrichment pathway analysis, DAVID keyword analysis, Ingenuity Pathway Analysis, and NextBio correlated compound analysis. The extrapolated data highlight the role of metabolism/metabolic pathways, the ER, immune response, and the unfolded protein response, each of which are key activities associated with tunicamycin exposure that were unrepresented or underrepresented in one or more of the analyses of the original "measured only" dataset. Furthermore, the inclusion of the extrapolated genes raised "tunicamycin" from third to first upstream regulator in Ingenuity Pathway Analysis and from sixth to second most correlated compound in NextBio analysis. Therefore, our case study suggests an approach to extend and enhance data from the S1500+ platform for improved insight into biological mechanisms and functional outcomes of diseases, drugs, and other perturbations.

From fatty hepatocytes to impaired bile flow: Matching model systems for liver biology and disease.

A large variety of model systems are used in hepatobiliary research. In this review, we aim to provide an overview of established and emerging models for specific research questions. We specifically discuss the value and limitations of these models for research on metabolic associated fatty liver disease (MAFLD), (previously named non-alcoholic fatty liver diseases/non-alcoholic steatohepatitis (NAFLD/NASH)) and cholestasis-related diseases such as primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC). The entire range of models is discussed varying from immortalized cell lines, mature or pluripotent stem cell-based models including organoids/spheroids, to animal models and human ex vivo models such as normothermic machine perfusion of livers and living liver slices. Finally, the pros and cons of each model are discussed as well as the need in the scientific community for continuous innovation in model development to better mimic the human (patho)physiology.

ATF6 Is a Critical Determinant of CHOP Dynamics during the Unfolded Protein Response.

The unfolded protein response (UPR) pathway senses unfolded proteins and regulates proteostasis and cell fate through activity of the transcription factors ATF4, ATF6, and XBP1 within a complex network of three main branches. Here, we investigated contributions of the three branches to UPR activity in single cells using microscopy-based quantification and dynamic modeling. BAC-GFP HepG2 reporter cell lines were exposed to tunicamycin, and activation of various UPR components was monitored for 24 h. We constructed a dynamic model to describe the adaptive UPR signaling network, for which incorporation of all three branches was required to match the data. Our calibrated model suggested that ATF6 shapes the early dynamics of pro-apoptotic CHOP. We confirmed this hypothesis by measurements beyond 24 h, by perturbing single siRNA knockdowns and by ATF6 measurements. Overall, our work indicates that ATF6 is an important regulator of CHOP, which in turn regulates cell fate decisions.

System Microscopy of Stress Response Pathways in Cholestasis Research.

Exposure to oxidative radical species and cytokine-mediated inflammatory stress are established contributors to hepatocyte cell death during cholestasis. Cellular counter measures against those stressors are called adaptive stress response pathways. While in early stages of the disease adaptive stress pathways protect the hepatocytes, in later stages during prolonged stressed conditions they fail. The quantitative imaging-based assessment of cellular stress response pathways using the HepG2 BAC-GFP response reporter platform is a powerful strategy to evaluate the impact of chemical substances and gene knockdown on activation of adaptive stress response pathways, hence allowing systematic screening for positive or negative influences on cholestasis progression. This protocol allows the application of a highly versatile screening tool for a systematic evaluation of the effect of compounds having cholestasis liability and affected genes during cholestatic injury on cellular adaptive stress pathway activation. The approach involves high-throughput live-cell visualization of GFP-tagged key proteins of the oxidative stress response/Nrf2 pathway and inflammatory cytokine signaling. Quantitative image analysis of temporal responses of individual cells is followed by informatics analysis. The overall practical approaches are discussed in this chapter.

Preclinical imaging methods for assessing the safety and efficacy of regenerative medicine therapies.

Regenerative medicine therapies hold enormous potential for a variety of currently incurable conditions with high unmet clinical need. Most progress in this field to date has been achieved with cell-based regenerative medicine therapies, with over a thousand clinical trials performed up to 2015. However, lack of adequate safety and efficacy data is currently limiting wider uptake of these therapies. To facilitate clinical translation, non-invasive in vivo imaging technologies that enable careful evaluation and characterisation of the administered cells and their effects on host tissues are critically required to evaluate their safety and efficacy in relevant preclinical models. This article reviews the most common imaging technologies available and how they can be applied to regenerative medicine research. We cover details of how each technology works, which cell labels are most appropriate for different applications, and the value of multi-modal imaging approaches to gain a comprehensive understanding of the responses to cell therapy in vivo.

Comprehensive Landscape of Nrf2 and p53 Pathway Activation Dynamics by Oxidative Stress and DNA Damage.

A quantitative dynamics pathway map of the Nrf2-mediated oxidative stress response and p53-related DNA damage response pathways as well as the cross-talk between these pathways has not systematically been defined. To allow the dynamic single cell evaluation of these pathways, we have used BAC-GFP recombineering to tag for each pathway's three key components: for the oxidative stress response, Keap1-GFP, Nrf2-GFP, and Srxn1-GFP; for the DNA damage response, 53bp1-GFP, p53-GFP, and p21-GFP. The dynamic activation of these individual components was assessed using quantitative high throughput confocal microscopy after treatment with a broad concentration range of diethyl maleate (DEM; to induce oxidative stress) and etoposide (to induce DNA damage). DEM caused a rapid activation of Nrf2, which returned to baseline levels at low concentrations but remained sustained at high concentrations. Srxn1-GFP induction and Keap1-GFP translocation to autophagosomes followed later, with upper boundaries reached at high concentrations, close to the onset of cell death. Etoposide caused rapid accumulation of 53bp1-GFP in DNA damage foci, which was later followed by the concentration dependent nuclear accumulation of p53-GFP and subsequent induction of p21-GFP. While etoposide caused activation of Srxn1-GFP, a modest activation of DNA damage reporters was observed for DEM at high concentrations. Interestingly, Nrf2 knockdown caused an inhibition of the DNA damage response at high concentrations of etoposide, while Keap1 knockdown caused an enhancement of the DNA damage response already at low concentrations of etoposide. Knockdown of p53 did not affect the oxidative stress response. Altogether, the current stress response landscapes provide insight in the time course responses of and cross-talk between oxidative stress and DNA-damage and defines the tipping points where cell injury may switch from adaptation to injury.

Quantitative high content imaging of cellular adaptive stress response pathways in toxicity for chemical safety assessment.

Over the past decade, major leaps forward have been made on the mechanistic understanding and identification of adaptive stress response landscapes underlying toxic insult using transcriptomics approaches. However, for predictive purposes of adverse outcome several major limitations in these approaches exist. First, the limited number of samples that can be analyzed reduces the in depth analysis of concentration-time course relationships for toxic stress responses. Second these transcriptomics analysis have been based on the whole cell population, thereby inevitably preventing single cell analysis. Third, transcriptomics is based on the transcript level, totally ignoring (post)translational regulation. We believe these limitations are circumvented with the application of high content analysis of relevant toxicant-induced adaptive stress signaling pathways using bacterial artificial chromosome (BAC) green fluorescent protein (GFP) reporter cell-based assays. The goal is to establish a platform that incorporates all adaptive stress pathways that are relevant for toxicity, with a focus on drug-induced liver injury. In addition, cellular stress responses typically follow cell perturbations at the subcellular organelle level. Therefore, we complement our reporter line panel with reporters for specific organelle morphometry and function. Here, we review the approaches of high content imaging of cellular adaptive stress responses to chemicals and the application in the mechanistic understanding and prediction of chemical toxicity at a systems toxicology level.