Analysis of the Peroxisome Proliferator-Activated Receptor-Β/Δ
Total Page:16
File Type:pdf, Size:1020Kb
Khozoie et al. BMC Genomics 2012, 13:665 http://www.biomedcentral.com/1471-2164/13/665 RESEARCH ARTICLE Open Access Analysis of the peroxisome proliferator-activated receptor-β/δ (PPARβ/δ) cistrome reveals novel co-regulatory role of ATF4 Combiz Khozoie1, Michael G Borland1,4, Bokai Zhu1, Songjoon Baek2, Sam John2,5, Gordon L Hager2, Yatrik M Shah3,6, Frank J Gonzalez3 and Jeffrey M Peters1* Abstract Background: The present study coupled expression profiling with chromatin immunoprecipitation sequencing (ChIP-seq) to examine peroxisome proliferator-activated receptor-β/δ (PPARβ/δ)-dependent regulation of gene expression in mouse keratinocytes, a cell type that expresses PPARβ/δ in high concentration. Results: Microarray analysis elucidated eight different types of regulation that modulated PPARβ/δ-dependent gene expression of 612 genes ranging from repression or activation without an exogenous ligand, repression or activation with an exogenous ligand, or a combination of these effects. Bioinformatic analysis of ChIP-seq data demonstrated promoter occupancy of PPARβ/δ for some of these genes, and also identified the presence of other transcription factor binding sites in close proximity to PPARβ/δ bound to chromatin. For some types of regulation, ATF4 is required for ligand-dependent induction of PPARβ/δ target genes. Conclusions: PPARβ/δ regulates constitutive expression of genes in keratinocytes, thus suggesting the presence of one or more endogenous ligands. The diversity in the types of gene regulation carried out by PPARβ/δ is consistent with dynamic binding and interactions with chromatin and indicates the presence of complex regulatory networks in cells expressing high levels of this nuclear receptor such as keratinocytes. Results from these studies are the first to demonstrate that differences in DNA binding of other transcription factors can directly influence the transcriptional activity of PPARβ/δ. Keywords: Peroxisome proliferator-activated receptor-β/δ, Gene expression, Keratinocytes Background PPARβ/δ for the prevention and/or treatment of diverse Peroxisome proliferator-activated receptor-β/δ (PPARβ/ diseases, including cancer, diabetes, metabolic syndrome, δ) is a ligand activated transcription factor with particu- and dyslipidemias [5-10,12]. larly high abundance in small intestine, colon, liver, and PPARβ/δ modulates cellular function by regulating skin [1-4]. In these tissues, PPARβ/δ likely has a consti- gene expression through several mechanisms. For ex- tutive physiological role, possibly modulated by the pres- ample, PPARβ/δ can interact and bind with other tran- ence of an endogenous ligand. This is consistent with scription factors, including NFκB, ERK5, and STAT3, the fact that PPARβ/δ is required to mediate the induc- and attenuate their signaling [5-10]. PPARβ/δ can also tion of terminal differentiation in epithelial cells [5-10]. repress gene expression by dynamically binding to chro- PPARβ/δ also attenuates inflammation and regulates matin in association with co-repressors [13,14]. The glucose and lipid homeostasis [5-11]. These important most commonly described mechanism by which PPARβ/ physiological roles underscore the potential for targeting δ was thought to regulate gene expression is that bind- ing of ligand to receptor induces a conformational * Correspondence: [email protected] change in the protein. This change in structure is 1 Department of Veterinary and Biomedical Sciences and The Center for accompanied by the release of co-repressors, heterodimer- Molecular Toxicology and Carcinogenesis, The Pennsylvania State University, University Park, PA 16802, USA ization with the retinoid X receptor (RXR), recruitment of Full list of author information is available at the end of the article © 2012 Khozoie et al.; licensee 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. Khozoie et al. BMC Genomics 2012, 13:665 Page 2 of 19 http://www.biomedcentral.com/1471-2164/13/665 co-activators and RNA polymerase II and increased tran- Results scription when this complex is bound to peroxisome PPARβ/δ-dependent regulation of genes in keratinocytes proliferator response elements (PPREs) proximal to tar- Six hundred and twelve genes were identified by get genes. However, more recent studies indicate that genome-wide expression profiling that were differentially nuclear receptors actually regulate gene expression regulated by either ligand, disruption of PPARβ/δ,or through dynamic and sometimes transient interactions both (Figure 1A, Table 1, Additional file 1: Table S2). with chromatin rather than through static complexes These genes were categorized into four major response occupying chromatin [15-18]. This suggests that the types: (I) repression without exogenous ligand (n=185), dynamic and transient occupancy of PPARβ/δ on chro- (II) activation without exogenous ligand (n=297), (III) matin must be considered when attempting to interpret activation with exogenous ligand (n=71), and (IV) receptor-DNA binding studies [15-18]. There are many repression with exogenous ligand (n=28). qPCR con- different classes of enzymes that can modify nucleosome firmed these four common types of PPARβ/δ-dependent location and/or chromatin structure that will allow for changes in gene expression detected by microarray ana- nuclear receptors to dynamically and transiently bind to lysis (Figure 2). Combined, these comprised 94.9% of the chromatin, and this binding can be markedly changed in genes identified that were differentially regulated by ei- response to the presence of endogenous and exogenous ther ligand, disruption of PPARβ/δ, or both. An add- ligands [15-18]. Further, there is evidence that the dy- itional four response types exhibiting combined namic interactions between agonist-activated nuclear responses were also observed: (V) repression without receptors with their co-regulators and DNA binding sites exogenous ligand and activation with exogenous ligand can lead to highly variable responses including differ- (n=12), (VI) activation without exogenous ligand and ences in: 1) the regions of chromatin occupied by the re- repression with exogenous ligand (n=9), (VII) activation ceptor, 2) expression of target genes and 3) the ultimate with and without exogenous ligand (n=1), and (VIII) biological effect. Conformational changes in a receptor repression with and without exogenous ligand (n=9) caused by different agonists, antagonists and partial ago- (Figure 1A, Table 1). The latter four response types only nists can cause differential recruitment of co-regulators comprised ~5% of the genes modulated by GW0742 resulting in alteration of the dynamics of transcriptional and/or disruption of PPARβ/δ. complexes and interactions with DNA binding sites. Characterization of the PPARβ/δ target genes and Thus, there are multiple levels of regulation by which their transcriptional responses was undertaken by func- PPARβ/δ can influence the expression of target genes. tional category enrichment analysis [27,28]. 50% of the Characterization of bona fide target genes directly enriched functional categories were common between regulated by PPARβ/δ in all tissues is incomplete. type I and type II responses (Figure 1B). Genes that Angiopoietin-like 4 (Angptl4) and adipocyte differentiation- regulate fatty acid metabolism were common between related protein (Adrp) are two genes that can be modu- type III and type V responses (Figure 1B). The level of lated by PPARβ/δ by direct transcriptional regulation. gene expression observed following ligand activation was That Angptl4 and Adrp are direct target genes is based compared to changes in gene expression observed by on analyses demonstrating: 1) that functional PPREs exist disruption of PPARβ/δ (Figure 1C). This analysis proximal to these genes, and 2) confirmed promoter revealed clustering for the eight different response types, occupancy of PPARβ/δ following ligand activation [19-21]. but there were differences in the magnitude of change In addition to these experimentally confirmed PPARβ/δ found for each of the response types (Figure 1C). The target genes, in silico screening based on genomic PPRE conventional view accounting for combined responses frequency predicted as many as 4000 to 5000 targets for involving activation/repression in both the presence and PPARs in the human genome [22-24] and a compar- absence of exogenous ligand (types V and VI, e.g. able number of binding sites for PPARs in some cells Lpcat3) suggests that a ligand-mediated switch occurs have been confirmed by chromatin immunoprecipitation between repression and activation (type V) or vice-versa sequencing (ChIP-seq) [25]. This suggests that many (type VI) [13]. For example, in the absence of PPARβ/δ, PPARβ/δ binding sites and target genes remain unidenti- expression of Lpcat3 is enhanced because PPARβ/δ fied. Moreover, recent studies have used next genera- represses expression, whereas ligand activation of PPARβ/δ tion sequencing to elucidate novel regulatory roles for increases expression of Lpcat3. This phenomenon was PPARβ/δ. For example, ChIP-seq was used to demon- also observed with type VI genes, consistent with a ligand- strate that PPARβ/δ and PPARγ can be exchanged on mediated release of an activating complex. However, it is target gene in adipocytes, following