Comparing in vitro human liver models to in vivo human liver using RNA-Seq.

The liver plays an important role in xenobiotic metabolism and represents a primary target for toxic substances. Many different in vitro cell models have been developed in the past decades. In this study, we used RNA-sequencing (RNA-Seq) to analyze the following human in vitro liver cell models in comparison to human liver tissue: cancer-derived cell lines (HepG2, HepaRG 3D), induced pluripotent stem cell-derived hepatocyte-like cells (iPSC-HLCs), cancerous human liver-derived assays (hPCLiS, human precision cut liver slices), non-cancerous human liver-derived assays (PHH, primary human hepatocytes) and 3D liver microtissues. First, using CellNet, we analyzed whether these liver in vitro cell models were indeed classified as liver, based on their baseline expression profile and gene regulatory networks (GRN). More comprehensive analyses using non-differentially expressed genes (non-DEGs) and differential transcript usage (DTU) were applied to assess the coverage for important liver pathways. Through different analyses, we noticed that 3D liver microtissues exhibited a high similarity with in vivo liver, in terms of CellNet (C/T score: 0.98), non-DEGs (10,363) and pathway coverage (highest for 19 out of 20 liver specific pathways shown) at the beginning of the incubation period (0 h) followed by a decrease during long-term incubation for 168 and 336 h. PHH also showed a high degree of similarity with human liver tissue and allowed stable conditions for a short-term cultivation period of 24 h. Using the same metrics, HepG2 cells illustrated the lowest similarity (C/T: 0.51, non-DEGs: 5623, and pathways coverage: least for 7 out of 20) with human liver tissue. The HepG2 are widely used in hepatotoxicity studies, however, due to their lower similarity, they should be used with caution. HepaRG models, iPSC-HLCs, and hPCLiS ranged clearly behind microtissues and PHH but showed higher similarity to human liver tissue than HepG2 cells. In conclusion, this study offers a resource of RNA-Seq data of several biological replicates of human liver cell models in vitro compared to human liver tissue.

Toxicity, recovery, and resilience in a 3D dopaminergic neuronal in vitro model exposed to rotenone.

To date, most in vitro toxicity testing has focused on acute effects of compounds at high concentrations. This testing strategy does not reflect real-life exposures, which might contribute to long-term disease outcome. We used a 3D-human dopaminergic in vitro LUHMES cell line model to determine whether effects of short-term rotenone exposure (100 nM, 24 h) are permanent or reversible. A decrease in complex I activity, ATP, mitochondrial diameter, and neurite outgrowth were observed acutely. After compound removal, complex I activity was still inhibited; however, ATP levels were increased, cells were electrically active and aggregates restored neurite outgrowth integrity and mitochondrial morphology. We identified significant transcriptomic changes after 24 h which were not present 7 days after wash-out. Our results suggest that testing short-term exposures in vitro may capture many acute effects which cells can overcome, missing adaptive processes, and long-term mechanisms. In addition, to study cellular resilience, cells were re-exposed to rotenone after wash-out and recovery period. Pre-exposed cells maintained higher metabolic activity than controls and presented a different expression pattern in genes previously shown to be altered by rotenone. NEF2L2, ATF4, and EAAC1 were downregulated upon single hit on day 14, but unchanged in pre-exposed aggregates. DAT and CASP3 were only altered after re-exposure to rotenone, while TYMS and MLF1IP were downregulated in both single-exposed and pre-exposed aggregates. In summary, our study shows that a human cell-based 3D model can be used to assess cellular adaptation, resilience, and long-term mechanisms relevant to neurodegenerative research.

A structure-activity relationship linking non-planar PCBs to functional deficits of neural crest cells: new roles for connexins.

Migration of neural crest cells (NCC) is a fundamental developmental process, and test methods to identify interfering toxicants have been developed. By examining cell function endpoints, as in the 'migration-inhibition of NCC (cMINC)' assay, a large number of toxicity mechanisms and protein targets can be covered. However, the key events that lead to the adverse effects of a given chemical or group of related compounds are hard to elucidate. To address this issue, we explored here, whether the establishment of two overlapping structure-activity relationships (SAR)-linking chemical structure on the one hand to a phenotypic test outcome, and on the other hand to a mechanistic endpoint-was useful as strategy to identify relevant toxicity mechanisms. For this purpose, we chose polychlorinated biphenyls (PCB) as a large group of related, but still toxicologically and physicochemically diverse structures. We obtained concentration-dependent data for 26 PCBs in the cMINC assay. Moreover, the test chemicals were evaluated by a new high-content imaging method for their effect on cellular re-distribution of connexin43 and for their capacity to inhibit gap junctions. Non-planar PCBs inhibited NCC migration. The potency (1-10 µM) correlated with the number of ortho-chlorine substituents; non-ortho-chloro (planar) PCBs were non-toxic. The toxicity to NCC partially correlated with gap junction inhibition, while it fully correlated (p < 0.0004) with connexin43 cellular re-distribution. Thus, our double-SAR strategy revealed a mechanistic step tightly linked to NCC toxicity of PCBs. Connexin43 patterns in NCC may be explored as a new endpoint relevant to developmental toxicity screening.

Integrate mechanistic evidence from new approach methodologies (NAMs) into a read-across assessment to characterise trends in shared mode of action.

Read-across approaches often remain inconclusive as they do not provide sufficient evidence on a common mode of action across the category members. This read-across case study on thirteen, structurally similar, branched aliphatic carboxylic acids investigates the concept of using human-based new approach methods, such as in vitro and in silico models, to demonstrate biological similarity. Five out of the thirteen analogues have preclinical in vivo studies. Three out of them induced lipid accumulation or hypertrophy in preclinical studies with repeated exposure, which leads to the read-across hypothesis that the analogues can potentially induce hepatic steatosis. To confirm the selection of analogues, the expression patterns of the induced differentially expressed genes (DEGs) were analysed in a human liver model. With increasing dose, the expression pattern within the tested analogues got more similar, which serves as a first indication of a common mode of action and suggests differences in the potency of the analogues. Hepatic steatosis is a well-known adverse outcome, for which over 55 adverse outcome pathways have been identified. The resulting adverse outcome pathway (AOP) network, comprised a total 43 MIEs/KEs and enabled the design of an in vitro testing battery. From the AOP network, ten MIEs, early and late KEs were tested to systematically investigate a common mode of action among the grouped compounds. The targeted testing of AOP specific MIE/KEs shows that biological activity in the category decreases with side chain length. A similar trend was evident in measuring liver alterations in zebra fish embryos. However, activation of single MIEs or early KEs at in vivo relevant doses did not necessarily progress to the late KE "lipid accumulation". KEs not related to the read-across hypothesis, testing for example general mitochondrial stress responses in liver cells, showed no trend or biological similarity. Testing scope is a key issue in the design of in vitro test batteries. The Dempster-Shafer decision theory predicted those analogues with in vivo reference data correctly using one human liver model or the CALUX reporter assays. The case study shows that the read-across hypothesis is the key element to designing the testing strategy. In the case of a good mechanistic understanding, an AOP facilitates the selection of reliable human in vitro models to demonstrate a common mode of action. Testing DEGs, MIEs and early KEs served to show biological similarity, whereas the late KEs become important for confirmation, as progression from MIEs to AO is not always guaranteed.

Development of a generic zebrafish embryo PBPK model and application to the developmental toxicity assessment of valproic acid analogs.

In order to better explain, predict, or extrapolate to humans the developmental toxicity effects of chemicals to zebrafish (Danio rerio) embryos, we developed a physiologically-based pharmacokinetic (PBPK) model designed to predict organ concentrations of neutral or ionizable chemicals, up to 120 h post-fertilization. Chemicals' distribution is modeled in the cells, lysosomes, and mitochondria of ten organs of the embryo. The model's partition coefficients are calculated with sub-models using physicochemical properties of the chemicals of interest. The model accounts for organ growth and changes in metabolic clearance with time. We compared ab initio model predictions to data obtained on culture medium and embryo concentrations of valproic acid (VPA) and nine analogs during continuous dosing under the OECD test guideline 236. We further improved the predictions by estimating metabolic clearance and partition coefficients from the data by Bayesian calibration. We also assessed the performance of the model at reproducing data published by Brox et al. (2016) on VPA and 16 other chemicals. We finally compared dose-response relationships calculated for mortality and malformations on the basis of predicted whole embryo concentrations versus those based on nominal water concentrations. The use of target organ concentrations substantially shifted the magnitude of dose-response parameters and the relative toxicity ranking of chemicals studied.