Evolutionary superscaffolding and chromosome anchoring to improve Anopheles genome assemblies.

BACKGROUND: New sequencing technologies have lowered financial barriers to whole genome sequencing, but resulting assemblies are often fragmented and far from 'finished'. Updating multi-scaffold drafts to chromosome-level status can be achieved through experimental mapping or re-sequencing efforts. Avoiding the costs associated with such approaches, comparative genomic analysis of gene order conservation (synteny) to predict scaffold neighbours (adjacencies) offers a potentially useful complementary method for improving draft assemblies. RESULTS: We evaluated and employed 3 gene synteny-based methods applied to 21 Anopheles mosquito assemblies to produce consensus sets of scaffold adjacencies. For subsets of the assemblies, we integrated these with additional supporting data to confirm and complement the synteny-based adjacencies: 6 with physical mapping data that anchor scaffolds to chromosome locations, 13 with paired-end RNA sequencing (RNAseq) data, and 3 with new assemblies based on re-scaffolding or long-read data. Our combined analyses produced 20 new superscaffolded assemblies with improved contiguities: 7 for which assignments of non-anchored scaffolds to chromosome arms span more than 75% of the assemblies, and a further 7 with chromosome anchoring including an 88% anchored Anopheles arabiensis assembly and, respectively, 73% and 84% anchored assemblies with comprehensively updated cytogenetic photomaps for Anopheles funestus and Anopheles stephensi. CONCLUSIONS: Experimental data from probe mapping, RNAseq, or long-read technologies, where available, all contribute to successful upgrading of draft assemblies. Our evaluations show that gene synteny-based computational methods represent a valuable alternative or complementary approach. Our improved Anopheles reference assemblies highlight the utility of applying comparative genomics approaches to improve community genomic resources.

Whole-genome sequence analysis of the Mycobacterium avium complex and proposal of the transfer of Mycobacterium yongonense to Mycobacterium intracellulare subsp. yongonense subsp. nov.

Bacterial whole-genome sequences contain informative features of their evolutionary pathways. Comparison of whole-genome sequences have become the method of choice for classification of prokaryotes, thus allowing the identification of bacteria from an evolutionary perspective, and providing data to resolve some current controversies. Currently, controversy exists about the assignment of members of the Mycobacterium avium complex, as is for the cases of Mycobacterium yongonense and 'Mycobacterium indicus pranii'. These two mycobacteria, closely related to Mycobacterium intracellulare on the basis of standard phenotypic and single gene-sequences comparisons, were not considered a member of such species on the basis on some particular differences displayed by a single strain. Whole-genome sequence comparison procedures, namely the average nucleotide identity and the genome distance, showed that those two mycobacteria should be considered members of the species M. intracellulare. The results were confirmed with other whole-genome comparison supplementary methods. According to the data provided, Mycobacterium yongonense and 'Mycobacterium indicus pranii' should be considered and renamed and included as members of M. intracellulare. This study highlights the problems caused when a novel species is accepted on the basis of a single strain, as was the case for M. yongonense. Based mainly on whole-genome sequence analysis, we conclude that M. yongonense should be reclassified as a subspecies of Mycobacterium intracellulareas Mycobacterium intracellularesubsp. yongonense and 'Mycobacterium indicus pranii' classified in the same subspecies as the type strain of Mycobacterium intracellulare and classified as Mycobacterium intracellularesubsp. intracellulare.

Organellar genomics: a useful tool to study evolutionary relationships and molecular evolution in Gracilariaceae (Rhodophyta).

Gracilariaceae has a worldwide distribution including numerous economically important species. We applied high-throughput sequencing to obtain organellar genomes (mitochondria and chloroplast) from 10 species of Gracilariaceae and, combined with published genomes, to infer phylogenies and compare genome architecture among species representing main lineages. We obtained similar topologies between chloroplast and mitochondrial genomes phylogenies. However, the chloroplast phylogeny was better resolved with full support. In this phylogeny, Melanthalia intermedia is sister to a monophyletic clade including Gracilaria and Gracilariopsis, which were both resolved as monophyletic genera. Mitochondrial and chloroplast genomes were highly conserved in gene synteny, and variation mainly occurred in regions where insertions of plasmid-derived sequences (PDS) were found. In mitochondrial genomes, PDS insertions were observed in two regions where the transcription direction changes: between the genes cob and trnL, and trnA and trnN. In chloroplast genomes, PDS insertions were in different positions, but generally found between psdD and rrs genes. Gracilariaceae is a good model system to study the impact of PDS in genome evolution due to the frequent presence of these insertions in organellar genomes. Furthermore, the bacterial leuC/leuD operon was found in chloroplast genomes of Gracilaria tenuistipitata, G. chilensis, and M. intermedia, and in extrachromosomal plasmid of G. vermiculophylla. Phylogenetic trees show two different origins of leuC/leuD: genes found in chloroplast and plasmid were placed with proteobacteria, and genes encoded in the nucleus were close to Viridiplantae and cyanobacteria.

GenFamClust: an accurate, synteny-aware and reliable homology inference algorithm.

BACKGROUND: Homology inference is pivotal to evolutionary biology and is primarily based on significant sequence similarity, which, in general, is a good indicator of homology. Algorithms have also been designed to utilize conservation in gene order as an indication of homologous regions. We have developed GenFamClust, a method based on quantification of both gene order conservation and sequence similarity. RESULTS: In this study, we validate GenFamClust by comparing it to well known homology inference algorithms on a synthetic dataset. We applied several popular clustering algorithms on homologs inferred by GenFamClust and other algorithms on a metazoan dataset and studied the outcomes. Accuracy, similarity, dependence, and other characteristics were investigated for gene families yielded by the clustering algorithms. GenFamClust was also applied to genes from a set of complete fungal genomes and gene families were inferred using clustering. The resulting gene families were compared with a manually curated gold standard of pillars from the Yeast Gene Order Browser. We found that the gene-order component of GenFamClust is simple, yet biologically realistic, and captures local synteny information for homologs. CONCLUSIONS: The study shows that GenFamClust is a more accurate, informed, and comprehensive pipeline to infer homologs and gene families than other commonly used homology and gene-family inference methods.

Comparative genomic of the teleost cathepsin B and H and involvement in bacterial induced immunity of miiuy croaker.

Cathepsins are a family of lysosomal proteases play different roles at physiological and pathological states and present in almost all animals as well as other organisms. Cathepsins B and H are both cysteine proteases of cathepsins. Cathepsin B and H have been studied playing parts in protein degradation/turnover, antigen presentation/processing and hormone maturation in mammals. However, little is known about the structures and functions of cathepsin B and H in teleosts. In the present study, we identified and characterized the full-length miiuy croaker (Miichthys miiuy) cathepsin B and H genes. The sequence analysis results showed that both cathepsin B and H contain the characteristics of papain family with a signal peptide, propeptide and mature peptide regions. The comparison of the genomic organizations and locations indicated the conserved synteny and mild evolution in the cathepsin B and H genes adjacent regions. In addition, the gene synteny analysis showed that miiuy croaker cathepsin B has a closer relationship to stickleback and fugu than to cave fish and zebrafish, and cathepsin H was most similar with the 2 subtype in tilapia and fugu. By phylogenetic analysis, miiuy croaker cathepsin B and H were all assigned to cysteine proteases, and with a close relationship to Salmo salar cathepsin B and Oplegnathus fasciatus cathepsin H, respectively. Quantitative real-time RT-PCR analysis results confirmed that cathepsin B and H genes expressed ubiquitously in all tested healthy tissues from miiuy croaker. Furthermore, up-regulated expression of the cathepsin B and H transcripts in liver, spleen and kidney after exposure upon Vibrio anguillarum suggested that they may play important roles in innate immune response and antigen processing of miiuy croaker.

Seawater acclimation and inositol monophosphatase isoform expression in the European eel (Anguilla anguilla) and Nile tilapia (Orechromis niloticus).

Inositol monophosphatase (IMPA) is responsible for the synthesis of inositol, a polyol that can function as an intracellular osmolyte helping re-establish cell volume when exposed to hypertonic environments. Some epithelial tissues in euryhaline teleosts such as the eel and tilapia encounter considerable hyperosmotic challenge when fish move from freshwater (FW) to seawater (SW) environments; however, the roles played by organic osmolytes, such as inositol, have yet to be determined. Syntenic analysis has indicated that, as a result of whole genome- and tandem-duplication events, up to six IMPA isoforms can exist within teleost genomes. Four isoforms are homologs of the mammalian IMPA1 gene, and two isoforms are homologs of the mammalian IMPA2 gene. Although the tissue-dependent isoform expression profiles of the teleost isoforms appear to be species-specific, it was primarily mRNA for the IMPA1.1 isoform that was upregulated in epithelial tissues after fish were transferred to SW (up to 16-fold in eel and 90-fold in tilapia). Although up-regulation of IMPA1.1 expression was evident in many tissues in the eel, more substantial increases in IMPA1.1 expression were found in tilapia tissues, where SW acclimation resulted in up to 2,000-fold increases in protein expression, 16-fold increases in enzyme activity and 15-fold increases in tissue inositol contents. Immunohistochemical studies indicated that the tissue and cellular distribution of IMPA1.1 protein differed slightly between eels and tilapia; however, in both species the basal epithelial cell layers within the skin and fin, and the branchial epithelium and interstitial cells within the kidney, exhibited high levels of IMPA1.1 protein expression.

Identification and functional characterization of interleukin-12 receptor beta 1 and 2 in grass carp (Ctenopharyngodon idella).

Interleukin 12 (IL-12) binds its receptor complex of IL-12 receptor beta 1 (IL-12Rß1) and IL-12Rß2 to transduce cellular signaling in mammals. In teleosts, the function of Il-12 is drawing increasing attention, but molecular and functional features of Il-12 receptors remain obscure. Especially, the existence of multiple Il-12 isoforms in some fish species elicits the requirement to clarify their receptors. In this study, we isolated three cDNA sequences as Il-12 receptor candidates from grass carp, entitled as grass carp Il-12rß1 (gcIl-12rß1), gcIl-12rß2a and gcIl-12rß2b. In silico analysis showed that gcIl-12rß1 and gcIl-12rß2a shared the conserved gene locus and similar structure characteristics with their orthologues of zebrafish, frog, chicken, mouse and human, respectively. However, the Il-12rß2b of grass carp and zebrafish was similar to IL-27Ra in non-fish species. Further locally installed BLAST and gene synteny analysis uncovered three gcIl-12 receptors being single copied genes. Tissue distribution assay revealed that gcil12rß1 and gcil12rß2a transcripts were predominantly expressed in head kidney, differing from the even distribution of gcil12rß2b transcripts in all detected tissues. Subsequently, the binding ability and antagonistic effects of recombinant extracellular region of gcIl-12rß1 with recombinant grass carp Il-12 (rgcIl-12) isoforms were explored, providing functional evidence of the newly cloned gcIl-12rß1 being genuine orthologues of mammalian IL-12Rß1. Moreover, our data showed that gcIl-12rß1 and gcIl-12rß2a but not gcIl-12rß1 and gcIl-12rß2b mediated the effects of rgcIl-12 isoforms on ifn-γ promoter activity, thereby revealing Il-12 receptor signaling in fish. These results identified grass carp Il-12 receptors, thereby advancing our understanding of Il-12 isoform signaling in fish.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search.

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

Evolution of the genes mediating phototransduction in rod and cone photoreceptors.

This paper reviews current knowledge of the evolution of the multiple genes encoding proteins that mediate the process of phototransduction in rod and cone photoreceptors of vertebrates. The approach primarily involves molecular phylogenetic analysis of phototransduction protein sequences, combined with analysis of the syntenic arrangement of the genes. At least 35 of these phototransduction genes appear to reside on no more than five paralogons - paralogous regions that each arose from a common ancestral region. Furthermore, it appears that such paralogs arose through quadruplication during the two rounds of genome duplication (2R WGD) that occurred in a chordate ancestor prior to the vertebrate radiation, probably around 600 millions years ago. For several components of the phototransduction cascade, it is shown that distinct isoforms already existed prior to WGD, with the likely implication that separate classes of scotopic and photopic photoreceptor cells had already evolved by that stage. The subsequent quadruplication of the entire genome then permitted the refinement of multiple distinct protein isoforms in rods and cones. A unified picture of the likely pattern and approximate timing of all the important gene duplications is synthesised, and the implications for our understanding of the evolution of rod and cone phototransduction are presented.

Contrasting Evolutionary Paths Among Indo-Pacific Pomacentrus Species Promoted by Extensive Pericentric Inversions and Genome Organization of Repetitive Sequences.

Pomacentrus (damselfishes) is one of the most characteristic groups of fishes in the Indo-Pacific coral reef. Its 77 described species exhibit a complex taxonomy with cryptic lineages across their extensive distribution. Periods of evolutionary divergences between them are very variable, and the cytogenetic events that followed their evolutionary diversification are largely unknown. In this respect, analyses of chromosomal divergence, within a phylogenetic perspective, are particularly informative regarding karyoevolutionary trends. As such, we conducted conventional cytogenetic and cytogenomic analyses in four Pomacentrus species (Pomacentrus similis, Pomacentrus auriventris, Pomacentrus moluccensis, and Pomacentrus cuneatus), through the mapping of repetitive DNA classes and transposable elements, including 18S rDNA, 5S rDNA, (CA)15, (GA)15, (CAA)10, Rex6, and U2 snDNA as markers. P. auriventris and P. similis, belonging to the Pomacentrus coelestis complex, have indistinguishable karyotypes (2n = 48; NF = 48), with a peculiar syntenic organization of ribosomal genes. On the other hand, P. moluccensis and P. cuneatus, belonging to another clade, exhibit very different karyotypes (2n = 48, NF = 86 and 92, respectively), with a large number of bi-armed chromosomes, where multiple pericentric inversions played a significant role in their karyotype organization. In this sense, different chromosomal pathways followed the phyletic diversification in the Pomacentrus genus, making possible the characterization of two well-contrasting species groups regarding their karyotype features. Despite this, pericentric inversions act as an effective postzygotic barrier in many organisms, which appear to be also the case for P. moluccensis and P. cuneatus; the extensive chromosomal similarities in the two species of P. coelestis complex suggest minor participation of chromosomal postzygotic barriers in the phyletic diversification of these species.

misMM: An Integrated Pipeline for Misassembly Detection Using Genotyping-by-Sequencing and Its Validation with BAC End Library Sequences and Gene Synteny.

As next-generation sequencing technologies have advanced, enormous amounts of whole-genome sequence information in various species have been released. However, it is still difficult to assemble the whole genome precisely, due to inherent limitations of short-read sequencing technologies. In particular, the complexities of plants are incomparable to those of microorganisms or animals because of whole-genome duplications, repeat insertions, and Numt insertions, etc. In this study, we describe a new method for detecting misassembly sequence regions of Brassica rapa with genotyping-by-sequencing, followed by MadMapper clustering. The misassembly candidate regions were cross-checked with BAC clone paired-ends library sequences that have been mapped to the reference genome. The results were further verified with gene synteny relations between Brassica rapa and Arabidopsis thaliana. We conclude that this method will help detect misassembly regions and be applicable to incompletely assembled reference genomes from a variety of species.

New Insights into the Diversity of the Genus Faecalibacterium.

Faecalibacterium prausnitzii is a commensal bacterium, ubiquitous in the gastrointestinal tracts of animals and humans. This species is a functionally important member of the microbiota and studies suggest it has an impact on the physiology and health of the host. F. prausnitzii is the only identified species in the genus Faecalibacterium, but a recent study clustered strains of this species in two different phylogroups. Here, we propose the existence of distinct species in this genus through the use of comparative genomics. Briefly, we performed analyses of 16S rRNA gene phylogeny, phylogenomics, whole genome Multi-Locus Sequence Typing (wgMLST), Average Nucleotide Identity (ANI), gene synteny, and pangenome to better elucidate the phylogenetic relationships among strains of Faecalibacterium. For this, we used 12 newly sequenced, assembled, and curated genomes of F. prausnitzii, which were isolated from feces of healthy volunteers from France and Australia, and combined these with published data from 5 strains downloaded from public databases. The phylogenetic analysis of the 16S rRNA sequences, together with the wgMLST profiles and a phylogenomic tree based on comparisons of genome similarity, all supported the clustering of Faecalibacterium strains in different genospecies. Additionally, the global analysis of gene synteny among all strains showed a highly fragmented profile, whereas the intra-cluster analyses revealed larger and more conserved collinear blocks. Finally, ANI analysis substantiated the presence of three distinct clusters-A, B, and C-composed of five, four, and four strains, respectively. The pangenome analysis of each cluster corroborated the classification of these clusters into three distinct species, each containing less variability than that found within the global pangenome of all strains. Here, we propose that comparison of pangenome subsets and their associated α values may be used as an alternative approach, together with ANI, in the in silico classification of new species. Altogether, our results provide evidence not only for the reconsideration of the phylogenetic and genomic relatedness among strains currently assigned to F. prausnitzii, but also the need for lineage (strain-based) differentiation of this taxon to better define how specific members might be associated with positive or negative host interactions.

The evolution of bacterial mechanosensitive channels.

Mechanosensitive channels are ubiquitous and highly studied. However, the evolution of the bacterial channels remains enigmatic. It can be argued that mechanosensitivity might be a feature of all membrane proteins with some becoming progressively less sensitive to membrane tension over the course of evolution. Bacteria and archaea exhibit two main classes of channels, MscS and MscL. Present day channels suggest that the evolution of MscL may be highly constrained, whereas MscS has undergone elaboration via gene fusion (and potentially gene fission) events to generate a diversity of channel structures. Some of these channel variants are constrained to a small number of genera or species. Some are only found in higher organisms. Only exceptionally have these diverse channels been investigated in any detail. In this review we consider both the processes that might have led to the evolved complexity but also some of the methods exploiting the explosion of genome sequences to understand (and/or track) their distribution. The role of MscS-related channels in calcium-mediated cell biology events is considered.

Molecular characterization of three IRF1 subfamily members reveals evolutionary significance of IRF11 in miiuy croaker.

The interferon regulatory factors IRF1 and IRF2 of the IRF1 subfamily play essential roles in immune responses against viruses. IRF11 is a novel IRF gene of the IRF1 subfamily; IRF11 genes share almost the same evolutionary distance with IRF1 and IRF2 genes. However, the structure and characteristics of IRF11 gene in fish have been rarely reported. In our study, IRF1, IRF2 and IRF11 genes were identified and characterized from miiuy croaker genome. Results showed that the IRF1, IRF2 and IRF11 genes contain the same domains; each of these genes is composed of conserved gene organizations and characterized by gene synteny with the orthologous genes. Interestingly, IRF11 was likely found only in fish (but not specific to teleost fish). Evolutionary analysis results showed that IRF1 gene in mammals, IRF2 and IRF11 gene in fish underwent positive selection. IRF1, IRF2 and IRF11 were expressed in a wide range of miiuy croaker tissues. These genes also exhibited the same expression patterns after miiuy croaker was infected with poly(I:C). Therefore, our data enhanced our understanding of the functions and evolution of IRF11 in fish.

Characterization of the duplicate L-SIGN and DC-SIGN genes in miiuy croaker and evolutionary analysis of L-SIGN in fishes.

Dendritic cell-specific ICAM-3-grabbing non-integrin (DC-SIGN/CD209) and liver/lymph node-specific ICAM-grabbing non-integrin (L-SIGN/CD299) which are homologues of DC-SIGN are important members in C-type lectin receptors family as key molecules to recognize and eliminate pathogens in the innate immune system. DC-SIGN and L-SIGN have become hot topics in recent studies which both served as cell adhesion and phagocytic pathogen recognition receptors in mammals. However, there have been almost no studies of DC-SIGN and L-SIGN structure and characters in fish, only DC-SIGN in the zebrafish had been studied. In our study, we identified and characterized the full-length miiuy croaker (Miichthys miiuy) DC-SIGN (mmDC-SIGN) and L-SIGN (mmL-SIGN) genes. The sequence analysis results showed that mmDC-SIGN and mmL-SIGN have the same domains with other vertebrates except primates, and share some conserved motifs in CRD among all the vertebrates which play a crucial role in interacting with Ca(2+) and for recognizing mannose-containing motifs. Gene synteny of DC-SIGN and L-SIGN were analyzed for the first time and gene synteny of L-SIGN was conserved among the five fishes. Interestingly, one gene next to L-SIGN from gene synteny had high similarity with L-SIGN gene that was described as L-SIGN-like in fish species. While only one L-SIGN gene existed in other vertebrates, two L-SIGN in fish may be in consequence of the fish-specific genome duplication to adapt the specific environment. The evolutionary analysis showed that the ancestral lineages of L-SIGN gene in fishes experienced purifying selection and the current lineages of L-SIGN gene in fishes underwent positive selection, indicating that the ancestral lineages and current lineages of L-SIGN gene in fishes underwent different evolutionary patterns. Both mmDC-SIGN and mmL-SIGN were expressed in all tested tissues and ubiquitously up-regulated in infected liver, spleen and kidney at different sampling time points, indicating that the mmDC-SIGN and mmL-SIGN participated in the immune response to defense against bacteria infection.

Comparative and evolutionary analysis of major peanut allergen gene families.

Peanut (Arachis hypogaea L.) causes one of the most serious food allergies. Peanut seed proteins, Arah1, Arah2, and Arah3, are considered to be among the most important peanut allergens. To gain insights into genome organization and evolution of allergen-encoding genes, approximately 617 kb from the genome of cultivated peanut and 215 kb from a wild relative were sequenced including three Arah1, one Arah2, eight Arah3, and two Arah6 gene family members. To assign polarity to differences between homoeologous regions in peanut, we used as outgroups the single orthologous regions in Medicago, Lotus, common bean, chickpea, and pigeonpea, which diverged from peanut about 50 Ma and have not undergone subsequent polyploidy. These regions were also compared with orthologs in many additional dicot plant species to help clarify the timing of evolutionary events. The lack of conservation of allergenic epitopes between species, and the fact that many different proteins can be allergenic, makes the identification of allergens across species by comparative studies difficult. The peanut allergen genes are interspersed with low-copy genes and transposable elements. Phylogenetic analyses revealed lineage-specific expansion and loss of low-copy genes between species and homoeologs. Arah1 syntenic regions are conserved in soybean, pigeonpea, tomato, grape, Lotus, and Arabidopsis, whereas Arah3 syntenic regions show genome rearrangements. We infer that tandem and segmental duplications led to the establishment of the Arah3 gene family. Our analysis indicates differences in conserved motifs in allergen proteins and in the promoter regions of the allergen-encoding genes. Phylogenetic analysis and genomic organization studies provide new insights into the evolution of the major peanut allergen-encoding genes.

Genomic organization and transcription of the medaka and zebrafish cellular retinol-binding protein (rbp) genes.

In this study, we examined the evolutionary trajectories and the common ancestor of medaka rbp genes by comparing them to the well-studied rbp/RBP genes from zebrafish and other vertebrates. We describe here gene structure, sequence identity, phylogenetic analysis and conserved gene synteny of medaka rbp genes and their putative proteins as well as the tissue-specific distribution of rbp transcripts in adult medaka and zebrafish. Medaka rbp genes consist of four exons separated by three introns that encode putative polypeptides of 134-138 amino acids, a genomic organization characteristic of rbp genes. Medaka Rbp sequences share highest sequence identity and similarity with their orthologs in vertebrates, and form a distinct clade with them in phylogenetic analysis. Conserved gene synteny was evident among medaka, zebrafish and human rbp/RBP genes, which provides compelling evidence that the medaka rbp1, rbp2a, rbp2b, rbp5, rbp7a and rbp7b genes arose from a common ancestor of vertebrates. Moreover, the duplicated rbp2 and rbp7 genes most likely exist owing to a whole-genome duplication (WGD) event specific to the teleost fish lineage. Selection pressure and the nonparametric relative rate test of the medaka and zebrafish duplicated rbp2 and rbp7 genes suggest that these duplicated genes are subjected to purifying selection and one paralog might have evolved at an accelerated rate compared to its sister duplicate since the WGD. The steady-state levels of medaka and zebrafish rbp1, rbp2a, rbp2b and rbp5 transcripts in various tissues suggest that medaka rbp1, rbp2a and rbp2b genes have retained the regulatory elements of an ancestral RBP1 and RBP2 genes, and the medaka rbp5 gene has acquired new function. Furthermore, the tissue-specific regulations of rbp7a and rbp7b genes have diverged markedly in medaka and zebrafish since the teleost-specific WGD.