Corpus
Pmc
Racine
Corpus
Checkpoint
Pmc
Racine
Pmc
Exploration
Accueil
Racine
Racine

Serveur d'exploration MERS

Attention, ce site est en cours de développement !
Attention, site généré par des moyens informatiques à partir de corpus bruts.
Les informations ne sont donc pas validées.

Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway

Identifieur interne : 000289 ( Pmc/Corpus ); précédent : 000288; suivant : 000290

Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway

Auteurs : Katherine Gurdziel ; Kyle R. Vogt ; Gary Schneider ; Neil Richards ; Deborah L. Gumucio

Source :

RBID : PMC:4765071

Abstract

Background

The Hedgehog (Hh) signaling pathway, acting through three homologous transcription factors (GLI1, GLI2, GLI3) in vertebrates, plays multiple roles in embryonic organ development and adult tissue homeostasis. At the level of the genome, GLI factors bind to specific motifs in enhancers, some of which are hundreds of kilobases removed from the gene promoter. These enhancers integrate the Hh signal in a context-specific manner to control the spatiotemporal pattern of target gene expression. Importantly, a number of genes that encode Hh pathway molecules are themselves targets of Hh signaling, allowing pathway regulation by an intricate balance of feed-back activation and inhibition. However, surprisingly few of the critical enhancer elements that control these pathway target genes have been identified despite the fact that such elements are central determinants of Hh signaling activity. Recently, ChIP studies have been carried out in multiple tissue contexts using mouse models carrying FLAG-tagged GLI proteins (GLIFLAG). Using these datasets, we tested whether a meta-analysis of GLI binding sites, coupled with a machine learning approach, could reveal genomic features that could be used to empirically identify Hh-regulated enhancers linked to loci of the Hh signaling pathway.

Results

A meta-analysis of four existing GLIFLAG datasets revealed a library of GLI binding motifs that was substantially more restricted than the potential sites predicted by previous in vitro binding studies. A machine learning method (kmer-SVM) was then applied to these datasets and enriched k-mers were identified that, when applied to the mouse genome, predicted as many as 37,000 potential Hh enhancers. For functional analysis, we selected nine regions which were annotated to putative Hh pathway molecules and found that seven exhibited GLI-dependent activity, indicating that they are directly regulated by Hh signaling (78 % success rate).

Conclusions

The results suggest that Hh enhancer regions share common sequence features. The kmer-SVM machine learning approach identifies those features and can successfully predict functional Hh regulatory regions in genomic DNA surrounding Hh pathway molecules and likely, other Hh targets. Additionally, the library of enriched GLI binding motifs that we have identified may allow improved identification of functional GLI binding sites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-016-0106-0) contains supplementary material, which is available to authorized users.


Url:
DOI: 10.1186/s12861-016-0106-0
PubMed: 26912062
PubMed Central: 4765071

Links to Exploration step

PMC:4765071

Le document en format XML

<record>
<TEI>
<teiHeader>
<fileDesc>
<titleStmt>
<title xml:lang="en">Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway</title>
<author>
<name sortKey="Gurdziel, Katherine" sort="Gurdziel, Katherine" uniqKey="Gurdziel K" first="Katherine" last="Gurdziel">Katherine Gurdziel</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="Aff2">Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Vogt, Kyle R" sort="Vogt, Kyle R" uniqKey="Vogt K" first="Kyle R." last="Vogt">Kyle R. Vogt</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Schneider, Gary" sort="Schneider, Gary" uniqKey="Schneider G" first="Gary" last="Schneider">Gary Schneider</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Richards, Neil" sort="Richards, Neil" uniqKey="Richards N" first="Neil" last="Richards">Neil Richards</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Gumucio, Deborah L" sort="Gumucio, Deborah L" uniqKey="Gumucio D" first="Deborah L." last="Gumucio">Deborah L. Gumucio</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt>
<idno type="wicri:source">PMC</idno>
<idno type="pmid">26912062</idno>
<idno type="pmc">4765071</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4765071</idno>
<idno type="RBID">PMC:4765071</idno>
<idno type="doi">10.1186/s12861-016-0106-0</idno>
<date when="2016">2016</date>
<idno type="wicri:Area/Pmc/Corpus">000289</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000289</idno>
</publicationStmt>
<sourceDesc>
<biblStruct>
<analytic>
<title xml:lang="en" level="a" type="main">Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway</title>
<author>
<name sortKey="Gurdziel, Katherine" sort="Gurdziel, Katherine" uniqKey="Gurdziel K" first="Katherine" last="Gurdziel">Katherine Gurdziel</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
<affiliation>
<nlm:aff id="Aff2">Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Vogt, Kyle R" sort="Vogt, Kyle R" uniqKey="Vogt K" first="Kyle R." last="Vogt">Kyle R. Vogt</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Schneider, Gary" sort="Schneider, Gary" uniqKey="Schneider G" first="Gary" last="Schneider">Gary Schneider</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Richards, Neil" sort="Richards, Neil" uniqKey="Richards N" first="Neil" last="Richards">Neil Richards</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
<author>
<name sortKey="Gumucio, Deborah L" sort="Gumucio, Deborah L" uniqKey="Gumucio D" first="Deborah L." last="Gumucio">Deborah L. Gumucio</name>
<affiliation>
<nlm:aff id="Aff1">Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</nlm:aff>
</affiliation>
</author>
</analytic>
<series>
<title level="j">BMC Developmental Biology</title>
<idno type="eISSN">1471-213X</idno>
<imprint>
<date when="2016">2016</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc>
<textClass></textClass>
</profileDesc>
</teiHeader>
<front>
<div type="abstract" xml:lang="en">
<sec>
<title>Background</title>
<p>The Hedgehog (Hh) signaling pathway, acting through three homologous transcription factors (GLI1, GLI2, GLI3) in vertebrates, plays multiple roles in embryonic organ development and adult tissue homeostasis. At the level of the genome, GLI factors bind to specific motifs in enhancers, some of which are hundreds of kilobases removed from the gene promoter. These enhancers integrate the Hh signal in a context-specific manner to control the spatiotemporal pattern of target gene expression. Importantly, a number of genes that encode Hh pathway molecules are themselves targets of Hh signaling, allowing pathway regulation by an intricate balance of feed-back activation and inhibition. However, surprisingly few of the critical enhancer elements that control these pathway target genes have been identified despite the fact that such elements are central determinants of Hh signaling activity. Recently, ChIP studies have been carried out in multiple tissue contexts using mouse models carrying FLAG-tagged GLI proteins (GLI
<sup>FLAG</sup>
). Using these datasets, we tested whether a meta-analysis of GLI binding sites, coupled with a machine learning approach, could reveal genomic features that could be used to empirically identify Hh-regulated enhancers linked to loci of the Hh signaling pathway.</p>
</sec>
<sec>
<title>Results</title>
<p>A meta-analysis of four existing GLI
<sup>FLAG</sup>
datasets revealed a library of GLI binding motifs that was substantially more restricted than the potential sites predicted by previous in vitro binding studies. A machine learning method (kmer-SVM) was then applied to these datasets and enriched k-mers were identified that, when applied to the mouse genome, predicted as many as 37,000 potential Hh enhancers. For functional analysis, we selected nine regions which were annotated to putative Hh pathway molecules and found that seven exhibited GLI-dependent activity, indicating that they are directly regulated by Hh signaling (78 % success rate).</p>
</sec>
<sec>
<title>Conclusions</title>
<p>The results suggest that Hh enhancer regions share common sequence features. The kmer-SVM machine learning approach identifies those features and can successfully predict functional Hh regulatory regions in genomic DNA surrounding Hh pathway molecules and likely, other Hh targets. Additionally, the library of enriched GLI binding motifs that we have identified may allow improved identification of functional GLI binding sites.</p>
</sec>
<sec>
<title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/s12861-016-0106-0) contains supplementary material, which is available to authorized users.</p>
</sec>
</div>
</front>
<back>
<div1 type="bibliography">
<listBibl>
<biblStruct>
<analytic>
<author>
<name sortKey="Briscoe, J" uniqKey="Briscoe J">J Briscoe</name>
</author>
<author>
<name sortKey="Therond, Pp" uniqKey="Therond P">PP Thérond</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Jiang, J" uniqKey="Jiang J">J Jiang</name>
</author>
<author>
<name sortKey="Hui, Cc" uniqKey="Hui C">CC Hui</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Mcmahon, Ap" uniqKey="Mcmahon A">AP McMahon</name>
</author>
<author>
<name sortKey="Ingham, Pw" uniqKey="Ingham P">PW Ingham</name>
</author>
<author>
<name sortKey="Tabin, Cj" uniqKey="Tabin C">CJ Tabin</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Petrova, R" uniqKey="Petrova R">R Petrova</name>
</author>
<author>
<name sortKey="Joyner, Al" uniqKey="Joyner A">AL Joyner</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Varjosalo, M" uniqKey="Varjosalo M">M Varjosalo</name>
</author>
<author>
<name sortKey="Taipale, J" uniqKey="Taipale J">J Taipale</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Creanga, A" uniqKey="Creanga A">A Creanga</name>
</author>
<author>
<name sortKey="Glenn, Td" uniqKey="Glenn T">TD Glenn</name>
</author>
<author>
<name sortKey="Mann, Rk" uniqKey="Mann R">RK Mann</name>
</author>
<author>
<name sortKey="Saunders, Am" uniqKey="Saunders A">AM Saunders</name>
</author>
<author>
<name sortKey="Talbot, Ws" uniqKey="Talbot W">WS Talbot</name>
</author>
<author>
<name sortKey="Beachy, Pa" uniqKey="Beachy P">PA Beachy</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Jakobs, P" uniqKey="Jakobs P">P Jakobs</name>
</author>
<author>
<name sortKey="Exner, S" uniqKey="Exner S">S Exner</name>
</author>
<author>
<name sortKey="Schurmann, S" uniqKey="Schurmann S">S Schurmann</name>
</author>
<author>
<name sortKey="Pickhinke, U" uniqKey="Pickhinke U">U Pickhinke</name>
</author>
<author>
<name sortKey="Bandari, S" uniqKey="Bandari S">S Bandari</name>
</author>
<author>
<name sortKey="Ortmann, C" uniqKey="Ortmann C">C Ortmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Burke, R" uniqKey="Burke R">R Burke</name>
</author>
<author>
<name sortKey="Nellen, D" uniqKey="Nellen D">D Nellen</name>
</author>
<author>
<name sortKey="Bellotto, M" uniqKey="Bellotto M">M Bellotto</name>
</author>
<author>
<name sortKey="Hafen, E" uniqKey="Hafen E">E Hafen</name>
</author>
<author>
<name sortKey="Senti, Ka" uniqKey="Senti K">KA Senti</name>
</author>
<author>
<name sortKey="Dickson, Bj" uniqKey="Dickson B">BJ Dickson</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Jeong, J" uniqKey="Jeong J">J Jeong</name>
</author>
<author>
<name sortKey="Mcmahon, Ap" uniqKey="Mcmahon A">AP McMahon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Izzi, L" uniqKey="Izzi L">L Izzi</name>
</author>
<author>
<name sortKey="Levesque, M" uniqKey="Levesque M">M Levesque</name>
</author>
<author>
<name sortKey="Morin, S" uniqKey="Morin S">S Morin</name>
</author>
<author>
<name sortKey="Laniel, D" uniqKey="Laniel D">D Laniel</name>
</author>
<author>
<name sortKey="Wilkes, Bc" uniqKey="Wilkes B">BC Wilkes</name>
</author>
<author>
<name sortKey="Mille, F" uniqKey="Mille F">F Mille</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Allen, Bl" uniqKey="Allen B">BL Allen</name>
</author>
<author>
<name sortKey="Song, Jy" uniqKey="Song J">JY Song</name>
</author>
<author>
<name sortKey="Izzi, L" uniqKey="Izzi L">L Izzi</name>
</author>
<author>
<name sortKey="Althaus, Iw" uniqKey="Althaus I">IW Althaus</name>
</author>
<author>
<name sortKey="Kang, Js" uniqKey="Kang J">JS Kang</name>
</author>
<author>
<name sortKey="Charron, F" uniqKey="Charron F">F Charron</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Chuang, Pt" uniqKey="Chuang P">PT Chuang</name>
</author>
<author>
<name sortKey="Mcmahon, Ap" uniqKey="Mcmahon A">AP McMahon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Kim, J" uniqKey="Kim J">J Kim</name>
</author>
<author>
<name sortKey="Kato, M" uniqKey="Kato M">M Kato</name>
</author>
<author>
<name sortKey="Beachy, Pa" uniqKey="Beachy P">PA Beachy</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Goetz, Sc" uniqKey="Goetz S">SC Goetz</name>
</author>
<author>
<name sortKey="Anderson, Kv" uniqKey="Anderson K">KV Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Dai, P" uniqKey="Dai P">P Dai</name>
</author>
<author>
<name sortKey="Akimaru, H" uniqKey="Akimaru H">H Akimaru</name>
</author>
<author>
<name sortKey="Tanaka, Y" uniqKey="Tanaka Y">Y Tanaka</name>
</author>
<author>
<name sortKey="Maekawa, T" uniqKey="Maekawa T">T Maekawa</name>
</author>
<author>
<name sortKey="Nakafuku, M" uniqKey="Nakafuku M">M Nakafuku</name>
</author>
<author>
<name sortKey="Ishii, S" uniqKey="Ishii S">S Ishii</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Shin, Sh" uniqKey="Shin S">SH Shin</name>
</author>
<author>
<name sortKey="Kogerman, P" uniqKey="Kogerman P">P Kogerman</name>
</author>
<author>
<name sortKey="Lindstrom, E" uniqKey="Lindstrom E">E Lindstrom</name>
</author>
<author>
<name sortKey="Toftgard, R" uniqKey="Toftgard R">R Toftgard</name>
</author>
<author>
<name sortKey="Biesecker, Lg" uniqKey="Biesecker L">LG Biesecker</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Sasaki, H" uniqKey="Sasaki H">H Sasaki</name>
</author>
<author>
<name sortKey="Nishizaki, Y" uniqKey="Nishizaki Y">Y Nishizaki</name>
</author>
<author>
<name sortKey="Hui, C" uniqKey="Hui C">C Hui</name>
</author>
<author>
<name sortKey="Nakafuku, M" uniqKey="Nakafuku M">M Nakafuku</name>
</author>
<author>
<name sortKey="Kondoh, H" uniqKey="Kondoh H">H Kondoh</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Beachy, Pa" uniqKey="Beachy P">PA Beachy</name>
</author>
<author>
<name sortKey="Hymowitz, Sg" uniqKey="Hymowitz S">SG Hymowitz</name>
</author>
<author>
<name sortKey="Lazarus, Ra" uniqKey="Lazarus R">RA Lazarus</name>
</author>
<author>
<name sortKey="Leahy, Dj" uniqKey="Leahy D">DJ Leahy</name>
</author>
<author>
<name sortKey="Siebold, C" uniqKey="Siebold C">C Siebold</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Liem, Kf" uniqKey="Liem K">KF Liem</name>
</author>
<author>
<name sortKey="He, M" uniqKey="He M">M He</name>
</author>
<author>
<name sortKey="Ocbina, Pj" uniqKey="Ocbina P">PJ Ocbina</name>
</author>
<author>
<name sortKey="Anderson, Kv" uniqKey="Anderson K">KV Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Aglyamova, Gv" uniqKey="Aglyamova G">GV Aglyamova</name>
</author>
<author>
<name sortKey="Agarwala, S" uniqKey="Agarwala S">S Agarwala</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Tenzen, T" uniqKey="Tenzen T">T Tenzen</name>
</author>
<author>
<name sortKey="Allen, Bl" uniqKey="Allen B">BL Allen</name>
</author>
<author>
<name sortKey="Cole, F" uniqKey="Cole F">F Cole</name>
</author>
<author>
<name sortKey="Kang, Js" uniqKey="Kang J">JS Kang</name>
</author>
<author>
<name sortKey="Krauss, Rs" uniqKey="Krauss R">RS Krauss</name>
</author>
<author>
<name sortKey="Mcmahon, Ap" uniqKey="Mcmahon A">AP McMahon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Martinelli, Dc" uniqKey="Martinelli D">DC Martinelli</name>
</author>
<author>
<name sortKey="Fan, Cm" uniqKey="Fan C">CM Fan</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Rahnama, F" uniqKey="Rahnama F">F Rahnama</name>
</author>
<author>
<name sortKey="Toftgard, R" uniqKey="Toftgard R">R Toftgard</name>
</author>
<author>
<name sortKey="Zaphiropoulos, Pg" uniqKey="Zaphiropoulos P">PG Zaphiropoulos</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Agren, M" uniqKey="Agren M">M Agren</name>
</author>
<author>
<name sortKey="Kogerman, P" uniqKey="Kogerman P">P Kogerman</name>
</author>
<author>
<name sortKey="Kleman, Mi" uniqKey="Kleman M">MI Kleman</name>
</author>
<author>
<name sortKey="Wessling, M" uniqKey="Wessling M">M Wessling</name>
</author>
<author>
<name sortKey="Toftgard, R" uniqKey="Toftgard R">R Toftgard</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Peterson, Ka" uniqKey="Peterson K">KA Peterson</name>
</author>
<author>
<name sortKey="Nishi, Y" uniqKey="Nishi Y">Y Nishi</name>
</author>
<author>
<name sortKey="Ma, W" uniqKey="Ma W">W Ma</name>
</author>
<author>
<name sortKey="Vedenko, A" uniqKey="Vedenko A">A Vedenko</name>
</author>
<author>
<name sortKey="Shokri, L" uniqKey="Shokri L">L Shokri</name>
</author>
<author>
<name sortKey="Zhang, X" uniqKey="Zhang X">X Zhang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Vokes, Sa" uniqKey="Vokes S">SA Vokes</name>
</author>
<author>
<name sortKey="Ji, H" uniqKey="Ji H">H Ji</name>
</author>
<author>
<name sortKey="Mccuine, S" uniqKey="Mccuine S">S McCuine</name>
</author>
<author>
<name sortKey="Tenzen, T" uniqKey="Tenzen T">T Tenzen</name>
</author>
<author>
<name sortKey="Giles, S" uniqKey="Giles S">S Giles</name>
</author>
<author>
<name sortKey="Zhong, S" uniqKey="Zhong S">S Zhong</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Vokes, Sa" uniqKey="Vokes S">SA Vokes</name>
</author>
<author>
<name sortKey="Ji, H" uniqKey="Ji H">H Ji</name>
</author>
<author>
<name sortKey="Wong, Wh" uniqKey="Wong W">WH Wong</name>
</author>
<author>
<name sortKey="Mcmahon, Ap" uniqKey="Mcmahon A">AP McMahon</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Holtz, Am" uniqKey="Holtz A">AM Holtz</name>
</author>
<author>
<name sortKey="Peterson, Ka" uniqKey="Peterson K">KA Peterson</name>
</author>
<author>
<name sortKey="Nishi, Y" uniqKey="Nishi Y">Y Nishi</name>
</author>
<author>
<name sortKey="Morin, S" uniqKey="Morin S">S Morin</name>
</author>
<author>
<name sortKey="Song, Jy" uniqKey="Song J">JY Song</name>
</author>
<author>
<name sortKey="Charron, F" uniqKey="Charron F">F Charron</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lee, Ey" uniqKey="Lee E">EY Lee</name>
</author>
<author>
<name sortKey="Ji, H" uniqKey="Ji H">H Ji</name>
</author>
<author>
<name sortKey="Ouyang, Z" uniqKey="Ouyang Z">Z Ouyang</name>
</author>
<author>
<name sortKey="Zhou, B" uniqKey="Zhou B">B Zhou</name>
</author>
<author>
<name sortKey="Ma, W" uniqKey="Ma W">W Ma</name>
</author>
<author>
<name sortKey="Vokes, Sa" uniqKey="Vokes S">SA Vokes</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Fletez Brant, C" uniqKey="Fletez Brant C">C Fletez-Brant</name>
</author>
<author>
<name sortKey="Lee, D" uniqKey="Lee D">D Lee</name>
</author>
<author>
<name sortKey="Mccallion, As" uniqKey="Mccallion A">AS McCallion</name>
</author>
<author>
<name sortKey="Beer, Ma" uniqKey="Beer M">MA Beer</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Lee, D" uniqKey="Lee D">D Lee</name>
</author>
<author>
<name sortKey="Karchin, R" uniqKey="Karchin R">R Karchin</name>
</author>
<author>
<name sortKey="Beer, Ma" uniqKey="Beer M">MA Beer</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Spitz, F" uniqKey="Spitz F">F Spitz</name>
</author>
<author>
<name sortKey="Furlong, Ee" uniqKey="Furlong E">EE Furlong</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Slattery, M" uniqKey="Slattery M">M Slattery</name>
</author>
<author>
<name sortKey="Zhou, T" uniqKey="Zhou T">T Zhou</name>
</author>
<author>
<name sortKey="Yang, L" uniqKey="Yang L">L Yang</name>
</author>
<author>
<name sortKey="Dantas Machado, Ac" uniqKey="Dantas Machado A">AC Dantas Machado</name>
</author>
<author>
<name sortKey="Gordan, R" uniqKey="Gordan R">R Gordan</name>
</author>
<author>
<name sortKey="Rohs, R" uniqKey="Rohs R">R Rohs</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Hallikas, O" uniqKey="Hallikas O">O Hallikas</name>
</author>
<author>
<name sortKey="Palin, K" uniqKey="Palin K">K Palin</name>
</author>
<author>
<name sortKey="Sinjushina, N" uniqKey="Sinjushina N">N Sinjushina</name>
</author>
<author>
<name sortKey="Rautiainen, R" uniqKey="Rautiainen R">R Rautiainen</name>
</author>
<author>
<name sortKey="Partanen, J" uniqKey="Partanen J">J Partanen</name>
</author>
<author>
<name sortKey="Ukkonen, E" uniqKey="Ukkonen E">E Ukkonen</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Heinz, S" uniqKey="Heinz S">S Heinz</name>
</author>
<author>
<name sortKey="Benner, C" uniqKey="Benner C">C Benner</name>
</author>
<author>
<name sortKey="Spann, N" uniqKey="Spann N">N Spann</name>
</author>
<author>
<name sortKey="Bertolino, E" uniqKey="Bertolino E">E Bertolino</name>
</author>
<author>
<name sortKey="Lin, Yc" uniqKey="Lin Y">YC Lin</name>
</author>
<author>
<name sortKey="Laslo, P" uniqKey="Laslo P">P Laslo</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bailey, Tl" uniqKey="Bailey T">TL Bailey</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gupta, S" uniqKey="Gupta S">S Gupta</name>
</author>
<author>
<name sortKey="Stamatoyannopoulos, Ja" uniqKey="Stamatoyannopoulos J">JA Stamatoyannopoulos</name>
</author>
<author>
<name sortKey="Bailey, Tl" uniqKey="Bailey T">TL Bailey</name>
</author>
<author>
<name sortKey="Noble, Ws" uniqKey="Noble W">WS Noble</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Crooks, Ge" uniqKey="Crooks G">GE Crooks</name>
</author>
<author>
<name sortKey="Hon, G" uniqKey="Hon G">G Hon</name>
</author>
<author>
<name sortKey="Chandonia, Jm" uniqKey="Chandonia J">JM Chandonia</name>
</author>
<author>
<name sortKey="Brenner, Se" uniqKey="Brenner S">SE Brenner</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Winklmayr, M" uniqKey="Winklmayr M">M Winklmayr</name>
</author>
<author>
<name sortKey="Schmid, C" uniqKey="Schmid C">C Schmid</name>
</author>
<author>
<name sortKey="Laner Plamberger, S" uniqKey="Laner Plamberger S">S Laner-Plamberger</name>
</author>
<author>
<name sortKey="Kaser, A" uniqKey="Kaser A">A Kaser</name>
</author>
<author>
<name sortKey="Aberger, F" uniqKey="Aberger F">F Aberger</name>
</author>
<author>
<name sortKey="Eichberger, T" uniqKey="Eichberger T">T Eichberger</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Machanick, P" uniqKey="Machanick P">P Machanick</name>
</author>
<author>
<name sortKey="Bailey, Tl" uniqKey="Bailey T">TL Bailey</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Bailey, Tl" uniqKey="Bailey T">TL Bailey</name>
</author>
<author>
<name sortKey="Machanick, P" uniqKey="Machanick P">P Machanick</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Creyghton, Mp" uniqKey="Creyghton M">MP Creyghton</name>
</author>
<author>
<name sortKey="Cheng, Aw" uniqKey="Cheng A">AW Cheng</name>
</author>
<author>
<name sortKey="Welstead, Gg" uniqKey="Welstead G">GG Welstead</name>
</author>
<author>
<name sortKey="Kooistra, T" uniqKey="Kooistra T">T Kooistra</name>
</author>
<author>
<name sortKey="Carey, Bw" uniqKey="Carey B">BW Carey</name>
</author>
<author>
<name sortKey="Steine, Ej" uniqKey="Steine E">EJ Steine</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Heintzman, Nd" uniqKey="Heintzman N">ND Heintzman</name>
</author>
<author>
<name sortKey="Stuart, Rk" uniqKey="Stuart R">RK Stuart</name>
</author>
<author>
<name sortKey="Hon, G" uniqKey="Hon G">G Hon</name>
</author>
<author>
<name sortKey="Fu, Y" uniqKey="Fu Y">Y Fu</name>
</author>
<author>
<name sortKey="Ching, Cw" uniqKey="Ching C">CW Ching</name>
</author>
<author>
<name sortKey="Hawkins, Rd" uniqKey="Hawkins R">RD Hawkins</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Visel, A" uniqKey="Visel A">A Visel</name>
</author>
<author>
<name sortKey="Thaller, C" uniqKey="Thaller C">C Thaller</name>
</author>
<author>
<name sortKey="Eichele, G" uniqKey="Eichele G">G Eichele</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Diez Roux, G" uniqKey="Diez Roux G">G Diez-Roux</name>
</author>
<author>
<name sortKey="Banfi, S" uniqKey="Banfi S">S Banfi</name>
</author>
<author>
<name sortKey="Sultan, M" uniqKey="Sultan M">M Sultan</name>
</author>
<author>
<name sortKey="Geffers, L" uniqKey="Geffers L">L Geffers</name>
</author>
<author>
<name sortKey="Anand, S" uniqKey="Anand S">S Anand</name>
</author>
<author>
<name sortKey="Rozado, D" uniqKey="Rozado D">D Rozado</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Mclean, Cy" uniqKey="Mclean C">CY McLean</name>
</author>
<author>
<name sortKey="Bristor, D" uniqKey="Bristor D">D Bristor</name>
</author>
<author>
<name sortKey="Hiller, M" uniqKey="Hiller M">M Hiller</name>
</author>
<author>
<name sortKey="Clarke, Sl" uniqKey="Clarke S">SL Clarke</name>
</author>
<author>
<name sortKey="Schaar, Bt" uniqKey="Schaar B">BT Schaar</name>
</author>
<author>
<name sortKey="Lowe, Cb" uniqKey="Lowe C">CB Lowe</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Gerber, An" uniqKey="Gerber A">AN Gerber</name>
</author>
<author>
<name sortKey="Wilson, Cw" uniqKey="Wilson C">CW Wilson</name>
</author>
<author>
<name sortKey="Li, Yj" uniqKey="Li Y">YJ Li</name>
</author>
<author>
<name sortKey="Chuang, Pt" uniqKey="Chuang P">PT Chuang</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Zhao, M" uniqKey="Zhao M">M Zhao</name>
</author>
<author>
<name sortKey="Qiao, M" uniqKey="Qiao M">M Qiao</name>
</author>
<author>
<name sortKey="Harris, Se" uniqKey="Harris S">SE Harris</name>
</author>
<author>
<name sortKey="Chen, D" uniqKey="Chen D">D Chen</name>
</author>
<author>
<name sortKey="Oyajobi, Bo" uniqKey="Oyajobi B">BO Oyajobi</name>
</author>
<author>
<name sortKey="Mundy, Gr" uniqKey="Mundy G">GR Mundy</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Levi, G" uniqKey="Levi G">G Levi</name>
</author>
<author>
<name sortKey="Topilko, P" uniqKey="Topilko P">P Topilko</name>
</author>
<author>
<name sortKey="Schneider Maunoury, S" uniqKey="Schneider Maunoury S">S Schneider-Maunoury</name>
</author>
<author>
<name sortKey="Lasagna, M" uniqKey="Lasagna M">M Lasagna</name>
</author>
<author>
<name sortKey="Mantero, S" uniqKey="Mantero S">S Mantero</name>
</author>
<author>
<name sortKey="Cancedda, R" uniqKey="Cancedda R">R Cancedda</name>
</author>
</analytic>
</biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Nonchev, S" uniqKey="Nonchev S">S Nonchev</name>
</author>
<author>
<name sortKey="Maconochie, M" uniqKey="Maconochie M">M Maconochie</name>
</author>
<author>
<name sortKey="Vesque, C" uniqKey="Vesque C">C Vesque</name>
</author>
<author>
<name sortKey="Aparicio, S" uniqKey="Aparicio S">S Aparicio</name>
</author>
<author>
<name sortKey="Ariza Mcnaughton, L" uniqKey="Ariza Mcnaughton L">L Ariza-McNaughton</name>
</author>
<author>
<name sortKey="Manzanares, M" uniqKey="Manzanares M">M Manzanares</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct>
<analytic>
<author>
<name sortKey="Santos, N" uniqKey="Santos N">N Santos</name>
</author>
<author>
<name sortKey="Reiter, Jf" uniqKey="Reiter J">JF Reiter</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article">
<pmc-dir>properties open_access</pmc-dir>
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">BMC Dev Biol</journal-id>
<journal-id journal-id-type="iso-abbrev">BMC Dev. Biol</journal-id>
<journal-title-group>
<journal-title>BMC Developmental Biology</journal-title>
</journal-title-group>
<issn pub-type="epub">1471-213X</issn>
<publisher>
<publisher-name>BioMed Central</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="pmid">26912062</article-id>
<article-id pub-id-type="pmc">4765071</article-id>
<article-id pub-id-type="publisher-id">106</article-id>
<article-id pub-id-type="doi">10.1186/s12861-016-0106-0</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Gurdziel</surname>
<given-names>Katherine</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
<xref ref-type="aff" rid="Aff2"></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Vogt</surname>
<given-names>Kyle R.</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Schneider</surname>
<given-names>Gary</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Richards</surname>
<given-names>Neil</given-names>
</name>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Gumucio</surname>
<given-names>Deborah L.</given-names>
</name>
<address>
<phone>734-647-0172</phone>
<email>dgumucio@umich.edu</email>
</address>
<xref ref-type="aff" rid="Aff1"></xref>
</contrib>
<aff id="Aff1">
<label></label>
Department of Cell and Developmental Biology, University of Michigan, Ann Arbor, MI 48109 USA</aff>
<aff id="Aff2">
<label></label>
Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109 USA</aff>
</contrib-group>
<pub-date pub-type="epub">
<day>24</day>
<month>2</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="pmc-release">
<day>24</day>
<month>2</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="collection">
<year>2016</year>
</pub-date>
<volume>16</volume>
<elocation-id>4</elocation-id>
<history>
<date date-type="received">
<day>23</day>
<month>6</month>
<year>2015</year>
</date>
<date date-type="accepted">
<day>16</day>
<month>2</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-statement>© Gurdziel et al. 2016</copyright-statement>
<license license-type="OpenAccess">
<license-p>
<bold>Open Access</bold>
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (
<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/publicdomain/zero/1.0/">http://creativecommons.org/publicdomain/zero/1.0/</ext-link>
) applies to the data made available in this article, unless otherwise stated.</license-p>
</license>
</permissions>
<abstract id="Abs1">
<sec>
<title>Background</title>
<p>The Hedgehog (Hh) signaling pathway, acting through three homologous transcription factors (GLI1, GLI2, GLI3) in vertebrates, plays multiple roles in embryonic organ development and adult tissue homeostasis. At the level of the genome, GLI factors bind to specific motifs in enhancers, some of which are hundreds of kilobases removed from the gene promoter. These enhancers integrate the Hh signal in a context-specific manner to control the spatiotemporal pattern of target gene expression. Importantly, a number of genes that encode Hh pathway molecules are themselves targets of Hh signaling, allowing pathway regulation by an intricate balance of feed-back activation and inhibition. However, surprisingly few of the critical enhancer elements that control these pathway target genes have been identified despite the fact that such elements are central determinants of Hh signaling activity. Recently, ChIP studies have been carried out in multiple tissue contexts using mouse models carrying FLAG-tagged GLI proteins (GLI
<sup>FLAG</sup>
). Using these datasets, we tested whether a meta-analysis of GLI binding sites, coupled with a machine learning approach, could reveal genomic features that could be used to empirically identify Hh-regulated enhancers linked to loci of the Hh signaling pathway.</p>
</sec>
<sec>
<title>Results</title>
<p>A meta-analysis of four existing GLI
<sup>FLAG</sup>
datasets revealed a library of GLI binding motifs that was substantially more restricted than the potential sites predicted by previous in vitro binding studies. A machine learning method (kmer-SVM) was then applied to these datasets and enriched k-mers were identified that, when applied to the mouse genome, predicted as many as 37,000 potential Hh enhancers. For functional analysis, we selected nine regions which were annotated to putative Hh pathway molecules and found that seven exhibited GLI-dependent activity, indicating that they are directly regulated by Hh signaling (78 % success rate).</p>
</sec>
<sec>
<title>Conclusions</title>
<p>The results suggest that Hh enhancer regions share common sequence features. The kmer-SVM machine learning approach identifies those features and can successfully predict functional Hh regulatory regions in genomic DNA surrounding Hh pathway molecules and likely, other Hh targets. Additionally, the library of enriched GLI binding motifs that we have identified may allow improved identification of functional GLI binding sites.</p>
</sec>
<sec>
<title>Electronic supplementary material</title>
<p>The online version of this article (doi:10.1186/s12861-016-0106-0) contains supplementary material, which is available to authorized users.</p>
</sec>
</abstract>
<kwd-group xml:lang="en">
<title>Keywords</title>
<kwd>Hedgehog signaling</kwd>
<kwd>Enhancers</kwd>
<kwd>Machine learning</kwd>
<kwd>GLI</kwd>
</kwd-group>
<funding-group>
<award-group>
<funding-source>
<institution>University of Michigan NIH Training Program in Basic and Translational Digestive Science</institution>
</funding-source>
<award-id>T32 DK094775</award-id>
<principal-award-recipient>
<name>
<surname>Gurdziel</surname>
<given-names>Katherine</given-names>
</name>
</principal-award-recipient>
</award-group>
<award-group>
<funding-source>
<institution>NIDDK NIH HHS</institution>
</funding-source>
<award-id>P01 DK062041</award-id>
<principal-award-recipient>
<name>
<surname>Gumucio</surname>
<given-names>Deborah L.</given-names>
</name>
</principal-award-recipient>
</award-group>
</funding-group>
<custom-meta-group>
<custom-meta>
<meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2016</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="Sec1" sec-type="introduction">
<title>Background</title>
<p>The Hedgehog (Hh) signaling pathway is critical for embryonic organ development and adult tissue homeostasis across animal phyla [
<xref ref-type="bibr" rid="CR1">1</xref>
–
<xref ref-type="bibr" rid="CR4">4</xref>
]. In multiple tissue-specific settings, Hh signaling directs specific cell fate choices, controls tissue patterning and governs cell proliferation. In mammals, Hh signaling originates with any of three ligands (Sonic hedgehog (
<italic>Shh</italic>
), Indian hedgehog (
<italic>Ihh</italic>
) or Desert hedgehog (
<italic>Dhh</italic>
)) (for Review see [
<xref ref-type="bibr" rid="CR5">5</xref>
]). Release of the lipid modified HH ligand has been shown to be facilitated by SCUBE in conjunction with the transmembrane protein Dispatched [
<xref ref-type="bibr" rid="CR6">6</xref>
–
<xref ref-type="bibr" rid="CR9">9</xref>
]. Once discharged, HH ligands interact with the Patched (PTCH1 or PTCH2) receptor protein and with Hh-binding proteins BOC, CDO and GAS1 on target cells to relieve PTCH-dependent inhibition of the Smoothened (SMO) transmembrane protein [
<xref ref-type="bibr" rid="CR10">10</xref>
,
<xref ref-type="bibr" rid="CR11">11</xref>
]. HH ligands can also be sequestered by the Hedgehog-interacting protein (HHIP), which dampens signaling [
<xref ref-type="bibr" rid="CR12">12</xref>
].</p>
<p>Hh-mediated signal transduction culminates in the nucleus, with the binding of zinc-finger transcription factors (GLI1, GLI2, GLI3) to target gene sequences [
<xref ref-type="bibr" rid="CR13">13</xref>
,
<xref ref-type="bibr" rid="CR14">14</xref>
]. However, proteolytic processing determines whether the GLI proteins act as repressors or activators. GLI1, which is not processed, functions exclusively as a transcriptional activator and may act to amplify Hh signals [
<xref ref-type="bibr" rid="CR15">15</xref>
]. GLI2 and GLI3 can be converted to a repressor form in the absence of Hh ligand. In the presence of the Hh ligand, this processing is inhibited, allowing full-length GLI proteins to traffic to the nucleus and activate gene expression [
<xref ref-type="bibr" rid="CR15">15</xref>
–
<xref ref-type="bibr" rid="CR17">17</xref>
]. Processing of GLI proteins requires passage through the cilia [
<xref ref-type="bibr" rid="CR13">13</xref>
,
<xref ref-type="bibr" rid="CR18">18</xref>
]; the kinesin KIF7 helps to properly construct the cilium and is enriched at the cilium tip, along with GLI and SUFU (Suppressor of Fused) [
<xref ref-type="bibr" rid="CR19">19</xref>
].</p>
<p>The Hh signaling pathway is regulated by both positive and negative feedback. Indeed, a number of Hh pathway components, including
<italic>Boc</italic>
,
<italic>Cdo</italic>
,
<italic>Gas1</italic>
,
<italic>Gli1</italic>
,
<italic>Hhip</italic>
,
<italic>Ptch1</italic>
and
<italic>Ptch2</italic>
are thought to be direct transcriptional targets of Hh signaling in multiple tissue contexts [
<xref ref-type="bibr" rid="CR12">12</xref>
,
<xref ref-type="bibr" rid="CR15">15</xref>
,
<xref ref-type="bibr" rid="CR20">20</xref>
–
<xref ref-type="bibr" rid="CR29">29</xref>
]. Thus, an important aspect of Hh pathway self-regulation is integrated at the level of the enhancers that control response of the pathway target genes to local Hh signaling levels. However, despite the high functional conservation of this pathway, surprisingly little is known about the enhancer elements that control self-regulation in any organism.</p>
<p>One way to identify Hh target enhancers is to perform chromatin immunoprecipitation (ChIP). Genetically modified mouse models carrying inducible FLAG-tagged GLI proteins have allowed analysis of GLI binding sites in vivo in several different tissue contexts. Four in vivo GLI binding studies, including three ChIP-chip analyses [
<xref ref-type="bibr" rid="CR26">26</xref>
,
<xref ref-type="bibr" rid="CR27">27</xref>
,
<xref ref-type="bibr" rid="CR29">29</xref>
] and one ChIP-seq study [
<xref ref-type="bibr" rid="CR25">25</xref>
], have been carried out using these models. Interestingly, examination of all four datasets for common GLI binding sites that are annotated to Hh pathway molecules reveals only three such sites (in
<italic>Gli1</italic>
,
<italic>Ptch1</italic>
, and
<italic>Ptch2</italic>
loci [
<xref ref-type="bibr" rid="CR15">15</xref>
,
<xref ref-type="bibr" rid="CR24">24</xref>
,
<xref ref-type="bibr" rid="CR28">28</xref>
]) that are uniformly detectable. Several other established Hh pathway genes, including
<italic>Boc</italic>
,
<italic>Hhip</italic>
,
<italic>Gli2</italic>
, and
<italic>Hipk2</italic>
, appear to exhibit different GLI-bound genomic locations, depending on context, suggesting that each of these pathway components is regulated by multiple distinct genomic enhancers that have context-specific features.</p>
<p>Using ChIP studies on diverse tissues, it may be possible to eventually identify all of the multiple enhancers that control each target gene in every context. While a valuable goal, such analyses are currently expensive and time consuming and often technically challenging where the number of cells available for analysis is limiting, as in many developmental contexts. Importantly, computational methods can reveal sequence features that characterize enhancer activity. We therefore asked whether analysis of all existing GLI ChIP data could reveal common sequence features that might be used to empirically and globally predict functional enhancers
<italic>de novo</italic>
. A publicly available machine learning approach, kmer-SVM [
<xref ref-type="bibr" rid="CR30">30</xref>
], was used to predict novel Hh enhancer regions. This tool uses a support vector machine (SVM) to determine sequence features (k-mer frequencies) that identify positive genomic regulatory regions [
<xref ref-type="bibr" rid="CR31">31</xref>
]. SVMs are classifier algorithms that define a boundary between members of two different groups. Kmer-SVM calculates weights for sequence features that determine the effectiveness of that feature to distinguish between positive and negative regulatory regions. Once the features are determined, they can be used to identify novel enhancer regions not present in the original positive set. The strength of this approach is that it relies exclusively on short regions of DNA sequence (length 3–10 bp) which are in the size range of transcription factor binding sites (TFBS). Additionally, the organization of the k-mers within a sequence does not impact the score; this feature is consistent with the variable arrangement of TFBS in enhancers [
<xref ref-type="bibr" rid="CR32">32</xref>
,
<xref ref-type="bibr" rid="CR33">33</xref>
].</p>
<p>Using the kmer-SVM tool [
<xref ref-type="bibr" rid="CR30">30</xref>
], analysis of the four existing GLI binding datasets identified a set of k-mers that appeared to successfully predict potential GLI-regulated enhancers. Application of this set of k-mers to the mouse genome pinpointed over 37,000 potential enhancers. Several putative enhancers that were annotated to Hh pathway components were then tested for their ability to drive Hh-dependent activity in transfected cells. The functional significance of the GLI binding motifs (GBM) was also tested within each active enhancer by mutation. Of the nine predicted regulatory regions tested, seven (78 %) drove reporter expression in a GLI-dependent fashion. These findings substantially increase the number of functionally verified Hh enhancers found in Hh pathway molecules and validate the use of machine learning on ChIP data as a valuable tool to empirically predict likely Hh-dependent regulatory regions.</p>
</sec>
<sec id="Sec2" sec-type="results">
<title>Results and discussion</title>
<sec id="Sec3">
<title>Analysis of GLI
<sup>FLAG</sup>
datasets to identify likely in vivo GLI transcription factor binding motifs</title>
<p>A previous in vitro analysis of GLI transcription factor binding resulted in the identification of a set of likely binding sites for this factor [
<xref ref-type="bibr" rid="CR34">34</xref>
]. However, this spectrum of sites may not accurately represent GLI binding site preferences in vivo. To begin to examine this, we performed a meta-analysis of four existing GLI-ChIP datasets. All of these datasets utilize transgenic mice carrying FLAG tagged GLI1 (GLI1
<sup>FLAG</sup>
) or GLI3 (GLI3
<sup>FLAG</sup>
) in the ROSA26 locus, activated by Cre recombination, in four different tissue contexts: limb bud development (LD) [
<xref ref-type="bibr" rid="CR27">27</xref>
], cerebellum development (CD), medulloblastoma (MB) resulting from Hh signaling overexpression [
<xref ref-type="bibr" rid="CR29">29</xref>
], and neural progenitor cells (NP) [
<xref ref-type="bibr" rid="CR25">25</xref>
]. An additional study of neural progenitors [
<xref ref-type="bibr" rid="CR26">26</xref>
] was excluded from analysis since it contained a low number of significant peaks and mirrors the same experimental conditions as the NP dataset [
<xref ref-type="bibr" rid="CR25">25</xref>
]. For each of the datasets, the reported percentage of ChIP peak sequences with GLI binding motifs (GBM) was as follows: LD 55 %, CD 26 %, MB 46 %, and NP 91 %. However, the definition of GBM was not the same across all datasets: one study allowed only two mismatches from the consensus [
<xref ref-type="bibr" rid="CR27">27</xref>
] and others generated a GLI motif
<italic>de novo</italic>
based on the sequences of recovered peaks [
<xref ref-type="bibr" rid="CR25">25</xref>
,
<xref ref-type="bibr" rid="CR29">29</xref>
].</p>
<p>To collate the spectrum of GBM observed in all four datasets, we applied a
<italic>de novo</italic>
motif enrichment analysis to each dataset individually [
<xref ref-type="bibr" rid="CR35">35</xref>
]. Sequences that contained at least one site that matched the
<italic>de novo</italic>
motifs were removed from the dataset. The remaining sequences were analyzed for residual motifs that resembled a GBM using DREME [
<xref ref-type="bibr" rid="CR36">36</xref>
] and Tomtom [
<xref ref-type="bibr" rid="CR37">37</xref>
] (see
<xref rid="Sec10" ref-type="sec">Methods</xref>
). This resulted in 548 putative GBM (12-mers) (Additional file
<xref rid="MOESM1" ref-type="media">1</xref>
: Table S1), encompassing the range of GBM that are present in existing ChIP data. This set therefore represents a collection of likely genomic GLI binding sites, although some functional GLI binding sites in vivo could be absent from this set and some false positive sites may be included. Each 12-mer was classified as high confidence (HC), medium confidence (MC), or low confidence (LC) if it was found within sequences from all four datasets, two to three datasets, or one dataset, respectively. The sequence logos [
<xref ref-type="bibr" rid="CR38">38</xref>
] for each classification, provided in Fig. 
<xref rid="Fig1" ref-type="fig">1a</xref>
, show a nearly absolute representation of CCxC in positions 4–7 for all sites. Indeed, concordant (C and C or G and G) nucleotides at the 5
<sup>th</sup>
and 7
<sup>th</sup>
position were previously found to be required for GLI binding [
<xref ref-type="bibr" rid="CR39">39</xref>
]. Interestingly, for high confidence sites, there is no variation at 5 of the 12 positions, including the 5
<sup>th</sup>
and 7
<sup>th</sup>
positions (xGxCCxCxCxxx).
<fig id="Fig1">
<label>Fig. 1</label>
<caption>
<p>Definition of GLI binding motifs (GBM) and characterization of GLI
<sup>FLAG</sup>
datastets.
<bold>a</bold>
Sequence logos (Weblogo) of 12-mer GBM. High confidence (HC), medium confidence (MC), and low confidence (LC) GBM are found in all four (HC), three or two (MC) or one (LC) datasets.
<bold>b</bold>
Relative frequency of peaks containing one or more GBM in the GLI
<sup>FLAG</sup>
ChIP-chip LD (
<italic>red</italic>
), CD (
<italic>blue</italic>
), MB (
<italic>green</italic>
), and GLI
<sup>FLAG</sup>
ChIP-seq NP (
<italic>yellow</italic>
) datasets. A high proportion of sequences contain only one GBM.
<bold>c</bold>
Overlap of sequences identified by all four GLI
<sup>FLAG</sup>
datasets; only 26 individual peaks are found in all contexts</p>
</caption>
<graphic xlink:href="12861_2016_106_Fig1_HTML" id="MO1"></graphic>
</fig>
</p>
<p>Using the recommended matrix similarity score cutoff of 81 % overall matrix similarity to the optimal consensus GLI site as defined by the in vitro DNA binding assay [
<xref ref-type="bibr" rid="CR34">34</xref>
] results in 1,432,161 putative GLI TFBS across the mouse genome. This is substantially more than the 191,745 found using the new GLI library defined by the in vivo ChIP studies. However, several of the sites predicted by the in vitro binding studies do not contain the concordant (C and C or G and G) nucleotides at the 5
<sup>th</sup>
and 7
<sup>th</sup>
position. Thus, the newly generated GLI library (Fig. 
<xref rid="Fig1" ref-type="fig">1a</xref>
) may more accurately represent functional GBM.</p>
<p>According to this new library of 548 GBM, 41 % of LD, 27 % of CD, 32 % of MB, and 80 % of NP peaks contain putative GLI binding sites. For those sequences that contain a GBM, the vast majority contain only a single site (85.5 % LD, 90.7 % CD, 88.0 % MB, 89.8 % NP) (Fig. 
<xref rid="Fig1" ref-type="fig">1b</xref>
). The overlap of genomic binding regions among datasets is shown in Fig. 
<xref rid="Fig1" ref-type="fig">1c</xref>
; only 26 genomic coordinates are shared among all datasets (Additional file
<xref rid="MOESM2" ref-type="media">2</xref>
: Table S2). Since pathway components must respond to the Hh signal in all tissues, it might be expected that this common response would be integrated by a single enhancer. However, only three of the 26 shared regions are annotated to known Hh pathway components (
<italic>Gli1</italic>
,
<italic>Ptch1</italic>
,
<italic>Ptch2</italic>
). Indeed, for
<italic>Boc</italic>
and
<italic>Hhip</italic>
, distinct genomic GLI binding regions are found in different datasets. This result suggests that some of these genes may have multiple enhancers that work to transduce the Hh signal in different tissue contexts.</p>
<p>Given this apparent complexity in regulatory regions, we next asked whether the existing datasets of ChIP peaks might contain additional sequence information that could be used to predict the location of other Hh-responsive enhancers in the mouse genome. A machine learning approach was employed to test this question.</p>
</sec>
<sec id="Sec4">
<title>Assessment of kmer-SVM performance and prediction</title>
<sec id="Sec5">
<title>Kmer-SVM assessment of classification using GLI
<sup>FLAG</sup>
datasets</title>
<p>For each GLI
<sup>FLAG</sup>
dataset, only sequences with at least one GBM (wGBM, meaning
<italic>with GBM</italic>
) were used. This was done since a high proportion of the ChIP-chip datasets did not contain a putative GLI binding site as defined in the original papers (LD 55 %, CD 26 %, MB 46 %, and NP 91 %) or by our assessment (41 % of LD, 27 % of CD, 32 % of MB, and 80 % of NP peaks). Each individual dataset was submitted to kmer-SVM and the ability of each classifier to correctly label a candidate sequence as positive was assessed. Background sequences were randomly selected from the genome, but matched for GC content with the positive set.</p>
<p>Kmer-SVM randomly divides the data as follows: 80 % of the sequences are used as a training set and 20 % are used as a testing set. The ability of the classifier built with the training set to accurately identify the members of the remaining 20 % testing set is then assessed. This is repeated five times, each with a different random division of the data. Receiver operating characteristic (ROC) curves and precision recall curves (PRC) are used to assess the success of the classifier to correctly label regions in the testing set as positive (see
<xref rid="Sec10" ref-type="sec">Methods</xref>
).</p>
<p>ROC curves display the cumulative distribution of the true positive rate compared to the false positive rate. This characteristic assesses how well the classifier is able to label the positive sequences from the test set. The area under the curve was 0.898 for LDwGBM (Fig. 
<xref rid="Fig2" ref-type="fig">2a</xref>
), 0.856 for CDwGBM (Fig. 
<xref rid="Fig2" ref-type="fig">2b</xref>
), 0.862 for MBwGBM (Fig. 
<xref rid="Fig2" ref-type="fig">2c</xref>
) and 0.976 for NPwGBM (Fig. 
<xref rid="Fig2" ref-type="fig">2d</xref>
). Thus, the classifier performs best in LD and NP datasets.
<fig id="Fig2">
<label>Fig. 2</label>
<caption>
<p>Assessment of classification capability of kmer-SVM trained GLI
<sup>FLAG</sup>
datasets containing sequences with at least one GBM. For all curves, each dataset is randomly split into 80 % for training and 20 % for prediction and the prediction is repeated five times (
<italic>represented by individual lines</italic>
). Plots assess the likelihood that the specified classifier can successfully predict sequences that have at least one GBM as positive or negative.
<bold>a</bold>
-
<bold>d</bold>
ROC plots depicting true positive rates (TPR) and false positive rates (FPR). Area under the curve (AUC) scores as calculated by kmer-SVM are: 0.89 for LDwGBM (
<bold>a</bold>
), 0.85 for CDwGBM (
<bold>b</bold>
), 0.86 for MBwGBM (
<bold>c</bold>
) and 0.97 for NPwGBM (
<bold>d</bold>
) datasets.
<bold>e</bold>
-
<bold>h</bold>
Precision recall curves depicting the positive predictive value (PPV), calculated as true positive / (true positive + false positive), versus the TPR. AUC of 0.75 for LDwGBM (
<bold>e</bold>
) and 0.88 for NPwGBM (
<bold>h</bold>
) indicate reasonable confidence in the classification while AUC of 0.49 for CDwGBM (
<bold>f</bold>
) and 0.55 for MBwGBM (
<bold>g</bold>
) indicate a low probability that the region is correctly labeled when the sequence is classified as positive</p>
</caption>
<graphic xlink:href="12861_2016_106_Fig2_HTML" id="MO2"></graphic>
</fig>
</p>
<p>PRC displays the predictive value against the true positive rate and represents the accuracy of the labeling. The PRC plots indicated high values for LDwGBM (AUC = 0.753) and NPwGBM (AUC = 0.880) but low values for CDwGBM (AUC = 0.490) (Fig. 
<xref rid="Fig2" ref-type="fig">2f</xref>
) and MBwGBM (AUC = 0.546) (Fig. 
<xref rid="Fig2" ref-type="fig">2e</xref>
-
<xref rid="Fig2" ref-type="fig">h</xref>
). The ROC and PRC plots for LDwGBM and NPwGBM datasets suggested that the classifier sequence features used were able to distinguish between positive and negative groups with a low level of false labeling.</p>
</sec>
<sec id="Sec6">
<title>Predictions</title>
<p>The classifiers for LDwGBM and NPwGBM were then individually run on 600 bp of sequence centered on every GBM in the mouse genome (191,745, as determined using the new GBM from in vivo data, described above). Use of both the LDwGBM and NPwGBM datasets for prediction incorporated data from the GLI1
<sup>FLAG</sup>
(predominately activator) and GLI3
<sup>FLAG</sup>
(predominantly repressor) transcription factors in two diverse contexts (neuronal precursor and limb development).</p>
<p>The length of 600 bp was selected based on motif enrichment analysis of the LD and NP datasets using MEME-ChIP [
<xref ref-type="bibr" rid="CR40">40</xref>
] and Centrimo [
<xref ref-type="bibr" rid="CR41">41</xref>
]. This analysis showed that, within the ChIP-chip LD dataset, enrichment for the location of GLI motifs (green line) has a broad profile that spans 200 bp to either side of the midpoint (Additional file
<xref rid="MOESM3" ref-type="media">3</xref>
: Figure S1A). The GLI motif has a narrower profile in the NP data, a feature that is expected for ChIP-seq (Additional file
<xref rid="MOESM3" ref-type="media">3</xref>
: Figure S1B). The profile for the Sox motif (blue line), an established tissue specific GLI cofactor [
<xref ref-type="bibr" rid="CR25">25</xref>
], shows an enrichment peak that is centered around 200 bp on either side of the midpoint (Additional file
<xref rid="MOESM3" ref-type="media">3</xref>
: Figure S1B) and suggests that cofactors for Hh may reside outside of the immediate vicinity of a GLI binding site. We therefore used 600 bp to capture both common Hh features as well as potential context specific sequence.</p>
<p>For the LDwGBM classifier, scores ranged from −4.33 to 12.00 with 18.4 % of the 191,745 analyzed genomic regions scoring as positive (Score > 0). The NPwGBM results ranged from −2.54 to 5.48 with 5.7 % positive (Fig. 
<xref rid="Fig3" ref-type="fig">3a</xref>
; Additional file
<xref rid="MOESM4" ref-type="media">4</xref>
: Table S3). The categorization of a sequence is dependent on the sign of the score and the weight of the value is less important than the ranking. Overall, the correlation between scores for individual genomic regions calculated by the LDwGBM and NPwGBM classifiers is poor (0.68 Pearson) (Fig. 
<xref rid="Fig3" ref-type="fig">3b</xref>
). However, if only sequences with positive scores are considered, the correlation improves (0.85 Pearson). If scores are restricted to values indicating only the high confidence scores (calculated posterior probabilities = 1.0, Additional file
<xref rid="MOESM5" ref-type="media">5</xref>
: Figure S2), the values are very well correlated (0.96 Pearson). In total, 8627 genomic regions were predicted as Hh enhancers by both classifiers. Of those, 1198 regions (14 %) overlapped at least one peak in the four GLI
<sup>FLAG</sup>
datasets. Among high confidence scores (LDwGBM: 5951 ≥ 1, NPwGM: 547 ≥ 1) 528 genomic regions were shared between the two datasets and 187 of these (35 %) overlapped with peaks from at least one of the four GLI
<sup>FLAG</sup>
datasets. All of the scored regions are listed in Additional file
<xref rid="MOESM4" ref-type="media">4</xref>
: Table S3.
<fig id="Fig3">
<label>Fig. 3</label>
<caption>
<p>Assessment of genomic kmer-SVM predictions using classifiers trained on LDwGBM and NPwGBM datasets.
<bold>a</bold>
All genomic sequences matching the restricted 548 GBM 12-mers (wGMB) were identified and the 600 bp surrounding each GBM were assessed and scored using the kmer-SVM classifier that was trained on each of the two datasets.
<bold>b</bold>
Correlation plot depicting the relationship between LDwGBM and NPwGBM scores; scores >1 are highly correlated in the two datasets.
<bold>c</bold>
GLI motif generated from overlapping high weighted k-mers shared between LDwGBM and NPwGBM classifiers.
<bold>d</bold>
High weighted k-mers (identified by Tomtom) represented in either LDwGBM (Tcf and Zfp) or NPwGBM (Fox and Sox)</p>
</caption>
<graphic xlink:href="12861_2016_106_Fig3_HTML" id="MO3"></graphic>
</fig>
</p>
</sec>
<sec id="Sec7">
<title>Evaluation of predictions</title>
<p>To assess whether kmer-SVM predictions were likely to represent Hh enhancers, we examined whether the predicted genomic regions overlapped publically available open chromatin and enhancer histone marks in tissues that were likely to be transducing Hh signals. We first examined the DNaseI hypersensitive profile collected from mouse mesoderm at E11.5 at genomic regions that were predicted with high confidence as positive (1 ≤ score; posterior probability = 1) or negative (−1 ≥ score; posterior probability = 0) (Additional file
<xref rid="MOESM5" ref-type="media">5</xref>
: Figure S2). A significantly higher proportion of overlap was found with the predicted positive regions than with predicted negative regions (Z-Score = 2.8332;
<italic>p</italic>
-value < 0.05) (Additional file
<xref rid="MOESM6" ref-type="media">6</xref>
: Table S4). We also examined publically available H3K4me1, H3K4me3, and H3K27ac ENCODE data collected from heart and liver at E14.5. Monomethylated H3K4 (H3K4me1) and histone H3 acetyl Lys27 H3K27ac [
<xref ref-type="bibr" rid="CR42">42</xref>
] were used as enhancer markers while trimethylated H3K4 (H3K4me3) was expected to be depleted in enhancer regions [
<xref ref-type="bibr" rid="CR43">43</xref>
]. Although Hh signaling is active during early development of both tissues, available in situ analysis for GLI1 (GenePaint: EN1215) [
<xref ref-type="bibr" rid="CR44">44</xref>
,
<xref ref-type="bibr" rid="CR45">45</xref>
] shows GLI1 expression in liver but not heart at E14.5 (Additional file
<xref rid="MOESM7" ref-type="media">7</xref>
: Figure S3). Thus, we expected signals for both H3K4me1 (poised enhancer) and H3K27ac (active enhancer) to be enriched in the predicted positive regions in liver but not heart at this time point. Indeed this was the case: H3K4me1 (Z-Score = 2.5511;
<italic>p</italic>
-value < 0.01) and H3K27ac (Z-Score = 8.076;
<italic>p</italic>
-value <0.01), with no significant difference in H3K4me3 when predicted positive regions were compared to predicted negative regions. As expected, the heart data did not show enrichment for H3K4me1 or H3K27ac. Together, the results (summarized in Additional file
<xref rid="MOESM6" ref-type="media">6</xref>
: Table S4) are consistent with the conclusion that the kmer-SVM classification correctly identifies Hh enhancer regions.</p>
<p>Next, we evaluated the sequence features, or k-mers, that kmer-SVM identified as primary components of Hh enhancer regions. The weights of k-mers are calculated during the SVM training and reflect the contribution of the k-mer to categorization of a sequence. Weights can be positive or negative and the sum of the weights of iterative k-mers across a sequence comprise the overall score of that sequence. Not surprisingly, alignment of k-mers with high scoring weights shared between both datasets returned a motif that strongly resembles the GBM (Fig. 
<xref rid="Fig3" ref-type="fig">3c</xref>
). Unique high weighted k-mers that occurred in each individual dataset represented potential context specific features. An E-box motif was identified for the LDwGBM dataset while a Sox motif was returned for NPwGBM (Fig. 
<xref rid="Fig3" ref-type="fig">3d</xref>
). Negative weights that occurred in both datasets include AC and ACC repeats as well as other C rich sequences.</p>
</sec>
</sec>
<sec id="Sec8">
<title>Functional verification of GLI-dependent enhancer activity</title>
<p>Predicted genomic regions were annotated to the two nearest genes using GREAT [
<xref ref-type="bibr" rid="CR46">46</xref>
]. Because our goal was to identify enhancers for Hh pathway components, we selected a subset of predictions that were positive in both the LDwGBM and NPwGBM datasets and that were annotated to members of the GO:0007224 Smoothened signaling pathway gene set. Because Hh pathway components are required for active Hh signaling, we reasoned that enhancers annotated to these genes would be more likely to function in any tissue that transduces Hh signal. Therefore, high scoring regions annotated to different members of the GO:0007224 gene set that were readily cloned were functionally tested for enhancer activity. Two previously known Hh enhancers for
<italic>Ptch1</italic>
and
<italic>Ptch2</italic>
appeared on this list [
<xref ref-type="bibr" rid="CR24">24</xref>
,
<xref ref-type="bibr" rid="CR28">28</xref>
]. Interestingly, an established
<italic>Gli1</italic>
regulatory region was not predicted [
<xref ref-type="bibr" rid="CR15">15</xref>
]. The test set consisted of genomic regions annotated to
<italic>Boc</italic>
,
<italic>Gli3</italic>
,
<italic>Hhip</italic>
,
<italic>Hipk2</italic>
,
<italic>Ptch1</italic>
,
<italic>Scube1</italic>
,
<italic>Shh</italic>
, and
<italic>Tgfbr2</italic>
. An additional region, annotated to
<italic>Dpp6</italic>
(near
<italic>Shh</italic>
) was also tested (Table 
<xref rid="Tab1" ref-type="table">1</xref>
).
<table-wrap id="Tab1">
<label>Table 1</label>
<caption>
<p>Assessment of predicted Hh enhancer regions</p>
</caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th>Annotated gene</th>
<th>Genomic coordinates (mm9)</th>
<th>Hh responsive</th>
<th>LD</th>
<th>CD</th>
<th>MB</th>
<th>NP</th>
</tr>
</thead>
<tbody>
<tr>
<td>
<italic>Ptch2</italic>
</td>
<td>chr4:116,767,757-116,769,455</td>
<td>+</td>
<td>+</td>
<td>+</td>
<td>+</td>
<td>+</td>
</tr>
<tr>
<td>
<italic>Boc</italic>
</td>
<td>chr16:44,502,136-44,503,346</td>
<td>+</td>
<td>-</td>
<td>+</td>
<td>+</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Dpp6</italic>
</td>
<td>chr5:27,248,056-27,249,266</td>
<td>+</td>
<td>+</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Gli3</italic>
</td>
<td>chr13:15,764,694-15,765,904</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Hhip</italic>
</td>
<td>chr8:82,838,195-82,839,405</td>
<td>+</td>
<td>+</td>
<td>-</td>
<td>-</td>
<td>+</td>
</tr>
<tr>
<td>
<italic>Hipk2</italic>
</td>
<td>chr6:38,614,001-38,615,211</td>
<td>+</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Ptch1</italic>
</td>
<td>chr13:63669992-63671202</td>
<td>+</td>
<td>+</td>
<td>-</td>
<td>+</td>
<td>+</td>
</tr>
<tr>
<td>
<italic>Scube1</italic>
</td>
<td>chr15:83503053-83504263</td>
<td>+</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Shh</italic>
</td>
<td>chr5:28832033-28833243</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
<tr>
<td>
<italic>Tgfbr2</italic>
</td>
<td>chr9:116,151,184-116,152,394</td>
<td>+</td>
<td>-</td>
<td>-</td>
<td>-</td>
<td>-</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>Seven of the nine regions predicted to be GLI-driven enhancers were indeed determined to be Hh responsive and GLI binding site dependent in a cell culture assay. Overlap of the predicted regions with peaks from the GLI
<sup>FLAG</sup>
ChIP datasets (LD, CD, MB, NP) is indicated by the plus sign.
<italic>Boc</italic>
,
<italic>Hipk2</italic>
,
<italic>Scube1</italic>
and
<italic>Tgfbr2</italic>
were predicted by kmer-SVM and found to be positive, even though those regions do not overlap any of the peaks identified in the LD or NP datasets used to generate the kmer-SVM classifier</p>
</table-wrap-foot>
</table-wrap>
</p>
<p>The nine genomic regions were screened for Hh responsiveness in a cell culture assay that has been previously used to validate Hh enhancers [
<xref ref-type="bibr" rid="CR47">47</xref>
,
<xref ref-type="bibr" rid="CR48">48</xref>
]. In this assay, C2C12 cells are transiently transfected with constructs containing the putative Hh regulatory region cloned upstream of a minimal promoter that drives luciferase expression (see
<xref rid="Sec10" ref-type="sec">Methods</xref>
). To induce a Hh response, enhancer reporter constructs were co-transfected with a plasmid that drives constitutive expression of GLI1 in mammalian cells [
<xref ref-type="bibr" rid="CR26">26</xref>
]. For those candidate enhancers that demonstrated apparent Hh activation, GLI-dependent activity was confirmed by retesting after mutagenesis of the GLI binding sites (GKO). Either complete loss of enhancer activity or attenuation of response in GKO sequence was considered GLI-dependent. The established Hh enhancer region for
<italic>Ptch2</italic>
was used as a positive control [
<xref ref-type="bibr" rid="CR28">28</xref>
].</p>
<p>For the 9 regulatory regions annotated to Hh pathway component genes, 7 exhibited Hh activation that was directly dependent on a GLI binding site (Fig. 
<xref rid="Fig4" ref-type="fig">4</xref>
). The
<italic>Ptch2</italic>
positive control region showed a complete loss of Hh response after mutation of the GLI site (Fig. 
<xref rid="Fig4" ref-type="fig">4</xref>
) as did the regulatory regions annotated to
<italic>Hhip</italic>
,
<italic>Hipk2</italic>
,
<italic>Ptch1</italic>
, and
<italic>Scube1</italic>
. Regions annotated to
<italic>Boc</italic>
,
<italic>Dpp6</italic>
, and
<italic>Tgfbr2</italic>
showed a significant decrease in Hh activation upon GLI TFBS mutation, but not a complete loss of response. This suggests that additional regulatory inputs influence the activity of these enhancers. Neither the
<italic>Gli3</italic>
nor
<italic>Shh</italic>
region exhibited Hh dependent enhancer activity. However, we cannot rule out the possibility that these regions might be positive if examined in a different cellular context [
<xref ref-type="bibr" rid="CR26">26</xref>
].
<fig id="Fig4">
<label>Fig. 4</label>
<caption>
<p>Functional verification of GLI-dependent enhancer activity. Putative regulatory regions were cloned upstream of a minimal promoter co-transfected into C2C12 cells, along with a GLI1 expression vector. Relative activity is plotted (stimulated/basal). The Ptch2 region is used as a positive control. Novel regions annotated to the
<italic>Boc</italic>
,
<italic>Dpp6</italic>
,
<italic>Hhip</italic>
,
<italic>Hipk2</italic>
,
<italic>Ptch1</italic>
,
<italic>Scube1</italic>
, and
<italic>Tgfbr2</italic>
loci exhibit upregulation in response to GLI1 co-transfection (
<italic>green</italic>
). Dependence on GLI was assessed by mutating all putative GLI TFBS (GKO) within the sequence and retesting in the assay (
<italic>gray</italic>
)</p>
</caption>
<graphic xlink:href="12861_2016_106_Fig4_HTML" id="MO4"></graphic>
</fig>
</p>
<p>Two previous
<italic>in silico</italic>
methods have been described for the identification of Hh-regulated enhancers in vertebrates: Enhancer Element Locator (EEL) and Module Cluster Analysis (MCA). EEL analyzes the estimated energy of a single transcription factor binding event, as well as the possible interaction between adjacent, highly conserved transcription factor binding sites, to detect potential enhancers [
<xref ref-type="bibr" rid="CR34">34</xref>
]. In contrast, MCA utilizes Poisson modeling to determine the relative enrichment of binding sites in highly conserved, non-coding sequence and, thereby, identify putative enhancers [
<xref ref-type="bibr" rid="CR26">26</xref>
]. While both of these methods have had some success (~25 %) in detecting putative Hh-regulated enhancers, a disadvantage to these approaches is that the search is limited to regions of high sequence conservation and to regions close to promoters. In contrast, kmer-SVM approach used here employs a genome-wide empirical analysis to locate regions that contain sequence features predictive of Hh enhancer function. Though these predictions miss one of the three known Hh pathway enhancers [
<xref ref-type="bibr" rid="CR15">15</xref>
] indicating that the algorithm does not capture all Hh-driven enhancers, the high success rate (78 %) of the kmer-SVM predictions far exceed the previous prediction rates for EEL or MCA. It is important to note, however, that the predictions tested here are all for pathway genes, which may have a unique signature. It would be necessary to test additional enhancers in tissue-specific assay systems (e.g., transgenic mice) to determine the overall success of this method in identification of tissue-specific enhancers.</p>
<p>To examine the impact of k-mers that contribute to predictions of the positive regions, weights were plotted across each of the sequences. The
<italic>Ptch2</italic>
sequence (Fig. 
<xref rid="Fig5" ref-type="fig">5a</xref>
), a known enhancer region [
<xref ref-type="bibr" rid="CR28">28</xref>
], contained matching profiles for 8-mers predicted from LDwGBM (red) and NPwGBM (yellow) that are GC rich and similar to the GBM k-mers However, mutation of the single GLI TFBS (Fig. 
<xref rid="Fig5" ref-type="fig">5b</xref>
, green box) ablates the Hh response, indicating that the presence of this GBM is required to transduce Hh signaling. Mutation of the GBMs annotated by the green boxes for
<italic>Hhip</italic>
,
<italic>Hipk2</italic>
,
<italic>Ptch1</italic>
, and
<italic>Scube1</italic>
(Fig. 
<xref rid="Fig5" ref-type="fig">5e</xref>
,
<xref rid="Fig5" ref-type="fig">f</xref>
,
<xref rid="Fig5" ref-type="fig">g</xref>
,
<xref rid="Fig5" ref-type="fig">h</xref>
) is also sufficient to abrogate Hh signal transduction. For
<italic>Boc</italic>
,
<italic>Dpp6</italic>
, and
<italic>Tgfbr2</italic>
, which show enhancer activity that remains after ablation of the GBM (Fig. 
<xref rid="Fig5" ref-type="fig">5b</xref>
,
<xref rid="Fig5" ref-type="fig">c</xref>
,
<xref rid="Fig5" ref-type="fig">j</xref>
), there were no sequence characteristics that were indicative of a shared feature responsible for this remaining response to induction by GLI1. The Boc profile was the only one that contained a high weighted k-mer (annotated with an asterisk) that was unique to the LDwGBM (red) profile. This k-mer was similar to a Krox motif (Tomtom
<italic>p</italic>
-value < 0.004) [
<xref ref-type="bibr" rid="CR37">37</xref>
] and may be enriched in the LDwGBM dataset as a context specific transcription factor, since it has roles in limb development [
<xref ref-type="bibr" rid="CR49">49</xref>
,
<xref ref-type="bibr" rid="CR50">50</xref>
]. In general, most of the tested regions contain distributed high weighted 8-mers in addition to the central GBM and had profiles that contained consistent peaks in both LDwGM and NPwGBM datasets (Fig. 
<xref rid="Fig5" ref-type="fig">5</xref>
).
<fig id="Fig5">
<label>Fig. 5</label>
<caption>
<p>K-mer weights plotted across sequences that show enhancer activity. Diagrams were generated in UCSC Genome browser and show coordinate information for regions annotated to
<italic>Ptch2</italic>
(
<bold>a</bold>
),
<italic>Boc </italic>
(
<bold>b</bold>
),
<italic>Dpp6</italic>
(
<bold>c</bold>
),
<italic>Gli3</italic>
(
<bold>d</bold>
),
<italic> Hhip </italic>
(
<bold>e</bold>
),
<italic>Hipk2</italic>
(
<bold>f</bold>
),
<italic> Ptch1 </italic>
(
<bold>g</bold>
),
<italic>Scube1</italic>
(
<bold>h</bold>
),
<italic>Shh </italic>
(
<bold>i</bold>
) and
<italic>Tgfbr2</italic>
(
<bold>j</bold>
). Green boxes represent GBM. Weights for LDwGBM and NPwGBM are represented by the red and yellow lines, respectively. Refseq gene annotations are represented in blue. A putative Krox-20 TFBS (*) that has a high weight in the LDwGBM classifier but not the NPwGBM classifier occurs in the sequence annotated to
<italic>Boc</italic>
. Note that most sequences show weighted k-mers located several hundred bp from the central GBM, suggesting that sequence motifs that predict Hh enhancer activity may be distributed</p>
</caption>
<graphic xlink:href="12861_2016_106_Fig5_HTML" id="MO5"></graphic>
</fig>
</p>
</sec>
</sec>
<sec id="Sec9" sec-type="conclusion">
<title>Conclusions</title>
<p>We have utilized the kmer-SVM machine learning approach to examine four existing GLI ChIP databases and to generate classifiers that can empirically predict functional Hh enhancers from genomic DNA. The analysis was facilitated by a new GBS library generated from a meta-analysis of genomic regions identified in in vivo binding studies [
<xref ref-type="bibr" rid="CR25">25</xref>
–
<xref ref-type="bibr" rid="CR27">27</xref>
,
<xref ref-type="bibr" rid="CR29">29</xref>
]. When compared to the previous library derived from in vitro binding studies [
<xref ref-type="bibr" rid="CR34">34</xref>
], this new GBS library resulted in the identification of a subset (approximately 10 %) of potential GLI transcription factor binding sites across the mouse genome. Of nine predicted Hh target regulatory regions tested, seven were functionally verified as GLI-dependent. All of the tested regions were annotated to genes believed to be components of the Hh pathway and important determinants of the Hh response. Given the high success rate of Hh enhancer prediction in this small scale screen (78 %), it is quite possible that a large number of the other 37,000 predicted regions (Score > 0 in Additional file
<xref rid="MOESM4" ref-type="media">4</xref>
: Table S3) may harbor Hh enhancer activity.</p>
</sec>
<sec id="Sec10" sec-type="materials|methods">
<title>Methods</title>
<sec id="Sec11">
<title>Computing resources</title>
<p>Except where otherwise indicated, all computational steps were performed using custom Perl and R scripts.</p>
</sec>
<sec id="Sec12">
<title>Publically available datasets</title>
<p>Genomic analysis was conducted on chromosomes 1 to 19, X and Y of mouse build mm9. Mouse ENCODE data [
<xref ref-type="bibr" rid="CR51">51</xref>
] comprising open chromatin DNaseI data that was collected at embryonic day 11.5 in the mesoderm and histone (H3K4me1, H3K4me3, and H3K27ac) data collected from embryonic day 14.5 for heart and liver were downloaded from the UCSC genome repository (goldenPath).</p>
</sec>
<sec id="Sec13">
<title>Definition of putative GLI binding motifs</title>
<p>The library of putative GLI binding motifs (GBM) was compiled using
<italic>de novo</italic>
motif analysis [
<xref ref-type="bibr" rid="CR35">35</xref>
] on each of the individual GLI
<sup>FLAG</sup>
datasets iteratively. Sequences that contained a GBM were removed from the dataset and the remaining sequences were analyzed for enriched motifs using DREME [
<xref ref-type="bibr" rid="CR36">36</xref>
]. If Tomtom [
<xref ref-type="bibr" rid="CR37">37</xref>
] returned a GLI motif, the dataset was reanalyzed using HOMER [
<xref ref-type="bibr" rid="CR35">35</xref>
]. The process continued until no residual GBM remained enriched in the dataset. Confidence in the GBM was classified as high (HC) if it was shared across sequences from all four GLI
<sup>FLAG</sup>
datasets, medium (MC) if it was found in two or three datasets, and low (LC) if it only occurred in one.</p>
</sec>
<sec id="Sec14">
<title>kmer-SVM parameters and evaluation of classifiers</title>
<p>Training was run on the Beer lab webserver (
<ext-link ext-link-type="uri" xlink:href="http://kmersvm.beerlab.org/">http://kmersvm.beerlab.org/</ext-link>
), using a k-mer of length of 8.</p>
<p>Performance of the classifier built by kmer-SVM’s training algorithm was assessed using Receiver Operating Characteristic (ROC) and Precision-Recall curves (PRC) generated within the kmer-SVM program. True positive, true negative, false positive and false negative counts were generated by segregating the sets of positive and negative sequences into a training set (80 % of the sequences) and a testing set (the remaining 20 % of the sequences). Each member (individual sequence) of the testing set that is correctly annotated as positive increases the true positive count while an incorrect prediction of a positive sequence as negative increases the false negative count. ROC curves asses the sensitivity and specificity of the classifier output. A steep curve with a high area under the curve (AUC) indicates a high true positive rate and a low false positive rate. PRC evaluate the accuracy and relevance of the classifier output. A high AUC indicates that the results have a low false positive rate (high precision) and a low false negative rate. The trained SVM is evaluated by assessing its ability to classify the testing set correctly. The classifier was assessed five times by resetting members in the training set and testing set.</p>
</sec>
<sec id="Sec15">
<title>Cloning of putative enhancer regions</title>
<p>Putative enhancers were amplified from C57BL/6 genomic DNA (supplied by Jackson Laboratory) using template-specific PCR primers (Additional file
<xref rid="MOESM8" ref-type="media">8</xref>
: Table S4). A CACC extension was added to the end of one primer to facilitate directional cloning. PCR fragments were cloned into the pENTR/D-TOPO vector using the standard kit (Invitrogen) and then shuttled into the pGL3-Promoter luciferase vector (Promega) using the Gateway® cloning system (Invitrogen). QuikChange mutagenesis (Stratagene) was used to mutate putative GLI binding sites by replacing the C in the 6
<sup>th</sup>
position to a G.</p>
</sec>
<sec id="Sec16">
<title>Luciferase assay</title>
<p>C2C12 cells (35,000) were plated per well on 12-well plates (10 % fetal bovine serum treated with penicillin, streptomycin and glutamate). After 24 h, cells were transfected, using lipofectamine, with 400 ng of the construct containing the putative enhancer region plus either a control vector or GLI1 (in equal molecular weight).
<italic>Renilla</italic>
(Promega pRL-CMV) was also included to normalize transfection efficiency. After an additional 24 h, cell media was changed to no serum to promote ciliogenesis [
<xref ref-type="bibr" rid="CR52">52</xref>
]. Cell lysate was collected after 48 h and measured for luciferase activity using the Dual-Luciferase® Reporter Assay System (Promega) on a Perkin Elmer Wallac Victor3 1420 Multilabel Counter. Three experimental replicates were collected for each condition.</p>
</sec>
</sec>
</body>
<back>
<app-group>
<app id="App1">
<sec id="Sec17">
<title>Additional files</title>
<p>
<media position="anchor" xlink:href="12861_2016_106_MOESM1_ESM.xlsx" id="MOESM1">
<label>Additional file 1: Table S1.</label>
<caption>
<p>GLI binding motif 12-mers. Library of 12-mers enriched in GLI
<sup>FLAG</sup>
datasets. Motifs were considered high confidence (HC) if the 12-mer occurred in all four GLI
<sup>FLAG</sup>
datasets, medium confidence (MC) if it occurred in 2 or 3, or low confidence (LC) if it only occurred in one dataset. (XLSX 17 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM2_ESM.xlsx" id="MOESM2">
<label>Additional file 2: Table S2.</label>
<caption>
<p>Peak coordinates that overlap between all four GLI
<sup>FLAG</sup>
datasets. The 26 regions that have shared peaks across all GLI
<sup>FLAG</sup>
datasets. Peaks are annotated to nearest genes. (XLSX 9 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM3_ESM.tif" id="MOESM3">
<label>Additional file 3: Figure S1.</label>
<caption>
<p>Determination of sequence length buffer surrounding the GBM. Plots depicting the positional distribution of the best GLI motif (green) were generated by submitting 300 bp of sequence surrounding the center of each peak to Centrimo. (A) LDwGBM shows a broad profile for the best GBM, consistent with ChIP-chip data. (B) The profile for the ChIP-seq sequences from NPwGBM is more narrow and suggests that most of the GBM fall within 240 bp around the center of the peak. In neural precursor cells, the motif for the GLI cofactor, Sox, has a profile that contains a central apex plus two additional summits at a distance of 240 bp on either side of the peak. This suggests that context-specific TF binding may occur outside the central peak region. (TIF 15871 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM4_ESM.xlsx" id="MOESM4">
<label>Additional file 4: Table S3.</label>
<caption>
<p>GLI
<sup>FLAG</sup>
dataset kmer-SVM scores. kmer-SVM scores for LDwGBM and NPwGBM datasets. (XLSX 20714 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM5_ESM.tif" id="MOESM5">
<label>Additional file 5: Figure S2.</label>
<caption>
<p>Posterior probability of kmer-SVM scores. Plots depicting the posterior probabilities assigned to scores for both (A) LDwGBM and (B) NPwGBM datasets. The graphs indicate that scores above 1 have a high confidence of being Hh regulatory regions. (TIF 14263 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM6_ESM.xlsx" id="MOESM6">
<label>Additional file 6: Table S4.</label>
<caption>
<p>Overlap of predicted high confidence positive and negative regions with embryonic open chromatin. Tabulation of the number of genomic regions predicted by both LDwGBM and NPwGBM that are classified with high confidence as Hh enhancer regions or as nonregulatory regions that overlap with mesoderm DNaseI (E 11.5) or enhancer markers (E14.5). (XLSX 9 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM7_ESM.jpg" id="MOESM7">
<label>Additional file 7: Figure S3.</label>
<caption>
<p>Expression of GLI1 within E14.5 mouse embryo. In situ hybridization of GLI1 (image from genepaint.org, EN1215) showing active Hh signaling at E14.5 in liver but not heart. (JPG 216 kb)</p>
</caption>
</media>
<media position="anchor" xlink:href="12861_2016_106_MOESM8_ESM.xlsx" id="MOESM8">
<label>Additional file 8: Table S5.</label>
<caption>
<p>PCR primers for amplification of mouse genomic regions. Mouse genomic coordinates (mm9) for primer sequence used to amplify candidate regions. (XLSX 10 kb)</p>
</caption>
</media>
</p>
</sec>
</app>
</app-group>
<glossary>
<title>Abbreviations</title>
<def-list>
<def-item>
<term>E14.5</term>
<def>
<p>embryonic day 14.5 (mouse)</p>
</def>
</def-item>
<def-item>
<term>GBM</term>
<def>
<p>GLI binding motif</p>
</def>
</def-item>
<def-item>
<term>Hh</term>
<def>
<p>Hedgehog</p>
</def>
</def-item>
<def-item>
<term>HPC</term>
<def>
<p>Hedgehog pathway component</p>
</def>
</def-item>
<def-item>
<term>SVM</term>
<def>
<p>support vector machine</p>
</def>
</def-item>
<def-item>
<term>TFBS</term>
<def>
<p>transcription factor binding sites</p>
</def>
</def-item>
</def-list>
</glossary>
<fn-group>
<fn>
<p>
<bold>Competing interests</bold>
</p>
<p>The authors have no competing interests to declare.</p>
</fn>
<fn>
<p>
<bold>Authors’ contributions</bold>
</p>
<p>KG and DG conceived and designed the experiments. KG, GS, and KV performed the experiments. KG, KV, GS, and NR cloned and mutagenized constructs. KG analyzed the data. KG and DG wrote the manuscript. All authors reviewed the final manuscript. All authors read and approved the final manuscript.</p>
</fn>
</fn-group>
<ack>
<title>Acknowledgements</title>
<p>KG was supported by University of Michigan NIH Training Program in Basic and Translational Digestive Science (T32 DK094775). DG acknowledges support from P01 DK062041.</p>
</ack>
<ref-list id="Bib1">
<title>References</title>
<ref id="CR1">
<label>1.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Briscoe</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Thérond</surname>
<given-names>PP</given-names>
</name>
</person-group>
<article-title>The mechanisms of Hedgehog signalling and its roles in development and disease</article-title>
<source>Nat Rev Mol Cell Biol</source>
<year>2013</year>
<volume>14</volume>
<fpage>416</fpage>
<lpage>29</lpage>
<pub-id pub-id-type="doi">10.1038/nrm3598</pub-id>
<pub-id pub-id-type="pmid">23719536</pub-id>
</element-citation>
</ref>
<ref id="CR2">
<label>2.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Hui</surname>
<given-names>CC</given-names>
</name>
</person-group>
<article-title>Hedgehog signaling in development and cancer</article-title>
<source>Dev Cell</source>
<year>2008</year>
<volume>15</volume>
<fpage>801</fpage>
<lpage>12</lpage>
<pub-id pub-id-type="doi">10.1016/j.devcel.2008.11.010</pub-id>
<pub-id pub-id-type="pmid">19081070</pub-id>
</element-citation>
</ref>
<ref id="CR3">
<label>3.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>McMahon</surname>
<given-names>AP</given-names>
</name>
<name>
<surname>Ingham</surname>
<given-names>PW</given-names>
</name>
<name>
<surname>Tabin</surname>
<given-names>CJ</given-names>
</name>
</person-group>
<article-title>Developmental roles and clinical significance of hedgehog signaling</article-title>
<source>Curr Top Dev Biol</source>
<year>2003</year>
<volume>53</volume>
<fpage>1</fpage>
<lpage>114</lpage>
<pub-id pub-id-type="doi">10.1016/S0070-2153(03)53002-2</pub-id>
<pub-id pub-id-type="pmid">12509125</pub-id>
</element-citation>
</ref>
<ref id="CR4">
<label>4.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Petrova</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Joyner</surname>
<given-names>AL</given-names>
</name>
</person-group>
<article-title>Roles for Hedgehog signaling in adult organ homeostasis and repair</article-title>
<source>Development</source>
<year>2014</year>
<volume>141</volume>
<fpage>3445</fpage>
<lpage>57</lpage>
<pub-id pub-id-type="doi">10.1242/dev.083691</pub-id>
<pub-id pub-id-type="pmid">25183867</pub-id>
</element-citation>
</ref>
<ref id="CR5">
<label>5.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Varjosalo</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Taipale</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Hedgehog: functions and mechanisms</article-title>
<source>Genes Dev</source>
<year>2008</year>
<volume>22</volume>
<fpage>2454</fpage>
<lpage>72</lpage>
<pub-id pub-id-type="doi">10.1101/gad.1693608</pub-id>
<pub-id pub-id-type="pmid">18794343</pub-id>
</element-citation>
</ref>
<ref id="CR6">
<label>6.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Creanga</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Glenn</surname>
<given-names>TD</given-names>
</name>
<name>
<surname>Mann</surname>
<given-names>RK</given-names>
</name>
<name>
<surname>Saunders</surname>
<given-names>AM</given-names>
</name>
<name>
<surname>Talbot</surname>
<given-names>WS</given-names>
</name>
<name>
<surname>Beachy</surname>
<given-names>PA</given-names>
</name>
</person-group>
<article-title>Scube/You activity mediates release of dually lipid-modified Hedgehog signal in soluble form</article-title>
<source>Genes Dev</source>
<year>2012</year>
<volume>26</volume>
<fpage>1312</fpage>
<lpage>25</lpage>
<pub-id pub-id-type="doi">10.1101/gad.191866.112</pub-id>
<pub-id pub-id-type="pmid">22677548</pub-id>
</element-citation>
</ref>
<ref id="CR7">
<label>7.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jakobs</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Exner</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Schurmann</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Pickhinke</surname>
<given-names>U</given-names>
</name>
<name>
<surname>Bandari</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Ortmann</surname>
<given-names>C</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Scube2 enhances proteolytic Shh processing from the surface of Shh-producing cells</article-title>
<source>J Cell Sci</source>
<year>2014</year>
<volume>127</volume>
<fpage>1726</fpage>
<lpage>37</lpage>
<pub-id pub-id-type="doi">10.1242/jcs.137695</pub-id>
<pub-id pub-id-type="pmid">24522195</pub-id>
</element-citation>
</ref>
<ref id="CR8">
<label>8.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Burke</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Nellen</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Bellotto</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Hafen</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Senti</surname>
<given-names>KA</given-names>
</name>
<name>
<surname>Dickson</surname>
<given-names>BJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Dispatched, a novel sterol-sensing domain protein dedicated to the release of cholesterol-modified hedgehog from signaling cells</article-title>
<source>Cell</source>
<year>1999</year>
<volume>99</volume>
<fpage>803</fpage>
<lpage>15</lpage>
<pub-id pub-id-type="doi">10.1016/S0092-8674(00)81677-3</pub-id>
<pub-id pub-id-type="pmid">10619433</pub-id>
</element-citation>
</ref>
<ref id="CR9">
<label>9.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jeong</surname>
<given-names>J</given-names>
</name>
<name>
<surname>McMahon</surname>
<given-names>AP</given-names>
</name>
</person-group>
<article-title>Cholesterol modification of Hedgehog family proteins</article-title>
<source>J Clin Invest</source>
<year>2002</year>
<volume>110</volume>
<fpage>591</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1172/JCI0216506</pub-id>
<pub-id pub-id-type="pmid">12208857</pub-id>
</element-citation>
</ref>
<ref id="CR10">
<label>10.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Izzi</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Levesque</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Morin</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Laniel</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Wilkes</surname>
<given-names>BC</given-names>
</name>
<name>
<surname>Mille</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Boc and Gas1 each form distinct Shh receptor complexes with Ptch1 and are required for Shh-mediated cell proliferation</article-title>
<source>Dev Cell</source>
<year>2011</year>
<volume>20</volume>
<fpage>788</fpage>
<lpage>801</lpage>
<pub-id pub-id-type="doi">10.1016/j.devcel.2011.04.017</pub-id>
<pub-id pub-id-type="pmid">21664577</pub-id>
</element-citation>
</ref>
<ref id="CR11">
<label>11.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Allen</surname>
<given-names>BL</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>JY</given-names>
</name>
<name>
<surname>Izzi</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Althaus</surname>
<given-names>IW</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>JS</given-names>
</name>
<name>
<surname>Charron</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Overlapping roles and collective requirement for the coreceptors GAS1, CDO, and BOC in SHH pathway function</article-title>
<source>Dev Cell</source>
<year>2011</year>
<volume>20</volume>
<fpage>775</fpage>
<lpage>87</lpage>
<pub-id pub-id-type="doi">10.1016/j.devcel.2011.04.018</pub-id>
<pub-id pub-id-type="pmid">21664576</pub-id>
</element-citation>
</ref>
<ref id="CR12">
<label>12.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chuang</surname>
<given-names>PT</given-names>
</name>
<name>
<surname>McMahon</surname>
<given-names>AP</given-names>
</name>
</person-group>
<article-title>Vertebrate Hedgehog signalling modulated by induction of a Hedgehog-binding protein</article-title>
<source>Nature</source>
<year>1999</year>
<volume>397</volume>
<fpage>617</fpage>
<lpage>21</lpage>
<pub-id pub-id-type="doi">10.1038/17611</pub-id>
<pub-id pub-id-type="pmid">10050855</pub-id>
</element-citation>
</ref>
<ref id="CR13">
<label>13.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Kato</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Beachy</surname>
<given-names>PA</given-names>
</name>
</person-group>
<article-title>Gli2 trafficking links Hedgehog-dependent activation of Smoothened in the primary cilium to transcriptional activation in the nucleus</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2009</year>
<volume>106</volume>
<fpage>21666</fpage>
<lpage>71</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0912180106</pub-id>
<pub-id pub-id-type="pmid">19996169</pub-id>
</element-citation>
</ref>
<ref id="CR14">
<label>14.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goetz</surname>
<given-names>SC</given-names>
</name>
<name>
<surname>Anderson</surname>
<given-names>KV</given-names>
</name>
</person-group>
<article-title>The primary cilium: a signalling centre during vertebrate development</article-title>
<source>Nat Rev Genet</source>
<year>2010</year>
<volume>11</volume>
<fpage>331</fpage>
<lpage>44</lpage>
<pub-id pub-id-type="doi">10.1038/nrg2774</pub-id>
<pub-id pub-id-type="pmid">20395968</pub-id>
</element-citation>
</ref>
<ref id="CR15">
<label>15.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Akimaru</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Tanaka</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Maekawa</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Nakafuku</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Ishii</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Sonic Hedgehog-induced activation of the Gli1 promoter is mediated by GLI3</article-title>
<source>J Biol Chem</source>
<year>1999</year>
<volume>274</volume>
<fpage>8143</fpage>
<lpage>52</lpage>
<pub-id pub-id-type="doi">10.1074/jbc.274.12.8143</pub-id>
<pub-id pub-id-type="pmid">10075717</pub-id>
</element-citation>
</ref>
<ref id="CR16">
<label>16.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shin</surname>
<given-names>SH</given-names>
</name>
<name>
<surname>Kogerman</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Lindstrom</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Toftgard</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Biesecker</surname>
<given-names>LG</given-names>
</name>
</person-group>
<article-title>GLI3 mutations in human disorders mimic Drosophila cubitus interruptus protein functions and localization</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>1999</year>
<volume>96</volume>
<fpage>2880</fpage>
<lpage>4</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.96.6.2880</pub-id>
<pub-id pub-id-type="pmid">10077605</pub-id>
</element-citation>
</ref>
<ref id="CR17">
<label>17.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sasaki</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Nishizaki</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Hui</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Nakafuku</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Kondoh</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>Regulation of Gli2 and Gli3 activities by an amino-terminal repression domain: implication of Gli2 and Gli3 as primary mediators of Shh signaling</article-title>
<source>Development</source>
<year>1999</year>
<volume>126</volume>
<fpage>3915</fpage>
<lpage>24</lpage>
<pub-id pub-id-type="pmid">10433919</pub-id>
</element-citation>
</ref>
<ref id="CR18">
<label>18.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beachy</surname>
<given-names>PA</given-names>
</name>
<name>
<surname>Hymowitz</surname>
<given-names>SG</given-names>
</name>
<name>
<surname>Lazarus</surname>
<given-names>RA</given-names>
</name>
<name>
<surname>Leahy</surname>
<given-names>DJ</given-names>
</name>
<name>
<surname>Siebold</surname>
<given-names>C</given-names>
</name>
</person-group>
<article-title>Interactions between Hedgehog proteins and their binding partners come into view</article-title>
<source>Genes Dev</source>
<year>2010</year>
<volume>24</volume>
<fpage>2001</fpage>
<lpage>12</lpage>
<pub-id pub-id-type="doi">10.1101/gad.1951710</pub-id>
<pub-id pub-id-type="pmid">20844013</pub-id>
</element-citation>
</ref>
<ref id="CR19">
<label>19.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liem</surname>
<given-names>KF</given-names>
<suffix>Jr</suffix>
</name>
<name>
<surname>He</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Ocbina</surname>
<given-names>PJ</given-names>
</name>
<name>
<surname>Anderson</surname>
<given-names>KV</given-names>
</name>
</person-group>
<article-title>Mouse Kif7/Costal2 is a cilia-associated protein that regulates Sonic hedgehog signaling</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2009</year>
<volume>106</volume>
<fpage>13377</fpage>
<lpage>82</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.0906944106</pub-id>
<pub-id pub-id-type="pmid">19666503</pub-id>
</element-citation>
</ref>
<ref id="CR20">
<label>20.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aglyamova</surname>
<given-names>GV</given-names>
</name>
<name>
<surname>Agarwala</surname>
<given-names>S</given-names>
</name>
</person-group>
<article-title>Gene expression analysis of the hedgehog signaling cascade in the chick midbrain and spinal cord</article-title>
<source>Dev Dyn</source>
<year>2007</year>
<volume>236</volume>
<fpage>1363</fpage>
<lpage>73</lpage>
<pub-id pub-id-type="doi">10.1002/dvdy.21146</pub-id>
<pub-id pub-id-type="pmid">17436280</pub-id>
</element-citation>
</ref>
<ref id="CR21">
<label>21.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tenzen</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Allen</surname>
<given-names>BL</given-names>
</name>
<name>
<surname>Cole</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>JS</given-names>
</name>
<name>
<surname>Krauss</surname>
<given-names>RS</given-names>
</name>
<name>
<surname>McMahon</surname>
<given-names>AP</given-names>
</name>
</person-group>
<article-title>The cell surface membrane proteins Cdo and Boc are components and targets of the Hedgehog signaling pathway and feedback network in mice</article-title>
<source>Dev Cell</source>
<year>2006</year>
<volume>10</volume>
<fpage>647</fpage>
<lpage>56</lpage>
<pub-id pub-id-type="doi">10.1016/j.devcel.2006.04.004</pub-id>
<pub-id pub-id-type="pmid">16647304</pub-id>
</element-citation>
</ref>
<ref id="CR22">
<label>22.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Martinelli</surname>
<given-names>DC</given-names>
</name>
<name>
<surname>Fan</surname>
<given-names>CM</given-names>
</name>
</person-group>
<article-title>Gas1 extends the range of Hedgehog action by facilitating its signaling</article-title>
<source>Genes Dev</source>
<year>2007</year>
<volume>21</volume>
<fpage>1231</fpage>
<lpage>43</lpage>
<pub-id pub-id-type="doi">10.1101/gad.1546307</pub-id>
<pub-id pub-id-type="pmid">17504940</pub-id>
</element-citation>
</ref>
<ref id="CR23">
<label>23.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rahnama</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Toftgard</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Zaphiropoulos</surname>
<given-names>PG</given-names>
</name>
</person-group>
<article-title>Distinct roles of PTCH2 splice variants in Hedgehog signalling</article-title>
<source>Biochem J</source>
<year>2004</year>
<volume>378</volume>
<fpage>325</fpage>
<lpage>34</lpage>
<pub-id pub-id-type="doi">10.1042/bj20031200</pub-id>
<pub-id pub-id-type="pmid">14613484</pub-id>
</element-citation>
</ref>
<ref id="CR24">
<label>24.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Agren</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Kogerman</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Kleman</surname>
<given-names>MI</given-names>
</name>
<name>
<surname>Wessling</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Toftgard</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Expression of the PTCH1 tumor suppressor gene is regulated by alternative promoters and a single functional Gli-binding site</article-title>
<source>Gene</source>
<year>2004</year>
<volume>330</volume>
<fpage>101</fpage>
<lpage>14</lpage>
<pub-id pub-id-type="doi">10.1016/j.gene.2004.01.010</pub-id>
<pub-id pub-id-type="pmid">15087129</pub-id>
</element-citation>
</ref>
<ref id="CR25">
<label>25.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peterson</surname>
<given-names>KA</given-names>
</name>
<name>
<surname>Nishi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Vedenko</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Shokri</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Neural-specific Sox2 input and differential Gli-binding affinity provide context and positional information in Shh-directed neural patterning</article-title>
<source>Genes Dev</source>
<year>2012</year>
<volume>26</volume>
<fpage>2802</fpage>
<lpage>16</lpage>
<pub-id pub-id-type="doi">10.1101/gad.207142.112</pub-id>
<pub-id pub-id-type="pmid">23249739</pub-id>
</element-citation>
</ref>
<ref id="CR26">
<label>26.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vokes</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>H</given-names>
</name>
<name>
<surname>McCuine</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Tenzen</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Giles</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Zhong</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genomic characterization of Gli-activator targets in sonic hedgehog-mediated neural patterning</article-title>
<source>Development</source>
<year>2007</year>
<volume>134</volume>
<fpage>1977</fpage>
<lpage>89</lpage>
<pub-id pub-id-type="doi">10.1242/dev.001966</pub-id>
<pub-id pub-id-type="pmid">17442700</pub-id>
</element-citation>
</ref>
<ref id="CR27">
<label>27.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vokes</surname>
<given-names>SA</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>WH</given-names>
</name>
<name>
<surname>McMahon</surname>
<given-names>AP</given-names>
</name>
</person-group>
<article-title>A genome-scale analysis of the cis-regulatory circuitry underlying sonic hedgehog-mediated patterning of the mammalian limb</article-title>
<source>Genes Dev</source>
<year>2008</year>
<volume>22</volume>
<fpage>2651</fpage>
<lpage>63</lpage>
<pub-id pub-id-type="doi">10.1101/gad.1693008</pub-id>
<pub-id pub-id-type="pmid">18832070</pub-id>
</element-citation>
</ref>
<ref id="CR28">
<label>28.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holtz</surname>
<given-names>AM</given-names>
</name>
<name>
<surname>Peterson</surname>
<given-names>KA</given-names>
</name>
<name>
<surname>Nishi</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Morin</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>JY</given-names>
</name>
<name>
<surname>Charron</surname>
<given-names>F</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Essential role for ligand-dependent feedback antagonism of vertebrate hedgehog signaling by PTCH1, PTCH2 and HHIP1 during neural patterning</article-title>
<source>Development</source>
<year>2013</year>
<volume>140</volume>
<fpage>3423</fpage>
<lpage>34</lpage>
<pub-id pub-id-type="doi">10.1242/dev.095083</pub-id>
<pub-id pub-id-type="pmid">23900540</pub-id>
</element-citation>
</ref>
<ref id="CR29">
<label>29.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>EY</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Ouyang</surname>
<given-names>Z</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>W</given-names>
</name>
<name>
<surname>Vokes</surname>
<given-names>SA</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Hedgehog pathway-regulated gene networks in cerebellum development and tumorigenesis</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2010</year>
<volume>107</volume>
<fpage>9736</fpage>
<lpage>41</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.1004602107</pub-id>
<pub-id pub-id-type="pmid">20460306</pub-id>
</element-citation>
</ref>
<ref id="CR30">
<label>30.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fletez-Brant</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>D</given-names>
</name>
<name>
<surname>McCallion</surname>
<given-names>AS</given-names>
</name>
<name>
<surname>Beer</surname>
<given-names>MA</given-names>
</name>
</person-group>
<article-title>kmer-SVM: a web server for identifying predictive regulatory sequence features in genomic data sets</article-title>
<source>Nucleic Acids Res</source>
<year>2013</year>
<volume>41</volume>
<fpage>W544</fpage>
<lpage>56</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkt519</pub-id>
<pub-id pub-id-type="pmid">23771147</pub-id>
</element-citation>
</ref>
<ref id="CR31">
<label>31.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Karchin</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Beer</surname>
<given-names>MA</given-names>
</name>
</person-group>
<article-title>Discriminative prediction of mammalian enhancers from DNA sequence</article-title>
<source>Genome Res</source>
<year>2011</year>
<volume>21</volume>
<fpage>2167</fpage>
<lpage>80</lpage>
<pub-id pub-id-type="doi">10.1101/gr.121905.111</pub-id>
<pub-id pub-id-type="pmid">21875935</pub-id>
</element-citation>
</ref>
<ref id="CR32">
<label>32.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Spitz</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Furlong</surname>
<given-names>EE</given-names>
</name>
</person-group>
<article-title>Transcription factors: from enhancer binding to developmental control</article-title>
<source>Nat Rev Genet</source>
<year>2012</year>
<volume>13</volume>
<fpage>613</fpage>
<lpage>26</lpage>
<pub-id pub-id-type="doi">10.1038/nrg3207</pub-id>
<pub-id pub-id-type="pmid">22868264</pub-id>
</element-citation>
</ref>
<ref id="CR33">
<label>33.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Slattery</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Dantas Machado</surname>
<given-names>AC</given-names>
</name>
<name>
<surname>Gordan</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Rohs</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Absence of a simple code: how transcription factors read the genome</article-title>
<source>Trends Biochem Sci</source>
<year>2014</year>
<volume>39</volume>
<fpage>381</fpage>
<lpage>99</lpage>
<pub-id pub-id-type="doi">10.1016/j.tibs.2014.07.002</pub-id>
<pub-id pub-id-type="pmid">25129887</pub-id>
</element-citation>
</ref>
<ref id="CR34">
<label>34.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hallikas</surname>
<given-names>O</given-names>
</name>
<name>
<surname>Palin</surname>
<given-names>K</given-names>
</name>
<name>
<surname>Sinjushina</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Rautiainen</surname>
<given-names>R</given-names>
</name>
<name>
<surname>Partanen</surname>
<given-names>J</given-names>
</name>
<name>
<surname>Ukkonen</surname>
<given-names>E</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity</article-title>
<source>Cell</source>
<year>2006</year>
<volume>124</volume>
<fpage>1</fpage>
<lpage>13</lpage>
<pub-id pub-id-type="doi">10.1016/j.cell.2005.10.042</pub-id>
</element-citation>
</ref>
<ref id="CR35">
<label>35.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heinz</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Benner</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Spann</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Bertolino</surname>
<given-names>E</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>YC</given-names>
</name>
<name>
<surname>Laslo</surname>
<given-names>P</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities</article-title>
<source>Mol Cell</source>
<year>2010</year>
<volume>38</volume>
<fpage>576</fpage>
<lpage>89</lpage>
<pub-id pub-id-type="doi">10.1016/j.molcel.2010.05.004</pub-id>
<pub-id pub-id-type="pmid">20513432</pub-id>
</element-citation>
</ref>
<ref id="CR36">
<label>36.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bailey</surname>
<given-names>TL</given-names>
</name>
</person-group>
<article-title>DREME: motif discovery in transcription factor ChIP-seq data</article-title>
<source>Bioinformatics</source>
<year>2011</year>
<volume>27</volume>
<fpage>1653</fpage>
<lpage>9</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btr261</pub-id>
<pub-id pub-id-type="pmid">21543442</pub-id>
</element-citation>
</ref>
<ref id="CR37">
<label>37.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gupta</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Stamatoyannopoulos</surname>
<given-names>JA</given-names>
</name>
<name>
<surname>Bailey</surname>
<given-names>TL</given-names>
</name>
<name>
<surname>Noble</surname>
<given-names>WS</given-names>
</name>
</person-group>
<article-title>Quantifying similarity between motifs</article-title>
<source>Genome Biol</source>
<year>2007</year>
<volume>8</volume>
<fpage>R24</fpage>
<pub-id pub-id-type="doi">10.1186/gb-2007-8-2-r24</pub-id>
<pub-id pub-id-type="pmid">17324271</pub-id>
</element-citation>
</ref>
<ref id="CR38">
<label>38.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Crooks</surname>
<given-names>GE</given-names>
</name>
<name>
<surname>Hon</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Chandonia</surname>
<given-names>JM</given-names>
</name>
<name>
<surname>Brenner</surname>
<given-names>SE</given-names>
</name>
</person-group>
<article-title>WebLogo: a sequence logo generator</article-title>
<source>Genome Res</source>
<year>2004</year>
<volume>14</volume>
<fpage>1188</fpage>
<lpage>90</lpage>
<pub-id pub-id-type="doi">10.1101/gr.849004</pub-id>
<pub-id pub-id-type="pmid">15173120</pub-id>
</element-citation>
</ref>
<ref id="CR39">
<label>39.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Winklmayr</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Schmid</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Laner-Plamberger</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Kaser</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Aberger</surname>
<given-names>F</given-names>
</name>
<name>
<surname>Eichberger</surname>
<given-names>T</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Non-consensus GLI binding sites in Hedgehog target gene regulation</article-title>
<source>BMC Mol Biol</source>
<year>2010</year>
<volume>11</volume>
<fpage>2</fpage>
<pub-id pub-id-type="doi">10.1186/1471-2199-11-2</pub-id>
<pub-id pub-id-type="pmid">20070907</pub-id>
</element-citation>
</ref>
<ref id="CR40">
<label>40.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Machanick</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Bailey</surname>
<given-names>TL</given-names>
</name>
</person-group>
<article-title>MEME-ChIP: motif analysis of large DNA datasets</article-title>
<source>Bioinformatics</source>
<year>2011</year>
<volume>27</volume>
<fpage>1696</fpage>
<lpage>7</lpage>
<pub-id pub-id-type="doi">10.1093/bioinformatics/btr189</pub-id>
<pub-id pub-id-type="pmid">21486936</pub-id>
</element-citation>
</ref>
<ref id="CR41">
<label>41.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bailey</surname>
<given-names>TL</given-names>
</name>
<name>
<surname>Machanick</surname>
<given-names>P</given-names>
</name>
</person-group>
<article-title>Inferring direct DNA binding from ChIP-seq</article-title>
<source>Nucleic Acids Res</source>
<year>2012</year>
<volume>40</volume>
<fpage>e128</fpage>
<pub-id pub-id-type="doi">10.1093/nar/gks433</pub-id>
<pub-id pub-id-type="pmid">22610855</pub-id>
</element-citation>
</ref>
<ref id="CR42">
<label>42.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Creyghton</surname>
<given-names>MP</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>AW</given-names>
</name>
<name>
<surname>Welstead</surname>
<given-names>GG</given-names>
</name>
<name>
<surname>Kooistra</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Carey</surname>
<given-names>BW</given-names>
</name>
<name>
<surname>Steine</surname>
<given-names>EJ</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Histone H3K27ac separates active from poised enhancers and predicts developmental state</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>2010</year>
<volume>107</volume>
<fpage>21931</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.1016071107</pub-id>
<pub-id pub-id-type="pmid">21106759</pub-id>
</element-citation>
</ref>
<ref id="CR43">
<label>43.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heintzman</surname>
<given-names>ND</given-names>
</name>
<name>
<surname>Stuart</surname>
<given-names>RK</given-names>
</name>
<name>
<surname>Hon</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Ching</surname>
<given-names>CW</given-names>
</name>
<name>
<surname>Hawkins</surname>
<given-names>RD</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome</article-title>
<source>Nat Genet</source>
<year>2007</year>
<volume>39</volume>
<fpage>311</fpage>
<lpage>8</lpage>
<pub-id pub-id-type="doi">10.1038/ng1966</pub-id>
<pub-id pub-id-type="pmid">17277777</pub-id>
</element-citation>
</ref>
<ref id="CR44">
<label>44.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Visel</surname>
<given-names>A</given-names>
</name>
<name>
<surname>Thaller</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Eichele</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>GenePaint.org: an atlas of gene expression patterns in the mouse embryo</article-title>
<source>Nucleic Acids Res</source>
<year>2004</year>
<volume>32</volume>
<fpage>D552</fpage>
<lpage>6</lpage>
<pub-id pub-id-type="doi">10.1093/nar/gkh029</pub-id>
<pub-id pub-id-type="pmid">14681479</pub-id>
</element-citation>
</ref>
<ref id="CR45">
<label>45.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Diez-Roux</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Banfi</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Sultan</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Geffers</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Anand</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Rozado</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A high-resolution anatomical atlas of the transcriptome in the mouse embryo</article-title>
<source>PLoS Biol</source>
<year>2011</year>
<volume>9</volume>
<fpage>e1000582</fpage>
<pub-id pub-id-type="doi">10.1371/journal.pbio.1000582</pub-id>
<pub-id pub-id-type="pmid">21267068</pub-id>
</element-citation>
</ref>
<ref id="CR46">
<label>46.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>McLean</surname>
<given-names>CY</given-names>
</name>
<name>
<surname>Bristor</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Hiller</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Clarke</surname>
<given-names>SL</given-names>
</name>
<name>
<surname>Schaar</surname>
<given-names>BT</given-names>
</name>
<name>
<surname>Lowe</surname>
<given-names>CB</given-names>
</name>
<etal></etal>
</person-group>
<article-title>GREAT improves functional interpretation of cis-regulatory regions</article-title>
<source>Nat Biotechnol</source>
<year>2010</year>
<volume>28</volume>
<fpage>495</fpage>
<lpage>501</lpage>
<pub-id pub-id-type="doi">10.1038/nbt.1630</pub-id>
<pub-id pub-id-type="pmid">20436461</pub-id>
</element-citation>
</ref>
<ref id="CR47">
<label>47.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gerber</surname>
<given-names>AN</given-names>
</name>
<name>
<surname>Wilson</surname>
<given-names>CW</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>YJ</given-names>
</name>
<name>
<surname>Chuang</surname>
<given-names>PT</given-names>
</name>
</person-group>
<article-title>The hedgehog regulated oncogenes Gli1 and Gli2 block myoblast differentiation by inhibiting MyoD-mediated transcriptional activation</article-title>
<source>Oncogene</source>
<year>2007</year>
<volume>26</volume>
<fpage>1122</fpage>
<lpage>36</lpage>
<pub-id pub-id-type="doi">10.1038/sj.onc.1209891</pub-id>
<pub-id pub-id-type="pmid">16964293</pub-id>
</element-citation>
</ref>
<ref id="CR48">
<label>48.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Qiao</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Harris</surname>
<given-names>SE</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Oyajobi</surname>
<given-names>BO</given-names>
</name>
<name>
<surname>Mundy</surname>
<given-names>GR</given-names>
</name>
</person-group>
<article-title>The zinc finger transcription factor Gli2 mediates bone morphogenetic protein 2 expression in osteoblasts in response to hedgehog signaling</article-title>
<source>Mol Cell Biol</source>
<year>2006</year>
<volume>26</volume>
<fpage>6197</fpage>
<lpage>208</lpage>
<pub-id pub-id-type="doi">10.1128/MCB.02214-05</pub-id>
<pub-id pub-id-type="pmid">16880529</pub-id>
</element-citation>
</ref>
<ref id="CR49">
<label>49.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Levi</surname>
<given-names>G</given-names>
</name>
<name>
<surname>Topilko</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Schneider-Maunoury</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Lasagna</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Mantero</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Cancedda</surname>
<given-names>R</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Defective bone formation in Krox-20 mutant mice</article-title>
<source>Development</source>
<year>1996</year>
<volume>122</volume>
<fpage>113</fpage>
<lpage>20</lpage>
<pub-id pub-id-type="pmid">8565822</pub-id>
</element-citation>
</ref>
<ref id="CR50">
<label>50.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nonchev</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Maconochie</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Vesque</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Aparicio</surname>
<given-names>S</given-names>
</name>
<name>
<surname>Ariza-McNaughton</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Manzanares</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The conserved role of Krox-20 in directing Hox gene expression during vertebrate hindbrain segmentation</article-title>
<source>Proc Natl Acad Sci U S A</source>
<year>1996</year>
<volume>93</volume>
<fpage>9339</fpage>
<lpage>45</lpage>
<pub-id pub-id-type="doi">10.1073/pnas.93.18.9339</pub-id>
<pub-id pub-id-type="pmid">8790331</pub-id>
</element-citation>
</ref>
<ref id="CR51">
<label>51.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<collab>ENCODE Project Consortium</collab>
</person-group>
<article-title>An integrated encyclopedia of DNA elements in the human genome</article-title>
<source>Nature</source>
<year>2012</year>
<volume>489</volume>
<fpage>57</fpage>
<lpage>74</lpage>
<pub-id pub-id-type="doi">10.1038/nature11247</pub-id>
<pub-id pub-id-type="pmid">22955616</pub-id>
</element-citation>
</ref>
<ref id="CR52">
<label>52.</label>
<element-citation publication-type="journal">
<person-group person-group-type="author">
<name>
<surname>Santos</surname>
<given-names>N</given-names>
</name>
<name>
<surname>Reiter</surname>
<given-names>JF</given-names>
</name>
</person-group>
<article-title>Building it up and taking it down: the regulation of vertebrate ciliogenesis</article-title>
<source>Dev Dyn</source>
<year>2008</year>
<volume>237</volume>
<fpage>1972</fpage>
<lpage>81</lpage>
<pub-id pub-id-type="doi">10.1002/dvdy.21540</pub-id>
<pub-id pub-id-type="pmid">18435467</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>

Pour manipuler ce document sous Unix (Dilib)

EXPLOR_STEP=$WICRI_ROOT/Sante/explor/MersV1/Data/Pmc/Corpus

HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000289 | SxmlIndent | more

Ou

HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000289 | SxmlIndent | more

Pour mettre un lien sur cette page dans le réseau Wicri

{{Explor lien
   |wiki=    Sante
   |area=    MersV1
   |flux=    Pmc
   |étape=   Corpus
   |type=    RBID
   |clé=     PMC:4765071
   |texte=   Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway
}}

Pour générer des pages wiki

HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i   -Sk "pubmed:26912062" \
       | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd   \
       | NlmPubMed2Wicri -a MersV1
Wicri

This area was generated with Dilib version V0.6.33.
Data generation: Mon Apr 20 23:26:43 2020. Site generation: Sat Mar 27 09:06:09 2021