Kohji Mabuchi

MitoFish and MitoAnnotator: a mitochondrial genome database of fish with an accurate and automatic annotation pipeline.

Mitofish is a database of fish mitochondrial genomes (mitogenomes) that includes powerful and precise de novo annotations for mitogenome sequences. Fish occupy an important position in the evolution of vertebrates and the ecology of the hydrosphere, and mitogenomic sequence data have served as a rich source of information for resolving fish phylogenies and identifying new fish species. The importance of a mitogenomic database continues to grow at a rapid pace as massive amounts of mitogenomic data are generated with the advent of new sequencing technologies. A severe bottleneck seems likely to occur with regard to mitogenome annotation because of the overwhelming pace of data accumulation and the intrinsic difficulties in annotating sequences with degenerating transfer RNA structures, divergent start/stop codons of the coding elements, and the overlapping of adjacent elements. To ease this data backlog, we developed an annotation pipeline named MitoAnnotator. MitoAnnotator automatically annotates a fish mitogenome with a high degree of accuracy in approximately 5 min; thus, it is readily applicable to data sets of dozens of sequences. MitoFish also contains re-annotations of previously sequenced fish mitogenomes, enabling researchers to refer to them when they find annotations that are likely to be erroneous or while conducting comparative mitogenomic analyses. For users who need more information on the taxonomy, habitats, phenotypes, or life cycles of fish, MitoFish provides links to related databases. MitoFish and MitoAnnotator are freely available at http://mitofish.aori.u-tokyo.ac.jp/ (last accessed August 28, 2013); all of the data can be batch downloaded, and the annotation pipeline can be used via a web interface.

Mitogenomic evaluation of the historical biogeography of cichlids toward reliable dating of teleostean divergences

BackgroundRecent advances in DNA sequencing and computation offer the opportunity for reliable estimates of divergence times between organisms based on molecular data. Bayesian estimations of divergence times that do not assume the molecular clock use time constraints at multiple nodes, usually based on the fossil records, as major boundary conditions. However, the fossil records of bony fishes may not adequately provide effective time constraints at multiple nodes. We explored an alternative source of time constraints in teleostean phylogeny by evaluating a biogeographic hypothesis concerning freshwater fishes from the family Cichlidae (Perciformes: Labroidei).ResultsWe added new mitogenomic sequence data from six cichlid species and conducted phylogenetic analyses using a large mitogenomic data set. We found a reciprocal monophyly of African and Neotropical cichlids and their sister group relationship to some Malagasy taxa (Ptychochrominae sensu Sparks and Smith). All of these taxa clustered with a Malagasy + Indo/Sri Lankan clade (Etroplinae sensu Sparks and Smith). The results of the phylogenetic analyses and divergence time estimations between continental cichlid clades were much more congruent with Gondwanaland origin and Cretaceous vicariant divergences than with Cenozoic transmarine dispersal between major continents.ConclusionWe propose to add the biogeographic assumption of cichlid divergences by continental fragmentation as effective time constraints in dating teleostean divergence times. We conducted divergence time estimations among teleosts by incorporating these additional time constraints and achieved a considerable reduction in credibility intervals in the estimated divergence times.

Independent evolution of the specialized pharyngeal jaw apparatus in cichlid and labrid fishes

BackgroundFishes in the families Cichlidae and Labridae provide good probable examples of vertebrate adaptive radiations. Their spectacular trophic radiations have been widely assumed to be due to structural key innovation in pharyngeal jaw apparatus (PJA), but this idea has never been tested based on a reliable phylogeny. For the first step of evaluating the hypothesis, we investigated the phylogenetic positions of the components of the suborder Labroidei (including Pomacentridae and Embiotocidae in addition to Cichlidae and Labridae) within the Percomorpha, the most diversified (> 15,000 spp) crown clade of teleosts. We examined those based on 78 whole mitochondrial genome sequences (including 12 newly determined sequences) through partitioned Bayesian analyses with concatenated sequences (13,933 bp).ResultsThe resultant phylogenies indicated that the Labridae and the remaining three labroid families have diverged basally within the Percomorpha, and monophyly of the suborder was confidently rejected by statistical tests using Bayes factors.ConclusionThe resultant phylogenies indicated that the specified PJA evolved independently at least twice, once in Labridae and once in the common ancestor of the remaining three labroid families (including the Cichlidae). Because the independent evolution of pharyngeal jaws appears to have been followed by trophic radiations, we consider that our result supports, from the aspect of historical repeatability, the idea that the evolution of the specialized PJA provided these lineages with the morphological potential for their spectacular trophic radiations. The present result will provide a new framework for the study of functional morphology and genetic basis of their PJA.

Evolutionary history of Otophysi (Teleostei), a major clade of the modern freshwater fishes: Pangaean origin and Mesozoic radiation

BackgroundFreshwater 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 worlds 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.ResultsPartitioned 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.ConclusionsThe 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.

Interrelationships of Atherinomorpha (medakas, flyingfishes, killifishes, silversides, and their relatives): The first evidence based on whole mitogenome sequences

Series Atherinomorpha, with its plentiful number of species and highly diversified ecological and morphological characters, is the most successful fish group at the surface layer of the ocean and many freshwater habitats, comprising 1552 species classified into three orders, six suborders, 21 families, and 193 genera. The group includes one of the most important research model organisms, the medaka (Oryzias latipes), together with diverse fishes with morphological, physiological, and ecological specializations, such as highly developed pectoral fins to glide, self-fertilization, and live-bearing. In this study, we examined the whole mitochondrial genomes (mitogenomes) from 17 species representing all of the three orders and six suborders within Atherinomorpha, with data from 70 additional percomorph species as ingroups, and two non-percomorph outgroup species. We subjected the unambiguously aligned mitogenome sequences to partitioned maximum likelihood and Bayesian phylogenetic analyses. The resulting phylogenies recovered a monophyletic Atherinomorpha within the Percomorpha, and demonstrated its phylogenetic affinity to the percomorph fishes (including cichlids) spawning demersal eggs with filaments. This study, further, provided the first molecular evidence for the monophyly of the respective atherinomorph orders (Atheriniformes, Beloniformes, and Cyprinodontiformes) with high posterior probabilities and mostly high bootstrap values, providing an important basis for the future studies on the phylogeny and evolution of this diverse group.

Divergence time of the two regional medaka populations in Japan as a new time scale for comparative genomics of vertebrates

The southern and northern Japanese populations of the medaka fish provide useful tools to gain insights into the comparative genomics and speciation of vertebrates, because they can breed to produce healthy and fertile offspring despite their highly divergent genetic backgrounds compared with those of human–chimpanzee. Comparative genomics analysis has suggested that such large genetic differences between the two populations are caused by higher molecular evolutionary rates among the medakas than those of the hominids. The argument, however, was based on the assumption that the two Japanese populations diverged approximately at the same time (4.0–4.7 Myr ago) as the human–chimpanzee lineage (5.0–6.0 Myr ago). This can be misleading, because the divergence time of the two populations was calculated based on estimated, extremely higher molecular evolutionary rates of other fishes with an implicit assumption of a global molecular clock. Here we show that our estimate, based on a Bayesian relaxed molecular-clock analysis of whole mitogenome sequences from 72 ray-finned fishes (including 14 medakas), is about four times older than that of the previous study (18 Myr). This remarkably older estimate can be reconciled with the vicariant events of the Japanese archipelago, and the resulting rates of molecular evolution are almost identical between the medaka and hominid lineages. Our results further highlight the fact that reproductive isolation may not evolve despite a long period of geographical isolation.

Evolutionary origin of the Scombridae (tunas and mackerels): members of a paleogene adaptive radiation with 14 other pelagic fish families.

Uncertainties surrounding the evolutionary origin of the epipelagic fish family Scombridae (tunas and mackerels) are symptomatic of the difficulties in resolving suprafamilial relationships within Percomorpha, a hyperdiverse teleost radiation that contains approximately 17,000 species placed in 13 ill-defined orders and 269 families. Here we find that scombrids share a common ancestry with 14 families based on (i) bioinformatic analyses using partial mitochondrial and nuclear gene sequences from all percomorphs deposited in GenBank (10,733 sequences) and (ii) subsequent mitogenomic analysis based on 57 species from those targeted 15 families and 67 outgroup taxa. Morphological heterogeneity among these 15 families is so extraordinary that they have been placed in six different perciform suborders. However, members of the 15 families are either coastal or oceanic pelagic in their ecology with diverse modes of life, suggesting that they represent a previously undetected adaptive radiation in the pelagic realm. Time-calibrated phylogenies imply that scombrids originated from a deep-ocean ancestor and began to radiate after the end-Cretaceous when large predatory epipelagic fishes were selective victims of the Cretaceous-Paleogene mass extinction. We name this clade of open-ocean fishes containing Scombridae “Pelagia” in reference to the common habitat preference that links the 15 families.

Mitochondrial DNA analysis reveals cryptic large-scale invasion of non-native genotypes of common carp (Cyprinus carpio) in Japan.

Wild common carp (Cyprinus carpio) are probably suffering from biological invasions of conspecific domesticated strains. However, such invasions may be largely camouflaged by morphological similarities between introduced and native strains. We conducted a large survey of mitochondrial DNA sequences (complete D‐loop region) from 11 localities in Japan. From a total of 166 individuals, 28 haplotypes were determined to fit into six divergent clades. One of the six clades included 19 closely related haplotypes with moderate nucleotide differences; however, the remaining five clades each included either a single haplotype or two almost identical haplotypes. Phylogenetic analysis together with the previously published Eurasian haplotypes further demonstrated that the ‘monotypic’ clades were sisters to various Eurasian lineages, whereas the 19 related haplotypes formed a monophyletic group apart from the whole Eurasian clade. Given their monophyly and genetic diversity, the 19 related haplotypes were thought to originate from the Japanese native strain. Conversely, their phylogenetic affinities to Eurasian lineages and unnaturally low genetic diversities caused the haplotypes of the five monotypic clades to be considered as domesticated strains introduced from Eurasia. These hypotheses were supported by further evidences; i.e. the probable non‐native haplotypes were frequently found from Japanese domesticated strains, and the probable native population structure was rescued when the probable non‐native haplotypes were excluded from the analyses. This study revealed that almost half or more of the haplotypes in all of the locations studied originated from domesticated strains introduced from Eurasia.

Explosive Speciation of Takifugu: Another Use of Fugu as a Model System for Evolutionary Biology

Although the fugu Takifugu rubripes has attracted attention as a model organism for genomic studies because of its compact genome, it is not generally appreciated that there are approximately 25 closely related species with limited distributions in the waters of East Asia. We performed molecular phylogenetic analyses and constructed a time tree using whole mitochondrial genome sequences from 15 Takifugu species together with 10 outgroups to examine patterns of diversification. The resultant time tree showed that the modern Takifugu species underwent explosive speciation during the Pliocene 1.8-5.3 Ma, which is comparable with that of the Malawi cichlids and tropheine cichlids in Lake Tanganyika. Considering their limited distributions and remarkable variations in coloration, morphology, and behavior, the results of the present study strongly suggest that Takifugu species are strong candidates as a model system for evolutionary studies of speciation mechanisms in marine environments where few such organisms are available.

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 tRNAIle (I), tRNAGln (Q), and tRNAMet (M) genes in this order (IQM). However, in the mitogenome of the parrotfish, the tRNAMet gene was inserted between the tRNAIle and the tRNAGln genes, and the tRNAGln gene was followed by a putative tRNAMet pseudogene (ψM). Such a tRNA gene rearrangement including a pseudogene (IMQψM) 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 tRNAMet genes and the resulting pseudogenes demonstrated that the ancestral tRNAMet 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 tRNAMet 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 tRNAMet pseudogenes maintained clover-leaf secondary structures, with the remainder sustaining the clover-leaf structure in the „top half (TψC 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.