spacer gif spacer gif spacer gif spacer gif spacer gif
QUICK SEARCH:   [advanced]


spacer gif
     Home     Help     Feedback     Subscriptions     Archive     Search     Table of Contents    

First published online 25 June 2008
doi: 10.1242/dev.017384

Development 135, 2543-2553 (2008)
Published by The Company of Biologists 2008
This Article
Right arrow Summary Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Supplementary Material
Right arrow All Versions of this Article:
dev.017384v1
135/15/2543    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Menke, D. B.
Right arrow Articles by Kingsley, D. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Menke, D. B.
Right arrow Articles by Kingsley, D. M.
Social Bookmarking
Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?

Dual hindlimb control elements in the Tbx4 gene and region-specific control of bone size in vertebrate limbs

Douglas B. Menke, Catherine Guenther and David M. Kingsley*

Howard Hughes Medical Institute and Department of Developmental Biology, Stanford University, Stanford, CA 94305-5329, USA.

* Author for correspondence (e-mail: kingsley{at}cmgm.stanford.edu)

Accepted 19 May 2008


   SUMMARY
TOP
 SUMMARY
INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The Tbx4 transcription factor is crucial for normal hindlimb and vascular development, yet little is known about how its highly conserved expression patterns are generated. We have used comparative genomics and functional scanning in transgenic mice to identify a dispersed group of enhancers controlling Tbx4 expression in different tissues. Two independent enhancers control hindlimb expression, one located upstream and one downstream of the Tbx4 coding exons. These two enhancers, hindlimb enhancer A and hindlimb enhancer B (HLEA and HLEB), differ in their primary sequence, in their precise patterns of activity within the hindlimb, and in their degree of sequence conservation across animals. HLEB is highly conserved from fish to mammals. Although Tbx4 expression and hindlimb development occur at different axial levels in fish and mammals, HLEB cloned from either fish or mouse is capable of driving expression at the appropriate position of hindlimb development in mouse embryos. HLEA is highly conserved only in mammals. Deletion of HLEA from the endogenous mouse locus reduces expression of Tbx4 in the hindlimb during embryogenesis, bypasses the embryonic lethality of Tbx4-null mutations, and produces viable, fertile mice with characteristic changes in the size of bones in the hindlimb but not the forelimb. We speculate that dual hindlimb enhancers provide a flexible genomic mechanism for altering the strength and location of Tbx4 expression during normal development, making it possible to separately modify the size of forelimb and hindlimb bones during vertebrate evolution.

Key words: Tbx4, Hindlimb, cis-regulatory element, Enhancer


   INTRODUCTION
TOP
 SUMMARY
 INTRODUCTION
MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Tetrapod limbs show great morphological diversity (Flower, 1876Go; Mariani and Martin, 2003Go). Although a common pattern of skeletal structures can be seen in different species, the position of limbs, and the size, morphology and number of bones within limbs are extensively modified in organisms adapted to running, hopping, flying or swimming. Differences between forelimbs and hindlimbs of a single species can be just as striking as the differences exhibited among highly diverged species. For example, digit bones show striking elongation in flying animals owing to a local increase in skeletal growth in forelimbs but not hindlimbs. Conversely, pelvic structures are lost or greatly reduced in marine mammals, an alteration that affects the hindlimbs without disrupting normal forelimb development.

Numerous signaling and transcription factor pathways have now been identified that participate in limb development (Capdevila and Izpisua Belmonte, 2001Go; Mariani and Martin, 2003Go; Tickle, 2003Go). Most of these factors are expressed in both forelimbs and hindlimbs. However, particular members of the T-box family of DNA-binding proteins are expressed specifically in forelimbs or hindlimbs of many different animals (Chapman et al., 1996Go; Gibson-Brown et al., 1998Go; Logan et al., 1998Go; Tamura et al., 1999Go; Ruvinsky et al., 2000Go; Takabatake et al., 2000Go; Tanaka et al., 2002Go). Tbx5 is expressed in the forelimb field prior to forelimb outgrowth and continues to be expressed specifically in the forelimb mesenchyme as limb development proceeds. In humans, haploinsufficiency for TBX5 results in reduced and malformed forelimbs, as well as heart defects in Holt-Oram syndrome (Basson et al., 1997Go; Li et al., 1997Go; Agarwal et al., 2003Go; Rallis et al., 2003Go). Likewise, loss of Tbx5 function in mice and zebrafish completely prevents outgrowth of forelimbs and pectoral fins, respectively (Ahn et al., 2002Go; Garrity et al., 2002Go).

A related gene, Tbx4, is specifically expressed in the lateral plate mesoderm of the early hindlimb field and in the developing hindlimb bud (Chapman et al., 1996Go). Haploinsufficiency for TBX4 in humans results in the hindlimb-specific defects of small-patella syndrome (SPS), including incomplete ossification of the pelvis and small or reduced kneecaps (Bongers et al., 2004Go). Complete loss of mouse Tbx4 function also disrupts hindlimb development, although detailed studies have been difficult because of mid-gestation lethality due to vascular defects (Naiche and Papaioannou, 2003Go; Naiche and Papaioannou, 2007Go).

Although Tbx4 and Tbx5 are among the earliest markers of the prospective forelimb and hindlimb fields (Gibson-Brown et al., 1996Go), it is not yet clear how their forelimb or hindlimb specificity is established, or what causes the genes to turn on at different axial levels in animals whose forelimbs and hindlimbs form at different positions along the body. Previous studies suggest that the homeodomain transcription factor Pitx1 may act upstream of Tbx4. Pitx1 also shows hindlimb-specific expression, and can induce Tbx4 expression when overexpressed in chick forelimbs (Logan and Tabin, 1999Go). The related paralog, Pitx2, has similar trans-activating properties and is initially co-expressed with Pitx1 in the lateral plate mesoderm of the early hindlimb field. However, Tbx4 expression is reduced, but not eliminated, in Pitx1-/- or Pitx1-/- Pitx2+/- embryos (Lanctot et al., 1999Go; Szeto et al., 1999Go; Marcil et al., 2003Go). Although Tbx4 levels in Pitx1-/- Pitx2-/- embryos have not been examined owing to early embryonic lethality, Tbx4 expression is maintained in Pitx1-/- Pitx2+/- embryos at stages of hindlimb development at which Pitx2 is not expressed (Marcil et al., 2003Go). Together, these results suggest that Pitx1 and Pitx2 are important upstream regulators of Tbx4, and that additional factors might also contribute to hindlimb-specific expression of Tbx4.

Here we report a large-scale enhancer survey of the Tbx4 locus. Our transgenic and comparative studies show that two different cis- regulatory elements contribute to hindlimb expression. Targeted knockout of one of these elements bypasses the lethality seen in Tbx4-knockout mice, and reveals an essential role for this enhancer in regulating bone size in developing hindlimbs.


   MATERIALS AND METHODS
TOP
 SUMMARY
 INTRODUCTION
 MATERIALS AND METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
Transgenic mice
Transgenic mice were generated by pronuclear injections of FVB or C57BL6/CBA F2 embryos (Stanford Transgenic Research Facility and Xenogen Biosciences). BAC and plasmid DNAs were purified as previously described (DiLeone et al., 2000Go). Multiple independent embryos were scored for each transgene and only reproducible patterns were scored as significant (see Table S1 in the supplementary material).

Hsp68LacZ transgenes
Potential enhancer regions were amplified using primers (Table 1) containing NotI restriction sites and cloned into the NotI site of p5'-Not-HspLacZ (DiLeone et al., 1998Go). In cases where multiple tandem copies were desired, products were amplified with primers containing XbaI and SpeI sites, cloned into a modified pBS KS+ plasmid with two NotI sites flanking an XbaI site, excised with NotI and cloned into p5'-Not-HspLacZ.


View this table:
[in this window]
[in a new window]

  		Table 1. Hsp68LacZ transgene primers

 
BAC modifications
BACs RP24-376P4, RP24-84E15 and RP23-136J3 were identified by electronically searching mouse BAC end sequences (Zhao et al., 2001Go), and were modified by homologous recombination in bacteria (Lee et al., 2001Go). A Tbx4 targeting cassette was generated by amplifying 5' and 3' homology arms and cloning them into pIPTGfTet, which contains an IRES-βGeo cassette and an FRT-flanked tetracycline resistance gene (TetR) (Chandler et al., 2007Go). 5' homology primers: Tbx4-IBG1, 5'-AGCTTGCTAGCACATGTGGTCAGAGAGCATC-3' and Tbx4-IBG2, 5'-ATGAACTCGAGTAATGGTGGACTCGGCTCCT-3'. 3' homology primers: Tbx4-IBG3, 5'-TTGACACTAGTCCCTGGGTGGGATAGAAGTA-3' and Tbx4-IBG4, 5'-AGTTATGCGGCCGCAAGAAAGGCCCTCCCTGTG-3'. Modified BACs contained an IRES-βGeo cassette 245 bp downstream of the Tbx4 stop codon. The FRT-flanked TetR was removed by transient expression of FLPe recombinase.

BAC deletions were generated by amplifying TetR with primers possessing homology to regions flanking particular enhancers. HLEA: HLEA-DelF, 5'-CCATGGAGCTCCAGGCTGCTTGGGGGAGGAGGCCGAAGAGAGGGAACCCAAGATCTATGATTCCCTTTGTC-3' and HLEA-DelR, 5'-GGCTAGAACTGTGACTTCTCCAAGAGTCAACAGGCCCTAGACTGGACTCTTAGAGAATAGGAACTTCAAGCT-3'. HLEB: HLEB-DelF, 5'-AATGGCATTTTCCCTCACAAAGTCACACAAATAAGTTGGCTCAGGAGACTAGATCTATGATTCCCTTTGTC-3' and HLEB-DelR, 5'-CGAGCATGGGTTGTGTGCTCTGCATTCGCATGAGACATAACACCGTGTACTAGAGAATAGGAACTTCAAGCT-3'. TetR was subsequently removed from modified BACs by negative selection on fusaric acid after retargeting using homology arms identical to those used to create the initial TetR targeting cassette, but which lacked intervening TetR sequence (Yang et al., 1997Go). Retargeting fragments were generated from overlapping primers: HLEA-TetRDelF, 5'-CCATGGAGCTCCAGGCTGCTTGGGGGAGGAGGCCGAAGAGAGGGAACCCAAGAGTCCAGT-3' and HLEA-TetRDelR, 5'-GGCTAGAACTGTGACTTCTCCAAGAGTCAACAGGCCCTAGACTGGACTCTTGGGTTCCCT-3'; HLEB-TetRDelF, 5'-CCATGGAGCTCCAGGCTGCTTGGGGGAGGAGGCCGAAGAGAGGGAACCCAAGAGTCCAGT-3' and HLEB-TetRDelR, 5'-GGCTAGAACTGTGACTTCTCCAAGAGTCAACAGGCCCTAGACTGGACTCTTGGGTTCCCT-3'.

lacZ staining and in situ hybridization
Whole-mount lacZ staining and in situ hybridization were performed as described (Wilkinson, 1992Go; DiLeone et al., 1998Go). The Tbx4 riboprobe was generated by amplifying a 959 bp fragment from mouse embryonic hindlimb cDNA, cloning into pCR4Blunt-TOPO vector (Invitrogen), digesting with NotI, and transcribing with T3 RNA polymerase. Tbx4 primers: Tbx4-P1F, 5'-CAAGGAGTATCCCGTGATCT-3' and Tbx4-P2R, 5'-CACATTCTGAAATACCTTTCCATG-3'.

Quantification of lacZ staining
Forelimbs and hindlimbs of E12.5 embryos carrying transgenes pHLEA-768 and pPitx1Mut were measured for area of lacZ staining using ImageJ (Abramoff et al., 2004Go). Ratios of hindlimb to forelimb staining for different constructs were compared using ANCOVA in S-Plus (6.0) (Insightful Corporation).

Comparative sequence analysis
Sequences from Mus musculus, Homo sapiens, Felis catus, Canis familiaris, Bos taurus, Dasypus novemcinctus, Loxodonta africana, Monodelphis domestica, Ornithorhynchus anatinus, Gallus gallus, Anolis carolinensis, Danio rerio, Gasterosteus aculeatus and Callorhinchus milii were obtained from NCBI, aligned with Multi-LAGAN and analyzed with VISTA (Brudno et al., 2003Go; Frazer et al., 2004Go).

Generation of hindlimb enhancer-knockout mice
An HLEA targeting vector was generated by amplifying 5' and 3' homology arms from mouse strain 129P2 (5' arm: genome assembly mm9; Chr11: 85,687,031-85,689,553; 3' arm: Chr11: 85,690,623-85,696,175) and cloning into pLoxPNT flanking a floxed Neo cassette (Shalaby et al., 1995Go). The HLEA sequence with a 3' flanking loxP site was then inserted between the floxed Neo and the 3' homology arm to create targeting vector pHwHLEAKO. Homology arms were sequenced to verify the absence of mutations.

Homologous targeting was performed on ES cell line CGR8.8 (strain background 129P2) by the Stanford Transgenic Facility using standard protocols (Joyner, 1993Go; Townley et al., 1997Go). ES cell colonies were screened for correct targeting at the 5' end by PCR. PCR-positive clones were verified by Southern blot using a 3' probe. Chimeras were bred to C57BL/6J females resulting in germline transmission from two ES cell clones. Heterozygous mice were bred to mice carrying a CMV:Cre transgene (mouse strain DBA/2Tg-2.6I, maintained on a C57BL/6J background, provided by David Anderson, California Institute of Technology, Pasadena, CA) to generate the Tbx4{Delta}HLEA allele. Mice were maintained on a mixed background of 129P2 and C57BL/6J.


Figure 1
View larger version (95K):
[in this window]
[in a new window]

  		Fig. 1. A BAC enhancer scan across the mouse Tbx4 region. (A) Schematic of genomic region encompassing Tbx4. Arrows denote the direction of transcription of Tbx4 and flanking genes. Gray bars show BAC locations with the inserted IRES-βGeo cassette in blue. Mapped cis-regulatory intervals are shaded: yellow, hindlimb interval I; green, lung interval; pink, hindlimb interval II. (B,C) Whole-mount in situ hybridization for Tbx4 mRNA in E12.5 mouse embryos. (D-I) Whole-mount lacZ staining of transgenic embryos carrying RP24-376P4, RP24-84E15 or RP23-136J3. Dotted lines denote boundaries of limbs, umbilical cord (uc) and genital tubercle (gt). lu, lung.

 
Mice were genotyped using a three-primer multiplex PCR that produces a wild-type product of 164 bp and Tbx4{Delta}HLEA product of 263 bp: HLEA1F, 5'-AACCTGGCTGAAGACTCCTG-3'; HLEA2F, 5'-AGGACATGTTTCTGAACGAGC-3'; and HLEA1R, 5'-GCTGTCCGAGGAATGCCATG-3'.

Animals were handled using protocols approved by the Stanford University Animal Care and Use Committee.

Allele-specific expression
Levels of Tbx4 expression from wild-type and Tbx4{Delta}HLEA alleles were measured in E11.5 F1 hybrid embryos generated by crossing DBA/2J males and 129P2 or Tbx4{Delta}HLEA females. The second coding exon of Tbx4 was amplified from lung and hindlimb RNA using RT-PCR with forward primer 5'-CGGAGCAGACCATCGAGAA-3' and reverse primer (5'-biotinylated, HPLC-purified) 5'-TGCCTGCCTCGTGGAACTT-3'. Pyrosequencing was performed using primer 5'-CCATCGAGAACATCAA-3' to analyze the relative levels of mouse SNP rs29454648 derived from DBA/2J and 129P2 Tbx4 alleles (EpigenDx, Worcester, MA). Statistical analyses were performed by arcsine transforming the percentage of DBA/2J and 129P2 transcripts from each tissue sample and performing a two-tailed t-test assuming unequal variance.

Skeletal preparations and measurements
Skeletal preparations of 8-week-old male mice were prepared with Alcian Blue and Alizarin Red staining as described (Lufkin et al., 1992Go) and measured under a dissecting microscope with an eyepiece reticule (n=14 wild-type, n=15 Tbx4{Delta}HLEA heterozygous and n=15 Tbx4{Delta}HLEA homozygous animals). Body lengths and forelimb and hindlimb bone sizes were compared using a two-tailed t-test assuming unequal variation. No statistically significant differences were seen for body length or forelimb bone sizes between wild-type and mutant animals, and normalization of hindlimb measurements to body length or forelimb bone size produced similar results.


   RESULTS
TOP
 SUMMARY
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
A BAC scan across the Tbx4 locus localizes cis-regulatory elements
To screen large regions upstream and downstream of Tbx4 for limb control elements, we identified three mouse BACs that span a 344 kb region (Fig. 1A). Each BAC contains the entire coding region of Tbx4 plus varying amounts of flanking sequence. We inserted an IRES-βGeo reporter into the Tbx4 3' UTR of each clone, and generated transient transgenic embryos that were stained for β-galactosidase activity with X-Gal (lacZ staining).

The medial and distal BACs, RP24-84E15 and RP23-136J3, generated essentially identical patterns of lacZ staining in embryos at day 12.5 (E12.5) of development. Robust expression was observed in hindlimbs and in the lung, anterior genital tubercle and umbilical cord (Fig. 1F-I). This pattern closely resembles the endogenous pattern of Tbx4 expression, suggesting that these BACs contain most of the major Tbx4 cis-regulatory elements required at E12.5 (Fig. 1B,C) (Chapman et al., 1996Go). By contrast, the proximal BAC, RP24-376P4, exhibited staining in the lung and anterior genital tubercle, but minimal staining in the umbilical cord and incomplete staining in the hindlimb with no expression observed in the distal half of the autopod (Fig. 1D,E).

The similar lung and genital tubercle expression seen with all three BACs suggests that lung and genital tubercle regulatory elements are located in an 86 kb region shared by all three clones. We refer to this region as the lung interval (Fig. 1A). The weak hindlimb expression driven by the proximal BAC suggests that some hindlimb control elements map to the region contained within this BAC (hindlimb interval I, Fig. 1A). The much more complete hindlimb and umbilical cord staining observed with the medial and distal BACs suggests that a distal shared 59 kb region contains additional hindlimb elements, as well as umbilical cord elements (hindlimb interval II, Fig. 1A).

Isolation and characterization of Tbx4 enhancers
Previous studies suggest that evolutionarily conserved sequences often correspond to tissue-specific enhancers (Fortini and Rubin, 1990Go; Mortlock et al., 2003Go; Woolfe et al., 2005Go; Pennacchio et al., 2006Go). Alignment of mouse and human Tbx4 region sequences revealed numerous evolutionarily conserved regions (ECRs) with 70% or higher sequence identity over at least 100 bp (Fig. 2A). We cloned several of these ECR-containing regions in front of an Hsp68 minimal promoter and a lacZ reporter, and tested whether individual fragments were capable of driving reporter gene activity at consistent locations in transgenic mice.


Figure 2
View larger version (70K):
[in this window]
[in a new window]

  		Fig. 2. Comparative genomics reveals locations of candidate cis-regulatory elements. (A) Sequence comparison of mouse with human, opossum, chicken and zebrafish Tbx4. Colored bars (top) denote location of cis-regulatory intervals; blue boxes indicate position of Tbx2 and Tbx4 exons. All regions with >70% identity over a 100 bp window are colored: blue, exons; red, conserved intergenic sequence; pink, conserved intronic sequences). Block-shaded portions of plots indicate positions of confirmed enhancers (see B): yellow, hindlimb enhancer A (HLEA); green, lung enhancer; pink, hindlimb enhancer B (HLEB); gray, umbilical cord. (B) Genomic regions assayed for enhancer activity. Black bars represent fragments that did not drive consistent expression patterns at E12.5. (C-G) Whole-mount lacZ staining of E12.5 transgenic mouse embryos illustrating expression patterns seen with the different constructs. tr, trachea; lu, lung; uc, umbilical cord.

 
A 6.3 kb fragment within the hindlimb I interval drove reproducible lacZ staining in the hindlimbs of transient transgenic E12.5 embryos (Fig. 2B, transgene pDBM2, and data not shown). A 1073 bp subregion of this sequence showed high sequence identity between placental mammals and opossum, and exhibited identical enhancer activity to the full-length fragment (Fig. 2B,C, transgene pDBM7). We designated this 1073 bp fragment hindlimb enhancer A (HLEA). HLEA produced strong expression in the proximal portion of the hindlimb, but was generally excluded from the distal half of the autopod. Sections through hindlimbs showed that HLEA exhibits strong activity throughout the proximal limb mesenchyme at E12.5 (see Fig. S1 in the supplementary material). Lower and more-restricted enhancer activity was also observed in the forelimb (Fig. 2C).

Within the lung interval defined by BAC scanning, the third intron of Tbx4 contains the highest concentration of ECRs, including sequences with high conservation to chicken (Fig. 2A). A 5.5 kb segment from this region drove consistent expression in the developing lung and trachea of transgenic mouse embryos at E12.5 (Fig. 2B,D, transgene pDBM3).

The 59 kb hindlimb II interval showed numerous regions conserved between mouse and human. We designed four transgenes that together enabled us to test 94% of these conserved sequences for enhancer activity (Fig. 2A,B). Only one of the four transgenes, pDBM40, drove consistent reporter expression at E12.5, with expression clearly observed in the hindlimb, umbilical cord and anterior genital tubercle (Fig. 2E and data not shown). The hindlimb staining exhibited a clear posterior bias, which was evident both in whole-mount and sectioned material (Fig. 2E and see Fig. S1 in the supplementary material). Reporter expression was also reproducibly observed in scattered cells throughout the torso and facial region, but this pattern did not resemble endogenous Tbx4 expression.

A 3.5 kb subregion of the original 9.2 kb pDBM40 fragment drove hindlimb, genital tubercle and umbilical cord expression but had lost the scattered torso/facial activity (Fig. 2B,F, transgene pDBM45, and data not shown). To further define the hindlimb control sequences, we subdivided the 3.5 kb fragment into two parts. The 3' half (pDBM50) retained the genital tubercle and hindlimb enhancer activity, albeit at somewhat reduced levels (Fig. 2B,G). Neither half reproducibly drove umbilical cord expression, suggesting that both regions are required for umbilical enhancer activity. The 3' portion of pDBM50 contains a region with high sequence conservation from mice to cartilaginous fish (Fig. 3A). We cloned a 654 bp region from mouse that encompassed all the sequences conserved in fish and placed four tandem copies upstream of Hsp68LacZ (Fig. 2B, pDBM5). This construct drove extremely intense lacZ expression throughout the hindlimb and genital tubercle (Fig. 3B and data not shown). We designated this 654 bp fragment hindlimb enhancer B (HLEB).


Figure 3
View larger version (60K):
[in this window]
[in a new window]

  		Fig. 3. HLEB is conserved from mammals to cartilaginous fish. (A) VISTA plot (Frazer et al., 2004Go) comparing the 654 bp mouse HLEB sequence against sequences from human, opossum, platypus, chicken, lizard, frog, teleost fish (stickleback) and cartilaginous fish. Percentage sequence identity [y-axis, ranging from 40% to 100% (lower to upper thick line in each panel)] is shown in a sliding window across the region (x-axis, bp). Pink shading indicates regions of 20 bp or more that are ≥65% identical. (B,C) Whole-mount lacZ staining of E12.5 transgenic mouse embryos. (B) Hsp68LacZ transgene pDBM5 with four copies of mouse HLEB. (C) Hsp68LacZ transgene pDBM20 with four tandem copies of stickleback HLEB.

 
Because of the extreme conservation of HLEB, we also tested whether the orthologous sequence from fish exhibited similar enhancer activity to mouse HLEB when injected into mouse embryos. We cloned four copies of HLEB from a marine threespine stickleback, a fish with large spine-like pelvic fins (Bell and Foster, 1994Go), upstream of Hsp68LacZ. Stickleback HLEB is 56% identical to mouse HLEB over a 300 bp region. Mouse embryos that carried the stickleback HLEB transgene exhibited strong lacZ expression in the posterior portion of the hindlimb at E12.5 (Fig. 3C, pDBM20). However, enhancer activity dropped off dramatically at the boundary of the autopod. Of fifteen independent F0 transgenic embryos, eleven completely lacked expression in the autopod. In the remaining four embryos, lacZ staining extended into the autopod, but this staining was weaker than that in the proximal hindlimb (Fig. 3C and data not shown). At E9.5 and E10.5, expression was also posteriorly restricted, but to a lesser degree than at E12.5 (see Fig. S2 in the supplementary material). In addition, expression at E10.5 extended into the distal portion of the limb bud. Thus, fish HLEB also drives hindlimb expression, but with less autopod expression than seen with the mouse element.

Developmental time course of HLEA and HLEB
We examined the activity of mouse HLEA and HLEB at additional stages by generating stable transgenic lines carrying either the 1073 bp HLEA element or four copies of the 654 bp HLEB element upstream of Hsp68LacZ. Both enhancers were active in the hindlimb field before the onset of hindlimb bud outgrowth and continued to drive strong expression at E10.5 and E11.5 in patterns that strongly resembled endogenous Tbx4 mRNA expression (Fig. 4). HLEA demonstrated fairly uniform activity throughout the hindlimb, before losing activity in the distal half of the autopod at E12.5 (Fig. 2C and Fig. 4). HLEA demonstrated varying levels of activity in the forelimb, but in all embryos forelimb staining was less intense and less extensive than in hindlimbs. By comparison, HLEB was much more hindlimb-specific. Both enhancers were active in hindlimbs during late embryogenesis. However, by E17.5, HLEA activity was restricted to knee and ankle bones (see Fig. S3 in the supplementary material). By E16.5, HLEB activity was strongest in knees, including the distal femur, proximal tibia and patellar ligament (see Fig. S3 in the supplementary material).

HLEA and HLEB are essential for Tbx4 hindlimb expression
We tested the relative importance HLEA and HLEB by examining the effects of deleting them, individually and in combination, from BAC RP24-84E15 (Fig. 5A-C). Deletion of HLEA dramatically reduced hindlimb expression, although expression was still seen in the posterior region of the hindlimb (Fig. 5D,E). Expression in the lungs, genital tubercle and umbilical cord was unaltered. The remaining hindlimb expression seen from the HLEA-deleted BAC (which still retains HLEB) appeared weaker than that seen when HLEB itself was cloned immediately upstream of a heterologous promoter (Fig. 2F). This difference might be due to the fact that in the BAC clone, HLEB is located over 75 kb downstream of the Tbx4 promoter. Deletion of only HLEB from RP24-84E15 reduced expression in the most proximal and distal regions of the hindlimb, but the loss of expression was less dramatic than that seen with the deletion of HLEA (Fig. 5F,G).

When we deleted both HLEA and HLEB from RP24-84E15, hindlimb lacZ staining was undetectable or barely detectable in transgenic embryos (Fig. 5H,I). Expression in lung, genital tubercle and umbilical cord was not impaired, showing that HLEA and HLEB are required primarily for limb expression. We conclude that both HLEA and HLEB are required for robust hindlimb expression. Moreover, the patterns observed when HLEA and HLEB are deleted individually are not complementary, suggesting that these enhancers act synergistically rather than in a purely additive manner.

Mutation of a putative Pitx1 binding site reduces HLEA activity
Pitx1 expression in chicken forelimbs can induce Tbx4 (Logan and Tabin, 1999Go). We therefore examined mouse HLEA for Pitx1/Pitx2 consensus binding sites [TAA(T/G)C(C/T) (Lamonerie et al., 1996Go; Tremblay et al., 2000Go)]. Three sites were found in a 768 bp subregion of the original 1073 bp HLEA that was sufficient to drive strong hindlimb expression in transgenic embryos (Fig. 6C, transgene HLEA-768). One predicted binding site is conserved in 17 of 18 placental mammals; a second is perfectly conserved in placental mammals, marsupials and monotremes; and a third is found only in mouse (Fig. 6A and see Fig. S4 in the supplementary material).


Figure 4
View larger version (134K):
[in this window]
[in a new window]

  		Fig. 4. Developmental time course of HLEA and HLEB activity compared with the endogenous expression of Tbx4. (A-F) Whole-mount lacZ staining of mouse embryos from stable transgenic lines carrying HLEA transgene pDBM7 (A-C) or HLEB transgene pDBM5 (D-F) show that these enhancers are active in the early hindlimb field (A,D) and developing hindlimb bud (B,C,E,F). (G-I) Tbx4 mRNA as detected by whole-mount in situ hybridization shows similar patterns of expression. Embryonic stages as indicated.

 


Figure 5
View larger version (87K):
[in this window]
[in a new window]

  		Fig. 5. HLEA and HLEB are both required for robust hindlimb expression of Tbx4. (A) Schematic of mouse Tbx4 genomic region. BAC transgene RP24-84E15 and derivatives are represented by gray bars. Yellow ovals, HLEA; green circles, lung enhancer; gray circles, umbilical cord element; pink ovals, HLEB; blue box, IRES-βGeo. (B-I) Lateral (B,D,F,H) and ventral-posterior (C,E,G,I) views of whole-mount lacZ staining of E12.5 transgenic mouse embryos carrying BAC transgenes RP24-84E15 (B,C), {Delta}HLEA (D,E), {Delta}HLEB (F,G) or {Delta}HLEA{Delta}HLEB (H,I).

 
We mutated three bp in the single perfectly conserved Pitx1 binding site, creating a sequence that should be incapable of binding Pitx1/Pitx2 (from TAATCC to GACTAC, Fig. 6, transgene Pitx1Mut) (Lamonerie et al., 1996Go; Tremblay et al., 2000Go; Lamolet et al., 2001Go). The Pitx1Mut enhancer still drove lacZ expression in transgenic embryos. However, expression in hindlimbs appeared lower than that seen with the wild-type enhancer (Fig. 6C,D). Quantitation of lacZ staining areas confirmed a significant reduction in hindlimb relative to forelimb staining in Pitx1Mut-LacZ compared with control transgenic embryos (P=0.012) (Fig. 6B). Thus, the perfectly conserved Pitx1 binding site is required for full hindlimb activity, but additional sequences must also contribute to the remaining forelimb and hindlimb activity of HLEA. These additional sequences might include additional Pitx1 binding sites, or sites for additional factors that also influence limb expression.


Figure 6
View larger version (81K):
[in this window]
[in a new window]

  		Fig. 6. Mutation of a conserved putative Pitx1 binding site in HLEA results in reduced hindlimb enhancer activity. (A) VISTA plot comparing 768 bp subregion of mouse HLEA against human, cat, dog, cow, armadillo, elephant, opossum and platypus sequences. Horizontal panels show sequence identity across the region, ranging from 50-100%. Pink shading indicates regions of 20 bp or more that are ≥65% identical. Arrowheads indicate positions of putative Pitx1 binding sites, including one that is perfectly conserved (red arrowhead). The accompanying sequence alignment (beneath) shows perfectly conserved bases (asterisks), the sequence of the wild-type Pitx1 site (red), and the sequence of the mutated site in Pitx1Mut. (B) Plot of the area of lacZ staining in forelimbs and hindlimbs of transgenic mouse embryos carrying wild-type (HLEA-768) or mutant (Pitx1Mut) enhancers. Each data point represents an individual transgenic embryo. (C,D) Whole-mount lacZ staining of E12.5 transgenic mouse embryos carrying wild-type (C) and mutant (D) enhancers.

 
Knockout of HLEA results in reduced Tbx4 expression
To rigorously test the function of a putative Tbx4 hindlimb control region, we used homologous recombination to precisely excise HLEA from the endogenous mouse Tbx4 locus. We focused on HLEA because removal of this element from the BAC had the largest effect on hindlimb expression. We anticipated that removal of HLEA would generate a hypomorphic regulatory allele that would bypass the embryonic lethality observed with Tbx4-null embryos (Naiche and Papaioannou, 2003Go), and thereby allow us to both determine the functional role of HLEA and observe the effects of reduced Tbx4 expression on normal hindlimb development. We created mice carrying a floxed HLEA allele via homologous recombination in ES cells (Fig. 7A,B). We bred heterozygous mice carrying the floxed allele to a Cre deleter strain to generate the Tbx4{Delta}HLEA allele in which HLEA has been completely deleted from the locus, leaving a single loxP site in its place (Fig. 7A,C). We refer to this allele as {Delta}HLEA.

Animals heterozygous or homozygous for the {Delta}HLEA allele were viable and fertile and occurred at normal Mendelian ratios in crosses between heterozygous carriers (50 +/+, 98 +/{Delta}HLEA, 58 {Delta}HLEA/{Delta}HLEA). At E9.5, whole-mount in situ hybridization revealed no obvious differences in Tbx4 expression between {Delta}HLEA/{Delta}HLEA and wild-type embryos (see Fig. S5 in the supplementary material). However, by E10.5, {Delta}HLEA/{Delta}HLEA embryos had significantly less Tbx4 hindlimb bud expression (Fig. 7E-H). This altered expression was most apparent in the anterior portion of hindlimb buds, and was also observed at E11.5 and E12.5 (Fig. 7I,J and data not shown). In wild-type embryos, low levels of Tbx4 mRNA were routinely found in forelimb buds. This forelimb expression was reduced in {Delta}HLEA/{Delta}HLEA E10.5 and E11.5 forelimb buds, but by E12.5 forelimb expression was less affected (Fig. 7E,F and data not shown). As expected, umbilical cord and genital tubercle expression did not appear to be affected.

We used an allele-specific expression assay to more precisely quantify the impact of deleting HLEA on Tbx4 expression. The Tbx4 transcript contains SNPs that differ between 129P2 and DBA mouse strains. Since the {Delta}HLEA allele was generated on a 129P2 genetic background, we generated mice that carried one DBA allele and either the wild-type 129P2 or {Delta}HLEA 129P2 Tbx4 allele. We performed RT-PCR on E11.5 hindlimbs and lungs and quantified the relative levels of the DBA and 129P2 allele-specific SNPs using pyrosequencing. In both hindlimbs and lungs, the wild-type DBA and 129P2 alleles were expressed at equivalent levels (Fig. 7D). By contrast, the {Delta}HLEA allele was expressed at levels comparable to the DBA allele in the lungs but at ~3-fold lower levels in the hindlimbs (129/DBA, n=18; 129{Delta}HLEA/DBA, n=16; P=2.3x10-21).

Previous work showed that conditional inactivation of Tbx4 at E9.5 results in perturbed expression of Fgf8, Hand2 and Alx4 and in smaller hindlimbs (Naiche and Papaioannou, 2007Go). No differences in expression of these genes were seen in {Delta}HLEA mutant embryos by in situ hybridization (see Fig. S6 in the supplementary material).

Skeletal alterations in {Delta}HLEA mice
Adult animals homozygous for the HLEA deletion showed small, but significant, reductions in the size of multiple bones in the hindlimb, including the pelvis, femur, tibia and patella (Fig. 8A). The hindlimb bones of heterozygous adults were also affected and were intermediate in size between the wild type and homozygous mutants (see Table S2 in the supplementary material). In the feet of homozygotes, anterior digits were more severely reduced than were posterior digits, and distal-most elements were more reduced than were proximal elements (Fig. 8C and see Table S3 in the supplementary material).

The peroneal process on the head of the fibula was missing in all mutants, and there were fusions between tarsal bones of the ankle. Kneecaps were also displaced in 33% of animals (see Fig. S7 in the supplementary material). By contrast, no significant change was seen in average body length of wild-type and mutant mice, and length and width measurements for forelimb bones were similar, confirming that HLEA acts as a quantitative element for hindlimb but not forelimb growth during normal development (Fig. 8B).


Figure 7
View larger version (55K):
[in this window]
[in a new window]

  		Fig. 7. Targeted knockout of HLEA reduces expression of Tbx4 in hindlimbs. (A) Schematic showing homologous targeting of the Tbx4 locus in mouse ES cells to create Tbx4floxNeoHLEA, a floxed HLEA knock-in allele, and Tbx4{Delta}HLEA, where the floxed HLEA allele has been deleted using Cre recombinase. Gray bars, first and second exons of Tbx4; yellow ovals, HLEA; triangles, loxP sites; red arrows, genotyping primers; Neo, neomycin resistance cassette. (B) Southern blot of ES cell clones digested with MfeI shows a wild-type (Wt) band of 24 kb and a targeted Tbx4floxNeoHLEA band of 9.1 kb with an external 3' probe (see A). (C) Mice carrying the Tbx4floxNeoHLEA allele were crossed to a Cre deleter strain to generate the Tbx4{Delta}HLEA allele. Genotyping by PCR confirmed the presence of heterozygous and homozygous animals. (D) The relative expression levels of wild-type 129P2 and {Delta}HLEA 129P2 Tbx4 alleles were compared with a wild-type DBA Tbx4 allele in E11.5 hindlimbs and lungs. Black bars, expression ratio of wild-type 129P2 to DBA; gray bars, ratio of {Delta}HLEA 129P2 to DBA. (E-J) Whole-mount in situ hybridization for Tbx4 mRNA in wild-type and homozygous {Delta}HLEA knockout embryos. Lateral views of E10.5 embryos (E,F) and dorsal views of E10.5 hindlimb buds (G,H) and E12.5 hindlimbs (I,J). Arrowheads indicate anterior side of hindlimb buds.

 
To determine whether the size reduction observed in adults is apparent during embryogenesis, we also measured the widths of forelimbs and hindlimbs of E11.5 wild-type and homozygous mutant embryos. When normalized to crown-rump length, we found that hindlimbs of homozygous mutants were ~7% narrower (+/+, n=12; {Delta}HLEA/{Delta}HLEA, n=11; P=4.1x10-5), whereas forelimbs of the wild type and mutants did not differ significantly (P=0.20). These results show that small, but significant, effects on hindlimb growth begin at the limb bud stage, before overt differentiation and ossification of skeletal elements.


   DISCUSSION
TOP
 SUMMARY
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
Our analysis of cis-regulatory sequences controlling Tbx4 shows a dispersed group of enhancers spanning upstream, downstream and intronic regions of the locus, including two distinct regions controlling hindlimb expression. The hindlimb enhancers differ in their sequence, precise patterns of activity within limbs, position relative to the Tbx4 promoter region, and their degree of conservation. Our BAC deletion analysis shows that both enhancers are required for robust hindlimb expression, and suggests that these two enhancers might act synergistically to achieve this expression.

Conservation of cis-regulatory sequences from fish to mammals
Tbx4 expression in developing hindlimbs/fins has been faithfully conserved among most vertebrates with paired appendages. Are the factors regulating this expression also well conserved? HLEB is conserved at the sequence level from mice to cartilaginous fish. This indicates that despite dramatic alterations in the location and morphology of limbs in different animals, at least some mechanisms regulating hindlimb expression of Tbx4 have been maintained across ~500 million years (Blair and Hedges, 2005Go).

In order to determine the extent of HLEB functional conservation, we examined enhancer activity of the HLEB from a fish in the context of a developing mouse embryo. At E9.5, the HLEB from the threespine stickleback drives patterns of hindlimb-specific expression that closely resemble those of mouse HLEB (compare Fig. 4D with Fig. S2A in the supplementary material). The functional conservation of this enhancer is noteworthy given that sticklebacks, as well as other spiny-rayed teleosts, have pelvic fins that develop at a far more anterior position than those of other vertebrates (Nelson, 1994Go). Previous studies have shown Tbx4 is expressed at the level of the developing hindfin in sticklebacks and so has also undergone an anterior shift (Cole et al., 2003Go). The different axial levels of Tbx4 expression and hindlimb development could have occurred through alterations in the cis-regulatory sequences of Tbx4 in different animals or through shifts in the expression domains of upstream trans-acting factors. Since the stickleback HLEB drives expression in the hindlimb of mouse embryos, both the fish and mammal enhancers presumably respond to similar trans-acting factors. These results suggest that shifts in Tbx4 expression and hindlimb position in different animals occur by changes in the expression domains of shared upstream factors, rather than by modification of Tbx4 cis-acting sequences.

Whereas at early stages of mouse hindlimb development the stickleback HLEB drives expression extending into the distal-most portion of the mouse hindlimb bud, by E12.5 the activity of the stickleback enhancer is sharply decreased at the boundary of the autopod (Fig. 3C and see Fig. S2 in the supplementary material). Thus, during early phases of mouse hindlimb bud outgrowth, stickleback HLEB behaves much more like mouse HLEB than at later phases of hindlimb development. This functional difference in enhancer activity is interesting when considered in the context of fish fin and tetrapod limb development and morphology. The digits of mice and of other tetrapods are believed to be novel distal structures without clear homologs in fish, whereas sticklebacks have a spine-like pelvic fin formed from the fusion of enlarged dermal fin rays. The stickleback enhancer might lack crucial binding sites for autopod-specific or distally restricted trans-acting factors that are required to maintain activity in the distal portion of the mouse hindlimb at E12.5. This absence of crucial sites in stickleback HLEB could be due to a gain of transcription factor binding sites in the HLEB of tetrapods. Alternatively, sticklebacks, and possibly other teleosts, might have lost binding sites in HLEB that are important in maintaining distal enhancer activity. Given that the HLEB of elephant shark, a cartilaginous fish, demonstrates significantly better conservation to mammals than does the HLEB of sticklebacks and other teleosts (our unpublished observations), we favor the second possibility (Fig. 3). This loss in teleosts would also be consistent with recent studies that have shown that teleosts have lost the later phases of Hox gene expression in their fins, whereas this expression has been retained in more-basal bony fish and cartilaginous fish (Davis et al., 2007Go; Freitas et al., 2007Go).


Figure 8
View larger version (28K):
[in this window]
[in a new window]

  		Fig. 8. Homozygous {Delta}HLEA mice show reduced size of hindlimb bones. (A,B) Bar charts showing mean sizes of wild-type (black bars) and homozygous {Delta}HLEA mutant (gray bars) bones in hindlimbs (A) and forelimbs (B) of adult mice. Error bars indicate s.e.m. Red brackets in diagrams to right show positions of measurements. (C) Metatarsals and phalanges of digit rays 1 (red), 3 (blue) and 5 (gray) were measured in the feet of the wild type and homozygous mutants. Shown are percent decreases in the sizes of individual bones, or in the combined length of digit rays, in mutants compared with wild type. *P<0.05, **P<0.005, ***P<5x10-5.

 
A second hindlimb enhancer highly conserved in mammals
In addition to the HLEB enhancer shared with fish, mammals have a second enhancer highly conserved from placental mammals to monotremes. HLEA is located over 86 kb from HLEB, and shows no obvious sequence homology with HLEB, suggesting that HLEA did not arise by simple duplication of an ancestral element. It is possible that some relationship exists between the two enhancers that cannot be detected by simple sequence alignment. Previous studies have shown that the sequence of enhancer elements can change to the point that global alignment is lost, even though the divergent enhancers still drive similar expression patterns and may respond to similar upstream factors (Ludwig et al., 2005Go; Fisher et al., 2006Go). Several features of Tbx4 regulation, however, suggest that the two enhancers have somewhat different functions. First, the detailed expression patterns driven by HLEA and HLEB clearly differ, and this difference becomes more pronounced as expression is traced to later stages of development (see Fig. S3 in the supplementary material). Second, the patterns of residual expression following deletion of HLEA or HLEB from a BAC transgene also differ, confirming that some regions of the hindlimb are more dependent on one or the other of the two enhancers for Tbx4 expression (Fig. 5). Third, the morphological effects of deleting HLEA are not uniform across the hindlimb. Anterior digits are more affected than posterior digits, a morphological effect that parallels the loss of Tbx4 expression in the most anterior portion of the autopod when HLEA is deleted (compare Fig. 7E-J with Fig. 8C).

Deletion of HLEA and HLEB from a BAC transgene results in loss of Tbx4 expression in distinct hindlimb domains. Anterior expression is completely lost from hindlimbs when HLEA is deleted. By contrast, deletion of HLEB results in reduced expression in proximal and distal regions (Fig. 5). These patterns of residual expression are not complementary and suggest that these two enhancers interact synergistically. Consistent with our BAC results, deletion of HLEA from the endogenous locus (with retention of a single 34 bp loxP site at the corresponding location) results in reduced Tbx4 expression in the anterior portion of mutant hindlimbs. However, the loss of anterior expression is less severe than that observed with our HLEA-deleted BAC transgene (compare Fig. 5D with Fig. 7J), perhaps because additional cis-regulatory information is present in the endogenous locus that is absent from the BAC clone.

The factors acting upstream of HLEA and HLEB are not yet known. Previous studies have shown that ectopic Pitx1 expression is sufficient to trigger Tbx4 expression (Logan and Tabin, 1999Go). In addition, Pitx1-knockout mice show reduced, but not absent, expression of Tbx4 in the hindlimb, suggesting that Pitx1 is one of the endogenous factors that controls Tbx4 expression (Lanctot et al., 1999Go; Szeto et al., 1999Go). We found that mutation of a highly conserved Pitx1 binding site in HLEA reduces but does not eliminate hindlimb Tbx4 expression. The remaining enhancer activity might be due to additional binding sites for Pitx1, or to sites for additional factors that bind in the several-hundred-bp conserved HLEA enhancer region.

Tbx4 and hindlimb morphology
Early studies of Tbx4 focused on the possibility that Tbx4 expression helped specify forelimb versus hindlimb identity (Rodriguez-Esteban et al., 1999Go; Takeuchi et al., 1999Go). More recent studies have shown that Tbx4 expression is capable of rescuing forelimb outgrowth defects in a Tbx5 mutant mouse (Minguillon et al., 2005Go). The resulting limbs have morphological features characteristic of typical forelimbs, not hindlimbs, suggesting that Tbx4 functions in limb initiation, rather than limb identity. Conditional inactivation studies also show that early loss of Tbx4 leads to severe truncation of hindlimb outgrowth, whereas Tbx4 inactivation at later stages produces milder limb defects (Naiche and Papaioannou, 2007Go).

Our genetic knockout of HLEA provides a new method to bypass viability problems associated with loss-of-function mutations in the Tbx4 gene. Animals missing HLEA survive, are fertile, and have normal forelimb structures. Loss of HLEA leads to quantitative changes in the size of several different hindlimb bones, however, confirming that HLEA functions as a hindlimb enhancer and that Tbx4 is required for regional growth of specific skeletal structures in the vertebrate skeleton. The phenotypic effects of deleting HLEA are milder than the effects seen when disrupting all Tbx4 function in limbs beginning at E9.5-10.5 (Naiche and Papaioannou, 2007Go). The milder phenotypes seen in {Delta}HLEA mice might be due to remaining Tbx4 expression and function in limbs driven by HLEB, a possibility that can be tested in the future by similar knockouts of HLEB, or of both HLEA and HLEB, from the mouse Tbx4 locus.

Comparative studies have long suggested that vertebrates can independently control the size and shape of specific bones in forelimbs and hindlimbs (Flower, 1876Go; Hinchliffe and Johnson, 1980Go). Bats show striking elongation of digits in forelimbs but not hindlimbs. By contrast, kangaroos show robust development of hindlimb structures, and proportionately much shorter bones in forelimbs. Striking species-specific variation is also seen for individual bones within the hindlimb. The first digit ray of the human foot is unusually robust, giving rise to a characteristic big toe that supports a bipedal gait. By contrast, the anterior-most and posterior-most digits are nearly lost in horses, with weight-bearing transferred almost entirely onto a central hoofed digit.

The multiple independent hindlimb enhancers in the Tbx4 gene might provide a flexible genomic mechanism for influencing the quantitative size of skeletal elements in vertebrate limbs. Our genetic studies show that even heterozygous loss of one copy of HLEA can produce significant changes in the size of hindlimb bones (see Table S2 in the supplementary material). Fine-scale changes in the size of hindlimb skeletal elements can thus arise from quantitative changes in Tbx4 expression, and these expression levels are themselves controlled by at least two different hindlimb enhancer regions. Recent studies suggest that regulatory alterations in key developmental control genes, such as Shh, Pitx1 and Prx1 (Prrx1), play an important role in limb modifications in naturally occurring species (Sagai et al., 2004Go; Shapiro et al., 2004Go; Cretekos et al., 2008Go). In each case, null mutations in the corresponding gene are lethal, whereas regulatory mutations are viable and fertile. Our results show that mutations in limb-specific enhancers of Tbx4 can also provide a way to bypass the lethal, pleiotropic effects seen when the coding regions of the Tbx4 gene are disrupted. Based on the localized skeletal phenotypes seen in the current studies, we hypothesize that structural changes in HLEA and HLEB enhancer regions might be one of several factors that contribute to quantitative modifications in the size of hindlimb bones during vertebrate evolution. This possibility can now be tested by comparative sequencing of HLEA and HLEB in a range of organisms with different hindlimb morphologies, followed by functional tests in transgenic and knock-in mice (Cretekos et al., 2008Go).

Supplementary material
Supplementary material for this article is available at http://dev.biologists.org/cgi/content/full/135/15/2543/DC1


   ACKNOWLEDGMENTS
 
We thank Frank Chan, Hao Chen, Christie Ham, Craig Miller, Alex Pollen, Phil Reno and Brian Summers for comments on the manuscript, Felicity Jones for help with statistical analyses, and Frank Chan for graphics support. This work was supported by a Ruth L. Kirschstein National Research Service Award (F32HD048006, D.B.M.) and by a Center of Excellence in Genomic Science Grant (2P50HG002568, D.M.K.). D.M.K. is an investigator of the Howard Hughes Medical Institute.


   REFERENCES
TOP
 SUMMARY
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 


Abramoff, M. D., Magelhaes, P. J. and Ram, S. J. (2004). Image Processing with ImageJ. Biophotonics International 11,36 -42.
Agarwal, P., Wylie, J. N., Galceran, J., Arkhitko, O., Li, C., Deng, C., Grosschedl, R. and Bruneau, B. G. (2003). Tbx5 is essential for forelimb bud initiation following patterning of the limb field in the mouse embryo. Development 130,623 -633.[Abstract/Free Full Text]
Ahn, D. G., Kourakis, M. J., Rohde, L. A., Silver, L. M. and Ho, R. K. (2002). T-box gene tbx5 is essential for formation of the pectoral limb bud. Nature 417,754 -758.[CrossRef][Medline]
Basson, C. T., Bachinsky, D. R., Lin, R. C., Levi, T., Elkins, J. A., Soults, J., Grayzel, D., Kroumpouzou, E., Traill, T. A., Leblanc-Straceski, J. et al. (1997). Mutations in human TBX5 cause limb and cardiac malformation in Holt-Oram syndrome. Nat. Genet. 15,30 -35.[CrossRef][Medline]
Bell, M. A. and Foster, S. A. (1994). The Evolutionary Biology of the Threespine Stickleback. Oxford: Oxford University Press.
Blair, J. E. and Hedges, S. B. (2005). Molecular phylogeny and divergence times of deuterostome animals. Mol. Biol. Evol. 22,2275 -2284.[Abstract/Free Full Text]
Bongers, E. M., Duijf, P. H., van Beersum, S. E., Schoots, J., Van Kampen, A., Burckhardt, A., Hamel, B. C., Losan, F., Hoefsloot, L. H., Yntema, H. G. et al. (2004). Mutations in the human TBX4 gene cause small patella syndrome. Am. J. Hum. Genet. 74,1239 -1248.[CrossRef][Medline]
Brudno, M., Do, C. B., Cooper, G. M., Kim, M. F., Davydov, E., Green, E. D., Sidow, A. and Batzoglou, S. (2003). LAGAN and Multi-LAGAN: efficient tools for large-scale multiple alignment of genomic DNA. Genome Res. 13,721 -731.[Abstract/Free Full Text]
Capdevila, J. and Izpisua Belmonte, J. C. (2001). Patterning mechanisms controlling vertebrate limb development. Annu. Rev. Cell Dev. Biol. 17, 87-132.[CrossRef][Medline]
Chandler, R. L., Chandler, K. J., McFarland, K. A. and Mortlock, D. P. (2007). Bmp2 transcription in osteoblast progenitors is regulated by a distant 3' enhancer located 156.3 kilobases from the promoter. Mol. Cell. Biol. 27,2934 -2951.[Abstract/Free Full Text]
Chapman, D. L., Garvey, N., Hancock, S., Alexiou, M., Agulnik, S. I., Gibson-Brown, J. J., Cebra-Thomas, J., Bollag, R. J., Silver, L. M. and Papaioannou, V. E. (1996). Expression of the T-box family genes, Tbx1-Tbx5, during early mouse development. Dev. Dyn. 206,379 -390.[CrossRef][Medline]
Cole, N. J., Tanaka, M., Prescott, A. and Tickle, C. (2003). Expression of limb initiation genes and clues to the morphological diversification of threespine stickleback. Curr. Biol. 13,R951 -R952.[CrossRef][Medline]
Cretekos, C. J., Wang, Y., Green, E. D., Martin, J. F., Rasweiler, J. J., 4th and Behringer, R. R. (2008). Regulatory divergence modifies limb length between mammals. Genes Dev. 22,141 -151.[Abstract/Free Full Text]
Davis, M. C., Dahn, R. D. and Shubin, N. H. (2007). An autopodial-like pattern of Hox expression in the fins of a basal actinopterygian fish. Nature 447,473 -476.[CrossRef][Medline]
DiLeone, R. J., Russell, L. B. and Kingsley, D. M. (1998). An extensive 3' regulatory region controls expression of Bmp5 in specific anatomical structures of the mouse embryo. Genetics 148,401 -408.[Abstract/Free Full Text]
DiLeone, R. J., Marcus, G. A., Johnson, M. D. and Kingsley, D. M. (2000). Efficient studies of long-distance Bmp5 gene regulation using bacterial artificial chromosomes. Proc. Natl. Acad. Sci. USA 97,1612 -1617.[Abstract/Free Full Text]
Fisher, S., Grice, E. A., Vinton, R. M., Bessling, S. L. and McCallion, A. S. (2006). Conservation of RET regulatory function from human to zebrafish without sequence similarity. Science 312,276 -279.[Abstract/Free Full Text]
Flower, W. H. (1876). Osteology of the Mammalia. London: MacMillan.
Fortini, M. E. and Rubin, G. M. (1990). Analysis of cis-acting requirements of the Rh3 and Rh4 genes reveals a bipartite organization to rhodopsin promoters in Drosophila melanogaster. Genes Dev. 4,444 -463.[Abstract/Free Full Text]
Frazer, K. A., Pachter, L., Poliakov, A., Rubin, E. M. and Dubchak, I. (2004). VISTA: computational tools for comparative genomics. Nucleic Acids Res. 32,W273 -W279.[Abstract/Free Full Text]
Freitas, R., Zhang, G. and Cohn, M. J. (2007). Biphasic Hoxd gene expression in shark paired fins reveals an ancient origin of the distal limb domain. PLoS ONE 2, e754.[CrossRef]
Garrity, D. M., Childs, S. and Fishman, M. C. (2002). The heartstrings mutation in zebrafish causes heart/fin Tbx5 deficiency syndrome. Development 129,4635 -4645.[Abstract/Free Full Text]
Gibson-Brown, J. J., Agulnik, S. I., Chapman, D. L., Alexiou, M., Garvey, N., Silver, L. M. and Papaioannou, V. E. (1996). Evidence of a role for T-box genes in the evolution of limb morphogenesis and the specification of forelimb/hindlimb identity. Mech. Dev. 56,93 -101.[CrossRef][Medline]
Gibson-Brown, J. J., Agulnik, S. I., Silver, L. M. and Papaioannou, V. E. (1998). Expression of T-box genes Tbx2-Tbx5 during chick organogenesis. Mech. Dev. 74,165 -169.[CrossRef][Medline]
Hinchliffe, J. R. and Johnson, D. R. (1980). The Development of the Vertebrate Limb. Oxford: Clarendon Press.
Joyner, A. L. (1993). Gene Targeting. A Practical Approach. Oxford, UK: Oxford University Press.
Lamolet, B., Pulichino, A. M., Lamonerie, T., Gauthier, Y., Brue, T., Enjalbert, A. and Drouin, J. (2001). A pituitary cell-restricted T box factor, Tpit, activates POMC transcription in cooperation with Pitx homeoproteins. Cell 104,849 -859.[CrossRef][Medline]
Lamonerie, T., Tremblay, J. J., Lanctot, C., Therrien, M., Gauthier, Y. and Drouin, J. (1996). Ptx1, a bicoid-related homeo box transcription factor involved in transcription of the pro-opiomelanocortin gene. Genes Dev. 10,1284 -1295.[Abstract/Free Full Text]
Lanctot, C., Moreau, A., Chamberland, M., Tremblay, M. L. and Drouin, J. (1999). Hindlimb patterning and mandible development require the Ptx1 gene. Development 126,1805 -1810.[Abstract]
Lee, E. C., Yu, D., Martinez de Velasco, J., Tessarollo, L., Swing, D. A., Court, D. L., Jenkins, N. A. and Copeland, N. G. (2001). A highly efficient Escherichia coli-based chromosome engineering system adapted for recombinogenic targeting and subcloning of BAC DNA. Genomics 73,56 -65.[CrossRef][Medline]
Li, Q. Y., Newbury-Ecob, R. A., Terrett, J. A., Wilson, D. I., Curtis, A. R., Yi, C. H., Gebuhr, T., Bullen, P. J., Robson, S. C., Strachan, T. et al. (1997). Holt-Oram syndrome is caused by mutations in TBX5, a member of the Brachyury (T) gene family. Nat. Genet. 15,21 -29.[CrossRef][Medline]
Logan, M. and Tabin, C. J. (1999). Role of Pitx1 upstream of Tbx4 in specification of hindlimb identity. Science 283,1736 -1739.[Abstract/Free Full Text]
Logan, M., Simon, H. G. and Tabin, C. (1998). Differential regulation of T-box and homeobox transcription factors suggests roles in controlling chick limb-type identity. Development 125,2825 -2835.[Abstract]
Ludwig, M. Z., Palsson, A., Alekseeva, E., Bergman, C. M., Nathan, J. and Kreitman, M. (2005). Functional evolution of a cis-regulatory module. PLoS Biol. 3, e93.[CrossRef][Medline]
Lufkin, T., Mark, M., Hart, C. P., Dolle, P., LeMeur, M. and Chambon, P. (1992). Homeotic transformation of the occipital bones of the skull by ectopic expression of a homeobox gene. Nature 359,835 -841.[CrossRef][Medline]
Marcil, A., Dumontier, E., Chamberland, M., Camper, S. A. and Drouin, J. (2003). Pitx1 and Pitx2 are required for development of hindlimb buds. Development 130, 45-55.[Abstract/Free Full Text]
Mariani, F. V. and Martin, G. R. (2003). Deciphering skeletal patterning: clues from the limb. Nature 423,319 -325.[CrossRef][Medline]
Minguillon, C., Del Buono, J. and Logan, M. P. (2005). Tbx5 and Tbx4 are not sufficient to determine limb-specific morphologies but have common roles in initiating limb outgrowth. Dev. Cell 8,75 -84.[CrossRef][Medline]
Mortlock, D. P., Guenther, C. and Kingsley, D. M. (2003). A general approach for identifying distant regulatory elements applied to the Gdf6 gene. Genome Res. 13,2069 -2081.[Abstract/Free Full Text]
Naiche, L. A. and Papaioannou, V. E. (2003). Loss of Tbx4 blocks hindlimb development and affects vascularization and fusion of the allantois. Development 130,2681 -2693.[Abstract/Free Full Text]
Naiche, L. A. and Papaioannou, V. E. (2007). Tbx4 is not required for hindlimb identity or post-bud hindlimb outgrowth. Development 134,93 -103.[Abstract/Free Full Text]
Nelson, J. S. (1994). Fishes of the World. New York: Wiley-Interscience.
Pennacchio, L. A., Ahituv, N., Moses, A. M., Prabhakar, S., Nobrega, M. A., Shoukry, M., Minovitsky, S., Dubchak, I., Holt, A., Lewis, K. D. et al. (2006). In vivo enhancer analysis of human conserved non-coding sequences. Nature 444,499 -502.[CrossRef][Medline]
Rallis, C., Bruneau, B. G., Del Buono, J., Seidman, C. E., Seidman, J. G., Nissim, S., Tabin, C. J. and Logan, M. P. (2003). Tbx5 is required for forelimb bud formation and continued outgrowth. Development 130,2741 -2751.[Abstract/Free Full Text]
Rodriguez-Esteban, C., Tsukui, T., Yonei, S., Magallon, J., Tamura, K. and Izpisua Belmonte, J. C. (1999). The T-box genes Tbx4 and Tbx5 regulate limb outgrowth and identity. Nature 398,814 -818.[CrossRef][Medline]
Ruvinsky, I., Oates, A. C., Silver, L. M. and Ho, R. K. (2000). The evolution of paired appendages in vertebrates: T-box genes in the zebrafish. Dev. Genes Evol. 210, 82-91.[CrossRef][Medline]
Sagai, T., Masuya, H., Tamura, M., Shimizu, K., Yada, Y., Wakana, S., Gondo, Y., Noda, T. and Shiroishi, T. (2004). Phylogenetic conservation of a limb-specific, cis-acting regulator of Sonic hedgehog. Mamm. Genome 15, 23-34.[CrossRef][Medline]
Shalaby, F., Rossant, J., Yamaguchi, T. P., Gertsenstein, M., Wu, X. F., Breitman, M. L. and Schuh, A. C. (1995). Failure of blood-island formation and vasculogenesis in Flk-1-deficient mice. Nature 376,62 -66.[CrossRef][Medline]
Shapiro, M. D., Marks, M. E., Peichel, C. L., Blackman, B. K., Nereng, K. S., Jonsson, B., Schluter, D. and Kingsley, D. M. (2004). Genetic and developmental basis of evolutionary pelvic reduction in threespine sticklebacks. Nature 428,717 -723.[CrossRef][Medline]
Szeto, D. P., Rodriguez-Esteban, C., Ryan, A. K., O'Connell, S. M., Liu, F., Kioussi, C., Gleiberman, A. S., Izpisua-Belmonte, J. C. and Rosenfeld, M. G. (1999). Role of the Bicoid-related homeodomain factor Pitx1 in specifying hindlimb morphogenesis and pituitary development. Genes Dev. 13,484 -494.[Abstract/Free Full Text]
Takabatake, Y., Takabatake, T. and Takeshima, K. (2000). Conserved and divergent expression of T-box genes Tbx2-Tbx5 in Xenopus. Mech. Dev. 91,433 -437.[CrossRef][Medline]
Takeuchi, J. K., Koshiba-Takeuchi, K., Matsumoto, K., Vogel-Hopker, A., Naitoh-Matsuo, M., Ogura, K., Takahashi, N., Yasuda, K. and Ogura, T. (1999). Tbx5 and Tbx4 genes determine the wing/leg identity of limb buds. Nature 398,810 -814.[CrossRef][Medline]
Tamura, K., Yonei-Tamura, S. and Belmonte, J. C. (1999). Differential expression of Tbx4 and Tbx5 in Zebrafish fin buds. Mech. Dev. 87,181 -184.[CrossRef][Medline]
Tanaka, M., Munsterberg, A., Anderson, W. G., Prescott, A. R., Hazon, N. and Tickle, C. (2002). Fin development in a cartilaginous fish and the origin of vertebrate limbs. Nature 416,527 -531.[CrossRef][Medline]
Tickle, C. (2003). Patterning systems-from one end of the limb to the other. Dev. Cell 4, 449-458.[CrossRef][Medline]
Townley, D. J., Avery, B. J., Rosen, B. and Skarnes, W. C. (1997). Rapid sequence analysis of gene trap integrations to generate a resource of insertional mutations in mice. Genome Res. 7,293 -298.[Abstract/Free Full Text]
Tremblay, J. J., Goodyer, C. G. and Drouin, J. (2000). Transcriptional properties of Ptx1 and Ptx2 isoforms. Neuroendocrinology 71,277 -286.[CrossRef][Medline]
Wilkinson, D. G. (1992). Whole mount in situ hybridization of vertebrate embryos. In In Situ Hybridization, A Practical Approach (ed. D. G. Wilkinson), pp.75 -83. Oxford: IRL Press.
Woolfe, A., Goodson, M., Goode, D. K., Snell, P., McEwen, G. K., Vavouri, T., Smith, S. F., North, P., Callaway, H., Kelly, K. et al. (2005). Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 3, e7.[CrossRef][Medline]
Yang, X. W., Model, P. and Heintz, N. (1997). Homologous recombination based modification in Escherichia coli and germline transmission in transgenic mice of a bacterial artificial chromosome. Nat. Biotechnol. 15,859 -865.[CrossRef][Medline]
Zhao, S., Shatsman, S., Ayodeji, B., Geer, K., Tsegaye, G., Krol, M., Gebregeorgis, E., Shvartsbeyn, A., Russell, D., Overton, L. et al. (2001). Mouse BAC ends quality assessment and sequence analyses. Genome Res. 11,1736 -1745.[Abstract/Free Full Text]


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?



This Article
Right arrow Summary Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Supplementary Material
Right arrow All Versions of this Article:
dev.017384v1
135/15/2543    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Menke, D. B.
Right arrow Articles by Kingsley, D. M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Menke, D. B.
Right arrow Articles by Kingsley, D. M.
Social Bookmarking
Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?