CTCF looping is established during gastrulation in medaka embryos.

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.

The art of medaka genetics and genomics: what makes them so unique?

The medaka fish, Oryzias latipes, is an emerging vertebrate model and now has a high quality draft genome and a number of unique mutants. The long history of medaka research in Japan has provided medaka with unique features, which are complementary to other vertebrate models. A large collection of spontaneous mutants collected over a century, the presence of highly polymorphic inbred lines established over decades, and the recently completed genome sequence all give the medaka a big boost. This review focuses on the state of the art in medaka genetics and genomics, such as the first isolation of active transposons in vertebrates, the influence of chromatin structure on sequence variation, fine quantitative trait locus (QTL) analysis, and versatile mutants as human disease models.

Modular development of the teleost trunk along the dorsoventral axis and zic1/zic4 as selector genes in the dorsal module.

Teleost fish exhibit remarkable diversity in morphology, such as fins and coloration, particularly on the dorsal side. These structures are evolutionary adaptive because their back is highly visible to other individuals. However, owing to the late phenotypic appearance (from larva to adult) and lack of appropriate mutants, the genetic mechanisms that regulate these dorsoventrally asymmetric external patterns are largely unknown. To address this, we have analyzed the spontaneous medaka mutant Double anal fin (Da), which exhibits a mirror-image duplication of the ventral half across the lateral midline from larva to adult. Da is an enhancer mutant for zic1 and zic4 in which their expression in dorsal somites is lost. We show that the dorsoventral polarity in Da somites is lost and then demonstrate using transplantation techniques that somites and their derived tissues globally determine the multiple dorsal-specific characteristics of the body (fin morphology and pigmentation) from embryo to adult. Intriguingly, the zic1/zic4 expression in the wild type persists throughout life in the dorsal parts of somite derivatives, i.e. the myotome, dermis and vertebrae, forming a broad dorsal domain in the trunk. Comparative analysis further implies a central role for zic1/zic4 in morphological diversification of the teleost body. Taken together, we propose that the teleost trunk consists of dorsal/ventral developmental modules and that zic1/zic4 in somites function as selector genes in the dorsal module to regulate multiple dorsal morphologies.

Associations between nucleosome phasing, sequence asymmetry, and tissue-specific expression in a set of inbred Medaka species.

BACKGROUND: Transcription start sites (TSSs) with pronounced and phased nucleosome arrays downstream and nucleosome-depleted regions upstream of TSSs are observed in various species. RESULTS: We have characterized sequence variation and expression properties of this set of TSSs (which we call "Nucleocyclic TSSs") using germline and somatic cells of three medaka (Oryzias latipes) inbred isolates from different locations. We found nucleocyclic TSSs in medaka to be associated with higher gene expression and characterized by a clear boundary in sequence composition with potentially-nucleosome-destabilizing A/T-enrichment upstream (p < 10(-60)) and nucleosome- accommodating C/G-enrichment downstream (p < 10(-40)) that was highly conserved from an ancestor. A substantial genetic distance between the strains facilitated the in-depth analysis of patterns of fixed mutations, revealing a localization-specific equilibrium between the rates of distinct mutation categories that would serve to maintain the conserved sequence anisotropy around TSSs. Downstream of nucleocyclic TSSs, C to T, T to C, and other mutation rates on the sense strand increased around first nucleosome dyads and decreased around first linkers, which contrasted with genomewide mutational patterns around nucleosomes (p < 5 %). C to T rates are higher than G to A rates around nucleosome associated with germline nucleocyclic TSS sites (p < 5 %), potentially due to the asymmetric effect of transcription-coupled repair. CONCLUSIONS: Our results demonstrate an atypical evolutionary process surrounding nucleocyclic TSSs.

Genome-wide genetic variations are highly correlated with proximal DNA methylation patterns.

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.

The establishment of rotational polarity in the airway and ependymal cilia: analysis with a novel cilium motility mutant mouse.

The airway is covered by multicilia that beat in a metachronous manner toward the mouth to eliminate debris and infectious particles. Coordinated one-directional beating is an essential feature of multicilia in the airway to guarantee proper mucociliary clearance. Defects in ciliary motility lead to primary ciliary dyskinesia (PCD), with major symptoms including bronchitis and other chronic respiratory diseases. Recent work suggested that ciliary motility and planar polarity are required in the process of ciliary alignment that produces coordinated beating. However, the extent to which cilia motility is involved in this process in mammals has not yet been fully clarified. Here, to address the role of ciliary motility in the process of coordinated ciliary alignment, we analyzed Kintoun mice mutants (Ktu(-/-)). Ktu(-/-) exhibited typical phenotypes of PCD with complete loss of ciliary motility in trachea and another ciliated tissue, the brain ependyma. Immunohistochemistry using antibodies against axonemal dynein confirmed the loss of multiple axonemal dynein components in mutant cilia. Observation of cilia orientation based on basal foot directions revealed that ciliary motility was not required in the alignment of airway cilia, whereas a strong requirement was observed in brain ependymal cells. Thus we conclude that the involvement of ciliary motility in the establishment of coordinated ciliary alignment varies among tissues.

The medaka draft genome and insights into vertebrate genome evolution.

Teleosts comprise more than half of all vertebrate species and have adapted to a variety of marine and freshwater habitats. Their genome evolution and diversification are important subjects for the understanding of vertebrate evolution. Although draft genome sequences of two pufferfishes have been published, 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 and developmental genetics. In the assembled medaka genome (700 megabases), which is less than half of the zebrafish genome, we predicted 20,141 genes, including approximately 2,900 new genes, using 5'-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 approximately 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.

Visualization of Actin Cytoskeleton in Cellular Protrusions in Medaka Embryos.

Cellular protrusions are fundamental structures for a wide variety of cellular behaviors, such as cell migration, cell-cell interaction, and signal reception. Visualization of cellular protrusions in living cells can be achieved by labeling of cytoskeletal actin with genetically encoded fluorescent probes. Here, we describe a detailed experimental procedure to visualize cellular protrusions in medaka embryos, which consists of the following steps: preparation of Actin-Chromobody-GFP and α-bungarotoxin mRNAs for actin labeling and immobilization of the embryo, respectively; microinjection of the mRNAs into embryos in a mosaic fashion to sparsely label individual cells; removal of the hard chorion, which hampers observation; and visualization of cellular protrusions in the embryo with a confocal microscope. Overall, our protocol provides a simple method to reveal cellular protrusions in vivo by confocal microscopy.

Identification of an adverse outcome pathway (AOP) for chemical-induced craniofacial anomalies using the transgenic zebrafish model.

Craniofacial anomalies are one of the most frequent birth defects worldwide and are often caused by genetic and environmental factors such as pharmaceuticals and chemical agents. Although identifying adverse outcome pathways (AOPs) is a central issue for evaluating the teratogenicity, the AOP causing craniofacial anomalies has not been identified. Recently, zebrafish has gained interest as an emerging model for predicting teratogenicity because of high throughput, cost-effectiveness and availability of various tools for examining teratogenic mechanisms. Here, we established zebrafish sox10-EGFP reporter lines to visualize cranial neural crest cells (CNCCs) and have identified the AOPs for craniofacial anomalies. When we exposed the transgenic embryos to teratogens that were reported to cause craniofacial anomalies in mammals, CNCC migration and subsequent morphogenesis of the first pharyngeal arch were impaired at 24 hours post-fertilization. We also found that cell proliferation and apoptosis of the migratory CNCCs were disturbed, which would be key events of the AOP. From these results, we propose that our sox10-EGFP reporter lines serve as a valuable model for detecting craniofacial skeletal abnormalities, from early to late developmental stages. Given that the developmental process of CNCCs around this stage is highly conserved between zebrafish and mammals, our findings can be extrapolated to mammalian craniofacial development and thus help in predicting craniofacial anomalies in human.

Mechanically sensitive HSF1 is a key regulator of left-right symmetry breaking in zebrafish embryos.

The left-right symmetry breaking of vertebrate embryos requires nodal flow. However, the molecular mechanisms that mediate the asymmetric gene expression regulation under nodal flow remain elusive. Here, we report that heat shock factor 1 (HSF1) is asymmetrically activated in the Kupffer's vesicle of zebrafish embryos in the presence of nodal flow. Deficiency in HSF1 expression caused a significant situs inversus and disrupted gene expression asymmetry of nodal signaling proteins in zebrafish embryos. Further studies demonstrated that HSF1 is a mechanosensitive protein. The mechanical sensation ability of HSF1 is conserved in a variety of mechanical stimuli in different cell types. Moreover, cilia and Ca2+-Akt signaling axis are essential for the activation of HSF1 under mechanical stress in vitro and in vivo. Considering the conserved expression of HSF1 in organisms, these findings unveil a fundamental mechanism of gene expression regulation by mechanical clues during embryonic development and other physiological and pathological transformations.

Establishment and characterization of Roberts syndrome and SC phocomelia model medaka (Oryzias latipes).

Roberts syndrome and SC phocomelia (RBS/SC) are genetic autosomal recessive syndromes caused by establishment of cohesion 1 homolog 2 ( ESCO 2) mutation. RBS/SC appear to have a variety of clinical features, even with the same mutation of the ESCO2 gene. Here, we established and genetically characterized a medaka model of RBS/SC by reverse genetics. The RBS/SC model was screened from a mutant medaka library produced by the Targeting Induced Local Lesions in Genomes method. The medaka mutant carrying the homozygous mutation at R80S in the conserved region of ESCO2 exhibited clinical variety (i.e. developmental arrest with craniofacial and chromosomal abnormalities and embryonic lethality) as characterized in RBS/SC. Moreover, widespread apoptosis and downregulation of some gene expression, including notch1a, were detected in the R80S mutant. The R80S mutant is the animal model for RBS/SC and a valuable resource that provides the opportunity to extend knowledge of ESCO2. Downregulation of some gene expression in the R80S mutant is an important clue explaining non-correlation between genotype and phenotype in RBS/SC.

Wnt11 acts on dermomyotome cells to guide epaxial myotome morphogenesis.

The dorsal axial muscles, or epaxial muscles, are a fundamental structure covering the spinal cord and vertebrae, as well as mobilizing the vertebrate trunk. To date, mechanisms underlying the morphogenetic process shaping the epaxial myotome are largely unknown. To address this, we used the medaka zic1/zic4-enhancer mutant Double anal fin (Da), which exhibits ventralized dorsal trunk structures resulting in impaired epaxial myotome morphology and incomplete coverage over the neural tube. In wild type, dorsal dermomyotome (DM) cells reduce their proliferative activity after somitogenesis. Subsequently, a subset of DM cells, which does not differentiate into the myotome population, begins to form unique large protrusions extending dorsally to guide the epaxial myotome dorsally. In Da, by contrast, DM cells maintain the high proliferative activity and mainly form small protrusions. By combining RNA- and ChIP-sequencing analyses, we revealed direct targets of Zic1, which are specifically expressed in dorsal somites and involved in various aspects of development, such as cell migration, extracellular matrix organization, and cell-cell communication. Among these, we identified wnt11 as a crucial factor regulating both cell proliferation and protrusive activity of DM cells. We propose that dorsal extension of the epaxial myotome is guided by a non-myogenic subpopulation of DM cells and that wnt11 empowers the DM cells to drive the coverage of the neural tube by the epaxial myotome.

Induction of c-fos transcription in the medaka brain (Oryzias latipes) in response to mating stimuli.

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.

Occurrence of a short variant of the Tol2 transposable element in natural populations of the medaka fish.

Tol2 is a member of the hAT (hobo/Activator/Tam3) transposable element family, residing as 10-30 copies per diploid genome in the medaka fish. We previously reported that this element is highly homogeneous in structure at both the restriction map level and the nucleotide sequence level. It was, however, possible that there is variation of such a low frequency as not to have been detected in our previous surveys, in which samples from 12 geographical locations were used. In the present study, we first conducted searches of genome sequence databases of medaka, and found a 119-bp-long internal deletion. We then conducted a survey of samples from 58 locations for this deletion by performing PCR preceded by restriction enzyme digestion to increase the sensitivity to this deletion. We found that copies suffering this deletion have spread, or have been generated by multiple origins, in the northern-to-central part of mainland Japan. Thus, although the high homogeneity in structure is a distinct feature of Tol2, variation does exist at low frequencies in natural populations of medaka. The current status of Tol2 is expected to provide information with which results of future surveys can be compared for clarification of determinants of population dynamics of this DNA-based element.

UTGB/medaka: genomic resource database for medaka biology.

Medaka (Oryzias latipes) is a small egg-laying freshwater teleost native to East Asia that has become an excellent model system for developmental genetics and evolutionary biology. The draft medaka genome sequence (700 Mb) was reported in June 2007, and its substantial genomic resources have been opened to the public through the University of Tokyo Genome Browser Medaka (UTGB/medaka) database. This database provides basic genomic information, such as predicted genes, expressed sequence tags (ESTs), guanine/cytosine (GC) content, repeats and comparative genomics, as well as unique data resources including (i) 2473 genetic markers and experimentally confirmed PCR primers that amplify these markers, (ii) 142,414 bacterial artificial chromosome (BAC) and 217,344 fosmid end sequences that amount to 15.0- and 11.1-fold clone coverage of the entire genome, respectively, and were used for draft genome assembly, (iii) 16,519,460 single nucleotide polymorphisms (SNPs), and 2 859 905 insertions/deletions detected between two medaka inbred strain genomes and (iv) 841 235 5'-end serial analyses of gene-expression (SAGE) tags that identified 344 266 transcription start sites on the genome. UTGB/medaka is available at: http://medaka.utgenome.org/.

Characterization of teleost Mdga1 using a gene-trap approach in medaka (Oryzias latipes).

MAM domain containing glycosilphosphatidilinositol anchor 1 (MDGA1) is an IgCAM protein present in many vertebrate species including humans. In mammals, MDGA1 is expressed by a subset of neurons in the developing brain and thought to function in neural cell migration. We identified a fish ortholog of mdga1 by a gene-trap screen utilizing the Frog Prince transposon in medaka (Japanese killifish, Oryzias latipes). The gene-trap vector was inserted into an intronic region of mdga1 to form a chimeric protein with green fluorescent protein, allowing us to monitor mdga1 expression in vivo. Expression of medaka mdga1 was seen in various types of embryonic brain neurons, and specifically in neurons migrating toward their target sites, supporting the proposed function of MDGA1. We also isolated the closely related mdga2 gene, whose expression partially overlapped with that of mdga1. Despite the fact that the gene-trap event eliminated most of the functional domains of the Mdga1 protein, homozygous embryos developed normally without any morphological abnormality, suggesting a functional redundancy of Mdga1 with other related proteins. High sequential homology of MDGA proteins between medaka and other vertebrate species suggests an essential role of the MDGA gene family in brain development among the vertebrate phylum.

Rescue from oculocutaneous albinism type 4 using medaka slc45a2 cDNA driven by its own promoter.

Patients and vertebrate mutants with oculocutaneous albinism type 4 (OCA4) have mutations in the solute carrier family 45 member 2 (slc45a2) gene. However, there is no empirical evidence for this gene-phenotype relationship. There is a unique OCA4 mutant in medaka (b) that exhibits albinism only in the skin, but the mechanism underlying this phenotype is also unknown. In this study, we rescued medaka OCA4 phenotypes, in both the eyes and the skin, by micro-injection of an slc45a2-containing genomic fragment or slc45a2 cDNA driven by its own 0.9-kb promoter. We also identified a spontaneous nucleotide change of 339 bp in the promoter as the b mutation. There are multiple transcription start sites in medaka slc45a2, as in its human ortholog, and only the shortest and eye-specific mRNA is transcribed with the b mutation. Interestingly, we further revealed a conserved pyrimidine (Py)-rich sequence of approximately 10 bp in the promoter by medaka-pufferfish comparative genomics and verified that it plays an indispensable role for expression of slc45a2 in the skin. Further studies of the 0.9-kb promoter identified in this study should provide insights into the cis/trans-regulatory mechanisms underlying the ocular and cutaneous expression of slc45a2.

Distribution of complete and defective copies of the Tol1 transposable element in natural populations of the medaka fish Oryzias latipes.

DNA-based transposable elements are present in the genomes of various organisms, and generally occur in autonomous and nonautonomous forms, with a good correspondence to complete and defective copies, respectively. In vertebrates, however, the vast majority of DNA-based elements occur only in the nonautonomous form. Until now, the only clear exception known has been the Tol2 element of the medaka fish, which still causes mutations in genes of the host species. Here, we report another exception: the Tol1 element of the same species. This element was thought likely to be a "dead" element like the vast majority of vertebrate elements, but recent identification of an autonomous Tol1 copy in a laboratory medaka strain gave rise to the possibility that the element is still "alive" in medaka natural populations. We examined variation in the structure of Tol1 copies through genomic Southern blot analysis, and revealed that 10 of the 32 fish samples examined contained full-length Tol1 copies in their genomes. The frequency at which these copies occur among Tol1 copies is at most 0.5%, yet some of them still have the ability to produce a functional transposase. The medaka fish thus harbors two active DNA-based elements in its genome, and is in this respect unique among vertebrates.

Horizontal Boundary Cells, a Special Group of Somitic Cells, Play Crucial Roles in the Formation of Dorsoventral Compartments in Teleost Somite.

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.

Genetic analysis of craniofacial traits in the medaka.

Family and twin studies suggest that a substantial genetic component underlies individual differences in craniofacial morphology. In the current study, we quantified 444 craniofacial traits in 100 individuals from two inbred medaka (Oryzias latipes) strains, HNI and Hd-rR. Relative distances between defined landmarks were measured in digital images of the medaka head region. A total of 379 traits differed significantly between the two strains, indicating that many craniofacial traits are controlled by genetic factors. Of these, 89 traits were analyzed via interval mapping of 184 F(2) progeny from an intercross between HNI and Hd-rR. We identified quantitative trait loci for 66 craniofacial traits. The highest logarithm of the odds score was 6.2 for linkage group (LG) 9 and 11. Trait L33, which corresponds to the ratio of head length to head height at eye level, mapped to LG9; trait V15, which corresponds to the ratio of snout length to head width measured behind the eyes, mapped to LG11. Our initial results confirm the potential of the medaka as a model system for the genetic analysis of complex traits such as craniofacial morphology.