Species-specific strategies underlying conserved functions of metabolic transcription factors - PubMed (original) (raw)
Species-specific strategies underlying conserved functions of metabolic transcription factors
Raymond E Soccio et al. Mol Endocrinol. 2011 Apr.
Abstract
The winged helix protein FOXA2 and the nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ) are highly conserved, regionally expressed transcription factors (TFs) that regulate networks of genes controlling complex metabolic functions. Cistrome analysis for Foxa2 in mouse liver and PPARγ in mouse adipocytes has previously produced consensus-binding sites that are nearly identical to those used by the corresponding TFs in human cells. We report here that, despite the conservation of the canonical binding motif, the great majority of binding regions for FOXA2 in human liver and for PPARγ in human adipocytes are not in the orthologous locations corresponding to the mouse genome, and vice versa. Of note, TF binding can be absent in one species despite sequence conservation, including motifs that do support binding in the other species, demonstrating a major limitation of in silico binding site prediction. Whereas only approximately 10% of binding sites are conserved, gene-centric analysis reveals that about 50% of genes with nearby TF occupancy are shared across species for both hepatic FOXA2 and adipocyte PPARγ. Remarkably, for both TFs, many of the shared genes function in tissue-specific metabolic pathways, whereas species-unique genes fail to show enrichment for these pathways. Nonetheless, the species-unique genes, like the shared genes, showed the expected transcriptional regulation by the TFs in loss-of-function experiments. Thus, species-specific strategies underlie the biological functions of metabolic TFs that are highly conserved across mammalian species. Analysis of factor binding in multiple species may be necessary to distinguish apparent species-unique noise and reveal functionally relevant information.
Figures
Fig. 1.
Conservation of TF-binding regions. The overall conservation of binding regions approximates 10% between mouse and human. All FOXA2 (A) and PPARγ (B) binding regions are represented in Venn diagrams. The darker-colored pie slices of each circle show the regions that fail to convert to the other genome, whereas the intersection of circles shows the conservation of convertible regions based on centers within 500 bp. There are different numbers of conserved regions in each species, reflecting cases in which a single GLITR region in one species overlaps multiple regions in the other species. C, Conserved regions are more likely to be located in gene promoters. Each dataset of binding regions was divided as in the Venn diagrams, and the percentage of regions localized within 5 kb of the transcription start site of the nearest gene is shown. D, The strongest regions are more likely to be conserved. The percent conservation of the 100 strongest convertible binding regions for each dataset is shown. E, Conservation is found even among the weakest binding regions. All convertible binding regions were divided into quintiles based on strength, and the percent conservation is shown. The data are clearly skewed from the null hypothesis of 20% (the P value for the overall skew in each dataset is <10−10), but even the weakest quintile accounts for more than 10% of all conserved sites. *, P < 0.05; **, P < 0.001; ***, P <10−10 vs. all sites.
Fig. 2.
Motifs in TF-binding regions. A, The complete datasets of mouse and human binding regions for FOXA2 and PPARγ were searched using de novo motif analysis, and the top scoring motif is shown for each. These motifs were compared with the TRANSFAC database to find the most similar motif. B, Conserved and species-specific binding correlates with the occurrence of motifs in binding regions. The top human de novo motif for FOXA2 (left) and PPARγ (right) was used to interrogate the convertible binding regions as well as control datasets matched for length and guanine and cytosine content. The binding regions were divided into conserved and species-specific regions, and coordinates were mapped on the mouse (mm8) and human (hg18) genomes as indicated. *, P < 0.01; **, P <10−8; ***, P <10−70 for comparisons as indicated. #, P <10−70 vs. mouse-only regions mapped on human genome, ^, P <10−70 vs. human-only regions mapped on mouse genome. C, PhastCons scores for conserved and species-specific (nonconserved) binding. PhastCons scores were calculated based on placental mammal sequences and were calculated in human binding regions. Conserved binding regions (yellow) had high conservation scores, while non-conserved regions (green) had much lower conservation scores, but were still higher than all regions (red). Non-convertible (blue) and matched control (purple) regions showed background conservation.
Fig. 3.
Shared genes nearest to TF-binding regions. A, The sharing of nearest genes approximates 50%. The nearest gene was found to each binding region, and this list was filtered to unique genes with orthologs in the other species. The intersection of the Venn diagrams represents orthologous genes with nearest binding regions in both species. B, Most shared genes have only nonconserved nearest regions. For each list of genes nearest to factor-binding regions in one species, about 4% of genes did not have an ortholog in the other species (red). When there was an ortholog, it either had no nearest region (yellow), only nonconserved nearest regions (dark blue), or a conserved nearest region (light blue). C, Shared genes have more nearest binding regions. The genes are divided as in B with the average number of nearest regions shown. *, P <10−6 by permutation t test, compared with 106 random permutations of the data.
Fig. 4.
Transcriptional regulation of shared and species-unique target genes by FOXA2 and PPARγ. A, Microarray experiments were performed comparing liver gene expression in Foxa1/a2 double-mutant mice to control mice (FOXA2) and 3T3-L1 mouse adipocyte gene expression in PPARγ siRNA knockdown cells vs. control siRNA (PPARγ). All nonredundant genes on each microarray were divided into the four classes shown based on the presence or absence of nearest factor-binding regions in the mouse and human genomes. For each class, the percentage of genes with statistically significant more than 1.5-fold microarray down-regulation upon factor knockout/knockdown was determined. *, P < 10−9; **, _P_ <10−40; NS, _P_ > 0.1 vs. negative control (black bar, genes with binding regions in neither species). B, UCSC browser tracks on the mouse genome (mm8) are shown to illustrate several classes of genes, with tracks shown for the mouse and convertible human PPARγ-binding regions. The fold down-regulation upon PPARγ siRNA knockdown from the microarray experiment in panel A is shown for each gene.
Fig. 5.
Conservation of metabolic pathways regulated by TFs. A, DAVID/PANTHER biological pathway analysis of FOXA2 nearest genes that are human unique (red), mouse-unique (blue), or shared human (purple). B, DAVID/PANTHER biological pathway analysis of PPARγ nearest genes that are human unique (green), mouse unique (yellow), or shared human (light green). All gene lists were generated by the analysis in Fig. 3A, and all pathways with P < 0.01 are shown.
Similar articles
- Cell-specific determinants of peroxisome proliferator-activated receptor gamma function in adipocytes and macrophages.
Lefterova MI, Steger DJ, Zhuo D, Qatanani M, Mullican SE, Tuteja G, Manduchi E, Grant GR, Lazar MA. Lefterova MI, et al. Mol Cell Biol. 2010 May;30(9):2078-89. doi: 10.1128/MCB.01651-09. Epub 2010 Feb 22. Mol Cell Biol. 2010. PMID: 20176806 Free PMC article. - De-novo identification of PPARgamma/RXR binding sites and direct targets during adipogenesis.
Hamza MS, Pott S, Vega VB, Thomsen JS, Kandhadayar GS, Ng PW, Chiu KP, Pettersson S, Wei CL, Ruan Y, Liu ET. Hamza MS, et al. PLoS One. 2009;4(3):e4907. doi: 10.1371/journal.pone.0004907. Epub 2009 Mar 20. PLoS One. 2009. PMID: 19300518 Free PMC article. - Adipocyte-specific expression of murine resistin is mediated by synergism between peroxisome proliferator-activated receptor gamma and CCAAT/enhancer-binding proteins.
Tomaru T, Steger DJ, Lefterova MI, Schupp M, Lazar MA. Tomaru T, et al. J Biol Chem. 2009 Mar 6;284(10):6116-25. doi: 10.1074/jbc.M808407200. Epub 2009 Jan 5. J Biol Chem. 2009. PMID: 19126543 Free PMC article. - Identification of DNA-dependent protein kinase as a cofactor for the forkhead transcription factor FoxA2.
Nock A, Ascano JM, Jones T, Barrero MJ, Sugiyama N, Tomita M, Ishihama Y, Malik S. Nock A, et al. J Biol Chem. 2009 Jul 24;284(30):19915-26. doi: 10.1074/jbc.M109.016295. Epub 2009 May 28. J Biol Chem. 2009. PMID: 19478084 Free PMC article. - Genome-wide analyses of ChIP-Seq derived FOXA2 DNA occupancy in liver points to genetic networks underpinning multiple complex traits.
Johnson ME, Schug J, Wells AD, Kaestner KH, Grant SF. Johnson ME, et al. J Clin Endocrinol Metab. 2014 Aug;99(8):E1580-5. doi: 10.1210/jc.2013-4503. Epub 2014 May 30. J Clin Endocrinol Metab. 2014. PMID: 24878043 Free PMC article.
Cited by
- A single-cell CRISPRi platform for characterizing candidate genes relevant to metabolic disorders in human adipocytes.
Bielczyk-Maczynska E, Sharma D, Blencowe M, Saliba Gustafsson P, Gloudemans MJ, Yang X, Carcamo-Orive I, Wabitsch M, Svensson KJ, Park CY, Quertermous T, Knowles JW, Li J. Bielczyk-Maczynska E, et al. Am J Physiol Cell Physiol. 2023 Sep 1;325(3):C648-C660. doi: 10.1152/ajpcell.00148.2023. Epub 2023 Jul 24. Am J Physiol Cell Physiol. 2023. PMID: 37486064 Free PMC article. - Natural human genetic variation determines basal and inducible expression of PM20D1, an obesity-associated gene.
Benson KK, Hu W, Weller AH, Bennett AH, Chen ER, Khetarpal SA, Yoshino S, Bone WP, Wang L, Rabinowitz JD, Voight BF, Soccio RE. Benson KK, et al. Proc Natl Acad Sci U S A. 2019 Nov 12;116(46):23232-23242. doi: 10.1073/pnas.1913199116. Epub 2019 Oct 28. Proc Natl Acad Sci U S A. 2019. PMID: 31659023 Free PMC article. - Patient Adipose Stem Cell-Derived Adipocytes Reveal Genetic Variation that Predicts Antidiabetic Drug Response.
Hu W, Jiang C, Guan D, Dierickx P, Zhang R, Moscati A, Nadkarni GN, Steger DJ, Loos RJF, Hu C, Jia W, Soccio RE, Lazar MA. Hu W, et al. Cell Stem Cell. 2019 Feb 7;24(2):299-308.e6. doi: 10.1016/j.stem.2018.11.018. Epub 2019 Jan 10. Cell Stem Cell. 2019. PMID: 30639037 Free PMC article. - Diet-Induced Circadian Enhancer Remodeling Synchronizes Opposing Hepatic Lipid Metabolic Processes.
Guan D, Xiong Y, Borck PC, Jang C, Doulias PT, Papazyan R, Fang B, Jiang C, Zhang Y, Briggs ER, Hu W, Steger D, Ischiropoulos H, Rabinowitz JD, Lazar MA. Guan D, et al. Cell. 2018 Aug 9;174(4):831-842.e12. doi: 10.1016/j.cell.2018.06.031. Epub 2018 Jul 26. Cell. 2018. PMID: 30057115 Free PMC article. - PPARγ is dispensable for clear cell renal cell carcinoma progression.
Sanchez DJ, Steger DJ, Skuli N, Bansal A, Simon MC. Sanchez DJ, et al. Mol Metab. 2018 Aug;14:139-149. doi: 10.1016/j.molmet.2018.05.013. Epub 2018 May 21. Mol Metab. 2018. PMID: 29866440 Free PMC article.
References
- Wasserman WW, Palumbo M, Thompson W, Fickett JW, Lawrence CE. 2000. Human-mouse genome comparisons to locate regulatory sites. Nat Genet 26:225–228 - PubMed
- Bejerano G, Pheasant M, Makunin I, Stephen S, Kent WJ, Mattick JS, Haussler D. 2004. Ultraconserved elements in the human genome. Science 304:1321–1325 - PubMed
Publication types
MeSH terms
Substances
Grants and funding
- U19DK/HL/ES 62434/DK/NIDDK NIH HHS/United States
- R01 DK049780/DK/NIDDK NIH HHS/United States
- P01 DK049210/DK/NIDDK NIH HHS/United States
- T32 DK007314/DK/NIDDK NIH HHS/United States
- P30 DK019525/DK/NIDDK NIH HHS/United States
LinkOut - more resources
Full Text Sources
Miscellaneous