Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.     

qPCR array检测分析 C57BL/6小鼠胚胎后肢发育基因表达谱

目的:胚胎生育过程中因肢体发育异常造成的出生缺陷比率不低,其相关基因表达模式尚不明确。本实验通过建立实时定量PCR芯片(Real-time quantitative polymerasechain reaction array,qPCR array)检测方案,研究C57BL/6品系小鼠后肢发育相关基因的表达谱。方法:以同源异形盒基因家族(Hox)、Wnt5a、配对同源结构域基因(Pitx1)、成纤维生长因子(Fgf8)、音猬因子(Shh)等小鼠肢体发育相关的重要基因制作基因检测表达谱,以C57BL/6品系怀孕雌鼠为材料,取胚胎肢芽发育的四个关键时期(E10.5,E11.5,E12.5,E13.5)的胎鼠后肢,利用qPCR array方案检测表达谱中基因的相对表达水平差异。结果:通过已建立的qPCR array检测了C57BL/6品系小鼠胚胎后肢发育时期Hox家族、Wnt5a、Pitx1、Fgf8、Shh等基因的表达差异。以E10.5为对照,检测出在后肢发育时期基因呈三种表达模式,即Hoxb6、Hoxb8、Hoxc8、Hoxc9、Hoxc10、Hoxd9和Shh基因的表达水平呈上调;Hoxa11、Hoxa13、Hoxc12、Hoxc13、Hoxd13等基因表达出现下调;Hoxc9、Hoxc10、Hoxc11、Hoxd9、Hoxd12、Fgf8和Pitx1等基因的相对表达量呈先上调后下调的曲线表达模式,且有少部分基因在小鼠后肢发育时期表达水平无明显变化。结论:Hox家族、Wnt5a、Pitx1、Fgf8、Shh等基因在小鼠后肢发育时期表达,并且表达模式存在明显差异。     

Developmental genetic basis for the evolution of pelvic fin loss in the pufferfish Takifugu rubripes

Paired appendages were a key developmental innovation among vertebrates and they eventually evolved into limbs. Ancient developmental control systems for paired fins and limbs are broadly conserved among gnathostome vertebrates. Some lineages including whales, some salamanders, snakes, and many ray-fin fish, independently lost the pectoral, pelvic, or both appendages over evolutionary time. When different taxa independently evolve similar developmental morphologies, do they use the same molecular genetic mechanisms? To determine the developmental genetic basis for the evolution of pelvis loss in the pufferfish Takifugu rubripes (fugu), we isolated fugu orthologs of genes thought to be essential for limb development in tetrapods, including limb positioning (Hoxc6, Hoxd9), limb bud initiation (Pitx1, Tbx4, Tbx5), and limb bud outgrowth (Shh, Fgf10), and studied their expression patterns during fugu development. Results showed that bud outgrowth and initiation fail to occur in fugu, and that pelvis loss is associated with altered expression of Hoxd9a, which we show to be a marker for pelvic fin position in three-spine stickleback Gasterosteus aculeatus. These results rule out changes in appendage outgrowth and initiation genes as the earliest developmental defect in pufferfish pelvic fin loss and suggest that altered Hoxd9a expression in the lateral mesoderm may account for pelvis loss in fugu. This mechanism appears to be different from the mechanism for pelvic loss in stickleback, showing that different taxa can evolve similar phenotypes by different mechanisms.     

HOXD13 may play a role in idiopathic congenital clubfoot by regulating the expression of FHL1

Idiopathic congenital clubfoot (CCF) is a type of congenital foot malformation, and previous studies using the extended transmission disequilibrium testing (ETDT) analysis have confirmed the HOXD13 gene as the susceptible gene of CCF. This study aimed to verify whether Hoxd13 directly regulates skeletal muscle LIM protein 1 (Fhl1) expression during limb development in rat embryo, which will serve as a basis for comprehending etiological research of idiopathic CCF. Immunofluorescence staining was utilized to detect Hoxd13 and Fhl1 expression in 12.5d rat embryo, while luciferase assay, electrophoretic mobility shift assay (EMSA), and chromatin immunoprecipitation (ChIP) assay were used to confirm the interaction between the two genes. Both Hoxd13 and Fhl1 were expressed in the interdigital tissues of E12.5 rat embryo. Luciferase assay and EMSA identified a novel promoter region of Fhl1 that directly interacts with Hoxd13. ChIP of the Hoxd13-Fhl1 promoter complex from the developing limb confirmed that endogenous Hoxd13 interacts with this region. Thus, Hoxd13 directly regulates Fhl1 expression in rat embryo. These findings suggest that HOXD13 may regulate the expression of FHL1 in the development of idiopathic CCF.     

Comparative analysis of 3D expression patterns of transcription factor genes and digit fate maps in the developing chick wing

Hoxd13, Tbx2, Tbx3, Sall1 and Sall3 genes are candidates for encoding antero-posterior positional values in the developing chick wing and specifying digit identity. In order to build up a detailed profile of gene expression patterns in cell lineages that give rise to each of the digits over time, we compared 3 dimensional (3D) expression patterns of these genes during wing development and related them to digit fate maps. 3D gene expression data at stages 21, 24 and 27 spanning early bud to digital plate formation, captured from in situ hybridisation whole mounts using Optical Projection Tomography (OPT) were mapped to reference wing bud models. Grafts of wing bud tissue from GFP chicken embryos were used to fate map regions of the wing bud giving rise to each digit; 3D images of the grafts were captured using OPT and mapped on to the same models. Computational analysis of the combined computerised data revealed that Tbx2 and Tbx3 are expressed in digit 3 and 4 progenitors at all stages, consistent with encoding stable antero-posterior positional values established in the early bud; Hoxd13 and Sall1 expression is more dynamic, being associated with posterior digit 3 and 4 progenitors in the early bud but later becoming associated with anterior digit 2 progenitors in the digital plate. Sox9 expression in digit condensations lies within domains of digit progenitors defined by fate mapping; digit 3 condensations express Hoxd13 and Sall1, digit 4 condensations Hoxd13, Tbx3 and to a lesser extent Tbx2. Sall3 is only transiently expressed in digit 3 progenitors at stage 24 together with Sall1 and Hoxd13; then becomes excluded from the digital plate. These dynamic patterns of expression suggest that these genes may play different roles in digit identity either together or in combination at different stages including the digit condensation stage.     

The Relationship between Gene Network Structure and Expression Variation among Individuals and Species

Variation among individuals is a prerequisite of evolution by natural selection. As such, identifying the origins of variation is a fundamental goal of biology. We investigated the link between gene interactions and variation in gene expression among individuals and species using the mammalian limb as a model system. We first built interaction networks for key genes regulating early (outgrowth; E9.5–11) and late (expansion and elongation; E11-13) limb development in mouse. This resulted in an Early (ESN) and Late (LSN) Stage Network. Computational perturbations of these networks suggest that the ESN is more robust. We then quantified levels of the same key genes among mouse individuals and found that they vary less at earlier limb stages and that variation in gene expression is heritable. Finally, we quantified variation in gene expression levels among four mammals with divergent limbs (bat, opossum, mouse and pig) and found that levels vary less among species at earlier limb stages. We also found that variation in gene expression levels among individuals and species are correlated for earlier and later limb development. In conclusion, results are consistent with the robustness of the ESN buffering among-individual variation in gene expression levels early in mammalian limb development, and constraining the evolution of early limb development among mammalian species.     

The chick oligozeugodactyly (ozd) mutant lacks sonic hedgehog function in the limb

We have analyzed a new limb mutant in the chicken that we name oligozeugodactyly (ozd). The limbs of this mutant have a longitudinal postaxial defect, lacking the posterior element in the zeugopod (ulna/fibula) and all digits except digit 1 in the leg. Classical recombination experiments show that the limb mesoderm is the defective tissue layer in ozd limb buds. Molecular analysis revealed that the ozd limbs develop in the absence of Shh expression, while all other organs express Shh and develop normally. Neither Ptc1 nor Gli1 are detectable in mutant limb buds. However, Bmp2 and dHAND are expressed in the posterior wing and leg bud mesoderm, although at lower levels than in normal embryos. Activation of Hoxd11-13 occurs normally in ozd limbs but progressively declines with time. Phase III of expression is more affected than phase II, and expression is more severely affected in the more 5' genes. Interestingly, re-expression of Hoxd13 occurs at late stages in the distal mesoderm of ozd leg buds, correlating with formation of digit 1. Fgf8 and Fgf4 expression are initiated normally in the mutant AER but their expression is progressively downregulated in the anterior AER. Recombinant Shh protein or ZPA grafts restore normal pattern to ozd limbs; however, retinoic acid fails to induce Shh in ozd limb mesoderm. We conclude that Shh function is required for limb development distal to the elbow/knee joints, similar to the Shh(-/-) mouse. Accordingly we classify the limb skeletal elements as Shh dependent or independent, with the ulna/fibula and digits other than digit 1 in the leg being Shh dependent. Finally we propose that the ozd mutation is most likely a defect in a regulatory element that controls limb-specific expression of Shh.     

Restricted patterns of Hoxd10 and Hoxd11 set segmental differences in motoneuron subtype complement in the lumbosacral spinal cord

During normal vertebrate development, Hoxd10 and Hoxd11 are expressed by differentiating motoneurons in restricted patterns along the rostrocaudal axis of the lumbosacral (LS) spinal cord. To assess the roles of these genes in the attainment of motoneuron subtypes characteristic of LS subdomains, we examined subtype complement after overexpression of Hoxd10 or Hoxd11 in the embryonic chick LS cord and in a Hoxd10 loss-of-function mouse embryo. Data presented here provide evidence that Hoxd10 defines the position of the lateral motor column (LMC) as a whole and, in rostral LS segments, specifically promotes the development of motoneurons of the lateral subdivision of the lateral motor column (LMCl). In contrast, Hoxd11 appears to impart a caudal and medial LMC (LMCm) identity to some motoneurons and molecular profiles suggestive of a suppression of LMC development in others. We also provide evidence that Hoxd11 suppresses the expression of Hoxd10 and the retinoic acid synthetic enzyme, retinaldehyde dehydrogenase 2 (RALDH2). In a normal chick embryo, Hoxd10 and RALDH2 are expressed throughout the LS region at early stages of motoneuron differentiation but their levels decline in Hoxd11-expressing caudal LS segments that ultimately contain few LMCl motoneurons. We hypothesize that one of the roles played by Hoxd11 is to modulate Hoxd10 and local retinoic acid levels and thus, perhaps define the caudal boundaries of the LMC and its subtype complement.     

Isolation, genomic structure and developmental expression of Fgf8 in the short-tailed fruit bat, Carollia perspicillata

Fibroblast growth factor-8 (Fgf8) encodes a secreted protein which was initially identified as the factor responsible for androgen-dependant growth of mouse mammary carcinoma cells (Tanaka et al., 1992). Fgf8 has been subsequently implicated in the patterning and growth of the gastrulating embryo, paraxial mesoderm (somites), limbs, craniofacial tissues, central nervous system and other organ systems during the development of several vertebrate model animals. Consistent with these findings, Fgf8 is expressed in a complex and dynamic pattern during vertebrate embryogenesis. Here we report the isolation and characterization of a bat (Carollia perspicillata) Fgf8 orthologue. Compared with those of other model vertebrates, Carollia Fgf8 is conserved with respect to genomic structure, sequence and many domains of developmental expression pattern. Interestingly, the expression domain marking the apical ectodermal ridge of the developing limb shows a striking difference compared to that of mouse, consistent with evolutionary diversification of bat limb morphology.     

Suppression of polydactyly of the Gli3 mutant (extra toes) by deltaEF1 homozygous mutation

Digit patterning is established through multiple genetic interactions. Delta-crystallin enhancer/E2-box factor (deltaEF1) is a zinc finger and homeodomain containing repressor protein, and is expressed in the posterior half of the forelimb bud and in the entire hindlimb bud during the early stage of limb development. The 6EF1-deficient mutant mice display various skeletal abnormalities, among which inferior ossification and abnormal patterning of autopodial bones are similar to those observed in Hox and Gli gene mutants. Gli3 mutant mice, extra toes (Xt), exhibit pre-axial polydactyly losing the identity of digit I. It is demonstrated here that deltaEF1null(lacZ) homozygosity suppressed formation of the extra digit, uniquely of the hindlimb, in both Gli3XtJ heterozygous and homozygous mutants, but with no restoration of digit I identity. In Gli3XtJ mutants, the Hoxd13 expression domain was expanded more dramatically in homozygotes. In Gli3XtJ;deltaEF1null(lacZ) double homozygous mutants, Hoxd13 expression once expanded in Gli3XtJ homozygous mutant was reduced, more conspicuously in the hindlimbs, which may account for hindlimb-restricted suppression of formation of the extra digit. The data suggest the possibility that the extent of Hoxd13 expression along the distal margin of the limb bud is determinative in defining the digit number.     

Hox genes, digit identities and the theropod/bird transition

Vargas and Fallon (2005. J Exp Zool (Mol Dev Evol) 304B:86-90) propose that Hox gene expression patterns indicate that the most anterior digit in bird wings is homologous to digit 1 rather than to digit 2 in other amniotes. This interpretation is based on the presence of Hoxd13 expression in combination with the absence of Hoxd12 expression in the second digit condensation from which this digit develops (the first condensation is transiently present). This is a pattern that is similar to that in the developing digit 1 of the chicken foot and the mouse hand and foot. They have tested this new hypothesis by analysing Hoxd12 and Hoxd13 expression patterns in two polydactylous chicken mutants, Silkie and talpid2. They conclude that the data support the notion that the most anterior remaining digit of the bird wing is homologous to digit 1 in other amniotes either in a standard phylogenetic sense, or alternatively in a (limited) developmental sense in agreement with the Frameshift Hypothesis of Wagner and Gautier (1999, i.e., that the developmental pathway is homologous to the one that leads to a digit 1 identity in other amniotes, although it occurs in the second instead of the first digit condensation). We argue that the Hoxd12 and Hoxd13 expression patterns found for these and other limb mutants do not allow distinguishing between the hypothesis of Vargas and Fallon (2005. J Exp Zool (Mol Dev Evol) 304B:86-90) and the alternative one, i.e., the most anterior digit in bird wings is homologous to digit 2 in other amniotes, in a phylogenetic or developmental sense. Therefore, at the moment the data on limb mutants does not present a challenge to the hypothesis, based on other developmental data (Holmgren, 1955. Acta Zool 36:243-328; Hinchliffe, 1984. In: Hecht M, Ostrom JH, Viohl G, Wellnhofer P, editors. The beginnings of birds. Eichst?tt: Freunde des Jura-Museum. p 141-147; Burke and Feduccia, 1997. Science 278:666-668; Kundrát et al., 2002. J Exp Zool (Mol Dev Evol) 294B:151-159; Larsson and Wagner, 2002. J Exp Zool (Mol Dev Evol) 294B:146-151; Feduccia and Nowicki, 2002. Naturwissenschaften 89:391-393), that the digits of bird wings are homologous to digits 2,3,4 in amniotes. We recommend further testing of the hypothesis by comparing Hoxd expression patterns in different taxa.     

Genetic Interactions Between Shox2 and Hox Genes During the Regional Growth and Development of the Mouse Limb

The growth and development of the vertebrate limb relies on homeobox genes of the Hox and Shox families, with their independent mutation often giving dose-dependent effects. Here we investigate whether Shox2 and Hox genes function together during mouse limb development by modulating their relative dosage and examining the limb for nonadditive effects on growth. Using double mRNA fluorescence in situ hybridization (FISH) in single embryos, we first show that Shox2 and Hox genes have associated spatial expression dynamics, with Shox2 expression restricted to the proximal limb along with Hoxd9 and Hoxa11 expression, juxtaposing the distal expression of Hoxa13 and Hoxd13. By generating mice with all possible dosage combinations of mutant Shox2 alleles and HoxA/D cluster deletions, we then show that their coordinated proximal limb expression is critical to generate normally proportioned limb segments. These epistatic interactions tune limb length, where Shox2 underexpression enhances, and Shox2 overexpression suppresses, Hox-mutant phenotypes. Disruption of either Shox2 or Hox genes leads to a similar reduction in Runx2 expression in the developing humerus, suggesting their concerted action drives cartilage maturation during normal development. While we furthermore provide evidence that Hox gene function influences Shox2 expression, this regulation is limited in extent and is unlikely on its own to be a major explanation for their genetic interaction. Given the similar effect of human SHOX mutations on regional limb growth, Shox and Hox genes may generally function as genetic interaction partners during the growth and development of the proximal vertebrate limb.