Use of mitogenomic information in teleostean molecular phylogenetics: a tree-based exploration under the maximum-parsimony optimality criterion

The complete mitochondrial genome of the whiting, Merlangius merlangus and the haddock, Melanogrammus aeglefinus: a detailed genomic comparison among closely related species of the Gadidae family

We determined the first complete mitochondrial DNA (mtDNA) sequences for the whiting (Merlangius merlangus, family Gadidae, order Gadiformes) and the haddock (Melanogrammus aeglefinus, family Gadidae, order Gadiformes). The entire mitogenomes were amplified and sequenced by primer walking using newly designed specific internal primers. Lengths were 16,569 and 16,585 bases for whiting and haddock respectively, lengths which lie within the range of previously reported gadiform sequences from Atlantic cod (Gadus morhua, 16,696 bases) and walleye pollock (Theragra chalcogramma, 16,570 bases). Gene arrangement in both species conformed to the order seen in most vertebrate mitochondrial genomes. We identified a long intergenic spacer located between the tRNA(Thr) and tRNA(Pro) genes (of 100 and 70 bp long for whiting and haddock, respectively), as previously described for other species of the order Gadiformes. Using nucleotide and amino acid divergence data of four complete gadoid mitogenomes (M. merlangius, M. aeglefinus, G. morhua and T. chalcogramma), we examined in detail the relative mtDNA mutation patterns across genes and among Gadidae species and tested for the performance of each protein-coding, transfer RNA and ribosomal RNA gene in depicting the expected phylogeny among the four species, as compared with the whole genome dataset. This comparison may be particularly useful in phylogenetic analyses of such a diverse fish family, as well as for the understanding of the patterns of nucleotide substitution of the mtDNA at low levels of divergence.

Round and pointed-head grenadier fishes (Actinopterygii: Gadiformes) represent a single sister group: Evidence from the complete mitochondrial genome sequences

The gene order of mitochondrial genomes (mitogenomes) has been employed as a useful phylogenetic marker in various metazoan animals, because it may represent uniquely derived characters shared by members of monophyletic groups. During the course of molecular phylogenetic studies of the order Gadiformes (cods and their relatives) based on whole mitogenome sequences, we found that two deep-sea grenadiers (Squalogadus modificatus and Trachyrincus murrayi: family Macrouridae) revealed a unusually identical gene order (translocation of the tRNA(Leu (UUR))). Both are members of the same family, although their external morphologies differed so greatly (e.g., round vs. pointed head) that they have been placed in different subfamilies Macrouroidinae and Trachyrincinae, respectively. Additionally, we determined the whole mitogenome sequences of two other species, Bathygadus antrodes and Ventrifossa garmani, representing a total of four subfamilies currently recognized within Macrouridae. The latter two species also exhibited gene rearrangements, resulting in a total of three different patterns of unique gene order being observed in the four subfamilies. Partitioned Bayesian analysis was conducted using available whole mitogenome sequences from five macrourids plus five outgroups. The resultant trees clearly indicated that S. modificatus and T. murrayi formed a monophyletic group, having a sister relationship to other macrourids. Thus, monophyly of the two species with disparate head morphologies was corroborated by two different lines of evidence (nucleotide sequences and gene order). The overall topology of the present tree differed from any of the previously proposed, morphology-based phylogenetic hypotheses.

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.

Production and utilization of a high-density oligonucleotide microarray in channel catfish, Ictalurus punctatus

Background

Functional analysis of the catfish genome will be useful for the identification of genes controlling traits of economic importance, especially innate disease resistance. However, this species lacks a platform for global gene expression profiling, so we designed a first generation high-density oligonucleotide microarray platform based on channel catfish EST sequences. This platform was used to profile gene expression in catfish spleens 2 h, 4 h, 8 h and 24 h after injection of lipopolysaccharide (LPS).

Results

In the spleen samples, 138 genes were significantly induced or repressed greater than 2-fold by LPS treatment. Real-time RT-PCR was used to verify the microarray results for nine selected genes representing different expression levels. The results from real-time RT-PCR were positively correlated (R² = 0.87) with the results from the microarray.

Conclusion

The first generation channel catfish microarray provided several candidate genes useful for further evaluation of immune response mechanisms in this species. This research will help us to better understand recognition of LPS by host cells and the LPS-signalling pathway in fish.

Spatial Covariation of Mutation and Nonsynonymous Substitution Rates in Vertebrate Mitochondrial Genomes

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.

Novel relationships among ten fish model species revealed based on a phylogenomic analysis using ESTs

The power of comparative phylogenomic analyses also depends on the amount of data that are included in such studies. We used expressed sequence tags (ESTs) from fish model species as a proof of principle approach in order to test the reliability of using ESTs for phylogenetic inference. As expected, the robustness increases with the amount of sequences. Although some progress has been made in the elucidation of the phylogeny of teleosts, relationships among the main lineages of the derived fish (Euteleostei) remain poorly defined and are still debated. We performed a phylogenomic analysis of a set of 42 of orthologous genes from 10 available fish model systems from seven different orders (Salmoniformes, Siluriformes, Cypriniformes, Tetraodontiformes, Cyprinodontiformes, Beloniformes, and Perciformes) of euteleostean fish to estimate divergence times and evolutionary relationships among those lineages. All 10 fish species serve as models for developmental, aquaculture, genomic, and comparative genetic studies. The phylogenetic signal and the strength of the contribution of each of the 42 orthologous genes were estimated with randomly chosen data subsets. Our study revealed a molecular phylogeny of higher-level relationships of derived teleosts, which indicates that the use of multiple genes produces robust phylogenies, a finding that is expected to apply to other phylogenetic issues among distantly related taxa. Our phylogenomic analyses confirm that the euteleostean superorders Ostariophysi and Acanthopterygii are monophyletic and the Protacanthopterygii and Ostariophysi are sister clades. In addition, and contrary to the traditional phylogenetic hypothesis, our analyses determine that killifish (Cyprinodontiformes), medaka (Beloniformes), and cichlids (Perciformes) appear to be more closely related to each other than either of them is to pufferfish (Tetraodontiformes). All 10 lineages split before or during the fragmentation of the supercontinent Pangea in the Jurassic.

The complete mitochondrial genome of the helmet catfish Cranoglanis bouderius (Siluriformes: Cranoglanididae) and the phylogeny of otophysan fishes

The complete sequence of the 16,539 nucleotide mitochondrial genome from the single species of the catfish family Cranoglanididae, the helmet catfish Cranoglanis bouderius, was determined using the long and accurate polymerase chain reaction (LA PCR) method. The nucleotide sequences of C. bouderius mitochondrial DNA have been compared with those of three other catfish species in the same order. The contents of the C. bouderius mitochondrial genome are 13 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, and a non-coding control region, the gene order of which is identical to that observed in most other vertebrates. Phylogenetic analyses for 13 otophysan fishes were performed using Bayesian method based on the concatenated mtDNA protein-coding gene sequence and the individual protein-coding gene sequence data set. The competing otophysan topologies were then tested by using the approximately unbiased test, the Kishino-Hasegawa test, and the Shimodaira-Hasegawa test. The results show that the grouping ((((Characiformes, Gymnotiformes), Siluriformes), Cypriniformes), outgroup) is the most likely but there is no significant difference between this one and the other alternative hypotheses. In addition, the phylogenetic placement of the family Cranoglanididae among siluriform families was also discussed.

Molecular phylogeny of the stromateoid fishes (Teleostei: Perciformes) inferred from mitochondrial DNA sequences and compared with morphology-based hypotheses

The phylogenetic relationships among 21 species of stromateoid fishes, representing five families and 13 genera, were reconstructed using 3263bp of mitochondrial DNA sequences, including the posterior half of the 16S rRNA and entire COI and Cytb genes. The resultant molecular phylogenies were compared with previous phylogenetic hypotheses inferred from morphological characters. Molecular phylogenetic trees were constructed using the maximum parsimony, maximum likelihood, and Bayesian methods. All three methods resulted in well-resolved trees with most nodes being supported by moderate to high support values. In contrast to previous morphological analyses, which resulted in non-monophyly of Centrolophidae, all three methods utilized for the present molecular analyses supported the monophyly of Centrolophidae, as well as the reciprocal monophyly of the other stromateoid families, previous morphological hypotheses being rejected by the Templeton and Shimodaira-Hasegawa tests. In addition, the three methods indicated a sister-group relationship between Ariommatidae and Nomeidae. The position of Tetragonuridae was, however, incongruent between the MP method and the ML and Bayesian methods, being placed in the most basal position of Stromateoidei in the former, but occupying a sister relationship to Stromateidae in the latter. Comparison of the molecular phylogenies to previous morphological hypotheses suggested that evolutionary changes in morphological characters have not occurred equally among the stromateoid lineages, the evolution of the centrolophids not having been accompanied by appreciable morphological changes, whereas other stromateoids have undergone considerable morphological changes during their evolutionary history. The molecular phylogenies also shed some light on the evolutionary pattern of the pharyngeal sac, two of the four types of sac corresponding to two main lineages of Stromateoidei. Some taxonomic implications were also discussed.

Evolutionary Rates, Divergence Dates, and the Performance of Mitochondrial Genes in Bayesian Phylogenetic Analysis

The mitochondrial genome is one of the most frequently used loci in phylogenetic and phylogeographic analyses, and it is becoming increasingly possible to sequence and analyze this genome in its entirety from diverse taxa. However, sequencing the entire genome is not always desirable or feasible. Which genes should be selected to best infer the evolutionary history of the mitochondria within a group of organisms, and what properties of a gene determine its phylogenetic performance? The current study addresses these questions in a Bayesian phylogenetic framework with reference to a phylogeny of plethodontid and related salamanders derived from 27 complete mitochondrial genomes; this topology is corroborated by nuclear DNA and morphological data. Evolutionary rates for each mitochondrial gene and divergence dates for all nodes in the plethodontid mitochondrial genome phylogeny were estimated in both Bayesian and maximum likelihood frameworks using multiple fossil calibrations, multiple data partitions, and a clock-independent approach. Bayesian analyses of individual genes were performed, and the resulting trees compared against the reference topology. Ordinal logistic regression analysis of molecular evolution rate, gene length, and the G-shape parameter a demonstrated that slower rate of evolution and longer gene length both increased the probability that a gene would perform well phylogenetically. Estimated rates of molecular evolution vary 84-fold among different mitochondrial genes and different salamander lineages, and mean rates among genes vary 15-fold. Despite having conserved amino acid sequences, cox1, cox2, cox3, and cob have the fastest mean rates of nucleotide substitution, and the greatest variation in rates, whereas rrnS and rrnL have the slowest rates. Reasons underlying this rate variation are discussed, as is the extensive rate variation in cox1 in light of its proposed role in DNA barcoding.

Mitochondrial molecular clocks and the origin of the major Otocephalan clades (Pisces: Teleostei): A new insight

The Otocephala, a clade including ostariophysan and clupeomorph teleosts, represents about a quarter of total fish species diversity, with about 1000 genera and more than 7000 species. A series of recent papers have defended that the origin of this clade and of its major groups may be significantly older than the oldest fossils of each of these groups suggest. Some of these recent papers explicitly defend a Pangean origin for some otocephalan groups such as the Siluriformes or Cypriniformes. To know whether or not the otocephalans as a whole, and particularly the mainly freshwater, cosmopolitan otophysans could have originated before the splitting of the Pangean supercontinent is of extreme importance, since otophysan fishes are among the most useful animal groups for the determination of historical continental relationships. In the present work we examined divergence times for each major otocephalan group by an analysis of complete mtDNA sequences, in order to investigate if these divergence times support the hypotheses advanced in recent studies. The complete mtDNA sequences of nine representative non-otocephalan fish species and of twenty-one representative otocephalan species was compared. The present study is thus, among the studies dealing with molecular divergence times of teleosts, the one in which a greater number of otocephalan species are included. The divergence times obtained support that the major otocephalan groups had a much older origin than the oldest fossil records available for these groups suggest. The origin of the Otocephala is estimated as having occurred about 282 Mya, with the origin of the Otophysi being estimated at about 251 Mya.

Molecular phylogeny and stripe pattern evolution in the cardinalfish genus Apogon

Cardinalfishes of the genus Apogon (Apogonidae) are one of the most speciose (>200 species) and numerically dominant fishes in coral reefs. Although the genus is divided into 10 subgenera, more than 70% of the species are included in the subgenus Ostorhinchus, most having either horizontal or vertical lines on the body. The phylogenetic relationship among 32 species of subgenus Ostorhinchus and 11 species of four other subgenera of Apogon, based on mitochondrially encoded 12S and 16S ribosomal genes and intervening tRNA(Val) gene, were investigated, using two species of the apogonid genus Fowleria as outgroups. The analyses demonstrated that Ostorhinchus (the most speciose subgenus) was polyphyletic, comprising at least three lineages, Ostorhinchus I, II, and III. Ostorhinchus I included two species, A. (O.) amboinensis and A. (O.) sangiensis, being a sister group to subgenus Zoramia. Ostorhinchus II and III included species with horizontal and vertical lines on the body, respectively. The respective monophylies of the latter two groups, together with a molecular clock calibration, indicated that in the evolutionary history of the genus, basic stripe patterns evolved first (more than 20 million years BP), with subsequent pattern diversification and modification.

Molecular systematics of the gonorynchiform fishes (Teleostei) based on whole mitogenome sequences: Implications for higher-level relationships within the Otocephala

Although the order Gonorynchiformes includes only 31 species assigned to seven genera and four families, it exhibits a large variety of anatomical structures, making difficult the reconstruction of phylogenetic relationships among its representatives. Within the basal teleosts, the Gonorynchiformes belong to the Otocephala where they have been alternatively placed as the sister group of the Otophysi and of the Clupeiformes. In this context, we investigated the phylogeny of the Gonorynchiformes using whole mitogenome sequences from 40 species (six being newly determined for this study). Our taxonomic sampling included at least one species of each gonorynchiform genus and of each other major otocephalan lineage. Unambiguously aligned, concatenated mitogenomic sequences (excluding the ND6 gene and control region) were divided into five partitions (1st, 2nd, and 3rd codon positions, tRNA genes, and rRNA genes) and partitioned Bayesian analyses were conducted. The resultant phylogenetic trees were fully resolved, with most of the nodes well supported by the high posterior probabilities. As expected, the Otocephala were recovered as monophyletic. Within this group, the mitogenome data supported the monophyly of Alepocephaloidei, Gonorynchiformes, Otophysi, and Clupeiformes. The Gonorynchiformes and the Otophysi formed a sister group, rending the Ostariophysi monophyletic. This result conflicts with previous mitogenomic phylogenetic studies, in which a sister relationship was found between Clupeiformes and Gonorynchiformes. We discussed the possible causes of this incongruence. Within the Gonorynchiformes, the following original topology was found: (Gonorynchus (Chanos (Phractolaemus (Cromeria (Grasseichthys (Kneria, Parakneria)))))). We confirmed that the paedomorphic species Cromeria nilotica and Grasseichthys gabonensis belong to the family Kneriidae; however, the two species together did not form a monophyletic group. This result challenges the value of reductive or absent characters as synapomorphies in this group.

The mitochondrial genome of Indonesian coelacanth Latimeria menadoensis (Sarcopterygii: Coelacanthiformes) and divergence time estimation between the two coelacanths

We determined the whole mitochondrial genome sequence for Indonesian coelacanth Latimeria menadoensis. The genome content and organization were identical to that of typical vertebrates including Comoran coelacanth, Latimeria chalumnae. The overall nucleotide differences between the two species (excluding the control region) was 4.28%. The divergence time between the two species was estimated using whole mitochondrial genome data from the two coelacanths and 26 actinopterygians that represent major actinopterygian lineages plus an outgroup. Partitioned Bayesian analyses were conducted with the two data sets that comprised concatenated amino acid sequences from 12 protein-coding genes (excluding ND6 gene) and concatenated nucleotide sequences from 12 protein-coding genes (without 3rd codon positions), 22 transfer RNA genes, and two ribosomal RNA genes. The molecular clock analysis was also conducted with the concatenated amino acid sequences from the 12 protein-coding genes after removing faster or more slowly evolving sequences. Using the sarcopterygian-actinopterygian split as a calibration point (450 Mya), divergence time estimation between L. menadoensis and L. chalumnae fell in the range of 40-30 Mya, which is much older than those of the previous studies (<6.3 Mya). Assuming that the most recent ancestor of Latimeria was distributed continuously along the deep coasts of Africa through Eurasia, our estimate is in agreement with the hypothesis that the collision of India with Eurasia (50 Mya) and the subsequent siltation caused by the formation of major rivers resulted in a coelacanth habitat disjunction that allowed populations on either side of India to diverge.

Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish

For many genes, ray-finned fish (Actinopterygii) have two paralogous copies, where only one ortholog is present in tetrapods. The discovery of an additional, almost-complete set of Hox clusters in teleosts (zebrafish, pufferfish, medaka, and cichlid) but not in basal actinopterygian lineages ( Polypterus) led to the formulation of the fish-specific genome duplication hypothesis. The phylogenetic timing of this genome duplication during the evolution of ray-finned fish is unknown, since only a few species of basal fish lineages have been investigated so far. In this study, three nuclear genes ( fzd8, sox11, tyrosinase) were sequenced from sturgeons (Acipenseriformes), gars (Semionotiformes), bony tongues (Osteoglossomorpha), and a tenpounder (Elopomorpha). For these three genes, two copies have been described previously teleosts (e.g., zebrafish, pufferfish), but only one orthologous copy is found in tetrapods. Individual gene trees for these three genes and a concatenated dataset support the hypothesis that the fish-specific genome duplication event took place after the split of the Acipenseriformes and the Semionotiformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. If these three genes were duplicated during the proposed fish-specific genome duplication event, then this event separates the species-poor early-branching lineages from the species-rich teleost lineage. The additional number of genes resulting from this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish.

Evolution of the Mitochondrial Genome in Cephalochordata as Inferred from Complete Nucleotide Sequences from Two Epigonichthys Species

Complete mitochondrial (mt) DNA sequences of two lancelets, Epigonichthys maldivensis and E. lucayanus, were compared with those of two Branchiostoma lancelets and several deuterostomes previously surveyed. The mt-gene order of E. lucayanus was quite different from that of E. maldivensis, the latter being identical to the two Branchiostoma species. A remarkable genomic change in E. lucayanus mtDNA was an inversion, indicating the possibility of recombination of the mt-genome. Gene rearrangements, probably attributable to tandem genome duplications and subsequent random deletions, were observed in two parts. Short major unassignable sequences of the examined lancelets were regarded as a part of putative regulative elements, judging from some sequence similarity to the conserved sequence block (CSB) in mammalian mtDNA. The considerable mt-genome reorganization in E. lucayanus seemed to have affected the nucleotide substitution pattern, suggested by base composition analyses. The present analysis also suggested that AGR codons in lancelet mtDNA were likely to correspond to serine residue, rather than glycine. Furthermore, the AGG codon, so far reputed to be unassignable in lancelet mtDNA, was found twice in E. maldivensis, indicating the availability of all four AGN codons in some lancelets. This finding lends support to an alternative hypothesis regarding the evolutionary history of AGR-codon assignment in extant chordates, rather than that previously proposed. A molecular phylogenetic tree of the Epigonichthys and Branchiostoma species based on DNA sequences of the 13 mt-protein genes doubted the monophyly of the former genus, unlike the prevailing classification based on their different gonadal arrangements.

Gene rearrangements and evolution of tRNA pseudogenes in the mitochondrial genome of the parrotfish (Teleostei: Perciformes: Scaridae)

Genomic size of animal mitochondrial DNA is usually minimized over time. Thus, when regional duplications occur, they are followed by a rapid elimination of redundant material. In contrast to this general view, we report here long-sustained tRNA pseudogenes in the mitochondrial genome (mitogenome) of teleost fishes of the family Scaridae (parrotfishes). During the course of a molecular phylogenetic study of the suborder Labroidei, we determined the complete nucleotide sequence of the mitogenome for a parrotfish, Chlorurus sordidus, and found a gene rearrangement accompanied by a tRNA pseudogene. In the typical gene order of vertebrates, a tRNA-gene cluster between ND1 and ND2 genes includes tRNA(Ile) (I), tRNA(Gln) (Q), and tRNA(Met) (M) genes in this order (IQM). However, in the mitogenome of the parrotfish, the tRNA(Met) gene was inserted between the tRNA(Ile) and the tRNA(Gln) genes, and the tRNA(Gln) gene was followed by a putative tRNA(Met) pseudogene (psiM). Such a tRNA gene rearrangement including a pseudogene (IMQpsiM) was found in all of the 10 examined species, representing 7 of the 10 currently recognized scarid genera. All sister groups examined (20 species of Labridae and a single species of Odacidae) had the typical gene order of vertebrate mitogenomes. Phylogenetic analysis of the tRNA(Met) genes and the resulting pseudogenes demonstrated that the ancestral tRNA(Met) gene was duplicated in a common ancestor of the parrotfish. Based on the fossil record, these results indicate that the pseudogenes have survived at least 14 million years. Most of the vertebrate mitochondrial gene rearrangements involving the IQM region have held the tRNA(Met) gene just upstream of the ND2 gene, and even in a few exceptional cases, including the present ones, the tRNA pseudogenes have been found in that position. In addition, most of these tRNA(Met) pseudogenes maintained clover-leaf secondary structures, with the remainder sustaining the clover-leaf structure in the "top half (TpsiC and acceptor arms). Considering their potential secondary structures (holding "top halves" of the clover-leaf structures), locations within mitogenomes (flanking the 5' ends of the ND2 genes) and stabilities over time (survived at least 14 Myr), it is likely that the tRNA pseudogenes retain function as punctuation marks for mitochondrial ND2 mRNA processing.

Mitogenomic evidence for the monophyly of elopomorph fishes (Teleostei) and the evolutionary origin of the leptocephalus larva

The monophyly of Elopomorpha (eels and their relatives) has long been one of the most problematic issues in systematic ichthyology. Since established the Elopomorpha based on the existence of the leaf-like larval form, termed a leptocephalus, no one has corroborated their monophyly using character matrices derived from both morphological and molecular data during the last 30 years. We investigated their monophyly and interrelationships at the ordinal level using complete mitochondrial genomic (mitogenomic) data from 33 purposefully chosen species (data for nine species being newly determined during the study) that fully represent the major teleostean and elopomorph lineages. Partitioned Bayesian analyses were conducted with the two data sets that comprised concatenated nucleotide sequences from 12 protein-coding genes (with and without third codon positions), 22 transfer RNA genes, and two ribosomal RNA genes. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high statistical values. Mitogenomic data strongly supported the monophyly of Elopomorpha, indicating the validity of the leptocephalus as an elopomorph synapomorphy. The order Elopiformes occupied the most basal position in the elopomorph phylogeny, with the Albuliformes and a clade comprising the Anguilliformes and the Saccopharyngiformes forming a sister group. The most parsimonious reconstruction of the three previously recognized, distinct larval types of elopomorphs onto the molecular phylogeny revealed that one of the types (fork-tailed type) had originated as the common ancestor of the Elopomorpha, the other two (filament-tailed and round-tailed types) having diversified separately in two more derived major clades.

Nonsynonymous Site Heteroplasmy in Fish Mitochondrial DNA

Heteroplasmic nucleotide polymorphisms are rarely observed in wild animal mitochondrial DNA. The occurrence of such site heteroplasmy is expected to be extremely rare at nonsynonymous sites where the number of nucleotide substitutions per site is low due to functional constraints. This report deals with nonsynonymous mitochondrial heteroplasmy from two wild fish species, chum salmon and Japanese flounder. We detected an A/C nonsynonymous heteroplasmic site corresponding to putative amino acids, Ile or Met, in NADH dehydrogenase subunit-5 (ND5) region of chum salmon. The heteroplasmic site was at the 3rd position of 58th codon. As for Japanese flounder we detected a C/T nonsynonymous heteroplasmic site corresponding to putative amino acids, Leu or Pro, in ND4 region. The heteroplasmic site was at the 2nd position of 450th codon. We also verified heteroplasmy at these sites by sequencing cloned fragments.

Efficiently resolving the basal clades of a phylogenetic tree using Bayesian and parsimony approaches: a case study using mitogenomic data from 100 higher teleost fishes

Many phylogenetic analyses that include numerous terminals but few genes show high resolution and branch support for relatively recently diverged clades, but lack of resolution and/or support for "basal" clades of the tree. The various benefits of increased taxon and character sampling have been widely discussed in the literature, albeit primarily based on simulations rather than empirical data. In this study, we used a well-sampled gene-tree analysis (based on 100 mitochondrial genomes of higher teleost fishes) to test empirically the efficiency of different methods of data sampling and phylogenetic inference to "correctly" resolve the basal clades of a tree (based on congruence with the reference tree constructed using all 100 taxa and 7990 characters). By itself, increased character sampling was an inefficient method by which to decrease the likelihood of "incorrect" resolution (i.e., incongruence with the reference tree) for parsimony analyses. Although increased taxon sampling was a powerful approach to alleviate "incorrect" resolution for parsimony analyses, it had the general effect of increasing the number of, and support for, "incorrectly" resolved clades in the Bayesian analyses. For both the parsimony and Bayesian analyses, increased taxon sampling, by itself, was insufficient to help resolve the basal clades, making this sampling strategy ineffective for that purpose. For this empirical study, the most efficient of the six approaches considered to resolve the basal clades when adding nucleotides to a dataset that consists of a single gene sampled for a small, but representative, number of taxa, is to increase character sampling and analyze the characters using the Bayesian method.

Molecular phylogeny of the antitropical genus Pseudolabrus (Perciformes: Labridae): evidence for a Southern Hemisphere origin

The genus Pseudolabrus comprises 11 species of marine nearshore fishes which are antitropically distributed: two species occur in East Asia, the remaining nine species being distributed in the Southern Hemisphere, mainly in the temperate Pacific. The distributions of their closely allied genera, collectively called "pseudolabrines" are, however, restricted to the Australia-New Zealand region. The molecular phylogeny of six of the 11 Pseudolabrus species from both Hemispheres and four of the five other pseudolabrine genera was reconstructed from nucleotide sequence data from mitochondrial DNA 12S rRNA, tRNAVal, and 16S rRNA genes. Both parsimony and Bayesian analyses were performed. Results are not consistent with a previous phylogenetic hypothesis based on osteological data, particularly in the relationship between Pseudolabrus and Notolabrus, indicating a probable need for reviewing the status of Notolabrus (or the delimitations of both Pseudolabrus and Notolabrus). The two Northern Hemisphere species of Pseudolabrus were monophyletic and nested deep into the clade of the Southern Hemisphere pseudolabrines, which indicates that both pseudolabrines and Pseudolabrus originated in the Southern Hemisphere. A dispersal rather than vicariance explanation for the antitropical distribution of Pseudolabrus is more parsimonious given the number of dispersal events, extinctions, and evolutionary adaptations required under the phylogeny. Based on molecular clock calibrations, the transequatorial divergence was suggested to be early to mid Pliocene at the earliest.

The complete mitochondrial genome of the javeline goby Acanthogobius hasta (Perciformes, Gobiidae) and phylogenetic considerations

We isolated Acanthogobius hasta mitochondrial DNA by long-polymerase chain reaction (long-PCR) with conserved primers, and sequenced this mitogenome with primer walking. The resultant A. hasta mitochondrial DNA sequence was found to consist of 16,663 bp with a structural organization conserved relative to that of other fish. In this paper, we report the basic characteristics of the A. hasta mitochondrial genome including structural organization, base composition of rRNAs and the tRNAs and protein-encoding genes, and characteristics of mitochondrial tRNAs. These findings are applicable to molecular phylogenetics in the suborder Gobioidei.

Molecular phylogeny and possible scenario of ponyfish (Perciformes:Leiognathidae) evolution

The family Leiognathidae, commonly known as ponyfish or slip mouth, comprises three genera, each being characterized mainly by mouth morphology. To date, however, neither the phylogenetic relationships within the family nor monophyly of the genera has been tested. The phylogenetic relationships among 14 species of Leiognathidae, inferred from two protein coding mitochondrial genes (ND4 and 5), indicated monophyly of the studied species form genera Gazza and Secutor, and paraphyly of the genus Leiognathus, with L. equulus occupying a basal branch of the family. The relationships allowed phylogenetic analyses of mouthpart structures and light organ systems. The results suggested that the morphology of the upwardly and forwardly protractile mouth types (latter with canine-like teeth) are phylogenetically informative, and the downwardly protractile mouth type being ancestral in the family. The results also suggested that internal sexual dimorphism of the light organ system was present in the common ancestor of a sister clade to L. equulus, whereas external sexual dimorphism seems to have evolved subsequently in two monophyletic subgroups.

Phylogenetic relationships of glyptosternoid fishes (Siluriformes: Sisoridae) inferred from mitochondrial cytochrome b gene sequences

To explore phylogenetic relationships among glyptosternoid fishes, we determined nucleotide sequences of the complete mitochondrial cytochrome b gene region (1138 base pair). Thirteen species of glyptosternoid fishes and six species of non-glyptosternoids represent 10 sisorid genera were examined. Molecular phylogenetic trees were constructed using the maximum parsimony, minimum evolution, maximum likelihood, and Bayesian methods. Bayesian and maximum likelihood analyses support the monophyly of glyptosternoids, but our hypothesis of internal relationships differs from previous hypothesis. Results indicated that glyptosternoid is a monophyletic group and genera Glyptosternum and Exostoma are two basal species having a primitive position among it. Genera Euchiloglanis and Pareuchiloglanis form a sister-group. Then they form a sister-group with Pseudexostoma plus Oreoglanis. Our result also found that Pareuchiloglanis anteanalis might be considered as the synonyms of Parechiloglanis sinensis, and genus Euchiloglanis might have only one valid species, Euchiloglanis davidi.

Complete sequence and characterization of the channel catfish mitochondrial genome

In order to support analysis of channel catfish populations and genetic improvement programs, the channel catfish, Ictalurus punctatus, mitochondrial genome was completely sequenced and revealed gene structure and gene order common to vertebrates. Nucleotide sequence comparisons of cytochrome b (Cytb) and cytochrome c oxidase subunit 1 (COI) demonstrated genetic separation of the genera Ictalurus, Pylodictis and Ameiurus consistent with the taxonomic classification within Ictaluridae. The ictalurid Cytb nucleotide sequences were significantly different from a putative channel catfish Cytb sequence in GenBank. Genetic relationships based on mitochondrial DNA sequences indicated the value of channel catfish in genomic comparisons between teleosts. Pairwise alignment of DNA sequences revealed conservation of regulatory sequences in the D-loop region with other vertebrates. Analysis of D-loop sequences in commercial populations and a research strain revealed 28 polymorphic sites and 33 D-loop haplotypes. Sequence analysis revealed clustering of haplotypes within commercial farms and the USDA103 research line, but D-loop haplotypes were not sufficient to discriminate the USDA103 fish from commercial catfish.

How Meaningful Are Bayesian Support Values?

In this study, we used an empirical example based on 100 mitochondrial genomes from higher teleost fishes to compare the accuracy of parsimony-based jackknife values with Bayesian support values. Phylogenetic analyses of 366 partitions, using differential taxon and character sampling from the entire data matrix of 100 taxa and 7,990 characters, were performed for both phylogenetic methods. The tree topology and branch-support values from each partition were compared with the tree inferred from all taxa and characters. Using this approach, we quantified the accuracy of the branch-support values assigned by the jackknife and Bayesian methods, with respect to each of 15 basal clades. In comparing the jackknife and Bayesian methods, we found that (1) both measures of support differ significantly from an ideal support index; (2) the jackknife underestimated support values; (3) the Bayesian method consistently overestimated support; (4) the magnitude by which Bayesian values overestimate support exceeds the magnitude by which the jackknife underestimates support; and (5) both methods performed poorly when taxon sampling was increased and character sampling was not increases. These results indicate that (1) the higher Bayesian support values are inappropriate (in magnitude), and (2) Bayesian support values should not be interpreted as probabilities that clades are correctly resolved. We advocate the continued use of the relatively conservative bootstrap and jackknife approaches to estimating branch support rather than the more extreme overestimates provided by the Markov Chain Monte Carlo-based Bayesian methods.

Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani

The complete nucleotide sequence (16,650 bp) of the mitochondrial genome of the salamander Mertensiella luschani (Caudata, Amphibia) was determined. This molecule conforms to the consensus vertebrate mitochondrial gene order. However, it is characterized by a long non-coding intervening sequence with two 124-bp repeats between the tRNA(Thr) and tRNA(Pro) genes. The new sequence data were used to reconstruct a phylogeny of jawed vertebrates. Phylogenetic analyses of all mitochondrial protein-coding genes at the amino acid level recovered a robust vertebrate tree in which lungfishes are the closest living relatives of tetrapods, salamanders and frogs are grouped together to the exclusion of caecilians (the Batrachia hypothesis) in a monophyletic amphibian clade, turtles show diapsid affinities and are placed as sister group of crocodiles+birds, and the marsupials are grouped together with monotremes and basal to placental mammals. The deduced phylogeny was used to characterize the molecular evolution of vertebrate mitochondrial proteins. Amino acid frequencies were analyzed across the main lineages of jawed vertebrates, and leucine and cysteine were found to be the most and least abundant amino acids in mitochondrial proteins, respectively. Patterns of amino acid replacements were conserved among vertebrates. Overall, cartilaginous fishes showed the least variation in amino acid frequencies and replacements. Constancy of rates of evolution among the main lineages of jawed vertebrates was rejected.

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

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

Phylogenetic and genomic analysis of the complete mitochondrial DNA sequence of the spotted asparagus beetle Crioceris duodecimpunctata

We report the complete mitochondrial DNA sequence of the spotted asparagus beetle, Crioceris duodecimpunctata. The genome complement, gene order, and nucleotide composition of this beetle's mitochondrial genome were found to be typical of those reported for other insects. Unusual features of this genome include the substitution of UCU for GCU as the anticodon for tRNA(Ser), an unusual TpsiC loop for the tRNA(Ile) gene, and the identification of a putative ATT start codon for cox1. The utility of complete mitochondrial genome data for phylogenetic inference of the insect orders was tested, and compared to that of cox1 and combined mitochondrial ribosomal DNA sequences. Even though the number of insect orders represented by complete mitochondrial genomes is still limited, several well-established relationships are evident in the phylogenetic analysis of the complete sequences. Monophyly of the orders Diptera, Lepidoptera, and Coleoptera were consistently recovered. Monophyly of the Holometabola was also observed in some (though not all) analyses. The accumulation of complete mitochondrial sequences from a broader array of insect orders holds the promise of clarifying the early diversification of insects.

Molecular phylogenetic evidence refuting the hypothesis of Batoidea (rays and skates) as derived sharks

Early morphological studies regarding the evolutionary history of elasmobranchs suggested sharks and batoids (skates and rays) were respectively monophyletic. More modern morphological cladistic studies, however, have tended to suggest that batoids are derived sharks, closely related to sawsharks and angelsharks, a phylogenetic arrangement known as the Hypnosqualea hypothesis. Very few molecular studies addressing interordinal relationships of elasmobranchs have been published; the few that do exist, are very limited in terms of both taxon representation and/or aligned sequence positions, and are insufficient to answer the question of whether batoids are derived sharks. The purpose of this study was to address this issue with more complete taxon representation, concomitant with a reasonable number of aligned sequence positions. The data set included a 2.4-kb segment of the mitochondrial 12S rRNA-tRNA valine-16S rRNA locus, and in terms of taxa, representatives of two orders of Batoidea, at least one representative of all orders of sharks, and as an outgroup, the widely recognized sister group to elasmobranchs-Holocephali. The results provide the first convincing molecular evidence for shark monophyly and the rejection of the Hypnosqualea hypothesis. Our phylogenetic placement of batoids as a basal elasmobranch lineage means that much of the current thinking regarding the evolution of morphological and life history characteristics in elasmobranchs needs to be re-evaluated.

Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences

A recent preliminary study using complete mitochondrial DNA sequences from 48 species of teleosts has suggested that higher teleostean phylogenies should be reinvestigated on the basis of more intensive taxonomic sampling. As a second step towards the resolution of higher teleostean phylogenies, which have been described as the "(unresolved) bush at the top of the tree," we reanalyzed their relationships using mitogenomic data from 100 purposefully chosen species that fully represented all of the higher teleostean orders, except for the Batrachoidiformes. Unweighted and weighted maximum parsimony analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding 3rd codon positions) and 21 transfer RNA (tRNA) genes (stem regions only) from each species. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high statistical values. All major, comprehensive groups above ordinal level as currently defined in higher teleosts (with the exception of the Neoteleostei and several monotypic groups), such as the Eurypterygii, Ctenosquamata, Acanthomorpha, Paracanthopterygii, Acanthopterygii, and Percomorpha, appeared to be nonmonophyletic in the present tree. Such incongruities largely resulted from differences in the placement and/or limits of the orders Ateleopodiformes, Lampridiformes, Polymixiiformes, Ophidiiformes, Lophiiformes, Beryciformes, Stephanoberyciformes, and Zeiformes, long-standing problematic taxa in systematic ichthyology. Of these, the resulting phylogenetic positions of the Ophidiiformes and Lophiiformes were totally unexpected, because, although they have consistently been considered relatively primitive groups within higher teleosts (Paracanthopterygii), they were confidently placed within a crown group of teleosts, herein called the Percomorpha. It should be noted that many unexpected, but highly supported relationships were found within the Percomorpha, being highly promising for the next investigative step towards resolution of this remarkably diversified group of teleosts.

Basal actinopterygian relationships: a mitogenomic perspective on the phylogeny of the "ancient fish"

The basal actinopterygians comprise four major lineages (polypteriforms, acipenseriforms, lepisosteids, and Amia) and have been collectively called "ancient fish." We investigated the phylogeny of this group of fishes in relation to teleosts using mitochondrial genomic (mitogenomic) data, and compared this to the various alternative phylogenetic hypotheses that have been proposed previously. In addition to the previously determined complete mitochondrial DNA (mtDNA) sequences from 14 teleosts and two outgroups, we used newly determined mitogenomic sequences of 12 purposefully chosen species representing all the ancient fish lineages plus related teleosts. This data set comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding the ND6 gene and third codon positions) and 22 transfer RNA (tRNA) genes (stem regions only) and these data were subjected to maximum parsimony, maximum likelihood, and Bayesian analyses. The resultant trees from the three methods were well resolved and largely congruent, with most internal branches being supported by high statistical values. Mitogenomic data strongly supported not only the monophyly of the teleosts (osteoglossomorphs and above), but also a sister-group relationship between the teleosts and a clade comprising the acipenseriforms, lepisosteids, and Amia, with the polypteriforms occupying the most basal position in the actinopterygian phylogeny. Although the tree topology differed from any of the previously proposed hypotheses based on morphology, it exhibited congruence with a recently proposed novel hypothesis based on nuclear markers.

The complete DNA sequence of the mitochondrial genome of the self-fertilizing fish Rivulus marmoratus (Cyprinodontiformes, Rivulidae) and the first description of duplication of a control region in fish

We isolated Rivulus marmoratus mitochondrial DNA by long-polymerase chain reaction with conserved primers, and sequenced it with 36 sets of internal conserved primers, which were designed from the extensive sequence similarities of mitochondrial DNA from several fish species. The R. marmoratus mitochondrial DNA has 17,329 bp with a conserved structural organization compared to those of other fish. Rivulus marmoratus mitochondrial DNA also has two nearly identical control regions. The basic characteristics of the R. marmoratus mitochondrial genome are discussed.

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

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

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

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