Differential gene expression in mouse primary hepatocytes exposed to the peroxisome proliferator-activated receptor α agonists
BMC Bioinformatics volume 7, Article number: S18 (2006)
Fibrates are a unique hypolipidemic drugs that lower plasma triglyceride and cholesterol levels through their action as peroxisome proliferator-activated receptor alpha (PPARα) agonists. The activation of PPARα leads to a cascade of events that result in the pharmacological (hypolipidemic) and adverse (carcinogenic) effects in rodent liver.
To understand the molecular mechanisms responsible for the pleiotropic effects of PPARα agonists, we treated mouse primary hepatocytes with three PPARα agonists (bezafibrate, fenofibrate, and WY-14,643) at multiple concentrations (0, 10, 30, and 100 μM) for 24 hours. When primary hepatocytes were exposed to these agents, transactivation of PPARα was elevated as measured by luciferase assay. Global gene expression profiles in response to PPARα agonists were obtained by microarray analysis. Among differentially expressed genes (DEGs), there were 4, 8, and 21 genes commonly regulated by bezafibrate, fenofibrate, and WY-14,643 treatments across 3 doses, respectively, in a dose-dependent manner. Treatments with 100 μM of bezafibrate, fenofibrate, and WY-14,643 resulted in 151, 149, and 145 genes altered, respectively. Among them, 121 genes were commonly regulated by at least two drugs. Many genes are involved in fatty acid metabolism including oxidative reaction. Some of the gene changes were associated with production of reactive oxygen species, cell proliferation of peroxisomes, and hepatic disorders. In addition, 11 genes related to the development of liver cancer were observed.
Our results suggest that treatment of PPARα agonists results in the production of oxidative stress and increased peroxisome proliferation, thus providing a better understanding of mechanisms underlying PPARα agonist-induced hepatic disorders and hepatocarcinomas.
Peroxisome proliferators are structurally diverse chemicals that include industrial pollutants, plasticizers, herbicides, and lipid-lowering drugs. Fibrates including bezafibrate, clofibrate, fenofibrate, WY-14,643, and others are a unique class of hypolipidemic drugs. They function as agonists for peroxisome proliferator-activated receptor alpha (PPARα). PPARα is a transcriptional nuclear receptor and forms a heterodimer with another nuclear receptor, retinoid X receptor (RXR). PPARα/RXR heterodimer binds effectively to the peroxisome proliferator response elements (PPREs) located in promoters of various target genes and regulate the expression of genes involved in lipid metabolism and peroxisome proliferation [1–3]. PPRE consists of direct repeats of TGA/TCCT which is separated by a single nucleotide (DR1) . Fibrates reduce plasma triglyceride and cholesterol levels via the activation of PPARα, which is considered to be the result of induction of fatty acid catabolism in the liver. At the molecular level, fibrates bind to PPARα and increase the expression of genes that involved in fatty acid uptake (fatty acid binding protein, FABP), β-oxidation (acyl-CoA oxidase, ACOX), and ω-oxidation (cytochrome P450) [4–7]. This pharmacological effect of fibrates is responsible for the therapeutic utility, and this effect was observed in preclinical species and also in humans.
Along with the pharmacological effects of fibrates, toxic effects such as marked peroxisome proliferation, hepatomegaly and hepatocarcinoma are observed in rodents . It is accepted that the mode of action (hepatocarcinogenesis) is dependent upon sustained PPARα activation. This mode of action is supported by the observation that even a one year exposure to PPARα agonists was insufficient to cause an increase in the incidence of hepatic neoplasms in PPARα knock-out mice. In addition, peroxisome proliferation and gene expression regulated by PPARα were not remarkably altered. One hypothesis for the carcinogenic mechanism of action of PPARα agonists in rodent liver is based on their ability to elevate peroxisomal β-oxidation and microsomal ω-oxidation of fatty acids, resulting in the generation of hydrogen peroxide. This excess production of hydrogen peroxide results in the generation of reactive oxygen species (ROS) and oxidative stress . The induction of oxidative stress has been suggested as a common pathway for many non-genotoxic carcinogens to elicit their carcinogenicity . In addition, increased peroxisome proliferation in response to activation by PPARα agonists is associated with tumor formation in rodent liver . The combined effect of increased oxidative stress and increased cell proliferation in the rodents exposed to PPARα agonists likely underlies their carcinogenic potential. The precise mechanism of the hepatocarcinogenesis of PPARα agonists in rodents is not fully understood. Since a number of fibrates (e.g., bezafibrate and fenofibrate) are used as therapeutic agents, it is important to analyze the mechanism of liver toxic effects occurred in rodents so that we can better evaluate the safety of these drugs.
Microarray technology provides a comprehensive, rapid and efficient method for large scale profiling of gene expression changes in biological samples (e.g., treatment versus control, disease versus normal). The advantages of DNA microarray technology include the ability to analyze expression patterns of thousands of genes simultaneously. Other advantages include the ability to characterize relationships between genes and the changes in biological processes such as disease states, developmental stages and responses to drugs [10, 11]. This method has been successfully employed in identifying gene expression changes in cells, including both hepatic cell lines  and isolated hepatocytes , in response to various stimuli.
In this study, using mouse primary hepatocytes, we examined global gene expression profiles observed after treatment with several concentrations of three PPARα agonists (i.e. bezafibrate, fenofibrate and WY-14,643). The gene expression profiles showed increased expression of genes involved in fatty acid oxidation and metabolism as expected, and also showed regulation of many genes involved in the production of ROS and those that are associated with liver cancer development. This study also demonstrates the similarity of gene expression changes induced by three different PPARα agonists. In addition, this whole genome microarray analysis performed following in vitro administration of PPARα agonists indicates a plausible mechanism of hepatocarcinogenesis in the mouse liver and may help with the safety assessment of this class of agents.
PPARα agonist administration to primary hepatocytes to assess PPARα activity
In order to estimate whether PPARα activity was elevated by the addition of PPARα agonists, pHD(x3) luciferase plasmid containing three direct tandem copies of a PPRE binding site was used as a reporter plasmid. Primary hepatocytes were co-transfected with pHD(x3) luc reporter plasmid and pSG5-PPARα or pSG5-PPARα/pSG5-RXR expression vectors. Twenty-four hours after transfection, cells were treated with PPARα agonist at various concentrations (10–100 μM) for 24 h. Basal level of luciferase activity was observed (data not shown) when cells were transfected with pHD(x3) luc reporter plasmid and pGS5 empty vector. This could be due to the activity of endogenous PPARα agonists (i.e., fatty acid and their metabolites) [5, 14]. Figure 1 shows that mouse primary hepatocytes transfected with pHD(x3) luc reporter plasmid and pGS5 empty vector produced 1.2, 0.8, and 1.8 fold increased in luciferase activity after the treatment with 100 μM bezafibrate, fenofibrate, and WY-14,643, respectively. When cells were treated with the same concentration of drugs in the presence of PPARα expression plasmid, the induction of luciferase activities increased to 2.3, 3.4 and 5.0 fold. Moreover, in the presence of both PPARα and RXR expression plasmids, luciferase activities increased even further to 6.6, 7.1, and 9.9 fold. PPARα activation inductions by these three agonists demonstrate a does-dependent increase for 100 μM compare wtih 10 and 30 μM treatments (Figure 1).
Gene expression patterns associated with exposures of PPARα agonists
To understand the molecular mechanisms responsible for the pleiotropic effects of fibrates, mouse primary hepatocytes were treated with three PPARα agonists (Figure 1). RNA was isolated following 24 h treatment and the Agilent Whole Mouse Genome Microarray analysis was performed. Three replicate arrays corresponding to each treatment (10, 30 and 100 μM) of each drug (bezafibrate, fenofibrate and WY-14,643) were compared to the control replicates (DMSO) using Student t-test. A gene was considered to be significantly regulated by a drug if the fold change was greater than 1.5 and the P-value was less than 0.05. Based on these two criteria, there were 4, 26, and 151 genes that showed an altered expression by the treatments with 10, 30 and 100 μM bezafibrate, respectively. There were 9, 41 and 149 genes altered by 10, 30 and 100 μM fenofibrate, and 31, 52 and 145 genes altered by the treatment with 10, 30 and 100 μM WY-14,643. The Venn diagrams (Figure 2) represent the numbers of differentially expressed genes (DEGs) from three drug treatments at various concentrations. Among these DEGs, 4, 8, and 21 genes were altered in common by bezafibrate, fenofibrate and WY-14,643 treatments across 3 doses, respectively (Figure 2). These commonly regulated genes showed clear dose-dependent changes (Figure 3). Except for the down regulation of Ahsg by WY-14,643, all genes were identified as up-regulated. Expression of Pdk4 and Cte1 showed a particularly prominent induction (10–20 fold at 100 μM) by treatment of WY-14,643 (Fig. 3C), as did Fabp1 with 20–30 fold induction by 100 μM treatments of fenofibrate and WY-14,643 (Figs. 3B &3C).
Hierarchical clustering analysis with DEGs revealed that these genes were grouped together based on treatment doses rather than on specific drugs (Figure 4), indicating that a class effect was detectable. Taking into consideration the fact that most of the genes altered by low (10 μM) and middle (30 μM) dose treatments were also altered by high (100 μM) dose treatments, further analysis was focused on the genes regulated by the high dose treatments. Figure 5 represents the numbers of genes regulated by the 100 μM treatments of three drugs and the numbers of overlapping genes among three drug treatments. Treatments with 100 μM of bezafibrate, fenofibrate and WY-14,643 resulted in 151, 149 and 145 genes with altered expression, respectively. Among them, 61 genes were concordantly regulated by three drugs and 121 genes were regulated commonly by at least 2 drugs. The correlation of log2 fold changes based on 61 genes which were regulated commonly by three drugs was determined by pairwise comparisons (Table 1). Comparison of gene expression profiles resulted in relatively high correlations (0.94–0.97) between three drug treatments, indicating that gene expression patterns are identical at the 100 μM concentration despite any differences among the drugs.
Gene function analysis
Using Ingenuity Pathway Analysis, we conducted gene function analysis with these 121 commonly regulated genes. Not surprisingly, the largest categories of induced genes were those involved in lipid metabolism (49 genes), including oxidation, modification, and metabolism. Table 2 shows the genes involved in oxidation (19 genes), of which 10 were involved in β-oxidation. As expected, many genes directly regulated by PPARα including Cpt1 (carnitine palmitoyltransferase 1), Fabp1 (fatty acid binding protein 1), Acox1 (acyl-Coenzyme A oxidase 1, palmitoyl), and Ehhadh (enoyl-Coenzyme A, hydratase/3-hydroxyacyl Coenzyme A dehydrogenase) were identified in this study. Induction of these genes corroborates previous data, and also serves as the validation of our microarray experiments. Several novel genes including Acacb, Acsl1, Ech1, Hadhb, and Pdk4 that are involved in oxidation of lipid metabolism were responsive to the treatments of PPARα agonists. In addition, four genes (Aldh3a2, Apoc2, Cd36, and Slc25a10) associated with the production of ROS exhibited 1.5–3.0 fold up-regulation, and two genes (Acox1 and Pex11a) involved in the cellular proliferation of peroxisomes were 3-fold up-regulated.
Administration of fibrates to rodents results in hepatic diseases and hepatocarcinoma . For gene function analysis, we also concentrated on genes involved in these effects. There were six genes related to hepatic disorders (Table 3), and all were up-regulated by fibrate administration. Four of these genes are involved in oxidation of lipids. In addition, 11 genes classified as being related to the development of liver cancer had altered expression levels. These include liver fatty acid binding protein 1 (Fabp1), lymphocyte antigen 6 complex locus D (Ly6d), monoglyceride lipase (Mgll), and angiopoietin-like 4 (Angptl4), with a wide range of up-regulation (1.6–34.4 fold) after treatments with these three PPARα agonists.
Fibrates, members of peroxisome proliferators and agonists of PPARα, are used to treat hyperlipidemia by reducing plasma triglyceride and cholesterol levels via accelerating lipid metabolism. In rodents, administration of fibrates can induce hepatomegaly and hepatocarcinoma, possibly due to the induction of cell proliferation and increased oxidative stress. Examination of PPARα-deficient mice demonstrated that the activation of PPARα is required exclusively for mediating both pharmacological (hypolipidemic) and toxic (carcinogenic) responses of fibrate administration [15, 16]. However, mechanisms of fibrate-induced hepatocarcinoma development and the potential risk of use of these drugs to humans remain unclear. Examination of gene expression profiles is an important approach that may help us better understand PPARα-mediated pleiotropic effects.
In this study, microarray analysis was applied to generate a molecular portrait of gene expression in mouse primary hepatocytes exposed to fibrates (Figure 1). We treated mouse primary hepatocytes with three fibrates (bezafibrate, fenofibrate and WY-14,643) at multiple doses (0, 10, 30, and 100 μM). Although global gene analysis study was conducted in vitro , the design of this study (i.e., treatment with multiple fibrates at multiple level doses) permitted us to detect whether changes in gene expressions are a class effect (e.g., genes are commonly regulated by multiple drugs) and whether changes are dose-dependent as well. Indeed, the majority of genes regulated by low and middle doses were also identified in high dose treatments. For example, 4/4, 9/9, and 25/31 genes that were regulated by 10 μM treatments of bezafibrate, fenofibrate and WY-14,643, respectively, were also found to be regulated in 100 μM treatments (Figure 2). In addition, dose-response dependency in gene expressions was also observed for the genes commonly regulated at multiple doses (Figure 3). The dose-dependent expression levels of genes altered by PPARα agonists allowed us to assess biological activity of this class of agents.
PPARα agonists have a therapeutic role in the management of fatty acid metabolism through their effects on β-oxidation and lipid transport. The gene expression changes in common across PPARα agonists may indicate those genes are directly regulated by PPARα stimulation. In this study, we demonstrated that 121 DEGs were altered in common by at least two of the three PPARα agonists tested. The Ingenuity Pathway Analysis was used to analyze gene functions and to provide pathway annotations. Based on this analysis, many of these genes (49 genes) are involved in the oxidation of fatty acids (Table 2) as has been previously shown for this class of agents [4, 18]. For example, acyl-coA synthetase catalyzes the precursor step to β-oxidation (ligates CoA to a free fatty acid) and three members of the long chain acyl CoA synthestase family (Acsl1, Acsl4, and Acsl5) were increased. This observation is supported by the work of Schoonjans et al. who demonstrated that the expression of Acs is altered by fibrates and that there is a PPRE in the Acsl promoter . These findings also agree with those of Cornwall et al. who reported that the expression of Acsl was elevated in the liver of rats exposed to fenofibrate . The effects on the β-oxidation pathway also include the induction of the first enzyme of peroxisomal β-oxidation, acyl-CoA oxidase (Acox), as well as the next enzyme in the cascade, enoyl-CoA hydratase/3-hydroxyacyl-CoA dehydrogenase (Ehhadh). The identification of a large number of lipid metabolizing genes following exposure to several PPARα agonists is in concordance with the known biochemical and molecular effects of these hypolipidemic agents to regulate lipid metabolism.
PPARα agonists are also considered to be nongenotoxic carcinogens in rodents. Oxidative stress has been proposed as a common pathway for many non-genotoxic carcinogens . In the present study, 10 genes involved in fatty acid β-oxidation were up-regulated upon exposure to PPARα agonists (Table 2), which included Acox, the key enzyme of peroxisomal fatty acid β-oxidation system. The elevation of peroxisomal fatty acid β-oxidation such as occurs with PPARα agonist exposure in rodents results in the elevated generation of hydrogen peroxide . Substantial production of hydrogen peroxide causes oxidative stress and the induction of ROS. The increased ROS associated with elevated levels of Acox has been postulated to mediate the hepatocarcinogenesis resulting from PPARα exposure in rodents. We observed four genes (Aldh3a2, Apoc2, Cd36, and Slc25a10) associated with the production of ROS were up-regulated. A growing body of evidence indicates that Cd36 (CD36 antigen) is involved in the cytotoxicity associated with inflammation and is an essential mediator of the production of ROS . In addition, six genes classified as related to hepatic disorders were identified as being up-regulated (Table 3). These observations support the hypothesis that increased peroxisome proliferation results in oxidative stress, which may be due to the disproportionate increase in the level of oxidation versus antioxidation enzyme activities .
It is believed that the activation of PPARα and ensuing cascade effects are linked to both pharmacological and tumorigenic effects of PPARα agnoists . The carcinogenic response seems likely to be associated with both the induction of oxidative stress and the increased cell proliferation from peroxisome proliferation after treatment with these chemicals. In this study, we found that two genes (Acox1 and Pex11a) associated with cellular proliferation of peroxisomes were up-regulated about 3-fold. The level of peroxisomal biogenesis factor 11 (Pex11) correlates roughly with peroxisome abundance in the cell, and over-expression of Pex11 alone is sufficient to accelerate peroxisome division and to increase peroxisome abundance . It is thought that alteration in the balance between cell proliferation and apoptosis is causally related to the induction of liver tumors, and induced cell proliferation plays a key role in carcinogenesis in animals and humans [24, 25].
Based on the Ingenuity Pathway Analysis, 11 genes associated with liver cancer development were up-regulated by at least two PPARα agonists tested (Table 3). For example, Bnip3 (BCL2 19 kDa-interacting protein 1), a pro-apoptotic factors of the Bcl-2-family, has been previously shown to be up-regulated in malignant tumors . Diazepam binding inhibitor (Dbi), interacts with hepatocyte nuclear factor-4 α that transcriptionally regulates the genes involved in both lipid and glucose metabolism , was also increased. Previous studies indicated that Dbi levels are higher in hepatocellular carcinoma (HCC) patients , and the elevation of Dbi expression is useful in evaluating malignancy and in diagnostic approaches of tumors in liver tissue . Fatty acid-binding proteins (Fabps) are involved in lipid metabolism by intracellular transport of long-chain fatty acids. Liver fatty acid-binding protein (Fabp1) is demonstrated immunohistochemically in human hepatocellular malignancies, suggesting that its immunoreactivity is a candidate for the tumor marker in hepatic cell malignancies . Fatty acid synthase (Fasn) is the key enzyme of de novo fatty acid synthesis. The over-expression of Fasn is an early phenomenon presented in both hormonally and chemically induced rat hepatocarcinogenesis . Hypoxia inducible factor-1 alpha (Hif1a) regulates the expression of a myriad of genes involved in oxygen transport, glucose uptake, glycolysis and angiogenesis. The expression of Hif1a in HCC tissue is higher than that in paraneoplastic tissue or normal liver tissue, and Hif1a plays an important role in neovascularization in HCC . Lgals3 (lectin) has been demonstrated to be associated with assorted processes such as cell growth, tumor transformation and metastasis. It has been reported that Lgals3 expression was induced in cirrhotic liver and HCC, and that the expression of Lgals3 in proliferating cells possibly indicates an early neoplastic event . The plasminogen activation system, including PAI-1 (plasminogen activator inhibitor 1), plays a crucial role in the process of cancer. PAI-1 is increased in HCC, and its expression is related to the invasiveness, metastasis, and prognosis [34, 35]. Our findings support the observation that PPARα agonists increase proliferation of peroxisomes in rodent hepatocytes and alter lipid metabolism. In addition, the gene expression profiles indicate a number of leads toward understanding PPARα agonist-induced hepatocarcinogenesis in the mouse.
In summary, primary mouse hepatocytes were treated with various concentrations (10, 30, 100 μM) of three PPARα agonists (bezafibrate, fenofibrate, and WY-14,643) for 24 hr. Transactivation analysis indicated that these three agents activated the PPARα in a dose-dependent manner. Global gene expression analysis was performed on whole mouse genome arrays following exposure of the mouse hepatocyte cultures to PPARα agonists. Hierarchal clustering analysis of these gene expression profiles indicated that expression profiles of DEGs were clustered based on doses rather than specific drugs, indicating there is a common effect across this class of compound. Gene expression changes were detected in mouse hepatocytes following exposure to 100 μM bezafibrate (151 genes), fenofibrate (149 genes), and WY-14,643 (145 genes). The expression of 121 genes was changed in common by at least two of the PPARα agonists tested. Based on Ingenuity Pathway Analysis, many of these genes (49) function in lipid metabolism. An additional 11 genes were mapped to cancer associated functions. Clear dose-dependent changes in DEGs were determined based on magnitude of fold change. These results provide a better understanding of the underlying mechanisms of the hepatic effects of PPARα agonists in the mouse.
Materials and methods
Chemicals and cell treatments
Bezafibrate and fenofibrate were purchased from Sigma (St. Louis, MO). WY-14,643 was purchased from Chemsyn Science Laboratories (Lenexa, KS). All compounds were prepared as 1000 x stock solutions in dimethyl sulfoxide (DMSO) and added to cell cultures in final concentrations of 10–100 μM. The same amount of DMSO (0.1% v/v) was added to control cells. During the treatments, serum-free medium was used with supplements (see below).
Three 6–8 week-old C57/BL6 male mice were obtained from the breeding colony of the FDA's National Center for Toxicological Research. Mice were anesthetized with 1.5 ml/kg of nembutal sodium solution containing 50 mg/ml of pentobarbital sodium prior to undergoing liver perfusion. All animals used in this study were handled in accordance with the principles and guidelines prepared by the National Institutes of Health, USA. Mouse primary hepatocytes were isolated by a two-stage collagenase perfusion process according to the methods described by Seglen et al. and Kreamer et al. [36, 37]. Primary hepatocytes were suspended in L-15 medium containing 2 mg/ml BSA, 18 mM HEPES, 3 mg/ml proline, 1 mg/ml galactose, 0.1% insulin-transferrin-selenite, 10 ng/ml epidermal growth factor, 50 U/ml penicillin, and 50 μg/ml streptomycin. The cells were treated with drugs six hours after plating.
Plasmid transfection and reporter assays
Luciferase reporter plasmid, pHD(x3) luc, contains three direct tandem copies of PPRE binding site, was used as previously described . pSG5-PPARα and pSG5-RXR expression plasmids were obtained from Dr. Marek Michalak (University of Alberta, Canada). Plasmid DNAs were purified by column chromatography (Qiagen Inc., Valencia, CA). Primary hepatocytes were plated in L15 medium at 1 × 105 cells/per well for 6 wells, and transfections were carried out after cell attachment with FuGENE reagent (Roche Diagnostics, Indianapolis, IN). Briefly, 300 μl of L15 (no additives) containing 9 μl of FuGENE reagent was mixed with a total of 3 μg of plasmid DNA with or without pSG5-PPARα/pSG5-RXR, luciferase reporter plasmid, and pSVβ-gal (internal control). This mixture was added to cells for a 10 min incubation at room temperature. For induction, medium was replaced with fresh L15 without BSA after 12 h incubation with plasmids/FuGENE. The PPARα agonists were added for 24 h at appropriate concentrations. To harvest lysates for luciferase activity, hepatocytes were washed twice in PBS, and then lysed in 150 μl of 1x reporter lysis buffer (Promega Corp., Madison, WI). Luciferase activity was measured using the Luciferase Assay System (Promega). Luminescence was determined using an automatic luminometer, LumiTeum II (Harta Instruments, Gaithersburg, MD). β-galactosidase enzyme assay was carried out using β-galactosidase Enzyme Assay System (Promega). Luciferase activity was normalized against β-galactosidase activity from the same lysate. Each assay was performed in duplicates.
RNA isolation and quality control
Total RNA from cells was isolated using an RNeasy system (Qiagen). The yield of the extracted RNA was determined spectrophotometrically by measuring the optical density at 260 nm. The purity and quality of extracted RNA were evaluated using the RNA 6000 LabChip and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA). Only high quality RNA with RNA integrity numbers (RINs) greater than 7.5 were used for microarray experiments.
Preparation of labeled in vitro transcribed cRNA targets
All total RNA samples were labeled by direct incorporation of cyanine 3 or cyanine 5 dyes using the Agilent Low RNA Input Linear Amplification Kit (Santa Clara, CA). A 500 ng quantity of total RNA was input to each reaction. Labeled cRNAs were purified using the Qiagen RNeasy Mini kit, and were analyzed for quality and quantity using standard UV spectrometry and the Agilent Bioanalyzer.
Hybridization of labeled cRNA to microarrays and microarray imaging
Cyanine 3 labeled cRNAs were mixed with cyanine 5 labeled cRNAs for hybridization to microarrays. Each test sample was hybridized against its corresponding control sample to two microarrays in a dye-swapped pair. Agilent Whole Mouse Genome Microarrays (Santa Clara, CA) were hybridized using the Agilent Gene Expression Hybridization kit and washed using Agilent Gene Expression Wash Buffers according to the manufacturers' protocols. Hybridized microarrays were scanned using the Agilent DNA Microarray Scanner and data were extracted from images using the Agilent Feature Extraction (version 7.5) software using default settings.
Microarray data analysis
The Agilent Whole Mouse Genome Microarray that contains 43,790 probes was used to generate gene expression profiles for bezafibrate, fenofibrate, and WY-14,643 at three dose levels (10 μM, 30 μM and 100 μM) with three biological replicates (hepatocytes isolated from mouse A, B, C). Each treated sample was paired with a control (non-treated hepatocytes) using a dye swap experiment design, resulting in 54 arrays [3 chemicals × 3 doses × 3 animals × (2 dye swap)]. In addition, self-self hybridizations were also conducted for each of three controls with three technical replicates, resulting in nine additional arrays. In total, 63 hybridizations were performed for this study.
Linear & Lowess method consists of median scaling to 1000 for each channel per array with a follow up Lowess normalization. The parameters used in Lowess normalization were: smoothing factor = 0.2 and robustness iterations = 3. Low intensity (<500) spots were filtered out after normalization and a subset of 25,010 genes was generated for further data analysis. The DEGs were identified using a combination of Student t-test and fold change (FC). A gene was considered differentially expressed if P-value was less than 0.05 and the FC was greater than 1.5. Two lists of DEGs were obtained. One list was obtained by comparing the nine self-self hybridization arrays with those polarity+ arrays in which control and treatment samples were labeled with Cy3 and Cy5, respectively. The other list was obtained by comparing the nine self-self hybridization arrays with those polarity- arrays in which control and treatment samples were labeled with Cy5 and Cy3, respectively. The genes in common between these two lists of DEGs were considered as final DEGs and used for biological interpretation. The described analysis approach, in terms of self-self hybridization experiment design, will be reported somewhere else, where its comparative performance analysis was conducted.
Microarray data management and analysis were conducted using an FDA microarray software, ArrayTrack [39, 40]. ArrayTrack also provides functionality for the interpretation of gene expression data. For example, pathway analysis is based on the Pathway Library in ArrayTrack, which contains pathways from Kyoto Encyclopedia of Genes and Genomes (KEGG)  and PathArt (Jubilant Biosys Ltd., Columbia, MD). The Fisher Exact Test  is implemented in ArrayTrack to assess the statistical significance of identified pathways. Ingenuity Pathway Analysis software (Mountain View, CA) was also used for gene function and pathway analysis. S-Plus (Insightful Corp., Seattle, WA) was used in this study for the statistical calculation.
Schoonjans K, Staels B, Auwerx J: Role of the peroxisome proliferator-activated receptor (PPAR) in mediating the effects of fibrates and fatty acids on gene expression. J Lipid Res 1996, 37(5):907–925.
Schoonjans K, Staels B, Auwerx J: The peroxisome proliferator activated receptors (PPARS) and their effects on lipid metabolism and adipocyte differentiation. Biochim Biophys Acta 1996, 1302(2):93–109.
Corton JC, Anderson SP, Stauber A: Central role of peroxisome proliferator-activated receptors in the actions of peroxisome proliferators. Annu Rev Pharmacol Toxicol 2000, 40: 491–518. 10.1146/annurev.pharmtox.40.1.491
Tugwood JD, Issemann I, Anderson RG, Bundell KR, McPheat WL, Green S: The mouse peroxisome proliferator activated receptor recognizes a response element in the 5' flanking sequence of the rat acyl CoA oxidase gene. Embo J 1992, 11(2):433–439.
Dreyer C, Krey G, Keller H, Givel F, Helftenbein G, Wahli W: Control of the peroxisomal beta-oxidation pathway by a novel family of nuclear hormone receptors. Cell 1992, 68(5):879–887. 10.1016/0092-8674(92)90031-7
Muerhoff AS, Griffin KJ, Johnson EF: The peroxisome proliferator-activated receptor mediates the induction of CYP4A6, a cytochrome P450 fatty acid omega-hydroxylase, by clofibric acid. J Biol Chem 1992, 267(27):19051–19053.
Vanden Heuvel JP, Sterchele PF, Nesbit DJ, Peterson RE: Coordinate induction of acyl-CoA binding protein, fatty acid binding protein and peroxisomal beta-oxidation by peroxisome proliferators. Biochim Biophys Acta 1993, 1177(2):183–190. 10.1016/0167-4889(93)90039-R
Klaunig JE, Babich MA, Baetcke KP, Cook JC, Corton JC, David RM, DeLuca JG, Lai DY, McKee RH, Peters JM, et al.: PPARalpha agonist-induced rodent tumors: modes of action and human relevance. Crit Rev Toxicol 2003, 33(6):655–780.
Yeldandi AV, Rao MS, Reddy JK: Hydrogen peroxide generation in peroxisome proliferator-induced oncogenesis. Mutat Res 2000, 448(2):159–177.
Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 1995, 270(5235):467–470.
Lipshutz RJ, Fodor SP, Gingeras TR, Lockhart DJ: High density synthetic oligonucleotide arrays. Nat Genet 1999, 21(1 Suppl):20–24. 10.1038/4447
Burczynski ME, McMillian M, Ciervo J, Li L, Parker JB, Dunn RT 2nd, Hicken S, Farr S, Johnson MD: Toxicogenomics-based discrimination of toxic mechanism in HepG2 human hepatoma cells. Toxicol Sci 2000, 58(2):399–415. 10.1093/toxsci/58.2.399
Waring JF, Ciurlionis R, Jolly RA, Heindel M, Ulrich RG: Microarray analysis of hepatotoxins in vitro reveals a correlation between gene expression profiles and mechanisms of toxicity. Toxicol Lett 2001, 120(1–3):359–368. 10.1016/S0378-4274(01)00267-3
Issemann I, Green S: Activation of a member of the steroid hormone receptor superfamily by peroxisome proliferators. Nature 1990, 347(6294):645–650. 10.1038/347645a0
Lee SS, Pineau T, Drago J, Lee EJ, Owens JW, Kroetz DL, Fernandez-Salguero PM, Westphal H, Gonzalez FJ: Targeted disruption of the alpha isoform of the peroxisome proliferator-activated receptor gene in mice results in abolishment of the pleiotropic effects of peroxisome proliferators. Mol Cell Biol 1995, 15(6):3012–3022.
Peters JM, Aoyama T, Cattley RC, Nobumitsu U, Hashimoto T, Gonzalez FJ: Role of peroxisome proliferator-activated receptor alpha in altered cell cycle regulation in mouse liver. Carcinogenesis 1998, 19(11):1989–1994. 10.1093/carcin/19.11.1989
Vanden Heuvel JP, Kreder D, Belda B, Hannon DB, Nugent CA, Burns KA, Taylor MJ: Comprehensive analysis of gene expression in rat and human hepatoma cells exposed to the peroxisome proliferator WY14,643. Toxicol Appl Pharmacol 2003, 188(3):185–198. 10.1016/S0041-008X(03)00015-2
Bardot O, Aldridge TC, Latruffe N, Green S: PPAR-RXR heterodimer activates a peroxisome proliferator response element upstream of the bifunctional enzyme gene. Biochem Biophys Res Commun 1993, 192(1):37–45. 10.1006/bbrc.1993.1378
Schoonjans K, Watanabe M, Suzuki H, Mahfoudi A, Krey G, Wahli W, Grimaldi P, Staels B, Yamamoto T, Auwerx J: Induction of the acyl-coenzyme A synthetase gene by fibrates and fatty acids is mediated by a peroxisome proliferator response element in the C promoter. J Biol Chem 1995, 270(33):19269–19276. 10.1074/jbc.270.33.19269
Cornwell PD, De Souza AT, Ulrich RG: Profiling of hepatic gene expression in rats treated with fibric acid analogs. Mutat Res 2004, 549(1–2):131–145.
Klaunig JE, Xu Y, Isenberg JS, Bachowski S, Kolaja KL, Jiang J, Stevenson DE, Walborg EF Jr: The role of oxidative stress in chemical carcinogenesis. Environ Health Perspect 1998, 106(Suppl 1):289–295.
Cho S, Park EM, Febbraio M, Anrather J, Park L, Racchumi G, Silverstein RL, Iadecola C: The class B scavenger receptor CD36 mediates free radical production and tissue injury in cerebral ischemia. J Neurosci 2005, 25(10):2504–2512. 10.1523/JNEUROSCI.0035-05.2005
Li X, Gould SJ: The dynamin-like GTPase DLP1 is essential for peroxisome division and is recruited to peroxisomes in part by PEX11. J Biol Chem 2003, 278(19):17012–17020. 10.1074/jbc.M212031200
Khan SH, Sorof S: Liver fatty acid-binding protein: specific mediator of the mitogenesis induced by two classes of carcinogenic peroxisome proliferators. Proc Natl Acad Sci USA 1994, 91(3):848–852. 10.1073/pnas.91.3.848
Roberts RA, James NH, Hasmall SC, Holden PR, Lambe K, Macdonald N, West D, Woodyatt NJ, Whitcome D: Apoptosis and proliferation in nongenotoxic carcinogenesis: species differences and role of PPARalpha. Toxicol Lett 2000, 112–113: 49–57. 10.1016/S0378-4274(99)00243-X
Stepan H, Leo C, Purz S, Hockel M, Horn LC: Placental localization and expression of the cell death factors BNip3 and Nix in preeclampsia, intrauterine growth retardation and HELLP syndrome. Eur J Obstet Gynecol Reprod Biol 2005, 122(2):172–176. 10.1016/j.ejogrb.2005.01.017
Petrescu AD, Payne HR, Boedecker A, Chao H, Hertz R, Bar-Tana J, Schroeder F, Kier AB: Physical and functional interaction of Acyl-CoA-binding protein with hepatocyte nuclear factor-4 alpha. J Biol Chem 2003, 278(51):51813–51824. 10.1074/jbc.M303858200
Venturini I, Zeneroli ML, Corsi L, Baraldi C, Ferrarese C, Pecora N, Frigo M, Alho H, Farina F, Baraldi M: Diazepam binding inhibitor and total cholesterol plasma levels in cirrhosis and hepatocellular carcinoma. Regul Pept 1998, 74(1):31–34. 10.1016/S0167-0115(98)00013-5
Venturini I, Alho H, Podkletnova I, Corsi L, Rybnikova E, Pellicci R, Baraldi M, Pelto-Huikko M, Helen P, Zeneroli ML: Increased expression of peripheral benzodiazepine receptors and diazepam binding inhibitor in human tumors sited in the liver. Life Sci 1999, 65(21):2223–2231. 10.1016/S0024-3205(99)00487-7
Suzuki T, Watanabe K, Ono T: Immunohistochemical demonstration of liver fatty acid-binding protein in human hepatocellular malignancies. J Pathol 1990, 161(1):79–83. 10.1002/path.1711610113
Evert M, Schneider-Stock R, Dombrowski F: Overexpression of fatty acid synthase in chemically and hormonally induced hepatocarcinogenesis of the rat. Lab Invest 2005, 85(1):99–108.
Huang GW, Yang LY, Lu WQ: Expression of hypoxia-inducible factor 1alpha and vascular endothelial growth factor in hepatocellular carcinoma: Impact on neovascularization and survival. World J Gastroenterol 2005, 11(11):1705–1708.
Hsu DK, Dowling CA, Jeng KC, Chen JT, Yang RY, Liu FT: Galectin-3 expression is induced in cirrhotic liver and hepatocellular carcinoma. Int J Cancer 1999, 81(4):519–526. 10.1002/(SICI)1097-0215(19990517)81:4<519::AID-IJC3>3.0.CO;2-0
Itoh T, Hayashi Y, Kanamaru T, Morita Y, Suzuki S, Wang W, Zhou L, Rui JA, Yamamoto M, Kuroda Y, et al.: Clinical significance of urokinase-type plasminogen activator activity in hepatocellular carcinoma. J Gastroenterol Hepatol 2000, 15(4):422–430. 10.1046/j.1440-1746.2000.02150.x
Zheng Q, Tang ZY, Xue Q, Shi DR, Song HY, Tang HB: Invasion and metastasis of hepatocellular carcinoma in relation to urokinase-type plasminogen activator, its receptor and inhibitor. J Cancer Res Clin Oncol 2000, 126(11):641–646. 10.1007/s004320000146
Kreamer BL, Staecker JL, Sawada N, Sattler GL, Hsia MT, Pitot HC: Use of a low-speed, iso-density percoll centrifugation method to increase the viability of isolated rat hepatocyte preparations. In Vitro Cell Dev Biol 1986, 22(4):201–211.
Seglen PO: Preparation of isolated rat liver cells. Methods Cell Biol 1976, 13: 29–83.
Zhang B, Marcus SL, Miyata KS, Subramani S, Capone JP, Rachubinski RA: Characterization of protein-DNA interactions within the peroxisome proliferator-responsive element of the rat hydratase-dehydrogenase gene. J Biol Chem 1993, 268(17):12939–12945.
Tong W, Harris S, Cao X, Fang H, Shi L, Sun H, Fuscoe J, Harris A, Hong H, Xie Q, et al.: Development of public toxicogenomics software for microarray data management and analysis. Mutat Res 2004, 549(1–2):241–253.
Tong W, Cao X, Harris S, Sun H, Fang H, Fuscoe J, Harris A, Hong H, Xie Q, Perkins R, et al.: ArrayTrack--supporting toxicogenomic research at the U.S. Food and Drug Administration National Center for Toxicological Research. Environ Health Perspect 2003, 111(15):1819–1826.
Kanehisa M, Goto S, Kawashima S, Nakaya A: The KEGG databases at GenomeNet. Nucleic Acids Res 2002, 30(1):42–46. 10.1093/nar/30.1.42
Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S, et al.: GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol 2003, 4(4):R28. 10.1186/gb-2003-4-4-r28
Helpful discussions, comments offered by Dr James Fuscoe from NCTR are greatly appreciated.
The views presented in this article do not necessarily reflect those of the US Food and Drug Administration.
LG performed the analysis of microarray data and wrote the manuscript. HF, XHF and WT performed data analysis. LS and YD helped writing the manuscript. JC, AW, KM conducted the microarray experiment and generated the raw data. SD and EB performed hepatocytes isolation and RNA extraction.
About this article
Cite this article
Guo, L., Fang, H., Collins, J. et al. Differential gene expression in mouse primary hepatocytes exposed to the peroxisome proliferator-activated receptor α agonists. BMC Bioinformatics 7 (Suppl 2), S18 (2006). https://doi.org/10.1186/1471-2105-7-S2-S18
- Ingenuity Pathway Analysis
- Liver Fatty Acid Binding Protein
- Mouse Primary Hepatocyte