Evolutionary history of Otophysi (Teleostei), a major clade of the modern freshwater fishes: Pangaean origin and Mesozoic radiation
© Nakatani et al; licensee BioMed Central Ltd. 2011
Received: 29 January 2011
Accepted: 22 June 2011
Published: 22 June 2011
Freshwater harbors approximately 12,000 fish species accounting for 43% of the diversity of all modern fish. A single ancestral lineage evolved into about two-thirds of this enormous biodiversity (≈ 7900 spp.) and is currently distributed throughout the world's continents except Antarctica. Despite such remarkable species diversity and ubiquity, the evolutionary history of this major freshwater fish clade, Otophysi, remains largely unexplored. To gain insight into the history of otophysan diversification, we constructed a timetree based on whole mitogenome sequences across 110 species representing 55 of the 64 families.
Partitioned maximum likelihood analysis based on unambiguously aligned sequences (9923 bp) confidently recovered the monophyly of Otophysi and the two constituent subgroups (Cypriniformes and Characiphysi). The latter clade comprised three orders (Gymnotiformes, Characiformes, Siluriformes), and Gymnotiformes was sister to the latter two groups. One of the two suborders in Characiformes (Characoidei) was more closely related to Siluriformes than to its own suborder (Citharinoidei), rendering the characiforms paraphyletic. Although this novel relationship did not receive strong statistical support, it was supported by analyzing independent nuclear markers. A relaxed molecular clock Bayesian analysis of the divergence times and reconstruction of ancestral habitats on the timetree suggest a Pangaean origin and Mesozoic radiation of otophysans.
The present timetree demonstrates that survival of the ancestral lineages through the two consecutive mass extinctions on Pangaea, and subsequent radiations during the Jurassic through early Cretaceous shaped the modern familial diversity of otophysans. This evolutionary scenario is consistent with recent arguments based on biogeographic inferences and molecular divergence time estimates. No fossil otophysan, however, has been recorded before the Albian, the early Cretaceous 100-112 Ma, creating an over 100 million year time span without fossil evidence. This formidable ghost range partially reflects a genuine difference between the estimated ages of stem group origin (molecular divergence time) and crown group morphological diversification (fossil divergence time); the ghost range, however, would be filled with discoveries of older fossils that can be used as more reasonable time constraints as well as with developments of more realistic models that capture the rates of molecular sequences accurately.
Although freshwater lakes and rivers occupy a small portion of the Earth's surface (0.8%) and hold a negligible amount of the total water on Earth (0.01%), these ecosystems support an extraordinarily high proportion of the world's biodiversity, consisting of at least 100,000 species or nearly 6% of all described species . While this enormous biodiversity has been described through continued efforts by taxonomists, its origin and the history of the diversification on a global scale remain largely unexplored across diverse taxa. This is even true for well-studied taxa such as fishes, which account for the largest proportion of vertebrate diversity and is the taxon that controls the trophic structure of freshwater ecosystems .
A summary of character and taxon sampling in the previous molecular phylogenetic studies that include the four otophysan orders
Dimmick & Larson 
nc + mt rDNA
Ortí & Meyer 
Ortí & Meyer 
Saitoh et al. 
Lavoué et al. 
Peng et al. 
Li et al. 
10 nc genes
Poulsen et al. 
10 nc genes
In Otophysi, a great deal of attention has been paid to the biogeographic history that has shaped the current distribution patterns. Indeed, contrasting patterns in the geographic distributions of modern otophysans (Gondwanan vs. Laurasian vs. Pangaean distributions; Figure 1), the general acceptance of plate tectonics and continental drift as explanatory factors for dispersal across land connections and/or as a causal mechanism of vicariant speciation, and the presumed designation of the ostariophysans (Otophysi + Gonorynchiformes) as "primary freshwater fishes" (with dispersal across oceans being unlikely) together led many authors to propose evolutionary scenarios that attempt to identify a center of origin and dispersal routes through land connections based on various assumptions [4, 12, 27, 32–37]. Reconstruction of the history of otophysan diversification (e.g., the history of the modern familial diversification), however, remains unchallenged, apparently because of poor representation in the fossil record before the Cenozoic period [24, 38–41], a huge extant taxonomic diversity encompassing over 7943 species placed in 64 families and 1068 genera , and the absence of an adequate timescale for the phylogenies across major lineages (but see ). Indeed, in a review of the early radiation of teleosts, Arratia  stated "... the enormous radiation of some modern groups such as otophysans, atherinomorphs, perciforms, etc. is missing a historical framework."
To provide an overview of the history of modern otophysan diversification within the broad context of the evolutionary history of ray-finned fishes (Actinopterygii), we assembled whole mitochondrial genome (mitogenome) sequences from 66 otophysans (including 51 newly determined sequences), representing 55 of the 64 currently recognized families (86%; Figure 1). The 66 sequences were concatenated with those from 44 outgroup species for a total of 110 species and unambiguously aligned sequences (9923 bp excluding quickly saturated third codon positions) were subjected to phylogenetic analysis and a relaxed-clock Bayesian divergence time estimation. The resultant timetree suggests that the modern otophysan diversity has been shaped through the two consecutive mass extinction events on the Pangaea supercontinent and subsequent radiations during the Jurassic through the early Cretaceous.
List of the species used in this study
Scaphirhynchus cf. albus
Specimens and DNA extraction
A portion of epaxial musculature or pectoral fins (~0.25 g) from fresh specimens of each species was excised and the tissue was immediately preserved in 99.5% ethanol. Total genomic DNA from the ethanol-preserved tissue was extracted using DNeasy (Qiagen) and Aquapure genomic DNA isolation kit (Bio-Rad Laboratories, Inc.) in accordance with the respective manufacturer's protocols, or the standard phenol-chloroform method as described in Asahida et al. .
PCR and sequencing
Whole mitogenome sequences of the 51 otophysans (double asterisks in Table 2) were determined using a combination of long and short PCR methods developed by Miya and Nishida . Briefly, the mitogenomes of the 51 otophysans in their entirety were amplified using a long PCR technique  in two or three reactions. 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 sets of fish-versatile primers (and species-specific primers if necessary) were used for short PCRs to amplify contiguous overlapping segments of the entire mitogenome for each otophysan species. The short PCR products were purified using the Exosap-IT enzyme (GE Healthcare Bio-Sciences Corp.) and subsequently sequenced with dye-labeled terminators (BigDye terminator ver. 1.1/3.1; Applied Biosystems) and the primers used in the short PCRs. Sequencing reactions were conducted according to the manufacturer's instructions, followed by electrophoresis on an ABI Prism 377, 3100, or 3130 DNA sequencer (Applied Biosystems). A list of PCR primers used in this study is available from MNa upon request.
In some cases when multiple bands were amplified during short PCRs, we conducted subcloning using MinElute (Qiagen), pGEM-T Easy Vector Systems (Promega), and Z-competent E. coli (ZYMO Research), in accordance with the manufactures' protocols. To avoid PCR errors, we sequenced eight clones for each fragment using SP6 and T7 primers.
Sequence editing and alignment
Mitogenome sequences from the 66 otophysans were concatenated with the pre-aligned sequences used in Azuma et al.  in FASTA format and subjected to multiple alignment using MAFFT ver. 6.707 . The aligned sequences were imported into MacClade ver. 4.08  and the resulting gaps in the aligned sequences were manually removed to correctly reproduce the alignment used by Azuma et al. . All the resulting positions with gaps were removed, so the final data set consisted of 6904 positions from the first and second codon positions of the 12 protein-coding genes (excluding the ND6 gene because of its heterogeneous base composition and poor phylogenetic performance ), 1622 positions from the two rRNA genes, and 1397 positions from the 22 tRNA genes (total 9923 positions) (designated as 12nRTn: where 1, 2, R and T represent 1st codon position, 2nd codon position, rRNA gene and tRNA gene, respectively, and the subscript "n" denotes nucleotides). The third codon positions of the protein-coding genes were excluded from the data set because of the extremely high substitution rates (and the resulting multiple hits) and heterogeneous base composition as sources of systematic noise in phylogenetic analysis at this taxonomic level [49, 50] and overestimation of divergence time [51, 52]. The aligned sequences are available from TreeBase with the following URL (http://purl.org/phylo/treebase/phylows/study/TB2:S11469).
To investigate the relationships within the otophysans, we also created an additional four data sets that treated 12 protein-coding genes differently. The first three data sets considered only transversional changes in the first and/or third codon positions by converting purine (A/G) and pyrimidine (C/T) nucleotides to A and C, respectively (1r2nRnTn, 12n3rRTn, 1r2n3rRTn: where the subscript "r" denotes a modified RY-coding following Saitoh et al. ). The transitional changes in the first codon positions are somewhat saturated among distantly related taxa , and the first data set (1r2nRTn) was expected to reduce phylogenetic noise from the original data set (12nRTn). The second data set (12n3rRTn) added the RY-coded third codon positions to the original data set, which was expected to increase phylogenetic signals and was predominantly used to resolve interrelationships within the Cypriniformes [11, 12], one of the major otophysan clades (Figure 1). The third data set (1r2n3rRTn) removed transitional changes in the first codon positions from the second data set. The last data set converted protein-coding genes into amino acids (designated as 123aRTn) to explore the utility of these sequences in resolving otophysan taxa. Only 66 otophysans plus 12 outgroup species (4 gonorynchiforms + 3 clupeiforms + 3 alepocephaliforms + 2 anguillifoms) were used in these additional data sets to minimize the computation time.
Unambiguously aligned sequences were divided into three to five partitions depending on the data sets (three partitions in the 123aRTn data set, four partitions in the 12nRTn and 1r2nRTn data sets, and five partitions in the 12n3rRTn and 1r2n3rRTn data sets) and subjected to ML analysis. We used RAxML ver. 7.2.8  because it is the only ML-based software that can handle large data sets with data partitioning. A general time reversible model (GTR)  with sites following a discrete gamma distribution (Γ) and some sites invariable (I) was selected as the best model of nucleotide sequence evolution by Modeltest ver. 3.7  using the Akaike information criterion (AIC). For amino acid sequences, the MTREV model  with sites following a discrete gamma distribution (Γ) and some sites invariable (I) was used. We performed a rapid bootstrap (BS) analysis using this model (GTR + Γ + I) with 1000 replications (-f a option). This performs BS analysis using GTRCAT, which is a GTR approximation with optimization of individual per-site substitution rates and classification of these individual rates into a certain number of rate categories. After implementing the BS analysis, the program uses every fifth BS tree as a starting point for another ML search using the GTR + Γ + I 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 (BSPs) on the best-scoring ML tree.
Evaluation of alternative hypotheses
We manually created the constrained tree topologies with reference to the alternative hypotheses using MacClade and then performed RAxML analysis with each constraint using the -g option. We conducted fast bootstrapping 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 likelihood scores. Probabilities of alternative phylogenetic hypotheses were calculated using the likelihood-based approximately unbiased (AU) test  as implemented in CONSEL v.0.1k . The P-values from this test are calculated using the multi-scale bootstrap technique and are less biased than those of conventional methods  such as the BS probability (BSP) , the Kishino-Hasegawa (KH) test  and the Shimodaira-Hasegawa (SH) test .
Supplementary analysis using nuclear genes
DNA accession numbers of the nuclear genes from the four characiphysans
Divergence time estimation
Time constraints used for 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 (Paleogene) 
The †Parabarbus fossil (Cyprinidae) from the Ypresian (Paleogene) 
Gymnotiform fossil from the Late Miocene 
The Santanichthys (Characiformes) fossil from Albian (Cretaceous) 
The ariid fossil from Campanian (Cretaceous) 
The †Esteseox foxi fossil (Esociformes) from the Campanian (Cretaceous) 
The †Berycopsis fossil (Polymixiidae) from the Cenomanian (Cretaceous) 
The tetraodontiform fossil from the Cenomanian 
The estimated divergence time between Takifugu and Tetraodon 
U 95 L 85
The upper and lower bounds of separation between Madagascar and Indian 
U 145 L 112
The upper and lower bounds of separation between Indo-Madagascar landmass and Gondwanaland 
U 120 L 100
The upper and lower bounds of separation between African and South American landmasses 
The MCMCTREE program allows for minimum (lower) and maximum (upper) time constraints, and multiple calibration points have been argued to provide overall more realistic divergence time estimates . Therefore we sought to obtain optimal phylogenetic coverage of calibration points across our tree, although we could set maximum constraints based on the fossil records only for the six nodes (Table 4). Other than these six nodes, 17 additional nodes were reasonably chosen to constraint their minimum ages only (total 29 time constraints for 23 nodes; Table 4). A hard and softbound version of the program (MCMCTREE-HS) was used, so that probabilities of the true divergence time falling outside the minimum bounds are zero, but small but not zero for the maximum bounds . All time constraints are provided in units of 100 Ma (i.e., 1 = 100 Ma) because some of the model components in the Bayesian analysis are scale-variant . The calibration nodes with minimal bound only were set as L (tmin) and those with both minimal and maximal bounds were set as B (tmin, tmax). The former setting (L) assumes a heavy-tailed density (nearly a flat prior) based on a truncated Cauchy distribution of p = 0.1 and c = 1 as the default  ("standard minimum-age constraints" ). We did not manipulate the two shape parameters of the truncated Cauchy distribution because of insufficient information with which to specify meaningful prior distributions for most otophysan diversification times.
MCMC approximation with a burn-in period of 10,000 cycles was obtained, and every 50 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, we performed two replicate MCMC runs with two different random seeds for each analysis. MCMC samples from the two runs were combined after checking the distributions of parameter values using Tracer 1.5 (available from http://tree.bio.ed.ac.uk/software/tracer/). The number of samples (20,000) was large enough to reach effective sample sizes (ESSs >200) for all parameters estimated in this study.
To evaluate the effects of topological uncertainties on divergence time estimation, we also conducted dating analyses with four of the 14 alternative topologies (four topologies with the second best P value, the two worst P values, plus a hypothesis advocated by Fink and Fink ).
Tracing character evolution
The ancestral habitat was reconstructed on the timetree under a ML optimality criterion using Mesquite ver. 2.71 . The ML reconstruction methods found that the ancestral states maximizing the probability of the observed states would evolve under a stochastic model of evolution [74, 75]. The Mk1 model ("Markov k-state 1 parameter model"), a k-state generalization of the Jukes-Cantor model that corresponds to Lewis' Mk model , was used to trace the character evolution. Two character states were assigned to the terminal node: saltwater (character state 0) and freshwater (state 1).
Results and Discussion
Complete L-strand nucleotide sequences from the mitogenomes of the 51 species newly determined during this study were deposited in the DNA Data Bank of Japan (DDBJ), European Molecular Biology Laboratory (EMBL), and GenBank (Table 2). The genome content of the 51 species included two rRNA, 22 tRNA, and 13 protein-coding genes, plus the putative control region, as found in other vertebrates. Their gene arrangements were identical to the typical gene order of vertebrates.
Paraphyletic characiforms have been repeatedly recovered in previous molecular studies (Figure 2C-E) but in different manners. For example, Ortí and Meyer  assembled nuclear ependymin sequences (588 bp) from 22 otophysans to explore the phylogenetic utility of this gene and found that characoids were more closely related to Gymnotiformes than to a citharinoid (Distichodus) (Figure 2C). To investigate the characiform radiation, Ortí and Meyer  used partial mitochondrial rDNA sequences (870 bp) as phylogenetic markers and observed that two citharinoids (Distichodus and Citharinus) were more closely related to Siluriformes than to characoids (Figure 2D). Although analyses of whole mitogenome sequences by Saitoh et al.  (maximum parsimony analysis of 8096 bp) and Peng et al.  (6198 bp) did not include citharinoids, one of the two characoids were more closely related to either Siluriformes or Gymnotiformes than to its own order (Figure 2E). These relationships, however, received weak statistical support (<50% BSPs) and alternative tree topologies could not be rejected [4, 19, 29].
Results from AU-tests among 15 alternative tree topologies of the four major lineages derived from analysis of whole mitogenome sequences
Diff -ln L
1(Gym, (Cit, (Cha, Sil)))
2(Gym, (Sil, (Cit, Cha)))
3((Cit, Gym), (Cha, Sil))
4(Gym, (Cha, (Cit, Sil)))
5(Cha, (Sil, (Cit, Gym)))
6(Cit, (Gym, (Cha, Sil)))
7(Cit, (Cha, (Sil, Gym)))
8(Cha, (Cit, (Sil, Gym)))
9((Cit, Cha), (Sil, Gym)))
10(Cha, (Gym, (Cit, Sil)))
11(Sil, (Cha, (Cit, Gym)))
12(Sil, (Gym, (Cit, Cha)))
13(Cit, (Sil, (Cha, Gym)))
14((Cha, Gym), (Cit, Sil))
15(Sil, (Cit, (Cha, Gym)))
Results from AU-tests among 15 alternative tree topologies of the four major lineages derived from analysis of 10 nuclear gene sequences
Diff -ln L
1(Gym, (Cit, (Cha, Sil)))
2(Gym, (Cha, (Cit, Sil)))
3(Gym, (Sil, (Cit, Cha)))
4((Cit, Cha), (Sil, Gym)))
5(Cha, (Gym, (Cit, Sil)))
6(Cit, (Gym, (Cha, Sil)))
7((Cha, Gym), (Cit, Sil))
8(Cha, (Sil, (Cit, Gym)))
9(Sil, (Gym, (Cit, Cha)))
10(Cha, (Cit, (Sil, Gym)))
11(Cit, (Cha, (Sil, Gym)))
12(Sil, (Cit, (Cha, Gym)))
13(Cit, (Sil, (Cha, Gym)))
14((Cit, Gym), (Cha, Sil))
15(Sil, (Cha, (Cit, Gym)))
Interfamilial relationships within each of the four major lineages are outside the scope of this study because taxonomic sampling is still sparse and requires at least several distantly related genera from the same family to reconstruct more reliable trees. Nevertheless, brief comparisons with previous studies are useful to direct future research.
In Cypriniformes, monophyly of one of the two superfamilies (Cyprinoidea) was consistently recovered, while another superfamily (Cobitoidea) has been variously recovered as mono- or paraphyletic depending on the data treatment (most notably that of the third codon position) even within some of the same studies [11–17]. An enigmatic Psilorhynchus (placed in its own family Psilorhynchidae) is nested within the cyprinid lineage (Cyprinus + Danio; Figure 3), rendering Cyprinidae paraphyletic as reported previously [15, 17, 18].
In Gymnotiformes, a sister group relationship between two gymnotids (Electrophorus and Gymnotus) was confidently recovered (BSP = 90%), but other portions of the tree were weakly supported (BSPs ≤58%; Figure 3) and do not warrant further comment. Such ambiguous interfamilial relationships have been reported in previous studies based on mitochondrial rDNA sequences [23, 24].
In Citharinoidei, one of the two suborders of Characiformes, two distichodontid genera (Distichodus and Ichthyborus) formed a strongly supported monophyletic group (BSP = 100%), which was sister to a citharinid (Citharinus). These relationships have been consistently recovered in previous studies based on both morphological [79, 82] and molecular [19, 20, 83] analyses. Within Charachoidei, another suborder of Characiformes, basal relationships were poorly resolved and only four of the 19 internal branches that connect two species from the same families received the highest statistical support (100% BSPs; Figure 3). With the exception of such ambiguous basal relationships, the resulting phylogenies showed several similarities to those reported by Calcagnotto et al. , who analyzed four nuclear and two mitochondrial genes across 124 characiform taxa. In particular, we confirmed that the two African lineages represented by Hepsetidae (Hepsetus) and Alestidae (Phenacogrammus and Micralestes) did not form a monophyletic group, but were nested in the same clade with strictly Neotropical species, such as Hydrolycus (Cynodontidae) and Hoplias (Erythrinidae) (Figure 3).
In Siluriformes, monophyly of the two currently recognized suborders (Loricaroidei, Siluroidei ) was confirmed and the Loricaroidei was recovered as sister to other clades, followed by the divergence of Diplomystoidei and Siluroidei (BSP = 82%; Figure 3). These basal divergences are fully congruent with the results of a recent molecular phylogenetic study by Sullivan et al. , who analyzed two nuclear genes across 110 catfish species representing 36 of 37 families. Thus, the two different lines of evidence support the most basal Loricaroidei in the siluriform phylogenies. Note that previous morphological studies have consistently argued the most basal Diplomystoidei ([41, 84–89]), which was supported by a molecular study by Hardman , who used Diplomystes mesembrinus (Diplomystidae) as the only outgroup to root the tree. Within the Loricaroidei, a sister group relationship between Pterygoplichthys (Loricaridae) and Astroblepus (Astroblepidae) was strongly supported as in previous morphological and molecular studies [21, 89]. As in the Characoidei, basal relationships in the Siluroidei were poorly resolved and only five of the 23 internal branches received high statistical support (≥95% BSPs; Figure 3). Nevertheless we recovered some of the suprafamilial clades reported by Sullivan et al. , and list them below with genera used in this study (Figure 3): "Big Asia" (Pseudobagrus, Hara, Liobagrus), "Claroidea (Cla)" (Clarias, Heteropneustes), "Aspredinidae + Doradoidea (Asp + Dor)" (Bunocephalus, Amblydoras, Tatia, Tetranematichthys), "Ictaluroidea (Ict)" (Ictalurus, Cranoglanis), and "Big Africa" (Synodontis, Amphilius, Malapterurus, Auchenoglanis, Chrysichthys, Pareutropius).
Estimation of divergence time
MCMCTREE analyses of the divergence times based on the two data sets (12nRTn and 123 a ) with the assumption of independent rates (IR) across nodes yielded similar estimated node ages (see Additional file 1). Overall dating analysis based on the nucleotide data set (12nRTn) provided slightly older node ages (absolute differences in posterior means: 16.3 million years ± 21.5 SD) with consistently smaller 95% credible intervals (absolute differences: 30.3 million years ± 23.0 SD) than those from the amino acid data set (123 a ). Consequently the 95% credible intervals from both data sets greatly overlapped with each other across all nodes. These tendencies held true for the estimated node ages from the most recent common ancestors (MRCAs) of the five major otophysan lineages themselves (Cypriniformes, Gymnotiformes, Citharinoidei, Characoidei, Siluriformes) and those from MRCAs comprising subsets of these five lineages (total = nine nodes), with the absolute difference in posterior means of 17.8 million years ± 5.7 SD. Considering the long evolutionary history of actinopterygians (≈ 450 million years), the differences seem minor and the following description and discussion of the results were based on the nucleotide data set (12nRTn) for simplicity.
Comparisons of the divergence time estimates (Ma) for selected nodes from alternative tree topologies
Topology 1 (P= 0.702)
Topology 2 (P= 0.588)
Topology 14 (P= 0.009)
Topology 15 (P= 0.001)
Topology 9 (P= 0.140)
Comparisons of divergence time estimates (Ma) for basal actinopterygian nodes between the present and previous mitogenomic studies
This study (Topology 1)
Setiamarga et al.
Azuma et al.
Our age estimates, however, are remarkably older than those indicated in previous reports based on the fossil record [43, 68, 90–92]. For example, the oldest otocephalan fossil discovered to date was from the Late Jurassic 150 Ma , while our estimate was 265 Ma. The oldest otophysan fossil record dates back to the Albian (100-112 Ma) , while our estimate was 248 Ma. For all other major otophysan lineages, our molecular estimates (from 160 Ma in Cypriniformes to 220 Ma in Characiformes) are over 100 million years older than those of the fossil records used as minimum constraints (from 8 Ma in Gymnotiformes to 98 Ma in Characiformes; Table 4). Azuma et al.  also noted such large differences and plotted minimum time constraints based on the fossil record against molecular time estimates of the corresponding nodes. They found that four data points in the Paleozoic showed a fairly good 1:1 relationship, whereas other points in the Mesozoic were considerably below the line of a 1:1 relationship. They considered that these significant departures from the expected relationships for the Mesozoic fossils may reflect the fact that they do not really represent the oldest fossil for the corresponding lineages.
For such remarkably older molecular estimates, Benton and Ayala  noted four pervasive biases that may cause molecular dates to be too old: 1) the calibration dates that are too old based on previous molecular studies; 2) undetected rapidly evolving genes; 3) an ancestral polymorphism that was maintained through a long evolutionary period; and 4) asymmetric distributions of estimated times, with a constrained younger end but an unconstrained older end. As discussed by Azuma et al. , whose data set formed the basis of this study, the first factor is not the case in the present study because we did not use calibration dates based on previous molecular studies. The third factor would be applicable when the genomic regions used are under long-term balancing selection, but no mitochondrial genes have been reported to be under such selection. Neither the second factor nor the fourth factor are true in this study because quickly saturated third codon positions and the control region were excluded from the present analysis, and because each mitochondrial gene used here was tested and shown to perform well for dating vertebrate divergences . We agree with the arguments of Azuma et al.  as well as those of Benton and Ayala , who stated that "careful choice of genes may be a more appropriate strategy (than the larger data strategy), with a focus on long and fast-evolving (yet alignable) sequences" for reliable dating. Whole mitogenome sequences seem to accommodate such requirements.
Comparisons of divergence time estimates (Ma) for basal actinopterygian nodes between the present and previous mitogenomic studies
This study (Topology 1)
Peng et al.
Saitoh et al.
Santini et al.
As expected from the remarkably narrow 95% credible intervals in the estimated ages for Otocephala (149-153 Ma; Table 9), Santini et al.  set both minimum (149 Ma) and maximum (152 Ma) time constraints for the MRCA of Ostarioclupeomorpha (= Otocephala in this study). The upper bound (152 Ma) was chosen based on their strong belief that the MRCA of Otocephala was unlikely to be older than the oldest stem elopomorph (true eels and their relatives: †Anaethalion) from the late Kimmeridgian, Jurassic 152 Ma . Accordingly, they manipulated three priors of the lognormal distribution (offset = 149; mean = 0.1; SD = 0.6), which allowed only 5% of the MCMC samples to exceed the soft upper bound (152 Ma). This approach is called "phylogenetic bracketing" , which obtains not only minimum, but also maximum constraints on the timing of a branching event using the dates of the preceding and subsequent branching episodes [98, 99]. Donoghue and Benton  argued that phylogenetic bracketing may be problematic because it assumes that the branching events above and below a calibration more reliably capture the timing of the branching event in question than the estimated date of the calibration itself. For other clades within the Otocephala, Santini et al.  also set relatively narrow time constraints: Ostariophysi (125-149 Ma), Characiformes (68-100 Ma), and Siluriformes (73-83.5 Ma). Accordingly, their resulting estimates (Table 9) fit very well with the fossil records and little chance exists (<5%) that the MRCA of Otocephala is more than 152 Ma.
We recognize that the approach in Santini et al.  is superior in that it can lend more credence to the fossil record than the standard minimum-age constraints. Nevertheless, we did not provide upper (maximum) bounds for these nodes or manipulate two parameters of the truncated Cauchy distributions (= lognormal distributions in Santini et al. ) because insufficient information exists with which to specify meaningful prior distributions for most otophysan diversification times. Thus, the remarkable gaps between estimated ages from the two studies mostly reflect differences in priors on node ages (as the time constraints) rather than differences in software (MCMCTREE vs. BEAST) or data (mitogenomes vs. nuclear genes).
Some recent studies have argued that mitochondrial genes are likely to yield older node ages than nuclear genes. For example, Hurley et al.  investigated the relationships and divergence times of the basal actinopterygians (including teleosts) using four nuclear genes and whole mitogenomes. A relaxed molecular clock Bayesian analysis based on these two data sets provided node ages for the most recent common ancestor of teleosts with posterior means of 219 Ma (181-265 Ma) and 246 Ma (206-292 Ma) depending orthologues in the nuclear genes (only the former estimate shown in their figure 4) and 296 Ma (268-326 Ma) in the mitochondrial genome (data taken from their supplementary tables 3 and 7). Despite the overlap between these two results (95% credible intervals), Hurley et al.  concluded that the mitochondrial estimate of the most recent common ancestor of teleosts was 50-100 million years older than that of the nuclear genes. Based on this comparison, they further argued that the discrepancy between their nuclear and mitochondrial estimates may have been due to evolutionary rate differences between these two genomes. Note that their nuclear estimation did not include any species from the osteoglossomorphs (the putative most primitive teleosts; see Figure 3), which is likely to result in underestimation of the crown node age of teleosts with insufficient taxon sampling.
More recently, Brandley et al.  investigated intercontinental dispersal of Plestiodon (Eumeces) lizards based on analysis of a single mitochondrial gene (ND1) plus seven nuclear genes. They performed a relaxed molecular clock Bayesian divergence time estimation using unpartitioned and partitioned data sets across these genes and found that extreme saturation obscured the underlying rate of evolution in the mitochondrial gene, resulting in overestimation of the divergence times. Such overestimation in the mitochondrial gene was most pronounced in the unpartitioned data set and less so in the two partitioned data sets. Brandley et al.  agreed with Phillips , who demonstrated that genes or gene partitions that evolve at extremely high rates may accumulate so many hidden substitutions, making it difficult to estimate the underlying process that created the data. Although their arguments are convincing, comparisons made by Brandley et al.  were based on a single mitochondrial gene, including third codon positions. We used whole mitochondrial genomes with various substitution rates across genes , excluded quickly saturated third codon positions entirely, and partitioned the data sets.
Also the estimated ages may be too old simply because of the entire absence of a fossil record . The absence of fossils, however, should not be taken as evidence of absence, as discussed extensively by Diogo . In an article on the early radiation of teleosts, Arratia  stated "... we know almost nothing concerning the fossil record of most otophysans, of most living perciforms, atheriniforms, etc." Thus, to assume a lack of otophysan fossil record during the Mesozoic is natural if the group existed as stem groups during this period.
As Brown et al.  convincingly argued, although ample room for improvement exists on both sides of the "rock-clock" divide (e.g., accounting for ghost lineages in the fossil record and developing more realistic models of rate evolution for molecular genetic sequences), the consistent and conspicuous disagreement between these two sources of data more likely reflects a genuine difference between estimated ages of 1) stem group origins and 2) crown group morphological diversification, respectively.
Comments on trans-Atlantic clades
Although this study was not intended to address questions regarding historical biogeography of the otophysans, divergence times of some of the transcontinental sister group relationships, particularly those between Africa and South America in Characiformes and Siluriformes, warrant brief comments regarding their origins.
Based on a hypothesis presented by Backup (unpublished thesis ), Lundberg  examined the implications of applying a strict vicariance scenario to the three transcontinental sister group pairs in Characiformes and noted that most of the diversification among characiforms would have occurred prior to the continental breakup (i.e., on Gondwana). Our timetree corroborates this hypothesis because all of our divergence time estimates of the three transcontinental pairs (Citharinoidei vs. Characoidei = 220 Ma; two alestids Micralestes + Phenacogrammus vs. Hydrolycus = 162 Ma; Hepsetus vs. Hoplias = 160 Ma) fall well before the Gondwanan separation (100-120 Ma) . Note, however, that recent discoveries of marine or brackish characiforms (Santanichthys  and Sorbinicharax ) from the Cretaceous challenge such a simple vicariance hypothesis .
Recently Lundberg et al.  examined the relationships of a highly distinct freshwater catfish, Lacantunia enigmatica (Lacantuniidae) from Chipas, Mexico, and found that Lacantunia was derived from within a multifamily clade of African freshwater catfishes ("Big Africa" in Sullivan et al. ; see Figures 3, 6). Based on a maximum-age constraint of 144 Ma for the stem of the siluriform lineage (= siluriform + gymnotiform node) and seven additional minimum-age constraints for the internal nodes, they estimated the divergence time of Lacantunia using a relaxed-molecular clock method implemented in MULTIDIVTIME  and r8s . The resulting divergence estimates ranged from 83 to 87 Ma depending on the partitioning schemes with 95% credible intervals of 75-94 Ma. These ages fall after the separation of Africa and South America, and Lundberg et al.  proposed a new scenario based on an ancient intercontinental passage for explaining disjunct distributions of relictually endemic Lacantunia and its African sister clade. Also note that they acknowledge that their choice of 144 Ma as the maximum-age constraint is arbitrary, but conservatively informed by the fossil record of actinopterygians. Although we did not include Lacantunia in our analyses, the onset of diversification of the Big African clade is estimated to be 129 Ma (113-144 Ma) and divergence of Lacantunia from its stem should be much older. Thus our estimated ages did not require an ancient intercontinental passage as hypothesized by Lundberg et al. .
History of diversification
ML reconstruction of the ancestral habitats on the timetree suggests that a shift from character state 0 (saltwater) to 1 (freshwater) occurred between the crown node of Ostariophysi (P0 = 0.859; P1 = 0.142) and that of Otophysi (P0 = 0.055; P1 = 0.945), corresponding to the end-Permian (263-249 Ma) (Figure 6). Moreover, a similar habitat shift was also recovered in Gonorynchiformes, a sister clade of Otophysi (Figure 6). Based on these ancestral habitat reconstructions and the divergence time estimates with reference to Earth's history, we identified the following five global turning points in the history of otophysan diversification. Note, however, that the true ancestral habitats and dynamics of diversification cannot be inferred from time-calibrated phylogenies of the extant lineages alone ; one need to include the fossil record to fully understand the ancestral habitats and diversity dynamics of the otophysans.
First, the common ancestor of Otophysi presumably entered freshwater around 263-249 Ma. This period corresponds to the end-Permian when the largest mass extinction in Earth's history occurred, wiping out 96% of all species . The end-Permian mass extinction has been associated with a massive release of carbon gases into the atmosphere, causing a global greenhouse effect and abrupt climate warming . During the same time, the oceans are believed to have become anoxic worldwide [112, 113] and to have contained free hydrogen sulfide . These environmental perturbations greatly altered the marine ecosystems  and may have driven fish extinctions. Note that freshwater contains much higher levels of dissolved oxygen than saltwater for a given atmospheric concentration, and this difference can affect the biology of the organisms . The possible habitat shift (from salt- to freshwater) in the common ancestor of Otophysi should have a causal relationship with the survival of this lineage across the end-Permian. To our knowledge, however, no comparable examples of habitat shift in other groups of organisms exist that independently support this evolutionary scenario.
Second, the five major otophysan lineages (corresponding to common ancestors of the present-day orders or suborders) have successively originated from an ancestral lineage of the otophysans in a short time window during the Triassic (251-200 Ma). During this period all continents remained united as the supercontinent Pangaea, with a climate characterized by globally warm temperatures, extreme seasonality, and high aridity over much of the inland region . However, the climate was more moderate around the edges of the supercontinent. The regions around Pangaea had sufficient rainfall to produce vigorous forests along riverbanks , which possibly facilitated stepwise divergences of the five lineages on Pangaea. Note that the Pangaea was cut almost in half along the east-west axis by a huge embayment called the Tethys Seaway (Figure 6). Over time, as the Tethys Seaway expanded, it may have led to the first vicariant divergence between the common ancestors of Cypriniformes and Characiphysi as suggested by Saitoh et al. [4, 12].
Third, these five major otophysan lineages survived the end-Triassic mass extinction, which was one of the five largest extinctions in Earth's history in which 80% of all species became extinct . Widespread magmatic activity of the Central Atlantic Magmatic Province (CAMP) has been suggested to have caused this catastrophic event, and repeated release of SO2 gas, heavy metal emissions, and darkening were the main environmental stressors . Relatively long stem branches from the common ancestors of the five major lineages across the Triassic-Jurassic boundary suggest profound effects of the mass extinction and associated environmental perturbations on the patterns of otophysan diversification.
Fourth, extant otophysan lineages began to diversify in each of the five clades during the Jurassic and through the early Cretaceous (200-100 Ma) (Figure 6), providing a framework for the modern otophysan diversity. The onset of diversification depends on the clade (193-160 Ma), although all fell in the Jurassic (200-146 Ma). Note that familial diversification of the South American gymnotiforms, characoids, and loricaroids was recovered much earlier than that of cypriniforms, citharinoids and siluroids (Figure 6). The Jurassic was a time of particularly swift change due to the Pangaean breakup and the resulting development of new oceans and a gentle tropical climate over the formerly arid interiors of the supercontinent . Such environmental changes likely yielded numerous novel habitats, facilitating otophysan diversification throughout the Mesozoic until about 100 Ma.
Finally these modern otophysan lineages established during the Mesozoic survived the third mass extinction event at the end of the Cretaceous (Figure 6). Such "mass survival" of the modern lineages across the Cretaceous-Paleogene boundary has been noted for birds and mammals based on molecular evidence .
The present timetree suggests a Pangaean origin and Mesozoic radiations of the modern otophysans. This evolutionary scenario is in good agreement with recent biogeographic inferences. For example Diogo  broadly surveyed the higher-level phylogeny, biogeographic distribution, physiology, and ecology of catfishes, and suggested that these fishes originated at a time when some Pangaean connections still existed between Laurasia and Gondwana. In addition, Briggs  examined the phylogeny and geographic distribution patterns of ostariophysan fishes and proposed the Late Jurassic (160-150 Ma) origins of the cypriniforms and siluriforms, which are highly congruent with the estimated node ages for the most recent common ancestor of Cypriniformes (160 Ma) and Siluriformes (180 Ma) in this study.
Of course the evolutionary scenario presented here just represents a testable hypothesis and should be viewed with caution because the fossil record provides no direct evidence. Although Mesozoic freshwater deposits are geographically and stratigraphically patchy and the freshwater fishes from most of these deposits have not been fully documented across all continents [121–123], available information suggests that major components of the freshwater fishes were basal actinopterygians (non-teleost ray-finned fishes), chondrichthyans (sharks, rays, chimaeras) and sarcopterygians (lungfishes and coelacanths) [121–124], and few freshwater teleosts (e.g., the Late Triassic Jiangilichthys) are referable to extinct "pholidophoriforms," phylogenetically located outside the modern teleosts . Moreover the fossil record suggests that the modern teleostean lineages did not diversify until the Late Jurassic about 150 Ma . Thus, acceptance of the present evolutionary scenario requires the origin and survival of the ghost lineage from the Triassic through Late Jurassic for over 100 million years.
The timetree presented here indicates that survival of the ancestral lineages through the two consecutive mass extinctions on Pangaea and subsequent radiations during the Jurassic and through the early Cretaceous shaped the modern familial diversity of otophysans. The Pangaean origin and Mesozoic radiations of the modern otophysans are consistent with recent arguments based on biogeographic inferences [40, 41] and molecular divergence time estimates [37, 40]. No fossil otophysan, however, has been recorded before the Albian, the early Cretaceous 100-112 Ma , creating an over 100 million year time span without fossil evidence. This extremely large ghost range partially reflects a genuine difference between the estimated ages of stem group origin (molecular divergence time) and crown group morphological diversification (fossil divergence time) ; the ghost range, however, may be filled with future discoveries of older fossils that can be used as more reasonable time constraints as well as with the development of more realistic models that accurately capture the divergence rates of molecular sequences.
We sincerely thank J.G. Inoue, S. Kanao, S. Lavoué, R.M.S. Vega for providing materials, C. Oliveira for species' identification, Z. Yang for technical advice on MCMCTREE, and M. MacCampbell for reviewing an earlier version of the manuscript. G. Arratia kindly reviewed a portion of the revised manuscript, provided a large amount of information on the Mesozoic fishes, and pointed out inherent problems of our inferences from a paleontological view. Of course, the review does not imply her acceptance of our hypothesis on the otophysan diversification and any mistakes in logic or fact remain our responsibility. Our manuscript was improved following constructive and insightful comments from the three anonymous reviewers and the associate editor. Funding was provided by JSPS Research grants (17207007, 19207007, and 22370035).
- Dudgeon D, Arthington AH, Gessner MO, Kawabata ZI, Knowler DJ, Lévéque C, Naiman RJ, Prieur-Richard AH, Soto D, Stiassny MLJ: Freshwater biodiversity: importance, threats, status and conservation challenges. Biol Rev. 2005, 81 (02): 163-182.PubMed
- Schindler DE: Fish extinctions and ecosystem functioning in tropical ecosystems. Proc Natl Acad Sci USA. 2007, 104 (14): 5707-5708. 10.1073/pnas.0700426104.PubMedPubMed Central
- Nelson JS: Fishes of the world. 2006, Hoboken, NJ: John Wiley & Sons, 4
- Saitoh K, Miya M, Inoue JG, Ishiguro NB, Nishida M: Mitochondrial genomics of ostariophysan fishes: Perspectives on phylogeny and biogeography. J Mol Evol. 2003, 56 (4): 464-472. 10.1007/s00239-002-2417-y.PubMed
- Poulsen JY, Møller PDR, Lavoué S, Knudsen SW, Nishida M, Miya M: Higher and lower-level relationships of the deep-sea fish order Alepocephaliformes (Teleostei: Otocephala) inferred from whole mitogenome sequences. Biol J Linn Soc. 2009, 98 (4): 923-936. 10.1111/j.1095-8312.2009.01323.x.
- Lavoué S, Miya M, Inoue JG, Saitoh K, Ishiguro NB, Nishida M: Molecular systematics of the gonorynchiform fishes (Teleostei) based on whole mitogenome sequences: Implications for higher-level relationships within the Otocephala. Mol Phylogenet Evol. 2005, 37 (1): 165-177. 10.1016/j.ympev.2005.03.024.PubMed
- Lavoué S, Miya M, Poulsen JY, Møller PR, Nishida M: Monophyly, phylogenetic position and inter-familial relationships of the Alepocephaliformes (Teleostei) based on whole mitogenome sequences. Mol Phylogenet Evol. 2008, 47 (3): 1111-1121. 10.1016/j.ympev.2007.12.002.PubMed
- Lecointre G: Gonorynchiformes in the teleostean phylogeny: molecules and morphology used to investigate interrelationships of the Ostariophysi. Gonorynchiformes and ostariophysan relationships: a comprehensive review. Edited by: Grande T, Poyato-Ariza FJ, Diogo R. 2010, Enfield, NH: Science Publishers
- 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. Biol J Linn Soc. 2005, 85 (3): 289-306. 10.1111/j.1095-8312.2005.00483.x.
- Miya M, Takeshima H, Endo H, Ishiguro NB, Inoue JG, Mukai T, Satoh TP, Yamaguchi M, Kawaguchi A, Mabuchi K, et al: Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences. Mol Phylogenet Evol. 2003, 26 (1): 121-138. 10.1016/S1055-7903(02)00332-9.PubMed
- 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. J Mol Evol. 2006, 63 (6): 826-841. 10.1007/s00239-005-0293-y.PubMed
- Saitoh K, Sado T, Doosey MH, Bart HLJ, Inoue JG, Nishida M, Mayden RL, Nishida M, Miya M: Evidence from mitochondrial genomics supports the lower Mesozoic of South Asia as the time and place of basal divergence of cypriniform fishes (Actinopterygii: Ostariophysi). Zool J Linn Soc. 2011, 161: 633-662. 10.1111/j.1096-3642.2010.00651.x.
- Mayden RL, Tang KL, Wood RM, Chen W-J, Agnew MK, Conway KW, Yang L, Simons AM, Bart HL, Harris PM, et al: Inferring the Tree of Life of the order Cypriniformes, the earth's most diverse clade of freshwater fishes: implications of varied taxon and character sampling. J Systemat Evol. 2008, 46 (3): 424-438.
- Chen WJ, Miya M, Saitoh K, Mayden RL: Phylogenetic utility of two existing and four novel nuclear gene loci in reconstructing Tree of Life of ray-finned fishes: The order Cypriniformes (Ostariophysi) as a case study. Gene. 2008, 423 (2): 125-134. 10.1016/j.gene.2008.07.016.PubMed
- He S, Gu X, Mayden RL, Chen W-J, Conway KW, Chen Y: Phylogenetic position of the enigmatic genus Psilorhynchus (Ostariophysi: Cypriniformes): Evidence from the mitochondrial genome. Mol Phylogenet Evol. 2008, 47 (1): 419-425. 10.1016/j.ympev.2007.10.012.PubMed
- Mayden RL, Chen WJ, Bart HL, Doosey MH, Simons AM, Tang KL, Wood RM, Agnew MK, Yang L, Hirt MV, et al: Reconstructing the phylogenetic relationships of the earth's most diverse clade of freshwater fishes --order Cypriniformes (Actinopterygii: Ostariophysi): A case study using multiple nuclear loci and the mitochondrial genome. Mol Phylogenet Evol. 2009, 51 (3): 500-514. 10.1016/j.ympev.2008.12.015.PubMed
- Mayden RL, Chen WJ: The world's smallest vertebrate species of the Genus Paedocypris: A new family of freshwater fishes and the sister group to the world's most diverse clade of freshwater fishes (Teleostei: Cypriniformes). Mol Phylogenet Evol. 2010, 57 (1): 152-175. 10.1016/j.ympev.2010.04.008.PubMed
- Chen W-J, Mayden RL: Molecular systematics of the Cyprinoidea (Teleostei: Cypriniformes), the World's largest clade of freshwater fishes: further evidence from six nuclear genes. Mol Phylogenet Evol. 2009, 52: 544-549. 10.1016/j.ympev.2009.01.006.PubMed
- Ortí G, Meyer A: The radiation of characiform fishes and the limits of resolution of mitochondrial ribosomal DNA sequences. Syst Biol. 1997, 46 (1): 75-100. 10.1093/sysbio/46.1.75.PubMed
- Calcagnotto D, Schaefer SA, Desalle R: Relationships among characiform fishes inferred from analysis of nuclear and mitochondrial gene sequences. Mol Phylogenet Evol. 2005, 36 (1): 135-153. 10.1016/j.ympev.2005.01.004.PubMed
- Sullivan J, Lundberg J, Hardman M: A phylogenetic analysis of the major groups of catfishes (Teleostei: Siluriformes) using rag1 and rag2 nuclear gene sequences. Mol Phylogenet Evol. 2006, 41 (3): 636-662. 10.1016/j.ympev.2006.05.044.PubMed
- Hardman M: The phylogenetic relationships among non-diplomystid catfishes as inferred from mitochondrial cytochrome b sequences; the search for the ictalurid sister taxon (Otophysi: Siluriformes). Mol Phylogenet Evol. 2005, 37 (3): 700-720. 10.1016/j.ympev.2005.04.029.PubMed
- Alves-Gomes JA, Orti G, Haygood M, Heiligenberg W, Meyer A: Phylogenetic analysis of the South American electric fishes (order Gymnotiformes) and the evolution of their electrogenic system: a synthesis based on morphology, electrophysiology, and mitochondrial sequence data. Mol Biol Evol. 1995, 12: 298-318.PubMed
- Alves-Gomes JA: The mitochondrial phylogeny of the South American electric fish (Gymnotiformes) and an alternative hypothesis for the otophysan historical biogeography. Gonorynchiformes and ostriophysan relationships: A comprehensive review. Edited by: Grande T, Poyato-Ariza FJ, Diogo R. 2010, Enfield, NH: Science Publishers, 517-565.
- Hennig W: Phylogenetic systematics. 1996, Urbana, IL: University of Illinois Press
- Wiley EO: Phylogenetics: the theory and practice of phylogenetic systematics. 1981, New York: Wiley & Sons
- Fink SV, Fink WL: Interrelationships of the ostariophysan fishes (Teleostei). Zool J Linn Soc. 1981, 72 (4): 297-353. 10.1111/j.1096-3642.1981.tb01575.x.
- Dimmick WW, Larson A: A molecular and morphological perspective on the phylogenetic relationships of the otophysan fishes. Mol Phylogenet Evol. 1996, 6 (1): 120-133. 10.1006/mpev.1996.0064.PubMed
- Ortí G, Meyer A: Molecular evolution of ependymin and the phylogenetic resolution of early divergences among euteleost fishes. Mol Biol Evol. 1996, 13 (4): 556-573.PubMed
- Peng Z, He S, Wang J, Wang W, Diogo R: Mitochondrial molecular clocks and the origin of the major otocephalan clades (Pisces: Teleostei): a new insight. Gene. 2006, 370: 113-124.PubMed
- Li C, Lu G, Orti G: Optimal data partitioning and a test case for ray-finned fishes (Actinopterygii) based on ten nuclear loci. Syst Biol. 2008, 57 (4): 519-539. 10.1080/10635150802206883.PubMed
- Novacek MJ, Marshall LG: Early biogeographic history of ostariophysan fishes. Copeia. 1976, 1976 (1): 1-12. 10.2307/1443767.
- Briggs JC: The biogeography of otophysan fishes (Ostariophysi: Otophysi): a new appraisal. J Biogeogr. 2005, 32 (2): 287-294. 10.1111/j.1365-2699.2004.01170.x.
- Briggs JC: Ostariophysan zoogeography: an alternative hypothesis. Copeia. 1979, 1979: 111-118. 10.2307/1443735.
- Gayet M: Consideration sur la phylogenie et la paleobiographie des Ostariophysaries. Geobios. 1982, 6: 39-52.
- Banarescu P: Zoogeography of fresh waters. 1990, Wiesbaden, Germany: AULA-Verlag, 1:
- Lundberg JG: African-South American freshwater clades and continental drift: problems with a paradigm. Biological relationships between Africa and South America. Edited by: Goldblatt P. 1993, New Haven: Yale University Press, 156-199.
- Fara E, Gayet M, Taverne L: The Fossil record of Gonorynchiformes. Gonorynchiformes and ostariophysan relationships: A comprehensive review. Edited by: Grande T, Poyato-Ariza FJ, Diogo R. 2010, Enfield, NH: Science Publishers, 173-226.
- Malabarba MC, Malabarba LR: Biogeography of Characiformes: an evaluation of the available information of fossil and extant taxa. Origin and Phylogenetic Interrelationships of Teleosts. Edited by: Nelson JS, Schultze H-P, Wilson MVH. 2010, München: Verlag Dr. Friedrich Pfeil
- Conway KW, Hirt MV, Yang L, Mayden RL, Simons AM: Cypriniformes: systematics and paleontology. Origin and phylogenetic interrelationships of teleosts. Edited by: Nelson JS, Schultze H-P, Wilson MVH. 2010, München: Verlag Dr. Friedrich Pfeil
- Diogo R: Phylogeny, origin and biogeography of catfishes: support for a Pangean origin of 'modern teleosts' and reexamination of some Mesozoic Pangean connections between the Gondwanan and Laurasian supercontinents. Anim Biol. 2004, 54 (4): 331-351. 10.1163/1570756042729546.
- Arratia G: Mesozoic halecostomes and the early radiation of teleosts. Mesozoic fishes 3. Edited by: Arratia G, Tintori A. 2004, München, Germany: Verlag Dr. Friedrich Pfeil, 279-315.
- 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 Evol Biol. 2008, 8 (1): 215-10.1186/1471-2148-8-215.PubMedPubMed Central
- Asahida T, Kobayashi T, Saitoh K, Nakayama I: Tissue preservation and total DNA extraction from fish stored at ambient temperature using buffers containing high concentration of urea. Fish Sci. 1996, 62: 727-730.
- Miya M, Nishida M: Organization of the mitochondrial genome of a deep-sea fish, Gonostoma gracile (Teleostei: Stomiiformes): first example of transfer RNA gene rearrangements in bony fishes. Mar Biotechnol. 1999, 1 (5): 416-426. 10.1007/PL00011798.PubMed
- Cheng S, Chang SY, Gravitt P, Respess R: Long PCR. Nature. 1994, 369 (6482): 684-645. 10.1038/369684a0.PubMed
- Katoh K, Toh H: Recent developments in the MAFFT multiple sequence alignment program. Brief Bioinformatics. 2008, 9 (4): 286-298. 10.1093/bib/bbn013.PubMed
- Maddison WP, Maddison DR: MacClade Version 3. 2000, Sunderland: Sinauer Associates
- Miya M, Nishida M: Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion. Mol Phylogenet Evol. 2000, 17 (3): 437-455. 10.1006/mpev.2000.0839.PubMed
- Broughton RE: Phylogeny of teleosts based on mitochondrial genome sequences. Origin and phylogenetic interrelationships of teleosts. Edited by: Nelson JS, Shultze H-S, Wilson MVH. 2010, München, Germany: Verlag Dr. Friedrich Pfeil, 61-76.
- Benton MJ, Ayala FJ: Dating the Tree of Life. Science. 2003, 300 (5626): 1698-1700. 10.1126/science.1077795.PubMed
- Brandley MC, Wang Y, Guo X, de Oca ANM, Fería-Ortíz M, Hikida T, Ota H: Accommodating high rates of evolution in molecular dating methods: an example using inter-continental dispersal of Plestiodon (Eumeces) lizards. Syst Biol. 2011, 60 (1): 3-15. 10.1093/sysbio/syq045.PubMed
- Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.PubMed
- Yang Z: Estimating the pattern of nucleotide substitution. J Mol Evol. 1994, 39 (1): 105-111.PubMed
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.PubMed
- Adachi J, Hasegawa M: Model of amino acid substitution in proteins encoded by mitochondrial DNA. J Mol Evol. 1996, 42 (4): 459-468. 10.1007/BF02498640.PubMed
- Shimodaira H, Hasegawa M: CONSEL: for assessing the confidence of phylogenetic tree selection. Bioinformatics. 2001, 17 (12): 1246-1247. 10.1093/bioinformatics/17.12.1246.PubMed
- Shimodaira H: An approximately unbiased test of phylogenetic tree selection. Syst Biol. 2002, 51 (3): 492-508. 10.1080/10635150290069913.PubMed
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39 (4): 783-791. 10.2307/2408678.
- Kishino H, Hasegawa M: Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea. J Mol Evol. 1989, 29 (2): 170-179. 10.1007/BF02100115.PubMed
- Shimodaira H, Hasegawa M: Multiple comparisons of log-likelihoods with applications to phylogenetic inference. Mol Biol Evol. 1999, 16: 1114-1116.
- Kawahara R, Miya M, Mabuchi K, Near T, Nishida M: Stickleback phylogenies resolved: Evidence from mitochondrial genomes and 11 nuclear genes. Mol Phylogenet Evol. 2009, 50 (1): 401-404.PubMed
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.PubMed
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMedPubMed Central
- Rannala B, Yang Z: Inferring speciation times under an episodic molecular clock. Syst Biol. 2007, 56 (3): 453-466. 10.1080/10635150701420643.PubMed
- Zhong B, Yonezawa T, Zhong Y, Hasegawa M: Episodic evolution and adaptation of chloroplast genomes in ancestral grasses. PLoS ONE. 2009, 4 (4): e5297-10.1371/journal.pone.0005297.PubMedPubMed Central
- Yang Z, Rannala B: Bayesian phylogenetic inference using DNA sequences: a Markov chain Monte Carlo method. Mol Biol Evol. 1997, 14 (7): 717-724.PubMed
- Setiamarga DHE, Miya M, 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. Biol Lett. 2009, 5 (5): 81-86.
- Thorne JL, Kishino H: Divergence time and evolutionary rate estimation with multilocus data. Syst Biol. 2002, 51 (5): 689-702. 10.1080/10635150290102456.PubMed
- Benton MJ, Donoghue PCJ: Paleontological evidence to date the Tree of Life. Mol Biol Evol. 2006, 24 (1): 26-53. 10.1093/molbev/msl150.PubMed
- Yang Z, Rannala B: Bayesian estimation of species divergence times under a molecular clock using multiple fossil calibrations with soft bounds. Mol Biol Evol. 2006, 23 (1): 212-226.PubMed
- Brown JW, Rest JS, García-Moreno J, Sorenson MD, Mindell DP: Strong mitochondrial DNA support for a Cretaceous origin of modern avian lineages. BMC Biol. 2008, 6: 6-PubMedPubMed Central
- Mesquite: a modular system for evolutionary analysis Ver. 2.73. [http://mesquiteproject.org]
- Page M: The maximum likelihood approach to reconstructing ancestral character states of discrete characters on phylogenies. Syst Biol. 1999, 48 (3): 612-622. 10.1080/106351599260184.
- Schluter D, Price T, Mooers A, Ludwig D: Likelihood of ancestor states in adaptive radiation. Evolution. 1997, 51 (6): 1699-1711. 10.2307/2410994.
- Lewis PO: A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001, 50 (6): 913-925. 10.1080/106351501753462876.PubMed
- Lartillot N, Philippe H: A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. Mol Biol Evol. 2004, 21 (6): 1095-1109. 10.1093/molbev/msh112.PubMed
- Lartillot N, Lepage T, Blanquart S: PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics. 2009, 25 (17): 2286-2288. 10.1093/bioinformatics/btp368.PubMed
- Buckup PA: Relationships of the Characidiinae and phylogeny of characiform fishes (Teleostei: Ostariophysi). Phylogeny and classification of neotropical fishes. Edited by: Malabarba LR, Reis RE, Vari RP, Lucena ZMS, Lucena CAS. 1998, Porto Alegre, Brasil: Edipucrs, 123-144.
- Fink SV, Fink WL: Interrelationships of ostariophysan fishes (Teleostei). Interrelationships of fishes. Edited by: Stiassny MLJ, Parenti LR, Johnson GD. 1996, San Diego: Academic Press, 209-249.
- Dahdul WM: Review of the phylogenetic relationships and fossil record of Characiformes. Gonorynchiformes and ostriophysan relationships: A comprehensive review. Edited by: Grande T, Poyato-Ariza FJ, Diogo R. 2010, Enfield, New Hampshire: Science Publishers, 441-464.
- Vari RP: Anatomy, relationships and classification of the families Citharinidae and Distichodontidae (Pisces, Characoidea). Bull Brit Mus (Nat Hist) Zool Ser. 1979, 36: 261-344.
- Ortí G, Petry P, Porto JIR, Jégu M, Meyer A: Patterns of nucleotide change in mitochondrial ribosomal RNA genes and the phylogeny of piranhas. J Mol Evol. 1996, 42 (2): 169-182. 10.1007/BF02198843.PubMed
- Grande L: Redescription of Hypsidoris farsonensis (Teleostei: Siluriformes), with a reassessment of its phylogenetic relationships. J Vert Paleontol. 1987, 7: 24-54. 10.1080/02724634.1987.10011636.
- Mo T: Anatomy, relationships and systematics of the Bagridae (Teleostei: Siluroidei) with a hypothesis of siluroid phylogeny. Theses Zoology. 1991, 17: 1-216.
- Arratia G: Development and variation of the suspensorium of primitive catfishes (Teleostei: Ostariophysi) and their phylogenetic relationships. Bonn Zool Monogr. 1992, 32: 1-148.
- De Pinna MCC: Phylogenetic relationships of Neotropical Siluriformes: History, overview and synthesis of hypotheses. Phylogeny and Classification of Neotropical Fishes. Edited by: Malabarba LR, Reis RE, Vari RP, Lucena ZM, Lucena CAS. 1998, Porto Alegre, Brasil: Edipucrs, 279-330.
- He S, Gayet M, Meunier FJ: Phylogeny of the Amphiliidae (Teleostei: Siluriformes). Ann Sci Natur. 1999, 20 (4): 117-146.
- Diogo R, Peng Z: State of the art of siluriform higher-level phylogeny. Gonorynchiformes and ostriophysan relationships: A comprehensive review. Edited by: Grande T, Poyato-Ariza FJ, Diogo R. 2010, Enfield, New Hampshire: Science Publishers, 465-515.
- Inoue JG, Miya M, Venkatesh B, Nishida M: The mitochondrial genome of Indonesian coelacanth Latimeria menadoensis (Sarcopterygii: Coelacanthiformes) and divergence time estimation between the two coelacanths. Gene. 2005, 349: 227-235.PubMed
- Yamanoue Y, Miya M, Inoue JG, Matsuura K, Nishida M: The mitochondrial genome of spotted green pufferfish Tetraodon nigroviridis (Teleostei: Tetraodontiformes) and divergence time estimation among model organisms in fishes. Genes Genet Syst. 2006, 81 (1): 29-39. 10.1266/ggs.81.29.PubMed
- Inoue JG, Kumazawa Y, Miya M, Nishida M: The historical biogeography of the freshwater knifefishes using mitogenomic approaches: a Mesozoic origin of the Asian notopterids (Actinopterygii: Osteoglossomorpha). Mol Phylogenet Evol. 2009, 51 (3): 486-499. 10.1016/j.ympev.2009.01.020.PubMed
- Arratia G: Basal teleosts and teleostean phylogeny. Palaeo-Ichthyology. 1997, 7: 5-168.
- Filleul A, Maisey JG: Redescription of Santanichthys diasii (Otophysi, Characiformes) from the Albian of the Santana Formation and comments on its implications for otophysan relationships. Am Mus Novit. 2004, 3455: 1-21. 10.1206/0003-0082(2004)455<0001:ROSDOC>2.0.CO;2.
- Kumazawa Y, Azuma Y, Nishida M: Tempo of mitochondrial gene evolution: Can mitochondrial DNA be used to date old divergences?. Endocytobiosis Cell Research. 2004, 15: 136-142.
- Santini F, Harmon LJ, Carnevale G, Alfaro ME: Did genome duplication drive the origin of teleosts? A comparative study of diversification in ray-finned fishes. BMC Evol Biol. 2009, 9 (1): 194-10.1186/1471-2148-9-194.PubMedPubMed Central
- Donoghue PCJ, Benton MJ: Rocks and clocks: calibrating the Tree of Life using fossils and molecules. Trends Ecol Evol. 2007, 22 (8): 424-431. 10.1016/j.tree.2007.05.005.PubMed
- Reisz RR, Müller J: Molecular timescales and the fossil record: a paleontological perspective. Trends Genet. 2004, 20: 237-241. 10.1016/j.tig.2004.03.007.PubMed
- Müller J, Reisz RR: Four well-constrained calibration points from the vertebrate fossil record for molecular clock estimates. Bioessays. 2005, 1069-1075. 27PubMed
- Hurley I, Mueller R, Dunn K, Schmidt E, Friedman M, Ho R, Prince V, Yang Z, Thomas M, Coates M: A new time-scale for ray-finned fish evolution. Proc R Soc B. 2007, 274 (1609): 489-498. 10.1098/rspb.2006.3749.PubMedPubMed Central
- Phillips MJ: Branch-length estimation bias misleads molecular dating for a vertebrate mitochondrial phylogeny. Gene. 2009, 441 (1-2): 132-140. 10.1016/j.gene.2008.08.017.PubMed
- Miya M, Saitoh K, Wood R, Nishida M, Mayden RL: New primers for amplifying and sequencing the mitochondrial ND4/ND5 gene region of the Cypriniformes (Actinopterygii: Ostariophysi). Ichthyol Res. 2006, 53 (1): 75-81. 10.1007/s10228-005-0303-5.
- Buckup PA: The Characidiinae: a phylogenetic study of the South American darters and their relationships with other characiform fishes. Ph.D. thesis. 1991, Ann Harbor: The University of Michigan
- Storey BC: The role of mantle plumes in continental breakup: case histories from Gondwanaland. Nature. 1995, 377: 301-308. 10.1038/377301a0.
- Taverne L: Les poissons crétacés de Nardò. 16° Sorbinicharax verraesi gen. sp. nov. (Teleostei, Ostariophysi, Otophysi, Characiformes). Bollettino del Museo Civico di Storia Naturale di Verona. 2003, 27: 29-45.
- Otero O, Valentin X, Garcia G: Cretaceous characiform fishes (Teleostei: Ostariophysi) from Northern Tethys: description of new material from the Maastrichtian of Provence (Southern France) and palaeobiogeographic implications. Fishes and the break-up of Pangaea. Edited by: Cavin L, Longbottom A, Richter M. 2008, London: Geological Society, London, Special Publications, 295: 155-164.
- Lundberg JG, Sullivan JP, Rodiles-Hernández R, Hendrickson DA: Discovery of African roots for the Mesoamerican Chiapas catfish, Lacantunia enigmatica, requires an ancient intercontinental passage. Proc Acad Nat Sci Phila. 2007, 156 (1): 39-53. 10.1635/0097-3157(2007)156[39:DOARFT]2.0.CO;2.
- Sanderson MJ: r8s ver. 1.70. 2004
- Quental TB, Marshall CR: Diversity dynamics: molecular phylogenies need the fossil record. Trends Ecol Evol. 2010, 25 (8): 434-441. 10.1016/j.tree.2010.05.002.PubMed
- Benton MJ: Vertebrate palaeontology. 2005, Malden, MA: Blackwell, 3
- Benton MJ, Tverdokhlebov VP, Surkov MV: Ecosystem remodeling among vertebrates at the Permian-Triassic boundary in Russia. Nature. 2004, 432 (7013): 97-100. 10.1038/nature02950.PubMed
- Isozaki Y: Permo-Triassic boundary superanoxia and stratified superocean: records from lost deep sea. Science. 1997, 276 (5310): 235-238. 10.1126/science.276.5310.235.PubMed
- Sperling EA, Ingle JC: A Permian-Triassic boundary section at Quinn River Crossing, northwestern Nevada, and implications for the cause of the Early Triassic chert gap on the western Pangean margin. Bull Geol Soc Am. 2006, 118 (5-6): 733-746. 10.1130/B25803.1.
- Grice K, Cao C, Love GD, B ttcher ME, Twitchett RJ, Grosjean E, Summons RE, Turgeon SC, Dunning W, Jin Y: Photic zone euxinia during the Permian-Triassic superanoxic event. Science. 2005, 307 (5710): 706-709. 10.1126/science.1104323.PubMed
- Jin YG, Wang Y, Wang W, Shang QH, Cao CQ, Erwin DH: Pattern of marine mass extinction near the Permian-Triassic boundary in South China. Science. 2000, 289 (5478): 432-436. 10.1126/science.289.5478.432.PubMed
- Chapelle G, Peck LS: Polar gigantism dictated by oxygen availability. Nature. 1999, 399 (6732): 114-115.
- Huynh TT, Poulsen CJ: Rising atmospheric CO2 as a possible trigger for the end-Triassic mass extinction. Palaeogeogr, Palaeoclimatol, Palaeoecol. 2005, 217 (3-4): 223-242. 10.1016/j.palaeo.2004.12.004.
- Dixon D, Jenkins I, Moody RTJ, Zhuravlev AY: Atlas of the evolving earth. Volume 2. From the Devonian to the Cretaceous. 2001, New York: Macmillan Reference USA
- Guex J, Bartolini A, Atudorei V, Taylor D: High-resolution ammonite and carbon isotope stratigraphy across the Triassic-Jurassic boundary at New York Canyon (Nevada). Earth Planet Sci Lett. 2004, 225 (1-2): 29-41. 10.1016/j.epsl.2004.06.006.
- Cooper A, Penny D: Mass survival of birds across the Cretaceous-Tertiary boundary: molecular evidence. Science. 1997, 275: 1109-1113. 10.1126/science.275.5303.1109.PubMed
- López-Arbarello A: The record of Mesozoic fishes from Gondwana (excluding India and Madagascar). Mesozoic fishes 3. Edited by: Arratia G, Tintori A. 2004, München, Germany: Verlag Dr. Friedrich Pfeil, 597-624.
- Wilson MVH, Bruner JC, Arratia G, Tintori A: Mesozoic fish assemblages of North America. Mesozoic fishes 3. 2004, München, Germany: Verlag Dr. Friedrich Pfeil, 575-595.
- Chang MM, Miao D, Arratia G, Tintori A: An overview of Mesozoic fishes in Asia. Mesozoic fishes 3. Edited by: Arratia G, Tintori A. 2004, München, Germany: Verlag Dr. Friedrich Pfeil, 535-563.
- Milner ARC, Kirkland JI, Birthisel TA: The geographic distribution and biostratigraphy of Late Triassic-Early Jurassic freshwater fish faunas of the southwestern United States. New Mexico Museum of Natural History and Science Bulletin. 2006, 37: 522-529.
- Arratia G, Scasso R, Kiessling W: Late Jurassic fishes from Longing Gap, Antarctic Peninsula. J Vert Paleontol. 2004, 24 (1): 41-55. 10.1671/1952-4.
- Smith AG, Smith DG, Funnell BM: Atlas of Mesozoic and Cenozoic coastlines. 2004, Cambridge, UK: Cambridge University Press
- Ziegler A, Eshel G, Rees PMA, Rothfus T, Rowley D, Sunderlin D: Tracing the tropics across land and sea: Permian to present. Lethaia. 2003, 36 (3): 227-254. 10.1080/00241160310004657.
- Janvier P: Early vertebrates. 1996, Oxford, UK: Oxford University Press
- Zhu M, Yu X, Wang W, Zhao W, Jia L: A primitive fish provides key characters bearing on deep osteichthyan phylogeny. Nature. 2006, 441 (7089): 77-80. 10.1038/nature04563.PubMed
- Patterson C: Osteichthyes: Teleostei. The fossil record 2. Edited by: Benton MJ. 1993, London, UK: Chapman & Hall, 621-656.
- Gayet M, Meunier FJ: Premiére découverte de Gymnotiformes fossiles (Pisces, Ostariophysi) dans le Miocène supérieur de Bolivie. CR Acad Sci Paris. 1991, 313 (4): 471-476.
- Gayet M, Meunier FJ: Paleontology and palaeobiogeography of catfishes. Catfishes. Edited by: Arratia G, Kapoor BG, Chardon M, Diogo R. 2003, Enfield, NH: Science Publishers, 491-522.
- Wilson MVH, Brinkman DB, Neuman AG: Cretaceous Esocoidei (Teleostei): early radiation of the pikes in North American fresh waters. J Paleontol. 1992, 66: 839-846.
- Tyler JC, Sorbini L: New superfamily and three new families of tetraodontiform fishes from the Upper Cretaceous: the earliest and most morphologically primitive plectognaths. Smith Contr Paleobio. 1996, 82: 1-59.
- Benton MJ, Donoghue PCJ, Asher RJ: Calibrating and constraining molecular clocks. Timetree of life. Edited by: Hedges SB, Kumar S. 2009, Oxford, UK: Oxford University Press, 35-86.
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.