Skip Navigation


ToxSci Advance Access originally published online on June 30, 2007
Toxicological Sciences 2007 99(1):58-69; doi:10.1093/toxsci/kfm168
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
99/1/58    most recent
kfm168v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Michel, C.
Right arrow Articles by Boitier, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Michel, C.
Right arrow Articles by Boitier, E.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2007. Published by Oxford University Press on behalf of the Society of Toxicology. All rights reserved. For Permissions, please email: journals.permissions@oxfordjournals.org

Diethylnitrosamine Initiation Does Not Alter Clofibric Acid–Induced Hepatocarcinogenesis in the Rat

Cecile Michel*, Chantal Desdouets{dagger}, Mohamed Slaoui*, Kevin Robert Isaacs{ddagger}, Ruth Angela Roberts§ and Eric Boitier*,1

* Department of Drug Safety Evaluation, sanofi aventis R&D, Centre de Recherche de Vitry/Alfortville-Evry, 94403 Vitry sur Seine, France {dagger} INSERM, U567, Paris, France {ddagger} 14 Rossett Park Road, Harrogate, North Yorkshire HG2 9NP, United Kingdom § AstraZeneca, Safety Assessment, Macclesfield, SK10 4TG Cheshire, United Kingdom

1 To whom correspondence should be addressed at sanofi aventis R&D, Centre de Recherche de Vitry/Alfortville-Evry, Drug Safety Evaluation, 13 quai Jules Guesde, 94403 Vitry sur Seine, France. Fax: +33-1-58-93-81-97. E-mail: Eric.Boitier{at}sanofi-aventis.com.

Received April 5, 2007; accepted May 22, 2007


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Clofibric acid (CLO) is a nongenotoxic hepatocarcinogen in rodents that causes altered hepatocellular foci and/or neoplasms. Initiation by DNA-damaging agents such as diethylnitrosamine (DEN) accelerates focus and tumor appearance and could therefore significantly contribute to shortening of the regulatory 2-year rodent carcinogenicity bioassays. However, it is crucial to evaluate the histological and molecular impact of initiation with DEN on hepatocarcinogenesis promoted by CLO. Male F344 rats were given a single nonnecrogenic injection of DEN (0 or 30 mg/kg) followed by Control diet or CLO (5000 ppm) in diet for up to 20 months. Histopathology and gene expression profiling were performed in liver tumors and surrounding nontumoral liver tissues. The molecular signature of DEN was characterized and its histopathological and immunohistopathological effects on focus and tumor types were also determined. Although foci and tumors appeared earlier in the DEN + CLO–treated group compared to the group treated with CLO alone, DEN had little impact on gene expression in nontumoral tissues since the gene expression profiles were highly similar between Control and DEN-treated rats, and DEN + CLO- and CLO-treated rats. Finally, tumors obtained from DEN + CLO and CLO-treated groups displayed highly correlated gene expression profiles (r > 0.83, independently of the time-point). The pathways involved in tumor development revealed by Gene Ontology functional analysis are similar when driven either by spontaneous initiation or by a chemically induced initiation step. Our work described here may contribute to the design optimization of shorter preclinical tests for the evaluation of the nongenotoxic hepatocarcinogenic potential of drugs under development.

Key Words: peroxisome proliferators; hepatocarcinogenesis; initiation; promotion; toxicogenomics; Clofibrate; diethylnitrosamine.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Nongenotoxic hepatocarcinogenesis is a slow, multistep process, which can be accelerated by chemical or physical initiation of liver cells (Scherer and Emmelot, 1975Go). Although the concept of initiation is well understood, its molecular drivers remain unclear and there is still debate on the relative contributions of chemically induced initiation, oncogene activation, tumor suppressor gene inactivation, and how these render cells more sensitive to tumor promoters. This lack of clear understanding probably underpins the persistent absence of accepted short-term tests for carcinogenic potential that could replace the 2-year rodent carcinogenicity studies required for the registration of pharmaceutical agents. Different exploratory approaches are available today for shortening this evaluation: models such as the medium-term rat liver bioassay (Ito et al., 1988Go) or the use of gene expression profiles for early biomarker determination (Michel et al., 2005Go). In this paper, we aim to clarify the influence of the DNA-damaging initiating agent diethylnitrosamine (DEN) on the development of nongenotoxic hepatocarcinogenesis, particularly in terms of focus development and of the molecular mechanisms involved in order to evaluate whether short-term carcinogenicity bioassays could be established as potential replacement in vivo models for the current regulatory prerequisites.

To characterize the histopathological and molecular influence of initiation on subsequent tumor promotion, Fisher F344 rats were given a single ip dose of DEN, an initiator that binds to DNA and induces mutations (den Engelse et al., 1983Go). DEN was given at a nonnecrogenic dose (30 mg/kg) (Scherer and Emmelot, 1975Go), therefore only exhibiting initiating properties. The glutathione-S-transferase-pi form (GST-p) is a highly sensitive liver focus biomarker that is overexpressed following DEN administration, allowing to evaluate the initiation efficiency (Moore et al., 1987Go). Clofibric acid (CLO), the primary metabolite of the hypolipidemic drug Clofibrate and its pharmacologically active form belongs to the class of peroxisome proliferators (PP), which have been shown to be nongenotoxic hepatocarcinogenic tumor promoters in the rodent (Mochizuki et al., 1982Go; Reddy and Qureshi, 1979Go) via the activation of the nuclear receptor/transcription factor peroxisome proliferator–activated receptor-{alpha} (Peters et al., 1997Go). This justifies the use of transcriptional profiling for studying the hepatocarcinogenic process induced by PP. As di(2-ethylhexyl)phthalate and Wyeth-14,643 but not Clofibrate were shown to significantly increase the mutation frequency in liver (Boerrigter, 2004Go), CLO was selected as a pure tumor promoter.

Together with the pathological evaluation, a global molecular snapshot of CLO-induced gene expression modulation was performed by gene expression profiling. Livers from control and DEN-treated rats presented highly correlated gene expression profiles, as were those obtained from DEN + CLO- or CLO-treated rats. Finally, tumors obtained from DEN + CLO and CLO-treated groups also exhibited molecular signatures that were strongly correlated in terms of genes modulated and functional categories involved. These results show that chemically induced initiation accelerates the molecular mechanisms of tumor promotion without modifying it.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Study design.
Details of the chemicals used and of the study design have been described elsewhere (Michel et al., 2005Go). Briefly, 175 Fisher F344 rats were divided into four groups: after 1 week on basal diet, rats belonging to the DEN and DEN + CLO groups underwent ip injection of 30 mg/kg DEN, those from the Control and CLO groups received NaCl 0.9%. Twelve days after injection, rats belonging to the DEN + CLO and CLO groups received a diet containing 5000 ppm of CLO for up to 12 or 20 months, respectively. DEN and Control rats were fed on normal diet for up to 12 or 20 months, respectively. Series of five rats from each group were necropsied on days 18, 46, 102, 264, 377, and, for the CLO and Control groups, also 524, and 607 days after the injection of saline. They were fasted approximately 16 h before being anesthetized by ip injection of Pentobarbital and culled by exsanguination. Livers were immediately excised under sterile conditions. They were examined macroscopically and visible tumors were excised and stored at – 80°C for further RNA extraction and downstream gene expression profiling. The number of tumors obtained per time-point is detailed in Figure 1. Portions of nontumoral liver from all animals were collected in 10% neutral buffer formalin and embedded in paraffin for further Hematoxylin & Eosin (H&E) staining, histopathological examination, and immunohistological evaluation. Other portions were flash-frozen in liquid nitrogen for total RNA extraction. Animal experiments were conducted in accord with accepted standards (local and national regulations) of humane animal care.


Figure 1
View larger version (24K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 1. Experimental design of the study summarizing the number and type of samples per time-point.

 
Evaluation of preneoplastic and neoplastic changes.
Liver foci, hepatocellular adenomas, and carcinomas were identified according to published criteria (Goodman et al., 1994Go; Mohr, 1997Go). Liver foci were distinguished from adenomas by the absence of pronounced compression of the surrounding tissue and retention of lobular architecture. Foci were classified as tigroid, eosinophilic, clear cell, or basophilic, according to the predominant staining characteristics of the hepatocytes. The incidence of preneoplastic lesions was recorded using a linear semiquantitative size grade from 1 to 5 for each replicate of the group. The number and size of each type of focus were used to calculate a quantifier (Table 1).


View this table:
[in this window]
[in a new window]

 
TABLE 1 Incidence of Each Type of Focus and of Neoplasm Determined Microscopically during a Long-Term CLO-Induced Hepatocarcinogenesis Study

 
Liver sections from animals euthanized on days 46 and 377 were used for immunohistochemical characterization of GST-p (Biotrin, Dublin, Ireland). Briefly, sections were deparaffinized and blocked for endogenous peroxidase activity. Sections were incubated with an anti-GST-p antibody (Dako France SAS, Trappes, France) diluted at 1:300 for 1 h at room temperature, then incubated with peroxidase-conjugated goat antirabbit immunoglobulin. Slides were then incubated with peroxidase antiperoxidase conjugate, reacted with 3',3-diaminobenzidine tetrahydrochloride and 0.2% hematoxylin counter-stained. GST-p–positive foci were counted and the number of cells per focus was evaluated. For each time-point, data were pooled for the five rats from each group.

Total RNA extraction and processing.
The standard Affymetrix protocol used has been described elsewhere (Michel et al., 2005Go).

Gene expression data analysis.
Reproducibility and similarity of the DNA chips were evaluated on normalized intensity data with the in-house developed analysis GECKO platform (Theilhaber et al., 2004Go). First, the 1000 genes that varied most between conditions were taken for a principal component analysis (PCA), and the first three dimensions (these account for most of the variance between treatment groups) were inspected graphically to verify that chips of the same treatment group clustered together. All microarray processing and downstream statistical analysis were performed in Rosetta Resolver version 4 (Rosetta Biosoftware, Seattle, WA). Data for each chip were automatically background-corrected via a chip-specific error model producing expression profiles. These individual scan data processing step were then followed by additional procedures to correct interchip variability performed via the standard analysis pipelines implemented in Resolver (normalization). Mean intensity ratio profiles were calculated for a treatment group over a baseline using the five replicates per group (except for the tumor samples).

In order to evaluate the influence of the DEN-induced initiation step on the hepatocarcinogenic process, ratio analyses were performed:

(1) For nontumor liver tissues, significantly modulated genes were selected by a two-factor ANOVA (treatment and time, p < 1.10–5) on at least three time-points for CLO versus Control, and DEN + CLO versus DEN. Significantly modulated genes obtained in the two analyses were compared and used for downstream functional pathway mapping (see below).
(2) For tumor tissues:
  • Significantly modulated genes were selected by a two-factor ANOVA (treatment and time, p < 1.10–5) on at least two time-points for CLO tumors versus CLO nontumors and DEN + CLO tumors versus DEN + CLO nontumors. Results of the two analyses were compared.
  • Functional pathway mapping of tumors was carried out on genes significantly modulated by DEN + CLO or CLO treatment on day 377 (early tumors, with similar adenoma phenotype).

(3) To understand the molecular mechanisms involved in DEN initiation, signature of DEN-modulated genes was determined by a two-factor ANOVA (treatment and time, p < 1.10–5) on at least three time-points for DEN versus Control.

Clustering was then performed on log-transformed ratios, using the union of the gene lists identified by the two ANOVA analyses, either for comparison 1 or comparison 2a by an agglomerative algorithm (Toronen et al., 1999Go), using a Euclidean distance metric, and visualized as heat-maps. Green shades indicate a transcriptional downregulation of a gene, and red an upregulation, whereas black represents no modulation. The intensity of the shades represents the magnitude of the change in gene expression (Chiang et al., 2001Go). Pearson's correlation coefficients were calculated on the basis of linear regression between pairs of mean intensity ratio profiles representative of the various study conditions.

Gene lists were analyzed for biological processes with an in-house Gene Ontology (GO) annotation tool (Muller et al., 2005Go). The complete data set can be accessed at the following web site: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE2216.

Quantitative real-time PCR.
Portion of the complementary DNA synthesized during the Affymetrix labeling process was diluted 1/150 into RNAse-free water. The following genes, most statistically modulated by DEN (p < 1.10–11), were tested relative to rat ß2-microglobulin: the chaperonin containing t complex polypeptide 1 (TCP1), epoxide hydrolase 1 (EPHX1), cyclin G1 (CCNG1), O(6)-methylguanine-DNA methyltransferase (MGMT), aldehyde dehydrogenase 1A1 (ALDH1A1), and p-glycoprotein 1 (MDR1). Their primers were purchased as Assays-on-Demand (respectively, Rn00562030, Rn00563349, Rn00563907, Rn00563462, Rn00755490, and Rn00561753, Applied Biosystems, Foster City, CA). Quantitative real-time polymerization chain reaction (qRT-PCR) was performed using the ABI 7700 Sequence Detection System (Applied Biosystems) according to the manufacturer's instructions. The Taqman Universal PCR Master Mix was used. Amplification reactions were carried out using the following temperature profile: 50°C, 2 min; 95°C, 10 min (95°C, 15 s; 60°C, 1 min) for 40 cycles. Fluorescence emission was detected for each PCR cycle and the threshold cycle (Ct) values were determined. Induction or repression of a gene in a DEN-treated sample relative to Control was calculated as follows:

Fold change = log2(2–(ct treated – ct ctrl) Gene of interest/ 2– (ct treated – ct ctrl) ß2-microglobulin). Values were reported as an average of duplicate analyses.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
DEN Accelerates the Appearance of Foci and Tumors
The results presented in Table 1 show the influence of DEN initiation on the kinetics of appearance of each type of focus and of neoplastic lesion. In Controls, tigroid foci were the first type to appear on day 264 of the study and their incidence increased with time. They represented the most abundant focus type in this group. DEN injection led to an earlier appearance of all types of focus, with tigroid foci appearing on day 102. In addition, clear cell and basophilic foci appeared on day 264, eosinophilic foci appeared on day 377, instead of day 524 in Controls, respectively. No tumors were observed in either the Control or DEN groups within the study course (607 days).

CLO treatment had no influence on tigroid foci, as they were seen on day 264 in both the CLO and Control groups with similar incidence. Clear cell and eosinophilic foci appeared 264 and 377 days, respectively, after the beginning of treatment with CLO. Basophilic foci appeared on day 524 in Control and CLO-treated groups but their size was increased after CLO treatment and they were the most prominent focus type on day 607.

The influence of DEN initiation on tigroid foci was inflected by CLO promotion as they appeared later and in fewer number in the DEN + CLO group compared to the DEN group, although more prominent than in the Control or CLO groups. DEN and CLO had a synergistic effect on eosinophilic, clear cell and basophilic foci as they appeared earlier and were more prominent in the DEN + CLO group compared to all the other groups. A higher incidence and an earlier onset of hepatocellular neoplasms were also noted in the DEN + CLO group as compared to the CLO group (day 264 vs. day 377).

Observation of the number and size of GST-p–positive foci enabled the evaluation of the efficiency of the initiation step and to monitor the fate of this neoplastic marker throughout the hepatocarcinogenic process (Fig. 2). Control animals developed a few GST-p–positive foci on day 377, likely corresponding to those observed by H&E staining. Both number and size of foci increased over time for all other groups. On day 46, 13 GST-p–positive foci were observed in DEN animals. Their number and size increased with time (see Fig. 2b). No GST-p–positive foci were observed on day 46 in the CLO and DEN + CLO groups, but 15 and 35 foci were counted on day 377, respectively. Figure 2a shows that some foci of the DEN + CLO group stained for GST-p, whereas others did not. There was no significant difference in terms of number of GST-p–positive foci between DEN and DEN + CLO groups on day 377.


Figure 2
View larger version (126K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 2. Incidence and size of GST-p–positive foci during CLO-induced hepatocarcinogenesis. (a) Paraffin liver section from a DEN + CLO animal on day 377 showing patterns of GST-p immunoreactivity. (b) Number and size of GST-p–positive foci per group on days 46 and 377. *The number of foci of DEN-injected groups on day 377 is statistically different from Control and CLO-treated groups at p < 0.01. However, there was no statistically significant group effect on the distribution concerning the size.

 
As a conclusion, in the two-step carcinogenesis protocol, DEN injection accelerated the appearance of foci and tumors without detectably modifying their types (though more GST-p–positive foci were found in the DEN + CLO group than after a pure CLO promotion).

DEN + CLO and CLO Groups Exhibit Highly Correlated Liver Gene Expression Patterns
We performed a general comparison of the liver gene expression profiles obtained from the four study groups. Figure 3 represents a qualitative comparison of the raw gene expression data for each replicate of the four groups from the nontumoral or tumoral liver tissues for the 46-day and 377-day time-points. The scatter plot obtained from the PCA shows that the gene expression profiles from Control replicates cluster by time-point, suggesting good data reproducibility, as for the other groups. Three main clusters could be distinguished in this scatter plot. Whatever the time-point, gene expression profiles from the Controls clustered together with the DEN-treated animals (green circle). Second, CLO and DEN + CLO–treated animals clustered together (orange circle). Third, tumor samples clustered away (red circle) from the nontumoral liver tissues and showed a broader spread, suggesting a higher variability between the various tumor samples. The similarity of the raw gene expression data between the DEN + CLO and the CLO groups was consistent throughout the study (data not shown). This similarity led us to compare in further details the gene expression modulations observed during the hepatocarcinogenic process after a spontaneous or a chemically induced initiation.


Figure 3
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 3. Qualitative comparison of the raw liver gene expression profiles obtained from the five replicates of the four treatment groups using PCA. Each symbol of the graph represents a time-point, color-coded as follows: Controls (blue), DEN-injected rats (green), DEN + CLO–treated rats (pink), or CLO-treated rats (yellow) for the 46- and 377-day time-points. CLO-free groups clustered together in a green circle, CLO-treated groups clustered together in an orange circle. All tumors clustered together in a red circle.

 
Genes Significantly Modulated in the DEN + CLO and the CLO Groups are Similar
This part focuses on identifying the effect of a chemically induced initiation on the CLO-induced transcriptional regulation during the hepatocarcinogenic process. Statistically modulated genes identified in the CLO-treated nontumoral groups (CLO vs. Control and DEN + CLO vs. DEN pairs) throughout the study were compared. In the CLO group (vs. Control), 495 genes were selected as significantly modulated, whereas 553 genes were identified in the DEN + CLO group (vs. DEN). Table 2 showed that 391 genes were in common between the two groups. The heat-map in Figure 4a showed that a majority of the differentially expressed genes exhibited a similar trend of modulation between the two CLO-treated groups. Since the GeneChip used contains 15,923 genes and the ANOVA p value cut-off was 1.10–5, less than one false positive was expected in each data set. This similarity was quantified by the calculation of Pearson correlation coefficients on the significantly modulated expression ratios obtained from the different sample comparisons: for a given time-point, the correlation coefficient was above 0.90 when comparing (CLO vs. Control) and (DEN + CLO vs. DEN) profiles, whereas this coefficient slightly decreased down to 0.71 when comparing different time-points. In order to place these complex lists of modulated genes in the context of underlying biological responses, we performed a statistical analysis of the significantly overexpressed (enriched) GO biological processes so that the molecular pathways affected could be identified. As shown in Figure 5a the unsupervised GO mapping analysis of the differentially regulated genes revealed that there were a majority of GO biological processes in common between the CLO and DEN + CLO nontumoral groups, namely lipid/fatty acid metabolism, fatty acid beta-oxidation, acyl-CoA metabolism, tricarboxylic acid cycle, mitochondrial electron transport, lipoprotein metabolism, lipid transport, steroid metabolism/biosynthesis, hormone metabolism, purine nucleotide metabolism, vitamin metabolism, amino acid metabolism/biosynthesis, sulfur metabolism, aromatic compound metabolism, translation, and complement activation. Four additional GO biological processes were unique to CLO-treated nontumoral liver tissues: hexose metabolism, regulation of cellular biosynthesis, amine biosynthesis, and heterocycle metabolism. DEN + CLO–treated nontumoral liver tissues exhibited two unique GO biological processes: tricarboxylic acid cycle/intermediate metabolism and cholesterol metabolism. However, all the GO biological processes could be clustered in four main functional categories: lipid and energy metabolism, nucleotide and vitamin metabolism, protein biosynthesis, and complement activation. When analyzed at this functional level, CLO and DEN + CLO groups displayed strictly identical molecular pathways.


View this table:
[in this window]
[in a new window]

 
TABLE 2 Treatment Effect on Nontumoral and Tumoral Liver Tissues

 

Figure 4
View larger version (61K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 4. Heat-map visualization of the differentially regulated genes and matrix of correlation coefficients among ratio profiles as a function of the treatment time and type, in nontumor (a) and tumor liver tissues (b): (a) DEN + CLO versus DEN and CLO versus Control. Two two-way ANOVAs were performed, using treatment ("CLO vs. Control" pair and "DEN + CLO vs. DEN" pair) and time as two factors, and statistically modulated genes were selected for a p value p < 1.10–5 on at least three time-points. (b) Tumors versus nontumors for DEN + CLO and CLO. Two two-way ANOVA were performed using tissue type (tumors vs. its surrounding tissue) and time as two factors, and statistically modulated genes were selected for a p value p < 1.10–5 on at least two time-points for the tumor samples. Briefly, mean intensity ratio profiles were then grouped and clustered hierarchically by heat-map using the union of the genes identified in the previous ANOVAs. Finally, mean Pearson correlation coefficients were calculated between ratio profiles representative of the various treatment groups and time-points and tabulated.

 

Figure 5
View larger version (50K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 5. GO-based pathway mapping of the genes significantly modulated either by CLO only or DEN + CLO (a) in nontumoral liver tissues and (b) in liver tumors. Significantly modulated genes observed in the microarray analyses (Table 2) were automatically analyzed for significantly overrepresented (enriched) GO biological processes. Given all genes on the microarray as the background population and a list of selected genes as the sample, the expected number of genes (as from a random sample) can be calculated for each GO term. The ratio of the observed over the expected number of genes in a term is the enrichment factor of that term in a gene list (sample). Enrichment greater than one corresponds to overrepresentation, and a p value for observing this enrichment by chance can be computed using a one-tailed Fisher exact test. For the selection of an enriched biological process within a set of genes, a p value ≤ 0.05 was required and a difference between the numbers of observed and expected genes of at least 4.

 
As shown in Figure 4a, DEN + CLO and CLO liver tumors also presented a similar trend of gene expression modulation throughout the study, although a higher intragroup data variability was noticed compared to the nontumoral liver tissues. When comparing profiles of statistically modulated genes between DEN + CLO tumors (vs. DEN + CLO nontumors) and CLO tumors (vs. CLO nontumors) at any time-point, correlations were bigger than 0.83 (Fig. 4b). Table 2 summarizes the number of genes in common between the CLO and DEN + CLO tumors (233), the number of genes unique to CLO tumors (218) and the number of genes unique to DEN + CLO tumors (601) on day 377. The statistical GO mapping analysis of the genes significantly modulated in the tumors revealed only two GO biological processes in common between DEN + CLO and CLO liver tumors: fatty acid metabolism, and nucleobase, nucleoside, nucleotide, and nucleic acid transport. The nine GO biological processes specific to CLO tumors were cholesterol metabolism, phospholipid metabolism, hexose metabolism, electron transport, nitrogen compound metabolism, RNA processing, translation, and protein import into nucleus. As for the DEN + CLO tumors, 15 GO biological processes were unique: response to drug, cellular carbohydrate metabolism, generation of precursor metabolites and energy, amino acid and derivative metabolism, amine biosynthesis, aromatic compound metabolism, rRNA processing, regulation of translation, translation initiation, protein folding, protein import, nuclear transport, nucleo-cytoplasmic transport, nucleotide metabolism, and DNA repair. When integrating the GO biological processes at a higher level, five main functional categories were identified, three of which were in common between CLO and DEN + CLO tumors: lipid and energy metabolism, protein biosynthesis, and nucleo-cytoplasmic trafficking and nucleotide metabolism. These three categories encompassed 91% of the GO biological processes identified across CLO and DEN + CLO tumors. The two remaining functional categories, namely response to drug and DNA repair, were unique to DEN + CLO tumors. Taken together, these results indicate that CLO and DEN + CLO tumors share similar treatment-impacted molecular pathways.

DEN Modulates the Expression of a Restricted Number of Genes
The similarity of the results obtained with our approach when comparing the effect of spontaneous and chemically induced initiation on CLO promotion raised the question of the sensitivity of our analysis to evaluate the impact of a single injection of DEN on transcription regulation. The six most-regulated genes in the DEN versus Control analysis were EPHX1, ALDH1A1, MGMT, MDR1, and cyclin G1: upregulated; TCP1: downregulated. Their modulation rapidly returned to baseline on day 18 (Fig. 6a). Their DEN-induced expression modulation was confirmed by qRT-PCR, except for TCP1 (Fig. 6b). Some of them were also modulated throughout CLO treatment (Fig. 6a): e.g., cyclin G1, MGMT, and MDR1 were slightly but increasingly downregulated throughout the carcinogenic process. ALDH1A1 was highly upregulated during the first 3 months of CLO treatment, following or not DEN injection.


Figure 6
View larger version (23K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 6. Heat-map visualization of the genes differentially regulated by DEN and qRT-PCR confirmation. (a) Time series of the expression of the six most statistically significant genes modulated by DEN (vs. Control) identified by two-way ANOVA (treatment and time as two factors) and monitored in the CLO-treated groups. (b) Confirmation of the expression of the selected genes in DEN samples (vs. Control) determined by qRT-PCR and normalized using ß2-microglobulin expression.

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
The majority of published data on the initiation–promotion paradigm are more than 10 years old. Although initiation mechanisms are well known, the events leading to the promotion of these initiated cells remain unclear and no molecular data have been provided to validate the concept. Oncogene activation or tumor suppressor gene inactivation was postulated to render cells more sensitive to tumor promoters. However, creation and use of transgenic mice for these genes has complicated the picture and most of these models are insensitive to nongenotoxic promoter chemicals (Gulezian et al., 2000Go). Here, we present an analysis of the histological and molecular liver changes occurring during preneoplastic and neoplastic events that take place when promotion is induced by CLO administration after an initiation step, either spontaneous or induced by the DNA-damaging agent DEN. We characterized the molecular signature of the initiation step and evaluated its influence on CLO-promoted tumors in terms of histological changes and molecular pathways modulated. This allowed us to show that, although DEN initiation favored the incidence and the early onset of hepatocellular foci and tumors in rats, it did not modify the molecular processes taking place in the liver, both in the nontumoral and tumoral tissues, under the hepatocarcinogenic effect of CLO.

Spontaneously initiated cells have been shown to be absent or rare in young Control Fisher 344 rats, but increase with age (Harada et al., 1989Go; Popp et al., 1985Go). We observed H&E and GST-p–stained foci, mainly tigroid of type, at later time-points in Controls. DEN increased the number of all focus type as compared to Controls, as described in the literature (Kraupp-Grasl et al., 1991Go), and specifically induced an earlier onset of eosinophilic and basophilic foci during CLO treatment, reinforcing the specific role of these foci during tumor development (Bannasch et al., 1985Go; Cattley and Popp, 1989Go; Grasl-Kraupp et al., 1993Go; Weber and Bannasch, 1994Go). Our results showing that basophilic foci appeared late upon CLO treatment were in line with what Kraupp-Grasl et al. (1990) reported in rats treated with nafenopin (another PP) for 13 months. Immunostaining for GST-p showed a significant increase of GST-p–positive foci after DEN initiation as previously described (Bralet et al., 2002Go). However, there was no statistically significant effect of CLO treatment on the number of GST-p–positive foci when comparing DEN + CLO and DEN groups, or CLO and Control groups, respectively. The few GST-p–positive foci found in the CLO group were in line with the findings of Rao et al. (1986)Go, who found that 10% and 48% of foci obtained after a treatment with 0.1% Wy-14,643 and 0.025% ciprofibrate, respectively, were GST-p positive. Unlike what was shown by some authors (Grasl-Kraupp et al., 1993Go), our results showed that CLO did not lead to the disappearance of GST-p–positive cells. However, the only tumor present in an immunohistochemically stained liver portion (from the DEN + CLO group) was GST-p negative. This would be in favor of the specific promotion of GST-p–negative lesions by PP.

The use of a low nonnecrogenic dose of DEN in our study could explain why our results were not in agreement with other publications that showed CLO treatment induced mainly basophilic and GST-p–negative foci (Cattley et al., 1989Go; Grasl-Kraupp et al., 1993Go; Weber and Bannasch, 1994Go). Williams et al. (2000)Go recently showed that the molecular mechanisms observed with low or high exposures to the DNA-reactive carcinogens, 2-acetylaminofluorene, and DEN drastically differed. Indeed, contrary to previous assumptions, the authors demonstrated the existence of thresholds for the hepatocellular initiating effects of these carcinogens. They also attributed the exaggerated responses at high exposures to cytotoxicity and compensatory hepatocyte proliferation. These papers point out the importance of the protocol used with regard to the mechanism, and to the development of different types of preneoplastic and neoplastic lesions in rodent models. In our study, the injection of DEN was nonnecrogenic and was not followed by partial hepatectomy. The lack of regenerative cell proliferation following initiation could explain the low amount and late appearance of basophilic foci obtained after CLO promotion in our study.

We observed that CLO alone was sufficient to induce foci and further led to the appearance of tumors. As demonstrated by others, the hepatocarcinogenicity of CLO alone could be due to the promotion of either spontaneously initiated cells (Cattley et al., 1991Go) or foci initiated by CLO itself, although all the evaluations of the initiating properties of PP previously tested in two-stage models turned out negative (Cattley and Popp, 1989Go; Glauert and Clark, 1989Go).

In addition to histopathology, we evaluated the molecular mechanisms underlying the hepatocarcinogenic process induced by CLO, focusing particularly on the influence of the initiation step on the promotion process. Indeed, we showed that genes selected as statistically modulated by CLO treatment were similar in the CLO and the DEN + CLO nontumoral groups, i.e., after a spontaneous or a chemically induced initiation. The corresponding molecular signatures were strongly correlated. When integrated at a biological level, the genes from both groups could be mapped to identical functional categories. The most enriched functional category was lipid and energy metabolism and corresponded to the specific signature of PP (genes such as peroxisomal bifunctional enzyme, peroxisomal 3-oxoacyl-coenzyme A thiolase, stearoyl-coenzyme A desaturase and fatty acid-binding protein, see supplementary data), therefore confirming previous microarray data (Cherkaoui-Malki et al., 2001Go; Yamazaki et al., 2002Go).

At the tumor level, the raw data were more variable within each group (Fig. 4). However, when focusing on the ratios, the gene expression profiles observed for the DEN + CLO and the CLO tumors were quantitatively similar (r > 0.9). This higher intragroup variability could be related to the drift of the tumor phenotype observed during the tumor progression, shifting from the adenoma type to the carcinoma type (see Fig. 1). The variability in the tumor phenotype together with the uneven number of tumors sampled and analyzed across the study decreased the power of the statistical analysis performed. Consequently, for the intergroup comparison of the molecular processes involved in tumor development, we decided to cross-analyze early tumors of the same phenotype (adenoma, see Table 1) and selected the day 377 time-point. Unlike what was observed with the nontumoral tissues, a statistical GO mapping analysis of the genes significantly modulated in the tumors revealed only a very few number of GO biological processes in common between DEN + CLO and CLO liver tumors. The higher number of genes found statistically modulated in the DEN + CLO tumors (Table 2) are in line with the higher number of statistically significant GO biological processes identified in the analysis when compared with those revealed in the CLO tumors. Here again, this difference may be due to the difference in the number of tumors used to generate the gene lists (3 and 2, respectively) impacting the statistical power of our analyses. However, when focusing on the functional categories, the three functional pathways in common between CLO and DEN + CLO tumors (lipid and energy metabolism, protein biosynthesis, and nucleo-cytoplasmic trafficking and nucleotide metabolism) encompassed 91% of the GO biological processes identified across CLO tumors, be after a spontaneous or a chemically induced initiation. Taken together, these functional categories indicated an increase in the protein turn-over and lipid metabolic state of the tumors, typical of these highly metabolically active tissues (Hruban, 1979Go; Seemann et al., 2006Go) and in line with what is described in the literature in response to PP treatment (Leonard et al., 2006Go). The two functional categories specific of the DEN + CLO tumors, namely, response to drug and DNA repair, led us to investigate the molecular effect of DEN and how the genes modulated could interfere with the CLO tumorigenic process.

Gene expression analysis linked to DEN injection led to the identification of six genes specifically and transiently modulated following DEN injection: EPHX1, ALDH1A1, MGMT, MDR1, cyclin G1, and TCP1. These genes, considered as the molecular signature of DEN initiation, have already been described in the literature to be modulated by alkylating agents such as DEN (except TCP1) within the scope of a 2-week toxicogenomics study in the rat up to 14 days after toxicogenomic evaluation (Ellinger-Ziegelbauer et al., 2004Go). DEN administration accelerated the appearance and increased the number of hepatocellular foci and tumors of rats treated with CLO (Cattley and Popp, 1989Go; Oesterle and Deml, 1983Go; Solt and Farber, 1976Go) but there was no evidence for long-term effect of DEN initiation on the gene expression level. Indeed, no DEN molecular effects could be identified when comparing DEN + CLO to CLO groups at the transcriptional level. Therefore, the mobilization of the DNA repair machinery in the tumors seems to be rather more an indication that the tumor cells undergo major gene/chromosome modifications as a function of the cancer progression and try to counteract these alterations.

This study demonstrates that initiation by a nonnecrogenic dose of DEN accelerates CLO-induced hepatocarcinogenesis, but without altering the histological profiles nor inducing any long-term gene expression changes in the focus. In particular, events linked to DEN-induced initiation and leading to the earlier appearance of carcinogenic lesions did not modify the long-term transcriptional modulations induced by CLO. These results need to be validated on microdissected foci, as it is probable that the effect of DEN is restricted to a few initiated cells. Gene expression profiling of microdissected basophilic foci (Michel et al., 2003Go) should be of particular interest to validate the molecular mechanisms of the hepatocarcinogenic process, to identify novel cancer early biomarkers and to confirm that nonnecrogenic doses of DEN can be used to boost the hepatocarcinogenic process without modifying the molecular events leading to CLO-induced cancer. Finally, our work could serve as a basis for the development of shorter-term bioassays for the evaluation of the nongenotoxic hepatocarcinogenic potential of drugs under development.


    SUPPLEMENTARY DATA
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Supplementary data are available online at http://toxsci.oxfordjournals.org/.


    ACKNOWLEDGMENTS
 
We thank A. Benevaut for performing the long-term in vivo study and the General Toxicology team for their support, particularly S. Trognon for preparing the diet. We thank Alexandre Secq for the histopathological preparations. S. Nicolas is thanked for statistical analysis. We warmly thank J. Theilhaber and A. Müller for assistance in data analysis.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 SUPPLEMENTARY DATA
 REFERENCES
 
Bannasch P, Benner U, Enzmann H, Hacker HJ. Tigroid cell foci and neoplastic nodules in the liver of rats treated with a single dose of aflatoxin B1. Carcinogenesis (1985) 6:1641–1648.[Abstract/Free Full Text]

Boerrigter ME. Mutagenicity of the peroxisome proliferators clofibrate, Wyeth 14,643 and di-2-ethylhexyl phthalate in the lacZ plasmid-based transgenic mouse mutation assay. J. Carcinog. (2004) 3:1–7.[CrossRef][Medline]

Bralet MP, Pichard V, Ferry N. Demonstration of direct lineage between hepatocytes and hepatocellular carcinoma in diethylnitrosamine-treated rats. Hepatology (2002) 36:623–630.[CrossRef][ISI][Medline]

Cattley RC, Marsman DS, Popp JA. Failure of the peroxisome proliferator WY-14,643 to initiate growth-selectable foci in rat liver. Toxicology (1989) 56:1–7.[CrossRef][ISI][Medline]

Cattley RC, Marsman DS, Popp JA. Age-related susceptibility to the carcinogenic effect of the peroxisome proliferator WY-14,643 in rat liver. Carcinogenesis (1991) 12:469–473.[Abstract/Free Full Text]

Cattley RC, Popp JA. Differences between the promoting activities of the peroxisome proliferator WY-14,643 and phenobarbital in rat liver. Cancer Res. (1989) 49:3246–3251.[Abstract/Free Full Text]

Cherkaoui-Malki M, Meyer K, Cao WQ, Latruffe N, Yeldandi AV, Rao MS, Bradfield CA, Reddy JK. Identification of novel peroxisome proliferator-activated receptor alpha (PPARalpha) target genes in mouse liver using cDNA microarray analysis. Gene Expr. (2001) 9:291–304.[ISI][Medline]

Chiang DY, Brown PO, Eisen MB. Visualizing association between genome sequences and gene expression data using genome-mean expression profiles. Bioinformatics (2001) 17:49S–55S.[Abstract]

den Engelse L, Floot BG, de Brij RJ, Tates AD. The induction of chromosomal damage in rat hepatocytes and lymphocytes. II. Alkylation damage and repair of rat-liver DNA after diethylnitrosamine, dimethylnitrosamine and ethyl methanesulphonate in relation to clastogenic effects. Mutat. Res. (1983) 107:153–166.[ISI][Medline]

Ellinger-Ziegelbauer H, Stuart B, Wahle B, Bomann W, Ahr HJ. Characteristic expression profiles induced by genotoxic carcinogens in rat liver. Toxicol. Sci. (2004) 77:19–34.[Abstract/Free Full Text]

Glauert HP, Clark TD. Lack of initiating activity of the peroxisome proliferator ciprofibrate in two-stage hepatocarcinogenesis. Cancer Lett. (1989) 44:95–100.[CrossRef][ISI][Medline]

Goodman DG, Maronpot RR, Newmberne PM, Popp JA, Squire RA. Proliferative and selected other lesions in the liver of Rats - Gl-5. In: In: Standard System of Nomenculature and Diagnostic Criteria Guides for Toxicologic Pathology (1994) Washington, DC: STP/ARP/AFIP.

Grasl-Kraupp B, Waldhor T, Huber W, Schulte-Hermann R. Glutathione S-transferase isoenzyme patterns in different subtypes of enzyme-altered rat liver foci treated with the peroxisome proliferator nafenopin or with phenobarbital. Carcinogenesis (1993) 14:2407–2412.[Abstract/Free Full Text]

Gulezian D, Jacobson-Kram D, McCullough CB, Olson H, Recio L, Robinson D, Storer R, Tennant R, Ward JM, Neumann DA. Use of transgenic animals for carcinogenicity testing: Considerations and implications for risk assessment. Toxicol. Pathol. (2000) 28:482–499.[Abstract/Free Full Text]

Harada T, Maronpot RR, Morris RW, Boorman GA. Observations on altered hepatocellular foci in National Toxicology Program two-year carcinogenicity studies in rats. Toxicol. Pathol. (1989) 17:690–706.[Medline]

Hruban Z. Ultrastructure of hepatocellular tumors. J. Toxicol. Environ. Health (1979) 5:403–433.[ISI][Medline]

Ito N, Tsuda H, Tatematsu M, Inoue T, Tagawa Y, Aoki T, Uwagawa S, Kagawa M, Ogiso T, Masui T, et al. Enhancing effect of various hepatocarcinogens on induction of preneoplastic glutathione S-transferase placental form positive foci in rats—An approach for a new medium-term bioassay system. Carcinogenesis (1988) 9:387–394.[Abstract/Free Full Text]

Kraupp-Grasl B, Huber W, Putz B, Gerbracht U, Schulte-Hermann R. Tumor promotion by the peroxisome proliferator nafenopin involving a certain subtype of altered foci in rat liver. Cancer Res. (1990) 50:3701–3708.[Abstract/Free Full Text]

Kraupp-Grasl B, Huber W, Taper H, Schulte-Hermann R. Increased susceptibility of aged rats to hepatocarcinogenesis by the peroxisome proliferator nafenopin and the possible involvement of altered liver foci occurring spontaneously. Cancer Res. (1991) 51:666–671.[Abstract/Free Full Text]

Leonard JF, Courcol M, Mariet C, Charbonnier A, Boitier E, Duchesne M, Parker F, Genet B, Supatto F, Roberts R, et al. Proteomic characterization of the effects of clofibrate on protein expression in rat liver. Proteomics (2006) 6:1915–1933.[CrossRef][ISI][Medline]

Michel C, Desdouets C, Sacre-Salem B, Gautier JC, Roberts R, Boitier E. Liver gene expression profiles of rats treated with clofibric acid: Comparison of whole liver and laser capture microdissected liver. Am. J. Pathol. (2003) 163:2191–2199.[Abstract/Free Full Text]

Michel C, Roberts RA, Desdouets C, Isaacs Kevin R, Boitier E. Characterization of an acute molecular marker of non-genotoxic hepatocarcinogenesis in rodent by gene expression profiling in a long term Clofibric acid study. Chem. Res. Toxicol. (2005) 18:611–618.[CrossRef][ISI][Medline]

Mochizuki Y, Furukawa K, Sawada N. Effects of various concentrations of ethyl-alpha-p-chlorophenoxyisobutyrate (clofibrate) on diethylnitrosamine-induced hepatic tumorigenesis in the rat. Carcinogenesis (1982) 3:1027–1029.[Abstract/Free Full Text]

Mohr UCE. (1997) International classification of rodent tumours, part 1. In The Rat, Fascicle 10. Digestive System. IARC Scientific Publications, No. 122. IARC, Lyon.

Moore MA, Nakagawa K, Satoh K, Ishikawa T, Sato K. Single GST-P positive liver cells—Putative initiated hepatocytes. Carcinogenesis (1987) 8:483–486.[Abstract/Free Full Text]

Muller A, Boitier E, Hu T, Carr GJ, Le Fevre AC, Marchandeau JP, Flor M, Jefferson F, Aardema MJ, Thybaud V. Laboratory variability does not preclude identification of biological functions impacted by hydroxyurea. Environ. Mol. Mutagen. (2005) 46:221–235.[CrossRef][ISI][Medline]

Oesterle D, Deml E. Promoting effect of polychlorinated biphenyls on development of enzyme-altered islands in livers of weanling and adult rats. J. Cancer Res. Clin. Oncol. (1983) 105:141–147.[CrossRef][ISI][Medline]

Peters JM, Cattley RC, Gonzalez FJ. Role of PPAR alpha in the mechanism of action of the nongenotoxic carcinogen and peroxisome proliferator Wy-14,643. Carcinogenesis (1997) 18:2029–2033.[Abstract/Free Full Text]

Popp JA, Scortichini BH, Garvey LK. Quantitative evaluation of hepatic foci of cellular alteration occurring spontaneously in Fischer-344 rats. Fundam. Appl. Toxicol. (1985) 5:314–319.[CrossRef][ISI][Medline]

Rao MS, Tatematsu M, Subbarao V, Ito N, Reddy JK. Analysis of peroxisome proliferator-induced preneoplastic and neoplastic lesions of rat liver for placental form of glutathione S-transferase and gamma-glutamyltranspeptidase. Cancer Res. (1986) 46:5287–5290.[Abstract/Free Full Text]

Reddy JK, Qureshi SA. Tumorigenicity of the hypolipidaemic peroxisome proliferator ethyl-alpha-p-chlorophenoxyisobutyrate (clofibrate) in rats. Br. J. Cancer (1979) 40:476–482.[ISI][Medline]

Scherer E, Emmelot P. Kinetics of induction and growth of precancerous liver-cell foci, and liver tumour formation by diethylnitrosamine in the rat. Eur. J. Cancer (1975) 11:689–696.[ISI][Medline]

Seemann MD, Beck R, Ziegler S. In vivo tumor imaging in mice using a state-of-the-art clinical PET/CT in comparison with a small animal PET and a small animal CT. Technol. Cancer Res. Treat. (2006) 5:537–542.[ISI][Medline]

Solt D, Farber E. New principle for the analysis of chemical carcinogenesis. Nature (1976) 236:702–703.

Theilhaber J, Malanthara A, Ulyanov A, Cole J, Xu D, Nahf R, Heuer M, Brockel C, Bushnell S. Gecko: A complete large-scale gene expression analysis platform. BMC Bioinformatics (2004) 5:195.[CrossRef][Medline]

Toronen P, Kolehmainen M, Wong G, Castren E. Analysis of gene expression data using self-organizing maps. FEBS Lett. (1999) 451:142–146.[CrossRef][ISI][Medline]

Weber E, Bannasch P. Dose and time dependence of the cellular phenotype in rat hepatic preneoplasia and neoplasia induced by single oral exposures to N-nitrosomorpholine. Carcinogenesis (1994) 15:1219–1226.[Abstract/Free Full Text]

Williams GM, Iatropoulos MJ, Jeffrey AM. Mechanistic basis for nonlinearities and thresholds in rat liver carcinogenesis by the DNA-reactive carcinogens 2-acetylaminofluorene and diethylnitrosamine. Toxicol. Pathol. (2000) 28:388–395.[Abstract/Free Full Text]

Yamazaki K, Kuromitsu J, Tanaka I. Microarray analysis of gene expression changes in mouse liver induced by peroxisome proliferator-activated receptor alpha agonists. Biochem. Biophys. Res. Commun. (2002) 290:1114–1122.[CrossRef][ISI][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
99/1/58    most recent
kfm168v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Right arrow Disclaimer
Google Scholar
Right arrow Articles by Michel, C.
Right arrow Articles by Boitier, E.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Michel, C.
Right arrow Articles by Boitier, E.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?