Revisiting the evolution of bow-tie architecture in signaling networks.

Bow-tie architecture is a layered network structure that has a narrow middle layer with multiple inputs and outputs. Such structures are widely seen in the molecular networks in cells, suggesting that a universal evolutionary mechanism underlies the emergence of bow-tie architecture. The previous theoretical studies have implemented evolutionary simulations of the feedforward network to satisfy a given input-output goal and proposed that the bow-tie architecture emerges when the ideal input-output relation is given as a rank-deficient matrix with mutations in network link intensities in a multiplicative manner. Here, we report that the bow-tie network inevitably appears when the link intensities representing molecular interactions are small at the initial condition of the evolutionary simulation, regardless of the rank of the goal matrix. Our dynamical system analysis clarifies the mechanisms underlying the emergence of the bow-tie structure. Further, we demonstrate that the increase in the input-output matrix reduces the width of the middle layer, resulting in the emergence of bow-tie architecture, even when evolution starts from large link intensities. Our data suggest that bow-tie architecture emerges as a side effect of evolution rather than as a result of evolutionary adaptation.

Cell deformability drives fluid-to-fluid phase transition in active cell monolayers.

Cell deformability is an essential determinant for tissue-scale mechanical nature, such as fluidity and rigidity, and is thus crucial for tissue homeostasis and stable developmental processes. However, large-scale simulations of deformable cells have been restricted to those of polygonal-shaped cells, limiting our understanding of populations of arbitrarily deformable cells, such as mesenchymal, amoeboid cells, and nonconfluent epithelial cells. Here, we present an efficient approach for simulating large populations of nonpolygonally deformable cells with considerably higher computational efficiency than existing methods. Using the method, we demonstrate that the densely packed active cell population interacting via excluded volume interactions exhibits a fluid-to-fluid transition. An experimentally measurable index of topological defects, defined using the number of neighboring cells, is also proposed to characterize this transition. This study provides a flexible approach to tissue-scale cell population and a broader perspective on the biological fluid phases.

Evolution of hierarchy and irreversibility in theoretical cell differentiation model.

The process of cell differentiation in multicellular organisms is characterized by hierarchy and irreversibility in many cases. However, the conditions and selection pressures that give rise to these characteristics remain poorly understood. By using a mathematical model, here we show that the network of differentiation potency (differentiation diagram) becomes necessarily hierarchical and irreversible by increasing the number of terminally differentiated states under certain conditions. The mechanisms generating these characteristics are clarified using geometry in the cell state space. The results demonstrate that the hierarchical organization and irreversibility can manifest independently of direct selection pressures associated with these characteristics, instead they appear to evolve as byproducts of selective forces favoring a diversity of differentiated cell types. The study also provides a new perspective on the structure of gene regulatory networks that produce hierarchical and irreversible differentiation diagrams. These results indicate some constraints on cell differentiation, which are expected to provide a starting point for theoretical discussion of the implicit limits and directions of evolution in multicellular organisms.

A deep learning approach for morphological feature extraction based on variational auto-encoder: an application to mandible shape.

Shape measurements are crucial for evolutionary and developmental biology; however, they present difficulties in the objective and automatic quantification of arbitrary shapes. Conventional approaches are based on anatomically prominent landmarks, which require manual annotations by experts. Here, we develop a machine-learning approach by presenting morphological regulated variational AutoEncoder (Morpho-VAE), an image-based deep learning framework, to conduct landmark-free shape analysis. The proposed architecture combines the unsupervised and supervised learning models to reduce dimensionality by focusing on morphological features that distinguish data with different labels. We applied the method to primate mandible image data. The extracted morphological features reflected the characteristics of the families to which the organisms belonged, despite the absence of correlation between the extracted morphological features and phylogenetic distance. Furthermore, we demonstrated the reconstruction of missing segments from incomplete images. The proposed method provides a flexible and promising tool for analyzing a wide variety of image data of biological shapes even those with missing segments.

Intercellular exchange of Wnt ligands reduces cell population heterogeneity during embryogenesis.

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.

Rare-event sampling analysis uncovers the fitness landscape of the genetic code.

The genetic code refers to a rule that maps 64 codons to 20 amino acids. Nearly all organisms, with few exceptions, share the same genetic code, the standard genetic code (SGC). While it remains unclear why this universal code has arisen and been maintained during evolution, it may have been preserved under selection pressure. Theoretical studies comparing the SGC and numerically created hypothetical random genetic codes have suggested that the SGC has been subject to strong selection pressure for being robust against translation errors. However, these prior studies have searched for random genetic codes in only a small subspace of the possible code space due to limitations in computation time. Thus, how the genetic code has evolved, and the characteristics of the genetic code fitness landscape, remain unclear. By applying multicanonical Monte Carlo, an efficient rare-event sampling method, we efficiently sampled random codes from a much broader random ensemble of genetic codes than in previous studies, estimating that only one out of every 1020 random codes is more robust than the SGC. This estimate is significantly smaller than the previous estimate, one in a million. We also characterized the fitness landscape of the genetic code that has four major fitness peaks, one of which includes the SGC. Furthermore, genetic algorithm analysis revealed that evolution under such a multi-peaked fitness landscape could be strongly biased toward a narrow peak, in an evolutionary path-dependent manner.

Microtopographical guidance of macropinocytic signaling patches.

In fast-moving cells such as amoeba and immune cells, dendritic actin filaments are spatiotemporally regulated to shape large-scale plasma membrane protrusions. Despite their importance in migration, as well as in particle and liquid ingestion, how their dynamics are affected by micrometer-scale features of the contact surface is still poorly understood. Here, through quantitative image analysis of Dictyostelium on microfabricated surfaces, we show that there is a distinct mode of topographical guidance directed by the macropinocytic membrane cup. Unlike other topographical guidance known to date that depends on nanometer-scale curvature sensing protein or stress fibers, the macropinocytic membrane cup is driven by the Ras/PI3K/F-actin signaling patch and its dependency on the micrometer-scale topographical features, namely PI3K/F-actin-independent accumulation of Ras-GTP at the convex curved surface, PI3K-dependent patch propagation along the convex edge, and its actomyosin-dependent constriction at the concave edge. Mathematical model simulations demonstrate that the topographically dependent initiation, in combination with the mutually defining patch patterning and the membrane deformation, gives rise to the topographical guidance. Our results suggest that the macropinocytic cup is a self-enclosing structure that can support liquid ingestion by default; however, in the presence of structured surfaces, it is directed to faithfully trace bent and bifurcating ridges for particle ingestion and cell guidance.

Three-dimensional morphodynamic simulations of macropinocytic cups.

Macropinocytosis refers to the non-specific uptake of extracellular fluid, which plays ubiquitous roles in cell growth, immune surveillance, and virus entry. Despite its widespread occurrence, it remains unclear how its initial cup-shaped plasma membrane extensions form without any external solid support, as opposed to the process of particle uptake during phagocytosis. Here, by developing a computational framework that describes the coupling between the bistable reaction-diffusion processes of active signaling patches and membrane deformation, we demonstrated that the protrusive force localized to the edge of the patches can give rise to a self-enclosing cup structure, without further assumptions of local bending or contraction. Efficient uptake requires a balance among the patch size, magnitude of protrusive force, and cortical tension. Furthermore, our model exhibits cyclic cup formation, coexistence of multiple cups, and cup-splitting, indicating that these complex morphologies self-organize via a common mutually-dependent process of reaction-diffusion and membrane deformation.

Comparative mapping of crawling-cell morphodynamics in deep learning-based feature space.

Navigation of fast migrating cells such as amoeba Dictyostelium and immune cells are tightly associated with their morphologies that range from steady polarized forms that support high directionality to those more complex and variable when making frequent turns. Model simulations are essential for quantitative understanding of these features and their origins, however systematic comparisons with real data are underdeveloped. Here, by employing deep-learning-based feature extraction combined with phase-field modeling framework, we show that a low dimensional feature space for 2D migrating cell morphologies obtained from the shape stereotype of keratocytes, Dictyostelium and neutrophils can be fully mapped by an interlinked signaling network of cell-polarization and protrusion dynamics. Our analysis links the data-driven shape analysis to the underlying causalities by identifying key parameters critical for migratory morphologies both normal and aberrant under genetic and pharmacological perturbations. The results underscore the importance of deciphering self-organizing states and their interplay when characterizing morphological phenotypes.

Adaptation of metabolite leakiness leads to symbiotic chemical exchange and to a resilient microbial ecosystem.

Microbial communities display remarkable diversity, facilitated by the secretion of chemicals that can create new niches. However, it is unclear why cells often secrete even essential metabolites after evolution. Based on theoretical results indicating that cells can enhance their own growth rate by leaking even essential metabolites, we show that such "leaker" cells can establish an asymmetric form of mutualism with "consumer" cells that consume the leaked chemicals: the consumer cells benefit from the uptake of the secreted metabolites, while the leaker cells also benefit from such consumption, as it reduces the metabolite accumulation in the environment and thereby enables further secretion, resulting in frequency-dependent coexistence of multiple microbial species. As supported by extensive simulations, such symbiotic relationships generally evolve when each species has a complex reaction network and adapts its leakiness to optimize its own growth rate under crowded conditions and nutrient limitations. Accordingly, symbiotic ecosystems with diverse cell species that leak and exchange many metabolites with each other are shaped by cell-level adaptation of leakiness of metabolites. Moreover, the resultant ecosystems with entangled metabolite exchange are resilient against structural and environmental perturbations. Thus, we present a theory for the origin of resilient ecosystems with diverse microbes mediated by secretion and exchange of essential chemicals.

Advantage of Leakage of Essential Metabolites for Cells.

Microbial cells generally leak various metabolites including those necessary to grow. Why cells secrete even essential chemicals so often is, however, still unclear. Based on analytical and numerical calculations, we show that if the intracellular metabolism includes multibody (e.g., catalytic) reactions, leakage of essential metabolites can promote the leaking cell's growth. This advantage is typical for most metabolic networks via "flux control" and "growth-dilution" mechanisms, as a general consequence of the balance between synthesis and growth-induced dilution with autocatalytic reactions. We further argue that this advantage may lead to a novel form of symbiosis among diverse cells.

Embedding dual function into molecular motors through collective motion.

Protein motors, such as kinesins and dyneins, bind to a microtubule and travel along it in a specific direction. Previously, it was thought that the directionality for a given motor was constant in the absence of an external force. However, the directionality of the kinesin-5 Cin8 was recently found to change as the number of motors that bind to the same microtubule is increased. Here, we introduce a simple mechanical model of a microtubule-sliding assay in which multiple motors interact with the filament. We show that, due to the collective phenomenon, the directionality of the motor changes (e.g., from minus- to plus- end directionality), depending on the number of motors. This is induced by a large diffusive component in the directional walk and by the subsequent frustrated motor configuration, in which multiple motors pull the filament in opposite directions, similar to a game of tug-of-war. A possible role of the dual-directional motors for the mitotic spindle formation is also discussed. Our framework provides a general mechanism to embed two conflicting tasks into a single molecular machine, which works context-dependently.

Symbiotic Cell Differentiation and Cooperative Growth in Multicellular Aggregates.

As cells grow and divide under a given environment, they become crowded and resources are limited, as seen in bacterial biofilms and multicellular aggregates. These cells often show strong interactions through exchanging chemicals, as evident in quorum sensing, to achieve mutualism and division of labor. Here, to achieve stable division of labor, three characteristics are required. First, isogenous cells differentiate into several types. Second, this aggregate of distinct cell types shows better growth than that of isolated cells without interaction and differentiation, by achieving division of labor. Third, this cell aggregate is robust with respect to the number distribution of differentiated cell types. Indeed, theoretical studies have thus far considered how such cooperation is achieved when the ability of cell differentiation is presumed. Here, we address how cells acquire the ability of cell differentiation and division of labor simultaneously, which is also connected with the robustness of a cell society. For this purpose, we developed a dynamical-systems model of cells consisting of chemical components with intracellular catalytic reaction dynamics. The reactions convert external nutrients into internal components for cellular growth, and the divided cells interact through chemical diffusion. We found that cells sharing an identical catalytic network spontaneously differentiate via induction from cell-cell interactions, and then achieve division of labor, enabling a higher growth rate than that in the unicellular case. This symbiotic differentiation emerged for a class of reaction networks under the condition of nutrient limitation and strong cell-cell interactions. Then, robustness in the cell type distribution was achieved, while instability of collective growth could emerge even among the cooperative cells when the internal reserves of products were dominant. The present mechanism is simple and general as a natural consequence of interacting cells with limited resources, and is consistent with the observed behaviors and forms of several aggregates of unicellular organisms.

Motif analysis for small-number effects in chemical reaction dynamics.

The number of molecules involved in a cell or subcellular structure is sometimes rather small. In this situation, ordinary macroscopic-level fluctuations can be overwhelmed by non-negligible large fluctuations, which results in drastic changes in chemical-reaction dynamics and statistics compared to those observed under a macroscopic system (i.e., with a large number of molecules). In order to understand how salient changes emerge from fluctuations in molecular number, we here quantitatively define small-number effect by focusing on a "mesoscopic" level, in which the concentration distribution is distinguishable both from micro- and macroscopic ones and propose a criterion for determining whether or not such an effect can emerge in a given chemical reaction network. Using the proposed criterion, we systematically derive a list of motifs of chemical reaction networks that can show small-number effects, which includes motifs showing emergence of the power law and the bimodal distribution observable in a mesoscopic regime with respect to molecule number. The list of motifs provided herein is helpful in the search for candidates of biochemical reactions with a small-number effect for possible biological functions, as well as for designing a reaction system whose behavior can change drastically depending on molecule number, rather than concentration.

Theoretical analysis of discreteness-induced transition in autocatalytic reaction dynamics.

Transitions in the qualitative behavior of chemical reaction dynamics with a decrease in molecule number have attracted much attention. Here, a method based on a Markov process with a tridiagonal transition matrix is applied to the analysis of this transition in reaction dynamics. The transition to bistability due to the small-number effect and the mean switching time between the bistable states are analytically calculated in agreement with numerical simulations. In addition, a novel transition involving the reversal of the chemical reaction flow is found in the model under an external flow, and also in a three-component model. The generality of this transition and its correspondence to biological phenomena are also discussed.

Baldwin effect under multipeaked fitness landscapes: phenotypic fluctuation accelerates evolutionary rate.

Phenotypic fluctuations and plasticity can generally affect the course of evolution, a process known as the Baldwin effect. Several studies have recast this effect and claimed that phenotypic plasticity accelerates evolutionary rate (the Baldwin expediting effect); however, the validity of this claim is still controversial. In this study, we investigate the evolutionary population dynamics of a quantitative genetic model under a multipeaked fitness landscape, in order to evaluate the validity of the effect. We provide analytical expressions for the evolutionary rate and average population fitness. Our results indicate that under a multipeaked fitness landscape, phenotypic fluctuation always accelerates evolutionary rate, but it decreases the average fitness. As an extreme case of the trade-off between the rate of evolution and average fitness, phenotypic fluctuation is shown to accelerate the error catastrophe, in which a population fails to sustain a high-fitness peak. In the context of our findings, we discuss the role of phenotypic plasticity in adaptive evolution.

Multicanonical sampling of rare events in random matrices.

A method based on multicanonical Monte Carlo is applied to the calculation of large deviations in the largest eigenvalue of random matrices. The method is successfully tested with the Gaussian orthogonal ensemble, sparse random matrices, and matrices whose components are subject to uniform density. Specifically, the probability that all eigenvalues of a matrix are negative is estimated in these cases down to the values of ∼10(-200), a region where simple random sampling is ineffective. The method can be applied to any ensemble of matrices and used for sampling rare events characterized by any statistics.