A first generation physical map of the medaka genome in BACs essential for positional cloning and clone-by-clone based genomic sequencing

In order to realize the full potential of the medaka as a model system for developmental biology and genetics, characterized genomic resources need to be established, culminating in the sequence of the medaka genome. To facilitate the map-based cloning of genes underlying induced mutations and to provide templates for clone-based genomic sequencing, we have created a first-generation physical map of the medaka genome in bacterial artificial chromosome (BAC) clones. In particular, we exploited the synteny to the closely related genome of the pufferfish, Takifugu rubripes, by marker content mapping. As a first step, we clustered 103,144 public medaka EST sequences to obtain a set of 21,121 non-redundant sequence entities. Avoiding oversampling of gene-dense regions, 11,254 of EST clusters were successfully matched against the draft sequence of the fugu genome, and 2363 genes were selected for the BAC map project. We designed 35mer oligonucleotide probes from the selected genes and hybridized them against 64,500 BAC clones of strains Cab and Hd-rR, representing 14-fold coverage of the medaka genome. Our data set is further supplemented with 437 results generated from PCR-amplified inserts of medaka cDNA clones and BAC end-fragment markers. Our current, edited, first generation medaka BAC map consists of 902 map segments that cover about 74% of the medaka genome. The map contains 2721 markers. Of these, 2534 are from expressed sequences, equivalent to a non-redundant set of 2328 loci. The 934 markers (724 different) are anchored to the medaka genetic map. Thus, genetic map assignments provide immediate access to underlying clones and contigs, simplifying molecular access to candidate gene regions and their characterization.     

An in silico mining for simple sequence repeats from expressed sequence tags of zebrafish, medaka, Fundulus, and Xiphophorus

Teleost fish genome projects involving model species are resulting in a rapid accumulation of genomic and expressed DNA sequences in public databases. The expressed sequence tags (ESTs) collected in the databases can be mined for the analysis of both structural and functional genomics. In this study, we in silico analyzed 49,430 unigenes representing a total of 692,654 ESTs from four model fish for their potential use in developing simple sequence repeats (SSRs), or microsatellites. After bioinformatical mining, a total of 3,018 EST derived SSRs (EST-SSRs) were identified for 2,335 SSR containing ESTs (SSR-ESTs). The frequency of identified SSR-ESTs ranged from 1.5% for Xiphophorus to 7.3% for zebrafish. The dinucleotide repeat motif is the most abundant SSR, accounting for 47%, 52%, 64%, and 78% for medaka, Fundulus, zebrafish, and Xiphophorus, respectively. Simulation analysis suggests that a majority of these EST-SSRs have sufficient flanking sequences for polymerase chain reaction (PCR) primer design. Comparative DNA sequence analyses of SSR-ESTs identified several cross-species SSRs and sequences that may be used as cross-reference genes in comparative studies. For example, the flanking sequences of one SSR (CTG)n within the pituitary tumor-transforming gene (PTTG) 1 interacting protein (PTTGIP), showed conservation spanning the medaka, Fundulus, human, and mouse genomes. This study provides a large body of information on EST-SSRs that can be useful for the development of polymorphic markers, gene mapping, and comparative genome analysis. Functional analysis of these SSR-ESTs may reveal their role in metabolism and gene evolution of these model species.     

Immunoglobulin light chains in medaka (Oryzias latipes)

The gene segments encoding antibodies have been studied in many capacities and represent some of the best-characterized gene families in traditional animal disease models (mice and humans). To date, multiple immunoglobulin light chain (IgL) isotypes have been found in vertebrates and it is unclear as to which isotypes might be more primordial in nature. Sequence data emerging from an array of fish genome projects is a valuable resource for discerning complex multigene assemblages in this critical branch point of vertebrate phylogeny. Herein, we have analyzed the genomic organization of medaka (Oryzias latipes) IgL gene segments based on recently released genome data. The medaka IgL locus located on chromosome 11 contains at least three clusters of IgL gene segments comprised of multiple gene assemblages of the kappa light chain isotype. These data suggest that medaka IgL gene segments may undergo both intra- and inter-cluster rearrangements as a means to generate additional diversity. Alignments of expressed sequence tags to concordant gene segments which revealed each of the three IgL clusters are expressed. Collectively, these data provide a genomic framework for IgL genes in medaka and indicate that Ig diversity in this species is achieved from at least three distinct chromosomal regions.     

The genome size evolution of medaka (Oryzias latipes) and fugu (Takifugu rubripes)

Evolution of the genome size in eukaryotes is often affected by changes in the noncoding sequences, for which insertions and deletions (indels) of small nucleotide sequences and amplification of repetitive elements are considered responsible. In this study, we compared the genomic DNA sequences of two kinds of fish, medaka (Oryzias latipes) and fugu (Takifugu rubripes), which show two-fold difference in the genome size (800 Mb vs. 400 Mb). We selected a contiguous DNA sequence of 790 kb from the medaka chromosome LG22 (linkage group 22), and made a precise comparison with the sequence (387 kb) of the corresponding region of Takifugu. The sequence of 178 kb in total was aligned common between two fishes, and the remaining sequences (612 kb for medaka and 209 kb for fugu) were found abundant in various repetitive elements including many types of unclassified low copy repeats, all of which accounted for more than a half (54%) of the genome size difference. Furthermore, we identified a significant difference in the length ratio of the unaligned sequences that locate between the aligned sequences (USBAS), particularly after eliminating known repetitive elements. These USBAS with no repetitive elements (USBAS-nr) located within the intron and intergenic region. These results strongly indicated that amplification of repetitive elements and compilation of indels are major driving forces to facilitate changes in the genome size.     

The genome of the marine medaka Oryzias melastigma

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field‐based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole‐genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.     

3640 Unique EST Clusters from the Medaka Testis and Their Potential Use for Identifying Conserved Testicular Gene Expression in Fish and Mammals

Background

The fish medaka is the first vertebrate capable of full spermatogenesis in vitro from self-renewing spermatogonial stem cells to motile test-tube sperm. Precise staging and molecular dissection of this process has been hampered by the lack of suitable molecular markers.

Methodology and Principal Findings

We have generated a normalized medaka testis cDNA library and obtained 7040 high quality sequences representing 3641 unique gene clusters. Among these, 1197 unique clusters are homologous to known genes, and 2444 appear to be novel genes. Ontology analysis shows that the 1197 gene products are implicated in diverse molecular and cellular processes. These genes include markers for all major types of testicular somatic and germ cells. Furthermore, markers were identified for major spermatogenic stages ranging from spermatogonial stem cell self-renewal to meiosis entry, progression and completion. Intriguingly, the medaka testis expresses at least 13 homologs of the 33 mouse X-chromosomal genes that are enriched in the testis. More importantly, we show that key components of several signaling pathways known to be important for testicular function in mammals are well represented in the medaka testicular EST collection.

Conclusions/Significance

Medaka exhibits a considerable similarity in testicular gene expression to mammals. The medaka testicular EST collection we obtained has wide range coverage and will not only consolidate our knowledge on the comparative analysis of known genes'' functions in the testis but also provide a rich resource to dissect molecular events and mechanism of spermatogenesis in vivo and in vitro in medaka as an excellent vertebrate model.     

Comparative genomics of medaka and fugu

A small freshwater fish medaka (Oryzias latipes) has been one of the most attractive experimental systems for research in genetics and developmental biology. We have formed an international consortium Medaka Genome Initiative (MGI) to collect and share various information and resources on medaka. The MGI has set an ambitious goal aiming at the complete sequencing of the medaka genome and as a feasibility study we have begun sequencing one particular chromosome, linkage group 22 (LG22) of 22 Mb in size. Initial sequence analysis revealed unique features of the medaka genome in comparison to fugu genome.     

Hox gene complexity in medaka fish may Be similar to that in pufferfish rather than zebrafish.

Changes in number and the genomic organization of Hox genes have played an important role in metazoan body-plan evolution. They make cluster(s), and in vertebrates, each cluster contains different number of Hox genes that have been classified into 13 groups. There are 39 Hox genes in four clusters on different chromosomes in the mammalian genome. In the fish, while 31 Hox genes in four clusters have been identified in pufferfish Fugu rubripes, 47 Hox genes in seven clusters exist in the zebrafish Danio rerio. To estimate the evolutionary origin of Hox organization in ray-finned fishes, we searched for Hox genes in the medaka fish Oryzias latipes, with a taxon thought to be widely separated from those of pufferfish and zebrafish. We synthesized various mixed oligonucleotides that can work as group-specific primers for PCR, then cloned and sequenced amplified fragments. Numbers of Hox genes identified in the present study were 2 for group 1, 2 for group 2, 1 for group 3, 3 for group 4, 6 for groups 5-7, 2 for group 8, 4 for group 9, 3 for group 10, 1 for group 12, and 3 for group 13. The primers specific for group 11 did not function in this study. Thus, at least 27 Hox genes are present in medaka genome, suggesting that the Hox gene complexity of the medaka genome is similar to that of the pufferfish rather than the zebrafish.     

Maternal-zygotic medaka mutants for fgfr1 reveal its essential role in the migration of the axial mesoderm but not the lateral mesoderm

The medaka fish (Oryzias latipes) is an emerging model organism for which a variety of unique developmental mutants have now been generated. Our recent mutagenesis screening of the medaka identified headfish (hdf), a null mutant for fgf receptor 1 (fgfr1), which fails to develop structures in the trunk and tail. Despite its crucial role in early development, the functions of Fgfr1-mediated signaling have not yet been well characterized due to the complexity of the underlying ligand-receptor interactions. In our present study, we further elucidate the roles of this pathway in the medaka using the hdf (fgfr1) mutant. Because Fgfr1 is maternally supplied in fish, we first generated maternal-zygotic (MZ) mutants by transplanting homozygous hdf germ cells into sterile interspecific hybrids. Interestingly, the host hybrid fish recovered their fertility and produced donor-derived mutant progeny. The resulting MZ mutants also exhibited severe defects in their anterior head structures that are never observed in the corresponding zygotic mutants. A series of detailed analyses subsequently revealed that Fgfr1 is required for the anterior migration of the axial mesoderm, particularly the prechordal plate, in a cell-autonomous manner, but is not required for convergence movement of the lateral mesoderm. Furthermore, fgfr1 was found to be dispensable for initial mesoderm induction. The MZ hdf medaka mutant was thus found to be a valuable model system to analyze the precise role of fgfr1-mediated signaling in vertebrate early development.     

Duplicated Abd-B class genes in medaka hoxAa and hoxAb clusters exhibit differential expression patterns in pectoral fin buds

Hox genes form clusters. Invertebrates and Amphioxus have only one hox cluster, but in vertebrates, they are multiple, i.e., four in the basal teleost fish Polyodon and tetrapods (HoxA, B, C, D), but seven or eight in common teleosts. We earlier completely sequenced the entire hox gene loci in medaka fish, showing a total of 46 hox genes to be encoded in seven clusters (hoxAa, Ab, Ba, Bb, Ca, Da, Db). Among them, hoxAa, hoxAb and hoxDa clusters are presumed to be important for fin-to-limb evolution because of their key role in forelimb and pectoral fin development. In the present study, we compared genome organization and nucleotide sequences of the hoxAa and hoxAb clusters to these of tetrapod HoxA clusters, and found greater similarity in hoxAa case. We then analyzed expression of Abd-B family genes in the clusters. In the trunk, those from the hoxAa cluster, i.e., hoxA9a, hoxA10a, hoxA11a and hoxA13a, were expressed in a manner keeping the colinearity rule of the hox expression as those of tetrapods, while those from the hoxAb cluster, i.e., hoxA9b, hoxA10b, hoxA11b and hoxA13b, were not. In the pectoral fins, the hoxAa cluster was expressed in split domains and did not obey the rule. By contrast, those from the hoxAb and hoxDa clusters were expressed in a manner keeping the rule, i.e., an ancestral pattern similar to those of tetrapods. It is plausible that this differential expression of the two clusters is caused by changes occurred in global control regions after cluster duplications.     

Cis-regulatory properties of medaka synexpression groups

During embryogenesis, tissue specification is triggered by the expression of a unique combination of developmental genes and their expression in time and space is crucial for successful development. Synexpression groups are batteries of spatiotemporally co-expressed genes that act in shared biological processes through their coordinated expression. Although several synexpression groups have been described in numerous vertebrate species, the regulatory mechanisms that orchestrate their common complex expression pattern remain to be elucidated. Here we performed a pilot screen on 560 genes of the vertebrate model system medaka (Oryzias latipes) to systematically identify synexpression groups and investigate their regulatory properties by searching for common regulatory cues. We find that synexpression groups share DNA motifs that are arranged in various combinations into cis-regulatory modules that drive co-expression. In contrast to previous assumptions that these genes are located randomly in the genome, we discovered that genes belonging to the same synexpression group frequently occur in synexpression clusters in the genome. This work presents a first repertoire of synexpression group common signatures, a resource that will contribute to deciphering developmental gene regulatory networks.     

Draft genome of the medaka fish: A comprehensive resource for medaka developmental genetics and vertebrate evolutionary biology

The medaka Oryzias latipes is a small egg-laying freshwater teleost, and has become an excellent model system for developmental genetics and evolutionary biology. The medaka genome is relatively small in size, ∼800 Mb, and the genome sequencing project was recently completed by Japanese research groups, providing a high-quality draft genome sequence of the inbred Hd-rR strain of medaka. In this review, I present an overview of the medaka genome project including genome resources, followed by specific findings obtained with the medaka draft genome. In particular, I focus on the analysis that was done by taking advantage of the medaka system, such as the sex chromosome differentiation and the regional history of medaka species using single nucleotide polymorphisms as genomic markers.     

The DNA sequence of medaka chromosome LG22

We report the genomic DNA sequence of a single chromosome (linkage group 22; LG22) of the small teleost fish medaka (Oryzias latipes) as a first whole chromosome sequence from a non-mammalian vertebrate. The order and orientation of 633 protein-coding genes were deduced from 18,803,338 bp of DNA sequence, providing the opportunity to analyze chromosome evolution of vertebrate genomes by direct comparison with the human genome. The average number of genes in the "conserved gene cluster" (CGC), a strict definition of "synteny" at the sequence basis, between medaka and human was 1.6. These and other data suggest that approximately 38.8% of pair-wise gene relationships would have been broken from their common ancestor in the human and medaka lineages and further imply that approx 20,000 (15,520-23,280) breaks would have occurred from the entire genome of the common ancestor. These breaks were generated mainly by intra-chromosomal shufflings at a specific era in the vertebrate lineage. These precise comparative genomics allowed us to identify the pieces of ancient chromosomes of the common vertebrate ancestor and estimate chromosomal evolution in the vertebrate lineage.     

Production of a maternal-zygotic medaka mutant using hybrid sterility

Taking advantage of the characteristics that make hybrids between Japanese and Chinese medaka grow well, albeit sterile, we have developed a method of germ-line replacement in which these hybrids are used as hosts for the production of a maternal-zygotic mutant. The protocol is described herein.     

Geographic variation and diversity of the cytochrome b gene in Japanese wild populations of medaka, Oryzias latipes

We conducted a polymerase chain reaction--restriction fragment length polymorphism (PCR-RFLP) analysis of the mitochondrial cytochrome b gene to elucidate the detailed genetic population structure of Japanese wild populations of medaka, Oryzias latipes. The analysis of 1,225 specimens collected from 303 sites identified 67 mitotypes. Subsequently we determined the nucleotide sequences of the complete cytochrome b gene (1141-bp) to clarify the phylogenetic relationships among mitotypes. The phylogenetic tree based on nucleotide sequences indicated three major clades (A, B and C) that differed by 11.3-11.8%, corresponding to three clusters previously identified by RFLP analysis of entire mitochondrial DNAs. The geographic distribution of mitotypes in clades A and B was fully concordant with the Northern and Southern Populations defined by allozymes. Clade A could be subdivided into three subclades and clade B into eleven, with sequence divergences among subclades of 1.3-5.8%. Each distribution of mitotypes in subclades roughly corresponded to that of mtDNA haplotypes in subclusters previously identified. Mitotypes in clade C were found only in the Kanto district. The phylogenetic relationships and the estimated divergence times suggest that three Japanese clades originated from a common ancestor and were separated during the Pliocene, and that the regional differentiation of subclades was closely connected with the geological history of the Quaternary. This study has also demonstrated the possibility of artificial disturbance of natural distribution especially in the Kanto district and the superior efficacy of PCR-RFLP analysis as a simple method for detecting genetic variation and artificial gene flow of medaka.     

Mutation in the abcb7 gene causes abnormal iron and fatty acid metabolism in developing medaka fish

The medaka fish (Oryzias latipes) is an emerging model organism for which a variety of unique developmental mutants have now been generated. Our recent mutagenesis screening of the medaka isolated a unique mutant that develops a fatty liver at larval stages. Positional cloning identified the responsible gene as medaka abcb7. Abcb7, a mitochondrial ABC (ATP binding cassette) half-transporter, has been implicated in iron metabolism. Recently, human Abcb7 was found to be mutated in X-linked sideroblastic anemia with cerebellar ataxia (XLSA/A). The homozygous medaka mutant exhibits abnormal iron metabolism in erythrocytes and accumulation of lipid in the liver. Microarray and in situ hybridization analyses demonstrated that the expression of genes involved in iron and lipid metabolisms are both affected in the mutant liver, suggesting novel roles of Abcb7 in the development of physiologically functional liver. The medaka abcb7 mutant thus could provide insights into the pathogenesis of XLSA/A as well as the normal function of the gene.