Species-specific strategies underlying conserved functions of metabolic transcription factors.

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.

[1]  J. Stone,et al.  Rapid evolution of cis-regulatory sequences via local point mutations. , 2001, Molecular biology and evolution.

[2]  K. Lindblad-Toh,et al.  Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals , 2005, Nature.

[3]  M. Lazar,et al.  HIV Protease Inhibitor‐Specific Alterations in Human Adipocyte Differentiation and Metabolism , 2006, Obesity.

[4]  H. Stunnenberg,et al.  Genome-wide profiling of PPARgamma:RXR and RNA polymerase II occupancy reveals temporal activation of distinct metabolic pathways and changes in RXR dimer composition during adipogenesis. , 2008, Genes & development.

[5]  G. Bourque,et al.  Transposable elements have rewired the core regulatory network of human embryonic stem cells , 2010, Nature Genetics.

[6]  N. Patel,et al.  Functional analysis of eve stripe 2 enhancer evolution in Drosophila: rules governing conservation and change. , 1998, Development.

[7]  K. Kaestner,et al.  The Foxa family of transcription factors in development and metabolism , 2006, Cellular and Molecular Life Sciences CMLS.

[8]  Cory Y. McLean,et al.  GREAT improves functional interpretation of cis-regulatory regions , 2010, Nature Biotechnology.

[9]  Klaus H. Kaestner,et al.  The initiation of liver development is dependent on Foxa transcription factors , 2005, Nature.

[10]  J. Costas,et al.  Turnover of binding sites for transcription factors involved in early Drosophila development. , 2003, Gene.

[11]  Jacques van Helden,et al.  Evaluation of phylogenetic footprint discovery for predicting bacterial cis-regulatory elements and revealing their evolution , 2008, BMC Bioinformatics.

[12]  D. Gifford,et al.  Tissue-specific transcriptional regulation has diverged significantly between human and mouse , 2007, Nature Genetics.

[13]  David A. Nix,et al.  Large-Scale Turnover of Functional Transcription Factor Binding Sites in Drosophila , 2006, PLoS Comput. Biol..

[14]  Michael D. Wilson,et al.  Five-Vertebrate ChIP-seq Reveals the Evolutionary Dynamics of Transcription Factor Binding , 2010, Science.

[15]  Jonathan Schug,et al.  PPARgamma and C/EBP factors orchestrate adipocyte biology via adjacent binding on a genome-wide scale. , 2008, Genes & development.

[16]  M. Ludwig,et al.  Functional evolution of noncoding DNA. , 2002, Current opinion in genetics & development.

[17]  J. Fickett,et al.  Identification of regulatory regions which confer muscle-specific gene expression. , 1998, Journal of molecular biology.

[18]  D. Haussler,et al.  Ultraconserved Elements in the Human Genome , 2004, Science.

[19]  Shane T. Jensen,et al.  Cis-regulatory modules in the mammalian liver: composition depends on strength of Foxa2 consensus site , 2008, Nucleic acids research.

[20]  M. Blanchette,et al.  Discovery of regulatory elements by a computational method for phylogenetic footprinting. , 2002, Genome research.

[21]  S. Hannenhalli,et al.  Position and distance specificity are important determinants of cis-regulatory motifs in addition to evolutionary conservation , 2007, Nucleic acids research.

[22]  R. Durbin,et al.  Comparative analysis of noncoding regions of 77 orthologous mouse and human gene pairs. , 1999, Genome research.

[23]  G. Tuteja,et al.  Cell-Specific Determinants of Peroxisome Proliferator-Activated Receptor γ Function in Adipocytes and Macrophages , 2010, Molecular and Cellular Biology.

[24]  Uwe Ohler,et al.  Phylogenetic simulation of promoter evolution: estimation and modeling of binding site turnover events and assessment of their impact on alignment tools , 2007, Genome Biology.

[25]  Hans Ellegren,et al.  Evidence for turnover of functional noncoding DNA in mammalian genome evolution. , 2004, Genomics.

[26]  Brad T. Sherman,et al.  Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources , 2008, Nature Protocols.

[27]  Geetu Tuteja,et al.  Foxa1 and Foxa2 regulate bile duct development in mice. , 2009, The Journal of clinical investigation.

[28]  F. Jacob,et al.  Evolution and tinkering. , 1977, Science.

[29]  A. Clark,et al.  Evolution of transcription factor binding sites in Mammalian gene regulatory regions: conservation and turnover. , 2002, Molecular biology and evolution.

[30]  C. Lawrence,et al.  Human-mouse genome comparisons to locate regulatory sites , 2000, Nature Genetics.

[31]  M. Gerstein,et al.  Of mice and men: phylogenetic footprinting aids the discovery of regulatory elements , 2003, Journal of biology.

[32]  M. Lazar,et al.  New developments in adipogenesis , 2009, Trends in Endocrinology & Metabolism.

[33]  Paul Schliekelman,et al.  Statistical Methods in Bioinformatics: An Introduction , 2001 .

[34]  Gregory R. Grant,et al.  Statistical Methods in Bioinformatics , 2001 .

[35]  Matthew W. Hahn,et al.  The evolution of transcriptional regulation in eukaryotes. , 2003, Molecular biology and evolution.

[36]  K. Kaestner,et al.  Foxa2 integrates the transcriptional response of the hepatocyte to fasting. , 2005, Cell metabolism.

[37]  M. Kreitman,et al.  Functional Evolution of a cis-Regulatory Module , 2005, PLoS biology.

[38]  R. Tibshirani,et al.  Significance analysis of microarrays applied to the ionizing radiation response , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[39]  M. Wabitsch,et al.  Characterization of a human preadipocyte cell strain with high capacity for adipose differentiation , 2001, International Journal of Obesity.

[40]  William Stafford Noble,et al.  Quantifying similarity between motifs , 2007, Genome Biology.

[41]  J. Garcia-Fernández,et al.  Implications of duplicated cis-regulatory elements in the evolution of metazoans: the DDI model or how simplicity begets novelty. , 2009, Briefings in functional genomics & proteomics.

[42]  Thomas Zeng,et al.  Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing , 2008, Nucleic acids research.

[43]  N. Patel,et al.  Evidence for stabilizing selection in a eukaryotic enhancer element , 2000, Nature.

[44]  Eric S. Lander,et al.  Comparative Epigenomic Analysis of Murine and Human Adipogenesis , 2010, Cell.

[45]  G. Tuteja,et al.  Extracting transcription factor targets from ChIP-Seq data , 2009, Nucleic acids research.

[46]  S. Mandrup,et al.  PPARγ in adipocyte differentiation and metabolism – Novel insights from genome‐wide studies , 2010, FEBS letters.

[47]  D. Haussler,et al.  Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes. , 2005, Genome research.

[48]  D. Tautz Evolution of transcriptional regulation. , 2000, Current opinion in genetics & development.

[49]  J. Schug,et al.  Re-expression of GATA2 Cooperates with Peroxisome Proliferator-activated Receptor-γ Depletion to Revert the Adipocyte Phenotype* , 2009, Journal of Biological Chemistry.

[50]  G. K. Davis,et al.  Phenotypic robustness conferred by apparently redundant transcriptional enhancers , 2010, Nature.

[51]  Jun S. Liu,et al.  An algorithm for finding protein–DNA binding sites with applications to chromatin-immunoprecipitation microarray experiments , 2002, Nature Biotechnology.

[52]  D. Haussler,et al.  Evolution's cauldron: Duplication, deletion, and rearrangement in the mouse and human genomes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[53]  K. Kaestner,et al.  Dynamic regulation of Pdx1 enhancers by Foxa1 and Foxa2 is essential for pancreas development. , 2008, Genes & development.