Hox genes play a critical role in the development of the vertebrate axis and limbs, and previous studies have implicated them in the specification of positional identity, the control of growth, and the timing of differentiation. Axolotl limbs offer an opportunity to distinguish these alternatives because the sequence of skeletal differentiation is reversed along the anterior-posterior axis relative to that of other tetrapods. We report that during early limb development, expression patterns of HoxD genes in axolotls resemble those in amniotes and anuran amphibians. At later stages, the anterior boundary of Hoxd-11 expression is conserved with respect to morphological landmarks, but there is no anterior-distal expansion of the posterior domain of Hoxd-11 expression similar to that observed in mice and chicks. Since axolotls do not form an expanded paddle-like handplate prior to digit differentiation, we suggest that anterior expansion of expression in higher vertebrates is linked to the formation of the handplate, but is clearly not necessary for digit differentiation. We also show that the 5' HoxD genes are reexpressed during limb regeneration. The change in the expression pattern of Hoxd-11 during the course of regeneration is consistent with the hypothesis that the distal tip of the regenerate is specified first, followed by intercalation of intermediate levels of the pattern. Both Hoxd-8 and Hoxd-10 are expressed in non-regenerating wounds, but Hoxd-11 is specific for regeneration. It is also expressed in the posterior half of nerve-induced supernumerary outgrowths.

Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

To connect human biology to fish biomedical models, we sequenced the genome of spotted gar (Lepisosteus oculatus), whose lineage diverged from teleosts before teleost genome duplication (TGD). The slowly evolving gar genome has conserved in content and size many entire chromosomes from bony vertebrate ancestors. Gar bridges teleosts to tetrapods by illuminating the evolution of immunity, mineralization and development (mediated, for example, by Hox, ParaHox and microRNA genes). Numerous conserved noncoding elements (CNEs; often cis regulatory) undetectable in direct human-teleost comparisons become apparent using gar: functional studies uncovered conserved roles for such cryptic CNEs, facilitating annotation of sequences identified in human genome-wide association studies. Transcriptomic analyses showed that the sums of expression domains and expression levels for duplicated teleost genes often approximate the patterns and levels of expression for gar genes, consistent with subfunctionalization. The gar genome provides a resource for understanding evolution after genome duplication, the origin of vertebrate genomes and the function of human regulatory sequences.