BACKGROUND:Recently, microarray data analyses using functional pathway information, e.g., gene set enrichment analysis (GSEA) and significance analysis of function and expression (SAFE), have gained recognition as a way to identify biological pathways/processes associated with a phenotypic endpoint. In these analyses, a local statistic is used to assess the association between the expression level of a gene and the value of a phenotypic endpoint. Then these gene-specific local statistics are combined to evaluate association for pre-selected sets of genes. Commonly used local statistics include t-statistics for binary phenotypes and correlation coefficients that assume a linear or monotone relationship between a continuous phenotype and gene expression level. Methods applicable to continuous non-monotone relationships are needed. Furthermore, for multiple experimental categories, methods that combine multiple GSEA/SAFE analyses are needed. RESULTS:For continuous or ordinal phenotypic outcome, we propose to use as the local statistic the coefficient of multiple determination (i.e., the square of multiple correlation coefficient) R2 from fitting natural cubic spline models to the phenotype-expression relationship. Next, we incorporate this association measure into the GSEA/SAFE framework to identify significant gene sets. Unsigned local statistics, signed global statistics and one-sided p-values are used to reflect our inferential interest. Furthermore, we describe a procedure for inference across multiple GSEA/SAFE analyses. We illustrate our approach using gene expression and liver injury data from liver and blood samples from rats treated with eight hepatotoxicants under multiple time and dose combinations. We set out to identify biological pathways/processes associated with liver injury as manifested by increased blood levels of alanine transaminase in common for most of the eight compounds. Potential statistical dependency resulting from the experimental design is addressed in permutation based hypothesis testing. CONCLUSION:The proposed framework captures both linear and non-linear association between gene expression level and a phenotypic endpoint and thus can be viewed as extending the current GSEA/SAFE methodology. The framework for combining results from multiple GSEA/SAFE analyses is flexible to address practical inference interests. Our methods can be applied to microarray data with continuous phenotypes with multi-level design or the meta-analysis of multiple microarray data sets.

This report details the standardized experimental design and the different data streams that were collected (histopathology, clinical chemistry, hematology and gene expression from the target tissue (liver) and a bio-available tissue (blood)) after treatment with eight known hepatotoxicants (at multiple time points and doses with multiple biological replicates). The results of the study demonstrate the classification of histopathological differences, likely reflecting differences in mechanisms of cell-specific toxicity, using either liver tissue or blood transcriptomic data.

Background Some of the biochemical events that lead to necrosis of the liver are well-known. However, the pathogenesis of necrosis of the liver from exposure to hepatotoxicants is a complex biological response to the injury. We hypothesize that gene expression profiles can serve as a signature to predict the level of necrosis elicited by acute exposure of rats to a variety of hepatotoxicants and postulate that the expression profiles of the predictor genes in the signature can provide insight to some of the biological processes and molecular pathways that may be involved in the manifestation of necrosis of the rat liver. Results Rats were treated individually with one of seven known hepatotoxicants and were analyzed for gene expression by microarray. Liver samples were grouped by the level of necrosis exhibited in the tissue. Analysis of significantly differentially expressed genes between adjacent necrosis levels revealed that inflammation follows programmed cell death in response to the agents. Using a Random Forest classifier with feature selection, 21 informative genes were identified which achieved 90%, 80% and 60% prediction accuracies of necrosis against independent test data derived from the livers of rats exposed to acetaminophen, carbon tetrachloride, and allyl alcohol, respectively. Pathway and gene network analyses of the genes in the signature revealed several gene interactions suggestive of apoptosis as a process possibly involved in the manifestation of necrosis of the liver from exposure to the hepatotoxicants. Cytotoxic effects of TNF-α, as well as transcriptional regulation by JUN and TP53, and apoptosis-related genes possibly lead to necrosis. Conclusion The data analysis, gene selection and prediction approaches permitted grouping of the classes of rat liver samples exhibiting necrosis to improve the accuracy of predicting the level of necrosis as a phenotypic end-point observed from the exposure. The strategy, along with pathway analysis and gene network reconstruction, led to the identification of 1) expression profiles of genes as a signature of necrosis and 2) perturbed regulatory processes that exhibited biological relevance to the manifestation of necrosis from exposure of rat livers to the compendium of hepatotoxicants.

The S checkpoint response to ultraviolet radiation (UVC) that inhibits replicon initiation is dependent on the ATR and Chk1 kinases. Downstream effectors of this response, however, are not well characterized. Data reported here eliminated Cdc25A degradation and inhibition of Cdk2-cyclin E as intrinsic components of the UVC-induced pathway of inhibition of replicon initiation in human cells. A sublethal dose of UVC (1 J/m(2)), which selectively inhibits replicon initiation by 50%, failed to reduce the amount of Cdc25A protein or decrease Cdk2-cyclin E kinase activity. Cdc25A degradation was observed after irradiation with cytotoxic fluences of UVC, suggesting that severe inhibition of DNA chain elongation and activation of the replication checkpoint might be responsible for the UVC-induced degradation of Cdc25A. Another proposed effector of the S checkpoint is the Cdc7-Dbf4 complex. Dbf4 interacted weakly with Chk1 in vivo but was recognized as a substrate for Chk1-dependent phosphorylation in vitro. FLAG-Dbf4 formed complexes with endogenous Cdc7, and this interaction was stable in UVC-irradiated HeLa cells. Overexpression of FLAG- or Myc-tagged Dbf4 abrogated the S checkpoint response to UVC but not ionizing radiation. These findings implicate a Dbf4-dependent kinase as a possible target of the ATR- and Chk1-dependent S checkpoint response to UVC.

Liver diseases that induce nonuniform lesions often give rise to greatly varying histopathology results in needle biopsy samples from the same patient. This study examines whether gene expression analysis of such biopsies could provide a more representative picture of the overall condition of the liver. We utilized acetaminophen (APAP) as a model hepatotoxicant that gives a multifocal pattern of necrosis following toxic doses. Rats were treated with a single toxic or subtoxic dose of APAP and sacrificed 6, 24, or 48 hours after exposure. Left liver lobes were harvested, and both gene expression and histopathological analysis were performed on biopsy-sized samples. While histopathological evaluation of such small samples revealed significant sample to sample differences after toxic doses of APAP, gene expression analysis provided a very homogeneous picture and allowed clear distinction between subtoxic and toxic doses. The main biological function differentiating animals that received sub-toxic from those that had received toxic doses was an acute stress response at 6 hours and signs of energy depletion at later time points. Our results suggest that the use of genomic analysis of biopsy samples together with histopathological analysis could provide a more precise representation of the overall condition of a patient’s liver than histopathological evaluation alone.

The molecular basis linking folate deficiency to certain health conditions and developmental defects is not fully understood. We examined the consequences of folate deficiency on global gene expression by microarray and compared transcript levels in normal human fibroblast cells (GM03349) grown in folate-deficient and -sufficient medium. The largest represented groups from the selected genes functioned in cell signaling, the cytoskeleton and the extracellular matrix and included the Wnt pathway genes DKK1, WISP1 and WNT5A. Twelve selected genes were further validated by qRT-PCR. Analysis of six genes at 4, 7, 10 and 14 days indicated that the relative differences in transcript levels between folate-sufficient and -deficient cells increases with time. Transcripts for 7 of the 12 selected genes were detected in the human lymphoblast cell line GM02257, and of these, changes in 4 genes corresponded to the results with fibroblast cells. Fibroblast cells were treated with the compounds homocysteine, methotrexate and the MEK1/2 inhibitor U0126, and relative transcript levels of six genes were determined. U0126 caused changes that more closely mimicked those detected in folate-deficient cells. The response of the DKK1 and TAGLN gene promoters to folate deficiency and compounds was examined in NIH3T3 cells using luciferase reporter plasmids. Promoter activity for both genes was decreased by folate deficiency and methotrexate and unaffected by homocysteine. U0126 caused a decrease in DKK1 promoter activity at 50 μM and had no effect on TAGLN promoter activity. These findings suggest an alternative mechanism for how folate deficiency leads to changes in gene expression and altered cell function.

Gene expression profiling, metabolomic screens, and other high-dimensional methods have become an integral part of many biological investigations. To facilitate interpretation of these data, it is important to have detailed phenotypic data--including histopathology--to which these data can be associated, or anchored. However, as the amount of phenotypic data increases, associations within and across these data can be difficult to visualize and interpret. We have developed an approach for categorizing and clustering biologically related histopathological diagnoses to facilitate their visualization, thereby increasing the possibility of identifying associations and facilitating the comparison with other data streams. In this study, we utilize histopathological data generated as part of a standardized toxicogenomics compendium study to generate composite histopathological scores and to develop visualizations that facilitate biological insight. The validity of this approach is illustrated by the identification of transcripts that correlate with the pathology diagnoses that comprise the categories of "response to hepatocellular injury" and "repair." This approach is broadly applicable to studies in which histopathology is used to phenotypically anchor other data, and results in visualizations that facilitate biological interpretation and the identification of associations and relationships within the data.

Toxicogenomics provides the ability to examine in greater detail the underlying molecular events that precede and accompany toxicity, thus allowing prediction of adverse events at much earlier times compared to classical toxicological end points. Acetaminophen (APAP) is a pharmaceutical that has similar metabolic and toxic responses in rodents and humans. Recent gene expression profiling studies with APAP found an oxidative stress signature at a subtoxic dose that we hypothesized can be phenotypically anchored to conventional biomarkers of oxidative stress. Liver tissue was obtained from experimental animals used to generate microarray data, where male rats were given APAP at subtoxic (150 mg/kg) or overtly toxic (1500 and 2000 mg/kg) doses and sacrificed at 6, 24, or 48 h. Oxidative stress in liver was evaluated by a diverse panel of markers that included assessing expression of base excision repair (BER) genes, quantifying oxidative lesions in genomic DNA, and evaluating protein and lipid oxidation. A subtoxic dose of APAP produced significant accumulation of nitrotyrosine protein adducts. Both subtoxic and toxic doses caused a significant increase in 8-hydroxy-deoxyguanosine (8-OH-dG) as well as a significant decrease in glutathione (GSH) content. Only toxic doses of APAP significantly induced expression levels of BER genes. None of the doses examined resulted in a significant increase in the number of abasic sites or in the amount of lipid peroxidation. The accumulation of nitrotyrosine and 8-OH-dG adducts along with reduced GSH content in the liver phenotypically anchors the oxidative stress gene expression signature observed with a subtoxic dose of APAP, lending support to the validity of gene expression studies as a sensitive and biologically meaningful end point in toxicology.

Different degrees of a toxic response between and within the various lobes of the liver have been observed in rodents following treatment with acetaminophen. This study was designed to compare 2 sampling methods of the rat liver for gene-expression analysis. Ten male Fischer 344/N rats, 12–14 weeks of age, were treated with vehicle (0.5% aqueous ethyl cellulose) or acetaminophen (APAP, 1500 mg/kg) and sacrificed 24 hours following dose administration. Two representative sections were collected from the left liver lobe, stained with hematoxylin and eosin (H&E), and evaluated independently by 2 pathologists. The central core of the left lobe was cubed and frozen. Five random cubes were conserved, while the remaining left lobe core was pulverized. From each of the 10 animals, 2 random cubes and 2 samples from the homogeneous, pulverized samples were prepared for microarray analysis. Histopathologic evaluation revealed a variable response of centrilobular necrosis within the left lobe. Multiple methods used to analyze the microarray data indicated that sampling technique was not a major contributor to the variability observed in the gene expression data; however, only the powdered samples clustered for all animals, even those with disparate histopathologic results. Additionally, a powdered sample provided the advantages of a homogenous sample pool and the ability to use sample aliquots for other analyses to include proteomics, metabonomics, and other molecular techniques.

A critical component in the design of the Chemical Effects in Biological Systems (CEBS) Knowledgebase is a strategy to capture toxicogenomics study protocols and the toxicity endpoint data (clinical pathology and histopathology). A Study is generally an experiment carried out during a period of time for the purpose of obtaining data, and the Study Design Description captures the methods, timing, and organization of the Study. The CEBS Data Dictionary (CEBS-DD) has been designed to define and organize terms in an attempt to standardize nomenclature needed to describe a toxicogenomics Study in a structured yet intuitive format and provide a flexible means to describe a Study as conceptualized by the investigator. The CEBS-DD will organize and annotate information from a variety of sources, thereby facilitating the capture and display of toxicogenomics data in biological context in CEBS, i.e., associating molecular events detected in highly-parallel data with the toxicology/pathology phenotype as observed in the individual Study Subjects and linked to the experimental treatments. The CEBS-DD has been developed with a focus on acute toxicity studies, but with a design that will permit it to be extended to other areas of toxicology and biology with the addition of domain-specific terms. To illustrate the utility of the CEBS-DD, we present an example of integrating data from two proteomics and transcriptomics studies of the response to acute acetaminophen toxicity (A. N. Heinloth et al., 2004, Toxicol. Sci. 80, 193-202).

Oxidative stress and reactive oxygen species play a major role in both normal and pathophysiologic cellular processes. Although many cellular constituents can be damaged by oxidant exposure, cysteine thiol groups are among the most readily oxidized moieties found within cells. To avoid potentially irreversible cysteine thiol oxidation, cells have developed multiple antioxidant defenses to preserve these moieties. Among these defenses, protein S-glutathiolation has emerged as an important mechanism, both in the maintenance of thiol stability during oxidant exposure and as a rapid and efficient mechanism regulating protein activity and cellular metabolic pathways. Here we review the known molecular targets of S-glutathiolation, with emphasis on the varying molecular effects of S-glutathiolation on different proteins.

Toxicology and the practice of pathology are rapidly evolving in the postgenomic era. Observable treatment related changes have been the hallmark of toxicology studies. Toxicogenomics is a powerful new tool that may show gene and protein changes earlier and at treatment levels below the limits of detection of traditional measures of toxicity. It may also aid in the understanding of toxic mechanisms. It is important to remember that it is only a tool and will provide meaningful results only when properly applied. As is often the case with new experimental tools, the initial utilization is driven more by the technology than application to problem solving. Toxicogenomics is interdisciplinary in nature including at a minimum, pathology, toxicology, and genomics. Most studies will require the input from the disciplines of toxicology, pathology, molecular biology, bioinformatics, biochemistry, and others depending on the types of questions being asked.

Background: Oxidatively modified low density lipoprotein (OxLDL) induces apoptosis in vascular cells including macrophages, while NO exerts antiapoptotic effects. Here we studied the impact of nitric oxide (NO) on OxLDL-induced cytochrome c release, apoptosis, and expression of the proapoptotic p53 in macrophages. Methods: Human LDL was oxidised by Cu2+, and monocytes were prepared from human buffy coats. Differentiation to macrophages was achieved by culturing cells in the presence of human serum and was followed by detecting monocyte chemoattractant protein 1 (MCP-1) expression (RT-PCR). Cytochrome c release and p53 expression of macrophages were detected by immunoblotting, and apoptosis by visualisation of nuclear condensation. Results: OxLDL dose-dependently (50–200 μg/ml) induced cytochrome c release that was prevented by preincubation with the NO-donor S-nitrosoglutathione (GSNO) (100 μM) or with the cGMP analogue 8-br-cGMP (100 μM) for 15 h. In cells co-treated with GSNO and the soluble guanylate cyclase (sGC) inhibitor oxadialoquinoxalione (ODQ, 10 μM, 15 h), OxLDL-evoked cytochrome c release remained effective, indicating that NO acted via sGC-dependent cGMP formation. Parallel incubation of macrophages with 8-br-cGMP (100 μM) and ODQ (10 μM) for 15 h left the protective effect of 8-br-cGMP unaltered. Short pre-incubation (30 min) with GSNO or 8-br-cGMP was ineffective in preventing OxLDL-elicited cytochrome c release. Initiation of cytochrome c release in macrophages was paralleled by a dose-dependent accumulation of the proapoptotic factor p53, and by enhanced rate of nuclear condensation. Stabilisation of p53 was prevented by preincubation with the NO-donor GSNO or 8-br-cGMP, thus implying a downmodulatory effect of cGMP on pathways that upregulate the tumor suppressor p53. Conclusions: OxLDL induces cytochrome c release and apoptosis in human macrophages in close association with p53 accumulation. NO attenuates OxLDL-induced cytochrome c release and p53 accumulation via activation of sGC and cGMP formation. These effects may be of particular importance in arterial tissue with reduced NO activity.

We measured oxidative stress caused by lipoproteins in human umbilical vein endothelial cells and rabbit aorta. Oxidized lipoprotein(a) was the more potent stimulus over oxidized low density lipoprotein, and we believe that both influence the pathogenesis of atherosclerosis.