Search PubMed⌕ Search

Biomedical subjects

Richard S Paules

Publications and source records attributed to Richard S Paules.

At least 19 recordsLinked to original sources

Phenotypic anchoring of acetaminophen-induced oxidative stress with gene expression profiles in rat liver.

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.

Acetaminophen↗

Major carcinogenic pathways identified by gene expression analysis of peritoneal mesotheliomas following chemical treatment in F344 rats.

This study was performed to characterize the gene expression profile and to identify the major carcinogenic pathways involved in rat peritoneal mesothelioma (RPM) formation following treatment of Fischer 344 rats with o-nitrotoluene (o-NT) or bromochloracetic acid (BCA). Oligo arrays, with over 20,000 target genes, were used to evaluate o-NT- and BCA-induced RPMs, when compared to a non-transformed mesothelial cell line (Fred-PE). Analysis using Ingenuity Pathway Analysis software revealed 169 cancer-related genes that were categorized into binding activity, growth and proliferation, cell cycle progression, apoptosis, and invasion and metastasis. The microarray data were validated by positive correlation with quantitative real-time RT-PCR on 16 selected genes including igf1, tgfb3 and nov. Important carcinogenic pathways involved in RPM formation included insulin-like growth factor 1 (IGF-1), p38 MAPkinase, Wnt/beta-catenin and integrin signaling pathways. This study demonstrated that mesotheliomas in rats exposed to o-NT- and BCA were similar to mesotheliomas in humans, at least at the cellular and molecular level.

Acetates↗

ATM requirement in gene expression responses to ionizing radiation in human lymphoblasts and fibroblasts.

The heritable disorder ataxia telangiectasia (AT) is caused by mutations in the AT-mutated (ATM) gene with manifestations that include predisposition to lymphoproliferative cancers and hypersensitivity to ionizing radiation (IR). We investigated gene expression changes in response to IR in human lymphoblasts and fibroblasts from seven normal and seven AT-affected individuals. Both cell types displayed ATM-dependent gene expression changes after IR, with some responses shared and some responses varying with cell type and dose. Interestingly, after 5 Gy IR, lymphoblasts displayed ATM-independent responses not seen in the fibroblasts at this dose, which likely reflect signaling through ATM-related kinases, e.g., ATR, in the absence of ATM function.

Ataxia Telangiectasia↗

Profiles of global gene expression in ionizing-radiation-damaged human diploid fibroblasts reveal synchronization behind the G1 checkpoint in a G0-like state of quiescence.

Cell cycle arrest and stereotypic transcriptional responses to DNA damage induced by ionizing radiation (IR) were quantified in telomerase-expressing human diploid fibroblasts. Analysis of cytotoxicity demonstrated that 1.5 Gy IR inactivated colony formation by 40-45% in three fibroblast lines; this dose was used in all subsequent analyses. Fibroblasts exhibited > 90% arrest of progression from G2 to M at 2 hr post-IR and a similarly severe arrest of progression from G1 to S at 6 and 12 hr post-IR. Normal rates of DNA synthesis and mitosis 6 and 12 hr post-IR caused the S and M compartments to empty by > 70% at 24 hr. Global gene expression was analyzed in IR-treated cells. A microarray analysis algorithm, EPIG, identified nine IR-responsive patterns of gene expression that were common to the three fibroblast lines, including a dominant p53-dependent G1 checkpoint response. Many p53 target genes, such as CDKN1A, GADD45, BTG2, and PLK3, were significantly up-regulated at 2 hr post-IR. Many genes whose expression is regulated by E2F family transcription factors, including CDK2, CCNE1, CDC6, CDC2, MCM2, were significantly down-regulated at 24 hr post-IR. Numerous genes that participate in DNA metabolism were also markedly repressed in arrested fibroblasts apparently as a result of cell synchronization behind the G1 checkpoint. However, cluster and principal component analyses of gene expression revealed a profile 24 hr post-IR with similarity to that of G0 growth quiescence. The results reveal a highly stereotypic pattern of response to IR in human diploid fibroblasts that reflects primarily synchronization behind the G1 checkpoint but with prominent induction of additional markers of G0 quiescence such as GAS1.

Cell Line↗

Global gene expression associated with hepatocarcinogenesis in adult male mice induced by in utero arsenic exposure.

Our previous work has shown that exposure to inorganic arsenic in utero produces hepatocellular carcinoma (HCC) in adult male mice. To explore further the molecular mechanisms of transplacental arsenic hepatocarcinogenesis, we conducted a second arsenic transplacental carcinogenesis study and used a genomewide microarray to profile arsenic-induced aberrant gene expression more extensively. Briefly, pregnant C3H mice were given drinking water containing 85 ppm arsenic as sodium arsenite or unaltered water from days 8 to 18 of gestation. The incidence of HCC in adult male offspring was increased 4-fold and tumor multiplicity 3-fold after transplacental arsenic exposure. Samples of normal liver and liver tumors were taken at autopsy for genomic analysis. Arsenic exposure in utero resulted in significant alterations (p < 0.001) in the expression of 2,010 genes in arsenic-exposed liver samples and in the expression of 2,540 genes in arsenic-induced HCC. Ingenuity Pathway Analysis revealed that significant alterations in gene expression occurred in a number of biological networks, and Myc plays a critical role in one of the primary networks. Real-time reverse transcriptase-polymerase chain reaction and Western blot analysis of selected genes/proteins showed > 90% concordance. Arsenic-altered gene expression included activation of oncogenes and HCC biomarkers, and increased expression of cell proliferation-related genes, stress proteins, and insulin-like growth factors and genes involved in cell-cell communications. Liver feminization was evidenced by increased expression of estrogen-linked genes and altered expression of genes that encode gender-related metabolic enzymes. These novel findings are in agreement with the biology and histology of arsenic-induced HCC, thereby indicating that multiple genetic events are associated with transplacental arsenic hepatocarcinogenesis.

Animals↗

Standardizing global gene expression analysis between laboratories and across platforms.

To facilitate collaborative research efforts between multi-investigator teams using DNA microarrays, we identified sources of error and data variability between laboratories and across microarray platforms, and methods to accommodate this variability. RNA expression data were generated in seven laboratories, which compared two standard RNA samples using 12 microarray platforms. At least two standard microarray types (one spotted, one commercial) were used by all laboratories. Reproducibility for most platforms within any laboratory was typically good, but reproducibility between platforms and across laboratories was generally poor. Reproducibility between laboratories increased markedly when standardized protocols were implemented for RNA labeling, hybridization, microarray processing, data acquisition and data normalization. Reproducibility was highest when analysis was based on biological themes defined by enriched Gene Ontology (GO) categories. These findings indicate that microarray results can be comparable across multiple laboratories, especially when a common platform and set of procedures are used.

Gene Expression Profiling↗

Differential gene expression profiling in whole blood during acute systemic inflammation in lipopolysaccharide-treated rats.

Microarrays have been used to evaluate the expression of thousands of genes in various tissues. However, few studies have investigated the change in gene expression profiles in one of the most easily accessible tissues, whole blood. We utilized an acute inflammation model to investigate the possibility of using a cDNA microarray to measure the gene expression profile in the cells of whole blood. Blood was collected from male Sprague-Dawley rats at 2 and 6 h after treatment with 5 mg/kg (ip) LPS. Hematology showed marked neutrophilia accompanied by lymphopenia at both time points. TNF-alpha and IL-6 levels were markedly elevated at 2 h, indicating acute inflammation, but by 6 h the levels had declined. Total RNA was isolated from whole blood and hybridized to the National Institute of Environmental Health Sciences Rat Chip v.3.0. LPS treatment caused 226 and 180 genes to be differentially expressed at 2 and 6 h, respectively. Many of the differentially expressed genes are involved in inflammation and the acute phase response, but differential expression was also noted in genes involved in the cytoskeleton, cell adhesion, oxidative respiration, and transcription. Real-time RT-PCR confirmed the differential regulation of a representative subset of genes. Principal component analysis of gene expression discriminated between the acute inflammatory response apparent at 2 h and the observed recovery underway at 6 h. These studies indicate that, in whole blood, changes in gene expression profiles can be detected that are reflective of inflammation, despite the adaptive shifts in leukocyte populations that accompany such inflammatory processes.

Animals↗

Microarray data analysis of mouse neoplasia.

Microarray gene expression analysis offers great promise to help us understand the molecular events of experimental carcinogenesis, but have such promises been fulfilled? Studies of gene expression profiles of rodent are being published and demonstrate that yes, indeed, gene array data is furthering our understanding of tumor biology. Recent studies have identified differentially expressed genes in rodent mammary, colon, lung, and liver tumors. Although relatively few genes on the rodent arrays have been fully characterized, information has been generated to better identify signatures of histologic type and grade, understand invasion and metastasis, identify candidate biomarkers of early development, identify gene networks in carcinogenesis, understand responses to therapy, and decifer overlap with molecular events in human cancers. Data from mouse lung, mammary gland, and liver tumor studies are reviewed as examples of how to approach and interpret gene array data. Methods of gene array data analysis were also applied for discovery of genes involved in the regression of mouse liver tumors induced by chlordane, a nongenotoxic murine hepatocarcinogen. Promises are beginning to be fulfilled and it is clear that pathologists and toxicologists, in collaboration with molecular biologists, bioinformatists,and other scientists are making great strides in the design, analysis, and interpretation of microarray data for cancer studies.

Animals↗

Systematic variation normalization in microarray data to get gene expression comparison unbiased.

Normalization removes or minimizes the biases of systematic variation that exists in experimental data sets. This study presents a systematic variation normalization (SVN) procedure for removing systematic variation in two channel microarray gene expression data. Based on an analysis of how systematic variation contributes to variability in microarray data sets, our normalization procedure includes background subtraction determined from the distribution of pixel intensity values from each data acquisition channel and log conversion, linear or non-linear regression, restoration or transformation, and multiarray normalization. In the case when a non-linear regression is required, an empirical polynomial approximation approach is used. Either the high terminated points or their averaged values in the distributions of the pixel intensity values observed in control channels may be used for rescaling multiarray datasets. These pre-processing steps remove systematic variation in the data attributable to variability in microarray slides, assay-batches, the array process, or experimenters. Biologically meaningful comparisons of gene expression patterns between control and test channels or among multiple arrays are therefore unbiased using normalized but not unnormalized datasets.

Algorithms↗

Analysis of ATP-binding cassette transporter expression in drug-selected cell lines by a microarray dedicated to multidrug resistance.

Discovery of the multidrug resistance protein 1 (MDR1), an ATP-binding cassette (ABC) transporter able to transport many anticancer drugs, was a clinically relevant breakthrough in multidrug resistance research. Although the overexpression of ABC transporters such as P-glycoprotein/ABCB1, MRP1/ABCC1, and MXR/ABCG2 seems to be a major cause of failure in the treatment of cancer, acquired resistance to multiple anticancer drugs may also be multifactorial, involving alteration of detoxification processes, apoptosis, DNA repair, drug uptake, and overexpression of other ABC transporters. As a tool for the study of such phenomena, we designed and created a microarray platform, the ABC-ToxChip, to evaluate relative levels of transcriptional activation among genes involved in the various mechanisms of resistance. In the ABC-ToxChip, a comprehensive set of genes important in toxicological responses (represented by 2200 cDNA probes) is complemented with probes specifically matching ABC transporters as well as oligonucleotides representing 18,000 unique human genes. By comparing the transcriptional profiles of KB-3-1 and DU-145 parental cells with resistant derivatives selected in colchicine (KB-8-5), and 9-nitro-camptothecin (RCO.1), respectively, we demonstrate that ABC transporters (ABCB1/MDR1 and ABCC2/MRP2, respectively) show dramatic overexpression, whereas the glutathione S-transferase gene GST-Pi shows the strongest decrease in expression among the 20,000 genes studied. The results were confirmed by quantitative reverse transcription-polymerase chain reaction and immunohistochemistry. The custom-designed ABC-Tox microarray presented here will be helpful to elucidate mechanisms leading to anticancer drug resistance.

ATP-Binding Cassette Transporters↗

Gene expression profiling of rat livers reveals indicators of potential adverse effects.

This study tested the hypothesis that gene expression profiling can reveal indicators of subtle injury to the liver induced by a low dose of a substance that does not cause overt toxicity as defined by conventional criteria of toxicology (e.g., abnormal clinical chemistry and histopathology). For the purpose of this study we defined this low dose as subtoxic, i.e., a dose that elicits effects which are below the detection of conventional toxicological parameters. Acetaminophen (APAP) was selected as a model hepatotoxicant because (1) considerable information exists concerning the mechanism of APAP hepatotoxicity that can occur following high doses, (2) intoxication with APAP is the leading cause of emergency room visits involving acute liver failure within the United States, and (3) conventional clinical markers have poor predictive value. Rats treated with a single dose of 0, 50, 150, or 1500 mg/kg APAP were examined at 6, 24, or 48 h after exposure for conventional toxicological parameters and for gene expression alterations. Patterns of gene expression were found which indicated cellular energy loss as a consequence of APAP toxicity. Elements of these patterns were apparent even after exposure to subtoxic doses. With increasing dose, the magnitude of changes increased and additional members of the same biological pathways were differentially expressed. The energy loss suggested by gene expression changes was confirmed at the 1500 mg/kg dose exposure by measuring ATP levels. Only by ultrastructural examination could any indication of toxicity be identified after exposure to a subtoxic dose of APAP and that was occasional mitochondrial damage. In conclusion, this study provides evidence that supports the hypothesis that gene expression profiling may be a sensitive means of identifying indicators of potential adverse effects in the absence of the occurrence of overt toxicity.

Acetaminophen↗

Cross-site comparison of gene expression data reveals high similarity.

Consistency and coherence of gene expression data across multiple sites depends on several factors such as platform (oligo, cDNA, etc.), environmental conditions at each laboratory, and data quality. The Hepatotoxicity Working Group of the International Life Sciences Institute Health and Environmental Sciences Institute consortium on the application of genomics to mechanism-based risk assessment is investigating these factors by comparing high-density gene expression data sets generated on two sets of RNA from methapyrilene (MP) experiments conducted at Abbott Laboratories and Boehringer-Ingelheim Pharmaceuticals, Inc. using a single platform (Affymetrix Rat Genome U34A GeneChip) at seven different sites. This article focuses on the evaluation of data quality and statistical models that facilitate the comparison of such data sets at the probe level. We present methods for exploring and quantitatively assessing differences in the data, with the principal goal being the generation of lists of site-insensitive genes responsive to low and high doses of MP. A combination of numerical and graphical techniques reveals important patterns and partitions of variability in the data, including the magnitude of the site effects. Although the site effects are significantly large in the analysis results, they appear to be primarily additive and therefore can be adjusted in the statistical calculations in a way that does not bias conclusions regarding treatment differences.

Animals↗

Identification of distinct and common gene expression changes after oxidative stress and gamma and ultraviolet radiation.

The human genome is exposed to many different kinds of DNA-damaging agents. While most damage is detected and repaired through complex damage recognition and repair machineries, some damage has the potential to escape these mechanisms. Unrepaired DNA damage can give rise to alterations and mutations in the genome in an individual cell, which can result in malignant transformation, especially when critical genes are deregulated. In this study, we investigated gene expression changes in response to oxidative stress, gamma (gamma) radiation, and ultraviolet (UV) radiation and their potential implications in cancer development. Doses were selected for each of the three treatments, based on their ability to cause a similar G(1) checkpoint induction and slow down in early S-phase progression, as reflected by a comparable reduction in cyclin E-associated kinase activity of at least 75% in logarithmically growing human dermal diploid fibroblasts. To investigate gene expression changes, logarithmically growing dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 microM of tert-butyl hydroperoxide (t-butyl-OOH)), or UV radiation (UVC) (7.5 J/m(2)) and RNA was harvested 6 h after treatment. Gene expression was analyzed using the NIEHS Human ToxChip 2.0 with approximately 1901 cDNA clones representing known genes and expressed sequence tags (ESTs). We were able to identify common and distinct responses in dermal diploid fibroblasts to the three different stimuli used. Within our analysis, gene expression profiles in response to gamma radiation and oxidative stress appeared to be more similar than profiles expressed after UV radiation. Interestingly, equivalent cyclin E-associated kinase activity reduction with all the three treatments was associated with greater transcriptional changes after UV radiation than after gamma radiation and oxidative stress. While samples treated with UV radiation displayed modulations of their mitogen activated protein kinase (MAPK) pathway, gamma radiation had its major influence on cell-cycle progression in S-phase and mitosis. In addition, cell cultures from different individuals displayed significant differences in their gene expression responses to DNA damage.

Adult↗

Changes in global gene and protein expression during early mouse liver carcinogenesis induced by non-genotoxic model carcinogens oxazepam and Wyeth-14,643.

We hypothesized that the mouse liver tumor response to non-genotoxic carcinogens would involve some common early gene and protein expression changes that could ultimately be used to predict chemical hepatocarcinogenesis. In order to identify a panel of genes to test, we analyzed global differences in gene and protein expression in livers from B6C3F1 mice following dietary treatment with two rodent carcinogens, the benzodiazepine anti-anxiety drug oxazepam (2500 p.p.m.) and the hypolipidemic agent Wyeth (Wy)-14,643 (500 p.p.m.) compared with livers from untreated mice. Male mice were exposed for 2 weeks and 1, 3 or 6 months to oxazepam or Wy-14,643 in an age-matched study design. By histopathological evaluation, no liver preneoplastic foci or tumors were detected at 6 months in treated or control groups. By cDNA microarray analysis [NIEHS Mouse Chip (8700 genes); n = 3 individual livers/group, four hybridizations/sample], expression of 36 genes or 220 genes were changed relative to control livers following 6 months of oxazepam or Wy-14,643 treatment, respectively. To obtain a more comprehensive picture of gene/protein expression changes, we also conducted a proteomics study by 2D-gel electrophoresis followed by matrix assisted laser desorption/ionization-mass spectrometry on cytoplasmic, nuclear, and microsomal subcellular fractions of the same liver samples utilized for the cDNA microarray analysis. Real-time PCR, western blot analysis and immunohistochemistry were utilized for validation and to expand the results to other time points. Cyp2b20, growth arrest- and damage-inducible gene beta (Gadd45beta), tumor necrosis factor alpha-induced protein 2 and insulin-like growth factor binding protein 1 (Igfbp5) genes and proteins were upregulated by oxazepam, and Cyp2b20, Cyclin D1, proliferating cell nuclear antigen, Igfbp5, Gadd45beta and cell death-inducing DNA fragmentation factor alpha subunit-like effector A exhibited higher expression after Wy-14,643 treatment. Most of these genes/proteins were also deregulated at 2 weeks. There appeared to be more distinct than common changes in the expression of carcinogenesis-related genes/proteins between the two compounds, suggesting that the major carcinogenic pathways are different for these compounds and may be distinct for different chemical classes.

Animals↗

ATM-dependent and -independent gene expression changes in response to oxidative stress, gamma irradiation, and UV irradiation.

Ataxia telangiectasia (AT) is an autosomal recessive disorder characterized by progressive cerebellar degeneration, immunodeficiencies, telangiectasias, sensitivity to ionizing radiation, and high predisposition for malignancies. The ataxia telangiectasia mutated (ATM) gene encodes a protein (ATM) with serine/threonine kinase activity. DNA-double strand breaks are known to increase its kinase activity. While cells from individuals with AT are attenuated in their G(1)-, S- and G(2)-phase cell cycle checkpoint functions in response to gamma irradiation and oxidative stress, their response to UV irradiation appears to be equivalent to that of wild-type cells. In this study, we investigated changes in gene expression in response to gamma irradiation, oxidative stress, and UV irradiation, focusing on the dependence on ATM. Doses for all three treatments were selected that resulted in roughly an equivalent induction of a G(1) checkpoint response and inhibition of progression through S phase. To investigate gene expression changes, logarithmically growing wild-type and AT dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 micro M t-butyl-hydroperoxide), or UV radiation (7.5 J/m(2)), and RNA was harvested 6 h after treatment. Gene expression analysis was performed using the NIEHS Human ToxChip 2.0 with approximately 1900 cDNA clones representing known genes and ESTs. All three treatments resulted in distinct patterns of gene expression changes, as shown previously. ATM-dependent and ATM-independent components were detected within these patterns, as were novel indications of involvement of ATM in regulation of transcription factors such as SP1, AP1 and MTF1.

Algorithms↗

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

Assessment of the impact of xenobiotic exposure on human health and disease progression is complex. Knowledge of mode(s) of action, including mechanism(s) contributing to toxicity and disease progression, is valuable for evaluating compounds. Toxicogenomics, the subdiscipline which merges genomics with toxicology, holds the promise to contributing significantly toward the goal of elucidating mechanism(s) by studying genome-wide effects of xenobiotics. Global gene expression profiling, revolutionized by microarray technology and a crucial aspect of a toxicogenomic study, allows measuring transcriptional modulation of thousands of genes following exposure to a xenobiotic. We use our results from previous studies on compounds representing two different classes of xenobiotics (barbiturate and peroxisome proliferator) to discuss the application of computational approaches for analyzing microarray data to elucidate mechanism(s) underlying cellular responses to toxicants. In particular, our laboratory demonstrated that chemical-specific patterns of gene expression can be revealed using cDNA microarrays. Transcript profiling provides discrimination between classes of toxicants, as well as, genome-wide insight into mechanism(s) of toxicity and disease progression. Ultimately, the expectation is that novel approaches for predicting xenobiotic toxicity in humans will emerge from such information.

Computational Biology↗

ATR enforces the topoisomerase II-dependent G2 checkpoint through inhibition of Plk1 kinase.

An ATR-dependent G(2) checkpoint responds to inhibition of topoisomerase II and delays entry into mitosis by sustaining nuclear exclusion of cyclin B1-Cdk1 complexes. Here we report that induction of this checkpoint with ICRF-193, a topoisomerase II catalytic inhibitor that does not cause DNA damage, was associated with an ATR-dependent inhibition of polo-like kinase 1 (Plk1) kinase activity and a decrease in cyclin B1 phosphorylation. Expression of constitutively active Plk1 but not wild type Plk1 reversed ICRF-193-induced mitotic delay in HeLa cells, suggesting that Plk1 kinase activity is important for the checkpoint response to ICRF-193. G(2)/M synchronized normal human fibroblasts, when treated with ICRF-193, showed a decrease in cyclin B1 phosphorylation and Plk1 kinase activity despite high cyclin B1-Cdk1 kinase activity. G(2) fibroblasts that were treated with caffeine to override the checkpoint response to ICRF-193 displayed a high incidence of chromosomal aberrations. Taken together, these results suggest that ATR-dependent inhibition of Plk1 kinase activity may be one mechanism to regulate cyclin B1 phosphorylation and sustain nuclear exclusion during the G(2) checkpoint response to topoisomerase II inhibition. Moreover, the results demonstrate an important role for the topoisomerase II-dependent G(2) checkpoint in the preservation of human genomic stability.

Alleles↗

Computational selection of distinct class- and subclass-specific gene expression signatures.

In this investigation we used statistical methods to select genes with expression profiles that partition classes and subclasses of biological samples. Gene expression data corresponding to liver samples from rats treated for 24 h with an enzyme inducer (phenobarbital) or a peroxisome proliferator (clofibrate, gemfibrozil or Wyeth 14,643) were subjected to a modified Z-score test to identify gene outliers and a binomial distribution to reduce the probability of detecting genes as differentially expressed by chance. Hierarchical clustering of 238 statistically valid differentially expressed genes partitioned class-specific gene expression signatures into groups that clustered samples exposed to the enzyme inducer or to peroxisome proliferators. Using analysis of variance (ANOVA) and linear discriminant analysis methods we identified single genes as well as coupled gene expression profiles that separated the phenobarbital from the peroxisome proliferator treated samples and discerned the fibrate (gemfibrozil and clofibrate) subclass of peroxisome proliferators. A comparison of genes ranked by ANOVA with genes assessed as significant by mixedlinear models analysis [J. Comput. Biol. 8 (2001) 625] or ranked by information gain revealed good congruence with the top 10 genes from each statistical method in the contrast between phenobarbital and peroxisome proliferators expression profiles. We propose building upon a classification regimen comprised of analysis of replicate data, outlier diagnostics and gene selection procedures to utilize cDNA microarray data to categorize subclasses of samples exposed to pharmacologic agents.

Algorithms↗