Evolutionary history of anglerfishes (Teleostei: Lophiiformes): a mitogenomic perspective
© Miya et al; licensee BioMed Central Ltd. 2010
Received: 30 August 2009
Accepted: 23 February 2010
Published: 23 February 2010
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.
Partitioned maximum likelihood analysis confidently recovered monophyly for all of the higher taxa (including the order itself) with the exception of the Thaumatichthyidae (Lasiognathus was deeply nested within the Oneirodidae). The mitogenomic trees strongly support the most basal and an apical position of the Lophioidei and a clade comprising Chaunacoidei + Ceratioidei, respectively, although alternative phylogenetic positions of the remaining two suborders (Antennarioidei and Ogcocephaloidei) with respect to the above two lineages are statistically indistinguishable. While morphology-based intra-subordinal relationships for relatively shallow, benthic dwellers (Lophioidei, Antennarioidei, Ogcocephaloidei, Chaunacoidei) are either congruent with or statistically indistinguishable from the present mitogenomic tree, those of the principally deep-sea midwater dwellers (Ceratioidei) cannot be reconciled with the molecular phylogeny. A relaxed molecular-clock Bayesian analysis of the divergence times suggests that all of the subordinal diversifications have occurred during a relatively short time period between 100 and 130 Myr ago (early to mid Cretaceous).
The mitogenomic analyses revealed previously unappreciated phylogenetic relationships among the lophiiform suborders and ceratioid familes. Although the latter relationships cannot be reconciled with the earlier hypotheses based on morphology, we found that simple exclusion of the reductive or simplified characters can alleviate some of the conflict. The acquisition of novel features, such as male dwarfism, bioluminescent lures, and unique reproductive modes allowed the deep-sea ceratioids to diversify rapidly in a largely unexploited, food-poor bathypelagic zone (200-2000 m depth) relative to the other lophiiforms occurring in shallow coastal areas.
Diversity of the Lophiiformes
Within the higher teleosts, the Lophiiformes has traditionally been allied with toadfishes of the order Batrachoidiformes, based primarily on osteological characters of the cranium [5–7]. Following the publication of the seminal work on higher-level relationships of teleosts by Greenwood et al.  and the advent of cladistic theory , both groups have been placed in the Paracanthopterygii, a presumed sister-group of the more derived Acanthopterygii . Other than the Lophiiformes and Batrachoidiformes, the original Paracanthopterygii  included those groups of fishes thought to be relatively primitive in the higher teleosts, such as Polymixiiformes, Percopsiformes, Ophidiiformes, Gadiformes, Zeioidei, Zoarcoidei and Gobiesocoidei. Subsequently, the taxonomic contents of the Paracanthopterygii have undergone significant changes, being finally reduced to five core orders (Percopsiformes, Ophidiiformes, Gadiformes, Batrachoidiformes, Lophiiformes) in an attempt to make the group monophyletic , and this taxonomic proposal has been followed in many reference books [11–14]. Thus the paracanthopterygian Lophiiformes (and its close association with the Batrachoidiformes) has been a prevailing view in the ichthyological community despite the lack of convincing evidence [1, 15, 16].
Recent molecular phylogenetic studies, however, have repeatedly cast doubt on such a paracanthoperygian position of the Lophiiformes within the higher teleosts [17–27]. These studies based on nucleotide sequences from both whole mitogenomes and various nuclear genes have strongly suggested that lophiiforms are highly derived teleosts, deeply nested in one of the larger percomorph clades, and that they are closely related to various percomorphs, such as the Tetraodontiformes, Caproidei, Acanthuroidei, Chaetodontidae, Pomacanthidae, Ephippidae and Drepanidae, all of them showing no indications of close affinity with the Lophiiformes before the advent of molecular phylogenetics. Significantly a mitochondrial phylogenomic study by Miya et al.  demonstrated that the Batrachoidiformes was deeply nested within a different percomorph clade consisting of the Synbranchiformes and Indostomiidae and a sister-group relationship between the Lophiiformes and Batrachoidiformes was confidently rejected by the Bayesian analyses. These novel relationships, however, have not been reflected in the most recently published classification of fishes .
In addition to the lack of available material of numerous rare taxa, the evolutionary history of the lophiiform fishes has remained elusive because of poor representation in the fossil record (but see [34–38]). Recent developments in the molecular estimation of divergence times, however, have provided promising tools to introduce time scales for the phylogenetic trees , thereby offering new insights into evolutionary history that cannot be inferred by the fossil data alone. Among the most significant advances common to these new methods is a departure from the molecular clock assumption and the use of time constraints at multiple nodes for rate calibration, usually based on fossil record. In higher teleosts, however, including lophiiforms, the fossil record is scarce and fragmentary, and alternative calibration points based on biogeographic events have proven useful for divergence time estimation. Azuma et al.  recently found that estimated divergence times of cichlid fishes showed excellent agreement with the history of Gondwanian fragmentation, arguing that such biogeographic events can be used as effective time constraints in dating teleostean divergences, which may be useful for dating lophiiform divergence times.
To address questions regarding the subordinal and familial relationships and evolutionary history of the Lophiiformes, we assembled the whole mitochondrial genome sequences from the 39 lophiiform species (33 sequences newly-determined during this study), representing all of the five suborders and 17 of the 18 families. Unambiguously aligned sequences (14,611 bp) from those 39 species plus 38 outgroups (total 77 species) were subjected to partitioned maximum likelihood (ML) analysis using RAxML . The resulting tree topology was then used to estimate the divergence time of the Lophiiformes using a Bayesian relaxed molecular-clock method to infer their evolutionary history, and patterns and rates of diversifications.
List of species used in this study
Outgroup (38 spp.)
Ingroup (39 spp.)
Lophiomus setigerus b
Caulophryne jordani c
DNA extraction, PCR, and Sequencing
We excised a small piece of epaxial musculature or fin-ray (ca. 0.25 g) from fresh or ethanol-fixed specimens of each species and preserved them in 99.5% ethanol. We extracted total genomic DNA from the tissue using QIAamp or DNeasy (Qiagen) following the manufacturer's protocol. We amplified the mitogenomes of the 33 lophiiform species in their entirety using a long PCR technique . We basically used seven fish-versatile PCR primers for the long PCR in the following four combinations (for locations of these primers, see [43–46]): L2508-16S (5'-CTC GGC AAA CAT AAG CCT CGC CTG TTT ACC AAA AAC-3') + H12293-Leu (5'-TTG CAC CAA GAG TTT TTG GTT CCT AAG ACC-3'); L2508-16S + H15149-CYB (5'-GGT GGC KCC TCA GAA GGA CAT TTG KCC TCA-3'); L8343-Lys (5'-AGC GTT GGC CTT TTA AGC TAA WGA TWG GTG-3') + H1065-12S (5'-GGC ATA GTG GGG TAT CTA ATC CCA GTT TGT-3'); and L12321-Leu (5'-GGT CTT AGG AAC CAA AAA CTC TTG GTG CAA-3') + S-LA-16S-H (5'-TGC ACC ATT RGG ATG TCC TGA TCC AAC ATC-3'). When we failed to cover the entire mitogenomes with these primer pairs, we used an additional five long PCR primers specifically designed to amplify the lophiiform mitogenomes: H8319-ANG-Lys (5'-GKA GKC ACC AKT TTT TAG MTT AAA AGG C-3'); L7567-ANG-Asp (5'-ACG CTG TTK TGT CAA GGC ARR AYT GTG GGT-3'); L10054-ANG-Gly (5'-CAC CWG GTC TTG GTT WAA MTC CMA GGA AAG-3'); H15149-ANG-CYB (5'-AGG TTK GTG ATG ACK GTK GCK CCT CA-3'); and L14850-ANG-CYB (5'-AAT ATC TCG GTK TGG TGG AAY TTT GGK TC-3'). Long PCR reaction conditions followed Miya and Nishida . Dilution of the long PCR products with TE buffer (1:10 to 100 depending on the concentration of the long PCR products) served as templates for subsequent short PCRs.
Standard set of 24 short PCR primer pairs for lophiiforms
Sequence (5' to 3')
We purified double-stranded short PCR products using a Pre-Sequencing kit (USB) for direct cycle sequencing with dye-labeled terminators (Applied Biosystems). We performed all sequencing reactions according to the manufacture's instructions with the same primers as those for the short PCRs. We analyzed labeled fragments on model 373/377/3100/3130xl sequencers (Applied Biosystems).
Sequence editing and alignment
We edited each sequence electropherogram with EditView (ver. 1.01; Applied Biosystems) and concatenated the multiple sequences using AutoAssembler (ver. 2.1; Applied Biosystems). We carefully checked the concatenated sequences using DNASIS (ver. 3.2; Hitachi Software Engineering) and created a sequence file for each gene. We compared the sequence files among closely related species to minimize sequence errors. Genes (or a portion of genes) that we were unable to sequence owing to technical difficulties were coded as missing.
To check sensitivity of additional taxon sampling of a number of the lophiiforms to the results reported in Yamanoue et al. , we used their pre-aligned sequences as a basis for further alignment with the newly determined sequences from 33 lophiiforms. Yamanoue et al.  aligned 13 protein-coding, two rRNA, and 22 tRNA genes using ProAlign ver. 0.5  and they used only those positions with posterior probabilities ≥70%. An exception to this was the alignment of tRNA genes, for which Yamanoue et al.  modified the alignment on the basis of the secondary structure, estimated with DNASIS. They used all the stem regions even if the aligned sequences were <70% posterior probabilities. Because the aligned sequences of Yamanoue et al.  included several overlapping positions between the two open reading frames (ATPase 8/6, ND4L/4, and ND5/6), we excluded those positions from the downstream genes (ATPase 6, ND4, and ND6).
To combine pre-aligned sequences from Yamanoue et al.  with the new sequences, we rearranged the dataset of Yamanoue et al.  into typical gene order of vertebrates (beginning from tRNA-Phe) and saved it in a FASTA format. The 33 newly determined sequences in the same format were concatenated to the rearranged, pre-aligned sequences and the dataset was subjected to the multiple alignment using MAFFT ver. 6 . We imported the aligned sequences into MacClade ver. 4.08  and removed the redundant regions appeared as gaps with slight modifications by eye to correctly reproduce the aligned sequences used in Yamanoue et al. . All the resulting gap positions from the alignment were coded as missing.
We divided unambiguously aligned sequences into five partitions (first, second, third codon positions, rRNA and tRNA genes) 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. We converted nucleotides at the third codon positions into purine (A/G) with "R" and pyrimidine (C/T) with "Y" (RY-coding; Phillips and Penny ) to take only transversions into account in the phylogenetic analysis following the recommendation of Saitoh et al. . This coding effectively removes likely "noise" from the dataset , and avoids the apparent lack of signal by retaining all available positions in the dataset. The "R" and "Y" were further recoded with "A" and "C," respectively, to avoid unnecessary estimation of transitional changes during the calculations using RAxML with the exception of an outgroup species (Polymixia japonicus) for running the program (RAxML does not accept all "A/C" partitions). We also constructed an additional two datasets that treat quickly saturated third codon positions differently (with or without third codon positions) to check sensitivity of the datasets to the phylogenetic analysis. The three datasets are designated as follows: 1) RY-coding (12n3rRTn); 2) all positions included (123nRTn); and 3) third codon positions excluded (12nRTn).
We subjected the above datasets to the partitioned maximum-likelihood (ML) analysis using RAxML ver. 7.0.4 . A general time reversible model with sites following a discrete gamma distribution (GTR + Γ; the model recommended by the author) was used and a rapid bootstrap (BS) analysis was conducted with 500 replications (-f a option). This option performs BS analysis using GTRCAT, which is GTR approximation with optimization of individual per-site substitution rates, and classification of those individual rates into certain number of rate categories. After implementing the BS analysis, the program uses every fifth BS tree as a starting point to search for the ML tree using GTR + Γ model of sequence evolution and saves the top 10 best-scoring ML trees (fast ML searches). Finally RAxML calculates more correct likelihood scores (slow ML searches) for those 10 trees and puts BS probabilities on the best-scoring ML tree.
Testing alternative phylogenetic hypotheses
We considered that the best-scoring ML tree resulting from 12n3rRTn (RY-coding) dataset as the best estimate of phylogeny because this coding effectively removes likely "noise" from the dataset and avoids the apparent lack of signal by retaining all available positions in the dataset (see discussions in Saitoh et al. ). Alternative tree topologies were thus individually compared to the resulting best-scoring ML tree derived from the 12n3rRTn dataset using the likelihood-based approximately unbiased (AU) test as implemented in CONSEL . The p-values from this test are calculated using the multi-scale bootstrap technique and are less biased than those of the conventional methods such as the bootstrap probability (BP), the Kishino-Hasegawa (KH) test and the Shimodaira-Hasegawa (SH) test .
We manually created the constrained tree topologies with reference to the alternative hypotheses using MacClade and then performed the RAxML analysis with each constraint using the -g option. We conducted the fast bootstrappings with 100 replicates as described above and the resulting best-scoring ML tree was considered as the constrained ML tree. The constrained and unconstrained ML trees (best-scoring ML tree without constraint) were used to compute the per-site log likelihood scores for each tree using the -f g option in RAxML and the output was subjected to CONSEL analysis to calculate statistical significance of the differences in the likelihood scores.
Tracing character evolution
Male sexual parasitism has been found among only the Ceratioidei . Its evolution was reconstructed on the best-scoring ML tree derived from 12n3rRTn dataset under an ML optimality criterion using Mesquite ver. 2.6 . The ML reconstruction methods find the ancestral states that maximize the probability the observed states would evolve under a stochastic model of evolution [57, 58]. The Mk1 model ("Markov k-state 1 parameter model"), a k-state generalization of the Jukes-Cantor model that corresponds to Lewis's Mk model , was used to trace the character evolution. Four character states were assigned to the male sexual parasitism based on extensive observations made by Pietsch  and Pietsch and Orr : males never attach to females (character state 0); males attach temporarily (state 1); males are facultative parasites (state 2); and males are obligate parasites (state 3).
Divergence time estimation
Because lophiiforms are rarely represented in the fossil record [34–37], the age of divergence of the lophiiform clades cannot be established precisely based on paleontological data alone. Thus external calibration points should be used at multiple nodes to estimate the divergence times of the Lophiiformes correctly. To that end, we used the mitogenomic dataset of Azuma et al.  who extensively sampled actinopterygians from the base to the top of the tree. Significantly the dataset of Azuma et al.  includes 1) all major lineages of the basal actinopterygians whose fossils and their relative phylogenetic positions are more reliable than those of the higher teleosts; and 2) all continental cichlids whose divergences show excellent agreement with the history of Gondwanian fragmentations.
Mitogenome sequences from the 39 lophiiforms were concatenated with the pre-aligned sequences used in Azuma et al.  in a FASTA format and the dataset was subjected to multiple alignment using MAFFT ver. 6  as described above. The dataset comprises 6966 positions from first and second codon positions of the 12 protein-coding genes, 1673 positions from the two rRNA genes and 1407 positions from the 22 tRNA genes (total 10,046 positions). The third codon positions of the protein-coding genes were entirely excluded because of their extremely accelerated rates of changes that may cause a high level of homoplasy at this taxonomic scale  and overestimation of divergence times .
Ideally all node ages for the 39 lophiiform species can be estimated in a single step; however, recent studies demonstrated that dense taxon sampling in a particular lineage (as has been done for the Lophiiformes in this study) tend to lead to overestimation of its age compared to the rest of the tree ("node-density effect" [61, 62]). To avoid such unnecessary overestimation, we retained a minimum number of taxa from each suborder in proportion to the logarithms of the species' diversity (Table 1). We selected the most distantly related species from each suborder to estimate crown node ages as correctly as possible. The nine selected species (three species from the most species-rich Ceratioidei and two from the rest of four suborders) are shown in Table 2 with asterisks. The resulting dataset contains 54 species used in Azuma et al  plus nine lophiiforms, with the total number of species being 63.
We used a relaxed molecular-clock method for dating analysis developed by Thorne and Kishino  to estimate divergence times. This method accommodates unlinked rate variation across different loci ("partitions" in this study), allows the use of time constraints on multiple divergences, and uses a Bayesian MCMC approach to approximate the posterior distribution of divergence times and rates based on a single tree topology estimated from the other method (ML tree in this study). A series of application in the software package multidistribute (v9/25/2003) were used for these analyses.
Baseml in PAML ver. 3.14 was used to estimate model parameters for each partition separately under the F84 + Γ model of sequence evolution (the most parameter-rich model implemented in multidistribute). Based on the outputs from baseml, branch lengths and the variance-covariance matrix were estimated using estbranches in multidistribute for each partition. Finally multidivtime in multidistribute was used to perform Bayesian MCMC analyses to approximate the posterior distribution of substitution rates, divergence times, and 95% credible intervals. In this step, multidivtime uses estimated branch lengths and the variance-covariance matrices from all partitions without information from the aligned sequences.
MCMC approximation with a burnin period of 100,000 cycles was obtained and every 100 cycles was taken to create a total of 10,000 samples. To diagnose possible failure of the Markov chains to converge to their stationary distribution, at least two replicate MCMC runs were performed with two different random seeds for each analysis.
Application of multidivtime requires values for the mean of the prior distribution for the time separating the ingroup root from the present (rttm) and its standard deviation (rttmsd), and we set conservative estimates of 4.2 (= 420 Myr ago [Ma]) and 4.2 SD, respectively. The tip-root branch lengths were calculated using TreeStat v. 1.1 http://tree.bio.ed.ac.uk/software/treestat/ for all terminals and their average was divided by rttm (4.2) to estimate rate of the root node (rtrate) and its standard deviation (rtratesd), which were set to 0.074 and 0.074, respectively. The priors for the mean of the Brownian motion constant, brownmean and brownsd, were both set to 0.5, specifying a relatively flexible prior.
List of time constraints used in divergence time estimation
The minimum age for the basal split of bony fish based on the earliest known acanthodian remains from Late Ordovician 
The minimum age for the Sarcopterygii/Actinopterygii split
The †Moythomasia fossil (actinopteran) from the Givetian/Eifelian boundary 
The minimum age for the Polypteriformes/Actinopteri split
The †Cosmoptychius fossil (neopterygian or actinopteran) from Tournasian 
The †Protopsephurus fossil (Polyodontidae) from Hauterivian (Cretaceous) 
The †Brachydegma fossil (stem amiids) from Artinskian (Permian) 
The †Yanbiania fossil (Hiodontidae) from the Lower Cretaceous 
The †Laeliichthys fossil (Osteoglossidae) from the Aptian (Cretaceous) 
The †Anaethalion, †Elopsomolos, and †Eoprotelops fossil (Elopomorpha) from Kimmeridgian (Jurassic) 
The †Lebonichthys (Albulidae) fossil from the Cenomanian (Cretaceous) 
The Conger (Congridae) and Anguilla (Anguillidae) fossils from the Ypresian (Tertiary) 
The †Tischlingerichthys fossil (Ostariophysi) from Tithonian (Jurassic) 
The †Knightia fossil (Clupeidae) from the Thanetian (Tertiary) 
The †Parabarbus fossil (Cyprinidae) from the Ypresian (Tertiary) 
The †Esteseox foxi fossil (Esociformes) from the Campanian (Cretaceous) 
The †Berycopsis fossil (Polymixiidae) from the Cenomanian (Cretaceous) 
The pleuronectiform fossil from the Ypresian (Tertiary) 
The tetraodontiform fossil from the Cenomanian 
The estimated divergence time between Takifugu and Tetraodon 
U 95 L 85
U 145 L 112
U 120 L 100
The lophiid fossil from Lutetian (Eocene) 
The ogcocephalid fossil from Lutetian (Eocene) 
The ceratioid fossil from upper Mohnian 
Net diversification rates
where n is the final number of lineages (present-day species diversity; Table 1) and t is the time interval considered (stem-group age).
Results and discussion
In the following sections, we describe and discuss the mitogenomic phylogenies and evolutionary history of the Lophiiformes. Whole mitogenomic phylogenetic analysis has been extremely useful in illuminating new ideas of interrelationships of fishes in particular, and renewed morphological analysis of these proposed relationships has often provided additional morphological support to challenge prevailing ideas of evolutionary relationships [27, 65]. We acknowledge, however, a phylogeny derived from the whole mitogenome only represents the mtDNA genealogy and may not necessarily match the evolution of the species under analysis. Because of the lack of recombination, the entire molecule of mtDNA has one molecular history that may be unusual because of various factors . Incongruence is a recurring problem at both higher and lower phylogenetic levels [67–70]. As noted by many authors, a broad approach to illuminating and reconciling this incongruence is to analyze other genetic evidence, such as that provided by nuclear DNA.
We newly determined the complete (or nearly complete) L-strand nucleotide sequences for mitogenomes of the 33 lophiiform species during this study and the sequences have been deposited in DDBJ/EMBL/GenBank under the accession numbers of AB282826-56 and AP005977-8 (Table 2). For Brachionichthys hirsutus (AB282832), however, we were unable to sequence a region spanning from tRNA-Gly to the control region (approximately 40% of the complete sequence) owing to technical difficulties and degradation of the tissues. The genome contents include two rRNA, 22 tRNA, and 13 protein-coding genes, plus the putative control region(s), as found in other vertebrates, and most of the genes are encoded on the H-strand, except for the ND6 and eight tRNA genes.
Patterns of intergenic non-coding sequences (≥ 100 bp) and their lengths (bp) in 23 oneiroid species
Met/N2 (111); A6/C3 (473); Glu/CB (257)
Ala/Asn (255); Tyr/C1 (314); N6/Glu (489)
A6/C3 (682); N6/Glu (186)
Such gene rearrangements and patterns of insertion sequences have been employed as useful phylogenetic markers in other fishes, as well as various metazoan animals, because they may represent uniquely derived characters shared by members of monophyletic groups (for reviews, see  but see also ). These genomic features also have been demonstrated to be useful in delimiting unexpected monophyletic groups in some teleosts, such as congroid eels  and macrouroid cods . However, the distributions of these unique genomic features across ceratioid families (not within-families; Table 5) are incongruent with the inferred inter-familial relationships derived from the nucleotide sequences (see below), suggesting either independent acquisitions or a single gain followed by independent losses of such unique features in a parsimony framework. Details of the gene rearrangements and patterns of insertion sequences in the Ceratioidei will be discussed elsewhere.
Monophyly and phylogenetic position of the Lophiiformes
Concerning the sister-group relationships of the Lophiiformes, no morphological study has provided a view that departs significantly from the previous paracanthopterygian notion advocated by Rosen and Patterson  and subsequently modified by Patterson and Rosen . Both mitogenomic  and nuclear gene  phylogenies, however, have convincingly demonstrated a percomorph relationship for the Lophiiformes and nullified the hypothesis of common ancestry with the Batrachoidiformes. In fact, use of the two whole mitogenome sequences from the Batrachoidiformes as only outgroups to root the lophiiform phylogenies disrupted the monophyletic Antennarioidei at the most basal position (as in Shedlock et al. ), followed by divergence of the Lophioidei, Ogcocephaloidei, and a clade comprising the Chaunacoidei and Ceratioidei at the top of the tree (results not shown). These subordinal relationships, particularly the non-monophyletic and most basal position of the Antennarioidei, are similar to those reported by Shedlock et al.  who used the batrachoidiform sequence as an only outgroup to root their tree.
We therefore excluded those two batrachoidiform sequences in the present study, thereby revealing a sister-group relationship either with the Caproidei alone (12n3rRTn and 123nRTn datasets) or with the Caproidei plus Tetraodontiformes (12nRTn dataset), as shown also by Yamanoue et al. . Nevertheless all nodal support values for these relationships were less than 50% bootstrap probabilities (BPs) and addition of unsampled members of the Percoidei (particularly putative members of Clade H in Kawahara et al. ; Yagishita et al. ) may eventually alter this picture of sister-group relationship of the Lophiiformes. Recently Li et al.  used 10 nuclear genes to analyze higher-level relationships of the actinopterygians and the only included lophiiform (a lophiid Lophius gastrophysus) was recovered as a sister species of two tetraodontiforms (Takifugu rubripes and Tetraodon nigroviridis). Although their dataset did not include a caproid sequence, it does appear from these and the other studies mentioned above that the tetraodontiforms are close relatives of the lophiiforms, within the Percomorpha.
Monophyly and interrelationships of the five suborders
Statistical comparisons among 15 alternative tree topologies of the four more derived suborders using AU test
Diff -ln L
We acknowledge, however, that no morphological data supports a sister-group relationship between the Chaunacoidei and Ceratioidei ( but see ). Instead, morphological data have indicated monophyly of the Ogcocephaloidei plus Ceratioidei with relatively strong statistical support (BS = 94%; Bremer index = 4; see ) with the following three unambiguous synapomorphies: 1) the first epibranchial is simple and without ligamentous connection to the second epibranchial (character 43); 2) the third cephalic dorsal-fin spine and pterygiophore are absent (character 60); and 3) the posttemporal is fused to the cranium (character 63). However, all of these characters appear in the Ogcocephaloidei and Ceratioidei to represent simplified or reductive trends, which are perhaps more likely to have occurred convergently, and the resulting homoplasy may undermine the robustness of the phylogenetic hypotheses based on morphology . Future evaluation of homology of these anatomical features, exploration of new morphological characters, and addition of molecular data from other genes may help resolve the conflict between these two different sources of phylogenetic information (for related discussion on the relationships within the Ceratioidei, see below).
The Lophioidei contains a single family, the Lophiidae, with 25 species distributed among four genera  (Table 1). Caruso  presented the first cladogram of lophiid genera based on 19 morphological characters (Figure 4C), of which 12 showed derived states shared by two or three genera. The reconstructed cladogram indicated the most basal position of Sladenia, followed by the divergence of Lophiodes and Lophiomus plus Lophius in sequential step-wise fashion, relationships that are fully congruent with the mitogenomic phylogenies, with all internal branches of the latter supported by 100% BPs (Figures 5, 6).
The Antennarioidei contains four families with 53 species distributed among 17 genera (Table 1). Pietsch and Grobecker  presented a cladogram of familial relationships of the suborder based on seven synapomorphies (Figure 4D), in which the Brachionichthyidae occupies the most basal position, followed by the divergence of Lophichthyidae, with Tetrabrachiidae and Antennariidae forming a sister-group at the top of the tree . Although we were unable to collect tissue samples from the only member of the Lophichthyidae (Lophichthys boschmai), the mitogenomic tree is completely congruent with the morphology-based phylogeny (Figures 5, 6).
Within the Antennariidae, Antennarius striatus is recovered as the sister of a terminal clade that includes Histrio histrio and A. coccineus, thus rendering Antennarius paraphyletic. The Antennariidae is by far the largest family of the suborder, including 45 species in 12 genera, of which only three species and two genera are included here. While our coverage of the Antennariidae is poor, an on-going molecular study by one of us (RJA), based on both mitochondrial and nuclear genes and considerably more taxa (25 species and 10 genera), also results in paraphyly for Antennarius. Thus, more extensive taxon sampling within Antennarius as well as within other antennariid genera is not likely to alter the topology shown here.
The Ogcocephaloidei contains a single family with 68 species distributed among 10 genera (Table 1). Endo and Shinohara , while describing a new species of the genus Coelophrys, cladistically analyzed nine morphological characters (all previously used in ) from nine of the 10 genera. As expected from such a small number of characters, resolution of the resulting cladogram was poor at the two most basal nodes (Figure 4E) and Coelophrys - an unusually globose genus among the typically dorsoventrally flattened ogcocephaloids - was placed at the top of the tree (Figure 4E). The placement of Coelophrys and the more basal Halieutaea in the cladogram (Figure 4E) agree with the mitogenomic phylogenies (Figures 5, 6), but the placement of Malthopsis and Zalieutes do not. A statistical test finds no significant difference between the morphological cladogram (Figure 4E) and the mitogenomic phylogeny (Figure 5) (AU test, P = 0.182), perhaps owing to the poor resolution of the morphological cladogram and low taxon sampling in the molecular phylogenies. Again more extensive taxon sampling will be required to obtain a better picture of their relationships.
The Chaunacoidei contains a single family with about 14 species divided between two genera  (Table 1). While we successfully obtained tissue samples from three species of the more common Chaunax, those from the rare genus Chaunacops were unavailable. Thus we are unable to evaluate monophyly for each of the two genera and to investigate their relationships. There is no phylogenetic hypothesis for chaunacoids at present.
The Ceratioidei contains 11 families with 160 species distributed among 35 genera  (Table 1). The first attempt to resolve relationships among ceratioid taxa after the advent of cladistic method  was that of Bertelsen . He admitted, however, that most of the derived osteological characters shared by two or more families are reductive states or loss of parts, and similarities among such characters may in many cases represent convergent development. Nevertheless, Bertelsen  presented a cladogram of the ceratioid taxa (Figure 4F), stating that the tree should be regarded only as "a very schematic compilation of expressed view." He concluded that future studies on additional characters and as yet unknown taxa may bring answers to at least some of the many questions about their phylogenetic relationships.
More recently, Pietsch and Orr , with the advantage of more than 20 years of additional accumulated data since Bertelsen's attempt , coupled with a re-examination of all previously identified characters and analyses of new characters, presented the first computer-assisted cladistic analysis of relationships of ceratioid families and genera (Figure 4G). In that study, Pietsch and Orr  showed two trees: one based on 71 morphological characters applicable to metamorphosed females (Figure 4G), and another one based on 17 morphological characters applicable to metamorphosed males and larvae, in addition to the 71 characters extracted from females, for a total of 88 characters. The latter tree was poorly resolved and Pietsch and Orr  thus considered the former as the best estimate of relationships.
Our dataset includes 23 species in 20 genera from all 11 ceratioid families. Our preferred dataset (12n3rRTn: RY-coding) reproduces the most basal Caulophrynidae, followed by divergence of the Ceratiidae, Gigantactinidae, Thaumatichthyidae plus Linophrynidae, Neoceratiidae plus Centrophrynidae, Oneirodidae (including Lasiognathus; see below), Himantolophidae, and Melanocetidae plus Diceratiidae at the top of the tree in sequential step-wise manner (Figure 5). More basal relationships among the seven families up to a clade comprising the Neoceratiidae plus Centrophrynidae are poorly resolved, with all internal branches supported by <60% BPs. Different treatments of the 3rd codon positions even yield different tree topologies, collapsing the most basal clade within the Ceratioidei in a strict consensus tree (Figure 6).
Interrelationships among the most apical four families (Oneirodidae, Himantolophidae, Melanocetidae, and Diceratiidae), on the other hand, are more robust with all internal branches supported by 99-100% BPs except for a sister-group relationship between the Melanocetidae and Diceratiidae (77-97% BPs; Figure 6). Shedlock et al.  also recovered an identical tree topology among the first three families (Oneirodidae, Himantolophidae, Melanocetidae) based on short sequences from the mitochondrial 16S rRNA gene, although their dataset lacked a member of the Diceratiidae. Lasiognathus (Figure 2H), long placed in the Thaumatichthyidae (Figure 2I), is here deeply nested within the Oneirodidae, and shown as the sister species of the derived oneirodid Puck pinnata at the top of the clade with 100% BPs (Figures 5, 6). The placement of Lasiognathus and Thaumatichthys in separate families was considered by Bertelsen and Struhsaker  who compared the osteology and pointed out that Lasiognathus appears more closely related to the Oneirodidae, but they chose to separate the two genera into different families. As expected from the most derived position of Lasiognathus within the oneirodids with the highest BPs (100%), monophyly of the traditional Thaumatichthyidae (including Lasiognathus) is confidently rejected by AU test (diff -ln L = 500.1; P < 0.0000).
Incongruence between morphology and molecule
As pointed out by Pietsch and Orr , who compared their tree with the unpublished molecular phylogeny provided by M.M. (referred to as Miya unpubl. data), the morphological hypothesis (Figure 4G) bears very little resemblance to the mitogenomic phylogenies (Figures 5, 6). Statistical differences between the constrained and unconstrained ML trees are so large (diff -ln L = 793.9, P < 0.0000 for ; diff -ln L = 1308.8, P < 0.0000 for ) that we are unable to reconcile these competing hypotheses. In fact, among the clades that differ between the two analyses, non-homoplastic morphological characters support only one clade (Himantolophidae, Diceratiidae, and Melanocetidae) and that with only a single character (8: the condition of the ventromedial extensions of the frontals).
Such remarkable incongruence between morphological and molecular hypotheses of ceratioid relationships requires an explanation. Although additional sequence data from other portions of the genome (e.g., nuclear genes) should be analyzed to confirm molecular conclusions [67, 69], Hedges and Sibley  argued that, in such cases of incongruence, morphological evidence should also be reevaluated. Following Hedges and Sibley's argument  and Bertelsen's empirical comments  that reductive states or loss of parts and similarities among such characters may in many cases represent convergent development, we have reviewed all of the 71 characters from the metamorphosed females and found the following 18 characters that are reductive, simplified, or absent for derived states (with the exception of those characters showing complete congruence with the molecular phylogenies; e.g., only autapomorphies for single families): vomerine teeth absent (character 3); parietal absent (9); pterosphenoids absent (10); endopterygoid absent (16); interopercle extremely reduced (23); rostral cartilage absent (26); maxillae considerably reduced (29); thick anterior maxillomandibular ligament very much reduced or absent (30); dentaries simple (31); first pharyngobranchial absent (39); first epibranchial absent (42); first epibranchial simple, not bearing a medial process (43); third hypobranchial absent (45); branchial teeth absent on the first three ceratobranchial (46); ninth or lower-most ray in caudal fin reduced (52); cephalic dorsa-fin spine absent (60); posttemporal is fused to the cranium (63); and pelvic bones reduced (66).
The revised cladogram (Figure 7) also recovers a monophyletic group comprising the Himantolophidae, Melanocetidae and Diceratiidae that is supported by 100% BPs (Figures 5, 6). Pietsch and Orr  observed that these three families uniquely share a single non-homoplastic morphological character (ventromedial extensions of the frontal that make no contact with the parasphenoid). In addition to these ceratioid relationships, monophyly of the Ogcocephaloidei + Ceratioidei is collapsed to form a trichotomy of these two suborders plus Chaunacoidei. Thus simple exclusion of reduced or simplified characters from the morphological dataset yields a tree that can be better reconciled with the molecular phylogenies (Figures 5, 6). However, simply deleting all reductive characters may also be misleading, by running the risk of rejecting informative characters. Homology of reductive morphological characters is commonly evaluated by ontogenetic analysis, but in the case of ceratioids, very little ontogenetic material is available for analysis [2, 33]. Considerably more work will be needed to further reconcile these competing phylogenetic hypotheses.
Evolution of male sexual parasitism
Pietsch and Orr  further argued that the disjunct pattern of sexual parasitism within ceratioids appears to be the result of independent acquisition among the various lineages rather than a repeated loss of this attribute within the suborder. To support this argument, Pietsch and Orr  listed many differences in the precise nature of male-female attachment among the various taxa , to the extent of the most extreme possibility being an independent acquisition of sexual parasitism within families, such as the Ceratiidae (Ceratias vs. Cryptopsaras) and Linophrynidae (Haplophryne vs. Linophryne). If so, evolution of sexual parasitism has independently occurred as many as seven times within the suborder (= the number of green or blue circles at terminal nodes in Figure 8). Similarly, although our simple character coding does not take into account such differences in male-female attachment, our ML reconstruction suggests that acquisition of this attribute has occurred at least five times during ceratioid evolution. Shedlock et al.  also found a paraphyletic pattern of sexual parasitism within the suborder in their much smaller dataset and suggested that the plasticity of this unique life history trait among vertebrates is likely shaped by a dynamic relationship between localized population densities and the feasibility of maintaining mate choice at low effective population size in the expanse of the deep ocean. Of course, it may be possible that availability of more specimens from these rare organisms will shed a new light for evolution of the male sexual parasitism.
Divergence time estimation
As Carnevale and Pietsch  stated, fishes of the order Lophiiformes are very rare in the fossil record and all of the recorded ages fall in the Cenozoic from 7.6 to 40 Myr ago (for details, see Table 4). Assuming a sister group relationship of the Lophiiformes and Tetraodontiformes, however, the origin of the modern Lophiiformes can be dated to the deep Mesozoic, because an articulated fossil that is convincingly assignable to the modern Tetraodontiformes was discovered from the mid Cretaceous (Cenomanian) 98 Myr ago . This fossil lineage would have appeared well after the divergence of the common ancestor of the Lophiiformes and Tetraodontiformes. Fossils are useful only for minimum time constraints to estimate divergence times of the Lophiiformes, as generally acknowledged [40, 84].
Comparisons of divergence time estimates between the present study and previous studies
This study (Figure 9)
Azuma et al. a
Setiamarga et al. 
Sarcopterygii vs. Actinopterygii
Teleostei vs. Neopterygii
Euteleostei vs. Otocephala
Cyprinus vs. Danio
Acanthopterygii vs. Paracanthopterygii
Percomorpha vs. Berycomorpha
Oryzias vs. Tetraodontiformes
Oryzias vs. Cichlidae
Gasterosteus vs. Tetraodontidae
Takifugu vs. Tetraodon
Unique among principally bathypelagic ceratioids are three species of the genus Thaumatichthys (Thaumatichthyidae; Figure 2I) that are abyssal-benthic, presumably staying in deep-sea bottom (> 1000 m) and luring prey items with esca "inside the mouth" . If this unusual life style had been attained concurrently in the origin of the common ancestor of Thaumatichthys, it took about 33 Myr after leaving the bottom of the sea around the continental shelves and subsequently returning to that unique benthic life style at greater depths.
Patterns and rates of diversification
Per-clade net diversification rates (events per lineage per Myr) for the five suborders of the Lophiiformes
Number of species
Divergence time (Myr ago)
The mitogenomic analyses demonstrated previously unappreciated phylogenetic relationships among the lophiiform suborders and deep-sea ceratioid familes. Although the latter relationships cannot be reconciled with the earlier hypotheses based on morphology, we found that simple exclusion of the reductive or simplified characters can alleviate some of the conflict. Reconstruction of the male reproductive modes of the ceratioids on the resultant phylogeny revealed complex evolutionary patterns of the sexual parasitism in males. A relaxed molecular-clock Bayesian analysis of the divergence times suggests that all of the subordinal diversifications have occurred during a relatively short time period between 100 and 130 Myr ago (early to mid Cretaceous). Comparisons of per-clade net diversification rates among the five lophiiform suborders suggest that the acquisition of novel features, such as male dwarfism, bioluminescent lures, and unique reproductive modes allowed the deep-sea ceratioids to diversify rapidly in a largely unexploited, food-poor bathypelagic zone (200-2000 m depth) relative to the other lophiiforms occurring in shallow coastal areas along continental shelves.
This study would not have been possible without donation of the study materials, for which we would like to thank A. Bentley, J.H. Caruso, H. Endo, A. Graham, K.E. Hartel, M. McGrouther, T.T. Sutton, E.O. Wiley, and M. Yamaguchi. We also thank J.G. Inoue for his kind advice in divergence time estimation, Y. Yamanoue for providing the pre-aligned sequences used in Yamanoue et al. , and C.P. Kenaley and D.E. Stevenson for helpful discussions. We thank the following for allowing us to reproduce their photographs: J.H. Caruso, P. David, D.B. Grobecker, C.P. Kenaley, R. Kuiter, D. Shale, and E. A. Widder. This study was supported in part by Grants-in-Aid from the Ministry of Education, Culture, Sports, Science and Technology, Japan (12NP0201, 15380131, 17207007, and 19207007); and by the U.S. National Science Foundation Grant DEB-0314637, T.W. Pietsch, principal investigator.
- Pietsch TW: Lophiiformes: development and relationships. Ontogeny and systematics of fishes. Edited by: Moser HG, Richards WJ, Cohen DE, Fahay MP, Kendall MPJ, Richardson SL. 1984, Lawrence, Kansas: American Society of Ichthyologists and Herpetologists, 320-325. vol. Special Publication 1Google Scholar
- Pietsch TW: Oceanic anglerfishes: extraordinary diversity in the deep sea. 2009, Berkeley, California: University of California PressGoogle Scholar
- Pietsch TW, Grobecker DB: Frogfishes of the world: systematics, zoogeography, and behavioral ecology. 1987, Stanford, California: Stanford University PressGoogle Scholar
- Pietsch TW: Dimorphism, parasitism, and sex revisited: modes of reproduction among deep-sea ceratioid anglerfishes (Teleostei: Lophiiformes). Ichthyological Research. 2005, 52: 207-236. 10.1007/s10228-005-0286-2.View ArticleGoogle Scholar
- Gregory WK: Fish skulls: A study of the evolution of natural mechanisms. Transactions of American Philosophical Society. 1933, 23: 75-481.View ArticleGoogle Scholar
- Regan CT: The classification of the teleostean fishes of the order Pediculati. Annals and Magazine of Natural History, Series. 1912, 8: 277-289.Google Scholar
- Rosen DE, Patterson C: The structure and relationships of the paracanthopterygian fishes. Bulletin of the American Museum of Natural History. 1969, 141: 359-474.Google Scholar
- Greenwood PH, Rosen DE, Weitzman SH, Myers GS: Phyletic studies of teleostean fishes, with a provisional classification of living forms. Bulletin of the American Museum of Natural History. 1966, 131: 339-456.Google Scholar
- Hennig W: Phylogenetic systematics. 1966, Urbana, IL: University of Illinois PressGoogle Scholar
- Patterson C, Rosen DE: The Paracanthopterygii revisited: order and disorder. Papers on the systematics of gadiform fishes. Edited by: Cohen DM. 1989, Los Angeles, California: Natural History Museum of Los Angeles County, 32: 5-36.Google Scholar
- Helfman GS, Collette BB, Facey DE: The diversity of fishes. 1997, Massachusetts: Blackwell PublishingGoogle Scholar
- Helfman GS, Collette BB, Facey DE, Bowen BW: The diversity of fishes: biology, evolution and ecology. 2009, West Sussex, UK: John Wiley & Sons, 2Google Scholar
- Nakabo T: Fishes of Japan with pictorial keys to the species. 2000, Tokyo: Tokai University Press, 2Google Scholar
- Nelson JS: Fishes of the world. 2006, Hoboken, New Jersey: John Wiley & Sons, 4Google Scholar
- Nelson JS: Fishes of the World. 1994, New York: John Wiley & Sons, 3Google Scholar
- Rosen DE: An essay on euteleostean classification. American Museum Novitates. 1985, 2827: 1-57.Google Scholar
- Chen WJ, Bonillo C, Lecointre G: Repeatability of clades as a criterion of reliability: a case study for molecular phylogeny of Acanthomorpha (Teleostei) with larger number of taxa. Molecular Phylogenetics and Evolution. 2003, 26: 262-288. 10.1016/S1055-7903(02)00371-8.View ArticlePubMedGoogle Scholar
- Dettai A, Bailly N, Vignes-Lebbe R, Lecointre G: Metacanthomorpha: essay on a phylogeny-oriented database for morphology--the acanthomorph (Teleostei) example. Systematic Biology. 2004, 53: 822-834. 10.1080/10635150490522313.View ArticlePubMedGoogle Scholar
- Dettai A, Lecointre G: Further support for the clades obtained by multiple molecular phylogenies in the acanthomorph bush. Comptes Rendus Biologies. 2005, 328: 674-689. 10.1016/j.crvi.2005.04.002.View ArticlePubMedGoogle Scholar
- Holcroft NI: A molecular test of alternative hypotheses of tetraodontiform (Acanthomorpha: Tetraodontiformes) sister group relationships using data from the RAG1 gene. Molecular Phylogenetics and Evolution. 2004, 32: 749-760. 10.1016/j.ympev.2004.04.002.View ArticlePubMedGoogle Scholar
- Holcroft NI: A molecular analysis of the interrelationships of tetraodontiform fishes (Acanthomorpha: Tetraodontiformes). Molecular Phylogenetics and Evolution. 2005, 34: 525-544. 10.1016/j.ympev.2004.11.003.View ArticlePubMedGoogle Scholar
- Li B, Dettaï A, Cruaud C, Couloux A, Desoutter-Meniger M, Lecointre G: RNF213, a new nuclear marker for acanthomorph phylogeny. Molecular Phylogenetics and Evolution. 2009, 50: 345-363. 10.1016/j.ympev.2008.11.013.View ArticlePubMedGoogle Scholar
- Li C, Lu G, Orti G: Optimal data partitioning and a test case for ray-finned fishes (Actinopterygii) based on ten nuclear loci. Systematic Biology. 2008, 57: 519-539. 10.1080/10635150802206883.View ArticlePubMedGoogle Scholar
- Miya M, Holcroft NI, Satoh TP, Yamaguchi M, Nishida M, Wiley EO: Mitochondrial genome and a nuclear gene indicate a novel phylogenetic position of deep-sea tube-eye fish (Stylephoridae). Ichthyological Research. 2007, 54: 323-332. 10.1007/s10228-007-0408-0.View ArticleGoogle Scholar
- Miya M, Satoh TP, 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. Biological Journal of the Linnean Society. 2005, 85: 289-306. 10.1111/j.1095-8312.2005.00483.x.View ArticleGoogle 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. Molecular Phylogenetics and Evolution. 2003, 26: 121-138. 10.1016/S1055-7903(02)00332-9.View ArticlePubMedGoogle Scholar
- Yamanoue Y, Miya M, Matsuura K, Yagishita N, Mabuchi K, Sakai H, Katoh M, Nishida M: Phylogenetic position of tetraodontiform fishes within the higher teleosts: Bayesian inferences based on 44 whole mitochondrial genome sequences. Molecular Phylogenetics and Evolution. 2007, 45: 89-101. 10.1016/j.ympev.2007.03.008.View ArticlePubMedGoogle Scholar
- Pietsch TW: The osteology and relationships of the anglerfish genus Tetrabrachium, with comments on lophiiform classification. Fishery Bulletin. 1981, 79: 387-419.Google Scholar
- Shedlock AM, Pietsch TW, Haygood MG, Bentzen P, Hasegawa M: Molecular systematics and life history evolution of anglerfishes (Teleostei: Lophiiformes): Evidence from mitochondrial DNA. Steenstrupia. 2004, 28: 129-144.Google Scholar
- Caruso JH: The systematics and distribution of the lophiid anglerfishes: III. Intergeneric relationships. Copeia. 1985, 1985: 870-875. 10.2307/1445235.View ArticleGoogle Scholar
- Endo H, Shinohara G: A new batfish, Coelophrys bradburyae (Lophiiformes: Ogcocephalidae) from Japan, with comments on the evolutionary relationships of the genus. Ichthyological Research. 1999, 46: 359-365. 10.1007/BF02673978.View ArticleGoogle Scholar
- Bertelsen E: Ceratioidei: Development and relationships. Ontogeny and systematics of fishes. Edited by: Moser HG, Richards WJ, Cohen DE, Fahay MP, Kendall MPJ, Richardson SL. 1984, Lawrence, Kansas: American Society of Ichthyologists and Herpetologists, 325-334. vol. Special Publication 1Google Scholar
- Pietsch TW, Orr JW: Phylogenetic relationships of deep-sea anglerfishes of the suborder Ceratioidei (Teleostei: Lophiiformes) based on morphology. Copeia. 2007, 2007: 1-34. 10.1643/0045-8511(2007)7[1:PRODAO]2.0.CO;2.View ArticleGoogle Scholar
- Carnevale G, Pietsch TW: Filling the gap: a fossil frogfish, genus Antennarius (Teleostei, Lophiiformes, Antennariidae), from the Miocene of Algeria. Journal of Zoology. 2006, 270: 448-457. 10.1111/j.1469-7998.2006.00163.x.View ArticleGoogle Scholar
- Carnevale G, Pietsch TW: The deep-sea anglerfish genus Acentrophryne (Teleostei, Ceratioidei, Linophrynidae) in the Miocene of California. Journal of Vertebrate Paleontology. 2009, 29: 372-378. 10.1671/039.029.0232.View ArticleGoogle Scholar
- Carnevale G, Pietsch TW: An Eocene frogfish from Monte Bolca, Italy: The earliest known skeletal record for the family. Palaeontology. 2009, 52: 745-752. 10.1111/j.1475-4983.2009.00874.x.View ArticleGoogle Scholar
- Carnevale G, Pietsch TW: Eocene handfishes from Monte Bolca, with description of a new genus and species, and a phylogeny of the family Brachionichthyidae (Teleostei: Lophiiformes). Zoological Journal of the Linnean Society. 2009.Google Scholar
- Carnevale G, Pietsch TW, Takeuchi GT, Huddleston RW: Fossil ceratioid anglerfishes (Teleostei: Lophiiformes) from the Miocene of the Los Angeles Basin, California. Journal of Paleontology. 2008, 82: 996-1008. 10.1666/07-113.1.View ArticleGoogle Scholar
- Yang Z: Computational molecular evolution. 2006, New York: Oxford University Press, USAView ArticleGoogle Scholar
- Azuma Y, Kumazawa Y, Miya M, Mabuchi K, Nishida M: Mitogenomic evaluation of the historical biogeography of cichlids toward reliable dating of teleostean divergences. BMC Evolutionary Biology. 2008, 8: 215-10.1186/1471-2148-8-215.PubMed CentralView ArticlePubMedGoogle Scholar
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.View ArticlePubMedGoogle Scholar
- Cheng S, Higuchi R, Stoneking M: Complete mitochondrial genome amplification. Nature Genetics. 1994, 7: 350-351. 10.1038/ng0794-350.View ArticlePubMedGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of the Japanese sardine Sardinops melanostictus. Fisheries Science. 2000, 66: 924-932. 10.1046/j.1444-2906.2000.00148.x.View ArticleGoogle Scholar
- Inoue JG, Miya M, Tsukamoto K, Nishida M: Complete mitochondrial DNA sequence of the Japanese anchovy Engraulis japonicus. Fisheries Science. 2001, 67: 828-835. 10.1046/j.1444-2906.2001.00329.x.View ArticleGoogle Scholar
- Ishiguro N, Miya M, Nishida M: Complete mitochondrial DNA sequence of ayu Plecoglossus altivelis. Fisheries Science. 2001, 67: 474-481. 10.1046/j.1444-2906.2001.00283.x.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. Ichthyological Research. 2001, 48: 213-223. 10.1007/s10228-001-8139-0.View ArticleGoogle 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. Marine Biotechnology. 1999, 1: 416-426. 10.1007/PL00011798.View ArticlePubMedGoogle Scholar
- Löytynoja A, Milinkovitch MC: A hidden Markov model for progressive multiple alignment. Bioinformatics. 2003, 19: 1505-1513. 10.1093/bioinformatics/btg193.View ArticlePubMedGoogle Scholar
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Briefings in Bioinformatics. 2008, 9: 286-298. 10.1093/bib/bbn013.View ArticlePubMedGoogle Scholar
- Maddison DR, Maddison WP: MacClade Version 4. 2003, Sunderland, Massachusetts: Sinauer AssociatesGoogle Scholar
- Phillips MJ, Penny D: The root of the mammalian tree inferred from whole mitochondrial genomes. Molecular Phylogenetics and Evolution. 2003, 28: 171-185. 10.1016/S1055-7903(03)00057-5.View ArticlePubMedGoogle Scholar
- Saitoh K, Sado T, Mayden RL, Hanzawa N, Nakamura K, Nishida M, Miya M: Mitogenomic evolution and interrelationships of the Cypriniformes (Actinopterygii: Ostariophysi): the first evidence toward resolution of higher-level relationships of the world's largest freshwater fish clade based on 59 whole mitogenome sequences. Journal of Molecular Evolution. 2006, 63: 826-841. 10.1007/s00239-005-0293-y.View ArticlePubMedGoogle Scholar
- Miya M, Nishida M: Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion. Molecular Phylogenetics and Evolution. 2000, 17: 437-455. 10.1006/mpev.2000.0839.View ArticlePubMedGoogle Scholar
- Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Systematic Biology. 2002, 51: 492-508. 10.1080/10635150290069913.View ArticlePubMedGoogle Scholar
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17: 1246-1247. 10.1093/bioinformatics/17.12.1246.View ArticlePubMedGoogle Scholar
- Maddison WP, Maddison DR: Mesquite: a modular system for evolutionary analysis. Version 2.71. 2009, [http://mesquiteproject.org]Google Scholar
- Pagel M: The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies. Systematic Biology. 1999, 48: 612-622. 10.1080/106351599260184.View ArticleGoogle Scholar
- Schluter D, Price T, Mooers AO, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51: 1699-1711. 10.2307/2410994.View ArticleGoogle Scholar
- Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data. Systematic Biology. 2001, 50: 913-925. 10.1080/106351501753462876.View ArticlePubMedGoogle Scholar
- Benton MJ, Ayala FJ: Dating the tree of life. Science. 2003, 300: 1698-1700. 10.1126/science.1077795.View ArticlePubMedGoogle Scholar
- Linder HP, Hardy CR, Rutschmann F: Taxon sampling effects in molecular clock dating: an example from the African Restionaceae. Molecular Phylogenetics and Evolution. 2005, 35: 569-582. 10.1016/j.ympev.2004.12.006.View ArticlePubMedGoogle Scholar
- Yoder AD, Yang Z: Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context. Molecular Ecology. 2004, 13: 757-773. 10.1046/j.1365-294X.2004.02106.x.View ArticlePubMedGoogle Scholar
- Thorne JL, Kishino H: Divergence time and evolutionary rate estimation with multilocus data. Systematic Biology. 2002, 51: 689-702. 10.1080/10635150290102456.View ArticlePubMedGoogle Scholar
- Magallónn S, Sanderson MJ: Absolute diversification rates in angiosperm clades. Evolution. 2001, 55: 1762-1780.View ArticleGoogle Scholar
- Holcroft NI, Wiley EO: Acanthuroid relationships revisited: a new nuclear gene-based analysis that incorporates tetraodontiform representatives. Ichthyological Research. 2008, 55: 274-283. 10.1007/s10228-007-0026-x.View ArticleGoogle Scholar
- Ballard JW, Whitlock MC: The incomplete natural history of mitochondria. Molecular Ecology. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.View ArticlePubMedGoogle Scholar
- Ballard JWO, Chernoff B, James AC: Divergence of mitochondrial DNA is not corroborated by nuclear DNA, morphology, or behavior in Drosophila simulans. Evolution. 2002, 56: 527-545.View ArticlePubMedGoogle Scholar
- Cameron SL, Miller KB, D'Haese CA, Whiting MF, Barker SC: Mitochondrial genome data alone are not enough to unambiguously resolve the relationships of Entognatha, Insecta and Crustacea sensu lato (Arthropoda). Cladistics. 2004, 20: 534-557. 10.1111/j.1096-0031.2004.00040.x.View ArticleGoogle Scholar
- Leache AD, McGuire JA: Phylogenetic relationships of horned lizards (Phrynosoma) based on nuclear and mitochondrial data: Evidence for a misleading mitochondrial gene tree. Molecular Phylogenetics and Evolution. 2006, 39: 628-644. 10.1016/j.ympev.2005.12.016.View ArticlePubMedGoogle Scholar
- Springer MS, DeBry RW, Douady C, Amrine HM, Madsen O, de Jong WW, Stanhope MJ: Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction. Molecular Biology and Evolution. 2001, 18: 132-143.View ArticlePubMedGoogle Scholar
- Boore JL: Animal mitochondrial genomes. Nucleic Acids Research. 1999, 27: 1767-1780. 10.1093/nar/27.8.1767.PubMed CentralView ArticlePubMedGoogle Scholar
- Mindell DP, Sorenson MD, Dimcheff DE: Multiple independent origins of mitochondrial gene order in birds. Proceedings of the National Academy of Sciences. 1998, 95: 10693-10697. 10.1073/pnas.95.18.10693.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. Journal of Molecular Evolution. 2001, 52: 311-320.PubMedGoogle Scholar
- Satoh TP, Miya M, Endo H, Nishida M: Round and pointed-head grenadier fishes (Actinopterygii: Gadiformes) represent a single sister group: evidence from the complete mitochondrial genome sequences. Molecular Phylogenetics and Evolution. 2006, 40: 129-138. 10.1016/j.ympev.2006.02.014.View ArticlePubMedGoogle Scholar
- Kawahara R, Miya M, Mabuchi K, Lavoue S, Inoue JG, Satoh TP, Kawaguchi A, Nishida M: Interrelationships of the 11 gasterosteiform families (sticklebacks, pipefishes, and their relatives): A new perspective based on whole mitogenome sequences from 75 higher teleosts. Molecular Phylogenetics and Evolution. 2007, 46: 224-236. 10.1016/j.ympev.2007.07.009.View ArticlePubMedGoogle Scholar
- Yagishita N, Miya M, Yamanoue Y, Shirai SM, Nakayama K, Suzuki N, Satoh TP, Mabuchi K, Nishida M, Nakabo T: Mitogenomic evaluation of the unique facial nerve pattern as a phylogenetic marker within the percifom fishes (Teleostei: Percomorpha). Molecular Phylogenetics and Evolution. 2009, 53: 258-266. 10.1016/j.ympev.2009.06.009.View ArticlePubMedGoogle Scholar
- Waters JM, Saruwatari T, Kobayashi T, Oohara I, McDowall RM, Wallis GP: Phylogenetic placement of retropinnid fishes: data set incongruence can be reduced by using asymmetric character state transformation costs. Systematic Biology. 2002, 51: 432-449. 10.1080/10635150290069887.View ArticlePubMedGoogle Scholar
- Caruso JH: The systematics and distribution of the lophiid anglerfishes: I. A revision of the genus Lophiodes with the description of two new species. Copeia. 1981, 1981: 522-549. 10.2307/1444556.View ArticleGoogle Scholar
- Bradbury MG: The genera of batfishes (family Ogcocephalidae). Copeia. 1967, 1967: 399-422. 10.2307/1442130.View ArticleGoogle Scholar
- Bertelsen E, Struhsaker PJ: The ceratioid fishes of the genus Thaumatichthys: Osteology, relationships, distribution, and biology. Galathea Report. 1977, 14: 7-40.Google Scholar
- Hedges SB, Sibley CG: Molecule versus morphology in avian evolution: The case of the "pelecaniform" birds. Proceedings of the National Academy of Sciences. 1994, 91: 9861-9865. 10.1073/pnas.91.21.9861.View ArticleGoogle Scholar
- Bertelsen E: The ceratioid fishes: ontogeny, taxonomy, distribution and biology. Dana Report. 1951, 39: 1-276.Google Scholar
- Tyler JC, Sorbini L: New superfamily and three new families of tetraodontiform fishes from the Upper Cretaceous: the earliest and most morphologically primitive plectognaths. Smithsonian Contributions to Paleontology. 1996, 82: 1-59.View ArticleGoogle Scholar
- Setiamarga DHE, Miya M, Yamanoue Y, Azuma Y, Inoue JG, Ishiguro NB, Mabuchi K, Nishida M: Divergence time of the two regional medaka populations in Japan as a new time scale for comparative genomics of vertebrates. Biology Letters. 2009, 5: 812-816. 10.1098/rsbl.2009.0419.PubMed CentralView ArticlePubMedGoogle Scholar
- Smith AG, Smith DG, Funnell BM: Atlas of Mesozoic and Cenozoic coastlines. 1994, Cambridge: Cambridge University PressGoogle Scholar
- Storey BC: The role of mantle plumes in continental breakup: case histories from Gondwanaland. Nature. 1995, 377: 301-308. 10.1038/377301a0.View ArticleGoogle Scholar
- Minegishi Y, Aoyama J, Inoue JG, Miya M, Nishida M, Tsukamoto K: Molecular phylogeny and evolution of the freshwater eels genus Anguilla based on the whole mitochondrial genome sequences. Molecular Phylogenetics and Evolution. 2005, 34: 134-146. 10.1016/j.ympev.2004.09.003.View ArticlePubMedGoogle Scholar
- Janvier P: Early vertebrates. 1996, Oxford: Oxford University PressGoogle Scholar
- Zhu M, Yu X, Wang W, Zhao W, Jia L: A primitive fish provides key characters bearing on deep osteichthyan phylogeny. Nature. 2006, 441: 77-80. 10.1038/nature04563.View ArticlePubMedGoogle Scholar
- Hurley IA, Mueller RL, Dunn KA, Schmidt EJ, Friedman M, Ho RK, Prince VE, Yang Z, Thomas MG, Coates MI: A new time-scale for ray-finned fish evolution. Proceedings of the Royal Society B. 2007, 274: 489-498. 10.1098/rspb.2006.3749.PubMed CentralView ArticlePubMedGoogle Scholar
- Patterson C: Osteichthyes: Teleostei. The Fossil Record 2. Edited by: Benton MJ. 1993, London: Chapman & Hall, 621-656.Google Scholar
- Wilson MVH, Brinkman DB, Neuman AG: Cretaceous Esocoidei (Teleostei): early radiation of the pikes in North American fresh waters. Journal of Paleontology. 1992, 66: 839-846.Google Scholar
- Benton MJ, Donoghue PC: Paleontological evidence to date the tree of life. Molecular Biology and Evolution. 2007, 24: 26-53. 10.1093/molbev/msl150.View ArticlePubMedGoogle Scholar
- Masters JC, de Wit MJ, Asher RJ: Reconciling the origins of Africa, India and Madagascar with vertebrate dispersal scenarios. Folia Primatologica. 2006, 77: 399-418. 10.1159/000095388.View ArticleGoogle Scholar
- Blot J: La faune ichthyologique des gisements du Monte Bolca (Province de Vérone, Italie). Bulletin Muséum national d'Histoire naturelle Paris (4). 1980, 2: 339-396.Google 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.