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.

Isl1 mediates mesenchymal expansion in the developing external genitalia via regulation of Bmp4, Fgf10 and Wnt5a.

Genital malformations are among the most common human birth defects, and both genetic and environmental factors can contribute to these malformations. Development of the external genitalia in mammals relies on complex signaling networks, and disruption of these signaling pathways can lead to genital defects. Islet-1 (ISL1), a member of the LIM/Homeobox family of transcription factors, has been identified as a major susceptibility gene for classic bladder exstrophy in humans, a common form of the bladder exstrophy-epispadias complex (BEEC), and is implicated in a role in urinary tract development. We report that deletion of Isl1 from the genital mesenchyme in mice led to hypoplasia of the genital tubercle and prepuce, with an ectopic urethral opening and epispadias-like phenotype. These mice also developed hydroureter and hydronephrosis. Identification of ISL1 transcriptional targets via ChIP-Seq and expression analyses revealed that Isl1 regulates several important signaling pathways during embryonic genital development, including the BMP, WNT, and FGF cascades. An essential function of Isl1 during development of the external genitalia is to induce Bmp4-mediated apoptosis in the genital mesenchyme. Together, these studies demonstrate that Isl1 plays a critical role during development of the external genitalia and forms the basis for a greater understanding of the molecular mechanisms underlying the pathogenesis of BEEC and urinary tract defects in humans.

PITX1 promotes chondrogenesis and myogenesis in mouse hindlimbs through conserved regulatory targets.

The PITX1 transcription factor is expressed during hindlimb development, where it plays a critical role in directing hindlimb growth and the specification of hindlimb morphology. While it is known that PITX1 regulates hindlimb formation, in part, through activation of the Tbx4 gene, other transcriptional targets remain to be elucidated. We have used a combination of ChIP-seq and RNA-seq to investigate enhancer regions and target genes that are directly regulated by PITX1 in embryonic mouse hindlimbs. In addition, we have analyzed PITX1 binding sites in hindlimbs of Anolis lizards to identify ancient PITX1 regulatory targets. We find that PITX1-bound regions in both mouse and Anolis hindlimbs are strongly associated with genes implicated in limb and skeletal system development. Gene expression analyses reveal a large number of misexpressed genes in the hindlimbs of Pitx1-/- mouse embryos. By intersecting misexpressed genes with genes that have neighboring mouse PITX1 binding sites, we identified 440 candidate targets of PITX1. Of these candidates, 68 exhibit ultra-conserved PITX1 binding events that are shared between mouse and Anolis hindlimbs. Among the ancient targets of PITX1 are important regulators of cartilage and skeletal muscle development, including Sox9 and Six1. Our data suggest that PITX1 promotes chondrogenesis and myogenesis in the hindlimb by direct regulation of several key members of the cartilage and muscle transcriptional networks.

Appendages and gene regulatory networks: Lessons from the limbless.

Among squamate reptiles, dozens of lineages have independently evolved complete or partial limb reduction. This remarkable convergence of limbless and limb-reduced phenotypes provides multiple natural replicates of different ages to explore the evolution and development of the vertebrate limb and the gene regulatory network that controls its formation. The most successful and best known of the limb-reduced squamates are snakes, which evolved a limb-reduced body form more than 100 million years ago. Recent studies have revealed the unexpected finding that many ancient limb enhancers are conserved in the genomes of snakes. Analyses in limbed animals show that many of these limb enhancers are also active during development of the phallus, suggesting that these enhancers may have been retained in snakes due their importance in regulating transcription in the external genitalia. This hypothesis is substantiated by functional tests of snake enhancers, which demonstrate that snake enhancer elements have lost limb function while retaining genital enhancer function. The large degree of overlap in the gene regulatory networks deployed during limb and phallus development may act to constrain the divergence of shared gene network components and the evolution of appendage morphology. Future studies will reveal whether limb regulatory elements have undergone similar functional changes in other lineages of limb-reduced squamates.

Pitx1 broadly associates with limb enhancers and is enriched on hindlimb cis-regulatory elements.

Extensive functional analyses have demonstrated that the pituitary homeodomain transcription factor Pitx1 plays a critical role in specifying hindlimb morphology in vertebrates. However, much less is known regarding the target genes and cis-regulatory elements through which Pitx1 acts. Earlier studies suggested that the hindlimb transcription factors Tbx4, HoxC10, and HoxC11 might be transcriptional targets of Pitx1, but definitive evidence for direct regulatory interactions has been lacking. Using ChIP-Seq on embryonic mouse hindlimbs, we have pinpointed the genome-wide location of Pitx1 binding sites during mouse hindlimb development and identified potential gene targets for Pitx1. We determined that Pitx1 binding is significantly enriched near genes involved in limb morphogenesis, including Tbx4, HoxC10, and HoxC11. Notably, Pitx1 is bound to the previously identified HLEA and HLEB hindlimb enhancers of the Tbx4 gene and to a newly identified Tbx2 hindlimb enhancer. Moreover, Pitx1 binding is significantly enriched on hindlimb relative to forelimb-specific cis-regulatory features that are differentially marked by H3K27ac. However, our analysis revealed that Pitx1 also strongly associates with many functionally verified limb enhancers that exhibit similar levels of activity in the embryonic mesenchyme of forelimbs and hindlimbs. We speculate that Pitx1 influences hindlimb morphology both through the activation of hindlimb-specific enhancers as well as through the hindlimb-specific modulation of enhancers that are active in both sets of limbs.

Conserved regulation of hoxc11 by pitx1 in Anolis lizards.

Anolis lizards are an emerging model system for the study of limb development and evolution, but very little is known concerning the regulatory interactions that control limb patterning differences among Anolis species or what regulatory interactions are deeply conserved between Anolis and other tetrapod groups. Here we report the establishment of an embryonic limb micromass culture system that enables functional studies of forelimb and hindlimb gene regulatory networks in Anolis. Characterization of this culture system demonstrated that embryonic forelimb and hindlimb micromasses from different Anolis species are easy to sustain in culture for weeks, and the expression of forelimb and hindlimb-specific gene expression patterns are maintained for at least 8 days in culture. We tested the ability of this system to explore regulatory linkages between transcription factors and their putative target genes through the ectopic expression of a hindlimb-specific transcription factor, pitx1, in forelimb micromasses. We found that pitx1 expression in forelimb cells is sufficient to strongly induce the expression of hoxc11, a gene that normally exhibits hindlimb-restricted expression.

Somitogenesis in the anole lizard and alligator reveals evolutionary convergence and divergence in the amniote segmentation clock.

The axial skeleton is a defining feature of vertebrates and is patterned during somitogenesis. Cyclically expressed members of the notch and other signaling pathways, described as the 'segmentation clock', regulate the formation of somite boundaries. Comparisons among vertebrate model systems have revealed fundamental shifts in the regulation of expression among critical genes in the notch pathway. However, insights into the evolution of these expression differences have been limited by the lack of information from non-avian reptiles. We analyzed the segmentation clock of the first Lepidosaurian reptile sequenced, the green anole lizard, Anolis carolinensis, for comparison with avian and mammalian models. Using genomic sequence, RNA-Seq transcriptomic data, and in situ hybridization analysis of somite-stage embryos, we carried out comparative analyses of key genes and found that the anole segmentation clock displays features common to both amniote and anamniote vertebrates. Shared features with anamniotes, represented by Xenopus laevis and Danio rerio, include an absence of lunatic fringe (lfng) expression within the presomitic mesoderm (PSM), a hes6a gradient in the PSM not observed in the chicken or mouse, and EGF repeat structure of the divergent notch ligand, dll3. The anole and mouse share cycling expression of dll1 ligand in the PSM. To gain insight from an Archosaurian reptile, we analysed LFNG and DLL1 expressions in the American alligator. LFNG expression was absent in the alligator PSM, like the anole but unlike the chicken. In contrast, DLL1 expression does not cycle in the PSM of the alligator, similar to the chicken but unlike the anole. Thus, our analysis yields novel insights into features of the segmentation clock that are evolutionarily basal to amniotes versus those that are specific to mammals, Lepidosaurian reptiles, or Archosaurian reptiles.

Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species.

Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei - a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates.

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1. The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB, that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB, suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature.

Comparative Genomics Reveals Accelerated Evolution in Conserved Pathways during the Diversification of Anole Lizards.

Squamates include all lizards and snakes, and display some of the most diverse and extreme morphological adaptations among vertebrates. However, compared with birds and mammals, relatively few resources exist for comparative genomic analyses of squamates, hampering efforts to understand the molecular bases of phenotypic diversification in such a speciose clade. In particular, the ∼400 species of anole lizard represent an extensive squamate radiation. Here, we sequence and assemble the draft genomes of three anole species-Anolis frenatus, Anolis auratus, and Anolis apletophallus-for comparison with the available reference genome of Anolis carolinensis. Comparative analyses reveal a rapid background rate of molecular evolution consistent with a model of punctuated equilibrium, and strong purifying selection on functional genomic elements in anoles. We find evidence for accelerated evolution in genes involved in behavior, sensory perception, and reproduction, as well as in genes regulating limb bud development and hindlimb specification. Morphometric analyses of anole fore and hindlimbs corroborated these findings. We detect signatures of positive selection across several genes related to the development and regulation of the forebrain, hormones, and the iguanian lizard dewlap, suggesting molecular changes underlying behavioral adaptations known to reinforce species boundaries were a key component in the diversification of anole lizards.

Molecular shifts in limb identity underlie development of feathered feet in two domestic avian species.

Birds display remarkable diversity in the distribution and morphology of scales and feathers on their feet, yet the genetic and developmental mechanisms governing this diversity remain unknown. Domestic pigeons have striking variation in foot feathering within a single species, providing a tractable model to investigate the molecular basis of skin appendage differences. We found that feathered feet in pigeons result from a partial transformation from hindlimb to forelimb identity mediated by cis-regulatory changes in the genes encoding the hindlimb-specific transcription factor Pitx1 and forelimb-specific transcription factor Tbx5. We also found that ectopic expression of Tbx5 is associated with foot feathers in chickens, suggesting similar molecular pathways underlie phenotypic convergence between these two species. These results show how changes in expression of regional patterning genes can generate localized changes in organ fate and morphology, and provide viable molecular mechanisms for diversity in hindlimb scale and feather distribution.

Shared Enhancer Activity in the Limbs and Phallus and Functional Divergence of a Limb-Genital cis-Regulatory Element in Snakes.

The amniote phallus and limbs differ dramatically in their morphologies but share patterns of signaling and gene expression in early development. Thus far, the extent to which genital and limb transcriptional networks also share cis-regulatory elements has remained unexplored. We show that many limb enhancers are retained in snake genomes, suggesting that these elements may function in non-limb tissues. Consistent with this, our analysis of cis-regulatory activity in mice and Anolis lizards reveals that patterns of enhancer activity in embryonic limbs and genitalia overlap heavily. In mice, deletion of HLEB, an enhancer of Tbx4, produces defects in hindlimbs and genitalia, establishing the importance of this limb-genital enhancer for development of these different appendages. Further analyses demonstrate that the HLEB of snakes has lost hindlimb enhancer function while retaining genital activity. Our findings identify roles for Tbx4 in genital development and highlight deep similarities in cis-regulatory activity between limbs and genitalia.