Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

A correction for this article has been published in BMC Evolutionary Biology 2010, 10:209

Open Access Research article

The mitochondrial phylogeny of an ancient lineage of ray-finned fishes (Polypteridae) with implications for the evolution of body elongation, pelvic fin loss, and craniofacial morphology in Osteichthyes

Dai Suzuki1, Matthew C Brandley2 and Masayoshi Tokita13*

Author Affiliations

1 Department of Zoology, Graduate School of Science, Kyoto University, Sakyo, Kyoto, 606-8502 Japan

2 Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 06520-8105 USA

3 Graduate School of Life and Environmental Sciences, University of Tsukuba, Tsukuba, Ibaraki, 305-8572 Japan

For all author emails, please log on.

BMC Evolutionary Biology 2010, 10:21  doi:10.1186/1471-2148-10-21

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/10/21


Received:22 April 2009
Accepted:25 January 2010
Published:25 January 2010

© 2010 Suzuki et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The family Polypteridae, commonly known as "bichirs", is a lineage that diverged early in the evolutionary history of Actinopterygii (ray-finned fish), but has been the subject of far less evolutionary study than other members of that clade. Uncovering patterns of morphological change within Polypteridae provides an important opportunity to evaluate if the mechanisms underlying morphological evolution are shared among actinoptyerygians, and in fact, perhaps the entire osteichthyan (bony fish and tetrapods) tree of life. However, the greatest impediment to elucidating these patterns is the lack of a well-resolved, highly-supported phylogenetic tree of Polypteridae. In fact, the interrelationships of polypterid species have never been subject to molecular phylogenetic analysis. Here, we infer the first molecular phylogeny of bichirs, including all 12 recognized species and multiple subspecies using Bayesian analyses of 16S and cyt-b mtDNA. We use this mitochondrial phylogeny, ancestral state reconstruction, and geometric morphometrics to test whether patterns of morphological evolution, including the evolution of body elongation, pelvic fin reduction, and craniofacial morphology, are shared throughout the osteichthyan tree of life.

Results

Our molecular phylogeny reveals 1) a basal divergence between Erpetoichthys and Polypterus, 2) polyphyly of P. endlicheri and P. palmas, and thus 3) the current taxonomy of Polypteridae masks its underlying genetic diversity. Ancestral state reconstructions suggest that pelvic fins were lost independently in Erpetoichthys, and unambiguously estimate multiple independent derivations of body elongation and shortening. Our mitochondrial phylogeny suggested species that have lower jaw protrusion and up-righted orbit are closely related to each other, indicating a single transformation of craniofacial morphology.

Conclusion

The mitochondrial phylogeny of polypterid fish provides a strongly-supported phylogenetic framework for future comparative evolutionary, physiological, ecological, and genetic analyses. Indeed, ancestral reconstruction and geometric morphometric analyses revealed that the patterns of morphological evolution in Polypteridae are similar to those seen in other osteichthyans, thus implying the underlying genetic and developmental mechanisms responsible for those patterns were established early in the evolutionary history of Osteichthyes. We propose developmental and genetic mechanisms to be tested under the light of this new phylogenetic framework.

Background

Osteichthyans (bony fish and tetrapods; Fig. 1) have evolved remarkably diverse body plans since their initial radiation in the Late Silurian ~420 Mya [1,2]. It is therefore not surprising that most major extant lineages have been the subject of extensive evolutionary biology research. As a result, we know much about the evolutionary history and patterns of morphological evolution in osteichthyans, most notably teleost fish and tetrapods (amphibians, reptiles, and mammals). In the age of genomics and advanced molecular techniques, knowledge of these relationships and patterns has proven useful in uncovering the developmental and genetic mechanisms responsible for morphological diversity [e.g., [3-9]].

thumbnailFigure 1. Consensus view of the phylogeny of extant Osteichthyes, including the position of bichirs (Polypteridae), inferred from phylogenetic analyses of molecular and morphological data.

However, the same cannot be said for one lineage that diverged early in the evolutionary history of Actinopterygii (ray-finned fishes) - the family Polypteridae. Commonly known as "bichirs", this clade includes 12 extant described species (as well as multiple subspecies) that inhabit freshwater rivers and lakes of tropical Africa [10]. Little is known about polypterid evolution primarily because of their long history of phylogenetic and taxonomic confusion, likely owing to their unique collection of "primitive" (e.g., ganoid scales, cartilaginous skeleton, the intestine with a spiral valve) and derived (e.g., highly modified dorsal fins, pectoral fins with lobed base covered with scales, possession of only four gill arches) anatomical features [11-13]. However, recent morphological and molecular analyses have determined that bichirs are a basal lineage of Actinopterygii (ray-finned fishes; Fig. 1) [14-19].

Bichirs provide an attractive system to test whether patterns of morphological evolution are shared throughout the osteichthyan tree of life. If common patterns exist, it would suggest that the underlying genetic and developmental mechanisms responsible for those patterns were established early in the evolutionary history of Osteichthyes. For example, there exist multiple independent derivations of body elongation in modern Actinopterygian fish [e.g., [20-23]], amphibians [e.g., [24]], squamate reptiles [e.g., [25-27]], and numerous extinct tetrapod lineages [e.g., [28,29]]. The loss of pelvic fins and limbs has occurred multiple times over the course of osteichthyan evolution [22,27,30-32]. Furthermore, diverse changes in craniofacial morphology have been the subject of particularly extensive research [5,33-41].

However, the greatest impediment to elucidating the patterns of morphological change in Polypteridae is the lack of a well-resolved, highly-supported phylogenetic tree. In fact, the interrelationships of polypterid species have never been subject to molecular phylogenetic analysis. Those molecular studies that have included polypterid representatives did so in the context of determining "deep" osteichthyan and actinopterygian relationships (Fig. 1) or surveying Hox gene clusters [42].

Extant polypterids comprise two extant genera, Erpetoichthys and Polypterus. Boulenger [43] distinguished these two genera based on the extremely elongate body and absent pelvic fins in Erpetoichthys, and split Polypterus into two taxonomic groups based on the position of the mandible relative to the snout ("lower jaw protrusion" hereafter). Poll [44-46] defined five clusters of species and an ancestral species within extant polypterids based on anatomical features such as relative jaw length, the location and size of the eyes, width of the suboperculum, and proportion of the gular plate. An analysis of 15 anatomical characters suggested, among other things, a sister relationship between Erpetoichthys and P. weeksi, thus rendering Polypterus polyphyletic [47]. However, this study did not incorporate objective rooting or optimality criteria and subsequent parsimony reanalysis revealed that this phylogeny is largely unresolved [12]. These conflicting and inconclusive results warrant additional phylogenetic scrutiny using molecular data.

Here, we infer the first molecular phylogeny of Polypteridae using all 12 recognized species including all but three recognized subspecies to achieve three major objectives. First, we evaluate whether the current species taxonomy underestimates the underlying phylogenetic diversity. Secondly, we examine patterns of evolution in body elongation pelvic fin reduction, and craniofacial morphological evolution in this new phylogenetic context. Finally, we recommend developmental and genetic mechanisms to be tested in the future in conjuction with this new phylogenetic framework.

Methods

Phylogenetic analyses

We obtained fresh tissues from commercial dealers of 15 recognized species and subspecies of polypterids including Erpetoichthys calabaricus, Polypterus ansorgii, P. bichir bichir, P. b. lapradei, P. delhezi, P. endlicheri congicus, P. e. endlicheri, P. mokelembembe, P. ornatipinnis. P. palmas buettikoferi, P. p. polli, P. retropinnis, P. senegalus senegalus, P. teugelsi and P. weeksii (Table 1). One to four specimens of each taxon were included in the phylogenetic analyses. Although we included all recognized species of polypterids, we were unable to obtain tissue samples for three subspecies, P. bichir katangae, P. palmas palmas, and P. senegalus meridionalis (Table 1). We obtained DNA sequences of outgroup taxa (Sturgeon [Acipenser transmontanus], Bowfin [Amia calva], Coelacanth [Latimeria chalumnae], Spotted gar [Lepisosteus oculatus], Paddlefish [Polyodon spatula], and Lungfish [Protopterus dolloi]) from GenBank (Additional file 1).

Additional file 1. Information of DNA sequences analyzed in this study.

Format: PDF Size: 586KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Table 1. Number of vertebrae and condition of the lower jaw and pelvic fins in each species and subspecies of Polypteridae.

DNA samples were extracted following the same procedure used in Tokita et al. [48]. A part of the mitochondrial 16S rRNA and cytochrome b (cyt b) genes were amplified by PCR System GeneAmp 2700 (Applied Biosystems, Lincoln, USA) using an Ex Taq polymerase kit (Takara Shuzo Co., Ltd., Otsu, Japan) and the primers provided in Additional file 2. After initial denaturing for 5 min at 95°C, 25-30 cycles were performed with denaturing 1 min at 94°C, annealing 3 min at 50-60°C, and primer extension for 3 min at 72°C, followed by a final elongation of 7 min at 72°C. PCR products were purified by PEG/NaCl precipitation and then sequenced using a Big Dye Terminator Cycle Sequencing Ready Reactions Kit v1.0 and an ABI PRISM 377 and 3130 DNA Sequencer (Applied Biosystems, Lincoln, USA), using the primers provided in Additional file 2.

Additional file 2. List of primers used in this study.

Format: PDF Size: 298KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Because of conserved codon reading frames, the cyt b protein coding sequences could be unambiguously aligned by eye. The 16S rRNA was also aligned by eye, but nucleotide positions that could not be unambiguously aligned were excluded from phylogenetic analysis (data set available from the authors). In some regions of the 16S data, the ability to align the polypterid taxa was improved if we removed sequences for the outgroup taxa (and their corresponding data replaced with "?" in these regions). Because our primary goal is to test the interrelationships of polypterids, and not necessarily their relationship to other osteichthyans, we feel the exclusion of these data is justified. All DNA sequences were deposited in GenBank (Additional file 1). The final data set for subsequent phylogenetic analysis included 925 base pairs (bp) of 16S (256 parsimony informative characters) and 932 bp of cyt b (103, 29, and 292 parsimony informative characters for the first, second, and third codon positions), for a total of 1857 total characters.

We employed partitioned Bayesian analyses to infer the phylogenetic relationships among polypterids. Previous studies have demonstrated that applying separate models of nucleotide evolution to specific subsets of nucleotide data (i.e., "partitioned" or "mixed-model" analyses) improves phylogenetic reconstruction [49-51]. We therefore partitioned the data a priori into four total partitions, one for the 16S rRNA and three for each codon position of cyt b. We determined the appropriate models of sequence evolution from a pool of 16 models (K80, HKY, SYM, GTR with and without I and Γ) for each partition using a JC corrected neighbor-joining tree in PAUP* [52] and the Bayesian Information Criterion [BIC: [53]]. We subsequently conducted Bayesian phylogenetic analyses using the most appropriate models (GTR+I+Γ, GTR+I+Γ, HKY+I+Γ, and GTR+Γ for the 16S and three cyt b codon positions, respectively) in the parallel version of MrBayes v3.1.1 [54,55]. Each Bayesian analysis consisted of 107 generations, using four chains sampled every 1000 generations and default priors (substitution rates, Dirichlet [1, 1, 1, 1, 1, 1]; base frequencies, Dirichlet [1, 1, 1, 1]; gamma shape parameter, uniform [0, 200]; topologies, uniform; branch lengths, unconstrained exponential [λ = 10]).

To determine convergence of the Bayesian analyses, we constructed cumulative posterior probability plots for each clade using the cumulative function in AWTY [56]. Stationarity was assumed when the cumulative posterior probabilities of all clades stabilized. These plots indicated that excluding the first two million generations as burn-in was sufficient. To decrease the chance of reaching apparent stationarity on local optima, four separate analyses were performed. Posterior probability estimates for each clade were then compared between the four analyses using a scatter-plot created by the compare command in AWTY. Posterior probability estimates for clades were similar in all four analyses, and the results of the analyses were combined and the frequency of inferred relationships in these 32,000 trees represented estimated posterior probabilities of clades. Posterior probabilities ≥ 0.95 are considered statistically significant (i.e., "strong") clade support [57].

Morphometric analysis

To investigate shape differences in the head region among Polypterus species, we employed landmark-based geometric morphometrics [58,59]. This approach allows us to describe morphological transformation of certain organismic structures visually and compare morphological differences among taxa quantitatively [59] and has been widely used in contemporary evolutionary biology research [e.g., [34,40,41]]. A total of 13 taxa and 26 specimens caught in wild were included in this analysis. Photographs of both dorsal and lateral views were taken with a digital camera. The photo images of the specimens were imported into the tpsDig2 software [60]. Eight and ten landmarks were placed for dorsal and lateral views, respectively (Additional file 3). The generalized least squares Procrustes superimposition, which removes information unrelated to shape such as position, scale, and orientation, was carried out using Coordgen 6n [61]. Based on the coordinates after the Procrustes superimposition, the partial warp scores were computed and principal component analysis (PCA) was conducted by calculating PCs based on the covariance matrix derived from the partial warp scores using PCAgen6 h [62].

Additional file 3. Landmark definition for morphometric analysis. (A) Polypterus palmas buettikoferi in dorsal view, (B) P. endlicheri congicus in dorsal view, (C) P. p. buettikoferi in lateral view, (D) P. e. congicus in lateral view.

Format: PDF Size: 638KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Ancestral state reconstruction

To evaluate the evolution of body elongation and pelvic fin loss within polypterids, we reconstructed ancestral states using data from extant sampled taxa, the rooted mtDNA phylogram (using mean branch lengths from the posterior distribution) estimated from our Bayesian phylogenetic analyses (above), and Mesquite v2.6 [63]. We removed the non-polypterid outgroups and pruned the data set and trees to include a single exemplar of each species and subspecies. We collected vertebral number for each species from FishBase [64] and preserved specimens. Polypterid species typically exhibit variation in the number of vertebrae (Table 1). Because there is currently no statistically elegant way to incorporate this information into ancestral state reconstruction analyses (e.g., modeling vertebral number in each species as statistical distributions) we used the mean number of vertebrae and coded it as a continuous character [e.g., [65]]. Ancestral vertebrae number for each node was subsequently estimated using weighted squared-change parsimony [66]. Presence of pelvic fin was coded as discrete binary characters and the marginal probabilities of each state at each node were estimated with maximum likelihood [67,68] using the 1-parameter symmetrical Markov model [69]. Ancestral states with a marginal probability ≥ 0.95 were considered resolved with statistical significance.

Results and Discussion

Polypterid phylogeny and hidden diversity

Posterior probabilities (PP) for a vast majority of nodes are statistically significant (≥ 0.95) and the clades with probabilities less than 0.90 are relationships between populations of subspecies (Fig. 2). The analyses support a basal split between Erpetoichthys and Polypterus. In Polypterus, the recently described P. mokelembembe [69] is the sister taxon to all other members of the genus; although support for this relationship is very high (PP = 0.94), it is not significant. The remaining Polypterus are divided into two major clades, one of which (P. bichir + P. endlicheri + P ansorgii) is characterized by lower jaw protrusion (Table 1; see below).

thumbnailFigure 2. Molecular phylogeny of the extant polypterid species inferred from partitioned Bayesian analyses 16SrRNA and cyt b mitochondrial genes. Branch lengths are means of the posterior distribution. Numbers above or below the node indicate the Bayesian posterior probability that clade is correctly estimated given the model. Posterior probabilities less than 0.50 are not shown. Colors indicate groups defined in Fig. 4.

These results reinforce the findings of Boulenger [43] who recognized Erpetoichthys and Polypterus as distinct genera, in addition to two groups within Polypterus based on the condition of the lower jaw (P. mokelembembe was not yet described). The phylogenetic position of Erpetoichthys is not particularly surprising considering numerous biological aspects of this species. The skull of E. calabaricus possesses a variety of unique characteristics that are not found in Polypterus [12]. Furthermore, the ecology of E. calabaricus quite different than Polypterus as it inhabits coastal and estuarine areas unlike Polypterus species that primarily inhabit freshwater [10,70-72].

Our data support Boulenger's conclusion of a close relationship between P. weeksi and P. ornatipinnis [43]. However, the present analyses reveal phylogenetic relationships that are substantially different from other previous studies focusing on the detailed interrelationships of the species. Analyses of 15 morphological and anatomical traits [47] identified close relationships among P. ornatipinnis, P. weeksii, and Erpetoichthys calabaricus, and among P. senagalus, P. retropinnis, P. palmas (although subsequent parsimony reanalysis of these data inferred an almost completely unresolved tree [12]). However, neither of these relationships was inferred by our analyses.

Poll [44-46] recognized close relationships among P. ornatipinnis, P. weeksii, and P. delhezi and that between P. senagalus and P. retropinnis based on morphological characteristics (e.g., relative jaw length, the location and size of the eyes). Our results statistically reject these phylogenetic hypotheses (i.e., alternative hypotheses have a posterior probability > 0.95). Our phylogenetic analyses also statistically reject Poll's hypothesis [45] that P. ansorgii is the sister lineage to all other Polypterus instead, this position is represented by P. mokelembembe.

Our phylogenetic analyses also demonstrate that at least two species, as currently described, are not natural groups - P. endlicheri and P. palmas. Polypterus endlicheri is paraphyletic with respect to both P. bichir and P. ansorgii. Taxonomic changes are clearly needed and we recognize P. e. congicus and P. e. endlicheri as distinct species, P. congicus and P. endlicheri. Hanssens et al. [73] recognized three subspecies in P. palmas: P. p. palmas, P. p. polli, P. p. buettikoferi based on morphometrical inferences. In our phylogeny, the sister relationship between P. p. polli and P. delhezi is strongly supported (PP = 1.0) as well as the relationship between P. p. buettikoferi and P. teugelsi (PP = 1.0). Thus, the morphometric resemblance between P. p. buettikoferi and P. p. polli is clearly the outcome of evolutionary convergence. To understand the evolutionary history of P. palmas complex and make subsequent taxonomic changes, it is necessary to include P. p. palmas in future phylogenetic analyses.

These results imply that morphological characteristics such as body proportion and body coloration used in previous studies are not suitable for inferring phylogeny, possibly because states of these characters tend to be flexibly altered depending on the condition of the species' environment. Furthermore, these phylogenetic results not only have implications for morphological evolution within Polypteridae (below), but also reveal the hidden genetic diversity of this ancient clade of enigmatic fish. That at least two species are polyphyletic (P. endlicheri and P. palmas) is evidence that the current taxonomy of Polypterus underestimates its genetic diversity.

Perhaps more importantly, we find that P. mokelembembe, a species only recently described [74], represents the sister lineage to all other members of the genus. This result suggests that much of the evolutionary diversity of extant Polypteridae has yet to be fully discovered.

As both the 16S and cyt b mitochondrial genes are maternally inherited as a single unit, we acknowledge that our phylogeny specifically tracks the evolutionary history of the mitochondrial genome, and not necessarily the species history. However, we note that there are no obvious instances of recent introgression (e.g., [75]) as all of the named taxa are monophyletic. However, we cannot exclude the possibility that incomplete lineage sorting obscures deeper relationships (see [76] for a review), but we also note that this is unlikely given the rapid coalescence of mitochondrial DNA. Furthermore, we note that, although we have not sampled the entire mitochondrial genome for this analysis, the fact that almost every node in our phylogeny is significantly supported (posterior probability > 0.95) is evidence that we have indeed captured almost as much "phylogenetic signal" as is possible from the mitochondrial genome. Nonetheless, for clarity, we refer to the phylogeny as the "mitochondrial phylogeny" to distinguish it from the true (albeit unknowable) species tree, and acknowledge that a more confident species tree of polypterids awaits future analyses of many more independently evolving loci.

Evolution of body elongation

Our ancestral state reconstructions (Fig. 3) estimate numerous increases and decreases of vertebral number, and thus, changes in body elongation. However, the most striking transformation is seen in Erpetoichthys. Erpetoichthys possesses between 110 and 113 vertebrae (Table 1), but the vertebral number reconstruction for the most recent common ancestor (MRCA) of this genus and Polypterus is ~71. It is not known whether the already relatively high number of vertebrae reconstructed for the MRCA is an artifact of the extremely high number of vertebrae in Erpetoichthys biasing the estimation of the root, or evidence that extinct stem polypterids were already relatively elongate. Regardless, Erpetoichthys is much more elongate than the estimated common ancestor of crown polypterids.

thumbnailFigure 3. Extant states and result of ancestral state reconstructions of vertebral number using squared-change parsimony.

Anguilliformity (eel-like elongated body form), such as that seen in Erpetoichthys, is seen across diverse lineages of vertebrates including ray-finned fishes, caecilian amphibians, and squamate reptiles [21,23-27,77], and an increase in vertebral number is considered as the common primary factor underlying such body form change [22]. Developmentally, vertebrae emerge from 'somites' that form one after another, with antero-posterior direction in the embryo. They bud off from the anterior end of unsegmented 'presomitic mesoderm' and this budding is regulated by a segmentation clock represented by oscillation of 'clock' genes such as Lunatic fringe [reviewed in [78,79]]. Recently, it was found that the tempo of the segmentation clock is accelerated in snake embryos, "ticking" around four times faster relative to growth rate than in shorter-bodied animals like chick and mouse and consequently generating a large number of somites [80].

Erpetoichthys has around 100 preanal vertebrae, which is almost double the number of other polypterid species [[22], pers. obs.]. Although few embryological studies of Erpetoichthys have been conducted [81,82], the yolk-sac larva and apterolarva of this species appear to be much more elongated and have more somites than those of Polypterus. It is plausible that an increase of vertebrate number in Erpetoichthys was brought about by a similar developmental mechanism to that exploited in snake embryogenesis [80,83], although the precise molecular mechanisms that regulate tempo of segmentation clock in vertebrate embryos have not been identified yet. The present study implies that body elongation observed in Erpetoichthys seems to be brought about in an early phase of polypterid evolution, possibly by accelerating the tempo of the segmentation clock twice as fast as that in other polypterids.

Although the extreme degree of body elongation observed in Erpetoichthys is not found in extant Polypterus, two species: P. bichir and P. teugelsi show a relative increase of vertebrate number (Table 1; Fig. 3). Although the total number of vertebrae is over 60 in both species, the mechanisms by which they achieve this body elongation are quite different. In P. bichir, the number of preanal vertebrae is less than 50 like that in majority of Polypterus species but the number of postanal vertebrae increases into over 14 [22]. On the other hand, P. teugelsi shows relative elongation of the trunk region, possessing at least 54 preanal vertebrae (pers. obs.). Instead, this species has only 10 postanal vertebrae keeping the condition in majority of Polypterus species. Thus, body elongation seems to have been achieved independently in two Polypterus species by employing different developmental mechanism that functions in separate developmental modules [22]. We note that this pattern is surprisingly similar to that seen in squamate reptiles (lizards and snakes) in which body elongation is also achieved by two distinct mechanisms elongating either the trunk or the tail [24-26].

Pelvic fin loss

In vertebrates, the body elongation tends to co-occur with a loss or reduction of paired appendages [21,24-26,29,84,85]. Following this evolutionary trend, Erpetoichthys, whose body is extremely elongated, also lacks pelvic fins (Table 1). Maximum likelihood ancestral state reconstructions strongly support the conclusion that pelvic fins were present in the most recent common ancestor of Erpetoichthys and Polypterus (Marginal probability = 0.97; data not shown), and that they were secondarily lost in the Erpetoichthys lineage. Given that the loss of pelvic appendages is frequent in multiple osteichthyan lineages, it is likely that the genetic mechanism responsible for pelvic fin loss in Erpetoichthys is similar to other members of this group. Translocation of Hox gene expression in trunk paraxial and/or lateral plate mesoderm is known to be the potential mechanism underlying limb elimination in vertebrates [[86,87]; but see recent paper by Wolstering et al. [88]]. However, this is a less likely explanation for pelvic fin loss in Erpetoichthys because such alteration of developmental program might also obliterate any vestiges of the pelvic girdle yet Erpetoichthys retains a pair of pelvic girdle elements (pers. obs.). Degeneration of apical ectodermal ridge (AER) and lack of expression of AER-associated genes might account for hindlimb truncation in lizards [89,90] and snakes [86], but this has never been recorded in non-tetrapods. Similarly, failure in establishing the Zone of Polarizing Activity (ZPA) has only been implicated in hindlimb loss in cetaceans [91]. The most promising potential mechanism responsible for pelvic fin loss in Erpetoichthys is the elimination of expression of transcription factor-encoding genes such as Pitx1 and Hoxd9 in the pelvic region - the mechanism that causes pelvic fin loss in some teleost fishes [31,32,92]. Unfortunately, the pattern of pelvic fin reduction in embryogenesis of Erpetoichthys has not been described precisely [82] and thus, discrimination of these various hypotheses is not yet possible. However, future studies of these genetic and developmental mechanisms can do so within the phylogenetic framework constructed here.

Transformation of craniofacial morphology

The first three PCs of the geometric morphometrical analysis of both dorsal and lateral views of the head accounted for over 70% of the total shape variation. In the comparison of the dorsal views, PC1 explains 36% of the total variance and primarily describes positional change of the orbit. PC2 explains 21% of the total variance (Additional file 4). In the comparison of the lateral views, PC1 explains 51% of the total variance and describes positional change of the tip of the mandible, the orbit, and labial end. While, PC2 explaining 14% of the total variance describes positional change of the opercle.

Additional file 4. Thin-plate splines (TPS) of principal component 1 (PC1) and 2 (PC2) for both dorsal and ventral views of the head of Polypterus. Arrows indicate PC value plus positive 0.1 score. (A) PC1 in dorsal view, (B) PC2 in dorsal view, (C) PC1 in lateral view, (D) PC2 in lateral view. Numbers along with each plot indicate landmarks defined in Additional file 3.

Format: PDF Size: 129KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Our morphometric analysis revealed the trend of morphological change in Polypterus species where PC1 positively correlates with the size of the head (Fig. 4a and 4c). Phylogenetically closely-related P. ansorgii, P. bichir and P. endlicheri occupied the most extreme positions on the positive sides of both PC1 and PC2 axes. This suggests that these three species not only tend to have a rostrally-projected lower jaw (a condition already noted by previous authors [43,47]), but also a medially located and up-righted orbit, enlarged mouth, and rostroventral distortion of the opercle. Judging from our mitochondrial phylogeny, the character complex seen in these three species is a derived condition.

thumbnailFigure 4. Plots of principal component (PC) and centroid size (CS) for morphometric characters of Polypterus. a) CS and PC1 for dorsal view, b) PC1 and PC 2 for dorsal view, c) CS and PC1 for ventral view, d) PC1 and PC2 for ventral view. Grouping corresponds to the clade inferred from molecular phylogenetic analysis (Fig. 2): Blue = Polypterus ansorgei, P. bichir lapradei, P. endlicheri endlicheri, and P. e. congicus; Green = P. mokelembembe; Purple = P. ornatipinnis, P. retropinnis, and P. weeksii; Red = P. delhezi, P. palmas buettikoferi, P. p. polli, P. senegalus senegalus, and P. teugelsi.

Interestingly, P. delhezi displayed a similar tendency of morphological change in PC1 axis of both dorsal and lateral views, indicating that this species actually possesses an intermediate condition where neither the upper or lower jaw show pronounced rostral projection. Because P. delhezi is distantly related to the lower jaw projection clade (P. ansorgii, P. bichir and P. endlicheri), the similar cranial shape seen in two polypterid lineages has apparently evolved independently. In polypterid species whose snout projected more rostrally to the tip of mandible, the range of variation was broader than that in species with lower jaw protrusion (Fig. 4). Furthermore, the range was broadly overlapped between two main clades: the monophyletic group composed of P. ornatipinnis, P. weeksii, and P. retropinnis (purple) and that composed of P. delhezi, P. palmas, P. senegalus, and P. teugelsi (red), showing no lineage-specific pattern of craniofacial transformation.

There are multiple potentially productive avenues of future genetic and developmental research into the mechanisms of craniofacial development in Polypteridae. In vertebrates, morphology of the lower jaw is regulated by cranial neural crest cells, and several transcription factors in BMP and/or FGF signaling pathway appear to be involved in the process [4,38,93-97]. Although we do not know what change in developmental program caused substantial transformation of craniofacial morphology in P. ansorgii, P. bichir and P. endlicheri, the event that occurred in relatively early phase of polypterid evolution might have facilitate the shift into a new ecological niche, affecting their feeding strategy.

Conclusion

The discovery of patterns of morphological evolution among "deep" lineages of organisms permits further exploration of the general underlying genetic and developmental mechanisms that evolved early in evolutionary history. Our new molecular phylogeny of Polypteridae, in conjunction with ancestral state reconstruction and geometric morphometric analyses reveals that many of the patterns of morphological evolution seen in polypterids are shared with other actinopterygian, and indeed, osteichthyan lineages. Thus, this implies the underlying genetic and developmental mechanisms responsible for those general patterns were established early in the evolutionary history of Osteichthyes, or perhaps even older in the craniate tree of life. Moreover, these patterns suggest interesting avenues of future evolution and developmental research that are only possible in a phylogenetic framework.

Authors' contributions

MT designed the study. DS and MT collected data. DS and MT performed the morphometric analyses and MCB performed the phylogenetic analyses. MT wrote the paper with the assistance of DS and MCB. All authors read and approved the final manuscript.

Acknowledgements

We are grateful to Nori Satoh, Michio Hori, Naomichi Ogihara, Takushi Kishida, Wataru Yano, and Nobuto Kikukawa for technical support. We also thank Hiroshi Senou (Kanagawa Prefectural Museum of Natural History) and Patrice Pruvost (Muséum d'Histoire naturelle de Paris) for access to museum collections; Kenko Takemura for providing fish specimens; Toshiaki Igarashi for access to literature; Tsutomu Hikida, Thomas Near, Ainsley Seago, Helen Whitaker and two anonymous reviewers for providing valuable comments on the manuscript; Ainsley Seago and Jun Inoue for the illustrations in Fig. 1; Aquarium Shop Water House (Yokohama, Japan) for the pictures in Fig. 2; and Jim McGuire and the Museum of Vertebrate Zoology for providing computing facilities. This study was financially supported by a Grant-in-Aid for the biodiversity research of the 21st Century COE (A14) and a Grant-in-Aid for the biodiversity and evolutionary research of the Global COE (A06) from the Ministry of Education, Culture, Sports and Science and Technology, Japan to Kyoto University.

References

  1. Benton MJ: Vertebrate Paleontology. 3rd edition. Malden, MA: Blackwell Science; 2004.

  2. Botella H, Blom H, Dorka M, Ahlberg PE, Janvier P: Jaws and teeth of the earliest bony fishes.

    Nature 2007, 448:583-586. PubMed Abstract | Publisher Full Text OpenURL

  3. Shapiro MD, Hanken J, Rosenthal N: Developmental basis of evolutionary digit loss in the Australian lizard Hemiergis.

    J Exp Zool 2003, 297:48-56. Publisher Full Text OpenURL

  4. Albertson RC, Streelman JT, Kocher TD, Yelick PC: Integration and evolution of the cichlid mandible: the molecular basis of alternate feeding strategies.

    Proc Natl Acad Sci USA 2005, 102:16287-16292. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ: The calmodulin pathway and evolution of elongated beak morphology in Darwin's finches.

    Nature 2006, 442:563-567. PubMed Abstract | Publisher Full Text OpenURL

  6. Hoekstra HE, Hirschmann RJ, Bundey RA, Insel PA, Crossland JP: A single amino acid mutation contributes to adaptive beach mouse color pattern.

    Science 2006, 313:101-104. PubMed Abstract | Publisher Full Text OpenURL

  7. Steiner CC, Weber JN, Hoekstra HE: Adaptive variation in beach mice produced by two interacting pigmentation genes.

    PLoS Biol 2007, 5:e219. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Vonk FJ, Admiraal JF, Jackson K, Reshef R, de Bakker MA, Vanderschoot K, Berge I, van Atten M, Burgerhout E, Beck A, Mirtschin PJ, Kochva E, Witte F, Fry BG, Woods AE, Richardson MK: Evolutionary origin and development of snake fangs.

    Nature 2008, 454:630-633. PubMed Abstract | Publisher Full Text OpenURL

  9. Gross JB, Borowsky R, Tabin CJ: A Novel Role for Mc1r in the parallel evolution of depigmentation in independent populations of the cavefish Astyanax mexicanus.

    PLoS Genet 2009, 5:e1000326. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Nelson JS: Fishes of the World. 4th edition. New York: Wiley; 2006.

  11. Britz R, Johnson GD: On the homology of the posteriormost gill arch in polypterids (Cladistia, Actinopterygii).

    Zool J Linn Soc 2003, 138:495-503. Publisher Full Text OpenURL

  12. Claeson KM, Bemis WE, Hagadorn JW: New interpretations of the skull of a primitive bony fish Erpetoichthys calabaricus (Actinopterygii: Cladistia).

    J Morphol 2007, 268:1021-1039. PubMed Abstract | Publisher Full Text OpenURL

  13. Takeuchi M, Okabe M, Aizawa S: The Genus Polypterus (Bichirs). A fish group diverged at the stem of ray-finned fishes (Actinopterygii). In Emerging Model Organisms: A Laboratory Manual. Volume 1. New York: Cold Spring Harbor Laboratory Press; 2009::447-467. OpenURL

  14. Patterson C: Morphology and interrelationships of primitive actinopterygian fishes.

    Am Zool 1982, 22:241-259. OpenURL

  15. Gardiner BG, Schaeffer B: Interrelationships of lower actinopterygian fishes.

    Zool J Linn Soc 1989, 97:135-187. Publisher Full Text OpenURL

  16. Noack K, Zardoya R, Meyer A: The complete mitochondrial DNA sequence of the bichir (Polypterus ornatipinnis), a basal ray-finned fish: Ancient establishment of the consensus vetebrate gene order.

    Genetics 1996, 144:1165-1180. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Venkatesh B, Erdmann MV, Brenner S: Molecular synapomorphies resolve evolutionary relationships of extant jawed vertebrates.

    Proc Natl Acad Sci USA 2001, 98:11382-11387. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Inoue JG, Miya M, Tsukamoto K, Nishida M: Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish.".

    Mol Phylogenet Evol 2003, 26:110-120. PubMed Abstract | Publisher Full Text OpenURL

  19. Kikugawa K, Katoh K, Kuraku S, Sakurai H, Ishida O, Iwabe N, Miyata T: Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes.

    BMC Evol Biol 2004, 2:3. OpenURL

  20. Agnese JF, Teugels GG: Insight into the phylogeny of African Clariidae (Teleostei, Siluriformes): Implications for their body shape evolution, biogeography, and taxonomy.

    Mol Phylogenet Evol 2005, 36:546-553. PubMed Abstract | Publisher Full Text OpenURL

  21. Devaere S, Jansen G, Adriaens D, Weekers P: Phylogeny of the African representatives of the catfish family Clariidae (Teleostei, Siluriformes) based on a combined analysis: independent evolution towards anguilliformity.

    J Zool Syst Evol Res 2007, 45:214-229. Publisher Full Text OpenURL

  22. Ward AB, Brainerd EL: Evolution of axial patterning in elongate fishes.

    Biol J Linn Soc 2007, 90:97-116. Publisher Full Text OpenURL

  23. Yamada T, Sugiyama T, Tamaki N, Kawakita A, Kato M: Adaptive radiation of gobies in the interstitial habitats of gravel beaches accompanied by body elongation and excessive vertebral segmentation.

    BMC Evol Biol 2009, 9:145. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Wiens JJ, Hoverman JT: Digit reduction, body size, and paedomorphosis in salamanders.

    Evol Dev 2008, 10:449-463. PubMed Abstract | Publisher Full Text OpenURL

  25. Wiens JJ, Slingluff JL: How lizards turn into snakes: a phylogenetic analysis of body-form evolution in anguid lizards.

    Evolution 2001, 55:2303-2318. PubMed Abstract OpenURL

  26. Wiens JJ, Brandley MC, Reeder TW: Why does a trait evolve multiple times within a clade? Repeated evolution of snake-like body form in squamate reptiles.

    Evolution 2006, 60:123-141. PubMed Abstract OpenURL

  27. Brandley MC, Huelsenbeck JP, Wiens JJ: Rates and patterns in the evolution of snake-like body form in squamate reptiles: evidence for repeated re-evolution of lost digits and long-term persistence of intermediate body forms.

    Evolution 2008, 62:2042-2064. PubMed Abstract | Publisher Full Text OpenURL

  28. Carroll RL: The origin and early radiation of terrestrial vertebrates.

    J Paleontol 2001, 75:1202-1213. Publisher Full Text OpenURL

  29. Palci A, Caldwell MW: Vestigial forelimbs and axial elongation in a 95 million-year old non-snake squamate.

    J Verteb Paleontol 2007, 27:1-7. Publisher Full Text OpenURL

  30. Adriaens D, Devaere S, Teugels GG, Dekegel B, Verraes W: Intraspecific variation in limblessness in vertebrates: a unique example of microevolution.

    Biol J Linn Soc 2002, 75:367-377. Publisher Full Text OpenURL

  31. Shapiro MD, Marks ME, Peichel CL, Blackman BK, Nereng KS, Jónsson B, Schluter D, Kingsley DM: Genetic and developmental basis of evolutionary pelvic reduction in threespine sticklebacks.

    Nature 2004, 428:717-23. PubMed Abstract | Publisher Full Text OpenURL

  32. Tanaka M, Hale LA, Amores A, Yan YL, Cresko WA, Suzuki T, Postlethwait JH: Developmental genetic basis for the evolution of pelvic fin loss in the pufferfish Takifugu rubripes.

    Dev Biol 2005, 15:227-39. Publisher Full Text OpenURL

  33. Schaefer SA, Lauder GV: Testing historical hypotheses of morphological change: biomechanical decoupling in Loricarioid catfishes.

    Evolution 1996, 50:1661-1675. Publisher Full Text OpenURL

  34. Albertson RC, Kocher TD: Assessing morphological differences in an adaptive trait: a landmark-based morphometric approach.

    J Exp Zool 2001, 289:385-403. PubMed Abstract | Publisher Full Text OpenURL

  35. Okada N, Takagi Y, Seikai T, Tanaka M, Tagawa M: Asymmetrical development of bones and soft tissues during eye migration of metamorphosing Japanese flounder, Paralichthys olivaceus.

    Cell Tissue Res 2001, 304:59-66. PubMed Abstract | Publisher Full Text OpenURL

  36. Kijimoto T, Watanabe M, Fujimura K, Nakazawa M, Murakami Y, Kuratani S, Kohara Y, Gojobori T, Okada N: cimp1, a novel astacin family metalloproteinase gene from East African Cichlids, is differentially expressed between species during growth.

    Mol Biol Evol 2005, 22:1649-1660. PubMed Abstract | Publisher Full Text OpenURL

  37. Westneat MW, Alfaro ME, Wainwright PC, Bellwood D, Grubich JR, Fessler JL, Clements KD, Smith LL: Local phylogenetic divergence and global evolutionary convergence of skull function in reef fishes of the family Labridae.

    Proc R Soc B 2005, 272:993-1000. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Albertson RC, Kocher TD: Genetic and developmental basis of cichlid trophic diversity.

    Heredity 2006, 97:211-221. PubMed Abstract | Publisher Full Text OpenURL

  39. Francis AW Jr, Turingan RG: Morphological and biomechanical changes of the feeding apparatus in developing southern flounder, Paralichthys lethostigma.

    J Morphol 2008, 269:1169-1180. PubMed Abstract | Publisher Full Text OpenURL

  40. Sidlauskas B: Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach.

    Evolution 2008, 62:3135-3156. PubMed Abstract | Publisher Full Text OpenURL

  41. Cooper WJ, Westneat MW: Form and function of damselfish skulls: rapid and repeated evolution into a limited number of trophic niches.

    BMC Evol Biol 2009, 9:24. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  42. Chiu C, Dewar K, Wagner GP, Takahashi K, Ruddle F, Ledje C, Bartsch P, Scemama J-L, Stellwag E, Fried C, Prohaska SJ, Stadler PF, Amemiya CT: Bichir HoxA cluster sequence reveals surprising trends in ray-finned fish genomic evolution.

    Genome Res 2004, 14:11-17. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Boulenger G: Catalogue of the freshwater fishes of Africa in the British Museum (Natural History). London: Wheldon and Wesley; 1909.

  44. Poll M: Contribution â l'étude systématique des Polyptèridae (première partie).

    Revue de Zoologie et de Botanique Africaines 1941, 35:141-179. OpenURL

  45. Poll M: Les tendances évolutives des Polyptères d'après l'étude systématique des espèces.

    Annales de la Société Royale Zoologique de Belgique 1941, 72:157-173. OpenURL

  46. Poll M: Contribution â l'étude systématique des Polyptèridae (suite et fin).

    Revue de Zoologie et de Botanique Africaines 1942, 35:269-308. OpenURL

  47. Daget J, Desoutter M: Essai de classification cladistique des Polyptéridés (Pisces, Branchiopterygii).

    Bull Mus Natn Hist Nat Paris 1983, 4:661-674. OpenURL

  48. Tokita M, Okamoto T, Hikida T: Evolutionary history of African lungfish: a hypothesis from molecular phylogeny.

    Mol Phylogenet Evol 2005, 35:281-286. PubMed Abstract | Publisher Full Text OpenURL

  49. Castoe TA, Doan TM, Parkinson CL: Data partitions and complex models in Bayesian analysis: the phylogeny of gymnophthalmid lizards.

    Syst Biol 2004, 53:448-469. PubMed Abstract | Publisher Full Text OpenURL

  50. Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data.

    Syst Biol 2004, 5:47-67. Publisher Full Text OpenURL

  51. Brandley MC, Schmitz A, Reeder TW: Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards.

    Syst Biol 2005, 54:373-390. PubMed Abstract | Publisher Full Text OpenURL

  52. Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4.0b10. Sunderland, MA: Sinauer Associates; 2002. PubMed Abstract | Publisher Full Text OpenURL

  53. Schwarz G: Estimating the dimension of a model.

    Ann Statis 1978, 6:461-464. Publisher Full Text OpenURL

  54. Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models.

    Bioinformatics 2003, 19:1572-1574. PubMed Abstract | Publisher Full Text OpenURL

  55. Altekar G, Dwarkadas S, Huelsenbeck JP, Ronquist F: Parallel Metropolis-coupled Markov chain Monte Carlo for Bayesian phylogenetic inference.

    Bioinformatics 2004, 20:407-415. PubMed Abstract | Publisher Full Text OpenURL

  56. Nylander JAA, Wilgenbusch JC, Warren DL, Swofford DL: AWTY (are we there yet?): a system for graphical exploration of MCMC convergence in Bayesian phylogenetics.

    Bioinformatics 2008, 24:581-583. PubMed Abstract | Publisher Full Text OpenURL

  57. Huelsenbeck JP, Rannala B: Frequentist properties of Bayesian posterior probabilities of phylogenetic trees under simple and complex substitution models.

    Syst Biol 2004, 53:904-913. PubMed Abstract | Publisher Full Text OpenURL

  58. Marcus LF, Corti M, Loy A, Naylor GJP, Slice DE: Advances in Morphometrics, NATO ASI Series A: Life Sciences, vol. 284. New York: Plenum Press; 1996.

  59. Zeldich ML, Swiderski D, Sheets DH, Fink W: Geometric Morphometrics for Biologists: a Primer. Elsevier Academic Press, Amsterdam; 2004:443.

  60. Rohlf FJ: tpsDig2.10. [http://life.bio.sunysb.edu/morph/] webcite

    2006. PubMed Abstract OpenURL

  61. Sheets HD: Coordgen6 h in Integrated Morphometrics Package (IMP). [http://www3.canisius.edu/~sheets/moremorph.html] webcite

    Department of Physics, Canisius College; 2006.

  62. Sheets HD: Pcagen6 h in Integrated Morphometrics Package (IMP). [http://www3.canisius.edu/~sheets/moremorph.html] webcite

    Department of Physics, Canisius College; 2007.

  63. Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.6. [http://mesquiteproject.org] webcite

    2009. PubMed Abstract | Publisher Full Text OpenURL

  64. Froese R, Pauly D: FishBase. version (09/2008). [http://www.fishbase.org/search.php] webcite

    2008. PubMed Abstract | Publisher Full Text OpenURL

  65. Schmitz A, Brandley MC, Mausfeld P, Vences M, Glaw F, Nussbaum RA, Reeder TW: Opening the black box: phylogenetics and morphological evolution of the Malagasy fossorial lizards of the subfamily "Scincinae".

    Mol Phylogenet Evol 2005, 34:118-133. PubMed Abstract | Publisher Full Text OpenURL

  66. Maddison WP: Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree.

    Syst Zool 1991, 40:304-314. Publisher Full Text OpenURL

  67. Schluter D, Price T, Mooers AO, Ludwig D: Likelihood of ancestor states in adaptive radiation.

    Evolution 1997, 51:1699-1711. Publisher Full Text OpenURL

  68. Pagel M: The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies.

    Syst Biol 1999, 48:612-622. Publisher Full Text OpenURL

  69. Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data.

    Syst Biol 2001, 50:913-925. PubMed Abstract | Publisher Full Text OpenURL

  70. Myers GS: Fresh-water fishes and West Indian zoogeography.

    Annual Report of the Board of Regents of the Smithsonian Institution 1938, 92:339-364. OpenURL

  71. Myers GS: Salt-tolerance of fresh-water fish groups in relation to zoogeographical problems.

    Bijdragen tot de Dierkunde 1949, 28:315-322. OpenURL

  72. Gosse JP: Polypteridae.

    In Faune des poisons d'eaux douces et saumâtres d'Afrique de l'Ouest I. ORSTOM 28 Edited by Lévèque C, Paugy D, Teugels GG. 1990, 79-87. OpenURL

  73. Hanssens MM, Teugels GG, Audenaerde DFE: Subspecies in the Polypterus palmas complex (Brachiopterygii; Polypteridae) from West and Central Africa.

    Copeia 1995, 1995:694-705. Publisher Full Text OpenURL

  74. Schliewen UK, Schäfer F: Polypterus mokelembembe, a new species of bichir from the central Congo River basin (Actinopterygii: Cladistia: Polypteridae).

    Zootaxa 2006, 1129:23-36. OpenURL

  75. Bossu CM, Near TJ: Gene trees reveal repeated instances of mitochondrial DNA introgression in orangethroat darters (Percidae:Etheostoma).

    Syst Biol 2009, 58:114-129. Publisher Full Text OpenURL

  76. Funk DJ, Omland KE: Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA.

    Annu Rev Ecol Syst 2003, 34:397-423. Publisher Full Text OpenURL

  77. Richardson MK, Allen SP, Wright GM, Raynaud A, Hanken J: Somite number and vertebrate evolution.

    Development 1998, 125:151-160. PubMed Abstract | Publisher Full Text OpenURL

  78. Özbudak EM, Pourquié O: The vertebrate segmentation clock: the tip of the iceburg.

    Curr Opin Genet Dev 2008, 18:317-323. PubMed Abstract | Publisher Full Text OpenURL

  79. Vonk FJ, Richardson MK: Serpent clocks tick faster.

    Nature 2008, 454:282-283. PubMed Abstract | Publisher Full Text OpenURL

  80. Gomez C, Özbudak EM, Wunderlich J, Baumann D, Lewis J, Pourquié O: Control of segment number in vertebrate embryos.

    Nature 2008, 454:335-339. PubMed Abstract | Publisher Full Text OpenURL

  81. Clausen HS: Larvae of polypterid fish Erpetoichthys Smith.

    Nature 1956, 178:932-933. Publisher Full Text OpenURL

  82. Britz R, Bartsch P: On the reproduction and early development of Erpetoichthys calabaricus, Polypterus senegalus, and Polypterus ornatipinnis (Actinopterygii: Polypteridae).

    Ichthyol Expl Freshwaters 1998, 9:325-334. OpenURL

  83. Gomez C, Pourquié O: Developmental control of segment numbers in vertebrates.

    J Exp Zool B Mol Dev Evol 2009, 312:533-544. PubMed Abstract | Publisher Full Text OpenURL

  84. Gans C: Tetrapod limblessness: evolution and functional corollaries.

    Am Zool 1975, 15:455-467. OpenURL

  85. Lande R: Evolutionary mechanisms of limb loss in tetrapods.

    Evolution 1978, 32:73-92. Publisher Full Text OpenURL

  86. Cohn MJ, Tickle C: Developmental basis of limblessness and axial patterning in snakes.

    Nature 1999, 399:474-479. PubMed Abstract | Publisher Full Text OpenURL

  87. Sanger TJ, Gibson-Brown JJ: The developmental bases of limb reduction and body elongation in squamates.

    Evolution 2004, 58:2103-2106. PubMed Abstract OpenURL

  88. Woltering JM, Vonk FJ, Müller H, Bardine N, Tuduce IL, de Bakker MAG, Knöchel WP, Sirbu IO, Durston AJ, Richardson MK: Axial patterning in snakes and caecilians: evidence for an alternative interpretation of the Hox code.

    Dev Biol 2009, 332:82-89. PubMed Abstract | Publisher Full Text OpenURL

  89. Rahmani TM-Z: Morphogenesis of the rudimentary hind-limb of the glass snake (Ophisaurus apodus Pallas).

    J Embryol Exp Morphol 1974, 32:431-443. PubMed Abstract OpenURL

  90. Raynaud A, Kan P, Bouche G, Duprat AM: Fibroblast growth factors (FGF-2) and delayed involution of the posterior limbs of the slow-worm embryo (Anguis fragilis, L.).

    Comptes rendus de l'Académie des sciences. Série III, Sciences de la vie 1995, 318:573-578. OpenURL

  91. Thewissen JGM, Cohn MJ, Stevens LS, Bajpai S, Heyning J, Horton WE Jr: Developmental basis for hind-limb loss in dolphins and origin of the cetacean bodyplan.

    Proc Natl Acad Sci USA 2006, 103:8414-8418. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  92. Shapiro MD, Bell MA, Kingsley DM: Parallel genetic origins of pelvic reduction in vertebrates.

    Proc Natl Acad Sci USA 2006, 103:13753-13758. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  93. Stottmann RW, Anderson RM, Klingensmith J: The BMP antagonists Chordin and Noggin have essential but redundant roles in mouse mandibular outgrowth.

    Dev Biol 2001, 240:457-473. PubMed Abstract | Publisher Full Text OpenURL

  94. Hide T, Hatakeyama J, Kimura-Yoshida C, Tian E, Takeda N, Ushio Y, Shiroishi T, Aizawa S, Matsuo I: Genetic modifiers of otocephalic phenotypes in Otx2 heterozygous mutant mice.

    Development 2002, 129:4347-4357. PubMed Abstract | Publisher Full Text OpenURL

  95. Wilson J, Tucker AS: Fgf and Bmp signals repress the expression of Bapx1 in the mandibular mesenchyme and control the position of the developing jaw joint.

    Dev Biol 2004, 266:138-150. PubMed Abstract | Publisher Full Text OpenURL

  96. Eames BF, Schneider RA: The genesis of cartilage size and shape during development and evolution.

    Development 2008, 135:3947-3958. PubMed Abstract | Publisher Full Text OpenURL

  97. Sperber SM, Dawid IB: barx1 is necessary for ectomesenchyme proliferation and osteochondroprogenitor condensation in the zebrafish pharyngeal arches.

    Dev Biol 2008, 321:101-110. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL