Predictive toxicogenomics approaches reveal underlying molecular mechanisms of nongenotoxic carcinogenicity

Characteristics of genomic signatures derived using univariate methods and mechanistically anchored functional descriptors for predicting drug- and xenobiotic-induced nephrotoxicity

ABSTRACT The ideal toxicity biomarker is composed of the properties of prediction (is detected prior to traditional pathological signs of injury), accuracy (high sensitivity and specificity), and mechanistic relationships to the endpoint measured (biological relevance). Gene expression-based toxicity biomarkers ("signatures") have shown good predictive power and accuracy, but are difficult to interpret biologically. We have compared different statistical methods of feature selection with knowledge-based approaches, using GeneGo's database of canonical pathway maps, to generate gene sets for the classification of renal tubule toxicity. The gene set selection algorithms include four univariate analyses: t-statistics, fold-change, B-statistics, and RankProd, and their combination and overlap for the identification of differentially expressed probes. Enrichment analysis following the results of the four univariate analyses, Hotelling T-square test, and, finally out-of-bag selection, a variant of cross-validation, were used to identify canonical pathway maps-sets of genes coordinately involved in key biological processes-with classification power. Differentially expressed genes identified by the different statistical univariate analyses all generated reasonably performing classifiers of tubule toxicity. Maps identified by enrichment analysis or Hotelling T-square had lower classification power, but highlighted perturbed lipid homeostasis as a common discriminator of nephrotoxic treatments. The out-of-bag method yielded the best functionally integrated classifier. The map "ephrins signaling" performed comparably to a classifier derived using sparse linear programming, a machine learning algorithm, and represents a signaling network specifically involved in renal tubule development and integrity. Such functional descriptors of toxicity promise to better integrate predictive toxicogenomics with mechanistic analysis, facilitating the interpretation and risk assessment of predictive genomic investigations.

Toxicogenomics discrimination of potential hepatocarcinogenicity of non-genotoxic compounds in rat liver

Long-term carcinogenicity testing of a compound is exceedingly time-consuming and costly, and requires many test animals, whereas the Ames test, which is based on the assumption that any substance that is mutagenic may also exert carcinogenic potential, is useful as a short-term screening assay but has major drawbacks. Although, in fact, 90% of compounds that give a positive Ames test cause cancer in laboratory animals, a good proportion of compounds that give a negative Ames test are also carcinogens; that is, there is no good correlation between carcinogenicity and negative Ames test results. As an alternative to these two approaches, we have tried applying toxicogenomics to predict the carcinogenicity of a compound from the gene expression profile induced in vivo. To establish our model, male Sprague-Dawley rats were orally administered test compounds (12 hepatocarcinogens and 26 non-hepatocarcinogens) for 28 days. Analysis of liver gene expression data by Support Vector Machines (SVM) dividing compounds into 'for training' and 'for test' (20 cases assigned randomly) allowed a set of marker genes to be tested for prediction of hepatocarcinogenicity. The developed prediction model was then validated with reference to the concordance rate with training data and test data, and a good performance was obtained. We will have new gene expression data and continue the validation of our model.

Toxicogenomics: the challenges and opportunities to identify biomarkers, signatures and thresholds to support mode-of-action

Toxicogenomics (TGx) can be defined as the application of "omics" techniques to toxicology and risk assessment. By identifying molecular changes associated with toxicity, TGx data might assist hazard identification and investigate causes. Early technical challenges were evaluated and addressed by consortia (e.g. ISLI/HESI and the Microarray Quality Control consortium), which demonstrated that TGx gave reliable and reproducible information. The MAQC also produced "best practice on signature generation" after conducting an extensive evaluation of different methods on common datasets. Two findings of note were the need for methods that control batch variability, and that the predictive ability of a signature changes in concert with the variability of the endpoint. The key challenge remaining is data interpretation, because TGx can identify molecular changes that are causal, associated with or incidental to toxicity. Application of Bradford Hill's tests for causation, which are used to build mode of action (MOA) arguments, can produce reasonable hypotheses linking altered pathways to phenotypic changes. However, challenges in interpretation still remain: are all pathway changes equal, which are most important and plausibly linked to toxicity? Therefore the expert judgement of the toxicologist is still needed. There are theoretical reasons why consistent alterations across a metabolic pathway are important, but similar changes in signalling pathways may not alter information flow. At the molecular level thresholds may be due to the inherent properties of the regulatory network, for example switch-like behaviours from some network motifs (e.g. positive feedback) in the perturbed pathway leading to the toxicity. The application of systems biology methods to TGx data can generate hypotheses that explain why a threshold response exists. However, are we adequately trained to make these judgments? There is a need for collaborative efforts between regulators, industry and academia to properly define how these technologies can be applied using appropriate case-studies.

Predictive modeling of chemical hazard by integrating numerical descriptors of chemical structures and short-term toxicity assay data

Quantitative structure-activity relationship (QSAR) models are widely used for in silico prediction of in vivo toxicity of drug candidates or environmental chemicals, adding value to candidate selection in drug development or in a search for less hazardous and more sustainable alternatives for chemicals in commerce. The development of traditional QSAR models is enabled by numerical descriptors representing the inherent chemical properties that can be easily defined for any number of molecules; however, traditional QSAR models often have limited predictive power due to the lack of data and complexity of in vivo endpoints. Although it has been indeed difficult to obtain experimentally derived toxicity data on a large number of chemicals in the past, the results of quantitative in vitro screening of thousands of environmental chemicals in hundreds of experimental systems are now available and continue to accumulate. In addition, publicly accessible toxicogenomics data collected on hundreds of chemicals provide another dimension of molecular information that is potentially useful for predictive toxicity modeling. These new characteristics of molecular bioactivity arising from short-term biological assays, i.e., in vitro screening and/or in vivo toxicogenomics data can now be exploited in combination with chemical structural information to generate hybrid QSAR-like quantitative models to predict human toxicity and carcinogenicity. Using several case studies, we illustrate the benefits of a hybrid modeling approach, namely improvements in the accuracy of models, enhanced interpretation of the most predictive features, and expanded applicability domain for wider chemical space coverage.

Metabolomics-on-a-chip and predictive systems toxicology in microfluidic bioartificial organs

The world faces complex challenges for chemical hazard assessment. Microfluidic bioartificial organs enable the spatial and temporal control of cell growth and biochemistry, critical for organ-specific metabolic functions and particularly relevant to testing the metabolic dose-response signatures associated with both pharmaceutical and environmental toxicity. Here we present an approach combining a microfluidic system with (1)H NMR-based metabolomic footprinting, as a high-throughput small-molecule screening approach. We characterized the toxicity of several molecules: ammonia (NH(3)), an environmental pollutant leading to metabolic acidosis and liver and kidney toxicity; dimethylsulfoxide (DMSO), a free radical-scavenging solvent; and N-acetyl-para-aminophenol (APAP, or paracetamol), a hepatotoxic analgesic drug. We report organ-specific NH(3) dose-dependent metabolic responses in several microfluidic bioartificial organs (liver, kidney, and cocultures), as well as predictive (99% accuracy for NH(3) and 94% for APAP) compound-specific signatures. Our integration of microtechnology, cell culture in microfluidic biochips, and metabolic profiling opens the development of so-called "metabolomics-on-a-chip" assays in pharmaceutical and environmental toxicology.

DNA methylation in the rat livers induced by low dosage isoniazid treatment

As a first-line anti-tuberculosis drug, isoniazid has a serious adverse side effect: hepatotoxicity. Therefore, the assessment and monitoring of hepatotoxicity from isoniazid to prevent liver injury are great concerns. In this experiment, we compared the levels of ALT in plasma and DNA methylation. 30 male SD rats were allocated randomly into two groups, a control group and an isoniazid group, and treated, respectively, with pure water and isoniazid at low dosage (10 mg/(kg day)) for 42 days by oral gavage. Five rats per group were sacrificed after 14, 28, and 42 days of isoniazid treatment. The levels of methylation in the genome and LINE-1 were measured, in which hypomethylation in the whole genome and LINE-1 repetitive sequences was observed in the INH group during the later period of the experiment (the 42nd day) accompanied with pathological changes in the liver. Thus, our results suggest that low dose of isoniazid can induce liver injury and that level of DNA methylation may be a more sensitive marker for monitoring drug-induced hepatotoxicity than aminotransferase.

Elevated progesterone receptor membrane component 1/sigma-2 receptor levels in lung tumors and plasma from lung cancer patients

Cancer is one of the leading causes of death, and there is an urgent need for new biomarkers and therapeutic targets. The progesterone receptor membrane component 1 (Pgrmc1) protein is upregulated in multiple types of cancer, and Pgrmc1 is required for tumor cell proliferation, motility and tumor formation in vivo. Furthermore, a small molecule inhibitor of Pgrmc1 suppressed the growth of lung, breast and cervical cancer cell lines. Recently, Pgrmc1 was identified as the sigma-2 receptor, a putative type of opioid receptor, and sigma-2 receptors are induced in cancers. However, Pgrmc1 shares no homology with known opioid or hormone receptors but is related to cytochrome b(5), and Pgrmc1 binds to heme and has reducing activity. In this study, we have analyzed Pgrmc1 levels in clinical tumor samples from squamous cell lung cancers (SCLC) and lung adenocarcinomas compared to corresponding nonmalignant tissue. Pgrmc1 levels increased significantly (p ≤ 0.05) in 12/15 SCLC samples and was elevated in poorly differentiated tumors. Pgrmc1 was highly expressed in SCLC cell lines, and SCLC cell survival was inhibited by siRNA knockdown of Pgrmc1 or the Pgrmc1 inhibitor AG-205. In adenocarcinomas, 6/15 tumors significantly had elevated Pgrmc1 levels, which correlated with patient survival. Pgrmc1 localizes to secretory vesicles in cancer cells, and Pgrmc1 was secreted by lung cancer cells. Furthermore, Pgrmc1 was significantly elevated in the plasma of lung cancer patients compared to noncancer patients. Together, the results demonstrate that Pgrmc1 is a potential tumor and serum biomarker, as well as a therapeutic target, for lung cancer.

Integrating transcriptomics and metabonomics to unravel modes-of-action of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) in HepG2 cells

BACKGROUND

The integration of different 'omics' technologies has already been shown in several in vivo studies to offer a complementary insight into cellular responses to toxic challenges. Being interested in developing in vitro cellular models as alternative to animal-based toxicity assays, we hypothesize that combining transcriptomics and metabonomics data improves the understanding of molecular mechanisms underlying the effects caused by a toxic compound also in vitro in human cells. To test this hypothesis, and with the focus on non-genotoxic carcinogenesis as an endpoint of toxicity, in the present study, the human hepatocarcinoma cell line HepG2 was exposed to the well-known environmental carcinogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD).

RESULTS

Transcriptomics as well as metabonomics analyses demonstrated changes in TCDD-exposed HepG2 in common metabolic processes, e.g. amino acid metabolism, of which some of the changes only being confirmed if both 'omics' were integrated. In particular, this integrated analysis identified unique pathway maps involved in receptor-mediated mechanisms, such as the G-protein coupled receptor protein (GPCR) signaling pathway maps, in which the significantly up-regulated gene son of sevenless 1 (SOS1) seems to play an important role. SOS1 is an activator of several members of the RAS superfamily, a group of small GTPases known for their role in carcinogenesis.

CONCLUSIONS

The results presented here were not only comparable with other in vitro studies but also with in vivo studies. Moreover, new insights on the molecular responses caused by TCDD exposure were gained by the cross-omics analysis.

Comparative analysis of predictive models for nongenotoxic hepatocarcinogenicity using both toxicogenomics and quantitative structure-activity relationships

The primary testing strategy to identify nongenotoxic carcinogens largely relies on the 2-year rodent bioassay, which is time-consuming and labor-intensive. There is an increasing effort to develop alternative approaches to prioritize the chemicals for, supplement, or even replace the cancer bioassay. In silico approaches based on quantitative structure-activity relationships (QSAR) are rapid and inexpensive and thus have been investigated for such purposes. A slightly more expensive approach based on short-term animal studies with toxicogenomics (TGx) represents another attractive option for this application. Thus, the primary questions are how much better predictive performance using short-term TGx models can be achieved compared to that of QSAR models, and what length of exposure is sufficient for high quality prediction based on TGx. In this study, we developed predictive models for rodent liver carcinogenicity using gene expression data generated from short-term animal models at different time points and QSAR. The study was focused on the prediction of nongenotoxic carcinogenicity since the genotoxic chemicals can be inexpensively removed from further development using various in vitro assays individually or in combination. We identified 62 chemicals whose hepatocarcinogenic potential was available from the National Center for Toxicological Research liver cancer database (NCTRlcdb). The gene expression profiles of liver tissue obtained from rats treated with these chemicals at different time points (1 day, 3 days, and 5 days) are available from the Gene Expression Omnibus (GEO) database. Both TGx and QSAR models were developed on the basis of the same set of chemicals using the same modeling approach, a nearest-centroid method with a minimum redundancy and maximum relevancy-based feature selection with performance assessed using compound-based 5-fold cross-validation. We found that the TGx models outperformed QSAR in every aspect of modeling. For example, the TGx models' predictive accuracy (0.77, 0.77, and 0.82 for the 1-day, 3-day, and 5-day models, respectively) was much higher for an independent validation set than that of a QSAR model (0.55). Permutation tests confirmed the statistical significance of the model's prediction performance. The study concluded that a short-term 5-day TGx animal model holds the potential to predict nongenotoxic hepatocarcinogenicity.

Genetic toxicology in the 21st century: reflections and future directions

A symposium at the 40th anniversary of the Environmental Mutagen Society, held from October 24-28, 2009 in St. Louis, MO, surveyed the current status and future directions of genetic toxicology. This article summarizes the presentations and provides a perspective on the future. An abbreviated history is presented, highlighting the current standard battery of genotoxicity assays and persistent challenges. Application of computational toxicology to safety testing within a regulatory setting is discussed as a means for reducing the need for animal testing and human clinical trials, and current approaches and applications of in silico genotoxicity screening approaches across the pharmaceutical industry were surveyed and are reported here. The expanded use of toxicogenomics to illuminate mechanisms and bridge genotoxicity and carcinogenicity, and new public efforts to use high-throughput screening technologies to address lack of toxicity evaluation for the backlog of thousands of industrial chemicals in the environment are detailed. The Tox21 project involves coordinated efforts of four U.S. Government regulatory/research entities to use new and innovative assays to characterize key steps in toxicity pathways, including genotoxic and nongenotoxic mechanisms for carcinogenesis. Progress to date, highlighting preliminary test results from the National Toxicology Program is summarized. Finally, an overview is presented of ToxCast™, a related research program of the U.S. Environmental Protection Agency, using a broad array of high throughput and high content technologies for toxicity profiling of environmental chemicals, and computational toxicology modeling. Progress and challenges, including the pressing need to incorporate metabolic activation capability, are summarized.

The embryonic stem cell test combined with toxicogenomics as an alternative testing model for the assessment of developmental toxicity

One of the most studied in vitro alternative testing methods for identification of developmental toxicity is the embryonic stem cell test (EST). Although the EST has been formally validated, the applicability domain as well as the predictability of the model needs further study to allow successful implementation of the EST as an alternative testing method in regulatory toxicity testing. Genomics technologies have already provided a proof of principle of their value in identification of toxicants such as carcinogenic compounds. Also within the EST, gene expression profiling has shown its value in the identification of developmental toxicity and in the evaluation of factors critical for risk assessment, such as dose and time responses. It is expected that the implementation of genomics into the EST will provide a more detailed end point evaluation as compared to the classical morphological scoring of differentiation cultures. Therefore, genomics may contribute to improvement of the EST, both in terms of definition of its applicability domain as well as its predictive capacity. In the present review, we present the progress that has been made with regard to the prediction of developmental toxicity using the EST combined with transcriptomics. Furthermore, we discuss the developments of additional aspects required for further optimization of the EST, including kinetics, the use of human embryonic stem cells (ESC) and computational toxicology. Finally, the current and future use of the EST model for prediction of developmental toxicity in testing strategies and in regulatory toxicity evaluations is discussed.

Alternative (non-animal) methods for cosmetics testing: current status and future prospects-2010

The 7th amendment to the EU Cosmetics Directive prohibits to put animal-tested cosmetics on the market in Europe after 2013. In that context, the European Commission invited stakeholder bodies (industry, non-governmental organisations, EU Member States, and the Commission's Scientific Committee on Consumer Safety) to identify scientific experts in five toxicological areas, i.e. toxicokinetics, repeated dose toxicity, carcinogenicity, skin sensitisation, and reproductive toxicity for which the Directive foresees that the 2013 deadline could be further extended in case alternative and validated methods would not be available in time. The selected experts were asked to analyse the status and prospects of alternative methods and to provide a scientifically sound estimate of the time necessary to achieve full replacement of animal testing. In summary, the experts confirmed that it will take at least another 7-9 years for the replacement of the current in vivo animal tests used for the safety assessment of cosmetic ingredients for skin sensitisation. However, the experts were also of the opinion that alternative methods may be able to give hazard information, i.e. to differentiate between sensitisers and non-sensitisers, ahead of 2017. This would, however, not provide the complete picture of what is a safe exposure because the relative potency of a sensitiser would not be known. For toxicokinetics, the timeframe was 5-7 years to develop the models still lacking to predict lung absorption and renal/biliary excretion, and even longer to integrate the methods to fully replace the animal toxicokinetic models. For the systemic toxicological endpoints of repeated dose toxicity, carcinogenicity and reproductive toxicity, the time horizon for full replacement could not be estimated.

Concentration-dependent gene expression responses to flusilazole in embryonic stem cell differentiation cultures

The murine embryonic stem cell test (EST) is designed to evaluate developmental toxicity based on compound-induced inhibition of embryonic stem cell (ESC) differentiation into cardiomyocytes. The addition of transcriptomic evaluation within the EST may result in enhanced predictability and improved characterization of the applicability domain, therefore improving usage of the EST for regulatory testing strategies. Transcriptomic analyses assessing factors critical for risk assessment (i.e. dose) are needed to determine the value of transcriptomic evaluation in the EST. Here, using the developmentally toxic compound, flusilazole, we investigated the effect of compound concentration on gene expression regulation and toxicity prediction in ESC differentiation cultures. Cultures were exposed for 24 h to multiple concentrations of flusilazole (0.54-54 μM) and RNA was isolated. In addition, we sampled control cultures 0, 24, and 48 h to evaluate the transcriptomic status of the cultures across differentiation. Transcriptomic profiling identified a higher sensitivity of development-related processes as compared to cell division-related processes in flusilazole-exposed differentiation cultures. Furthermore, the sterol synthesis-related mode of action of flusilazole toxicity was detected. Principal component analysis using gene sets related to normal ESC differentiation was used to describe the dynamics of ESC differentiation, defined as the 'differentiation track'. The concentration-dependent effects on development were reflected in the significance of deviation of flusilazole-exposed cultures from this transcriptomic-based differentiation track. Thus, the detection of developmental toxicity in EST using transcriptomics was shown to be compound concentration-dependent. This study provides further insight into the possible application of transcriptomics in the EST as an improved alternative model system for developmental toxicity testing.

Chemical carcinogenesis

Understanding the relationship of chemicals to carcinogenesis has progressed significantly since the initial observations of Hill and Pott in the 1700's. Distinguishing between DNA-reactive chemicals and those which increase cancer risk by increasing cell proliferation has been a major breakthrough in delineating overall mechanisms. Competing processes for activation versus inactivation of chemicals occur at many levels, including metabolism, DNA repair, and cellular repair processes. These processes can be affected by other agents to decrease carcinogenesis (chemoprevention). Increasing knowledge of the multiple steps of carcinogenesis is leading to improved methods for screening chemicals for carcinogenic activity and for regulatory decision making. Improvements in assessment of modes of action involved in animal and in vitro models have led to more rational approaches to assessing relevance to humans. The advent of genomics and high-throughput technologies have contributed to investigations of mechanisms and is beginning to impact development of better methods for screening chemicals. Based on developments in basic research, epidemiology, and astute clinical observations, the major risk factors and etiologic agents have been identified for a majority of cancers, which is beginning to lead to methods to decrease cancer incidence overall and identify targets for early detection and treatment.

Progesterone receptor membrane component 1 inhibits the activity of drug-metabolizing cytochromes P450 and binds to cytochrome P450 reductase

Progesterone receptor membrane component 1 (PGRMC1) has been shown to interact with several cytochromes P450 (P450s) and to activate enzymatic activity of P450s involved in sterol biosynthesis. We analyzed the interactions of PGRMC1 with the drug-metabolizing P450s, CYP2C2, CYP2C8, and CYP3A4, in transfected cells. Based on coimmunoprecipitation assays, PGRMC1 bound efficiently to all three P450s, and binding to the catalytic cytoplasmic domain of CYP2C2 was much more efficient than to a chimera containing only the N-terminal transmembrane domain. Down-regulation of PGRMC1 expression levels in human embryonic kidney 293 and HepG2 cell lines stably expressing PGRMC1-specific small interfering RNA had no effect on the endoplasmic reticulum localization and expression levels of P450s, whereas enzymatic activities of CYP2C2, CYP2C8, and CYP3A4 were slightly higher in PGRMC1-deficient cells. Cotransfection of cells with P450s and PGRMC1 resulted in PGRMC1 concentration-dependent inhibition of the P450 activities, and this inhibition was partially reversed by increased expression of the P450 reductase (CPR). In contrast, CYP51 activity was decreased by down-regulation of PGRMC1 and expression of PGRMC1 in the PGRMC1-deficient cells increased CYP51 activity. In cells cotransfected with CPR and PGRMC1, strong binding of CPR to PGRMC1 was observed; however, in the presence of CYP2C2, interaction of PGRMC1 with CPR was significantly reduced, suggesting that CYP2C2 competes with CPR for binding to PGRMC1. These data show that in contrast to sterol synthesizing P450, PGRMC1 is not required for the activities of several drug-metabolizing P450s, and its overexpression inhibits those P450 activities. Furthermore, PGRMC1 binds to CPR, which may influence P450 activity.

The enhanced value of combining conventional and "omics" analyses in early assessment of drug-induced hepatobiliary injury

The InnoMed PredTox consortium was formed to evaluate whether conventional preclinical safety assessment can be significantly enhanced by incorporation of molecular profiling ("omics") technologies. In short-term toxicological studies in rats, transcriptomics, proteomics and metabolomics data were collected and analyzed in relation to routine clinical chemistry and histopathology. Four of the sixteen hepato- and/or nephrotoxicants given to rats for 1, 3, or 14days at two dose levels induced similar histopathological effects. These were characterized by bile duct necrosis and hyperplasia and/or increased bilirubin and cholestasis, in addition to hepatocyte necrosis and regeneration, hepatocyte hypertrophy, and hepatic inflammation. Combined analysis of liver transcriptomics data from these studies revealed common gene expression changes which allowed the development of a potential sequence of events on a mechanistic level in accordance with classical endpoint observations. This included genes implicated in early stress responses, regenerative processes, inflammation with inflammatory cell immigration, fibrotic processes, and cholestasis encompassing deregulation of certain membrane transporters. Furthermore, a preliminary classification analysis using transcriptomics data suggested that prediction of cholestasis may be possible based on gene expression changes seen at earlier time-points. Targeted bile acid analysis, based on LC-MS metabonomics data demonstrating increased levels of conjugated or unconjugated bile acids in response to individual compounds, did not provide earlier detection of toxicity as compared to conventional parameters, but may allow distinction of different types of hepatobiliary toxicity. Overall, liver transcriptomics data delivered mechanistic and molecular details in addition to the classical endpoint observations which were further enhanced by targeted bile acid analysis using LC/MS metabonomics.

Pgrmc1 (progesterone receptor membrane component 1) associates with epidermal growth factor receptor and regulates erlotinib sensitivity

Tumorigenesis requires the concerted action of multiple pathways, including pathways that stimulate proliferation and metabolism. Epidermal growth factor receptor (EGFR) is a transmembrane receptor-tyrosine kinase that is associated with cancer progression, and the EGFR inhibitors erlotinib/tarceva and tyrphostin/AG-1478 are potent anti-cancer therapeutics. Pgrmc1 (progesterone receptor membrane component 1) is a cytochrome b(5)-related protein that is up-regulated in tumors and promotes cancer growth. Pgrmc1 and its homologues have been implicated in cell signaling, and we show here that Pgrmc1 increases susceptibility to AG-1478 and erlotinib, increases plasma membrane EGFR levels, and co-precipitates with EGFR. Pgrmc1 co-localizes with EGFR in cytoplasmic vesicles and co-fractionates with EGFR in high density microsomes. The findings have therapeutic potential because a Pgrmc1 small molecule ligand, which inhibits growth in a variety of cancer cell types, de-stabilized EGFR in multiple tumor cell lines. EGFR is one of the most potent receptor-tyrosine kinases driving tumorigenesis, and our data support a role for Pgrmc1 in promoting several cancer phenotypes at least in part by binding EGFR and stabilizing plasma membrane pools of the receptor.

An elastic network model to identify characteristic stress response genes

Exposing eukaryotic cells to a toxic compound and subsequent gene expression profiling may allow the prediction of selected toxic effects based on changes in gene expression. This objective is complicated by the observation that compounds with different modes of toxicity cause similar changes in gene expression and that a global stress response affects many genes. We developed an elastic network model of global stress response with nodes representing genes which are connected by edges of graded coexpression. The expression of only few genes have to be known to model the global stress response of all but a few atypical responder genes. Those required genes and the atypical response genes are shown to be good biomarker for tox predictions. In total, 138 experiments and 13 different compounds were used to train models for different toxicity classes. The deduced biomarkers were shown to be biologically plausible. A neural network was trained to predict the toxic effects of compounds from profiling experiments. On a validation data set of 189 experiments with 16 different compounds the accuracy of the predictions was assessed: 14 out of 16 compounds have been classified correctly. Derivation of model based biomarkers through the elastic network approach can naturally be extended to other areas beyond toxicology since subtle signals against a broad response background are common in biological studies.

Alterations in mRNA expression of steroid receptors and heat shock proteins in the liver of rare minnow (Grobiocypris rarus) exposed to atrazine and p,p'-DDE

The chaperon role of heat shock proteins (HSPs) throughout the life cycle of steroid receptors have been demonstrated in vitro. However, the actions of HSPs in steroid pathways in animals especially in fish were unclear. In this study, sexually mature rare minnows (Gobiocypris rarus) were exposed to typical endocrine disruptors (atrazine and p,p'-DDE). Hypertrophy of hepatocytes at the 333 microg/l atrazine treatment and atrophy of hepatocytes in all p,p'-DDE treatments were observed. The expression of liver hsp70 and hsp90 in atrazine treatments were significantly up-regulated. Moreover, remarkable increases in the expression of androgen receptor (ar) and estrogen receptor (er) were observed, while alterations of the glucoticorcoid receptor (gr) expression was not significant in atrazine exposed fish. The expression of ar, er, gr, hsp70 and hsp90 were significantly suppressed following p,p'-DDE exposure. These results demonstrate that the expression of hsp70 and hsp90 is altered along with changes of steroid receptors in vivo. Therefore, the results are consistent with the possibility that in fish heat shock proteins (HSPs) play chaperon roles for the steroid receptors in vivo, which also concurs with previous in vitro mammalian studies.

The Japanese toxicogenomics project: application of toxicogenomics

Biotechnology advances have provided novel methods for the risk assessment of chemicals. The application of microarray technologies to toxicology, known as toxicogenomics, is becoming an accepted approach for identifying chemicals with potential safety problems. Gene expression profiling is expected to identify the mechanisms that underlie the potential toxicity of chemicals. This technology has also been applied to identify biomarkers of toxicity to predict potential hazardous chemicals. Ultimately, toxicogenomics is expected to aid in risk assessment. The following discussion explores potential applications and features of the Japanese Toxicogenomics Project.

Liverbeads: a practical and relevant in vitro model for gene induction investigations

Cryopreserved rat hepatocytes entrapped within an alginate matrix, commercially available as Liverbeads, were evaluated for their relevance as a screening tool for gene induction in vitro, using quantitative real-time reverse transcriptase-polymerase chain reaction. They were treated with the reference compounds beta-naphthoflavone (BNF), phenobarbital (PB), pregnenolone 16alpha-carbonitrile (PCN), and clofibric acid (CLO) and analyzed for mRNA levels of Cyp1a1, Cyp2b1, Cyp3a1, Cyp4a1, Ugt1a6, and Ugt2b1. In addition, for PB and PCN, the results were compared with those obtained in rat liver in vivo. For each inducer, the gene induction profiles obtained with the Liverbeads in vitro model were time- and dose-dependent. The in vitro gene expression profiles confirmed the corresponding known P450 and UGT induction by each reference compound. In particular, the most strongly induced genes were Cyp1a1 by BNF, Cyp2b1 by PB, Cyp3a1 and Ugt2b1 by PCN, and Cyp4a1 and Cyp2b1 by CLO. Other genes investigated were also induced by the reference compounds, but the expression levels were lower, and increases were seen only after prolonged treatment. In particular, Ugt1a6 and Cyp2b1 were increased by BNF, Cyp1a1, Cyp3a1, and Ugt2b1 by PB, and Cyp3a1 and Ugt2b1 by CLO. All of these results correlated well with published in vitro data and our in vivo data. In conclusion, our results suggest that Liverbeads is a relevant and useful in vitro screening tool for determining gene induction profiles of new molecules. In addition, because Liverbeads from different species are available, this tool offers the possibility to conduct interspecies comparisons.

Ginkgo biloba extract induces gene expression changes in xenobiotics metabolism and the Myc-centered network

The use of herbal dietary supplements in the United States is rapidly growing, and it is crucial that the quality and safety of these preparations be ensured. To date, it is still a challenge to determine the mechanisms of toxicity induced by mixtures containing many chemical components, such as herbal dietary supplements. We previously proposed that analyses of the gene expression profiles using microarrays in the livers of rodents treated with herbal dietary supplements is a potentially practical approach for understanding the mechanism of toxicity. In this study, we utilized microarrays to analyze gene expression changes in the livers of male B6C3F1 mice administered Ginkgo biloba leaf extract (GBE) by gavage for 2 years, and to determine pathways and mechanisms associated with GBE treatments. Analysis of 31,802 genes revealed that there were 129, 289, and 2,011 genes significantly changed in the 200, 600, and 2,000 mg/kg treatment groups, respectively, when compared with control animals. Drug metabolizing genes were significantly altered in response to GBE treatments. Pathway and network analyses were applied to investigate the gene relationships, functional clustering, and mechanisms involved in GBE exposure. These analyses indicate alteration in the expression of genes coding for drug metabolizing enzymes, the NRF2-mediated oxidative stress response pathway, and the Myc gene-centered network named "cell cycle, cellular movement, and cancer" were found. These results indicate that Ginkgo biloba-related drug metabolizing enzymes may cause herb-drug interactions and contribute to hepatotoxicity. In addition, the outcomes of pathway and network analysis may be used to elucidate the toxic mechanisms of Ginkgo biloba.

State-of-the-art genomics approaches in toxicology

Genomics may be an effective tool in decreasing the lengthy drug development process and reducing compound attrition. It can generate specific gene expression profiles induced by chemicals that can be linked to dose and response. Toxicogenomics can identify sensitive biomarkers of early deleterious effects, distinguish genotoxic from non-genotoxic carcinogens and can provide information on the mechanism of action. It can help bridge in vitro to in vivo findings and provide context for preclinical data and thus address human health risks. Issues and shortcomings that still need to be resolved or improved for efficient incorporation of genomics in drug development and environmental toxicology research include data analysis, data interpretation tools and accessible data repositories. In addition, implementation of toxicogenomics in early screening or drug discovery phases and effective use of this information by project teams remains a challenge.

Atrazine affects kidney and adrenal hormones (AHs) related genes expressions of rare minnow (Gobiocypris rarus)

Atrazine, one of the most widely used herbicides, has been proved to interfere with sexual hormones. However few studies have considered the effects of atrazine on adrenal hormones (AH). In this study, rare minnow (Gobiocypris rarus) was exposed to 0, 3, 10, 33, 100 and 333microg/l atrazine for 28 days. The histopathology of kidney and gill was examined and the expressions of AHs-related genes including Na(+),K(+)-ATPase, glucocorticoid receptor (gr), heat shock protein 70 (hsp70), and heat shock protein 90 (hsp90) in kidney and gill were quantitatively determined. Histopathological observation revealed obvious lesions in gill including hyperplasia, necrosis in epithelium region, aneurysm and lamellar fusion at concentrations as low as 10microg/l. The observed lesions in kidney included extensive expansion in the lumen, degenerative and necrotic changes of the tubular epithelia, shrinkage of the glomerulus as well as increase of the Bowman's space at concentrations as low as 10microg/l. The expressions of Na(+),K(+)-ATPase, gr, hsp70 and hsp90 in the kidney of females were significantly decreased at all concentrations. For males, the expressions of hsp90 in the kidney of all treated groups were significantly down-regulated, while gr at all concentrations and hsp70 at 10, 33, 100microg/l were significantly up-regulated. However in the gill, the expressions of these genes were not significantly different from the control. These results indicated that exposure to atrazine caused impairments of kidney and gill of fish at environmental related concentrations. Histopathological lesions could partly attribute to the changes of the expressions of AHs-related genes in kidney. We concluded also that atrazine is a potential AHs-disruptor and AHs-related genes in kidney of fish could be used as sensitive molecular biomarkers.

Modelling the genesis and treatment of cancer: the potential role of physiologically based pharmacodynamics

Physiologically based modelling of pharmacodynamics/toxicodynamics requires an a priori knowledge on the underlying mechanisms causing toxicity or causing the disease. In the context of cancer, the objective of the expert meeting was to discuss the molecular understanding of the disease, modelling approaches used so far to describe the process, preclinical models of cancer treatment and to evaluate modelling approaches developed based on improved knowledge. Molecular events in cancerogenesis can be detected using 'omics' technology, a tool applied in experimental carcinogenesis, but also for diagnostics and prognosis. The molecular understanding forms the basis for new drugs, for example targeting protein kinases specifically expressed in cancer. At present, empirical preclinical models of tumour growth are in great use as the development of physiological models is cost and resource intensive. Although a major challenge in PKPD modelling in oncology patients is the complexity of the system, based in part on preclinical models, successful models have been constructed describing the mechanism of action and providing a tool to establish levels of biomarker associated with efficacy and assisting in defining biologically effective dose range selection for first dose in man. To follow the concentration in the tumour compartment enables to link kinetics and dynamics. In order to obtain a reliable model of tumour growth dynamics and drug effects, specific aspects of the modelling of the concentration-effect relationship in cancer treatment that need to be accounted for include: the physiological/circadian rhythms of the cell cycle; the treatment with combinations and the need to optimally choose appropriate combinations of the multiple agents to study; and the schedule dependence of the response in the clinical situation.

Evaluation of Possible Carcinogenic Risk to Humans Based on Liver Tumors in Rodent Assays

The two-year rodent bioassay remains the mainstay for carcinogenicity testing, although numerous difficulties have been identified. Fundamentally, a chemical can increase the risk of cancer (1) by damaging DNA directly (DNA reactive) or (2) indirectly by increasing the number of DNA replications (non–DNA reactive). Mechanistic research has identified numerous precursor lesions in the sequence of key events necessary for neoplasia development. Based on these concepts, the author has proposed a short-term (thirteen-week) assay for screening for carcinogenic potential based on a mode of action analysis and on readily available, identifiable preneoplastic changes. A screening assay that detects all potential rodent hepatocarcinogens has been previously identified ( Toxicol Pathol32 [2004], 393–401) including increased liver weight, hepatocellular necrosis, hypertrophy, and cytomegaly. Labeling index for DNA replication might supply additional support. These markers have high sensitivity but low specificity. However, most chemicals can be appropriately classified as to their mode(s) of action for hepatocarcinogenesis with follow-up mechanistic studies, and a rational evaluation of their relevance to humans can be made. A similar process can be envisioned for other tissues for evaluation for carcinogenic potential.

Progesterone receptor membrane component 1 (Pgrmc1): a heme-1 domain protein that promotes tumorigenesis and is inhibited by a small molecule

Tumorigenesis requires the concerted action of multiple pathways, including pathways that stimulate proliferation and increase metabolism. Progesterone receptor membrane component 1 (Pgrmc1) is related to cytochrome b5, binds to heme, and is associated with DNA damage resistance and apoptotic suppression. Pgrmc1 is induced by carcinogens, including dioxin, and is up-regulated in multiple types of cancer. In the present study, we found that Pgrmc1 increased in vivo tumor growth, anchorage-independent growth, and migration. Pgrmc1 also promoted proliferation in the absence of serum in A549 non-small cell lung cancer cells but enhanced proliferation regardless of serum concentration in MDA-MB-468 breast cancer cells. Pgrmc1 promotes cholesterol synthesis and binds to Insig (insulin-induced gene), Scap (sterol regulatory element binding protein cleavage activating protein), and P450 proteins, but Pgrmc1 did not affect cholesterol synthesis in lung cancer cells. Pgrmc1 is also associated with progesterone signaling and plasminogen activator inhibitor (PAI1) RNA binding protein, but neither progesterone activity nor PAI1 transcript levels were altered in Pgrmc1-knockdown lung cancer cells. Pgrmc1 homologues bind to aryl ligands identified in an in silico screen, and we have found that a Pgrmc1 ligand induced cell death in a Pgrmc1-specific manner in multiple breast and lung tumor cell lines. Our data support a role for Pgrmc1 in promoting cancer-associated phenotypes and provide a therapeutic approach for targeting Pgrmc1 with a small molecule in lung and breast cancer.

Discrimination of carcinogens by hepatic transcript profiling in rats following 28-day administration

This study aimed at discriminating carcinogens on the basis of hepatic transcript profiling in the rats administrated with a variety of carcinogens and non-carcinogens. We conducted 28-day toxicity tests in male F344 rats with 47 carcinogens and 26 non-carcinogens, and then investigated periodically the hepatic gene expression profiles using custom microarrays. By hierarchical cluster analysis based on significantly altered genes, carcinogens were clustered into three major groups (Group 1 to 3). The formation of these groups was not affected by the gene sets used as well as the administration period, indicating that the grouping of carcinogens was universal independent of the conditions of both statistical analysis and toxicity testing. Seventeen carcinogens belonging to Group 1 were composed of mainly rat hepatocarcinogens, most of them being mutagenic ones. Group 2 was formed by three subgroups, which were composed of 23 carcinogens exhibiting distinct properties in terms of genotoxicity and target tissues, namely nonmutagenic hepatocarcinogens, and mutagenic and nonmutagenic carcinogens both of which are targeted to other tissues. Group 3 contained 6 carcinogens including 4 estrogenic substances, implying the group of estrogenic carcinogens. Gene network analyses revealed that the significantly altered genes in Group 1 included Bax, Tnfrsf6, Btg2, Mgmt and Abcb1b, suggesting that p53-mediated signaling pathway involved in early pathologic alterations associated with preceding mutagenic carcinogenesis. Thus, the common transcriptional signatures for each group might reflect the early molecular events of carcinogenesis and hence would enable us to identify the biomarker genes, and then to develop a new assay for carcinogenesis prediction.

Cancer risk assessment approaches at the FDA/CDER: Is the era of the 2-year bioassay drawing to a close?

Determining the carcinogenic potential of materials to which humans have significant exposures is an important, complex, and imperfect exercise. Not only are the methods for such determinations protracted and expensive and use large numbers of animals, extrapolation of data from such studies to human risk is imprecise. With improved understanding of oncogene activation and tumor suppressor gene inactivation, a number of animal models have been developed to dramatically reduce latency for chemically induced cancers and has led to the development and use of shorter carcinogenicity assays. Recent studies by a number of investigators suggest that specific gene signature patterns seen after short-term exposure of rats to test chemicals can predict long-term outcomes in cancer bioassays with relatively high accuracy. In addition, a recent survey performed by PhRMA member companies examined two hundred drug years to determine whether histological biomarkers seen at the end of a six- or twelve-month toxicology study in rats can predict the outcome of a two-year carcinogenicity study. With only a handful of exceptions, chronic studies appear capable of predicting effects at the end of two years with good accuracy. It is hoped that the combination of results from transgenic mouse assays and six-month rat studies will soon supplant the need for most two-year bioassays.

Effects of exposure to acetochlor on the expression of thyroid hormone related genes in larval and adult rare minnow (Gobiocypris rarus)

Acetochlor is a commonly used herbicide in agricultural environments, which can accelerate T(3)-induced metamorphosis in amphibians. However, effects of acetochlor on the thyroid system in fish are still unclear at environmentally relevant concentrations. In this study, expression profiles of thyroid hormone receptor alpha (tralpha), deiodinase (d1 and d2), malic enzyme (me), and sodium iodide symporter (nis) genes were determined in larval and adult rare minnow (Gobiocypris rarus) after exposure to different levels of acetochlor (20, 200, and 2000ng/l) for 21 days, respectively. Furthermore, plasma thyroid hormones (THs) levels in adults were evaluated simultaneously. The results showed that d1, d2, me, and nis mRNA levels were significantly down-regulated in the larvae. Similar expression changes of these genes in female brains were found. However, in adult livers, these gene expressions had no significantly effects, except the tralpha mRNA level significantly up-regulated. These results indicated that exposure to acetochlor could result in tissue-specific alternative expression of TH-related genes in adults. Moreover, the expression of d2 and me showed a positive correlation with plasma T(4) levels in female brains. Therefore, larval development and adult brain of rare minnow could be affected by acetochlor at environmentally relevant concentrations.

Erratum to "Prediction of non-genotoxic carcinogenesis in rats using changes in gene expression following acute dosing"

Carcinogenicity of chemicals can currently only be evaluated in 2-year rodent bioassays. Therefore, the development of early biomarkers for carcinogenesis would result in substantial savings in time and expense. The current study investigates whether early changes in gene expression may be developed as markers for cancer. Animals were treated for 1 or 5 days with either non-genotoxic carcinogens or non-carcinogens and gene expression was analyzed by quantitative PCR (qPCR).We tested two gene signatures previously reported to detect non-genotoxic carcinogens. Using one gene signature it was confirmed that 3/3 nongenotoxic carcinogens and 2/2 non-carcinogens are correctly identified with data from 1 or 5 days of dosing. In contrast an alternative signature correctly identified 0/3 and 2/3 nongenotoxic carcinogens at 1 and 5 days of treatment, respectively and 2/2 non-carcinogens at both time-points. Additionally, we evaluated a novel panel of putative biomarker genes, from the literature, many of which have roles in cell growth and division, including myc, cdc2 and mcm6. These genes were significantly induced by non-genotoxic carcinogens and not by non-carcinogens. Using the average fold-induction across this panel, 2/3 non-genotoxic carcinogens were detected on both day 1 and day 5. These data support the idea that acute changes in gene expression may provide biomarkers for non-genotoxic carcinogenesis but also highlight interesting differences in the sensitivities of distinct gene signatures.

Global liver proteomics of rats exposed for 5 days to phenobarbital identifies changes associated with cancer and with CYP metabolism

A global proteomics approach was applied to model the hepatic response elicited by the toxicologically well-characterized xenobiotic phenobarbital (PB), a prototypical inducer of hepatic xenobiotic metabolizing enzymes and a well-known nongenotoxic liver carcinogen in rats. Differential detergent fractionation two-dimensional liquid chromatography electrospray ionization tandem mass spectrometry and systems biology modeling were used to identify alterations in toxicologically relevant hepatic molecular functions and biological processes in the livers of rats following a 5-day exposure to PB at 80 mg/kg/day or a vehicle control. Of the 3342 proteins identified, expression of 121 (3.6% of the total proteins) was significantly increased and 127 (3.8%) significantly decreased in the PB group compared to controls. The greatest increase was seen for cytochrome P450 (CYP) 2B2 (167-fold). All proteins with statistically significant differences from control were then analyzed using both Gene Ontology (GO) and Ingenuity Pathways Analysis (IPA, 5.0 IPA-Tox) for cellular location, function, network connectivity, and possible disease processes, especially as they relate to CYP-mediated metabolism and nongenotoxic carcinogenesis mechanisms. The GO results suggested that PB's mechanism of nongenotoxic carcinogenesis involves both increased xenobiotic metabolism, especially induction of the 2B subfamily of CYP enzymes, and increased cell cycle activity. Apoptosis, however, also increased, perhaps, as an attempt to counter the rising cancer threat. Of the IPA-mapped proteins, 41 have functions which are procarcinogenic and 14 anticarcinogenic according to the hypothesized nongenotoxic mechanism of imbalance between apoptosis and cellular proliferation. Twenty-two additional IPA nodes can be classified as procarcinogenic by the competing theory of increased metabolism resulting in the formation of reactive oxygen species. Since the systems biology modeling corresponded well to PB effects previously elucidated via more traditional methods, the global proteomic approach is proposed as a new screening methodology that can be incorporated into future toxicological studies.

Dose- and time-dependent effects of phenobarbital on gene expression profiling in human hepatoma HepaRG cells

Phenobarbital (PB) induces or represses a wide spectrum of genes in rodent liver. Much less is known about its effects in human liver. We used pangenomic cDNA microarrays to analyze concentration- and time-dependent gene expression profile changes induced by PB in the well-differentiated human HepaRG cell line. Changes in gene expression profiles clustered at specific concentration ranges and treatment times. The number of correctly annotated genes significantly modulated by at least three different PB concentration ranges (spanning 0.5 to 3.2 mM) at 20 h exposure amounted to 77 and 128 genes (p< or =0.01) at 2- and 1.8-fold filter changes, respectively. At low concentrations (0.5 and 1 mM), PB-responsive genes included the well-recognized CAR- and PXR-dependent responsive cytochromes P450 (CYP2B6, CYP3A4), sulfotransferase 2A1 and plasma transporters (ABCB1, ABCC2), as well as a number of genes critically involved in various metabolic pathways, including lipid (CYP4A11, CYP4F3), vitamin D (CYP24A1) and bile (CYP7A1 and CYP8B1) metabolism. At concentrations of 3.2 mM or higher after 20 h, and especially 48 h, increased cytotoxic effects were associated with disregulation of numerous genes related to oxidative stress, DNA repair and apoptosis. Primary human hepatocyte cultures were also exposed to 1 and 3.2 mM PB for 20 h and the changes were comparable to those found in HepaRG cells treated under the same conditions. Taken altogether, our data provide further evidence that HepaRG cells closely resemble primary human hepatocytes and provide new information on the effects of PB in human liver. These data also emphasize the importance of investigating dose- and time-dependent effects of chemicals when using toxicogenomic approaches.

Evolving concepts in liver tissue modeling and implications for in vitro toxicology

The development of human cell models stably expressing functional properties of the in vivo cells they are derived from for predicting toxicity of chemicals is a major challenge. For mimicking the liver, a major target of toxic chemicals, primary hepatocytes represent the most pertinent model. Their use is limited by interdonor functional variability and early phenotypic changes although their lifespan can be extended not only by culturing in a 2D dimension under sophisticated conditions but also by the use of synthetic and natural scaffolds as 3D supporting templates that allow cells to have a more stable microenvironment. Hepatocytes derived from stem cells could be the most appropriate alternative but up to now only liver progenitors/hepatoblasts are obtained in vitro. A few hepatocyte cell lines have retained a variable set of liver-specific functions. Among them are the human hepatoma HepaRG cells that express drug metabolism capacity at levels close to those found in primary hepatocytes making them a suitable model for both acute and chronic toxicity studies. New screening strategies are now proposed based on miniaturized and automated systems; they include the use of microfluidic chips and cell chips coupled with high content imaging analysis. Toxicogenomics technologies (particularly toxicotranscriptomics) have emerged as promising in vitro approaches for better identification and discrimination of cellular responses to chemicals. They should allow to discriminate compounds on the basis of the identification of a set of markers and/specific signaling pathways.