Utilizing toxicogenomic data to understand chemical mechanism of action in risk assessment

Functional differences and similarities in activated peripheral blood mononuclear cells by lipopolysaccharide or phytohemagglutinin stimulation between human and cynomolgus monkeys

Background

The monkey is a primary species used in toxicological research. However, the failures of preclinical studies to predict a life-threatening “cytokine storm”, which, for instance, rapidly occurred in six healthy volunteers with the CD28 superagonist monoclonal antibody (mAb) TGN1412 in the first-in-human phase I clinical trial, have emphasized a need to clarify the differences between human and monkey immune systems.

Methods

In the present study, we analyzed and compared the lymphocyte proliferation, cytokine secretion, and gene expression profiles after phytohemagglutinin (PHA) and lipopolysaccharide (LPS) stimulation of peripheral blood mononuclear cells (PBMCs) from three healthy humans and cynomolgus monkeys (Macaca fascicularis).

Results

The results derived from comparison with the corresponding control groups showed that PHA in humans induced a stronger proliferation and wider range of cytokine secretion, along with a greater number of differently expressed genes (DEGs), than when PHA was applied in cynomolgus monkeys. The significant upregulation of genes involved in the mitotic cell cycle, including cyclin B2, TOP2A, TYMS, and CEP55, was observed in human PBMCs with PHA stimulation, while only infrequent or slight upregulation occurred in cynomolgus monkey PBMCs, which may be one of the reasons for a stronger response to PHA in humans. In contrast to PHA, LPS in both species induced a similar proliferation ratio, cytokine profile, and DEG count, suggesting that human and cynomolgus monkeys have a similar response intensity for innate immune responses. Furthermore, 38 and 20 overlapped genes under PHA and LPS stimulation, respectively, were found in both species. These overlapped DEGs were associated with the same biological functions, including DNA replication, mitosis, immune response, chemotaxis, and inflammatory response. Thus, these results might reflect the highly conserved signatures of immune responses to PHA/LPS stimulation across the primates. Moreover, there were some differences in antigen processing and presentation, and the interferon gamma (INF-γ)–mediated signaling pathway in these species detected by gene expression profile study.

Conclusions

In conclusion, this is the first study to compare data on the responses of PBMCs to PHA and LPS in humans versus cynomolgus monkeys, and these findings may provide crucial insights into translating non-human primate (NHP) studies into human trials.

Evaluation of seeds ethanolic extracts of Triplaris gardneriana Wedd. using in vitro and in vivo toxicological methods

Triplaris gardneriana Wedd. is a tree used in folk medicine to treat venereal diseases and inflammation as well as a source of biological compounds with antioxidant capacity. In order to assess the safety of these bioactive compounds, the present study aimed to determine the toxicity of an ethanolic extract of T. gardneriana, (EETg). Toxicological tests included hemolytic activity, toxicity toward the brine shrimp Artemia, cytotoxicity against breast cancer cells (MCF7) and acute oral toxicity in rodents. In addition, toxicogenomics techniques were used to determine genome expression in MCF7 cells exposed to EETg. The results showed that the extract exhibits approximately 60% of hemolytic activity at the highest tested concentration (64 µg/ml) and toxicity against nauplii of Artemia sp. (LC₅₀ of 67.85 µg/ml). Further, EETg appears to be cytotoxic to MCF7 (cell viability reduced to 40% at 250 µg/ml after 24 hr). Genomic data demonstrated differential expression of 14 genes. Data analysis indicated possible altered pathways (e.g., xenobiotic metabolism), possible adverse health risks (e.g., hepatotoxicity), and drugs with similar gene expression profile (e.g., antimicrobials). The investigation provides important information on potentially adverse aspects of EETg, which need to be considered prior to the therapeutic utilization of this plant.Abbreviations: EETg: ethanolic extract of T. gardneriana seeds; MCF7: michigan cancer foundation-7 which refers to a human breast cell line (adenocarcinoma); NGS: next-generation sequencing; edgeR: empirical analysis of digital gene expression data in R; Consensus: consensus path database; FDR: false discovery rate; NCBI: national center for biotechnology information; KEGG: kyoto encyclopedia of genes and genomes; Ingenuity: ingenuity pathway analysis software; CMAP: connectivity map; OECD: organization for economic co-operation and development; HL-60: human promyelocytic leukemia cells; PC3: prostate cancer cells.

Next generation sequencing data for use in risk assessment

Next generation sequencing (NGS) represents several powerful platforms that have revolutionized RNA and DNA analysis. The parallel sequencing of millions of DNA molecules can provide mechanistic insights into toxicology and provide new avenues for biomarker discovery with growing relevance for risk assessment. The evolution of NGS technologies has improved over the last decade with increased sensitivity and accuracy to foster new biomarker assays from tissue, blood and other biofluids. NGS sequencing technologies can identify transcriptional changes and genomic targets with base pair precision in response to chemical exposure. Further, there are several exciting movements within the toxicology community that incorporate NGS platforms into new strategies for more rapid toxicological characterizations. These include the Tox21 in vitro high throughput transcriptomic screening program, development of organotypic spheroids, alternative animal models, mining archival tissues, liquid biopsy and epigenomics. This review will describe NGS-based technologies, demonstrate how they can be used as tools for target discovery in tissue and blood, and suggest how they might be applied for risk assessment.

Variation in transcriptional responses to copper exposure across Daphnia pulex lineages

Copper pollution is pervasive in aquatic habitats and is particularly harmful to invertebrates sensitive to environmental changes such as Daphnia pulex. Mechanisms of toxicity and tolerance to copper are not well understood. We used RNA-sequencing to investigate these mechanisms in three genetically distinct D. pulex clonal lineages with different histories of copper exposure. Upregulated genes after copper exposure were enriched with Gene Ontology (GO) categories involved in digestion, molting and growth, whereas downregulated genes after copper exposure were enriched in the metal-regulatory system, immune response and epigenetic modifications. The three D. pulex clones in our study show largely similar transcriptional patterns in response to copper, with only a total of twenty genes differentially expressed in a single clonal lineages. We also detected lower relative expression of some genes known to be important for copper tolerance, metallothionein and glutathione-S-transferase, in a sensitive lineage sampled from an uncontaminated habitat. Daphnia-specific genes (without orthologs outside the genus) and Daphnia-specific duplications (genes duplicated in the Daphnia lineage) were overrepresented in differentially expressed genes, highlighting an important role for newly emerged genes in tolerating environmental stressors. The results indicate that the D. pulex lineages tested in this study generally respond to copper stress using the same major pathways, but that the more resistant clone with previous copper exposure might be better able to regulate key genes. This finding highlights the important nuances in gene expression among clones, shaped by historical exposure and influencing copper tolerance.

Screening Strategies and Methods for Better Off-Target Liability Prediction and Identification of Small-Molecule Pharmaceuticals

Pharmaceutical discovery and development is a long and expensive process that, unfortunately, still results in a low success rate, with drug safety continuing to be a major impedance. Improved safety screening strategies and methods are needed to more effectively fill this critical gap. Recent advances in informatics are now making it possible to manage bigger data sets and integrate multiple sources of screening data in a manner that can potentially improve the selection of higher-quality drug candidates. Integrated screening paradigms have become the norm in Pharma, both in discovery screening and in the identification of off-target toxicity mechanisms during later-stage development. Furthermore, advances in computational methods are making in silico screens more relevant and suggest that they may represent a feasible option for augmenting the current screening paradigm. This paper outlines several fundamental methods of the current drug screening processes across Pharma and emerging techniques/technologies that promise to improve molecule selection. In addition, the authors discuss integrated screening strategies and provide examples of advanced screening paradigms.

Transcriptomic study of the toxic mechanism triggered by beauvericin in Jurkat cells

Beauvericin (BEA), an ionophoric cyclic hexadepsipeptide mycotoxin, is able to increase oxidative stress by altering membrane ion permeability and uncoupling oxidative phosphorylation. A toxicogenomic study was performed to investigate gene expression changes triggered by BEA exposure (1.5, 3 and 5 μM; 24 h) in Jurkat cells through RNA-sequencing and differential gene expression analysis. Perturbed gene expression was observed in a concentration dependent manner, with 43 differentially expressed genes (DEGs) overlapped in the three studied concentrations. Gene ontology (GO) analysis showed several biological processes related to electron transport chain, oxidative phosphorylation, and cellular respiration significantly altered. Molecular functions linked to mitochondrial respiratory chain and oxidoreductase activity were over-represented (q-value < 0.01). Pathway analysis revealed oxidative phosphorylation and electron transport chain as the most significantly altered pathways in all studied doses (z-score > 1.96; adj p-value < 0.05). 77 genes involved in the respiratory chain were significantly down-regulated at least at one dose. Moreover, 21 genes related to apoptosis and programmed cell death, and 12 genes related to caspase activity were significantly altered, mainly affecting initiator caspases 8, 9 and 10. The results demonstrated BEA-induced mitochondrial damage affecting the respiratory chain, and pointing to apoptosis through the caspase cascade in human lymphoblastic T cells.

Transcriptional profiling reveals gene expression changes associated with inflammation and cell proliferation following short-term inhalation exposure to copper oxide nanoparticles

Our recent studies revealed a dose-dependent proinflammatory response to copper oxide nanoparticles (CuO NPs) in rats following short-term inhalation exposure for five consecutive days. Here transcriptomics approaches were applied using the same model to assess global gene expression in lung tissues obtained 1 day post-exposure and after a recovery period of 22 days from rats exposed to clean air or 6 hour equivalent doses of 3.3 mg m^-3 (low dose) and 13.2 mg m^-3 (high dose). Microarray analyses yielded about 1000 differentially expressed genes in the high-dose group and 200 in low-dose compared to the clean air control group, and less than 20 after the recovery period. Pathway analysis indicated cell proliferation/survival and inflammation as the main processes triggered by exposure to CuO NPs. We did not find significant perturbations of pathways related to oxidative stress. Upregulation of epithelial cell transforming protein 2 (Ect2), a known oncogene, was noted and ECT2 protein was upregulated in the lungs of exposed animals. Proliferation of alveolar epithelial cells was demonstrated based on Ki67 expression. The gene encoding monocyte chemoattractant protein 1 (or CCL2) was also upregulated and this was confirmed by immunohistochemistry. However, no aberrant DNA methylation of inflammation-associated genes was observed. In conclusion, we have found that inhalation of CuO NPs in rats causes upregulation of the oncoprotein ECT2 and the chemokine CCL2 and other proinflammatory markers as well as proliferation in bronchoalveolar epithelium after a short-term inhalation exposure. Thus, pathways known to be associated with neoplastic processes and inflammation were affected in this model.

Integration of the TGx-28.65 genomic biomarker with the flow cytometry micronucleus test to assess the genotoxicity of disperse orange and 1,2,4-benzenetriol in human TK6 cells

In vitro gene expression signatures to predict toxicological responses can provide mechanistic context for regulatory testing. We previously developed the TGx-28.65 genomic biomarker from a database of gene expression profiles derived from human TK6 cells exposed to 28 well-known compounds. The biomarker comprises 65 genes that can classify chemicals as DNA damaging or non-DNA damaging. In this study, we applied the TGx-28.65 genomic biomarker in parallel with the in vitro micronucleus (MN) assay to determine if two chemicals of regulatory interest at Health Canada, disperse orange (DO: the orange azo dye 3-[[4-[(4-Nitrophenyl)azo]phenyl] benzylamino]propanenitrile) and 1,2,4-benzenetriol (BT: a metabolite of benzene) are genotoxic or non-genotoxic. Both chemicals caused dose-dependent declines in relative survival and increases in apoptosis. A strong significant increase in MN induction was observed for all concentrations of BT; the top two concentrations of DO also caused a statistically significant increase in MN, but these increases were <2-fold above controls. TGx-28.65 analysis classified BT as genotoxic at all three concentrations and DO as genotoxic at the mid and high concentrations. Thus, although DO only caused a small increase in MN, this response was sufficient to induce a cellular DNA damage response. Benchmark dose modeling confirmed that BT is much more potent than DO. The results strongly suggest that follow-up work is required to assess whether DO and BT are also genotoxic in vivo. This is particularly important for DO, which may require metabolic activation by bacterial gut flora to fully induce its genotoxic potential. Our previously published data and this proof of concept study suggest that the TGx-28.65 genomic biomarker has the potential to add significant value to existing approaches used to assess genotoxicity.

Barriers to the use of toxicogenomics data in human health risk assessment: A survey of Canadian risk assessors

Regulatory agencies worldwide need to modernize human health risk assessment (HHRA) to meet challenges of the 21st century. Toxicogenomics is at the core of this improvement. Today, however, the use of toxicogenomics data in HHRA is very limited. The purpose of this survey was to identify barriers to the application of toxicogenomics data in HHRA by human health risk assessors. An online survey targeting Canadian risk assessors gathered information on their knowledge and perception of toxicogenomics, their current and future inclusion of toxicogenomics data in HHRA, and barriers to the use of such data. Twenty-nine (29) participants completed a questionnaire after 2 months of solicitation. The results show that the application of toxicogenomics data in Canada is marginal, with 85% of respondents reporting that they never or rarely used such data. Knowledge of toxicogenomics by Canadian risk assessors is also limited: about two-thirds of respondents (68%) were not at all or only slightly familiar with the concept. Lack of guidelines for toxicogenomics data interpretation, data quality assessment and on their use in HHRA, were found to be major barriers. In conclusion, there is a need for interventions aimed at facilitating the use of toxicogenomics data in HHRA, when available.

Systems Biology to Support Nanomaterial Grouping

The assessment of potential health risks of engineered nanomaterials (ENMs) is a challenging task due to the high number and great variety of already existing and newly emerging ENMs. Reliable grouping or categorization of ENMs with respect to hazards could help to facilitate prioritization and decision making for regulatory purposes. The development of grouping criteria, however, requires a broad and comprehensive data basis. A promising platform addressing this challenge is the systems biology approach. The different areas of systems biology, most prominently transcriptomics, proteomics and metabolomics, each of which provide a wealth of data that can be used to reveal novel biomarkers and biological pathways involved in the mode-of-action of ENMs. Combining such data with classical toxicological data would enable a more comprehensive understanding and hence might lead to more powerful and reliable prediction models. Physico-chemical data provide crucial information on the ENMs and need to be integrated, too. Overall statistical analysis should reveal robust grouping and categorization criteria and may ultimately help to identify meaningful biomarkers and biological pathways that sufficiently characterize the corresponding ENM subgroups. This chapter aims to give an overview on the different systems biology technologies and their current applications in the field of nanotoxicology, as well as to identify the existing challenges.

Transcriptomic responses in the oral cavity of F344 rats and B6C3F1 mice following exposure to Cr(VI): Implications for risk assessment

Exposure to hexavalent chromium [Cr(VI)] in drinking water was previously reported to increase oral tumor incidence in F344 rats. To investigate the mode of action for these tumors, transcriptomic profiles in oral mucosa samples of F344 rats and B6C3F1 mice were analyzed following exposure to 0.1-180 ppm Cr(VI) for 7 or 90 days. In rats, genome-wide microarray analyses identified no significantly differentially expressed genes (DEGs) at either time point. In mice, 14 and 1 DEGs were respectively identified after 7 and 90 days of exposure. Therefore, relaxed statistical criteria were employed to identify potential DEGs (pDEGs), followed by high-throughput benchmark dose modeling to identify responsive pDEGs for pathway enrichment analysis. This identified 288 and 168 pDEGs in the rat oral mucosa, of which only 20 and 7 showed evidence of dose-response. No significant pathway enrichment was obtained with either pDEG or dose-responsive pDEG lists. Similar results were obtained in mice. These analyses indicate a negligible transcriptional response in the oral mucosa of both species. Comparison of the total number of gene changes in the oral mucosa of rats and mice with responses in the duodenum of animals from the same study demonstrated remarkable dose-response concordance across tissues and species as a function of tissue chromium concentration. The low chromium levels in the oral mucosa and negligible transcript response are consistent with an absence of tissue lesions. These findings are used to compare the merits of linear and nonlinear approaches for deriving toxicity criteria based on the oral tumors in rats. Environ. Mol. Mutagen. 57:706-716, 2016. © 2016 The Authors. Environmental and Molecular Mutagenesis Published by Wiley Periodicals, Inc.

Detection of phase specificity of in vivo germ cell mutagens in an in vitro germ cell system

In vivo tests for male reproductive genotoxicity are time consuming, resource-intensive and their use should be minimised according to the principles of the 3Rs. Accordingly, we investigated the effects in vitro, of a variety of known, phase-specific germ cell mutagens, i.e., pre-meiotic, meiotic, and post-meiotic genotoxins, on rat spermatogenic cell types separated using Staput unit-gravity velocity sedimentation, evaluating DNA damage using the Comet assay. N-ethyl-N-nitrosourea (ENU), N-methyl-N-nitrosourea (MNU) (spermatogenic phase), 6-mercaptopurine (6-MP) and 5-bromo-2'-deoxy-uridine (5-BrdU) (meiotic phase), methyl methanesulphonate (MMS) and ethyl methanesulphonate (EMS) (post-meiotic phase) were selected for use as they are potent male rodent, germ cell mutagens in vivo. DNA damage was detected directly using the Comet assay and indirectly using the TUNEL assay. Treatment of the isolated cells with ENU and MNU produced the greatest concentration-related increase in DNA damage in spermatogonia. Spermatocytes were most sensitive to 6-MP and 5-BrdU while spermatids were particularly susceptible to MMS and EMS. Increases were found when measuring both Olive tail moment (OTM) and% tail DNA, but the greatest changes were in OTM. Parallel results were found with the TUNEL assay, which showed highly significant, concentration dependent effects of all these genotoxins on spermatogonia, spermatocytes and spermatids in the same way as for DNA damage. The specific effects of these chemicals on different germ cell types matches those produced in vivo. This approach therefore shows potential for use in the detection of male germ cell genotoxicity and could contribute to the reduction of the use of animals in such toxicity assays.

Nano-risk Science: application of toxicogenomics in an adverse outcome pathway framework for risk assessment of multi-walled carbon nanotubes

BACKGROUND

A diverse class of engineered nanomaterials (ENMs) exhibiting a wide array of physical-chemical properties that are associated with toxicological effects in experimental animals is in commercial use. However, an integrated framework for human health risk assessment (HHRA) of ENMs has yet to be established. Rodent 2-year cancer bioassays, clinical chemistry, and histopathological endpoints are still considered the 'gold standard' for detecting substance-induced toxicity in animal models. However, the use of data derived from alternative toxicological tools, such as genome-wide expression profiling and in vitro high-throughput assays, are gaining acceptance by the regulatory community for hazard identification and for understanding the underlying mode-of-action. Here, we conducted a case study to evaluate the application of global gene expression data in deriving pathway-based points of departure (PODs) for multi-walled carbon nanotube (MWCNT)-induced lung fibrosis, a non-cancer endpoint of regulatory importance.

METHODS

Gene expression profiles from the lungs of mice exposed to three individual MWCNTs with different physical-chemical properties were used within the framework of an adverse outcome pathway (AOP) for lung fibrosis to identify key biological events linking MWCNT exposure to lung fibrosis. Significantly perturbed pathways were categorized along the key events described in the AOP. Benchmark doses (BMDs) were calculated for each perturbed pathway and were used to derive transcriptional BMDs for each MWCNT.

RESULTS

Similar biological pathways were perturbed by the different MWCNT types across the doses and post-exposure time points studied. The pathway BMD values showed a time-dependent trend, with lower BMDs for pathways perturbed at the earlier post-exposure time points (24 h, 3d). The transcriptional BMDs were compared to the apical BMDs derived by the National Institute for Occupational Safety and Health (NIOSH) using alveolar septal thickness and fibrotic lesions endpoints. We found that regardless of the type of MWCNT, the BMD values for pathways associated with fibrosis were 14.0-30.4 μg/mouse, which are comparable to the BMDs derived by NIOSH for MWCNT-induced lung fibrotic lesions (21.0-27.1 μg/mouse).

CONCLUSIONS

The results demonstrate that transcriptomic data can be used to as an effective mechanism-based method to derive acceptable levels of exposure to nanomaterials in product development when epidemiological data are unavailable.

Application of in silico and in vitro methods in the development of adverse outcome pathway constructs in wildlife

There is a long history of using both in silico and in vitro methods to predict adverse effects in humans and environmental species where toxicity data are lacking. Currently, there is a great deal of interest in applying these methods to the development of so-called 'adverse outcome pathway' (AOP) constructs. The AOP approach provides a framework for organizing information at the chemical and biological level, allowing evidence from both in silico and in vitro studies to be rationally combined to fill gaps in knowledge concerning toxicological events. Fundamental to this new paradigm is a greater understanding of the mechanisms of toxicity and, in particular, where these mechanisms may be conserved across taxa, such as between model animals and related wild species. This presents an opportunity to make predictions across diverse species, where empirical data are unlikely to become available as is the case for most species of wildlife.

Development of a toxicogenomics signature for genotoxicity using a dose-optimization and informatics strategy in human cells

The development of in vitro molecular biomarkers to accurately predict toxicological effects has become a priority to advance testing strategies for human health risk assessment. The application of in vitro transcriptomic biomarkers promises increased throughput as well as a reduction in animal use. However, the existing protocols for predictive transcriptional signatures do not establish appropriate guidelines for dose selection or account for the fact that toxic agents may have pleiotropic effects. Therefore, comparison of transcriptome profiles across agents and studies has been difficult. Here we present a dataset of transcriptional profiles for TK6 cells exposed to a battery of well-characterized genotoxic and nongenotoxic chemicals. The experimental conditions applied a new dose optimization protocol that was based on evaluating expression changes in several well-characterized stress-response genes using quantitative real-time PCR in preliminary dose-finding studies. The subsequent microarray-based transcriptomic analyses at the optimized dose revealed responses to the test chemicals that were typically complex, often exhibiting substantial overlap in the transcriptional responses between a variety of the agents making analysis challenging. Using the nearest shrunken centroids method we identified a panel of 65 genes that could accurately classify toxicants as genotoxic or nongenotoxic. To validate the 65-gene panel as a genomic biomarker of genotoxicity, the gene expression profiles of an additional three well-characterized model agents were analyzed and a case study demonstrating the practical application of this genomic biomarker-based approach in risk assessment was performed to demonstrate its utility in genotoxicity risk assessment.

Approaches and perspectives to toxicogenetics and toxicogenomics

Toxicology is one of the scientific disciplines that has most<br />evolved in recent years due to scientific and technological<br />advances that have created a deeper understanding of the<br />genetic and molecular basis for appreciative variability in<br />toxic response from one person to another. The application<br />of this knowledge in toxicology is known as toxicogenetics<br />and toxicogenomics. The latter is the discipline that studies<br />the genomic response of organisms exposed to chemical<br />agents, including drugs, environmental pollutants, food<br />additives, and other commonly used chemical products.<br />The use of emerging omic technologies, such as genomics,<br />transcriptomics, proteomics, metabolomics and bioinformatics<br />techniques, permits the analysis of many variants of genes<br />simultaneously in an organism exposed to toxic agents in order<br />to search for genes susceptible to damage, to detect patterns<br />and mechanisms of toxicity, and determine specific profiles<br />of gene expression that give origin to biomarkers of exposure<br />and risk. This constitutes predictive toxicology.

Technical guide for applications of gene expression profiling in human health risk assessment of environmental chemicals

Toxicogenomics promises to be an important part of future human health risk assessment of environmental chemicals. The application of gene expression profiles (e.g., for hazard identification, chemical prioritization, chemical grouping, mode of action discovery, and quantitative analysis of response) is growing in the literature, but their use in formal risk assessment by regulatory agencies is relatively infrequent. Although additional validations for specific applications are required, gene expression data can be of immediate use for increasing confidence in chemical evaluations. We believe that a primary reason for the current lack of integration is the limited practical guidance available for risk assessment specialists with limited experience in genomics. The present manuscript provides basic information on gene expression profiling, along with guidance on evaluating the quality of genomic experiments and data, and interpretation of results presented in the form of heat maps, pathway analyses and other common approaches. Moreover, potential ways to integrate information from gene expression experiments into current risk assessment are presented using published studies as examples. The primary objective of this work is to facilitate integration of gene expression data into human health risk assessments of environmental chemicals.

A gene signature for gold nanoparticle-exposed human cell lines

A unique four-gene signature for AuNP exposure was identified using the cDNA microarray and evaluated by qPCR and biological assays in mammalian cell lines.

Proteomics for systems toxicology

Current toxicology studies frequently lack measurements at molecular resolution to enable a more mechanism-based and predictive toxicological assessment. Recently, a systems toxicology assessment framework has been proposed, which combines conventional toxicological assessment strategies with system-wide measurement methods and computational analysis approaches from the field of systems biology. Proteomic measurements are an integral component of this integrative strategy because protein alterations closely mirror biological effects, such as biological stress responses or global tissue alterations. Here, we provide an overview of the technical foundations and highlight select applications of proteomics for systems toxicology studies. With a focus on mass spectrometry-based proteomics, we summarize the experimental methods for quantitative proteomics and describe the computational approaches used to derive biological/mechanistic insights from these datasets. To illustrate how proteomics has been successfully employed to address mechanistic questions in toxicology, we summarized several case studies. Overall, we provide the technical and conceptual foundation for the integration of proteomic measurements in a more comprehensive systems toxicology assessment framework. We conclude that, owing to the critical importance of protein-level measurements and recent technological advances, proteomics will be an integral part of integrative systems toxicology approaches in the future.

Cancer biology, toxicology and alternative methods development go hand-in-hand

Toxicological research faces the challenge of integrating knowledge from diverse fields and novel technological developments generally in the biological and medical sciences. We discuss herein the fact that the multiple facets of cancer research, including discovery related to mechanisms, treatment and diagnosis, overlap many up and coming interest areas in toxicology, including the need for improved methods and analysis tools. Common to both disciplines, in vitro and in silico methods serve as alternative investigation routes to animal studies. Knowledge on cancer development helps in understanding the relevance of chemical toxicity studies in cell models, and many bioinformatics-based cancer biomarker discovery tools are also applicable to computational toxicology. Robotics-aided, cell-based, high-throughput screening, microscale immunostaining techniques and gene expression profiling analyses are common tools in cancer research, and when sequentially combined, form a tiered approach to structured safety evaluation of thousands of environmental agents, novel chemicals or engineered nanomaterials. Comprehensive tumour data collections in databases have been translated into clinically useful data, and this concept serves as template for computer-driven evaluation of toxicity data into meaningful results. Future 'cancer research-inspired knowledge management' of toxicological data will aid the translation of basic discovery results and chemicals- and materials-testing data to information relevant to human health and environmental safety.

An approach for integrating toxicogenomic data in risk assessment: the dibutyl phthalate case study

An approach for evaluating and integrating genomic data in chemical risk assessment was developed based on the lessons learned from performing a case study for the chemical dibutyl phthalate. A case study prototype approach was first developed in accordance with EPA guidance and recommendations of the scientific community. Dibutyl phthalate (DBP) was selected for the case study exercise. The scoping phase of the dibutyl phthalate case study was conducted by considering the available DBP genomic data, taken together with the entire data set, for whether they could inform various risk assessment aspects, such as toxicodynamics, toxicokinetics, and dose-response. A description of weighing the available dibutyl phthalate data set for utility in risk assessment provides an example for considering genomic data for future chemical assessments. As a result of conducting the scoping process, two questions--Do the DBP toxicogenomic data inform 1) the mechanisms or modes of action?, and 2) the interspecies differences in toxicodynamics?--were selected to focus the case study exercise. Principles of the general approach include considering the genomics data in conjunction with all other data to determine their ability to inform the various qualitative and/or quantitative aspects of risk assessment, and evaluating the relationship between the available genomic and toxicity outcome data with respect to study comparability and phenotypic anchoring. Based on experience from the DBP case study, recommendations and a general approach for integrating genomic data in chemical assessment were developed to advance the broader effort to utilize 21st century data in risk assessment.

Assessment of the mode of action underlying development of rodent small intestinal tumors following oral exposure to hexavalent chromium and relevance to humans

Abstract Chronic exposure to high concentrations of hexavalent chromium (Cr(VI)) in drinking water causes intestinal adenomas and carcinomas in mice, but not in rats. Cr(VI) causes damage to intestinal villi and crypt hyperplasia in mice after only one week of exposure. After two years of exposure, intestinal damage and crypt hyperplasia are evident in mice (but not rats), as are intestinal tumors. Although Cr(VI) has genotoxic properties, these findings suggest that intestinal tumors in mice arise as a result of chronic mucosal injury. To better understand the mode of action (MOA) of Cr(VI) in the intestine, a 90-day drinking water study was conducted to collect histological, biochemical, toxicogenomic and pharmacokinetic data in intestinal tissues. Using MOA analyses and human relevance frameworks proposed by national and international regulatory agencies, the weight of evidence supports a cytotoxic MOA with the following key events: (a) absorption of Cr(VI) from the intestinal lumen, (b) toxicity to intestinal villi, (c) crypt regenerative hyperplasia and (d) clonal expansion of mutations within the crypt stem cells, resulting in late onset tumorigenesis. This article summarizes the data supporting each key event in the MOA, as well as data that argue against a mutagenic MOA for Cr(VI)-induced intestinal tumors.

Profiling environmental chemicals for activity in the antioxidant response element signaling pathway using a high throughput screening approach

BACKGROUND

Oxidative stress has been implicated in the pathogenesis of a variety of diseases ranging from cancer to neurodegeneration, highlighting the need to identify chemicals that can induce this effect. The antioxidant response element (ARE) signaling pathway plays an important role in the amelioration of oxidative stress. Thus, assays that detect the up-regulation of this pathway could be useful for identifying chemicals that induce oxidative stress.

OBJECTIVES

We used cell-based reporter methods and informatics tools to efficiently screen a large collection of environmental chemicals and identify compounds that induce oxidative stress.

METHODS

We utilized two cell-based ARE assay reporters, β-lactamase and luciferase, to screen a U.S. National Toxicology Program 1,408-compound library (NTP 1408, which contains 1,340 unique compounds) for their ability to induce oxidative stress in HepG2 cells using quantitative high throughput screening (qHTS).

RESULTS

Roughly 3% (34 of 1,340) of the unique compounds demonstrated activity across both cell-based assays. Based on biological activity and structure-activity relationship profiles, we selected 50 compounds for retesting in the two ARE assays and in an additional follow-up assay that employed a mutated ARE linked to β-lactamase. Using this strategy, we identified 30 compounds that demonstrated activity in the ARE-bla and ARE-luc assays and were able to determine structural features conferring compound activity across assays.

CONCLUSIONS

Our results support the robustness of using two different cell-based approaches for identifying compounds that induce ARE signaling. Together, these methods are useful for prioritizing chemicals for further in-depth mechanism-based toxicity testing.

Review of current and "omics" methods for assessing the toxicity (genotoxicity, teratogenicity and nephrotoxicity) of herbal medicines and mushrooms

ETHNOPHARMACOLOGICAL RELEVANCE

The increasing use of traditional herbal medicines around the world requires more scientific evidence for their putative harmlessness. To this end, a plethora of methods exist, more or less satisfying. In this post-genome era, recent reviews are however scarce, not only on the use of new "omics" methods (transcriptomics, proteomics, metabonomics) for genotoxicity, teratogenicity, and nephrotoxicity assessment, but also on conventional ones.

METHODS

The present work aims (i) to review conventional methods used to assess genotoxicity, teratogenicity and nephrotoxicity of medicinal plants and mushrooms; (ii) to report recent progress in the use of "omics" technologies in this field; (iii) to underline advantages and limitations of promising methods; and lastly (iv) to suggest ways whereby the genotoxicity, teratogenicity, and nephrotoxicity assessment of traditional herbal medicines could be more predictive.

RESULTS

Literature and safety reports show that structural alerts, in silico and classical in vitro and in vivo predictive methods are often used. The current trend to develop "omics" technologies to assess genotoxicity, teratogenicity and nephrotoxicity is promising but most often relies on methods that are still not standardized and validated.

CONCLUSION

Hence, it is critical that toxicologists in industry, regulatory agencies and academic institutions develop a consensus, based on rigorous methods, about the reliability and interpretation of endpoints. It will also be important to regulate the integration of conventional methods for toxicity assessments with new "omics" technologies.

Comparison of hepatocarcinogen-induced gene expression profiles in conventional primary rat hepatocytes with in vivo rat liver

At present, substantial efforts are focused on the development of in vitro assays coupled with "omics" technologies for the identification of carcinogenic substances as an alternative to the classical 2-year rodent carcinogenicity bioassay. A prerequisite for the eventual regulatory acceptance of such assays, however, is the in vivo relevance of the observed in vitro findings. In the current study, hepatocarcinogen-induced gene expression profiles generated after the exposure of conventional cultures of primary rat hepatocytes to three non-genotoxic carcinogens (methapyrilene hydrochloride, piperonyl butoxide, and Wy-14643), three genotoxic carcinogens (aflatoxin B1, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, and 2-nitrofluorene), and two non-carcinogens (nifedipine and clonidine) are compared with previously obtained in vivo data after oral administration for up to 14 days of the same hepatocarcinogens to rats. In addition to the comparison of deregulated genes and functions per compound between in vivo and in vitro models, the major discriminating cellular pathways found in vivo in livers of exposed rats were examined for deregulation in vitro. Further, in vivo-derived gene signatures for the identification of genotoxic versus non-genotoxic carcinogens are used to classify in vitro-tested hepatocarcinogens and non-carcinogens. In the primary hepatocyte cultures, two out of the three tested genotoxic carcinogens mimicked the in vivo-relevant DNA damage response and were correctly assessed. Exposure to the non-genotoxic hepatocarcinogens, however, triggered a relatively weak response in the in vitro system, with no clear similarities to in vivo. This study contributes to the further optimization of toxicogenomics predictive tools when applied in in vitro settings.

Current understanding of the mechanism of benzene-induced leukemia in humans: implications for risk assessment

Benzene causes acute myeloid leukemia and probably other hematological malignancies. As benzene also causes hematotoxicity even in workers exposed to levels below the US permissible occupational exposure limit of 1 part per million, further assessment of the health risks associated with its exposure, particularly at low levels, is needed. Here, we describe the probable mechanism by which benzene induces leukemia involving the targeting of critical genes and pathways through the induction of genetic, chromosomal or epigenetic abnormalities and genomic instability, in a hematopoietic stem cell (HSC); stromal cell dysregulation; apoptosis of HSCs and stromal cells and altered proliferation and differentiation of HSCs. These effects modulated by benzene-induced oxidative stress, aryl hydrocarbon receptor dysregulation and reduced immunosurveillance, lead to the generation of leukemic stem cells and subsequent clonal evolution to leukemia. A mode of action (MOA) approach to the risk assessment of benzene was recently proposed. This approach is limited, however, by the challenges of defining a simple stochastic MOA of benzene-induced leukemogenesis and of identifying relevant and quantifiable parameters associated with potential key events. An alternative risk assessment approach is the application of toxicogenomics and systems biology in human populations, animals and in vitro models of the HSC stem cell niche, exposed to a range of levels of benzene. These approaches will inform our understanding of the mechanisms of benzene toxicity and identify additional biomarkers of exposure, early effect and susceptibility useful for risk assessment.

Use of genomic data in risk assessment case study: II. Evaluation of the dibutyl phthalate toxicogenomic data set

An evaluation of the toxicogenomic data set for dibutyl phthalate (DBP) and male reproductive developmental effects was performed as part of a larger case study to test an approach for incorporating genomic data in risk assessment. The DBP toxicogenomic data set is composed of nine in vivo studies from the published literature that exposed rats to DBP during gestation and evaluated gene expression changes in testes or Wolffian ducts of male fetuses. The exercise focused on qualitative evaluation, based on a lack of available dose-response data, of the DBP toxicogenomic data set to postulate modes and mechanisms of action for the male reproductive developmental outcomes, which occur in the lower dose range. A weight-of-evidence evaluation was performed on the eight DBP toxicogenomic studies of the rat testis at the gene and pathway levels. The results showed relatively strong evidence of DBP-induced downregulation of genes in the steroidogenesis pathway and lipid/sterol/cholesterol transport pathway as well as effects on immediate early gene/growth/differentiation, transcription, peroxisome proliferator-activated receptor signaling and apoptosis pathways in the testis. Since two established modes of action (MOAs), reduced fetal testicular testosterone production and Insl3 gene expression, explain some but not all of the testis effects observed in rats after in utero DBP exposure, other MOAs are likely to be operative. A reanalysis of one DBP microarray study identified additional pathways within cell signaling, metabolism, hormone, disease, and cell adhesion biological processes. These putative new pathways may be associated with DBP effects on the testes that are currently unexplained. This case study on DBP identified data gaps and research needs for the use of toxicogenomic data in risk assessment. Furthermore, this study demonstrated an approach for evaluating toxicogenomic data in human health risk assessment that could be applied to future chemicals.

Use of genomic data in risk assessment case study: I. Evaluation of the dibutyl phthalate male reproductive development toxicity data set

A case study was conducted, using dibutyl phthalate (DBP), to explore an approach to using toxicogenomic data in risk assessment. The toxicity and toxicogenomic data sets relative to DBP-related male reproductive developmental outcomes were considered conjointly to derive information about mode and mechanism of action. In this manuscript, we describe the case study evaluation of the toxicological database for DBP, focusing on identifying the full spectrum of male reproductive developmental effects. The data were assessed to 1) evaluate low dose and low incidence findings and 2) identify male reproductive toxicity endpoints without well-established modes of action (MOAs). These efforts led to the characterization of data gaps and research needs for the toxicity and toxicogenomic studies in a risk assessment context. Further, the identification of endpoints with unexplained MOAs in the toxicity data set was useful in the subsequent evaluation of the mechanistic information that the toxicogenomic data set evaluation could provide. The extensive analysis of the toxicology data set within the MOA context provided a resource of information for DBP in attempts to hypothesize MOAs (for endpoints without a well-established MOA) and to phenotypically anchor toxicogenomic and other mechanistic data both to toxicity endpoints and to available toxicogenomic data. This case study serves as an example of the steps that can be taken to develop a toxicological data source for a risk assessment, both in general and especially for risk assessments that include toxicogenomic data.