Conserved cis-regulatory regions in a large genomic landscape control SHH and BMP-regulated Gremlin1expression in mouse limb buds
© Zuniga et al.; licensee BioMed Central Ltd. 2012
Received: 28 March 2012
Accepted: 12 July 2012
Published: 13 August 2012
Mouse limb bud is a prime model to study the regulatory interactions that control vertebrate organogenesis. Major aspects of limb bud development are controlled by feedback loops that define a self-regulatory signalling system. The SHH/GREM1/AER-FGF feedback loop forms the core of this signalling system that operates between the posterior mesenchymal organiser and the ectodermal signalling centre. The BMP antagonist Gremlin1 (GREM1) is a critical node in this system, whose dynamic expression is controlled by BMP, SHH, and FGF signalling and key to normal progression of limb bud development. Previous analysis identified a distant cis-regulatory landscape within the neighbouring Formin1 (Fmn1) locus that is required for Grem1 expression, reminiscent of the genomic landscapes controlling HoxD and Shh expression in limb buds.
Three highly conserved regions (HMCO1-3) were identified within the previously defined critical genomic region and tested for their ability to regulate Grem1 expression in mouse limb buds. Using a combination of BAC and conventional transgenic approaches, a 9 kb region located ~70 kb downstream of the Grem1 transcription unit was identified. This region, termed Grem1 Regulatory Sequence 1 (GRS1), is able to recapitulate major aspects of Grem1 expression, as it drives expression of a LacZ reporter into the posterior and, to a lesser extent, in the distal-anterior mesenchyme. Crossing the GRS1 transgene into embryos with alterations in the SHH and BMP pathways established that GRS1 depends on SHH and is modulated by BMP signalling, i.e. integrates inputs from these pathways. Chromatin immunoprecipitation revealed interaction of endogenous GLI3 proteins with the core cis-regulatory elements in the GRS1 region. As GLI3 is a mediator of SHH signal transduction, these results indicated that SHH directly controls Grem1 expression through the GRS1 region. Finally, all cis-regulatory regions within the Grem1 genomic landscape locate to the DNAse I hypersensitive sites identified in this genomic region by the ENCODE consortium.
This study establishes that distant cis-regulatory regions scattered through a larger genomic landscape control the highly dynamic expression of Grem1, which is key to normal progression of mouse limb bud development.
During embryonic development, the spatio-temporal expression of transcriptional regulators and morphogenetic signals is tightly controlled. In fact, different types of congenital malformations are caused by mutations affecting specific cis-regulatory and non-coding genomic regions that alter the expression of genes or gene clusters in specific tissues (reviewed in ref. [1–3]). Genome-wide functional mapping revealed many of the large genomic landscapes that control the expression of developmental regulator genes in a global and tissue-specific manner . This study and others reveal that the cis-regulatory regions controlling gene expression in a particular tissue are often located several hundred kilobases (kb) or even megabases (Mb) up- or downstream of the transcriptional start sites in other loci, and may control the expression of several genes [5–7]. Recently, it has been shown how multiple cis-regulatory regions interact to control the expression of the 5’HoxD gene complex in the presumptive digit territory. The transcriptionally active part of the 5’HoxD gene cluster forms a so-called regulatory archipelago in which dispersed cis-regulatory elements co-operate to control gene expression in the distal limb bud by interacting with proximal promoters over large distances . In fact, regulatory landscapes with distant and dispersed cis-regulatory regions seem to be a recurring theme in the dynamic spatio-temporal regulation of genes that are essential for limb bud development (reviewed in ref. [2, 9, 10], see also below).
Vertebrate limb bud development is controlled by interactions between two main signalling centres, the apical ectodermal ridge (AER) and the zone of polarizing activity (ZPA) located in the posterior limb mesenchyme. These two signalling centres interact as part of a self-regulatory feedback signalling system involving several signalling pathways . In particular, ZPA cells produce the Sonic Hedgehog (SHH) signal, which together with GREMLIN1-mediated antagonism of Bone Morphogenetic Proteins (BMPs) in the posterior-distal limb bud mesenchyme propagates Fibroblast Growth Factor (FGF) signalling in the AER [11–14]. This SHH/GREM1/AER-FGF feedback-signalling loop promotes distal progression of limb bud outgrowth and formation of the autopod which gives rise to carpals and digit rays (reviewed in ref. ). SHH expression in the posterior limb bud mesenchyme is controlled by a cis-regulatory region located about 800 kb upstream of the transcriptional start site. Deletion of this ZPA regulatory region (ZRS) results in loss of distal limb structures , while point mutations within the ZRS result in anterior ectopic Shh expression and duplications of the thumb and/or anterior digits in different species [17, 18]. In cells actively transcribing Shh, the ZRS loops out of its chromosomal territory to the Shh transcription unit, which reveals how distant cis-regulatory elements control gene expression . Recently, several transcriptional regulators have been identified that control Shh expression by directly interacting with the distant ZRS region. These include HOX proteins, the bHLH transcription factor HAND2, and ETS transcription factors, providing a glimpse of the complex transcriptional regulation of Shh in the posterior limb bud mesenchyme [20–22].
SHH signalling is transduced by the GLI family of transcriptional regulators and inhibits the constitutive processing of the full-length GLI3 activator to its repressor form GLI3R (see e.g. ref. ). Of the three GLI transcription factors expressed during limb bud development, only Gli3 is essential on its own (see e.g. ref. ). The inactivation of Gli3 alters morphogenetic feedback signalling and results in formation of additional anterior digits [23, 25–27]. In particular, GLI3 restrains the proliferation of mesenchymal progenitors in the anterior limb bud mesenchyme and promotes initiation of digit ray chondrogenesis by directly repressing Grem1 expression during handplate formation .
We have previously shown that the expression of Grem1 in the limb bud mesenchyme is controlled by a large genomic landscape downstream of its transcription unit . In fact, several of the classical limb deformity mutations that disrupt distal limb bud development and formation of metanephric kidneys in the mouse are caused by deletions affecting this cis-regulatory landscape rather than directly altering the Grem1 gene products [7, 13, 28]. Molecular and genetic analysis in the mouse identified a 70 kb genomic region located downstream of the Grem1 transcription unit within the neighbouring Formin1 (Fmn1) gene that is required for Grem1 expression in the limb bud mesenchyme [7, 29]. Detailed analysis of this Grem1 genomic landscape revealed similarities with the global control region (GCR) that controls the expression of 5’HoxD genes in the limb bud mesenchyme , but did not reveal the structural nature and transacting factors and/or signalling pathways controlling these cis-regulatory regions.
To gain further insights into the Grem1 landscape, the potential cis-regulatory activities of the three highest evolutionarily conserved genomic regions within the 70 kb Grem1 critical region were analysed. In addition to its expression in the posterior mesenchyme, emphasis was also given to the coordinated distal-anterior expansion of Grem1 expression, which is key to orderly progression of limb bud development [14, 30]. Combining BAC with conventional transgenic approaches, we identified a 9 kb genomic region that is able to recapitulate major but not all aspects of the dynamic spatio-temporal Grem1 expression in the limb bud mesenchyme. This 9 kb region was termed Grem1 Regulatory Sequence 1 (GRS1), and contains a core sequence that is essential to express a LacZ transgene under control of a ß-globin minimal promoter in a Grem1-like pattern. The GRS1 region drives gene expression into the posterior and subsequently distal-anterior mesenchyme, i.e. reproduces aspects of the distal-anterior expansion of Grem1 expression during progression of limb bud development. By crossing GRS1 transgenic mice into different mutant contexts, we establish that the GRS1 region is controlled by inputs from both the SHH and BMP signalling pathways in limb buds. In addition, chromatin immunoprecipitation (ChIP) shows that GLI3 interacts with the conserved HMCO1 sequences in the GRS1 region. The functionally relevant cis-regulatory regions identified by the present and two previous studies [31, 32] map to the endogenous DNAse I hypersensitive sites within the Grem1 genomic landscape recently identified by the ENCODE consortium [33, 34]. This indicates that the Grem1 cis-regulatory landscape is composed of at least five active cis-regulatory regions that control the spatio-temporal expression of Grem1 in mouse embryos and limb buds.
Results and Discussion
Identification of conserved limb bud regulatory regions within the Grem1genomic landscape
Spatio-temporal activity of the GRS1region during mouse limb bud development
To analyze comparatively the activity of the GRS1 region with respect to the spatio-temporal regulation of Grem1 expression, transgenic mouse strains expressing the GRS1-ßglob-LacZ reporter construct were established (Figure 4). Seven transgenic founders were obtained. In three independent strains, LacZ was expressed in the posterior fore- and hindlimb bud mesenchyme in a pattern comparable to the transgenic founder embryos (Figure 3B and data not shown). Two of these three transgenic strains were analyzed in detail and used to study the spatio-temporal LacZ distribution (Figure 4B, C and Additional file 4). Initial analysis revealed that ß-galactosidase activity of the LacZ transgene was detected earlier in one strain and is significantly higher than in the other strain (Figure 4B, compare to Figure 4C). As levels did not change during subsequent generations (data not shown), this is likely due to differences in transgene copy number and/or integration site. Therefore, the copy number of both strains was determined by real-time qPCR analysis (Additional file 4). This analysis revealed that ten copies of the GRS1-ßglob-LacZ transgene were integrated into the genome of the strain expressing higher LacZ levels (GRS1 tg10 ), while only two copies were detected in the GRS1 tg2 strain (Additional file 4). In GRS1 tg10/+ embryos, ß-galactosidase activity is first detected at ~ E10.25 in the posterior forelimb bud mesenchyme and continuously increases until ~ E11.75 (Figure 4B). From ~ E11.0 onwards, scattered positive cells were detected in the anterior mesenchyme, and this anterior expression increased during distal limb bud outgrowth to form a crescent in the distal-anterior autopod (right panels, Figure 4B). By E11.75, ß-galactosidase activity was rather variable, such that the anterior crescent separated from the posterior domain in some forelimb buds. In GRS1 tg2/+ embryos, the spatio-temporal pattern is similar, but expression levels are significantly reduced due to the lower transgene copy number, resulting in detection of ß-galactosidase activity from only ~ E10.5 onwards (Figure 4C, arrowhead). In both strains, GRS1 drove expression of the ßglob-LacZ reporter specifically in the limb bud mesenchyme, despite some low ß-galactosidase activity detected in the developing eyes of GRS1 tg10/+ transgenic embryos (Additional file 4). This analysis reveals the robust nature of the GRS1 region, which functions in positive regulation of Grem1 expression in the limb bud mesenchyme. The posterior domain of ß-galactosidase activity at E10.5 (Figure 4B) is comparable to the Grem1 transcript distribution in wild-type limb buds at earlier stages (E9.5-E10.25, Figure 4D). This temporal delay is likely due to postponed transcriptional activation and/or up-regulation of the GRS1 transgene as the establishment of a posterior LacZ expression domain is also only apparent at E10.5 in GRS1 tg10/+ transgenic limb buds (Additional file 5). Therefore, additional cis-regulatory regions likely control the temporally correct early onset of Grem1 expression. Similar delays in the onset of LacZ reporter gene expression have been previously observed by analyzing the ZRS cis-regulatory region that controls Shh expression in mouse limb buds . During distal progression of limb bud development, endogenous Grem1 expression expands distal-anterior within the developing handplate and begins to fade by E11.75 (Figure 4D), due to termination of the SHH/GREM1/AER-FGF feedback loop and GLI3-mediated repression in the anterior limb bud mesenchyme [11, 26, 39, 40]. In contrast, ß-galactosidase activity remains high in the posterior limb bud mesenchyme of GRS1 transgenic embryos and the distal-anterior expansion of its expression is significantly delayed and only occurs as transcription of the endogenous Grem1 locus is starting to terminate (E11.5 onwards Figure 4B, C, compare to Figure 4D).
The GRS1 cis-regulatory region integrates inputs from the SHH and BMP signalling pathways
BMPs control Grem1 activation in the posterior limb bud mesenchyme and directly modulate its expression as part of the self-regulatory SHH/GREM1/AER-FGF feedback signalling system. Therefore, a conditional Smad4 loss-of-function allele  in combination with the Prx1-Cre recombinase strain  was used to inactivate Smad4 (Smad4 Δ/ΔM ) and thereby canonical BMP signal transduction in the limb bud mesenchyme. In Smad4 Δ/ΔM mutant limb buds, endogenous Grem1 expression is down-regulated but not lost, while Shh expression remains (Benazet et al., submitted). In Smad4 Δ/ΔM GRS1 tg10/+ forelimb buds, ß-galactosidase activity appears normal in the posterior mesenchyme, while distal-anterior expansion fails to occur (Figure 6C, compare to Figure 4B). This indicates that GRS1 activity is modulated by SMAD4-mediated signal transduction during the progression of limb bud development. These results indicate that also the distal-anterior expansion of Grem1 expression depends on BMP activity and agree with the observation that genetic lowering of Bmp4 results in a global reduction of Grem1 expression in the limb bud mesenchyme .
Finally, in Grem1-deficient embryos, BMP activity is increased due to reduced BMP antagonism, which results in up-regulation of non-functional Grem1 transcript . As BMP signal transduction modulates GRS1 activity during advanced limb bud development (Figure 6C), we also determined the potential influence of Grem1 deficiency on GRS1 activity. In Grem1 Δ/Δ GRS1 tg10/+ E10.5 limb buds ß-galactosidase activity is initially similar to GRS1 tg10/+ limb buds (Figure 6D, compare to Figure 4B). Subsequently, ß-galactosidase activity is increased and by E11.75 the entire Grem1 Δ/Δ GRS1 tg10/+ limb bud is positive (Figure 6D). These alterations of ß-galactosidase activity in Smad4 Δ/ΔM and Grem1 Δ/Δ forelimb buds reveal that GRS1 activity is extensively modulated by changes in BMP activity, in particular also the anterior expansion of its expression. Taken together, this genetic analysis (Figure 6) shows that GRS1 activity critically depends on SHH and is modulated extensively by BMP signal transduction. This analysis identifies the GRS1, located ~70 kb downstream of the Grem1 transcription unit as a cis-regulatory region that integrates inputs by both SHH and BMP signal transduction. However, as previous analysis provided good evidence that Grem1 is activated by mesenchymal BMP signalling upstream of establishing the SHH/FGF feedback loop , these studies point to the existence of additional unknown BMP response regions that activate Grem1 expression in the posterior limb bud mesenchyme.
Endogenous GLI3 proteins are part of the chromatin complexes interacting with specific parts of the GRS1region in mouse limb buds
Spatio-temporal Grem1 expression is regulated by the interaction of multiple cis-regulatory regions far downstream of the transcription unit
Previous genetic analysis showed that Grem1 transcription in limb buds is initiated by BMP signalling and up-regulated under the influence of SHH and AER-FGF signaling as limb bud development progresses . Finally, Grem1 expression is terminated concomitantly with the initiation of digit formation by high levels of FGF signal transduction and GLI3-mediated repression in the anterior mouse limb bud mesenchyme [26, 40]. This dynamic spatio-temporal regulation of Grem1 expression indicates that the activity of different cis-regulatory elements may change over time. The spatio-temporal activity of the GRS1 transgene (Figure 4B) indicates that it is primarily regulated by the SHH/GREM1/AER-FGF feedback signaling system as limb bud development progresses. In particular, its activity is first apparent when the SHH/GREM1/AER-FGF feedback is already established, and the expected termination does not occur, as ß-galactosidase activity remains in the posterior mesenchyme after the endogenous Grem1 transcripts have been down-regulated (Figure 4B, D). These results indicate that FGF-mediated termination of Grem1 expression and the underlying self-regulatory feedback signalling system  does not occur by FGF signal transduction impacting the GRS1 region. In addition, neither the 70 kb critical region of the Grem1 landscape (Grem1-LacZ, Figure 2A) nor transgenes derived from this region (Figures. 234 and ref. ) are able to accurately recapitulate the entire spatio-temporal distribution of Grem1 transcript in the limb bud mesenchyme. Therefore, either the interactions among these elements or additional as yet unknown cis-regulatory regions located outside the critical region must also participate in the regulation of Grem1 expression. Indeed, there is circumstantial evidence for the latter, as deletion of the GC-rich Fmn1 coding exon 9 (located ~200 kb downstream of the Grem1 transcription unit) results in non-complementation with a Grem1 null allele, causing a hypomorphic Grem1 limb skeletal phenotype, and significantly reduced Grem1 expression (ref.  and A.Z and R.Z., unpublished data).
In this study, we identify GRS1 as a distant cis-regulatory region that encompasses the highly conserved HMCO1 region and show that it controls important aspects of Grem1 expression in the posterior and distal-anterior limb bud mesenchyme. GRS1 activity depends critically on SHH signalling and endogenous GLI3 protein complexes interact with specific regions within HMCO1 in wild-type limb buds. In addition, the anterior expansion of GRS1 activity is also regulated by BMP signal transduction. These results together with previous studies reveal the large genomic architecture composed of several distant cis-regulatory regions that control the highly dynamic Grem1 expression in response to signalling inputs from the BMP, SHH and FGF pathways. It is likely that the interactions among these dispersed cis-regulatory regions provide the dynamic Grem1 expression in limb buds with the necessary robustness, but at the same time allow for evolutionary plasticity of its expression .
All the experiments were conducted in strict accordance with Swiss law following the 3R principles and the conduct defined by the Basel Declaration. All studies involving mice were classified as grade zero by the Animal Welfare and Ethics Commission of both cantons of Basel and Argovia, which implies minimal suffering.
The Shh loss-of-function allele , the Gli3 extra-toes-J allele , the Grem1 null allele  and the GRS1 tg/+ alleles (this study) were maintained in an NMRI background, while the Smad4 conditional allele (Smad4 flox , ) and the Prx1-Cre transgene  were maintained in a C57BL/6 J background.
BAC constructs were engineered using ET cloning as previously described . The deletions ΔHMCO1, ΔHMCO2 and ΔHMCO3 were induced in the Grem1-LacZ BAC clone using a Kanamycin resistance selection cassette. All primer sequences are available upon request. See Additional file 8 for the genomic coordinates of HMCO1-3.
An expression cassette was constructed in pKS-Bluescript using LacZeocin reporter under control of the minimal human ß-globin promoter sequences . Various combinations of HMCO sequences were inserted downstream into this cassette. The HMCO1, -2 and -3 core sequences were amplified by PCR directly from mouse genomic DNA (primers available on request). The 9 kb GRS1 region was initially subcloned as a 9.6 kb Bgl2 fragment from BAC RP23-113 H17 (chr2: 113,611,499-113,847,299 (mm10). BacPac Resources, Children’s Hospital, Oakland, USA). All constructs were linearized with Ksp1 before microinjection.
Generation of transgenic founder embryos and transgenic strains
BAC and ßglob-LacZ transgenic constructs were injected into the pronucleus of fertilized mouse eggs. Several founder embryos for each construct were scored for ß-galactosidase activity (Additional file 3). For the GRS1 ßglob-LacZ transgene, seven independent transgenic mouse strains were established by crossing founders into the NMRI background. Transgene copy numbers were determined by real time qPCR  using the Bio-Rad CFX96 Real-Time PCR System in combination with the iQ SYBR Green Supermix (Bio-Rad). For each mouse, 20 ng of genomic DNA were analysed in triplicate. The primers to amplify the DBH genomic region were used for normalization and non-transgenic littermate DNA was used as control. The normalized control levels were set to 2 as the primers also amplified both alleles of the endogenous HMCO1 regions. The 2-ΔΔCt formula was used to calculate and normalize transgene copy numbers. In addition to three different F1 males, three embryos from the F2 and three embryos from the F3 generation were analysed for the two transgenic strains shown to ascertain stable transmission of the transgenes. Primers sequences are available on request.
Detection of ß-galactosidase activity
Embryos were isolated in ice-cold PBS and staged according to somites numbers, then fixed in 1% formaldehyde, 0.2% glutaraldehyde, 0.02% NP40, 0.01% sodium deoxycholate in 1x PBS for 20–30 minutes at 4°C. Subsequently, embryos were washed three times in 1x PBS for 5 minutes at room temperature and incubated overnight at 37°C and in the dark in 1 mg/mL X-Gal, 0.25 mM K3Fe(CN6), 0.25 mM K4Fe(CN6), 0.01% NP-40, 0.4 mM MgCl2 and 1% sodium deoxycholate to detect ß-galactosidase activity, which colours cells blue. To stop the reaction, embryos were washed three times in 1x PBS for 5 minutes each at room temperature.
Chromatin Immunoprecipitation (ChIP)
Forelimbs and hindlimbs of 10 wild-type or Gli3 Xt/Xt embryos at E11.5 were dissected and processed for ChIP as described  using a polyclonal anti-GLI3 antibody (#2676, Genentech, ). To compute the level of enrichment of a given region, the Ct values of both input and ChIP samples were compared with those of a negative control amplicon located in the mouse ß-actin locus . All results (mean ± SD) were obtained by analysing three completely independent experiments per genotype. The significance of all differences was assessed using the two-tailed, non-parametric Mann–Whitney test. The coordinates of the relevant amplification are shown in Additional file 7.
The authors are grateful to A. Offinger and her staff for excellent animal care. L. d'Amato and Christina Torres are thanked for technical assistance. S. Scales (Genentech) provided the GLI3 antibodies. We thank group members for critical input on the manuscript. The ENCODE Group at University of Washington (J. Stamatoyannopoulos and coworkers) are acknowledged for the genome-wide DNase I hypersensitive site mapping data made available freely through the ENCODE Data Coordination Center at UCSC (release date on March 29, 2012). NM received a EU Marie Curie postdoctoral fellowship. This research is supported by grants of the Swiss National Science Foundation (grants no. 31003A_130803 to RZ and 310000_122558 to AZ), a EU reintegration grant (PERG-GA-2009-246576 to JLR) and the University of Basel. The authors declare no conflict of interest.
- Gordon CT, Tan TY, Benko S, Fitzpatrick D, Lyonnet S, et al: Long-range regulation at the SOX9 locus in development and disease. Journal of medical genetics. 2009, 46: 649-656. 10.1136/jmg.2009.068361.View ArticlePubMedGoogle Scholar
- Vandermeer JE, Ahituv N: cis-regulatory mutations are a genetic cause of human limb malformations. Dev Dyn. 2011, 240: 920-930. 10.1002/dvdy.22535.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuniga A, Probst S, Zeller R: The molecular basis of human congenital limb malformations. 2012, Wiley Interdisciplinary Reviews: Developmental BiologyGoogle Scholar
- Ruf S, Symmons O, Uslu VV, Dolle D, Hot C, et al: Large-scale analysis of the regulatory architecture of the mouse genome with a transposon-associated sensor. Nat Genet. 2011, 43: 379-386. 10.1038/ng.790.View ArticlePubMedGoogle Scholar
- Lettice LA, Horikoshi T, Heaney SJ, Van Baren MJ, Van Der Linde HC, et al: Disruption of a long-range cis-acting regulator for Shh causes preaxial polydactyly. Proc Natl Acad Sci U S A. 2002, 99: 7548-7553. 10.1073/pnas.112212199.PubMed CentralView ArticlePubMedGoogle Scholar
- Spitz F, Gonzalez F, Duboule D: A global control region defines a chromosomal regulatory landscape containing the HoxD cluster. Cell. 2003, 113: 405-417. 10.1016/S0092-8674(03)00310-6.View ArticlePubMedGoogle Scholar
- Zuniga A, Michos O, Spitz F, Haramis AP, Panman L, et al: Mouse limb deformity mutations disrupt a global control region within the large regulatory landscape required for Gremlin expression. Genes Dev. 2004, 18: 1553-1564. 10.1101/gad.299904.PubMed CentralView ArticlePubMedGoogle Scholar
- Montavon T, Soshnikova N, Mascrez B, Joye E, Thevenet L, et al: A regulatory archipelago controls Hox genes transcription in digits. Cell. 2011, 147: 1132-1145. 10.1016/j.cell.2011.10.023.View ArticlePubMedGoogle Scholar
- Spitz F: Control of vertebrate Hox clusters by remote and global cis-acting regulatory sequences. Adv Exp Med Biol. 2010, 689: 63-78. 10.1007/978-1-4419-6673-5_4.View ArticlePubMedGoogle Scholar
- Zeller R, Zuniga A: Shh and Gremlin1 chromosomal landscapes in development and disease. Curr Opin Genet Dev. 2007, 17: 428-434. 10.1016/j.gde.2007.07.006.View ArticlePubMedGoogle Scholar
- Benazet JD, Bischofberger M, Tiecke E, Goncalves A, Martin JF, et al: A self-regulatory system of interlinked signaling feedback loops controls mouse limb patterning. Science. 2009, 323: 1050-1053. 10.1126/science.1168755.View ArticlePubMedGoogle Scholar
- Capdevila J, Tsukui T, Rodriquez Esteban C, Zappavigna V, Izpisua Belmonte JC: Control of vertebrate limb outgrowth by the proximal factor Meis2 and distal antagonism of BMPs by Gremlin. Mol Cell. 1999, 4: 839-849. 10.1016/S1097-2765(00)80393-7.View ArticlePubMedGoogle Scholar
- Michos O, Panman L, Vintersten K, Beier K, Zeller R, et al: Gremlin-mediated BMP antagonism induces the epithelial-mesenchymal feedback signaling controlling metanephric kidney and limb organogenesis. Development. 2004, 131: 3401-3410. 10.1242/dev.01251.View ArticlePubMedGoogle Scholar
- Zuniga A, Haramis AP, McMahon AP, Zeller R: Signal relay by BMP antagonism controls the SHH/FGF4 feedback loop in vertebrate limb buds. Nature. 1999, 401: 598-602. 10.1038/44157.View ArticlePubMedGoogle Scholar
- Zeller R, Lopez-Rios J, Zuniga A: Vertebrate limb bud development: moving towards integrative analysis of organogenesis. Nat Rev Genet. 2009, 10: 845-858. 10.1038/nrg2681.View ArticlePubMedGoogle Scholar
- Sagai T, Hosoya M, Mizushina Y, Tamura M, Shiroishi T: Elimination of a long-range cis-regulatory module causes complete loss of limb-specific Shh expression and truncation of the mouse limb. Development. 2005, 132: 797-803. 10.1242/dev.01613.View ArticlePubMedGoogle Scholar
- Lettice LA, Hill AE, Devenney PS, Hill RE: Point mutations in a distant sonic hedgehog cis-regulator generate a variable regulatory output responsible for preaxial polydactyly. Hum Mol Genet. 2008, 17: 978-985.View ArticlePubMedGoogle Scholar
- Maas SA, Suzuki T, Fallon JF: Identification of spontaneous mutations within the long-range limb-specific Sonic hedgehog enhancer (ZRS) that alter Sonic hedgehog expression in the chicken limb mutants oligozeugodactyly and silkie breed. Dev Dyn. 2011, 240: 1212-1222. 10.1002/dvdy.22634.PubMed CentralView ArticlePubMedGoogle Scholar
- Amano T, Sagai T, Tanabe H, Mizushina Y, Nakazawa H, et al: Chromosomal dynamics at the Shh locus: limb bud-specific differential regulation of competence and active transcription. Dev Cell. 2009, 16: 47-57. 10.1016/j.devcel.2008.11.011.View ArticlePubMedGoogle Scholar
- Capellini TD, Di Giacomo G, Salsi V, Brendolan A, Ferretti E, et al: Pbx1/Pbx2 requirement for distal limb patterning is mediated by the hierarchical control of Hox gene spatial distribution and Shh expression. Development. 2006, 133: 2263-2273. 10.1242/dev.02395.View ArticlePubMedGoogle Scholar
- Galli A, Robay D, Osterwalder M, Bao X, Benazet JD, et al: Distinct roles of Hand2 in initiating polarity and posterior Shh expression during the onset of mouse limb bud development. PLoS Genet. 2010, 6: e1000901-10.1371/journal.pgen.1000901.PubMed CentralView ArticlePubMedGoogle Scholar
- Lettice LA, Williamson I, Wiltshire JH, Peluso S, Devenney PS, et al: Opposing Functions of the ETS Factor Family Define Shh Spatial Expression in Limb Buds and Underlie Polydactyly. Dev Cell. 2012, 22: 459-467. 10.1016/j.devcel.2011.12.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang B, Fallon JF, Beachy PA: Hedgehog-Regulated Processing of Gli3 Produces an Anterior/Posterior Repressor gradient in the Developing Vertebrate Limb. Cell. 2000, 100: 423-434. 10.1016/S0092-8674(00)80678-9.View ArticlePubMedGoogle Scholar
- Mo R, Freer AM, Zinyk DL, Crackower MA, Michaud J, et al: Specific and redundant functions of Gli2 and Gli3 zinc finger genes in skeletal patterning and development. Development. 1997, 124: 113-123.PubMedGoogle Scholar
- Litingtung Y, Dahn RD, Li Y, Fallon JF, Chiang C: Shh and Gli3 are dispensable for limb skeleton formation but regulate digit number and identity. Nature. 2002, 418: 979-983. 10.1038/nature01033.View ArticlePubMedGoogle Scholar
- Lopez-Rios J, Speziale D, Robay D, Scotti M, Osterwalder M, et al: GLI3 Constrains Digit Number by Controlling Both Progenitor Proliferation and BMP-Dependent Exit to Chondrogenesis. Developmental Cell. 2012, 22: 837-848. 10.1016/j.devcel.2012.01.006.PubMed CentralView ArticlePubMedGoogle Scholar
- te Welscher P, Zuniga A, Kuijper S, Drenth T, Goedemans HJ, et al: Progression of Vertebrate Limb Development through SHH-Mediated Counteraction of GLI3. Science. 2002, 298: 827-830. 10.1126/science.1075620.View ArticlePubMedGoogle Scholar
- Khokha MK, Hsu D, Brunet LJ, Dionne MS, Harland RM: Gremlin is the BMP antagonist required for maintenance of Shh and Fgf signals during limb patterning. Nat Genet. 2003, 34: 303-307. 10.1038/ng1178.View ArticlePubMedGoogle Scholar
- Zhou F, Leder P, Zuniga A, Dettenhofer M: Formin1 disruption confers oligodactylism and alters Bmp signaling. Hum Mol Genet. 2009, 18: 2472-2482. 10.1093/hmg/ddp185.PubMed CentralView ArticlePubMedGoogle Scholar
- Panman L, Galli A, Lagarde N, Michos O, Soete G, et al: Differential regulation of gene expression in the digit forming area of the mouse limb bud by SHH and gremlin 1/FGF-mediated epithelial-mesenchymal signalling. Development. 2006, 133: 3419-3428. 10.1242/dev.02529.View ArticlePubMedGoogle Scholar
- Soshnikova N, Montavon T, Leleu M, Galjart N, Duboule D: Functional analysis of CTCF during mammalian limb development. Dev Cell. 2010, 19: 819-830. 10.1016/j.devcel.2010.11.009.View ArticlePubMedGoogle Scholar
- Vokes SA, Ji H, Wong WH, McMahon AP: A genome-scale analysis of the cis-regulatory circuitry underlying sonic hedgehog-mediated patterning of the mammalian limb. Genes Dev. 2008, 22: 2651-2663. 10.1101/gad.1693008.PubMed CentralView ArticlePubMedGoogle Scholar
- Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447: 799-816. 10.1038/nature05874.View ArticlePubMedGoogle Scholar
- Myers RM, Stamatoyannopoulos J, Snyder M, Dunham I, Hardison RC, et al: A user's guide to the encyclopedia of DNA elements (ENCODE). PLoS biology. 2011, 9: e1001046-10.1371/journal.pbio.1001046.View ArticleGoogle Scholar
- Pearl EJ, Barker D, Day RC, Beck CW: Identification of genes associated with regenerative success of Xenopus laevis hindlimbs. BMC Dev Biol. 2008, 8: 66-10.1186/1471-213X-8-66.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicoli S, Gilardelli CN, Pozzoli O, Presta M, Cotelli F: Regulated expression pattern of gremlin during zebrafish development. Gene expression patterns : GEP. 2005, 5: 539-544. 10.1016/j.modgep.2004.11.001.View ArticlePubMedGoogle Scholar
- Boisvert CA, Mark-Kurik E, Ahlberg PE: The pectoral fin of Panderichthys and the origin of digits. Nature. 2008, 456: 636-638. 10.1038/nature07339.View ArticlePubMedGoogle Scholar
- Lettice LA, Heaney SJ, Purdie LA, Li L, de Beer P, et al: A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly. Hum Mol Genet. 2003, 12: 1725-1735. 10.1093/hmg/ddg180.View ArticlePubMedGoogle Scholar
- Scherz PJ, Harfe BD, McMahon AP, Tabin CJ: The limb bud Shh-Fgf feedback loop is terminated by expansion of former ZPA cells. Science. 2004, 305: 396-399. 10.1126/science.1096966.View ArticlePubMedGoogle Scholar
- Verheyden JM, Sun X: An Fgf/Gremlin inhibitory feedback loop triggers termination of limb bud outgrowth. Nature. 2008, 454: 638-641. 10.1038/nature07085.PubMed CentralView ArticlePubMedGoogle Scholar
- Bandyopadhyay A, Tsuji K, Cox K, Harfe BD, Rosen V, et al: Genetic Analysis of the Roles of BMP2, BMP4, and BMP7 in Limb Patterning and Skeletogenesis. PLoS Genet. 2006, 2: e216-10.1371/journal.pgen.0020216.PubMed CentralView ArticlePubMedGoogle Scholar
- Tzchori I, Day TF, Carolan PJ, Zhao Y, Wassif CA, et al: LIM homeobox transcription factors integrate signaling events that control three-dimensional limb patterning and growth. Development. 2009, 136: 1375-1385. 10.1242/dev.026476.PubMed CentralView ArticlePubMedGoogle Scholar
- Nissim S, Hasso SM, Fallon JF, Tabin CJ: Regulation of Gremlin expression in the posterior limb bud. Dev Biol. 2006, 299: 12-21. 10.1016/j.ydbio.2006.05.026.View ArticlePubMedGoogle Scholar
- te Welscher P, Fernandez-Teran M, Ros MA, Zeller R: Mutual genetic antagonism involving GLI3 and dHAND prepatterns the vertebrate limb bud mesenchyme prior to SHH signaling. Genes Dev. 2002, 16: 421-426. 10.1101/gad.219202.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang X, Li C, Herrera PL, Deng CX: Generation of Smad4/Dpc4 conditional knockout mice. Genesis. 2002, 32: 80-81. 10.1002/gene.10029.View ArticlePubMedGoogle Scholar
- Logan M, Martin JF, Nagy A, Lobe C, Olson EN, et al: Expression of Cre recombinase in the developing mouse limb bud driven by a Prxl enhancer. Genesis. 2002, 33: 77-80. 10.1002/gene.10092.View ArticlePubMedGoogle Scholar
- Wen X, Lai CK, Evangelista M, Hongo JA, de Sauvage FJ, et al: Kinetics of hedgehog-dependent full-length Gli3 accumulation in primary cilia and subsequent degradation. Molecular and cellular biology. 2010, 30: 1910-1922. 10.1128/MCB.01089-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Ong CT, Corces VG: Enhancer function: new insights into the regulation of tissue-specific gene expression. Nature reviews Genetics. 2011, 12: 283-293.PubMed CentralView ArticlePubMedGoogle Scholar
- Zeller R: The temporal dynamics of vertebrate limb development, teratogenesis and evolution. Curr Opin Genet Dev. 2010, 20: 384-390. 10.1016/j.gde.2010.04.014.View ArticlePubMedGoogle Scholar
- St-Jacques B, Dassule HR, Karavanova I, Botchkarev VA, Li J, et al: Sonic hedgehog signaling is essential for hair development. Curr Biol. 1998, 8: 1058-1068. 10.1016/S0960-9822(98)70443-9.View ArticlePubMedGoogle Scholar
- Maynard TM, Jain MD, Balmer CW, LaMantia AS: High-resolution mapping of the Gli3 mutation extra-toes reveals a 51.5-kb deletion. Mamm Genome. 2002, 13: 58-61. 10.1007/s00335-001-2115-X.View ArticlePubMedGoogle Scholar
- Morgan BA, Conlon FL, Manzanares M, Millar JB, Kanuga N, et al: Transposon tools for recombinant DNA manipulation: characterization of transcriptional regulators from yeast, Xenopus, and mouse. Proc Natl Acad Sci U S A. 1996, 93: 2801-2806. 10.1073/pnas.93.7.2801.PubMed CentralView ArticlePubMedGoogle Scholar
- Tesson L, Heslan JM, Menoret S, Anegon I: Rapid and accurate determination of zygosity in transgenic animals by real-time quantitative PCR. Transgenic research. 2002, 11: 43-48. 10.1023/A:1013928600442.View ArticlePubMedGoogle Scholar
- Ovcharenko I, Nobrega MA, Loots GG, Stubbs L: ECR Browser: a tool for visualizing and accessing data from comparisons of multiple vertebrate genomes. Nucleic acids research. 2004, 32: W280-W286. 10.1093/nar/gkh355.PubMed CentralView ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG: ClustalW and ClustalX version 2. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.