Koichi Fujimoto

The Onset of Collective Behavior in Social Amoebae

All Together Now In the social amoeba Dictyostelium discoideum, periodic synthesis and release of cyclic AMP (cAMP) guides the cellular aggregation required to form fruiting bodies. It has been unclear whether the initiation of this behavior is owing to synchronization of autonomously oscillating cells or whether individual cells remain nonoscillatory unless the entire population becomes oscillatory. Gregor et al. (p. 1021, published online 22 April; see the Perspective by Prindle and Hasty) used live-cell imaging to show that cAMP pulses originate from a specific location in space and that individual cells move in and out of these signaling centers. The observations suggest that oscillations do not originate from autonomous activities of specialized cells. However, individual cells do display stochastic cAMP-pulsing below a threshold external concentration of cAMP, and the generation of synchronized oscillations could only be modeled accurately when this random pulsing was taken into account. Stochastic pulsing of individual cells plays a critical role in initiating cyclic adenosine monophosphate pulses. In the social amoebae Dictyostelium discoideum, periodic synthesis and release of extracellular cyclic adenosine 3′,5′-monophosphate (cAMP) guide cell aggregation and commitment to form fruiting bodies. It is unclear whether these oscillations are an intrinsic property of individual cells or if they exist only as a population-level phenomenon. Here, we showed by live-cell imaging of intact cell populations that pulses originate from a discrete location despite constant exchange of cells to and from the region. In a perfusion chamber, both isolated single cells and cell populations switched from quiescence to rhythmic activity depending on the concentration of extracellular cAMP. A quantitative analysis showed that stochastic pulsing of individual cells below the threshold concentration of extracellular cAMP plays a critical role in the onset of collective behavior.

Cross talking of network motifs in gene regulation that generates temporal pulses and spatial stripes.

Gene regulatory networks contain several substructures called network motifs, which frequently exist throughout the networks. One of such motifs found in Escherichia coli, Saccharomyces cerevisiae, and Drosophila melanogaster is the feed‐forward loop, in which an effector regulates its target by a direct regulatory interaction and an indirect interaction mediated by another gene product. Here, we theoretically analyze the behavior of networks that contain feed‐forward loops cross talking to each other. In response to levels of the effecter, such networks can generate multiple rise‐and‐fall temporal expression profiles and spatial stripes, which are typically observed in developmental processes. The mechanism to generate these responses reveals the way of inferring the regulatory pathways from experimental results. Our database study of gene regulatory networks indicates that most feed‐forward loops actually cross talk. We discuss how the feed‐forward loops and their cross talks can play important roles in morphogenesis.

Network Evolution of Body Plans

One of the major goals in evolutionary developmental biology is to understand the relationship between gene regulatory networks and the diverse morphologies and their functionalities. Are the diversities solely triggered by random events, or are they inevitable outcomes of an interplay between evolving gene networks and natural selection? Segmentation in arthropod embryogenesis represents a well-known example of body plan diversity. Striped patterns of gene expression that lead to the future body segments appear simultaneously or sequentially in long and short germ-band development, respectively. Moreover, a combination of both is found in intermediate germ-band development. Regulatory genes relevant for stripe formation are evolutionarily conserved among arthropods, therefore the differences in the observed traits are thought to have originated from how the genes are wired. To reveal the basic differences in the network structure, we have numerically evolved hundreds of gene regulatory networks that produce striped patterns of gene expression. By analyzing the topologies of the generated networks, we show that the characteristics of stripe formation in long and short germ-band development are determined by Feed-Forward Loops (FFLs) and negative Feed-Back Loops (FBLs) respectively, and those of intermediate germ-band development are determined by the interconnections between FFL and negative FBL. Network architectures, gene expression patterns and knockout responses exhibited by the artificially evolved networks agree with those reported in the fly Drosophila melanogaster and the beetle Tribolium castaneum. For other arthropod species, principal network architectures that remain largely unknown are predicted. Our results suggest that the emergence of the three modes of body segmentation in arthropods is an inherent property of the evolving networks.

Inter‐origin cooperativity of geminin action establishes an all‐or‐none switch for replication origin licensing

In metazoans, geminin functions as a molecular switch for preventing re‐replication of chromosomal DNA. Geminin binds to and inhibits Cdt1, which is required for replication origin licensing, but little is known about the mechanisms underlying geminin’s all‐or‐none action in licensing inhibition. Using Xenopus egg extract, we found that the all‐or‐none activity correlated with the formation of Cdt1 foci on chromatin, suggesting that multiple Cdt1–geminin complexes on origins cooperatively inhibit licensing. Based on experimental identification of licensing intermediates targeted by geminin and Cdt1, we developed a mathematical model of the licensing process. The model involves positive feedback owing to the cooperative action of geminin at neighboring origins and accurately accounts for the licensing activity mediated by geminin and Cdt1 in the extracts. The model also predicts that such cooperativity leads to clustering of licensing‐inhibited origins, an idea that is supported by the experimentally measured distribution of inter‐origin distances. We propose that geminin inhibits licensing through an inter‐origin interaction, ensuring strict and coordinated control of multiple replication origins on chromosomes.

Collective oscillations in developing cells: Insights from simple systems

From hormonal secretion to gene expression, multicellular dynamics are rich in oscillatory regulation. When organized in space and time, periodic cell–cell signaling can give rise to long‐range coordination of gene expression and cell movement in tissues. Lack of synchrony of the oscillations on the other hand can serve as a source of initial divergence of cell fate in stem cells. How properties of individual cells can account for collective rhythmic behaviors at the tissue level remains elusive in most cases. Recently, studies in chemical reactions, synthetic gene circuits, yeast and social amoeba Dictyostelium have greatly enhanced our view of collective oscillations in cell populations. From these relatively simple systems, a unified view of how excitable and oscillatory regulations could be tuned and coupled to give rise to tissue‐level oscillations is emerging. The review focuses on recent progress in cyclic adenosine monophosphate oscillations in Dictyostelium and highlights similarities and differences with other systems. We will see that the autonomy of single‐cell level oscillations and different ways in which cells are coupled influence how group‐level information can be encoded in collective oscillations.

A Design Principle of Group-level Decision Making in Cell Populations

Populations of cells often switch states as a group to cope with environmental changes such as nutrient availability and cell density. Although the gene circuits that underlie the switches are well understood at the level of single cells, the ways in which such circuits work in concert among many cells to support group-level switches are not fully explored. Experimental studies of microbial quorum sensing show that group-level changes in cellular states occur in either a graded or an all-or-none fashion. Here, we show through numerical simulations and mathematical analysis that these behaviors generally originate from two distinct forms of bistability. The choice of bistability is uniquely determined by a dimensionless parameter that compares the synthesis and the transport of the inducing molecules. The role of the parameter is universal, such that it not only applies to the autoinducing circuits typically found in bacteria but also to the more complex gene circuits involved in transmembrane receptor signaling. Furthermore, in gene circuits with negative feedback, the same dimensionless parameter determines the coherence of group-level transitions from quiescence to a rhythmic state. The set of biochemical parameters in bacterial quorum-sensing circuits appear to be tuned so that the cells can use either type of transition. The design principle identified here serves as the basis for the analysis and control of cellular collective decision making.

A developmental basis for stochasticity in floral organ numbers

Stochasticity ubiquitously inevitably appears at all levels from molecular traits to multicellular, morphological traits. Intrinsic stochasticity in biochemical reactions underlies the typical intercellular distributions of chemical concentrations, e.g., morphogen gradients, which can give rise to stochastic morphogenesis. While the universal statistics and mechanisms underlying the stochasticity at the biochemical level have been widely analyzed, those at the morphological level have not. Such morphological stochasticity is found in foral organ numbers. Although the floral organ number is a hallmark of floral species, it can distribute stochastically even within an individual plant. The probability distribution of the floral organ number within a population is usually asymmetric, i.e., it is more likely to increase rather than decrease from the modal value, or vice versa. We combined field observations, statistical analysis, and mathematical modeling to study the developmental basis of the variation in floral organ numbers among 50 species mainly from Ranunculaceae and several other families from core eudicots. We compared six hypothetical mechanisms and found that a modified error function reproduced much of the asymmetric variation found in eudicot floral organ numbers. The error function is derived from mathematical modeling of floral organ positioning, and its parameters represent measurable distances in the floral bud morphologies. The model predicts two developmental sources of the organ-number distributions: stochastic shifts in the expression boundaries of homeotic genes and a semi-concentric (whorled-type) organ arrangement. Other models species- or organ-specifically reproduced different types of distributions that reflect different developmental processes. The organ-number variation could be an indicator of stochasticity in organ fate determination and organ positioning.

Bifurcation cascade as chaotic itinerancy with multiple time scales

A coupled chaotic system with a variety of time scales is studied. Under a certain condition, it is shown that a change in fast dynamics can influence slow dynamics with a huge time-scale difference, successively through propagation of correlation over elements. This propagation is given by bifurcation cascade, for which three conditions are found to be essential: strong correlation, bifurcation of fast element dynamics by the change of its parameter, and marginal stability. By using coupled Lorenz equation with multiple time scales, it is shown that chaotic itinerancy (CI) observed there leads for the system to satisfy the three conditions, and the bifurcation cascade occurs. Through the analysis of the CI, characteristic properties of the bifurcation cascade, asymmetry in propagation with regards to the time scale, and the universality of the results are discussed, with possible relevance to biological memory.

A dynamical phyllotaxis model to determine floral organ number.

How organisms determine particular organ numbers is a fundamental key to the development of precise body structures; however, the developmental mechanisms underlying organ-number determination are unclear. In many eudicot plants, the primordia of sepals and petals (the floral organs) first arise sequentially at the edge of a circular, undifferentiated region called the floral meristem, and later transition into a concentric arrangement called a whorl, which includes four or five organs. The properties controlling the transition to whorls comprising particular numbers of organs is little explored. We propose a development-based model of floral organ-number determination, improving upon earlier models of plant phyllotaxis that assumed two developmental processes: the sequential initiation of primordia in the least crowded space around the meristem and the constant growth of the tip of the stem. By introducing mutual repulsion among primordia into the growth process, we numerically and analytically show that the whorled arrangement emerges spontaneously from the sequential initiation of primordia. Moreover, by allowing the strength of the inhibition exerted by each primordium to decrease as the primordium ages, we show that pentamerous whorls, in which the angular and radial positions of the primordia are consistent with those observed in sepal and petal primordia in Silene coeli-rosa, Caryophyllaceae, become the dominant arrangement. The organ number within the outmost whorl, corresponding to the sepals, takes a value of four or five in a much wider parameter space than that in which it takes a value of six or seven. These results suggest that mutual repulsion among primordia during growth and a temporal decrease in the strength of the inhibition during initiation are required for the development of the tetramerous and pentamerous whorls common in eudicots.

Robustness of Self-Organizing Chemoattractant Field Arising from Precise Pulse Induction of Its Breakdown Enzyme: A Single-Cell Level Analysis of PDE Expression in Dictyostelium

The oscillation of chemoattractant cyclic AMP (cAMP) in Dictyostelium discoideum is a collective phenomenon that occurs when the basal level of extracellular cAMP exceeds a threshold and invokes cooperative mutual excitation of cAMP synthesis and secretion. For pulses to be relayed from cell to cell repetitively, secreted cAMP must be cleared and brought down to the subthreshold level. One of the main determinants of the oscillatory behavior is thus how much extracellular cAMP is degraded by extracellular phosphodiesterase (PDE). To date, the exact nature of PDE gene regulation remains elusive. Here, we performed live imaging analysis of mRNA transcripts for pdsA--the gene encoding extracellular PDE. Our analysis revealed that pdsA is upregulated during the rising phase of cAMP oscillations. Furthermore, by analyzing isolated cells, we show that expression of pdsA is strictly dependent on the presence of extracellular cAMP. pdsA is induced only at ∼1 nM extracellular cAMP, which is almost identical to the threshold concentration for the cAMP relay response. The observed precise regulation of PDE expression together with degradation of extracellular cAMP by PDE form a dual positive and negative feedback circuit, and model analysis shows that this sets the cAMP level near the threshold concentration for the cAMP relay response for a wide range of adenylyl cyclase activity. The overlap of the thresholds could allow oscillations of chemoattractant cAMP to self-organize at various starving conditions, making its development robust to fluctuations in its environment.