Genomic interrogation of mechanism(s) underlying cellular responses to toxicants

Expression of immune, antioxidant and stress related genes in different organs of common carp exposed to indoxacarb

The aim of the present study was to investigate the effects of chronic exposure of common carp (Cyprinus carpio) to indoxacarb on immune, antioxidant and stress gene expression. After 21 days exposure to 0, 0.75, 1.5 and 3 ppm indoxacarb, expression of IL-1β, IL-8, IL-10, TNF-α, IFN-γ, SOD, CAT, HSP70, IGF-I and IGF-II were assessed in liver, kidney and gills. In general, exposure to low concentration of indoxacarb increased inflammatory cytokine gene expression (IL-1β, IL-8, IL-10, TNF-α and IFN-γ) and inhibits inflammatory cytokines' expression at higher concentrations. The assessment of antioxidant gene expression (SOD and CAT) in different organs indicate that they were increased by low concentrations of indoxacarb to deal with primary oxidative situation. However, higher concentrations of indoxacarb caused reduction in oxidative gene expression. IGF genes expression in liver significantly increased at a concentration of 0.75 ppm treatment, then it decreased at 1.5 ppm indoxacarb and increased again by increasing in the indoxacarb concentration to 3 ppm. The expression of HSP70 in kidney showed a significant elevation in 0.75 and 1.5 ppm treatments compared with 3 ppm treatment and the control group. The expression of this gene in liver was significantly increased in 1.5 and 3 ppm treatments. The same pattern of expression was also observed in gill. Overall, indoxacarb exposure affects common carp health at transcription levels. Changes in the genes expression generally suggest that indoxacarb exposure led to interference in inflammation, oxidative stress and tissue damage.

Crosstalk between Receptor and Non-receptor Mediated Chemical Modes of Action in Rat Livers Converges through a Dysregulated Gene Expression Network at Tumor Suppressor Tp53

Chemicals, toxicants, and environmental stressors mediate their biologic effect through specific modes of action (MOAs). These encompass key molecular events that lead to changes in the expression of genes within regulatory pathways. Elucidating shared biologic processes and overlapping gene networks will help to better understand the toxicologic effects on biological systems. In this study we used a weighted network analysis of gene expression data from the livers of male Sprague-Dawley rats exposed to chemicals that elicit their effects through receptor-mediated MOAs (aryl hydrocarbon receptor, orphan nuclear hormone receptor, or peroxisome proliferator-activated receptor-α) or non-receptor-mediated MOAs (cytotoxicity or DNA damage). Four gene networks were highly preserved and statistically significant in each of the two MOA classes. Three of the four networks contain genes that enrich for immunity and defense. However, many canonical pathways related to an immune response were activated from exposure to the non-receptor-mediated MOA chemicals and deactivated from exposure to the receptor-mediated MOA chemicals. The top gene network contains a module with 33 genes including tumor suppressor TP53 as the central hub which was slightly up-regulated in gene expression compared to control. Although, there is crosstalk between the two MOA classes of chemicals at the TP53 gene network, more than half of the genes are dysregulated in opposite directions. For example, Thromboxane A Synthase 1 (Tbxas1), a cytochrome P450 protein coding gene regulated by Tp53, is down-regulated by exposure to the receptor-mediated chemicals but up-regulated by the non-receptor-mediated chemicals. The regulation of gene expression by the chemicals within MOA classes was consistent despite varying alanine transaminase and aspartate aminotransferase liver enzyme measurements. These results suggest that overlap in toxicologic pathways by chemicals with different MOAs provides common mechanisms for discordant regulation of gene expression within molecular networks.

Investigating the Mechanistic Basis for Hepatic Toxicity Induced by an Experimental Chemokine Receptor 5 (CCR5) Antagonist Using a Compendium of Gene Expression Profiles

A compendium of hepatic gene expression signatures was used to identify a mechanistic basis for the hepatic toxicity of an experimental CCR5 antagonist (MrkA). Development of MrkA, a potential HIV therapeutic, was discontinued due to hepatotoxicity in preclinical studies. Rats were treated with MrkA at 3 dose levels (50, 250, and 500 mg/kg) for 1, 3, or 7 days. Hepatic toxicity (vacuolation, consistent with steatosis, and elevated serum transaminase levels) was observed at 250 and 500 mg/kg, but not at 50 mg/kg. Hepatic gene expression profiles were compared to a compendium of hepatic expression profiles. MrkA was similar to 3 β-oxidation inhibitors (valproate, cyclopropane carboxylate, pivalate), 8 PPARα agonists (fenofibrate, bezafibrate and 6 fibrate analogues), and 3 other diverse compounds (diethylnitrosamine, microcystin LR & actinomycin D). These data indicate MrkA to be a mitochondrial inhibitor, and activation of PPARα-regulated transcription was thought to be due to an accumulation of endogenous ligands. While mitochondrial inhibition was likely responsible for steatosis, canonical pathway analysis revealed that progression to liver injury may be mediated by activation of the innate immune system primarily through NF-kB pathways. These results demonstrate the utility of a gene expression response compendium in developing transcriptional biomarkers and identifying the mechanistic basis for toxicity.

Challenges for using quantitative PCR test batteries as a TIE-type approach to identify metal exposure in benthic invertebrates

The epibenthic amphipod Melita plumulosa shows unique gene expression profiles when exposed to different contaminants. We hypothesized that specific changes in transcript abundance could be used in a battery of quantitative polymerase chain reaction (qPCR) assays as a toxicity identification evaluation (TIE)-like approach to identify the most relevant stressor in field-contaminated sediments. To test this hypothesis, seven candidate transcriptomic markers were selected, and their specificity following metal exposure was confirmed. The performance of these markers across different levels of added metals was verified. The ability of these transcripts to act as markers was tested by exposing amphipods to metal-contaminated field-collected sediments and measuring changes in transcript abundance via qPCR. For two of the three sediments tested, at least some of the transcriptomic patterns matched our predictions, suggesting that they would be effective in helping to identify metal exposure in field sediments. However, following exposure to the third sediment, transcriptomic patterns were unlike our predictions. These results suggest that the seven transcripts may be insufficient to discern individual contaminants from complex mixtures and that microarray or RNA-Seq global gene expression profiles may be more effective for TIE. Changes in transcriptomics based on laboratory exposures to single compounds should be carefully validated before the results are used to analyze mixtures.

Approaches and perspectives to toxicogenetics and toxicogenomics

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

Dissolved and particulate copper exposure induces differing gene expression profiles and mechanisms of toxicity in the deposit feeding amphipod Melita plumulosa

Uptake of metals via ingestion is an important route of exposure for many invertebrates, and it has been suggested that the toxic response to metals accumulated via food differs from that of metals accumulated via the dissolved phase. To test this hypothesis, the deposit-feeding epibenthic amphipod Melita plumulosa was exposed to nontoxic or reproductively toxic concentrations of copper via the overlying water, via ingestion of sediment, or via a combination of the two. Rates of copper uptake from the two exposure routes were predicted using a biokinetic model. Gene expression profiles were measured via microarray analysis and confirmed via quantitative polymerase chain reaction. Differences in expression profiles were related to the exposure route more than to individual or combined rates of copper uptake. Chitinase and digestive protease transcript expression levels correlated to the copper uptake rate from sediment, rather than from the dissolved phase or combined total uptake rate. Overall, this study supports the hypothesis that metals accumulated via ingestion have a different mode of toxic action than metals taken up from water. Consequently, guidelines that only consider dissolved metal exposure, including equilibrium-partitioning-based guidelines, may underestimate the potential effects from deposited or resuspended metal-contaminated sediments.

Assessing mechanisms of toxicant response in the amphipod Melita plumulosa through transcriptomic profiling

This study describes the function of transcripts with altered abundance in the epibenthic amphipod, Melita plumulosa, following whole-sediment exposure to a series of common environmental contaminants. M. plumulosa were exposed for 48 h to sediments spiked and equilibrated with the following contaminants at concentrations predicted to cause sublethal effects to reproduction: porewater ammonia 30 mg L(-1); bifenthrin at 100 μg kg(-1); fipronil at 50 μg kg(-1); 0.6% diesel; 0.3% crude oil; 250 mg Cu kg(-1); 400 mg Ni kg(-1); and 400 mg Zn kg(-1). RNA was extracted and hybridized against a custom Agilent microarray developed for this species. Although the microarray represented a partial transcriptome and not all features on the array could be annotated, unique transcriptomic profiles were generated for each of the contaminant exposures. Hierarchical clustering grouped the expression profiles together by contaminant class, with copper and zinc, the petroleum products and nickel, and the pesticides each forming a distinct cluster. Many of the transcriptional changes observed were consistent with patterns previously described in other crustaceans. The changes in the transcriptome demonstrated that contaminant exposure caused changes in digestive function, growth and moulting, and the cytoskeleton following metal exposure, whereas exposure to petroleum products caused changes in carbohydrate metabolism, xenobiotic metabolism and hormone cycling. Functional analysis of these gene expression profiles can provide a better understanding of modes of toxic action and permits the prediction of mixture effects within contaminated ecosystems.

Molecular interactions of nanomaterials and organisms: defining biomarkers for toxicity and high-throughput screening using traditional and next-generation sequencing approaches

The expression of molecular pathways in an organism provides a clue as to the potential impacts of exposure to nanomaterials.

Selection of reliable biomarkers from PCR array analyses using relative distance computational model: methodology and proof-of-concept study

It is increasingly evident about the difficulty to monitor chemical exposure through biomarkers as almost all the biomarkers so far proposed are not specific for any individual chemical. In this proof-of-concept study, adult male zebrafish (Danio rerio) were exposed to 5 or 25 µg/L 17β-estradiol (E2), 100 µg/L lindane, 5 nM 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) or 15 mg/L arsenic for 96 h, and the expression profiles of 59 genes involved in 7 pathways plus 2 well characterized biomarker genes, vtg1 (vitellogenin1) and cyp1a1 (cytochrome P450 1A1), were examined. Relative distance (RD) computational model was developed to screen favorable genes and generate appropriate gene sets for the differentiation of chemicals/concentrations selected. Our results demonstrated that the known biomarker genes were not always good candidates for the differentiation of pair of chemicals/concentrations, and other genes had higher potentials in some cases. Furthermore, the differentiation of 5 chemicals/concentrations examined were attainable using expression data of various gene sets, and the best combination was the set consisting of 50 genes; however, as few as two genes (e.g. vtg1 and hspa5 [heat shock protein 5]) were sufficient to differentiate the five chemical/concentration groups in the present test. These observations suggest that multi-parameter arrays should be more reliable for biomonitoring of chemical exposure than traditional biomarkers, and the RD computational model provides an effective tool for the selection of parameters and generation of parameter sets.

Gene expression analysis of so called asian dust extracts in human acute myeloid leukemia cells

As the frequency and the intensity of so called Asian dust (AD) events have increased, public concerns about the adverse health effects has spiked sharply over the last two decades. Despite the recent reports on the correlation between AD events and the risk for cardiovascular and respiratory disease, the nature of the toxicity and the degree of the risk are yet largely unknown. In the present study, we investigated the effects of the dichloromethane extract of AD (AD-X) and that of urban dust (NAD-X) collected during a non-AD period on gene expression in HL-60 cells using Illumina Sentrix HumanRef-8 Expression BeadChips. Global changes in gene expression were analyzed after 24 h of incubation with 50 or 100 μg/ml AD-X and NAD-X. By one-way analysis of variance (p < 0.05) and Benjamini-Hochberg multiple testing correction for false discovery rate of the results, 573 and 297 genes were identified as AD-X- and NAD-X-responsive, respectively. The genes were classified into three groups by Venn diagram analysis of their expression profile, i.e., 290 AD-X-specific, 14 NAD-X-specific, and 283 overlapping genes. Quantitative realtime PCR confirmed the changes in the expression levels of the selected genes. The expression patterns of five genes, namely SORL1, RABEPK, DDIT4, AZU1, and NUDT1 differed significantly between the two groups. Following rigorous validation process, these genes may provide information in developing biomarker for AD exposure.

Comparative genomic hybridization (CGH) in genotoxicology

In the past two decades comparative genomic hybridization (CGH) and array CGH have become crucial and indispensable tools in clinical diagnostics. Initially developed for the genome-wide screening of chromosomal imbalances in tumor cells, CGH as well as array CGH have also been employed in genotoxicology and most recently in toxicogenomics. The latter methodology allows a multi-endpoint analysis of how genes and proteins react to toxic agents revealing molecular mechanisms of toxicology. This chapter provides a background on the use of CGH and array CGH in the context of genotoxicology as well as a protocol for conventional CGH to understand the basic principles of CGH. Array CGH is still cost intensive and requires suitable analytical algorithms but might become the dominating assay in the future when more companies provide a large variety of different commercial DNA arrays/chips leading to lower costs for array CGH equipment as well as consumables such as DNA chips. As the amount of data generated with microarrays exponentially grows, the demand for powerful adaptive algorithms for analysis, competent databases, as well as a sound regulatory framework will also increase. Nevertheless, chromosomal and array CGH are being demonstrated to be effective tools for investigating copy number changes/variations in the whole genome, DNA expression patterns, as well as loss of heterozygosity after a genotoxic impact. This will lead to new insights into affected genes and the underlying structures of regulatory and signaling pathways in genotoxicology and could conclusively identify yet unknown harmful toxicants.

Comparison of toxicity and transcriptomic profiles in a diatom exposed to oil, dispersants, dispersed oil

Dispersants are commonly used to mitigate the impact of oil spills, however, the ecological cost associated with their use is uncertain. The toxicity of weathered oil, dispersed weathered oil, and the hydrocarbon-based dispersant Slickgone NS(®), to the diatom Phaeodactylum tricornutum has been examined using standardized toxicity tests. The assumption that most toxicity occurs via narcosis was tested by measuring membrane damage in diatoms after exposure to one of the petroleum products. The mode of toxic action was determined using microarray-based gene expression profiling in diatoms after exposure to one of the petroleum products. The diatoms were found to be much more sensitive to dispersants than to the water accommodated fraction (WAF), and more sensitive to the chemically enhanced WAF (CEWAF) than to either the WAF itself or the dispersants. Exposure to dispersants and CEWAF caused membrane damage, while exposure to WAF did not. The gene expression profiles resulting from exposure to all three petroleum mixtures were highly similar, suggesting a similar mode of action for these compounds. The observed toxicity bore no relationship to PAH concentrations in the water column or to total petroleum hydrocarbon (TPH), suggesting that an undescribed component of the oil was causing toxicity. Taken together, these results suggest that the use of dispersants to clean up oil spills will dramatically increase the oil toxicity to diatoms, and may have implications for ecological processes such as the timing of blooms necessary for recruitment.

mRNA levels in control rat liver display strain-specific, hereditary, and AHR-dependent components

Rat is a major model organism in toxicogenomics and pharmacogenomics. Hepatic mRNA profiles after treatment with xenobiotic chemicals are used to predict and understand drug toxicity and mechanisms. Surprisingly, neither inter- and intra-strain variability of mRNA abundances in control rats nor the heritability of rat mRNA abundances yet been established. We address these issues by studying five populations: the popular Sprague-Dawley strain, sub-strains of Long-Evans and Wistar rats, and two lines derived from crosses between the Long-Evans and Wistar sub-strains. Using three independent techniques--variance analysis, linear modelling, and unsupervised pattern recognition--we characterize extensive intra- and inter-strain variability in mRNA levels. We find that both sources of variability are non-random and are enriched for specific functional groups. Specific transcription-factor binding-sites are enriched in their promoter regions and these genes occur in "islands" scattered throughout the rat genome. Using the two lines generated by crossbreeding we tested heritability of hepatic mRNA levels: the majority of rat genes appear to exhibit directional genetics, with only a few interacting loci. Finally, a comparison of inter-strain heterogeneity between mouse and rat orthologs shows more heterogeneity in rats than mice; thus rat and mouse heterogeneity are uncorrelated. Our results establish that control hepatic mRNA levels are relatively homogeneous within rat strains but highly variable between strains. This variability may be related to increased activity of specific transcription-factors and has clear functional consequences. Future studies may take advantage of this phenomenon by surveying panels of rat strains.

Promise and progress in environmental genomics: a status report on the applications of gene expression-based microarray studies in ecologically relevant fish species

The advent of any new technology is typically met with great excitement. So it was a few years ago, when the combination of advances in sequencing technology and the development of microarray technology made measurements of global gene expression in ecologically relevant species possible. Many of the review papers published around that time promised that these new technologies would revolutionize environmental biology as they had revolutionized medicine and related fields. A few years have passed since these technological advancements have been made, and the use of microarray studies in non-model fish species has been adopted in many laboratories internationally. Has the relatively widespread adoption of this technology really revolutionized the fields of environmental biology, including ecotoxicology, aquaculture and ecology, as promised? Or have these studies merely become a novelty and a potential distraction for scientists addressing environmentally relevant questions? In this review, the promises made in early review papers, in particular about the advances that the use of microarrays would enable, are summarized; these claims are compared to the results of recent studies to determine whether the forecasted changes have materialized. Some applications, as discussed in the paper, have been realized and have led to advances in their field, others are still under development.

Mapping drug physico-chemical features to pathway activity reveals molecular networks linked to toxicity outcome

The identification of predictive biomarkers is at the core of modern toxicology. So far, a number of approaches have been proposed. These rely on statistical inference of toxicity response from either compound features (i.e., QSAR), in vitro cell based assays or molecular profiling of target tissues (i.e., expression profiling). Although these approaches have already shown the potential of predictive toxicology, we still do not have a systematic approach to model the interaction between chemical features, molecular networks and toxicity outcome. Here, we describe a computational strategy designed to address this important need. Its application to a model of renal tubular degeneration has revealed a link between physico-chemical features and signalling components controlling cell communication pathways, which in turn are differentially modulated in response to toxic chemicals. Overall, our findings are consistent with the existence of a general toxicity mechanism operating in synergy with more specific single-target based mode of actions (MOAs) and provide a general framework for the development of an integrative approach to predictive toxicology.

Involvement of multiple molecular events in pyrogallol-induced hepatotoxicity and silymarin-mediated protection: evidence from gene expression profiles

In this study, the involvement of various molecular events in pyrogallol-mediated hepatotoxicity was deciphered by differential mRNA transcription profiles of control and pyrogallol treated mice liver. The modulatory effects of silymarin on pyrogallol-induced differentially expressed transcripts were also looked into. Swiss albino mice were treated with or without pyrogallol. In some set of experiments, mice were also treated with silymarin 2 h prior to pyrogallol. Total RNA was isolated from liver and polyadenylated RNA was reverse-transcribed into Cye 3 or Cye 5 labelled cDNA. Equal amounts of labelled cDNA from two different groups were mixed and hybridized with mouse 15 k array. The hybridized arrays were scanned, analyzed and the expression level of each transcript was calculated. The differential expression was validated by quantitative real time polymerase chain reaction. Comparative transcription pattern showed an alteration in the expression of 183 transcripts (150 up-regulated and 33 down-regulated) associated with oxidative stress, cell cycle, cytoskeletal network, cell-cell adhesion, extra-cellular matrix, inflammation, apoptosis, cell-signaling and intermediary metabolism in pyrogallol-exposed liver and silymarin pre-treatment modulated the expression of many of these transcripts. Results obtained thus suggest that pyrogallol induces multiple molecular events leading to hepatotoxicity and silymarin effectively counteracts pyrogallol-mediated alterations.

The anti-tumour agent, cisplatin, and its clinically ineffective isomer, transplatin, produce unique gene expression profiles in human cells

Cisplatin is a DNA-damaging anti-cancer agent that is widely used to treat a range of tumour types. Despite its clinical success, cisplatin treatment is still associated with a number of dose-limiting toxic side effects. The purpose of this study was to clarify the molecular events that are important in the anti-tumour activity of cisplatin, using gene expression profiling techniques. Currently, our incomplete understanding of this drug's mechanism of action hinders the development of more efficient and less harmful cisplatin-based chemotherapeutics. In this study the effect of cisplatin on gene expression in human foreskin fibroblasts has been investigated using human 19K oligonucleotide microarrays. In addition its clinically inactive isomer, transplatin, was also tested. Dualfluor microarray experiments comparing treated and untreated cells were performed in quadruplicate. Cisplatin treatment was shown to significantly up- or down-regulate a consistent subset of genes. Many of these genes responded similarly to treatment with transplatin, the therapeutically inactive isomer of cisplatin. However, a smaller proportion of these transcripts underwent differential expression changes in response to the two isomers. Some of these genes may constitute part of the DNA damage response induced by cisplatin that is critical for its anti-tumour activity. Ultimately, the identification of gene expression responses unique to clinically active compounds, like cisplatin, could thus greatly benefit the design and development of improved chemotherapeutics.

Expression of copper-responsive genes in HepG2 cells

The hypothesis that copper modulates the activity of intracellular signal transduction pathways to affect transcription, which ultimately disrupts normal development was investigated. Preliminary analysis of transcriptomes from HepG2 cells exposed to copper for 4 and 24 h identified 19 and 7 up-regulated genes (twofold; p <or= 0.05), respectively. Among the up-regulated genes, several have been previously reported to be responsive to metals or oxidative stress. Differentially expressed genes were grouped by the functional categories based on gene ontology (GO). Significantly enriched GO categories (p < 0.01) included copper ion homeostasis, cadmium and copper ion binding, and heme oxygenase and oxidoreductase activities. Real-time RT-PCR confirmed the effect of copper on the levels of MT2A, HSPA1A, CYP1A1 and HMOX1 expression.

Temporal changes in gene expression in rainbow trout exposed to ethynyl estradiol

We examined changes in the genomic response during continuous exposure to the xenoestrogen ethynyl estradiol. Isogenic rainbow trout Oncorhynchus mykiss were exposed to nominal concentrations of 100 ng/L ethynyl estradiol (EE2) for a period of 3 weeks. At fixed time points within the exposure, fish were euthanized, livers harvested and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Salmonid array (GRASP project, University of Victoria, Canada) spotted with 16,000 cDNAs. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up and down regulated genes, and to determine gene clustering patterns that can be used as "expression signatures". Gene ontology was determined using the annotation available from the GRASP website. Our analysis indicates each exposure time period generated specific gene expression profiles. Changes in gene expression were best understood by grouping genes by their gene expression profiles rather than examining fold change at a particular time point. Many of the genes commonly used as biomarkers of exposure to xenoestrogens were not induced initially and did not have gene expression profiles typical of the majority of genes with altered expression.

A dolphin peripheral blood leukocyte cDNA microarray for studies of immune function and stress reactions

A microarray focused on stress response and immune function genes of the bottlenosed dolphin has been developed. Random expressed sequence tags (ESTs) were isolated and sequenced from two dolphin peripheral blood leukocyte (PBL) cDNA libraries biased towards T- and B-cell gene expression by stimulation with IL-2 and LPS, respectively. A total of 2784 clones were sequenced and contig analysis yielded 1343 unigenes (archived and annotated at ). In addition, 52 dolphin genes known to be important in innate and adaptive immune function and stress responses of terrestrial mammals were specifically targeted, cloned and added to the unigene collection. The set of dolphin sequences printed on a cDNA microarray comprised the 1343 unigenes, the 52 targeted genes and 2305 randomly selected (but unsequenced) EST clones. This set was printed in duplicate spots, side by side, and in two replicates per slide, such that the total number of features per microarray slide was 19,200, including controls. The dolphin arrays were validated and transcriptomic profiles were generated using PBL from a wild dolphin, a captive dolphin and dolphin skin cells. The results demonstrate that the array is a reproducible and informative tool for assessing differential gene expression in dolphin PBL and in other tissues.

Genotoxicants target distinct molecular networks in neonatal neurons

BACKGROUND

Exposure of the brain to environmental agents during critical periods of neuronal development is considered a key factor underlying many neurologic disorders.

OBJECTIVES

In this study we examined the influence of genotoxicants on cerebellar function during early development by measuring global gene expression changes.

METHODS

We measured global gene expression in immature cerebellar neurons (i.e., granule cells) after treatment with two distinct alkylating agents, methylazoxymethanol (MAM) and nitrogen mustard (HN2). Granule cell cultures were treated for 24 hr with MAM (10-1,000 microM) or HN2 (0.1-20 microM) and examined for cell viability, DNA damage, and markers of apoptosis.

RESULTS

Neuronal viability was significantly reduced (p < 0.01) at concentrations > 500 microM for MAM and > 1.0 microM for HN2; this correlated with an increase in both DNA damage and markers of apoptosis. Neuronal cultures treated with sublethal concentrations of MAM (100 microM) or HN2 (1.0 microM) were then examined for gene expression using large-scale mouse cDNA microarrays (27,648). Gene expression results revealed that a) global gene expression was predominantly up-regulated by both genotoxicants; b) the number of down-regulated genes was approximately 3-fold greater for HN2 than for MAM; and c) distinct classes of molecules were influenced by MAM (i.e, neuronal differentiation, the stress and immune response, and signal transduction) and HN2 (i.e, protein synthesis and apoptosis).

CONCLUSIONS

These studies demonstrate that individual genotoxicants induce distinct gene expression signatures. Further study of these molecular networks may explain the variable response of the developing brain to different types of environmental genotoxicants.

DNA microarrays in herbal drug research

Natural products are gaining increased applications in drug discovery and development. Being chemically diverse they are able to modulate several targets simultaneously in a complex system. Analysis of gene expression becomes necessary for better understanding of molecular mechanisms. Conventional strategies for expression profiling are optimized for single gene analysis. DNA microarrays serve as suitable high throughput tool for simultaneous analysis of multiple genes. Major practical applicability of DNA microarrays remains in DNA mutation and polymorphism analysis. This review highlights applications of DNA microarrays in pharmacodynamics, pharmacogenomics, toxicogenomics and quality control of herbal drugs and extracts.

Prediction of carcinogenic potential by a toxicogenomic approach using rat hepatoma cells

The long-term rodent bioassay is the standard method to predict the carcinogenic hazard of chemicals for humans. However, this assay is costly, and the results take at least two years to produce. In the present study, we conducted gene expression profiling of cultured cells exposed to carcinogenic chemicals with the aim of providing a basis for rapid and reliable prediction of carcinogenicity using microarray technology. We selected 39 chemicals, including 17 rat hepatocarcinogens and eight compounds demonstrating carcinogenicity in organs other than the liver. The remaining 14 were non-carcinogens. When rat hepatoma cells (MH1C1) were treated with the chemicals for 3 days at a non-toxic dose, analysis of gene expression changes with our in-house microarray allowed a set of genes to be identified differentiating hepatocarcinogens from non-carcinogens, and all carcinogens from non-carcinogens, by statistical methods. Moreover, optimization of the two gene sets for classification with an SVM and LOO-CV resulted in selection of 39 genes. The highest predictivity was achieved with 207 genes for differentiation between non-hepatocarcinogens and non-carcinogens. The overlap between the two selected gene sets encompassed 26 genes. This gene set contained significant genes for prediction of carcinogenicity, with a concordance of 84.6% by LOO-CV SVM. Using nine external samples, correct prediction of carcinogenicity by SVM was 88.9%. These results indicate that short-term bioassay systems for carcinogenicity using gene expression profiling in hepatoma cells have great promise.

NSAID-induced acute phase response is due to increased intestinal permeability and characterized by early and consistent alterations in hepatic gene expression

Toxicogenomics using a reference database can provide a better understanding and prediction of toxicity, largely by creating biomarkers that tie gene expression to actual pathology events. During the course of building a toxicogenomic database, an observation was made that a number of non-steroidal anti-inflammatory compounds (NSAIDs) at supra-pharmacologic doses induced an acute phase response (APR) and displayed hepatic gene expression patterns similar to that of intravenous lipopolysaccharide (LPS). Since NSAIDs are known to cause injury along the gastrointestinal tract, it has been suggested that NSAIDs increase intestinal permeability, allowing LPS and/or bacteria into the systemic circulation and stimulating an APR detectable in the liver. A short term study was subsequently conducted examining the effects of aspirin, indomethacin, ibuprofen, and rofecoxib to rats and a variety of endpoints were examined that included serum levels of inflammatory cytokines, histologic evaluation, and hepatic gene expression. Both indomethacin and ibuprofen injured the gastrointestinal tract, induced an APR, and increased serum levels of LPS, while rofecoxib and aspirin did not affect the GI tract or induce an APR. In treatments that eventually showed a systemic inflammatory response, hepatic expression of many inflammatory genes was noted as early as 6 hours after treatment well before alterations in traditional clinical pathology markers were detected. This finding led to the creation of a hepatic gene expression biomarker of APR that was effectively shown to be an early identifier of imminent inflammatory injury. In terms of the relative gastrointestinal safety and the NSAIDs studied, an important safety distinction can be made between the presumptive efficacious dose and the APR-inducing dose for indomethacin (1-2-fold), ibuprofen (5-fold), and rofecoxib (approximately 250-fold). Our data support the notion that NSAID-induced intestinal injury results in leakage of commensural bacteria and/or LPS into the circulation, provoking a systemic inflammatory response and that hepatic gene expression-based biomarkers can be used as early and sensitive biomarkers of APR onset.

Gene expression patterns in rainbow trout, Oncorhynchus mykiss, exposed to a suite of model toxicants

The increased availability and use of DNA microarrays has allowed the characterization of gene expression patterns associated with exposure to different toxicants. An important question is whether toxicant induced changes in gene expression in fish are sufficiently diverse to allow for identification of specific modes of action and/or specific contaminants. In theory, each class of toxicant may generate a gene expression profile unique to its mode of toxic action. In this study, isogenic (cloned) rainbow trout Oncorhynchus mykiss were exposed to sublethal levels of a series of model toxicants with varying modes of action, including ethynylestradiol (xeno-estrogen), 2,2,4,4'-tetrabromodiphenyl ether (BDE-47, thyroid active), diquat (oxidant stressor), chromium VI, and benzo[a]pyrene (BaP) for a period of 1-3 weeks. An additional experiment measured trenbolone (anabolic steroid; model androgen) induced gene expression changes in sexually mature female trout. Following exposure, fish were euthanized, livers removed and RNA extracted. Fluorescently labeled cDNA were generated and hybridized against a commercially available Atlantic Salmon/Trout array (GRASP project, University of Victoria) spotted with 16,000 cDNA's. The slides were scanned to measure abundance of a given transcript in each sample relative to controls. Data were analyzed via Genespring (Silicon Genetics) to identify a list of up- and downregulated genes, as well as to determine gene clustering patterns that can be used as "expression signatures". The results indicate each toxicant exposure caused between 64 and 222 genes to be significantly altered in expression. Most genes exhibiting altered expression responded to only one of the toxicants and relatively few were co-expressed in multiple treatments. For example, BaP and Diquat, both of which exert toxicity via oxidative stress, upregulated 28 of the same genes, of over 100 genes altered by either treatment. Other genes associated with steroidogenesis, p450 and estrogen responsive genes appear to be useful for selectively identifying toxicant mode of action in fish, suggesting a link between gene expression profile and mode of toxicity. Our array results showed good agreement with quantitative real time polymerase chain reaction (qRT PCR), which demonstrates that the arrays are an accurate measure of gene expression. The specificity of the gene expression profile in response to a model toxicant, the link between genes with altered expression and mode of toxic action, and the consistency between array and qRT PCR results all suggest that cDNA microarrays have the potential to screen environmental contaminants for biomarkers and mode of toxic action.

Human placenta: a human organ for developmental toxicology research and biomonitoring

Pregnant mothers are exposed to a wide variety of foreign chemicals. This exposure is most commonly due to maternal medication, lifestyle factors, such as smoking, drug abuse, and alcohol consumption, or occupational and environmental sources. Foreign compounds may interfere with placental functions at many levels e.g. signaling, production and release of hormones and enzymes, transport of nutrients and waste products, implantation, cellular growth and maturation, and finally, at the terminal phase of placental life, i.e. delivery. Placental responses may also be due to pharmaco-/toxicodynamic responses to foreign chemicals, e.g. hypoxia. On the other hand, placental xenobiotic-metabolizing enzymes can detoxify or activate foreign chemicals, and transporters either enhance or prevent cellular accumulation and transfer across the placenta. The understanding of what xenobiotics do to the placenta and what the placenta does to the xenobiotics should provide the basis for the use of placenta as a tool to investigate and predict some aspects of developmental toxicity. This review aims to give an update of the fate and behavior of xenobiotics in the placenta from the viewpoint of xenobiotic-metabolizing enzymes and transporters. Their response levels will be described according to gestational status and methods used. The effects of foreign chemicals on placental metabolizing enzymes will be discussed. Also, interactions in the transporter protein level will be covered. The role of the placenta in contributing to developmental effects and fetotoxicity will be examined. The toxicological effects of maternal medications, smoking, and environmental exposures (dioxins, pesticides) as well as some possibilities for biomonitoring will be highlighted.

Differentiating mechanisms of toxicity using global gene expression analysis in Saccharomyces cerevisiae

Genotoxic stress triggers a variety of biological responses including the transcriptional activation of genes regulating DNA repair, cell survival and cell death. Genomic approaches, which monitor gene expressions across large numbers of genes, can serve as a powerful tool for exploring mechanisms of toxicity. Here, using five different agents, we investigated whether the analysis of genome-wide expression profiles in Saccharomyces cerevisiae could provide insights into mechanisms of genotoxicity versus cytotoxicity. To differentiate the genotoxic stress-associated expression signatures from that of a general cytotoxic stress, we compared gene expression profiles following the treatment with DNA-reactive (cisplatin, MMS, bleomycin) and DNA non-reactive (ethanol and sodium chloride) compounds. Although each of the tested chemicals produced a distinct gene expression profile, we were able to identify a gene expression signature consisting of a relatively small number of biologically relevant genes capable of differentiating genotoxic and cytotoxic stress. The gene set includes such upregulated genes as HUG1, ECM4 and previously uncharacterized gene, YLR297W in the genotoxic and GAP1, CGR1 in the cytotoxic group. Our results indicate the potential of gene expression profile analysis for elucidating mechanism of action of genotoxic agents.

Development and application of a brain-specific cDNA microarray for effect evaluation of neuro-active pharmaceuticals in zebrafish (Danio rerio)

The environmental fate and ecotoxicological effect of pharmaceuticals are poorly understood, and standardized tests to detect and evaluate their potential effects in the environment are not available. We developed a zebrafish brain-specific microarray containing 682 neurologically relevant cDNA-fragments. To investigate the applicability of this microarray for studying neurotoxic modes-of-action and impact assessment of neuro-active pharmaceuticals in zebrafish, chlorpromazine was used as a model compound. After exposure to chlorpromazine (75 microg/L) for 2, 4, 14 and 28 days or control treatment RNA was extracted from brains of males and females. Fluorescently labeled cDNA was prepared and hybridized to the custom microarray. In total, 56 genes were differentially expressed in brains of male and/or female zebrafish, of which most genes were down-regulated. A clear difference in response to chlorpromazine exposure between males and females was observed with exposure time as well as in functional classes of affected genes. The presented study is one of the first reports on molecular effects of human neuro-active pharmaceuticals in aquatic non-target organisms. This new genomic tool successfully detected gene expression effects of exposure to chlorpromazine in the brain of zebrafish. Reported gene expression effects are found to be consistent with literature data for other laboratory animals.

Using genome-wide transcriptional profiling to elucidate small-molecule mechanism

Transcriptional profiling with DNA microarrays can be used to measure the genome-wide transcriptional response to small molecules. Recent progress in the analysis of gene-expression data has relied on the generation of databases of profiles documenting the transcriptional effects of various compound treatments and genetic perturbations. A positive correlation between the transcriptional response induced by a novel small molecule and a database profile can provide insight into the molecule's mechanism. Transcriptional profiling can also be used to assess a small molecule's specificity for its target and to facilitate analysis of pathways downstream of the target.

Ecotoxicogenomics: linkages between exposure and effects in assessing risks of aquatic contaminants to fish

Understanding the biological effects of exposures to chemicals in the environment relies on classical methods and emerging technologies in the areas of genomics, proteomics, and metabonomics. Linkages between the historical and newer toxicological tools are currently being developed in order to predict and assess risk. Being able to classify chemicals and other stressors based on effects they have at the molecular, tissue, and organismal levels helps define a systems biology approach to development of streamlined, cost-effective, and comprehensive testing approaches for evaluating environmental hazards. The challenges of the individual technologies and the combinations of tools for ecotoxicogenomics are discussed in application to aquatic toxicology with a particular emphasis on fish testing.

Gene expression profiling in chronic copper overload reveals upregulation ofPrnpandApp

The level at which copper becomes toxic is not clear. Several studies have indicated that copper causes oxidative stress; however, most have tested very high levels of copper exposure. We currently have only a limited understanding of the protective systems that operate in cells chronically exposed to copper. Additionally, the limits of homeostatic regulation are not known, making it difficult to define the milder effects of copper excess. Furthermore, a robust assay to facilitate the diagnosis of copper excess and to distinguish mild, moderate, and severe copper overload is needed. To address these issues, we have investigated the effects on steady-state gene expression of chronic copper overload in a cell culture model system using cDNA microarrays. For this study we utilized cells from genetic models of copper overload: fibroblast cells from two mouse mutants, C57BL/6- Atp7aMobrand C57BL/6- Atp7aModap. These cell lines accumulate copper to abnormally high levels in normal culture media due to a defect in copper export from the cell. We identified 12 differentially expressed genes in common using our outlier identification methods. Surprisingly, our results show no evidence of oxidative stress in the copper-loaded cells. In addition, candidate components perhaps responsible for a copper-specific homeostatic response are identified. The genes that encode for the prion protein and the amyloid-β precursor protein, two known copper-binding proteins, are upregulated in both cell lines.

Toxicogenomics in predictive toxicology in drug development

The goal of toxicology is the assessment of possible risk to man. An emerging technology with the potential to have a major impact on risk assessment is toxicogenomics. In this review, we provide an overview of the many possibilities for toxicogenomics including technology platforms, data interpretation, and regulatory perspective and we give examples of toxicogenomics investigations. Toxicogenomics is a powerful tool for compound classification, for mechanistic studies, and for the detection of toxicity markers. Thus, toxicogenomics helps in the extrapolation of findings across species and increases predictability. Biomarkers are valuable in the evaluation of compounds at earlier development phases, improving clinical candidate selection. Caution regarding the interpretation of the results is still necessary. Nevertheless, toxicogenomics will accelerate preclinical safety assessments and improve the prediction of toxic liabilities, as well as of potential risk accumulation for drug-drug or drug-disease interactions.

Toxicogenomics and systems toxicology: aims and prospects

Toxicogenomics combines transcript, protein and metabolite profiling with conventional toxicology to investigate the interaction between genes and environmental stress in disease causation. The patterns of altered molecular expression that are caused by specific exposures or disease outcomes have revealed how several toxicants act and cause disease. Despite these success stories, the field faces noteworthy challenges in discriminating the molecular basis of toxicity. We argue that toxicology is gradually evolving into a systems toxicology that will eventually allow us to describe all the toxicological interactions that occur within a living system under stress and use our knowledge of toxicogenomic responses in one species to predict the modes-of-action of similar agents in other species.

Temporal changes in gene expression in the liver of male plaice (Pleuronectes platessa) in response to exposure to ethynyl oestradiol analysed by macroarray and Real-Time PCR

Suppression subtractive hybridisation (SSH) was used to generate cDNA libraries representing genes differentially-expressed in liver from male plaice (Pleuronectes platessa) exposed to ethynyl oestradiol (EE2). BLAST analysis and alignments of the clones with database sequence suggested at least three vitellogenin (VTG) genes and three zona radiata protein (ZRP) genes were represented. Clones with unique sequence (62 up-, 13 down-regulated) were arrayed as probes on nylon membranes to investigate temporal expression of oestrogen-responsive genes in experimental animals. Arrays were hybridised with radiolabelled cDNAs prepared from hepatic mRNA from animals treated with EE2 for various times upto 21 days and from treated animals transferred to clean water for upto a further 31 days. By day 21 of treatment 11 out of 17 probes from unidentified genes, 21/22 VTG, 13/14 ZRP, 2/2 liver aspartic proteinase (LAP) and 8/10 other gene sequences were induced by EE2 exposure. Of the down-regulated sequences, only three showed significant, decreased expression and these encode cytochrome b and two with cryptic functions. Based on the pattern of temporal response the up-regulated probes fell into two classes. Pattern A reached maximum expression by day 16 of exposure and then declined prior to removal of EE2 at 21 days. Pattern B genes reached maximal expression between day 16 and 22, declining only after removal of EE2. Independent investigation of the expression patterns of selected probes using quantitative Real-Time PCR reproduced the distinctive patterns. The results indicate a previously unrecognised mechanism for oestrogenic toxicity in which there is a selective down-regulation of some egg proteins, potentially diminishing the quality of eggs and this may contribute to reproductive failure described elsewhere.

Functional genomics of UV radiation responses in human cells

The gene expression responses of MCF-7, a p53 wild-type (wt) human cell line, were monitored by cDNA microarray hybridization after exposure to different wavelengths of UV irradiation. Equitoxic doses of UVA, UVB, and UVC radiation were used to reduce survival to 37%. The effects of suramin, a signal pathway inhibitor, on the gene expression responses to the three UV wavelengths were also compared in this model system. UVB radiation triggered the broadest gene expression responses, and 172 genes were found to be consistently responsive in at least two-thirds of independent UVB experiments. These UVB radiation-responsive genes encode proteins with diverse cellular roles including cell cycle control, DNA repair, signaling, transcription, protein synthesis, protein degradation, and RNA metabolism. The set of UVB-responsive genes included most of the genes responding to an equitoxic dose of UVC radiation, plus additional genes that were not strongly triggered by UVC radiation. There was also some overlap with genes responding to an equitoxic dose of UVA radiation, although responses to this lower energy UV radiation were overall weaker. Signaling through growth factor receptors and other cytokine receptors was shown to have a major role in mediating UV radiation stress responses, as suramin, which inhibits such receptors, attenuated responses to UV radiation in nearly all the cases. Inhibition by suramin was greater for UVC than for UVB irradiation. This probably reflects the more prominent role in UVB damage response of signaling by reactive oxygen species, which would not be affected by suramin. Our results with suramin demonstrate the power of cDNA microarray hybridization to illuminate the global effects of a pharmacologic inhibitor on cell signaling.

Response of rainbow trout transcriptome to model chemical contaminants

We used high-density cDNA microarray in studies of responses of rainbow trout fry at sublethal ranges of beta-naphthoflavone, cadmium, carbon tetrachloride, and pyrene. The differentially expressed genes were grouped by the functional categories of Gene Ontology. Significantly different response to the studied compounds was shown by a number of classes, such as cell cycle, apoptosis, signal transduction, oxidative stress, subcellular and extracellular structures, protein biosynthesis, and modification. Cluster analysis separated responses to the contaminants at low and medium doses, whereas at high levels the adaptive reactions were masked with general unspecific response to toxicity. We found enhanced expression of many mitochondrial proteins as well as genes involved in metabolism of metal ions and protein biosynthesis. In parallel, genes related to stress and immune response, signal transduction, and nucleotide metabolism were down-regulated. We performed computer-assisted analyses of Medline abstracts retrieved for each compound, which helped us to indicate the expected and novel findings.

In situ synthesis of oligonucleotide microarrays

This contribution presents a brief overall look of the methods for the preparation of various types of DNA microarrays and a thorough examination of the methods for in situ synthesis of oligonucleotide microarrays.

Effect of serum cholesterol on the mRNA content of amyloid precursor protein in rat livers

Genes that showed mRNA content profiles, which correlated with serum concentrations of total cholesterol (T.CHO), were screened from the microarray data of phenobarbital (PB)- or clofibrate (CLO)-treated rat livers, and the correlation was evaluated based on Spearman's correlation coefficient. Many genes involved in the cholesterol or bile acid metabolism were highly correlated such as UDP-glucuronosyltransferase-21, apolipoprotein A-I and cMOAT. The mRNA content of the amyloid precursor protein (APP) showed the 5th highest correlation among the 8799 probes in the Affymetrix Rat Genome U34 Array. In the livers of rats fed a high-cholesterol (1%) diet for 33 days, serum T.CHO levels increased by 4.6-fold, and the hepatic APP mRNA content also increased by 1.9-fold compared to the control group. These data suggest that the hepatic APP mRNA content was affected by serum T.CHO, and that hepatic APP was involved in cholesterol metabolism in rat livers.

The pain of antisense: in vivo application of antisense oligonucleotides for functional genomics in pain and analgesia

As the genomic revolution continues to evolve, there is an increasing demand for efficient and reliable tools for functional characterization of individual gene products. Antisense oligonucleotide-mediated knockdown has been used successfully as a functional genomics tool in animal models of pain and analgesia yet skepticism regarding the validity and utility of antisense technology remains. Contributing to this uncertainty are the lack of systematic studies exploring antisense oligonucleotide use in vivo and the many technical and methodological challenges intrinsic to the method. This article reviews the contributions of antisense oligonucleotide-based studies to the field of pain and analgesia and the general principles of antisense technology. A special emphasis is placed on technical issues surrounding the successful application of antisense oligonucleotides in vivo, including sequence selection, antisense oligonucleotide chemistry, DNA controls, route of administration, uptake, dose-dependence, time-course and adequate evaluation of knockdown.

What is the relevance of bioinformatics to pharmacology?

Although bioinformatics achieved prominence because of its central role in genome data storage, management and analysis, its focus has shifted as the life sciences exploit these data. In pharmacology, genomic, transcriptomic and proteomic data are being used in the quest for drugs that fulfill unmet medical needs, are disease modifying or curative and are more effective and safer than current drugs. Bioinformatics is used in drug target identification and validation and in the development of biomarkers and toxicogenomic and pharmacogenomic tools to maximize the therapeutic benefit of drugs. Now that the 'parts list' of cellular signalling pathways is available, integrated computational and experimental programmes are being developed, with the goal of enabling in silico pharmacology by linking the genome, transcriptome and proteome to cellular pathophysiology.