Axial patterning in snakes and caecilians: evidence for an alternative interpretation of the Hox code

Patterns of growth in the presacral vertebral column of the leopard gecko (Eublepharis macularius)

Postnatal growth patterns within the vertebral column may be informative about body proportions and regionalization. We measured femur length, lengths of all pre-sacral vertebrae, and lengths of intervertebral spaces, from radiographs of a series of 21 Eublepharis macularius, raised under standard conditions and covering most of the ontogenetic body size range. Vertebrae were grouped into cervical, sternal, and dorsal compartments, and lengths of adjacent pairs of vertebrae were summed before analysis. Femur length was included as an index of body size. Principal component analysis of the variance-covariance matrix of these data was used to investigate scaling among them. PC1 explained 94.19% of total variance, interpreted as the variance due to body size. PC1 differed significantly from the hypothetical isometric vector, indicating overall allometry. The atlas and axis vertebrae displayed strong negative allometry; the remainder of the vertebral pairs exhibited weak negative allometry, isometry or positive allometry. PC1 explained a markedly smaller amount of variance for the vertebral pairs of the cervical compartment than for the remainder of the vertebral pairs, with the exception of the final pair. The relative standard deviations of the eigenvalues from the PCAs of the three vertebral compartments indicated that the vertebrae of the cervical compartment were less strongly integrated by scaling than were the sternal or dorsal vertebrae, which did not differ greatly between themselves in their strong integration, suggesting that the growth of the cervical vertebrae is constrained by the mechanical requirements of the head. Regionalization of the remainder of the vertebral column is less clearly defined but may be associated with wave form propagation incident upon locomotion, and by locomotory changes occasioned by tail autotomy and regeneration. Femur length exhibits negative allometry relative to individual vertebral pairs and to vertebral column length, suggesting a change in locomotor requirements over the ontogenetic size range.

The skeletal ontogeny of Astatotilapia burtoni - a direct-developing model system for the evolution and development of the teleost body plan

Background

The experimental approach to the evolution and development of the vertebrate skeleton has to a large extent relied on “direct-developing” amniote model organisms, such as the mouse and the chicken. These organisms can however only be partially informative where it concerns secondarily lost features or anatomical novelties not present in their lineages. The widely used anamniotes Xenopus and zebrafish are “indirect-developing” organisms that proceed through an extended time as free-living larvae, before adopting many aspects of their adult morphology, complicating experiments at these stages, and increasing the risk for lethal pleiotropic effects using genetic strategies.

Results

Here, we provide a detailed description of the development of the osteology of the African mouthbrooding cichlid Astatotilapia burtoni, primarily focusing on the trunk (spinal column, ribs and epicentrals) and the appendicular skeleton (pectoral, pelvic, dorsal, anal, caudal fins and scales), and to a lesser extent on the cranium. We show that this species has an extremely “direct” mode of development, attains an adult body plan within 2 weeks after fertilization while living off its yolk supply only, and does not pass through a prolonged larval period.

Conclusions

As husbandry of this species is easy, generation time is short, and the species is amenable to genetic targeting strategies through microinjection, we suggest that the use of this direct-developing cichlid will provide a valuable model system for the study of the vertebrate body plan, particularly where it concerns the evolution and development of fish or teleost specific traits. Based on our results we comment on the development of the homocercal caudal fin, on shared ontogenetic patterns between pectoral and pelvic girdles, and on the evolution of fin spines as novelty in acanthomorph fishes. We discuss the differences between “direct” and “indirect” developing actinopterygians using a comparison between zebrafish and A. burtoni development.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

Deep time perspective on turtle neck evolution: chasing the Hox code by vertebral morphology

The unparalleled ability of turtle neck retraction is possible in three different modes, which characterize stem turtles, living side-necked (Pleurodira), and hidden-necked (Cryptodira) turtles, respectively. Despite the conservatism in vertebral count among turtles, there is significant functional and morphological regionalization in the cervical vertebral column. Since Hox genes play a fundamental role in determining the differentiation in vertebra morphology and based on our reconstruction of evolutionary genetics in deep time, we hypothesize genetic differences among the turtle groups and between turtles and other land vertebrates. We correlated anterior Hox gene expression and the quantifiable shape of the vertebrae to investigate the morphological modularity in the neck across living and extinct turtles. This permitted the reconstruction of the hypothetical ancestral Hox code pattern of the whole turtle clade. The scenario of the evolution of axial patterning in turtles indicates shifts in the spatial expression of HoxA-5 in relation to the reduction of cervical ribs in modern turtles and of HoxB-5 linked with a lower morphological differentiation between the anterior cervical vertebrae observed in cryptodirans. By comparison with the mammalian pattern, we illustrate how the fixed count of eight cervical vertebrae in turtles resulted from the emergence of the unique turtle shell.

Somitogenesis and Axial Development in Reptiles

Among amniote vertebrates, reptiles display the greatest variation in axial skeleton morphology. Only recently have they been used in gene expression studies of somitogenesis , challenging previous assumptions about the segmentation clock and axial patterning. An increasing number of reptile genomes and transcriptomes are becoming available as next-generation sequencing becomes more affordable. Information regarding gene sequence and structure can be used to design and synthesize labeled riboprobes by in vitro transcription for gene expression analysis by in situ hybridization, thus, enabling the characterization of spatial and temporal expression patterns of genes involved in somitogenesis, a topic of great interest within evolutionary developmental studies of vertebrates.

Evolutionary trajectories of snake genes and genomes revealed by comparative analyses of five-pacer viper

Snakes have numerous features distinctive from other tetrapods and a rich history of genome evolution that is still obscure. Here, we report the high-quality genome of the five-pacer viper, Deinagkistrodon acutus, and comparative analyses with other representative snake and lizard genomes. We map the evolutionary trajectories of transposable elements (TEs), developmental genes and sex chromosomes onto the snake phylogeny. TEs exhibit dynamic lineage-specific expansion, and many viper TEs show brain-specific gene expression along with their nearby genes. We detect signatures of adaptive evolution in olfactory, venom and thermal-sensing genes and also functional degeneration of genes associated with vision and hearing. Lineage-specific relaxation of functional constraints on respective Hox and Tbx limb-patterning genes supports fossil evidence for a successive loss of forelimbs then hindlimbs during snake evolution. Finally, we infer that the ZW sex chromosome pair had undergone at least three recombination suppression events in the ancestor of advanced snakes. These results altogether forge a framework for our deep understanding into snakes' history of molecular evolution.

Snakes have many characteristics that distinguish them from their relatives. Here, Yin et al. sequence the genome of the five-pacer viper, Deinagkistrodon acutus, and use comparative genomic analyses to elucidate the evolution of transposable elements, developmental genes and sex chromosomes in snakes.

Differential growth of body segments explains ontogenetic shifts in organ position for the Diamondback Water Snake (Nerodia rhombifer)

As snakes grow, their organs move anteriorly relative to body size. We explored a developmental explanation for the ontogenetic shift in the relative position of internal organs for snakes using the Diamondback Water Snake (Nerodia rhombifer (Hallowell, 1852)). With age, this water snake’s heart, liver, small intestine, and right kidney move anteriorly by 2.5–5.0 percentage points of snout–vent length. The number of precaudal vertebrae did not vary due to size or sex. The anterior edge of the heart, liver, small intestine, and right kidney were typically aligned within a span of 4–8 vertebrae that likewise did not differ as a function of size or sex. Snakes exhibited a positive relationship between the number of precaudal vertebrae and the vertebra number aligned with each organ. Total length, centrum length, centrum width, ball width, height, and mass of eight vertebrae sampled at consistent vertebral number revealed that vertebrae in the middle region of the body grow at a greater rate than vertebrae at the anterior or distal ends of the body. For N. rhombifer, the observed forward shift in relative organ positions is the product of regional differences in the growth of body segments. Predictably, these differences arise from a developmental program generated by the differential expression of Hox genes.

Of Mice and Snakes: A Tail of Oct4

The vertebrate axial skeleton comprises regions of specialized vertebrae, which vary in length between lineages. Aires et al. (2016) uncover a key role for Oct4 in determining trunk length in mice. Additionally, a heterochronic shift in Oct4 expression may underlie the extreme elongation of the trunk in snakes.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

New insights into the vertebral Hox code of archosaurs

Variation in axial formulae (i.e., number and identity of vertebrae) is an important feature in the evolution of vertebrates. Vertebrae at different axial positions exhibit a region-specific morphology. Key determinants for the establishment of particular vertebral shapes are the highly conserved Hox genes. Here, we analyzed Hox gene expression in the presacral vertebral column in the Nile crocodile in order to complement and extend a previous examination in the alligator and thus establish a Hox code for the axial skeleton of crocodilians in general. The newly determined expression of HoxA-4, C-5, B-7, and B-8 all revealed a crocodilian-specific pattern. HoxA-4 and HoxC-5 characterize cervical morphologies and the latter furthermore is associated with the position of the forelimb relative to the axial skeleton. HoxB-7 and HoxB-8 map exclusively to the dorsal vertebral region. The resulting expression patterns of these two Hox genes is the first description of their exact expression in the archosaurian embryo. Our comparative analyses of the Hox code in several amniote taxa provide new evidence that evolutionary differences in the axial skeleton correspond to changes in Hox gene expression domains. We detect two general processes: (i) expansion of a Hox gene's expression domain as well as (ii) a shift of gene expression. We infer that the ancestral archosaur Hox code may have resembled that of the crocodile. In association with the evolution of morphological traits, it may have been modified to patterns that can be observed in birds.

Hox genes and evolution

Hox proteins are a deeply conserved group of transcription factors originally defined for their critical roles in governing segmental identity along the antero-posterior (AP) axis in Drosophila. Over the last 30 years, numerous data generated in evolutionarily diverse taxa have clearly shown that changes in the expression patterns of these genes are closely associated with the regionalization of the AP axis, suggesting that Hox genes have played a critical role in the evolution of novel body plans within Bilateria. Despite this deep functional conservation and the importance of these genes in AP patterning, key questions remain regarding many aspects of Hox biology. In this commentary, we highlight recent reports that have provided novel insight into the origins of the mammalian Hox cluster, the role of Hox genes in the generation of a limbless body plan, and a novel putative mechanism in which Hox genes may encode specificity along the AP axis. Although the data discussed here offer a fresh perspective, it is clear that there is still much to learn about Hox biology and the roles it has played in the evolution of the Bilaterian body plan.

Molecular evolution of HoxA13 and the multiple origins of limbless morphologies in amphibians and reptiles

Developmental processes and their results, morphological characters, are inherited through transmission of genes regulating development. While there is ample evidence that cis-regulatory elements tend to be modular, with sequence segments dedicated to different roles, the situation for proteins is less clear, being particularly complex for transcription factors with multiple functions. Some motifs mediating protein-protein interactions may be exclusive to particular developmental roles, but it is also possible that motifs are mostly shared among different processes. Here we focus on HoxA13, a protein essential for limb development. We asked whether the HoxA13 amino acid sequence evolved similarly in three limbless clades: Gymnophiona, Amphisbaenia and Serpentes. We explored variation in ω (dN/dS) using a maximum-likelihood framework and HoxA13sequences from 47 species. Comparisons of evolutionary models provided low ω global values and no evidence that HoxA13 experienced relaxed selection in limbless clades. Branch-site models failed to detect evidence for positive selection acting on any site along branches of Amphisbaena and Gymnophiona, while three sites were identified in Serpentes. Examination of alignments did not reveal consistent sequence differences between limbed and limbless species. We conclude that HoxA13 has no modules exclusive to limb development, which may be explained by its involvement in multiple developmental processes.

Tetrapod axial evolution and developmental constraints; Empirical underpinning by a mouse model

The tetrapod vertebral column has become increasingly complex during evolution as an adaptation to a terrestrial life. At the same time, the evolution of the vertebral formula became subject to developmental constraints acting on the size of the cervical and thoraco-lumbar regions. In the course of our studies concerning the evolution of Hox gene regulation, we produced a transgenic mouse model expressing fish Hox genes, which displayed a reduced number of thoraco-lumbar vertebrae and concurrent sacral homeotic transformations. Here, we analyze this mutant stock and conclude that the ancestral, pre-tetrapodial Hox code already possessed the capacity to induce vertebrae with sacral characteristics. This suggests that alterations in the interpretation of the Hox code may have participated to the evolution of this region in tetrapods, along with potential modifications of the HOX proteins themselves. With its reduced vertebral number, this mouse stock violates a previously described developmental constraint, which applies to the thoraco-lumbar region. The resulting offset between motor neuron morphology, vertebral patterning and the relative positioning of hind limbs illustrates that the precise orchestration of the Hox-clock in parallel with other ontogenetic pathways places constraints on the evolvability of the body plan.

•

A transgenic mouse line expressing fish Hox genes has anterior homeotic transformations.

•

Fish Hox genes are capable of inducing tetrapod specific vertebral characters.

•

A sacral Hox-code influences adult hindlimb position, yet not the position of limb budding.

Correlation between Hox code and vertebral morphology in archosaurs

The relationship between developmental genes and phenotypic variation is of central interest in evolutionary biology. An excellent example is the role of Hox genes in the anteroposterior regionalization of the vertebral column in vertebrates. Archosaurs (crocodiles, dinosaurs including birds) are highly variable both in vertebral morphology and number. Nevertheless, functionally equivalent Hox genes are active in the axial skeleton during embryonic development, indicating that the morphological variation across taxa is likely owing to modifications in the pattern of Hox gene expression. By using geometric morphometrics, we demonstrate a correlation between vertebral Hox code and quantifiable vertebral morphology in modern archosaurs, in which the boundaries between morphological subgroups of vertebrae can be linked to anterior Hox gene expression boundaries. Our findings reveal homologous units of cervical vertebrae in modern archosaurs, each with their specific Hox gene pattern, enabling us to trace these homologies in the extinct sauropodomorph dinosaurs, a group with highly variable vertebral counts. Based on the quantifiable vertebral morphology, this allows us to infer the underlying genetic mechanisms in vertebral evolution in fossils, which represents not only an important case study, but will lead to a better understanding of the origin of morphological disparity in recent archosaur vertebral columns.

Distinct patterns of desynchronized limb regression in malagasy scincine lizards (squamata, scincidae)

Scincine lizards in Madagascar form an endemic clade of about 60 species exhibiting a variety of ecomorphological adaptations. Several subclades have adapted to burrowing and convergently regressed their limbs and eyes, resulting in a variety of partial and completely limbless morphologies among extant taxa. However, patterns of limb regression in these taxa have not been studied in detail. Here we fill this gap in knowledge by providing a phylogenetic analysis of DNA sequences of three mitochondrial and four nuclear gene fragments in an extended sampling of Malagasy skinks, and microtomographic analyses of osteology of various burrowing taxa adapted to sand substrate. Based on our data we propose to (i) consider Sirenoscincus Sakata & Hikida, 2003, as junior synonym of Voeltzkowia Boettger, 1893; (ii) resurrect the genus name Grandidierina Mocquard, 1894, for four species previously included in Voeltzkowia; and (iii) consider Androngo Brygoo, 1982, as junior synonym of Pygomeles Grandidier, 1867. By supporting the clade consisting of the limbless Voeltzkowia mira and the forelimb-only taxa V. mobydick and V. yamagishii, our data indicate that full regression of limbs and eyes occurred in parallel twice in the genus Voeltzkowia (as hitherto defined) that we consider as a sand-swimming ecomorph: in the Voeltzkowia clade sensu stricto the regression first affected the hindlimbs and subsequently the forelimbs, whereas the Grandidierina clade first regressed the forelimbs and subsequently the hindlimbs following the pattern prevalent in squamates. Timetree reconstructions for the Malagasy Scincidae contain a substantial amount of uncertainty due to the absence of suitable primary fossil calibrations. However, our preliminary reconstructions suggest rapid limb regression in Malagasy scincids with an estimated maximal duration of 6 MYr for a complete regression in Paracontias, and 4 and 8 MYr respectively for complete regression of forelimbs in Grandidierina and hindlimbs in Voeltzkowia.

Evolution of the snake body form reveals homoplasy in amniote Hox gene function

Hox genes regulate regionalization of the axial skeleton in vertebrates, and changes in their expression have been proposed to be a fundamental mechanism driving the evolution of new body forms. The origin of the snake-like body form, with its deregionalized pre-cloacal axial skeleton, has been explained as either homogenization of Hox gene expression domains, or retention of standard vertebrate Hox domains with alteration of downstream expression that suppresses development of distinct regions. Both models assume a highly regionalized ancestor, but the extent of deregionalization of the primaxial domain (vertebrae, dorsal ribs) of the skeleton in snake-like body forms has never been analysed. Here we combine geometric morphometrics and maximum-likelihood analysis to show that the pre-cloacal primaxial domain of elongate, limb-reduced lizards and snakes is not deregionalized compared with limbed taxa, and that the phylogenetic structure of primaxial morphology in reptiles does not support a loss of regionalization in the evolution of snakes. We demonstrate that morphometric regional boundaries correspond to mapped gene expression domains in snakes, suggesting that their primaxial domain is patterned by a normally functional Hox code. Comparison of primaxial osteology in fossil and modern amniotes with Hox gene distributions within Amniota indicates that a functional, sequentially expressed Hox code patterned a subtle morphological gradient along the anterior-posterior axis in stem members of amniote clades and extant lizards, including snakes. The highly regionalized skeletons of extant archosaurs and mammals result from independent evolution in the Hox code and do not represent ancestral conditions for clades with snake-like body forms. The developmental origin of snakes is best explained by decoupling of the primaxial and abaxial domains and by increases in somite number, not by changes in the function of primaxial Hox genes.

Hox gene regulation and timing in embryogenesis

Hox genes are critical regulators of embryonic development in bilaterian animals. They exhibit a unique mode of transcriptional regulation where the position of the genes along the chromosome corresponds to the time and place of their expression during development. The sequential temporal activation of these genes in the primitive streak helps determining their subsequent pattern of expression along the anterior-posterior axis of the embryo, yet the precise correspondence between these two collinear processes is not fully understood. In addition, vertebrate Hox genes evolved similar modes of regulation along secondary body axes, such as the developing limbs. We review the current understanding of the mechanisms operating during activation, maintenance and silencing of Hox gene expression in these various contexts, and discuss the evolutionary significance of their genomic organization.

Conservation and divergence of regulatory strategies at Hox Loci and the origin of tetrapod digits

During development, expression of the Hoxa and Hoxd genes in zebrafish fins and mouse limbs are regulated via a conserved chromatin structure. However, zebrafish lack certain regulatory elements required to produce digits, revealing that radials—the fin's bony elements—are likely not homologous to tetrapod digits.

The evolution of tetrapod limbs from fish fins enabled the conquest of land by vertebrates and thus represents a key step in evolution. Despite the use of comparative gene expression analyses, critical aspects of this transformation remain controversial, in particular the origin of digits. Hoxa and Hoxd genes are essential for the specification of the different limb segments and their functional abrogation leads to large truncations of the appendages. Here we show that the selective transcription of mouse Hoxa genes in proximal and distal limbs is related to a bimodal higher order chromatin structure, similar to that reported for Hoxd genes, thus revealing a generic regulatory strategy implemented by both gene clusters during limb development. We found the same bimodal chromatin architecture in fish embryos, indicating that the regulatory mechanism used to pattern tetrapod limbs may predate the divergence between fish and tetrapods. However, when assessed in mice, both fish regulatory landscapes triggered transcription in proximal rather than distal limb territories, supporting an evolutionary scenario whereby digits arose as tetrapod novelties through genetic retrofitting of preexisting regulatory landscapes. We discuss the possibility to consider regulatory circuitries, rather than expression patterns, as essential parameters to define evolutionary synapomorphies.

Our upper limbs differ from fish fins, notably by their subdivision into arm and hand regions, which are separated by a complex articulation, the wrist. The development of this anatomy is associated with two distinct waves of expression of the Hoxa and Hoxd genes during development. Would such a shared expression pattern be sufficient to infer homology between fish fins and mouse limbs? We investigated this question here, looking at whether the two phases of Hox gene transcription that are observed during tetrapod limb development also occur during zebrafish fin development. We find the answer is “not quite.” For although the mechanisms that regulate the expression of Hoxa and Hoxd are comparable between zebrafish fins and mouse limbs, when the genomic regions that regulate Hox gene expression in fish fins are introduced into transgenic mice, they trigger Hox gene expression in only the proximal limb segment (the segment nearest the body) and not in the presumptive digits. We conclude that although fish have the Hox regulatory toolkit to produce digits, this potential is not utilized as it is in tetrapods, and as a result we propose that fin radials—the bony elements of fins—are not homologous to tetrapod digits.

Evolution of transcriptional enhancers and animal diversity

Deciphering the genetic bases that drive animal diversity is one of the major challenges of modern biology. Although four decades ago it was proposed that animal evolution was mainly driven by changes in cis-regulatory DNA elements controlling gene expression rather than in protein-coding sequences, only now are powerful bioinformatics and experimental approaches available to accelerate studies into how the evolution of transcriptional enhancers contributes to novel forms and functions. In the introduction to this Theme Issue, we start by defining the general properties of transcriptional enhancers, such as modularity and the coexistence of tight sequence conservation with transcription factor-binding site shuffling as different mechanisms that maintain the enhancer grammar over evolutionary time. We discuss past and current methods used to identify cell-type-specific enhancers and provide examples of how enhancers originate de novo, change and are lost in particular lineages. We then focus in the central part of this Theme Issue on analysing examples of how the molecular evolution of enhancers may change form and function. Throughout this introduction, we present the main findings of the articles, reviews and perspectives contributed to this Theme Issue that together illustrate some of the great advances and current frontiers in the field.

Role of a polymorphism in a Hox/Pax-responsive enhancer in the evolution of the vertebrate spine

Patterning of the vertebrate skeleton requires the coordinated activity of Hox genes. In particular, Hox10 proteins are essential to set the transition from thoracic to lumbar vertebrae because of their rib-repressing activity. In snakes, however, the thoracic region extends well into Hox10-expressing areas of the embryo, suggesting that these proteins are unable to block rib formation. Here, we show that this is not a result of the loss of rib-repressing properties by the snake proteins, but rather to a single base pair change in a Hox/Paired box (Pax)-responsive enhancer, which prevents the binding of Hox proteins. This polymorphism is also found in Paenungulata, such as elephants and manatees, which have extended rib cages. In vivo, this modified enhancer failed to respond to Hox10 activity, supporting its role in the extension of rib cages. In contrast, the enhancer could still interact with Hoxb6 and Pax3 to promote rib formation. These results suggest that a polymorphism in the Hox/Pax-responsive enhancer may have played a role in the evolution of the vertebrate spine by differently modulating its response to rib-suppressing and rib-promoting Hox proteins.

Retinoic acid signaling in spinal cord development

Retinoic acid (RA) is an important signaling molecule mediating intercellular communication through vertebrate development. Here, we present and discuss recent information on the roles of the RA signaling pathway in spinal cord development. RA is an important player in the patterning and definition of the spinal cord territory from very early stages of development, even before the appearance of the neural plate and further serves a role in the patterning of the spinal cord both along the dorsoventral and anteroposterior axes, particularly in the promotion of neuronal differentiation. It is thus required to establish a variety of neuronal cell types at specific positions of the spinal cord. The main goal of this review is to gather information from vertebrate models, including fish, frogs, chicken and mice, and to put this information in a comparative context in an effort to visualize how the RA pathway was incorporated into the evolving vertebrate spinal cord and to identify mechanisms that are both common and different in the various vertebrate models. In doing so, we try to reconstruct how spinal cord development has been regulated by the RA signaling cascade through vertebrate diversification, highlighting areas which require further studies to obtain a better understanding of the evolutionary events that shaped this structure in the vertebrate lineage.

Salamander Hox clusters contain repetitive DNA and expanded non-coding regions: a typical Hox structure for non-mammalian tetrapod vertebrates?

Hox genes encode transcription factors that regulate embryonic and post-embryonic developmental processes. The expression of Hox genes is regulated in part by the tight, spatial arrangement of conserved coding and non-coding sequences. The potential for evolutionary changes in Hox cluster structure is thought to be low among vertebrates; however, recent studies of a few non-mammalian taxa suggest greater variation than originally thought. Using next generation sequencing of large genomic fragments (>100 kb) from the red spotted newt (Notophthalamus viridescens), we found that the arrangement of Hox cluster genes was conserved relative to orthologous regions from other vertebrates, but the length of introns and intergenic regions varied. In particular, the distance between hoxd13 and hoxd11 is longer in newt than orthologous regions from vertebrate species with expanded Hox clusters and is predicted to exceed the length of the entire HoxD clusters (hoxd13–hoxd4) of humans, mice, and frogs. Many repetitive DNA sequences were identified for newt Hox clusters, including an enrichment of DNA transposon-like sequences relative to non-coding genomic fragments. Our results suggest that Hox cluster expansion and transposon accumulation are common features of non-mammalian tetrapod vertebrates.

From lizard to snake; behind the evolution of an extreme body plan

The elongated, snake-like skeleton, as it has convergently evolved in numerous reptilian and amphibian lineages, is from a developmental biologist’s point of view amongst the most fascinating anatomical peculiarities in the animal kingdom. This type of body plan is characterized by a greatly increased number of vertebrae, a reduction of skeletal regionalization along the primary body axis and loss of the limbs. Recent studies conducted on both mouse and snakes now hint at how changes inside the gene regulatory circuitries of the Hox genes and the somitogenesis clock likely underlie these striking departures from standard tetrapod morphology, suggesting scenarios by which snakes and other elongated species may have evolved from more ordinarily bodied ancestors.

Developmental diversity of amphibians

The current model amphibian, Xenopus laevis, develops rapidly in water to a tadpole which metamorphoses into a frog. Many amphibians deviate from the X. laevis developmental pattern. Among other adaptations, their embryos develop in foam nests on land or in pouches on their mother's back or on a leaf guarded by a parent. The diversity of developmental patterns includes multinucleated oogenesis, lack of RNA localization, huge non-pigmented eggs, and asynchronous, irregular early cleavages. Variations in patterns of gastrulation highlight the modularity of this critical developmental period. Many species have eliminated the larva or tadpole and directly develop to the adult. The wealth of developmental diversity among amphibians coupled with the wealth of mechanistic information from X. laevis permit comparisons that provide deeper insights into developmental processes.

Primitive duplicate Hox clusters in the European eel's genome

The enigmatic life cycle and elongated body of the European eel (Anguilla anguilla L., 1758) have long motivated scientific enquiry. Recently, eel research has gained in urgency, as the population has dwindled to the point of critical endangerment. We have assembled a draft genome in order to facilitate advances in all provinces of eel biology. Here, we use the genome to investigate the eel's complement of the Hox developmental transcription factors. We show that unlike any other teleost fish, the eel retains fully populated, duplicate Hox clusters, which originated at the teleost-specific genome duplication. Using mRNA-sequencing and in situ hybridizations, we demonstrate that all copies are expressed in early embryos. Theories of vertebrate evolution predict that the retention of functional, duplicate Hox genes can give rise to additional developmental complexity, which is not immediately apparent in the adult. However, the key morphological innovation elsewhere in the eel's life history coincides with the evolutionary origin of its Hox repertoire.

Evolving possibilities: Post-embryonic axial elongation in salamanders with biphasic (Eurcyea cirrigera, E. longicauda, E. quadridigitata) and paedomorphic life cycles (Eurycea nana and Ambystoma mexicanum). Submitted Acta Zoologica

Typically the number of vertebrae an organism will have post-embryonically is determined during embryogenesis via the development of paired somites. Our research investigates the phenomenon of post-embryonic vertebral addition in salamander tails. We describe body and tail growth, and patterns of postsacral vertebral addition and elongation in context with caudal morphology for four plethodontids (Eurycea), and one ambystomatid. Eurycea nana and A. mexicanum have paedomorphic life cycles; E. cirrigera, E. longicauda and E. quadridigitata are biphasic. Specimens were collected, borrowed and/or purchased, and cleared and stained for bone and cartilage. Data collected include snout-vent length (SVL), tail length (TL), vertebral counts and centrum lengths. Eurycea species with biphasic life cycles had TLs that surpassed SVL following metamorphosis. Tails in paedomorphic species elongated but rarely exceeded body length. Larger TLs were associated with more vertebrae and longer vertebrae in all species. We observed that rates of postsacral vertebral addition varied little amongst species. Regional variation along the tail becomes prominent following metamorphosis in biphasic developers. In all species vertebrae in the posterior one-half of the tail taper towards the tip. We suggest a developmental link might exist between the ability to continually add vertebrae and regeneration in salamanders.

Analysis of cervical ribs in a series of human fetuses

In humans, an increasing body of evidence has linked the frequency of cervical ribs to stillbirths, other malformations and early childhood cancers. However, the frequency of cervical ribs in a putatively healthy fetal population is not sufficiently known to assess the actual medical risks of these prenatal findings. We therefore analyzed the presence of skeletal anomalies in a series of 199 electively aborted fetuses, which were whole-mount stained with alizarin red specific for skeletal tissues. Results show that approximately 40% of the fetuses had cervical ribs, even though external congenital abnormalities such as craniofacial and limb defects were absent. A literature overview indicates that the observed frequency of cervical ribs is comparable to results previously obtained for deceased fetuses with no or minor congenital anomalies, and higher than expected for healthy fetuses. This unexpected result can probably in part be explained by a higher detection rate of small cervical ribs when using alizarin red staining instead of radiographs. Additionally, studies in the literature suggest that the size of a cervical rib may indicate the severity of abnormalities, but this possibility requires further research. Anomalies of the axial skeleton are known to be caused by a disturbance of early development, which alters Hox gene expression, but in this study the origin of the stress could not be verified as maternal medical data were not available. The co-occurrence of rudimentary or absent 12th ribs in 23.6% of the cases with cervical ribs indicates that in approximately 8% of the fetuses a homeotic shift occurred over a larger part of the vertebral column. This suggests that the expression of multiple Hox genes may have been affected in these fetuses. Together, the high incidence of cervical ribs and also their co-occurrence with rudimentary or absent 12th ribs suggests that there may have been a disturbance of early development such that the studied fetuses are probably not informative about the general population. Future studies determining the frequency of cervical ribs in a more healthy fetal population are therefore needed to evaluate their potential as an indicator of medical risks.

Evolution of gene regulatory networks controlling body plan development

Evolutionary change in animal morphology results from alteration of the functional organization of the gene regulatory networks (GRNs) that control development of the body plan. A major mechanism of evolutionary change in GRN structure is alteration of cis-regulatory modules that determine regulatory gene expression. Here we consider the causes and consequences of GRN evolution. Although some GRN subcircuits are of great antiquity, other aspects are highly flexible and thus in any given genome more recent. This mosaic view of the evolution of GRN structure explains major aspects of evolutionary process, such as hierarchical phylogeny and discontinuities of paleontological change.

Breaking evolutionary and pleiotropic constraints in mammals: On sloths, manatees and homeotic mutations

Background

Mammals as a rule have seven cervical vertebrae, except for sloths and manatees. Bateson proposed that the change in the number of cervical vertebrae in sloths is due to homeotic transformations. A recent hypothesis proposes that the number of cervical vertebrae in sloths is unchanged and that instead the derived pattern is due to abnormal primaxial/abaxial patterning.

Results

We test the detailed predictions derived from both hypotheses for the skeletal patterns in sloths and manatees for both hypotheses. We find strong support for Bateson's homeosis hypothesis. The observed vertebral and rib patterns cannot be explained by changes in primaxial/abaxial patterning. Vertebral patterns in sloths and manatees are similar to those in mice and humans with abnormal numbers of cervical vertebrae: incomplete and asymmetric homeotic transformations are common and associated with skeletal abnormalities. In sloths the homeotic vertebral shift involves a large part of the vertebral column. As such, similarity is greatest with mice mutant for genes upstream of Hox.

Conclusions

We found no skeletal abnormalities in specimens of sister taxa with a normal number of cervical vertebrae. However, we always found such abnormalities in conspecifics with an abnormal number, as in many of the investigated dugongs. These findings strongly support the hypothesis that the evolutionary constraints on changes of the number of cervical vertebrae in mammals is due to deleterious pleitropic effects. We hypothesize that in sloths and manatees low metabolic and activity rates severely reduce the usual stabilizing selection, allowing the breaking of the pleiotropic constraints. This probably also applies to dugongs, although to a lesser extent.

Hox expression in the American alligator and evolution of archosaurian axial patterning

The avian body plan has undergone many modifications, most associated with adaptation to flight and bipedal walking. Some of these modifications may be owing to avian-specific changes in the embryonic Hox expression code. Here, we have examined Hox expression in alligator, the closest living relative of birds, and an archosaur with a more conservative body plan. Two differences in Hox expression between chick, alligator, and other tetrapods correlate with aspects of alligator or bird-specific skeletal morphology. First, absence of a thoracic subdomain of Hoxc-8 expression in alligator correlates with morphological adaptations in crocodilian thoracic segments. Second, Hoxa-5, a gene required to pattern the cervical-thoracic transition, shows unique patterns of expression in chick, alligator, and mouse, correlating with species-specific morphological patterning of this region. Given that cervical vertebral morphologies evolved independently in the bird and mammalian lineages, the underlying developmental mechanisms, including refinement of Hox expression domains, may be distinct.

A general scenario of Hox gene inventory variation among major sarcopterygian lineages

BACKGROUND

Hox genes are known to play a key role in shaping the body plan of metazoans. Evolutionary dynamics of these genes is therefore essential in explaining patterns of evolutionary diversity. Among extant sarcopterygians comprising both lobe-finned fishes and tetrapods, our knowledge of the Hox genes and clusters has largely been restricted in several model organisms such as frogs, birds and mammals. Some evolutionary gaps still exist, especially for those groups with derived body morphology or occupying key positions on the tree of life, hindering our understanding of how Hox gene inventory varied along the sarcopterygian lineage.

RESULTS

We determined the Hox gene inventory for six sarcopterygian groups: lungfishes, caecilians, salamanders, snakes, turtles and crocodiles by comprehensive PCR survey and genome walking. Variable Hox genes in each of the six sarcopterygian group representatives, compared to the human Hox gene inventory, were further validated for their presence/absence by PCR survey in a number of related species representing a broad evolutionary coverage of the group. Turtles, crocodiles, birds and placental mammals possess the same 39 Hox genes. HoxD12 is absent in snakes, amphibians and probably lungfishes. HoxB13 is lost in frogs and caecilians. Lobe-finned fishes, amphibians and squamate reptiles possess HoxC3. HoxC1 is only present in caecilians and lobe-finned fishes. Similar to coelacanths, lungfishes also possess HoxA14, which is only found in lobe-finned fishes to date. Our Hox gene variation data favor the lungfish-tetrapod, turtle-archosaur and frog-salamander relationships and imply that the loss of HoxD12 is not directly related to digit reduction.

CONCLUSIONS

Our newly determined Hox inventory data provide a more complete scenario for evolutionary dynamics of Hox genes along the sarcopterygian lineage. Limbless, worm-like caecilians and snakes possess similar Hox gene inventories to animals with less derived body morphology, suggesting changes to their body morphology are likely due to other modifications rather than changes to Hox gene numbers. Furthermore, our results provide basis for future sequencing of the entire Hox clusters of these animals.

Hox genes and regional patterning of the vertebrate body plan

Several decades have passed since the discovery of Hox genes in the fruit fly Drosophila melanogaster. Their unique ability to regulate morphologies along the anteroposterior (AP) axis (Lewis, 1978) earned them well-deserved attention as important regulators of embryonic development. Phenotypes due to loss- and gain-of-function mutations in mouse Hox genes have revealed that the spatio-temporally controlled expression of these genes is critical for the correct morphogenesis of embryonic axial structures. Here, we review recent novel insight into the modalities of Hox protein function in imparting specific identity to anatomical regions of the vertebral column, and in controlling the emergence of these tissues concomitantly with providing them with axial identity. The control of these functions must have been intimately linked to the shaping of the body plan during evolution.