Loss and Re-emergence of Legs in Snakes by Modular Evolution of Sonic hedgehog and HOXD Enhancers

Ancient developmental genes underlie evolutionary novelties in walking fish

A critical question in biology is how new traits evolve, but studying this in wild animals remains challenging. Here, we probe the genetic basis of trait gain in sea robin fish, which have evolved specialized leg-like appendages for locomotion and digging along the ocean floor. We use genome sequencing, transcriptional profiling, and interspecific hybrid analysis to explore the molecular and developmental basis of leg formation. We identified the ancient, conserved transcription factor tbx3a as a major determinant of sensory leg development. Genome editing confirms that tbx3a is required for normal leg formation in sea robins, and for formation of enlarged central nervous system lobes, sensory papillae, and adult digging behavior. Our study establishes sea robins as a model organism for studying the evolution of major trait gain and illustrates how ancient developmental control genes can underlie novel organ formation.

Tgfbr1 controls developmental plasticity between the hindlimb and external genitalia by remodeling their regulatory landscape

The hindlimb and external genitalia of present-day tetrapods are thought to derive from an ancestral common primordium that evolved to generate a wide diversity of structures adapted for efficient locomotion and mating in the ecological niche occupied by the species. We show that despite long evolutionary distance from the ancestral condition, the early primordium of the mouse external genitalia preserved the capacity to take hindlimb fates. In the absence of Tgfbr1, the pericloacal mesoderm generates an extra pair of hindlimbs at the expense of the external genitalia. It has been shown that the hindlimb and the genital primordia share many of their key regulatory factors. Tgfbr1 controls the response to those factors by modulating the accessibility status of regulatory elements that control the gene regulatory networks leading to the formation of genital or hindlimb structures. Our work uncovers a remarkable tissue plasticity with potential implications in the evolution of the hindlimb/genital area of tetrapods, and identifies an additional mechanism for Tgfbr1 activity that might also contribute to the control of other physiological or pathological processes.

Resolving the mechanisms underlying epithelial-to-mesenchymal transition of the lateral plate mesoderm

Formation of the vertebrate limb buds begins with a localized epithelial-to-mesenchymal transition (EMT) of the somatic lateral plate mesoderm (LPM). While the processes that drive proliferation and outgrowth of the limb mesenchyme are well established, the fundamental mechanisms that precede this process and initiate EMT are less understood. In this review, we outline putative drivers of EMT of the LPM, drawing from analyses across a range of vertebrates and developmental models. We detail the expression patterns of key EMT transcriptional regulators in the somatic LPM of the presumptive limb fields, and their potential role in producing a mesenchymal cell fate. These include a putative cooperative role between the EMT inducers PRRX1 and TWIST1, supported by evidence in zebrafish and chicken models but unconfirmed data from mice. As such, additional functional data are required to definitively determine the mechanisms that initiate and drive EMT of the somatic LPM, a critical transition preceding formation of the limb bud mesenchyme.

Identification and characterization of Hand2 upstream genomic enhancers active in developing stomach and limbs

BACKGROUND

The bHLH transcription factor HAND2 plays important roles in the development of the embryonic heart, face, limbs, and sympathetic and enteric nervous systems. To define how and when HAND2 regulates these developmental systems, requires understanding the transcriptional regulation of Hand2.

RESULTS

Remarkably, Hand2 is flanked by an extensive upstream gene desert containing a potentially diverse enhancer landscape. Here, we screened the regulatory interval 200 kb proximal to Hand2 for putative enhancers using evolutionary conservation and histone marks in Hand2-expressing tissues. H3K27ac signatures across embryonic tissues pointed to only two putative enhancer regions showing deep sequence conservation. Assessment of the transcriptional enhancer potential of these elements using transgenic reporter lines uncovered distinct in vivo enhancer activities in embryonic stomach and limb mesenchyme, respectively. Activity of the identified stomach enhancer was restricted to the developing antrum and showed expression within the smooth muscle and enteric neurons. Surprisingly, the activity pattern of the limb enhancer did not overlap Hand2 mRNA but consistently yielded a defined subectodermal anterior expression pattern within multiple transgenic lines.

CONCLUSIONS

Together, these results start to uncover the diverse regulatory potential inherent to the Hand2 upstream regulatory interval.

Trema orientale (L.) Blume: A review of its taxonomy, traditional uses, phytochemistry, pharmacological activities and domestication potential

Trema orientale (L.) Blume is an important medicinal plant with multiple applications for treating several disease conditions. This study compiled published data on botanical, traditional uses, phytochemistry, pharmacology, and antimicrobials, coupled with discussing the conservation and domestication potential of T. orientale. Data were sourced from databases such as Google Scholar, PubMed, Scopus, Elsevier Plants of the World Online (Kew Science), Global Biodiversity Information Facility (GBIF), and World Flora Online (WFO), using key search terms: Trema orientale or orientalis, phytochemistry pharmacology, taxonomy, and domestication with Boolean operators to include and exclude articles for the review. The review indicated that molecular studies have shown that T. orientale is closely related to a sister group of Cannabis through plastome phylogenetic evidence which accounts for its transfer from Ulmaceae to the Cannabaceae family. T. orientale is distributed across several African countries and has recently been assessed as the Least Concern by the IUCN Red List of Threatened Species. Nevertheless, deforestation continues to pose an extinction risk to their population. Currently, 31 compounds have been isolated from different parts of T. orientale justifying many traditional uses accredited to it. T. orientale is considered a dose-dependent safe remedy for the treatment of infectious diseases, cancer, cardiovascular diseases and other disease conditions ascribed to it except for its continuous application. This review underscores the domestication potential of T. orientale including evidence of molecular markers, soil seed banks, and promising outcomes of germination experimentations. This, therefore, presents significant gains toward sustainable utilization of Trema orientale.

Embryonic development of the pelvic girdle and hindlimb skeletal elements in Anilius scytale (Linnaeus, 1758) (Serpentes: Aniliidae)

Anilius scytale is the sister lineage of all other alethinophidian snakes. Morphology of the hind limb complex in adult A. scytale (Aniliidae) has been documented. We herein, for the first time, describe the embryology of the skeletal elements of its hind limb and pelvic girdle and contextualize the evolution of these structures. We identified pregnant females of A. scytale in the Herpetology Collection of the Museu Paraense Emílio Goeldi and separated 40 embryos. The embryos were sequentially staged using external and internal anatomy, collectively comprising a developmental series representing six stages. We cleared-stained one specimen of stages 31, 34, 36, and 37. Using the embryological information gleaned from A. scytale, we reinterpret evidence relating to the ossification of the pelvis and hindlimbs. In A. scytale hindlimb buds develop as transient structures that developed before Stage 30 and regresses in subsequent stages. There is no external or internal evidence of the forelimb or scapular girdle. From Stage 31 onwards the ischium, pubis, ilium, femur and zeugopodial cartilages are visible. Pubis and femur ossify towards the end of embryonic life, and cloacal spurs do not develop in the embryo. Skeletal elements of the hindlimb and pelvic girdle develop initially in the ventral zone of the cloaca-tail region. In subsequent stages the hindlimb and pelvic girdle elements migrate dorsally, with the pubis/ischium positioned medial to the ribs. A similar process may be associated with the achievement of the condition of the pelvic girdle in adults of scolecophidians, pythonids and boids.

Reversibility of digit loss revisited: Limb diversification in Bachia lizards (gymnophthalmidae)

Strict interpretations of the Dollo's Law lead to postulation that trait loss is irreversible and organisms never recover ancestral phenotypes. Dollo, however, admitted the possibility of reversals in trait loss when predicted differences between reversed (derived) and ancestral forms. Phenotypic signatures from reversals are expected, as the historical context of a reversal in trait loss differs from the initial setting where the trait originally evolved. This article combines morphological and molecular information for Bachia scolecoides to discuss phenotypic and genetic patterns established during processes that reversed digit loss in Gymnophthalmidae (also termed microteiid lizards). Results suggest that pathways leading to the derived tetradactyl state of B. scolecoides comprise particularities in their origin and associated processes. Autopodial bones of B. scolecoides lack digit identity, and muscle anatomy is very similar between manus and pes. Gymnophthalmidae sequence patterns in the limb-specific sonic hedgehog enhancer (ZRS) suggest that regulation of shh expression did not degenerate in Bachia, given the prediction of similar motifs despite mutations specific to Bachia. Persistence of developmental mechanisms might explain intermittent character expression leading to reversals of digit loss, as ZRS signaling pathways remain active during the development of at least one pair of appendices in Bachia, especially if some precursors persisted at early stages. Patterns of ZRS sequences suggest that irreversibility of trait loss might be lineage-specific (restricted to Gymnophthalmini) and contingent to the type of signature established. These results provide insights regarding possible mechanisms that may allow reactivation of developmental programs in specific regions of the embryo.

The genetic basis of novel trait gain in walking fish

A major goal in biology is to understand how organisms evolve novel traits. Multiple studies have identified genes contributing to regressive evolution, the loss of structures that existed in a recent ancestor. However, fewer examples exist for genes underlying constructive evolution, the gain of novel structures and capabilities in lineages that previously lacked them. Sea robins are fish that have evolved enlarged pectoral fins, six mobile locomotory fin rays (legs) and six novel macroscopic lobes in the central nervous system (CNS) that innervate the corresponding legs. Here, we establish successful husbandry and use a combination of transcriptomics, CRISPR-Cas9 editing, and behavioral assays to identify key transcription factors that are required for leg formation and function in sea robins. We also generate hybrids between two sea robin species with distinct leg morphologies and use allele-specific expression analysis and gene editing to explore the genetic basis of species-specific trait diversity, including a novel sensory gain of function. Collectively, our study establishes sea robins as a new model for studying the genetic basis of novel organ formation, and demonstrates a crucial role for the conserved limb gene tbx3a in the evolution of chemosensory legs in walking fish.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Increased chromatin accessibility promotes the evolution of a transcriptional silencer in Drosophila

The loss of discrete morphological traits, the most common evolutionary transition, is typically driven by changes in developmental gene expression. Mutations accumulating in regulatory elements of these genes can disrupt DNA binding sites for transcription factors patterning their spatial expression, or delete entire enhancers. Regulatory elements, however, may be silenced through changes in chromatin accessibility or the emergence of repressive elements. Here, we show that increased chromatin accessibility at the gene yellow, combined with the gain of a repressor site, underlies the loss of a wing spot pigmentation pattern in a Drosophila species. The gain of accessibility of this repressive element is regulated by E93, a transcription factor governing the progress of metamorphosis. This convoluted evolutionary scenario contrasts with the parsimonious mutational paths generally envisioned and often documented for morphological losses. It illustrates how evolutionary changes in chromatin accessibility may directly contribute to morphological diversification.

The role of timing in the development and evolution of the limb

The term heterochrony was coined to describe changes in the timing of developmental processes relative to an ancestral state. Limb development is a well-suited system to address the contribution of heterochrony to morphological evolution. We illustrate how timing mechanisms have been used to establish the correct pattern of the limb and provide cases where natural variations in timing have led to changes in limb morphology.

Evolutionary genetics of flipper forelimb and hindlimb loss from limb development-related genes in cetaceans

Background

Cetacean hindlimbs were lost and their forelimb changed into flippers characterized by webbed digits and hyperphalangy, thus allowing them to adapt to a completely aquatic environment. However, the underlying molecular mechanism behind cetacean limb development remains poorly understood.

Results

In the present study, we explored the evolution of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with that of other mammals. TBX5, a forelimb specific expression gene, was identified to have been under accelerated evolution in the ancestral branches of cetaceans. In addition, 32 cetacean-specific changes were examined in the SHH signaling network (SHH, PTCH1, TBX5, BMPs and SMO), within which mutations could yield webbed digits or an additional phalange. These findings thus suggest that the SHH signaling network regulates cetacean flipper formation. By contrast, the regulatory activity of the SHH gene enhancer—ZRS in cetaceans—was significantly lower than in mice, which is consistent with the cessation of SHH gene expression in the hindlimb bud during cetacean embryonic development. It was suggested that the decreased SHH activity regulated by enhancer ZRS might be one of the reasons for hindlimb degeneration in cetaceans. Interestingly, a parallel / convergent site (D42G) and a rapidly evolving CNE were identified in marine mammals in FGF10 and GREM1, respectively, and shown to be essential to restrict limb bud size; this is molecular evidence explaining the convergence of flipper-forelimb and shortening or degeneration of hindlimbs in marine mammals.

Conclusions

We did evolutionary analyses of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with those of other mammals to provide novel insights into the molecular basis of flipper forelimb and hindlimb loss in cetaceans.

Gene expression changes during the evolution of the tetrapod limb

Major changes in the vertebrate anatomy have preceded the conquest of land by the members of this taxon, and continuous changes in limb shape and use have occurred during the later radiation of tetrapods. While the main, conserved mechanisms of limb development have been discerned over the past century using a combination of classical embryological and molecular methods, only recent advances made it possible to identify and study the regulatory changes that have contributed to the evolution of the tetrapod appendage. These advances include the expansion of the model repertoire from traditional genetic model species to non-conventional ones, a proliferation of predictive mathematical models that describe gene interactions, an explosion in genomic data and the development of high-throughput methodologies. These revolutionary innovations make it possible to identify specific mutations that are behind specific transitions in limb evolution. Also, as we continue to apply them to more and more extant species, we can expect to gain a fine-grained view of this evolutionary transition that has been so consequential for our species as well.

How vertebrates got their bite

A newly discovered enhancer region may have allowed vertebrates to evolve the ability to open and close their jaws.

Developmental basis of evolutionary lung loss in plethodontid salamanders

One or more members of four living amphibian clades have independently dispensed with pulmonary respiration and lack lungs, but little is known of the developmental basis of lung loss in any taxon. We use morphological, molecular, and experimental approaches to examine the Plethodontidae, a dominant family of salamanders, all of which are lungless as adults. We confirm an early anecdotal report that plethodontids complete early stages of lung morphogenesis: Transient embryonic lung primordia form but regress by apoptosis before hatching. Initiation of pulmonary development coincides with expression of the lung-specification gene Wnt2b in adjacent mesoderm, and the lung rudiment expresses pulmonary markers Nkx2.1 and Sox9. Lung developmental-genetic pathways are at least partially conserved despite the absence of functional adult lungs for at least 25 and possibly exceeding 60 million years. Adult lung loss appears associated with altered expression of signaling molecules that mediate later stages of tracheal and pulmonary development.

The pZRS non-coding regulatory mutation resulting in triphalangeal thumb-polysyndactyly syndrome changes the pattern of local interactions

Herein, we report on a large Polish family presenting with a classical triphalangeal thumb-polysyndactyly syndrome (TPT-PS). This rare congenital limb anomaly is generally caused by microduplications encompassing the Sonic Hedgehog (SHH) limb enhancer, termed the zone of polarizing activity (ZPA) regulatory sequence (ZRS). Recently, a pathogenic variant in the pre-ZRS (pZRS), a conserved sequence located near the ZRS, has been described in a TPT-PS Dutch family. We performed targeted ZRS sequencing, array comparative genomic hybridization, and whole-exome sequencing. Next, we sequenced the recently described pZRS region. Finally, we performed a circular chromatin conformation capture-sequencing (4C-seq) assay on skin fibroblasts of one affected family member and control samples to examine potential alterations in the SHH regulatory domain and functionally characterize the identified variant. We found that all affected individuals shared a recently identified pathogenic point mutation in the pZRS region: NC_000007.14:g.156792782C>G (GRCh38/hg38), which is the same as in the Dutch family. The results of 4C-seq experiments revealed increased interactions within the whole SHH regulatory domain (SHH-LMBR1 TAD) in the patient compared to controls. Our study expands the number of TPT-PS families carrying a pathogenic alteration of the pZRS and underlines the importance of routine pZRS sequencing in the genetic diagnostics of patients with TPT-PS or similar phenotypes. The pathogenic mutation causative for TPT-PS in our patient gave rise to increased interactions within the SHH regulatory domain in yet unknown mechanism.

Macroevolutionary Analyses Provide New Evidence of Phasmid Wings Evolution as a Reversible Process

The concept that complex ancestral traits can never be recovered after their loss is still widely accepted, despite phylogenetic and molecular approaches suggest instances where phenotypes may have been lost throughout the evolutionary history of a clade and subsequently reverted back in derived lineages. One of the first and most notable examples of such a process is wing evolution in phasmids; this polyneopteran order of insects, which comprises stick and leaf insects, has played a central role in initiating a long-standing debate on the topic. In this study, a novel and comprehensive time tree including over 300 Phasmatodea species is used as a framework for investigating wing evolutionary patterns in the clade. Despite accounting for several possible biases and sources of uncertainty, macroevolutionary analyses consistently revealed multiple reversals to winged states taking place after their loss, and reversibility is coupled with higher species diversification rates. Our findings support a loss of or reduction in wings that occurred in the lineage leading to the extant phasmid most recent common ancestor, and brachyptery is inferred to be an unstable state unless co-opted for nonaerodynamic adaptations. We also explored how different assumptions of wing reversals probability could impact their inference: we found that until reversals are assumed to be over 30 times more unlikely than losses, they are consistently inferred despite uncertainty in tree and model parameters. Our findings demonstrate that wing evolution is a reversible and dynamic process in phasmids and contribute to our understanding of complex trait evolution.

Snake-like limb loss in a Carboniferous amniote

Among living tetrapods, many lineages have converged on a snake-like body plan, where extreme axial elongation is accompanied by reduction or loss of paired limbs. However, when and how this adaptive body plan first evolved in amniotes remains poorly understood. Here, we provide insights into this question by reporting on a new taxon of molgophid recumbirostran, Nagini mazonense gen. et sp. nov., from the Francis Creek Shale (309-307 million years ago) of Illinois, United States, that exhibits extreme axial elongation and corresponding limb reduction. The molgophid lacks entirely the forelimb and pectoral girdle, thus representing the earliest occurrence of complete loss of a limb in a taxon recovered phylogenetically within amniotes. This forelimb-first limb reduction is consistent with the pattern of limb reduction that is seen in modern snakes and contrasts with the hindlimb-first reduction process found in many other tetrapod groups. Our findings suggest that a snake-like limb-reduction mechanism may be operating more broadly across the amniote tree.

Sonic hedgehog is Essential for Proximal-Distal Outgrowth of the Limb Bud in Salamanders

The developing forelimb has been a foundational model to understand how specified progenitor cells integrate genetic information to produce the tetrapod limb bauplan. Although the reigning hypothesis is that all tetrapods develop limbs in a similar manner, recent work suggests that urodeles have evolved a derived mode of limb dvelopment. Here, we demonstrate through pharmacological and genetic inactivation of Sonic hedgehog (Shh) signaling in axolotls that Shh directs expansion and survival of limb progenitor cells in addition to patterning the limb across the proximodistal and antero-posterior axis. In contrast to inactivation of Shh in mouse or chick embryos where a humerus, radius, and single digit develop, Shh crispant axolotls completely lack forelimbs. In rescuing limb development by implanting SHH-N protein beads into the nascent limb field of Shh crispants, we show that the limb field is specified in the absence of Shh and that hedgehog pathway activation is required to initiate proximodistal outgrowth. When our results are examined alongside other derived aspects of salamander limb development and placed in a phylogenetic context, a new hypothesis emerges whereby the ability for cells at an amputation plane to activate morphogenesis and regenerate a limb may have evolved uniquely in urodeles.

Long-read genome sequencing provides molecular insights into scavenging and societal complexity in spotted hyena Crocuta crocuta

The spotted hyena (Crocuta crocuta) is a large and unique terrestrial carnivore. It is a particularly fascinating species due to its distinct phenotypic traits, especially its complex social structure and scavenging lifestyle, with associated high dietary exposure to microbial pathogens. However, the underlying molecular mechanisms related to these phenotypes remain elusive. Here, we sequenced and assembled a high-quality long-read genome of the spotted hyena, with a contig N50 length of ∼13.75 Mb. Based on comparative genomics, immunoglobulin family members (e.g., IGKV4-1) showed significant adaptive duplications in the spotted hyena and striped hyena. Furthermore, immune-related genes (e.g., CD8A, LAG3, and TLR3) experienced species-specific positive selection in the spotted hyena lineage. These results suggest that immune tolerance between the spotted hyena and closely related striped hyena has undergone adaptive divergence to cope with prolonged dietary exposure to microbial pathogens from scavenging. Furthermore, we provided the potential genetic insights underlying social complexity, hinting at social behavior and cognition. Specifically, the RECNE-associated genes (e.g., UGP2 and ACTR2) in the spotted hyena genome are involved in regulation of social communication. Taken together, our genomic analyses provide molecular insights into the scavenging lifestyle and societal complexity of spotted hyenas.

Theories, laws, and models in evo-devo

Evolutionary developmental biology (evo-devo) is the study of the evolution of developmental mechanisms. Here, I review some of the theories, models, and laws in evo-devo, past and present. Nineteenth-century evo-devo was dominated by recapitulation theory and archetypes. It also gave us germ layer theory, the vertebral theory of the skull, floral organs as modified leaves, and the "inverted invertebrate" theory, among others. Newer theories and models include the frameshift theory, the genetic toolkit for development, the ABC model of flower development, the developmental hourglass, the zootype, Urbilateria, and the hox code. Some of these new theories show the influence of archetypes and recapitulation. Interestingly, recent studies support the old "primordial leaf," "inverted invertebrate," and "segmented head" theories. Furthermore, von Baer's first three laws may now need to be rehabilitated, and the hourglass model modified, in view of what Abzhanov has pointed out about the maternal-zygotic transition. There are many supposed "laws" of evo-devo but I argue that these are merely generalizations about trends in particular lineages. I argue that the "body plan" is an archetype, and is often used in such a way that it lacks any scientific meaning. Looking to the future, one challenge for evo-devo will be to develop new theories and models to accommodate the wealth of new data from high-throughput sequencing, including single-cell sequencing. One step in this direction is the use of sophisticated in silico analyses, as in the "transcriptomic hourglass" models.

Comparative development of limb musculature in phylogenetically and ecologically divergent lizards

BACKGROUND

Squamate reptiles (lizards, snakes, and amphisbaenians) exhibit incredible diversity in their locomotion, behavior, morphology, and ecological breadth. Although they often are used as models of locomotor diversity, surprisingly little attention has been given to muscle development in squamate reptiles. In fact, the most detailed examination was conducted almost 80 years ago and solely focused on the proximal limb regions. Herein, we present forelimb and hindlimb muscle morphogenesis data for three lizard species with different locomotion and feeding strategies: the desert grassland whiptail lizard, the central bearded dragon, and the veiled chameleon. This study fills critical gaps in our understanding of muscle morphogenesis in squamate reptiles and presents a comparative and temporospatial analysis of muscle development.

RESULTS

Our results reveal a conserved pattern of early muscle development among lizards with different adult morphologies and ecologies. The variations that exist are concentrated in distal regions, particularly the specialized autopodia of chameleons, where differentiation of muscles associated with the digits is delayed.

CONCLUSIONS

The chameleon autopod provides an example of major evolutionary modifications to the skeleton with only minor disruption of the conserved order and pattern of limb muscle development. This robustness of muscle patterning facilitates the evolution of extreme yet functional phenotypes.

(Non)Parallel developmental mechanisms in vertebrate appendage reduction and loss

Appendages have been reduced or lost hundreds of times during vertebrate evolution. This phenotypic convergence may be underlain by shared or different molecular mechanisms in distantly related vertebrate clades. To investigate, we reviewed the developmental and evolutionary literature of appendage reduction and loss in more than a dozen vertebrate genera from fish to mammals. We found that appendage reduction and loss was nearly always driven by modified gene expression as opposed to changes in coding sequences. Moreover, expression of the same genes was repeatedly modified across vertebrate taxa. However, the specific mechanisms by which expression was modified were rarely shared. The multiple routes to appendage reduction and loss suggest that adaptive loss of function phenotypes might arise routinely through changes in expression of key developmental genes.

Spatial regulation by multiple Gremlin1 enhancers provides digit development with cis-regulatory robustness and evolutionary plasticity

Precise cis-regulatory control of gene expression is essential for normal embryogenesis and tissue development. The BMP antagonist Gremlin1 (Grem1) is a key node in the signalling system that coordinately controls limb bud development. Here, we use mouse reverse genetics to identify the enhancers in the Grem1 genomic landscape and the underlying cis-regulatory logics that orchestrate the spatio-temporal Grem1 expression dynamics during limb bud development. We establish that transcript levels are controlled in an additive manner while spatial regulation requires synergistic interactions among multiple enhancers. Disrupting these interactions shows that altered spatial regulation rather than reduced Grem1 transcript levels prefigures digit fusions and loss. Two of the enhancers are evolutionary ancient and highly conserved from basal fishes to mammals. Analysing these enhancers from different species reveal the substantial spatial plasticity in Grem1 regulation in tetrapods and basal fishes, which provides insights into the fin-to-limb transition and evolutionary diversification of pentadactyl limbs.

The BMP antagonist Gremlin1 balances BMP and SHH signalling, endowing limb bud development with robustness. Here, the authors identify enhancers controlling Grem1 levels in an additive, and spatial regulation in a synergistic manner, providing digit patterning with cis-regulatory robustness and evolutionary plasticity.

Picking up the threads: Comparative osteology and associated cartilaginous elements for members of the genus Trilepida Hedges, 2011 (Serpentes, Leptotyphlopidae) with new insights on the Epictinae systematics

The threadsnakes of the family Leptotyphlopidae have been historically neglected in terms of their natural history, ecology, systematics, and morphology. Given the relevance of morphological data for resolving systematic, evolutionary, and functional issues, we aimed to provide a detailed comparative description of osteology and associated cartilaginous elements for members of the genus Trilepida. Data were obtained through high-resolution computed tomography images, cleared and stained specimens, and radiography images of a total of 47 specimens and 12 species. Both cranial and axial osteology characters exhibited a relevant degree of intraspecific variation regarding qualitative and quantitative data associated with skull and vertebrae foramina and the shape of bony sutures and processes. The high representativeness of examined species and specimens allows us to provide a comprehensive discussion on the inter- and intraspecific osteological variation, as well as a compelling osteological diagnosis for the entire genus. Trilepida spp. differ from all Epictinae by the presence of the following combination of characters: paired nasals, fused supraoccipitals (distinct from parietal, prootics, and otooccipitals), a single (fused) parietal without a dorsal fontanelle, and the basioccipital participating in the foramen magnum (except in Trilepida nicefori). Our results reinforce the need for integration of detailed anatomical traits to usually conserved external morphological data to provide accurate diagnostic features for Epictinae. In addition, new phylogenetic hypotheses or even taxonomic re-allocations may broadly benefit from these detailed comparative studies.

Selection on phalanx development in the evolution of the bird wing

The frameshift hypothesis is a widely-accepted model of bird wing evolution. This hypothesis postulates a shift in positional values, or molecular-developmental identity, that caused a change in digit phenotype. The hypothesis synthesised developmental and palaeontological data on wing digit homology. The 'most anterior digit' (MAD) hypothesis presents an alternative view based on changes in transcriptional regulation in the limb. The molecular evidence for both hypotheses is that the most anterior digit expresses Hoxd13 but not Hoxd11 and Hoxd12. This digit I 'signature' is thought to characterise all amniotes. Here, we studied Hoxd expression patterns in a phylogenetic sample of 18 amniotes. Instead of a conserved molecular signature in digit I, we find wide variation of Hoxd11, Hoxd12 and Hoxd13 expression in digit I. Patterns of apoptosis, and Sox9 expression, a marker of the phalanx-forming region, suggest that phalanges were lost from wing digit IV because of early arrest of the phalanx-forming region followed by cell death. Finally, we show that multiple amniote lineages lost phalanges with no frameshift. Our findings suggest that the bird wing evolved by targeted loss of phalanges under selection. Consistent with our view, some recent phylogenies based on dinosaur fossils eliminate the need to postulate a frameshift in the first place. We suggest that the phenotype of the Archaeopteryx lithographica wing is also consistent with phalanx loss. More broadly, our results support a gradualist model of evolution based on tinkering with developmental gene expression.

Over-expression of the bottlenose dolphin Hoxd13 gene in zebrafish provides new insights into the cetacean flipper formation

Cetaceans have evolved elongated soft-tissue flipper with digits made of hyperphalangy. Cetaceans were found to have 2-3 more alanine residues in Hoxd13 than other mammals, which were suggested to be related to their flipper. However, how Hoxd13 regulates other genes and induces hyperphalangy in cetaceans remain poorly understood. Here, we overexpressed the bottlenose dolphin Hoxd13 in zebrafish (Danio rerio). Combined with transcriptome data and evolutionary analyses, our results revealed that the Wingless/Integrated (Wnt) and Hedgehog signaling pathways and multiple genes might regulate hyperphalangy development in cetaceans. Meanwhile, the Notch and mitogen-activated protein kinase (Mapk) signaling pathways and Fibroblast growth factor receptor 1 (Fgfr1) are probably correlated with interdigital tissues retained in the cetacean flipper. In conclusion, this is the first study to use a transgenic zebrafish to explore the molecular evolution of Hoxd13 in cetaceans, and it provides new insights into cetacean flipper formation.

Amplicon genome fishing (AGF): a rapid and efficient method for sequencing target cis-regulatory regions in nonmodel organisms

Cis-regulatory sequences play a crucial role in regulating gene expression and are evolutionary hot spots that drive phenotypic divergence among organisms. Sequencing some cis-regulatory regions of interest in many different species is common in comparative genetic studies. For nonmodel organisms lacking genomic data, genome walking is often the preferred method for this type of application. However, applying genome walking will be laborious and time-consuming when the number of cis-regulatory regions and species to be analyzed is large. In this study, we propose a novel method called amplicon genome fishing (AGF), which can isolate and sequence cis-regulatory regions of interest for any organism. The main idea of the AGF method is to use fragments amplified from the target cis-regulatory regions as enrichment baits to capture and sequence the whole target cis-regulatory regions from genomic library pools. Unlike genome walking, the AGF method is based on hybridization capture and high-throughput sequencing, which makes this method rapid and efficient for projects where some cis-regulatory regions have to be sequenced for many species. We used human amplicons as capture baits and successfully sequenced five target enhancer regions of Homo sapiens, Mus musculus, Gallus gallus, and Xenopus tropicalis, proving the feasibility and repeatability of AGF. To show the utility of the AGF method in real studies, we used it to sequence the ZRS enhancer, a cis-regulatory region associated with the limb loss of snakes, for twenty-three vertebrate species (includes many limbless species never sequenced before). The newly obtained ZRS sequences provide new perspectives into the relationship between the ZRS enhancer's evolution and limb loss in major tetrapod lineages.

On the sequence heterochrony of cranial ossification of bats in light of Haeckel's recapitulation theory

Haeckel's recapitulation theory has been a controversial topic in evolutionary biology. However, we have seen some recent cases applying Haeckel's view to interpret the interspecific variation of prenatal ontogeny. To revisit the validity of Haeckel's recapitulation theory, we take bats that have undergone drastic morphological changes and possess a characteristic ecology as a case study. All members of Rhinolophoidea and Yangochiroptera can generate an ultrasonic pulse from the larynx to interpret surrounding objects (laryngeal echolocation) whereas Pteropodidae lacks such ability. It is known that the petrosal bone is particularly derived in shape and expanded in laryngeal echolocators. If Haeckel's recapitulation theory holds, the formation of this derived trait should occur later than those of other bones. Therefore, we compared the prenatal ossification timing of the petrosal in 15 bat species and five outgroup species. We found that the ossification of the petrosal is accelerated in laryngeal echolocators while it is the last bone to ossify in non-laryngeal echolocating bats and non-volant mammals, which runs counter to the prediction generated by Haeckel's recapitulation theory. We point out the evolutionarily labile nature of trait developmental timing and emphasize that Haeckel's recapitulation theory does not hold in many cases. We caution that generating predictions on ancestral conditions and evolutionary history leading from Haeckel's recapitulation theory is not well supported.

Embryonic evidence uncovers convergent origins of laryngeal echolocation in bats

Bats are the second-most speciose group of mammals, comprising 20% of species diversity today. Their global explosion, representing one of the greatest adaptive radiations in mammalian history, is largely attributed to their ability of laryngeal echolocation and powered flight, which enabled them to conquer the night sky, a vast and hitherto unoccupied ecological niche. While there is consensus that powered flight evolved only once in the lineage, whether laryngeal echolocation has a single origin in bats or evolved multiple times independently remains disputed. Here, we present developmental evidence in support of laryngeal echolocation having multiple origins in bats. This is consistent with a non-echolocating bat ancestor and independent gain of echolocation in Yinpterochiroptera and Yangochiroptera, as well as the gain of primitive echolocation in the bat ancestor, followed by convergent evolution of laryngeal echolocation in Yinpterochiroptera and Yangochiroptera, with loss of primitive echolocation in pteropodids. Our comparative embryological investigations found that there is no developmental difference in the hearing apparatus between non-laryngeal echolocating bats (pteropodids) and terrestrial non-bat mammals. In contrast, the echolocation system is developed heterotopically and heterochronically in the two phylogenetically distant laryngeal echolocating bats (rhinolophoids and yangochiropterans), providing the first embryological evidence that the echolocation system evolved independently in these bats.

Diversification of the vertebrate limb: sequencing the events

Naturalists leading up to the early 20th century were captivated by the diversity of limb form and function and described its development in a variety of species. The advent of discoveries in genetics followed by molecular biology led to focused efforts in few 'model' species, namely mouse and chicken, to understand conserved mechanisms of limb axis specification and development of the musculoskeletal system. 'Non-traditional' species largely fell by the wayside until their recent resurgence into the spotlight with advances in next-generation sequencing technologies (NGS). In this review, we focus on how the use of NGS has provided insights into the development, loss, and diversification of amniote limbs. Coupled with advances in chromatin interrogation techniques and functional tests in vivo, NGS is opening possibilities to understand the genetic mechanisms that govern the remarkable radiation of vertebrate limb form and function.

Congenital Malformations in Sea Turtles: Puzzling Interplay between Genes and Environment

Simple Summary

Congenital malformations can lead to embryonic mortality in many species, and sea turtles are no exception. Genetic and/or environmental alterations occur during early development in the embryo, and may produce aberrant phenotypes, many of which are incompatible with life. Causes of malformations are multifactorial; genetic factors may include mutations, chromosomal aberrations, and inbreeding effects, whereas non-genetic factors may include nutrition, hyperthermia, low moisture, radiation, and contamination. It is possible to monitor and control some of these factors (such as temperature and humidity) in nesting beaches, and toxic compounds in feeding areas, which can be transferred to the embryo through their lipophilic properties. In this review, we describe possible causes of different types of malformations observed in sea turtle embryos, as well as some actions that may help reduce embryonic mortality.

Abstract

The completion of embryonic development depends, in part, on the interplay between genetic factors and environmental conditions, and any alteration during development may affect embryonic genetic and epigenetic regulatory pathways leading to congenital malformations, which are mostly incompatible with life. Oviparous reptiles, such as sea turtles, that produce numerous eggs in a clutch that is buried on the beach provide an opportunity to study embryonic mortality associated with malformations that occur at different times during development, or that prevent the hatchling from emerging from the nest. In sea turtles, the presence of congenital malformations frequently leads to mortality. A few years ago, a detailed study was performed on external congenital malformations in three species of sea turtles from the Mexican Pacific and Caribbean coasts, the hawksbill turtle, Eretmochelys imbricata (n = 23,559 eggs), the green turtle, Chelonia mydas (n = 17,690 eggs), and the olive ridley, Lepidochelys olivacea (n = 20,257 eggs), finding 63 types of congenital malformations, of which 38 were new reports. Of the three species, the olive ridley showed a higher incidence of severe anomalies in the craniofacial region (49%), indicating alterations of early developmental pathways; however, several malformations were also observed in the body, including defects in the carapace (45%) and limbs (33%), as well as pigmentation disorders (20%), indicating that deviations occurred during the middle and later stages of development. Although intrinsic factors (i.e., genetic mutations or epigenetic modifications) are difficult to monitor in the field, some environmental factors (such as the incubation temperature, humidity, and probably the status of feeding areas) are, to some extent, less difficult to monitor and/or control. In this review, we describe the aetiology of different malformations observed in sea turtle embryos, and provide some actions that can reduce embryonic mortality.

A farewell to arms and legs: a review of limb reduction in squamates

Elongated snake-like bodies associated with limb reduction have evolved multiple times throughout vertebrate history. Limb-reduced squamates (lizards and snakes) account for the vast majority of these morphological transformations, and thus have great potential for revealing macroevolutionary transitions and modes of body-shape transformation. Here we present a comprehensive review on limb reduction, in which we examine and discuss research on these dramatic morphological transitions. Historically, there have been several approaches to the study of squamate limb reduction: (i) definitions of general anatomical principles of snake-like body shapes, expressed as varying relationships between body parts and morphometric measurements; (ii) framing of limb reduction from an evolutionary perspective using morphological comparisons; (iii) defining developmental mechanisms involved in the ontogeny of limb-reduced forms, and their genetic basis; (iv) reconstructions of the evolutionary history of limb-reduced lineages using phylogenetic comparative methods; (v) studies of functional and biomechanical aspects of limb-reduced body shapes; and (vi) studies of ecological and biogeographical correlates of limb reduction. For each of these approaches, we highlight their importance in advancing our understanding, as well as their weaknesses and limitations. Lastly, we provide suggestions to stimulate further studies, in which we underscore the necessity of widening the scope of analyses, and of bringing together different perspectives in order to understand better these morphological transitions and their evolution. In particular, we emphasise the importance of investigating and comparing the internal morphology of limb-reduced lizards in contrast to external morphology, which will be the first step in gaining a deeper insight into body-shape variation.

Regulation of vertebrate forelimb development and wing reduction in the flightless emu

The vertebrate limb is a dynamic structure which has evolved into many diverse forms to facilitate complex behavioral adaptations. The principle molecular and cellular processes that underlie development of the vertebrate limb are well characterized. However, how these processes are altered to drive differential limb development between vertebrates is less well understood. Several vertebrate models are being utilized to determine the developmental basis of differential limb morphogenesis, though these typically focus on later patterning of the established limb bud and may not represent the complete developmental trajectory. Particularly, heterochronic limb development can occur prior to limb outgrowth and patterning but receives little attention. This review summarizes the genetic regulation of vertebrate forelimb diversity, with particular focus on wing reduction in the flightless emu as a model for examining limb heterochrony. These studies highlight that wing reduction is complex, with heterochronic cellular and genetic events influencing the major stages of limb development. Together, these studies provide a broader picture of how different limb morphologies may be established during development.

Modes of genetic adaptations underlying functional innovations in the rumen

The rumen is the hallmark organ of ruminants and hosts a diverse ecosystem of microorganisms that facilitates efficient digestion of plant fibers. We analyzed 897 transcriptomes from three Cetartiodactyla lineages: ruminants, camels and cetaceans, as well as data from ruminant comparative genomics and functional assays to explore the genetic basis of rumen functional innovations. We identified genes with relatively high expression in the rumen, of which many appeared to be recruited from other tissues. These genes show functional enrichment in ketone body metabolism, regulation of microbial community, and epithelium absorption, which are the most prominent biological processes involved in rumen innovations. Several modes of genetic change underlying rumen functional innovations were uncovered, including coding mutations, genes newly evolved, and changes of regulatory elements. We validated that the key ketogenesis rate-limiting gene (HMGCS2) with five ruminant-specific mutations was under positive selection and exhibits higher synthesis activity than those of other mammals. Two newly evolved genes (LYZ1 and DEFB1) are resistant to Gram-positive bacteria and thereby may regulate microbial community equilibrium. Furthermore, we confirmed that the changes of regulatory elements accounted for the majority of rumen gene recruitment. These results greatly improve our understanding of rumen evolution and organ evo-devo in general.

Gene regulatory landscape of the sonic hedgehog locus in embryonic development

The organs of vertebrate species display a wide variety of morphology. A remaining challenge in evolutionary developmental biology is to elucidate how vertebrate lineages acquire distinct morphological features. Developmental programs are driven by spatiotemporal regulation of gene expression controlled by hundreds of thousands of cis-regulatory elements. Changes in the regulatory elements caused by the introduction of genetic variants can confer regulatory innovation that may underlie morphological novelties. Recent advances in sequencing technology have revealed a number of potential regulatory variants that can alter gene expression patterns. However, a limited number of studies demonstrate causal dependence between genetic and morphological changes. Regulation of Shh expression is a good model to understand how multiple regulatory elements organize tissue-specific gene expression patterns. This model also provides insights into how evolution of molecular traits, such as gene regulatory networks, lead to phenotypic novelty.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

Regeneration enhancers: Starting a journey to unravel regulatory events in tissue regeneration

Regeneration, an ability to replace lost body parts, is widespread across animal species. While mammals poorly regenerate most tissues, teleost fish and urodele amphibians possess remarkable regenerative capacity. Earlier work demonstrated that genes driving regeneration are evolutionarily conserved, indicating that a key factor in diverse tissue regeneration is not the presence or absence of regeneration-driving genes but the mechanisms controlling activation of these genes after injury. Thus, understanding the regulatory events of tissue regeneration could provide the means for unlocking latent capacities for tissue regeneration. After injury, cells undergo extensive epigenetic changes to establish new transcriptional programs for tissue regeneration. Gene transcription in eukaryotes is a complicated process that requires specific interactions between trans-acting regulators and cis-regulatory DNA elements. Among cis-regulatory elements, enhancers are essential to control precise gene expression. Recently, multiple regeneration/injury-associated enhancers have been identified in several model organisms. In this review, we highlight recently discovered regeneration/injury enhancers and their specific characteristics. We also discuss how abnormal regulation of regeneration enhancers influences animal development and physiology. Investigation of regeneration enhancers potentially allows us to begin understanding the fundamental biology of tissue regeneration and inspires new solutions for manipulating regenerative ability.

The homologies and evolutionary reduction of the pelvis and hindlimbs in snakes, with the first report of ossified pelvic vestiges in an anomalepidid (Liotyphlops beui)

We report the first example of ossified pelvic vestiges in an anomalepidid snake, Liotyplophs beui, and provide a review of the diversity of limb and pelvic elements within Serpentes. We trace the evolution, homology and reduction of the pelvic elements and hindlimbs from the oldest known snakes through to living forms. Evolutionary analysis of the pelvic and limb data shows that the most recent common ancestor of all living snakes (Serpentes) most probably retained all three pelvic elements and rudimentary hindlimbs (femoral spurs). Subsequently, there have been multiple losses of ossified pelvic and hindlimb elements and regaining of ossified pelvic elements. Reduction of the pelvis has followed different routes in the two primary groups of living snakes (scolecophidians and alethinophidians). The single remaining rod-like element in some scolecophidians is the ischium, whereas the single remaining rod-like element in many basal alethinophidians is the pubis. Notably, many basal alethinophidians share a distinctive configuration of cloacal spur (claw), femur and a sizeable pubis, which is likely to be related functionally to the widespread use of the hindlimbs in mating and courtship, rather than the usual representation of the bones as non-functional vestiges.

REforge Associates Transcription Factor Binding Site Divergence in Regulatory Elements with Phenotypic Differences between Species

Elucidating the genomic determinants of morphological differences between species is key to understanding how morphological diversity evolved. While differences in cis-regulatory elements are an important genetic source for morphological evolution, it remains challenging to identify regulatory elements involved in phenotypic differences. Here, we present Regulatory Element forward genomics (REforge), a computational approach that detects associations between transcription factor binding site divergence in putative regulatory elements and phenotypic differences between species. By simulating regulatory element evolution in silico, we show that this approach has substantial power to detect such associations. To validate REforge on real data, we used known binding motifs for eye-related transcription factors and identified significant binding site divergence in vision-impaired subterranean mammals in 1% of all conserved noncoding elements. We show that these genomic regions are significantly enriched in regulatory elements that are specifically active in mouse eye tissues, and that several of them are located near genes, which are required for eye development and photoreceptor function and are implicated in human eye disorders. Thus, our genome-wide screen detects widespread divergence of eye-regulatory elements and highlights regulatory regions that likely contributed to eye degeneration in subterranean mammals. REforge has broad applicability to detect regulatory elements that could be involved in many other phenotypes, which will help to reveal the genomic basis of morphological diversity.

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.

How a growing organismal perspective is adding new depth to integrative studies of morphological evolution

Over the past half century, the field of Evolutionary Developmental Biology, or Evo-devo, has integrated diverse fields of biology into a more synthetic understanding of morphological diversity. This has resulted in numerous insights into how development can evolve and reciprocally influence morphological evolution, as well as generated several novel theoretical areas. Although comparative by default, there remains a great gap in our understanding of adaptive morphological diversification and how developmental mechanisms influence the shape and pattern of phenotypic variation. Herein we highlight areas of research that are in the process of filling this void, and areas, if investigated more fully, that will add new insights into the diversification of morphology. At the centre of our discussion is an explicit awareness of organismal biology. Here we discuss an organismal framework that is supported by three distinct pillars. First, there is a need for Evo-devo to adopt a high-resolution phylogenetic approach in the study of morphological variation and its developmental underpinnings. Secondly, we propose that to understand the dynamic nature of morphological evolution, investigators need to give more explicit attention to the processes that generate evolutionarily relevant variation at the population level. Finally, we emphasize the need to address more thoroughly the processes that structure variation at micro- and macroevolutionary scales including modularity, morphological integration, constraint, and plasticity. We illustrate the power of these three pillars using numerous examples from both invertebrates and vertebrates to emphasize that many of these approaches are already present within the field, but have yet to be formally integrated into many research programs. We feel that the most exciting new insights will come where the traditional experimental approaches to Evo-devo are integrated more thoroughly with the principles of this organismal framework.

Reassessing the Role of Hox Genes during Vertebrate Development and Evolution

Since their discovery Hox genes have been at the core of the established models explaining the development and evolution of the vertebrate body plan as well as its paired appendages. Recent work brought new light to their role in the patterning processes along the main body axis. These studies show that Hox genes do not control the basic layout of the vertebrate body plan but carry out region-specific patterning instructions loaded on the derivatives of axial progenitors by Hox-independent processes. Furthermore, the finding that Hox clusters are embedded in functional chromatin domains, which critically impacts their expression, has significantly altered our understanding of the mechanisms of Hox gene regulation. This new conceptual framework has broadened our understanding of both limb development and the evolution of vertebrate paired appendages.