Design principles of concentration-dependent transcriptome deviations in drug-exposed differentiating stem cells

Risk assessment of parabens in a transcriptomics-based in vitro test

Parabens have been used for decades as preservatives in food, drugs and cosmetics. The majority however, were banned in 2009 and 2014 leaving only methyl-, ethyl-, propyl-, and butyl-derivates available for subsequent use. Methyl- and propylparaben have been extensively tested in vivo, with no resulting evidence for developmental and reproductive toxicity (DART). In contrast, ethylparaben has not yet been tested for DART in animal experiments, and it is currently debated if additional animal studies are warranted. In order to perform a comparison of the four currently approved parabens, we used a previously established in vitro test based on human induced pluripotent stem cells (iPSC) that are exposed to test substances during their differentiation to neuroectodermal cells. EC50 values for cytotoxicity were 906 μM, 698 μM, 216 μM and 63 μM for methyl-, ethyl-, propyl- and butylparaben, respectively, demonstrating that cytotoxicity increases with increasing alkyl chain length. Genome-wide analysis demonstrated that FDR-adjusted significant gene expression changes occurred only at cytotoxic or close to cytotoxic concentrations, for example 1720 differentially expressed genes (DEG) at 1000 μM ethylparaben, 1 DEG at 316 μM, and no DEG at 100 μM or lower concentrations. The highest concentration of ethylparaben that did not induce any cytotoxicity nor DEG was 1670-fold above the highest concentration reported in biomonitoring studies (60 nM ethylparaben in cord blood). In conclusion, cytotoxicity and gene expression alterations of ethylparaben occurred at concentrations of approximately three orders of magnitude above human blood concentrations; moreover, the substance fitted well into a scenario where toxicity increases with the alkyl chain length, and gene expression changes only occur at cytotoxic or close to cytotoxic concentrations. Therefore, no evidence was obtained suggesting that in vivo DART with ethylparaben would lead to different results as the methyl- or propyl derivates.

Performance assessment and economic analysis of a human Liver-Chip for predictive toxicology

BACKGROUND

Conventional preclinical models often miss drug toxicities, meaning the harm these drugs pose to humans is only realized in clinical trials or when they make it to market. This has caused the pharmaceutical industry to waste considerable time and resources developing drugs destined to fail. Organ-on-a-Chip technology has the potential improve success in drug development pipelines, as it can recapitulate organ-level pathophysiology and clinical responses; however, systematic and quantitative evaluations of Organ-Chips' predictive value have not yet been reported.

METHODS

870 Liver-Chips were analyzed to determine their ability to predict drug-induced liver injury caused by small molecules identified as benchmarks by the Innovation and Quality consortium, who has published guidelines defining criteria for qualifying preclinical models. An economic analysis was also performed to measure the value Liver-Chips could offer if they were broadly adopted in supporting toxicity-related decisions as part of preclinical development workflows.

RESULTS

Here, we show that the Liver-Chip met the qualification guidelines across a blinded set of 27 known hepatotoxic and non-toxic drugs with a sensitivity of 87% and a specificity of 100%. We also show that this level of performance could generate over $3 billion annually for the pharmaceutical industry through increased small-molecule R&D productivity.

CONCLUSIONS

The results of this study show how incorporating predictive Organ-Chips into drug development workflows could substantially improve drug discovery and development, allowing manufacturers to bring safer, more effective medicines to market in less time and at lower costs.

Neuronal differentiation pathways and compound-induced developmental neurotoxicity in the human neural progenitor cell test (hNPT) revealed by RNA-seq

There is an increased awareness that the use of animals for compound-induced developmental neurotoxicity (DNT) testing has limitations. Animal-free innovations, especially the ones based on human stem cell-based models are pivotal in studying DNT since they can mimic processes relevant to human brain development. Here we present the human neural progenitor test (hNPT), a 10-day protocol in which neural progenitor cells differentiate into a neuron-astrocyte co-culture. The study aimed to characterise differentiation over time and to find neurodevelopmental processes sensitive to compound exposure using transcriptomics. 3992 genes regulated in unexposed control cultures (p ≤ 0.001, log2FC ≥ 1) showed Gene Ontology (GO-) term enrichment for neuronal and glial differentiation, neurite extension, synaptogenesis, and synaptic transmission. Exposure to known or suspected DNT compounds (acrylamide, chlorpyrifos, fluoxetine, methyl mercury, or valproic acid) at concentrations resulting in 95% cell viability each regulated unique combinations of GO-terms relating to neural progenitor proliferation, neuronal and glial differentiation, axon development, synaptogenesis, synaptic transmission, and apoptosis. Investigation of the GO-terms 'neuron apoptotic process' and 'axon development' revealed common genes that were responsive across compounds, and might be used as biomarkers for DNT. The GO-term 'synaptic signalling', on the contrary, whilst also responsive to all compounds tested, showed little overlap in gene expression regulation patterns between the conditions. This GO-term may articulate compound-specific effects that may be relevant for revealing differences in mechanism of toxicity. Given its focus on neural progenitor cell to mature multilineage neuronal cell maturation and its detailed molecular readout based on gene expression analysis, hNPT might have added value as a tool for neurodevelopmental toxicity testing in vitro. Further assessment of DNT-specific biomarkers that represent these processes needs further studies.

Optimization of the TeraTox assay for preclinical teratogenicity assessment

Current animal-free methods to assess teratogenicity of drugs under development still deliver high numbers of false negatives. To improve the sensitivity of human teratogenicity prediction, we characterized the TeraTox test, a newly developed multilineage differentiation assay using 3D human-induced pluripotent stem cells. TeraTox produces primary output concentration-dependent cytotoxicity and altered gene expression induced by each test compound. These data are fed into an interpretable machine-learning model to perform prediction, which relates to the concentration-dependent human teratogenicity potential of drug candidates. We applied TeraTox to profile 33 approved pharmaceuticals and 12 proprietary drug candidates with known in vivo data. Comparing TeraTox predictions with known human or animal toxicity, we report an accuracy of 69% (specificity: 53%, sensitivity: 79%). TeraTox performed better than 2 quantitative structure-activity relationship models and had a higher sensitivity than the murine embryonic stem cell test (accuracy: 58%, specificity: 76%, and sensitivity: 46%) run in the same laboratory. The overall prediction accuracy could be further improved by combining TeraTox and mouse embryonic stem cell test results. Furthermore, patterns of altered gene expression revealed by TeraTox may help grouping toxicologically similar compounds and possibly deducing common modes of action. The TeraTox assay and the dataset described here therefore represent a new tool and a valuable resource for drug teratogenicity assessment.

Classification of Developmental Toxicants in a Human iPSC Transcriptomics-Based Test

Despite the progress made in developmental toxicology, there is a great need for in vitro tests that identify developmental toxicants in relation to human oral doses and blood concentrations. In the present study, we established the hiPSC-based UKK2 in vitro test and analyzed genome-wide expression profiles of 23 known teratogens and 16 non-teratogens. Compounds were analyzed at the maximal plasma concentration (C_max) and at 20-fold C_max for a 24 h incubation period in three independent experiments. Based on the 1000 probe sets with the highest variance and including information on cytotoxicity, penalized logistic regression with leave-one-out cross-validation was used to classify the compounds as test-positive or test-negative, reaching an area under the curve (AUC), accuracy, sensitivity, and specificity of 0.96, 0.92, 0.96, and 0.88, respectively. Omitting the cytotoxicity information reduced the test performance to an AUC of 0.94, an accuracy of 0.79, and a sensitivity of 0.74. A second method, which used the number of significantly deregulated probe sets to classify the compounds, resulted in a specificity of 1; however, the AUC (0.90), accuracy (0.90), and sensitivity (0.83) were inferior compared to those of the logistic regression-based procedure. Finally, no increased performance was achieved when the high test concentrations (20-fold C_max) were used, in comparison to testing within the realistic clinical range (1-fold C_max). In conclusion, although further optimization is required, for example, by including additional readouts and cell systems that model different developmental processes, the UKK2-test in its present form can support the early discovery-phase detection of human developmental toxicants.

A novel human stem cell-based biomarker assay for in vitro assessment of developmental toxicity

BACKGROUND

Testing for developmental toxicity according to the current regulatory guidelines requires large numbers of animals, making these tests very resource intensive, time-consuming, and ethically debatable. Over the past decades, several alternative in vitro assays have been developed, but these often suffered from low predictability and the inability to provide a mechanistic understanding of developmental toxicity.

METHODS

To identify embryotoxic compounds, we developed a human induced pluripotent stem cells (hiPSCs)-based biomarker assay. The assay is based on the differentiation of hiPSCs into functional cardiomyocytes and hepatocytes. Proper stem cell differentiation is investigated by morphological profiling and assessment of time-dependent expression patterns of cell-specific biomarkers. In this system, a decrease in the expression of the biomarker genes and morphology disruption of the differentiated cells following compound treatment indicated teratogenicity.

RESULTS

The hiPSCs-based biomarker assay was validated with 21 well-established in vivo animal teratogenic and non-teratogenic compounds during cardiomyocyte and hepatocyte differentiation. The in vivo teratogenic compounds (e.g., thalidomide and valproic acid) markedly disrupted morphology, functionality, and the expression pattern of the biomarker genes in either one or both cell types. Non-teratogenic chemicals generally had no effect on the morphology of differentiated cells, nor on the expression of the biomarker genes. Compared to the in vivo classification, the assay achieved high accuracy (91%), sensitivity (91%), and specificity (90%).

CONCLUSION

The assay, which we named ReproTracker®, is a state-of-the-art in vitro method that can identify the teratogenicity potential of new pharmaceuticals and chemicals and signify the outcome of in vivo test systems.

OUP accepted manuscript

In vitro models of the peripheral nervous system would benefit from further refinements to better support studies on neuropathies. In particular, the assessment of pain-related signals is still difficult in human cell cultures. Here, we harnessed induced pluripotent stem cells (iPSCs) to generate peripheral sensory neurons enriched in nociceptors. The objective was to generate a culture system with signaling endpoints suitable for pharmacological and toxicological studies. Neurons generated by conventional differentiation protocols expressed moderate levels of P2X3 purinergic receptors and only low levels of TRPV1 capsaicin receptors, when maturation time was kept to the upper practically useful limit of 6 weeks. As alternative approach, we generated cells with an inducible NGN1 transgene. Ectopic expression of this transcription factor during a defined time window of differentiation resulted in highly enriched nociceptor cultures, as determined by functional (P2X3 and TRPV1 receptors) and immunocytochemical phenotyping, complemented by extensive transcriptome profiling. Single cell recordings of Ca²⁺-indicator fluorescence from >9000 cells were used to establish the “fraction of reactive cells” in a stimulated population as experimental endpoint, that appeared robust, transparent and quantifiable. To provide an example of application to biomedical studies, functional consequences of prolonged exposure to the chemotherapeutic drug oxaliplatin were examined at non-cytotoxic concentrations. We found (i) neuronal (allodynia-like) hypersensitivity to otherwise non-activating mechanical stimulation that could be blocked by modulators of voltage-gated sodium channels; (ii) hyper-responsiveness to TRPV1 receptor stimulation. These findings and several other measured functional alterations indicate that the model is suitable for pharmacological and toxicological studies related to peripheral neuropathies.

RNA modifications as emerging therapeutic targets

The field of epitranscriptome, posttranscriptional modifications to RNAs, is still growing up and has presented substantial evidences for the role of RNA modifications in diseases. In terms of new drug development, RNA modifications have a great promise for therapy. For example, more than 170 type of modifications exist in various types of RNAs. Regulatory genes and their roles in critical biological process have been identified and they are associated with several diseases. Current data, for example, identification of inhibitors targeting RNA modifications regulatory genes, strongly support the idea that RNA modifications have potential as emerging therapeutic targets. Therefore, in this review, RNA modifications and regulatory genes were comprehensively documented in terms of drug development by summarizing the findings from previous studies. It was discussed how RNA modifications or regulatory genes can be targeted by altering molecular mechanisms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > RNA Editing and Modification.

Mapping the cellular response to electron transport chain inhibitors reveals selective signaling networks triggered by mitochondrial perturbation

Mitochondrial perturbation is a key event in chemical-induced organ toxicities that is incompletely understood. Here, we studied how electron transport chain (ETC) complex I, II, or III (CI, CII and CIII) inhibitors affect mitochondrial functionality, stress response activation, and cell viability using a combination of high-content imaging and TempO-Seq in HepG2 hepatocyte cells. CI and CIII inhibitors perturbed mitochondrial membrane potential (MMP) and mitochondrial and cellular ATP levels in a concentration- and time-dependent fashion and, under conditions preventing a switch to glycolysis attenuated cell viability, whereas CII inhibitors had no effect. TempO-Seq analysis of changes in mRNA expression pointed to a shared cellular response to CI and CIII inhibition. First, to define specific ETC inhibition responses, a gene set responsive toward ETC inhibition (and not to genotoxic, oxidative, or endoplasmic reticulum stress) was identified using targeted TempO-Seq in HepG2. Silencing of one of these genes, NOS3, exacerbated the impact of CI and CIII inhibitors on cell viability, indicating its functional implication in cellular responses to mitochondrial stress. Then by monitoring dynamic responses to ETC inhibition using a HepG2 GFP reporter panel for different classes of stress response pathways and applying pathway and gene network analysis to TempO-Seq data, we looked for downstream cellular events of ETC inhibition and identified the amino acid response (AAR) as being triggered in HepG2 by ETC inhibition. Through in silico approaches we provide evidence indicating that a similar AAR is associated with exposure to mitochondrial toxicants in primary human hepatocytes. Altogether, we (i) unravel quantitative, time- and concentration-resolved cellular responses to mitochondrial perturbation, (ii) identify a gene set associated with adaptation to exposure to active ETC inhibitors, and (iii) show that ER stress and an AAR accompany ETC inhibition in HepG2 and primary hepatocytes.

Evaluation of postmortem microarray data in bipolar disorder using traditional data comparison and artificial intelligence reveals novel gene targets

Large-scale microarray studies on post-mortem brain tissues have been utilized to investigate the complex molecular pathology of bipolar disorder. However, a major challenge in characterizing the dysregulation of gene expression in patients with bipolar disorder includes the lack of convergence between different studies, limiting comprehensive understanding from individual results. In this study, we aimed to identify genes that are both validated in published literature and are important classification features of unsupervised machine learning analysis of Stanley Brain Bank microarray database, followed by augmented intelligence method to identify distinct patient molecular subgroups. Through combining traditional literature approaches and machine learning, we identified TBL1XR1, SMARCA2, and CHMP5 to be replicated in 3 of the 4 studies included our analysis. The expression of these genes segregated unique subgroups of patients with bipolar disorder. Our study suggests the involvement of PPARγ pathway regulation in patients with bipolar disorder.

Reproducibility and robustness of high-throughput S1500+ transcriptomics on primary rat hepatocytes for chemical-induced hepatotoxicity assessment

•

TempO-Seq assays of rat hepatocytes in collagen sandwich are highly reproducible.

•

Gene expression analysis shows S1500+ is representative of the whole transcriptome.

•

Connectivity mapping shows consistency between TempO-Seq and Affymetrix data.

•

Gene set enrichment shows consistency between S1500+ and the whole transcriptome.

•

Gene set enrichment using hallmark gene sets informs hepatotoxicity.

Cell-based in vitro models coupled with high-throughput transcriptomics (HTTr) are increasingly utilized as alternative methods to animal-based toxicity testing. Here, using a panel of 14 chemicals with different risks of human drug-induced liver injury (DILI) and two dosing concentrations, we evaluated an HTTr platform comprised of collagen sandwich primary rat hepatocyte culture and the TempO-Seq surrogate S1500+ (ST) assay. First, the HTTr platform was found to exhibit high reproducibility between technical and biological replicates (r greater than 0.85). Connectivity mapping analysis further demonstrated a high level of inter-platform reproducibility between TempO-Seq data and Affymetrix GeneChip data from the Open TG-GATES project. Second, the TempO-Seq ST assay was shown to be a robust surrogate to the whole transcriptome (WT) assay in capturing chemical-induced changes in gene expression, as evident from correlation analysis, PCA and unsupervised hierarchical clustering. Gene set enrichment analysis (GSEA) using the Hallmark gene set collection also demonstrated consistency in enrichment scores between ST and WT assays. Lastly, unsupervised hierarchical clustering of hallmark enrichment scores broadly divided the samples into hepatotoxic, intermediate, and non-hepatotoxic groups. Xenobiotic metabolism, bile acid metabolism, apoptosis, p53 pathway, and coagulation were found to be the key hallmarks driving the clustering. Taken together, our results established the reproducibility and performance of collagen sandwich culture in combination with TempO-Seq S1500+ assay, and demonstrated the utility of GSEA using the hallmark gene set collection to identify potential hepatotoxicants for further validation.

A human stem cell-derived test system for agents modifying neuronal N-methyl-D-aspartate-type glutamate receptor Ca2+-signalling

Methods to assess neuronal receptor functions are needed in toxicology and for drug development. Human-based test systems that allow studies on glutamate signalling are still scarce. To address this issue, we developed and characterized pluripotent stem cell (PSC)-based neural cultures capable of forming a functional network. Starting from a stably proliferating neuroepithelial stem cell (NESC) population, we generate “mixed cortical cultures” (MCC) within 24 days. Characterization by immunocytochemistry, gene expression profiling and functional tests (multi-electrode arrays) showed that MCC contain various functional neurotransmitter receptors, and in particular, the N-methyl-d-aspartate subtype of ionotropic glutamate receptors (NMDA-R). As this important receptor is found neither on conventional neural cell lines nor on most stem cell-derived neurons, we focused here on the characterization of rapid glutamate-triggered Ca²⁺ signalling. Changes of the intracellular free calcium ion concentration ([Ca²⁺]_i) were measured by fluorescent imaging as the main endpoint, and a method to evaluate and quantify signals in hundreds of cells at the same time was developed. We observed responses to glutamate in the low µM range. MCC responded to kainate and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA), and a subpopulation of 50% had functional NMDA-R. The receptor was modulated by Mg²⁺, Zn²⁺ and Pb²⁺ in the expected ways, and various toxicologically relevant agonists (quinolinic acid, ibotenic acid, domoic acid) triggered [Ca²⁺]_i responses in MCC. Antagonists, such as phencyclidine, ketamine and dextromethorphan, were also readily identified. Thus, the MCC developed here may fill an important gap in the panel of test systems available to characterize the effects of chemicals on neurotransmitter receptors.

Kinetic modeling of stem cell transcriptome dynamics to identify regulatory modules of normal and disturbed neuroectodermal differentiation

Thousands of transcriptome data sets are available, but approaches for their use in dynamic cell response modelling are few, especially for processes affected simultaneously by two orthogonal influencing variables. We approached this problem for neuroepithelial development of human pluripotent stem cells (differentiation variable), in the presence or absence of valproic acid (signaling variable). Using few basic assumptions (sequential differentiation states of cells; discrete on/off states for individual genes in these states), and time-resolved transcriptome data, a comprehensive model of spontaneous and perturbed gene expression dynamics was developed. The model made reliable predictions (average correlation of 0.85 between predicted and subsequently tested expression values). Even regulations predicted to be non-monotonic were successfully validated by PCR in new sets of experiments. Transient patterns of gene regulation were identified from model predictions. They pointed towards activation of Wnt signaling as a candidate pathway leading to a redirection of differentiation away from neuroepithelial cells towards neural crest. Intervention experiments, using a Wnt/beta-catenin antagonist, led to a phenotypic rescue of this disturbed differentiation. Thus, our broadly applicable model allows the analysis of transcriptome changes in complex time/perturbation matrices.

Development of a neural rosette formation assay (RoFA) to identify neurodevelopmental toxicants and to characterize their transcriptome disturbances

The first in vitro tests for developmental toxicity made use of rodent cells. Newer teratology tests, e.g. developed during the ESNATS project, use human cells and measure mechanistic endpoints (such as transcriptome changes). However, the toxicological implications of mechanistic parameters are hard to judge, without functional/morphological endpoints. To address this issue, we developed a new version of the human stem cell-based test STOP-tox_(UKN). For this purpose, the capacity of the cells to self-organize to neural rosettes was assessed as functional endpoint: pluripotent stem cells were allowed to differentiate into neuroepithelial cells for 6 days in the presence or absence of toxicants. Then, both transcriptome changes were measured (standard STOP-tox_(UKN)) and cells were allowed to form rosettes. After optimization of staining methods, an imaging algorithm for rosette quantification was implemented and used for an automated rosette formation assay (RoFA). Neural tube toxicants (like valproic acid), which are known to disturb human development at stages when rosette-forming cells are present, were used as positive controls. Established toxicants led to distinctly different tissue organization and differentiation stages. RoFA outcome and transcript changes largely correlated concerning (1) the concentration-dependence, (2) the time dependence, and (3) the set of positive hits identified amongst 24 potential toxicants. Using such comparative data, a prediction model for the RoFA was developed. The comparative analysis was also used to identify gene dysregulations that are particularly predictive for disturbed rosette formation. This 'RoFA predictor gene set' may be used for a simplified and less costly setup of the STOP-tox_(UKN) assay.

Current Availability of Stem Cell-Based In Vitro Methods for Developmental Neurotoxicity (DNT) Testing

There is evidence that chemical exposure during development can cause irreversible impairments of the human developing nervous system. Therefore, testing compounds for their developmentally neurotoxic potential has high priority for different stakeholders: academia, industry, and regulatory bodies. Due to the resource-intensity of current developmental neurotoxicity (DNT) in vivo guidelines, alternative methods that are scientifically valid and have a high predictivity for humans are especially desired by regulators. Here, we review availability of stem-/progenitor cell-based in vitro methods for DNT evaluation that is based on the concept of neurodevelopmental process assessment. These test methods are assembled into a DNT in vitro testing battery. Gaps in this testing battery addressing research needs are also pointed out.

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 (C_max) 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 (EC₁₀) yielded better metrics than higher toxicity thresholds (EC₅₀); (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 C_max were best suited to estimate human toxicity. With a support vector machine-based classifier, using EC₁₀ and C_max 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.

Road Map for Development of Stem Cell-Based Alternative Test Methods

Much progress has been made in establishing strategies for differentiation of induced human pluripotent stem cells (hiPSCs). However, differentiated hiPSCs are not yet routinely used for prediction of toxicity. Here, limiting factors are summarised and possibilities for improvement are discussed, with a focus on hepatocytes, cardiomyocytes, tubular epithelial cells, and developmental toxicity. Moreover, we make recommendations for further fine-tuning of differentiation protocols for hiPSCs to hepatocyte-like cells by comparing individual steps of currently available protocols to the mechanisms occurring during embryonic development. A road map is proposed to facilitate test system development, including a description of the most useful performance metrics.

Assessment of stem cell differentiation based on genome-wide expression profiles

In recent years, protocols have been established to differentiate stem and precursor cells into more mature cell types. However, progress in this field has been hampered by difficulties to assess the differentiation status of stem cell-derived cells in an unbiased manner. Here, we present an analysis pipeline based on published data and methods to quantify the degree of differentiation and to identify transcriptional control factors explaining differences from the intended target cells or tissues. The pipeline requires RNA-Seq or gene array data of the stem cell starting population, derived 'mature' cells and primary target cells or tissue. It consists of a principal component analysis to represent global expression changes and to identify possible problems of the dataset that require special attention, such as: batch effects; clustering techniques to identify gene groups with similar features; over-representation analysis to characterize biological motifs and transcriptional control factors of the identified gene clusters; and metagenes as well as gene regulatory networks for quantitative cell-type assessment and identification of influential transcription factors. Possibilities and limitations of the analysis pipeline are illustrated using the example of human embryonic stem cell and human induced pluripotent cells to generate 'hepatocyte-like cells'. The pipeline quantifies the degree of incomplete differentiation as well as remaining stemness and identifies unwanted features, such as colon- and fibroblast-associated gene clusters that are absent in real hepatocytes but typically induced by currently available differentiation protocols. Finally, transcription factors responsible for incomplete and unwanted differentiation are identified. The proposed method is widely applicable and allows an unbiased and quantitative assessment of stem cell-derived cells.This article is part of the theme issue 'Designer human tissue: coming to a lab near you'.

Comparison of base-line and chemical-induced transcriptomic responses in HepaRG and RPTEC/TERT1 cells using TempO-Seq

The utilisation of genome-wide transcriptomics has played a pivotal role in advancing the field of toxicology, allowing the mapping of transcriptional signatures to chemical exposures. These activities have uncovered several transcriptionally regulated pathways that can be utilised for assessing the perturbation impact of a chemical and also the identification of toxic mode of action. However, current transcriptomic platforms are not very amenable to high-throughput workflows due to, high cost, complexities in sample preparation and relatively complex bioinformatic analysis. Thus, transcriptomic investigations are usually limited in dose and time dimensions and are, therefore, not optimal for implementation in risk assessment workflows. In this study, we investigated a new cost-effective, transcriptomic assay, TempO-Seq, which alleviates the aforementioned limitations. This technique was evaluated in a 6-compound screen, utilising differentiated kidney (RPTEC/TERT1) and liver (HepaRG) cells and compared to non-transcriptomic label-free sensitive endpoints of chemical-induced disturbances, namely phase contrast morphology, xCELLigence and glycolysis. Non-proliferating cell monolayers were exposed to six sub-lethal concentrations of each compound for 24 h. The results show that utilising a 2839 gene panel, it is possible to discriminate basal tissue-specific signatures, generate dose-response relationships and to discriminate compound-specific and cell type-specific responses. This study also reiterates previous findings that chemical-induced transcriptomic alterations occur prior to cytotoxicity and that transcriptomics provides in depth mechanistic information of the effects of chemicals on cellular transcriptional responses. TempO-Seq is a robust transcriptomic platform that is well suited for in vitro toxicity experiments.

A transcriptome comparison of time-matched developing human, mouse and rat neural progenitor cells reveals human uniqueness

It is widely accepted that human brain development has unique features that cannot be represented by rodents. Obvious reasons are the evolutionary distance and divergent physiology. This might lead to false predictions when rodents are used for safety or pharmacological efficacy studies. For a better translation of animal-based research to the human situation, human in vitro systems might be useful. In this study, we characterize developing neural progenitor cells from prenatal human and time-matched rat and mouse brains by analyzing the changes in their transcriptome profile during neural differentiation. Moreover, we identify hub molecules that regulate neurodevelopmental processes like migration and differentiation. Consequences of modulation of three of those hubs on these processes were studied in a species-specific context. We found that although the gene expression profiles of the three species largely differ qualitatively and quantitatively, they cluster in similar GO terms like cell migration, gliogenesis, neurogenesis or development of multicellular organism. Pharmacological modulation of the identified hub molecules triggered species-specific cellular responses. This study underlines the importance of understanding species differences on the molecular level and advocates the use of human based in vitro models for pharmacological and toxicological research.

Recommendation on test readiness criteria for new approach methods in toxicology: Exemplified for developmental neurotoxicity

Multiple non-animal-based test methods have never been formally validated. In order to use such new approach methods (NAMs) in a regulatory context, criteria to define their readiness are necessary. The field of developmental neurotoxicity (DNT) testing is used to exemplify the application of readiness criteria. The costs and number of untested chemicals are overwhelming for in vivo DNT testing. Thus, there is a need for inexpensive, high-throughput NAMs, to obtain initial information on potential hazards, and to allow prioritization for further testing. A background on the regulatory and scientific status of DNT testing is provided showing different types of test readiness levels, depending on the intended use of data from NAMs. Readiness criteria, compiled during a stakeholder workshop, uniting scientists from academia, industry and regulatory authorities are presented. An important step beyond the listing of criteria, was the suggestion for a preliminary scoring scheme. On this basis a (semi)-quantitative analysis process was assembled on test readiness of 17 NAMs with respect to various uses (e.g. prioritization/screening, risk assessment). The scoring results suggest that several assays are currently at high readiness levels. Therefore, suggestions are made on how DNT NAMs may be assembled into an integrated approach to testing and assessment (IATA). In parallel, the testing state in these assays was compiled for more than 1000 compounds. Finally, a vision is presented on how further NAM development may be guided by knowledge of signaling pathways necessary for brain development, DNT pathophysiology, and relevant adverse outcome pathways (AOP).

New animal-free concepts and test methods for developmental toxicity and peripheral neurotoxicity

The complex toxicological fields of repeat dose organ toxicity (RDT) and developmental and reproductive toxicity (DART) still require new concepts and approaches to achieve a fully animal-free safety assessment of chemicals. One novel approach is the generation of relevant human cell types from pluripotent stem cells, and the use of such cells for the establishment of phenotypic test methods. Due to their broad endpoints, such tests capture multiple types of toxicants, i.e. they are a readout for the activation of many adverse outcome pathways (AOPs). The 2016 Lush Science Prize was awarded for the development of one such assay, the PeriTox test, which uses human peripheral neurons generated from stem cells. The assay endpoints measure various cell functions, and these give information on the potential neurotoxicity and developmental neurotoxicity hazard of test compounds. The PeriTox test method has a high predictivity and sensitivity for peripheral neurotoxicants, and thus addresses the inherent challenges in pesticide testing and drug development. Data from the test can be obtained quickly and at a relatively high-throughput, and thus, the assay has the potential to replace animal-based safety assessment during early product development or for screening potential environmental toxicants.

Whole genome microarray analysis of neural progenitor C17.2 cells during differentiation and validation of 30 neural mRNA biomarkers for estimation of developmental neurotoxicity

Despite its high relevance, developmental neurotoxicity (DNT) is one of the least studied forms of toxicity. Current guidelines for DNT testing are based on in vivo testing and they require extensive resources. Transcriptomic approaches using relevant in vitro models have been suggested as a useful tool for identifying possible DNT-generating compounds. In this study, we performed whole genome microarray analysis on the murine progenitor cell line C17.2 following 5 and 10 days of differentiation. We identified 30 genes that are strongly associated with neural differentiation. The C17.2 cell line can be differentiated into a co-culture of both neurons and neuroglial cells, giving a more relevant picture of the brain than using neuronal cells alone. Among the most highly upregulated genes were genes involved in neurogenesis (CHRDL1), axonal guidance (BMP4), neuronal connectivity (PLXDC2), axonogenesis (RTN4R) and astrocyte differentiation (S100B). The 30 biomarkers were further validated by exposure to non-cytotoxic concentrations of two DNT-inducing compounds (valproic acid and methylmercury) and one neurotoxic chemical possessing a possible DNT activity (acrylamide). Twenty-eight of the 30 biomarkers were altered by at least one of the neurotoxic substances, proving the importance of these biomarkers during differentiation. These results suggest that gene expression profiling using a predefined set of biomarkers could be used as a sensitive tool for initial DNT screening of chemicals. Using a predefined set of mRNA biomarkers, instead of the whole genome, makes this model affordable and high-throughput. The use of such models could help speed up the initial screening of substances, possibly indicating alerts that need to be further studied in more sophisticated models.

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 µ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 cardiomyocyte-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.

Can Valproic Acid Regulate Neurogenesis from Nestin+ Cells in the Adult Midbrain?

Degeneration of dopamine (DA) neurons in the substantia nigra pars compacta (SNc) causes the motor symptoms (e.g. tremor, muscle rigidity, bradykinesia, postural instability) of Parkinson's disease (PD). It is generally agreed that replacing these neurons will provide better motor symptom relief and fewer side effects than current pharmacotherapies. One potential approach to this is up-regulating endogenous DA neurogenesis in SNc. In the present study, we conducted bioinformatics analyses to identify signalling pathways that control expression of Pax6 and Msx1 genes, which have been identified as potentially important neurogenic regulators in the adult midbrain. From this Valproic acid (VPA) was identified as a regulator of these pathways, and we tested VPA for its ability to regulate midbrain neurogenesis in adult mice. VPA was infused directly into the midbrain of adult NesCreER^T2/R26eYFP mice using osmotic pumps attached to implanted cannula. These mice enable permanent eYFP+ labelling of adult Nestin-expressing neural precursor cells and their progeny/ontogeny. VPA did not affect the number of eYFP+ midbrain cells, but significantly reduced the number of Pax6+, Pax6+/NeuN+, eYFP+/NeuN+ and eYFP-/NeuN+ cells. However, this reduction in NeuN expression was probably via VPA's Histone de-acetylase inhibitory properties rather than reduced neuronal differentiation by eYFP + cells. We conclude that Pax6 and Msx1 are not viable targets for regulating neurogenesis in the adult midbrain.

Impairment of human neural crest cell migration by prolonged exposure to interferon-beta

Human cell-based toxicological assays have been used successfully to detect known toxicants, and to distinguish them from negative controls. However, there is at present little experience on how to deal with hits from screens of compounds with yet unknown hazard. As a case study to this issue, we characterized human interferon-beta (IFNβ) as potential developmental toxicant affecting neural crest cells (NCC). The protein was identified as a hit during a screen of clinically used drugs in the ‘migration inhibition of neural crest’ (MINC) assay. Concentration–response studies in the MINC combined with immunocytochemistry and mRNA quantification of cellular markers showed that IFNβ inhibited NCC migration at concentrations as low as 20 pM. The effective concentrations found here correspond to levels found in human plasma, and they were neither cytostatic nor cytotoxic nor did they did they affect the differentiation state and overall phenotype of NCC. Data from two other migration assays confirmed that picomolar concentration of IFNβ reduced the motility of NCC, while other interferons were less potent. The activation of JAK kinase by IFNβ, as suggested by bioinformatics analysis of the transcriptome changes, was confirmed by biochemical methods. The degree and duration of pathway activation correlated with the extent of migration inhibition, and pharmacological block of this signaling pathway before, or up to 6 h after exposure to the cytokine prevented the effects of IFNβ on migration. Thus, the reduction of vital functions of human NCC is a hitherto unknown potential hazard of endogenous or pharmacologically applied interferons.

Stem Cell Transcriptome Responses and Corresponding Biomarkers That Indicate the Transition from Adaptive Responses to Cytotoxicity

Analysis of transcriptome changes has become an established method to characterize the reaction of cells to toxicants. Such experiments are mostly performed at compound concentrations close to the cytotoxicity threshold. At present, little information is available on concentration-dependent features of transcriptome changes, in particular, at the transition from noncytotoxic concentrations to conditions that are associated with cell death. Thus, it is unclear in how far cell death confounds the results of transcriptome studies. To explore this gap of knowledge, we treated pluripotent stem cells differentiating to human neuroepithelial cells (UKN1 assay) for short periods (48 h) with increasing concentrations of valproic acid (VPA) and methyl mercury (MeHg), two compounds with vastly different modes of action. We developed various visualization tools to describe cellular responses, and the overall response was classified as "tolerance" (minor transcriptome changes), "functional adaptation" (moderate/strong transcriptome responses, but no cytotoxicity), and "degeneration". The latter two conditions were compared, using various statistical approaches. We identified (i) genes regulated at cytotoxic, but not at noncytotoxic, concentrations and (ii) KEGG pathways, gene ontology term groups, and superordinate biological processes that were only regulated at cytotoxic concentrations. The consensus markers and processes found after 48 h treatment were then overlaid with those found after prolonged (6 days) treatment. The study highlights the importance of careful concentration selection and of controlling viability for transcriptome studies. Moreover, it allowed identification of 39 candidate "biomarkers of cytotoxicity". These could serve to provide alerts that data sets of interest may have been affected by cell death in the model system studied.

Stem cell toxicology: a powerful tool to assess pollution effects on human health

Environmental pollution is a global problem; the lack of comprehensive toxicological assessments may lead to increased health risks. To fully understand the health effects of pollution, it is paramount to implement fast, efficient and specific toxicity screening that relies on human models rather than on time-consuming, expensive and often inaccurate tests involving live animals. Human stem cell toxicology represents a valid alternative to traditional toxicity assays because it takes advantage of the ability of stem cells to differentiate into multiple cell types and tissues of the human body. Thus, this branch of toxicology provides a possibility to assess cellular, embryonic, developmental, reproductive and functional toxicity in vitro within a single system highly relevant to human physiology. In this review, we describe the development, performance and future perspectives of stem cell toxicology, with an emphasis on how it can meet the increasing challenges posed by environmental pollution in the modern world.

In vitro acute and developmental neurotoxicity screening: an overview of cellular platforms and high-throughput technical possibilities

Neurotoxicity and developmental neurotoxicity are important issues of chemical hazard assessment. Since the interpretation of animal data and their extrapolation to man is challenging, and the amount of substances with information gaps exceeds present animal testing capacities, there is a big demand for in vitro tests to provide initial information and to prioritize for further evaluation. During the last decade, many in vitro tests emerged. These are based on animal cells, human tumour cell lines, primary cells, immortalized cell lines, embryonic stem cells, or induced pluripotent stem cells. They differ in their read-outs and range from simple viability assays to complex functional endpoints such as neural crest cell migration. Monitoring of toxicological effects on differentiation often requires multiomics approaches, while the acute disturbance of neuronal functions may be analysed by assessing electrophysiological features. Extrapolation from in vitro data to humans requires a deep understanding of the test system biology, of the endpoints used, and of the applicability domains of the tests. Moreover, it is important that these be combined in the right way to assess toxicity. Therefore, knowledge on the advantages and disadvantages of all cellular platforms, endpoints, and analytical methods is essential when establishing in vitro test systems for different aspects of neurotoxicity. The elements of a test, and their evaluation, are discussed here in the context of comprehensive prediction of potential hazardous effects of a compound. We summarize the main cellular characteristics underlying neurotoxicity, present an overview of cellular platforms and read-out combinations assessing distinct parts of acute and developmental neurotoxicology, and highlight especially the use of stem cell-based test systems to close gaps in the available battery of tests.