Various animals, including humans, have been suggested to perform Bayesian inferences to handle noisy, time-varying external information. In performing Bayesian inference by the brain, the prior distribution must be acquired and represented by sampling noisy external inputs. However, the mechanism by which neural activities represent such distributions has not yet been elucidated. Our findings reveal that networks with modular structures, composed of fast and slow modules, are adept at representing this prior distribution, enabling more accurate Bayesian inferences. Specifically, the modular network that consists of a main module connected with input and output layers and a sub-module with slower neural activity connected only with the main module outperformed networks with uniform time scales. Prior information was represented specifically by the slow sub-module, which could integrate observed signals over an appropriate period and represent input means and variances. Accordingly, the neural network could effectively predict the time-varying inputs. Furthermore, by training the time scales of neurons starting from networks with uniform time scales and without modular structure, the above slow-fast modular network structure and the division of roles in which prior knowledge is selectively represented in the slow sub-modules spontaneously emerged. These results explain how the prior distribution for Bayesian inference is represented in the brain, provide insight into the relevance of modular structure with time scale hierarchy to information processing, and elucidate the significance of brain areas with slower time scales.

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Cell polarity regulates the orientation of the cytoskeleton members that directs intracellular transport for cargo-like organelles, using chemical gradients sustained by ATP or GTP hydrolysis. However, how cargo transports are directly mediated by chemical gradients remains unknown. We previously proposed a physical mechanism that enables directed movement of cargos, referred to as chemophoresis. According to the mechanism, a cargo with reaction sites is subjected to a chemophoresis force in the direction of the increased concentration. Based on this, we introduce an extended model, the chemophoresis engine, as a general mechanism of cargo motion, which transforms chemical free energy into directed motion through the catalytic ATP hydrolysis. We applied the engine to plasmid motion in a ParABS system to demonstrate the self-organization system for directed plasmid movement and pattern dynamics of ParA-ATP concentration, thereby explaining plasmid equi-positioning and pole-to-pole oscillation observed in bacterial cells and in vitro experiments. We mathematically show the existence and stability of the plasmid-surfing pattern, which allows the cargo-directed motion through the symmetry-breaking transition of the ParA-ATP spatiotemporal pattern. We also quantitatively demonstrate that the chemophoresis engine can work even under in vivo conditions. Finally, we discuss the chemophoresis engine as one of the general mechanisms of hydrolysis-driven intracellular transport.

Robustness and plasticity are essential features that allow biological systems to cope with complex and variable environments. In a constant environment, robustness, i.e., insensitivity of phenotypes, is expected to increase, whereas plasticity, i.e., the changeability of phenotypes, tends to diminish. Under a variable environment, existence of plasticity will be relevant. The robustness and plasticity, on the other hand, are related to phenotypic variances. As phenotypic variances decrease with the increase in robustness to perturbations, they are expected to decrease through the evolution. However, in nature, phenotypic fluctuation is preserved to a certain degree. One possible cause for this is environmental variation, where one of the most important "environmental" factors will be inter-species interactions. As a first step toward investigating phenotypic fluctuation in response to an inter-species interaction, we present the study of a simple two-species system that comprises hosts and parasites. Hosts are expected to evolve to achieve a phenotype that optimizes fitness. Then, the robustness of the corresponding phenotype will be increased by reducing phenotypic fluctuations. Conversely, plasticity tends to evolve to avoid certain phenotypes that are attacked by parasites. By using a dynamic model of gene expression for the host, we investigate the evolution of the genotype-phenotype map and of phenotypic variances. If the host-parasite interaction is weak, the fittest phenotype of the host evolves to reduce phenotypic variances. In contrast, if there exists a sufficient degree of interaction, the phenotypic variances of hosts increase to escape parasite attacks. For the latter case, we found two strategies: if the noise in the stochastic gene expression is below a certain threshold, the phenotypic variance increases via genetic diversification, whereas above this threshold, it is increased mediated by noise-induced phenotypic fluctuation. We examine how the increase in the phenotypic variances caused by parasite interactions influences the growth rate of a single host, and observed a trade-off between the two. Our results help elucidate the roles played by noise and genetic mutations in the evolution of phenotypic fluctuation and robustness in response to host-parasite interactions.

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.

Proteins in cellular environments are highly susceptible. Local perturbations to any residue can be sensed by other spatially distal residues in the protein molecule, showing long-range correlations in the native dynamics of proteins. The long-range correlations of proteins contribute to many biological processes such as allostery, catalysis, and transportation. Revealing the structural origin of such long-range correlations is of great significance in understanding the design principle of biologically functional proteins. In this work, based on a large set of globular proteins determined by X-ray crystallography, by conducting normal mode analysis with the elastic network models, we demonstrate that such long-range correlations are encoded in the native topology of the proteins. To understand how native topology defines the structure and the dynamics of the proteins, we conduct scaling analysis on the size dependence of the slowest vibration mode, average path length, and modularity. Our results quantitatively describe how native proteins balance between order and disorder, showing both dense packing and fractal topology. It is suggested that the balance between stability and flexibility acts as an evolutionary constraint for proteins at different sizes. Overall, our result not only gives a new perspective bridging the protein structure and its dynamics but also reveals a universal principle in the evolution of proteins at all different sizes.

The emergence of replicases that can replicate themselves is a central issue in the origin of life. Recent experiments suggest that such replicases can be realized if an RNA polymerase ribozyme is divided into fragments short enough to be replicable by the ribozyme and if these fragments self-assemble into a functional ribozyme. However, the continued self-replication of such replicases requires that the production of every essential fragment be balanced and sustained. Here, we use mathematical modeling to investigate whether and under what conditions fragmented replicases achieve continued self-replication. We first show that under a simple batch condition, the replicases fail to display continued self-replication owing to positive feedback inherent in these replicases. This positive feedback inevitably biases replication toward a subset of fragments, so that the replicases eventually fail to sustain the production of all essential fragments. We then show that this inherent instability can be resolved by small rates of random content exchange between loose compartments (i.e., horizontal transfer). In this case, the balanced production of all fragments is achieved through negative frequency-dependent selection operating in the population dynamics of compartments. The horizontal transfer also ensures the presence of all essential fragments in each compartment, sustaining self-replication. Taken together, our results underline compartmentalization and horizontal transfer in the origin of the first self-replicating replicases.

A most important property of biochemical systems is robustness. Static robustness, e.g., homeostasis, is the insensitivity of a state against perturbations, whereas dynamics robustness, e.g., homeorhesis, is the insensitivity of a dynamic process. In contrast to the extensively studied static robustness, dynamics robustness, i.e., how a system creates an invariant temporal profile against perturbations, is little explored despite transient dynamics being crucial for cellular fates and are reported to be robust experimentally. For example, the duration of a stimulus elicits different phenotypic responses, and signaling networks process and encode temporal information. Hence, robustness in time courses will be necessary for functional biochemical networks. Based on dynamical systems theory, we uncovered a general mechanism to achieve dynamics robustness. Using a three-stage linear signaling cascade as an example, we found that the temporal profiles and response duration post-stimulus is robust to perturbations against certain parameters. Then analyzing the linearized model, we elucidated the criteria of when signaling cascades will display dynamics robustness. We found that changes in the upstream modules are masked in the cascade, and that the response duration is mainly controlled by the rate-limiting module and organization of the cascade's kinetics. Specifically, we found two necessary conditions for dynamics robustness in signaling cascades: 1) Constraint on the rate-limiting process: The phosphatase activity in the perturbed module is not the slowest. 2) Constraints on the initial conditions: The kinase activity needs to be fast enough such that each module is saturated even with fast phosphatase activity and upstream changes are attenuated. We discussed the relevance of such robustness to several biological examples and the validity of the above conditions therein. Given the applicability of dynamics robustness to a variety of systems, it will provide a general basis for how biological systems function dynamically.

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.

Embryonic stem cells exhibit pluripotency: they can differentiate into all types of somatic cells. Pluripotent genes such as Oct4 and Nanog are activated in the pluripotent state, and their expression decreases during cell differentiation. Inversely, expression of differentiation genes such as Gata6 and Gata4 is promoted during differentiation. The gene regulatory network controlling the expression of these genes has been described, and slower-scale epigenetic modifications have been uncovered. Although the differentiation of pluripotent stem cells is normally irreversible, reprogramming of cells can be experimentally manipulated to regain pluripotency via overexpression of certain genes. Despite these experimental advances, the dynamics and mechanisms of differentiation and reprogramming are not yet fully understood. Based on recent experimental findings, we constructed a simple gene regulatory network including pluripotent and differentiation genes, and we demonstrated the existence of pluripotent and differentiated states from the resultant dynamical-systems model. Two differentiation mechanisms, interaction-induced switching from an expression oscillatory state and noise-assisted transition between bistable stationary states, were tested in the model. The former was found to be relevant to the differentiation process. We also introduced variables representing epigenetic modifications, which controlled the threshold for gene expression. By assuming positive feedback between expression levels and the epigenetic variables, we observed differentiation in expression dynamics. Additionally, with numerical reprogramming experiments for differentiated cells, we showed that pluripotency was recovered in cells by imposing overexpression of two pluripotent genes and external factors to control expression of differentiation genes. Interestingly, these factors were consistent with the four Yamanaka factors, Oct4, Sox2, Klf4, and Myc, which were necessary for the establishment of induced pluripotent stem cells. These results, based on a gene regulatory network and expression dynamics, contribute to our wider understanding of pluripotency, differentiation, and reprogramming of cells, and they provide a fresh viewpoint on robustness and control during development.

Cellular memory, which allows cells to retain information from their environment, is important for a variety of cellular functions, such as adaptation to external stimuli, cell differentiation, and synaptic plasticity. Although posttranslational modifications have received much attention as a source of cellular memory, the mechanisms directing such alterations have not been fully uncovered. It may be possible to embed memory in multiple stable states in dynamical systems governing modifications. However, several experiments on modifications of proteins suggest long-term relaxation depending on experienced external conditions, without explicit switches over multi-stable states. As an alternative to a multistability memory scheme, we propose "kinetic memory" for epigenetic cellular memory, in which memory is stored as a slow-relaxation process far from a stable fixed state. Information from previous environmental exposure is retained as the long-term maintenance of a cellular state, rather than switches over fixed states. To demonstrate this kinetic memory, we study several models in which multimeric proteins undergo catalytic modifications (e.g., phosphorylation and methylation), and find that a slow relaxation process of the modification state, logarithmic in time, appears when the concentration of a catalyst (enzyme) involved in the modification reactions is lower than that of the substrates. Sharp transitions from a normal fast-relaxation phase into this slow-relaxation phase are revealed, and explained by enzyme-limited competition among modification reactions. The slow-relaxation process is confirmed by simulations of several models of catalytic reactions of protein modifications, and it enables the memorization of external stimuli, as its time course depends crucially on the history of the stimuli. This kinetic memory provides novel insight into a broad class of cellular memory and functions. In particular, applications for long-term potentiation are discussed, including dynamic modifications of calcium-calmodulin kinase II and cAMP-response element-binding protein essential for synaptic plasticity.

Cells generally adapt to environmental changes by first exhibiting an immediate response and then gradually returning to their original state to achieve homeostasis. Although simple network motifs consisting of a few genes have been shown to exhibit such adaptive dynamics, they do not reflect the complexity of real cells, where the expression of a large number of genes activates or represses other genes, permitting adaptive behaviors. Here, we investigated the responses of gene regulatory networks containing many genes that have undergone numerical evolution to achieve high fitness due to the adaptive response of only a single target gene; this single target gene responds to changes in external inputs and later returns to basal levels. Despite setting a single target, most genes showed adaptive responses after evolution. Such adaptive dynamics were not due to common motifs within a few genes; even without such motifs, almost all genes showed adaptation, albeit sometimes partial adaptation, in the sense that expression levels did not always return to original levels. The genes split into two groups: genes in the first group exhibited an initial increase in expression and then returned to basal levels, while genes in the second group exhibited the opposite changes in expression. From this model, genes in the first group received positive input from other genes within the first group, but negative input from genes in the second group, and vice versa. Thus, the adaptation dynamics of genes from both groups were consolidated. This cooperative adaptive behavior was commonly observed if the number of genes involved was larger than the order of ten. These results have implications in the collective responses of gene expression networks in microarray measurements of yeast Saccharomyces cerevisiae and the significance to the biological homeostasis of systems with many components.

Recent experimental measurements have demonstrated that spontaneous neural activity in the absence of explicit external stimuli has remarkable spatiotemporal structure. This spontaneous activity has also been shown to play a key role in the response to external stimuli. To better understand this role, we proposed a viewpoint, "memories-as-bifurcations," that differs from the traditional "memories-as-attractors" viewpoint. Memory recall from the memories-as-bifurcations viewpoint occurs when the spontaneous neural activity is changed to an appropriate output activity upon application of an input, known as a bifurcation in dynamical systems theory, wherein the input modifies the flow structure of the neural dynamics. Learning, then, is a process that helps create neural dynamical systems such that a target output pattern is generated as an attractor upon a given input. Based on this novel viewpoint, we introduce in this paper an associative memory model with a sequential learning process. Using a simple hebbian-type learning, the model is able to memorize a large number of input/output mappings. The neural dynamics shaped through the learning exhibit different bifurcations to make the requested targets stable upon an increase in the input, and the neural activity in the absence of input shows chaotic dynamics with occasional approaches to the memorized target patterns. These results suggest that these dynamics facilitate the bifurcations to each target attractor upon application of the corresponding input, which thus increases the capacity for learning. This theoretical finding about the behavior of the spontaneous neural activity is consistent with recent experimental observations in which the neural activity without stimuli wanders among patterns evoked by previously applied signals. In addition, the neural networks shaped by learning properly reflect the correlations of input and target-output patterns in a similar manner to those designed in our previous study.

Precise temporal coordination of gene expression is crucial for many developmental processes. One central question in developmental biology is how such coordinated expression patterns are robustly controlled. During embryonic development of the Drosophila central nervous system, neural stem cells called neuroblasts express a group of genes in a definite order, which leads to the diversity of cell types. We produced all possible regulatory networks of these genes and examined their expression dynamics numerically. From the analysis, we identified requisite regulations and predicted an unknown factor to reproduce known expression profiles caused by loss-of-function or overexpression of the genes in vivo, as well as in the wild type. Following this, we evaluated the stability of the actual Drosophila network for sequential expression. This network shows the highest robustness against parameter variations and gene expression fluctuations among the possible networks that reproduce the expression profiles. We propose a regulatory module composed of three types of regulations that is responsible for precise sequential expression. This study suggests that the Drosophila network for sequential expression has evolved to generate the robust temporal expression for neuronal specification.

How can a microorganism adapt to a variety of environmental conditions despite the existence of a limited number of signal transduction mechanisms? We show that for any growing cells whose gene expression fluctuate stochastically, the adaptive cellular state is inevitably selected by noise, even without a specific signal transduction network for it. In general, changes in protein concentration in a cell are given by its synthesis minus dilution and degradation, both of which are proportional to the rate of cell growth. In an adaptive state with a higher growth speed, both terms are large and balanced. Under the presence of noise in gene expression, the adaptive state is less affected by stochasticity since both the synthesis and dilution terms are large, while for a nonadaptive state both the terms are smaller so that cells are easily kicked out of the original state by noise. Hence, escape time from a cellular state and the cellular growth rate are negatively correlated. This leads to a selection of adaptive states with higher growth rates, and model simulations confirm this selection to take place in general. The results suggest a general form of adaptation that has never been brought to light--a process that requires no specific mechanisms for sensory adaptation. The present scheme may help explain a wide range of cellular adaptive responses including the metabolic flux optimization for maximal cell growth.

In the past few years, modeling and querying video databases have been a subject of extensive research to develop tools for effective search of videos. In this paper, we present a hierarchal approach to model videos at three levels, object level (OL), frame level (FL), and shot level (SL). The model captures the visual features of individual objects at OL, visual-spatio-temporal (VST) relationships between objects at FL, and time-varying visual features and time-varying VST relationships at SL. We call the combination of the time-varying visual features and the time-varying VST relationships a Content trajectory which is used to represent and index a shot. A novel query interface that allows users to describe the time-varying contents of complex video shots such as those of skiers, soccer players, etc., by sketch and feature specification is presented. Our experimental results prove the effectiveness of modeling and querying shots using the content trajectory approach.

Feature-based similarity retrieval has become an important research issue in multimedia database systems. The features of multimedia data are useful for discriminating between multimedia objects (e.g., documents, images, video, music score, etc.). For example, images are represented by their color histograms, texture vectors, and shape descriptors. A feature vector is a vector that represents a set of features, and are usually high-dimensional data. The performance of conventional multidimensional data structures (e.g., R-tree family K-D-B tree, grid file, TV-tree) tends to deteriorate as the number of dimensions of feature vectors increases. The R*-tree is the most successful variant of the R-tree. We propose a SOM-based R*-tree as a new indexing method for high-dimensional feature vectors. The SOM-based R*-tree combines SOM and R*-tree to achieve search performance more scalable to high dimensionalities. Self-organizing maps (SOMs) provide mapping from high-dimensional feature vectors onto a two-dimensional space. The mapping preserves the topology of the feature vectors. The map is called a topological feature map, and preserves the mutual relationships (similarity) in the feature spaces of input data, clustering mutually similar feature vectors in neighboring nodes. We experimentally compare the retrieval time cost of a SOM-based R*-tree with that of an SOM and an R*-tree using color feature vectors extracted from 40,000 images.

Proposes a new parallel algorithm for computing path expressions, named the "parallel cascade semi-join" (PCSJ) algorithm. Moreover, a new scheduling strategy called the "right-deep zigzag tree" is designed to further improve the performance of the PCSJ algorithm. The experiments have been implemented in a distributed and parallel NOW (network of workstations) environment. The results show that the PCSJ algorithm outperforms two other parallel algorithms [the parallel forward pointer chasing (PFPC) algorithm and the index-splitting parallel algorithm (IndexSplit)] when computing path expressions with restrictive predicates, and that the right-deep zigzag tree scheduling strategy has a better performance than the right-deep tree scheduling strategy.

In order to support processing of highly concurrent transactions at low price on distributed database server systems, a network of workstations (NOW) is used as the hardware environment. In order to share resources that are distributed among different sites of NOW and support the database functionalities, persistent distributed shared virtual memory (PDSVM) is implemented on ShusseUo, an object database system developed by Kyushu University, Japan. In ShusseUo, all workstations cooperate to perform jobs submitted by database applications. Each job consists of several transactions. These transactions are executed on PDSVM and the cost of each transaction varies according to the workstation on which the transaction runs. We present the cost-based transaction coordinator (CTC) algorithm. In CTC, the load information of a transaction is collected automatically while the transaction is running, and it is fed back when the transaction is committed. In CTC, each transaction is coordinated to a certain workstation based on its cost as calculated using the fed back information and the distribution information of the database. The algorithm is evaluated in terms of the TPC-C benchmark. The benchmark result is presented and analyzed.

Earth science (ES) applications handle very large geo-spatial data sets and interactive response time is required by its query processing. Spatial selection is one of the very important basic operations for geo-spatial databases. It retrieves all the objects that intersect with a given point or rectangle. We present a novel approach for the parallel processing of spatial selection of very large geo-spatial databases using partitioned parallelism. To evaluate this approach, we use the Extended Sequoia 2000 benchmark, which has real world data and real queries. In addition, we use an actual object database management system, ShusseUo, which we developed previously. The experimental results of parallel processing of spatial selection show good speed-up.

This article lists fourteen open problems in artificial life, each of which is a grand challenge requiring a major advance on a fundamental issue for its solution. Each problem is briefly explained, and, where deemed helpful, some promising paths to its solution are indicated.

By introducing a dynamical system model of a multicellular system, it is shown that an organism with a variety of differentiated cell types and a complex pattern emerges through cell-cell interactions even without postulating any elaborate control mechanism. Such an organism is found to maintain a larger growth speed as an ensemble, by achieving a cooperative use of resources, than do simple homogeneous cells, which behave "selfishly." This suggests that the emergence of multicellular organisms with complex organization is a necessity in evolution. According to our theoretical model, there initially appear multipotent stem cells, which undergo stochastic differentiation to other cell types. With development and differentiation, both the chemical diversity and the complexity of intra-cellular dynamics are decreased, as a general consequence of our system. Robustness of the developmental process is also confirmed.

Evolution of genetic codes is studied as change in the choice of enzymes that are used to synthesize amino acids from the genetic information of nucleic acids. We propose the following theory: the differentiation of physiological states of a cell allows for a choice of enzymes, and this choice is later fixed genetically through evolution. To demonstrate this theory, a dynamical systems model consisting of the concentrations of metabolites, enzymes, amino acyl tRNA synthetase, and tRNA - amino acid complexes in a cell is introduced and studied numerically. It is shown that the biochemical states of cells are differentiated by cell-cell interactions, and each differentiated type starts to use a different synthetase. Through the mutation of genes, this difference in the genetic code is amplified and stabilized. The relevance of this theory to the evolution of non-universal genetic code in mitochondria is suggested. The present theory is based on our recent theory of isologous symbiotic speciation, which is briefly reviewed. According to the theory, phenotypes of organisms are first differentiated into distinct types through the interaction and developmental dynamics, even though they have identical genotypes; later, with mutation in the genotype, the genotype also differentiates into discrete types, while maintaining the "symbiotic" relationship between the types. Relevance of the theory to natural as well as artificial evolution is discussed.

Though parallel database systems have been extensively studied, as far as we know, the parallel algorithms of R-tree proposed so far are limited to one workstation with multiprocessors or multi disks, where a parallel sorting algorithm or concurrent I/O is used to improve the performance. For the searching of R-trees, multiple search paths from the root to leaves are traversed sequentially. This sequential traverse can be transformed into multiple parallel traverses based on multiple search paths, where the query is divided into subqueries which can be executed concurrently. Aiming at parallel I/O and CPU operations, we introduce a parallel R-tree search algorithm running on distributed shared virtual memory (DSVM), especially on Shusseuo which is an ODBMS providing global persistent object management on persistent DSVM. The related problems are discussed and the evaluations are made based on Shusseuo. Experimental results show that optimal performance can be reached in dealing with large volumes of data.

To understand an architecture of a living system, "Tile Automaton" is introduced as an abstract model of chemical reaction of molecules scattered over a space. The model consists of tiles of various shapes that stand for molecules. The chemical reaction, induced by the collisions of tiles, is represented by the change of the tile shapes. The rules for reaction are deterministic, and the evolution of the system strongly depends on mutual spatial relationship among tiles. The evolution often leads to self-organization of a "factory," a set of tiles that produces tiles continuously and keeps its structure. Several interesting phenomena, such as a deformation or a division of a factory, are also observed. It is proposed that the formation of the factory is due to the interference between different aspects of tiles - the shape and the motion. The concept of "entanglement" is introduced as a mechanism of living systems.

The origin of multicellular organisms and the mechanism of development in cell societies are studied by choosing a model with intracellular biochemical dynamics allowing for oscillations, cell-cell interaction through diffusive chemicals on a two-dimensional grid, and state-dependent cell adhesion. Cells differentiate due to a dynamical instability, as described by our "isologous diversification" theory. A fixed spatial pattern of differentiated cells emerges, where spatial information is sustained by cell-cell interactions. This pattern is robust against perturbations. With an adequate cell adhesion force, active cells are release that form the seed of a new generation of multicellular organisms, accompanied by death of the original multicellular unit as a halting state. It is shown that the emergence of multicellular organisms with differentiation, regulation, and life cycle is not an accidental event, but a natural consequence in a system of replicating cells with growth.

This paper presents the transaction management in a high performance distributed object storage system WAKASHI. Unlike other systems that use centralized client/server architecture and other conventional buffer management for distributed persistent object management, WAKASHI is based on symmetric peer-peer architecture and employs memory-mapping and distributed shared virtual memory techniques. Several novel techniques of transaction management for WAKASHI are developed. First, a multi-threaded transaction manager offers "multi-threaded connection" so that data control and transaction operations can be performed in parallel manner. Secondly, a concurrency control mechanism supports transparent page-level locks to reduce the complexity of user programs and locking overhead. Thirdly, a "compact commit" method is proposed to minimize the communication cost by reducing the amount of data and the number of connections. Fourthly, a redo-only recovery method is implemented by "shadowed cache" method to minimize the logging cost, and to allow fast recovery and system restart. Moreover, the system offers "hierarchical" control to support nested transactions. A performance evaluation by the OO7 benchmark is presented.

A nonzero-sum three-person coalition game is presented to study the evolution of complexity and diversity in cooperation, where the population dynamics of players with strategies is given according to their scores in the iterated game and mutations. Two types of differentiation emerge initially: a biased one to classes and a temporal one to change their roles for coalition. Rules to change the hands are self-organized in a society through evolution. The coevolution of diversity and complexity of strategies and interactions (or communications) are found at later stages of the simulation. Relevance of our results to the biological society is briefly discussed.

The relevance of chaos to evolution is discussed in the context of the origin and maintenance of diversity and complexity. Evolution to the edge of chaos is demonstrated in an imitation game. As an origin of diversity, dynamic clustering of identical chaotic elements, globally coupled each to the other, is briefly reviewed. The clustering is extended to nonlinear dynamics on hypercubic lattices, which enables us to construct a self-organizing genetic algorithm. A mechanism of maintenance of diversity, “homeochaos,” is given in an ecological system with interaction among many species. Homeochaos provides a dynamic stability sustained by high-dimensional weak chaos. A novel mechanism of cell differentiation is presented, based on dynamic clustering. Here, a new concept—“open chaos”—is proposed for the instability in a dynamical system with growing degrees of freedom. It is suggested that studies based on interacting chaotic elements can replace both top-down and bottom-up approaches.