cDNA microarray gene expression profiling of hydroxyurea, paclitaxel, and p-anisidine, genotoxic compounds with differing tumorigenicity results

HbF Levels in Sickle Cell Disease Are Associated with Proportion of Circulating Hematopoietic Stem and Progenitor Cells and CC-Chemokines

The concentration of circulating hematopoietic stem and progenitor cells has not been studied longitudinally. Here, we report that the proportions of Lin-CD34+38- hematopoietic multipotent cells (HMCs) and of Lin-CD34+CD38+ hematopoietic progenitors cells (HPCs) are highly variable between individuals but stable over long periods of time, in both healthy individuals and sickle cell disease (SCD) patients. This suggests that these proportions are regulated by genetic polymorphisms or by epigenetic mechanisms. We also report that in SCD patients treated with hydroxyurea, the proportions of circulating HMCs and HPCs show a strong positive and negative correlation with fetal hemoglobin (HbF) levels, respectively. Titration of 65 cytokines revealed that the plasma concentration of chemokines CCL2, CCL11, CCL17, CCL24, CCL27, and PDGF-BB were highly correlated with the proportion of HMCs and HPCs and that a subset of these cytokines were also correlated with HbF levels. A linear model based on four of these chemokines could explain 80% of the variability in the proportion of circulating HMCs between individuals. The proportion of circulating HMCs and HPCs and the concentration of these chemokines might therefore become useful biomarkers for HbF response to HU in SCD patients. Such markers might become increasingly clinically relevant, as alternative treatment modalities for SCD are becoming available.

Variation in Gamma-Globin Expression before and after Induction with Hydroxyurea Associated with BCL11A, KLF1 and TAL1

The molecular mechanisms governing γ-globin expression in a subset of fetal hemoglobin (α2γ2: HbF) expressing red blood cells (F-cells) and the mechanisms underlying the variability of response to hydroxyurea induced γ-globin expression in the treatment of sickle cell disease are not completely understood. Here we analyzed intra-person clonal populations of basophilic erythroblasts (baso-Es) derived from bone marrow common myeloid progenitors in serum free cultures and report the level of fetal hemoglobin production in F-cells negatively correlates with expression of BCL11A, KLF1 and TAL1. We then examined the effects of hydroxyurea on these three transcription factors and conclude that a successful induction of γ-globin includes a reduction in BCL11A, KLF1 and TAL1 expression. These data suggests that expression changes in this transcription factor network modulate γ-globin expression in F-cells during steady state erythropoiesis and after induction with hydroxyurea.

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.

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: transcription profiling for toxicology assessment

Toxicogenomics, the application of transcription profiling to toxicology, has been widely used for elucidating the molecular and cellular actions of chemicals and other environmental stressors on biological systems, predicting toxicity before any functional damages, and classification of known or new toxicants based on signatures of gene expression. The success of a toxicogenomics study depends upon close collaboration among experts in different fields, including a toxicologist or biologist, a bioinformatician, statistician, physician and, sometimes, mathematician. This review is focused on toxicogenomics studies, including transcription profiling technology, experimental design, significant gene extraction, toxicological results interpretation, potential pathway identification, database input and the applications of toxicogenomics in various fields of toxicological study.

Microarray analyses of hypoxia-regulated genes in an aryl hydrocarbon receptor nuclear translocator (Arnt)-dependent manner

We investigated hypoxia-inducible factor (HIF)-dependent changes in the expression of 5592 genes in response to hypoxia (0.1% O(2), 16 h) by performing cDNA microarray analyses of mouse hepa1c1c7 and BpRc1 cells. BpRc1 cells are a hepa1c1c7 variant defective in HIF-beta/aryl hydrocarbon receptor nuclear translocator (Arnt), and are therefore unable to induce HIF target genes in response to hypoxia. By comparing hepa1c1c7 cells with BpRc1 cells, we were able to investigate hypoxia-regulated gene expression as well as the role played by HIF in regulating the hypoxic-dependent response of gene expression. This study identified 50 hypoxia-induced genes and 36 hypoxia-repressed genes. Quantitative PCR analysis of nine genes confirmed our ability to accurately analyze changes in hypoxia-induced gene expression by microarray analysis. By comparing quantitative PCR analyses of these nine genes in BpRc1 and hepa1c1c7 cells, we determined that eight of the nine hypoxia-induced genes are Arnt dependent. Additional quantitative PCR analyses of eight hypoxia-repressed genes confirmed, with a 50% probability, that microarray analysis was able to predict hypoxia-repressed gene expression. Only two of the four confirmed genes were found to be repressed in an Arnt-dependent manner. Collectively, six of these 13 genes (46.2% probability) showed a pattern of expression consistent with the microarray analysis with regard to Arnt dependence. Finally, we investigated the HIF-1alpha dependence of these 13 genes by quantitative PCR analysis in HIF-1alpha knockdown 3T3-L1 cells. These analyses identified novel hypoxia-regulated genes and confirmed the role of Arnt and HIF-1alpha in regulating their expression. These results identify additional HIF target genes and provide a more complete understanding of hypoxia signaling.

Application of toxicogenomics to genetic toxicology risk assessment

Based on the assumption that compounds having similar toxic modes of action induce specific gene expression changes, the toxicity of unknown compounds can be predicted after comparison of their molecular fingerprints with those obtained with compounds of known toxicity. These predictive models will therefore rely on the characterization of marker genes. Toxicogenomics (TGX) also provides mechanistic insight into the mode of toxicity, and can therefore be used as an adjunct to the standard battery of genotoxicity tests. Promising results, highlighting the ability of TGX to differentiate genotoxic from non-genotoxic carcinogens, as well as DNA-reactive from non-DNA reactive genotoxins, have been reported. Additional data suggested the possibility of ranking genotoxins according to the nature of their interactions with DNA. This new approach could contribute to the improvement of risk assessment. TGX could be applied as a follow-up testing strategy in case of positive in vitro genotoxicity findings, and could contribute to improve our ability to identify the molecular mechanism of action and to possibly better assess dose-response curves. TGX has been found to be less sensitive than the standard genotoxicity end-points, probably because it measures the whole cell population response, when compared with standard tests designed to detect rare events in a small number of cells. Further validation will be needed (1) to better link the profiles obtained with TGX to the established genotoxicity end-points, (2) to improve the gene annotation tools, and (3) to standardise study design and data analysis and to better evaluate the impact of variability between platforms and bioinformatics approaches.

DNA microarray technology in toxicogenomics of aquatic models: methods and applications

Toxicogenomics represents the merging of toxicology with genomics and bioinformatics to investigate biological functions of genome in response to environmental contaminants. Aquatic species have traditionally been used as models in toxicology to characterize the actions of environmental stresses. Recent completion of the DNA sequencing for several fish species has spurred the development of DNA microarrays allowing investigators access to toxicogenomic approaches. However, since microarray technology is thus far limited to only a few aquatic species and derivation of biological meaning from microarray data is highly dependent on statistical arguments, the full potential of microarray in aquatic species research has yet to be realized. Herein we review some of the issues related to construction, probe design, statistical and bioinformatical data analyses, and current applications of DNA microarrays. As a model a recently developed medaka (Oryzias latipes) oligonucleotide microarray was described to highlight some of the issues related to array technology and its application in aquatic species exposed to hypoxia. Although there are known non-biological variations present in microarray data, it remains unquestionable that array technology will have a great impact on aquatic toxicology. Microarray applications in aquatic toxicogenomics will range from the discovery of diagnostic biomarkers, to establishment of stress-specific signatures and molecular pathways hallmarking the adaptation to new environmental conditions.

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.

Animal carcinogenicity studies: 3. Alternatives to the bioassay

Conventional animal carcinogenicity tests take around three years to design, conduct and interpret. Consequently, only a tiny fraction of the thousands of industrial chemicals currently in use have been tested for carcinogenicity. Despite the costs of hundreds of millions of dollars and millions of skilled personnel hours, as well as millions of animal lives, several investigations have revealed that animal carcinogenicity data lack human specificity (i.e. the ability to identify human non-carcinogens), which severely limits the human predictivity of the bioassay. This is due to the scientific inadequacies of many carcinogenicity bioassays, and numerous serious biological obstacles, which render profoundly difficult any attempts to accurately extrapolate animal data in order to predict carcinogenic hazards to humans. Proposed modifications to the conventional bioassays have included the elimination of mice as a second species, and the use of genetically-altered or neonatal mice, decreased study durations, initiation-promotion models, the greater incorporation of toxicokinetic and toxicodynamic assessments, structure-activity relationship (computerised) systems, in vitro assays, cDNA microarrays for detecting changes in gene expression, limited human clinical trials, and epidemiological research. The potential advantages of non-animal assays when compared to bioassays include the superior human specificity of the results, substantially reduced time-frames, and greatly reduced demands on financial, personnel and animal resources. Inexplicably, however, the regulatory agencies have been frustratingly slow to adopt alternative protocols. In order to decrease the enormous cost of cancer to society, a substantial redirection of resources away from excessively slow and resource-intensive rodent bioassays, into the further development and implementation of non-animal assays, is both strongly justified and urgently required.

RUNX3 suppresses gastric epithelial cell growth by inducing p21(WAF1/Cip1) expression in cooperation with transforming growth factor {beta}-activated SMAD

RUNX3 has been suggested to be a tumor suppressor of gastric cancer. The gastric mucosa of the Runx3-null mouse develops hyperplasia due to enhanced proliferation and suppressed apoptosis accompanied by a decreased sensitivity to transforming growth factor beta1 (TGF-beta1). It is known that TGF-beta1 induces cell growth arrest by activating CDKN1A (p21(WAF1)(/Cip1)), which encodes a cyclin-dependent kinase inhibitor, and this signaling cascade is considered to be a tumor suppressor pathway. However, the lineage-specific transcription factor that cooperates with SMADs to induce p21 expression is not known. Here we show that RUNX3 is required for the TGF-beta-dependent induction of p21 expression in stomach epithelial cells. Overexpression of RUNX3 potentiates TGF-beta-dependent endogenous p21 induction. In cooperation with SMADs, RUNX3 synergistically activates the p21 promoter. In contrast, RUNX3-R122C, a mutation identified in a gastric cancer patient, abolished the ability to activate the p21 promoter or cooperate with SMADs. Furthermore, areas in mouse and human gastric epithelium where RUNX3 is expressed coincided with those where p21 is expressed. Our results suggest that at least part of the tumor suppressor activity of RUNX3 is associated with its ability to induce p21 expression.

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.

Systems-ADME/Tox: resources and network approaches

The increasing cost of drug development is partially due to our failure to identify undesirable compounds at an early enough stage of development. The application of higher throughput screening methods have resulted in the generation of very large datasets from cells in vitro or from in vivo experiments following the treatment with drugs or known toxins. In recent years the development of systems biology, databases and pathway software has enabled the analysis of the high-throughput data in the context of the whole cell. One of the latest technology paradigms to be applied alongside the existing in vitro and computational models for absorption, distribution, metabolism, excretion and toxicology (ADME/Tox) involves the integration of complex multidimensional datasets, termed toxicogenomics. The goal is to provide a more complete understanding of the effects a molecule might have on the entire biological system. However, due to the sheer complexity of this data it may be necessary to apply one or more different types of computational approaches that have as yet not been fully utilized in this field. The present review describes the data generated currently and introduces computational approaches as a component of ADME/Tox. These methods include network algorithms and manually curated databases of interactions that have been separately classified under systems biology methods. The integration of these disparate tools will result in systems-ADME/Tox and it is important to understand exactly what data resources and technologies are available and applicable. Examples of networks derived with important drug transporters and drug metabolizing enzymes are provided to demonstrate the network technologies.

Identification of potential biomarkers of genotoxicity and carcinogenicity in L5178Y mouse lymphoma cells by cDNA microarray analysis

In the present study, cDNA microarray analyses were performed with mouse cDNA chips in order to evaluate similarities and differences in the gene expression profiles for compounds differing in their genotoxic and carcinogenic potential. Eight test substances were evaluated, two each from four classes of compounds: genotoxic carcinogens (1,2-dibromoethane and glycidol), genotoxic noncarcinogens (8-hydroxyquinoline and emodin), nongenotoxic carcinogens (methyl carbamate and o-nitrotoluene), and nongenotoxic noncarcinogens (D-mannitol and 1,2-dichlorobenzene). Quadruplicate hybridization experiments were performed in order to identify a set of genes with significant expression changes for these four classes of substances. Twelve genes were consistently altered more than twofold by the genotoxic noncarcinogens while four genes were consistently regulated by the nongenotoxic carcinogens. One gene (Trp63) was identified whose expression was upregulated by all four genotoxic substances regardless of the presence or absence of carcinogenicity; this finding, however, was not confirmed by quantitative real-time RT-PCR. RT-PCR did confirm the change in expression of 9 of 15 genes (60%) identified by microarray analysis. Interestingly, the downregulated genes were least likely to be validated by real-time RT-PCR. Those genes showing more than a twofold change in expression level in response to at least one substance were further analyzed with hierarchical clustering after category assignment of each gene according to its main cellular function. Clustering revealed differences in the gene expression profiles between the genotoxic and nongenotoxic substances for genes involved in cell cycle control, the stress response, and the immune response. However, no clustering specific to all four carcinogenic substances was observed in any of the functional categories. Taken together, these results suggest that gene expression profiling in mouse lymphoma cells can provide valuable information for the evaluation of potential genotoxicity but may have limitations in predicting carcinogenicity.

Applications of microarrays with toxicologically relevant genes (tox genes) for the evaluation of chemical toxicants in Sprague Dawley rats in vivo and human hepatocytes in vitro

Microarrays with toxicologically relevant genes (tox genes) have been developed in our laboratory for toxicogenomics studies in rat, dog and man. The genes were chosen using published information as well as a discovery process for genes responsive to toxic treatments using transcription profiling experiments conducted with rats and dogs. In addition to published information human tox genes were derived from rat tox genes based on gene homology. Using the microarray with rat-specific tox genes, a database containing gene expression, histopathology, and clinical chemistry findings has been generated for 89 compounds. Analysis of the database indicates that treatment with toxic compounds induces specific gene expression patterns. Dose- and time-dependent response relationships in gene expression were observed for treatment with toxic compounds. Gene expression at 24h was found to correlate well with organ toxicity observed at 72 h. Mining of the database led to the selection of specific groups of genes (predictive gene sets) whose expression patterns are predictive of organ toxicity with a high degree of accuracy (approximately 90%). The data also provide insight on toxic mechanism and gene regulation pathways. For instance, carbon tetrachloride and chloroform treatments were found to decrease the expression of the cytochrome P450 isoform 3A1 gene while enhancing the expression of the multiple drug resistance gene MDR1 in liver, clearly demonstrating that the CYP3A1 and MDR1 genes were not co-regulated as postulated by some researchers. This approach, the use of gene expression as an endpoint to define organ toxicity, is extended to the definition of human drug toxicity using primary human hepatocytes as a test system. Preliminary results demonstrate that the toxic drug, troglitazone, can be clearly distinguished from the less toxic analogues, rosiglitazone and pioglitazone based on their effects on tox gene expression in human hepatocytes. Our results with both rats in vivo and human hepatocytes in vitro suggest that microarrays with toxicologically relevant genes can be used routinely for the evaluation of chemical toxicity.