Teleosts comprise more than half of all vertebrate species and have adapted to a variety of marine and freshwater habitats 1 . Their genome evolution and diversification are important subjects for the understanding of vertebrate evolution. Although draft genome sequences of two pufferfishes have been published 2,3 , analysis of more fish genomes is desirable. Here we report a high-quality draft genome sequence of a small egg-laying freshwater teleost, medaka (Oryzias latipes). Medaka is native to East Asia and an excellent model system for a wide range of biology, including ecotoxicology, carcinogenesis, sex determination 4-6 and developmental genetics 7 . In the assembled medaka genome (700 megabases), which is less than half of the zebrafish genome, we predicted 20,141 genes, including 2,900 new genes, using 59-end serial analysis of gene expression tag information. We found single nucleotide polymorphisms (SNPs) at an average rate of 3.42% between the two inbred strains derived from two regional populations; this is the highest SNP rate seen in any vertebrate species. Analyses based on the dense SNP information show a strict genetic separation of 4 million years (Myr) between the two populations, and suggest that differential selective pressures acted on specific gene categories. Four-way comparisons with the human, pufferfish (Tetraodon), zebrafish and medaka genomes revealed that eight major interchromosomal rearrangements took place in a remarkably short period of 50 Myr after the whole-genome duplication event in the teleost ancestor and afterwards, intriguingly, the medaka genome preserved its ancestral karyotype for more than 300 Myr.We applied the whole-genome shotgun approach to an inbred strain, , derived from the southern Japanese population, as the main target. A total of 13.8 million reads amounting to approximately 10.6-fold genome coverage were obtained from the shotgun plasmid, fosmid and bacterial artificial chromosome (BAC) libraries. A newly developed RAMEN assembler was used to process the shotgun reads to generate contigs and scaffolds. The N50 values (50% of nucleotides in an assembly are in scaffolds-or contigs-longer than or equal to the N50 value) are ,1.41 megabases (Mb) for scaffolds and ,9.8 kilobases (Kb) for contigs. The total length of the contigs reached 700.4 Mb, which, from now on, we refer to as the medaka genome size.To construct ultracontigs, the scaffolds were integrated with the medaka genetic map by using SNP markers. For this purpose, we further obtained about 2.8-fold coverage of shotgun reads from another inbred strain HNI (refs 9, 10), which is derived from the northern Japanese population. The reads were assembled by RAMEN to scaffolds covering 648 Mb. Aligning the HNI contigs with the HdrR genome using BLASTZ 11 , we identified 16.4 million SNPs as well as 1.40 million insertions and 1.45 million deletions in non-repetitive regions (Supplementary Table 2). We selected 2,401 SNPs and genetically mapped them onto medaka chromosomes using a backcross panel between the...

Among mammals, only eutherians and marsupials are viviparous and have genomic imprinting that leads to parent-of-origin-specific differential gene expression. We used comparative analysis to investigate the origin of genomic imprinting in mammals. PEG10 (paternally expressed 10) is a retrotransposon-derived imprinted gene that has an essential role for the formation of the placenta of the mouse. Here, we show that an orthologue of PEG10 exists in another therian mammal, the marsupial tammar wallaby (Macropus eugenii), but not in a prototherian mammal, the egg-laying platypus (Ornithorhynchus anatinus), suggesting its close relationship to the origin of placentation in therian mammals. We have discovered a hitherto missing link of the imprinting mechanism between eutherians and marsupials because tammar PEG10 is the first example of a differentially methylated region (DMR) associated with genomic imprinting in marsupials. Surprisingly, the marsupial DMR was strictly limited to the 5′ region of PEG10, unlike the eutherian DMR, which covers the promoter regions of both PEG10 and the adjacent imprinted gene SGCE. These results not only demonstrate a common origin of the DMR-associated imprinting mechanism in therian mammals but provide the first demonstration that DMR-associated genomic imprinting in eutherians can originate from the repression of exogenous DNA sequences and/or retrotransposons by DNA methylation.

BackgroundThe Solanaceae family includes several economically important vegetable crops. The tomato (Solanum lycopersicum) is regarded as a model plant of the Solanaceae family. Recently, a number of tomato resources have been developed in parallel with the ongoing tomato genome sequencing project. In particular, a miniature cultivar, Micro-Tom, is regarded as a model system in tomato genomics, and a number of genomics resources in the Micro-Tom-background, such as ESTs and mutagenized lines, have been established by an international alliance.ResultsTo accelerate the progress in tomato genomics, we developed a collection of fully-sequenced 13,227 Micro-Tom full-length cDNAs. By checking redundant sequences, coding sequences, and chimeric sequences, a set of 11,502 non-redundant full-length cDNAs (nrFLcDNAs) was generated. Analysis of untranslated regions demonstrated that tomato has longer 5'- and 3'-untranslated regions than most other plants but rice. Classification of functions of proteins predicted from the coding sequences demonstrated that nrFLcDNAs covered a broad range of functions. A comparison of nrFLcDNAs with genes of sixteen plants facilitated the identification of tomato genes that are not found in other plants, most of which did not have known protein domains. Mapping of the nrFLcDNAs onto currently available tomato genome sequences facilitated prediction of exon-intron structure. Introns of tomato genes were longer than those of Arabidopsis and rice. According to a comparison of exon sequences between the nrFLcDNAs and the tomato genome sequences, the frequency of nucleotide mismatch in exons between Micro-Tom and the genome-sequencing cultivar (Heinz 1706) was estimated to be 0.061%.ConclusionThe collection of Micro-Tom nrFLcDNAs generated in this study will serve as a valuable genomic tool for plant biologists to bridge the gap between basic and applied studies. The nrFLcDNA sequences will help annotation of the tomato whole-genome sequence and aid in tomato functional genomics and molecular breeding. Full-length cDNA sequences and their annotations are provided in the database KaFTom http://www.pgb.kazusa.or.jp/kaftom/ via the website of the National Bioresource Project Tomato http://tomato.nbrp.jp.

JapanCommonly used classical inbred mouse strains have mosaic genomes with sequences from different subspecific origins. Their genomes are derived predominantly from the Western European subspecies Mus musculus domesticus, with the remaining sequences derived mostly from the Japanese subspecies Mus musculus molossinus. However, it remains unknown how this intersubspecific genome introgression occurred during the establishment of classical inbred strains. In this study, we resequenced the genomes of two M. m. molossinus-derived inbred strains, MSM/Ms and JF1/Ms. MSM/Ms originated from Japanese wild mice, and the ancestry of JF1/Ms was originally found in Europe and then transferred to Japan. We compared the characteristics of these sequences to those of the C57BL/6J reference sequence and the recent data sets from the resequencing of 17 inbred strains in the Mouse Genome Project (MGP), and the results unequivocally show that genome introgression from M. m. molossinus into M. m. domesticus provided the primary framework for the mosaic genomes of classical inbred strains. Furthermore, the genomes of C57BL/6J and other classical inbred strains have long consecutive segments with extremely high similarity (>99.998%) to the JF1/Ms strain. In the early 20th century, Japanese waltzing mice with a morphological phenotype resembling that of JF1/Ms mice were often crossed with European fancy mice for early studies of ''Mendelism,'' which suggests that the ancestor of the extant JF1/Ms strain provided the origin of the M. m. molossinus genome in classical inbred strains and largely contributed to its intersubspecific genome diversity.[Supplemental material is available for this article.]Classical inbred strains of mice were established in America in the early 20th century from European-derived fancy mice that were reared as pets (Morse 1981;Beck et al. 2000). It is well established that extant classical inbred strains are hybrids between multiple subspecies of Mus musculus, and they have a mosaic genome architecture comprised of sequences originating from different subspecies (Wade et al. 2002;Frazer et al. 2007;Yang et al. 2007). Recent high-resolution single nucleotide polymorphism (SNP) genotyping of wild-caught mice and a comparison of these sequences to those of classical inbred strains revealed that classical inbred strains are derived from a relatively small pool of fancy mice with limited haplotype diversity; their genomes are overwhelmingly derived from the Western European subspecies M. musculus domesticus, with the remaining sequences mostly derived from the Japanese subspecies M. musculus molossinus (Yang et al. 2011).Although detailed genealogical information is of great help for discovering the genes responsible for specific phenotypes, very little is known about how the intersubspecific genome introgression from M. m. molossinus into M. m. domesticus occurred during the establishment of classical inbred strains. In this study, we traced the ancestries of classical inbred strains, focusing on the contribution ...

BackgroundBody coloration is an ecologically important trait that is often involved in prey-predator interactions through mimicry and crypsis. Although this subject has attracted the interest of biologists and the general public, our scientific knowledge on the subject remains fragmentary. In the caterpillar of the swallowtail butterfly Papilio xuthus, spectacular changes in the color pattern are observed; the insect mimics bird droppings (mimetic pattern) as a young larva, and switches to a green camouflage coloration (cryptic pattern) in the final instar. Despite the wide variety and significance of larval color patterns, few studies have been conducted at a molecular level compared with the number of studies on adult butterfly wing patterns.ResultsTo obtain a catalog of genes involved in larval mimetic and cryptic pattern formation, we constructed expressed sequence tag (EST) libraries of larval epidermis for P. xuthus, and P. polytes that contained 20,736 and 5,376 clones, respectively, representing one of the largest collections available in butterflies. A comparison with silkworm epidermal EST information revealed the high expression of putative blue and yellow pigment-binding proteins in Papilio species. We also designed a microarray from the EST dataset information, analyzed more than five stages each for six markings, and confirmed spatial expression patterns by whole-mount in situ hybridization. Hence, we succeeded in elucidating many novel marking-specific genes for mimetic and cryptic pattern formation, including pigment-binding protein genes, the melanin-associated gene yellow-h3, the ecdysteroid synthesis enzyme gene 3-dehydroecdysone 3b-reductase, and Papilio-specific genes. We also found many cuticular protein genes with marking specificity that may be associated with the unique surface nanostructure of the markings. Furthermore, we identified two transcription factors, spalt and ecdysteroid signal-related E75, as genes expressed in larval eyespot markings. This finding suggests that E75 is a strong candidate mediator of the hormone-dependent coordination of larval pattern formation.ConclusionsThis study is one of the most comprehensive molecular analyses of complicated morphological features, and it will serve as a new resource for studying insect mimetic and cryptic pattern formation in general. The wide variety of marking-associated genes (both regulatory and structural genes) identified by our screening indicates that a similar strategy will be effective for understanding other complex traits.

Recent studies have revealed that a cilium-generated liquid flow in the node has a crucial role in the establishment of the left-right (LR) axis in the mouse. In fish, Kupffer's vesicle (KV), a teleost-specific spherical organ attached to the tail region, is known to have an equivalent role to the mouse node during LR axis formation. However, at present, there has been no report of an asymmetric gene expressed in KV under the control of fluid flow. Here we report the earliest asymmetric gene in teleost KV, medaka charon, and its regulation. Charon is a member of the Cerberus/DAN family of proteins, first identified in zebrafish. Although zebrafish charon was reported to be symmetrically expressed in KV, medaka charon displays asymmetric expression with more intense expression on the right side. This asymmetric expression was found to be regulated by KV flow because symmetric and up-regulated charon expression was observed in flow-defective embryos with immotile cilia or disrupted KV. Taken together, medaka charon is a reliable gene marker for LR asymmetry in KV and thus, will be useful for the analysis of the early steps downstream of the fluid flow.

The medaka is becoming an attractive model organism for the study of vertebrate early development and organogenesis and large-scale mutagenesis projects that are aimed at creating developmentally defective mutants are now being conducted by several groups in Japan. To strengthen the study of medaka developmental genetics, we have conducted a large-scale isolation of ESTs from medaka embryos and developed tools that facilitate mutant analysis. In this study, we have characterized a total of 132,082 sequences from both ends of cloned insert cDNAs from libraries generated at different stages of medaka embryo development. Clustering analysis with 3-prime sequences finally identified a total of 12,429 clusters. As a pilot analysis, 924 clusters were subjected to in situ hybridization to determine the spatial localization of their transcripts. Using EST sequence data generated in the present study, a 60-mer oligonucleotide microarray with 8,091 unigenes (Medaka Microarray 8K) was constructed and tested for its usefulness in expression profiling. Furthermore, we have developed a rapid and reliable mutant mapping system using a set of mapped EST markers (M-marker 2003) that covers the entire medaka genome. These resources will accelerate medaka mutant analyses and make an important contribution to the medaka genome project.

The proinflammatory cytokine interleukin 1β (IL-1β) induces prostaglandin E2 (PGE2) production via upregulation of cyclooxygenase-2 (COX-2) expression in synovial fibroblasts. This effect of IL-1β is involved in osteoarthritis. We investigated MAPK signaling pathways in IL-1β-induced COX-2 expression in feline synovial fibroblasts. In the presence of MAPK inhibitors, IL-1β-induced COX-2 expression and PGE2 release were both attenuated. IL-1β induced the phosphorylation of p38, JNK, MEK, and ERK1/2. A JNK inhibitor prevented not only JNK phosphorylation but also MEK and ERK1/2 phosphorylation in IL-1β-stimulated cells, but MEK and ERK1/2 inhibitors had no effect on JNK phosphorylation. A p38 inhibitor prevented p38 phosphorylation, but had no effect on MEK, ERK1/2, and JNK phosphorylation. MEK, ERK1/2, and JNK inhibitors had no effect on p38 phosphorylation. We also observed that in IL-1β-treated cells, phosphorylated MEK, ERK1/2, and JNK were co-precipitated with anti-phospho-MEK, ERK1/2, and JNK antibodies. The silencing of JNK1 in siRNA-transfected fibroblasts prevented IL-1β to induce phosphorylation of MEK and ERK1/2 and COX-2 mRNA expression. These observations suggest that JNK1 phosphorylation is necessary for the activation of the MEK/ERK1/2 pathway and the subsequent COX-2 expression for PGE2 release, and p38 independently contributes to the IL-1β effect in synovial fibroblasts.