Partial penetrance facilitates developmental evolution in bacteria
Development normally occurs similarly in all individuals within an isogenic population, but mutations often affect the fates of individual organisms differently. This phenomenon, known as partial penetrance, has been observed in diverse developmental systems. However, it remains unclear how the underlying genetic network specifies the set of possible alternative fates and how the relative frequencies of these fates evolve. We examined this question in the context of a simple developmental network; Bacillus subtilis sporulation. Mutations in an intercompartmental signalling process generate a set of discrete alternative fates not observed in wild-type cells, including rare formation of two viable 'twin' spores, rather than one within a single cell. Using genetic perturbations and time-lapse microscopy we identify the genetic control over twin penetrance and suggest a mechanistic model for twin formation. The results suggest a potential pathway for developmental evolution between mono-sporulation and twin-sporulation through states of intermediate twin penetrance. Furthermore, time-lapse microscopy of twin sporulation in wild-type Clostridium oceanicum shows a strong resemblance to twin sporulation in these B. subtilis mutants. Together the results suggest that noise can facilitate developmental evolution by enabling the initial expression of discrete morphological traits at low penetrance, and allowing their stabilization by gradual adjustment of genetic parameters. A simple mathematical model capture some of the essential requirements for such an evolution to occur.Network Evolution of Morphogenesis.
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. 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.Biodiversity, Endosymbiosis and Evolution
In nature, organisms are living in close association with surrounding physical environment as well as other organisms. Thus, each of the organisms is regarded as a component of the ecosystem. Considering the diverse microbial community found inside the organisms in general, each of the organisms can also be regarded as constituting an ecosystem.A Dynamic Model for Irreversible Differentiation in a Stem Cell System: Chaos Hypothesis
During normal development, cells undergo a unidirectional course of differentiation that progressively decreases the number of cell types they can potentially become. Pluripotent cells can differentiate into any of the cell types, but terminally differentiated cells cannot differentiate any more. A fundamental problem in stem cell biology is how to characterize the difference in cellular states between undifferentiated stem cells and terminally differentiated cells. In this study, to answer the raised question, we have carried out a dynamical systems modeling of cells accounting for the loss of multipotency resulting from stem cell differentiation triggered by the cell-cell interactions. We used an abstract cell model in which intra-cellular states are controlled by gene regulatory network and cell-cell interaction. Using Genetic Algorithm (GA) simulations, we screened regulatory networks that enable to generate a cell population with heterogeneous cell types by the cell-cell interactions. After screening regulatory networks generating cell type heterogeneity, we found that intra-cellular gene expression dynamics of cells that have potential to differentiate (i.e. stem cells) generally show irregular (chaotic) oscillation, whereas such complex oscillation is lost for determined cell types. Also, such stem cells express larger number of genes in comparison with differentiated cells in which lineage-specific genes strongly expressed. Based on the simulation results we proposed following hypothesis: i) Expression dynamics in stem cells show irregular chaotic oscillation and such oscillation disappear in differentiated cells. ii) The oscillatory dynamics in stem cells involves changes of many genes, whose number (effective degrees of freedom in the sense of dynamical systems) decreases as the differentiation progresses. iii) Robustness in the population ratio of several cell types is maintained by regulated differentiation of stem cells with irregular oscillatory dynamics. The universal features in regulatory networks generating cell type heterogeneity and experimental results supporting this hypothesis will also be discussed.Surfaces of Selective Value Revisited, Again
Imagine a protein undergoing evolution driven by positive selection through the sequential replacement of individual amino acids, in a haploid population sufficiently large that random genetic drift can be ignored. After a single amino acid substitution at each of n distinct sites, the protein will have passed through an evolutionary pathway consisting of n successive alleles, each of which increased the fitness of the organism by some increment.Multipotency and Cell fate decision on the Epigenetic landscape: From Metaphor to Molecules and Mathematical Model
Each cell type in the metazoan body is characterized by its specific genome-wide profile of gene expression - which now can be approximated by measuring the transcriptome. How do cells, given the unfathomable combinatorial possibilities in the expression of ten thousands of genes, establish the correct transcriptome that specifies their fate with such stunning reliability and generates the discrete cell types? Obviously, the answer lies in the gene regulatory network that orchestrates the improbable yet inevitable establishment of each of the cell type-specific gene expression configurations. Conrad Waddington, not knowing about gene networks, proposed in the 1940s a now famous metaphor that captures the dynamics of cell differentiation into discrete and robust entities: The 'Epigenetic Landscape'. But what is the molecular and mathematical basis of this metaphor? In this talk, I will explain some intuitive theoretical concepts that may offer a formal basis for Waddington's landscape, and then, present experimental evidence in support of the ideas (i) that cell fates are "attractors" of the network (= Waddington's "valleys") and (ii) that cell fate decisions correspond to bifurcation-driven destabilization of the stem cell attractors (Waddington's "watersheds"). Finally I will show how transcriptome-wide gene expression fluctuations are at the core of multi-potency, allowing cells to vacillate in the state of indeterminacy before making fate decisions in response to external signals.Constructing self-replication system in liposome
To date, most knowledge regarding biological systems has been obtained by analyzing existing biological systems, such as cells or viruses. However, "analysis" is only a partial approach to understanding a system, which can be compensated by the constructive approach; the constructive approach is a method to construct a biological system from defined components. During the construction process, we can improve our understanding of the biological systems. We are focusing on a self-replication system of genetic information, which is a fundamental characteristic of biological systems.Shuffled gene networks and liar paradoxes
Synthetic biology aims to understand biological systems through reconstruction. We tested the limits of reconstruction by "rewiring" a large-scale network - the transcription system in Escherichia coli. By constructing 598 recombinations of promoters (including regulatory regions) with different transcription or sigma-factor genes, new network connections were added over the wild-type genetic background. This revealed both the robustness of the underlying network and its propensity to evolve new properties. In the process of shuffling bacterial transcription networks, many new network motifs were created, involving both positive and negative feedback loops. Commonly, such networks are illustrated with topology diagrams - where pointed or blunt arrows indicate network interactions. Because the same topology can have drastically different behaviours in different dimensions, it is important to emphasise that arrow diagrams are not sufficient to understand biological systems. Unless the dimensions of time and space are considered explicitly, these 'snapshot' diagrams contain inherent contradictions and resemble the classical Greek philosophy problems known as liar paradoxes.Macroscopic relationship among plasticity, robustness, evolvablity, and phenotypic fluctuations: Waddinton's legacy revisited under the spirit of Einstein
Characterization of biological plasticity, robustness, and evolvability in terms of dynamical systems and statistical analysis is an important issue in Complex Systems Biology.Adaptive Response of a Gene Network to Environmental Changes by Fitness-induced Attractor Selection
Cells switch between various stable genetic programs (attractors) to accommodate environmental conditions. Signal transduction machineries efficiently convey environmental changes to the gene regulation apparatus in order to express the appropriate genetic program. However, since the number of environmental conditions is much larger than that of available genetic programs so that the cell may utilize the same genetic program for a large set of conditions, it may not have evolved a signaling pathway for every environmental condition, notably, those that are rarely encountered. Here we show that in the absence of signal transduction, switching to the appropriate attractor state expressing the genes that afford adaptation to the external condition can occur. In a synthetic bistable gene switch in Escherichia coli in which mutually inhibitory operons govern the expression of two genes required in two alternative nutritional environments, cells reliably selected the "adaptive attractor" driven by gene expression noise. A mathematical model suggests that the "non-adaptive attractor" is avoided because in unfavorable conditions, cellular activity is lower which suppresses mRNA metabolism, leading to larger fluctuations in gene expression. This in turn renders the non-adaptive state less stable. Although attractor selection is not as efficient as signal transduction via a dedicated cascade, it is simple and robust and may represent a primordial mechanism for adaptive responses that preceded the evolution of signaling cascades for the frequently encountered environmental changes.Resistance in multi-drug environments
The application of antibiotics is hindered by a well-known "catch-22": The use of a drug promotes the emergence and spread of drug-resistant mutants that ultimately render it ineffective. While certain combination therapies are known to be more effective than single drugs, the impact of such treatments on the evolution of drug resistance is unclear. In particular, very little is known about how the evolution of resistance is affected by the nature of the interactions - synergy or antagonism - between the drugs. I will describe a combined theoretical-experimental approach to understand drug interactions and their effect on the evolution of resistance in bacteria. Our results indicate that synergistic drug pairs, typically preferred in clinical settings, can accelerate bacterial adaptation. Drug antagonism, on the other hand, can generate selection against resistant bacteria, thereby slowing down the evolution of resistance. These results suggest a tradeoff in drug-combination therapy between immediate inhibition of growth and long-term inhibition of resistance.Exploring the fitness landscape of the protein translation system
The protein translation reaction, one of the most important regulators of cell behavior, involves the interactions of a large number of components, and has been studied extensively because of its importance in the cell. By reconstruction of an Escherichia coli-based in vitro translation system with protein components highly purified on an individual basis, it has been demonstrated that 36 enzymes and ribosomes are sufficient to carry out protein translation [1]. Thus, the minimum number of components required for the reaction has been identified. Systems in which the components interact with each other are nonlinear, which makes it difficult to link the component concentrations and the synthesis activity, and therefore the relationship between the two remains to be elucidated. In a system composed of n components, the activity of the system (activity) is written as a function of the component concentrations (c1, c2, ..., cn):A "pattern language" for evolution and development of animal form
Ancient animals arose more than half a billion years ago from unicellular organisms with sophisticated genetic repertoires. Unexpectedly from the neo-Darwinian perspective, most of the genes of the metazoan "developmental-genetic toolkit," the products of which coordinate multicellular development, were present in these single-celled forms. The rapidity with which animal developmental mechanisms emerged during the Precambrian-Cambrian transition suggests that more than just incremental adaptive evolution was involved. Here I will consider the role played by a core set of "dynamical patterning modules" (DPMs) in the origination, development and evolution of multicellular animals. DPMs are functional modules that consist of the products of certain toolkit genes in association with physical processes and effects characteristic of chemically and mechanically excitable systems of the "mesoscale" (i.e., linear dimension ~0.1-10 mm). Once cellular life achieved this spatial scale by the most basic DPM, cell-cell adhesion, a variety of physical forces and effects came into play, including cohesion, viscoelasticity, diffusion, spatiotemporal heterogeneity based on lateral inhibition, and global synchronization of oscillatory dynamics. I will show how toolkit gene products and pathways that pre-existed the metazoans acquired novel morphogenetic functions simply by virtue of the change in scale and context inherent to multicellularity. DPMs, acting singly and in combination with each other, will be seen to constitute a "pattern language" capable of generating all metazoan body plans and organ forms. The stable developmental trajectories ("developmental programs") and morphological phenotypes of modern organisms, in this view, arose by stabilizing selection that routinized the generation of morphological motifs which were originally manifestations of the physical properties of multicellular aggregates. With respect to the origination of multicellular form, this perspective diverges from the tenets of the Modern Synthesis by relinquishing uniformitarianism, incrementalism, and a central role for adaptive directional selection.A cell-density dependent transition to collective behavior in Dictyostelium
Coherent dynamics in populations of microorganisms such as collective migration and synchronized oscillations often emerge as a result of communication among individual cells via cell-cell signaling. Recent studies in chemical oscillators suggest a possibility that a dynamical density-sensing maybe at work whereby the coupling between the excitable cells gives rise to synchronized oscillations whose frequency encodes cell-density. Although the scheme provides a robust means to encode cell-density information as a frequency in the periodic intracellular signaling, whether such mechanism is at work remain elusive. The chemoattractant signaling in social amoebae Dictyostelium discoideum is one of the best known examples of synchronized oscillations. Chemoattractant cAMP is synthesized and secreted periodically and serves as a cue that directs chemotaxis of individual cells. Our recent live-cell FRET imaging revealed the very onset of oscillations and sequential changes in the frequency at the single cell level resolution. Together with a mathematical modeling, we demonstrate that the cells prefer inputs of a certain resonant frequency and that this frequency is determined by cell-intrinsic feedback. In the talk, I will focus on how such intracellular dynamics gives rise to the transition between decoupled and collective states of the periodic cAMP signaling. We show from a series of perfusion experiments that the transition is abrupt and the frequency of the oscillations depends strictly on cell density. Similar quorum-sensing schemes are also thought to be at work in glycolytic oscillations in yeast and maybe prevalent in other organisms.A Road to Artificial Cell
INTRODUCTIONMathematical structure of the evolution of cell differentiation in clonal cell groups
Multicellular organisms are composed of a number of cell types which are produced through cell differentiation processes from clonal cell populations or single fertilized cells. Since multicellularity has evolved several times independently in different lineages, universality in these evolutions should be taken into consideration. In particular, here we investigate the evolutionary establishment of cell differentiation in clonal cell groups. How the evolutionary processes are characterized mathematically? Are they classified into several types based on the characterization? In this study, we attempt to classify the evolution of cell differentiation by bifurcation theory. Suppose a dynamical system which (potentially) controls the cell differentiation, the evolutionary establishment and/or change of differentiation process can be characterized by the bifurcation of the dynamical system. In order to test the idea, we performed evolutionary simulations of differentiation for cell groups with gene regulatory network models which describe cell state dynamics. Clonal cell groups that develop into same cell types were subjected to selection pressure to produce two cell types whose population ratio is controlled to become about 1:1. We then found two clearly distinguishable evolutionary pathways; i) a function to produce two cell types was achieved first, and then a function to control the population ratio was achieved, ii) both functions were achieved simultaneously. The results indicate that the evolutionary relation of the two cellular functions noted above depends on the situation. To verify that these evolutionary pathways are two distinct processes, we identified the bifurcations (mathematical description for qualitative changes in the behavior of a dynamical system) involved in these pathways. The analysis reveals that the above two pathways correspond to two different bifurcation types of cellular systems; i.e. these pathways are characterized by qualitatively different mathematical structures. Our study suggests the flexibility of evolution of differentiation, and proposes the concept of applying the bifurcation analysis for characterization of flexible evolutionary pathways.On the relationship between robustness and evolutionary innovation
Mutational robustness and evolvability, a system's ability to produce heritable phenotypic variation, harbour a paradoxical tension. On one hand, high robustness implies low production of heritable phenotypic variation. On the other hand, both experimental and computational analyses of neutral networks indicate that robustness enhances evolvability. I here resolve this tension using a molecular study system. To resolve the tension, one must distinguish between robustness of a genotype and a phenotype. I confirm that genotype (sequence) robustness and evolvability share an antagonistic relationship. In stark contrast, phenotype (structure) robustness promotes structure evolvability. A consequence is that finite populations of sequences with a robust phenotype can access large amounts of phenotypic variation while spreading through a neutral network. Population-level processes and phenotypes rather than individual sequences are key to understand the relationship between robustness and evolvability. My observations may apply to other genetic systems where many connected genotypes produce the same phenotypes.Bacterial persistence against antibiotics
Bacterial persistence is an epigenetic phenomenon in which a subset of cells in an isogenic population shows phenotypic tolerance to bactericidal drugs and is able to survive prolonged exposure to the drugs. We analyzed a persistence phenomenon in Mycobacterium smegmatis generated on exposure to the antibiotic, isoniazid (INH), by time-lapse single cell microscopy utilizing a microfluidic device. Our results demonstrate that persistence trait in M. smegmatis against INH is independent of the growth state of the bacteria prior to the drug exposure, which counteracts the widely accepted hypothesis for the INH persistence of mycobacteria that persisters are pre-existing non-growing cells selected on exposure to the drug. Furthermore, on INH exposure they continue to grow and divide at suppressed rates while the killing by the drug is ongoing; persistence phase, where a number of surviving cells is nearly unchanged at the population-level, is in fact quite dynamic at the single-cell level and emerges by balancing the division and killing rates in a population.Phenotypic fluctuation in experimental evolution
Recent studies on stochastic gene expression of individual cells have uncovered the large degree of phenotypic fluctuation, thus leading substantial phenotypic diversity even under identical environments. What does the phenotypic fluctuation affect survival in Darwinian evolution? We have conducted evolution experiments of random mutation on GFP gene in E.coli and selection based on expressed green fluorescence as phenotype in two different methods. In each cycle of the first method, the mutant gene with the highest "average" green fluorescence was selected in excluding the effect of cell-to-cell variation. As expected, the fitness (the average green fluorescence) of the mutant selected genotype increased while the cell-to-cell variation of each of the selected genotypes became smaller. But, interestingly, the linear correlation between the evolution rate and the cell-to-cell variation in the phenotype was found (Sato K, et al. (2003)). In contrast to the above selection based on the "average" phenotype, in the second selection method where 0.2% of individual cells were selected for higher green fluorescence, some of the selected genotypes showed a low average green fluorescence but larger cell-to-cell variations in their isogenic populations. Because the cells of each of the selected mutants, even with their low averages, occasionally could show higher green fluorescence by phenotypic fluctuation, some of the cells survived under the strong selection of the top 0.2% in the phenotype distribution. These mutant genotypes of larger fluctuation appeared independently on the different branches of the phylogenic tree. In addition to the average phenotype change by genetic mutation, the observed increase in phenotypic fluctuation acts as an evolutionary strategy to produce an extreme phenotype under severe selective environments (Ito Y, et al. (2009)).