Genes are thought to be transmitted to offspring by random fertilization of a small number of oocytes with numerous spermatozoa. Here we analyzed the dynamics of male germline transmission by genetic marking and transplantation of spermatogonial stem cells (SSCs). We found that offspring deriving from a small number of specific SSCs appear within a limited time. Interestingly, the same SSC clones reappear later with an average functional lifespan of ∼124.4 days. Cyclic offspring production from SSCs was not caused by changes in SSC self-renewal activity because lineage-tracing analyses suggested that all SSCs actively proliferated. Selection appears to occur during the differentiating spermatogonia stage, when extensive apoptosis was observed. The pattern of germline transmission could be predicted using a mathematical model in which SSCs repeat cycles of transient spermatogenic burst and refractory periods. Thus, spermatogenesis is a regulated process whereby specific SSC clones are repeatedly recruited for fertilization with long-term cycles.

A small number of offspring are born from the numerous sperm generated from spermatogonial stem cells (SSCs). However, little is known regarding the rules and molecular mechanisms that govern germline transmission patterns. Here we report that the Trp53 tumor suppressor gene limits germline genetic diversity via Cdkn1a. Trp53-deficient SSCs outcompeted wild-type (WT) SSCs and produced significantly more progeny after co-transplantation into infertile mice. Lentivirus-mediated transgenerational lineage analysis showed that offspring bearing the same virus integration were repeatedly born in a non-random pattern from WT SSCs. However, SSCs lacking Trp53 or Cdkn1a sired transgenic offspring in random patterns with increased genetic diversity. Apoptosis of KIT+ differentiating germ cells was reduced in Trp53- or Cdkn1a-deficient mice. Reduced CDKN1A expression in Trp53-deficient spermatogonia suggested that Cdkn1a limits genetic diversity by supporting apoptosis of syncytial spermatogonial clones. Therefore, the TRP53-CDKN1A pathway regulates tumorigenesis and the germline transmission pattern.

Extracellular signal-regulated kinase (ERK) is a key effector of many growth signalling pathways. In this study, we visualise epidermal ERK activity in living mice using an ERK FRET biosensor. Under steady-state conditions, the epidermis occasionally revealed bursts of ERK activation patterns where ERK activity radially propagated from cell to cell. The frequency of this spatial propagation of radial ERK activity distribution (SPREAD) correlated with the rate of epidermal cell division. SPREADs and proliferation were stimulated by 12-O-tetradecanoylphorbol 13-acetate (TPA) in a manner dependent on EGF receptors and their cognate ligands. At the wounded skin, ERK activation propagated as trigger wave in parallel to the wound edge, suggesting that ERK activation propagation can be superimposed. Furthermore, by visualising the cell cycle, we found that SPREADs were associated with G2/M cell cycle progression. Our results provide new insights into how cell proliferation and transient ERK activity are synchronised in a living tissue.

A long-standing question in neurodevelopment is how neurons develop a single axon and multiple dendrites from common immature neurites. Long-range inhibitory signaling from the growing axon is hypothesized to prevent outgrowth of other immature neurites and to differentiate them into dendrites, but the existence and nature of this inhibitory signaling remains unknown. Here, we demonstrate that axonal growth triggered by neurotrophin-3 remotely inhibits neurite outgrowth through long-range Ca2+ waves, which are delivered from the growing axon to the cell body. These Ca2+ waves increase RhoA activity in the cell body through calcium/calmodulin-dependent protein kinase I. Optogenetic control of Rho-kinase combined with computational modeling reveals that active Rho-kinase diffuses to growing other immature neurites and inhibits their outgrowth. Mechanistically, calmodulin-dependent protein kinase I phosphorylates a RhoA-specific GEF, GEF-H1, whose phosphorylation enhances its GEF activity. Thus, our results reveal that long-range inhibitory signaling mediated by Ca2+ wave is responsible for neuronal polarization.Emerging evidence suggests that gut microbiota influences immune function in the brain and may play a role in neurological diseases. Here, the authors offer in vivo evidence from a Drosophila model that supports a role for gut microbiota in modulating the progression of Alzheimer's disease.

During somite segmentation, clock genes oscillate within the posterior presomitic mesoderm (PSM). The temporal information ties up with the posteriorly moving FGF gradient, leading to the formation of a presumptive somite within the PSM. We previously investigated Erk activity downstream of FGF signaling by collecting stained zebrafish embryos, and discovered that the steep gradient of Erk activity was generated in the PSM, and the Erk activity border regularly shifted in a stepwise manner. However, since these interpretations come from static analyses, we needed to firmly confirm them by applying an analysis that has higher spatiotemporal resolutions. Here we developed a live imaging system for Erk activity in zebrafish embryos, using a Förster resonance energy transfer (FRET)-based Erk biosensor. With this system, we firmly showed that Erk activity exhibits stepwise regression within the PSM. Although our static analyses could not detect the stepwise pattern of Erk activity in clock-deficient embryos, our system revealed that, in clock-deficient embryos, the stepwise regression of Erk activity occurs at an irregular timing, eventually leading to formation of irregularly-sized somites. Therefore, our system overcame the limitation of static analyses and revealed that clock-dependent spatiotemporal regulation of Erk is required for proper somitogenesis in zebrafish.

The biophysical framework of collective cell migration has been extensively investigated in recent years; however, it remains elusive how chemical inputs from neighboring cells are integrated to coordinate the collective movement. Here, we provide evidence that propagation waves of extracellular signal-related kinase (ERK) mitogen-activated protein kinase activation determine the direction of the collective cell migration. A wound-healing assay of Mardin-Darby canine kidney (MDCK) epithelial cells revealed two distinct types of ERK activation wave, a "tidal wave" from the wound, and a self-organized "spontaneous wave" in regions distant from the wound. In both cases, MDCK cells collectively migrated against the direction of the ERK activation wave. The inhibition of ERK activation propagation suppressed collective cell migration. An ERK activation wave spatiotemporally controlled actomyosin contraction and cell density. Furthermore, an optogenetic ERK activation wave reproduced the collective cell migration. These data provide new mechanistic insight into how cells sense the direction of collective cell migration.

Myosin light chain (MLC) phosphorylation plays important roles in various cellular functions such as cellular morphogenesis, motility, and smooth muscle contraction. MLC phosphorylation is determined by the balance between activities of Rho-associated kinase (Rho-kinase) and myosin phosphatase. An impaired balance between Rho-kinase and myosin phosphatase activities induces the abnormal sustained phosphorylation of MLC, which contributes to the pathogenesis of certain vascular diseases, such as vasospasm and hypertension. However, the dynamic principle of the system underlying the regulation of MLC phosphorylation remains to be clarified. Here, to elucidate this dynamic principle whereby Rho-kinase regulates MLC phosphorylation, we developed a mathematical model based on the behavior of thrombin-dependent MLC phosphorylation, which is regulated by the Rho-kinase signaling network. Through analyzing our mathematical model, we predict that MLC phosphorylation and myosin phosphatase activity exhibit bistability, and that a novel signaling pathway leading to the auto-activation of myosin phosphatase is required for the regulatory system of MLC phosphorylation. In addition, on the basis of experimental data, we propose that the auto-activation pathway of myosin phosphatase occurs in vivo. These results indicate that bistability of myosin phosphatase activity is responsible for the bistability of MLC phosphorylation, and the sustained phosphorylation of MLC is attributed to this feature of bistability.

We developed two new FRET imaging measures for intramolecular FRET biosensors, called linearly proportional (LP) and highly contrasting (HC) measures, which can be easily calculated by the fluorescence intensities of donor and acceptor as a ratio between their weighted sums. As an alternative to the conventional ratiometric measure, which non-linearly depends on the fraction of active molecule, we first developed the LP measure, which is linearly proportional to the fraction of active molecules. The LP measure inherently unmixes bleed-through signals and is robust against fluorescence noise. By extending the LP measure, we furthermore designed the HC measure, which provides highly contrasting images of the molecular activity, more than the ratiometric measure. In addition to their advantages, these measures are insensitive to the biosensor expression level, which is a fundamental property of the ratiometric measure. Using artificial data and FRET imaging data, we showed that the LP measure effectively represents the fraction of active molecules and that the HC measure improves visual interpretability by providing high contrast images of molecular activity. Therefore, the LP and HC measures allow us to gain more quantitative and qualitative insights from FRET imaging than the ratiometric measure.

During development, positional information directs cells to specific fates, leading them to differentiate with their own transcriptomes and express specific behaviors and functions. However, the mechanisms underlying these processes in a genome-wide view remain ambiguous, partly because the single-cell transcriptomic data of early developing embryos containing accurate spatial and lineage information are still lacking. Here, we report a single-cell transcriptome atlas of Drosophila gastrulae, divided into 77 transcriptomically distinct clusters. We find that the expression profiles of plasma-membrane-related genes, but not those of transcription-factor genes, represent each germ layer, supporting the nonequivalent contribution of each transcription-factor mRNA level to effector gene expression profiles at the transcriptome level. We also reconstruct the spatial expression patterns of all genes at the single-cell stripe level as the smallest unit. This atlas is an important resource for the genome-wide understanding of the mechanisms by which genes cooperatively orchestrate Drosophila gastrulation.

We propose a new computation-based approach for elucidating how signaling molecules are decoded in cell migration. In this approach, we performed FRET time-lapse imaging of Rac1 and Cdc42, members of Rho GTPases which are responsible for cell motility, and quantitatively identified the response functions that describe the conversion from the molecular activities to the morphological changes. Based on the identified response functions, we clarified the profiles of how the morphology spatiotemporally changes in response to local and transient activation of Rac1 and Cdc42, and found that Rac1 and Cdc42 activation triggers laterally propagating membrane protrusion. The response functions were also endowed with property of differentiator, which is beneficial for maintaining sensitivity under adaptation to the mean level of input. Using the response function, we could predict the morphological change from molecular activity, and its predictive performance provides a new quantitative measure of how much the Rho GTPases participate in the cell migration. Interestingly, we discovered distinct predictive performance of Rac1 and Cdc42 depending on the migration modes, indicating that Rac1 and Cdc42 contribute to persistent and random migration, respectively. Thus, our proposed predictive approach enabled us to uncover the hidden information processing rules of Rho GTPases in the cell migration.

Decoding spatial transcriptomes from single-cell RNA sequencing (scRNA-seq) data has become a fundamental technique for understanding multicellular systems; however, existing computational methods lack both accuracy and biological interpretability due to their model-free frameworks. Here, we introduce Perler, a model-based method to integrate scRNA-seq data with reference in situ hybridization (ISH) data. To calibrate differences between these datasets, we develop a biologically interpretable model that uses generative linear mapping based on a Gaussian mixture model using the Expectation-Maximization algorithm. Perler accurately predicts the spatial gene expression of Drosophila embryos, zebrafish embryos, mammalian liver, and mouse visual cortex from scRNA-seq data. Furthermore, the reconstructed transcriptomes do not over-fit the ISH data and preserved the timing information of the scRNA-seq data. These results demonstrate the generalizability of Perler for dataset integration, thereby providing a biologically interpretable framework for accurate reconstruction of spatial transcriptomes in any multicellular system.

Wnt signaling is required to maintain bipotent progenitors for neural and paraxial mesoderm cells, the neuromesodermal progenitor (NMP) cells that reside in the epiblast and tailbud. Since epiblast/tailbud cells receive Wnt ligands produced by one another, this exchange may average out the heterogeneity of Wnt signaling levels among these cells. Here, we examined this possibility by replacing endogenous Wnt3a with a receptor-fused form that activates signaling in producing cells, but not in neighboring cells. Mutant mouse embryos show a unique phenotype in which maintenance of many NMP cells is impaired, although some cells persist for long periods. The epiblast cell population of these embryos increases heterogeneity in Wnt signaling levels as embryogenesis progresses and are sensitive to retinoic acid, an endogenous antagonist of NMP maintenance. Thus, mutual intercellular exchange of Wnt ligands in the epiblast cell population reduces heterogeneity and achieves robustness to environmental stress.

The reproducibility of embryonic development is remarkable, although molecular processes are intrinsically stochastic at the single-cell level. How the multicellular system resists the inevitable noise to acquire developmental reproducibility constitutes a fundamental question in developmental biology. Toward this end, we focused on vertebrate somitogenesis as a representative system, because somites are repeatedly reproduced within a single embryo whereas such reproducibility is lost in segmentation clock gene-deficient embryos. However, the effect of noise on developmental reproducibility has not been fully investigated, because of the technical difficulty in manipulating the noise intensity in experiments. In this study, we developed a computational model of ERK-mediated somitogenesis, in which bistable ERK activity is regulated by an FGF gradient, cell-cell communication, and the segmentation clock, subject to the intrinsic noise. The model simulation generated our previous in vivo observation that the ERK activity was distributed in a step-like gradient in the presomitic mesoderm, and its boundary was posteriorly shifted by the clock in a stepwise manner, leading to regular somite formation. Here, we showed that this somite regularity was robustly maintained against the noise. Removing the clock from the model predicted that the stepwise shift of the ERK activity occurs at irregular timing with irregular distance owing to the noise, resulting in somite size variation. This model prediction was recently confirmed by live imaging of ERK activity in zebrafish embryos. Through theoretical analysis, we presented a mechanism by which the clock reduces the inherent somite irregularity observed in clock-deficient embryos. Therefore, this study indicates a novel role of the segmentation clock in noise-resistant developmental reproducibility.

The immune system discriminates between harmful and harmless antigens based on past experiences; however, the underlying mechanism is largely unknown. From the viewpoint of machine learning, the learning system predicts the observation and updates the prediction based on prediction error, a process known as "predictive coding." Here, we modeled the population dynamics of T cells by adopting the concept of predictive coding; conventional and regulatory T cells predict the antigen concentration and excessive immune response, respectively. Their prediction error signals, possibly via cytokines, induce their differentiation to memory T cells. Through numerical simulations, we found that the immune system identifies antigen risks depending on the concentration and input rapidness of the antigen. Further, our model reproduced history-dependent discrimination, as in allergy onset and subsequent therapy. Taken together, this study provided a novel framework to improve our understanding of how the immune system adaptively learns the risks of diverse antigens.

To obtain more of a particular uncertain reward, animals must learn to actively overcome the lack of reward and adjust behavior to obtain it again. The neural mechanisms underlying such coping with reward omission remain unclear. Here, we developed a task in rats to monitor active behavioral switch toward the next reward after no reward. We found that some dopamine neurons in the ventral tegmental area exhibited increased responses to unexpected reward omission and decreased responses to unexpected reward, following the opposite responses of the well-known dopamine neurons that signal reward prediction error (RPE). The dopamine increase reflected in the nucleus accumbens correlated with behavioral adjustment to actively overcome unexpected no reward. We propose that these responses signal error to actively cope with lack of expected reward. The dopamine error signal thus cooperates with the RPE signal, enabling adaptive and robust pursuit of uncertain reward to ultimately obtain more reward.