Joshua F. Robinson

Discriminating classes of developmental toxicants using gene expression profiling in the embryonic stem cell test

The embryonic stem cell test (EST) has been shown to be a promising in vitro method for the prediction of developmental toxicity. In our previous studies, we demonstrated that the implementation of gene expression analysis in the EST may further improve the identification of developmental toxicants. In the present study, we investigated if gene expression profiling could be used to discriminate compound classes with distinct modes of action (MoA) using the EST protocol. Gene expression data of our previous study were used and were analyzed of embryonic stem cell (ESC) differentiation cultures exposed to six compounds belonging to two classes with distinct MoA, namely phthalates and triazoles. We used three approaches to study class-characteristic gene regulation that may be useful for discrimination of compound classes. First, at the individual gene level, gene signatures characteristic for each class were identified that successfully discriminated both classes using principal component analysis. Second, at the functional level, enriched gene ontology (GO) biological processes showed their usefulness for class discrimination via hierarchical clustering. Third, two previously identified gene sets, which we designed to predict developmental toxicity, appeared successful in separating phthalate from triazole compounds. In summary, we established the possibility to discriminate between compound classes in the EST system using three different specific transcriptomics-based approaches. Differential gene expression information specific for the class of compound under study may be employed to optimize prioritization of compounds within that class for further testing.

Time-response evaluation by transcriptomics of methylmercury effects on neural differentiation of murine embryonic stem cells

Current globally harmonized Organisation for Economic Co-operation and Development (OECD) animal test guidelines for developmental toxicity require high numbers of experimental animals. To reduce animal use in this field, alternative developmental toxicity assays are highly desirable. We previously developed a dynamic in vitro model for screening effects of possible neurodevelopmental toxicants, using neural cell differentiation of pluripotent murine embryonic stem cells. To further mechanistically characterize the mouse neural embryonic stem cell test (ESTn) and to improve detection of possible neurodevelopmental toxicants, gene expression patterns were studied describing neural cell differentiation over time, as well as the impact on gene expression of exposure to the well-known neurotoxicant methylmercury (MeHg). A transcriptomics study was performed to examine whole-genome expression changes during the first 7 days of the cell differentiation protocol. Specific gene clusters were identified and enrichment analysis of Gene Ontology (GO) terms and gene sets derived from literature was performed using DAVID and T-profiler. Over time, a decrease of blastocyst and trophectoderm GO terms was observed, which included well-characterized pluripotency genes. Furthermore, an increase in the range of neural development-related GO terms, such as neuron differentiation and the wnt pathway, was observed. Analysis of gene expression using principle component analysis showed a time-dependent track in untreated cells, describing the process of neural differentiation. Furthermore, MeHg was shown to induce deviation from the predefined differentiation track. The compound inhibited general development GO terms and induced neural GO terms over time. This system appears promising for studying compound effects on neural differentiation in a mechanistic approach.

A system-based comparison of gene expression reveals alterations in oxidative stress, disruption of ubiquitin-proteasome system and altered cell cycle regulation after exposure to cadmium and methylmercury in mouse embryonic fibroblast.

Environmental and occupational exposures to heavy metals such as methylmercury (MeHg) and cadmium (Cd) pose significant health risks to humans, including neurotoxicity. The underlying mechanisms of their toxicity, however, remain to be fully characterized. Our previous studies with Cd and MeHg have demonstrated that the perturbation of the ubiquitin-proteasome system (UPS) was associated with metal-induced cytotoxicity and apoptosis. We conducted a microarray-based gene expression analysis to compare metal-altered gene expression patterns with a classical proteasome inhibitor, MG132 (0.5 microM), to determine whether the disruption of the UPS is a critical mechanism of metal-induced toxicity. We treated mouse embryonic fibroblast cells at doses of MeHg (2.5 microM) and Cd (5.0 microM) for 24 h. The doses selected were based on the neutral red-based cell viability assay where initial statistically significant decreases in variability were detected. Following normalization of the array data, we employed multilevel analysis tools to explore the data, including group comparisons, cluster analysis, gene annotations analysis (gene ontology analysis), and pathway analysis using GenMAPP and Ingenuity Pathway Analysis (IPA). Using these integrated approaches, we identified significant gene expression changes across treatments within the UPS (Uchl1 and Ube2c), antioxidant and phase II enzymes (Gsta2, Gsta4, and Noq1), and genes involved in cell cycle regulation pathways (ccnb1, cdc2a, and cdc25c). Furthermore, pathway analysis revealed significant alterations in genes implicated in Parkinsons disease pathogenesis following metal exposure. This study suggests that these pathways play a critical role in the development of adverse effects associated with metal exposures.

Cadmium-Induced Differential Toxicogenomic Response in Resistant and Sensitive Mouse Strains Undergoing Neurulation

Common inbred mouse strains, such as the C57BL/6 (C57) and the SWV, display differences in sensitivity to environmental teratogens during gestation. For example, the C57 is more sensitive than the SWV to cadmium (Cd) exposure during neurulation, inducing a higher incidence of neural tube defects (NTDs). Here, we report, using Cd as a model teratogen, the first large scale toxicogenomic study to compare teratogen-induced gene expression alterations in C57 and SWV embryos undergoing neurulation, identifying toxicogenomic responses that associate with developmental toxicity and differential sensitivity. Using a systems-based toxicogenomic approach, comparing Cd-exposed and control C57 and SWV embryos (12- and 24-h postinjection [p.i.] [gestational day 8.0, ip]), we examined differentially expressed genes at multiple levels (biological process, pathway, gene) using Gene Ontology (GO) analysis, pathway mapping and cross-scatter plots. In both C57 and SWV embryos, we observed several gene expression alterations linked with cell cycle-related classifications, however, only in the C57 we observed upregulation of p53-dependent mediators Ccng1 and Pmaip1, previously associated with cell cycle arrest, apoptosis and NTD formation. In addition, we also identified a greater reduction in expression of nervous system development-related genes (e.g., Zic1, En2, Neurog1, Elavl4, Metrn, Nr2f1, Nr2f2) in the C57 compared to the SWV (12-h p.i.). In summary, our results indicate that differences in Cd-induced gene expression profiles between NTD resistant and sensitive strains within enriched biological processes (including developmental and cell cycle-related categories) associate with increased sensitivity to developmental toxicity as determined by observations of increased NTD formation, mortality (resorptions) and reduced fetal growth. Such observations may provide more detailed and useful mechanistic clues for identification of differences in life-stage specific teratogenic response.

Transcriptomic concentration-response evaluation of valproic acid, cyproconazole, and hexaconazole in the neural embryonic stem cell test (ESTn).

Alternative developmental toxicity assays are urgently needed to reduce animal use in regulatory developmental toxicology. We previously designed an in vitro murine neural embryonic stem cell test (ESTn) as a model for neurodevelopmental toxicity testing (Theunissen et al., 2010). Toxicogenomic approaches have been suggested for incorporation into the ESTn to further increase predictivity and to provide mechanistic insights. Therefore, in this study, using a transcriptomic approach, we investigated the concentration-dependent effects of three known (neuro) developmental toxicants, two triazoles, cyproconazole (CYP) and hexaconazole (HEX), and the anticonvulsant valproic acid (VPA). Compound effects on gene expression during neural differentiation and corresponding regulated gene ontology (GO) terms were identified after 24 h of exposure in relation to morphological changes on day 11 of culture. Concentration-dependent responses on individual gene expression and on biological processes were determined for each compound, providing information on mechanism and concentration-response characteristics. All compounds caused enrichment of the embryonic development process. CYP and VPA but not HEX significantly enriched the neuron development process. Furthermore, specific responses for triazole compounds and VPA were observed within the GO-term sterol metabolic process. The incorporation of transcriptomics in the ESTn was shown to enable detection of effects, which precede morphological changes and provide a more sensitive measure of concentration-dependent effects as compared with classical morphological assessments. Furthermore, mechanistic insight can be instrumental in the extrapolation of effects in the ESTn to human hazard assessment.

Embryotoxicant-specific transcriptomic responses in rat postimplantation whole-embryo culture.

Rat postimplantation whole-embryo culture (WEC) is a promising alternative test for the assessment of developmental toxicity. Toxicogenomic-based approaches may improve the predictive ability of the WEC model by providing a means to identify compound-specific mechanistic responses associated with embryotoxicity in vivo. Furthermore, alterations in gene expression may serve as a sensitive, objective, and robust marker, which precedes the observation of classical developmental toxicity endpoints in time. In this study, in combination with morphological developmental assessments, we studied transcriptomic responses associated with four distinct teratogens (caffeine [CAF], methylmercury [MM], monobutyl phthalate, and methoxyacetic acid) after 4 h of exposure, well before apparent embryotoxicity in WEC. We evaluated gene expression changes associated with similar levels of induced morphological embryotoxicity for each teratogen (as determined by total morphological score), evaluating for functional enrichment and quantitative changes in response. Concentrations selected for each of the four teratogens used induced a number of common effects on embryonic development (neural tube closure and optic/otic system). Despite inducing common morphological effects, our analysis suggests limited overlap in terms of toxicogenomic response at the gene expression level and at the level of biological processes across all four test chemicals. Many unique responses associated with each chemical correlated with previously hypothesized modes of developmental toxicity. For example, alterations in developmental signaling and cholesterol metabolism were observed with MM and CAF, respectively. This initial study suggests that distinct chemically induced toxicogenomic responses precede morphological effects in WEC and that these responses are relevant with mechanisms of toxicity previously observed in vivo.

Toxicogenomic profiling in maternal and fetal rodent brains following gestational exposure to chlorpyrifos

Considering the wide variety of effects that have been reported to occur in the developmental neurotoxicity of chlorpyrifos (CP) and the lack of consensus on their dependence of brain acetylcholinesterase (AChE) activity inhibition, we applied microarray technology to explore dose-dependent alterations in transcriptional response in the fetal and maternal C57BL/6 mouse brain after daily gestational exposure (days 6 to 17) to CP (2, 4, 10, 12 or 15 mg/kg, sc). We identified significantly altered genes across doses and assessed for overrepresentation of Gene Ontology (GO) biological processes and KEGG pathways. We further clustered genes based on their expression profiles across doses and repeated the GO/pathways analysis for each cluster. The dose-effect relationship of CP on gene expression, both at the gene and pathway levels was non-monotonic and not necessarily related to brain AChE inhibition. The largest impact was observed in the 10mg/kg dose group which was also the LOAEL for brain AChE inhibition. In the maternal brain, lower doses (4 mg/kg) influenced GO categories and pathways such as cell adhesion, behavior, lipid metabolism, long-term potentiation, nervous system development, neurogenesis, synaptic transmission. In the fetal brain, lower doses (2 and/or 4 mg/kg) significantly altered cell division, translation, transmission of nerve impulse, chromatin modification, long-term potentiation. In addition, some genes involved in nervous system development and signaling were shown to be specifically influenced by these lower CP doses. Our approach was sensitive and reflected the diversity of responses known to be disrupted by CP and highlighted possible additional consequences of CP neurotoxicity, such as disturbance of the ubiquitin proteasome system.

A Comparison of Gene Expression Responses in Rat Whole Embryo Culture and In Vivo: Time-Dependent Retinoic Acid-Induced Teratogenic Response

The whole embryo culture (WEC) model serves as a potential alternative for classical in vivo developmental toxicity testing. In the WEC, cultured rat embryos are exposed during neurulation and early organogenesis and evaluated for morphological effects. Toxicogenomic-based approaches may improve the predictive ability of WEC by providing molecular-based markers associated with chemical exposure, which can be compared across multiple parameters (e.g., exposure duration, developmental time, experimental model). Additionally, comparisons between in vitro and in vivo models may identify objective relevant molecular responses linked with developmental toxicity endpoints in vivo. In this study, using a transcriptomic approach, we compared all-trans retinoic acid (RA)-exposed and nonexposed Wistar rat embryos derived using WEC (RA, 0.5 μg/ml) or in vivo (RA, 50 mg/kg, oral gavage) to identify overlapping and nonoverlapping effects of RA on RNA expression in parallel with morphological changes. Across six time points (gestational day 10 + 2-48 h), we observed strong similarities in RA response at the gene (directionality, significance) and functional (e.g., embryonic development, cell differentiation) level which associated with RA-induced adverse morphological effects, including growth reduction as well as alterations in neural tube, limb, branchial, and mandible development. We observed differences between models in the timing of RA-induced effects on genes related to embryonic development and RA metabolism. These observations on the gene expression level were associated with specific differential morphological outcomes. This study supports the use of WEC to examine compound-induced molecular responses relative to in vivo and, furthermore, assists in defining the applicability domain of the WEC in determining complementary windows of sensitivity for developmental toxicological investigations.

Arsenic- and Cadmium-Induced Toxicogenomic Response in Mouse Embryos Undergoing Neurulation

Arsenic (As) and cadmium (Cd) are well-characterized teratogens in animal models inducing embryotoxicity and neural tube defects (NTDs) when exposed during neurulation. Toxicological research is needed to resolve the specific biological processes and associated molecular pathways underlying metal-induced toxicity during this timeframe in gestational development. In this study, we investigated the dose-dependent effects of As and Cd on gene expression in C57BL/6J mouse embryos exposed in utero during neurulation (GD8) to identify significantly altered genes and corresponding biological processes associated with embryotoxicity. We quantitatively examined the toxicogenomic dose-response relationship at the gene level. Our results suggest that As and Cd induce dose-dependent gene expression alterations representing shared (cell cycle, response to UV, glutathione metabolism, RNA processing) and unique (alcohol/sugar metabolism) biological processes, which serve as robust indicators of metal-induced developmental toxicity and indicate underlying embryotoxic effects. Our observations also correlate well with previously identified impacts of As and Cd on specific genes associated with metal-induced toxicity (Cdkn1a, Mt1). In summary, we have identified in a quantitative manner As and Cd induced dose-dependent effects on gene expression in mouse embryos during a peak window of sensitivity to embryotoxicity and NTDs in the sensitive C57BL/6J strain.

Comparison of MeHg-induced toxicogenomic responses across in vivo and in vitro models used in developmental toxicology.

Toxicogenomic evaluations may improve toxicity prediction of in vitro-based developmental models, such as whole embryo culture (WEC) and embryonic stem cells (ESC), by providing a robust mechanistic marker which can be linked with responses associated with developmental toxicity in vivo. While promising in theory, toxicogenomic comparisons between in vivo and in vitro models are complex due to inherent differences in model characteristics and experimental design. Determining factors which influence these global comparisons are critical in the identification of reliable mechanistic-based markers of developmental toxicity. In this study, we compared available toxicogenomic data assessing the impact of the known teratogen, methylmercury (MeHg) across a diverse set of in vitro and in vivo models to investigate the impact of experimental variables (i.e. model, dose, time) on our comparative assessments. We evaluated common and unique aspects at both the functional (Gene Ontology) and gene level of MeHg-induced response. At the functional level, we observed stronger similarity in MeHg-response between mouse embryos exposed in utero (2 studies), ESC, and WEC as compared to liver, brain and mouse embryonic fibroblast MeHg studies. These findings were strongly correlated to the presence of a MeHg-induced developmentally related gene signature. In addition, we identified specific MeHg-induced gene expression alterations associated with developmental signaling and heart development across WEC, ESC and in vivo systems. However, the significance of overlap between studies was highly dependent on traditional experimental variables (i.e. dose, time). In summary, we identify promising examples of unique gene expression responses which show in vitro-in vivo similarities supporting the relevance of in vitro developmental models for predicting in vivo developmental toxicity.