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.

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.

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.

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.

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.

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.

ZMYND10 functions in a chaperone relay during axonemal dynein assembly.

Molecular chaperones promote the folding and macromolecular assembly of a diverse set of 'client' proteins. How ubiquitous chaperone machineries direct their activities towards specific sets of substrates is unclear. Through the use of mouse genetics, imaging and quantitative proteomics we uncover that ZMYND10 is a novel co-chaperone that confers specificity for the FKBP8-HSP90 chaperone complex towards axonemal dynein clients required for cilia motility. Loss of ZMYND10 perturbs the chaperoning of axonemal dynein heavy chains, triggering broader degradation of dynein motor subunits. We show that pharmacological inhibition of FKBP8 phenocopies dynein motor instability associated with the loss of ZMYND10 in airway cells and that human disease-causing variants of ZMYND10 disrupt its ability to act as an FKBP8-HSP90 co-chaperone. Our study indicates that primary ciliary dyskinesia (PCD), caused by mutations in dynein assembly factors disrupting cytoplasmic pre-assembly of axonemal dynein motors, should be considered a cell-type specific protein-misfolding disease.

Complete fusion of a transposon and herpesvirus created the Teratorn mobile element in medaka fish.

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.

Effect of Stress-Free Therapy on immune system: Induction of Interleukin 10 expression in lymphocytes through activation of CD19(+) CD24(hi) CD38(hi) regulatory B Cells.

BACKGROUND AND AIMS: Mild thermal treatment with "Pinpoint Plantar Long-wavelength Infrared Light Irradiation (PP-LILI)" named as Stress-Free Therapy(®) increases peripheral-deep body temperature and blood flow, and improves multiple disorders including hyperpiesia, type II diabetes and cardiovascular patients. Immunomodulatory effects of PP-LILI were investigated. MATERIALS AND METHODS: Seven healthy individuals and 4 people with underlying medical condition (UMC) participated in this study. Participants were given PP-LILI stimuli twice a week over 3 weeks and followed with placebo stimuli over 3 weeks. This set of sessions was repeated 3 times. For analyses, fresh peripheral mononuclear cells from participants were stained with fluorescencedye conjugated monoclonal antibodies and changes in populational compositions and IL-10 expression levels were observed by flow cytometry. RESULTS: Distinct expression of IL-10 in lymphocytes was induced by PP-LILI from the second session in the healthy individuals. This induction was terminated during the following placebo sessions. PP-LILI induced activation of CD19(+) CD24(hi) CD38(hi) regulatory B cells in every session prior to induce the IL-10 in major lymphocytes. Activated regulatory B cells in the individuals with UMC decreased as same levels of healthy individuals after second PP-LILI session and re-activated with the stimuli. Significant population changes in neither regulatory T cells nor proinflammatory IL-17A expressing CD4(+) T cells were observed. CONCLUSIONS: PP-LILI is a potent immunomodulatory inducer that activates regulatory B cells and consequent IL-10 expression in lymphocytes. Moreover, its stimulatory intervals down-regulate the higher activation of regulatory B cells and lymphocyte's IL-10 expression occurred by UMC to the healthy people's level.

Unexpected link between polyketide synthase and calcium carbonate biomineralization.

INTRODUCTION: Calcium carbonate biominerals participate in diverse physiological functions. Despite intensive studies, little is known about how mineralization is initiated in organisms. RESULTS: We analyzed the medaka spontaneous mutant, ha, defective in otolith (calcareous ear stone) formation. ha lacks a trigger for otolith mineralization, and the causative gene was found to encode polyketide synthase (pks), a multifunctional enzyme mainly found in bacteria, fungi, and plant. Subsequent experiments demonstrate that the products of medaka PKS, most likely polyketides or their derivatives, act as nucleation facilitators in otolith mineralization. The generality of this novel PKS function is supported by the essential role of echinoderm PKS in calcareous skeleton formation together with the presence of PKSs in a much wider range of animals from coral to vertebrates. CONCLUSION: The present study first links PKS to biomineralization and provides a genetic cue for biogeochemistry of carbon and calcium cycles.

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.

Genetic control of startle behavior in medaka fish.

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.

A neural mechanism underlying mating preferences for familiar individuals in medaka fish.

Social familiarity affects mating preference among various vertebrates. Here, we show that visual contact of a potential mating partner before mating (visual familiarization) enhances female preference for the familiarized male, but not for an unfamiliarized male, in medaka fish. Terminal-nerve gonadotropin-releasing hormone 3 (TN-GnRH3) neurons, an extrahypothalamic neuromodulatory system, function as a gate for activating mating preferences based on familiarity. Basal levels of TN-GnRH3 neuronal activity suppress female receptivity for any male (default mode). Visual familiarization facilitates TN-GnRH3 neuron activity (preference mode), which correlates with female preference for the familiarized male. GnRH3 peptides, which are synthesized specifically in TN-GnRH3 neurons, are required for the mode-switching via self-facilitation. Our study demonstrates the central neural mechanisms underlying the regulation of medaka female mating preference based on visual social familiarity.

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.

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.

Trunk exoskeleton in teleosts is mesodermal in origin.

The vertebrate mineralized skeleton is known to have first emerged as an exoskeleton that extensively covered the fossil jawless fish. The evolutionary origin of this exoskeleton has long been attributed to the emergence of the neural crest, but experimental evaluation for this is still poor. Here we determine the embryonic origin of scales and fin rays of medaka (teleost trunk exoskeletons) by applying long-term cell labelling methods, and demonstrate that both tissues are mesodermal in origin. Neural crest cells, however, fail to contribute to these tissues. This result suggests that the trunk neural crest has no skeletogenic capability in fish, instead highlighting the dominant role of the mesoderm in the evolution of the trunk skeleton. This further implies that the role of the neural crest in skeletogenesis has been predominant in the cephalic region from the early stage of vertebrate evolution.

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.

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.

The medaka zic1/zic4 mutant provides molecular insights into teleost caudal fin evolution.

Teleosts have an asymmetrical caudal fin skeleton formed by the upward bending of the caudal-most portion of the body axis, the ural region. This homocercal type of caudal fin ensures powerful and complex locomotion and is regarded as one of the most important innovations for teleosts during adaptive radiation in an aquatic environment. However, the mechanisms that create asymmetric caudal fin remain largely unknown. The spontaneous medaka (teleost fish) mutant, Double anal fin (Da), exhibits a unique symmetrical caudal skeleton that resembles the diphycercal type seen in Polypterus and Coelacanth. We performed a detailed analysis of the Da mutant to obtain molecular insight into caudal fin morphogenesis. We first demonstrate that a large transposon, inserted into the enhancer region of the zic1 and zic4 genes (zic1/zic4) in Da, is associated with the mesoderm-specific loss of their transcription. We then show that zic1/zic4 are strongly expressed in the dorsal part of the ural mesenchyme and thereby induce asymmetric caudal fin development in wild-type embryos, whereas their expression is lost in Da. Comparative analysis further indicates that the dorsal mesoderm expression of zic1/zic4 is conserved in teleosts, highlighting the crucial role of zic1/zic4 in caudal fin morphogenesis.

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.