Immediate-early genes (IEGs) are useful for mapping active brain regions in various vertebrates. Here we identified a c-fos homologue gene in medaka and demonstrated that the amounts of c-fos transcripts and proteins in the medaka brain increased in relation to an artificially evoked seizure, suggesting that the homologue gene has the characteristics of IEGs, which are used as markers of neural activity. Next, quantitative reverse-transcription-polymerase chain reaction revealed that female mating behaviors upregulated c-fos transcription in some brain regions including the telencephalon, optic tectum, and cerebellum. In addition, we performed in situ hybridization with a c-fos intron probe to detect the de novo synthesis of c-fos transcripts and confirmed induction of c-fos transcription in these brain regions after mating. This is the first report of IEG induction in response to mating stimuli in teleost fish. Our results indicated that c-fos expression was induced in response to behavioral stimuli in the medaka brain and that medaka c-fos could be a useful marker of neural activity.

Establishment of robust gene expression boundary is crucial for creating elaborate morphology during development. However, mechanisms underlying boundary formation have been extensively studied only in a few model systems. We examined the establishment of zic1/zic4-expression boundary demarcating dorsoventral boundary of the entire trunk of medaka fish (Oryzias latipes) and identified a subgroup of dermomyotomal cells called horizontal boundary cells (HBCs) as crucial players for the boundary formation. Embryological and genetic analyses demonstrated that HBCs play crucial roles in the two major events of the process, i.e., refinement and maintenance. In the refinement, HBCs could serve as a chemical barrier against Wnts from the neural tube by expressing Hhip. At later stages, HBCs participate in the maintenance of the boundary by differentiating into the horizontal myoseptum physically inhibiting cell mixing across the boundary. These findings reveal the mechanisms underlying the dorsoventral boundary in the teleost trunk by specialized boundary cells.

Mobile genetic elements (e.g., transposable elements and viruses) display significant diversity with various life cycles, but how novel elements emerge remains obscure. Here, we report a giant (180-kb long) transposon, Teratorn, originally identified in the genome of medaka, Oryzias latipes. Teratorn belongs to the piggyBac superfamily and retains the transposition activity. Remarkably, Teratorn is largely derived from a herpesvirus of the Alloherpesviridae family that could infect fish and amphibians. Genomic survey of Teratorn-like elements reveals that some of them exist as a fused form between piggyBac transposon and herpesvirus genome in teleosts, implying the generality of transposon-herpesvirus fusion. We propose that Teratorn was created by a unique fusion of DNA transposon and herpesvirus, leading to life cycle shift. Our study supports the idea that recombination is the key event in generation of novel mobile genetic elements. Teratorn is a large mobile genetic element originally identified in the small teleost fish medaka. Here, the authors show that Teratorn is derived from the fusion of a piggyBac superfamily DNA transposon and an alloherpesvirus and that it is widely found across teleost fish.

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110-200 million years of evolution.

Medaka (Oryzias latipes) is a small freshwater teleost that provides an excellent developmental genetic model complementary to zebrafish. Our recent mutagenesis screening using medaka identified headfish (hdf) which is characterized by the absence of trunk and tail structures with nearly normal head including the midbrain-hindbrain boundary (MHB). Positional-candidate cloning revealed that the hdf mutation causes a functionally null form of Fgfr1. The fgfr1hdf is thus the first fgf receptor mutant in fish. Although FGF signaling has been implicated in mesoderm induction, mesoderm is induced normally in the fgfr1hdf mutant, but subsequently, mutant embryos fail to maintain the mesoderm, leading to defects in mesoderm derivatives, especially in trunk and tail. Furthermore, we found that morpholino knockdown of medaka fgf8 resulted in a phenotype identical to the fgfr1hdf mutant, suggesting that like its mouse counterpart, Fgf8 is a major ligand for Fgfr1 in medaka early embryogenesis. Intriguingly, Fgf8 and Fgfr1 in zebrafish are also suggested to form a major ligand-receptor pair, but their function is much diverged, as the zebrafish fgfr1 morphant and zebrafish fgf8 mutant acerebellar (ace) only fail to develop the MHB, but develop nearly unaffected trunk and tail. These results provide evidence that teleost fish have evolved divergent functions of Fgf8-Fgfr1 while maintaining the ligand-receptor relationships. Comparative analysis using different fish is thus invaluable for shedding light on evolutionary diversification of gene function.

5-methyl-cytosines at CpG sites frequently mutate into thymines, accounting for a large proportion of spontaneous point mutations. The repair system would leave substantial numbers of errors in neighboring regions if the synthesis of erased gaps around deaminated 5-methyl-cytosines is error-prone. Indeed, we identified an unexpected genome-wide role of the CpG methylation state as a major determinant of proximal natural genetic variation. Specifically, 507 Mbp (∼18%) of the human genome was within 10 bp of a CpG site; in these regions, the single nucleotide polymorphism (SNP) rate significantly increased by ∼50% (P < 10(-566) by a two-proportion z-test) if the neighboring CpG sites are methylated. To reconfirm this finding in another vertebrate, we compared six single-base resolution methylomes in two inbred medaka (Oryzias latipes) strains with sufficient genetic divergence (3.4%). We found that the SNP rate also increased by ∼50% (P < 10(-2170)), and the substitution rates in all dinucleotides increased simultaneously (P < 10(-441)) around methylated CpG sites. In the hypomethylated regions, the "CGCG" motif was significantly enriched (P < 10(-680)) and evolutionarily conserved (P = ∼ 0.203%), and slow CpG deamination rather than fast CpG gain was seen, indicating a possible role of CGCG as a candidate cis-element for the hypomethylation state. In regions that were hypermethylated in germline-like tissues but were hypomethylated in somatic liver cells, the SNP rate was significantly smaller than that in hypomethylated regions in both tissues, suggesting a positive selective pressure during DNA methylation reprogramming. This is the first report of findings showing that the CpG methylation state is significantly correlated with the characteristics of evolutionary change in neighboring DNA.

Genetic polymorphisms are thought to generate intraspecific behavioral diversities, both within and among populations. The mechanisms underlying genetic control of behavioral properties, however, remain unclear in wild-type vertebrates, including humans. To explore this issue, we used diverse inbred strains of medaka fish (Oryzias latipes) established from the same and different local populations. Medaka exhibit a startle response to a visual stimulus (extinction of illumination) by rapidly bending their bodies (C-start) 20-ms after the stimulus presentation. We measured the rates of the response to repeated stimuli (1-s interval, 40 times) among four inbred strains, HNI-I, HNI-II, HO5, and Hd-rR-II1, and quantified two properties of the startle response: sensitivity (response rate to the first stimulus) and attenuation of the response probability with repeated stimulus presentation. Among the four strains, the greatest differences in these properties were detected between HNI-II and Hd-rR-II1. HNI-II exhibited high sensitivity (approximately 80%) and no attenuation, while Hd-rR-II1 exhibited low sensitivity (approximately 50%) and almost complete attenuation after only five stimulus presentations. Our findings suggested behavioral diversity of the startle response within a local population as well as among different populations. Linkage analysis with F2 progeny between HNI-II and Hd-rR-II1 detected quantitative trait loci (QTL) highly related to attenuation, but not to sensitivity, with a maximum logarithm of odds score of 11.82 on linkage group 16. The three genotypes (homozygous for HNI-II and Hd-rR-II1 alleles, and heterozygous) at the marker nearest the QTL correlated with attenuation. Our findings are the first to suggest that a single genomic region might be sufficient to generate individual differences in startle behavior between wild-type strains. Further identification of genetic polymorphisms that define the behavioral trait will contribute to our understanding of the neural mechanisms underlying behavioral diversity, allowing us to investigate the adaptive significance of intraspecific behavioral polymorphisms of the startle response.