Collaborative studies in toxicogenomics in rodent liver in JEMS·MMS; a useful application of principal component analysis on toxicogenomics
Genes and Environment volume 38, Article number: 15 (2016)
Toxicogenomics is a rapidly developing discipline focused on the elucidation of the molecular and cellular effects of chemicals on biological systems. As a collaborative study group of Toxicogenomics/JEMS·MMS, we conducted studies on hepatocarcinogens in rodent liver in which 100 candidate marker genes were selected to discriminate genotoxic hepatocarcinogens from non-genotoxic hepatocarcinogens. Differential gene expression induced by 13 chemicals were examined using DNA microarray and quantitative real-time PCR (qPCR), including eight genotoxic hepatocarcinogens [o-aminoazotoluene, chrysene, dibenzo[a,l]pyrene, diethylnitrosamine (DEN), 7,12-dimethylbenz[a]anthracene, dimethylnitrosamine, dipropylnitrosamine and ethylnitrosourea (ENU)], four non-genotoxic hepatocarcinogens [carbon tetrachloride, di(2-ethylhexyl)phthalate (DEHP), phenobarbital and trichloroethylene] and a non-genotoxic non-hepatocarcinogen [ethanol]. Using qPCR, 30 key genes were extracted from mouse livers at 4 h and 28 days following dose-dependent gene expression alteration induced by DEN and ENU: the most significant changes in gene expression were observed at 4 h. Next, we selected key point times at 4 and 48 h from changes in time-dependent gene expression during the acute phase following administration of chrysene by qPCR. We successfully showed discrimination of eight genotoxic hepatocarcinogens [2-acetylaminofluorene, 2,4-diaminotoluene, diisopropanolnitrosamine, 4-dimethylaminoazobenzene, 4-(methylnitsosamino)-1-(3-pyridyl)-1-butanone, N-nitrosomorpholine, quinoline and urethane] from four non-genotoxic hepatocarcinogens [1,4-dichlorobenzene, dichlorodiphenyltrichloroethane, DEHP and furan] using qPCR and principal component analysis. Additionally, we successfully identified two rat genotoxic hepatocarcinogens [DEN and 2,6-dinitrotoluene] from a nongenotoxic-hepatocarcinogen [DEHP] and a non-genotoxic non-hepatocarcinogen [phenacetin] at 4 and 48 h. The subsequent gene pathway analysis by Ingenuity Pathway Analysis extracted the DNA damage response, resulting from the signal transduction of a p53-class mediator leading to the induction of apoptosis. The present review of these studies suggests that application of principal component analysis on the gene expression profile in rodent liver during the acute phase is useful to predict genotoxic hepatocarcinogens in comparison to non-genotoxic hepatocarcinogens and/or non-carcinogenic hepatotoxins.
Recently, a radical overhaul of toxicological test protocols has been proposed [1–4]. For example, Hartung wrote that after several productive decades, in which a patchwork of testing approaches was formed, fewer and fewer of the latest scientific development were incorporated . Caiment et al.  wrote that one of the main challenges of toxicology is the accurate prediction of compound carcinogenicity. The default test model for assessing chemical carcinogenicity, the 2-year rodent cancer bioassay, is currently criticized because of its limited specificity. With increased societal attention and new legislation against animal testing, toxicologists urgently need an alternative to the current rodent bioassays for chemical cancer risk assessment. In the beginning of the 21st century, toxicogenomics approaches proposed to use global high-throughput technologies (transcriptomics) to study the toxic effect of compounds on a biological system.
For risk assessment purposes, there is a general agreement that the chemicals acting through genotoxic and non-genotoxic mechanisms of carcinogenesis should be distinguished . Mathijs et al. hypothesized that genotoxic and non-genotoxic carcinogens induce distinct gene expression profiles, which may therefore be used to classify the mechanisms of compounds as either genotoxic carcinogens or non-genotoxic carcinogens . DNA microarray, which is a powerful technology for characterizing gene expression on a genome-wide scale , developed toxicogenomics. Quantitative real-time PCR (qPCR) is the field standard for measuring gene expression and is the most sensitive technique for the detection and quantification of mRNA target .
In the present study, we summarize our collaborative studies in toxicogenomics. We first selected about 100 candidate marker genes to discriminate mouse genotoxic hepatocarcinogens from non-genotoxic hepatocarcinogens by DNA microarrays, which were next quantified by qPCR . We extracted about 30 key genes from dose responses in gene expression  and selected key point times at the beginning and end of the acute phase (4 and 48 h) . We successfully showed the discrimination of genotoxic and non-genotoxic hepatocarcinogens in mouse liver  and rat liver  by qPCR and the application of principal component analysis (PCA) at 4 and 48 h after administration of hepatocarcinogens. The subsequent gene pathway analysis by Ingenuity Pathway Analysis extracted the DNA damage response, resulting from signal transduction by a p53-class mediator leading to the induction of apoptosis. Application of PCA was useful to discriminate genotoxic hepatocarcinogens from non-genotoxic and/or non-genotoxic non-hepatocarcinogens on rodent liver.
Selection of genes by DNA microarray and quantified by real-time PCR
In our preliminary study, we examined differential gene expression of 13 chemicals including eight genotoxic hepatocarcinogens [o-aminoazotoluene, chrysene, dibenzo[a,l]pyrene, DEN, 7,12-dimethylbenz[a]anthracene, dimethylnitrosamine, dipropylnitrosamine, and ENU], four non-genotoxic hepatocarcinogens [carbon tetrachloride, DEHP, phenobarbital, and trichloroethylene], and a non-genotoxic non-hepatocarcinogen (for mouse) [ethanol] using DNA microarray (Affymetrix GeneChip Mu74A V2 and in-house microarray) in mouse liver at 4 h and up to 28 days following a single intraperitoneal administration to groups of five 9-week-old B6C3F1 male mice. The cDNA was prepared with total RNA combined from pooled livers. After preliminary DNA microarray data were generated, results were confirmed by qPCR. We identified about 100 candidate genes to discriminate the genotoxic hepatocarcinogens from the non-genotoxic hepatocarcinogens. The results were published in part  and registered to the GEO database (GEO accession GSE33248). The changes in gene expression at 4 h were much greater than at 20 h, 14 days, and 28 days. We used qPCR in continual studies.
Dose-dependent alterations in gene expression at 4 h and 28 days
We examined the dose-dependent gene expression changes in candidate marker genes from our previous studies in mouse liver treated with two N-nitroso genotoxic hepatocarcinogens to extract key genes, and reported the results of 51 genes determined by qPCR . DEN at doses of 3, 9, 27, and 80 mg/kg body weight (bw) (LD50: 200 mg/kg bw, oral) or ENU at doses of 6, 17, 50, and 150 mg/kg bw (LD50: 200 mg/kg bw, intraperitoneally) were administered to groups of five 9-week-old B6C3F1 male mice, and the livers were dissected after 4 h and 28 days. Control mice received sterile water. The cDNA was prepared with total RNA from pooled livers and qPCR relative quantitative values were normalized using the Gapdh housekeeping gene. A total of 32 genes exhibited a dose response either via increased or decreased expression at least once at 4 or 48 h by DEN or ENU. At 4 h, as shown in Fig. 1 (Fig. 2 in ), 26 genes showed an obvious dose-dependent increase in gene expression by DEN [Aen (Isg20l1), Bax, Btg2, Ccng1, Ccng2, Cdkn1a, Cyp4a10, Cyp21a1, Fos, Gadd45b, Gdf15, Hmox1, Hspb1, Hspb2, Igfbp1, Jun, Mbd1, Mdm2, Myc, Net1, Plk2, Pmm1, Ppp1r3c, Rad52, Rcan1, and Tubb4b (Tubb2c)] from over 2-fold to 64-fold. Two genes [Cyp1a2 and Glul] showed a dose-dependent decrease in the DEN-treated sample at 4 h. ENU exhibited similar results except for a few genes (Fabp5 and Hist1h1c), although the increase in gene expression to ENU was generally weaker than to DEN. At 28 days, DEN induced a dose-dependent increase, between 2- and 4-fold, in four genes [Btg2, Cdkn1a, Cyp21a1, and Gdf15], and a dose-dependent decrease in Igfbp1 by less than 0.3-fold. ENU exhibited similar results except for the genes Casp1, Gstk1, Hspab1, and Ung. Only Gdf15 displayed a dose-dependent increase in expression on day 28 for both carcinogens. In addition, gene networks were analyzed using Ingenuity Pathway Analysis (IPA, http://www.ingenuity.com/products/ipa), a web-based software application for the analysis, integration, and interpretation of data derived from ‘omics experiments’ such as our qPCR data. Five gene networks were extracted by IPA: Network 1 consisted of genes related to cancer and cell-cycle arrest, such as Bax, Btg2, Ccng1, Cdkn1a, Gadd45b, Gdf15, Hspb1, Hspb2, Mdm2, Plk2, and Pmm1; Network 2 comprised cell cycle, DNA replication and recombination, repair, and cell death genes, such as Ccng2, Cyp1a2, Cyp4a10, Cyp21a1, Gdf15, Ppp1r3c, Rcan1, and Tubb4b (Tubb2c).
Time-course changes in gene expression at the acute stage within 48 h
We previously noticed that changes in gene expression were greater at 4 h, while reports on changes in the gene expression profile in rodent liver at the acute stage in the first 48 h after administration of a hepatocarcinogen were limited. We therefore selected key point times at 4 and 48 h from changes in time-dependent gene expression in mouse liver during the acute phase between 4 and 48 h after administration of chrysene, a polycyclic aromatic hydrocarbon (PAH) and genotoxic hepatocarcinogen, as determined by qPCR . Chrysene (100 mg/kg bw) was injected intraperitoneally into groups of three 9-week-old B6C3F1 male mice, and 4, 16, 20, 24, and 48 h later, livers were dissected and processed for gene expression. The cDNA was prepared with total RNA from each individual liver, and the amount of each gene was quantified by qPCR. We reported the results from 50 genes, 35 of which exhibited statistically significant increases at least once within 48 h after exposure to chrysene (Table 1). Fifteen genes [Bhlhe40, Btg2, Casp4, Ccng2, Cdkn1a, Crp, Cyp1a1, Cyp1a2, Fkbp5, Gadd45b, Gadd45g, Hmox1, Igfbp1, Lcn2, and Ly6a] at 4 h, six genes at 16 h, seven genes at 20 h, seven genes at 24 h, and 10 genes [Bhlhe40, Ccnf, Cyp1a1, Cyp1a2, Ephx1, Hhex, Hmox1, Rcan1, Tubb2a, and Tubb4b] at 48 h showed statistically significant increases of more than 2-fold. No significant decreases in gene expression were observed in this study. IPA at 4 h revealed that 7 genes [Btg2, Ccng2, Cdkn1a, Gadd45b, Gadd45g, Phlda3, and Mdm2] of 18 genes, which showed statistically significant increases, were associated with cancer, cell cycle, cell death and survival, and cellular growth and proliferation. The expression-increased genes from 16 to 48 h were associated with various biological processes including cancer. Cyp1a1 and Cyp1a2 showed remarkably consistent increases in gene expression during 4–48 h. These two genes are associated with toxin metabolism, the oxidation-reduction process, and the induction by carcinogenic polycyclic aromatic hydrocarbons as reported previously . We noticed that the greatest characteristic differences between 4 and 48 h were with 11 genes [Ly6a, Gadd45g, Igfbp1, Lcn2, Casp4, Cdkn1a, Btg2, Ccng2, Fkbp5, Crp, and Gadd45b], which differentially exhibited a statistically significant increase more than 2-fold at 4 h, and six genes [Tubb2a, Ephx1, Hhex, Ccnf, Rcan1, and Tubb4b] differentially showed a statistically significant increase more than 2-fold at 48 h.
Discrimination of genotoxic and non-genotoxic hepatocarcinogens at 4 and 48 h in mouse liver
We next successfully showed the discrimination of eight genotoxic hepatocarcinogens from four non-genotoxic hepatocarcinogens at 4 and 48 h in mouse liver by qPCR and statistical analysis using the Dunnett’s test, Welch’s t-test, and PCA . Eight genotoxic hepatocarcinogens, 2-acetylaminofuluorene (300 mg/kg bw), 2,4-diaminotoluene (200 mg/kg bw), diisopropanolnitrosamine (500 mg/kg bw), 4-dimethylaminoazobenzene (100 mg/kg bw), 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (250 mg/kg bw), N-nitrosomorpholine (32 mg/kg bw), quinoline (100 mg/kg bw), and urethane (1000 mg/kg bw) and four non-genotoxic hepatocarcinogens, 1,4-dichlorobenzene (1000 mg/kg bw), dichlorodiphenyltrichloroethane (50 mg/kg bw), DEHP (2000 mg/kg bw), and furan (30 mg/kg bw) were injected intraperitoneally into groups of five 9-week-old B6C3F1 males, livers were collected at 4 and 48 h later, and processed for gene expression. The cDNA was prepared with total RNA from each individual liver, and the gene expression was quantified by qPCR. Control mice received a solvent diluent, either saline or olive oil. We reported the results from 35 genes: 34 genes [Aen, Bax, Bhlhe40, Btg2, Ccnf, Ccng1, Cdkn1a, Cyp1a2, Ddit4, Ddit4l, Egfr, Ephx1, Gadd45b, Gdf15, HistH1, Hmox1, Hspb1, Igfbp1, Jun, Lrp1, Ly6a, Mbd1, Mdm2, Phlda3, Plk2, Pml, Pmm1, Ppp1r3c, Psma3, Rad52, Rcan1, St3gal5, Trp53, and Tubb4b (Tubb2c)] showed statistically significant changes in their gene expression, at least once at 4 h and/or 48 h, as computed by the Dunnett’s test using the Gapdh gene to normalize the data. The statistical significance between the genotoxic and non-genotoxic hepatocarcinogens for each gene was assessed by the Welch’s test at 4 and 48 h after chemical administration. Different sets of 17 genes [Aen, Bax, Btg2, Ccng1, Cdkn1a, Egfr, Gdf15, Hist1h1c, Jun, Lrp1, Mbd1, Mdm2, Phlda3, Plk2, Pml, Ppp1r3c, and Tubb4b (Tubb2c)] at 4 h, and 19 genes [Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Ddit4l, Ephx1, Gadd45b, Gdf15, Lrp1, Ly6a, Mdm2, Phlda3, Plk2, Pmm1, Ppp1r3c, St3gal5, and Tubb4b (Tubb2c)] at 48 h showed a statistical significance between the genotoxic and non-genotoxic hepatocarcinogens, as analyzed by the Welch’s test.
Fourteen genes [Aen, Bax, Cdkn1a, Mdmd2, Btg2, Ccng1, Ddit4, Gdf15, Hist1h1c, Hmox1, Hspb1, Phlda3, Plk2, and Pml] identified in this study have been reported to be directly associated with Trp53. Among these, 11 genes [Aen, Bax, Btg2, Ccng1, Cdkn1a, Gdf15, Hist1h1c, Mdm2, Phlda3, Plk2, and Pml] showed a statistical significance between the genotoxic and non-genotoxic hepatocarcinogens analyzed by the Welch’s t-test at 4 and/or 48 h. Seven major biological processes were extracted from the Gene Ontology analysis (Gene Ontology Consortium: geneontology.org), which were apoptosis, cell cycle and proliferation, DNA damage and repair, oncogenes, and tumor suppression. IPA suggested the DNA damage response pathway resulting from signal transduction by a p53-class mediator was likely leading to the induction of apoptosis. Although we did not observe a significant increase more than 2-fold in Trp53 expression, it was reported that after exposure to DNA-damaging agents, and other stress stimuli, p53 protein was stabilized and activated by a series of post-translational modifications that freed it from MDM2, a ubiquitination ligase responsible for its ubiquitination prior to proteasome degradation .
Discrimination of the gene expression profile between the genotoxic and nongenotoxic hepatocarcinogens was achieved by statistical analysis using PCA.
Useful application of PCA on gene expression profile to discriminate genotoxic and non-genotoxic hepatocarcinogens
We performed a statistical analysis using a logarithmic (log2) transformation of the data to stabilize the variance. PCA is a classic statistical procedure and is recently increasingly applied to biological data. PCA involves a mathematical procedure that transforms a number of possibly correlated variables into a smaller number of uncorrelated variables called “principal components”. The first principal component (PC1) accounts for as much of the variability in the data as possible, and each succeeding component accounts for as much as of the remaining variability as possible.
The mathematical formula of PC1 (z1) for 4 h is presented as the following:
where a1p is the eigenvector and x is the canonical logarithmic (log2)-transformed gene ratios (exp/cont). PCA was performed using the PCA programs in GeneSpringGX11.0.1 (Agilent Technologies, Santa Clara, CA, USA). Initially, PCA was applied to all 34 logarithmic (log2) transformed ratios (exp/cont), and subsequently tried with various candidate gene sets until the optimal discrimination was observed. The candidate genes were selected primarily using the Welch’s t-test from the results at 4 h, 48 h, and a combination of both times . PCA can be also performed using a free software R (https://cloud.r-project.org/).
We selected specific genes to obtain an optimal separation between genotoxic hepatocarcinogens and non-genotoxic hepatocarcinogens using PCA. Seven genes [Btg2, Ccnf, Ccng1, Lrp1, Mbd1, Phlda3, and Tubb4b (Tubb2c)] were used for PCA at 4 h, 12 genes [Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Gdf15, Lrp1, Mbd1, Phlda3, Plk2, and Tubb4b] at 48 h, and a combination of both time points (Fig. 2).
Differentiation between genotoxic and non-genotoxic hepatocarcinogens at 4 and 48 h in rat liver
Finally, we examined hepatocarcinogens in rat liver, and showed successful differentiation of two genotoxic hepatocarcinogens [DEN and 2,6-dinitrotoluene] from a nongenotoxic-hepatocarcinogen [DEHP], and a non-genotoxic non-hepatocarcinogen [phenacetin] at 4 and 48 h by qPCR and PCA . Candidate genes were selected from the data generated in mice. Two genotoxic hepatocarcinogens: DEN (12.5, 25, and 50 mg/kg bw) and 2,6-dinitrotoluene (125 and 250 mg/kb bw), a non-genotoxic hepatocarcinogen: DEHP (1000 and 2000 mg/kg bw), and a non-genotoxic non-hepatocarcinogen: phenacetin (500 and 1000 mg/kg bw) were examined in liver samples from groups of four 4-week-old F344 male mice at 4 and 48 h after a single oral administration of a chemical. Control rats received a solvent of sterile water or olive oil. The cDNA was prepared with total RNA from each individual liver. We reported results from 33 genes: 32 genes [Aen, Bax, Btg2, Ccnf, Ccng1, Cdkn1a, Cyp21a1, Cyp4a1, Ddit4l, Egfr, Ephx1, Gadd45b, Gadd45g, Gdf15, Hhex, Hmox1, Hspb1, Igfbp1, Jun, Lpp, Ly6al, Mdm2, Myc, Net1, Phlda3, Plk2, Pml, Pmm1, Rcan1, Tnf, Tp53, and Tubb4b (Tubb2c)] exhibited statistically significant changes in expression according to statistical analysis using the Williams’ test and the Dunnett’s test; and a normalized gene, Gapdh. The changes appeared to be greater at 4 h than at 48 h. Statistical analysis via PCA successfully differentiated the genotoxic hepatocarcinogens from the nongenotoxic hepatocarcinogen and non-genotoxic non-hepatocarcinogen at 4 h based on 16 genes [Ccnf, Ccng1, Cy4a10, Ddit4l, Egfr, Gadd45g, Gdf15, Hspb1, Igfbp1, Jun, Myc, Net1, Phlda3, Pml, Rcan1, and Tubb4b (Tubb2c)], and 48 h based on 10 genes [Aen, Ccng1, Cdkn1a, Cyp21a1, Cyp4a10, Gdf15, Igfbp1, Mdm2, Phlda3, and Pmm1] (Fig. 3). Eight major biological processes were extracted from a Gene Ontology analysis: apoptosis, cell cycle and proliferation, DNA damage and repair, oxidative stress, oncogenes, and tumor suppression. IPA suggested the DNA damage response, which signals through a Tp53-mediated pathway and leads to the induction of apoptosis: 24 genes are associated with Tp53 directly or indirectly (Fig. 4). This study showed that mouse candidate marker genes are applicable to rats for the differentiation of the genotoxic hepatocarcinogens from the non-genotoxic hepatocarcinogens examined in this study.
Recently, a new toxicogenomics tool for hepatocarcinogenicity evaluation of drug candidates in rodents (mainly rats) was reported: ToxDBScan (http://www.ra.cs.uni-tuebingen.de/software/ToxDBScan/) , which is a web tool offering a quick and easy similarity screening of new drug candidates against two large-scale public databases, which contain expression profiles for substances with known carcinogenic profiles: TG-GATEs (http://toxico.nibiohn.go.jp/english/)  and DrugMatrix (https://ntp.niehs.nih.gov/drugmatrix/) . TG-GATEs contains DNA microarray data on 170 chemicals, primarily medicinal compounds. DrugMatrix contains toxicogenomic profiles (DNA microarray data) for 638 different compounds. These compounds include US Food and Drug Administration-approved drugs, drugs approved in Europe and Japan, withdrawn drugs, drugs in preclinical and clinical studies, biochemical standards, and industrial and environmental toxicants. Although these large databases based on DNA microarrays were prepared, the number of published papers on toxicogenomics by DNA microarrays and qPCR in rodent liver or liver cells was not as expected.
Since its first application to toxicogenomics in 2003, PCA is a classic statistical technique that is recently increasingly applied to biological data. Previously, we successfully applied PCA to human lung cancer cell lines [19, 20]. Successful discrimination was performed in some toxicogenomics studies, such as hepatocarcinogens against non-carcinogens in rat liver , and carcinogenic PAHs against non-carcinogenic PAHs in HepG2 cells . However, the number of publications using PCA in toxicogenomics is still limited. We are now trying to apply this type of analysis on selected key genes to rodent liver gene expression profiles that have been described previously (unpublished).
Additionally, the involvement of next-generation sequencing (NGS) technology for the study of toxicogenomics is now being introduced [23–25]. Jiang et al. reported that NGS technologies, in comparison to microarray-based technologies, may overcome the current limitations, and are promising for the development of predictive models in the near future . Maslov et al.  suggested that the NGS era is well underway; new methods have been developed to directly analyze genetic material in a genome-wide manner with single nucleotide resolution. Moreover, there is no dependency on any particular gene or cell line, and the genetic material derived from any cell or tissue can be analyzed. This makes NGS-based mutagenicity assays particularly suitable for use in genetic toxicology. As toxicology continues to develop, we expect that testing methods will continue to change in concert with increased knowledge and understanding.
In the present review, we summarize our toxicogenomics collaborative studies. We selected and quantified by qPCR candidate marker genes to discriminate mouse genotoxic hepatocarcinogens from non-genotoxic hepatocarcinogens examined by DNA microarrays. We determined 30 key genes by dose responses in mouse liver gene expression induced by DEN and ENU at 4 h and 28 days, and extracted key times between 4 and 48 h from time-course studies during the acute phase induced by chrysene. Finally, we successfully showed the discrimination in mouse liver of eight genotoxic hepatocarcinogens [2-acetylaminofuluorene, 2,4-diaminotoluene, diisopropanolnitrosamine, 4-dimethylaminoazobenzene, 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone, N-nitrosomorpholine, quinoline, and urethane] from four non-genotoxic hepatocarcinogens [1,4-dichlorobenzene, dichlorodiphenyltrichloroethane, DEHP, and furan] and in rat liver two genotoxic hepatocarcinogens [diethylnitrosamine and 2,6-dinitrotoluene] from a non-genotoxic hepatocarcinogen [DEHP] and a non-genotoxic and non-hepatocarcinogen [phenacetin] determined by qPCR and PCA at 4 and 48 h after administration of chemicals. The subsequent gene pathway studies extracted the DNA damage response, resulting from signal transduction by a p53-class mediator leading to the induction of apoptosis. These studies suggest that application of PCA in the study of toxicogenomics is useful to discriminate genotoxic hepatocarcinogens from non-genotoxic hepatocarcinogens and/or non-hepatocarcinogens in rodent liver.
All animal experiments in the original papers [9–13] were conducted in accordance with the NIH Guide for Care and Use of Laboratory Animals and approved by the Animal Care and Use Committee at the Biosafety Research Center, Foods, Drugs, and Pesticides (applicant: NM, RD-07-005; approval: No. 07–066, in 2007) and the Animal Care and Use Committee of the Mitsubishi Chemical Medience Corp (applicant: HS; approval: 2007–0138).
ingenuity pathway analysis
the Japanese environmental mutagen society
mammalian mutagenicity study group
next generation sequencing
polycyclic aromatic hydrocarbon
principal component analysis
quantitative real-time PCR
Hartung T. Toxicology for the twenty-first century. Nature. 2009;460:208–12.
Ward JM. Value of rodent carcinogenesis bioassay. Toxicol Appl Pharmacol. 2008;226:212.
Knight A, Bailey J, Balcombe J. Which drugs cause cancer? For and against: cancer bioassays. BMJ. 2005;331:E389–91.
Caiment F, Tsamou M, Jennen D, Kleinjans J. Assessing compound carcinogenicity in vitro using connectivity mapping. Carcinogenesis. 2014;35:201–7.
Bol HM, Huici-Montagud A. Strategy of the scientific committee on occupational exposure limits (SCOEL) in the derivation of occupational exposure limits for carcinogens and mutagens. Arc Toxicol. 2008;82:61–4.
Mathijs K, Brauers KJ, Jennen DG, Boorsma A, van Herwijnen MH, Gottschalk RW, et al. Discrimination for genotoxic and nongenotoxic carcinogens by gene expression profiling in primary mouse hepatocytes improves with exposure time. Toxicol Sci. 2009;112:374–84.
Wu Z. A review of statistical methods for preprocessing oligonucleotide microarrays. Stat Methods Med Res. 2009;18:533–41.
Nardon E, Donada M, Bonin S, Dotti I, Stanta G. Higher random oligo concentration improves reverse transcription yield of cDNA from bioptic tissues and quantitative RT-PCR reliability. Exp Mol Pahtol. 2009;87:146–51.
Watanabe T, Tobe K, Nakachi Y, Kondoh Y, Nakajima M, Hamada S, et al. Differential gene expression induced by two N-nitroso carcinogens, phenobarbital and ethanol in mouse liver examined with oligonucleotide microarray and quantitative real-time PCR. Gene Env. 2007;29:115–27.
Watanabe T, Tanaka G, Hamada S, Namiki C, Suzuki T, Nakajima M, et al. Dose-dependent alterations in gene expression in mouse liver induced by diethylnitrosamine and ethylnitrosourea and determined by quantitative real-time PCR. Mutat Res. 2009;673:9–20.
Sakurai M, Watanabe T, Suzuki T, Furihata C. Time-course comparison of gene expression profiles induced by the genotoxic hepatocarcinogen, chrysene, in the mouse liver. Gene Env. 2014;36:54–64.
Watanabe T, Suzuki T, Natsume M, Nakajima M, Narumi K, Hamada S, et al. Discrimination of genotoxic and non-genotoxic hepatocarcinogens by statistical analysis based on gene expression profiling in the mouse liver as determined by quantitative real-time PCR. Mutat Res. 2012;747:164–75.
Suenaga K, Takasawa H, Watanabe T, Wako Y, Suzuki T, Hamada S, et al. Differential gene expression profiling between genotoxic and non-genotoxic hepatocarcinogens in young rat liver determined by quantitative real-time PCR and principal component analysis. Mutat Res. 2013;751:73–83.
Shimada T, Inoue K, Suzuki Y, Kawai T, Azuma E, Nakajima T, et al. Arylhydrocarbon receptor-dependent induction of liver and lung cytochromes P450 1A1, 1A2, and 1B1 by polycyclic aromatic hydrocarbons and polychlorinated biphenyls in genetically engineered C57BL/6 J mice. Carcinogenesis. 2002;23:1199–207.
Lavin MF, Gueven N. The complexity of p53 stabilization and activation. Cell Death Differ. 2006;13:941–50.
Römer M, Backert L, Eichner J, Zell A. ToxDBScan: large-scale similarity screening of toxicological databases for drug candidates. Int J Mol Sci. 2014;15:19037–55.
Igarashi Y, Nakatsu N, Yamashita T, Ono A, Ohno Y, Urushidani T, et al. Open TG-GATEs: a large-scale toxicogenomics database. Nucleic Acids Res. 2015;43(Database issue):D921–7.
Ganter B, Snyder RD, Halbert DN, Lee MD. Toxicogenomics in drug discovery and development: mechanistic analysis of compound/class-dependent effects using the DrugMatrix database. Pharmacogenomics. 2006;7:1025–44.
Oshima Y, Shinzawa H, Takenaka T, Furihata C, Sato H. Discrimination analysis of human lung cancer cells associated with histological type and malignancy using Raman spectroscopy. J Biomed Opt. 2010;15:017009.
Watanabe T, Miura T, Degawa Y, Fujita Y, Inoue M, et al. Comparison of lung cancer cell lines representing four histopathological subtypes with gene expression profiling using quantitative real-time PCR. Cancer Cell Int. 2010;10:2.
Uehara T, Minowa Y, Morikawa Y, Kondo C, Maruyama T, Kato I, et al. Prediction model of potential hepatocarcinogenicity of rat hepatocarcinogens using a large-scale toxicogenomics database. Toxicol Appl Pharmacol. 2011;255:297–306.
Song MK, Song M, Choi HS, Kim YJ, Park YK, et al. Identification of molecular signatures predicting the carcinogenicity of polycyclic aromatic hydrocarbons (PAHs). Toxicol Lett. 2012;212:18–28.
Jiang J, Wolters JE, van Breda SG, Kleinjans JC, de Kok TM. Development of novel tools for the in vitro investigation of drug-induced liver injury. Expert Opin Drug Metab Toxicol. 2015;11:1523–37.
Maslove AY, Quispe-Tintaya W, Gorbacheva T, White RR, Vijg J. High-throughput sequencing in mutation detection: a new generation of genotoxicity tests? Mutat Res. 2015;776:136–43.
Ning B, Su Z, Mei N, Hong H, Deng H, Shi L, et al. Toxicogenomics and cancer susceptibility: advances with next-generation sequencing. J Environ Sci Health C Environ Carcinog Ecotoxicol Rev. 2014;32:121–58.
We thank all participants in the present collaborative studies. Original works were partly supported by KAKENHI (18310047) (C. Furihata, T. Watanabe, and T. Suzuki) and a High-Tech Research Center Project of Private Universities with a matching fund subsidy from the Japanese Ministry of Education, Culture, Sports, Science and Technology (C Furihata).
The authors declare that they have no competing interests.
FC, ST, HS and NM proposed the present collaborative studies. FC, WT, ST, HS and NM participated in experiments. All authors participated in rearing and anatomy of animals. FC, WT and ST performed DNA microarray and qPCR experiments and analysis of results. FC and ST wrote the present review and all authors read and approved the final manuscript.
About this article
Cite this article
Furihata, C., Watanabe, T., Suzuki, T. et al. Collaborative studies in toxicogenomics in rodent liver in JEMS·MMS; a useful application of principal component analysis on toxicogenomics. Genes and Environ 38, 15 (2016). https://doi.org/10.1186/s41021-016-0041-0
- Rodent liver
- DNA microarray
- Quantitative real-time PCR
- Principal component analysis
- Gene network