Wataru Iwasaki

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.

MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species

We developed a set of universal PCR primers (MiFish-U/E) for metabarcoding environmental DNA (eDNA) from fishes. Primers were designed using aligned whole mitochondrial genome (mitogenome) sequences from 880 species, supplemented by partial mitogenome sequences from 160 elasmobranchs (sharks and rays). The primers target a hypervariable region of the 12S rRNA gene (163–185 bp), which contains sufficient information to identify fishes to taxonomic family, genus and species except for some closely related congeners. To test versatility of the primers across a diverse range of fishes, we sampled eDNA from four tanks in the Okinawa Churaumi Aquarium with known species compositions, prepared dual-indexed libraries and performed paired-end sequencing of the region using high-throughput next-generation sequencing technologies. Out of the 180 marine fish species contained in the four tanks with reference sequences in a custom database, we detected 168 species (93.3%) distributed across 59 families and 123 genera. These fishes are not only taxonomically diverse, ranging from sharks and rays to higher teleosts, but are also greatly varied in their ecology, including both pelagic and benthic species living in shallow coastal to deep waters. We also sampled natural seawaters around coral reefs near the aquarium and detected 93 fish species using this approach. Of the 93 species, 64 were not detected in the four aquarium tanks, rendering the total number of species detected to 232 (from 70 families and 152 genera). The metabarcoding approach presented here is non-invasive, more efficient, more cost-effective and more sensitive than the traditional survey methods. It has the potential to serve as an alternative (or complementary) tool for biodiversity monitoring that revolutionizes natural resource management and ecological studies of fish communities on larger spatial and temporal scales.

Functional characterization of flavobacteria rhodopsins reveals a unique class of light-driven chloride pump in bacteria

Significance Ion-translocating, light-activated membrane proteins known as rhodopsins are found in all three domains of life. Proton-pumping rhodopsins, such as proteorhodopsin, are known to be broadly distributed in marine bacteria. The first known sodium-pumping rhodopsin was recently described in marine flavobacterium. We report the discovery and characterization of a unique type of light-activated ion-translocating rhodopsin that translocates chloride ions into the cell and is evolutionarily distinct from the other known rhodopsin chloride pump, halorhodopsin, found in haloarchaea. Our data show that rhodopsins with different ion specificities have evolved independently in marine bacteria, with individual strains containing as many as three functionally different rhodopsins. Light-activated, ion-pumping rhodopsins are broadly distributed among many different bacteria and archaea inhabiting the photic zone of aquatic environments. Bacterial proton- or sodium-translocating rhodopsins can convert light energy into a chemiosmotic force that can be converted into cellular biochemical energy, and thus represent a widespread alternative form of photoheterotrophy. Here we report that the genome of the marine flavobacterium Nonlabens marinus S1-08T encodes three different types of rhodopsins: Nonlabens marinus rhodopsin 1 (NM-R1), Nonlabens marinus rhodopsin 2 (NM-R2), and Nonlabens marinus rhodopsin 3 (NM-R3). Our functional analysis demonstrated that NM-R1 and NM-R2 are light-driven outward-translocating H+ and Na+ pumps, respectively. Functional analyses further revealed that the light-activated NM-R3 rhodopsin pumps Cl− ions into the cell, representing the first chloride-pumping rhodopsin uncovered in a marine bacterium. Phylogenetic analysis revealed that NM-R3 belongs to a distinct phylogenetic lineage quite distant from archaeal inward Cl−-pumping rhodopsins like halorhodopsin, suggesting that different types of chloride-pumping rhodopsins have evolved independently within marine bacterial lineages. Taken together, our data suggest that similar to haloarchaea, a considerable variety of rhodopsin types with different ion specificities have evolved in marine bacteria, with individual marine strains containing as many as three functionally different rhodopsins.

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.

CLOCK-Controlled Polyphonic Regulation of Circadian Rhythms through Canonical and Noncanonical E-Boxes

ABSTRACT In mammalian circadian clockwork, the CLOCK-BMAL1 complex binds to DNA enhancers of target genes and drives circadian oscillation of transcription. Here we identified 7,978 CLOCK-binding sites in mouse liver by chromatin immunoprecipitation-sequencing (ChIP-Seq), and a newly developed bioinformatics method, motif centrality analysis of ChIP-Seq (MOCCS), revealed a genome-wide distribution of previously unappreciated noncanonical E-boxes targeted by CLOCK. In vitro promoter assays showed that CACGNG, CACGTT, and CATG(T/C)G are functional CLOCK-binding motifs. Furthermore, we extensively revealed rhythmically expressed genes by poly(A)-tailed RNA-Seq and identified 1,629 CLOCK target genes within 11,926 genes expressed in the liver. Our analysis also revealed rhythmically expressed genes that have no apparent CLOCK-binding site, indicating the importance of indirect transcriptional and posttranscriptional regulations. Indirect transcriptional regulation is represented by rhythmic expression of CLOCK-regulated transcription factors, such as Krüppel-like factors (KLFs). Indirect posttranscriptional regulation involves rhythmic microRNAs that were identified by small-RNA-Seq. Collectively, CLOCK-dependent direct transactivation through multiple E-boxes and indirect regulations polyphonically orchestrate dynamic circadian outputs.

BioHackathon series in 2011 and 2012: penetration of ontology and linked data in life science domains

The application of semantic technologies to the integration of biological data and the interoperability of bioinformatics analysis and visualization tools has been the common theme of a series of annual BioHackathons hosted in Japan for the past five years. Here we provide a review of the activities and outcomes from the BioHackathons held in 2011 in Kyoto and 2012 in Toyama. In order to efficiently implement semantic technologies in the life sciences, participants formed various sub-groups and worked on the following topics: Resource Description Framework (RDF) models for specific domains, text mining of the literature, ontology development, essential metadata for biological databases, platforms to enable efficient Semantic Web technology development and interoperability, and the development of applications for Semantic Web data. In this review, we briefly introduce the themes covered by these sub-groups. The observations made, conclusions drawn, and software development projects that emerged from these activities are discussed.

CapR: revealing structural specificities of RNA-binding protein target recognition using CLIP-seq data

RNA-binding proteins (RBPs) bind to their target RNA molecules by recognizing specific RNA sequences and structural contexts. The development of CLIP-seq and related protocols has made it possible to exhaustively identify RNA fragments that bind to RBPs. However, no efficient bioinformatics method exists to reveal the structural specificities of RBP–RNA interactions using these data. We present CapR, an efficient algorithm that calculates the probability that each RNA base position is located within each secondary structural context. Using CapR, we demonstrate that several RBPs bind to their target RNA molecules under specific structural contexts. CapR is available at https://sites.google.com/site/fukunagatsu/software/capr.

Rapid Pathway Evolution Facilitated by Horizontal Gene Transfers across Prokaryotic Lineages

The evolutionary history of biological pathways is of general interest, especially in this post-genomic era, because it may provide clues for understanding how complex systems encoded on genomes have been organized. To explain how pathways can evolve de novo, some noteworthy models have been proposed. However, direct reconstruction of pathway evolutionary history both on a genomic scale and at the depth of the tree of life has suffered from artificial effects in estimating the gene content of ancestral species. Recently, we developed an algorithm that effectively reconstructs gene-content evolution without these artificial effects, and we applied it to this problem. The carefully reconstructed history, which was based on the metabolic pathways of 160 prokaryotic species, confirmed that pathways have grown beyond the random acquisition of individual genes. Pathway acquisition took place quickly, probably eliminating the difficulty in holding genes during the course of the pathway evolution. This rapid evolution was due to massive horizontal gene transfers as gene groups, some of which were possibly operon transfers, which would convey existing pathways but not be able to generate novel pathways. To this end, we analyzed how these pathways originally appeared and found that the original acquisition of pathways occurred more contemporaneously than expected across different phylogenetic clades. As a possible model to explain this observation, we propose that novel pathway evolution may be facilitated by bidirectional horizontal gene transfers in prokaryotic communities. Such a model would complement existing pathway evolution models.

Methylocaldum marinum sp. nov., a thermotolerant, methane-oxidizing bacterium isolated from marine sediments, and emended description of the genus Methylocaldum.

An aerobic, methane-oxidizing bacterium (strain S8(T)) was isolated from marine sediments in Kagoshima Bay, Japan. Phylogenetic analysis based on 16S rRNA gene sequences indicated that this strain is closely related to members of the genus Methylocaldum (97.6-97.9 % similarity) within the class Gammaproteobacteria. Strain S8(T) was a Gram-staining-negative, non-motile, coccoid or short rod-shaped organism. The temperature range for growth of strain S8(T) was 20-47 °C (optimum growth at 36 °C). It required NaCl (>0.5 %), tolerated up to 5 % NaCl and utilized methane and methanol. The major cellular fatty acid and major respiratory quinone were C16 : 0 and 18-methylene ubiquinone 8, respectively. The DNA G+C content was 59.7 mol%. Strain S8(T) possessed mmoX, which encodes soluble methane monooxygenase, as well as pmoA, which encodes the particulate methane monooxygenase. On the basis of this morphological, physiological, biochemical and genetic information, the first marine species in the genus Methylocaldum is proposed, with the name Methylocaldum marinum sp. nov. The type strain is S8(T) ( = NBRC 109686(T) = DSM 27392(T)). An emended description of the genus Methylocaldum is also provided.

From β- to α-Proteobacteria: The Origin and Evolution of Rhizobial Nodulation Genes nodIJ

Although many α- and some β-proteobacterial species are symbiotic with legumes, the evolutionary origin of nitrogen-fixing nodulation remains unclear. We examined α- and β-proteobacteria whose genomes were sequenced using large-scale phylogenetic profiling and revealed the evolutionary origin of two nodulation genes. These genes, nodI and nodJ (nodIJ), play key roles in the secretion of Nod factors, which are recognized by legumes during nodulation. We found that only the nodulating β-proteobacteria, including the novel strains isolated in this study, possess both nodIJ and their paralogous genes (DRA-ATPase/permease genes). Contrary to the widely accepted scenario of the α-proteobacterial origin of rhizobia, our exhaustive phylogenetic analysis showed that the entire nodIJ clade is included in the clade of Burkholderiaceae DRA-ATPase/permease genes, that is, the nodIJ genes originated from gene duplication in a lineage of the β-proteobacterial family. After duplication, the evolutionary rates of nodIJ were significantly accelerated relative to those of homologous genes, which is consistent with their novel function in nodulation. The likelihood analyses suggest that this accelerated evolution is not associated with changes in either nonsynonymous/synonymous substitution rates or transition/transversion rates, but rather, in the GC content. Although the low GC content of the nodulation genes has been assumed to reflect past horizontal transfer events from donor rhizobial genomes with low GC content, no rhizobial genome with such low GC content has yet been found. Our results encourage a reconsideration of the origin of nodulation and suggest new perspectives on the role of the GC content of bacterial genes in functional adaptation.