Actin flows mediate a universal coupling between cell speed and cell persistence.

Cell movement has essential functions in development, immunity, and cancer. Various cell migration patterns have been reported, but no general rule has emerged so far. Here, we show on the basis of experimental data in vitro and in vivo that cell persistence, which quantifies the straightness of trajectories, is robustly coupled to cell migration speed. We suggest that this universal coupling constitutes a generic law of cell migration, which originates in the advection of polarity cues by an actin cytoskeleton undergoing flows at the cellular scale. Our analysis relies on a theoretical model that we validate by measuring the persistence of cells upon modulation of actin flow speeds and upon optogenetic manipulation of the binding of an actin regulator to actin filaments. Beyond the quantitative prediction of the coupling, the model yields a generic phase diagram of cellular trajectories, which recapitulates the full range of observed migration patterns.

Sexually dimorphic dynamics of the microtubule network in medaka (Oryzias latipes) germ cells.

Gametogenesis is the process through which germ cells differentiate into sexually dimorphic gametes, eggs and sperm. In the teleost fish medaka (Oryzias latipes), a germ cell-intrinsic sex determinant, foxl3, triggers germline feminization by activating two genetic pathways that regulate folliculogenesis and meiosis. Here, we identified a pathway involving a dome-shaped microtubule structure that may be the basis of oocyte polarity. This structure was first established in primordial germ cells in both sexes, but was maintained only during oogenesis and was destabilized in differentiating spermatogonia under the influence of Sertoli cells expressing dmrt1. Although foxl3 was dispensable for this pathway, dazl was involved in the persistence of the microtubule dome at the time of gonocyte development. In addition, disruption of the microtubule dome caused dispersal of bucky ball RNA, suggesting the structure may be prerequisite for the Balbiani body. Collectively, the present findings provide mechanistic insight into the establishment of sex-specific polarity through the formation of a microtubule structure in germ cells, as well as clarifying the genetic pathways implementing oocyte-specific characteristics.

The Hox code responsible for the patterning of the anterior vertebrae in zebrafish.

The vertebral column is a characteristic structure of vertebrates. Genetic studies in mice have shown that Hox-mediated patterning plays a key role in specifying discrete anatomical regions of the vertebral column. Expression pattern analyses in several vertebrate embryos have provided correlative evidence that the anterior boundaries of Hox expression coincide with distinct anatomical vertebrae. However, because functional analyses have been limited to mice, it remains unclear which Hox genes actually function in vertebral patterning in other vertebrates. In this study, various zebrafish Hox mutants were generated for loss-of-function phenotypic analysis to functionally decipher the Hox code responsible for the zebrafish anterior vertebrae between the occipital and thoracic vertebrae. We found that Hox genes in HoxB- and HoxC-related clusters participate in regulating the morphology of the zebrafish anterior vertebrae. In addition, medaka hoxc6a was found to be responsible for anterior vertebral identity, as in zebrafish. Based on phenotypic similarities with Hoxc6 knockout mice, our results suggest that the Hox patterning system, including at least Hoxc6, may have been functionally established in the vertebral patterning of the common ancestor of ray-finned and lobe-finned fishes.

The decision of male medaka to mate or fight depends on two complementary androgen signaling pathways.

Adult male animals typically court and attempt to mate with females, while attacking other males. Emerging evidence from mice indicates that neurons expressing the estrogen receptor ESR1 in behaviorally relevant brain regions play a central role in mediating these mutually exclusive behavioral responses to conspecifics. However, the findings in mice are unlikely to apply to vertebrates in general because, in many species other than rodents and some birds, androgens-rather than estrogens-have been implicated in male behaviors. Here, we report that male medaka (Oryzias latipes) lacking one of the two androgen receptor subtypes (Ara) are less aggressive toward other males and instead actively court them, while those lacking the other subtype (Arb) are less motivated to mate with females and conversely attack them. These findings indicate that, in male medaka, the Ara- and Arb-mediated androgen signaling pathways facilitate appropriate behavioral responses, while simultaneously suppressing inappropriate responses, to males and females, respectively. Notably, males lacking either receptor retain the ability to discriminate the sex of conspecifics, suggesting a defect in the subsequent decision-making process to mate or fight. We further show that Ara and Arb are expressed in intermingled but largely distinct populations of neurons, and stimulate the expression of different behaviorally relevant genes including galanin and vasotocin, respectively. Collectively, our results demonstrate that male teleosts make adaptive decisions to mate or fight as a result of the activation of one of two complementary androgen signaling pathways, depending on the sex of the conspecific that they encounter.

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins.

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

Role of trafficking protein particle complex 2 in medaka development.

The skeletal dysplasia spondyloepiphyseal dysplasia tarda (SEDT) is caused by mutations in the TRAPPC2 gene, which encodes Sedlin, a component of the trafficking protein particle (TRAPP) complex that we have shown previously to be required for the export of type II collagen (Col2) from the endoplasmic reticulum. No vertebrate model for SEDT has been generated thus far. To address this gap, we generated a Sedlin knockout animal by mutating the orthologous TRAPPC2 gene (olSedl) of Oryzias latipes (medaka) fish. OlSedl deficiency leads to embryonic defects, short size, diminished skeletal ossification and altered Col2 production and secretion, resembling human defects observed in SEDT patients. Moreover, SEDT knock-out animals display photoreceptor degeneration and gut morphogenesis defects, suggesting a key role for Sedlin in the development of these organs. Thus, by studying Sedlin function in vivo, we provide evidence for a mechanistic link between TRAPPC2-mediated membrane trafficking, Col2 export, and developmental disorders.

Incomplete erasure of histone marks during epigenetic reprogramming in medaka early development.

Epigenetic modifications undergo drastic erasure and reestablishment after fertilization. This reprogramming is required for proper embryonic development and cell differentiation. In mammals, some histone modifications are not completely reprogrammed and play critical roles in later development. In contrast, in nonmammalian vertebrates, most histone modifications are thought to be more intensively erased and reestablished by the stage of zygotic genome activation (ZGA). However, histone modifications that escape reprogramming in nonmammalian vertebrates and their potential functional roles remain unknown. Here, we quantitatively and comprehensively analyzed histone modification dynamics during epigenetic reprogramming in Japanese killifish, medaka (Oryzias latipes) embryos. Our data revealed that H3K27ac, H3K27me3, and H3K9me3 escape complete reprogramming, whereas H3K4 methylation is completely erased during cleavage stage. Furthermore, we experimentally showed the functional roles of such retained modifications at early stages: (i) H3K27ac premarks promoters during the cleavage stage, and inhibition of histone acetyltransferases disrupts proper patterning of H3K4 and H3K27 methylation at CpG-dense promoters, but does not affect chromatin accessibility after ZGA; (ii) H3K9me3 is globally erased but specifically retained at telomeric regions, which is required for maintenance of genomic stability during the cleavage stage. These results expand the understanding of diversity and conservation of reprogramming in vertebrates, and unveil previously uncharacterized functions of histone modifications retained during epigenetic reprogramming.

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish.

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

De novo PAM generation to reach initially inaccessible target sites for base editing.

Base editing by CRISPR crucially depends on the presence of a protospacer adjacent motif (PAM) at the correct distance from the editing site. Here, we present and validate an efficient one-shot approach termed 'inception' that expands the editing range. This is achieved by sequential, combinatorial base editing: de novo generated synonymous, non-synonymous or intronic PAM sites facilitate subsequent base editing at nucleotide positions that were initially inaccessible, further opening the targeting range of highly precise editing approaches. We demonstrate the applicability of the inception concept in medaka (Oryzias latipes) in three settings: loss of function, by introducing a pre-termination STOP codon in the open reading frame of oca2; locally confined multi-codon changes to generate allelic variants with different phenotypic severity in kcnh6a; and the removal of a splice acceptor site by targeting intronic sequences of rx3. Using sequentially acting base editors in the described combinatorial approach expands the number of accessible target sites by 65% on average. This allows the use of well-established tools with NGG PAM recognition for the establishment of thus far unreachable disease models, for hypomorphic allele studies and for efficient targeted mechanistic investigations in a precise and predictable manner.

Genetic and functional insights into the fractal structure of the heart.

The inner surfaces of the human heart are covered by a complex network of muscular strands that is thought to be a remnant of embryonic development1,2. The function of these trabeculae in adults and their genetic architecture are unknown. Here we performed a genome-wide association study to investigate image-derived phenotypes of trabeculae using the fractal analysis of trabecular morphology in 18,096 participants of the UK Biobank. We identified 16 significant loci that contain genes associated with haemodynamic phenotypes and regulation of cytoskeletal arborization3,4. Using biomechanical simulations and observational data from human participants, we demonstrate that trabecular morphology is an important determinant of cardiac performance. Through genetic association studies with cardiac disease phenotypes and Mendelian randomization, we find a causal relationship between trabecular morphology and risk of cardiovascular disease. These findings suggest a previously unknown role for myocardial trabeculae in the function of the adult heart, identify conserved pathways that regulate structural complexity and reveal the influence of the myocardial trabeculae on susceptibility to cardiovascular disease.

Looking for a needle in a haystack: de novo phenotypic target identification reveals Hippo pathway-mediated miR-202 regulation of egg production.

Understanding microRNA (miRNA) functions has been hampered by major difficulties in identifying their biological target(s). Currently, the main limitation is the lack of a suitable strategy to identify biologically relevant targets among a high number of putative targets. Here we provide a proof of concept of successful de novo (i.e. without prior knowledge of its identity) miRNA phenotypic target (i.e. target whose de-repression contributes to the phenotypic outcomes) identification from RNA-seq data. Using the medaka mir-202 knock-out (KO) model in which inactivation leads to a major organism-level reproductive phenotype, including reduced egg production, we introduced novel criteria including limited fold-change in KO and low interindividual variability in gene expression to reduce the list of 2853 putative targets to a short list of 5. We selected tead3b, a member of the evolutionarily-conserved Hippo pathway, known to regulate ovarian functions, due to its remarkably strong and evolutionarily conserved binding affinity for miR-202-5p. Deleting the miR-202-5p binding site in the 3' UTR of tead3b, but not of other Hippo pathway members sav1 and vgll4b, triggered a reduced egg production phenotype. This is one of the few successful examples of de novo functional assignment of a miRNA phenotypic target in vivo in vertebrates.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model.

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

A transcriptional program underlying the circannual rhythms of gonadal development in medaka.

To cope with seasonal environmental changes, organisms have evolved approximately 1-y endogenous circannual clocks. These circannual clocks regulate various physiological properties and behaviors such as reproduction, hibernation, migration, and molting, thus providing organisms with adaptive advantages. Although several hypotheses have been proposed, the genes that regulate circannual rhythms and the underlying mechanisms controlling long-term circannual clocks remain unknown in any organism. Here, we show a transcriptional program underlying the circannual clock in medaka fish (Oryzias latipes). We monitored the seasonal reproductive rhythms of medaka kept under natural outdoor conditions for 2 y. Linear regression analysis suggested that seasonal changes in reproductive activity were predominantly determined by an endogenous program. Medaka hypothalamic and pituitary transcriptomes were obtained monthly over 2 y and daily on all equinoxes and solstices. Analysis identified 3,341 seasonally oscillating genes and 1,381 daily oscillating genes. We then examined the existence of circannual rhythms in medaka via maintaining them under constant photoperiodic conditions. Medaka exhibited approximately 6-mo free-running circannual rhythms under constant conditions, and monthly transcriptomes under constant conditions identified 518 circannual genes. Gene ontology analysis of circannual genes highlighted the enrichment of genes related to cell proliferation and differentiation. Altogether, our findings support the "histogenesis hypothesis" that postulates the involvement of tissue remodeling in circannual time-keeping.

The TBC1D31/praja2 complex controls primary ciliogenesis through PKA-directed OFD1 ubiquitylation.

The primary cilium is a microtubule-based sensory organelle that dynamically links signalling pathways to cell differentiation, growth, and development. Genetic defects of primary cilia are responsible for genetic disorders known as ciliopathies. Orofacial digital type I syndrome (OFDI) is an X-linked congenital ciliopathy caused by mutations in the OFD1 gene and characterized by malformations of the face, oral cavity, digits and, in the majority of cases, polycystic kidney disease. OFD1 plays a key role in cilium biogenesis. However, the impact of signalling pathways and the role of the ubiquitin-proteasome system (UPS) in the control of OFD1 stability remain unknown. Here, we identify a novel complex assembled at centrosomes by TBC1D31, including the E3 ubiquitin ligase praja2, protein kinase A (PKA), and OFD1. We show that TBC1D31 is essential for ciliogenesis. Mechanistically, upon G-protein-coupled receptor (GPCR)-cAMP stimulation, PKA phosphorylates OFD1 at ser735, thus promoting OFD1 proteolysis through the praja2-UPS circuitry. This pathway is essential for ciliogenesis. In addition, a non-phosphorylatable OFD1 mutant dramatically affects cilium morphology and dynamics. Consistent with a role of the TBC1D31/praja2/OFD1 axis in ciliogenesis, alteration of this molecular network impairs ciliogenesis in vivo in Medaka fish, resulting in developmental defects. Our findings reveal a multifunctional transduction unit at the centrosome that links GPCR signalling to ubiquitylation and proteolysis of the ciliopathy protein OFD1, with important implications on cilium biology and development. Derangement of this control mechanism may underpin human genetic disorders.

A non-disruptive and efficient knock-in approach allows fate tracing of resident osteoblast progenitors during repair of vertebral lesions in medaka.

During bone development and repair, osteoblasts are recruited to bone deposition sites. To identify the origin of recruited osteoblasts, cell lineage tracing using Cre/loxP recombination is commonly used. However, a confounding factor is the use of transgenic Cre drivers that do not accurately recapitulate endogenous gene expression or the use of knock-in Cre drivers that alter endogenous protein activity or levels. Here, we describe a CRISPR/Cas9 homology-directed repair knock-in approach that allows efficient generation of Cre drivers controlled by the endogenous gene promoter. In addition, a self-cleaving peptide preserves the reading frame of the endogenous protein. Using this approach, we generated col10a1p2a-CreERT2 knock-in medaka and show that tamoxifen-inducible CreERT2 efficiently recombined loxP sites in col10a1 cells. Similar knock-in efficiencies were obtained when two unrelated loci (osr1 and col2a1a) were targeted. Using live imaging, we traced the fate of col10a1 osteoblast progenitors during bone lesion repair in the medaka vertebral column. We show that col10a1 cells at neural arches represent a mobilizable cellular source for bone repair. Together, our study describes a previously unreported strategy for precise cell lineage tracing via efficient and non-disruptive knock-in of Cre.

Neuropeptide FF indirectly affects testicular morphogenesis and functions in medaka.

Testicular morphogenesis and functions are considered to be under the control of neural and endocrine systems. However, the available literature is mainly limited to mammals, and it remains unclear how they are regulated in teleost species. Here, we demonstrated that neuropeptide FF (NPFF) in the brain is responsible for the follicle-stimulating hormone expression in the pituitary, which facilitates the testicular morphogenesis and androgen synthesis, and subsequently contributes to successful spermatogenesis. The present findings give us important insights into the neuroendocrine regulatory mechanisms underlying the testicular morphogenesis and functions in teleosts.

Spatio-temporal control of targeted gene expression in combination with CRISPR/Cas and Tet-On systems in Medaka.

Spatial and temporal control of transgene expression is a powerful approach to understand gene functions in specific cells and tissues. The Tet-On system is a robust tool for controlling transgene expression spatially and temporally; however, few studies have examined whether this system can be applied to postembryonic stages of Medaka (Oryzias latipes) or other fishes. Here, we first improved a basal promoter sequence on the donor vector for a nonhomologous end joining (NHEJ)-based knock-in (KI) system. Next, using transgenic Medaka for establishing the Tet-On system by KI, we demonstrated that doxycycline administration for four or more days by feeding can be a stable and efficient method to achieve expression of the transduced reporter gene in adult fish. From these analyses, we propose an optimized approach for a spatio-temporal gene-expression system in the adult stage of Medaka and other small fishes.

Sarcomere level mechanics of the fast skeletal muscle of the medaka fish larva.

The medaka fish (Oryzias latipes) is a vertebrate model used in developmental biology and genetics. Here we explore its suitability as a model for investigating the molecular mechanisms of human myopathies caused by mutations in sarcomeric proteins. To this end, the relevant mechanical parameters of the intact skeletal muscle of wild-type medaka are determined using the transparent tail at larval stage 40. Tails were mounted at sarcomere length of 2.1 µm in a thermoregulated trough containing physiological solution. Tetanic contractions were elicited at physiological temperature (10°C-30°C) by electrical stimulation, and sarcomere length changes were recorded with nanometer-microsecond resolution during both isometric and isotonic contractions with a striation follower. The force output has been normalized for the actual fraction of the cross section of the tail occupied by the myofilament lattice, as established with transmission electron microscopy (TEM), and then for the actual density of myofilaments, as established with X-ray diffraction. Under these conditions, the mechanical performance of the contracting muscle of the wild-type larva can be defined at the level of the half-thick filament, where ∼300 myosin motors work in parallel as a collective motor, allowing a detailed comparison with the established performance of the skeletal muscle of different vertebrates. The results of this study point out that the medaka fish larva is a suitable model for the investigation of the genotype/phenotype correlations and therapeutic possibilities in skeletal muscle diseases caused by mutations in sarcomeric proteins.NEW & NOTEWORTHY The suitability of the medaka fish as a model for investigating the molecular mechanisms of human myopathies caused by mutations of sarcomeric proteins is tested by combining structural analysis and sarcomere-level mechanics of the skeletal muscle of the tail of medaka larva. The mechanical performance of the medaka muscle, scaled at the level of the myosin-containing thick filament, together with its reduced genome duplication makes this model unique for investigations of the genotype/phenotype correlations in human myopathies.

Egg envelope formation of medaka Oryzias latipes requires ZP proteins originating from both the liver and ovary.

Teleost oocytes are surrounded by a structure called chorion or egg envelopes, which is composed of zona pellucida (ZP) proteins. As a result of the gene duplication in teleost, the expression site of the zp genes, coding the major component protein of egg envelopes, changed from the ovary to the maternal liver. In Euteleostei, there are three liver-expressed zp genes, named choriogenin (chg) h, chg hm, and chg l, and the composition of the egg envelope is mostly made up of these Chgs. In addition, ovary-expressed zp genes are also conserved in the medaka genomes, and their proteins have also been found to be minor components of the egg envelopes. However, the specific role of liver-expressed versus ovary-expressed zp genes was unclear. In the present study, we showed that ovary-synthesized ZP proteins first form the base layer of the egg envelope and then Chgs polymerize inwardly to thicken the egg envelope. To analyze the effects of dysfunction of the chg gene, we generated some chg knockout medaka. All knockout females failed to produce normally fertilized eggs by the natural spawning. The egg envelopes lacking Chgs were significantly thinner, but layers formed by ZP proteins synthesized in the ovary were found in the thin egg envelope of knockout as well as wildtype eggs. These results suggest that the ovary-expressed zp gene is well conserved in all teleosts, including those species in which liver-derived ZP proteins are the major component, because it is essential for the initiation of egg envelope formation.

Diversification processes of teleost intron-less opsin genes.

Opsins are universal photosensitive proteins in animals. Vertebrates have a variety of opsin genes for visual and non-visual photoreceptions. Analysis of the gene structures shows that most opsin genes have introns in their coding regions. However, teleosts exceptionally have several intron-less opsin genes that are presumed to have been duplicated by an RNA-based gene duplication mechanism, retroduplication. Among these retrogenes, we focused on the Opn4 (melanopsin) gene responsible for non-image-forming photoreception. Many teleosts have five Opn4 genes including one intron-less gene, which is speculated to have been formed from a parental intron-containing gene in the Actinopterygii. In this study, to reveal the evolutionary history of Opn4 genes, we analyzed them in teleost (zebrafish and medaka) and non-teleost (bichir, sturgeon, and gar) fishes. Our synteny analysis suggests that the intron-less Opn4 gene emerged by retroduplication after the branching of the bichir lineage. In addition, our biochemical and histochemical analyses showed that, in the teleost lineage, the newly acquired intron-less Opn4 gene became abundantly used without substantial changes in the molecular properties of the Opn4 protein. This stepwise evolutionary model of Opn4 genes is quite similar to that of rhodopsin genes in the Actinopterygii. The unique acquisition of rhodopsin and Opn4 retrogenes would have contributed to the diversification of the opsin gene repertoires in the Actinopterygii and the adaptation of teleosts to various aquatic environments.