Skip to main content

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



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.


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.


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.


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., [39]].

Figure 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 [1113]. However, recent morphological and molecular analyses have determined that bichirs are a basal lineage of Actinopterygii (ray-finned fishes; Fig. 1) [1419].

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., [2023]], amphibians [e.g., [24]], squamate reptiles [e.g., [2527]], 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, 3032]. Furthermore, diverse changes in craniofacial morphology have been the subject of particularly extensive research [5, 3341].

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 [4446] 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.


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).

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.

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 [4951]. 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].

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).

Figure 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, 7072].

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 [4446] 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.

Figure 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, 2327, 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 [2426].

Pelvic fin loss

In vertebrates, the body elongation tends to co-occur with a loss or reduction of paired appendages [21, 2426, 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.

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.

Figure 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, 9397]. 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.


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.


  1. 1.

    Benton MJ: Vertebrate Paleontology. 2004, Malden, MA: Blackwell Science, 3

    Google Scholar 

  2. 2.

    Botella H, Blom H, Dorka M, Ahlberg PE, Janvier P: Jaws and teeth of the earliest bony fishes. Nature. 2007, 448: 583-586. 10.1038/nature05989.

    Article  CAS  PubMed  Google Scholar 

  3. 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. 10.1002/jez.b.19.

    Article  Google Scholar 

  4. 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. 10.1073/pnas.0506649102.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  5. 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. 10.1038/nature04843.

    Article  CAS  PubMed  Google Scholar 

  6. 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. 10.1126/science.1126121.

    Article  CAS  PubMed  Google Scholar 

  7. 7.

    Steiner CC, Weber JN, Hoekstra HE: Adaptive variation in beach mice produced by two interacting pigmentation genes. PLoS Biol. 2007, 5: e219-10.1371/journal.pbio.0050219.

    PubMed Central  Article  PubMed  Google Scholar 

  8. 8.

    Vonk FJ, Admiraal JF, Jackson K, Reshef R, de Bakker MA, Vanderschoot K, Berge van den 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. 10.1038/nature07178.

    Article  CAS  PubMed  Google Scholar 

  9. 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-10.1371/journal.pgen.1000326.

    PubMed Central  Article  PubMed  Google Scholar 

  10. 10.

    Nelson JS: Fishes of the World. 2006, New York: Wiley, 4

    Google Scholar 

  11. 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. 10.1046/j.1096-3642.2003.t01-1-00067.x.

    Article  Google Scholar 

  12. 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. 10.1002/jmor.10567.

    Article  PubMed  Google Scholar 

  13. 13.

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

    Google Scholar 

  14. 14.

    Patterson C: Morphology and interrelationships of primitive actinopterygian fishes. Am Zool. 1982, 22: 241-259.

    Article  Google Scholar 

  15. 15.

    Gardiner BG, Schaeffer B: Interrelationships of lower actinopterygian fishes. Zool J Linn Soc. 1989, 97: 135-187. 10.1111/j.1096-3642.1989.tb00550.x.

    Article  Google Scholar 

  16. 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 Central  CAS  PubMed  Google Scholar 

  17. 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. 10.1073/pnas.201415598.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  18. 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. 10.1016/S1055-7903(02)00331-7.

    Article  CAS  PubMed  Google Scholar 

  19. 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-

    Article  Google Scholar 

  20. 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. 10.1016/j.ympev.2005.03.028.

    Article  PubMed  Google Scholar 

  21. 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. 10.1111/j.1439-0469.2006.00399.x.

    Article  Google Scholar 

  22. 22.

    Ward AB, Brainerd EL: Evolution of axial patterning in elongate fishes. Biol J Linn Soc. 2007, 90: 97-116. 10.1111/j.1095-8312.2007.00714.x.

    Article  Google Scholar 

  23. 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-10.1186/1471-2148-9-145.

    PubMed Central  Article  PubMed  Google Scholar 

  24. 24.

    Wiens JJ, Hoverman JT: Digit reduction, body size, and paedomorphosis in salamanders. Evol Dev. 2008, 10: 449-463. 10.1111/j.1525-142X.2008.00256.x.

    Article  PubMed  Google Scholar 

  25. 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.

    Article  CAS  PubMed  Google Scholar 

  26. 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  Google Scholar 

  27. 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. 10.1111/j.1558-5646.2008.00430.x.

    Article  PubMed  Google Scholar 

  28. 28.

    Carroll RL: The origin and early radiation of terrestrial vertebrates. J Paleontol. 2001, 75: 1202-1213. 10.1666/0022-3360(2001)075<1202:TOAERO>2.0.CO;2.

    Article  Google Scholar 

  29. 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. 10.1671/0272-4634(2007)27[1:VFAAEI]2.0.CO;2.

    Article  Google Scholar 

  30. 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. 10.1111/j.1095-8312.2002.tb02078.x.

    Article  Google Scholar 

  31. 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. 10.1038/nature02415.

    Article  CAS  PubMed  Google Scholar 

  32. 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. 10.1016/j.ydbio.2005.02.016.

    Article  Google Scholar 

  33. 33.

    Schaefer SA, Lauder GV: Testing historical hypotheses of morphological change: biomechanical decoupling in Loricarioid catfishes. Evolution. 1996, 50: 1661-1675. 10.2307/2410902.

    Article  Google Scholar 

  34. 34.

    Albertson RC, Kocher TD: Assessing morphological differences in an adaptive trait: a landmark-based morphometric approach. J Exp Zool. 2001, 289: 385-403. 10.1002/jez.1020.

    Article  CAS  PubMed  Google Scholar 

  35. 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. 10.1007/s004410100353.

    Article  CAS  PubMed  Google Scholar 

  36. 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. 10.1093/molbev/msi159.

    Article  CAS  PubMed  Google Scholar 

  37. 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. 10.1098/rspb.2004.3013.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  38. 38.

    Albertson RC, Kocher TD: Genetic and developmental basis of cichlid trophic diversity. Heredity. 2006, 97: 211-221. 10.1038/sj.hdy.6800864.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Francis AW, Turingan RG: Morphological and biomechanical changes of the feeding apparatus in developing southern flounder, Paralichthys lethostigma. J Morphol. 2008, 269: 1169-1180. 10.1002/jmor.10631.

    Article  PubMed  Google Scholar 

  40. 40.

    Sidlauskas B: Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach. Evolution. 2008, 62: 3135-3156. 10.1111/j.1558-5646.2008.00519.x.

    Article  PubMed  Google Scholar 

  41. 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-10.1186/1471-2148-9-24.

    PubMed Central  Article  PubMed  Google Scholar 

  42. 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. 10.1101/gr.1712904.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  43. 43.

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

    Google Scholar 

  44. 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.

    Google Scholar 

  45. 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.

    Google Scholar 

  46. 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.

    Google Scholar 

  47. 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.

    Google Scholar 

  48. 48.

    Tokita M, Okamoto T, Hikida T: Evolutionary history of African lungfish: a hypothesis from molecular phylogeny. Mol Phylogenet Evol. 2005, 35: 281-286. 10.1016/j.ympev.2004.11.025.

    Article  CAS  PubMed  Google Scholar 

  49. 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. 10.1080/10635150490445797.

    Article  PubMed  Google Scholar 

  50. 50.

    Nylander JAA, Ronquist F, Huelsenbeck JP, Nieves-Aldrey JL: Bayesian phylogenetic analysis of combined data. Syst Biol. 2004, 5: 47-67. 10.1080/10635150490264699.

    Article  Google Scholar 

  51. 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. 10.1080/10635150590946808.

    Article  PubMed  Google Scholar 

  52. 52.

    Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4.0b10. 2002, Sunderland, MA: Sinauer Associates

    Google Scholar 

  53. 53.

    Schwarz G: Estimating the dimension of a model. Ann Statis. 1978, 6: 461-464. 10.1214/aos/1176344136.

    Article  Google Scholar 

  54. 54.

    Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.

    Article  CAS  PubMed  Google Scholar 

  55. 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. 10.1093/bioinformatics/btg427.

    Article  CAS  PubMed  Google Scholar 

  56. 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. 10.1093/bioinformatics/btm388.

    Article  CAS  PubMed  Google Scholar 

  57. 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. 10.1080/10635150490522629.

    Article  PubMed  Google Scholar 

  58. 58.

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

    Chapter  Google Scholar 

  59. 59.

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

    Google Scholar 

  60. 60.

    Rohlf FJ: tpsDig2.10. 2006, []

    Google Scholar 

  61. 61.

    Sheets HD: Coordgen6 h in Integrated Morphometrics Package (IMP). 2006, Department of Physics, Canisius College, []

    Google Scholar 

  62. 62.

    Sheets HD: Pcagen6 h in Integrated Morphometrics Package (IMP). 2007, Department of Physics, Canisius College, []

    Google Scholar 

  63. 63.

    Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.6. 2009, []

    Google Scholar 

  64. 64.

    Froese R, Pauly D: FishBase. version (09/2008). 2008, []

    Google Scholar 

  65. 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. 10.1016/j.ympev.2004.08.016.

    Article  CAS  PubMed  Google Scholar 

  66. 66.

    Maddison WP: Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree. Syst Zool. 1991, 40: 304-314. 10.2307/2992324.

    Article  Google Scholar 

  67. 67.

    Schluter D, Price T, Mooers AO, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51: 1699-1711. 10.2307/2410994.

    Article  Google Scholar 

  68. 68.

    Pagel M: The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies. Syst Biol. 1999, 48: 612-622. 10.1080/106351599260184.

    Article  Google Scholar 

  69. 69.

    Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50: 913-925. 10.1080/106351501753462876.

    Article  CAS  PubMed  Google Scholar 

  70. 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.

    Google Scholar 

  71. 71.

    Myers GS: Salt-tolerance of fresh-water fish groups in relation to zoogeographical problems. Bijdragen tot de Dierkunde. 1949, 28: 315-322.

    Google Scholar 

  72. 72.

    Gosse JP: Polypteridae. 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.

    Google Scholar 

  73. 73.

    Hanssens MM, Teugels GG, Audenaerde Thys van den DFE: Subspecies in the Polypterus palmas complex (Brachiopterygii; Polypteridae) from West and Central Africa. Copeia. 1995, 1995: 694-705. 10.2307/1446767.

    Article  Google Scholar 

  74. 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.

    Google Scholar 

  75. 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. 10.1093/sysbio/syp014.

    Article  CAS  PubMed  Google Scholar 

  76. 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. 10.1146/annurev.ecolsys.34.011802.132421.

    Article  Google Scholar 

  77. 77.

    Richardson MK, Allen SP, Wright GM, Raynaud A, Hanken J: Somite number and vertebrate evolution. Development. 1998, 125: 151-160.

    CAS  PubMed  Google Scholar 

  78. 78.

    Özbudak EM, Pourquié O: The vertebrate segmentation clock: the tip of the iceburg. Curr Opin Genet Dev. 2008, 18: 317-323. 10.1016/j.gde.2008.06.007.

    Article  PubMed  Google Scholar 

  79. 79.

    Vonk FJ, Richardson MK: Serpent clocks tick faster. Nature. 2008, 454: 282-283. 10.1038/454282a.

    Article  CAS  PubMed  Google Scholar 

  80. 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. 10.1038/nature07020.

    Article  CAS  PubMed  Google Scholar 

  81. 81.

    Clausen HS: Larvae of polypterid fish Erpetoichthys Smith. Nature. 1956, 178: 932-933. 10.1038/178932b0.

    Article  Google Scholar 

  82. 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.

    Google Scholar 

  83. 83.

    Gomez C, Pourquié O: Developmental control of segment numbers in vertebrates. J Exp Zool B Mol Dev Evol. 2009, 312: 533-544. 10.1002/jez.b.21305.

    PubMed Central  Article  PubMed  Google Scholar 

  84. 84.

    Gans C: Tetrapod limblessness: evolution and functional corollaries. Am Zool. 1975, 15: 455-467.

    Article  Google Scholar 

  85. 85.

    Lande R: Evolutionary mechanisms of limb loss in tetrapods. Evolution. 1978, 32: 73-92. 10.2307/2407411.

    Article  Google Scholar 

  86. 86.

    Cohn MJ, Tickle C: Developmental basis of limblessness and axial patterning in snakes. Nature. 1999, 399: 474-479. 10.1038/20944.

    Article  CAS  PubMed  Google Scholar 

  87. 87.

    Sanger TJ, Gibson-Brown JJ: The developmental bases of limb reduction and body elongation in squamates. Evolution. 2004, 58: 2103-2106.

    Article  PubMed  Google Scholar 

  88. 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. 10.1016/j.ydbio.2009.04.031.

    Article  CAS  PubMed  Google Scholar 

  89. 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.

    CAS  PubMed  Google Scholar 

  90. 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.

    CAS  PubMed  Google Scholar 

  91. 91.

    Thewissen JGM, Cohn MJ, Stevens LS, Bajpai S, Heyning J, Horton WE: Developmental basis for hind-limb loss in dolphins and origin of the cetacean bodyplan. Proc Natl Acad Sci USA. 2006, 103: 8414-8418. 10.1073/pnas.0602920103.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  92. 92.

    Shapiro MD, Bell MA, Kingsley DM: Parallel genetic origins of pelvic reduction in vertebrates. Proc Natl Acad Sci USA. 2006, 103: 13753-13758. 10.1073/pnas.0604706103.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  93. 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. 10.1006/dbio.2001.0479.

    Article  CAS  PubMed  Google Scholar 

  94. 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.

    CAS  PubMed  Google Scholar 

  95. 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. 10.1016/j.ydbio.2003.10.012.

    Article  CAS  PubMed  Google Scholar 

  96. 96.

    Eames BF, Schneider RA: The genesis of cartilage size and shape during development and evolution. Development. 2008, 135: 3947-3958. 10.1242/dev.023309.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  97. 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. 10.1016/j.ydbio.2008.06.004.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

Download references


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.

Author information



Corresponding author

Correspondence to Masayoshi Tokita.

Additional information

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.

Electronic supplementary material

Additional file 1: Information of DNA sequences analyzed in this study. (PDF 586 KB)

Additional file 2: List of primers used in this study. (PDF 298 KB)

Landmark definition for morphometric analysis

Additional file 3: . (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. (PDF 638 KB)

Thin-plate splines (TPS) of principal component 1 (PC1) and 2 (PC2) for both dorsal and ventral views of the head of Polypterus

Additional file 4: . 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. (PDF 129 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Suzuki, D., Brandley, M.C. & Tokita, M. 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. BMC Evol Biol 10, 21 (2010).

Download citation


  • Much Recent Common Ancestor
  • Phylogenetic Framework
  • Vertebral Number
  • Ancestral State Reconstruction
  • Craniofacial Morphology