Analysis of the heat shock response in mouse liver reveals transcriptional dependence on the nuclear receptor peroxisome proliferator-activated receptor α (PPARα)
© Vallanat et al; licensee BioMed Central Ltd. 2010
Received: 17 February 2009
Accepted: 7 January 2010
Published: 7 January 2010
The nuclear receptor peroxisome proliferator-activated receptor alpha (PPARα) regulates responses to chemical or physical stress in part by altering expression of genes involved in proteome maintenance. Many of these genes are also transcriptionally regulated by heat shock (HS) through activation by HS factor-1 (HSF1). We hypothesized that there are interactions on a genetic level between PPARα and the HS response mediated by HSF1.
Wild-type and PPARα-null mice were exposed to HS, the PPARα agonist WY-14,643 (WY), or both; gene and protein expression was examined in the livers of the mice 4 or 24 hrs after HS. Gene expression profiling identified a number of Hsp family members that were altered similarly in both mouse strains. However, most of the targets of HS did not overlap between strains. A subset of genes was shown by microarray and RT-PCR to be regulated by HS in a PPARα-dependent manner. HS also down-regulated a large set of mitochondrial genes specifically in PPARα-null mice that are known targets of PPARγ co-activator-1 (PGC-1) family members. Pretreatment of PPARα-null mice with WY increased expression of PGC-1β and target genes and prevented the down-regulation of the mitochondrial genes by HS. A comparison of HS genes regulated in our dataset with those identified in wild-type and HSF1-null mouse embryonic fibroblasts indicated that although many HS genes are regulated independently of both PPARα and HSF1, a number require both factors for HS responsiveness.
These findings demonstrate that the PPARα genotype has a dramatic effect on the transcriptional targets of HS and support an expanded role for PPARα in the regulation of proteome maintenance genes after exposure to diverse forms of environmental stress including HS.
Physiological and chemical stresses produce disease states in which proteins are damaged or misfolded in part through increases in oxidative stress. Many endogenous pathways are engaged in restoring cellular homeostasis, including stabilization of unfolded proteins to prevent aggregation and removing damaged or excess proteins through proteolysis. Stabilization of unfolded proteins is performed by molecular chaperones that assist in the folding of nascent polypeptides. Many genes encoding for chaperones exhibit increased expression after exposure to a wide variety of stimuli including chemical exposure or increased temperatures and are thus called heat shock (HS) protein (Hsp) genes [1–3]. These proteins play key roles in a number of human diseases . Expression of some Hsp is essential for cellular survival under physical or chemical exposure conditions that increase oxidative stress [5, 6]. Regulation of the Hsp genes by heat or chemical-induced oxidative stress is controlled in part by HS factor 1 (HSF1), activated under conditions in which the level of unfolded proteins increase [1–3]. Microarray studies of mouse embryonic fibroblasts from wild-type and HSF1-null mice or a human cervical carcinoma cell line have shown that HSF1 controls only a subset of the genes altered by HS [7, 8] indicating that other inducible pathways play roles in regulating the Hsp genes.
The nuclear receptor peroxisome proliferator-activated receptor α (PPARα) is one of three PPAR subtypes that regulate lipid and glucose homeostasis, tissue growth and inflammation after exposure to a large class of structurally heterogeneous industrial and pharmaceutical chemicals called peroxisome proliferator chemicals (PPC) [9, 10]. The PPARα subtype plays a key role in mediating the effects of hypolipidemic and xenobiotic PPC in liver, kidney, heart and skin. Exposure to PPC leads to regulation of a large number of genes including up-regulation of those involved in lipid homeostasis and down-regulation of inflammatory genes that are generally abolished in PPC-exposed PPARα-null mice [11–13].
There is compelling evidence that PPARα protects tissues from chemical-induced oxidative stress (reviewed in ). Prior exposure of rats and mice to PPC protects the liver from damage by cytotoxic agents that induce oxidative stress . The hypolipidemic drug and PPC, clofibrate, protects the liver from damage by the cytotoxicant acetaminophen in wild-type but not PPARα-null mice . Compared to wild-type mice, untreated PPARα-null mice or primary hepatocytes isolated from PPARα-null mice were more sensitive to carbon tetrachloride-, paraquat- or cadmium-induced toxicity . The beneficial effects of caloric restriction in protecting the liver from cytotoxicant-induced liver injury were shown to depend on PPARα . Additional studies have shown that PPARα plays a positive role in recovering from partial hepatectomy . In the kidney, PPARα-null mice were more sensitive to damage after ischemia-reperfusion injury, and prior exposure of wild-type mice to PPC reduces the injury [19, 20]. Our previous microarray studies identified an overlap in the genes regulated by the PPARα agonist WY-14,643 (WY) and those regulated by HS through HSF1 . The proteome maintenance genes included those involved in protein folding (e.g., Hsp genes) as well as ubiquitin-dependent and -independent proteolytic processing through the proteosome (e.g., Psm genes). Altered regulation of these genes by PPC could help to explain why PPC exposure through PPARα helps to protect tissues from environmental stressors.
In the present study, we hypothesized that there were interactions on a genetic level between PPARα and the HS response mediated by HSF1. We dissected the contribution of PPARα by examining gene expression changes in the livers of wild-type and PPARα-null mice after HS. We found that a number of genes regulated by HS were dependent on PPARα for regulation. Furthermore, we showed that PPARα-null mice were more sensitive to the transcriptional effects of HS and exhibited a remarkably different transcriptional response compared to wild-type mice. Our findings suggest PPARα is a major regulator of stress responses in the liver.
Animals, exposure to heat and WY
All animal studies were carried out at CIIT Centers for Health Research, Research Triangle Park, NC. Studies utilized wild-type and PPARα-null male mice 9-12 weeks of age on a mixed SV129/C57BL/6 background. These mice have been described previously . The mice were originally obtained from Dr. Frank Gonzalez, NCI, National Institutes of Health to establish a breeding colony at CIIT. Control and treated mice were provided with NIH-07 rodent chow (Zeigler Brothers, Gardeners, PA) and deionized, filtered water ad libitum. Lighting was on a 12-hr light/dark cycle. Mice (n = 5-6 per group) were fed a control diet or a diet containing the PPARα agonist WY-14,643 (ChemSyn Science, Lenexa, KS) (WY) in the diet (500 ppm) for 1 week followed by a 40-min heat stress at 42°C or held at room temperature in wire rack cages. Mice were sacrificed 4 hrs, 24 hrs and 3 days after the heat stress.
Heat in the exposure chamber was generated by a combination of the building reheat system, an in-line air heater and two small surface heaters. The in-line heater (Riley Equipment Company, Inc., Houston, TX) was set to approximately 300°F. The heated air traveled through a 2 foot by 2 inch inlet tube constructed of PVC, stainless steel and Teflon into the bottom of the H1000 inhalation chamber. A small diverter was placed in the bottom of the H1000 inhalation chamber to prevent a chimney effect in the chamber and increasing the efficiency of the heat distribution. The stainless steel animal cage unit was placed in the center of the H1000 inhalation chamber and a feces/urine catch pan was placed in the chamber. A small (1500 watt) heater was placed on either side of the H1000 inhalation chamber approximately 4 inches from the surface. These heaters heated the stainless steel sides and heated the surface enough to allow the internal temperature of the chamber to reach approximately 42°C. The temperature was monitored in three locations of the H1000 inhalation chamber: front door above the animal cage unit, back door below the animal cage unit and in the center using temperature probes (Rotronic Hygrometer Series 1200, Rotronic Instrument Corp., Huntington, NY).
Portions of the livers were rapidly snap-frozen in liquid nitrogen and stored at -70°C until analysis. Where appropriate, slices of liver were fixed in 10% neutral buffered formalin for 48 h, transferred to 70% ethanol, and embedded in paraffin; 5-μm sections were cut and mounted on slides and stained with H&E. H&E-stained liver sections were examined by light microscopy. Clinical chemistry markers in the serum from animals were measured by Antech Diagnostics, Charlotte, NC. All animal studies were conducted under federal guidelines for the use and care of laboratory animals and were approved by Institutional Animal Care and Use Committees.
RNA Isolation and Analysis of Gene Expression
Liver RNA was isolated using a modified guanidium isothiocyanate method (TRIzol®; Invitrogen) and was further purified using silica membrane spin columns (RNeasy®; Qiagen, Valencia, CA). RNA integrity was assessed by the RNA 6000 LabChip® kit using a 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). Global gene expression changes were examined using the Affymetrix platform. Gene expression changes were assessed in the livers from three mice in each of the 8 treatment groups in wild-type and PPARα-null mice using a total of 24 chips. Gene expression in each animal was assayed on a separate chip. Biotin-labeled cRNA was produced from 15 μg total RNA using an Affymetrix "one-way" labeling kit. Total cRNA was quantified using a Nano-Drop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE) and evaluated for quality after fragmentation on a 2100 Bioanalyzer. Following overnight hybridization at 45°C to Affymetrix U74Av2 GeneChips in an Affymetrix Model 640 GeneChip hybridization oven, the arrays were washed and stained using an Affymetrix 450 fluidics station as recommended by the manufacturer and scanned on an Affymetrix Model 3000 scanner. After scanning, raw data (Affymetrix .cel files) were obtained using Affymetrix GeneChip Operating Software (version 1.4). All of the Affymetrix (Santa Clara, CA) .cel files were analyzed by Bioconductor SimpleAffy to assess data quality .
Data (.cel files) was analyzed and statistically filtered using Rosetta Resolver® version 7.1 software (Rosetta Inpharmatics, Kirkland, WA). Statistically significant genes were identified using one-way ANOVA with a false discovery rate (Benjamini-Hochberg test) of ≤ 0.05 followed by a post-hoc test (Scheffe) for significance. Significant transcripts were evaluated for relevance to canonical pathways and biological functions using Ingenuity Pathways Analysis (Ingenuity Systems, http://www.ingenuity.com). Heat maps were generated using Eisen Lab Cluster and Treeview software http://rana.lbl.gov/EisenSoftware.htm.
Gene Set Enrichment Analysis (GSEA; http://www.broad.mit.edu/gsea/) was used to evaluate whether a pre-defined set of genes showed statistically significant, concordant differences between two biological states . GSEA generates a list of genes from submitted .cel files and tests whether, within that queried list of genes, there is a statistically significant enrichment or not of pre-defined groups of genes, or "gene sets". Gene sets examined through GSEA include canonical, metabolic and signaling pathways as well as groups of genes previously identified and validated to be up- or down-regulated when cells are given a particular stimulus. This type of analysis can sometimes detect more subtle changes present in the data. Expression profiles were submitted to GSEA using default settings and searched for enriched gene sets. Gene set C2 (curated gene sets) includes genes from online pathway databases, publications in PubMed, and knowledge from domain experts. Gene set C3 (motif gene sets) contains genes that share a transcription factor (TF) binding site defined in TRANSFAC®, a database for transcription factors and their genomic binding sites. Additional .cel files were analyzed using GSEA including those from the livers of wild-type and PPARα-null mice exposed to WY (400 uL of a 10 mg/mL solution/day) for 6 hr or 5 days .
A description of the HS and WY exposure microarray experiment is available through Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov/geo/, as accession number GSE14,869.
We categorized the genes regulated by HS in wild-type or PPARα-null mice based on whether they were altered by HS in mouse embryonic fibroblasts (MEF) from wild-type or HSF1-null mice . From this dataset, genes which fell into three groups were identified including 1) HSF1-dependent regulation by HS, 2) HSF1-independent regulation by HS and 3) altered regulation between wild-type and HSF1-null MEFs in the absence of HS. Genes were identified using procedures outlined in .
Real Time RT-PCR Analysis
The levels of gene expression were quantified using real time RT-PCR analysis. Real-time reverse transcriptase-PCR was performed as follows. After DNase treatment, total RNA was quantified (Ribogreen®, Molecular Probes, Inc., Eugene, OR) and diluted with water. Fifty ng of RNA and PCR reagents were aliquoted into 96-well plates using an ABI Prism™ 6700 Automated Nucleic Acid Workstation (Applied Biosystems, Foster City, CA) and subjected to real-time quantitative PCR (TaqMan®, Applied Biosystems) using gene-specific primers (Additional File 1) and fluorescently labeled probes (Molecular Probes) designed by the Primer Express® software (Applied Biosystems). Amplification curves were generated using the ABI Prism™ 7900 HT Sequence Detection System (Applied Biosystems). Expression relative to vehicle control animals was determined after normalizing to the ribosomal 18S gene. There were four animals per treatment group, and each sample was analyzed in duplicate. Variability is expressed as standard error of the mean. Means and S.E. (n = 4) for RT-PCR data were calculated by Student's t test. The level of significance was set at p ≤ 0.05.
Western Blot Analysis
Liver lysates were prepared in 250 mM sucrose, 10 mM Tris-HCl, pH 7.4, and 1 mM EDTA with protease inhibitors (0.2 mM phenylmethylsulphonyl fluoride, 0.1% aprotinin, 1 μg/ml pepstatin, 1 μg/ml leupeptin) as previously described . Fifty μg hepatocyte whole cell lysate was subjected to 12% sodium dodecyl sulfate - polyacrylamide gel electrophoresis followed by transfer to nitrocellulose membranes. Immunoblots were developed using primary antibodies against acyl-CoA oxidase (ACO) (a gift from S. Alexson, Huddinge University Hospital, Huddinge, Sweden), HS proteins (Santa Cruz Biotechnology, Santa Cruz, CA; StressGen, Victoria, B.C., Canada) or CYP4A (GenTest, Waltham, MA) and appropriate secondary antibodies conjugated with horseradish peroxidase (Santa Cruz Biotechnology) in the presence of chemiluminescent substrate ECL (Amersham, Piscataway, NJ). Blots were quantitated using Gel-Pro (MediaCybernetics, Silver Spring, MD). Most antibodies recognized only one major band of the expected size. Antibodies to TCP1η routinely gave 2 bands of ~60 kDa and ~40 kDa, both of which were elevated after HS, WY or both treatments. In this study we report the levels of ACO-B protein, the 52 kDa processed form of the full-length ACO protein (ACO-A) .
Results and Discussion
Identification of PPARα-dependent heat shock-regulated genes
The gene expression pattern after treatment was dramatically different in PPARα-null mice. HS, WY alone or WY+HS treatments in PPARα-null mice altered 774, 6, or 68 genes, respectively (Figure 1B).
The HS genes were categorized based on their expression behavior in the two strains after HS and WY exposure (Figure 2B). Genes in class I and II were altered by either HS alone (class I) or by HS and WY (class II) in a PPARα-dependent manner. The class I genes included FK506 binding protein 4 (Fkbp4) and GrpE-like 2, mitochondrial gene (Grpel2). Also included was ubiquitin carboxy-terminal hydrolase L1 (Uchl1) identified using less stringent microarray analysis cut-offs but confirmed by RT-PCR (Figure 2C). The class II genes included Dnaja1 (Hsp40) and a tubulin family member (Tubb2a). RT-PCR revealed concordance with array expression for Tubb2a; Dnaja1 exhibited weak induction by HS in PPARα-null mice. Genes in class III were altered only by HS and were partially, if not completely independent of PPARα, including stress-inducible phosphoprotein-1 (Stip1) in which increases in expression after HS in PPARα-null mice were detected by RT-PCR. Genes in class IV were induced by HS independent of PPARα and also exhibited PPARα-dependent WY induction. These genes included a number of chaperones Hsp90aa1, Hspe1, Hspa4l, Hspd1 and Hsph1. The expression pattern of Hsp90aa1 and Hsph1 was confirmed by RT-PCR. Class V genes were those that exhibited altered expression only in PPARα-null mice after HS but not WY and included two metallothionein genes (Mt1, Mt2). Class VI genes were those that exhibited altered regulation after HS in PPARα-null mice but not after HS+WY indicating suppression of HS transcriptional effects by WY. These genes included the up-regulation of serum amyloid protein family members (Saa2, Saa3) and down-regulation of a large number of mitochondrial genes (described below). Overall, the data is consistent with a subset of 1) HS-regulated genes (Class I and II) that are dependent on PPARα for altered expression, 2) HS-regulated genes (Class III and IV) that are independent of PPARα and 3) genes responsive to HS only in the absence of PPARα (class V and VI).
Heat shock activates different biological responses in wild-type and PPARα-null mice
We identified sets of genes that exhibited altered regulation in the two strains using Gene Set Enrichment Analysis (GSEA). GSEA is an analysis method that evaluates the expression of biological pathways on defined gene sets, as an alternative to examining individual genes, to assist in identifying significant biological changes in microarray data sets . Consistent with HS effects on proteome maintenance genes, HS in wild-type mice increased the expression of genes involved in protein synthesis (HSA03010_RIBOSOME, RIBOSOMAL_PROTEINS, TRNA_SYNTHETASES) and degradation (HSA03050_PROTEASOME, PROTEASOME_DEGRADATION, PROTEASOMEPATHWAY, PROTEASOME) (Additional File 3). There were no gene sets that exhibited significant overlap with the down-regulated HS wild-type genes.
HS in PPARα-null mice altered the expression of genes in a different set of biological categories (Additional File 4) including the up-regulation of genes associated with mouse liver tumors (LEE_DENA_UP, LEE_MYC_E2F1_UP, LEE_E2F1_UP, LEE_MYC_TGFA_UP, LEE_ACOX1_UP, LEE_CIP_UP, LEE_MYC_TGFA_DN), muscle contraction (STRIATED_MUSCLE_CONTRACTION) and the aging kidney (AGEING_KIDNEY_UP).
There were a large number of genes down-regulated by HS in PPARα-null mice associated with mitochondrial biogenesis and function regulated by the PPARγ Co-activator 1 (PGC-1) family which includes PGC-1α, PGC-1β and PRC (MOOTHA_VOXPHOS, PGC, HUMAN_MITODB_6_2002, MITOCHONDRIA) (Additional File 5). There was also down-regulation of genes involved in mouse and human liver tumors (LEE_MYC_E2F1_DN, LEE_DENA_DN, LEE_MYC_TGFA_DN, LEE_CIP_DN, HCC_SURVIVAL_GOOD_VS_POOR_UP) and down-regulation of genes associated with amino acid and fatty acid metabolism (VALINE_LEUCINE_AND_ISOLEUCINE_DEGRADATION, HSA00280_VALINE_LEUCINE_AND_ISOLEUCINE_DEGRADATION, HSA00380_TRYPTOPHAN_METABOLISM, HSA00071_FATTY_ACID_METABOLISM). The fatty acid metabolism genes included those that are regulated by PPARα upon exposure to PPC (discussed below).
GSEA was also used to identify candidate transcription factors responsible for altered regulation of the HS genes. A number of stress-inducible transcription factors are likely involved in the altered expression of HS genes in wild-type mice. Consistent with well-characterized mechanisms of HS gene regulation, GSEA twice identified HSF1 as a candidate transcription factor of the HS genes (Additional File 6). Two sets of genes positively regulated by HS overlapped with those regulated by X-box binding protein (XBP1). XBP1 is a key component of the endoplasmic reticulum (ER) stress response . Accumulation of unfolded proteins in the ER activates the unfolded protein response (UPR), resulting in transcriptional induction of ER chaperones and proteases in part through activation of XBP1. A number of genes associated with XBP1 were either up- (HOXC4, STAT3, PAK1IP1, SGK, VDP, SEC61A1, SFRS2, UBQLN1, MORF4L2) or down-regulated (COL1A2, ATBF1, POLG, GRIA3) in wild-type mice. HS genes linked to up-regulation of nuclear respiratory factor 1 (NRF1) were also identified (data not shown). NRF1 is activated by increased levels of heme and regulates a large number of genes involved in mitochondrial biogenesis and oxidative stress that overlap with those regulated by PGC-1α or PGC-1β[30, 31].
In contrast, the genes regulated by HS in PPARα-null mice overlapped with serum response factor (SRF) (one set significantly and three approaching significance) (Additional File 7). SRF is a ubiquitously expressed transcription factor that plays a role in regulating early response genes to various stimuli including healing after tissue injury . SRF plays a prominent role in regulating the expression of heart muscle genes  and may explain increases in muscle contraction genes (ACTA1, ACTA2, DES, MYL1, MYL2, MYL3, MYL9, MYOM1, TNNC2, TNNT1, TNNT3, TPM2, TPM4) after HS. Two gene sets that approached significance were regulated by HSF1. In summary, these results indicate that the PPARα genotype has a dramatic effect on the biological pathways that are affected by HS.
PPARα-independent changes in triglyceride levels and fatty acid oxidation protein expression by heat shock
These findings indicate that HS elicits complex changes in fatty acid metabolism gene and protein expression leading to alterations in triglyceride levels, some of which are independent of PPARα. The increases in triglyceride levels in wild-type mice by HS at 4 hrs may be attributed in part to decreases in fatty acid oxidation genes; the decreases in triglyceride levels in PPARα-null mice are correlated with increases in Cyp4a and ACO-B in the liver or kidney. The decreases in triglyceride levels in wild-type mice 24 hrs after HS cannot be explained by the decreases in ACO and Cyp4a in liver and kidney. Thus, additional mechanisms that control triglyceride levels are likely at play including effects on mitochondrial fatty acid oxidation and require further study. In summary, 24 hrs after HS the decreases in triglyceride levels were PPARα-independent whereas decreases in glucose levels were PPARα-dependent.
Chaperone protein expression after heat shock and WY treatment
We assessed protein expression in the kidneys from the animals described in this study. Tcp1η, although not affected by HS in liver was induced up to 8-fold by HS in kidneys from both strains, indicating tissue-specific induction by HS (Figure 5D). Hsp86 but not Hsp25, Hsp70, and Hsp110 was induced by HS in the kidney (Figure 5D and data not shown). After WY exposure Hsp86, but not Tcp1η, Hsp25, Hsp70, or Hsp110 was weakly induced in the kidneys of wild-type, but not PPARα-null mice.
Most of the Hsp genes/proteins that were induced by WY (Hsp25, Hsp70, Hsp86, Hsp110, TCP1) carry out protein folding in the cytoplasm in a coordinated manner. For example, Hsp70 and Hsp110 cooperate to sequentially fold proteins in the cytoplasm . In contrast, we did not observe any mitochondrial chaperones such as mtHsp70 altered on the transcriptional level; a number of chaperones which operate in the endoplasmic reticulum including BiP/Grp78  and Grp94 (our study) were down-regulated. An earlier report showed that the Grp94 gene and protein were increased after exposure to WY or nafenopin . Overall, our studies indicate that WY coordinately regulates a number of chaperones which operate in the cytoplasm.
What might the increased levels of chaperone proteins be doing in the liver after WY exposure? Increased levels of chaperones might allow tight control of the inducibility of PPARα. Many nuclear receptors interact with chaperone proteins including the ones induced by WY in our studies . PPARα interacts with Hsp72  and is inhibited by Hsp90 . Thus, induction of Hsp90 and other family members may dampen the PPARα transcriptional response. Additionally, Hsp induction may help support the increases in protein synthesis associated with hepatocyte replication after PPC exposure. In addition, increased expression of Hsp family members has been associated with protection from apoptosis [42, 43] and PPC, at least under acute exposure conditions decrease basal levels of apoptosis . PPC exposure leads to dramatic increases in peroxisome size and number and increased expression of TCP1 subunits may be important for proper protein insertion into the peroxisomal membrane .
It is generally appreciated that increased expression of Hsps by mild HS can protect from subsequent challenges by more severe treatments . Our microarray and protein expression analysis showed that WY can induce the gene and protein expression of a number of chaperones. These same chaperones including Hsp25/27 [46, 47], Hsp70 [45, 48], Hsp86/90 [49–51] and Hsp105/110 [36, 52, 53] when over-expressed can protect cells from various stressors. Our results identify a number of potential candidate genes/proteins that could protect cells from environmental stressors by pretreatment with PPARα agonists.
Global decreases in mitochondrial gene expression in PPARα-null mice after heat shock are prevented by WY pretreatment
To identify candidates involved in resistance to HS-induced decreases in mitochondrial gene expression, we searched for genes regulated by WY independently of PPARα. Our previous experiments indicated that WY had minimal effects in the livers of PPARα-null mice, altering only up to ~8% of the total number of genes compared to wild-type mice . These results were based on transcript profiles generated using the Affymetrix gene chip U74Av2 which contains ~12,000 genes. Using full-genome arrays containing ~45,000 genes (Affymetrix MOE430_2 arrays), we examined gene expression in PPARα-null mice after exposure to WY for 6 hrs or 5 days. WY had minimal effects on gene expression after 6 hrs in PPARα-null mice (1 gene change) but at 5 days, 92 genes were significantly altered in PPARα-null mice. A number of the changes would affect intermediary metabolism and mitochondrial gene expression, the most striking of these changes being increases in PGC-1β (Ppargc1b) (~2-fold, p = 0.01; Figure 6C). (This gene is not included on the U74Av2 array.) Genes known to be regulated by PGC-1β included Cpt1[55–57] and a number of others involved in fatty acid oxidation including acetyl-Coenzyme A acyltransferase 1B (Acaa1b), aldehyde dehydrogenase family 3, subfamily A2 (Aldh3a2), acetyl-Coenzyme A acyltransferase 1A ( Acaa1a), enoyl coenzyme A hydratase 1, peroxisomal (Ech1) and peroxisomal delta3, delta2-enoyl-Coenzyme A isomerase (Peci) (Figure 6D).
PGC-1β is a master controller of hepatic energy homeostasis in part through regulation of mitochondrial oxidative metabolism genes. Overproduction of PGC-1β leads to increases in mitochondrial volume density in C2C12 myotubes , and defects in PGC-1β expression have been associated with global decreases in mitochondrial gene expression in the skeletal muscle of diabetics . PGC-1β physically interacts with transcription factors involved in mitochondrial biogenesis and lipid metabolism including Nrf-1, PPARα and PPARγ [55, 59, 60]. In trans-activation assays, WY activates PPARβ and PPARγ, albeit to a lesser extent than PPARα [61, 62]. Increases in PGC-1β co-activator function in the presence of the weak agonist function of WY may lead to increased expression of fatty acid oxidation genes regulated by PPARβ or PPARγ similar to what we have observed with the PPC perfluorooctanoic acid in PPARα-null mice . These results indicate that WY can alter the expression and activity of effectors of mitochondrial biogenesis and fatty acid oxidation, preventing HS from decreasing those gene batteries in PPARα-null mice. This protective effect of WY is reminiscent of thermotolerance induced by prior treatment with a low temperature HS .
Hsp gene regulation by PPARα and HSF1
The genes fell into 4 groups based on dependence on PPARα and HSF1 (Figure 7B). The largest groups were those that were PPARα-independent and were about equally split between those that were dependent or independent of HSF1. PPARα-dependent genes included those that were HSF1-dependent (Cacybp, Dnaja1, Fkbp4) or independent (Apoe, Irak1). This comparison reveals that although some genes were altered by HS independent of both PPARα and HSF1, most genes required one or both factors for HS regulation.
An important question is why does PPARα expression/activation determine in part the profile of HS genes? We hypothesize that PPARα is required for the induction of a subset of the HS genes. Class I and II genes (Figure 2B) require PPARα for HS regulation indicating that they could possess functional PPREs. Genes have been identified that possess HSEs that bind HSF1 but that were not HS-inducible. This provides evidence that other cis-acting elements or interacting factors were required for HS inducibility . A linkage between HSE and PPRE in the promoters of HS genes was found as binding sites for PPAR family members were significantly increased within the 1000 bp upstream of the genes that were up- or down-regulated by HS in a human epithelial cell line . Follow-up PPAR-HSF1 co-immunoprecipitation studies would reveal the nature of the interactions. Our analysis demonstrated major differences in the target genes of HS between wild-type and PPARα-null strains indicating that in the absence of PPARα, major changes in chromatin accessibility to HSF1 and other factors leads to greater sensitivity to HS alteration.
How are the HS genes regulated by PPC through PPARα? Given that Hsp gene expression is controlled in part by heat shock factor 1 (HSF1), one possibility is that the increases in Hsp gene expression are secondary to increases in the expression and activity of HSF1. However, we did not observe changes in Hsf1 expression in our transcript profiling studies and earlier studies showed that HSF1 and HSF2 expression and binding to HSE was not altered by WY exposure in the rat liver . To help determine whether regulation of Hsp gene expression is direct or indirect, we examined their promoters and found that only a few genes possess a putative PPRE(s) (data not shown) indicating indirect activation.
Many Hsp genes may be regulated indirectly through increases in oxidative stress. There is abundant evidence for the increased expression of chaperone gene expression under conditions that also increase oxidative stress [5, 45]. PPC exposure leads to increases in oxidative stress and lipid peroxidation mediated through increased activities of enzymes that generate reactive oxygen species . Although direct evidence that Hsp genes are activated by PPC-mediated increases in oxidative stress is lacking, the absolute increases in expression of some Hsp genes was higher in mice nullizygous for Nrf2, a transcription factor activated by oxidative stress that regulates genes that decrease oxidative stress. Thus, in the absence of Nrf2, increased levels of oxidative stress may have contributed to the greater increases in the Hsp genes by WY . However, we cannot rule out the induction of a HSP1 cofactor secondary to PPARα induction." In summary, PPARα may regulate the Hsp genes in a species-specific manner secondarily to increases in oxidative stress. Further work is needed to confirm this hypothesis.
Our microarray analysis of HS in the livers of mice uncovered a dramatic influence of PPARα on the overall gene expression pattern. There were remarkable differences in the transcriptional response to HS between wild-type and PPARα-null mice. A number of HS genes required PPARα for induction. Many mitochondrial genes were uniquely down-regulated in PPARα-null but not wild-type mice. Pretreatment of PPARα-null mice with WY prevented the down-regulation of these genes possibly through the increased expression of PGC-1β and regulated genes. Our findings support the hypothesis that PPARα in the liver acts as an integrator of chemical and physical stress signals manifested in the regulation of stress gene responses in the liver.
heat shock elements
heat shock factor
heat shock proteins
polymerase chain reaction
peroxisome proliferator-activated receptor
peroxisome proliferator chemical
peroxisome proliferator response element
We thank Drs. David Dix and Kevin Kregel for guidance in performing the animal exposures and western blots, Mr. Artie James, the CIIT Animal Care and Necropsy and Histology Units for assistance in performing these studies, Dr. Dennis House for assistance with statistical analyses, Mr. Lawrence Yoon for performing the RT-PCR analysis, Mr. John Havel for graphics support and Drs. David Dix, Stephen Edwards and Susan Hester for critical review of the manuscript.
The information in this document has been funded in part by the U.S. Environmental Protection Agency. It has been subjected to review by the National Health and Environmental Effects Research Laboratory and approved for publication. Approval does not signify that the contents reflect the views of the Agency, nor does mention of trade names or commercial products constitute endorsement or recommendation for use.
- Christians ES, Benjamin IJ: Heat shock response: lessons from mouse knockouts. Handbook of experimental pharmacology. 2006, 139-152. full_text. 172
- Pirkkala L, Nykanen P, Sistonen L: Roles of the heat shock transcription factors in regulation of the heat shock response and beyond. Faseb J. 2001, 15 (7): 1118-1131. 10.1096/fj00-0294rev.PubMedView Article
- Shamovsky I, Nudler E: New insights into the mechanism of heat shock response activation. Cell Mol Life Sci. 2008, 65 (6): 855-861. 10.1007/s00018-008-7458-y.PubMedView Article
- Slavotinek AM, Biesecker LG: Unfolding the role of chaperones and chaperonins in human disease. Trends Genet. 2001, 17 (9): 528-535. 10.1016/S0168-9525(01)02413-1.PubMedView Article
- Akerfelt M, Trouillet D, Mezger V, Sistonen L: Heat shock factors at a crossroad between stress and development. Annals of the New York Academy of Sciences. 2007, 1113: 15-27. 10.1196/annals.1391.005.PubMedView Article
- Arya R, Mallik M, Lakhotia SC: Heat shock genes - integrating cell survival and death. Journal of biosciences. 2007, 32 (3): 595-610. 10.1007/s12038-007-0059-3.PubMedView Article
- Page TJ, Sikder D, Yang L, Pluta L, Wolfinger RD, Kodadek T, Thomas RS: Genome-wide analysis of human HSF1 signaling reveals a transcriptional program linked to cellular adaptation and survival. Molecular bioSystems. 2006, 2 (12): 627-639. 10.1039/b606129j.PubMedView Article
- Trinklein ND, Murray JI, Hartman SJ, Botstein D, Myers RM: The role of heat shock transcription factor 1 in the genome-wide regulation of the mammalian heat shock response. Molecular biology of the cell. 2004, 15 (3): 1254-1261. 10.1091/mbc.E03-10-0738.PubMed CentralPubMedView Article
- Corton JC, Lapinskas PJ, Gonzalez FJ: Central role of PPARalpha in the mechanism of action of hepatocarcinogenic peroxisome proliferators. Mutation research. 2000, 448 (2): 139-151.PubMedView Article
- Klaunig JE, Babich MA, Baetcke KP, Cook JC, Corton JC, David RM, DeLuca JG, Lai DY, McKee RH, Peters JM: PPARalpha agonist-induced rodent tumors: modes of action and human relevance. Critical reviews in toxicology. 2003, 33 (6): 655-780. 10.1080/713608372.PubMedView Article
- Anderson SP, Dunn C, Laughter A, Yoon L, Swanson C, Stulnig TM, Steffensen KR, Chandraratna RA, Gustafsson JA, Corton JC: Overlapping transcriptional programs regulated by the nuclear receptors peroxisome proliferator-activated receptor alpha, retinoid × receptor, and liver × receptor in mouse liver. Molecular pharmacology. 2004, 66 (6): 1440-1452. 10.1124/mol.104.005496.PubMedView Article
- Anderson SP, Howroyd P, Liu J, Qian X, Bahnemann R, Swanson C, Kwak MK, Kensler TW, Corton JC: The transcriptional response to a peroxisome proliferator-activated receptor alpha agonist includes increased expression of proteome maintenance genes. The Journal of biological chemistry. 2004, 279 (50): 52390-52398. 10.1074/jbc.M409347200.PubMedView Article
- Peraza MA, Burdick AD, Marin HE, Gonzalez FJ, Peters JM: The toxicology of ligands for peroxisome proliferator-activated receptors (PPAR). Toxicol Sci. 2006, 90 (2): 269-295. 10.1093/toxsci/kfj062.PubMedView Article
- Corton JC, Brown-Borg HM: Peroxisome proliferator-activated receptor gamma coactivator 1 in caloric restriction and other models of longevity. The journals of gerontology. 2005, 60 (12): 1494-1509.PubMedView Article
- Mehendale HM: PPAR-alpha: a key to the mechanism of hepatoprotection by clofibrate. Toxicol Sci. 2000, 57 (2): 187-190. 10.1093/toxsci/57.2.187.PubMedView Article
- Chen C, Hennig GE, Whiteley HE, Corton JC, Manautou JE: Peroxisome proliferator-activated receptor alpha-null mice lack resistance to acetaminophen hepatotoxicity following clofibrate exposure. Toxicol Sci. 2000, 57 (2): 338-344. 10.1093/toxsci/57.2.338.PubMedView Article
- Corton JC, Apte U, Anderson SP, Limaye P, Yoon L, Latendresse J, Dunn C, Everitt JI, Voss KA, Swanson C: Mimetics of caloric restriction include agonists of lipid-activated nuclear receptors. The Journal of biological chemistry. 2004, 279 (44): 46204-46212. 10.1074/jbc.M406739200.PubMedView Article
- Anderson SP, Yoon L, Richard EB, Dunn CS, Cattley RC, Corton JC: Delayed liver regeneration in peroxisome proliferator-activated receptor-alpha-null mice. Hepatology (Baltimore, Md). 2002, 36 (3): 544-554. 10.1053/jhep.2002.35276.View Article
- Portilla D, Dai G, Peters JM, Gonzalez FJ, Crew MD, Proia AD: Etomoxir-induced PPARalpha-modulated enzymes protect during acute renal failure. American journal of physiology. 2000, 278 (4): F667-675.PubMed
- Sivarajah A, Chatterjee PK, Hattori Y, Brown PA, Stewart KN, Todorovic Z, Mota-Filipe H, Thiemermann C: Agonists of peroxisome-proliferator activated receptor-alpha (clofibrate and WY14643) reduce renal ischemia/reperfusion injury in the rat. Med Sci Monit. 2002, 8 (12): BR532-539.PubMed
- 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. Molecular and cellular biology. 1995, 15 (6): 3012-3022.PubMed CentralPubMedView Article
- Wilson CL, Miller CJ: Simpleaffy: a BioConductor package for Affymetrix Quality Control and data analysis. Bioinformatics (Oxford, England). 2005, 21 (18): 3683-3685. 10.1093/bioinformatics/bti605.View Article
- Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proceedings of the National Academy of Sciences of the United States of America. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.PubMed CentralPubMedView Article
- Sanderson LM, de Groot PJ, Hooiveld GJ, Koppen A, Kalkhoven E, Muller M, Kersten S: Effect of synthetic dietary triglycerides: a novel research paradigm for nutrigenomics. PLoS ONE. 2008, 3 (2): e1681-10.1371/journal.pone.0001681.PubMed CentralPubMedView Article
- Jonnalagadda S, Srinivasan R: Principal components analysis based methodology to identify differentially expressed genes in time-course microarray data. BMC bioinformatics. 2008, 9: 267-10.1186/1471-2105-9-267.PubMed CentralPubMedView Article
- Fan LQ, Coley J, Miller RT, Cattley RC, Corton JC: Opposing mechanisms of NADPH-cytochrome P450 oxidoreductase regulation by peroxisome proliferators. Biochemical pharmacology. 2003, 65 (6): 949-959. 10.1016/S0006-2952(03)00004-2.PubMedView Article
- Corton JC, Fan LQ, Brown S, Anderson SP, Bocos C, Cattley RC, Mode A, Gustafsson JA: Down-regulation of cytochrome P450 2C family members and positive acute-phase response gene expression by peroxisome proliferator chemicals. Molecular pharmacology. 1998, 54 (3): 463-473.PubMed
- Dinh HK, Zhao B, Schuschereba ST, Merrill G, Bowman PD: Gene expression profiling of the response to thermal injury in human cells. Physiological genomics. 2001, 7 (1): 3-13.PubMed
- Yoshida H: Unconventional splicing of XBP-1 mRNA in the unfolded protein response. Antioxidants & redox signaling. 2007, 9 (12): 2323-2333.View Article
- Patti ME, Butte AJ, Crunkhorn S, Cusi K, Berria R, Kashyap S, Miyazaki Y, Kohane I, Costello M, Saccone R: Coordinated reduction of genes of oxidative metabolism in humans with insulin resistance and diabetes: Potential role of PGC1 and NRF1. Proceedings of the National Academy of Sciences of the United States of America. 2003, 100 (14): 8466-8471. 10.1073/pnas.1032913100.PubMed CentralPubMedView Article
- Spiegelman BM: Transcriptional control of energy homeostasis through the PGC1 coactivators. Novartis Foundation symposium. 2007, 286: 3-6. full_text. discusssion 6-12, 162-163, 196-203PubMedView Article
- Chai J, Tarnawski AS: Serum response factor: discovery, biochemistry, biological roles and implications for tissue injury healing. J Physiol Pharmacol. 2002, 53 (2): 147-157.PubMed
- Miano JM, Long X, Fujiwara K: Serum response factor: master regulator of the actin cytoskeleton and contractile apparatus. Am J Physiol Cell Physiol. 2007, 292 (1): C70-81. 10.1152/ajpcell.00386.2006.PubMedView Article
- Cho MC, Lee K, Paik SG, Yoon DY: Peroxisome Proliferators-Activated Receptor (PPAR) Modulators and Metabolic Disorders. PPAR research. 2008, 2008: 679137-10.1155/2008/679137.PubMed CentralPubMedView Article
- Fievet C, Fruchart JC, Staels B: PPARalpha and PPARgamma dual agonists for the treatment of type 2 diabetes and the metabolic syndrome. Current opinion in pharmacology. 2006, 6 (6): 606-614. 10.1016/j.coph.2006.06.009.PubMedView Article
- Easton DP, Kaneko Y, Subjeck JR: The hsp110 and Grp1 70 stress proteins: newly recognized relatives of the Hsp70s. Cell stress & chaperones. 2000, 5 (4): 276-290.View Article
- Motojima K, Goto S: Rat liver BiP/GRP78 is down-regulated by a peroxisome-proliferator, clofibrate. FEBS letters. 1992, 308 (2): 207-210. 10.1016/0014-5793(92)81275-Q.PubMedView Article
- Macdonald N, Barrow K, Tonge R, Davison M, Roberts RA, Chevalier S: PPARalpha-dependent alteration of GRP94 expression in mouse hepatocytes. Biochemical and biophysical research communications. 2000, 277 (3): 699-704. 10.1006/bbrc.2000.3741.PubMedView Article
- Morimoto RI: Dynamic remodeling of transcription complexes by molecular chaperones. Cell. 2002, 110 (3): 281-284. 10.1016/S0092-8674(02)00860-7.PubMedView Article
- Huang Q, Alvares K, Chu R, Bradfield CA, Reddy JK: Association of peroxisome proliferator-activated receptor and Hsp72. The Journal of biological chemistry. 1994, 269 (11): 8493-8497.PubMed
- Sumanasekera WK, Tien ES, Davis JW, Turpey R, Perdew GH, Vanden Heuvel JP: Heat shock protein-90 (Hsp90) acts as a repressor of peroxisome proliferator-activated receptor-alpha (PPARalpha) and PPARbeta activity. Biochemistry. 2003, 42 (36): 10726-10735. 10.1021/bi0347353.PubMedView Article
- Beere HM: Stressed to death: regulation of apoptotic signaling pathways by the heat shock proteins. Sci STKE. 2001, 2001 (93): RE1-10.1126/stke.2001.93.re1.PubMed
- Latchman DS: Heat shock proteins and cardiac protection. Cardiovascular research. 2001, 51 (4): 637-646. 10.1016/S0008-6363(01)00354-6.PubMedView Article
- Pause B, Diestelkotter P, Heid H, Just WW: Cytosolic factors mediate protein insertion into the peroxisomal membrane. FEBS letters. 1997, 414 (1): 95-98. 10.1016/S0014-5793(97)00975-7.PubMedView Article
- Benjamin IJ, McMillan DR: Stress (heat shock) proteins: molecular chaperones in cardiovascular biology and disease. Circulation research. 1998, 83 (2): 117-132.PubMedView Article
- Martin JL, Mestril R, Hilal-Dandan R, Brunton LL, Dillmann WH: Small heat shock proteins and protection against ischemic injury in cardiac myocytes. Circulation. 1997, 96 (12): 4343-4348.PubMedView Article
- Park YM, Han MY, Blackburn RV, Lee YJ: Overexpression of HSP25 reduces the level of TNF alpha-induced oxidative DNA damage biomarker, 8-hydroxy-2'-deoxyguanosine, in L929 cells. Journal of cellular physiology. 1998, 174 (1): 27-34. 10.1002/(SICI)1097-4652(199801)174:1<27::AID-JCP4>3.0.CO;2-I.PubMedView Article
- Tolson JK, Dix DJ, Voellmy RW, Roberts SM: Increased hepatotoxicity of acetaminophen in Hsp70i knockout mice. Toxicology and applied pharmacology. 2006, 210 (1-2): 157-162. 10.1016/j.taap.2005.10.001.PubMedView Article
- Conconi M, Petropoulos I, Emod I, Turlin E, Biville F, Friguet B: Protection from oxidative inactivation of the 20S proteasome by heat-shock protein 90. The Biochemical journal. 1998, 333 (Pt 2): 407-415.PubMed CentralPubMedView Article
- Cumming DV, Heads RJ, Watson A, Latchman DS, Yellon DM: Differential protection of primary rat cardiocytes by transfection of specific heat stress proteins. Journal of molecular and cellular cardiology. 1996, 28 (12): 2343-2349. 10.1006/jmcc.1996.0227.PubMedView Article
- Lee MW, Park SC, Chae HS, Bach JH, Lee HJ, Lee SH, Kang YK, Kim KY, Lee WB, Kim SS: The protective role of HSP90 against 3-hydroxykynurenine-induced neuronal apoptosis. Biochemical and biophysical research communications. 2001, 284 (2): 261-267. 10.1006/bbrc.2001.4938.PubMedView Article
- Hatayama T, Yamagishi N, Minobe E, Sakai K: Role of hsp105 in protection against stress-induced apoptosis in neuronal PC12 cells. Biochemical and biophysical research communications. 2001, 288 (3): 528-534. 10.1006/bbrc.2001.5802.PubMedView Article
- Oh HJ, Chen X, Subjeck JR: Hsp110 protects heat-denatured proteins and confers cellular thermoresistance. The Journal of biological chemistry. 1997, 272 (50): 31636-31640. 10.1074/jbc.272.50.31636.PubMedView Article
- Rosen MB, Lee JS, Ren H, Vallanat B, Liu J, Waalkes MP, Abbott BD, Lau C, Corton JC: Toxicogenomic dissection of the perfluorooctanoic acid transcript profile in mouse liver: evidence for the involvement of nuclear receptors PPAR alpha and CAR. Toxicol Sci. 2008, 103 (1): 46-56. 10.1093/toxsci/kfn025.PubMedView Article
- Lin J, Tarr PT, Yang R, Rhee J, Puigserver P, Newgard CB, Spiegelman BM: PGC-1beta in the regulation of hepatic glucose and energy metabolism. The Journal of biological chemistry. 2003, 278 (33): 30843-30848. 10.1074/jbc.M303643200.PubMedView Article
- Sadana P, Zhang Y, Song S, Cook GA, Elam MB, Park EA: Regulation of carnitine palmitoyltransferase I (CPT-Ialpha) gene expression by the peroxisome proliferator activated receptor gamma coactivator (PGC-1) isoforms. Molecular and cellular endocrinology. 2007, 267 (1-2): 6-16. 10.1016/j.mce.2006.11.012.PubMed CentralPubMedView Article
- Wolfrum C, Stoffel M: Coactivation of Foxa2 through Pgc-1beta promotes liver fatty acid oxidation and triglyceride/VLDL secretion. Cell metabolism. 2006, 3 (2): 99-110. 10.1016/j.cmet.2006.01.001.PubMedView Article
- St-Pierre J, Lin J, Krauss S, Tarr PT, Yang R, Newgard CB, Spiegelman BM: Bioenergetic analysis of peroxisome proliferator-activated receptor gamma coactivators 1alpha and 1beta (PGC-1alpha and PGC-1beta) in muscle cells. The Journal of biological chemistry. 2003, 278 (29): 26597-26603. 10.1074/jbc.M301850200.PubMedView Article
- Kressler D, Schreiber SN, Knutti D, Kralli A: The PGC-1-related protein PERC is a selective coactivator of estrogen receptor alpha. The Journal of biological chemistry. 2002, 277 (16): 13918-13925. 10.1074/jbc.M201134200.PubMedView Article
- Lin J, Puigserver P, Donovan J, Tarr P, Spiegelman BM: Peroxisome proliferator-activated receptor gamma coactivator 1beta (PGC-1beta ), a novel PGC-1-related transcription coactivator associated with host cell factor. The Journal of biological chemistry. 2002, 277 (3): 1645-1648. 10.1074/jbc.C100631200.PubMedView Article
- Corton JC, Lapinskas PJ: Peroxisome proliferator-activated receptors: mediators of phthalate ester-induced effects in the male reproductive tract?. Toxicol Sci. 2005, 83 (1): 4-17. 10.1093/toxsci/kfi011.PubMedView Article
- Lapinskas PJ, Brown S, Leesnitzer LM, Blanchard S, Swanson C, Cattley RC, Corton JC: Role of PPARalpha in mediating the effects of phthalates and metabolites in the liver. Toxicology. 2005, 207 (1): 149-163. 10.1016/j.tox.2004.09.008.PubMedView Article
- Laramie JM, Chung TP, Brownstein B, Stormo GD, Cobb JP: Transcriptional profiles of human epithelial cells in response to heat: computational evidence for novel heat shock proteins. Shock (Augusta, Ga). 2008, 29 (5): 623-630. 10.1097/SHK.0b013e318157f33c.
- O'Brien ML, Cunningham ML, Spear BT, Glauert HP: Peroxisome proliferators do not activate the transcription factors AP-1, early growth response-1, or heat shock factors 1 and 2 in rats or hamsters. Toxicol Sci. 2002, 69 (1): 139-48. 10.1093/toxsci/69.1.139.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.