Axial homeosis and appendicular skeleton defects in mice with a targeted disruption of hoxd-11

Evolution of vertebral numbers in primates, with a focus on hominoids and the last common ancestor of hominins and panins

The primate vertebral column has been extensively studied, with a particular focus on hominoid primates and the last common ancestor of humans and chimpanzees. The number of vertebrae in hominoids-up to and including the last common ancestor of humans and chimpanzees-is subject to considerable debate. However, few formal ancestral state reconstructions exist, and none include a broad sample of primates or account for the correlated evolution of the vertebral column. Here, we conduct an ancestral state reconstruction using a model of evolution that accounts for both homeotic (changes of one type of vertebra to another) and meristic (addition or loss of a vertebra) changes. Our results suggest that ancestral primates were characterized by 29 precaudal vertebrae, with the most common formula being seven cervical, 13 thoracic, six lumbar, and three sacral vertebrae. Extant hominoids evolved tail loss and a reduced lumbar column via sacralization (homeotic transition at the last lumbar vertebra). Our results also indicate that the ancestral hylobatid had seven cervical, 13 thoracic, five lumbar, and four sacral vertebrae, and the ancestral hominid had seven cervical, 13 thoracic, four lumbar, and five sacral vertebrae. The last common ancestor of humans and chimpanzees likely either retained this ancestral hominid formula or was characterized by an additional sacral vertebra, possibly acquired through a homeotic shift at the sacrococcygeal border. Our results support the 'short-back' model of hominin vertebral evolution, which postulates that hominins evolved from an ancestor with an African ape-like numerical composition of the vertebral column.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Whole-genome resequencing to investigate the determinants of the multi-lumbar vertebrae trait in sheep

Multi-lumbar vertebrae trait is a beneficial mutation that can significantly improve livestock meat production. However, the genetic basis of the multi-lumbar vertebrae in sheep is still unclear. Here, we analysed the number of lumbar vertebrae of Duolang sheep and found three different traits of lumbar vertebrae number. Compared with the normal sheep, the length and weight of animal carcass from the multi-lumbar vertebrae sheep increased by 2.21 cm and 0.78 kg, respectively. We performed high-throughput genome resequencing on multi-lumbar vertebrae (n = 18) and normal (n = 11) Duolang sheep and obtained a total of more than 528.87 GB data. We found that the most significantly selective region were located in the 49.68-49.74 MB of chromosome 4 by selective-sweep analysis. We annotated this region and found that it contains SFRP4 which is known to regulate bone development. We further used the PCR-SSCP technology to detect the single nucleotide polymorphism (SNP) of the putative candidate SFRP4 and found that the two SNPs (rs600370085:C > T and rs415133338: A > G) of this gene were significantly associated with the multi-lumbar vertebrae of Duolang sheep. Our study indicates that the SFRP4 may be a potential major gene that affects the number of lumbar vertebrae in Duolang sheep, and has the potential to be utilized for sheep breeding in the future.

Genetic basis for an evolutionary shift from ancestral preaxial to postaxial limb polarity in non-urodele vertebrates

In most tetrapod vertebrates, limb skeletal progenitors condense with postaxial dominance. Posterior elements (such as ulna and fibula) appear prior to their anterior counterparts (radius and tibia), followed by digit-appearance order with continuing postaxial polarity. The only exceptions are urodele amphibians (salamanders), whose limb elements develop with preaxial polarity and who are also notable for their unique ability to regenerate complete limbs as adults. The mechanistic basis for this preaxial dominance has remained an enigma and has even been proposed to relate to the acquisition of novel genes involved in regeneration. However, recent fossil evidence suggests that preaxial polarity represents an ancestral rather than derived state. Here, we report that 5'Hoxd (Hoxd11-d13) gene deletion in mouse is atavistic and uncovers an underlying preaxial polarity in mammalian limb formation. We demonstrate this shift from postaxial to preaxial dominance in mouse results from excess Gli3 repressor (Gli3R) activity due to the loss of 5'Hoxd-Gli3 antagonism and is associated with cell-cycle changes promoting precocious cell-cycle exit in the anterior limb bud. We further show that Gli3 knockdown in axolotl results in a shift to postaxial dominant limb skeleton formation, as well as expanded paddle-shaped limb-bud morphology and ensuing polydactyly. Evolutionary changes in Gli3R activity level, which also played a key role in the fin-to-limb transition, appear to be fundamental to the shift from preaxial to postaxial polarity in formation of the tetrapod limb skeleton.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Homeotic change in segment identity derives the human vertebral formula from a chimpanzee-like one

OBJECTIVES

One of the most contentious issues in paleoanthropology is the nature of the last common ancestor of humans and our closest living relatives, chimpanzees and bonobos (panins). The numerical composition of the vertebral column has featured prominently, with multiple models predicting distinct patterns of evolution and contexts from which bipedalism evolved. Here, we study total numbers of vertebrae from a large sample of hominoids to quantify variation in and patterns of regional and total numbers of vertebrae in hominoids.

MATERIALS AND METHODS

We compile and study a large sample (N = 893) of hominoid vertebral formulae (numbers of cervical, thoracic, lumbar, sacral, caudal segments in each specimen) and analyze full vertebral formulae, total numbers of vertebrae, and super-regional numbers of vertebrae: presacral (cervical, thoracic, lumbar) vertebrae and sacrococcygeal vertebrae. We quantify within- and between-taxon variation using heterogeneity and similarity measures derived from population genetics.

RESULTS

We find that humans are most similar to African apes in total and super-regional numbers of vertebrae. Additionally, our analyses demonstrate that selection for bipedalism reduced variation in numbers of vertebrae relative to other hominoids.

DISCUSSION

The only proposed ancestral vertebral configuration for the last common ancestor of hominins and panins that is consistent with our results is the modal formula demonstrated by chimpanzees and bonobos (7 cervical-13 thoracic-4 lumbar-6 sacral-3 coccygeal). Hox gene expression boundaries suggest that a rostral shift in Hox10/Hox11-mediated complexes could produce the human modal formula from the proposal ancestral and panin modal formula.

HoxD transcription factors define monosynaptic sensory-motor specificity in the developing spinal cord

The specificity of monosynaptic connections between proprioceptive sensory neurons and their recipient spinal motor neurons depends on multiple factors, including motor neuron positioning and dendrite morphology, axon projection patterns of proprioceptive sensory neurons in the spinal cord, and the ligand-receptor molecules involved in cell-to-cell recognition. However, with few exceptions, the transcription factors engaged in this process are poorly characterized. Here, we show that members of the HoxD family of transcription factors play a crucial role in the specificity of monosynaptic sensory-motor connections. Mice lacking Hoxd9, Hoxd10 and Hoxd11 exhibit defects in locomotion but have no obvious defects in motor neuron positioning or dendrite morphology through the medio-lateral and rostro-caudal axes. However, we found that quadriceps motor neurons in these mice show aberrant axon development and receive inappropriate inputs from proprioceptive sensory axons innervating the obturator muscle. These genetic studies demonstrate that the HoxD transcription factors play an integral role in the synaptic specificity of monosynaptic sensory-motor connections in the developing spinal cord.

Latent developmental potential to form limb-like skeletal structures in zebrafish

Changes in appendage structure underlie key transitions in vertebrate evolution. Addition of skeletal elements along the proximal-distal axis facilitated critical transformations, including the fin-to-limb transition that permitted generation of diverse modes of locomotion. Here, we identify zebrafish mutants that form supernumerary long bones in their pectoral fins. These new bones integrate into musculature, form joints, and articulate with neighboring elements. This phenotype is caused by activating mutations in previously unrecognized regulators of appendage patterning, vav2 and waslb, that function in a common pathway. This pathway is required for appendage development across vertebrates, and loss of Wasl in mice causes defects similar to those seen in murine Hox mutants. Concordantly, formation of supernumerary bones requires Hox11 function, and mutations in the vav2/wasl pathway drive enhanced expression of hoxa11b, indicating developmental homology with the forearm. Our findings reveal a latent, limb-like pattern ability in fins that is activated by simple genetic perturbation.

Accelerated Evolution of Limb-Related Gene Hoxd11 in the Common Ancestor of Cetaceans and Ruminants (Cetruminantia)

Reduced numbers of carpal and tarsal bones (wrist and ankle joints) are extensively observed in the clade of Cetacea and Ruminantia (Cetruminantia). Homebox D11 (Hoxd11) is one of the important genes required for limb development in mammals. Mutations in Hoxd11 can lead to defects in particular bones of limbs, including carpus and tarsus. To test whether evolutionary changes in Hoxd11 underlie the loss of these bones in Cetruminantia, we sequenced and analyzed Hoxd11 coding sequences and compared them with other 5′ HoxA and HoxD genes in a taxonomic coverage of Cetacea, Ruminantia and other mammalian relatives. Statistical tests on the Hoxd11 sequences found an accelerated evolution in the common ancestor of cetaceans and ruminants, which coincided with the reduction of carpal and tarsal bones in this clade. Five amino acid substitutions (G222S, G227A, G229S, A240T and G261V) and one amino acid deletion (G254Del) occurred in this lineage. In contrast, other 5′ HoxA and HoxD genes do not show this same evolutionary pattern, but instead display a highly conserved pattern of evolution in this lineage. Accelerated evolution of Hoxd11, but not other 5′ HoxA and HoxD genes, is probably related to the reduction of the carpal and tarsal bones in Cetruminantia. Moreover, we found two amino acid substitutions (G110S and D223N) in Hoxd11 that are unique to the lineage of Cetacea, which coincided with hindlimb loss in the common ancestor of cetaceans. Our results give molecular evidence of Hoxd11 adaptive evolution in cetaceans and ruminants, which could be correlated with limb morphological adaptation.

Bone morphology is regulated modularly by global and regional genetic programs

Bone protrusions provide stable anchoring sites for ligaments and tendons and define the unique morphology of each long bone. Despite their importance, the mechanism by which superstructures are patterned is unknown. Here, we identify components of the genetic program that control the patterning of Sox9 ⁺/Scx ⁺ superstructure progenitors in mouse and show that this program includes both global and regional regulatory modules. Using light-sheet fluorescence microscopy combined with genetic lineage labeling, we mapped the broad contribution of the Sox9 ⁺/Scx ⁺ progenitors to the formation of bone superstructures. Then, by combining literature-based evidence, comparative transcriptomic analysis and genetic mouse models, we identified Gli3 as a global regulator of superstructure patterning, whereas Pbx1, Pbx2, Hoxa11 and Hoxd11 act as proximal and distal regulators, respectively. Moreover, by demonstrating a dose-dependent pattern regulation in Gli3 and Pbx1 compound mutations, we show that the global and regional regulatory modules work in a coordinated manner. Collectively, our results provide strong evidence for genetic regulation of superstructure patterning, which further supports the notion that long bone development is a modular process.This article has an associated 'The people behind the papers' interview.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Disruption of Hox9,10,11 function results in cellular level lineage infidelity in the kidney

Hox genes are important regulators of development. The 39 mammalian Hox genes have considerable functional overlap, greatly confounding their study. In this report, we generated mice with multiple combinations of paralogous and flanking Abd-B Hox gene mutations to investigate functional redundancies in kidney development. The resulting mice developed a number of kidney abnormalities, including hypoplasia, agenesis, and severe cysts, with distinct Hox functions observed in early metanephric kidney formation and nephron progenitor maintenance. Most surprising, however, was that extensive removal of Hox shared function in these kidneys resulted in cellular level lineage infidelity. Strikingly, mutant nephron tubules consisted of intermixed cells with proximal tubule, loop of Henle, and collecting duct identities, with some single cells expressing markers associated with more than one nephron segment. These results indicate that Hox genes are required for proper lineage selection/maintenance and full repression of genes involved in cell fate restriction in the developing kidney.

Reduced Abd-B Hox function during kidney development results in lineage infidelity

Hox genes can function as key drivers of segment identity, with Hox mutations in Drosophila often resulting in dramatic homeotic transformations. In addition, however, they can serve other essential functions. In mammals, the study of Hox gene roles in development is complicated by the presence of four Hox clusters with a total of 39 genes showing extensive functional overlap. In this study, in order to better understand shared core Hox functions, we examined kidney development in mice with frameshift mutations of multiple Abd-B type Hox genes. The resulting phenotypes included dramatically reduced branching morphogenesis of the ureteric bud, premature depletion of nephron progenitors and abnormal development of the stromal compartment. Most unexpected, however, we also observed a cellular level lineage infidelity in nephron segments. Scattered cells within the proximal tubules, for example, expressed genes normally expressed only in collecting ducts. Multiple combinations of inappropriate nephron segment specific marker expression were found. In some cases, cells within a tubule showed incorrect identity, while in other cases cells showed ambiguous character, with simultaneous expression of genes associated with more than one nephron segment. These results give evidence that Hox genes have an overlapping core function at the cellular level in driving and/or maintaining correct differentiation decisions.

Proximate cause, anatomical correlates, and obstetrical implication of a supernumerary lumbar vertebra in humans

OBJECTIVES

Three issues are considered on variation in number of presacral vertebrae (PSV) in humans: (1) sexual difference in number of PSV, (2) inactivation of Hoxd-11 gene as etiology for a supernumerary lumbar vertebra, and (3) anatomical correlates of a supernumerary lumbar vertebra, including lumbar-sacral nearthrosis, and pelvic size.

MATERIALS AND METHODS

Sample was 407 skeletonized females and 1,318 males from United States; ages at death were 20 to 49 years. Two subsamples of males were used: (1) 98 with modal numbers of cervical, thoracic, lumbar, and sacral vertebrae (PSV = 24) and (2) 45 with a supernumerary lumbar vertebra but modal numbers for other vertebral segments (PSV = 25). Measurements were taken of ulna, second metacarpal, vertebrae, femur, and pelvis; presence of lumbar-sacral nearthrosis was observed.

RESULTS

Although 90% of females and males have 24 PSV, females have higher frequency of 23 PSV and males have higher frequency of 25 PSV. Compared to males with 24 PSV, males with 25 PSV and supernumerary lumbar vertebra show (1) no difference in anatomies associated with inactivation of Hoxd-11, and (2) higher frequency of lumbar-sacral nearthrosis and smaller pelvic inlet circumference.

DISCUSSION

Sexual difference in number of PSV may be due to tempo of somite formation and Hox gene activation. Hypothesis is not supported that a supernumerary lumbar vertebra is due to inactivation of Hoxd-11. The presence of a supernumerary lumbar vertebra is associated with small pelvic inlet circumference, which can be obstetrically disadvantageous.

Sacral anatomy of the phytosaur Smilosuchus adamanensis, with implications for pelvic girdle evolution among Archosauriformes

The sacrum - consisting of those vertebrae that articulate with the ilia - is the exclusive skeletal connection between the hindlimbs and axial skeleton in tetrapods. Therefore, the morphology of this portion of the vertebral column plays a major role in the evolution of terrestrial locomotion. Whereas most extant reptiles only possess the two plesiomorphic sacral vertebrae, additional vertebrae have been incorporated into the sacrum multiple times independently among early-diverging archosaurian (crocodylians + birds) clades. Phytosauria was a diverse, abundant, and cosmopolitan clade of archosauriforms throughout the Late Triassic, but postcrania of this clade are rarely described and few species-level taxonomic placements of phytosaurian postcranial material are available, potentially hampering knowledge of morphological disparity in the postcranial skeleton among phytosaurs. Here, we describe the sacrum of Smilosuchus adamanensis, a phytosaur recovered from the Upper Triassic Chinle Formation of Arizona. This sacrum consists of the two primordial sacral vertebrae, but has a vertebra incorporated from the trunk into the sacrum (= a dorsosacral) and is therefore the first Late Triassic phytosaur and one of the first non-archosaurian archosauromorphs to be described with more than two sacral vertebrae. Our interpretation of this element as a dorsosacral is justified by the lateral extent of the dorsosacral ribs, clear surfaces of articulation between the distal ends of the dorsosacral ribs and the first primordial sacral ribs, and the scar on the medial surface of each ilium for articulation with each dorsosacral rib. Additionally, we provide the first detailed description of the vertebral junction formed by two anteriorly projecting flanges on the first primordial sacral ribs and their corresponding facets on the centrum of the dorsosacral. Computed tomographic (CT) imaging reveals that the two primordial sacrals are not co-ossified and that the dorsosacral morphology of this specimen is not the result of obvious pathology. We place this incorporation of a trunk vertebra into the phytosaurian sacrum in a broader evolutionary context, with this shift in vertebral identity occurring at least seven times independently among Triassic archosauriforms, including at least three times in early crocodylian-line archosaurs and at least four times among bird-line archosaurs. Additionally, anteriorly projecting flanges of sacral ribs which articulate with the anterior-adjacent centrum have evolved several times in archosauriforms, and we interpret 'shared' sacral ribs (= a sacral rib that articulates with two adjacent sacral centra more or less equally) present in some archosaurian clades as a more extreme example of this morphology. In extant taxa the highly conserved Hox gene family plays a central role in the patterning of the axial skeleton, especially vertebral identity; therefore, the independent incorporation of a trunk vertebra into the sacrum across multiple archosauriform lineages may suggest a homologous underlying developmental mechanism for this evolutionary trend.

A genome-wide association study identifies a genomic region for the polycerate phenotype in sheep (Ovis aries)

Horns are a cranial appendage found exclusively in Bovidae, and play important roles in accessing resources and mates. In sheep (Ovies aries), horns vary from polled to six-horned, and human have been selecting polled animals in farming and breeding. Here, we conducted a genome-wide association study on 24 two-horned versus 22 four-horned phenotypes in a native Chinese breed of Sishui Fur sheep. Together with linkage disequilibrium (LD) analyses and haplotype-based association tests, we identified a genomic region comprising 132.0–133.1 Mb on chromosome 2 that contained the top 10 SNPs (including 4 significant SNPs) and 5 most significant haplotypes associated with the polycerate phenotype. In humans and mice, this genomic region contains the HOXD gene cluster and adjacent functional genes EVX2 and KIAA1715, which have a close association with the formation of limbs and genital buds. Our results provide new insights into the genetic basis underlying variable numbers of horns and represent a new resource for use in sheep genetics and breeding.

Developmental Patterning as a Quantitative Trait: Genetic Modulation of the Hoxb6 Mutant Skeletal Phenotype

The process of patterning along the anterior-posterior axis in vertebrates is highly conserved. The function of Hox genes in the axis patterning process is particularly well documented for bone development in the vertebral column and the limbs. We here show that Hoxb6, in skeletal elements at the cervico-thoracic junction, controls multiple independent aspects of skeletal pattern, implicating discrete developmental pathways as substrates for this transcription factor. In addition, we demonstrate that Hoxb6 function is subject to modulation by genetic factors. These results establish Hox-controlled skeletal pattern as a quantitative trait modulated by gene-gene interactions, and provide evidence that distinct modifiers influence the function of conserved developmental genes in fundamental patterning processes.

Hox11 genes regulate postnatal longitudinal bone growth and growth plate proliferation

Hox genes are critical regulators of skeletal development and Hox9-13 paralogs, specifically, are necessary for appendicular development along the proximal to distal axis. Loss of function of both Hoxa11 and Hoxd11 results in severe malformation of the forelimb zeugopod. In the radius and ulna of these mutants, chondrocyte development is perturbed, growth plates are not established, and skeletal growth and maturation fails. In compound mutants in which one of the four Hox11 alleles remains wild-type, establishment of a growth plate is preserved and embryos develop normally through newborn stages, however, skeletal phenotypes become evident postnatally. During postnatal development, the radial and ulnar growth rate slows compared to wild-type controls and terminal bone length is reduced. Growth plate height is decreased in mutants and premature growth plate senescence occurs along with abnormally high levels of chondrocyte proliferation in the reserve and proliferative zones. Compound mutants additionally develop an abnormal curvature of the radius, which causes significant distortion of the carpal elements. The progressive bowing of the radius appears to result from physical constraint caused by the disproportionately slower growth of the ulna than the radius. Collectively, these data are consistent with premature depletion of forelimb zeugopod progenitor cells in the growth plate of Hox11 compound mutants, and demonstrate a continued function for Hox genes in postnatal bone growth and patterning.

Mesenchymal Hox6 function is required for mouse pancreatic endocrine cell differentiation

Despite significant advances in our understanding of pancreatic endocrine cell development, the function of the pancreatic mesodermal niche in this process is poorly understood. Here we report a novel role for mouse Hox6 genes in pancreatic organogenesis. Hox6 genes are expressed exclusively in the mesoderm of the developing pancreas. Genetic loss of all three Hox6 paralogs (Hoxa6, Hoxb6 and Hoxc6) leads to a dramatic loss of endoderm-derived endocrine cells, including insulin-secreting β-cells, and to mild delays and disruptions in pancreatic branching and exocrine differentiation. Ngn3-expressing pan-endocrine progenitor cells are specified normally in Hox6 mutant pancreata, but fail to mature into hormone-producing cells. Reduced expression of Wnt5a is observed in mutant pancreatic mesenchyme, leading to subsequent loss of expression of the crucial Wnt inhibitors Sfrp3 and Dkk1 in endocrine progenitor cells. These results reveal a key role for Hox6 genes in establishing Wnt mesenchymal-epithelial crosstalk in pancreatic development.

DNA-Methylation Patterns in Trisomy 21 Using Cells from Monozygotic Twins

DNA methylation is essential in mammalian development. We have hypothesized that methylation differences induced by trisomy 21 (T21) contribute to the phenotypic characteristics and heterogeneity in Down syndrome (DS). In order to determine the methylation differences in T21 without interference of the interindividual genomic variation, we have used fetal skin fibroblasts from monozygotic (MZ) twins discordant for T21. We also used skin fibroblasts from MZ twins concordant for T21, normal MZ twins without T21, and unrelated normal and T21 individuals. Reduced Representation Bisulfite Sequencing (RRBS) revealed 35 differentially methylated promoter regions (DMRs) (Absolute methylation differences = 25%, FDR < 0.001) in MZ twins discordant for T21 that have also been observed in comparison between unrelated normal and T21 individuals. The identified DMRs are enriched for genes involved in embryonic organ morphogenesis (FDR = 1.60 e -03) and include genes of the HOXB and HOXD clusters. These DMRs are maintained in iPS cells generated from this twin pair and are correlated with the gene expression changes. We have also observed an increase in DNA methylation level in the T21 methylome compared to the normal euploid methylome. This observation is concordant with the up regulation of DNA methyltransferase enzymes (DNMT3B and DNMT3L) and down regulation of DNA demethylation enzymes (TET2 and TET3) observed in the iPSC of the T21 versus normal twin. Altogether, the results of this study highlight the epigenetic effects of the extra chromosome 21 in T21 on loci outside of this chromosome that are relevant to DS associated phenotypes.

Independent regulation of vertebral number and vertebral identity by microRNA-196 paralogs

The Hox genes play a central role in patterning the embryonic anterior-to-posterior axis. An important function of Hox activity in vertebrates is the specification of different vertebral morphologies, with an additional role in axis elongation emerging. The miR-196 family of microRNAs (miRNAs) are predicted to extensively target Hox 3' UTRs, although the full extent to which miR-196 regulates Hox expression dynamics and influences mammalian development remains to be elucidated. Here we used an extensive allelic series of mouse knockouts to show that the miR-196 family of miRNAs is essential both for properly patterning vertebral identity at different axial levels and for modulating the total number of vertebrae. All three miR-196 paralogs, 196a1, 196a2, and 196b, act redundantly to pattern the midthoracic region, whereas 196a2 and 196b have an additive role in controlling the number of rib-bearing vertebra and positioning of the sacrum. Independent of this, 196a1, 196a2, and 196b act redundantly to constrain total vertebral number. Loss of miR-196 leads to a collective up-regulation of numerous trunk Hox target genes with a concomitant delay in activation of caudal Hox genes, which are proposed to signal the end of axis extension. Additionally, we identified altered molecular signatures associated with the Wnt, Fgf, and Notch/segmentation pathways and demonstrate that miR-196 has the potential to regulate Wnt activity by multiple mechanisms. By feeding into, and thereby integrating, multiple genetic networks controlling vertebral number and identity, miR-196 is a critical player defining axial formulae.

Key pathways regulated by HoxA9,10,11/HoxD9,10,11 during limb development

Background

The 39 mammalian Hox genes show problematic patterns of functional overlap. In order to more fully define the developmental roles of Hox genes it is necessary to remove multiple combinations of paralogous and flanking genes. In addition, the downstream molecular pathways regulated by Hox genes during limb development remain incompletely delineated.

Results

In this report we examine limb development in mice with frameshift mutations in six Hox genes, Hoxa9,10,11 and Hoxd9,10,11. The mice were made with a novel recombineering method that allows the simultaneous targeting of frameshift mutations into multiple flanking genes. The Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant mice show a reduced ulna and radius that is more severe than seen in Hoxa11^−/−/Hoxd11^−/− mice, indicating a minor role for the flanking Hox9,10 genes in zeugopod development, as well as their primary function in stylopod development. The mutant mice also show severe reduction of Shh expression in the zone of polarizing activity, and decreased Fgf8 expression in the apical ectodermal ridge, thereby better defining the roles of these specific Hox genes in the regulation of critical signaling centers during limb development. Importantly, we also used laser capture microdissection coupled with RNA-Seq to characterize the gene expression programs in wild type and mutant limbs. Resting, proliferative and hypertrophic compartments of E15.5 forelimb zeugopods were examined. The results provide an RNA-Seq characterization of the progression of gene expression patterns during normal endochondral bone formation. In addition of the Hox mutants showed strongly altered expression of Pknox2, Zfp467, Gdf5, Bmpr1b, Dkk3, Igf1, Hand2, Shox2, Runx3, Bmp7 and Lef1, all of which have been previously shown to play important roles in bone formation.

Conclusions

The recombineering based frameshift mutation of the six flanking and paralogous Hoxa9,10,11 and Hoxd9,10,11 genes provides a resource for the analysis of their overlapping functions. Analysis of the Hoxa9,10,11^−/−/Hoxd9,10,11^−/− mutant limbs confirms and extends the results of previous studies using mice with Hox mutations in single paralogous groups or with entire Hox cluster deletions. The RNA-Seq analysis of specific compartments of the normal and mutant limbs defines the multiple key perturbed pathways downstream of these Hox genes.

Electronic supplementary material

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

Suppressor of Fused Is Required for Determining Digit Number and Identity via Gli3/Fgfs/Gremlin

The anterior-posterior patterning of the vertebrate limb bud requires closely coordinated signaling interactions, including Sonic Hedgehog (Shh)-mediated counteraction of the Gli3 transcription factor in the distal and posterior mesenchyme of the limb bud. Suppressor of Fused (Sufu), an intracellular negative regulator of Shh signaling via Gli2 and Gli3, is implicated in early development of the mouse limb bud. However, how Sufu is involved in the genetic regulation of limb bud patterning still remains elusive. In this study, we show that the conditional deletion of Sufu in the mesenchyme of the early limb bud results in polydactyly with loss of digit identity and supernumerary bones in the wrist and the ankle. These pattern alterations are associated with anterior expansion of HoxD genes located at the 5’ end of the cluster. By focusing on gene expression analysis of Shh/Gremlin1/Fgf signaling critical for the establishment and maintenance of anterior-posterior patterning, we show that early response to loss of Sufu involves anterior prolongation of Fgf4 and Fgf8 expression in the apical ectodermal ridge at E10.5. We also reveal the anterior activation of Shh-dependent posterior markers Ptc1, Gli1 and Gremlin in limb buds lacking Sufu. Furthermore, we find that loss of Sufu leads to attenuated levels of repressor Gli2 and repressor Gli3 in the early limb bud. Moreover, expression of Hand2 is activated in the entire limb bud at the early outgrowth stage in the mutant lacking Sufu. Thus, we provide evidence that Sufu is involved in the genetic network that restricts the posterior expression of Gli2/3/Hand2 and Gremlin/Fgf in limb bud patterning.

Insights in the etiopathology of galactosyltransferase II (GalT-II) deficiency from transcriptome-wide expression profiling of skin fibroblasts of two sisters with compound heterozygosity for two novel B3GALT6 mutations

Mutations in B3GALT6, encoding the galactosyltransferase II (GalT-II) involved in the synthesis of the glycosaminoglycan (GAG) linkage region of proteoglycans (PGs), have recently been associated with a spectrum of connective tissue disorders, including spondyloepimetaphyseal dysplasia with joint laxity type 1 (SEMDJL1) and Ehlers–Danlos-like syndrome. Here, we report on two sisters compound heterozygous for two novel B3GALT6 mutations that presented with severe short stature and progressive kyphoscoliosis, joint hypermobility and laxity, hyperextensible skin, platyspondyly, short ilia, and elbow malalignment. Microarray-based transcriptome analysis revealed the differential expression of several genes encoding extracellular matrix (ECM) structural components, including COMP, SPP1, COL5A1, and COL15A1, enzymes involved in GAG synthesis and in ECM remodeling, such as CSGALNACT1, CHPF, LOXL3, and STEAP4, signaling transduction molecules of the TGFβ/BMP pathway, i.e., GDF6, GDF15, and BMPER, and transcription factors of the HOX and LIM families implicated in skeletal and limb development. Immunofluorescence analyses confirmed the down-regulated expression of some of these genes, in particular of the cartilage oligomeric matrix protein and osteopontin, encoded by COMP and SPP1, respectively, and showed the predominant reduction and disassembly of the heparan sulfate specific GAGs, as well as of the PG perlecan and type III and V collagens. The key role of GalT-II in GAG synthesis and the crucial biological functions of PGs are consistent with the perturbation of many physiological functions that are critical for the correct architecture and homeostasis of various connective tissues, including skin, bone, cartilage, tendons, and ligaments, and generates the wide phenotypic spectrum of GalT-II-deficient patients.

•

Clinical features/molecular characterization of two patients with spondyloepimetaphyseal dysplasia with joint laxity type 1

•

Identification of two novel B3GALT6 mutations

•

First report of transcriptome-wide gene expression profiling on GalT-II-deficient fibroblasts

•

Immunofluorescence studies of several ECM structural components in GalT-II-deficient cells

•

Enlargement of the knowledge on the GalT-II deficiency’s molecular pathogenesis

The pisiform growth plate is lost in humans and supports a role for Hox in growth plate formation

The human pisiform is a small, nodular, although functionally significant, bone of the wrist. In most other mammals, including apes and Australopithecus afarensis, pisiforms are elongate. An underappreciated fact is that the typical mammalian pisiform forms from two ossification centers. We hypothesize that: (i) the presence of a secondary ossification center in mammalian pisiforms indicates the existence of a growth plate; and (ii) human pisiform reduction results from growth plate loss. To address these hypotheses, we surveyed African ape pisiform ossification and confirmed the presence of a late-forming secondary ossification center in chimpanzees and gorillas. Identification of the initial ossification center occurs substantially earlier in apes relative to humans, raising questions concerning the homology of the human pisiform and the two mammalian ossification centers. Second, we conducted histological and immunohistochemical analyses of pisiform ossification in mice. We confirm the presence of two ossification centers separated by organized columnar and hypertrophic chondrocyte zones. Flattened chondrocytes were highly mitotic, indicating the presence of a growth plate. Hox genes have been proposed to play a fundamental role in growth plate patterning. The existence of a pisiform growth plate presents an interesting test case for the association between Hox expression and growth plate formation, and could explain the severe effects on the pisiform observed in Hoxa11 and Hoxd11 knockout mice. Consistent with this hypothesis, we show that Hoxd11 is expressed adjacent to the pisiform in late-stage embryonic mouse limbs supporting a role for Hox genes in growth plate specification. This raises questions concerning the mechanisms regulating Hox expression in the developing carpus.

Long bone development requires a threshold of Hox function

The Hoxd(Del(11-13)) mutant is one of the animal models for human synpolydactyly, characterized by short and syndactylous digits. Here we have characterized in detail the cartilage and bone defects in these mutants. We report two distinct phenotypes: (i) a delay and change in pattern of chondrocyte maturation of metacarpals/metatarsals and (ii) formation of a poor and not centrally positioned primary ossification center in the proximal-intermediate phalanx. In the metacarpals of Hoxd(Del(11-13)) mutants, ossification occurs postnataly, in the absence of significant Ihh expression and without the establishment of growth plates, following patterns similar to those of short bones. The strong downregulation in Ihh expression is associated with a corresponding increase of the repressor form of Gli3. To evaluate the contribution of this alteration to the phenotype, we generated double Hoxd(Del(11-13));Gli3 homozygous mutants. Intriguingly, these double mutants showed a complete rescue of the phenotype in metatarsals but only partial phenotypic rescue in metacarpals. Our results support Hox genes being required in a dose-dependent manner for long bone cartilage maturation and suggest that and excess of Gli3R mediates a significant part of the Hoxd(Del(11-13)) chondrogenic phenotype.

Unique expression patterns of multiple key genes associated with the evolution of mammalian flight

Bats are the only mammals capable of true flight. Critical adaptations for flight include a pair of dramatically elongated hands with broad wing membranes. To study the molecular mechanisms of bat wing evolution, we perform genomewide mRNA sequencing and in situ hybridization for embryonic bat limbs. We identify seven key genes that display unique expression patterns in embryonic bat wings and feet, compared with mouse fore- and hindlimbs. The expression of all 5'HoxD genes (Hoxd9-13) and Tbx3, six known crucial transcription factors for limb and digit development, is extremely high and prolonged in the elongating wing area. The expression of Fam5c, a tumour suppressor, in bat limbs is bat-specific and significantly high in all short digit regions (the thumb and foot digits). These results suggest multiple genetic changes occurred independently during the evolution of bat wings to elongate the hand digits, promote membrane growth and keep other digits short. Our findings also indicate that the evolution of limb morphology depends on the complex integration of multiple gene regulatory networks and biological processes that control digit formation and identity, chondrogenesis, and interdigital regression or retention.

Human HOX gene disorders

The Hox genes are an evolutionarily conserved family of genes, which encode a class of important transcription factors that function in numerous developmental processes. Following their initial discovery, a substantial amount of information has been gained regarding the roles Hox genes play in various physiologic and pathologic processes. These processes range from a central role in anterior-posterior patterning of the developing embryo to roles in oncogenesis that are yet to be fully elucidated. In vertebrates there are a total of 39 Hox genes divided into 4 separate clusters. Of these, mutations in 10 Hox genes have been found to cause human disorders with significant variation in their inheritance patterns, penetrance, expressivity and mechanism of pathogenesis. This review aims to describe the various phenotypes caused by germline mutation in these 10 Hox genes that cause a human phenotype, with specific emphasis paid to the genotypic and phenotypic differences between allelic disorders. As clinical whole exome and genome sequencing is increasingly utilized in the future, we predict that additional Hox gene mutations will likely be identified to cause distinct human phenotypes. As the known human phenotypes closely resemble gene-specific murine models, we also review the homozygous loss-of-function mouse phenotypes for the 29 Hox genes without a known human disease. This review will aid clinicians in identifying and caring for patients affected with a known Hox gene disorder and help recognize the potential for novel mutations in patients with phenotypes informed by mouse knockout studies.

Nonsyndromic brachydactyly type D and type E mapped to 7p15 in healthy children and adults from the Jirel ethnic group in eastern Nepal

OBJECTIVES

There is phenotypic overlap between Brachydactyly Type D (BDD) and Brachydactyly Type E (BDE) that suggests a possible common underlying etiology. We seek to understand the genetic underpinnings of, and relationship between, these skeletal anomalies.

METHODS

The Jirel ethnic group of eastern Nepal participates in various genetic epidemiologic studies, including those in which hand-wrist radiographs have been taken to examine skeletal development. Nearly 2,130 individuals (969 males; 1,161 females) were phenotyped for BDD/BDE. Of these, 1,722 individuals (773 males; 949 females) were genotyped for 371 STR markers spanning the autosomal genome. Variance components-based linkage analysis was used to conduct a genome-wide linkage scan for QTL influencing the BDD/BDE phenotype.

RESULTS

BDD was present in 3.55%, and BDE was present in 0.39%, of the study sample. Because of the phenotypic overlap between two traits, affecteds of either type were considered as affected by a single combined phenotype (BDD/BDE) having a prevalence of 3.94%. The additive genetic heritability of BDD/BDE was highly significant (h(2) ± SE = 0.89 ± 0.13; P = 1.7 × 10(-11) ). Significant linkage of BDD/BDE was found to markers on chromosome 7p21-7p14 (peak LOD score = 3.74 at 7p15 between markers D7S493 and D7S516).

CONCLUSIONS

Possible positional candidate genes in the one-lod support interval of this QTL include TWIST and the HOXA1-A13 cluster. This is the first study to report significant linkage results for BDD/BDE using a large extended pedigree, and the first to suggest that mutations in TWIST and/or the HOXA1-A13 cluster may contribute to these specific skeletal anomalies.

Direct activation of a mouse Hoxd11 axial expression enhancer by Gdf11/Smad signalling

A Hoxd11/lacZ reporter, expressed with a Hoxd11-like axial expression pattern in transgenic mouse embryos, is stimulated in tailbud fragments when cultured in presence of Gdf11, a TGF-β growth/differentiation factor. The same construct is also stimulated by Gdf11 when transiently transfected into cultures of HepG2 cells. Stimulation of the reporter in HepG2 cells is enhanced where it contains only the 332 bp Hoxd11 enhancer region VIII upstream or downstream of a luciferase or lacZ reporter. This enhancer contains three elements conserved from fish to mice, one of which has the sequence of a Smad3/4 binding element. Mutation of this motif inhibits the ability of Gdf11 to enhance reporter activity in the HepG2 cell assay. Chromatin immunoprecipitation experiments show direct evidence of Smad2/3 protein binding to the Hoxd11 region VIII enhancer. The action of Gdf11 upon Hoxd11 in HepG2 cells is inhibited, at least in part, by SIS3, a specific inhibitor of Smad3. SIS3 also produces partial inhibition of Hoxd11/lacZ expression in cultured transgenic tailbuds, indicating that Smad3 may play a similar role in the embryonic expression of Hoxd11. Transgenic mouse experiments show that the Smad binding motif is essential for the axial expression of Hoxd11/lacZ reporter in the embryo tailbud, posterior mesoderm and neurectoderm.

Recombineering-based dissection of flanking and paralogous Hox gene functions in mouse reproductive tracts

Hox genes are key regulators of development. In mammals, the study of these genes is greatly confounded by their large number, overlapping functions and interspersed shared enhancers. Here, we describe the use of a novel recombineering strategy to introduce simultaneous frameshift mutations into the flanking Hoxa9, Hoxa10 and Hoxa11 genes, as well as their paralogs on the HoxD cluster. The resulting Hoxa9,10,11 mutant mice displayed dramatic synergistic homeotic transformations of the reproductive tracts, with the uterus anteriorized towards oviduct and the vas deferens anteriorized towards epididymis. The Hoxa9,10,11 mutant mice also provided a genetic setting that allowed the discovery of Hoxd9,10,11 redundant reproductive tract patterning function. Both shared and distinct Hox functions were defined. Hoxd9,10,11 play a crucial role in the regulation of uterine immune function. Non-coding non-polyadenylated RNAs were among the key Hox targets, with dramatic downregulation in mutants. We observed Hox cross-regulation of transcription and splicing. In addition, we observed a surprising anti-dogmatic apparent posteriorization of the uterine epithelium. In caudal regions of the uterus, the normal simple columnar epithelium flanking the lumen was replaced by a pseudostratified transitional epithelium, normally found near the more posterior cervix. These results identify novel molecular functions of Hox genes in the development of the male and female reproductive tracts.

Copy number analysis of 413 isolated talipes equinovarus patients suggests role for transcriptional regulators of early limb development

Talipes equinovarus is one of the most common congenital musculoskeletal anomalies and has a worldwide incidence of 1 in 1000 births. A genetic predisposition to talipes equinovarus is evidenced by the high concordance rate in twin studies and the increased risk to first-degree relatives. Despite the frequency of isolated talipes equinovarus and the strong evidence of a genetic basis for the disorder, few causative genes have been identified. To identify rare and/or recurrent copy number variants, we performed a genome-wide screen for deletions and duplications in 413 isolated talipes equinovarus patients using the Affymetrix 6.0 array. Segregation analysis within families and gene expression in mouse E12.5 limb buds were used to determine the significance of copy number variants. We identified 74 rare, gene-containing copy number variants that were present in talipes equinovarus probands and not present in 759 controls or in the Database of Genomic Variants. The overall frequency of copy number variants was similar between talipes equinovarus patients compared with controls. Twelve rare copy number variants segregate with talipes equinovarus in multiplex pedigrees, and contain the developmentally expressed transcription factors and transcriptional regulators PITX1, TBX4, HOXC13, UTX, CHD (chromodomain protein)1, and RIPPLY2. Although our results do not support a major role for recurrent copy number variations in the etiology of isolated talipes equinovarus, they do suggest a role for genes involved in early embryonic patterning in some families that can now be tested with large-scale sequencing methods.

Forward genetics uncovers Transmembrane protein 107 as a novel factor required for ciliogenesis and Sonic hedgehog signaling

Cilia are dynamic organelles that are essential for a vast array of developmental patterning events, including left-right specification, skeletal formation, neural development, and organogenesis. Despite recent advances in understanding cilia form and function, many key ciliogenesis components have yet to be identified. By using a forward genetics approach, we isolated a novel mutant allele (schlei) of the mouse Transmembrane protein 107 (Tmem107) gene, which we show here is critical for cilia formation and embryonic patterning. Tmem107 is required for normal Sonic hedgehog (Shh) signaling in the neural tube and acts in combination with Gli2 and Gli3 to pattern ventral and intermediate neuronal cell types. schlei mutants also form extra digits, and we demonstrate that Tmem107 acts in the Shh pathway to determine digit number, but not identity, by regulating a subset of Shh target genes. Phenotypically, schlei mutants share several features with other cilia mutants; however, spatial restriction of mutant phenotypes and lack of left-right patterning defects in schlei animals suggest differential requirements for Tmem107 in cilia formation in distinct tissues. Also, in contrast to mutants with complete loss of cilia, schlei mutants retain some function of both Gli activator and repressor forms. Together, these studies identify a previously unknown regulator of ciliogenesis and provide insight into how ciliary factors affect Shh signaling and cilia biogenesis in distinct tissues.

RNA-Seq defines novel genes, RNA processing patterns and enhancer maps for the early stages of nephrogenesis: Hox supergenes

During kidney development the cap mesenchyme progenitor cells both self renew and differentiate into nephrons. The balance between renewal and differentiation determines the final nephron count, which is of considerable medical importance. An important goal is to create a precise genetic definition of the early differentiation of cap mesenchyme progenitors. We used RNA-Seq to transcriptional profile the cap mesenchyme progenitors and their first epithelial derivative, the renal vesicles. The results provide a global view of the changing gene expression program during this key period, defining expression levels for all transcription factors, growth factors, and receptors. The RNA-Seq was performed using two different biochemistries, with one examining only polyadenylated RNA and the other total RNA. This allowed the analysis of noncanonical transcripts, which for many genes were more abundant than standard exonic RNAs. Since a large fraction of enhancers are now known to be transcribed the results also provide global maps of potential enhancers. Further, the RNA-Seq data defined hundreds of novel splice patterns and large numbers of new genes. Particularly striking was the extensive sense/antisense transcription and changing RNA processing complexities of the Hox clusters.

A function for all posterior Hoxd genes during digit development?

BACKGROUND

Four posterior Hoxd genes, from Hoxd13 to Hoxd10, are collectively regulated during the development of tetrapod digits. Besides the well-documented role of Hoxd13, the function of the neighboring genes has been difficult to evaluate due to the close genetic linkage and potential regulatory interferences. We used a combination of five small nested deletions in cis, involving from two to four consecutive genes of the Hoxd13 to Hoxd9 loci, in mice, to evaluate their combined functional importance.

RESULTS

We show that deletions leading to a gain of function of Hoxd13, via regulatory re-allocation, generate abnormal phenotypes, in agreement with the dominant negative role of this gene. We also show that Hoxd10, Hoxd11, and Hoxd12 all seem to play a genuine role in digit development, though less compelling than that of Hoxd13. In contrast, the nearby Hoxd9 contributed no measurable function in digits.

CONCLUSIONS

We conclude that a slight and transient deregulation of Hoxd13 expression can readily affect the relative lengths of limb segments and that all posterior Hoxd genes likely contribute to the final limb morphology. We discuss the difficulty to clearly assess the functional share of individual genes within such a gene family, where closely located neighbors, coding for homologous proteins, are regulated by a unique circuitry and all contribute to shape the distal parts of our appendages.

Shox2 regulates progression through chondrogenesis in the mouse proximal limb

In humans, loss of SHOX gene function is responsible for the mesomelic short stature characteristic of Turner syndrome, Leri-Weill dyschondrosteosis, and Langer dysplasia. In a mouse model of SHOX deficiency, Prrx1-Cre-driven limb-specific deletion of the paralogous gene Shox2 results in severe rhizomelia. In this study, we show that Col2a1-Cre-driven deletion of Shox2 in developing chondrocytes also results in shortening of the stylopodial skeleton (i.e. humerus, femur) and that this rhizomelia is due to precocious chondrocyte maturation and hypertrophy. We demonstrate, using the micromass culture model system, that increased BMP activity triggers accelerated maturation and hypertrophy in Col2a1-Cre Shox2 mutant chondrocytes and we confirm in vivo that elevated transcript levels and expanded expression domains of Bmp2 and 4 are associated with premature formation of the hypertrophic zone in mutant humeri. In micromass cultures of Prrx1-Cre Shox2 mutant limb cells, we find that Shox2 deletion in undifferentiated mesenchymal cells results in increased BMP activity that enhances early chondrogenesis, but is insufficient to provoke chondrocyte maturation and hypertrophy. Similarly, shRNA-mediated Shox2 knockdown in multipotent C3H10T1/2 cells and primary mouse bone marrow mesenchymal stem cells results in spontaneous chondrogenesis in the absence of chondrostimulation, but again fails to induce progression through the later stages of chondrogenic differentiation. Importantly, exogenous BMP supplementation can overcome the block to maturation and hypertrophy caused by Shox2 depletion prior to overt chondrogenesis. Thus, we provide evidence that Shox2 regulates progression through chondrogenesis at two distinct stages – the onset of early differentiation and the transition to maturation and hypertrophy.

Homeobox genes d11-d13 and a13 control mouse autopod cortical bone and joint formation

The molecular mechanisms that govern bone and joint formation are complex, involving an integrated network of signaling pathways and gene regulators. We investigated the role of Hox genes, which are known to specify individual segments of the skeleton, in the formation of autopod limb bones (i.e., the hands and feet) using the mouse mutant synpolydactyly homolog (spdh), which encodes a polyalanine expansion in Hoxd13. We found that no cortical bone was formed in the autopod in spdh/spdh mice; instead, these bones underwent trabecular ossification after birth. Spdh/spdh metacarpals acquired an ovoid shape and developed ectopic joints, indicating a loss of long bone characteristics and thus a transformation of metacarpals into carpal bones. The perichondrium of spdh/spdh mice showed abnormal morphology and decreased expression of Runt-related transcription factor 2 (Runx2), which was identified as a direct Hoxd13 transcriptional target. Hoxd11-/-Hoxd12-/-Hoxd13-/- triple-knockout mice and Hoxd13-/-Hoxa13+/- mice exhibited similar but less severe defects, suggesting that these Hox genes have similar and complementary functions and that the spdh allele acts as a dominant negative. This effect was shown to be due to sequestration of other polyalanine-containing transcription factors by the mutant Hoxd13 in the cytoplasm, leading to their degradation. These data indicate that Hox genes not only regulate patterning but also directly influence bone formation and the ossification pattern of bones, in part via Runx2.

Hox11 paralogous genes are required for formation of wrist and ankle joints and articular surface organization

Limb skeletal elements are connected by distinct synovial joints, but the mechanisms regulating joint formation, diversity, and organization remain unclear. Previous studies showed that Hox11 mouse mutants have severe developmental defects in radius and ulna and tibia and fibula, but wrist and ankle joint formation and characteristics were not examined in detail. We now find that E11.5 and E12.5 triple Hox11aaccdd mutants exhibit a significant reduction in prospective carpal and tarsal mesenchyme. Although the mesenchyme became segmented into individual carpal and tarsal skeletal elements with further development, the elements were ill defined and the more proximal elements (radiale, ulnare, talus, and calcaneous) actually underwent involution and/or fusion. Wild-type carpal and tarsal elements displayed a thick articulating superficial zone at their outer perimeter that expressed genes typical of developing joint interzones and articulating cells, including Gdf5, Erg, Gli3, collagen IIA, and lubricin, and defined each element anatomically. In mutant wrists and ankles, the superficial zone around each element was thin and ill defined, and expression of several of those genes was low and often interrupted. These and other data provide novel and clear evidence that Hox11 paralogous genes regulate wrist and ankle joint organization and are essential for establishing carpal and tarsal element boundary and maintaining their articulating surface tissue.

Control of vertebrate Hox clusters by remote and global cis-acting regulatory sequences

Despite apparently shared structural organisation and functional roles, vertebrate Hox genes are controlled by regulatory mechanisms rather distinct from those of the prototypic Drosophila Antennapedia (ANT-C) and Bithorax (BX-C) Complexes. If individual regulatory modules have been shown to recapitulate specific Hox expression patterns, other experimental studies underscore that vertebrate Hox clusters are controlled in many of their functions in a global manner, through distinct mechanisms. We will discuss the different models that have been proposed to account for these global regulatory modes. In this context, the studies of the regulation of the HoxD complex during limb development highlighted the role of global regulatory elements and the different mechanisms associated to transform a structural organisation into distinct temporal and spatial expression domains. We will further discuss how these mechanisms may have benefited from the structure of the vertebrate homeotic clusters and reciprocally contribute to shape their evolution towards an increased level of organisation and compaction.