BACKGROUND
understanding the complex relationship between genotype and phenotype requires an integrative and interdisciplinary biological framework  <cit> .
loci known to be involved in development, morphogenesis and in the pathogenesis of congenital malformations have been identified through forward genetic approaches, genomic mapping or genome wide association studies .
concurrent with studying how malformations arise, cellular and genetic mechanisms have been uncovered which have direct consequences on our understanding of the processes governing normal development  <cit> .
examining the development of innovations and novelties  <cit>  in natural groups provides an understanding of how body plans have been molded by natural selection  <cit>  as optimal phenotypes through adaptive evolution.
limbs and digits were key innovations in the evolution and diversification of tetrapods .
while much has been learned about gene function during morphogenesis and differentiation of the tetrapod limb  <cit> , most studies examining non-traditional model organisms have focused on limb reduction or loss of distal elements of the autopodium .
within reptilia, the order squamata  comprises more than  <dig> lineages representing independent limb reduction toward a snake-like body form  <cit> .
however, one of the most dramatic limb skeletal modifications is found in the hands and feet of the family chamaeleonidae, the chameleons .
chamaeleonidae limbs have been modified from the generalized terrestrial tetrapod plan to exhibit an architecture highly adapted for an arboreal lifestyle.
this was facilitated through evolution of a midline autopodial cleft  and two opposable syndactylous bundles  of digits that are highly mobile .
in addition, the proximal skeleton of the autopodium  was modified through a reduction in the number of bone elements.
furthermore, those elements that remain are enlarged and form a ball-and-socket joint between the autopodium and forearm .
this specialized joint allows for greater rotation of the wrist and ankle, which is important while climbing.
limb modifications such as these facilitated an increase in precision, security and mobility in the three-dimensional chameleon arboreal environment.
while other lizards are also arboreal , chameleons differ by using zygodactyly and gripping of branches rather than depending on claws and epidermal modifications .fig.
<dig> chameleon autopodia.
a the veiled chameleon  has a laterally compressed body, prehensile tail, turreted and independent eyes, a projectile tongue and zygodactylous cleft autopdia.
b chameleon hand digits i, ii, iii are syndactylous medially with digits iv, v are bound laterally.
feet  have digits i, ii syndactylous medially while digits iii, iv, v are bound laterally.
d the desert grassland whiptail lizard  was used as for comparison due to its ‘typical’ terrestrial lizard body plan
few studies to date have focused on the development of chameleon autopodial architecture  <cit> , leaving unanswered the question of how the highly conserved tetrapod pentadactyl autopodium  <cit>  was modified during chameleon evolution.
to address this issue, we have developed and utilized the veiled chameleon , which offers key advantages over other lizard and chameleon species  <cit> .
one of the most important and distinctive developmental traits in c. calyptratus is their slow rate of embryogenesis , which facilitates the examination of fine, progressive, morphological changes  <cit> .
chameleons provide an important example of how developmental mechanisms may be modified through heterochrony and constrained by body size to produce highly adaptive morphologies.
with respect to the limb, a proximal reduction in the number of carpal and tarsal elements in adult chameleons with concurrent distal cleft formation of autopodia have facilitated the evolution of complex zygodactylous grasping in a three dimensional arboreal environment.
we collected c. calyptratus embryos at various time points to study gene expression, morphogenesis, and skeletal and histological development of autopodial elements.
we also examined  <dig> of the currently recognized  <dig> genera of chameleons and  <dig> outgroup taxa to understand the evolution and diversity of carpal and tarsal elements across this phenotypically derived clade of arboreal lizards.
a reversal to an increased number of carpal and tarsal elements together with increased independence of elements, appears to have played a significant role in the evolution of ‘true chameleons’.
such a mechanism could provide more components by which to increase wrist flexion, a biomechanical advantage in a complex arboreal environment.
methods
embryos
veiled chameleon eggs were collected at oviposition in the stowers institute for medical research’s reptiles and aquatics facility and at the la sierra university animal facility.
eggs were subsequently incubated at a constant 26 °c in deli cups with moist vermiculite.
embryos were harvested in cold 1× phosphate buffered saline  at varying stages, fixed overnight in depc-treated 4 % paraformaldehyde in 1xpbs at 4 °c and later dehydrated to absolute methanol before storing at −20 °c  <cit> .
dapi helicon imaging
whole embryos were collected in cold 1x phosphate buffered saline  and fixed overnight at room temperature with 1 μl dapi, dilactate   per 1 ml of pbs on a roller   <cit> .
embryos were placed flat in a 50 ml petri dish filled with pbs and imaged on a leica mz fliii with a carl zeiss color axiocam hrc using axiocam software.
the uv filter was used to visualize dapi fluorescence .
settings were changed to black and white in the software to generate a grayscale image of the embryo.
z-stack images were taken and subsequently processed with helicon focus  <dig>  software , which functions to create a projection of all images into a single high focus photo.
this process is analogous to doing a dapi  z-stack on a confocal, and is useful for generating high resolution images of embryos too large for visualization on a confocal microscope.
this method also overcomes some of the limitations of sem imaging in that the embryos never have to be dehydrated, can be kept at room temperature in fixative indefinitely, can be used on museum specimens, is non-destructive , and also provides simultaneous tissue resolution below the differentiating skin.
a 50 ml plastic petri dish was half-filled with 1 % agarose and left to solidify.
troughs in the agarose were created to accommodate embryos in various orientations and provide stability during imaging .
rna collection, cdna synthesis, probe development and in situ hybridization
embryos between  <dig> and  <dig> dpo  were collected in cold depc-treated 1xpbs.
limbs were pooled together and total rna was extracted with the rneasy mini kit .
complementary dna  was derived from messenger rna using the superscript® first-strand synthesis system for rt-pcr .
degenerate primers were designed for genes of interest.
gotaq® colorless mastermix  was used for pcr amplification.
reaction products were run on a 1 % agarose gel and bands of expected sizes were excised, purified using the qiaex ii gel extraction kit , and stored at −20 °c.
amplified fragments were ligated into pgem®-t easy vectors .
dna was sequenced using sp <dig> and t <dig> primers, and checked via blast  to determine the direction of the insert.
plasmids were digested using the appropriate restriction enzymes for the development of sense and antisense in situ hybridization probes using the dig rna labeling kit  .
degenerate primers used include: for sonic hedgehog  shhf: 5’-ccggcttcgactgggtntayta-3’, shhr: 5’-catgggcggtcagtggngcrtangc-3’.
for fibroblast growth factor  <dig>  fgf8f: 5’-ctmrtscgsacctaccarct-3’, fgf8r: 5’-ggsarsckcttcatgaagtg-3’.
for gremlin  <dig>  grem1f: 5’-caaggmtcagcmcaatgaytc-3’, grem1r: 5’-yttcttgggcttgcagaag-3’.
genbank accession numbers for the above amplified fragments are .
for in situ hybridization, embryos were collected in cold depc-treated 1x pbs, and fixed in 4 % paraformaldehyde overnight at 4 °c on a rocker.
specimens were then dehydrated through a graded series of methanol/pbs-depc into 100 % methanol and stored at −20 °c.
whole mount in situ hybridization was performed as previously described  <cit> .
micro-computed tomography
comparative limb osteological data was obtained by isolating limbs from adult specimens followed by scanning with a scanco vivact <dig> micro-computed tomography scanner.
dicom data was visualized using amira  <dig>  .
scanning electron microscopy
embryos were fixed overnight in  <dig>  % gluteraldehyde/2 % paraformaldehyde in 1xpbs, rinsed in water , incubated in 2 % osmium tetroxide , rinsed again in water , and dehydrated through a graded series of ethanol  for 20 min each.
critical point drying was conducted in tousimis samdri- <dig> critical point dryer.
scanning was accomplished on a hitachi tm- <dig> tabletop microscope.
alcian/alizarin skeletal preparations
staining of bone with alizarin red  and cartilage with alcian blue 8gx  was adapted from walker and kimmel  <cit> .
initially, embryos were fixed in absolute ethanol  and dehydrated/stored in absolute etoh .
specimens were stained overnight in alizarin red/alcian blue solution at room temperature with rocking.
on the following day, specimens were rinsed in water, and, depending on the size of the specimen, were incubated in 4 %, 2 %, 1 % or  <dig>  % koh until sufficient soft tissue digestion had occurred.
specimens were then cleared in 25 % glycerol in  <dig>  % koh, followed by 50 % glycerol in  <dig>  % koh before being post-fixed in 2 % formaldehyde in 50 % glycerol/1xpbs, which is also the permanent storage media for specimens.
adult museum specimen catalogue numbers  are included in the images in additional file 1: figure s <dig>  with each species represented by its unique tiff file labeled a-n. specimen measurements can be found in additional file  <dig>
h&e paraffin histology
embryos were fixed in 4 % pfa overnight, dehydrated through an ethanol series, cleared in xylene, and embedded in paraffin.
8–10 μm sections cut onto slides were deparaffinized with xylene, re-hydrated through an ethanol series into water and counterstained with hematoxylin and eosin.
whole embryo culture
eggs were incubated at 26 °c for 100 days and candled to check for viable embryos.
healthy eggs limit light passing through the egg  and also exhibit membrane vasculature, whereas infertile or unhealthy eggs appear quite clear .
pin-hole punctures with forcep tips allowed for very clear and viscous albumin/chorioallantoic fluid to leak onto the egg shell, while infertile eggs burst, releasing a creamy white yolk .
shell opening ruptures the chorioallantoic membrane , however the embryo remains within the amnion embedded within the yolk sac.
the amnion can be seen as a remnant of membranes externally .
the supine position of the embryo is associated with the constraint imposed by the umbilical arteries and veins leading toward this thin transparent membrane where the yolk sac and allantoic membranes both diverge .
care must be taken to cut through the clear membrane without severing the umbilical veins/arteries.
the egg yolk was cut to form a ‘blanket’ below the embryo while retaining the remaining chorioallantoic membrane and yolk sac vasculature .
embryos were then placed in organ culture media within a 6-well plate .
embryos were cultured  in an exo-terra thermoelectric reptile incubator  for 4 days at 26 °c under high humidity.
media was changed daily.
media for culture
culture media consisted of  <dig>  ml of fetal bovine serum ,  <dig>  ml dmem/f- <dig> 50/ <dig> with l-glutamine , 10 μl of 50 mg/ml ascorbic acid for final 50 μm/ml concentration .
200 μl of pen-strep  was added to limit bacterial contamination.
nile blue staining
embryos were incubated in 1: <dig> dilution of a  <dig>  % stock solution of nile blue in culture media for 40 min in a 26 °c incubator, followed by two 30 min washes in pbs with rocking at 4 °c prior to photographing.
ldn <dig>  treatment
dorsomorphin  was added to a final concentration of 15 μm in culture media .
animal research
all procedures were formally approved by the institutional animal care and use committee at the stowers institute for medical research and the animal care and use committee at la sierra university.
no wild caught animals were used in this study.
all embryos studied were produced in house.
RESULTS
interdigital clefting is independent of a robust aer
to begin to understand the mechanisms governing the modified hand and foot morphology of chameleons, we compared autopodial development in the veiled chameleon to the desert grassland whiptail lizard , which serves as a reference species due to its overall ‘typical’ lizard limb and body plan.
in lizards, it takes approximately 10 days for limb development to progress from a round digital plate to the typical pentadactyl pattern which coincides with the onset of interdigital cell death  .
in contrast, autopodial morphogenesis in c. calyptratus lasts approximately 14– <dig> days from the rounded paddle stage through cleft completion.
phenotypic divergence between these species commences at the late round digital plate stage when distal flattening occurs in the chameleon forelimb and hindlimb .
subsequent anteroposterior expansion of the flattened region delimits the clefting domain such that forelimb and hindlimb digits differ in their syndactylous clusters, with digits grouped as    and   , respectively .fig.
<dig> comparative morphogenesis of lizard autopodia.
in aspidoscelis uniparens, the digital plate progresses through a round shape with the pentadactyl digital rays developing within and visible through the thin tissue .
subsequent thickening of digital rays and thinning of the interdigital mesenchyme  ultimately leads to loss of the interdigital tissue through apoptosis .
in the chameleon hand, a semi round digital plate develops  with subsequent distal flattening.
clefting initiates at this distally flattened location .
two syndactylous digit bundles are approximately 180° apart in the anterior-posterior plane.
the a. uniparens foot also passes through a semi round digital plate stage  which subsequently presents digit differentiation with subsequent interdigital cell death .
the chameleon foot undergoes the same stages as in the hand , but terminates with a different complement of digits in syndactyly .
dpo = days post oviposition; a = anterior; p = posterior
current knowledge of cleft formation in autopodia comes from examination of avian and mammalian systems as well as linkage studies in humans, with all work supporting a failure to maintain the integrity of the apical ectodermal ridge , particularly along the distal midline .
loss of the aer inhibits not only distal outgrowth of the limb but also leads to the loss or splitting of the digital rays.
to determine whether or not aer destabilization is the primary factor driving cleft formation, we examined the aer of a. uniparens in order to first understand how a ‘typical’ aer is structured in pentadactyl lizard autopodia .
at  <dig> dpo , a robust distal ectodermal thickening  is present at the dorso-ventral boundary of the forelimbs and hindlimbs.
at  <dig> dpo  the aer begins to narrow as the digital rays differentiate within the underlying mesenchyme.
by 20– <dig> dpo  the aer further narrows prior to the initiation of interdigital cell death.
by comparison, in chameleon autopodia, the aer is thicker and more prominent .
interestingly, upon distal flattening of the digital plate , commencement of distal clefting in chameleon interdigital tissue is visible despite the retainment of a very robust aer.
syndactylous digit clusters diverge toward the anterior and posterior poles, increasing the angle between clusters and further expanding the cleft.
the proximodistal size of the autopodium remains unchanged during these stages.
overall, autopodial shape shifts from a distally flattened rounded paddle to one that is rectangular shaped .
it is interesting to note that as the aer is maintained, it appears slightly thicker in the region overlying the clefting mesenchyme , which may be compensatory to maintain interdigital tissue proliferation during this midline mesenchymal expansion in concert with digit bundle realignment prior to cleft formation.fig.
<dig> sem of aspidoscelis uniparens hand and foot morphogenesis.
aspidoscelis uniparens hands and feet undergo the same stages of morphogenesis as presented in fig.
<dig> in dorsal view , respectively.
in distal view , the hand shows an apical ectodermal ridge  which is initially very robust as a distal ectodermal thickening and is present at the midline between the dorsal and ventral half of the developing limb.
at later stages, the aer narrows significantly as it spans a greater anterior-posterior domain.
in the foot , we see a similar situation with a very straight and stereotypical aer which is initially robust and later thins
fig.
<dig> sem of chamaeleo calyptratus hand and foot morphogenesis.
early stages of morphogenesis show that the chameleon hand develops as a round digital plate  which subsequently develops a distal flattening .
in dorsal view, these stages present a very robust aer, relatively larger than that present in a. uniparens  while in distal view the aer is seen as having a greater thickness along the dorsoventral midline .
significantly, during the stage at which the distal autopodium begins to flatten , the distal aer is no longer the stereotypical a-p flattened ectodermal thickening but is instead arched ventrally .
surprisingly, despite having a robust distal aer, proximal mesenchymal cleft formation has already begun .
at later stages , both the forelimb and hindlimb expand the cleft while maintaining the aer quite robust until the thickness tapers and narrows at later stages
autopodial patterning is normal despite a divergent skeleton
chameleons retain the pentadactyl tetrapod autopodial plan but differ significantly from other tetrapods due to the presence of ectrodactyly and a highly modified mesopodial skeleton.
to establish whether early changes in limb patterning influence subsequent development of the novel hand and foot morphology in this group, we characterized the spatiotemporal activity of the major signaling pathways involved in limb development as previously determined in murine and avian species.
we initially examined the spatiotemporal activity of sonic hedgehog .
shh is known to play a central role in regulating anteroposterior patterning of the digits as well as proliferation of the autopodial mesenchyme which gives rise to the digits.
furthermore, shh plays a feedback role with grem <dig> in the limb mesenchyme and with fgf <dig> in the aer to maintain distal limb outgrowth  <cit> .
to date, shh expression has only been examined during squamate limb development in lizard limbs, and it was shown that the temporal activity of shh expression in the posterior zone of polarizing activity  plays a role in determining the number of skeletal elements in the digits of element-reduced species in the scincid genus hemiergis.
in this case, shorter expression periods were observed in association with phalangeal and digit reduced taxa  <cit> .
in chameleons, the tetrapod pentadactyl plan is conserved.
the phalangeal formulae for the hands and feet of c. calyptratus are only slightly modified   <cit> .
therefore, we hypothesized that the spatiotemporal activity of shh would also be conserved in chameleons.
shh is in fact expressed in the expected posterior domain demarcating the zpa in veiled chameleons, as it is in the pentadactyl hemiergis  <cit> .
furthermore, shh is expressed in c. calyptratus throughout early limb bud development while the digital plate is rounded , but is subsequently downregulated upon distal flattening of the digital plates .fig.
<dig>
c. calyptratus autopodial sonic hedgehog  expression.
wholemount in situ hybridization indicates that shh is expressed in a similar pattern to other amniotes along the posterior margin of the developing limb within the zone of polarizing activity  in the hands  and feet .
significantly, the expression of shh is lost upon distal flattening of the limb bud in both the hands  and feet
the aer forms at the dorsoventral boundary along the distal edge of the limb bud and in mice is the site of expression for members of the fibroblast growth factor family   <cit> .
in particular, fgf <dig> is necessary and sufficient for distal outgrowth and patterning of the growing limb  <cit> .
temporal removal of the aer at various developmental stages allows for only more proximal elements to form  <cit> .
using fgf <dig> as a candidate marker for aer integrity , we find that the aer maintains normal fgf <dig> expression along the anteroposterior distal boundary across all digits as would be expected for a pentadactyl limb .
however, fgf <dig> continues to be expressed in the aer covering the distal boundary of the expanding interdigital domain which will form a cleft during anteroposterior movements of the two syndactylous digit bundles .
ultimately, fgf <dig> expression diminishes in association with interdigital mesenchymal regression .fig.
<dig>
c. calyptratus fibroblast growth factor  <dig>  expression in the apical ectodermal ridge.
in wholemount in situ hybridization, the apical ectodermal ridge expresses fgf <dig> in the highly conserved pattern as seen in tetrapods  along the distal anteroposterior boundary of the limb.
fgf <dig> expression is maintained in the aer at stages where mesenchymal clefting of the autopodia has already commenced in both the hands and feet, though at a reduced amount prior to complete cleft formation
clefting is due to increased cell death
recent work exploring digit number reduction  <cit>  supports an important role for cell survival in sculpting adult autopodial elements.
the presence of macrophages in interdigital cleft mesenchyme has led to the suggestion that apoptosis is a primary driver of the derived chameleon autopodium  <cit> .
this idea warrants further validation, as mature digit formation generally requires apoptotic removal of interdigital mesenchyme in all known tetrapod limbs  <cit> .
macrophage associated or non-associated cell death would therefore be expected to be localized to the clefting mesenchyme in most tetrapods, with the exception of taxa with interdigital webbing where cell death is inhibited .
we used nile blue staining as evidence for cell death in whole embryo cultured chameleons at stages during clefting .
while the presence of macrophages is common in the aer to aid its maintenance and extension  <cit> , we also observed increased cell death in the mesenchyme in the distal margin of the autopod.
at the onset of interdigital cell death, we observed differential nile blue staining between the major clefting domain along the distal midline relative to the interdigital mesenchyme of digits bound in syndactyly .
this reflects the comparative retention of mesenchyme between digits, compared to the cleft.
at later stages of embryogenesis, the aer covering the cleft is no longer stained by nile blue .
in contrast, anterior and posterior syndactylous domains retain nile blue staining in the aer of the hand  and foot , highlighting the requirement for an aer in mesenchymal tissue maintenance.
in ducks and bats, gremlin <dig> , an inhibitor of bmp signaling, is expressed in the interdigital mesenchyme promoting the retention of mesenchymal tissue, resulting in webbing  <cit> .
chameleons initially exhibit expression of grem <dig> in the developing limb bud mesenchyme anterior to shh in the zpa  <cit> , which is typical for limbs during outgrowth.
however, grem <dig> is subsequently downregulated in interdigital regions where syndactyly occurs .fig.
<dig> nile blue staining for cell death during autopodial morphogenesis.
during autopodial morphogenesis, nile blue granules are concentrated along the distal periphery of the developing limb within the mesenchyme as well as delimiting the apical ectodermal ridge in both the hands and feet .
the aer is also highlighted as ventrally arched in distal view  and is also clearly dorsoventrally thinner at stages when cleft formation is underway .
during final stages of cleft morphogenesis, nile blue staining is concentrated within the interdigital mesenchyme , with a significantly greater concentration in the interdigital mesenchyme between digits iii and iv in the hand  and ii and iii in the foot .
reduced nile blue is present in the tissue within syndactylous digits.
nile blue is visible in the aer on the anterior and posterior regions of the limb while absent in the cleft domain in the hand
fig.
<dig>
c. calyptratus gremlin  <dig>  expression during autopodial morphogenesis.
grem <dig> is expressed in the mesenchyme of the autopodium anterior to the zpa where a feedback loop is known to be present between grem <dig> and shh, bmps, and fgf <dig> .
grem <dig> appears to decrease its anterior expression during digital plate formation in the forelimb  and hindlimb .
during distal flattening and cleft formation, grem <dig> does not appear ectopically between syndactylous digits
members of the bone morphogenetic protein  family play diverse roles in limb morphogenesis.
during early limb development bmps are initially required for formation of the aer  <cit> , but are then negatively regulated by grem <dig> in the limb mesenchyme to allow for activation of fgfs present in the aer .
later during autopodial morphogenesis, bmp <dig>   <dig>  and  <dig> are expressed within interdigital mesenchyme prior to and during the process of programmed cell death where soft tissue is removed to delineate each digit  <cit>  .
inhibiting bmp signaling within interdigital mesenchyme through ectopic expression of bmp inhibitors such as noggin or gremlin or by dominant negative bmp receptors leads to a loss of programmed cell death in chick  <cit>  and mouse .
to functionally test the role of active bmp signaling in establishing the apoptotic program in association with cleft formation in the chameleon limb, we cultured whole embryos in vitro in the presence of ldn <dig>  a pharmacological molecule that inhibits bmp signaling by blocking smad activity .
two groups of embryos were cultured for 4 days.
one group termed ‘younger’ were collected at the distally flattened limb paddle stage of development prior to the onset of interdigital cell death .
the ‘older’ group was collected after clefting had commenced .
for ‘older’ embryos cultured at the onset of clefting, tissue loss continued  despite exposure to dorsomorphin .
however, if dorsomorphin exposure occured prior to the initiation of cell death, clefting was prevented .
however, digital ray development was still slightly impaired, and the digits also appeared smaller relative to controls .
smaller digital rays in bmp inhibited mouse limbs are due to direct inhibition of sox <dig> by fgf, and a delay in chondrification  <cit> .
blocking bmp activity maintains digital ray cells in an undifferentiated state due to continued fgf expression in the aer  <cit> .
treating mouse embryonic limbs with dorsomorphin during the period of interdigital cell death led to a reduction of smad phosphorylation in interdigital mesenchyme but did not affect apoptosis  <cit> .
however, treating mouse embryos with dorsomorphin at earlier stages prior to initiation of apoptosis, did inhibit interdigital cell death.
previous studies genetically perturbing bmp signaling in chicken and mouse embryos has also led to comparable limb phenotypes  <cit> .
thus, our study of cleft formation in chameleon embryos reveals that blocking bmp signaling can inhibit interdigital cell death in a time dependent fashion prior to the onset of programmed cell death.fig.
<dig> disruption of bmp signaling during autopod development.
dorsomorphin  at 15 μm was used to disrupt bmp signaling during whole embryo in vitro culture for 4 days.
in the ‘younger’ group of embryos  collected prior to distal flattening of the digital plate, treatment led to impaired cleft formation.
asterisks highlight reductions in cleft formation relative to controls.
b older embryos collected during the onset of clefting  had no effect from the ldn <dig> treatment over a 4 day growth period
cryptic diversity of mesopodial skeletal element number in c. calyptratus embryos
limbs and fins come in a variety of shapes and sizes, with variable numbers of skeletal elements  <cit> .
despite divergent morphologies, all autopodia are dorsoventrally flattened structures, as can be seen by juxtaposing the hands and feet of the evolutionarily distant and morphologically divergent lizard and mouse .
in instances where clefting occurs in mammals , the autopodium is bifurcated along the distal midline  and results in missing elements or a reconfigured architecture  <cit> .fig.
<dig> tetrapod autopodia are dorsoventrally flattened, except in chameleons.
chameleon autopodia  while retaining the pentadactyl complement of digits, not only have a distal cleft separating the digits iii and iv in the hand  and ii and iii in the foot  but also develop a rosette shaped cluster of metacarpals and metatarsals around a central enlarged distal carpal element.
in species  with congenital distal clefting , clefting does not alter the dorsoventrally flattened morphology
despite having a relatively normal complement of phalanges and a pentadactyl number of digits, chameleons have an extremely modified carpus and tarsus .
in addition, the proximal skeleton of the autopodium  is reduced in terms of the number of bone elements.
furthermore, those that remain are enlarged and form a ball-and-socket joint between the autopodium and forelimb .
architecturally, clefting and syndactyly are necessary for grasping, but these two phenotypic novelties are not sufficient to explain the ball-and-socket mesopodial skeleton.
thus, wrist and ankle bone homologies remain contentious relative to outgroup lizards and other tetrapods  <cit> .
previous investigations of the potential homology of derived mesopodial skeletal elements in chameleons relied on late gestation stage trioceros hoehnelii and bradypodion pumilum embryos .
similarly, comparisons to outgroup lepidosaurs also relied on late gestation stage embryos and were also based solely on the number of visible ossification centers in the mesopodia.
consequently, distal carpals 2– <dig> were interpreted to have fused to make a large central element with the proximal forelimb elements comprising the radiale, ulnare and pisiform, while in the hindlimb only distal tarsal  <dig> and a ‘tarsale proximale’ were present  <cit> .
to better understand autopodial evolution and structural diversity across chameleons, we generated a skeletal ontogenetic series of chamaeleo calyptratus in order to: 1) identify the developmental sequence of appearance of autopodial skeletal elements, 2) homologize the derived and simplified elements of the chameleon autopodium, and 3) understand what developmental events occur in the formation of c. calyptratus autopodia.
in our skeletal analysis of chamaeleo calyptratus, we identified  <dig> rows of skeletogenesis during development of the carpals :  <dig> proximal carpal elements  and  <dig> distal elements .
tarsal elements also developed from two rows :  <dig> proximal elements  and  <dig> distal elements .
interestingly, the intermedium appears to condense in association with different mesopodial elements, from the distal medial radius in the hand and from the distal medial fibula in the foot .
thus, c. calyptratus retains more mesopodial elements than previously acknowledged, with only the medial and lateral centrale of lepidosauria failing to develop additional file  <dig>  centralia are generally the last elements to form ontogenetically  <cit>  and are consistent with a heterochronic terminal truncation.
a clear series of proximal condensations is not present in the outgroup taxa aspidoscelis uniparens and pogona vitticeps, as observed by the increased and more diffuse alcian blue staining present across this domain .fig.
<dig> chondrogenic differentiation in chameleon carpus.
during distal flattening and expansion , the digital rays are already present and are lightly stained by alcian blue.
the first elements of the mesopodium to form as cartilage are the fi, dc <dig>  and metacarpals 3– <dig>  distal carpal  <dig> and  <dig> form next  with later appearance of the intermedium, metacarpal  <dig>  and ultimately the radiale and pisiform .
all embryos were at  <dig> dpo from clutchmates
fig.
<dig> chondrogenic differentiation in chameleon tarsus.
during distal flattening expansion, the first elements to appear are the ulnare and distal tarsal  <dig>  followed by distal tarsal  <dig> and metatarsals 3– <dig> and the intermedium .
metatarsal  <dig> is followed along with an appearance of the tibiale .
all embryos were at  <dig> dpo from clutchmates
fig.
<dig> outgroup lizard species autopodial chondrification.
desert grassland whiptail lizard  and bearded dragon  chondrification patterns show post-axial dominance in elements appearance along the metapterygial axis.
a. uniparens  hands show early formation of the ulnare, distal carpals  <dig> and  <dig>  metacarpals 2– <dig> and a diffuse staining in the area of the radiale, lateral centrale and intermedium a. this empty domain later fills with a large medial and lateral centrale b. the foot develops from less elements and completes with a fused metatarsal  <dig> and distal tarsal  <dig> and a large distal tarsal  <dig> c-d. a large proximal tarsal element comprising the fibulare, intermedium and tibiale + lateral centrale develops in the ankle.
in p. vitticeps
e-j, a very similar sequence of cartilage elements in the hand  develop relative to a. uniparens with the added intermediate stage of seeing the medial carpale segment from the proximal end of distal carpal  <dig>  while the lateral centrale shifts distally adjacent to dc <dig> and dc <dig>  the tarsus is more difficult to follow due to a small sample size and diffuse staining, but the proximal tarsal element is at least formed by the tibiale, intermedium and fibulare  with centralia not distinguishable
during early embryogenesis,  <dig> carpal elements and  <dig> tarsal elements form in chamaeleo calyptratus autopodia.
the hands and feet at  <dig> dpo are dorsoventrally flattened, the plesiomorphic condition for tetrapods.
during later stages of embryogenesis, fusion of proximal elements and a reorganization of the distal limb skeleton gives way to a distal cleft and zygodactyly.
in the adult manus of c. calyptratus, the  <dig> proximal elements of the first row are retained, while only  <dig> distal carpal  elements are present, dc <dig>  , dc <dig>  in the pes, the fibulare, intermedium, and tibiale fuse, leading to a single ossification center comprising the proximal tarsal element , while distally, the distal tarsals , dt <dig>  dt <dig>  dt <dig> remain as independent bony elements.
while observing the formation of the carpals, the ulnare was noted to undergo a segmentation event from its ventrolateral surface at  <dig> dpo .
subsequent analysis with later stage embryos provided evidence that this particular element later ossified as the pisiform bone.
the pisiform is a large ossified element on the posterior margin of the manus in leaf mimic chameleons .
the pisiform increases in size  as one samples higher taxa across the chameleon phylogeny  <cit> .
thus, the origin of the pisiform as a sesamoid is not supported by our study .
both the pisometacarpal ligament and the tendon of the flexor carpi ulnaris are visible with separate insertions, not enclosing the pisiform within a tendon as expected for a sesamoid in adult chameleons.
this conclusion contrasts with a recent study which also recognized the pisiform’s origin from the ulnare, but declared it as a sesamoid  <cit> .
the bony element of the posterior tarsus is widely considered a sesamoid but was not apparent until after hatching in our study of c. calyptratus .
this element clearly originates within a tendon  and supports the effect of biomechanical stress in inducing the skeletogenic program of sesamoids  <cit> .
additional file 5: figure s5a-c illustrates the segmentation of the pisiform from the ventrolateral wall of the ulnare while transverse histological sections show the pisiform segmenting from the ulnare with a shared outer condensed layer of mesenchymal cells, the prospective perichondrium .fig.
<dig> contentious skeletal elements of the chameleon mesopodium.
while broadly considered to be a sesamoid, we see the pisiform segmenting from the ventrolateral margin of the ulnare, underscoring its origin from the primary skeleton .
in ‘true chameleons’ we see what appear to be ectopic and possibly de novo elements on the posterior tarsus.
b the two distal elements are actually the tibiale and the intermedium which become shifted distally toward distal tarsals by the enlarged distal epiphysis of the tibia.
c despite the highly modified autopodium, chameleons retain distal tarsal 5
broad mesopodial skeletal diversity across chamaeleonidae
to understand how skeletal fusions and heterochronic changes during skeletogenesis in the mesopodium of chameleons have been modified by natural selection, due to the need to adapt to an arboreal lifestyle, we sampled species representing  <dig> of the currently recognized  <dig> genera of extant chameleons  <cit> .
palleon was examined from an embedded 3d microct isosurface reconstruction in glaw  <cit> .
despite an overall conserved body plan, chameleon body size ranges from the large tree adapted furcifer oustaleti and calumma parsonii, at approximately 68 cm total length, to the recently described brookesia micra, whose maximum total length is less than 30 mm  <cit> , a striking  <dig> fold difference.
within the family chamaeleonidae,  <dig> genera diverged early  in madagascar with all other genera forming the sister clade.
the two malagasy genera and two basal genera found within the sister clade are known as the leaf chameleons as they are leaf mimics, whose niche encompasses ground cover and low shrubs in madagascar and africa  <cit> .
classically, morphological descriptions of chameleon hand and foot skeletons came from two derived and highly arboreal genera, trioceros and chamaeleo  <cit> , thus biasing our knowledge of morphological diversity of autopodia across chameleons.
this has inadvertently led to an underappreciation of skeletal element number with respect to carpals and tarsals in chameleons generally.
our comparative analysis illustrates that the most basal lineages of chameleons  maintain the least number of independent carpal and tarsal elements as adults .
in the forelimb mesopodium, the lowest number of independent skeletal elements are found in the genera brookesia, rieppeleon and rhampholeon.
members of these genera have a mesopodial formula of 3/ <dig>  representing a /).
all chameleons examined lack the middle row comprised of two centralia which are found in outgroup lepidosaurs.
an isosurface microct 3d rendering of palleon lolontany grossly presents with the same autopodial morphology as brookesia  <cit> .
the foot mesopodial skeleton was greatly simplified in these three genera.
the most reduced architecture is found in rieppeleon brevicaudatus, having a single distal tarsal and two proximal elements of which one was a sesamoid within the fibularis brevis tendon, whereas the other was a fused tibiale + intermedium + fibulare.
surprisingly, in rieppeleon the proximal face of all metatarsals fuse with adjacent metatarsals to form a single crown surrounding the distal tarsal element .
brookesia and palleon were found to be the only taxa sampled lacking the fibularis brevis sesamoid ossification.
this may be representative of basal chameleons  <cit>  and/or locomotory mode.fig.
<dig> histological examination of mesopodial mesenchymal and cartilaginous condensations.
to complement the wholemount alcian cartilage stains, standard h&e paraffin sections were made to look for condensations which may not have shown due to a failure to reach the chondrification stage.
in the carpals, we were able to find additional support for elements such as the intermedium, radiale and ulnare .
in addition, we were able to identify a condensation which appears to be small and transient as distal carpal  <dig>  which does not chondrify and has a significant contribution of mesenchyme from the preaxial side of the limb.
in the tarsals, we see a significant amount of shared mesenchymal cells between the distal tarsal  <dig>  fibulare, intermedium and the tibiale , with the tibiale and intermedium always in closer association and smaller
from this study we can conclude that increased body length in chameleons >20 cm accompanies a corresponding increase in the number of carpal elements, from 3/ <dig>  in leaf mimic chameleons, to 3/ <dig>  3/ <dig>  4/ <dig> and 4/ <dig> in the phylogenetically derived ‘true chameleons’.
hindlimb mesopodial skeletal elements also increase from 1/ <dig>  [proximal  <dig> rows with sesamoid in parentheses) in leaf mimics, through 2/ <dig> , 2/ <dig> , 3/ <dig> , with variation in chamaeleo spanning 2/ <dig> , 1/ <dig> , and 2/ <dig> .
all but brookesia and palleon in this study retain a tarsal sesamoid.
as in other lepidosaurs  <cit> , the foot skeleton is always composed of fewer elements than the manus.
in the hindlimb mesopodium, we observed 1– <dig> elements which appeared to be ectopic or de novo structures in bradypodion, chamaeleo, furcifer and trioceros.
knowing that in c. calyptratus the proximal tarsal is initially composed of condensations for the tibiale, intermedium and fibulare , we examined cartilage development in a  <dig> dpo c. calyptratus embryo  to follow the possible origin of these ectopic elements.
at this stage, the distal epiphysis of the tibia has grown extensively in the distal direction and displaced the cartilages of the intermedium and the tibiale toward the metapodial skeleton.
in adults , the tibiale, intermedium and fibulare fuse, forming a projection medially on the proximal tarsal element.
the tibia is able to redirect the placement of the tibiale and intermedium from their proximal position.
using these results as a model, the distal mesopodial element in bradypodion would be the tibiale and the proximal tarsal element would be comprised of a fused intermedium and fibulare.
across the genus chamaeleo, 0– <dig> ectopic elements are present distally while  <dig> are visible in both furcifer and trioceros.
using the distal epiphyseal growth of the tibia as a model for proximal tarsal elements as seen in c. calyptratus, it is evident that the distal elements in the former two genera are due to a more distal tibiale and a more proximal intermedium, with the proximal tarsal element composed of only the fibulare.
thus, distal displacement of the intermedium and tibiale is due to an accelerated growth heterochrony in the epiphysis as well as a loss of fusion for the three proximal tarsal elements which generally produce a single proximal tarsal.
the proximal tarsal element in different chameleon taxa is also shown to have different compositions and are thus not homologous in our study.
functional and phylogenetic constraints on mesopodial skeletal elements
carpal and tarsal tunnels of the wrist and ankle, respectively, guide tendons, nerves and blood vessels from the zeugopod into the autopodial domain and are involved in flexion of the more distal skeletal elements.
chameleons, unlike other tetrapods, have an autopodial skeleton with an approximately  <dig> degree perpendicular fixed angle to the proximodistal axis of the appendicular skeleton .
in association with a highly modified ‘ball-and-socket’ architecture of mesopodial elements, tissues extending through the tunnels must not only be enclosed through retinacula but must also be inhibited from sliding laterally and out of the tunnel by digits  <dig> and  <dig> during wrist extension, the position of which is fixed in chameleons  <cit> .
to compensate for the perpendicular orientation of the autopodium in chameleons, wrist elements of the hand  have enlarged and grown ventromedially to laterally border the tunnel connective tissues.
furthermore, the large distal carpal element  develops a ventral concavity to funnel the connective tissues to the palmar side toward the ventral aspect of the metacarpals .
in lepidosaurs, metatarsal  <dig> and distal tarsal  <dig> are fused  <cit>  .
thus, a hooked ‘5th metatarsal’ serves as a synapomorphy for this clade.
in the chameleon foot, ‘mt5’ and mt <dig> bilaterally border the connective tissues through the tarsal tunnel to the plantar surface of the foot.
basal species exhibit a single proximal tarsal element and a large distal tarsal element.
phylogenetically derived genera  enclose the tarsal tunnel connective tissues by increasing the relative size of the fibularis brevis tendon while also making the tibiale and intermedium independent from the proximal tarsal element.
this changes the composition of the proximal tarsal element across taxa in exchange for increased ventromedial support of the tarsal tunnel connective tissues .
discussion
chameleon autopodia have undergone extensive modification for optimized locomotion in a complex three-dimensional arboreal habitat.
here we have taken advantage of the slowly developing veiled chameleon   <cit>  to characterize the morphogenesis of the squamate limb and study various aspects of morphogenesis of the hands and feet.
this has aided our understanding of the diversification of autopodial morphologies across  <dig> of the  <dig> currently recognized chameleon genera and how the derived autopodial morphology of chameleons emerged within squamate reptiles.
chameleons exhibit a typical pentadactyl morphology, but also present with a midline autopodial cleft  and two opposable syndactylous bundles of highly mobile digits.
in addition, the proximal skeleton of the autopodium in chameleons was modified via a reduction in the number of bone elements.
due to the presence of  <dig> digits in chameleons, the ancestral condition for tetrapods, we expected a conserved pattern of expression for both fgf <dig> and shh in the developing autopodia of chameleons.
indeed, fgf <dig> was observed in the apical ectodermal ridge  while shh was present in the zone of polarizing activity  from the onset of limb bud outgrowth through termination of distal limb outgrowth and morphogenesis.
all previous studies of autopodial ectrodactyly and clefting have concluded this phenotype as having occurred due to a loss of integrity of the aer.
in contrast, we show that in chameleons, the aer remains a robust ectodermal thickening and expresses fgf <dig> during proximal mesenchymal clefting.
interdigital cell death, also an ancestral condition for tetrapods, was mainly observed in chameleon autopodia between the digit pairs undergoing clefting while adjacent digits retained interdigital mesenchyme .
interestingly, inhibition of bmp signaling prior to the initiation of interdigital cell death prevented cleft formation.
however, upon the commencement of cleft formation, blocking bmp signaling has no effect on interdigital cell death.
cell death between syndactylous digits was present in chameleons, but at significantly lower levels relative to the clefting mesenchyme.
the syndactyly observed in chameleons may differ to that observed in duck feet and bat wings due to the absence of grem <dig>  a bmp antagonist, in interdigitial tissue.
while previous studies credited interdigital cell death along the autopodial midline for the modified chameleon limb morphology, we highlight the previously unstudied and underappreciated skeletal diversity present in the mesopodial architecture across chamaeleonidae.
while retaining the ancestral pentadactyl morphology together with the typical  phalangeal formulae, the chameleon wrists have been highly modified to have a reduced number of skeletal elements with subsequent enlargement of those elements that are maintained.
proximal carpals and tarsals of digits  <dig> and  <dig>  respectively, of the hands and feet also make contact with each other.
we have shown that during embryogenesis in c. calyptratus, the carpals comprise  <dig> skeletal elements and the tarsals  <dig>  each of which originate as mesenchymal condensations, chondrify, ossify or fuse with adjacent elements at various stages of skeletogenesis.
we have used this model from c. calyptratus to homologize carpal and tarsal structures in  <dig> of the currently recognized  <dig> chameleon genera and found that chameleon research has significantly underestimated mesopodial skeletal diversity.
the earliest diverged genera palleon + brookesia, and rieppeleon are not only the smallest members of the family, but also show the greatest degree of skeletal number reduction in the wrist and ankle.
more phylogenetically derived chameleon genera from this study chamaeleo, trioceros, furcifer  show the greatest number of independent skeletal elements present in the wrist and ankle within the family.
this pattern also most closely resembles the skeletal complement found in outgroup lepidosaurs.
a phylogenetic trend of gigantism appears to potentially correlate with the ‘re-evolution’ of mesopodial skeletal elements, as this is associated with an increased number of cells from which to build mesenchymal condensations in the mesopodia.
an increase in mesopodial elements in ‘true chameleons’ also changed the resting posture of the hand and foot of chameleons, resembling a more typical lizard position where the two syndactylous bundles are oriented more anteriorly along the midline but are constrained due to the retainment of a cleft .
chameleons therefore provide an important example of how heterochrony can modify developmental patterns while being constrained by functional requirements and possibly body size.
thus, extrinsic and intrinsic factors are selected by nature to produce highly adaptive morphologies.
loss of carpal and tarsal elements increases the angle across the cleft zone, with the greatest angle being present in the most basal and smaller genera.
more derived chameleon genera have an independent intermedium in the ankle, a structure which has not been described in any extant terrestrial amniote since the origin of the astragalus and calcaneum from the amphibian skeletal complement.
in the wrist, segmentation takes place along the ventrolateral wall of the ulnare and gives rise to the pisiform bone.
in contrast to the sesamoid bone of the fibularis brevis tendon of the hindlimb, the pisiform appears to originate with the primary endochondral skeleton and not within a tendon.
the family chamaeleonidae is thought to have originated through the heterochronic process of paedomorphosis  <cit> , with adults resembling late gestation embryos of more distant taxa  <cit> .
malagasy species belonging to palleon and brookesia range in total length from 2–10 cm  <cit> .
here we have identified that early diverged and basal genera  exhibit a dramatic reduction of independent carpal and tarsal elements in adults despite maintaining a normal pentadactyl number of digits and phalanges.
this is in contrast to salamanders of equivalent size which exhibit digit loss and mesopodial element reduction .
while we only sampled  <dig> of the  <dig> extant genera , an ascending phylogenetic trend of increased body size appears to correlate with an increase  in mesopodial element number to a complement close to that of outgroup lepidosaurs.
more components in the wrist and ankle may have facilitated increased wrist flexion and provided a biomechanical advantage which has allowed more derived chameleons to leave the ground cover and low bushes for the trees.
clefting of the autopodium via interdigital cell death is the ancestral condition in the tetrapod limb.
in chameleons, interdigital tissue is however retained between fingers neighboring the midline cleft.
this represents the derived condition.
in phylogenetically more derived chameleon lineages with ‘re-evolved’ mesopodial skeletal elements, the ancestral tetrapod condition where all digits are present along the distal midline of the limb  is almost achieved.
while clefting separates the syndactylous bundles of digits, it is clear from the transiently bent chondrified digital rays of the alcian prepared skeletons  that there is either a musculoskeletal component driving the split, pulling the digit bundles toward the anteroposterior poles, or there may be a proliferative phase in the mesenchyme prior to cleft formation by apoptosis.
mesopodial elements have not received as much attention in the development literature when compared to studies of digit loss and phalangeal reduction.
despite the three major subdomains of the limb  being highly conserved across tetrapods, species-specific differences in the morphometric dimensions of long bone elements of the zeugopodia may elicit a direct effect on carpal and tarsal elements.
broader elements  would provide a wider  dimension of initial mesopodial mesenchymal pioneer cells for the more distal limb structures.
we hypothesize that altered proximal limb width, especially in the zeugopodial domain, may affect chameleon evolution as a narrowing of this domain would ultimately provide a reduction in the distal mesenchymal cells available for building the mesopodial and acropodial skeleton.
this model may also have been important for the evolution of digit and carpal reduced taxa such as hemiergis  <cit>  as they, like chameleons, have qualitatively thinner, elongated limbs with a smaller distal radius in zeugopodial skeletal elements.
larger taxa such as mammals, which lack a bottleneck in cell number for the mesopodial element and acropodial modification, employ genetic mechanisms for modifying adult phenotypes  <cit> .
in contrast, small ectotherms appear to exploit the cell number constraint model and typically exhibit element size and element number variation.
within this framework, a smaller mesopodial mesenchymal domain would be expected to give rise to a reduced number of carpal/tarsal elements as well as fewer cells able to contribute to digit  <dig>  notably, while two different studies  <cit>  have demonstrated a temporal reduction in shh expression in association with digit reduction, the tetradactyl pattern in these studies was different.
in the mouse, the digits were identified as having the character of digits i, ii, iv, v, whereas the digit reduced lizard retained digits ii, iii, iv, v together with a reduction in mesopodial elements.
digit i loss in hemiergis is consistent with our hypothesis that this phenotype is associated with a smaller initial pool of mesopodial mesenchymal cells.
fusion and reduction are more prevalent in organisms of smaller size, as was seen in this study, and this may simply be due to fewer mesenchymal cells being available for individual condensations to form or reach a sufficient density required for establishment of the chondrogenic pathway  <cit> .
mesopodial elements do not ‘re-evolve’ de novo
dollo’s law on the irreversibility of complex traits has recently been challenged in taxa where autopodial elements appear to have been re-evolved  <cit> .
phylogenetic trends in chameleon autopodial elements are not due to de novo ‘re-evolution’ of skeletal elements, but rather the elaboration of skeletal condensations already present.
in members of the genera chamaeleo, furcifer and trioceros sampled, mesopodial elements appear to be present but are underdeveloped or fused to adjacent condensations or skeletal elements, character states which can be easily missed when a detailed developmental series is not available for study.
mammalian species appear to reduce autopodial elements through genetically controlled cell death .
while a reduced period of shh expression in the zpa may also decrease the population size of mesenchymal cells available for autopodial development  <cit> , narrow limbs such as those in hemiergis lizards may have played an equally important role in reducing the cellular pool leading to element reduction  <cit> .
small mesenchymal population pools may also be a possible constraint on other members of the family scincidae such as brachymeles  <cit>  and larutia  <cit> .
more laboratory and field studies are needed to test this.
differential fusion between mesopodial condensations in chameleons blurs their skeletal homologies with outgroup lepidosaurs, especially if the topographic position of an element of interest or the number of ossification centers present are the only criteria considered.
herein we show that interspecific variability exists in not only the number of elements present  but also that different combinations of mesopodial elements fuse to give rise to the adult carpal and tarsal complement along with distal migration of other elements .
caution must be exercised when interpreting skeletal morphology using only dry skeletal specimens and microct imaging, as neither allows for observation of cartilaginous tissue  <cit> .
we have observed that fused condensations do not necessarily have multiple ossification centers, as individual condensations remain modular  and arrest their skeletogenic program at the mesenchymal, chondrogenic or ossification stage.fig.
<dig>
chamaeleo calyptratus atuopodial skeletogenic condensations.
despite a very specialized and element reduced adult chameleon autopodial skeleton, most ancestral skeletal elements remained present through ontogeny with the exception of the centralia.
between the manus and pes, the preaxial elements were always significantly delayed  in appearance.
the intermedium appeared to have a ‘dual’ origin, with its condensation occurring closer to the radius in the forelimb and the fibulare in the hindlimb.
we found that elements condensed and completed skeletogenesis .
there was a high tendency for fusion amongst proximal elements of the foot in c. calyptratus despite delayed appearance of the tibiale.
distal carpal  <dig> and distal tarsal  <dig> were seen to condense as mesenchyme but later fuse to adjacent elements while we were unable to find any condensation for distal tarsal  <dig>  distal carpal  <dig> remains cartilaginous.
the pisiform was found to be derived through segmentation of the ulnare and the fibularis brevis tendon to enclose a sesamoid
fig.
<dig> evolutionary trends in chameleon autopodial morphology.
outgroup comparison was conducted with the extant rhynchocephalian sphenodon punctatus  and  <dig> outgroups within squamata ], agamidae ).
fusion of proximal tarsal elements arose during amniote evolution by formation of the astragalus  and the calcaneum  .
aspidoscelis  and pogona  both show an increase in fusion of mesopodial elements relative to sphenodon .
the small and early diverged chameleons , palleon) show the distal clefting of the autopodia between the forelimb digits iii and iv and hindlimb ii and iii characteristic of chameleons.
these two genera also present the greatest modified proximal autopodial skeleton within chameleons.
the larger rieppeleon and rhampholeon ) presents a sesamoid in the fibularis brevis tendon of the heel which is maintained in all higher chameleons.
bradypodion , chameleo , furcifer  and trioceros   sampled in this paper showed an ascending phylogenetic trend toward increased body size and an associated increase in skeletal elements present in the mesopodium, leading to a wrist with increased flexion capability as well as being the first amniote group to form an independently ossified tibiale and intermedium which is only known to be present in amphibians.
thus, in early chameleon lineages there appears to be a skeletal bottleneck relative to outgroup squamates with a subsequent phylogenetic trend of ‘re-evolving’ mesopodial elements present in outgroup taxa.
the red gradient shows a trend toward reduction of elements from outgroup taxa to chameleons while the red gradient shows a relative increase in mesopodial skeletal elements in advanced chameleon genera studied.
genera in “<genus>” are currently recognized taxa but were not examined in this study.
letters within figure correspond to the species with the same letters in additional file 1: figure s1
fig.
<dig> phylogenetic reversal to a ‘typical’ lizard conformation of the autopodia.
despite retaining the cleft, an increase in mesopodial elements leads to increase wrist flexion in higher chameleons.
concomitant with the increase in the number of carpal and tarsal elements, we see a bifurcation where the syndactylous bundles are almost 180° separate in the early diverged lineages , rieppeleon , and rhampholeon ) with the angle decreasing in bradypodion .
‘true chameleon’ syndactylous bundles shift toward the distal midline as is present in the basal squamate plan .
while chamaeleo  show an anteromedial transition of digit bundles, furcifer  and trioceros , the taxa with greater number of carpal and tarsal elements in this study, exemplify this the best through a greatly reduced angle between syndactylous bundles as digits move closer together along the distal midline.
fl: refers to forelimb with hl: referring to hindlimb.
the numbers associated with each lineage on the phylogeny correlate to the same order of illustrations from fig.
<dig> and represent formulae for number of elements present in the mesopodium in association with a representative image for autopodial morphology of a preserved specimen.
fl: x/x refers to the number of elements in the two mesopodial rows while fl: x/x refers to the two mesopodial rows with the number in parentheses representing the presence or absence of the sesamoid of the fibularis brevis.
no centralia are present in chameleons, which would have made the standard  <dig> rows of mesopodial elements.
genera in “<>” are currently recognized taxa, but were specimens were not obtained in this study
homology and evolution of the astragalus and tarsal bones
a proximal reduction in the number of carpal and tarsal elements in adult chameleons, relative to embryos and outgroup taxa, together with distal clefting, have allowed for evolution of the novel chameleon hand and foot which are used for complex grasping in an arboreal three dimensional environment.
the progressive addition  of skeletal elements in higher taxa may have facilitated increased flexion and a greater ability to grasp, which correlates with chameleons invading the arboreal niche  <cit>  where only mammals with autopodial modifications such as marsupials have evolved to inhabit  <cit> .
chameleons are thought to have invaded madagascar through two independent events  <cit> .
the initial invasion by the ancestor of brookesia and palleon evolved in situ, and was followed by a secondary invasion of the ancestor of calumma and furcifer by dispersal from africa.
calumma + furcifer  are the sister lineage to trioceros + kinyongia , which collectively comprise the four most phylogenetically derived extant chameleon genera  <cit> .
from our studies of the mesopodial architecture of these african and malagasy clades, we can infer that the dispersed ancestor of calumma + furcifer to madagascar was already mesopodially ‘pre-adapted’ to be an advanced climber due the larger complement of skeletal elements in the wrist and ankle.
brookesia and palleon were already in madagascar as forest floor dwellers.
given the advanced nature of the wrists in calumma, furcifer, trioceros and kinyongia, such increased mesopodial flexion has allowed them to comprise ~57 % of the  <dig> extant chameleon species  <cit> .
innovations associated with the evolution of terrestrial cursorial ability in amniotes involved modifications to the ankle.
evolution of the astragalus has elicited much discussion with recent models attempting to ascribe it an accurate homology with respect to tarsal skeletal elements found in anamniotes, generally involving 1– <dig> ossification centers .
our work reveals a high degree of variability in the composition of the chameleon and squamate astragalus prior to the generally accepted co-ossification step following mesenchymal condensation.
chameleons of the genera chamaeleo, furcifer and trioceros as shown here, are currently the only known amniotes to have successfully ‘re-evolved’ an independent intermedium and tibiale.
a free intermedium and tibiale was lost in amniotes during the tetrapod progression onto land which selected for innovation of the astragalus and calcaneum  <cit> .
reptiles are considered to not have a tibiale .
this ‘re-evolution’ of an ancestral series of tarsal elements may be due to a loss in selective pressure toward the typical terrestrial mode of tetrapod locomotion, in combination with a shift toward arboreality.
within this framework, different amniote lineages should be expected to develop an astragalus through the fusion of tarsal elements at different points through ontogeny.
although each element is modular, the distinct stage of skeletogenesis prior to and after fusion may vary across squamata .
in summary, examining the development of innovations and novelties in natural groups such as chameleons can further our understanding of how body plans are molded by natural selection as optimal phenotypes through adaptive evolution.
CONCLUSIONS
studying organisms in the lab with unique morphologies complements biomedical studies of normal development and congenital malformations while also providing insight into the limits  and novel deployment of molecular genetic pathways during embryonic development.
veiled chameleons  present a highly derived body plan within tetrapods due to their specialization for living in an arboreal environment while also presenting a very slow rate of embryonic growth .
developmental analysis of the specialized and bifurcate autopodia of c. calyptratus using histological, in situ hybridization and in vivo techniques show that the specialized mesopodial skeleton  forms despite the expression of highly conserved patterning genes during limb outgrowth.
our detailed skeletal developmental sequence present an unexpected doubling in the number of recognized wrist and ankle elements which has allowed us to homologize mesopodial elements across  <dig> of the  <dig> recognized chameleon genera.
our data also present a correlation between body size and skeletal diversity of the wrist and ankle during chameleon evolution with phylogenetically higher taxa exhibiting an increase in mesopodial skeletal element number, larger body size, and increased climbing ability.
exploring development in a slow developing taxon has also allowed us to observe the fine progressive development and segmentation of the pisiform from the ulnare in the hand and differential fusion of proximal tarsal elements in the ankle.
additional files
additional file 1:
figure s <dig>  comparative morphology of autopodia in this study.
dorsoventrally flattened autopodia in outgroup taxa  were drawn in dorsal view only  for both forelimb and hindlimb autopodia.
d-n are chameleon illustrations from taxa cleared and stained in this study to look for ossified and chondrified elements comprising the autopodial architecture and are illustrated in both anterior/dorsal and posterior/ventral views.
a) sphenodon punctatus , b) aspidoscelis uniparens , c) pogona vitticeps , d) brookesia stumpffi, e) rieppeleon brevicaudatus, f) rhampholeon boulengeri, g) bradypodion pumilum, h) chamaeleo dilepis, i) chamaeleo calyptratus, j) chamaeleo chameleon, k) furcifer lateralis, l) furcifer pardalis, m) trioceros ellioti, n) trioceros hohnelii.
additional file 2:
figure s <dig>  forelimb condensations of a. uniparens in h&e stained sections.
a large contribution of preaxial mesenchymal cells contributing to digit 1/distal carpal  <dig> are visible in these sections as well as the development of an independent lateral centrale and intermedium at an older stage.
the medial centrale will later appear as a medial segmentation from the distal carpal  <dig>
additional file 3:
figure s <dig>  pisiform versus sesamoid identity in the chameleon autopodium.
the top row  shows the forelimb with the pisiform bone which we have identified in this study as segmenting from the ulnare.
this bony element is not present within a tendon and instead has a proximal insertion of the tendon of the flexor carpi ulnaris with a more distal pisometacarpal ligament.
the lower row  shows the hindlimb with a large sesamoid present within the tendon of the fibularis brevis which inserts onto the proximal margin of the ‘5th metatarsal.’ the species used for this image was trioceros ellioti .
additional file 4:
figure s <dig>  whole embryo culture of c. calyptratus.
details for this protocol are provided in the article text in the methods section and correlate with labeled figures.
additional file 5:
figure s <dig>  pisiform formation in the manus of c. calyptratus.
pisiform development from the ulnare is shown in posteroventral , lateral , and dorsal view  in whole alizarin/alcian skeletal preparations.
transverse h&e stained paraffin sections show, from distal to proximal mesopodium, the shared mesenchymal cells between the ulnare and the pisiform, supporting the ventrolateral segmentation of the pisiform from the ulnare .
additional file 6:
figure s <dig>  autopodia are dorsoventrally flattened divergent tetrapods such as the common lab mouse  and the desert grassland whiptail lizard  have divergent hands  and feet .
despite architectural differences, they remain dorsoventrally flattened.
additional file 7:
chameleon specimen measurements.
abbreviations
gremgremlin
bmpbone morphogenetic protein
shhsonic hedgehog
fgffibroblast growth factor
icdinterdigital cell death
aerapical ectodermal ridge
zpazone of polarizing activity
dpodays post oviposition
pteproximal tarsal element
svlsnout-vent length
tltail length
flforelimb
hlhindlimb
roman numerals i-vrepresent digit 1- digit 5
raradiale
rradius
uulna
ululnare
ppisiform
iintermedium
ttibia
titibiale
ffibula
fifibulare
dcdistal carpal
dtdistal tarsal
mcmetacarpal
mtmetatarsal
lacelateral centrale
mecemedial centrale
competing interests
the authors declare that they have no competing interests.
authors’ contributions
red and pat conceived and designed the study.
red carried out the laboratory work and data collection.
red and pat drafted, read and approved the final manuscript.

 understanding the complex relationship between genotype and phenotype requires an integrative and interdisciplinary biological framework. loci known to be involved in development, morphogenesis and in the pathogenesis of congenital malformations have been identified through forward genetic approaches, genomic mapping or genome wide association studies. concurrent with studying how malformations arise, cellular and genetic mechanisms have been uncovered which have direct consequences on our understanding of the processes governing normal development. examining the development of innovations and novelties in natural groups provides an understanding of how body plans have been molded by natural selection as optimal phenotypes through adaptive evolution. limbs and digits were key innovations in the evolution and diversification of tetrapods. while much has been learned about gene function during morphogenesis and differentiation of the tetrapod limb, most studies examining non-traditional model organisms have focused on limb reduction or loss of distal elements of the autopodium. within reptilia, the order squamata comprises more than lineages representing independent limb reduction toward a snake-like body form. however, one of the most dramatic limb skeletal modifications is found in the hands and feet of the family chamaeleonidae, the chameleons. chamaeleonidae limbs have been modified from the generalized terrestrial tetrapod plan to exhibit an architecture highly adapted for an arboreal lifestyle. this was facilitated through evolution of a midline autopodial cleft and two opposable syndactylous bundles of digits that are highly mobile. in addition, the proximal skeleton of the autopodium was modified through a reduction in the number of bone elements. furthermore, those elements that remain are enlarged and form a ball-and-socket joint between the autopodium and forearm. this specialized joint allows for greater rotation of the wrist and ankle, which is important while climbing. limb modifications such as these facilitated an increase in precision, security and mobility in the three-dimensional chameleon arboreal environment. while other lizards are also arboreal, chameleons differ by using zygodactyly and gripping of branches rather than depending on claws and epidermal modifications.fig. chameleon autopodia. a the veiled chameleon has a laterally compressed body, prehensile tail, turreted and independent eyes, a projectile tongue and zygodactylous cleft autopdia. b chameleon hand digits i, ii, iii are syndactylous medially with digits iv, v are bound laterally. feet have digits i, ii syndactylous medially while digits iii, iv, v are bound laterally. d the desert grassland whiptail lizard was used as for comparison due to its ‘typical’ terrestrial lizard body plan few studies to date have focused on the development of chameleon autopodial architecture, leaving unanswered the question of how the highly conserved tetrapod pentadactyl autopodium was modified during chameleon evolution. to address this issue, we have developed and utilized the veiled chameleon, which offers key advantages over other lizard and chameleon species. one of the most important and distinctive developmental traits in c. calyptratus is their slow rate of embryogenesis, which facilitates the examination of fine, progressive, morphological changes. chameleons provide an important example of how developmental mechanisms may be modified through heterochrony and constrained by body size to produce highly adaptive morphologies. with respect to the limb, a proximal reduction in the number of carpal and tarsal elements in adult chameleons with concurrent distal cleft formation of autopodia have facilitated the evolution of complex zygodactylous grasping in a three dimensional arboreal environment. we collected c. calyptratus embryos at various time points to study gene expression, morphogenesis, and skeletal and histological development of autopodial elements. we also examined of the currently recognized genera of chameleons and outgroup taxa to understand the evolution and diversity of carpal and tarsal elements across this phenotypically derived clade of arboreal lizards. a reversal to an increased number of carpal and tarsal elements together with increased independence of elements, appears to have played a significant role in the evolution of ‘true chameleons’. such a mechanism could provide more components by which to increase wrist flexion, a biomechanical advantage in a complex arboreal environment. methods embryos veiled chameleon eggs were collected at oviposition in the stowers institute for medical research’s reptiles and aquatics facility and at the la sierra university animal facility. eggs were subsequently incubated at a constant 26 °c in deli cups with moist vermiculite. embryos were harvested in cold 1× phosphate buffered saline at varying stages, fixed overnight in depc-treated 4 % paraformaldehyde in 1xpbs at 4 °c and later dehydrated to absolute methanol before storing at −20 °c. dapi helicon imaging whole embryos were collected in cold 1x phosphate buffered saline and fixed overnight at room temperature with 1 μl dapi, dilactate per 1 ml of pbs on a roller. embryos were placed flat in a 50 ml petri dish filled with pbs and imaged on a leica mz fliii with a carl zeiss color axiocam hrc using axiocam software. the uv filter was used to visualize dapi fluorescence. settings were changed to black and white in the software to generate a grayscale image of the embryo. z-stack images were taken and subsequently processed with helicon focus software, which functions to create a projection of all images into a single high focus photo. this process is analogous to doing a dapi z-stack on a confocal, and is useful for generating high resolution images of embryos too large for visualization on a confocal microscope. this method also overcomes some of the limitations of sem imaging in that the embryos never have to be dehydrated, can be kept at room temperature in fixative indefinitely, can be used on museum specimens, is non-destructive, and also provides simultaneous tissue resolution below the differentiating skin. a 50 ml plastic petri dish was half-filled with 1 % agarose and left to solidify. troughs in the agarose were created to accommodate embryos in various orientations and provide stability during imaging. rna collection, cdna synthesis, probe development and in situ hybridization embryos between and dpo were collected in cold depc-treated 1xpbs. limbs were pooled together and total rna was extracted with the rneasy mini kit. complementary dna was derived from messenger rna using the superscript® first-strand synthesis system for rt-pcr. degenerate primers were designed for genes of interest. gotaq® colorless mastermix was used for pcr amplification. reaction products were run on a 1 % agarose gel and bands of expected sizes were excised, purified using the qiaex ii gel extraction kit, and stored at −20 °c. amplified fragments were ligated into pgem®-t easy vectors. dna was sequenced using sp and t primers, and checked via blast to determine the direction of the insert. plasmids were digested using the appropriate restriction enzymes for the development of sense and antisense in situ hybridization probes using the dig rna labeling kit. degenerate primers used include: for sonic hedgehog shhf: 5’-ccggcttcgactgggtntayta-3’, shhr: 5’-catgggcggtcagtggngcrtangc-3’. for fibroblast growth factor fgf8f: 5’-ctmrtscgsacctaccarct-3’, fgf8r: 5’-ggsarsckcttcatgaagtg-3’. for gremlin grem1f: 5’-caaggmtcagcmcaatgaytc-3’, grem1r: 5’-yttcttgggcttgcagaag-3’. genbank accession numbers for the above amplified fragments are. for in situ hybridization, embryos were collected in cold depc-treated 1x pbs, and fixed in 4 % paraformaldehyde overnight at 4 °c on a rocker. specimens were then dehydrated through a graded series of methanol/pbs-depc into 100 % methanol and stored at −20 °c. whole mount in situ hybridization was performed as previously described. micro-computed tomography comparative limb osteological data was obtained by isolating limbs from adult specimens followed by scanning with a scanco vivact micro-computed tomography scanner. dicom data was visualized using amira. scanning electron microscopy embryos were fixed overnight in  % gluteraldehyde/2 % paraformaldehyde in 1xpbs, rinsed in water, incubated in 2 % osmium tetroxide, rinsed again in water, and dehydrated through a graded series of ethanol for 20 min each. critical point drying was conducted in tousimis samdri- critical point dryer. scanning was accomplished on a hitachi tm- tabletop microscope. alcian/alizarin skeletal preparations staining of bone with alizarin red and cartilage with alcian blue 8gx was adapted from walker and kimmel. initially, embryos were fixed in absolute ethanol and dehydrated/stored in absolute etoh. specimens were stained overnight in alizarin red/alcian blue solution at room temperature with rocking. on the following day, specimens were rinsed in water, and, depending on the size of the specimen, were incubated in 4 %, 2 %, 1 % or  % koh until sufficient soft tissue digestion had occurred. specimens were then cleared in 25 % glycerol in  % koh, followed by 50 % glycerol in  % koh before being post-fixed in 2 % formaldehyde in 50 % glycerol/1xpbs, which is also the permanent storage media for specimens. adult museum specimen catalogue numbers are included in the images in additional file 1: figure s with each species represented by its unique tiff file labeled a-n. specimen measurements can be found in additional file h&e paraffin histology embryos were fixed in 4 % pfa overnight, dehydrated through an ethanol series, cleared in xylene, and embedded in paraffin. 8–10 μm sections cut onto slides were deparaffinized with xylene, re-hydrated through an ethanol series into water and counterstained with hematoxylin and eosin. whole embryo culture eggs were incubated at 26 °c for 100 days and candled to check for viable embryos. healthy eggs limit light passing through the egg and also exhibit membrane vasculature, whereas infertile or unhealthy eggs appear quite clear. pin-hole punctures with forcep tips allowed for very clear and viscous albumin/chorioallantoic fluid to leak onto the egg shell, while infertile eggs burst, releasing a creamy white yolk. shell opening ruptures the chorioallantoic membrane, however the embryo remains within the amnion embedded within the yolk sac. the amnion can be seen as a remnant of membranes externally. the supine position of the embryo is associated with the constraint imposed by the umbilical arteries and veins leading toward this thin transparent membrane where the yolk sac and allantoic membranes both diverge. care must be taken to cut through the clear membrane without severing the umbilical veins/arteries. the egg yolk was cut to form a ‘blanket’ below the embryo while retaining the remaining chorioallantoic membrane and yolk sac vasculature. embryos were then placed in organ culture media within a 6-well plate. embryos were cultured in an exo-terra thermoelectric reptile incubator for 4 days at 26 °c under high humidity. media was changed daily. media for culture culture media consisted of  ml of fetal bovine serum,  ml dmem/f- 50/ with l-glutamine, 10 μl of 50 mg/ml ascorbic acid for final 50 μm/ml concentration. 200 μl of pen-strep was added to limit bacterial contamination. nile blue staining embryos were incubated in 1: dilution of a  % stock solution of nile blue in culture media for 40 min in a 26 °c incubator, followed by two 30 min washes in pbs with rocking at 4 °c prior to photographing. ldn treatment dorsomorphin was added to a final concentration of 15 μm in culture media. animal research all procedures were formally approved by the institutional animal care and use committee at the stowers institute for medical research and the animal care and use committee at la sierra university. no wild caught animals were used in this study. all embryos studied were produced in house.  interdigital clefting is independent of a robust aer to begin to understand the mechanisms governing the modified hand and foot morphology of chameleons, we compared autopodial development in the veiled chameleon to the desert grassland whiptail lizard, which serves as a reference species due to its overall ‘typical’ lizard limb and body plan. in lizards, it takes approximately 10 days for limb development to progress from a round digital plate to the typical pentadactyl pattern which coincides with the onset of interdigital cell death. in contrast, autopodial morphogenesis in c. calyptratus lasts approximately 14– days from the rounded paddle stage through cleft completion. phenotypic divergence between these species commences at the late round digital plate stage when distal flattening occurs in the chameleon forelimb and hindlimb. subsequent anteroposterior expansion of the flattened region delimits the clefting domain such that forelimb and hindlimb digits differ in their syndactylous clusters, with digits grouped as and, respectively.fig. comparative morphogenesis of lizard autopodia. in aspidoscelis uniparens, the digital plate progresses through a round shape with the pentadactyl digital rays developing within and visible through the thin tissue. subsequent thickening of digital rays and thinning of the interdigital mesenchyme ultimately leads to loss of the interdigital tissue through apoptosis. in the chameleon hand, a semi round digital plate develops with subsequent distal flattening. clefting initiates at this distally flattened location. two syndactylous digit bundles are approximately 180° apart in the anterior-posterior plane. the a. uniparens foot also passes through a semi round digital plate stage which subsequently presents digit differentiation with subsequent interdigital cell death. the chameleon foot undergoes the same stages as in the hand, but terminates with a different complement of digits in syndactyly. dpo = days post oviposition; a = anterior; p = posterior current knowledge of cleft formation in autopodia comes from examination of avian and mammalian systems as well as linkage studies in humans, with all work supporting a failure to maintain the integrity of the apical ectodermal ridge, particularly along the distal midline. loss of the aer inhibits not only distal outgrowth of the limb but also leads to the loss or splitting of the digital rays. to determine whether or not aer destabilization is the primary factor driving cleft formation, we examined the aer of a. uniparens in order to first understand how a ‘typical’ aer is structured in pentadactyl lizard autopodia. at dpo, a robust distal ectodermal thickening is present at the dorso-ventral boundary of the forelimbs and hindlimbs. at dpo the aer begins to narrow as the digital rays differentiate within the underlying mesenchyme. by 20– dpo the aer further narrows prior to the initiation of interdigital cell death. by comparison, in chameleon autopodia, the aer is thicker and more prominent. interestingly, upon distal flattening of the digital plate, commencement of distal clefting in chameleon interdigital tissue is visible despite the retainment of a very robust aer. syndactylous digit clusters diverge toward the anterior and posterior poles, increasing the angle between clusters and further expanding the cleft. the proximodistal size of the autopodium remains unchanged during these stages. overall, autopodial shape shifts from a distally flattened rounded paddle to one that is rectangular shaped. it is interesting to note that as the aer is maintained, it appears slightly thicker in the region overlying the clefting mesenchyme, which may be compensatory to maintain interdigital tissue proliferation during this midline mesenchymal expansion in concert with digit bundle realignment prior to cleft formation.fig. sem of aspidoscelis uniparens hand and foot morphogenesis. aspidoscelis uniparens hands and feet undergo the same stages of morphogenesis as presented in fig.  in dorsal view, respectively. in distal view, the hand shows an apical ectodermal ridge which is initially very robust as a distal ectodermal thickening and is present at the midline between the dorsal and ventral half of the developing limb. at later stages, the aer narrows significantly as it spans a greater anterior-posterior domain. in the foot, we see a similar situation with a very straight and stereotypical aer which is initially robust and later thins fig. sem of chamaeleo calyptratus hand and foot morphogenesis. early stages of morphogenesis show that the chameleon hand develops as a round digital plate which subsequently develops a distal flattening. in dorsal view, these stages present a very robust aer, relatively larger than that present in a. uniparens while in distal view the aer is seen as having a greater thickness along the dorsoventral midline. significantly, during the stage at which the distal autopodium begins to flatten, the distal aer is no longer the stereotypical a-p flattened ectodermal thickening but is instead arched ventrally. surprisingly, despite having a robust distal aer, proximal mesenchymal cleft formation has already begun. at later stages, both the forelimb and hindlimb expand the cleft while maintaining the aer quite robust until the thickness tapers and narrows at later stages autopodial patterning is normal despite a divergent skeleton chameleons retain the pentadactyl tetrapod autopodial plan but differ significantly from other tetrapods due to the presence of ectrodactyly and a highly modified mesopodial skeleton. to establish whether early changes in limb patterning influence subsequent development of the novel hand and foot morphology in this group, we characterized the spatiotemporal activity of the major signaling pathways involved in limb development as previously determined in murine and avian species. we initially examined the spatiotemporal activity of sonic hedgehog. shh is known to play a central role in regulating anteroposterior patterning of the digits as well as proliferation of the autopodial mesenchyme which gives rise to the digits. furthermore, shh plays a feedback role with grem in the limb mesenchyme and with fgf in the aer to maintain distal limb outgrowth. to date, shh expression has only been examined during squamate limb development in lizard limbs, and it was shown that the temporal activity of shh expression in the posterior zone of polarizing activity plays a role in determining the number of skeletal elements in the digits of element-reduced species in the scincid genus hemiergis. in this case, shorter expression periods were observed in association with phalangeal and digit reduced taxa. in chameleons, the tetrapod pentadactyl plan is conserved. the phalangeal formulae for the hands and feet of c. calyptratus are only slightly modified. therefore, we hypothesized that the spatiotemporal activity of shh would also be conserved in chameleons. shh is in fact expressed in the expected posterior domain demarcating the zpa in veiled chameleons, as it is in the pentadactyl hemiergis. furthermore, shh is expressed in c. calyptratus throughout early limb bud development while the digital plate is rounded, but is subsequently downregulated upon distal flattening of the digital plates.fig. c. calyptratus autopodial sonic hedgehog expression. wholemount in situ hybridization indicates that shh is expressed in a similar pattern to other amniotes along the posterior margin of the developing limb within the zone of polarizing activity in the hands and feet. significantly, the expression of shh is lost upon distal flattening of the limb bud in both the hands and feet the aer forms at the dorsoventral boundary along the distal edge of the limb bud and in mice is the site of expression for members of the fibroblast growth factor family. in particular, fgf is necessary and sufficient for distal outgrowth and patterning of the growing limb. temporal removal of the aer at various developmental stages allows for only more proximal elements to form. using fgf as a candidate marker for aer integrity, we find that the aer maintains normal fgf expression along the anteroposterior distal boundary across all digits as would be expected for a pentadactyl limb. however, fgf continues to be expressed in the aer covering the distal boundary of the expanding interdigital domain which will form a cleft during anteroposterior movements of the two syndactylous digit bundles. ultimately, fgf expression diminishes in association with interdigital mesenchymal regression.fig. c. calyptratus fibroblast growth factor expression in the apical ectodermal ridge. in wholemount in situ hybridization, the apical ectodermal ridge expresses fgf in the highly conserved pattern as seen in tetrapods along the distal anteroposterior boundary of the limb. fgf expression is maintained in the aer at stages where mesenchymal clefting of the autopodia has already commenced in both the hands and feet, though at a reduced amount prior to complete cleft formation clefting is due to increased cell death recent work exploring digit number reduction supports an important role for cell survival in sculpting adult autopodial elements. the presence of macrophages in interdigital cleft mesenchyme has led to the suggestion that apoptosis is a primary driver of the derived chameleon autopodium. this idea warrants further validation, as mature digit formation generally requires apoptotic removal of interdigital mesenchyme in all known tetrapod limbs. macrophage associated or non-associated cell death would therefore be expected to be localized to the clefting mesenchyme in most tetrapods, with the exception of taxa with interdigital webbing where cell death is inhibited. we used nile blue staining as evidence for cell death in whole embryo cultured chameleons at stages during clefting. while the presence of macrophages is common in the aer to aid its maintenance and extension, we also observed increased cell death in the mesenchyme in the distal margin of the autopod. at the onset of interdigital cell death, we observed differential nile blue staining between the major clefting domain along the distal midline relative to the interdigital mesenchyme of digits bound in syndactyly. this reflects the comparative retention of mesenchyme between digits, compared to the cleft. at later stages of embryogenesis, the aer covering the cleft is no longer stained by nile blue. in contrast, anterior and posterior syndactylous domains retain nile blue staining in the aer of the hand and foot, highlighting the requirement for an aer in mesenchymal tissue maintenance. in ducks and bats, gremlin, an inhibitor of bmp signaling, is expressed in the interdigital mesenchyme promoting the retention of mesenchymal tissue, resulting in webbing. chameleons initially exhibit expression of grem in the developing limb bud mesenchyme anterior to shh in the zpa, which is typical for limbs during outgrowth. however, grem is subsequently downregulated in interdigital regions where syndactyly occurs.fig. nile blue staining for cell death during autopodial morphogenesis. during autopodial morphogenesis, nile blue granules are concentrated along the distal periphery of the developing limb within the mesenchyme as well as delimiting the apical ectodermal ridge in both the hands and feet. the aer is also highlighted as ventrally arched in distal view and is also clearly dorsoventrally thinner at stages when cleft formation is underway. during final stages of cleft morphogenesis, nile blue staining is concentrated within the interdigital mesenchyme, with a significantly greater concentration in the interdigital mesenchyme between digits iii and iv in the hand and ii and iii in the foot. reduced nile blue is present in the tissue within syndactylous digits. nile blue is visible in the aer on the anterior and posterior regions of the limb while absent in the cleft domain in the hand fig. c. calyptratus gremlin expression during autopodial morphogenesis. grem is expressed in the mesenchyme of the autopodium anterior to the zpa where a feedback loop is known to be present between grem and shh, bmps, and fgf. grem appears to decrease its anterior expression during digital plate formation in the forelimb and hindlimb. during distal flattening and cleft formation, grem does not appear ectopically between syndactylous digits members of the bone morphogenetic protein family play diverse roles in limb morphogenesis. during early limb development bmps are initially required for formation of the aer, but are then negatively regulated by grem in the limb mesenchyme to allow for activation of fgfs present in the aer. later during autopodial morphogenesis, bmp and are expressed within interdigital mesenchyme prior to and during the process of programmed cell death where soft tissue is removed to delineate each digit. inhibiting bmp signaling within interdigital mesenchyme through ectopic expression of bmp inhibitors such as noggin or gremlin or by dominant negative bmp receptors leads to a loss of programmed cell death in chick and mouse. to functionally test the role of active bmp signaling in establishing the apoptotic program in association with cleft formation in the chameleon limb, we cultured whole embryos in vitro in the presence of ldn a pharmacological molecule that inhibits bmp signaling by blocking smad activity. two groups of embryos were cultured for 4 days. one group termed ‘younger’ were collected at the distally flattened limb paddle stage of development prior to the onset of interdigital cell death. the ‘older’ group was collected after clefting had commenced. for ‘older’ embryos cultured at the onset of clefting, tissue loss continued despite exposure to dorsomorphin. however, if dorsomorphin exposure occured prior to the initiation of cell death, clefting was prevented. however, digital ray development was still slightly impaired, and the digits also appeared smaller relative to controls. smaller digital rays in bmp inhibited mouse limbs are due to direct inhibition of sox by fgf, and a delay in chondrification. blocking bmp activity maintains digital ray cells in an undifferentiated state due to continued fgf expression in the aer. treating mouse embryonic limbs with dorsomorphin during the period of interdigital cell death led to a reduction of smad phosphorylation in interdigital mesenchyme but did not affect apoptosis. however, treating mouse embryos with dorsomorphin at earlier stages prior to initiation of apoptosis, did inhibit interdigital cell death. previous studies genetically perturbing bmp signaling in chicken and mouse embryos has also led to comparable limb phenotypes. thus, our study of cleft formation in chameleon embryos reveals that blocking bmp signaling can inhibit interdigital cell death in a time dependent fashion prior to the onset of programmed cell death.fig. disruption of bmp signaling during autopod development. dorsomorphin at 15 μm was used to disrupt bmp signaling during whole embryo in vitro culture for 4 days. in the ‘younger’ group of embryos collected prior to distal flattening of the digital plate, treatment led to impaired cleft formation. asterisks highlight reductions in cleft formation relative to controls. b older embryos collected during the onset of clefting had no effect from the ldn treatment over a 4 day growth period cryptic diversity of mesopodial skeletal element number in c. calyptratus embryos limbs and fins come in a variety of shapes and sizes, with variable numbers of skeletal elements. despite divergent morphologies, all autopodia are dorsoventrally flattened structures, as can be seen by juxtaposing the hands and feet of the evolutionarily distant and morphologically divergent lizard and mouse. in instances where clefting occurs in mammals, the autopodium is bifurcated along the distal midline and results in missing elements or a reconfigured architecture.fig. tetrapod autopodia are dorsoventrally flattened, except in chameleons. chameleon autopodia while retaining the pentadactyl complement of digits, not only have a distal cleft separating the digits iii and iv in the hand and ii and iii in the foot but also develop a rosette shaped cluster of metacarpals and metatarsals around a central enlarged distal carpal element. in species with congenital distal clefting, clefting does not alter the dorsoventrally flattened morphology despite having a relatively normal complement of phalanges and a pentadactyl number of digits, chameleons have an extremely modified carpus and tarsus. in addition, the proximal skeleton of the autopodium is reduced in terms of the number of bone elements. furthermore, those that remain are enlarged and form a ball-and-socket joint between the autopodium and forelimb. architecturally, clefting and syndactyly are necessary for grasping, but these two phenotypic novelties are not sufficient to explain the ball-and-socket mesopodial skeleton. thus, wrist and ankle bone homologies remain contentious relative to outgroup lizards and other tetrapods. previous investigations of the potential homology of derived mesopodial skeletal elements in chameleons relied on late gestation stage trioceros hoehnelii and bradypodion pumilum embryos. similarly, comparisons to outgroup lepidosaurs also relied on late gestation stage embryos and were also based solely on the number of visible ossification centers in the mesopodia. consequently, distal carpals 2– were interpreted to have fused to make a large central element with the proximal forelimb elements comprising the radiale, ulnare and pisiform, while in the hindlimb only distal tarsal and a ‘tarsale proximale’ were present. to better understand autopodial evolution and structural diversity across chameleons, we generated a skeletal ontogenetic series of chamaeleo calyptratus in order to: 1) identify the developmental sequence of appearance of autopodial skeletal elements, 2) homologize the derived and simplified elements of the chameleon autopodium, and 3) understand what developmental events occur in the formation of c. calyptratus autopodia. in our skeletal analysis of chamaeleo calyptratus, we identified rows of skeletogenesis during development of the carpals: proximal carpal elements and distal elements. tarsal elements also developed from two rows: proximal elements and distal elements. interestingly, the intermedium appears to condense in association with different mesopodial elements, from the distal medial radius in the hand and from the distal medial fibula in the foot. thus, c. calyptratus retains more mesopodial elements than previously acknowledged, with only the medial and lateral centrale of lepidosauria failing to develop additional file centralia are generally the last elements to form ontogenetically and are consistent with a heterochronic terminal truncation. a clear series of proximal condensations is not present in the outgroup taxa aspidoscelis uniparens and pogona vitticeps, as observed by the increased and more diffuse alcian blue staining present across this domain.fig. chondrogenic differentiation in chameleon carpus. during distal flattening and expansion, the digital rays are already present and are lightly stained by alcian blue. the first elements of the mesopodium to form as cartilage are the fi, dc and metacarpals 3– distal carpal and form next with later appearance of the intermedium, metacarpal and ultimately the radiale and pisiform. all embryos were at dpo from clutchmates fig. chondrogenic differentiation in chameleon tarsus. during distal flattening expansion, the first elements to appear are the ulnare and distal tarsal followed by distal tarsal and metatarsals 3– and the intermedium. metatarsal is followed along with an appearance of the tibiale. all embryos were at dpo from clutchmates fig. outgroup lizard species autopodial chondrification. desert grassland whiptail lizard and bearded dragon chondrification patterns show post-axial dominance in elements appearance along the metapterygial axis. a. uniparens hands show early formation of the ulnare, distal carpals and metacarpals 2– and a diffuse staining in the area of the radiale, lateral centrale and intermedium a. this empty domain later fills with a large medial and lateral centrale b. the foot develops from less elements and completes with a fused metatarsal and distal tarsal and a large distal tarsal c-d. a large proximal tarsal element comprising the fibulare, intermedium and tibiale + lateral centrale develops in the ankle. in p. vitticeps e-j, a very similar sequence of cartilage elements in the hand develop relative to a. uniparens with the added intermediate stage of seeing the medial carpale segment from the proximal end of distal carpal while the lateral centrale shifts distally adjacent to dc and dc the tarsus is more difficult to follow due to a small sample size and diffuse staining, but the proximal tarsal element is at least formed by the tibiale, intermedium and fibulare with centralia not distinguishable during early embryogenesis, carpal elements and tarsal elements form in chamaeleo calyptratus autopodia. the hands and feet at dpo are dorsoventrally flattened, the plesiomorphic condition for tetrapods. during later stages of embryogenesis, fusion of proximal elements and a reorganization of the distal limb skeleton gives way to a distal cleft and zygodactyly. in the adult manus of c. calyptratus, the proximal elements of the first row are retained, while only distal carpal elements are present, dc, dc in the pes, the fibulare, intermedium, and tibiale fuse, leading to a single ossification center comprising the proximal tarsal element, while distally, the distal tarsals, dt dt dt remain as independent bony elements. while observing the formation of the carpals, the ulnare was noted to undergo a segmentation event from its ventrolateral surface at dpo. subsequent analysis with later stage embryos provided evidence that this particular element later ossified as the pisiform bone. the pisiform is a large ossified element on the posterior margin of the manus in leaf mimic chameleons. the pisiform increases in size as one samples higher taxa across the chameleon phylogeny. thus, the origin of the pisiform as a sesamoid is not supported by our study. both the pisometacarpal ligament and the tendon of the flexor carpi ulnaris are visible with separate insertions, not enclosing the pisiform within a tendon as expected for a sesamoid in adult chameleons. this conclusion contrasts with a recent study which also recognized the pisiform’s origin from the ulnare, but declared it as a sesamoid. the bony element of the posterior tarsus is widely considered a sesamoid but was not apparent until after hatching in our study of c. calyptratus. this element clearly originates within a tendon and supports the effect of biomechanical stress in inducing the skeletogenic program of sesamoids. additional file 5: figure s5a-c illustrates the segmentation of the pisiform from the ventrolateral wall of the ulnare while transverse histological sections show the pisiform segmenting from the ulnare with a shared outer condensed layer of mesenchymal cells, the prospective perichondrium.fig. contentious skeletal elements of the chameleon mesopodium. while broadly considered to be a sesamoid, we see the pisiform segmenting from the ventrolateral margin of the ulnare, underscoring its origin from the primary skeleton. in ‘true chameleons’ we see what appear to be ectopic and possibly de novo elements on the posterior tarsus. b the two distal elements are actually the tibiale and the intermedium which become shifted distally toward distal tarsals by the enlarged distal epiphysis of the tibia. c despite the highly modified autopodium, chameleons retain distal tarsal 5 broad mesopodial skeletal diversity across chamaeleonidae to understand how skeletal fusions and heterochronic changes during skeletogenesis in the mesopodium of chameleons have been modified by natural selection, due to the need to adapt to an arboreal lifestyle, we sampled species representing of the currently recognized genera of extant chameleons. palleon was examined from an embedded 3d microct isosurface reconstruction in glaw. despite an overall conserved body plan, chameleon body size ranges from the large tree adapted furcifer oustaleti and calumma parsonii, at approximately 68 cm total length, to the recently described brookesia micra, whose maximum total length is less than 30 mm, a striking fold difference. within the family chamaeleonidae, genera diverged early in madagascar with all other genera forming the sister clade. the two malagasy genera and two basal genera found within the sister clade are known as the leaf chameleons as they are leaf mimics, whose niche encompasses ground cover and low shrubs in madagascar and africa. classically, morphological descriptions of chameleon hand and foot skeletons came from two derived and highly arboreal genera, trioceros and chamaeleo, thus biasing our knowledge of morphological diversity of autopodia across chameleons. this has inadvertently led to an underappreciation of skeletal element number with respect to carpals and tarsals in chameleons generally. our comparative analysis illustrates that the most basal lineages of chameleons maintain the least number of independent carpal and tarsal elements as adults. in the forelimb mesopodium, the lowest number of independent skeletal elements are found in the genera brookesia, rieppeleon and rhampholeon. members of these genera have a mesopodial formula of 3/ representing a /). all chameleons examined lack the middle row comprised of two centralia which are found in outgroup lepidosaurs. an isosurface microct 3d rendering of palleon lolontany grossly presents with the same autopodial morphology as brookesia. the foot mesopodial skeleton was greatly simplified in these three genera. the most reduced architecture is found in rieppeleon brevicaudatus, having a single distal tarsal and two proximal elements of which one was a sesamoid within the fibularis brevis tendon, whereas the other was a fused tibiale + intermedium + fibulare. surprisingly, in rieppeleon the proximal face of all metatarsals fuse with adjacent metatarsals to form a single crown surrounding the distal tarsal element. brookesia and palleon were found to be the only taxa sampled lacking the fibularis brevis sesamoid ossification. this may be representative of basal chameleons and/or locomotory mode.fig. histological examination of mesopodial mesenchymal and cartilaginous condensations. to complement the wholemount alcian cartilage stains, standard h&e paraffin sections were made to look for condensations which may not have shown due to a failure to reach the chondrification stage. in the carpals, we were able to find additional support for elements such as the intermedium, radiale and ulnare. in addition, we were able to identify a condensation which appears to be small and transient as distal carpal which does not chondrify and has a significant contribution of mesenchyme from the preaxial side of the limb. in the tarsals, we see a significant amount of shared mesenchymal cells between the distal tarsal fibulare, intermedium and the tibiale, with the tibiale and intermedium always in closer association and smaller from this study we can conclude that increased body length in chameleons >20 cm accompanies a corresponding increase in the number of carpal elements, from 3/ in leaf mimic chameleons, to 3/ 3/ 4/ and 4/ in the phylogenetically derived ‘true chameleons’. hindlimb mesopodial skeletal elements also increase from 1/ [proximal rows with sesamoid in parentheses) in leaf mimics, through 2/, 2/, 3/, with variation in chamaeleo spanning 2/, 1/, and 2/. all but brookesia and palleon in this study retain a tarsal sesamoid. as in other lepidosaurs, the foot skeleton is always composed of fewer elements than the manus. in the hindlimb mesopodium, we observed 1– elements which appeared to be ectopic or de novo structures in bradypodion, chamaeleo, furcifer and trioceros. knowing that in c. calyptratus the proximal tarsal is initially composed of condensations for the tibiale, intermedium and fibulare, we examined cartilage development in a dpo c. calyptratus embryo to follow the possible origin of these ectopic elements. at this stage, the distal epiphysis of the tibia has grown extensively in the distal direction and displaced the cartilages of the intermedium and the tibiale toward the metapodial skeleton. in adults, the tibiale, intermedium and fibulare fuse, forming a projection medially on the proximal tarsal element. the tibia is able to redirect the placement of the tibiale and intermedium from their proximal position. using these results as a model, the distal mesopodial element in bradypodion would be the tibiale and the proximal tarsal element would be comprised of a fused intermedium and fibulare. across the genus chamaeleo, 0– ectopic elements are present distally while are visible in both furcifer and trioceros. using the distal epiphyseal growth of the tibia as a model for proximal tarsal elements as seen in c. calyptratus, it is evident that the distal elements in the former two genera are due to a more distal tibiale and a more proximal intermedium, with the proximal tarsal element composed of only the fibulare. thus, distal displacement of the intermedium and tibiale is due to an accelerated growth heterochrony in the epiphysis as well as a loss of fusion for the three proximal tarsal elements which generally produce a single proximal tarsal. the proximal tarsal element in different chameleon taxa is also shown to have different compositions and are thus not homologous in our study. functional and phylogenetic constraints on mesopodial skeletal elements carpal and tarsal tunnels of the wrist and ankle, respectively, guide tendons, nerves and blood vessels from the zeugopod into the autopodial domain and are involved in flexion of the more distal skeletal elements. chameleons, unlike other tetrapods, have an autopodial skeleton with an approximately degree perpendicular fixed angle to the proximodistal axis of the appendicular skeleton. in association with a highly modified ‘ball-and-socket’ architecture of mesopodial elements, tissues extending through the tunnels must not only be enclosed through retinacula but must also be inhibited from sliding laterally and out of the tunnel by digits and during wrist extension, the position of which is fixed in chameleons. to compensate for the perpendicular orientation of the autopodium in chameleons, wrist elements of the hand have enlarged and grown ventromedially to laterally border the tunnel connective tissues. furthermore, the large distal carpal element develops a ventral concavity to funnel the connective tissues to the palmar side toward the ventral aspect of the metacarpals. in lepidosaurs, metatarsal and distal tarsal are fused. thus, a hooked ‘5th metatarsal’ serves as a synapomorphy for this clade. in the chameleon foot, ‘mt5’ and mt bilaterally border the connective tissues through the tarsal tunnel to the plantar surface of the foot. basal species exhibit a single proximal tarsal element and a large distal tarsal element. phylogenetically derived genera enclose the tarsal tunnel connective tissues by increasing the relative size of the fibularis brevis tendon while also making the tibiale and intermedium independent from the proximal tarsal element. this changes the composition of the proximal tarsal element across taxa in exchange for increased ventromedial support of the tarsal tunnel connective tissues. discussion chameleon autopodia have undergone extensive modification for optimized locomotion in a complex three-dimensional arboreal habitat. here we have taken advantage of the slowly developing veiled chameleon to characterize the morphogenesis of the squamate limb and study various aspects of morphogenesis of the hands and feet. this has aided our understanding of the diversification of autopodial morphologies across of the currently recognized chameleon genera and how the derived autopodial morphology of chameleons emerged within squamate reptiles. chameleons exhibit a typical pentadactyl morphology, but also present with a midline autopodial cleft and two opposable syndactylous bundles of highly mobile digits. in addition, the proximal skeleton of the autopodium in chameleons was modified via a reduction in the number of bone elements. due to the presence of digits in chameleons, the ancestral condition for tetrapods, we expected a conserved pattern of expression for both fgf and shh in the developing autopodia of chameleons. indeed, fgf was observed in the apical ectodermal ridge while shh was present in the zone of polarizing activity from the onset of limb bud outgrowth through termination of distal limb outgrowth and morphogenesis. all previous studies of autopodial ectrodactyly and clefting have concluded this phenotype as having occurred due to a loss of integrity of the aer. in contrast, we show that in chameleons, the aer remains a robust ectodermal thickening and expresses fgf during proximal mesenchymal clefting. interdigital cell death, also an ancestral condition for tetrapods, was mainly observed in chameleon autopodia between the digit pairs undergoing clefting while adjacent digits retained interdigital mesenchyme. interestingly, inhibition of bmp signaling prior to the initiation of interdigital cell death prevented cleft formation. however, upon the commencement of cleft formation, blocking bmp signaling has no effect on interdigital cell death. cell death between syndactylous digits was present in chameleons, but at significantly lower levels relative to the clefting mesenchyme. the syndactyly observed in chameleons may differ to that observed in duck feet and bat wings due to the absence of grem a bmp antagonist, in interdigitial tissue. while previous studies credited interdigital cell death along the autopodial midline for the modified chameleon limb morphology, we highlight the previously unstudied and underappreciated skeletal diversity present in the mesopodial architecture across chamaeleonidae. while retaining the ancestral pentadactyl morphology together with the typical phalangeal formulae, the chameleon wrists have been highly modified to have a reduced number of skeletal elements with subsequent enlargement of those elements that are maintained. proximal carpals and tarsals of digits and respectively, of the hands and feet also make contact with each other. we have shown that during embryogenesis in c. calyptratus, the carpals comprise skeletal elements and the tarsals each of which originate as mesenchymal condensations, chondrify, ossify or fuse with adjacent elements at various stages of skeletogenesis. we have used this model from c. calyptratus to homologize carpal and tarsal structures in of the currently recognized chameleon genera and found that chameleon research has significantly underestimated mesopodial skeletal diversity. the earliest diverged genera palleon + brookesia, and rieppeleon are not only the smallest members of the family, but also show the greatest degree of skeletal number reduction in the wrist and ankle. more phylogenetically derived chameleon genera from this study chamaeleo, trioceros, furcifer show the greatest number of independent skeletal elements present in the wrist and ankle within the family. this pattern also most closely resembles the skeletal complement found in outgroup lepidosaurs. a phylogenetic trend of gigantism appears to potentially correlate with the ‘re-evolution’ of mesopodial skeletal elements, as this is associated with an increased number of cells from which to build mesenchymal condensations in the mesopodia. an increase in mesopodial elements in ‘true chameleons’ also changed the resting posture of the hand and foot of chameleons, resembling a more typical lizard position where the two syndactylous bundles are oriented more anteriorly along the midline but are constrained due to the retainment of a cleft. chameleons therefore provide an important example of how heterochrony can modify developmental patterns while being constrained by functional requirements and possibly body size. thus, extrinsic and intrinsic factors are selected by nature to produce highly adaptive morphologies. loss of carpal and tarsal elements increases the angle across the cleft zone, with the greatest angle being present in the most basal and smaller genera. more derived chameleon genera have an independent intermedium in the ankle, a structure which has not been described in any extant terrestrial amniote since the origin of the astragalus and calcaneum from the amphibian skeletal complement. in the wrist, segmentation takes place along the ventrolateral wall of the ulnare and gives rise to the pisiform bone. in contrast to the sesamoid bone of the fibularis brevis tendon of the hindlimb, the pisiform appears to originate with the primary endochondral skeleton and not within a tendon. the family chamaeleonidae is thought to have originated through the heterochronic process of paedomorphosis, with adults resembling late gestation embryos of more distant taxa. malagasy species belonging to palleon and brookesia range in total length from 2–10 cm. here we have identified that early diverged and basal genera exhibit a dramatic reduction of independent carpal and tarsal elements in adults despite maintaining a normal pentadactyl number of digits and phalanges. this is in contrast to salamanders of equivalent size which exhibit digit loss and mesopodial element reduction. while we only sampled of the extant genera, an ascending phylogenetic trend of increased body size appears to correlate with an increase in mesopodial element number to a complement close to that of outgroup lepidosaurs. more components in the wrist and ankle may have facilitated increased wrist flexion and provided a biomechanical advantage which has allowed more derived chameleons to leave the ground cover and low bushes for the trees. clefting of the autopodium via interdigital cell death is the ancestral condition in the tetrapod limb. in chameleons, interdigital tissue is however retained between fingers neighboring the midline cleft. this represents the derived condition. in phylogenetically more derived chameleon lineages with ‘re-evolved’ mesopodial skeletal elements, the ancestral tetrapod condition where all digits are present along the distal midline of the limb is almost achieved. while clefting separates the syndactylous bundles of digits, it is clear from the transiently bent chondrified digital rays of the alcian prepared skeletons that there is either a musculoskeletal component driving the split, pulling the digit bundles toward the anteroposterior poles, or there may be a proliferative phase in the mesenchyme prior to cleft formation by apoptosis. mesopodial elements have not received as much attention in the development literature when compared to studies of digit loss and phalangeal reduction. despite the three major subdomains of the limb being highly conserved across tetrapods, species-specific differences in the morphometric dimensions of long bone elements of the zeugopodia may elicit a direct effect on carpal and tarsal elements. broader elements would provide a wider dimension of initial mesopodial mesenchymal pioneer cells for the more distal limb structures. we hypothesize that altered proximal limb width, especially in the zeugopodial domain, may affect chameleon evolution as a narrowing of this domain would ultimately provide a reduction in the distal mesenchymal cells available for building the mesopodial and acropodial skeleton. this model may also have been important for the evolution of digit and carpal reduced taxa such as hemiergis as they, like chameleons, have qualitatively thinner, elongated limbs with a smaller distal radius in zeugopodial skeletal elements. larger taxa such as mammals, which lack a bottleneck in cell number for the mesopodial element and acropodial modification, employ genetic mechanisms for modifying adult phenotypes. in contrast, small ectotherms appear to exploit the cell number constraint model and typically exhibit element size and element number variation. within this framework, a smaller mesopodial mesenchymal domain would be expected to give rise to a reduced number of carpal/tarsal elements as well as fewer cells able to contribute to digit notably, while two different studies have demonstrated a temporal reduction in shh expression in association with digit reduction, the tetradactyl pattern in these studies was different. in the mouse, the digits were identified as having the character of digits i, ii, iv, v, whereas the digit reduced lizard retained digits ii, iii, iv, v together with a reduction in mesopodial elements. digit i loss in hemiergis is consistent with our hypothesis that this phenotype is associated with a smaller initial pool of mesopodial mesenchymal cells. fusion and reduction are more prevalent in organisms of smaller size, as was seen in this study, and this may simply be due to fewer mesenchymal cells being available for individual condensations to form or reach a sufficient density required for establishment of the chondrogenic pathway. mesopodial elements do not ‘re-evolve’ de novo dollo’s law on the irreversibility of complex traits has recently been challenged in taxa where autopodial elements appear to have been re-evolved. phylogenetic trends in chameleon autopodial elements are not due to de novo ‘re-evolution’ of skeletal elements, but rather the elaboration of skeletal condensations already present. in members of the genera chamaeleo, furcifer and trioceros sampled, mesopodial elements appear to be present but are underdeveloped or fused to adjacent condensations or skeletal elements, character states which can be easily missed when a detailed developmental series is not available for study. mammalian species appear to reduce autopodial elements through genetically controlled cell death. while a reduced period of shh expression in the zpa may also decrease the population size of mesenchymal cells available for autopodial development, narrow limbs such as those in hemiergis lizards may have played an equally important role in reducing the cellular pool leading to element reduction. small mesenchymal population pools may also be a possible constraint on other members of the family scincidae such as brachymeles and larutia. more laboratory and field studies are needed to test this. differential fusion between mesopodial condensations in chameleons blurs their skeletal homologies with outgroup lepidosaurs, especially if the topographic position of an element of interest or the number of ossification centers present are the only criteria considered. herein we show that interspecific variability exists in not only the number of elements present but also that different combinations of mesopodial elements fuse to give rise to the adult carpal and tarsal complement along with distal migration of other elements. caution must be exercised when interpreting skeletal morphology using only dry skeletal specimens and microct imaging, as neither allows for observation of cartilaginous tissue. we have observed that fused condensations do not necessarily have multiple ossification centers, as individual condensations remain modular and arrest their skeletogenic program at the mesenchymal, chondrogenic or ossification stage.fig. chamaeleo calyptratus atuopodial skeletogenic condensations. despite a very specialized and element reduced adult chameleon autopodial skeleton, most ancestral skeletal elements remained present through ontogeny with the exception of the centralia. between the manus and pes, the preaxial elements were always significantly delayed in appearance. the intermedium appeared to have a ‘dual’ origin, with its condensation occurring closer to the radius in the forelimb and the fibulare in the hindlimb. we found that elements condensed and completed skeletogenesis. there was a high tendency for fusion amongst proximal elements of the foot in c. calyptratus despite delayed appearance of the tibiale. distal carpal and distal tarsal were seen to condense as mesenchyme but later fuse to adjacent elements while we were unable to find any condensation for distal tarsal distal carpal remains cartilaginous. the pisiform was found to be derived through segmentation of the ulnare and the fibularis brevis tendon to enclose a sesamoid fig. evolutionary trends in chameleon autopodial morphology. outgroup comparison was conducted with the extant rhynchocephalian sphenodon punctatus and outgroups within squamata ], agamidae ). fusion of proximal tarsal elements arose during amniote evolution by formation of the astragalus and the calcaneum. aspidoscelis and pogona both show an increase in fusion of mesopodial elements relative to sphenodon. the small and early diverged chameleons, palleon) show the distal clefting of the autopodia between the forelimb digits iii and iv and hindlimb ii and iii characteristic of chameleons. these two genera also present the greatest modified proximal autopodial skeleton within chameleons. the larger rieppeleon and rhampholeon ) presents a sesamoid in the fibularis brevis tendon of the heel which is maintained in all higher chameleons. bradypodion, chameleo, furcifer and trioceros sampled in this paper showed an ascending phylogenetic trend toward increased body size and an associated increase in skeletal elements present in the mesopodium, leading to a wrist with increased flexion capability as well as being the first amniote group to form an independently ossified tibiale and intermedium which is only known to be present in amphibians. thus, in early chameleon lineages there appears to be a skeletal bottleneck relative to outgroup squamates with a subsequent phylogenetic trend of ‘re-evolving’ mesopodial elements present in outgroup taxa. the red gradient shows a trend toward reduction of elements from outgroup taxa to chameleons while the red gradient shows a relative increase in mesopodial skeletal elements in advanced chameleon genera studied. genera in “<genus>” are currently recognized taxa but were not examined in this study. letters within figure correspond to the species with the same letters in additional file 1: figure s1 fig. phylogenetic reversal to a ‘typical’ lizard conformation of the autopodia. despite retaining the cleft, an increase in mesopodial elements leads to increase wrist flexion in higher chameleons. concomitant with the increase in the number of carpal and tarsal elements, we see a bifurcation where the syndactylous bundles are almost 180° separate in the early diverged lineages, rieppeleon, and rhampholeon ) with the angle decreasing in bradypodion. ‘true chameleon’ syndactylous bundles shift toward the distal midline as is present in the basal squamate plan. while chamaeleo show an anteromedial transition of digit bundles, furcifer and trioceros, the taxa with greater number of carpal and tarsal elements in this study, exemplify this the best through a greatly reduced angle between syndactylous bundles as digits move closer together along the distal midline. fl: refers to forelimb with hl: referring to hindlimb. the numbers associated with each lineage on the phylogeny correlate to the same order of illustrations from fig.  and represent formulae for number of elements present in the mesopodium in association with a representative image for autopodial morphology of a preserved specimen. fl: x/x refers to the number of elements in the two mesopodial rows while fl: x/x refers to the two mesopodial rows with the number in parentheses representing the presence or absence of the sesamoid of the fibularis brevis. no centralia are present in chameleons, which would have made the standard rows of mesopodial elements. genera in “<>” are currently recognized taxa, but were specimens were not obtained in this study homology and evolution of the astragalus and tarsal bones a proximal reduction in the number of carpal and tarsal elements in adult chameleons, relative to embryos and outgroup taxa, together with distal clefting, have allowed for evolution of the novel chameleon hand and foot which are used for complex grasping in an arboreal three dimensional environment. the progressive addition of skeletal elements in higher taxa may have facilitated increased flexion and a greater ability to grasp, which correlates with chameleons invading the arboreal niche where only mammals with autopodial modifications such as marsupials have evolved to inhabit. chameleons are thought to have invaded madagascar through two independent events. the initial invasion by the ancestor of brookesia and palleon evolved in situ, and was followed by a secondary invasion of the ancestor of calumma and furcifer by dispersal from africa. calumma + furcifer are the sister lineage to trioceros + kinyongia, which collectively comprise the four most phylogenetically derived extant chameleon genera. from our studies of the mesopodial architecture of these african and malagasy clades, we can infer that the dispersed ancestor of calumma + furcifer to madagascar was already mesopodially ‘pre-adapted’ to be an advanced climber due the larger complement of skeletal elements in the wrist and ankle. brookesia and palleon were already in madagascar as forest floor dwellers. given the advanced nature of the wrists in calumma, furcifer, trioceros and kinyongia, such increased mesopodial flexion has allowed them to comprise ~57 % of the extant chameleon species. innovations associated with the evolution of terrestrial cursorial ability in amniotes involved modifications to the ankle. evolution of the astragalus has elicited much discussion with recent models attempting to ascribe it an accurate homology with respect to tarsal skeletal elements found in anamniotes, generally involving 1– ossification centers. our work reveals a high degree of variability in the composition of the chameleon and squamate astragalus prior to the generally accepted co-ossification step following mesenchymal condensation. chameleons of the genera chamaeleo, furcifer and trioceros as shown here, are currently the only known amniotes to have successfully ‘re-evolved’ an independent intermedium and tibiale. a free intermedium and tibiale was lost in amniotes during the tetrapod progression onto land which selected for innovation of the astragalus and calcaneum. reptiles are considered to not have a tibiale. this ‘re-evolution’ of an ancestral series of tarsal elements may be due to a loss in selective pressure toward the typical terrestrial mode of tetrapod locomotion, in combination with a shift toward arboreality. within this framework, different amniote lineages should be expected to develop an astragalus through the fusion of tarsal elements at different points through ontogeny. although each element is modular, the distinct stage of skeletogenesis prior to and after fusion may vary across squamata. in summary, examining the development of innovations and novelties in natural groups such as chameleons can further our understanding of how body plans are molded by natural selection as optimal phenotypes through adaptive evolution.  studying organisms in the lab with unique morphologies complements biomedical studies of normal development and congenital malformations while also providing insight into the limits and novel deployment of molecular genetic pathways during embryonic development. veiled chameleons present a highly derived body plan within tetrapods due to their specialization for living in an arboreal environment while also presenting a very slow rate of embryonic growth. developmental analysis of the specialized and bifurcate autopodia of c. calyptratus using histological, in situ hybridization and in vivo techniques show that the specialized mesopodial skeleton forms despite the expression of highly conserved patterning genes during limb outgrowth. our detailed skeletal developmental sequence present an unexpected doubling in the number of recognized wrist and ankle elements which has allowed us to homologize mesopodial elements across of the recognized chameleon genera. our data also present a correlation between body size and skeletal diversity of the wrist and ankle during chameleon evolution with phylogenetically higher taxa exhibiting an increase in mesopodial skeletal element number, larger body size, and increased climbing ability. exploring development in a slow developing taxon has also allowed us to observe the fine progressive development and segmentation of the pisiform from the ulnare in the hand and differential fusion of proximal tarsal elements in the ankle. additional files additional file 1: figure s comparative morphology of autopodia in this study. dorsoventrally flattened autopodia in outgroup taxa were drawn in dorsal view only for both forelimb and hindlimb autopodia. d-n are chameleon illustrations from taxa cleared and stained in this study to look for ossified and chondrified elements comprising the autopodial architecture and are illustrated in both anterior/dorsal and posterior/ventral views. a) sphenodon punctatus, b) aspidoscelis uniparens, c) pogona vitticeps, d) brookesia stumpffi, e) rieppeleon brevicaudatus, f) rhampholeon boulengeri, g) bradypodion pumilum, h) chamaeleo dilepis, i) chamaeleo calyptratus, j) chamaeleo chameleon, k) furcifer lateralis, l) furcifer pardalis, m) trioceros ellioti, n) trioceros hohnelii. additional file 2: figure s forelimb condensations of a. uniparens in h&e stained sections. a large contribution of preaxial mesenchymal cells contributing to digit 1/distal carpal are visible in these sections as well as the development of an independent lateral centrale and intermedium at an older stage. the medial centrale will later appear as a medial segmentation from the distal carpal additional file 3: figure s pisiform versus sesamoid identity in the chameleon autopodium. the top row shows the forelimb with the pisiform bone which we have identified in this study as segmenting from the ulnare. this bony element is not present within a tendon and instead has a proximal insertion of the tendon of the flexor carpi ulnaris with a more distal pisometacarpal ligament. the lower row shows the hindlimb with a large sesamoid present within the tendon of the fibularis brevis which inserts onto the proximal margin of the ‘5th metatarsal.’ the species used for this image was trioceros ellioti. additional file 4: figure s whole embryo culture of c. calyptratus. details for this protocol are provided in the article text in the methods section and correlate with labeled figures. additional file 5: figure s pisiform formation in the manus of c. calyptratus. pisiform development from the ulnare is shown in posteroventral, lateral, and dorsal view in whole alizarin/alcian skeletal preparations. transverse h&e stained paraffin sections show, from distal to proximal mesopodium, the shared mesenchymal cells between the ulnare and the pisiform, supporting the ventrolateral segmentation of the pisiform from the ulnare. additional file 6: figure s autopodia are dorsoventrally flattened divergent tetrapods such as the common lab mouse and the desert grassland whiptail lizard have divergent hands and feet. despite architectural differences, they remain dorsoventrally flattened. additional file 7: chameleon specimen measurements. abbreviations gremgremlin bmpbone morphogenetic protein shhsonic hedgehog fgffibroblast growth factor icdinterdigital cell death aerapical ectodermal ridge zpazone of polarizing activity dpodays post oviposition pteproximal tarsal element svlsnout-vent length tltail length flforelimb hlhindlimb roman numerals i-vrepresent digit 1- digit 5 raradiale rradius uulna ululnare ppisiform iintermedium ttibia titibiale ffibula fifibulare dcdistal carpal dtdistal tarsal mcmetacarpal mtmetatarsal lacelateral centrale mecemedial centrale competing interests the authors declare that they have no competing interests. authors’ contributions red and pat conceived and designed the study. red carried out the laboratory work and data collection. red and pat drafted, read and approved the final manuscript. 
to better understand the homology of chameleon carpal and tarsal elements, the process of syndactyly, cleft formation, and how modification of the mesopodial skeleton has played a role in the evolution and diversification of chameleons, we have studied the veiled chameleon. phylogenetically derived chameleons exhibit a ‘typical’ outgroup complement of mesopodial elements, with twice the number of currently recognized carpal and tarsal elements considered for this clade. our study uncovered a previously underappreciated degree of mesopodial skeletal diversity in chameleons. we analysed limb patterning and morphogenesis through in situ hybridization, in vitro whole embryo culture and pharmacological perturbation, scoring for apoptosis, clefting, and skeletogenesis. body size may have played a role in the characteristic mesopodial skeletal architecture of chameleons by constraining deployment of the skeletogenic program in the smaller and earliest diverged and basal taxa. our study challenges the ‘re-evolution’ of osteological features by showing that ‘re-evolving’ a ‘lost’ feature de novo may instead be due to so called ‘missing structures’ being present but underdeveloped and/or fused to other adjacent elements whose independence may be re-established under changes in adaptive selective pressure. interdigital cell death therefore may be an ancestral characteristic of the autopodium, however syndactyly is an evolutionary novelty. one of the most distinctive traits found within chamaeleonidae is their split/cleft autopodia and the simplified and divergent morphology of the mesopodial skeleton.
 one of the most distinctive traits found within chamaeleonidae is their split/cleft autopodia and the simplified and divergent morphology of the mesopodial skeleton. these anatomical characteristics have facilitated the adaptive radiation of chameleons to arboreal niches. to better understand the homology of chameleon carpal and tarsal elements, the process of syndactyly, cleft formation, and how modification of the mesopodial skeleton has played a role in the evolution and diversification of chameleons, we have studied the veiled chameleon. we analysed limb patterning and morphogenesis through in situ hybridization, in vitro whole embryo culture and pharmacological perturbation, scoring for apoptosis, clefting, and skeletogenesis. furthermore, we framed our data within a phylogenetic context by performing comparative skeletal analyses in of the currently recognized genera of extant chameleons.  our study uncovered a previously underappreciated degree of mesopodial skeletal diversity in chameleons. phylogenetically derived chameleons exhibit a ‘typical’ outgroup complement of mesopodial elements, with twice the number of currently recognized carpal and tarsal elements considered for this clade. in contrast to avians and rodents, mesenchymal clefting in chameleons commences in spite of the maintenance of a robust apical ectodermal ridge. furthermore, bmp signaling appears to be important for cleft initiation but not for maintenance of apoptosis. interdigital cell death therefore may be an ancestral characteristic of the autopodium, however syndactyly is an evolutionary novelty. in addition, we find that the pisiform segments from the ulnare and that chameleons lack an astragalus-calcaneum complex typical of amniotes and have evolved an ankle architecture convergent with amphibians in phylogenetically higher chameleons.  our data underscores the importance of comparative and phylogenetic approaches when studying development. body size may have played a role in the characteristic mesopodial skeletal architecture of chameleons by constraining deployment of the skeletogenic program in the smaller and earliest diverged and basal taxa. our study challenges the ‘re-evolution’ of osteological features by showing that ‘re-evolving’ a ‘lost’ feature de novo may instead be due to so called ‘missing structures’ being present but underdeveloped and/or fused to other adjacent elements whose independence may be re-established under changes in adaptive selective pressure. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. issue-copyright-statement© the author 2015 