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

ECM-derived biomaterials for regulating tissue multicellularity and maturation

Recent breakthroughs in developing human-relevant organotypic models led to the building of highly resemblant tissue constructs that hold immense potential for transplantation, drug screening, and disease modeling. Despite the progress in fine-tuning stem cell multilineage differentiation in highly controlled spatiotemporal conditions and hosting microenvironments, 3D models still experience naive and incomplete morphogenesis. In particular, existing systems and induction protocols fail to maintain stem cell long-term potency, induce high tissue-level multicellularity, or drive the maturity of stem cell-derived 3D models to levels seen in their in vivo counterparts. In this review, we highlight the use of extracellular matrix (ECM)-derived biomaterials in providing stem cell niche-mimicking microenvironment capable of preserving stem cell long-term potency and inducing spatial and region-specific differentiation. We also examine the maturation of different 3D models, including organoids, encapsulated in ECM biomaterials and provide looking-forward perspectives on employing ECM biomaterials in building more innovative, transplantable, and functional organs.

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.

Capturing time-dependent activation of genes and stress-response pathways using transcriptomics in iPSC-derived renal proximal tubule cells

Transcriptomic analysis is a powerful method in the utilization of New Approach Methods (NAMs) for identifying mechanisms of toxicity and application to hazard characterization. With this regard, mapping toxicological events to time of exposure would be helpful to characterize early events. Here, we investigated time-dependent changes in gene expression levels in iPSC-derived renal proximal tubular-like cells (PTL) treated with five diverse compounds using TempO-Seq transcriptomics with the aims to evaluate the application of PTL for toxicity prediction and to report on temporal effects for the activation of cellular stress response pathways. PTL were treated with either 50 μM amiodarone, 10 μM sodium arsenate, 5 nM rotenone, or 300 nM tunicamycin over a temporal time course between 1 and 24 h. The TGFβ-type I receptor kinase inhibitor GW788388 (1 μM) was used as a negative control. Pathway analysis revealed the induction of key stress-response pathways, including Nrf2 oxidative stress response, unfolding protein response, and metal stress response. Early response genes per pathway were identified much earlier than 24 h and included HMOX1, ATF3, DDIT3, and several MT1 isotypes. GW788388 did not induce any genes within the stress response pathways above, but showed deregulation of genes involved in TGFβ inhibition, including downregulation of CYP24A1 and SERPINE1 and upregulation of WT1. This study highlights the application of iPSC-derived renal cells for prediction of cellular toxicity and sheds new light on the temporal and early effects of key genes that are involved in cellular stress response pathways.

Analyses of Transcriptomics Cell Signalling for Pre-Screening Applications in the Integrated Approach for Testing and Assessment of Non-Genotoxic Carcinogens

With recent rapid advancement of methodological tools, mechanistic understanding of biological processes leading to carcinogenesis is expanding. New approach methodologies such as transcriptomics can inform on non-genotoxic mechanisms of chemical carcinogens and can be developed for regulatory applications. The Organisation for the Economic Cooperation and Development (OECD) expert group developing an Integrated Approach to the Testing and Assessment (IATA) of Non-Genotoxic Carcinogens (NGTxC) is reviewing the possible assays to be integrated therein. In this context, we review the application of transcriptomics approaches suitable for pre-screening gene expression changes associated with phenotypic alterations that underlie the carcinogenic processes for subsequent prioritisation of downstream test methods appropriate to specific key events of non-genotoxic carcinogenesis. Using case studies, we evaluate the potential of gene expression analyses especially in relation to breast cancer, to identify the most relevant approaches that could be utilised as (pre-) screening tools, for example Gene Set Enrichment Analysis (GSEA). We also consider how to address the challenges to integrate gene panels and transcriptomic assays into the IATA, highlighting the pivotal omics markers identified for assay measurement in the IATA key events of inflammation, immune response, mitogenic signalling and cell injury.

A comprehensive transcriptomic comparison of hepatocyte model systems improves selection of models for experimental use

The myriad of available hepatocyte in vitro models provides researchers the possibility to select hepatocyte-like cells (HLCs) for specific research goals. However, direct comparison of hepatocyte models is currently challenging. We systematically searched the literature and compared different HLCs, but reported functions were limited to a small subset of hepatic functions. To enable a more comprehensive comparison, we developed an algorithm to compare transcriptomic data across studies that tested HLCs derived from hepatocytes, biliary cells, fibroblasts, and pluripotent stem cells, alongside primary human hepatocytes (PHHs). This revealed that no HLC covered the complete hepatic transcriptome, highlighting the importance of HLC selection. HLCs derived from hepatocytes had the highest transcriptional resemblance to PHHs regardless of the protocol, whereas the quality of fibroblasts and PSC derived HLCs varied depending on the protocol used. Finally, we developed and validated a web application (HLCompR) enabling comparison for specific pathways and addition of new HLCs. In conclusion, our comprehensive transcriptomic comparison of HLCs allows selection of HLCs for specific research questions and can guide improvements in culturing conditions.

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.

Susceptibility to develop cancer varies among people, partially due to differences in genetic background. Ideally, healthy human-derived cells are used to investigate intracellular signaling pathways and their interindividual variability contributing to cancer susceptibility. Because cells from healthy human tissue are difficult to obtain and culture for periods longer than a few days, cell lines are often used as substitute. However, it is unclear to what extent signaling dynamics in cell lines represent dynamics in healthy human tissue. We asked whether we could reproduce interindividual variability in DNA damage response gene expression in a set of 50 human liver cell donors. Therefore, we built a mathematical model that simulates temporal expression dynamics of the DNA damage response in the HepG2 liver cell line upon chemical activation and used the simulations to create virtual donors. Our virtual donors displayed similar relations between genes as the samples from human donors, provided that we adjusted the strength of specific molecular interactions. Thus, our approach can be used to examine the applicability of widely used cell systems to healthy human tissue in terms of their dynamic responses.

Mapping the dynamics of Nrf2 antioxidant and NFκB inflammatory responses by soft electrophilic chemicals in human liver cells defines the transition from adaptive to adverse responses

A comprehensive understanding of the dynamic activation and crosstalk between different cellular stress response pathways that drive cell adversity is crucial in chemical safety assessment. Various chemicals have electrophilic properties that drive cell injury responses in particular oxidative stress signaling and inflammatory signaling. Here we used bacterial artificial chromosome-based GFP cellular stress reporters with live cell confocal imaging, to systematically monitor the differential modulation of the dynamics of stress pathway activation by six different soft electrophiles: sulforaphane, andrographolide, diethyl maleate, CDDO-Me, ethacrynic acid and tert-butyl hydroquinone. The various soft electrophiles showed differential potency and dynamics of Nrf2 activation and nuclear translocation. These differences in Nrf2 dynamics correlated with distinct activation pattern of Nrf2 downstream targets SRNX1 and HMOX1. All soft electrophiles caused a strong dose dependent suppression of a cytokine-induced NFĸB response represented by suppression of NFĸB nuclear oscillation and inhibition of the downstream target gene activation A20 and ICAM1, which followed the potency of Nrf2 modulation but occurred at higher concentration close to saturation of Nrf2 activation. RNAi-based depletion of RelA resulted in a prolonged presence of Nrf2 in the nucleus after soft electrophile treatment; depletion of Nrf2 caused the induction of NFĸB signaling and activation of its downstream targets A20 and ICAM1. A systematic transcriptome analysis confirmed these effects by soft electrophiles on Nrf2 and NFκB signaling crosstalk in human induced-pluripotent stem cell-derived hepatocyte-like cells. Altogether our data indicate that modulation of Nrf2 by soft electrophiles may have consequences for efficient inflammatory signaling.

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.