Comprehensive analysis of DNA methylation and gene expression of rat liver in a 2-stage hepatocarcinogenesis model.

Recent studies have shown that epigenetic alterations correlate with carcinogenesis in various tissues. Identification of these alterations might help characterize the early stages of carcinogenesis. We comprehensively analyzed DNA methylation and gene expression in livers obtained from rats exposed to nitrosodiethylamine (DEN) followed by a promoter of hepatic carcinogenesis, phenobarbital (PB). The combination of DEN and PB induced marked increases in number and area of glutathione S-transferase-placental form (GST-P)-positive foci in the liver. In the liver of rats that received 30 mg/kg of DEN, pathway analysis revealed alterations of common genes in terms of gene expression and DNA methylation, and that these alterations were related to immune responses. Hierarchical clustering analysis of the expression of common genes from public data obtained through the Toxicogenomics Project-Genomics Assisted Toxicity Evaluation system (TG-GATEs) showed that carcinogenic compounds clustered together. MBD-seq and GeneChip analysis indicated that major histocompatibility complex class Ib gene RT1-CE5, which has an important role in antigen presentation, was hypomethylated around the promoter region and specifically induced in the livers of DEN-treated rats. Further, immunohistochemical analysis indicated that the co-localization of GST-P and protein homologous to RT1-CE5 was present at the foci of some regions. These results suggest that common genes were altered in terms of both DNA methylation and expression in livers, with preneoplastic foci indicating carcinogenic potential, and that immune responses are involved in early carcinogenesis. In conclusion, the present study identified a specific profile of DNA methylation and gene expression in livers with preneoplastic foci. Early epigenetic perturbations of immune responses might correlate with the early stages of hepatocarcinogenesis.

Detection of initiating potential of non-genotoxic carcinogens in a two-stage hepatocarcinogenesis study in rats.

We previously reported a toxicogenomics-based prediction model for hepatocarcinogens in which the expression patterns of signature genes following repeated doses of either genotoxic or non genotoxic compounds were similar. Based on the results of our prediction model, we hypothesized that repeated doses of non-genotoxic carcinogens might have initiating potential. Here, we conducted a two stage hepatocarcinogenesis study in rats exposed to the initiating agent nitrosodiethylamine (DEN), and hepatotoxic compounds thioacetamide (TAA), methapyrilene (MP) and acetaminophen (APAP) for 1-2 weeks followed by the liver tumor promoter phenobarbital (PB). The duration of initial treatment was determined based on positive results from our prediction model. Combined treatment of 3 or 30 mg/kg of genotoxic DEN and PB induced marked increases in altered hepatocellular foci and a DEN dose-dependent increase in the number and area of glutathione S-transferase-placental form (GST-P)-positive foci. A low number of altered hepatocellular foci were also observed in rats treated with TAA at a dose of 45 mg/kg.MP at a dose of 100 mg/kg induced a very low number of foci, but APAP did not. Hierarchical clustering analysis using gene expression data revealed that 2-week treatment with TAA at a dose of 30 mg/kg and MP at 45 mg/kg induced specific expression of DNA damage-related genes, similar to 1-week treatment with DEN at a dose of 30 mg/kg. These results suggest that TAA and MP induce DNA damage, which partially supports our hypothesis. Although this study does not indicate whether tumor growth in response to these compounds can be assessed in this model, our results suggest that cumulative treatment with non genotoxic TAA might have initiating potential in the liver.

miRNA expression atlas in male rat.

MicroRNAs (miRNAs) are small (~22 nucleotide) noncoding RNAs that play pivotal roles in regulation of gene expression. The value of miRNAs as circulating biomarkers is now broadly recognized; such tissue-specific biomarkers can be used to monitor tissue injury and several pathophysiological conditions in organs. In addition, miRNA profiles of normal organs and tissues are important for obtaining a better understanding of the source of modulated miRNAs in blood and how those modulations reflect various physiological and toxicological conditions. This work was aimed at creating an miRNA atlas in rats, as part of a collaborative effort with the Toxicogenomics Informatics Project in Japan (TGP2). We analyzed genome-wide miRNA profiles of 55 different organs and tissues obtained from normal male rats using miRNA arrays. The work presented herein represents a comprehensive dataset derived from normal samples profiled in a single study. Here we present the whole dataset with miRNA profiles of multiple organs, as well as precise information on experimental procedures and organ-specific miRNAs identified in this dataset.

Gene expression profiling in rat liver treated with various hepatotoxic-compounds inducing coagulopathy.

A large-scale transcriptome database of rat liver (TG-GATEs) has been established by the Toxicogenomics Project in Japan. In the present study, we focused on 8 hepatotoxic compounds within TG-GATEs, i.e., clofibrate, omeprazole, ethionine, thioacetamide, benzbromarone, propylthiouracil, Wy-14,643 and amiodarone, which induced coagulation abnormalities. Aspirin was selected as a reference compound that directly causes coagulation abnormality, but not through liver toxicity. In blood chemical examinations, for all the coagulopathic compounds there was little elevation of aspartate aminotransferase (AST) and/or alanine aminotransferase (ALT), suggesting no severe cell death by treatment with the compounds. We extracted 344 probe sets from the data for these 8 typical drugs, which induced this phenotype at any time from 3 to 28 days of repeated administration. Principal component analysis using these probe sets clearly separated dose- and time-dependent clusters of the treated groups from their controls, except aspirin and propylthiouracil, both of which were considered to cause coagulopathy not due to their hepatotoxicity but due to their direct effects on the blood coagulation system. Reviewing the extracted genes, changes in lipid metabolism were found to be dominant. Genes related to blood coagulation were generally down-regulated by these drugs except that vitamin K epoxide reductase complex subunit 1 (Vkorc1) like 1, a paralogous gene of Vkorc1, was up-regulated. As expected, expression changes of these genes were least prominent in aspirin or propylthiouracil-treated liver. We concluded that these probe sets could be a good starting point in developing mechanism-based biomarkers for diagnosis or prognosis of hepatotoxicity-related coagulation abnormalities in the early stage of drug development.

A toxicogenomics approach for early assessment of potential non-genotoxic hepatocarcinogenicity of chemicals in rats.

For assessing carcinogenicity in animals, it is difficult and costly, an alternative strategy has been desired. We explored the possibility of applying a toxicogenomics approach by using comprehensive gene expression data in rat liver treated with various compounds. As prototypic non-genotoxic hepatocarcinogens, thioacetamide (TAA) and methapyrilene (MP) were selected and 349 commonly changed genes were extracted by statistical analysis. Taking both compounds as positive with six compounds, acetaminophen, aspirin, phenylbutazone, rifampicin, alpha-naphthylisothiocyanate, and amiodarone as negative, prediction analysis of microarray (PAM) was performed. By training and 10-fold cross validation, a classifier containing 112 probe sets that gave an overall success rate of 95% was obtained. The validity of the present discriminator was checked for 30 chemicals. The PAM score showed characteristic time-dependent increases by treatment with several non-genotoxic hepatocarcinogens, including TAA, MP, coumarin, ethionine and WY-14643, while almost all of the non-carcinogenic samples were correctly predicted. Measurement of hepatic glutathione content suggested that MP and TAA cause glutathione depletion followed by a protective increase, but the protective response is exhausted during repeated administration. Therefore, the presently obtained PAM classifier could predict potential non-genotoxic hepatocarcinogenesis within 24 h after single dose and the inevitable pseudo-positives could be eliminated by checking data of repeated administrations up to 28 days. Tests for carcinogenicity using rats takes at least 2 years, while the present work suggests the possibility of lowering the time to 28 days with high precision, at least for a category of non-genotoxic hepatocarcinogens causing oxidative stress.

Gene expression profiling of methapyrilene-induced hepatotoxicity in rat.

The present study was conducted as a model case of the toxicogenomics approach for analyzing toxicological mechanisms and toxicity assessments in the early stage of drug development by comparing with classical toxicology data. Methapyrilene (MP) 100 mg/kg produced obvious histopathological changes in liver of rats by single or repeated dose up to 28 days with significant elevation of ALT and AST. In the middle dose groups (30 mg/kg MP), no apparent changes were noted in blood biochemical data by single dosing or repeated dosing up to one week, and no obvious histopathological changes were observed except a slight hypertrophy in the hepatocytes. Comprehensive gene expression changes were analyzed using Affymetrix GeneChip and differentially expressed probe sets were statistically extracted. These contained many genes related to "glutathione metabolism", "apoptosis", "MAPK signaling pathway" and "regulation of cell cycle", which were all thought to be involved in the development of presently observed phenotypes. In the high dose groups, TGP1 scores (developed in our system in order to overview the responsiveness of drugs to multiple marker gene lists) for these categories were markedly increased from the early time point after single dose and kept their high expression throughout the repeated dose period. In the middle dose groups, the increment of the scores were noted not only at the time points when apparent pathological changes emerged, but also at the earlier stage of repeated dosing and even after single dosing. We conclude that toxicogenomics would enable a more sensitive assessment at the earlier time point than classical toxicology evaluation.

Gene expression profiling of rat liver treated with serum triglyceride-decreasing compounds.

We have constructed a large-scale transcriptome database of rat liver treated with various drugs. In an effort to identify a biomarker for interpretation of plasma triglyceride (TG) decrease, we extracted 218 probe sets of rat hepatic genes from data of 15 drugs that decreased the plasma TG level but differentially affected food consumption. Pathway and gene ontology analysis revealed that the genes belong to amino acid metabolism, lipid metabolism and xenobiotics metabolism. Principal component analysis (PCA) showed that 12 out of 15 compounds were separated in the direction of PC1, and these 12 were separated in the direction of PC2, according to their hepatic gene expression profiles. It was found that genes with either large or small eigenvector values in principal component PC 2 were those reported to be regulated by peroxisome proliferator-activated receptor (PPAR)alpha or constitutive androstane receptor (CAR), respectively. In fact, WY-14,643, clofibrate, gemfibrozil and benzbromarone, reported to be PPARalpha activators, distributed to the former, whereas propylthiouracil, omeprazole, phenobarbital, thioacetamide, methapyrilene, sulfasalazine and coumarin did to the latter. We conclude that these identified 218 probe sets could be a useful source of biomarkers for classification of plasma TG decrease, based on the mechanisms involving PPARalpha and CAR.

Identification of glutathione depletion-responsive genes using phorone-treated rat liver.

To identify candidate biomarker gene sets to evaluate the potential risk of chemical-induced glutathione depletion in livers, we conducted microarray analysis on rat livers administered with phorone (40, 120 and 400 mg/kg), a prototypical glutathione depletor. Hepatic glutathione content was measured and glutathione depletion-responsive gene probe sets (GSH probe sets) were identified using Affymetrix Rat Genome 230 2.0 GeneChip by the following procedure. First, probe sets, whose signal values were inversely correlated with hepatic glutathione content throughout the experimental period, were statistically identified. Next, probe sets, whose average signal values were greater than 1.5-fold compared to those of controls 3 hr after phorone treatment, were selected. Finally, probe sets without unique Entrez Gene ID were removed, ending up with 161 probe sets in total. The usefulness of the identified GSH probe sets was verified by a toxicogenomics database. It was shown that signal profiles of the GSH probe sets in rats treated with bromobenzene were strongly altered compared with other chemicals. Focusing on bromobenzene, time-course profiles of hepatic glutathione content and gene expression revealed that the change in gene expression profile was marked after the bromobenzene treatment, whereas hepatic glutathione content had recovered after initial acute depletion, suggesting that the gene expression profile did not reflect the hepatic glutathione content itself, but rather reflects a perturbation of glutathione homeostasis. The identified GSH probe sets would be useful for detecting glutathione-depleting risk of chemicals from microarray data.

Utilization of a one-dimensional score for surveying chemical-induced changes in expression levels of multiple biomarker gene sets using a large-scale toxicogenomics database.

A large-scale toxicogenomcis database has now been constructed in the Toxicogenomics Project in Japan (TGP). To facilitate the analytical procedures for such large-scale microarray data, we developed a simple one-dimensional score, named TGP1 which expresses the trend of the changes in expression of biomarker genes as a whole. To evaluate the usefulness of the TGP1 score, microarray data of rat liver and rat hepatocytes deposited in the TGP database were scored for three biomarker gene sets, i.e., carcinogenesis-related, PPARalpha-regulated and glutathione depletion-related gene sets. The TGP1 scoring system gave reasonable results, i.e., the scores for carcinogenesis-related genes were high in omeprazole-, chlorpromazine-, hexachlorobenzene-, sulfasalazine- and Wy-14,643-treated rat livers, that for PPARalpha-regulated genes were high in clofibrate-, Wy-14,643-, gemfibrozil-, benzbromarone- and aspirin-treated rat livers as well as rat hepatocytes, and for glutathione deficiency-related genes were high in omeprazole-, bromobenzene-, acetaminophen- and coumarin-treated rat liver. We concluded that the TGP1 score is useful for surveying the expression changes in multiple biomarker gene sets for a large-scale toxicogenomics database, which would reduce the time of doing conventional multivariate statistical analysis. In addition, the TGP1 score can be applied to screening of compatible biomarker gene sets between rat liver and rat hepatocytes, like PPARalpha-regulated gene sets, which will allow us to develop an appropriate in vitro system for drug safety assessment in vivo.