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 (Cmax) and at 20-fold Cmax 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 Cmax) were used, in comparison to testing within the realistic clinical range (1-fold Cmax). 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.

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.

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 Ca2+ signalling. Changes of the intracellular free calcium ion concentration ([Ca2+]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 Mg2+, Zn2+ and Pb2+ in the expected ways, and various toxicologically relevant agonists (quinolinic acid, ibotenic acid, domoic acid) triggered [Ca2+]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.

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.

The EU-ToxRisk method documentation, data processing and chemical testing pipeline for the regulatory use of new approach methods.

Hazard assessment, based on new approach methods (NAM), requires the use of batteries of assays, where individual tests may be contributed by different laboratories. A unified strategy for such collaborative testing is presented. It details all procedures required to allow test information to be usable for integrated hazard assessment, strategic project decisions and/or for regulatory purposes. The EU-ToxRisk project developed a strategy to provide regulatorily valid data, and exemplified this using a panel of > 20 assays (with > 50 individual endpoints), each exposed to 19 well-known test compounds (e.g. rotenone, colchicine, mercury, paracetamol, rifampicine, paraquat, taxol). Examples of strategy implementation are provided for all aspects required to ensure data validity: (i) documentation of test methods in a publicly accessible database; (ii) deposition of standard operating procedures (SOP) at the European Union DB-ALM repository; (iii) test readiness scoring accoding to defined criteria; (iv) disclosure of the pipeline for data processing; (v) link of uncertainty measures and metadata to the data; (vi) definition of test chemicals, their handling and their behavior in test media; (vii) specification of the test purpose and overall evaluation plans. Moreover, data generation was exemplified by providing results from 25 reporter assays. A complete evaluation of the entire test battery will be described elsewhere. A major learning from the retrospective analysis of this large testing project was the need for thorough definitions of the above strategy aspects, ideally in form of a study pre-registration, to allow adequate interpretation of the data and to ensure overall scientific/toxicological validity.

The DEK oncoprotein is a Su(var) that is essential to heterochromatin integrity.

Heterochromatin integrity is crucial for genome stability and regulation of gene expression, but the factors involved in mammalian heterochromatin biology are only incompletely understood. Here we identify the oncoprotein DEK, an abundant nuclear protein with a previously enigmatic in vivo function, as a Suppressor of Variegation [Su(var)] that is crucial to global heterochromatin integrity. We show that DEK interacts directly with Heterochromatin Protein 1 α (HP1α) and markedly enhances its binding to trimethylated H3K9 (H3K9me3), which is key for maintaining heterochromatic regions. Loss of Dek in Drosophila leads to a Su(var) phenotype and global reduction in heterochromatin. Thus, these findings show that DEK is a key factor in maintaining the balance between heterochromatin and euchromatin in vivo.

Stage-specific metabolic features of differentiating neurons: Implications for toxicant sensitivity.

Developmental neurotoxicity (DNT) may be induced when chemicals disturb a key neurodevelopmental process, and many tests focus on this type of toxicity. Alternatively, DNT may occur when chemicals are cytotoxic only during a specific neurodevelopmental stage. The toxicant sensitivity is affected by the expression of toxicant targets and by resilience factors. Although cellular metabolism plays an important role, little is known how it changes during human neurogenesis, and how potential alterations affect toxicant sensitivity of mature vs. immature neurons. We used immature (d0) and mature (d6) LUHMES cells (dopaminergic human neurons) to provide initial answers to these questions. Transcriptome profiling and characterization of energy metabolism suggested a switch from predominantly glycolytic energy generation to a more pronounced contribution of the tricarboxylic acid cycle (TCA) during neuronal maturation. Therefore, we used pulsed stable isotope-resolved metabolomics (pSIRM) to determine intracellular metabolite pool sizes (concentrations), and isotopically non-stationary 13C-metabolic flux analysis (INST 13C-MFA) to calculate metabolic fluxes. We found that d0 cells mainly use glutamine to fuel the TCA. Furthermore, they rely on extracellular pyruvate to allow continuous growth. This metabolic situation does not allow for mitochondrial or glycolytic spare capacity, i.e. the ability to adapt energy generation to altered needs. Accordingly, neuronal precursor cells displayed a higher sensitivity to several mitochondrial toxicants than mature neurons differentiated from them. In summary, this study shows that precursor cells lose their glutamine dependency during differentiation while they gain flexibility of energy generation and thereby increase their resistance to low concentrations of mitochondrial toxicants.

Major changes of cell function and toxicant sensitivity in cultured cells undergoing mild, quasi-natural genetic drift.

Genomic drift affects the functional properties of cell lines, and the reproducibility of data from in vitro studies. While chromosomal aberrations and mutations in single pivotal genes are well explored, little is known about effects of minor, possibly pleiotropic, genome changes. We addressed this question for the human dopaminergic neuronal precursor cell line LUHMES by comparing two subpopulations (SP) maintained either at the American-Type-Culture-Collection (ATCC) or by the original provider (UKN). Drastic differences in susceptibility towards the specific dopaminergic toxicant 1-methyl-4-phenylpyridinium (MPP+) were observed. Whole-genome sequencing was performed to identify underlying genetic differences. While both SP had normal chromosome structures, they displayed about 70 differences on the level of amino acid changing events. Some of these differences were confirmed biochemically, but none offered a direct explanation for the altered toxicant sensitivity pattern. As second approach, markers known to be relevant for the intended use of the cells were specifically tested. The "ATCC" cells rapidly down-regulated the dopamine-transporter and tyrosine-hydroxylase after differentiation, while "UKN" cells maintained functional levels. As the respective genes were not altered themselves, we conclude that polygenic complex upstream changes can have drastic effects on biochemical features and toxicological responses of relatively similar SP of cells.

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.

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.

Combination of multiple neural crest migration assays to identify environmental toxicants from a proof-of-concept chemical library.

Many in vitro tests have been developed to screen for potential neurotoxicity. However, only few cell function-based tests have been used for comparative screening, and thus experience is scarce on how to confirm and evaluate screening hits. We addressed these questions for the neural crest cell migration test (cMINC). After an initial screen, a hit follow-up strategy was devised. A library of 75 compounds plus internal controls (NTP80-list), assembled by the National Toxicology Program of the USA (NTP) was used. It contained some known classes of (developmental) neurotoxic compounds. The primary screen yielded 23 confirmed hits, which comprised ten flame retardants, seven pesticides and six drug-like compounds. Comparison of concentration-response curves for migration and viability showed that all hits were specific. The extent to which migration was inhibited was 25-90%, and two organochlorine pesticides (DDT, heptachlor) were most efficient. In the second part of this study, (1) the cMINC assay was repeated under conditions that prevent proliferation; (2) a transwell migration assay was used as a different type of migration assay; (3) cells were traced to assess cell speed. Some toxicants had largely varying effects between assays, but each hit was confirmed in at least one additional test. This comparative study allows an estimate on how confidently the primary hits from a cell function-based screen can be considered as toxicants disturbing a key neurodevelopmental process. Testing of the NTP80-list in more assays will be highly interesting to assemble a test battery and to build prediction models for developmental toxicity.

Definition of transcriptome-based indices for quantitative characterization of chemically disturbed stem cell development: introduction of the STOP-Toxukn and STOP-Toxukk tests.

Stem cell-based in vitro test systems can recapitulate specific phases of human development. In the UKK test system, human pluripotent stem cells (hPSCs) randomly differentiate into cells of the three germ layers and their derivatives. In the UKN1 test system, hPSCs differentiate into early neural precursor cells. During the normal differentiation period (14 days) of the UKK system, 570 genes [849 probe sets (PSs)] were regulated >fivefold; in the UKN1 system (6 days), 879 genes (1238 PSs) were regulated. We refer to these genes as 'developmental genes'. In the present study, we used genome-wide expression data of 12 test substances in the UKK and UKN1 test systems to understand the basic principles of how chemicals interfere with the spontaneous transcriptional development in both test systems. The set of test compounds included six histone deacetylase inhibitors (HDACis), six mercury-containing compounds ('mercurials') and thalidomide. All compounds were tested at the maximum non-cytotoxic concentration, while valproic acid and thalidomide were additionally tested over a wide range of concentrations. In total, 242 genes (252 PSs) in the UKK test system and 793 genes (1092 PSs) in the UKN1 test system were deregulated by the 12 test compounds. We identified sets of 'diagnostic genes' appropriate for the identification of the influence of HDACis or mercurials. Test compounds that interfered with the expression of developmental genes usually antagonized their spontaneous development, meaning that up-regulated developmental genes were suppressed and developmental genes whose expression normally decreases were induced. The fraction of compromised developmental genes varied widely between the test compounds, and it reached up to 60 %. To quantitatively describe disturbed development on a genome-wide basis, we recommend a concept of two indices, 'developmental potency' (D p) and 'developmental index' (D i), whereby D p is the fraction of all developmental genes that are up- or down-regulated by a test compound, and D i is the ratio of overrepresentation of developmental genes among all genes deregulated by a test compound. The use of D i makes hazard identification more sensitive because some compounds compromise the expression of only a relatively small number of genes but have a high propensity to deregulate developmental genes specifically, resulting in a low D p but a high D i. In conclusion, the concept based on the indices D p and D i offers the possibility to quantitatively express the propensity of test compounds to interfere with normal development.

Adverse outcome pathways: opportunities, limitations and open questions.

Adverse outcome pathways (AOPs) are a recent toxicological construct that connects, in a formalized, transparent and quality-controlled way, mechanistic information to apical endpoints for regulatory purposes. AOP links a molecular initiating event (MIE) to the adverse outcome (AO) via key events (KE), in a way specified by key event relationships (KER). Although this approach to formalize mechanistic toxicological information only started in 2010, over 200 AOPs have already been established. At this stage, new requirements arise, such as the need for harmonization and re-assessment, for continuous updating, as well as for alerting about pitfalls, misuses and limits of applicability. In this review, the history of the AOP concept and its most prominent strengths are discussed, including the advantages of a formalized approach, the systematic collection of weight of evidence, the linkage of mechanisms to apical end points, the examination of the plausibility of epidemiological data, the identification of critical knowledge gaps and the design of mechanistic test methods. To prepare the ground for a broadened and appropriate use of AOPs, some widespread misconceptions are explained. Moreover, potential weaknesses and shortcomings of the current AOP rule set are addressed (1) to facilitate the discussion on its further evolution and (2) to better define appropriate vs. less suitable application areas. Exemplary toxicological studies are presented to discuss the linearity assumptions of AOP, the management of event modifiers and compensatory mechanisms, and whether a separation of toxicodynamics from toxicokinetics including metabolism is possible in the framework of pathway plasticity. Suggestions on how to compromise between different needs of AOP stakeholders have been added. A clear definition of open questions and limitations is provided to encourage further progress in the field.

Identification of transcriptome signatures and biomarkers specific for potential developmental toxicants inhibiting human neural crest cell migration.

The in vitro test battery of the European research consortium ESNATS ('novel stem cell-based test systems') has been used to screen for potential human developmental toxicants. As part of this effort, the migration of neural crest (MINC) assay has been used to evaluate chemical effects on neural crest function. It identified some drug-like compounds in addition to known environmental toxicants. The hits included the HSP90 inhibitor geldanamycin, the chemotherapeutic arsenic trioxide, the flame-retardant PBDE-99, the pesticide triadimefon and the histone deacetylase inhibitors valproic acid and trichostatin A. Transcriptome changes triggered by these substances in human neural crest cells were recorded and analysed here to answer three questions: (1) can toxicants be individually identified based on their transcript profile; (2) how can the toxicity pattern reflected by transcript changes be compacted/dimensionality-reduced for practical regulatory use; (3) how can a reduced set of biomarkers be selected for large-scale follow-up? Transcript profiling allowed clear separation of different toxicants and the identification of toxicant types in a blinded test study. We also developed a diagrammatic system to visualize and compare toxicity patterns of a group of chemicals by giving a quantitative overview of altered superordinate biological processes (e.g. activation of KEGG pathways or overrepresentation of gene ontology terms). The transcript data were mined for potential markers of toxicity, and 39 transcripts were selected to either indicate general developmental toxicity or distinguish compounds with different modes-of-action in read-across. In summary, we found inclusion of transcriptome data to largely increase the information from the MINC phenotypic test.

A transcriptome-based classifier to identify developmental toxicants by stem cell testing: design, validation and optimization for histone deacetylase inhibitors.

Test systems to identify developmental toxicants are urgently needed. A combination of human stem cell technology and transcriptome analysis was to provide a proof of concept that toxicants with a related mode of action can be identified and grouped for read-across. We chose a test system of developmental toxicity, related to the generation of neuroectoderm from pluripotent stem cells (UKN1), and exposed cells for 6 days to the histone deacetylase inhibitors (HDACi) valproic acid, trichostatin A, vorinostat, belinostat, panobinostat and entinostat. To provide insight into their toxic action, we identified HDACi consensus genes, assigned them to superordinate biological processes and mapped them to a human transcription factor network constructed from hundreds of transcriptome data sets. We also tested a heterogeneous group of 'mercurials' (methylmercury, thimerosal, mercury(II)chloride, mercury(II)bromide, 4-chloromercuribenzoic acid, phenylmercuric acid). Microarray data were compared at the highest non-cytotoxic concentration for all 12 toxicants. A support vector machine (SVM)-based classifier predicted all HDACi correctly. For validation, the classifier was applied to legacy data sets of HDACi, and for each exposure situation, the SVM predictions correlated with the developmental toxicity. Finally, optimization of the classifier based on 100 probe sets showed that eight genes (F2RL2, TFAP2B, EDNRA, FOXD3, SIX3, MT1E, ETS1 and LHX2) are sufficient to separate HDACi from mercurials. Our data demonstrate how human stem cells and transcriptome analysis can be combined for mechanistic grouping and prediction of toxicants. Extension of this concept to mechanisms beyond HDACi would allow prediction of human developmental toxicity hazard of unknown compounds with the UKN1 test system.

Epigenetic changes and disturbed neural development in a human embryonic stem cell-based model relating to the fetal valproate syndrome.

Exposure to the antiepileptic drug valproic acid (VPA) during gestation causes neurofunctional and anatomic deficits in later life. At present, there are little human data on how early neural development is affected by chemicals. We used human embryonic stem cells, differentiating to neuroectodermal precursors, as a model to investigate the modes of action of VPA. Microarray expression profiling, qPCR of specific marker genes, immunostaining and the expression of green fluorescent protein under the control of the promoter of the canonical neural precursor cell marker HES5 were used as readouts. Exposure to VPA resulted in distorted marker gene expression, characterized by a relative increase in NANOG and OCT4 and a reduction in PAX6. A similar response pattern was observed with trichostatin A, a potent and specific histone deacetylase inhibitor (HDACi), but not with several other toxicants. Differentiation markers were disturbed by prolonged, but not by acute treatment with HDACi, and the strongest disturbance of differentiation was observed by toxicant exposure during early neural fate decision. The increased acetylation of histones observed in the presence of HDACi may explain the up-regulation of some genes. However, to understand the down-regulation of PAX6 and the overall complex transcript changes, we examined further epigenetic markers. Alterations in the methylation of lysines 4 and 27 of histone H3 were detected in the promoter region of PAX6 and OCT4. The changes in these activating and silencing histone marks provide a more general mechanistic rational for the regulation of developmentally important genes at non-cytotoxic drug concentrations.

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

Information on design principles governing transcriptome changes upon transition from safe to hazardous drug concentrations or from tolerated to cytotoxic drug levels are important for the application of toxicogenomics data in developmental toxicology. Here, we tested the effect of eight concentrations of valproic acid (VPA; 25-1000 µM) in an assay that recapitulates the development of human embryonic stem cells to neuroectoderm. Cells were exposed to the drug during the entire differentiation process, and the number of differentially regulated genes increased continuously over the concentration range from zero to about 3000. We identified overrepresented transcription factor binding sites (TFBS) as well as superordinate cell biological processes, and we developed a gene ontology (GO) activation profiler, as well as a two-dimensional teratogenicity index. Analysis of the transcriptome data set by the above biostatistical and systems biology approaches yielded the following insights: (i) tolerated (≤25 µM), deregulated/teratogenic (150-550 µM), and cytotoxic (≥800 µM) concentrations could be differentiated. (ii) Biological signatures related to the mode of action of VPA, such as protein acetylation, developmental changes, and cell migration, emerged from the teratogenic concentrations range. (iii) Cytotoxicity was not accompanied by signatures of newly emerging canonical cell death/stress indicators, but by catabolism and decreased expression of cell cycle associated genes. (iv) Most, but not all of the GO groups and TFBS seen at the highest concentrations were already overrepresented at 350-450 µM. (v) The teratogenicity index reflected this behavior, and thus differed strongly from cytotoxicity. Our findings suggest the use of the highest noncytotoxic drug concentration for gene array toxicogenomics studies, as higher concentrations possibly yield wrong information on the mode of action, and lower drug levels result in decreased gene expression changes and thus a reduced power of the study.

Chromatin dynamics during epigenetic reprogramming in the mouse germ line.

A unique feature of the germ cell lineage is the generation of totipotency. A critical event in this context is DNA demethylation and the erasure of parental imprints in mouse primordial germ cells (PGCs) on embryonic day 11.5 (E11.5) after they enter into the developing gonads. Little is yet known about the mechanism involved, except that it is apparently an active process. We have examined the associated changes in the chromatin to gain further insights into this reprogramming event. Here we show that the chromatin changes occur in two steps. The first changes in nascent PGCs at E8.5 establish a distinctive chromatin signature that is reminiscent of pluripotency. Next, when PGCs are residing in the gonads, major changes occur in nuclear architecture accompanied by an extensive erasure of several histone modifications and exchange of histone variants. Furthermore, the histone chaperones HIRA and NAP-1 (NAP111), which are implicated in histone exchange, accumulate in PGC nuclei undergoing reprogramming. We therefore suggest that the mechanism of histone replacement is critical for these chromatin rearrangements to occur. The marked chromatin changes are intimately linked with genome-wide DNA demethylation. On the basis of the timing of the observed events, we propose that if DNA demethylation entails a DNA repair-based mechanism, the evident histone replacement would represent a repair-induced response event rather than being a prerequisite.