The Dlx genes as clues to vertebrate genomics and craniofacial evolution

Transcriptome and open chromatin analysis reveals the process of myocardial cell development and key pathogenic target proteins in Long QT syndrome type 7

OBJECTIVE

Long QT syndrome type 7 (Andersen-Tawil syndrome, ATS), which is caused by KCNJ2 gene mutation, often leads to ventricular arrhythmia, periodic paralysis and skeletal malformations. The development, differentiation and electrophysiological maturation of cardiomyocytes (CMs) changes promote the pathophysiology of Long QT syndrome type 7(LQT7). We aimed to specifically reproduce the ATS disease phenotype and study the pathogenic mechanism.

METHODS AND RESULTS

We established a cardiac cell model derived from human induced pluripotent stem cells (hiPSCs) to the phenotypes and electrophysiological function, and the establishment of a human myocardial cell model that specifically reproduces the symptoms of ATS provides a reliable platform for exploring the mechanism of this disease or potential drugs. The spontaneous pulsation rate of myocardial cells in the mutation group was significantly lower than that in the repair CRISPR group, the action potential duration was prolonged, and the Kir2.1 current of the inward rectifier potassium ion channel was decreased, which is consistent with the clinical symptoms of ATS patients. Only ZNF528, a chromatin-accessible TF related to pathogenicity, was continuously regulated beginning from the cardiac mesodermal precursor cell stage (day 4), and continued to be expressed at low levels, which was identified by WGCNA method and verified with ATAC-seq data in the mutation group. Subsequently, it indicated that seven pathways were downregulated (all p < 0.05) by used single sample Gene Set Enrichment Analysis to evaluate the overall regulation of potassium-related pathways enriched in the transcriptome and proteome of late mature CMs. Among them, the three pathways (GO: 0008076, GO: 1990573 and GO: 0030007) containing the mutated gene KCNJ2 is involved that are related to the whole process by which a potassium ion enters the cell via the inward rectifier potassium channel to exert its effect were inhibited. The other four pathways are related to regulation of the potassium transmembrane pathway and sodium:potassium exchange ATPase (p < 0.05). ZNF528 small interfering (si)-RNA was applied to hiPSC-derived cardiomyocytes for CRISPR group to explore changes in potassium ion currents and growth and development related target protein levels that affect disease phenotype. Three consistently downregulated proteins (KCNJ2, CTTN and ATP1B1) associated with pathogenicity were verificated through correlation and intersection analysis.

CONCLUSION

This study uncovers TFs and target proteins related to electrophysiology and developmental pathogenicity in ATS myocardial cells, obtaining novel targets for potential therapeutic candidate development that does not rely on gene editing.

DLX gene expression in the developing chick pharyngeal arches and relationship to endothelin signaling and avian jaw patterning

BACKGROUND

A hinged jaw that articulates with the skull base is a striking feature of the vertebrate head and has been greatly modified between, and within, vertebrate classes. Genes belonging to the DLX homeobox family are conserved mediators of local signaling pathways that distinguish the dorsal and ventral aspects of the first pharyngeal arch. Specifically, a subset of DLX genes are expressed in the cranial neural crest-derived mandibular ectomesenchyme in response to ventral endothelin signaling, an important step that confers the first arch with maxillary and mandibular identities. Downstream targets of DLX genes then execute the morphogenetic processes that lead to functional jaws. Identifying lineage-specific variations in DLX gene expression and the regulatory networks downstream of DLX action is necessary to understand how different kinds of jaws evolved.

RESULTS

Here, we describe and compare the expression of all six DLX genes in the chick pharyngeal arches, focusing on the period of active patterning in the first arch. Disruption of endothelin signaling results in the down-regulation of ventral-specific DLX genes and confirms their functional role in avian jaw patterning.

CONCLUSIONS

This expression resource will be important for comparative embryology and for identifying synexpression groups of DLX-regulated genes in the chick.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

Ectomesenchymal Six1 controls mandibular skeleton formation

Craniofacial development requires intricate cooperation between multiple transcription factors and signaling pathways. Six1 is a critical transcription factor regulating craniofacial development. However, the exact function of Six1 during craniofacial development remains elusive. In this study, we investigated the role of Six1 in mandible development using a Six1 knockout mouse model (Six1 ^-/- ) and a cranial neural crest-specific, Six1 conditional knockout mouse model (Six1 ^f/f ; Wnt1-Cre). The Six1 ^-/- mice exhibited multiple craniofacial deformities, including severe microsomia, high-arched palate, and uvula deformity. Notably, the Six1 ^f/f ; Wnt1-Cre mice recapitulate the microsomia phenotype of Six1 ^-/- mice, thus demonstrating that the expression of Six1 in ectomesenchyme is critical for mandible development. We further showed that the knockout of Six1 led to abnormal expression of osteogenic genes within the mandible. Moreover, the knockdown of Six1 in C3H10 T1/2 cells reduced their osteogenic capacity in vitro. Using RNA-seq, we showed that both the loss of Six1 in the E18.5 mandible and Six1 knockdown in C3H10 T1/2 led to the dysregulation of genes involved in embryonic skeletal development. In particular, we showed that Six1 binds to the promoter of Bmp4, Fat4, Fgf18, and Fgfr2, and promotes their transcription. Collectively, our results suggest that Six1 plays a critical role in regulating mandibular skeleton formation during mouse embryogenesis.

The pseudobranch of jawed vertebrates is a mandibular arch-derived gill

The pseudobranch is a gill-like epithelial elaboration that sits behind the jaw of most fishes. This structure was classically regarded as a vestige of the ancestral gill arch-like condition of the gnathostome jaw. However, more recently, hypotheses of jaw evolution by transformation of a gill arch have been challenged, and the pseudobranch has alternatively been considered a specialised derivative of the second (hyoid) pharyngeal arch. Here, we demonstrate in the skate (Leucoraja erinacea) that the pseudobranch does, in fact, derive from the mandibular arch, and that it shares gene expression features and cell types with gills. We also show that the skate mandibular arch pseudobranch is supported by a spiracular cartilage that is patterned by a shh-expressing epithelial signalling centre. This closely parallels the condition seen in the gill arches, where cartilaginous appendages called branchial rays, which support the respiratory lamellae of the gills, are patterned by a shh-expressing gill arch epithelial ridge. Together with similar discoveries in zebrafish, our findings support serial homology of the pseudobranch and gills, and an ancestral origin of gill arch-like anatomical features from the gnathostome mandibular arch.

Summary: The skate pseudobranch is a gill serial homologue and reveals the ancestral gill arch-like nature of the jawed vertebrate mandibular arch.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Morphological stasis masks ecologically divergent coral species on tropical reefs

Coral reefs are the epitome of species diversity, yet the number of described scleractinian coral species, the framework-builders of coral reefs, remains moderate by comparison. DNA sequencing studies are rapidly challenging this notion by exposing a wealth of undescribed diversity, but the evolutionary and ecological significance of this diversity remains largely unclear. Here, we present an annotated genome for one of the most ubiquitous corals in the Indo-Pacific (Pachyseris speciosa) and uncover, through a comprehensive genomic and phenotypic assessment, that it comprises morphologically indistinguishable but ecologically divergent lineages. Demographic modeling based on whole-genome resequencing indicated that morphological crypsis (across micro- and macromorphological traits) was due to ancient morphological stasis rather than recent divergence. Although the lineages occur sympatrically across shallow and mesophotic habitats, extensive genotyping using a rapid molecular assay revealed differentiation of their ecological distributions. Leveraging "common garden" conditions facilitated by the overlapping distributions, we assessed physiological and quantitative skeletal traits and demonstrated concurrent phenotypic differentiation. Lastly, spawning observations of genotyped colonies highlighted the potential role of temporal reproductive isolation in the limited admixture, with consistent genomic signatures in genes related to morphogenesis and reproduction. Overall, our findings demonstrate the presence of ecologically and phenotypically divergent coral species without substantial morphological differentiation and provide new leads into the potential mechanisms facilitating such divergence. More broadly, they indicate that our current taxonomic framework for reef-building corals may be scratching the surface of the ecologically relevant diversity on coral reefs, consequently limiting our ability to protect or restore this diversity effectively.

DLX5 and HOXC8 enhance the chondrogenic differentiation potential of stem cells from apical papilla via LINC01013

BACKGROUND

Mesenchymal stem cell (MSC)-based cartilage tissue regeneration is a treatment with great potential. How to enhance the MSC chondrogenic differentiation is a key issue involved in cartilage formation. In the present study, we seek to expound the phenotypes and mechanisms of DLX5 in chondrogenic differentiation function in MSCs.

METHODS

Stem cells from apical papilla (SCAPs) were used. The Alcian Blue staining, pellet culture system, and cell transplantation in rabbit knee cartilage defect were used to evaluate the chondrogenic differentiation function of MSCs. Western blot, real-time RT-PCR, and ChIP assays were used to evaluate the molecular mechanisms.

RESULTS

DLX5 and HOXC8 expressions were upregulated during chondrogenic differentiation. In vitro results showed that DLX5 and HOXC8 enhanced the expression of chondrogenic markers including collagen II (COL2), collagen V (COL5), and sex-determining region Y box protein 9 (SOX9) and promoted the chondrogenic differentiation and the formation of cartilage clumps in the pellet culture system. Mechanically, DLX5 and HOXC8 formed protein complexes and negatively regulated the LncRNA, LINC01013, via directly binding its promoter. In vivo transplantation experiment showed that DLX5 and HOXC8 could restore the cartilage defect in the rabbit knee model. In addition, knock-down of LINC01013 enhanced the chondrogenic differentiation of SCAPs.

CONCLUSIONS

In conclusion, DLX5 and HOXC8 enhance the chondrogenic differentiation abilities of SCAPs by negatively regulating LINC01013 in SCAPs, and provided the potential target for promoting cartilage tissue regeneration.

Identification and validation of prognosis-related DLX5 methylation as an epigenetic driver in myeloid neoplasms

The deregulated DLX gene family members DLX1/2/3/4/5/6 (DLXs) caused by DNA methylation has been demonstrated in various cancers with therapeutic target value. However, the potential role of DLXs methylation in myeloid neoplasms such as acute myeloid leukemia (AML) and myelodysplastic syndromes (MDS) remains to be elucidated. Clinical significance of DLXs methylation/expression was analyzed in patient with AML and MDS. The functional roles of DLXs were determined in vitro. In the identification stage, we found that lower DLX5 expression was correlated with prognosis in AML among all DLXs analyzed by The Cancer Genome Atlas datasets. In the validation stage, we revealed that reduced DLX5 expression was frequently occurred, and was also correlated with promoter hypermethylation in AML evaluated by targeted bisulfite sequencing. Epigenetic studies also showed that DLX5 promoter DNA methylation was associated with its expression. By quantitative polymerase chain reaction, we also validated that DLX5 hypermethylation was frequent event in both AML and MDS, and also correlated with MDS transformation to leukemia. Moreover, DLX5 hypermethylation was associated with lower rate of complete remission and shorter time of leukemia‐free/overall survival, and was also confirmed by Logistic/Cox regression analysis. Functional studies revealed the antiproliferative and pro‐apoptotic effects of DLX5 in MDS‐derived AML cell‐line SKM‐1. Finally, bioinformatics analysis demonstrated that DLX5 functioned in leukemogenesis may be through the association with PI3K/Akt signaling pathway. Collectively, our findings demonstrated that DLX5 methylation, negatively correlated DLX5 expression, was a potential prognostic and predictive indicator in patients with AML and MDS, which could also act as an epigenetic driver in myeloid neoplasms.

Development of hypobranchial muscles with special reference to the evolution of the vertebrate neck

BACKGROUND

The extant vertebrates include cyclostomes (lamprey and hagfish) and crown gnathostomes (jawed vertebrates), but there are various anatomical disparities between these two groups. Conspicuous in the gnathostomes is the neck, which occupies the interfacial domain between the head and trunk, including the occipital part of the cranium, the shoulder girdle, and the cucullaris and hypobranchial muscles (HBMs). Of these, HBMs originate from occipital somites to form the ventral pharyngeal and neck musculature in gnathostomes. Cyclostomes also have HBMs on the ventral pharynx, but lack the other neck elements, including the occipital region, the pectoral girdle, and cucullaris muscles. These anatomical differences raise questions about the evolution of the neck in vertebrates.

RESULTS

In this study, we observed developing HBMs as a basis for comparison between the two groups and show that the arrangement of the head-trunk interface in gnathostomes is distinct from that of lampreys. Our comparative analyses reveal that, although HBM precursors initially pass through the lateral side of the pericardium in both groups, the relative positions of the pericardium withrespect to the pharyngeal arches differ between the two, resulting in diverse trajectories of HBMs in gnathostomes and lampreys.

CONCLUSIONS

We suggest that a heterotopic rearrangement of early embryonic components, including the pericardium and pharyngeal arches, may have played a fundamental role in establishing the gnathostome HBMs, which would also have served as the basis for neck formation in the jawed vertebrate lineage.

Flightless-I governs cell fate by recruiting the SUMO isopeptidase SENP3 to distinct HOX genes

Background

Despite recent studies on the role of ubiquitin-related SUMO modifier in cell fate decisions, our understanding on precise molecular mechanisms of these processes is limited. Previously, we established that the SUMO isopeptidase SENP3 regulates chromatin assembly of the MLL1/2 histone methyltransferase complex at distinct HOX genes, including the osteogenic master regulator DLX3. A comprehensive mechanism that regulates SENP3 transcriptional function was not understood.

Results

Here, we identified flightless-I homolog (FLII), a member of the gelsolin family of actin-remodeling proteins, as a novel regulator of SENP3. We demonstrate that FLII is associated with SENP3 and the MLL1/2 complex. We further show that FLII determines SENP3 recruitment and MLL1/2 complex assembly on the DLX3 gene. Consequently, FLII is indispensible for H3K4 methylation and proper loading of active RNA polymerase II at this gene locus. Most importantly, FLII-mediated SENP3 regulation governs osteogenic differentiation of human mesenchymal stem cells.

Conclusion

Altogether, these data reveal a crucial functional interconnection of FLII with the sumoylation machinery that converges on epigenetic regulation and cell fate determination.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0122-8) contains supplementary material, which is available to authorized users.

Developmental transcriptional regulation by SUMOylation, an evolving field

SUMOylation is a reversible post-translational protein modification that affects the intracellular localization, stability, activity, and interactions of its protein targets. The SUMOylation pathway influences several nuclear and cytoplasmic processes. The expression of many genes, in particular those involved in development is finely tuned in space and time by several groups of proteins. There is growing evidence that transcriptional regulation mechanisms involve direct SUMOylation of transcriptional-related proteins such as initiation and elongation factors, and subunits of chromatin modifier and remodeling complexes originally described as members of the trithorax and Polycomb groups in Drosophila. Therefore, it is being unveiled that SUMOylation has a role in both, gene silencing and gene activation mechanisms. The goal of this review is to discuss the information on how SUMO modification in components of these multi-subunit complexes may have an effect in genome architecture and function and, therefore, in the regulation of gene expression in time and space.

Sp7/Osterix Is Restricted to Bone-Forming Vertebrates where It Acts as a Dlx Co-factor in Osteoblast Specification

In extant species, bone formation is restricted to vertebrate species. Sp7/Osterix is a key transcriptional determinant of bone-secreting osteoblasts. We performed Sp7 chromatin immunoprecipitation sequencing analysis identifying a large set of predicted osteoblast enhancers and validated a subset of these in cell culture and transgenic mouse assays. Sp family members bind GC-rich target sequences through their zinc finger domain. Several lines of evidence suggest that Sp7 acts differently, engaging osteoblast targets in Dlx-containing regulatory complexes bound to AT-rich motifs. Amino acid differences in the Sp7 zinc finger domain reduce Sp7's affinity for the Sp family consensus GC-box target; Dlx5 binding maps to this domain of Sp7. The data support a model in which Dlx recruitment of Sp7 to osteoblast enhancers underlies Sp7-directed osteoblast specification. Because an Sp7-like zinc finger variant is restricted to vertebrates, the emergence of an Sp7 member within the Sp family was likely closely coupled to the evolution of bone-forming vertebrates.

Evidence from cyclostomes for complex regionalization of the ancestral vertebrate brain

The vertebrate brain is highly complex, but its evolutionary origin remains elusive. Because of the absence of certain developmental domains generally marked by the expression of regulatory genes, the embryonic brain of the lamprey, a jawless vertebrate, had been regarded as representing a less complex, ancestral state of the vertebrate brain. Specifically, the absence of a Hedgehog- and Nkx2.1-positive domain in the lamprey subpallium was thought to be similar to mouse mutants in which the suppression of Nkx2-1 leads to a loss of the medial ganglionic eminence. Here we show that the brain of the inshore hagfish (Eptatretus burgeri), another cyclostome group, develops domains equivalent to the medial ganglionic eminence and rhombic lip, resembling the gnathostome brain. Moreover, further investigation of lamprey larvae revealed that these domains are also present, ruling out the possibility of convergent evolution between hagfish and gnathostomes. Thus, brain regionalization as seen in crown gnathostomes is not an evolutionary innovation of this group, but dates back to the latest vertebrate ancestor before the divergence of cyclostomes and gnathostomes more than 500 million years ago.

Osteogenesis Imperfecta Type I Caused by COL1A1 Deletions

Osteogenesis imperfecta (OI) type I is usually caused by COL1A1 stop or frameshift mutations, leading to COL1A1 haploinsufficiency. Here we report on 12 individuals from 5 families who had OI type I due to an unusual cause—heterozygous deletions of the entire COL1A1 gene. The deletions were initially detected by semiconductor-based sequencing of genomic DNA and confirmed by quantitative PCR. Array comparative genomic hybridization in DNA of the index patient in each family showed that deletion size varied from 18.5 kb to 2.23 Mb between families, encompassing between 1 and 47 genes (COL1A1 included). The skeletal phenotype of the affected individuals was similar to that of patients with haploinsufficiency caused by COL1A1 stop or frameshift mutations. However, one individual with a deletion that included also DLX3 and DLX4 had tooth discoloration and bone fragility. Three individuals from 2 families had deletions that included also CACNA1G, and these individuals had learning disabilities. These features are not usually observed in COL1A1 haploinsufficiency, but are in accordance with previously described individuals in whom deletions included the same genes. In summary, we found deletions of COL1A1 in 5 out of 161 families (3 %) with OI type I that were evaluated. Deletions encompassing not only COL1A1 but also neighboring genes can lead to contiguous gene syndromes that may include dental involvement and learning disability.

Molecular phylogeny of four homeobox genes from the purple sea star Pisaster ochraceus

Homeobox genes cloned from the purple sea star Pisaster ochraceus (Phylum Echinodermata/Class Asteroidea) were used along with related sequences available from members of other representative animal phyla to generate molecular phylogenies for Distal-less/Dlx, Hox5, Hox7, and Hox9/10 homeobox genes. Phylogenetic relationships were inferred based on the predicted 60 amino acid homeodomain, using amino acid (AA) and nucleotide (NT) models as well as the recently developed codon substitution models of sequence evolution. The resulting phylogenetic trees were mostly congruent with the consensus species-tree, grouping these newly identified genes with those isolated from other Asteroidea. This analysis also allowed a preliminary comparison of the performance of codon models with that of NT and AA evolutionary models in the inference of homeobox phylogeny. We found that, overall, the NT models displayed low reliability in recovering major clades at the Superphylum/Phylum level, and that codon models were slightly more dependable than AA models. Remarkably, in the majority of cases, codon substitution models seemed to outperform both AA and NT models at both the Class level and homeobox paralogy-group level of classification.

Drug discovery for alopecia: gone today, hair tomorrow

INTRODUCTION

Hair loss or alopecia affects the majority of the population at some time in their life, and increasingly, sufferers are demanding treatment. Three main types of alopecia (androgenic [AGA], areata [AA] and chemotherapy-induced [CIA]) are very different, and have their own laboratory models and separate drug-discovery efforts.

AREAS COVERED

In this article, the authors review the biology of hair, hair follicle (HF) cycling, stem cells and signaling pathways. AGA, due to dihydrotesterone, is treated by 5-α reductase inhibitors, androgen receptor blockers and ATP-sensitive potassium channel-openers. AA, which involves attack by CD8(+)NK group 2D-positive (NKG2D(+)) T cells, is treated with immunosuppressives, biologics and JAK inhibitors. Meanwhile, CIA is treated by apoptosis inhibitors, cytokines and topical immunotherapy.

EXPERT OPINION

The desire to treat alopecia with an easy topical preparation is expected to grow with time, particularly with an increasing aging population. The discovery of epidermal stem cells in the HF has given new life to the search for a cure for baldness. Drug discovery efforts are being increasingly centered on these stem cells, boosting the hair cycle and reversing miniaturization of HF. Better understanding of the molecular mechanisms underlying the immune attack in AA will yield new drugs. New discoveries in HF neogenesis and low-level light therapy will undoubtedly have a role to play.

A gene expression map of the larval Xenopus laevis head reveals developmental changes underlying the evolution of new skeletal elements

The morphology of the vertebrate head skeleton is highly plastic, with the number, size, shape, and position of its components varying dramatically between groups. While this evolutionary flexibility has been key to vertebrate success, its developmental and genetic bases are poorly understood. The larval head skeleton of the frog Xenopus laevis possesses a unique combination of ancestral tetrapod features and anuran-specific novelties. We built a detailed gene expression map of the head mesenchyme in X. laevis during early larval development, focusing on transcription factor families with known functions in vertebrate head skeleton development. This map was then compared to homologous gene expression in zebrafish, mouse, and shark embryos to identify conserved and evolutionarily flexible aspects of vertebrate head skeleton development. While we observed broad conservation of gene expression between X. laevis and other gnathostomes, we also identified several divergent features that correlate to lineage-specific novelties. We noted a conspicuous change in dlx1/2 and emx2 expression in the second pharyngeal arch, presaging the differentiation of the reduced dorsal hyoid arch skeletal element typical of modern anamniote tetrapods. In the first pharyngeal arch we observed a shift in the expression of the joint inhibitor barx1, and new expression of the joint marker gdf5, shortly before skeletal differentiation. This suggests that the anuran-specific infrarostral cartilage evolved by partitioning of Meckel's cartilage with a new paired joint. Taken together, these comparisons support a model in which early patterning mechanisms divide the vertebrate head mesenchyme into a highly conserved set of skeletal precursor populations. While subtle changes in this early patterning system can affect skeletal element size, they do not appear to underlie the evolution of new joints or cartilages. In contrast, later expression of the genes that regulate skeletal element differentiation can be clearly linked to the evolution of novel skeletal elements. We posit that changes in the expression of downstream regulators of skeletal differentiation, like barx1 and gdf5, is one mechanism by which head skeletal element number and articulation are altered during evolution.