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.

Modular control of vertebrate axis segmentation in time and space.

How the timing of development is linked to organismal size is a longstanding question. Although numerous studies have reported a correlation of temporal and spatial traits, the developmental or selective constraints underlying this link remain largely unexplored. We address this question by studying the periodic process of embryonic axis segmentation in-vivo in Oryzias fish. Interspecies comparisons reveal that the timing of segmentation correlates to segment, tissue and organismal size. Segment size in turn scales according to tissue and organism size. To probe for underlying causes, we genetically hybridised two closely related species. Quantitative analysis in ~600 phenotypically diverse F2 embryos reveals a decoupling of timing from size control, while spatial scaling is preserved. Using developmental quantitative trait loci (devQTL) mapping we identify distinct genetic loci linked to either the control of segmentation timing or tissue size. This study demonstrates that a developmental constraint mechanism underlies spatial scaling of axis segmentation, while its spatial and temporal control are dissociable modules.

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.

Cell cycle length governs heterochromatin reprogramming during early development in non-mammalian vertebrates.

Heterochromatin marks such as H3K9me3 undergo global erasure and re-establishment after fertilization, and the proper reprogramming of H3K9me3 is essential for early development. Despite the widely conserved dynamics of heterochromatin reprogramming in invertebrates and non-mammalian vertebrates, previous studies have shown that the underlying mechanisms may differ between species. Here, we investigate the molecular mechanism of H3K9me3 dynamics in medaka (Japanese killifish, Oryzias latipes) as a non-mammalian vertebrate model, and show that rapid cell cycle during cleavage stages causes DNA replication-dependent passive erasure of H3K9me3. We also find that cell cycle slowing, toward the mid-blastula transition, permits increasing nuclear accumulation of H3K9me3 histone methyltransferase Setdb1, leading to the onset of H3K9me3 re-accumulation. We further demonstrate that cell cycle length in early development also governs H3K9me3 reprogramming in zebrafish and Xenopus laevis. Together with the previous studies in invertebrates, we propose that a cell cycle length-dependent mechanism for both global erasure and re-accumulation of H3K9me3 is conserved among rapid-cleavage species of non-mammalian vertebrates and invertebrates such as Drosophila, C. elegans, Xenopus and teleost fish.

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.