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

In a previous study based on 100 whole mitochondrial genome (mitogenome) sequences, we sought to provide a new perspective on the ordinal relationships of higher ray-finned fish (Actinopterygii). The study left unexplored the phylogenetic position of toadfishes (order Batrachoidiformes), as data were unavailable owing to technical difficulties. In the present study, we successfully determined mitogenomic sequences for two toadfish species ( Batrachomoeus trispinosus and Porichthys myriaster ) and found that the difficulties resulted from unusual gene arrangements and associated repetitive non-coding sequences. Unambiguously aligned, concatenated mitogenomic sequences (13 461 bp) from 102 higher actinopterygians (excluding the ND6 gene and control region) were divided into five partitions (1st, 2nd and 3rd codon positions of the protein-coding genes, tRNA genes and rRNA genes) and partitioned Bayesian analyses were conducted. The resultant phylogenies strongly suggest that the toadfishes are not members of relatively primitive higher actinopterygians (Paracanthopterygii), but belong to a crown group of actinopterygians (Percomorpha), as was demonstrated for ophidiiform eels (Ophidiiformes) and anglerfishes (Lophiiformes) in the previous study. We propose revised limits of major unranked categories for higher actinopterygians and a new name (Berycomorpha) for a clade comprising two reciprocally paraphyletic orders (Beryciformes and Stephanoberyciformes) based on the present mitogenomic phylogenies.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 289–306.     

Mitogenomic exploration of higher teleostean phylogenies: a case study for moderate-scale evolutionary genomics with 38 newly determined complete mitochondrial DNA sequences

Although adequate resolution of higher-level relationships of organisms apparently requires longer DNA sequences than those currently being analyzed, limitations of time and resources present difficulties in obtaining such sequences from many taxa. For fishes, these difficulties have been overcome by the development of a PCR-based approach for sequencing the complete mitochondrial genome (mitogenome), which employs a long PCR technique and many fish-versatile PCR primers. In addition, recent studies have demonstrated that such mitogenomic data are useful and decisive in resolving persistent controversies over higher-level relationships of teleosts. As a first step toward resolution of higher teleostean relationships, which have been described as the "(unresolved) bush at the top of the tree," we investigated relationships using mitogenomic data from 48 purposefully chosen teleosts, of which those from 38 were newly determined during the present study (a total of 632,315 bp), using the above method. Maximum-parsimony and maximum-likelihood analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and third codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) from the 48 species. The resultant two trees from the two methods were well resolved and largely congruent, with many internal branches supported by high statistical values. The tree topologies themselves, however, exhibited considerable variation from the previous morphology-based cladistic hypotheses, with most of the latter being confidently rejected by the mitogenomic data. Such incongruence resulted largely from the phylogenetic positions or limits of long-standing problematic taxa, which were quite unexpected from previous morphological and molecular analyses. We concluded that the present study provided a basis of and guidelines for future investigations of teleostean evolutionary mitogenomics and that purposeful higher-density taxonomic sampling, subsequent sequencing efforts, and phylogenetic analyses of their mitogenomes may be decisive in resolving persistent controversies over higher-level relationships of teleosts, the most diversified group of all vertebrates, comprising over 23,500 extant species.     

Phylogenetic position of tetraodontiform fishes within the higher teleosts: Bayesian inferences based on 44 whole mitochondrial genome sequences

Tetraodontiformes includes approximately 350 species assigned to nine families, sharing several reduced morphological features of higher teleosts. The order has been accepted as a monophyletic group by many authors, although several alternative hypotheses exist regarding its phylogenetic position within the higher teleosts. To date, acanthuroids, zeiforms, and lophiiforms have been proposed as sister-groups of the tetraodontiforms. The monophyly and sister-group status was investigated using whole mitochondrial genome (mitogenome) sequences from 44 purposefully-chosen species (26 sequences newly-determined during the study) that fully represent the major tetraodontiform lineages plus all the groups that have been hypothesized as being close relatives. Partitioned Bayesian analyses were conducted with the three datasets that comprised concatenated nucleotide sequences from 13 protein-coding genes (with and without, or with RY-coding, 3rd codon positions), plus 22 transfer RNA and two ribosomal RNA genes. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high posterior probabilities. Mitogenomic data strongly supported the monophyly of tetraodontiform fishes, placing them as a sister-group of either Lophiiformes plus Caproidei or Caproidei only. The sister-group relationship between Acanthuroidei and Tetraodontiformes was statistically rejected using Bayes factors. These results were confirmed by a reanalysis of the previously published nuclear RAG1 gene sequences using the Bayesian method. Within the Tetraodontiformes, however, monophylies of the three superfamilies were not recovered and further taxonomic sampling and subsequent efforts should clarify these relationships.     

Mitogenomic circumscription of a novel percomorph fish clade mainly comprising “Syngnathoidei” (Teleostei)

Percomorpha, comprising about 60% of modern teleost fishes, has been described as the “(unresolved) bush at the top” of the tree, with its intrarelationships still being ambiguous owing to huge diversity (> 15,000 species). Recent molecular phylogenetic studies based on extensive taxon and character sampling, however, have revealed a number of unexpected clades of Percomorpha, and one of which is composed of Syngnathoidei (seahorses, pipefishes, and their relatives) plus several groups distributed across three different orders. To circumscribe the clade more definitely, we sampled several candidate taxa with reference to the previous studies and newly determined whole mitochondrial genome (mitogenome) sequences for 16 percomorph species across syngnathoids, dactylopterids, and their putatively closely-related fishes (Mullidae, Callionymoidei, Malacanthidae). Unambiguously aligned sequences (13,872 bp) from those 16 species plus 78 percomorphs and two outgroups (total 96 species) were subjected to partitioned Bayesian and maximum likelihood analyses. The resulting trees revealed a highly supported clade comprising seven families in Syngnathoidei (Gasterosteiformes), Dactylopteridae (Scorpaeniformes), Mullidae in Percoidei and two families in Callionymoidei (Perciformes). We herein proposed to call this clade “Syngnathiformes” following the latest nuclear DNA studies with some revisions on the included families.     

Mitogenomic evaluation of the unique facial nerve pattern as a phylogenetic marker within the percifom fishes (Teleostei: Percomorpha)

Percomorpha has been described as the “(unresolved) bush at the top” of the teleostean phylogenies and its intrarelationships are intrinsically difficult to solve because of its huge diversity (>15,000 spp.) and ill-defined higher taxa. Patterns of facial nerves, such as those of the ramus lateralis accessorius (RLA), have been considered as one of the candidate characters to delimit a monophyletic group within the percomorphs. Six families of the suborder Percoidei (Arripidae, Dichistiidae, Kyphosidae, Terapontidae, Kuhliidae, and Oplegnathidae) and suborder Stromateoidei (including six families) share the unique pattern 10 of RLA and it has been suggested that those fishes form a monophyletic group across the two perciform suborders. To evaluate the usefulness of the RLA pattern 10 as a phylogenetic marker within the percomorphs, we newly determined whole mitochondrial genome (mitogenome) sequences for the 13 species having RLA pattern 10 and their putatively, closely-related species (5 spp.). Unambiguously aligned sequences (14,263 bp) from those 18 species plus 50 percomrphs and two outgroups (total 70 species) were subjected to partitioned maximum likelihood and Bayesian analyses. The resulting trees clearly indicated that there were at least two independent origins of the unique facial nerve pattern: one in a common ancestor of Kyphosidae, Terapontidae, Kuhliidae, and Oplegnathidae and another one in that of the percoid Arripidae and Stromateoidei. Thus further detailed anatomical studies are needed to clarify the homology of this character between the two lineages. It should be noted that the latter two taxa (Arripidae and Stromateoidei) formed an unexpected, highly-supported monophyletic group together with Scombridae and possibly Chiasmodontidae and Bramidae, all lacking RLA pattern 10 (the former two are members of other perciform suborders Scombroidei and Trachinoidei, respectively). This novel, trans-subordinal clade has never been suggested by any morphological studies, although they share a common ecological characteristic, dwelling in the pelagic realm and often associated with long-distance migrations.     

A study on teleost phylogeny using specific antisera

The phylogenetic position of 26 teleost species has been studied by the use of a dot blot test based on eight rabbit antisera generated against serum components from representatives from various teleost orders. The isolated positions of Anguilliformes, Cypriniformes, Salmoniformes and Clupeiformes were confirmed. Lophiiformes showed some affinities for Pleuronectiformes and less for Gadiformes.     

A mitogenomic perspective on the basal teleostean phylogeny: resolving higher-level relationships with longer DNA sequences

A recent study demonstrated that mitochondrial genomic (mitogenomic) data comprising nucleotide sequences from the concatenated protein-coding (no 3rd codon positions) plus transfer RNA (stem regions only) genes reproduced the expected phylogeny of teleosts with high statistical support. We reexamined the interrelationships of the five major, basal teleostean lineages (Osteoglossomorpha, Elopomorpha, Clupeomorpha, Ostariophysi, and Protacanthopterygii; given various rankings) using mitogenomic data for which five alternative phylogenetic hypotheses have been previously proposed on the basis of both morphological and molecular analyses. In addition to previously determined complete mitochondrial DNA (mtDNA) sequences from eight basal teleosts and two outgroups, we determined the complete mtDNA sequences (excluding a portion of the control region) for two, purposefully chosen species of Osteoglossomorpha (Osteoglossum bicirrhosum and Pantodon buchholzi), and the data were subjected to maximumparsimony and maximum-likelihood analyses. The resultant tree topologies from the two methods were congruent, although they differed from any of the previously proposed hypotheses. Furthermore, the mitogenomic data confidently rejected all of these hypotheses with high statistical significance.     

Phylogenomic analysis of kinetoplastids supports that trypanosomatids arose from within bodonids

Kinetoplastids are a large group of free-living and parasitic eukaryotic flagellates, including the medically important trypanosomatids (e.g., Trypanosoma and Leishmania) and the widespread free-living and parasitic bodonids. Small subunit rRNA- and conserved protein-based phylogenies support the division of kinetoplastids into five orders (Prokinetoplastida, Neobodonida, Parabodonida, Eubodonida, and Trypanosomatida), but they produce incongruent results regarding their relative branching order, in particular for the position of the Trypanosomatida. In general, small subunit rRNA tends to support their early emergence, whereas protein phylogenies most often support a more recent origin from within bodonids. In order to resolve this question through a phylogenomic approach, we carried out massive parallel sequencing of cDNA from representatives of three bodonid orders (Bodo saltans -Eubodonida-, Procryptobia sorokini -Parabodonida-, and Rhynchomonas nasuta -Neobodonida-). We identified 64 well-conserved proteins shared by these species, four trypanosomatids, and two closely related outgroup species (Euglena gracilis and Diplonema papillatum). Phylogenetic analysis of a concatenated data set yielded a strongly supported tree showing the late emergence of trypanosomatids as a sister group of the Eubodonida. In addition, we identified homologues of proteins involved in trypanosomatid mitochondrial mRNA editing in the three bodonid species, suggesting that editing may be widespread in kinetoplastids. Comparison of expressed sequences from mitochondrial genes showed variability at U positions, in agreement with the existence of editing activity in the three bodonid orders most closely related to trypanosomatids (Neobodonida, Parabodonida, and Eubodonida). Mitochondrial mRNA editing appears to be an ancient phenomenon in kinetoplastids.     

Evolutionary history of anglerfishes (Teleostei: Lophiiformes): a mitogenomic perspective

Background  

The teleost order Lophiiformes, commonly known as the anglerfishes, contains a diverse array of marine fishes, ranging from benthic shallow-water dwellers to highly modified deep-sea midwater species. They comprise 321 living species placed in 68 genera, 18 families and 5 suborders, but approximately half of the species diversity is occupied by deep-sea ceratioids distributed among 11 families. The evolutionary origins of such remarkable habitat and species diversity, however, remain elusive because of the lack of fresh material for a majority of the deep-sea ceratioids and incompleteness of the fossil record across all of the Lophiiformes. To obtain a comprehensive picture of the phylogeny and evolutionary history of the anglerfishes, we assembled whole mitochondrial genome (mitogenome) sequences from 39 lophiiforms (33 newly determined during this study) representing all five suborders and 17 of the 18 families. Sequences of 77 higher teleosts including the 39 lophiiform sequences were unambiguously aligned and subjected to phylogenetic analysis and divergence time estimation.     

Species tree inference in a recent radiation of orioles (Genus Icterus): multiple markers and methods reveal cytonuclear discordance in the northern oriole group

Recent computational advances provide novel opportunities to infer species trees based on multiple independent loci. Thus, single gene trees no longer need suffice as proxies for species phylogenies. Several methods have been developed to deal with the challenges posed by incomplete and stochastic lineage sorting. In this study, we employed four Bayesian methods to infer the phylogeny of a clade of 11 recently diverged oriole species within the genus Icterus. We obtained well-resolved and mostly congruent phylogenies using a set of seven unlinked nuclear intron loci and sampling multiple individuals per species. Most notably, Bayesian concordance analysis generally agreed well with concatenation; the two methods agreed fully on eight of nine nodes. The coalescent-based method ^∗BEAST further supported six of these eight nodes. The fourth method used, BEST, failed to converge despite exhaustive efforts to optimize the tree search. Overall, the results obtained by new species tree methods and concatenation generally corroborate our findings from previous analyses and data sets. However, we found striking disagreement between mitochondrial and nuclear DNA involving relationships within the northern oriole group. Our results highlight the danger of reliance on mtDNA alone for phylogenetic inference. We demonstrate that in spite of low variability and incomplete lineage sorting, multiple nuclear loci can produce largely congruent phylogenies based on multiple species tree methods, even for very closely-related species.     

Basal euteleostean relationships: a mitogenomic perspective on the phylogenetic reality of the "Protacanthopterygii"

Higher-level relationships of the basal Euteleostei (=Protacanthopterygii) are so complex and controversial that at least nine different morphology-based phylogenetic hypotheses have been proposed during the last 30 years. Relationships of the Protacanthopterygii were investigated using mitochondrial genomic (mitogenomic) data from 34 purposefully chosen species (data for 12 species being newly determined during the study) that fully represented major basal euteleostean lineages and some basal teleosts plus neoteleosts as outgroups. Unweighted and weighted maximum parsimony (MP) and maximum likelihood (ML) analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and 3rd codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) from the 34 species. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high statistical values. Monophyly of the protacanthopterygians was confidently rejected by the mitogenomic data. Of the five major monophyletic groups that received high statistical support within the protacanthopterygians, a clade comprising members of the alepocephaloids was unexpectedly nested within the Otocephala, sister-group of the euteleosts. The remaining four major monophyletic groups, on the other hand, occupied phylogenetic positions intermediate between the otocephalans and neoteleosts, with a clade comprising esociforms + salmoniforms being more basal to the argentinoids and osmeroids. Although interrelationships of the latter two clades (argentinoids and osmeroids) with the neoteleosts remained ambiguous, the present results indicated explicitly that the protacanthopterygians as currently defined merely represent a collective, polyphyletic group of the basal euteleosts, located between the basal teleosts (elopomorphs and below) and neoteleosts (stomiiforms and above).     

Analysis of the Distribution of Bootstrap Tree Lengths Using the Maximum Parsimony Method

The bootstrap is an important tool for estimating the confidence interval of monophyletic groups within phylogenies. Although bootstrap analyses are used in most evolutionary studies, there is no clear consensus as how best to interpret bootstrap probability values. To study further the bootstrap method, nine small subunit ribosomal DNA (SSU rDNA) data sets were submitted to bootstrapped maximum parsimony (MP) analyses using unweighted and weighted sequence positions. Analyses of the lengths (i.e., parsimony steps) of the bootstrap trees show that the shape and mean of the bootstrap tree distribution may provide important insights into the evolutionary signal within the sequence data. With complex phylogenies containing nodes defined by short internal branches (multifurcations), the mean of the bootstrap tree distribution may differ by 2 standard deviations from the length of the best tree found from the original data set. Weighting sequence positions significantly increases the bootstrap values at internal nodes. There may, however, be strong bootstrap support for conflicting species groupings among different data sets. This phenomenon appears to result from a correlation between the topology of the tree used to create the weights and the topology of the bootstrap consensus tree inferred from the MP analysis of these weighted data. The analyses also show that characteristics of the bootstrap tree distribution (e.g., skewness) may be used to choose between alternative weighting schemes for phylogenetic analyses.     

Evolutionary relationships among heterokont algae (the autotrophic stramenopiles) based on combined analyses of small and large subunit ribosomal RNA

In order to study the phylogenetic relationships within the stramenopiles, and particularly among the heterokont algae, we have determined complete or nearly complete large-subunit ribosomal RNA sequences for different species of raphidophytes, phaeophytes, xanthophytes, chrysophytes, synurophytes and pinguiophytes. With the small- and large-subunit ribosomal RNA sequences of representatives for nearly all known groups of heterokont algae, phylogenetic trees were constructed from a concatenated alignment of both ribosomal RNAs, including more than 5,000 positions. By using different tree construction methods, inferred phylogenies showed phaeophytes and xanthophytes as sister taxa, as well as the pelagophytes and dictyochophytes, and the chrysophytes/synurophytes and eustigmatophytes. All these relationships are highly supported by bootstrap analysis. However, apart from these sister group relationships, very few other internodes are well resolved and most groups of heterokont algae seem to have diverged within a relatively short time frame.     

The mitochondrial genome of the pufferfish,Fugu rubripes,and ordinal teleostean relationships

The small nuclear genome of the pufferfish, Fugu rubripes (order Tetraodontiformes), makes this species highly interesting for genome research. In order to establish the phylogenetic position of the Tetraodontiformes relative to other teleostean orders that might also have a reduced nuclear genome size, we have sequenced the mitochondrial (mt) genome of the pufferfish. The gene order, nucleotide composition and evolutionary rate of the mt genome of the fugu correspond to those of other teleosts. This suggests that the evolution of this genome has not been affected by the processes that led to the dramatic reduction of the size of the nuclear genome of the fugu. The phylogenetic analyses, which were based on the concatenated amino acid sequences of twelve protein-coding mt genes, placed the fugu among the percomorphs. The affinities between the Tetraodontiformes and either the Perciformes or the Zeiformes were limited, however. The common notion of a separate euteleostean clade remained unsupported. The analyses did not support the traditional systematic understanding that the Clupeiformes constitute a basal teleostean lineage. In addition the findings strongly suggest that three teleostean orders, the Perciformes, Zeiformes and Scorpaeniformes, are paraphyletic.     

Gene trees, species trees, and morphology converge on a similar phylogeny of living gars (Actinopterygii: Holostei: Lepisosteidae), an ancient clade of ray-finned fishes

Extant gars represent the remaining members of a formerly diverse assemblage of ancient ray-finned fishes and have been the subject of multiple phylogenetic analyses using morphological data. Here, we present the first hypothesis of phylogenetic relationships among living gar species based on molecular data, through the examination of gene tree heterogeneity and coalescent species tree analyses of a portion of one mitochondrial (COI) and seven nuclear (ENC1, myh6, plagl2, S7 ribosomal protein intron 1, sreb2, tbr1, and zic1) genes. Individual gene trees displayed varying degrees of resolution with regards to species-level relationships, and the gene trees inferred from COI and the S7 intron were the only two that were completely resolved. Coalescent species tree analyses of nuclear genes resulted in a well-resolved and strongly supported phylogenetic tree of living gar species, for which Bayesian posterior node support was further improved by the inclusion of the mitochondrial gene. Species-level relationships among gars inferred from our molecular data set were highly congruent with previously published morphological phylogenies, with the exception of the placement of two species, Lepisosteus osseus and L. platostomus. Re-examination of the character coding used by previous authors provided partial resolution of this topological discordance, resulting in broad concordance in the phylogenies inferred from individual genes, the coalescent species tree analysis, and morphology. The completely resolved phylogeny inferred from the molecular data set with strong Bayesian posterior support at all nodes provided insights into the potential for introgressive hybridization and patterns of allopatric speciation in the evolutionary history of living gars, as well as a solid foundation for future examinations of functional diversification and evolutionary stasis in a "living fossil" lineage.     

Phylogenetic relationships among ascomycetes: evidence from an RNA polymerse II subunit

In an effort to establish a suitable alternative to the widely used 18S rRNA system for molecular systematics of fungi, we examined the nuclear gene RPB2, encoding the second largest subunit of RNA polymerase II. Because RPB2 is a single-copy gene of large size with a modest rate of evolutionary change, it provides good phylogenetic resolution of Ascomycota. While the RPB2 and 18S rDNA phylogenies were highly congruent, the RPB2 phylogeny did result in much higher bootstrap support for all the deeper branches within the orders and for several branches between orders of the Ascomycota. There are several strongly supported phylogenetic conclusions. The Ascomycota is composed of three major lineages: Archiascomycetes, Saccharomycetales, and Euascomycetes. Within the Euascomycetes, plectomycetes, and pyrenomycetes are monophyletic groups, and the Pleosporales and Dothideales are distinct sister groups within the Loculoascomycetes. We confirm the placement of Neolecta within the Archiascomycetes, suggesting that fruiting body formation and forcible discharge of ascospores were characters gained early in the evolution of the Ascomycota. These findings show that a slowly evolving protein-coding gene such as RPB2 is useful for diagnosing phylogenetic relationships among fungi.     

Central nervous system myelin of teleosts: comparative electrophoretic analysis of its proteins by staining and immunoblotting

CNS myelin was isolated from 24 teleostean fishes and the proteins were analyzed by staining and immunoblotting. All species showed a 36 K protein, two or more glycosylated hydrophobic intermediate protein (IP) components and several myelin basic protein bands (BP). The 36 K protein was specific for teleostean fishes. The IP and BP components displayed substantial variations in their proportions as well as in molecular sizes when comparing the different teleosts. This contrasts with CNS myelin proteins which appear more stable in terrestrial vertebrates.     

DNA barcoding Australia's fish species

Two hundred and seven species of fish, mostly Australian marine fish, were sequenced (barcoded) for a 655 bp region of the mitochondrial cytochrome oxidase subunit I gene (cox1). Most species were represented by multiple specimens, and 754 sequences were generated. The GC content of the 143 species of teleosts was higher than the 61 species of sharks and rays (47.1% versus 42.2%), largely due to a higher GC content of codon position 3 in the former (41.1% versus 29.9%). Rays had higher GC than sharks (44.7% versus 41.0%), again largely due to higher GC in the 3rd codon position in the former (36.3% versus 26.8%). Average within-species, genus, family, order and class Kimura two parameter (K2P) distances were 0.39%, 9.93%, 15.46%, 22.18% and 23.27%, respectively. All species could be differentiated by their cox1 sequence, although single individuals of each of two species had haplotypes characteristic of a congener. Although DNA barcoding aims to develop species identification systems, some phylogenetic signal was apparent in the data. In the neighbour-joining tree for all 754 sequences, four major clusters were apparent: chimaerids, rays, sharks and teleosts. Species within genera invariably clustered, and generally so did genera within families. Three taxonomic groups-dogfishes of the genus Squalus, flatheads of the family Platycephalidae, and tunas of the genus Thunnus-were examined more closely. The clades revealed after bootstrapping generally corresponded well with expectations. Individuals from operational taxonomic units designated as Squalus species B through F formed individual clades, supporting morphological evidence for each of these being separate species. We conclude that cox1 sequencing, or 'barcoding', can be used to identify fish species.     

The Infrabranchial Musculature and Its Bearing on the Phylogeny of Percomorph Fishes (Osteichthyes: Teleostei)

The muscles serving the ventral portion of the gill arches ( = infrabranchial musculature) are poorly known in bony fishes. A comparative analysis of the infrabranchial muscles in the major percomorph lineages reveals a large amount of phylogenetically-relevant information. Characters derived from this anatomical system are identified and discussed in light of current hypotheses of phylogenetic relationships among percomorphs. New evidence supports a sister-group relationship between the Batrachoidiformes and Lophiiformes and between the Callionymoidei and Gobiesocoidei. Investigated data also corroborate the existence of two monophyletic groups, one including the Pristolepididae, Badidae, and Nandidae, and a second clade consisting of all non-amarsipid stromateiforms. New synapomorphies are proposed for the Atherinomorphae, Blenniiformes, Lophiiformes, Scombroidei (including Sphyraenidae), and Gobiiformes. Within the latter order, the Rhyacichthyidae and Odontobutidae are supported as the successive sister families of all remaining gobiiforms. The present analysis further confirms the validity of infrabranchial musculature characters previously proposed to support the grouping of the Mugiliformes with the Atherinomorphae and the monophyly of the Labriformes with the possible inclusion of the Pholidichthyiformes. Interestingly, most hypotheses of relationships supported by the infrabranchial musculature have been advanced by preceding anatomists on the basis of distinct data sources, but were never recovered in recent molecular phylogenies. These conflicts clearly indicate the current unsatisfactory resolution of the higher-level phylogeny of percomorphs.