Independent evolution of the specialized pharyngeal jaw apparatus in cichlid and labrid fishes
© Mabuchi et al; licensee BioMed Central Ltd. 2007
Received: 01 August 2006
Accepted: 30 January 2007
Published: 30 January 2007
Fishes in the families Cichlidae and Labridae provide good probable examples of vertebrate adaptive radiations. Their spectacular trophic radiations have been widely assumed to be due to structural key innovation in pharyngeal jaw apparatus (PJA), but this idea has never been tested based on a reliable phylogeny. For the first step of evaluating the hypothesis, we investigated the phylogenetic positions of the components of the suborder Labroidei (including Pomacentridae and Embiotocidae in addition to Cichlidae and Labridae) within the Percomorpha, the most diversified (> 15,000 spp) crown clade of teleosts. We examined those based on 78 whole mitochondrial genome sequences (including 12 newly determined sequences) through partitioned Bayesian analyses with concatenated sequences (13,933 bp).
The resultant phylogenies indicated that the Labridae and the remaining three labroid families have diverged basally within the Percomorpha, and monophyly of the suborder was confidently rejected by statistical tests using Bayes factors.
The resultant phylogenies indicated that the specified PJA evolved independently at least twice, once in Labridae and once in the common ancestor of the remaining three labroid families (including the Cichlidae). Because the independent evolution of pharyngeal jaws appears to have been followed by trophic radiations, we consider that our result supports, from the aspect of historical repeatability, the idea that the evolution of the specialized PJA provided these lineages with the morphological potential for their spectacular trophic radiations. The present result will provide a new framework for the study of functional morphology and genetic basis of their PJA.
Although this scenario is one of the most well-known examples of an evolutionary key innovation (cited in Futuyma's  "Evolutionary Biology", the most widely-used textbook in the field), it has never been evaluated adequately. Unfortunately, the key innovation hypothesis on the specialized PJA has been proposed being coupled with a taxonomic hypothesis that has widely been accepted ever since: the family Labridae and Cichlidae have been classified into the suborder Labroidei, together with Pomacentridae and Embiotocidae, based on the sharing of the specialized PJA [10, 12]. According to Stiassny & Jensen , the specialized PJA has seven morphological features including the above three major ones. None of them is, however, unique to labroid fishes, although there is no other group known to have them all . Moreover, there is no independent, corroborative evidence for a monophyletic Labroidei outside of the specialized PJA [12, 13]. This situation has led some evolutionary biologists to doubt not only the taxonomic hypothesis but also the key innovation hypothesis (e.g. ). The doubt about the latter hypothesis is a little rash, however, because appropriateness of the evolutionary hypothesis essentially does not depend on that of the appended taxonomic hypothesis.
A discussion of the specialized PJA as a key innovation, however, hinges, first of all, upon a phylogenetic hypothesis independent of pharyngeal characters as below. The widely-accepted taxonomic hypothesis (monophyletic Labroidei) assumed that the evolution of the specialized PJA occurred only once during the percomorph evolution, a situation that has prevented the famous evolutionary example from being examined from the aspect of historical repeatability. The traditional taxonomic hypothesis based on the sharing of the specialized PJA, however, needs reconsideration. Phylogeny independent of the PJA may come to demonstrate multiple origins of the specialized PJA followed by trophic radiations, which support the key innovation hypothesis at least from the aspect of the historical repeatability. The hypothesis, however, must be tested statistically by comparing the diversity of the clades with the trait to that of the sister taxa without that trait [15, 16] when the true sister taxa are found in a future study.
For the first step of evaluating the key innovation hypothesis on the specialized PJA, we examined labroid monophyly based on 78 whole mitochondrial genome sequences from a variety of the Perciformes and related other ordinal taxa (collectively called "Percomorpha"). Usefulness of whole mitochondrial genome sequences (ca. 16,000 bp) for resolving higher-level phylogenetic relationships of teleost has been demonstrated in several studies [e.g. [17–19]]. Although non-monophyly of the labroid families have been suggested in previous studies [14, 20, 21], limited sequence data used in these studies (2,222 bp at longest) have precluded them from drawing definite conclusions, and no statistical test has been conducted for corroborating that hypothesis. Our analysis used the whole mitochondrial genome sequences (>13,000 bp per species), and confirmed the independent origins of the two diversified families by statistical tests using Bayes factors.
List of species used in this study, with DDBJ/EMBL/GenBank Accession Numbers. Classifications follow Nelson .
Percomorpha sensu Miya et al. 
Thunnus thynnus thynnus
In the resultant trees, every family of the suborder Labroidei (Cichlidae, Pomacentridae, Embiotocidae, and Labridae) was monophyletic with 100% PPs, and moreover, the family Labridae (clade A) was phylogenetically distant from the remaining three labroid families forming clade B. Monophyly of the Labroidei comprised of all the four families was confidently rejected by a statistical test using Bayes factor (2 ln = 540.28). The family Pseudochromidae (represented by Labracinus cyclophthalma, underlined in Fig. 3), which was sometimes thought to be a possible sister to the Labridae based on larval morphology , was placed not close to the Labridae (clade A), but more closely related to the clade B.
Overall relationships among the percomorph fishes examined here were quite similar to those obtained by unweighted and weighted maximum-parsimony analyses in Miya et al. , except for the newly added fishes, such as labroids. Although the phylogenies recovered many interesting inter-subordinal relationships that were not congruent with the traditional classification, sparse taxon sampling outside the Labroidei precludes further discussions.
The resultant trees and the statistical test with Bayes factor demonstrated that the Labridae and the remaining three labroid families including Cichlidae have independent origins within the Percomorpha. Based on molecular analyses, various relationships among the labroid families have been so far proposed [14, 20, 21]. Our analysis demonstrated entirely different relationships from the various ones with high statistical support, whereas none of the previous relationships was supported with high statistical value, probably owing to short sequences (2,222 bp at longest). The present analyses were conducted based on the long nucleotide sequences more than six times the lengths of those so far used (13,933 bp).
In spite of the high statistical supports, we must keep in mind that the sister relationships obtained here are highly tentative because of the sparse taxon sampling (76 species from 53 families) compared to the tremendous taxonomic diversity of the Percomorpha (ca. 14,000 species in 251 families; calculated from Nelson ). That is to say, although the present phylogeny demonstrated that the three labroid families, Cichlidae, Pomacentridae, and Embiotocidae, formed a well-supported monophyletic group (clade B) without non-labroid taxa, there remains a possibility that some non-labroid taxa would break into the group when more extensive taxon sampling is conducted. In the same way, there remains not a few possibilities that some yet-to-be-analyzed taxa (ca. 200 families left) would become sister group of the "labroid" taxa. To evaluate these possibilities, we need to conduct analyses with more extensive taxon sampling. In doing so, we would be able to evaluate "key innovation" hypothesis in a more rigorous manner using statistical test [15, 16].
Implications of the revised labroid phylogeny
The present result supports, from the aspect of historical repeatability, the famous (but never evaluated) evolutionary hypothesis that the dramatic labroid radiations in trophic ecology are due to evolution of specialized pharyngeal jaw apparatus (PJA). Based on the present phylogeny, a single origin of the specialized PJA in a common ancestor of the suborder Labroidei is unlikely, because such an evolutionary scenario requires subsequent losses of the complex structures in a wide variety of percomorphs. Rather, it is plausible that the specialized PJA has evolved twice independently in the Percomorpha history, in a common ancestor of the Labridae (clade A) and that of the Cichlidae, Pomacentridae, and Embiotocidae (clade B). Interestingly, evolutionary radiations in jaw forms and head shape (see Fig. 1) have occurred within both of the lineages with specialized PJA: one in the Labridae, and the other in the Cichlidae. The independent occurrence of such evolutionary succession (evolution of specialized PJA followed by trophic radiation) appears to support the idea that the evolution of the specialized PJA provided these lineages with the morphological "potential" for trophic radiations.
Despite their specialized PJA, the Pomacentridae and Embiotocidae seem not to have experienced clear trophic radiation. This may reflect the context-dependence of a key trait's effect on diversification . There are many potential factors that could influence how a key trait affects diversification. De Queiroz  categorized these factors roughly into three types: (1) other taxa, (2) other traits of the group itself, and (3) the physical environment. Evolutionary radiations in African Great Lake cichlids are supposedly facilitated by colonization of novel habitats (newly created lakes) without other competitive taxa . The Pomacentridae and Embiotocidae, both of which inhabit coastal marine waters, may not have had such opportunity to colonize novel habitats without competitors.
Evolution of the specialized "labroid" PJA
It is interesting fact that the complex "labroid" PJA evolved multiply in the Percomorpha history (other than the two lineages recognized here, two fossil labroid families have been recognized so far from the lower Middle Eocene [ca. 50 million years ago] of Monte Bolca, Italy, together with the oldest fossil records of the Labridae and Pomacentridae [26–29]). Is the acquisition of the "successful" PJA not so difficult for percomorph fishes, like parallel evolution of trophic morphologies in African Great Lake cichlids ? According to Stiassny & Jensen , typical labroid PJA has seven features (three of them shown in Fig. 2). Interestingly, none of them is unique to labroid fishes (; e.g. fused lower pharyngeal jaw is also shown in several families of Perciformes and members of Beloniformes), though there is no other group known to have them all . In addition to it, it seems that such features can be accomplished by a very simple change in ontogenetic mechanism (e.g. fusion of lower pharyngeal jaws into one unit) . The labroid PJA is therefore considered to have been achieved through the unique combination of several commonplace and slight morphological modifications. It will be worthwhile to study the molecular mechanisms underlying the morphological modifications, which would reveal molecular basis for the probable key innovations. According to Hulsey et al. , more than ten genes have been so far recognized which putatively influence teleost pharyngeal jaw.
The present result will provide a new view point for functional morphology of the unique organ (the "labroid" PJA). Although many works on the trophic apparatus have been so far published [e.g. [8–10, 12, 32–35]], comparison between labrid and other labroid fishes (including cichlids) as evolutionary independent lineages has never been conducted. Careful comparisons, particularly based on more extensive phylogenies including labroid and non-labroid fishes possessing "partly" specialized PJA (e.g. members of Beloniformes), will shed light on overlooked evolutionary implications.
The phylogenetic analyses of whole mitochondrial DNAs from various labroid and non-labroid "percomorph" fishes revealed that the Labridae and the remaining three labroid families have diverged basally within Percomorpha, indicating that the specified "labroid" PJA evolved independently at least twice, once in Labridae and once in the common ancestor of the remaining three labroid families (including the Cichlidae). Because both of the evolution appear to be followed by trophic radiations, we consider that our result supports the idea that the evolution of the specialized PJA provided these lineages with the morphological potential for their spectacular trophic radiations.
The present result will provide a new view point for functional morphology of the unique organ: although many works on the specialized PJA have been so far published, comparison between labrid and other labroid fishes (including cichlids) as evolutionary independent lineages has never been conducted. Careful comparisons, particularly based on more extensive phylogenies, will shed light on overlooked evolutionary implications.
To examine monophyly of suborder Labroidei, we used all four labroid families (Labridae, Cichlidae, Pomacentridae, and Embiotocidae) and a possible sister family for Labridae (Pseudochromidae) in our analysis. In addition, we used all available whole mitochondrial genome (mitogenome) sequences from the Percomorpha (sensu Miya et al. ) in the database (NCBI Organelle Genome Resources)  in the analyses. We chose two basal acanthomorph (but not percomorph) taxa, Polymixia japonica and Beryx splendens, as collective outgroups to root the trees. A list of taxa examined in this study is provided in Table 1 along with DDBJ/EMBL/GenBank accession.
DNA extraction, PCR, and sequencing
We amplified whole mitogenome sequences for the 11 labroids plus a single non-labroid species (Labracinus cyclophthalma: Pseudochromidae) using a long PCR technique . We used six fish-versatile long PCR primers in the following four combinations (S-LA-16S-L + S-LA-16S-H; L2508-16S + H1065-12S; L2508-16S + H12293-Leu; L12321-Leu + S-LA-16S-H; for locations and sequences of these primers, see Miya & Nishida [17, 38], Inoue et al. [39, 40], Ishiguro et al. , Kawaguchi et al. ) so as to amplify the entire mitogenome in two reactions. Long PCR reaction conditions followed Miya & Nishida . Dilution of the long PCR products was with sterile water (1:10–100) served for subsequent uses as short PCR templates.
We used a total of 162 fish-versatile PCR primers in various combinations to amplify contiguous, overlapping segments of the entire mitogenome for each of the 12 species (for locations and sequences of the primers, see Miya & Nishida [17, 38], Inoue et al. [39, 40, 43–45]; Ishiguro et al. , Kawaguchi et al. ). Short PCR reaction conditions followed Miya & Nishida .
Double-stranded short PCR products, purified using a Pre-Sequencing kit (USB), subsequently served for direct cycle sequencing with dye-labeled terminators (Applied Biosystems and Amersham Pharmacia) with the same primers for the short PCRs. All sequencing reactions were performed according to the manufacture's instructions. Labeled fragments were analyzed on model 373/377/3100 sequencers (Applied Biosystems).
Sequences editing and alignment
The sequence electropherograms were edited with the computer programs EDITVIEW (Version 1.01; Applied Biosystems). AUTOASSEMBLER (Version 2.1; Applied Biosystems) and DNASIS (Version 3.2; Hitachi Software Engineering) were used to concatenate the consensus mitogenomic sequences. Then, sequences were exported to phylogenetic software programs. To build our character matrix, we combined the 12 newly determined sequences with the 66 previously published sequences.
For each individual protein-coding gene, we manually aligned the sequences for the 78 species, with respect to the translated amino acid sequence, using MACCLADE . All stop codons and gaps were excluded from the subsequent phylogenetic analyses, as well as ambiguously aligned regions. The ND6 gene was also excluded because of its heterogeneous base composition and consistently poor phylogenetic performance . The 22 tRNA genes were also aligned manually, and ambiguous and gap alignment were excluded from the phylogenetic analyses. The 12S and 16S rRNA sequences were aligned using the software PROALIGN (Version 0.5) , and default setting parameters. Regions with posterior probabilities of ≤ 50% were excluded from the subsequent phylogenetic analyses.
When using a mitogenomic dataset for a particular taxonomic sampling, Simmons & Miya  empirically demonstrated that Bayesian analysis [49–52] is the most efficient character-based method for accurately reconstructing phylogeny. Following their recommendations, we used this method for constructing the labroid phylogeny.
Partitioned Bayesian phylogenetic analyses were performed with MRBAYES (Version 3.04 b)  for three different character matrices. The first matrix (dataset #1, 13,933 positions) includes concatenated nucleotide sequences from 12 protein-coding genes (10,758 positions), 22 transfer RNA genes (1,406 positions) and the two ribosomal RNA genes (1,769 positions). The second (dataset #2, 13,933 positions) includes the same set of characters with 3rd codon positions of the protein-coding genes converted into purine (R) and pyrimidine (Y), and the third matrix (dataset #3, 10,347 positions) with the 3rd codon positions excluded. Assuming that functional constraints on sequence evolution are more similar within codon positions (or types of molecules) across genes than across codon positions (or types of molecules) within genes, five partitions (1st, 2nd and 3rd codon positions of protein-coding genes, tRNA genes and rRNA genes) were set for the datasets #1 and #2, and four partitions (no third codon positions included) were done for dataset #3.
The general time reversible model with some sites assumed to be invariable and with variable sites assumed to follow a discrete gamma distribution [ref. ; GTR + I + Γ] was used for model of sequence evolution, as it was selected as best-fitting model with MRMODELTEST (Version 2)  for each partition except for positions with the RY-coding. We assumed that all of the model parameters were unlinked and rate multipliers were variable across partitions. For 3rd codon positions in the dataset #2, we used arbitrarily "A" and "C" instead of "R" and "Y" and set a single rate category (lset nst = 1) instead of six (lset nst = 6) to allow the program to estimate only transversional changes between purine (R) and pyrimidine (Y) nucleotides. MrBayes unnecessarily estimated transitional changes (A ↔ G or C ↔ T) when RY-coding was employed, which simply imposed fixed nucleotide compositions on R (A + G) and Y (C + T) during the calculation.
For each of the three matrices, two independent Bayesian analyses with Markov chain Monte Carlo (MCMC) process of 3,000,000 generations were performed. Parameter values and trees were sampled every 1,000 generations, and the samples before the convergence of the Markov chain were discarded for each run. The remaining samples from the two independent run were combined into a single file with a total of 5,200, 4,700 and 5,520 phylogenetic trees, respectively. The combined tree files were then imported into PAUP (Version 4.0b10)  to compute the 50% majority rule consensus trees. The percentages for the branches in the consensus trees represent the Bayesian posterior probabilities, which are the rough equivalent of a maximum likelihood search with bootstrapping .
Testing an alternative phylogenetic hypothesis
We tested an alternative hypothesis for labroid phylogeny using Bayes factors. Constrained Bayesian trees (imposing labroid monophyly on the trees) were estimated with MRBAYES (Version 3.1.2)  and harmonic means of log likelihood scores (calculated using "sump" command) were compared between constrained and unconstrained trees. Two independent partitioned Bayesian analyses with MCMC process of 3,000,000 generations were performed based on the dataset #2 (RY coding), which effectively removes likely "noise" from the dataset and avoids the apparent lack of signal by retaining all available positions in the dataset. Parameters and trees were sampled in the same way as described above, being combined into a single file with a total of 5,200 trees. The alternative trees were then compared using a Bayesian approach with Bayes factors , the traditional criterion of 2 ln Bayes factor of ≥ 10 being used as very strong evidence against the alternative hypothesis .
We wish to thank Australian Museum and Tokyo Sea Life Park for providing tissues. We also thank the following individuals and institutions for their help in collecting the specimens: J. Paxton, M. McGrouther, T. Mukai, and Iriomote Station, Tropical Biosphere Research Center, University of the Ryukyus. We are grateful to D. R. Bellwood for comments on an early version of the manuscript. This study was supported by Research Fellowships of Japan Society for the Promotion of Science for Young Scientists (No. 72503), and Grants-in-Aid from the Ministry of Education, Culture, Sports, Science, and Technology, Japan (12NP0201 and 15380131).
- Westneat MW, Alfaro ME: Phylogenetic relationships and evolutionary history of the reef fish family Labridae. Mol Phylogenet Evol. 2005, 36: 370-390.View ArticlePubMedGoogle Scholar
- Meyer A, Kocher TD, Basasibwaki P, Wilson AC: Monophyletic origin of Lake Victoria cichlid fishes suggested by mitochondrial DNA sequences. Nature. 1990, 347: 550-553.View ArticlePubMedGoogle Scholar
- Westneat MW, Alfaro ME, Wainwright PC, Bellwood DR, 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 CentralView ArticlePubMedGoogle Scholar
- FishBase. [http://www.fishbase.org]
- Wainwright PC, Bellwood DR, Westneat MW, Grubich JR, Hoey AS: A functional morphospace for the skull of labrid fishes: patterns of diversity in a complex biomechanical system. Biol J Linn Soc. 2004, 82: 1-25.View ArticleGoogle Scholar
- Westneat MW: A biomechanical model for analysis of muscle force, power output and lower jaw motion in fishes. J Theoret Biol. 2003, 223: 269-281.View ArticleGoogle Scholar
- Westneat MW: Evolution of levers and linkages in the feeding mechanisms of fishes. Integr Comp Biol. 2004, 44: 378-389.View ArticlePubMedGoogle Scholar
- Wainwright PC: Functional morphology of the pharyngeal jaw apparatus. Biomechanics of fishes. Edited by: Shadwick R, Lauder GV. 2005, Chicago: Elsevier, 77-101.View ArticleGoogle Scholar
- Liem KF: Evolutionary strategies and morphological innovations: cichlid pharyngeal jaws. Syst Zool. 1973, 22: 425-441.View ArticleGoogle Scholar
- Kaufman LS, Liem KF: Fishes of the suborder Labroidei (Pisces: Perciformes): phylogeny, ecology and evolutionary significance. Breviora. 1982, 472: 1-19.Google Scholar
- Futuyma DJ: Evolutionary biology. 1998, Sunderland: Sinauer, 3Google Scholar
- Stiassny MLJ, Jensen JS: Labroid interrelationships morphological complexity, key innovations, and the study of comparative diversity. Bull Mus Comp Zool. 1987, 151: 269-319.Google Scholar
- Johnson GD: Percomorph phylogeny: progress and problems. Bull Mar Sci. 1993, 52: 3-28.Google Scholar
- Streelman JT, Karl SA: Reconstructing labroid evolution with single-copy nuclear DNA. Proc R Soc B. 1997, 264: 1011-1020.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanderson MJ, Donoghue MJ: Reconstructing shifts in diversification rates on phylogenetic trees. Trends Ecol Evol. 1996, 11: 15-20.View ArticlePubMedGoogle Scholar
- Slowinski JB, Guyer C: Testing the stochasticity of patterns of organismal diversity: an improved null model. Am Nat. 1989, 134: 907-921.View ArticleGoogle Scholar
- Miya M, Nishida M: Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion. Mol Phylogenet Evol. 2000, 17: 437-455.View ArticlePubMedGoogle Scholar
- Miya M, Takeshima H, Endo H, Ishiguro NB, Inoue JG, Mukai T, Satoh TP, Yamaguchi M, Kawaguchi A, Mabuchi K, Shirai SM, Nishida M: Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences. Mol Phylogenet Evol. 2003, 26: 121-138.View ArticlePubMedGoogle Scholar
- Miya M, Satoh PS, Nishida M: The phylogenetic position of toadfishes (order Batrachoidiformes) in the higher ray-finned fish as inferred from partitioned Bayesian analysis of 102 whole mitochondrial genome sequences. Biol J Linn Soc. 2005, 85: 289-306.View ArticleGoogle Scholar
- Sparks JS, Smith WL: Phylogeny and biogeography of cichlid fishes (Teleostei: Perciformes: Cichlidae). Cladistics. 2004, 20: 501-517.View ArticleGoogle Scholar
- Dettai A, Lecointre G: Further support for the clades obtained by multiple molecular phylogenies in the acanthomorph bush. C R Biologies. 2005, 328: 674-689.View ArticlePubMedGoogle Scholar
- Richards WJ, Leis JM: Labroidei: development and relationships. Ontogeny and Systematics of fishes. Edited by: Moser HG. 1984, Gainesville: American Society of Ichthyologists and Herpetologists, 542-547.Google Scholar
- Nelson JS: Fishes of the world. 1994, New York: John Wiley & Sons, 3Google Scholar
- De Queiroz A: Contingent predictability in evolution: key traits and diversification. Syst Biol. 2002, 51: 917-929.View ArticlePubMedGoogle Scholar
- Schluter D: The ecology of adaptive radiation. 2000, Oxford: Oxford University PressGoogle Scholar
- Sorbini L, Boscaini E, Bannikov AF: On the morphology and systematics of the Eocene fish genus Tortonesia Sorbini from Bolca. Studi e Ricerche sui Giacimenti Terziari di Bolca. 1990, 4: 115-132.Google Scholar
- Bellwood DR: A new Eocene fossil fish, Sorbinia caudopunctata Gen. et Sp. nov., from Monte Bolca, Italy. Boll Mus Civ di St Nat Verona. 1995, 19: 461-474.Google Scholar
- Bellwood DR, Sorbini L: A review of the fossil record of the Pomacentridae (Teleostei: Labroidei) with a description of a new genus and species from the Eocene of Monte Bolca, Italy. Zool J Linn Soc. 1996, 117: 159-174.View ArticleGoogle Scholar
- Bellwood DR: Fossil pharyngognath fishes from Monte Bolca, Italy, with a description of a new pomacentrid genus and species. Studi e Ricerche sui Giacimenti Terziari di Bolca. 1999, 8: 207-217.Google Scholar
- Kocher TD, Conroy JA, McKaye KR, Stauffer JR: Similar morphologies of cichlid fish in Lakes Tanganyika and Malawi are due to convergence. Mol Phylogenet Evol. 1993, 2: 158-165.View ArticlePubMedGoogle Scholar
- Hulsey CD, Fraser GJ, Streelman JT: Evolution and development of complex biomechanical systems: 300 million years of fish jaws. Zebrafish. 2005, 2: 243-257.View ArticlePubMedGoogle Scholar
- Liem KF, Greenwood PH: A functional approach to the phylogeny of the pharyngognath teleosts. Amer Zool. 1981, 21: 83-101.View ArticleGoogle Scholar
- Galis F, Drucker EG: Pharyngeal biting mechanics in centrarchids and cichlids: insights into a key evolutionary innovation. J Evol Biol. 1996, 9: 641-670.View ArticleGoogle Scholar
- Wainwright PC, Bellwood DR: Ecomorphology of feeding in coral reef fishes. Coral reef fishes. Dynamics and diversity in a complex ecosystem. Edited by: Sale PF. 2002, Orlando: Academic Press, 33-55.View ArticleGoogle Scholar
- Hulsey CD: Function of a key morphological innovation: fusion of the cichlid pharyngeal jaw. Proc R Soc B. 2006, 273: 669-675.View ArticleGoogle Scholar
- NCBI Organelle Genome Resources. [http://www.ncbi.nlm.nih.gov/genomes/ORGANELLES/7898.html]
- Cheng S, Chang SY, Gravitt P, Respess R: Long PCR. Nature. 1994, 369: 684-685.View ArticlePubMedGoogle Scholar
- Miya M, Nishida M: Organization of the mitochondrial genome of a deep-sea fish Gonostoma gracile (Teleostei: Stomiiformes): first example of transfer RNA gene rearrangements in bony fishes. Mar Biotechnol. 1999, 1: 416-426.View ArticlePubMedGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of the Japanese sardine, Sardinops melanostictus. Fish Sci. 2000, 66: 924-932.View ArticleGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: A mitogenomic perspective on the basal teleostean phylogeny: resolving higher-level relationships with longer DNA sequences. Mol Phylogenet Evol. 2001, 20: 275-285.View ArticlePubMedGoogle Scholar
- Ishiguro N, Miya M, Nishida M: Complete mitochondrial DNA sequence of ayu, Plecoglossus altivelis. Fish Sci. 2001, 67: 474-481.View ArticleGoogle Scholar
- Kawaguchi A, Miya M, Nishida M: Complete mitochondrial DNA sequence of Aulopus japonicus (Teleostei: Aulopiformes), a basal Eurypterygii: longer DNA sequences and higher-level relationships. Ichthyol Res. 2001, 48: 213-223.View ArticleGoogle Scholar
- Inoue JG, Miya M, Aoyama J, Ishikawa S, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of the Japanese eel, Anguilla japonica. Fish Sci. 2001, 67: 118-125.View ArticleGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of Conger myriaster (Teleostei: Anguilliformes): novel gene order for vertebrate mitochondrial genomes and the phylogenetic implications for anguilliform families. J Mol Evol. 2001, 52: 311-320.PubMedGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of the Japanese anchovy, Engraulis japonicus. Fish Sci. 2001, 67: 828-835.View ArticleGoogle Scholar
- Maddison DR, Maddison WP: MACCLADE: Analysis of Phylogeny and Character Evolution, version 4.03. 2000, Sunderland: SinauerGoogle Scholar
- Löytynoja A, Milinkovitch MC: A hidden Markov model for progressive multiple alignment. Bioinformatics. 2003, 19: 1505-1513.View ArticlePubMedGoogle Scholar
- Simmons MP, Miya M: Efficiently resolving the basal clades of a phylogenetic tree using Bayesian and parsimony approach: a case study using mitogenomic data from 100 higher teleost fishes. Mol Phylogenet Evol. 2004, 31: 351-362.View ArticlePubMedGoogle Scholar
- Rannala B, Yang Z: Probability distribution of molecular evolutionary trees: a new method of phylogenetic inference. J Mol Evol. 1996, 43: 304-311.View ArticlePubMedGoogle Scholar
- Yang Z, Rannala B: Bayesian phylogenetic inference using DNA sequences: a Markov Chain Monte Carlo method. Mol Biol Evol. 1997, 14: 717-724.View ArticlePubMedGoogle Scholar
- Larget B, Simon DL: Markov Chain Monte Carlo algorithms for the Bayesian analysis of phylogenetic trees. Mol Biol Evol. 1999, 16: 750-759.View ArticleGoogle Scholar
- Mau B, Newton MA, Larget B: Bayesian phylogenetic inference via Markov Chain Monte Carlo methods. Biometrics. 1999, 55: 1-12.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian Inference of PhylogeneticTrees. Bioinformatics. 2001, 17: 754-755.View ArticlePubMedGoogle Scholar
- Yang Z: Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J Mol Evol. 1994, 39: 306-314.View ArticlePubMedGoogle Scholar
- Johan Nylander. [http://www.abc.se/~nylander/]
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods), version 4.0b10. 2002, Sunderland: Sinauer, [http://paup.csit.fsu.edu]Google Scholar
- Brandley MC, Schmitz A, Reeder TW: Partitioned Bayesian analyses, partition choice, and the phylogenetic relationships of scincid lizards. Syst Biol. 2005, 54: 373-390.View ArticlePubMedGoogle Scholar
- Kass RE, Raftery AE: Bayes factors. J Am Stat Assoc. 1995, 90: 773-795.View ArticleGoogle Scholar
- Fryer G, Iles TD: The cichlid fishes of the Great Lakes of Africa, their biology and evolution. 1972, London: Oliver & BoydGoogle Scholar
This article is published under license to 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.