Gene networks and transcription factor motifs defining the differentiation of stem cells into hepatocyte-like cells

Background & AimsThe differentiation of stem cells to hepatocyte-like cells (HLC) offers the perspective of unlimited supply of human hepatocytes. However, the degree of differentiation of HLC remains controversial. To obtain an unbiased characterization, we performed a transcriptomic study with HLC derived from human embryonic and induced stem cells (ESC, hiPSC) from three different laboratories.MethodsGenome-wide gene expression profiles of ESC and HLC were compared to freshly isolated and up to 14days cultivated primary human hepatocytes. Gene networks representing successful and failed hepatocyte differentiation, and the transcription factors involved in their regulation were identified.ResultsGene regulatory network analysis demonstrated that HLC represent a mixed cell type with features of liver, intestine, fibroblast and stem cells. The “unwanted” intestinal features were associated with KLF5 and CDX2 transcriptional networks. Cluster analysis identified highly correlated groups of genes associated with mature liver functions (n=1057) and downregulated proliferation associated genes (n=1562) that approach levels of primary hepatocytes. However, three further clusters containing 447, 101, and 505 genes failed to reach levels of hepatocytes. Key TF of two of these clusters include SOX11, FOXQ1, and YBX3. The third unsuccessful cluster, controlled by HNF1, CAR, FXR, and PXR, strongly overlaps with genes repressed in cultivated hepatocytes compared to freshly isolated hepatocytes, suggesting that current in vitro conditions lack stimuli required to maintain gene expression in hepatocytes, which consequently also explains a corresponding deficiency of HLC.ConclusionsThe present gene regulatory network approach identifies key transcription factors which require modulation to improve HLC differentiation

Prediction of human drug-induced liver injury (DILI) in relation to oral doses and blood concentrations

Drug-induced liver injury (DILI) cannot be accurately predicted by animal models. In addition, currently available in vitro methods do not allow for the estimation of hepatotoxic doses or the determination of an acceptable daily intake (ADI). To overcome this limitation, an in vitro/in silico method was established that predicts the risk of human DILI in relation to oral doses and blood concentrations. This method can be used to estimate DILI risk if the maximal blood concentration (Cmax) of the test compound is known. Moreover, an ADI can be estimated even for compounds without information on blood concentrations. To systematically optimize the in vitro system, two novel test performance metrics were introduced, the toxicity separation index (TSI) which quantifies how well a test differentiates between hepatotoxic and non-hepatotoxic compounds, and the toxicity estimation index (TEI) which measures how well hepatotoxic blood concentrations in vivo can be estimated. In vitro test performance was optimized for a training set of 28 compounds, based on TSI and TEI, demonstrating that (1) concentrations where cytotoxicity first becomes evident in vitro (EC10) yielded better metrics than higher toxicity thresholds (EC50); (2) compound incubation for 48 h was better than 24 h, with no further improvement of TSI after 7 days incubation; (3) metrics were moderately improved by adding gene expression to the test battery; (4) evaluation of pharmacokinetic parameters demonstrated that total blood compound concentrations and the 95%-population-based percentile of Cmax were best suited to estimate human toxicity. With a support vector machine-based classifier, using EC10 and Cmax as variables, the cross-validated sensitivity, specificity and accuracy for hepatotoxicity prediction were 100, 88 and 93%, respectively. Concentrations in the culture medium allowed extrapolation to blood concentrations in vivo that are associated with a specific probability of hepatotoxicity and the corresponding oral doses were obtained by reverse modeling. Application of this in vitro/in silico method to the rat hepatotoxicant pulegone resulted in an ADI that was similar to values previously established based on animal experiments. In conclusion, the proposed method links oral doses and blood concentrations of test compounds to the probability of hepatotoxicity

Detoxication Strategy of Epoxide Hydrolase—The Basis for a Novel Threshold for Definable Genotoxic Carcinogens

From our recent work on the three-dimensional structure of epoxide hydrolases we theoretically deduced the likelihood of a two-step catalytic mechanism that we and others have subsequently experimentally confirmed. Analysis of the rate of the two steps by us and by others show that the first step—responsible for removal of the reactive epoxide from the system—works extraordinarily fast (typically three orders of magnitude faster than the second step), sucking up the epoxide like a sponge. Regeneration of the free enzyme (the second step of the catalytic mechanism) is slow. This becomes a toxicological problem only at doses of the epoxide that titrate the enzyme out. Our genotoxicity work shows that indeed this generates a practical threshold below which no genotoxicity is observed. This shows that—contrary to old dogma—practical thresholds exist for definable genotoxic carcinogens

Sequestration of biological reactive intermediates by trapping as covalent enzyme-intermediate complex

One important class of biological reactive intermediates arising in the course of human xenobiotic metabolism are arene and alkene oxides. The major safeguard against the potential genotoxic effects of these compounds is the microsomal epoxide hydrolase (mEH). This enzyme has a broad substrate specificity but--on the first sight--seems to be inadequately suited for this protection task due to its low turnover number with most of its substrates. The recent progress in the understanding of the mechanism of enzymatic epoxide hydrolysis has shed new light on this apparent dilemma: Epoxide hydrolases convert their substrates via the intermediate formation of a covalent enzyme-substrate complex, and it has been shown that the formation of the intermediate proceeds by orders of magnitudes faster than the subsequent hydrolysis, i.e. the formation of the terminal product. Thus, the enzyme acts like a molecular sponge by binding and inactivating the dangerous metabolite very fast while the subsequent product release is considerably slower, and quantification of the latter heavily underestimates the speed of detoxification. Usually, the slow enzyme regeneration does not pose a problem, since the mEH is highly abundant in human liver, the organ with the highest capacity to metabolically generate epoxides. Computer simulation provides evidence that the high amount of mEH enzyme is crucial for the control of the steady-state level of a substrate epoxide and can keep it extremely low. Once the mEH is titrated out under conditions of extraordinarily high epoxide concentration, the epoxide steady-state level steeply rises, leading to a sudden burst of the genotoxic effect. This prediction of the computer simulation is in perfect agreement with our experimental work. V79 Chinese Hamster cells that we have genetically engineered to express human mEH at about the same level as that observed in human liver are well protected from any measurable genotoxic effect of the model compound styrene oxide (STO) up to an apparent threshold level of 100 microM in the cell culture medium. In V79 cells that do not express mEH, STO triggers the formation of DNA strand breaks in a dose-dependent manner with no apparent threshold. Above 100 microM, the genotoxic effect of STO in the mEH-expressing cell line parallels the one in the parental cell line

A liver digital twin for in silico testing of cellular and inter-cellular mechanisms in regeneration after drug-induced damage

Summary: This communication presents a mathematical mechanism-based model of the regenerating liver after drug-induced pericentral lobule damage resolving tissue microarchitecture. The consequence of alternative hypotheses about the interplay of different cell types on regeneration was simulated. Regeneration dynamics has been quantified by the size of the damage-induced dead cell area, the hepatocyte density and the spatial-temporal profile of the different cell types. We use deviations of observed trajectories from the simulated system to identify branching points, at which the systems behavior cannot be explained by the underlying set of hypotheses anymore. Our procedure reflects a successful strategy for generating a fully digital liver twin that, among others, permits to test perturbations from the molecular up to the tissue scale. The model simulations are complementing current knowledge on liver regeneration by identifying gaps in mechanistic relationships and guiding the system toward the most informative (lacking) parameters that can be experimentally addressed

Hepatobiliary Thyroid Hormone Deficiency Impacts Bile Acid Hydrophilicity and Aquaporins in Cholestatic C57BL/6J Mice

Women are more prone to develop either hypothyroidism or cholesterol gallstones than men. However, a male predominance in cholesterol gallstones under hypothyroidism was reported. Recently, a novel pathogenic link between thyroid hormone (TH) deficiency and cholesterol gallstones has been described in male mice. Here, we investigate if TH deficiency impacts cholesterol gallstone formation in females by the same mechanism. Three-month-old C57BL/6J mice were randomly divided into a control, a TH deficient, a lithogenic, and a lithogenic + TH deficient group and diet-treated for two, four, and six weeks. Gallstone prevalence, liver function tests, bile composition, hepatic gene expression, and gallbladder aquaporin expression and localization were investigated. Cholesterol gallstones were observed in lithogenic + TH deficient but not lithogenic only female mice. Diminished hydrophilicity of primary bile acids due to decreased gene expression of hepatic detoxification phase II enzymes was observed. A sex-specific expression and localization of hepatobiliary aquaporins involved in transcellular water and glycerol permeability was observed under TH deficient and lithogenic conditions. TH deficiency promotes cholesterol gallstone formation in female C57BL/6J mice by the same mechanism as observed in males. However, cholesterol gallstone prevalence was lower in female than male C57BL/6J mice. Interestingly, the sex-specific expression and localization of hepatobiliary aquaporins could protect female C57BL/6J mice to cholestasis and could reduce biliary water transport in male C57BL/6J mice possibly contributing to the sex-dependent cholesterol gallstone prevalence under TH deficiency

Cell death mechanisms of the anti-cancer drug etoposide on human cardiomyocytes isolated from pluripotent stem cells

Etoposide (ETP) and anthracyclines are applied for wide anti-cancer treatments. However, the ETP-induced cardiotoxicity remains to be a major safety issue and the underlying cardiotoxic mechanisms are not well understood. This study is aiming to unravel the cardiotoxicity profile of ETP in comparison to anthracyclines using physiologically relevant human pluripotent stem cell-derived cardiomyocytes (hPSC-CMs). Using xCELLigence real-time cell analyser (RTCA), we found that single high dose of ETP induces irreversible increase in hPSC-CMs beating rate and decrease in beating amplitude. We also identified 58 deregulated genes consisting of 33 upregulated and 25 downregulated genes in hPSC-CMs after ETP treatment. Gene ontology (GO) and pathway analysis showed that most upregulated genes are enriched in GO categories like positive regulation of apoptotic process, regulation of cell death, and mitochondria organization, whereas most downregulated genes were enriched in GO categories like cytoskeletal organization, muscle contraction, and Ca2+ ion homeostasis. Moreover, we also found upregulation in 5 miRNAs (has-miR-486-3p, has-miR-34c-5p, has-miR-4423-3p, has-miR-182-5p, and has-miR-139-5p) which play role in muscle contraction, arginine and proline metabolism, and hypertrophic cardiomyopathy (HCM). Immunostaining and transmission electron microscopy also confirmed the cytoskeletal and mitochondrial damage in hPSC-CMs treated with ETP, as well as noticeable alterations in intracellular calcium handling and mitochondrial membrane potential were also observed. The apoptosis inhibitor, Pifithrin-alpha, found to protect hPSC-CMs from ETP-induced cardiotoxicity, whereas hPSC-CMs treated with ferroptosis inhibitor, Liproxstatin-1, showed significant recovery in hPSC-CMs functional properties like beating rate and amplitude after ETP treatment. We suggest that the damage to mitochondria is a major contributing factor involved in ETP-induced cardiotoxicity and the activation of the p53-mediated ferroptosis pathway by ETP is likely the critical pathway in ETP-induced cardiotoxicity. We also conclude that the genomic biomarkers identified in this study will significantly contribute to develop and predict potential cardiotoxic effects of novel anti-cancer drugs in vitro

Unique Metabolic Features of Stem Cells, Cardiomyocytes, and Their Progenitors

Recently, growing attention has been directed toward stem cell metabolism, with the key observation that the plasticity of stem cells also reflects the plasticity of their energy substrate metabolism. There seems to be a clear link between the self-renewal state of stem cells, in which cells proliferate without differentiation, and the activity of specific metabolic pathways. Differentiation is accompanied by a shift from anaerobic glycolysis to mitochondrial respiration. This metabolic switch of differentiating stem cells is required to cover the energy demands of the different organ-specific cell types. Among other metabolic signatures, amino acid and carbohydrate metabolism is most prominent in undifferentiated embryonic stem cells, whereas the fatty acid metabolic signature is unique in cardiomyocytes derived from embryonic stem cells. Identifying the specific metabolic pathways involved in pluripotency and differentiation is critical for further progress in the field of developmental biology and regenerative medicine. The recently generated knowledge on metabolic key processes may help to generate mature stem cell-derived somatic cells for therapeutic applications without the requirement of genetic manipulation. In the present review, the literature about metabolic features of stem cells and their cardiovascular cell derivatives as well as the specific metabolic gene signatures differentiating between stem and differentiated cells are summarized and discussed

Functional cardiotoxicity assessment of cosmetic compounds using human-induced pluripotent stem cell-derived cardiomyocytes

There is a large demand of a human relevant in vitro test system suitable for assessing the cardiotoxic potential of cosmetic ingredients and other chemicals. Using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), we have already established an in vitro cardiotoxicity assay and identified genomic biomarkers of anthracycline-induced cardiotoxicity in our previous work. Here, five cosmetic ingredients were studied by the new hiPSC-CMs test; kojic acid (KJA), triclosan (TS), triclocarban (TCC), 2,7-naphthalenediol (NPT), and basic red 51 (BR51) based on cytotoxicity as well as ATP assays, beating rate, and genomic biomarkers to determine the lowest observed effect concentration (LOEC) and no observed effect concentration (NOEC). The LOEC for beating rate were 400, 10, 3, >400, and 3 mu M for KJA, TS, TCC, NPT, and BR51, respectively. The corresponding concentrations for cytotoxicity or ATP depletion were similar, with the exception of TS and TCC, where the cardiomyo-cyte-beating assay showed positive results at non-cytotoxic concentrations. Functional analysis also showed that the individual compounds caused different effects on hiPSC-CMs. While exposure to KJA, TS, TCC, and BR51 induced significant arrhythmic beating, NPT slightly decreased cell viability, but did not influence beating. Gene expression studies showed that TS and NPT caused down-regulation of cytoskeletal and cardiac ion homeostasis genes. Moreover, TS and NPT deregulated genomic biomarkers known to be affected also by anthracyclines. The present study demonstrates that hiPSC-CMs can be used to determine LOECs and NOECs in vitro, which can be compared to human blood concentrations to determine margins of exposure. Our in vitro assay, which so far has been tested with several anthracyclines and cosmetics, still requires validation by larger numbers of positive and negative controls, before it can be recommended for routine analysis