Gene expression profile in the livers of rats orally administered ethinylestradiol for 28 days using a microarray technique

Dietary exposure of 17-alpha ethinylestradiol modulates physiological endpoints and gene signaling pathways in female largemouth bass (Micropterus salmoides)

17Alpha-ethinylestradiol (EE2), used for birth control in humans, is a potent estrogen that is found in wastewater at low concentrations (ng/l). EE2 has the ability to interfere with the endocrine system of fish, affecting reproduction which can result in population level effects. The objective of this study was to determine if dietary exposure to EE2 would alter gene expression patterns and key pathways in the liver and ovary and whether these could be associated with reproductive endpoints in female largemouth bass during egg development. Female LMB received 70ng EE2/g feed (administered at 1% of body weight) for 60 days. EE2 dietary exposure significantly reduced plasma vitellogenin concentrations by 70%. Hepatosomatic and gonadosomatic indices were also decreased with EE2 feeding by 38.5% and 40%, respectively. Transcriptomic profiling revealed that there were more changes in steady state mRNA levels in the liver compared to the ovary. Genes associated with reproduction were differentially expressed, such as vitellogenin in the liver and aromatase in the gonad. In addition, a set of genes related with oxidative stress (e.g. glutathione reductase and glutathione peroxidase) were identified as altered in the liver and genes associated with the immune system (e.g. complement component 1, and macrophage-inducible C-type lectin) were altered in the gonad. In a follow-up study with 0.2ng EE2/g feed for 60 days, similar phenotypic and gene expression changes were observed that support these findings with the higher concentrations. This study provides new insights into how dietary exposure to EE2 interferes with endocrine signaling pathways in female LMB during a critical period of reproductive oogenesis.

Toxicogenomic biomarkers for liver toxicity

Toxicogenomics (TGx) is a widely used technique in the preclinical stage of drug development to investigate the molecular mechanisms of toxicity. A number of candidate TGx biomarkers have now been identified and are utilized for both assessing and predicting toxicities. Further accumulation of novel TGx biomarkers will lead to more efficient, appropriate and cost effective drug risk assessment, reinforcing the paradigm of the conventional toxicology system with a more profound understanding of the molecular mechanisms of drug-induced toxicity. In this paper, we overview some practical strategies as well as obstacles for identifying and utilizing TGx biomarkers based on microarray analysis. Since clinical hepatotoxicity is one of the major causes of drug development attrition, the liver has been the best documented target organ for TGx studies to date, and we therefore focused on information from liver TGx studies. In this review, we summarize the current resources in the literature in regard to TGx studies of the liver, from which toxicologists could extract potential TGx biomarker gene sets for better hepatotoxicity risk assessment.

Gene expression analysis of estrogenic compounds in the liver of common carp (Cyprinus carpio) using a custom cDNA microarray

Exposure to a variety of compounds with estrogenic activity has been shown to interfere with normal developmental and reproductive processes in various vertebrate species. The aim of this study was to determine the transcriptional profile of the natural estrogen, 17 beta-estradiol, and three synthetic estrogenic compounds (4-nonylphenol, bisphenol A, ethinylestradiol) in the liver of common carp, using a custom cDNA microarray. For that purpose, fish were aqueously exposed to three concentrations of each chemical for 24 or 96 h. Microarray analysis revealed that a total of 185 different gene transcripts were differentially expressed following exposure to at least one of the estrogen(-like) concentrations. We were able to identify a common set of 28 gene fragments, whose expression was significantly modified in the same way by the three xenoestrogens and 17 beta-estradiol. Although several of these gene expression effects corroborated past literature data, we also discovered some novel target genes of (xeno)estrogen exposure, providing interesting insights into the molecular basis of estrogenic effects. In addition, each of the four compounds induced gene expression changes that were not, or only partially, shared by the other chemicals, suggesting that not all chemicals with estrogenic activity act alike. These results demonstrate the potential of our custom Cyprinus carpio microarray to detect common estrogen-like activity as well as to identify unique compound-associated effects of (estrogenic) endocrine disruptors in fish.

Biomarkers for exposure to estrogenic compounds: gene expression analysis in zebrafish (Danio rerio)

Gene expression analyses in male zebrafish (Danio rerio) were carried out using microarray technique and quantitative polymerase chain reaction. Genes responding to the exposure to 17beta-estradiol, bisphenol A and genistein were identified, among them genes involved in metabolism, reproductional and developmental processes. Threshold levels of 17beta-estradiol (200 ng/L), bisphenol A (2000 microg/L), and genistein (5000 microg/L) for the upregulation of the vtg1 gene in short-time exposures (11 days) were determined by qPCR. 14k microarrays were used to generate complete lists of genes regulated by these estrogenic compounds. For this purpose, liver samples from 10 exposed zebrafish and 10 controls were processed. In this case the expressions of 211 genes were significantly regulated by 17beta-estradiol, 47 by bisphenol A and 231 by genistein. Furthermore, it is shown that fish exposed to 17beta-estradiol and genistein have similarities in their gene expression patterns, whereas bisphenol A apparently affected gene expression in a different way. Only genes coding for egg-yolk precursor protein vitellogenin were found to be regulated by all three compounds, which shows that these genes are the only suitable markers for exposure to different estrogenic compounds. The regulated genes were assigned to gene ontology classes. All three estrogenic compounds regulated genes mainly involved in primary and cellular metabolism, but genistein regulated several genes involved in cell cycle-regulation and bisphenol A several genes involved in protein biosynthesis. Genistein also upregulated the expression of four eggshell proteins, which can be used as biomarkers for exposure to this chemical.

o,p'-DDT elicits PXR/CAR-, not ER-, mediated responses in the immature ovariectomized rat liver

Technical-grade dichlorodiphenyltrichloroethane (DDT) is an agricultural pesticide and malarial vector control agent that has been designated a potential human hepatocarcinogen. The o,p'-enantiomer exhibits estrogenic activity that has been associated with the carcinogenicity of DDT. The temporal and dose-dependent hepatic estrogenicity of o,p'-DDT was investigated using complementary DNA microarrays in immature ovariectomized Sprague-Dawley rats with complementary histopathology and tissue-level analysis. Animals were gavaged with 300 mg/kg o,p'-DDT either once or once daily for 3 consecutive days. Liver samples were examined 2, 4, 8, 12, 18, or 24 h after a single dose or following three daily doses. For dose-response studies, a single dose of 3, 10, 30, 100, or 300 mg/kg body weight o,p'-DTT was administered for 3 consecutive days. Genes associated with drug metabolism (Cyp2b2 and Cyp3a2), the nuclear receptors constitutive androstane receptor (CAR) and pregnane X receptor (PXR), cell proliferation (Ccnd1, Ccnb1, Ccnb2, and Stmn1), and oxidative stress (Gclm and Hmox1) were significantly induced. Cyp2b2 exhibited dose-dependent regulation and was significantly induced across all time points, while cell proliferation- and oxidative stress-related genes exhibited transient induction. The induction of Cyp2b2 and Cyp3a2 mRNA levels suggest PXR/CAR activation, consistent with expression of genes associated with oxidative stress. Few genes known to be estrogen receptor (ER) regulated were differentially expressed when compared to the hepatic gene expression profile elicited by ethynyl estradiol in immature ovariectomized C57BL/6 mice using the same study design and analysis methods. These data indicate that o,p'-DDT elicits PXR/CAR-, not ER-, mediated gene expression in the rat liver. Based on the species-specific differences in CAR regulation, the extrapolation of rodent DDT hepatocarcinogenicity to humans warrants further investigation.

Identification of a transcriptional fingerprint of estrogen exposure in rainbow trout liver

The goal of this study was to identify a set of hepatic genes regulated by ligand-dependent activation of the estrogen receptor in juvenile rainbow trout (Oncorhynchus mykiss). A custom rainbow trout oligo DNA microarray, which contains probes targeting approximately 1450 genes relevant to carcinogenesis, toxicology, endocrinology, and stress physiology was utilized to identify transcriptional fingerprints of in vivo dietary exposure to 17 beta-estradiol (E2), tamoxifen (TAM), estradiol + tamoxifen (E2 + TAM), diethylstilbestrol (DES), dehydroepiandrosterone (DHEA), dihydrotestosterone (DHT), and cortisol (CORT). Estrogen exposure altered the expression of up to 49 genes involved in reproduction, immune response, cell growth, transcriptional regulation, protein synthesis and modification, drug metabolism, redox regulation, and signal transduction. E2, DES, and DHEA regulated 18 genes in common, mostly those associated with vitellogenesis, cell proliferation, and signal transduction. Interestingly, DHEA uniquely regulated several complement component genes of importance to immune response. While the effect of TAM on E2-induced changes in gene expression was mostly antagonistic, TAM alone increased expression of VTG1 and other genes associated with egg development and immune response. Few genes responded to CORT treatment, and DHT significantly altered expression of only one gene targeted by the OSUrbt array. Hierarchical cluster and principal components analyses revealed distinct patterns of gene expression corresponding to estrogens and non-estrogens, though unique patterns could also be detected for individual chemicals. A set of estrogen-responsive genes has been identified that can serve as a biomarker of environmental exposure to xenoestrogens.

The Evolving Role of Drug Metabolism in Drug Discovery and Development

Drug metabolism in pharmaceutical research has traditionally focused on the well-defined aspects of absorption, distribution, metabolism and excretion, commonly-referred to ADME properties of a compound, particularly in the areas of metabolite identification, identification of drug metabolizing enzymes (DMEs) and associated metabolic pathways, and reaction mechanisms. This traditional emphasis was in part due to the limited scope of understanding and the unavailability of in vitro and in vivo tools with which to evaluate more complex properties and processes. However, advances over the past decade in separate but related fields such as pharmacogenetics, pharmacogenomics and drug transporters, have dramatically shifted the drug metabolism paradigm. For example, knowledge of the genetics and genomics of DMEs allows us to better understand and predict enzyme regulation and its effects on exogenous (pharmacokinetics) and endogenous pathways as well as biochemical processes (pharmacology). Advances in the transporter area have provided unprecedented insights into the role of transporter proteins in absorption, distribution, metabolism and excretion of drugs and their consequences with respect to clinical drug-drug and drug-endogenous substance interactions, toxicity and interindividual variability in pharmacokinetics. It is therefore essential that individuals involved in modern pharmaceutical research embrace a fully integrated approach and understanding of drug metabolism as is currently practiced. The intent of this review is to reexamine drug metabolism with respect to the traditional as well as current practices, with particular emphasis on the critical aspects of integrating chemistry and biology in the interpretation and application of metabolism data in pharmaceutical research.

Cross-species comparative toxicogenomics as an aid to safety assessment

Cross-species comparative toxicogenomics has the potential for improving the understanding of the different responses of animal models to toxicants at a molecular level. This understanding could then lead to a more accurate extrapolation of the risk posed by these toxicants to humans. Cross-species comparative studies have been carried out at the genomic sequence level and using microarrays to examine changes in global mRNA profiles. However, these studies face considerable bioinformatic challenges in terms of identifying which genes are truly orthologous across species. The resources to analyse such studies, in the context of such orthologues, beg improvement. Finally, the experimental design of such studies needs to be carefully considered to make their results fully interpretable. These issues are discussed, along with the current state-of-the-art cross-species comparative toxicogenomics in this review.

Gene expression responses of European flounder (Platichthys flesus) to 17-beta estradiol

Male European flounder (Platichthys flesus) were intraperitoneally injected with 10mg/kg 17-beta estradiol and tissues taken from individuals over a timecourse of 16 days. The GENIPOL P. flesus cDNA microarray was employed to detect hepatic gene expression differences between fish treated with estradiol and saline controls. Known biomarkers of estrogen exposure, choriogenin L and vitellogenins, showed sustained induction over the time-course. Among 175 identified clones showing sustained statistically significant induction or repression, those associated with the Gene Ontology terms mitochondria, amino acid synthesis, ubiquitination and apoptosis were included amongst those induced while those associated with immune function, electron transport, cell signalling and protein phosphorylation were repressed. Thus, we show the gene expression response of an environmentally relevant fish species to a high dose of an estrogenic endocrine disruptor and also report the sequencing of a further 2121 flounder ESTs.

Novel Microarray-Based Method for Estimating Exposure to Ionizing Radiation

Accurate estimation of the dose of ionizing radiation to which individuals have been exposed is critical for therapeutic treatment. We investigated whether gene expression profiles could be used to evaluate the dose received, thereby serving as a biological dosimeter. We used cDNA microarrays to monitor changes in gene expression profiles induced by ionizing radiation in mouse total blood. The subsets of genes best characterizing each dose were identified by resampling the original data set and calculating the intersection of the dose signatures. This analytical strategy minimizes the impact of potential genetic/epigenetic variation between mice and overcomes the bias in gene selection inherent to microarray technology. The significance of the identified signatures was evaluated by monitoring the type I error rate by in silico negative control simulation. Based on the distribution of the mean ratios of the selected probes, we were able to identify transcription profiles giving 83% to 100% correct estimation of the dose received by test mice, demonstrating that the selected probes could be used to determine the dose of radiation to which the animals had been exposed. This method could potentially be generalized to determine the level of exposure to other toxins and could be used to develop new related clinical applications.

Hepatic gene expression profiling using Genechips in zebrafish exposed to 17alpha-ethynylestradiol

Genomic, proteomic, and metabolomic technologies continue to receive increasing interest from environmental toxicologists. This interest is due to the great potential of these technologies to identify detailed modes of action and to provide assistance in the evaluation of a contaminant's risk to aquatic organisms. Our experimental model is the zebrafish (Danio rerio) exposed to reference endocrine disrupting compounds in order to investigate compound-induced changes in gene transcript profiles. Adult, female zebrafish were exposed to 0, 15, 40, and 100ng/L of 17alpha-ethynylestradiol (EE2) and concentration and time-dependent changes in hepatic gene expression were examined using Affymetrix GeneChip Zebrafish Genome Microarrays. At 24, 48, and 168h, fish were sacrificed and liver mRNA was extracted for gene expression analysis (24 and 168h only). In an effort to link gene expression changes to effects on higher levels of biological organization, body and ovary weights were measured and blood was collected for measurement of plasma steroid hormones (17beta-estradiol (E2), testosterone (T)) and vitellogenin (VTG) using ELISA. EE2 exposure significantly affected gene expression, GSI, E2, T, and VTG. We observed 1622 genes that were significantly affected (p< or =0.001) in a concentration-dependent manner by EE2 exposure at either 24 or 168h. Gene ontology (GO) analysis revealed that EE2 exposure affected genes involved in hormone metabolism, vitamin A metabolism, steroid binding, sterol metabolism, and cell growth. Plasma VTG was significantly increased at 24, 48, and 168h (p< or =0.05) at 40 and 100ng/L and at 15ng/L at 168h. E2 and T were significantly reduced following EE2 exposure at 48 and 168h. GSI was decreased in a concentration-dependent manner at 168h. In this study, we identified genes involved in a variety of biological processes that have the potential to be used as markers of exposure to estrogenic substances. Future work will evaluate the use of these genes in zebrafish exposed to weak estrogens to determine if these genes are indicative of exposure to estrogens with varying potencies.

Risk Assessment of Endocrine Disrupters: The Role of Toxicological Studies

Endocrine disrupters (ED) represent a good challenge for experimental toxicology. In order to deal with several critical points relevant to risk assessment: (a) ED may induce long-term effects upon exposure in susceptible developmental phases, including postnatal life up to puberty; thus, efforts are required to refine testing strategies, for example, by supporting the two-generation rodent study with a comprehensive in vitro/in vivo screening battery; (b) due to the regulatory role of endocrine homeostasis, mechanisms of endocrine disruption may impact on immune, neurobehavioral, and reproductive development, as well as on susceptibility to cancer; (c) the potential multiple exposure to ED with common targets through diet and/or living environment calls for the development of models to understand mechanisms of interactions and effects of mixtures; and (d) last but not least, ED may interact with a number of factors related to differential vulnerability of individuals or population subgroups, including the intake of nutrients or bioactive food components. Besides reducing the chance for noxious chemicals to enter our life, toxicological research on mechanisms may also lead to the definition of possible biomarkers of exposure, effect, and susceptibility that may be further exploited in human health surveillance.

Analysis of dose-response effects on gene expression data with comparison of two microarray platforms

MOTIVATION

The problems of analyzing dose effects on gene expression are gaining attention in biomedical research. A specific challenge is to detect genes with expression levels that change according to dose levels in a non-random manner, but nonetheless may be considered as potential biomarkers.

METHOD

We are among the first to formally apply a tool that uses an isotonic (monotonic) regression approach to this area of study. We introduce a test statistic to select genes with significant dose-response expression in a monotonic fashion based on a permutation procedure. We then compare the results with those achieved from the application of a likelihood ratio-based test.

RESULTS

We apply the isotonic regression approach to a study of gene expression in the RKO colon carcinoma cell line in response to varying dosage levels of the chemotherapeutic agent 5-fluorouracil. A feature of both Affymetrix and printed 75mer oligomer cDNA arrays produced from the same samples provides an opportunity to compare the two microarray platforms.

AVAILABILITY

Statistical software S-plus Code to implement the method is available from the authors.

CONTACT

kcoombes@mdanderson.org

Systems-ADME/Tox: resources and network approaches

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

Molecular pathological analysis on the mechanism of liver carcinogenesis in dicyclanil-treated mice

To investigate the mechanism of hepatocarcinogenesis due to dicyclanil (DC), an insect growth regulator for sheep, histopathological and molecular biological analyses were performed in the liver of male ICR mice fed on a diet containing 1500 ppm of DC for 2 weeks (Experiment I; Exp. I). In gene expression analyses using a large-scale cDNA microarray and RT-PCR, fluctuations of expressions of metabolism-/oxidation-/reduction-related genes, such as CYP1A, aldehyde dehydrogenase family 1 subfamily A1 (Aldh1a1), and thioredoxin reductase 1 (Txnrd1), were predominantly observed in the liver of the DC-treated group. In Experiment II (Exp. II), small-scale and metabolism/oxidative stress-specific cDNA microarray, real-time RT-PCR, and measurement of NF-kappaB protein were performed in the mice liver using a two-stage hepatocarcinogenesis model, in which the male ICR mice were fed on a diet containing 1500 ppm of DC for 7 weeks after a single injection of dimethylnitrosamine (DMN). These mice were subjected to two-thirds partial hepatectomy (PH) at week 3. During histopathological examinations, a remarkable increase in gamma-glutamyltransferase-positive cells was observed in the DMN+DC+PH group. During the microarray and PCR analyses, the metabolism and oxidative stress-related genes, such as Cyp1a, P450 oxidoreductase (Por), and thioredoxin reductase 1 (Txnrd1); a few DNA damage/repair genes, such as 8-oxoguanine DNA-glycosylase 1 (Ogg1); and growth arrest and DNA-damage-inducible 45 alpha (Gadd45a), were fluctuated in this group, together with a slight increase in the concentration of activated NF-kappaB. These results suggest that DNA damages due to oxidative stress may be involved in the mechanism of DC-induced hepatocarcinogenesis in mice.

Dietary exposure to soy or whey proteins alters colonic global gene expression profiles during rat colon tumorigenesis

BACKGROUND

We previously reported that lifetime consumption of soy proteins or whey proteins reduced the incidence of azoxymethane (AOM)-induced colon tumors in rats. To obtain insights into these effects, global gene expression profiles of colons from rats with lifetime ingestion of casein (CAS, control diet), soy protein isolate (SPI), and whey protein hydrolysate (WPH) diets were determined.

RESULTS

Male Sprague Dawley rats, fed one of the three purified diets, were studied at 40 weeks after AOM injection and when tumors had developed in some animals of each group. Total RNA, purified from non-tumor tissue within the proximal half of each colon, was used to prepare biotinylated probes, which were hybridized to Affymetrix RG_U34A rat microarrays containing probes sets for 8799 rat genes. Microarray data were analyzed using DMT (Affymetrix), SAM (Stanford) and pair-wise comparisons. Differentially expressed genes (SPI and/or WPH vs. CAS) were found. We identified 31 induced and 49 repressed genes in the proximal colons of the SPI-fed group and 44 induced and 119 repressed genes in the proximal colons of the WPH-fed group, relative to CAS. Hierarchical clustering identified the co-induction or co-repression of multiple genes by SPI and WPH. The differential expression of I-FABP (2.92-, 3.97-fold down-regulated in SPI and WPH fed rats; P = 0.023, P = 0.01, respectively), cyclin D1 (1.61-, 2.42-fold down-regulated in SPI and WPH fed rats; P = 0.033, P = 0.001, respectively), and the c-neu proto-oncogene (2.46-, 4.10-fold down-regulated in SPI and WPH fed rats; P < 0.001, P < 0.001, respectively) mRNAs were confirmed by real-time quantitative RT-PCR. SPI and WPH affected colonic neuro-endocrine gene expression: peptide YY (PYY) and glucagon mRNAs were down-regulated in WPH fed rats, whereas somatostatin mRNA and corresponding circulating protein levels, were enhanced by SPI and WPH.

CONCLUSIONS

The identification of transcripts co- or differentially-regulated by SPI and WPH diets suggests common as well as unique anti-tumorigenesis mechanisms of action which may involve growth factor, neuroendocrine and immune system genes. SPI and WPH induction of somatostatin, a known anti-proliferative agent for colon cancer cells, would inhibit tumorigenesis.