Comparison of gene expression changes induced in mouse and human cells treated with direct-acting mutagens

Mitotic syndicates Aurora Kinase B (AURKB) and mitotic arrest deficient 2 like 2 (MAD2L2) in cohorts of DNA damage response (DDR) and tumorigenesis

Aurora Kinase B (AURKB) and Mitotic Arrest Deficient 2 Like 2 (MAD2L2) are emerging anticancer therapeutic targets. AURKB and MAD2L2 are the least well studied members of their protein families, compared to AURKA and MAD2L1. Both AURKB and MAD2L2 play a critical role in mitosis, cell cycle checkpoint, DNA damage response (DDR) and normal physiological processes. However, the oncogenic roles of AURKB and MAD2L2 in tumorigenesis and genomic instability have also been reported. DDR acts as an arbitrator for cell fate by either repairing the damage or directing the cell to self-destruction. While there is strong evidence of interphase DDR, evidence of mitotic DDR is just emerging and remains largely unelucidated. To date, inhibitors of the DDR components show effective anti-cancer roles. Contrarily, long-term resistance towards drugs that target only one DDR target is becoming a challenge. Targeting interactions between protein-protein or protein-DNA holds prominent therapeutic potential. Both AURKB and MAD2L2 play critical roles in the success of mitosis and their emerging roles in mitotic DDR cannot be ignored. Small molecule inhibitors for AURKB are in clinical trials. A few lead compounds towards MAD2L2 inhibition have been discovered. Targeting mitotic DDR components and their interaction is emerging as a potent next generation anti-cancer therapeutic target. This can be done by developing small molecule inhibitors for AURKB and MAD2L2, thereby targeting DDR components as anti-cancer therapeutic targets and/or targeting mitotic DDR. This review focuses on AURKB and MAD2L2 prospective synergy to deregulate the p53 DDR pathway and promote favourable conditions for uncontrolled cell proliferation.

TGx-DDI, a Transcriptomic Biomarker for Genotoxicity Hazard Assessment of Pharmaceuticals and Environmental Chemicals

Genotoxicity testing is an essential component of the safety assessment paradigm required by regulatory agencies world-wide for analysis of drug candidates, and environmental and industrial chemicals. Current genotoxicity testing batteries feature a high incidence of irrelevant positive findings—particularly for in vitro chromosomal damage (CD) assays. The risk management of compounds with positive in vitro findings is a major challenge and requires complex, time consuming, and costly follow-up strategies including animal testing. Thus, regulators are urgently in need of new testing approaches to meet legislated mandates. Using machine learning, we identified a set of transcripts that responds predictably to DNA-damage in human cells that we refer to as the TGx-DDI biomarker, which was originally referred to as TGx-28.65. We proposed to use this biomarker in conjunction with current genotoxicity testing batteries to differentiate compounds with irrelevant “false” positive findings in the in vitro CD assays from true DNA damaging agents (i.e., for de-risking agents that are clastogenic in vitro but not in vivo). We validated the performance of the TGx-DDI biomarker to identify true DNA damaging agents, assessed intra- and inter- laboratory reproducibility, and cross-platform performance. Recently, to augment the application of this biomarker, we developed a high-throughput cell-based genotoxicity testing system using the NanoString nCounter® technology. Here, we review the status of TGx-DDI development, its integration in the genotoxicity testing paradigm, and progress to date in its qualification at the US Food and Drug Administration (FDA) as a drug development tool. If successfully validated and implemented, the TGx-DDI biomarker assay is expected to significantly augment the current strategy for the assessment of genotoxic hazards for drugs and chemicals.

Naked mole rat cells display more efficient excision repair than mouse cells

Naked mole rat (NMR) is the long-lived and tumor-resistant rodent. NMRs possess multiple adaptations that may contribute to longevity and cancer-resistance. However, whether NMRs have more efficient DNA repair have not been directly tested. Here we compared base excision repair (BER) and nucleotide excision repair (NER) systems in extracts from NMR and mouse fibroblasts after UVC irradiation. Transcript levels of the key repair enzymes demonstrated in most cases higher inducibility in the mouse vs the NMR cells. Ratios of repair enzymes activities in the extracts somewhat varied depending on post-irradiation time. NMR cell extracts were 2–3-fold more efficient at removing the bulky lesions, 1.5–3-fold more efficient at removing uracil, and about 1.4-fold more efficient at cleaving the AP-site than the mouse cells, while DNA polymerase activities being as a whole higher in the mouse demonstrate different patterns of product distribution. The level of poly(ADP-ribose) synthesis was 1.4–1.8-fold higher in the NMR cells. Furthermore, NMR cell extracts displayed higher binding of PARP1 to DNA probes containing apurinic/apyrimidinic site or photo-reactive DNA lesions. Cumulatively, our results suggest that the NMR has more efficient excision repair systems than the mouse, which may contribute to longevity and cancer resistance of this species.

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.

Transcriptional regulation of human DNA repair genes following genotoxic stress: trigger mechanisms, inducible responses and genotoxic adaptation

DNA repair is the first barrier in the defense against genotoxic stress. In recent years, mechanisms that recognize DNA damage and activate DNA repair functions through transcriptional upregulation and post-translational modification were the focus of intensive research. Most DNA repair pathways are complex, involving many proteins working in discrete consecutive steps. Therefore, their balanced expression is important for avoiding erroneous repair that might result from excessive base removal and DNA cleavage. Amelioration of DNA repair requires both a fine-tuned system of lesion recognition and transcription factors that regulate repair genes in a balanced way. Transcriptional upregulation of DNA repair genes by genotoxic stress is counteracted by DNA damage that blocks transcription. Therefore, induction of DNA repair resulting in an adaptive response is only visible through a narrow window of dose. Here, we review transcriptional regulation of DNA repair genes in normal and cancer cells and describe mechanisms of promoter activation following genotoxic exposures through environmental carcinogens and anticancer drugs. The data available to date indicate that 25 DNA repair genes are subject to regulation following genotoxic stress in rodent and human cells, but for only a few of them, the data are solid as to the mechanism, homeostatic regulation and involvement in an adaptive response to genotoxic stress.

Validation of microarray data in human lymphoblasts shows a role of the ubiquitin-proteasome system and NF-kB in the pathogenesis of Down syndrome

BACKGROUND

Down syndrome (DS) is a complex disorder caused by the trisomy of either the entire, or a critical region of chromosome 21 (21q22.1-22.3). Despite representing the most common cause of mental retardation, the molecular bases of the syndrome are still largely unknown.

METHODS

To better understand the pathogenesis of DS, we analyzed the genome-wide transcription profiles of lymphoblastoid cell lines (LCLs) from six DS and six euploid individuals and investigated differential gene expression and pathway deregulation associated with trisomy 21. Connectivity map and PASS-assisted exploration were used to identify compounds whose molecular signatures counteracted those of DS lymphoblasts and to predict their therapeutic potential. An experimental validation in DS LCLs and fetal fibroblasts was performed for the most deregulated GO categories, i.e. the ubiquitin mediated proteolysis and the NF-kB cascade.

RESULTS

We show, for the first time, that the level of protein ubiquitination is reduced in human DS cell lines and that proteasome activity is increased in both basal conditions and oxidative microenvironment. We also provide the first evidence that NF-kB transcription levels, a paradigm of gene expression control by ubiquitin-mediated degradation, is impaired in DS due to reduced IkB-alfa ubiquitination, increased NF-kB inhibitor (IkB-alfa) and reduced p65 nuclear fraction. Finally, the DSCR1/DYRK1A/NFAT genes were analysed. In human DS LCLs, we confirmed the presence of increased protein levels of DSCR1 and DYRK1A, and showed that the levels of the transcription factor NFATc2 were decreased in DS along with a reduction of its nuclear translocation upon induction of calcium fluxes.

CONCLUSIONS

The present work offers new perspectives to better understand the pathogenesis of DS and suggests a rationale for innovative approaches to treat some pathological conditions associated to DS.

Characterization of formaldehyde's genotoxic mode of action by gene expression analysis in TK6 cells

Gene expression analysis has been established as a tool for the characterization of genotoxic mechanisms of chemical mutagens. It has been suggested that expression analysis is capable of distinguishing compounds that cause DNA damage from those that interfere with mitotic spindle function. Formaldehyde (FA) is known to be a DNA-reactive substance which mainly induces chromosomal damage in cultured mammalian cells. However, there has been concern that FA might also induce leukemia-specific aneuploidies, although recent cytogenetic studies excluded a relevant aneugenic potential of FA. We now investigated whether gene expression profiling can be used as a molecular tool to further characterize FA's genotoxic mode of action and to differentiate between clastogenic and aneugenic activity. TK6 cells were exposed to FA for 4 and 24 h, and changes in gene expression were analyzed using a whole-genome human microarray. Results were compared to the expression profiles of two DNA-damaging clastogens (methyl methanesulfonate and ethyl methanesulfonate) and two aneugens (colcemid and vincristine). The genotoxic activity of FA, MMS and EMS under these conditions was confirmed by comet assay experiments. The gene expression profiles indicated that clastogens and aneugens induce discriminable gene expression patterns. Exposure of TK6 cells to FA led to a discrete gene expression pattern, and all toxicogenomics analyses revealed a closer relationship of FA with clastogens than with aneugens.

MondoA is highly overexpressed in acute lymphoblastic leukemia cells and modulates their metabolism, differentiation and survival

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. To identify novel candidates for targeted therapy, we performed a comprehensive transcriptome analysis identifying MondoA (MLXIP) - a transcription factor regulating glycolysis - to be overexpressed in ALL compared to normal tissues. Using microarray-profiling, gene-set enrichment analysis, RNA interference and functional assays we show that MondoA overexpression increases glucose catabolism and maintains a more immature phenotype, which is associated with enhanced survival and clonogenicity of leukemia cells. These data point to an important contribution of MondoA to leukemia aggressiveness and make MondoA a potential candidate for targeted treatment of ALL.

Effect of chemical mutagens and carcinogens on gene expression profiles in human TK6 cells

Characterization of toxicogenomic signatures of carcinogen exposure holds significant promise for mechanistic and predictive toxicology. In vitro transcriptomic studies allow the comparison of the response to chemicals with diverse mode of actions under controlled experimental conditions. We conducted an in vitro study in TK6 cells to characterize gene expression signatures of exposure to 15 genotoxic carcinogens frequently used in European industries. We also examined the dose-responsive changes in gene expression, and perturbation of biochemical pathways in response to these carcinogens. TK6 cells were exposed at 3 dose levels for 24 h with and without S9 human metabolic mix. Since S9 had an impact on gene expression (885 genes), we analyzed the gene expression data from cells cultures incubated with S9 and without S9 independently. The ribosome pathway was affected by all chemical-dose combinations. However in general, no similar gene expression was observed among carcinogens. Further, pathways, i.e. cell cycle, DNA repair mechanisms, RNA degradation, that were common within sets of chemical-dose combination were suggested by clustergram. Linear trends in dose–response of gene expression were observed for Trichloroethylene, Benz[a]anthracene, Epichlorohydrin, Benzene, and Hydroquinone. The significantly altered genes were involved in the regulation of (anti-) apoptosis, maintenance of cell survival, tumor necrosis factor-related pathways and immune response, in agreement with several other studies. Similarly in S9+ cultures, Benz[a]pyrene, Styrene and Trichloroethylene each modified over 1000 genes at high concentrations. Our findings expand our understanding of the transcriptomic response to genotoxic carcinogens, revealing the alteration of diverse sets of genes and pathways involved in cellular homeostasis and cell cycle control.

Identification of biomarkers of chemically induced hepatocarcinogenesis in rasH2 mice by toxicogenomic analysis

Toxicogenomic approaches have been applied to chemical-induced heptocarcinogenesis rodent models for the identification of biomarkers of early-stage hepatocarcinogenesis and to help clarify the underlying carcinogenic mechanisms in the liver. In this study, we used toxiciogenomic methods to identify candidate biomarker genes associated with hepatocarcinogenesis in rasH2 mice. Blood chemical, histopathologic, and gene expression analyses of the livers of rasH2 mice were performed 7 and 91 days after the administration of the genotoxic hepatocarcinogens 2-acetylaminofluorene (AAF) and diethylnitrosoamine (DEN), the genotoxic carcinogen melphalan (Mel), and the nongenotoxic noncarcinogen 1-naphthylisothiocynate (ANIT). Histopathologic lesions and a rise in accompanying serum marker levels were found in the DEN-treated rasH2 mice, whereas no neoplastic lesions were observed in the rasH2 mice. However, biological functional analysis using Ingenuity Pathways Analysis (IPA) software revealed that genes with comparable molecular and cellular functions were similarly deregulated in the AAF- and DEN-treated rasH2 mice. We selected 68 significantly deregulated genes that represented a hepatocarcinogen-specific signature; these genes were commonly deregulated in both the AAF- and DEN-treated rasH2 mice on days 7 and 91. Hierarchical clustering analysis indicated that the expression patterns of the selected genes in the hepatocarcinogen (AAF and DEN) groups were distinctive from the patterns in the control, Mel, and ANIT groups. Biomarker filter analysis using IPA software suggested that 28 of the 68 signature genes represent promising candidate biomarkers of cancer. Quantitative real-time PCR analysis confirmed that the deregulated genes, which exhibited sustained up- and down-regulation up to day 91, are likely involved in early-stage hepatocarcinogenesis. In summary, the common and significant gene expression changes induced by AAF and DEN may reflect early molecular events associated with hepatocarcinogenesis, and these "signature" genes may be useful as biomarkers of hepatocarcinogenesis in mice.

Study of gene expression profiles in TK6 human cells exposed to DNA-oxidizing agents

During the last decade, there has been clear progress in using threshold in risk assessment but its acceptance by scientists is still under debate. Contrary to indirect DNA-damaging agents, DNA-reactive agents have been assumed to have a non-threshold mode of action, as they directly induce DNA lesions that potentially can be converted into mutations. However, in recent years there is a growing number of data establishing threshold doses even for these DNA-reactive compounds. Indeed, there are several defence and repair mechanisms that provide protection and that may be responsible for genotoxic thresholds. In this context, we recently showed that DNA-oxidizing agents exhibit a thresholded dose-response in vitro with respect to chromosomal alterations. We have hypothesized the involvement of different cellular responses whose nature and efficiency depend on the stress level. The aim of this study was to develop a more complete understanding of these underlying mechanisms. We investigated global gene expression profiles of human lymphoblastoid TK6 cells after exposure to potassium bromate and hydrogen peroxide (via glucose oxidase). Cells were treated for 1h and mRNAs were isolated either immediately at the end of the treatment or after a 23-h recovery period. Our results showed that cells have developed elaborate cellular responses to oxidative stress in order to maintain genomic integrity. Many of altered genes were redox-sensitive transcription factors such as p53, NF-kappaB, AP-1 and Nrf2. Their downstream target genes and signalling pathways were subsequently activated leading mainly to the induction of antioxidant defenses, inflammation, cell cycle arrest, DNA repair and cell death. Overall, our study allowed the identification of key events involved in the thresholded response observed after DNA-oxidizing agents exposure and shows the usefulness of the combination of standard in vitro genotoxicity assays with gene expression profiling technology to determine modes of action, particularly for critical risk assessment.

Effects of ionizing radiation on self-renewal and pluripotency of human embryonic stem cells

Human embryonic stem cells (hESC) present a novel platform for in vitro investigation of the early embryonic cellular response to ionizing radiation. Thus far, no study has analyzed the genome-wide transcriptional response to ionizing radiation in hESCs, nor has any study assessed their ability to form teratomas, the definitive test of pluripotency. In this study, we use microarrays to analyze the global gene expression changes in hESCs after low-dose (0.4 Gy), medium-dose (2 Gy), and high-dose (4 Gy) irradiation. We identify genes and pathways at each radiation dose that are involved in cell death, p53 signaling, cell cycling, cancer, embryonic and organ development, and others. Using Gene Set Enrichment Analysis, we also show that the expression of a comprehensive set of core embryonic transcription factors is not altered by radiation at any dose. Transplantation of irradiated hESCs to immune-deficient mice results in teratoma formation from hESCs irradiated at all doses, definitive proof of pluripotency. Further, using a bioluminescence imaging technique, we have found that irradiation causes hESCs to initially die after transplantation, but the surviving cells quickly recover by 2 weeks to levels similar to control. To conclude, we show that similar to somatic cells, irradiated hESCs suffer significant death and apoptosis after irradiation. However, they continue to remain pluripotent and are able to form all three embryonic germ layers. Studies such as this will help define the limits for radiation exposure for pregnant women and also radiotracer reporter probes for tracking cellular regenerative therapies.

Differential gene expression in primary human skin keratinocytes and fibroblasts in response to ionizing radiation

Although skin is usually exposed during human exposures to ionizing radiation, there have been no thorough examinations of the transcriptional response of skin fibroblasts and keratinocytes to radiation. The transcriptional response of quiescent primary fibroblasts and keratinocytes exposed to from 10 cGy to 5 Gy and collected 4 h after treatment was examined. RNA was isolated and examined by microarray analysis for changes in the levels of gene expression. Exposure to ionizing radiation altered the expression of 279 genes across both cell types. Changes in RNA expression could be arranged into three main categories: (1) changes in keratinocytes but not in fibroblasts, (2) changes in fibroblasts but not in keratinocytes, and (3) changes in both. All of these changes were primarily of p53 target genes. Similar radiation-induced changes were induced in immortalized fibroblasts or keratinocytes. In separate experiments, protein was collected and analyzed by Western blotting for expression of proteins observed in microarray experiments to be overexpressed at the mRNA level. Both Q-PCR and Western blot analysis experiments validated these transcription changes. Our results are consistent with changes in the expression of p53 target genes as indicating the magnitude of cell responses to ionizing radiation.

Characterization and interlaboratory comparison of a gene expression signature for differentiating genotoxic mechanisms

The genotoxicity testing battery is highly sensitive for detection of chemical carcinogens. However, it features a low specificity and provides only limited mechanistic information required for risk assessment of positive findings. This is especially important in case of positive findings in the in vitro chromosome damage assays, because chromosome damage may be also induced secondarily to cell death. An increasing body of evidence indicates that toxicogenomic analysis of cellular stress responses provides an insight into mechanisms of action of genotoxicants. To evaluate the utility of such a toxicogenomic analysis we evaluated gene expression profiles of TK6 cells treated with four model genotoxic agents using a targeted high density real-time PCR approach in a multilaboratory project coordinated by the Health and Environmental Sciences Institute Committee on the Application of Genomics in Mechanism-based Risk Assessment. We show that this gene profiling technology produced reproducible data across laboratories allowing us to conclude that expression analysis of a relevant gene set is capable of distinguishing compounds that cause DNA adducts or double strand breaks from those that interfere with mitotic spindle function or that cause chromosome damage as a consequence of cytotoxicity. Furthermore, our data suggest that the gene expression profiles at early time points are most likely to provide information relevant to mechanisms of genotoxic damage and that larger gene expression arrays will likely provide richer information for differentiating molecular mechanisms of action of genotoxicants. Although more compounds need to be tested to identify a robust molecular signature, this study confirms the potential of toxicogenomic analysis for investigation of genotoxic mechanisms.

Toxicogenomics: overview and potential applications for the study of non-covalent DNA interacting chemicals

Non-covalent DNA interacting agents, DNA-groove binding chemicals and DNA intercalators, are generally considered less cytotoxic than agents producing covalent DNA adducts and other DNA damage. Although the impact of non-covalent compound-DNA interactions on convoluted molecular and biochemical pathways is not well characterized, the most prominent effects include DNA conformational and related structural perturbations, interference with normal DNA protein interactions, such as topoisomerases, as well as effects on mitochondrial DNA and function. The cellular responses to such perturbations would be expected to include changes in transcription of a variety of genes. The emerging field of toxicogenomics seeks to exploit gene responses to define expression profiling signatures for various types of drugs and toxicants, and to provide mechanistic insight into their cellular effects. There are a variety of examples whereby different classes of genotoxicants and non-genotoxic agents can be distinguished by gene expression profiling using functional genomics approaches, which survey global transcriptional responses. In this review, we will discuss the promises and precautions in the use of functional genomics approaches to characterize stress agents including non-covalent DNA interacting agents.

Gene expression profile analysis: an emerging approach to investigate mechanisms of genotoxicity

The response to stress triggers transcriptional activation of genes involved in cell survival and/or cell death. Thus, the monitoring of gene expression levels in large gene sets or whole genomes in response to various agents (toxicogenomics) has been proposed as a tool for investigating mechanisms of toxicity. Although standard in vitro genetic toxicity testing provides relatively simple and accurate hazard detection, interpretation of positive findings, i.e., in vitro chromosome aberrations, in terms of relevant risk to humans is difficult, due to the limited insight into the underlying mechanisms. Therefore, the development of experimental approaches capable of differentiating a wide range of genotoxic mechanisms is expected to significantly improve risk assessment. The goal of this review is to summarize current developments in toxicogenomic analysis of genotoxic stress, and to provide a perspective on the application of gene expression profile analysis in genetic toxicology.

Bystander effect: biological endpoints and microarray analysis

In cell populations exposed to ionizing radiation, the biological effects occur in a much larger proportion of cells than are estimated to be traversed by radiation. It has been suggested that irradiated cells are capable of providing signals to the neighboring unirradiated cells resulting in damage to these cells. This phenomenon is termed the bystander effect. The bystander effect induces persistent, long-term, transmissible changes that result in delayed death and neoplastic transformation. Because the bystander effect is relevant to carcinogenesis, it could have significant implications for risk estimation for radiation exposure. The nature of the bystander effect signal and how it impacts the unirradiated cells remains to be elucidated. Examination of the changes in gene expression could provide clues to understanding the bystander effect and could define the signaling pathways involved in sustaining damage to these cells. The microarray technology serves as a tool to gain insight into the molecular pathways leading to bystander effect. Using medium from irradiated normal human diploid lung fibroblasts as a model system we examined gene expression alterations in bystander cells. The microarray data revealed that the radiation-induced gene expression profile in irradiated cells is different from unirradiated bystander cells suggesting that the pathways leading to biological effects in the bystander cells are different from the directly irradiated cells. The genes known to be responsive to ionizing radiation were observed in irradiated cells. Several genes were upregulated in cells receiving media from irradiated cells. Surprisingly no genes were found to be downregulated in these cells. A number of genes belonging to extracellular signaling, growth factors and several receptors were identified in bystander cells. Interestingly 15 genes involved in the cell communication processes were found to be upregulated. The induction of receptors and the cell communication processes in bystander cells receiving media from irradiated cells supports the active involvement of these processes in inducing bystander effect.

Comparison of the expression profiles induced by genotoxic and nongenotoxic carcinogens in rat liver

Application of recently developed gene expression techniques using microarrays in toxicological studies (toxicogenomics) facilitate the interpretation of a toxic compound's mode of action and may also allow the prediction of selected toxic effects based on gene expression changes. In order to test this hypothesis, we investigated whether carcinogens at doses known to induce liver tumors in the 2-year rat bioassay deregulate characteristic sets of genes in a short term in vivo study and whether these deregulated genes represent defined biological pathways. Male Wistar rats were dosed with the four nongenotoxic hepatocarcinogens methapyrilene (MPy, 60 mg/kg/day), diethylstilbestrol (DES, 10 mg/kg/day), Wy-14643 (Wy, 60 mg/kg/day), and piperonylbutoxide (PBO, 1200 mg/kg/day). After 1, 3, 7, and 14 days, the livers were taken for histopathological evaluation and for analysis of the gene expression profiles on Affymetrix RG_U34A arrays. The expression profile of the four nongenotoxic carcinogens were compared to the profiles of the four genotoxic carcinogens 2-nitrofluorene (2-NF), dimethylnitrosamine (DMN), 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), and aflatoxin B1 (AB1) from a similar study reported previously. By using statistical and clustering tools characteristically deregulated genes were extracted and functionally classified. Distinct cellular pathways were affected by the nongenotoxic carcinogens compared to the genotoxic carcinogens which at least partly correlated with the two-stage model of carcinogenesis. Characteristic to genotoxic carcinogens were a DNA damage response and the activation of proliferative and survival signaling. Nongenotoxic carcinogens showed responses to oxidative DNA or protein damage, as well as cell cycle progression and signs of regeneration. Many of the gene alterations found with the nongenotoxic carcinogens imply compound-specific mechanisms. Although neither a single gene nor a single pathway will be sufficient to discriminate the two classes of carcinogens, it became evident that combinations of pathway-associated gene expression profiles may be used to predict a genotoxic or nongenotoxic carcinogenic potential of a compound in short-term studies.